搜索

x

留言板

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

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

Jaynes-Cummings晶格模型和Rabi晶格模型的量子相变

尤冰凌 刘雪莹 成书杰 王晨 高先龙

引用本文:
Citation:

Jaynes-Cummings晶格模型和Rabi晶格模型的量子相变

尤冰凌, 刘雪莹, 成书杰, 王晨, 高先龙

The quantum phase transition in the Jaynes-Cummings lattice model and the Rabi lattice model

You Bing-Ling, Liu Xue-Ying, Cheng Shu-Jie, Wang Chen, Gao Xian-Long
PDF
HTML
导出引用
  • 采用平均场近似的方法, 分别研究了Jaynes-Cummings晶格模型和Rabi晶格模型的量子相变: Mott绝缘体相-超流体相量子相变, 探索了光的聚束-反聚束行为, 研究了Kerr非线性作用对量子相变与光子统计特征的影响. 研究结果表明, 在Rabi晶格模型中二能级原子和光子相互作用强度g和格点之间光子跃迁强度J的增大会使晶格体系从Mott绝缘体相向超流体相转变, 同时, 光子统计行为由聚束转变为反聚束, 而Kerr非线性强度的增大抑制了Mott绝缘体相-超流体相相变, 但促进了光子聚束与反聚束之间的转变.
    We use the mean field approximation method to study the quantum phase transitions of the Jaynes-Cummings lattice model and the Rabi lattice model. The effective Hamiltonians are obtained for the JC and Rabi model including the Kerr nonlinear term. Numerically we diagonalized the Hamiltonian matrix and calculated the superfluidity order parameter and the two-photon correlation function by solving the iteration equations.We have explored the Mott insulating-superfluid quantum phase transition, the bunching-antibunching behavior of light, and the effect of Kerr nonlinear term on the quantum phase transition and photon statistical characteristics. Our results show that in the JC lattice model, by increasing J, a quantum phase transition takes place and the system is driven to a superfluid phase. The phase boundaries of the Mott lobes are N-dependent. However the photon will always be in a bunching statistical behavior irrelevant of the coupling strength between the two-level atom and the phonton and the nonlinear Kerr effect.In the Rabi lattice model, the anti-rotating wave term breaks Mott-lobe structure of the phase diagram and the increase of the two-level atom and photon interaction strength g and the photon transition strength J between the lattices drive the system from the Mott insulating phase to the superfluid phase. The photon statistical behavior changes from the bunching to the antibunching one when considering the anti-rotating wave term, which is important in the strongly coupled systems. Most interestingly, the increase of the Kerr nonlinear coefficient will inhibit the Mott insulating phase-superfluid phase transition, but favor the superfluid phase and the transition from the bunching to anti-bunching statistics.
      通信作者: 高先龙, gaoxl@zjnu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11774316, 11704093)资助的课题
      Corresponding author: Gao Xian-Long, gaoxl@zjnu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 11774316, 11704093).
    [1]

    Hartmann M J, Brandao F G S L, Plenio M B 2006 Nat. Phys. 2 849Google Scholar

    [2]

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

    [3]

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

    [4]

    Rossini D, Fazio R 2007 Phys. Rev. Lett. 99 186401Google Scholar

    [5]

    Aichhorn M, Hohenadler M, Tahan C, Littlewood P B 2008 Phys. Rev. Lett. 100 216401Google Scholar

    [6]

    Na N, Utsunomiya S, Tian L, Yamamoto Y 2008 Phys. Rev. A 77 031803Google Scholar

    [7]

    Carusotto I, Gerace D, Türeci H E, De Liberato S, Ciuti C, and Imamoğlu A 2009 Phys. Rev. Lett. 103 033601Google Scholar

    [8]

    Schmidt S, Blatter G 2009 Phys. Rev. Lett. 103 086403Google Scholar

    [9]

    Koch J, Le Hur K 2009 Phys. Rev. A 80 023811Google Scholar

    [10]

    Pippan P, Evertz H G, Hohenadler M 2009 Phys. Rev. A 80 033612Google Scholar

    [11]

    Ferretti S, Andreani L C, Türeci H E, Gerace D 2010 Phys. Rev. A 82 013841Google Scholar

    [12]

    Umucalilar R O, Carusotto I 2011 Phys. Rev. A 84 043804Google Scholar

    [13]

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

    [14]

    Tian L, Carmichael H J 1992 Phys. Rev. A 46 R6801Google Scholar

    [15]

    Imamoḡlu A, Schmidt H, Woods G, Deutsch M 1997 Phys. Rev. Lett. 79 1467Google Scholar

    [16]

    Rebic S, Tan S M, Parkins A S, Walls D F 1999 J. Opt. B 1 490Google Scholar

    [17]

    Bloch I, Dalibard J, Zwerger W 2008 Rev. Mod. Phys. 80 885Google Scholar

    [18]

    Schmidt S, Koch J 2013 Ann. Phys. 525 395Google Scholar

    [19]

    Greiner M, Mandel O, Esslinger T, Hänsch T W, Bloch I 2002 Nature 415 39

    [20]

    Lundqvist S, Nilsson N B 1989 Physics of Low-dimensional Systems (Sweden: World Scientific) pp89−95

    [21]

    Fisher M P, Weichman P B, Grinstein G, Fisher D S 1989 Phys. Rev. B 40 546Google Scholar

    [22]

    van Oosten D, van Der Straten P, Stoof H T C 2001 Phys. Rev. A 63 053601Google Scholar

    [23]

    van Oosten D, van Der Straten P, Stoof H T C 2003 Phys. Rev. A 67 033606Google Scholar

    [24]

    Sheshadri K, Krishnamurthy H R, Pandit R, Ramakrishnan T V 1993 Europhys. Lett. 22 257Google Scholar

    [25]

    Krauth W, Trivedi N 1991 Europhys. Lett. 14 627Google Scholar

    [26]

    Krauth W, Trivedi N, Ceperley D 1991 Phys. Rev. Lett. 67 2307Google Scholar

    [27]

    Xie Z W, Liu W M 2004 Phys. Rev. A 70 045602Google Scholar

    [28]

    Albus A, Illuminati F, Eisert J 2003 Phys. Rev. A 68 023606Google Scholar

    [29]

    Lewenstein M, Santos L, Baranov M A, Fehrmann H 2004 Phys. Rev. Lett. 92 050401Google Scholar

    [30]

    Illuminati F, Albus A 2004 Phys. Rev. Lett. 93 090406Google Scholar

    [31]

    Cramer M, Eisert J, Illuminati F 2004 Phys. Rev. Lett. 93 190405Google Scholar

    [32]

    Fehrmann H, Baranov M A, Damski B, Lewenstein M, Santos L 2004 Opt. Commun. 243 23

    [33]

    Littlewood P B, Eastham P R, Keeling J M J, Marchetti F M, Simons B D, Szymanska M H 2004 J. Phys. Condens. Matter 16 S3597Google Scholar

    [34]

    He Y, Zhu X, Mihalache D, Liu J, Chen Z 2012 Phys. Rev. A 85 013831Google Scholar

    [35]

    Eguchi K, Takagi Y, Nakagawa T, Yokoyama T 2012 Phys. Rev. B 85 174415Google Scholar

    [36]

    Vitali D, Fortunato M, Tombesi P 2000 Phys. Rev. Lett. 85 445Google Scholar

    [37]

    Angelakis D G, Dai L, Kwek L C 2010 Europhys. Lett. 91 10003Google Scholar

    [38]

    Patargias N, Bartzis V, Jannussis A 1995 Phys. Scr. 52 554Google Scholar

    [39]

    Bu S P, Zhang G F, Liu J, Chen Z Y 2008 Phys. Scr. 78 065008Google Scholar

    [40]

    Cordero S, Récamier J 2011 J. Phys. B: At. Mol. Opt. Phys. 44 135502Google Scholar

    [41]

    Schmidt H, Imamoğlu A 1996 Opt. Lett. 21 1936Google Scholar

    [42]

    Harris S E, Hau L V 1999 Phys. Rev. Lett. 82 4611Google Scholar

    [43]

    Niu Y, Gong S 2006 Phys. Rev. A 73 053811Google Scholar

    [44]

    Glauber R J 1963 Phys. Rev. 130 2529Google Scholar

    [45]

    Gomes C B C, Almeida F A G, Souza A M C 2016 Phys. Lett. A 38 1799

  • 图 1  (a), (b)平均场近似下, 不同晶格模型关于超流序参量$ A = \braket{a} $$ J \text {-} g $相图 (a) JC晶格模型; (b) Rabi 晶格模型. 横坐标为格点之间的光子跃迁强度$ J $, 纵坐标为二能级原子和光子相互作用强度$ g $, 横纵坐标的单位为$ \omega_0 $, 颜色条表示超流序参量$ A = \braket{a} $的大小. 深蓝色表示Mott绝缘相, 浅黄色表示超流体相. 其他参量取值为: $ \mathop{\omega_{0} = \omega_{1}} = 1 $, 光子截断数$ N = 20 $. (c), (d)对于不同的$ J $, 不同晶格模型的超流序参量$ A $$ g $变化的图像 (c) JC晶格模型; (d) Rabi晶格模型

    Fig. 1.  (a), (b) Under the mean field approximation, the $ J \text {-} g $ phase diagram of different lattice models with respect to the superfluid order parameter $ A = \braket{a} $: (a) JC lattice model; (b) Rabi lattice model. The abscissa is the photon transition intensity $ J $ between the lattice, the ordinate is the two-level atom and photon interaction strength $ g $, the unit of the abscissa and the ordinate is $ \omega_0 $, and the color bar represents the value of the superfluid order parameter $ A = \braket{a} $. Dark blue indicates Mott insulating phase, and light yellow indicates superfluid phase. Other parameters are taken as $ \mathop{\omega_{0} = \omega_{1}} = 1 $, and the number of the photon truncation $ N = 20 $. (c), (d) For different $ J $, the superfluid order parameter $ A $ of different lattice models varies with $ g $: (c) JC lattice model; (d) Rabi lattice model.

    图 2  平均场近似下, 不同晶格模型关于二阶关联函数$ g^{2}(0) $$J\text {-} g$相图 (a) JC晶格模型; (b) Rabi晶格模型. 横坐标为格点之间的光子跃迁强度$ J $, 纵坐标为二能级原子和光子相互作用强度$ g $, 横纵坐标的单位为$ \omega_0 $, 颜色条表示二阶关联函数$ g^{2}(0) $ 的值. 其他参量取值为: $ \mathop{\omega_{0} = \omega_{1}} = 1 $, 光子截断数$ N = 20 $

    Fig. 2.  Under the mean field approximation, the $ J \text {-} g $ phase diagram of different lattice models with respect to the second-order correlation function $ g^{2}(0) $: (a) JC lattice model; (b) Rabi lattice model. The abscissa is the photon transition intensity $ J $ between the lattice, the ordinate is the two-level atom and photon interaction strength $ g $, the unit of the abscissa and the ordinate is $ \omega_0 $, the color bar is represented by the value of the second-order correlation function $ g^{2}(0) $. $ \mathop{\omega_{0} = \omega_{1}} = 1 $, and the number of photon truncation $ N = 20 $.

    图 3  Kerr效应下不同晶格模型关于超流序参量$ A = \braket{a} $$ \kappa \text {-} g $相图 (a) JC晶格模型; (b) Rabi 晶格模型. 横坐标为Kerr非线性强度$ \kappa $, 纵坐标为二能级原子和光子相互作用强度$ g $, 横纵坐标的单位为$ \omega_0 $, 颜色条表示超流序参量$ A $的大小. 其他参量取值为: $ \mathop{\omega_{0} = \omega_{1}} = 1 $, $ J = 0.05 $, 光子截断数$ N = 20 $

    Fig. 3.  The $ \kappa \text {-} g $ phase diagram of different lattice models under the Kerr effect with respect to the superfluid order parameter $ A = \braket{a} $: (a) JC lattice model; (b) Rabi lattice model. The abscissa is the Kerr nonlinear intensity $ \kappa $, the ordinate is the two-level atom and photon interaction strength $ g $, the unit of the abscissa and the ordinate is $ \omega_0 $, and the color bar represents the value of the superfluid order parameter $ A $. Other parameters are taken as $ \mathop{\omega_{0} = \omega_{1}} = 1 $, $ J = 0.05 $, and the number of photon truncation $ N = 20 $.

    图 4  Kerr效应下不同晶格模型关于二阶关联函数$ g^{2}(0) $$ \kappa \text {-} g $相图 (a) JC晶格模型; (b) Rabi晶格模型. 横坐标为Kerr非线性强度$ \kappa $, 纵坐标为二能级原子和光子相互作用强度$ g $, 横纵坐标的单位为$ \omega_0 $, 颜色条表示二阶关联函数$ g^{2}(0) $. 其他参量取值为: $ \mathop{\omega_{0} = \omega_{1}} = 1 $, $ J = 0.05 $, 光子截断数$ N = 20 $

    Fig. 4.  The $ \kappa \text {-} g $ phase diagram of different lattice models under the Kerr effect with respect to the second-order correlation function $ g^{2}(0) $: (a) JC lattice model; (b) Rabi lattice model. The abscissa is the Kerr nonlinear intensity $ \kappa $, the ordinate is the two-level atom and photon interaction strength $ g $, the unit of the abscissa and the ordinate is $ \omega_0 $, and the color bar represents the value of second-order correlation function $ g^{2}(0) $. Other parameters are taken as $ \mathop{\omega_{0} = \omega_{1}} = 1 $, $ J = 0.05 $, and the number of photon truncation $ N = 20 $.

  • [1]

    Hartmann M J, Brandao F G S L, Plenio M B 2006 Nat. Phys. 2 849Google Scholar

    [2]

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

    [3]

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

    [4]

    Rossini D, Fazio R 2007 Phys. Rev. Lett. 99 186401Google Scholar

    [5]

    Aichhorn M, Hohenadler M, Tahan C, Littlewood P B 2008 Phys. Rev. Lett. 100 216401Google Scholar

    [6]

    Na N, Utsunomiya S, Tian L, Yamamoto Y 2008 Phys. Rev. A 77 031803Google Scholar

    [7]

    Carusotto I, Gerace D, Türeci H E, De Liberato S, Ciuti C, and Imamoğlu A 2009 Phys. Rev. Lett. 103 033601Google Scholar

    [8]

    Schmidt S, Blatter G 2009 Phys. Rev. Lett. 103 086403Google Scholar

    [9]

    Koch J, Le Hur K 2009 Phys. Rev. A 80 023811Google Scholar

    [10]

    Pippan P, Evertz H G, Hohenadler M 2009 Phys. Rev. A 80 033612Google Scholar

    [11]

    Ferretti S, Andreani L C, Türeci H E, Gerace D 2010 Phys. Rev. A 82 013841Google Scholar

    [12]

    Umucalilar R O, Carusotto I 2011 Phys. Rev. A 84 043804Google Scholar

    [13]

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

    [14]

    Tian L, Carmichael H J 1992 Phys. Rev. A 46 R6801Google Scholar

    [15]

    Imamoḡlu A, Schmidt H, Woods G, Deutsch M 1997 Phys. Rev. Lett. 79 1467Google Scholar

    [16]

    Rebic S, Tan S M, Parkins A S, Walls D F 1999 J. Opt. B 1 490Google Scholar

    [17]

    Bloch I, Dalibard J, Zwerger W 2008 Rev. Mod. Phys. 80 885Google Scholar

    [18]

    Schmidt S, Koch J 2013 Ann. Phys. 525 395Google Scholar

    [19]

    Greiner M, Mandel O, Esslinger T, Hänsch T W, Bloch I 2002 Nature 415 39

    [20]

    Lundqvist S, Nilsson N B 1989 Physics of Low-dimensional Systems (Sweden: World Scientific) pp89−95

    [21]

    Fisher M P, Weichman P B, Grinstein G, Fisher D S 1989 Phys. Rev. B 40 546Google Scholar

    [22]

    van Oosten D, van Der Straten P, Stoof H T C 2001 Phys. Rev. A 63 053601Google Scholar

    [23]

    van Oosten D, van Der Straten P, Stoof H T C 2003 Phys. Rev. A 67 033606Google Scholar

    [24]

    Sheshadri K, Krishnamurthy H R, Pandit R, Ramakrishnan T V 1993 Europhys. Lett. 22 257Google Scholar

    [25]

    Krauth W, Trivedi N 1991 Europhys. Lett. 14 627Google Scholar

    [26]

    Krauth W, Trivedi N, Ceperley D 1991 Phys. Rev. Lett. 67 2307Google Scholar

    [27]

    Xie Z W, Liu W M 2004 Phys. Rev. A 70 045602Google Scholar

    [28]

    Albus A, Illuminati F, Eisert J 2003 Phys. Rev. A 68 023606Google Scholar

    [29]

    Lewenstein M, Santos L, Baranov M A, Fehrmann H 2004 Phys. Rev. Lett. 92 050401Google Scholar

    [30]

    Illuminati F, Albus A 2004 Phys. Rev. Lett. 93 090406Google Scholar

    [31]

    Cramer M, Eisert J, Illuminati F 2004 Phys. Rev. Lett. 93 190405Google Scholar

    [32]

    Fehrmann H, Baranov M A, Damski B, Lewenstein M, Santos L 2004 Opt. Commun. 243 23

    [33]

    Littlewood P B, Eastham P R, Keeling J M J, Marchetti F M, Simons B D, Szymanska M H 2004 J. Phys. Condens. Matter 16 S3597Google Scholar

    [34]

    He Y, Zhu X, Mihalache D, Liu J, Chen Z 2012 Phys. Rev. A 85 013831Google Scholar

    [35]

    Eguchi K, Takagi Y, Nakagawa T, Yokoyama T 2012 Phys. Rev. B 85 174415Google Scholar

    [36]

    Vitali D, Fortunato M, Tombesi P 2000 Phys. Rev. Lett. 85 445Google Scholar

    [37]

    Angelakis D G, Dai L, Kwek L C 2010 Europhys. Lett. 91 10003Google Scholar

    [38]

    Patargias N, Bartzis V, Jannussis A 1995 Phys. Scr. 52 554Google Scholar

    [39]

    Bu S P, Zhang G F, Liu J, Chen Z Y 2008 Phys. Scr. 78 065008Google Scholar

    [40]

    Cordero S, Récamier J 2011 J. Phys. B: At. Mol. Opt. Phys. 44 135502Google Scholar

    [41]

    Schmidt H, Imamoğlu A 1996 Opt. Lett. 21 1936Google Scholar

    [42]

    Harris S E, Hau L V 1999 Phys. Rev. Lett. 82 4611Google Scholar

    [43]

    Niu Y, Gong S 2006 Phys. Rev. A 73 053811Google Scholar

    [44]

    Glauber R J 1963 Phys. Rev. 130 2529Google Scholar

    [45]

    Gomes C B C, Almeida F A G, Souza A M C 2016 Phys. Lett. A 38 1799

  • [1] 孙振辉, 胡丽贞, 徐玉良, 孔祥木. 准一维混合自旋(1/2,5/2)Ising-XXZ模型的量子相干和互信息. 物理学报, 2023, 0(0): . doi: 10.7498/aps.72.20230381
    [2] 陈西浩, 夏继宏, 李孟辉, 翟福强, 朱广宇. 自旋-1/2量子罗盘链的量子相与相变. 物理学报, 2022, 71(3): 030302. doi: 10.7498/aps.71.20211433
    [3] 陈西浩, 夏继宏, 李孟辉, 翟福强, 朱广宇. 自旋-1/2量子罗盘链的量子相与相变. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211433
    [4] 周晓凡, 樊景涛, 陈刚, 贾锁堂. 光学腔中一维玻色-哈伯德模型的奇异超固相. 物理学报, 2021, 70(19): 193701. doi: 10.7498/aps.70.20210778
    [5] 保安. 各向异性ruby晶格中费米子体系的Mott相变. 物理学报, 2021, 70(23): 230305. doi: 10.7498/aps.70.20210963
    [6] 文凯, 王良伟, 周方, 陈良超, 王鹏军, 孟增明, 张靖. 超冷87Rb原子在二维光晶格中Mott绝缘态的实验实现. 物理学报, 2020, 69(19): 193201. doi: 10.7498/aps.69.20200513
    [7] 陈爱民, 刘东昌, 段佳, 王洪雷, 相春环, 苏耀恒. 含有Dzyaloshinskii-Moriya相互作用的自旋1键交替海森伯模型的量子相变和拓扑序标度. 物理学报, 2020, 69(9): 090302. doi: 10.7498/aps.69.20191773
    [8] 伊天成, 丁悦然, 任杰, 王艺敏, 尤文龙. 具有Dzyaloshinskii-Moriya相互作用的XY模型的量子相干性. 物理学报, 2018, 67(14): 140303. doi: 10.7498/aps.67.20172755
    [9] 陈西浩, 王秀娟. 一维扩展量子罗盘模型的拓扑序和量子相变. 物理学报, 2018, 67(19): 190301. doi: 10.7498/aps.67.20180855
    [10] 宋加丽, 钟鸣, 童培庆. 横场中具有周期性各向异性的一维XY模型的量子相变. 物理学报, 2017, 66(18): 180302. doi: 10.7498/aps.66.180302
    [11] 赵红霞, 赵晖, 陈宇光, 鄢永红. 一维扩展离子Hubbard模型的相图研究. 物理学报, 2015, 64(10): 107101. doi: 10.7498/aps.64.107101
    [12] 俞立先, 梁奇锋, 汪丽蓉, 朱士群. 双模Dicke模型的一级量子相变. 物理学报, 2014, 63(13): 134204. doi: 10.7498/aps.63.134204
    [13] 赵建辉, 王海涛. 应用多尺度纠缠重整化算法研究量子自旋系统的量子相变和基态纠缠. 物理学报, 2012, 61(21): 210502. doi: 10.7498/aps.61.210502
    [14] 赵建辉. 应用约化密度保真度确定自旋为1的一维量子 Blume-Capel模型的基态相图. 物理学报, 2012, 61(22): 220501. doi: 10.7498/aps.61.220501
    [15] 石筑一, 童 红, 石筑亚, 张春梅, 赵行知, 倪绍勇. 转动诱发原子核量子相变的一种可能途径. 物理学报, 2007, 56(3): 1329-1333. doi: 10.7498/aps.56.1329
    [16] 李 园, 李 刚, 张玉驰, 王晓勇, 王军民, 张天才. 计数率和分辨时间对光场统计性质测量的影响——单探测器直接测量的实验分析. 物理学报, 2006, 55(11): 5779-5783. doi: 10.7498/aps.55.5779
    [17] 王晓波, 黄 涛, 邵军虎, 降雨强, 肖连团, 贾锁堂. 利用光子统计概率测量单分子光子源的信号背景比. 物理学报, 2005, 54(2): 617-621. doi: 10.7498/aps.54.617
    [18] 梁文青, 储开芹, 张智明, 谢绳武. 超冷V型三能级原子注入的微波激射:原子相干性对腔场光子统计的影响. 物理学报, 2001, 50(12): 2345-2355. doi: 10.7498/aps.50.2345
    [19] 顾樵. Jaynes-Cummings模型的量子统计性质. 物理学报, 1989, 38(5): 735-744. doi: 10.7498/aps.38.735
    [20] 李孝申. 量子统计的级联三能级Jaynes-Cummings模型. 物理学报, 1985, 34(6): 833-840. doi: 10.7498/aps.34.833
