搜索

文章查询

x

留言板

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

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

Fabry-Perot腔与光学参量放大复合系统中实现可调谐的非常规光子阻塞

李宏 张斯淇 郭明 李美萱 宋立军

Fabry-Perot腔与光学参量放大复合系统中实现可调谐的非常规光子阻塞

李宏, 张斯淇, 郭明, 李美萱, 宋立军
PDF
HTML
导出引用
导出核心图
  • 本文提出在Fabry-Perot腔和光学参量放大复合系统中实现非常规光子阻塞效应. 此系统包含可调谐的复合型驱动强度相位, 用二阶关联函数描述光子统计性质, 数值模拟不同参数下的光子阻塞效应, 研究发现通过调节复合型驱动强度相位可以控制非常规光子阻塞. 在弱驱动条件下, 计算得到了强光子反聚束的最优化条件, 并给出了二阶关联函数解析式, 研究发现数值模拟结果与解析结果相符合. 研究结果为光子阻塞的相干操作提供了平台, 在量子信息处理和量子光学器件等方面具有潜在的应用前景.
      通信作者: 宋立军, songlj@jlenu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11347013)、吉林省教育厅科技规划项目(批准号: JJKH20190764KJ)、吉林工程技术师范学院博士科研启动经费专项(批准号: BSKJ201825)和吉林工程技术师范学院校级一般项目(批准号: XYB201820)资助的课题.
    [1]

    Imamoglu A, Schmidt H, Woods G, Deutsch M 1997 Phys. Rev. Lett. 79 1467

    [2]

    Liew T C H, Savona V 2010 Phys. Rev. Lett. 104 183601

    [3]

    Cao C, Mi S C, Wang T, Zhang R, Wang C 2016 IEEE J. Quantum Electron. 52 7000205

    [4]

    Cao C, Mi S C, Gao Y P, He L Y, Yang D, Wang T J, Zhang R, Wang C 2016 Sci. Rep. 6 22920

    [5]

    Cao Cong, Chen Xi, Duan Y W, Fan L, Zhang R, Wang T J, Wang C 2017 Optik 130 659

    [6]

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

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

    [7]

    廖庆洪, 叶杨, 李红珍, 周南润 2018 物理学报 67 40302

    Liao Q H, Ye Y, Li H Z, Zhou N R 2018 Acta Phys. Sin. 67 40302

    [8]

    Birnbaum K M, Boca A, Miller R, Boozer A D, Northup T E, Kimble H J 2005 Nature 436 87

    [9]

    Greentree A D, Tahan C, Cole J H, Hollenberg L C L 2006 Nat. Phys. 2 856

    [10]

    Angelakis D G, Santos M F, Bose S 2007 Phys. Rev. A 76 031805

    [11]

    Shen H Z, Zhou Y H, Yi X X 2015 Phys. Rev. A 91 063808

    [12]

    Shen H Z, Zhou Y H, Yi X X 2014 Phys. Rev. A 90 023849

    [13]

    Irvine W T M, Hennessy K, Bouwmeester D 2006 Phys. Rev. Lett. 96 057405

    [14]

    Zhou Y H, Shen H Z, Yi X X 2015 Phys. Rev. A 92 023838

    [15]

    Shen H Z, Zhou Y H, Liu H D, Wang G C, Yi X X 2015 Opt. Express 23 32835

    [16]

    Zhou Y H, Zhang S S, Shen H Z, Yi X X 2017 Opt. Lett. 42 1289

    [17]

    Shen H Z, Shang C, Zhou Y H, Yi X X 2018 Phys. Rev. A 98 023856

    [18]

    Shen H Z, Xu S, Zhou Y H, Wang G C, Yi X X 2018 J. Phys. B 51 035503

    [19]

    Zhou Y H, Shen H Z, Zhang X Y, Yi X X 2018 Phys. Rev. A 97 043819

    [20]

    Su S L, Tian Y Z, Shen H Z, Zang H P, Liang E J, Zhang S 2017 Phys. Rev. A 96 042335

    [21]

    Su S L, Gao Y, Liang E J, Zhang S 2017 Phys. Rev. A 95 022319

    [22]

    Su S L, Liang E J, Zhang S, Wen J J, Sun L L, Jin Z, Zhu A D 2016 Phys. Rev. A 93 012306

    [23]

    Zhou Y H, Shen H Z, Shao X Q, Yi X X 2016 Opt. Express 24 17332

    [24]

    Tang J, Geng W, Xu X 2015 Sci. Rep. 5 9252

    [25]

    Majumdar A, Bajcsy M, Rundquist A, Vuckovic J 2012 Phys. Rev. Lett. 108 183601

    [26]

    Zhang W, Yu Z, Liu Y, Peng Y 2014 Phys. Rev. A 8 043832

    [27]

    Flayac H, Savona V 2016 Phys. Rev. A 94 013815

    [28]

    Gerace D, Savona V 2014 Phys. Rev. A 89 031803

    [29]

    Lemonde M A, Didier N, Clerk A A 2014 Phys. Rev. A 90 063824

    [30]

    Xu X W, Li Y J 2013 J. Phys. B 46 035502

    [31]

    Wicz A, Li H R, Miranoao J Q, Nori F, Jing H 2018 Phys. Rev. Lett. 121 153601

    [32]

    石海泉, 谢智强, 徐勋卫, 刘念华 2018 物理学报 67 044203

    Shi H Q, Xie Z Q, Xu X W, Liu N H 2018 Acta Phys. Sin. 67 044203

    [33]

    Sarma B, Sarma A K 2017 Phys. Rev. A 96 053827

  • 图 1  (a) 用激光抽运OPA, 在腔内产生参量放大的腔结构示意图; (b) 量子干涉系统的跃迁路径

    Fig. 1.  (a) Schematic diagram of the cavity setup with an OPA which is pumped by a laser to produce parametric amplification in the cavity; (b) transition paths of the system for quantum interference.

    图 3  等时二阶关联函数$\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$数值结果随OPA非线性增益G和相位$\theta $的等高线图${\varOmega / {\kappa = 0.01}}$, ${\varDelta _a} = 1$, $\phi = 0.5\;{\rm{rad}}$

    Fig. 3.  Contour plot of the second-order correlation functions $\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$ vs. the nonlinear gain G of the OPA and phase $\theta $. Other parameters are ${\varOmega / {\kappa = 0.01}}$, ${\varDelta _a} = 1$, $\phi = 0.5\;{\rm{rad}}$.

    图 2  等时二阶关联函数${g^{\left(2\right)}}(0)$随OPA非线性增益G的变化${\varOmega / {\kappa = 0.01}}$, $\theta = - 0.0341{\text{π}}$, ${\varDelta _a} = 1$

    Fig. 2.  Variation curves of the zero-time-delay second-order correlation function ${g^{\left(2\right)}}(0)$ with the nonlinear gain G of the OPA. Other parameters are ${\varOmega / {\kappa = 0.01}}$, $\theta = - 0.0341{\text{π}}$, ${\varDelta _a} = 1$.

    图 4  等时二阶关联函数$\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$数值结果随相位$\theta $$\phi $的等高线图${\varOmega / {\kappa = 0.01}}$, ${\varDelta _a} = 1$, $G = {G_{\rm opt}}$

    Fig. 4.  Contour plot of the second-order correlation functions $\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$ vs. the phase $\theta $ and $\phi $. Other parameters are ${\varOmega / {\kappa = 0.01}}$, ${\varDelta _a} = 1$, $G = {G_{\rm opt}}$.

    图 5  等时二阶关联函数$\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$数值结果随OPA非线性增益$G$和相位$\phi $变化的等高线图${\varOmega / {\kappa = 0.01}}$, ${\varDelta _a} = 1$, $G = {G_{\rm opt}}$

    Fig. 5.  Contour plot of the second-order correlation functions $\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$ vs. the nonlinear gain of the optical parametric amplifier $G$ and phase $\phi $. Other parameters are ${\varOmega / {\kappa = 0.01}}$, ${\varDelta _a} = 1$, $G = {G_{\rm opt}}$.

    图 6  (a) 二阶关联函数${g^{\left(2\right)}}(0)$随OPA非线性增益$G$的变化, 其中蓝色实线由数值求解方程(3)得出, 红色菱形由(13)和(15)式解析得出; 其他参数为${\varDelta _a} = 1, \; \varOmega /\kappa = 0.01,$ $\phi = 0.5\;{\rm{rad}}, \; \theta = {\theta _{\rm opt}}$; (b) 二阶关联函数${g^{\left(2\right)}}(0)$随相位$\theta $的变化, 蓝色实线由数值求解方程(3)得出, 红色圆形由(13)和(15)式解析得出; 其他参数为${\varDelta _a} = 1, \;\varOmega /\kappa = 0.01,$ $\phi = 0.5\;{\rm{rad}}, \;G = {G_{\rm opt}}$

    Fig. 6.  (a) The second-order correlation functions ${g^{\left(2\right)}}(0)$ vs. the nonlinear gain of the optical parametric amplifier $G$; the blue solid line indicates the numerical results by numerically solving Eq. (3) and the red diamond corresponds to the analytical results of Eq. (13) and Eq. (15); other parameters are ${\varDelta _a} = 1, \; \varOmega /\kappa = 0.01, \; \phi = 0.5\;{\rm{rad}}, \; \theta = {\theta _{\rm opt}}$; (b) the second-order correlation functions ${g^{\left(2\right)}}(0)$ vs. the phase $\theta $; the blue solid line indicates the numerical results by numerically solving Eq. (3) and the red diamond corresponds to the analytical results of Eq. (13) and Eq. (15). Other parameters are ${\varDelta _a} = 1, \;\varOmega /\kappa = 0.01, \;\phi = 0.5\;{\rm{rad}},$ G = Gopt.

  • [1]

    Imamoglu A, Schmidt H, Woods G, Deutsch M 1997 Phys. Rev. Lett. 79 1467

    [2]

    Liew T C H, Savona V 2010 Phys. Rev. Lett. 104 183601

    [3]

    Cao C, Mi S C, Wang T, Zhang R, Wang C 2016 IEEE J. Quantum Electron. 52 7000205

    [4]

    Cao C, Mi S C, Gao Y P, He L Y, Yang D, Wang T J, Zhang R, Wang C 2016 Sci. Rep. 6 22920

    [5]

    Cao Cong, Chen Xi, Duan Y W, Fan L, Zhang R, Wang T J, Wang C 2017 Optik 130 659

    [6]

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

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

    [7]

    廖庆洪, 叶杨, 李红珍, 周南润 2018 物理学报 67 40302

    Liao Q H, Ye Y, Li H Z, Zhou N R 2018 Acta Phys. Sin. 67 40302

    [8]

    Birnbaum K M, Boca A, Miller R, Boozer A D, Northup T E, Kimble H J 2005 Nature 436 87

    [9]

    Greentree A D, Tahan C, Cole J H, Hollenberg L C L 2006 Nat. Phys. 2 856

    [10]

    Angelakis D G, Santos M F, Bose S 2007 Phys. Rev. A 76 031805

    [11]

    Shen H Z, Zhou Y H, Yi X X 2015 Phys. Rev. A 91 063808

    [12]

    Shen H Z, Zhou Y H, Yi X X 2014 Phys. Rev. A 90 023849

    [13]

    Irvine W T M, Hennessy K, Bouwmeester D 2006 Phys. Rev. Lett. 96 057405

    [14]

    Zhou Y H, Shen H Z, Yi X X 2015 Phys. Rev. A 92 023838

    [15]

    Shen H Z, Zhou Y H, Liu H D, Wang G C, Yi X X 2015 Opt. Express 23 32835

    [16]

    Zhou Y H, Zhang S S, Shen H Z, Yi X X 2017 Opt. Lett. 42 1289

    [17]

    Shen H Z, Shang C, Zhou Y H, Yi X X 2018 Phys. Rev. A 98 023856

    [18]

    Shen H Z, Xu S, Zhou Y H, Wang G C, Yi X X 2018 J. Phys. B 51 035503

    [19]

    Zhou Y H, Shen H Z, Zhang X Y, Yi X X 2018 Phys. Rev. A 97 043819

    [20]

    Su S L, Tian Y Z, Shen H Z, Zang H P, Liang E J, Zhang S 2017 Phys. Rev. A 96 042335

    [21]

    Su S L, Gao Y, Liang E J, Zhang S 2017 Phys. Rev. A 95 022319

    [22]

    Su S L, Liang E J, Zhang S, Wen J J, Sun L L, Jin Z, Zhu A D 2016 Phys. Rev. A 93 012306

    [23]

    Zhou Y H, Shen H Z, Shao X Q, Yi X X 2016 Opt. Express 24 17332

    [24]

    Tang J, Geng W, Xu X 2015 Sci. Rep. 5 9252

    [25]

    Majumdar A, Bajcsy M, Rundquist A, Vuckovic J 2012 Phys. Rev. Lett. 108 183601

    [26]

    Zhang W, Yu Z, Liu Y, Peng Y 2014 Phys. Rev. A 8 043832

    [27]

    Flayac H, Savona V 2016 Phys. Rev. A 94 013815

    [28]

    Gerace D, Savona V 2014 Phys. Rev. A 89 031803

    [29]

    Lemonde M A, Didier N, Clerk A A 2014 Phys. Rev. A 90 063824

    [30]

    Xu X W, Li Y J 2013 J. Phys. B 46 035502

    [31]

    Wicz A, Li H R, Miranoao J Q, Nori F, Jing H 2018 Phys. Rev. Lett. 121 153601

    [32]

    石海泉, 谢智强, 徐勋卫, 刘念华 2018 物理学报 67 044203

    Shi H Q, Xie Z Q, Xu X W, Liu N H 2018 Acta Phys. Sin. 67 044203

    [33]

    Sarma B, Sarma A K 2017 Phys. Rev. A 96 053827

  • [1] 王晶晶, 何博, 于波, 刘岩, 王晓波, 肖连团, 贾锁堂. 单光子调制锁定Fabry-Perot腔. 物理学报, 2012, 61(20): 204203. doi: 10.7498/aps.61.204203
    [2] 白岩, 赵卫疆, 任德明, 曲彦臣, 刘闯, 袁晋鹤, 钱黎明, 陈振雷. 基于有源Fabry-Perot腔的激光脉冲延时自外差研究. 物理学报, 2012, 61(9): 094218. doi: 10.7498/aps.61.094218
    [3] 程桂平, 柯莎莎, 李高翔, 张立辉. 光腔中两原子共振荧光的相干性质. 物理学报, 2007, 56(2): 830-836. doi: 10.7498/aps.56.830
    [4] 岳宏卫, 阎少林, 周铁戈, 谢清连, 游峰, 王争, 何明, 赵新杰, 方兰, 杨扬, 王福音, 陶薇薇. 嵌入Fabry-Perot谐振腔的高温超导双晶约瑟夫森结的毫米波辐照特性研究. 物理学报, 2010, 59(2): 1282-1287. doi: 10.7498/aps.59.1282
    [5] 王争, 赵新杰, 周铁戈, 岳宏卫, 阎少林, 何明. 嵌入到Fabry-Perot谐振腔的双晶约瑟夫森结阵列的阻抗匹配和相位锁定研究. 物理学报, 2010, 59(5): 3481-3487. doi: 10.7498/aps.59.3481
    [6] 刘明礼, 刘仁杰, 田相龙, 王亚伟, 雷海娜. Fabry-Perot腔谐振对横电波激励下亚波长一维金属光栅的异常透射性的作用. 物理学报, 2011, 60(2): 024217. doi: 10.7498/aps.60.024217
    [7] 詹达三. 完全相干场的二阶关联函数的分解性质. 物理学报, 1979, 164(1): 117-120. doi: 10.7498/aps.28.117
    [8] 甘琛利, 张彦鹏, 余孝军, 聂志强, 李 岭, 宋建平, 葛 浩, 姜 彤, 张相臣, 卢克清. 基于双光子不对称色锁二阶随机关联的阿秒极化拍研究. 物理学报, 2007, 56(5): 2670-2677. doi: 10.7498/aps.56.2670
    [9] 许 鸥, 鲁韶华, 简水生. 用于单频光纤激光器的光纤光栅双腔Fabry-Perot结构传输谱特性理论研究. 物理学报, 2008, 57(10): 6404-6411. doi: 10.7498/aps.57.6404
    [10] 张晨, 曹祥玉, 高军, 李思佳, 郑月军. 一种基于共享孔径Fabry-Perot谐振腔结构的宽带高增益磁电偶极子微带天线. 物理学报, 2016, 65(13): 134205. doi: 10.7498/aps.65.134205
    [11] 韩笑纯, 黄靖正, 方晨, 曾贵华. 群速度色散对于纠缠光场二阶关联函数影响的研究. 物理学报, 2015, 64(7): 070301. doi: 10.7498/aps.64.070301
    [12] 程金光. 高压调控的磁性量子临界点和非常规超导电性. 物理学报, 2017, 66(3): 037401. doi: 10.7498/aps.66.037401
    [13] 赵国栋, 杨亚利, 任伟. 钙钛矿型氧化物非常规铁电研究进展. 物理学报, 2018, 67(15): 157504. doi: 10.7498/aps.67.20180936
    [14] 周可余, 叶辉, 甄红宇, 尹伊, 沈伟东. 基于压电聚合物薄膜可调谐Fabry-Perot滤波器的研究. 物理学报, 2010, 59(1): 365-369. doi: 10.7498/aps.59.365
    [15] 李永贵, 张洪钧, 杨君慧, 高存秀. 混合型非线性Fabry-Perot标准具的光学双稳特性. 物理学报, 1982, 31(4): 446-459. doi: 10.7498/aps.31.446
    [16] 王建波, 钱进, 刘忠有, 陆祖良, 黄璐, 杨雁, 殷聪, 李同保. 计算电容中Fabry-Perot干涉仪测量位移的相位修正方法. 物理学报, 2016, 65(11): 110601. doi: 10.7498/aps.65.110601
    [17] 俞建宁, 张永祥, 孔贵芹. 振动筛系统的两类余维三分岔与非常规混沌演化. 物理学报, 2008, 57(10): 6182-6187. doi: 10.7498/aps.57.6182
    [18] 胡顺星, 魏合理, 谭 锟, 胡欢陵, 汪少林, 苏 嘉, 赵培涛, 曹开法. 基于三级Fabry-Perot标准具的纯转动拉曼测温激光雷达. 物理学报, 2008, 57(6): 3941-3947. doi: 10.7498/aps.57.3941
    [19] 钟权德. 内含等离子体的Fabry-Perot干涉仪的光学开关和双稳现象. 物理学报, 1985, 34(2): 182-187. doi: 10.7498/aps.34.182
    [20] 曾志文, 刘海涛, 张斯文. 基于Fabry-Perot模型设计亚波长金属狭缝阵列光学异常透射折射率传感器. 物理学报, 2012, 61(20): 200701. doi: 10.7498/aps.61.200701
  • 引用本文:
    Citation:
计量
  • 文章访问数:  157
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-01-27
  • 修回日期:  2019-04-04
  • 上网日期:  2019-08-16
  • 刊出日期:  2019-06-01

Fabry-Perot腔与光学参量放大复合系统中实现可调谐的非常规光子阻塞

  • 1. 吉林工程技术师范学院量子信息技术交叉学科研究院, 长春 130052
  • 2. 吉林省量子信息技术工程实验室, 长春 130052
  • 通信作者: 宋立军, songlj@jlenu.edu.cn
    基金项目: 国家自然科学基金(批准号: 11347013)、吉林省教育厅科技规划项目(批准号: JJKH20190764KJ)、吉林工程技术师范学院博士科研启动经费专项(批准号: BSKJ201825)和吉林工程技术师范学院校级一般项目(批准号: XYB201820)资助的课题.

摘要: 本文提出在Fabry-Perot腔和光学参量放大复合系统中实现非常规光子阻塞效应. 此系统包含可调谐的复合型驱动强度相位, 用二阶关联函数描述光子统计性质, 数值模拟不同参数下的光子阻塞效应, 研究发现通过调节复合型驱动强度相位可以控制非常规光子阻塞. 在弱驱动条件下, 计算得到了强光子反聚束的最优化条件, 并给出了二阶关联函数解析式, 研究发现数值模拟结果与解析结果相符合. 研究结果为光子阻塞的相干操作提供了平台, 在量子信息处理和量子光学器件等方面具有潜在的应用前景.

