Nuclear dynamics and particle production near threshold energies in heavy-ion collisions

NUCLEAR PHYSICS AND INTERDISCIPLINARY RESEARCH

Nuclear dynamics and particle production near threshold energies in heavy-ion collisions

Zhao-Qing Feng

Nuclear Science and Techniques

Vol.29, No.3

Article number 40

Published in print 01 Mar 2018

Available online 20 Feb 2018

DOI：10.1007/s41365-018-0379-z

45101

Recent progress of the quantum molecular dynamics model for describing the dynamics of heavy-ion collisions is viewed, in particular the nuclear fragmentation, isospin physics, particle production and in-medium effect, hadron-induced nuclear reactions, hypernucleus etc. The neck fragmentation in Fermi-energy heavy-ion collisions is investigated for extracting the symmetry energy at subsaturation densities. The isospin effects, in-medium properties, and the behavior of high-density symmetry energy in medium and high energy heavy-ion collisions are thoroughly discussed. The hypernuclide dynamics formed in heavy-ion collisions and in hadron induced reactions is analyzed and addressed in the future experiments at the High-Intensity heavy-ion Accelerator Facility (HIAF).

LQMD transport modelHeavy-ion reactionsHadron induced nuclear reactionsSymmetry energyStrangeness production

1 Introduction

The nuclear dynamics in heavy-ion collisions and in hadron induced reactions at the medium and high energies is complicated and associated with the many-body interaction, collision dynamics, fermionic nature etc. Over the past several decades, two sorts of models based on transport theories have been developed for investigating the hot and dense nuclear matter created in heavy-ion collisions, i.e., the quantum molecular dynamics and Boltzmann Uehling Uhlenbeck (BUU) transport models [1-4]. The liquid-gas phase transition formed in dilute nuclear matter and the quark-gluon plasma (QGP) phase transition in dense matter attracted much attention and have obtained progress both in theories and in experiments [5-11]. The isospin degree of freedom in the nuclear matter plays a significant role on the phase transition.

The nuclear equation of state (EOS) is a basic quantity for describing the nuclear matter with the density and temperature. The EOS in the isospin asymmetric nuclear matter is usually expressed with the energy per nucleon as E(ρ, T, δ)=E(ρ, T, 0)+E_sym(ρ,T)δ²+ O(δ²) in terms of baryon density, ρ, isospin asymmetry, δ=(ρ_n-ρ_p)/(ρ_n+ρ_p), and temperature, T. For treating the properties of a nucleus and cold dense matter, the temperature is usually not taken into account. The high-density symmetry energy E_sym is a key ingredient in understanding the stellar structure, the cooling of protoneutron stars, the nucleosynthesis during supernova explosion of massive stars, and even the binary neutron star merger [12, 13]. The symmetry energy at suprasaturation densities can be constrained from the observables emitted from the high-density regime in nucleus-nucleus collisions. Up to now, the high-density behavior of symmetry energy is not well understood. With operating and constructing the new generation radioactive beam facilities in the world, such as the Cooling Storage Ring (Institute of Modern Physics, Lanzhou, China), Facility for Antiproton and Ion Research (GSI, Darmstadt, Germany), RIKEN (Japan), SPIRAL2 (GANIL, Caen, France), and Facility for Rare Isotope Beams (MSU, East Lansing, USA), the isospin physics is to attract much attention both in theories and in experiments. Besides nucleonic observables, particles produced in heavy-ion collisions would be preferable probes for extracting the information of high-density phase diagram. Kaons as probing the high-density EOS were proposed for the first time [14]. The available experimental data from the KaoS collaboration for K⁺ production favored a soft EOS at high baryon densities associated with transport model calculations [15-18]. Similar structure for Λ production on the EOS was found in Ref. [19]. The ratios of isospin particles produced in heavy-ion collisions such as π^-/π⁺, K⁰/K⁺, ∑^-/∑⁺ etc. [20-26], the yields of η etc [27, 28], and the flow difference of isospin particles [29, 30] have been proposed as sensitive probes for extracting the high-density behavior of the nuclear symmetry energy.

Properties of hadrons in dense nuclear matter is related to the Quantum Chromodynamics (QCD) phase structure, in particular the chiral symmetry restoration, phase transition from hadronic to partonic degrees of freedom, hypernucleus formation, nuclear EOS etc [31-33]. Nucleus-nucleus collisions provide possibilities for exploring the in-medium properties of hadrons, which are related to the interaction potentials and production cross sections in nuclear reactions. The energy spectra of particles are distorted by the surrounding nucleons. It has obtained progress in extracting the in-medium properties of hadrons in dense nuclear matter, in particular for strange particles K, $\bar{K}$ , Λ, and ∑ [34-37]. Studies of hypernuclei have attracted much attention over the past several decades. The interested topics related to hypernuclei are the hyperon-nucleon and hyperon-hyperon interactions, opening a new horizon with strangeness (three-dimensional nuclear chart) in nuclear physics and probing the in-medium properties of hadrons and the inner structure of a nucleus [31, 38]. Moreover, hyperons as essential ingredients impact the EOS at high-baryon densities. The strangeness constitution in dense matter softens the high-density EOS, and consequently decreases the mass of neutron stars[39-41]. Since the first observation of Λ-hypernuclide in nuclear multifragmentation reactions induced by cosmic rays in the 1950s [42], a remarkable progress has been obtained in producing hypernuclides via different reaction mechanisms, such as hadron (pion, kaon, proton) induced reactions, bombarding a fixed target with high-energy photons or electrons, and fragmentation reactions with high energy heavy-ion collisions. Experimental collaborations of nuclear physics facilities in the world, e.g., PANDA [43], FOPI/CBM and Super-FRS/NUSTAR [44] at FAIR (GSI, Germany), STAR at RHIC (BNL,USA) [45], ALICE at LHC (CERN) [46], NICA (Dubna, Russia) [47], J-PARC (Japan) [48], and HIAF (IMP, China) [49] have started or planned to investigate hypernuclei and their properties. In the laboratories, strangeness nuclear physics is to be concentrated on the isospin degree of freedom (neutron-rich/proton-rich hypernuclei), multiple or anti-strangeness (s=-2 or s=1), high-density hypermatter production etc. A more localized energy deposition enables the secondary collisions available for producing hyperons in antiproton-nucleus collisions. The cold quark-gluon plasma (QGP) could be formed from the annihilation of antiprotons on nuclei [50, 51], and the searching is still undergoing in experiments.

In this article, I will review the recent progress on the theoretical treatment of nuclear dynamics in heavy-ion collisions. The symmetry energy, strange particles and in-medium effects, and hyperfragment formation are to be discussed. In Sect.2, I will give a transport model description for simulating the nuclear dynamics essentially based on the quantum molecular dynamics (QMD). The nuclear dynamics and isospin effect in heavy-ion collisions are shown in Sect. 3. The hadron induced reactions and hypernuclide production are investigated in Sect. 4. Summary and perspective on the medium-energy nuclear dynamics are presented in Sect. 5.

2 The quantum molecular dynamics model and its modification

In the QMD transport model, the ith nucleon is represented with a Gaussian wave packet [3, 4]

ψ_{i} (r, t) = \frac{1}{{(2 π σ_{r}^{2})}^{3 / 4}} \exp [- \frac{{(r - r_{i} (t))}^{2}}{4 σ_{r}^{2}}] \cdot \exp (\frac{i p_{i} (t) \cdot r}{ℏ}) .

(1)

Here ri(t), and pi(t) are the centers of the wave packet in the coordinate and momentum space, respectively. The σr is the width of the Gaussian wave packet. The N-body wave function is assumed as the direct product of the coherent states $Φ (r, t) = \prod_{i} ψ_{i} (r, t)$ , and the anti-symmetrization is neglected. Therefore, the fermionic nature in the dynamical evolution is lost. After performing Wigner transformation, the phase-space density is obtained as

f (r, p, t) = \sum_{i} f_{i} (r, p, t),

(2)

\begin{array}{l} f_{i} (r, p, t) = \frac{8}{h^{3}} \exp [- \frac{{(r - r_{i} (t))}^{2}}{2 σ_{r}^{2}}] \\ \times \exp [- \frac{2 (p - p_{i} (t {))}^{2} \cdot σ_{r}^{2}}{ℏ^{2}}] . \end{array}

(3)

The density distributions $ρ (r, t) = \sum_{i} ρ_{i} (r, t)$ in the coordinate and $g (p, t) = \sum_{_{i}} g_{i} (p, t)$ in the momentum space can be easily evaluated with the phase-space density as

ρ_{i} (r, t) = \frac{1}{{(2 π)}^{3 / 2} σ_{r}^{3}} \exp [- \frac{{(r - r_{i} (t))}^{2}}{2 σ_{r}^{2}}],

(4)

g_{i} (p, t) = (8 π)^{3 / 2} {(σ_{r} / h)}^{3} \exp [- \frac{2 (p - p_{i} (t {))}^{2} \cdot σ_{r}^{2}}{ℏ^{2}}] .

(5)

The QMD transport model was improved by different groups, such as implementing the Pauli potential into the mean-field potential [52], extending to the relativistic heavy-ion collisions (TuQMD, UrQMD) [36, 53], introducing the cooling method in the initialization and the time-dependent wave-packet width (EQMD) [54, 55], distinguishing the isospin degree of freedom (IQMD) [56, 57], taking the phase-space constraint approach into the dynamical evolution and nucleon-nucleon collisions (CoMD) [58]. Recently, the ImQMD model by the CIAE group was proposed for describing the fusion-fission, quasifission dynamics, multinucleon transfer, symmetry energy etc [59-63]. The mean-field potential in the ImQMD model is derived selfconsistently from the Skyrme energy-density function. The shell effect in the fusion reactions was investigated within the framework of the QMD model by the Lanzhou group and a phenomenological formula was proposed [64, 65]. Recently, the giant dipole resonance formed in heavy-ion collisions was investigated by SINAP group for probing the α clustering structure in a nucleus and the EOS within the QMD model [66, 67].

