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

登录注册

应《网络安全法》要求，自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 ]

回复此楼

» 猜你喜欢

求调剂一志愿武汉理工大学材料工程（085601）已经有7人回复
一志愿C9材料与化工专业总分300求调剂已经有4人回复
321求调剂已经有4人回复
307求调剂已经有4人回复
293求调剂已经有6人回复
085602 289分求调剂已经有3人回复
一志愿北京化工大学材料与化工 264分各科过A区国家线已经有3人回复
277分求调剂，跨调材料已经有3人回复
材料292调剂已经有6人回复
材料与化工考研调剂已经有6人回复

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

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

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

[考研] 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

24小时热门版块排行榜

forestwolf9291

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

» 猜你喜欢

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

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

dxf_dxf_dxf

buct2010

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

查看全部散金贴