祝锦 彭肖宇 罗思远 肖万成 何列 刘雨晨 罗凤姣 肖敏 王晓冬

（南华大学 核科学技术学院 衡阳 421001）

重离子超导同步加速器（Nuclotron-based Ion Collider fAcility，NICA）是俄罗斯杜布纳联合核子研究所（Joint Institute for Nuclear Research，JINR）正在建造的新型对撞机。它能提供4～11 GeV 能量范围内的Au-Au 重离子碰撞，以此来寻找高重子密度区夸克物质和强子物质的混合相，以及研究量子色动力学（Quantum Chromodynamics，QCD）物质的相变临界点［1］。NICA 的重离子实验计划十分有利于把此前的相关研究拓展到高重子密度区。尽管第一性原理计算无法用于该能量域，但是理论预期、高密度下的QCD相图可能具有丰富的结构［2］。

NICA加速环上共设置了两个对撞点，两个对撞点上的探测器可以同时工作，多功能探测谱仪（Mult-Purpose Detector，MPD）是位于其中一个对撞点上的探测器［3］，可用于研究重离子碰撞中高温高密核物质的特性，特别是用于研究夸克禁闭、手征对称破缺、核物质相变及相变临界点等科学问题。电磁量能器（Electromagnetic Calorimeter，Ecal）作为MPD 上重要的子探测器，其主要用于探测10 MeV至几个GeV 能域内电子和光子的能量、时间、位置信息［4］。为了满足实验的要求，MPD-Ecal 必须具备莫里哀半径较小的致密工作物质和最小的电磁级联簇射重叠的优点，并且能在MPD 螺线管的磁场（B=0.5 T）中运行。为了避免粒子入射角度对探测器的性能造成影响和减少探测器的探测死区，Ecal 采用物理性能优异的“Shashlik”型抽样量能器，并且采取圆桶状全对心的设计方案［5］，将其模块制作为楔形。

“Shashlik”型电磁量能器凭借稳定的性能被广泛应用于多个大科学装置上，例如PHENIX［6］、KOPIO［7］、LHCB［8］、T2K［9］等实验。能量分辨率、时间分辨率、位置分辨率是评估电磁量能器性能的重要依据。MPD-Ecal 的技术报告（Technical Design Report，TDR）［10］中要求探测器在3 GeV 电子入射条件下，量能器必须满足好于5%的能量分辨率、小于1 ns的时间分辨率（最好能够接近100 ps）、优异的位置分辨率的条件。目前，Ecal的整体结构尺寸［11-12］、吸收体和反射材料的选取［13-14］、波长位移光纤的选择以及光子输运过程［15］已经得到了充分的研究。本工作基于Geant4 软件［16］模拟了多个物理参数对量能器的能量分辨率、时间分辨率的影响，包括电子入射点位置、tower 中闪烁体层数、入射粒子能量和类别、光纤端面抛光度等，以及使用电荷重心法对7×7阵列组合的tower 的位置分辨率进行了研究。模拟所得到的能量分辨率、时间分辨率、位置分辨率参数将为后续的实验提供参考。

1 电磁量能器结构

JINR对MPD-Ecal的整体设计结构如图1所示，内径为3.42 m，外径为4.226 m，长度为6 m，重量为60 t。整个Ecal由2 400个module组成，根据倾斜角度不同，分为8 个不同区域的module［17］。中方总负责768个module的研制与测试，分别由清华大学、山东大学、复旦大学、南华大学共同完成。南华大学负责研制图1中第一区域的77个module，每个module由16 个tower 通过8×2 的排列方式粘接在一起组成。tower 经过铣切加工呈棱台形，向底部逐渐变小，端部面积为33 mm×33 mm，长度为415.5 mm。每个tower中穿插16根由可乐丽（Kuraray）公司生产的Y-11（200）［18］型号的波长位移光纤（Wavelength-Shifted Fibers，WLSF），光纤由三层不同材料组成。16 根光纤在tower 的末端集束成为一个5.5 mm×5.5 mm的水平端面，距离光纤端面约0.3 mm处耦合一个大小为6 mm×6 mm的硅光电倍增器（SiPM）以读出WLSF的光信号，后接电子学系统，获取入射粒子的能量、时间、位置等物理信息的数据。

