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

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 如何正确使用VASP计算晶体的弹性常数？
391/1返回列表
查看: 3984  |  回复: 38
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 2OLEGEND 的 34 个金币 ，回帖就立即获得 2 个金币，每人有 1 次机会

2OLEGEND

新虫 (正式写手)

[交流] 如何正确使用VASP计算晶体的弹性常数？

各位老师，本人最近在使用vasp计算petn晶体的弹性常数。步骤如下，首先优化晶体，优化的incar和kpoints如下：
incar
system=petn
ismear=0
sigma=0.05
encut=1000
ediff=0.0001
ediffg=0.001
prec=accurate
potim=0.0100
lcharg=false
lwave=false
ibrion=2
isif=3
nfree=4
nsw=100
lreal=auto

kpoints
k
0
gamma
   3   3   4
0.0  0.0  0.0

优化结束后（56步优化结束），将优化得到的contcar文件重新命名为poscar文件（计算弹性常数用）。使用下列incar作为计算，kpoints保持不变：
system=petn
ismear=0
sigma=0.05
encut=1000
ediff=0.0001
ediffg=0.001
prec=accurate
potim=0.01
lcharg=false
lwave=false
ibrion=6
isif=4
nfree=4
nsw=1
lreal=auto
为了测得更多结果，更改encut的值为：400,500,600,700,800,1000，其他的所有都保持不变，为什么计算的弹性常数差距如此大：

如何正确使用VASP计算晶体的弹性常数？
encut=400.png


如何正确使用VASP计算晶体的弹性常数？-1
encut=500.png


如何正确使用VASP计算晶体的弹性常数？-2
encut=600.png


如何正确使用VASP计算晶体的弹性常数？-3
encut=700.png


如何正确使用VASP计算晶体的弹性常数？-4
encut=800.png


如何正确使用VASP计算晶体的弹性常数？-5
encut=1000.png
回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼2019-01-22 19:14:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bobcanada

至尊木虫 (文坛精英)
★
2OLEGEND(金币+2): 谢谢参与
ding
回复此楼
9楼2019-01-23 04:58:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuqiang68

至尊木虫 (文坛精英)
★
2OLEGEND(金币+2): 谢谢参与
加油!
回复此楼
10楼2019-01-23 06:46:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhm_list

木虫 (著名写手)
★
2OLEGEND(金币+2): 谢谢参与
也学习学习

发自小木虫IOS客户端
回复此楼
13楼2019-01-23 10:55:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

漫天飘雪

版主 (知名作家)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
2OLEGEND(金币+2): 谢谢参与
2OLEGEND: 金币+10, 回复对问题有帮助 2019-01-25 12:22:07
kpoints
k
0
gamma
   3   3   4
0.0  0.0  0.0

k点是不是有点少了？
赞一下 回复此楼
15楼2019-01-23 14:40:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bobcanada

至尊木虫 (文坛精英)
ding
回复此楼
17楼2019-01-24 10:59:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2OLEGEND

新虫 (正式写手)
引用回帖:
15楼: Originally posted by 漫天飘雪 at 2019-01-23 14:40:56
kpoints
k
0
gamma
   3   3   4
0.0  0.0  0.0

k点是不是有点少了？

我又重新测试了几种K点，改成446，仍然是这样的结果，搞得我都不晓得哪个是正确的，文献中用最多的ENCUT=1000，但是他们不是计算弹性常数的。
赞一下 回复此楼
20楼2019-01-25 12:23:10
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

123dingyc

新虫 (知名作家)
★
2OLEGEND(金币+2): 谢谢参与
继续测试算

发自小木虫IOS客户端
回复此楼
23楼2019-01-26 11:48:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

598878157

版主 (文学泰斗)
★
2OLEGEND(金币+2): 谢谢参与
棵棵

发自小木虫Android客户端
回复此楼
26楼2019-01-26 18:41:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zuiaixuexi

至尊木虫 (文坛精英)
★
2OLEGEND(金币+2): 谢谢参与
祝你成功
回复此楼
29楼2019-01-28 20:20:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (著名写手)
★
2OLEGEND(金币+2): 谢谢参与
本帖仅楼主可见
赞一下
35楼2019-02-03 17:08:03
 已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

2358952545

铜虫 (正式写手)
加油，

发自小木虫Android客户端
回复此楼
36楼2019-02-04 10:52:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

千寻与猫

新虫 (小有名气)
★
2OLEGEND(金币+2): 谢谢参与
引用回帖:
20楼: Originally posted by 2OLEGEND at 2019-01-25 12:23:10
我又重新测试了几种K点，改成446，仍然是这样的结果，搞得我都不晓得哪个是正确的，文献中用最多的ENCUT=1000，但是他们不是计算弹性常数的。...

据我了解，用应变-应力的办法直接算弹性常数，对精度要求很高
赞一下 回复此楼
37楼2020-04-22 17:12:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

千寻与猫

新虫 (小有名气)
引用回帖:
20楼: Originally posted by 2OLEGEND at 2019-01-25 12:23:10
我又重新测试了几种K点，改成446，仍然是这样的结果，搞得我都不晓得哪个是正确的，文献中用最多的ENCUT=1000，但是他们不是计算弹性常数的。...

据我了解，用应变-应力的办法直接算弹性常数，对精度要求很高
赞一下 回复此楼
38楼2020-04-22 17:12:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

长风和秋雨

新虫 (初入文坛)
★
2OLEGEND(金币+2): 谢谢参与
我的也是，不同K点弹性常数都不一样

