PT对称势下自旋轨道耦合玻色-爱因斯坦凝聚体的能带结构.pdf

西 北 师 范 大 学 学 报(自然科学版)第6 0卷2 0 2 4年第3期 J o u r n a l o f N o r t h w e s t N o r m a l U n i v e r s i t y(N a t u r a l S c i e n c e)V o l.6 0 2 0 2 4 N o.3 D O I:1 0.1 6 7 8 3/j.c n k i.n w n u z.2 0 2 4.0 3.0 0 5收稿日期:2 0 2 3 1 2 1 0;修改稿收到日期:2 0 2 4 0 2 0 5基金项目:国家自然科学基金资助项目(1 2 0 6 5 0 2 2)作者简介:邵凯花(1 9 9 7),女,甘肃天水人,硕士研究生.主要研究方向为非线性物理.E m a i l:2 8 5 9 8 5 4 4 6 0q q.c o m*通信联系人,男,教授,博士,博士研究生导师.主要研究方向为非线性物理.E m a i l:s h i y r n w n u.e d u.c nP T对称势下自旋轨道耦合玻色-爱因斯坦凝聚体的能带结构邵凯花,马金萍,王青青,席保龙,石玉仁*(西北师范大学 物理与电子工程学院,甘肃 兰州 7 3 0 0 7 0)摘要:数值研究P T对称势下具有自旋轨道耦合玻色-爱因斯坦凝聚体的能带结构.用傅立叶配置法求解布洛赫能带,分析P T对称晶格势的深度、虚部相对大小、自旋轨道耦合强度和拉比耦合强度对能带结构的影响.结果表明,第一、二能带的上下边沿随晶格深度的增加而上升,但能带宽度变窄;当P T对称晶格势的虚部大于临界值时,第一和第二能带发生重叠;随着自旋轨道耦合强度的增大,第一和第二能带的上下边界均减小;拉比耦合强度对能带结构也具有重要影响.关键词:玻色-爱因斯坦凝聚;P T对称势;自旋轨道耦合;能带结构中图分类号:O 1 4 5 文献标志码:A 文章编号:1 0 0 1-9 8 8(2 0 2 4)0 3-0 0 4 2-0 6B a n d s t r u c t u r e s i n s p i n-o r b i t c o u p l e d B o s e-E i n s t e i nc o n d e n s a t e w i t h P T s y mm e t r i c p o t e n t i a lS HAO K a i-h u a,MA J i n-p i n g,WANG Q i n g-q i n g,X I B a o-L o n g,S H I Y u-r e n(C o l l e g e o f P h y s i c s a n d E l e c t r o n i c E n g i n e e r i n g,N o r t h w e s t N o r m a l U n i v e r s i t y,L a n z h o u 7 3 0 0 7 0,G a n s u,C h i n a)A b s t r a c t:T h e b a n d s t r u c t u r e s i n s p i n-o r b i t c o u p l e d B o s e-E i n s t e i n c o n d e n s a t e w i t h P T s y mm e t r i c p o t e n t i a l a r e i n v e s t i g a t e d n u m e r i c a l l y.T h e B l o c h b a n d s a r e o b t a i n e d b y t h e F o u r i e r c o l l o c a t i o n m e t h o d a n d t h e e f f e c t s o f t h e d e p t h o f P T s y mm e t r i c l a t t i c e p o t e n t i a l,t h e m a g n i t u d e o f i m a g i n a r y p a r t,t h e s p i n-o r b i t c o u p l e d s t r e n g t h a n d t h e R a b i c o u p l i n g s t r e n g t h o n b a n d s t r u c t u r e s a r e a n a l y z e d.T h e r e s u l t s i n d i c a t e t h a t t h e u p p e r a n d l o w e r b o u n d a r i e s o f t h e f i r s t a n d s e c o n d b a n d s i n c r e a s e w i t h t h e i n c r e a s e o f l a t t i c e d e p t h,b u t t h e w i d t h o f b a n d b e c o m e s n a r r o w e r.T h e f i r s t a n d s e c o n d b a n d s o v e r l a p w h e n t h e i m a g i n a r y p a r t o f t h e P T s y mm e t r i c l a t t i c e p o t e n t i a l i s l a r g e r t h a n t h e c r i t i c a l v a l u e.F o r t h e f i r s t a n d s e c o n d b a n d s,t h e i r u p p e r a n d l o w e r b o u n d a r i e s d e c r e a s e w i t h t h e i n c r e a s e o f s p i n-o r b i t c o u p l e d s t r e n g t h.I t i s a l s o f o u n d t h a t t h e R a b b i c o u p l i n g s t r e n g t h a l s o h a s a n i m p o r t a n t e f f e c t o n t h e b a n d s t r u c t u r e.K e y w o r d s:B o s e-E i n s t e i n c o n d e n s a t e;P T s y mm e t r i c p o t e n t i a l;s p i n-o r b i t c o u p l i n g;b a n d s t r u c t u r e24 2 0 2 4年第3期 邵凯花等:P T对称势下自旋轨道耦合玻色-爱因斯坦凝聚体的能带结构 2 0 2 4 N o.3B a n d s t r u c t u r e s i n s p i n-o r b i t c o u p l e d B o s e-E i n s t e i n c o n d e n s a t e w i t h P T s y mm e t r i c p o t e n t i a l 玻 色-爱 因 斯 坦 凝 聚(B o s e-E i n s t e i n c o n d e n-s a t e,B E C)是指系统的温度低于临界温度时,所有原子 都 聚 集 在 最 低 能 级 上 的 宏 观 量 子 现 象1.1 9 9 5年,C o r n e l l等2通过激光冷却和射频蒸发冷却技术首次在1 7 0 n K的低温下获得了8 7R b原子的B E C.同年,K e t t e r l e等3利用塞曼减速技术冷却的原子束系统实现了2 3N a原子的B E C.之后,科学家们陆续实现了各种原子的B E C,使得其在物理研究中获得了重要的地位.在超冷原子气体中,自旋轨道耦合(S p i n-o r b i t c o u p l i n g,S O C)的加入为探索新的量子现象提供了平台.自2 0 1 1年实验上成功实现了带有S O C的8 7R b B E C4以来,激发了人们对量子相变、拓扑激发以及自旋电子器件等的兴趣.实验上B E C被外势束缚,如光晶格势5、谐振子势6、同心耦合环形阱7、径向周期阱8等.在光晶格周期势中,非线性相互作用会明显改变布洛赫带的结构,使能带在第一布里渊区的边界处呈现出燕尾形状的环状结构8.并且,晶格势和S O C的结合极大地扩展了冷原子的可调性,为研究B E C的相关物理特性 提 供 了 很 好 的 平 台,例 如B l o c h 振 荡9、L a n d a u-Z e n e r 隧穿1 0、局域化等.宇称-时间(p a r i t y-t i m e,P T)对称是物理学的热门研究领域.1 9 9 8年,B e n d e r等1 1提出了一个具有P T对称的非厄米哈密顿量仍然具有完全实数的能谱.同时,非厄米系统中的P T对称势平衡了系统的增益和损耗.实验中,增益可通过使用原子激光器将原子注入冷凝物来实现1 2,损耗采用激光束将原子激发到激发态,然后通过光子反冲将它们从凝聚物中喷射出来而实现1 3.随着研究的进展,许多研究表明P T对称性在各个领域具有丰富的研究内容.2 0 2 2年,Q i n等1 4的研究表明,在P T对称伪自旋1/2系统中,与动量相关的平衡增益和损耗可以合成虚S O C,可以改变系统的能谱.Z h a n g等1 5通过引入虚势,从理论上研究了一维非厄米S O C S u-S c h r i e f f e r-H e e g e r模 型 的 能 带 结构.本文考虑一准二维S O C-B E C系统,用傅立叶配置法求解在P T对称晶格势下的线性能带结构,详细研究系统中各物理参量(如晶格的深度、自旋轨道耦合强度、拉比耦合强度等)对能带结构的影响.1 理论模型与方法1.1 理论模型将一准二维S O C-B E C束缚在外势中,经无量纲处理后,在平均场近似下系统的动力学行为可用如下耦合G r o s s-P i t a e v s k i i(G P)方程组描述1 6i t=-2+V+g1 12+g1 22+i -+,i t=-2+V+g2 12+g2 22+i+,(1)其中,=(r,t)表示波函数;r=(x,y),g1 1,g2 2和g1 2分别表征原子种内和种间相互作用强度,其具体表达式请参看文献1 6;2=2 x2+2 y2是二维拉普 拉斯算符;和分 别 表 示S O C强度和拉比耦合强度,=/xi/y.V为P T对称晶格势,可写为1 7V(x,y)=VR(x,y)+iVI(x,y),(2)满足V(x,y)=V*(-x,-y),即V(x,y)的实部VR为偶函数,虚部VI为奇函数.这里我们选择P T对称晶格势的实部和虚部分别为1 8VR(x,y)=V0(c o s2x+c o s2y),(3)VI(x,y)=V0W0(s i n(2x)+s i n(2y),(4)其中,V0表示晶格深度,W0为P T对称晶格势虚部的相对大小.图1为该周期晶格势的实部和虚部的等值线图,其在空间呈现明显的周期结构.1.2 线性能带结构若寻找方程(1)形式为(x,y,t)=(x,y)e-it的定态解,则化学势和(x,y)满足=(-2+V+g1 12+g1 22)+i-+,=(-2+V+g1 22+g2 22)+i+.(5)由于(x,y)和V(x,y)均为周期函数,其有界解可用布洛赫定理(x,y)=ei kr(x,y)表示,其中波矢k=(kx,ky).若(x,y)为小量,则方程(5)中非线性项可以忽略,得到本征方程34西 北 师 范 大 学 学 报(自然科学版)第6 0卷 J o u r n a l o f N o r t h w e s t N o r m a l U n i v e r s i t y(N a t u r a l S c i e n c e)V o l.6 0(k)=-(+i k)2+V+i-+,(k)=-(+i k)2+V+i+.(6)特征值(k)的集合构成布洛赫带(B l o c h b a n d s).图1 二维周期晶格势F i g 1 C o n t o u r p l o t s o f t h e t w o-d i m e n s i o n a l l a t t i c e p o t e n t i a l1.3 傅立叶配置法一般情况下,很难得到本征方程(6)的精确解析解,可考虑对其进行数值求解,如有限差分法或傅立叶配置法等.傅立叶配置法是一种高精度的数值方法,在此做一简要介绍.首先,将无界域(x,y)R2截断为有界域:-Lx/2xLx/2,-Ly/2yLy/2,其中Lx和Ly分别为x和y方向的区间宽度.然后,在此有界域上定义内积f,g=f(x,y)g*(x,y)dxdy,这里“*”表示复共轭.设(x,y)=+m=-+n=-a,m nei(n k0 x+m k0yy),V(x,y)=+m=-+n=-bm nei(n k0 x+m k0yy),其中k0 x=2Lx,k0y=2Ly,bm n=1LxLyLx/2-Lx/2Ly/2-Ly/2V(x,y)ei(n k0 x+m k0yy)dxdy.代入本征方程(6),并对求和指标进行截断后,得am n=(-(in k0 x)2-(im k0y)2+2n kxk0 x+2m kyk0y+(k2x+k2y)am n+Nyl=-Ny+Nxj=-Nxbm-l,n-jai j+(c1-n k0 x+i mk0y)am n,am n=(-(in k0 x)2-(im k0y)2+2n kxk0 x+2m kyk0y+(k2x+k2y)am n+Nyl=-Ny+Nxj=-Nxbm-l,n-jai j+(c2-n k0 x+i mk0y)am n.n=-Nx,-1,0,1,Nx;m=-Ny,-1,0,1,Ny.其中c1=-kx+i ky,c2=-kx-i ky.此时本征方程(6)便转化为离散的矩阵形式AX=X,其中本征矢X=am nam n ,矩阵A是(2Nx+1)(2Ny+1)阶稠密方阵.随着Nx和Ny的增加,会导致存储量和计算量急剧增加.在后面计算中,取Ny=Nx=2 0.作为特例,图2给出了V0=-4,W0=0.1,=0.5时的布洛赫能带结构.从图3中可以看到一系列带隙将布洛赫能带分隔开来.通常将位于最低布洛赫能带下面的带隙称为半无限带隙,其它带隙依次称为第一带隙、第二带隙等等.这些带隙在系统中起着非常重要的作用.图3中,半无限带隙的 区 域 约 为-5.4 2,第 一 带 隙 区 域 约 为0.4 8 4-4.4 7.图3 a和3 b分别是图2 a沿kx和ky方向的投影图.从图3可以看出,布洛赫能带沿kx方向关于kx=0不对称,而沿ky方向关于ky=0对称.44 2 0 2 4年第3期 邵凯花等:P T对称势下自旋轨道耦合玻色-爱因斯坦凝聚体的能带结构 2 0 2 4 N o.3B a n d s t r u c t u r e s i n s p i n-o r b i t c o u p l e d B o s e-E i n s t e i n c o n d e n s a t e w i t h P T s y mm e t r i c p o t e n t i a l图2 P T对称光晶格势中B E C的线性布洛赫能带结构(V0=-4,W0=0.1,=0.5)F i g 2 L i n e a r B l o c h b a n d s t r u c t u r e s o f B E Ci n P T s y mm e t r i c l a t t i c e p o t e n t i a l(V0=-4,W0=0.1,=0.5)2 数值结果及分析下面我们系统地计算晶格深度V0,P T对称晶格势的虚部W0,S O C强度以及拉比耦合强度的变化对能带结构的影响.2.1 P T对称晶格势对能带结构的影响数值结果表明,P T对称晶格势的深度V0与W0对能带结构有着明显的影响.图4给出了能带结构随 V0的变化,所选参数为W0=0.1,=0.5.从图4可以看出,当晶格深度V0很小时,能带结构出现了交叠,此时只存在半无限带隙,第一带隙和第二带隙均没有打开.随着晶格深度V0的增加,第一能带的下边沿和上边沿均上升,但能带宽度变窄.当V0增大至2.3 7时,第一带隙开始出现且其宽度随V0的增大而增大.当V0大于2.4 7时,第二带隙开始出现,且逐渐变宽.图5给出了V0=-4,=0.5时能带结构随W0的变化.从图5可以看到,随着W0的增大,半无限带隙的上边沿逐渐增大.当W00.4时,图3 线性布洛赫能带结构的轴向侧视图,阴影区域对应于线性能带F i g 3 A x i a l s i d e v i e w o f l i n e a r B l o c h b a n ds t r u c t u r e.T h e s h a d e d a r e a s c o r r e s p o n d t ol i n e a r e n e r g y b a n d s图4 能带结构随晶格深度V0的变化F i g 4 T h e v a r i a t i o n o f b a n d s t r u c t u r e w i t h l a t t i c e d e p t h时,第一能带和第二能带的宽度变化并不明显,但第一带隙和第二带隙的宽度则随W0的增大而明显变窄.当W0增大至临界值0.4时,所有能带出现交叠,此时只存在半无限带隙.2.2 自旋轨道耦合强度对能带结构的影响自旋轨道耦合强度也是影响能带结构的重要参数之一.图6给出了能带结构随的变化,所选54西 北 师 范 大 学 学 报(自然科学版)第6 0卷 J o u r n a l o f N o r t h w e s t N o r m a l U n i v e r s i t y(N a t u r a l S c i e n c e)V o l.6 0 参数为V0=-4,W0=0.1,=0.5.从图6可以看出,随着的增大,第一能带和第二能带的上边沿均逐渐减小,能带整体向下移动,其宽度也逐渐变窄.但当增加到一定值时,它们的宽度变化不再明显.在此过程中,第一带隙和第二带隙始终存在.图5 能带结构随P T对称晶格势虚部W0的变化F i g 5 T h e v a r i a t i o n o f b a n d s t r u c t u r e w i t h t h e i m a g i n a r yp a r t o f P T s y mm e t r i c l a t t i c e p o t e n t i a l图6 能带结构随S O C强度的变化F i g 6 T h e v a r i a t i o n o f b a n d s t r u c t u r e s w i t h t h eS O C s t r e n g t h2.3 拉比耦合对能带结构的影响图7给出了V0=-4,W0=0.1,=0.5时拉比耦合强度对能带结构的影响.从图7可以看出,在所选参数条件下,第一能带的上下边沿均随的增大而减小,且两边沿大致呈平行状.值得注意的是,当0.3 2时,系统产生图7所示第一带隙,且其宽度随着的增大而增大.但当足够大时,第一带隙的宽度变化不再明显.图7中所示第二带隙的变化则与此有点相反.随着的增加,第二带隙的宽度明显变窄.当0.9 6时,第二带隙消失.该结果表明拉比耦合强度对于系统的带隙结构也有着重要的影响.在实验中,可通过调节其值控制和改变系统的能带结构.图7 能带结构随拉比耦合强度的变化F i g 7 T h e v a r i a t i o n o f b a n d s t r u c t u r e s w i t h t h eR a b b i c o u p l i n g s t r e n g t h3 结论数值研究了P T对称晶格势 下准二维S O C-B E C中的线性能带结构.在平均场近似下,B E C的动力学行为可用耦合G P方程组来描述.用傅立叶配置法数值求解了不同参数下的能带结构.结果表明,随着晶格深度V0的增加,第一能带的上下边沿均上升,但能带宽度变窄.P T对称势的虚部W0存在临界值,当W0小于临界值时,第一能带和第二能带的宽度变化并不明显;W0大于临界值时所有能带出现交叠.随着S O C强度的增大,第一能带和第二能带的上边沿均逐渐减小,能带整体向下移动,其宽度也逐渐变窄.系统的能带和带隙结构也受拉比耦合强度的影响.故实验中可通过调节这些参数来改变能带结构.参考文献:1 漆伟,李明,郭晓刚.三体作用下玻色-爱因斯坦凝聚体中的P T-对称孤子J.西北师范大学学报(自然科学版),2 0 2 2,5 8(6):6 6.2 AN D E R S ON M H,E N S HE R J R,MA T THEWS M R,e t a l.O b s e r v a t i o n o f B o s e-E i n s t e i n 64 2 0 2 4年第3期 邵凯花等:P T对称势下自旋轨道耦合玻色-爱因斯坦凝聚体的能带结构 2 0 2 4 N o.3B a n d s t r u c t u r e s i n s p i n-o r b i t c o u p l e d B o s e-E i n s t e i n c o n d e n s a t e w i t h P T s y mm e t r i c p o t e n t i a lc o n d e n s a t i o n i n a d i l u t e a t o m i c v a p o rJ.P h y s R e p,1 9 9 5,2 6 9(5 2 2 1):1 9 8.3 梁青青,周小燕,李好财.有限温度下含三体相一经作用的玻色-爱因斯坦凝聚体的集体激发特性J.西北师 范 大 学 学 报(自 然 科 学 版),2 0 2 3,5 9(4):5 4.4 L I N Y J,G A R C I A K J,S P I E LMAN I B.S p i n-o r b i t-c o u p l e d B o s e-E i n s t e i n c o n d e n s a t e sJ.N a t u r e,2 0 1 1,4 7 1:8 3.5 KUHN S,J U D D T E.T r a n s p o r t o f d i p o l a r B o s e-E i n s t e i n c o n d e n s a t e s i n a o n e-d i m e n s i o n a l o p t i c a l l a t t i c eJ.P h y s R e v A,2 0 1 3,8 7(2):0 2 3 6 0 8.6 XU X Q,HAN J H.S p i n-o r b i t c o u p l e d B o s e-E i n s t e i n c o n d e n s a t e u n d e r r o t a t i o nJ.P h y s R e v L e t t,2 0 1 1,1 0 7(2 0):2 0 0 4 0 1.7 MA L E T F,KAVOU L AK I S G M,R E I MANN S M.M i x t u r e s o f B o s e g a s e s c o n f i n e d i n c o n c e n t r i c a l l y c o u p l e d a n n u l a r t r a p sJ.P h y s R e v A,2 0 1 2,8 1(1):0 1 3 6 3 0.8 WAN G Y J,WE N L,GUO H,e t a l.S p i n-o r b i t-c o u p l e d B o s e-E i n s t e i n c o n d e n s a t e s i n r a d i a l l y p e r i o d i c p o t e n t i a l sJ.J P h y S o c J p n,2 0 1 9,8 8(2):0 2 4 0 0 5.9 MA CHHO LM M,P E TH I C K C J,S M I TH H.B a n d s t r u c t u r e,e l e m e n t a r y e x c i t a t i o n s,a n d s t a b i l i t y o f a B o s e-E i n s t e i n c o n d e n s a t e i n a p e r i o d i c p o t e n t i a lJ.P h y s R e v A,2 0 0 3,6 7(5):0 5 3 6 1 3.1 0 J I W,Z HAN G K,Z HAN G W,e t a l.B l o c h o s c i l l a t i o n s o f s p i n-o r b i t-c o u p l e d c o l d a t o m s i n a n o p t i c a l l a t t i c e a n d s p i n-c u r r e n t g e n e r a t i o nJ.P h y s R e v A,2 0 1 9,9 9(2):0 2 3 6 0 4.1 1 KONO TO P V V,K E V R E K I D I S P G,S A L E R NO M.L a n d a u-Z e n e r t u n n e l i n g o f B o s e-E i n s t e i n c o n d e n s a t e s i n a n o p t i c a l l a t t i c eJ.P h y s R e v A,2 0 0 5,7 2(2):0 2 3 6 1 1.1 2 B E N D E R C M,B O E T T CHE R S.R e a l s p e c t r a i n n o n-H e r m i t i a n H a m i l t o n i a n s h a v i n g P T s y mm e t r yJ.P h y s R e v L e t t,1 9 9 8,8 0(2 4):5 2 4 3.1 3 R O B I N S N P,F I A L C,J E P P E S E N M,e t a l.A p u m p e d a t o m l a s e rJ.N a t P h y s,2 0 0 8,4(9):7 3 1.1 4 L I J,HA R T E R A K,L I U J,e t a l.O b s e r v a t i o n o f p a r i t y-t i m e s y mm e t r y b r e a k i n g t r a n s i t i o n s i n a d i s s i p a t i v e F l o q u e t s y s t e m o f u l t r a c o l d a t o m sJ.N a t C o mm,2 0 1 9,1 0(1):8 5 5.1 5 Q I N J L,Z HOU L,D ONG G J.I m a g i n a r y s p i n-o r b i t a l c o u p l i n g i n p a r i t y-t i m e s y mm e t r i c s y s t e m s w i t h m o m e n t u m-d e p e n d e n t g a i n a n d l o s sJ.N e w J P h y s,2 0 2 2,2 4(6):0 6 3 0 2 5.1 6 WAN G Q Q,TU P,MA J P,e t a l.V e c t o r g a p s o l i t o n s o f s p i n-o r b i t-c o u p l e d B o s e-E i n s t e i n c o n d e n s a t e i n s q u a r e o p t i c a l l a t t i c eJ.J P h y s B,2 0 2 4,5 7(6):0 6 5 2 0 1.1 7 ME NG H J,WANG J,S H I Y R,e t a l.V e c t o r g a p s o l i t o n s o f s p i n-o r b i t-c o u p l e d B o s e-E i n s t e i n c o n d e n s a t e i n h o n e y c o m b o p t i c a l l a t t i c e sJ.P h y s R e v E,2 0 2 3,1 0 8(3):0 3 4 2 1 5.1 8 S CHWA R Z L,C A R T A R I U S H,MU S S L I MAN I Z H,e t a l.V o r t i c e s i n B o s e-E i n s t e i n c o n d e n s a t e s w i t h P T-s y mm e t r i c g a i n a n d l o s sJ.P h y s R e v A,2 0 1 7,9 5(5):0 5 3 6 1 3.1 9 X I E J N,S U Z K,CHE N W C,e t a l.D e f e c t s o l i t o n s i n t w o-d i m e n s i o n a l p h o t o n i c l a t t i c e s w i t h p a r i t y-t i m e s y mm e t r yJ.O p t C o mm,2 0 1 4,3 1 3(0 0 3 0-4 0 1 8):1 3 9.(责任编辑 孙对兄)74