搜索

文章查询

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

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

谷开慧 严冬 张孟龙 殷景志 付长宝

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

谷开慧, 严冬, 张孟龙, 殷景志, 付长宝
PDF
HTML
导出引用
导出核心图
  • 随着纳米科技以及半导体技术的迅猛发展, 光力诱导透明、快慢光和光存储以及其他在光力系统中发现的量子光学和非线性光学效应成为人们目前研究的热点. 本文将薄膜腔光力系统同被束缚在腔中的二能级冷原子系综相耦合, 通过直接在薄膜振子上引入弱辅助驱动场来研究该原子辅助光力系统中原子和相位对量子相干性质及其快慢光的调控. 经过分析发现, 通过改变辅助驱动场的强度可直接实现对光力诱导透明窗口深度的调控, 通过改变辅助场与探测场之间的相位差, 可实现输出的探测场在“吸收”、“透明”和“增益”之间相互转换, 进而对弱探测场进行动态调控实现光开关. 与此同时, 还发现系统的群延迟时间随相位差的改变呈周期性变化. 通过调节相位差及原子数, 不但可以改变群延迟时间, 还可实现快慢光之间的相互转换.
      通信作者: 严冬, ydbest@126.com ; 殷景志, yjz886666@163.com
    • 基金项目: 教育部“春晖计划”项目(批准号: Z2017030)和吉林省自然科学基金(批准号: 2016286, GH16102)资助的课题.
    [1]

    Agarwal G S, Huang S M 2010 Phys. Rev. A 81 041803(R)

    [2]

    Safavi-Naeini A H, Mayer Alegre T P, Chan J, Eichenfield M, Winger M, Lin Q, Hill J T, Chang D E, Painter O 2011 Nature 472 69

    [3]

    Chang D E, Safavi-Naeini A H, Hafezi M, Painter O 2011 New J. Phys. 13 023003

    [4]

    Gu K H, Yan X B, Zhang Y, Fu C B, Liu Y M, Wang X, Wu J H 2015 Opt. Commun. 338 569

    [5]

    Chen B, Jiang C, Zhu K D 2011 Phys. Rev. A 83 055803

    [6]

    Jiang C, Liu H X, Cui Y S, Li X W, Chen G B, Chen B 2013 Opt. Express 21 12165

    [7]

    Tarhan D, Huang S, Müstecaplıoğlu Ö E 2013 Phys. Rev. A 87 013824

    [8]

    Fiore V, Yang Y, Kuzyk M C, Barbour R, Tian L, Wang H 2011 Phys. Rev. Lett. 107 133601

    [9]

    Tian L, Wang H L 2010 Phys. Rev. A 82 053806

    [10]

    Farman F, Bahrampour A R 2015 Phys. Rev. A 91 033828

    [11]

    McGee S A, Meiser D, Regal C A, Lehnert K W, Holland M J 2013 Phys. Rev. A 87 053818

    [12]

    Qu K, Agarwal G S 2013 Phys. Rev. A 87 031802(R)

    [13]

    Yan X B, Cui C L, Gu K H, Tian X D, Fu C B, Wu J H 2014 Opt. Express 22 4886

    [14]

    Barzanjeh S, Naderi M H, Soltanolkotabi M 2011 Phys. Rev. A 84 023803

    [15]

    Chan J, Alegr T P M, Safavi-Naeini A H, Hill J F, Krause A, Gröblacher S, Aspelmeyer M, Painter O 2011 Nature 478 89

    [16]

    Rabl P, Genes C, Hammerer K, Aspelmeyer M 2009 Phys. Rev. A 80 063819

    [17]

    Mari A, Eisert J 2009 Phys. Rev. Lett. 103 213603

    [18]

    Gu W J, Li G X 2013 Phys. Rev. A 88 013835

    [19]

    Mancini S, Giovannetti V, Vitali D, Tombesi P 2002 Phys. Rev. Lett. 88 120401

    [20]

    Hartmann M J, Plenio M B 2008 Phys. Rev. Lett. 101 200503

    [21]

    Mazzola L, Paternostro M 2011 Sci. Rep. 1 199

    [22]

    Yan X B 2017 Phys. Rev. A 96 053831

    [23]

    张秀龙, 鲍倩倩, 杨明珠, 田雪松 2018 物理学报 67 104203

    Zhang X L, Bao Q Q, Yang M Z, Tian X S 2018 Acta Phys. Sin. 67 104203

    [24]

    Stannigel K, Rabl P, Sorensen A S, Zoller P, Lukin M D 2010 Phys. Rev. Lett. 105 220501

    [25]

    Li J J, Zhu K D 2011 J. Appl. Phys. 110 114308

    [26]

    Agarwal G S, Huang S M 2014 New J. Phys. 16 033023

    [27]

    刘欢, 曹士英, 孟飞, 林百科, 方占军 2015 物理学报 64 094204

    Liu H, Cao S Y, Meng F, Lin B K, Fang Z J 2015 Acta Phys. Sin. 64 094204

    [28]

    Chang Y, Shi T, Liu Y X, Sun C P, Nori F 2011 Phys. Rev. A 83 063826

    [29]

    Fu C B, Yan X B, Gu K H, Cui C L, Wu J H, Fu T D 2013 Phys. Rev. A 87 053841

    [30]

    韩明, 谷开慧, 刘一谋, 张岩, 王晓畅, 田雪冬, 付长宝, 崔淬砺 2014 物理学报 63 094206

    Han M, Gu K H, Liu Y M, Zhang Y, Wang X C, Tian X D, Fu C B, Cui C L 2014 Acta Phys. Sin. 63 094206

    [31]

    Yan D, Wang Z H, Ren C N, Gao H, Li Y, Wu J H 2015 Phys. Rev. A 91 023813

    [32]

    He Q Y, Ficek Z 2014 Phys. Rev. A 89 022332

    [33]

    Genes C, Vitali D, Tombesi P 2008 Phys. Rev. A 77 050307

    [34]

    Genes C, Ritsch H, Drewsen M 2009 Phys. Rev. A 80 061803

    [35]

    Yi Z, Li G X, Wu S P, Yang Y P 2014 Opt. Express 22 20060

    [36]

    Walls D F, Milburn G J 1944 Quantum Optics (Berlin: Springer) p296

  • 图 1  由单模FP腔以及束缚在其中的N个全同二能级87Rb冷原子系综和中间由弱辅助驱动场${\varepsilon _{\rm{f}}}$驱动的振动频率为${\omega _{\rm{m}}}$的薄膜振子构成的光力学系统, 该系统的探测场和驱动场分别为${\varepsilon _{\rm{p}}}$${\varepsilon _{\rm{d}}}$

    Fig. 1.  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}}}$.

    图 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部分

    Fig. 2.  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$$\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部分

    Fig. 3.  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  $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部分

    Fig. 4.  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  ${\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

    Fig. 5.  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  $\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部分

    Fig. 6.  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.

  • [1]

    Agarwal G S, Huang S M 2010 Phys. Rev. A 81 041803(R)

    [2]

    Safavi-Naeini A H, Mayer Alegre T P, Chan J, Eichenfield M, Winger M, Lin Q, Hill J T, Chang D E, Painter O 2011 Nature 472 69

    [3]

    Chang D E, Safavi-Naeini A H, Hafezi M, Painter O 2011 New J. Phys. 13 023003

    [4]

    Gu K H, Yan X B, Zhang Y, Fu C B, Liu Y M, Wang X, Wu J H 2015 Opt. Commun. 338 569

    [5]

    Chen B, Jiang C, Zhu K D 2011 Phys. Rev. A 83 055803

    [6]

    Jiang C, Liu H X, Cui Y S, Li X W, Chen G B, Chen B 2013 Opt. Express 21 12165

    [7]

    Tarhan D, Huang S, Müstecaplıoğlu Ö E 2013 Phys. Rev. A 87 013824

    [8]

    Fiore V, Yang Y, Kuzyk M C, Barbour R, Tian L, Wang H 2011 Phys. Rev. Lett. 107 133601

    [9]

    Tian L, Wang H L 2010 Phys. Rev. A 82 053806

    [10]

    Farman F, Bahrampour A R 2015 Phys. Rev. A 91 033828

    [11]

    McGee S A, Meiser D, Regal C A, Lehnert K W, Holland M J 2013 Phys. Rev. A 87 053818

    [12]

    Qu K, Agarwal G S 2013 Phys. Rev. A 87 031802(R)

    [13]

    Yan X B, Cui C L, Gu K H, Tian X D, Fu C B, Wu J H 2014 Opt. Express 22 4886

    [14]

    Barzanjeh S, Naderi M H, Soltanolkotabi M 2011 Phys. Rev. A 84 023803

    [15]

    Chan J, Alegr T P M, Safavi-Naeini A H, Hill J F, Krause A, Gröblacher S, Aspelmeyer M, Painter O 2011 Nature 478 89

    [16]

    Rabl P, Genes C, Hammerer K, Aspelmeyer M 2009 Phys. Rev. A 80 063819

    [17]

    Mari A, Eisert J 2009 Phys. Rev. Lett. 103 213603

    [18]

    Gu W J, Li G X 2013 Phys. Rev. A 88 013835

    [19]

    Mancini S, Giovannetti V, Vitali D, Tombesi P 2002 Phys. Rev. Lett. 88 120401

    [20]

    Hartmann M J, Plenio M B 2008 Phys. Rev. Lett. 101 200503

    [21]

    Mazzola L, Paternostro M 2011 Sci. Rep. 1 199

    [22]

    Yan X B 2017 Phys. Rev. A 96 053831

    [23]

    张秀龙, 鲍倩倩, 杨明珠, 田雪松 2018 物理学报 67 104203

    Zhang X L, Bao Q Q, Yang M Z, Tian X S 2018 Acta Phys. Sin. 67 104203

    [24]

    Stannigel K, Rabl P, Sorensen A S, Zoller P, Lukin M D 2010 Phys. Rev. Lett. 105 220501

    [25]

    Li J J, Zhu K D 2011 J. Appl. Phys. 110 114308

    [26]

    Agarwal G S, Huang S M 2014 New J. Phys. 16 033023

    [27]

    刘欢, 曹士英, 孟飞, 林百科, 方占军 2015 物理学报 64 094204

    Liu H, Cao S Y, Meng F, Lin B K, Fang Z J 2015 Acta Phys. Sin. 64 094204

    [28]

    Chang Y, Shi T, Liu Y X, Sun C P, Nori F 2011 Phys. Rev. A 83 063826

    [29]

    Fu C B, Yan X B, Gu K H, Cui C L, Wu J H, Fu T D 2013 Phys. Rev. A 87 053841

    [30]

    韩明, 谷开慧, 刘一谋, 张岩, 王晓畅, 田雪冬, 付长宝, 崔淬砺 2014 物理学报 63 094206

    Han M, Gu K H, Liu Y M, Zhang Y, Wang X C, Tian X D, Fu C B, Cui C L 2014 Acta Phys. Sin. 63 094206

    [31]

    Yan D, Wang Z H, Ren C N, Gao H, Li Y, Wu J H 2015 Phys. Rev. A 91 023813

    [32]

    He Q Y, Ficek Z 2014 Phys. Rev. A 89 022332

    [33]

    Genes C, Vitali D, Tombesi P 2008 Phys. Rev. A 77 050307

    [34]

    Genes C, Ritsch H, Drewsen M 2009 Phys. Rev. A 80 061803

    [35]

    Yi Z, Li G X, Wu S P, Yang Y P 2014 Opt. Express 22 20060

    [36]

    Walls D F, Milburn G J 1944 Quantum Optics (Berlin: Springer) p296

  • [1] 严晓波, 杨柳, 田雪冬, 刘一谋, 张岩. 参量放大器腔中光力诱导透明与本征模劈裂性质. 物理学报, 2014, 63(20): 204201. doi: 10.7498/aps.63.204201
    [2] 王甫, 王智, 吴重庆, 刘国栋, 毛雅亚, 孙振超, 李强. 掺铒光纤中方波信号高次谐波的快慢光特性. 物理学报, 2015, 64(24): 244205. doi: 10.7498/aps.64.244205
    [3] 郑狄, 潘炜, 闫连山, 罗斌, 邹喜华, 江宁, 马雅男. 基于一种优化的梳状布里渊增益谱实现对任意周期信号的零展宽快慢光. 物理学报, 2010, 59(2): 1040-1046. doi: 10.7498/aps.59.1040
    [4] 赵建朋, 罗斌, 潘炜, 闫连山, 朱宏娜, 邹喜华, 叶佳. 光纤参量放大增益谱边带快慢光特性研究. 物理学报, 2014, 63(4): 044203. doi: 10.7498/aps.63.044203
    [5] 邓瑞婕, 闫智辉, 贾晓军. 基于电磁诱导透明机制的压缩光场量子存储. 物理学报, 2017, 66(7): 074201. doi: 10.7498/aps.66.074201
    [6] 付静, 刘万芳, 赵玉杰. 电磁诱导光透明过程中的Wigner-Yanse偏振信息. 物理学报, 2013, 62(17): 170302. doi: 10.7498/aps.62.170302
    [7] 刘岩, 张文明, 仲作阳, 彭志科, 孟光. 光梯度力驱动纳谐振器的非线性动力学特性研究. 物理学报, 2014, 63(2): 026201. doi: 10.7498/aps.63.026201
    [8] 林建潇, 吴九汇, 刘爱群, 陈喆, 雷浩. 光梯度力驱动的纳米硅基光开关. 物理学报, 2015, 64(15): 154209. doi: 10.7498/aps.64.154209
    [9] 张鹏飞, 李刚, 张玉驰, 杨榕灿, 郭龑强, 王军民, 张天才. 光致原子解吸附对冷原子磁光阱装载的动力学研究. 物理学报, 2010, 59(9): 6423-6429. doi: 10.7498/aps.59.6423
    [10] 李高芳, 马国宏, 马红, 初凤红, 崔昊杨, 刘伟景, 宋小军, 江友华, 黄志明, 褚君浩. 光抽运太赫兹探测技术研究ZnSe的光致载流子动力学特性. 物理学报, 2016, 65(24): 247201. doi: 10.7498/aps.65.247201
    [11] 张秀龙, 鲍倩倩, 杨明珠, 田雪松. 双腔光力学系统中输出光场纠缠特性的研究. 物理学报, 2018, 67(10): 104203. doi: 10.7498/aps.67.20172467
    [12] 李森, 李浩珍, 许静平, 朱成杰, 羊亚平. 基于腔光力学系统的全光三极管的压缩特性. 物理学报, 2019, 68(17): 174202. doi: 10.7498/aps.68.20190078
    [13] 朱琦, 潘佰良, 陈立, 王亚娟, 张迅懿. 光泵铯蒸气激光的动力学模型. 物理学报, 2010, 59(3): 1797-1801. doi: 10.7498/aps.59.1797
    [14] 陈雪, 刘晓威, 张可烨, 袁春华, 张卫平. 腔光力学系统中的量子测量. 物理学报, 2015, 64(16): 164211. doi: 10.7498/aps.64.164211
    [15] 任志君, 吴琼, 周卫东, 吴根柱, 施逸乐. 空间诱导产生艾里-贝塞尔光弹研究. 物理学报, 2012, 61(17): 174207. doi: 10.7498/aps.61.174207
    [16] 佘卫龙, 何穗荣, 汪河洲, 余振新, 莫党. 热自聚焦诱导光折变非对称自散焦. 物理学报, 1996, 45(12): 2022-2026. doi: 10.7498/aps.45.2022
    [17] 伍小平, 胡耿军, 李静, 龙潜, 陶陶, 张恭轩. 时域有限差分法数值仿真单光镊中微球受到的光阱力. 物理学报, 2011, 60(3): 030301. doi: 10.7498/aps.60.030301
    [18] 邓长发, 燕少安, 王冬, 彭金峰, 郑学军. 基于导电原子力显微镜的单根GaN纳米带光调控力电耦合性能. 物理学报, 2019, 68(23): 237304. doi: 10.7498/aps.68.20191097
    [19] 姚 鸣, 朱卡的, 袁晓忠, 蒋逸文, 吴卓杰. 声子辅助的电磁感应透明和超慢光效应的研究. 物理学报, 2006, 55(4): 1769-1773. doi: 10.7498/aps.55.1769
    [20] 张斌, 张浩佳, 杨秋红, 陆神洲. α-Al2O3透明陶瓷的发光及热释光特性. 物理学报, 2010, 59(2): 1333-1337. doi: 10.7498/aps.59.1333
  • 引用本文:
    Citation:
计量
  • 文章访问数:  164
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-07-25
  • 修回日期:  2018-10-23
  • 上网日期:  2019-03-23
  • 刊出日期:  2019-03-01

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

  • 1. 吉林大学电子科学与工程学院, 长春 130012
  • 2. 长春理工大学光电信息学院光电科学分院, 长春 130114
  • 3. 长春大学理学院, 材料设计与量子模拟实验室, 长春 130022
  • 4. 通化师范学院物理学院, 通化 134000
  • 通信作者: 严冬, ydbest@126.com ; 殷景志, yjz886666@163.com
    基金项目: 教育部“春晖计划”项目(批准号: Z2017030)和吉林省自然科学基金(批准号: 2016286, GH16102)资助的课题.

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

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}}}$

      Figure 1.  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部分

      Figure 2.  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部分

      Figure 3.  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部分

      Figure 4.  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

      Figure 5.  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部分

      Figure 6.  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窗口深度的调控. 当保持辅助场与探测场振幅比不变, 只改变它们之间的相位差时发现, 相位差的取值也直接影响着系统的光力耦合效应, 使得仅通过调节相位就可以实现输出的探测场在“吸收”、“透明”和“增益”之间相互转换, 从而改变整个系统的吸收和色散性质. 与此同时, 发现群延迟时间随相位差的改变呈周期性变化. 通过调节相位差及原子数, 可以改变群延迟时间, 使其表现出慢光或快光效应以及快慢光之间的转换. 利用该系统不但可以通过相位的调节来操控光在“吸收”和“透明”之间相互转换从而实现全光开关, 同时还可以通过光在“吸收”、“透明”和“增益”之间相互转换实现光电倍增, 因而该系统在量子光学器件方面有着潜在的应用价值.

参考文献 (36)

目录

    /

    返回文章
    返回