1.   引　言

孤子作为一种非线性波, 因其独特的传播性质及潜在的应用价值, 已成为非线性科学研究领域的重要研究课题之一. 孤子也是自然界中的一种普遍的非线性现象, 并广泛地存在于各种非线性介质中, 如水波、等离子体、粒子物理、分子生物学及光纤等^[1]. 特别地, 随着玻色-爱因斯坦凝聚(BEC)和简并费米气体的实验实现, 大量的研究结果展示, 超冷原子气体中也存在物质波孤子现象, 实验上已经相继观察到了物质波亮孤子、暗孤子及涡旋孤子等非线性现象^[2-13]. 由于物质波孤子在相干原子光学、原子干涉仪及原子输运等领域中存在着潜在的应用价值, 研究超冷原子气体中的孤子动力学性质也成为了近几十年的热点研究课题之一.

近年来, 人造自旋-轨道耦合在超冷原子气体中的实验实现, 也为探索规范场中孤子的动力学性质提供了平台^[14-17]. 一方面, 自旋-轨道耦合使得体系单粒子基态在有限动量处简并^[18,19], 自旋-轨道耦合在BEC中将导致许多新奇的静态孤子, 例如条纹孤子和分数涡旋能隙孤子等^[20-42]. 另一方面, 在多组分BEC中, 孤子可看作是具有赝自旋的粒子, 自旋-轨道耦合将孤子自旋与质心耦合在一起, 使得孤子的自旋旋进将影响其质心运动. 例如, 孤子自旋周期性翻转将提供一个周期性的力去驱使孤子质心做周期性振荡^[43-46]. 由此可见, 自旋-轨道耦合为孤子的宏观量子调控提供了新的手段.

实验上在超冷原子气体中产生人造自旋-轨道耦合时, 不同分量间的能量差及Raman激光的频率差将产生一个有效的线性塞曼劈裂^[14-17], 它使得单粒子能谱的对称结构被破坏, 并导致许多新奇的量子态, 如BEC中的极化平面波态、超冷费米气体中的拓扑超流态和Majorana费米子等^[47-51]. 除此以外, 甚至在没有自旋-轨道耦合的旋量BEC中, 塞曼劈裂也会导致一些新奇的拓扑孤子态^[52]. 因此, 在具有自旋-轨道耦合的BEC中, 线性塞曼劈裂将对孤子的动力学性质产生明显的影响.

本文以一维自旋-轨道耦合双组份BEC为研究对象, 深入研究线性塞曼劈裂对亮孤子动力学性质的影响. 假设体系的原子相互作用守恒SU(2)对称性, 并取双曲正弦函数作为亮孤子的试探波函数, 本文首先利用变分法导出了试探波函数中的变分参数(未知参数)随时间演化所满足的欧拉-拉格朗日方程; 然后求解欧拉-拉格朗日方程的固定点解, 在自旋-轨道耦合强度较弱时, 发现了两个近似的静态亮孤子解; 进一步对这两个静态孤子做线性稳定性分析, 获得了一个零能的Goldstone激发模和一个谐振激发模, 前者对应于平移对称性的破缺, 后者的振荡频率与线性塞曼劈裂的强度有关; 最后, 通过求解欧拉-拉格朗日方程, 获得了变分参数的含时精确解, 并发现线性塞曼劈裂明显地影响孤子的运动速度和振荡周期. 这些变分计算结果与基于Gross-Pitaevskii (GP)方程的直接数值模拟结果相吻合.

2.   理论模型

考虑沿$z$方向的一维均匀自旋-轨道耦合双组分BEC^[14-17]. 在由同一种原子的两种不同超精细态所形成的双组分BEC的实验中, 由于两个分量的种内原子相互作用强度和种间原子相互作用强度通常比较接近, 本文将假设两者相等, 即原子间的相互作用守恒SU(2)对称性. 因此, 在平均场近似下, 系统的动力学性质可用如下的无量纲化GP方程描述：

