自从20世纪英国物理学家狄拉克预言正电子开始，直至近期RHIC-STAR合作组发现反氦-4原子核^[1]，反物质原子核和引起其的核反应机制研究一直是核物理学的前沿课题。迄今为止，反物质相关物理仍有若干未解之谜，如宇宙中正反物质不对称问题、反粒子-粒子相互作用等。目前，人们利用高能核-核碰撞可以产生大量的反质子，再通过冷却储存环收集反质子进行核反应^[2]，欧洲核子中心利用低能反质子束流开展散裂反应、延迟裂变、反质子原子谱学等与反质子相关的研究，并寻找冷的夸克胶子等离子体物质^[3]。反质子引起的核反应在强子核介质效应、中子皮厚度、奇异超核、不稳定原子核结构等方面有重要应用^[4]。反质子引起的核反应动力学涉及大量反应道，如湮灭反应、电荷交换反应、弹性碰撞与非弹性碰撞等。反质子与核子湮灭主要反应产物是π介子，随着入射能量的升高也会出现超子，如反质子湮灭产生$ \overline{\varLambda}\varLambda $阈值1.439 GeV/c^[5]。因此，高能反质子束流可以产生超子并形成超核。迄今为止，已有一些模型被提出来研究反质子引起的核反应机制，如核内级联模型(Intra-Nuclear Cascade, INC)^[6]、统计碎裂模型(Statistics Multifragmentation Model, SMM)^[7]、微观输运模型等^[8−9]。

本文采用兰州量子分子动力学模型(Lanzhou Quantum Molecular Dynamics transport model, LQMD)研究入射动量为5 GeV/c反质子轰击$ ^{58}{\rm{Ni}} $的粒子产生与碎片分布。采用Skyrme能量密度泛函和相对论协变密度泛函(Relativistic Mean-Field, RMF)分别描述核子的动力学演化，而反质子、介子、超子等演化采用光学势和库仑相互作用势，并对计算结果做了分析。

LQMD模型中考虑了重子-重子、介子-重子和介子-介子碰撞以及质量2 GeV以下共振态的产生和演化，用于质子、反质子和介子引起的散裂反应和重离子碰撞动力学研究，在核物质对称能、强子的介质效应、结团产生、超核形成等方面取得了重要成果。本工作中我们将研究反质子引起的核反应机制，分析在不同平均场演化下反应末态产物的影响。

在LQMD模型中核子的演化可以用泊松括号表示：

根据泊松括号的定义$\{A,\, B\} = \sum\limits_{k}^{{}}{{}}\left(\frac{\partial{A}}{\partial{p_k^{}}}\boldsymbol\cdot\frac{\partial{B}}{\partial{q_k^{}}}-\frac{\partial{A}}{\partial{q_k^{}}}\boldsymbol\cdot\frac{\partial{B}}{\partial{p_k^{}}}\right)$，上式可以表示为

其中：$ p_i^{} $和$ q_i^{} $分别表示反应中第$ i $个核子的动量坐标和空间坐标；H为哈密顿量。在RQMD中，由于一个N体的系统具有4N个空间坐标和4N个动量坐标，但是在相空间中，N体系统只有3N个空间坐标和3N个动量坐标，一共6N个自由度，所以需要降低2个维度将相对论8N减少到相空间6N个自由度，采用质量壳条件和时间约束条件可以将8N个自由度减少至6N个自由度^[10]，此时核子的态可以表示为：$ \phi_i^{} = {p^*_i}^2 - {M^*_i}^2 $。而哈密顿量可以表示为多个态的线性求和^[11]：

其中拉格朗日乘子$ \lambda _i^{}(\tau) $在每个时间步长中由必须始终满足2N个约束的条件确定^[12]。这样上面的核子演化方程可以写成：

在LQMD模型中，单个核子利用高斯波包来表示^[13−14]：

其中：$ {\boldsymbol{r}}_i^{}(t) $和$ {\boldsymbol{p}}_i^{}(t) $为相空间中6维时间相关的参数；L为坐标空间高斯波包宽度的平方。全系统波函数可表示为N体单核子波函数的直积：

由于LQMD模型是半经典框架下的模型，因此需要对核子的量子力学波函数做Wigner变换，得到其在经典力学中的相空间Wigner密度分布函数^[15]：

对上面的相空间分布函数动量和坐标求积分别可得核子的单核子空间密度和动量密度：

