低能自旋波在超导电子掺杂铜氧化物中的涌现

背景与学术渊源

起源与学术渊源

本文所解决的问题源于理解高温超导性的长期探索，特别是在一类被称为铜氧化物的材料中。数十年的凝聚态物理研究已经确立了超导性与磁性在这些材料中的关键联系。一个核心的谜团源于观察到电子掺杂（n型）和空穴掺杂（p型）铜氧化物相图的显著不对称性。虽然两类材料都能实现超导，但n型铜氧化物，如Nd$_{1.85}$Ce$_{0.15}$CuO$_{4-\delta}$（NCCO），表现出一种特殊的先决条件：它们在合成后必须经过还原退火处理才能实现超导。未退火的n型铜氧化物保持反铁磁性且不超导。

这种还原退火对材料结构的确切影响及其随后对磁性和超导性质的影响，一直是激烈争论但尚未达成明确共识的主题。先前对退火后的超导NCCO的研究，在其磁涨落谱中识别出了一个“自旋赝隙”。然而，由于未对未退火的、非超导NCCO的低能自旋动力学进行专门研究，理解上仍然存在显著的空白。这种缺乏直接比较的状况，使得关于还原退火如何特异性地影响这些低能自旋动力学以及进而影响超导性的涌现，留下了一个悬而未决的问题。本文旨在通过直接比较退火前后NCCO中的磁激发，来填补这一空白。

直观的领域术语

为了帮助零基础的读者理解核心概念，这里列出了一些论文中的专业术语，并用日常类比进行解释：

• 铜氧化物 (Cuprates): 想象它们是一类特殊的陶瓷材料，如同高级的陶器，它们在冷却到特定温度时，以在没有任何电阻的情况下完美导电（超导性）而闻名。它们之所以特别有趣，是因为与已知其他超导体相比，它们能在相对“高”（尽管仍然非常冷）的温度下实现这一点。
• 反铁磁莫特绝缘体 (Antiferromagnetic Mott Insulators): 想象一个微小的棋盘，每个格子都放着一个微型磁铁（电子的自旋）。在反铁磁材料中，这些磁铁呈交替图案——北-南，北-南——完美对齐，因此它们的整体磁效应相互抵消。而“莫特绝缘体”意味着，尽管具有这些磁性，电子却被固定在原地，无法自由移动以导电，从而使材料成为绝缘体。
• 自旋赝隙 (Spin Pseudogap): 想象一个音乐厅，音乐家（磁激发）通常会演奏从非常低到非常高的完整音域。一个“自旋赝隙”就像是低音区的一个临时“静默区”。它不是完全的寂静，但这些低频声音明显受到抑制或缺失，在可用的音乐范围内造成了一个“间隙”。
• 还原退火 (Reductive Annealing): 想象你烤了一个蛋糕（合成了材料），但它不太对劲。“还原退火”就像是把蛋糕放回一个特殊烤箱，在受控气氛下（可能含氧量较低）进行“再烘烤”，以微妙地改变其内部结构和性质，使其变得完美。对于这些材料来说，这是一个关键的“再烘烤”步骤，用于去除一些氧原子并解锁其超导能力。
• 自旋波 (Spin Waves): 想象一个体育场里的人们在做“波浪式起立”。这里不是人站起来坐下，而是想象材料中电子的微小磁铁（自旋）集体倾斜并通过结构传播。这些“自旋波”是磁能量在材料中传播的方式，研究它们有助于我们理解其磁行为。

符号表

符号	描述	类型
$T_c$	超导性的临界温度	参数
$\hbar\omega$	能量转移（通常与被测激发的能量有关，如自旋波）	变量
$\chi''(\omega)$	动态磁化率（衡量材料磁化强度在振荡磁场作用下改变的难易程度的指标）	变量
$Q$	反铁磁序化波矢量（描述磁序的空间周期性和方向）	参数
$E_{gap}$	自旋赝隙起始能量（自旋赝隙，即低能磁激发被抑制的能量，开始出现的能量）	参数

问题定义与约束

核心问题表述与困境

本文所解决的核心问题是精确定义还原退火在电子掺杂铜氧化物，特别是Nd$_{1.85}$Ce$_{0.15}$CuO$_{4-\delta}$中超导性涌现中的作用。

输入/当前状态 是Nd$_{1.85}$Ce$_{0.15}$CuO$_{4-\delta}$的未退火晶体。在此状态下，该材料是一种非超导反铁磁莫特绝缘体，在其磁涨落谱中表现出大的自旋赝隙。先前的研究并未充分研究这些未退火样品的低能自旋动力学。