计量
  • 文章访问数:  3954
  • PDF下载量:  269
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-12-06
  • 修回日期:  2021-01-04
  • 上网日期:  2021-05-06
  • 刊出日期:  2021-05-20

Jaynes-Cummings晶格模型和Rabi晶格模型的量子相变

  • 浙江师范大学物理系, 金华 321004
  • 通信作者: 高先龙, gaoxl@zjnu.edu.cn
    基金项目: 国家自然科学基金(批准号: 11774316, 11704093)资助的课题

摘要: 采用平均场近似的方法, 分别研究了Jaynes-Cummings晶格模型和Rabi晶格模型的量子相变: Mott绝缘体相-超流体相量子相变, 探索了光的聚束-反聚束行为, 研究了Kerr非线性作用对量子相变与光子统计特征的影响. 研究结果表明, 在Rabi晶格模型中二能级原子和光子相互作用强度g和格点之间光子跃迁强度J的增大会使晶格体系从Mott绝缘体相向超流体相转变, 同时, 光子统计行为由聚束转变为反聚束, 而Kerr非线性强度的增大抑制了Mott绝缘体相-超流体相相变, 但促进了光子聚束与反聚束之间的转变.

English Abstract

    • 光与物质之间的作用是自然界中最基本的相互作用过程之一, 是理解多种物理现象的基础. Jaynes-Cummings (JC)模型和Rabi模型是描述二能级系统与光场相互作用的核心模型. 随着科技的发展, 人们对原子和固态系统与光相互作用的控制水平不断提高, 在光晶格系统上实现了一系列新的量子物质的多体效应[1-12]. 在这些系统中, 光和物质之间的相互作用有着非常重要的作用, 例如由单模的电磁谐振器和二能级系统相互作用形成的基本腔量子电动力学系统. 原子与光子之间的耦合会导致有效的非线性光子-光子排斥. Birnbaum等[13]研究表明这种光子排斥会阻止一个以上的光子进入空腔. 光子阻塞效应[14-16]与腔内光子隧穿之间的相互作用导致腔量子电动力学系统耦合形成晶格, 形成大型的Bose-Hubbard模型[17]. 最近, 有研究者提出了一种光晶格中存在JC模型且腔之间存在光子跃迁项的模型[1,8,18], 称之为JC晶格模型, 它描述了原子-光子之间的强耦合与腔间光子跃迁的竞争. 已经有一些文章使用JC晶格模型研究从Mott绝缘体相到超流体相的量子相图[19,20]问题. 同样, 也可以在Rabi晶格模型中进行相应的研究.

      近几年, 平均场近似[2,6]为从Mott绝缘体相到超流体相的相变研究提供了方法, 这与Bose-Hubbard模型[21]中的研究类似. 与蒙特卡罗方法[22,23]相比, 平均场近似可以定性地描述大的系统[24-26], 同时, 平均场近似的方法已经扩展到偶极玻色子[27]和玻色子-费米子原子混合系统[28-32], 以及激子-极化子凝聚[33]等众多领域. 由于可以在光晶格中观察到非线性光学效应[34,35], 且这种效应具有产生纠缠态的可能, 因此这种效应引起了人们的极大兴趣, 特别是在量子信息领域[36,37]. 有一些实验已经对涉及Kerr效应的JC模型的特性进行了研究. 特别是探索了由外部激光源驱动非线性腔发出光子的统计特性[38-40]. 与此同时, 控制Kerr非线性的方法也已经得到了解决[41-43], 因此进一步研究Kerr非线性对Rabi晶格模型的影响具有重要的理论意义.

      本文采用平均场近似的方法, 研究了JC晶格模型和Rabi晶格模型中Mott绝缘体相到超流体相的量子相图, 以及从聚束到反聚束的转变行为. 并在此基础上考虑Kerr非线性项, 分析Kerr效应对JC晶格模型和Rabi晶格模型中光量子相变与统计特征的影响.

    • Rabi模型描述了一个二能级系统和光场之间的基本相互作用,

      $ H_{\rm{Rabi}} = \omega_{0}a^{\dagger}a+\frac{\omega_{1}}{2}{{\sigma}}_{z} +g({a}^{\dagger}+{a})({{\sigma}}_-+{{\sigma}}_+), $

      这里, $ \omega_{0} $是光场的频率, $ \omega_{1} $是二能级原子的能级差, $ g $是光子与二能级原子的耦合强度, $ a^{\dagger}(a) $为光子的产生(湮灭)算符, $ {{\sigma}}_{z} $是二能级原子的泡利矩阵, $ {{\sigma}}_{\pm} $是泡利矩阵的升降算符.

      Rabi模型的相互作用项$ ({a}^{\dagger}+{a})({{\sigma}}_{-}+{{\sigma}}_{+}) $可以写成旋波项$ a^{\dagger}{{\sigma}}_{-}+a{{\sigma}}_{+} $和反旋波项$ a{{\sigma}}_{-}+a^{\dagger}{{\sigma}}_{+} $之和, 其中$ a^{\dagger}{{\sigma}}_{-} $表示释放一个光子的同时二能级原子从激发态跃迁到基态, $ a{{\sigma}}_{+} $表示吸收一个光子的同时二能级原子从基态跃迁到激发态, $ a{{\sigma}}_{-} $表示吸收一个光子的同时二能级原子从激发态跃迁到基态, $ a^{\dagger}{{\sigma}}_{+} $表示释放一个光子的同时二能级原子从基态跃迁到激发态. 在光场频率和原子频率近共振、弱耦合的近似条件下, 反旋波项随时间演化会迅速衰减到零, 可以被忽略, 即在旋波近似下, (1)式变为

      $ H_{\rm{JC}} = \omega_{0}a^{\dagger}a+\frac{\omega_{1}}{2}{{\sigma}}_{z} +g(a^{\dagger}{{\sigma}}_-+a{{\sigma}}_+), $

      即为JC模型.

      如果考虑在晶格中每个格点上都有一个Rabi模型, 则得到Rabi晶格模型, 其哈密顿量为

      $ \begin{split} {\cal{H}} =\;& J\mathop \sum \limits_{\langle i,j\rangle } a_i^{\dagger}{a_j} + \mathop \sum \limits_i \Big[{\omega _0}a_i^{\dagger}{a_i} + \frac{{{\omega _1}}}{2}{{{\sigma}} _z} \\ &+ g(a_i^{\dagger} + {a_i})({{{\sigma}} _ - } + {{{\sigma}} _ + })\Big],\end{split}$

      在这里, 第一项是格点与格点之间的光子跃迁项, $ J $是任意两格点的跃迁能量.

      如果考虑在晶格中每个格点上的Rabi模型都有Kerr非线性项, 哈密顿量则为

      $ \begin{split} {{\cal{H}}_{{\rm{Kerr}}}} =\;& J\mathop \sum \limits_{\langle i,j\rangle } a_i^{\dagger}{a_j} + \mathop \sum \limits_i \Big[{\omega _0}a_i^{\dagger}{a_i} + \frac{{{\omega _1}}}{2}{{{\sigma}} _z} \\ & + g(a_i^{\dagger} + {{\hat a}_i})({{{\sigma}} _ - } + {{{\sigma}} _ + }) + \kappa a_i^{\dagger}a_i^{\dagger}{a_i}{a_i}\Big]\\ = \;&J\mathop \sum \limits_{\langle i,j\rangle } a_i^{\dagger}{a_j} + \mathop \sum \limits_i H_{{\rm{Rabi}}}^i + \kappa \mathop \sum \limits_i a_i^{\dagger}a_i^{\dagger}{a_i}{a_i}, \end{split}$

      这里的$ \kappa $是Kerr非线性强度, $ H_{\rm{Rabi}}^{i} $代表第$ i $个格点上的Rabi模型的哈密顿量.

    • Greentree等[2]已经用平均场理论讨论过JC晶格模型, 他们揭示了从极化子的超流体相到Mott绝缘相的二阶相变. 这里用同样的方法研究Rabi晶格模型.

      引入一个超流序参量

      $ A = \braket{a}, $

      来描述从Mott绝缘体相到超流体相的相变, $ \braket{\cdots} $表示对物理量求基态的期望值, 将跃迁项中的$ a_{i}^{\dagger} $代换为$ (a_{i}^{\dagger}-\braket{a_{i}^{\dagger}})+\braket{a_{i}^{\dagger}} $, $ a_{j} $ 代换为$ (a_{j}-\braket{a_{j}})+\braket{a_{j}} $:

      $ \begin{split} {\cal{H}} = \;&J \sum\limits_{\langle i,j\rangle }[(a_{i}^{\dagger}-\langle{a_{i}^{\dagger}}\rangle)+\langle{a_{i}^{\dagger}}\rangle][(a_{j}-\langle{a_{j}}\rangle)+\langle{a_{j}}\rangle]\\ & +\sum\limits_{i}H_{\rm{Rabi}}^{i}\\ = \;&J\sum\limits_{\langle i,j\rangle }[(a_{i}^{\dagger}-\langle{a_{i}^{\dagger}}\rangle)\langle{a_{j}}\rangle +(a_{j}-\langle{a_{j}}\rangle)\langle{a_{i}^{\dagger}}\rangle\\ &+\langle{a_{i}^{\dagger}}\rangle\langle{a_{j}}\rangle +(a_{i}^{\dagger}-\langle{a_{i}^{\dagger}}\rangle)(a_{j}-\langle{a_{j}}\rangle)]+\sum\limits_{i}H_{\rm{Rabi}}^{i}\\ \thickapprox \;& J\sum\limits_{\langle i,j\rangle }[A(a_{i}^{\dagger}+a_{j})-A^{2}]+\sum\limits_{i}H_{\rm{Rabi}}^{i}, \\[-17pt] \end{split} $

      其中$ (a_{i}^{\dagger}-\braket{a_{i}^{\dagger}})(a_{j}-\braket{a_{j}}) $为高阶涨落项, 在平均场近似下可以忽略.

      在一维晶格链中, 一个格点可以向相邻的两个格点跃迁, 所以, 有效哈密顿量可以写成:

      $ \begin{split} {\cal{H}}_{\rm{eff}} = \;&2JA\underset{i}{\sum}(a_{i}^{\dagger}+a_{i})+\underset{i}{\sum}\Big[\omega_{0}a_{i}^{\dagger}a_{i}+\frac{\omega_{1}}{2}{{\sigma}}_{z}\\ &+g(a_{i}^{\dagger}+a_{i})({{\sigma}}_-+{{\sigma}}_+)\Big]+{\rm {const}}. \end{split} $

      考虑矩阵直和, 有效哈密顿量可以分解为 ${\cal{H}}_{\rm{eff}} = $$ {\displaystyle\sum\nolimits_i}{\cal{H}}_{\rm{eff}}^{i}$. 其中, 第$ i $个格点的有效哈密顿分量为

      $ \begin{split} {\cal{H}}_{\rm{eff}}^{i} =\;& 2JA(a_{i}^{\dagger}+a_{i})+g({{\sigma}}_- +{{\sigma}}_+)(a_{i}^{\dagger}+a_{i}) \\ & +\omega_{0}a_{i}^{\dagger}a_{i} +\frac{\omega_{1}}{2}{{\sigma}}_{z}. \end{split} $

      如果考虑Kerr非线性项, 则有效哈密顿量为

      $ \begin{split} {\cal{H}}_{\rm{Kerr,eff}}^{i} =\;& 2JA(a_{i}^{\dagger}+a_{i}) +g({{\sigma}}_-+{{\sigma}}_+)(a_{i}^{\dagger}+a_{i}) \\ &+\omega_{0}a_{i}^{\dagger}a_{i}+\frac{\omega_{1}}{2}{{\sigma}}_{z} +{\kappa}a_{i}^{\dagger}a_{i}^{\dagger}a_{i}a_{i}, \end{split} $

      这里的$ \kappa $是Kerr非线性强度.

      由此, 可以根据有效哈密顿量写出对应的哈密顿量矩阵, 从而数值计算出相图.

    • 光的量子相变是指当控制光子之间的有效相互作用的参数发生变化时, 腔中光子的量子态发生质的改变, 所以本文将通过不同的物理变量来观察光的量子相变. 在数值计算上主要采用迭代法求得超流序参量$ A = \braket{a} $. 先为$ A $设定一个初始值, 得到一个有效哈密顿量, 通过求解相应的薛定谔方程求出定态波函数, 并由此计算得出下一个$ A = \braket{a} $, 如此迭代直到得到一个收敛的$ A $值. 利用超流序参量$ A = \braket{a} $、腔内光子数的基态期望值$ N = \braket{a^{\dagger}a} $和二阶关联函数$ g^{2}(0) $就可以描述Rabi晶格模型的量子相变.

      首先研究两种晶格模型基态下的光量子相变行为. 图1(a)图1(b)分别给出了JC晶格模型和Rabi晶格模型的$ J \text {-}g $相图, 即Mott绝缘体相-超流体相转变的相图, 其中$ g $为二能级原子和光子之间的相互作用强度, $ J $为格点之间的光子跃迁强度. 图中深蓝色区域$ A = 0 $, 表示Mott绝缘体相, 而$ A > 0 $的区域为超流体相. 根据哈密顿量, JC晶格模型相图(图1(a))中特殊的平台型图像(深蓝色区域)对应着不同的单格点内的总激发数, 并且此时的超流序参数$ A = 0 $. 而在Rabi晶格模型中, 由于考虑了反旋转波项, 平均场近似后的单格点约化系统的总激发数即使在Mott相中也不再守恒. 因此, 为了更加直观地讨论$ J $$ g $对量子相变的作用, 在图1(c)图1(d)中对JC晶格模型和Rabi晶格模型的相变进行了分析. 图1(c)为对于不同的$ J $, JC晶格模型的超流序参量$ A $$ g $变化的图像, 在$ J $很小时, 超流序参量$ A $不会随$ g $的增大而增大. 但当$ J $$ g $增大达到一定值时, $ A $会随$ g $的增大而增大. 图1(d)为对于不同的$ J $, Rabi晶格模型的超流序参量$ A $$ g $变化的图像, 可以看出, 当$ g $达到一定值时, $ A $$ g $的增大而增大. 图1(c)图1(d)说明光子跃迁强度$ J $对超流体相的产生具有促进作用. 对于JC晶格模型, 在光子跃迁强度$ J $较大的时候, 相互作用强度$ g $才会对超流体相的产生有一定的促进作用; 对于Rabi晶格模型, 由于反旋波项的存在, 相互作用强度$ g $会促进Mott绝缘相到超流体相的转变.

      图  1  (a), (b)平均场近似下, 不同晶格模型关于超流序参量$ A = \braket{a} $$ J \text {-} g $相图 (a) JC晶格模型; (b) Rabi 晶格模型. 横坐标为格点之间的光子跃迁强度$ J $, 纵坐标为二能级原子和光子相互作用强度$ g $, 横纵坐标的单位为$ \omega_0 $, 颜色条表示超流序参量$ A = \braket{a} $的大小. 深蓝色表示Mott绝缘相, 浅黄色表示超流体相. 其他参量取值为: $ \mathop{\omega_{0} = \omega_{1}} = 1 $, 光子截断数$ N = 20 $. (c), (d)对于不同的$ J $, 不同晶格模型的超流序参量$ A $$ g $变化的图像 (c) JC晶格模型; (d) Rabi晶格模型

      Figure 1.  (a), (b) Under the mean field approximation, the $ J \text {-} g $ phase diagram of different lattice models with respect to the superfluid order parameter $ A = \braket{a} $: (a) JC lattice model; (b) Rabi lattice model. The abscissa is the photon transition intensity $ J $ between the lattice, the ordinate is the two-level atom and photon interaction strength $ g $, the unit of the abscissa and the ordinate is $ \omega_0 $, and the color bar represents the value of the superfluid order parameter $ A = \braket{a} $. Dark blue indicates Mott insulating phase, and light yellow indicates superfluid phase. Other parameters are taken as $ \mathop{\omega_{0} = \omega_{1}} = 1 $, and the number of the photon truncation $ N = 20 $. (c), (d) For different $ J $, the superfluid order parameter $ A $ of different lattice models varies with $ g $: (c) JC lattice model; (d) Rabi lattice model.

      然后, 研究两类晶格体系中的光子统计行为, 利用Glauber[44]于1963年在量子光学的研究中提出的无延时两光子关联函数:

      $ g^{2}(0) = \frac{\braket{a^{\dagger}a^{\dagger}aa}}{\braket{a^{\dagger}a}^{2}}, $

      它表示测量第一个光子时测量到第二个光子的相对概率, 描述了两个连续光子测量之间的相关性. 在量子光学中, 光子聚束和反聚束表征了腔内光子的两种可区分状态. 光子聚束表示更多的光子被检测到是相互靠近, 而不是相互远离, 所以光子有呈束状分布的趋势, 这时$ g^{2}(0) > 1 $. 光子反聚束表示更多的光子单独存在而不是相互靠近, 这时光子没有呈束状分布的趋势, 这意味着要检测到第二个光子具有反相关效果, 即$ g^{2}(0) < 1 $. 对于相干态, 两次连续光子检测之间没有相关性, 光子呈随机分布, 此时$ g^{2}(0) = 1 $.

      图2为二阶关联函数$ g^{2}(0) $$J\text {-} g$相图, $ g^{2}(0) $$ J $$ g $的增大而增大. 黑线为$ g^{2}(0) = 1 $, $ g^{2}(0) > 1 $对应于光子聚束相, $ g^{2}(0) < 1 $对应于光子反聚束相. 可以看到, JC晶格模型没有光子聚束行为, 一直处于反聚束状态. 而Rabi晶格模型在相互作用强度$ g $较小的时候处于光子聚束, 随着$ g $的增大, Rabi晶格模型逐渐从光子聚束到反聚束变化. 同时, 随着光子跃迁强度$ J $的增大, Rabi晶格模型也逐渐从光子聚束到反聚束变化. 这说明相互作用强度$ g $和光子跃迁强度$ J $都会影响Rabi晶格模型光子聚束向反聚束的转变, 其性质与JC晶格模型有显著差异.

      图  2  平均场近似下, 不同晶格模型关于二阶关联函数$ g^{2}(0) $$J\text {-} g$相图 (a) JC晶格模型; (b) Rabi晶格模型. 横坐标为格点之间的光子跃迁强度$ J $, 纵坐标为二能级原子和光子相互作用强度$ g $, 横纵坐标的单位为$ \omega_0 $, 颜色条表示二阶关联函数$ g^{2}(0) $ 的值. 其他参量取值为: $ \mathop{\omega_{0} = \omega_{1}} = 1 $, 光子截断数$ N = 20 $

      Figure 2.  Under the mean field approximation, the $ J \text {-} g $ phase diagram of different lattice models with respect to the second-order correlation function $ g^{2}(0) $: (a) JC lattice model; (b) Rabi lattice model. The abscissa is the photon transition intensity $ J $ between the lattice, the ordinate is the two-level atom and photon interaction strength $ g $, the unit of the abscissa and the ordinate is $ \omega_0 $, the color bar is represented by the value of the second-order correlation function $ g^{2}(0) $. $ \mathop{\omega_{0} = \omega_{1}} = 1 $, and the number of photon truncation $ N = 20 $.

      接下来, 研究Kerr非线性项对Rabi晶格模型和JC晶格模型的量子相变以及对聚束和反聚束转变行为的影响.

      图3(a)图3(b)分别为JC晶格模型和Rabi晶格模型的Mott绝缘体相-超流体相相图, 刻画了二阶关联函数随二能级原子和光子之间的相互作用强度$ g $和Kerr非线性强度$ \kappa $的变化情况. 图3(a)中的Mott绝缘相被分为三个部分, 这与图1(a)中JC晶格模型的平台型结构有关, 可以验证对于不同的$ J $, Mott绝缘相可以分成不同的区域. 从图3(a)可以看出, 序参量$ A $随着Kerr非线性强度$ \kappa $的增大而减小, 即Kerr非线性项对超流体相有抑制作用. 这是因为JC晶格模型中非线性项增加了每个格点中光子的排斥, 从而导致系统趋于Mott相[2]. 图3(a)中淡蓝色直线区域与图1(a)$ N = 0 $$ N = 1 $之间的超流区域对应, 即在跃迁强度$ J $较小的时候, Kerr非线性项对JC晶格在$ g = 1 $$ N = 0 $$ N = 1 $之间的超流区域没有影响, 这与文献[45]的结果一致.

      图  3  Kerr效应下不同晶格模型关于超流序参量$ A = \braket{a} $$ \kappa \text {-} g $相图 (a) JC晶格模型; (b) Rabi 晶格模型. 横坐标为Kerr非线性强度$ \kappa $, 纵坐标为二能级原子和光子相互作用强度$ g $, 横纵坐标的单位为$ \omega_0 $, 颜色条表示超流序参量$ A $的大小. 其他参量取值为: $ \mathop{\omega_{0} = \omega_{1}} = 1 $, $ J = 0.05 $, 光子截断数$ N = 20 $

      Figure 3.  The $ \kappa \text {-} g $ phase diagram of different lattice models under the Kerr effect with respect to the superfluid order parameter $ A = \braket{a} $: (a) JC lattice model; (b) Rabi lattice model. The abscissa is the Kerr nonlinear intensity $ \kappa $, the ordinate is the two-level atom and photon interaction strength $ g $, the unit of the abscissa and the ordinate is $ \omega_0 $, and the color bar represents the value of the superfluid order parameter $ A $. Other parameters are taken as $ \mathop{\omega_{0} = \omega_{1}} = 1 $, $ J = 0.05 $, and the number of photon truncation $ N = 20 $.

      图4(a)4(b)中, 研究了 JC晶格模型没有聚束到反聚束的转变, 以及Rabi晶格模型随$ g $$ \kappa $的变化从聚束到反聚束的转变. 可以看出, 二阶关联函数$ g^{2}(0) $随着Kerr非线性强度$ \kappa $的增大而减小, 说明Kerr非线性强度$ \kappa $会促进光子反聚束的产生.

      图  4  Kerr效应下不同晶格模型关于二阶关联函数$ g^{2}(0) $$ \kappa \text {-} g $相图 (a) JC晶格模型; (b) Rabi晶格模型. 横坐标为Kerr非线性强度$ \kappa $, 纵坐标为二能级原子和光子相互作用强度$ g $, 横纵坐标的单位为$ \omega_0 $, 颜色条表示二阶关联函数$ g^{2}(0) $. 其他参量取值为: $ \mathop{\omega_{0} = \omega_{1}} = 1 $, $ J = 0.05 $, 光子截断数$ N = 20 $

      Figure 4.  The $ \kappa \text {-} g $ phase diagram of different lattice models under the Kerr effect with respect to the second-order correlation function $ g^{2}(0) $: (a) JC lattice model; (b) Rabi lattice model. The abscissa is the Kerr nonlinear intensity $ \kappa $, the ordinate is the two-level atom and photon interaction strength $ g $, the unit of the abscissa and the ordinate is $ \omega_0 $, and the color bar represents the value of second-order correlation function $ g^{2}(0) $. Other parameters are taken as $ \mathop{\omega_{0} = \omega_{1}} = 1 $, $ J = 0.05 $, and the number of photon truncation $ N = 20 $.

    • 本文采用平均场近似法研究了JC晶格模型和Rabi晶格模型, 得到了JC晶格模型和Rabi晶格模型的有效哈密顿量, 再用迭代法求得超流序参量, 然后对得到的物理图像进行分析, 从而研究了由二能级原子和光子耦合的晶格系统的性质.

      研究发现, JC晶格模型只存在Mott绝缘体相-超流体相量子相变, 而Rabi晶格模型存在Mott绝缘体相-超流体相的量子相变和光子聚束-反聚束的转变, 这是反旋波项影响的结果. 同时, 分别讨论了光子跃迁强度$ J $和相互作用强度$ g $对超流序参量$ A $的影响, 发现在Rabi晶格模型中, 二能级原子与光子之间的相互作用强度$ g $和格点之间的光子跃迁强度$ J $的增大, 会促进模型从Mott绝缘体相到超流体相的量子相变和光子聚束到反聚束的转变. 同时, 考虑了Kerr非线性项对JC晶格模型和Rabi晶格模型的影响, 发现在Rabi晶格模型中, Kerr非线性项抑制Mott绝缘体相向超流体相的转变, 但是促进光子聚束向反聚束转变, 这说明Kerr非线性项可以在光子之间产生有效的光子-光子排斥作用.

      这些研究结果凸显了二能级系统与光腔的相互作用强度$ g $和光子跃迁强度$ J $对JC晶格模型和Rabi晶格模型的重要性, 以及Kerr非线性项对Rabi晶格模型具有重要影响, 为研究这些模型的量子相变提供了理论方法, 对实验研究JC晶格模型和Rabi晶格模型有一定的指导意义.

参考文献 (45)

目录

    /

    返回文章
    返回