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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于BCC-Fe的弹性常数的能量-应变计算
481/1返回列表
查看: 5275  |  回复: 47
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价34次，作者dengc09增加金币 25.8 个

dengc09

禁虫 (正式写手)

[资源] 关于BCC-Fe的弹性常数的能量-应变计算

最近看到很多关于计算弹性常数的帖子，特别是BCC-Fe的例子比较典型，我在此大概算了一下，结果是a0=2.835埃，C11=262GPa，C44=111GPa，C12=158GPa。vasp计算文件和origin拟合曲线见附件。用的是侯博士的脚本，取五个点拟合，大家可以参考一下，你们的主要问题可能还是晶格常数拟合没有达到要求。当然，我的2.835也只保留到3位小数，不过基本上可以了，不说了，上附件。
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : vasp.rar
    • 2014-06-08 20:04:24, 1.1 M
  • 附件 2 : origin.rar
    • 2014-06-08 20:04:34, 74.01 K

» 收录本帖的淘帖专辑推荐

第一性原理杂 分子动力学

» 本帖已获得的红花（最新10朵）

dahui2013

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2014-06-08 20:05:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dengc09

禁虫 (正式写手)
本帖内容被屏蔽
2楼2014-06-08 20:13:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fdd096030079

新虫 (小有名气)

★ 一星级,一般

Vasp现在可以直接计算弹性常数了，不需要用侯博士的脚本。
在INCAR里设置：IBRION=6，ISIF=3，NFREE=4，即可。在OUTCAR里即可找到弹性常数矩阵。
我算的bcc_Fe的弹性常数跟楼主差不多，也补充下数据,供学习计算弹性常数的同学们参考：
C11=267Gpa,C12=145Gpa,C44=100Gpa.
另外我的晶格常数小数点后只有两位有效数字，取得是2.84A。
赞一下 回复此楼
11楼2014-06-25 08:46:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dahui2013

铜虫 (小有名气)
送红花一朵
我看到你在进行relax和static时候用的ISMEAR不同我也是和你取的一样，不过我ISIF=4，进行了离子弛豫，VASP手册上有这样一句话for the calculation of DOS and very accurate total energy calculation (no relaxation in metals) use the tetrahedron method (ISMEAR=-5),我计算的是金属体系，我这样是不是错了，还有个问题想请教下你，在做结构优化的时候是不是必须要进行离子弛豫？？希望能回复！！
赞一下 回复此楼
17楼2014-08-19 11:23:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyu0315

新虫 (小有名气)
系系诶些
回复此楼
24楼2015-04-27 21:50:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

123dingyc

新虫 (知名作家)
引用回帖:
2楼: Originally posted by dengc09 at 2014-06-08 20:13:11
不好意思，漏掉了晶格常数优化的文件，补上

有人私下问我，还是解释一句吧，根据能量-应变那个公式，计算Cij只要看拟合得到的二次项系数B2，对于文件1，B2=V0C11/2，文件2，B2=V0C44/2，文件3，B2=V0（C11+C12） ...

其他晶系的怎么拟合

发自小木虫Android客户端
赞一下 回复此楼
30楼2015-12-12 10:45:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

123dingyc

新虫 (知名作家)
引用回帖:
30楼: Originally posted by 123dingyc at 2015-12-12 10:45:43
其他晶系的怎么拟合
...

用什么公式

发自小木虫Android客户端
回复此楼
31楼2015-12-12 10:46:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

蓝小茉Q

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by dengc09 at 2014-06-08 20:13:11
不好意思，漏掉了晶格常数优化的文件，补上

有人私下问我，还是解释一句吧，根据能量-应变那个公式，计算Cij只要看拟合得到的二次项系数B2，对于文件1，B2=V0C11/2，文件2，B2=V0C44/2，文件3，B2=V0（C11+C12） ...

楼主，你的vasp 中2和3的defvector是一毛一样，输入参数啥的也是一毛一样，应变加错了吧，SUMMARY居然不一样
赞一下 回复此楼
37楼2016-06-05 07:37:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
请问楼主，我在计算cij的时候 E-S曲线的最低点不在0，是什么原因？
赞一下 回复此楼
38楼2016-08-03 19:37:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

JP1993

金虫 (小有名气)

★★★ 三星级,支持鼓励

引用回帖:
30楼: Originally posted by 123dingyc at 2015-12-12 10:45:43
其他晶系的怎么拟合
...

其他的可以根据所属晶系来判断，可以查下候老师的文献，他讲了立方晶系和六角晶系等耦合关系

