1.   引　言

由于在空间光通信、量子纠缠、粒子操控等领域的应用, 光学涡旋已成为近年来的研究热点. 与普通光束相比, 涡旋光束可以携带轨道角动量. 拉盖尔-高斯(L-G)光束就是一种典型的涡旋光束, 其具有螺旋相位结构和环形场强分布, 光束中心的强度为零^[1]. L-G光束沿其传播方向携带每光子$\ell \hbar$的轨道角动量, $\ell$表示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}}{\ell _{\text{1}}}$和${{ - 2}}{\ell _{\text{2}}}$的拓扑荷; 处于中间位置的旋转腔镜是完全反射的, 可以给反射光增加$+2{\ell _{\text{1}}}$和$+2{\ell _{\text{2}}}$的拓扑荷, 并可以在支架上绕腔轴$z$旋转(角平衡位置${\phi _{{0}}}=0$); 同时两个腔场通过光纤线性耦合在一起. 两个频率均为${\omega _{\text{c}}}$, 振幅分别为${\varepsilon _{\text{c}}}$和${\varepsilon _{\text{d}}}$的强驱动场分别从两侧注入系统驱动频率为${\omega _0}$的L-G腔模. 利用一束频率为${\omega _{\text{p}}}$、振幅为${\varepsilon _{\text{L}}}$(${\varepsilon _{\text{R}}}$)的弱探测场从系统左侧(右侧)入射, 从而检查该光旋转系统相对于左侧(右侧)探测光的响应特性. 在驱动场频率${\omega _{\text{c}}}$的旋转框架下, 系统的哈密顿量($\hbar$ = 1)可写为

$$\begin{split} \;& H = {\varDelta _{\text{c}}}(c_1^\dagger {c_1} + c_2^\dagger {c_2}) + \frac{{L_z^2}}{{2I}} + \frac{1}{2}I\omega _\phi ^2{\phi ^2} \\ &~~+ {g_{\phi 1}}c_1^\dagger {c_1}\phi + {g_{\phi 2}}c_2^\dagger {c_2}\phi + J(c_1^\dagger {c_2} + c_2^\dagger {c_1}) \\ &~~ + {\text{i}}({\varepsilon _{\text{c}}}c_1^\dagger - \varepsilon _{\text{c}}^ * {c_1})+ {\text{i}}({\varepsilon _{\text{d}}}c_2^\dagger - \varepsilon _{\text{d}}^ * {c_2}) \\ &~~ + {\text{i}}{\varepsilon _{\text{L}}}(c_1^\dagger {{\text{e}}^{ - {\text{i}}\delta t}} - {c_1}{{\text{e}}^{{\text{i}}\delta t}}) + {\text{i}}{\varepsilon _{\text{R}}}(c_2^\dagger {{\text{e}}^{ - {\text{i}}\delta t}} - {c_2}{{\text{e}}^{{\text{i}}\delta t}}) \text{, } \end{split}$$

(1)

式中, 第1项表示L-G腔模的自由哈密顿量, 其中${c_1} (c_1^\dagger )$和${c_2}(c_2^\dagger )$为两个腔的湮灭(产生)算符, ${\varDelta _{\text{c}}} = {\omega _0} - {\omega _{\text{c}}}$为驱动场与腔模的失谐; 第2项和第3项给出了旋转镜的自由哈密顿量, ${\omega _\phi }$, ${L_z}$和$\phi$分别为旋转镜的角频率、角动量和角位移, $I = M{R^2}/2$为旋转镜的转动惯量, 其中$M$和$R$分别为旋转镜的质量和半径; 第4项和第5项表示L-G腔模与旋转镜的光旋转耦合, 单光子耦合强度分别为${g_{\phi 1}} = c{\ell _1}/L$和${g_{\phi 2}} = c{\ell _2}/L$, 其中L为两个腔的长度; 第6项表示两个腔模的线性耦合, 其中$J$为线性耦合强度; 最后4项分别表示强驱动场和弱探测场对腔模的驱动, $\delta = {\omega _{\text{p}}} - {\omega _{\text{c}}}$为探测场与驱动场之间的失谐.

