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

旋量玻色-爱因斯坦凝聚体拓扑性质的研究进展

本站小编 Free考研考试/2021-12-29

摘要:实现玻色-爱因斯坦凝聚的原子大多具备内部自旋自由度, 在光势阱下原子内部自旋被解冻, 从而使原子可以凝聚到各个超精细量子态上, 形成旋量玻色-爱因斯坦凝聚体. 灵活的自旋自由度成为体系相关的动力学变量, 可以使体系出现新奇的拓扑量子态, 如自旋畴壁、涡旋、磁单极子、斯格明子等. 本文综述了旋量玻色-爱因斯坦凝聚的实验和理论研究, 旋量玻色-爱因斯坦凝聚体中拓扑缺陷的种类, 以及两分量、三分量玻色-爱因斯坦凝聚体中拓扑缺陷的研究进展.
关键词: 旋量玻色-爱因斯坦凝聚体/
自旋畴壁/
涡旋/
斯格明子/
磁单极子

English Abstract


--> --> -->
对于一个多粒子的玻色系统, 当体系的温度极低, 且原子间的距离足够靠近时, 大量原子将会凝聚在动量为零的最低能态上, 从而出现一种具有宏观量子特性的简并态, 即玻色-爱因斯坦凝聚(Bose-Einstein condensate, BEC)[1-5]. 当BEC被限制在磁势阱(如四极阱(Paul)[6], Ioffe-Pritvhard 阱(IPT)[7], 时间轨道势阱(time orbiting potential, TOP)[8]等)中时, 原子处于弱场束缚态, 内部自旋自由度被冻结, 此时凝聚体性质可以通过标量形式的序参量来描述[9,10]. 当BEC被限制在光势阱[11,12]中时, 内部自旋自由度被释放, 其磁化性质由自旋相互作用决定, 这样的体系称为旋量BEC[13-16].
实验方面, 1998年美国麻省理工学院(MIT)的Stamper-Kurn研究组[16]在沿磁势阱的中心轴射入红外线光束形成的光偶极阱中, 首次实现了金属23Na原子系统的BEC, 通过这一突破性进展把在磁势阱中冻结的自旋自由度解放出来, 在冷原子物理领域中开辟了旋量BEC研究的崭新领域. 实验采用光势阱捕获和冷却钠原子, 当温度达到1—2 μK, 原子密度数约为1 × 1014 cm–3时, 转移到磁势阱中经由射频诱导蒸发进一步冷却, 在 F = 1, mF = –1的电子基态上产生含有(5—10) × 106个原子的凝聚体. 随着红外线激光束能量的提升, 在稳定的磁势阱中控制凝聚体进入光势阱, 然后关闭磁势阱, 让凝聚体处于全光阱中. 图1为光势阱中F = 1旋量BEC所有精细态上钠原子的分布, 从图1(b)可以清晰地看到凝聚体分布在三个超精细态上, 金属钠原子的自旋自由度被充分释放.
图 1 光势阱中F = 1 23Na凝聚体的超精细态[16]. (a) 250 ms时光势阱中钠原子的吸收图像; (b) 340 ms时光势阱中钠原子的吸收图像
Figure1. Optical trapping of 23Na condensates in all F = 1 hyperfine states: shown are absorption images after (a) 250 ms and (b) 340 ms of optical confinement.

2001年, Barrett等[17]在两束CO2激光束(功率为12 W, 波长为10.6 μm)交叉会聚形成的全光阱中得到了F = 1的87Rb原子的旋量BEC. 实验中首先将低速的铷原子装载到磁光阱中, 接着逐步降低光强, 关闭磁场, 然后把原子抽运到F = 1, mF = 1态, 最后打开CO2激光束, 使原子装载到光势阱中. 整个过程中是通过在2 s时间内逐步降低光强, 即减少阱深来实现蒸发冷却的. 图2为实验得到的Stern-Gerlach梯度磁场中自由膨胀10 ms后的原子吸收图像, 从下到上分别是F = 1, mF = –1, 0, 1凝聚体的三个分量.
图 2 铷原子云在Stern-Gerlach梯度磁场中自由膨胀10 ms后的吸收图像[17]. 从下到上分别是F = 1, mF = (–1, 0, 1)凝聚体的三个分量
Figure2. Absorptive image of Rb atomic cloud after 10 ms free expansion in a Stern-Gerlach magnetic field gradient. Three distinct components are observed corresponding to F = 1, mF = (–1, 0, 1) spin projections from bottom to top, respectively.

