A charge-flow instability in plasmas with charge fluctuations

技術要約：磁化プラズマにおける電荷流不安定性（Charge-Flow Instability）

問題提起
標準的な磁気流体力学（MHD）およびそのカイラル拡張は、通常、準中性条件または急速な電荷遮蔽を仮定しており、バルク流における電荷揺らぎ（$\mu \neq 0$）の動的な効果を無視している。カイラルMHDの文献では、カイラル化学ポテンシャル $\mu_5$ によって駆動されるカイラル磁気効果（CME）やカイラル渦効果（CVE）が広く扱われているが、標準的な化学ポテンシャル $\mu$（フェルミオンと反フェルミオンの非対称性を定量化するもの）とバルクプラズマ速度 $\mathbf{v}$ との相互作用については、体系的な分析が行われてこなかった。具体的には、準中性からの逸脱が生じた際に発生する、$\mu \mathbf{v}$ に比例する電流の存在が、安定性解析から漏れていた。本論文では、この「電荷流（C-flow）」項が、既知の静電不安定性（ダイオクトロン不安定性など）とは異なる、磁化プラスマにおける線形不安定性を誘発するかどうかを調査する。

手法
著者らは、電荷を持つフェルミオンのプラズマに対する一流体近似に基づき、彼らの付随する研究 [29] から導出された完全なカイラルMHD方程式を用いた理論的および数値的なアプローチを採用している。

• 理論的導出: 著者らは、静的で均質な背景磁場 $\mathbf{B}_0$ とゼロの背景速度を持つ系に対して、支配方程式を線形化した。小さな化学ポテンシャル（$\mu, \mu_5$）および亜相対論的な速度を仮定している。解析では、一般化されたオームの法則（式5）における、$\mu \mathbf{v}$ に比例する起電力（EMF）として機能するC-flow項を孤立させている。
• 分散解析: 並行伝播（$\mathbf{k} \parallel \mathbf{B}_0$）を仮定し、C-flowの効果を分離するためにカイラル化学ポテンシャルを $\mu_5 = 0$ と設定することで、分散関係（式14）を導出した。著者らは、成長率を決定するために固有モードを解析し、不安定性の条件を特定した。
• 数値シミュレーション: 線形理論の検証および非線形進化の探索のため、Pencil Codeを用いた直接数値シミュレーションを実施した。著者らは、CMEの寄与を除去するために $\mu_5 = 0$ とした簡約化されたカイラルMHD方程式（式2–7）を解いている。シミュレーションは、速度場と磁場に対するガウス型ランダム初期場、および $\mu$ に対する正弦波状の初期プロファイルを用いて行われた。研究には、均質な $\mu$ を用いた実行と、平均がゼロで長距離の揺らぎを持つ $\mu$ を用いた実行が含まれている。

主要な貢献

1. 新しい不安定性の特定: 本論文は、C-flow電流（$\propto \mu \mathbf{v}$）によって駆動される線形不安定性を特定し、これをC-flow不安定性と命名した。このメカニズムは、$\mu_5$ を必要とせず $\mu_5 = 0$ でも動作するため、カイラルプラズマ不安定性とは明確に異なる。
2. 解析的分散関係: 著者らは分散関係（式14）と、この不安定性の最大成長率を導出した。C-flow項が、非エルミート線形演算子の特徴である、分散関係の平方根における虚数成分から生じることを示している。
3. スケーリング則: 小さな背景磁場（$B_0$）および単位磁気プラントル数（$Pm = 1$）において、最大成長率は$|\omega_i|{\text{max}} = \frac{C{\text{flow}} k_\mu |\mu| B_0}{8}$ と導出される。注目すべきは、C-flow項が誘導方程式に入るために $\eta$ が必要であるにもかかわらず、この成長率が磁気拡散率 $\eta$ に依存しないことである。
4. 数値的確認: 直接数値シミュレーションにより、理論的な成長率とスケーリング則が、自由パラメータなしで高い精度（数パーセント以内の誤差）で一致することが確認された。

結果

• 線形領域: シミュレーションは、並行伝播において系が線形に不安定であることを確認している。支配的なモードは特定の波数 $k_{\text{max}}$ に落ち着き、解析的な予測と一致する一定の割合で成長する。成長率は、C-flow係数（$C_{\text{flow}}$）、背景磁場（$B_0$）、および化学ポテンシャルの振幅（$\mu$）に対して線形にスケールする。
• ゼロ平均の揺らぎ: 初期化学ポテンシャル $\mu$ の体積平均がゼロであっても（Run M）、この不安定性は持続する。このシナリオでは、線形段階において平均 $\mu$ は一定に保たれるが、不安定性が磁場と速度場の成長を駆動する。
• 非線形進化と暴走: カイラルダイナモが、全カイラリティと磁気ヘリシティの保存（$\mu_5$ から磁気ヘリシティへの変換）によって飽和するのに対し、C-flow不安定性は一流体近似において $\mu$ の枯渇メカニズムを欠いている。その結果、摂動が背景場と同程度になると、磁場と速度場が有限時間内に爆発する「暴走」挙動を示す。
• 飽和の不在: 著者らは、一流体モデルにおいて飽和が見られないことを観察している。これは、動的な電場（変位電流）の無視と準中性条件の仮定に起因すると考えている。彼らは、電荷分離の静電エネルギーとランダウ振動を考慮した適切な二流体記述を用いることが、飽和を記述するために必要であると示唆している。

意義
本論文は、C-flow不安定性が、標準的およびカイラルMHDにおいてこれまで考慮されてこなかったメカニズムであり、平均ゼロの電荷揺らぎによって引き起こされ得ることを主張している。この不安定性に関連する駆動電流は、CVE電流を超えることがあり、特定のシナリオではCME電流に匹敵することもある。したがって、著者らは、凝縮系物理学から初期宇宙の磁気生成、原始中性子星に至るまで、$\mu \neq 0$ であるプラズマの包括的な解析において、この項を含める必要があると論じている。本研究は、カイラル不安定性と電荷流不安定性の根本的な違いを浮き彫りにしている。前者はカイラリティ保存を通じて飽和するが、後者は現在の近似においては有限時間の特異点へと至る。このことは、飽和の物理が、より複雑な二流体記述を必要とする縦方向の電子的応答によって支配されていることを示唆している。