为了方便计算, 设定$b = \sqrt {\dfrac{{I{\omega _\phi }}}{2}} \Big(\phi + \dfrac{{\text{i}}}{{I{\omega _\phi }}}{L_z}\Big)$和${b^\dagger } = \sqrt {\dfrac{{I{\omega _\phi }}}{2}} \Big(\phi - \dfrac{{\text{i}}}{{I{\omega _\phi }}}{L_z}\Big)$分别表示转动腔镜的湮灭(产生)算符, 系统的哈密顿量可以写为

$$\begin{split} H = \;&{\varDelta _{\text{c}}}(c_1^\dagger {c_1} + c_2^\dagger {c_2}) + {\omega _\phi }{b^\dagger }b \\ &+ ({g_{\phi 1}}c_1^\dagger {c_1} + {g_{\phi 2}}c_2^\dagger {c_2})({b^\dagger } + b)\\ & + J(c_1^\dagger {c_2} + c_2^\dagger {c_1}) + {\text{i}}({\varepsilon _{\text{c}}}c_1^\dagger - {\varepsilon _{\text{c}}}^ * {c_1})\\ & + {\text{i}}({\varepsilon _{\text{d}}}c_2^\dagger - \varepsilon _{\text{d}}^ * {c_2}) + {\text{i}}{\varepsilon _{\text{L}}}(c_1^\dagger {{\text{e}}^{ - {\text{i}}\delta t}} - {c_1}{{\text{e}}^{{\text{i}}\delta t}}) \\ &+ {\text{i}}{\varepsilon _{\text{R}}}(c_2^\dagger {{\text{e}}^{ - {\text{i}}\delta t}} - {c_2}{{\text{e}}^{{\text{i}}\delta t}}) .\end{split}$$

(2)

由系统的哈密顿量, 可以写出描述系统的动力学的相关算符的海森伯-郎之万方程:

$$\begin{split} &{{\dot c}_{1}}=- \Big[\text{i}{\varDelta }_{\text{c}}+\frac{{\kappa }_{1}}{2}+\text{i}{g}_{\phi 1}({b}^{†}+b)\Big]{c}_{1}\\ &\qquad+{\varepsilon }_{\text{c}}+{\varepsilon }_{\text{L}}{\text{e}}^{-\text{i}\delta t}-\text{i}J{c}_{2} \text{, } \\ & {{\dot c}_{2}}=-\Big[\text{i}{\varDelta }_{\text{c}} + \frac{{\kappa }_{2}}{2}+\text{i}{g}_{\phi 2}({b}^{†}+b)\Big]{c}_{2}\\ &\qquad+{\varepsilon }_{\text{d}}+{\varepsilon }_{\text{R}}{e}^{-\text{i}\delta t}-\text{i}J{c}_{1} \text{, } \\ &\dot{b}=-\text{i}{\omega }_{m}b-\frac{{\gamma }_{\phi }}{2}b-\text{i}({g}_{\phi 1}{c}_{1}^{†}{c}_{1}+{g}_{\phi 2}{c}_{\text{2}}^{†}{c}_{\text{2}}) \text{, } \end{split}$$

(3)

其中${\kappa _j}$和${\gamma _\phi }$分别为两腔和转动腔镜的弛豫速率. 在乘积算符的平均值满足分解假设$\left\langle {b{c_j}} \right\rangle = \left\langle b \right\rangle \left\langle {{c_j}} \right\rangle$且忽略探测场的情况下, 可得到算符的稳态解:

$$\begin{split} &{b_{\text{s}}} = \frac{{ - {\text{i}}({g_{\phi 1}}{{\left| {{c_{1{\text{s}}}}} \right|}^2} + {g_{\phi 2}}{{\left| {{c_{2{\text{s}}}}} \right|}^2})}}{{\dfrac{{{\gamma _\phi }}}{2} + {\text{i}}{\omega _\phi }}} \text{, } \\ &{c_{1{\text{s}}}} = \frac{{\left(\dfrac{{{\kappa _2}}}{2} + {\text{i}}{\varDelta _2}\right){\varepsilon _{\text{c}}} - {\text{i}}J{\varepsilon _{\text{d}}}}}{{{J^2} + \left(\dfrac{{{\kappa _1}}}{2} + {\text{i}}{\varDelta _1}\right)\left(\dfrac{{{\kappa _2}}}{2} + {\text{i}}{\varDelta _2}\right)}} \text{, } \\ &{c_{{\text{2s}}}} = \frac{{\left(\dfrac{{{\kappa _1}}}{2} + {\text{i}}{\varDelta _1}\right){\varepsilon _{\text{d}}} - {\text{i}}J{\varepsilon _{\text{c}}}}}{{{J^2} + \left(\dfrac{{{\kappa _1}}}{2} + {\text{i}}{\varDelta _1}\right)\left(\dfrac{{{\kappa _2}}}{2} + {\text{i}}{\varDelta _2}\right)}} \text{, }\end{split}$$

