摘要:量子光学中的Rabi模型描述了一个与量子谐振子耦合的两能级系统, 当耦合强度与振子频率相当时, 会产生丰富的物理现象. 本文研究了谐波势阱中具有拉曼诱导自旋轨道耦合的Bose气体, 通过将受限系统映射为Rabi模型, 引入量子光学中的平移Fock态利用变分方法求解了系统的本征能态和基态, 发现左右平移Fock态的奇宇称叠加态能量低于平移态的能量，并分别研究了粒子在动量和坐标空间的动力学Zitterbewegung振荡特性以及原子极化的动力学, 给出了一个直观清晰的物理图像, 与相关实验的结果定性一致. 本文的研究结果有助于进一步研究量子光学领域目前难以实现的深度强耦合参数区域的量子Rabi模型, 对冷原子物理的研究也有一些借鉴和启发.
关键词: Rabi模型/
自旋轨道耦合/
Zitterbewegung振荡/
平移Fock态/
变分法

English Abstract

全文HTML

--> --> -->

自旋轨道耦合(spin-orbit coupling, SOC)是粒子自旋内禀自由度与其外部运动自由度之间的相互作用, 在凝聚态物理许多重要的现象中扮演着重要的角色^[1—3]. 近几年来, 冷原子物理学的一个重要进展是实现了光诱导合成规范场的中性原子的SOC^[4—7], 目前在实验上可以在玻色子和费米子超冷原子气体中实现各种SOC的哈密顿量^[8—12]. 虽然大多理论工作都集中在Rashba耦合的SOC的均匀系统上, 但是谐振子束缚势下具有SOC的冷原子系统也受到关注和研究^[13—17].
量子Rabi模型^[18,19]是量子光学中重要的基础模型之一, 模型形式简单, 但精确求解并不容易. 直到2011年, Braak^[20]提出Rabi模型具有Z₂对称性保证了其可积性, 因而可以通过求解超越方程得到系统的精确能谱. 该工作重新激起了人们对量子Rabi模型研究的兴趣^[21—24]. Rabi模型被广泛应用于不同的物理领域, 包括量子光学^[25]、量子信息^[26]、凝聚态物理^[27]和腔量子电动力学^[28—31], 也为研究受限空间中的冷原子系统提供了模型和方法.
在量子光学中考虑原子的动能大大超过了相互作用能, 通常采用绝热近似不考虑原子质心动能. 然而, 利用激光冷却原子运动降低原子动能可以制备出1 K量级的冷原子气体, 这个能量与腔量子电动力学实验的相互作用能相当, 从哈密顿量中排除动能项不再合理. 本文研究一维(1D)谐波势阱中具有SOC的Bose气体(考虑稀薄原子气体, 忽略相互作用), 采用量子光学中的方法求解系统的本征能态及可观测物理量的动力学演化, 与目前相关的实验结果定性一致^[11,12].

三维空间谐振子束缚势中具有SOC的单原子Bose气体的哈密顿量形式为

${\hat H_{3{\rm{d}}}} = \frac{{{{\hat P}^2}}}{{2m}} + {\hat V_{{\rm{SO}}}} + \hat V(\hat r),$

其中$\hat V\left( {\hat r} \right) = m\left( {\omega _x^2{x^{\rm{2}}} + \omega _y^2y{}^{\rm{2}} + \omega _z^2{z^{\rm{2}}}} \right){\rm{/2}}$是束缚势, m为原子质量, $\hat P$为原子动量; ${\hat V_{{\rm{SO}}}}$为SOC项, 对于x方向具有SOC两组分的Bose气体(视为赝自旋1/2)可以写为^[4,30]

${\hat V_{{\rm{SO}}}} = 2\hbar {\kappa _r}{\hat p_x}{\hat \sigma _z} + \frac{{\hbar \varOmega }}{2}{\hat \sigma _x} + \frac{{\hbar \delta }}{2}{\hat \sigma _z},$

