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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【原创】晶格优化中ISIF取值的探究（六方MgO为例）
51/1返回列表
查看: 3096  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

奥古拉斯

无虫 (小有名气)

[交流] 【原创】晶格优化中ISIF取值的探究（六方MgO为例）已有6人参与

一、目的：搞清楚用vasp进行晶格优化中ISIF值的重要性
二、方法：采用分步优化和一次性优化两种方法探究
三、重点参数说明：
       ISIF：2－>relax ions?/Y; || change cell shape?/N; || change cell volume?/N
                 3－>relax ions?/Y; || change cell shape?/Y; || change cell volume?/Y
                 4－>relax ions?/Y; || change cell shape?/Y; || change cell volume?/N
                 5－>relax ions?/N; || change cell shape?/Y; || change cell volume?/N
                 6－>relax ions?/N; || change cell shape?/Y; || change cell volume?/Y
                 7－>relax ions?/N; || change cell shape?/N; || change cell volume?/Y
     POSCAR：
a
0.0              -1.0     0.0
0.8660254037844   0.5     0.0
0.0               0.0    1.6230529595
2
Direct
0.6666666666666667    0.3333333333333333 0.750
0.3333333333333333    0.6666666666666667 0.250

四、实验步骤：
A、采用分步优化法：
      算法原理：
      a.先进行一次驰豫优化计算
      b.将CONTCAR拷贝成POSCAR
      c.再进行一次静态计算，求出总能和体积，拟合状态方程，得出晶格常 数

     I、取ISIF=2时：（＃上下对应。上行为系列缩放系数；下行为相应的总能；所取CONTAR＆volume为缩放系数取3.71时的结构情况）
  2.91              3.01               3.11              3.21              3.31              3.41              3.51              3.61 3.71   -0.09768706 -2.60381322  -2.62609278 -2.62612787 -2.62612791 -2.60064488  -2.59725747 -2.59725252 -2.597253
volume of cell :       71.78
CONTCAR：
3.710000000000000
     0.0000000000000000   -1.0000000000000000    0.0000000000000000
     0.8660254037844000    0.5000000000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.6230529595000000
                              
    II、取ISIF＝3时：
2.91              3.01               3.11              3.21              3.31              3.41              3.51              3.61 3.71
-1.47880875 -3.56749270  -3.57829654 -3.57830700 -3.57830701 -3.55534231 -3.55344206 -3.55344601 -3.553446
volume of cell :       43.12
CONTCAR：
3.710000000000000
     0.0000000000000000   -0.8480309609621116    0.0000000000000000
     0.7344163553888973    0.4240154804810558    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.3556928126902820
                                       
    III、取ISIF＝4时：
2.91              3.01               3.11              3.21              3.31              3.41              3.51              3.61 3.71
-0.24918897 -2.60258815 -2.63108804 -2.63118786 -2.63118806 -2.60762910 -2.60458644 -2.60459724 -2.604597
volume of cell :       71.78        
CONTCAR：
3.710000000000000
     0.0000000000000001   -1.0553127846071310    0.0000000000000000
     0.9139276804085778    0.5276563923035652    0.0000000000000000
     0.0000000000000001    0.0000000000000000    1.4573715661012910

  IV、取ISIF＝5时：
2.91              3.01               3.11              3.21              3.31              3.41              3.51              3.61 3.71
-0.24918897 -2.60258815 -2.63108804 -2.63118786 -2.63118806 -2.60762910 -2.60458644 -2.60459724 -2.604597
volume of cell :       71.78
CONTCAR：
3.710000000000000
     0.0000000000000001   -1.0553127846071310    0.0000000000000000
     0.9139276804085778    0.5276563923035652    0.0000000000000000
     0.0000000000000001    0.0000000000000000    1.4573715661012910

V、取ISIF＝6时：
2.91              3.01               3.11              3.21              3.31              3.41              3.51              3.61 3.71
-1.47880875 -3.56749270 -3.57829654 -3.57830700 -3.57830701 -3.55534231 -3.55344206 -3.55344601 -3.553446
volume of cell :       43.12
CONTCAR:
3.710000000000000
     0.0000000000000000   -0.8480309609621116    0.0000000000000000
     0.7344163553888973    0.4240154804810558    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.3556928126902820

VI、取ISIF＝7时：
2.91              3.01               3.11              3.21              3.31              3.41              3.51              3.61 3.71
-1.28296288 -3.56234825 -3.57418319 -3.57419641 -3.57419642 -3.55112710 -3.54922372 -3.54922856 -3.549229
volume of cell :       43.07
CONTCAR：
3.710000000000000
     0.0000000000000000   -0.8434465875014347    0.0000000000000000
     0.7304461715115043    0.4217232937507173    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.3689584800243790
