| 3 | 1/1 | 返回列表 |
| 查看: 1293 | 回復(fù): 2 | |||
| 【懸賞金幣】回答本帖問題,作者大力2010將贈(zèng)送您 10 個(gè)金幣 | |||
[求助]
請(qǐng)教高手,微擾密度泛函方法(DFPT)計(jì)算聲子譜,phonopy得不到band.yaml文件,咋辦 已有1人參與
|
|||
|
VASP結(jié)合Phonopy計(jì)算聲子譜,微擾密度泛函方法(DFPT),計(jì)算結(jié)束后運(yùn)行phonopy --fc vasprun.xml得到FORCE_CONSTANTS文件,準(zhǔn)備band.conf文件, 運(yùn)行phonopy --dim 2 2 2 -c POSCAR-unitcell band.conf得不到band.yaml文件。出現(xiàn)提示: "band.conf" was read as phonopy configuration file. Crystal structure was read from "POSCAR-unitcell". Unit of length: angstrom Band structure mode Settings: Force constants symmetrization: on Supercell: [2 2 2] Spacegroup: Pbcm (57) Use -v option to watch primitive cell, unit cell, and supercell structures. Force constants are read from "FORCE_CONSTANTS". Number of atoms in supercell (224) is not consistent with the matrix shape of force constants (28, 28) read from FORCE_CONSTANTS. Please carefully check DIM, FORCE_CONSTANTS, and POSCAR-unitcell. |
專家顧問 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +83 |
|
我是參照找個(gè)教程 ( force constant 文件是在計(jì)算結(jié)束后運(yùn)行phonopy --fc vasprun.xml得到的 ): 微擾密度泛函方法(DFPT) 將SPOSCAR復(fù)制為POSCAR,INCAR設(shè)置如下 Global Parameters ISTART = 0 ISPIN = 1 ICHARG = 2 LREAL = .FALSE. ENCUT = 500 PREC = A LWAVE = .FALSE. LCHARG = .FALSE. ADDGRID= .TRUE. Electronic Relaxation ISMEAR = 0 SIGMA = 0.05 NELM = 90 NELMIN = 6 EDIFF = 1E-08 GGA = PE ALGO = N Ionic Relaxation NSW = 1 IBRION = 8 EDIFFG = -0.01 KPOINTS為 K-Spacing Value to Generate K-Mesh: 0.040 0 Gamma 4 4 4 0.0 0.0 0.0 計(jì)算結(jié)束后運(yùn)行phonopy --fc vasprun.xml得到FORCE_CONSTANTS文件。 準(zhǔn)備band.conf文件,如下 ATOM_NAME = Si DIM = 2 2 2 BAND = 0.0 0.0 0.0 0.5 0.5 0.0, 0.5 0.5 1.0 0.0 0.0 0.0, 0.0 0.0 0.0 0.5 0.5 0.5 BAND_POINTS = 51 BAND_LABELS = $\Gamma$ X, X $\Gamma$, $\Gamma$ L BAND_CONNECTION = .TRUE. FORCE_CONSTANTS = READ FC_SYMMETRY = .TRUE. 運(yùn)行phonopy --dim 2 2 2 -c POSCAR-unitcell band.conf得到band.yaml文件。 運(yùn)行phonopy-bandplot --gnuplot > band.dat得到可用Origin畫圖的band.dat文件。 注意:如果想要直接出圖,可以運(yùn)行phonopy -c POSCAR-unitcell band.conf -p -s,之后會(huì)得到band.pdf文件等,里邊繪制好了聲子譜。 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281電子信息求調(diào)劑 +4 | jhtfeybgj 2026-03-02 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
[考研] 一志愿985材料與化工 326分求調(diào)劑 +3 | Hz795795 2026-03-04 | 3/150 |
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
[考研] 307求調(diào)劑 +4 | xuyuan111 2026-03-04 | 5/250 |
|
|
[考研] 085600 材料與化工 298 +14 | 小西笑嘻嘻 2026-03-03 | 14/700 |
|
|
[考研] 085601 材料工程 320 +6 | 和樂瑤 2026-03-03 | 6/300 |
|
|
[考博] 26申博 +5 | north, 2026-02-28 | 5/250 |
|
|
[論文投稿]
EST拒稿重投
5+3
|
15102603076 2026-03-02 | 3/150 |
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5 | Qwertyuop 2026-03-01 | 5/250 |
|
|
[考研] 26考研報(bào)考西工大材料308分求調(diào)劑 +4 | weizhong123 2026-03-01 | 5/250 |
|
|
[考研] 材料類考研調(diào)劑 +6 | gemmgemm 2026-03-01 | 7/350 |
|
|
[考研] 298求調(diào)劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
[考研] 材料284求調(diào)劑,一志愿鄭州大學(xué)英一數(shù)二專碩 +15 | 想上岸的土撥鼠 2026-02-28 | 15/750 |
|
|
[考研] 295求調(diào)劑 +8 | 19171856320 2026-02-28 | 8/400 |
|
|
[考研] 化工專碩342,一志愿大連理工大學(xué),求調(diào)劑 +6 | kyf化工 2026-02-28 | 7/350 |
|
|
[考研] 274求調(diào)劑 +3 | cgyzqwn 2026-03-01 | 7/350 |
|
|
[考研] 275求調(diào)劑 +3 | L-xin? 2026-03-01 | 6/300 |
|
|
[基金申請(qǐng)] 成果系統(tǒng)訪問量大,請(qǐng)一小時(shí)后再嘗試。---NSFC啥時(shí)候好哦,已經(jīng)兩天這樣了 +4 | NSFC2026我來了 2026-02-28 | 4/200 |
|