发自小木虫Android客户端
赞一下 回复此楼
39楼2023-06-08 01:18:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
syhorchid2楼
2019-01-22 19:18   回复  
2OLEGEND(金币+2): 谢谢参与
Sun_10833楼
2019-01-22 19:21   回复  
2OLEGEND(金币+2): 谢谢参与
1 发自小木虫Android客户端
飞行鸟4楼
2019-01-22 19:53   回复  
2OLEGEND(金币+2): 谢谢参与
， 发自小木虫Android客户端
youngen5楼
2019-01-22 20:07   回复  
2OLEGEND(金币+2): 谢谢参与
time886楼
2019-01-22 20:21   回复  
2OLEGEND(金币+2): 谢谢参与
tzynew7楼
2019-01-22 21:28   回复  
2OLEGEND(金币+2): 谢谢参与
nono20098楼
2019-01-22 21:47   回复  
2OLEGEND(金币+2): 谢谢参与
·
Biotage11楼
2019-01-23 09:19   回复  
2OLEGEND(金币+2): 谢谢参与
sz_kebaoyuan12楼
2019-01-23 09:24   回复  
2OLEGEND(金币+2): 谢谢参与
slytjiaofei14楼
2019-01-23 11:54   回复  
2OLEGEND(金币+2): 谢谢参与
happyever8316楼
2019-01-23 14:54   回复  
2OLEGEND(金币+2): 谢谢参与
祝福 发自小木虫IOS客户端
jiaoxg18楼
2019-01-24 16:22   回复  
2OLEGEND(金币+2): 谢谢参与
tzynew19楼
2019-01-25 12:06   回复  
胡小虾21楼
2019-01-26 10:33   回复  
2OLEGEND(金币+2): 谢谢参与
祝福
ckxbear22楼
2019-01-26 11:13   回复  
2OLEGEND(金币+2): 谢谢参与
youngen24楼
2019-01-26 12:54   回复  
xautldhldh25楼
2019-01-26 18:05   回复  
2OLEGEND(金币+2): 谢谢参与
235895254527楼
2019-01-27 09:06   回复  
2OLEGEND(金币+2): 谢谢参与
发自小木虫Android客户端
paperhunter28楼
2019-01-28 09:22   回复  
2OLEGEND(金币+2): 谢谢参与
。
youngen30楼
2019-01-28 20:27   回复  
youngen31楼
2019-01-28 20:28   回复  
xautldhldh32楼
2019-01-28 21:42   回复  
slytjiaofei33楼
2019-01-29 09:54   回复  
纳米材料75334楼
2019-01-29 17:21   回复  
2OLEGEND(金币+2): 谢谢参与
相关版块跳转 我要订阅楼主 2OLEGEND 的主题更新
391/1返回列表
提示: 如果您在30分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复(可上传附件) SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 普通院校药学硕士，做合成的，感觉找不到工作 +15 pom戴墨镜 2024-04-24 25/1250 2024-04-28 22:03 by eieie
[考博] 申博成果界定是根据Jcr分区还是中科院分区 +6 我属驴核动力驴 2024-04-22 7/350 2024-04-28 20:21 by J584326
[硕博家园] 科研好难啊，想退学去打工。 +17 byywnyl 2024-04-27 19/950 2024-04-28 17:48 by Totorop
[分析] C18色谱柱可以分析单糖，二糖，低聚糖（3-7）吗 +3 备安网公 2024-04-26 3/150 2024-04-28 17:26 by eieie
[论文投稿] with editor 两个月了，什么原因？ +9 yiersan9 2024-04-24 22/1100 2024-04-28 17:12 by ds521chun
[考博] 真的好想读博！ +17 wangzhe_bs 2024-04-22 25/1250 2024-04-28 17:10 by wangzhe_bs
[考博] 25年博士申请 +8 Changzixuan 2024-04-25 14/700 2024-04-28 16:10 by Napoleonsky
[论文投稿] 环境类微生物英文论文投稿 50+3 lyyyyds0 2024-04-26 5/250 2024-04-28 14:42 by nono2009
[有机交流] 以2氟甲苯为原料做2氟4溴甲苯 +3 高树桩 2024-04-24 7/350 2024-04-28 14:37 by 高树桩
[有机交流] 锂盐的合成 4+3 liurui1203 2024-04-27 3/150 2024-04-28 10:45 by 591950582
[教师之家] 大学直属学院卸任的副院长退休后还享受副处级的养老待遇吗？ +6 苏东坡二世 2024-04-27 6/300 2024-04-28 08:35 by ou0551
[论文投稿] AC审稿 +3 马winner 2024-04-26 7/350 2024-04-27 23:33 by topedit
[基金申请] 优青已经开始了么？ +3 majunge000 2024-04-27 3/150 2024-04-27 18:51 by spring_1987
[考博] 材料方向24博士申请/一作SCI三篇 +3 白天不碰 2024-04-24 5/250 2024-04-27 17:41 by 安塔瓦拉多
[有机交流] 环肽的合成 +3 徐来不惊 2024-04-25 6/300 2024-04-27 13:54 by xxzhangci
[基金申请] 两类问题算是白选了~ +8 jurkat.1640 2024-04-23 13/650 2024-04-27 12:03 by 淀粉搬运工
[考研] 学硕专硕 +5 小蜗牛* 2024-04-26 5/250 2024-04-26 16:43 by 鱼翔浅底1
[考博] 申博求助 +4 dskabdh 2024-04-24 11/550 2024-04-26 15:54 by dskabdh
[教师之家] 某种做法不行。说过几遍了。同学还那样做。再那样做就给低分 +4 河西夜郎 2024-04-24 4/200 2024-04-26 08:51 by Quakerbird
[考博] 博士招生 +4 zx179 2024-04-24 7/350 2024-04-24 15:01 by H考研成功
信息提示
关闭
请填处理意见
关闭