| 查看: 1361 | 回复: 2 | ||
yi315730581木虫 (正式写手)
|
[求助]
vasp计算体弹性模量出问题
|
|
将文件写好之后，开始运行，summary直接宣示所设置的体积 -235 -237 -239 -241 -243 -245 -247 -249 -251 后面的数据不出来就直接计算完成。 summary的完整数据应该是 -246 3 F= -.24039430E+03 E0= -.24039430E+03 d E =-.299303E-03 -248 3 F= -.24050611E+03 E0= -.24050611E+03 d E =-.151022E-03 -250 1 F= -.24058649E+03 E0= -.24058649E+03 d E =-.240586E+03 -252 1 F= -.24063860E+03 E0= -.24063860E+03 d E =-.240639E+03 -254 1 F= -.24066267E+03 E0= -.24066267E+03 d E =-.240663E+03 -256 1 F= -.24065903E+03 E0= -.24065903E+03 d E =-.240659E+03 -258 1 F= -.24063015E+03 E0= -.24063015E+03 d E =-.240630E+03 -260 3 F= -.24057665E+03 E0= -.24057665E+03 d E =-.211674E-03 -262 3 F= -.24049992E+03 E0= -.24049992E+03 d E =-.307053E-03 问题到底出在哪里，哪位高人帮忙解释下，感激不尽  |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有62人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于TiN超晶胞掺杂的问题
已经有11人回复
- VASP结构优化出错
已经有9人回复
- 无机非金属版资源索引
已经有63人回复
- 加压后的体系怎么计算弹性常数？（或其它性质）
已经有6人回复
- 【求助】关于拟合得到体弹性模量&自旋极化&LDA+U等问题，望大家指点指点
已经有4人回复
- 【求助】非零压下弹性常数计算的初始结构
已经有9人回复
- 【求助】侯老师的vasp入门的晶格优化部分没看懂，问一下
已经有25人回复
- 【讨论】vasp结构优化之前的测试
已经有11人回复
- 【讨论】交流讨论下用VASP计算不同温度下材料的热力学性能
已经有19人回复
2楼2012-03-26 15:58:22
yi315730581
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2560.8
- 红花: 1
- 帖子: 662
- 在线: 82.4小时
- 虫号: 1011583
- 注册: 2010-05-05
- 性别: GG
- 专业: 金属结构材料
|
您好： INCAR system=TiO2-Anatase-gga ISIF=4 NSW=15 ISTART=0 IBRION=2 ISMEAR=-5 !0 ! for less K-point, set 0. try -5 (tetrahedron method valid for NKPT>=4) with multi-KPOINTS NFREE =15 NELM =120 NELMIN= 4 NELMDL = -12 ENCUT=520 NBANDS=160 Prec=Normal ROPT=-5E-4 -5E-4 -5E-4 ! * Making the calculations faster IALGO = 48 ! RMM-DIIS algorithm for electrons LREAL =Auto ! evaluate projection operators in real space NSIM = 4 ! blocked algorithm update, four bands at a time LELF=TURE ICHARG=1 K-Points 0 Monkhorst-Pack 4 4 4 0 0 0 那个summary的脚本是指run文件么 for a in -235 -237 -239 -241 -243 -245 -247 -249 -251 do echo "a= $a" cat >POSCAR < system=Sm2O3-mc $a ! opt v=243.40 0.0000000000000000 4.9314122559415781 4.9314122559415781 4.9314122559415781 0.0000000000000000 4.9314122559415781 4.9800985103559263 4.9800985103559263 -0.0486862544143462 3 5 16 Direct 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.5000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.5000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.1289041758412512 0.1289041758412512 0.1132874724762462 0.6330511441617183 0.1177105896269768 0.1161871220495867 0.1177105896269768 0.6330511441617183 0.1161871220495867 0.1189283512638254 0.1189283512638254 0.6392864472229567 0.6330511441617183 0.6330511441617183 0.1161871220495867 0.1189283512638254 0.6228568502493927 0.6392864472229567 0.6228568502493927 0.1189283512638254 0.6392864472229567 0.6235366114367832 0.6235366114367832 0.6293901656896503 0.3764633885632168 0.3764633885632168 0.3706098343103497 0.8810716487361746 0.3771431497506073 0.3607135527770435 0.3771431497506073 0.8810716487361746 0.3607135527770435 0.3669488558382817 0.3669488558382817 0.8838128779504132 0.8810716487361746 0.8810716487361746 0.3607135527770435 0.3669488558382817 0.8822894103730230 0.8838128779504132 0.8822894103730230 0.3669488558382817 0.8838128779504132 0.8710958241587488 0.8710958241587488 0.8867125275237537 ! vasp E=`tail -1 OSZICAR` echo $a $E >>SUMMARY done 我这个是一个超胞体系，而且在晶胞优化当中也会出现问题，r31能计算，但是r32不能进行下去。本人也是初步接触计算这些东西，望您不吝赐教，谢谢 |
3楼2012-03-26 20:30:18
