形成与自组织12 分钟

自旋链中的局部扰动,怎样分解成能分开传播的激发?

局部改变一对自旋,可以产生两个能够分开的集体边界;它们的运动对象和量子数不必与最初操作一一对应。本文用双畴壁有效模型追踪传播与束缚,计算各向同性链的两自旋子连续谱范围,再比较YbAlO3在低温有序与高温散射中的表现。谱的连续宽度、寿命展宽、动态相干和纠缠见证分别讨论,并用一个可精确求解的热态例子说明,见证未达门槛不等于没有纠缠。

局部改了两处方向,为什么后来可以看见两个远离的对象

一串相互作用的自旋,并不一定把局部扰动当成一个整体往前搬。最初只改变一小块,随后真正容易传播的对象,可能是这块改变与周围背景之间的两个边界。

自旋子是这种集体激发的一个重要例子。它的“分数化”,说的是有效激发携带的量子数,并不是把电子、原子或一个实体小磁针从中间切开。

本篇先用一个能逐步计算的极限看清图景,再进入近各向同性的自旋1/2链和真实中子散射。最后讨论一个前沿问题:为什么某些动态谱结构能在相对交换能标很高的温度下继续存在,而另一种量子指标却已经很小。12

先看一个容易画清的极限:两个畴壁

设相邻自旋倾向于反向排列,形成上下交替的背景。现在把一对相邻方向,从“下、上”改为“上、下”。中间一小段与原来的背景错开了相位,它的左右两端出现两个边界,称为畴壁。

强各向异性极限下的基底图景。局部改变相邻两个方向后,出现两个不符合原交替模式的边界。标出的半单位为按相邻两点平均定义的自旋投影;这不是近各向同性链真实量子基态的完整图片。
强各向异性极限下的基底图景。局部改变相邻两个方向后,出现两个不符合原交替模式的边界。标出的半单位为按相邻两点平均定义的自旋投影;这不是近各向同性链真实量子基态的完整图片。

如果两种交替排列具有相同的体能量,那么把中间那段延长,并不会不断增加“每一段都站错队”的代价;主要代价可以集中在两个边界上。

加入能交换相邻自旋方向的量子项后,边界有机会移动。局部扰动于是能够分解成两个可传播的集体对象,原来的原子仍然留在晶格位置。

这里还可以具体说明“半单位”。取自旋投影为$\pm1/2$,对相邻两点求平均。在完全交替的背景中,这个平均为零;在两个相同方向构成的边界处,则分别为$+1/2$或$-1/2$。本例两个边界携带相反投影,总改变仍为零。

这个平均不会凭空改变总的自旋投影:固定链端后,把相邻平均相加,仍得到同样的总改变。半单位集中在边界,是这套背景与归一化下的守恒量分配。

这是一个便于理解的强各向异性极限。近各向同性Heisenberg链的基态本身具有强量子涨落,其自旋1/2激发由多体量子理论定义,不能只用一张经典交替图代替。1

把传播算出来:分离是否需要越付越多能量

本库建立了一个双畴壁有效模型。基底记录左、右畴壁的位置;相邻自旋交换,使每个畴壁以两个格距为一步跳动。先保留这一子空间内的最低阶运动,不包含额外畴壁的产生、湮灭和基态的虚涨落修正。

然后比较两种情形:

一是分离代价不随距离增加;二是加入一个与间距$r-l$成正比的能量项:

$$ V(r-l)=\kappa(r-l). $$

第二种情况可以表示周围环境偏爱原来的交替排列。中间反向区域越长,与这种偏好不一致的区域越大,两个边界便受到拉回一起的趋势。

本库有效模型中局部相邻翻转后的量子演化。色值是左右两个畴壁位置概率之和,空间总和为2,使用对数色标;每次跨两格的跃迁保留固定键位奇偶性,形成图中竖条。下图比较平均分离距离。两个模型的其他参数相同,时间使用hbar/Jxy单位;概率图不是单次粒子轨迹。
本库有效模型中局部相邻翻转后的量子演化。色值是左右两个畴壁位置概率之和,空间总和为2,使用对数色标;每次跨两格的跃迁保留固定键位奇偶性,形成图中竖条。下图比较平均分离距离。两个模型的其他参数相同,时间使用hbar/Jxy单位;概率图不是单次粒子轨迹。

