| 查看: 4872 | 回复: 32 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
llh2010至尊木虫 (著名写手)
|
[求助]
ATAT中gensqs工具如何改变取代浓度
|
||||
| 请问各位高手们，使用ATAT软件包的gensqs工具产生sqs结构时，如何改变取代的浓度？例如对于一个bcc结构，体心上的一种原子不变，而想让立方体八个顶点上的另一种原子有一个或几个被其它原子取代，请问应该如何写lat.in和conc.in文件。 |
回复此楼» 收录本帖的淘帖专辑推荐
 学习 | VASP | 我的分享 |
» 猜你喜欢
323求调剂
已经有6人回复
-
一志愿北京化工大学 070300 学硕 336分 求调剂
已经有4人回复
-
352求调剂
已经有3人回复
-
一志愿东华大学化学070300，求调剂
已经有8人回复
-
277材料科学与工程080500求调剂
已经有7人回复
-
317求调剂
已经有18人回复
-
293求调剂
已经有5人回复
-
280分求调剂 一志愿085802
已经有7人回复
-
0854电子信息求调剂
已经有3人回复
-
263求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 金刚石工具的切削原理？
已经有7人回复
- 用过冠醚的朋友们, 你们一般往有机溶剂中溶盐都用多大浓度啊?
已经有3人回复
- 表面活性剂F127和P123在水中，随浓度变化，胶束形态改变的具体数据。
已经有14人回复
- 血红细胞的分离与固定所用的PBS浓度以及固定液
已经有5人回复
- 旋光仪计算氨基酸及对映体的浓度
已经有3人回复
- 有谁知道二硫化硒SeS2的浓度积常数Ksp值
已经有13人回复
- 【求助】树脂分子在乙醇中随着浓度的降低而变浑析出的反常现象解释
已经有3人回复
linggang87
银虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 48.4
- 散金: 60
- 红花: 1
- 帖子: 198
- 在线: 108.8小时
- 虫号: 635564
- 注册: 2008-10-24
- 专业: 金属材料的制备科学与跨学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★
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产生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
linggang87
银虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 48.4
- 散金: 60
- 红花: 1
- 帖子: 198
- 在线: 108.8小时
- 虫号: 635564
- 注册: 2008-10-24
- 专业: 金属材料的制备科学与跨学
3楼2012-06-21 23:55:49
llh2010
至尊木虫 (著名写手)
- 应助: 34 (小学生)
- 金币: 10007.8
- 散金: 14
- 红花: 11
- 帖子: 1614
- 在线: 350.1小时
- 虫号: 1471268
- 注册: 2011-11-01
- 性别: GG
- 专业: 原子和分子物理
5楼2012-06-22 09:40:27
llh2010
至尊木虫 (著名写手)
- 应助: 34 (小学生)
- 金币: 10007.8
- 散金: 14
- 红花: 11
- 帖子: 1614
- 在线: 350.1小时
- 虫号: 1471268
- 注册: 2011-11-01
- 性别: GG
- 专业: 原子和分子物理
6楼2012-06-22 10:01:05
