| 查看: 3705 | 回复: 15 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】求证晶格常数的优化 已有6人参与
|
|||
|
各位虫友，你们好！ 近日看到侯老师关于优化晶格常数的INCAR文件（Page22）如下： SYSTEM=Al-fcc ENCUT=250 ISTAR=0；ICAHRG=2 ISMEAR=-5 PREC=Accurate 因为NSW默认为0, IBRION为-1(原子不动)，ISIF为0 按照这样我做了个Pt的晶格常数计算，EtVo.dat文件如下： 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 19.8800 -6.934469 这是怎么搞的？请各位帮忙解答下，谢谢！ run_a0脚本如下，与侯老师给的几乎一样： #!/bin/sh rm WAVECAR for i in 3.5 3.6 3.7 3.8 3.9 4.0 4.1 4.2 4.3 do cat > POSCAR < fcc $i 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 1 Direct 0.0 0.0 0.0 ! echo "a=$i"; vaspall E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f\n", $5}'` V=`grep "volume" OUTCAR | tail -1 | awk '{printf "%12.4f\n", $5}'` echo $V $E >> EtVo.dat done [ Last edited by lugang911 on 2010-6-21 at 15:41 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
考博求助
已经有3人回复
lvjian8596
木虫 (小有名气)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 0.014
- 金币: 3531.7
- 红花: 2
- 帖子: 296
- 在线: 209.1小时
- 虫号: 545009
- 注册: 2008-04-13
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
15楼2010-07-12 14:32:27
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
4楼2010-06-21 16:03:14
6楼2010-06-21 16:27:24
condensed
木虫 (著名写手)
- 应助: 21 (小学生)
- 金币: 3223.5
- 散金: 44
- 红花: 3
- 帖子: 1002
- 在线: 208.4小时
- 虫号: 292567
- 注册: 2006-11-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+1):谢谢回复 2010-06-21 20:09:45
gavinliu7390(金币+2):非常感谢详细指导！ 2010-06-24 17:38:14
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+1):谢谢回复 2010-06-21 20:09:45
gavinliu7390(金币+2):非常感谢详细指导！ 2010-06-24 17:38:14
|
ISMEAR=-5 你用这个结构优化是不太好的. -5是tetrahedral的取点方式,有点是计算总能准但是计算原子受力的时候不准.最好用1来做.而且NSW=0只会做一步SCF,根本不能达到优化的目的. 一般来说我自己会用 NSW = 50 IBION = 2 ISIF = 3 用的CG的优化方式同时优化所有的晶格参数. [ Last edited by condensed on 2010-6-21 at 16:51 ] |
7楼2010-06-21 16:48:23
