핵분열에서 바륨과 크세논 동위 원소 수율의 상세 분석

초록

이 논문은 4차원 랭게뱅 모델과 Fourier‑over‑Spheroid 형상 파라미터를 이용해 다양한 액티늄 핵의 자발 및 중성자 유도 핵분열에서 바륨·크세논 동위 원소 수율을 계산하고, 평가된 실험 데이터와 비교한다. 평균 전하 분배와 주요 분열 채널의 평균 중성자 수는 잘 재현되지만, 계산된 분포는 꼬리 부분에서 급격히 감소해 실험값보다 좁게 나타나는 체계적인 편차가 확인된다.

상세 요약

본 연구는 4차원(Langevin) 동역학 모델에 Fourier‑over‑Spheroid 형상 파라미터화를 도입하여, 분열 과정에서 발생하는 복합적인 잠재 에너지와 마찰·확산 항을 정밀하게 기술한다. 모델은 4개의 자유도(두 축 비율, 비대칭성, 그리고 길이 변형)를 동시에 추적함으로써, 전통적인 2차원 또는 3차원 접근법이 놓치는 미세한 형상 변화를 포착한다. 특히 전하 분할 Z와 중성자 수 N의 상관관계를 실시간으로 계산해, 특정 Z에 고정된 동위 원소 사슬 Y(N_f) 를 직접 얻을 수 있다.

벤치마크 대상으로는 ^244,246Cm 및 ^250Cf의 자발 분열, 그리고 ^229Th, ^235U, ^239Pu, ^249Cf의 열중성자 및 14 MeV 중성자 유도 분열을 선택하였다. 평가된 데이터베이스(ENDF/B‑VIII, JEFF‑3.3 등)와 비교했을 때, 주요 분열 채널(예: 1‑2 MeV 이하의 저에너지 채널)에서 평균 전하 분할과 중성자 수 최대값(N_max)은 0.5 ~ 1 % 이내의 오차로 재현된다. 이는 모델이 평균적인 잠재 에너지 최소점과 전하-중성자 비율을 정확히 잡아낸다는 강점을 보여준다.

그러나 전반적인 동위 원소 분포 폭(σ)은 실험값보다 현저히 좁게 계산된다. 특히 무거운 파편(바륨 사슬)에서 꼬리 부분이 급격히 감소하는 경향이 두드러진다. 이는 두 가지 원인으로 해석된다. 첫째, 마찰 계수와 확산 텐서의 파라미터화가 실제 핵 내부의 복잡한 열역학적 변동성을 충분히 반영하지 못했을 가능성이다. 둘째, 형상 파라미터 공간에서의 경계 조건(예: 스핑클링 형태에서의 전하 분할 제한)이 과도하게 제약을 가해, 비대칭적인 극단 형상으로의 전이 확률을 억제했을 것으로 보인다.

또한, 무거운 파편과 가벼운 파편 사이의 상관관계(예: Z_heavy + Z_light = Z_CN)와 중성자 수 보존(N_heavy + N_light + ν = N_CN)도 모델이 전반적으로 만족하지만, 특정 채널에서 ν(방출 중성자 수)의 분산이 과소평가되어 전체 분포 폭이 축소되는 현상이 관찰된다. 이러한 결과는 향후 마찰·확산 파라미터의 에너지 의존성 도입과, 형상 파라미터의 고차원 확장을 통해 개선될 여지가 있음을 시사한다.

요약하면, 4D Langevin‑Fourier‑over‑Spheroid 모델은 평균적인 분열 특성(전하 분할, 중성자 수 최대값) 재현에 성공했으며, 특히 중성자 수 최대값이 핵종마다 일관되게 나타나는 점은 모델의 물리적 타당성을 뒷받침한다. 그러나 분포 폭과 꼬리 부분의 체계적인 과소평가 문제는 마찰·확산 메커니즘과 형상 파라미터의 제한 조건을 재검토해야 함을 강조한다. 향후 이러한 개선을 통해 핵연료 사이클, 방사성 폐기물 관리, 핵무기 비확산 등 실용적인 분야에서 요구되는 고정밀 동위 원소 수율 예측이 가능해질 것이다.