tower 结构如图2（a）所示，每个tower 由上下两片厚度分别为7 mm和8 mm的底座及211片塑料闪烁体片（厚度1.5 mm）和210 片铅片（厚度0.4 mm）交错堆叠组成，其横截面积均为40 mm×40 mm，并且铅片的两侧均喷涂有二氧化钛作为反射层。闪烁体片与铅片均留有18个直径为1.6 mm的小孔，其中按照4×4 排列（间距9.6 mm）的16 个小孔用于安装WLSF，剩余的2 个用于安置固定钢丝。同时，为了辅助定位，闪烁体设计成“LEGO”体结构。其表面有4个直径为2 mm且高度为1.25 mm的圆柱突起，背面是直径为2.1 mm 且深度为0.9 mm 的凹槽，相应的在铅片相同位置有4 个直径为2.2 mm 的圆孔。南华大学所研制的1 232个tower的闪烁体层数分布如图2（b）所示，闪烁体层数范围在207～217，其中闪烁体层数为211的tower数目最多。

图2 电磁量能器单元 (a) Tower结构示意图，(b) 闪烁体层数分布Fig.2 Electromagnetic calorimeter unit (a) Schematic diagram of Tower structure, (b) Distribution of scintillator layer

2 电磁量能器模拟建模

本工作基于Geant4.10.07.p03版本，对影响电磁量能器性能的多个参数进行模拟研究。粒子源分别设置为单能电子和天然宇宙射线源（Cosmic Ray，CRY），表1为模拟中使用的塑料闪烁体与反射材料的主要性能参数。程序中所建立的module 模型如图3（a）所示，模型中的反射材料均为TiO2。闪烁体选取的是北京高能科迪科技有限公司生产的HNDS2型塑料闪烁体［19］，该塑料闪烁体的发光衰减曲线如图3（b）所示。

图3 电磁量能器工作原理 (a) Module模型，(b) 发光衰减曲线，(c) 电磁簇射简化模型，(d) 光子传输过程Fig.3 Working principle of electromagnetic calorimeter (a) Module model, (b) Luminescence decay curve, (c) Simplified model of electromagnetic shower, (d) Photon transmission process

电磁簇射过程即高能电子在tower 中经过一段路程后发生轫致辐射，产生次级电子和轫致辐射光子，后者在经过一段路程后发生电子对效应，又会转换为正负电子对。此过程交替出现，直至粒子能量下降到临界值以下时，簇射粒子停止增殖。对于电磁簇射的发展，通常使用辐射长度（X0）进行描述。X0表示电子能量减少至原来的1/e所走过的射程，图3（c）所示为电磁簇射简化模型。tower 的长度为415.5 mm，等效于电子经过11.8X0。塑料闪烁体内沉积的能量使闪烁体发光（波长为440 nm），产生的光子大部分直接进入WLSF，小部分经过反射层反射后进入WLSF，光子进入后将以470 nm 的波长再次发射，被光纤一端的SiPM收集，图3（d）表述了光子在塑料闪烁体与光纤中的传输过程。

3 电磁量能器模拟

3.1 电子入射位置对module能量分辨率的影响

module 的能量分辨率表征了它对相近能量的辐射粒子的分辨本领，探测器的能量分辨率越小，其对粒子能量分辨的精度越高，蒙特卡罗模拟中通常将能量的相对涨落σE/E作为量能器的能量分辨率的表征［20］。模拟设置3 GeV 的单能电子从不同tower的中心径直入射，来研究入射位置对module能量沉积和能量分辨率的影响，对module 一侧的tower 进行编号如图4（a）所示，一共有8 个入射点。粒子入射位置对module 的闪烁体内能量沉积影响结果如图4（b）所示，2～7号tower中心位置入射时，相比于1号和8 号位置入射，闪烁体内的能量沉积较多。对于闪烁体层数为211 层的tower，随着电子入射位置由module 边缘向中心移动，闪烁体内能量沉积由718 MeV增加至758 MeV。图4（c）为电子入射位置对module 能量分辨率的影响，当入射位置为2～7 号tower 中心位置时，module 能量分辨率较好，且可好于4.35%。入射位置接近module 边缘时，module 周围存在较大的横向能量泄漏，导致module的闪烁体内能量沉积少，使module 的能量分辨率变差。然而，MPD-Ecal 是桶型结构，因此不存在大量横向能量泄漏，所以推断根据Ecal 总体结构模拟给出的能量分辨率要好于4.35%。