其中核子采用固定的波包宽度$ \sigma_r^{} = 2 \; {\mathrm{fm}}$，$ L = \sigma_r^2 $，在RMF下还需要考虑相对论带来的尺缩效应，计算两核子的平均密度时还必须加上洛伦兹因子：

其中：$ M_i^*，M_j^* $是两核子的相对论质量；$ p_i^*，p_j^* $是两核子相对论动量，单位皆为MeV/c；$ {\boldsymbol{r}}_i^{}-{\boldsymbol{r}}_j^{} $是两核子实验室系下的距离。关于核子的哈密顿量具体形式将在后面平均场部分详细介绍。

Skyrme平均场作用力由Skyrme能量密度泛函出发推导得到，本次工作不考虑动量相关项，核子的哈密顿量可以表示为

$ \sqrt{p_i^2 + m_i^2} $是相对论能量，其中$ p_i^{} $和$ m_i^{} $分别表示第i个核子的动量和质量。$ U_{{\mathrm{int}}}^{} $是有效相互作用势，它包括库仑势$ U_{{\mathrm{Coul}}}^{} $与局域势$ U_{{\mathrm{loc}}}^{} $。库仑势可以表示为

其中：$ e_i^{} $和$ e_j^{} $是两核子的电荷数；$ r_{ij}^{} $是两核子距离；erf是误差函数，$ L $是两核子平均高斯波包宽度的平方。而局域相互作用势$ U_{{\mathrm{loc}}}^{} = \int V_{{\mathrm{loc}}}^{}[\rho(r)]{\rm{d}}r $，能量密度泛函$ V_{{\mathrm{loc}}}^{} $又表示为

分别是两体项、三体项、对称能项、表面能与表面对称能项，这里的$ \rho_{\mathrm{n}}^{} $、$ \rho_{\mathrm{p}}^{} $和$ \rho $分别表示中子、质子和核子的总密度，$ t_{{\textit{z}}_i^{}}^{}t_{{\textit{z}}_j^{}}^{} $是核子的同位旋第三分量。参数$\alpha,\; \beta,\;\gamma,\;g_{{\mathrm{sur}}}^{},\;g_{{\mathrm{sur}}}^{{\mathrm{iso}}},\;C_{{\mathrm{sym}}}^{}$和$ \rho_0^{} $的值分别取−226.5, 173.7 MeV, 1.309，23，−2.7 MeV·fm²，38 MeV和0.16 fm^−3[16]。

相对论平均场作用力是采用介子交换实现的也就是介子场标量介子和矢量介子分别给核子之间提供中距吸引与短距排斥^[17]。引入σ介子的非线性自相互作用来合理地降低不可压缩性^[18]，研究过程中我们还考虑了同位旋矢量矢量ρ介子和同位旋矢量标量δ介子^[19−20]，系统拉氏量如下：

这里M_N = 938 MeV是核子在自由空间的质量，g_i(σ，ω，ρ，δ)是四种介子作用的耦合常数，m_i(σ，ω，ρ，δ)是四种介子的质量。$ g_2^{} $和$ g_3^{} $是σ介子非线性自相互作用的耦合常数，$ F_{\mu v}^{} = \partial_\mu^{}\omega_v^{} - \partial_v^{}\omega_\mu^{} $和$ {{\boldsymbol{B}}}_{\mu v}^{} = \partial_\mu^{} {{\boldsymbol{b}}}_v^{} - \partial_v^{} {{\boldsymbol{b}}}_\mu^{} $分别表示矢量介子ω和ρ的全反对称张量。介子与核子的运动方程可以用欧拉-拉格朗日等式得到，利用相对论平均场近似^[21−22]、无海近似以及局域密度近似，此时每种介子场的克莱尔戈登方程左边微分项皆为零，最终核子的运动方程可以表示为

其中的ρ和$ \rho _S^{} $分别是重子密度和重子标量密度，$ \rho _3^{} = \rho _{\mathrm{p}}^{}-\rho _{\mathrm{n}}^{} $，$ \rho _{S3}^{} =\rho _{S{\mathrm{p}}}^{}-\rho _{S{\mathrm{n}}}^{} $。在相对论近似下，系统的能量密度泛函可以写成：

这里的$ p_F^i $是核子费米动量(i = p，n)，$ M_i^{} = M_N^{}- g_\sigma^{}\sigma \pm g_\delta^{}\delta_3^{} $是核子的有效质量(质子符号为−，中子符号为+)。

