应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 用VASP-DFPT+phonopy计算石墨烯声子谱求助
381/4返回列表1234下一页
查看: 11086  |  回复: 37
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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 数学资料专题集

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-01-31 21:26:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

chuanghua304

禁虫 (职业作家)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
conperint(ljw4010代发): 金币+40, 楼主奖赏 2015-02-04 11:20:04
ljw4010: 应助指数+1 2015-02-05 16:04:53
本帖内容被屏蔽

» 本帖已获得的红花（最新10朵）

conperint
3楼2015-02-01 11:27:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

conperint

金虫 (正式写手)
引用回帖:
3楼: Originally posted by chuanghua304 at 2015-02-01 11:27:03
首先看文献用的是LDA还是PBE？SIGMA = 0.02太小了，用0.1就行了。还有就是看人家都到底是扩的多大的包。

谢谢回复。

（1）文献里没说扩的多大胞，这个也就没法比较了。我又重新扩了6*6*1的胞，算出的声子谱只有3支声学波，光学波完全没有。

（2）文献里泛函用的是LDA，我用的PBE。我用LDA重新算了一下声子谱。画出来的Gamma点LO和TO的频率反而变得更小了，只有1000cm-1多一点，与文献的数据1586cm-1差的更多了，很奇怪。。。不知道是什么原因？

LDA泛函画出的声子谱
用VASP-DFPT+phonopy计算石墨烯声子谱求助-2
赞一下(1人) 回复此楼
4楼2015-02-01 14:44:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuanghua304

禁虫 (职业作家)
本帖内容被屏蔽
6楼2015-02-01 15:01:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
conperint(ljw4010代发): 金币+10, 铁柱奖赏 2015-02-04 11:20:40
ljw4010: 抱歉打字快了，上为帖主奖赏 2015-02-04 11:21:40
ljw4010: 应助指数+1 2015-02-05 16:05:04
石墨烯为半金属，用MP型展宽力常数算得更准确。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

conperint
9楼2015-02-01 18:00:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

conperint

金虫 (正式写手)
我刚才又仔细看了一下，只有LO和TO这两支光学波在Gamma点附近不一样，其他地方两幅图基本完全一样。
赞一下 回复此楼
2楼2015-01-31 21:33:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

conperint

金虫 (正式写手)
送红花一朵
引用回帖:
3楼: Originally posted by chuanghua304 at 2015-02-01 11:27:03
首先看文献用的是LDA还是PBE？SIGMA = 0.02太小了，用0.1就行了。还有就是看人家都到底是扩的多大的包。

你这不是应助贴，我没法给你评分
赞一下 回复此楼
5楼2015-02-01 14:44:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

conperint

金虫 (正式写手)
引用回帖:
6楼: Originally posted by chuanghua304 at 2015-02-01 15:01:07
没事，LDA和PBE计算的频率有些差别很正常。另外你可以用IBRION=6算声子，其他和IBRION=8是一样的。

好的，我试试IBRION=6,看看。另外，我刚刚试着讲SIGMA改为0.1，其他的不便，算出来的频率为1520cm-1，跟文献值接近很多了！看来这个参数影响还是蛮大的
赞一下 回复此楼
7楼2015-02-01 16:45:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

conperint

金虫 (正式写手)
引用回帖:
6楼: Originally posted by chuanghua304 at 2015-02-01 15:01:07
没事，LDA和PBE计算的频率有些差别很正常。另外你可以用IBRION=6算声子，其他和IBRION=8是一样的。

我还想请教一下，phonopy里建议IBRION=8，如果IBRION=6，算出力常数以后，仍然使用phonopy照同样的步骤去画声子谱吗？
赞一下 回复此楼
8楼2015-02-01 16:49:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

conperint

金虫 (正式写手)
引用回帖:
9楼: Originally posted by liqizuiyang at 2015-02-01 18:00:26
石墨烯为半金属，用MP型展宽力常数算得更准确。

好的，谢谢，我试试看
赞一下 回复此楼
10楼2015-02-01 18:43:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 conperint 的主题更新
381/4返回列表1234下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