通过Python理解Jones矩阵

Jones矩阵

Jones向量

假设光波沿z轴传播，那么其三个方向的电场分量可以表示为

$$\begin{array}{rl}{E}_x& =A_x\cos (\omega t-\mathbf{k}\mathbf{r}+{\phi }_x)\\ E_y& =A_t\cos (\omega t-\mathbf{k}\mathbf{r}+{\phi }_y)\\ E_z& =0\end{array}$$

由于传播方向为0，故可以通过一个二维向量来表示

$$\left[\begin{array}{c}{E}_x\\ E_y\end{array}\right]=\left[\begin{array}{c}\text{Re}[A_x\exp (\text{j}{\phi }_x)\exp (\text{j}\omega t-\text{j}\mathbf{k}\mathbf{r})]\\ \text{Re}[A_y\exp (\text{j}{\phi }_y)\exp (\text{j}\omega t-\text{j}\mathbf{k}\mathbf{r})]\end{array}\right]$$

则记

$$\hat{J}=\left[\begin{array}{c}{A}_x\exp (\text{j}{\phi }_x)\\ A_y\exp (\text{j}{\phi }_y)\end{array}\right]==\left[\begin{array}{c}{J}_x\\ J_y\end{array}\right]$$

此即琼斯向量，其中包含了偏振态的全部信息，而且附赠了光波的振幅项$\mathbf{A}$。如果只关心偏振量，则可对$\hat{J}$进行归一化，则令

$$\tan \psi =\frac{A_y}{A_x}\exp [\text{j}({\phi }_y-{\phi }_x)]$$

可表示为

$$\hat{J}=\left[\begin{array}{c}\cos \psi \\ \sin \psi \end{array}\right]$$

当${\phi }_x={\phi }_y$时，表示x和y方向不存在相位差，此时$J_x,J_y$均为实数，描述的是线偏振光；当存在相位差时，琼斯向量中出现虚数，表示圆偏振光。

对于偏振光$\hat{J}$，与其正交的偏振向量可以表示为

$${\hat{J}}^{\ast }=\left[\begin{array}{c}-\sin \psi \\ \cos \psi \end{array}\right]$$

二者符合关系

$$\hat{J}\cdot {\hat{J}}^{\ast }=0$$

常见的琼斯向量有

向量	说明
${\left[\begin{array}{cc}1& 0\end{array}\right]}^T$	x方向线偏振
${\left[\begin{array}{cc}0& 1\end{array}\right]}^T$	y方向线偏振
$\frac{1}{\sqrt{2}}{\left[\begin{array}{cc}1& 1\end{array}\right]}^T$	45°偏振
$\frac{1}{\sqrt{2}}{\left[\begin{array}{cc}1& -1\end{array}\right]}^T$	-45°偏振
$\frac{1}{\sqrt{2}}{\left[\begin{array}{cc}\cos \theta & \sin \theta \end{array}\right]}^T$	$\theta$偏振
$\frac{1}{\sqrt{2}}{\left[\begin{array}{cc}1& -\text{i}\end{array}\right]}^T$	右旋圆偏振
$\frac{1}{\sqrt{2}}{\left[\begin{array}{cc}1& \text{i}\end{array}\right]}^T$	左旋圆偏振

Jones矩阵

能够保证二维列向量形状不变的运算有无穷多种，但最符合我们直觉的一定是$2\times 2$矩阵。好在这种矩阵已经可以提供足够多的运算，从而满足我们描述偏振变化的需求。

先考虑最简单的情况，一束偏振光经过x轴的线偏振片之后会有什么后果？无论入射光如何，其出射光必然遵循$[A0{]}^T$的形式，故对应的Jones矩阵可以表示为

$$\left[\begin{array}{cc}1& 0\\ 0& 0\end{array}\right]$$

对于与x轴夹角为$\Psi$的偏振片，其表达形式为

$$\frac{1}{2}\left[\begin{array}{cc}{\cos }^2\Psi & \cos \Psi \sin \Psi \\ \sin \Psi \cos \Psi & {\sin }^2\Psi \end{array}\right]$$

常见偏振片的Jones矩阵可以表示为

Jonoes矩阵	说明
$\left[\begin{array}{cc}1& 0\\ 0& 0\end{array}\right]$	x方向线偏振片
$\left[\begin{array}{cc}0& 0\\ 0& 1\end{array}\right]$	y方向线偏振片
$\frac{1}{2}\left[\begin{array}{cc}1& 1\\ 1& 1\end{array}\right]$	45°偏振片
$\frac{1}{2}\left[\begin{array}{cc}1& -1\\ -1& 1\end{array}\right]$	-45°偏振片
$\frac{1}{2}\left[\begin{array}{cc}1& \text{i}\\ -\text{i}& 1\end{array}\right]$	右旋偏振片
$\frac{1}{2}\left[\begin{array}{cc}1& -\text{i}\\ \text{i}& 1\end{array}\right]$	左旋偏振片

光在通过波片之后，会在不同方向产生差异性的相位延迟，对于与x轴角度为$\Psi$，相位差为$\Phi$的波片，其Jones矩阵为

