Monte Carlo으로 핵분열 관측값 변동성 완전 분석

초록

본 논문은 입자수 보존을 투영한 Bogoliubov 진공 상태에서 핵분열 조각들의 다양한 물리량 확률분포함수(pdf)를 Monte Carlo 샘플링으로 추정하는 방법을 제시한다. 위치·스핀 좌표에서 핵자 구성을 무작위로 생성하고, 일·이체 상호작용을 이용해 총운동량·전하·형상 변동 등을 계산한다. 실험적으로 중요한 총운동에너지·각운동량 등도 평균적인 평균장(mean‑field) 수준에서 이미 큰 변동을 보임을 확인한다.

상세 분석

이 연구는 에너지 밀도 함수(EDF) 기반의 평균장 이론에서 발생하는 양자 요동을 정량화하기 위해, 입자수 보존을 구현한 Bogoliubov 진공을 직접 샘플링하는 새로운 Monte Carlo 프레임워크를 구축하였다. 핵자 각각을 (r, σ) 좌표와 스핀 투영으로 기술하고, 정규화된 확률밀도 p(r₁…r_N)=|⟨r₁…r_N|ψ⟩|²/N! 를 기반으로 Metropolis‑Hastings 알고리즘을 적용한다. 여기서 핵자 하나씩을 작은 공간적 변동(σ_space≈1 fm)과 스핀 플립(P_flip≈10 %)으로 제안함으로써 높은 수용률을 유지하고, 충분한 burn‑in과 jump 간격을 통해 독립 표본을 확보한다.

관측가능 연산자는 위치·스핀 표현에서 대각형(diagonal)이어야 하며, 이는 핵자 밀도·전하·운동량 등 1‑body 및 2‑body 연산자를 포함한다. 대각형 연산자는 샘플링된 구성에서 직접 계산된 커널 O(r₁…r_N) 로 표현되며, 그 확률분포는 바로 물리량 측정값의 pdf와 일치한다. 따라서 평균값뿐 아니라 전 구간의 분포, 높은 모멘트까지 자연스럽게 얻을 수 있다.

Bogoliubov 진공의 파동함수 절대값 제곱은 Z 행렬식(det Z) 형태로 표현되며, Z_ij=Z(r_i,r_j)−Z(r_j,r_i) 로 정의된다. 여기서 Z(r,r′)=∑_k v_k u_k φ_k(r) φ̄_k(r′) 로 전개되며, φ_k는 정준 기저의 단일입자 파동함수이다. 계산 효율을 위해 정준 상태를 10⁻⁴ 입자 오차 기준으로 절단하고, 격자화된 박스(L_x≈32 fm, dx≈0.5 fm) 위에서 φ_k를 상수화하여 근사한다. 이러한 근사는 다중극자 모멘트와 전하 분포에 미치는 영향을 최소화한다.

수치 실험에서는 가벼운 ²⁰Ne와 중핵 ²⁵²Cf의 스키션 전후 상태를 대상으로 검증하였다. ²⁰Ne에서는 600여 개 자유도, ²⁵²Cf에서는 6000여 개 자유도를 가진 고차원 샘플링이 수행되었으며, 수천에서 수만 개의 독립 표본을 확보해 pdf 수렴성을 확인했다. 결과는 조각 형상·전하·각운동량·상호작용 토크 등 다양한 물리량이 평균장 수준에서도 상당한 변동성을 보이며, 실험적으로 관측되는 폭넓은 분포의 상당 부분을 설명한다는 점을 시사한다.

이 방법은 기존의 투사 기법이 다루기 어려웠던 총운동에너지(TKE)와 같은 2‑body 관측값의 전 확률분포를 제공함으로써, 향후 fission fragment의 동역학 모델링과 실험 데이터 해석에 중요한 도구가 될 전망이다. 또한, NucleoScope라는 오픈소스 구현을 통해 커뮤니티가 손쉽게 접근하고 확장할 수 있도록 하였다.