用VASP+VTST计算Raman光谱（抛砖引玉）

最新推荐文章于 2024-09-20 20:14:50 发布

zesuzesuzesu

最新推荐文章于 2024-09-20 20:14:50 发布

阅读量6.7k

点赞数 7

分类专栏： vasp 文章标签：量子力学

本文链接：https://blog.csdn.net/zesuzesuzesu/article/details/120360014

版权

本文详细介绍了如何利用VASP和VTST工具进行Raman光谱的计算，包括所需软件、预设条件、计算步骤和结果处理。通过 phonon 计算、dymmatrix.py 脚本的应用，以及vasp_raman.py代码的调整，最终得到vasp_raman.dat文件，实现与实验数据的对比。作者呼吁有经验的专家提供指导，探讨ENCUT、温度等因素对结果的影响，并提及了使用vasp+phonopy+phono3py的另一种计算方法。

摘要由CSDN通过智能技术生成

鉴于目前没有找到很好的，很完整的教程由A到Z介绍如何用VASP计算Raman光谱。可能有干货的大神都留着自己用

摸索了一会儿找到了以下其中一个☝️方法，如果大神有补充或者纠正，欢迎留言

pre-requests:

VASP 和 VTST 安装好，没有的可以参照一下这篇文章
将会用到VTST scripts里面这两个code：dymselsph.pl, dymmatrix.py
会用最基本的vasp
安装好python
把github上的raman-sc下载下来，用python3.6的可以用知乎上这位大神的code。或者其实就是把原作者的code，vasp_raman.py里面的print都加上括号就可以用了
最重要的是根据原作者issues里面的这个问答，加上：

except:
    from lxml import etree
    doc = etree.parse('vasprun.xml')
    values = [c.text for c in doc.xpath("/modeling/calculation/varray"[3].getchildren()]



# 修改以后如下
if "MACROSCOPIC STATIC DIELECTRIC TENSOR" in line:
    outcar_fh.readline()
    try:
        epsilon.append([float(x) for x in outcar_fh.readline().split()])
        epsilon.append([float(x) for x in outcar_fh.readline().split()])
        epsilon.append([float(x) for x in outcar_fh.readline().split()])
    except:
        #以下是预防在OUTCAR里找不到MACROSCOPIC STATIC DIELECTRIC TENSOR
        #或者有时候OUTCAR里出现了一连串超大数字和负号， 把两个数连起来了。例如
        #1287356918273.125234-1234