输出/目标状态 是在材料从其未退火、非超导状态转变为通过还原退火诱导的超导状态的过程中，全面理解磁激发谱，特别是自旋赝隙，是如何演变的。最终目标是建立材料缺陷、磁性与高温超导性涌现之间的直接、机制性联系。

本文试图弥合的确切缺失环节或数学鸿沟在于未退火、非超导电子掺杂铜氧化物中低能自旋动力学的详细表征，以及与退火后的超导对应物的直接比较。虽然已知退火后的超导NCCO中存在自旋赝隙，但未退火状态下自旋涨落的行为，以及退火如何精确地改变它们以实现超导性，仍然是一个悬而未决的问题。本文旨在量化自旋赝隙的变化，并将其与缺陷的“修复”以及长波自旋波的发展联系起来。

困扰先前研究人员的痛苦的权衡或困境源于铜氧化物中超导性与磁性之间复杂的相互作用，特别是自旋赝隙的作用。普遍的理解是超导性打开了自旋赝隙。然而，本研究揭示了一个显著的困境：未退火的、非超导样品已经表现出比退火后的超导样品（在2 K时起始于2 $\pm$ 0.6 meV）更明显甚至更大的自旋赝隙（在2 K时起始于10 $\pm$ 0.5 meV）。这一观察直接挑战了超导性是打开赝隙唯一原因的普遍观念。相反，它表明诱导超导性的还原退火实际上减小了自旋赝隙。这就产生了一个悖论，即“改进”（超导性）与较小的赝隙相关联，而不是其涌现。

约束与失效模式

电子掺杂铜氧化物中低能自旋波和超导性涌现的问题，由于以下几个严峻的现实约束而变得异常困难：

• 物理与材料约束：

  • n型铜氧化物的化学复杂性： 电子掺杂铜氧化物如NCCO比p型铜氧化物在化学上更复杂，因为它们在合成后需要还原退火处理才能实现超导。这一必要性立即打破了铜氧化物物理中常考虑的电子-空穴对称性。
  • 未退火状态下的缺陷： 未退火的NCCO晶体本身就存在缺陷，包含各种扰动局部离子势并充当散射中心的缺陷（例如，间隙的顶端氧原子、CuO$_{2}$平面中的氧空位或Cu空位）。这些缺陷被认为会抵消超导性。这些缺陷的确切性质和位置，以及退火如何影响它们，一直是长期争论而无共识的问题。
  • 反铁磁基态： 未退火材料以反铁磁基态开始，超导性仅在退火后出现，这暗示了磁序与超导性之间微妙的平衡和竞争。
  • 晶体电场干扰： 非弹性中子散射测量受到Nd晶体电场能级（尤其是在15 meV附近）的强烈干扰，严重受限。这使得在约14 meV以上的能量下可靠测量非弹性信号成为不可能，限制了自旋动力学可观测的能量范围。
  • 样品变异性： 确保未退火和退火状态之间直接比较的样品一致性至关重要。作者通过将一个最优掺杂的晶体分成两部分来解决这个问题，一部分进行退火，另一部分保持未退火状态。这最大限度地减少了样品变异性带来的不确定性，这是材料科学中常见的挑战。
  • 孪晶畴： 还原退火后的超导样品表现出额外的、旋转45°的孪晶畴，有效地减少了用于某些测量的样品质量。这使数据分析复杂化，并可能影响信号质量。

• 计算与数据驱动约束：

  • 微弱磁信号： Cu$^{2+}$ (S = 1/2) 的磁矩很小，导致磁信号微弱，难以检测。这需要使用具有更大分辨率体积的热三轴谱仪来积分这些微弱信号。
  • 归一化挑战： 为了确保不同样品和测量之间的可比性，必须仔细进行强度归一化，通常需要额外的声学声子测量，这增加了实验程序的复杂性。
  • 信号检测的统计严谨性： 确定磁信号是否存在（例如，q扫描中的高斯峰）需要稳健的统计方法，例如使用指定置信区间（p = 0.05）的Wilks定理，以避免未经证实的结论或误解。
  • 绝对单位转换有限： 中子散射实验所需的大晶体通常太大，无法用于标准SQUID磁力计将其磁化强度转换为绝对单位，导致测量结果以每克磁化强度报告，这可能阻碍与某些理论模型的直接定量比较。

为什么选择此方法

选择的必然性

本研究的核心在于揭示电子掺杂铜氧化物中磁性与超导性之间错综复杂的相互作用，特别是聚焦于磁激发——即自旋波——的动态行为以及神秘的自旋赝隙。要真正理解这些现象，我们需要一种能够直接“看到”这些微小磁涨落的工具，而不仅仅是它们的宏观效应。这时，非弹性中子散射，特别是使用热中子三轴谱仪，就不仅仅是一个绝佳的选择，甚至可以说是唯一可行的实验方法。

