24小时热门版块排行榜    

查看: 4078  |  回复: 38
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 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
回复此楼

» 猜你喜欢

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

查看全部散金贴

已阅   回复此楼   关注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的回帖
查看全部 39 个回答

zhm_list

木虫 (著名写手)



2OLEGEND(金币+2): 谢谢参与
也学习学习

发自小木虫IOS客户端
13楼2019-01-23 10:55:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
Sun_10833楼
2019-01-22 19:21   回复  
2OLEGEND(金币+2): 谢谢参与
1 发自小木虫Android客户端
飞行鸟4楼
2019-01-22 19:53   回复  
2OLEGEND(金币+2): 谢谢参与
发自小木虫Android客户端
提示: 如果您在30分钟内回复过其他散金贴,则可能无法领取此贴金币
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 第一次投SCI,一审给了Revise +4 慎独的小花卷 2024-05-13 8/400 2024-05-14 01:21 by yifeng11
[教师之家] 学生家长私下联系老师修改成绩不成,唆使19名学生联名要求复核成绩 +20 sjtu2012 2024-05-11 22/1100 2024-05-13 23:44 by fysh2023
[教师之家] 加上“青年”两个字,意义就变了 +6 zylfront 2024-05-13 7/350 2024-05-13 23:01 by peterrjp
[有机交流] 爬大板 +4 abcde! 2024-05-12 8/400 2024-05-13 22:58 by 爱一回伤一回
[硕博家园] 哈工大硕博招生! +4 nailooo 2024-05-12 5/250 2024-05-13 22:04 by yuanjijoy
[教师之家] 南京林业大学老师首聘期考核不过,服药自杀 +11 babu2015 2024-05-12 11/550 2024-05-13 16:27 by jjsfb
[考研] 求考研导师(24年底我才考,材料类,提前联系) +4 笑笑宝公主 2024-05-12 5/250 2024-05-13 16:10 by 356169183
[教师之家] 相比过去,现在高校导师水平越来越高(可见招聘条件),研究生越来越差 +10 苏东坡二世 2024-05-11 14/700 2024-05-13 10:07 by yanglp
[教师之家] 请问有没有想要联合培养学生的老师? +5 吴小懒 2024-05-08 6/300 2024-05-13 09:22 by Xiaolin81
[考博] 准研三生态毒理学,爱化学生物新能源碳中和,电池等,25博士岗位求捞啦! +4 wangchensi 2024-05-10 14/700 2024-05-13 00:26 by 大林挺好
[基金申请] F06送审了吗? +4 昂首走 2024-05-11 7/350 2024-05-13 00:07 by blueearth171
[基金申请] 国社科系统提交后,下载申请书,里面的思路图不显示,是什么问题 +3 十三画福将 2024-05-09 3/150 2024-05-12 23:01 by 老虎当猫养
[硕博家园] 求助,请问我是否该转专业 +7 ygcbl 2024-05-12 9/450 2024-05-12 22:22 by tfang
[催化] 齿球形催化剂的尺寸 +3 anndy1971 2024-05-08 5/250 2024-05-11 23:57 by 596699273
[基金申请] 是这样的? +17 xiangsu121 2024-05-10 19/950 2024-05-11 11:43 by zqwcr
[硕博家园] 科研不是打打杀杀,科研就是人情世故 +10 杞天大圣 2024-05-08 10/500 2024-05-11 09:05 by 贪吃fish
[硕博家园] 化斋 求解惑 +5 tyl1111111 2024-05-08 7/350 2024-05-11 08:27 by chemhua
[电化学] 硫化锂正极TEM测试如何制样? +3 发仔子呀 2024-05-07 3/150 2024-05-10 13:14 by 姜尚真。
[考博] 25光催化申博 +3 wjtab 2024-05-08 5/250 2024-05-09 08:26 by 安塔瓦拉多
[基金申请] 想问下,F口信息的函评是不是还没开始? +3 wowhow 2024-05-07 5/250 2024-05-08 12:31 by cuicyge
信息提示
请填处理意见