| 查看: 1234 | 回复: 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 问题到底出在哪里，哪位高人帮忙解释下，感激不尽  |
回复此楼» 猜你喜欢
青椒八年已不青，大家都被折磨成啥样了？
已经有11人回复
-
限项规定
已经有9人回复
-
免疫学博士有名额，速联系
已经有4人回复
-
交叉科学部支持青年基金，对三无青椒是个机会吗？
已经有5人回复
-
国家基金申请书模板内插入图片不可调整大小?
已经有6人回复
-
国家级人才课题组招收2026年入学博士
已经有5人回复
-
Fe3O4@SiO2合成
已经有6人回复
-
青年基金C终止
已经有4人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于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