24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1146  |  回复: 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个字符以上)
信息提示
请填处理意见