[논문 리뷰] CEM2k - Recent Developments in CEM

S. G. Mashnik, Arnold J. Sierk|ArXiv.org|2000. 11. 17.

Nuclear reactor physics and engineering참고 문헌 11인용 수 28

한 줄 요약

이 논문은 고에너지 프로톤 및 중입자 충돌에서 잔류 뉴클리드 수율을 보다 정확히 예측하기 위해 개선된 캐스케이드-익스시톤 모델(Cascade-Exciton Model)인 CEM2k를 제시한다. 이 모델은 정밀한 핵 준위 밀도, 향상된 피온 흡수 및 파울리 금지 원리 적용, 정확한 운동량-에너지 보존, 업데이트된 결합 에너지 등을 포함하여 최근 GSI의 역-역학 실험 데이터와의 일치도를 크게 향상시켰다. 특히 800–1000 MeV/핵자에서 197Au 및 208Pb에 대해 뛰어난 예측 성능을 보였다.

ABSTRACT

Recent developments of the Cascade-Exciton Model (CEM) of nuclear reactions are briefly described. These changes are motivated by new data on isotope production measured recently in "reverse kinematics" at GSI for interactions of 208-Pb and 238-U at 1 GeV/nucleon and 197-Au at 800 MeV/nucleon with liquid 1-H. This study leads us to CEM2k, which is a new version of the CEM code that is still under development. The increased accuracy and predictive power of the code CEM2k are shown by several examples. Further necessary work is outlined.

연구 동기 및 목표

• 중간 및 고에너지 핵 반응에 대한 캐스케이드-익스시톤 모델(Cascade-Exciton Model, CEM)의 예측 정확도를 향상시키기 위해.
• 최근의 고정밀도 GSI 데이터(역-역학 실험에서의 동위원소 생성)와 기존 모델(LAHET, YIELDX 등) 간의 괴리 문제를 해결하기 위해.
• 가속기 구동 시스템에서의 입자 방출 스펙트럼과 잔류 뉴클리드 수율을 보다 정확히 기술하기 위해 새로운 코드 버전인 CEM2k를 개발하기 위해.
• 예비 평형 상태와 증발 단계의 처리를 정교화하고, 특히 캐스케이드에서 예비 평형 붕괴로의 전이 조건 기준을 개선하기 위해.
• 우라늄 같은 작동성 우라늄 동위체(예: 238U)에서처럼 붕괴 단면적이 지배적인 경우의 붕괴 모델링의 한계를 해결하기 위해.

제안 방법

• 이전 버전에서 사용하던 고정된 $\mathcal{P} = 0.3$ 임계값을 대체하기 위해, 내부 핵 캐스케이드(INC)에서 예비 평형 단계로의 전이 기준으로 새로운 기준인 $\mathcal{P} = |(W_{\text{opt.mod.}} - W_{\text{opt.exp.}})/W_{\text{opt.exp.}}|$ 를 도입하였다.
• 향상된 기본 단면적, 업데이트된 핵 질량, $Q$-값, 쌍화 및 결합 에너지, 보정된 준위 밀도 체계를 통합하였다.
• 이전 버전에서 사용하던 통계적 보존 방식을 대체하여, 각 몬테카를로 사건에 대해 정확한 운동량-에너지 보존을 적용하였다.
• 캐스케이드 단계 동안 일정한 7 MeV 근사치를 사용하는 대신 실제 중성자 결합 에너지를 사용하였다.
• 방출 폭 계산에서 반동 없음 근사를 제거하고, 방출 입자의 감소 질량을 적용하였다.
• 핵 외부에서 빠른 캐스케이드 중성자로부터 복합 입자(예: 알파, 리튬)의 응집을 도입하였다.

실험 결과

연구 질문

• RQ1 최근 GSI의 역-역학 반응에서의 동위원소 생성에 대한 고정밀도 실험 데이터를 더 잘 재현하기 위해 CEM 모델은 어떻게 개선될 수 있는가?
• RQ2 캐스케이드에서 예비 평형 상태로의 전이 조건 기준을 어떻게 수정하면 실험 데이터와의 일치도를 높일 수 있는가?
• RQ3 LAHET 및 YIELDX와 같은 기존 코드들이 800–1000 MeV/핵자에서 197Au 및 238U 반응에서 측정된 동위원소 수율을 설명하지 못하는 이유는 무엇인가?
• RQ4 붕괴 및 중성자-붕괴 폭 비율($\Gamma_n/\Gamma_f$)의 처리 방식이 238U와 같은 작동성 우라늄 타겟의 예측에 어떻게 영향을 미치는가?
• RQ5 향상된 보존 법칙과 현실적인 결합 에너지 적용이 CEM2k의 모델 정확도에 어느 정도 기여하는가?

주요 결과

• CEM2k는 800 MeV/핵자에서 197Au + p 반응에 대해 GSI 데이터와 뚜렷한 개선된 일치도를 보였으며, 특히 수은에서 텅스텐, lutetium에서 neodymium에 이르는 동위원소 분포에서 뛰어난 성능을 보였다.
• 1 GeV/핵자에서 208Pb에 대해 CEM2k는 CEM97나 LAHET-Bertini보다 측정된 동위원소 수율을 더 정확히 재현하였으며, 특히 중성자 과잉 영역에서 뛰어난 성능을 보였다.
• 1 GeV/핵자에서 238U + p 반응에 대해 중성자 부족 영역에서는 데이터를 과소평가하고 있어 붕괴 조각 생성 모델링이 부족함을 시사한다.
• 1 GeV에서 p + 238U 반응에서 붕괴가 총 반응 단면적의 약 80%를 차지한다(1788 mb 중 1390 mb), 이는 $\Gamma_n/\Gamma_f$ 가 모델 예측에 있어 핵심적인 역할을 한다는 것을 시사한다.
• 238U 데이터는 중성자 부족 핵으로 향하는 동위원소 수율이 강한 지수 감소를 보이며, 이는 주로 붕괴 기여 때문임을 일관되게 나타낸다.
• 초기 계산 결과에서는 1 GeV 이상에서 붕괴 단면적이 과대 예측되는 경향이 있었으며, 특히 Hg와의 p 반응에서 두드러져 붕괴 모델의 추가 보완이 필요함을 시사한다.