24小时热门版块排行榜    

查看: 5881  |  回复: 47
【奖励】 本帖被评价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朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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的回帖
简单回复
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   回复  
五星好评  顶一下,感谢分享!
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   回复  
五星好评  顶一下,感谢分享!
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   回复  
五星好评  顶一下,感谢分享!
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   回复  
五星好评  顶一下,感谢分享!
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   回复  
顶一下,感谢分享!
2019-09-07 10:41   回复  
五星好评  顶一下,感谢分享!
相关版块跳转 我要订阅楼主 dengc09 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +37 医学老男孩 2026-08-13 86/4300 2026-08-18 00:18 by longsantaizi
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
请填处理意见