方程中$\psi_{s}\left(z, t\right)$为描述两个分量的动力学性质的波函数($s = \uparrow, \downarrow$代表两个不同的分量), 且满足归一化条件$\displaystyle\sum_{s = \uparrow, \downarrow}\int_{-\infty}^{+\infty}\vert\psi_{s}\vert^2 {\rm d}z = 1$, $z$和$t$分别为空间坐标和时间; $n = \vert \psi_{\uparrow}\vert^2+\vert \psi_{\downarrow}\vert^2$为总密度分布; $\varepsilon$为线性塞曼劈裂的强度; $\varOmega$和$k_{\rm R}$分别表示Raman激光强度和自旋-轨道耦合强度; $g<0$代表原子间的吸引相互作用强度. 当自旋-轨道耦合存在时, GP方程((1a), (1b))不可积^[25], 本文利用变分近似方法解析研究孤子的动力学性质^[1]. 在变分近似方法中, 体系所对应的拉格朗日量为

其中$\ast$代表复共轭.

对于SU(2)自旋对称的均匀双组分BEC, 假设BEC凝聚在单个准动量为$k$的态上, 则可以利用双曲正弦函数作为亮孤子的试探波函数, 并取两个分量孤子的宽度$\eta^{-1}$、质心坐标$\langle z\rangle$和波矢$k$相等, 即

其中变分参数$\theta$, $\eta$, $\langle z\rangle$, $k$和$\varphi_s$都是时间的函数; $\varphi_{s}$代表两个组分中亮孤子的相位; $\theta$可描述两个亮孤子的振幅比. 将波函数(3)代入拉格朗日量(2)式, 并定义$\varphi_{\pm} = \left(\varphi_{\uparrow}\pm \varphi_{\downarrow}\right)/2$及$\beta = k_{\rm R} k+\varepsilon$, 则

然后利用欧拉-拉格朗日方程$\dfrac{\partial \mathcal{L}}{\partial A}-\dfrac{{\rm d}}{{\rm d}t}\left( \dfrac{\partial \mathcal{L}}{ \partial \dot{A}}\right) =$0 ($A$为变分参数$\theta$, $k$, $\langle z\rangle$, $\varphi_{\pm}$及$\eta$, 且$\dot{A} = {{\rm d}A}/{{\rm d}t}$), 可得变分参数随时间演化的运动方程

方程(5a)和(5b)说明孤子的宽度仅由原子间的相互作用决定, 且两个孤子的动量守恒. 方程(5c)—(5e)表明, 自旋-轨道耦合将$\theta$、$\langle z\rangle$及$\varphi_{-}$耦合在一起, 并且线性塞曼劈裂也影响这些参数的动力学演化, 这将导致孤子展现出有趣的动力学特征.

3.   静态孤子

