[논문 리뷰] Anisotropic transport and early dynamical impact of Cosmic Rays around Supernova remnants

Philipp Girichidis, Thorsten Naab외 2명|arXiv (Cornell University)|2014. 06. 18.

Solar and Space Plasma Dynamics인용 수 4

한 줄 요약

이 논문은 코스믹 레이즈(CRs)를 별도의 유체로 모델링한 새로운 구현 방식을 제시한다. 이 방식은 에너지에 따라 변화하는 비대칭 확산, 유동, 압축 손실을 고려하며, 자기장에 따라 수직 및 수평으로 방향성 있는 확산을 모델링한다. 연구 결과, CRs는 주로 자기장에 수직인 방향으로 열값의 약 100배에 이르는 강력하고 비대칭적인 가속도를 유도하며, 1 GeV에서 피크를 이룬다. 이는 80 pc 규모에서 약 150 kyr 동안의 시뮬레이션에서, 거듭제곱 법칙 스펙트럼보다 더 큰 역학적 영향을 미친다.

ABSTRACT

We present a novel implementation of cosmic rays (CR) in the magneto-hydrodynamic code FLASH. CRs are described as separate fluids with different energies. CR advection, energy dependent anisotropic diffusion with respect to the magnetic field and adiabatic losses to follow the evolution of spectra are taken into account. We present a first study of the transport and immediate (~150 kyr) dynamical impact of CRs on the turbulent magnetised interstellar medium around supernova remnants on scales up to 80 pc. CR diffusion quickly leads to an efficient acceleration of low-density gas (mainly perpendicular to the magnetic field) with accelerations up to two orders of magnitude above the thermal values. Peaked (at 1 GeV) CR injection spectra have a stronger impact on the dynamics than power-law spectra. For realistic magnetic field configurations low energy CRs (with smaller diffusion coefficients) distribute anisotropically with large spatial variations of a factor of ten and more. Adiabatic losses can change the local spectra perceptibly but do not have an integral effect on the dynamics at the spatial and temporal scales considered here. We discuss the potential global impact of CRs and anisotropic transport on the dynamical structure of the ISM and also detail the limitations of the model.

연구 동기 및 목표

• 자기장역학유체역학(MHD) 코드에 물리적으로 일관된 코스믹 레이즈(CRs) 모델을 개발하고 구현하여, 에너지에 따라 변화하는 운반 특성을 고려한다.
• 초기 역학적 영향(약 150 kyr 이내)을 분석한다. 특히 초신성잔해(SNRs) 주변의 난류적이고 자기장이 있는 은하간매질(ISM)에서 80 pc 이내 스케일에서 CRs의 영향을 연구한다.
• CR 확산의 비대칭성, 스펙트럼 형태(피크형 대비 거듭제곱 법칙형), 그리고 압축 손실이 CR 압력 기울기와 기체 가속도에 미치는 영향을 평가한다.
• 스펙트럼 진화 과정에서 CR 에너지 운반의 수치적 안정성, 특히 압축 손실의 가역성과 총 에너지 보존 여부를 평가한다.

제안 방법

• 에너지에 따라 변화하는 운반 특성을 고려하여, 다중 밴드 에너지 분포를 사용해 스펙트럼 진화를 해상도를 높여 모델링한 별도의 유체로 CRs를 구현한다.
• 유체 방정식을 통해 CR 유동을 모델링하고, 자기장 선과 수직 방향으로 별도의 계수를 가진 비대칭 확산을 도입하며, 입자 에너지에 따라 달라진다.
• 에너지 밴드 간의 유량 기반 방법을 사용해 압축 손실을 포함하며, 기계 정밀도 수준에서 에너지 보존을 보장한다.
• 유한체적 MHD 해소기(FLASH)를 사용하고 자기장 진화를 자동으로 일관되게 처리하여, 난류적 ISM에서 SNR 팽창과 CR 피드백을 시뮬레이션한다.
• CR 스펙트럼에 대해 조각별로 일정한 분포 함수를 적용하고, 압축 흐름 하에서 압축 손실의 가역성을 시험한다.
• 단일 세포 초신성 폭발 및 정현파 압축을 포함한 통제된 수치 시험을 수행하여 에너지 보존 및 스펙트럼 진화의 정당성을 검증한다.

실험 결과

연구 질문

• RQ1 자기장에 따라 강하게 방향성 있는 확산(자기장 선 방향으로는 강하고 수직 방향으로는 약한)이 CRs와 주변 ISM 간의 역학적 결합에 어떤 영향을 미치는가?
• RQ2 초기 단계(약 150 kyr 이내)에서 1 GeV에서 피크를 이룬 스펙트럼과 거듭제곱 법칙 스펙트럼 간에, 기체 가속도에 대한 CRs의 상대적 역학적 영향은 어떠한가?
• RQ3 압축 손실이 국소적 CR 스펙트럼을 얼마나 변화시키며, 80 pc 규모에서 시간이 지남에 따라 ISM로의 전체 에너지 전달에 어떤 영향을 미치는가?
• RQ4 스펙트럼 해상도와 에너지 유량의 가역성에 대한 수치적 제약이 CR 에너지 운반 및 스펙트럼 진화의 정확도에 어떤 영향을 미치는가?
• RQ5 CR 압력 기울기가 초신성잔해 충격면 근처의 저밀도 기체 영역에서 초기 비열역학적 가속도를 유도하는 데 어떤 역할을 하는가?

주요 결과

• 비대칭 CR 확산은 주로 자기장 선에 수직인 방향으로 열값보다 약 두 개의 지수 차수 높은 기체 가속도를 유도한다.
• 1 GeV에서 피크를 이룬 CR의 주입 스펙트럼은 거듭제곱 법칙 스펙트럼보다 더 강력한 ISM에 대한 역학적 영향을 미친다. 이는 더 강한 압력 기울기 때문이며.
• 저에너지 CRs는 매우 비대칭적인 공간 분포를 보이며, 그 변화 폭이 10배 이상에 이르며, 이는 그들의 작은 확산 계수에 기인한다.
• 압축 손실은 국소적 CR 스펙트럼을 눈에 띄게 변화시키지만, 연구된 공간 및 시간 스케일에서 ISM 전체에 대한 CR의 전반적 역학적 영향에는 크게 영향을 주지 않는다.
• 수치적 구현은 총 CR 에너지를 기계 정밀도 수준에서 보존하지만, 10⁻¹² s⁻¹ 이상의 압축 흐름 속도에서는 압축 손실의 가역성이 붕괴되며, 이는 실제로 SNR 충격에서 관측되는 물리적 범위에 가깝다.
• 단일 세포 초신성 주입 시뮬레이션 결과, 내부 충격 구조를 해상도로 표현하지 않으면 스펙트럼 경화 현상이 발생하지 않으며, 이는 정확한 CR 스펙트럼 진화를 위해 충격 구조를 해상도로 표현하는 것이 중요함을 시사한다.