English Abstract

    • 单光子源在量子光学领域至关重要, 理想的单光子源在量子信息和量子通信领域意义重大. 1997年, Imamoglu等[1]首次提出光子阻塞的概念, 光子阻塞效应是实现单光子源的方法之一. 2010年, Liew和Savona[2]在弱Kerr非线性条件下仍获得强的反聚束效应, 首次提出非常规光子阻塞的概念. 近年来, 腔光力系统在理论和实验上都取得了快速的发展[3-7], 为量子现象和光学效应研究提供了新途径, 尤其在非线性方面为光子阻塞机制实现单光子源奠定了基础. Birnbaum等[8]首次在腔量子电动力学(quantum electro-dynamics, QED)系统中成功观测到光子阻塞效应以来, 人们分别在耦合空腔阵列[9,10]、一维光波导[11]、量子腔耦合系统[12]、光学系统[13-22]和电路QED系统[23]中观测到光子阻塞现象. 在实现非常规光子阻塞方面同样进展迅速, 如耦合单量子点腔系统[24]、有损双模纳米谐振腔系统[25,26]、三阶非线性耦合的单模腔[11,15,27]、二阶非线性耦合腔[10,28]、高斯压缩态[29]和耦合光机械系统[30]等均观测到了非常规光子阻塞现象. 2018年, Wicz等[31]利用非互易量子光学首次揭示了依赖方向的光子阻塞效应. 同年, 石海泉等[32]在多模光力系统中研究了非传统声子阻塞效应.

      本文研究基于包含简并光学参量放大(optical parametric amplifier, OPA)的Fabry-Perot腔内实现非常规光子阻塞的可能性. 尽管相关系统与内容已经有人研究, 但本文包含复合型驱动强度的相位, 并解析给出包含此相位的最优化条件. 通过数值模拟, 验证了优化条件的有效性. 研究发现复合型驱动的相位对非常规光子阻塞效应有着显著的影响.

    • 考虑典型的光力系统: 在Fabry-Perot腔中包含简并的OPA装置, 如图1(a)所示, 系统的哈密顿量可以写为[33]

      $\begin{split}{H_0} =\; & {\omega _a}{a^† }a + \varOmega ({a^† }{{\rm{e}}^{{\rm{i}}\phi + {\rm{i}}{\omega _{\rm{l}}}t}} + a{{\rm{e}}^{ - {\rm{i}}\phi - {\rm{i}}{\omega _{\rm{l}}}t}}) \\ &+ {\rm{i}}G({a^† }^2{{\rm{e}}^{{\rm{i}}\theta + 2{\rm{i}}{\omega _{\rm{l}}}t}} - {a^2}{{\rm{e}}^{ - {\rm{i}}\theta - 2{\rm{i}}{\omega _{\rm{l}}}t}}),\end{split}$

      其中$a({a^† })$是光模的降(升)算符, ${\omega _a}$是腔场的共振频率, ${\omega _{\rm{l}}}$是驱动激光频率, $\varOmega {{\rm{e}}^{{\rm{i}}\phi }}$是复合型驱动强度, $G$是OPA的非线性增益, $\theta $是外场驱动OPA的相位.

      假设驱动强度$\varOmega $非常小, 通过$U = \exp \left[ {{\rm{i}}{\omega _{\rm{l}}}t{a^† }a} \right]$旋转变换, 有效哈密顿量变为

      ${H_0} = {\varDelta _a}{a^† }a + \varOmega ({a^† }{{\rm{e}}^{{\rm{i}}\phi }} + a{{\rm e}^{ - {\rm{i}}\phi }}) + {\rm{i}}G({a^† }^2{{\rm{e}}^{{\rm{i}}\theta }} - {a^2}{{\rm{e}}^{ - {\rm{i}}\theta }}), $

      其中${\varDelta _a} = {\omega _a} - {\omega _{\rm{l}}}$表示腔的失谐量. 在下面的研究中, 主要探索相位$\phi $对光子阻塞的影响.

      图  1  (a) 用激光抽运OPA, 在腔内产生参量放大的腔结构示意图; (b) 量子干涉系统的跃迁路径

      Figure 1.  (a) Schematic diagram of the cavity setup with an OPA which is pumped by a laser to produce parametric amplification in the cavity; (b) transition paths of the system for quantum interference.

      用等时二阶关联函数分析光子阻塞效应, 其表达式如下:

      ${g^{\left( 2 \right)}}(0) = \frac{{\left\langle {{a^† }(t){a^† }(t)a(t)a(t)} \right\rangle }}{{{{\left\langle {{a^† }(t)a(t)} \right\rangle }^2}}}. $

      二阶关联函数${g^{\left( 2 \right)}}\left( 0 \right) > 1$表示光子存在聚束效应, 会极大提高腔内双光子存在的概率, 反之${g^{\left( 2 \right)}}\left( 0 \right) < 1$表示光子存在反聚束效应, 会有效抑制腔内双光子存在的概率. 如果二阶关联函数${g^{\left( 2 \right)}}\left( 0 \right) \to 0$, 表示系统处于完全光子阻塞机制下, 腔内同时出现两个光子的概率趋近于零.

      考虑到系统的耗散, 系统动力学演化过程可以由如下主方程描述:

      $\dot \rho = - {\rm{i}}[{H_0}, \rho ] + \frac{\kappa }{2}(2a\rho {a^† } - {a^† }a\rho - \rho {a^† }a), $

      其中$\kappa $表示腔的耗散率. 通过数值求解主方程得到稳态解

      ${g^{\left(2\right)}}(0) = \frac{{{\rm Tr}(\rho {a^† }{a^† }aa)}}{{{{[{\rm Tr}(\rho {a^† }a)]}^2}}}.$

    • 本文所有计算结果都是基于弱驱动条件, 令${\varOmega / {\kappa = 0.01}}$, 数值模拟结果如图2图3所示. 为方便起见, 将耗散率$\kappa $归一化. 在图2中, 展示了关联函数数值计算结果随相位$\phi $$G$的变化, 结果表明对于不同的$G$存在光子反聚束效应. 相位$\phi $分别取0.5, 0.8和1.2 rad, 数值结果表明$\phi = 0.5\;{\rm{rad}}$时对应优化的强反聚束效应. 为了得到对应反聚束效应的优化参数, 图3展示了$\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$随OPA非线性增益$G$和相位$\theta $的等高线图, 其他参数为$\phi = 0.5\;{\rm{rad}}$, ${\varDelta _a} = 1$. 由图3可见, 在一小区域内${g^{\left(2\right)}}(0)\ll 1$, 这区域内选定参数可实现强的反聚束效应.

      图  3  等时二阶关联函数$\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$数值结果随OPA非线性增益G和相位$\theta $的等高线图${\varOmega / {\kappa = 0.01}}$, ${\varDelta _a} = 1$, $\phi = 0.5\;{\rm{rad}}$

      Figure 3.  Contour plot of the second-order correlation functions $\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$ vs. the nonlinear gain G of the OPA and phase $\theta $. Other parameters are ${\varOmega / {\kappa = 0.01}}$, ${\varDelta _a} = 1$, $\phi = 0.5\;{\rm{rad}}$.

      图  2  等时二阶关联函数${g^{\left(2\right)}}(0)$随OPA非线性增益G的变化${\varOmega / {\kappa = 0.01}}$, $\theta = - 0.0341{\text{π}}$, ${\varDelta _a} = 1$

      Figure 2.  Variation curves of the zero-time-delay second-order correlation function ${g^{\left(2\right)}}(0)$ with the nonlinear gain G of the OPA. Other parameters are ${\varOmega / {\kappa = 0.01}}$, $\theta = - 0.0341{\text{π}}$, ${\varDelta _a} = 1$.

    • 系统演化过程用Fock态表示, 假设初始时刻处在$\left| 0 \right\rangle $态上, 系统含时演化的态$\left| \varPsi \right\rangle $

      $\left| \varPsi \right\rangle = {C_0}\left| 0 \right\rangle + {C_1}\left| 1 \right\rangle + {C_2}\left| 2 \right\rangle, $

      其中${C_m}\;(m = 0, 1, 2)$表示量子态的概率幅, 通过求解薛定谔方程可得到${C_m}$, 考虑到系统的耗散情况, 此时薛定谔方程为

      ${\rm{i}}\frac{{{\rm{d}}\left| \varPsi \right\rangle }}{{{\rm{d}}t}} = H'\left| \varPsi \right\rangle, $

      有效非厄米哈密顿量

      $H' = {\varDelta '_a}{a^† }a + \varOmega ({a^† }{{\rm{e}}^{{\rm{i}}\phi }} + a{{\rm{e}}^{ - {\rm{i}}\phi }}) + {\rm{i}}G({a^† }^2{{\rm{e}}^{{\rm{i}}\theta }} - {a^2}{{\rm{e}}^{ - {\rm{i}}\theta }}), $

      其中${\varDelta '_a} = {\varDelta _a} - {{{\rm{i}}\kappa } / 2}$. 把态$\left| \varPsi \right\rangle $和哈密顿量$H'$代入到(7)式中, 得到概率幅所满足的方程:

      $\begin{split} &{\rm{i}}{{\dot C}_0} = \varOmega {C_1}{{\rm{e}}^{ - {\rm{i}}\phi }} - \sqrt 2 {\rm{i}}G{\rm{e}}{}^{ - {\rm{i}}\theta }{C_2}, \\ &{\rm{i}}{{\dot C}_1} = \varOmega {C_0}{{\rm{e}}^{{\rm{i}}\phi }} + {{\varDelta '}_a}{C_1} + \sqrt 2 \varOmega {\rm{e}}{}^{ - {\rm{i}}\phi }{C_2}, \\ &{\rm{i}}{{\dot C}_2} = \sqrt 2 {\rm{i}}G{\rm{e}}{}^{{\rm{i}}\theta }{C_0} + \sqrt 2 \varOmega {\rm{e}}{}^{{\rm{i}}\phi }{C_1} + 2{{\varDelta '}_a}{C_2}. \end{split} $

      通过对系数耦合方程组的求解, 可以得到稳态解. 当态$\left| 2 \right\rangle $等于0时, 系统可达到完全光子阻塞效应, 在此条件下可以解出光子阻塞的最优化条件. 因此, 在(9)式中, 令${C_2} = 0$, 在弱驱动条件下, 方程(9)中第一个式子总是近似满足, 所以只需考虑最后两个方程进行计算:

      $\begin{split} & 0 = \varOmega {C_0}{{\rm e}^{{\rm{i}}\phi }} + {{\varDelta '}_a}{C_1}, \\ & 0 = {\rm{i}}G{{\rm{e}}^{i\theta }}{C_0} + \varOmega {{\rm{e}}^{{\rm{i}}\phi }}{C_1}. \end{split} $

      为了保证${C_0}$${C_1}$有非奇异解, 可解析给出优化条件

      ${G_{\rm opt}} = \frac{{2{\varOmega ^2}\cos 2\phi \sqrt {{{\tan }^2}2\phi + 1} }}{{\sqrt {4\varDelta _a^2 + {\kappa ^2}} }}$

      ${\theta _{\rm opt}} = \arctan \left[ {\frac{{\kappa \tan 2\phi - 2{\varDelta _a}}}{{\kappa + 2{\varDelta _a}\tan 2\phi }}} \right], $

      其中“opt”表示$G$$\theta $的优化解. 值得一提的是, 这些条件取决于腔失谐、驱动激光振幅和复合驱动强度的相位. 由于最优条件与驱动OPA的抽运场参数相对应, 所以这些参数可以通过调节OPA抽运场来控制. 当最优条件(11)和(12)式同时满足时, 可获得较强的反聚束效应.

      把方程(6)代入方程(3), 在弱耦合条件下, 态的概率幅满足${C_0} \gg {C_1} \gg C{}_2$, 此时得到等时二阶关联函数

      ${g^{\left(2\right)}}(0) \approx \frac{{2{{\left| {{C_2}} \right|}^2}}}{{{{\left| {{C_1}} \right|}^4}}}, $

      在弱抽运极限下, 基态布居数近似为1, 其他能级布居数微乎其微, 可忽略不计, 在这种情况下, 方程(9)变为

      $ \begin{split} & 0 = \varOmega {{\rm{e}}^{{\rm{i}}\phi }} + {{\varDelta '}_a}{C_1} + \sqrt 2 \varOmega {{\rm{e}}^{ - {\rm{i}}\phi }}{C_2}, \\ & 0 = \sqrt 2 {\rm{i}}\varOmega {{\rm{e}}^{i\theta }} + \sqrt 2 \varOmega {C_1}{{\rm{e}}^{{\rm{i}}\phi }} + 2{{\varDelta '}_a}{C_2}. \end{split} $

      求解(14)式得到

      $ \begin{split} & {C_1} = \frac{{\varOmega [{\rm{i}}G{{\rm{e}}^{{\rm{i}}\left( {\theta - \phi } \right)}} - {{\varDelta '}_a}{{\rm{e}}^{{\rm{i}}\phi }}]}}{{{{\varDelta '}_a}^2 - {\varOmega ^2}}}, \\ &{C_2} = - \frac{{{{\varDelta '}_a}{C_1} + \varOmega {{\rm{e}}^{{\rm{i}}\phi }}}}{{\sqrt 2 \varOmega {{\rm{e}}^{ - {\rm{i}}\phi }}}}. \end{split} $

      通过数值求解方程(3)模拟了$\lg\left[ {{g^{\left(2\right)}}\left(0\right)} \right]$随外场驱动OPA的相位$\theta $和复合型驱动相位$\phi $变化的等高线图, 结果如图4所示. 其中$G$满足(11)式, 即$G = {G_{\rm opt}}$. 其中红色虚线由(11)式画出, 研究发现非常规光子阻塞发生的地方正好是$G$取最优值的地方. 同样地, 在图5中, 模拟了$\lg\left[ {{g^{\left(2\right)}}\left(0\right)} \right]$随非线性增益$G$和复合型驱动相位$\phi $变化的等高线图, 其中$\theta $满足(12)式, 即$\theta = {\theta _{\rm opt}}$. 其中红色虚线部分由(12)式画出, 发现非常规光子阻塞发生的地方正好是$\theta $取最优值的地方. 由此可知, 光子的统计性质可以通过改变复合型驱动强度的相位$\phi $、OPA的非线性增益$G$和外场驱动OPA的相位$\theta $来调节. 系统的跃迁路径如图1(b)所示, 双光子激发共有2条跃迁路径, 一条是由驱动场激发从$\left| 0 \right\rangle $态到$\left| 1 \right\rangle $态, 然后从$\left| 1 \right\rangle $态再到$\left| 2 \right\rangle $态; 另一条是OPA作用直接从$\left| 0 \right\rangle $态到$\left| 2 \right\rangle $态. 当满足阻塞优化条件时, 这两条路径上的光子发生量子干涉相消, 干涉的结果为光子不能占据$\left| 2 \right\rangle $态, 所以发生强反聚束效应, 即产生光子阻塞效应.

      图  4  等时二阶关联函数$\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$数值结果随相位$\theta $$\phi $的等高线图${\varOmega / {\kappa = 0.01}}$, ${\varDelta _a} = 1$, $G = {G_{\rm opt}}$

      Figure 4.  Contour plot of the second-order correlation functions $\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$ vs. the phase $\theta $ and $\phi $. Other parameters are ${\varOmega / {\kappa = 0.01}}$, ${\varDelta _a} = 1$, $G = {G_{\rm opt}}$.

      图  5  等时二阶关联函数$\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$数值结果随OPA非线性增益$G$和相位$\phi $变化的等高线图${\varOmega / {\kappa = 0.01}}$, ${\varDelta _a} = 1$, $G = {G_{\rm opt}}$

      Figure 5.  Contour plot of the second-order correlation functions $\lg \left[ {{g^{\left(2\right)}}\left(0\right)} \right]$ vs. the nonlinear gain of the optical parametric amplifier $G$ and phase $\phi $. Other parameters are ${\varOmega / {\kappa = 0.01}}$, ${\varDelta _a} = 1$, $G = {G_{\rm opt}}$.

      图6中, 分别数值模拟了二阶关联函数${g^{\left(2\right)}}(0)$随非线性增益$G$和相位$\theta $的变化. 图6(a)为模拟得到的${g^{\left(2\right)}}(0)$$G$的变化曲线, 蓝色实线是由数值求解方程(3)得到的, 红色菱形是由(13)和(15)式得到的解析结果, 其中$\theta = {\theta _{\rm opt}}$$\phi = 0.5\;{\rm{rad}}$; 图6(b)为模拟得到的${g^{\left(2\right)}}(0)$$\theta $的变化曲线, 蓝色实线是由数值求解方程(3)得到的, 红色圆形是由(13)和(15)式得到的解析结果, 其中$G = {G_{\rm opt}}$$\phi = 0.5\;{\rm{rad}}$. 可以看出数值模拟与解析结果相符合, 说明了解析结果的正确性.

      图  6  (a) 二阶关联函数${g^{\left(2\right)}}(0)$随OPA非线性增益$G$的变化, 其中蓝色实线由数值求解方程(3)得出, 红色菱形由(13)和(15)式解析得出; 其他参数为${\varDelta _a} = 1, \; \varOmega /\kappa = 0.01,$ $\phi = 0.5\;{\rm{rad}}, \; \theta = {\theta _{\rm opt}}$; (b) 二阶关联函数${g^{\left(2\right)}}(0)$随相位$\theta $的变化, 蓝色实线由数值求解方程(3)得出, 红色圆形由(13)和(15)式解析得出; 其他参数为${\varDelta _a} = 1, \;\varOmega /\kappa = 0.01,$ $\phi = 0.5\;{\rm{rad}}, \;G = {G_{\rm opt}}$

      Figure 6.  (a) The second-order correlation functions ${g^{\left(2\right)}}(0)$ vs. the nonlinear gain of the optical parametric amplifier $G$; the blue solid line indicates the numerical results by numerically solving Eq. (3) and the red diamond corresponds to the analytical results of Eq. (13) and Eq. (15); other parameters are ${\varDelta _a} = 1, \; \varOmega /\kappa = 0.01, \; \phi = 0.5\;{\rm{rad}}, \; \theta = {\theta _{\rm opt}}$; (b) the second-order correlation functions ${g^{\left(2\right)}}(0)$ vs. the phase $\theta $; the blue solid line indicates the numerical results by numerically solving Eq. (3) and the red diamond corresponds to the analytical results of Eq. (13) and Eq. (15). Other parameters are ${\varDelta _a} = 1, \;\varOmega /\kappa = 0.01, \;\phi = 0.5\;{\rm{rad}},$ G = Gopt.

    • 本文研究了在Fabry-Perot腔和OPA复合系统中实现非常规光子阻塞效应, 给出了光子阻塞出现的最优化条件. 研究发现, 可以通过调整复合驱动强度中的相位来实现光子反聚束效应. 本文考虑了弱驱动和弱非线性条件, 通过选择最优解, 从数值和解析两方面论证了系统的强反聚束现象, 发现数值模拟与解析结果是一致的. 本研究为精确控制光子阻塞提供了方案, 并为制备优良单光子源提供了理论基础.

参考文献 (33)

目录

    /

    返回文章
    返回