Effects of nucleon-nucleon short-range correlation on fragment generation in isotopic nuclear reactions

The saturation character of nuclear matter reveals a balance between long-range attractive forces and short-range repulsive forces within the nuclear interaction. This short-range aspect of nuclear forces engenders correlations among nucleon pairs, especially neutron-proton (n-p) pairs, which give rise to nucleons in the ground state populating the phase space with momenta exceeding the Fermi momentum [1−5]. The high-momentum tail (HMT) phenomenon, resulting from nucleon-nucleon short-range correlations (SRCs), has been substantiated through proton removal experiments using high-energy electrons or protons [6−11]. Recent experiments with high momentum transfer indicate that nuclei in the ground state can form nucleon pairs characterized by high relative momentum and low center-of-mass momentum [12−16]. Such findings reveal that approximately 20% of nucleons in the ground state exist in these paired states, focusing interest on these intricate nucleon correlations [17−21]. Furthermore, the HMT profile typically exhibits a $ C/k^4 $ tail across all nuclei, ranging from deuterons to much heavier nuclei [22−25]. In the nucleon momentum distribution within the HMT, the ratio of n-p SRC to p-p or n-n SRC is approximately 18 in heavier nuclei [26, 27], but significantly lower in lighter nuclei (A<5) [28]. Consequently, the nucleon component in the HMT is markedly isospin dependent, indicating distinct structures in the high-momentum wave function, particularly in lighter nuclei.

Using the isospin- and momentum-dependent transport model, one can study the isospin effect in intermediate-energy heavy-ion collisions (HICs) through some isospin-sensitive observables, such as nucleonic collective flows, the free neutron to proton ratio ($ n/p $), and the $ \pi^{-}/\pi^{+} $ ratio [29−32]. They all have been shown to be more or less sensitive to the HMT [33−35]. And the quantum molecular dynamics transport model can better explain the multi-fragmentation and the production of fragments in intermediate-energy HICs [36−38]. Since the effects of nucleon-nucleon SRC on fragments and the light clusters, e.g., triton and $ ^3 {\rm{He}}$ are studied with central $ ^{112} {\rm{Sn}}$ + $ ^{112} {\rm{Sn}}$ and $ ^{124} {\rm{Sn}}$ + $ ^{124} {\rm{Sn}}$ collisions at $ E_{\text{beam}} = 120 $ MeV/u. Before this work, we first constructed the HMT of the initializing nucleus in phase space. Afterwards, the gas-phase nucleons sensitive to the high-momentum tail is studied in the transport model with the HMT of the initializing nucleus, and compared with experimental data. Simultaneously, as a verification of the constructed HMT, it offers certain evidentiary support for the subsequent study of fragments and the light clusters.

In the QMD-like models, the wave function for each nucleon is represented by Gaussian wave packet [39, 40] as follows

where $ {\bf{r}}_{i}(t) $ and $ {\bf{p}}_{i}(t) $ are the centers of the $ i^{th} $ nucleon in the coordinate and momentum space, respectively. $ \sigma_{r} $ is the width of the Gaussian wave packet depending on the mass number of the nucleus in $ \sigma_{r} = 0.92+0.08A^{1/3} $ fm [41]. After performing a Wigner transformation for Eq. (1), the Wigner density is obtained as,

with

The $ i^{th} $ nucleon density distribution in the coordinate and momentum space are

and

respectively, and the total density distributions are the sum of all nucleons in reaction systems.

In the LQMD model, the temporal evolutions of the baryons (nucleons and resonances) and mesons in the reaction system under the self-consistently generated mean-field are governed by Hamilton’s equations of motion [42−46]. On the basis of the Skyrme interactions, isospin-, density-, and momentum-dependent Hamiltonians were constructed. The Hamiltonian of baryons consists of the relativistic energy, effective interaction potential, and momentum related components. The effective interaction potential is composed of the Coulomb potential and the local interactions.

