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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】一定压强下的弹性常数计算方法
261/3返回列表123下一页
查看: 4871  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 312求调剂 +3 陌宸希 2026-03-16 4/200 2026-03-16 15:06 by peike
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 材料与化工（0856）304求B区调剂 +7 邱gl 2026-03-10 11/550 2026-03-14 12:18 by 邱gl
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 一志愿郑大070303，338分，求调剂 +4 dadawaf 2026-03-10 5/250 2026-03-14 01:20 by lsw010101
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 0703求调剂 +7 jtyq001 2026-03-10 7/350 2026-03-14 01:06 by JourneyLucky
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 08食品或轻工求调剂，本科发表3篇sci一区top论文，一志愿南师大食品科学与工程 +3 我是一个兵， 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[基金申请] 提交后的基金本子，已让学校撤回了，可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
信息提示
关闭
请填处理意见
关闭