Vortices in Two-Dimensional Chiral Superfluids

本論文は、Bogoliubov-de Gennes理論を用いて多重量子化渦を持つ二次元カイラル超流体の軌道角運動量を調査し、角運動量がBEC領域では普遍的な公式に従う一方で、スペクトルの非対称性と対を組んでいないフェルミオンのためにBCS領域では著しく抑制されること、そしてその減少の程度が特定のペアリング対称性と渦の渦度（渦度）に依存することを明らかにしている。

要約

超流動体を、フェルミ粒子（粒子）が完璧に調和して動くためにペアを組む、巨大で目に見えないダンスフロアだと想像してみてください。 「カイラル」超流動体では、これらのペアはただ踊るだけでなく、特定の方向に回転しています。まるで、同期した一列のダンサーたちが全員で時計回りに回転しているようなものです。この論文は、このダンスフロアに「ひねり」や「渦（ボルテックス）」、つまりダンサーたちが中心点の周りを回転する渦巻きを導入したときに何が起こるのかを調査しています。

著者である Yan He と Wenxing Nie は、シンプルかつトリッキーな問いを投げかけています。「もしこのダンスフロアを回転させたら、システム全体の総『スピン』（軌道角運動量、または OAM）はどれくらいになるのか？」 ということです。

以下に、日常的な比喩を用いた彼らの研究結果の解説をまとめます。

1. 二つのダンススタイル：「簡単な」方法 vs 「難しい」方法

この論文では、ダンサーたちの二つの異なるレジーム（状態）について考察しています。

• BECレジーム（「密接な」ダンス）: ダンサーたちが手を取り合い、非常に固く結びついて、一つの単一のユニットとして振る舞っている状態を想像してください。この状態では、数学は単純です。もし、強さ $k$ の渦があり、ダンサーたちが自然に強さ $\nu$ で回転しているなら、部屋の総スピンは、予想通り「（$k + \nu$）× ダンサーの数」となります。これは完璧で、予測可能な計算です。
• BCSレジーム（「緩い」ダンス）: 今度は、ダンサーたちが手を緩く繋ぎ、かろうじてつながっているだけの状態を想像してください。彼らはより独立しています。この状態では、物事は複雑になります。論文によれば、総スピンは上述の「完璧な」数値よりも少なくなることがよくあります。

2. 消えたスピンの謎

なぜ「緩い」ダンスではスピンが消えてしまうのでしょうか？ 著者たちは スペクトル非対称性（Spectral Asymmetry）（またはスペクトルフロー）という概念を用いて説明しています。

ダンサーたちのエネルギー準位を、一連の階段だと考えてください。完璧な世界では、一人のダンサーが一段上がるごとに、別の誰かが一段下がることでバランスが保たれます。しかし、これらの渦を持つ超流動体では、この階段が乱れてしまいます。

• ペアを失ったフェルミ粒子: これらはパートナーを失ったダンサーたちです。彼らはグループと一緒に回転する代わりに、逆方向に回転します。
• 打ち消し合い: これらの「はぐれもの」のダンサーたちは逆方向に回転するため、ペアになったダンサーたちの前向きのスピンを打ち消してしまいます。これが、総スピンが減少する理由です。

3. 様々な種類の「ひねり」（渦）

論文では、二つの主要な変数、すなわち「ペアリングの強さ（p波、d波など）」と「渦の強さ（単一のひねりか、複数のひねりか）」をテストしています。

• 「完璧な」ひねり（単一のひねり、p波）:
  ダンサーたちがシンプルな「p波」のダンス（一度回転する）をしており、渦が単一のひれ（$k=1$）である場合、システムは美しく機能します。たとえ「緩い」ダンスのレジームであっても、総スピンは完璧なままです。「はぐれもの」のダンサーが現れて何も打ち消すことはありません。

  • しかし、 ひねりの中にさらなるひねりがあります。もし渦が反対方向（$k=-1$）に回転する場合、総スピンはゼロになります。しかし、論文は、たとえ総計がゼロであっても、スピンの分布は複雑であることを指摘しています。それは、部屋の半分が左に回り、もう半分が右に回っており、全体としては打ち消し合っているものの、局所的には動きが活発で、静かな部屋とは異なる状態にあるようなものです。

