删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

原子辅助光力系统中快慢光的量子调控

本站小编 Free考研考试/2021-12-29

摘要:随着纳米科技以及半导体技术的迅猛发展, 光力诱导透明、快慢光和光存储以及其他在光力系统中发现的量子光学和非线性光学效应成为人们目前研究的热点. 本文将薄膜腔光力系统同被束缚在腔中的二能级冷原子系综相耦合, 通过直接在薄膜振子上引入弱辅助驱动场来研究该原子辅助光力系统中原子和相位对量子相干性质及其快慢光的调控. 经过分析发现, 通过改变辅助驱动场的强度可直接实现对光力诱导透明窗口深度的调控, 通过改变辅助场与探测场之间的相位差, 可实现输出的探测场在“吸收”、“透明”和“增益”之间相互转换, 进而对弱探测场进行动态调控实现光开关. 与此同时, 还发现系统的群延迟时间随相位差的改变呈周期性变化. 通过调节相位差及原子数, 不但可以改变群延迟时间, 还可实现快慢光之间的相互转换.
关键词: 光力学/
光力诱导透明/
快慢光

English Abstract


--> --> -->
随着纳米技术和半导体工艺的迅猛发展, 人们越来越关注光力学这个新领域, 它为人们在高精密测量、量子信息处理以及量子基本原理的验证等方面打开了一扇大门. 近年来, 光力系统实现了一些经典的量子光学现象, 例如Agarwal和Huang[1]在Fabry-Perot (FP) 腔光力系统中发现了类似原子的电磁诱导透明现象(electromagnetically induced transparency, EIT), 即光力诱导透明现象(optomechanically induced transparency, OMIT), 同原子中的电磁诱导透明一样, 光力诱导透明效应不仅促进了量子光学的发展, 同时也在快慢光[2-7]、光存储[8-11]和光开关[12,13]等量子信息处理方面有着广泛的应用. 此外, 在光力系统中人们还实现了机械振子的基态冷却[14,15], 在宏观物体上观察到相干态[16]、压缩态[17,18]和纠缠态[19-23]等量子现象. 然而, 目前大部分典型腔光力学系统仍然只涉及两个自由度之间的耦合, 超过两个自由度耦合的多元光力学系统较少受到研究者关注.
最近在杂化的、多模的光力学系统中发现了一些奇异的量子现象, 例如量子态转移[24]、声子诱导透明[25]、相干完美吸收及相干完美透射[13,26,27]等, 这些研究成果在量子信息处理方面有着潜在的应用价值, 因此成为当前光力学的一个研究趋势. 同时在原子辅助光力系统中, Chang等[28]、Fu等[29]和韩明等[30]发现了电磁诱导透明的多稳现象, Yan等[31]则在与FP腔耦合的具有偶极-偶极相互作用的里德伯超级原子系综中发现了双稳现象, 另外还观察到了EPR (Einstein-Podolsky-Rosen)态的产生[32]以及基态冷却[33-35]等量子现象. 但是在原子辅助光力系统中通过相位调控快慢光的研究则较少见. 很明显这种杂化光力学系统中, 原子-腔-薄膜振子三模耦合, 原子作为辅助物质扮演着一个非常重要的角色, 它同光力系统的相互作用通过腔模联系起来, 控制原子可改变光力学腔中的腔场, 进而影响整个杂化光力系统的量子相干性, 从而实现对快慢光的量子调控, 因此原子辅助光力系统比常见的光力系统表现出了更多的量子操控自由度.
本文研究在弱辅助驱动场直接激发薄膜振子的原子辅助光力系统中, 原子和相位对OMIT及快慢光的调控. 首先讨论如何通过改变辅助场与探测场的相位差及振幅比来实现对OMIT的控制, 同时实现输出的探测光在“吸收”、“透明”和“增益”之间相互转换, 进而对弱探测光进行动态调控实现光开关. 其次讨论如何通过调节相位差及原子数来改变群延迟时间, 使其表现出慢光或快光效应以及快慢光之间的转换. 本文的研究有助于该系统应用于存储光脉冲、全关开关、光电倍增管等量子光机械器件的研制工作.
考虑如图1所示的原子辅助光力系统, 它由线性耦合的薄膜腔光力系统和被束缚在其中频率为${\omega _{\rm{a}}}$$N$个全同二能级87Rb超冷原子系综构成, 其中二能级的量子态$\left| 0 \right\rangle $$\left| 1 \right\rangle $分别对应87Rb原子的$D1$线的量子态$\left| {5{S_{1/2}}} \right\rangle $$\left| {5{P_{1/2}}} \right\rangle $. 选择的薄膜振子的本征频率为${\omega _{\rm{m}}}$, 弛豫速率为${\gamma _{\rm{m}}}$, 其力学品质因子为$Q=\omega_{\rm{m}} /\gamma_{\rm{m}}$. 如图1, 整个系统由外加的一个频率为${\omega _{\rm{d}}}$的强驱动场和一个频率为${\omega _{\rm{p}}}$的弱探测场驱动, 而薄膜振子则是直接由弱辅助驱动场${\varepsilon _{\rm{f}}}$驱动, 这三个外加相干场之间的频率假设满足${\omega _{\rm{p}}} - {\omega _{\rm{d}}} = {\omega _{\rm{f}}}$关系.
图 1 由单模FP腔以及束缚在其中的N个全同二能级87Rb冷原子系综和中间由弱辅助驱动场${\varepsilon _{\rm{f}}}$驱动的振动频率为${\omega _{\rm{m}}}$的薄膜振子构成的光力学系统, 该系统的探测场和驱动场分别为${\varepsilon _{\rm{p}}}$${\varepsilon _{\rm{d}}}$
Figure1. Schematic diagram of an optomechanical cavity containing N identical two-level cold 87Rb atoms with two fixed-end mirrors of equal reflectivity, which is driven by a strong coupling field ${\varepsilon _{\rm{d}}}$, a weak auxiliary drive field ${\varepsilon _{\rm{f}}}$and probed by a weak field ${\varepsilon _{\rm{p}}}$.

