-
针对临近空间高超声速飞行器目标探测与识别研究的需求, 开展了高超声速飞行器非均匀等离子体电磁散射特性模拟测量研究. 利用弹道靶设备发射高超声速类HTV2模型形成模拟的超高速复杂外形目标, 弹道靶高精度阴影成像系统和雷达测量系统分别测量高超声速类HTV2模型姿态、全目标C波段/X波段电磁散射特性, 获得了不同实验条件下模型全目标雷达散射截面积(RCS)等实验数据. 研究结果表明: 在不同实验状态下, 包覆等离子体鞘套的高超声速类HTV2模型同一测量波段的RCS差别超过1个数量级, 模型姿态角对包覆等离子体鞘套的高超声速类HTV2模型RCS影响较大, 最大相差1个多数量级; 在给定的实验条件下, 模型尾迹C波段RCS远小于包覆等离子体鞘套的模型RCS, 模型尾迹X波段RCS显著增强; 高超声速类HTV2模型全目标C波段电磁散射能量主要分布在模型及其绕流区域, X波段电磁散射能量主要分布在模型及其绕流区域和等离子体尾迹区域. 根据弹道靶实验条件, 开展了包覆等离子体鞘套的高超声速类HTV2模型电磁散射特性数值仿真, 仿真结果与实验结果之间的最大误差小于4 dB, 验证了本文提出的非均匀等离子体包覆目标电磁散射特性建模方法的有效性.According to the requirements for target detection and recognition of hypervelocity vehicles in near space, the simulation and measurement of corresponding electromagnetic scattering characteristics of non-uniform plasma generated by hypervelocity targets are conducted. A numerical calculation method with dynamic plasma parameters is developed and hypervelocity HTV2-like models launched by the ballistic ranges are used to simulate complex shape target flying at a hypervelocity velocity. The high-precision shadow imaging systems and radar measurement systems of the ballistic range are used to measure the model postures, the electromagnetic scattering characteristics of the whole targets and their flow fields in both C band and X band. The experimental measurement results of the radar cross section (RCS) of the models and their flow fields under different experimental conditions are obtained. The results show that the numerical simulation methods of unsteady high-temperature ionized air flow can be used to simulate the unsteady thermal chemical flow fields around the head and body of the simplified hypervelocity HTV2-like flight models. The electromagnetic scattering characteristics of the models and their plasma sheath differ by more than one order of magnitude under different experimental conditions. The total RCS of the model’s wake in the C band is much smaller than that of the model, and the total RCS of the model’s wake in the X band is significantly enhanced. The attitude angles of the models have great influence on their electromagnetic scattering characteristics and their RCSs with different attitude angles can differ by one order of magnitude. The electromagnetic scattering energy of the model in the C band is distributed mainly around the targets and their flow fields surrounding them. The X band electromagnetic scattering energy of the model is distributed mainly in the regions around the targets, surrounding flow fields and the wake flow fields. According to the experimental conditions of the ballistic range, the numerical simulation analyses of the electromagnetic scattering characteristics of the models and flow fields around them are carried out, and the maximum error between simulation and experimental results is less than 4 dB, verifying the effectiveness of the modeling methods of simulating electromagnetic scattering characteristics of non-uniform plasma coated targets.
-
Keywords:
- hypervelocity /
- model /
- plasma sheath /
- ballistic range
1. 引 言
由于在空间光通信、量子纠缠、粒子操控等领域的应用, 光学涡旋已成为近年来的研究热点. 与普通光束相比, 涡旋光束可以携带轨道角动量. 拉盖尔-高斯(L-G)光束就是一种典型的涡旋光束, 其具有螺旋相位结构和环形场强分布, 光束中心的强度为零[1]. L-G光束沿其传播方向携带每光子
ℓℏ 的轨道角动量,ℓ 表示L-G光场的拓扑荷数[2].实验表明, 通过交换轨道角动量, L-G光束可以对物体施加扭矩, 包括微观吸附颗粒[3]、介观玻色-爱因斯坦凝聚[4]及亚微米布朗粒子[5]. 2007年, 基于此性质, 拉盖尔-高斯旋转腔系统被提出[6], 该系统由两个用作腔镜的螺旋相位元件构成, 分别为一个固定腔镜和一个围绕腔轴旋转的旋转腔镜(可视为扭摆). 在这种光旋转系统中, 腔模是L-G光, 它携带轨道角动量, 可以与螺旋相位元件交换轨道角动量. 近期人们对此类系统的兴趣逐渐增加, 并做了一些相关研究, 包括旋转腔镜的冷却[6,7]、光场轨道角动量的检测[8,9]、二阶边带效应[10]和纠缠[11,12]. 然而, 据我们所知, 关于L-G光束在拉盖尔-高斯旋转腔中的非互易传输现象至今还未被报道过. 由于光旋转相互作用[13]和光机械相互作用[14,15]之间存在类似的哈密顿量, 光机械相互作用产生的许多效应[16-29], 在拉盖尔-高斯旋转腔中通过与L-G光场轨道角动量的交换也会发生[6-12]. 这为利用光的轨道角动量而不是动量提供了可能性.
光的传播通常是互易的, 不过实现光学非互易对于构建非对称量子网络极为重要. 光学非互易性描述的是光场从一个方向通过光学系统但是不能沿原路返回的特性. 光学非互易是实现定向放大器、隔离器、循环器的基础[30], 并且需要打破时间反演对称性. 传统上, 非互易传输依赖于施加的偏置磁场 [31], 但这些传统器件通常体积庞大, 与超低损耗超导电路不兼容, 并且需要相当大的磁场. 近年来, 有人利用辐射-压力诱导光力耦合来打破时间反转对称性, 实现光的非互易效应[32-46]. 这一领域已经取得了重要成果, 包括非互易传输和放大[32-34]、非互易单光子效应[35]、非互易慢光[37]等. 文献[38]提出了一种基于具有光增益的光机械系统来实现光场和微波场之间光定向放大的方案, 发现放大方向可以通过调节有效光机耦合之间的相对相位(非互易相位差[46])来控制. 文献[43]给出了由两个光学膜和一个机械模组成的三模光机械系统中的光学非互易响应, 并证明了光学非互易响应是通过调节光机械耦合速率之间的相位差来打破系统的时间反转对称性来实现的. 然而, 通过光旋转耦合来实现携带轨道角动量的涡旋光束的非互易传输还有待进一步研究, 且之前工作在利用光力耦合来实现普通光束的非互易传输时, 都是利用光学模式与机械模式交换动量来实现, 而非交换角动量.
本文通过在单个拉盖尔-高斯旋转腔中额外加一个固定腔镜来构建双拉盖尔-高斯旋转腔, 并通过光纤将两个腔模线性耦合, 并且用两个强驱动场驱动不同的腔模, 进而研究涡旋光束的非互易传输现象. 将以往通过光与机械振子交换动量来实现普通光束的非互易传输推广到光与旋转腔镜交换角动量来实现涡旋光束的非互易传输. 研究发现, 非互易相位差可以决定涡旋光束非互易传输的产生及方向; 在一定的拓扑荷比值下, 该系统可以实现完美的非互易性; 系统耗散对传输振幅的变化有比较大的影响; 此外, 通过调节非互易相位差, 在该系统中可以实现非互易的慢光效应. 在此系统中的涡旋光束传输的非互易性, 来源于光旋转耦合相互作用和腔模耦合相互作用形成的两条路径之间发生的量子干涉效应.
2. 理论模型与主要公式
如图1所示, 本方案考虑了双拉盖尔-高斯旋转腔光旋转系统, 其中3个腔镜均为螺旋相位元件, 其中两个输入腔镜部分透明且刚性固定, 不改变透过它们光束的轨道角动量而会给反射光束增加
−2ℓ1 和−2ℓ2 的拓扑荷; 处于中间位置的旋转腔镜是完全反射的, 可以给反射光增加+2ℓ1 和+2ℓ2 的拓扑荷, 并可以在支架上绕腔轴z 旋转(角平衡位置ϕ0=0 ); 同时两个腔场通过光纤线性耦合在一起. 两个频率均为ωc , 振幅分别为εc 和εd 的强驱动场分别从两侧注入系统驱动频率为ω0 的L-G腔模. 利用一束频率为ωp 、振幅为εL (εR )的弱探测场从系统左侧(右侧)入射, 从而检查该光旋转系统相对于左侧(右侧)探测光的响应特性. 在驱动场频率ωc 的旋转框架下, 系统的哈密顿量(ℏ = 1)可写为图 1 双拉盖尔-高斯旋转腔示意图, 两L-G腔场通过光旋转相互作用与中间的旋转腔镜耦合, 振幅为 和εc (εd 和εL )的强驱动场(弱探测场)分别从两侧入射系统来驱动L-G腔模εR 和c1 , 同时用光纤将两L-G腔模线性耦合. 旋转腔镜的平衡位置为c2 , 在扭力作用下的角位移用ϕ0 角表示ϕ Fig. 1. Schematic diagram of double Laguerre-Gaussian (L-G) rotational-cavity. The two L-G cavity modes are coupled with a rotating cavity mirror in the middle via the optical rotation interaction. Two strong pump fields (weak probe fields) with amplitudes andεc (εd andεL ) are incident on the system from both sides to drive the L-G cavity modesεR andc1 , and the two L-G cavity modes are linearly coupled with an optical fiber. The equilibrium position of the rotational mirror isc2 , and the angular displacement is indicated by angleϕ0 under the action of the torsion.ϕ H=Δc(c†1c1+c†2c2)+L2z2I+12Iω2ϕϕ2 +gϕ1c†1c1ϕ+gϕ2c†2c2ϕ+J(c†1c2+c†2c1) +i(εcc†1−ε∗cc1)+i(εdc†2−ε∗dc2) +iεL(c†1e−iδt−c1eiδt)+iεR(c†2e−iδt−c2eiδt), (1) 式中, 第1项表示L-G腔模的自由哈密顿量, 其中
c1(c†1) 和c2(c†2) 为两个腔的湮灭(产生)算符,Δc=ω0−ωc 为驱动场与腔模的失谐; 第2项和第3项给出了旋转镜的自由哈密顿量,ωϕ ,Lz 和ϕ 分别为旋转镜的角频率、角动量和角位移,I=MR2/2 为旋转镜的转动惯量, 其中M 和R 分别为旋转镜的质量和半径; 第4项和第5项表示L-G腔模与旋转镜的光旋转耦合, 单光子耦合强度分别为gϕ1=cℓ1/L 和gϕ2=cℓ2/L , 其中L为两个腔的长度; 第6项表示两个腔模的线性耦合, 其中J 为线性耦合强度; 最后4项分别表示强驱动场和弱探测场对腔模的驱动,δ=ωp−ωc 为探测场与驱动场之间的失谐.为了方便计算, 设定
b=√Iωϕ2(ϕ+iIωϕLz) 和b†=√Iωϕ2(ϕ−iIωϕLz) 分别表示转动腔镜的湮灭(产生)算符, 系统的哈密顿量可以写为H=Δc(c†1c1+c†2c2)+ωϕb†b+(gϕ1c†1c1+gϕ2c†2c2)(b†+b)+J(c†1c2+c†2c1)+i(εcc†1−εc∗c1)+i(εdc†2−ε∗dc2)+iεL(c†1e−iδt−c1eiδt)+iεR(c†2e−iδt−c2eiδt). (2) 由系统的哈密顿量, 可以写出描述系统的动力学的相关算符的海森伯-郎之万方程:
˙c1=−[iΔc+κ12+igϕ1(b†+b)]c1+εc+εLe−iδt−iJc2, ˙c2=−[iΔc+κ22+igϕ2(b†+b)]c2+εd+εRe−iδt−iJc1, ˙b=−iωmb−γϕ2b−i(gϕ1c†1c1+gϕ2c†2c2), (3) 其中
κj 和γϕ 分别为两腔和转动腔镜的弛豫速率. 在乘积算符的平均值满足分解假设⟨bcj⟩=⟨b⟩⟨cj⟩ 且忽略探测场的情况下, 可得到算符的稳态解:bs=−i(gϕ1|c1s|2+gϕ2|c2s|2)γϕ2+iωϕ, c1s=(κ22+iΔ2)εc−iJεdJ2+(κ12+iΔ1)(κ22+iΔ2), c2s=(κ12+iΔ1)εd−iJεcJ2+(κ12+iΔ1)(κ22+iΔ2), (4) 其中
Δ1,2=Δc+gϕ1,ϕ2(bs+b*s) 为L-G腔模和驱动场之间的有效失谐.在强驱动场的驱动下, 每个算符可以由其平均值和量子涨落的和表示, 即
cj=cjs+δcj ,b=bs+δb . 将cj 和b 的表达式代入方程(3), 可以得到线性化的海森伯-郎之万方程:δ˙c1=−(κ12+iΔ1)δc1−igϕ1c1s(δb†+δb)+εLe−iδt−iJδc2, δ˙c2=−(κ22+iΔ2)δc2−igϕ2c2s(δb†+δb)+εRe−iδt−iJδc1, δ˙b=−(γϕ2+iωϕ)δb−igϕ1(c1sδc†1+c∗1sδc1)−igϕ2(c2sδc†2+c∗2sδc2), (5) 采用变换
δb→δbe−iωϕt ,δc1→δc1e−iΔ1t ,δc2→δc2e−iΔ2t , 则方程(5)为δ˙c1=−κ12δc1−iG1[δb†ei(ωϕ+Δ1)t+δbe−i(ωϕ−Δ1)t]+εLe−i(δ−Δ1)t−iJδc2ei(Δ1−Δ2)t, δ˙c2=−κ22δc2−iG2eiθ[δb†ei(ωϕ+Δ2)t+δbe−i(ωϕ−Δ2)t]+εRe−i(δ−Δ2)t−iJδc1ei(Δ2−Δ1)t, δ˙b=−γϕ2δb−iG1[δc1ei(ωϕ−Δ1)t+δc†1ei(ωϕ+Δ1)t]−iG2[e−iθδc2ei(ωϕ−Δ2)t+eiθδc†2ei(ωϕ+Δ2)t], (6) 其中
G1 ,G2 为两个涡旋腔场与旋转腔镜之间的有效光旋转耦合强度,G1=gϕ1c1s ,G2=gϕ2c2se−iθ ,gϕ1c1s 和gϕ2c2s 之间的非互易相位差θ 可以通过驱动场εc 和εd 的相对相位来调节.如果每个强驱动场都处于力学红边带(
Δ1≈Δ2≈ωϕ ), 同时旋转镜频率ωϕ 远大于有效耦合强度Gj , 方程(6)可以化简为δ˙c1=−κ12δc1−iG1δb−iJδc2+εLe−ixt, δ˙c2=−κ22δc2−iG2eiθδb−iJδc1+εRe−ixt, δ˙b=−γϕ2δb−iG1δc1−iG2e−iθδc2, (7) 其中
x=δ−ωϕ . 由方程(7)的具体形式, 可以假设方程(7)的解具有下面的形式:δc1=δc1+e−ixt+δc1−eixt, δc2=δc2+e−ixt+δc2−eixt, δb=δb+e−ixt+δb−eixt. (8) δc1+=(8G22+2γxκ2x)εL−(8G1G2e−iθ+4iJγx)εRΩ−16iG1G2Jcosθ, δc2+=(8G21+2γxκ1x)εR−(8G1G2eiθ+4iJγx)εLΩ−16iG1G2Jcosθ, b+=(4iG1κ2x+8G2Je−iθ)εL+(8G1J+4iG2κ1xe−iθ)εR16iG1G2Jcosθ−Ω, (9) 其中
γx=γ−2ix ,κjx=κj−2ix ,Ω=4G22κ1x+4G21κ2x+4J2γx+κ1xκ2xγx .为了研究涡旋光束传输的非互易性, 需要得到
εoutL 和εoutR . 根据输入输出关系[47,48]得εoutL+εinLe−ixt=√κ1δc1, εoutR+εinRe−ixt=√κ2δc2, (10) 其中
εinL, R=εL, R/√κ1κ2 , 输出场具有δc1=δc1+e−ixt+δc1−eixt 和δc2=δc2+e−ixt+δc2−eixt 的形式, 由方程(9)和方程(10)得εoutL + =√κ1δc1+−εL/√κ1, εoutR + =√κ2δc2+−εR/√κ2, εoutL−=0, εoutR −=0. (11) 本文定义
T12 (T21 )为入射腔c1(c2) 的探测场通过系统后输出时的传输振幅, 根据(9)式和(11)式计算得T12=TL→R=|εoutR+εinL|εinR=0=−√κ1κ2(8G1G2eiθ+4iJγx)Ω−16iG1G2Jcosθ, T21=TR→L=|εoutL+εinR|εinL=0=−√κ1κ2(8G1G2e−iθ+4iJγx)Ω−16iG1G2Jcosθ. (12) 3. 涡旋光束传输非互易性
本节将详细说明在双拉盖尔-高斯旋转腔系统中非互易相位差
θ 、涡旋光束所携带拓扑荷ℓ 和系统耗散对传输振幅的影响. 首先讨论非互易相位差θ 对传输振幅的影响.图2给出传输振幅
T12 (红色实线)和T21 (蓝色虚线)在不同非互易相位差θ 下随标准化失谐x/κ1 的演化. 可以看出, 相位差θ 对涡旋光束传输有比较大的影响, 在θ = 0时, 探测场在两个方向上的传输是互易的, 但当θ≠0 探测场的传输会呈现明显的非互易性. 观察图2(b)—图2(e)可以发现, 非互易现象会出现在共振点及附近, 还可以发现相位差的正负可决定非互易的方向(θ<0 时,T21>T12 ; 而θ>0 时,T21<T12 ). 当θ=±π /4 时, 非互易性在x=−0.4κ 处呈现的最明显. 当θ=±π /2 时, 系统可以在共振处(x=0 )呈现完美的涡旋光束非互易传输(T12=0 时T21=1 , 或T12=1 时T21=0 ). 由(12)式也可以看出, 当θ=nπ 时,T12=T21 , 涡旋光束传输呈现互易性, 而当θ≠nπ 时,T12≠T21 , 即涡旋光束的传输是非互易的. 实际上, 这里的涡旋光束的非互易传输主要源于光子不同传输路径之间的量子干涉作用. 为了更清楚地展示更多非互易相位差对传输振幅的影响, 图2(f)给出了传输振幅T12 (红色实线)和T21 (蓝色虚线)随非互易相位差θ 的演化.图 2 传输振幅 (红色实线)和T12 (蓝色虚线)在不同非互易相位差T21 下随标准化失谐θ 的演化 (a)x/κ1 ; (b)θ=0 ; (c)θ=−π /4 ; (d)θ=π /4 ; (e)θ=−π /2 . (f)在标准化失谐为零(x = 0)时, 传输振幅θ=π /2 (红色实线)和T12 (蓝色虚线)随非互易相位差T21 的演化. 其他参数为θ ,κ1=1 ,κ2=1 ,γ=1 ,Gi=√κiγ/2 .J=√κ1κ2/2 Fig. 2. Transmission amplitudes (red solid line) andT12 (blue dotted line) versus normalized detuningT21 under different nonreciprocal phase difference: (a)x/κ1 ; (b)θ=0 ; (c)θ=−π /4 ; (d)θ=π /4 ; (e)θ=−π /2 . (f) Transmission amplitudesθ=π /2 (red solid line) andT12 (blue dotted line) versus nonreciprocal phase difference withT21 . Other parameters arex=0 ,κ1=1 ,κ2=1 ,γ=1 ,Gi=√κiγ/2 .J=√κ1κ2/2 接下来讨论系统耗散对涡旋光束传输非互易性的影响. 为不失一般性, 取非互易相位差
θ=π /2 . 图3(a)和图3(b)给出了传输振幅T12 和T21 在旋转腔镜的不同耗散率γ 下随标准化失谐x/κ1 的演化. 可以看出在共振点(x = 0)附近, 随着旋转镜耗散率γ 增加, 从腔c1 到腔c2 的传输振幅T12 = 1 保持不变, 而从腔c2 到腔c1 的传输振幅T21 呈现出明显的变化, 显示出先降低后增加的趋势. 当旋转镜耗散率γ = κ1 = 1 时, 涡旋光束的传输非互易性达到最佳(T12=1 ,T21=0 , 在共振点处).图 3 (a)传输振幅 和(b)传输振幅T12 作为标准化失谐T21 和旋转镜与腔x/κ1 耗散比c1 的函数; (c)传输振幅γ/κ1 和(d)传输振幅T12 作为标准化失谐T21 和腔x/κ1 与腔c2 耗散比c1 的函数. 其他参数为κ2/κ1 ,θ=π /2 ,κ1=1 G1=G2=J=1/2 Fig. 3. (a) Transmission amplitude and (b) transmission amplitudeT12 as a function of the normalized detuningT21 and the dissipation ratiox/κ1 of rotating mirror to cavityγ/κ1 ; (c) transmission amplitudesc1 and (d) transmission amplitudesT12 as a function of the normalized detuningT21 and the dissipation ratiox/κ1 of cavityκ2/κ1 to cavityc2 . Other parameters arec1 ,θ=π /2 ,κ1=1 .G1=G2=J=1/2 图3(c)和图3(d)给出了传输振幅
T12 和T21 在腔c2 的不同耗散率下随标准化失谐x/κ1 的演化. 可以看出, 在共振点(x = 0)附近, 随着腔c2 耗散率κ2 增加, 从腔c2 到腔c1 的传输振幅T21 = 0 保持不变, 而从腔c1 到腔c2 的传输振幅T12 呈现逐渐减小的趋势. 当腔c2 耗散率κ2 = κ1 = 1 时, 涡旋光束传输的非互易性达到最佳(T12=1 ,T21=0 , 在共振点处). 所以对于涡旋光束的非互易传输来说, 腔场c2 的耗散不是越低越好.接下来讨论两个涡旋光场所携带拓扑荷比值
ℓ2/ℓ1 对传输振幅的影响. 通常, 涡旋光场所携带拓扑荷实际上会影响旋转腔镜与L-G腔模的耦合强度. 为了方便, 这里通过调节驱动场振幅εc 和εd 的大小从而使得拓扑荷比ℓ2/ℓ1 与有效耦合强度比G2/G1 相等. 图4(a)—图4(f)给出了在标准化失谐x = 0时传输振幅T12 和T21 在不同的非互易相位差下随两个涡旋光场所携带拓扑荷比ℓ2/ℓ1 的演化. 可以发现, 当θ=0 和θ=π 时, 拓扑荷的比值对传输振幅T12 和T21 的影响是相同的, 即系统对涡旋光束的传输呈现互易性. 当θ=±π /4 和θ=±π /2 时, 系统对涡旋光束传输呈现明显的非互易性. 当θ=±π /4 时, 随着拓扑荷比值的增加, 传输振幅T12 和T21 的变化趋势是相同的, 且是偶对称, 但相同拓扑荷比值所对应的传输振幅T12 和T21 是不同的, 在拓扑荷比值为±1 处的非互易性是最好的. 当θ=±π /2 时, 随着拓扑荷比值的增加, 传输振幅T12 仍呈现关于拓扑荷比值等于0 处的对称性, 但此时传输振幅T21 保持为0 不变. 在拓扑荷比值为±1 处, 系统对涡旋光束传输可以呈现完美的非互易性(T12=0 时T21=1 , 或T12=1 时T21=0 ). 此外, 拓扑荷的正负代表涡旋光矢量的旋转方向, 根据涡旋光矢量的旋转方向可以将涡旋光束分为左旋涡旋光束和右旋涡旋光束. 从图4(d)和图4(e)可以了解到涡旋光束想要呈现完美的非互易性与涡旋光束的左、右旋无关, 只要保证|ℓ2/ℓ1| = 1 .图 4 传输振幅 (红色圆圈)和T12 (蓝色圆圈)在不同非互易相位差T21 下随两个涡旋光场所携带拓扑荷之比θ 的演化 (a)ℓ2/ℓ1 ; (b)θ=0 ; (c)θ=−π /4 ; (d)θ=π /4 ; (e)θ=−π /2 ; (f)θ=π /2 . 其他参数为θ=π ,κ1=1 ,κ2=1 ,γ=1 J=2G1G2/γ Fig. 4. Transmission amplitudes (red circle) andT12 (blue circle) with ratio of topological charges carried by two vortex optical fieldsT21 under different nonreciprocal phase difference: (a)ℓ2/ℓ1 ; (b)θ=0 ; (c)θ=−π /4 ; (d)θ=π /4 ; (e)θ=−π /2 ; (f)θ=π /2 . Other parameters areθ=π ,κ1=1 ,κ2=1 ,γ=1 .J=2G1G2/γ 4. 非互易慢光效应
在腔系统中, 光的群延迟是描述光群速度的一个重要指标, 首先介绍群延迟的定义, 即[49,50]
τij=dΘ21(12)dωp, (13) 其中
Θ21(12) 是在频率为ωp 下输出场c1→c2 (c2→c1 )的相位, 群延迟τij<0 对应快光,τij>0 对应慢光.图5(a)—图5(c)给出了针对涡旋光束所对应的群延迟
τ12 和τ21 在不同的非互易相位差下随两个涡旋光场所携带拓扑荷比ℓ2/ℓ1 的演化. 在θ=0 时群速度是互易的且关于拓扑荷比值偶对称. 在θ=±π /2 时, 系统可以呈现明显的非互易慢光效应, 且相位的变化可以影响非互易慢光的群延迟的变换. 随着拓扑荷比值的增加,τ12 和τ21 均呈现先增大后减小的趋势且关于拓扑荷比值偶对称, 但群延迟的差值(τ12−τ21 或τ21−τ12 )一直保持不变.图 5 群延迟 (红色圆圈)和τ12 (蓝色圆圈)在不同非互易相位差τ21 下随两个涡旋光场所携带拓扑荷比θ 的演化 (a)ℓ2/ℓ1 ; (b)θ=0 ; (c)θ=−π /2 . 其他参数为θ=π /2 ,κ1=1 ,κ2=1 ,γ=1 J=2G1G2/γ Fig. 5. Group delay (red circle) andτ12 (blue circle) with ratio of topological charges carried by two vortex optical fieldsτ21 under different nonreciprocal phase difference: (a)ℓ2/ℓ1 ; (b)θ=0 ; (c)θ=−π /2 . Other parameters areθ=π /2 ,κ1=1 ,κ2=1 ,γ=1 .J=2G1G2/γ 5. 实验可行性分析
最后, 根据目前的实验进展, 对本文使用的系统方案的可行性进行了讨论. 本文讨论的双拉盖尔-高斯旋转腔系统是在文献[6]提出的由两个螺旋相位元件所构成的拉盖尔-高斯旋转腔的基础上, 添加一个固定的螺旋相位元件所实现的. 通过三个螺旋相位元件作为腔镜来组成两个光腔, 并且用光纤将两个光腔线性耦合. 其中, 两个输入腔镜FM1和 FM2为部分透明且被刚性固定, 不会改变透射光束的轨道角动量但会给反射光束增加
−2ℓ1 和−2ℓ2 的拓扑荷, 这一效应已经被实验观测到[51]. 另外, 处于中间位置的旋转腔镜RM是完全反射的, 可以给反射光增加+2ℓ1 和+2ℓ2 的拓扑荷[6]. 随着纳米技术的发展, 螺旋相位元件可以利用超精密点金刚石车削机床直接加工铝盘表面得到[6,52]. 利用螺旋相位元件的反射和透射可以改变激光光束的方位结构, 进而改变其携带的拓扑荷的值. 拓扑荷ℓ 依赖于螺旋相位元件的结构, 通过将螺旋相位元件的方位角坡道划分为离散的阶跃并控制阶跃高度和旋向, 可以设计出具有特定拓扑荷ℓ 的涡旋光场[6,53]. 最近实验已经可以通过螺旋相位元件产生携带拓扑荷高达1000的L-G光束[54]. 对于本工作, 若想实现理想的涡旋光隔离器, 就需要在本系统中实现涡旋光束的完美非互易性, 即涡旋光束的单向传输. 在理论分析中得到, 若要呈现完美的非互易涡旋光束传输, 需要保证|ℓ2/ℓ1| = 1 , 因此这里的涡旋光束所携带的拓扑荷需是可以调节的. 在实验上, 可以通过调整旋转腔镜两侧的阶跃高度和旋向来实现. 另外, 理论分析中得到非互易相位差对涡旋光束的非互易传输有着重要影响, 而非互易相位差的改变在实验上可以通过调节两个强驱动激光场εc 和εd 的相对相位来实现[45]. 值得强调的是, 本文所构建的系统在实验实现时需要较高精细度的光腔, 在光腔的精细度比较低时, 在该系统所要探究相关现象会不明显[53].6. 结 论
本文研究了涡旋光束在双拉盖尔-高斯旋转腔系统中的非互易传输. 将以往大多在光机械系统中研究的普通光束的非互易性推广到在光旋转系统中研究携带轨道角动量的涡旋光束传输的非互易性. 系统中, 光旋转耦合作用以及光纤的线性耦合作用提供的两条路径, 从而可以产生量子干涉. 首先, 研究了非互易相位差对传输振幅的影响, 发现相位差可以决定涡旋光束非互易性的发生及方向. 然后, 在相位差
θ=π /2 时, 研究了系统耗散对传输振幅的影响, 分析了涡旋光束实现完美的非互易传输时系统的耗散应该满足的条件. 接下来, 分析了涡旋光束携带的拓扑荷比值对系统非互易性的影响, 探究了实现涡旋光束完美的非互易传输时拓扑荷比值应该满足的条件, 发现涡旋光的左、右旋不会对传输产生影响. 拓扑荷对光传输产生影响, 主要是因为涡旋光束携带的拓扑荷数会影响旋转腔镜与腔模的耦合强度. 最后, 分析了该系统产生的非互易的慢光效应. 这些研究成果可用于设计针对携带轨道角动量的涡旋光束的理想光隔离器, 有望应用于光通信等领域.[1] Akey N D, Cross A E 1970 NASA-TN-D-5615 [1970-02-01]
[2] Dirsa E F 1960 Pro. IRE 48 703
Google Scholar
[3] Golden K E, Pridmore D C, Stewart G E 1970 NASA 19710011626 (Washington: NASA LANGLEY Research Center)
[4] Wood G E, Asmar S W, Rebold T A 1997 TDA Progress Report 42 131
[5] Bachynski M P, Gibbs B W 1970 NASA 19710011649 (Washington: NASA LANGLEY Research Center)
[6] Sotnikov V L, Leboeuf J N, Mudaliar S 2010 IEEE Trans. Plasma Sci. 38 2208
Google Scholar
[7] Usui H, Yamashita F, Matsumoto H 1999 Adv. Space Res. 24 1069
Google Scholar
[8] Mather D E, Pasqual J M, Sillence J P, Lewis P 2005 AIAA/CIRA 13th International Space Planes and Hypersonics Systems and Technology Capua, Italy, May 16–20, 2005 p3443
[9] 周超, 张小宽, 张晨新, 吴国成 2014 现代雷达 36 83
Google Scholar
Zhou C, Zhang X K, Zhang C X, Wu G C 2014 Modern Radar 36 83
Google Scholar
[10] 吴巍, 刘方, 钟建林, 王国宏 2019 电波科学学报 34 610
Wu W, Liu F, Zhong J L, Wang G H 2019 Chin. J. Radio. Sci. 34 610
[11] 金铭, 韦笑, 吴洋, 张羽淮, 余西龙 2015 物理学报 64 205205
Google Scholar
Jin M, Wei X, Wu Y, Zhang Y H, Yu X L 2015 Acta Phys. Sin. 64 205205
Google Scholar
[12] 马平, 石安华, 杨益兼, 于哲峰, 黄洁 2015 强激光与粒子束 27 073201
Google Scholar
Ma P, Shi A H, Yang Y J, Yu Z F, Huang J 2015 High Power Laser Part. Beams 27 073201
Google Scholar
[13] 马平, 石安华, 杨益兼, 于哲峰, 梁世昌, 黄洁 2017 物理学报 66 102401
Google Scholar
Ma P, Shi A H, Yang Y J, Yu Z F, Liang S C, Huang J 2017 Acta Phys. Sin. 66 102401
Google Scholar
[14] 邾继贵, 于之靖 2012 视觉测量原理与方法 (北京: 机械工业出版社) 第4678页
Zhu J G, Yu Z J 2012 Principles and Methods of Visual Measurements (Beijing: China Machine Press) pp46–78 (in Chinese)
[15] Sun H Y, Cui Z W, Wang J J, Han Y P, Shi P 2018 Phys. Plasmas 25 063514
Google Scholar
[16] Sun H Y, Wang J J, Han Y P, Cui Z W, Sun P, Shi X W, Zhao W J 2018 Int. J. Antennas Propag. 1 14
[17] 孙浩宇 2018 博士学位论文 (西安: 西安电子科技大学)
Sun H Y 2018 Ph. D. Dissertation (Xi’an: Xidian University) (in Chinese)
[18] 艾夏 2013 博士学位论文 (西安: 西安电子科技大学)
Ai X 2013 Ph. D. Dissertation (Xi’an: Xidian University) (in Chinese)
[19] 陈安涛 2019 博士学位论文 (西安: 西安电子科技大学)
Chen A T 2019 Ph. D. Dissertation (Xi’an: Xidian University) (in Chinese)
[20] 葛德彪, 闫玉波 2005 电磁波时域有限差分方法 (西安: 西安电子科技大学出版社) 第88—89页
Ge D B, Yan Y B 2005 Finite-Difference Time-Domain Method for Electromagnetic Waves (Xi’an: Xidian University Press) pp88–89 (in Chinese)
[21] Ai X, Han Y, Li C Y, Shi X W 2011 Prog. Electromagn. Res. Lett. 22 83
Google Scholar
[22] Ai X, Han Y, Chen Z, Shi X W 2011 Prog. Electromagn. Res. M. 18 143
Google Scholar
[23] Ai X, Tian Y, Han Y P, Shi X W, Li W T 2013 J. Quant. Spectrosc. Radiat Transfer 124 28
Google Scholar
[24] 陈伟, 郭立新, 李江挺, 淡荔 2017 物理学报 66 084102
Google Scholar
Chen W, Guo L X, Li J T, Dan L 2017 Acta Phys. Sin. 66 084102
Google Scholar
-
表 1 弹道靶高超声速球模型全目标RCS实验测量结果与数值计算结果对比
Table 1. Comparison between measurement and numerical results of the RCS of the ball models flying at hypervelocity.
实验状态 球模型及其等离子体鞘套套RCS/dBsm 尾迹RCS/dBsm 压力/kPa 速度/(km·s–1) 测量实验值 计算结果 实验测量值 计算结果 4.2 5.0 –31.40 –31.31 –52.10 –50.63 表 2 高超声速类HTV2模型C波段全目标RCS实验测量结果
Table 2. Measurement results of the C band full target RCS of the simplified HTV2 models flying at hypervelocity.
实验状态 模型姿态 高超声速类HTV2模型全目标RCS实测结果/dBsm 速度/(km·s–1) 压力/kPa 俯仰角/(°) 偏航角/(°) 滚转角/(°) 模型全目标总RCS 本体及绕流场RCS 尾迹RCS 5.0 7.6 22.43 5.18 103.65 5.91 5.90 –18.83 5.0 8.8 –1.85 14.79 –107.27 –4.19 –4.27 –21.57 5.0 10.4 –7.15 –2.10 69.40 –2.61 –2.72 –18.75 5.0 15.3 6.02 –10.1 –109.07 1.90 1.87 –19.79 表 3 高超声速类HTV2模型X波段全目标RCS实验测量结果
Table 3. Measurement results of the X band full target RCS of the simplified HTV2 models flying at hypervelocity.
实验状态 模型姿态 高超声速类HTV2模型全目标RCS实测结果/dBsm 速度/(km·s–1) 压力/kPa 俯仰角/(°) 偏航角/(°) 滚转角/(°) 模型全目标总RCS 本体及绕流场RCS 尾迹RCS 4.0 8.5 –17.00 23.00 6.80 –8.75 –8.76 –35.56 5.0 6.6 7.50 –11.00 108.28 1.05 1.05 –30.70 5.0 8.2 18.93 17.11 109.95 –2.92 –3.44 –12.43 5.0 11.6 4.03 –3.25 107.73 –3.63 –3.94 –15.23 5.0 14.8 –7.00 9.10 70.93 –5.99 –9.36 –8.68 表 4 高超声速类HTV2模型不同波段RCS的FDTD方法数值模拟结果与弹道靶实验测量结果的对比
Table 4. Various bands RCS comparisons between FDTD simulations and experiment results of the simplified HTV2 models flying at hypervelocity.
测量波段 实验状态 姿态角/(°) 模型及等离子体鞘套总RCS 压力/kPa 速度/(km·s–1) 俯仰角 偏航角 滚转角 测量/dBsm 计算/dBsm 误差/dB X 8.2 5.0 18.93 17.11 109.95 –3.30 –1.51 1.79 C 10.4 5.0 –7.15 –2.10 69.4 0.73 3.14 2.41 X 11.6 5.0 4.03 –3.25 107.73 –3.94 –5.20 1.26 X 14.8 5.0 –7.00 9.10 70.93 –9.26 –6.01 3.35 -
[1] Akey N D, Cross A E 1970 NASA-TN-D-5615 [1970-02-01]
[2] Dirsa E F 1960 Pro. IRE 48 703
Google Scholar
[3] Golden K E, Pridmore D C, Stewart G E 1970 NASA 19710011626 (Washington: NASA LANGLEY Research Center)
[4] Wood G E, Asmar S W, Rebold T A 1997 TDA Progress Report 42 131
[5] Bachynski M P, Gibbs B W 1970 NASA 19710011649 (Washington: NASA LANGLEY Research Center)
[6] Sotnikov V L, Leboeuf J N, Mudaliar S 2010 IEEE Trans. Plasma Sci. 38 2208
Google Scholar
[7] Usui H, Yamashita F, Matsumoto H 1999 Adv. Space Res. 24 1069
Google Scholar
[8] Mather D E, Pasqual J M, Sillence J P, Lewis P 2005 AIAA/CIRA 13th International Space Planes and Hypersonics Systems and Technology Capua, Italy, May 16–20, 2005 p3443
[9] 周超, 张小宽, 张晨新, 吴国成 2014 现代雷达 36 83
Google Scholar
Zhou C, Zhang X K, Zhang C X, Wu G C 2014 Modern Radar 36 83
Google Scholar
[10] 吴巍, 刘方, 钟建林, 王国宏 2019 电波科学学报 34 610
Wu W, Liu F, Zhong J L, Wang G H 2019 Chin. J. Radio. Sci. 34 610
[11] 金铭, 韦笑, 吴洋, 张羽淮, 余西龙 2015 物理学报 64 205205
Google Scholar
Jin M, Wei X, Wu Y, Zhang Y H, Yu X L 2015 Acta Phys. Sin. 64 205205
Google Scholar
[12] 马平, 石安华, 杨益兼, 于哲峰, 黄洁 2015 强激光与粒子束 27 073201
Google Scholar
Ma P, Shi A H, Yang Y J, Yu Z F, Huang J 2015 High Power Laser Part. Beams 27 073201
Google Scholar
[13] 马平, 石安华, 杨益兼, 于哲峰, 梁世昌, 黄洁 2017 物理学报 66 102401
Google Scholar
Ma P, Shi A H, Yang Y J, Yu Z F, Liang S C, Huang J 2017 Acta Phys. Sin. 66 102401
Google Scholar
[14] 邾继贵, 于之靖 2012 视觉测量原理与方法 (北京: 机械工业出版社) 第4678页
Zhu J G, Yu Z J 2012 Principles and Methods of Visual Measurements (Beijing: China Machine Press) pp46–78 (in Chinese)
[15] Sun H Y, Cui Z W, Wang J J, Han Y P, Shi P 2018 Phys. Plasmas 25 063514
Google Scholar
[16] Sun H Y, Wang J J, Han Y P, Cui Z W, Sun P, Shi X W, Zhao W J 2018 Int. J. Antennas Propag. 1 14
[17] 孙浩宇 2018 博士学位论文 (西安: 西安电子科技大学)
Sun H Y 2018 Ph. D. Dissertation (Xi’an: Xidian University) (in Chinese)
[18] 艾夏 2013 博士学位论文 (西安: 西安电子科技大学)
Ai X 2013 Ph. D. Dissertation (Xi’an: Xidian University) (in Chinese)
[19] 陈安涛 2019 博士学位论文 (西安: 西安电子科技大学)
Chen A T 2019 Ph. D. Dissertation (Xi’an: Xidian University) (in Chinese)
[20] 葛德彪, 闫玉波 2005 电磁波时域有限差分方法 (西安: 西安电子科技大学出版社) 第88—89页
Ge D B, Yan Y B 2005 Finite-Difference Time-Domain Method for Electromagnetic Waves (Xi’an: Xidian University Press) pp88–89 (in Chinese)
[21] Ai X, Han Y, Li C Y, Shi X W 2011 Prog. Electromagn. Res. Lett. 22 83
Google Scholar
[22] Ai X, Han Y, Chen Z, Shi X W 2011 Prog. Electromagn. Res. M. 18 143
Google Scholar
[23] Ai X, Tian Y, Han Y P, Shi X W, Li W T 2013 J. Quant. Spectrosc. Radiat Transfer 124 28
Google Scholar
[24] 陈伟, 郭立新, 李江挺, 淡荔 2017 物理学报 66 084102
Google Scholar
Chen W, Guo L X, Li J T, Dan L 2017 Acta Phys. Sin. 66 084102
Google Scholar
计量
- 文章访问数: 6557
- PDF下载量: 104