| 查看: 4694 | 回复: 32 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
llh2010至尊木虫 (著名写手)
|
[求助]
ATAT中gensqs工具如何改变取代浓度
|
|
| 请问各位高手们，使用ATAT软件包的gensqs工具产生sqs结构时，如何改变取代的浓度？例如对于一个bcc结构，体心上的一种原子不变，而想让立方体八个顶点上的另一种原子有一个或几个被其它原子取代，请问应该如何写lat.in和conc.in文件。 |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有24人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
tt-0-8
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 676.8
- 帖子: 193
- 在线: 127.2小时
- 虫号: 1315323
- 注册: 2011-06-04
- 专业: 金属非晶态、准晶和纳米晶
23楼2012-11-11 11:59:00
linggang87
银虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 48.4
- 散金: 60
- 红花: 1
- 帖子: 198
- 在线: 108.8小时
- 虫号: 635564
- 注册: 2008-10-24
- 专业: 金属材料的制备科学与跨学
3楼2012-06-21 23:55:49
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
llh2010
至尊木虫 (著名写手)
- 应助: 34 (小学生)
- 金币: 10007.8
- 散金: 14
- 红花: 11
- 帖子: 1614
- 在线: 350.1小时
- 虫号: 1471268
- 注册: 2011-11-01
- 性别: GG
- 专业: 原子和分子物理
5楼2012-06-22 09:40:27