从横场伊辛链到量子E8可积模型*

王骁 杨家豪 吴建达2)

1) (上海交通大学李政道研究所,上海 201210)

2) (上海交通大学物理与天文学院,上海 200240)

本文综述了近年来关于横场伊辛链模型及量子 E8可积模型研究工作的一系列理论及实验进展.在针对横场伊辛链的研究工作中,理论上发现格林艾森比率 (热或磁胀系数与比热之比) 在此模型中独特的奇异量子临界行为,并实质性地扩展了能涌现横场伊辛链普适类的微观模型.这些理论进展成功地推进了一系列合作实验在准一维反铁磁材料BaCo2V2O8及SrCo2V2O8中首次实现横场伊辛链普适类.在针对量子 E8可积模型的研究工作中,理论上严格计算了该系统的低温局域动力学行为及在零转移动量之时的动力学结构因子,并在动力学结构因子的连续谱区域得到了级联幂次发散的奇性谱边激发.这些理论进展在结合详细的量子临界标度行为分析及大规模iTEBD 数值计算之后,成功促成了包括太赫兹谱学测量、非弹性中子散射、核磁共振等系列实验,在BaCo2V2O8中首次实现量子 E8 可积模型.量子 E8可积模型的物理实现为在真实材料中研究量子可积系统的物理开拓了新的边界.这一系列关于横场伊辛链及量子 E8可积模型的研究进展为量子不可积系统的研究带来了新的契机,并将启发凝聚态系统、冷原子系统、统计场论和共形场论等相关方向的研究.

1 引言

量子相变指在零温下通过调节参数,强关联多体系统发生的相变行为[1].在量子临界点(quantum critical point,QCP) 附近,强烈的量子涨落往往会涌现出奇异的量子临界行为[2-7].对这些量子临界行为的深入理解将有助于解答长期悬而未决的物理问题[8-16],包括非常规高温超导体、非费米液体等[17-19],并有助于在材料中实现量子自旋液体[20,21],因此吸引着理论和实验的持续研究.量子临界行为能被临界指数所刻画,每一组量子临界指数对应于一种普适类,并往往伴随相应的低能有效量子场论[1].即使微观起源不同,属于同一普适类的模型在QCP 附近也会展示出非常一致的物理,如标度规律、关联行为及动力学激发等[1].因此,对普适类的深入研究可以使人们同时理解QCP 附近大量具有该普适类的微观模型,而这意味着将实际系统调节向QCP 附近是实现并研究普适类之必要.然而,由于真实材料内在的复杂性,只有极少数的QCP 可以被实验实现,即使一些形式上并不复杂的微观模型,如横场伊辛链模型[22](transverse field Ising chain model,TFIC模型,见方程(1)),实验上也一直难以找到候选材料实现并研究其对应的普适类.

在QCP 处对体系的微扰可能会进一步涌现出丰富的物理,如涌现对称性和涌现可积性(涌现对称性指量子系统在低能物理中的对称性大于完整的哈密顿量的对称性.涌现可积性指不可积系统在某些特定的参数区间成为可积系统,产生了可积性;或可积系统在恰当扰动下又呈现出新的可积性)[3-6,23-26].自1989 年起Zamolodchikov等[23,24,27-29]发展了微扰共形场论(perturbed conformal field theory,pCFT).他们通过在不同普适类的临界点处施加微扰,发现了一系列能够由例外李代数(exceptional Lie algebra)描述的可积模型[24].比如,QCP 处的TFIC 模型在纵向磁场的微扰下会涌现出量子E8可积模型[23,24,27-31].这类涌现现象同样具有普适性,一方面不依赖于微观模型的具体细节,另一方面拓展了量子可积模型范围,有利于推动实验在实际材料中对可积模型进行研究.

