Scattering characteristics of Hermite-Gaussian beam on anisotropic cylinder

近些年来，高斯波束在各向异性材料中传播特性一直受到持续性的关注，该研究不仅具有重要的理论价值，而且一些成果已广泛应用于雷达散射截面、光信号处理、微波器件、天线罩、光纤的优化设计和微带天线等领域。1972年，ALEXOPOULOS等人^[1]研究了无限长非均匀介质圆柱对高斯波束片的散射。LOCK^[2]研究了无限长均匀圆柱体对入射聚焦高斯光束的散射问题。GUO等人^[3]提出一种有效而精确的递推算法来计算高斯波束垂直入射无限长多层圆柱的散射问题，并重点分析了圆柱的彩虹效应。随着广义洛伦兹-米理论的发展，HUANG和ZHANG等人^[4-6]通过将圆柱散射场、内场和入射的高斯光束用适当的圆柱矢量波函数展开，应用电磁场边界条件，解析地解决了高斯光束在单轴各向异性圆柱体上的散射特性。CHEN和ZHANG等人^[7]提出一种精确的半解析方法来计算回旋各向异性圆柱体对在轴高斯光束的散射问题。高斯波束的研究及其性质已经非常成熟^[8-14]，只有少量文献中研究了其它高斯波束的性质^[15-18]。

厄米-高斯波束是由高斯波束演化而来，具有重要的理论价值。作者在前人基础上，重点研究了厄米-高斯波束，首先讨论了电磁媒质的本构关系和分类，其次讨论了厄米-高斯波束在单轴各向异性圆柱中的散射特性，最后给出了数值结果和结论。

任意电磁媒质的本构关系可表达为：

式中，D是电位移，E是电场强度，B是磁感应强度，H是磁场强度，$\boldsymbol{\overline{\overline \xi }} $和$\boldsymbol{\overline{\overline \zeta }}$为手征参量张量，$\boldsymbol{\overline{\overline \varepsilon }}$和$ \boldsymbol{\overline{\overline \mu }}$分别为媒质的介电常数张量和磁导率张量。当媒质的磁导率为标量，介电常数为张量时，称媒质为电各向异性媒质，反之为磁各向异性媒质。在直角坐标系Oxyz中，电各向异性媒质中的介电常数可表示为$\overline{\overline \varepsilon } =\hat{x}\hat{x}{{\varepsilon }_{1}}+\hat{y}\hat{y}{{\varepsilon }_{2}}+\hat{z}\hat{z}{{\varepsilon }_{3}}$，而D与E的关系可表示为：

式中，下标x, y, z表示电场的3个方向分量; $ \hat{x}, \hat{y}, \hat{z}$为电场的3个方向的单位矢量。若ε₁，ε₂，ε₃全部相等，称为各向同性媒质; 若其中有两个相等，称为单轴各向异性媒质; 若3个对角元素均不相等，则称为双轴各向异性媒质。本文中主要研究单轴各向异性媒质ε₁=ε₂的情形。

如图 1所示^[4]，在直角坐标系Oxyz中有一个半径为r₀的无限长单轴各向异性圆柱体，高斯波束在自由空间且沿直角坐标系Ox′y′z′的z′轴正方向传播，高斯波束的束腰半径为w₀，束腰中心与圆心O′重合。波束的传播方向与z轴正方向夹角为β，在O′x′y′z′中点O坐标为(x₀, y₀, z₀)。随时间变化的部分规定为exp(-iωt)。

Figure 1.  The incidence of a Gaussian beam on an uniaxial anisotropic cylinder

正如参考文献[4]中所描述的，散射场用圆柱矢量波函数展开的表达式为：

同样参考文献[4]中的单轴各向异性圆柱内部电磁场表示为：