在反质子诱导的核反应中，介子的产生主要通过反质子-核子湮灭反应得到。以及反质子-核子非弹性散射产生共振态$\Delta(1232)，{\mathrm{N}}^*(1440)，{\mathrm{N}}^*(1535)$衰变也可以得到π介子^[23]。奇异粒子主要由反质子-核子湮灭产生，反质子-核子的湮灭反应使用赝标量和矢量介子的SU(3)对称的统计模型描述的^[24]，同时还考虑了几个介子的最终态混合可能^[9]。此外LQMD模型还考虑了电荷交换反应(Charge-exchange reaction, CEX)、弹性散射(Elsatic scattering, EL)和非弹性散射，包含的反应道如下^[5]：

其中B表示$ {\mathrm{N}}, \; \Delta, \; {\mathrm{N}}^* $等核子与共振态，Y表示包括Λ，Σ等在内的超子，$ {\mathrm{K}} $表示${\mathrm{K}}^0,\;{\mathrm{K}}^+$，$ \overline{{\mathrm{K}}} $表示$\overline{{\mathrm{K}}^0},\;{\mathrm{K}}^-$介子，每个反应道的截面均采用实验数据参数化或外推得到。

对于核子我们按照上面的两种平均场进行演化，但对于碰撞产生的介子、超子以及入射的反质子，我们采用的是相对论平均场或手征有效场论来描述，它的哈密顿量表示为^[25−26]

即由三部分组成：库仑相互作用势、核子能量和光学势。其中库仑相互作用为

$ N $和$ N_{\mathrm{B}}^{} $分别为反应体系介子、超子或者反质子总数和重子的总数，光学势$ V_{{\mathrm{opt}}}^{}({{\boldsymbol{p}}_{\boldsymbol{i}}^{}},\, \rho_i^{}) = \omega({{\boldsymbol{p}}_{\boldsymbol{i}}^{}},\, \rho_i^{}) - \sqrt{{\boldsymbol{p}}^2_i+m^2_i} $，$ \omega({{\boldsymbol{p}}_i^{}},\, \rho_i^{}) $是介子、超子和反质子在核介质中的能量。超子在核介质中的能量由色散关系可以表示为^[27]

这里自能标量$\varSigma^{\varLambda,\, \varSigma}_{\mathrm{S}} = 2\varSigma^{\mathrm{N}}_{\mathrm{S}}/3,\; \varSigma^{\varXi}_{\mathrm{S}} = 1\varSigma^{\mathrm{N}}_{\mathrm{S}}/3$和自能矢量$\varSigma^{\varLambda,\, \varSigma}_{\mathrm{V}} = 2\varSigma^{\mathrm{N}}_{\mathrm{S}}/3,\; \varSigma^{\varXi}_V = 1\varSigma^{\mathrm{N}}_{\mathrm{S}}/3$均可由核子的自能标量$ \varSigma^{\mathrm{N}}_{\mathrm{S}} $及矢量$ \varSigma^{\mathrm{N}}_{\mathrm{S}} $得到，由此得到的$ \varLambda $与核子的光学势在饱和密度处($ \rho_0^{} = 0.16 \; {\boldsymbol{}}{{\rm{fm}}^{-3}}$，动量为0)为−32 MeV。

反核子的能量由G宇称变换计算得到^[28]：

反质子的自能标量$ \varSigma^{\overline{{\mathrm{N}}}}_{\mathrm{S}} = \xi\varSigma^{\mathrm{N}}_{\mathrm{S}} $与自能矢量$ \varSigma^{\overline{{\mathrm{N}}}}_{\mathrm{V}} = -\xi \varSigma^{\mathrm{N}}_{\mathrm{V}} $，其中$ \xi $取0.25，由反质子-核子散射拟合得到^[28]，由此得到反质子与核子的光学势在饱和密度下为−164 MeV。

