Collective flow and nuclear stopping in heavy ion collisions in Fermi energy domain

Special Section on International Workshop on Nuclear Dynamics in Heavy-Ion Reactions (IWND2018)

Collective flow and nuclear stopping in heavy ion collisions in Fermi energy domain

Peng-Cheng Li，

Yong-Jia Wang，

Qing-Feng Li，

Hong-Fei Zhang

Nuclear Science and Techniques

Vol.29, No.12

Article number 177

Published in print 01 Dec 2018

Available online 17 Nov 2018

DOI：10.1007/s41365-018-0510-1

65102

The effects of the in-medium nucleon-nucleon (NN) elastic cross section on the observables in heavy ion collisions in the Fermi energy domain are investigated within the framework of the ultrarelativistic quantum molecular dynamics model. The results simulated using medium correction factors of $F = σ_{N N}^{in - medium} / σ_{N N}^{free} = 0.2, 0.3, 0.5$ and the density- and momentum-dependent factor obtained from the FU3FP1 parametrization are compared with the FOPI and INDRA experimental data. It is found that the calculations using the correction factors $F$ = 0.2 and 0.5 reproduce the experimental data (i.e., collective flow and nuclear stopping) at 40 and 150 MeV/nucleon, respectively. Calculations with the FU3FP1 parametrization can best fit these experimental data. These conclusions can be confirmed in both ¹⁹⁷Au+¹⁹⁷Au and ¹²⁹Xe+¹²⁰Sn.

Nucleon-nucleon elastic cross sectionHeavy ion collisionsNuclear stoppingCollective flow.

1 Introduction

The properties of hot, dense nuclear matter are usually extracted by comparing the observables in heavy ion collision (HIC) experiments with the corresponding results from transport model simulations. The most popular transport methods for studying HICs at low and intermediate energies are the quantum molecular dynamics (QMD) model [1], the Boltzmann (Vlasov)–Uehling–Uhlenbeck (BUU or VUU) model [2], and their improved versions [3-5]. The mean-field potential and nucleon-nucleon cross section (NNCS) are two essential components of these models. The mean-field potential in transport models has been studied extensively [6-8]. Regarding the NNCS, it is well-known that compared with the free NNCS, the in-medium NNCS is suppressed [9-24]; however, the degree of suppression is still uncertain and requires further improvement. Thus, the in-medium NNCS has been studied by many groups using different methods [25-34].

Recently, a systematic experimental study of nuclear stopping in central collisions at intermediate energies pointed out that the correction factor $F = σ_{N N}^{in - medium} / σ_{N N}^{free}$ is deduced to be 0.16±0.04 at 40 MeV/nucleon and 0.5±0.06 at 100 MeV/nucleon [35]. New opportunities emerged to study the in-medium NNCS using transport models based on these stopping data [36-39]. There are many parameterized in-medium NNCSs adopted by various models. For example, $σ_{N N}^{in - medium} = (1 - η ρ / ρ_{0}) σ_{N N}^{free}$ with η=0.2 has been adopted in the BUU, relativistic BUU, and improved QMD models [10, 18]. In the pBUU model, the in-medium NNCS is independent of the density, energy, and isospin, and reads as $σ_{N N}^{in - medium} = 0.85 ρ^{- 2 / 3} / \tanh (\frac{σ^{free}}{0.85 ρ^{- 2 / 3}})$ [21]. In the isospin-dependent BUU model and the Lanzhou QMD model [40], the NNCS in the nuclear medium is corrected by a factor of $F = σ_{N N}^{in - medium} / σ_{N N}^{free} = (μ_{N N}^{*} / μ_{N N})^{2}$ . Here μNN and $μ_{N N}^{*}$ are the κ masses of the colliding nucleon partners in free space and in the nuclear medium, respectively [16, 22, 41]. In our previous works, we found that the collective flow and nuclear stopping can be reproduced using the parameterized medium correction factor $F$ (which depends on the density and momentum) for the free NNCS within the ultrarelativistic QMD (UrQMD) model [42-44]. It is desirable to evaluate the differences between the $F$ values extracted from experimental data [35] and the parameterized one adopted currently in the UrQMD model.

This paper is organized as follows. In Sec. 2, we review the UrQMD model with the improved potential. Sec. 3 presents an investigation of the effects of the in-medium NNCS on the observables of free protons and hydrogen isotopes in HICs in the Fermi energy domain. Finally, a summary is given in Sec. 4.

2 Model description and observables

