偶极玻色-爱因斯坦凝聚体中的各向异性耗散*

赵珊珊 贺丽 余增强†

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

2)(山西大学物理电子工程学院,太原 030006)

针对偶极相互作用的玻色-爱因斯坦凝聚体,解析计算了点状杂质沿平行极化轴和垂直极化轴运动的能量耗散率,证明了在超流临界速度更大的方向上耗散率也更高.该结论为最近在162Dy原子气体中观测到的实验现象提供了理论支持.对于一般的运动方向,给出了耗散率在高速极限下以及临界速度附近的渐近形式.结合数值计算的结果,论证了耗散率随方向角的变化总是表现出与临界速度一致的各向异性.

1 引 言

无阻尼运动是超流现象的基本特征之一.对于一个在超流体中运动的大质量杂质而言,仅当其速度大于超流体的朗道临界速度时才会产生能量耗散,临界速度 vc的取值由体系的元激发能谱决定[1,2].在各向异性超流体中,vc的大小与运动方向有关,不同方向上的耗散行为也会表现出差异.

超冷原子气体作为高度可调控的量子多体系统为超流性的研究提供了一个理想平台.实验上通常采用搅动激光束的方法模拟杂质的运动,并通过观察原子云的加热效应测定玻色-爱因斯坦凝聚体(BEC)或费米超流体的朗道临界速度[3−7].最近,德国Wenzel等[8]在具有偶极相互作用的162Dy原子BEC中首次观测到了各向异性的临界速度.实验上发现,相较于垂直极化轴的运动而言,平行极化轴方向的临界速度更大,耗散率也更高.

本文将通过理论计算对上述实验现象进行解释.利用线性响应理论和Bogoliubov近似,考察了在偶极BEC中点状杂质沿不同方向运动的能量耗散率.对于实验上观测的两个主轴方向,解析地推导出耗散率随约化速度变化的函数形式,从理论上证明了Wenzel等[8]发现的各向异性关系.对于偏离主轴的运动方向,我们给出了耗散率在高速极限下以及 vc附近的渐近形式,并结合数值计算的结果,论证了耗散率与临界速度随方向角的变化趋势是普遍一致的.

2 模型和理论框架

考虑零温下一个大质量的点状杂质以速度v在N 个磁性原子(如 Cr,Dy,Er等[9−11])形成的BEC中运动.针对实验的情形,设BEC原子的磁偶极矩被外场极化在z方向.杂质的运动方向与极化轴之间的夹角记为θv.采用接触势模型描述杂质与磁性原子之间的相互作用,其哈密顿量为表示相互作用强度,R0为杂质在t=0时刻的位置,ri为第i个原子的坐标.

在BEC原子的不断碰撞下,杂质的动能一旦出现减损,则意味着耗散的产生.根据线性响应理论,能量耗散率(即单位时间内的平均能量变化)可写为[12,13]

其中 Eq为BEC中动量为 q的准粒子激发能,q=|q|为动量的大小,n和m分别为玻色子的密度和质量,ℏ 为约化普朗克常数.由于磁性原子间同时具有各向同性的s波相互作用和各向异性的偶极相互作用,因此耗散率P不仅依赖于杂质的速度大小v,也依赖于运动方向的极角 θv.

偶极BEC的元激发能谱可以在Bogoliubov近似下解得[14,15]

其中 θq为动量q与极化轴之间的夹角,为长波极限下声子激发的声速.在平行极化轴方向上,声速取最大值在垂直极化轴方向上,声速取最小值这里 c0表示仅考虑s波相互作用时的声速,ϵdd是衡量偶极相互作用相对强度的无量纲参数.本文将限于讨论 c0＞0 且 0＜ ϵdd＜1的情形,以保证元激发能谱是稳定的.

由(1)式可知,产生耗散的必要条件是体系中出现满足关系式 Eq=q·v 的准粒子激发,这意味着杂质的运动速度必须大于BEC的临界速度 vc.与各向同性体系不同的是,在各向异性超流体中,vc并非简单地由运动方向的元激发能谱决定,朗道提出的临界速度公式应做适当的推广[13,16].对于偶极BEC,已有理论工作给出[13]