近年来,吴建达及其合作者在理论和实验的共同努力下,在强关联材料BaCo2V2O8(BCVO)和SrCo2V2O8(SCVO) 中首次成功实现TFIC 模型的普适类[22,32-35],并在BCVO 材料中首次实现量子E8可积模型[30,36,37].这一系列进展不仅是量子可积磁性理论和实验上的突破[2],也展示了对量子临界行为、普适类、涌现对称性和涌现可积性研究的重要意义,更开拓了于真实材料中实现并研究量子可积模型的重要方向.本文按照历史的脉络,从理论和实验两个方面综述了这一系列的进展.全文组织如下:第2节分别回顾TFIC 模型和横场海森伯-伊辛链的量子临界行为及其普适类;第3节介绍证实BCVO 材料和SCVO 材料中TFIC 模型普适类的实验进展;第4节概括量子E8模型的理论背景和进展;第5节叙述量子E8模型在材料中的实验实现;第6节为总结与展望.

2 伊辛普适类

2.1 横场伊辛链

TFIC 模型的哈密顿量为[1,22]

其中J及gJ分别为铁磁及横场与自旋间的相互作用强度.=/2,α=x,y,z表示在格点i处的自旋 1/2算符,是泡利矩阵.该模型是一精确可解模型,利用约当—魏格纳变换 (Jordan-Wigner transformation)和波格留波夫 (Bogoliubov)变换可以对角化该哈密顿量并求解出其激发谱[1,22,32]:

1970 年,Pfeuty[22]首次研究了TFIC 模型的量子相变.1984 年,Zamolodchikov 等[3,30,38]证明了当g=gc时TFIC 模型的低能物理可以被中心荷 1/2 的共形场论 (central charge 1/2 conformal field theory,CFTc=1/2)所描述.在标度极限下[30],TFIC 模型中的自旋算符及分别对应于CFTc=1/2中自旋密度算符σ(x) 及能量密度算符ε(x).当|x|→0 时,其两点关联函数具有标度行为[31,38]

σ(x)和ε(x)对应的标度维度分别是2Δσ=1/8和2Δε=1[30,31,38],它们刻画了TFIC 模型的普适类,同时也刻画了二维经典伊辛模型的普适类[39](以下简称该普适类为伊辛普适类).然而,自TFIC模型提出后的近半个世纪里,具有伊辛普适类的强关联材料却一直没有被确认存在,这是因为一方面实验上受限于实验精度,另一方面缺乏有效的理论指导去分辨并证实观测到的量子临界行为是否属于伊辛普适类.比如,方程(3)所提及的标度行为是等时短程关联,实验上不具有可操作性,无法进行测量.2018 年,吴建达与合作者[32]关于伊辛普适类格林艾森比率 (Grüneisen ratio)量子临界标度行为的理论突破为实验提供了一个可靠的判定标准.这里先简要归纳伊辛普适类在QCP 附近的标度规律,这些标度规律皆可作为判断伊辛普适类的依据[1,22,32].首先,QCP 附近系统的磁化Mz满足Mz≈|g-gc|δ,其中δ=1/8.其次,在QCP 附近考虑自旋关联函数

该关联随格点间距呈指数衰减,即Cij≈e-|i-j|/ξ,其中ξ被称为关联长度,满足标度律 1/ξ≈|gc-g|ν,ν=1.此时系统激发的能隙满足标度律Δ≈|gc-g|zν,z=1 被称为动力学临界指数[1,32,35].

由于临界指数也可以直接用于分析热力学量临界行为,因此在实验方面,可以通过测量相应热力学量的量子临界行为考察相关材料的普适类性质.比如考虑和磁热效应相关的热力学量格林艾森比率:

其中H是可控外场;cH,αH,M和S分别为摩尔比热、磁膨胀系数、每摩尔磁化强度和熵.在TFIC模型中,相关物理量均可被解析地求解出来[32].零温极限下,格林艾森比率在横场趋于TFIC 模型的QCP 时趋于发散而在量子临界区内,趋于QCP 时,格林艾森比率则趋于常数[32].此奇异的量子临界行为一方面表明了TFIC 模型的QCP 是本性奇点,另一方面则可作为TFIC 模型中伊辛普适类的独特特征,为实验探测伊辛普适类提供有效的理论指导.值得一提的是,格林艾森比率的发散行为意味着熵极大值会出现在QCP 附近[40],这在相关材料的量子临界测量中得到确认,如Sr3Ru2O7和CeCu6-xAux等[41,42].于是,基于上述理论,在实验中可以通过测量相关材料格林艾森比率的量子临界行为,判断该材料的普适类是否是伊辛普适类.

2.2 横场海森伯-伊辛链

接下来总结横场海森伯 -伊辛链在外加横场之下的可能量子相变及其普适类特征,首先考察外加均匀横场之下的海森伯-伊辛链:

这里Hx是自旋和横场的耦合系数.ε通常取0—1之间,当ε=0,1 时该模型分别对应TFIC和横场海森伯链,于是称该模型为横场海森伯-伊辛链[43,44].此前已经证明了该模型在加横场后会存在伊辛普适类[44].2019 年,吴建达及合作者[35]利用大规模infinite time-evolving block decimation (iTEBD)数值计算证明了引入垂直于所加均匀横场的交错横场后,该模型在量子相变点处依旧展示了伊辛普适类.考虑一般的带有交错横场的海森伯-伊辛链哈密顿量

图1 TFIC 模型相图.g 为横场参数,gc=1 为QCP.在相变点左侧,青蓝色实线代表零温下的铁磁相;在相变点右侧,黄色实线代表零温下的顺磁相.相图中的两斜虚线为低温下各个不同无序区域的过渡边界,而上方虚线则是量子临界区域到经典区域的过渡边界Fig.1.A phase diagram of TFIC,where g labels transverse field,and gc=1 is the QCP.The cyan solid line represents for ferromagnetic phase at zero temperature on the left of the critical point,while on the right the yellow solid line stands for a paramagnetic phase at zero temperature.The two tilted dotted lines illustrate the crossovers of different disorder phases in low temperature region,while the dotted line above shows a crossover from quantum critical region to classical region.

利用iTEBD 数值方法,在gx-gy-ε三维参数空间内可以精确计算该有效模型中展现伊辛普适类物理的QCPs,所有的QCPs 组成了gx-gy-ε参数空间内一量子临界曲面,如图2[35]所示.

图2 由iTEBD 计算得出的针对模型方程(7)的量子临界曲面 (图片经文献[35]允许转载,版权归2019 IOP Publishing Ltd 所有)Fig.2.Quantum critical surface calculated by iTEBD algorithm for Eq.(7).(Reprinted with permission from Ref.[35].Copyright 2019 IOP Publishing Ltd.).

该理论工作进一步考察了带有交错横场的海森伯-伊辛链存在沿z方向的4 周期磁场微扰项的情形,发现这一微扰项仅轻微移动QCP 的位置而不改变相应的伊辛普适类.计算结果见图3[35].

