갤륨 이상 현상의 재평가: 자체 일관 전자 파동함수 활용

초록

본 논문은 기존의 전자 파동함수 근사법을 버리고, Dirac‑Coulomb 방정식을 수치적으로 풀어 바운드와 연속 전자 상태를 정확히 계산한다. 이를 통해 71Ga → 71Ge 중성자 캡처 단면을 재평가하고, 최신 실험(GALLEX, SAGE, BEST) 데이터를 이용해 갤륨 이상 현상의 통계적 유의성을 새롭게 추정한다. 결과는 스테릴 중성미자 해석에 필요한 파라미터 공간을 크게 축소한다.

상세 분석

이 연구는 갤륨‑71(71Ga) 핵을 이용한 전자 캡처와 중성미자 흡수 반응의 상세한 이론적 계산을 수행한다. 기존 문헌에서는 전자 파동함수를 전자기장에 대한 일차 근사(LO)로 처리해 전자 밀도와 페르미 함수(Fermi function)를 단순화했지만, 이러한 접근은 핵 내부에서 전자 파동함수의 급격한 변화를 충분히 반영하지 못한다는 한계가 있다. 저자들은 RADIAL 패키지를 기반으로 Dirac‑Hartree‑Fock‑Slater(DHFS) 방정식을 풀어, 바운드 전자(특히 1s 궤도)의 정확한 큰·작은 성분 g₁ₛ(r), f₁ₛ(r)와 연속 전자 파동함수 g_κ(r), f_κ(r)를 얻었다. 핵 전위 V_DHFS(r)는 핵 전하 분포(2pF 모델, R_ch = 4.032 fm)와 전자 구름, 교환 상호작용을 포함하며, Latter tail 보정을 통해 장거리 거동을 보정한다. 교환 상수 C_ex은 Slater(3/2)와 Kohn‑Sham(1) 사이의 평균값 1.25를 채택하고, 이에 대한 시스템틱 불확실성을 추가하였다.

핵심 결과는 두 가지이다. 첫째, 1s 전자 밀도 4π|ψ_be,1s(r₀)|²는 r₀ = 0 및 r₀ = (5/3)R_ch에서 기존 계산과 0.2 % 이내의 일치성을 보이며, RADIAL와 GRASP(전완전 Dirac 해석기) 사이의 차이는 1 % 미만이다. 이는 전자 밀도 계산에 있어 선택된 전위와 교환 모델이 충분히 정확함을 의미한다. 둘째, 일반화된 페르미 함수 F(E_e,Z,r) = (|f₁(r)|²+|g_{−1}(r)|²)/(2p_e²) 를 직접 계산함으로써, 전통적인 F₀·L₀·U·S 형태와 비교했을 때 동일한 값이 얻어짐을 확인했다. 특히 핵 표면이 아닌 원점에서 평가한 경우, 전자 스크리닝과 유한 핵 크기 보정이 자연스럽게 포함된다.

이러한 정밀한 전자 파동함수와 페르미 함수는 상세균형 원리(detailed balance)를 통해 전자 캡처 반감기 t₁/₂와 중성미자 흡수 단면 σ_gs를 연결한다. 저자들은 최신 전자 캡처 데이터(71Ge의 t₁/₂, 전자 캡처 비율 P_K, P_L, P_M 등)와 최신 Q값을 사용해 f_EC와 B_GT를 재계산하였다. 결과적으로 σ_gs는 기존 LO 기반 값보다 약 2–3 % 감소했으며, 이는 BEST, GALLEX, SAGE에서 보고된 관측/예측 비율(R)과의 차이를 약간 완화한다. 그러나 통계적 분석에 따르면, 전체 데이터셋을 합친 갤륨 이상 현상의 유의성은 여전히 3.8σ 수준으로 남아 있다.

스테릴 중성미자 해석에서는 두 자유 파라미터(질량 차이 Δm²₁₄, 혼합각 sin²2θ₁₄)를 고려한다. 재계산된 σ_gs를 사용해 전역 χ² 분석을 수행한 결과, 기존에 제시된 허용 영역보다 약 30 % 축소된 영역이 남았다. 특히 Δm²₁₄ ≈ 1.3 eV², sin²2θ₁₄ ≈ 0.1 근처가 최적점으로 나타났으며, 이 영역은 다른 단거리 실험(NEOS, DANSS 등)과의 텐션을 완화한다.

요약하면, 전자 파동함수의 전자기적 왜곡을 정확히 다루는 것이 갤륨 이상 현상의 이론적 불확실성을 크게 줄였으며, 스테릴 중성미자 가설에 대한 실험적 제약을 더욱 강력히 만든다. 향후 연구는 핵 구조 효과(예: 핵 변형, 전이 강도)와 고차 전자 상호작용(진공 편극, 자기 상호작용)까지 포함해 오차 범위를 더욱 축소하는 것이 필요하다.