这里只考虑单模腔场的情况, 则系统总的哈密顿量为
$H = {H_0} + {H_{\rm{i}}} + {H_{\rm{d}}},$
其中第一项${H_0}$为系统的自由哈密顿量
${H_0} = \hbar {\omega _{\rm{c}}}{c^\dagger }c + \hbar {\omega _{\rm{a}}}{a^\dagger }a + \hbar {\omega _{\rm{m}}}{b^\dagger }b,$
式中右侧三项依次分别代表腔场的能量、原子的能量和薄膜振子的能量. 其中${\omega _{\rm{c}}}$为单模腔场的本征频率, $c$${c^\dagger }$分别表示腔中光子的湮灭和产生算符, 它们之间满足对易关系$[ c, {c^\dagger }] = 1$. ${\omega _{\rm{a}}}$表示原子在能级$\left| 0 \right\rangle $$\left| 1 \right\rangle $之间的跃迁频率, 在弱激发及原子数$N \gg 1$条件下, $a = {1 / {\sqrt N }}\Sigma _{i = 1}^N|0{\rangle _{ii}}\langle 1|$表示集体跃迁算符[28,29], 满足对易关系$[a, {a^\dagger }] = 1$. 由于中间的薄膜振子在其平衡位置做微小振动, 所以可用量子谐振子模型来描述其量子性质, 写成声子产生算符($b^\dagger$)和湮灭算符($b$)的形式, 满足对易关系$[b,b^\dagger]=1$. 第二项${H_{\rm{i}}}$为原子-腔和腔-薄膜振子相互作用哈密顿量
${H_{\rm{i}}} = \hbar {g_{{\rm{ac}}}}\sqrt N (c{a^\dagger } + {c^\dagger }a) - \hbar {g_{{\rm{mc}}}}{c^\dagger }c({b^\dagger } + b),$
式中右侧两项分别代表腔场同原子系综之间的相互作用和腔场与薄膜机械振子之间的相互作用. 其中${g_{{\rm{ac}}}} = \mu \sqrt {{\omega _{\rm{c}}}/(2\hbar V{\varepsilon _0})} $ 表示腔场与原子之间相互作用的耦合系数, $\mu $表示原子跃迁的电偶极矩, $V$表示腔的体积, ${\varepsilon _0}$表示真空介电常数. 光力耦合系数用${g_{{\rm{mc}}}}{{ = }}{{{\omega _{\rm{c}}}} / l}\sqrt {{\hbar / {(2m{\omega _{\rm{m}}})}}} $表示. (1)式中的第三项${H_{\rm{d}}}$为三个外加输入场与腔场和薄膜振子相互作用哈密顿量
$\begin{split} {H_{\rm{d}}} =\;& {\rm{i}}\hbar {\varepsilon _{\rm{d}}}({c^\dagger }{{\rm{e}}^{ - {\rm{i}}{\omega _{\rm{d}}}t}} - c{{\rm{e}}^{{\rm{i}}{\omega _{\rm{d}}}t}}) \\& + {\rm{i}}\hbar ({c^\dagger }{{\varepsilon}_{\rm{p}}}{{\rm{e}}^{ - {\rm{i}}{\omega _{\rm{p}}}t}} - c\varepsilon _{\rm{p}}^*{{\rm{e}}^{{\rm{i}}{\omega _{\rm{p}}}t}}) \\&+{\rm{ i}}\hbar ({b^\dagger }{\varepsilon _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}{\omega _{\rm{f}}}t}} - b\varepsilon _{\rm{f}}^{\rm{*}}{{\rm{e}}^{{\rm{i}}{\omega _{\rm{f}}}t}}), \end{split}$
式中右侧三项分别代表腔外输入的强驱动场和弱探测场与腔场之间的相互作用以及弱辅助驱动场与薄膜振子之间的相互作用. 其中${\varepsilon _{\rm{d}}} = \sqrt {2\kappa {P_{\rm{d}}}/(\hbar {\omega _{\rm{d}}})} $, $\left| {{\varepsilon _{\rm{p}}}} \right| = \sqrt {2\kappa {P_{\rm{p}}}/(\hbar {\omega _{\rm{p}}})} $$\left| {{\varepsilon _{\rm{f}}}} \right| = \sqrt {2\kappa {P_{\rm{f}}}/(\hbar {\omega _{\rm{f}}})} $分别代表驱动场、探测场和辅助驱动场的输入振幅, ${P_{\rm{d}}}$, ${P_{\rm{p}}}$${P_{\rm{f}}}$分别代表驱动场、探测场和辅助驱动场的输入功率.
在得出系统总的哈密顿量之后, 通过海森伯-郎之万方程研究体系的动力学性质,
$\dot a = - [{\gamma _{\rm{a}}} + {\rm{i}}({\omega _{\rm{a}}} - {\omega _{\rm{d}}})]a + {\rm{i}}{g_{{\rm{ac}}}}\sqrt N c + \sqrt {{\gamma _{\rm{a}}}} {a_{{\rm{in}}}},\tag{5a}$
$\dot b = - ({\gamma _{\rm{m}}} + {\rm{i}}{\omega _{\rm{m}}})b + {\rm{i}}{g_{{\rm{mc}}}}{c^\dagger }c + {\varepsilon _f}{{\rm e}^{ - {\rm{i}}{\omega _{\rm{f}}}t}} + \sqrt {{\gamma _{\rm{m}}}} {b_{{\rm{in}}}}, \tag{5b}$
$\begin{split} \dot c = & - [2\kappa + {\rm{i}}({\omega _{\rm{c}}} - {\omega _{\rm{d}}} - {g_{{\rm{mc}}}}(b + {b^\dagger }))]c \\& + {\varepsilon _{\rm{d}}} + {\varepsilon _{\rm{p}}}{{\rm{e}}^{ - {\rm{i}}\delta t}} + {\rm{i}}{g_{{\rm{ac}}}}\sqrt N a + \sqrt {2\kappa } c _{{\rm{in}}}, \end{split}\tag{5c}$
其中$\delta = {\omega _{\rm{p}}} - {\omega _{\rm{d}}} = {\omega _{\rm{f}}}$是探测场与驱动场之间的失谐, ${\gamma _{\rm{a}}}$为原子的弛豫速率, $\gamma_{\rm{m}}$为薄膜振子的弛豫速率, $\kappa $为腔场的弛豫速率. ${a_{{\rm{in}}}}$为作用在原子上的平均值为零的量子噪声, ${b_{{\rm{in}}}}$为作用在可移动薄膜振子上平均值为零的热噪声, ${c _{\rm{in}}}$为由输入场引入的平均值为零的量子真空噪声. 即
$\left\langle {{a_{{\rm{in}}}}} \right\rangle = \left\langle {{b_{{\rm{in}}}}} \right\rangle = \left\langle {{c_{{\rm{in}}}}} \right\rangle = 0.$
为了获得光力系统算符的平均值方程, 这里采用算符因式分解假设$\left\langle {{c^\dagger }c} \right\rangle = \left\langle {{c^\dagger }} \right\rangle \left\langle c \right\rangle $, $\left\langle {{b^\dagger }c} \right\rangle = \left\langle {{b^\dagger }} \right\rangle $$\left\langle c \right\rangle$$\left\langle {bc} \right\rangle = \left\langle b \right\rangle \left\langle c \right\rangle $, 则在稳态情况下, 算符在平均值附近涨落的关联与相应算符平均值乘积之比可忽略.
由此可得到系统算符的平均值方程为:
$\left\langle {\dot a} \right\rangle = - ({\gamma _{\rm{a}}} + {\rm{i}}{\varDelta _{\rm{a}}})\left\langle a \right\rangle + {\rm{i}}{g_{{\rm{ac}}}}\sqrt N \left\langle c \right\rangle ,\tag{7a}$
$\big\langle {\dot b} \big\rangle = - ({\gamma _{\rm{m}}} + {\rm{i}}{\omega _{\rm{m}}})\left\langle b \right\rangle + {\rm{i}}{g_{{\rm{mc}}}}\left\langle {{c^\dagger }} \right\rangle \left\langle c \right\rangle + {\varepsilon _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}\delta t}},\tag{7b}$
$\begin{split} \left\langle {\dot c} \right\rangle =\; & - \left[2\kappa + {\rm{i}}\left[{\varDelta _0} - {g_{{\rm{mc}}}}\left(\left\langle {{b^\dagger }} \right\rangle + \left\langle b \right\rangle \right)\right]\right]\left\langle c \right\rangle \\ &+ {\rm{i}}{g_{{\rm{ac}}}}\sqrt N \left\langle a \right\rangle + {\varepsilon _{\rm{d}}} + {\varepsilon _{\rm{p}}}{{\rm{e}}^{ - {\rm{i}}\delta t}}, \end{split} \tag{7c}$
其中${\varDelta _{\rm{a}}} = {\omega _{\rm{a}}} - {\omega _{\rm{d}}}$表示原子与驱动场之间的失谐, ${\varDelta _0} = {\omega _c} - {\omega _{\rm{d}}}$表示腔场与驱动场之间的失谐. 当不考虑含时间因子的探测场和辅助驱动场时, 令(7a)—(7c)式左侧为零, 可以得出系统相关算符的稳态平均值:
${a_{\rm{s}}} = \frac{{{\rm{i}}{g_{{\rm{ac}}}}\sqrt N }}{{{\gamma _{\rm{a}}} + {\rm{i}}{\varDelta _{\rm{a}}}}}{c_{\rm{s}}},\tag{8a}$
${b_{\rm{s}}} = \frac{{{\rm{i}}{g_{{\rm{mc}}}}|{c_{\rm{s}}}{|^2}}}{{{\gamma _{\rm{m}}} + {\rm{i}}{\omega _{\rm{m}}}}},\tag{8b}$
${c_{\rm{s}}} = \frac{{{\varepsilon _d}}}{{2\kappa + {\rm{i}}\varDelta + \dfrac{{g_{{\rm{ac}}}^{\rm{2}}N}}{{{\gamma _{\rm{a}}} + {\rm{i}}{\varDelta _{\rm{a}}}}}}},\tag{8c}$
其中
$\varDelta = {\omega _{\rm{c}}} - {\omega _{\rm{d}}} - {g_{{\rm{mc}}}}(b_{\rm{s}}^{\rm{*}} + {b_{\rm{s}}}) = {\varDelta _0} - {g_{{\rm{mc}}}}(b_{\rm{s}}^{\rm{*}} + {b_{\rm{s}}}),\tag{8d}$
表示腔场的有效失谐, 包含光压效应.
本文主要研究该原子辅助光力系统对探测场的线性响应, 所以在$\left| {{\varepsilon _{\rm{p}}}} \right| \ll \left| {{\varepsilon _{\rm{d}}}} \right|, \left| {{\varepsilon _{\rm{f}}}} \right| \ll \left| {{\varepsilon _{\rm{d}}}} \right|$情况下, 把系统算符写成各自平均值及其微小涨落的和:
$a = {a_{\rm{s}}} + \delta a, \quad {\rm{ }}b = {b_{\rm{s}}} + \delta b, \quad {\rm{ }}c = {c_{\rm{s}}} + \delta c.$
将(9)式代入(7)式中, 并且只保留线性项, 忽略非线性项$\delta {c^\dagger }\delta c$, $\delta b\delta c$, $\delta {b^\dagger }\delta c$和平均值为零的量子热噪声项, 可得涨落算符平均值满足的方程为:
$\langle \delta \dot a\rangle = - ({\gamma _{\rm{a}}} + {\rm{i}}{\varDelta _{\rm{a}}})\langle \delta a\rangle + {\rm{i}}{g_{{\rm{ac}}}}\sqrt N \langle \;\delta c\rangle ,\tag{10a}$
$\langle \delta \dot b\rangle = - ({\gamma _{\rm{m}}} + {\rm{i}}{\omega _{\rm{m}}})\langle \delta b\rangle + {\rm{i}}{g_{{\rm{mc}}}}(c_{\rm{s}}^{\rm{*}}\langle \delta c\rangle + {c_{\rm{s}}}\langle \delta {c^\dagger }\rangle ) + {\varepsilon _{\rm{f}}}{{\rm{e}}^{ - {\rm{i}}\delta t}}, \tag{10b}$
$\begin{split}\langle \delta \dot c\rangle =\; & - (2\kappa + {\rm{i}}\varDelta )\langle \delta c\rangle + {\rm{i}}{g_{{\rm{mc}}}}{c_{\rm{s}}}(\langle \delta b\rangle + \langle \delta {b^\dagger }\rangle ) \\ & + {\rm{i}}{g_{{\rm{ac}}}}\sqrt N \langle \delta a\rangle + {\varepsilon _{\rm{p}}}{{\rm{e}}^{ - {\rm{i}}\delta t}}.\end{split}\tag{10c}$
根据(10)式, 可以假设其解具有如下形式:
$\langle \delta s\rangle = \delta {s_ + }{\varepsilon _{\rm{p}}}{{\rm{e}}^{ - {\rm{i}}\delta t}} + \delta {s_ - }\varepsilon _{\rm{p}}^{\rm{*}}{{\rm{e}}^{{\rm{i}}\delta t}},$
其中$s = a, b, c$. 把(11)式代入(10)式, 比较等式两边含${\varepsilon _{\rm{p}}}{{\rm{e}}^{ - {\rm{i}}\delta t}}$$\varepsilon _{\rm{p}}^{\rm{*}}{{\rm{e}}^{ - {\rm{i}}\delta t}}$项的系数可得:

$ \delta {b_ + } = \frac{{ - g_{{\rm{mc}}}^{\rm{2}}|{c_{\rm{s}}}{|^2}({B_ - } - {B_ + })[{\rm{i}}{g_{{\rm{mc}}}}({R_ - } - {R_ + }){B_ - }{c_{\rm{s}}}^ * + {B_ - }{B_ + }{R_ - }Y{{\rm{e}}^{{\rm{i}}\varPhi }}]}}{{{B_ - }{B_ + }{R_ + }\left[ {{B_ - }{B_ + }{R_ - }{R_ + } + g_{{\rm{mc}}}^{\rm{2}}|{c_{\rm{s}}}{|^2}({B_ - } - {B_ + })({R_ - } - {R_ + })} \right]}} +\frac{{{\rm{i}}{g_{{\rm{mc}}}}{c_s}^ * + {B_ + }Y{{\rm{e}}^{ - {\rm{i}}\varPhi }}}}{{{B_ + }{R_ + }}},\tag{12a}$
$\delta {c_ + } = \frac{{{B_ - }{R_ - }{R_ + } - g_{{\rm{mc}}}^2|{c_{\rm{s}}}{|^2}({R_ - } - {R_ + })}}{{{B_ - }{B_ + }{R_ - }{R_ + } + g_{{\rm{mc}}}^{\rm{2}}|{c_s}{|^2}({B_ - } - {B_ + })({R_ - } - {R_ + })}} +\frac{{{\rm{i}}{g_{{\rm{mc}}}}{c_{\rm{s}}}{B_ - }{R_ - }Y{{\rm{e}}^{ - {\rm{i}}\varPhi }}}}{{{B_ - }{B_ + }{R_ - }{R_ + } + g_{{\rm{mc}}}^{\rm{2}}|{c_{\rm{s}}}{|^2}({B_ - } - {B_ + })({R_ - } - {R_ + })}}, \tag{12b}$

