[논문 리뷰] Diffusion of relativistic charged particles and field lines in isotropic turbulence: II. Analytical models

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

Solar and Space Plasma Dynamics인용 수 4

한 줄 요약

이 논문은 등방성 난류에서 수직 우주선 운반에 대한 분석 모델과 시뮬레이션 간 오랫동안 지속된 괴리 문제를 해결한다. 세 단계 분석 모델(초기 입자 운반, 필드선 확산, 입자 비상관)을 도입하여, 수직 확산 계수의 비표준 강성 스케일링의 근본 원인으로 하위확산 필드선 운반을 규명한다. 시뮬레이션과 모델 간의 뛰어난 일치는 강성 전역에서 관찰되며, 특히 기존 이론이 실패하는 저강성 영역에서도 성립한다.

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⁻³)에 도달한다.
• 수직 확산 계수를 세 단계 과정으로 모델링한다: (1) 필드선 沿 운반, (2) 난류 자기장 구조에 의해 지배되는 필드선 확산, (3) 입자가 필드선에서 최종적으로 비상관됨.
• Corrsin의 가정과 등방성 Kolmogorov 난류 조건 하에서, σ_x²(z)와 d_FL(z)의 연립 미분방정식을 풀이하는 반분석적 ODE 접근법을 사용하여 런닝 필드선 확산 계수를 모델링한다.
• 외부 스케일 L_max와 상관 길이 L_c ≈ L_max/5를 갖는 Kolmogorov 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)가 초기에 선형보다 빠르게 증가하다가 최고점에 도달하고 진동을 반복한 후 점차 수렴하는 것으로 규명된다.
• 필드선 ODE 기반 모델과 입자 비상관을 통합한 분석 모델은 모든 강성과 난류 수준에서 시뮬레이션된 확산 계수를 매우 정확하게 재현한다.
• ODE 모델은 적분 내 간섭 항(cos(k cosθ z))으로 인해 발생하는 런닝 필드선 확산 계수의 진동적 초과반응을 포착하며, 구동성에서 확산성으로의 전이를 재현한다.
• 소규모 조정 인자(η=0.2일 경우 0.8, η=0.5일 경우 0.56)를 적용함으로써 시뮬레이션과 분석 모델 간의 일치가 향상되며, 이는 진동 위상과 주기의 약간의 불일치에도 불구하고 모델의 강건성을 시사한다.