| 5 | 1/1 | 返回列表 |
| 查看: 1880 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
澤潤東方金蟲 (正式寫手)
|
[求助]
關(guān)于縮放系數(shù)的優(yōu)化 已有1人參與
|
|
|
關(guān)于VASP計算我些問題想請教各位。 結(jié)構(gòu)弛豫:空間群是Pnma,三個晶格常數(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ǔ)開始優(yōu)化,采用GGA-PBE泛函。步驟如下: 1.從0.98-1.8之間選取了二十多個數(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)行了十多個測試。 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)弛豫,測試了K點(diǎn)、Ecut。最后確定下來結(jié)構(gòu)和Ecut。 之后的靜待自洽、能帶以及態(tài)密度的計算都是步驟的結(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 - - 我的問題如下: a. 在此基礎(chǔ)上所做的能帶和態(tài)密度可信度高么 ? b. 如果不可信,應(yīng)如何從新開始結(jié)構(gòu)參數(shù)的優(yōu)化 ? c. 出現(xiàn)計算誤差的原因有哪些 ? |

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

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
[考研] 331求調(diào)劑(0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[基金申請] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[基金申請]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|