电磁飞环的缓扩散传播特性*
2020-08-29王任刘胜张默然王秉中
王任 刘胜 张默然 王秉中
(电子科技大学应用物理研究所, 成都 610054)
1 引 言
无源麦克斯韦方程组是描述自由空间电磁场的基本方程, 平面波、X 波、贝塞尔波束、高斯脉冲和艾里脉冲等是麦克斯韦方程组的典型解[1−5]. 这些解对应的电磁场是典型的横波, 即传播方向的电场和磁场分量均为0. 而且, 这些解具有空时可分离特性, 即可以写成独立空间函数和独立时间函数的乘积[6,7]. 然而, Brittingham[8]和Ziolkowski 等[9,10]发现, 麦克斯韦方程组有一类有限能量解和上述典型解具有明显区别, 其空间函数和时间函数因具有内在联系而无法分离, 被称为电磁定向能量脉冲串(electromagnetic directed energy pulse trains,EDEPTs). 随后, 学者们发现了一系列EDEPTs,包括调制功率谱脉冲(modified power spectrum pulses)[9]、方位角依赖脉冲 (pulses with azimuthal dependence)[11]以及单周期电磁脉冲(single-cycle electromagnetic pulses)[12,13]等. 这些有限能量脉冲形态各异, 但都是麦克斯韦方程组的严格解.
在单周期电磁脉冲中, 有一种脉冲因其奇特的形态引起了学者们的关注, 即电磁飞环(flying electromagnetic toroid, FET), 也称飞翔的电磁甜 甜 圈 (flying electromagnetic doughnut)或 聚焦的电磁甜甜圈(focused electromagnetic doughnut)[12−15]. 近期研究表明, 电磁飞环具有空时不可分离特性、宽频谱特性、有限能量非衍射特性及能量回流特性[14−17]. 尤其引人注目的是, 电磁飞环具有超环面拓扑结构, 该拓扑结构具有与电磁飞环传播方向平行的场分量, 即纵向场分量[14−16]. 因为超环面拓扑结构和纵向场分量的存在, 电磁飞环在粒子加速和特殊场激励方面具有重要潜在应用[12−15].2017 年, Raybould 等[18]发现电磁飞环可以在介质球内同时激励起电偶极子模式和超环偶极子(toroidal dipole)模 式, 并 由 此 激 励 起 著 名 的Anapoles 非辐射模式, 为研究场与物质的相互作用提供了有效工具[19,20]. 作为应用和产生电磁飞环的基础, 对其特性的研究尤为重要. 如前文所述,学者们已经研究了电磁飞环的基本传播特性, 发现了电磁飞环具有一个“聚焦范围”, 在此范围内具有很多新奇特性[12−17].
本文通过理论计算对电磁飞环的传播特性进行了深入研究. 通过讨论电磁飞环横向和纵向场分量的空间分布和频谱特征, 发现电磁飞环在“聚焦范围”内保持超环面拓扑结构传播的同时, 也会缓慢向外发散, 即电磁飞环具有缓扩散传播特性. 本文的研究为电磁飞环的应用提供了理论基础.
2 电磁飞环的缓扩散传播特性
电磁飞环具有横磁(TM)和横电(TE)两种形式. TM 电磁飞环的磁场在与传播方向垂直的面上呈环形, 其电场绕磁场旋转, 整体呈现一种超环面拓扑结构, 如图1 所示. 在柱坐标系下, 沿z方向传播的TM 电磁飞环的电场和磁场可以表示为[12]:
式中,τ=z −ct,σ=z+ct,q1表示电磁飞环脉冲的有效波长,q2表示电磁飞环的聚焦范围,q1⩽q2. 有效波长指的是电磁飞环在传播方向上所覆盖的空间范围, 聚焦范围指的是电磁飞环传播过程中发散速度较慢的空间范围. 当|z|>q2时, 电磁飞环以与高斯脉冲相同的方式传播, 其有效波长为q1, 瑞利(Rayleigh)长度为q2. TM 电磁飞环的电场既有横向场分量Eρ, 又有纵向场分量Ez, 而其磁场只有横向场分量Hφ. TE 电磁飞环与TM 电磁飞环拓扑结构相同, 电场和磁场位置互换. 从(1)式—(3)式可知, 根据方程的实部和虚部可以分别构建两种不同的脉冲. 这两种脉冲都是麦克斯韦方程组的精确解. 实部对应的脉冲具有周期电场和单周期磁场, 虚部对应的脉冲具有单周期电场和周期磁场[14]. 因此, 我们将(1)式—(3)式实部和虚部对应的脉冲分别称为周期脉冲和单周期脉冲. 由于古伊相移作用, 单周期脉冲和周期脉冲会在传播过程中逐渐相互转化[12,13]. 不失一般性地, 在本文中我们以(1)式—(3)式的实部为例,对TM 电磁飞环的纵向场分量和横向场分量进行理论研究. 通过对方程中的时间变量t赋不同的值, 可以计算出电磁飞环在传播到不同位置时的场分布、拓扑结构及频谱分布, 从而在理论上准确研究电磁飞环的传播特性.
图1 TM 电磁飞环示意图Fig. 1. Schematic of TM FET.
考虑参数为q2= 100q1的TM 电磁飞环, 根据(1)式—(3)式画出其在t= 0 时刻的场分布, 如图2 所示.电磁飞环的各场分量在垂直于传播方向(z方向)的面上(xy面)都是旋转对称的. 纵向电场Ez在z= 0 附近沿+z方向, 在r= 0 处(即x= 0,y= 0)场值最大. 随着r增大, 场值逐渐减小, 而后电场方向反转, 变为沿–z方向. 在z= 0 平面,横向电场Er与横向磁场Hq场值均为0. 无论z取何值,Er与Hq在r= 0 处的场值均为零, 随着r的增大, 场值逐渐增加到最大值, 再逐渐减小为零.
在t= 0 时刻, 上述电磁飞环的纵向电场Ez的最大值位于直线r= 0 上, 横向电场Er的最大值位于柱面r= 4.5q1上, 这两个位置上的场分量如图3 所示. 由图可知, 在直线r= 0 上,Ez的最大值位于z= 0 处, 此处Er恒为零. 在柱面r=4.5q1上, 横向电场Er的最大值和最小值分居z= 0两侧, 即横向电场Er在z= 0 处改变方向; 纵向电场Ez在柱面r= 4.5q1上的最大值仍位于z= 0处, 即Er的0 值位置.
图2 电磁飞环在t = 0 时刻的场分布 (a) y = 0 平面Ez; (b) y = 0 平面Er; (c) y = 0 平面Hq; (d) z = 0 平面Ez; (e) z = –q1 平面Er; (f) z = –q1 平面HqFig. 2. Field distribution of the FET when t = 0: (a) Ez on the y = 0 plane; (b) Er on the y = 0 plane; (c) Hq on the y = 0 plane;(d) Ez on the z = 0 plane; (e) Er on the z = –q1 plane; (f) Hq on the z = –q1 plane.
图3 在t = 0 时刻的z 方向上的场分布 (a) 直线r = 0 处; (b) 柱面r = 4.5q1 上直线(x = 4.5q1, y = 0)处Fig. 3. Field distribution along z direction when t = 0: (a) on the line r = 0; (b) on the line (x = 4.5q1, y = 0) of cylindrical surface r = 4.5q1.
图2和图3 讨论了t= 0 时刻电磁飞环的拓扑特征. 值得说明的是, 电磁飞环的特性是随着传播过程不断变化的. 当上述TM 电磁飞环沿+z方向传播时, 其中心(即纵向场Ez最大值所在的位置)以光速向+z方向移动. 为了研究电磁飞环传播过程中结构的发散情况, 我们根据(1)式—(3)式理论计算了电磁飞环传播到各位置处的场, 并提取了每个传播状态下电磁飞环中心平面上的场值. 随着电磁飞环的传播, 即(1)式—(3)式中t 的增加,纵向场Ez和横向电场Er的变化如图4 所示. 图4的纵坐标为传播过程中电磁飞环中心在z 轴上的位置, 色值表示电磁飞环中心移动时纵向场Ez和横向电场Er最大值所在的水平面(xy 面)上的场分布. 从图中可以看出, 在电磁飞环传播过程中,纵向场Ez的最大值始终位于直线r = 0 上, 且随着传播过程逐渐减小. 当电磁飞环的中心位于z = 0处时, 横向场分量Er的最大值在r = 4.5q1处, 随着电磁飞环的传播, Er的最大值逐渐向外移动. 当电磁飞环传播到z = q2= 100q1时, 横向场分量Er的最大值移动至大约r = 10q1处. 这表明, 即使在一般认为的聚焦区域q2(或称非衍射区域)内[12−15], 电磁飞环也是不断向外扩散的, 只是扩散的范围较小, 即电磁飞环具有缓扩散传播特性. 以电磁飞环中心位于z = 0 面上时其 Eρ分量最大值的 ρ 坐标为基准(扩散0%), 电磁飞环传输距离(z 坐标, 以q1为单位)与环半径( ρ 坐标, 以q1为单位)的扩散关系如表1 所示. 值得说明的是, 电磁飞环在z < q2范围内虽然会缓慢扩散, 但其各场分量的相对位置并没有发生改变, 也就是说其仍然能够维持超环面拓扑结构.
图4 电磁飞环在传播过程中场分布的演化 (a) 纵向电场Ez; (b)横向电场Er, 黑线表示电磁飞环传播到不同位置时横向电场最大值所在位置Fig. 4. Evolution of the field distribution of FET: (a) Longitudinal electric field Ez; (b) transverse electric field Er, the black line indicates the position of maximum transverse electric field when the FET propagates to different positions.
为了更清晰地展示电磁飞环的缓扩散传输特性, 我们分别使用线性方程和对数方程对其传播到不同位置时横向电场最大值所在的位置轨迹进行了拟合,z
表1 电磁飞环传输距离z 与环半径r 的扩散关系Table 1. Relation between propagation distance and toroid radius of FET.
图5 电磁飞环传播到不同位置时横向电场最大值所在位置 (a) z < q2; (b) z > q2Fig. 5. Position of maximum transverse electric field when the FET propagates to different positions: (a) z < q2; (b) z > q2.
前文已经从时域讨论了电磁飞环的拓扑结构和传播特性, 接下来展示电磁飞环的频域特性. 根据(1)式—(3)式, 在空间中各点提取电磁飞环传播过程中经过该点的场值, 即获得该点对应的电磁飞环的时域信号, 然后对该信号做傅里叶变换, 即可得到电磁飞环在该位置处的频谱. 在z= 0 平面上记录电磁飞环各分量对应的时域场并做傅里叶变换, 得到的频谱分布如图6 所示. 由图可知, 纵向电场Ez的频谱在所有频点的最大值均位于r=0 处. 随着r的增大, 纵向场分量Ez各频点的谱值先减小后增加, 在频谱图上形成一条谱值为零的带状区域. 横向电场Er的频谱在r= 0 处为零, 随着r的增大, 各频点的值先增大后减小. 比较电磁飞环的横向电场和纵向电场的空间分布可以知道,各频点的电场均在r= 0 附近沿z方向穿过电磁飞环的中心, 然后沿r方向逐渐变为横向场, 再经由电磁飞环外边缘(即r较大的位置)演变为沿z方向的纵向场, 由此形成电场闭环, 即形成图1所示的超环结构.
图6 电磁 飞 环 经 过z = 0 平 面 时的 归 一 化 频 谱分 布(a)纵向电场Ez; (b)横向电场ErFig. 6. Normalized spectrum distribution when the FET propagates through the z = 0 plane: (a) Longitudinal electric field Ez; (b) transverse electric field Er.
从图6(a)可以, 纵向电场Ez的频谱在各个频点的最大值的位置都为r= 0, 那么横向电场Er频谱在各个频点的最大值的位置是否保持不变呢? 答案是否定的. 在上述q2= 100q1的TM 电磁飞环传播过程中, 计算垂直于传播方向的每个面的频谱, 并记录各个频点最大值的坐标, 如图7 所示.图7 中最外侧曲线对应频率为c/15q1的最大值的坐标, 越往内侧的曲线频率越高, 最内侧曲线对应频率为c/q1的最大值的坐标. 每一个频点的最大值坐标都随着z增大(即随着电磁飞环的传播)逐渐向外移动, 并且频率越低的曲线向外偏转的幅度越大. 这表明即使在z 图7 横向电场Er 各个频点最大值的位置曲线Fig. 7. Position curves of maximum value of each frequency of the transverse electric field Er. 本文通过理论计算揭示了电磁飞环的缓扩散传播特性. 电磁飞环的超环面拓扑结构在传输过程中的稳定性可以为粒子加速、特殊多极子激励等应用提供有力支持[12,18], 而本文的研究为这些应用提供了重要参考. 例如, 在粒子加速应用中, 要精确计算场与粒子的作用力, 电磁飞环的缓扩散传输特性是必须加以考虑的; 在激励Anapole 等模式时,因为缓扩散传输特性的存在, 电磁飞环起始位置与目标激励区域的距离将影响场激励的效果. 当前,关于电磁飞环特性的理论研究日趋完善, 要在实际系统中应用电磁飞环的诸多新奇特性, 则必须将电磁飞环实际产生出来. 因此, 如何产生电磁飞环将成为后续研究需要关注的重点. 感谢英国南安普顿大学Nikolay I. Zheludev 教授的指导.
3 结 论