其中${\kappa _r}$是双光子反冲动量, $\varOmega $是拉曼耦合强度, $\delta $是双光子失谐, $\hat \sigma $为泡利矩阵算符. 三维哈密顿量${\hat H_{{\rm{3d}}}} = \hat H(x) + {\hat H_{{\rm{2d}}}}(y,z)$可以通过维度塌缩方法变为一维系统, 这样一维系统的哈密顿量为

$\hat H = \frac{{\hat p_x^2}}{{2m}} + \frac{1}{2}m{\omega ^2}{x^2} + {\hat V_{{\rm{SO}}}}.$

将${\omega _x}$简记为$\omega $. (3)式的模型看似简单, 除了特殊情况$\omega = 0$或$\varOmega = {\rm{0}}$之外, 该模型不容易解析求解.

为了将谐波势阱中的SOC模型变换到量子Rabi模型, 将原子质心运动量子化用${a^\dagger }(a)$产生(湮灭)算符来表示, 方程(3)改写为

$\hat H = \hbar \omega \left({a^\dagger }a + \frac{1}{2}\right) + \hbar \delta {\hat \sigma _z}/2 + \hbar \varOmega {\hat \sigma _x}/2 + {\rm{i}}\hbar \lambda ({a^\dagger } - a){\hat \sigma _z},$

其中${p_x} = {\rm{i}}\sqrt {m\hbar \omega /2} ({a^\dagger } - a)$和$\lambda = {k_r}\sqrt {2m\hbar \omega } $. 略去常数项哈密顿量简化为

$\hat H = \hbar \omega {a^\dagger }a + \hbar \delta {\hat \sigma _z}/2 + \hbar \varOmega {\hat \sigma _x}/2 + {\rm{i}}\hbar \lambda ({a^\dagger } - a){\hat \sigma _z}.$

做一旋转变换${a^\dagger } \to {a^\dagger }{{\rm{e}}^{{\rm{i}}\theta }},\;a \to a{{\rm{e}}^{ - {\rm{i}}\theta }}\left( {\theta = {\text{π}}/2} \right)$, 方程(5)变为标准的Rabi模型

${\hat H_{\rm{R}}} = \hbar \omega {a^\dagger }a + \hbar \delta {\hat \sigma _z}/2 + \hbar \varOmega {\hat \sigma _x}/2 + \hbar \lambda ({a^\dagger } + a){\hat \sigma _z}.$

量子光学中Rabi模型描述了量子化光场与二能级原子内态的耦合, 此处则描述了原子赝自旋(此后简称为量子比特)与质心动量的耦合. 这样就可以利用量子光学中的方法来求解该系统. 当$\delta = {\rm{0}}$时, 系统进一步简化为

$\hat H = \frac{{\hbar \varOmega }}{2}{\hat \sigma _x} + \hbar \omega {a^\dagger }a + \hbar \lambda {\hat \sigma _z}({a^\dagger } + a),$

其中$\varOmega $是原子内部能级的能量差, $\omega $为谐振子束缚势频率, $\lambda $为自旋-轨道耦合强度.

系统的哈密顿量(7)式在${\hat \sigma _z}$表象中可以表示为

$\begin{split}H = & \frac{{\hbar \varOmega }}{2}\left( {\left| { + z} \right\rangle \left\langle { - z} \right| \!+\! \left| { - z} \right\rangle \left\langle { + z} \right|} \right) \!+\! {H_{\rm{L}}}\left| { + z} \right\rangle \left\langle { + z} \right| \\ & + {H_{\rm{R}}}\left| { - z} \right\rangle \left\langle { - z} \right| - {{\hbar {\lambda ^2}}}/{\omega },\end{split}$

其中${H_{\rm{L}}}$和${H_{\rm{R}}}$分别是依赖于量子比特的振子左右平移哈密顿量^[32],

${H_{{\rm{L,R}}}} = \hbar \omega \left( {{a^\dagger } \pm \lambda {\rm{/}}\omega } \right)\left( {a \pm \lambda {\rm{/}}\omega } \right),$

