Phonopy-Spectroscopy 计算Si拉曼
按位移生成、读取真实介电响应、后处理三步操作，展示张量、活动量和差分步长。
阅读本篇目录
非共振拉曼需要模式位移引起的电子极化率变化。Phonopy-Spectroscopy的三个步骤分别为 -d生成位移、-r读取介电响应、-p计算活动量和光谱。
下载完整光谱实例包。在新目录解压后操作；包内包含真实输入、已有DFT响应、25个操作条件、原生输出、数据来源说明和绘图脚本。
1. 先看本例结构与模式#
Si原胞2原子，Γ点三个光学模简并。这里在实际简并子空间中选用实数笛卡尔基：两个Si沿x、y、z分别反向运动。脚本先检查该基与实际动力学矩阵的本征子空间一致，再写入模式文件；频率和力常数保持原值。
参考原胞：
Si
1.0000000000000000
0.0000000000000000 2.7344629562984917 2.7344629562984917
2.7344629562984917 0.0000000000000000 2.7344629562984917
2.7344629562984917 2.7344629562984917 0.0000000000000000
Si
2
Direct
0.0000000000000000 0.0000000000000000 -0.0000000000000000
0.2500000000000000 0.2500000000000000 0.2500000000000000
简并模式可以换基，单个模式的偏振响应要结合所选基底解释。
2. 生成位移#
bash prepare.sh
bash run.sh si-displacements
程序在 work/si-displacements/ 生成六份POSCAR，每条光学模式有负、正两个结构，没有零位移结构。底层命令：
phonopy-raman -d --bands "4 5 6" --amplitude 0.01 --step-type max_r
max_r表示每个结构中最大原子位移0.01 Å；模式4/5/6的序号从1开始。本例每个Si都移动0.01 Å，方向相反。它对应约0.0749473 Å√amu的质量加权模式步长。
3. 对位移结构计算电子介电张量#
实例包已提供在集群完成的真实VASP响应。自己的新结构应逐个固定构型计算，不要让位移结构重新优化回参考位置。
本例完整INCAR
SYSTEM = Si finite-field Raman reference
ENCUT = 550
PREC = Accurate
EDIFF = 1E-6
NELM = 200
ALGO = Normal
ISMEAR = 0
SIGMA = 0.05
ISYM = 0
NSW = 0
IBRION = -1
LREAL = .FALSE.
LASPH = .TRUE.
LWAVE = .FALSE.
LCHARG = .FALSE.
LEPSILON = .TRUE.
NBANDS = 16
实际KPOINTS：
Gamma
0
Gamma
6 6 6
0 0 0
关键是 LEPSILON、固定原子以及相同的电子设置。读取电子介电响应时使用包括局域场效应的矩阵，不能把离子贡献也混入这里的差分。已有DFT准备方法见VASP非共振拉曼。
4. 读取响应并画谱#
bash run.sh si-read
bash run.sh si-raman-L5
读取入口先逐份核对生成POSCAR与真实DFT结构的周期几何，再按正确顺序将原始OUTCAR传给 phonopy-raman -r。原生读取过程会更新 Raman.yaml，运行器保留读取前的副本。
-p输出 Raman-Tensors.yaml、Raman-PeakTable.dat 和 Raman-Spectrum.dat/png。
5. 结果和单位#
三个光学模式频率为502.798114 cm⁻¹，活动量分别为2687.553003、2687.573587、2687.552659 Å⁴/amu。张量按 R=V/(4π) × Δε/ΔQ 得到，物理单位为 Ų/√amu。
当前原生张量YAML的 activity_units文字把质量平方根的位置写得不准确；本文按质量加权定义解释数值。标量峰表的Å⁴/amu单位与上述定义一致。
6. 差分步长检查#
bash run.sh si-step02
这个条件只计算x方向0.02 Å的既有响应：活动量2718.456011 Å⁴/amu，相比0.01 Å变化约1.15%。改变位移幅度后必须有对应的DFT响应，不能只改YAML步长后复用不匹配的OUTCAR。
这些量属于静态电子响应下的非共振拉曼活动量，还未包含完整的激光频率、玻色占据和实验几何系数。