图4 入射位置对module性能的影响 (a) 8个不同电子入射位置，(b) 入射位置对能量沉积的影响，(c) 入射位置对能量分辨率的影响Fig.4 Influence of incident position on module performance (a) Eight different electron incident positions, (b) Influence of incident position on energy deposition, (c) Influence of incident position on energy resolution

3.2 闪烁体层数对module能量分辨率的影响

为了研究闪烁体层数对module 能量分辨率的影响，模拟中设置3 GeV 电子入射不同闪烁体层数的module，入射点位置为4 号tower 中心。基于tower 长度为415.5 mm、tower 中的闪烁体层数为207～217 片、闪烁体厚度为（1.50±0.01） mm、铅片厚度为（0.39±0.04） mm条件下，tower中的闪烁体层数与材料（铅片、闪烁体片）厚度之间的对应关系如图5（a）所示。图5（b）为闪烁体层数对module 的能量沉积的影响结果，闪烁体层数与能量沉积之间具有良好的线性关系。闪烁体的层数每增加1 层，闪烁体内的能量沉积增加9 MeV，铅片内的能量沉积减少20 MeV。图5（c）展示了闪烁体层数对module的能量分辨率的影响，随着闪烁体的层数增加，module能量分辨率变差。当闪烁体的厚度为1.49 mm、层数为209 时，module 能量分辨率达到最佳，可好于4.15%。Tower 层数在207～217 区间内，其组成的module 能量分辨率均能好于4.79%，模拟结果好于TDR（Ecal）中5%的性能要求。

图5 闪烁体层数对module性能的影响 (a) 闪烁体层数与铅片厚度的关系，(b) 闪烁体层数对能量沉积的影响，(c) 闪烁体层数对能量分辨率的影响Fig.5 Effect of scintillator layer number on module performance (a) Scintillator layer vs. thickness of lead, (b) Effect of scintillator layer on energy deposition, (c) Effect of scintillator layer on energy resolution

以上的模拟过程中同时改变了闪烁体的层数和铅片厚度两个参数，为了进一步确定层数和铅片厚度对module 能量分辨率的影响。在模拟参数为1.5 mm 厚闪烁体片、0.025 mm 厚反射层的条件下，对不同层数的tower的能量分辨率进行模拟研究；以及对相同层数的tower，通过改变铅片的厚度来研究铅片厚度对tower 的能量分辨率的影响。模拟结果如图6 所示，闪烁体层数或铅片厚度的增加都能使module 能量分辨率变好。这是因为闪烁体层数和铅片厚度增加都能增长tower的等效辐射长度，有利于tower 中的电磁簇射事件的发展。闪烁体层数在208～216 区间内时，铅片厚度比闪烁体层数对module 能量分辨率的影响要大。因此，在tower 长度不变的条件下增加闪烁体层数，module的能量分辨率变差。

图6 材料参数与能量分辨率的关系 (a) 闪烁体层数与能量分辨率的关系，(b) 铅片厚度与能量分辨率的关系Fig.6 Material parameters vs. energy resolution (a) Relationship between the number of scintillator layers and energy resolution,(b) Relationship between the thickness of lead and energy resolution

3.3 闪烁体层数对tower时间分辨率的影响

Tower 的时间分辨率表征其对不同时刻到达探测器的辐射粒子的区分能力，时间分辨率越小，所探测到粒子的时间信息越准确。图7（a）为模拟方案，设置电子径直入射tower 中心，在tower 的前方放置两个大小为40 mm×40 mm 的塑料闪烁体（A、B）作为触发探测器。在模拟中的定时方法借助了过零定时［21］的定义，使用信号的峰值时刻t0作为SiPM探测到光电子的时刻。SiPM 探测到的光电子的时间信息分布服从朗道分布，使用root 程序［22］对其进行朗道拟合，再对拟合曲线寻峰获取峰值时刻t0，如图7（b）所示。再通过式（1）计算可得tower 的时间分辨率σT（t），式中，σT（t）-T（a）、σT（t）-T（b）、σT（a）-T（b）为探测器两两之间的时间差分布的标准差。

图7 时间分辨率模拟 (a) 模拟方案示意图，(b) SiPM探测到光电子的时间分布Fig.7 Time resolution simulation (a) Schematic diagram of the simulation scenario, (b) Time distribution of SiPM detected photoelectrons