In this work, I will give a description of the quantum molecular dynamics model developed by the Lanzhou group (LQMD). The dynamics of the resonances (Δ(1232), N*(1440), N*(1535), etc), hyperons (Λ, ∑, Ξ, Ω), and mesons (π, η, K, $\bar{K}$ , ρ, ω) are described via hadron-hadron collisions, Pauli-blocking, decays of resonances, mean-field potentials, and corrections on threshold energies of elementary cross sections [68, 69]. Besides the hadron-hadron collisions, we have further included the annihilation channels, charge-exchange reaction, elastic and inelastic collisions in antinucleon-nucleon collisions for understanding antiproton induced reactions [70]. The LQMD model has been applied to investigate the medium and high energy heavy-ion collisions, antiproton induced reactions, meson-nucleus dynamics etc.

In the LQMD model, the time evolutions of the baryons (nucleons and resonances) and mesons in the reaction system are governed by Hamilton’s equations of motion under the self-consistently generated mean-field potential as

{\dot{p}}_{i} = - \frac{\partial H}{\partial r_{i}}, {\dot{r}}_{i} = \frac{\partial H}{\partial p_{i}} .

(6)

The Hamiltonian of baryons consists of the relativistic energy, the effective potential, and the momentum dependent interaction. The effective potential is composed of the Coulomb interaction and the local potential

U_{int} = U_{Coul} + U_{loc} .

(7)

The Coulomb interaction potential is written as

U_{Coul} = \frac{1}{2} \sum_{i, j, j \neq i} \frac{e_{i} e_{j}}{r_{i j}} e r f (r_{i j} / 2 σ_{r}),

(8)

where the ej is the charged number including protons and charged resonances. The rij=|ri-rj| is the relative distance of two charged particles.

The local potential can be derived from the energy-density functional as

U_{loc} = \int V_{loc} (ρ (r)) d r,

(9)

which is constructed with the 2-body and 3-body Skyrme-type interaction as

U_{loc}^{2} = \frac{a_{2}}{2} \sum_{i, j, j \neq i} ρ_{i} (r) δ (r - r^{'}) ρ_{j} (r^{'}) d r d r^{'},

(10)

and

\begin{array}{l} U_{loc}^{3} = \frac{a_{3}}{6} \sum_{i, j, k, j \neq i, j \neq k, i \neq k} ρ_{i} (r) ρ_{j} (r^{'}) ρ_{k} (r^{''}) δ (r - r^{'}) \\ \times δ (r - r^{''}) d r d r^{'} d r^{''} . \end{array}

(11)

In the model, the 2-body nucleon-nucleon (NN) interaction includes the attractive potential, symmetry energy, surface term and momentum dependent interaction. The 3-body interaction is embodied in the short-range repulsive force and an approximation treatment is taken in evaluating Eq. (11) [4]. The energy-density functional is expressed by

\begin{array}{l} V_{loc} (ρ) = \frac{α}{2} \frac{ρ^{2}}{ρ_{0}} + \frac{β}{1 + γ} \frac{ρ^{1 + γ}}{ρ_{0}^{γ}} + E_{sym}^{loc} (ρ) ρ δ^{2} \\ + \frac{g_{sur}}{2 ρ_{0}} {(\nabla ρ)}^{2} + \frac{g_{sur}^{iso}}{2 ρ_{0}} {[\nabla (ρ_{n} - ρ_{p})]}^{2}, \end{array}

(12)

where the ρ_n, ρ_p, and ρ=ρ_n+ρ_p are the neutron, proton, and total densities, respectively, and the δ=ρ_n-ρ_p)/(ρ_n+ρ_p) is the isospin asymmetry. The surface coefficients g_sur and $g_{sur}^{iso}$ are taken as 23 MeV fm² and -2.7 MeV fm², respectively. The bulk parameters α, β, and γ are readjusted after inclusion of the momentum term in order to reproduce the compression modulus of symmetric nuclear matter and the binding energy of isospin symmetric nuclear matter at saturation density. The $E_{sym}^{loc}$ is the local part of the symmetry energy, which can be adjusted to mimic predictions of the symmetry energy calculated by microscopical or phenomenological many-body theories and has two-type forms as follows:

E_{sym}^{loc} (ρ) = \frac{1}{2} C_{sym} {(ρ / ρ_{0})}^{γ_{s}},

(13)

and

E_{sym}^{loc} (ρ) = a_{sym} (ρ / ρ_{0}) + b_{sym} {(ρ / ρ_{0})}^{2} .

(14)

The parameters C_sym, a_sym, and b_sym are taken as the values of 52.5 MeV, 43 MeV, -16.75 MeV and 23.52 MeV, 32.41 MeV, -20.65 MeV corresponding to the mass splittings of $m_{n}^{*} > m_{p}^{*}$ and $m_{n}^{*} < m_{p}^{*}$ , respectively. The parameter C_sym is taken as the value of 38 MeV for the case without the mass splitting. The values of γ_s=0.5, 1., 2. have the soft, linear and hard symmetry energy, respectively, and the last formula corresponding to a supersoft symmetry energy.

A Skyrme-type momentum-dependent potential is used in the LQMD model [71]

\begin{array}{l} U_{mom} = \frac{1}{2 ρ_{0}} \sum_{i, j, j \neq i} \sum_{τ, τ^{'}} C_{τ, τ^{'}} δ_{τ, τ_{i}} δ_{τ^{'}, τ_{j}} \int \int \int d p d p' d r \\ \times f_{i} (r, p, t) {[\ln (\in {(p - p^{'})}^{2} + 1]}^{2} f_{j} (r, p', t) . \end{array}

(15)

Here, Cτ,τ=C_mom(1+x), Cτ,τ’=C_mom(1-x) (τ≠τ’) and the isospin symbols τ(τ’) represent a proton or neutron. The parameters C_mom and ϵ were determined by fitting the real part of the optical potential as a function of incident energy from the proton-nucleus elastic scattering data. The parameter x is the strength of the isospin splitting of nucleon effective mass, e.g., the value of -0.65 leading to the splitting of $m_{n}^{*} > m_{p}^{*}$ in nuclear medium. Both sets of the parameters with momentum dependent and independent interactions give the same compression modulus of K=230 MeV for isospin symmetric nuclear matter, as shown in Table 1.

The parameters and properties of isospin symmetric EoS used in the LQMD model at the density of 0.16 fm^-3.

Parameters	α (MeV)	β (MeV)	γ	C_mom (MeV)	ϵ (c²/MeV²)	$m_{\infty}^{*} / m$	K_∞ (MeV)
PAR1	-215.7	142.4	1.322	1.76	5×10^-4	0.75	230
PAR2	-226.5	173.7	1.309	0.	0.	1.	230

The symmetry energy per nucleon in the LQMD model is composed of three parts, namely the kinetic energy, the local part, and the momentum dependence of the potential energy as

E_{sym} (ρ) = \frac{1}{3} \frac{ℏ^{2}}{2 m} {(\frac{3}{2} π^{2} ρ)}^{2 / 3} + E_{sym}^{loc} (ρ) + E_{sym}^{mom} (ρ) .

(16)

After an expansion to second order around the normal density, the symmetry energy can be expressed as

E_{sym} (ρ) \approx E_{sym} (ρ_{0}) + \frac{L}{3} (\frac{ρ - ρ_{0}}{ρ_{0}}) + \frac{K_{sym}}{18} {(\frac{ρ - ρ_{0}}{ρ_{0}})}^{2}

(17)

in terms of a slope parameter of $L = 3 ρ_{0} (\partial E_{sym} / \partial ρ {) |}_{ρ = ρ_{0}}$ and a curvature parameter of $K_{sym} = 9 ρ_{0}^{2} (\partial^{2} E_{sym} / \partial ρ^{2} {) |}_{ρ = ρ_{0}}$ . The values of slope parameters are 203.7 MeV, 124.9 MeV, 85.6 MeV, and 74.7 MeV for the hard, linear, soft, and supersoft symmetry energies, respectively, and the corresponding 448 MeV, -24.5 MeV, -83.5 MeV, and -326 MeV for the curvature parameters. Figure 1 is a comparison of the isospin-symmetric nuclear EOS and nuclear symmetry energy. The EOS of different compression modulus and the symmetry energy of different stiffness cross at the saturation density with the values of -16 MeV and 31.5 MeV, respectively. The difference is pronounced at the high density. The observables from the reaction zone in heavy-ion collisions can extract the high-density EOS.

Figure 1:

(Color online) Density dependence of the symmetric nuclear equation of state at different compression modulus and symmetry energies of supersoft, soft, linear, and hard cases.

The single-particle optical potential can be obtained from the density functional as [72]

\begin{array}{l} U_{τ} (ρ, δ, p) = \frac{1}{ρ_{0}} C_{τ, τ} \int d p' f_{τ} (r, p) {[\ln (ε {(p - p^{'})}^{2} + 1)]}^{2} \\ + \frac{1}{ρ_{0}} C_{τ, τ^{'}} \int d p' f_{τ^{'}} (r, p) {[\ln (\in {(p - p^{'})}^{2} + 1)]}^{2} \\ + α \frac{ρ}{ρ_{0}} + β \frac{ρ^{γ}}{ρ_{0}^{γ}} + E_{sym}^{loc} (ρ) δ^{2} + \frac{\partial E_{sym}^{loc} (ρ)}{\partial ρ} ρ δ^{2} \\ + E_{sym}^{l o c} (ρ) ρ \frac{\partial δ^{2}}{\partial ρ_{τ}} . \end{array}

(18)