The UrQMD model is based on the same principles as the QMD model. In the UrQMD model, more than 55 different baryon and 32 different meson species, as well as the corresponding antiparticle and isospin-projected states, are considered. The model has been extensively and successfully used to study the nuclear reactions of p + p, p + A, and A + A systems within a large range of beam energies, from the low-energy regime of the INDRA/GSI experiments up to the highest energies presently available at LHC/CERN. In addition, the other two main features unique to the collision term of the UrQMD model are the unique collision time for each individual collision and the two-step particle production process. More discussions can be found in recent works on the transport model comparison project, i.e., Refs. [3, 4].

At intermediate energies, the potential used in the UrQMD model depends on the momentum and density [1, 42, 45], and reads as

U = α \cdot (\frac{ρ}{ρ_{0}}) + β \cdot {(\frac{ρ}{ρ_{0}})}^{γ} + t_{md} \ln^{2} [1 + a_{md} {(p_{i} - p_{j})}^{2}] \frac{ρ}{ρ_{0}} .

(1)

Here α = -393 MeV, β = 320 MeV, γ = 1.14, t_md = 1.57 MeV, and a_md=0.0005 MeV^-2; these values yield an incompressibility K₀ of 200 MeV for isospin-symmetric nuclear matter. In order to better describe the recent experimental data, the surface, surface asymmetry energy, and symmetry energy terms obtained from the Skyrme potential energy density functional have been further introduced into the present version [44, 46, 47]:

\begin{array}{l} u_{Skyrme} & = u_{sur} + u_{sur,iso} + u_{sym} \\ = \frac{g_{sur}}{2 ρ_{0}} {(\nabla ρ)}^{2} + \frac{g_{sur,iso}}{2 ρ_{0}} {[\nabla (ρ_{n} - ρ_{p})]}^{2} \\ + (A_{sym} \frac{ρ^{2}}{ρ_{0}} + B_{sym} \frac{ρ^{η + 1}}{ρ_{0}^{η}} + C_{sym} \frac{ρ^{8 / 3}}{ρ_{0}^{5 / 3}}) δ^{2} . \end{array}

(2)

\frac{g_{sur}}{2} = \frac{1}{64} (9 t_{1} - 5 t_{2} - 4 x_{2} t_{2}) ρ_{0},

(3)

\frac{g_{sur,iso}}{2} = - \frac{1}{64} [3 t_{1} (2 x_{1} + 1) + t_{2} (2 x_{2} + 1)] ρ_{0},

(4)

A_{sym} = - \frac{t_{0}}{4} (x_{0} + 1 / 2) ρ_{0},

(5)

B_{sym} = - \frac{t_{3}}{24} (x_{3} + 1 / 2) ρ_{0}^{η},

(6)

C_{sym} = \frac{1}{24} {(\frac{3 π^{2}}{2})}^{2 / 3} ρ_{0}^{5 / 3} Θ_{sym},

(7)

where Θ_sym=3t₁x₁-t₂(4+5x₂) [44]. In this work, the SV-sym34 force, which gives g_sur=18.2 MeV fm², g_sur,iso=8.9 MeV fm², A_sym=20.3 MeV, B_sym=14.4 MeV, C_sym=-9.2 MeV, and the slope parameter of the symmetry energy L=80.95 MeV, is chosen.

In the UrQMD model, the in-medium NN elastic cross section is the product of the free cross sections with a medium correction factor $F$ (ρ,p), which is momentum- and density-dependent. In the UrQMD model, the free proton-neutron cross sections are larger than the free neutron-neutron (proton-proton) cross sections at low beam energy, which is consistent with experimental data. In our previous work [17, 33, 43], we found that the medium correction factor $F$ depends on the density, momentum, and isospin. However, the isospin dependence is relatively weak; see, e.g., [48, 49].

σ_{N N}^{in - medium} = F (ρ, p) * σ_{N N}^{free},

(8)

F (ρ, p) = \frac{λ + (1 - λ) e^{- ρ / ρ_{0} / ζ} - f_{0}}{1 + {(p_{N N} / p_{0})}^{κ}} + f_{0} .

(9)

Here λ=1/3, ζ=1/3, f₀=1, p₀=0.425 GeV/c, and κ=5, which corresponds to the FU3FP1 parametrization used in Ref. [43]. Moreover, $F$ (ρ,p) is set to unity when pNN is larger than 1 GeV/c.

The in-medium correction factor $F$ (ρ,p) obtained from the FU3FP1 parametrization as functions of both the reduced density ρ/ρ₀ and momentum p=pNN is shown in Fig.1. It is seen that the values deduced from the experiment [35] can be covered by the values of the in-medium correction factor $F$ (ρ,p) obtained from the FU3FP1 parametrization. In order to compare the effect of the in-medium NN cross section on the collective flow and nuclear stopping, three fixed values, $F$ = 0.2, 0.3, and 0.5, are further adopted in this work.

