Probing neutron-proton effective mass splitting using nuclear stopping and isospin mix in heavy-ion collisions in GeV energy region

Probing neutron-proton effective mass splitting using nuclear stopping and isospin mix in heavy-ion collisions in GeV energy region 1 corresp first-author Corresponding author, zhangfan@mail.bnu.edu.cn Corresponding author, zhangfan@mail.bnu.edu.cn zhangfan@mail.bnu.edu.cn https://orcid.org/0000-0003-3642-6639 2 Department of Electronic Information and Physics, Changzhi University, Changzhi 046011, China Department of Electronic Information and Physics, Changzhi University, Changzhi 046011, China Sino-French Institute of Nuclear Engineering and Technology, Sun Yat-sen University, Zhuhai 519082, China Sino-French Institute of Nuclear Engineering and Technology, Sun Yat-sen University, Zhuhai 519082, China The ramifications of the effective mass splitting on the nuclear stopping and isospin tracer during heavy-ion collisions within the gigaelectron volt energy region are studied using an isospin-dependent quantum molecular dynamics model. Three isotope probes, i.e., a proton, deuteron, and triton, are used to calculate the nuclear stopping. Compared to the <math id="M1"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> case, the <math id="M2"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> parameter results in a stronger stopping for protons but a weaker stopping for tritons. The calculations of the isospin tracer show that the <math id="M3"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> parameter results in a higher isospin mix than the <math id="M4"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> parameter. The rapidity and impact parameter dependences of the isospin tracer are also studied. A constraining of the effective mass splitting using the free nucleons with high rapidity and in a central rather than peripheral collision is suggested. Neutron-proton effective mass splitting Nuclear stopping Isospin mix 1 INTRODUCTION 1 Despite wide investigations in both nuclear physics and nuclear astrophysics, the nuclear potential between nucleons is not fully understood [ 1 1 - 4 4 1-4 ]. There are three degrees of freedom for the nuclear potential, namely, the density ρ , momentum p, and isospin asymmetry δ = ( ρ n - ρ p )/( ρ n + ρ p ). For the density and isospin dependence, many theoretical and experimental studies have been conducted during the past decades [ 5 5 - 7 7 5-7 ]. Readers are encouraged to refer to previous reviews on the accepted acknowledge regarding two degrees of freedom [ 8 8 , 9 9 ]. The momentum dependence of the nuclear potential has been recently emphasized [ 10 10 - 12 12 10-12 ]. For example, it was found that the momentum degree of freedom is requisite in extracting the compressibility of nuclear matter [ 18 18 ]. The momentum dependence of the nuclear potential between a proton and nucleus has been measured and parameterized within an energy range of 10 MeV < E kin 1 GeV [ 13 13 - 17 17 13-17 ]. Apart from studying the momentum-dependent nuclear potential, the gradient of momentum, i.e., d U /d p , has also been studied. This is the concept of the effective mass m * [ 19 19 ]. An isoscalar effective mass of m * / m = 0.7 ± 0.05 at normal density and the Fermi momentum have been wildly used in transport models [ 20 20 ]. However, the concept of effective mass splitting is applied to describe the difference in the effective masses between a neutron and proton in asymmetric nuclear matter [ 21 21 - 24 24 21-24 ]. Knowledge regarding neutron-proton (np) effective mass splitting is essential for understanding neutron stars and exotic nuclei [ 25 25 , 26 26 ]. Unfortunately, np effective mass splitting is not fully understood [ 27 27 - 30 30 27-30 ]. Some theories have predicted <math id="M5"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> in neutron-rich matter [ 31 31 - 33 33 31-33 ], whereas other theories have made an opposite prediction [ 34 34 - 36 36 34-36 ]. Effective mass splitting plays an important role in the transport simulation of heavy-ion collisions [ 37 37 ]. There are two main parts in a transport simulation, i.e., the evolution in the mean field and nucleon-nucleon (NN) scattering. The effective mass splitting results in different mean fields of the proton and neutron. Thus, in an energy region of approximately 100 MeV/nucleon, the calculations of the isospin observable in heavy-ion collisions will depend on the effective mass splitting applied in the model. Using the double n/p spectral ratios and isospin diffusion, the role of the effective mass splitting on the mean field has been investigated [ 37 37 - 39 39 37-39 ]. However, few studies have reported the influence of the effective mass splitting on NN scattering, which is dominant in heavy-ion collisions at an energy level of gigaelectron volts [ 40 40 - 42 42 40-42 ]. It has been indicated that the effective mass provides an in-medium correction of the NN cross sections [ 43 43 , 44 44 ]. The in-medium factors related to the nucleon effective mass have been widely used in transport models [ 45 45 - 47 47 45-47 ]. Owing to the effective mass splitting, a splitting of the in-medium factors will occur between the neutron-neutron (nn) and proton-proton (pp) cross sections [ 48 48 ]. Thus, an observable isospin that is sensitive to the NN scattering may depend on the effective mass splitting. The effects of the in-medium NN cross sections on the observable stopping and the ratio of free protons have been found in 96 Zr( 96 Ru) + 96 Ru( 96 Zr) collisions at 400 MeV/nucleon [ 49 49 ]. In this study, an attempt was made to determine the effect of the effective mass splitting on the nuclear stopping and isospin mix in heavy-ion collisions at incident energies of 150 to 1,500 MeV/nucleon. The remainder of this paper is organized as follows. In Sect. 2, the proposed method is described. In Sect. 3, the results are presented and discussed. In Sect. 4, some concluding remarks are given. 2 THEORETICAL FRAMEWORK 1 In this study, a similar theoretical framework as that applied in a previous study is used [ 50 50 ]. Readers can refer to the references for the isospin-dependent quantum molecular dynamics (IQMD) model. For brevity, the nuclear potential energy and in-medium corrections to the NN cross sections are discussed here. The nuclear potential energy is given by the following: <math display="block" id="M6"><mtable columnalign="left"><mtr><mtd><mi>U</mi><mo>=</mo><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mrow><mo>[</mo><mrow><mfrac><mi>α</mi><mn>2</mn></mfrac><mfrac><mrow><msup><mi>ρ</mi><mn>2</mn></msup></mrow><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mi>β</mi><mrow><mi>γ</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mfrac><mrow><msup><mi>ρ</mi><mrow><mi>γ</mi><mo>+</mo><mn>1</mn></mrow></msup></mrow><mrow><msubsup><mi>ρ</mi><mn>0</mn><mi>γ</mi></msubsup></mrow></mfrac><mo>+</mo><mfrac><mrow><msub><mi>g</mi><mrow><mtext>sur</mtext></mrow></msub></mrow><mn>2</mn></mfrac><mfrac><mrow><msup><mrow><mo stretchy="false">(</mo><mo>∇</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>]</mo></mrow><mtext>d</mtext><mstyle mathvariant="bold" mathsize="normal"><mi>r</mi></mstyle></mtd></mtr><mtr><mtd><mo>+</mo><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mfrac><mrow><msub><mi>C</mi><mrow><mi>s</mi><mn>,</mn><mi>p</mi></mrow></msub></mrow><mn>2</mn></mfrac><msup><mrow><mo stretchy="false">(</mo><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><msub><mi>γ</mi><mi>i</mi></msub></mrow></msup><mi>ρ</mi><msup><mi>δ</mi><mn>2</mn></msup><mtext>d</mtext><mstyle mathvariant="bold" mathsize="normal"><mi>r</mi></mstyle></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mi>τ</mi></munder></mstyle><mi>C</mi><mo stretchy="false">(</mo><mn>1</mn><mo>+</mo><mi>x</mi><mo stretchy="false">)</mo><mstyle displaystyle="true"><mrow><mo>∭</mo></mrow></mstyle><mtext>d</mtext><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mtext>d</mtext><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mo>'</mo><mtext>d</mtext><mstyle mathvariant="bold" mathsize="normal"><mi>r</mi></mstyle><msub><mi>f</mi><mi>τ</mi></msub><mo stretchy="false">(</mo><mstyle mathvariant="bold" mathsize="normal"><mi>r</mi></mstyle><mn>,</mn><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>×</mo><msup><mstyle mathsize="140%" displaystyle="true"><mrow><mi>ln</mi></mrow></mstyle><mn>2</mn></msup><mo stretchy="false">[</mo><mi>ϵ</mi><msup><mrow><mo stretchy="false">(</mo><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mo>−</mo><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mo>'</mo><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>1</mn><mo stretchy="false">]</mo><msub><mi>f</mi><mi>τ</mi></msub><mo stretchy="false">(</mo><mstyle mathvariant="bold" mathsize="normal"><mi>r</mi></mstyle><mn>,</mn><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mo>'</mo><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>+</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mi>τ</mi></munder></mstyle><mi>C</mi><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mi>x</mi><mo stretchy="false">)</mo><mstyle displaystyle="true"><mrow><mo>∭</mo></mrow></mstyle><mtext>d</mtext><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mtext>d</mtext><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mo>'</mo><mtext>d</mtext><mstyle mathvariant="bold" mathsize="normal"><mi>r</mi></mstyle><msub><mi>f</mi><mi>τ</mi></msub><mo stretchy="false">(</mo><mstyle mathvariant="bold" mathsize="normal"><mi>r</mi></mstyle><mn>,</mn><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>×</mo><msup><mstyle mathsize="140%" displaystyle="true"><mrow><mi>ln</mi></mrow></mstyle><mn>2</mn></msup><mo stretchy="false">[</mo><mi>ϵ</mi><msup><mrow><mo stretchy="false">(</mo><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mo>−</mo><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mo>'</mo><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>1</mn><mo stretchy="false">]</mo><msub><mi>f</mi><mrow><mo>−</mo><mi>τ</mi></mrow></msub><mo stretchy="false">(</mo><mstyle mathvariant="bold" mathsize="normal"><mi>r</mi></mstyle><mn>,</mn><mstyle mathvariant="bold" mathsize="normal"><mi>p</mi></mstyle><mo>'</mo><mo stretchy="false">)</mo><mn>.</mn></mtd></mtr></mtable></math> Here, ρ is the nucleon density, δ = ( ρ n - ρ p )/( ρ n + ρ p ) is the asymmetry denoted by the neutron density ρ n and proton density ρ p , and f τ (r, p) is the phase-space density, with τ = -1/2 for a proton and τ = 1/2 for a neutron. The summation includes both a neutron and proton. The parameters α = -130 MeV, β = 59 MeV, γ = 2.09, g sur = 130 MeV, C = 1.57 MeV, and ϵ = 0.0005 c 2 /MeV 2 provide a hard equation of state ( K = 380 MeV) for symmetric nuclear matter [ 51 51 ]. In fact, the commonly accepted value of the compressibility is approximately K = 240 MeV. In [ 17 17 ], it has been shown that the calculation of the nuclear stopping using a tiff compressibility of K = 380 MeV with a momentum-dependent interaction agrees more with the data than that for a soft case ( K = 200 MeV) with a momentum dependent interaction. The effective mass splitting rather than the compressibility is focused on, and thus a compressibility of K = 380 MeV is still used. The effective mass splitting depends on parameter x . The value x = -0.65 is applied to provide <math id="M7"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> in neutron rich nuclear matter, and x = 0.65 for <math id="M8"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> . The parameters C s , p = 23.5 MeV and γ i = 1.0 are used for x = 0.65, whereas C s , p = 52.5 MeV and γ i = 1.0 for x = -0.65 to give the same density dependence of the symmetry energy. These parameters have been applied to calculate the nuclear collective flows in the HICs [ 35 35 , 36 36 ]. Refer to Fig. 1 Fig. 1 in [ 35 35 ] for the symmetry energy as a function of density, and Fig. 1 Fig. 1 in [ 36 36 ] for the lane potential as a function of momentum provided by the parameters. Particularly in asymmetric nuclear matter with δ = 0.2 at the saturation density, <math id="M9"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> results in an effective mass of <math id="M10"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>/</mo><mi>m</mi></mrow></math> = 0.72 and <math id="M11"><mrow><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup><mo>/</mo><mi>m</mi></mrow></math> = 0.68, whereas <math id="M12"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> results in effective masses of <math id="M13"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>/</mo><mi>m</mi></mrow></math> = 0.68 and <math id="M14"><mrow><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup><mo>/</mo><mi>m</mi></mrow></math> = 0.72. It is worth noting that the integrals of momenta in Eq.(1) are solved by a summation in the code. Meanwhile, the phase-space densities f τ are reduced to the nucleon densities, which only depend on the coordinate r. In this case, the difference in the algorithms between the momentum dependence and two body interactions is only the intensities of the interaction between two nucleons. The former intensity depends on the momenta of two nucleons, whereas the latter intensity is <math id="M15"><mrow><mfrac><mi>α</mi><mrow><mn>2</mn><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac></mrow></math> . In-medium corrections to NN scattering cross sections are needed in the calculations. An isospin independent parametrization, expressed as <math id="M16"><mrow><msubsup><mi>f</mi><mrow><mtext>NN</mtext></mrow><mrow><mtext>med</mtext></mrow></msubsup><mo>=</mo><mn>1</mn><mo>−</mo><mn>0.2</mn><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac></mrow></math> , has been used in [ 52 52 ]. In this study, the effective mass of each nucleon is calculated from the momentum-dependent potential. Next, the in-medium factor of an NN cross section for i th and j th nucleons is calculated [ 43 43 - 46 46 43-46 ] <math display="block" id="M17"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><msubsup><mi>f</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mtext>med</mtext></mrow></msubsup><mo>=</mo><mfrac><mrow><msub><mi>σ</mi><mrow><mtext>med</mtext></mrow></msub></mrow><mrow><msub><mi>σ</mi><mrow><mtext>free</mtext></mrow></msub></mrow></mfrac><mo>=</mo><msup><mrow><mrow><mo>[</mo><mrow><mfrac><mrow><msubsup><mi>m</mi><mi>i</mi><mo>*</mo></msubsup><msubsup><mi>m</mi><mi>j</mi><mo>*</mo></msubsup><mo>/</mo><mo stretchy="false">(</mo><msubsup><mi>m</mi><mi>i</mi><mo>*</mo></msubsup><mo>+</mo><msubsup><mi>m</mi><mi>j</mi><mo>*</mo></msubsup><mo stretchy="false">)</mo></mrow><mrow><msub><mi>m</mi><mi>i</mi></msub><msub><mi>m</mi><mi>j</mi></msub><mo>/</mo><mo stretchy="false">(</mo><msub><mi>m</mi><mi>i</mi></msub><mo>+</mo><msub><mi>m</mi><mi>j</mi></msub><mo stretchy="false">)</mo></mrow></mfrac></mrow><mo>]</mo></mrow></mrow><mn>2</mn></msup><mn>,</mn></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><msubsup><mi>m</mi><mi>i</mi><mo>*</mo></msubsup><mo>=</mo><msup><mrow><mrow><mo>[</mo><mrow><mfrac><mn>1</mn><mrow><msub><mi>m</mi><mi>i</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>∂</mo><mi>U</mi></mrow><mrow><msub><mi>p</mi><mi>i</mi></msub><mo>∂</mo><msub><mi>p</mi><mi>i</mi></msub></mrow></mfrac></mrow><mo>]</mo></mrow></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><mn>.</mn></mrow></mtd></mtr></mtable></mrow></math> The in-medium factor <math id="M18"><mrow><msubsup><mi>f</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mtext>med</mtext></mrow></msubsup></mrow></math> is an isospin and is energy dependent [see Eq. (1)]. Nucleon colliding pairs with E lab = 400 MeV are selected from all simulated events at t = 20 fm/ c and then applied to calculate the in-medium correction factors to show the isospin dependence. The results for the central reaction of 96 Zr + 96 Zr at 400 MeV/nucleon are shown in Fig. 1 Fig. 1 . Two panels correspond to the effective mass splitting of <math id="M19"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> and <math id="M20"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> . The dash-dotted curve is the isospin independent parametrization <math id="M21"><mrow><mn>1</mn><mo>−</mo><mn>0.2</mn><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac></mrow></math> . The other three curves are the in-medium factors of nn, np, and pp cross sections extracted from the effective masses. The values of both cases are shown to be similar at subnormal and normal densities. However, in the density region above 1.5 ρ 0 , the in-medium factors extracted from the effective masses are globally larger than in the parameterized case. It can be seen that the free-space cross sections of the nn and pp scattering are the same. However, the in-medium nn cross sections are smaller than those for the pp scattering when using the parameters providing <math id="M22"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> . Oppositely, <math id="M23"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> leads to larger in-medium nn cross sections than the pp case. Finally, it should be noted that the Pauli blocking algorithm is different compared to that in [ 50 50 ]. The phase-space density <math id="M24"><mrow><msub><msup><mi>f</mi><mo>′</mo></msup><mi>i</mi></msub></mrow></math> is calculated after the scattering. It is then interpreted as a blocking probability. The probability of (1 - <math id="M25"><mrow><msub><msup><mi>f</mi><mo>′</mo></msup><mn>1</mn></msub></mrow></math> )(1 - <math id="M26"><mrow><msub><msup><mi>f</mi><mo>′</mo></msup><mn>2</mn></msub></mrow></math> ) is applied to allow the scattering. 3 RESULTS AND DISCUSSIONS 1 It was determined that nuclear stopping can be influenced by the in-medium effect of NN scattering [ 49 49 ]. The in-medium factor of NN scattering has been calculated from the effective mass splitting in this study [see Eq. (2)] . Thus, the calculation of nuclear stopping is expected to be affected by the effective mass splitting. The stopping of observable varxz calculated from the proton, deuteron, and triton rapidity is applied. The definition of the observable varxz and the data are referenced from [ 40 40 ], <math display="block" id="M27"><mrow><mi>v</mi><mi>a</mi><mi>r</mi><mi>x</mi><mi>z</mi><mo>=</mo><mfrac><mrow><mo>〈</mo><msubsup><mi>y</mi><mi>x</mi><mn>2</mn></msubsup><mo>〉</mo></mrow><mrow><mo>〈</mo><msubsup><mi>y</mi><mi>z</mi><mn>2</mn></msubsup><mo>〉</mo></mrow></mfrac></mrow></math> where y x and y z are the transverse and longitudinal rapidity, respectively. The averages are obtained for the proton, deuteron, and triton. Compared to global nuclear stopping, the stopping observable for the isotopes can provide further insight into the isospin property. As an example, Fig. 2 Fig. 2 shows the excitation function of the stopping of observable varxz for the proton, deuteron, and triton produced in central Au + Au collisions ( b < 2 fm). The values of varxz for three isotopes show a similar rise and fall with the maximum at near 400 MeV/nucleon, which has also been found for the global stopping [ 41 41 ]. The maximum nuclear stopping, which is related to the competition between the mean field, Pauli blocking, and inelastic collisions, is not perfectly reproduced by the model. The phenomenological method for the Pauli blocking may be responsible for this. However, the isotope dependence of the observable varxz , which is guided by the dotted line, is reproduced well by the model. For the same incident energy, more stopping occurs for a proton, less occurs for a triton, and the average amount occurs for a deuteron. It should be noted that there are two data points deviating from the systematic nature of isotope dependence (shown in the dashed circle). The model does not reproduce such abnormal increases. In Fig. 3 Fig. 3 , the stopping of observable varxz is shown as a function of the system mass. The systems 40 Ca + 40 Ca, 58 Ni + 58 Ni, 96 Ru + 96 Ru, 129 Xe + 131 I, and 197 Au + 197 Au, and incident energies of 150 to 1500 MeV/nucleon, are considered. For the FOPI data, the central events are sorted based on the ratio of total kinetic energy between the transverse and longitudinal directions. As indicated by the FOPI collaboration, the dependence of the system size on the nuclear stopping is observable for studying the transparency of the nuclei. The data support the transparency of the nuclei and incident energies involved. All three isotope probes show that lighter nuclei have higher transparency. Moreover, the dependence of the system size is more obvious in a higher incident energy region (see Fig. 24 in [ 40 40 ]). Two aspects may be responsible for the system size dependence of nuclear stopping. First, a peripheral collision has more transparency than a central collision. The central collision is defined by b < 0.15 b max , where the maximum b max = 1.15( <math id="M30"><mrow><msubsup><mtext>A</mtext><mtext>p</mtext><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup><mo>+</mo><msubsup><mtext>A</mtext><mtext>t</mtext><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msubsup></mrow></math> ) fm is approximate. The system mass and incident energy dependencies of b max may influence the centralities of the sorted events. Second, the surfaces of the nuclei are more transparent than the core with normal density. For light nuclei, the proportion of the surface is larger than the core case. Note the fact that the mean field, i.e., the long-range effect, dominates the nuclear stopping at a low incident energy. The long-range effect allows the core to stop more of the surface. This may be responsible for the weak system size dependence of the nuclear stopping at low energy. Regarding the effect of the effective mass splitting, note the global difference between the solid and dashed curves in Fig. 2 Fig. 2 . Owing to more pp scatterings provided by the <math id="M34"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> case than the <math id="M35"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> case (referring to Fig. 1 Fig. 1 ), the <math id="M36"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> results in a stronger stopping for the proton, as shown in Fig. 2(a) Fig. 2(a) . Through a similar inference to nn scatterings, a stronger stopping for neutron and neutron-rich fragments is expected when <math id="M37"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> is considered in the calculation. Similar effects can be observed for all reaction systems in Fig. 3 Fig. 3 . For the proton probe, as shown in Fig. 3(a) Fig. 3(a) , the squares representing the <math id="M38"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> case are globally at the top of the triangles representing the <math id="M39"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> case. The effect globally for the triton probe can be seen in Fig. 3(c) Fig. 3(c) . However, the effect is too weak compared to the error bars in the data. To show this effect more clearly, the observable differences in the proton and triton are shown in Fig. 4 Fig. 4 . Once again, the squares representing the <math id="M40"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> case are globally at the top of the triangles representing the <math id="M41"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> case. However, it cannot be concluded which case reproduces the data better because the errors in the data are too large compared to the effect of the effective mass splitting. Another observation related to the amount of rapidity is proposed in [ 42 42 ]. The relative abundance of the free protons was adopted to trace the isospin mixing between the projectile and target nucleons <math display="block" id="M42"><mrow><msub><mi>R</mi><mi>Z</mi></msub><mo>=</mo><mfrac><mrow><mn>2</mn><msup><mi>Z</mi><mrow><mtext>mix</mtext></mrow></msup><mo>−</mo><msup><mi>Z</mi><mrow><mtext>Zr</mtext></mrow></msup><mo>−</mo><msup><mi>Z</mi><mrow><mtext>Ru</mtext></mrow></msup></mrow><mrow><msup><mi>Z</mi><mrow><mtext>Zr</mtext></mrow></msup><mo>−</mo><msup><mi>Z</mi><mrow><mtext>Ru</mtext></mrow></msup></mrow></mfrac><mn>,</mn></mrow></math> where Z mix , Z Zr , and Z Ru are the number of free protons in reactions of 96 Ru + 96 Zr (or 96 Zr + 96 Ru), 96 Zr + 96 Zr, and 96 Ru + 96 Ru, respectively. In the case of full equilibrium in the mixed reactions 96 Ru + 96 Zr and 96 Zr + 96 Ru, the value of R Z is zero. The centrality dependence of R Z for both 96 Ru + 96 Zr and 96 Zr + 96 Ru reactions at 400 MeV/nucleon has been determined using FOPI collaboration, for which a rapidity range of -0.75 < y 0 < -0.25 is considered [ 42 42 ]. The proton and t/ 3 He are applied as the probes, as shown in Figs. 5 Figs. 5 and 6 6 . It can be seen that the non-zero R Z for 96 Ru + 96 Zr and 96 Zr + 96 Ru appear to be symmetric relative to the axis R Z = 0, supporting an image of incomplete isospin equilibration. Moreover, the difference between two mixed reactions increases slightly with the impact parameter. Compared to those in the peripheral collisions, the temperatures and densities are higher in a central collision. Hence, there is less mixing in the peripheral collisions and more mixing in the central collisions. Using the IQMD model, the values of R Z are calculated using the free protons in the rapidity ranges of -0.75 < y 0 < -0.25 (shown in Fig. 5(a) Fig. 5(a) ) and -1.25 < y 0 < -0.75 (shown in Fig. 5(b) Fig. 5(b) ), respectively. The center of each square or triangle is the average of those five simulated groups, whereas the error bar refers to the standard deviation. In Figs. 5(a) Figs. 5(a) and 6(a) 6(a) , the effect of the effective mass splitting can be seen. However, the effective mass splitting cannot be constrained because the calculations for <math id="M43"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> and <math id="M44"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> are similar and generally agree with the data. Interestingly, most values for <math id="M45"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> are closer to zero than the <math id="M46"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> case, except for b = 4.0 fm. The divergences in R Z between the two calculations are more obvious for a larger rapidity (-1.25 < y 0 < -0.75) and at a smaller centrality. This indicates that the observable R Z for fast protons in the central collisions is sensitive to the effective mass splitting. The rapidity dependence of R Z for central 96 Ru + 96 Zr reactions at 400 MeV/nucleon was determined using FOPI collaboration [ 42 42 ], as shown in the circles in Fig. 7 Fig. 7 . The data for negative rapidity are available. Because of the reaction symmetry, the data are shown for 96 Zr + 96 Ru in a positive rapidity region. The data show approximately zero values of R Z at midrapidity ( y 0 0), with positive values on the y 0 < 0 side and negative values on the y 0 > 0 side. The positive values of R Z indicate that the protons in this rapidity cell are predominantly populated from the Zr target, whereas the negative values indicate they are populated from a Ru projectile [ 42 42 ]. To investigate the sensitivity of R Z to the effective mass splitting, simulations with <math id="M49"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> and <math id="M50"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> are conducted, respectively. Here, the number of events of each group is 300,000. In fact, the symmetry energy is also expected to affect the isospin mixings, but there is no difference between when a soft or stiff symmetry energy is applied to the calculation. In fact, during the compression phase, with the nucleon density approaching approximately 2 ρ 0 , soft parameters provide a stronger repulsion of protons than stiff parameters because the colliding system is neutron-rich. However, in a stiff central collision at 400 MeV/nucleon, the NN scatterings are violent. The similar values of R Z for soft and stiff symmetry energies support the hypothesis that the effect of the local symmetry potential to an isospin tracer is reduced by the violent NN scatterings. The calculations for <math id="M51"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> and <math id="M52"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> are shown as squares and triangles in Fig. 7 Fig. 7 , respectively. It can be seen that the results for <math id="M53"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> show more isospin stoppings, and even an over balance, whereas <math id="M54"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> is more appropriate for reproducing the data. Our results in Fig. 1 Fig. 1 , together with the calculations in [ 45 45 , 46 46 ], have shown that <math id="M55"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> results in a larger in-medium factor under more scattering than <math id="M56"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> . This could be the dynamic reason for more proton mixing for <math id="M57"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> than for <math id="M58"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> because violent NN scattering, rather than the mean field, is responsible for the isospin mixing. Figure 7 Figure 7 shows the robust effect of the effective mass splitting on the isospin tracer for the central collision and high rapidity region. However, other effects on the isospin tracer, such as those of the compressibility and Pauli blocking, should be investigated before using the isospin tracer data to constrain the effective mass splitting. 4 CONCLUSION 1 In summary, the effects of the effective mass splitting on the nuclear stopping during heavy-ion collisions within a gigaelectron volt energy region were studied within the IQMD model. Three isotope probes, i.e., proton, deuteron, and triton, were used to calculate the nuclear stopping. Compared to the <math id="M59"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> case, the <math id="M60"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> parameter results in a stronger stopping for a proton but a weaker stopping for triton. The opposite influences of the effective mass splitting on the proton-proton and neutron-neutron scatterings may be responsible for the above phenomenon. However, it is worth mentioning that, not only the number, but also the angular distribution of the scattering, influence the nuclear stopping. The effect of the effective mass splitting on the angular distribution was further studied. The effect of the effective mass splitting on the isospin tracer was also considered. The calculations of the isospin tracer show that the <math id="M61"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo>></mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> parameter results in more isospin mixing than the case of <math id="M62"><mrow><msubsup><mi>m</mi><mtext>n</mtext><mo>*</mo></msubsup><mo><</mo><msubsup><mi>m</mi><mtext>p</mtext><mo>*</mo></msubsup></mrow></math> . In fact, the effects of the effective mass splitting on the nuclear stopping and isospin tracer were investigated in [ 49 49 ], using the central 96 Ru + 96 Zr collision at 400 MeV/nucleon. Compared to [ 49 49 ], the innovations in this study are as follows. The effect of the effective mass splitting on the nuclear stopping is ubiquitous in 40 Ca + 40 Ca, 58 Ni + 58 Ni, 96 Ru + 96 Ru, 129 Xe + 131 I, and 197 Au + 197 Au at incident energies of 150, 250, 400, 1,000, and 1,500 MeV/nucleon. The system mass dependence shows that the effect is stronger in a lighter system. However, the effect is quite weak compared to the error in the data, although the differences in the varxz values between a proton and triton are used to intensify the effect. Concerning the isospin mixing, the impact parameter and the dependence on the rapidity were studied using a proton and the ratio between 3 H and 3 He as the probes. It was found that the effect of the effective mass splitting on an isospin tracer is the strongest for a central collision and in a high rapidity region. However, other effects on the isospin tracer, such as those of the compressibility and Pauli blocking, should be investigated before using the isospin tracer data to constrain the effective mass splitting. References H. A. Bethe , Supernova mechanisms . Rev. Mod. Phys . 62 , 801 ( 1990 ). https://doi.org/10.1103/RevModPhys.62.801 https://doi.org/10.1103/RevModPhys.62.801 Supernova mechanisms A.W. Steiner , M. Prakash , J.M. Lattimer et al ., Isospin asymmetry in nuclei and neutron stars . Phys. Rep . 411 , 325 ( 2005 ). https://doi.org/10.1016/j.physrep.2005.02.004 https://doi.org/10.1016/j.physrep.2005.02.004 Isospin asymmetry in nuclei and neutron stars P. Danielewicz , R. Lacey , W.G. Lynch , Determination of the equation of state of dense matter . Science 298 , 1592 ( 2002 ). https://doi.org/10.1126/science.1078070 https://doi.org/10.1126/science.1078070 Determination of the equation of state of dense matter H.T. Janka , Explosion mechanisms of core-collapse supernovae . Annu. Rev. Nucl. Part. Sci . 62 , 407 ( 2012 ). https://doi.org/10.1146/annurev-nucl-102711-094901 https://doi.org/10.1146/annurev-nucl-102711-094901 Explosion mechanisms of core-collapse supernovae T.Z. Yan , S. Li , Y.N. Wang et al ., Yield ratios and directed flows of light particles from proton-rich nuclei-induced collisions . Nucl. Sci. Tech . 30 , 15 ( 2019 ). https://doi.org/10.1007/s41365-018-0534-6 https://doi.org/10.1007/s41365-018-0534-6 Yield ratios and directed flows of light particles from proton-rich nuclei-induced collisions Byungsik Hong , Jung Keun Ahn , Gyeonghwan Bak et al ., Development of large acceptance multi-purpose spectrometer in Korea for symmetry energy . Nucl. Sci. Tech . 29 , 171 ( 2018 ). https://doi.org/10.1007/s41365-018-0507-9 https://doi.org/10.1007/s41365-018-0507-9 Development of large acceptance multi-purpose spectrometer in Korea for symmetry energy B.-J. Cai , L.-W. Chen , Constraints on the skewness coefficient of symmetric nuclear matter within the nonlinear relativistic mean field model . Nucl. Sci. Tech . 28 , 185 ( 2017 ). https://doi.org/10.1007/s41365-017-0329-1 https://doi.org/10.1007/s41365-017-0329-1 Constraints on the skewness coefficient of symmetric nuclear matter within the nonlinear relativistic mean field model C. Fuchs and H.H. Wolter , Modelization of the EOS . Eur. Phys. J. A 30 , 5 ( 2006 ). https://doi.org/10.1140/epja/i2005-10313-x https://doi.org/10.1140/epja/i2005-10313-x Modelization of the EOS B. A. Li , L. W. Chen , and C. M. Ko , Recent progress and new challenges in isospin physics with heavy-ion reactions . Phys. Rep . 464 , 113 ( 2008 ). https://doi.org/10.1016/j.physrep.2008.04.005 https://doi.org/10.1016/j.physrep.2008.04.005 Recent progress and new challenges in isospin physics with heavy-ion reactions V. Baran , M. Colonna , V. Greco et al ., Reaction dynamics with exotic nuclei . Phys. Rep . 410 , 335 ( 2005 ). https://doi.org/10.1016/j.physrep.2004.12.004 https://doi.org/10.1016/j.physrep.2004.12.004 . Reaction dynamics with exotic nuclei L.L. Li , Z.H. Li , E.G. Zhao , et al ., Isospin splitting of the nucleon-nucleus optical potential . Phys. Rev. C 80 , 064607 ( 2009 ). https://doi.org/10.1103/PhysRevC.80.064607 https://doi.org/10.1103/PhysRevC.80.064607 Isospin splitting of the nucleon-nucleus optical potential Q. B. Shen , Y. L. Han , and H. R Guo , Isospin dependent nucleon-nucleus optical potential with Skyrme interactions . Phys. Rev. C 80 , 024604 ( 2009 ). https://doi.org/10.1103/PhysRevC.80.024604 https://doi.org/10.1103/PhysRevC.80.024604 Isospin dependent nucleon-nucleus optical potential with Skyrme interactions J.P. Jeuhenne , A. Lejeune , and C. Mahaux , Many-body theory of nuclear matter . Phys. Rep . 25 , 85 ( 1976 ). https://doi.org/10.1016/0370-1573(76)90017-X https://doi.org/10.1016/0370-1573(76)90017-X Many-body theory of nuclear matter L.G. Arnold , B.C. Clark , E.D. Cooper , et al ., Energy dependence of the p - 40Ca optical potential: A Dirac equation perspective . Phys. Rev. C 25 , 936 ( 1982 ). https://doi.org/10.1103/PhysRevC.25.936 https://doi.org/10.1103/PhysRevC.25.936 Energy dependence of the p - 40Ca optical potential: A Dirac equation perspective S. Hama , B.C. Clark , E.D. Cooper et al ., Global Dirac optical potentials for elastic proton scattering from heavy nuclei . Phys. Rev. C 41 , 2737 ( 1990 ). https://doi.org/10.1103/PhysRevC.41.2737 https://doi.org/10.1103/PhysRevC.41.2737 Global Dirac optical potentials for elastic proton scattering from heavy nuclei C. Hartnack and J. Aichelin , New parametrization of the optical potential . Phys. Rev. C 49 , 2801 ( 1994 ). https://doi.org/10.1103/PhysRevC.49.2801 https://doi.org/10.1103/PhysRevC.49.2801 New parametrization of the optical potential G. Q. Zhang , Y. G. Ma , X. G. Cao et al ., Unified description of nuclear stopping in central heavy-ion collisions from 10 MeV to 1.2 GeV . Phys. Rev. C 84 , 034612 ( 2011 ). https://doi.org/10.1103/PhysRevC.84.034612 https://doi.org/10.1103/PhysRevC.84.034612 Unified description of nuclear stopping in central heavy-ion collisions from 10 MeV to 1.2 GeV J. Aichelin , A. Rosenhauer , G. Peilert et al ., Importance of momentum-dependent interactions for the extraction of the nuclear equation of state from high-energy heavy-ion collisions . Phys. Rev. Lett . 58 , 1926 ( 1987 ). https://doi.org/10.1103/PhysRevLett.58.1926 https://doi.org/10.1103/PhysRevLett.58.1926 Importance of momentum-dependent interactions for the extraction of the nuclear equation of state from high-energy heavy-ion collisions E.N.E. van Dalen , C. Fuchs , and A. Faessler , Effective nucleon masses in symmetric and asymmetric nuclear matter . Phys. Rev. Lett . 95 , 022302 ( 2005 ). https://doi.org/10.1103/PhysRevLett.95.022302 https://doi.org/10.1103/PhysRevLett.95.022302 Effective nucleon masses in symmetric and asymmetric nuclear matter M. Jaminon and C. Mahaux , Effective masses in relativistic approaches to the nucleon-nucleus mean field . Phys.Rev.C 40 , 354 ( 1989 ). https://doi.org/10.1103/PhysRevC.40.354 https://doi.org/10.1103/PhysRevC.40.354 Effective masses in relativistic approaches to the nucleon-nucleus mean field H.Y. Kong , J. Xu , L.W. Chen et al ., Constraining simultaneously nuclear symmetry energy and neutron-proton effective mass splitting with nucleus giant resonances using a dynamical approach . Phys. Rev. C 95 , 034324 ( 2017 ). https://doi.org/10.1103/PhysRevC.95.034324 https://doi.org/10.1103/PhysRevC.95.034324 Constraining simultaneously nuclear symmetry energy and neutron-proton effective mass splitting with nucleus giant resonances using a dynamical approach B. A. Li , X. Han , What can we learn about the neutron-proton effective mass splitting from constraints on the density dependence of nuclear symmetry energy around normal density . Phys. Lett. B 727 , 276 ( 2013 ). https://doi.org/10.1016/j.physletb.2013.10.006 https://doi.org/10.1016/j.physletb.2013.10.006 What can we learn about the neutron-proton effective mass splitting from constraints on the density dependence of nuclear symmetry energy around normal density C. B. Das , S.D. Gupta , C. Gale et al ., Momentum dependence of symmetry potential in asymmetric nuclear matter for transport model calculations . Phys. Rev. C 67 , 034611 ( 2003 ). https://doi.org/10.1103/PhysRevC.67.034611 https://doi.org/10.1103/PhysRevC.67.034611 Momentum dependence of symmetry potential in asymmetric nuclear matter for transport model calculations Ch.C. Moustakidis , Temperature and momentum dependence of single-particle properties in hot asymmetric nuclear matter . Phys. Rev. C 78 , 054323 ( 2008 ). https://doi.org/10.1103/PhysRevC.78.054323 https://doi.org/10.1103/PhysRevC.78.054323 Temperature and momentum dependence of single-particle properties in hot asymmetric nuclear matter J. W. Negele and K. Yazaki , Mean free path in a nucleus . Phys. Rev. Lett . 62 , 71 ( 1981 ). https://doi.org/10.1103/PhysRevLett.47.71 https://doi.org/10.1103/PhysRevLett.47.71 Mean free path in a nucleus B. A. Li , Constraining the neutron-proton effective mass splitting in neutron-rich matter . Phys. Rev. C 69 , 064602 ( 2004 ). https://doi.org/10.1103/PhysRevC.69.064602 https://doi.org/10.1103/PhysRevC.69.064602 Constraining the neutron-proton effective mass splitting in neutron-rich matter X.H. Li , W.J. Guo , B.A. Li , et al ., Neutron-proton effective mass splitting in neutron-rich matter at normal density from analyzing nucleon Cnucleus scattering data within an isospin dependent optical model . Phys. Lett. B 743 , 408 ( 2015 ). https://doi.org/10.1016/j.physletb.2015.03.005 https://doi.org/10.1016/j.physletb.2015.03.005 Neutron-proton effective mass splitting in neutron-rich matter at normal density from analyzing nucleon Cnucleus scattering data within an isospin dependent optical model Z. Q. Feng , Effective mass splitting of neutron and proton and isospin emission in heavy-ion collisions . Nucl. Phys. A 878 , 3 ( 2012 ). https://doi.org/10.1016/j.nuclphysa.2012.01.014 https://doi.org/10.1016/j.nuclphysa.2012.01.014 Effective mass splitting of neutron and proton and isospin emission in heavy-ion collisions C. Xu , B.A. Li , and L. W. Chen , Symmetry energy, its density slope, and neutron-proton effective mass splitting at normal density extracted from global nucleon optical potentials . Phys. Rev. C 82 , 054607 ( 2010 ). https://doi.org/10.1103/PhysRevC.82.054607 https://doi.org/10.1103/PhysRevC.82.054607 Symmetry energy, its density slope, and neutron-proton effective mass splitting at normal density extracted from global nucleon optical potentials Z. Q. Feng , Momentum dependence of the symmetry potential and its influence on nuclear reactions . Phys. Rev. C 84 , 024610 ( 2011 ). https://doi.org/10.1103/PhysRevC.84.024610 https://doi.org/10.1103/PhysRevC.84.024610 Momentum dependence of the symmetry potential and its influence on nuclear reactions E.N.E. van Dalen , C. Fuchs , and A. Faessler , Momentum, density, and isospin dependence of symmetric and asymmetric nuclear matter properties . Phys. Rev. C 72 , 065803 ( 2005 ). https://doi.org/10.1103/PhysRevC.72.065803 https://doi.org/10.1103/PhysRevC.72.065803 Momentum, density, and isospin dependence of symmetric and asymmetric nuclear matter properties W. Zuo , L.G. Cao , B.A. Li et al ., Isospin splitting of the nucleon mean field . Phys. Rev. C 60 , 024605 ( 2005 ). https://doi.org/10.1103/PhysRevC.72.014005 https://doi.org/10.1103/PhysRevC.72.014005 Isospin splitting of the nucleon mean field L. Ou , Z. Li , Y. Zhang et al ., Effect of the splitting of the neutron and proton effective masses on the nuclear symmetry energy at finite temperatures . Phys. Lett. B 697 , 246 ( 2011 ). https://doi.org/10.1016/j.physletb.2011.01.062 https://doi.org/10.1016/j.physletb.2011.01.062 Effect of the splitting of the neutron and proton effective masses on the nuclear symmetry energy at finite temperatures J. Rizzo , M. Colonna , M. Di Toro et al ., Transport properties of isospin effective mass splitting . Nucl. Phys. A 732 , 202 ( 2004 ). https://doi.org/10.1016/j.nuclphysa.2003.11.057 https://doi.org/10.1016/j.nuclphysa.2003.11.057 Transport properties of isospin effective mass splitting Z. Q. Feng , Transverse emission of isospin ratios as a probe of high-density symmetry energy in isotopic nuclear reactions . Phys. Lett. B 707 , 83 ( 2012 ). https://doi.org/10.1016/j.physletb.2011.12.001 https://doi.org/10.1016/j.physletb.2011.12.001 Transverse emission of isospin ratios as a probe of high-density symmetry energy in isotopic nuclear reactions W. J. Xie and F. S. Zhang , Nuclear collective flows as a probe to the neutron Cproton effective mass splitting . Phys. Lett. B 735 , 250 ( 2014 ). https://doi.org/10.1016/j.physletb.2014.06.050 https://doi.org/10.1016/j.physletb.2014.06.050 Nuclear collective flows as a probe to the neutron Cproton effective mass splitting Y.X. Zhang , M.B. Tsang , Z. Li et al ., Constraints on nucleon effective mass splitting with heavy ion collisions . Phys. Lett. B 732 , 186 ( 2014 ). https://doi.org/10.1016/j.physletb.2014.03.030 https://doi.org/10.1016/j.physletb.2014.03.030 Constraints on nucleon effective mass splitting with heavy ion collisions L. W. Chen , C. M. Ko , and B. A Li , Determination of the stiffness of the nuclear symmetry energy from isospin diffusion . Phys. Rev. Lett . 94 , 032701 ( 2005 ). https://doi.org/10.1103/PhysRevLett.94.032701 https://doi.org/10.1103/PhysRevLett.94.032701 Determination of the stiffness of the nuclear symmetry energy from isospin diffusion H.Y. Kong , Y. Xia , J. Xu et al ., Reexamination of the neutron-to-proton-ratio puzzle in intermediate-energy heavy-ion collisions . Phys. Rev. C 91 , 047601 ( 2015 ). https://doi.org/10.1103/PhysRevC.91.047601 https://doi.org/10.1103/PhysRevC.91.047601 Reexamination of the neutron-to-proton-ratio puzzle in intermediate-energy heavy-ion collisions W. Reisdorf , A. Andronic , R. Averbeck et al ., (FOPI Collaboration), Systematics of central heavy ion collisions in the 1A GeV regime . Nucl. Phys. A 848 , 366 ( 2010 ). https://doi.org/10.1016/j.nuclphysa.2010.09.008 https://doi.org/10.1016/j.nuclphysa.2010.09.008 (FOPI Collaboration), Systematics of central heavy ion collisions in the 1A GeV regime A. Andronic , J. Łukasik1 , W. Reisdorf et al ., Systematics of stopping and flow in Au+ Au collisions . Eur. Phys. J. A 30 , 31 ( 2006 ). https://doi.org/10.1140/epja/i2006-10101-2 https://doi.org/10.1140/epja/i2006-10101-2 Systematics of stopping and flow in Au+ Au collisions F. Rami , Y. Leifels , B. De Schauenburg et al ., (FOPI Collaboration), Isospin tracing: A probe of nonequilibrium in central heavy-ion collisions . Phys. Rev. Lett . 84 , 1120 ( 2000 ). https://doi.org/10.1103/PhysRevLett.84.1120 https://doi.org/10.1103/PhysRevLett.84.1120 (FOPI Collaboration), Isospin tracing: A probe of nonequilibrium in central heavy-ion collisions D. Persram and C. Gale , Elliptic flow in intermediate energy heavy ion collisions and in-medium effects . Phys. Rev. C 65 , 064611 ( 2002 ). https://doi.org/10.1103/PhysRevC.65.064611 https://doi.org/10.1103/PhysRevC.65.064611 Elliptic flow in intermediate energy heavy ion collisions and in-medium effects V. R. Pandharipande and S. C. Pieper , Nuclear transparency to intermediate-energy nucleons from (e,e’p) reactions . Phys. Rev. C 45 , 791 ( 1991 ). https://doi.org/10.1103/PhysRevC.45.791 https://doi.org/10.1103/PhysRevC.45.791 Nuclear transparency to intermediate-energy nucleons from (e,e’p) reactions B. A. Li and L. W. Chen , Nucleon-nucleon cross sections in neutron-rich matter and isospin transport in heavy-ion reactions at intermediate energies . Phys. Rev. C 72 , 064611 ( 2005 ). https://doi.org/10.1103/PhysRevC.72.064611 https://doi.org/10.1103/PhysRevC.72.064611 Nucleon-nucleon cross sections in neutron-rich matter and isospin transport in heavy-ion reactions at intermediate energies Z. Q. Feng , Nuclear in-medium effects and collective flows in heavy-ion collisions at intermediate energies . Phys. Rev. C 85 , 014604 ( 2012 ). https://doi.org/10.1103/PhysRevC.85.014604 https://doi.org/10.1103/PhysRevC.85.014604 Nuclear in-medium effects and collective flows in heavy-ion collisions at intermediate energies Z. Q. Feng , Nuclear dynamics and particle production near threshold energies in heavy-ion collisions . Nucl. Sci. Tech . 29 , 40 ( 2018 ). https://doi.org/10.1007/s41365-018-0379-z https://doi.org/10.1007/s41365-018-0379-z Nuclear dynamics and particle production near threshold energies in heavy-ion collisions P.-C. Li , Y.-J. Wang , Q.-F. Li , et al ., Collective flow and nuclear stopping in heavy ion collisions in Fermi energy domain . Nucl. Sci. Tech . 29 , 177 ( 2018 ). https://doi.org/10.1007/s41365-018-0510-1 https://doi.org/10.1007/s41365-018-0510-1 Collective flow and nuclear stopping in heavy ion collisions in Fermi energy domain J. Su , C.Y. Huang , W.J. Xie , et al ., Effects of in-medium nucleon-nucleon cross sections on stopping observable and ratio of free protons in heavy-ion collisions at 400 MeV/nucleon . Eur. Phys. J. A 52 , 207 ( 2016 ). https://doi.org/10.1140/epja/i2016-16207-x https://doi.org/10.1140/epja/i2016-16207-x Effects of in-medium nucleon-nucleon cross sections on stopping observable and ratio of free protons in heavy-ion collisions at 400 MeV/nucleon F. Zhang and J. Su , A dynamical description of the 136Xe+ p spallation at 1000 MeV/nucleon . Chinese Physics C 43 , 024103 ( 2019 ). https://doi.org/10.1088/1674-1137/43/2/024103 https://doi.org/10.1088/1674-1137/43/2/024103 A dynamical description of the 136Xe+ p spallation at 1000 MeV/nucleon 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. Rep . 202 , 233 ( 1991 ). https://doi.org/10.1016/0370-1573(91)90094-3 https://doi.org/10.1016/0370-1573(91)90094-3 Quantum molecular dynamics a dynamical microscopic n-body approach to investigate fragment formation and the nuclear equation of state in heavy ion collisions G. D. Westfall , W. Bauer , D. Craig et al ., Mass dependence of the disappearance of flow in nuclear collisions . Phys. Rev. Lett . 71 , 1986 ( 1993 ). https://doi.org/10.1103/PhysRevLett.71.1986 https://doi.org/10.1103/PhysRevLett.71.1986 Mass dependence of the disappearance of flow in nuclear collisions 10.1007/s41365-020-00787-6 This work was supported by the National Natural Science Foundation of China (Nos. 11905018 and 11875328), and the Scientific and Technological Innovation Programs of Higher Education Institutions of Shanxi Province, China (No. 2019L0908). 2020-03-16 2020-05-30 2020-06-02 2020-07-29 2020-08-01 2021-04-10 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2020 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-020-00787-6 http://dx.doi.org/10.1007/s41365-020-00787-6 . 2020 Nuclear Science and Techniques 8 31