$$\left[\begin{array}{cc}{e}^{-\text{i}\frac{\Phi }{2}}{\cos }^2\Psi +e^{\text{i}\frac{\Phi }{2}}{\sin }^2\Psi & -\text{i}\sin \frac{\Phi }{2}\sin 2\Psi \\ -\text{i}\sin \frac{\Phi }{2}\sin 2\Psi & e^{\text{i}\frac{\Phi }{2}{\sin }^2\Psi }+e^{\text{i}\frac{\Phi }{2}}{\cos }^2\Psi \end{array}\right]$$

一些常用的光学矩阵如下

Jonoes矩阵	说明
$\left[\begin{array}{cc}{e}^{-\text{i}\frac{\Phi }{2}}& 0\\ 0& e^{\text{i}\frac{\Phi }{2}}\end{array}\right]$	x轴方向相位差为$\Phi$的波片
$\left[\begin{array}{cc}-\text{i}& 0\\ 0& \text{i}\end{array}\right]$	x轴方向半波片
$\frac{1}{\sqrt{2}}\left[\begin{array}{cc}1-\text{i}& 0\\ 0& 1+\text{i}\end{array}\right]$	x轴方向$\lambda /4$波片

我们平时使用的偏振器件，多半放在可旋转的镜架上，对偏振器件进行旋转操作的矩阵不出意外地表示为

$$\left[\begin{array}{cc}\cos \theta & -\sin \theta \\ \sin \theta & \cos \theta \end{array}\right]$$

Jones矩阵的表示

为了对Jones矩阵所对应的偏振状态进行绘制，我们需要进一步理解Jones矩阵中每个值所对应的物理概念。如果将Jones矩阵写成虚数形式，可以表示为

$$J=\left[\begin{array}{c}{a}_x+i{b}_x\\ a_y+i{b}_y\end{array}\right]$$

其模矩阵表示x轴与y轴的相对强度，可以表示为

$$\left[\begin{array}{c}\sqrt{a_x^2+b_x^2}\\ \sqrt{a_y^2+b_y^2}\end{array}\right]$$

其俯角矩阵表示x轴与y轴的俯角则表示相位延迟，可以表示为

$$\left[\begin{array}{c}\arctan \frac{b_x}{a_x}\\ \arctan \frac{b_y}{a_y}\end{array}\right]$$

所以，对于右旋圆偏振光$\frac{1}{\sqrt{2}}{\left[\begin{array}{cc}1& -\text{i}\end{array}\right]}^T$而言，其x方向幅角为0，y方向幅角为$-\frac{\pi }{2}$，二者相差$\frac{\pi }{2}$，故为右旋圆偏振光。

那么，对于任意Jones矩阵，我们可以很方便地绘制其对应的偏振图像。方便起见，我们只对1.064um光波在z轴方向传播5um这段距离进行采样，然后画出

def drawJones(J=[1,-2j]):
    J = np.array(J).reshape(2,1)
    # 设1.064um的光在z方向传播5um后，沿传播方向看去的振幅分布
    z = np.arange(0,5,0.01)
    k = 2*np.pi/1.064

    # 将J改写成x和y方向的振幅序列
    J = np.abs(J)*np.cos(k*z+np.angle(J))
    
    fig = plt.figure()
    ax1 = fig.add_subplot(121)
    ax1.plot(J[0],J[1])
    ax2 = fig.add_subplot(122,projection='3d')
    ax2.plot3D(z,J[0],J[1])
    
    plt.show()

if __name__ == "__main__":
    drawJones()

其图像为

现在光路中引入$\lambda /4$波片，令其旋转一周，观察一下线偏振光会有怎样的变化，设波片与x轴所成夹角为$\Psi$，则其对应的Jones矩阵为

$$\left[\begin{array}{cc}1-\cos 2\Psi & -\text{i}\sin 2\Psi \\ -\text{i}\sin 2\Psi & 1+\cos 2\Psi \end{array}\right]$$

代码为

def drawQuaterPlate():
    quater = lambda psi : np.array(
        [[1-np.cos(2*psi),-1j*np.sin(2*psi)],
         [-1j*np.sin(2*psi),1+np.cos(2*psi)]])
    
    z = np.arange(0,5,0.01)
    k = 2*np.pi/1.064
    # 初始光波为x方向线偏振光
    J0 = np.array([1,0]).reshape(2,1)

    fig = plt.figure()
    ax = fig.add_subplot(111,xlim=(-2,2),ylim=(-2,2))
    ax.grid()

    line, = ax.plot([],[],lw=0.2)
    time_text = ax.text(0.1,0.9,'',transform=ax.transAxes)


    thetas = np.linspace(0,np.pi*2,80)

    def init():
        line.set_data([],[])
        time_text.set_text("")
        return line, time_text
    
    def animate(theta):
        J = quater(theta)@J0    #经过波片后的Jones矩阵
        arrJ = np.abs(J)*np.cos(k*z+np.angle(J)) #采样后的振幅
        line.set_data(arrJ[0],arrJ[1])
        time_text.set_text("angle:"+str(theta))
        return line, time_text

    ani = animation.FuncAnimation(fig, animate, thetas,
        interval=100, init_func=init)

    ani.save('polor.gif',writer='imagemagick')
    plt.show()

最终得到偏振情况与波片角度的关系

可见，当角度合适的时候，$\lambda /4$波片会将线偏振光变为圆偏振光。