| 查看: 1397 | 回复: 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 问题到底出在哪里，哪位高人帮忙解释下，感激不尽  |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有12人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有20人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有16人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于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
