摘要:光学非互易性对建立量子网络和现代通讯技术都是不可或缺的. 本文研究了双腔光力学系统在蓝失谐驱动下如何实现完美的光学非互易性. 研究发现此系统中的光学非互易性来源于系统中的光力相互作用和腔模线性耦合相互作用之间的量子干涉效应. 在应用光力学输入输出关系得出输出光场表达式后, 给出了在此系统中实现完美光学非互易性的条件以及影响非互易谱线宽度的因素. 另外还发现当系统参数(耗散速率)一定时, 可以存在两套耦合强度来实现完美的光学非互易性. 最后利用劳斯-霍尔维茨(Routh-Hurwitz)稳定性判据给出了系统的稳定条件.
关键词: 腔光力学/
光学非互易性/
光频隔离器

English Abstract

全文HTML

--> --> -->

非互易光学器件在光源和接收器调换位置后可以使光信号表现出不同的传输特性. 由于其可以抑制多余的信号, 因此在量子信号处理和量子通信中有着重要的应用. 例如, 在量子超导电路中它们可以保护信号源不被读取器件发出的噪声干扰^[1]. 为实现光学非互易性, 时间反演对称性破缺是必须的. 传统的非互易性光学器件都是依赖 强磁场去实现时间反演对称性破缺^[2]. 然而由于需要较强磁场, 这些传统器件体积往往较大, 不便于微型化和集成化. 近年来由于纳米技术的进步, 使微纳系统中光学现象得到广泛研究^[3-7]. 腔光力学系统中的光辐射压力可以使系统呈现出各种有趣的量子现象. 例如, 腔光力学系统中的量子纠缠^[8-16], 力学振子的基态冷却^[17-21], 光力诱导透明^[22-26]以及非线性效应^[26-33]和声子阻塞^[34]等量子现象. 最近, 人们意识到光力耦合相互作用也可以产生光学非互易传输现象. 例如通过光力相互作用可以产生非互易光学反应在理论上被预言^[35-37], 并在实验上得到证实^[38-42]. 并在理论上指出如果采用适当的驱动场, 以力学模为中介的两个腔模之间的态转换可以是非互易的^[43-45], 以及通过光力耦合相互作用可以实现信号非互易放大现象^[46-48]. 在文献[49,50]中, 理论上给出了通过光力相互作用可以实现非互易光子阻塞效应, 以及在文献[51]中, 作者理论上指出通过光力相互作用可以实现非互易慢光. 另外, 在文献[52,53]中, 作者理论上预言了通过光力耦合可以实现声子环形器和热二极管. 然而在大部分文献当中人们常常采用红失谐的驱动场和$ \pm\dfrac{{\text{π}}}{2} $的非互易相位差去实现光学非互易性.
本文研究了在蓝失谐驱动下, 在双腔光力系统(如图1所示)中如何实现光场的非互易传输. 在此模型中, Li等^[48]利用力学驱动机制实现了光的非互易放大, 并在红失谐驱动下, 利用$ \pm\dfrac{{\text{π}}}{2} $的非互易相位差去实现光学非互易性^[43]. 实际上, 此系统中的光学非互易性源于光力耦合和腔模线性耦合的共同作用, 使从不同路径传输的光信号之间产生干涉效应. 本文根据此物理机理并由腔光力学中标准的光场输入输出关系, 得到了实现完美的非互易光传输条件. 研究发现, 在系统中各耗散速率一定的情况下, 会有两套耦合强度可以实现光学非互易传输, 并且即使在非互易相位差不为$ \pm\dfrac{{\text{π}}}{2} $时系统依然可以实现完美光学非互易性. 最后根据劳斯-霍尔维茨(Routh-Hurwitz)稳态判据给出了系统在蓝失谐驱动下的稳定条件. 这些研究结果有望能应用于在光力系统中实现光频隔离器、非互易态转换等量子信息处理过程.
图 1 双腔光力学系统示意图, 两光学腔通过光力相互作用与一个力学振子相耦合, 振幅为$\varepsilon_{\rm c}$和 $\varepsilon_{\rm d}$($\varepsilon_{\rm L}$和$\varepsilon_{\rm R}$)的强耦合场 (探测场)分别从左右两侧驱动腔模$c_{1}$和$c_{2}$, 同时两腔模之间存在线性耦合相互作用J
Figure1. A two-cavity optomechanical system with a mechanical resonator interacted with two cavities. Two strong coupling fields (probe fields) with amplitudes $\varepsilon_{\rm c}$ and $\varepsilon_{\rm d}$ ($\varepsilon_{\rm L}$ and $\varepsilon _{\rm R}$) are used to drive cavity $c_{1}$ and $c_{2}$ respectively. Meanwhile, the two cavities are linearly coupled to each other with coupling strength J

