搜索

x

留言板

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

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

基于四波混频过程和线性分束器产生四组份纠缠

余胜 刘焕章 刘胜帅 荆杰泰

基于四波混频过程和线性分束器产生四组份纠缠

余胜, 刘焕章, 刘胜帅, 荆杰泰
PDF
HTML
导出引用
导出核心图
  • 多组份纠缠是量子信息处理的重要资源, 它的产生通常涉及到许多复杂的线性和非线性过程. 本文从理论上提出了一种利用两个独立的四波混频过程和线性分束器产生真正的四组份纠缠的方案, 其中, 线性分束器的作用是将两个独立的四波混频过程联系起来. 首先应用部分转置正定判据研究了强度增益对四组份纠缠的影响, 结果表明, 在整个增益区域内都存在真正的四组份纠缠, 并且随着强度增益的增加, 纠缠也在增强. 然后研究了线性分束器的透射率对四组份纠缠的影响, 发现只要线性分束器的透射率不为0或1, 该系统也可以产生真正的四组份纠缠. 最后, 通过研究该系统可能存在的三组份纠缠和两组份纠缠来揭示该系统的纠缠结构. 本文理论结果为实验上利用原子系综四波混频过程产生真正的四组份纠缠提供了可靠的方案.
      通信作者: 刘胜帅, 852729711@qq.com ; 荆杰泰, jtjing@phy.ecnu.edu.cn
    • 基金项目: 国家级-基于人工表面等离激元的功能集成型辐射器研究(11874155, 91436211, 11374104, 10974057)
    [1]

    Greenberger D M, Horne M A, Zeilinger A 1993 Phys. Today 46 22

    [2]

    Lloyd S, Braunstein S 1999 Phys. Rev. Lett. 82 1784

    [3]

    Braunstein S L, van Loock P 2005 Rev. Mod. Phys. 77 513

    [4]

    Kimble H J 2008 Nature 453 1023

    [5]

    Bruβ D, Macchiavello C 2011 Phys. Rev. A 83 052313

    [6]

    Jing J T, Zhang J, Yan Y, Zhao F G, Xie C D, Peng K C 2003 Phys. Rev. Lett. 90 167903

    [7]

    Aoki T, Takei N, Yonezawa H, Wakui K, Hiraoka T, Furusawa A, van Loock P 2003 Phys. Rev. Lett. 91 080404

    [8]

    Yonezawa H, Aoki T, Furasawa A 2004 Nature 431 430

    [9]

    van Loock P, Braunstein S L 2000 Phys. Rev. Lett. 84 3482

    [10]

    Coelho A S, Barbosa F A S, Cassemiro K N, Villar A S, Martinelli M, Nussenzveig P 2009 Science 326 823

    [11]

    Armstrong S, Wang M, Teh R Y, Gong Q, He Q, Janousek J, Bachor H A, Reid M D, Lam P K 2015 Nat. Phys. 11 167

    [12]

    Cassermiro K N, Villar A S 2008 Phys. Rev. A 77 022311

    [13]

    Yokoyama S, Ukai R, Armstrong S C, Sornphiphatphong C, Kaji T, Suzuki S, Yoshikawa J I, Yonezawa H, Menicucci N C, Furusawa A 2013 Nat. Photonics 7 982

    [14]

    Roslund J, De Araujo R M, Jiang S, Fabre C, Treps N 2014 Nat. Photonics 8 109

    [15]

    Gerke S, Sperling J, Vogel W, Cai Y, Roslund J, Treps N, Fabre C 2015 Phys. Rev. Lett. 114 050501

    [16]

    Chen M, Menicucci N C, Pfister O 2014 Phys. Rev. Lett. 112 120505

    [17]

    Qin Z Z, Cao L M, Wang H L, Marino A M, Zhang W P, Jing J T 2014 Phys. Rev. Lett. 113 023602

    [18]

    Cao L M, Qi J, Du J J, Jing J T 2017 Phys. Rev. A 95 023803

    [19]

    Liu S S, Wang H L, Jing J T 2018 Phys. Rev. A 97 043846

    [20]

    Wang W, Cao L M, Lou Y B, Du J J, Jing J T 2018 Appl. Phys. Lett. 112 034101

    [21]

    Cao L M, Wang W, Lou Y B, Du J J, Jing J T 2018 Appl. Phys. Lett. 112 251102

    [22]

    Wang H L, Fabre C, Jing J T 2017 Phys. Rev. A 95 051802

    [23]

    Wang H L, Zheng Z, Wang Y X, Jing J T 2016 Opt. Express 24 23459

    [24]

    Lv S C, Jing J T 2017 Phys. Rev. A 96 043873

    [25]

    Lv S C, Jing J T 2018 Opt. Commun. 424 63

    [26]

    Abdisa G, Ahmed I, Wang X X, Liu Z C, Wang H X, Zhang Y P 2016 Phys. Rev. A 94 023849

    [27]

    Li C B, Jiang Z H, Zhang Y, Zhang Z Y, Wen F, Chen H X, Zhang Y P, Xiao M 2017 Phys. Rev. Appl. 7 014023

    [28]

    Zhang D, Li C B, Zhang Z Y, Zhang Y Q, Zhang Y P, Xiao M 2017 Phys. Rev. A 96 043847

    [29]

    Horodecki R, Horodecki P, Horodecki M, Horodecki K 2009 Rev. Mod. Phys. 81 865

    [30]

    Giedke G, Kraus B, Lewenstein M, Cirac J I 2001 Phys. Rev. Lett. 87 167904

    [31]

    Duan L M, Giedke G, Cirac J I, Zoller P 2000 Phys. Rev. Lett. 84 2722

    [32]

    Serafini A, Adesso G, Illuminati F 2005 Phys. Rev. A 71 032349

    [33]

    Peres A 1996 Phys. Rev. Lett. 77 1413

    [34]

    Simon R 2000 Phys. Rev. Lett. 84 2726

    [35]

    Werner R F, Wolf M M 2001 Phys. Rev. Lett. 86 3658

    [36]

    McCormick C F, Boyer V, Arimondo E, Lett P D 2007 Opt. Lett. 32 178

  • 图 1  产生四组份纠缠的简化图及铷-85 D1线的双Λ能级结构 (a) C0和C2是真空态注入, Pr0和Pr2是相干态注入; C1$ \rm Pr'_{1} $是第一个四波混频过程产生的孪生光束, C3$ \rm Pr'_{3} $是第二个四波混频过程产生的孪生光束; 光束Pr1和Pr3是光束$ \rm Pr'_{1} $$ \rm Pr'_{3} $经过线性分束器混合后产生的; (b)铷-85 D1线的双Λ能级结构, Δδ分别表示单光子失谐和双光子失谐

    Fig. 1.  A simplified diagram of quadripartite entanglement and an energy level diagram of rubidium-85: (a) C0 and C2 are vacuum states, Pr0 and Pr2 are coherent states; C1 and $ \rm Pr'_{1} $ are the twin beams generated by the first four-wave mixing process, C3 and $ \rm Pr'_{3} $ are the twin beams generated by the second four-wave mixing process; Pr1 and Pr3 are produced by mixing beams $ \rm Pr'_{1} $ and $ \rm Pr'_{3} $ through a linear beam splitter; (b) the double Λ energy level structure of D1 line in rubidium-85, Δ and δ represent one-photon detuning and two-photon detuning respectively.

    图 2  四种1 × 3情形的最小辛本征值v, 其为强度增益G1G2的函数 (a) C1被部分转置; (b) Pr1被部分转置; (c) Pr3被部分转置; (d) C3被部分转置

    Fig. 2.  The smallest symplectic eigenvalue v of all 1 × 3 scenarios, as a function of the power gains G1 and G2: (a) C1 is partially transposed; (b) Pr1 is partially transposed; (c) Pr3 is partially transposed; (d) C3 is partially transposed.

    图 3  三种2 × 2情形的最小辛本征值v, 其为强度增益G1G2的函数 (a) C1和Pr1被部分转置; (b) C1和Pr3被部分转置; (c) C1和C3被部分转置

    Fig. 3.  The smallest symplectic eigenvalues v of all 2 × 2 scenarios, as a function of the power gains G1 and G2: (a) C1 and Pr1 arepartially transposed; (b) C1 and Pr3 are partially transposed; (c) C1 and C3 are partially transposed.

    图 4  线性分束器的透射率η对四组份态的最小辛本征值v的影响 (a) C1被部分转置; (b) Pr1被部分转置; (c) Pr3被部分转置; (d) C3被部分转置; (e) C1和Pr1被部分转置; (f) C1和Pr3被部分转置; (g) C1和C3被部分转置

    Fig. 4.  Effect of the transmissivity of the linear beam splitter on the quadripartite entanglement of the system: (a) C1 is partially transposed; (b) Pr1 is partially transposed; (c) Pr3 is partially transposed; (d) C3 is partially transposed; (e) C1 and Pr1 are partially transposed; (f) C1 and Pr3 are partially transposed; (g) C1 and C3 are partially transposed.

    图 5  四个三组份态的最小辛本征值v, 其为强度增益G1G2的函数 (a)−(c)是由C1, Pr1和Pr3组成的三组份态的最小辛本征值v; (d)−(f)是由C1, Pr1和C3组成的三组份态的最小辛本征值v; (g)−(i)是由C1, Pr3和C3组成的三组份态的最小辛本征值v; (j)−(l)是由Pr1, Pr3和C3组成的三组份态的最小辛本征值v

    Fig. 5.  The smallest symplectic eigenvalues v of all tripartite states as a function of power gains G1 and G2: (a)−(c) The smallest symplectic eigenvalues v of tripartite state composed of C1, Pr1 and Pr3; (d)−(f) the smallest symplectic eigenvalues v of tripartite state composed of C1, Pr1 and C3; (g)−(i) the smallest symplectic eigenvalues v of tripartite state composed of C1, Pr3 and C3; (j)−(l) the smallest symplectic eigenvalues v of tripartite state composed of Pr1, Pr3 and C3.

    图 6  六种两组份态的最小辛本征值v, 其为强度增益G1G2的函数 (a) 由C1和Pr1组成的两组份态的最小辛本征值v; (b) 由C1和Pr3组成的两组份态的最小辛本征值v; (c) 由C1和C3组成的两组份态的最小辛本征值v; (d) 由Pr1和Pr3组成的两组份态的最小辛本征值v; (e) 由Pr1和C3组成的两组份态的最小辛本征值v; (f) 由Pr3和C3组成的两组份的最小辛本征值v

    Fig. 6.  The smallest symplectic eigenvalues v of all bipartite states as a function of power gains G1 and G2: (a) The smallest symplectic eigenvalues v of bipartite state composed of C1 and Pr1; (b) the smallest symplectic eigenvalues v of bipartite state composed of C1 and Pr3; (c) the smallest symplectic eigenvalues v of bipartite state composed of C1 and C3; (d) the smallest symplectic eigenvalues v of bipartite state composed of Pr1 and Pr3; (e) the smallest symplectic eigenvalues v of bipartite state composed of Pr1 and C3; (f) the smallest symplectic eigenvalues v of bipartite state composed of Pr3 and C3.

    表 1  四组份态的七种二分形式

    Table 1.  Seven partitions of quadripartite state.

    数目二分形式数目二分形式
    11|23422|134
    33|12444|123
    512|34613|24
    714|23
    下载: 导出CSV
  • [1]

    Greenberger D M, Horne M A, Zeilinger A 1993 Phys. Today 46 22

    [2]

    Lloyd S, Braunstein S 1999 Phys. Rev. Lett. 82 1784

    [3]

    Braunstein S L, van Loock P 2005 Rev. Mod. Phys. 77 513

    [4]

    Kimble H J 2008 Nature 453 1023

    [5]

    Bruβ D, Macchiavello C 2011 Phys. Rev. A 83 052313

    [6]

    Jing J T, Zhang J, Yan Y, Zhao F G, Xie C D, Peng K C 2003 Phys. Rev. Lett. 90 167903

    [7]

    Aoki T, Takei N, Yonezawa H, Wakui K, Hiraoka T, Furusawa A, van Loock P 2003 Phys. Rev. Lett. 91 080404

    [8]

    Yonezawa H, Aoki T, Furasawa A 2004 Nature 431 430

    [9]

    van Loock P, Braunstein S L 2000 Phys. Rev. Lett. 84 3482

    [10]

    Coelho A S, Barbosa F A S, Cassemiro K N, Villar A S, Martinelli M, Nussenzveig P 2009 Science 326 823

    [11]

    Armstrong S, Wang M, Teh R Y, Gong Q, He Q, Janousek J, Bachor H A, Reid M D, Lam P K 2015 Nat. Phys. 11 167

    [12]

    Cassermiro K N, Villar A S 2008 Phys. Rev. A 77 022311

    [13]

    Yokoyama S, Ukai R, Armstrong S C, Sornphiphatphong C, Kaji T, Suzuki S, Yoshikawa J I, Yonezawa H, Menicucci N C, Furusawa A 2013 Nat. Photonics 7 982

    [14]

    Roslund J, De Araujo R M, Jiang S, Fabre C, Treps N 2014 Nat. Photonics 8 109

    [15]

    Gerke S, Sperling J, Vogel W, Cai Y, Roslund J, Treps N, Fabre C 2015 Phys. Rev. Lett. 114 050501

    [16]

    Chen M, Menicucci N C, Pfister O 2014 Phys. Rev. Lett. 112 120505

    [17]

    Qin Z Z, Cao L M, Wang H L, Marino A M, Zhang W P, Jing J T 2014 Phys. Rev. Lett. 113 023602

    [18]

    Cao L M, Qi J, Du J J, Jing J T 2017 Phys. Rev. A 95 023803

    [19]

    Liu S S, Wang H L, Jing J T 2018 Phys. Rev. A 97 043846

    [20]

    Wang W, Cao L M, Lou Y B, Du J J, Jing J T 2018 Appl. Phys. Lett. 112 034101

    [21]

    Cao L M, Wang W, Lou Y B, Du J J, Jing J T 2018 Appl. Phys. Lett. 112 251102

    [22]

    Wang H L, Fabre C, Jing J T 2017 Phys. Rev. A 95 051802

    [23]

    Wang H L, Zheng Z, Wang Y X, Jing J T 2016 Opt. Express 24 23459

    [24]

    Lv S C, Jing J T 2017 Phys. Rev. A 96 043873

    [25]

    Lv S C, Jing J T 2018 Opt. Commun. 424 63

    [26]

    Abdisa G, Ahmed I, Wang X X, Liu Z C, Wang H X, Zhang Y P 2016 Phys. Rev. A 94 023849

    [27]

    Li C B, Jiang Z H, Zhang Y, Zhang Z Y, Wen F, Chen H X, Zhang Y P, Xiao M 2017 Phys. Rev. Appl. 7 014023

    [28]

    Zhang D, Li C B, Zhang Z Y, Zhang Y Q, Zhang Y P, Xiao M 2017 Phys. Rev. A 96 043847

    [29]

    Horodecki R, Horodecki P, Horodecki M, Horodecki K 2009 Rev. Mod. Phys. 81 865

    [30]

    Giedke G, Kraus B, Lewenstein M, Cirac J I 2001 Phys. Rev. Lett. 87 167904

    [31]

    Duan L M, Giedke G, Cirac J I, Zoller P 2000 Phys. Rev. Lett. 84 2722

    [32]

    Serafini A, Adesso G, Illuminati F 2005 Phys. Rev. A 71 032349

    [33]

    Peres A 1996 Phys. Rev. Lett. 77 1413

    [34]

    Simon R 2000 Phys. Rev. Lett. 84 2726

    [35]

    Werner R F, Wolf M M 2001 Phys. Rev. Lett. 86 3658

    [36]

    McCormick C F, Boyer V, Arimondo E, Lett P D 2007 Opt. Lett. 32 178

  • [1] 杨荣国, 张超霞, 李妮, 张静, 郜江瑞. 级联四波混频系统中纠缠增强的量子操控. 物理学报, 2019, 68(9): 094205. doi: 10.7498/aps.68.20181837
    [2] 杨 磊, 李小英, 王宝善. 利用光纤中自发四波混频产生纠缠光子的实验装置. 物理学报, 2008, 57(8): 4933-4940. doi: 10.7498/aps.57.4933
    [3] 仲银银, 潘晓州, 荆杰泰. 级联四波混频相干反馈控制系统量子纠缠特性. 物理学报, 2020, 69(13): 130301. doi: 10.7498/aps.69.20200042
    [4] 李述标, 武保剑, 文峰, 韩瑞. 高非线性光纤中四波混频的磁控机理研究. 物理学报, 2013, 62(2): 024213. doi: 10.7498/aps.62.024213
    [5] 黄德修, 张新亮, 李培丽. 基于PolSK调制的四波混频型超快全光译码器. 物理学报, 2009, 58(3): 1785-1792. doi: 10.7498/aps.58.1785
    [6] 曹亚敏, 武保剑, 万峰, 邱昆. 四波混频光相位运算器原理及其噪声性能研究. 物理学报, 2018, 67(9): 094208. doi: 10.7498/aps.67.20172638
    [7] 李培丽, 张新亮, 陈 俊, 黄黎蓉, 黄德修. 基于环行腔激光器四波混频型可调谐波长转换的理论研究. 物理学报, 2005, 54(3): 1222-1228. doi: 10.7498/aps.54.1222
    [8] 李建设, 李曙光, 赵原源, 刘强, 范振凯, 王光耀. 在单零色散微结构光纤中一次抽运同时发生两组四波混频的实验观察. 物理学报, 2016, 65(21): 214201. doi: 10.7498/aps.65.214201
    [9] 胡明列, 王清月, 栗岩峰, 王 专, 柴 路, 张伟力. 飞秒激光在双折射微结构光纤中模式控制的四波混频效应的实验研究. 物理学报, 2005, 54(9): 4411-4415. doi: 10.7498/aps.54.4411
    [10] 贾新鸿, 钟东洲, 陈海涛, 王 飞. 基于λ/4相移分布反馈半导体激光器四波混频的THz波长转换特性研究. 物理学报, 2007, 56(5): 2637-2646. doi: 10.7498/aps.56.2637
    [11] 苗向蕊, 高士明, 高 莹. 基于光纤四波混频效应的新型组播方法. 物理学报, 2008, 57(12): 7699-7704. doi: 10.7498/aps.57.7699
    [12] 刘瑞, 於亚飞, 张智明. 利用冷原子系综制备窄线宽三光子频率纠缠态. 物理学报, 2014, 63(14): 144203. doi: 10.7498/aps.63.144203
    [13] 孙 江, 左战春, 米 辛, 俞祖和, 吴令安, 傅盘铭. 引入量子干涉的双光子共振非简并四波混频. 物理学报, 2005, 54(1): 149-154. doi: 10.7498/aps.54.149
    [14] 孙 江, 郭庆林, 王英龙, 怀素芳, 王 颖, 左战春, 傅盘铭. 应用双光子共振非简并四波混频测量Ba原子里德伯态. 物理学报, 2006, 55(1): 221-225. doi: 10.7498/aps.55.221
    [15] 朱成禹, 吕志伟, 何伟明, 巴德欣, 王雨雷, 高 玮, 董永康. 布里渊增强四波混频时域特性的理论研究. 物理学报, 2007, 56(1): 229-235. doi: 10.7498/aps.56.229
    [16] 邓 莉, 孙真荣, 林位株, 文锦辉. 亚10 fs激光脉冲产生中的受激拉曼散射与四波混频效应. 物理学报, 2008, 57(12): 7668-7673. doi: 10.7498/aps.57.7668
    [17] 李博, 谭中伟, 张晓兴. 利用交叉相位调制和四波混频制作的时间透镜的仿真分析. 物理学报, 2012, 61(1): 014203. doi: 10.7498/aps.61.014203
    [18] 万峰, 武保剑, 曹亚敏, 王瑜浩, 文峰, 邱昆. 空频复用光纤中四波混频过程的解析分析方法. 物理学报, 2019, 68(11): 114207. doi: 10.7498/aps.68.20182129
    [19] 孙 江, 刘 霞, 牛金艳, 米 辛, 姜 谦, 吴令安, 傅盘铭. 布里渊增强非简并四波混频. 物理学报, 2008, 57(8): 4991-4994. doi: 10.7498/aps.57.4991
    [20] 王彦斌, 侯静, 陆启生, 彭杨, 陈子伦, 熊春乐. 长脉冲抽运光子晶体光纤四波混频和超连续谱的理论研究. 物理学报, 2011, 60(1): 014201. doi: 10.7498/aps.60.014201
  • 引用本文:
    Citation:
