24小时热门版块排行榜

返回列表

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

lzlgcdc

金虫 (正式写手)

应助: 2 (幼儿园)
金币: 677.6
散金: 17
红花: 2
帖子: 824
在线: 318.8小时
虫号: 1940713
注册: 2012-08-15

[求助] EV曲线计算问题已有2人参与

各位大侠求教一个问题在做EV曲线时 INCAR是以什么标准设置的参数？
我现在对一个结构在ISIF=3下优化过但是想算一下体积模量所以需要重算这时候怎么简便点呢？能不能就是做EV曲线的时候只改变晶格参数不改变原子位置静态计算得到的能量画出EV曲线？

感激啊

回复此楼

1楼2014-03-09 19:58:49

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

jpchou

专家顾问 (著名写手)

专家经验: +20
1ST强帖: 5
应助: 231 (大学生)
金币: 2342
散金: 4367
红花: 103
帖子: 1061
在线: 757.3小时
虫号: 585915
注册: 2008-08-08
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
lzlgcdc: 金币+5, ★★★★★最佳答案 2014-03-10 13:38:25
fzx2008: 金币+5, 谢谢指导 2014-03-11 09:09:52

手冊要多讀，把它當工具書?碛�
http://cms.mpi.univie.ac.at/vasp/guide/node112.html

ISIF    是否計算是否計算原子是否允許 cell 形狀是否允許 cell體積是否允許
   受力 (force) ?埩�(tensor) 移動 (relax) 改變 (shape) 改變 (volume)
------------------------------------------------------------------------------------------------------------------------------------------
0      yes    no yes      no      no
1      yes      trace only    yes      no      no
2      yes      yes      yes      no      no
3      yes      yes      yes      yes      yes
4      yes      yes      yes      yes      no
5      yes      yes      no      yes      no
6      yes      yes      no      yes      yes
7      yes      yes      no      no      yes
------------------------------------------------------------------------------------------------------------------------------------------
?倪@表可以看到
前兩項只影響 OUTCAR 內是否寫出 force and stress tensor
與計算本質無關 (且第一項全都是 YES, 根本沒得選)，後三項才是重點
常用的 ISIF = 3 所有欄位全都是 YES，重點是 cell 的形狀體積都是允許變化
另一個常用的是 ISIF=2，它是 cell 的形狀體積都是 "不" 允許變化
研究目的不同，選取不同 ISIF

回到你的問題

在做EV曲线时 INCAR是以什么标准设置的参数？
> 跟你做結構優化(ISIF=3)時的INCAR機乎一模一樣
唯一差別就是將 ISIF 改為 2

我现在对一个结构在ISIF=3下  优化过但是想算一下体积模量所以需要重算  这时候怎么简便点呢？
> "体积模量"...... 是什麼? bulk modules?
如果你指的是 bulk modules~ 那就是要重算, 就我所知應該沒有簡便的方法
做法就是你下面所說的內容

能不能就是做EV曲线的时候只改变晶格参数不改变原子位置  静态计算得到的能量  画出EV曲线？
> 這不是"能不能"的問題，EV曲線就是這樣做

EV 曲線就是 E(energy 能量)對V(volume 體積) 的作圖
通常 x 軸為體積，y軸為能量
體積主要是由什麼決定? --> lattice constant
所以 EV 曲線作法就是
  改變晶格常數 --> 固定原子位置 --> 不同晶格常數可以得到不同能量 --> 作出 EV 曲線

關於2樓 jghe 大俠所說的
"ISIF=3 是全优化。做EV曲线，要用ISIF=4 固定体积优化，然后拟合E-V。可以从你ISIF=3优化好的结构开始左右各取5个点，优化，算能量。"
小弟我�?在是忍不住想講些什麼，希望您能見�?....

1.  ISIF=4 與 ISIF=2 的差別是差在 4 是允許晶胞形狀變化，而 2 是不允許
   在高溫或高壓下晶格的相是會改變的 (phase transition)
   著名的就是 BCC-HCP or BCC-FCC phase transition
   若今天研究的主題是某個元素在高壓下 BCC 與 FCC 兩種結構穩定度比較
   此時對此元素在 BCC 與 FCC 下的 EV 曲線 (假設此元素高壓下 BCC 不穩定)
   由於是"高壓"，若是允許晶胞形狀改變 (ISIF=4)
   因為此元素高壓下 BCC 不穩定，則 BCC 作出?淼� EV 曲線，其每一個點都不會是 BCC structure (BCC shape)
   就會得到不正確的結果
2.  補充第一點
   EV 曲線是要  "在某個瞬間的體積" 與 "此瞬間體積之能量" 作圖
   也就是說 EV 曲線下每一個 "點" 的計算都只會有一個 ionic loop
   (non-SCF calculation，或叫 "單點能計算" single point calculation)
   VASP 是不允許做任何 optimization
   若是 ISIF=4，則VASP計算時就會對 structure shape 做 optimization
   這是不對的
3.  銜接第二點
   後面您所提到的 "可以从你ISIF=3优化好的结构开始左右各取5个点，优化，算能量"
   整句話是正確的，除了 "优化" 之外~
   如同我第二點所說，結構的 "优化" 在 EV 曲線計算下是不允許的