 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1296 | 回復(fù): 2 | |||
| 【懸賞金幣】回答本帖問題，作者大力2010將贈送您 10 個(gè)金幣 | |||
[求助]
請教高手，微擾密度泛函方法(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. |
回復(fù)此樓» 猜你喜歡
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有19人回復(fù)
-
085600材料與化工，一志愿廣州985，求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料學(xué)學(xué)碩308分/本科東北大學(xué)/一志愿西北工業(yè)大學(xué)/
已經(jīng)有8人回復(fù)
-
334求調(diào)劑
已經(jīng)有4人回復(fù)
-
0856材料與化工求調(diào)劑！
已經(jīng)有8人回復(fù)
-
一志愿武理085601專碩347分求調(diào)劑
已經(jīng)有4人回復(fù)
-
招調(diào)劑
已經(jīng)有6人回復(fù)
-
材料與化工，291，求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑，學(xué)校研究所都可以，材料與化工267分
已經(jīng)有4人回復(fù)
-
332材料求調(diào)劑
已經(jīng)有4人回復(fù)
dxcharlary
專家顧問 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +83 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 251 (大學(xué)生)
- 金幣: 10414.6
- 散金: 155
- 紅花: 65
- 帖子: 3137
- 在線: 757.4小時(shí)
- 蟲號: 936718
- 注冊: 2010-01-04
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
- 管轄: 第一性原理
2樓2024-02-11 20:55:39
|
我是參照找個(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，之后會得到band.pdf文件等，里邊繪制好了聲子譜。 |
3樓2024-02-21 10:37:25
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 大力2010 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招調(diào)劑 +4 | 帆船哥 2026-03-04 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 沒上岸的看過來 +3 | tangxiaotian 2026-03-01 | 5/250 |
|
|
[考研] 304分材料專碩求調(diào)劑 +11 | qiuzhigril 2026-03-03 | 14/700 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-05 | 5/250 |
|
|
|
[考研] 一志愿985材料與化工 326分求調(diào)劑 +3 | Hz795795 2026-03-04 | 3/150 |
|
|
|
[考研] 320材料與化工，求調(diào)劑 +6 | 鶴遨予卿 2026-03-04 | 8/400 |
|
|
|
[考研] 307求調(diào)劑 +4 | xuyuan111 2026-03-04 | 5/250 |
|
|
|
[考研] 學(xué)碩材料275調(diào)劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|
|
[考研] 本科太原理工采礦工程，求調(diào)劑 +3 | onlx 2026-03-01 | 3/150 |
|
|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
 5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
|
[考研] 291求調(diào)劑 +4 | Afy123456 2026-03-03 | 7/350 |
|
|
|
[考研] 環(huán)境調(diào)劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
|
[考研] 324求調(diào)劑 +4 | wxz2 2026-03-03 | 5/250 |
|
|
|
[考研] 材料工程274求調(diào)劑 +5 | Lilithan 2026-03-01 | 5/250 |
|
|
|
[考博] 博士自薦 +4 | kkluvs 2026-02-28 | 5/250 |
|
|
|
[考研] 一志愿中石油（華東）本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-02 | 3/150 |
|
|
|
[考研] 一志愿華南理工大學(xué)材料與化工326分，求調(diào)劑 +3 | wujinrui1 2026-02-28 | 3/150 |
|
|
|
[基金申請] 此成果不能導(dǎo)入原因：元數(shù)據(jù)必填信息不完整，可 進(jìn)行補(bǔ)充。 +4 | Kittylucky 2026-03-02 | 5/250 |
|
|
|
[考博] 26申博 +4 | 想申博！ 2026-02-26 | 6/300 |
|