超导电流为什么也会有方向偏好？

零电阻，并没有消除状态的条件

假设一个器件在正向电流达到某个值时开始出现电压，而反向需要更大的电流才出现电压。把同样大小的电流分别向两个方向加上去，就可能得到不同的结果：一个方向已经耗散，另一个方向仍处在低耗散状态。

这里的方向性通常只在一个有限区间内成立。电流很小时，两边都可能保持超导；电流过大时，两边都可能失去原来的低电压状态。“超导二极管”并不意味着所有大小的电流都只能向一个方向通过。1

超导状态由许多电子协同形成，能够支持无直流电压的持续电流，但这种状态有自己的稳定条件。电流会改变凝聚态的相位梯度，也会对磁通涡旋施力，改变弱连接中的相位差。当某个维持条件先被破坏，就可能出现电压。

问题因此变得具体：同一器件在正反两个方向，究竟是哪一项条件先到了极限？

先把“临界电流”拆开

几种常见电流读数，在不同论文中可能都被宽泛地称为临界值，但它们对应的过程不同。

相同方向的开关值与恢复值也可能不同，这叫回滞。它本身不等于正反方向不对称。另一方面，即使曲线上有很明显的方向差，也未必直接量到了材料的退配对极限。16

为免符号混乱，下文用$I_c^+$和$I_c^-$分别表示正反向阈值的正数幅度，并定义：

它衡量两边相差多少，不是输出能量除以输入能量。论文常把它称为“二极管效率”，阅读时仍应记住这个实际定义。不同实验的$I_c$取法不同，直接比较它们的最大$\eta$也需要相应条件。

相位关系怎样产生不同门槛

约瑟夫森结是两个超导部分之间的弱连接。即使中间有很薄的绝缘层或一段正常金属，两侧也可以通过量子耦合传递超电流。影响电流的一个核心量，是两侧凝聚态之间的规范不变相位差$\phi$；有磁场时，它还包含矢势的贡献。2

最简单的关系是$I(\phi)=I_0\sin\phi$。正向最大值与反向最小值的幅度一样，所以没有这种阈值方向性。

若把整条曲线平移成$I_0\sin(\phi+\phi_0)$，正反极值仍然一样。平移改变了相位位置，却没有改变上下两端的高度。

甚至“出现更高次谐波”也不自动够用。本篇比较三个无量纲例子：

实际求极值后，第三种关系的正向门槛为0.8、反向门槛为1.2，因而$\eta=-20\%$；前两种均没有阈值差异。第三种的导数还能写成$\cos\phi(1-0.8\sin\phi)$，括号始终为正，所以这两个极值可以直接解析核对。

这些函数展示的是相位关系的后果，没有预先解释材料为什么产生某个谐波。要完成机制解释，还需要找出哪些物理环节改变了曲线。

另外，$I_C(0)$不为零，并不意味着一个没有外部驱动的孤立器件会凭空持续输出电流。相位能够自行调整到允许的平衡位置；外界怎样固定相位、磁通或电流，属于问题的一部分。

对称性给出约束，还要找到生成机制

一种对称变换如果能在保持同一实验条件的同时，把正向状态变成对应的反向状态，就会限制两者的差异。

Wang等的分析把这件事写成自由能关系。若在固定磁场$B$下存在常数$\theta$，使

那么由$I=(2e/\hbar)\partial F/\partial\phi$可以推出成对的相反电流，以及相等的正反极值。这比只检查某一个熟悉的空间或时间对称性更直接。2

因此，打破空间反演和时间反演，通常只是排除了一些禁止条件。剩余的组合对称、相位平移，乃至特定参数下的抵消，仍可能保留相等的极值。证明“允许不同”，没有完成“为什么实际不同”的工作。

这些自由能约束针对相应的静态问题。若实验主要测的是耗散运动怎样停止，就需要把阻尼与历史放进来，不能机械套用同一个判据。

电流反过来改变相位：一个能走通的反馈例子

2024年Chen等研究薄膜约瑟夫森器件时，发现电流本身对相位分布的影响十分重要。薄膜中的超流有动能；电流变化对应动能储存的变化，可以用动能电感描述。在超流刚度较小的薄膜中，它能显著影响相位，即使电流自己产生的磁场还很弱。3

论文用两个并联的小结说明机制。两条通道有不同强度，外加磁通使它们的相位错开，电流又反过来增加相位差。此时，必须同时求电流和它造成的相位改变。

将两条通道的电流记作$I_1$、$I_2$，它们之间的相位错开量记作$\chi$，这个关系可以写为：

其中$f$是外加磁通贡献的相位，$k$控制电流对相位的反馈。最后一行又会改变前两行，所以不能先固定相位图样，再把两条电流简单相加了事。

本篇按附件的两个隐式方程实际求解，选取通道强度1.1和0.9、反馈系数0.4。这里的参数均为无量纲示意值，选择在唯一解容易验证的范围内，没有拟合真实器件。