本文研究了一个双腔光力学系统, 如图1所示, 左右两个光学腔与中间一个力学振子通过光力相互作用耦合; $ {c}_{i} $($ \omega_{0} $)和$ b $($ \omega_{m} $)分别表示光学腔$ i $和力学振子的湮灭算符(本征频率); $ \kappa_{i} $和$ \gamma $分别表示光腔i和力学振子的弛豫速率. 两个频率均为$ \omega_{\rm c} $($ \omega_{\rm p} $)、振幅分别为$ \varepsilon_{\rm c} $和$ \varepsilon_{\rm d} $ ($ \varepsilon_{\rm L} $和$ \varepsilon_{\rm R} $)的强驱动场(弱探测场) 分别从左右两侧射入并驱动腔模$ c_{1} $和$ c_{2} $. 同时左右两腔之间由线性相互作用$ \hbar J(c_{1}^{?}c_{2}+c_{2}^{?}c_{1}) $相耦合, J为线性耦合强度. 在相对驱动场频率$ \omega_{\rm c} $做旋转后, 系统的哈密顿量($ \hbar = 1 $)可写为

$ \begin{split} H =& \varDelta_{\rm c}(c_{1}^{?}c_{1}+c_{2}^{?}c_{2})+\omega_{\rm{m}} b^{?}b+(g_{1}c_{1}^{?}c_{1} \\ &+g_{2}c_{2}^{?}c_{2})(b^{? }+b) +J(c_{1}^{?}c_{2}+c_{2}^{?}c_{1})\\ &+{\rm i}(\varepsilon_{\rm c}c_{1}^{? }-\varepsilon^{\ast}_{\rm c}c_{1})+{\rm i}(\varepsilon_{\rm d}c_{2}^{?}-\varepsilon ^{\ast}_{\rm d}c_{2})\\ &+{\rm i}\varepsilon_{\rm L}(c_{1}^{?}{\rm e}^{-{\rm i}\varDelta t}-c_{1}{\rm e}^{{\rm i}\varDelta t}) \\ & +{\rm i}\varepsilon_{\rm R}(c_{2}^{?}{\rm e}^{-{\rm i}\varDelta t}-c_{2}{\rm e}^{{\rm i}\varDelta t}), \end{split} $

其中$ \varDelta_{\rm c} = \omega_{0}-\omega_{\rm c} $ ($ \varDelta = \omega_{\rm p}-\omega_{\rm c} $)为腔模(探测场)与驱动场之间的失谐, $ g_{i} $为光学腔i与力学振子之间的单光子耦合常数. 实际上, 此三模光力耦合系统在实验上是可行的, 如在法布里-珀罗(Fabry-Pérot)腔中加入力学膜的实验装置, 见文献[54-56].
根据海森伯-郎之万方程, 由系统哈密顿量(1)式可得系统相关算符的运动方程为:

$ \begin{split} \dot{c}_{1} =& - \Big[{\rm i}\varDelta_{\rm c}+\dfrac{\kappa_{1}}{2}+{\rm i}g_{1}(b^{?} +b)\Big]c_{1}+\varepsilon_{\rm c} \\ & +\varepsilon_{\rm L}{\rm e}^{-{\rm i} \varDelta t}-{\rm i}Jc_{2},\\ \dot{c}_{2} =& - \Big[{\rm i}\varDelta_{\rm c}+\dfrac{\kappa_{2}}{2}+{\rm i}g_{2}(b^{?} +b)\Big]c_{2}+\varepsilon_{\rm d}\\ & +\varepsilon_{\rm R}{\rm e}^{-{\rm i} \varDelta t}-{\rm i}Jc_{1},\\ \dot{b} =& -{\rm i}\omega_{\rm m}b-\dfrac{\gamma}{2}b-{\rm i}(g_{1}c_{1}^{?}c_{1}+g_{2}c_{2}^{?} c_{2}). \end{split} $

在没有探测场时, 根据假设$ \langle bc_{i}\rangle = \langle b\rangle\langle c_{i}\rangle $, 可以得出各算符的稳态平均值为:

$ \begin{split} &\langle b\rangle = b_{\rm s} = \dfrac{-{\rm i}(g_{1}\left\vert c_{1{\rm s}}\right\vert ^{2}+g_{2}\left\vert c_{2{\rm s}}\right\vert ^{2})}{\dfrac{\gamma}{2}+{\rm i}\omega_{m}},\\ &\left\langle c_{1}\right\rangle = c_{1{\rm s}} = \dfrac{\left(\dfrac{\kappa_{2}}{2}+{\rm i}\varDelta_{2}\right)\varepsilon_{\rm c}-{\rm i}J\varepsilon_{\rm d}}{J^{2}+\left(\dfrac{\kappa_{1}}{2}+{\rm i}\varDelta_{1}\right)\left(\dfrac{\kappa_{2}}{2}+{\rm i}\varDelta_{2}\right)},\\ & \left\langle c_{2}\right\rangle = c_{2{\rm s}} = \dfrac{\left(\dfrac{\kappa_{1}}{2}+{\rm i}\varDelta_{1}\right)\varepsilon_{\rm d}-{\rm i}J\varepsilon_{\rm c}}{J^{2}+\left(\dfrac{\kappa_{1}}{2}+{\rm i}\varDelta_{1}\right)\left(\dfrac{\kappa_{2}}{2}+{\rm i}\varDelta_{2}\right)}, \end{split} $

其中, $ \varDelta_{1,2} = \varDelta_{\rm c}+g_{1,2}(b_{\rm s}+b_{\rm s}^{\ast}) $表示腔模和驱动场间的有效失谐. 当探测场不为零时, 可以把算符写成这个平均值和一个小的涨落值之和, 即$ b = b_{\rm s}+{\text{δ}} b $, $ c_{1} = c_{1{\rm s}}+{\text{δ}} c_{1} $, $ c_{2} = c_{2{\rm s}}+{\text{δ}} c_{2} $. 把算符的这种相加形式代入方程(2)中, 仅保留涨落的线性项, 并采用相互作用表象, 即$ {\text{δ}} b\rightarrow{\text{δ}} b{\rm e}^{-{\rm i}\omega_{\rm m}t} $, $ {\text{δ}} c_{1}\rightarrow{\text{δ}} c_{1}{\rm e}^{-{\rm i}\varDelta_{1}t} $, $ {\text{δ}} c_{2}\rightarrow{\text{δ}} c_{2}{\rm e}^{-{\rm i}\varDelta_{2}t} $, 则可以得到线性化的海森伯-郎之万方程为:

$ \begin{split} {\text{δ}}\dot{c}_{1} = & -\dfrac{\kappa_{1}}{2}{\text{δ}} c_{1}-{\rm i}G_{1}({\text{δ}} b^{?}{\rm e}^{{\rm i}(\omega_{\rm m}+\varDelta_{1})t}+{\text{δ}} b{\rm e}^{-{\rm i}(\omega_{\rm m}-\varDelta_{1} )t})\\ &+\varepsilon_{\rm L}{\rm e}^{-{\rm i}(\varDelta-\varDelta_{1})t}-{\rm i}J{\text{δ}} c_{2}{\rm e}^{{\rm i}(\varDelta_{1}-\varDelta_{2})t},\\ {\text{δ}}\dot{c}_{2} = & -\dfrac{\kappa_{2}}{2}{\text{δ}} c_{2} \!-\! {\rm i}G_{2}{\rm e}^{{\rm i}\theta }({\text{δ}} b^{?}{\rm e}^{{\rm i}(\omega_{\rm m} \!+\! \varDelta_{2})t} \!+\! {\text{δ}} b{\rm e}^{-{\rm i}(\omega_{\rm m} \!-\! \varDelta_{2})t})\\ &+\varepsilon_{\rm R}{\rm e}^{-{\rm i}(\varDelta-\varDelta_{2})t}-{\rm i}J{\text{δ}} c_{1}{\rm e}^{{\rm i}(\varDelta_{2}-\varDelta_{1})t},\\ {\text{δ}}\dot{b} = & -\dfrac{\gamma}{2}{\text{δ}} b-{\rm i}G_{1}({\text{δ}} c_{1} {\rm e}^{{\rm i}(\omega_{\rm m}-\varDelta_{1})t}+{\text{δ}} c^{?}_{1}{\rm e}^{{\rm i}(\omega_{\rm m}+\varDelta_{1} )t})\\ & -{\rm i}G_{2}({\rm e}^{-{\rm i}\theta}{\text{δ}} c_{2}{\rm e}^{{\rm i}(\omega_{\rm m}-\varDelta_{2})t}+{\rm e}^{{\rm i}\theta }{\text{δ}} c^{?}_{2}{\rm e}^{{\rm i}(\varDelta_{2}+\omega_{m})t}), \end{split} $