(3) 式~(6)式中，E_s, H_s表示散射场中电场和磁场，E_w, H_w表示圆柱内部的电场和磁场，λ=k₀sinζ，h=k₀cosζ，k₀为自由空间中的波数，ζ为圆柱矢量波与坐标轴z方向的夹角; E₀为电场振幅，η₀为波阻抗; m_mλq^(j), n_mλq^(j), l_mλq^(j)是圆柱矢量波函数, j=1, 2, 3分别对应三类贝塞尔函数; 而α_m(ζ)，β_m(ζ)，F_m1(ζ)和F_m2(ζ)是待求的未知系数。可以令a₁²=ω²ε₁μ₀，a₂²=ω²ε₃μ₀，ε₀和μ₀是自由空间中的介电常数和磁导率，则式中其它参量为：k₁=a₁，k₂=[a₁²a₂²+(a₁²-a₂²)k₀²cos²ζ]^1/2/a₁，λ₁=a₁²－k₀²cos²ζ，λ₂=a₂a₁^-1a₁²－k₀²cos²ζ，α_{1, e}(ζ)=1，β_{1, e}(ζ)=γ_1，e(ζ)=α_2，e(ζ)=0，$ {{\beta }_{2, \text{e}}}(\zeta )=-\text{i}\times \frac{a_{1}^{2}{{a}_{2}}}{\sqrt{(a_{1}^{2}-k_{0}^{2}{{\cos }^{2}}\zeta )[a_{1}^{2}a_{2}^{2}+(a_{1}^{2}-a_{2}^{2})k_{0}^{2}{{\cos }^{2}}\zeta ]}}$, $ {{\gamma }_{2, \text{e}}}(\zeta )=-\frac{a_{1}^{2}-a_{2}^{2}}{a_{1}^{2}}\ \frac{{{a}_{1}}{{a}_{2}}{{k}_{0}}\cos \zeta \sqrt{a_{1}^{2}-k_{0}^{2}{{\cos }^{2}}\zeta }}{a_{1}^{2}a_{2}^{2}+(a_{1}^{2}-a_{2}^{2})k_{0}^{2}{{\cos }^{2}}\zeta }$

电磁场的切向分量在r=r₀的边界连续，则边界条件可以表示为：

式中，E_i和H_i表示入射电磁波束的电场和磁场。将(3)式~(6)式代入到(7)式中，边界条件可以写为：

遵循投影法的一般理论步骤，分别在(8)式和(9)式两边点乘$ \hat{z}{{\text{e}}^{-\text{i}{{h}_{1}}z}}{{\text{e}}^{-\text{i}{m}'\varphi }}$和$ \hat{\varphi }{{\text{e}}^{-\text{i}{{h}_{1}}z}}{{\text{e}}^{-\text{i}{m}'\varphi }}$，然后在圆柱面进行积分，可以得到未知的展开系数与E_i和H_i的关系式：

式中，ξ=λr₀，ξ₁=λ₁r₀和ξ₂=λ₂r₀。

DAVIS和BARTON给出了一种计算方法^[19-20]，高斯波束TEM₀₀^(y′)或TE模在坐标系O′x′y′z′中电磁分量中各阶近似描述可表示为：

式中，s=1/(kw₀)，ψ₀(ξ, η, ζ)=iQexp(－iρ²Q)，Q=1/(i－2ζ)，ρ²=ξ²+η²。

根据电磁场理论中的对偶关系，即在(14)式~(19)式中做如下替换:E→－H，H→E，ε₀→μ₀，μ₀→ε₀，得到的电磁场分量仍然近似满足麦克斯韦方程组，是高斯波束的TEM₀₀^(x′)或TM模。

正如参考文献[19]中提到的，厄米-高斯波束的各种模式可以通过求解TEM₀₀^(y′)或TEM₀₀^(x′)的偏导数得到，即：

式中，ξ和η为无量纲的参量：ξ=x′/w₀，η=y′/w₀，w₀为高斯波束束腰半径。

为了得到厄米-高斯波束的具体表达式，以TEM₀₁^(y′)为例，即在(20)式中令m=0, n=1可得:

将(14)式~(19)式代入(22)式，即可求得TEM₀₁^(y′)模：

(23) 式~(28)式即为厄米-高斯波束的具体表达式。类似可得出厄米-高斯波束的TEM₁₀^(x′)模。

