Angular and Kinetic Properties of Scission Neutrons within Time-dependent Density Functional Theory

본 연구는 시간 의존 밀도 범함수 이론을 활용하여 절단 중성자(scission neutrons)가 $^{235}\mathrm{U}$, $^{239}\mathrm{Pu}$, 그리고 $^{252}\mathrm{Cf}$ 핵분열에서 즉각 핵분열 중성자 스펙트럼의 유의미한 고에너지 성분을 구성함을 입증하며, 이들의 포함이 기존의 증발 전용 모델에서 관찰된 고에너지 수율의 체계적 과소평가를 해결한다는 직접적인 증거를 제공한다.

핵심

거대하고 불안정한 풍선(원자핵)이 갑자기 두 개로 툭 끊어지는 모습을 상상해 보세요. 이것이 바로 핵분열입니다. 수십 년 동안 과학자들은 이런 현상이 일 발생할 때, 결과물인 두 조각(파편)이 엄청난 속도로 튕겨 나가며 중성자라고 불리는 아주 작은 입자들을 내뿜는다는 사실을 알고 있었습니다.

오랫동안 물리학자들은 이 중성자들이 마치 뜨거운 커피 잔에서 피어오르는 김처럼, 나중에 "증발"되는 것이라고 생각했습니다. 그들은 파편들이 완전히 형성되어 안정적으로 움직이는 과정에서 중성자를 놓아주는 것이라고 가정했습니다.

하지만 이 새로운 논문은 일부 중성자들이 바로 풍선이 끊어지는 그 순간에 "차여져서" 나온다는 것을 시사합니다. 이들을 **절단 중성자(scission neutrons)**라고 부릅니다. 이들은 나중에 식어가는 파편으로부터 나오는 것이 아니라, 분열의 찰나적인 혼돈 속에서 태어납니다.

연구진이 이 "끊어지는 순간"의 중성자들을 찾아낸 증거를 아주 쉽게 설명하면 다음과 같습니다.

1. 슈퍼컴퓨터 시뮬레이션

분열 중에 어떤 일이 일어나는지 보기 위해, 과학자들은 현미경을 사용하는 대신 **시간 의존 밀도 범함수 이론(TDDFT)**이라는 이론을 사용하여 이 사건의 영화를 실행하는 슈퍼컴퓨터를 사용했습니다.

이것은 원자들이 춤추며 부서지는 모습을 시뮬레이션하는 고속 3D 비디오 게임과 같습니다. 이전 버전의 이 "게임"에서는 가상 세계가 너무 작았습니다. 중성자들이 과학자들이 정확히 얼마나 빠른 속도로 어느 방향으로 날아가는지 파악하기도 전에 화면 가장자리에 부딪혀 버렸기 때문입니다.

이번 연구에서 그들은 훨씬 더 큰 가상 세계(이전보다 약 3배 더 큰)를 구축했습니다. 이를 통해 중성자들이 충분히 날아가서 안착할 수 있는 공간을 확보했고, "벽"이 데이터에 영향을 주지 않는 상태에서 중성자의 속도와 방향을 정확하게 측정할 수 있었습니다.

2. "속도 제한"의 발견

명확한 시야를 확보한 후, 연구진은 특정 각도(주로 분열 방향의 옆쪽이나 약간 뒤쪽)로 날아가는 중성자들을 관찰했습니다. 그들은 놀라운 사실을 발견했습니다.

• "출입 금지" 구역: 에너지 수준이 낮은(약 150만~200만 전자볼트 미만) 절단 중성자는 존재하지 않습니다. 마치 속도 제한이 있는 것처럼, "끊어지는 순간"의 중성자 중 느린 것은 허용되지 않습니다.
• 고속 집단: 대신, 이 중성자들은 모두 빠릅니다. 이들은 특정 높은 속도(3~3.5 MeV) 근처에 모여 있으며, 그 이후에는 훨씬 더 빠른 입자들로 이어지는 긴 꼬리를 형성합니다.

이것은 다이빙대에서 뛰어내리는 사람들의 무리와 같습니다. "증발된" 중성자들이 나중에 여유롭게 수영장 데크를 걷어차는 사람이라면, "절단" 중사들은 다이빙대가 부서지는 바로 그 순간에 격렬하게 튕겨 나가는 사람들입니다. 다이빙대에서 튕겨 나간 사람들은 항상 빠르게 움직이며, 그 사건으로부터 느리게 움직이는 사람은 결코 볼 수 없습니다.

3. "사라진" 에너지의 미스터리 해결

과학자들은 수년 동안 컴퓨터 모델을 실제 실험 결과와 일치시키기 위해 노력해 왔습니다. 그들에게는 한 가지 문제가 있었습니다.

