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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【原创】晶格优化中ISIF取值的探究（六方MgO为例）
51/1返回列表
查看: 3306  |  回复: 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的回帖

ufo991

铁虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
楼主为何比较缩放系数为3.71时候的各参量呢，为么不比较能量最低时候的呢？而一次性优化ISIF=3用的是实验值，这样有可比性吗？结论是咋得出来的呢，没看明白。
赞一下(1人) 回复此楼
9楼2011-07-19 09:56:04
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工学调剂 +6 用户573181 2026-03-20 10/500 2026-03-23 11:04 by 用户573181
[考研] 298求调剂 +8 上岸6666@ 2026-03-20 8/400 2026-03-23 11:02 by laoshidan
[考研] 0703 调剂 +3 我可以上岸的对� 2026-03-16 6/300 2026-03-23 08:52 by 醉在风里
[考研] 石河子大学（211、双一流）硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 310求调剂 +4 baibai1314 2026-03-16 4/200 2026-03-22 20:19 by edmund7
[考研] 生物学071000 329分求调剂 +5 我爱生物生物爱� 2026-03-17 5/250 2026-03-22 16:42 by tcx007
[考研] 生物学调剂 +5 Surekei 2026-03-21 5/250 2026-03-22 14:39 by tcx007
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 材料与化工（0856）304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 材料工程专硕 348分求调剂 +3 冬辞. 2026-03-17 5/250 2026-03-21 18:47 by 学员8dgXkO
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 299求调剂 +4 某某某某位 2026-03-21 4/200 2026-03-21 16:30 by barlinike
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
信息提示
关闭
请填处理意见
关闭