대기 중 우주선 양성자 이온화 수율 함수 Y의 afxy‑코드 기반 파라미터화

초록

본 연구는 시뮬레이션 데이터에 역비선형 문제 해결 기법인 afxy‑코드를 적용하여, 500 MeV부터 1 TeV까지의 1차 양성자에 대한 대기 이온화 수율 함수 Y를 간단한 유리함수 형태로 파라미터화한다. 제시된 식은 대기 전층(지표 ~ 상층)에서의 이온쌍 생성량을 직접 계산할 수 있게 하며, 태양 활동·지자기 변동에 따른 우주선 플럭스 변화와 기후·구름 형성 연구에 활용될 수 있다.

상세 분석

이 논문은 우주선 유도 대기 이온화의 핵심 변수인 이온화 수율 함수 Y(x,E)를 정량화하기 위해, 기존 Monte‑Carlo 시뮬레이션 결과를 역문제로 전환하고 afxy‑code( “analyze f(x)=y” )를 이용해 비선형 시스템 f(x)=y 를 해석한다. afxy‑code는 정규화된 Gauss‑Newton 반복식 (5)을 기반으로 하며, 프로세스와 문제 자체에 대한 정규화 파라미터 ε_k, α_k 를 자동 조정한다. 특히, ill‑conditioned 상황에서 L∞ 노름 기반 자동 정규화식 (6)과 특이값 기반 L2 정규화식 (7‑8)을 도입해 수렴성을 확보한다.

수치 해법은 로그 변환과 변수 스케일링을 병행함으로써, 10^17 수준의 조건 수를 가진 시스템에서도 100 회 이하의 반복으로 χ²≈0.001 수준의 적합을 달성한다. 파라미터화 모델은 2차 다항식의 분자·분모 형태인 Y(h)= (a h² + b h + c)/(d h² + e h + f) 로 설정했으며, h는 대기 깊이(g cm⁻²)이다. 각 에너지(500 MeV, 1 GeV, 5 GeV, 10 GeV, 50 GeV, 100 GeV, 500 GeV, 1 TeV)에 대해 a‑f 값을 최적화하였다.

결과적으로, 모든 에너지 구간에서 Pfotzer 최대(대기 약 15 km 고도)와 저층 대기에서의 Y값을 정확히 재현했으며, 특히 1 GeV와 5 GeV 구간에서는 시뮬레이션 데이터와 거의 일치한다. 50 GeV·100 GeV 구간에서는 Pfotzer 최대에서 약간의 과대평가가 있었지만 전체적인 오차는 미미하다. 파라미터는 간단한 지수·멱함수로 추가 보간이 가능해, 다양한 대기 모델과 지자기 위도 λ_m 에 대한 D(E,λ_m)와 결합해 식 (3) 형태의 이온쌍 생성률 q(h,λ_m) 를 손쉽게 계산할 수 있다.

이러한 분석은 기존에 복잡한 Monte‑Carlo 시뮬레이션을 매번 수행해야 했던 문제를 해소하고, 실시간 혹은 장기 기후 모델에 우주선 이온화 효과를 삽입하는 데 큰 장점을 제공한다. 또한, afxy‑code가 비선형 역문제 해결에 강인한 수렴성을 보이며, 다양한 물리‑수학 문제에 적용 가능함을 시연한다.