24小时热门版块排行榜    

查看: 4757  |  回复: 25
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xh512

金虫 (正式写手)

[交流] 【求助】一定压强下的弹性常数计算方法 已有10人参与

我计算的思路是这样的:
1。改变晶胞的体积,由BM方程拟合出该体积下的压强。
2。选出拟合出来的某个压强下的晶胞体积,此体积的晶胞对应某个压强,然后按侯博士的方法,施加应变,拟合二次曲线。

可是给一定压强下的晶胞,施加应变后,能量-应变关系基本就是一条单调递减的直线,指向平衡体积的方向。能量选用的是VASP输出的能量,没有考虑压强进去。
请问各位在计算高压下的弹性常数时,是怎么操作的呢?谢谢大家了。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2):谢谢指教 2010-11-13 14:47:28
引用回帖:
Originally posted by xh512 at 2010-10-25 16:20:37:
我计算的思路是这样的:
1。改变晶胞的体积,由BM方程拟合出该体积下的压强。
2。选出拟合出来的某个压强下的晶胞体积,此体积的晶胞对应某个压强,然后按侯博士的方法,施加应变,拟合二次曲线。

可是给一定 ...

最近也一直学这个, 得到多位虫友帮助, 思路如下:
(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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 26 个回答

valenhou001

至尊木虫 (职业作家)


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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

虾米学习

新虫 (初入文坛)

求助:babaleo


小木虫(金币+0.5):给个红包,谢谢回帖交流
我最近也在算这个,分别设置PSTRESS=50kb,100kb ,150kb,用vasp5.2计算,得到的弹性常数一样,不知道怎么回事,请教一下各位高手!
还有请教侯老师,我想计算加了单轴方向的应力,计算弹性常数应该怎样实现,看了不少文献,没弄明白。
谢谢!!
5楼2010-11-15 14:47:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by 虾米学习 at 2010-11-15 14:47:31:
我最近也在算这个,分别设置PSTRESS=50kb,100kb ,150kb,用vasp5.2计算,得到的弹性常数一样,不知道怎么回事,请教一下各位高手!
还有请教侯老师,我想计算加了单轴方向的应力,计算弹性常数应该怎样实现,看 ...

我只有4.6,没用过5.2版. 帮你顶一下, 盼望高手出招.
6楼2010-11-16 09:10:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见