Figure 1:

(Color online) In-medium correction coefficient

F

obtained from the parametrization FU3FP1, which is density- and momentum-dependent.

3 Results and discussion

3.1 Charge distribution

In this work, fragments are recognized by the isospin-dependent minimum spanning tree (iso-MST) method. Nucleons with relative distances smaller than $R_{0}^{p p} = 2.8$ fm, $R_{0}^{n n} = R_{0}^{n p} = 3.8$ fm and relative momenta smaller than P₀ =0.25 GeV/c are considered to belong to the same cluster. The charge distribution of Au + Au collisions at a beam energy of 150 MeV/nucleon is displayed in Fig.2. The distribution of charged fragments decreased exponentially with increasing charge number, and this behavior can be reproduced by both the simulations with FU3FP1 and those with the fixed in-medium correction factors. For large fragments (Z ≥ 6), calculations with $F$ =0.5 seem to be closer to the experimental data. However, one cannot conclude that $F$ =0.5 is preferable to the others because the charge distribution will also be affected by unconstrained (or less constrained) parameters in the transport model; e.g., the chosen coalescence model parameters will strongly influence the charge distribution but not the collective flow and nuclear stopping [44, 50-57].

Figure 2:

(Color online) Charge multiplicity distribution in ¹⁹⁷Au+¹⁹⁷Au central collisions at beam energy of 150 MeV/nucleon. The experimental data are taken from Ref. [58]. The calculations performed with the in-medium correction factors of

F

=0.2 (circles),

F

=0.3 (up-facing triangles), and

F

=0.5 (down-facing triangles) are compared with calculations using the FU3FP1 parametrization of the in-medium NNCS (squares).

3.2 Collective flow

One can obtain the directed (v₁) and elliptic (v₂) flows from the Fourier expansion of the azimuthal distribution of detected particles [59, 60]; they read as

v_{1} \equiv 〈 \cos (ϕ) 〉 = 〈 \frac{p_{x}}{p_{t}} 〉,

(10)

v_{2} \equiv 〈 \cos (2 ϕ) 〉 = 〈 \frac{p_{x}^{2} - p_{y}^{2}}{p_{t}^{2}} 〉, p_{t} = \sqrt{p_{x}^{2} + p_{y}^{2}} .

(11)

The reduced rapidity (y_z/y_pro) dependence of the directed (v₁) and elliptic (v₂) flows of free protons from ¹⁹⁷Au + ¹⁹⁷Au collisions at 40, 100, and 150 MeV/nucleon with different $F$ is shown in Fig.3. One can see that both v₁ and v₂ obtained with different $F$ are different, and the differences become more evident with increasing beam energy. The reason is that with increasing beam energy, the collision term plays a more important role. The v₁ and v₂ values obtained with the FU3FP1 and $F$ =0.3 are very close to each other. In addition, with a larger $F$ , which means a smaller reduction in the cross section, one can obtain a larger slope for v₁ and a more negative v₂ at mid-rapidity, y_z/y_pro = 0. Nucleons are known to be more likely to undergo a bounce-off (positive v₁ slope) motion and squeeze-out (negative v₂) pattern because of the increasing collision number.

Figure 3:

(Color online) Rapidity distribution of the directed flow v₁ (upper, b = 2–5.5 fm) and elliptic flow v₂ (lower, b = 5.5–7.5 fm) of free protons from Au + Au collisions at E_lab = 40 (left panels), 100 (middle panels), and 150 MeV/nucleon (right panels).

Figure 4 shows the v₁ slope and v₂ at y_z/y_pro = 0 for hydrogen isotopes calculated with different $F$ values in comparison with the FOPI and INDRA experimental data. The chosen impact parameters are b = 2–5.5 fm and b = 5.5–7.5 fm for v₁ and v₂, respectively. The v₁ slope is obtained assuming $v_{1} (y_{0}) = v_{11} y_{0} + v_{13} y_{0}^{3} + c$ in the range of |y₀|=|yz/y_pro| < 0.4. With increasing E_lab, the value of the v₁ slope increases, and the value of v₂ decreases. The results obtained with FU3FP1 and $F$ =0.3 are very close to each other. Furthermore, the correction factors $F$ = 0.2 and 0.5 are required to fit the experimental data at beam energies of 40 and 150 MeV/nucleon, respectively. In addition, calculations with the FU3FP1 parametrization can simultaneously reproduce both v₁ and v₂ well.

