24小时热门版块排行榜    

查看: 2013  |  回复: 6

zll201101

金虫 (正式写手)

[求助] VASP+Phon计算声子谱,如何引入LO_TO劈裂?已有1人参与

初学Vasp,现在想计算LiH的声子谱。目前已经计算出了Born有效电荷,该如何利用Phon引入Gamma点的LO-TO劈裂?
必需用phon+QE吗?还是直接在输入文件INPHON中加入某个指令?
如果用现有结果,利用Phonopy计算,该如何操作?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (职业作家)

本帖仅楼主可见

» 本帖已获得的红花(最新10朵)

2楼2018-08-20 19:48:26
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

上一个回答的是的。这个phonopy手册上有详细的介绍。

https://atztogo.github.io/phonopy/vasp.html
见该网页上Non-analytical term correction (Optional)之后的部分。

由于phonopy版本的不同,可能该网页提到的命令比如phonopy-vasp-born会略有变化。

» 本帖已获得的红花(最新10朵)

3楼2018-08-21 13:29:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zll201101

金虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by zyl雨田 at 2018-08-20 19:48:26
计算的时候incar里面多一个参数,就是考虑波恩有效电荷的那个参数

谢谢!我看到vasp手册上有四种计算BORN EFFECTIVE CHARGES 的方法,对于同一个例子,不同方法计算出的结果有差异,该怎样选择?
还有,引入LO-TO劈裂时,只用phon+vasp能算吗?

发自小木虫Android客户端
4楼2018-08-21 21:11:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zll201101

金虫 (正式写手)

送红花一朵
引用回帖:
3楼: Originally posted by valenhou001 at 2018-08-21 13:29:44
上一个回答的是的。这个phonopy手册上有详细的介绍。

https://atztogo.github.io/phonopy/vasp.html
见该网页上Non-analytical term correction (Optional)之后的部分。

由于phonopy版本的不同,可能该网页提 ...

谢谢!我看到vasp手册上有四种计算BORN EFFECTIVE CHARGES 的方法,对于同一个例子,不同方法计算出的结果有差异,该怎样选择?
还有,引入LO-TO劈裂时,只用phon+vasp能算吗?目前phonopy有问题,不能正常使用。

发自小木虫Android客户端
5楼2018-08-21 21:15:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (职业作家)

本帖仅楼主可见
6楼2018-08-21 21:32:56
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zll201101: 金币+20, ★★★★★最佳答案 2018-08-24 08:08:23
引用回帖:
5楼: Originally posted by zll201101 at 2018-08-21 22:15:28
谢谢!我看到vasp手册上有四种计算BORN EFFECTIVE CHARGES 的方法,对于同一个例子,不同方法计算出的结果有差异,该怎样选择?
还有,引入LO-TO劈裂时,只用phon+vasp能算吗?目前phonopy有问题,不能正常使用。
...

phon和phonopy是两个不同的软件。所用的原理基本上都是一样的,都是通过small displacement的方法来做。前者是fotran代码的,目前好像还不直接支持加入LO-TO,需要自己手动的写代码加入。

你看下这篇文章上对各种计算声子谱程序的综述和比较:
https://www.nature.com/articles/npjcompumats20166

BORN EFFECTIVE CHARGES 的计算选用那种计算方法最简便的就可以了(比如直接在INCAR里加入LEPSILON=.TRUE.和IBRION = 8以及其他相关的参数)。各种方法算出来的结果会有细微的差别,比较正常。
7楼2018-08-23 22:07:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zll201101 的主题更新
信息提示
请填处理意见