24小时热门版块排行榜    

查看: 4295  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注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的回帖
查看全部 3 个回答

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 321求调剂 +4 Ymlll 2026-03-24 4/200 2026-03-24 14:44 by sprinining
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 【双一流院校新能源、环境材料,材料加工与模拟招收大量调剂】 +4 Higraduate 2026-03-22 7/350 2026-03-24 11:23 by 种大树
[考研] 一志愿哈工大,085400,320,求调剂 +3 gdlf9999 2026-03-24 3/150 2026-03-24 10:08 by 搏击518
[考研] 284求调剂 +10 Zhao anqi 2026-03-22 10/500 2026-03-24 00:08 by Equinoxhua
[考研] 327求调剂 +5 prayer13 2026-03-23 5/250 2026-03-23 22:11 by 星空星月
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 一志愿武理材料工程348求调剂 +6  ̄^ ̄゜汗 2026-03-19 9/450 2026-03-23 19:53 by pswait
[考研] 招08考数学 +6 laoshidan 2026-03-20 14/700 2026-03-23 14:37 by 15614359529
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 一志愿东华大学化学070300,求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
信息提示
请填处理意见