其中
$\begin{split}& {B_ + } = 2\kappa - {\rm{i}}\delta + {\rm{i}}\varDelta + \frac{{g_{{\rm{ac}}}^2N}}{{{\gamma _{\rm{a}}} - {\rm{i}}\delta + {\rm{i}}{\varDelta _{\rm{a}}}}}, \\& {B_ - } = 2\kappa - {\rm{i}}\delta - {\rm{i}}\varDelta + \frac{{g_{{\rm{ac}}}^{\rm{2}}N}}{{{\gamma _{\rm{a}}} - {\rm{i}}\delta - {\rm{i}}{\varDelta _{\rm{a}}}}}, \\& {R_ + } = {\gamma _{\rm{m}}} - {\rm{i}}\delta + {\rm{i}}{\omega _{\rm{m}}}{\rm{, }}\\& {R_ - } = {\gamma _{\rm{m}}} - {\rm{i}}\delta - {\rm{i}}{\omega _{\rm{m}}}, \\& Y = \Big|\frac{{{\varepsilon _{\rm{f}}}}}{{{\varepsilon _{\rm{p}}}}}\Big|,~~\varPhi = {\phi _{\rm{f}}} - {\phi _{\rm{p}}},\end{split}$
其中$Y$为辅助驱动场同探测场振幅大小之比, $\varPhi $为它们之间的相位差. 由(12b)式$\delta {c_ + }$的表达式可以看出, 第一项主要体现原子对光力学系统的影响, 第二项主要体现辅助驱动场对原子辅助光力学系统的影响.
相应地把输出的透射场${\varepsilon _{\rm{T}}}$和反射场${\varepsilon _{\rm{R}}}$的表达式也写成与(11)式相同的形式,
${\varepsilon _{{\rm{Tout}}}}(t) = {\varepsilon _{\rm{T}}}{\varepsilon _{\rm{p}}}{{\rm{e}}^{ - {\rm{i}}\delta t}} + {\varepsilon _{{\rm{T}} - }}\varepsilon _{\rm{p}}^ * {{\rm{e}}^{{\rm{i}}\delta t}},\tag{14a}$
${\varepsilon _{{\rm{Rout}}}}(t) = {\varepsilon _{\rm{R}}}{\varepsilon _{\rm{p}}}{{\rm{e}}^{ - {\rm{i}}\delta t}} + {\varepsilon _{{\rm{R}} - }}\varepsilon _{\rm{p}}^ * {{\rm{e}}^{{\rm{i}}\delta t}}\tag{14b}.$
这里需要指出的是, (14a)式和(14b)式中的${\varepsilon _{\rm{T}}}$${\varepsilon _{\rm{R}}}$表述的是输出场对探测场${\omega _{\rm{p}}}$的响应, 而${\varepsilon _{{\rm{T}} - }}$${\varepsilon _{{\rm{R}} - }}$描述的则是输出场对非线性四波混频过程产生的新频率$2{\omega _{\rm{d}}} - {\omega _{\rm{p}}}$的响应.
由腔的输入输出关系[36]
${\varepsilon _{{\rm{in}}}}(t) + {\varepsilon _{{\rm{out}}}}(t) = 2\kappa \left\langle {\delta c} \right\rangle , $
得到
$\begin{split}& {\varepsilon _{{\rm{Tout}}}}(t) = 2\kappa (\delta {c_ + }{\varepsilon _{\rm{p}}}{{\rm{e}}^{ - {\rm{i}}\delta t}} + \delta {c_ - }\varepsilon _{\rm{p}}^{\rm{*}}{{\rm{e}}^{{\rm{i}}\delta t}}), \\& {\varepsilon _{{\rm{Rout}}}}(t) + {\varepsilon _{\rm{p}}}{{\rm{e}}^{ - {\rm{i}}\delta t}} = 2\kappa (\delta {c_ + }{\varepsilon _{\rm{p}}}{{\rm{e}}^{ - {\rm{i}}\delta t}} + \delta {c_ - }\varepsilon _{\rm{p}}^{\rm{*}}{{\rm{e}}^{{\rm{i}}\delta t}}).\end{split}$
由此可得
$\begin{split}& {\varepsilon _{\rm{T}}} = 2\kappa \delta {c_ + }, ~~ {\varepsilon _{\rm{R}}} = 2\kappa \delta {c_ + } - 1,\end{split}$
这里只考虑${\varepsilon _{\rm{T}}} = 2\kappa \delta {c_ + }$${\varepsilon _{\rm{R}}} = 2\kappa {\rm{\delta}} {c_ + } - 1$是因为本文只关心输出场对探测场的响应. 这里${\varepsilon _{\rm{T}}}$的实部和虚部分别描述透射的探测场的吸收和色散性质, ${\varepsilon _{\rm{R}}}$的实部和虚部则分别描述反射的探测场的吸收和色散性质.
相应地, 根据(17)式可得到群延迟时间[7]为:
$\begin{split}& {\tau _{\rm T}} = \frac{{\partial {\phi _{\rm{T}}}}}{{\partial {\omega _{\rm{p}}}}} = {\rm{Im}}\left[ {\frac{1}{{{\varepsilon _{\rm{T}}}}}\frac{{\partial {\varepsilon _{\rm{T}}}}}{{\partial {\omega _{\rm{p}}}}}} \right],\\& {\tau _{\rm{R}}} = \frac{{\partial {\phi _{\rm{R}}}}}{{\partial {\omega _{\mathop{\rm p}\nolimits} }}} = {\rm{Im}}\left[ {\frac{1}{{{\varepsilon _{\rm{R}}}}}\frac{{\partial {\varepsilon _{\rm{R}}}}}{{\partial {\omega _{\rm{p}}}}}} \right].\end{split}$
群延迟时间如果大于零, 则系统表现出慢光现象, 群延迟时间如果小于零, 则为快光现象.
以下主要讨论原子辅助光力学腔中相位对OMIT及快慢光现象的影响. 本文采用的是实验上真实可行的实验参数[29,33], 并满足劳斯-胡尔维茨稳态判定准则[36], 确保系统稳态的稳定性. 各参数具体值为: $l \!=\! 1$ mm, $m \!=\! 10$ ng, ${\omega _{\rm{m}}}/(2{\text{π}}) \!=\! 10$ MHz, ${\gamma _{\rm{m}}}/(2{\text{π}}) = 100$ Hz, $\kappa = {\omega _{\rm{m}}}/10$, ${\varDelta _{\rm{a}}} = {\omega _{\rm{m}}}$, ${P_{\rm{d}}} = 1$ mW, $\lambda \! = \! 2{\text{π}}c/{\omega _{\rm{d}}} \! = \! 794.98$ nm, ${\gamma _{\rm{a}}}/2{\text{π}} \! = \! 2.875$ MHz, ${g_{{\rm{ac}}}}\sqrt N \!= $$3 \times 2{\text{π}}$ MHz.
图2所示为当考虑外加辅助场与探测场同相位时($\varPhi = 0$), 振幅比$Y$对杂化系统吸收和色散性质的影响. 如图2(a)所示, 在没有辅助驱动场的情况下($Y = 0$${\varepsilon _{\rm{f}}} = 0$时), 在共振$\delta = {\omega _{\rm{m}}}$处观察到了由原子-腔-薄膜振子三模耦合形成的杂化的EIT和OMIT. 如图2(b)所示, 当辅助场与机械振子耦合时($Y = 0.05$), 在共振$\delta = {\omega _{\rm{m}}}$处OMIT窗口的深度变浅, 这是因为外加辅助场对薄膜振子直接激发, 而辅助场的频率与薄膜振子频率共振${\omega _{\rm{f}}} = {\omega _{\rm{m}}} = \delta $, 使得薄膜振子在$\delta = {\omega _{\rm m}}$处发生共振吸收, 部分抵消了OMIT现象, 从而使得OMIT窗口深度变浅. 在图2(c)中发现, 随着$Y$的增加, OMIT窗口逐渐消失并形成了吸收峰, 这是因为随着辅助场强度的增加, 薄膜振子的共振吸收也随之增强, 从而导致光力效应被完全抵消掉, 产生吸收峰. 由以上的分析可知, 在辅助场与探测场同相位时, 改变辅助场同探测场振幅比$Y$可直接影响系统光力耦合效果, 从而实现对OMIT窗口深度的调控.
图 2 $\operatorname{Re} ({\varepsilon _{\rm{T}}})$(黑色实线)和${\rm{Im}}({\varepsilon _{\rm{T}}})$(红色虚线)随频率失谐$\delta /{\omega _{\rm{m}}}$的变化曲线 (a) $Y = 0$; $\varPhi = 0$; (b) $Y = 0.05$; $\varPhi = 0$; (c) $Y = 0.2$; $\varPhi = 0$; 其他参数值见文中第4部分
Figure2. The real (black line) and the imaginary (red dotted line) parts of ${\varepsilon _{\rm{T}}}$ as a function of $\delta /{\omega _{\rm{m}}}$with (a) $Y = 0$; $\varPhi = 0$; (b) $Y = 0.05$; $\varPhi = 0$; (c) $Y = 0.2$; $\varPhi = 0$. Relevant parameters are the same as those in Sec. 4.

