应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】算VASP弹性常数时的能量-应变出错
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）
191/2返回列表12下一页
查看: 4060  |  回复: 18
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 收录本帖的淘帖专辑推荐

弹性常数【自用】

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2010-08-25 12:23:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

thinkllh

铁杆木虫 (正式写手)
★
小木虫: 金币+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 ]
赞一下(3人) 回复此楼
Thinkwiththeworld.
7楼2012-06-14 09:18:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢建议 2010-08-26 09:02:45
建议初始优化设置ISIF=3，一般这样优化出来的晶格参数就很接近能量最小值了
赞一下(3人) 回复此楼
心之所向无惧无悔
3楼2010-08-25 23:38:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

wangle_xq

银虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢讨论 2010-08-26 09:02:36
不是原点对称，意味着没有找到能量最低点。原因么，还是晶格优化没到位。按照你现在的图，选取0.0025位置的晶格参数应该就可以了。

个人观点，欢迎讨论。
赞一下(2人) 回复此楼
网络日志 http://blog.163.com/wangle_xq/ 关注计算材料科学的应用，以及相关的计算机技术。文章以转贴为主，少量原创。欢迎参观
2楼2010-08-25 22:27:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)
谢谢讨论和建议. 看来我是需要再重新优化.
赞一下 回复此楼
4楼2010-08-26 09:41:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ufo991

铁虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
1860789楼: Originally posted by babaleo at 2010-08-25 12:23:37
用侯老师的脚本练习算弹性常数时出错, 以C33为例,
(1) 用侯老师文中INCAR, 开始时设ISIF=4 优化, 优化后的POSCAR拷成OLDPOS
(2)加e3=delta, 改defvector.f 中strain(i) 并编译,以此来算对应的基矢
(3)用侯老师的 ...

你好，请问你这个问题解决了嘛？是因为结构没优化好的原因嘛？我计算的一正交结构在有些形变的时候是对称分布，有些不是。如C11时最低能量点出现在detal=0.003时候，C22时出现在detal=-0.003时，能否请教一下？
赞一下 回复此楼
5楼2012-06-11 09:11:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)
引用回帖:
2458418楼: Originally posted by ufo991 at 2012-06-11 09:11:01
你好，请问你这个问题解决了嘛？是因为结构没优化好的原因嘛？我计算的一正交结构在有些形变的时候是对称分布，有些不是。如C11时最低能量点出现在detal=0.003时候，C22时出现在detal=-0.003时，能否请教一下？...

可能是结构没优化好的原因. 不过后来看了其他帖, 就开始并一直用应力应变方法算了.
赞一下 回复此楼
6楼2012-06-14 08:14:58
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mousekingadv

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
2458715楼: Originally posted by thinkllh at 2012-06-14 09:18:56
这个是有步骤的：
1. 计算体系E-V关系，每个体积V下先采用ISIF=4，保持体积不变，优化晶格形状和离子位置，ISMEAR取0或1（对金属），优化完之后要增加一步静态能量计算，一般ISMEAR=-5.
2. 拟合E-V关系，找到能量 ...

谢谢分享！虽然我不做应力-应变优化，还是很有帮助
赞一下 回复此楼
8楼2012-06-15 20:22:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fjfuzhou

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
2458715楼: Originally posted by thinkllh at 2012-06-14 09:18:56
这个是有步骤的：
1. 计算体系E-V关系，每个体积V下先采用ISIF=4，保持体积不变，优化晶格形状和离子位置，ISMEAR取0或1（对金属），优化完之后要增加一步静态能量计算，一般ISMEAR=-5.
2. 拟合E-V关系，找到能量 ...

你的意思就是先进性E-V计算，得到最小体积，然后在最小体积下设置ISIF=4？         -->得到contcar,然后转化成为OLDPOS?
赞一下 回复此楼
9楼2012-11-14 11:41:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
http://muchong.com/bbs/viewthread.php?tid=6328686麻烦您给解答下！！！！十分感激！！！
赞一下 回复此楼
人，要靠自己！！！
10楼2013-09-07 22:08:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 babaleo 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