[논문 리뷰] Three dimensional modeling of CR propagation

Daniele Gaggero, Luca Maccione외 3명|arXiv (Cornell University)|2013. 06. 28.

Astrophysics and Cosmic Phenomena인용 수 3

한 줄 요약

이 논문은 은하계에서 태양계에서의 우주선(CR) 전파를 위한 DRAGON 코드에 중대한 업그레이드를 도입하여, 정규 자기장에 수직 및 평행한 방향으로 위치 및 레지드리티(에너지)에 따라 변화하는 확산 계수를 고려한 완전한 3차원 이방성 확산 모델링을 제시한다. 이는 처음으로 태양의 위치가 간성 사이 지역에 있음과 더불어 에너지 손실이 증가함으로써 전자 스펙트럼이 자연스럽게 기울어지며, 충격 가속 이론과 호환되는 2.38의 주요 주입 지수를 사용하여 AMS-02 및 PAMELA 데이터와 일관된 피팅이 가능하다는 것을 보여준다.

ABSTRACT

We present here a major upgrade of DRAGON, a numerical package that computes the propagation of a wide set of CR species from both astrophysical and exotic origin in the Galaxy in a wide energy range from tens of MeV to tens of TeV. DRAGON takes into account all relevant processes in particular diffusion, convection, reacceleration, fragmentation and energy losses. For the first time, we present a full 3D version of DRAGON with anisotropic position-dependent diffusion. In this version, the propagation is calculated within a 3D cartesian grid and the user is able to implement realistic and structured three dimensional source, gas and regular magnetic field distributions. Moreover, it is possible to specify an arbitrary function of position and rigidity for the diffusion coefficients in the parallel and perpendicular direction to the regular magnetic field of the Galaxy. The code opens many new possibilities in the study of CR physics. In particular, we can study for the first time the impact of the spiral arm structure on the leptonic spectra: taking into account the fact that we live in an interarm region, far from most sources, we obtain - due to increased energy losses - a steeper electron spectrum compared to the assumption of a smooth source term. We discuss the implications of these results on our understanding on leptonic spectra and we briefly mention future studies that can be performed with our new 3D code.

연구 동기 및 목표

• 실제 은하계 구조와 이방성 확산을 고려한 3차원, 이방성 우주선 전파 모델을 개발하기 위해.
• 이전 모델들이 등방성 확산과 매끄러운 원천 분포를 가정하는 데서 비롯하는 한계를 해결하기 위해.
• 태양이 간성 사이 지역에 위치함으로써 증가된 에너지 손실이 전자 및 양성자 스펙트럼에 미치는 영향을 조사하기 위해.
• 물리적으로 타당한 주입 스펙트럼을 사용하여 AMS-02 및 PAMELA의 고정밀도 전자 및 양성자 데이터를 재현하기 위해.
• 향후 태양계 내 이방성 확산과 그가 태양계 및 국소 간성간 매질 내 우주선 운반에 미치는 영향을 연구할 수 있도록 하기 위해.

제안 방법

• 코드는 은하 자기장에 대해 수직 및 평행 방향으로 위치 및 레지드리티에 따라 변화하는 확산 계수를 가진 3차원 우주선 운반 방정식을 카르테시안 좌표계에서 해석한다.
• CR 원천, 간성간 가스, 정규 자기장에 대해 임의의 3차원 분포를 허용하여 나선형 패턴과 국소 구조를 현실적으로 모델링할 수 있다.
• 확산 텐서는 $ D_{\parallel} $, $ D_{\bot} $ 및 자기장 단위 벡터로 구성되며, 공간적 및 에너지적 의존성이 분리되어 있다.
• 운반 방정식은 원천, 대류, 재가속, 에너지 손실 및 분열을 포함하여 연산자 분할과 유한차분 방법을 사용해 수치적으로 해석한다.
• 태양 조절은 HelioProp 코드를 통해 포함되며, 전하에 의존하는 이동과 시간에 따라 변화하는 태양계 매개변수를 고려한다.
• 순수한 수직 및 수평 확산에 대해 그린 함수를 사용하여 모델의 수치 정확도를 검증한다.

실험 결과

연구 질문

• RQ1 태양이 간성 사이 지역에 위치함으로써 증가된 에너지 손실로 인해 예측된 전자 스펙트럼에 어떤 영향을 미치는가?
• RQ2 3차원 원천 구조와 에너지 손실을 적절히 모델링할 경우, $ \gamma_0 = 2.38 $의 단일 주요 전자 주입 스펙트럼이 AMS-02 및 PAMELA 데이터와 일치할 수 있는가?
• RQ3 이방성 확산이 국소 은하계 환경에서 고에너지 우주선 전자 스펙트럼을 형성하는 데 어떤 역할을 하는가?
• RQ4 실제 3차원 원천 분포를 고려할 때, 우주선의 관측된 양성자 분수 상승을 하드 스펙트럼의 추가 성분 없이 설명할 수 있는가?
• RQ5 국소 기포와 나선형 패턴의 과밀도와 같은 국소 구조는 고에너지 렙톤의 전파에 어떤 영향을 미치는가?

주요 결과

• 주어진 전자 주입 지수 $ \gamma_0 = 2.38 $를 사용하여 모델은 AMS-02의 양성자 분수와 PAMELA의 전자 데이터를 잘 재현하며, 충격 가속 이론과 일치한다.
• 관측된 고에너지 영역에서 전자 스펙트럼의 기울어짐은 하드 스펙트럼의 추가 성분이 필요 없이 태양이 간성 사이 지역에 위치해 있고 에너지 손실이 증가함으로써 자연스럽게 설명된다.
• 10 테라전자볼트(10 TeV)의 截止 에너지와 $ \gamma_0 = 1.7 $를 가진 제2차 성분은 양성자 분수에 좋은 피팅을 제공하며, 1 테라전자볼트(1 TeV)의 截止 에너지도 수용 가능한 피팅을 제공한다.
• 200 기로전자볼트(200 GeV) 이상의 페르미-LAT 전자 스펙트럼은 모델에서 과소 예측되며, 벨라 초신성 잔해와 같은 근접한 국소 원천의 기여 가능성을 시사한다.
• 단일 근접 전자 가속기에서 $ E_{\text{cut}} = 1~\text{TeV} $ 이며 $ \gamma = 2.1 $일 경우, 페르미-LAT 및 AMS-02 데이터에 대한 통합 피팅이 양호하게 이루어지며, 이는 국소 원천이 고에너지 초과를 설명할 수 있음을 시사한다.