图3所示为外加辅助场与探测场振幅比$Y = 0.2$时, 相位差$\varPhi $对系统吸收和色散性质的影响. 如图3(a)所示, 辅助场与探测场同相位, 即$\varPhi = 0$时, 在共振$\delta = {\omega _{\rm{m}}}$处OMIT消失形成了吸收峰. 但是随着相位差增加到$\varPhi = {\text{π}}/4$(图3(b)), 却发现在共振处出现了吸收为零的现象. 随着相位差进一步增加到$\varPhi = {\text{π}}/2$$\varPhi = 3{\text{π}}/4$时(如图3(c)图3(d)所示), 在共振$\delta = {\omega _{\rm{m}}}$处则出现了吸收增益现象. 随着相位差$\varPhi $的再增加(如图3(e)图3(h)所示), 也发现了由“增益”到“透明”再到“吸收”的过程, 图3(e)图3(h)的趋势与图3(a)图3(d)的趋势相同, 呈现一定的周期性变化. 实际上, 从(12a)式就可以看出, $\varPhi $的取值直接影响$\delta b_+$的大小, 也就是说辅助场与探测场的相位差直接影响着系统的光力耦合效应, 使得通过调节相位可以实现输出探测场在“吸收”、“透明”和“增益”之间相互转换, 从而改变整个系统的吸收和色散性质.
图 3 $Y = 0.2$$\operatorname{Re} ({\varepsilon _{\rm{T}}})$(黑色实线)和${\rm{Im}}({\varepsilon _{\rm{T}}})$(红色虚线)随频率失谐$\delta /{\omega _{\rm{m}}}$的变化曲线 (a) $\varPhi = 0$; (b) $\varPhi = {\text{π}}/4$; (c) $\varPhi = {\text{π}}/2$; (d) $\varPhi = 3{\text{π}}/4$; (e) $\varPhi = {\text{π}}$; (f) $\varPhi = 5{\text{π}} /4$; (g) $\varPhi = 6{\text{π}}/4$; (h) $\varPhi = 7{\text{π}}/4$; 其他参数值见文中第4部分
Figure3. The real (black line) and the imaginary (red dotted line) parts of ${\varepsilon _{\rm{T}}}$ as a function of $\delta /{\omega _{\rm{m}}}$with $Y = 0.2$: (a) $\varPhi = 0$; (b) $\varPhi = {\text{π}}/4$; (c) $\varPhi = {\text{π}}/2$; (d) $\varPhi = 3{\text{π}}/4$; (e) $\varPhi = {\text{π}} $; (f) $\varPhi = 5{\text{π}}/4$; (g) $\varPhi = 6{\text{π}}/4$; (h) $\varPhi = 7{\text{π}}/4$. Other parameters are the same as those in Sec. 4.

