| 查看: 1249 | 回复: 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 问题到底出在哪里，哪位高人帮忙解释下，感激不尽  |
回复此楼» 猜你喜欢
材料调剂
已经有10人回复
-
一志愿哈尔滨工业大学085600英一数二337分求调剂
已经有5人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有11人回复
-
求调剂，一志愿南京航空航天大学 ，080500材料科学与工程学硕
已经有9人回复
-
本科211，293分请求调剂
已经有4人回复
-
求调剂
已经有3人回复
-
085400电子信息319求调剂（接受跨专业调剂）
已经有5人回复
-
268求调剂
已经有9人回复
-
一志愿双非085502，267分，过四级求调剂
已经有3人回复
-
一志愿085404，总分291，四级已过，求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于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
