冷原子物理中的一维少体问题*

刘彦霞 张云波

1）(山西大学,理论物理研究所,太原 030006）

2）(中国科学院物理研究所,北京 100190)

作为构成量子多体系统的基本单元,一维少体系统的研究不仅可以在理论上为多体系统的量子关联及动力学等性质提供更为基本的理解,也可以为实验上制备多体系统提供更加方便和功能更加全面的方法.本文回顾了冷原子物理中一维少体系统最新的实验和理论进展.首先介绍了少体实验中实现的谐振子势阱中确定原子数的精确制备,亚稳态势阱和双阱系统中原子的隧穿,以及强相互作用下等效自旋链的实验结果.然后深度解析了理论研究方面,特别是基于精确可解模型的一些重要结果,包括亚稳态势阱中相互作用原子的隧穿概率,以及相应实验上常见势阱的能谱分析、密度分布、隧穿动力学以及强相互作用极限下的有效自旋链模型等.

1 引 言

少体物理几乎在所有的物理分支中都起着重要作用,比如天体力学中的三体开普勒问题以及量子力学中把氦原子中的核子看做是质点的三体问题,一般情况下三体问题很难退化到一个有效的单体问题.这里简单介绍一些囚禁在势阱中的玻色子或费米子的少体问题,这些都是研究少体物理的理想模型.这些系统中有很多是可以解析处理的,尽管存在很强的相互作用和很多的自由度,但是在某些情况下,可以转换成单体问题来求解.

如今,简并的玻色和费米气体在世界各地的实验室都可以产生.大多数实验小组都以碱金属原子作为研究对象,如 Li,Na,K,Rb,Cs,这些原子的共同特点是最外层只有一个价电子.在超冷温度下,当两体相互作用的有效力程远小于德布罗意波长时,原子间的相互作用主要是以各向同性的s波碰撞为主.s波碰撞的效应可以通过s波散射长度描述,在很多情况下成为用来描述系统的唯一微观参数[1,2].一个经典的例子就是平均场Gross-Pitaevskii方程[3−6]只依赖于粒子数 N 和s波散射长度.Gross-Pitaevskii方程将多体问题转换成一个有效的单体问题,可以非常准确地描述弱相互作用稀薄玻色凝聚体的特性.除了碱金属原子,实验上也可以实现价电子结构更复杂的原子的凝聚,如 Ca,He,Cr,Yb,Sr.例如两个 Cr原子间的碰撞是各向异性和长程的[7,8].这些特性会引起超出s波的物理,两原子间的相互作用除了依赖s波散射长度外还依赖偶极相互作用的强度.

实验上囚禁少体系统可以通过所谓的微势阱装载少数个粒子来实现[9].一系列的少体系统可以通过将多个原子囚禁在很深的光晶格中来实现,其中不同格点间的隧穿被抑制,近邻相互作用可以忽略[10,11].对于冷原子系统,原子间的相互作用可以通过外磁场在Feshbach共振附近调节[12].超冷碱金属原子间的Feshbach共振现象可以用两通道模型解释.两个原子最外面的两个电子的相互作用可以由一个比较浅的电子自旋三重态和一个比较深的单重态的Born-Oppenheimer势能曲线来描述.超精细相互作用将单态和三重态的势能耦合,形成一个闭通道和一个开通道.当闭通道的束缚能量和开通道入口处的能量相同时就发生Feshbach共振,散射长度趋于发散,1998年首次在超冷原子气体中观察到Feshbach共振[13−15],为理论和实验对比铺平了道路.

一个特别有趣的少体现象是Efimov效应.早在 1970年,俄罗斯核物理学家 Efimov[16,17]就发现在无穷大s-波散射长度时,自由空间的两体系统只存在一个零能束缚态,没有其他弱的束缚态.但是当加入第三个粒子以后,就会出现令人惊讶的结果,相应的三体系统会出现无穷多个弱的束缚三体态,且束缚能有分立标度不变性[18,19](这里假设只有两体相互作用).Efimov效应激发了很多理论和实验的研究[20],直到现在仍然是少体物理中的重要研究课题.在氦的三聚物、核物理、还有其他系统中都有三体Efimov效应和关联的Efimov共振的研究.然而,它的第一个实验验证是2006年由Kraemer等[21]在铯原子气体中完成的.从更广的范围看,Efimov态是 Borromean环的一种,这一类物态的标志特征是只有三个粒子会形成束缚态环链,而其中任意两个都不互相束缚.这是一个纯粹由量子力学导致的效应,有丰富而深刻的物理内涵.自从 Efimov效应提出以后,推广到大系统已经有很多进展,但这仍然是一个开放问题.在自由空间或者有外势的情况下N体系统(N>3）的能谱如何依赖于三体参数,四体、五体和更多体系统的研究可能会给出答案.

本文主要回顾最近超冷原子实验在准一维少体系统中的突破性进展及相关的理论研究[22,23].海德堡Zürn等首先将两个费米子囚禁在光势阱中,研究了强相互作用下费米子的费米化,通过施加一个磁场梯度,在排斥相互作用下研究了单粒子的隧穿[24],在吸引相互作用下研究了配对粒子的隧穿[25].然后通过精确地逐个增加原子个数,测量单个杂质粒子与背景气体的相互作用可以用来研究从少体到多体的跨越,观察到了费米海的形成[26].该实验组接着实现了两个原子在双势阱间的量子隧穿,得到相互作用强度对双势阱的本征态以及占据概率的影响[27],在一维的谐振子势阱中准确制备由自旋1/2费米子构成的反铁磁海森堡自旋链[28].Kaufman等[29]用光镊子将单个原子冷却到它的振动基态,随后扩展实验技术将两个原子冷却到双阱的基态,通过两个粒子的干涉观察粒子的统计效应[30].

2 一维冷原子少体实验

目前超冷气体的研究大部分关注的是多体系统的性质,在热力学极限下粒子数的涨落和散粒噪声一样,不会明显改变热力学量,例如温度、压强等.与此相反,少体系统非常依赖于粒子的数目.目前的实验技术已经可以很精确地控制少体系统中粒子的个数.海德堡Jochim实验小组从2011年开始对一维冷原子物理中的少体问题做了一系列关键实验,包括相互作用原子在亚稳态和双阱模型中的隧穿[24−28],从少体到多体的转变,以及一维自旋链的量子磁学.他们用光偶极阱和磁场梯度将少数6Li费米原子制备到基态,精确度可以达到93±2%.粒子囚禁在纵横比为1:10的三维雪茄状的势阱中,各方向可近似认为是谐振子势阱,此系统可以作为一维系统来处理.文献[31]中给出的解析解可以很好地描述这种势阱中两个相互作用原子的基态.一维光势阱加上梯度磁场可以将势阱的一边打开而形成一个亚稳态势阱,是研究隧穿动力学的重要模型之一.原子间的相互作用可以通过Feshbach共振[12]和束缚诱导共振[32]来调节.一维的耦合常数g可以通过三维散射长度和横向谐振子势的特征长度来确定,这里是普朗克常数,是两个质量为m的原子的约化质量.耦合常数为

2.1 原子的费米化及单粒子隧穿

当相互作用强度趋于无穷时,两个可区分的粒子的能量和波函数的模平方与两个无相互作用的全同费米子相同,这被称为原子的费米化.实验上通过直接比较强相互作用极限的两种可区分费米子和相同势阱中的两个全同费米子来观察这种现象.文献[24]将两个6Li费米原子制备到最低的两个Zeeman子能级的超精细态,分别标记为和.对于可区分粒子系统和全同系统分别用和和表示.由于费米子的反对称性,全同系统感受不到s波相互作用,对于可区分粒子系统的相互作用则可以通过磁Feshbach共振来调节.实验上可以测量其中一个原子隧穿出势阱的时间随相互作用强度g的变化规律,通过拟合留在阱中的平均粒子数可定义一个隧穿时间常数.实验发现随着磁场的增加，两个可区分原子的减少了两个数量级,而全同费米子的保持为常数,二者在 CIR处相交.对6Li原子,该 CIR 共振发生在磁场为 (783.4±0.4）G()处.也就是说在相互作用强度时,态和态的能量,波函数密度分布及隧穿时间常数都是相同的.文献 [33]通过WKB(Wenzel-Kramers-Brillouin)近似给出了隧穿的准粒子理论,很好地解释了实验的数据点,揭示了波函数在束缚诱导共振点处的费米化特点.

2.2 原子的配对隧穿

原子费米化的研究是将两个可区分的原子制备到排斥相互作用区域,该区域系统中的粒子只会出现单粒子隧穿.当原子处在吸引相互作用区域时,会发生两个原子同时隧穿的现象[25].测量过程如下:首先将两个原子制备到基态,一个原子处在态,另一个处在态,然后突然加上磁场梯度,使势阱的一边倾斜形成势垒,这时粒子会跑出去,在某个固定的时间迅速去掉磁场梯度,将势阱恢复,然后测量势阱里面的原子个数,重复多次,计算某一时刻的势阱出现两个原子、单个原子、以及零个原子的概率.通过这些概率可以得到平均粒子数随时间的演化.

为了描述两个粒子的配对隧穿,实验中用一个简单的模型来描述,如图1所示.对于两个原子的隧穿有两个过程:一个过程是配对隧穿,即两个原子同时离开势阱,发生这个过程的速率记为.另一个过程是原子按次序离开势阱,第一个粒子先离开势阱,第二个粒子处在势阱的非微扰基态,这个过程第一个粒子隧穿的速率记成.这个速率主要由隧穿势垒的高度决定,同时依赖于两个原子的相互作用能.对于第二个原子,没有相互作用的影响,因此离开势阱的速率是无相互作用系统的.考虑实验上有限的准备精确度以及通过改变磁场调节相互作用强度会影响原子自旋的取向,隧穿速率可以修正成为一个依赖于自旋的结果.实验结果表明在相互作用强度g>–0.59时是没有配对隧穿的.在强吸引相互作用区域 g<–0.64,粒子配对隧穿占主导地位,单粒子的隧穿速率很小,可以认为两个粒子绑在一起形成了束缚态,相互作用可看做是其内部的自由度.文献[33,34]通过WKB近似理论计算给出了隧穿率和其随相互作用强度的变化.

图1 双原子的两种隧穿过程:单原子次序隧穿及两原子配对隧穿.本图摘自参考文献[25]Fig.1.The loss processes include two tunneling processes of two atoms out of a metastable potential:subsequent singleparticle tunneling and direct pair tunneling (Reproduced with permission from Ref.[25]).

2.3 双阱中的两个费米子

实验上双阱可以由两束激光通过高分辨率的物镜聚焦产生,通过声光偏转器独立地控制两束激光的强度和位置,可以调节双阱的隧穿耦合系数J和两阱之间的偏差.双阱是构成光晶格的基本单元,文献[27]介绍了半满填充时Fermi-Hubbard模型的基本构建块的实现,即在双阱中态和态原子形成的自旋单态.在Hubbard模型体系,两个粒子的空间波函数可以由基矢展开,其中和分别表示单粒子在左阱和右阱的基态.通过对角化对称双阱在双模近似下的哈密顿量,可以得到四个本征态,其本征能量分别是,这里c 是依赖隧穿能J和相互作用强度U的,并且这两个态都是满足交换对称的.文献[27]给出了实验结果:基态a的双占据概率随着排斥相互作用的增加而增加,第二激发态c则相反,双占据几率随着相互作用的增加而减小.b和d这两个态的占据概率则不依赖于相互作用强度.d是交换反对称的态,其占据概率不随相互作用强度改变.但是b是交换对称的态,实验中并没有给出这个态的占据概率随相互作用强度的改变.另一方面,从这几个态的宇称对称性来看,a和c是偶宇称态,b和d是奇宇称态.双模近似下b是一个两粒子的NOON态,是量子力学两体的最大纠缠态.

为了研究两个粒子在双阱中的隧穿动力学.实验上首先将两个原子制备到态,然后迅速将势垒降低到可以使原子在双阱中来回隧穿.为了观察动力学的结果,在演化的不同时刻,迅速增加势垒的高度使原子在双阱中的空间分布锁定,然后通过荧光谱测量每个阱中原子的个数,重复测量多次得到每一时刻的占据概率.实验可以给出不同参数下占据概率随时间的振荡图像,也可以将原子制备到系统的本征态,如基态、激发态,探测两个粒子在不同阱的概率和两个粒子在同一个阱的概率,得到相互作用强度对双势阱的本征态以及占据概率的影响.

2.4 少体到多体: 观察费米海的形成

知道一个物理体系什么时候可以具有宏观特性并且可以用多体理论很好地描述是非常困难的.一维冷原子的少体实验[26]通过精确控制单杂质系统中全同费米子的数目来研究从少体到多体的跨越.首先制备一个 N+1 个粒子的系统,其中 1 个杂质粒子和N个囚禁在细长光偶极阱基态的全同粒子之间存在排斥相互作用.实验上将超冷费米子原子制备到两个精细态来实现这样的系统,这里杂质粒子处在态,全同粒子都处在态.为了探测这个系统,实验上主要测量杂质和全同粒子间的相互作用能随多数粒子数的变化.相互作用能的测量可以通过RF(radio frequency)脉冲改变杂质粒子的内态来实现.如果没有全同粒子的存在,转变发生在频率,对应杂质粒子初末态超精细能级差.对于有N个全同粒子存在的情况,杂质粒子和全同粒子之间的相互作用引起一个的频率移动.对于给定的相互作用强度,相互作用能为.

每个全同粒子的添加都会增加与杂质粒子相互作用的原子数目,因此随着粒子数的增加,√相互作用能是增加的.对于弱相互作用,满足.当时,是发散的.因此,通过自然能量标度重新定义相互作用能,无量纲的相互作用能为,这里是全同原子的费米能.由于只考虑相互作用能,因此忽略零点动能,这意味着对于谐振子的情况.为了补偿由于增加粒子引起的密度改变,引入无量纲的排斥相互作用参数>0.为了判断系统是否达到多体极限,实验的数据要和两种极限情况比较,一种是N=1 的情况,另一种是的情况.前一种情况的相互作用能可以通过两相互作用原子的精确解得到,后一种情况是多粒子极限,即单个杂质粒子浸入到无穷多个粒子形成的费米海中.对于=0,相互作用能为0.对于,系统达到费米化极限,杂质粒子和N个全同费米子相互作用的能量与N+1个没有相互作用的全同粒子的能量相同,因此这种情况的相互作用能与费米能量是相同的.N=1 和两种极限情况在时结果是一致的.比较理论和实验的结果,可以发现时的相互作用能和无穷多个粒子情况符合很好,也就是说,就已经体现了多体的特性.

2.5 一维势阱中的反铁磁Heisenberg自旋链

实验[28]报导了个自旋向上的粒子和个自旋向下的粒子在一维势阱中如何实现Heisenberg自旋链.实验中研究的系统有.这些体系的基态在强相互作用区域都表现出反铁磁态的特征.实验上通过两个独立的测量来判断系统是否处于反铁磁态:利用隧穿技术探测自旋链最外边粒子自旋的方向,以及利用波函数在单粒子能级的投影探测自旋向下粒子在体系中的空间波函数.

在准一维的系统中可以通过磁Feshbach共振和束缚诱导共振,将不同组分间的一维相互作用强度调到共振区域,即费米化区域,此时相同自旋组分的散射是被禁止的.在费米化区域,原子之间是不可穿透的,因此原子在势纵轴方向的排序是固定,这样就可以通过隧穿确定最外边原子的自旋方向.倾斜势阱,最外边的原子会隧穿出去.通过探测势阱中剩余自旋向上原子的个数,来判断隧穿出去的原子的自旋,重复多次测量,计算概率,与理论结果比较.例如,对于系统 (2,1）,其基态为,自旋向下粒子隧穿出去的概率为.在实验中由于要加磁场梯度,因此实验结果要考虑外加磁场梯度使外势不对称引起的变化.通过以上的测量计算就可以判断是否处在反铁磁态.对于强相互作用系统,相对坐标在杂质处会产生一个尖峰,这个尖峰会导致更高能量的能级占据.一个态空间波函数越对称,尖峰就会越多.因此单粒子能级占据数的分布直接反应了系统自旋的结构.不引入外势的情况下,两分量间的强相互作用就可以使自旋链达到稳定,这使得在非格点模型中制备自旋链成为可能,也为磁学性质的研究提供了一个新的平台.

3 一维冷原子少体可解模型

量子力学中只有屈指可数的几个问题是可以严格求解的,其中应用最广的例子之一是谐振子势阱中的单个粒子.考虑相互作用的两体问题可精确求解的则更少,特别地,氢原子中由于库仑相互作用只和相对坐标有关,引入相对坐标和质心坐标可将动能部分分离变量,描述质心运动的方程是一个自由粒子的能量本征方程,而描述相对运动部分的方程形式上与单体波动方程完全一样,两体问题于是简化为单体问题.实验上使用激光冷却和蒸发冷却技术实现的简并原子气体需要将原子囚禁于外势场中,这些囚禁系统中的相互作用所起的作用与自由空间中的超冷碰撞不同.在自由空间中,通过微分散射截面和碰撞相移因子来探测两原子间的相互作用,而在存在囚禁势的情况下这些物理量都没有了意义,因为囚禁系统中无法构造入射态和出射态的渐进行为.相反,在囚禁势中相互有作用的多原子系统的特征在于其离散能谱以及能级间的跃迁矩阵元.下面简单介绍几种一维囚禁势阱中严格可解的相互作用少原子体系的离散能谱及相关的动力学结果,该体系的哈密顿量为

3.1 Bethe Ansatz方法简介

对于多组分系统,粒子具有内部自旋自由度,S 矩阵通常是算符,满足Yang-Baxter方程

该方程是Yang[39]和Baxter[40]在研究一维势相互作用费米气体模型和二维格子统计模型时分别提出的.其含义是系统从初态到末态有两种途径,这两种散射过程是等价的.(4）式是周期性边界条件系统可积的条件.根据Yang-Baxter方程,得到转移矩阵,这里Lax算符形式与散射矩阵相同:

3.2 谐振子势阱中的两粒子问题

1998年Busch等[31]给出了两个冷原子在谐振子势中的精确波函数和解析解.为了避免处理真实原子间相互作用的困难,在稀薄原子气体中通常采用零程的点状势将原子间相互作用近似为接触相互作用,一维系统中只含有简单的项.质量为m的两个原子在频率为的一维谐振子势中的运动由(2）式中描述的哈密顿量,其中

图2 三种相互作用的两原子系统的能量E与相互作用强度g的关系:(a)谐振子势阱中的两原子,黑实线和红虚线分别表示相对运动奇宇称和偶宇称波函数所对应的能量;(b)周期边界条件下两异核原子,黑实线和红虚线分别表示质量相等两原子和质量比为2.175情况下对应的能量;(c)双势阱中相互作用两原子的奇宇称态,一维无限深方势阱中心为一个强度为d=0.5的 势垒劈开,黑实线和红虚线分别表示准动量为实数的原子本√征态和准动量为复数的分子态对应的能量.这里能量的单位分别是,相互作用强度g的单位分别是Fig.2.Energy spectrum of three types of interacting two-atom system:(a)Two atoms in a harmonic oscillator potential.Black solid lines and red dashed lines are odd parity and even parity energy level respectively;(b)two heteronuclear atoms in a ring trap.Black solid lines and red dashed lines are energy levels for equal mass and mass ratio =2.175respectively;(c)two atoms in a -split hard-wall double well.Five lowest odd parity levels for barrier height d=0.5.Black solid lines and red dashed lines are the bound states for atoms with real-valued quasimomentum and the molecule states with complex-valued quasimomentum respectively.

Rubeni等[75]用变分法对一维谐振子势阱中两相互作用费米子的基态波函数进行了研究,但其试探波函数存在缺陷,即在边界处一阶导数不连续,Liu等[76]改进了该变分法方案得到了关于动量的一个限制条件,该条件类似于 Bethe ansatz方程,但是与边界条件无关,结合变分法和波函数及其一阶导数的连续性条件得到的基态能量和精确结果符合,给出了更为合理的基态波函数,可以作为两原子亚稳态隧穿的初始波函数.

3.3 周期边界条件下两异核原子的解析解

冷原子混合物体系是不同质量原子的混合,实验 上 已 经 成 功 制 备 了40K-6Li,40K-87Rb,6Li-23Na等异核系统[77–79].异核费米子的超流性质、配对机制以及异核分子的形成等新现象得到广泛研究[80–82].在异核系统中,一个新的自由度——质量比表征了原子混合物体系的特性,诱导了异核模型中新的物理效应.这里介绍一种精确求解周期边界条件下两体相互作用异核系统的精确求解方法[83],其关键在于较为普适的波函数假设.

具有接触相互作用的两个异核原子在周期边界条件下可以看作原子囚禁在一维环状的势阱中,在哈密顿量(2）式中取

异核系统的动量分布与全同粒子系统区别很大,因为系统不再有交换对称性.质量比为时的动量分布表明[83],对于全同粒子系统,每个粒子的动量分布相同,动量位于或者0的概率相同;对于异核系统,质量大的粒子以为中心,质量小的粒子则集中在0处,从而保证总体能量更低.

3.4 双原子在亚稳态势阱的隧穿

冷原子少体实验的准确制备使得在少体系统中研究隧穿问题成为可能.隧穿理论的两个主要模型包括亚稳态势阱和双势阱.在这两种势阱中研究隧穿已有很长的历史,最新进展包括隧穿的干涉效应以及共振隧穿等[24,25,84].亚稳态势阱的量子隧穿衰变是早已熟知的量子现象,如核衰变、量子成核等,这些模型都是在单粒子图像上建立的.这里主要介绍相互作用的双原子体系在亚稳态势阱的隧穿问题,哈密顿量中的外势场如下

图3 (a)衰减时间随磁场的变化.黑色虚线和实线分别是WKB和考虑准粒子波动方程的修正的结果.绿色线表示通过微扰论考虑非谐项贡献的结果;(b)相互作用强度和隧穿能量随磁场的变化.本图摘自参考文献[33]Fig.3.(a)Decay time vs magnetic field B.The points with error bars are the experimental data[24],the dashed and solid lines are,respectively,the WKB and QPWF predictions. The green light gray lines include the perturbation theory correction to the tunneling energy ;(b) interaction strength (red gray curve) and tunneling energy(black and green light gray curve)vs B(Reproduced with permission from Ref.[33]).

3.5 双阱中少体玻色系统的基态性质及隧穿动力学

双阱是研究超导比特Josephson振荡[86]和玻色爱因斯坦凝聚中自囚禁现象[87–89]的典型模型.受实验启发,对于双阱的研究,目前大部分的工作都是基于双模近似[90,91],或者是数值求解[92,93,94],没有严格解的结果.Murphy等[95]研究中间带有势垒的一维简谐势阱中一对超冷玻色原子的基态和低能激发态性质.通过计算约化单粒子密度矩阵、动量分布和双粒子纠缠来完整表征基态性质,得到了无限排斥相互作用Tonks-Girardeau极限的几个解析表达式,并通过数值解处理了有限粒子间的相互作用,这里空间坐标的离散化是通过分立变量表象(DVR)做到的.双阱中的对相互作用形成了周期势中多体系统的基本构件,势垒的添加使得系统变得不平庸,由于中间势垒的存在,势垒两边的粒子可以发生隧穿,人们在双阱中发现了约瑟夫森振荡和自束缚现象等.文献[96]研究了束缚在中间带有势垒的无限深方势阱中的少体玻色气体的基态性质.利用玻色子与费米子的对应关系构建了Tonks-Girardeau气体的基态波函数,研究了势垒两边的关联特性以及这种特性和粒子数目宇称的关系,并且运用精确对角化方法计算了在任意势垒高度和原子间相互作用强度下系统基态的密度分布、占据数分布、动量分布等.计算表明动量次峰峰值随着势垒的升高而增大,随着原子间相互作用强度的增大而减小.该模型甚至可以被推广到玻色费米混合物系统[97],在动量分布中次峰的出现显示了势垒两边有着显著的关联.

Liu和Zhang[98]给出了双阱中的BA严格解,在严格解的基础上研究了其隧穿动力学问题.考虑N个质量为m的玻色子在中间加势垒的一维无限深方势阱中运动,(2）式的系统哈密顿量中

区别于(3）式中的BA波函数,这里引入额外的求和指标n,其物理意义为位于势垒左边的原子个数.用下标n来标记这一块(block)区域,(17）式中对n的求和包含所有可能的 N+1个块.另外表示粒子向右或向左运动,在势垒处粒子发生反射.因此对于固定的块区域n,开边界条件对于左边界和右边界是不同的,在左阱的粒子()只能感受到左边的边界条件,在右阱的粒子()只能感受到右边的边界条件0.于是存在两个反射矩阵和分别在左右边界上反射.

至于两个粒子的散射,注意到仅当相邻两个粒子的坐标位于双阱的同一侧时才能发生散射.由相互作用条件导致波函数相对坐标的一阶导数不连续,可以得到散射矩阵.根据位于中心处势垒的阶跃条件连接块区域和n,可以得到第n个粒子和势垒的散射矩阵.存在势垒时准动量满足的BA类型的方程为

准动量要同时满足(18）式和(19）两式,通过数值求解得到N=2即两粒子在双阱中的能谱图.图2(c)展示了d=0和2的能谱随着相互作用强度g的变化.在无势垒的情况下奇宇称和偶宇称的解析解都存在,除g=0外该方法只找到了奇宇称的解.对于 g>0 的情况,准动量的解都是实数,对于 g<0 的情况,准动量有两种解,实数解对应亚稳态,互为共轭的复数解对应分子态.

精确解模型得到的能谱可用来研究两个原子在双阱中的布局数及隧穿动力学.随着排斥相互作用的增加,双占据的概率减小,单占据概率则增加,强相互作用时双占据和单占据的概率都接近50%.将两个原子制备在一个无限高势垒的NOON态(数值上取 d=300),然后将势垒的强度突然降低到 d=0.5,原子在两阱间隧穿.图 4 给出了两个原子在同一个阱中的占据概率随时间的演化.随着排斥相互作用强度从弱到强,双占据概率振荡过程的振幅逐渐增强.强排斥相互作用使得振荡不再完整,双占据概率不能再100%回到初始NOON态.图4(b)展示了相互作用强度分别为g=1.5和20的两体密度分布在一个周期内的变化.

3.6 强相互作用系统的等效自旋链

强相互作用系统有很丰富的物理,但通常情况下求解是非常困难的.然而有一种情况除外,即相互作用趋于无穷大时.对于无穷大相互作用极限的情况,不论是玻色子还是费米子都表现出无自旋费米子的行为,每个粒子占据不同的位置,不可穿透,也就是说具有了费米化的特点,因此相互作用趋于无穷时的系统是容易处理的.对于强相互作用区域,相互作用强度为有限值的连续系统可以等效成自旋链模型,在实验上通过自旋1/2的费米子系统已经实现了海森堡自旋链.一个重要的区别在于外势的不同会导致格点间的耦合强度的非均匀性.比如对强相互作用极限下的自旋1/2系统,无穷大相互作用附近的微扰计算可以将连续模型的哈密顿量等效为一个自旋链模型

图4 (a)两个原子在同一个阱中的占据概率随时间的演化图.黑实线,蓝虚线,红点线分别表示不同的相互作用强度g=0.5,1.5,20.在 t=0 时刻势垒高度突然从 d=300 降到 d=0.5.(b)g=1.5 和 20 时一个周期里不同时刻的密度分布图,这里坐标 和 的单位是 .本图摘自参考文献 [98]Fig.4.(a)Tunneling dynamics of the occupation probability of finding both atoms in the same well for g=0.5(black solid line),g=1.5(blue dashed line),and g=20(red dotted line).The barrier is abruptly lowered from a height d=300to0.5at time t= 0.(b)the two-body density functions at different times t,for g= 1.5 and 20,respectively.Here the coordinates and are in units of (Reproduced with permission from Ref.[98]).

4 结 论

量子少体问题是当前研究的一个热点问题.少体系统的研究对于理解微观现象有很大的帮助,也为研究多体系统奠定了基础.随着冷原子实验迅速发展,实验上的精确操控,少体系统的研究成为目前关注的热点.本文简单介绍了实验和理论两方面的进展.由于目前冷原子实验的精确操控,少体实验可以准确地观察到一些重要的物理现象,比如强相互作用时的费米化,格点模型中的量子相变,以及少体到多体转变等等.这些基本物理现象在少体实验中的实现让我们看到了少体研究的重要性.对于少体的理论研究不同于多体物理,基于精确可解模型的一些重要结果,包括亚稳态势阱中相互作用原子的衰减率以及双势阱中的相互作用原子的隧穿动力学研究,使得该研究领域富有挑战性又对理解单体到多体系统的物理性质提供了过渡.

感谢陈星博士在关于异核系统精确解方面的讨论,以及昌茂林同学在论文撰写方面提供的帮助.