우주선 전파 파라미터의 체계적 불확실성 B/C 비율로 델타 1분의 3 혹은 1분의 2를 맞출 수 있을까

초록

B/C 비율을 이용해 우주선 전파 모델의 확산 지수 델타를 추정하면, 입력 가스 밀도·수소 비율·원천 스펙트럼·확산 계수의 저에너지 형태·핵분열 단면적 등 여러 변수에 따라 결과가 크게 달라진다. 저에너지에서의 불확실성이 통계적 오차보다 크며, 대류와 재가속을 동시에 포함한 모델은 델타 > 0.6를 요구해 이론적 기대값(1/3, 1/2)과 충돌한다. 반면 대류를 제외하고 저에너지 확산 계수를 급격히 상승시키면 델타 ≈ 1/3–1/2를 얻을 수 있다. 최종적으로는 고에너지(TeV/n) B/C 측정이나 다른 2차/1차 비율과의 결합이 필요하다.

상세 요약

이 논문은 현재까지 발표된 B/C 비율 기반 전파 파라미터 추정값이 0.3에서 0.8까지 넓은 범위에 걸쳐 산재하는 현상을 체계적인 입력 변수 변화로 설명하고자 한다. 저자들은 USINE 전파 코드와 Minuit 최소화 알고리즘을 이용해, (1) 은하 원반의 가스 밀도와 수소 비율, (2) 원천 스펙트럼의 파워‑로우 형태, (3) 확산 계수 D(R)=β^η D0 (R/R0)^δ 의 저에너지 지수 η, (4) 핵분열 단면적 데이터베이스의 선택·보정, 네 가지 요소를 독립적으로 변동시켰다. 각 경우마다 B/C 데이터에 대한 χ² 최소값을 찾고, 최적 파라미터(δ, D0, Vc, Va 등)의 통계적 불확실성과 시스템적 변동을 비교하였다.

주요 결과는 다음과 같다. 첫째, 가스 밀도와 수소 비율을 10 % 정도 바꾸어도 δ는 0.05–0.07 정도 변동한다. 둘째, 원천 스펙트럼의 저에너지 경사(γ0) 변화는 δ에 0.03–0.05 정도 영향을 미친다. 셋째, η를 -0.5에서 +1.0 사이로 조정하면 δ는 0.2 이상 크게 변한다. 특히 η>0인 경우, 저에너지에서 확산이 억제되어 B/C가 상승하고, 이를 보정하기 위해 모델은 큰 δ와 강한 대류 풍속 Vc≈15–20 km s⁻¹를 요구한다. 넷째, 핵분열 단면적을 최신 TALYS 기반 데이터와 전통적인 Webber‑Koch 데이터 사이에서 교체하면 δ가 0.04–0.06 정도 달라진다. 이러한 시스템적 변동은 통계적 오차(Δδ≈0.02)보다 훨씬 크다.

대류와 재가속을 동시에 포함한 모델(VC>0, Va>30 km s⁻¹)은 χ² 최소값이 가장 낮아 통계적으로 선호된다. 그러나 이 경우 δ는 0.65–0.80 범위에 머물며, Kolmogorov(δ=1/3)이나 Kraichnan(δ=1/2) 난류 이론과는 크게 차이 난다. 반면 대류를 없애고 η를 양수(≈0.8)로 크게 설정하면, δ≈0.35–0.45를 얻을 수 있다. 이때 확산 계수는 저에너지에서 급격히 상승하는 형태를 띠며, 이는 현재 관측된 B/C의 저에너지 플랫함을 설명한다. 하지만 이러한 “강한 저에너지 상승” 모델은 물리적 근거가 약하고, 재가속 파라미터 Va가 비현실적으로 낮아지는 경향이 있다.

결론적으로, 저에너지 영역(~ GeV/n)에서의 입력 불확실성이 δ 추정에 지배적인 영향을 미치며, 현재 B/C 데이터만으로는 이론적 기대값을 확정하기 어렵다. 저에너지 모델링을 개선하거나, TeV/n 수준의 고에너지 B/C 측정, 혹은 B/O, sub‑Fe/Fe 등 다른 2차/1차 비율과의 동시 분석이 필요하다. 이러한 추가 데이터가 확보되면, 대류·재가속·저에너지 확산 형태 사이의 파라미터 상쇄 효과를 깨고, δ를 이론적 값에 근접하게 결정할 수 있을 것으로 기대된다.