| 查看: 3341 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP |
» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有123人回复
-
今天基金会出结果吗？20260819
已经有13人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于优化晶格常数
已经有4人回复
- 用LDA赝势优化出来的晶格常数一定比实验值小吗？谢谢！
已经有10人回复
- 晶格常数优化
已经有7人回复
- 【求助】请教手动进行晶格常数的优化
已经有4人回复
- 【求助】新虫求助！哪位前辈指点一下用vasp做空位缺陷的流程
已经有7人回复
- 【求助】用CASTEP计算晶体时结构问题
已经有10人回复
- 【求助】晶格常数
已经有6人回复
- 【求助】用murn.f来拟合晶格常数
已经有4人回复
- 【求助】镍110表面优化失败的原因在那？
已经有23人回复
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
2楼2010-06-21 15:53:27
3楼2010-06-21 15:54:37
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
5楼2010-06-21 16:10:10
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
锐利的碎片
木虫 (正式写手)
star watcher
- 1ST强帖: 1
- 应助: 136 (高中生)
- 金币: 3637.1
- 散金: 252
- 红花: 22
- 帖子: 988
- 在线: 1224.9小时
- 虫号: 961933
- 注册: 2010-03-05
- 专业: 凝聚态物性 II ：电子结构
8楼2010-06-21 16:53:10
condensed
木虫 (著名写手)
- 应助: 21 (小学生)
- 金币: 3223.5
- 散金: 44
- 红花: 3
- 帖子: 1002
- 在线: 208.4小时
- 虫号: 292567
- 注册: 2006-11-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2010-06-21 16:54:59
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
10楼2010-06-21 20:10:58
