| 查看: 1238 | 回复: 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 问题到底出在哪里，哪位高人帮忙解释下，感激不尽  |
回复此楼» 猜你喜欢
298求调剂
已经有3人回复
-
博士自荐
已经有7人回复
-
博士推荐
已经有4人回复
-
求环氧树脂研发1名
已经有10人回复
-
280求调剂
已经有5人回复
-
什么是人一生最重要的？
已经有10人回复
-
面上可以超过30页吧？
已经有13人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有13人回复
-
版面费该交吗
已经有17人回复
-
【博士招生】太原理工大学2026化工博士
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于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