其中$ G_{1} = g_{1}c_{1{\rm s}} $, $ G_{2} = g_{2}c_{2{\rm s}}{\rm e}^{-{\rm i}\theta} $. 由等式(3)可知, 通过调节驱动场$ \varepsilon_{\rm c} $和$ \varepsilon_{\rm d} $可以有效调节光力耦合$ g_{1}c_{1{\rm s}} $和$ g_{2}c_{2{\rm s}} $之间的非互易相位差θ (即调节驱动场$ \varepsilon_{\rm c} $和$ \varepsilon_{\rm d} $的强度和相位可以使$ c_{1{\rm s}} $为实数, 而此时$ c_{2{\rm s}} $的辐角便是非互易相位差θ). 为简化, 本文只讨论相等耦合$ G_{1} = G_{2} = G $和$ \kappa_{1} = \kappa_{2} = \kappa $, 并且设$ G(J) > 0 $. 本文讨论蓝失谐驱动, 即$ (\varDelta_{1}\approx\varDelta_{2}\approx-\omega_{\rm m}) $, 并假设力学振子频率$ \omega_{\rm m} $远大于耦合强度G, 则方程(4)可以化简为

$ \begin{array}{l} \dot{{{{v}}}} = {{M}}{{{v}}}+{{{v}}}_{\rm {in}}, \end{array} $

其中$ {{{v}}} \!=\! \{{\text{δ}} c_{1},{\text{δ}} c_{2},{\text{δ}} b^{?}\}^{\rm T}\!$, ${{{v}}}_{\rm {in}} \!=\! \{\varepsilon_{\rm L}{\rm e}^{-{\rm i}xt},\varepsilon_{\rm R}{\rm e}^{-{\rm i}xt},0\}^{\rm T} $ ($ x = \varDelta+\omega_{\rm m} $), 以及矩阵

$ \begin{array}{l} {{M}} = \left( \begin{array} [c]{ccc}-\dfrac{\kappa_{1}}{2} & -{\rm i}J & -{\rm i}G\\ -{\rm i}J & -\dfrac{\kappa_{2}}{2} & -{\rm i}G{\rm e}^{{\rm i}\theta}\\ {\rm i}G & {\rm i}G{\rm e}^{-{\rm i}\theta} & -\dfrac{\gamma}{2} \end{array} \right). \end{array} $

由$ {\nu}_{\rm {in}} $的形式可以假设方程(5)的解具有${\text{δ}} s =$$ {\text{δ}} s_{+}{\rm e}^{-{\rm i}xt}+{\text{δ}} s_{-}{\rm e}^{{\rm i}xt} $ ($ s = b, c_{1}, c_{2} $)的形式, 经计算可得

$ \begin{split} & {\text{δ}} b_{-}^{\ast} = \dfrac{4G[({\rm i}\kappa_{x}+2J{\rm e}^{-{\rm i}\theta})\varepsilon _{\rm L}+(2J+{\rm i}\kappa_{x}{\rm e}^{-{\rm i}\theta})\varepsilon_{\rm R}]}{4J^{2}\gamma_{x}-8G^{2}\kappa_{x}+\kappa_{x}^{2}\gamma_{x}+16{\rm i}G^{2}J\cos\theta},\\ & {\text{δ}} c_{1+} = \dfrac{(2\gamma_{x}\kappa_{x}-8G^{2})\varepsilon_{\rm L}+(8G^{2}{\rm e}^{-{\rm i}\theta}-4{\rm i}J\gamma_{x})\varepsilon_{\rm R}}{4J^{2}\gamma_{x}-8G^{2}\kappa _{x}+\kappa_{x}^{2}\gamma_{x}+16{\rm i}G^{2}J\cos\theta},\\ & {\text{δ}} c_{2+} = \dfrac{(8G^{2}{\rm e}^{{\rm i}\theta}-4{\rm i}J\gamma_{x})\varepsilon_{\rm L}+(2\gamma _{x}\kappa_{x}-8G^{2})\varepsilon_{\rm R}}{4J^{2}\gamma_{x}-8G^{2}\kappa _{x}+\kappa_{x}^{2}\gamma_{x}+16{\rm i}G^{2}J\cos\theta}, \end{split} $

其中$ \gamma_{x} = \gamma-2{\rm i}x $, $ \kappa_{x} = \kappa-2{\rm i}x $, $ \delta s_{-} = 0 $.
系统产生光学非互易性的物理根源是时间反演对称发生破缺, 这点也可从(5)式看出, 运动方程(5)式对应的系统等效哈密顿量为$ H_{\rm {eff}} = G{\text{δ}} c_{1}{\text{δ}} b+$$G{\rm e}^{-{\rm i}\theta}{\text{δ}} c_{2}{\text{δ}} b+J{\text{δ}} c_{1}{\text{δ}} c_{2}^{?}+{\rm {H.C.}} $. 当$ \theta\neq n{\text{π}} $($ n $为整数)时, 时间反演算符T与等效哈密顿算符$ H_{\rm {eff}} $不对易, 即$ [T,H]\neq0 $. 为研究系统的光学非互易性, 首先必须要求出系统左右两侧的输出光场$ \varepsilon^{\rm {out}}_{\rm L} $和 $ \varepsilon^{\rm {out}}_{\rm R} $. 输出场可由光力学中的输入输出关系^[57,58]得出, 即

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

