#LAMMPS计算二阶力矩阵与声子谱

1.安装phonolammps

pip install phonolammps

2.准备lammps输入文件
本文用graphene作为例子

• in.lammps LAMMPS输入文件

units    metal
atom_style  atomic
dimension 3
boundary p p p
read_data     graphene.data #读取结构
pair_style    nep nep.txt  #选择自己的势函数
pair_coeff    * *
#pair_style    airebo 3.0
#pair_coeff    * * CH.airebo C

• graphene.data 结构文件
• nep.txt 势函数文件

3.利用phonolammps生成二阶力矩阵

phonolammps in.lammps -c POSCAR --dim  3 3 1

运行完之后会生成
FORCE_CONSTANTS
这个就是二阶力矩阵

安装lammps python接口
运行如下代码：
make serial mode=shlib
经过较长时间的编译，生成了liblammps.so文件。再运行
make install-python
则安装成功。
在python环境下运行
from lammps import lammps
成功。

5.利用vaspkit生成高对称点

vaspkit #运行vaspkit 输入305 2生成二维材料的高对称点
305 
2

cp KPATH.phonopy band.conf
vi band.conf

NPOINTS = 501
DIM =  2 2 1
BAND = 0.000000 0.000000 0.000000 0.500000 -0.500000 0.000000 0.625456 -0.374544 0.000000 0.500000 0.000000 0.000000 0.374544 0.374544 0.000000 0.000000 0.500000 0.000000 0.000000 0.000000 0.000000
BAND_LABELS = $\Gamma$ X H$_1$ C H Y $\Gamma$

MP = 30 30 30
TETRAHEDRON = .TRUE.
#PDOS = 1 2 3 4 5 6 7 8 9 10 11 12
BAND_CONNECTION = .TRUE.
FORCE_CONSTANTS = READ

# FORCE_SETS = READ
# IRREPS = 0  0  0
# SHOW_IRREPS = .TRUE.
# LITTLE_COGROUP = .TRUE.

phonopy -c POSCAR band.conf --dim=" 3 3 1" -p -s #绘图

文件下会产生band-dos.pdf
也可以自已拿数据画

phonopy-bandplot --gnuplot > phonon.out #导出数据自己画