24小时热门版块排行榜    

查看: 1700  |  回复: 4

维特的日记

木虫 (小有名气)

[求助] VASP 已有2人参与

大家好,本人最近在学习VASP相关的计算。看侯柱锋的那个学习手册,学习优化晶格常数。但是我发现里面一个优化六方晶格的Mg(有两个参数a,c)的例子是不是有问题的。其中他把ISIF设为了5,这样一来晶格形状就允许改变了,通过拟合BM方程得到了体积也不能计算出晶格常数a吧,因为晶格基矢也变了。

手册的做法是先固定c/a,改变a的取值,利用BM方程拟合得到平衡态的体积,进而求出a。之后固定a,改变c/a,再次得到一系列值,同样利用BM方程拟合得到c。是这样吗?但是例子使用了ISIF=5,允许晶格形状就允许改变,这样晶格基矢也变了。如何利用平衡态的体积得到a?

希望各位网友能解答一下这个问题。。。

VASP
1.png


VASP-1
2.png
回复此楼

» 猜你喜欢

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

漫天飘雪

版主 (知名作家)

认真做事,踏实做人

【答案】应助回帖

感谢参与,应助指数 +1
ISIF=3,让程序自动找mini就可以了
http://blog.sciencenet.cn/u/pfliu89
2楼2016-12-21 15:31:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

维特的日记

木虫 (小有名气)

引用回帖:
2楼: Originally posted by 漫天飘雪 at 2016-12-21 15:31:54
ISIF=3,让程序自动找mini就可以了

嗯,我知道ISIF为3的时候直接一次性优化可以的,而且大多数情况会总这种一次性优化的方法但是我还是想知道用BM方程应该如何处理。

发自小木虫IOS客户端
3楼2016-12-21 17:09:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiedeng

铜虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★
漫天飘雪: 金币+5, 谢谢交流 2016-12-28 01:14:24
我觉得是这样的:
例子中用ISIF=5,即只变晶格结构,体积和离子位置不变;这样用bash生成一系列的体积下的relax结果。找到能量最小值,定下体积,然后在此体积下取ISIF=4,找到结构和离子位置的最小值。

ISIF=3的话,我看VASP Manual里是不建议的,不如用脚本来算。
如果我做的话,第一步直接用ISIF=4来算,也许结果会更精确。不知道侯博士是不是因为计算量的原因才把离子位置先固定住的。
wwwwww
4楼2016-12-28 00:28:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

维特的日记

木虫 (小有名气)

引用回帖:
4楼: Originally posted by xiedeng at 2016-12-28 00:28:25
我觉得是这样的:
例子中用ISIF=5,即只变晶格结构,体积和离子位置不变;这样用bash生成一系列的体积下的relax结果。找到能量最小值,定下体积,然后在此体积下取ISIF=4,找到结构和离子位置的最小值。

ISIF=3 ...

谢谢您,但是ISIF为5时,晶格基矢也会变,基矢变化后你如果确定平衡体积下的长度,所以你不可能知道晶格常数

发自小木虫IOS客户端
5楼2016-12-28 09:09:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 维特的日记 的主题更新
信息提示
请填处理意见