通过设$\dfrac{{\rm d}\theta}{{\rm d}t} = \dfrac{{\rm d} \langle z\rangle}{{\rm d}t} = \dfrac{{\rm d}\varphi_{-}}{{\rm d}t} = 0$, 首先求解欧拉-拉格朗日方程(5a)—(5e)的固定点解(用符号“$\sim$”作为上标标记), 它或许对应于静态孤子解^[1,53]. 从方程(5c)和(5d)中可看出, $\tilde{\varphi}_{-} = n\dfrac{{\text{π}}}{2}$ ($n$为整数)和$\tilde{\langle z\rangle} =$任意常数, 后者说明静态孤子具有平移对称性. 而固定点解$\tilde{\theta}$和$\tilde{k}$满足非线性方程组$\tilde{k} \!=\! k_{{\rm R}}\cos(2\tilde{\theta})$和$0 \!=\! k_{\rm R}^2\cos(2\tilde{\theta})+\varepsilon+(-1)^{n}\varOmega\cot(2\tilde{\theta})$. 由于$n$取奇数和偶数时, $\tilde{\theta}$和$\tilde{k}$的解仅仅相差一个负号, 因此本文考虑$n = 0$的简单情形. 图1给出了固定点解$\tilde{\theta}$和$\tilde{k}$随线性塞曼劈裂强度的变化, 曲线颜色代表不同的解. 从图1(a)和(b)中可看出, 当$\varOmega/k_{\rm R}^2<1$时, 线性塞曼劈裂存在着一个临界值$\varepsilon_{{\rm c}}$, 它随Raman耦合强度增加而单调递减(见图1(e)). 如果$\varepsilon<\varepsilon_{\rm c}$, $\tilde{\theta}$和$\tilde{k}$分别具有四个不同的解, 否则$\tilde{\theta}$和$\tilde{k}$分别具有两个不同的解. 而对于$\varOmega/k_{\rm R}^2>1$, 不论线性塞曼劈裂强度取何值, $\tilde{\theta}$和$\tilde{k}$都只具有两个不同的解(见图1(c)和(d)). 特别的, 如图1(a)— (d)所示, 归咎于线性塞曼劈裂的存在, 这些固定点所对应的孤子总是具有非零的动量, 并且两个分量间的粒子数不相等, 即

图 1  (a)，(b)$\varOmega/k_{\rm R}^2 = 0.5$时固定点解$\tilde{\theta}$和$\tilde{k}$随线性塞曼劈裂$\varepsilon$的变化; (c)，(d)$\varOmega/k_{\rm R}^2 = 1.5$时固定点解$\tilde{\theta}$和$\tilde{k}$随线性塞曼劈裂$\varepsilon$的变化; (e)$\varOmega/k_{\rm R}^2<1$时临界值$\varepsilon_{\rm c}$随$\varOmega$的变化

Fig. 1.  (a) and (b) show the $\tilde{\theta}$ and $\tilde{k}$ change with $\varepsilon$ for $\varOmega/k_{\rm R}^2 = 0.5$; (c) and (d) display $\tilde{\theta}$ and $\tilde{k}$ change with $\varepsilon$ for $\varOmega/k_{\rm R}^2 = 1.5$; (e) shows the critical value $\varepsilon_{\rm c}$ versus $\varOmega$ for $\varOmega/k_{{\rm R}}^2<1$.

下载: 全尺寸图片 幻灯片

$\displaystyle\int_{-\infty}^{+\infty} \left(\vert\psi_{\uparrow}\vert^2-\vert\psi_{\downarrow}\vert^2\right){\rm d}z = -\tilde{k}/k_{\rm R}\neq 0.$

变分法仅仅是一种近似方法, 有必要将欧拉-朗格朗日方程(5a)—(5e)的固定点解与GP方程的数值解作对比^[1]. 一方面, 可利用虚时演化方法求解GP方程的静态孤子数值解, 其结果展示在图2(a)和(b)中. 对于弱的自旋-轨道耦合($\varOmega/k_{\rm R}^2\gg 1$), 图2(a)展示欧拉-朗格朗日方程(5a)—(5e)的固定点解与GP方程的静态孤子数值解相一致. 然而, 对于强自旋-轨道耦合($\varOmega/k_{\rm R}^2\ll 1$)的情况, 两者存在着明显的差别(见图2(b)), 图2(a), (b)中变分静态孤子解分别取自于图1(d)和1(a)中的蓝色和红色曲线. 因此, 在弱自旋-轨道耦合情况下, 变分法能产生一个较好的静态孤子近似解. 另一方面, 也可以将欧拉-朗格朗日方程(5a)—(5e)的固定点解作为初始条件去数值求解含时GP方程, 如果孤子在含时演化中能够保持其初始波形而不运动, 则欧拉-拉格朗日方程(5a)—(5e)的固定点解可认作是GP方程的静态孤子近似解. 这些含时演化结果展示在图2(c)—(f)中, 从中可以看出, 对于弱自旋-轨道耦合的情况($\varOmega/k_{\rm R}^2\gg 1$), 尽管线性塞曼劈裂导致初始孤子具有一个有限的动量, 但是这些孤子总是能保持其初始的波形而静止在初始位置(见图2(c)，(d)). 相反, 对于强自旋-轨道耦合的情形($\varOmega/k_{\rm R}^2\ll 1$), GP方程的含时数值演化结果表明, 孤子将偏离初始位置, 它不仅沿着$z$方向线性运动, 而且运动过程中还会出现振荡运动(见图2(e)，(f)).