根据前文反应道的介绍，反质子诱导的核反应中存在着非弹性散射。非弹性散射会使得核子激发形成共振态，共振态继续衰变，也会产生大量的粒子，这些粒子又会引发一连串的级联散射。所以研究这些粒子的产生，有利于理解反应体系所经历的复杂动力学过程，特别是各个微观散射反应道之间的关联。本次工作，反应的碰撞参数b取的是0~r₀ fm(r₀是靶核半径)，囊括了对心碰撞与擦边碰撞，核物质的饱和密度取$ \rho_0^{} = 0.16\; {{\rm{fm}}^{-3}} $。相对论平均场我们采用的参数是set2^[21]，此套参数与前面Skyrme平均场所选参数计算所得的对称能软硬接近，可以更好地比较对称核物质在不同平均场下动力学机制的差异。如图1所示，入射动量为5 GeV/c，$ \overline{{\mathrm{p}}} + ^{58}{\rm{Ni}} $反应系统下，π介子、K介子、Λ和Σ超子数量随时间演化的关系。可以发现在反应初期就有大量的介子产生，在反应发生的短时间内迅速产生。随着反应的进行，π介子的产生逐渐平缓，反质子湮灭给靶核带来区域的能量冗余，使得部分散射后的核子处于共振态，共振态会随时间衰变产生π介子，同时π介子与基态核子发生散射也会产生共振态，导致反质子反应结束后π介子会小幅度增加直至共振态产生与衰变完全平衡，K介子的产生在反应初期的情况与π介子相似，然而K⁺和K⁰介子能量较低时与其他核子基本只发生弹性散射，在反应末期数量几乎不再变化，K⁻和$ \overline{{\mathrm{K}}}^0 $介子的奇异数为1，继续散射会与核子发生电荷交换和奇异数交换从而产生超子，具体反应道如：$\pi {\mathrm{N}} \rightarrow \Delta,\; \overline{{\mathrm{K}}}{\mathrm{N}} \rightarrow\pi + {\mathrm{Y}}$等。同时当反核子动能达到实验给出的动能阈值时将产生一对相反的超子，反应道为$ \overline{{\mathrm{N}}}{\mathrm{N}} \rightarrow \overline{\Lambda} \Lambda $，一对相反的夸克湮灭产生一对相反奇异夸克$ {\mathrm{u}}\overline{{\mathrm{u}}}({\mathrm{d}}\overline{{\mathrm{d}}})\rightarrow {\mathrm{s}} \overline{{\mathrm{s}}} $，反应过程遵循电荷守恒与奇异数守恒。因此反应产生的介子与核子发生散射也会产生超子$\pi (\eta, \, \rho,\, \omega) + {\mathrm{N}} \rightarrow {\mathrm{K}} + {\mathrm{Y}},\;\overline{{\mathrm{K}}} + {\mathrm{N}} \rightarrow \pi + {\mathrm{Y}}$。随着湮灭反应和K⁻和$ \overline{{\mathrm{K}}}^0 $介子的电荷、奇异数交换结束超子不再增加，超子的反应截面在mb量级。同时超子本身会发生衰变，通过弱相互作用力超子会衰变成介子或者其他超子$ \Lambda\rightarrow {\mathrm{N}} + \pi, \; \Sigma\rightarrow {\mathrm{N}} + \pi $，它们的平均寿命约为10⁻¹⁰ s($ \varSigma^0 $通过电弱相互作用衰变，$\Sigma^0\rightarrow\Lambda + \gamma$，平均寿命约为10⁻²⁰ s)，因此超子就有强产生弱衰变的特点，在以fm/c为单位的时间系下相对稳定，如图1中反应产生的K⁺和K⁰介子和超子数目在反应初期迅速增加，最后达到平衡。在不同的平均场下所有粒子反应道截面相同，因此产生粒子的多重性和动力学演化相似。

图  1  入射动量5 GeV/c下$ \overline{{\mathrm{p}}}+ ^{58}{\rm{Ni}} $反应(a) π介子多重数，(b) K介子多重数和(c) 超子多重性的时间演化(在线彩图)

5 GeV/c的反质子打中质量靶核发生的是散裂反应，湮灭反应的存在使得反质子能给靶核内部带来更大的能量沉积。如前文所说，能量沉积会使周围核子激发，同时在反质子入射方向上会发生核子之间级联散射，根据动量守恒，部分核子或团簇从靶核中沿入射方向飞出。在这个过程中，由于超子与核子之间的光学势是排斥的，所以产生的超子会飞离反应区，然后被核子结团或者核碎片捕获吸收形成超核。研究碎片的质量分布有利于评估反质子在核内的能量释放，而超核的质量分布可以更好地理解超子的产生以及超核的形成，图2给出了不同平均场下Λ超核碎片和普通核子碎片的分布，横坐标表示质量数或者电荷数。其中，碎片与结团的识别均采用最小树法(Minimum Spanning Tree, MST)，MST的参数r₀=3.5 fm，p₀=0.2 GeV/c。

图  2  入射动量为5 GeV/c下$ \overline{{\mathrm{p}}}+ ^{58}{\rm{Ni}} $反应产生核碎片和超核在不同平均场下(a) 质量和(b)电荷分布(在线彩图)