图8（a）展示了在3 GeV 单能电子入射条件下，闪烁体的层数对SiPM探测到光电子数和tower时间分辨率的影响。随着闪烁体层数增加，探测到的光电子数目增多，tower的时间分辨率变好。闪烁体厚度为1.51 mm、层数为215时，tower时间分辨率可达最好，且小于96 ps。Tower层数在207～217区间内，其时间分辨率均能小于118 ps，模拟结果好于TDR对tower 的时间分辨率小于1 ns 的要求。模拟得到的时间分辨率是tower在理想模型下的理论性能，该结果可以作为性能的参考指标，从而对探测器进行结构上的改进或优化。但是，模拟中并未考虑电子学器件对粒子时间信息的测量影响，所以与实验测量值存在一定差异。

图8 时间分辨率模拟结果 (a) 闪烁体层数对tower时间分辨率的影响，(b) 入射电子能量对tower时间分辨率的影响Fig.8 Simulation results of time resolution (a) Effect of scintillator layers on tower time resolution, (b) Effect of incident electron energy on tower time resolution

使用不同能量的电子径直入射tower中心，入射电子能量对SiPM探测到光电子数和tower时间分辨率的影响结果如图8（b）所示。随着电子能量增加，SiPM 探测到光电子数目增多，tower 的时间分辨率变好。探测到的光电子数目越多，SiPM对信号的峰值时刻分辨越精确。因此入射电子能量在0.1～1 GeV阶段，随着入射电子能量增加，探测到的光电子数目增多，tower时间分辨率变好趋势明显。当入射电子能量为1 GeV 时，探测到的光电子数目为2 903，此时SiPM 已经能够很好分辨信号的峰值位置。因此在1～3.5 GeV 阶段，随着入射电子能量增加，探测到的光电子数目增多，tower 时间分辨率变好的趋势较缓慢。对于闪烁体层数为211层的单个tower，在3.5 GeV电子入射条件下，其时间分辨率小于110 ps。

3.4 光纤端面抛光度对tower时间分辨率的影响

NICA 实验中桶型电磁量能器采取的是单端读出方案，桶外层的光纤端部按照tower各自集束后耦合SiPM，桶内壁的光纤端面则通过涂覆反射材料，以提高每个tower的光子输出。由于SiPM是直接耦合在光纤端部，所以耦合端的截面需要尽可能光滑以免损伤SiPM。涂覆反射材料的非耦合端面的抛光度根据加工工艺不同，则有多种选择。端面反射层的抛光度会影响SiPM光电子的收集效率，进而影响tower的时间分辨率，因此本小节研究了涂覆反射材料的光纤端面不同抛光度对SiPM 探测到光电子数和tower 时间分辨率的影响。蒙特卡罗模拟中抛光度为镜面反射在反射中所占的比值，0 表示光纤端面的光子全发生漫反射，1 表示光纤端面的光子全发生镜面反射。

粒子源设置为3 GeV 单能电子，端面抛光度对探测到光电子数的影响如图9（a）所示，对tower 时间分辨率的影响结果如图9（b）所示。光纤端面抛光度越高，SiPM探测到的光电子数越多，tower的时间分辨率越好。抛光度为0 时，探测到的光电子数目为8 331个，tower的时间分辨率小于126 ps；抛光度为1时，探测到的光电子数目为9 484个，tower时间分辨率小于103 ps。两者的结果存在差异，是因为光子在波移光纤内传输时，漫反射会使光子在光纤内的反射次数增多，使光子在传输过程中损失增大，探测到的光电子数少，导致tower 时间分辨率变差。当端面抛光度为0.6时，光纤端面涂抹反射材料相比于不涂抹反射材料的tower，SiPM 探测到光电子数提升了约42%。

图9 抛光度对探测器性能的影响 (a) 光纤端面抛光度对SiPM探测到的光电子数目的影响，(b) 光纤端面抛光度对tower时间分辨率的影响Fig.9 Effect of polishing on detector performance (a) Influence of optical fiber end face polishing on the number of photoelectrons detected by SiPM, (b) Influence of optical fiber end face polishing on the time resolution of tower

3.5 天然宇宙射线模拟tower时间分辨率

