| 查看: 5010 | 回复: 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  |
回复此楼» 收录本帖的淘帖专辑推荐
 弹性常数【自用】 |
» 猜你喜欢
能否退出参与的面上项目解除限项
已经有14人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
-
放榜前的不淡定
已经有17人回复
-
科研孤儿太难了
已经有19人回复
-
明天应该可查了！？
已经有4人回复
-
什么时候开奖？
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 弹性常数的计算
已经有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
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
2楼2010-08-25 22:27:05
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
4楼2010-08-26 09:41:20
ufo991
铁虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 金币: 1148.3
- 散金: 14
- 红花: 1
- 帖子: 358
- 在线: 104.6小时
- 虫号: 671678
- 注册: 2008-12-10
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
5楼2012-06-11 09:11:01
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
6楼2012-06-14 08:14:58
mousekingadv
金虫 (小有名气)
- 应助: 51 (初中生)
- 金币: 1141.8
- 红花: 1
- 帖子: 140
- 在线: 28.7小时
- 虫号: 1721002
- 注册: 2012-03-27
- 专业: 理论和计算化学
8楼2012-06-15 20:22:09
9楼2012-11-14 11:41:44
tonglijia
木虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 1150.6
- 散金: 4
- 红花: 13
- 帖子: 519
- 在线: 305.9小时
- 虫号: 1776353
- 注册: 2012-04-24
- 性别: GG
- 专业: 金属材料的制备科学与跨学
10楼2013-09-07 22:08:23
20