Phonopy-Spectroscopy 拉曼偏振操作

偏振拉曼需要指定表面、光的偏振和转角。同一组三重简并模式，在不同实验几何下可以出现完全不同的峰强。

下载完整光谱实例包。在新目录解压后操作；包内包含真实输入、已有DFT响应、25个操作条件、原生输出、数据来源说明和绘图脚本。

1. 先准备实际拉曼张量#

在新目录依次执行：

bash prepare.sh
bash raman.sh
bash polarized.sh

前两步完成Si模式和响应读取；最后一步生成平行、交叉及指定角度的结果。

2. 表面指数采用哪个晶胞#

本例使用FCC原胞基矢，配套参数为 --surface "0 1 1"，其法向对应常规立方晶胞的[100]。不能直接将原胞Miller(011)称作常规晶胞(011)。脚本用倒易晶格独立核对这个映射。

偏振向量在仪器坐标系中定义：表面法向旋转到入射轴方向，平行偏振取相同的入射/出射电场方向，交叉偏振取正交方向。

3. 使用配套偏振入口#

需要手工指定参数时，在包含已读取 Raman.yaml 的目录中加载本专题环境，再调用包内入口：

python /你的实例包路径/polarized_raman.py \
  --surface "0 1 1" --theta "0 45 90" --scattered-pol parallel

完整测试条件通过短入口选择：

bash run.sh si-polar-parallel
bash run.sh si-polar-cross
bash run.sh si-polar-angles

完整转角条件从0°到360°、每5°一个点；指定角度条件为0、22.5、45、67.5、90°。不要在已经由 polarized.sh生成同名结果的目录重复运行这些条件，换新目录复现。

4. 变换与投影#

配套入口使用标准二阶张量变换：

T′ = R · T · Rᵀ
I ∝ |e_inᵀ · T′ · e_out|²

固定源码提交的原始偏振分支采用了不同的旋转乘法；本例通过配套入口在当前进程中使用上式，程序安装文件保持原样。入口同时整理角度CSV的模式列名，输出Mode4/5/6。

这一几何下，主要模式的平行与交叉响应交替增强。配套分析同时检查张量迹、Frobenius范数和独立投影，避免只凭曲线看起来平滑就判断变换正确。

5. 读输出#

head work/si-polar-parallel/Raman-Pol-IntensityTheta.csv

第一列角度，后面三列对应Mode4/5/6。指定角度条件另输出按角度排列的峰表和旋转后张量；完整转角条件还输出二维谱线图。

固定几何投影与非偏振活动量采用不同的几何系数，不直接比较两者峰高。弱小分量还会受有限差分与电子响应数值精度影响。