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

本論文は、絶対零度における超流動ヘリウム4中を運動する物体によるロトン放出および渦核生成の臨界速度に圧力と障害物サイズがどのように影響するかを解析するために一般化された非局所グロス・ピタエフスキーモデルを用いた数値研究を提示し、これら圧力依存メカニズムを同時に扱う最初の理論的枠組みを確立するものである。

要約

摩擦を一切伴わず、まるで指の間をすり抜ける幽霊のように流れる液体の世界を想像してみてください。これが**超流動ヘリウム（He II）**です。これはヘリウムを絶対零度に近い温度まで冷却したときにのみ現れる、物質の特別な状態です。

この論文は、まるでハイテクなシミュレーション実験室のようであり、著者たちはこの幽霊のような液体を試験にかけました。彼らが目指したのは、小さな物体（微小な気泡や帯電粒子など）をこの超流動体の中を通そうとしたときに何が起こるかを理解することでした。具体的には、液体がその完全な流れを「壊す」前に、物体をどの程度の速度まで押し進めることができるのかを知りたかったのです。

以下に、彼らの発見の物語を、簡単な概念に分解して紹介します。

1. 液体が壊れる二つの方法

超流動ヘリウムの中に物体を押し進めると、液体は水のように単に「引きずられる」わけではありません。ある特定の速度限界（臨界速度と呼ばれるもの）に達すると、液体は二つの明確な反応を示します。

• 「ロトン」の爆発: 液体を穏やかな海だと想像してください。物体を押し進めすぎると、単に波を立てるだけでなく、突然ロトンと呼ばれる微小でエネルギーに満ちた粒子の群れが生まれます。まるで液体が突然、百万個の小さなエネルギーの火花に砕け散ると決めたかのようです。これは特定の速度で起こります。
• 「渦」の渦巻き: さらに速く押し進める（あるいは物体が十分に大きい）と、液体は回転し始めます。それは量子渦と呼ばれる微小なマイクロレベルの竜巻を作り出します。これらは物体に絡みつき、引きずり下ろす小さな渦巻きのようなものです。

この論文の主な目的は、その「火花」（ロトン）と「渦巻き」（渦）のどちらを誘発するために、どの程度の速度が必要かを正確に突き止めることでした。

2. 圧力鍋実験

著者たちは、液体を一つの圧力だけで観察したわけではありません。真空（0 バール）から、固体の岩に変わる点（約 25 バール）まで、ヘリウムを締め付けていく過程で何が起こるかをシミュレーションしました。

彼らは、**「一般化された非局所グロス・ピタエフスキーモデル」**と呼ばれる特別な数学的モデルを使用しました。これはスーパーに正確なビデオゲームエンジンとして機能します。このエンジンは、標準的な物理方程式では通常見逃されてしまう奇妙な「ロトン」の振る舞いを含め、ヘリウム原子の実際で複雑な挙動を模倣するようにプログラムされていました。

3. 大きな発見：締め付けがルールを変える

彼らが発見したことを、簡単な比喩を使って説明します。

あなたが人混みの中を走ろうとしていると想像してください。

• 低圧（緩やかな人混み）: 人々は広がっています。彼らは遠く離れているため、暴動（ロトンの生成）を起こすのは実際にはかなり困難です。しかし、十分に速く走れば、誰かを転ばせて、転び続ける人々の連鎖反応（渦）を起こすかもしれません。
• 高圧（ぎゅうぎゅうの人混み）: 人々は肩を並べて詰め込まれています。今や、彼らが非常に近接しているため、暴動（ロトン）を起こすのははるかに容易です。しかし、人混みが非常に密度が高く硬いため、回転に抵抗し、転び続ける人々の連鎖反応（渦）を起こすのは困難になります。

結果:

• ロトン速度: ヘリウムを締め付ける（圧力を上げる）につれて、その「火花」（ロトン）を作るために必要な速度は低下しました。流れを壊すために、それほど速く走る必要がなくなったのです。
• 渦速度: ヘリウムを締め付けるにつれて、「渦巻き」（渦）を作るために必要な速度は上昇しました。液体を回転させるには、はるかに速く走る必要があります。

4. 検出のための「絶妙なポイント」

これにより、興味深い隙間が生まれます。高圧下では、「火花」（ロトン）を「渦巻き」（渦）を作ることなく生成できる、広範な速度範囲が存在します。

過去、科学者たちは、しばしば厄介な渦に隠れてしまうロトンを研究することに苦労していました。著者たちは、ヘリウムを高圧に締め付けることで、ロトンが単独で現れる「クリーンな」環境を作り出すことができ、それによりロトンの研究がはるかに容易になると提案しています。

5. サイズが重要

この論文は、液体の中を移動する物体の大きさにも注目しました。

• 微小な物体（単一のイオンのようなもの）: これらは非常に敏感です。彼らはまず「ロトン限界」に達します。
• 大きな物体（大きな円盤のようなもの）: これらはロトンに対してあまり敏感ではありません。圧力に関係なく、彼らはまず「渦限界」に達する傾向があります。

まとめ