(4)

其中${\varDelta _{1, 2}} = {\varDelta _{\text{c}}} + {g_{\phi 1, \phi 2}}({b_{\text{s}}} + b_{\text{s}}^{\text{*}})$为L-G腔模和驱动场之间的有效失谐.

在强驱动场的驱动下, 每个算符可以由其平均值和量子涨落的和表示, 即${c_j} = {c_{j{\text{s}}}} + {\text{δ}}{c_j}$, $b = {b_{\text{s}}} + {\text{δ}}b$. 将${c_j}$和$b$的表达式代入方程(3), 可以得到线性化的海森伯-郎之万方程:

$$\begin{split} & {\text{δ}\dot{c}_{1}}=-\left(\frac{{\kappa }_{1}}{2}+\text{i}{\varDelta }_{1}\right)\text{δ}{c}_{1}-\text{i}{g}_{\phi 1}{c}_{1s}(\text{δ}{b}^{†}+\text{δ}b)+{\varepsilon }_{\text{L}}{e}^{-\text{i}\delta t}-\text{i}J\text{δ}{c}_{2} \text{, } \\ & {\text{δ}\dot{c}_{2}}=-\left(\frac{{\kappa }_{2}}{2}+\text{i}{\varDelta }_{2}\right)\text{δ}{c}_{2}-\text{i}{g}_{\phi 2}{c}_{\text{2s}}(\text{δ}{b}^{†}+\text{δ}b)+{\varepsilon }_{\text{R}}{\text{e}}^{-\text{i}\delta t}-\text{i}J\text{δ}{c}_{1} \text{, } \\ & {\text{δ}\dot b}=-\left(\frac{{\gamma }_{\phi }}{2}+\text{i}{\omega }_{\phi }\right)\text{δ}b-\text{i}{g}_{\phi 1}({c}_{1\text{s}}\text{δ}{c}_{1}^{†}+{c}_{1\text{s}}^{\ast }\text{δ}{c}_{1})-\text{i}{g}_{\phi 2}({c}_{2\text{s}}\text{δ}{c}_{2}^{†}+{c}_{2\text{s}}^{\ast }\text{δ}{c}_{2}) \text{, } \end{split}$$

(5)

采用变换${\text{δ}}b \to {\text{δ}}b{{\text{e}}^{ - {\text{i}}{\omega _\phi }t}}$, ${\text{δ}}{c_1} \to {\text{δ}}{c_1}{{\text{e}}^{ - {\text{i}}{\varDelta _1}t}}$, ${\text{δ}}{c_2} \to {\text{δ}}{c_2}{{e} ^{ - {\text{i}}{\varDelta _2}t}}$, 则方程(5)为

$$\begin{split} & {\text{δ}\dot{c}_{1}}=-\frac{{\kappa }_{1}}{2}\text{δ}{c}_{1}-\text{i}{G}_{1}\left[\text{δ}{b}^{†}{\text{e}}^{\text{i}({\omega }_{\phi }+{\varDelta }_{1})t}+\text{δ}b{\text{e}}^{-\text{i}({\omega }_{\phi }-{\varDelta }_{1})t}\right]+{\varepsilon }_{\text{L}}{\text{e}}^{-\text{i}(\delta -{\varDelta }_{1})t}-\text{i}J\text{δ}{c}_{2}{\text{e}}^{\text{i}({\varDelta }_{1}-{\varDelta }_{2})t} \text{, } \\ & {\text{δ}\dot{c}_{2}}=-\frac{{\kappa }_{2}}{2}\text{δ}{c}_{2}-\text{i}{G}_{2}{\mathrm{e}}^{\text{i}\theta }\left[\text{δ}{b}^{†}{\mathrm{e}}^{\text{i}({\omega }_{\phi }+{\varDelta }_{2})t}+\text{δ}b{\mathrm{e}}^{-\text{i}({\omega }_{\phi }-{\varDelta }_{2})t}\right]+{\varepsilon }_{\text{R}}{\mathrm{e}}^{-\text{i}(\delta -{\varDelta }_{2})t}-\text{i}J\text{δ}{c}_{1}{\mathrm{e}}^{\text{i}({\varDelta }_{2}-{\varDelta }_{1})t} \text{, } \\ & {\text{δ}\dot b}=-\frac{{\gamma }_{\phi }}{2}\text{δ}b-\text{i}{G}_{1}\left[\text{δ}{c}_{1}{\text{e}}^{\text{i}({\omega }_{\phi }-{\varDelta }_{1})t}+\text{δ}{c}_{1}^{†}{\text{e}}^{\text{i}({\omega }_{\phi }+{\varDelta }_{1})t}\right]-\text{i}{G}_{2}\left[{\text{e}}^{-\text{i}\theta }\text{δ}{c}_{2}{\text{e}}^{\text{i}({\omega }_{\phi }-{\varDelta }_{2})t}+{\text{e}}^{\text{i}\theta }\text{δ}{c}_{2}^{†}{\text{e}}^{\text{i}({\omega }_{\phi }+{\varDelta }_{2})t}\right] \text{, }\end{split}$$

