1. 引言

在科学研究领域，偏振光已作为一门独立的光学学科。在生物、化学、医疗领域有着广泛的应用。一个普通的发光体是由大量的发光原子或分子组成的，每个发光原子每次发射一个持续时间为微秒级或更短的近似平面波的波列，每一个原子或分子都是一个随机的自发辐射过程，彼此没有联系，从而导致合成的振动方向光场矢量各不相同，相关时间大大缩短了。光波在介质的表面上发生了反射、透射或是散射时，其偏振态的改变与物质自身的性质有关，如目标的表面结构，材料属性等。光在物体内部传输时，在一些物体外部施加磁场、应力时，也会出现不同的偏振态形式。偏振态能够提供比非偏振技术更丰富的信息。根据光波偏振特性来研究目标识别的偏振特性，最初是应用于天体表面物质特性研究和空间天体探测方面，在天体物理学中，空间被测目标也有一定偏振特性。因此，偏振态可以识别目标的位置和材料属性，偏振遥感成像与分析技术已应用于大气和云层中溶胶状态的探测、资源调查等领域 [1]。对生物组织病变前后反馈的光偏振参数进行对比和分析，可以快速、准确的给病理信息。在医疗诊断方面，圆偏振光的传播方式和螺旋型的蛋白质等分子结构相吻合，透过人体时吸收效果较好，光偏振检测技术在影像处理方面也有应用 [2] [3] [4]，对偏振参数和光谱参数进行融合计算，在恶劣的环境下可以增强图像边缘信息，了解探测目标散射特点后，被侦查目标可以通过雷达识别，可最大限度的对敌方军事目标进行毁伤，并对敌方更有效的军事打击 [5] [6]。2010年，吕鑫，蒋卓勤 [7] 等人在现有的方案上推导出了新的斯托克斯计算公式，提出了一种能对高速变化的偏振态进行实时测量的新方案。近几年来，随着激光技术的发展，光学的应用产品越来越多，其中偏振光作为一种重要的测量手段更受人们的欢迎 [8] [9] [10]，可以在复杂的背景环境中检测出实用的信息。通过高速偏振态测量系统的理论计算，偏振态探测技术的研究对光通信、光传感等领域，有着重要的指导意义。

2. 理论计算

如果用偏振态测量装置测量自然景物反射光的偏振特性时，假设进入测量系统的光子集合表示为：

$A\left(u_1,u_2,u_3,u_4,\cdots \right)$ (1)

其中 $u_n$ 表示状态为

$u_n=a_n\exp \left(j{\phi }_{nx}\right)\hat{x}+b_n\exp \left(j{\phi }_{ny}\right)\hat{y}$ (2)

$u_n$ 为总的光子数量，其中 $\hat{x}$ 和 $\hat{y}$ 分别表示沿x轴和y轴的单位矢量，并且 $a_n^2+b_n^2=1$。

在大多数偏振文献中，实现偏振态测量的方案很多，图1~图4为一种典型的偏振态测量方法 [11]。

图1的四路输出测量系统的分光器是等光强分束器、并且偏振无关(保偏分束器)，因此入射光子经过分束器后，光子的状态保持原有的信息，对于四路分束器来讲，根据归一化要求，可得到各路光子状态的归一化系数为0.5。对于 ${\stackrel{\to }{u}}_n$ 状态的光子，在分束器各路的状态可以表示为：

${u^{\prime }}_n=0.5a_n\exp \left(j{\phi }_{nx}\right)\hat{x}+0.5b_n\exp \left(j{\phi }_{ny}\right)\hat{y}$ (3)

在图1的偏振测量装置中，根据(3)式，可得第1路光强为：

$i_{n1}=0.25\left(a_n^2+b_n^2\right)$ (4)

当测量系统入射的光子数为 $c_n$ 时，则第1路的光强为：

$I_{n1}=c_n{i}_{n1}=0.25c_n\left(a_n^2+b_n^2\right)$ (5)

在第2路偏振片之后，光子状态为：

${{\stackrel{\to }{u}}^{\prime }}_{n2}=0.5a_n\exp \left(j{\phi }_{nx}\right)\hat{x}$ (6)

则第二路的光强为：

$i_{n2}=0.25a_n^2$ (7)

入射 $u_n$ 状态的光子数为 $c_n$ 时，则第2路的光强为：

$I_{n2}=0.25c_n{a}_n^2$ (8)

对于第3路，根据偏振片透光轴与x轴夹角为 $\theta$，可得到偏振片后的光子状态为：

${u^{″}}_{n3}=0.5\left[a_n\cos \theta \exp \left(j{\phi }_{nx}\right)+b_n\sin \theta \exp \left(j{\phi }_{ny}\right)\right]{\hat{r}}_{\theta }$ (9)

在(9)式中， ${\hat{r}}_{\theta }$ 表示沿偏振片透光轴方向的单位矢量，其表达式为

${\hat{r}}_{\theta }=\cos \theta \hat{x}+\sin \theta \hat{y}$ (10)

