Modeling the Effect of the Heliospheric Magnetic Field on Cosmic Ray Muon Shadows

기술 요약: 헬리오스pheric 자기장이 우주선 뮤온 그림자에 미치는 영향 모델링

문제 제기
우주선 (CR) 은 주로 양성자와 알파 입자로 구성되어 있으며, 태양계를 통과할 때 전자기장에 의해 편향됩니다. 지구 자기장 (GMF) 은 잘 매핑되어 있지만, 헬리오스pheric 자기장 (HMF) 은 복잡하고 모델링하기 어렵습니다. 태양과 달은 입사하는 우주선을 흡수하여 지구에서 관측 가능한 우주선 플럭스 내 "그림자"를 생성합니다. 그러나 이러한 그림자는 정적이지 않으며, HMF 와 GMF 에 의해 이동하고 왜곡됩니다. MINOS+ 협력단은 태양 주기 (태양 극소기, 태양 극대기, 그리고 전체 13 년 기간) 의 다양한 위상 동안 이러한 태양 그림자를 위치시키기 위해 MINOS 원거리 검출기의 뮤온 데이터를 활용했습니다. 다루어진 핵심 문제는 관측된 태양 뮤온 그림자의 이동과 왜곡이 HMF 의 표준 파커 나선 모델과 GMF 의 쌍극자 모델과 일관성이 있는지, 그리고 그렇지 않다면 이론과 관측을 조화시키기 위해 모델이나 매개변수에 어떤 수정이 필요한지입니다.

방법론
본 연구는 소우단 광산 지하 2070 미터 수환당 (water equivalent) 깊이에 위치한 MINOS 원거리 검출기에서 수집된 2003 년 7 월부터 2016 년 7 월까지의 13 년간 뮤온 데이터를 활용했습니다. 분석은 입사 천정각이 80°미만이고 궤적 재구성을 위한 충분한 트랙 길이를 가진 단일 뮤온 사건에 초점을 맞추었습니다.

1. 그림자 위치 확인: 연구팀은 그림자를 위치시키기 위해 $\Lambda$ 확률 테스트 (로그 가능도 비) 를 사용했습니다. 이는 태양 근처에서 관측된 입자 수 ($N_{obs}$) 를 이론적 평탄 배경 ($N_{back}$) 과 이론적 그림자 템플릿 ($N_{th}$) 과 비교하는 과정을 포함했습니다. 템플릿은 태양의 각지름 (0.52°) 과 암석 상부 덮개에서의 산란 효과를 고려하였습니다.
2. 통계적 검증: 몬테카를로 시뮬레이션을 사용하여 무작위 배경 히스토그램을 생성하여 관측된 그림자 위치의 유의성을 결정했습니다. 무작위 요동이 그림자를 모방할 확률을 계산하여 높은 신뢰 수준 (전체 데이터셋의 경우 최대 6$\sigma$) 을 도출했습니다.
3. 계산 모델링: 그림자 이동을 해석하기 위해 저자는 검출기에서 태양까지 우주선의 시간 역행 궤적을 시뮬레이션했습니다.
   • 장 모델: 시뮬레이션은 HMF 에 대한 파커 나선 모델과 GMF 에 대한 쌍극자 모델을 중첩시켰습니다. 파커 나선 매개변수에는 일정한 태양풍 속도 ($u_{sowi} = 400$ km/s) 와 태양 회전 주파수가 포함되었습니다. GMF 는 쌍극자 모멘트를 사용하여 모델링되고 태양 중심 좌표계로 변환되었습니다.
   • 입자 역학: 로런츠 힘 방정식은 4 차 룽게 - 쿠타 방식을 사용하여 적분되었습니다. 입자들은 상대론적 입자로 처리되었으며, 에너지는 멱함수 분포 ($E^{-2.75}$) 를 따랐습니다.
   • 시나리오: 시뮬레이션은 다양한 자기장 극성 비율 (50/50, 33/67, 67/33) 과 두 가지 그림자 템플릿 유형 (단순 원형 그림자와 실제 $\Lambda$ 확률 히스토그램에서 유도된 "왜곡된" 그림자) 을 테스트했습니다.
   • 매개변수 변화: 1 차 우주선 에너지 차단값 ($E_0$) 은 모델이 관측된 그림자 위치와 일치하는지 확인하기 위해 자유 매개변수로 처리되었습니다.

주요 기여