(6)

其中${G_1}$, ${G_2}$为两个涡旋腔场与旋转腔镜之间的有效光旋转耦合强度, ${G_1} = {g_{\phi 1}}{c_{1{\text{s}}}}$, ${G_2} = {g_{\phi 2}}{c_{2{\text{s}}}}{{\text{e}}^{ - {\text{i}}\theta }}$, ${g_{\phi 1}}{c_{1{\text{s}}}}$和${g_{\phi 2}}{c_{2{\text{s}}}}$之间的非互易相位差$\theta$可以通过驱动场${\varepsilon _{\text{c}}}$和${\varepsilon _{\text{d}}}$的相对相位来调节.

如果每个强驱动场都处于力学红边带(${\varDelta _1} \approx {\varDelta _2} \approx {\omega _\phi }$), 同时旋转镜频率${\omega _\phi }$远大于有效耦合强度${G_j}$, 方程(6)可以化简为

$$\begin{split} & {\text{δ}\dot{c}_{1}}=-\frac{{\kappa }_{1}}{2}\text{δ}{c}_{1}-\text{i}{G}_{1}\text{δ}b-\text{i}J\text{δ}{c}_{2}+{\varepsilon }_{\text{L}}{\text{e}}^{-\text{i}xt} \text{, } \\ & {\text{δ}\dot{c}_{2}}=-\frac{{\kappa }_{2}}{2}\text{δ}{c}_{2}-\text{i}{G}_{2}{e}^{\text{i}\theta }\text{δ}b-\text{i}J\text{δ}{c}_{1}+{\varepsilon }_{\text{R}}{\text{e}}^{-\text{i}xt} \text{, } \\ & {\text{δ}\dot b}=-\frac{{\gamma }_{\phi }}{2}\text{δ}b-\text{i}{G}_{1}\text{δ}{c}_{1}-\text{i}{G}_{2}{\text{e}}^{-\text{i}\theta }\text{δ}{c}_{2} \text{, } \end{split}$$

(7)

其中$x = \delta - {\omega _\phi }$. 由方程(7)的具体形式, 可以假设方程(7)的解具有下面的形式:

$$\begin{split} &{\text{δ}}{c_1} = {\text{δ}}{c_{1 + }}{{\text{e}}^{ - {\text{i}}xt}} + {\text{δ}}{c_{1 - }}{{\text{e}}^{{\text{i}}xt}} \text{, }\\ &{\text{δ}}{c_2} = {\text{δ}}{c_{2 + }}{{\text{e}}^{ - {\text{i}}xt}} + {\text{δ}}{c_{2 - }}{{\text{e}}^{{\text{i}}xt}} \text{, } \\ &{\text{δ}}b = {\text{δ}}{b_ + }{{\text{e}}^{ - {\text{i}}xt}} + {\text{δ}}{b_ - }{{\text{e}}^{{\text{i}}xt}} . \end{split}$$

(8)

将方程(8)代入方程(7)得