$\left| { \pm z} \right\rangle $为${\hat \sigma _z}$的本征态. 如果忽略$\varOmega $项, 则能量本征态就是与量子比特相关的振子的平移Fock态($\left| { + z} \right\rangle $对应左平移态, $\left| { - z} \right\rangle $对应右平移态). 通常人们将$\varOmega $项视为微扰(绝热近似), 但本文主要研究$\varOmega {\rm{/}}\omega > 1$的情况. 下面我们采用试探波函数$\left| {{\psi _0}} \right\rangle $方法, 将振子平移量和量子比特内态的相干叠加参数分别作为变分参数:

$\left| {{\psi _0}} \right\rangle = \left| \alpha \right\rangle \otimes \left( {\cos \frac{\theta }{2}\left| { + x} \right\rangle + \sin \frac{\theta }{2}\left| { - x} \right\rangle } \right),$

其中$\left| \alpha \right\rangle $为相干态, $\left| { \pm x} \right\rangle $为${\hat \sigma _x}$的本征态. 在坐标和动量表象中表示为

$\begin{split}& \left\langle {q|\alpha } \right\rangle = {\left( {\dfrac{{m\omega }}{{{\text{π}}\hbar }}} \right)^{{\rm{1/4}}}}\exp \left[ { - {{\left( {\sqrt {\dfrac{{m\omega }}{{2\hbar }}} q - \alpha } \right)}^2}} \right],\\& \left\langle {p|\alpha } \right\rangle \!=\! {\left( {{\text{π}}m\hbar \omega } \right)^{ - {\rm{1/4}}}}\exp \bigg[ { - {{\bigg( {\sqrt {\dfrac{1}{{2m\hbar \omega }}} p \!+\! {\rm{i}}\alpha } \bigg)}^2} - {\alpha ^2}} \Bigg],\end{split}$

其中m是原子的质量. 假定$\alpha $为实数, 则能量期待值表示为

${E_0} \!=\! \left\langle {{\psi _0}} \right|H\left| {{\psi _0}} \right\rangle \! =\! \frac{{\hbar \varOmega }}{2}\cos \theta \!+\! 2\hbar \lambda \alpha \sin \theta \!+\! \hbar \omega {\alpha ^2}.$

故此, 可以通过对含参数$\alpha $和$\theta $的能量泛函最小化以求得基态能量. 根据二元函数极值判据可知当${{{\rm{4}}{\lambda ^{\rm{2}}}}/{\left( {\varOmega \omega } \right)}} < 1$时, 系统能量只有一个局域最小值, 在$\alpha = {\rm{0}}$和$\theta = {\rm{0}}$时取得; 但如果条件

$\frac{1}{\varepsilon } \equiv \frac{{4{\lambda ^2}}}{{\varOmega \omega }} > 1$

成立, 则能量有两个局域极小值, 分别由$\alpha = \pm {\alpha _0}$和$\theta = \pm {\theta _0}$给出, 且有

$ \cos {\theta _0} = - {{\varOmega \omega }}{{/(4{\lambda ^2)}}}, $

$ {\alpha _0}= - \frac{\lambda }{\omega }\sqrt {1 - {{\left( {\frac{{\varOmega \omega }}{{4{\lambda ^2}}}} \right)}^2}}.$

此后我们假设条件(13)始终成立, 另外我们假定选择${\theta _{\rm{0}}}$使得${\rm{sin}}{\theta _0} \geqslant 0$, 这使得${\alpha _{\rm{0}}} \leqslant {\rm{0}}$, 分别对应两个简并态$\left| {{\psi _{0,{\rm{L}}}}} \right\rangle = \left| {{\alpha _0}} \right\rangle \left( {{\rm{cos}}\dfrac{{{\theta _0}}}{2}\left| { + x} \right\rangle + {\rm{sin}}\dfrac{{{\theta _0}}}{2}\left| { - x} \right\rangle } \right)$(左平移)和$\left| {{\psi _{0,{\rm{R}}}}} \right\rangle = \left| { - {\alpha _0}} \right\rangle \left( {{\rm{cos}}\dfrac{{{\theta _0}}}{2}\left| { + x} \right\rangle - {\rm{sin}}\dfrac{{{\theta _0}}}{2}\left| { - x} \right\rangle } \right)$(右平移), 相应的能量为