图 2  (a)和(b)分别展示$k_{\rm R} = 0.2\varOmega$和$k_{\rm R} = 1.5\varOmega$时, 变分静态孤子解(圆圈)与GP方程（2）静态孤子的数值解(实线)的对比, 其他参数取值为$\varepsilon = 0.3$, $\varOmega = 0.5$及$g = -10$; (c)—(f)分别为(a)和(b)中的变分静态孤子解作为初始条件在含时GP方程中的动力学演化

Fig. 2.  (a), (b) show the comparisons between the variationally predicted stationary soliton solutions (circles) and the numerical solutions (solid lines) of stationary solitons of GP equation (2) for $k_{\rm R}=0.2\varOmega$ and $k_{\rm R}=1.5\varOmega$ with $\varOmega=0.5$, respectively. The other parameters are $\varepsilon=0.3$ and $g=-10$；(c)−(f) are the dynamical evolutions of solitons in time-dependent GP simulations by using the variationally predicted stationary soliton solutions in (a) and (b) as initial wave functions, respectively

下载: 全尺寸图片 幻灯片

利用线性稳定性理论可进一步分析这些静态孤子在微扰下的稳定性及激发模式^[1,53]. 设变分参数$A(t) = \tilde{A}+{\text{δ}} {\rm e}^{{\rm i}\omega t}$, 其中${\text{δ}} A$代表在微扰下变分参数相对于其固定点解$\tilde{A}$的偏离, $\omega$为激发的本征频率. 将$A(t) = \tilde{A}+{\text{δ}} A {\rm e}^{{\rm i}\omega t}$代入欧拉-拉格朗日方程, 并保留至${\text{δ}} A$的一阶项, 可得矩阵方程:

解该矩阵方程, 可得到本征频率$\omega_1 = \omega_2 = 0$和$\omega_{\pm} = \pm\dfrac{2\varOmega}{\sin\left(2\tilde{\theta}\right)}$, 以及对应的本征矢量${ V}_{1} =$$\left(0, 0, 0, 0\right)^{\rm T},$ ${ V}_{2} \!=\! \left(0, 1, 0, 0\right)^{\rm T}$及${ V}_{\pm} \!=\! \left(0, \dfrac{k_{\rm R}}{\varOmega}\sin^3(2\tilde{\theta})\right.,$$\left.\pm {\rm i}\sin(2\tilde{\theta})\right)^{\rm T}$, 其中上标“${\rm T}$”表示转置. 根据线性稳定性分析理论可知^[1], 所有4个本征频率都为实数, 说明静态孤子在微扰下是动力学稳定的. 与本征矢量${ V}_1$相对应的零能模表明, 孤子在扰动下将保持不变. 而与本征矢量${ V}_2$相对应的零能模就是所谓的Goldstone模. 由于本征矢量${ V}_2$中的${\text{δ}} \langle z\rangle\neq 0$, 而${\text{δ}} k$, ${\text{δ}} \theta$及${\text{δ}} \varphi_{-}$均为0, 说明孤子质心坐标在扰动下将偏离其平衡位置而以零频率振动, 即孤子在扰动下将以固定速度做线性运动, 其平移对称性被破缺^[54]. 对于频率为$\omega_{\pm}$的谐振模, 孤子在扰动下将以频率$\omega_{\pm}$而振荡, 且$\omega_{\pm}$也与线性塞曼劈裂有关(见图3， 图3(a), (b)中的$\tilde{\theta}$分别取图1(a)和(c)中的数据).