The local interaction potential is derived from the energy-density functional in the form $ U_{\mathrm{loc}} = \int V_\mathrm{loc}(\rho({\bf{r}})) d {\bf{r}} $. The functional $ \int V_{\mathrm{loc}} $ can be written as:

where $ \rho_{\mathrm{n}} $, $ \rho_{\mathrm{p}} $ and $ \rho = \rho_{\mathrm{n}} + \rho_{\mathrm{p}} $ are the neutron, proton, and total densities, respectively, and $ \delta = (\rho_{\mathrm{n}}-\rho_{\mathrm{p}})/(\rho_{\mathrm{n}}+\rho_{\mathrm{p}}) $ is the isospin asymmetry. The parameters α, β, γ, $ g_{\mathrm{sur}} $, $ g_{\mathrm{sur}}^{\mathrm{iso}} $ and $ \rho_{0} $ are set to -215.7 MeV, 142.4 MeV, 1.322 MeV $ \mathrm{fm}^{2} $, -2.7 MeV $ \mathrm{fm}^{2} $, and 0.16 $ \mathrm{fm}^{-3} $, respectively. $ E_{\mathrm{sym}}^{\mathrm{loc}}(\rho) = 26.25\left(\rho / \rho_0\right)^{0.5} $ is the local part of the symmetry energy, which leads to soft-symmetry energy and the effective mass splitting of $ m_\mathrm{n}^* > m_\mathrm{p}^* $ to obtain the symmetry energy of 31.5 MeV at saturation density. In the LQMD model, a Skyrme-type momentum-dependent potential is used [42, 43].

in which $ C_{\tau, \tau} = C_{\mathrm{mom}}(1+x), $ $ C_{\tau,\tau^{\prime}} = C_{\mathrm{mom}}(1-x)\left(\tau \neq \tau^{\prime}\right) $ and the isospin symbols $ \tau(\tau^{\prime}) $ represent the proton or neutron. The parameters $ C_{\mathrm{mom}} $ and $ \epsilon $ were determined by fitting the real part of the optical potential as a function of the incident energy from the proton-nucleus elastic scattering data, and the obtained values of $ C_{\mathrm{mom}} $ and $ \epsilon $ were 1.76 MeV and 500 $ \mathrm{c}^{2}/\mathrm{GeV}^{2} $, respectively. Thus, the effective mass of the nuclear medium at saturation density is $ m^*/m = 0.75 $. The parameter x is the strength of the isospin splitting, for which a value of -0.65 is adopted in this study, and the mass splitting is $ m_\mathrm{n}^* > m_\mathrm{p}^* $ in the nuclear medium [47].

The initial density distributions of nucleons in the projectile and target come from the density distributions of stable nuclei selected by the self-evolution of nucleons in the LQMD model. Each nucleus is divided into many spherical shells centered around its center of mass. By using the local Thomas-Fermi approximation in each shell of radius r, the local Fermi momenta of neutrons and the protons in each shell are calculated according to

where $ \rho(r) $ is given by Eq. (4). The Fermi momentum obtained based on Eq. (8) can provide nuclear momentum space information by using the Monte Carlo sampling method.

We use a nucleon momentum distribution with a high-momentum tail reaching $ \lambda k_{F_{n,p}} = 2.0 k_{F_{n,p}} $ [48], where $ \lambda = k_{max}/k_{F} $ is the high momentum cutoff parameter and $ k_{F_{n,p}} $ is the neutron or proton Fermi momentum [49, 50]. According to the n-p dominance model [26, 27], 20% of nucleons with equal numbers of neutrons and protons are inside the HMT. The momentum distribution for nucleons with high-momentum tail in initialization reads

and

keeping 20% fraction of total nucleons in the HMT, i.e.,

with the normalization condition

In the above equations, the parameters $ C_{1} $ and $ C_{2} $ in Eq. (10) are determined automatically from the above equations. The calculation of momentum distribution also uses the form of Gaussian wave packet in Eq. (5).

