24小时热门版块排行榜    

查看: 4068  |  回复: 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

新虫 (正式写手)


引用回帖:
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的回帖
相关版块跳转 我要订阅楼主 2OLEGEND 的主题更新
提示: 如果您在30分钟内回复过其他散金贴,则可能无法领取此贴金币
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 24/25申博,科研助理找导师 +6 cc毕业了吗 2024-05-12 6/300 2024-05-12 21:36 by 申博小能手
[基金申请] NSFC函评几月几日截止提交评审意见? +9 瞬息宇宙 2024-05-11 9/450 2024-05-12 19:53 by 20081002
[硕博家园] 求助,请问我是否该转专业 +6 ygcbl 2024-05-12 8/400 2024-05-12 19:44 by 考研er
[基金申请] 十个本子来自一个高校,六个本子来自一个导师的学生 +20 babu2015 2024-05-10 21/1050 2024-05-12 18:23 by saint0link
[基金申请] 基金委治打招呼顽疾越治越严重 +28 zzahkj 2024-05-10 45/2250 2024-05-12 16:42 by wolfgangHugh
[有机交流] 长烷基链氰基水解成羧酸 求助 5+5 Sheerio 2024-05-05 11/550 2024-05-12 16:11 by zyp0009928
[考博] 准研三生态毒理学,爱化学生物新能源碳中和,电池等,25博士岗位求捞啦! +3 wangchensi 2024-05-10 13/650 2024-05-12 15:57 by wangchensi
[基金申请] 基金函评会不会跨学部的送啊 +7 cls512 2024-05-09 7/350 2024-05-12 15:36 by 一路向东
[教师之家] 学生家长私下联系老师修改成绩不成,唆使19名学生联名要求复核成绩 +16 sjtu2012 2024-05-11 18/900 2024-05-12 15:32 by 会飞的猪157
[找工作] 985博士毕业,前几天通过了河南省科学院下面一个所的面试,值得去吗? +12 wendao_2016 2024-05-07 18/900 2024-05-12 14:26 by 安小樱
[考博] 矿大本-中科大硕(推免)-24年科研助理求职-方向不限 +6 Kelaizhang 2024-05-08 8/400 2024-05-12 12:50 by Kelaizhang
[基金申请] 啥也不想干 +10 lr小机灵鬼 2024-05-10 10/500 2024-05-12 09:11 by tfang
[基金申请] 面上项目代表作必须有通讯作者文章吗? +21 zlxnnx111 2024-05-09 33/1650 2024-05-12 00:40 by 星火12
[催化] 齿球形催化剂的尺寸 +3 anndy1971 2024-05-08 5/250 2024-05-11 23:57 by 596699273
[有机交流] 常见的固体有机酸由哪些呢? 40+4 whw818 2024-05-09 8/400 2024-05-11 18:45 by mapenggao
[论文投稿] 期刊 International Journal of Fuzzy Systems 稿费问题求助 20+4 一眼灬清泉 2024-05-09 7/350 2024-05-11 10:15 by allen123412
[基金申请] 信息口青基送审了吗 +5 WeichaoDing 2024-05-08 7/350 2024-05-10 19:05 by pinyuan33
[论文投稿] 论文投稿这种情况希望大吗? +3 woshichuxuan 2024-05-08 6/300 2024-05-09 17:27 by topedit
[考博] 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
信息提示
请填处理意见