如果我们从凝聚态物理学的角度考虑，传统的测量方法，如标准磁测量（例如图1所示的SQUID测量），可以告诉我们材料是否超导或反铁磁。然而，它们对磁激发本身的能量和动量是盲目的。它们提供宏观图景，但无法提供微观动力学信息。作者需要探测自旋波如何涌现，其能谱如何随退火而变化，以及缺陷如何影响它们。这需要一种能够同时分辨能量转移（图3-6中的“$\hbar\omega$”轴）和动量转移（图2中的“q”扫描）的技术。中子散射在这方面表现出色，因为中子直接与原子磁矩相互作用，使我们能够绘制出自旋激发谱。没有这种直接探测，自旋赝隙的细微变化以及潜在的自旋动力学，而这些正是本文发现的核心，将完全隐藏起来。

比较优势

除了能够测量自旋动力学之外，非弹性中子散射还提供了定性优势，使其在该特定问题上具有压倒性优势。本文强调了几个结构性优势：

首先，它直接测量自旋波的态密度（如第7页所述并如图7所示）。这意味着它不仅仅是探测一个磁信号，而是提供了一个详细的能量图景，显示了这些自旋激发可能存在的位置。这对于理解自旋赝隙至关重要，自旋赝隙本质上是这些低能自旋态的耗尽。

其次，该技术允许进行精确的动量分辨测量（q扫描，图2）。自旋波是集体激发，在材料中传播，其行为取决于其波长（或动量）。能够绘制出这种“色散”是表征它们的根本。其他技术可能推断出某些方面，但中子散射提供了直接、明确的图像。

第三，作者选择的是热中子三轴谱仪而非冷中子谱仪（第9页）。这是基于信号性质的深思熟虑的决定。热中子因其“更大的分辨率体积”，更适合“积分源于Cu$^{2+}$ (S = 1/2) 小磁矩的微弱信号”。这意味着它们能够有效地捕捉到铜离子发出的微弱磁信号，而铜离子是铜氧化物超导性的关键参与者，即使这些信号难以检测。这一选择直接解决了处理磁矩小、信号可能微弱的材料的实验挑战。

与约束的契合

所选的非弹性中子散射方法与问题的固有约束和要求完美契合，形成了科学问题与实验工具的真正“结合”。

1. 理解自旋波的“涌现”和“演变”： 本文的目标是观察低能自旋波如何涌现以及自旋赝隙如何随还原退火而演变。中子散射直接提供了这些自旋波的能量和动量谱，允许在未退火（非超导）和退火（超导）状态之间进行直接比较（图3、4、5）。这种直接比较，使用来自同一晶体的样品，最大限度地减少了样品变异性，这是可靠结果的关键约束。
2. 聚焦低能动力学： 问题特别针对“低能自旋波”。中子谱仪被调谐到测量2 meV至13 meV范围内的能量转移（第9页），精确捕捉了自旋赝隙表现出来的相关低能区域。
3. 探测缺陷及其影响： 核心假设围绕着退火去除的缺陷如何影响磁性和超导性。通过测量退火之前和之后的自旋动力学，该方法直接揭示了这些缺陷对自旋波谱和赝隙的影响（图7对此进行了精彩说明）。
4. 微弱磁信号： 铜氧化物涉及相对较小的Cu$^{2+}$离子磁矩。选择具有更大分辨率体积的热中子谱仪，是为了响应检测这些“微弱信号”的需求（第9页），确保即使是细微的磁涨落也能被可靠地测量。

替代方案的排除

虽然本文没有明确将讨论框架化为对其他流行机器学习方法的“排除”（因为它是一项基础物理实验），但它间接证明了为什么其他实验技术对于其主要目标来说是不够的。

