| 查看: 12266 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理相关文档 | 第一性原理和电化学 | phonopy | 数学资料专题集 |
» 猜你喜欢
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有7人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于VASP+phonopy计算声子过程的讨论
已经有6人回复
- vasp计算声子谱已得FORCE_SETS 计算声子谱的时候出现以下错误提示 请求各位帮忙
已经有4人回复
- vasp压力下计算band、dos、声子谱的问题
已经有19人回复
- vasp5.2可以用linear response theory方法计算声子谱，声子态密度吗？？
已经有15人回复
- Vasp+phonopy 声子谱出现离奇虚频！
已经有25人回复
- 关于声子谱计算热膨胀系数的学位论文
已经有3人回复
- VASP+PHON算声子谱
已经有9人回复
- 声子谱计算耗内存 想把两台自己连起来一起算可以吗？
已经有8人回复
- 关于vasp超胞中LREAL如何设置
已经有7人回复
- 怎么用VASP计算声子谱
已经有4人回复
- QE-PWscf计算声子谱流程。
已经有10人回复
- 用VASP计算声子谱时出现错误
已经有3人回复
- 求助pwscf高手，改变k点之后声子谱出现很大虚频
已经有10人回复
- VASP+phonopy声子谱的问题
已经有18人回复
- 单层和多层石墨电子谱和声子谱研究
已经有66人回复
- 采用vasp ，phononpy或frophon计算声子谱的具体方法是什么？
已经有5人回复
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朵）
conperint3楼2015-02-01 11:27:03
4楼2015-02-01 14:44:15
chuanghua304
禁虫 (职业作家)
|
本帖内容被屏蔽 |
6楼2015-02-01 15:01:07
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2015-02-01 18:00:26
2楼2015-01-31 21:33:12
5楼2015-02-01 14:44:54
7楼2015-02-01 16:45:54
8楼2015-02-01 16:49:24
10楼2015-02-01 18:43:54