发自小木虫Android客户端
赞一下 回复此楼
43楼2016-10-17 16:35:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxtp123

新虫 (初入文坛)
引用回帖:
11楼: Originally posted by fdd096030079 at 2014-06-25 08:46:41
Vasp现在可以直接计算弹性常数了，不需要用侯博士的脚本。
在INCAR里设置：IBRION=6，ISIF=3，NFREE=4，即可。在OUTCAR里即可找到弹性常数矩阵。
我算的bcc_Fe的弹性常数跟楼主差不多，也补充下数据,供学习计算弹 ...

isif=3，警惕的体积不是改变了吗，算出来的弹性常数矩阵里面并不符合立方晶系只有三个晶格常数的规律啊
赞一下 回复此楼
47楼2018-04-01 20:01:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
hw458887923楼
2014-06-09 03:46   回复  
五星好评  顶一下，感谢分享！
Quan.4楼
2014-06-09 06:14   回复  
五星好评  顶一下，感谢分享！
dalongmao5楼
2014-06-09 11:09   回复  
五星好评  
qing1086楼
2014-06-09 15:44   回复  
五星好评  顶一下，感谢分享！
小小静7楼
2014-06-09 21:30   回复  
五星好评  顶一下，感谢分享！
boneboy8楼
2014-06-10 11:02   回复  
五星好评  顶一下，感谢分享！
6759324069楼
2014-06-10 14:59   回复  
五星好评  顶一下，感谢分享！
灿烂的幸福10楼
2014-06-11 09:23   回复  
三星好评  顶一下，感谢分享！
raoqinglei12楼
2014-07-02 17:02   回复  
三星好评  顶一下，感谢分享！
tonyngch13楼
2014-07-24 09:23   回复  
一般  顶一下，感谢分享！
yansunning14楼
2014-08-10 03:30   回复  
五星好评  顶一下，感谢分享！
dahui201315楼
2014-08-19 08:53   回复  
五星好评  顶一下，感谢分享！
dahui201316楼
2014-08-19 09:05   回复  
顶一下，感谢分享！
wangzg929218楼
2014-08-26 10:51   回复  
五星好评  顶一下，感谢分享！
excellencer19楼
2014-09-04 00:07   回复  
五星好评  顶一下，感谢分享！
songbone20楼
2014-10-12 16:44   回复  
五星好评  顶一下，感谢分享！
skylve21楼
2015-01-15 20:27   回复  
五星好评  顶一下，感谢分享！
阿雏22楼
2015-04-24 08:23   回复  
五星好评  顶一下，感谢分享！
zhaoyu031523楼
2015-04-27 21:49   回复  
五星好评  顶一下，感谢分享！
iylni16325楼
2015-06-01 11:56   回复  
五星好评  顶一下，感谢分享！
qing10826楼
2015-07-11 08:26   回复  
顶一下，感谢分享！
jing043127楼
2015-10-22 20:34   回复  
五星好评  顶一下，感谢分享！
fangchaoge28楼
2015-12-03 11:20   回复  
五星好评  顶一下，感谢分享！
jimuchun29楼
2015-12-11 10:17   回复  
五星好评  顶一下，感谢分享！
酱油x酱油32楼
2015-12-20 16:29   回复  
五星好评  顶一下，感谢分享！
LoogCat33楼
2016-03-05 21:26   回复  
五星好评  顶一下，感谢分享！
begin201634楼
2016-03-24 18:04   回复  
五星好评  顶一下，感谢分享！
begin201635楼
2016-03-24 18:08   回复  
顶一下，感谢分享！
kervinzhao36楼
2016-06-01 22:00   回复  
五星好评  顶一下，感谢分享！
sbingyi39楼
2016-09-03 13:08   回复  
五星好评  顶一下，感谢分享！
江南小毛羊40楼
2016-10-02 10:59   回复  
五星好评  顶一下，感谢分享！
ECHOAJMMY41楼
2016-10-17 16:06   回复  
五星好评  顶一下，感谢分享！
ECHOAJMMY42楼
2016-10-17 16:07   回复  
顶一下，感谢分享！
但丁44楼
2017-01-16 20:55   回复  
五星好评  顶一下，感谢分享！
meizhixin45楼
2017-03-01 15:55   回复  
五星好评  顶一下，感谢分享！
meizhixin46楼
2017-03-01 15:56   回复  
顶一下，感谢分享！
zhishihaoyu48楼
2019-09-07 10:41   回复  
五星好评  顶一下，感谢分享！
相关版块跳转 我要订阅楼主 dengc09 的主题更新
481/1返回列表
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