[论文解读] CEM2k - Recent Developments in CEM

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

Nuclear reactor physics and engineering参考文献 11被引用 28

一句话总结

本文提出 CEM2k，即核反应级联-爆发模型（Cascade-Exciton Model）的高级版本，旨在改进高能质子与重离子碰撞后残余核素产额的预测。该模型引入了更精确的核态密度、改进的π介子吸收与泡利不相容效应、精确的能量-动量守恒，以及更新的结合能，显著提升了与近期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.

研究动机与目标

• 提升级联-爆发模型（CEM）在中高能核反应中的预测准确性。
• 解决现有模型（如 LAHET、YIELDX）与新近高精度 GSI 实验数据在逆向动力学实验中同位素产额方面的差异。
• 开发新版本代码 CEM2k，以更准确描述加速器驱动系统中粒子发射谱与残余核素产额。
• 改进预平衡与蒸发阶段的处理，特别是从级联到预平衡衰变的过渡判据。
• 解决裂变建模中的局限性，特别是对 238U 等锕系元素，其裂变截面在总反应截面中占主导地位。

提出的方法

• 采用新过渡判据，从核内级联（INC）阶段进入预平衡阶段，使用参数 $\mathcal{P} = |(W_{\text{opt.mod.}} - W_{\text{opt.exp.}})/W_{\text{opt.exp.}}|$，取代旧版本中固定的 $\mathcal{P} = 0.3$ 阈值。
• 引入改进的基底截面、更新的核质量、$Q$-值、配对与结合能，以及校正后的态密度系统学。
• 对每个蒙特卡罗事件实施精确的能量-动量守恒，取代旧版本中使用的统计守恒方法。
• 在级联阶段使用真实的核子结合能，而非恒定的 7 MeV 近似值。
• 在发射宽度计算中，用发射粒子的约化质量替代无反冲近似。
• 在核外引入快级联核子的复合粒子（如 α、锂）形成过程。

实验结果

研究问题

• RQ1 如何改进 CEM 模型，以更准确地再现近期 GSI 实验中逆向动力学反应的高精度同位素产额数据？
• RQ2 对级联到预平衡过渡判据的哪些修改能提升模型与实验数据的一致性？
• RQ3 为何现有代码（如 LAHET 和 YIELDX）在 800–1000 MeV/nucleon 条件下无法准确描述 197Au 和 238U 反应中的实测同位素产额？
• RQ4 对 238U 等锕系靶核，裂变与中子至裂变宽度比（$\Gamma_n/\Gamma_f$）的处理方式如何影响模型预测结果？
• RQ5 改进的守恒律与真实的结合能处理在多大程度上提升了 CEM2k 模型的准确性？

主要发现

• 在 800 MeV/nucleon 条件下，CEM2k 对 197Au + p 反应的 GSI 实验数据表现出显著改进，尤其在汞到铪以及镥到钕的同位素分布区域。
• 在 1 GeV/nucleon 条件下，CEM2k 对 208Pb 的实测同位素产额预测精度优于 CEM97 与 LAHET-Bertini，特别是在富中子区域。
• 在 1 GeV/nucleon 条件下，CEM2k 对 238U + p 反应的中子缺乏区域数据存在低估，表明裂变碎片产额的建模仍不足。
• 在 1 GeV 条件下，p + 238U 反应中裂变占总反应截面的约 80%（1390 mb / 1788 mb），凸显 $\Gamma_n/\Gamma_f$ 在模型预测中的关键作用。
• 238U 实验数据揭示，同位素产额随中子缺乏核素的增加呈强指数下降趋势，与裂变主导的贡献一致。
• 初步计算显示，在 1 GeV 以上能量区域，特别是 p + Hg 反应中，裂变截面被高估，表明裂变模型仍需进一步优化。