| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1878 | 回復(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ì)算誤差的原因有哪些 ? |

金蟲(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)題。 |

至尊木蟲(chóng) (職業(yè)作家)

| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求材料調(diào)劑 +8 | 隔壁陳先生 2026-03-12 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +7 | Liwangman 2026-03-15 | 7/350 |
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +7 | rare12345 2026-03-18 | 7/350 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機(jī)器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
[考研] 考研調(diào)劑 +4 | 芬達(dá)46 2026-03-12 | 4/200 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|