删除或更新信息，请邮件至freekaoyan#163.com(#换成@)

线性塞曼劈裂对自旋-轨道耦合玻色-爱因斯坦凝聚体中亮孤子动力学的影响

本站小编 Free考研考试/2021-12-29
摘要:利用变分近似及基于Gross-Pitaevskii方程的直接数值模拟方法, 研究了自旋-轨道耦合玻色-爱因斯坦凝聚体中线性塞曼劈裂对亮孤子动力学的影响, 发现线性塞曼劈裂将导致体系具有两个携带有限动量的静态孤子, 以及它们在微扰下存在一个零能的Goldstone激发模和一个频率与线性塞曼劈裂有关的谐振激发模. 同时给出了描述孤子运动的质心坐标表达式, 发现线性塞曼劈裂明显影响孤子的运动速度和振荡周期.
关键词: 孤子/
玻色-爱因斯坦凝聚/
Gross-Pitaevskii方程/
自旋-轨道耦合

English Abstract


全文HTML

--> --> -->
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方程描述：
${\rm i}\frac{\partial \psi_{\uparrow}}{\partial t}= -\frac{1}{2}·\frac{\partial^2 \psi_{\uparrow} }{\partial z^2}- {\rm i} k_{\rm R} \frac{\partial \psi_{\uparrow}}{\partial z}+\varepsilon \psi_{\uparrow}+\varOmega\psi_{\downarrow}+g n\psi_{\uparrow},\tag{1a}$
${\rm i}\frac{\partial \psi_{\downarrow}}{\partial t}=-\frac{1}{2}·\frac{\partial^2 \psi_{\downarrow} }{\partial z^2}+{\rm i} k_{\rm R}\frac{\partial \psi_{\downarrow}}{\partial z}-\varepsilon \psi_{\downarrow}+\varOmega\psi_{\uparrow}+g n\psi_{\downarrow},\tag{1b}$
方程中$ \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]. 在变分近似方法中, 体系所对应的拉格朗日量为
$\begin{split} & \mathcal{L}(t) \\= & \int_{-\infty}^{+\infty}\bigg[\frac{{\rm i}}{2}\sum_{s = \uparrow, \downarrow}\left(\psi_{s}^{\ast}\frac{\partial \psi_{s}}{\partial t}-\psi_{s}\frac{\partial \psi_{s}^{\ast}}{\partial t}\right)+\!\frac{1}{2}\psi_{\uparrow}^{\ast}\frac{\partial^2 \psi_{\uparrow} }{\partial z^2} \\ & +\frac{1}{2}\psi_{\downarrow}^{\ast}\frac{\partial^2 \psi_{\downarrow} }{\partial z^2}+{\rm i} k_{\rm R}\left(\!\psi_{\uparrow}^{\ast}\frac{\partial \psi_{\uparrow}}{\partial z}-\psi_{\downarrow}^{\ast}\frac{\partial \psi_{\downarrow}}{\partial z}\!\right)-\varepsilon\big(\vert\psi_{\uparrow}\vert^2\\ & -\vert\psi_{\downarrow}\vert^2\big)-\varOmega\left(\psi_{\uparrow}^{\ast}\psi_{\downarrow}+\psi_{\downarrow}^{\ast}\psi_{\uparrow}\right)-\frac{g}{2}n^4\bigg]{\rm d}z, \end{split} $
其中$ \ast $代表复共轭.
对于SU(2)自旋对称的均匀双组分BEC, 假设BEC凝聚在单个准动量为$ k $的态上, 则可以利用双曲正弦函数作为亮孤子的试探波函数, 并取两个分量孤子的宽度$ \eta^{-1} $、质心坐标$ \langle z\rangle $和波矢$ k $相等, 即
$\left( \begin{aligned} \psi_{\uparrow} \\ \psi_{\downarrow}\end{aligned} \right) \!=\! \sqrt{\frac{\eta}{2}}\left( \begin{aligned} \sin\theta{\rm e}^{{\rm i}(k z+\varphi_{\uparrow})} \\ \cos\theta {\rm e}^{{\rm i}(k z+\varphi_{\downarrow})} \\ \end{aligned} \right)\mbox{sech}\left[\eta \left(z-\langle z\rangle\right)\right], $
其中变分参数$ \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 $, 则
$\begin{split}\mathcal{L}\left(t\right) \!=\! & -\langle z\rangle\frac{{\rm d} k}{{\rm d}t}\!+\!\cos\left(2\theta\right)\frac{{\rm d}\varphi_{-}}{{\rm d}t} \!-\!\frac{{\rm d}\varphi_+}{{\rm d}t} \!+\!\beta\cos\left(2\theta\right) \\ & -\varOmega\sin\left(2\theta\right)\cos\left(2\varphi_{-}\right) \!-\!\frac{\eta\left(\eta \!+\! g\right)}{6} \!-\! \frac{1}{2}k^2.\end{split}$
然后利用欧拉-拉格朗日方程$ \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} $), 可得变分参数随时间演化的运动方程
$\eta =-g/2,\tag{5a}$
${{\rm d}k}/{{\rm d}t}= 0,\tag{5b}$
$\frac{{\rm d}\langle z\rangle }{{\rm d}t} = k- k_{{\rm R}}\cos \left( 2\theta \right),\tag{5c}$
$\frac{{\rm d}\theta }{{\rm d}t}= -\varOmega \sin \left( 2\varphi _{-}\right),\tag{5d}$
$ \frac{{\rm d}\varphi _{-}}{{\rm d}t}= -\beta-\varOmega \cot \left( 2\theta \right) \cos \left( 2\varphi _{-}\right).\tag{5e}$
方程(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 $的变化
Figure1. (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方程中的动力学演化
Figure2. (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 $的一阶项, 可得矩阵方程:
$\left( \begin{aligned} &\;\; {\rm i}\omega \quad\quad\;\; 0 \quad\quad\quad\quad 0 \quad\quad\quad\quad\quad 0 \\ & -1\quad\quad {\rm i}\omega \quad-2k_{\rm R}\sin \left( 2\tilde{\theta}\right) \quad\;\; 0 \\ &\;\; 0 \quad\quad\;\;\;\; 0 \quad\quad\quad\quad {\rm i}\omega \quad\quad\quad \quad 2\varOmega \\ & k_{\rm R} \quad\quad\;\;\; 0 \quad -\frac{2\varOmega }{\sin ^{2}\left( 2\tilde{\theta}\right) } \quad \quad\quad {\rm i}\omega\end{aligned}\right)\left( \begin{aligned} {\text{δ}} k \;\;\\ {\text{δ}}\langle z\rangle \\ {\text{δ}} \theta \;\;\\ {\text{δ}} \varphi_{-}\end{aligned}\right) = 0.\tag{6}$
解该矩阵方程, 可得到本征频率$ \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 $的变化
Figure3. 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) $. 结合欧拉-拉格朗日方程, 可导出孤子自旋满足的运动方程
$\frac{{\rm d} S_x}{{\rm d}t}=-2\beta S_y, \tag{7a}$
$\frac{{\rm d} S_y}{{\rm d}t}= 2\beta S_x-2\varOmega S_z,\tag{7b}$
$\frac{{\rm d} S_z}{{\rm d}t}= 2\varOmega S_y.\tag{7c}$
该方程为常系数线性微分方程组, 其精确解为
$S_x =\frac{a\beta}{\varpi}\sin\left(2\varpi t-\phi\right)+\frac{\varOmega c}{\varpi^2},\tag{8a}$
$S_y = a \sin\left(2\varpi t+\phi\right), \tag{8b}$
$S_z = -\frac{a\varOmega}{\varpi}\sin\left(2\varpi t-\phi\right)+\frac{\beta c}{\varpi^2},\tag{8c}$
其中
$ \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 $的精确解代入, 并积分可得
$ \begin{split} \left \langle z\right \rangle\left(t\right) =\; & \frac{k_{\rm R}\varOmega}{2 \varpi^{2}}\left[2{S_{y, 0}\sin }^{2}\left( \varpi t\right)-b\sin\left(2\varpi t\right) \right]\\ & + \left(\frac{k_{\rm R}\beta c}{\varpi^2}+k\right) t, \end{split}\tag{9} $
其中已假设孤子质心的初始值为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) $随线性塞曼劈裂的变化
Figure4. (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方程的直接数值模拟相吻合.
相关话题/运动 孤子 实验 物质 线性
  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效，10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材，产品类型包含电子书、题库、全套资料以及视频，无论您是考研复习、考证刷题，还是考前冲刺等，不同类型的产品可满足您学习上的不同需求。 ...
    考试优惠券 本站小编 Free壹佰分学习网 2022-09-19
  • W波段分布作用速调管的设计和实验研究
    摘要:基于运动学理论、感应电流定理和电荷守恒定律,分析了分布作用谐振腔的渡越时间效应,推导了各个谐振腔工作于${\text{π}}$模的电子注与微波之间的能量转换系数、电子负载电导和电子负载电纳,计算结果显示采用分布作用谐振腔有利于提高速调管的工作效率.利用三维电磁仿真软件,设计了一款工作于W波段的 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 量子微波制备方法与实验研究进展
    摘要:量子微波信号既保留了经典微波信号的空间远距离传播能力,又具有非经典的量子特性,为微波频段量子通信、量子导航及量子雷达等基于大尺度动态空间环境无线传输的量子信息技术提供了可资利用的重要信号源.按照腔量子电动力学系统、超导电路量子电动力学系统和腔–光(电)–力学系统三大类型实验平台,归纳、分析了微 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 正弦波沟槽对湍流边界层相干结构影响的TR-PIV实验研究
    摘要:利用高时间分辨率粒子图像测速(time-resolvedparticleimagevelocimetry,TR-PIV)技术,在不同雷诺数下对光滑壁面和二维顺流向、三维正弦波(two/threedimensional,2D/3D)沟槽壁面湍流边界层流场进行了实验测量,从不同沟槽对湍流边界层相干 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 激光加载下金属锡材料微喷颗粒与低密度泡沫混合实验研究
    摘要:金属材料的微喷是冲击加载下金属表面发生的一种动态破碎现象,微喷研究在很多领域都具有重要意义,包括惯性约束聚变(ICF)和烟火制造等.由于激光实验特有的优势,近几年国内外开展了很多利用强激光驱动冲击加载研究材料微喷过程的实验.利用泡沫材料对微喷颗粒进行静态软回收虽然可以获得颗粒的形态分布、颗粒尺 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • Kiselev黑洞的热力学性质和物质吸积特性
    摘要:本文考虑带有黑洞视界和宇宙视界的Kiselev时空.研究以黑洞视界和宇宙视界为边界的系统的热力学性质.统一地给出了两个系统的热力学第一定律;在黑洞视界半径远小于宇宙视界半径的情况下,近似地计算了通过宇宙视界和黑洞视界的热能.然后,探讨Kiselev时空的物质吸积特性.在吸积能量密度正比于背景能 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 光纤1/f 热噪声的实验研究
    摘要:光纤热噪声是限制光纤传感、测量系统性能的最终因素.但是低频区域呈1/f谱特性的光纤热噪声的形成机制迄今仍然存在争论.实验研究了光纤1/f热噪声水平与光纤内杂质离子浓度和光纤施加张力的关系,验证了这类热噪声来源于光纤内部的机械耗散引起的长度自发抖动,符合热机械噪声的理论假设.关键词:光纤光学/热 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 实验优化设计Sr<sub>2</sub>MgSi<sub>2</sub>O<sub>7</sub>:Eu<sup>2
    摘要:为了得到最长有效余辉时间的Sr2MgSi2O7:Eu2+,Dy3+荧光粉,应用二次通用旋转组合设计对实验进行全程优化,建立了稀土离子掺杂浓度Eu2+,Dy3+和有效余辉时间的二元二次回归方程模型,应用遗传算法计算得到有效余辉时间的理论最大值.采用高温固相法合成了最优掺杂浓度Sr2MgSi2O7 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 氦离子显微镜对钨中氦行为的实验研究
    摘要:针对热核聚变面向等离子体钨材料中氦泡形成、演变以及机理研究的需求,克服目前常用离子注入、电子扫描显微镜和透射电子显微镜等离线研究手段存在的不足,提出氦离子显微镜对钨中氦的上述行为原位实时在线研究方法.借助氦离子显微镜的离子注入、显微成像和聚焦离子束纳米加工功能,它可以提供能量为0.5—35ke ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 类KBe<sub>2</sub>BO<sub>3</sub>F<sub>2</sub>结构硼酸盐深紫外非线性光学材料的研究进展
    摘要:利用非线性光学(NLO)晶体材料和变频技术,可以把波长范围有限的激光光源扩展到紫外、深紫外区,这已成为深紫外光源的热点研究方向.然而,目前限制深紫外全固态激光器发展和应用的关键问题是缺乏能够在该波段进行频率转换并且产业化应用的NLO晶体材料.因此,该领域的各国科学家都在积极探索并发展新一代的深 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 微结构气体探测器中紫外激光束的信号和指向精度实验研究
    摘要:在气体探测器研究中,利用266nm紫外激光的双光子电离物理机制使气体电离产生可测量的信号,是一种重要的标定方法.随着微结构气体探测器(MPGD)的不断发展,用紫外激光标定来实现较高精度位置分辨率成为了一种研究需求,对此有两个关键技术问题需要解决:实验研究激光可测信号大小以及激光指向精度.分析和 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29

Free考研考试FreeKaoYan.Com
欢迎来到Free考研考试，"为实现人生的Free而奋斗"

© 2020 FreeKaoYan! . All rights reserved.