24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gleerat

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3926.9
散金: 111
红花: 2
帖子: 322
在线: 102.6小时
虫号: 860161
注册: 2009-09-29
性别: GG
专业: 凝聚态物性 II ：电子结构

[交流] 【求助成功】晶格优化中，murn.f的使用方法

在晶格参数优化中，拟合Murn状态方程需要用到murn.f编译生成的可执行文件，使用的输入文件为inp.m，其内容和格式为：

引用回帖:

2
0.25
6.00 7.45 50
8
6.0849 -.62120891E+01
6.2739 -.64553428E+01
6.4629 -.65246916E+01
6.6518 -.64694519E+01
6.8408 -.63284020E+01
7.0298 -.61271095E+01
7.2188 -.58843994E+01
7.4077 -.56183090E+01
-----------------------------------------

第二行是原胞与晶胞的的体积之比。
这里的晶胞是否就是在优化过程中使用的POSCAR所构建的晶胞，即如果使用原胞构建POSCAR的话，这里的比值就应该是1？

对于VASP来说，使用惯用晶胞优化得出的数据再次构建原胞而进行计算的话，对计算有什么影响么？也就是说使用原胞还是使用惯用晶胞优化对结果有无影响？

在静态自洽计算中使用的kmesh的大小与非自洽中的密集kmesh的大小有什么关系，也就是说满足条件的非自洽计算的kmesh至少要比静态计算的kmesh大多少？

[ Last edited by gleerat on 2010-11-13 at 13:01 ]

回复此楼

» 猜你喜欢

有没有人能给点建议已经有5人回复
假如你的研究生提出不合理要求已经有12人回复
实验室接单子已经有7人回复
全日制（定向）博士已经有5人回复
萌生出自己或许不适合搞科研的想法，现在跑or等等看？已经有4人回复
Materials Today Chemistry审稿周期已经有4人回复
参与限项已经有3人回复
对氯苯硼酸纯化已经有3人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复

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

环境样品地下水中细菌总DNA的提取方法已经有23人回复
水分子建模方法已经有6人回复
**有图有真相**求助-日立F-4500荧光光谱仪--模式选择？？？已经有20人回复
晶界处的氧气和晶格中的氧空位如何表征已经有8人回复
关于硼氢化钠的使用方法已经有18人回复
在做CASTP模拟时，建立晶格需要晶格常数，如何获取呢？已经有10人回复
利用什么方法可以确定掺杂元素究竟进入晶格里面没有？已经有18人回复
CVD方法制graphene，选substrate有什么标准啊？为啥选的是Cu，Ni。与晶格有什么关系？已经有9人回复
急！在vasp中使用murn.f程序拟合晶格常数的问题已经有18人回复
求问用X-ray diffraction 能定量的测定出晶体的晶格参数吗？已经有5人回复
晶格常数优化已经有7人回复
【讨论】使用药典上的分析方法，是否还需要做方法验证？已经有20人回复
【求助】用murn.f来拟合晶格常数已经有4人回复
【原创】晶格优化中ISIF取值的探究（六方MgO为例）已经有9人回复
【求助】si晶格常数已经有7人回复

明察、慎思、笃行

1楼 2010-08-24 17:45:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

gleerat

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3926.9
散金: 111
红花: 2
帖子: 322
在线: 102.6小时
虫号: 860161
注册: 2009-09-29
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

Originally posted by 锐利的碎片 at 2010-08-25 20:25:30:

最低点体积在param.out里面

在候柱峰的《入门指南》中说道：

引用回帖:

采用状态方程拟合得到平衡状态下的体积、体弹模量。
在该体积下，重复上面第二步和第三步，得到平衡状态下的其他的晶胞参数。这一步也就是：在得到了E(V)曲线以后，通过状态方程拟合得到平衡状态下的体积，计算出上面脚本中变量$i的值，并改变$i的循环值，再运行runcell计算一次，得到其他的结构参数c和位置u。

1、这里说到“计算出上面脚本中变量$i的值”，因为极有可能平衡体积不是计算中给定的$i值对应的体积，而晶胞形状也改变了，那么如何计算平衡状态下的$i值？
2、“改变$i的循环值”，是说“原来优化的是a，现在改为优化b”这个意思么？
3、在EVPAI.OUT和EVPAP.OUT中两列数据是对应着么？怎么同样是能量，在EVPAP.OUT第一列数据（能量）明显要比EVPAI.OUT中的（能量）大的很多，而且EVPAI.OUT的两列数据（能量和体积）明显不在EVPAP.OUT范围内？在PARAM.OUT中出现的平衡体积也不在EVPAI.OUT中出现？

引用回帖:

CaFe2As2

Murnaghan EOS
Murnaghan F D, Am. J. Math. 49, p235 (1937)

(Default units are atomic: Hartree, Bohr etc.)

V0             =          0.8581027941E-05
E0             =          -155.8020874
B0             =          0.3000079390E-02
B0'             =          0.1000005891E-04

B0 (GPa)       =          88.26536818

[ Last edited by gleerat on 2010-8-25 at 21:58 ]