The HMT is strongly isospin dependent and the isotopes $ ^{112} {\rm{Sn}}$ and $ ^{124} {\rm{Sn}}$ have been applied to the study of isospin effects. Therefore, the HMT of $ ^{112} {\rm{Sn}}$ and $ ^{124} {\rm{Sn}}$ have been constructed as shown in Fig. 1. The nucleon momentum distributions with and without HMT are compared. Compared to the ideal gas case, protons in the neutron-rich nucleus $ ^{124} {\rm{Sn}}$ have a higher probability than neutrons of possessing momenta exceeding the nuclear without HMT momentum [51]. This phenomenon arises from the n-p dominance model, which dictates that equal numbers of neutrons and protons must populate the HMT.

Figure 1.  (Color online) Comparison between nucleon momentum distribution n(k) of $^{112}{\rm{Sn}}$ and $^{124}{\rm{Sn}}$ with and without HMT. The normalization condition $\int_{0}^{\lambda k_{F}}n(k)k^{2}dk = 1.$

Since the performance of nucleons is relevant to the isospin-dependent HMT, the observables should be kinetic-energy spectra or momentum spectra with isospin correlation. The central collisions of the isotopic nuclear systems A = $ ^{112} {\rm{Sn}}$ + $ ^{112} {\rm{Sn}}$ and B = $ ^{124} {\rm{Sn}}$ + $ ^{124} {\rm{Sn}}$ are simulated at $ E_{beam} $ = 120 MeV/u and an impact parameter $ b = 0-3 $ fm to investigate the effects of the HMT on the kinetic-energy spectra of neutron-to-proton ratios ($ R_{n/p} $) from the yields of gas-phase fragments (nucleons, hydrogen and helium isotopes) in Fig. 2. These particles, emitted perpendicular to the beam direction with a cut at the polar angle of $ 70^{\circ}<\theta_{\mathrm{c.m.}}<110^{\circ} $ $ (\cos\theta_{\mathrm{c.m.}} = p_{z}/ $ $ \sqrt{p_{x}^{2}+p_{y}^{2}+p_{z}^{2}}) $, are analyzed. The nucleons in the gas phase of the experimental data are obtained by combining the free nucleons with those bound in light isotopes with $ 1 < A < 5 $ [52, 53]. As shown in Fig. 2, it can be observed that the value of the gas-phase $ R_{n/p} $ without the HMT is obviously higher than that with the HMT, both in $ ^{112} {\rm{Sn}}$ + $ ^{112} {\rm{Sn}}$ and $ ^{124} {\rm{Sn}}$ + $ ^{124} {\rm{Sn}}$ collisions. The HMT is induced by SRCs and the n-p pairs are predominant. This predominance implies that in the high-kinetic energy region, neutrons and protons are more prone to appear as correlated pairs, resulting in reduced $ R_{n/p} $. The reduction in the $ n/p $ ratio will tend to approach 1, and the extent of the reduction will become smaller as the isospin asymmetry of the reaction system decreases. Therefore, the HMT effect in neutron-rich reaction systems is more noticeable. It can also be seen that in both $ ^{112} {\rm{Sn}}$ + $ ^{112} {\rm{Sn}}$ and $ ^{124} {\rm{Sn}}$ + $ ^{124} {\rm{Sn}}$ collisions under the influence of soft symmetry energy and the effective mass splitting of $ m_\mathrm{n}^{\ast} > m_\mathrm{p}^{\ast} $, the gas-phase $ R_{n/p} $ with the HMT is closer to the experimental data than that without the HMT. On the other hand, the bump structure around the energy of 20 MeV/u comes from the competition of free nucleons and light fragments to $ R_{n/p} $. The light fragments with Z $ \leq $ 2 are mainly produced at low kinetic energy and have smaller $ n/p $ ratios in comparison to the free nucleons [44].

