光学薄膜中的驻波场计算

最新推荐文章于 2024-10-04 12:40:24 发布

旋算仿真工作室

最新推荐文章于 2024-10-04 12:40:24 发布

阅读量418

点赞数 4

文章标签： matlab python 算法

本文链接：https://blog.csdn.net/2401_84906130/article/details/139907923

版权

激光在光路中传输时，反射镜会吸收激光能量，形成不均匀的温升分布，并产生热弹性应力\应变，甚至导致光学薄膜损伤，最终影响系统光学性能。在研究高功率激光辐照反射镜片表面的热形变分布时，通常会忽略其表面镀有的高反射膜系，我们采用菲涅尔系数矩阵法计算光学薄膜电场分布，有利于评估高功率光学系统性能。

反射镜中各层薄膜内的吸收、散射损耗与光强度在各层薄膜内的分布密切相关，当两列振幅相同的相干波各沿相反方向传播时，叠加后便产生驻波、反射镜膜层中的入射波和反射波正好能满足上述条件于是形成驻波场。当薄膜处在电磁波场中时，膜层之间可通过电矢量以矩阵方程式相联系，故采用矩阵法可求解出膜层内的电场分布。

图1 界面处的电场

图1表示两种折射率分别 $N_i$ 和 $N_j$ 所构成的界面，考虑到界面上两侧连续的边界条件，可得到如下方程：

$E_{i}^{-} =r_{ij}E_{i}^{+}+t_{ji}E_{j}^{-}\\ E_{j}^{+} =t_{ij}E_{i}^{+}+r_{ji}E_{j}^{-}$

(1)

式中， $r_{ij}$ 、 $r_{ji}$ 和 $t_{ij}$ 、 $t_{ji}$ 分别为界面上的菲涅尔反射系数和透射系数。由于 $r_{ij}=-r_{ji}$ ， $r_{ij}^2+t_{ij}t_{ji}=1$ ，通过代入上式可得如下矩阵：

$\left[\begin{array}{c}E_{j}^{+} \\ E_{j}^{-}\end{array}\right]=\frac{t_{i j}}{1-r_{i j}^{2}}\left[\begin{array}{cc}1& -r_{j j} \\ -r_{i j} &1\end{array}\right]\left[\begin{array}{l}E_{i}^{+} \\ E_{j}^{-}\end{array}\right]$

(2)

而膜层得逆矩阵可表示为

$M^{\prime}=\left[\begin{array}{cc}e^{-i \delta_{j}} &0\\0& e^{i \delta_{j}}\end{array}\right]$

(3)

将各膜层的入射波和反射波电场在不同界面两侧得关系用图2进行表示，可以得到膜层中任意界面上正向场和反向场，表示成矩阵形式的电矢量为

$\begin{array}{l}{\left[\begin{array}{l}E_{j}^{+} \\ E_{j}^{-}\end{array}\right]=\frac{t_{1} t_{2} \cdots t_{j}}{\left(1-r_{1}^{2}\right)\left(1-r_{2}^{2}\right) \cdots\left(1-r_{j}^{2}\right)}\left[\begin{array}{cc}e^{-i \delta_{j-1}} & -r_{j} e^{i \delta_{j-1}} \\ -r_{j} e^{-i \delta_{j-1}} & e^{i \delta_{j-1}}\end{array}\right]} \\ \cdots\left[\begin{array}{cc}e^{-i \delta_{1}} & -r_{2} e^{i \delta_{1}} \\ -r_{2} e^{-i \delta_{1}} & e^{i \delta_{1}}\end{array}\right]\left[\begin{array}{cc}1& -r_{1} \\ -r_{1} &1\end{array}\right]\left[\begin{array}{c}E_{0}^{+} \\ E_{0}^{-}\end{array}\right]\end{array}$