B：采用一次优化法：
与上面方法不同的是INCAR中取ISIF＝3，取定的实验缩放系数a＝3.21，计算的精度选得更高，且仅做一次优化。
结果为：
TOTEN  =        -3.549568 eV
volume of cell :       43.12
CONTCAR：
Mg-hex                                 
    3.210000000000000   
     0.0000000000000000   -0.9801232359458133    0.0000000000000000
     0.8488116211684474    0.4900616179729067    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.5669819270266370
   2
Direct
  0.6666666666666643  0.3333333333333357  0.7500000000000000
  0.3333333333333357  0.6666666666666643  0.2500000000000000

结论：
1、特殊位置的粒子随晶格形状的改变而变化，不论允许驰豫与否。
2、不同ISIF的取值对晶格优化影响不大，本例中缩放系数均在3.21~3.31之范围内
3、采用分步优化法得出的优化晶格比较可靠
回复此楼

» 收录本帖的淘帖专辑推荐

VASP 计算-vasp VASP相关主题 科研
vasp

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2010-06-07 17:14:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
大概清楚了；）
非常感谢；）


我只知道pwscf可以做物态方程拟合
但也仅限于集中特殊的晶包(比如立方，六方)
赞一下(1人) 回复此楼
6楼2010-06-09 09:28:15
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

xirainbow

木虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiaohunhun(金币+1):谢谢 2010-06-08 23:39:11
做的非常好；）
谢谢；）

有个小问题：）
    “II、取ISIF＝3”，和“B：采用一次优化法”
都是ISIF=3.这两个有什么区别呢？
还有，为什么有的volume of cell 是七十多，有的四十多呢？
谢谢；）
赞一下(2人) 回复此楼
2楼2010-06-08 14:47:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奥古拉斯

无虫 (小有名气)
★ ★
xiaohunhun(金币+2):谢谢 2010-06-08 23:39:26
引用回帖:
Originally posted by xirainbow at 2010-06-08 14:47:21:
做的非常好；）
谢谢；）

有个小问题：）
    “II、取ISIF＝3”，和“B：采用一次优化法”
都是ISIF=3.这两个有什么区别呢？
还有，为什么有的volume of cell 是七十多，有的四十多呢？
谢谢；）

谢谢鼓励

1. Q:  “II、取ISIF＝3”，和“B：采用一次优化法”
        都是ISIF=3.这两个有什么区别呢？
    A:“II、取ISIF＝3”，这样取定ISIF后，会对所有的缩放系数试探值下的结构全优化，也就是在每次优化过程中离子可以驰豫，晶格形状允许改变。这样就得到了一系列的V～E数对，然后进行状态方程的拟合，从而得到最稳定的结构情况；而“B：采用一次优化法”，仅仅是取定了实验上得到的缩放系数，在这个基础上进行结构的全优化，而且仅作这一次运算，所得到的稳定结构情况在输出文件CONTCAR中。
2.Q:，为什么有的volume of cell 是七十多，有的四十多呢?
   A:看I和II这两种情况，I是ISIF＝2，II是ISIF＝3，然后你看再前面的ISIF参数说明，就会知道ISIF＝2在优化过程中要保持c/a，以及volume不变，而ISIF＝3则相反；同时我已经说明了我列出的volume of cell都是在缩放系数为3.71时的情况。据此可以这样看，ISIF＝2求出的结构是ISIF＝3的初始结构，经过可变晶格的优化计算，当然会得到与初始结构完全不同的结构。
呵呵，希望我解释清楚了
  
赞一下(2人) 回复此楼
3楼2010-06-08 17:16:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiaohunhun(金币+2):谢谢 2010-06-08 23:39:36
非常非常感谢您给出如此详细的说明；）

但是对于第一个小问题，我依然有些疑问：）
既然ISIF都是=3，那么vasp将对所有的原子位置和晶包参数进行优化。
那么A和B中的ISIF=3的时，在相同的输入结构的前提下，他们的计算过程和结果有差别么？

另：2.91              3.01               3.11              3.21              3.31              3.41              3.51              3.61 3.71
以上的这几个参数，是您人工指定的，还是vasp在优化时，自己产生的中间结构？
如果是vasp自己产生的，您如何保证A中的ISIF=不同数值的时候，都取这几个参数？
再：将V～E进行状态方程的拟合。vasp有提供这个功能么？还请指教；）

[ Last edited by xirainbow on 2010-6-8 at 20:36 ]
赞一下(2人) 回复此楼
4楼2010-06-08 20:27:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