在相同的反应系统中，不同的平均场下模型给出的碎片分布趋势是一致的，在低质量和高质量的碎片数量比较多，在中质量碎片区形成低谷。因为在随机且均匀的碰撞参数下，从擦边碰撞到中心碰撞都会发生，其中大质量碎片主要是由擦边碰撞产生的，碰撞时反质子会捕获远核半径附近的核子或与其发生弹性与非弹性散射，形成一个小核与一个近靶质量的大核；中质量碎片主要是由近中心碰撞或中心碰撞产生的，中质量碎片描述中心碰撞时核子的碎裂程度，低质量碎片在擦边碰撞和中心碰撞都会产生，它们由反应后动量较大的核子从结团中脱出形成。Λ超核碎片与普通核碎片相比，在近靶质量上的分布差异是由于擦边碰撞产生的超子会受到核子的排斥势向核外飞出，难以被其他核子捕获形成超核。然而，从图中可以发现RMF相比于Skyrme平均场，中质量碎片分布有着较大的差异，说明反应后不同的平均场给核子带来的效应不完全相同。从高质量碎片来看，两者平均场都能保持稳定，然而RMF下在中质量碎片上数量更多，大结团更容易碎裂成中小结团。

实验上，反质子诱导的核反应为核物质产生奇异粒子提供了良好的环境。介子和超子主要是通过反质子湮灭反应道$ \overline{{\mathrm{p}}}{\mathrm{N}} \rightarrow {\rm{annihilation}}(\pi, \; {\mathrm{K}}, \; \Lambda, \; \Sigma, \; \Xi) $产生，它们的演化过程会受到反应区周围的核子强相互作用与多重散射过程的影响。介子、奇异粒子的快度分布、动能不变谱和集体运动等与高温的核物质状态方程有关。快度是表征核子运动速度的物理量，采用快度描述利于坐标洛伦兹变换下核子之间速度或者动量的关系，图3给出了π介子、K介子、Λ和Σ超子的快度分布，可以发现每种核子的分布峰值都向右偏移了一些，同时随着核子速度或者动量增加，核子数量在迅速减少。峰值的偏移是因为反质子诱导的核反应中反应区的动量受到反质子入射能量影响，在湮灭反应过程遵循动量守恒，所以质心的动量方向与入射方向相同，且会随着入射能量的增加峰值动量也会增加。本次工作分析了所有的事件数据，包括没有产生超核的事件，即反质子与靶核子发生弹性散射、非弹性散射，或者完全湮灭放出了γ光子的情况。Λ，Σ超子的中心快度的产生截面在mb量级与前文的时间演化的截面相符，图3中超子的截面偏离快度峰值后迅速减少，超子质量比核子大，由于一对超子的质量较大，根据能量守恒湮灭产生的超子拥有的动量不会太高，介子质量比较小，快度分布则比较广。此外，不同产物的快度分布在不同的平均场下几乎完全相同，说明产物在RMF与Skyrme平均场下的动力学演化是相当的。

图  3  入射动量5 GeV/c下$ \overline{p}+ ^{58}{\rm{Ni}} $的产物快度分布(在线彩图)

动能谱的情况与快度分布相似，不同的是动能不变谱以动能来描述反应产物在平均场下的运动状态，随着入射能量的增加不变谱的截面分布会变得平缓^[29]，因为入射能的增加使得产物得到的平均动能增加，拥有高动能的产物截面也随之增加。图4给出了π介子和Λ+Σ超子的不变谱，随着动能的增加，产物的不变谱截面都在减小。对于比较轻π介子，不变谱的分布比较平缓，高动能的数目比较多；而质量较大的超子，不变谱呈指数迅速减少。此外与产物快度分布相似，平均场对它们分布影响不明显。

图  4  入射动量5 GeV/c下$ \overline{p}+ ^{58}{\rm{Ni}} $的介子和超子不变能谱(在线彩图)

我们基于LQMD模型研究了在不同平均场下入射动量为5 GeV/c反质子轰击$ ^{58}{\rm{Ni}} $粒子产生相空间分布和碎裂机制。结果显示，在相对论平均场下中等质量碎片产额比Skyrme能量密度泛函高出约三倍，反映出更加明显的涨落效应。而湮灭产生的粒子相空间分布不会受到核子平均场的影响，依赖于反质子与核子湮灭反应截面和粒子-核子散射。反质子湮灭产物主要为π介子，奇异粒子产额比π介子低一个数量级。反质子引起的核反应有利于产生“类靶”区域超核，截面在毫靶量级，超核碎片的截面比核碎片截面低约两个数量级。另一方面，反质子引起的核反应在研究强子介质效应方面也具有一定优势。