强关联电子体系二维相干光谱的理论研究评述*

李子龙 万源3)

1) (中国科学院物理研究所,凝聚态理论与材料计算重点实验室,北京 100190)

2) (中国科学院大学物理科学学院,北京 100049)

3) (松山湖材料实验室,东莞 523808)

非线性相干光谱是一种测量材料非线性光学响应的谱学手段.相比线性光谱,非线性相干光谱具有多个时间变量,能够提供材料的更多信息.作为非线性相干光谱的代表之一,二维相干光谱在众多领域取得了丰硕的成果.在化学、生物学等领域,二维相干光谱已经展现出相比线性谱学的优越性.随着太赫兹技术的发展,二维相干光谱在强关联体系也显现出巨大潜力,相关的理论与实验工作正在开展.本文简要概括了强关联体系二维相干光谱的发展现状,介绍了二维相干光谱的基本概念与理论工具,分析了二维相干光谱的主要特点,并重点总结了我们研究组近几年关于强关联体系二维相干光谱的理论研究工作.

1 研究现状

凝聚态物理的进步与谱学手段的发展息息相关.例如,中子散射、扫描隧道光谱以及角分辨电子能谱的发展拓宽了我们对材料性质的认识.在强关联电子学领域,二维相干光谱作为一种新兴的谱学手段受到越来越多的重视.二维相干光谱利用多束相干电磁波来探测材料的非线性响应,因而在频域上产生一个二维的光谱来表征材料的非线性响应[1,2,3].

相比传统的线性光谱,二维相干光谱不仅能得到材料中光学激发的信息,也能揭示激发之间的联系.在红外以及高频波段,二维相干光谱已经被广泛应用于化学、生物学等领域,用来刻画化学反应以及生命活动的动力学过程[1-4].太赫兹技术的发展使这一谱学手段达到了强关联电子系统的能量窗口.太赫兹波段的频率范围为 0.1—10×1012Hz,对应的能量范围大致为 1meV—0.1 eV.高温超导体的赝能隙和超导转变温度对应的能量为几十个mV[5],几种量子磁性材料的磁激发对应的能量为几个mV[6-8],典型拓扑绝缘体的能隙也为几十个mV[9].这些强关联体系的特征能量尺度在太赫兹的波段范围,因而可以用太赫兹二维相干光谱进行研究.目前,这一谱学手段被应用在量子阱[10]、反铁磁磁体[11]、电子玻璃[12]等强关联体系中.

与不断涌现的实验成果相比,相关的理论仍然不成熟.首先,人们并不清楚强关联体系的非线性响应是否与体系的性质存在直截了当的关系;其次,目前大多数理论工作仍然使用能级的图像来描述非线性响应过程.这样一套图像将非线性响应与量子系统在时间维度上的干涉现象联系起来,为二维相干光谱提供了一个直观的物理解释.然而,量子多体系统具有空间上的广延性.因此,多体系统的非线性响应,除了时间维度上存在干涉效应外,还可能存在空间维度上的干涉效应.因此,针对强关联体系的非线性光谱研究需要一套新的物理图像.

理论工作者们利用数值和解析手段计算了一些重要强关联模型的非线性响应,并通过联系非线性光谱学中已有的分析手段,对这些模型的非线性响应光谱做出解释.与此同时,重整化、玻色化以及规范理论等多体物理的分析工具被引入并与原有的双边费曼图等方法结合,从而给出强关联体系二维相干光谱更深刻的理解.

Wan 与Armitage[13]从理论上研究了一种典型的强关联自旋系统—一维横场伊辛(Ising)模型.自旋算符作用到基态激发两个分数激发—“自旋子”.线性光谱是不同动量“自旋子”的连续谱,体系的退相干时间等信息隐藏在谱内.在不引入耗散的前提下,退相干时间是无穷大,因而线性谱中的展宽与耗散无关.二维相干光谱通过光子回波将引起展宽的失相过程消除,因而重现了“自旋子”的尖峰结构.在引入耗散过程后,退相干时间有限,二维光谱的尖峰结构出现展宽.这说明二维光谱将耗散与非耗散过程引起的展宽区分开来,从而获得相比线性光谱关于激发更多的信息.虽然横场伊辛模型是一维模型,但由于分数激发是自旋液体中普遍存在的激发,因而他们认为二维相干光谱可以用来分析更高维的、具有连续谱结构的自旋液体模型.

