 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3804 | 回復(fù): 12 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||
葵cc_木蟲 (小有名氣)
|
[求助]
表面體系的結(jié)構(gòu)優(yōu)化 已有3人參與
|
|||
|
各位前輩，本人剛接觸VASP，現(xiàn)在想用它來計算表面體系，ENCUT、K-mesh、SIGMA都優(yōu)化好了，接下來想進(jìn)行結(jié)構(gòu)優(yōu)化。 看了侯老師的教程，關(guān)于表面體系的結(jié)構(gòu)優(yōu)化這一塊，候老師寫到“在對表面體系的結(jié)構(gòu)進(jìn)行優(yōu)化時，主要是對原子的位置進(jìn)行優(yōu)化，而不再對原胞（Slab模型得到的）大小進(jìn)行優(yōu)化，一般采用的是Selective Dynamics（也就是有選擇性的位置弛豫）�！� 我計算的體系是2x2x3的超胞，對照侯老師的教程進(jìn)行了Selective Dynamics，計算結(jié)果出現(xiàn)ZBRENT: fatal error in bracketing please rerun with smaller EDIFF, or copy CONTCAR to POSCAR and continue，重新設(shè)置之后也沒有什么效果。 問題： （1）表面體系的結(jié)構(gòu)優(yōu)化具體是怎么做的？有經(jīng)驗的前輩能分享一下你們的例子嗎？包括具體的步驟和設(shè)置，不勝感激~~~~ （2）不需要優(yōu)化晶格常數(shù)，而需要優(yōu)化原子位置又是為什么？ （3）像以上的Selective Dynamics有選擇性的位置弛豫，又是怎么判斷該體系的原子已經(jīng)弛豫到了最佳位置。 謝謝~~~ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 表面 | 光催化及vasp計算 |
» 猜你喜歡
283分材料與化工求調(diào)劑
已經(jīng)有12人回復(fù)
-
349求調(diào)劑
已經(jīng)有7人回復(fù)
-
296材料專碩求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿上海海洋大學(xué)083200食品學(xué)碩，求調(diào)劑，接受其他專業(yè)
已經(jīng)有7人回復(fù)
-
材料340分調(diào)劑
已經(jīng)有8人回復(fù)
-
316求調(diào)劑
已經(jīng)有13人回復(fù)
-
321求調(diào)劑
已經(jīng)有11人回復(fù)
-
315求調(diào)劑
已經(jīng)有6人回復(fù)
-
330分求調(diào)劑
已經(jīng)有13人回復(fù)
-
化工調(diào)劑303分，過四級
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 這樣的優(yōu)化問題有哪些方法可以做？
已經(jīng)有11人回復(fù)
- vasp可以做體系激發(fā)態(tài)的結(jié)構(gòu)優(yōu)化和計算嗎？
已經(jīng)有7人回復(fù)
- 表面吸附體系的結(jié)構(gòu)優(yōu)化（弛豫）無法收斂 求分析求指導(dǎo)！
已經(jīng)有5人回復(fù)
- Pwscf 固定原子（表面體系）算聲子譜
已經(jīng)有4人回復(fù)
- MS Discover 動力學(xué)求助
已經(jīng)有10人回復(fù)
- 如何用VASP計算體系的介電常數(shù)
已經(jīng)有24人回復(fù)
- 收斂測試和結(jié)構(gòu)優(yōu)化
已經(jīng)有10人回復(fù)
- vasp表面計算的結(jié)構(gòu)優(yōu)化問題
已經(jīng)有3人回復(fù)
- vasp表面計算的結(jié)構(gòu)優(yōu)化問題
已經(jīng)有7人回復(fù)
- 請教一下那款免費軟件的結(jié)構(gòu)優(yōu)化功能比較好用？
已經(jīng)有7人回復(fù)
- 計算Cl-在Fe表面吸附問題時體系加電荷方法
已經(jīng)有27人回復(fù)
- 兩個表面之間作用能或原子吸附能的計算疑惑
已經(jīng)有9人回復(fù)
- 【活動or資源】vasp晶體結(jié)構(gòu)優(yōu)化
已經(jīng)有142人回復(fù)
- 【求助】金紅石結(jié)構(gòu)TiO2（110）表面體系能帶結(jié)構(gòu)計算得到的結(jié)果只有一半
已經(jīng)有5人回復(fù)
- 【求助完畢】vasp-slab-邊界附近原子頂位吸附小分子的結(jié)構(gòu)優(yōu)化很難收斂嗎？
已經(jīng)有5人回復(fù)
- 【求助】結(jié)構(gòu)優(yōu)化時如何將原子位置固定在原胞內(nèi)
已經(jīng)有9人回復(fù)
- 【原創(chuàng)】關(guān)于結(jié)構(gòu)弛豫設(shè)置及收斂判據(jù)選擇的一點體會
已經(jīng)有140人回復(fù)
- 【求助】請教分子動力學(xué)模擬
已經(jīng)有8人回復(fù)
- 【求助】如何用VASP 計算體系的表面能
已經(jīng)有11人回復(fù)
9樓2014-04-25 10:15:04
2樓2014-04-24 17:04:03
mywai520
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 339 (大學(xué)生)
- 金幣: 10425.6
- 紅花: 45
- 帖子: 1664
- 在線: 296.9小時
- 蟲號: 265142
- 注冊: 2006-07-09
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2014-04-24 23:59:39
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗: +20 - 1ST強(qiáng)帖: 5
- 應(yīng)助: 231 (大學(xué)生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時
- 蟲號: 585915
- 注冊: 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
葵cc_: 金幣+5, ★有幫助 2014-04-25 10:32:06
fzx2008: 金幣+5, 謝謝指導(dǎo)！辛苦了 2014-04-28 19:22:39
葵cc_: 金幣+5, ★有幫助 2014-04-25 10:32:06
fzx2008: 金幣+5, 謝謝指導(dǎo)！辛苦了 2014-04-28 19:22:39
|
我想~應(yīng)該是對於 POSCAR 內(nèi)容無法與所學(xué)的固態(tài)物理無法對應(yīng) 才造成迷惑 POSCAR內(nèi)主要分成三個部份 第二行可以當(dāng)成是 lattice constant --> A 第三~五行是 lattice vector --> VA "Direct" 這個關(guān)鍵字以下所有內(nèi)容都是 atomic position --> p 無論分子(molecule, cluster) or 表面(surface) or 塊材(bulk)系統(tǒng) ~ 皆是由這三個參數(shù)所組成 所以每顆原子在 cartesian coordinate 座標(biāo)可簡單寫成 A*VA*p 以上，都與固態(tài)所學(xué)是相對應(yīng)的 你所謂的 "對原胞大小進(jìn)行優(yōu)化" 計算中主要控制參數(shù)就是INCAR內(nèi)中的 ISIF=3 而 "在對表面體系的結(jié)構(gòu)進(jìn)行優(yōu)化時，主要是對原子的位置進(jìn)行優(yōu)化" 對應(yīng)的就是計算中 INCAR 內(nèi)的參數(shù) ISIF=2 請仔細(xì)讀 ISIF 參數(shù)說明 http://cms.mpi.univie.ac.at/vasp/guide/node112.html 簡單說 ISIF 就是控制計算過程中 VA & p 是否允許變化(或是說，是否固定) 所謂的結(jié)構(gòu)優(yōu)化就是優(yōu)化這兩個參數(shù) ISIF=3 是 VA & p 都允許優(yōu)化 而 ISIF=2 是 VA 固定， p 允許優(yōu)化 塊材系統(tǒng)使用 ISIF=3 我想你應(yīng)該是可以理解 但~ 可不可以用 ISIF=2 呢? 答案是可以的，這個我也猜想你應(yīng)該是知道的，所以就不詳述 表面系統(tǒng) (由原子與真空組成，真空層大多數(shù)會使用約 10 A) 一般~ 學(xué)習(xí)過程都會教你對於真空系統(tǒng)，就是用 ISIF=2 也就是只有 p 允許優(yōu)化 但~ 可不可以跟塊材一樣，使用 ISIF=3 去優(yōu)化呢? 計算上~當(dāng)然是可以的 如果你計算資源夠多，時間上也很空閒的話，你可以試試看用 ISIF=3 去優(yōu)化真空層約 10A 的表面結(jié)構(gòu) 你就會發(fā)現(xiàn) 你的真空層會一直被吃掉 原本 10A 的真空，一直算一直算，可能會變成 7~8 A 你再把 7~8 A 的 CONTCAR 拿?砝^續(xù)算下去~ VASP 還是會乖乖的一直算一直算~然後可能真空變成 5 A 若是在努力的算下去，最後表面系統(tǒng)可能就會變回成塊材 理由很簡單 這是因為 VASP 計算中，結(jié)構(gòu)是週期性的呈現(xiàn) 建構(gòu) 10 A 真空的表面系統(tǒng) 對表面系統(tǒng)而言 其頭上與腳下距離 10A 的地方都有一個長的跟它一樣的表面存在著 週期表內(nèi)所有元素都是傾向變成塊材 (塊材一定比表面穩(wěn)定) 所以這樣的計算最後一定回到塊材態(tài) 這就是為什麼 "在對表面體系的結(jié)構(gòu)進(jìn)行優(yōu)化時，主要是對原子的位置進(jìn)行優(yōu)化" 其中一個原因 |
4樓2014-04-25 02:26:28
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 330分求調(diào)劑 +13 | qzenlc 2026-03-29 | 13/650 |
|
|---|---|---|---|---|
|
|
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +16 | 睿08 2026-03-27 | 17/850 |
|
|
[考研] 292分，材料與化工，申請調(diào)劑 +19 | 程晴之 2026-04-01 | 23/1150 |
|
|
|
[考研] 324求調(diào)劑 +5 | 想上學(xué)求調(diào) 2026-04-01 | 6/300 |
|
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +14 | c297914 2026-03-28 | 14/700 |
|
|
|
[考研] 279求調(diào)劑 +6 | 學(xué)而思兮知 2026-04-01 | 6/300 |
|
|
|
[考研] 材料調(diào)劑 +12 | 一樣YWY 2026-04-01 | 12/600 |
|
|
|
[考研] 一志愿346上海大學(xué)生物學(xué) +3 | 上海大學(xué)346調(diào)劑 2026-04-01 | 3/150 |
|
|
|
[考研] 一志愿北交大材料工程，總分358 +4 | cs0106 2026-04-01 | 4/200 |
|
|
|
[考研] 08工科275分求調(diào)劑 +13 | AaAa7420 2026-03-31 | 13/650 |
|
|
|
[考研] 11408 321分求調(diào)劑 +3 | huchun12138 2026-03-30 | 4/200 |
|
|
|
[考研] 建環(huán)，能源，土木老師路過看一看�。�！ +4 | 嘿嘿uu 2026-04-01 | 4/200 |
|
|
|
[考研] 070300求調(diào)劑306分 +5 | 26要上岸 2026-03-27 | 5/250 |
|
|
|
[考研] 一志愿武理材料工程302調(diào)劑環(huán)化或化工 +15 | Doleres 2026-03-31 | 16/800 |
|
|
|
[考研] 085601英二數(shù)二求調(diào)劑 總分325 +4 | 余航航 2026-03-31 | 4/200 |
|
|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +4 | Yam. 2026-03-30 | 5/250 |
|
|
|
[考研] 334求調(diào)劑 +7 | Trying] 2026-03-31 | 7/350 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-29 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +16 | kikiki7 2026-03-30 | 16/800 |
|
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|