Here τ≠τ’, ∂δ²/∂ρn=4∂δ_ρ/ρ² and ∂δ²/∂ρ_p=-4δρ_n/ρ₂. The nucleon effective (Landau) mass in nuclear matter of isospin asymmetry δ=ρ_n-ρ_p)/(ρ_n+ρ_p) with ρ_n and ρ_p being the neutron and proton density, respectively, is calculated through the potential as $m_{τ}^{*} = m_{τ} / (1 + \frac{m_{τ}}{| p |} | \frac{d U_{τ}}{d p} |)$ with the free mass mτ at Fermi momentum p=p_F. Figure 2 is the single-particle potentials with the mass splittings of $m_{n}^{*} > m_{p}^{*}$ and $m_{n}^{*} > m_{p}^{*}$ . It should be noticed that the neutron and proton optical potentials are across with increasing the momentum for the mass splitting of $m_{n}^{*} > m_{p}^{*}$ . The case of $m_{n}^{*} > m_{p}^{*}$ enlarges the difference in the neutron-rich matter. The variance of optical potential influences the dynamical evolution of nucleons and the kinetic energy spectra of the neutron to proton ratios are changed [30]. It was found that the neutron/proton ratio at high transverse momenta (kinetic energies) and elliptic flow difference between neutron and proton sensitively are sensitive probes for constraining the isospin splitting. The isospin splitting of nucleon effective mass in neutron-rich nuclear matter has been planning in the RIKEN-SAMURAI experiments.

Figure 2:

(Color online) Momentum dependence of single-nucleon optical potential at the saturation density for isospin symmetric matter with the solid lines (δ=0) and neutron-rich matter (δ=0.2) with the mass splittings of

m_{n}^{*} > m_{p}^{*}

(left panel) and

m_{n}^{*} < m_{p}^{*}

(right panel), respectively.

The dynamics of mesons, hyperons, and antiprotons is influenced by the mean-field potentials. The Hamiltonian is constructed as follows

H_{ν} = \sum_{i = 1}^{N_{ν}} (V_{i}^{C o u l} + ω_{ν} (p_{i}, ρ_{i})) .

(19)

Here, ν denotes the hadron species, and the Coulomb potential is evaluated by point-like charged particles as

V_{i}^{C o u l} = \sum_{j = 1}^{N_{h}} \frac{e_{i} e_{j}}{r_{i j}} .

(20)

The Nν and Nh are the total numbers of the same species and all charged particles. The in-medium energy ων(pi,ρi) is related to the hadron momentum and baryon density, which is usually evaluated from the scalar and vector potential components to hadron self-energies.

The energy of pion in the nuclear medium is composed of the isoscalar and isovector contributions as follows

ω_{π} (p_{i}, ρ_{i}) = ω_{isoscalar} (p_{i}, ρ_{i}) + C_{π} τ_{z} δ {(ρ / ρ_{0})}^{γ_{π}} .

(21)

The coefficient $C_{π} = ρ_{0} ℏ^{3} / (4 f_{π}^{2}) = 36$ MeV, and the isospin quantities are taken as τz= 1, 0, and -1 for π^-, π⁰, and π⁺, respectively. The isospin asymmetry δ=(ρ_n-ρ_p)/(ρ_n+ρ_p) and the quantity, γπ, adjusts the isospin splitting of the pion optical potential. Usually, we take the γπ = 2 having the obvious difference of π^--nucleon and π⁺-nucleon potentials in the dense neutron-rich matter. The isoscalar pion energy in the nuclear medium is calculated with the Δ-hole model [69, 73, 74]. The in-medium dispersion relation consists of a pion branch (smaller value) and a Δ-hole (larger value) branch. The π-like energy in nuclear medium decreases with increasing the baryon density. However, the Δ-like energy is just the opposite. The dispersion relation reads

\begin{array}{l} ω_{isoscalar} (p_{i}, ρ_{i}) = S_{π} (p_{i}, ρ_{i}) ω_{π - like} (p_{i}, ρ_{i}) + \\ S_{Δ} (p_{i}, ρ_{i}) ω_{Δ - like} (p_{i}, ρ_{i}) . \end{array}

(22)

The probabilities of the pion component satisfy the relation

S_{π} (p_{i}, ρ_{i}) + S_{Δ} (p_{i}, ρ_{i}) = 1

(23)

The value of the probability is determined from the pion self-energy as [75]

S (p_{i}, ρ_{i}) = \frac{1}{1 - \partial Π (ω) / \partial ω^{2}},

(24)

where the pion self-energy is given by

