24小时热门版块排行榜    

查看: 1593  |  回复: 6

jinwentao

新虫 (初入文坛)

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

候柱峰老师教程讲用murn.x <inp.m> out.m得到out.m。
但六方晶胞有ac两个参数,我感觉不太好写inp.m。只能得到体积与能量的EtVo.dat.
这应该怎么拟合?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

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的回帖

zique

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by jinwentao at 2014-04-16 09:10:24
POSCAR第二行可以写体积?不是晶格参数嘛?...

可以,不过写体积时前面要加负号
6楼2014-04-16 10:30:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zym765319199

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
按照候柱峰的讲义,是一个一个拟合的,先a后c
7楼2016-07-27 09:13:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jinwentao 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +3 哈哈157349 2026-03-21 3/150 2026-03-24 14:11 by zhyzzh
[考研] 生物学学硕求调剂 +6 小羊睡着了? 2026-03-23 6/300 2026-03-24 13:22 by 李常安
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 298-一志愿中国农业大学-求调剂 +11 手机用户 2026-03-17 12/600 2026-03-23 23:51 by 热情沙漠
[考研] 一志愿国科过程所081700,274求调剂 +3 三水研0水立方 2026-03-23 3/150 2026-03-23 23:11 by MajorWen
[考研] 一志愿211 初试270分 求调剂 +4 谷雨上岸 2026-03-23 5/250 2026-03-23 21:18 by 不惑可乐
[考研] 352求调剂 +3 大米饭! 2026-03-22 3/150 2026-03-22 23:28 by king123!
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 313求调剂 +4 肆叁贰壹22 2026-03-19 4/200 2026-03-21 17:33 by ColorlessPI
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
信息提示
请填处理意见