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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】一定压强下的弹性常数计算方法
261/3返回列表123下一页
查看: 4870  |  回复: 25
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xh512

金虫 (正式写手)

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

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

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

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

vasp压力

» 猜你喜欢

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

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的回帖

sg18408926

至尊木虫 (著名写手)
不错，顶一个
回复此楼
4楼2010-11-15 11:57:34
已阅   回复此楼   关注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的回帖

babaleo

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

我只有4.6,没用过5.2版. 帮你顶一下, 盼望高手出招.
赞一下(1人) 回复此楼
6楼2010-11-16 09:10:59
已阅   回复此楼   关注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的回帖

虾米学习

新虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
谢谢，我再试试，还有其他意见么？
期盼中。。。。
赞一下(1人) 回复此楼
8楼2010-11-17 23:15:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengshiquan

铁杆木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by babaleo at 2010-11-13 11:59:00:
最近也一直学这个, 得到多位虫友帮助, 思路如下:
(1)高精度ISIF=3优化结构,压强通过PSTRESS设定
(2)根据材料的对称性, 选定几个应变,比如 (100000), 等,每个选定应变方式下取6-8个不同幅度应变值优化原子位置 ...

不知道虫友是如何设定几个应变,比如 (100000)的，需要在那个文件中设置
另外你说的对6-8个不同幅度应变值下的得到的不同 F-Cij 进行拟和 就可以得到Cij ，用什么拟合，是BM方程吗
赞一下(1人) 回复此楼
9楼2011-03-31 10:14:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 谢谢指点 2011-04-01 08:22:44
引用回帖:
Originally posted by fengshiquan at 2011-03-31 10:14:22:
不知道虫友是如何设定几个应变,比如 (100000)的，需要在那个文件中设置
另外你说的对6-8个不同幅度应变值下的得到的不同 F-Cij 进行拟和 就可以得到Cij ，用什么拟合，是BM方程吗

搜一下侯博士的指南,(100000) 是来说明应变(e1,e2,e3,e4,e5,e6)里的e1不为零, 其他五个都是零, e1取几个值算一下,就可以stresss=cij*strain 拟合了. 文献里看到别人有的就算e=-0.01 和0.01下的stress, 然后直线拟合.
至于加那几个,  . 可以根据对称性和stress(i)=cij*stress(j)随你自己选.比如说cubic, 只用(1001000)就可以了.
赞一下(2人) 回复此楼
10楼2011-04-01 07:29:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xh512 的主题更新
261/3返回列表123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 306求调剂 +4 唐薏薏 2026-03-09 4/200 2026-03-14 01:19 by JourneyLucky
[考研] 0703求调剂 +7 jtyq001 2026-03-10 7/350 2026-03-14 01:06 by JourneyLucky
[考研] 材料工程专硕，一志愿中国矿业大学，总分314，求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
[考研] 一志愿：武汉理工，材料工程，英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
关闭
请填处理意见
关闭