本工作拟搭建一套基于宇宙射线的电磁量能器性能测试平台，为了更好地优化测试条件和性能，在模拟中利用了天然缪子源程序包［23-24］（CRY，平均能量为3～4 GeV）研究单个tower（211 层）的时间分辨率。研究方案如图10（a）所示，tower置于探测器A、B 之间，三个探测器中两两做触发。图10 中（b）、（c）、（d）、（e）为模拟结果，计算得到tower 的时间分辨率为185 ps。宇宙射线测试tower 相比于3 GeV单能电子测试，tower 时间分辨率变差了73 ps。造成该结果是因为缪子质量约为电子的207 倍，具有极强的穿透能力；与电子在闪烁体中通过轫致辐射损失能量引发电磁簇射过程不同，缪子穿过物体通过电磁相互作用损失能量，穿越闪烁体时损失能量使闪烁体原子发生激发—退激作用产生荧光。因此缪子穿过tower时在塑料闪烁体内沉积的能量少，导致探测到的光电子数目少，tower 的时间分辨率变差。

图10 宇宙线测试 (a)宇宙线测试模拟方案，(b) 光电子产额，(c) Tower与A的时间差分布，(d) Tower与B的时间差分布，(e) A与B的时间差分布Fig.10 Cosmic ray test (a) Simulation scheme for cosmic ray testing, (b) Photoelectric yield, (c) Distribution of time difference between Tower and A, (d) Distribution of time difference between Tower and B, (e) Distribution of time difference between A and B

3.6 电磁量能器的位置分辨率模拟

Ecal的位置分辨率表征了探测器对粒子入射位置的分辨能力，一般通过位置重建的方法来获取粒子入射位置，位置分辨率越小，探测器对粒子位置重建的准确度越高。Ecal 是由大量tower 组成的桶型电磁量能器，本工作将tower组成的阵列近似看作量能器的组成结构，以此来模拟Ecal的位置分辨率，图11（a）为7×7 阵列组合tower 在横向上的排列方式。在3 GeV电子径直入射tower阵列的正中心条件下，tower 阵列组合的闪烁体内能量沉积分布图如图11（b）所示，能量沉积主要集中在中心入射的tower内。由于tower横向面积小，当电子入射到某个tower时，相邻区域的tower也会响应，且距离入射点位置越远的tower 能量沉积越少。本工作利用电荷重心法对入射tower 及周围8 个tower 的能量沉积进行分析，假设各个方位具有相同权重，在此假设条件下重建了入射粒子的入射位置。式（2）是电荷重心法的计算方法，其中Ei为水平方向上第i个tower 的能量沉积，Xi为第i个tower 的水平中心位置，Yi为第i个tower的垂直中心位置。

图11 位置分辨率模拟 (a) 7×7的tower组合，(b) 阵列tower的能量沉积分布Fig.11 Coordinate resolution simulation (a) 7×7 tower combination, (b) Energy deposition distribution of array towers

由于最外层tower 的边缘没有tower 分布，无法对该区域内的射入粒子进行位置重建。设置3 GeV电子随机打入中间区域（2≤i≤6，2≤j≤6）内的tower，利用电荷重心法得到了水平方向重建位置与入射位置的对应关系如图12（a）所示。其中，点为利用电荷重心法计算出的电子位置，直线为入射粒子实际位置，二者的差别表现出粒子重建位置与入射位置存在有规律的波动偏差。3 GeV 电子打入tower 正中心时，入射点周围3×3 区域内水平方向上的能量沉积分布如图12（b）所示，入射点附近10 mm 范围内沉积了约80%的能量。因此，当入射点位置由中心向两边位移的幅度较小时（Δx<10 mm），不会较大影响到每个tower内的能量沉积，则通过电荷重心法得到的重建位置不会发生较大改变，这将导致距离中心点越远重建精度越差。随着位移幅度的增大（10 mm<Δx<20 mm），靠近移动方向tower内的能量沉积将会快速增加，此时探测器的位置重建精度随着位置幅度增大逐渐恢复。因此当入射点位置由tower中心向边缘移动时，量能器的位置重建误差会出现先增大后减小的正弦波动现象，造成该现象的主要原因来源于位置重建方法本身。为了消除正弦波动特性的影响，使用正弦函数对其进行拟合如图12（c）所示，利用拟合后的函数曲线对重建位置进行校正，从而减少由电子入射位置导致的重建位置偏差。校正后的结果如图12（d）所示，与图12（a）相比，它的重建位置与入射位置具有更好的线性关系。

图12 电子重建位置校正 (a) 电子位置重建，(b) 入射点周围能量沉积分布，(c) 重建位置误差与入射位置的关系，(d) 校正后的电子重建位置Fig.12 Correction of the reconstructed position of the electrons (a) Electronic position reconstruction, (b) Energy deposition distribution around the point of incidence, (c) Reconstruction position error vs. incident position, (d) Calibrated electronically

