fix ID group-ID phonon N Noutput Nwait map_file prefix keyword values
## 比如 每10step测量一组格林函数，测了每5000组后输出一次动力学矩阵，
## 等待50w个步长后开始测量
fix 1 all phonon 10 5000 500000 GAMMA EAM0D nasr 100- N -> 每隔这么多时间步长测量一次格林函数
- Noutput -> 每测量这么多次后输出一次动力学矩阵
- Mwait -> 等待这么多个时间步长后开始测量
- map_file -> file or GAMMA
在lammps的安装文件夹下面有examples，在PACKAGES/phonon里。整个MD计算需要这4个文件：
- in.EAM3D -> lammps 输入脚本
- data.pos -> 结构文件，一般是超胞结构
- map.in -> 超胞中原子的序号对应情况
- cuu3.eam -> Cu的力场文件（eam/opt力场）
最后得到的结果
