다중 스케일 자기장과 대류가 만든 태양 디스크 감마선

초록

이 논문은 태양 광구와 최상부 대류층에서 대류가 만들어내는 다중 스케일 자기장 구조와 알프벤 파동 난류가 은하계 우주선(GCR)의 반사와 감마선 생성에 미치는 영향을 정량적으로 모델링한다. 큰 규모의 필라멘트형 자기장이 전체 감마선 스펙트럼을 형성하고, 30 km 이하 규모의 알프벤 난류가 100 GeV 이하에서 추가 억제를 일으킨다. 모델은 1 GeV–1 TeV 구간에서 관측된 Fermi‑LAT와 HAWC 스펙트럼 기울기를 성공적으로 재현하지만, 절대 플럭스는 관측치보다 2–5배 낮다.

상세 분석

본 연구는 기존의 두 가지 접근법—(1) 저층 대기에서의 정적 압력‑자기장 평형 모델, (2) 전지구적 PFSS와 같은 대규모 코르오날 모델—의 한계를 인식하고, 대류에 의해 형성되는 다중 스케일 자기장을 직접 시뮬레이션한다. 저층(‑200 km ~ 12 Mm) 영역을 세 부분으로 구분한다. 첫 번째는 수백 킬로미터 규모의 필라멘트형 네트워크 구조로, 플럭스 튜브와 플럭스 시트가 교차하며 ‘역와인병’ 형태의 확장‑수축을 보인다. 두 번째는 100 km 이하의 작은 스케일 필라멘트, 이는 대류 플럭스 압출에 의해 고강도(≈ kG) 영역이 간격을 두고 배열되는 현상이다. 세 번째는 대류 흐름이 생성하는 알프벤 파동이 상향·하향으로 전파되고, 반사된 파동과의 비선형 상호작용을 통해 30 km~1 km 규모까지 에너지를 cascade하는 난류이다. 연구팀은 RMHD(Reduced MHD) 방식을 이용해 7개의 개방 필드 라인에 대해 파동 전파·반사·감쇠를 3차원 MHD 시뮬레이션으로 구현하고, Grad‑Shafranov 방정식으로 평균 자기장을 정적 구조로 계산한 뒤, 두 결과를 중첩해 복합 자기장을 구성한다.

입자 추적에서는 전자기학적 로렌츠 힘을 완전하게 포함하고, GCR(주로 p와 He)의 입자 궤적을 Monte‑Carlo 방식으로 10⁶개 이상 시뮬레이션한다. 각 입자는 위에서 정의한 복합 자기장과 알프벤 난류에 의해 발생하는 피치각 산란을 경험한다. 특히, 난류 수준을 파라미터 δB/B₀로 조절해 0.1, 0.3, 1.0 등 다양한 경우를 탐색했으며, δB/B₀≈0.3(표준 난류)에서 100 GeV 이하 감마선 플럭스가 관측보다 약 3배 억제되는 것을 확인했다. 이는 알프벤 난류가 입자의 반사 확률을 감소시키고, 고에너지 입자는 더 깊은 층으로 침투해 효율적인 파이온 생성이 일어나게 함을 의미한다.

감마선 생성은 Kelner et al.(2006)와 Kafexhiu et al.(2014)의 최신 파이온·중성 파이온 생산 모델을 사용해 0.1 GeV부터 10 TeV까지 전 범위에서 계산하였다. 결과 스펙트럼은 dNγ/dEγ ∝ E⁻²·⁴ (1 GeV–100 GeV)와 dNγ/dEγ ∝ E⁻³·¹ (100 GeV–1 TeV) 형태로, Fermi‑LAT가 보고한 E⁻²·²와 HAWC가 보고한 E⁻³·⁰ 사이에 위치한다. 스펙트럼 기울기는 필라멘트형 네트워크가 주도하는 ‘거울 효과’와 알프벤 난류에 의한 ‘산란 억제’가 결합된 결과이며, 특히 100 GeV 이하에서 난류가 기울기를 크게 가파르게 만든다.

절대 플럭스가 관측보다 2–5배 낮은 원인으로는(1) 모델에 포함되지 않은 태양 활동 주기별 자기장 강화, (2) 높은 고도(> 5 Mm)에서의 추가적인 난류, (3) GCR 스펙트럼의 지역적 변동성을 들 수 있다. 저자들은 향후 고해상도 관측(예: DKIST)과 연계한 자기장 매핑, 그리고 전자·양성자 혼합 효과를 포함한 전자기적 전이 모델을 제안한다.

전반적으로, 이 연구는 ‘다중 스케일 자기장 + 알프벤 난류’라는 새로운 프레임워크를 제시함으로써, 태양 디스크 감마선이 저층 대류와 직접 연결된 물리적 메커니즘을 밝히고, GCR 전파 연구에 새로운 관측적 도구를 제공한다는 점에서 큰 의미를 가진다.