在平行和垂直极化轴方向上,vc分别达到最大和最小,其取值与运动方向的BEC声速相等.而在其他方向上,临界速度则小于运动方向的声速.

(1)式—(3)式构成了本文理论计算的出发点.以下为了对比不同方向的耗散行为,我们把P重新表达为约化速度 u ≡v/vc的函数.于是,在任意方向上,只有当 u ＞1 时,耗散率才不为 0.

3 主轴方向的能量耗散率

首先来考察实验观测的两个主轴方向.

当杂质沿平行极化轴方向运动时(θv=0 或π),(1)式中的被积函数仅与q和θq有关,代入(2)式中 Eq的表达式并完成积分,即可得到纵向耗散率为

显然,由于元激发能谱的各向异性,P‖与P⊥的形式也是不同的.易于验证,当偶极相互作用参数ϵdd趋于0时,(4)式和(5)式均回到各向同性体系的结果[17,18].

通过比较 P‖和P⊥的大小很容易解释实验上观察到的各向异性耗散现象[8].由声速关系c‖＞ c⊥可知,因此得到

(6)式表明,对于任意 u ＞1 的 约化速度,在 vc更大的主轴方向上(即平行极化轴方向上)能量耗散率更高.这个关系式为 Wenzel等[8]最近在162Dy 原子BEC中观测到的实验现象提供了简明、直观的理论依据.

4 能量耗散率的各向异性

现在进一步考察在一般方向上的耗散行为,并对前述各向异性规律做更为普遍的推广.

对于给定的约化速度u,可以用 ∂ P/∂θv来表征能量耗散随运动方向的变化.当耗散率与临界速度对方向角的依赖具有一致的变化趋势时,∂P/∂θv与∂ vc/∂θv符号相同.这意味着,若将 P 写为 vc和u的函数,必然有

下面就从三个方面来论证在偶极BEC中(7)式关系总是成立的.

首先,考虑 u ≫1 时的高耗散情形.此时耗散率的主要贡献来自于大动量的准粒子激发,能谱的各向异性不再重要.将(1)式中的 Eq作大q展开并保留至头阶项,可以推出耗散率在高速极限下的渐近形式为 P=Av3.把该结果重新用 vc和u表示,于是得到

图1 沿不同方向的能量耗散率 P在(a)高耗散区间和(b)低耗散区间随约化速度u的变化曲线.离散点为数值计算的结果,(a)和(b)中的虚线分别对应于(8)式和(9)式给出的渐近表达式.从上至下,各曲线的临界速度分别为vc=1.41c0,1.22c0,1.07c0,0.89c0,0.71c0.两图采用的图例相同.偶极相互作用参数取为ϵdd=0.5Fig.1.Energy dissipation rate P as a function of reduced velocity u along different directions in(a)high dissipation regime and(b)low dissipation regime.Discrete symbols are numerical results,and dashed lines in(a)and(b)correspond to the asymptotic expressions(8)and(9),respectively.From up to bottom,the critical velocities for each line are given by vc=1.41c0,1.21c0,1.07c0,0.89c0,and 0.71c0.The two plots use the same legend.The relative strength of dipolar interaction is set as ϵdd=0.5 .

显然,不等式关系(7)式在此情形下是成立的.

其次,再来考虑u略大于1时的低耗散情形.此时耗散率的主要贡献来自于小动量的准粒子,且激发动量的极角集中在θc附近,θc对应于v=vc时产生声子激发的传播方向.将(1)式中的被积函数在q=0和θq=θc处作展开,经过一系列较为繁琐的推导(具体步骤见附录A),最终得出

对于给定的 u,(9)式中的P是vc的增函数,因此在临界速度附近,不等式(7)式也是成立的.

