| 24小時(shí)熱門版塊排行榜 |
| 1 | 1/1 | 返回列表 |
| 查看: 714 | 回復(fù): 0 | |||
[求助]
計(jì)算體模量時(shí),晶格常數(shù)該如何變化?
|
|
若是從模擬實(shí)驗(yàn)的角度出發(fā),體積的變化不應(yīng)該超過(guò)1%,也就是晶格常數(shù)的取值以 0.002 A 增加。 但是,按照VASP手冊(cè)中的說(shuō)法,為了消除 Puley stress 的影響,體積的變化要大一些,有的教程推薦使用 0.1 A 的間隔。 在我的計(jì)算中,當(dāng)使用 0.002 A 時(shí),同一 K-point 取值下的一系列計(jì)算結(jié)果可以很好的擬合出 Birch-Murnaghan 方程。但是計(jì)算出的體模量并不能隨著 K-point 取值的增加而穩(wěn)定/收斂。 如果使用 0.1 A 時(shí),則甚至連同一 K-point 取值下的一系列計(jì)算都不能很好的擬合出狀態(tài)方程。 請(qǐng)問(wèn),按照大家的計(jì)算經(jīng)驗(yàn),應(yīng)該如何選擇?謝謝! P.S. 我剛開(kāi)始使用 VASP,現(xiàn)在主要進(jìn)行一些小的計(jì)算上手。上面的結(jié)果來(lái)自于計(jì)算 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ā)表 | |
|---|---|---|---|---|
|
[考研]
|
小羊36 2026-03-30 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿西交288化工專碩求調(diào)劑 +4 | 好運(yùn)好運(yùn)接接 2026-03-24 | 4/200 |
|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +6 | 材料學(xué)257求調(diào)劑 2026-03-28 | 7/350 |
|
|
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +4 | 睿08 2026-03-27 | 4/200 |
|
|
[考研] 哈爾濱工業(yè)大學(xué)材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +9 | @taotao 2026-03-29 | 9/450 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
[考研] 367求調(diào)劑 +5 | 芋泥啵! 2026-03-28 | 5/250 |
|
|
[考研] 359求調(diào)劑 +5 | 王了個(gè)楠 2026-03-25 | 5/250 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 297求調(diào)劑 +17 | 田洪有 2026-03-26 | 18/900 |
|
|
[考研] 292求調(diào)劑 +13 | 是妍子也是研子 2026-03-30 | 13/650 |
|
|
[考研] 材料專碩 085600求調(diào)劑 +7 | BBQ233 2026-03-30 | 7/350 |
|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
[考研] 化學(xué)308分調(diào)劑 +7 | 你好明天你好 2026-03-30 | 8/400 |
|
|
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +9 | 天天好運(yùn)來(lái)上岸?/a> 2026-03-24 | 10/500 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|