24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1356  |  回复: 2
【悬赏金币】回答本帖问题,作者大力2010将赠送您 10 个金币

大力2010

金虫 (著名写手)

[求助] 请教高手,微扰密度泛函方法(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.
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

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.

这不是写的很清楚吗?你计算的force constant文件要用超胞的POSCAR
2楼2024-02-11 20:55:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大力2010

金虫 (著名写手)

引用回帖:
2楼: Originally posted by dxcharlary at 2024-02-11 20:55:39
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.
...

我是参照找个教程 ( 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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 大力2010 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +5 否极泰来2026 2026-04-15 6/300 2026-04-16 09:43 by 学员JpLReM
[考研] 急需调剂 +5 绝不放弃22 2026-04-15 5/250 2026-04-15 23:11 by Tide man
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 一志愿沪9,生物学326求调剂 +7 刘墨墨 2026-04-15 7/350 2026-04-15 18:51 by 浮云166
[考研] 材料299专硕求调剂 +17 +21 2026-04-09 17/850 2026-04-15 17:20 by zs92450
[考研] 367求调剂 +11 hffQAQ 2026-04-09 11/550 2026-04-14 17:48 by lhj2009
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[考研] 0854调剂 +10 长弓傲 2026-04-11 11/550 2026-04-13 10:38 by wp06
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[教师之家] 山东双非院校考核超级无底线,领导幸灾乐祸,教师遭殃恐 +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 调剂 +6 青灯不负 2026-04-09 6/300 2026-04-11 20:35 by dongdian1
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 机械专硕270求调剂,接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 0854调剂 +8 950824he@ 2026-04-09 8/400 2026-04-11 10:11 by zhq0425
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 本科西工大 0856 324求调剂 +10 wysyjs25 2026-04-09 11/550 2026-04-10 08:37 by 5268321
信息提示
请填处理意见