图4所示为辅助驱动场存在时, 原子数$N$对系统吸收和色散性质的影响. 在图4(a)图4(c)中, 令振幅比$Y = 0.2$, 相位差$\varPhi = {\text{π}}/4$, 可观察到随着原子数$N$的增加, 系统的吸收增益也随之增大. 相反, 在图4(e)图4(f)中, 令$Y = 0.2$, $\varPhi = 6{\text{π}}/4$, 可以观察到随着原子数$N$的增加, 系统的吸收增益反而逐渐减小. 也就是说在不同的相位下, 原子对系统吸收和色散性质的影响不同. 由此可知, 除了可通过相位$\varPhi $来操控系统的“吸收”、“透明”和“增益”以外, 还可通过改变原子数$N$来调控整个系统的吸收和色散性质.
图 4 $Y = 0.2$$\operatorname{Re} ({\varepsilon _{\rm{T}}})$(黑色实线)和${\rm{Im}}({\varepsilon _{\rm{T}}})$(红色虚线)随频率失谐$\delta /{\omega _{\rm{m}}}$的变化曲线 (a) $\varPhi = {\text{π}}/4$, $g\sqrt N = 3 \times 2{\text{π}}$ MHz; (b) $\varPhi = {\text{π}}/4$, $g\sqrt N = 6 \times 2{\text{π}}$ MHz; (c) $\varPhi = {\text{π}}/4$, $g\sqrt N = 9 \times 2{\text{π}}$ MHz; (d) $\varPhi = 6{\text{π}}/4$, $g\sqrt N = 3 \times 2{\text{π}}$ MHz; (e) $\varPhi = 6{\text{π}}/4$, $g\sqrt N = 6 \times 2{\text{π}}$ MHz; (f) $\varPhi = 6{\text{π}}/4$, $g\sqrt N = 9 \times 2{\text{π}}$ MHz; 其他参数值见文中第4部分
Figure4. The real (black line) and the imaginary (red dotted line) parts of ${\varepsilon _{\rm{T}}}$ as a function of $\delta /{\omega _{\rm m}}$with $Y = 0.2$: (a) $\varPhi = {\text{π}}/4$, $g\sqrt N = 3 \times 2{\text{π}}$ MHz; (b) $\varPhi = {\text{π}}/4$, $g\sqrt N = 6 \times 2{\text{π}}$ MHz; (c) $\varPhi = {\text{π}}/4$, $g\sqrt N = 9 \times 2{\text{π}}$ MHz; (d) $\varPhi = 6{\text{π}}/4$, $g\sqrt N = 3 \times 2{\text{π}}$ MHz; (e) $\varPhi = 6{\text{π}}/4$, $g\sqrt N = 6 \times {{2{\text{π}} }}$ MHz; (f) $\varPhi = 6{\text{π}}/4$, $g\sqrt N = 9 \times {{2{\text{π}} }}$ MHz. Other parameters are the same as those in Sec. 4.

