아리랑 입자 흐름의 방위 변조와 지자기 효과

초록

ARGO‑YBJ 실험에서 대기 중 2차 입자들의 지자기 라렌츠 힘에 의해 발생하는 방위(azimuth) 분포의 비균일성을 조사하였다. 데이터와 시뮬레이션을 통해 1차와 2차 조화 성분을 분리하고, 1차는 순수 지자기 효과, 2차는 지자기와 검출기 구조 비대칭이 결합된 결과임을 확인했다.

상세 분석

본 연구는 고도 4300 m에 위치한 ARGO‑YBJ 전천공 입자 검출기가 지자기장에 의해 2차 입자들의 측면 분포가 비대칭적으로 늘어나는 현상을 어떻게 감지하고 정량화하는지를 상세히 다룬다. 먼저, 라렌츠 힘에 의해 전자·양전자가 평균 고도 h 와 에너지 Eₑ 에 따라 횡방향으로 이동하는 평균 변위 d 를 식(1)로 정의하고, 이 변위가 샤워 풋프린트의 밀도를 감소시켜 트리거 효율을 각 방위 φ 에 따라 변조시킨다. 이론적으로는 방위 분포 dN/dφ 가 상수 N₀ 에 1차와 2차 코사인 항을 더한 형태(식 2)로 나타나며, 1차 진폭 g₁ 은 sin²θ 에 비례하고 위상 φ₁ 은 지자기 방위 −φ_B 와 일치한다. 2차 진폭 g₂ 도 동일한 sin²θ 의 의존성을 갖지만, 검출기 비대칭(패드 밀도 차이)으로 인해 추가적인 180° 주기의 변조가 겹칠 수 있다.

데이터는 2010년 10월 12‑14일 3일간 수집된 1억 3000만 건의 이벤트 중, 코어가 중앙 40 × 40 m² 내에 있고 θ < 60° 인 1300만 건을 선택하였다. 타이밍 보정과 프리모듈레이션 보정을 거친 뒤 방위 히스토그램을 얻었으며, 식 2에 대한 최소제곱 피팅 결과 g₁ = 1.521 % ± 0.012 %, φ₁ = −72.4° ± 0.47°, g₂ = 0.587 % ± 0.012 %, φ₂ = −86.5° ± 0.60°를 얻었다. g₁ 은 기대값보다 약 3배 큰 반면 φ₁ 은 −φ_B와 거의 일치한다. 반면 φ₂ 는 −φ_B와 −90° 사이에 위치해 검출기 비대칭을 시사한다.

다음으로 θ 에 따른 g₁, g₂ 의 변화를 조사했으며, g₁ 은 sin²θ 에 정확히 비례함을 확인했다(그림 3). g₂ 는 기대와 달리 일정하지 않아, 식 6으로 두 개의 2차 항을 분리하였다. 첫 번째는 지자기 효과(2B), 두 번째는 검출기 구조 효과(2A)로 가정하고, 세 개의 θ 구간(α, β, γ)에서 동시에 피팅하였다. 결과는 k₁ = 2.094 % ± 0.016, k₂ᴮ = 0.87 % ± 0.32, k₁/k₂ᴮ = 2.41 ± 0.90(예상 ≈ 1.9)이며, φ₁ = −71.86°와 φ₂ᴬ = −92.5° ± 2.9°(≈ −90°)가 얻어졌다. 이는 검출기 x‑축과 y‑축의 패드 밀도 차이(1.54 vs 1.76 pads/m)가 트리거 효율에 방위 의존성을 부여한다는 물리적 해석과 일치한다.

시뮬레이션에서는 CORSIKA와 GEANT4 기반 코드를 이용해 θ = 45°, E = 3 TeV 인 양성자 빔을 9가지 방위와 3가지 자기장 세기(0, B, 2B)로 전파하였다. 결과는 sinχ 가 큰 방위에서 트리거 효율이 최대 ~ 30% 감소함을 보여준다(표 2). 또한, 전체적인 프로톤·헬륨 샤워 시뮬레이션(10 GeV–10 PeV)에서도 방위 비균일성이 재현되었으나 통계가 부족해 정확한 진폭 추정이 어려웠다.

결론적으로, ARGO‑YBJ에서 관측된 방위 변조는 1차 코사인 항이 순수 지자기 라렌츠 힘에 의해, 2차 항은 지자기와 검출기 비대칭이 복합적으로 작용한 결과임을 명확히 규명하였다. 이는 향후 EAS 배열의 절대 포인팅 정확도 검증 및 달 그림자 분석에 중요한 교정 요소가 될 수 있다.