2. 偏振光学基础：光的偏振态描述、琼斯矩阵与穆勒矩阵、偏振测量原理

做应力双折射，说白了就是在跟光的偏振态打交道。你连光怎么振、怎么描述、怎么测都搞不清楚，后面分析数据就是瞎蒙。这一章我把偏振光学里最核心的东西拎出来，全是干活用得上的。

2.1 光的偏振态——光矢量到底在怎么动？

光是一种横波，电场矢量 E 的振动方向垂直于传播方向。这个振动方向怎么随时间变化，就是偏振态。我个人习惯把偏振态分成三类：

• 线偏振光：电场矢量只在一个固定方向上振动。比如你透过偏振片看天空，转一下角度，亮度会变，那就是线偏振光在起作用。
• 圆偏振光：电场矢量端点画出一个圆。左旋还是右旋？用右手定则，拇指指向光传播方向，四指弯曲方向跟电场旋转方向一致，就是右旋。我刚开始学的时候老搞反，后来干脆记住：对着光源看，顺时针转就是右旋。
• 椭圆偏振光：最一般的情况。线偏振和圆偏振都是它的特例。应力双折射样品出来的光，绝大多数是椭圆偏振光。

2.2 琼斯矩阵——处理完全偏振光的利器

琼斯矩阵用一个 2×1 的列向量表示光的偏振态，用 2×2 的矩阵表示光学元件。我特别喜欢它的简洁：

琼斯矢量：E = [Ex; Ey] = [Ax * exp(iφx); Ay * exp(iφy)]

琼斯矩阵：J_out = M * J_in

举个例子，水平线偏振光通过一个快轴水平的 1/4 波片，会变成右旋圆偏振光。用琼斯矩阵算一下：

J_in = [1; 0]                     // 水平线偏振
M_qwp = [1, 0; 0, -i]            // 快轴水平，相位延迟 π/2
J_out = M_qwp * J_in = [1; -i]    // 右旋圆偏振

嗯，这里要注意：琼斯矩阵只能处理完全偏振光。如果你遇到部分偏振光或者非偏振光，就得请出穆勒矩阵了。

2.3 穆勒矩阵——通吃所有偏振态

穆勒矩阵用 4×1 的斯托克斯矢量表示光，用 4×4 的矩阵表示元件。斯托克斯矢量的四个参数是：

参数	含义	物理意义
S₀	总光强	I₀
S₁	水平 vs 垂直线偏振分量	I₀·cos(2χ)·cos(2ψ)
S₂	+45° vs -45° 线偏振分量	I₀·cos(2χ)·sin(2ψ)
S₃	右旋 vs 左旋圆偏振分量	I₀·sin(2χ)

其中 ψ 是椭圆方位角，χ 是椭圆率角。你想想看，有了这四个数，任何偏振态都能唯一确定。

2.4 偏振测量原理——怎么把偏振态测出来？

测量偏振态，核心思路就一句话：让光通过一系列已知的偏振元件，测出光强变化，反推出入射光的偏振态。

常用的方法有：

1. 旋转波片法：让 1/4 波片匀速旋转，后面放一个固定偏振片，探测器记录光强随时间的变化曲线。然后做傅里叶分析，提取出斯托克斯参数。我在项目中用过这个方法，优点是光路简单，缺点是旋转机构有机械误差。
2. 液晶可变延迟器法：用电控液晶代替旋转波片，没有运动部件，速度快。我记得第一次用这个方案时，被液晶的响应时间和温度漂移折腾得不轻。
3. 分振幅法：把入射光分成四路，每路经过不同的偏振分析器，同时测量四个光强，直接算出斯托克斯参数。实时性好，但光路复杂，对校准要求高。

下面这张图展示了旋转波片法偏振测量系统的核心逻辑：

探测器记录的光强 I(θ) 与波片旋转角 θ 的关系是：

I(θ) = 0.5 * [S₀ + S₁·cos²(2θ) + S₂·sin(2θ)·cos(2θ) - S₃·sin(2θ)]

对测得的 I(θ) 做傅里叶级数展开，提取出直流分量、2θ 分量和 4θ 分量的系数，就能反解出 S₀、S₁、S₂、S₃。我在项目里一般采集 200 个角度点的数据，然后用最小二乘法拟合，精度能到 0.5%。

2.5 偏振测量中的误差源

做偏振测量，误差来源比你想象的多。我踩过的坑列出来给你参考：

误差来源	典型影响	我的应对方法
波片相位延迟不准	圆偏振分量 S₃ 偏差	用标准偏振片+1/4波片做系统校准
偏振片消光比不足	线偏振分量 S₁、S₂ 偏差	选用消光比 > 10⁵ 的偏振片
探测器非线性	所有分量等比例偏差	用中性密度滤光片做线性度标定
环境杂散光	S₀ 偏大，其他分量信噪比下降	加光阑、做暗室测量、用锁相放大器
温度漂移	液晶/波片延迟量变化	恒温控制，或者实时监测温度做补偿

2.6 从偏振测量到应力双折射

好了，现在你能测出样品的穆勒矩阵了。怎么从中提取应力双折射信息？

对于无退偏的透明样品，穆勒矩阵可以分解为：

M = M_ret * M_diat

其中：
M_ret —— 相位延迟矩阵（包含双折射信息）
M_diat —— 二向色性矩阵（包含吸收各向异性信息）

对于应力双折射测量，我们主要关注 M_ret。从 M_ret 可以提取出两个关键参数：

• 相位延迟量 δ：光通过样品后，两个正交偏振分量的相位差。δ = 2π·Δn·d / λ，其中 Δn 是双折射率，d 是样品厚度，λ 是波长。
• 快轴方位角 φ：快轴在样品平面内的方向。

这两个参数，就是应力双折射分析的基础。下一章我会详细讲怎么从穆勒矩阵的测量数据中，稳定地提取出 δ 和 φ。