• SQUID磁测量： 如图1所示，SQUID测量用于确认宏观磁状态（反铁磁 vs. 超导迈斯纳效应）。然而，它们无法提供关于动态自旋激发谱或自旋波能量和动量的信息。对于自旋赝隙的详细理解，SQUID数据是必要的表征，但不是主要的研究工具。
• 角分辨光电子能谱（ARPES）： 本文提到ARPES提供了“对我们中子测量推断出的配对相互作用动量依赖性的额外视角”（第8页）。ARPES探测的是电子激发谱。虽然对于理解电子结构和配对至关重要，但它不能直接测量磁性自旋波及其动力学，而这正是本研究的重点。它是一种互补技术，而非替代品。
• 共振非弹性X射线散射（RIXS）： 一项先前的RIXS研究被引用（参考文献54），以支持作者的模型，但有一个关键区别。该RIXS研究关注的是“高能磁激发谱……能量范围从约100 meV到1 eV”。然而，本文的作者特别关注“低能（即长波长）自旋波”（第8页）。虽然RIXS可以探测磁激发，但其在引用的工作中的应用是在不同于本文核心低能现象的能量尺度上。
• 冷中子三轴谱仪： 这可能是对所选方法的一种替代变体最直接的“排除”。作者明确指出，热中子三轴谱仪比冷中子三轴谱仪更受青睐，“因为其分辨率体积更大，使我们能够积分源于Cu$^{2+}$ (S = 1/2) 小磁矩的微弱信号，以更好地验证磁涨落的存在或不存在”（第9页）。这是选择最适合特定、微弱磁信号检测的中子能量范围的明确技术理由。

因此，选择非弹性中子散射，特别是热中子变体，是一个经过深思熟虑的决定，由直接、高分辨率探测这些复杂材料中低能自旋动力学的需求所决定。

Figure 1. shows the magnetization measurements of the two crystal pieces from the same growth, of which one has been reductively annealed. Note that here the mag- netization is simply given as magnetization per gram of crystal, as the large crystals needed for our neutron scat- tering experiments are generally too large for the SQUID magnetometer to convert into absolute units. The an- nealed sample displays a clear negative magnetization at low temperatures, indicative of the Meissner effect, with an onset temperature of the superconducting transition at Tc = 23 K. In contrast, the as-grown sample shows a flat magnetization curve, with only a slight increase at low temperatures. This is typical of an antiferromag- netic response and clearly differs from the sharp super- conducting transition. The insert shows the tetragonal crystal structure, I4/mmm for both annealed and as- grown, optimally doped NCCO with lattice parameters a = b = 3.957 ˚A and c = 12.075 ˚A.25

数学与逻辑机制

主方程

坦白说，本文主要是一项使用中子谱仪的凝聚态物理实验研究，它并没有呈现一个单一的“主方程”，例如驱动底层物理现象或计算模型的客观函数、常微分方程/随机微分方程或复杂的变换逻辑。相反，作者精心测量和分析的核心数学可观测量是动态磁化率 $\chi''(\omega)$。该量对于理解磁激发和自旋赝隙至关重要。

文章指出，“积分强度被转换为动态磁化率 $\chi''(\omega)$”（第4页），并且该量随后用于通过使用误差函数进行拟合来确定自旋赝隙起始能量（第4页，参考补充说明4）。虽然 $\chi''(\omega)$ 的完整理论定义或用于拟合的误差函数的精确形式并未在正文中明确列出，但 $\chi''(\omega)$ 是被研究行为的基本量。

在非弹性中子散射的背景下，动态磁化率 $\chi''(\mathbf{Q}, \omega)$ 与动态结构因子 $S(\mathbf{Q}, \omega)$ 相关，后者与测量的中子散射截面成正比。作者对动量转移 $\mathbf{Q}$ 进行积分以获得 $\chi''(\omega)$。动态磁化率的虚部作为响应函数的一般形式可以认为是：

$$ \chi''(\omega) = \int d\mathbf{Q} \, \chi''(\mathbf{Q}, \omega) $$

其中 $\chi''(\mathbf{Q}, \omega)$ 是动量和能量分辨的动态磁化率。该积分代表了在给定能量 $\omega$ 下，所有相关动量转移 $\mathbf{Q}$ 的总磁耗散响应。

逐项解剖

让我们来剖析核心可观测量 $\chi''(\omega)$ 以及与其测量和解释相关的概念。

• $\chi''$（动态磁化率，虚部）：

  1. 数学定义： $\chi''$ 是复数动态磁化率 $\chi(\mathbf{Q}, \omega) = \chi'(\mathbf{Q}, \omega) + i\chi''(\mathbf{Q}, \omega)$ 的虚部。它描述了材料对振荡磁场的磁响应的耗散部分。在中子散射的背景下，根据涨落-耗散定理，尤其是在低温和低能情况下，它与动态结构因子 $S(\mathbf{Q}, \omega)$ 成正比。作者对动量 $\mathbf{Q}$ 进行积分以获得 $\chi''(\omega)$。
  2. 物理/逻辑作用： 该项是材料中磁激发，特别是自旋波的直接探针。非零的 $\chi''(\omega)$ 表明材料可以通过激发自旋波从入射中子那里吸收能量。其幅度反映了这些激发的密度和强度。在低能量下 $\chi''(\omega)$ 的抑制表明存在“自旋赝隙”，意味着没有（或很少有）低能磁激发可用。
  3. 为何积分（隐含）： 作者对动量分辨的磁化率 $\chi''(\mathbf{Q}, \omega)$ 在 $\mathbf{Q}$ 上进行积分以获得 $\chi''(\omega)$。该积分的使用是因为自旋赝隙是在反铁磁序化波矢量周围一系列动量转移上观察到的现象。对 $\mathbf{Q}$ 的积分将多维数据压缩为能量 $\omega$ 下的动态磁化率的动量平均图景，有效地给出了自旋激发的态密度，这对于表征赝隙至关重要。积分而非求和的选择反映了块体材料中动量空间的连续性。

