| 1 | 1/1 | 返回列表 |
| 查看: 711 | 回復(fù): 0 | ||
[求助]
計算體模量時,晶格常數(shù)該如何變化?
|
|
若是從模擬實(shí)驗(yàn)的角度出發(fā),體積的變化不應(yīng)該超過1%,也就是晶格常數(shù)的取值以 0.002 A 增加。 但是,按照VASP手冊中的說法,為了消除 Puley stress 的影響,體積的變化要大一些,有的教程推薦使用 0.1 A 的間隔。 在我的計算中,當(dāng)使用 0.002 A 時,同一 K-point 取值下的一系列計算結(jié)果可以很好的擬合出 Birch-Murnaghan 方程。但是計算出的體模量并不能隨著 K-point 取值的增加而穩(wěn)定/收斂。 如果使用 0.1 A 時,則甚至連同一 K-point 取值下的一系列計算都不能很好的擬合出狀態(tài)方程。 請問,按照大家的計算經(jīng)驗(yàn),應(yīng)該如何選擇?謝謝! P.S. 我剛開始使用 VASP,現(xiàn)在主要進(jìn)行一些小的計算上手。上面的結(jié)果來自于計算 NiFe 合金。 具體的輸入文件如下: #! /bin/bash for i in `seq -w 3 1 15` do mkdir $i cd $i for j in `seq -w 3.0 0.1 4.0` do mkdir $j cd $j cat > KPOINTS << EOF automatic k-mesh generation 0 G $i $i $i 0 0 0 EOF cat > POSCAR << EOF fcc Ni_3_Fe_1 $j 1 0 0 0 1 0 0 0 1 3 1 c 0 0 0 0 0.5 0.5 0.5 0 0.5 0.5 0.5 0 EOF cat > INCAR << EOF SYSTEM = Ni_3_Fe_1 ISTART = 0 job: 0-new 1-cont 2-samecut ICHARG = 2 charge: 1-file 2-atom 10-const ENCUT = 350 eV ISMEAR = 2 SIGMA = 0.1 LREAL = .FALSE. PREC = Normal NGX = 16 NGY = 16 NGZ = 16 ISPIN = 2 MAGMOM = 1 1 1 5 NPAR = 2 IBRION = 2 NSW = 300 EOF |
找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +4 | 申子申申 2026-03-19 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +5 | Ncdx123456 2026-03-19 | 5/250 |
|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
[考研] 281求調(diào)劑(0805) +11 | 煙汐憶海 2026-03-16 | 21/1050 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 287求調(diào)劑 +3 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|