Pressure and Size Dependence of Roton Emission and Vortex Creation by Moving Objects in He~II in $T \to 0$ Limit: Generalized Nonlocal Gross-Pitaevskii Model

本文采用广义非局域 Gross-Pitaevskii 模型进行数值研究，分析了在零温下压力与障碍物尺寸如何影响运动物体在超流氦-4 中诱发旋子发射和涡旋成核的临界速度，这标志着首个能够同时处理这些压力依赖机制的理论框架。

要点

想象一个液体可以毫无摩擦地流动的世界，就像幽灵从你的指缝间滑过。这就是超流氦（He II），一种仅在氦被冷却至接近绝对零度时才会存在的特殊物质状态。

本文就像一个高科技模拟实验室，作者们将这种幽灵般的液体置于考验之中。他们想要了解，当你推动一个小物体（如微小气泡或带电粒子）穿过这种超流体时会发生什么。具体来说，他们想知道：在液体“打破”其完美流动之前，你最多能推多快？

以下是他们发现的简化解说：

1. 液体破裂的两种方式

当你推动一个物体穿过超流氦时，液体并不会像水那样仅仅受到“拖拽”。相反，一旦你达到某个速度极限（称为临界速度），它会以两种截然不同的方式做出反应：

• “旋子”爆发： 把液体想象成一片平静的海洋。如果你推动物体的速度太快，你产生的不仅仅是波浪；你会突然制造出一群微小的、高能的粒子，称为旋子（rotons）。这就像液体突然决定碎裂成百万个微小的、充满能量的火花。这种情况发生在特定的速度下。
• “涡旋”漩涡： 如果你推得更快（或者物体足够大），液体开始旋转。它会形成微小的、微观的龙卷风，称为量子涡旋（quantum vortices）。这些就像附着在物体上的微小漩涡，将其向下拖拽。

本文的主要目标是弄清楚，究竟需要多快的速度才能触发“火花”（旋子）， versus 触发“漩涡”（涡旋）。

2. 高压锅实验

作者们并没有只观察单一压力下的液体。他们模拟了随着氦被挤压得越来越紧所发生的情况，从真空（0 巴）一直增加到它会变成固体岩石的程度（约 25 巴）。

他们使用了一种特殊的数学模型（“广义非局域 Gross-Pitaevskii 模型”），它就像一个超精确的视频游戏引擎。这个引擎被编程以模仿氦原子的真实复杂行为，包括那些标准物理方程通常忽略的奇怪“旋子”行为。

3. 重大发现：挤压改变规则

以下是他们的发现，使用一个简单的类比：

想象你正试图穿过一群人的拥挤人群。

• 在低压下（松散的人群）： 人们分散开来。实际上很难引发骚乱（产生旋子），因为他们相距甚远。但是，如果你跑得足够快，你可能会绊倒某人，引发一连串人们摔倒的连锁反应（涡旋）。
• 在高压下（拥挤的人群）： 人们肩并肩地挤在一起。现在，引发骚乱（旋子）要容易得多，因为他们靠得太近了。然而，引发一连串人们摔倒的连锁反应（涡旋）变得更困难，因为人群如此密集且僵硬，以至于抗拒旋转。

结果：

• 旋子速度： 随着他们挤压氦（增加压力），产生那些“火花”（旋子）所需的速度下降了。你不必跑得那么快就能打破流动。
• 涡旋速度： 随着他们挤压氦，产生“漩涡”（涡旋）所需的速度上升了。你必须跑得更快才能让液体旋转。

4. 探测的“甜蜜点”

这创造了一个迷人的间隙。在高压下，存在一个广泛的速度范围，在这个范围内你可以产生“火花”（旋子）而不产生“漩涡”（涡旋）。

过去，科学家们难以研究旋子，因为它们经常被混乱的漩涡所掩盖。作者们建议，通过将氦挤压到高压，我们可以创造一个“干净”的环境，让旋子独自出现，从而使它们更容易被研究。

5. 尺寸很重要

本文还考察了穿过液体运动的物体的大小。

• 微小物体（如单个离子）： 它们非常敏感。它们首先达到“旋子极限”。
• 大物体（如大圆盘）： 它们对旋子不太敏感。无论压力如何，它们往往首先达到“涡旋极限”。

总结

作者们构建了一个数字显微镜来观察受压下的超流氦。他们发现，挤压氦使得产生能量火花（旋子）变得更容易，但产生旋转漩涡（涡旋）变得更困难。

这解释了为什么过去的实验在不同压力下观察到了不同的行为，并表明如果我们想要研究神秘的“旋子”粒子，我们应该在高压下进行实验，在那里液体更有可能在不被漩涡搞乱的情况下揭示其秘密。

注：作者承认他们的模拟是在二维（世界的扁平切片）中完成的，因为在完整的三维中进行计算过于繁重，但他们认为他们发现的物理规律在真实的三维世界中同样成立。

技术摘要

技术摘要：零温极限下运动物体在 He II 中发射旋子和产生涡旋的压力与尺寸依赖性

