| 查看: 2026 | 回复: 6 | ||
zll201101金虫 (正式写手)
|
[求助]
VASP+Phon计算声子谱，如何引入LO_TO劈裂？ 已有1人参与
|
|
初学Vasp，现在想计算LiH的声子谱。目前已经计算出了Born有效电荷，该如何利用Phon引入Gamma点的LO-TO劈裂？ 必需用phon+QE吗？还是直接在输入文件INPHON中加入某个指令？ 如果用现有结果，利用Ｐｈｏｎｏｐｙ计算，该如何操作？ |
回复此楼» 猜你喜欢
回收溶剂求助
已经有7人回复
-
职称评审没过，求安慰
已经有40人回复
-
硝基苯如何除去
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
求助文献
已经有3人回复
2楼2018-08-20 19:48:26
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
|
上一个回答的是的。这个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
zll201101
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4883.7
- 散金: 56
- 红花: 2
- 帖子: 784
- 在线: 120.5小时
- 虫号: 3371917
- 注册: 2014-08-18
- 专业: 凝聚态物性 II ：电子结构
送红花一朵 |
谢谢！我看到vasp手册上有四种计算BORN EFFECTIVE CHARGES 的方法，对于同一个例子，不同方法计算出的结果有差异，该怎样选择？ 还有，引入LO-TO劈裂时，只用phon+vasp能算吗？ 发自小木虫Android客户端 |
4楼2018-08-21 21:11:27
zll201101
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4883.7
- 散金: 56
- 红花: 2
- 帖子: 784
- 在线: 120.5小时
- 虫号: 3371917
- 注册: 2014-08-18
- 专业: 凝聚态物性 II ：电子结构
送红花一朵 |
谢谢！我看到vasp手册上有四种计算BORN EFFECTIVE CHARGES 的方法，对于同一个例子，不同方法计算出的结果有差异，该怎样选择？ 还有，引入LO-TO劈裂时，只用phon+vasp能算吗？目前phonopy有问题，不能正常使用。 发自小木虫Android客户端 |
5楼2018-08-21 21:15:28
6楼2018-08-21 21:32:56
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zll201101: 金币+20, ★★★★★最佳答案 2018-08-24 08:08:23
zll201101: 金币+20, ★★★★★最佳答案 2018-08-24 08:08:23
|
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