| 查看: 5086 | 回复: 18 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
babaleo木虫 (正式写手)
|
[交流]
【求助】算VASP弹性常数时的能量-应变出错 已有12人参与
|
|||
|
用侯老师的脚本练习算弹性常数时出错, 以C33为例, (1) 用侯老师文中INCAR, 开始时设ISIF=4 优化, 优化后的POSCAR拷成OLDPOS (2)加e3=delta, 改defvector.f 中strain(i) 并编译,以此来算对应的基矢 (3)用侯老师的脚本算应变区间-0.012 到0.012(步长0.003)的弛豫和能量, 然后拟合能量-应变曲线. 可是得到的能量-应变曲线怎么不是原点对称的呢?哪一步出错了?恳请大虾指点. 图中红色是计算结果,蓝色是拟合deltaE/V=3C33^2/2  |
回复此楼» 收录本帖的淘帖专辑推荐
 弹性常数【自用】 |
» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 弹性常数的计算
已经有6人回复
- VASP计算TiO2表面能问题 求助
已经有8人回复
- 弹性常数原子单位转换
已经有3人回复
- VASP直接计算弹性常数出错
已经有6人回复
- vasp计算应力--应变曲线时应力值发生突变的问题
已经有8人回复
- 求助：VASP单点能计算没有找到能量最低点
已经有6人回复
- 计算弹性常数时的问题
已经有10人回复
- 弹性常数的计算
已经有12人回复
- 在VASP4.6中如何利用应力应变关系求弹性常数
已经有4人回复
- VASP 计算E-V曲线求助。。
已经有10人回复
- VASP直接计算弹性常数时报错：Very bad news! internal error in subroutine IBZKPT:
已经有5人回复
- 求助VASP计算弹性常数的问题
已经有5人回复
- VASP5.2 直接计算弹性常数
已经有52人回复
- 求助 利用弹性常数计算德拜温度求解
已经有55人回复
- 弹性常数计算出错
已经有5人回复
- 计算不同晶系的弹性常数
已经有260人回复
- 【求助】castep计算弹性常数问题
已经有14人回复
- 【求助】一定压强下的弹性常数计算方法
已经有25人回复
- 【求助】怎样用VASP计算反应速率常数
已经有8人回复
- 【讨论】用vasp计算过弹性常数的请进！！更新中！！
已经有52人回复
thinkllh
铁杆木虫 (正式写手)
- 应助: 34 (小学生)
- 金币: 6569.1
- 散金: 50
- 红花: 1
- 帖子: 538
- 在线: 157小时
- 虫号: 715800
- 注册: 2009-03-05
- 性别: GG
- 专业: 极端条件下使用的金属材料
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
这个是有步骤的： 1. 计算体系E-V关系，每个体积V下先采用ISIF=4，保持体积不变，优化晶格形状和离子位置，ISMEAR取0或1（对金属），优化完之后要增加一步静态能量计算，一般ISMEAR=-5. 2. 拟合E-V关系，找到能量最小点的体积（推荐采用Gibbs2程序计算，如果不知道，请自行搜索，本版有我多个相关回复）。 3. 在最小体积下设置ISIF=4，优化晶格形状和离子位置，一遍不放心就两遍吧，最终的CONTCAR即为最稳定构型。 4. 基于上述稳定构型，采用应力-应变方法或应力-应变能方法，后者可根据侯老师教程或直接在VASP中设置IBRION=6计算，采用侯老师的方法是变形之后先优化再计算静态能，精度更高。 Note： 1. 直接采用ISIF=3得到的CONTCAR是不能用的，尤其是较复杂的体系，不信可以grep volume vasprun.xml，一般最后几步的体积都在个位数上震荡；采用上述方法得到的最小能量对应体积精度一般在小数点之后，  。2. 凡是想要得到比较精确的稳定结构的，建议都采用上述1-3步逐步施行。 3. 以上是本人失败了N次之后总结出来的，这样得到的应力-应变能曲线对称性很高，特此分享，以帮助各位少走弯路。 [ Last edited by thinkllh on 2012-6-14 at 09:22 ] |
7楼2012-06-14 09:18:56
2楼2010-08-25 22:27:05
sunyang1988
至尊木虫 (著名写手)
- 1ST强帖: 1
- 应助: 24 (小学生)
- 贵宾: 0.984
- 金币: 9397.7
- 散金: 1593
- 红花: 23
- 沙发: 1
- 帖子: 1744
- 在线: 1255.8小时
- 虫号: 965316
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2010-08-25 23:38:47
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
4楼2010-08-26 09:41:20
10