24小时热门版块排行榜    

查看: 968  |  回复: 2

rnchen

金虫 (小有名气)

[求助] GULP 计算热膨胀系数遇到的问题

各位好,
最近开始学习gulp, 试了一下算一个钙钛矿的热膨涨系数
input 文件写的很简单,keyword 是 opti conp free dos
然后就是设置cell 和frac
在然后是势能参数,用的是bush.lib

然后计算不同温度下的lattice parameter.
发现的问题是, 自由能并不随着温度单调变化,温度高的时候, 例如大于1000K, 自由能随温度没有明显规律。
还有就是,到温度高以后,lattice parameter几乎不变了。 怀疑是不是哪里设置了一个默认的上界?因为对这个软件还不是很了解,也很少有论坛,所以才在这里问。
如果哪位大侠肯赐教,在这里先谢过了!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rnchen

金虫 (小有名气)

自己发现了原因。 温度高一点的时候, 比如超过2/3 melting point. quasiharmonic approximation 不再适用, GULP 不再适用。
温度高的时候输出文件里有警告信息 too many failed attempts to optimise.
所以optimization根本没有达到最小值

原因很大一部分是选取的势能参数
2楼2013-10-21 16:55:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uy_naiT

铜虫 (小有名气)

楼主 您好! 请问一下 gulp里input文件是怎么写得呢?在哪里有相关的命令解读文档啊(类似于lammps的manul)

发自小木虫Android客户端
3楼2018-12-09 21:22:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 rnchen 的主题更新
信息提示
请填处理意见