24小时热门版块排行榜    

查看: 4291  |  回复: 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的回帖

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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 306求0703调剂一志愿华中师范 +7 纸鱼ly 2026-03-21 8/400 2026-03-23 23:31 by chixmc
[考研] 一志愿河北工业大学0817化工278分求调剂 +7 jhybd 2026-03-23 11/550 2026-03-23 23:29 by 呆呆师姐
[考研] 085600材料与化工调剂 +7 A-哆啦Z梦 2026-03-23 12/600 2026-03-23 23:16 by 星空星月
[考研] 材料专业求调剂 +11 hanamiko 2026-03-18 11/550 2026-03-23 23:12 by peike
[考研] 361求调剂 +3 Glack 2026-03-22 3/150 2026-03-23 22:03 by fuyu_
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 311求调剂 +6 冬十三 2026-03-18 6/300 2026-03-22 20:18 by edmund7
[考研] 材料与化工085600,总分304,本科有两篇sci参与,求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 一志愿华中农业071010,总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 寻找调剂 +4 倔强芒? 2026-03-21 4/200 2026-03-22 16:14 by 木托莫露露
[考研] 286分人工智能专业请求调剂愿意跨考! +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 295求调剂 +4 一志愿京区211 2026-03-18 6/300 2026-03-20 23:41 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
信息提示
请填处理意见