24小时热门版块排行榜    

查看: 4008  |  回复: 14

lzlgcdc

金虫 (正式写手)

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

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

感激啊
回复此楼

» 收录本帖的淘帖专辑推荐

soft study

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
lzlgcdc: 金币+5, ★★★★★最佳答案 2014-03-10 13:38:25
fzx2008: 金币+5, 谢谢指导 2014-03-11 09:09:52
???????x?????????????????
http://cms.mpi.univie.ac.at/vasp/guide/node112.html

ISIF           ??????                ??????                ?????????S        cell ?Π???????S                cell?w?e??????S
           ???? (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
?c????|?o?P (?????????? YES, ?????]???x)?????????????c
????? ISIF = 3 ???Й?λ????? YES?????c?? cell ???Π??w?e???????S???
???????????? ISIF=2?????? cell ???Π??w?e???? "??" ???S???
?о?????????x???? ISIF

????????}

????EV?????   INCAR???????????????????  
> ???????Y??????(ISIF=3)?r??INCAR?C???????
   Ψ???e????? ISIF ??? 2

??????????????ISIF=3??  ?????    ?????????????????    ???????????  ?????????????? ??
> "??????"...... ????N? bulk modules?
   ?????????? bulk modules~ ??????????, ???????????]?к???????  
   ?????????????????f?????

??????????EV???????????侧????? ????????λ??  ???????????????  ????EV?????  
> ?@????"?????"????}??EV?????????@????


EV ???????? E(energy ????)??V(volume ?w?e) ?????D
??? x ?S???w?e??y?S??????
?w?e?????????N?Q??? --> lattice constant
???? EV ????????????
  ???????? --> ??????λ?? --> ????????????????????? --> ???? EV ????



?P?2?? jghe ??b???f??
"ISIF=3 ??????????EV????????ISIF=4 ?????????????????E-V?????????ISIF=3????????????????5?????????????????"
С????????????????vЩ??N????????????....

1.  ISIF=4 ?c ISIF=2 ???e????? 4 ?????S?????Π???????? 2 ??????S
     ??????????????????????? (phase transition)
     ????????? BCC-HCP or BCC-FCC phase transition
     ???????о??????}?????????????? BCC ?c FCC ??N?Y??????????^
     ??r????????? BCC ?c FCC ??? EV ???? (???O????????? BCC ??????)
     ?????"???"?????????S?????Π??? (ISIF=4)
     ??????????? BCC ?????????t BCC ???????? EV ?????????????c???????? BCC structure (BCC shape)
     ???????????_??Y??
2.  ?a?????c
     EV ???????  "???????g???w?e" ?c "????g?w?e?????" ???D
     ??????f EV ?????????? "?c" ????????????? ionic loop
     (non-SCF calculation????? "???c?????" single point calculation)
     VASP ??????S???κ? optimization
     ???? ISIF=4???tVASP???r????? structure shape ?? optimization
     ?@???????
3.  ?????c
     ???????????? "???????ISIF=3????????????????5???????????????"
     ??????????_??????? "???" ???~
     ????????c???f???Y???? "???" ?? EV ???????????????S??
5楼2014-03-10 04:41:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

★ ★ ★
fzx2008: 金币+3, 谢谢指导 2014-03-11 09:25:56
???????:
6?: Originally posted by jghe at 2014-03-10 10:53:09
???????????????
????????????????????????????????????????
VASP???????????????????????????????????????EV?????????????????????棨???????????????????? ...

"VASP???????????????????????????????????????EV?????????????????????棨?????????????????????????????????ISIF=2??ISIF=4 ?????????????????????????б??... ISIF=2???к??????????????????????????????????"
> ???v????????_?????@????????N "????" ?? ISIF=4 ?????
   ???ε?ξ??]??N??????~????????f?? ISIF=2 ?? 4 ??????
   ???}??????~?????????????????б....?@?????Q???????????y
   ?? ISIF=2 ??_????к???????? (???????????h?x?????B?r)
   ?????@?????? EV curve ???????
   ?F???@???}???????? EV curve ???
   ??????????l??б?Y???? EV curve, ????????21?c?????g?????B?????????
   ???????侧???????^?L???^???????????????(????????g???B)
   ??r???????? ISIF=4?????????c??????????????? ISIF=2 ?????
   ?Π???????????????????
   ???????c??Y??????????g?????B??Y??????????Y????
   ?@??? EV curve .... ?????????



    BTW......?o???......... ????N?
12楼2014-03-10 16:32:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

jghe

金虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
lzlgcdc: 金币+5, ★★★★★最佳答案 2014-03-10 02:07:56
ISIF=3 是全优化。做EV曲线,要用ISIF=4 固定体积优化,然后拟合E-V。可以从你ISIF=3优化好的结构开始左右各取5个点,优化,算能量。
2楼2014-03-10 00:08:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzlgcdc

金虫 (正式写手)

引用回帖:
2楼: Originally posted by jghe at 2014-03-10 00:08:29
ISIF=3 是全优化。做EV曲线,要用ISIF=4 固定体积优化,然后拟合E-V。可以从你ISIF=3优化好的结构开始左右各取5个点,优化,算能量。