计算结果显示：有外加磁通、两侧强度不同且保留反馈时，正反极值可以不同；交换两侧强度会反转偏好。保持几何不对称但去掉反馈，或保留反馈但让两侧相同，则恢复相等的极值。

这里新增的解释是：偏置电流改变空间分布，空间分布再改变相位，整体电流关系因此超出了简单正弦。它连接了前一节抽象的谐波条件与一个具体的器件过程。

这仍只是两通道约化模型。所选参数使相位方程的迭代具有收缩性，能验证唯一解；更强反馈、多稳态或连续空间器件，需要额外判断分支的可达性和稳定性，不能直接套用本图。

为什么局部磁成像比一条电阻曲线多告诉我们一些事

Chen等用钻石尖端中的单个氮空位中心测量局部磁场，再反演电流分布。论文报告两个器件；主要比较7 K和4 K下的电流、磁场与空间图样。3

这项测量的直接信号是两种偏置条件下的磁场差，电流箭头则需要反演。作者比较了傅里叶、正则化与机器学习等重建方法；测量窗口以外的电流、边界延拓和空间分辨率仍属于反演条件。图像提供了额外证据，但不是一台直接拍到每个电子的相机。

在4 K下，输运测得正反阈值不同；局部图样也显示不对称。更有意义的是，即使仍处在输运的低电阻区，改变偏置也能切换不同的约瑟夫森涡旋配置。单看电阻，会漏掉这些内部状态的变化。

作者测得电流自场相对外场较小，并结合薄膜电感与模拟，支持电流诱导相位的重要作用。自场、结宽变化、界面透明度不均匀等解释仍需逐项区分；“存在不均匀”本身不能决定不均匀来自哪里。零电阻区内观察到不同状态，也没有量出整个切换过程的零能耗。

本篇取得作者公开数据包，将4 K图4a中的221个磁场点逐行重算，再按原表重绘：

这张图说明，方向偏好随外场变化，不能只用一个最大值概括。小电流区域的分母较小，同样大小的阈值误差会造成更大的比例摆动。221个点是场扫描数据，不是221个独立器件的复现。

附件另把7 K下的结果描述为未观察到二极管效应。不过其公开源表存在列名、部分电流符号和比例不一致的问题，本篇已登记，未用这些数值证明“效应严格为零”。本篇关于“几何不对称不一定产生阈值差异”的判断，由前面的解析与反馈对照支撑，不依赖这一较强的零效应表述。

涡旋路线：材料仍能配对，电压却可以先出现

磁通涡旋是超导状态中的一种空间结构，周围相位发生绕转，核心附近的超导受到抑制。电流对涡旋施加洛伦兹力；若涡旋开始穿过样品，持续的相位变化就能产生电压。因此，测到电压并不总意味着整个材料已经同时失去配对。4

2023年Gutfreund等在EuS/Nb双层器件中，把局部磁成像与输运放在一起比较。铁磁层的磁化会在边缘产生杂散磁场，超导层为屏蔽它而产生电流。这些电流与外加输运电流叠加后，能够在一个方向帮助涡旋通过，在另一个方向抑制通过。

关键对照是改变磁化方向。横向磁化翻转时，可见的涡旋通道换到另一侧，输运不对称也翻转；在净磁化较小的条件或磁化转到平行电流方向时，这种差异明显减弱。空间信号、磁化条件与输运读数相互对应，使边缘屏蔽电流的解释比仅凭一条非对称曲线更有依据。

成像平均了大量快速涡旋运动，不能把它看作单个涡旋的连续高速录像。实验还施加了垂直磁场以提供涡旋；不能把结果搬到完全无磁通的器件条件下。

附件进一步说明，主图的临界值取自平滑后的$dV/dI=20$ mΩ门槛；换成2 μV电压门槛，趋势相同但散点更多。作者明确承认，这种做法主要描述曲线不对称，并不必然精确定位宏观磁通流动的唯一物理起点。这项限制没有抹掉空间与磁化对照的价值，却限制了临界值的具体含义。4

2026年的层状模型：改变排列，能改变运动路径

Li等在2026年用时变Ginzburg–Landau模型研究Nb、V、Ta层状结构，把材料参数差异、磁场、电流、涡旋运动与热方程连在一起。5

在其选定结构中，电流反向改变洛伦兹力方向，各层的超导状态、涡旋占据和穿层过程随之不同。某一方向更容易出现剧烈涡旋运动与快速升温，因此更早转入高电阻状态。

例如，原文在4 K、归一化磁场为0.08时，一个三周期结构的模拟给出正反电流密度阈值幅度约0.0090和0.0065，按差与和计算为16.1%。这些是论文自己的无量纲电流密度。另一个交换层序、归一化磁场为0.10的例子，两侧都约为0.0061，方向差异在其分辨率下消失。本篇重算了原文五组打印数值的比例，但没有重新运行这些空间模拟。

