| 查看: 3977 | 回复: 14 | ||||
[求助]
EV曲线计算问题已有2人参与
|
||||
|
各位大侠 求教一个问题 在做EV曲线时 INCAR是以什么标准设置的参数？ 我现在对一个结构在ISIF=3下 优化过 但是想算一下体积模量 所以需要重算 这时候怎么简便点呢 ？ 能不能就是做EV曲线的时候只改变晶格参数 不改变原子位置 静态计算得到的能量 画出EV曲线？ 感激啊 |
回复此楼» 收录本帖的淘帖专辑推荐
 soft study |
» 猜你喜欢
导师想让我从独立一作变成了共一第一
已经有8人回复
-
博士读完未来一定会好吗
已经有23人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
bulk modulus 在vasp下如何计算？
已经有4人回复- 关于matlab两个自变量、一个因变量曲线拟合问题
已经有3人回复
- 用castep计算单斜晶系弹性常数的问题
已经有7人回复
- 如何计算反应活化能
已经有7人回复
- 如何通过UV吸收曲线计算direct bandgap和indirect-bandgap
已经有8人回复
- vasp计算应力--应变曲线时应力值发生突变的问题
已经有8人回复
- 请教一个充放电曲线计算比电容问题
已经有9人回复
- 求助origin曲线拟合的一个问题
已经有6人回复
- 掺杂后直接计算能量，能说明什么问题？
已经有6人回复
- 使用peakfit软件 曲线拟合和分峰 的问题求助
已经有7人回复
- 请教各位计算牛人有关HOMO/LUMO能极差的问题，万急！
已经有29人回复
- 求问为什么使用DMOL3模块优化石墨烯结构的时候没有能量曲线？
已经有13人回复
- 循环伏安曲线基础问题--氧化峰，还原峰，氧化电位，还原电位及计算HOMO，LUMO的方法
已经有6人回复
- vasp计算中，如何选择正确的K点
已经有7人回复
- 【求助】如何计算应变对应的应力值
已经有5人回复
- 有关origin曲线拟合的问题
已经有7人回复
- 【求助】截断能优化该选择曲线的最低点还是继续优化以致稳定？
已经有24人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
- 【求助】放电曲线计算比能量和比功率的问题
已经有6人回复
- 【求助】计算出来的导带位置和价带位置和实验有差别吗？怎么修正？谢谢了！
已经有9人回复
- 【求助】如何解决自旋极化计算时的难以收敛问题？
已经有16人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
lzlgcdc: 金币+5, ★★★★★最佳答案 2014-03-10 13:38:25
fzx2008: 金币+5, 谢谢指导 2014-03-11 09:09:52
感谢参与，应助指数 +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 曲線計算下是不允許的 |
5楼2014-03-10 04:41:41
★ ★ ★
fzx2008: 金币+3, 谢谢指导 2014-03-11 09:25:56
fzx2008: 金币+3, 谢谢指导 2014-03-11 09:25:56
|
"VASP在优化的时候会考虑对称性的，如果有两个相，就要做两个EV曲线，在高压下就会观察到曲线交叉（相变）。你所说的大概是比较简单的体系，如立方，ISIF=2和ISIF=4 是一样的。对于四方，正交，单斜，... ISIF=2会有很大的晶格（各向异性）应力，因为你没有优化晶格。" > 你講得完全正確，但這也就是為什麼 "不能" 用 ISIF=4 的原因 簡單的單晶沒什麼好討論的~就如同你說的 ISIF=2 或 4 是一樣的 問題就是在~若是四方，正交，單斜....這�?對稱性相對低的系統 用 ISIF=2 的確是會有很大的晶格應力 (當晶格常數已經遠離最穩定態時) 可是這就是做 EV curve 所需要的 現在這問題不是在討論 EV curve 嬤? 假如我做一條單斜結構的 EV curve, 總共用了21點，中間是最穩態，能量最低 左右兩端其晶格常數已經過長或過短，因此能量會很高(相對於中間最穩態) 此時，若是用 ISIF=4，左右兩端點，其能量肯定比用 ISIF=2 ?淼牡� 形狀被改變，晶格角度就跟著改變 左右兩端點的結構已經跟中間穩定態的結構是不同一個結構了 這樣的 EV curve .... 可以接受嬤? BTW......靜水壓......... 是什麼? |
12楼2014-03-10 16:32:44
2楼2014-03-10 00:08:29
3楼2014-03-10 02:08:58
4楼2014-03-10 02:49:50
★ ★ ★ ★ ★
fzx2008: 金币+5, 谢谢指导 2014-03-11 09:23:27
fzx2008: 金币+5, 谢谢指导 2014-03-11 09:23:27
|
没关系，交流而已！ 我完全不同意你的观点，也可能我们所讨论的问题不一样。 VASP在优化的时候会考虑对称性的，如果有两个相，就要做两个EV曲线，在高压下就会观察到曲线交叉（相变）。你所说的大概是比较简单的体系，如立方，ISIF=2和ISIF=4 是一样的。对于四方，正交，单斜，... ISIF=2会有很大的晶格（各向异性）应力，因为你没有优化晶格。 我不明白你所说的瞬间体积，瞬间能量。你是说动态过程，如爆炸过程？对于我的理解，模拟静水压必须用ISIF=4，因为晶体受到的压力是各向同性的 （所谓静水压）。如你用ISIF=2，你用的不同a，b，c （但是相同的体积），得到的结果肯定是不同的。 |
6楼2014-03-10 10:53:09
★ ★ ★
fzx2008: 金币+3, 谢谢回帖！ 2014-03-11 09:24:00
fzx2008: 金币+3, 谢谢回帖！ 2014-03-11 09:24:00
|
针对5楼说的“3. 銜接第二點 後面您所提到的 "可以从你ISIF=3优化好的结构开始左右各取5个点，优化，算能量" 整句話是正確的，除了 "优化" 之外~ 如同我第二點所說，結構的 "优化" 在 EV 曲線計算下是不允許的” 根据侯老师的经验手册，是要做优化的，优化完后计算静态能量。以下是手册中的部分内容 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
8楼2014-03-10 11:37:06
9楼2014-03-10 13:35:48
10楼2014-03-10 13:39:19