$$\begin{split} &{\text{δ}}{c_{1 + }} = \frac{{(8G_2^2 + 2{\gamma _x}{\kappa _{2x}}){\varepsilon _{\text{L}}} - (8{G_1}{G_2}{{\text{e}}^{ - {\text{i}}\theta }} + 4{\text{i}}J{\gamma _x}){\varepsilon _{\text{R}}}}}{{\varOmega - 16{\text{i}}{G_1}{G_2}J\cos \theta }} \text{, } \\ &{\text{δ}}{c_{2 + }} = \frac{{(8G_1^2 + 2{\gamma _x}{\kappa _{1x}}){\varepsilon _{\text{R}}} - (8{G_1}{G_2}{{\text{e}}^{{\text{i}}\theta }} + 4{\text{i}}J{\gamma _x}){\varepsilon _{\text{L}}}}}{{\varOmega - 16{\text{i}}{G_1}{G_2}J\cos \theta }} \text{, } \\ &{b_ + } = \frac{{(4{\text{i}}{G_1}{\kappa _{2x}} + 8{G_2}J{{\text{e}}^{ - {\text{i}}\theta }}){\varepsilon _{\text{L}}} + (8{G_1}J + 4{\text{i}}{G_2}{\kappa _{1x}}{{\text{e}}^{ - {\text{i}}\theta }}){\varepsilon _{\text{R}}}}}{{16{\text{i}}{G_1}{G_2}J\cos \theta - \varOmega }} \text{, } \end{split}$$

(9)

其中${\gamma _x} = \gamma - 2{\text{i}}x$, ${\kappa _{jx}} = {\kappa _j} - 2{\text{i}}x$, $\varOmega = 4 G_2^2{\kappa _{1 x}} + 4 G_1^2{\kappa _{2 x}} + 4{J^2}{\gamma _x} + {\kappa _{1 x}}{\kappa _{2 x}}{\gamma _x}$.

为了研究涡旋光束传输的非互易性, 需要得到$\varepsilon _{\text{L}}^{{\text{out}}}$和$\varepsilon _{\text{R}}^{{\text{out}}}$. 根据输入输出关系^[47,48]得

$$\begin{split} &{\varepsilon }_{\text{L}}^{\text{out}}+{\varepsilon }_{\text{L}}^{\text{in}}{\text{e}}^{-\text{i}xt}=\sqrt{{\kappa }_{1}}\text{δ}{c}_{1}\text{, }\\ &{\varepsilon }_{\text{R}}^{\text{out}}+{\varepsilon }_{\text{R}}^{\text{in}}{\text{e}}^{-\text{i}xt}=\sqrt{{\kappa }_{2}}\text{δ}{c}_{2}\text{, }\end{split}$$

(10)

其中$\varepsilon _{{\text{L, R}}}^{{\text{in}}} = {\varepsilon _{{\text{L, R}}}}/\sqrt {{\kappa _{\text{1}}}{\kappa _{\text{2}}}}$, 输出场具有${\text{δ}}{c_1} = {\text{δ}}{c_{1 + }}{{\text{e}}^{ - {\text{i}}xt}} + {\text{δ}}{c_{1 - }}{{\text{e}}^{{\text{i}}xt}}$和${\text{δ}}{c_2} = {\text{δ}}{c_{2 + }}{{\text{e}}^{ - {\text{i}}xt}} + {\text{δ}}{c_{2 - }}{{\text{e}}^{{\text{i}}xt}}$的形式, 由方程(9)和方程(10)得

$$\begin{split} &\varepsilon _{{\text{L + }}}^{{\text{out}}} = \sqrt {{\kappa _1}} {\text{δ}}{c_{1 + }} - {\varepsilon _{\text{L}}}/\sqrt {{\kappa _{\text{1}}}} \text{, } \\ &\varepsilon _{{\text{R + }}}^{{\text{out}}} = \sqrt {{\kappa _2}} {\text{δ}}{c_{2 + }} - {\varepsilon _{\text{R}}}/\sqrt {{\kappa _{\text{2}}}} \text{, }\\ &\varepsilon _{{\text{L}}-}^{{\text{out}}} = 0 \text{, }~~ \varepsilon _{{\text{R }}-}^{{\text{out}}} = 0 . \end{split}$$

(11)

本文定义${T_{12}}$(${T_{21}}$)为入射腔${c_1}({c_2})$的探测场通过系统后输出时的传输振幅, 根据(9)式和(11)式计算得