$E\left( {{\alpha _{\rm{0}}},{\theta _{\rm{0}}}} \right) = - \frac{{\hbar {\lambda ^{\rm{2}}}}}{\omega } - \frac{{\hbar {\varOmega ^{\rm{2}}}\omega }}{{{\rm{16}}{\lambda ^{\rm{2}}}}} = - \frac{{\hbar {\lambda ^{\rm{2}}}}}{\omega }\left( {{\rm{1}} + {\varepsilon ^{\rm{2}}}} \right),$

此处$\varepsilon=\varOmega \omega/(4\lambda^2) $. 这两个简并态均不是宇称算符$\hat P = \exp [{\rm i}{\text{π}}({a^ + }a + {\sigma _x}/$$2 - 1/2)]$的本征态^[33], 通过相干叠加可以得到具有确定宇称的量子态

$\left| {{\varPhi _{ \pm {\rm{,0}}}}} \right\rangle = \frac{{\rm{1}}}{{\sqrt {\rm{2}} {\eta _{ \pm {\rm{,0}}}}}}\left( {\left| {{\psi _{{\rm{0,}}{\kern 1pt} {\rm{L}}}}} \right\rangle \pm \left| {{\psi _{0,{\rm{R}}}}} \right\rangle } \right),$

这里${\eta _{ \pm ,0}} = \sqrt {1 \mp \varepsilon {{\rm{e}}^{ - 2\alpha _0^2}}} $是归一化系数, $\left| {{\psi _{{\rm{0,}}\,{\rm{L}}}}} \right\rangle \,$和$\left| {{\psi _{0,{\rm{R}}}}} \right\rangle $在一般不正交(见(20), (22)式). 计算证明可得若(13)式成立, 就有$\left\langle {{\varPhi _{ - ,0}}} \right|H\left| {{\varPhi _{ - ,0}}} \right\rangle < E\left( {{\alpha _0},{\theta _0}} \right)$成立(见图1), 因此, $\left| {{\varPhi _{ - ,0}}} \right\rangle $是系统基态更好的近似. 为保证波函数具有确定的宇称, 定义
图 1 简并量子态$\left| {{\psi _{N{\rm{,L}}\left( {\rm{R}} \right)}}} \right\rangle $能量${E_{N,{\rm{L/R}}}}$与左右平移奇宇称叠加态$\left| {{\psi _{ - ,N}}} \right\rangle $能量${E_{ - ,N}}$随SO耦合强度$\lambda $的变化　可见$N = 0$叠加态$\left| {{\psi _{ - ,0}}} \right\rangle $能量最低, 更接近基态; 而对于激发态$N \ne 0$, 二者能量随参数变化出现交叉; 相关参数取值为$\varOmega \; = \;{\rm{1}}{\rm{.4}}\omega $, 与文献[19]精确解的结果基本一致
Figure1. The energies of degenerate quantum states $\left| {{\psi _{N{\rm{,L}}\left( {\rm{R}} \right)}}} \right\rangle $ and the superposition state of odd parity $\left| {{\psi _{ - ,N}}} \right\rangle $ of left(right)-displaced number states varies as the spin-orbit coupling strength $\lambda $. It is seen that for $N = 0$, the superposition state has the lowest energy which is the best approximation for the ground state in our interest. And for the cases of $N \ne 0$, the energies of the two quantum states have pitchforks.The relevant parameters is Ω=1.4 and the results are in agreement with those in Ref.[19].

$\left| {{\varPhi _{ \pm ,N}}} \right\rangle = \frac{1}{{\sqrt 2 {\eta _{ \pm ,N}}}}\left( {\left| {{\psi _{N,{\rm{L}}}}} \right\rangle \pm \left| {{\psi _{N,{\rm{R}}}}} \right\rangle } \right),$

这里左右平移态分别为:

