一条自旋链中的三种激发区间
这篇论文要回答的问题是:在S = 1/2海森堡反铁磁链化合物铜苯甲酸盐的同一块高质量单晶中,能否同时存在几种不同的低能激发,核磁共振(NMR)又能否把它们区分开。简短的回答是可以:作者报告了0.2 T下的自旋子特征、0.4 K以下一个额外的无能隙弛豫通道(他们将其归因于类随机单态行为),以及1 T以上由磁场诱导的自旋能隙(他们将其认定为孤子激发)。
这种区分之所以重要,是因为这些区间在温度上重叠,但在磁场上并不重叠。自旋子的证据来自0.2 T下1.4至0.03 K之间几乎不随温度变化的自旋-晶格弛豫率1/T1,以及约2 K以下按T^(-1/2)上升的自旋-自旋弛豫率1/T2G。类随机单态通道表现为0.4 K以下出现的第二个慢得多的弛豫分量。孤子区间则表现为2至10 T下约2 K以下奈特位移的上翘,并用热激活能隙进行拟合。
材料的选择本身就是论证的一部分。铜苯甲酸盐即使发生有序,其有序温度也在0.8 mK或更低,因此在30 mK以上进行的实验应当看到一条近乎纯粹的一维链。正是这种一维性,使得本征的链激发和由无序驱动的低温行为可以被当作彼此独立的贡献来讨论,而不是混成一个模糊的响应。
平坦的弛豫率为何指向自旋子
1/T1衡量核自旋与周围电子自旋交换能量的快慢,因此它对低能磁激发的密度和性质敏感。在自旋子图像中,链的spin-1/2激发具有线性色散,标度关系1/T1正比于T^(eta-1)在eta = 1时给出不随温度变化的弛豫率。因此,0.2 T下观察到的平台是自旋子的核心证据,而不只是弛豫率变慢。
第二个自旋子指标是1/T2G,它反映邻近自旋产生的局域场分布。它在2 K以下拟合出的T^(-1/2)上翘与自旋子的预期行为一致,也与此前在Sr2CuO3中报道的结果一致。作者还指出,在4 K以上,1/T1遵循随T线性变化的标度关系,与零波矢处的自旋涨落一致,这属于与低温平台不同的区间。
一个预期中的特征并未出现。低温下1/T1的预期对数修正没有被观察到;相反,1/T1f在0.3 K以下出现一个微弱的下弯,并在0.5 T和1 T下增强。作者把这个下弯归因于磁场诱导效应,而不是纯链本身,因此自旋子的结论专门绑定在最低磁场的数据上。
两个弛豫通道把无序与纯链分开
0.2 K以下,核磁化恢复曲线不再能用单指数拟合。作者采用包含快、慢两个T1分量的双指数形式,并检验了这种双分量行为是本征的,而不是谱线重叠造成的假象。随后,逆拉普拉斯变换分析给出每个温度下的弛豫率分布P(1/T1),使两个通道以彼此分离的分布显现出来。
快分量在受到磁场诱导抑制之前大体不随温度变化,与自旋子通道相符。慢分量约小两个数量级,在0.4 K以下出现,其速率随降温线性减小。线性的1/T1属于无能隙行为,纯海森堡链模型不会产生这种行为,因此作者把它归因于受淬火无序影响的区域。
慢分量的相对谱权重w2在0.03至0.4 K范围内按-ln(T/T0)对数增长,拟合得到T0 = 0.40(2) K、c = 25(1)。在0.03 K时,w2达到65%并继续上升,作者将其解读为向另一种基态的渡越,而不是在某个确定温度发生的相变。
发散的奈特位移支持类随机单态解读
奈特位移Kn反映质子感受到的局域磁场,因此它追踪链的均匀自旋磁化率。在0.2 T下,Kn被经验性地拟合为Kn正比于T^(-1) ln^2(Omega/T),其中Omega约为10(1) K,接近约18 K的交换能标J。作者以此作为向类随机单态行为渡越的证据,并强调该函数形式不同于常规杂质情形——后者通常表现为有能隙。
他们还从低频H(1)谱分量中提取w2,并与弛豫分析得到的w2比较。两者都随降温上升,也都遵循对数形式,但拟合出的T0值不同。作者用这一差异论证T0是一个较宽的渡越温度,而不是相变温度。
这一结果附带两点提醒。Kn的拟合只在一个很窄的温度窗口内进行,作者表示其确切函数形式仍需进一步验证。他们还将0.2 K以下Kn的下弯视为谱线展宽和信噪比下降造成的可能假象,而不是物理特征。
磁场诱导的能隙指向孤子
在2至10 T的磁场下,Kn在约2 K以下上翘,并在最低温度饱和。作者用热激活形式Kn = a exp(-Delta/kBT) + b拟合,其中Delta是激发能隙,b来自交错场造成的零温自旋倾斜效应。由这些拟合得到的能隙在1至10 T范围内随磁场单调增大。
磁场的依赖关系是关键判别证据。对NMR数据的幂律拟合给出Delta正比于J(h/J)^0.8(1),其中h是有效交错场;而孤子激发的理论曲线是Delta正比于J(h/J)^(2/3)。作者把自己的能隙值与比热及非弹性中子散射的结果比较,并将整体行为描述为与交错磁化诱导的孤子激发一致。
同样的带隙图像也出现在P5位点的1/T1中,那里的数据用带隙激发谱拟合。作者把低能孤子与可积的sine-Gordon描述联系起来,并指出该模型的耦合常数随磁场变化,他们认为这为研究其他理论预言的激发(如呼吸子)提供了机会。
类随机单态背后的无序有多强
论文没有对晶体进行有意掺杂,因此作者认为无序很可能是畴类型的:由大片更干净的区域隔开的小块更脏区域,可能来自生长或冷却过程中的应力,或苯甲酸基团转动冻结的不平衡。他们用三个独立的能量尺度估计无序水平:类随机单态体积分数在T0约0.4 K以下开始非零、0.8 K以下展宽上翘,以及无能隙行为在0.5 T以下磁场中保持稳定(对应约0.4 K)。
结合约18 K的J,这些尺度指向约0.02至0.05的Delta J/J无序强度。作者把这一弱无序水平与BaCu2(Si0.5Ge0.5)2O7中报道的更强Ge掺杂无序作对比。他们还对双峰耦合分布的无序链做了量子蒙特卡洛模拟,其中10%的键带有Delta J = 0.03,用来模拟这种弱无序链的关联长度。
作者用随机单态图像解释该材料此前的一个观测:0.3 K以下热导率的抑制。按他们的解释,随机分布的单态内部的局域自旋激发会阻碍巡游自旋子,而在严格一维链中,渗流阈值恰好位于零无序处,因此即使很弱的无序也能中断输运。
随机单态图像适用的边界
随机单态情形定义在零温和零场极限,而这里的关键证据却是在0.2 T下采集的。作者直接回应了这一矛盾:磁场会使长距离单态极化,但也会降低交换耦合与塞曼能量相当的配对的三重态能隙,使Delta接近g muB H的配对变得有效无能隙。他们认为,这一组合在低场下保留了有限的低能态密度。
作者指出,在0.5 T以上,塞曼效应会部分破坏类随机单态,通过相分离造成带隙与无能隙激发的共存。这就是类随机单态的结论锚定在0.2 T数据、而不延伸到用于孤子分析的更高磁场的原因。
作者还标出了尚未测量的内容。他们呼吁未来直接测量关联长度及其确切函数形式,并且不能完全排除稀有缺陷附近链间耦合的局域改变。因此,类随机单态的认定依赖弛豫、谱权重和奈特位移三方面特征的组合,而不是对单态配对的直接测量。
区分各区间的是磁场轴而非单一温度
阅读这篇论文最有用的方式,是把它当作一张随磁场变化的图谱,而不是关于铜苯甲酸盐的单一结论。在0.2 T下,恒定的1/T1平台和1/T2G的T^(-1/2)上升标识出自旋子。在同一低场的0.4 K以下,1/T1正比于T的慢弛豫通道和对数增长的w2标识出类随机单态贡献,其在0.03 K时达到65%的谱权重。在1 T以上,激活式的奈特位移上翘和大致按H^(2/3)增长的能隙标识出孤子。
真正重要的比较是实测能隙指数与理论指数之间:NMR拟合给出0.8(1),而理论预言为2/3,作者把孤子的认定表述为与该预言一致,而非精确证实。同样,类随机单态由三个特征支持,但被描述为有限温度下的渡越,T0被视为较宽的渡越尺度而非转变温度。
对于想判断这些结果能延伸多远的读者,论文自身的边界就是磁场轴:类随机单态行为被论证在0.5 T以下仍然存在,孤子区间在1至10 T得到确立,自旋子平台在0.2 T最干净。其独特贡献在于:同一种材料在毫开尔文温度下经NMR局域探测,同时展现出这三种区间,并使无序得以与链的本征响应分开。
本文引用的文献
铜苯甲酸盐自旋链化合物中的自旋子、孤子和类随机单态行为
在低至0.2 T的磁场和低至0.03 K的温度下对铜苯甲酸盐进行质子核磁共振,分辨出1/T1中的自旋子平台、0.4 K以下谱权重对数增长的慢速无能隙弛豫通道,以及1 T以上大致按H^(2/3)标度的磁场诱导能隙。