著者たちは、圧力下での超流動ヘリウムを観察するためのデジタル顕微鏡を構築しました。彼らは、ヘリウムを締め付けることで、エネルギーの火花（ロトン）を作りやすくする一方、回転する渦巻き（渦）を作るのは難しくなることを発見しました。

これは、過去の実験が異なる圧力で異なる挙動を見た理由を説明し、もし私たちが神秘的な「ロトン」粒子を研究したいのであれば、液体が渦巻きで汚れることなく秘密を明かす可能性が高い高圧下で実験を行うべきであることを示唆しています。

注：著者らは、完全な 3 次元で行うと計算量が重すぎるため、シミュレーションは 2 次元（世界の平らなスライス）で行われたことを認めています。しかし、彼らは発見した物理法則が実際の 3 次元の世界でも真実であると信じています。

技術概要

技術的サマリー：極低温（T→0）極限における He II 中を運動する物体によるロトン放出および渦の生成の圧力・サイズ依存性

問題提起
ヘリウム II（He II）中を物体が運動する際、超流動性の崩壊は主に 2 つのメカニズムによって支配される。すなわち、ロトン（素励起）の放出と、量子化渦の核生成である。ロトン放出の臨界速度（$v_{rot}^{crit}$）および渦生成の臨界速度（$v_{v}^{crit}$）の圧力依存性は、特にランカスター研究グループが負イオンを用いて実験的に調査されてきたが、これらの現象は歴史的に統一的な理論的記述を欠いてきた。実験的観測によれば、低圧力領域では渦の核生成がイオンの速度を制限するが、高圧力領域（約 10〜12 バール以上）ではロトン放出が制限要因となる。さらに、ロトン放出の臨界速度は圧力の増加とともに低下するのに対し、渦の核生成の臨界速度は圧力の増加とともに上昇する。標準的なグロス・ピタエフスキー方程式（GPE）などの既存の理論的枠組みは、厳密に単調な励起スペクトルしか持たず、特徴的な「ロトン」極小値を持つランドウ分散曲線を再現できないため、ロトン放出を捉えることができない。一方、微視的理論は、しばしば単一の枠組み内で両方のメカニズムを記述することに苦慮している。

手法
著者らは、絶対零度（T→0）における超流動$^4$He のダイナミクスをシミュレートするために、一般化された非局所グロス・ピタエフスキー（gnlGPE）モデルを採用した。このモデルは複素秩序パラメータ$\psi(x,t)$によって定義され、非局所的な原子間相互作用ポテンシャル$V(r)$および高次の密度補正項を組み込んでいる。

1. モデル較正: 物理的な精度を圧力範囲（0 から 24 バール、凝固圧力である約 25 バールに近づく）全体にわたって確保するため、著者らは相互作用ポテンシャルに固定された解析的形式を用いない。代わりに、Godfrin ら（2021）の実験データを用いて一般化されたボゴリューボフ分散関係を逆算する。これにより、モデルはフォノン、マクソン、ロトンの各分枝およびピタエフスキー高原を含む、圧力依存性を持つ正確な分散曲線を再現できる。
2. シミュレーション設定: 本研究は主に 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$バールで行われる。障害物のサイズは$D = 10\xi$から$D = 1000\xi$まで変化する。特に、圧力とともに変化する負イオンのサイズ（0 バールで約 34.5 オングストロームから 24 バールで約 21.7 オングストローム）に注意が払われる。

主要な結果

• 臨界速度の圧力依存性: シミュレーションは実験的に観測された相反する傾向を確認する。圧力が上昇するにつれて、ロトン放出の臨界マッハ数（$M_{rot}^{crit}$）は低下し、0 バールで約 0.247 から 24 バールで約 0.128 まで減少する。一方、渦の核生成の臨界マッハ数（$M_{v}^{crit}$）は圧力の増加とともに上昇する。
• サイズ依存性: 小さな障害物（負イオンのサイズに相当し、$D < 100$ nm）の場合、臨界速度は 2 つの閾値（$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 を利用することで、本研究はランドウ基準（ロトン放出）と流体力学的な渦の核生成の間のギャップを成功裏に埋めている。

本論文は、3D シミュレーションの計算複雑性（特に 3D における渦核生成のための最小シードの決定に関連して）により研究が 2D に制限されていることを強調しているが、2D の結果は実験的観測と整合する基礎物理の定性的および半定量的な説明を提供すると述べている。著者らは、臨界速度の数値が実験データと比較可能であることを指摘し、実験的に導出された分散関係を用いて相互作用ポテンシャルを制約するアプローチの有効性を検証している。この研究は、臨界圧力$P_{crit} \approx 10-12$バールでイオンの運動の性質が変化する理由に対する物理的に動機付けられた説明を提供し、それを$v_{rot}^{crit}$および$v_{v}^{crit}$の圧力依存曲線の交差に帰因している。

著者らは 3D への含意については控えめであり、2D シミュレーションが一般的な傾向を説明する一方で、3D における渦核生成の具体的なトポロジー（例えば、イオン・渦ループ対リング）は今後の調査を必要とする未解決の問題であると認めている。また、イオンを完全な硬球としてモデル化することは簡略化であるが、圧力の傾向は頑健であると予想されるとも述べている。