| 3 | 1/1 | 返回列表 |
| 查看: 1298 | 回復(fù): 2 | |||
| 【懸賞金幣】回答本帖問題,作者大力2010將贈送您 10 個金幣 | |||
[求助]
請教高手,微擾密度泛函方法(DFPT)計算聲子譜,phonopy得不到band.yaml文件,咋辦 已有1人參與
|
|||
|
VASP結(jié)合Phonopy計算聲子譜,微擾密度泛函方法(DFPT),計算結(jié)束后運行phonopy --fc vasprun.xml得到FORCE_CONSTANTS文件,準(zhǔn)備band.conf文件, 運行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)驗: +83 |
|
我是參照找個教程 ( force constant 文件是在計算結(jié)束后運行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 計算結(jié)束后運行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. 運行phonopy --dim 2 2 2 -c POSCAR-unitcell band.conf得到band.yaml文件。 運行phonopy-bandplot --gnuplot > band.dat得到可用Origin畫圖的band.dat文件。 注意:如果想要直接出圖,可以運行phonopy -c POSCAR-unitcell band.conf -p -s,之后會得到band.pdf文件等,里邊繪制好了聲子譜。 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] A區(qū)一本交叉課題組,低分調(diào)劑,招收機械電子信息通信等交叉方向 +19 | lisimayy 2026-03-04 | 27/1350 |
|
|---|---|---|---|---|
|
[考研] 中科大材料299求調(diào)劑 +4 | DAIjiayo 2026-03-05 | 6/300 |
|
|
[考研] 接收調(diào)劑 +19 | 津萌津萌 2026-03-02 | 29/1450 |
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|
|
[考研] 334求調(diào)劑 +3 | Trying] 2026-03-03 | 4/200 |
|
|
[考研] 301求調(diào)劑 +4 | 李LJR 2026-03-04 | 4/200 |
|
|
[考研] 293求調(diào)劑 +3 | 是樂渝哇 2026-03-04 | 3/150 |
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
[考博] 26申博-目前4篇SCI一作 +3 | chen_2024 2026-03-02 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | G212 2026-03-03 | 5/250 |
|
|
[考研] 環(huán)境調(diào)劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
[考研] 化工專碩調(diào)劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|
|
[考研] 一志愿華中科技大學(xué),化學(xué)專業(yè)344分,求調(diào)劑 +6 | 邢xing1 2026-03-02 | 6/300 |
|
|
[考研] 268求調(diào)劑 +10 | 簡單點0 2026-03-02 | 14/700 |
|
|
[考研]
|
好好好1233 2026-02-28 | 16/800 |
|
|
[考研] 求調(diào)劑 +7 | repeatt?t 2026-02-28 | 7/350 |
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
[考研] 265分求調(diào)劑不調(diào)專業(yè)和學(xué)校有行學(xué)上就 +6 | 禮堂丁真258 2026-02-28 | 9/450 |
|
|
[考研] 328求調(diào)劑 +3 | aaadim 2026-03-01 | 5/250 |
|
|
[考研] 311求調(diào)劑 +6 | 亭亭亭01 2026-03-01 | 6/300 |
|