• 그림자 특성화: 본 논문은 전체 13 년 데이터셋, 태양 극소기 (2007–2009), 그리고 태양 극대기 (2012–2014) 의 세 가지 명확한 시간 구간에서 높은 통계적 유의성을 가진 태양 뮤온 그림자를 성공적으로 위치시켰습니다.
• 템플릿 개발: 저자들은 완벽한 원형 그림자에 대한 가정을 넘어, 관측된 $\Lambda$ 확률 분포에 기반한 "왜곡된" 그림자 샘플링 템플릿을 생성하는 방법을 도입했습니다.
• 모델 테스트: 본 연구는 파커 나선 모델을 관측 데이터에 대해 엄격하게 테스트하여 시뮬레이션된 그림자 이동과 관측된 그림자 이동 사이의 일치 정도를 적경 ($\Delta RA$) 과 적위 ($\Delta DEC$) 에서 정량화했습니다.
• 에너지 스펙트럼 분석: 본 논문은 그림자 위치가 1 차 우주선 에너지 스펙트럼에 얼마나 민감한지, 특히 저에너지 차단값 ($E_0$) 의 영향을 테스트하여 조사했습니다.

결과

• 그림자 이동: 관측된 그림자는 태양 주기 위상에 따라 뚜렷한 이동을 보였습니다. 태양 극소기 그림자는 태양 극대기 그림자와 다르게 이동했으며, 전체 데이터셋은 이러한 이동들의 조합을 나타냈습니다.
• 모델 일관성:
  • 파커 나선 모델은 태양 극소기 데이터와 가장 잘 일치했습니다.
  • 이 모델은 태양 극대기 데이터, 특히 $\Delta DEC$ 이동과 관련하여 가장 일치하지 않았습니다.
  • 왜곡된 그림자 템플릿은 태양 극대기와 전체 데이터 기간에서 더 잘 수행된 반면, 단순한 원형 템플릿은 태양 극소기 데이터에서 더 효과적이었습니다.
• 극성 효과: 시뮬레이션은 자기장 극성 비율이 그림자의 $\Delta DEC$ 위치에 직접적인 영향을 미친다는 것을 확인했습니다. 데이터는 특히 태양 극대기 동안 양의 극성이 지배적인 HMF 를 선호하는 것으로 나타났습니다.
• 에너지 차단 불일치: 파커 나선 모델을 사용하여 시뮬레이션된 그림자 평균을 이론적 중심 (0,0) 에 정렬시키기 위해, 1 차 우주선 에너지 차단값 ($E_0$) 을 22.8 $\pm$ 1.2 TeV로 조정해야 했습니다. 이는 실험적으로 결정된 약 10 TeV 값보다 현저히 높습니다 (2~3 배).
• 천정각 보정: 에너지 차단값의 천정각 의존성에 대한 보정이 적용되었습니다. 이는 관측된 뮤온 에너지 분포와 스펙트럼을 더 잘 일치시키도록 조정되었지만, 관측된 $\Delta RA$ 이동을 완전히 설명하기에는 부족했습니다.

의의 및 결론
본 논문은 우주선 태양 뮤온 그림자가 관측 가능하며 HMF 구조를 탐구하는 데 사용될 수 있다고 결론지었습니다. 결과는 태양 극소기 그림자가 표준 파커 나선 모델과 쌍극자 GMF 의 조합과 가장 일관됨을 나타냅니다. 그러나 모델은 비현실적으로 높은 에너지 차단값을 도입하지 않고는 태양 극대기 그림자와 전체 데이터셋을 완전히 재현하지 못합니다.

저자들은 이 불일치가 파커 나선 모델이 태양 극대기 동안 존재하는 조건에 대한 과도한 단순화임을 시사한다고 주장합니다. 그들은 중성 전류층 기울기 ("발레리나 모델"), 충격 전면, 그리고 복잡한 극성 변화와 같은 모델링되지 않은 요인들이 관측된 이동에 기여할 가능성이 있다고 제안합니다. 본 연구는 편향의 일반적인 방향이 파커 나선과 일치하지만, 태양 극대기 동안의 HMF 의 상세한 구조는 단순 나선 모델이 제공하는 것보다 더 복잡한 모델링이 필요함을 확립합니다. 이 작업은 HMF 모델을 제약하는 데 뮤온 그림자 관측의 유용성을 보여주지만, 태양 주기 변동을 설명하기 위해 더 정교한 시뮬레이션이 필요함을 강조합니다.