$\begin{split}& \left| {{\psi _{N,{\rm{L}}}}} \right\rangle \!=\! \hat D\left[ {{\alpha _0}} \right]\left| N \right\rangle \otimes \left( {\cos \frac{{{\theta _0}}}{2}\left| { + x} \right\rangle + \sin \frac{{{\theta _0}}}{2}\left| { - x} \right\rangle } \right), \\ & \left| {{\psi _{N,{\rm{R}}}}} \right\rangle \!=\! \hat D\left[ { - {\alpha _0}} \right]\left| N \right\rangle \!\otimes\! \left( {\cos \frac{{{\theta _0}}}{2}\left| { + x} \right\rangle \!-\! \sin \frac{{{\theta _0}}}{2}\left| { - x} \right\rangle } \right).\end{split} $

归一化系数为

${\eta _{ \pm ,N}} \!=\! {\bigg[ {1 \!\mp\! \varepsilon {{\rm{e}}^{ - 2\left( {{{{\lambda ^2}}/{\omega _{}^2}}} \right)\left( {1 - {\varepsilon ^2}} \right)}}{L_N}\Big( {\frac{{4{\lambda ^2}}}{{\omega _{}^2}}\left( {1 \!-\! {\varepsilon ^2}} \right)} \Big)} \bigg]^{{1/2}}}.$

在(19)式中, $\hat D\left[ \alpha \right] \equiv \exp \left[ {\alpha \left( {{a^ + } - a} \right)} \right]$是平移算符, $\hat D\left[ { \pm \alpha } \right]\left| N \right\rangle $是平移Fock态, 简记为$ \left| { \pm\alpha ,N} \right\rangle $, 左右平移Fock态的内积可用Laguerre多项式给出:

$\left\langle { - \alpha ,N|\alpha ,N} \right\rangle = {{\rm{e}}^{ - 2{\alpha ^2}}}{L_N}(4{\alpha ^2}),$

平移Fock态$\left| {{\psi _{N,{\rm{L}}\left( {\rm{R}} \right)}}} \right\rangle $与$\left| {{\psi _{0,{\rm{L}}\left( {\rm{R}} \right)}}} \right\rangle $源于相同的变分计算, 对于不同声子数N均具有相同的${\alpha _0}$和${\theta _0}$, 并且对于不同N , $\left| {{\psi _{N,{\rm{L/R}}}}} \right\rangle $彼此正交, 但左右平移振子态并不正交; 对于相同的声子数且当N较大时, 内积以${\rm{1/}}\left( {{N^{{1/4}}}{\alpha ^{{1/2}}}} \right)$衰变. $\left| {{\psi _{N,{\rm{L}}}}} \right\rangle $和$\left| {{\psi _{N,{\rm{R}}}}} \right\rangle $在能量上简并, 考虑宇称自然构造出(18)式的叠加态. 容易看出, $\left| {{\varPhi _{ \pm ,N}}} \right\rangle $为一完备非正交的基. $\left| {{\varPhi _{ \pm ,N}}} \right\rangle $的宇称是$ \pm {\left( { - 1} \right)^N}$, 并且相反宇称态正交, 相同宇称态一般不正交; 对于大声子数N和M, 交叠积分$\left\langle { - \alpha ,N|\alpha ,M} \right\rangle $至少以${\alpha ^{ - 1/2}}$的速度衰减. 考虑(15)式, 在极限$\lambda \gg \omega $的参数区域可以得到近似的正交基. 此外交叠积分也为因子$\cos {\theta _0} = - \varepsilon $所抑制, 对于大的耦合常数λ积分也是小量. 这些条件允许我们用期望值近似作为H的能量本征值, 也就是哈密顿量在近似正交基(18)式下的对角元. 对于足够大的${\alpha _0}$非对角元很小. 近似能量${E_{ \pm ,N}} = \left\langle {{\varPhi _{ \pm ,N}}} \right|H\left| {{\varPhi _{ \pm ,N}}} \right\rangle $表示为

