| 查看: 4768 | 回复: 25 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】一定压强下的弹性常数计算方法 已有10人参与
|
|||
|
我计算的思路是这样的： 1。改变晶胞的体积，由BM方程拟合出该体积下的压强。 2。选出拟合出来的某个压强下的晶胞体积，此体积的晶胞对应某个压强，然后按侯博士的方法，施加应变，拟合二次曲线。 可是给一定压强下的晶胞，施加应变后，能量-应变关系基本就是一条单调递减的直线，指向平衡体积的方向。能量选用的是VASP输出的能量，没有考虑压强进去。 请问各位在计算高压下的弹性常数时，是怎么操作的呢？谢谢大家了。 |
回复此楼» 猜你喜欢
求个博导看看
已经有17人回复
-
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
gzqdyouxia(金币+2): 谢谢提示 2011-04-01 12:01:57
小木虫(金币+0.5):给个红包，谢谢回帖交流
gzqdyouxia(金币+2): 谢谢提示 2011-04-01 12:01:57
|
加应变(e1,e2,e3,e4,e5,e6)通过 POSCAR来表示, 就是侯博士指南里的”defvector.f”那个程序. 不加应变优化后的CONTCAR拷成OLDPOS, 加注原子数, 执行那个程序可以生成不同应变后的POSCAR, 然后ISIF=2优化后再静态算, 就得到能量和应力了. 这一步不管是用应变-能量法, 或者应力-应变法都是一样的. 比如Orthorombic, 有9个非零 Cij, 加应变(00e00e)的话, s1 C11 C12 C13 0 0 0 0 s2 C12 C22 C23 0 0 0 0 s3 = C13 C23 C33 0 0 0 X e s4 0 0 0 C44 0 0 0 s5 0 0 0 0 C55 0 0 s6 0 0 0 0 0 C66 e 对应力s1, 从上面矩阵s1=c13xe, 加不同应变, 算出不同应力, 比如, strain magnitute1=-0.006, Stress s1== -30.775258 strain magnitute2=-0.0036, stress s1=-30.376274 strain magnitute3=-0.0012, Stress s1=-29.925712 strain magnitute4=0.0012, Stress s1=-29.504233 strain magnitute5=0.0036, Stress s1=-29.063 strain magnitute6=0.006, Stress s1=-28.6 对上面应力-应变值 拟和, 拟和直线斜率就是c13班(c13=181.4GPa).同理,用s2, s3, s6可得到 c23,c33, c66. |
12楼2011-04-01 11:35:40
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
★
xh512(金币+10):谢谢侯老师 2010-10-26 21:11:17
zzy870720z(金币+1):谢谢指教 2010-11-15 12:05:50
xh512(金币+10):谢谢侯老师 2010-10-26 21:11:17
zzy870720z(金币+1):谢谢指教 2010-11-15 12:05:50
|
参考这篇文献： J. Phys.: Condens. Matter 21 (2009) 115401 http://iopscience.iop.org/0953-8984/21/11/115401/ |
2楼2010-10-26 20:56:16
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢指教 2010-11-13 14:47:28
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2):谢谢指教 2010-11-13 14:47:28
|
最近也一直学这个, 得到多位虫友帮助, 思路如下: (1)高精度ISIF=3优化结构,压强通过PSTRESS设定 (2)根据材料的对称性, 选定几个应变,比如 (100000), 等,每个选定应变方式下取6-8个不同幅度应变值优化原子位置然后算静态能(参考侯博士教程) (3)从静态能结果OUTCAR里找出各应变幅度对应的应力值,根据材料对称性, 决定每个应力对应的应变, 比如说(1000000) 中F1=C11*应变, 等等. $>grep "in kB" OUTCAR 会显示 F1, F2, F3, F4, F5 F6 (单位kbar ) (4)对6-8个不同幅度应变值下的得到的不同 F-Cij 进行拟和 就可以得到Cij |
3楼2010-11-13 11:59:00
5楼2010-11-15 14:47:31