Theory of frozen flux in a narrow uniform superconducting strip after cooling in a small magnetic field

本文通过推导并求解一个动态平衡方程，以量化热激活和冷却速率如何决定在小磁场下穿过转变温度时，窄超导带中的冻结温度及由此产生的剩余涡旋密度。

要点

想象一下，一条超导带就像一条长而窄的走廊。在走廊内部，一些被称为“涡旋”的微小磁粒子想要居住其中。走廊的边缘（带子的边缘）以及一种被称为“迈斯纳效应”的特殊力量，创造了一个凹凸不平的能量景观。你可以把这个景观想象成一系列的山丘和山谷。

当带子处于高温时，这些涡旋充满能量且躁动不安。它们可以轻易地翻越“山丘”（能量势垒）进入走廊，或者从走廊中逃逸出来。随着带子冷却，涡旋失去能量。最终，这些山丘变得太高，令它们无法攀爬，于是它们就被困住了。

这篇由 Alexei E. Kosselev 撰写的论文，研究了在磁场中随温度降低，这些涡旋究竟在何时以及如何被困住（或“冻结”）的。以下是使用日常类比对研究结果进行的解析：

1. 设置：一条狭窄的走廊

该研究关注的是极薄、极窄的超导材料带。在这些窄带中，物理机制比宽带要简单得多。阻挡涡旋进入的“山丘”是由带子本身的几何形状产生的。

• 最小排斥场 ($H_e$)： 想象一种如此微弱的磁场强度，以至于“山丘”变得极高，导致没有任何涡旋能够进入。这是理论上的极限，此时带子是完全洁净的。
• 现实情况： 在真实的实验中，即使磁场强度高于这个理论极限，科学家们经常仍会观察到被捕获的涡旋。这篇论文探讨的问题是：为什么？

2. 与时间的赛跑：冷却过程

问题的关键在于冷却。

• 平衡态： 如果你能以无限慢的速度冷却带子，涡旋将有充足的时间找到完美的平衡。如果磁场太强，它们会离开走廊；如果磁场恰到好处，它们则会留下。
• 冻结现象 (The Freeze-Out)： 在现实世界中，我们以特定的速度冷却物体。随着温度下降，“山丘”变得越来越陡峭，涡旋也变得越来越迟缓。在某个临界点，涡旋变得如此笨拙，以至于即便“理想状态”显示它们应该离开，它们也无法在足够快的时间内爬过山丘逃脱。
• 冻结温度 ($T_{fr}$)： 这是指涡旋停止逃离并被困住的具体时刻（温度）。论文精确计算了这种情况何时发生。

3. “冻结”机制

作者描述了一种“动态平衡”。把它想象成走廊里一个繁忙的门口：

• 进入： 涡旋试图跳进来。
• 退出： 涡旋试图跳出去。
• 平衡： 在高温下，人们（涡旋）在进出之间快速往返。走廊内部的人数根据外部的拥挤程度保持稳定。
• 锁定： 随着温度下降，“出口门”变得极其难以开启。内部的涡旋无法出去。同时，“入口门”也变得难以开启，但已经在里面的涡旋现在被困住了。
• 结果： 被捕获的涡旋数量停止变化，并保持在一个固定的数值，尽管“理想”的数量应该是零。这就是“冻结磁通量”。

4. 用通俗语言解释的核心发现

• 它发生在非常接近“熔化”点的位置： 涡旋并不是在带子很冷的时候冻结的；它们是在带子刚刚开始进入超导状态时（非常接近转变温度）发生冻结的。
• “对数”因子： 论文发现，冻结发生的温度比随机热噪声通常起作用的温度略高。这是一个微小的差异，但在数学上具有显著意义（被称为“大的对数因子”）。
• 速度至关重要： 如果你更慢地冷却带子，涡旋就有更多时间逃脱，因此它们会在更低的温度下冻结，捕获的涡旋也会更少。如果你冷却得更快，它们会更早被捕获，留下的涡旋也会更多。
• 磁场是一个开关： 捕获的磁通量量取决于磁场强度的变化。
  • 仅略高于最小极限（$H_e$）时，捕获的涡旋数量极少（几乎为零）。
  • 当你稍微增加磁场，捕获的涡旋数量就会爆炸式增长（急剧增加）。
  • 由于这种剧烈的增长，科学家可以定义一个“有效排斥场”。这是指被捕获的涡旋强度足以被仪器检测到的磁场强度。

5. 为什么现实实验与理论不同

论文解释了一个常见的谜题：实验通常显示，需要比简单数学预测的强得多的磁场，才能使带子达到“洁净”状态（无涡旋）。

• 解释： 数学模型假设走廊是完美光滑、均匀的。而真实的带子存在凸起、划痕和杂质（不均匀性）。
• 影响： 这些缺陷可以充当“陷阱”，即使在低磁场下也能将涡旋固定在原处。这使得看起来像是带子捕获了比预期更多的磁通量，从而将“有效”排斥场推向了更高的数值。

总结

这篇论文提供了一个数学“配方”，用于预测在冷却过程中，会有多少磁涡旋被困在狭窄的超导带中。它解释了涡旋被捕获并非因为磁场太强，而是因为带子冷却得太快，导致涡旋无法及时逃离能量势垒。这种“冻结”现象发生在材料进入超导态的温度附近，且捕获的磁通量量取决于冷却速度和精确的磁场强度。

技术摘要

技术摘要：窄均匀超导带中冻结通量的理论研究

