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

24小时热门版块排行榜    

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

sungjen

银虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你这个只是在探究对晶格常数的影响，还有就是你一次性优化的晶格常数怎么只有3.21  和3.71差别这么大
赞一下 回复此楼
希望有一天成为大神
10楼2013-12-04 11:23:17
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085404 293求调剂 +9 勇远库爱314 2026-04-08 9/450 2026-04-12 02:24 by 秋豆菜芽
[考研] 280求调剂 +13 wzzz王 2026-04-09 13/650 2026-04-12 00:31 by 勇攀高峰0126
[考研] 085600材料与化工329分求调剂 +16 叶zilin 2026-04-10 16/800 2026-04-11 11:04 by may_新宇
[考研] 300分求调剂 （085501机械专硕，本科扬大） +8 xu@841019 2026-04-11 8/400 2026-04-11 10:46 by qingpingzhu
[考研] 还有化工二轮调剂的学校吗 5+14 化工人999 2026-04-09 48/2400 2026-04-11 10:27 by 89436494
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 298求调剂 +9 钉叮咚冬瓜 2026-04-07 11/550 2026-04-11 09:35 by zhq0425
[考研] 22408 366分，本科211，一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 中科院总分315求调剂 +8 lallalh 2026-04-09 8/400 2026-04-10 19:30 by dick_runner
[考研] 273求调剂 +51 麦小叮当 2026-04-06 58/2900 2026-04-10 15:54 by jiajinhpu
[考研] 0702物理学学硕299求调剂 +6 祁柒连 2026-04-06 6/300 2026-04-10 11:10 by Roomoo
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
[考研] 考研调剂-材料类-284 +28 想换手机不想解� 2026-04-08 28/1400 2026-04-09 20:08 by 倒数321?
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 349学科化学045106求调剂，化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 材料工程322 +18 哈哈哈吼吼吼哈 2026-04-07 19/950 2026-04-09 10:44 by cymywx
[考研] 334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 09:28 by wdyheheeh
[考研] 338求调剂 +5 小猪红色 678 2026-04-06 6/300 2026-04-07 21:18 by 乔哒哒哒
[考研] 专硕085403，291分，有两篇专利，一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
信息提示
关闭
请填处理意见
关闭