[논문 리뷰] Computation of diffusive shock acceleration using stochastic differential equations

L. Grémillet, J. G. Kirk|arXiv (Cornell University)|1999. 05. 14.

Cosmology and Gravitation Theories참고 문헌 1인용 수 6

한 줄 요약

이 논문은 천체 플라즈마에서 확산 충격 가속(DSA)을 모델링하기 위해 확률적 미분 방정식(SDE) 방법을 제시한다. 이는 비상대론적 충격파에서 입자 운반을 정확하게 시뮬레이션할 수 있도록 유체 속도 기울기의 암묵적 계산을 가능하게 한다. 이 방법은 평탄하거나 역전된 싱크로트론 스펙트럼을 성공적으로 재현하며, 특히 충격 간 거리가 10–100개의 확산 길이에 해당하는 시스템에서 가속-손실 평형으로 인한 피크 형성 특징까지 재현한다.

ABSTRACT

The present work considers diffusive shock acceleration at non-relativistic shocks using a system of stochastic differential equations (SDE) equivalent to the Fokker-Planck equation. We compute approximate solutions of the transport of cosmic particles at shock fronts with a SDE numerical scheme. The momentum gain is given by implicit calculations of the fluid velocity gradients using a linear interpolation between two consecutive time steps. We validate our procedure in the case of single shock acceleration with different shock thickness, with or without synchrotron losses. A comparative discussion of implicit and explicit schemes for different shock thickness shows that implicit calculations extend the range of applicability of SDE schemes to infinitely thin 1D shocks. The method is then applied to multiple shock acceleration for a system of identical shocks which free parameters are the distance between two consecutive shocks, the synchrotron losses time and the escape time of the particles. The stationary distribution reproduces quite well the flat differential logarithm energy distribution produced by multiple shock effect, and also the piling-up effect due synchrotron losses at a momentum where they equilibrate the acceleration rate. At higher momenta particle losses dominate and the spectrum drops. The competition between acceleration and loss effects leads to a pile-up shaped distribution which appears to be effective only in a restrict range of inter-shock distances of 10-100 diffusion lengths. We finally compute the optically thin synchrotron spectrum produced such periodic pattern which can explain flat and/or inverted spectra observed in Flat Radio spectrum Quasars and in the galactic centre.

연구 동기 및 목표

• 확률적 미분 방정식을 사용한 확산 충격 가속을 시뮬레이션하기 위한 강력한 수치적 방법을 개발하기 위해.
• 유체 속도 기울기의 암묵적 시간 적분을 통해 SDE 기법의 적용 범위를 임의로 두꺼운 충격파에까지 확장하기 위해.
• 다중 충격, 싱크로트론 손실, 입자 탈출의 병합 효과가 고에너지 우주선 스펙트럼에 미치는 영향을 조사하기 위해.
• 편평하거나 역전된 라디오 스펙트럼을 관측한 천체 물리적 천체, 예를 들어 편평한 스펙트럼 라디오 퀘이사와 은하 중심에서 관측된 스펙트럼을 재현하기 위해.
• 특히 충격 간 거리에 해당하는 매개변수 영역을 규명하여, 가속-손실 평형에 의해 발생하는 스펙트럼 피크 형성이 관측 가능한 조건을 도출하기 위해.

제안 방법

• 입자 운동량과 충격에서의 공간 확산을 기술하는 포커-플랑크 운반 방정식과 등가인 확률적 미분 방정식 시스템을 수립한다.
• 연속적인 시간 단계 사이의 선형 보간을 통해 유체 속도 기울기의 암묵적 시간 적분을 사용하여 얇은 충격파에 대한 해를 안정화시킨다.
• 변동 가능한 충격 두께와 공간적 간격을 가진 단일 및 다중 충격 구성을 대상으로 SDE 기법을 적용한다.
• 자기장과 입자 에너지에 기반한 운동량에 의존하는 에너지 손실률을 도입하여 싱크로트론 손실을 포함한다.
• 시뮬레이션 프레임워크 내에서 운동량과 시간에 의존하는 탈출 시간 상수를 통해 입자 탈출을 모델링한다.
• 관측 결과와 비교하기 위해 최종 입자 에너지 분포와 해당되는 광학적으로 얇은 싱크로트론 복사 스펙트럼을 계산한다.

실험 결과

연구 질문

• RQ1 암묵적 SDE 기법이 임의로 얇은 충격파에서 확산 충격 가속을 정확하게 모델링할 수 있는가?
• RQ2 싱크로트론 손실과 입자 탈출은 다중 충격 시스템에서 고에너지 우주선 스펙트럼의 형태에 어떤 영향을 미치는가?
• RQ3 어느 정도의 충격 간 거리 범위에서 가속-손실 평형에 의해 발생하는 관측 가능한 스펙트럼 피크 형성이 발생하는가?
• RQ4 SDE 방법이 편평한 스펙트럼 라디오 퀘이사와 은하 중심에서 관측된 편평하거나 역전된 라디오 스펙트럼을 재현할 수 있는가?
• RQ5 충격 두께와 손실 시간 상수가 최종 입자 분포 함수 형성에 어떤 역할을 하는가?

주요 결과

• 암묵적 SDE 기법은 무한히 얇은 1차원 충격파에 대해서도 안정적이고 정확한 계산을 가능하게 하여 SDE 기법의 적용 가능 범위를 확장한다.
• 정적인 입자 에너지 분포는 싱크로트론 손실이 가속률과 균형을 이루는 고운동량 영역에서 특징적인 피크 형성을 보이며, 더 높은 에너지에서는 주로 손실이 지배적이므로 급격한 감소가 나타난다.
• 피크 형성 특징은 충격 간 거리가 10에서 100개의 확산 길이에 해당하는 제한된 범위 내에서만 관측 가능하며, 이를 초월하면 효과가 약해진다.
• 계산된 광학적으로 얇은 싱크로트론 복사 스펙트럼은 편평하거나 역전된 스펙트럼을 재현하며, 이는 다중 충격 모델이 타당한 설명임을 지지한다.
• 이 방법은 가속, 싱크로트론 손실, 탈출 간의 경쟁을 성공적으로 포착하였으며, 손실 시간 상수와 충격 간격이 관측된 스펙트럼 특징 형성에 핵심적인 역할을 한다는 것을 보여준다.
• 단일 충격에서 싱크로트론 손실 유무에 따라 분석적 해와의 검증을 통해 SDE 수치적 접근의 정확성과 신뢰성을 확인하였다.