对校正后的重建位置偏差分布进行高斯拟合，所得到的分布宽度σX即Ecal 的位置分辨率。设置3 GeV电子随机打入tower上，闪烁体层数对探测器位置分辨率的影响如图13（a）所示。随着tower 中闪烁体层数的增加，Ecal的位置分辨率变差，但是变化的趋势非常缓慢，这是因为闪烁体层数增加虽然增加了闪烁体内的能量沉积，但是使tower内的总能量沉积减少了，导致电子位置重建误差变大。当闪烁体的厚度为1.49 mm、层数为209 时，Ecal 位置分辨率达到最佳为1.84 mm。在兼顾电磁量能器的各项性能后，tower 中的闪烁体层数最佳为211 层、厚度为1.5 mm。此时，Ecal的能量分辨率好于4.35%，时间分辨率小于112 ps，位置分辨率好于2 mm。

图13 位置分辨率模拟结果 (a) 闪烁体层数对位置分辨率的影响，(b) 入射电子能量对位置分辨率的影响Fig.13 Simulation results of coordinate resolution (a) Effect of scintillator layer number on coordinate resolution, (b) Effect of incident electron energy on coordinate resolution

为了研究不同能量下Ecal 的位置分辨能力，设置不同能量的电子随机打入到tower上，以此来模拟入射电子能量对Ecal 位置分辨率的影响，其结果如图13（b）所示。随着电子能量提高，Ecal的位置分辨率变好。电子能量在0～1 GeV，随着能量增加，Ecal的位置分辨率变好趋势快；在1～3.5 GeV 阶段，Ecal的位置分辨率变好趋势缓慢。造成该结果是因为随着电子能量增加，周围的tower 内能量沉积增多，有利于提高位置重建精度。当电子能量达到1 GeV以后，周围tower 沉积的能量足以用于对电子位置重建。此时，Ecal 的位置分辨率随电子能量变化的趋势减缓。入射电子能量为3.5 GeV 时，探测器的位置分辨率好于1.85 mm。

4 结语

本工作基于Geant4 软件，使用单能电子模拟研究了多个参数对电磁量能器性能的影响，包括电子入射点位置和能量、闪烁体和铅片的层数和厚度、光纤端面抛光度。在tower长度为415.5 mm限定条件下，当3 GeV 电子沿着module 中心区域入射时，module 的横向能量泄漏最少，且随着tower 中闪烁体层数增加，module 能量分辨率变差，tower 时间分辨率变好，Ecal 位置分辨率变差。在兼顾能量分辨率、时间分辨率、位置分辨率后，tower中闪烁体层数最佳为211 层。闪烁体层数在207～217 范围内，module能量分辨率均能好于4.79%、tower时间分辨率小于118 ps，Ecal 位置分辨率好于2.2 mm。随着入射电子能量增加，量能器的时间分辨率和位置分辨率变好。当光纤端面抛光度为0.6时，光纤端面涂覆反射层与不涂覆相比，SiPM探测到的光电子数目增加了42%。随着光纤端面抛光度提高，探测到的光电子数目增加，tower时间分辨率变好，抛光度为1时，时间分辨率小于103 ps。为了给后续宇宙线测试实验提供一些参考，本工作还使用了天然宇宙射线程序包模拟了单个tower的时间分辨率，得到的结果小于185 ps。模拟得到的结果是tower 在理想模型下的理论性能，该结果可以作为性能的参考指标，从而对探测器进行结构上的改进或优化。但是，模拟中并未考虑电子学器件对粒子时间信息的测量影响，所以将会与实验测量值存在一定差异。

作者贡献声明祝锦负责研究的提出及设计、数据的收集和整理、文章的起草和最终版本的修订；彭肖宇负责模拟程序设计；罗思远负责数据的收集；肖万成负责数据的收集；何列负责数据的验证与核实；刘雨晨负责模拟结果可视化；罗凤姣负责研究课题的指导；肖敏负责研究课题的指导；王晓冬负责最终版本的修订、项目的监督和管理。

1Taranenko A. Status of the mega-science project NICA[J].Journal of Physics: Conference Series, 2020, 1685(1):012021. DOI: 10.1088/1742-6596/1685/1/012021.