Fund Project:Project supported by the National Natural Science Foundation of China (Grant Nos. 11774316, 11704093).
Received Date:06 December 2020
Accepted Date:04 January 2021
Available Online:06 May 2021
Published Online:20 May 2021
Abstract: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. Keywords:Jaynes-Cummings lattice model/ Rabi lattice model/ quantum phase transition/ photon statistics
4.数值结果及讨论光的量子相变是指当控制光子之间的有效相互作用的参数发生变化时, 腔中光子的量子态发生质的改变, 所以本文将通过不同的物理变量来观察光的量子相变. 在数值计算上主要采用迭代法求得超流序参量$ 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晶格模型 Figure1. (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.
它表示测量第一个光子时测量到第二个光子的相对概率, 描述了两个连续光子测量之间的相关性. 在量子光学中, 光子聚束和反聚束表征了腔内光子的两种可区分状态. 光子聚束表示更多的光子被检测到是相互靠近, 而不是相互远离, 所以光子有呈束状分布的趋势, 这时$ 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 $ Figure2. 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 $ Figure3. 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 $ Figure4. 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 $.
5.结 论本文采用平均场近似法研究了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晶格模型有一定的指导意义.