大侠  从你ISIF=3优化好的结构开始左右各取5个点 这句不太理解
能再解释下吗   再优化的话NSW设置多少合适呢?
3楼2014-03-10 02:08:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzlgcdc

金虫 (正式写手)

引用回帖:
2楼: Originally posted by jghe at 2014-03-10 00:08:29
ISIF=3 是全优化。做EV曲线,要用ISIF=4 固定体积优化,然后拟合E-V。可以从你ISIF=3优化好的结构开始左右各取5个点,优化,算能量。

大侠  我是想着 ISIF=3优化后得到了较为准确的结构   然后做EV曲线时候  晶格缩放后  继续优化  大概跑完也是在ISIF=3的那种结构
大侠意思这样是不对的   是不是说改变晶格参数后  再优化  材料的结构都变化了
4楼2014-03-10 02:49:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)

★ ★ ★ ★ ★
fzx2008: 金币+5, 谢谢指导 2014-03-11 09:23:27
???????:
5?: Originally posted by jpchou at 2014-03-10 04:41:41
???????x?????????????????
http://cms.mpi.univie.ac.at/vasp/guide/node112.html

ISIF           ??????                ??????                ?????????S        cell ?Π???????S                cell?w?e??????S
           ???? (force)        ????(tensor)        ??? ( ...

???????????????
????????????????????????????????????????
VASP???????????????????????????????????????EV?????????????????????棨?????????????????????????????????ISIF=2??ISIF=4 ?????????????????????????б??... ISIF=2???к??????????????????????????????????

?????????????????????????????????????????????籬????????????????????????????ISIF=4?????????????????????????? ????ν??????????????ISIF=2?????????a??b??c ???????????????????????????????????
6楼2014-03-10 10:53:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tianyali5665

新虫 (初入文坛)

★ ★ ★
fzx2008: 金币+3, 谢谢回帖! 2014-03-11 09:24:00
???????:
5?: Originally posted by jpchou at 2014-03-10 04:41:41
???????x?????????????????
http://cms.mpi.univie.ac.at/vasp/guide/node112.html

ISIF           ??????                ??????                ?????????S        cell ?Π???????S                cell?w?e??????S
           ???? (force)        ????(tensor)        ??? ( ...

???5?????3.  ?????c
     ???????????? "???????ISIF=3????????????????5???????????????"
      ??????????_??????? "???" ???~
      ????????c???f???Y???? "???" ?? EV ???????????????S???
?????????????????????????????????????????????????????е????????
b) ?????????????????????????????????ISIF ??????????ISIF ????????
2??4 ??5??
ISIF ?????????????????????????????????????????ISIF ??????????С???
???????????????????????????c/a ??b/a ?????
c) ???????????CONTCAR ??????POSCAR????????ξ???????
d) ??a ????10 ???????????????????????????????????????Volume-Etot ?????
7楼2014-03-10 11:29:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tianyali5665

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by jghe at 2014-03-10 10:53:09
没关系,交流而已!
我完全不同意你的观点,也可能我们所讨论的问题不一样。
VASP在优化的时候会考虑对称性的,如果有两个相,就要做两个EV曲线,在高压下就会观察到曲线交叉(相变)。你所说的大概是比较简单的 ...

我不明白的一点是,在侯老师的教程中,“以六角结构的Mg为例计算E-V关系,在结构优化的时候用了ISIF=5,即不改变元胞的体积,不进行原子位置的弛豫,但改变元胞的形状”,
疑问1
在计算E-V数据的时候,晶格常数要在实验晶格常数附近改变,这样晶格常数改变了,元胞的体积不就改变了吗?
疑问2
为什么要改变元胞的形状呢?
8楼2014-03-10 11:37:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzlgcdc

金虫 (正式写手)

在EV计算中   到底应该优化后静态计算E呢,还是直接改变晶格参数后 直接静态计算E呢?

我因为是做过一次优化了  得到较为准确的POSCAR    我想着  改变晶格参数   材料的结构也不应该变吧……所以就用改变晶格参数后直接静态计算了E       不知道这种想法是不是错误的     请楼上的各位大侠讲解一下  
感激不尽
9楼2014-03-10 13:35:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzlgcdc

金虫 (正式写手)

???????:
5?: Originally posted by jpchou at 2014-03-10 04:41:41
???????x?????????????????
http://cms.mpi.univie.ac.at/vasp/guide/node112.html

ISIF           ??????                ??????                ?????????S        cell ?Π???????S                cell?w?e??????S
           ???? (force)        ????(tensor)        ??? ( ...

????  лл?????????????????    ??????????????????  
???????е?????  ????
??EV??????   ?????????????????E???????????侧??????? ?????????E???

???????????????????  ??????????POSCAR    ??????  ??侧?????   ?????????????ɡ???????????侧??????????????????E       ????????????????????     ???????λ???????????  
?м?????
10楼2014-03-10 13:39:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lzlgcdc 的主题更新
信息提示
请填处理意见