Figure 2.  (Color online)Effects of the HMT on the kinetic-energy spectra of neutron to proton ratios from the yields of gas-phase fragments (nucleons, hydrogen and helium isotopes) [panels (a) and (b)] in $ ^{112} {\rm{Sn}}$ + $ ^{112} {\rm{Sn}}$ and $ ^{124} {\rm{Sn}}$ + $ ^{124} {\rm{Sn}}$ reactions at the beam energy $ E_{beam} $ = 120 MeV/u in the transverse direction $ 70^{\circ}<\theta_{\mathrm{c.m.}}<110^{\circ} $.

The kinetic energy spectra of the gas-phase double neutron-to-proton ratios ($ DR_{n/p} = R_{n/p}(B)/R_{n/p}(A) $) has been calculated at $ E_{beam} $ = 120 MeV/u in the transverse direction $ 70^{\circ}<\theta_{\mathrm{c.m.}}<110^{\circ} $ with and without the HMT as shown in Fig. 3. However, the $ DR_{n/p} $ whether with or without the HMT is less sensitive to the experimental data since the HMT has the same effect on neutron-rich and -deficient systems. Therefore, the $ DR_{n/p} $ serves as a good experimental observable for other isospin effects, such as symmetry energy and nucleon effective mass splitting, because it can be considered without taking into account the HMT, neutron measurement efficiency, and the influence of Coulomb energy.

Figure 3.  (Color online)Effects of the HMT on the kinetic-energy spectra of the double neutron to proton ratios of gas-phase nucleons in collisions of $^{124}{\rm{Sn}}$ + $^{124}{\rm{Sn}}$ over $^{112}{\rm{Sn}}$ + $^{112}{\rm{Sn}}$ at the beam energy $E_{beam}$ = 120 MeV/u in the transverse direction $70^{\circ}<\theta_{\mathrm{c.m.}}<110^{\circ}$.

Among gas-phase nucleons, light clusters such as triton and $ ^{3} {\rm{He}}$ constitute essential components. Therefore, analyzing the triton-to-$ ^{3} {\rm{He}}$ ratio enables us to further explore the isospin-dependent effects within gas-phase nucleons and establish an understanding of the dynamics of light clusters under the influence of HMT. Shown in Fig. 4 is the triton-to-$ ^{3} {\rm{He}}$ ratio as a function of the total momentum in the peripheral collisions of 120 MeV/u $ ^{112} {\rm{Sn}}$ + $ ^{112} {\rm{Sn}}$ and $ ^{124} {\rm{Sn}}$ + $ ^{124} {\rm{Sn}}$. With respect to HMT, it is seen that the triton-to-$ ^{3} {\rm{He}}$ ratio is sensitive to HMT. Moreover, the value of the triton-to-$ ^{3} {\rm{He}}$ ratio with HMT is much smaller than that without HMT, especially at high momenta. This is due to the fact that the HMT arises from SRCs, where n-p pairs dominate. This dominance suggests that neutrons and protons in the high-momentum region are more likely to form correlated pairs, leading to a decrease in the triton-to-$ ^{3} {\rm{He}}$ ratio. For the peripheral collisions, since more (less) energetic nucleons are from the HMT (nucleon-nucleon collisions), the larger effects of the HMT on the triton-to-$ ^{3} {\rm{He}}$ ratio are expected.

Figure 4.  (Color online)Effects of the HMT on the triton to $^{3}{\rm{He}}$ ratio as a function of total momentum in the Sn+Sn reactions at 120 MeV/u with different reaction systems.

Light clusters, such as triton and $ ^{3} {\rm{He}}$, offer insights into isospin dynamics under the influence of the HMT. Extending this analysis, it turns to heavier fragments whose production involves more complex interactions and structural dynamics. Figure 5 shows multiplicity distributions of fragments produced in HICs with and without HMT in the peripheral collisions of $ ^{112} {\rm{Sn}}$ + $ ^{112} {\rm{Sn}}$ and $ ^{124} {\rm{Sn}}$ + $ ^{124} {\rm{Sn}}$. It can be seen that the yields of fragments decrease with the increasing fragment charge whether with or without the HMT and are not affected by the different isotope reaction systems. However, it is clear that there are more yields of very light fragments and less yields of heavier fragments with the HMT compared to without the HMT. This is because the very light fragments are mostly pre-equilibrium emissions, the momentum correlations among nucleons increase the coalescence probability among nucleons, thus more light fragments are produced. The heavier fragments are mainly deexcited remnants of the reaction, thus, heavier fragments are produced with smaller nucleon momenta (i.e., without the HMT).