随后,Choi 等[14]进行了Kitaev 自旋液体的二维相干光谱的理论研究.利用Kitaev 模型的严格可解性,他们计算了体系的三阶非线性响应,发现一个自旋算符可以从体系的基态中激发一个马约拉纳费米子与两个量子涡旋(vison),但这两种激发在二维光谱中的特征不同,因而可以被区分出来.其中,二维光谱的对角线方向尖锐的峰代表了马约拉纳费米子的激发,这对实验上认证马约拉纳费米子有指导意义.

作为工作的延续,Nandkishore 等[15]对具有能隙的自旋液体进行了进一步研究.他们考虑了几种典型的有能隙的自旋液体模型,如伊辛模型、二维toric code 模型、Haah’s code 模型以及X-cube 模型.对于常规自旋液体,各类激发在线性光谱中无法表现出显著的特征.但他们通过理论计算发现,这些激发在二维相干光谱上都有明显特征.对于非常规自旋液体,第I 类分形相(X-cube)的激发在线性光谱以及二维相干光谱中都表现出明显尖峰特征.与之相比,第II 类分形相(Haah’s code)需要更高阶的非线性响应来获取体系激发的特征,这也对实验提出了更大的挑战.

Parameswaran和Gopalakrishnan[16]从理论上研究了无序量子磁体的非线性响应.以一维无序量子横场伊辛模型为例,他们通过实空间重整化技术,说明在相变点附近该模型可以用一系列孤立二能级系统描述.紧接着,他们人为地引入声子与体系中二能级系统的相互作用,从而使体系的激发存在寿命.在这种情况下,他们考虑了对应于“泵浦-探测(PP)”技术的非线性响应函数.通过分析他们发现,相比线性光谱,二维相干光谱不仅能够将弛豫过程与退相干过程区分开,也能够得到体系的态密度和重整化后的磁矩等物理信息.

Phuc和Trung[17]将二维相干光谱运用到一个典型的电子体系—双带哈伯德模型中.通过模拟计算,他们发现在没有占位能的情况下,二维光谱的峰都分布在对角线上.他们将这种特征归因于无相互作用的准激发.随着占位能增加,非对角峰开始显现,体系中相互作用的大小正比于非对角峰的强度;在体系进入莫特相后,二维相干光谱表现为一个原点位置单一的峰,这对应于基态到激发态的跃迁.此外,他们进一步发现,二维相干光谱也可以用来分析强关联体系中的含时演化问题.

可以发现,以上几个工作都局限于有能隙的体系.这些体系在经过一定处理后,都可以直接利用分析原子、分子体系非线性相干光谱的方法,因而无法说明二维相干光谱在强关联体系中的新特性.为了进一步说明二维相干光谱在强关联体系的巨大潜力,Li 等[18]讨论了一种典型的无能隙系统—拉廷格(Luttinger)液体.他们发现拉廷格液体的二维相干光谱中仍然存在光子回波信号.与其他光子回波产生的机理不同,他们将这一现象归因于一种特殊的时空间干涉过程.第1 个自旋算符产生两个分数激发,第2和第3 个自旋算符分别改变两个分数激发的手性,最终二者汇聚于原点并被湮灭掉.在这个过程中,两个分数激发的世界线犹如被一面透镜聚焦于焦点,因而这一过程又被命名为“透镜”效应.当体系存在色散或者耗散,这种“透镜”效应不再完美,因而光子回波信号随时间衰减.这一工作将非线性光谱的时间干涉图像扩展为时空间干涉图像,对于理解强关联体系的非线性光谱有启发性.

在接下来的部分中,本文将介绍如下内容:第2节以二能级系统为例,简要介绍二维相干光谱的含义与方法;第3节和第4节以本研究组发表的两篇工作为例,说明强关联电子体系二维相干光谱的新特点.其中,第3节主要介绍二维相干光谱在有能隙体系—横场伊辛模型的应用[13],第4节主要介绍二维相干光谱在无能隙体系—拉廷格液体的应用[18].

2 二维相干光谱基本概念—以二能级系统为例[1,2,12]

2.1 响应理论

考虑一个由二能级模型描述的分子体系:

其中,σz是泡利矩阵,Δ是二能级系统的激发能.希尔伯特空间的任何一个态都可以由|0〉和|1〉展开.假设交变电场的作用是使得二能级系统跃迁,哈密顿量可以加入一项含时相互作用:

其中,E(t)是电磁场的电场分量,µ对应于分子体系中的偶极子算符,µ=µ10σx.假设体系初始处于基态,其含时演化的薛定谔方程为

假设含时相互作用W(t) 很弱,电场表现为δ函数的形式:E(t)=ϵδ(t).将ϕ(t)=c0(t)|0〉+c1(t)|1〉代入(3)式O(ϵ),得到:

在与电磁场相互作用后,密度矩阵的形式为

密度矩阵此后自由演化,也就是

从而得到任意时刻体系中电偶极矩的值:

由麦克斯韦方程知,周期运动的电偶极矩产生其运动频率的电磁场,因而(7)式是一可测量.响应的强度与电偶极子的跃迁矩阵元和电场强度的乘积成比例,Hamm 等[19]仔细讨论了不同体系电偶极矩的估计方法以及线性和非线性光谱实验需要的光场强度.以上得出了线性响应理论的主要结果,但与此同时,物理图像在计算中并不清晰,因而需要使用一套更恰当的描述方法.

从密度矩阵的刘维尔方程出发,

在相互作用绘景中,密度矩阵表达为

式中下标 I 代表相互作用绘景.将(9)式展开到电场的一阶,得到:

从而得到电偶极矩的时间演化:

将电场的具体形式E(t)=ϵδ(t) 代入(11)式再次得到(7)式,这证明了两种描述方法的等价性.更高阶的响应函数可以通过保留(9)式更高阶的项得到:

因而,电偶极矩的更高阶近似为

2.2 刘维尔路径

2.1节给出了线性和非线性响应的具体形式,但仍然缺乏清晰的物理图像.为此,将(13)式展开并按项分析其中的物理含义.首先,以线性响应为例,将〈µ(t)〉(1)分解为两项:

仍然假设体系初始处于基态,即ρ(t0)=|0〉〈0|.R1(t,t1)代表密度矩阵的时间演化,可以用如下刘维尔路径描述:

以上的刘维尔路径描述了这样的过程:二能级体系在t1时刻与电磁场相互作用,密度矩阵由|0〉〈0|激发到|1〉〈0|,而后在H0下自由演化并积累相位因子,最终测量过程使|1〉〈0|重新回到|0〉〈0|,对该矩阵取迹即得到R1(t,t1) 的值.时间上的自由演化贡献了刘维尔路径的相位,最终的响应函数是所有刘维尔路径的求和(图1).Mukamel 等[1]将之类比为一种时间维度上的干涉行为.这种分析方法优点在于演化过程的物理图像清晰,可以容易地为体系添加衰减过程.体系的衰减过程分为两种:一种是激发态向基态自发弛豫,其特征时间为T1;另一种是体系的退相干过程,其特征时间为T2.在加入这两种过程后,R1(t,t1) 可以自然写作:

图1 线性响应的刘维尔路径.R1对应的路径为图中的红色实线Fig.1.Liouville paths of linear responce.R1 is illustrated by the solid red line.

在考虑进衰减过程后,线性响应的电偶极矩为

若系统是由激发能Δ各异的二能级体系组成的系综,由(17)式知,各二能级系统的电偶极矩正比于sin(Δt),因而大小不同、正负各异.系综平均使电偶极矩互相抵消,这种过程被称为“失相”.例如,假设Δ服从正态分布Δ～(Δ0,σ2),电偶极矩在系综平均后得到:“失相”过程贡献了衰减因子

以上的分析可以自然推广到非线性响应.由于(2)式具有粒子空穴对称性,因而二阶非线性响应为0.三阶非线性响应是领头阶的非线性响应.由(11)式知,体系的三阶非线性响应为

考虑图2所示的脉冲序列,E(t)=ϵ1δ(t)+ϵ2δ(t-t1)+ϵ3δ(t-t1-t2)).扣除(19)式中体系仅与一个脉冲和两个脉冲相互作用的部分从而保留体系与3 个光脉冲都相互作用的部分.用代表这部分贡献.

图2 一个典型的脉冲序列Fig.2.A prototypical pulse sequence.

式中取迹的部分可以展开成8 项,每一项分别代表一条刘维尔路径:

