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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】一定压强下的弹性常数计算方法
261/3返回列表123下一页
查看: 4869  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂，求各位老师收留 +8 秋有木北 2026-03-14 8/400 2026-03-16 15:21 by 哦哦123
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 8/400 2026-03-16 13:42 by Dream007008
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 321求调剂 +4 大米饭！ 2026-03-15 4/200 2026-03-16 08:41 by Linda Hu
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 土木第一志愿276求调剂，科研和技能十分丰富，求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 一志愿河海大学085900土木水利专硕279求调剂不挑专业 +4 SunWwWwWw 2026-03-10 8/400 2026-03-13 02:23 by SunWwWwWw
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
信息提示
关闭
请填处理意见
关闭