| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1876 | 回復(fù): 4 | ||
澤潤(rùn)東方金蟲(chóng) (正式寫(xiě)手)
|
[求助]
關(guān)于縮放系數(shù)的優(yōu)化 已有1人參與
|
|
關(guān)于VASP計(jì)算我些問(wèn)題想請(qǐng)教各位。 結(jié)構(gòu)弛豫:空間群是Pnma,三個(gè)晶格常數(shù)互不相等。 實(shí)驗(yàn)方面的結(jié)構(gòu)數(shù)據(jù)如下: **************************** a = 0.3846 nm, b = 1.1268 nm, c = 1.1321 nm, α = 90 °, β = 90 °, γ = 90 ° 。 **************************** 以實(shí)驗(yàn)結(jié)果為基礎(chǔ)開(kāi)始優(yōu)化,采用GGA-PBE泛函。步驟如下: 1.從0.98-1.8之間選取了二十多個(gè)數(shù)值作為縮放系數(shù)進(jìn)行了一系列的結(jié)構(gòu)弛豫。 2.選擇步驟1中能量最低的結(jié)果,作為初始結(jié)構(gòu),首先固定a方向晶格常數(shù),保持b/a和c/a不變(ISIF=4),進(jìn)行晶格常數(shù)優(yōu)化,但是保持a,b,c之間的夾角不變;其次再調(diào)整a方向的常數(shù),重復(fù)上面的步驟。反復(fù)進(jìn)行了十多個(gè)測(cè)試。 3.最后選擇步驟2中得到的結(jié)構(gòu)作為最穩(wěn)定的結(jié)構(gòu)。如下: **************************** direct lattice vectors 3.915449903 0.000000000 0.000000000 0.000000000 11.390399806 0.000000000 0.000000000 0.000000000 11.492100194 **************************** 4.在步驟3的基礎(chǔ)上重新做了結(jié)構(gòu)弛豫,測(cè)試了K點(diǎn)、Ecut。最后確定下來(lái)結(jié)構(gòu)和Ecut。 之后的靜待自洽、能帶以及態(tài)密度的計(jì)算都是步驟的結(jié)構(gòu)為基礎(chǔ)進(jìn)行的。 目前我比較了一下我的結(jié)構(gòu)跟之前發(fā)表的理論值與實(shí)驗(yàn)值。誤差稍微有些大: present EV-GGA(已發(fā)表的理論值) Exp1 Exp2 Exp3 a(angstr) 3.87 3.84 3.85 3.83 3.84 b(angstr) 11.21 11.22 11.27 11.23 11.22 c(angstr) 11.96 11.30 11.32 11.29 11.31 Vol(angstr3) 519.8 - 490.62 - - 我的問(wèn)題如下: a. 在此基礎(chǔ)上所做的能帶和態(tài)密度可信度高么 ? b. 如果不可信,應(yīng)如何從新開(kāi)始結(jié)構(gòu)參數(shù)的優(yōu)化 ? c. 出現(xiàn)計(jì)算誤差的原因有哪些 ? |

用戶注銷(xiāo) (職業(yè)作家)
金蟲(chóng) (正式寫(xiě)手)
|
最近按照侯老師的建議做了一下。取在實(shí)驗(yàn)值的左右取10來(lái)個(gè)體積值,進(jìn)行ISIF=4的計(jì)算,得到體積---總能的10來(lái)對(duì)數(shù)據(jù),然后狀態(tài)方程擬合,得到一個(gè)能量最低所對(duì)應(yīng) 的體積值(即平衡態(tài)下的體積值),然后以此體積值作為基礎(chǔ),計(jì)算出一個(gè)初始的a, b, c,再對(duì)此體積下結(jié)構(gòu)進(jìn)行ISIF=4計(jì)算。 但是計(jì)算出來(lái)的結(jié)果與實(shí)驗(yàn)值相差還是很大。但是跟2013年JPCC上面發(fā)表的一篇文章的結(jié)果比較接近。我想這多半不是我的方法的問(wèn)題。 |

用戶注銷(xiāo) (職業(yè)作家)
至尊木蟲(chóng) (職業(yè)作家)

| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求材料調(diào)劑 +8 | 隔壁陳先生 2026-03-12 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料工程專(zhuān)碩330分求調(diào)劑 +5 | 小材化本科 2026-03-18 | 5/250 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 材料與化工專(zhuān)碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 一志愿,福州大學(xué)材料專(zhuān)碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[論文投稿]
投稿問(wèn)題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|