24小时热门版块排行榜    

查看: 1622  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +16 羊腰板 2026-08-21 17/850 2026-08-23 23:02 by tianxiaochun
[基金申请] 让我中一个面上吧! +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 范进中举一文的中心思想 +3 炎黄贵胄 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见