图 3  $\varOmega/k_{\rm R}^2 = 0.5$ (a)和$\varOmega/k_{\rm R}^2 = 1.5$(b)时, 频率$\omega_{\pm} = \pm 2\varOmega/\sin\left(2\tilde{\theta}\right)$随线性塞曼劈裂强度$\varepsilon$的变化

Fig. 3.  The frequency $\omega_{\pm}=\pm 2\varOmega/\sin\left(2\tilde{\theta}\right)$ changes with $\varepsilon$ for $\varOmega/k_{\rm R}^2=0.5$ and $\varOmega/k_{\rm R}^2=1.5$ in (a) and (b), respectively.

下载: 全尺寸图片 幻灯片

4.   孤子的运动

通过求解含时欧拉-拉格朗日方程的解, 可研究线性塞曼劈裂对孤子的运动的影响. 由于自旋-轨道耦合将变分参数$\langle z\rangle$、$\theta$及$\varphi_{-}$非线性地耦合在一起, 很难直接求解欧拉-拉格朗日方程的精确解. 为此, 首先引入满足条件$\displaystyle\sum_{s = \uparrow, \downarrow}\vert\chi_{s}\vert^2 = 1$的复值旋量$\chi_{s} = \psi_{s}/\sqrt{n}$, 并定义亮孤子自旋

$S_x = \chi_{\uparrow}\chi_{\downarrow}^{\ast}+\chi_{\uparrow}^{\ast}\chi_{\downarrow} = \sin\left(2\theta\right)\cos\left(2\varphi_{-}\right)$,

$S_y = {\rm i}\left(\chi_{\uparrow}\chi_{\downarrow}^{\ast}-\chi_{\uparrow}^{\ast}\chi_{\downarrow}\right) = -\sin\left(2\theta\right)\sin\left(2\varphi_{-}\right),$

及$S_z = \vert\chi_{\uparrow}\vert^2-\vert\chi_{\downarrow}\vert^2 = -\cos\left(2\theta\right)$. 结合欧拉-拉格朗日方程, 可导出孤子自旋满足的运动方程

该方程为常系数线性微分方程组, 其精确解为

其中

$\varpi = \sqrt{\varOmega^2+\beta^2}$,

$\phi = \arctan\left(b/S_{y, 0}\right)$,

$a = \sqrt{b^2+S_{y, 0}^2}$,

$b = \left(\beta S_{x, 0}-\varOmega S_{z, 0}\right)/\omega ,$

$c = \varOmega S_{x, 0}+\beta S_{z, 0}$, 下标“0”标记$S_{x, y, z}$的初始值, 它们由变分参数$\theta$和$\varphi_{-}$的初始值$\theta_0$和$\varphi_{-, 0}$决定. 从解(8a)—(8c)式中可反解得出$\theta$和$\varphi_{-}$的解, 即

$\theta\left(t\right) = \dfrac{1}{2}\arccos\left(S_z\right)$及$\varphi_{-} = -\dfrac{1}{2}\arctan\left(\frac{S_y}{S_x}\right)$.

对于孤子质心$\langle z\rangle$, 由于欧拉-拉格朗日方程(5c)可写作$\dfrac{{\rm d} \langle z\rangle}{{\rm d}t} = k+k_{\rm R} S_z$, 则将$S_z$的精确解代入, 并积分可得

其中已假设孤子质心的初始值为0. 从精确解中可以看出, 在自旋-轨道耦合、Raman耦合及线性塞曼劈裂的作用下, 孤子的质心运动是周期振荡及线性运动的叠加, 其振荡频率$2\varpi$及线性运动速度$v = k_{\rm R}\beta c/\varpi^2+k$均与线性塞曼劈裂强度有关.