问题陈述
当物体在超流体氦 II（He II）中运动时，超流性的破坏主要由两种机制主导：旋子（基本激发）的发射和量子化涡旋的成核。尽管旋子发射的临界速度（$v_{rot}^{crit}$）和涡旋产生的临界速度（$v_{v}^{crit}$）的压力依赖性已通过实验进行了研究——特别是兰开斯特小组利用负离子进行的研究——但这些现象在历史上缺乏统一的理论描述。实验观察表明，在低压下，涡旋成核限制了离子速度，而在高压下（高于 $\approx 10-12$ bar），旋子发射成为限制因素。此外，旋子发射的临界速度随压力增加而降低，而涡旋成核的临界速度则随压力增加而升高。现有的理论框架，如标准 Gross-Pitaevskii 方程（GPE），无法捕捉旋子发射，因为它们具有严格单调的激发谱，无法重现具有特征性“旋子”极小值的朗道色散曲线。相反，微观理论往往难以在单一框架内描述这两种机制。

方法论
作者采用广义非局部 Gross-Pitaevskii（gnlGPE）模型来模拟零温（$T \to 0$）下超流体 $^4$He 的动力学。该模型由复序参量 $\psi(x,t)$ 定义，并包含非局部原子间相互作用势 $V(r)$ 和高阶密度修正项。

1. 模型校准：为了确保在一系列压力（0 至 24 bar，接近约 25 bar 的固化压力）下的物理准确性，作者未对相互作用势使用固定的解析形式。相反，他们利用 Godfrin 等人（2021）的实验数据，通过反演广义 Bogoliubov 色散关系来构建模型。这使得模型能够重现精确的压力依赖性色散曲线，包括声子、最大子和旋子分支，以及 Pitaevskii 平台。
2. 模拟设置：研究主要在二维（2D）近似下进行。作者模拟了一个在静止 He II 中运动的圆形障碍物（圆盘）。该障碍物被建模为具有高斯分布的外部势，能够完全耗尽超流体密度。
3. 临界速度确定：
   • 旋子发射（$v_{rot}^{crit}$）：通过分析稳态解的分岔图确定。作者使用牛顿 - 拉夫逊法耦合双共轭梯度稳定法，寻找存在稳定稳态解的最大马赫数（$M = U/c$）。超过这一点，系统变得不稳定，旋子被发射。
   • 涡旋成核（$v_{v}^{crit}$）：通过动态方法确定。障碍物被加速直至达到目标速度，随后系统演化约 $10\xi/c$ 的时间尺度（其中 $\xi$ 为愈合长度），以观察是否成核了量子化涡旋。
4. 参数：模拟涵盖压力 $P = 0, 2, 5, 10,$ 和 $24$ bar。障碍物尺寸从 $D = 10\xi$ 变化至 $D = 1000\xi$。特别关注了随压力变化的负离子尺寸（从 0 bar 时的 $\approx 34.5$ Å 到 24 bar 时的 $\approx 21.7$ Å）。

关键结果

• 临界速度的压力依赖性：模拟证实了实验观察到的相反趋势。随着压力增加，旋子发射的临界马赫数（$M_{rot}^{crit}$）降低，从 0 bar 时的 $\approx 0.247$ 降至 24 bar 时的 $\approx 0.128$。相反，涡旋成核的临界马赫数（$M_{v}^{crit}$）随压力增加而升高。
• 尺寸依赖性：对于小障碍物（与负离子尺寸相当，$D < 100$ nm），临界速度由两个阈值（$v_{rot}^{crit}$ 或 $v_{v}^{crit}$）中较低者决定。在低压下，$v_{v}^{crit}$ 是限制因素；在高压下，$v_{rot}^{crit}$ 是限制因素。随着障碍物尺寸增大（$D \ge 1000\xi$），所有压力下的临界速度汇聚到一条由涡旋成核主导的单一曲线上，这与标准 GPE 的预测相似。
• “间隙”：一个重要的发现是，随着压力增加，$v_{rot}^{crit}$ 和 $v_{v}^{crit}$ 之间的差距扩大。在高压下，存在一个速度窗口，在此窗口内旋子被发射，但涡旋成核尚未发生。这表明在高压下可以隔离并实验观察旋子发射，而无需立即面对涡旋生成的复杂情况。
• 可视化：密度可视化证实，在中等马赫数下（例如 $M=0.28$），旋子在所有压力下均被发射，而涡旋成核仅发生在更高的马赫数下（$M \approx 0.3$），且随着压力升高，成核阈值向更高值移动。

意义与主张
作者声称，这项工作代表了首次在一个理论框架内分析旋子发射和涡旋成核及其相关临界速度的压力依赖性。通过利用广义非局部 GPE，该研究成功弥合了朗道判据（旋子发射）与流体动力学涡旋成核之间的差距。

论文强调，尽管由于三维模拟的计算复杂性（特别是关于确定三维中涡旋成核的最小种子），研究仅限于二维，但二维结果提供了与实验观察一致的、关于底层物理的定性和半定量解释。作者指出，他们计算出的临界速度数值与实验数据相当，验证了使用实验导出的色散关系来约束相互作用势的方法。这项工作为离子运动性质在临界压力 $P_{crit} \approx 10-12$ bar 处发生变化的原因提供了基于物理的解释，将其归因于 $v_{rot}^{crit}$ 和 $v_{v}^{crit}$ 的压力依赖性曲线的交叉。

作者对三维影响保持谦逊，承认虽然二维模拟解释了总体趋势，但三维中涡旋成核的具体拓扑结构（例如离子 - 涡旋环与涡旋环）仍然是一个有待未来研究解决的开放性问题。他们还指出，将离子建模为完美的硬球是一种简化，但预计压力趋势将保持稳健。