图3 (a)交错磁化 Mz(g)的计算结果.蓝线和红线分别代表了带有和不带有4 周期项的数据.这两条曲线可以分别用Mz=0.524(gc1-g)0.126和Mz=0.530(gc2-g)0.128来拟合,其中 gc1=0.1454,gc2=0.1456,在误差精度范围内均可得到临界指数为 1/8的结论.内嵌图是利用对数坐标轴画出的 Mz的标度行为.(b)半对数坐标轴下纠缠熵 SEE(l)和链长 l 的关系,两种情况都符合伊辛普适类的 1/2中心荷.(c)当 g=0.1448时,自旋两点关联函数和距离 i-j间的函数关系.内嵌图展示了 ln C(i-j) 在长程时与 i-j成正比.(d)关联长度的倒数和g的函数关系,两种情况都符合伊辛普适类中关联长度指数 ν=1 的结论(图片经文献[35]允许转载,版权归2019 IOP Publishing Ltd 所有)Fig.3.(a) iTEBD results for staggered magnetization Mz(g).The blue line and red line represent iTEBD data with and without the four periodic term,respectively.The two curves can be fitted with Mz=0.524(gc1-g)0.126and Mz=0.530(gc2-g)0.128,where gc1=0.1454,gc2=0.1456,and critical exponent δ=1/8 is obtained within error bar for both cases.Inset shows scaling behavior of Mzin log-logplot.(b) Entanglement entropy SEE(l)versus length lin a semi-log plot,both fall into central charge c=1/2of TFIC universality.(c) Spin-spin correlation function versus distance i-jat g=0.1448.Inset shows ln C(i-j)being proportional to i-jin long range region.(d) The inverse of correlation length in terms of g,both agree with correlation length exponent ν=1.(Reprinted with permission from Ref.[35].Cop yright 2019 IOP Publishing Ltd.).

综上所述,该工作利用iTEBD 数值算法证明了在带有交错横场的海森伯-伊辛链中通过合理调整所加磁场,可以到达伊辛普适类的量子临界区,从而实现伊辛普适类的物理.而由于BCVO和SCVO 材料均可被海森伯-伊辛链所描述[33-37],理论上预期这两种材料在横场调节之下可以实现伊辛普适类.第3节中将叙述在BCVO和SCVO材料中实验实现伊辛普适类的一系列工作.实验测量所得到的BCVO 材料中格林艾森比率的量子临界行为完美符合了伊辛普适类中该物理量量子临界行为的理论结果,从而为BCVO 在外加强横场之下涌现伊辛普适类提供了坚实的判定依据[33].而利用核磁共振(NMR)实验辅以数值算法,通过对SCVO 材料及其有效模型进行研究,我们证实了SCVO 材料在弱横场下展现了伊辛普适类[34].

3 BCVO和SCVO 材料中的伊辛普适类

利用上述伊辛普适类的相关理论,吴建达及合作者[32-37]对BCVO和SCVO 材料展开测量并证实了这两类准一维材料在外加横场之下能够实现伊辛普适类.BCVO和SCVO 是能被海森伯-伊辛链有效描述的准一维反铁磁材料,在外加磁场下一般性的有效哈密顿量为[43]

其中µB是玻尔磁子,H0为外加磁场,代表g张量:

考虑加横向磁场H0=H0ex,由于g张量的非对角元,此时在Y轴和Z轴会诱导出垂直于该磁场的有效磁场,因此有[45]

其中

而ξ轴、ψ轴和ζ轴代表将实验室的XYZ坐标轴绕Y轴转θ角,绕Z轴转ϕ角后得到的材料中的局域坐标轴.CoO6的晶格结构、螺旋链结构及坐标轴间的关系如图4所示.

图4 (a) CoO6八面体的结构;(b) CoO6四周期螺旋链的单元结构及局域坐标系ξψζ;(c)局域坐标系ξψζ 与实验室坐标系XYZ 的关系Fig.4.(a) Structure of CoO6octahedron;(b) a unit structure of CoO6four-period screw chain and the local coordinate ξψζ ;(c) relation between the local coordinate ξψζ and the lab coordinate XYZ.

对于外加磁场沿[110]和[100]两种情况,ϕ1分别取 0°和45°.以BCVO 为例,若对该材料施加沿[110]的磁场,则gxy=gxz=0,可以避免诱导出有效交错横场.在低温下对该材料进行磁热效应测量,可以得到随温度和横场变化的格林艾森比率数据.对该数据的标度行为分析完全符合理论上对TFIC 模型格林艾森比率量子临界行为的预期,从而确认该材料QCP 在=40 T 附近,并在该点附近涌现出伊辛普适类(图5)[33].

