基于蒙特卡洛模拟评估一种测量电荷改变截面的实验技术

《Particles》：Evaluation of an Experimental Technique for Measuring Charge-Changing Cross-Sections via Monte Carlo Simulations

　　

测量电荷改变截面（charge-changing cross-section）作为一种确定质子半径的方法已被发展起来，尤其适用于不稳定、短寿命原子核。为获得精确的电荷半径，此类截面必须以高精度测量。然而实验方法中存在复杂性，导致在确定精确核半径时引入不确定度。因此，需要能够描述带电粒子相互作用复杂物理过程的良好模型，以验证实验方法并估算系统不确定度的贡献。GEANT4是一款广泛用于描述从原子物理到宇宙线能区重离子碰撞中相互作用的蒙特卡洛模拟程序。研究人员利用10,12C在不同次级靶上的电荷改变反应实验测量数据，在GEANT4中尽可能真实地重现了实验探测器几何、束流剖面与探测器配置。研究人员将实验获取能谱与相应GEANT4模拟结果比较，以验证对所测能谱的解释并评估系统效应；同时推导次级粒子产额比并与GEANT4结果对照。

研究背景方面，放射性离子束电荷改变截面（σ_cc）测量是探究远离稳定谷不稳定核结构的有力工具，该观测量直接关联核内质子半径分布，可揭示中子晕、中子皮等奇异现象。提取可靠核结构信息要求截面测量不确定度达到或优于1%，这使得对系统效应、探测器响应及反应动力学的精细理解至关重要。现有研究表明，高能区（>1.5 GeV/u）反应截面模型总体趋势尚可，但中能区（300–500 MeV/u）σ_cc能量依赖显著、模型预测分歧大，且缺乏该能区实验数据。GSI的BARB计划旨在探索正电子发射碳氧同位素用于重离子治疗，本研究聚焦其中实验方法验证而非治疗应用。论文发表于《Particles》。

关键技术方法上，研究人员依托GSI的FRS（FRagment Separator）在线碎片分离器，以470–480 MeV/u¹²C初级束轰击Be靶产生¹⁰C次级束，在F2与F4焦平面分别布置反应靶与探测阵列；采用透射法（transmission technique）由靶有无条件下的同电荷态出射粒子计数比推σ_cc；利用LISE++模拟空间分离同位素，以MUSIC（Multiple-Sampling Ionization Chamber）、塑料闪烁体（SCI）、TPC（Time Projection Chamber）、BaF₂晶体望远镜作粒子鉴别与能量测量；构建GEANT4 11.4.0仿真，采用QBBC_ABLA物理列表（QBBC处理强子级联、ABLA处理退激发），植入实测束流相空间与探测器几何及响应参数，但不含电子学数字化。

研究结果部分，第二章实验方法指出，FRS标准消色差模式下，F2处经磁刚度选¹⁰C、±20 mm狭缝压制¹¹C，上游塑料闪烁体定Z；F2下游MUSIC+SCI鉴别出射电荷，F4处用ΔE-TOF-MUSIC作逐事件A/Q与Z识别，BaF₂测残余动能。第三章GEANT4模拟实验复现F2（TPC 24×6×7 cm³、MUSIC有源22×10×30 cm³、SCI厚0.49 cm）与F4（六边形塑料闪烁体、BaF₂六棱柱）布局，束流均值方差取自实测，电磁过程走Bethe–Bloch，光学光子产额等按参数植入但截面提取仅依赖能损与电荷鉴别。第四章结果与讨论中，图5（F2的¹⁰C束MUSIC-SCI二维ΔE关联）显示靶入时碳主线下方出现B(Z=5)、Be(Z=4)碎裂带，靶出时极弱，模拟重复该结构且分辨率更优；图6（F4的¹²C束SCI-ΔE vs BaF₂ E）同样见Z=6主线及轻核对角带。研究人员用两种计数法（下游探测器围线计数、靶后紧邻计数）从模拟提σ_cc：¹⁰C@F2得637 mb与638 mb，¹²C@F4得1232 mb与1233 mb，差异≈0.1%，小于目标精度，证明透射法受输运与探测器效应引入的系统差极小。次级粒子产额比方面，Z≥6不变粒子占比与Z=5（硼）占比实验与模拟吻合；¹²C入射时Z=5约占总碎裂7%，¹⁰C入射时仅约0.2%。

讨论与结论部分翻译：综合F2的MUSIC-SCI与F4的塑料-BaF₂关联实验分析及模拟可见，未反应束与碎裂产物鉴别清晰一致；靶入构型中碎裂组分显著证实源自靶内核相互作用，靶出与靶入共有展宽归于束流本征、空气隙/探测器物质作用及探测器响应。GEANT4以与实验相同流程分析并与模拟真值比对，偏差<0.1%，故该方法适于<1%精度σ_cc测量；次级粒子产额比亦与GEANT4相符。该研究意义在于，为中能区放射性束σ_cc实验提供经模拟背书的系统不确定度量化框架，支撑后续由σ_cc推质子半径及检验晕核结构的可信度。