其中$ \varepsilon_{\rm {L,R}}^{\rm {in}} = \varepsilon_{\rm {L,R}}/\sqrt{\kappa} $. 输出场也具有${\text{δ}} s = {\text{δ}} s_{+}{\rm e}^{-{\rm i}xt}+$$ {\text{δ}} s_{-}{\rm e}^{{\rm i}xt} $的形式, 由方程(7), (8)可得

$ \begin{split} & \varepsilon^{\rm {out}}_{\rm L+} = \sqrt{\kappa}{\text{δ}} c_{1+}-\varepsilon_{\rm L}/\sqrt{\kappa }\\ & \varepsilon^{\rm {out}}_{\rm R+} = \sqrt{\kappa}{\text{δ}} c_{2+}-\varepsilon_{\rm R}/\sqrt{\kappa} \end{split} $

和$ \varepsilon^{\rm {out}}_{\rm L-} = \varepsilon^{\rm {out}}_{\rm R-} = 0 $.

当系统呈现出完美的光学非互易性时, 传输振幅 $ T_{i\rightarrow j} $ ($ i,j =\rm{ L,R} $)应满足

$ \begin{align} T_{{\rm L}\rightarrow {\rm R}} = \left\vert \dfrac{\varepsilon_{\rm R}^{\rm {out}}}{\varepsilon _{\rm L}^{\rm {in}}}\right\vert _{\varepsilon_{\rm R}^{\rm {in}} = 0} = 1,\;T_{{\rm R}\rightarrow {\rm L}} = \left\vert \dfrac{\varepsilon_{\rm L}^{\rm {out}}}{\varepsilon_{\rm R}^{\rm {in}}}\right\vert _{\varepsilon_{\rm L}^{\rm {in}} = 0} = 0, \end{align} $

或者

$ \begin{align} T_{{\rm L}\rightarrow {\rm R}} = \left\vert \dfrac{\varepsilon_{\rm R}^{\rm {out}}}{\varepsilon _{\rm L}^{\rm {in}}}\right\vert _{\varepsilon_{\rm R}^{\rm {in}} = 0} = 0,\;T_{{\rm R}\rightarrow {\rm L}} = \left\vert \dfrac{\varepsilon_{\rm L}^{\rm {out}}}{\varepsilon_{\rm R}^{\rm {in}}}\right\vert _{\varepsilon_{\rm L}^{\rm {in}} = 0} = 1. \end{align} $

完美非互易性就意味着信号可以从系统的一侧完全传输到另一侧, 而另一侧的信号却一点也不可以传输过来. (10)式和(11)式代表光频隔离的两个不同的方向. 本文只讨论(10)式, 因为对(11)式的讨论是类似的. 下标$ \varepsilon_{\rm {R/L}}^{\rm {in}} = 0 $表示没有信号从右侧/左侧输入. 我们将忽略这些下标, 因为一般来说完美非互易性只讨论单侧输入的情况, 并且为方便将$ T_{i\rightarrow j} $简写为$ T_{ij} $.
由(7)式和(9)式可得输出场为

$ \begin{split} & \dfrac{\varepsilon_{\rm R+}^{\rm {out}}}{\varepsilon_{\rm L}^{\rm {in}}} = \dfrac{\kappa (8G^{2}{\rm e}^{{\rm i}\theta}-4{\rm i}J\gamma_{x})}{4J^{2}\gamma_{x}-8G^{2}\kappa_{x}+\kappa_{x}^{2}\gamma_{x}+16{\rm i}G^{2}J\cos\theta},\\ & \dfrac{\varepsilon_{\rm L+}^{\rm {out}}}{\varepsilon_{\rm R}^{\rm {in}}} = \dfrac{\kappa (8G^{2}{\rm e}^{-{\rm i}\theta}-4{\rm i}J\gamma_{x})}{4J^{2}\gamma_{x}-8G^{2}\kappa_{x}+\kappa_{x}^{2}\gamma_{x}+16{\rm i}G^{2}J\cos\theta}. \end{split} $

