24小时热门版块排行榜    

查看: 4077  |  回复: 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的回帖

2358952545

铜虫 (正式写手)


36楼2019-02-04 10:52:46
已阅   回复此楼   关注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分钟内回复过其他散金贴,则可能无法领取此贴金币
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 加上“青年”两个字,意义就变了 +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 爱一回伤一回
[基金申请] 12个本子最多让给1a2b +5 地球e村长 2024-05-13 5/250 2024-05-13 22:29 by lucky_my2010
[基金申请] 河南工业大学过年开学后,就不让用国基金报账,这合法合规吗? 100+8 知己而已 2024-05-09 34/1700 2024-05-13 21:28 by 鱼翔浅底1
[基金申请] 啥也不想干 +12 lr小机灵鬼 2024-05-10 12/600 2024-05-13 18:14 by zwp9308
[基金申请] 面上项目代表作必须有通讯作者文章吗? +22 zlxnnx111 2024-05-09 35/1750 2024-05-13 18:07 by 星火12
[考博] 上海大学2024年秋季博士招生 +10 94872300 2024-05-08 13/650 2024-05-13 17:00 by lmlm199903
[基金申请] 函评这一关首先干掉一大半,今年卷上天了 (金币+5) +27 mycaas 2024-05-07 38/1900 2024-05-13 16:32 by hhs666
[考研] 求考研导师(24年底我才考,材料类,提前联系) +4 笑笑宝公主 2024-05-12 5/250 2024-05-13 16:10 by 356169183
[考博] 24/25申博,科研助理找导师 +9 cc毕业了吗 2024-05-12 9/450 2024-05-13 14:31 by zxl_1105
[论文投稿] Journal of Electrical Engineering&Technology Reviews Completed 快一周了 10+4 qweasd12345 2024-05-12 5/250 2024-05-13 10:00 by bear2007
[考博] 准研三生态毒理学,爱化学生物新能源碳中和,电池等,25博士岗位求捞啦! +4 wangchensi 2024-05-10 14/700 2024-05-13 00:26 by 大林挺好
[基金申请] 国社科系统提交后,下载申请书,里面的思路图不显示,是什么问题 +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
[基金申请] 科研之友老是给我发消息 +6 问君611 2024-05-11 8/400 2024-05-12 17:24 by 淀粉搬运工
[有机交流] 常见的固体有机酸由哪些呢? 40+4 whw818 2024-05-09 8/400 2024-05-11 18:45 by mapenggao
[基金申请] 是这样的? +17 xiangsu121 2024-05-10 19/950 2024-05-11 11:43 by zqwcr
[硕博家园] 化斋 求解惑 +5 tyl1111111 2024-05-08 7/350 2024-05-11 08:27 by chemhua
[基金申请] 是不是面上基金获得资助最多限三项,大家怎么看 +20 bailu929 2024-05-07 21/1050 2024-05-09 20:15 by 星火12
[基金申请] 为自己的基金祈福 +7 阿呆不呆 2024-05-07 7/350 2024-05-08 17:14 by yh3271126
信息提示
请填处理意见