$$\begin{split} {T}_{12}=\;&{T}_{\rm{L\to R}}={\left|\frac{{\varepsilon }_{\text{R+}}^{\text{out}}}{{\varepsilon }_{\text{L}}^{\text{in}}}\right|}_{{\varepsilon }_{\text{R}}^{\text{in}}=0}\\ =\;&\frac{-\sqrt{{\kappa }_{1}{\kappa }_{2}}(8{G}_{1}{G}_{2}{\text{e}}^{\text{i}\theta }+4\text{i}J{\gamma }_{x})}{\varOmega -16\text{i}{G}_{1}{G}_{2}J\mathrm{cos}\theta } \text{, } \\ {T}_{21}=\;&{T}_{\rm{R\to L}}={\left|\frac{{\varepsilon }_{\text{L+}}^{\text{out}}}{{\varepsilon }_{\text{R}}^{\text{in}}}\right|}_{{\varepsilon }_{\text{L}}^{\text{in}}=0}\\ =\;&\frac{-\sqrt{{\kappa }_{1}{\kappa }_{2}}(8{G}_{1}{G}_{2}{\text{e}}^{-\text{i}\theta }+4\text{i}J{\gamma }_{x})}{\varOmega -16\text{i}{G}_{1}{G}_{2}J\mathrm{cos}\theta } . \end{split}$$

(12)

3.   涡旋光束传输非互易性

本节将详细说明在双拉盖尔-高斯旋转腔系统中非互易相位差$\theta$、涡旋光束所携带拓扑荷$\ell$和系统耗散对传输振幅的影响. 首先讨论非互易相位差$\theta$对传输振幅的影响.

图2给出传输振幅${T_{12}}$(红色实线)和${T_{21}}$(蓝色虚线)在不同非互易相位差$\theta$下随标准化失谐$x/{\kappa _1}$的演化. 可以看出, 相位差$\theta$对涡旋光束传输有比较大的影响, 在$\theta$ = 0时, 探测场在两个方向上的传输是互易的, 但当$\theta \ne 0$探测场的传输会呈现明显的非互易性. 观察图2(b)—图2(e)可以发现, 非互易现象会出现在共振点及附近, 还可以发现相位差的正负可决定非互易的方向($\theta < 0$时, ${T_{21}} > {T_{12}}$; 而$\theta > 0$时, ${T_{21}} < {T_{12}}$). 当$\theta = \pm {\text{π }}/4$时, 非互易性在$x = - 0.4\kappa$处呈现的最明显. 当$\theta = \pm {\text{π }}/2$时, 系统可以在共振处($x = 0$)呈现完美的涡旋光束非互易传输(${T_{12}} = 0$时${T_{21}} = 1$, 或${T_{12}} = 1$时${T_{21}} = 0$). 由(12)式也可以看出, 当$\theta = n{\text{π }}$时, ${T_{12}} = {T_{21}}$, 涡旋光束传输呈现互易性, 而当$\theta \ne n{\text{π }}$时, ${T_{12}} \ne {T_{21}}$, 即涡旋光束的传输是非互易的. 实际上, 这里的涡旋光束的非互易传输主要源于光子不同传输路径之间的量子干涉作用. 为了更清楚地展示更多非互易相位差对传输振幅的影响, 图2(f)给出了传输振幅${T_{12}}$(红色实线)和${T_{21}}$(蓝色虚线)随非互易相位差$\theta$的演化.

接下来讨论系统耗散对涡旋光束传输非互易性的影响. 为不失一般性, 取非互易相位差$\theta = {\text{π }}/2$. 图3(a)和图3(b)给出了传输振幅${T_{12}}$和${T_{21}}$在旋转腔镜的不同耗散率$\gamma$下随标准化失谐$x/{\kappa _1}$的演化. 可以看出在共振点(x = 0)附近, 随着旋转镜耗散率$\gamma$增加, 从腔${c_1}$到腔${c_2}$的传输振幅${T_{12}}{\text{ = }}1$保持不变, 而从腔${c_2}$到腔${c_1}$的传输振幅${T_{21}}$呈现出明显的变化, 显示出先降低后增加的趋势. 当旋转镜耗散率$\gamma {\text{ = }}{\kappa _{\text{1}}}{\text{ = }}1$时, 涡旋光束的传输非互易性达到最佳(${T_{12}} = {\text{1}}$, ${T_{21}} = {\text{0}}$, 在共振点处).