由(12)式可以看出, 当非互易相位差$ \theta = n{\text{π}} $(n为整数)时, 两输出场相等, 这说明光子传输是互易的. 而当$ \theta\neq n{\text{π}} $时, 两输出场不再相等, 即系统呈现出光学非互易性. 由(12)式可看出系统的非互易性来源于光力耦合相互作用G和腔模线性耦合相互作用J之间的量子相干效应. 由(12)式可得出(10)式成立时失谐x和线性耦合强度J必须满足的条件为:

$ \begin{split} & x = -\dfrac{\gamma\cot\theta}{2},\\ & J = -\dfrac{2G^{2}\sin\theta}{\gamma}. \end{split} $

首先讨论$ \theta = -\dfrac{{\text{π}}}{2} $的情况, 此时输出场(12)式变为:

$ \begin{split} & \dfrac{\varepsilon_{\rm R+}^{\rm {out}}}{\varepsilon_{\rm L}^{\rm {in}}} = \dfrac{-4{\rm i}\kappa (2G^{2}+J\gamma)}{4J^{2}\gamma-8G^{2}\kappa+\kappa^{2}\gamma},\\ & \dfrac{\varepsilon_{\rm L+}^{\rm {out}}}{\varepsilon_{\rm R}^{\rm {in}}} = \dfrac{-4{\rm i}\kappa (J\gamma-2G^{2})}{4J^{2}\gamma-8G^{2}\kappa+\kappa^{2}\gamma}. \end{split} $

把(13)式中的J代入(14)式中, 得到

$ \begin{split} & \dfrac{\varepsilon_{\rm R+}^{\rm {out}}}{\varepsilon_{\rm L}^{\rm {in}}} = \dfrac{-16{\rm i}G^{2}\kappa\gamma}{(4G^{2}-\kappa\gamma)^{2}},\\ & \dfrac{\varepsilon_{\rm L+}^{\rm {out}}}{\varepsilon_{\rm R}^{\rm {in}}} = 0. \end{split} $

由(10)式和(15)式, 可以得出系统出现完美非互易性时, 耦合强度G和J必须满足:

$ \begin{split} & G = G_{\pm} = \dfrac{1}{2}\sqrt{(3\pm2\sqrt{2})\kappa\gamma},\\ & J = J_{\pm} = \left(\dfrac{3}{2}\pm\sqrt{2}\right)\kappa. \end{split} $

由于系统处于蓝失谐驱动下, 在某些条件下系统会不稳定. 系统要稳定, 矩阵M (见(6)式)的本征值一定具有负实部. 根据劳斯-霍尔维茨(Routh-Hurwitz)稳态判据^[59], 可以得出具体的稳定条件如下:

$ \begin{align} 16G^{4}\kappa-4G^{2}\gamma^{2}(\gamma+\kappa)+\gamma^{2}\kappa(\gamma +\kappa)^{2} > 0. \end{align} $

系统所有的参数必须满足(17)式. 由(17)式可知, 当耦合$ G = G_{-} $时, 系统始终都是稳定的, 而当$ G = G_{+} $时, 系统参数只有满足以下条件时系统才是稳定的:

$ \gamma < \dfrac{1}{4}\left(\sqrt{2}-3+\sqrt{59+42\sqrt{2}}\right)\kappa\approx2.32\kappa. $

图2给出了在$ \theta = -\dfrac{{\text{π}}}{2} $, $ G = G_{+} $和$ J = J_{+} $(见(16)式)时, 在力学振子弛豫速率$ \gamma $不同取值下, 传输振幅$ T_{\rm {LR}} $(红线)和$ T_{\rm {RL}} $(黑线)随着标准化失谐$ x/\kappa $的变化曲线. 从图2可知, $ \theta = -\dfrac{{\text{π}}}{2} $时传输谱线左右对称、有多个峰值并且完美光学非互易性出现在共振频率位置($ x = 0 $). 通过计算可知, 在$ \gamma/\kappa\ll1 $时非互易传输谱线线宽为$ \Delta\omega\propto\gamma $, 就是为什么当$ \gamma $很小时传输谱线非常狭窄的原因(见图2(a)和图2(b)). 而当力学振子耗散速率$ \gamma $增大时, 谱线宽度也随之增加, 并且会出现传输增幅增益(幅值大于1)的现象(见图2(c)和图2(d)), 这是由于系统是处于蓝失谐驱动. 具体可由稳态条件(18)式看出, 当$ \gamma\gtrsim2.32\kappa $时, 系统会不稳定, 这就是为什么当$ \gamma $增大时系统会出现增益现象, 当$ \gamma/\kappa $增加到2时, 传输振幅$ T_{\rm {LR}} $和$ T_{\rm {RL}} $的最大值均已远远大于1(见图2(d)).
图 2 传输振幅$T_{\rm {LR}}$(红线)和$T_{\rm {RL}}$(黑线)在不同力学振子耗散速率下随着标准化失谐$x/\kappa$的变化曲线　(a) $\gamma/\kappa$=1/100; (b) $\gamma/\kappa$=1/10; (c) $\gamma/\kappa$=1; (d) $\gamma/\kappa$=2; 其他参数: $\theta=-\dfrac{{\text{π}}}{2}$, $G =G_{+}$和$J=J_{+}$(见(16)式)
Figure2. Transmission amplitudes $T_{\rm {LR}}$ (red line) and $T_{\rm {RL}}$ (black line) are plotted vs normalized detuning $x/\kappa$ for different cavity damping rate: (a) $\gamma/\kappa$=1/100; (b) $\gamma/\kappa$=1/10; (c) $\gamma/\kappa$=1; (d) $\gamma/\kappa$=2. Other parameters: $\theta=-\dfrac{{\text{π}}}{2}$, $G =G_{+}$ and $J=J_{+}$ according to Eq. (16)