• $\omega$（能量转移）：

  1. 数学定义： $\omega$ 代表激发的角频率，通过 $\Delta E = \hbar\omega$ 与能量转移相关，其中 $\hbar$ 是约化普朗克常数。
  2. 物理/逻辑作用： 该变量决定了被探测的磁激发的能量尺度。通过改变 $\omega$（或 $\hbar\omega$），实验者可以绘制出自旋波的能量谱。在特定 $\omega$ 值下信号的存在或不存在直接揭示了磁涨落的能量图景，从而可以识别间隙或峰值。
  3. 为何是连续变量： 非弹性散射中的能量转移是一个连续变量，反映了固体中可能的激发的连续谱。因此，$\chi''$ 自然地表示为 $\omega$ 的函数，允许进行详细的光谱分析，而不是离散点。

逐步流程

想象一个单一的抽象中子，代表数据流通过这个实验和分析管道：

1. 初始状态： 我们的抽象中子，携带特定的初始能量（$E_i$）和动量（$\mathbf{k}_i$），接近 Nd$_{1.85}$Ce$_{0.15}$CuO$_{4-\delta}$ 样品。
2. 相互作用与散射： 中子与样品内的磁矩发生相互作用。如果存在一个能量为 $\hbar\omega$ 和动量为 $\mathbf{Q}$ 的磁激发（如自旋波），中子可以产生或湮灭该激发。然后，我们的中子散射出去，以新的最终能量（$E_f$）和动量（$\mathbf{k}_f$）出现。
3. 能量和动量转移计算： 探测器测量 $E_f$ 和 $\mathbf{k}_f$。由此，精确确定了这次单次散射事件的能量转移 $\hbar\omega = E_i - E_f$ 和动量转移 $\mathbf{Q} = \mathbf{k}_i - \mathbf{k}_f$。
4. 原始数据累积： 这个过程重复了数百万次。特定范围 $\mathbf{Q}$ 和 $\omega$ 内散射中子的数量（计数）被累积起来，形成一个原始强度图，例如图2所示的q扫描。该图显示了磁激发在动量-能量空间中的位置。
5. 归一化： 在进行有意义的比较之前，这些原始计数需要进行归一化。这包括除以解释仪器效率、样品体积和其他实验条件的因子，通常是通过与已知标准（如声学声子信号）进行比较（第9页）。这确保了测量的强度真正反映了样品的内在磁特性，而不是实验伪影。
6. 动量积分： 归一化后的强度（与 $\chi''(\mathbf{Q}, \omega)$ 成正比）然后在相关范围的动量转移 $\mathbf{Q}$ 上进行积分（例如，围绕反铁磁序化矢量）。该积分将多维数据压缩为能量 $\omega$ 下的动态磁化率的一维谱，代表了每个能量 $\omega$ 下的总磁响应。这是图3和图5所示的数据。
7. 自旋赝隙确定： 最后，对 $\chi''(\omega)$ 谱进行分析。作者使用误差函数拟合该曲线。从该拟合中提取参数，特别是起始能量（$E_{gap}$）。如果 $\chi''(\omega)$ 在低能量下显著抑制，则识别出自旋赝隙，并量化其能量尺度（$E_{gap}$）。整个过程对未退火和退火样品在不同温度下重复进行，从而可以对它们在不同条件下的磁激发谱进行直接比较。

优化动力学

本文中的“优化动力学”可以从两个互补的角度来理解：材料的物理转变和数据拟合的分析过程。