图5 (a) BCVO 材料在加[110]方向磁场后的相图,利用从不同起始温度出发的绝热磁热测量得出.红色区域内直到 20 T 处标记了一有限温的三维序.一维QCP 出现在 =40 T 附近,上方黄色区域为量子临界区.(b),(c)以磁场和温度作为自变量的格林艾森比率实验数据.在图(b)中,大于QCP 时不同温度的数据展现了 ΓB≈(B-)-1 的发散行为.相对应地,当场强在临界场强附近时,低温下格林艾森比率 ΓB(T) 趋于收敛(图片经文献[33]允许转载,版权归2018 American Physical Society 所有)Fig.5.(a) Phase diagram of BCVO material with a transverse field along [110],obtained from adiabatic magnetocaloric-effect measurements starting from different temperatures.Red area labels a finite-temperature three dimensional (3D) order until 20 T.The one dimensional (1D) QCP appears around =40 T,while the yellow area above labels quantum critical region.(b),(c) Experimental data of Grüneisen ratio in terms of magnetic field and temperature.In panel (b) data of different temperature above the QCP shows divergent behavior of ΓB≈(B-)-1.Correspondingly,with the field being around the critical field,Grüneisen ratio converges in low temperature.(Reprinted with permission from Ref.[33].Copyright 2018 American Physical Society).

然而,40 T 左右的量子临界磁场过于昂贵,难以从此出发继续深入研究与伊辛普适类相关的物理.此后不久,吴建达及合作者[34-37]在对SCVO和BCVO 的进一步研究中发现,沿[100](或[010])方向施加一横向磁场所诱导出的沿[010](或[100])的有效交错横场能极大降低达到QCP 所需要的磁场强度,相关的数值分析已在2.2节中提及(图2)[35].这里介绍吴建达及合作者对施加沿[100]方向磁场的SCVO 材料的伊辛普适类QCP 的实验进展.此时,该材料的有效哈密顿量为[34]

其中J≈7.5meV,各向异性因子ε=0.47.gx≈2.79为旋磁比率,诱导出的有效交错横场Hy≈0.4H,四周期场Hz≈0.14H.利用NMR 实验可以测量此时的QCP 及其普适类,测量结果见图6[34].在SCVO 中,沿[100]方向调节磁场会涌现出两个QCPs,记为Hc1和Hc2.其中,Hc1≈7.03 T 代表三维序的终点,而Hc2≈7.7 T 对应伊辛普适类的量子临界点[34].

图6 外加[100]方向磁场的SCVO 材料相图.图中AFM,QCR和QD 分别代表反铁磁相、量子临界区和量子无序区.两个QCPs分别为 Hc1≈7.03 T和Hc2≈7.7 T.蓝色线可以被曲线 ～|Hc2-H| 所拟合,揭示了该QCP 对应伊辛普适类.内嵌图是两QCPs 附近拟合曲线的放大 (图片经文献[34]允许转载,版权归2019 American Physical Society 所有)Fig.6.Phase diagram of SCVO with transverse field along[100].The AFM,QCR and QD in the figure corresponds to anti-ferromagnetic phase,quantum critical region and quantum disordered region,respectively.The two QCPs are Hc1≈7.03 T and Hc2≈7.7 T.The blue line can be fitted with ～|Hc2-H|,implying the Ising universality of the 1D QCP.Inset is an enlarged figure of the fitting curves near two QCPs.(Reprinted with permission from Ref.[34].Copyright 2019 American Physical Society).

这篇工作表明,对于能被海森伯-伊辛链有效描述的磁性材料而言,若外加横向磁场时能诱导出新的平面内交错磁场,这些诱导场能够有效压制达到临界磁场所需要的的场强大小,从而在一个较小的磁场环境下实现伊辛普适类.

本节总结了在BCVO和SCVO 中实现伊辛普适类的实验工作[33,34].在实现伊辛普适类之后,笔者及合作者结合理论工作进展[30,46],推动了量子E8可积模型(以下简称E8模型)的物理实现.下文将首先介绍与该模型相关的理论背景及研究进展[30,46-50],再详细描述理论及实验共同合作促成的该奇异模型在BCVO 材料中的实现.

