[논문 리뷰] Diffusion of relativistic charged particles and field lines in isotropic turbulence: I. Numerical simulations

Marco Kuhlen, Vo Hong Minh Phan외 1명|arXiv (Cornell University)|2022. 11. 10.

Solar and Space Plasma Dynamics인용 수 4

한 줄 요약

이 연구는 등방성 콜모고로프 turbulence에서 상대론적 입자에 대한 고해상도 시험입자 시뮬레이션을 제시하며, 이는 수평 평균 자유경로 λ⊥이 중간 강성도에서 (rg/Lc)^0.5 비례함을 밝혀내지만, 매우 낮은 강성도에서는 (rg/Lc)^1/3 비례로 전이됨을 보여주며, λ∥는 모든 강성도에서 일관되게 (rg/Lc)^1/3 비례함을 확인한다. 이 비보편적 스케일링은 이전의 가정을 도전하며 은하 내 우주선 운반에 중요한 함의를 지닌다.

ABSTRACT

The transport of non-thermal particles across a large-scale magnetic field in the presence of magnetised turbulence has been a long-standing issue in high-energy astrophysics. Of particular interest is the dependence of the parallel and perpendicular mean free paths $λ_{\parallel}$ and $λ_{\perp}$ on rigidity $\mathcal{R}$. We have revisited this important issue with a view to applications from the transport of Galactic cosmic rays to acceleration at astrophysical shocks. We have run test particle simulations of cosmic ray transport in synthetic, isotropic Kolmogorov turbulence at unprecedentedly low reduced rigidities $ g/\Lc \simeq 10^{-4}$, corresponding to $\mathcal{R} \simeq 10 \, ext{TV}$ for a turbulent magnetic field of $\Brms = 4 \, μ ext{G}$ and correlation length $\Lc = 30 \, ext{pc}$. Extracting the (asymptotic) parallel and perpendicular mean free paths $λ_{\parallel}$ and $λ_{\perp}$, we have found $λ_{\parallel} \propto ( g/\Lc)^{1/3}$ as expected for a Kolmogorov turbulence spectrum. In contrast, $λ_{\perp}$ has a faster dependence on $ g/\Lc$ for $10^{-2} \lesssim g/\Lc \lesssim 1$, but for $ g/\Lc \ll 10^{-2}$, also $λ_{\perp} \propto ( g/\Lc)^{1/3}$. Our results have important implications for the transport of Galactic cosmic rays.

연구 동기 및 목표

• 등방성 난류에서 상대론적 우주선의 수평 평균 자유경로 λ⊥의 강성 의존성에 대한 오랫동안 지속된 불확실성을 해결하기 위해.
• 낮은 강성도에서 수평 확산 스케일링이 (λ∥ ∝ (rg/Lc)^1/3)와 다를 수 있는지 테스트하기 위해.
• 우주선 확산 또는 유동에 의한 탈출이 우세해지는 조건을 규명하기 위해.
• 난류 자기장에서 자기장선 및 입자 확산의 분석 모델에 대한 수치 기준을 제공하기 위해.

제안 방법

• 합성적 등방성 콜모고로프 난류에서 수치적 시험입자 시뮬레이션을 수행하였으며, 스펙트럼은 E(k) ∝ k^{-5/3}의 거듭제곱 법칙을 따른다.
• 시뮬레이션은 B_rms = 4 μG 및 Lc = 30 pc 조건에서 강성도 ≈10 TV인 매우 낮은 감소된 강성도 rg/Lc ≈ 10^{-4}를 사용하였다.
• 시간 평균 제곱 이동을 이용하여 입자 궤적에서 최종 상태의 수평 및 수직 평균 자유경로 λ∥ 및 λ⊥을 추출하였다.
• 시뮬레이션은 자기장선 따라의 입자 확산과 자기장선의 무작위 이동을 모두 고려하였으며, 후자는 수평 운반에 기여하였다.
• 결과는 준선형 이론 및 비선형 운반 모델과 비교되었으며, 특히 90도 문제와 공진 너비 확장에 초점을 맞추었다.
• 추가 분석을 통해 경쟁적 탈출 메커니즘을 평가하였으며, 이는 수평 확산, 속도 U로 인한 유동 탈출, 복사 냉각 손실이 포함되었다.

실험 결과

연구 질문

• RQ1 등방성 콜모고로프 난류에서 수평 평균 자유경로 λ⊥은 감소된 강성도 rg/Lc에 따라 어떻게 스케일링되는가?
• RQ2 수평 확산 스케일링 λ⊥ ∝ (rg/Lc)^1/3은 보편적인가, 아니면 낮은 강성도에서 변화하는가?
• RQ3 다양한 강성도에서 우주선 탈출에 있어 수평 확산, 유동, 냉각 손실의 상대적 중요도는 무엇인가?
• RQ4 결과는 준선형 및 비선형 운반 이론의 예측과 어떻게 비교되는가?

주요 결과

• 모든 강성도에서 수평 평균 자유경로 λ∥는 (rg/Lc)^1/3 비례로 스케일링되며, 콜모고로프 난류에 대한 이론적 예측과 일치한다.
• 중간 강성도에서 (rg/Lc ≲ 1) 수평 자유경로 λ⊥은 (rg/Lc)^0.5 비례로 스케일링되며, 단순 스케일링보다 더 강한 수평 확산을 나타낸다.
• 매우 낮은 강성도에서 (rg/Lc ≪ 1) 수평 자유경로 λ⊥은 다시 (rg/Lc)^1/3 비례로 전이되며, 수평 스케일링과 일치한다.
• 수평 및 수직 확산 계수의 비율 κ⊥/κ∥는 rg/Lc ≤ 3×10^{-3} 범위에서 약 2×10^{-2}로 일정하게 유지되며, 넓은 에너지 범위에서 중요한 수평 운반을 나타낸다.
• λ⊥ ∝ (rg/Lc)^1/3를 가정할 경우, rg/Lc ≈ 3×10^{-6} 이하의 강성도에서는 유동이 확산보다 우세하지만, λ⊥ ∝ (rg/Lc)^1/2일 경우 이 임계값은 약 3×10^{-5}로 증가한다.
• 고강성도에서는 냉각 손실이 확산 손실보다 우세하며, λ⊥ ∝ (rg/Lc)^1/3일 경우 전이점은 rg/Lc ≈ 3.1×10^{-10}에서 발생한다. 이는 (rg/Lc)^1/2 스케일링보다 약 200배 낮은 값이다.