Π = p_{i}^{2} \frac{χ}{1 - g' χ},

(25)

with the Migdal parameter being g′∼0.6 and

χ = - \frac{8}{9} {(\frac{f_{Δ}}{m_{π}})}^{2} \frac{ω_{Δ} ρ ℏ^{3}}{ω_{Δ}^{2} - ω^{2}} \exp (- 2 p_{i}^{2} / b^{2}) .

(26)

In $ω_{Δ} = \sqrt{m_{Δ}^{2} + p_{i}^{2}} - m_{N}$ , the mπ, m_N and m_Δ are the masses of pion, nucleon, and delta, respectively. The πNΔ has a coupling constant of f_Δ 2 and a cutoff factor of b 7 mπ. Two eigenvalues of ω_π-like and ω_Δ-like are obtained from the pion dispersion relation as

ω^{2} = p_{i}^{2} + m_{π}^{2} + Π (ω) .

(27)

The Hamiltonian of η is composed of

H_{η} = \sum_{i = 1}^{N_{η}} (\sqrt{m_{η}^{2} + p_{i}^{2}} + V_{η}^{opt} (p_{i}, ρ_{i})) .

(28)

The eta optical potential is evaluated from the dispersion relation based on the chiral perturbation theory [76] as

ω_{η} (p_{i}, ρ_{i}) = \sqrt{(m_{η}^{2} - a_{η} ρ_{s}) {(1 + b_{η} ρ_{s})}^{- 1} + p_{i}^{2}}

(29)

with $a_{η} = ℏ^{3} \frac{Σ_{η N}}{f_{π}^{2}}$ and $b_{η} = ℏ^{2} \frac{κ}{f_{π}^{2}}$ . The pion decay constant was fπ=92.4 MeV, ∑_ηN =280 MeV, and κ =0.4 fm. The optical potential is given by $V_{η}^{opt} (p_{i}, ρ_{i}) = ω_{η} (p_{i}, ρ_{i}) - \sqrt{m_{η}^{2} + p_{i}^{2}}$ with the eta mass mη=547 MeV. The value of $V_{η}^{opt} =$ -94 MeV is obtained with zero momentum and a saturation density ρ=ρ0.

The energies of kaon and anti-kaon in nuclear medium are calculated with the chiral Lagrangian approach as follows [69, 77]

\begin{array}{l} ω_{K} (p_{i}, ρ_{i}) = [m_{K}^{2} + p_{i}^{2} - a_{K} ρ_{i}^{S} - τ_{3} c_{K} ρ_{i 3}^{S} \\ + {(b_{K} ρ_{i} + τ_{3} d_{K} ρ_{i}_{3})}^{2}]^{1 / 2} + b_{K} ρ_{i} \\ + τ_{3} d_{K} ρ_{i 3} \end{array}

(30)

and

\begin{array}{l} ω_{\bar{K}} (p_{i}, ρ_{i}) = [m_{\bar{K}}^{2} + p_{i}^{2} - a_{\bar{K}} ρ_{i}^{S} - τ_{3} c_{K} ρ_{i 3}^{S} \\ + {(b_{K} ρ_{i} + τ_{3} d_{K} ρ_{i 3})}^{2}]^{1 / 2} - b_{K} ρ_{i} \\ - τ_{3} d_{K} ρ_{i 3}, \end{array}

(31)

respectively. Here the $b_{K} = 3 ℏ^{3} / (8 f_{π}^{2}) \approx 0.337$ GeVfm³, the aK and $a_{\bar{K}}$ are 0.18 GeV²fm³ and 0.31 GeV²fm³. The τ₃=1 and -1 for the isospin pair K⁺( ${\bar{K}}^{0}$ ) and K⁰(K^-), respectively. The parameters cK=0.0298 GeV²fm³ and dK=0.111 GeVfm³ determine the isospin splitting of kaons in neutron-rich nuclear matter. The optical potential of kaon is derived from the in-medium energy as $V_{opt} (p, ρ) = ω (p, ρ) - \sqrt{p^{2} + m_{K}^{2}}$ . The optical potentials being the values of 31 MeV, 25 MeV, -96 MeV, and 104 MeV at saturation density with the isospin asymmetry of 0.2 are obtained for K⁺, K⁰, K^-, and ${\bar{K}}^{0}$ in nuclear medium, respectively. The values of $m_{K}^{*} / m_{K}$ =1.056 and $m_{\bar{K}}^{*} / m_{\bar{K}}$ = 0.797 at normal baryon density are concluded with the parameters in isospin symmetric nuclear matter. The effective mass is used to evaluate the threshold energy for kaon and antikaon production, e.g., the threshold energy in the pion-baryon collisions $\sqrt{s_{th}} = m_{Y}^{*} + m_{K}^{*}$ .

The in-medium dispersion relation for hyperons reads as

ω_{Y} (p_{i}, ρ_{i}) = \sqrt{{(m_{Y} + Σ_{S}^{Y})}^{2} + p_{i}^{2}} + Σ_{V}^{Y} .

(32)

The hyperon self-energies are evaluated on the basis of the light-quark counting rules, i.e., Λ and ∑ being assumed to be two thirds of nucleon self-energies, and the Ξ self-energy being one third of nucleon’s ones. Namely, for hyperons $Σ_{S}^{Λ, Σ} = 2 Σ_{S}^{N} / 3$ , $Σ_{V}^{Λ, Σ} = 2 Σ_{V}^{N} / 3$ , $Σ_{S}^{Ξ} = Σ_{S}^{N} / 3$ , and $Σ_{V}^{Ξ} = Σ_{V}^{N} / 3$ . The antibaryon energy is computed from the G-parity transformation of baryon potential as

ω_{\bar{B}} (p_{i}, ρ_{i}) = \sqrt{{(m_{\bar{B}} + Σ_{S}^{\bar{B}})}^{2} + p_{i}^{2}} + Σ_{V}^{\bar{B}}

(33)

with $Σ_{S}^{\bar{B}} = Σ_{S}^{B}$ and $Σ_{V}^{\bar{B}} = Σ_{V}^{\bar{B}}$ . The nuclear scalar $Σ_{S}^{N}$ and vector $Σ_{V}^{N}$ self-energies are computed from the well-known relativistic mean-field model with the NL3 parameter ( $g_{σ N}^{2}$ =80.8, $g_{ω N}^{2}$ =155 and $g_{ρ N}^{2}$ =20). The optical potential of baryon or antibaryon is derived from the in-medium energy as $V_{o p t} (p, ρ) = ω (p, ρ) - \sqrt{p^{2} + m^{2}}$ . The relativistic self-energies are used for the construction of hyperon and antibaryon potentials only. A very deep antiproton-nucleus potential is obtained with the G-parity approach being $V_{o p t} (p = 0, ρ = ρ_{0}) = - 655$ MeV. From fitting the antiproton-nucleus scattering [78] and the real part of phenomenological antinucleon-nucleon optical potential [79], a factor, ξ, is introduced in order to moderately evaluate the optical potential as $Σ_{S}^{\bar{N}} = ξ Σ_{S}^{N}$ and $Σ_{V}^{\bar{N}} = - ξ Σ_{V}^{N}$ with ξ=0.25, which leads to the strength of $V_{\bar{N}} = - 164$ MeV at the normal nuclear density ρ₀=0.16 fm^-3. It should be noted that the scaling approach violates the fundamental G-symmetry. The antihyperon potentials exhibit strong attractive interaction in the nuclear medium, e.g., the strengths at saturation density being -436 MeV and -218 MeV for $\bar{Λ}$ and $\bar{Ξ}$ , respectively. The optical potentials will affect the dynamics of hyperons, consequently the production of hypernuclei. Figure 3 is the optical potentials for pions, etas, kaons, and hyperons. The empirical value extracted from hypernucleus experiments [31, 80] can be well reproduced by the optical potential. The same approach is used for the ∑ propagation in nuclear medium. The optical potentials are of importance in binding a hyperon into a nucleus.

Figure 3:

Density dependence of the optical potentials for pions, etas, kaons, and hyperons in dense nuclear matter. The empirical value of Λ potential extracted from hypernucleus experiments is denoted by the solid circle [80].

The probability of two-particle or three-particle collision is sampled with a Monte Carlo procedure, in which the scattering of two particles is determined by a geometrical minimum distance criterion, $d \leq \sqrt{0.1 σ_{tot} / π}$ fm, weighted by the Pauli blocking of the final states [2, 81]. Here, the total NN cross section, σ_tot, in mb is the sum of the elastic and all inelastic baryon-baryon collisions. The channel probability is evaluated by comparing the channel to the total cross section as P_ch=σ_ch/σ_tot. The choice of reaction channel is randomly performed by the weight of the channel probability. The total, elastic and inelastic NN cross sections are parameterized in accordance with the available experimental data [82], as shown in Fig. 4. One expects that the in-medium NN elastic collisions are reduced in comparison with the free-space ones. The in-medium elastic cross section is scaled according to the effective mass that was used in the BUU model [83]. The in-medium elastic cross section is scaled according to the effective mass through $σ_{NN}^{medium} = (μ_{NN}^{*} / μ_{NN})^{2} σ_{NN}^{free}$ with the $μ_{NN}^{*}$ and $μ_{NN}$ being the reduced masses of colliding nucleon pairs in the medium and in the free space [84], respectively. Figure 5 is a comparison of the scaling factor as functions of baryon density, nucleon momentum, and isospin asymmetry for the different mass splittings of $m_{n}^{*} > m_{p}^{*}$ and $m_{n}^{*} < m_{p}^{*}$ , respectively. It is interesting to notice that a splitting of nn (neutron-neutron), pn (proton-neutron), and pp (proton-proton) exists because of the mixing of neutron and proton effective masses. The mass splitting results in an opposite splitting of the in-medium cross sections of nn and pp channels. The NN elastic cross section impacts the collective flows in heavy-ion collisions, in particular the balance energy [72]. It is found that the in-medium elastic cross sections play a significant role in isospin emission and result in a flatter distribution for transverse flows and elliptic flows of free nucleons compared with the in-vacuum ones. However, the flow difference of neutrons and protons is independent on the in-medium modifications of NN cross sections and more sensitive to the stiffness of symmetry energy.

Figure 4:

Nucleon-nucleon total, elastic, and inelastic cross sections in free space parameterized in the LQMD model as a function of nucleon incident energy in the laboratory frame.

Figure 5:

In-medium nucleon-nucleon elastic cross section scaled by nucleon effective mass as functions of baryon density (δ=0.2), nucleon momentum (δ=0.2), and isospin asymmetry (ρ=ρ₀).

The inelastic NN collisions are complicated and a number of reaction channels are opened and coupled to each other. The primary products in NN collisions are the resonances and their decay processes, such as Δ(1232), N^*(1440), N^*(1535) etc. We have included the reaction channels as follows:

\begin{array}{l} N N \leftrightarrow N Δ, N N \leftrightarrow N N^{*}, N N \leftrightarrow Δ Δ, \\ Δ \leftrightarrow N π, N^{*} \leftrightarrow N π, N N \leftrightarrow N N π (s - s t a t e), \\ N^{*} (1535) \leftrightarrow N η . \end{array}

(34)

Here, hadron-hadron collisions take place as two-body processes and three-body (s-state pion production) reactions. The momentum-dependent decay widths are used for the resonances of Δ(1232) and N^*(1440) [85, 86]. We have taken a constant width of Γ=150 MeV for the N^*(1535) decay. Elastic scattering of nucleon-resonance (NR→ NR) and resonance-resonance (RR→ RR) collisions and inelastic collisions of nucleon-resonance (NR→ NN, NR→ NR′) and resonance-resonance (RR→ NN, RR→ RR′, R and R′ being different resonances) have been included in the model. A parameterized cross section is used in the LQMD model for the channel of NΔ→NN by fitting the available experimental data [87, 88] similar to the GiBUU model [89]. Calculations based on relativistic Dirac-Brueckner also favor the decrease trend of the Δ absorption cross section with increasing energy in nuclear medium at low densities (less than 2ρ₀) [90]. We used the parametrized cross sections calculated by the one-boson exchange model [86] for resonance production and the absorption of N^* with the detailed balancing principle. Figure 6 is the energy dependence of resonance production in proton-proton collisions. It is obvious that the Δ(1232) production is dominant in the medium-energy heavy-ion collisions. The mass of resonance is usually sampled with the Breit-Wigner (BW) distribution as

Figure 6:

Production and inverse absorption cross sections of the resonances Δ(1232), N^*(1440) and N^*(1535). The experimental data of the channel NΔ→NN are shown for comparison [87, 88].

p (m_{r}) = \frac{1}{1 + 4 [(m_{r} - m_{r}^{0}) / Γ]^{2}}

(35)

with the peak mass, $m_{r}^{0}$ , and the decay width, Γ. Figure 7 is the mass distributions of the resonances Δ(1232) and N^*(1440) in the ¹⁹⁷Au+¹⁹⁷Au collisions. A modified formula proposed by Pawel [91] is also shown for the Δ(1232) mass sampling at the incident energy of 0.4 GeV/nucleon. The distribution is slightly moved towards the lower mass region with the Pawel’s formula in comparison with the standard BW form.

Figure 7:

Mass distributions of the resonances Δ(1232) and N^*(1440) in the ¹⁹⁷Au+¹⁹⁷Au reaction at the incident energies of 400, 1000, and 1500 A MeV, respectively.

The cross section of pion-nucleon scattering is evaluated with the Breit-Wigner formula as the form of

σ_{π N \to R} (\sqrt{s}) = σ_{max} {(| p_{0} / p |)}^{2} \frac{0.25 Γ^{2} (p)}{0.25 Γ^{2} (p) + {(\sqrt{s} - m_{0})}^{2}},

(36)

where the p and p₀ are the momenta of pions at the energies of $\sqrt{s}$ and m₀, respectively, and m₀ being the centroid of resonance mass, e.g., 1.232 GeV, 1.44 GeV, and 1.535 GeV for Δ(1232), N^*(1440), and N^*(1535), respectively. The maximum cross section, σ_max, is taken from fitting the total cross sections of the available experimental data. For example, 200 mb, 133.3 mb, and 66.7 mb for π⁺+p→ Δ⁺⁺ (π^-+n→ Δ^-), π⁰+p→ Δ⁺ (π⁰+n→ Δ⁰) and π^-+p→ Δ⁰ (π⁺+n→ Δ⁺), respectively. And 24 mb, 12 mb, 32 mb, 16 mb for π^-+p→ N^*0(1440) (π⁺+n→ N^{* +}(1440)), π⁰+p→ N^{* +}(1440) (π⁰+n→ N^{* 0}(1440)), π^-+p→ N^{* 0}(1535) (π⁺+n→ N^{* +}(1535)), and π⁰+p→ N^{* +}(1535) (π⁰+n→ N^{* 0}(1535)), respectively. Figure 8 is a comparison of the elastic and total cross sections in the pion-nucleon collisions with the available data from the PDG collaboration [92]. The total cross sections include the sum of three resonances and the contributions of strangeness production. It is obvious that more resonances are needed to be implemented at the pion momentum above 0.5 GeV/c. The spectra can be reproduced nicely with the resonance approach, in particular in the domain of Δ resonance momenta (0.298 GeV/c). The pion-nucleon scattering dominates the energy deposition in the pion-induced reactions, which contributes to the fragmentation of target nuclei [93].

Figure 8:

(Color online) The elastic (upper panel) and total (down panel) pion-nucleon scattering cross sections contributed from different resonances. The available data are taken from the PDG collaboration [92].

The strangeness and vector mesons (ρ, ω) are created in inelastic hadron-hadron collisions without intermediate resonances. These are the 2-body and 3-body reaction channels as follows:

\begin{array}{l} B B \to B Y K, B B \to B B K \bar{K}, B π (η) \to Y K, Y K \to B π, \\ B π \to N K \bar{K}, Y π \to B \bar{K}, B \bar{K} \to Y π, Y N \to \bar{K} N N, \\ N N \to N N ρ, N N \to N N ω . \end{array}

(37)

Here the B stands for (N, △, N^*) and Y(Λ, ∑), K(K⁰, K⁺), and $\bar{K}$ ( ${\bar{K}}^{0}$ , K^-). The parameterized cross sections of each isospin channel BB → BYK [94] are used in the calculation. We take the parametrizations of the channels Bπ → YK [95, 96] besides the Nπ → Λ K reaction [97]. The results are close to the experimental data at near threshold energies. The cross section of antikaon production in inelastic hadron-hadron collisions is taken as the same form of the parametrization used in the hadron string dynamics (HSD) calculations [98]. Furthermore, the elastic scattering and strangeness-exchange reaction between strangeness and baryons have been considered through the channels of KB → KB, YB → YB, and $\bar{K} B \to \bar{K} B$ , and we use the parametrizations in Ref. [99]. The charge-exchange reactions between the KN → KN and YN → YN channels are included by using the same cross sections with the elastic scattering, such as K⁰p→ K⁺n, K⁺n→ K⁰p etc [100]. The cross sections for ρ and ω production are taken from the fitting of experimental data [101]. Figure 9 is the elementary cross section of different channels associated with the strangeness reactions. It is obvious that the Δ-nucleon collisions are dominant in the production of strange particles. The large absorption cross section in the K^--nucleon collisions retards the K^- emission in near-threshold energy heavy-ion collisions, but enhances the hyperon production. The pion induced reactions have considerable contributions on the strangeness production, roughly, one third to one fourth K⁺ yields being from the pion channels in heavy-ion collisions around the threshold energies [68].

Figure 9:

Elementary cross sections in the production of strange particles.

In order to describe the antiproton-nucleus collisions we have further included the annihilation channels, charge-exchange reaction (CEX), elastic (EL), and inelastic scattering as follows [70, 102]:

\begin{array}{l} \bar{B} B \to annihilation (π, η, ρ, ω, K, \bar{K}, η', K^{*}, {\bar{K}}^{*}, ϕ), \\ \bar{B} B \to \bar{B} B (CEX,EL), \bar{N} N \leftrightarrow \bar{N} Δ (\bar{Δ} N), \bar{B} B \to \bar{Y} Y . \end{array}

(38)

Here the B stands for nucleon and Δ(1232), Y(Λ, ∑, Ξ), K(K⁰, K⁺), and $\bar{K}$ ( $\bar{K^{0}}$ , K^-). The overline of B (Y) means its an antiparticle. The cross sections of these channels are based on the parametrization or extrapolation from available experimental data [103]. The annihilation dynamics in antibaryon-baryon collisions is described by a statistical model with SU(3) symmetry inclusion of all pseudoscalar and vector mesons [104]. Pions are the dominant products in the antiproton annihilation on a nucleus. The energy released in the antiproton-nucleon annihilation is mainly deposited in the target nucleus via the pion-nucleon collisions. Hyperons are mainly produced via strangeness exchange reactions in collisions of antikaons and nucleons, which have smaller relative momentum and could be easily captured by the residue nuclei to form hypernuclei. The nuclear fragmentation, hypernuclide formation, particle production, unexpected neutron/proton ratio and isospin effect in antiproton induced reactions have been thoroughly investigated within the LQMD model [102, 105, 106].

3 Nuclear dynamics and isospin effect in heavy-ion collisions

The LQMD transport model has been developed for many issues in nuclear dynamics, i.e., the neck fragmentation and isospin diffusion in quasi-fission dynamics from Fermi-energy heavy-ion collisions, high-density symmetry energy, in-medium properties of hadrons etc. Motivation of the future experiments is discussed on the topics of hypernuclide production, strangeness dynamics, isospin physics, and nuclear phase diagrams in heavy-ion collisions and hadron induced reactions. Parts of the results are presented in the article.

3.1 Neck fragmentation in isotopic nuclear reactions

Heavy-ion collisions at the Fermi energies (10-100A MeV) attract much attention on several topical issues in nuclear physics, i.e., spinodal multifragmentation, liquid-gas phase transition, properties of highly excited nuclei, symmetry energy in the domain of subnormal densities etc [5, 6, 107-112]. The isospin dynamics in the Fermi-energy heavy-ion collisions is related to several interesting topics, e.g., in-medium nucleon-nucleon cross sections, cluster formation and correlation, isotopic fragment distribution, the symmetry energy in dilute nuclear matter etc. Different mechanisms can coexist and are correlated in heavy-ion collisions at the Fermi energies, such as projectile or target fragmentation, neck emission, preequilibrium emission of light clusters (complex particles), fission of heavy fragments, multifragmentation etc, in which the isospin dependent nucleon-nucleon potential dominates the dynamical processes. The time scales from dynamical and pre-equilibrium emissions to statistical decay of excited systems at equilibrium and the isospin effect in neck fragmentation were investigated experimentally [113, 114]. The LQMD transport model is used to describe the nuclear dynamics in Fermi-energy heavy-ion collisions.

The composite system formed in Fermi energy heavy-ion collisions is highly excited with an excitation energy up to several tens of MeV/nucleon. The heated system is unstable and fragmentation even multifragmentation takes place. The fragmentation dynamics in the Fermi-energy heavy-ion collisions is described by the LQMD model. The primary fragments are constructed in phase space with a coalescence model, in which nucleons at freeze-out are considered to belong to one cluster with the relative momentum smaller than P₀ and with the relative distance smaller than R₀ (here P₀ = 200 MeV/c and R₀ = 3 fm). At the freeze-out, the primary fragments are highly excited. The de-excitation of the fragments is described within the GEMINI code [115]. Particles produced from the neck fragmentations in heavy-ion collisions could be probes of the low-density phase diagram, which are constrained within the midrapidities (|y/y_proj|<0.3) in semicentral nuclear collisions. Figure 10 is the kinetic energy spectra of neutron/proton ratios from the yields of free nucleons and ’gas-phase’ nucleons (nucleons, hydrogen, and helium isotopes) from the neck fragmentations in the ¹¹²Sn+⁵⁸Ni and ¹²⁴Sn+⁶⁴Ni reactions at a beam energy of 35 MeV/nucleon with the different symmetry energies [116]. A larger value of the n/p ratio with softening symmetry energy is found, in particular at the kinetic energies above the Fermi energy (36 MeV), which is caused from the fact that more repulsive force is enforced on neutrons in dilute nuclear matter. Calculations are performed without the momentum dependent interaction, and the set of parameter of PAR2 is taken.

Figure 10:

(Color online) Kinetic energy spectra of neutron/proton ratios from the yields of free nucleons ((a) and (b)) and ’gas-phase’ nucleons (nucleons, hydrogen, and helium isotopes) ((c) and (d)) from neck fragmentations in the ¹¹²Sn+⁵⁸Ni and ¹²⁴Sn+⁶⁴Ni reactions at the fermi energy of 35 MeV/nucleon within the collision centralities of 6-8 fm.

Besides the fast nucleons being probes of symmetry energy in the dilute matter, the neutron to proton ratio of light intermediate mass fragments (IMFs) could be sensitive to the stiffness of symmetry energy because of the isospin migration from the neck fragmentation [117]. The light IMFs (Z≤10) are measured by the CHIMERA detector at the INFN-LNS Superconducting Cyclotron of Catania (Italy) and emitted preferentially towards the midrapidity domain on a short timescale within 50 fm/c, which can not be entirely described through the decay of the excited projectile-like (PLF) and target-like (TLF) fragments [118]. We constrained the particles emitted from the neck fragmentation within the rapidity range of |y/y_proj|<0.3. Figure 11 is the average neutron to proton ratio of light IMFs in the isotopic reactions of ¹¹²Sn+⁵⁸Ni, ¹²⁴Sn+⁶⁴Ni, ¹¹²Sn+¹¹²Sn, and ¹²⁴Sn+¹²⁴Sn. The isospin effect is pronounced for the light isotopes, i.e., lithium, beryllium, and boron. The soft symmetry energy leads to a larger n/p ratio of light IMFs. It is caused from the fact that the heavy fragments tend to the neutron-rich side with the hard symmetry energy. However, the statistical decay leads to the disappearance of the isospin effect in heavy fragments. The light IMFs from the neck fragmentation are emitted preferentially towards the midrapidity domain on a short time scale in comparison to PLFs and TLFs. The isospin ratios depend on the stiffness of the symmetry energy and the effects increase with softening of the symmetry energy, in particular in neutron-rich nuclear reactions. It should be mentioned that the hard symmetry energy leads to the neutron-rich fragment formation in the PLF and TLF regime.

Figure 11:

(Color online) The average neutron to proton ratio of IMFs emitted within the rapidity range of |y/y_proj|<0.3 as a function of charged number with the different symmetry energies. The data from the CHIMERA detector [118] are shown for comparison.

The spectra of the isospin ratios are influenced by both the symmetry energy and the isospin splitting of nucleon effective mass, in particular at the high kinetic energies, which confuse the extraction of the density dependence of symmetry energy. On the other hand, the n/p ratios are also influenced by the Coulomb potential and the detector efficiencies of protons and neutrons in experiments. To eliminate the uncertainties, the double ratios of two isotopic systems would be nice probes for constraining the isospin splitting of nucleon effective mass and the symmetry energy beyond the saturation density from the experimental data. Figure 12 is the double ratio spectra in collisions of ¹²⁴Sn+¹²⁴Sn over ¹¹²Sn+¹¹²Sn. The influence of the isospin splitting with a linear symmetry energy (left panels) and the stiffness of symmetry energy with the mass splitting of $m_{n}^{*} > m_{p}^{*}$ (right panels) on the spectra is compared with the new data from the MSU-NSCL (National Superconducting Cyclotron Laboratory) [119]. It should be noticed that the effective mass splitting of the neutron and proton in the nuclear medium is pronounced in the kinetic energies above 30 MeV/nucleon, and the splitting of $m_{n}^{*} < m_{p}^{*}$ is nicely consistent with the available data at the beam energy of 120 MeV/nucleon. However, the difference of both the mass splittings on the spectra is very small at the beam energy of 50 MeV/nucleon. The symmetry energy effect is obvious and appears at the kinetic energy above 20 MeV/nucleon. A soft symmetry energy (γ_s=0.5) is constrained from both of the incident energies [120]. It is concluded that the mass splitting from the isovector momentum dependent interaction potential is negligible in the Fermi-energy heavy-ion collisions. The stiffness of symmetry energy is a dominant quantity on the double ratio spectra. However, the n/p ratios in high-energy heavy-ion collisions are influenced by both the symmetry energy and isospin splitting of effective mass. One should pay attention to some observables independent on the isospin splitting for constraining the density dependence of symmetry energy, such as the charged pion spectra from the high-density domain [121]. It should be noticed that the fragmentation reactions have been investigated within the antisymmetrized molecular dynamics model for extracting the temperature from the isobaric yield ratios [122].

Figure 12:

(Color online) The double neutron to proton ratios of gas-phase particles in collisions of ¹²⁴Sn+¹²⁴Sn over ¹¹²Sn+¹¹²Sn at the incident energy of 50 MeV/nucleon (upper panels) and 120 MeV/nucleon (lower panels) with different isospin splitting of effective mass and different stiffness of symmetry energy. The available data were measured at NSCL [119].

3.2 Pseudoscalar meson production and in-medium effects in heavy-ion collisions

Hadronic matter is created in high-energy heavy-ion collisions and exists in the compact stars, such as neutron stars. The hadron-hadron interaction in the superdense matter is complicated and varying with the baryon density. The hadron-nucleon potential impacts the ingredients of hadrons in the compact stars. In nuclear reactions, i.e., heavy-ion collisions, antiproton (proton) induced reactions etc, the production and phase-space distribution of particles were modified in the nuclear medium. The isospin dependence of the potential and the corrections on threshold energies influence the ratios of isospin particles. Consequently, the extraction of high-density symmetry energy is to be moved from particle production. On the other hand, the yields of particles and bound fragments such as hypernuclides, kaonic nucleus, and antiprotonic nucleus are related to the potential of particles in the nuclear medium. Dynamics of pseudoscalar mesons (π, η, K, and $\bar{K}$ ) and hyperons (Λ and ∑) produced in heavy-ion collisions near threshold energies has been investigated within the LQMD transport model. The in-medium modifications on particle production in dense nuclear matter are implemented in the model through corrections on the elementary cross sections and by inclusion of the meson (hyperon)-nucleon potentials, in which the isospin effects are considered. Pion is the lightest meson in the nature and can be produced in nucleon-nucleon collisions with a lower threshold energy e.g., E_th=289 MeV for π^±. The stochastic processes in heavy-ion collisions accelerate part of nucleons, which make the particle production at sub-threshold energies. Dynamics of pions produced in heavy-ion collisions have been investigated both in theories and in experiments, which is related to the issues of the high-density symmetry energy and the pion-nucleon interaction in dense matter. The production rate and density profile of π, η and K⁺ produced in central ¹⁹⁷Au+¹⁹⁷Au collisions at the incident energies of 0.4 and 1.5 A GeV are shown in Fig. 13. The reabsorption process retards the pion production towards the subnormal densities. The maximal yields correspond to the densities of 1.0ρ₀ and 1.6ρ₀ with and without the pion-nucleon (πN) scattering, respectively. The average density in the production pion is evaluated by $〈 ρ 〉 = \int ρ d N / N_{total}$ with N_total being the number of total pions. The average density are 1.07ρ₀ and 1.44ρ₀ for the cases with and without the πN absorption reactions at the energy of 0.4 A GeV, respectively. Therefore, pions produced in heavy-ion collisions provide the information of the phase diagram around saturation densities. The production of η is similar to the pion dynamics. However, K⁺ is created within a narrower time range at suprasaturation densities. The average density of the total kaons is around 2ρ₀ at the energy of 1.5 A GeV. One might evaluate the π^-/π⁺ ratio squeezed out the reaction zone for extracting the high-density symmetry energy.

Figure 13:

(Color online) Evolution of production rate and density profiles of particles in collisions of ¹⁹⁷Au+¹⁹⁷Au.

Dynamics of pseudoscalar mesons in heavy-ion collisions near threshold energies provide the isospin effect of nucleon force and the in-medium properties of u, d, and s quark in dense matter. The strange particles are produced at the supra-saturation densities and emitted at the early stage of the phase diagram when compared with the pion production. The strangeness exchange process retards the saturation of Λ and ∑ production [123]. Higher incident energy and central collisions enhance the domain of the high-density phase diagram, hence leading to the increase of the strangeness multiplicity. A larger high-density region of compressed nuclear matter is formed in heavy collision systems when compared with the light systems, which increases the strangeness production. The total multiplicity of particles produced in central ⁴⁰Ca+⁴⁰Ca collisions are calculated and shown in Fig. 14. The η and strange particles (K, $\bar{K}$ , Λ, ∑) are enhanced owing to the π induced reactions, such as π N→ N * (1535), π N→ KY and $π N \to N K \bar{K}$ . The available data from the FOPI collaboration for pions [124] and from theTAPS collaboration for etas [125] can be well reproduced. It should be mentioned that the inclusion of the KN potential in the model leads to about a 30% reduction of the total kaon yields in the subthreshold domain [68]. The competition of the π (η)N and KN potentials results in a bit of increase of kaon production because of the contribution of π (η)N→ KY.

Figure 14:

Excitation functions of particles produced in central ⁴⁰Ca+⁴⁰Ca collisions. The available data from FOPI collaboration for pions [124] and from TAPS collaboration for etas [125] are shown for comparison.

To investigate the high-density symmetry, the isospin observables emitted from the high-density reaction zone are expected. However, the dynamics of particles produced in the high-density reaction domain is distorted in the nuclear medium. The in-medium corrections on particle production in heavy-ion collisions are the elementary cross section modifications and the dynamical evolutions by the mean-field potentials [69]. The interaction potential of the strange particle and nucleon is of significance in the formation of a hypernucleus, the core structure of a neutron star, etc. However, it is not well understood up to now, in particular, in the dense nuclear matter. Consequently, the KN potential reduces the kaon production at midrapidities and at high transverse momenta. However, an opposite contribution of the ηN potential is obtained because of the attractive interaction of eta and nucleon in the nuclear medium. The hyperon-nucleon interaction is negligible for Λ dynamics. Figure 15 is the inclusive spectra of π⁰, η, K⁰, ${\bar{K}}^{0}$ , and neutral hyperons (Λ+∑⁰) in ¹⁹⁷Au+¹⁹⁷Au collisions. Similar to the transverse momentum spectra, the eta-nucleon, kaon-nucleon, and hyperon-nucleon potentials weakly impact the particle emission. However, the interaction of pions (antikaons) and nucleons contributes to the pion (antikaon) dynamics, i.e., increasing the antikaon production at low kinetic energies. More pronounced information of the in-medium effects has been investigated from the K^-/K⁺ spectrum in heavy-ion collisions. A deeply attractive K^-N potential, being the value of -110±15 MeV was obtained at saturation density, and the weakly repulsive K⁺N potential has been concluded to be 25±10 MeV from the heavy-ion collisions [8, 34, 36, 37]. The in-medium modifications influence the kaon production and dynamical emissions in phase space, i.e., inclusive spectrum, collective flows etc. We compared the effects of kaon(antikaon)-nucleon potentials in the ¹²C+¹²C reaction as a function of kinetic energy and in the collisions of ¹⁹⁷Au+¹⁹⁷Au and proton on ¹²C and ¹⁹⁷Au versus the transverse mass ( $m_{t} = \sqrt{p_{t}^{2} + m_{0}^{2}}$ with p_t being the transverse momentum and the mass of kaon (antikaon) being m₀) as shown in Fig. 16 [126]. The experimental data from the KaoS collaboration [127-129] can be nicely reproduced within the inclusion of the mean-field potentials for kaons and antikaons, in which a weakly repulsive KN potential of the order of 28 MeV and a deeply attractive $\bar{K} N$ potential of -100 MeV at normal nuclear density are used in the LQMD model. The attractive $\bar{K} N$ potential reduces the threshold energies associated with enhancing K^- production. However, the KN potential leads to an opposite contribution for the K⁺ emission. The value of the K^-/K⁺ ratio is more sensitive to the optical potentials because of the opposite interactions of kaon and antikaon in nuclear medium. The K^--bound state and kaon condensation in dense baryonic matter are being searched in experiments.

Figure 15:

Inclusive spectra of the neutral particles produced in ¹⁹⁷Au+¹⁹⁷Au collisions and in-medium effect on the particle production [69].

Figure 16:

The ratio of K^-/K⁺ as a function of transverse mass (kinetic energy) in collisions of ¹²C+¹²C and protons on ¹²C and ¹⁹⁷Au at the beam energies of 1.8A GeV and 2.5 GeV, respectively. The reaction of ¹⁹⁷Au+¹⁹⁷Au at the energy of 1.5A GeV [126]. The experimental data are taken from the KaoS collaboration [127-129]

The high-density matter larger than normal nuclear density can be formed in high-energy heavy-ion collisions, where strange particles are produced. The influence of the in-medium potential by distinguishing the isospin effect and the stiffness of symmetry energy on strangeness production can be observed from the emission of isospin pairs. Figure 17 is a comparison of the energy dependence of the yield ratios of K⁰/K⁺ and ∑^-/∑⁺ produced in central ¹⁹⁷Au+¹⁹⁷Au collisions. It should be noticed that the KN potential increases the K⁰/K⁺ ratio by distinguishing isospin effect and also enlarges the influence of symmetry energy on isospin ratio in comparison to the in-vacuum and the nonisospin cases, in particular in the domain of deep subthreshold energies. The effect of symmetry energy on the isospin ratios disappears at high incident energy. The inclusion of the KN potential leads to a reduction of the total kaon yields. The contribution of LF is very small on the kaon production, but changes the kaon distribution in momentum space [68, 100]. Calculations of the double ratio excitations of K⁰/K⁺ taken from the two systems of ⁹⁶Ru+⁹⁶Ru and ⁹⁶Zr+⁹⁶Zr present that the symmetry energy effect is very weak over the whole energy range [100]. In analogue to the isospin ratio of kaons, the effect of symmetry energy on the ∑^-/∑⁺ excitation function is pronounced in the region of the deep threshold energies. Overall, a hard symmetry energy leads to a lower isospin ratio in neutron-rich nuclear collisions. The ∑^-/∑⁺ ratio is larger than the value of K⁰/K⁺ in the whole energy regime. Production of ∑ increases dramatically with the incident energy and the isospin effects of the yields for charged ∑ are obvious with decreasing the energy. Despite that part of sigma can be absorbed in nuclear matter through the reaction of $π Σ \to \bar{K} B$ , most ∑ are produced in the high-density region in heavy-ion collisions, which can be used to probe the information of the high-density phase diagram of nuclear matter. The KN and ∑N potentials are of importance on the kaon and ∑ emissions in phase space. The isospin ratio of K⁰/K⁺ depends on the kaon potential at high transverse momenta, and a flat spectrum appears after inclusion of the KN potential. The K⁰/K⁺ and ∑^-/∑⁺ ratios of neutron-rich heavy systems in the domain of subthreshold energies are sensitive to the stiffness of nuclear symmetry energy, which are created in the high-density domain in heavy-ion collisions and might be the promising probes to extract the high-density information of symmetry energy.

Figure 17:

Excitation functions of the isospin ratios of K⁰/K⁺ and ∑^-/∑⁺ in the ¹⁹⁷Au + ¹⁹⁷Au reaction for the hard and soft symmetry energies, respectively.

3.3 Hypernucleus formation in heavy-ion collisions

The formation mechanism of fragments with strangeness in high-energy heavy-ion collisions has been investigated within the LQMD transport approach combined with the statistical model (GEMINI). The production and dynamics of hyperons is described within the LQMD model. A coalescence model is used for constructing the primary fragments and the hyperon capture by residual nucleons. The combined approach is used to describe the formation of hyperfragments. The production of hypernuclei is associated to the hyperon production, hyperon-nucleon and hyperon-hyperon interactions, capture of hyperons by nucleonic fragments, and decay of excited hyper-fragments. The investigation of hypernucleus properties is an essential way for extracting the in-medium information of hyperons. To form the hyperfragments in heavy-ion collisions, the incident energy is chosen to be enough for creating hyperons, but not too high that the hyperons are captured by surrounding nucleons. At the near threshold energies, the reaction channels of BB → BYK and π(η)B → KY dominate the hyperon production. Usually, the hyperons are created in the domain of the dense nuclear medium. The hyperon-nucleon potential impacts the hyperon dynamics and hyperfragment formation. Calculations from a statistical model showed that the beam energy regime of 3-5 GeV/nucleon is available for producing hypernuclei [130]. The production and structure studies of neutron-rich and even double-strangeness hypernuclides have been planned at the high intensity heavy-ion facility (HIAF) in the future. Figure 18 is the rapidity and kinetic energy spectra of Λ-hypernuclides with the charged numbers of Z≥2 in collisions of ¹¹²Sn+¹¹²Sn, ¹²⁴Sn +¹²⁴Sn, and ¹³²Sn +¹²⁴Sn at the incident energy of 2 GeV/nucleon and within the collision centrality bin of 0-8 fm in the center of mass frame. The hyperfragments are formed within a narrower rapidity domain and at the less kinetic energies in comparison to hyperons in Fig. 15. The phase-space distributions are available for arranging the detector setup in experiments. The isotopic and charge distributions of Λ-hyperfragments are also calculated, as shown in Fig. 19. It is obvious that the neutron-rich systems are favorable for producing the more neutron-rich hyperfragments and the production rate decreases drastically with the charged number. The medium and heavy hyperfragments might be created with the peripheral collisions, but have very small cross sections owing to the less production number of hyperons. The strangeness physics opens a new window in nuclear physics, besides the hypernuclide properties, and also are the main ingredients in the dense baryonic matter, such as neutron stars etc. An accurate description of the hypernucleus formation in heavy-ion collisions is still in progress, i.e., the quantal clustering hyperfragments, the multi-mode decays of excited hyperfragments, the formation of exotic hypernucleus ( $_{Λ}^{3} n$ , $_{Λ}^{5} n$ ) etc.

Figure 18:

Rapidity and kinetic energy distributions of hypernuclides with the charged numbers of Z≥2 in the isotopic reactions.

Figure 19:

Isotopic and charge distributions of Λ-hypernuclides produced in the ^132,124,112Sn + ¹¹²Sn/¹²⁴Sn collisions at the incident energy of 2A GeV.

4 Nuclear dynamics induced by hadrons

The nuclear matter formed in heavy-ion collisions is compressed, and the nuclear density varies with the time-space evolutions and the incident energy. The spinodal and chemical instabilities complicate the nuclear fragmentation and particle emission in heavy-ion collisions. However, the nuclear reactions induced by hadrons have advantages on many topical issues, i.e., the in-medium properties of hadrons around the saturation density, the highly excited nucleus, the production of hypernucleus etc.

4.1 Particle production in proton-nucleus collisions

The nuclear dynamics induced by protons has many interesting aspects, in particular, the spallation reactions on heavy targets, the in-medium properties of hadrons, the hypernucleus formation etc. Microscopic description on the issues is quite necessary. It has been found that the strangeness production is strongly suppressed in proton-induced reactions in comparison to heavy-ion collisions [126]. The kaon-nucleon and antikaon-nucleon potentials change the structures of rapidity and transverse momentum distributions and the inclusive spectra in proton-nucleus collisions. Figure 20 is the temporal evolution of pions, resonances, η, and strange particles produced in collisions of a proton on ⁴⁰Ca at the momentum of 3 GeV/c. The pions are dominant products and the secondary collisions of π(η)N → KY are negligible for the hyperon production. Once the particles are created, they go away from the surrounding nucleons rapidly and the shadowing effect is very weak in the proton induced reactions. The strange particles are mainly from the direct process of the channel NN(Δ) → NYK. The phase-space distributions of particles produced in proton-nucleus collisions are distorted by the surrounding nucleons although the secondary collisions are slight. Therefore, it provides a good opportunity for studying the optical potentials around the normal nuclear density and the exotic bound states with the proton induced reactions, i.e., πN, ηN, YN etc.

Figure 20:

Temporal evolution of particle production in collisions of a proton on ⁴⁰Ca at the momentum of 3 GeV/c.

It has been found that the kaon- (antikaon-) nucleon potential plays a significant role on the strangeness production and dynamical emission in phase space, which reduces (enhances) the kaon (antikaon) yields and is pronounced when the incident energy is close to the threshold values. Although there are very low hyperon yields produced in proton-nucleus collisions, the hyperfragments could be formed because the hyperons are created inside the target nucleus. The fragments and hyperfragments are calculated, as shown in Fig. 21, with the 5 GeV/c protons on the targets of ⁴⁰Ca, ¹²⁴Sn, and ¹⁹⁷Au, respectively. Roughly, the Λ-hyperfragments are produced with the yields of the 4-order magnitude lower than the nucleonic fragments. The fragments tend to be formed in the target-like region. It provides the possibility for studying the delayed fission of heavy hypernuclide with the high-energy protons.

Figure 21:

Production of nucleonic fragments and hyperfragments in collisions of a proton on ⁴⁰Ca, ¹²⁴Sn, and ¹⁹⁷Au at the momentum of 5 GeV/c, respectively.

4.2 Meson induced nuclear reactions

The nuclear dynamics induced by charged mesons has attracted much attention on the aspects of the charge-exchange reactions, nuclear fragmentation, hypernuclide production etc. On the other hand, the meson-nucleus collisions have advantages in investigating the energy dissipation mechanism, the meson-nucleon interaction, and resonance properties in nuclear medium around the saturation density. Figure 22 is the comparison of the nucleonic fragments, hyperfragments, and free Λ in the K^-+⁴⁰Ca reaction at an incident momentum of 2 GeV/c. A similar structure of hyperfragments and free Λ is found. The strangeness exchange reaction of K^-N → π Y dominates the hyperon production. Once the Λ is created in the nucleus, a large probability is captured to form a hypernucleus. The free Λ hyperons and hyperfragments have the narrower rapidity regime and lower kinetic energies. The nuclear fragmentation is still pronounced owing to the decay of the highly excited target nucleus.

Figure 22:

Rapidity and kinetic energy distributions of nucleonic fragments, hyperfragments and free Λ in the K^-+⁴⁰Ca reaction at an incident momentum of 2 GeV/c.

Nuclear reactions induced by pions provide the opportunity to study the pion-nucleon interaction, the charge-exchange reactions, decay modes of highly excited nucleus, hypernucleus formation, and the in-medium properties of Δ resonance. Recently, the bound state of the neutron-rich Λ-hypernucleus, $Λ_{6}^{6}$ H, was studied with the double charge-exchange reaction ⁶Li(π^-,K⁺)X at a beam momentum of 1.2 GeV/c at J-PARC [131]. The nuclear fragmentations and the charge exchange reactions in pion-nucleus collisions near the Δ(1232) resonance energies has been investigated within the LQMD model [93]. It is found that the fragments tend to be formed along the β stability line and the fission reactions take place for heavy targets. The relative motion energy is deposited in the nucleus via the pion-nucleon collisions for the pion induced reactions. The dissipated energy weakly depends on the incident pion energy, which is available through the hypernuclide formation. Figure 23 is the distributions of fragment and hyperfragment produced in collisions of π⁺ on ¹⁹⁷Au at the momentum of 3 GeV/c. The medium fragments are from the fission of target-like nuclei. It is noticed that the Λ- hyperfragments are formed in the target-like region.

Figure 23:

Nucleonic fragment and hyperfragment production in collisions of π⁺ on ¹⁹⁷Au at the momentum of 3 GeV/c.

The hyperon does not obey the Pauli principle and can occupy the same quantum state with nucleons. On the other hand, the attractive hyperon-nucleon interaction enables the binding of the very neutron-rich hypernucleus and even the two-neutron-Λ or four-neutron-Λ state. The isotopic distributions of hyperfragments in collisions of K^- on ⁴⁰Ca at the momentum of 2 GeV/c are shown in Fig. 24. The maximal yields are located around the neutron number of N=19. A comparison of nucleonic fragments and hyperfragments is shown in Fig. 25. The decay of primary hyperfragments enables the nucleon emissions and lightens the hypernuclei. The hyperfragment yields are reduced by the two order magnitude than the nucleonic fragments, which are available for hypernucleus formation with the strangeness exchange reactions in comparison to heavy-ion collisions and π induced reactions.

Figure 24:

Isotopic distributions of hyperfragments in collisions of K^- on ⁴⁰Ca at the momentum of 2 GeV/c.

Figure 25:

Comparison of nucleonic fragments and hyperfragments in collisions of K^- on ⁴⁰Ca at the momentum of 2 GeV/c.

4.3 Hypernucleus formation in antiproton-nucleus collisions

Since the first evidence of antiprotons was found in 1955 at Berkeley in collisions of protons on copper at the energy of 6.2 GeV by Chamberlain, Segrè, Wiegand, and Ypsilantis [132], the secondary beams of antiprotons were produced at many laboratories, such as CERN, BNL, KEK etc [133-136]. The stochastic cooling method provides the possibility for storing the antiprotons produced in proton-nucleus collisions. The particles W^± and Z⁰ were found for the first time with the high energy protons colliding the stored antiprotons at CERN [137, 138]. On the other hand, the antiproton-nucleus collisions are motivated by many interesting issues, i.e., charmonium physics, strangeness physics, antiprotonic atom, symmetry, in-medium properties of hadrons, cold quark-gluon plasma, highly excited nucleus etc [50, 139]. Recently, the antiproton-antiproton interaction was investigated by the STAR collaboration in relativistic heavy-ion collisions [140]. The low-energy antiprotons usually annihilate at the nucleus surface because of the large absorption cross section. The huge annihilation energy is available for producing the 2-6 pions. The subsequent processes are complicated and also associated with the multiple pion-nucleon interaction, which result in the fragmentation of the target nucleus and the preequilibrium emissions of complex particles. The localized energy is deposited in the nucleus with an excitation energy of several hundreds of MeV. The hot nucleus proceeds to the explosive decay via multifragmentation process or the sequential particle evaporation. On the other hand, the collisions of the antiproton and secondary particles with surrounding nucleons lead to the pre-equilibrium particle emissions, which are related to the scattering cross sections of each reaction channels, antiproton-nucleon interaction, particle-nucleon potentials, and density profile of target nucleus. The unexpected large neutron yields produced by stopped antiprotons in nuclei were reported in the LEAR experiments [141].

The dynamics of the antiproton-nucleus collisions is more complicated in comparison to hadron (proton, π, K etc) induced reactions and to heavy-ion collisions, in which the particles produced in the annihilation of the antiproton in a nucleus are coupled to the subsequent collisions with surrounding nucleons. The dynamics of antiproton-nucleus collisions are complicated, which is associated with the mean-field potentials of hadrons in nuclear medium, and also coupled to a number of reaction channels, i.e., the annihilation channels, charge-exchange reaction, and elastic and inelastic collisions. A more localized energy deposition enables the secondary collisions to be available for producing hyperons. Hyperons produced in antiproton induced reactions can be captured in the potential of nucleon fragments to form hypernuclei. The dynamics of antiproton-nucleus collisions is complicated, which is associated with the mean-field potentials of hadrons in nuclear medium, and also with a number of reaction channels, i.e., the annihilation channels, charge-exchange reaction, and elastic and inelastic collisions. The larger yields of strange particles in antiproton induced reactions are favorable to form hypernuclei in comparison to proton-nucleus and heavy-ion collisions. To understand the nuclear dynamics induced by antiprotons, several approaches have been proposed, such as the intranuclear cascade (INC) model [142], kinetic approach [143], Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) transport model [144, 145], Statistical Multifragmentation Model (SMM) [146], and the Lanzhou quantum molecular dynamics (LQMD) approach [102, 147]. Particle production in antiproton induced reactions is significant in understanding the contributions of different reaction channels and the energy dissipation into the nucleus. In Fig. 26 the temporal evolutions of pions, kaons, antikaons, hyperons, and antihyperons in the reaction of $\bar{p}$ +⁴⁰Ca at an incident momentum of 4 GeV/c are shown. It is shown that the kaons are emitted immediately after the annihilation in collisions of antibaryons and baryons. The secondary collisions of pions and antikaons on nucleons retard the emissions and even a bump appears for antikaons in the reaction dynamics, i.e., π N→ Δ, $\bar{K} N \to π Y$ etc, which lead to the production of hyperons. At the considered momentum above its threshold energy, e.g., the reaction $\bar{N} N \to \bar{Λ} Λ$ (p_threshold=1.439 GeV/c), production of hyperons are attributed from the direct reaction (annihilation and creation of quark pairs, $u \bar{u} (d \bar{d}) \to s \bar{s}$ ) and also from the secondary collisions after annihilations in antibaryon-baryon collisions, i.e., meson induced reactions, π (η, ρ, ω) N→ KY, and the strangeness exchange reaction $\bar{K} N \to π Y$ . It is noticed that the antikaon yields are reduced for the heavy targets owing to more strong reabsorption reactions. However, the hyperon production is more pronounced with increasing the target mass, which is available for producing hypernucleus.

Figure 26:

(Color online) Temporal evolution of particle production in collisions of

\bar{p}

on ⁴⁰Ca at incident momentum of 4 GeV/c.

The fragmentation process of the target nucleus is induced by an antiproton undergoing the explosive process (fast stage) in which the preequilibrium nucleons (light clusters) are emitted or the multifragments are produced after collisions between baryons and nucleons, and the decay process (slow stage) of the highly excited nucleus after the relative motion energy is deposited via the meson-nucleon collisions. The decay mechanism is determined by the excitation energy, i.e., the particle evaporation or fission dominating at low excitation energies (1-2 MeV/nucleon). The system is broken via multifragmentation when the local energy is close to the binding energy. More information on the hyperfragment formation in antiproton induced reactions is pronounced from the mass and charged number distributions. Direct production of hypernuclei with strangeness s=-2 (double Λ- hypernucleus) in heavy-ion collisions or by proton induced reactions are difficult because of very small cross sections, in particular for the heavy-mass region. Properties of the hypernuclei would be significant in understanding the Λ-Λ and $\bar{Λ}$ -nucleon interactions, which are not well understood up to now. The (K^-, K⁺) reactions are used for producing the double hypernucleus, $_{Λ Λ}^{6}$ He [148]. The antiproton-nucleus collisions would be a chance for producing the s=-2 and s=1 hypernuclei. Figure 27 is the mass and charged number distributions of hyperfragments with strangeness s=-1 (_ΛX and _∑±X) and s=-2 (_ΛΛX) in collisions of $\bar{p}$ + ⁶³Cu and $\bar{p}$ + ¹⁸¹Ta at the same of incident momentum of 4 GeV/c. The Λ-hyperfragments spread the whole isotope range with the larger yields. The maximal cross sections for s=-1 and s=-2 hyperfragments are at the levels of 1 mb and 0.01 mb, respectively. The lower production yields of $\bar{Λ}$ -hyperfragments at the level of 1 μb are found [102]. The hypernuclear physics with the antiproton beams is planned at the FAIR-PANDA (Antiproton Annihilation at Darmstadt) in the near future.

Figure 27:

(Color online) Hyperfragments with strangeness s=-1 (_ΛX and _∑±X) and s=-2 (_ΛΛX) spectra as functions of atomic number (left panel) and charged number (right panel) in collisions of

\bar{p}

on ⁶³Cu and ¹⁸¹Ta at incident momentum of 4 GeV/c.

5 Summary

In this article, I reviewed the recent progress on the nuclear dynamics in heavy-ion collisions and in hadron induced reactions, in particular, with the topics of symmetry energy, in-medium properties of hadrons, fragmentation reactions, and hypernucleus formation being stressed. The possible challenges and treatments in simulating nucleus-nucleus collisions with transport models are discussed. Further experiments are motivated at the facilities in the world such as CSR, HIAF, FRIB, Spiral2, FAIR etc. Within the framework of the LQMD transport model, the high-density symmetry energy, decay mode of highly excited nucleus, hyperon-nucleon and hyperon-hyperon interactions in dense nuclear matter, and hypernuclide production with double strangeness are thoroughly investigated.

The isospin effect in Fermi-energy heavy-ion collisions is pronounced, in particular from the neck fragmentation. The ratios of neutron to proton yields from free nucleons and from light IMFs are sensitive to the stiffness of symmetry energy. The soft symmetry energy at the domain of subsaturation densities is constrained from the CHIMERA and MSU data. The mass splitting of neutron and proton in the nuclear medium is obvious at the high beam energy. The optical potentials of pions, etas, and kaons influence the structures of invariant spectra. The hadron production is coupled to the reaction channels and the transportation is related in the mean-field potentials. The phase-space distributions of particles produced in hadron-hadron collisions are influenced by the potentials. A weakly repulsive KN potential roughly being 25 MeV and the attractive $\bar{K} N$ potential being -100 MeV at the normal density are obtained in comparison to the KaoS data. The pion-nucleon potential is not well constrained up to now, which impacts the pion energy spectra and associated with the properties of Δ(1232) in the nuclear medium. Parts of pions and etas are created at the subsaturation densities because of the reabsorption process. The rapidity and azimuthal angle cuts are necessary to extract the high-density symmetry energy from the yield ratios of π^-/π⁺. Further experiments on the pion measurements are expected, such as HIRFL-CSR at Lanzhou, RIKEN-SAMURAI etc.

The nucleonic fragments and hyperfragments are investigated within the LQMD transport model and are associated with the GEMINI statistical code. Hyperfragments are formed in the region of the narrower rapidity and lower kinetic energy in comparison with nucleonic fragments and hyperons. The hyperfragment yields formed in heavy-ion and proton (pion)-nucleus collisions are 10^-4 times smaller in comparison to the fragments. The nuclear reactions induced by K^- and antiprotons are available for producing the Λ-hypernuclide. Moreover, the heavy-ion collisions, in particular with the neutron-rich beams, have an advantage for creating the neutron-rich hypernuclide. Experiments on the strangeness physics and hypernucleus properties are planned at the HIAF in the future.

References

[1]

E.A. Uehling, G.E. Uhlenbeck.

Transport phenomena in Einstein-Bose and Fermi-Dirac gases

. Phys Rev, 1933, 43: 552. doi: 10.1103/PhysRev.43.552