4 量子E8可积模型

1990 年前后,Zamolodchikov[23]和Fateev 等[24]证明了对应于某一特定中心荷的CFT,若它被一相关元场所微扰,该CFT 对应的哈密顿量会成为能够被例外李代数所描述的可积模型[16,27-31].特别地,对于中心荷c=1-6/[p(p+1)]的CFT 哈密顿量Hp,其被相关元场Φp微扰的哈密顿量为

当p=6,4,3时,该微扰哈密顿量可以分别被E6,E7,E8例外李代数所描述,对应的系统分别包含有6,7,8种准粒子激发[24,51-57].在QCP处,TFIC模型是中心荷c=1/2 的CFT[3],在标度极限下,自旋密度算符为相关元场[30,31],在其微扰下的哈密顿量

成为一个激发谱能够被E8例外李代数所描述的量子可积模型[16,27-38],其对应的格点模型为被一纵向磁场所微扰的QCP 处(g=gc=1)的TFIC 模型[30,36,37]:

E8模型包含有8 种不同的准粒子激发.其中最轻的准粒子质量m1=Ch8/15,C=4.4049···[58].而次轻准粒子质量m2与m1满足黄金分割,即m2=2 cos(π/5)m1≈1.618m1.其他准粒子的质量都可以通 过m1及m2精确地得到.E8模型相图如图7所示.该相图给出了E8准粒子激发与伊辛普适类的联系:在TFIC 模型的伊辛普适类处,QCP 附近的物理可以用中心荷 1/2 的共形场论描述.此时纵向磁场的微扰提供了一个微弱的禁闭效应,引起一系列完美共振的介子态,从而产生8 种E8模型的准粒子激发[59,60].因此纵向磁场较小的区间正是可以涌现出量子E8模型的区间,如图7中蓝色区域所示.而这些准粒子间的精确质量关系亦见于图7.

图7 E8模型相图.图中蓝色实线代表了 E8模型涌现的区域,其中 g和hz分别对应于横场参数 g和方程(17)中的hz.其余E8准粒子的质量和m1和m2的关系列于右上,其中 m1为最轻准粒子的质量Fig.7.A phase diagram of E8 model.The blue solid line in the figure implies region of E8 model emerging,whereg and hzcorrespond to transverse parameter gand hz in Eq.(17) respectively.Relations of other E8 particles'masses to m1 and m2 are listed in the top right-hand corner,where m1 is mass of the lightest particle.

关于FF 的一般性讨论可参考文献[30,31,61-64].而对于E8模型FF 的求解则可见于文献[30,47],这里不再赘述.

在方程(20)中,通过引入不同粒子数态,依粒子数的不同,总体的贡献可分解为不同的通道,即单粒子通道、两粒子通道、三粒子通道和四粒子通道[30].由于质量最大的E8粒子的质量约在 4.8m1左右,而计算表明随着能量增加谱权重急剧衰减,因此选取DSF直至ω=5m1.以Dyy(ω)为例,其DSF 谱如图8所示.这一解析计算的结果展示了丰富的物理:首先,对于单粒子通道,所有的谱均为共振峰,位置分别对应了E8模型中8 种准粒子的质量.随着能量的增大共振峰的强度不断衰减.这是由于在粒子质量增大的同时,单粒子形式因子的绝对数值也在不断减小,而单粒子通道的谱权重的峰值正比于|F|2/m,因此质量越大的单粒子对应的峰值越小.其次,对于两粒子通道,含不同类型粒子a和b的通道谱存在谱边发散行为,位于这两粒子的质量之和边界ω0=ma+mb处.这是由两粒子态密度的发散行为所导致的,即范·霍夫奇性 (Van Hove’s singularity)[30].两相同类型的粒子组合却不存在这种发散行为,这是由于它们的FF 特殊结构恰好能够消除由态密度导致的发散.最后,三粒子和四粒子通道的谱贡献不存在上述的发散行为,且二者行为相似,贡献都比较小.这是由于随着引入粒子数的增多,FF的绝对大小在不断衰减,同时由于相空间的增大,奇性行为被消去.从图8的解析数据可以预期,随着能量的不断升高,高能区的DSF 谱会不断衰减,其行为相对来说较为平凡.将解析计算得出的DSF谱线以实验测量的能量分辨率作展宽后[30],便可以与实验进行对比以验证具体材料中是否实现E8模型.