计算中,初始间距为2个格距。在所展示的12个时间单位内,无额外分离代价时,平均间距增加到约26.0;加入$\kappa=0.15$后,间距先增加再回落,在这个窗口内最大约9.3,末时约5.9。

本库用101点和141点链比较了这个时间窗口,平均间距结果相符;还从小系统中的实际自旋交换矩阵元,独立核对了有效模型的投影形式。概率归一与总自旋投影也保持正确。

这展示的是明确有效模型内的分离与约束,不能拿这些自选时间、距离直接预测真实材料。真实链还可能存在更多激发、不同相互作用和额外散射。

“有约束”也不要求两个激发完全静止。它们可以作为有关联的对象继续运动;区别在于,彼此分离是否要不断支付更高代价。

为什么固定动量下会出现一整片能量连续谱

如果一个探针把总动量和总能量交给两个激发,两者之间可以有很多种分法。同一个总动量,不必对应唯一能量。

对于零温、各向同性自旋1/2 Heisenberg反铁磁链,两自旋子部分的典型上下边界为:

$$ \epsilon_{\mathrm{lower}}(q) =\frac{\pi J}{2}|\sin(qd)|, $$
$$ \epsilon_{\mathrm{upper}}(q) =\pi J\left|\sin\frac{qd}{2}\right|. $$

$J$是交换能量,$d$是相邻链自旋的间距。这些边界可以从单自旋子的色散和两个动量的允许分配关系求出;本库也逐点做了数值极值核对。2

零温各向同性链的两自旋子运动学范围。阴影不是散射强度预测。使用J=0.21 meV作为材料参照;YbAlO3的链间距为晶体c轴晶格常数的一半,所以qd=πL。这个模型条件与前面的强各向异性畴壁演示不同。
零温各向同性链的两自旋子运动学范围。阴影不是散射强度预测。使用J=0.21 meV作为材料参照;YbAlO3的链间距为晶体c轴晶格常数的一半,所以qd=πL。这个模型条件与前面的强各向异性畴壁演示不同。

因此,一片宽的连续谱,并不意味着每个组成它的激发都寿命很短。两个本身比较清晰的激发,也会因为内部动量有很多种分配,组成一个宽谱。

另一方面,普通激发的有限寿命、温度、无序和仪器分辨率也能让谱变宽。只看到“宽”,不足以识别自旋子。更有区分力的是完整的动量—能量结构、边界怎样移动,以及改变温度或相互作用后出现什么新结构。

两自旋子边界也不是完整响应绝对不能越过的围墙。更高数目的自旋子态仍有贡献,其中一部分可以落在这一区域之外;有限温度还会改变可参与跃迁的状态。

真实材料里,为什么单离子很像Ising,交换却可以接近Heisenberg

YbAlO3中的低能自由度来自Yb离子的晶场双重态,可以用有效自旋1/2描述。这是低能模型,不等于把整个离子的自旋轨道结构都消去。

2019年的研究指出,单离子的磁响应具有很强的方向性,但同一链内的交换仍包含量子自旋翻转项,并且接近Heisenberg形式。磁矩对外场的方向性,与交换相互作用的各向异性,不是同一件事。1

这个区别直接影响实验能看见什么。该材料的中子散射主要耦合到纵向自旋分量,对横向分量的灵敏度低得多。前面的相邻翻转示例选择总自旋投影改变为零的制备,帮助理解成对边界;它仍然没有重建真实散射算符的全部矩阵元。

在约0.8至0.9 K以下,较弱的链间作用使材料建立磁有序。邻链的有序对一条链产生交错有效场,使中间反向区域的代价随长度增加,自旋子因此受到约束。

2019年实验在有序转变上下,观察到谱从无可分辨能隙的连续响应,转向带明显能隙和束缚激发特征的响应;高能连续部分仍然可以保留。与交错场模型的比较,把“分离代价改变”与实际谱结构联系起来。1

外场也不能笼统地当作一种“让自旋子更难分开”的旋钮。均匀场会改变磁化和磁相,交错场则在交替位置偏爱相反方向;场的方向、空间形式以及它怎样改变链间有序,都需要分别考虑。

2025年的前沿观察:高温留下的上边界

Kish等人在2025年的研究中,测量了YbAlO3从2 K到100 K的非弹性中子散射,并与有限温度DMRG计算比较。2