1. 物理系统的“更新”（还原退火）：
   本研究中的主要“更新”机制是应用于Nd$_{1.85}$Ce$_{0.15}$CuO$_{4-\delta}$样品的还原退火过程。这不是一种算法学习过程，而是一种改变材料状态的物理转变。

   • 初始状态（未退火）： 未退火样品具有高浓度的缺陷（例如，间隙氧原子或Cu空位），这些缺陷充当自旋波的散射中心。这导致反铁磁斑块破碎，并产生大的自旋赝隙，抑制低能自旋涨落。
   • “学习”/“更新”机制： 还原退火涉及在还原气氛中加热样品。该过程通过去除过量的氧原子或允许Cu原子迁移并填充空位来“修复”缺陷。这类似于系统“学习”变得更有序、电阻更低。
   • “损失景观”（比喻）： 可以将材料的能量景观看作是比喻。未退火状态，连同其缺陷，可能代表一个能量较高、稳定性较低的构型。退火驱动系统趋向于一个能量较低、更有序、超导的状态。这里的“梯度”是趋向平衡的热力学驱动力。
   • 收敛状态（退火/超导）： 退火后的样品表现出较少的缺陷，允许形成更大的反铁磁斑块和长波自旋波。这导致自旋赝隙显著减小，并出现超导性。材料已“收敛”到一个对超导性更优化的物理状态。

2. 数据分析“收敛”（曲线拟合）：
   在数据分析中，采用标准的曲线拟合技术从测量的 $\chi''(\omega)$ 谱中提取定量信息。

   • 目标： 目标是找到最佳拟合参数（例如，对于高斯函数或误差函数，$E_{gap}$、峰值幅度、宽度），这些参数能够描述实验数据。
   • 损失函数： 这种拟合的常见损失函数是实验数据点与所选拟合模型（例如，高斯函数用于峰值，误差函数用于赝隙起始）之间残差平方和（SSR）。目标是最小化此损失。
   • 梯度行为： 优化算法（例如，Levenberg-Marquardt、梯度下降）迭代地调整拟合参数。“梯度”函数相对于每个参数的梯度指示了向最小值的最陡下降方向。
   • 损失景观： 参数空间形成一个“损失景观”，其中损失函数的值会变化。算法在该景观中导航，寻找全局最小值，这对应于最佳拟合参数集。
   • 收敛： 当参数和损失函数的变化低于预定容差时，拟合过程“收敛”，表明已找到局部（理想情况下是全局）最小值。这提供了自旋赝隙（$E_{gap}$）和其他光谱特征的定量度量，从而可以对不同样品和条件进行精确比较。作者使用统计检验，如Wilks定理，来评估其拟合的显著性并确保结论的稳健性。

结果、局限性与结论

实验设计与基线

为了严格研究还原退火对电子掺杂铜氧化物中磁激发和超导性的影响，作者采用了一种以直接比较为中心、细致的实验设计。他们方法的核心是使用一个单一的最优掺杂Nd$_{1.85}$Ce$_{0.15}$CuO$_{4-\delta}$晶体，然后将其分成两部分。其中一半经过了诱导超导性的还原退火处理，而另一半则保持其未退火的非超导状态。这一策略对于最大限度地减少样品间的变异性至关重要，确保任何观察到的差异都可以直接归因于退火过程。

本研究中的“受害者”（基线模型）本质上是该材料本身的未退火、非超导状态。退火效应的第一个明确证据来自磁测量。使用Quantum Design MPMS-XL SQUID磁力计，在1.8 K至50 K的温度范围内，以10 Oe的施加场进行了零场冷（ZFC）测量。退火样品在低温下明显显示出负磁化强度，这是迈斯纳效应的标志，超导转变起始温度为$T_c = 23$ K（图1）。

相比之下，未退火样品表现出平坦的磁化强度曲线，仅在低温下略有增加，这符合反铁磁响应的特征，并证实了其非超导性质。这确立了两种样品之间的基本差异。

探测磁激发的首要实验技术是非弹性中子散射，在ANSTO TAIPAN仪器和ILL IN20仪器上进行。选择这些热中子三轴谱仪是因为它们具有更大的分辨率体积，这有利于检测Cu$^{2+}$ ($S = 1/2$) 小磁矩产生的微弱信号。通过X射线和中子劳厄衍射的组合，对两个晶体样品在(h, k, 0)平面进行了仔细对齐。实验包括在反铁磁布拉格点(h, 1-h, 0)（对于 $h = 0.5$）周围进行对角q扫描，能量转移（$\hbar\omega$）范围从2 meV到13 meV。在两个关键温度下进行测量：1.9 K（基准温度，远低于退火样品的 $T_c$）和27 K（高于退火样品的 $T_c$）。为了追踪自旋赝隙的温度依赖性，在固定能量转移 $\hbar\omega = 2$ meV 和 $\hbar\omega = 8$ meV 下，从2 K到55 K测量了磁信号。