图8 Dyy(ω)的解析计算结果 (a)单粒子通道的谱权重贡献,8 个峰的位置依次对应 E8模型8 种粒子的质量,其中 m1=1 ;(b)—(d)二粒子、三粒子和四粒子通道的谱权重贡献,图中分别标出了谱线和峰的具体贡献来源(图片经文献[30]允许转载,版权归2021 American Physical Society 所有)Fig.8.(a) Spectrum contribution from single particle channel,the eight delta peaks one-by-one correspond to the masses of the eight particles in the quantum E8 model with m1=1 ;(b)—(d) contributions from two,three and four particle channels.The detailed contributions of different peaks and curves are highlighted in the sub-figures.(Reprinted with permission from Ref.[30].Copyright 2021 American Physical Society).

5 BCVO 材料中的E8激发

2010 年,牛津大学的Coldea 实验组[65]利用非弹性中子散射技术(inelastic neutron scattering,INS)对准一维磁性材料CoNb2O6展开了自旋动力学谱的测量.其测量结果提供了量子E8模型可在材料中实现的初步实验证据,这体现在质量最轻的两个粒子质量比满足黄金分割.然而,受限于实验分辨率,其测量结果未能得到较为精确的连续谱,无法就连续谱区域的相关的E8物理做出有效的比对,从而未能就该材料是否能完整实现量子E8模型给出确定性的回答.之后,2020 年对该材料的太赫兹 (THz) 谱测量[66]得到了详细且精确的连续谱,实验所展示的连续谱区域特征和解析上得到的E8谱存有明显的偏差,因此该材料是否能实现量子E8模型有待进一步深入仔细的研究[67,68].接下来聚焦于 BCVO 材料中实现量子E8模型的工作.

前文所述的E8模型对应的格点模型为铁磁相互作用,而若想在反铁磁相互作用下实现E8模型,需要将方程(17)中的纵场转为交错纵场,这在凝聚态系统实验中几乎不可能通过外加手段实现,因此只能另辟蹊径.考虑沿[010]方向的外加磁场,材料的有效哈密顿量为[34-37]

图9 (a) BCVO 的晶格结构.(b) 加沿[010]方向磁场后BCVO 的相图,QCP 为 ,隐藏在三维序中,其一维量子临界行为可在三维序外测得.棕色圆圈代表由NMR 实验测出的尼尔(Néel)温度 TN.=(10.4±0.1)T,=(4.7±0.3)T.E8模型在蓝色缎带区涌现.(c) 8 种 E8粒子的质量以 m1为单位沿能量轴分布(图片经文献[37]允许转载,版权归2021 American Physical Society 所有)Fig.9.(a) Crystal structure of BCVO.(b) Phase diagram of BCVO with a transverse field along [010]direction.QCP is at ,hidden in the 3D order.Its quantum criticality can be measured outside the 3D ordering phase.Brown circles represent for Néel temperature TN measured by NMR experiments.=(10.4±0.1)T, =(4.7±0.3)T.E8 model emerges in the blue ribbon area.(c) Masses of eight E8 particles along energy axis in unit of m1.(Reprinted with permission from Ref.[37].Copyright 2021 American Physical Society).