$\begin{split} {E_{ \pm ,N}} & \!=\! \frac{\hbar }{{1 \mp \varepsilon {{\rm{e}}^{ - 2\alpha _0^2}}{L_N}\left( {4\alpha _0^2} \right)}}\left[ { - \frac{\varOmega }{2}\varepsilon \!+\! N\omega \!-\! \frac{{{\lambda ^2}}}{\omega }\left( {1 \!-\! {\varepsilon ^2}} \right) \!\pm\! \left( {\frac{\varOmega }{2} \!-\! N\omega \varepsilon \!+\! \frac{{{\lambda ^2}}}{\omega }\varepsilon \left( {1 \!-\! {\varepsilon ^2}} \right)} \right){{\rm{e}}^{ - 2\alpha _0^2}}{L_N}\left( {4\alpha _0^2} \right)} \right] \\ & = \frac{\hbar }{{1 \mp \varepsilon {{\rm{e}}^{ - 2\alpha _0^2}}{L_N}\left( {4\alpha _0^2} \right)}}\left[ { - \frac{\varOmega }{2}\left( {\varepsilon \mp {{\rm{e}}^{ - 2\alpha _0^2}}{L_N}\left( {4\alpha _0^2} \right)} \right)} \right] + \hbar N\omega - \hbar \frac{{{\lambda ^2}}}{\omega }\left( {1 - {\varepsilon ^2}} \right).\end{split} $

这个近似结果解有助于我们直观地理解系统的动力学.

考虑系统的动力学演化, 对于稀薄原子气体可以忽略原子间相互作用. 我们取初态$\left| {{\psi _{{\rm{0,L}}}}} \right\rangle $, 在初始时刻时开启拉曼诱导的SOC, 则系统的初始波函数为

$\begin{split}&~~~~ \left| {\varPsi \left( {t = 0} \right)} \right\rangle = \left| {{\psi _{0,{\rm L}}}} \right\rangle \\ &= \left| {{\alpha _0}} \right\rangle \left( {\cos \frac{{{\theta _0}}}{2}\left| { + x} \right\rangle + \sin \frac{{{\theta _0}}}{2}\left| { - x} \right\rangle } \right) \\ &\simeq \frac{1}{{\sqrt 2 }}\left( {\left| {{\varPhi _{ - ,0}}} \right\rangle + \left| {{\varPhi _{ + ,0}}} \right\rangle } \right),\end{split}$

其时间演化近似为

$\left| {\varPsi \left( t \right)} \right\rangle = \frac{{{{\rm{e}}^{{{ - {\rm{i}}{E_{ - ,0}}t}/\hbar }}}}}{{\sqrt 2 }}\left( {\left| {{\varPhi _{ - ,0}}} \right\rangle + {{\rm{e}}^{ - {\rm{i}}\vartriangle \omega t}}\left| {{\varPhi _{ + ,0}}} \right\rangle } \right),$

其中$\Delta\omega $是频率差,

$\begin{split}\Delta \omega & = \frac{{{E_{ + ,0}} - {E_{ - ,0}}}}{\hbar } = \frac{{\varOmega \left( {1 - {\varepsilon ^2}} \right){{\rm{e}}^{ - 2\alpha _0^2}}}}{{1 - {\varepsilon ^2}{{\rm{e}}^{ - 4\alpha _0^2}}}} \\ &\simeq \varOmega \left( {1 - {\varepsilon ^2}} \right){{\rm{e}}^{ - 2\alpha _0^2}},\end{split}$