对于给定的初始条件, 我们数值求解GP方程(1a), (1b), 并将GP方程的数值解与变分精确解做对比. 为了研究线性塞曼劈裂对孤子运动的影响, 本文考虑初始条件$\theta_{0} = \dfrac{{\text{π}}}{4}$及$\varphi_{-, 0} = 0$, 并假设BEC的凝聚动量$k = 0$. 如果线性塞曼劈裂为零, 变分精确解为$\theta\left(t\right) = \dfrac{{\text{π}}}{4}$, $\varphi_{-} = 0$及$\langle z\rangle\left(t\right) = 0$, 说明孤子并不会运动, 与GP方程的数值模拟结果一致(见图4(a)—(c)). 图4(a)—(f)其他参数取值为$g = -10$, $\eta = -g/2$, $k_{\rm R} = \sqrt{\varOmega/2}$及$\varOmega = 0.5$. 然而, 当线性塞曼劈裂不为零时, 孤子质心的变分精确解为

图 4  (a)—(f)初始值为$\theta_{0} = \dfrac{{\text{π}}}{4}$, $\varphi_{-, 0} = 0$及$k = 0$的孤子动力学演化　(a)—(c) $\varepsilon = 0$, (d)—(f) $\varepsilon = 0.35$; 孤子振荡周期$T(g)$及速度$v(h)$随线性塞曼劈裂的变化

Fig. 4.  (a)−(f) show the dynamical evolutions of initially balanced solitons with $\theta_0=\dfrac{{\text{π}}}{4}$, $k=0$ and $\varphi_{-, 0}=0$ in GP simulations, $\varepsilon=0$ in (a)−(c), and $\varepsilon=0.35$ in (d)−(f); oscillation period $T(g)$ and moving velocity $v(h)$ of solitons change with the linear Zeeman splitting.

下载: 全尺寸图片 幻灯片

$\langle z\rangle \left(t\right) = -\dfrac{k_{\rm R}\varOmega\varepsilon}{2\left(\varOmega^2+\varepsilon^2\right)^{3/2}}\sin\left(2\sqrt{\varOmega^2+\varepsilon^2}t\right)+\dfrac{k_{\rm R}\varOmega\varepsilon}{\varOmega^2+\varepsilon^2}t,$

如图4(d)—(f)所示, 孤子沿$z$方向线性运动的同时将做振荡运动, 振荡周期$T$及线性运动速度$v$与线性塞曼劈裂的强度有关(见图4(g)，(h), 其中$k_{\rm R} = \sqrt{\varOmega/2}$, $\varOmega = 0.5$). 此外, 我们也选择了其他初始条件进行数值模拟, 并发现线性塞曼劈裂将影响孤子的运动速度及振荡频率, 且变分精确解与GP方程的数值模拟结果相吻合.

5.   结　论

本文研究了线性塞曼劈裂对一维自旋-轨道耦合双组份BEC中亮孤子的动力学性质的影响. 通过选择双曲正弦函数作为亮孤子的变分试探波函数, 可运用变分法导出变分参数随时间演化所满足的欧拉-拉格朗日方程. 求解不含时的欧拉-拉格朗日方程, 在弱自旋-轨道耦合情况下, 获得了两个近似的静态孤子解. 基于线性稳定性分析, 进一步发现了一个零能的Goldstone激发模和一个谐振激发模, 前者对应于在外界扰动下静态孤子的平移对称性破缺, 后者表明静态孤子在外界扰动下将做谐振运动, 其谐振频率也与线性塞曼劈裂有关. 最终, 通过求解含时欧拉-拉格朗日方程, 描述孤子质心运动的精确变分解被获得, 并发现线性塞曼劈裂将明显地影响孤子的运动速度和谐振周期. 所有这些变分计算结果都与GP方程的直接数值模拟相吻合.