版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

forestwolf9291

木虫 (著名写手)

应助: 7 (幼儿园)
金币: 2702.3
帖子: 1048
在线: 289.7小时
虫号: 695237

[交流] 关于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 ]

回复此楼

» 猜你喜欢

有没有人能给点建议已经有5人回复
假如你的研究生提出不合理要求已经有12人回复
实验室接单子已经有7人回复
全日制（定向）博士已经有5人回复
萌生出自己或许不适合搞科研的想法，现在跑or等等看？已经有4人回复
Materials Today Chemistry审稿周期已经有4人回复
参与限项已经有3人回复
对氯苯硼酸纯化已经有3人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复

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

一些凝聚态专业的知识和关于vasp的知识已经有83人回复
关于vasp收敛和优化已经有24人回复
关于VASP一些琐事的求助（K点，吸附优化等）已经有12人回复
求助关于vasp磁矩设置的说明书的详细理解已经有12人回复
关于vasp计算自旋轨道耦合的问题已经有4人回复
求助成功：VASP进行体积弛豫时用ISIF＝7，如何看弛豫后的结果已经有7人回复
新手求助~关于VASP位置优化的问题已经有6人回复
【求助】vasp中COOP如何理解及输出文件如何处理？已经有12人回复
【求助】vasp优化时选择了isif=3 怎么限制晶格的一个维度呢？已经有16人回复
【求助】关于vasp进行分子动力学模拟的求助已经有12人回复
【求助】简单问题----关于vasp晶格优化已经有7人回复
【求助】关于VASP中Kpoint设置问题已经有6人回复
【求助】关于vasp态密度计算已经有10人回复
【求助】vasp中关于电场单位的问题已经有12人回复
【原创】晶格优化中ISIF取值的探究（六方MgO为例）已经有9人回复
【讨论】用vasp计算过弹性常数的请进！！更新中！！已经有52人回复

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

查看全部散金贴

1楼 2012-04-28 11:03:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dxf_dxf_dxf

银虫 (小有名气)

1ST强帖: 4
应助: 15 (小学生)
金币: 269.1
帖子: 104
在线: 29.1小时
虫号: 1449737

★ ★ ★ ★ ★ ★
小木虫: 金币+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-----------------------

引用回帖:

---------------------ISIF=1-----------------------

引用回帖:

---------------------ISIF=2-----------------------

引用回帖:

---------------------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-----------------------

引用回帖:

---------------------ISIF=5-----------------------

引用回帖:

---------------------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

金虫 (正式写手)

1ST强帖: 10
应助: 27 (小学生)
金币: 2328.6
帖子: 351
在线: 231.2小时
虫号: 1129027

★ ★ ★ ★ ★
小木虫: 金币+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的回帖

查看全部 3 个回答

24小时热门版块排行榜

forestwolf9291

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

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到: 查看全部散金贴

dxf_dxf_dxf

buct2010

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

查看全部散金贴