这种材料的交换能量约为0.21 meV,对应$J/k_B\approx2.44$ K。因此,40 K已经约为该交换尺度的16倍,100 K约为41倍;这里的“高温”指相对于链内交换能量,并不是高于所有原子能级的尺度。

随着温度上升,下边界逐渐模糊。与此同时,上边界的动量依赖仍然可以辨认。作者用理想链的有限温度谱作为参照,再加入仪器分辨率和额外洛伦兹展宽,讨论环境所限制的有效寿命。

“未分辨出额外展宽”不等于完全没有相互作用或寿命无限。它表示,在该仪器和分析框架下,额外宽度还没有可靠地超过可辨识范围。理想链计算本身已经包含了一部分有限温度多体效应。

原文据此给出了约几十个格距的有效相干长度,并在40 K附近报告超过约35个间距的量级。按$d=c/2\approx0.366$ nm换算,35个间距约为12.8 nm。这个数字来自谱分析和分辨率条件,不是直接用显微镜量到的单个自旋子行程。

作者还给出了一个更简化的谱形分析,得到约13.4个格距的有效长度。两者对多自旋子贡献与外部展宽的划分不同,不能把它们当作同一个无模型尺子的两次读数。

在更高温度,额外展宽开始变得可分辨;其温度依赖与约20 meV的激活能相容。晶场激发参与是作者提出的解释。这个能标远高于交换$J$,有助于理解低能集体动力学为什么能坚持到很高的$T/J$;但相近的激活尺度本身不能替代对具体散射通道的独立验证。

时空图中的传播,和逐颗追踪有什么区别

中子散射记录的是动量和能量相关的响应。把相应谱转换到位置和时间,可以显示一个局部扰动的响应怎样扩展,以及早期较快传播与晚期输运之间的关系。

2025年论文用这种方法重建了时空响应,并与拟合和DMRG结果比较。高温下仍能看到与传播速度有关的结构;作者还讨论了较长时段的超扩散迹象。2

这些图是相关函数和线性响应的重建,并非对两颗自旋子连续拍摄。转换还受到测量窗口、背景、分辨率和数值延拓的影响。

本库前面展示的畴壁概率图,则来自一个选定初态的有效模型演化。两种图都与传播有关,但初态、算符、温度和读数定义不同,不能只因外形相似就当成同一项实验复现。

纠缠见证很小,为什么不能直接说没有纠缠

2025年论文还讨论了量子Fisher信息(QFI)。它与针对某个生成元的状态灵敏度有关;在合适的归一化和局部算符条件下,超过阈值可以证明存在相应的多体纠缠。3

这是一个单向的保证:足够大时能给出认证,低于阈值时则可能只是这项见证没有检测出来。它也不是所有动态相干性质的统一刻度。

用一个只有两个自旋的模型,就能把这个逻辑算清。取反铁磁相互作用$H=J\mathbf S_1\cdot\mathbf S_2$,选择生成元$A=S_1^z-S_2^z$,计算热态的归一化QFI。

在$k_BT/J=0.7$时,本库得到$F_Q/2\approx0.543$,低于见证门槛1。可是,对同一个密度矩阵做部分转置,会得到约$-0.0817$的负本征值。可分离态的部分转置必须非负,因此这里确实有纠缠。

精确两自旋热态例子。阴影区内,这项QFI见证没有越过门槛,但部分转置仍能证明纠缠。数值矩阵计算与闭式表达式相核。这个小模型只用于说明见证的逻辑,不是自旋子模型,也不给YbAlO3的纠缠温度作预测。
精确两自旋热态例子。阴影区内,这项QFI见证没有越过门槛,但部分转置仍能证明纠缠。数值矩阵计算与闭式表达式相核。这个小模型只用于说明见证的逻辑,不是自旋子模型,也不给YbAlO3的纠缠温度作预测。

本例QFI见证在约$k_BT/J=0.477$以上不再越过门槛,而纠缠直到约0.910才消失。两者之间就是不能把“没见证到”写成“没有”的具体区间。

对YbAlO3,原文高温QFI未能提供相应纠缠认证,并不自动否定谱中看到的动态结构;同样,谱中仍有相干特征,也不能补写为已经证明高温下存在任意形式的多体纠缠。