• 「混沌とした」ひねり（複数のひねり、または複雑なダンス）:
  もし、渦を二回以上回転させたり（$|k| \ge 2$）、あるいはダンサーたちがより複雑なダンス（自然に二回転するd波のような）を行ったりすると、「はぐれもの」のダンサーが現れます。

  • 複数のひねり（$|k| \ge 2$）: 「はぐれもの」のダンサーたちは、渦の中心部（コア）に集まります。彼らの逆方向のスピンは中程度ですが、コアの大きさに依存します。
  • 複雑なダンス（$\nu \ge 2$）: 「はぐれもの」のダンサーたちは、部屋の壁（端）の近くに集まります。彼らの逆方向のスピンは鋭く、顕著です。

4. 「逆流（カウンターフロー）」の驚き

最も興味深い発見の一つは、**逆流（counter-flows）**の存在です。
メインのダンスフロアが時計回りに回転していると想像してください。論文では、特定の複雑なシナリオにおいて、反時計回りに回転する小さなポケット（集団）が存在することを発見しました。

• 強力な渦の中心では、一部のダンサーが後ろ向きに回転しています。
• 壁の近くでは、他のダンサーが後ろ向きに回転しています。
  これら「後ろ向きに回転するポケット」こそが、先述の「ペアを失ったフェルミ粒子」です。彼らはブレーキとして機能し、システムの総スピンを減少させます。

まとめ

この論文は本質的に次のように述べています。

1. 単純さは予測可能: ダンスが単純で、ひねりが単純であれば、総スピンは計算通りになります。
2. 複雑さが混沌を生む: より多くのひねりを加えたり、ダンスをより複雑にしたりすると、「はぐれもの」のダンサーが現れます。
3. はぐれものがスピンを打ち消す: これらのペアになっていないダンサーたちは逆方向に回転し、システムの総スピンを減少させます。
4. 場所が重要: ひねりが強いか、あるいはダンスが複雑かによって、これらの「はぐれもの」は渦の中心に隠れるか、あるいは壁の近くに隠れます。

著者たちは新しい機械や医療用途を提案したわけではありません。彼らは単に、回転させたときにこれらの量子的なダンサーがどのように振る舞うかを正確に描き出し、「完璧な」スピンは非常に特定の、単純な条件下でのみ起こることを証明したのです。

技術概要

技術要約：二次元カイラル超流動体における渦

問題提起
本研究は、$(p_x + ip_y)^\nu$ 波のペアリング対称性を特徴とする二次元カイラル超流動体（SF）における軌道角運動量（OAM）、$L_z$ を調査するものである。具体的には、ゼロ温度下のディスク上において、軸対称な多重量子化渦（MQV）であり、渦度 $k$ を持つ系を対象としている。中心となる問題は「軌道角運動量パラドックス」であり、これは異なる理論的近似が、全OAMに対して相反する推定値を与えるというものである。素朴な図式では $L_z = \nu N/2$（ここで $N$ は全フェルミオン数、$\nu$ はペアリング次数）が示唆されるが、フェルミ液体理論および近年の Bogoliubov-de Gennes (BdG) 研究は、弱結合のバルディーン・クーパー・シュリーファー（BCS）領域において、スペクトル非対称性と非対のフェルミオンによる大幅な抑制が生じることを示唆している。本論文は、カイラルなペアリングと渦の循環の共存が、全OAMおよびその空間分布にどのように影響するかを明らかにすることを目的としており、特に、ペアリング対称性（$\nu=1$ の $p+ip$、$\nu=2$の$d+id$）が異なる条件下での、単一量子化渦（$|k|=1$）とMQV（$|k|>1$）の区別を行う。

手法
著者らは、無限円形ウェル（鏡面壁を持つディスク）内に閉じ込められた二次元カイラル・フェルミオン超流動体の枠組みとして、Bogoliubov-de Gennes (BdG) 理論を用いている。

• ハミルトニアン: 系は、対称化されたペアリング演算子 $\hat{\Delta} \sim (\hat{p}_x + i\hat{p}_y)^\nu$ を持つBCS平均場ハミルトニアンによって記述される。
• 渦の構成: 渦度 $k$ を持つ渦をディスクの中心に配置する。ギャップ関数は $\Delta(r) = \Delta_0 \tanh(r/\xi)e^{ik\theta}$ とモデル化される（$\xi$ はコヒーレンス長）。
• 数値解法: ベッセル関数の基底（円形ウェルにおける運動エネルギー演算子の固有関数）を用いて、準粒子波動関数（$u_n, v_n$）を展開することで、BdG方程式を解く。これにより、問題は行列固有値問題へと帰着される。
• 観測量: 著者らは、粒子密度 $n(r)$、方位角質量流、および局所的なOAM分布 $L_z(r)$ を計算する。全OAMは、基底状態の期待値から導出される。
• スペクトル非対称性: 主要な解析ツールとして、角運動量チャネル $l$ ごとの正のエネルギー固有値と負のエネルギー固有値の差として定義されるスペクトルフロー（またはスペクトル非対称性）$\eta_l$ を用いる。著者らは、全OAMの「完全な」値からの偏差を求めるために、一般化されたBogoliubov変換を利用し、これをスペクトル非対称性と関連付けている。