2003年, Gustavson等[18]研究了F = 2的23Na旋量BEC, 观察到了数秒寿命的精细态(F = 2, mF = –2), 且原子密度数为每立方厘米1014 . Chang[19]、Schmaljohann[20]、Kuwamoto[21]等观察到了F = 2的87Rb旋量BEC, 并研究了其动力学性质. 2012年, Pasquiou等[22]也从实验上观察到了F = 3的52Cr旋量BEC, 并研究了其热力学性质, 表明在极低磁场下, 当自旋自由度被热激活时, 凝聚体的临界温度会降低.
特别地, 在2011年, 美国国家标准与技术研究院的Spielman[23]研究组在87Rb旋量BEC中两个内部自旋态之间通过双光子拉曼过程实现了等效的自旋-轨道耦合, 这是人类历史上首次通过人工合成磁场实现了BEC的自旋-轨道耦合效应, 成为近年来冷原子物理中的一个研究热点[24-27]. 2012年, 中国科技大学潘建伟[28]研究组也实现了BEC中的一维自旋-轨道耦合. 此外, 山西大学[29]和麻省理工[30]的研究组, 分别用40K和6Li实现了费米冷原子气体中的自旋-轨道耦合. 量子气体中自旋-轨道耦合的实现, 开辟了研究冷原子物理的新方向, 如光与原子相互作用对旋量BEC拓扑结构的影响, 自旋涡旋阵列的产生, 以及各种新奇拓扑缺陷的产生. 科研工作者不仅详细研究了自旋-轨道耦合的BEC的性质[31-34], 还在实验上成功模拟了Zitterbewegung(狄拉克颤动)[35,36]、自旋霍尔效应[37]等现象, 并提出了研究量子霍尔效应[38]、反常量子自旋霍尔效应[39]、费米气体中的拓扑量子相变[40]等物理问题的方案.
理论方面, 1998年Ho[41]和Ohmi[42]独立进行了旋量BEC的理论研究. 他们推广了Gross-Pitaevskii方程, 在平均场理论下研究了自旋F = 1的BEC的基态结构和拓扑缺陷的性质. 同年, Law等[43]利用量子光学中的代数方法研究了其自旋混合动力学行为. 接着, Koashi[44,45]、Ciobanu[46]和Semenoff[47]等研究了F = 2的旋量BEC的基态结构、自旋相关、磁效应、激发谱及相位等一系列理论性质. Santos[48]、Diener[49]、Makela[50]等对F = 3的旋量BEC也进行了理论研究. 2007年, Yip[51]研究了自旋F ≤ 4的旋量BEC的基态及其对称性.
BEC作为一种新颖的物质形态, 是物理学前沿的研究热点之一, 特别是实验上实现了旋量BEC、分子的BEC、费米凝聚体、自旋-轨道耦合的超冷原子气体以来, 又掀起了新一轮的研究热潮. 本文综述了旋量BEC的实验和理论研究, 旋量BEC中产生的拓扑缺陷的种类, 如自旋畴壁、涡旋、磁单极子、skymions、扭结, 并结合我们的工作着重介绍了两分量和三分量旋量BEC中拓扑缺陷的研究进展.
2
2.1.旋量玻色-爱因斯坦凝聚体的哈密顿量表示
-->自旋为f的BEC的序参量在旋转下有2f + 1个分量, 且这些分量是随着时间空间变化的, 携带关于超流和磁性的信息, 能形成自旋流及质量流, 产生丰富的自旋结构[52].
考虑自旋f的全同玻色子体系, 取$ {{\hat \psi}_m}({{r}}) $为相应的场算符, 假设满足对易关系:
$\left[ {{{\hat \psi }_m}\left( {{r}} \right),{\hat \psi}_n^{\dagger} \left( {{{r}}'} \right)} \right] = {\delta _{nm}}\delta \left( {{{r}} - {{r}}'} \right),$
$\left[ {{{\hat \psi }_n}\left( {{r}} \right),{{\hat \psi }_m}\left( {{{r}}'} \right)} \right] = 0,$
$\left[ {{\hat \psi}_m^{\dagger} \left( {{r}} \right),{\hat \psi}_n^{\dagger} \left( {{r}} \right)} \right] = 0.$
考虑两个自旋为f的玻色子的耦合作用, 其总自旋F为0, 2, …, 2f, 耦合态对应的湮灭算符为
$ \begin{split}{\hat A_{{\rm{FM}}}}\left( {{{r}},{{r'}}} \right) = & \mathop \sum \limits_{{m_1},{m_2} = - f}^f \langle F,M|f,{m_1};f,{m_2} \rangle \\ & {\hat \psi}_{{m_1}}^{}\left( {{r}} \right)\psi _{{m_2}}^{}\left( {{{r}}'} \right).\end{split}$
哈密顿量具有标量性, 在旋转下保持不变, 进一步得到相互作用哈密顿量:
$ \begin{split} {\hat V^F} =\, & \frac{1}{2} \mathop \int {\rm{d}}{{r}}\mathop \int \limits {\rm{d}}{{r}}'{v^{\left( F \right)}}\left( {{{r}} - {{r}}'} \right)\\ & \times\mathop \sum \limits_{M = - F}^F \hat A_{{\rm{FM}}}^{\dagger} \left( {r,r'} \right){\hat A_{{\rm{FM}}}}\left( {{{r}},{{r}}'} \right), \end{split}$
其中${v^{\left( F \right)}}\left( {{{r}} - {{r}}'} \right)$表示相互作用势, 且
$\sum _{{F}} \sum_{{{M}} = - F}^F \hat A_{{\rm{FM}}}^{\dagger} ( {{r}},{{r'}} ){\hat A_{{\rm{FM}}}} \left({{{r}},{{r'}}} \right)\! =\! :\hat n ( {{r}} )\hat n ( {{{r'}}} ):,$
${\rm{\hat n}}\left( r \right) \equiv \mathop \sum \limits_{m = - f}^f {\hat \psi}_m^{\dagger} \left( {{r}} \right){{\hat \psi}_m}\left( {{r}} \right),$
其中${{\hat n}}\left( {{r}} \right)$是粒子数密度算符, ::表示取正规序. 当相互作用势的量级与原子平均间隔相比忽略不计时, 可以将${v^{\left( F \right)}}$近似表示为δ函数的形式:
${v^{\left( F \right)}} = {g_{\rm{F}}}\delta \left( {{r}} \right),$
${g_{\rm{F}}} = \frac{{4{\text{π}}{\hbar ^2}}}{m}{a_{\rm{F}}},$
其中gF描述两个粒子间相互作用强度, aF指S波散射长度. 由此哈密顿量可以表示为
${\hat V^F} = \frac{{{g_{\rm{F}}}}}{2}\int {\rm{d}}{{r}}\mathop \sum \limits_{M = - F}^F \hat A_{{\rm{FM}}}^{\dagger} \left( {{r}} \right){\hat A_{{\rm{FM}}}}\left( {{r}} \right),$
其中${\hat A_{{\rm{FM}}}}\left( {{r}} \right)$如(4)式所示
$ \langle 0,0|f,{m_1};f,{m_2} \rangle = {\delta _{{m_1} + {m_2},0}}\frac{{\left( { - 1} \right){f^{ - {m_1}}}}}{{\sqrt {2f + 1} }},$
${\hat A_{00}}\left( {{{r}},{{r'}}} \right) \!=\! \frac{1}{{\sqrt {2f \!+\! 1} }}\!\!\mathop \sum \limits_{ - f}^f \!{\left( { - 1} \right)^{f - m}}{{\hat \psi}_m}\left( {{r}} \right){{\hat \psi}_{ - m}}\left( {{{r}}'} \right).$
对于自旋f = 1时的相互作用哈密顿量, 此时两个玻色子碰撞的总自旋F只能是0和2, 相应的相互作用哈密顿量可表示为
${\hat V_0} = \frac{{{g_0}}}{2}\int{{\rm{d}}^3}{{r}}\hat A_{00}^{\dagger} \left( {{r}} \right){\hat A_{00}}\left( {{r}} \right),$
$\begin{split}{\hat V_2} \, & = \frac{{{g_2}}}{2}\int {{\rm{d}}^3}{{r}}\mathop \sum \limits_{M = - 2}^2 \hat A_{2M}^{\dagger} \left( {{r}} \right){\hat A_{2M}}\left( {{r}} \right) \\ & = \frac{{{g_2}}}{2}\int {{\rm{d}}^3}\left[ {:{{\hat n}^2}\left( {{r}} \right): - \hat A_{00}^{\dagger} \left( {{r}} \right){{\hat A}_{00}}\left( {{r}} \right)} \right],\end{split}$
联合两式得到总的相互作用哈密顿量为
$\hat V = \mathop \int \limits {{\rm{d}}^3}\left[ {\frac{{{g_2}}}{2}:{{\hat n}^2}\left( {{r}} \right): + \frac{{{g_0} - {g_2}}}{2}\hat A_{00}^{\dagger} \left( {{r}} \right){{\hat A}_{00}}\left( {{r}} \right)} \right],$
将自旋1情况的算符恒等式$:{{\hat F}}\left( {{r}} \right) \cdot {{\hat F}}\left( {{r}} \right): +$$ 3\hat A_{00}^{\dagger} \left( {{{r}}, {{r'}}} \right){\hat A_{00}}\left( {{{r}}, {{r'}}} \right) \!=: \hat n\left( {{r}} \right)\hat n\left( {{{r'}}} \right):$代入上式, 得到:
$\hat V = \frac{1}{2}\int {{\rm{d}}^3}\left[ {{c_0}:{{\hat n}^2}\left( {{r}} \right): + {c_1}:{{{{\hat F}}}^2}\left( {{r}} \right):} \right],$
其中
${c_0} = \frac{{{g_0} + 2{g_2}}}{3},$
${c_1} = \frac{{{g_2} - {g_0}}}{3},$
同样地, 可以得到f = 2时的相互作用哈密顿量:
$\begin{split} {\hat V} =\, & \frac{1}{2}\int {{\rm{d}}^3}\left[{c_0}:{{\hat n}^2}\left( {{r}} \right): + {c_1}:{{{{\hat F}}}^2}\left( {{r}} \right):\right. \\ & +\left. {c_2}\hat A_{00}^{\dagger} \left( {{r}} \right){{\hat A}_{00}}\left( {{r}} \right)\right], \end{split}$
其中
${c_0} = ({{4{g_2} + 3{g_4}}})/{7},$
${c_1} = ({{{g_4} - {g_2}}})/{7},$
${c_2} = ({{7{g_0} - 10{g_2} + 3{g_4}}})/{7},$
f = 3时的相互作用哈密顿量为
$\begin{split}\hat V = \, &\frac{1}{2} \int{{\rm{d}}^3}\bigg[{c_0}:{{\hat n}^2}\left( {{r}} \right): + {c_1}:{{{{\hat F}}}^2}\left( {{r}} \right): \\& + {c_2}\hat A_{00}^{\dagger} \left( {{r}} \right){{\hat A}_{00}}\left( {{r}} \right) \\& + {c_3}\mathop \sum \limits_{M = - 2}^2 \hat A_{2M}^{\dagger} \left( {{r}} \right){{\hat A}_{2M}}\left( {{r}} \right) \bigg],\end{split}$
其中
${c_0} = ({{9{g_4} + 2{g_6}}})/{{11}},$
${c_1} = ({{{g_6} - {g_4}}})/{{11}},$
${c_2} = ({{11{g_0} - 21{g_4} + 10{g_6}}})/{{11}},$
${c_3} = (11g_0 - 18g_4 + 7g_6)/11.$

2
2.2.旋量玻色-爱因斯坦凝聚体的拓扑荷
-->对于各种形式的拓扑缺陷主要通过计算拓扑荷来判定. 线缺陷形式的涡旋可通过观察相位图中的奇异点和缠绕数来识别, 对于点缺陷, 如磁单极子, 可以通过计算拓扑荷来确定是整数磁单极子还是分数磁单极子[53]. 在两分量BEC系统中, 原子可以是两种不同的元素[54], 也可以是同种元素的两种同位素[55], 或者是同种原子的不同超精细态[56], 该体系可用贋自旋${1}/{2}$系统来描述[57,58].
贋自旋${{S}} = {\varPsi ^{\dagger} }\sigma \varPsi /{\left| \varPsi \right|^2}$, 其分量形式为
$\left\{ {\begin{aligned}& {{S_{{x}}} = \frac{{2\left| {{\psi _1}} \right|\left| {{\psi _2}} \right|}}{{{{\left| \varPsi \right|}^2}}}\cos \left( {{\theta _1} - {\theta _2}} \right)}\\& {{S_{{y}}} = - \frac{{2\left| {{\psi _1}} \right|\left| {{\psi _2}} \right|}}{{{{\left| \varPsi \right|}^2}}}\sin \left( {{\theta _1} - {{\rm{\theta }}_2}} \right)}\\& {{S_{{z}}} = \frac{{{{\left| {{\psi _1}} \right|}^2} - {{\left| {{\psi _2}} \right|}^2}}}{{{{\left| \varPsi \right|}^2}}}}\end{aligned}} \right.,$
其中ψ1, ψ2, θ1, θ2分别为两个分量的波函数和相位. 局部贋自旋的三个分量可表示为
$\left\{ {\begin{aligned}& {{S_{{x}}} = \sin \theta \cos \varphi }\\& {{S_{{y}}} = \sin \theta \sin \varphi }\\& {{S_{{z}}} = \cos \theta }\end{aligned}} \right.,$
其中θφ分别为极角和方位角. 比较两式可知贋自旋的极角只与两分量的相对密度有关, 即贋自旋的Sz分量由两组分的密度差决定, 而Sx, Sy分量由相对密度和相对相位共同决定, 表明两分量BEC不同的相对密度和相对相位会导致不同结构的自旋纹理, 可以用拓扑荷密度q(r)来描述其拓扑结构的空间分布:
$q\left( r \right) = \frac{1}{{8{\text{π}} }}{\varepsilon ^{ij}}{\bf{S}} \cdot {\partial _i}{{S}} \times {\partial _j}{{S}},$
对全空间积分得拓扑荷:
$\begin{split}{\cal Q}\left( {{S_{{x}}},{S_{{y}}},{S_{{z}}}} \right) & = \int \dfrac{1}{{8{\text{π}}}}{\varepsilon ^{ij}}{{S}} \cdot {\partial _i}{{S}} \times {\partial _j}{{S}}{\rm{d}}{{r}}\\ & = \dfrac{1}{{4{\text{π}}}}\int {\left| {\begin{array}{*{20}{c}}{{S_{{x}}}}&{{S_{{y}}}}&{{S_z}}\\{\dfrac{{\partial {{{S}}_{{x}}}}}{{\partial x}}}&{\dfrac{{\partial {S_{{y}}}}}{{\partial x}}}&{\dfrac{{\partial {S_{\rm{z}}}}}{{\partial x}}}\\{\dfrac{{\partial {S_{{x}}}}}{{\partial y}}}&{\dfrac{{\partial {S_{{y}}}}}{{\partial y}}}&{\dfrac{{\partial {S_{{z}}}}}{{\partial y}}}\end{array}} \right|{\rm{d}}{{r}}}\\ &= - \dfrac{1}{{4{\text{π}}}}\int {\left| {\begin{array}{*{20}{c}}{{S_{{x}}}}&{{S_{{z}}}}&{{S_y}}\\{\dfrac{{\partial {{{S}}_{{x}}}}}{{\partial x}}}&{\dfrac{{\partial {S_{{z}}}}}{{\partial x}}}&{\dfrac{{\partial {S_{{y}}}}}{{\partial x}}}\\{\dfrac{{\partial {S_{{x}}}}}{{\partial y}}}&{\dfrac{{\partial {S_{{z}}}}}{{\partial y}}}&{\dfrac{{\partial {S_{{y}}}}}{{\partial y}}}\end{array}} \right|{\rm{d}}{{r}}} \\ & = - {\cal Q}\left( {{S_{{x}}},{S_z},{S_{{y}}}} \right).\\[-14pt]\end{split}$
从(31)式可以看出, 任意交换自旋密度矢量的三个分量Sx, Sy, Sz中的两个分量, 或者任意改变三个分量中其中一个分量的正负号, 将会出现不同结构的自旋纹理, 但其拓扑荷密度q(r)和拓扑荷的绝对值$\left| {\cal Q} \right|$不会变化.
拓扑缺陷的独特性是在弱微扰下能稳定存在, 在空间自由连续地变换而不改变其自身性质. 这种稳定性由表征序参量流形的离散拓扑量子数保护. 拓扑缺陷的分类主要依靠同伦论, 描述不同的拓扑缺陷会出现不同的序参量流形. 表1总结了利用同伦群对拓扑缺陷进行分类的结果. 旋量BEC有丰富的序参量流形, 因而可产生不同种类的拓扑缺陷, 如整数和分数涡旋[59-64], 非阿贝尔涡旋[65],’ tHooft-Polyakov磁单极子[66], 狄拉克磁单极子[67-68], skyrmions[69-71]和扭结[72].
πn缺陷孤子
π0磁畴壁暗孤子
π1涡旋非奇异磁畴壁
π2磁单极二维skyrmions
π3skyrmions, 扭结
π4瞬子


表1同伦群描述的拓扑缺陷结构
Table1.Topological defect structures described by homotopy groups.

2
3.1.自旋畴壁
-->畴壁可由零阶同伦群${{\text{π}}_0}\left( {\cal R} \right)$表征. 若${{\text{π}}_0}\left( {\cal R} \right) = 0$说明体系是相连的, 若${{\text{π}} _0}\left( {\cal R} \right) = 1$, 说明体系被分为两个非连续区域. 以两分量赝自旋${1}/{2}$的序参量来描述, 其赝自旋表示为(28)式. 图3表示在赝自旋表象中自旋密度Sx, Sy, Sz的分布情况. 第一组分分布在势阱的左边, 标记为自旋朝上Sz = 1(图中红色表示); 第二组分在强烈的排斥作用下分布在第一组分的低密度区域, 即势阱的右边, 标记为自旋朝下Sz = –1(图中蓝色表示). 两个分量的相分离导致在赝自旋表象中形成了两个自旋畴, 在两个自旋畴的界面上形成了自旋畴壁. 在这个区域自旋既不朝上也不朝下(即$\left| {{S_{{z}}}} \right| \ne 1$), 而是在x方向上有了分量. 从左边自旋畴到右边自旋畴的过程中, 畴壁上的自旋沿着x方向翻转, 赝自旋在x-y平面的投影均指向x轴的正方向, 形成了如图3(d)中的经典奈尔型畴壁[73].
图 3 赝自旋密度Sz, Sx, Sy的空间分布[73] (a)?(c)表示旋转角频率为0; (d)自旋纹理投影到x-y平面内的矢量表示
Figure3. The pseudospin density distribution for (a) Sz, (b) Sx and (c) Sy for Ω = 0; (d) the vectorial representation of the spin texture projected onto the x-y plane.

2
3.2.涡 旋
-->涡旋可由一阶同伦群${{\text{π}} _1}\left( {\cal R} \right)$表征, 描述一个从实空间回路到序参量流形的映射. 它属于线缺陷, 常见于超流和液晶位错中. 对于S波超导体、液4He以及旋量BEC, 绕涡旋线的质量环流是量子化的, 其涡旋由缠绕数nw表征:
$\mathop \oint \nolimits_C {{{v}}_{\rm{s}}} \cdot {\rm{d}}{{l}} = {n_{\rm{w}}}\kappa,$
积分沿着闭合路径C进行, 其中$\kappa \equiv {h}/{M}$, ${{{v}}_{\rm{s}}}( = {\hbar }/{M}\nabla \phi )$为超流速度. 对于标量序参量, 其序参量流形为${\cal R} = {\rm{U}}\left( 1 \right)$, 其基本群为整数加群: ${{\text{π}} _1}\left( {U\left( 1 \right)} \right) \cong \mathbb{Z}$. ${{\text{π}} _1}\left( {U\left( 1 \right)} \right)$的群元与缠绕数nw的值是一一对应的[74].
对于自旋1的铁磁BEC, 序参量流形为SO(3), 其基本群为${{\text{π}} _1}\left( {SO\left( 3 \right)} \right) \cong {\mathbb{Z}_2} = \left\{ {0, 1} \right\}$, 可存在Mermin-Ho涡旋, 对应序参量为:
$\left( {\begin{array}{*{20}{c}}{{\psi _1}}\\{{\psi _0}}\\{{\psi _{ - 1}}}\end{array}} \right) = \sqrt n \left( {\begin{array}{*{20}{c}}{{\rm{co}}{{\rm{s}}^2}\dfrac{\beta }{2}}\\{\sqrt 2 {{\rm{e}}^{{\rm{i}}\phi }}\sin \dfrac{\beta }{2}\cos \dfrac{\beta }{2}}\\{{{\rm{e}}^{2{\rm{i}}\phi }}{\rm{si}}{{\rm{n}}^2}\dfrac{\beta }{2}}\end{array}} \right),$
其中β是方位角, $\phi $表示极坐标下极角. Mermin-Ho涡旋在中心处的方位角β = 0, 在边界处对应的方位角$\beta = {{\text{π}}}/{2}$. 因此, 空间内自旋方向在中心处是垂直的, 边界处是水平的. 这种涡旋的缠绕数组合为(0, 1, 2), 每个组分的密度分布都是轴对称的, 势阱中心被ψ1组分占据, ψ0组分被推到外部, ψ–1组分被排挤到最边缘区域.
另一种涡旋结构为极核涡旋, 对应序参量为
$\left( {\begin{array}{*{20}{c}}{\begin{array}{*{20}{c}}{{\psi _1}}\\{{\psi _0}}\end{array}}\\{{\psi _{ - 1}}}\end{array}} \right) = \sqrt n \left( {\begin{array}{*{20}{c}}{{{\rm{e}}^{ - {\rm{i}}\phi }}{\rm{co}}{{\rm{s}}^2}\dfrac{\beta }{2}}\\{\dfrac{1}{{\sqrt 2 }}\sin \beta }\\{{{\rm{e}}^{{\rm{i}}\phi }}{\rm{si}}{{\rm{n}}^2}\dfrac{\beta }{2}}\end{array}} \right).$
这种涡旋的缠绕数满足组合(1, 0, –1). 势阱中心被ψ0组分占据, ψ1ψ–1组分被推到外部区域. 不同于Mermin-Ho涡旋中心处的纵向磁化, 极核涡旋代表手征对称性的自发破缺.
对于涡旋而言, 相位的改变量是量子化的, 应是2π的整数倍[75]. 在两分量BEC系统中, 通过光学方法可以实现半量子化涡旋(half-quantized vortice), 描述两分量系统中, 组分一涡旋的相位改变量是2π, 组分二相位改变量为0的情况. 相对于整个系统, 每个组分的相位改变量为π, 因此称之为半量子化涡旋或Alice涡旋[76]. 2017年Liu C F和Liu W M[77]利用变分方法, 得到了自旋-轨道耦合和外磁场联合作用下自旋1BEC中无核半量子涡旋可能的稳定解. 计算表明, 各向同性自旋-轨道耦合提供了稳定无核半量子涡旋和相应的奇异半量子自旋纹理激发的局部能量最小. 各向同性自旋-轨道耦合和垂直磁场的联合约束是获得这种激发的关键因素, 研究还提供了半量子涡旋的稳定相图和最可能的尺寸. 图4为自旋1BEC中半量子涡旋近似解和相应的奇异自旋纹理. 在自旋2BEC的单轴向列相和四面体循环相中还存在非阿贝尔涡旋, 其迷向群都是阿贝尔群[78]. 非阿贝尔涡旋的特点在于其涡旋的生成元是不能互换的, 因此当两个涡旋碰撞时, 它们不能合并或者通过对方, 而是形成一个rung涡旋粘在一起.
图 4 自旋1 BEC中半量子涡旋的近似解和相应的奇异自旋纹理[77] (a)和(b)对应$\left| {F = 1, {{{m}}_{\rm F}} = 0} \right\rangle $$\left| {F = 1, {{{m}}_{\rm F}} = - 1} \right\rangle $分量的密度; (c)和(d)是对应的相; (e)为半量子涡旋的分布; (f)|S|自旋密度; (g) |S|自旋密度分布; (h)自旋纹理; (i)拓扑荷密度$q\left( {x, y} \right)$
Figure4. Approximate half-quantum vortex solution in the spin-1 BEC and the corresponding singular spin texture: (a) and (b) are the densities of the $\left| {F = 1, {m_{\rm F}} = 0} \right\rangle $$\left| {F = 1, {m_F} = - 1} \right\rangle $ components, respectively; (c) and (d) are the corresponding phases; (e) shows the profile of the half-quantum vortex; (f) spin density|S|; (g) the profile of the spin density|S|; (h) spin texture; (i) topological charge density $q\left( {x, y} \right)$.

2
3.3.磁单极子
-->磁单极子可由二阶同伦群${{\text{π}} _2}\left( {\cal R} \right)$表征, 描述从实空间一个球到序参量流形的映射, 为二维skyrmion结构, 常见的两种二维skyrmions构型如图5. 在铁磁BEC系统中可以产生狄拉克磁单极子, 其铁磁相的序参量为
图 5 两种常见的二维skyrmions的矢量场构型[79] (a) 豪猪型skyrmion; (b) 螺旋型skyrmion
Figure5. Two common vector field configurations of two-dimensional skyrmions: (a) The hedgehog type skyrmion; (b) the spiral type skyrmion.

$\begin{split} \left( {\begin{array}{*{20}{c}}{{\psi _1}}\\{{\psi _0}}\\{{\psi _{ - 1}}}\end{array}} \right) & = \sqrt n {{\rm{e}}^{{\rm{i}}\theta }}U\left( {\alpha,\beta,\gamma } \right)\left( {\begin{array}{*{20}{c}}1\\0\\0\end{array}} \right) \\ & = \sqrt n {{\rm{e}}^{{\rm{i}}\left( {\theta - \gamma } \right)}}\left( {\begin{array}{*{20}{c}}{{{\rm{e}}^{ - {\rm{i}}\alpha }}{\rm{co}}{{\rm{s}}^2}\dfrac{\beta }{2}}\\{\dfrac{1}{{\sqrt 2 }}\sin \beta }\\{{{\rm{e}}^{{\rm{i}}\alpha }}{\rm{si}}{{\rm{n}}^2}\dfrac{\beta }{2}}\end{array}} \right).\end{split}$
这样的点缺陷在铁磁BEC中三个自旋态的分布分别为: 在ψ1分量表现为双量子数涡旋线, 在ψ0分量表现为单量子数涡旋线, 在ψ–1分量表现为孤子态, 没有涡旋线.
自1931年狄拉克磁单极子被首次提出以来, 便受到了广泛关注, 尽管到目前为止仍然没有探测到真实的磁单极子, 科研工作者在不同领域也取得了较大的进展, 例如固态物理中的自旋冰材料, 便提供了一种人造磁单极子环境. 最具有突破性进展的是研究者在铁磁BEC中实验上实现了人造磁单极子[80,81].
2
3.4.三维skyrimion
-->三维skyrimion可由三阶同伦群${{\text{π}} _3}\left( {\cal R} \right)$表征, 这是一个扩展到整个三维实空间的拓扑结构, 类似于粒子的拓扑孤子, 所有自旋有序排列, 在有限空间下自旋会发生反转. 考虑一个铁磁系统, 对应的序参量具有$SO(3)$对称性. 它是在位置r处序参量通过一个向量$ \varOmega $表示, 向量的方向和梯度描述了自旋的方向和旋转的角度. 给定条件f(∞) = 0, 保证序参量均匀且在无限远处等于ζFM = (1, 0, 0)T. 给定另一条件f(0) = 2πn保证skyrmion属于同伦分类. 通过酋变换得到序参量表示:
$\left( {\begin{array}{*{20}{c}}{{\psi _1}}\\{{\psi _0}}\\{{\psi _{ - 1}}}\end{array}} \right) = \left( \begin{array}{c}{\left( {\cos \dfrac{{f\left( r \right)n}}{2} - {\rm{i}}\cos \theta \sin \dfrac{{f\left( r \right)n}}{2}} \right)^2} \\ - \sqrt 2 {\rm{i}}\left( {\cos \dfrac{{f\left( r \right)n}}{2} - {\rm{i}}\cos \theta \sin \dfrac{{f\left( r \right)n}}{2}} \right)\sin \dfrac{{f\left( r \right)n}}{2}\sin \theta {{\rm{e}}^{{\rm{i}}\phi }}\\ - {\rm{si}}{{\rm{n}}^2}\dfrac{{f\left( r \right)n}}{2}{\rm{si}}{{\rm{n}}^2}\theta {{\rm{e}}^{2{\rm{i}}\phi }}\end{array} \right),$
其中$\hat f = \left( {{{\hat f}_{{x}}}, {{\hat f}_{{y}}}, {{\hat f}_{{z}}}} \right)$, ($r, \theta, \phi $)是位置r处的极坐标.
Skyrmion在凝聚态的许多体系中发挥了作用, 例如液氦、量子霍尔体系、液晶以及螺旋铁磁系统, 目前在旋量BEC中的skyrimon已经成为研究热点. 图6是研究者在两分量BEC中通过非阿贝尔规范场诱导观察到的三维skyrmions[82,83].
图 6 稳定的三维skyrmions在x-yz-x平面的空间分布[83] (a)中的箭头和颜色分别表示贋自旋方向和OP的U(1)相分布. 彩图(b)和(c)分别表示$\left| {{\varPsi _ \uparrow }\left( {\rm{r}} \right)} \right|$$\left| {{\Psi _ \downarrow }\left( {\rm{r}} \right)} \right|$的振幅
Figure6. The spatial profile of the stable 3D skyrmions in the x-y and z-x planes: The arrows and their colors in (a) indicate the pseudospin direction and the U(1) phase of the OP, respectively; the color maps of (b) and (c) give the amplitudes $\left| {{\varPsi _ \uparrow }\left( {\rm{r}} \right)} \right|$ and $\left| {{\varPsi _ \downarrow }\left( {\rm{r}} \right)} \right|$, respectively.

2
3.5.扭 结
-->扭结也可由三阶同伦群${{\text{π}} _3}\left( {\cal R} \right)$表征, 表现为一个闭合回路互相嵌套的结构, 描述从三维球空间S 3S2的映射. 不同于其他拓扑缺陷, 如涡旋、磁单极子和三维skyrmion通过缠绕数来描述, 扭结是通过连接数或Hopf不变量来描述. Hopf映射是从S3S2, 在S2中一个点的原像会形成一个无结的回路. 假设在BEC系统中给定一个闭合路径C1, 自旋矢量固定在一个给定的方向, 另一个回路C2沿着自旋不同的方向. 如果两个回路C1C2连接一次, 那么连接数就为1. 连接数可正可负, 决定于两个回路的相对取向. 扭结自Faddeev和Niemi提出在三维经典场论中可以稳定孤子的形式存在以来, 引起了物理学家极大的研究兴趣. 与三维skyrmion的情况一样, S3域是通过设置一个边界条件来给定的, 即在空间无穷大的各个方向上, 序参量的值是相同的. 考虑自旋1的极性相, 序参量流形为${\cal{R}^{{\rm{ploar}}}} \cong \left( {S_r^2 \times U{{\left( 1 \right)}_\phi }} \right)/{\left( {{\mathbb{Z}_2}} \right)_{r, \phi }}$, 这里U(1)和${\mathbb{Z}_2}$对一维空间以上的同伦群无对称贡献, 因此, 得到${{\text{π}} _3}{{\cal R}^{{\rm{ploar}}}} \cong {{\text{π}} _3}\left( {{S^2}} \right) \cong \mathbb{Z}$. 对应的拓扑荷, 也即Hopf荷${\cal Q} \in \mathbb{Z}$[84].
通过操纵外部磁场, 可以在自旋1BEC中形成扭结. 在外部磁场存在时, 线性塞曼效应引起${{\hat d}}$的拉莫尔进动, 而二次塞曼效应${{\hat d}}$往往平行与磁场. 假设一个处于光势阱的BEC, 其${{\hat d}}$ = (0 0 1)T, 在z方向施加均匀磁场, 然后突然关闭均匀磁场, 打开四极场, 由于线性塞曼效应, ${{\hat d}}$开始围绕磁场旋转. 此时${{\hat d}}$作为时间的函数运动, 从而导致扭结的形成. 图7显示了四极场作用下球形光势阱中扭结产生的动力学过程.
图 7 四极场作用下球形光势阱中扭结产生的动力学过程[85]. 上一行表示${{\hat d}} = {\left( {0, 0, - 1} \right)^{\rm{T}}}$${{\hat d}}$ = (1, 0, 0)T的图像快照, 下一行表示x-y平面上m = –1分量的密度截面
Figure7. Dynamics of the creation of knots in a spherical optical trap under a quadrupole magnetic field. Snapshots of the preimages of ${{\hat d}}$ = (0, 0, –1)T and ${{\hat d}}$ = (1, 0, 0)T(top), and the cross sections of the density for the m = –1 components on the xy plane (bottom).

2
4.1.两分量旋量玻色-爱因斯坦凝聚体中的拓扑缺陷研究
-->对于两分量BEC[86], 由于组分内相互作用与组分间相互作用两者的竞争导致体系出现丰富而有趣的拓扑结构, 如上文介绍的涡旋、skyrmion、磁单极子和量子扭结等, 现已成为实验研究的理想平台. 实验方面, 1999年Matthews等[87]在两分量BEC中产生了量子涡旋. 2001年Anderson等[88]采用两束正交的探测光观察到双组分的BEC中暗孤子受动力学不稳定的影响衰变为涡旋环. 2016年Hall等[89]演示了在旋量BEC中扭结孤子的实验创建和检测, 所观察到的扭结纹理对应于三阶同伦群, 如图8.
图 8 扭结孤子的结构及其产生方法[89] (a)和(b)为扭结形成之前和形成过程中磁感应线的示意图, 绿色椭圆为对应的凝聚体; (c)和(d)显示扭结形成时, 最初的z方向的向列相矢量(黑色箭头)沿着局部磁场(青色线)的方向进动, 以实现最终的结构(彩色箭头). 灰色虚线表示dz = 0, 白线表示孤子核(dz = –1), 深灰色线表示体积V (dz = 1)的边界; (e)表示实空间中扭结孤子的构型及其与S2中向列矢量${{\hat d}}$的关系
Figure8. Structure of the knot soliton and the method of its creation: Schematic magnetic field lines before (a) and during (b) the knot formation, with respect to the condensate (green ellipse); (c), (d) as the knot is tied, the initially z-pointing nematic vector (black arrows) precesses about the direction of the local magnetic field (cyan lines) to achieve the final configuration (coloured arrows); the dashed grey line shows where dz = 0, the white line indicates the soliton core (dz = –1), and the dark grey line defines the boundary of the volume V (dz= 1); (e) the knot soliton configuration in real space and its relation to the nematic vector ${{\hat d}}$ in S2 (inset).

理论方面, 1999年Williams和Holland[90]展示了在BEC中选择性地产生具有不同角动量超流涡旋的方法, 包括求解具有强耦合原子态的双组分凝聚体的时间相关运动方程. 2002年Battye等[91]证实了在两分量BEC中存在稳定的skyrmion结构. 同年Martikainen等[92]理论研究了在两分量BEC中产生磁单极子的方法, 并证明磁单极子的产生并不局限于反铁磁自旋凝聚, 同时研究了这种磁单极子的膨胀探测, 以及势阱中位移磁单极子的动力学. 2004年Kasamatsu和Tsubota[93]通过数值积分耦合的Gross-Pitaevskii方程, 研究了双组分BEC在轴对称势阱中由调制不稳定性引起的多畴壁形成动力学. 2010年Wang等[94]在无外势的自旋-轨道耦合两分量BEC中发现了平面波相和条纹相, 在考虑外势的情况下, 体系将出现新的量子态, 如分数涡旋和涡旋格子[95]以及skyrmion格子[96].
接着, 人们研究了旋转势下自旋-轨道耦合两分量BEC的基态性质. 2011年Xu和Han[97]在旋转势下自旋-轨道耦合两分量BEC体系中发现了对称排列的涡旋列和中心伴有巨skyrmion的三角涡旋格子. 2011年Zhou等[98]探讨了具有旋转和自旋-轨道耦合两分量BEC的基态性质, 发现了半量子数涡旋格子结构. 2012年Liu等[99]总结了skyrmions的类型, 并通过研究自旋-轨道耦合的两分量BEC的随机投影Gross-Pitaevskii方程, 发现自旋-轨道耦合能够诱导系统产生具有两个Sz极值的环形-双曲状的skyrmions. 图9为总结的skyrmions类型.
图 9 Skyrmions的类型(λ = 0.5)[99] (a)?(h)表示自旋矢量的模式: (a)径向-向外skyrmion, (b)径向-向内skyrmion, (c)环形skyrmion, (d)双曲skyrmion, (e)双曲-径向向外skyrmion, (f)双曲-径向向内skyrmion, (g)环形-双曲skyrmion-I, (h)环形-双曲skyrmion-II
Figure9. Configuration of the skyrmion where λ = 0.5: The (a)?(h) figures indicate the mode of the spin vectors: (a) radial-out skyrmion, (b) radial-in skyrmion, (c) circular skyrmion, (d) hyperbolic skyrmion, (e) hyperbolic-radial(out) skyrmion, (f) hyperbolic-radial (in) skyrmion, (g) circular-hyperbolic skyrmion-I, and (h) circular-hyperbolic skyrmion-II[99].

2014年Wang等[100]研究了准二维旋转双分量BEC在改变自旋-轨道耦合和旋转频率强度的情况下呈现出的各种丰富的基态结构, 表明各向异性自旋-轨道耦合引起的不同基态相之间的转变明显不同于各向同性的情况. Fetter[101]用时间相关的拉格朗日形式和变分函数研究了自旋-轨道耦合的BEC中双组分涡旋的动力学. 2016年Sakaguchi和Umeda[102]通过数值模拟和变分方法, 研究了Rashba型自旋-轨道耦合的双组分BEC的Gross-Pitaevskii方程. 发现当不存在相互作用时, 多量子涡旋态成为谐波势中的基态. 当引力相互作用较强时, 多量子涡旋态在方位角方向呈现调制不稳定性, 出现孤子态. 当排斥相互作用较强时, 形成中心为多量子涡旋的涡旋晶格态, 且涡旋晶格态近似于多量子涡旋态的线性组合.
2017年Sakaguchi[103]通过数值分析研究了双组分BEC中自旋-轨道耦合条纹和半涡旋物质波孤子在组分间线性混合(Rabi耦合)作用下, 以耦合结构振荡和穿梭运动形式产生宏观量子效应的可能性. 研究得到在一维系统中, 本征振荡表现为条纹孤子在空间上的偶数分量和奇数分量之间的翻转, 而在二维系统中则表现为半涡旋孤子的零涡旋和涡旋分量之间的周期性跃迁. 同年, Wang等[104]研究了环形势阱中自旋-轨道耦合两分量BEC的旋涡态和自旋纹理, 系统地讨论了旋转、自旋-轨道耦合和原子间相互作用对系统基态涡旋结构和自旋纹理的影响. 特别是当旋转频率固定在临界值以上时, 各向同性自旋-轨道耦合的增强使每个分量中有一个可见的涡旋链, 在中心区域伴随着一个隐藏的巨涡旋和一个(或几个)隐藏的涡旋链. Kato等[105]研究了具有Rashba自旋-轨道耦合双组分BEC中涡旋-反涡旋对的动力学性质, 并且发现旋涡-反涡旋对的速度远小于无自旋-轨道耦合的速度且存在稳态, 而两个具有相同环流的涡旋相互移动或联合会形成一个静止状态.
2018年Shi等[106]研究了旋转非对称势阱中具有Dreselhaus自旋-轨道耦合(DSOC)的两分量BEC中的拓扑缺陷和自旋纹理. 结果表明, 对于不含自旋-轨道耦合的初始混合凝聚体, 旋转频率的增强可导致系统的结构相变. 在各向同性DSOC存在的情况下, 该系统维持涡流对、Anderson-Toulouse无核涡流、圆形涡旋片和组合涡旋结构. 特别地, 当旋转频率固定在径向势阱频率之上时, 强DSOC导致了由多层可见涡旋链、隐藏涡旋链和隐藏巨涡旋流组成的特殊拓扑结构. Li和Liu[107]利用Gross-Pitaevskii方程研究了自旋-轨道耦合作用和梯度磁场对旋转两分量BEC基态的影响. 研究结果表明, 在梯度磁场中, 随着自旋-轨道耦合强度增大, 基态结构由skyrmion格子逐渐过渡为沿着对角线方向排列的skyrmion列. 当自旋-轨道耦合强度和旋转频率都小的情况, 磁场梯度的增强可导致基态由平面波相转变为half-skyrmion; 当自旋-轨道耦合强度和旋转频率都大的情况, 梯度磁场可诱导hidden涡旋的产生. 图10为不同自旋-轨道耦合强度下梯度磁场中两分量87RbBEC基态粒子数密度分布和相位分布.
图 10 不同自旋-轨道耦合强度下梯度磁场中两分量87RbBEC基态粒子数密度分布(第1、2列)和相位分布(第3、4列)[107] (a)?(d)的${\tilde {\rm{\kappa}} }$值分别为0, 0.2, 0.8, 2
Figure10. Particle number densities (the first and second columns) and phase distributions (the third and fourth columns) of ground state of the two-component BEC of 87Rb for the different spin-orbit coupling strengths: the parameters of ${\tilde {\rm{\kappa}} }$ in (a)?(d) are 0, 0.2, 0.8, 2, respectively[107].

2
4.2.三分量旋量玻色-爱因斯坦凝聚体中的拓扑缺陷研究
-->对于自旋F = 1的旋量三分量BEC, 原子可能占据的塞曼态有$\left| {1, 1} \right\rangle $, $\left| {1, 0} \right\rangle $, $\left| {1, -1} \right\rangle $三种, 体系内会出现两种类型基态相——磁相和晶列相, 依赖于自旋无关相互作用和自旋相关相互作用.
实验上2001年Raghavan等[108]通过静态直流磁场在三分量旋量BEC中产生了暗孤子和涡旋结构. 2002年Ogawa等[109]采用四极磁场将自旋1BEC限制在Ioffe-Pritchard势阱中, 发现在凝聚体中产生了涡旋. 2006年Itin等[110]研究了自旋1BEC在一对载流线和偏置磁场控制的双磁阱中产生涡旋的几种机制, 图11为凝聚体快速分裂过程中涡旋的动力学形成, 显示动态涡旋存在于凝聚体的所有分量中, 在ψ–1分量中占99%以上, 在ψ0分量中动态涡旋和拓扑涡旋共存. 2014年, Ray等[80,81]在自旋1BEC中借助梯度磁场, 观察到了狄拉克磁单极子, 图12为狄拉克磁单极子的实验产生过程, 显示了凝聚体中不同自旋组分中的粒子密度随BZ, f的减小而减小.
图 11 涡旋的动力学形成[110]. 涡旋形成于凝聚体的所有分量中, 在ψ–1分量中占99%以上, 在ψ0分量中动态涡旋和拓扑涡旋共存
Figure11. Dynamical formation of vortices: vortices are formed in all components, more than 99% of total population is in ψ1 component. In the ψ0 component, dynamical and topological vortices coexist[110].

图 12 狄拉克磁单极子的实验产生[80] (a)?(f)每一行都包含单个凝聚体的图像. 最左边的列显示了三种自旋状态$\left\{ {\left| 1 \right\rangle, \left| 0 \right\rangle, \left| { - 1} \right\rangle } \right\}$沿水平轴的柱状密度彩色图像; 最右边三列显示沿纵轴拍摄的图像
Figure12. Experimental creation of Dirac monopoles. Each row (a)?(f) contains images of an individual condensate. The leftmost column shows colour composite images of the column densities taken along the horizontal axis for the three spin states $\left\{ {\left| 1 \right\rangle, \left| 0 \right\rangle, \left| { - 1} \right\rangle } \right\}$; The rightmost three columns show images taken along the vertical axis[80].

理论上2002年Isoshima和Machida[111]在旋转外势下, 采用推广的Bogoliubov理论研究了自旋1BEC中各种轴对称涡旋的稳定结构. 同年Mizushima等[112]通过求解广义Gross-Pitaevskii方程研究了旋转条件下铁磁自旋1BEC中的Mermin-Ho和Anderson-Toulouse无核涡旋, 得到在铁磁情况下, Mermin-Ho涡旋是稳定的. 2006年Saito等[113,114], 以及2009年Turner[115]研究了铁磁自旋1BEC中的极核涡旋, 发现系统核心由mF = 0原子填充, 围绕涡旋有mF = 1分量的质量流, 以及mF = –1分量的相等但自旋相反的质量流, 这就产生了净自旋流, 但没有净质量流. 2007年Mottonen等[116]研究了非旋转条件下在三维抛物势与Ioffe-Pritchard场中的自旋1铁磁态旋量BEC中无核涡旋的能量和动力学稳定性.
2008年Ji等[117]研究了自旋1钠原子BEC中半量子涡旋的动态生成. 模拟结果表明, 在外加脉冲磁捕获势的情况下, 旋转光势阱中可以同时产生独立的半量子涡旋和涡旋晶格, 同时还发现, 一个明显的周期调制的自旋密度波空间结构总是嵌入在方形半量子涡旋晶格中. 2012年Liu和Liu[118]研究了旋转和快速淬火的自旋1BEC中由自旋-轨道耦合引起的半skyrmion激发, 给出了描述半skyrmion自旋矢量的三个表达式. 结果表明, 半skyrmion激发依赖于自旋-轨道耦合和旋转的结合, 当自旋-轨道耦合强度和旋转频率均大于某些临界值时, 半skyrmion由一个或几个圆包围中心, 形成径向晶格, 即使在强铁磁(或反铁磁)凝聚体中也会发生. 图13为旋转频率对23Na旋量BEC自旋纹理的影响.
图 13 旋转频率对23Na旋量BEC的影响[118], 其中${\mu _{j, 0}}\left( {j = 0, \pm 1} \right) = 3.6\;\hbar {\rm{\omega }}$, ${\text{μ}} = 25\;\hbar {\rm{\omega }}$, κx = κy = κz = 1, ${a_0} = 50\;{a_{\rm{B}}}$, and a2 = 55 aB (a) Ω = 0; (b) Ω = 0.2 ω; (c) Ω = 0.5 ω. 第四列显示了相应的自旋纹理和涡旋的位置
Figure13. The effect of rotation frequency for spinor BEC of 23Na with ${\mu _{j, 0}}\left( {j = 0, \pm 1} \right) = 3.6\;\hbar {\rm{\omega }}$, ${\rm{\mu }} = 25\;\hbar {\rm{\omega }}$, κx = κy = κz = 1, ${a_0} = 50\;{a_{\rm{B}}}$, and a2 = 55 aB: (a) Ω = 0; (b) Ω = 0.2 ω; (c) Ω = 0.5 ω. The fourth column shows the corresponding spin textures and the positions of the vortices[118].

2013年Liu等[119]利用阻尼映射Gross-Pitaevkii方程, 研究了二维体系中自旋-轨道耦合的23Na自旋1BEC中的涡旋斑图, 研究发现较弱的自旋-轨道耦合可以完全破坏不考虑自旋-轨道耦合情况下出现的周期性涡旋晶格; 在自旋-轨道耦合较强的情况下, 各自旋态的涡旋易形成涡旋组, 并绕凝聚体中心形成花瓣状涡旋斑图. 2014年Song等[120]利用精确对角化和平均场理论研究了弱相互作用的自旋-轨道耦合自旋1玻色气体在外谐波势阱中的碎裂问题, 研究发现这种碎裂倾向源于总角动量守恒, 且受自旋-轨道耦合强度和自旋相关相互作用的影响. Lovegrove等[121]将铁磁态无核涡旋通过相位植入法, 在极化态凝聚体中得到了混合态下稳定的无核涡旋. 2015年Zhao等[122]在自旋1BEC中解析得到了两种不均匀的自旋畴构型, 它们分别由正二次塞曼效应和负二次塞曼效应所致. 分析表明, 二次塞曼效应可以诱导自旋畴的动态相变, 其符号可以影响自旋模式的拓扑结构. 2016年, Gautam和Adhikari[123]对零磁化强度的自旋1和自旋2旋量BEC中的分数涡旋进行分类, 并利用精确的数值解和拉格朗日变分近似研究了准二维旋量BEC中涡旋的静力学和动力学性质.
2017年Liu等[124]通过虚实演化方法研究了具有面内四极磁场自旋1的旋量BEC的基态结构. 研究发现, 面内四极磁场和旋转双重作用可导致中央Mermin-Ho涡旋的产生; 随着磁场梯度增强, Mermin-Ho涡旋周围环绕的涡旋趋向对称化排布; 在四极磁场下, 密度相互作用和自旋交换相互作用作为体系的调控参数, 可以控制Mermin-Ho涡旋周围的涡旋数目. Li等[125]研究了铁磁自旋1BEC中由自旋-轨道耦合引起具有极核旋涡的狄拉克磁单极子, 随着自旋-轨道耦合强度的增加, 具有极核涡旋的狄拉克磁单极子可以转化为正方形晶格. 在自旋-轨道耦合确定的情况下, 增大相互作用强度可引起从具有极核涡旋的狄拉克磁单极子向具有Mermin-Ho涡旋的循环相变. 图14显示了具有Mermin-Ho涡旋的磁单极子的形成.
图 14 具有Mermin-Ho涡旋的磁单极子[125] (a)等值面的粒子数密度; (b)粒子数密度等深线段(y ≤ 0), 节点线(Dirac线)的位置用红色箭头突出显示; (c) z=0平面上的位相分布. 单涡旋(mF = 0)和双涡旋(mF = –1)具有相同的环流, 由红圈突出显示
Figure14. The monopoles with the Mermin-Ho vortex: (a) Isosurface of particle densities; (b) segments of isosurface of particle densities (y ≤ 0). the position of the nodal line (Dirac string) is highlighted by the red arrow; (c) phase distributions in the z = 0 planes. the single vortex (mF = 0) and double vortex (mF = –1) have the same circulations, as highlighted by the red circles[125].

1998年Ketterle研究组首次在自旋为1的23Na原子系统中实现了BEC, 为冷原子物理开辟了旋量BEC的研究领域. 尤其是在超冷原子气体中人造自旋-轨道耦合的实验实现, 为研究拓扑量子态提供了一个理想的实验平台. 本文综述了旋量BEC的实验和理论研究, 旋量BEC中产生的拓扑缺陷的种类, 如自旋畴壁、涡旋、磁单极子、skymion、扭结, 着重介绍了两分量和三分量旋量BEC中拓扑缺陷的研究进展. 今后的工作可以推广到具有高自旋BEC体系, 以及不同自旋-轨道耦合形式的玻色气体等, 如Rashba自旋-轨道耦合和旋转势作用下的铁磁自旋2BEC中新奇的拓扑结构研究. 此外, 在冷原子平台上研究具有长程相互作用的拓扑缺陷也是未来的一个方向, 以及通过研究自旋-轨道耦合旋量BEC的动力学行为, 并从非平衡过程中来观察一些拓扑结构, 如冯卡门涡街和量子扭结, 也是非常有意义的工作.
相关话题/结构 实验 系统 空间 涡旋

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 离散可积系统: 多维相容性
    摘要:对比已有完善而系统理论的微分方程领域,差分方程理论尚处于发展之中.近年来离散可积理论的进展,带来了差分方程理论的革命.多维相容性是伴随离散可积系统研究出现的新的概念,作为对离散可积性的一种理解,提供了构造离散可积系统的B?cklund变换、Lax对和精确解的工具.本文旨在综述多维相容性的概念及 ...
    本站小编 Free考研考试 2021-12-29
  • 可积谐振系统中的极端波事件研究进展
    摘要:从微观角度上讲,单个极端异常波事件可视为可积模型方程的时空局域有理函数解.本文主要讨论了三类典型的可积谐振相互作用模型(即长波短波谐振方程,三波谐振相互作用方程,非线性薛定谔和麦克斯韦-布洛赫方程)的基阶Peregrine异常波解及其相关研究进展;明确指出了这些基阶异常波解形式具有普适性,可推 ...
    本站小编 Free考研考试 2021-12-29
  • 可积系统多孤子解的全反演对称表达式
    摘要:多孤子解是非线性数学物理系统的基本激发模式.文献中存在各种类型的表达式,如广田(Hirota)形式,朗斯基(Wronskian)或双朗斯基形式和法夫(Phaffian)形式.最近在多地系统的研究中,我们发现使用一种全新但等价的形式具有极为简洁和方便的优点.本文主要综述多种类型可积非线性系统的多 ...
    本站小编 Free考研考试 2021-12-29
  • 二维材料/铁电异质结构的研究进展
    摘要:二维材料是一类具有原子层厚度的层状材料,拥有独特的电学、磁学、光学和力学性能.以石墨烯和过渡金属硫族化合物为代表的二维材料展现出迁移率高、能带可调、可见光透过率高等特点,是近年来微纳科学领域的前沿热点.将二维材料与各种功能材料,如SiO2绝缘体、半导体、金属、有机化合物等结合,可以深化和拓宽二 ...
    本站小编 Free考研考试 2021-12-29
  • 基于密度泛函理论的外电场下C<sub>5</sub>F<sub>10</sub>O的结构及其激发特性
    摘要:采用密度泛函(DFT)B3LYP/6-311g(d)对C5F10O分子进行几何结构优化,研究外加电场(0-0.03a.u.,1a.u.=5.142×1011V/m)对分子的几何结构、能量、前线轨道能级、红外光谱的影响.在相同基组下,采用TD-DFT方法计算和分析C5F10O的轨道成分和激发特性 ...
    本站小编 Free考研考试 2021-12-29
  • 彗差和球差对涡旋光束斜程传输特性的影响
    摘要:利用涡旋光束作为空间光通信载波可以大大提高数据传输的容量,因此,研究涡旋光束在大气湍流中的传输具有重要意义.涡旋光束在大气湍流中传输时会产生光束漂移,进而影响通信系统的性能.本文基于多相位屏和傅里叶变换的方法,研究了带有彗差和球差的涡旋光束在大气湍流中传输时的光束漂移特性.结果表明,涡旋光束在 ...
    本站小编 Free考研考试 2021-12-29
  • 具有内参考热补偿功能的三层膜结构微球腔折射率传感器
    摘要:光学微腔在高灵敏度传感中有着重要的应用前景,而在传感中热漂移是制约其走向实用的重要因素.本文提出了一种镀有三层膜结构的微球腔,可以在实现高灵敏度折射率传感的同时,具备内参考热补偿功能.该结构由内到外分别涂覆折射率为高、低、高的薄膜,内外两高折射率层可以分别支持各自的回音壁模式,称之为内层模式和 ...
    本站小编 Free考研考试 2021-12-29
  • 准Λ型四能级系统选择反射光谱
    摘要:运用密度矩阵理论研究了气固界面准$\Lambda$型四能级原子系统的非线性选择反射光谱.基于刘维尔方程给出了一阶近似条件下探测光场的解析式.在探测场为弱场时,分析了信号场拉比频率、失谐量和耦合场失谐量对反射光谱线型的影响.数值模拟表明:信号场参与产生的选择反射峰线宽可以利用信号场拉比频率进行调 ...
    本站小编 Free考研考试 2021-12-29
  • He离子辐照对石墨烯微观结构及电学性能的影响
    摘要:本文采用5.4keV不同剂量的He离子辐照单层石墨烯,利用X射线光电子能谱(XPS)、拉曼光谱(Raman)和半导体参数分析仪表征辐照前后石墨烯微观结构和电学性能变化.研究结果表明:随着辐照剂量增大,单层石墨烯的缺陷密度逐渐增加,当辐照剂量增至1.6×1013He+/cm2,石墨烯开始由纳米晶 ...
    本站小编 Free考研考试 2021-12-29
  • 隔板对流系统的热流特性及热量输入与传递特性
    摘要:采用DNS方法对隔板对流装置进行模拟计算,研究系统中热流特性以及热量输入与传递特性.讨论了热流的纵向和横向输运特性,在此基础上对传热通道和狭缝区域的热通量以及对应底板外界输入热通量进行了定量化分析.研究结果表明,通道中低温流体向下冲击底板而后转入水平运动,流入狭缝区域并不断被加热,只进行水平的 ...
    本站小编 Free考研考试 2021-12-29