24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 10920  |  回复: 37
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

conperint

金虫 (正式写手)

[求助] 用VASP-DFPT+phonopy计算石墨烯声子谱求助

本人用VASP-DFPT+phonopy计算声子谱,下图是我计算得到的石墨烯的声子谱
用VASP-DFPT+phonopy计算石墨烯声子谱求助

下图是一篇PRB文章里给出的石墨烯的声子谱(PHYSICAL REVIEW B 77, 125401 2008)
用VASP-DFPT+phonopy计算石墨烯声子谱求助-1

两个图形状差不多,也都有一点虚频(很小)。但是仔细观察可以发现,PRB文章算出的声子谱在Gamma点的频率为1586cm-1,而我算出的值为1300cm-1,其他点频率值也都不一样。我不清楚是什么原因,还请高手指导一下。(我测试了一下phonopy给的NaCl的例子,没有问题,会不会是计算二维材料声子谱跟块体材料声子谱的计算有不一样的地方呢?)

我在这里把我算声子谱的过程贴出来给大家分享一下:
(1)优化石墨烯晶格,我用的GGA-PBE,能量判据为10-4eV,力收敛到0.01eV/A。优化得到石墨烯晶格常数为2.4673A。这一步应该没有问题。

(2)按照phonopy给的提示,扩胞,我扩了4*4的胞,a,b方向长度都接近10A,总共32个C原子。利用该超胞计算力常数。K点设定6*6*1。计算力常数的INCAR文件贴出来如下:
SYSTEM = graphene
PREC = Accurate
ISTART = 0
ICHARG = 2
LREAL = .FALSE.
LWAVE = .FALSE.
LCHARG = .FALSE.

ENCUT = 550
EDIFF = 1.0E-8
GGA = PE
NELM = 60
NELMIN = 4

NSW = 1
IBRION = 8
ISIF = 2
POTIM = 0.2
ISMEAR = 0 !Gaussian
SIGMA = 0.02
IALGO = 38
ADDGRID = .TRUE.

(3)计算完力常数以后,根据phonopy里给出的提示,运行命令:phonopy --fc vasprun.xml,得到文件FORCE_CONSTANTS。

(4)运行命令:phonopy --dim="4 4 1" -c POSCAR-unitcell band.conf 得到band.yaml文件,然后采用phonopy提供的scirpts程序bandplot得到数据,导入origin画声子谱

band.conf文件设置如下:
ATOM_NAME = C
DIM = 4 4 1
PRIMITIVE_AXIS = 0.25 0.0 0.0 0.0 0.25 0.0 0.0 0.0 1.0
BAND = 0.0 0.0 0.0 0.08333333 0.08333333 0.0 0.125 0.0 0.0 0.0 0.0 0.0
BAND_POINTS = 200
FORCE_CONSTANTS = READ
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档 第一性原理和电化学 phonopy 数学资料专题集

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

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

bjwang

铁杆木虫 (正式写手)

我尝试计算石墨烯纳米带的声子谱,4x1x1扩胞,各位帮忙看看我的band.conf文件哪里有问题,画出图和文献不对(计算方法用的和2D 石墨烯的类似,石墨烯的和文献一样,我怀疑是我的band.conf文件有问题),谢谢!
ATOM_NAME = H C
DIM =4 1 1
PRIMITIVE_AXIS = 0.25 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0
BAND = 0.0 0.0 0.0 0.125 0.0 0.0
BAND_POINTS = 200
FORCE_CONSTANTS = READ
It is a shame for the soul to be first to give way in this life,when thy body does not give way.
26楼2016-04-07 10:49:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 conperint 的主题更新
信息提示
请填处理意见