Figure 4:

(Color online) Flow excitation functions at mid-rapidity for hydrogen isotopes from ¹⁹⁷Au+¹⁹⁷Au collisions. The results with four in-medium correction factors are represented by different symbols and lines as indicated. The FOPI experimental data (open stars) and INDRA experimental data (solid stars) are taken from Ref. [58].

To understand more clearly the influence of the medium correction factor on the NNCS, Fig.5 shows the proton energy spectra of ¹⁹⁷Au+¹⁹⁷Au collisions at beam energies of 40 and 150 MeV/nucleon. Overall, the energy spectra drop exponentially with energy and can also be affected by changing the medium correction factor. The energy spectrum obtained with $F = σ_{N N}^{in - medium} / σ_{N N}^{free} = 0.2$ has a steeper slope than the others because the number of collisions is smallest in this case.

Figure 5:

(Color online) Kinetic energy distribution of free protons from ¹⁹⁷Au+¹⁹⁷Au central collisions at beam energies of 40 (left panels) and 150 (right panels) MeV/nucleon.

3.3 Nuclear stopping

The nuclear stopping power is another important observable that characterizes the transparency of colliding nuclei [15, 18, 19 43]. In this work, we calculated the two quantities RE and vartl for the same reaction. The former was proposed by the FOPI collaboration [59] and is defined as the ratio of the variances of the transverse and longitudinal rapidity distributions:

v a r t l = \frac{〈 y_{x}^{2} 〉}{〈 y_{z}^{2} 〉}, 〈 y_{x, z}^{2} 〉 = \frac{\sum (y_{x, z}^{2} N_{y_{x, z}})}{\sum N_{y_{x, z}}} .

(12)

The other quantity, R_E, was proposed by the INDRA collaboration and is defined as the ratio of the transverse and parallel energies:

R_{E} = \frac{\sum E_{⊥}}{2 \sum E_{∥}},

(13)

where E_⊥ (E_‖) is the transverse (parallel) kinetic energy of particles in the center-of-mass system [61].

The beam energy dependence of the degree of nuclear stopping for free protons in central HICs is displayed in Fig.6. Similar to the results shown in Fig.4, the degree of nuclear stopping can be reproduced well with $F$ ≈0.3 at E_lab = 40 MeV/nucleon. However, $F$ ≈0.5 is needed at E_lab = 150 MeV/nucleon. With increasing beam energy, a nuclear medium with higher density and momentum can be created; thus, it is reasonable that the reduction factors at different beam energies are different. With increasing beam energy, the difference in R_E between $F$ =0.3 and FU3FP1 increases steadily. The experimental data and the trend in stopping power with increasing beam energy can be reproduced reasonably well by the calculations with FU3FP1. In contrast, the other calculations remain almost constant because the reduced factor in the FU3FP1 depends on both the density and momentum. In addition, these conclusions can be confirmed in both ¹⁹⁷Au+¹⁹⁷Au and ¹²⁹Xe+¹²⁰Sn.

Figure 6:

(Color online) Comparison of R_E (solid symbols) and vartl (open symbols) for free protons from central ¹⁹⁷Au+¹⁹⁷Au (left panels) and ¹²⁹Xe+¹²⁰Sn (right panels) collisions with four different medium correction factors (dashed lines with different symbols) to the FOPI (open stars) and INDRA (solid stars) experimental data [35, 62].

4 Summary

The effects of the in-medium nucleon-nucleon cross section on the collective flow and nuclear stopping in the Fermi energy domain are investigated within the UrQMD model, in which medium correction factors of $F = σ_{N N}^{in - medium} / σ_{N N}^{free} = 0.2, 0.3, 0.5$ and the parametrization factor FU3FP1 (which is density- and momentum- dependent) are considered. It is demonstrated that to best fit both the collective flow and nuclear stopping data, $F$ ≈ 0.2 and 0.5 are required at beam energies of 40 and 150 MeV/nucleon, respectively. These results are consistent with the results obtained by the INDRA collaboration. The calculations with the FU3FP1 parametrization reproduce both the directed and elliptic flow data reasonably well. The slightly increased stopping power with increasing beam energy can also be found in the calculations with the FU3FP1 parametrization, whereas the calculations with $F$ = 0.2, 0.3, and 0.5 yield roughly flat energy-dependent.

References

[1]

J. Aichelin,

"Quantum" molecular dynamics-a dynamical microscopic n-body approach to investigate fragment formation and the nuclear equation of state in heavy ion collisions

. Phys. Rept. 202, 233 (1991). doi: 10.1016/0370-1573(91)90094-3