主な貢献と結果
本研究は、ペアリング次数 $\nu$、渦の渦度 $k$、および結合領域（BEC vs BCS）に応じて、OAMがどのように異なる挙動を示すかを明らかにしている。

1. BEC領域: ボース・アインシュタイン凝縮（BEC）領域（$\mu/E_F < 0$）では、エネルギー・スペクトルは完全にギャップが開いており、インギャップ状態が存在しない。その結果、すべての $l$ に対してスペクトルフロー $\eta_l$ はゼロとなる。全OAMは「完全な」値に従う：
   $$L_z = \frac{(k + \nu)N}{2}$$
   これは任意の整数 $\nu$ および $k$ に対して成立し、すべてのクーパー対が $k+\nu$ の角運動量を寄与することを示している。

2. **BCS領域 ($p+ip$かつ$k=\pm 1$):** 単一量子化渦（$k=\pm 1$）を持つカイラル$p+ip$超流動体（$\nu=1$）の場合、スペクトルフローは依然としてゼロである。全OAMは完全な値 $L_z = (k+\nu)N/2$ を保持する。

   • $k=1$ のとき、$L_z = N$ である。
   • $k=-1$（反渦）のとき、$L_z = 0$ である。全OAMはゼロであるが、空間分布 $L_z(r)$ は非自明である。すなわち、反渦コアからの負の寄与が、$p+ip$ ペアリングによる正のエッジ電流の寄与と正確に相殺される。

3. BCS領域 ($\nu \ge 2$ または $|k| \ge 2$): ペアリング次数が $p$ 波より高い場合（$\nu \ge 2$）、あるいは渦度が1より大きい場合（$|k| \ge 2$）には、顕著な偏差が生じる。

   • スペクトルフロー: これらの場合、インギャップ状態（エッジモード）が現れ、非ゼロのスペクトルフロー $\eta_l$ を引き起こす。
   • OAMの抑制: 全OAMは「完全な」値から著しく減少する。この抑制は、基底状態における非対のフェルミオンが、クーパー対とは逆向きのOAMを持つことに起因する。
   • 抑制のメカニズム:
     • $\nu \ge 2$（例：$d+id$で$k=1$）の場合、抑制は系の境界付近のエッジモードに関連している。
     • $|k| \ge 2$（例：$p+ip$で$k=2$）の場合、抑制は中程度であり、コア径に依存する。これはインギャップ・ボルテックスモードに起因する。
   • カウンターフロー: 非対フェルミオンの存在は、インギャップ・ボルテックスモード（$|k| \ge 2$ の場合）または境界付近（$\nu \ge 2$ の場合）における「カウンターフロー（負のOAM電流）」として現れる。これらのカウンターフローが、全OAMの減少を説明している。

意義と主張
本論文は、カイラルなペアリングにおける渦誘起の循環と、マクロなOAMを決定する相互作用を解明したと主張している。著者らは、「完全な」OAMの値は、システムにスペクトル非対称性（すなわち、フェルミ準位を横切るインギャップ・エッジモード）が存在しない場合にのみ堅牢であることを示している。

• BCS領域におけるOAMの減少は、スペクトルフローおよびBdGハミルトニアンの基底状態における非対フェルミオンの存在と根本的に結びついていることを確立した。
• 本研究は、高次のペアリング（$\nu \ge 2$）と高渦度（$|k| \ge 2$）の両方のシナリオにおけるOAM抑制について、共通のメカニズム、すなわちスペクトル非対称性が非対フェルミオンを引き起こすという点に注目し、統一的な説明を提供している。
• 著者らは、新しい実験装置や将来の応用を提案するのではなく、カイラル超流動体における渦のトポロジカルな性質（スペクトルフロー、エッジモード）に関する理論的な明確化を提供しており、これは軌道角運動量パラドックスに関するこれまでの知見を支持するものである。