VASP非共振拉曼张量
用硅光学模位移前后的介电响应，提取静态拉曼张量与相对活性。
阅读本篇目录
振动频率告诉我们谱峰可能出现在哪里，拉曼活性还要看振动怎样改变体系的电子极化响应。本例用Si原胞在Γ点的光学振动，展示“声子→正负位移→介电响应→拉曼张量”的流程。
1. 先准备振动模式#
沿用晶体声子教程的实际Si力常数。原胞含2个Si、体积40.892732 ų，Γ点三个光学模简并，频率约502.798 cm⁻¹。
在这个立方晶体的简并子空间中，可选两个Si沿x、y、z方向分别反向运动作为三个模式。简并空间可以换基，不能将其中一个基底方向视作唯一不可改变的本征矢。
本例Si质量28.0855 amu。若两个原子分别移动+δ和−δ，则对应的质量加权模式坐标Q=δ√(2m)。忽略质量归一化会使不同元素或不同原胞的强度单位无法比较。
2. 对正负位移计算电子介电张量#
每个方向采用δ=0.01 Å，另外对x方向补δ=0.02 Å以检查差分步长。所有位移结构固定原子，计算电子响应，不再优化回原来的结构。
LEPSILON = .TRUE.
IBRION = -1
NSW = 0
EDIFF = 1E-6
ENCUT = 550
本例6×6×6 Gamma网格、PBE PAW64、NBANDS=16。使用VASP 6.6.1 CPU分支vasp/6.6.1-8173m，经Slurm在兼容的V100宿主CPU上以4个MPI进程运行。它是CPU微扰响应计算，不作为GPU拉曼加速示例；CPU版使用oneAPI 2024.2；进一步构建信息与极化和介电教程一致。
在OUTCAR读取包含局域场效应的MACROSCOPIC STATIC DIELECTRIC TENSOR。本次未位移参考的对角分量约16.4115，三个方向一致；该值对应本页较粗的采样，未宣称是收敛的Si介电常数。VASP微扰介电响应
3. 从响应差分得到模式张量#
采用每原胞极化率体积形式α=V(ε−I)/(4π)，得到：
R = ∂α/∂Q
≈ V/(4π) × [ε(+δ)−ε(−δ)] / [2δ√(2m)]
Q以Å√amu计，R以Ų/√amu计。脚本保留原始张量反对称部分的数值大小，再使用对称部分计算非共振活性。
沿x反向运动时，主要响应出现在yz、zy分量；另外两个模式为循环置换。这展示了模式和偏振方向之间的关系，不能只凭频率相同就认为各实验偏振下强度都相同。
4. 活性与差分步长#
取a′=Tr(R)/3，γ′²为张量各向异性组合，常用取向平均活动量为S=45a′²+7γ′²。它还没有包含激光频率、温度因子、折射与实验几何等全部散射强度因素。
| 模式方向 | 每原子位移δ / Å | 频率 / cm⁻¹ | S / Å⁴·amu⁻¹ |
|---|---|---|---|
| x | 0.010 | 502.798 | 2687.553 |
| x | 0.020 | 502.798 | 2718.456 |
| y | 0.010 | 502.798 | 2687.573 |
| z | 0.010 | 502.798 | 2687.552 |
三个δ=0.01 Å模式给出接近的活性；x方向把位移加倍后，活性变化约1.15%。这个对照检查差分步长的影响，不能替代k网格和电子响应收敛检查。
图中的谱峰是把三个简并模式相加后，使用人为高斯展宽展示；展宽不是本次计算得到的声子寿命或实验线宽。
5. 适用范围#
本页使用静态电子响应，属于非共振近似。若激光能量接近电子跃迁、需要共振拉曼或复杂电子–声子耦合，不能直接把这套静态结果当作完整光谱。金属、极性材料以及二维体系还需分别确认响应和体积归一化。
输入与提取#
按README提交参考和8个位移目录，再运行python extract_raman.py。脚本读取OUTCAR张量，结合附件中已实算的Γ点振动参照，输出RAMAN-DATA.json。换结构时必须重新取得同一晶胞的振动模式，不能沿用本例的频率和反向位移。