图5所示为辅助驱动场与探测场相位差$\varPhi $对与色散性质紧密相关的群延迟时间的影响, 其中黑色实线对应透射的群延迟时间${\tau _{\rm{T}}}$, 红色虚线对应反射的群延迟时间${\tau _{\rm{R}}}$. 如图5(a)图5(h)所示, 透射群延迟时间${\tau _{\rm{T}}}$和反射群延迟时间${\tau _{\rm{R}}}$随着相位差$\varPhi $的增加, 出现了大于零及小于零的情况, 这意味着输出场包含着慢光和快光传播, 并伴有一定的周期性规律. 同时在$\varPhi = {\text{π}}/4$$\varPhi = 5{\text{π}}/4$时还观察到了快慢光同时存在的现象. 从前面的分析可知, 相位差$\varPhi $的取值影响整个系统的吸收和色散性质, 从而直接影响群延迟时间${\tau _{\rm{T}}}$${\tau _{\rm{R}}}$的取值. 因而, 通过调节相位差$\varPhi $可使原子辅助光力系统表现出慢光或快光效应, 并可实现快慢光之间的转换.
图 5 ${\tau _{\rm{T}}}$(黑色实线)和${\tau _{\rm{R}}}$(红色虚线)随频率失谐$\delta /{\omega _{\rm{m}}}$的变化曲线: (a) $\varPhi = 0$; (b) $\varPhi = {\text{π}}/4$; (c) $\varPhi = {\text{π}}/2$; (d) $\varPhi = 3{\text{π}}/4$; (e) $\varPhi = {\text{π}}$; (f) $\varPhi = 5{\text{π}} /4$; (g) $\varPhi = 6{\text{π}}/4$; (h) $\varPhi = 7{{{\text{π}} }}/4$; 其他参数取值同图3
Figure5. The ${\tau _{\rm{T}}}$ (black line) and the ${\tau _{\rm{R}}}$ (red dotted line) as a function of $\delta /{\omega _{\rm{m}}}$with (a) $\varPhi = 0$; (b) $\varPhi = {\text{π}}/4$; (c) $\varPhi = {\text{π}}/2$; (d) $\varPhi = 3{\text{π}}/4$; (e) $\varPhi = {\text{π}}$; (f) $\varPhi = 5{\text{π}}/4$; (g) $\varPhi = {{6{\text{π}} /}}4$; (h) $\varPhi = 7{\text{π}}/4$. Other parameters are the same as those in Fig.3.

图6(a)图6(c)中, 给出了在共振处$\delta = {\omega _{\rm{m}}}$, 随着原子数$N$的增加, 透射的群延迟时间${\tau _{\rm{T}}}$和反射的群延迟时间${\tau _{\rm{R}}}$随相位差$\varPhi /{\text{π}}$的变化情况. 可以看出, 黑色实线对应的${\tau _{\rm{T}}}$和红色虚线对应的${\tau _{\rm{R}}}$随相位差$\varPhi /{\text{π}}$呈周期性变化. 在图6(a)中反射和透射群延迟时间全部大于零, 表现出明显的慢光效应, 但随着原子数$N$的增加, 在图6(b)图6(c)中可以看出${\tau _{\rm{T}}}$${\tau _{\rm{R}}}$的取值逐渐减小, 其中${\tau _{\rm{R}}}$还出现了负值, 即表现出了快光现象. 出现这种情况的原因是原子数的增加, 改变了光力学腔中的腔场以及耦合强度, 进而影响了整个杂化光力系统的光学特性, 从而能够改变群延迟时间. 通过上述分析发现, 除了可以通过相位$\varPhi $对快慢光进行周期性调制外, 还可以通过控制原子数$N$来实现快慢光之间的转换.
图 6 $\delta = {\omega _{\rm{m}}}$${\tau _{\rm{T}}}$(黑色实线)和${\tau _{\rm{R}}}$(红色虚线)随相位差$\varPhi /{\text{π}}$和振幅比$Y$的变化曲线: (a) $Y = 0.2$, $g\sqrt N = 3 \times 2{\text{π}}$ MHz; (b) $Y = 0.2$, $g\sqrt N = 6 \times 2{\text{π}}$ MHz; (c) $Y = 0.2$, $g\sqrt N = 9 \times 2{\text{π}}$ MHz; 其他参数取值见文中第4部分
Figure6. The ${\tau _{\rm{T}}}$ (black line) and the ${\tau _{\rm{R}}}$ (red dotted line) as a function of $\varPhi /{\text{π}}$ and Y with $\delta = {\omega _{\rm{m}}}$: (a) $Y = 0.2$, $g\sqrt N = 3 \times 2{\text{π}}$ MHz; (b) $Y = 0.2$, $g\sqrt N = 6 \times 2{\text{π}}$ MHz; (c) $Y = 0.2$, $g\sqrt N = 9 \times 2{\text{π}}$ MHz. Other parameters are the same as those in Sec. IV.