其中,µ0,µ1,µ2,µ3分别是µI(0),µI(t1),µI(t1+t2),µI(t1+t2+t3)的简写.不失一般性,在这里仅给出R1的分析过程:

沿用线性光谱的分析方法,以上刘维尔路径(图3)可做如下解释.体系在 0 时刻与电磁场相互作用,由基态|0〉〈0|激发到|0〉〈1|,此后在H0作用下自由演化积累相位eiΔt1并由此产生失相,同时退相干过程引入一衰减因子.在t1时刻,体系再次与电磁场相互作用并被激发到激发态|1〉〈1|.而后在激发态下演化,并由于体系弛豫引进一个衰减因子.在t1+t2时刻体系与电磁场相互作用,密度矩阵由|1〉〈1|转变为|1〉〈0|.之后自由演化,积累一个与失相过程符号相反的相位,此过程又被称为回相过程.与此同时,退相干过程产生另一个衰减因子.当t1=t3时,失相与回相的相位完全抵消,因而由失相引起的信号衰减现象被移除,这个效应又被称为“光子回波”.最终体系在t1+t2+t3时刻被测量,再次回到基态.这里的“光子回波”与Hahn[20]首先提出并由Carr和Purcell[21]改进的“自旋回波”过程类似,由 Abella 等[22]首次观测到.综合以上可得:

图3 非线性响应的刘维尔路径.R1由图中红色实线表示Fig.3.Liouville paths of non-linear responce.R1 is illustrated by the solid red line.

其中〈···〉Δ代表对整个系综做平均.其他各项均可以通过这种方法得到:

因为三阶非线性响应有3 个时间指标,可以选取任意两个进行二维相干光谱研究.由(23)式和(24)式可知,不同的衰减过程对应于不同的时间指标.退相干引起的信号衰减正比于t1+t3,自发弛豫带来的信号衰减正比于t2,失相带来的信号衰减正比于|t1-t3|.通过选取特定的时间通道,可以得到不同物理过程的特征时间.与之相比,线性光谱无法将3 种物理过程区分开.除此之外,也可以采用双边费曼图的办法来分析非线性响应过程,具体可参见文献[1,2].

作为回顾,本部分主要介绍了非线性响应理论,并通过分析刘维尔路径将响应函数中的时间变量与不同的物理过程联系起来.下面两部分将介绍这些分析技巧在具体例子中的应用.

3 一维横场伊辛模型的二维相干光谱[13]

Wan 等[13]考虑了一个典型的强关联体系—一维横场伊辛模型[23].其哈密顿量为

为了与(25)式中的σz区分开,在这里使用τz代表泡利矩阵.x方向的总磁矩在安德森准自旋表象下表达为

其中 sinθk=Δk/λk.若体系与x方向偏振的光脉冲耦合,那么哈密顿量写作

由于不同动量的安德森准自旋的算符之间互相对易,因而描述不同动量的二能级体系互相独立.该模型是二能级体系的系综,这也正是第2节中考虑的情形.但要注意,在第2节中偶极子算符仅有非对角元,与之相比,既有非对角元也有对角元.但这种不同只增加了刘维尔路径的数量,第2节中的分析方法仍然适用.

沿着这样的分析,详细计算线性响应、二阶非线性响应以及三阶非线性响应等各阶响应函数.线性响应的结果为

退相干衰减与失相导致的衰减都与t有关,因而无法被区分.在频率空间,线性谱描述了不同动量“自旋子”的连续谱,而体系退相干时间的信息都隐藏在连续谱内(图4(a)和图4(f)).

由于(28)式不具有粒子-空穴对称性,因而二阶非线性响应非零,其具体形式为

其中右侧第一项不随τ演化.在频域空间,这对应于wt=0并沿wτ方向分布的条带(图4(b)).条带随着T1的减小而增宽,因而可以来表征体系的弛豫时间.右侧第二项随着t+τ振荡,这说明体系经历失相过程,由此带来的衰减与退相干引起的衰减混合在一起无法区分,在频域空间对应于第一象限弥散的信号(图4(g)).

