2. 偏振的数学描述

在偏振态制备与测量系统中,仅用“水平偏振”“圆偏振”这类定性说法无法支撑器件选型、光路计算和测量反演。工程上必须把偏振态写成可运算的数学对象。本节围绕两套最常用的描述体系展开:琼斯矢量适合描述完全偏振光在光学元件中的传播与变换,斯托克斯参数适合描述可含部分偏振的实际光束并直接对应可测光强,而庞加莱球则把斯托克斯参数映射为几何位置,使偏振态的演化、转换和误差变得直观可读。

一、琼斯矢量:用复振幅描述完全偏振态

先约定坐标系:光沿 +z 方向传播,x 轴为水平方向,y 轴为垂直方向,观察者沿传播方向看(即迎着光看)。任意完全偏振光的横向电场可写成两个正交分量的叠加:

E = [Ex; Ey] = [Ax·e^(jφx); Ay·e^(jφy)]

提取公共相位后,琼斯矢量通常归一化为:

J = [cosθ; sinθ·e^(jδ)],δ = φy − φx

其中 θ 决定两分量的振幅比,δ 是 y 相对 x 的相位差。几个基准态可直接写出:水平偏振 J_H = [1; 0],垂直偏振 J_V = [0; 1],45° 线偏振 J_45 = (1/√2)[1; 1],右旋圆偏振 J_R = (1/√2)[1; −j],左旋圆偏振 J_L = (1/√2)[1; +j]。注意圆偏振的符号约定与“迎着光看时旋转方向”相关,不同教材可能相反,工程文档中必须写明约定,本文统一采用上述形式。

琼斯矢量的价值在于变换是线性的:偏振元件用 2×2 琼斯矩阵表示,输出态 J_out = M·J_in。例如理想水平偏振片矩阵为 [[1,0],[0,0]],快轴在水平方向的四分之一波片矩阵为 [[1,0],[0,j]]。若输入 45° 线偏振 (1/√2)[1;1],经该波片后得到 (1/√2)[1;j],即左旋圆偏振——这正是“线偏振→圆偏振”制备的典型计算路径。实际使用时要注意:波片矩阵必须建立在波片快轴与 x 轴重合的前提下,若快轴旋转角为 α,需用旋转矩阵 R(−α)·M·R(α) 变换,否则相位差会算错。

二、斯托克斯参数:用可测光强描述任意偏振态

琼斯矢量只适用于完全偏振光,而实际激光器输出、经粗糙表面散射或长距离光纤传输后的光往往含非偏振成分。斯托克斯参数用四个实数量描述:

S0 = Ix + Iy           总光强
S1 = Ix − Iy           水平/垂直线偏振差
S2 = I45 − I−45         45°/135° 线偏振差
S3 = IR − IL           右旋/左旋圆偏振差

对完全偏振光,S0² = S1² + S2² + S3²;对部分偏振光,S0² > S1² + S2² + S3²,定义偏振度 DoP = √(S1²+S2²+S3²)/S0。工程上常把 S1、S2、S3 归一化为 s1、s2、s3(除以 S0),完全偏振态对应单位球面上的点。

斯托克斯参数的最大优势是每一项都对应一次可直接测量的光强组合。典型测量配置为:在探测器前依次插入水平偏振片、45° 偏振片、右旋圆偏振分析器(四分之一波片快轴与水平成 45° 再加水平偏振片),分别记录 I_H、I_V、I_45、I_−45、I_R、I_L,再按上式组合。例如测得 I_H = 0.80 mW、I_V = 0.20 mW、I_45 = 0.65 mW、I_−45 = 0.35 mW、I_R = 0.55 mW、I_L = 0.45 mW(示例数据),则 S0 = 1.00 mW,S1 = 0.60,S2 = 0.30,S3 = 0.10,DoP ≈ 0.68,说明这是一束部分偏振光,且线偏振成分占主导、方向偏向水平。

三、庞加莱球:把偏振态放到球面上

将归一化斯托克斯参数 (s1, s2, s3) 作为三维直角坐标,就得到庞加莱球:球心为完全非偏振态,球面为完全偏振态,球内为部分偏振态,半径即 DoP。赤道上 s3 = 0,对应所有线偏振态,s1 轴正方向为水平偏振、负方向为垂直偏振,s2 轴正方向为 45° 线偏振;北极 s3 = +1 为右旋圆偏振,南极 s3 = −1 为左旋圆偏振。波片的作用相当于绕某条过球心的轴旋转球面上的点:四分之一波片把赤道上的线偏振点转到极点,半波片则把点绕轴翻转 180°。这一几何图像让“波片快轴角度决定旋转轴”变得一目了然,也便于判断补偿方案。

庞加莱球:斯托克斯参数与偏振态的几何映射 s₁ s₃ s₂ H(水平) V(垂直) R(右旋圆) L(左旋圆) 45° 线偏振 球心:非偏振 部分偏振态(球内,半径=DoP) 两套描述的对应 琼斯矢量 J = [cosθ; sinθ·e^(jδ)] → 仅描述完全偏振光 → 线性变换 J_out = M·J_in 斯托克斯 S = (S₀,S₁,S₂,S₃) → 可含部分偏振,DoP 可算 → 每项对应可测光强组合 完全偏振时二者可互转: s₁ = cos2θ,s₂ = sin2θ·cosδ s₃ = sin2θ·sinδ 波片作用 = 球面上绕轴旋转

四、相互转换与工程使用要点

完全偏振态下,琼斯矢量与斯托克斯参数可以严格互转。由 J = [cosθ; sinθ·e^(jδ)] 可得 s1 = cos2θ,s2 = sin2θ·cosδ,s3 = sin2θ·sinδ;反向由 (s1, s2, s3) 求 θ = 0.5·arccos(s1),δ = atan2(s3, s2)。例如 J = (1/√2)[1; j] 对应 θ = 45°、δ = 90°,于是 s1 = 0,s2 = 0,s3 = 1,即北极右旋圆偏振,与几何图像一致。若测得 s1 = 0.6、s2 = 0.3、s3 = 0.1(示例),可先算 DoP ≈ 0.68,再归一化得 (0.88, 0.44, 0.15),据此定位球面上的点并反推主导偏振方向。

实际系统设计中有三点值得强调。第一,测量前必须标定分析器的角度与波片相位延迟:偏振片透光轴角度误差 1° 会引入约 3% 的 S1/S2 串扰,四分之一波片延迟偏离 90° 会直接污染 S3。第二,符号与坐标系约定要贯穿全链路:从琼斯矩阵建模、波片安装角度到斯托克斯反演,若圆偏振手性约定不一致,S3 会整体反号,导致补偿方向判断错误。第三,部分偏振光不能直接套用琼斯矩阵,应改用穆勒矩阵与斯托克斯矢量,或先通过偏振分析仪分解出完全偏振分量再处理。

快速自查:若一组斯托克斯参数算得 S0² 明显大于 S1²+S2²+S3²,说明光束含非偏振成分,此时用琼斯矢量描述会丢失信息;若接近相等,则可放心用琼斯矩阵做器件级光路计算,并用庞加莱球核对旋转结果。

把两套描述配合使用,是偏振系统设计的常见做法:用琼斯矢量完成波片、偏振分束器等器件的定量计算与制备方案设计,用斯托克斯参数和庞加莱球完成测量反演、偏振度评估与误差可视化。二者通过上述转换公式衔接,构成从“设计”到“验证”的完整数学闭环。