计量
  • 文章访问数:  595
  • PDF下载量:  46
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-01-06
  • 修回日期:  2020-02-24
  • 刊出日期:  2020-05-05

基于四波混频过程和线性分束器产生四组份纠缠

  • 1. 华东师范大学, 精密光谱科学与技术国家重点实验室, 上海 200062
  • 2. 山西大学, 极端光学协同创新中心, 太原 030006
  • 通信作者: 刘胜帅, 852729711@qq.com ; 荆杰泰, jtjing@phy.ecnu.edu.cn
    基金项目: 国家级-基于人工表面等离激元的功能集成型辐射器研究(11874155, 91436211, 11374104, 10974057)

摘要: 多组份纠缠是量子信息处理的重要资源, 它的产生通常涉及到许多复杂的线性和非线性过程. 本文从理论上提出了一种利用两个独立的四波混频过程和线性分束器产生真正的四组份纠缠的方案, 其中, 线性分束器的作用是将两个独立的四波混频过程联系起来. 首先应用部分转置正定判据研究了强度增益对四组份纠缠的影响, 结果表明, 在整个增益区域内都存在真正的四组份纠缠, 并且随着强度增益的增加, 纠缠也在增强. 然后研究了线性分束器的透射率对四组份纠缠的影响, 发现只要线性分束器的透射率不为0或1, 该系统也可以产生真正的四组份纠缠. 最后, 通过研究该系统可能存在的三组份纠缠和两组份纠缠来揭示该系统的纠缠结构. 本文理论结果为实验上利用原子系综四波混频过程产生真正的四组份纠缠提供了可靠的方案.

