| 查看: 4068 | 回复: 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  |
回复此楼» 收录本帖的淘帖专辑推荐
 弹性常数【自用】 |
» 猜你喜欢
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
孩子确诊有中度注意力缺陷
已经有12人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
弹性常数的计算
已经有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计算过弹性常数的请进！！更新中！！
已经有51人回复
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