当第3路的偏振片方位角 $\theta ={45}^{\circ }$ 时，透过偏振片的光子状态为：

${u^{″}}_{n3}=0.25\sqrt{2}\left[a_n\exp \left(j{\phi }_{nx}\right)+b_n\exp \left(j{\phi }_{ny}\right)\right]{\hat{r}}_{\theta }$ (11)

第3路偏振片后的光强为：

$i_{n3}=0.125\left[a_n^2+b_n^2+2a_n{b}_n\cos {\delta }_n\right]$ (12)

在(12)式中 ${\delta }_n={\phi }_{ny}-{\phi }_{nx}$。

当入射光子数为 $c_n$ 时，则偏振片后光强为：

$I_{n3}=c_n{i}_{n3}=0.125c_n\left[a_n^2+b_n^2+2a_n{b}_n\cos {\delta }_n\right]$ (13)

对于第4路，波片后的光子状态为：

${{\stackrel{\to }{u}}^{\prime }}_{n4}=0.5a_n\exp \left(j{\phi }_{nx}\right)\hat{x}+0.5b_n\exp \left[j\left({\phi }_{ny}+0.5\pi \right)\right]\hat{y}$ (14)

可得到偏振片后的光子状态为：

${u^{″}}_{n4}=\left\{0.5a_n\cos \theta \exp \left(j{\phi }_{nx}\right)+0.5b_n\sin \theta \exp \left[j\left({\phi }_{ny}+0.5\pi \right)\right]\right\}{\hat{r}}_{\theta }$ (15)

当第4路的偏振片方位角 $\theta ={45}^{\circ }$ 时，偏振片后的光子状态为：

${u^{″}}_{n4}=0.25\sqrt{2}\left\{a_n\exp \left(j{\phi }_{nx}\right)+b_n\exp \left[j\left({\phi }_{ny}+0.5\pi \right)\right]\right\}{\hat{r}}_{\theta }$ (16)

此时光强为：

$\begin{array}{c}{i}_{n4}=0.125\left[a_n^2+b_n^2+2a_n{b}_n\cos \left({\phi }_{ny}+0.5\pi -{\phi }_{nx}\right)\right]\\ =0.125\left[a_n^2+b_n^2-2a_n{b}_n\sin \left({\phi }_{ny}-{\phi }_{nx}\right)\right]\\ =0.125\left[a_n^2+b_n^2-2a_n{b}_n\sin {\delta }_n\right]\end{array}$ (17)

当入射光子数为 $c_n$ 时，则加偏振片后光强为：

$I_{n4}=0.125c_n\left[a_n^2+b_n^2-2a_n{b}_n\sin {\delta }_n\right]$ (18)

因为所有的光子状态相同，通过 $I_{n1}$ 、 $I_{n2}$ 、 $I_{n3}$ 和 $I_{n4}$ 的结果得出光子集合态的斯托克斯参量。并给出每个斯托克斯参量的计算公式。

$s_{n0}=I_{n1}=0.25c_n\left(a_n^2+b_n^2\right)$ (19)

$s_{n1}=2I_{n2}-I_{n1}=0.5c_n{a}_n^2-0.25c_n\left(a_n^2+b_n^2\right)=0.25c_n\left(a_n^2-b_n^2\right)$ (20)

$s_{n2}=2I_{n3}-I_{n1}=0.25c_n\left(a_n^2+b_n^2+2a_n{b}_n\cos {\delta }_n\right)-0.25c_n\left(a_n^2+b_n^2\right)=0.5c_n{a}_n{b}_n\cos {\delta }_n$ (21)

$s_{n3}=I_{n1}-2I_{n4}=0.25c_n\left(a_n^2+b_n^2\right)-0.25c_n\left(a_n^2+b_n^2-2a_n{b}_n\sin {\delta }_n\right)=0.5c_n{a}_n{b}_n\sin {\delta }_n$ (22)

3. 光子集合与Stokes关系

如果进入偏振光测量系统的光子集合为 $A\left(c_1,c_2,c_3,\cdots \right)$ 的情况时，各光探测器输出电流分别为：

$I_1={\displaystyle {\sum }_n{I}_{n1}}=0.25{\displaystyle {\sum }_n{c}_n\left(a_n^2+b_n^2\right)}$ (23)

$I_2={\displaystyle {\sum }_n{I}_n{}_2}=0.25{\displaystyle {\sum }_n{c}_n{a}_n^2}$ (24)

$I_3={\displaystyle {\sum }_n{I}_{n3}}=0.125{\displaystyle {\sum }_n{c}_n\left[a_n^2+b_n^2+2a_n{b}_n\cos {\delta }_n\right]}$ (25)

$I_4={\displaystyle {\sum }_n{I}_{n4}}=0.125{\displaystyle {\sum }_n{c}_n\left[a_n^2+b_n^2-2a_n{b}_n\sin {\delta }_n\right]}$ (26)

(22)式至(25)式探测器输出信号 $I_1$ 、 $I_2$ 、 $I_3$ 和 $I_4$ 为光子集合的贡献。此时系统给出的斯托克斯参量为

