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

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

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 的主题更新

返回列表

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

[考研] 266分，求材料相关专业调剂 +6	哇呼哼呼哼 2026-03-30	8/400	2026-03-30 23:48 by 哇呼哼呼哼
[考研] 085701环境工程求调剂 +11	多久上课 2026-03-27	12/600	2026-03-30 21:21 by 研究僧导导
[考研] 11408总分309，一志愿东南大学求调剂，不挑专业 +5	天赋带到THU 2026-03-29	6/300	2026-03-30 20:49 by dick_runner
[考研] 抱歉 +3	田洪有 2026-03-30	3/150	2026-03-30 19:11 by 迷糊CCPs
[考研] 334分一志愿武理材料求调剂 +16	李李不服输 2026-03-26	16/800	2026-03-30 17:47 by wangjy2002
[考研] 求调剂 +7	青春裁为三截 2026-03-29	7/350	2026-03-30 13:14 by laoshidan
[考研] 311求调剂 +10	lin0039 2026-03-26	10/500	2026-03-30 10:26 by herarysara
[考研] 一志愿南航 335分 \| 0856 \| GPA 4.07 \| 有科研经历 +8	cccchenso 2026-03-29	8/400	2026-03-29 23:53 by 我是小康
[考研] 275求调剂 +15	Micky11223 2026-03-25	20/1000	2026-03-29 20:44 by 唐沐儿
[考研] 329求调剂 +10	钮恩雪 2026-03-25	10/500	2026-03-29 13:32 by peike
[考研] 316求调剂 +7	江辞666 2026-03-26	7/350	2026-03-28 21:28 by sanrepian
[考研] 311（085601）求调剂 +4	liziyeyeye 2026-03-28	4/200	2026-03-28 18:50 by 535743368
[考研] 299求调剂 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐儿
[考研] 求调剂推荐材料 304 +15	荷包蛋hyj 2026-03-26	15/750	2026-03-28 04:13 by fmesaito
[考研] 求调剂 +4	零八# 2026-03-27	4/200	2026-03-27 18:07 by yu221
[考研] 复试调剂，一志愿南农083200食品科学与工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麦当当
[考研] 化学调剂一志愿上海交通大学336分-本科上海211 +4	小鱼爱有机 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 打过很多竞赛，085406控制工程300分，求调剂 +3	askeladz 2026-03-26	3/150	2026-03-26 09:08 by 给你你注意休息
[考研] 一志愿吉林大学材料与化工303分求调剂 +4	为学666 2026-03-24	4/200	2026-03-25 11:27 by BruceLiu320
[考研] 277分求调剂，跨调材料 +3	考研调剂lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky

24小时热门版块排行榜

forestwolf9291

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

» 猜你喜欢

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

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

buct2010

dxf_dxf_dxf

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

查看全部散金贴