在(25)式中忽略了${{\rm{e}}^{ - 2\alpha _0^2}}$的高阶幂次. 可以看出在初始时刻$t = 0$时, 初态动量分布主要位于左侧(振子相干态$\left| {{\alpha _0}} \right\rangle $); 而在时刻$t = {{\text{π}}/{\Delta\omega }}$时, 动量分布主要位于右侧(振子相干态$\left| { - {\alpha _0}} \right\rangle $); 在时刻$t = {{\text{π}}/({2\Delta \omega })}$, 原子动量概率分布呈双峰分布, 对应于两个相干态的叠加, 这是标准的隧穿运动, 与经典双势阱完全类似.
图2和图3分别给出了粒子在动量空间和坐标空间概率分布的动态特性, 由(23)式的近似值计算得出, 在这里我们取$\varOmega \; = \;{\rm{3}}\omega $和$ \lambda \; = \;{\rm{2}}\omega $. 可以清楚地看到原子质心动量和空间位置分布的特征隧穿行为, 即所谓 Zitterbewegung振荡.
图 2 原子动量分布概率的粗粒动力学演化 (3D, 左侧; 2D, 右侧)　相关参数取值为$\varOmega \; = \;{\rm{3}}\omega $, $ \lambda \; = \;{\rm{2}}\omega $, 初态为$\varPsi (t = 0) = {\psi _{0,{\rm L}}}$, 动量$ \tilde p = \sqrt {1/m\hbar \omega } p $
Figure2. The coarse dynamics evulution of momentum distribution of single particle (left for 3D; right for 2D) with $\varOmega \; = \;{\rm{3}}\omega $ and $ \lambda \; = \;{\rm{2}}\omega $. The initial state is set as $\varPsi (t = 0) = {\psi _{0,{\rm{L}}}}$. Momentum $\tilde p$ is defined by $ \tilde p = \sqrt {1/m\hbar \omega } p $.

图 3 原子空间位置分布概率的粗粒动力学演化(3D, 左侧; 2D, 右侧)　 相关参数取值及初态同图2, 位置$ \tilde q = \sqrt {m\omega /\hbar } q $
Figure3. The coarse dynamics evolution of position distribution of single particle (left for 3D; right for 2D) with the same parameters and the initial state in Fig. 2 and $ \tilde q = \sqrt {m\omega /\hbar } q $.

另外两组分原子布居差${\sigma _z}$的期望值$\left\langle {{\sigma _z}} \right\rangle = \sin {\theta _0}\cos \left( {\Delta \omega t} \right)$描述了原子的极化率的动力学. 图4显示了原子极化$\left\langle {{\sigma _z}} \right\rangle $随时间的演化, 可以看到$\left\langle {{\sigma _z}} \right\rangle $在1和–1之间周期振荡.
图 4 原子极化$\left\langle {{\sigma _z}} \right\rangle $随时间演化初态为$\varPsi \left( {t = 0} \right) = {\psi _{0,{\rm{L}}}}$, 参数取值为$\varOmega \; = \;{\rm{3}}\omega $和$ \lambda \; = \;{\rm{2}}\omega $, 时间以因子${{2{\text{π}}}/{\Delta \omega }}$标度
Figure4. Time evolution of $\left\langle {{\sigma _z}} \right\rangle $ with the initial state being $\varPsi \left( {t = 0} \right) = {\psi _{0,{\rm{L}}}}$ and the parameters $\varOmega \; = \;{\rm{3}}\omega $ and $ \lambda \; = \;{\rm{2}}\omega $. The time is scaled by the tunneling period $2{\text{π}}/\Delta\omega $.

综上, 我们求解了谐波势阱中拉曼诱导自旋轨道耦合的Bose气体, 通过将系统完全映射到量子Rabi模型, 将辐射场变为声子场, 运用量子光学中平移Fock态的方法得到了强耦合区域谐波势阱中自旋轨道耦合的Bose气体模型的基态解及系统的动力学演化, 直观地给出了原子质心空间坐标和动量及原子极化随时间的振荡图像, 与相关的实验结果定性相符.
传统量子光学中二能级系统与振子系统的耦合强度受到很大限制, 但在本系统中原子自旋轨道耦合强度可以通过Raman耦合来调节, 冷原子的质心动能很小, 不能像量子光学中惯常采用绝热近似忽略掉原子的质心动能, 这使得本模型科学合理. 冷原子系统具有良好的可调控性, 通过改变束缚势阱的频率以及Raman激光的波长和强度, 可以实现量子光学中Rabi模型目前无法达到的参数区域—深度强耦合区域. 本文的研究也为自旋轨道耦合的冷原子系统提供了一个新的方法和视角.