图3给出了在$ \theta = -\dfrac{{\text{π}}}{2} $, $ G = G_{-} $和$ J = J_{-} $(见(16)式)时, 传输振幅$ T_{\rm {LR}} $(红线)和$ T_{\rm {RL}} $(黑线)在不同力学振子耗散速率下, 随着标准化失谐$ x/\kappa $的变化曲线. 从图3可知, 此时无论力学振子耗散速率取何值, 传输谱线$ T_{\rm {LR}} $都只有一个峰值, 并且峰值出现在共振频率的位置 ($ x = 0 $). 当$ \gamma/\kappa\ll1 $(见图3(a)和图3(b))时非互易传输谱线非常狭窄, 这是因为此时谱线线宽为$ \Delta\omega\propto\gamma $. 随着力学振子耗散速率的增加, 传输谱线的宽度也随之不断增加(见图3(c)), 同时传输谱线$ T_{\rm {RL}} $(黑线)的最大值则不断降低. 但现在即使力学振子耗散速率$ \gamma\gg\kappa $(见图3(d))时也不再出现增益现象(因为此时系统不会出现非稳情况).
图 3 传输振幅$T_{\rm {LR}}$(红线)和$T_{\rm {RL}}$(黑线)在不同力学振子耗散速率下随着标准化失谐$x/\kappa$的变化曲线　(a) $\gamma/\kappa$=1/100; (b) $\gamma/\kappa$=1/10; (c) $\gamma/\kappa$=1; (d) $\gamma/\kappa$=10. 其他参数: $\theta=-\dfrac{{\text{π}}}{2}$, $G =G_{-}$和$J=J_{-}$(见(16)式)
Figure3. Transmission amplitudes $T_{\rm {LR}}$ (red line) and $T_{\rm {RL}}$ (black line) are plotted vs normalized detuning $x/\kappa$ for different mechancial damping rate: (a) $\gamma/\kappa$=1/100; (b) $\gamma/\kappa$=1/10; (c) $\gamma/\kappa$=1; (d) $\gamma/\kappa$=10. Other parameters: $\theta=-\dfrac{{\text{π}}}{2}$, $G =G_{-}$ and $J=J_{-}$ according to Eq. (16)

由(13)式可知, 当$ \theta\neq\pm\dfrac{{\text{π}}}{2} $时, 非互易性不会出现在$ x\neq0 $的位置, 也就是说, 当$ \theta\neq\pm\dfrac{{\text{π}}}{2} $时, 非互易谱线一定是非对称的图形. 下面以$ \theta = -\dfrac{{\text{π}}}{4} $和$ -\dfrac{3{\text{π}}}{4} $为例讨论 $ \theta\neq\pm\dfrac{{\text{π}}}{2} $时传输谱线的性质. 由(12)式和(13)式可得, 当$ \theta = -\dfrac{{\text{π}}}{4} $和$ -\dfrac{3{\text{π}}}{4} $时, 若要实现完美的非互易性, 耦合强度G和J必须满足

$ \begin{split} & G = G_{\pm} = \dfrac{1}{2}\sqrt{\gamma\left[3\kappa-\gamma\pm\sqrt{(\kappa-\gamma )(\gamma+7\kappa)}\right]},\\ & J = J_{\pm} = \dfrac{\sqrt{2}}{4}\left[3\kappa-\gamma\pm\sqrt{(\kappa-\gamma )(\gamma+7\kappa)}\right]. \end{split} $