本文主要研究了在薄膜振子上耦合弱辅助场的原子辅助光力系统中相位对OMIT及快慢光的调控. 当没有辅助场与薄膜振子耦合时, 在原子-腔-薄膜振子三模耦合的光力系统中观察到了杂化EIT和OMIT. 当辅助场直接激发薄膜振子且与探测场同相位时发现, 在共振处辅助场使薄膜振子产生共振吸收, 从而抑制了OMIT. 所以本文通过改变辅助场同探测场的相对振幅比, 直接影响了腔场与薄膜振子的光力耦合效果, 从而实现了对OMIT窗口深度的调控. 当保持辅助场与探测场振幅比不变, 只改变它们之间的相位差时发现, 相位差的取值也直接影响着系统的光力耦合效应, 使得仅通过调节相位就可以实现输出的探测场在“吸收”、“透明”和“增益”之间相互转换, 从而改变整个系统的吸收和色散性质. 与此同时, 发现群延迟时间随相位差的改变呈周期性变化. 通过调节相位差及原子数, 可以改变群延迟时间, 使其表现出慢光或快光效应以及快慢光之间的转换. 利用该系统不但可以通过相位的调节来操控光在“吸收”和“透明”之间相互转换从而实现全光开关, 同时还可以通过光在“吸收”、“透明”和“增益”之间相互转换实现光电倍增, 因而该系统在量子光学器件方面有着潜在的应用价值.
相关话题/系统 力学 光学 机械 控制

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 层状氧化钼的电子结构、磁和光学性质第一原理研究
    摘要:按照基于自旋密度泛函理论的赝势平面波第一原理计算方法,理论研究了两种层堆叠结构氧化钼(正交和单斜MoO3)的电子结构、磁性和光学特性,探讨其作为电致变色材料或电磁材料在光电子器件中的技术应用.采用先进的半局域GGA-PW91和非局域HSE06交换相关泛函精确计算晶体结构和带隙宽度.计算得出较低 ...
    本站小编 Free考研考试 2021-12-29
  • 库的量子关联相干辅助系统能量提取的研究
    摘要:基于单模微腔与二能级原子系综(库)构成的混合动力学模型,探索了非平衡库中量子关联相干(quantumcorrelatedcoherence,QCC)[TanKC,etal.2016Phys.Rev.A94,022329])对系统动力学的影响.推导了量子关联相干库下系统演化的动力学方程.借助于含 ...
    本站小编 Free考研考试 2021-12-29
  • 水中受激拉曼散射的能量增强及受激布里渊散射的光学抑制
    摘要:为提高液体介质中受激拉曼散射的输出能量,提出了通过温度调控来抑制受激布里渊散射的方法,设计了532nm多纵模宽带脉冲激光泵浦的受激拉曼散射发生系统,测量了不同温度下水中前向受激拉曼散射及后向受激布里渊散射的输出能量,分析了水温、泵浦激光线宽及热散焦效应对受激拉曼散射输出能量影响的物理机制.实验 ...
    本站小编 Free考研考试 2021-12-29
  • 动力学淬火过程中的不动点及衍生拓扑现象
    摘要:本文对近两年来有关淬火动力学过程中拓扑现象的研究做简要综述.这些动力学拓扑现象被动力学过程中的衍生拓扑不变量保护,与淬火前后体系的拓扑性质有密切关系.基于人工量子模拟平台的高度可控性,已在诸如超冷原子、超导量子比特、核磁共振、线性光学等众多物理体系中,通过对人工拓扑体系动力学过程的调控,观测到 ...
    本站小编 Free考研考试 2021-12-29
  • 光学微腔中一维费米气的磁性关联特性
    摘要:对于准一维两组分费米气与光学微腔耦合的系统,证明了微腔光子的超辐射可以驱动原子系统的磁性转变,该磁性转变与原子的失谐以及费米子的填充数密切相关.对于无相互作用原子气,在超辐射相区内平均场近似合理.基于该近似,分析了不同的填充和失谐情况下体系的静态自旋结构因子,由此刻画出腔光子协助的磁性关联转变 ...
    本站小编 Free考研考试 2021-12-29
  • 光晶格中超冷原子系统的磁激发
    摘要:囚禁在光学晶格中的旋量凝聚体由于其长的相干性和可调控性,使其成为时下热点的多比特量子计算的潜在候选载体,清楚地了解该体系的自旋和磁性的产生和调控就显得尤为重要.本文主要从理论上回顾了光晶格原子自旋链的磁性的由来和操控手段.从激光冷却原子出发,制备旋量玻色-爱因斯坦凝聚体,并装载进光晶格,最后实 ...
    本站小编 Free考研考试 2021-12-29
  • 超冷费米气体的膨胀动力学研究新进展
    摘要:多体系统的非平衡动力学演化是当前物理学中最具挑战性的问题之一.超冷量子费米原子气体具有较强的可控性,是研究多体非平衡动力学的理想系统,可以用来模拟和理解大爆炸后的早期宇宙、重离子碰撞中产生的夸克-胶子以及核物理等动力学.一般多体系统演化是非常复杂的,往往需要利用对称性来研究.利用Feshbac ...
    本站小编 Free考研考试 2021-12-29
  • 强场非次序双电离中再碰撞动力学的强度依赖
    摘要:利用三维经典系综模型系统地研究了不同强度线偏振激光脉冲驱动下He原子的非次序双电离.结果表明在非次序双电离中回碰电子的返回次数、两电子的碰撞距离和电子对的关联特性都强烈地依赖于激光强度.对于750nm,随着激光强度的增加,单次返回诱导的非次序双电离事件逐渐减少,而多次返回事件的比例显著增加.对 ...
    本站小编 Free考研考试 2021-12-29
  • 基于赝热光照明的单发光学散斑成像
    摘要:散射介质对光的散射是当前限制光学成像深度或距离的一个严重的问题.本文首先数值模拟比较了光透过随机散射介质成像研究中常用的基于光学记忆效应(memoryeffect,ME)和自相关(autocorrelation,AC)方法的HIO&ER算法和乒乓(Ping-Pang,PP)算法的优缺点.通过对 ...
    本站小编 Free考研考试 2021-12-29
  • 飞秒激光场中原子所受光学偶极力研究
    摘要:通过求解全波矢布洛赫方程研究了两能级原子与飞秒超快激光脉冲的相互作用过程,计算了不同拉比频率取值下原子所受光学偶极力和粒子数布居随时间的演化情况,分析了光场失谐量对光学势分布情况的影响.研究发现:由飞秒激光场产生的横向光力的时间平均值并不等于零,而是随着拉比频率的增加呈现振荡的增大趋势;纵向光 ...
    本站小编 Free考研考试 2021-12-29