图 2代表的是厄米-高斯波束TEM₀₁^(y′)的强度分布图, 图 3代表的是厄米-高斯波束TEM₁₀^(x′)的强度分布图。图中横纵坐标ξ，η均是无量纲的参量。

Figure 2.  The intensity distribution of the TEM₁₀^(y′) mode Hermite-Gaussian beam

Figure 3.  The intensity distribution of the TEM₁₀^(x′) mode Hermite-Gaussian beam

入射的波束，其E_i和H_i可由(23)式~(28)式得到，代入到(10)式~(13)式，可以得到未知展开系数α_m(ζ)，β_m(ζ), F_m1(ζ)和F_m2(ζ)构成的线性方程组，从而求出这些未知系数。求出这些系数后再代入到(3)式~(6)式，进而可以求出散射场和内场。

定义归一化内场和近场的强度分布如下：

式中，E_{w, r}, E_{w, φ}, E_{w, z}分别为圆柱内部场中电场的3个分量; E_{i, r}, E_{i, φ}, E_{i, z}分别为入射场中电场的3个分量; E_{s, r}, E_{s, φ}, E_{s, z}分别为散射场中电场的3个分量。

对于高斯波束入射的情形，使用的参量与模型是：单轴各向异性圆柱a₁=3k₀，a₂=2k₀, 高斯波束的束腰半径w₀为2倍入射高斯波束的波长，圆柱的半径为5倍入射高斯波束的波长，入射角β=π/4，z₀=0。图 4表示高斯波束TE模通过单轴各向异性圆柱的归一化内场和近场。其中对于x(wavelength)轴上的范围，5~15表示入射场，-5~5表示圆柱的内场，-15~-5表示透射场，图中色柱表示的物理量是无量纲的。

Figure 4.  The normalized internal-field and near-field of a TE mode Gaussian beam through an uniaxial anisotropic cylinder

通过与已有方法及结果^[4]比较可以发现，两者实现了很好的吻合，这在很大程度上验证了作者方法的正确性。

使用相同的参量和模型，可得厄米-高斯波束入射的情形，图 5和图 6分别表示厄米高斯波束TEM₁₀^(x′)模和TEM₀₁^(y′)通过单轴各向异性圆柱的归一化内场和近场。

Figure 5.  The normalized internal-field and near-field of a TEM₁₀^(x′) mode Hermite-Gaussian beam through an uniaxial anisotropic cylinder

Figure 6.  The normalized internal-field and near-field of a TEM₁₀^(y′) mode Hermite-Gaussian beam through an uniaxial anisotropic cylinder

从图 5可以看出，TEM₁₀^(x′)模式厄米-高斯波束入射单轴各向异性圆柱时的反射场强度很弱。圆柱类似于凸透镜，有一个会聚作用，所以波束在通过圆柱后的近场强度明显增强，由入射波和反射波叠加而成的驻波现象也同样在图中表现的非常明显。在图 6中，波束透过圆柱后内部场强度逐渐增强，同样也有驻波现象。比较两图形，发现在相同的情况下TEM₁₀^(x′)模式厄米-高斯波束通过单轴各向异性圆柱时的近场强度比TEM₀₁^(y′)模式的强，而反射场强度弱。另一个值得注意的现象是，TEM₁₀^(x′)模式的厄米-高斯波束在通过单轴各向异性圆柱时有一个明显的折射现象，而在TEM₀₁^(y′)模式厄米-高斯波束入射时则表现的不明显。

主要基于单轴各向异性圆柱对厄米-高斯波束的散射特性进行研究。应用电磁场边界条件和投影法，精确半解析地得到了单轴各向异性圆柱对厄米-高斯波束散射特性和内场以及近场的归一化强度分布图，分析对比了两种不同的厄米高斯波束入射情形，发现波束透过圆柱后都有驻波现象。在相同的情况下TEM₁₀^(x′)模式厄米-高斯波束通过单轴各向异性圆柱时的近场强度比TEM₀₁^(y′)模式的强，而前者在通过单轴各向异性圆柱时折射现象更明显。