| 查看: 1570 | 回复: 2 | |||
| 【悬赏金币】回答本帖问题，作者大力2010将赠送您 10 个金币 | |||
[求助]
请教高手，微扰密度泛函方法(DFPT)计算声子谱，phonopy得不到band.yaml文件，咋办 已有1人参与
|
|||
|
VASP结合Phonopy计算声子谱，微扰密度泛函方法(DFPT)，计算结束后运行phonopy --fc vasprun.xml得到FORCE_CONSTANTS文件，准备band.conf文件， 运行phonopy --dim 2 2 2 -c POSCAR-unitcell band.conf得不到band.yaml文件。出现提示： "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. |
回复此楼» 猜你喜欢
准备明年的基金了
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有15人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有13人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有19人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有19人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有15人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有19人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有24人回复
dxcharlary
专家顾问 (职业作家)
-
专家经验: +83 - 1ST强帖: 2
- 应助: 251 (大学生)
- 金币: 10505.9
- 散金: 155
- 红花: 65
- 帖子: 3141
- 在线: 758.7小时
- 虫号: 936718
- 注册: 2010-01-04
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
2楼2024-02-11 20:55:39
|
我是参照找个教程 （ force constant 文件是在计算结束后运行phonopy --fc vasprun.xml得到的 ）： 微扰密度泛函方法(DFPT) 将SPOSCAR复制为POSCAR，INCAR设置如下 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 计算结束后运行phonopy --fc vasprun.xml得到FORCE_CONSTANTS文件。 准备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楼2024-02-21 10:37:25
50