由(19)式可知, 当$ \theta = -\dfrac{{\text{π}}}{4} $和$ -\dfrac{3{\text{π}}}{4} $时, 要实现完美的非互易性, 则系统耗散必须满足$ \kappa\geqslant \gamma $. 图4给出了在力学振子耗散速率$ \gamma/\kappa = 10^{-3} $时, 传输振幅$ T_{\rm {LR}} $(红线)和$ T_{\rm {RL}} $(黑线)在不同非互易相位差θ和耦合强度G时 随着标准化失谐$ x/\gamma $的变化曲线. 其中耦合强度$ G_{\pm} $由(19)式给出. 由图4可知, 传输谱线$ T_{\rm {LR}} $(红线)和$ T_{\rm {RL}} $(黑线)在失谐$ x = \pm\dfrac{1}{2}\gamma $处确实会出现完美光学非互易性, 并且均会出现增益现象(增益是因为系统处于蓝失谐驱动). 虽然传输谱线$ T_{\rm {LR}} $和$ T_{\rm {RL}} $自身不是左右对称的, 但对于耦合$ G = G_{+} $与$ G = G_{-} $的谱线关于完美非互易位置轴对称(当$ \theta = -\dfrac{{\text{π}}}{4} $时, 对称轴$ x = \dfrac{1}{2}\gamma $, 如图4(a)和图4(b)所示; 当$ \theta = -\dfrac{3{\text{π}}}{4} $时, 对称轴$ x = -\dfrac{1}{2}\gamma $, 见图4(c)和图4(d)), 并且在整个非互易传输频谱上有$ T_{\rm {LR}} > T_{\rm {RL}} $.
图 4 传输振幅$T_{\rm {LR}}$(红线)和$T_{\rm {RL}}$(黑线)在不同非互易相位差θ和耦合强度G时随着标准化失谐$x/\gamma$的变化曲线　(a) $\theta=-\dfrac{{\text{π}}}{4}$和$G=G_{-}$; (b) $\theta=-\dfrac{{\text{π}}}{4}$和$G=G_{+}$; (c) $\theta=-\dfrac{3{\text{π}}}{4}$和$G=G_{-}$; (d) $\theta=-\dfrac{3{\text{π}}}{4}$和$G=G_{+}$; 其他参数: $\gamma/\kappa=10^{-3}$, $J=J_{\pm}$和$G=G_{\pm}$(见(19)式)
Figure4. Transmission amplitudes $T_{\rm {LR}}$ (red line) and $T_{\rm {RL}}$ (black line) are plotted vs normalized detuning $x/\gamma$ for different nonreciprocal phase θ and coupling strength G: (a) $\theta=-\dfrac{{\text{π}}}{4}$ and $G=G_{-}$; (b) $\theta=-\dfrac{{\text{π}}}{4}$ and $G=G_{+}$; (c) $\theta=-\dfrac{3{\text{π}}}{4}$ and $G=G_{-}$; (d) $\theta=-\dfrac{3{\text{π}}}{4}$ and $G=G_{+}$. Other parameters: $\gamma/\kappa=10^{-3}$, coupling strengths $J=J_{\pm}$ and $G=G_{\pm}$ according to Eq. (19)

本文研究了双腔光力学系统在蓝失谐驱动下的光学非互易性. 由系统中的光力耦合相互作用G和腔模线性耦合相互作用J之间的量子相干效应, 在某些条件下, 可以使系统呈现出完美的光学非互易现象. 首先研究了非互易相位差$ \theta =$$ -\dfrac{{\text{π}}}{2} $的情况, 研究发现, 当系统中各耗散速率(力学耗散速率γ和腔模耗散速率κ)一定的情况下, 会有两组耦合强度($ G = G_{\pm} $和$ J = J_{\pm} $)均可使系统出现完美非互易性. 由于系统处于蓝失谐驱动下, 会使系统出现非稳现象, 根据劳斯-霍尔维茨(Routh-Hurwitz)稳态判据我们给出了系统的稳定条件, 这种非稳现象也表现为非互易传输谱线会出现增益现象(谱线幅值大于1). 我们还发现当$ \gamma\ll\kappa $时, 非互易传输谱线的线宽$ \Delta\omega\propto\gamma $, 即当力学振子耗散速率很小时, 非互易传输谱线将会变得很狭窄. 最后研究了更一般的非互易相位差的情况$\left(\theta = -\dfrac{{\text{π}}}{4}\right.$和$\left. -\dfrac{3{\text{π}}}{4}\right)$, 并给出了实现完美非互易传输的必要条件. 这些研究结果有望能应用于光力系统中量子态转换、非互易传输等量子信息处理过程.