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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (6762)
>考研 (1585)
>导师招生 (1517)
>文献求助 (543)
>虫友互识 (354)
>考博 (232)
>休闲灌水 (214)
>基金申请 (208)
>硕博家园 (178)
>招聘信息布告栏 (138)
>博后之家 (125)
>论文投稿 (62)
>绿色求助(高悬赏) (51)
>找工作 (49)
>SciFinder/Reaxys (33)
>教师之家 (27)

返回列表

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 ]

回复此楼

» 猜你喜欢

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

一些凝聚态专业的知识和关于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的回帖

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的回帖

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的回帖

相关版块跳转我要订阅楼主 forestwolf9291 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 341求调剂(一志愿湖南大学070300) +4	番茄头--- 2026-03-22	5/250	2026-03-23 21:59 by fuyu_
[考研] 269求调剂 +4	我想读研11 2026-03-23	4/200	2026-03-23 21:25 by pswait
[考研] 一志愿重庆大学085700资源与环境，总分308求调剂 +7	墨墨漠 2026-03-23	8/400	2026-03-23 20:36 by Creta
[考研] 材料292调剂 +3	橘颂思美人 2026-03-23	3/150	2026-03-23 19:36 by pswait
[考研] 336化工调剂 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 一志愿东华大学化学070300，求调剂 +7	2117205181 2026-03-21	8/400	2026-03-22 22:55 by chixmc
[考研] 310求调剂 +4	baibai1314 2026-03-16	4/200	2026-03-22 20:19 by edmund7
[考研] 一志愿华中农业071010，总分320求调剂 +5	困困困困坤坤 2026-03-20	6/300	2026-03-22 17:41 by hxsm
[考研] 324求调剂 +6	lucky呀呀呀鸭 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 319求调剂 +4	小力气珂珂 2026-03-20	4/200	2026-03-22 15:53 by ColorlessPI
[考研] 354求调剂 +7	Tyoumou 2026-03-18	10/500	2026-03-22 11:11 by 人来盛
[考研] 初试 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 0805 316求调剂 +3	大雪深藏 2026-03-18	3/150	2026-03-21 18:55 by 学员8dgXkO
[考研] 材料工程专硕 348分求调剂 +3	冬辞. 2026-03-17	5/250	2026-03-21 18:47 by 学员8dgXkO
[考研] 311求调剂 +3	勇敢的小吴 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 一志愿华南师大 070300（化学）304分求调剂 +3	0703武芊慧雪304 2026-03-18	3/150	2026-03-21 00:48 by JourneyLucky
[考研] 22408 344分求调剂一志愿华电计算机技术 +4	solanXXX 2026-03-20	4/200	2026-03-20 23:49 by alg094825
[考研] 0703化学调剂 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 收复试调剂生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考博] 26博士申请 +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 轻松不少随

24小时热门版块排行榜

forestwolf9291

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

» 猜你喜欢

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

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

buct2010

dxf_dxf_dxf

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

查看全部散金贴