[논문 리뷰] Monte-Carlo simulation for fragment mass and kinetic energy distributions from neutron induced fission of 235U

M. Montoya, E. Saettone외 1명|ArXiv.org|2007. 09. 07.

Nuclear Physics and Applications참고 문헌 2인용 수 4

한 줄 요약

이 연구는 열성자 유도 $^{235}$U 열분열에서 파편 질량 및 운동에너지 분포를 몬테카를로 시뮬레이션을 통해 모델링하며, 중성자 방출 및 초기 파편 수율과 운동에너지의 변동성을 고려한다. 시뮬레이션은 최종 파편 운동에너지 표준편차의 실험적으로 관측된 두 가지 넓어짐을 성공적으로 재현한다: 질량 $m \approx 109$에서의 넓어짐과 $m \approx 125$에서의 넓어짐으로, 각각 중성자 방출 역학과 초기 파편 평균 운동에너지 곡선의 비연속성에 기인한다.

ABSTRACT

The mass and kinetic energy distribution of nuclear fragments from thermal neutron induced fission of 235U have been studied using a Monte-Carlo simulation. Besides reproducing the pronounced broadening on the standard deviation of the final fragment kinetic energy distribution $σ_{e}(m)$ around the mass number m = 109, our simulation also produces a second broadening around m = 125, that is in agreement with the experimental data obtained by Belhafaf et al. These results are consequence of the characteristics of the neutron emission, the variation in the primary fragment mean kinetic energy and the yield as a function of the mass.

연구 동기 및 목표

• 열성자 유도 $^{235}$U 열분열에서 최종 파편 운동에너지 표준편차 ($\sigma_e(m)$)의 실험적으로 관측된 넓어짐의 기원을 이해하는 것.
• Belhafaf 등이 관측한 $m \approx 125$에서의 넓어짐이 초기 파편 성질에 기인하는지 아니면 분열 후 역학에 기인하는지 조사하는 것.
• 중성자 방출을 통해 초기 파편 분포(질량 수율 $Y(A)$, 운동에너지 $E(A)$)에서 최종 파편 분포($Y(m)$, $e(m)$)로의 변환을 모델링하는 것.
• 중성자 방출 특성과 초기 파편 운동에너지 및 수율의 변동성이 최종 운동에너지 분포 형성에 미치는 역할을 규명하는 것.

제안 방법

• 초기 파편 질량 및 운동에너지 분포를 중성자 방출을 거쳐 최종 파편 분포로 전파하기 위해 몬테카를로 시뮬레이션을 사용하였다.
• 초기 파편 운동에너지 $E(A)$는 질량에 따라 달라지는 평균과 표준편차를 가진 가우시안 분포로 모델링되었다.
• 중성자 방출은 실험 데이터에서 유도한 질량에 따라 달라지는 평균 중성자 수 $\bar{\nu}(A)$를 사용하여 시뮬레이션되었으며, 이는 매개변수 $\alpha$와 $\beta$를 통해 조정되었다.
• 최종 파편 운동에너지 $e(m)$는 에너지 보존 법칙을 적용하여 계산되었으며, 중성자 방출에 의한 운동에너지 손실를 고려했다.
• 최종 질량 $m$에 대한 최종 운동에너지 분포의 표준편차 $\sigma_e(m)$를 계산하여 실험 데이터와 비교하였다.
• 감도 분 析를 위해 매개변수 $\alpha$와 $\beta$를 변화시켜 $\sigma_e(m)$의 넓어짐의 형태와 크기에 미치는 영향을 평가하였다.

실험 결과

연구 질문

• RQ1 왜 $^{235}$U 열분열에서 최종 파편 운동에너지 분포의 표준편차($\sigma_e(m)$)가 $m \approx 109$에서 뚜렷한 넓어짐을 보이는가?
• RQ2 왜 $\sigma_e(m)$에서 관측된 $m \approx 125$에서의 두 번째 넓어짐은 이전의 시뮬레이션으로는 재현되지 않았는가?
• RQ3 중성자 방출 특성과 초기 파편 평균 운동에너지 $\bar{E}(A)$의 변동성이 $\sigma_e(m)$의 관측된 넓어짐에 어느 정도 기여하는가?
• RQ4 간단한 모델이 중성자 방출과 초기 파편 수율/에너지 변동성을 포함하여 전체 실험적 $\sigma_e(m)$ 곡선, 특히 $m \approx 125$ 특징까지 재현할 수 있는가?
• RQ5 질량 $m \approx 130$ 근처에 세 번째 넓어짐이 존재하는 증거가 있으며, 만약 그렇다면 이를 설명할 수 있는 물리적 메커니즘은 무엇인가?

주요 결과

• 시뮬레이션은 $m \approx 109$에서 관측된 뚜렷한 $\sigma_e(m)$ 넓어짐을 성공적으로 재현하였으며, 이는 중성자 방출과 초기 파편 운동에너지 분포의 형태에 기인함을 확인하였다.
• 모델은 $m \approx 125$에서의 두 번째 넓어짐도 재현하였으며, 이는 $A = 126$와 $A = 125$ 사이에서 초기 파편 평균 운동에너지 곡선 $\bar{E}(A)$의 비연속성에 기인함을 확인하였다. 이는 실험 데이터와 일치한다.
• $\sigma_e(m)$의 넓어짐의 형태와 크기는 중성자 방출 모델의 매개변수 $\alpha$와 $\beta$에 민감하며, 실험과의 최적 일치를 위해 $\alpha = 0.62$와 $\beta = 0.35$로 설정되었다.
• 최종 파편에 대한 시뮬레이션된 평균 중성자 수 $\bar{\nu}(m)$는 초기 파편의 $\bar{\nu}(A)$와 비교해 질량 이동에 따라 이동함을 보였다.
• $m \approx 125$에서의 넓어짐은 초기 운동에너지 분포의 내재적 넓어짐 때문이 아니라, 중성자 방출과 $Y(A)$ 및 $\bar{E}(A)$의 기울기 변화 간의 상호작용에 기인한다.
• 시뮬레이션은 $m \approx 130$ 근처에 가능한 세 번째 넓어짐이 존재할 수 있음을 시사하며, 이는 $m = 129$에서 실험적 $\sigma_e = 3.9$ MeV 값을 외삽하여 지원된다.