值得继续研究的机制预期是：层的排列、单层厚度、磁场与散热条件共同改变运动路径。这里的“交换层序”也不是简单把原器件整体镜像翻转；文中具体交换了Nb与V的位置，改变了相邻材料与边界组合。

这项工作有几个重要边界。它是带热反馈的二维模型，所用单层远厚于早期1 nm超晶格实验；按其长度单位与参数表换算，90个相干长度对应约3.51 μm，不能当作120 nm原实验的精确重建。热学参数还被简化为统一的现象学系数。

作者还用自由能分布解释涡旋的偏好方向。理解时要进一步分清：不同材料背景的凝聚自由能，与在指定背景中增加或移动一个涡旋的能量代价，并非同一个量。空间演化提供了该模型中的机制线索，但不能仅凭不同材料自由能颜色的高低，就为所有层状器件确定一个普遍涡旋势能排序。

模型中稳定磁通流动区间很短，也不自动证明整流器件总能耗很低。该模型本身就包含涡旋运动、焦耳热与快速升温；实际性能还要计算脉冲过程、负载、散热和制冷成本。

恢复过程也能有方向性，而且原因可能在耗散

Trahms等在2023年研究单磁原子约瑟夫森结，得到另一种有区分力的结果：加入Cr或Mn原子后，最大的正反差异出现在恢复低电压状态的电流上，开关电流的差异小得多。两种原子的恢复不对称方向还相反。6

这类小结不能只看静态相位关系。相位运动还受到电容储能、准粒子电流、环境阻尼和热噪声影响。减小电流时，已经运动起来的相位可能继续运动，直到能量输入不足以抵消损失，才重新被困在某个状态。

可以用一个地形类比帮助理解：开始越过山坡，取决于坡形、推动和扰动；已经跑起来之后何时停住，还取决于运动中损失多少能量。类比中的地形对应相位关系，摩擦对应阻尼，它们在真实器件里由不同物理量决定。

原实验测到磁原子引入的亚能隙谱峰在正反偏置下具有不同强度。作者将相关准粒子电流带入含噪声和频率依赖阻尼的模型，在保持其他主要参数相同的比较下，再现了恢复电流不对称的方向。改用非对称相位关系、保持耗散部分对称，得到的却主要是开关电流的不对称，与实验的主特征不合。

这个对照说明，必须同时解释“哪种阈值变了”。它并没有取消静态超电流的对称性约束，因为它着重研究的是不同的动力学过程。

该器件的低电压分支还有残余相位扩散，不能严格称为理想零电压。开关和恢复都具有随机性；论文按连续100次扫描汇总，并检查尖端漂移、导电性变化和仪器偏置。扫描次数也不等于独立器件数。6

能整流，也要有能量收支

回到本篇第三种相位关系，另加一个零噪声、无电容的阻尼模型：

这里使用无量纲量，平均相位速度对应平均电压。在允许静态解的区间，相位能够停住；超过相应极值时，相位持续绕转，出现平均电压。本篇对一个周期的运动时间积分，并用直接时间积分交叉核对。

幅度为1时，正向超过0.8的门槛，反向还没超过1.2的门槛。若以很慢的方波在两者之间切换，就可以产生非零平均电压。但有电压的半周期也在消耗输入功率；方向性没有绕过能量守恒。

实际器件中的涡旋、恢复过程、热反馈和有限切换速度会增加结构。本例只展示相位关系与耗散的基本连接，没有给出某台装置的工程效率。

怎样把相似现象分出不同机制

较有用的研究路线，是让证据逐层连接：先确定哪种电流读数有差异，再检查局部状态如何变化，最后比较改变条件后的预期。

• 若主要解释是相位关系与电感反馈，应同时比较正反极值、局部电流分布，以及改变通道强度和反馈后的结果。
• 若主要解释是涡旋路径，应追踪磁化、边缘屏蔽、入射与退出条件，并检查空间信号是否与电压起点对应。
• 若主要解释是恢复阻尼，应同时比较开关和恢复分布、准粒子谱、扫描历史及环境阻抗。
• 若提出一种器件性能优势，还应测输入能量、完整切换过程、负载和散热，而不能只引用阈值差异比例。

这些机制可以共存。一个器件同时具有材料不对称、边缘不均匀和热反馈时，单一曲线未必足以给每项贡献分账。

超导电流的方向偏好，来自结构、状态与运动过程对两个方向给出了不同条件。 对称性帮助识别哪些差异被禁止，具体模型说明差异怎样产生，而区分性测量决定这套解释在某台真实器件中能走多远。

可回查：4 K原表与重算比例、三种相位曲线、两结反馈计算、电压与功率、原文TDGL打印值比例、计算结果及来源表缺口。