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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » murn.x怎么拟合含a和c两个参数的晶胞？
51/1返回列表
查看: 1512  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jinwentao

新虫 (初入文坛)

[交流] murn.x怎么拟合含a和c两个参数的晶胞？ 已有3人参与

候柱峰老师教程讲用murn.x <inp.m> out.m得到out.m。
但六方晶胞有ac两个参数，我感觉不太好写inp.m。只能得到体积与能量的EtVo.dat.
这应该怎么拟合？
回复此楼

» 猜你喜欢

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

jinwentao

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by 伤心春草 at 2014-04-10 21:12:25
不知道侯老师的表述有歧义还是那个软件不太合适，我没用那个算，用的别的eos软件拟合的~
这样得到的最优体积，再解出晶格常数。此时的晶格常数不算最优，但体积已达最优，再做一次结构迟豫以及静态计算就可以了。。 ...

我想我后来明白那个murn.x了，意思跟你的差不多。就是固定晶胞形状，a就对应了一个体积，再弛豫。反复得到几个a-E关系，然后用murn.x得到最优的体积对应的a。在POSCAR中写入那个固定形状的晶胞和a,再弛豫，就得到了这个体积下的最优的晶胞。

abc三个参数不等也没关系，先固定abc的形状然后随便取a，a对应一个体积V。把V/a^3这个比值写入inp.m（比如立方晶系就是1），murn.x可以通过这个比值由a得到V。取十几个a值计算得到inp.m，murn.x拟合得到最优a，写入a，弛豫，就可以了。
赞一下 回复此楼
5楼2014-04-16 09:23:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

zique

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
如果你用的是vasp，在用murn.x拟合得到平衡体积（例如是50）后，利用原来的POSCAR，但把第二行改为-50（注意有负号），在INCAR中设ISIF=4或5，进行一次计算。在得到的CONTCAR中有最终优化后的a, b, c晶格常数
赞一下 回复此楼
2楼2014-04-09 11:14:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伤心春草

铜虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
不知道侯老师的表述有歧义还是那个软件不太合适，我没用那个算，用的别的eos软件拟合的~
这样得到的最优体积，再解出晶格常数。此时的晶格常数不算最优，但体积已达最优，再做一次结构迟豫以及静态计算就可以了。。
不知道这个思维对不对。
我还有个问题就是，abc三个参数不等的体系应该怎么优化
赞一下 回复此楼
且喜且忧，且歌且行
3楼2014-04-10 21:12:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinwentao

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by zique at 2014-04-09 11:14:26
如果你用的是vasp，在用murn.x拟合得到平衡体积（例如是50）后，利用原来的POSCAR，但把第二行改为-50（注意有负号），在INCAR中设ISIF=4或5，进行一次计算。在得到的CONTCAR中有最终优化后的a, b, c晶格常数

POSCAR第二行可以写体积？不是晶格参数嘛？
赞一下 回复此楼
4楼2014-04-16 09:10:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