Figure 5.  (Color online)Effects of the HMT on the multiplicities of fragments as a function of fragment charge in the Sn+Sn reactions at 120 MeV/u with different reaction systems.

Due to the dominant isospin dependence of short-range correlation pairs, in addition to the multiplicity distribution, the influence of the HMT should also be reflected in average neutron to proton ratio $ \langle N/Z\rangle $ of the fragments. Figure 6 shows the $ \langle N/Z\rangle $ as a function of the fragment charge number $ A_{frag} $ for the $ ^{112} {\rm{Sn}}$ + $ ^{112} {\rm{Sn}}$ reactions (left panel) and $ ^{124} {\rm{Sn}}$ + $ ^{124} {\rm{Sn}}$ reactions (right panel) at the incident beam energy of 120 MeV/u with and without the HMT in semiperipheral collision. By the way, to show more clearly the isospin effects of produced isobars, we analyze the average neutron to proton ratio $ \langle N/Z\rangle $ of produced isobars, which is defined as [54]

Figure 6.  (Color online)Effects of the HMT on the average neutron to proton ratio $\langle N/Z\rangle$ as a function of charge number $A_{frag}$ in the Sn+Sn reactions at 120 MeV/u with different reaction systems.

where $ \sum Y\left({ }^A X_Z\right) $ stands for the yield of the isobar $ {}^AX_Z $ of the nuclear mass number A, the summation is over all different isobars. It is obvious that the $ A_{frag} $ distribution of average neutron to proton ratio $ \langle N/Z\rangle $ is sensitive to the HMT especially for heavier fragments from Fig. 6. Moreover, we can see that $ \langle N/Z\rangle $ with the HMT decreases more sharply than that without the HMT, especially for heavier isobars. This is understandable since the HMT causes de-excitation of remnants easily and usually more neutrons are emitted for neutron-rich system. Thus, the neutron-rich $ ^{124} {\rm{Sn}}$ + $ ^{124} {\rm{Sn}}$ reaction system is more sensitive to the HMT than the $ ^{112} {\rm{Sn}}$ + $ ^{112} {\rm{Sn}}$ reaction system in the $ A_{frag} $ distribution of $ \langle N/Z\rangle $.

In this work, the effects of nucleon-nucleon SRCs lead to the HMT are studied in isotopic nuclear reactions at $ E_{beam} = $ 120 MeV/u using the LQMD transport model. The HMT is incorporated into the initialization of nucleons of the isotopes of $ ^{112} {\rm{Sn}}$ and $ ^{124} {\rm{Sn}}$ in phase space. The analysis of single and double ratios of gas phase neutron and proton spectra shows that HMT affects single ratios but not double ratios in 120 MeV/u $ ^{112} {\rm{Sn}}$ + $ ^{112} {\rm{Sn}}$ and $ ^{124} {\rm{Sn}}$ + $ ^{124} {\rm{Sn}}$ reactions. Therefore, the double ratios can serve as a good experimental observable for other isospin effects. The ratio of triton-to-$ ^3 {\rm{He}}$, like the single ratios of gas-phase neutrons and protons, is influenced by the HMT, with its effects being more pronounced in neutron-rich reaction systems. The single ratios with the HMT are consistent with the experimental data, but the double ratios whether with or without the HMT can fit better to the experimental data. Finally, the isospin-dependent HMT significantly alters the fragment multiplicity distribution and the average neutron-to-proton ratio of produced isobars. Interestingly, the average neutron-to-proton ratio of heavier isobars influenced by HMT is also more pronounced in neutron-rich systems. The production of heavier fragments main originates from the de-excitation of remnants, which is influenced by the HMT. This differs from the production mechanism of light clusters.