本篇使用标准密度矩阵定义和附录中的归一化来理解见证。原文正文与附录有一处QFI前因子不一致,单激发与双激发的寿命展宽也需要按定义区分;相关问题已在来源记录保留,本库没有据此重新宣称测得该材料的绝对QFI或单粒子寿命。

哪些证据让“集体激发”成为有用的物理对象

一个有效粒子的价值,在于它能把原来复杂的集体运动组织成可计算、可检验的描述:携带哪些量、怎样传播、哪些相互作用会使它束缚或失去清晰的谱特征。

在本题中,成对产生与分离、两激发动量分配形成的连续谱、链间有序改变分离代价,以及高温谱与理论之间的比较,共同支持了自旋子描述的用途。仅凭一个宽峰或一张图,还不足以走完这条推论链。

前沿问题是,这些有效对象在真实材料中能保持到什么温度、受哪些外部自由度限制,以及改变环境后是否还能承担具体功能。谱分析中的相干长度提供了重要线索,实际信息编码、操控与读取的能力还需要另外验证。

可回查:计算结果、畴壁分离随时间变化、选定时刻的概率与自旋投影、两自旋子边界、两自旋热态见证。

来源与依据

  1. S1 原始稀土链结构、磁响应与中子谱/模型比较

    Tomonaga–Luttinger liquid behavior and spinon confinement in YbAlO3

    L. S. Wu, S. E. Nikitin, Z. Wang, W. Zhu, C. D. Batista, A. M. Tsvelik et al. · 2019

    Nature Communications10,698,2019-02-11;完整主文及已取得SI · 相关全文已核

    Results16–37、Fig2/3、有效Hamiltonian及Discussion74–76、Methods79–82

    核单离子Ising样g响应不等于交换Ising;有效自旋在Kramers双重态中,真实离子J=7/2不直接等同S1/2。INS以纵向DeltaSz0为主;0.88K链间有序产生交错有效场和束缚,均匀外场作用不同。低温有能隙/束缚模仍有高能连续响应;未重跑原64点DMRG或拟合全部相图。SI取得,未将其全部段落宣称已独立复核。

  2. S2 原始2–100K中子谱、有限温度DMRG与寿命分析

    High-temperature quantum coherence of spinons in a rare-earth spin chain

    Lazar L. Kish, Andreas Weichselbaum, Daniel M. Pajerowski, Andrei T. Savici, Andrey Podlesnyak, Leonid Vasylechko, Alexei Tsvelik, Robert Konik, Igor A. Zaliznyak · 2025

    Nature Communications16,6594,2025-07-17;完整主文/SI · 相关全文已核

    全主线Results/Discussion/Methods;原Fig1/3/4实际查看;SI归一化S3–S11、卷积/展宽、IIE S40–S45与S7图注

    J.21meV、d=c/2=3.6565Å,qd=πL;零温二自旋子边界非全谱绝对支撑,更高自旋子可在外侧。未分辨额外Γ非无限寿命;Γ还依DMRG已包含的多体效应/分辨率处理,简化谱形13.4格距与主分析约35分解不同。35格距换算12.8nm只复算长度单位,未重拟合Γ。原正文QFI Eq2写1/(4π),SI S10及s3用4/π,前因子不一致;SI S42两激发宽度为单激发2倍,不混单/双宽与FWHM/HWHM。正文关于ξ0的上限/Arrhenius写法不作本库拟合公式。原实空间图是谱的傅里叶响应且乘T/J,不是单粒子录像;原DMRG有64点边界、线性预测/延拓。20meV激活与CEF解释相容但不唯一,未由相干谱补写量子门/任意纠缠能力。

  3. S3 原始QFI与动态响应/纠缠见证理论

    Measuring multipartite entanglement through dynamic susceptibilities

    Philipp Hauke, Markus Heyl, Luca Tagliacozzo, Peter Zoller · 2016

    Nature Physics12,778–782;Strathprints公开作者PDF · 相关全文已核

    前两页:纯态4Var、混态谱表达Eq2、单位局部谱宽的fQ>m条件、动态响应Eq4证明

    实际核所列定义/证明,未复核所有QPT应用。见证是充分判据,未越界不证明可分离。热平衡下QFI与χ''关系有4/π前因子;本库dimer用密度矩阵定义和N=2归一,不重算材料绝对QFI。

最后更新 2026-09-10

来源