图10 三维反铁磁序外NMR 实验测量结果 (a) 在不同横场下测得自旋-晶格弛豫率 1/T1关于温度的函数.下箭头标记出决定TN的 1/T1的峰.(b)拟合 1/T1得到能隙,温度区间为6—12 K.具体细节见文献[37].QCP 位于 (4.7±0.3) T.内嵌图是对数坐标下放大的低温区数据.灰色直线区域内的量子临界行为可以拟合为 1/T1～T-0.75 (图片经文献[37]允许转载,版权归2021 American Physical Society 所有)Fig.10.NMR experiment result outside the 3D anti-ferromagnetic ordering phase.(a) Spin-lattice relaxation rate 1/T1 with different transverse fields as a function of temperature.The down arrows label 1/T1 peaks that determining TN.(b) Energy gap by fitting 1/T1,with temperature being from 6 K to 12 K.The details can be found in Ref.[37].QCP is at (4.7±0.3) T.Inset:a log-logplot of enlarged data in low-temperature region.The quantum critical behavior around the gray line area can be fitted by 1/T1～T-0.75.(Reprinted with permission from Ref.[37].Copyright 2021 American Physical Society).

图11 (a) Q=(002)QCP 处INS 测量数据,场强 H=4.7T,温度 0.4 K.蓝色方块及误差棒为实验测量数据,黑线是高斯函数的拟合.红色竖线对应于8 种 E8模型的粒子.其他的峰来自于多粒子通道的连续谱贡献以及布里渊区折叠效应(标记于 F1和F2处).(b) 解析计算的 Dxx(ω),数据经过 0.08m1洛伦兹展宽.m1=1.2 meV.红色线和黑色线分别代表单粒子激发贡献及引入多粒子激发贡献的总和.(c) iTEBD 数值计算有效模型在转移动量 q=0 处的激发谱,其中蓝线和黑线代表不包含及包含折叠效应的数据.(d)—(g) 单粒子到四粒子激发通道的贡献,ijkl代表激发来自于 mimjmkml 通道(图片经文献[37]允许转载,版权归2021 American Physical Society 所有)Fig.11.(a) INS data obtained from Q=(002)at QCP with H=4.7T,temperature 0.4 K.The blue diamonds with error bars are experimental data,the black curve is fitted with Gaussian functions.The red vertical lines correspond to eight particles ofE8model.Other peaks are contributed from multi-particle channels and zone-folding effect (labelled at F1 and F2).(b) Analytical result of Dxx(ω)with a 0.08m1 Lorentzian boradening.m1=1.2 meV.Red curve and black curve are single particle spectrum and total spectrum including multi-particle excitations respectively.(c) Effective model spectrum of iTEBD calculation at q=0,where blue curve and black curve are results of without and with zone-folding effect.(d)—(g) single particle-four particle excited channels' contributions,where ijklrepresents for mimjmkml channel.(Reprinted with permission from Ref.[37].Copyright 2021 American Physical Society).

6 结论

本文综述了在BCVO和SCVO 材料中首次发现伊辛普适类并在BCVO 材料中首次实现E8模型激发的研究进展.叙述了便于判断伊辛普适类的格林艾森比率的奇异量子临界行为.随后介绍了横场海森伯-伊辛链并叙述了数值验证其具有伊辛普适类的理论进展.通过将数值方法和实验相结合,确认了BCVO和SCVO 材料在外加横场之下会涌现出伊辛普适类.在此基础上,介绍了E8模型主要的理论进展及在理论与实验的共同推动下于BCVO 材料中实现E8模型的一系列工作.

这些研究进展充分展示了普适类的强大力量,其所蕴含的迷人的物理现象持续吸引着人们对量子临界物理及相关普适类物理的进一步的探索[2].我们期待未来通过更多的凝聚态理论和实验物理学家的合作,能在实际材料中实现更多非平凡的普适类,并能通过深入研究这些奇异系统中的非平凡物理激发,最终实现物理操纵这些奇异的激发并期望能够服务于新颖量子器件的设计与开发.