问题陈述
本文研究了在恒定垂直磁场（$H$）中，长而窄且均匀的超导带在穿过其转变温度（$T_c$）时产生的残余捕获磁通量问题。这一现象对于超导电子器件（如 SQUID、transmon 量子比特和单光子探测器）的性能至关重要。虽然实验成像已经证实，在低于某一“有效通量驱逐场”（$H_{exp}$）的条件下不存在捕获通量，但该实验值通常比理论上的最小通量驱逐场（$H_e$）高出 3–4 倍。在动态过程中存在一个根本性的理论空白：具体而言，是什么决定了冷却后的冻结通量密度？系统是如何从平衡态过渡到冻结态的？以往的工作缺乏对这种冻结态形成的定量分析。

方法论
作者采用动态平衡法来模拟冷却过程中涡旋密度（$n$）的演化。分析重点在于带宽 $W$ 远小于 Pearl 屏蔽长度（$\Lambda$）的机制，这使得解析处理涡旋能量剖面成为可能。

1. 能量景观： 研究利用了窄带中已知的涡旋能量剖面 $E(X)$，该剖面包含了核能量、与带边缘的相互作用以及 Meissner 屏蔽电流的贡献。该剖面定义了涡旋退出的能量势垒 $U_{out}$ 和涡旋进入的能量势垒 $U_{in}$。
2. 动态平衡方程： 涡旋密度的随时间演化由描述热激活跳跃的速率方程控制：
   $$\frac{dn}{dt} = -n \frac{1}{\tau} e^{-U_{out}/T} + n_\xi \frac{1}{\tau} e^{-U_{in}/T}$$
   其中 $\tau$ 是 Kramers 尝试时间，$n_\xi$ 是与最大可能线性密度相关的前指数因子。
3. 冷却方案： 假设系统以恒定速率 $R = -T_c^{-1} dT/dt$ 进行冷却。通过求解动态方程来确定涡旋密度何时脱离与快速变化的温度的平衡状态。
4. 冻结条件： “冻结”通量的定义并非是一个突发的过程，而是指弛豫时间（随能量势垒呈指数级变化）变得过长，以至于无法追踪不断变化的平衡密度。通过令激活速率等于冷却诱导的进入 Boltzmann 因子的变化率，推导出冻结时间 $t_{fr}$ 和冻结温度 $T_{fr}$。

主要贡献与结果

• 冻结温度（$T_{fr}$）的推导： 本文推导出了约化冻结温度 $\varepsilon_{fr} = 1 - T_{fr}/T_c$ 的定量表达式。分析表明，$\varepsilon_{fr}$ 比转变过程中的涨落宽度（$\varepsilon_f$）高出一个巨大的对数因子。至关重要的是，当施加磁场 $H$ 接近最小通量驱逐场 $H_e$（即能量势垒消失处）时，$\varepsilon_{fr}$ 会迅速增加；且随着冷却速率降低，$\varepsilon_{fr}$ 呈对数增长。
• 冻结通量密度（$n_{fr}$）： 获得了一个关于冻结涡旋密度的解析表达式。对于任何 $H > H_e$ 的情况，该密度均为有限值，但在 $H_e$ 附近呈指数级减小。该密度表现出极强的磁场依赖性，随着 $H$ 的增加而剧烈上升。
• 标度律： 尽管涉及复杂的材料参数（相干长度 $\xi$、穿透深度 $\lambda$、粘度），但比值 $\varepsilon_{fr}/\varepsilon_f$ 和 $n_{fr}/n_\xi$ 仅取决于两个约化参数：无量纲冷却速率（$R\tau_0/\varepsilon_f$）以及带宽与涨落宽度处相干长度的比值（$2w\sqrt{\varepsilon_f}/\xi_{GL}$）。
• 有效通量驱逐场（$H_{eff}$）： 本文建议将有效通量驱逐场定义为计算出的冻结通量密度达到可检测水平（例如每条带长度一个涡旋）时的磁场。理论预测，$H_{eff}$ 总是大于理论上的 $H_e$。
• 特定材料的估算： 该理论被应用于铌（Nb）和铌钛氮化物（NbTiN）带。
  • 对于 Nb 带，计算出的 $H_{eff}$ 接近于交叉场 $H_s$（即单涡旋近似失效的场），但仍显著低于实验值（$H_{exp}$）。作者认为这种差异可能源于大规模的样品不均匀性，而非理想带模型的失效。
  • 对于 NbTiN 带，计算出的 $H_{eff}$ 与实验数据基本吻合，虽略高。较大的伦敦穿透深度和较小的相干长度使得 NbTiN 具有更宽的涨落机制，并表现出与 Nb 不同的标度行为。

意义与主张
本文声称提供了第一个关于在固定磁场下冷却的窄超导带中冻结通量密度的明确定量预测。通过求解动态平衡方程，作者确立了最终状态是由热激活与冷却速率之间的竞争决定的，这一过程发生在非常接近 $T_c$ 但处于弱涨落机制内的温度下。

其主要意义在于解释了“有效通量驱逐场”的起源。论文指出，实验观察到的 $H_{exp}$ 不是一个基本的热力学极限，而是一个由检测灵敏度和特定冷却历史定义的动力学阈值。该模型成功捕捉到了捕获通量的陡峭磁场依赖性，这表明在某些材料（如 Nb）中，理论上的 $H_e$ 与实验上的 $H_{exp}$ 之间存在的巨大差异，可能是由于样品不均匀性（钉扎）引起的，而非窄带理论本身的失效。这项工作为解释实验通量捕获数据以及基于材料参数和冷却速率定义有效驱逐场提供了一个框架。