| 查看: 5014 | 回复: 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  |
回复此楼» 收录本帖的淘帖专辑推荐
 弹性常数【自用】 |
» 猜你喜欢
能否退出参与的面上项目解除限项
已经有24人回复
-
号外号外
已经有8人回复
-
filecode，4个jtjc了
已经有17人回复
-
人气不行了
已经有9人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有5人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有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
