-
光与物质之间的作用是自然界中最基本的相互作用过程之一, 是理解多种物理现象的基础. 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晶格模型有一定的指导意义.
-
采用平均场近似的方法, 分别研究了Jaynes-Cummings晶格模型和Rabi晶格模型的量子相变: Mott绝缘体相-超流体相量子相变, 探索了光的聚束-反聚束行为, 研究了Kerr非线性作用对量子相变与光子统计特征的影响. 研究结果表明, 在Rabi晶格模型中二能级原子和光子相互作用强度g和格点之间光子跃迁强度J的增大会使晶格体系从Mott绝缘体相向超流体相转变, 同时, 光子统计行为由聚束转变为反聚束, 而Kerr非线性强度的增大抑制了Mott绝缘体相-超流体相相变, 但促进了光子聚束与反聚束之间的转变.
-
关键词:
- Jaynes-Cummings晶格模型 /
- Rabi晶格模型 /
- 量子相变 /
- 光子统计
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
[1] Hartmann M J, Brandao F G S L, Plenio M B 2006 Nat. Phys. 2 849
Google Scholar
[2] Greentree A D, Tahan C, Cole J H, Hollenberg L C 2006 Nat. Phys. 2 856
Google Scholar
[3] Angelakis D G, Santos M F, Bose S 2007 Phys. Rev. A 76 031805
Google Scholar
[4] Rossini D, Fazio R 2007 Phys. Rev. Lett. 99 186401
Google Scholar
[5] Aichhorn M, Hohenadler M, Tahan C, Littlewood P B 2008 Phys. Rev. Lett. 100 216401
Google Scholar
[6] Na N, Utsunomiya S, Tian L, Yamamoto Y 2008 Phys. Rev. A 77 031803
Google Scholar
[7] Carusotto I, Gerace D, Türeci H E, De Liberato S, Ciuti C, and Imamoğlu A 2009 Phys. Rev. Lett. 103 033601
Google Scholar
[8] Schmidt S, Blatter G 2009 Phys. Rev. Lett. 103 086403
Google Scholar
[9] Koch J, Le Hur K 2009 Phys. Rev. A 80 023811
Google Scholar
[10] Pippan P, Evertz H G, Hohenadler M 2009 Phys. Rev. A 80 033612
Google Scholar
[11] Ferretti S, Andreani L C, Türeci H E, Gerace D 2010 Phys. Rev. A 82 013841
Google Scholar
[12] Umucalilar R O, Carusotto I 2011 Phys. Rev. A 84 043804
Google 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 R6801
Google Scholar
[15] Imamoḡlu A, Schmidt H, Woods G, Deutsch M 1997 Phys. Rev. Lett. 79 1467
Google Scholar
[16] Rebic S, Tan S M, Parkins A S, Walls D F 1999 J. Opt. B 1 490
Google Scholar
[17] Bloch I, Dalibard J, Zwerger W 2008 Rev. Mod. Phys. 80 885
Google Scholar
[18] Schmidt S, Koch J 2013 Ann. Phys. 525 395
Google 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 546
Google Scholar
[22] van Oosten D, van Der Straten P, Stoof H T C 2001 Phys. Rev. A 63 053601
Google Scholar
[23] van Oosten D, van Der Straten P, Stoof H T C 2003 Phys. Rev. A 67 033606
Google Scholar
[24] Sheshadri K, Krishnamurthy H R, Pandit R, Ramakrishnan T V 1993 Europhys. Lett. 22 257
Google Scholar
[25] Krauth W, Trivedi N 1991 Europhys. Lett. 14 627
Google Scholar
[26] Krauth W, Trivedi N, Ceperley D 1991 Phys. Rev. Lett. 67 2307
Google Scholar
[27] Xie Z W, Liu W M 2004 Phys. Rev. A 70 045602
Google Scholar
[28] Albus A, Illuminati F, Eisert J 2003 Phys. Rev. A 68 023606
Google Scholar
[29] Lewenstein M, Santos L, Baranov M A, Fehrmann H 2004 Phys. Rev. Lett. 92 050401
Google Scholar
[30] Illuminati F, Albus A 2004 Phys. Rev. Lett. 93 090406
Google Scholar
[31] Cramer M, Eisert J, Illuminati F 2004 Phys. Rev. Lett. 93 190405
Google 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 S3597
Google Scholar
[34] He Y, Zhu X, Mihalache D, Liu J, Chen Z 2012 Phys. Rev. A 85 013831
Google Scholar
[35] Eguchi K, Takagi Y, Nakagawa T, Yokoyama T 2012 Phys. Rev. B 85 174415
Google Scholar
[36] Vitali D, Fortunato M, Tombesi P 2000 Phys. Rev. Lett. 85 445
Google Scholar
[37] Angelakis D G, Dai L, Kwek L C 2010 Europhys. Lett. 91 10003
Google Scholar
[38] Patargias N, Bartzis V, Jannussis A 1995 Phys. Scr. 52 554
Google Scholar
[39] Bu S P, Zhang G F, Liu J, Chen Z Y 2008 Phys. Scr. 78 065008
Google Scholar
[40] Cordero S, Récamier J 2011 J. Phys. B: At. Mol. Opt. Phys. 44 135502
Google Scholar
[41] Schmidt H, Imamoğlu A 1996 Opt. Lett. 21 1936
Google Scholar
[42] Harris S E, Hau L V 1999 Phys. Rev. Lett. 82 4611
Google Scholar
[43] Niu Y, Gong S 2006 Phys. Rev. A 73 053811
Google Scholar
[44] Glauber R J 1963 Phys. Rev. 130 2529
Google 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 849
Google Scholar
[2] Greentree A D, Tahan C, Cole J H, Hollenberg L C 2006 Nat. Phys. 2 856
Google Scholar
[3] Angelakis D G, Santos M F, Bose S 2007 Phys. Rev. A 76 031805
Google Scholar
[4] Rossini D, Fazio R 2007 Phys. Rev. Lett. 99 186401
Google Scholar
[5] Aichhorn M, Hohenadler M, Tahan C, Littlewood P B 2008 Phys. Rev. Lett. 100 216401
Google Scholar
[6] Na N, Utsunomiya S, Tian L, Yamamoto Y 2008 Phys. Rev. A 77 031803
Google Scholar
[7] Carusotto I, Gerace D, Türeci H E, De Liberato S, Ciuti C, and Imamoğlu A 2009 Phys. Rev. Lett. 103 033601
Google Scholar
[8] Schmidt S, Blatter G 2009 Phys. Rev. Lett. 103 086403
Google Scholar
[9] Koch J, Le Hur K 2009 Phys. Rev. A 80 023811
Google Scholar
[10] Pippan P, Evertz H G, Hohenadler M 2009 Phys. Rev. A 80 033612
Google Scholar
[11] Ferretti S, Andreani L C, Türeci H E, Gerace D 2010 Phys. Rev. A 82 013841
Google Scholar
[12] Umucalilar R O, Carusotto I 2011 Phys. Rev. A 84 043804
Google 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 R6801
Google Scholar
[15] Imamoḡlu A, Schmidt H, Woods G, Deutsch M 1997 Phys. Rev. Lett. 79 1467
Google Scholar
[16] Rebic S, Tan S M, Parkins A S, Walls D F 1999 J. Opt. B 1 490
Google Scholar
[17] Bloch I, Dalibard J, Zwerger W 2008 Rev. Mod. Phys. 80 885
Google Scholar
[18] Schmidt S, Koch J 2013 Ann. Phys. 525 395
Google 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 546
Google Scholar
[22] van Oosten D, van Der Straten P, Stoof H T C 2001 Phys. Rev. A 63 053601
Google Scholar
[23] van Oosten D, van Der Straten P, Stoof H T C 2003 Phys. Rev. A 67 033606
Google Scholar
[24] Sheshadri K, Krishnamurthy H R, Pandit R, Ramakrishnan T V 1993 Europhys. Lett. 22 257
Google Scholar
[25] Krauth W, Trivedi N 1991 Europhys. Lett. 14 627
Google Scholar
[26] Krauth W, Trivedi N, Ceperley D 1991 Phys. Rev. Lett. 67 2307
Google Scholar
[27] Xie Z W, Liu W M 2004 Phys. Rev. A 70 045602
Google Scholar
[28] Albus A, Illuminati F, Eisert J 2003 Phys. Rev. A 68 023606
Google Scholar
[29] Lewenstein M, Santos L, Baranov M A, Fehrmann H 2004 Phys. Rev. Lett. 92 050401
Google Scholar
[30] Illuminati F, Albus A 2004 Phys. Rev. Lett. 93 090406
Google Scholar
[31] Cramer M, Eisert J, Illuminati F 2004 Phys. Rev. Lett. 93 190405
Google 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 S3597
Google Scholar
[34] He Y, Zhu X, Mihalache D, Liu J, Chen Z 2012 Phys. Rev. A 85 013831
Google Scholar
[35] Eguchi K, Takagi Y, Nakagawa T, Yokoyama T 2012 Phys. Rev. B 85 174415
Google Scholar
[36] Vitali D, Fortunato M, Tombesi P 2000 Phys. Rev. Lett. 85 445
Google Scholar
[37] Angelakis D G, Dai L, Kwek L C 2010 Europhys. Lett. 91 10003
Google Scholar
[38] Patargias N, Bartzis V, Jannussis A 1995 Phys. Scr. 52 554
Google Scholar
[39] Bu S P, Zhang G F, Liu J, Chen Z Y 2008 Phys. Scr. 78 065008
Google Scholar
[40] Cordero S, Récamier J 2011 J. Phys. B: At. Mol. Opt. Phys. 44 135502
Google Scholar
[41] Schmidt H, Imamoğlu A 1996 Opt. Lett. 21 1936
Google Scholar
[42] Harris S E, Hau L V 1999 Phys. Rev. Lett. 82 4611
Google Scholar
[43] Niu Y, Gong S 2006 Phys. Rev. A 73 053811
Google Scholar
[44] Glauber R J 1963 Phys. Rev. 130 2529
Google Scholar
[45] Gomes C B C, Almeida F A G, Souza A M C 2016 Phys. Lett. A 38 1799
计量
- 文章访问数: 3954
- PDF下载量: 269
- 被引次数: 0