实验设计的一个关键方面是中子散射强度的归一化。为了确保未退火和退火样品之间的直接可比性，强度被归一化到声学声子扫描（例如，在(2,0,0)处）。自旋赝隙起始能量（$E_{gap}$）通过使用误差函数拟合动态磁化率 $\chi''(\omega)$ 来确定。对于2 K时的未退火样品，由于未观察到明显的饱和，$E_{gap}$估计为2 K和27 K数据点重叠的第一个数据点。这种严谨的方法允许对不同条件下的磁激发谱进行直接、定量的比较。

证据证明了什么

本文提供的证据有力地叙述了电子掺杂铜氧化物中缺陷、磁性与超导性之间的相互作用，直接挑战了一些传统理解。最初的磁测量（图1）明确表明，还原退火成功地将未退火的反铁磁Nd$_{1.85}$Ce$_{0.15}$CuO$_{4-\delta}$晶体转变为具有23 K $T_c$ 的超导体，如迈斯纳效应所证实的。这为后续的磁学研究确立了该材料的两种不同状态。

核心机制通过对同一晶体退火前后自旋激发谱的直接比较得到了无情地证明。对于未退火的、非超导样品，中子散射数据显示在27 K时 $\hbar\omega = 6$ meV 处存在清晰的磁响应峰（图2）。然而，在2 K时，该信号显著消失，表明低能磁激发受到强烈抑制。更引人注目的是，2 K时未退火样品的动态磁化率 $\chi''(\omega)$ 显示出明显的自旋赝隙，从约10 meV逐渐演变到4 meV，起始点为10 $\pm$ 0.5 meV（图3a）。

温度依赖性进一步凸显了这一点：在低于约40 K的所有温度下，8 meV处的能量涨落比2 meV处的更强（图5a，图6），证实了未退火状态下低能涨落的抑制。

相比之下，退火后的超导样品表现出截然不同的行为。在2 K时（处于超导状态），它仅显示出2 $\pm$ 0.6 meV 的小自旋赝隙（图3b）。在 $T_c$ 以上，未观察到自旋赝隙。对2 K和27 K之间动态磁化率变化量（$\Delta\chi''$）的分析（图4）进一步强调了这种差异：退火样品显示出自旋赝隙急剧关闭，起始点为3.0 $\pm$ 0.1 meV。2 meV的涨落一直占主导地位，直到它们在约5 K以下被间隙化（图5b，图6），这与未退火样品的行为明显不同。

这是明确的、不容否认的证据：诱导超导性的还原退火，同时减小了自旋赝隙，从未退火状态的大约10 meV减小到超导状态的约3 meV。这直接挑战了超导性打开自旋赝隙的普遍观点。相反，作者提出，未退火样品中的大赝隙是由缺陷引起的，这些缺陷使CuO$_{2}$平面破碎，从而抑制了长波自旋波。退火“修复”了这些缺陷，允许形成更长波长的自旋波并占据更低的能量态，这反过来减小了自旋赝隙（图7）。

这提供了一个直接、来之不易的联系，将材料缺陷、磁关联和高温超导性涌现联系起来。此外，弹性散射测量证实，与未退火样品相比，超导样品表现出反铁磁序被抑制，这与竞争序的广泛理解一致。

局限性与未来方向

尽管这项研究为缺陷和还原退火在电子掺杂铜氧化物中形成自旋赝隙和实现超导性中的作用提供了令人信服的证据，但它也突显了几项局限性，并为未来的研究开辟了丰富的途径。

一个显著的局限性是由于Nd晶体电场能级（约15 meV）的干扰，无法测量超过14 meV的非弹性信号。这限制了对更高能量磁激发谱的全面理解，而这些谱可能包含关于配对机制或其他磁现象的进一步线索。

也许最关键的局限性是作者自己承认的，即关于还原退火对材料结构的确切化学后果的“长期争论”和缺乏共识。虽然本文提出了一个缺陷修复模型，但它指出文献中存在相互矛盾的假设，关于退火是主要减少顶端氧缺陷、产生平面氧空位，还是“修复”Cu位点。在没有对确切的原子级变化有明确理解的情况下，尽管强烈暗示，但“缺陷修复”与观察到的自旋动力学之间的直接联系仍然在某种程度上是概念性的。

另一个讨论点是，与一些p型铜氧化物甚至相关的n型铜氧化物（如PLCCO）中观察到的情况不同，该超导样品中缺乏清晰的共振峰。这引发了关于共振峰作为超导性标志的普遍性问题，并表明Nd$_{1.85}$Ce$_{0.15}$CuO$_{4-\delta}$中的磁激发谱可能具有独特的特征。