图4 横场伊辛模型铁磁相[h/(h+J)=0.3]的一维和二维光谱.从上到下,(a)-(e)无耗散的结果(1/T1,2=0);(f)-(j)有耗散的结果(1/T1,2=0.2(J +h));(k)-(o)添加无序后的结果.从左到右,每列分别是(ω),以及(t,t+τ),×(t,t+τ,t+τ),(t,t,t+τ)的傅里叶变换,以及沿着黑色箭头方向的信号轮廓(本图来自文献[13])Fig.4.One dimensional (1D) and two dimensional (2D) spectra in the FM phase [ h/(h+J)=0.3 ]of the TFIC.From the top to bottom,the rows show (a)-(e) the case with no dissipation (1/T1,2=0),(f)-(j) with dissipation (1/T1,2=0.2(J +h)),and(k)-(o) with quenched disorder.From the left to right,the columns show,respectively,(ω),and the FTs of (t,t+τ),(t,t+τ,t+τ),(t,t,t+τ),and its profile along a cut indicated by the black arrow.(This figure is reprinted from ref.[13]).

三阶非线性响应函数荷载更多信息.对于图5(a)表示的脉冲序列,非线性响应函数表达为

图5 脉冲序列 (a) (t3,t3+t2,t3+t2+t1)对应的三脉冲过程以及“自旋子”回波过程 对应的刘维尔路径;(b) 作为三脉冲极限的两脉冲序列下的三阶响应(本图来自文献[13])χ(3)Fig.5.Pulse sequences:(a) Three-pulse process associated with (t3,t3+t2,t3+t2+t1).The spinon echo process that produces the rephasing signal is also shown.(b) The χ(3) terms measured in the two-pulse setup are special limits of the threepulse process.(This figure is reprinted from ref.[13]).

其中

Wan 等将刘维尔路径分为4 类,前3 类中相因子没有抵消,所描述的刘维尔路径中没有回波过程.与之相比,第4 类中与t和τ相关的相因子符号相反,互相抵消,说明相应的刘维尔路径中存在回波过程.当t=τ时,失相和回相过程产生的相位完全抵消,这是第2节中提及的“光子回波”过程.由于体系中的准激发是自旋子,因而文献[13]中又称之为“自旋子”回波.

“光子回波”能够敏感地探测体系中的耗散效应.考虑一特殊的时间通道,对应于t3→t,t2→0,t1→τ(图5(b)).在此时间通道中,第4 类刘维尔路径存在回波过程,所以在频谱上位于第二、四象限,而其他3 类刘维尔路径均处于一、三象限或坐标轴附近.如此一来,二维相干光谱将“自旋子回波”同其他信号独立出来.进一步地,这种回波信号可以将退相干与失相引起的衰减区分开.具体来说,在频谱中,尖锐的峰分布在频谱的第四象限的对角线并组成一条条带,这是回波信号的重要特征.在不引入退相干过程时,这些峰沿反对角线方向的峰宽为零(图4(d)和图4(e)).在加入退相干过程后,反对角线的峰出现展宽(图4(i)和图4(j)),并且展宽与退相干过程的特征时间T2有关.与之相比,其他时间通道并不包含回波信号.比如考虑(t,t+τ,t+τ) (图5(b)),此时第3 类和第4 类刘维尔路径的相因子只和t有关.在二维相干光谱中,对应的信号呈分布在ωτ=0的一条条带.在不引入激发态向基态弛豫过程时,该条带沿ωτ方向峰宽为0 (图4(c)).在引入这样的过程后,沿ωτ方向出现展宽(图4(h)),并且展宽与弛豫时间T1有关.这些结果再一次印证了(23)式和(24)式后的分析.

除此之外,Wan 等还讨论了系统无序对结果的影响,作者将相互作用常数J和横场h设置为位置相关Jn和hn,并假定其满足正态分布.由于模型仍然可积,“自旋子”依然是体系的严格本征态,因而二维相干光谱中回波信号的反对角线方向宽度为零(图4(m)—图4(o)).这说明这种无序不影响二维相干光谱的主要结构,无法为体系引入耗散过程.

4 拉廷格液体的二维相干光谱[18]

第3节讨论了一个有能隙的强关联体系,但本质上该模型可以映射为二能级体系系综,进而使用双边费曼图等技巧进行分析.这样的体系在本质上并未脱离针对分子体系的非线性光谱的分析框架.为了讨论更一般的情形,Li 等[18]考虑了S=1/2 的XXZ模型(图6(a),(b)),这一模型在某些参数空间的激发是无能隙的,自然不能简化为二能级体系系综:

图6 (a) 法拉第构型示意图.磁场沿 z 方向.3 个圆偏振短光脉冲通过自旋模型,传播方向平行于 z 方向.第1 个光脉冲为右旋偏振,第2和第3 个光脉冲为左旋偏振.第1和第2 个光脉冲的时间间隔为 τ,第2和第3 个光脉冲的时间间隔为 tw,第三个光脉冲和测量时间的时间间隔为 t.(b) t≈τ 时,光子回波信号出现(本图来自文献[18])Fig.6.(a) The Faraday configuration.A magnetic field B is applied in the z axis.Three short electromagnetic pulses with circular polarizations pass through the S=1/2 spin chain.The propagation direction is parallel with the spin z axis.The first pulse is right-handed,whereas the second and the third are left-handed.The time delay between the first and the second pulse is denoted by τ,the second and the third by tw,and the third pulse and the time of detection by t.(b) When t≈τ,photon echo appears.(This figure is reprinted from ref.[18]).

其中 0,-1,-2,-3是时空坐标(0,0),(-t1,-x1),(-t2,-x2),(-t3,-x3)的简写.为了得到(35)式的表达式,首先利用约当-维格纳变换将XXZ模型映射为费米子模型[26].在费米子的表象下,自旋算符不是局域算符,因而(35)式仍难得到封闭表达式.在自旋液体相,进一步将费米子模型玻色化[27-30],忽略掉模型中重整化群意义下的无关项后,玻色子哈密顿量由拉廷格自旋液体描述:

其中的费米速度v和拉廷格参数K由模型的微观参数J⊥,Jz以及外磁场B共同决定.场算符的对易关系为

哈密顿量同样可以通过手征玻色场表达:

其中手征玻色场定义为

ϕL仅是x+=vt+x的函数,ϕR仅是x-=vt-x的函数.

在玻色子的表象下,自旋算符由一类特殊的顶点算子表达:

其中m是由外磁场B诱导的磁密度,x是格点j的空间坐标.注意到在铁磁情形下,自旋算符在玻色场表象中与磁化密度无关,因而沿z方向的外场不改变铁磁二维相干光谱的行为.与之相比,反铁磁情形的自旋算符在玻色场表象下与磁化密度有关,哈密顿量中的塞曼项会明显改变反铁磁二维相干光谱的行为.

将(40)式代入(35)式中,利用场论的相关技巧,得到(35)式的具体表达式.将结果代入(34)式并作适当近似,得到三阶响应函数χ+--+(t,tw,τ)的渐近形式.对于铁磁情形,他们发现χ+--+(t,tw,τ)的渐近形式为

其中,Δ=2K+1/(2K).相比铁磁情形,反铁磁χ(3)在时域是复的,这是因为磁化密度影响了自旋算符在玻色场中的表达.同样,数值积分的结果再一次确认了回波信号的存在(图8(a)—图8(d)).与横场伊辛模型不同,XXZ模型不能写成二能级系综的形式,因而二维相干光谱中的回波信号需要其他解释.

图8 (a),(b) 以 πTt, πTτ为自变量,反铁磁链的三阶非线性响应的实部和虚部.固定 πTtw=1,拉廷格参数 K=1,磁化密度 2mu/T=1.15.(c),(d) 分别是二维相干光谱的实部和虚部(本图来自文献[18])Fig.8.(a),(b) The real and imaginary parts of Nonlinear magnetic susceptibility of an antiferromagnetic chain as function of πTtand πTτ.The waiting time πTtw=1.The Luttinger parameter is K=1.The magnetization density2mu/T=1.15.(c),(d) The real and imaginary parts of the two-dimensional spectrum.(This figure is reprinted from Ref.[18]).

按照玻色化方法,磁化密度算符写作:

由基本的对易关系(37)式以及手征玻色场定义可知:

该刘维尔路径用“自旋子”产生湮灭图像理解如下:首先在 (-x1,-t1) 由S+(-1) 产生一对“自旋子”,在演化时间t后,这对“自旋子”在空间的位置分别为-x1-vt1(左行)以及-x1+vt1(右行).在 0 时刻,在原点位置由S-(0) 产生一对反自旋子.如图7(a),“自旋子”与“反自旋子”无法互相抵消,并且“自旋子”与“反自旋子”的距离随t1增加而增大,初态演化后与末态的交叠积分减小,因而该刘维尔路径随t1衰减.