$S_0={\displaystyle {\sum }_n{I}_{n1}}=0.25{\displaystyle {\sum }_n{c}_n\left(a_n^2+b_n^2\right)}={\displaystyle {\sum }_n{s}_{n0}}$ (27)

$S_1={\displaystyle {\sum }_n\left(2I_{n2}-I_{n1}\right)}=0.25{\displaystyle {\sum }_n{c}_n\left(a_n^2-b_n^2\right)}={\displaystyle {\sum }_n{s}_{n1}}$ (28)

$S_2={\displaystyle {\sum }_n\left(2I_{n3}-I_{n1}\right)}=0.5{\displaystyle {\sum }_n{c}_n{a}_n{b}_n\cos {\delta }_n}={\displaystyle {\sum }_n{s}_{n2}}$ (29)

$S_3={\displaystyle {\sum }_n\left(I_{n1}-2I_{n4}\right)}=0.5{\displaystyle {\sum }_n{c}_n{a}_n{b}_n\sin {\delta }_n}={\displaystyle {\sum }_n{s}_{n3}}$ (30)

由于系统输入的是光子集合，无法利用(26)~(29)结果得出集合中光子的斯托克斯参量。

4. 结论

根据以上推导分析，可分两种情况进行解释不同集合光子的偏振态的叠加。假设第一种情况是 $c_1$ 个 $u_1=a_1\exp \left(j{\phi }_{1x}\right)\hat{x}+b_1\exp \left(j{\phi }_{1y}\right)\hat{y}$ 状态的光子和 $c_2$ 个 $u_2=a_2\exp \left(j{\phi }_{2x}\right)\hat{x}+b_2\exp \left(j{\phi }_{2y}\right)\hat{y}$ 状态的光子，利用(26)-(29)式得出斯托克斯参量为：

${S^{\prime }}_0={\displaystyle {\sum }_n{s}_{n0}}=0.25c_1\left(a_1^2+b_1^2\right)+0.25c_1\left(a_2^2+b_2^2\right)$ (31)

${S^{\prime }}_1={\displaystyle {\sum }_n{s}_{n1}}=0.25c_1\left(a_1^2-b_1^2\right)+0.25c_1\left(a_2^2-b_2^2\right)$ (32)

${S^{\prime }}_2={\displaystyle {\sum }_n{s}_{n2}}=0.5c_1{a}_1{b}_1\cos {\delta }_1+0.5c_2{a}_2{b}_2\cos {\delta }_2$ (33)

${S^{\prime }}_3={\displaystyle {\sum }_n{s}_{n3}}=0.5c_1{a}_1{b}_1\sin {\delta }_1+0.5c_2{a}_2{b}_2\sin {\delta }_2$ (34)

如第二种情况是 $c_3$ 个 $u_3=a_3\exp \left(j{\phi }_{3x}\right)\hat{x}+b_3\exp \left(j{\phi }_{3y}\right)\hat{y}$ 状态的光子和 $c_4$ 个 $u_4=a_4\exp \left(j{\phi }_{4x}\right)\hat{x}+b_4\exp \left(j{\phi }_{4y}\right)\hat{y}$ 状态的光子，利用(26)~(29)式得出斯托克斯参量为：

${S^{″}}_0={\displaystyle {\sum }_n{s}_{n0}}=0.25c_3\left(a_3^2+b_3^2\right)+0.25c_4\left(a_4^2+b_4^2\right)$ (35)

${S^{″}}_1={\displaystyle {\sum }_n{s}_{n1}}=0.25c_3\left(a_3^2-b_3^2\right)+0.25c_4\left(a_4^2-b_4^2\right)$ (36)

${S^{″}}_2={\displaystyle {\sum }_n{s}_{n2}}=0.5c_3{a}_3{b}_3\cos {\delta }_3+0.5c_4{a}_4{b}_4\cos {\delta }_4$ (37)

${S^{″}}_3={\displaystyle {\sum }_n{s}_{n3}}=0.5c_3{a}_3{b}_3\sin {\delta }_3+0.5c_4{a}_4{b}_4\sin {\delta }_4$ (38)

针对光学偏振态的测试需求，对光子集合偏振态的理论研究，推导出光子集合与斯托克斯矢量的关系以及相应的计算公式，得到光子集合偏振态的斯托克斯四个矢量参数，如果适当选择 $c_1~c_4$ 、 $a_1~a_4$ 、 $b_1~b_4$ 、 ${\delta }_1~{\delta }_4$ 这些参数，可以使得 ${S^{″}}_0={S^{\prime }}_0$ 、 ${S^{″}}_1={S^{\prime }}_1$ 、 ${S^{″}}_2={S^{\prime }}_2$ 和 ${S^{″}}_3={S^{\prime }}_3$。虽然光子集合偏振态的构成有所不同，但是测量装置给出的结果却是相同的，这就是目前偏振态测量装置的特性。理论推导出光子集合偏振态与斯托克斯矢量的函数关系，对今后的偏振态测量技术有一定的指导意义。

NOTES

^*通讯作者。

参考文献