此外，本文指出，先前关于退火NCCO中磁关联长度（$\xi$）的研究发现它比未退火NCCO小，这似乎与暗示更大、更少破碎斑块的“晶格修复”模型相矛盾。尽管作者将这种差异归因于散射方法和积分能量范围的不同，但这突显了跨不同实验技术比较结果的复杂性，以及对更统一理论框架的需求。

展望未来，出现了几个可以进一步发展和演变这些发现的讨论主题：

1. 阐明退火的原子尺度机制： 未来的研究应侧重于明确解决还原退火引起的化学和结构变化。先进的实验技术，如原子分辨率扫描透射电子显微镜（STEM）、X射线吸收光谱（XAS）或核磁共振（NMR），可以提供氧空位产生、顶端氧去除或Cu位点重构的直接证据。这将为“缺陷修复”模型提供坚实的基础。

2. 工程化缺陷景观： 如果缺陷确实会破碎自旋链并抑制低能自旋波，我们能否有意地工程化特定的缺陷类型和密度来调节磁性能并可能增强超导性？这可能涉及受控掺杂策略或超越简单退火的合成后处理。

3. 重新审视自旋涨落介导的配对： 研究结果支持自旋涨落介导的配对机制，退火后自旋赝隙的减小与超导性相关。这挑战了“更大的自旋赝隙总是更有益”的观点。需要进一步的理论建模，可能结合特定的缺陷结构及其对自旋波色散的影响，才能完全理解这些改变的自旋涨落如何促成配对。与角分辨光电子能谱（ARPES）数据的更深入比较，正如作者所建议的，也可以为配对相互作用的动量依赖性提供关键见解。

4. 磁信号的普遍性： 在该NCCO样品中缺乏清晰的共振峰，值得进一步研究。这是材料特有的特征，还是取决于掺杂、测量条件或自旋涨落的确切性质？跨越更广泛的电子掺杂和空穴掺杂铜氧化物使用一致方法学的比较研究，可以帮助确定哪些磁特征真正普遍适用于高温超导性。

5. 关联长度差异的弥合： 一项在相同样品（未退火 vs. 退火）上使用能量积分和非弹性中子散射方法直接比较磁关联长度的研究，可以帮助调和文献中明显的矛盾。这将阐明缺陷如何影响静态磁序和动态自旋涨落。

通过解决这些问题，我们可以更全面、更细致地理解结构缺陷、磁激发与高温超导性涌现之间错综复杂的关系，最终指导改进超导材料的设计。

Figure 1. Magnetization measurements as a function of temperature. Zero-field cooled (ZFC) measurement at 10 Oe applied field, for the as-grown and reductively annealed, superconducting Nd1.85Ce0.15CuO4–δsingle crystals, depicted in blue pen- tagons and orange triangles, respectively. The criti- cal temperature Tc is defined as the onset tempera- ture of superconductivity. Insert: crystal structure of Nd1.85Ce0.15CuO4–δ,24,25 with Cu, O and Nd depicted in blue, red and green, respectively. The 15% Ce doping on the Nd site is denoted as a pink slice on the green Nd atoms Figure 3. Dynamic susceptibility χ′′(ω), as a function of energy transfer. a) as-grown sample. b) annealed, superconducting sample. The black out- lined points indicate 3-point scans, while colored out- lined points indicate q-scans. Error bars represent the fitting error of the area under the Gaussian signal. For the 3-point scans, error bars are determined as outlined in Supplementary Note 2. The solid lines are fits to the response following Supplementary Note 4. The dashed lines are drawn as guide to the eye, while the colored ver- tical dotted lines are the estimate of the spin pseudogap onset with the faded area representing the uncertainty Figure 7. Schematic illustrating how the size of the antiferromagnetic patches influences the spin waves allowed in the system. Left column (a-d) rep- resents the as-grown sample, the right column (e-h) rep- resents the annealed, superconducting sample. a) and e) show the antiferromagnetic structure in each case, with the structure composed of smaller patches created by defects (black circles), that are still weakly antiferro- magnetically interacting. In the annealed sample, the undisturbed patches are larger. b) and f) show how the patches restrict the spin waves above a certain wave- length. By having larger patches, more low-energy states are occupied, minimizing the energy spin pseudogap, as illustrated in c) and g). This is more quantitatively ex- pressed as a (partial) suppression of the spin wave density of states (DoS) at low energies, seen in d) and h)

与其他领域的同构性

结构骨架

本文阐述的核心机制是：结构缺陷（缺陷）如何破坏有序系统，抑制低能集体激发，以及它们的去除如何恢复这些激发，从而产生新的宏观性质。