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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】一定压强下的弹性常数计算方法
51/1返回列表
查看: 5447  |  回复: 25
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xh512

金虫 (正式写手)

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

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

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

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

vasp压力

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-10-25 16:20:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

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

不会是PSTRESS设置问题吧, 重新执行 "PSTRESS=50" 和 "PSTRESS=100" "PSTRESS=150", 然后检查三个压力下的优化后的晶格常数, 如果晶格常数有变化, , 弹性常数也应该变吧.
赞一下(2人) 回复此楼
7楼2010-11-16 09:45:59
已阅   回复此楼   关注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/
赞一下(1人) 回复此楼
2楼2010-10-26 20:56:16
已阅   回复此楼   关注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
赞一下(2人) 回复此楼
3楼2010-11-13 11:59:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

虾米学习

新虫 (初入文坛)

求助：babaleo

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
我最近也在算这个，分别设置PSTRESS=50kb，100kb ，150kb，用vasp5.2计算，得到的弹性常数一样，不知道怎么回事，请教一下各位高手！
还有请教侯老师，我想计算加了单轴方向的应力，计算弹性常数应该怎样实现，看了不少文献，没弄明白。
谢谢！！
赞一下(1人) 回复此楼
5楼2010-11-15 14:47:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 26 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！ +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[有机交流] 同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？ +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
关闭
请填处理意见
关闭