| 查看: 10914 | 回复: 37 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
用VASP-DFPT+phonopy计算石墨烯声子谱求助
|
|||
本人用VASP-DFPT+phonopy计算声子谱，下图是我计算得到的石墨烯的声子谱 下图是一篇PRB文章里给出的石墨烯的声子谱（PHYSICAL REVIEW B 77, 125401 2008）  两个图形状差不多，也都有一点虚频（很小）。但是仔细观察可以发现，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 |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
孩子确诊有中度注意力缺陷
已经有6人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
送红花一朵
 |
11楼2015-02-01 18:44:08
2楼2015-01-31 21:33:12
chuanghua304
禁虫 (职业作家)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
conperint(ljw4010代发): 金币+40, 楼主奖赏 2015-02-04 11:20:04
ljw4010: 应助指数+1 2015-02-05 16:04:53
conperint(ljw4010代发): 金币+40, 楼主奖赏 2015-02-04 11:20:04
ljw4010: 应助指数+1 2015-02-05 16:04:53
|
本帖内容被屏蔽 |
» 本帖已获得的红花（最新10朵）
3楼2015-02-01 11:27:03
4楼2015-02-01 14:44:15