REAS3 지오신크로트론 모델 개정 구현과 라디오 방출 예측

초록

**
REAS3은 기존 REAS2 코드의 결함을 보완하여 입자 수 변동에 의한 전자기 복사를 ‘엔드‑포인트’ 방식으로 구현하였다. 이를 통해 순수 기하학적 방사뿐 아니라 전하 과잉에 의한 방사도 시뮬레이션에 포함시켰으며, 펄스 형태가 단극에서 양극으로 바뀌고 주파수 스펙트럼이 0 Hz에서 0으로 수렴한다. LOPES 실험 데이터와의 비교에서도 측면 분포가 크게 개선되어 REAS3이 실제 공기 샤워 라디오 방출을 보다 정확히 재현함을 확인하였다.

**

상세 분석

**
본 논문은 REAS2에서 누락된 ‘입자 수 변동에 의한 전자기 복사’를 보완하기 위해 엔드‑포인트(formulation) 접근법을 도입한 점이 가장 큰 혁신이다. 기존 REAS2는 입자 궤적의 연속 구간에서만 가속도에 의한 복사를 계산했으며, 입자가 생성·소멸되는 순간(시작·종료점)에서 발생하는 급격한 속도 변화는 무시되었다. 이는 실제 공기 샤워에서 전하 과잉(electron excess)이 시간에 따라 변하면서 발생하는 라디오 신호를 과소평가하게 만든다.

엔드‑포인트 방식은 입자 궤적을 ‘직선 구간 + 급격한 꺾임(kink)’으로 분해하고, 각 꺾임에서 속도 변화 Δβ를 이용해 복사 전기장을 시간 적분 형태로 계산한다(식 1). 이때 시작점·종료점은 하나의 속도가 0인 특수한 꺾임으로 취급되어, 전하 생성·소멸에 따른 복사가 자연스럽게 포함된다. 이렇게 하면 전자·양성자 등 모든 입자에 대해 동일한 수식이 적용되므로 모델이 완전히 보편적이며, 별도의 ‘전하 과잉 전용’ 모듈이 필요 없게 된다.

시뮬레이션 결과는 두드러진 변화를 보여준다. 첫째, 펄스 형태가 단극에서 양극(bipolar)으로 바뀌어, 복사가 유한한 시간·공간 영역에서만 발생한다는 물리적 기대와 일치한다. 둘째, 주파수 스펙트럼이 저주파(≈0 Hz)에서 0으로 수렴하여, MGMR 모델과의 차이가 크게 감소한다. 셋째, 전하 과잉에 의한 방사는 지자기장과 무관하게 발생하므로, 관측자 위치에 따라 동쪽‑서쪽 비대칭이 나타난다. 이는 REAS3이 순수한 v × B(geomagnetic) 방사만을 고려하던 REAS2와 근본적으로 다른 점이다.

또한, LOPES 실험 데이터와의 비교에서 REAS3은 측면(라테럴) 분포의 기울기가 REAS2보다 평탄해 실제 관측값과 거의 일치한다. 특히, 기존 모델이 과소평가하던 ‘플랫 라테럴 분포’를 성공적으로 재현했으며, 이는 전하 과잉 복사가 외부 전자기장에 의존하지 않는 추가적인 방사원을 제공하기 때문이다. 마지막으로, REAS3은 자유 파라미터가 전혀 없으며, 입력으로 CORSIKA에서 제공되는 입자 분포와 샤워 매개변수만을 사용한다는 점에서 물리적 일관성이 높다.

**