English Abstract

    • 多组份纠缠不仅在检验基本量子效应方面具有重要作用[1], 而且在量子通信以及量子计算领域也有着广泛的应用[2-5]. 在连续变量量子系统中, 已经有大量产生多组份纠缠的方案在理论上提出并在实验上实现[6-12]. 在这些方案中, 比较典型的方法是利用光学参量振荡器和线性分束器网络产生的独立单模压缩光束来生成连续变量多组份纠缠[9]. 此外, 另一种有效且简单的方法是在时域[13]和频域[14-16]中使用复用技术产生大规模连续变量簇态. 最近, 我们课题组已经在实验上证明了热铷原子蒸气池中的级联四波混频过程可以产生多个量子关联光束[17-22], 并在理论上提出了产生真正的三组份纠缠[23]和四组份纠缠[24,25]的方案, 这种级联是基于铷原子的双Λ能级结构. 不同于铷原子池的级联, 类原子系统能级级联的多波混频过程是基于类原子系统的N型能级结构, 这种级联不仅可以产生关联光束[26-28], 而且关联和压缩可以通过修饰场引起的相对非线性相移来控制[27]. 检验多组份量子纠缠时, 有很多判据可供选择[29-35]. 本文选择部分转置正定判据[34,35].

      由于强的非线性过程、空间多模性和可产生空间分离的非经典光束等优点, 热铷原子蒸气池中的四波混频过程是产生连续变量多组份纠缠的一项有前景的技术[36]. 本文提出了一个基于两个独立的四波混频过程和一个线性分束器来产生真正的四组份纠缠的方案, 其中, 线性分束器的作用是将两个独立的四波混频过程联系起来. 然后, 应用部分转置正定判据从理论上研究了系统的纠缠特性对强度增益和线性分束器的透射率的依赖性, 证实了该系统存在真正的四组份纠缠. 此外, 为了更好地解释系统的纠缠结构, 研究了该系统中可能存在的三组份纠缠和两组份纠缠.

    • 图1(a)为由线性分束器连接的两个独立的四波混频过程的结构图. 该方案可以产生四组份纠缠的关键是引入了线性分束器, 它将两个独立的四波混频过程联系起来. 图1(b)为铷-85 D1线的双Λ能级结构, 其中, 非线性相互作用强度极大地依赖于单光子失谐Δ和双光子失谐δ.

      图  1  产生四组份纠缠的简化图及铷-85 D1线的双Λ能级结构 (a) C0和C2是真空态注入, Pr0和Pr2是相干态注入; C1$ \rm Pr'_{1} $是第一个四波混频过程产生的孪生光束, C3$ \rm Pr'_{3} $是第二个四波混频过程产生的孪生光束; 光束Pr1和Pr3是光束$ \rm Pr'_{1} $$ \rm Pr'_{3} $经过线性分束器混合后产生的; (b)铷-85 D1线的双Λ能级结构, Δδ分别表示单光子失谐和双光子失谐

      Figure 1.  A simplified diagram of quadripartite entanglement and an energy level diagram of rubidium-85: (a) C0 and C2 are vacuum states, Pr0 and Pr2 are coherent states; C1 and $ \rm Pr'_{1} $ are the twin beams generated by the first four-wave mixing process, C3 and $ \rm Pr'_{3} $ are the twin beams generated by the second four-wave mixing process; Pr1 and Pr3 are produced by mixing beams $ \rm Pr'_{1} $ and $ \rm Pr'_{3} $ through a linear beam splitter; (b) the double Λ energy level structure of D1 line in rubidium-85, Δ and δ represent one-photon detuning and two-photon detuning respectively.

      在四波混频过程中, 湮灭的两个抽运光子会同时转换成一个探针光子和一个共轭光子. 在这里, 分别用$\hat a$$\hat b$标记探针光子和共轭光子, 相互作用强度用ε表示, $\hat H$表示相互作用哈密顿量. 在无耗尽及经典抽运近似条件下, 相互作用哈密顿量$\hat H$可以表示为

      $\hat H = {\rm{i}}\hbar \varepsilon {\hat b^\dagger }{\hat a^\dagger } + {\rm{h}}.{\rm{c}}.,$

      其中h.c.表示厄米共轭项.

      为了方便, 我们分别用1, 2, 3和4表示光束C1, Pr1, Pr3和C3. 由(1)式可得, 系统的输入输出关系可以写为

      $\begin{split} {{\hat a}_1} =\;& \sqrt {{G_1}} {{\hat a}_{{{\rm{C}}_0}}} + \sqrt {{g_1}} \hat a_{{\rm{P}}{{\rm{r}}_0}}^\dagger ,\\ {{\hat a}_2} =\;& \sqrt {\eta {G_1}} {{\hat a}_{{\rm{P}}{{\rm{r}}_0}}} + \sqrt {\eta {g_1}} \hat a_{{{\rm{C}}_0}}^\dagger \\ & + \sqrt {\left( {1 - \eta } \right){G_2}} {{\hat a}_{{\rm{P}}{{\rm{r}}_2}}} + \sqrt {\left( {1 - \eta } \right){g_2}} \hat a_{{{\rm{C}}_2}}^\dagger ,\\ {{\hat a}_3} =\;& - \sqrt {(1 - \eta ){G_1}} {{\hat a}_{{\rm{P}}{{\rm{r}}_0}}} - \sqrt {\left( {1 - \eta } \right){g_1}} \hat a_{{{\rm{C}}_0}}^\dagger \\ & + \sqrt {\eta {G_2}} {{\hat a}_{{\rm{P}}{{\rm{r}}_2}}} + \sqrt {\eta {g_2}} \hat a_{{{\rm{C}}_2}}^\dagger ,\\ {{\hat a}_4} =\;& \sqrt {{G_2}} {{\hat a}_{{{\rm{C}}_2}}} + \sqrt {{g_2}} \hat a_{{\rm{P}}{{\rm{r}}_2}}^\dagger , \end{split}$

      其中, ${\hat a_{{\rm{P}}{{\rm{r}}_{\rm{0}}}}}$${\hat a_{\rm P{r_2}}}$为相干态注入, ${\hat a_{{{\rm{C}}_{\rm{0}}}}}$${\hat a_{{{\rm{C}}_{\rm{2}}}}}$是真空态注入; G1G2分别是第一个四波混频过程和第二个四波混频过程的强度增益, 并且有Gkgk = 1 (k = 1, 2); η代表线性分束器的透射率. 为了方便分析, 假设光束$ \rm Pr'_{1} $$ \rm Pr'_{3} $的相位差ψ为0.

      在这里, 分别将正交振幅算符和正交相位算符定义为${\hat X_j} = (\hat a_j^\dagger + {\hat a_j}), {\hat Y_j} = {\rm{i}}(\hat a_j^\dagger - {\hat a_j})$, 其中$\hat a_j^\dagger $${\hat a_j}$是模式j的产生算符和湮灭算符. 正交振幅算符和正交相位算符的对易关系可以由产生算符和湮灭算符的对易关系推导得到, 它们的对易关系为$\left[ {{{\hat a}_j}, \hat a_j^\dagger } \right] = 1, \left[ {{{\hat X}_j}, {{\hat Y}_j}} \right] = 2{\rm{i}}$.

      基于(2)式和上述分析, 本方案中的正交振幅算符和正交相位算符的输入输出关系可以表示为

      $\left( \begin{aligned} {{\hat X}_1}\\ {{\hat X}_2}\\ {{\hat X}_3}\\ {{\hat X}_4} \end{aligned} \right) = \left( \begin{array}{ccccccc} \sqrt {{g_1}} &\sqrt {{G_1}} & 0 & 0\\ \sqrt {\eta {G_1}} & \sqrt {\eta {g_1}} & \sqrt {(1 - \eta ){G_2}} & \sqrt {(1 - \eta ){g_2}} \\ - \sqrt {(1 - \eta ){G_1}} & - \sqrt {(1 - \eta ){g_1}} & \sqrt {\eta {G_2}} & \sqrt {\eta {g_2}} \\ 0 & 0 & \sqrt {{g_2}} & \sqrt {{G_2}} \end{array} \right)\;\left( \begin{aligned} {{\hat X}_{{\rm{P}}{{\rm{r}}_{\rm{0}}}}}\\ {{\hat X}_{{{\rm{C}}_{\rm{0}}}}}\\ {{\hat X}_{{\rm{P}}{{\rm{r}}_{\rm{2}}}}}\\ {{\hat X}_{{{\rm{C}}_{\rm{2}}}}} \end{aligned} \right)$

      $\left( \begin{aligned} {{\hat Y}_1}\\ {{\hat Y}_2}\\ {{\hat Y}_3}\\ {{\hat Y}_4} \end{aligned} \right) = \left( \begin{array}{ccccccc} - \sqrt {{g_1}} &\sqrt {{G_1}} & 0& 0\\ \sqrt {\eta {G_1}} &- \sqrt {\eta {g_1}}&\sqrt {(1 - \eta ){G_2}}&- \sqrt {(1 - \eta ){g_2}} \\ - \sqrt {(1 - \eta ){G_1}}&\sqrt {(1 - \eta ){g_1}}&\sqrt {\eta {G_2}} & - \sqrt {\eta {g_2}} \\ 0&0& - \sqrt {{g_2}} &\sqrt {{G_2}} \end{array} \right)\;\left( \begin{aligned} {{\hat Y}_{{\rm{P}}{{\rm{r}}_{\rm{0}}}}}\\ {{\hat Y}_{{{\rm{C}}_0}}}\\ {{\hat Y}_{{\rm{P}}{{\rm{r}}_{\rm{2}}}}}\\ {{\hat Y}_{{{\rm{C}}_2}}} \end{aligned} \right).$

      在本方案中, 两个四波混频过程的强度增益G1G2是相互独立的. 在这里分别用

      $\left\langle {{{\hat X}_j}{{\hat X}_k}} \right\rangle = \frac{{\left\langle {{{\hat X}_j}{{\hat X}_k} + {{\hat X}_k}{{\hat X}_j}} \right\rangle }}{2} - \left\langle {{{\hat X}_j}} \right\rangle \left\langle {{{\hat X}_k}} \right\rangle $

      $\begin{split} & \left\langle {{{\hat Y}_j}{{\hat Y}_k}} \right\rangle = \frac{{\left\langle {{{\hat Y}_j}{{\hat Y}_k} + {{\hat Y}_k}{{\hat Y}_j}} \right\rangle }}{2} - \left\langle {{{\hat Y}_j}} \right\rangle \left\langle {{{\hat Y}_k}} \right\rangle \\ & \qquad \quad \left( {j, \;k = 1, 2, 3, 4} \right)\end{split}$

      来表示输出场的正交振幅算符和正交相位算符的协方差, 并且在j = k情况下时, 其代表的是相应的正交算符的方差. 此外, 正交振幅算符和正交相位算符相互正交[3], 即$\left\langle {\hat X\hat Y} \right\rangle = \left\langle {\hat Y\hat X} \right\rangle = {\rm{0}}$. 因此, 所有光束的协方差如下:

      $\begin{split} & \left\langle {\hat X_1^2} \right\rangle = \left\langle {\hat Y_1^2} \right\rangle = {G_1} + {{{g}}_1},\\ & \left\langle {\hat X_2^2} \right\rangle = \left\langle {\hat Y_2^2} \right\rangle = \eta ({G_1} + {{{g}}_1}) + (1 - \eta )({G_2} + {{{g}}_2}),\\ & \left\langle {\hat X_3^2} \right\rangle = \left\langle {\hat Y_3^2} \right\rangle = (1 - \eta )({G_1} + {{{g}}_1}) + \eta ({G_2} + {{{g}}_2}),\\ & \left\langle {\hat X_4^2} \right\rangle = \left\langle {\hat Y_4^2} \right\rangle = {G_2} + {g_2};\\[-15pt] \end{split}$

      $\begin{split} & \left\langle {{{\hat X}_1}{{\hat X}_2}} \right\rangle = - \left\langle {{{\hat Y}_1}{{\hat Y}_2}} \right\rangle = 2\sqrt {\eta {G_1}{g_1}} ,\\ & \left\langle {{{\hat X}_1}{{\hat X}_3}} \right\rangle = - \left\langle {{{\hat Y}_1}{{\hat Y}_3}} \right\rangle = - 2\sqrt {(1 - \eta ){G_1}{g_1}} ,\\ & \left\langle {{{\hat X}_1}{{\hat X}_4}} \right\rangle = \left\langle {{{\hat Y}_1}{{\hat Y}_4}} \right\rangle = 0,\\ & \left\langle {{{\hat X}_2}{{\hat X}_3}} \right\rangle \!=\! \left\langle {{{\hat Y}_2}{{\hat Y}_3}} \right\rangle \!=\! - \sqrt {\eta (1 - \eta )} ({G_1} \!+\! {{{g}}_1} \!-\! {G_2} \!-\! {{{g}}_2}),\\ & \left\langle {{{\hat X}_2}{{\hat X}_4}} \right\rangle = - \left\langle {{{\hat Y}_2}{{\hat Y}_4}} \right\rangle = 2\sqrt {(1 - \eta ){G_2}{g_2}} ,\\ & \left\langle {{{\hat X}_3}{{\hat X}_4}} \right\rangle = - \left\langle {{{\hat Y}_3}{{\hat Y}_4}} \right\rangle = 2\sqrt {\eta {G_2}{g_2}}.\\[-15pt] \end{split}$

      这一节推导了由线性分束器连接的两个独立的四波混频过程的输出场表达式, 接下来将重点研究该系统产生真正的四组份纠缠的可能性以及系统的纠缠结构.

    • 这一节用部分转置正定判据[34,35]来检验系统的四组份纠缠, 该判据是使用部分转置协方差矩阵σ的辛本征值来表示系统的纠缠特性, 如果所有的最小辛本征值v都小于1, 则表示系统是不可分离的. 首先研究了这个方案产生四组份纠缠的可能性. 在四组份情形中, 需要检验七种可能的两分形式, 其中包括四种1 × 3形式和三种2 × 2形式, 这七种二分形式如表1所列.

      数目二分形式数目二分形式
      11|23422|134
      33|12444|123
      512|34613|24
      714|23

      表 1  四组份态的七种二分形式

      Table 1.  Seven partitions of quadripartite state.

      为了方便, 假设线性分束器的透射率η = 0.5. 这里首先研究了四种1 × 3形式, 这种形式描述了每束光是如何与其他三束光纠缠的. 图2为四种1 × 3形式的最小辛本征值v, 它们都是G1G2的函数, 其中图2(a)—(d)分别是部分转置应用于光束C1, Pr1, Pr3和C3时的最小辛本征值v. 在图2中, 对于任何 G1 > 1, G2 > 1, 最小辛本征值v都小于1, 这表明每束光与剩余的三束光都存在纠缠, 并且纠缠随强度增益G1G2的增大而增强. 此外, 可以看出, 图2(b)图2(c)的最小辛本征值v是相同的, 这是系统的对称性导致的. 而当部分转置分别应用于光束C1和C3时, 其最小辛本征值v与另外一个四波混频过程的强度增益无关, 这说明它们的纠缠来源于对应的四波混频过程, 而与另外一个四波混频过程无关.

      图  2  四种1 × 3情形的最小辛本征值v, 其为强度增益G1G2的函数 (a) C1被部分转置; (b) Pr1被部分转置; (c) Pr3被部分转置; (d) C3被部分转置

      Figure 2.  The smallest symplectic eigenvalue v of all 1 × 3 scenarios, as a function of the power gains G1 and G2: (a) C1 is partially transposed; (b) Pr1 is partially transposed; (c) Pr3 is partially transposed; (d) C3 is partially transposed.

      同样地, 也研究了三种2 × 2形式, 这种形式是将部分转置应用于四束光中的任意两束光. 图3为三种2 × 2形式的最小辛本征值v, 其为强度增益G1G2的函数, 其中图3(a)(c)分别是部分转置应用于光束C1和Pr1, C1和Pr3以及C1和C3时的最小辛本征值v. 由图3可得, 所有的最小辛本征值v都小于1, 这表明在四组份系统中, 每对光束都与其他光束纠缠, 并且强度增益G1G2越大, 纠缠越强. 值得注意的是, 当部分转置应用于一对光束C1C3时, 其最小辛本征值v与强度增益G1无关, 并且只受强度增益G2的影响. 换句话说, 它们的纠缠来源于第二个四波混频过程, 而与第一个四波混频过程无关.

      图  3  三种2 × 2情形的最小辛本征值v, 其为强度增益G1G2的函数 (a) C1和Pr1被部分转置; (b) C1和Pr3被部分转置; (c) C1和C3被部分转置

      Figure 3.  The smallest symplectic eigenvalues v of all 2 × 2 scenarios, as a function of the power gains G1 and G2: (a) C1 and Pr1 arepartially transposed; (b) C1 and Pr3 are partially transposed; (c) C1 and C3 are partially transposed.

      基于上述的分析可知, 在线性分束器的透射率η = 0.5以及强度增益G1 > 1, G2 > 1时, 系统所有的最小辛本征值v都小于1, 并且随着强度增益G1G2增加, 纠缠也在增强, 这表明本文系统可以产生真正的四组份纠缠.

      在研究了强度增益G1G2对四组份纠缠的影响之后, 研究线性分束器的透射率η对四组份纠缠的影响也十分必要. 在这里, 为了研究透射率η对四组份纠缠的影响, 假设两个四波混频过程的强度增益G1 = G2 = 3. 图4中所有的最小辛本征值v都是透射率η的函数, 其中图4(a)(g)分别是部分转置应用于光束C1, Pr1, Pr3, C3以及C1和Pr1, C1和Pr3, C1和C3时的最小辛本征值v. 由图4可得, 在透射率η ≠ 0或η ≠ 1时, 系统存在真正的四组份纠缠. 此外, 有些最小辛本征值v与透射率η无关, 有些则与其紧密相关. 这是因为线性分束器仅连接了两束光, 这会导致部分最小辛本征值v与线性分束器的透射率η相关.

      图  4  线性分束器的透射率η对四组份态的最小辛本征值v的影响 (a) C1被部分转置; (b) Pr1被部分转置; (c) Pr3被部分转置; (d) C3被部分转置; (e) C1和Pr1被部分转置; (f) C1和Pr3被部分转置; (g) C1和C3被部分转置

      Figure 4.  Effect of the transmissivity of the linear beam splitter on the quadripartite entanglement of the system: (a) C1 is partially transposed; (b) Pr1 is partially transposed; (c) Pr3 is partially transposed; (d) C3 is partially transposed; (e) C1 and Pr1 are partially transposed; (f) C1 and Pr3 are partially transposed; (g) C1 and C3 are partially transposed.

    • 到目前为止, 我们已经证明该方案可以产生真正的四组份纠缠. 同时, 研究系统所产生的三组份纠缠和两组份纠缠将有助于进一步了解该系统的纠缠结构. 在这里, 假设线性分束器的透射率η = 0.5. 对于三组份态的纠缠特性, 我们是通过去除系统所产生的四束光中的一束光进行研究的. 很明显, 这里有四种可能的三组份态, 并且每种三组份态有三种1 × 2形式, 因此, 总共需要检验12个1 × 2形式的最小辛本征值v. 图5给出了所有三组份态的最小辛本征值v, 其中图5(a)(c)是由光束C1, Pr1 和Pr3组成的三组份态的最小辛本征值v, 图5(d)(f)是由光束C1, Pr1和C3组成的三组份态的最小辛本征值v, 图5(g)(i)是由光束C1, Pr3和C3组成的三组份态的最小辛本征值v, 图5(j)(l)是由光束Pr1, Pr3和C3组成的三组份态的最小辛本征值v. 由图5可得, 在G1 > 1, G2 > 1时, 所有三组份态的最小辛本征值v都小于1, 并且强度增益G1G2越大, 纠缠越强, 这表明输出场的四束光在去除其中的任意一束光之后仍存在真正的三组份纠缠. 此外, 由于系统的对称性, 三组份态(C1, Pr1, C3)和(C1, Pr3, C3)的最小辛本征值v是一样的. 而三组份态(C1, Pr1, Pr3)和(Pr1, Pr3, C3)的最小辛本征值v却是不一样的, 其最小辛本征值v分别受到两个独立的四波混频过程的强度增益影响.

      图  5  四个三组份态的最小辛本征值v, 其为强度增益G1G2的函数 (a)−(c)是由C1, Pr1和Pr3组成的三组份态的最小辛本征值v; (d)−(f)是由C1, Pr1和C3组成的三组份态的最小辛本征值v; (g)−(i)是由C1, Pr3和C3组成的三组份态的最小辛本征值v; (j)−(l)是由Pr1, Pr3和C3组成的三组份态的最小辛本征值v

      Figure 5.  The smallest symplectic eigenvalues v of all tripartite states as a function of power gains G1 and G2: (a)−(c) The smallest symplectic eigenvalues v of tripartite state composed of C1, Pr1 and Pr3; (d)−(f) the smallest symplectic eigenvalues v of tripartite state composed of C1, Pr1 and C3; (g)−(i) the smallest symplectic eigenvalues v of tripartite state composed of C1, Pr3 and C3; (j)−(l) the smallest symplectic eigenvalues v of tripartite state composed of Pr1, Pr3 and C3.

    • 接下来将通过去除系统输出场中的任意两束光来研究剩余两束光的纠缠特性. 这里需要检验六种不同的两组份态的纠缠特性. 图6为所有的两组份态的最小辛本征值v, 其中图6(a)(f)分别是由C1和Pr1, C1和Pr3, C1和C3, Pr1和Pr3, Pr1和C3以及Pr3和C3组成的两组份态的最小辛本征值v. 如图6所示, 有些两组份态(例如由C1和C3组成的两组份态)总不存在纠缠, 这是由于两个四波混频过程是相互独立而导致的. 而有些两组份态(例如由C1和Pr1组成的两组份态)的纠缠特性则随着强度增益G1G2的增大而逐渐减小直至消失, 这是因为随着强度增益G1G2的增大, 线性分束器引入的额外噪声迅速增大, 从而导致纠缠消失.

      图  6  六种两组份态的最小辛本征值v, 其为强度增益G1G2的函数 (a) 由C1和Pr1组成的两组份态的最小辛本征值v; (b) 由C1和Pr3组成的两组份态的最小辛本征值v; (c) 由C1和C3组成的两组份态的最小辛本征值v; (d) 由Pr1和Pr3组成的两组份态的最小辛本征值v; (e) 由Pr1和C3组成的两组份态的最小辛本征值v; (f) 由Pr3和C3组成的两组份的最小辛本征值v

      Figure 6.  The smallest symplectic eigenvalues v of all bipartite states as a function of power gains G1 and G2: (a) The smallest symplectic eigenvalues v of bipartite state composed of C1 and Pr1; (b) the smallest symplectic eigenvalues v of bipartite state composed of C1 and Pr3; (c) the smallest symplectic eigenvalues v of bipartite state composed of C1 and C3; (d) the smallest symplectic eigenvalues v of bipartite state composed of Pr1 and Pr3; (e) the smallest symplectic eigenvalues v of bipartite state composed of Pr1 and C3; (f) the smallest symplectic eigenvalues v of bipartite state composed of Pr3 and C3.

      由于线性分束器的引入, 该系统存在真正的四组份纠缠和三组份纠缠, 并且随着强度增益G1G2的增大, 纠缠也会逐渐增强. 至于系统两组份态的纠缠特性, 有些两组份态总是可分离的. 而有些两组份态的纠缠的存在与否却是有条件的, 即它们的纠缠特性受强度增益G1G2的影响, 随着两个四波混频过程的强度增益G1G2的增大, 这些两组份纠缠都会逐渐减小直至消失.

    • 本文从理论上提出了一种利用两个独立的四波混频过程和线性分束器产生真正的四组份纠缠的方案, 两个独立的四波混频过程由线性分束器连接. 首先应用部分转置正定判据研究了强度增益对四组份纠缠的影响, 结果表明在整个增益区域内都存在四组份纠缠, 并且强度增益越大, 纠缠程度越强, 这表明本文系统能够生成真正的四组份纠缠. 然后研究了线性分束器的透射率对系统的四组份纠缠的影响, 发现只要线性分束器的透射率不为0或1, 系统也存在真正的四组份纠缠. 最后, 通过研究该系统中可能存在的三组份纠缠和两组份纠缠揭示了系统的纠缠结构. 我们的理论研究为实验上利用原子系综四波混频过程产生真正的四组份纠缠提供了一个有效的方案.

参考文献 (36)

目录

    /

    返回文章
    返回