24小时热门版块排行榜    

查看: 1594  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 300求调剂,材料科学英一数二 +3 leaflight 2026-03-24 3/150 2026-03-24 15:10 by 冠c哥
[考研] 293求调剂 +6 加一一九 2026-03-24 6/300 2026-03-24 14:29 by JourneyLucky
[考研] 085602 289分求调剂 +3 WWW西西弗斯 2026-03-24 3/150 2026-03-24 14:26 by JourneyLucky
[考研] 一志愿哈工大,085400,320,求调剂 +3 gdlf9999 2026-03-24 3/150 2026-03-24 10:08 by 搏击518
[考研] 资源与环境 调剂申请(333分) +6 holy J 2026-03-21 6/300 2026-03-24 09:51 by hengsmith
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗? +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 求材料,环境专业调剂 +3 18567500178 2026-03-18 3/150 2026-03-23 23:50 by 热情沙漠
[考研] 环境学硕288求调剂 +8 皮皮皮123456 2026-03-22 8/400 2026-03-23 23:47 by 热情沙漠
[考研] 材料专硕英一数二306 +8 z1z2z3879 2026-03-18 8/400 2026-03-23 20:49 by baobaoye
[考研] 263求调剂 +6 yqdszhdap- 2026-03-22 9/450 2026-03-23 12:57 by yqdszhdap-
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 289求调剂 +7 怀瑾握瑜l 2026-03-20 7/350 2026-03-22 15:57 by ColorlessPI
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 求助 +5 梦里的无言 2026-03-21 6/300 2026-03-21 17:51 by 学员8dgXkO
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
信息提示
请填处理意见