• 기존 모델: 만약 "증기"(증발된 중성자)만을 계산한다면, 컴퓨터 모델은 고에너지 중성자가 너무 적다고 예측합니다. 이는 작은 컵으로 양동이를 채우려는데, 양동이가 필요로 하는 물의 양이 컵이 제공할 수 있는 양보다 계속 많은 것과 같습니다.
• 새로운 모델: 연구진이 "끊어지는 순간"의 중성자(그들이 큰 시뮬레이션에서 발견한 것)를 "증기" 중성자에 더했을 때, 수학적 계산이 마침내 맞아떨어졌습니다. 결합된 모델은 우라늄과 칼리포늄에 대해 측정된 실제 실험의 고에너지 데이터를 완벽하게 일치시켰습니다.

4. 이것이 중요한 이유

이것은 순수하게 미시적인 이론(무언가가 존재한다고 단순히 추측하거나 가정하지 않는 이론)이 이 "끊어지는 순간"의 중성자들을 예측하고 그것이 실재함을 증명한 첫 사례이기에 매우 중요합니다.

• 이전에는: 과학자들은 수학적 계산이 맞지 않기 때문에 이 중성자들이 존재할 것이라고 추측해야만 했습니다.
• 이제는: 컴퓨터 시뮬레이션이 별도의 지시 없이도 이 중성자들을 자연스럽게 생성해 냈습니다. 이는 단순히 기상 보고를 보고 폭풍이 올 것이라고 가정하는 것이 아니라, 구름의 움직임을 관찰하여 폭풍을 예측하는 것과 같습니다.

핵심 요약

논문은 원자가 분열될 때, 중성자의 작은 부분(전체의 약 6%~10%)이 분열의 격렬한 순간에 태어난다고 결론짓습니다. 이 중성자들은 특정 각도 내에서 항상 빠르고 결코 느리지 않다는 점에서 구별됩니다.

데이터에서 이 "지문"을 찾아냄으로써, 연구진은 마침내 "끊어지는 순간"의 중성자를 "증기" 중성자와 분리해 냈으며, 핵분열이 실제로 어떻게 일어나는지에 대한 더 명확하고 정확한 그림을 제시했습니다. 이는 물질을 결합하고 해체하는 근본적인 힘에 대한 우리의 이해를 정교화하는 데 도움을 줍니다.

기술 요약

기술 요약: 시변 밀도 범함수 이론(TDDFT) 내 절단 중성자(Scission Neutrons)의 각도 및 운동적 특성

문제 제기
핵분열의 발견으로부터 거의 한 세기가 지났음에도 불구하고, 특히 외부 안장점(outer saddle point)에서 절단(scission)으로의 전이와 그에 따른 1차 핵분열 파편(FFs) 형성 과정에 대한 완전한 미시적 묘사는 여전히 난제로 남아 있다. 현재의 통계적 증발 코드는 모든 즉발 중성자가 완전히 가속된 파편으로부터 등방적으로 방출된다고 가정한다. 그러나 이러한 프레임워크는 목(neck)의 파열로 알려진 고도의 비평형 역학적 과정 중에 방출되는 중성자인 절단 중성자(SNs)를 설명할 수 없다. SN은 1939년부터 가설로 제기되어 다양한 모델을 통해 연구되어 왔으나, 즉발 중성자 관측량만으로는 절단 시점에 방출된 중성자와 나중에 증발된 중성자를 구별할 수 없기 때문에 실험적으로 해결되지 않은 상태로 남아 있다. 기존의 시간 의존 밀도 범함수 이론(TDDFT)을 이용한 예측 시도는 시뮬레이션 영역의 제한으로 인해, 중성자가 계산 경계에 도달하기 전에 각도 및 운동적 특성을 동시에 결정하는 데 한계가 있었다.

방법론
본 연구는 $^{235}\text{U}(n_{th}, f)$, $^{239}\text{Pu}(n_{th}, f)$, 그리고 $^{252}\text{Cf}(sf)$에서의 SN 방출을 조사하기 위해 초유체 시스템으로 확장된 TDDFT를 기반으로 한 완전한 미시적 프레임워크를 채택하였다.

