 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 4897 | 回復(fù): 32 | |||||
llh2010至尊木蟲 (著名寫手)
|
[求助]
ATAT中g(shù)ensqs工具如何改變?nèi)〈鷿舛?
|
||||
| 請問各位高手們，使用ATAT軟件包的gensqs工具產(chǎn)生sqs結(jié)構(gòu)時，如何改變?nèi)〈臐舛�？例如對于一個bcc結(jié)構(gòu)，體心上的一種原子不變，而想讓立方體八個頂點上的另一種原子有一個或幾個被其它原子取代，請問應(yīng)該如何寫lat.in和conc.in文件。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 學(xué)習(xí) | VASP | 我的分享 |
» 猜你喜歡
286求調(diào)劑
已經(jīng)有9人回復(fù)
-
一志愿南昌大學(xué)324求調(diào)劑
已經(jīng)有6人回復(fù)
-
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可
已經(jīng)有10人回復(fù)
-
291求調(diào)劑
已經(jīng)有5人回復(fù)
-
308求調(diào)劑
已經(jīng)有7人回復(fù)
-
085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富
已經(jīng)有4人回復(fù)
-
壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！
已經(jīng)有3人回復(fù)
-
一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑
已經(jīng)有4人回復(fù)
-
083000學(xué)碩274求調(diào)劑
已經(jīng)有6人回復(fù)
-
304求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 金剛石工具的切削原理？
已經(jīng)有7人回復(fù)
- 用過冠醚的朋友們, 你們一般往有機溶劑中溶鹽都用多大濃度啊?
已經(jīng)有3人回復(fù)
- 表面活性劑F127和P123在水中，隨濃度變化，膠束形態(tài)改變的具體數(shù)據(jù)。
已經(jīng)有14人回復(fù)
- 血紅細胞的分離與固定所用的PBS濃度以及固定液
已經(jīng)有5人回復(fù)
- 旋光儀計算氨基酸及對映體的濃度
已經(jīng)有3人回復(fù)
- 有誰知道二硫化硒SeS2的濃度積常數(shù)Ksp值
已經(jīng)有13人回復(fù)
- 【求助】樹脂分子在乙醇中隨著濃度的降低而變渾析出的反常現(xiàn)象解釋
已經(jīng)有3人回復(fù)
linggang87
銀蟲 (小有名氣)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 48.4
- 散金: 60
- 紅花: 1
- 帖子: 198
- 在線: 108.8小時
- 蟲號: 635564
- 注冊: 2008-10-24
- 專業(yè): 金屬材料的制備科學(xué)與跨學(xué)
3樓2012-06-21 23:55:49
linggang87
銀蟲 (小有名氣)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 48.4
- 散金: 60
- 紅花: 1
- 帖子: 198
- 在線: 108.8小時
- 蟲號: 635564
- 注冊: 2008-10-24
- 專業(yè): 金屬材料的制備科學(xué)與跨學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★
fzx2008: 金幣+3, 謝謝分享 2012-06-22 08:59:52
llh2010: 金幣+5, ★★★很有幫助 2012-06-22 09:41:26
fzx2008: 金幣+3, 謝謝分享 2012-06-22 08:59:52
llh2010: 金幣+5, ★★★很有幫助 2012-06-22 09:41:26
|
我看到過一個人寫的ATAT產(chǎn)生SQS的筆記，貼在這里，你可以看看： SQS “Tutorial” Amy Bengtson 1 June 2007 Here are my notes on setting up SQS (with many thanks from Yueh-Lin for teaching me). Setting up the most random SQS cells can be subtle, and my notes aren’t perfect, so please make sure you double check everything here. I make special note of the files you may want to double check with Dane. For more explanations, please see the atat manual (section 6.1.8): http://www.its.caltech.edu/~avdw/atat/ http://www.its.caltech.edu/~avdw/atat/manual.pdf I am using CoFe as an example since it is an SQS structure I have set up in the past. In order to get sqs structure for 50% Fe in BCC CoFe: 1. Create a lat.in file. This is the structure from which you want to get the disorded structure. lat.in is of this format: 2.810643 0.000000 0.000000 0.000000 2.810643 0.000000 0.000000 0.000000 2.810643 0.500000 0.500000 0.500000 0.500000 -0.500000 0.500000 0.500000 0.500000 -0.500000 0.000000 0.000000 0.000000 Co,Fe This lat.in is for primitive BCC with one atom. If your primitive cell has more than 1 atom, then continue to list all atomic positions after 0.000000 0.000000 0.000000 Co,Fe. 2. Make conc.in. This tells ATAT what concentration you want to create. conc.in is of the same format as lat.in, only with all atomic positions listed. 2.810643 0.000000 0.000000 0.000000 2.810643 0.000000 0.000000 0.000000 2.810643 -1.000000 0.000000 0.000000 0.000000 0.000000 1.000000 0.000000 1.000000 0.000000 -1.000000 1.000000 1.000000 Fe -0.500000 0.500000 0.500000 Co 3. Now create the clusters by using this command: corrdump -l lat.in -2=maxradius -clus The maxradius is the length of the longest pair desired. A reasonable starting guess is: -2=3 4. Find the correlation of the clusters you created in step 3. Here, maxradius should be the same number as in 3. corrdump -noe -2=maxradius -rnd -s=conc.in > tcorr.out 5. Generate the sqs structure (disorded structures) using this command: gensqs -n=Natoms > sqs.out Where Natoms = the size of the cell you want to create. sqs.out will contain all of the sqs structures, there may be more than 1! The structures will be in the formation of str.out (like conc.in – some conversions will be needed to make them into POSCARs – see step *). 6. In many cases you will find you have many structures in sqs.out. The problem is how to choose which structures to run. You can do two things: a. Redo everything above with a larger maxradius – this will give more correlations to match. Or you can add another radius to match: corrdump -l lat.in -2=maxradius -3=another_radius –clus corrdump -noe -2=maxradius -3=another_radius -rnd -s=conc.in > tcorr.out The “-2” matches pairs, the “-3” matches triplets. b. Rank the sqs structures that you have in sqs.out to find the “most random” structures. To do this: i. Find the correlation on the sqs structures in sqs.out: corrdump -noe -2=LargerRadius -s=sqs.out > tcorr_final.out You want LargerRadius to be larger than maxradius because you are trying find the correlations beyond the original correlation given in tcorr.out. Usually you want to rank based on the first 3 pairs. The columns in tcorr_final.out are: [point correlation] [1st pair] [2nd pair] [3rd pair] [4th pair] … If you don’t have at least 4 columns in your tcorr_final.out, then you need to increase LargerRadius until you have at least 4 columns. Each row corresponds to an sqs structure in sqs.out ii. Find the target correlation out to the LargerRadius corrdump -noe -2= LargerRadius -rnd -s=sqs.out > tcorr_finalRND.out iii. Open tcorr_final.out in excel and rank the structures based on most random. 7. Once you have found which structures from sqs.out you are going to use, split them up and make each structure into its own str.out file. 8. Create a vasp.wrap file that gives the VASP input information. Here is an example: [INCAR] SYSTEM = CoFe ENCUT=455 ISPIND = 2 #makes spin-polarized calc. possible ISPIN = 2 #does spin-polarized calc. MAGMOM= 5 5 5 5 5 5 5 5 ISTART = 0 INIWAV = 1 NSW = 191 IBRION = 2 ISIF = 3 ISMEAR = 1 SIGMA = 0.2 PREC = Accurate LWAVE = .FALSE. KPPRA = 17576 KSCHEME=Monkorst-Pack DOGGA SUBATOM = s/Co/Co/g SUBATOM = s/Fe/Fe_pv/g 9. In the directory with sqs.out, create the VASP input files using this command: runstruct_vasp -nr 10. Run as you normally do. If you make a mistake and need to start over, remove all of these files: rm clusters.out corrdump.log gensqs_0_1.stat sqscell.out sqs* sym.out tcorr.out Important!!! These are the files you should double check with Dane to make sure they were set up correctly. o lat.in o conc,in o tcorr.out o rms ranking (show the Excel sheet where you rank tcorr_final.out) o And confirm that maxradius is a reasonable value. |
4樓2012-06-22 00:17:13
tt-0-8
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 676.8
- 帖子: 193
- 在線: 127.2小時
- 蟲號: 1315323
- 注冊: 2011-06-04
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
2樓2012-05-11 16:01:15
llh2010
至尊木蟲 (著名寫手)
- 應(yīng)助: 34 (小學(xué)生)
- 金幣: 10007.8
- 散金: 14
- 紅花: 11
- 帖子: 1614
- 在線: 350.1小時
- 蟲號: 1471268
- 注冊: 2011-11-01
- 性別: GG
- 專業(yè): 原子和分子物理
5樓2012-06-22 09:40:27
llh2010
至尊木蟲 (著名寫手)
- 應(yīng)助: 34 (小學(xué)生)
- 金幣: 10007.8
- 散金: 14
- 紅花: 11
- 帖子: 1614
- 在線: 350.1小時
- 蟲號: 1471268
- 注冊: 2011-11-01
- 性別: GG
- 專業(yè): 原子和分子物理
★ ★
fzx2008: 金幣+2, 鼓勵交流 2012-06-22 10:30:20
fzx2008: 金幣+2, 鼓勵交流 2012-06-22 10:30:20
|
周期性結(jié)構(gòu)中八個頂角是等價的，在隨機排布中一個單胞中的八個頂角可能排不同的原子。隨機排布的晶體很多情況下只具有短程有序。 你說的方法我也看到其他人做過，這種方法計算出來的結(jié)果不是很準(zhǔn)，原因是這種方法建出來的超胞具有周期性，雖然在一個超胞中是可以隨機取代，但是這個超胞之外的其它超胞就都和它一樣了。而實際情況是所有的超胞都是隨機取代的，和所建的超胞可能相同，也可能不相同。為了克服這個問題，往往需要建很大的超胞，以便達到統(tǒng)計效果，但是問題又來了，機器可能算不動這么大的超胞。 |
6樓2012-06-22 10:01:05
linggang87
銀蟲 (小有名氣)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 48.4
- 散金: 60
- 紅花: 1
- 帖子: 198
- 在線: 108.8小時
- 蟲號: 635564
- 注冊: 2008-10-24
- 專業(yè): 金屬材料的制備科學(xué)與跨學(xué)
★ ★
franch: 金幣+2, 謝謝回帖交流。 2012-06-25 17:29:51
franch: 金幣+2, 謝謝回帖交流。 2012-06-25 17:29:51
|
看來樓主可能不是用超晶胞方法材料計算方向的。SQS是基于超晶胞方法，用相對較小的晶胞去模擬無序情況，這要求SQS超胞內(nèi)各個原子的周圍環(huán)境‘盡可能’和無序時相近，而且誤差可以衡量。 超晶胞方法進行晶體學(xué)計算中，周期性是必須的，而且如果你用的超胞能夠模擬無序情況了，這個超胞之外的和它等價的那些自然也可以模擬無序。 超晶胞方法模擬不同濃度時是有限制的，比如2%濃度，你只能建一個50個原子的胞，取代其中一個原子，進行模擬，你不能取代半個原子。對于這個限制VCA 和CPA模擬無序就有優(yōu)勢了，濃度可以隨意，不過VCA和CPA都不能進行原子位置的弛豫，有時候這是很大的劣勢。 |
» 本帖已獲得的紅花（最新10朵）
llh20107樓2012-06-22 19:58:45
llh2010
至尊木蟲 (著名寫手)
- 應(yīng)助: 34 (小學(xué)生)
- 金幣: 10007.8
- 散金: 14
- 紅花: 11
- 帖子: 1614
- 在線: 350.1小時
- 蟲號: 1471268
- 注冊: 2011-11-01
- 性別: GG
- 專業(yè): 原子和分子物理
8樓2012-06-23 14:26:45
souledge
專家顧問 (著名寫手)
-
專家經(jīng)驗: +174 - 1ST強帖: 19
- 應(yīng)助: 303 (大學(xué)生)
- 金幣: 11306.4
- 散金: 1123
- 紅花: 108
- 帖子: 2814
- 在線: 916.1小時
- 蟲號: 941486
- 注冊: 2010-01-12
- 專業(yè): 結(jié)構(gòu)陶瓷
- 管轄: 第一性原理
9樓2012-06-25 13:04:22
tt-0-8
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 676.8
- 帖子: 193
- 在線: 127.2小時
- 蟲號: 1315323
- 注冊: 2011-06-04
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
10樓2012-06-25 16:37:42
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 llh2010 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +9 | PolarBear11 2026-03-26 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +8 | 新時代材料 2026-03-27 | 8/400 |
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +10 | 路癡小琪 2026-03-27 | 10/500 |
|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
|
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 10/500 |
|
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
|
[考研] 總分293求調(diào)劑 +6 | 加一一九 2026-03-25 | 8/400 |
|
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。，。，。，。i 2026-03-23 | 6/300 |
|
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
|
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[考研] 277分求調(diào)劑，跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|