图3(c)和图3(d)给出了传输振幅${T_{12}}$和${T_{21}}$在腔${c_2}$的不同耗散率下随标准化失谐$x/{\kappa _1}$的演化. 可以看出, 在共振点(x = 0)附近, 随着腔${c_2}$耗散率${\kappa _2}$增加, 从腔${c_2}$到腔${c_1}$的传输振幅${T_{21}}{\text{ = 0}}$保持不变, 而从腔${c_1}$到腔${c_2}$的传输振幅${T_{12}}$呈现逐渐减小的趋势. 当腔${c_2}$耗散率${\kappa _2}{\text{ = }}{\kappa _{\text{1}}}{\text{ = }}1$时, 涡旋光束传输的非互易性达到最佳(${T_{12}} = {\text{1}}$, ${T_{21}} = {\text{0}}$, 在共振点处). 所以对于涡旋光束的非互易传输来说, 腔场${c_2}$的耗散不是越低越好.

接下来讨论两个涡旋光场所携带拓扑荷比值${\ell _2}/{\ell _1}$对传输振幅的影响. 通常, 涡旋光场所携带拓扑荷实际上会影响旋转腔镜与L-G腔模的耦合强度. 为了方便, 这里通过调节驱动场振幅${\varepsilon _{\text{c}}}$和${\varepsilon _{\text{d}}}$的大小从而使得拓扑荷比${\ell _2}/{\ell _1}$与有效耦合强度比${G_2}/{G_1}$相等. 图4(a)—图4(f)给出了在标准化失谐x = 0时传输振幅${T_{12}}$和${T_{21}}$在不同的非互易相位差下随两个涡旋光场所携带拓扑荷比${\ell _2}/{\ell _1}$的演化. 可以发现, 当$\theta = 0$和$\theta = {\text{π }}$时, 拓扑荷的比值对传输振幅${T_{12}}$和${T_{21}}$的影响是相同的, 即系统对涡旋光束的传输呈现互易性. 当$\theta = \pm {\text{π }}/4$和$\theta = \pm {\text{π }}/2$时, 系统对涡旋光束传输呈现明显的非互易性. 当$\theta = \pm {\text{π }}/4$时, 随着拓扑荷比值的增加, 传输振幅${T_{12}}$和${T_{21}}$的变化趋势是相同的, 且是偶对称, 但相同拓扑荷比值所对应的传输振幅${T_{12}}$和${T_{21}}$是不同的, 在拓扑荷比值为$\pm 1$处的非互易性是最好的. 当$\theta = \pm {\text{π }}/2$时, 随着拓扑荷比值的增加, 传输振幅${T_{12}}$仍呈现关于拓扑荷比值等于$0$处的对称性, 但此时传输振幅${T_{21}}$保持为$0$不变. 在拓扑荷比值为$\pm 1$处, 系统对涡旋光束传输可以呈现完美的非互易性(${T_{12}} = 0$时${T_{21}} = 1$, 或${T_{12}} = 1$时${T_{21}} = 0$). 此外, 拓扑荷的正负代表涡旋光矢量的旋转方向, 根据涡旋光矢量的旋转方向可以将涡旋光束分为左旋涡旋光束和右旋涡旋光束. 从图4(d)和图4(e)可以了解到涡旋光束想要呈现完美的非互易性与涡旋光束的左、右旋无关, 只要保证$\left| {{\ell _2}/{\ell _1}} \right|{\text{ = }}1$.

4.   非互易慢光效应

在腔系统中, 光的群延迟是描述光群速度的一个重要指标, 首先介绍群延迟的定义, 即^[49,50]

$${\tau _{ij}} = \frac{{{\text{d}}{\varTheta _{21(12)}}}}{{{\text{d}}{\omega _{\text{p}}}}} ,$$

(13)

其中${\varTheta _{21(12)}}$是在频率为${\omega _{\text{p}}}$下输出场${c_1} \to {c_2}$(${c_2} \to {c_1}$)的相位, 群延迟${\tau _{ij}} < 0$对应快光, ${\tau _{ij}} > 0$对应慢光.

图5(a)—图5(c)给出了针对涡旋光束所对应的群延迟${\tau _{12}}$和${\tau _{21}}$在不同的非互易相位差下随两个涡旋光场所携带拓扑荷比${\ell _{\text{2}}}{\text{/}}{\ell _{\text{1}}}$的演化. 在$\theta = 0$时群速度是互易的且关于拓扑荷比值偶对称. 在$\theta = \pm {\text{π }}/2$时, 系统可以呈现明显的非互易慢光效应, 且相位的变化可以影响非互易慢光的群延迟的变换. 随着拓扑荷比值的增加, ${\tau _{12}}$和${\tau _{21}}$均呈现先增大后减小的趋势且关于拓扑荷比值偶对称, 但群延迟的差值(${\tau _{12}}-{\tau _{21}}$或${\tau _{21}}-{\tau _{12}}$)一直保持不变.

