| 查看: 10967 | 回复: 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 | 数学资料专题集 |
» 猜你喜欢
博士申请都是内定的吗？
已经有6人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有12人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
关于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
|
本帖内容被屏蔽 |
conperint3楼2015-02-01 11:27:03
2楼2015-01-31 21:33:12
4楼2015-02-01 14:44:15
5楼2015-02-01 14:44:54