最后,对于中等u值、中等耗散的情形,(1)式中的积分无法得到简单的解析结果,需要借助数值计算验证耗散率与临界速度的相关性.图1给出了杂质沿不同方向运动时耗散率P随约化速度u的变化曲线.可以看出,在 vc更大的方向上耗散率也更高,这表明不等式(7)式对于任意 u ＞1 的速度区间都是成立的.在 u ≫1和u →1+两种极限情形下,数值结果与前面推导出的渐近表达式符合得很好.

根据以上分析,我们得出结论: 在偶极BEC中,耗散率与临界速度随方向角的变化趋势总是一致的.这个结论将(6)式描述的各向异性关系推广为更加普遍的形式,也为实验未观测的非主轴方向上的耗散现象提供了新的理论预言.

5 结 论

在前文的计算中,假定了每个元激发过程只会产生一个准粒子,而没有考虑同时产生多个准粒子的元激发过程.可以证明,激发单个准粒子的临界速度一定小于等于激发多个准粒子的临界速度(具体证明见附录B).因此,当杂质速度由小及大逐渐增加时,单粒子激发必然先于多粒子激发出现.另一方面,在Bogoliubov近似下,密度扰动引起的元激发过程自然地分成了两类: 单粒子激发对应于原子从零动量态到非零动量态的跃迁,多粒子激发对应于原子在非零动量态之间的跃迁.由于占据零动量态的凝聚原子远远多于非凝聚原子,所以多粒子激发对于耗散的贡献是可以忽略不计的.

综上所述,本文运用线性响应理论,全面考察了偶极BEC中点状杂质沿不同方向运动的能量耗散率.在两个主轴方向上,解析地得到了耗散率随约化速度变化的函数形式,证明了平行极化轴方向的耗散率比垂直极化轴方向的更高,从而定性地解释了实验上观测到的各向异性耗散现象.此外,我们还普遍地论证了,耗散率与临界速度随方向角的变化趋势总是一致的,这为进一步的实验研究提供了新的理论预言.本文的结论对于探索其他冷原子体系(如自旋-轨道耦合BEC[19−21])的各向异性超流现象也有一定的借鉴意义.

附录A

考虑杂质沿偏离主轴的方向运动.根据体系的轴对称性,不妨取速度v在x-y平面内投影的方位角φv=0.将q·v=qv(sinθqcosφqsinθv+cosθqcosθv)代入(1)式可算出φq的积分,再对q和θq做适当的变量替换后,得到

对于一般情形,(A1)式中的积分并不能完整地解析算出.下面重点讨论u略大于1时耗散率的渐近形式.

当杂质速度刚好达到vc时,触发耗散的准粒子为q→0的长波声子,其传播方向的极角θc满足关系由于激发条件 Eq=q ·v 的等式两端随v和q连续变化,故当杂质速度比 vc略大时,贡献耗散的准粒子只出现在小动量区间,且动量的极角集中在θc附近.据此分析,可以将函数展开为

其中 γc=cosθc,展开系数 λ0,λ1,λ2,λ3由泰勒公式给出

λ1表 达式中的因子sgn(γv)意为取γv的正负号.

易于验证,当 u →1 时,λ2和λ3均 小于零,这表明在vc附近函数可简化为一个抛物面.在此近似下,耗散率可由(A1)式解析地算出:

将(A4)式—(A7)式中的展开系数代入(A8)式,并将所得结果保留至 u -1 的平方阶,即得到(9)式给出的渐近形式.

附录B

触发耗散的元激发过程可以产生一个或多个准粒子.元激发的能量、动量以及杂质速度之间满足关系式Eq+q·v=0 .对于多粒子激发过程,Eq和q 分别表示全体准粒子的总能量和总动量.根据这个关系,要同时激发j个准粒子,杂质速度必须达到阈值

将该结果与(B1)式对比,即可推断,产生多粒子激发的阈值速度一定大于产生单粒子激发的阈值速度,

这意味着触发耗散的临界速度总是由单粒子激发过程决定的,即

上述结论作为朗道超流判据的推广,同时适用于各向同性体系与各向异性体系.对于偶极BEC,将Bogoliubov近似下的元激发能谱代入(B5)式,即得到(3)式给出的临界速度[13].