5.   实验可行性分析

最后, 根据目前的实验进展, 对本文使用的系统方案的可行性进行了讨论. 本文讨论的双拉盖尔-高斯旋转腔系统是在文献[6]提出的由两个螺旋相位元件所构成的拉盖尔-高斯旋转腔的基础上, 添加一个固定的螺旋相位元件所实现的. 通过三个螺旋相位元件作为腔镜来组成两个光腔, 并且用光纤将两个光腔线性耦合. 其中, 两个输入腔镜FM1和 FM2为部分透明且被刚性固定, 不会改变透射光束的轨道角动量但会给反射光束增加${{ - 2}}{\ell _{\text{1}}}$和${{ - 2}}{\ell _{\text{2}}}$的拓扑荷, 这一效应已经被实验观测到^[51]. 另外, 处于中间位置的旋转腔镜RM是完全反射的, 可以给反射光增加$+2{\ell _{\text{1}}}$和$+2{\ell _{\text{2}}}$的拓扑荷^[6]. 随着纳米技术的发展, 螺旋相位元件可以利用超精密点金刚石车削机床直接加工铝盘表面得到^[6,52]. 利用螺旋相位元件的反射和透射可以改变激光光束的方位结构, 进而改变其携带的拓扑荷的值. 拓扑荷$\ell$依赖于螺旋相位元件的结构, 通过将螺旋相位元件的方位角坡道划分为离散的阶跃并控制阶跃高度和旋向, 可以设计出具有特定拓扑荷$\ell$的涡旋光场^[6,53]. 最近实验已经可以通过螺旋相位元件产生携带拓扑荷高达1000的L-G光束^[54]. 对于本工作, 若想实现理想的涡旋光隔离器, 就需要在本系统中实现涡旋光束的完美非互易性, 即涡旋光束的单向传输. 在理论分析中得到, 若要呈现完美的非互易涡旋光束传输, 需要保证$\left| {{\ell _2}/{\ell _1}} \right|{\text{ = }}1$, 因此这里的涡旋光束所携带的拓扑荷需是可以调节的. 在实验上, 可以通过调整旋转腔镜两侧的阶跃高度和旋向来实现. 另外, 理论分析中得到非互易相位差对涡旋光束的非互易传输有着重要影响, 而非互易相位差的改变在实验上可以通过调节两个强驱动激光场${\varepsilon _{\text{c}}}$和${\varepsilon _{\text{d}}}$的相对相位来实现^[45]. 值得强调的是, 本文所构建的系统在实验实现时需要较高精细度的光腔, 在光腔的精细度比较低时, 在该系统所要探究相关现象会不明显^[53].

6.   结　论

本文研究了涡旋光束在双拉盖尔-高斯旋转腔系统中的非互易传输. 将以往大多在光机械系统中研究的普通光束的非互易性推广到在光旋转系统中研究携带轨道角动量的涡旋光束传输的非互易性. 系统中, 光旋转耦合作用以及光纤的线性耦合作用提供的两条路径, 从而可以产生量子干涉. 首先, 研究了非互易相位差对传输振幅的影响, 发现相位差可以决定涡旋光束非互易性的发生及方向. 然后, 在相位差$\theta = {\text{π }}/2$时, 研究了系统耗散对传输振幅的影响, 分析了涡旋光束实现完美的非互易传输时系统的耗散应该满足的条件. 接下来, 分析了涡旋光束携带的拓扑荷比值对系统非互易性的影响, 探究了实现涡旋光束完美的非互易传输时拓扑荷比值应该满足的条件, 发现涡旋光的左、右旋不会对传输产生影响. 拓扑荷对光传输产生影响, 主要是因为涡旋光束携带的拓扑荷数会影响旋转腔镜与腔模的耦合强度. 最后, 分析了该系统产生的非互易的慢光效应. 这些研究成果可用于设计针对携带轨道角动量的涡旋光束的理想光隔离器, 有望应用于光通信等领域.