• 시뮬레이션 설정: 저자들은 SkM* 에너지 밀도 범함수(EDF)를 사용하는 축대칭 조화 진동자 솔버인 AxialHOHFB를 활용하였다. 핵심적인 발전 사항은 이전 연구들보다 약 3.2배 더 큰 시뮬레이션 영역(반축 $Z_{max} = 70$ fm, $R_{max} = 55$ fm인 타원체)을 사용했다는 점이다.
• 2단계 적분: 계산적 제약을 관리하면서 정확도를 유지하기 위해 2단계 시간 적분 절차를 사용하였다. 시스템을 작은 기저($B_1$)에서 정적 구속 헤세-보골리우프(HFB) 해로부터 절단 순간까지 진화시킨 후, 중성자가 경계에 도관할 때까지 진화를 계속하기 위해 더 큰 기저($B_2$)로 변환하였다.
• SN 특성 추출: 집단 흐름 에너지 밀도와 내적 운동 에너지 밀도의 비율이 임계값 $\eta = 3$을 초과하는 영역 $V^{(SN)}_\eta$를 정의함으로써 SN을 식별하였다. 이 기준은 외부 집단 흐름에 의해 지배되는 중성자를 분리해 낸다.
• 각도 및 운동 분석: 이 영역 내에서 국소 속도장과 중성자당 집단 흐름 에너지를 계산하였다. 연구는 경계 반사가 일어나기 전 SN 수율이 포화되는 각도 구간 $110^\circ \leq \theta \leq 150^\circ$ (경량 파편 방향 기준)에 집중하였다.
• 실험과의 비교: 추출된 SN 운동 에너지 스펙트럼을 증발 중성자(EN) 성분에 대한 맥스웰 분포 증발 모델과 결합하였다. 증발 모델의 파라미터(파편 온도 $T$ 및 중성자 다중도 비율 $\bar{\nu}_L/\bar{\nu}_H$)는 오직 저에너지 실험 데이터($E \leq E_{thr}$)에 의해서만 제약되었다. 이 결합 모델은 $^{239}\text{Pu}$ 및 $^{252}\text{Cf}$의 고에너지 즉발 핵분열 중성자 스펙트럼(PFNS) 데이터와 비교되었다.

주요 결과

• SN 스펙트럼 특성: 세 가지 시스템 모두에서 추출된 SN 스펙트럼은 뚜렷한 임계값 거동을 보인다. SN은 특정 에너지 임계값($E_{thr}$) 미만에서는 존재하지 않으며, 이 값은 약 1.5–2 MeV(유도 분열의 경우 $\approx 2$ MeV, 자발 분열의 경우 $\approx 1.5$ MeV)이다. 스펙트럼은 3–3.5 MeV 근처에서 정점을 찍으며, 약 18 MeV까지 뻗어 있는 지수적 고에너지 꼬리(tail)를 특징으로 한다.
• 수율 추정: $^{239}\text{Pu}(n_{th}, f)$의 경우, 총 SN 수율은 $0.19 \pm 0.01$ (전체 즉발 중성자의 약 6.6%)로 계산되었다. $^{252}\text{Cf}(sf)$의 경우, 수율은 $0.39 \pm 0.02$ (약 10.4%)이다.
• 모델 검증: 계산된 SN 성분을 (저에너지 데이터에 의해서만 제약된) 증발 모델에 더했을 때, 결과적인 총 PFNS는 $^{239}\text{Pu}$와 $^{252}\text{Cf}$ 모두에서 측정된 고에너지 수율을 정확하게 재현한다. 반면, 증발 전용 모델은 95% 신뢰 구간 내에서도 고에너지 수율을 체계적으로 과소평가한다.
• 각도 안정성: 본 연구는 신뢰할 수 있는 점근적 특성이 $110^\circ \leq \theta \leq 150^\circ$의 각도 범위 내에서만 추출될 수 있음을 확인하였다. 이 범위를 벗어나면 과도적 간섭 효과와 경계 반사로 인해 현재의 시뮬레이션 영역 내에서 수율이 포화되지 않는다.

의의 및 주장
본 논문은 완전한 미시적 프레임워크를 사용하여 기존에 측정된 즉발 핵분열 중성자 스펙트럼 내에서 절단 중성자의 징후를 최초로 식별했다고 주장한다. TDDFT가 사전 가정 없이 자연스럽게 SN을 생성한다는 것과, 이러한 SN이 증발 전용 모델과 고에너지 실험 데이터 사이의 불일치를 해결한다는 것을 입증함으로써, 본 연구는 즉발 핵분열 내에 무시할 수 없는 SN 성분이 존재한다는 직접적인 증거를 제시한다.

저자들은 이 작업이 표준 통계적 증발 개념을 넘어 핵분열 중성자 스펙트럼을 해석하기 위한 미시적 토대를 구축한다고 상정한다. 또한, 특정 실험적 징후를 식별하였다: 에너지 임계값($E_{thr}$) 아래에서 SN과 EN 기여를 분리하는 것이다. 이는 지정된 각도 범위에서의 저에너지 측정이 증발 성분에 대한 가장 강력한 제약 조건을 제공하는 반면, 고에너지 꼬리는 SN 기여의 주요 지표 역할을 함을 의미한다. 본 연구는 현재의 시뮬레이션이 특정 각도 범위로 제한되어 있지만, 방법론적으로 거의 한 세기 동안 이론적으로 가설로만 존재했던 절단 중성자의 측정 가능한 징후를 성공적으로 분리해 냈다고 결론짓는다.