图7 (a) 以 πTt,πTτ为自变量,铁磁链的三阶非线性响应.固定 πTtw=1,拉廷格参数 K=1.(b),(c) 分别是图(a)中数据傅里叶变换后二维光谱的实部和虚部(本图来自文献[18])Fig.7.(a) Nonlinear magnetic susceptibility of a ferromagnetic chain as function of πTtand πTτ.The waiting time πTtw=1.The Luttinger parameter is K=1.(b),(c) The real and imaginary parts of two dimensional spectrum,obtained by Fourier transforming the data of panel (a).(This figure is reprinted from Ref.[18]).

刘维尔路径中产生的“自旋子”与“反自旋子”全部抵消掉(图9(c)),所以这一刘维尔路径不随t衰减.Li 等将这一过程命名为“透镜效应”,这一效应正是回波信号的起源.Li 等进一步检验了这一物理图像,发现“透镜效应”对应的路径对响应函数贡献最大(图9(d)).反铁磁也可做同样分析,与铁磁不同,产生一对“自旋子”和一对“劳弗林(Laughlin)准粒子”[31].“透镜效应”在反铁磁情形中依然存在,因而二维相干光谱中依然有回波信号.

图9 (a) 两点关联函数中的“自旋子”产生湮灭过程.实线代表“自旋子”的动力学过程.虚线代表“反自旋子”的动力学过程.(b) 四点关联函数中的“自旋子”产生湮灭过程.(c) 四点关联函数中的“透镜效应”构型.(d) 在图(c)阴影部分的行为(本图来自文献[18])Fig.9.(a) The spinon creation/annihilation process in two-point correlation function.Solid and dashed lines represent dynamical processes of spinon and antispinon respectively.(b) The spinon creation/annihilation process in four-point correlation function.(c) The“Lensing”configuration in four-point correlation function.(d) The behavior of the shaded area in panel (c).(This figure is reprinted from ref.[18]).

当玻色场存在色散时,空间中的“自旋子”波包随时间扩散开.“透镜效应”说明,波包扩散导致“光子回波”衰减.具体来说,τ很大时,由S+(-3)产生的左行自旋子的波包在-t2时刻会非常平坦.这一扩散开的波包无法与由S-(-2) 产生的左行反自旋子互相抵消.这使得“透镜效应”不再完美,从而导致“光子回波”信号的衰减.联系第2节的讨论,热弛豫引起的信号衰减仅与tw有关,而色散引起的“光子回波”衰减仅与t和τ有关.二维相干光谱将两种衰减方式区分开来.

我们已经注意到回波现象对于研究强关联电子体系的重要作用.按照Mukamel[1]的物理图像,二维相干光谱中的回波信号来源于不同刘维尔路径的时间维的干涉作用.“透镜效应”将这一图像进一步发展.回波现象实际上描述了准激发的波包之间的时空间干涉作用.由电磁场激发的准激发以波包的形式在空间运动,不同波包之间互相干涉.当满足干涉最大条件时(准激发的波包完全重叠),二维相干光谱对应的响应强度最大.这一图像自然解释了XXZ模型中的回波现象并预言了色散对回波的影响.这种色散导致回波衰减的效应难以被传统的分析方法如双边费曼图所描述,因而这种新图像涵盖了更多物理内容.

5 总结与展望

本文通过本研究组的两个工作,试图说明二维相干光谱在强关联电子体系的重大潜力.这些理论工作为非线性谱学提供了新的概念与方法,对未来的研究具有很大启发性.与此同时,强关联体系的二维相干光谱作为一个理论上崭新的方向,也充满未知与挑战.首先,目前的理论工作局限于一些可以严格求解的解析模型,并没有一套成熟的微扰方法对一般体系进行研究.此外,大部分理论工作集中于强关联系统的磁响应,相比电响应,磁响应实验难度更大[19].最后,仍然没有一个确定性的实验验证理论预言的结果.未知也意味着可能性,能否通过非平衡场论[32]发展出一套微扰方法? 能否将这样的方法应用到超导等体系中获得诸如库珀对的寿命等信息? 能否利用一些半经典方法[33-35]理解二维相干光谱? 我们相信,未来对这些问题的研究会揭示强关联电子体系中更丰富的物理信息.

感谢与Peter Armitage和Masaki Oshikawa 的有益合作.