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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于VASP中ISIF取值的理解
31/1返回列表
查看: 4418  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

forestwolf9291

木虫 (著名写手)

[交流] 关于VASP中ISIF取值的理解

结构优化是下面一切计算的关键。VASP中队结构优化都是通过ISIF来控制的。
ISIF控制主要通过(1) 离子迟豫；(2) 形状改变； (3)体积改变 来实现。
(1). 离子迟豫：就是离子坐标相对基矢做变化，即Direct坐标变化；
(2). 形状改变：就三个基矢的大小和相对夹角的变化；
(3).体积改变：是指三个基矢按相同比例增减；
以上理解没有问题吧？
下面以金刚石硅为例，原胞含有两个原子。
(1). 如果仅允许离子迟豫，则就是这两个原子间相对位置的变化
(2). 如果仅允许形状改变，那原本三个基矢大小相等和夹角相同的关系就有可能被
打破
(3). 仅允许体积改变就是三个基矢的大小同比例增减。
如果体系是面心立方等，则原胞就一个原子，离子迟豫就没意义了。
不知上述理解有多少不妥，望指点！

[ Last edited by forestwolf9291 on 2012-4-28 at 18:57 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-04-28 11:03:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)
★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
uuv2010: 金币+2, 专家考核, 3q 2012-04-28 18:59:36
forestwolf9291: 金币+2, 谢谢关注。我的帖子的1,2,3是问题编号，已改为（1）（2）（3）。不好意思误会了。 2012-04-28 18:59:45
不是这样的，ISIF=3晶胞体积改变，原子也会弛豫的，形状也会改变
赞一下 回复此楼
2楼2012-04-28 18:18:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxf_dxf_dxf

银虫 (小有名气)
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
forestwolf9291: 金币+5, 谢谢，我也在测试。但没你搞得详细。再次感谢！ 2012-05-03 19:54:55
测试：ISIF的设定值对vasp几何优化结果的影响
体系：fcc-Si
输入文件：
见[url=http://muchong.com/bbs/viewthread.php?tid=4458443&fpage=1]http://muchong.com/bbs/viewthread.php?tid=4458443&fpage=1[/rul]5楼
测试语句：
CODE:
grep "lattice vectors" -A2 OUTCAR;grep "POSITION" -A3 OUTCAR

测试用脚本:
CODE:
#! /bin/bash
rm WAVECAR
rm ISIF-Geo.txt
for (( i = 0; i < 8 ; i++ ));do
cat >INCAR < general:
System = diamond Si
ISMEAR = 0; SIGMA = 0.1;
ENMAX  =  240
IBRION=2; ISIF=$i ; NSW=15
EDIFF  = 0.1E-04
EDIFFG = -0.01
!
echo "ENCUT= $i" ;
mpirun -np 8 $HOME/vasp/bin/vasp52mpi.x
echo "---------------------ISIF=$i-----------------------">>ISIF-Geo.txt
grep "lattice vectors" -A2 OUTCAR >> ISIF-Geo.txt
grep "POSITION" -A3 OUTCAR >> ISIF-Geo.txt
done

结果
---------------------ISIF=0-----------------------
引用回帖:
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.81250      4.81250      4.81250         0.000000      0.000000      0.000000
      0.68750      0.68750      0.68750         0.000000      0.000000      0.000000

---------------------ISIF=1-----------------------
引用回帖:
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.81250      4.81250      4.81250         0.000000      0.000000      0.000000
      0.68750      0.68750      0.68750         0.000000      0.000000      0.000000

---------------------ISIF=2-----------------------
引用回帖:
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.81250      4.81250      4.81250         0.000000      0.000000      0.000000
      0.68750      0.68750      0.68750         0.000000      0.000000      0.000000

---------------------ISIF=3-----------------------
引用回帖:
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.721745648  2.721745648    -0.183705630  0.183705630  0.183705630
     2.721745648  0.000000000  2.721745648     0.183705630 -0.183705630  0.183705630
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.732455265  2.732455265    -0.182985612  0.182985612  0.182985612
     2.732455265  0.000000000  2.732455265     0.182985612 -0.182985612  0.182985612
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.81250      4.81250      4.81250         0.000000      0.000000      0.000000
      0.68750      0.68750      0.68750         0.000000      0.000000      0.000000
--
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.76305      4.76305      4.76305         0.000000      0.000000      0.000000
      0.68044      0.68044      0.68044         0.000000      0.000000      0.000000
--
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.78180      4.78180      4.78180         0.000000      0.000000      0.000000
      0.68311      0.68311      0.68311         0.000000      0.000000      0.000000

---------------------ISIF=4-----------------------
引用回帖:
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.81250      4.81250      4.81250         0.000000      0.000000      0.000000
      0.68750      0.68750      0.68750         0.000000      0.000000      0.000000

---------------------ISIF=5-----------------------
引用回帖:
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.81250      4.81250      4.81250         0.000000      0.000000      0.000000
      0.68750      0.68750      0.68750         0.000000      0.000000      0.000000

---------------------ISIF=6-----------------------
引用回帖:
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.721744920  2.721744920    -0.183705680  0.183705680  0.183705680
     2.721744920  0.000000000  2.721744920     0.183705680 -0.183705680  0.183705680
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.732455491  2.732455491    -0.182985597  0.182985597  0.182985597
     2.732455491  0.000000000  2.732455491     0.182985597 -0.182985597  0.182985597
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.81250      4.81250      4.81250         0.000000      0.000000      0.000000
      0.68750      0.68750      0.68750         0.000000      0.000000      0.000000
--
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.76305      4.76305      4.76305         0.000000      0.000000      0.000000
      0.68044      0.68044      0.68044         0.000000      0.000000      0.000000
--
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.78180      4.78180      4.78180         0.000000      0.000000      0.000000
      0.68311      0.68311      0.68311         0.000000      0.000000      0.000000

---------------------ISIF=7-----------------------
引用回帖:
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.750000000  2.750000000    -0.181818182  0.181818182  0.181818182
     2.750000000  0.000000000  2.750000000     0.181818182 -0.181818182  0.181818182
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.721771879  2.721771879    -0.183703860  0.183703860  0.183703860
     2.721771879  0.000000000  2.721771879     0.183703860 -0.183703860  0.183703860
--
      direct lattice vectors                 reciprocal lattice vectors
     0.000000000  2.732449693  2.732449693    -0.182985986  0.182985986  0.182985986
     2.732449693  0.000000000  2.732449693     0.182985986 -0.182985986  0.182985986
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.81250      4.81250      4.81250         0.000000      0.000000      0.000000
      0.68750      0.68750      0.68750         0.000000      0.000000      0.000000
--
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.76310      4.76310      4.76310         0.000000      0.000000      0.000000
      0.68044      0.68044      0.68044         0.000000      0.000000      0.000000
--
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      4.78179      4.78179      4.78179         0.000000      0.000000      0.000000
      0.68311      0.68311      0.68311         0.000000      0.000000      0.000000
赞一下 回复此楼
3楼2012-05-03 19:09:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 forestwolf9291 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 wsYGJnURYvm6 2026-08-31 3/150 2026-09-01 01:52 by o3EMxHUAQDtS
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 01:30 by o3EMxHUAQDtS
[考博] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 01:30 by o3EMxHUAQDtS
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了，不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:23 by o3EMxHUAQDtS
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
关闭
请填处理意见
关闭