[論文レビュー] Diffusion of relativistic charged particles and field lines in isotropic turbulence: II. Analytical models

Marco Kuhlen, Vo Hong Minh Phan+1 more|arXiv (Cornell University)|Nov 10, 2022

Solar and Space Plasma Dynamics被引用数 4

ひとこと要約

この論文は、等方的乱流における垂直方向宇宙線輸送に関して、長年の理論モデルとシミュレーションの不一致を解消するため、3段階の解析的モデルを導入した。初期粒子輸送（磁場線に沿った輸送）、磁場線拡散、粒子の磁場線からの相関消失の段階を含む。非標準的な剛性スケーリングの原因は、磁場線輸送の準拡散的挙動に起因し、低剛性領域（標準理論が失敗する領域）を含め、広い剛性範囲でシミュレーションとモデルの一致が良好に得られた。

ABSTRACT

The transport of high-energy particles in the presence of small-scale, turbulent magnetic fields is a long-standing issue in astrophysics. Analytical theories on transport perpendicular to the large-scale magnetic field disagree with numerical simulations at rigidities where the particles' gyroradii are slightly smaller than the correlation length of turbulence. At the same time, extending the numerical simulations to lower rigidities has proven computationally prohibitive. We present an analytical model for the perpendicular transport, based on (1) initial particle transport along field lines, (2) the transport of field lines and (3) the eventual decorrelation of particles from field lines. Transport parallel to the large-scale field is governed by pitch-angle scattering and so for times larger than the inverse pitch-angle diffusion coefficient, particles spatially diffuse in the parallel direction. Our results suggest that perpendicular diffusion occurs when particles have displaced in the perpendicular direction by a few correlation lengths of turbulence. We have tested the analytical theory by running a large suite of test particle simulations at unprecedentedly low rigidities, making extensive use of graphical processing units (GPUs). Our numerical results exhibit a non-standard rigidity-dependence for the perpendicular diffusion coefficient at intermediate rigidities. At the lowest rigidities, the standard rigidity-dependence is recovered. The simulated diffusion coefficients are nicely reproduced by our analytical model. We have traced the non-standard rigidity-dependence to a subdiffusive phase in the field line transport. Our study confirms our understanding of the escape of cosmic rays from the Galactic halo and its rigidity-dependence. [abridged]

研究の動機と目的

• 中程度の剛性領域における垂直方向宇宙線輸送に関して、解析的モデルと数値的シミュレーションの間の継続的な不一致を解消すること。
• 標準理論が予測しないが、シミュレーションで観測される垂直拡散係数における非標準的剛性スケーリングの原因を特定すること。
• 広い剛性範囲（特に低剛性領域）において、シミュレーション結果を正確に再現できる解析的モデルの構築。
• 磁場線輸送と宇宙線の垂直拡散の物理的関係、特に準拡散的磁場線運動の役割を明確化すること。
• 太陽系および銀河間媒体に代表される乱れた宇宙プラズマにおける宇宙線輸送を物理的に一貫した枠組みでモデル化すること。

提案手法

• GPU加速計算を用いた大規模テスト粒子シミュレーションにより、前例のない低剛性領域（r_g/L_c ≲ 10⁻³）に到達し、従来のシミュレーション限界を超えた。
• 垂直拡散係数を3段階のプロセスとしてモデル化する：(1) 磁場線に沿った初期球面的輸送、(2) 乱流磁場構造に支配される磁場線拡散、(3) 粒子が磁場線から最終的に相関を失う過程。
• コルシンの仮定と等方的コルゴロフ乱流を前提に、σ_x²(z) と d_FL(z) の連立常微分方程式を半解析的に解き、走査的磁場線拡散係数をモデル化する。
• 外スケール L_max と相関長 L_c ≈ L_max/5 を持つコルゴロフ k⁻⁵/³ スペクトルを乱流磁場に適用する。
• 粒子速度がアルベール速度を超えるものと仮定し、静磁場近似を適用し、乱流磁場実現と粒子初期方向のアンサンブル平均を取る。
• 異なる乱流度（η = 0.2, 0.5）と剛性において、走査的および漸近的拡散係数を比較することで、解析的モデルの妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 なぜ標準的な解析的モデルは、中程度の剛性領域で数値的シミュレーションで観測される剛性依存的垂直拡散係数を再現できないのか？
• RQ2 等方的乱流における垂直拡散係数の非標準的剛性スケーリングの背後にある物理的メカニズムは何か？
• RQ3 磁場線輸送は時間とともにどのように変化するのか？また、異常な粒子拡散を説明できるような準拡散的挙動を示すのか？
• RQ4 粒子の垂直拡散係数を、磁場線拡散と粒子の相関消失に分解することで、どの程度正確にモデル化できるのか？
• RQ5 磁場線ダイナミクスに基づく解析的モデルは、標準理論が適用できない低剛性領域を含め、全剛性範囲でシミュレーション結果を再現できるか？

主な発見

• シミュレーションされた垂直拡散係数は、中程度の剛性領域（r_g/L_c ≲ 1）で非標準的な剛性依存性を示し、λ_⊥/λ_∥ は高剛性領域に達するまで増加し、その後急激に低下する。
• 最低剛性領域（r_g/L_c ≲ 10⁻³）では、標準的な剛性スケーリング（λ_⊥ ∝ r_g²）が回復し、理論的期待と整合的である。
• 非標準的剛性依存性は、磁場線輸送における準拡散的段階に起因し、d_FL(z) が初期に線形より速く増加し、過剰応答と振動を経て、定常値に到達する。
• 磁場線常微分方程式と粒子の相関消失に基づく解析的モデルは、全剛性範囲および全乱流度において、シミュレーションされた拡散係数を高い精度で再現した。
• 常微分方程式モデルは、積分における干渉項（cos(k cosθ z)）に起因する、走査的磁場線拡散係数の振動的過剰応答を捉え、球面的から拡散的磁場線運動への遷移を再現した。
• わずかな周期・位相差のずれを補正するための小さなスケーリング係数（η=0.2 では 0.8、η=0.5 では 0.56）を適用することで、シミュレーションと解析的モデルの一致が向上し、モデルの頑健性が示された。