Isospin effects on intermediate mass fragments at intermediate energy-heavy ion collisions

Isospin effects on intermediate mass fragments at intermediate energy-heavy ion collisions 1 first-author 1 1 2 corresp zhyx@ciae.ac.cn zhyx@ciae.ac.cn zhyx@ciae.ac.cn https://orcid.org/0000-0003-1860-4745 China Institute of Atomic Energy, P. O. Box 275(18), Beijing 102413, China China Institute of Atomic Energy, P. O. Box 275(18), Beijing 102413, China Department of Physics and Technology, Guangxi normal University, Guilin 540101, China Department of Physics and Technology, Guangxi normal University, Guilin 540101, China In this study, we investigated the isospin properties of intermediate mass fragments (IMFs) for the central collisions of 112,124 Sn+ 112,124 Sn at a beam energy of 50 MeV per nucleon using an improved quantum molecular dynamics model (ImQMD) coupled with a sequential decay model (GEMINI). Three observables were analyzed: 1) the average center-of-mass kinetic energy per nucleon <math id="M1"><mrow><mo>〈</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo>/</mo><mi>A</mi><mo>〉</mo></mrow></math> of fragments as a function of their charge number Z ; 2) the average neutron number to proton number ratio ( <math id="M2"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> ) of fragments with a given charge number Z as a function of their center-of-mass kinetic energy per nucleon ( E c.m. / A ); and 3) the average total neutron number to total proton number ratio ( <math id="M3"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> ) and double ratio ( DR ( N/Z )) of IMFs with Z =3-8 as a function of their center-of-mass kinetic energy per nucleon E c.m. / A . Our calculations revealed that the sensitivity of the isospin properties of IMFs relative to the stiffness of the symmetry energy remains even after sequential decay. By comparing the calculations of <math id="M4"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> and DR ( N/Z ) with the data, it was found that the soft symmetry energy, i.e., γ =0.5, is favored. Symmetry energy Intermediate mass fragments Isospin effects Introduction 1 The Equation of State (EOS) of asymmetric nuclear matter can be written approximately as <math display="block" id="M5"><mrow><mi>E</mi><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>δ</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>E</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>δ</mi><mo>=</mo><mn>0</mn><mo stretchy="false">)</mo><mo>+</mo><mi>S</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><msup><mi>δ</mi><mn>2</mn></msup><mo>+</mo><mi>O</mi><mrow><mo>(</mo><mrow><msup><mi>δ</mi><mn>4</mn></msup></mrow><mo>)</mo></mrow><mo>+</mo><mo>…</mo><mn>,</mn></mrow></math> where δ =( ρ n - ρ p )/( ρ n + ρ p ), which is related to different areas of nuclear physics, such as the nuclear structure of ground-state nuclei, the dynamics of neutron-rich heavy ion collisions (HICs), the physics of neutron stars, etc. [ 1 1 , 2 2 ]. To date, the isoscalar component of the EOS has been well-constrained by the kaon condensation and the flow data of HICs [ 3 3 ]. However, the symmetry energy, particularly the density dependence of the symmetry energy S ( ρ ), still has large uncertainties away from the normal density. Heavy-ion collisions with neutron-rich nuclei in facilities, such as the Heavy Ion Research Facility in Lanzhou (HIRFL/Lanzhou), Facility for Rare Isotope Beams at Michigan States University (FRIB/MSU), the Radioactive Isotope Beam Factory (RIBF/RIKEN), will produce rare isotope beams to study the dynamical evolution of neutron-rich nuclear systems and obtain information on the density dependence of the symmetry energy by comparing the acquired data with the transport model simulations [ 4 4 , 5 5 ]. Many theoretical and experimental studies have been conducted to constrain the density dependence of the symmetry energy via intermediate-energy HICs [ 2 2 , 6 6 - 20 20 6-20 ]. In addition, Zheng et al. proposed a novel method for studying the symmetry energy [ 21 21 ], which may impose additional restriction on the constraints in the future. There are also a lot of studies that aim to constrain the density dependence of symmetry energy at the suprasaturation density by analyzing the properties of neutron stars, such as the mass-radius relationship and tidal deformability [ 23 23 - 34 34 23-34 ], subsequent to the publication of GW170817 events data [ 23 23 , 24 24 ]. By analyzing the combination of the observables of the neutron skin of 208 Pb, the collective flow, isospin diffusion, or π - / π + ratios of HICs, and properties of neutron stars, the relatively tight constraints on the symmetry energy were inferred from 0.3 ρ 0 to 3 ρ 0 [ 14 14 , 17 17 , 12 12 , 22 22 ]. In a recent review paper by Li et al.[ 35 35 ], the constraints on the symmetry energy were compiled and yielded the average values and their variances. For the slope of the symmetry energy, i.e. <math id="M6"><mrow><mi>L</mi><mo>=</mo><mn>3</mn><msub><mi>ρ</mi><mn>0</mn></msub><mfrac><mrow><mo>∂</mo><mi>S</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><mi>ρ</mi></mrow></mfrac><msub><mo>|</mo><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></msub></mrow></math> , the value is 57.7± 19 MeV, and the symmetry energy at twice the normal density S (2 ρ 0 ) is 51± 13 MeV. However, there is tension around the symmetry energy constraints following the Lead Radius Experiment (PREX-II) reporting of a model-independent extraction of the neutron skin of 208 Pb and Δ r np =0.29± 0.07 fm [ 36 36 ]. Yue et al. [ 17 17 ] analyzed newly published neutron skin data with ground-state properties and giant monopole resonances of finite nuclei, the properties of neutron stars for PSR J0740+6620, the mass-radius of PSR J0030+045, and the neutron star tidal deformability extracted from the GW170817 event to obtain the constraints of symmetry energy in a large density region. Their calculations suggested that S ( ρ 0 )=34.5± 1.5 MeV and L =85.5 ± 22.2 MeV. The values obtained are consistent with the constraints of the symmetry energy from the S π RIT data [ 13 13 ] for Sn+Sn at 270A MeV, but are less than the values obtained for the analysis of PREX-II using a specific class of relativistic energy density functional [ 36 36 ], in which L =106 ± 37 MeV and S 0 =38.7± 4.7 MeV. Their constraints on the symmetry energy also overlap with the symmetry energy extracted from the charge radius of 54 Ni, where 21 MeV < L < 88 MeV [ 37 37 ], and the combination of astrophysical data with PREX-II and chiral effective field theory, where <math id="M7"><mrow><mi>L</mi><msubsup><mrow><mo>=</mo><mn>53</mn></mrow><mrow><mo>−</mo><mn>15</mn></mrow><mrow><mo>+</mo><mn>14</mn></mrow></msubsup></mrow></math> MeV [ 38 38 ]. Currently, the big challenge in constraining the density dependence of the symmetry energy with HICs is to understand the model dependence and reduce the uncertainties of the symmetry energy constraints. The understanding of the model dependence has been achieved in the last 10 years via the transport model evaluation project, and significant progress has been made in the treatment of nucleon-nucleon collisions and the mean-field potential [ 39 39 - 43 43 39-43 ]. The reduction of the uncertainties of symmetry energy constraints is still underway, which can be realized by simultaneously describing multi-observables using transport models. The data for isospin-sensitive observables from intermediate-energy HICs, such as isoscaling [ 44 44 ], double neutron to proton ratio [ 45 45 ], isospin diffusion [ 46 46 ], rapidity dependence of the isobaric ratio [ 10 10 ], and neutron-excess to proton ratio [ 11 11 ], have been used for comparison with the predictions from different transport models [ 4 4 , 5 5 ]. Among these transport models, the ImQ MD model successfully describes the neutron-to-proton yield ratio (n/p ratio), isospin diffusion, and isospin transport ratio as a function of rapidity, and constraints on the symmetry energy at the subsaturation density have been achieved [ 5 5 ]. Except for the aforementioned isospin-sensitive observables, the neutron number divided by the proton number of IMFs was considered to be complementary information of the neutron-to-proton yield ratio and can be used to extract information related to the symmetry energy [ 6 6 ]. This can be explained by the multi-fragmentation mechanism where the copious production of IMFs and its isospin contents are closely related to the dynamic mechanism of the HICs [ 47 47 - 49 49 47-49 ]. The measurement of IMFs has been performed at the National Superconducting Conducting Cyclotron Laboratory at Michigan State University [ 50 50 ], but there are few models for the theoretical analysis of the isospin effects of IMFs and their use to constrain the symmetry energy. Thus, it is worth to check whether the ImQMD model can simultaneously describe the isospin properties of IMFs. In this study, we focus on the properties of IMFs ( Z =3-8) for the central collision ( b =2 fm) of 112,124 Sn+ 112,124 Sn at E beam =50 MeV/u using the ImQMD model coupled with the GEMINI model. Furthermore, we compared the calculated results to the data for learning the form of the symmetry energy and the reaction dynamics. Theoretical model and symmetry energy 1 The version of ImQMD used in this study was the same as that used in Ref. [ 5 5 ], i.e., ImQMD05, in which the isospin-independent momentum-dependent interaction was utilized. To simulate low-intermediate energy HICs, the basic properties of the initial nuclei are important, and they are guaranteed by requiring that the binding energy and root-mean-square (rms) radius of the sampled nuclei fall within the range of the experimental data with specific uncertainties. For the standard initialization treatment, the initial nuclei remain stable within 200 fm/ c , which is sufficient for Sn+Sn at 50 MeV/u because the reaction dynamics are completed within 100 fm/ c [ 51 51 ]. Additional details regarding the initialization of the ImQMD model can be found in Refs. [ 51 51 - 55 55 51-55 ]. The nucleonic mean field potential is derived from the Skyrme-type energy density functional, and the energy density of the symmetry potential is as follows: <math display="block" id="M8"><mrow><msubsup><mi>w</mi><mrow><mtext>asy</mtext></mrow><mrow><mtext>pot</mtext></mrow></msubsup><mo>=</mo><mfrac><mrow><msub><mi>C</mi><mtext>s</mtext></msub></mrow><mn>2</mn></mfrac><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mi>γ</mi></msup><msup><mi>δ</mi><mn>2</mn></msup><mi>ρ</mi><mn>.</mn></mrow></math> where ρ is the nucleon density. C s is the symmetry potential parameter, C s =35.2 MeV, and γ is a parameter related to the density dependence of the symmetry energy. Based on these interaction parameters, we can obtain the symmetry energy per nucleon for cold nuclear matter, as follows: <math display="block" id="M9"><mrow><mi>S</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mfrac><mrow><msup><mi>ℏ</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><mi>m</mi></mrow></mfrac><msubsup><mi>ρ</mi><mn>0</mn><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msubsup><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mn>3</mn><msup><mi>π</mi><mn>2</mn></msup></mrow><mn>2</mn></mfrac><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>+</mo><mfrac><mrow><msub><mi>C</mi><mi>s</mi></msub></mrow><mn>2</mn></mfrac><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mi>ρ</mi><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mi>γ</mi></msup><mn>,</mn></mrow></math> where m denotes the nucleon mass. The corresponding values of the symmetry energy S ( ρ 0 ) and the slope of the symmetry energy <math id="M10"><mrow><mi>L</mi><mo>=</mo><mn>3</mn><msub><mi>ρ</mi><mn>0</mn></msub><mfrac><mrow><mo>∂</mo><mi>S</mi><mo stretchy="false">(</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><mi>ρ</mi></mrow></mfrac><msub><mo>|</mo><mrow><msub><mi>ρ</mi><mn>0</mn></msub></mrow></msub></mrow></math> at a normal density are listed in Table 1 Table 1 . For this particular parameterization, the symmetry energy increases with a decrease of γ at a subsaturation density, whereas the opposite is true for supranormal density. 1 1 Para. 1 1 S 0 1 1 L 1 1 S (2 ρ 0 ) 1 1 γ =0.5 1 1 30 1 1 51 1 1 44.4 1 1 γ =2.0 1 1 30 1 1 130 1 1 90.0 The ImQMD simulations were terminated at 400 fm/ c when the systems reached “dynamical freeze-out”. The primary fragments were recognized using the isospin-dependent minimum spanning tree method [ 51 51 ], and the excitation energy was calculated in the rest frame of the fragment by subtracting the ground state energy from the energy of the fragment in its rest frame. It should be noted that the fragments produced during the ImQMD calculations are primary, and they may have very neutron-rich or very neutron-poor “unphysical” fragments. The probability of generation of these “unphysical” primary fragments is tiny in ImQMD simulations, and their ground state energies are approximately obtained using the liquid-drop mass formula. The sequential decay of the primary fragments was performed using the GEMINI code [ 56 56 , 57 57 ], which describes the de-excitation process of the compound nucleus and facilitates binary decay, including light-particle evaporation and symmetric fission. The decay continues until the particle emission is energetically forbidden or impossible owing to competition with γ -ray emission. Results and discussions 1 Before studying the isospin properties of the IMFs, we focused on the cluster formation mechanism of HICs. Figure 1 Figure 1 shows the time evolution of the multiplicities of emitted nucleons (top panels) and fragments (bottom panels), for the collisions of 124 Sn+ 124 Sn and 112 Sn+ 112 Sn. It is evident that the nucleons are emitted prior to IMFs with Z =3-8, and IMFs with Z =3-8 prior to fragments with Z =9-20. For example, for the reaction system 112 Sn+ 112 Sn, the emission of nucleons begins from t ∼50 fm/ c , but the emission of IMFs with Z =3-8 occurs at approximately t ∼100 fm/ c and Z =9-20 fragments at approximately t ∼140 fm/ c . Such a dynamic process results in a decrease in the fragment velocity with the fragment charge number. After t ∼250 fm/ c , the primary fragments are completely formed. Ideally, the system should be in “dynamical freeze-out,” and its multiplicities should not change with time. However, it is difficult for the QMD type Hamiltonian to bind the fragments for a long time, and the pseudo-emission of nucleons occurs at the late stage of the simulations of HICs. Therefore, the multiplicities of fragments with Z =3-8 or 9-20 decrease slightly, and the magnitude depends on the stiffness of the symmetry energy. To understand the influence of different symmetry potentials on the reaction dynamics, we investigated the time evolution of the multiplicities of nucleons and fragments for γ =0.5 (solid lines) and γ =2.0 (dashed lines). The results show that the multiplicities of the nucleons and the fragments depend on γ and the isospin asymmetry of the reaction system. The multiplicities of emitted neutrons calculated with γ =0.5 are greater than those with γ =2.0, and the effect is the opposite for the multiplicities of protons. Consequently, it is expected that the neutron-to-proton yield ratio would be sensitive to the stiffness of the symmetry energy. Numerous studies have investigated for constraining the symmetry energy using n/p ratios [ 45 45 , 10 10 , 54 54 , 58 58 , 59 59 , 39 39 ]. For the multiplicities of fragments with Z =3-8 (or 9-20), the calculations show that the multiplicities of IMFs obtained with γ =0.5 are smaller than those with γ =2.0. This is because a higher density is reached in the overlap region and more compressional energy is stored for the soft isospin asymmetric nuclear equation of the state ( γ =0.5) compared to that of the stiff isospin asymmetric nuclear equation of state ( γ =2.0) in the simulations. Thus, the reaction system simulated with γ =0.5 disintegrates into more light-charged fragments compared to the case of γ =2.0. Owing to the conservation of the nucleon number in the reaction system, the multiplicities for fragments with Z =3-8 (or 9-20) for the γ =0.5 case are smaller than those for γ =2.0. Considering the reaction dynamics, it is expected that the size of the fragments and their kinetic energies should be correlated to the stiffness of the symmetry energy. In Fig. 2 Fig. 2 , the average center-of-mass kinetic energy per nucleon <math id="M11"><mrow><mo>〈</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo>/</mo><mi>A</mi><mo>〉</mo></mrow></math> is plotted for the primary fragments as a function of their charge number Z for 124 Sn+ 124 Sn (left panel),and 112 Sn+ 112 Sn (right panel). The solid red lines represent the calculated results for γ =0.5, and the dashed red lines are the results obtained for γ =2.0. For both symmetry energies, the values of <math id="M12"><mrow><mo>〈</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo>/</mo><mi>A</mi><mo>〉</mo></mrow></math> decrease with an increase in Z owing to the aforementioned reaction dynamics. It is interesting to note that <math id="M13"><mrow><mo>〈</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo>/</mo><mi>A</mi><mo>〉</mo></mrow></math> as a function of Z is also sensitive to the density dependence of the symmetry energy, and the values of <math id="M14"><mrow><mo>〈</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo>/</mo><mi>A</mi><mo>〉</mo></mrow></math> obtained for γ =0.5 are larger than those obtained for γ =2.0 when Z > 3. This can be understood based on the reaction dynamics shown in Fig. 1 Fig. 1 , wherein the multiplicities of the IMFs are larger for γ =2.0 compared to the case of γ =0.5. If there is an energy equipartition about fragments, one can expect each IMF to carry less kinetic energy for the case of γ =2.0 compared to the case of γ =0.5. The difference between the values of <math id="M15"><mrow><mo>〈</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo>/</mo><mi>A</mi><mo>〉</mo></mrow></math> obtained for γ =0.5 and γ =2.0 increase with an increase in Z , and the difference exceeds 10% at Z =8. However, the measured fragments in the experiments cold, which indicate that the hot primary fragments should be deexcited. In this study, the decay of the primary fragments created in the ImQMD simulations was performed using the GEMINI code [ 56 56 , 57 57 ] using their default parameter sets. The blue lines with symbols represent the results for the cold fragments. For Z =3 and 4, the values of <math id="M16"><mrow><mo>〈</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo>〉</mo><mo>/</mo><mi>A</mi></mrow></math> are reduced by ∼ 7% compared to the results obtained for the primary fragments. This reduction is caused by the decay of primary fragments, which converts the fragments with higher mass or charge numbers to those with lower mass or charge numbers. For fragments with Z ≥ 6, the effects of decay on <math id="M17"><mrow><mo>〈</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo>/</mo><mi>A</mi><mo>〉</mo></mrow></math> are weak, and the sensitivities of <math id="M18"><mrow><mo>〈</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo>/</mo><mi>A</mi><mo>〉</mo></mrow></math> to the stiffness of the symmetry energy remain. Consequently, <math id="M19"><mrow><mo>〈</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo>/</mo><mi>A</mi><mo>〉</mo></mrow></math> as a function of Z of the IMFs can also be used to understand the fragmentation mechanism and to probe the density dependence of the symmetry energy. We will now discuss the isospin contents of the IMFs in more detail. The <math id="M20"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> of the IMFs was investigated, which is defined as <math display="block" id="M21"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi><mo>=</mo><mfrac><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mi>i</mi></msub><mrow><msub><mi>N</mi><mi>i</mi></msub><mo>×</mo><mi>Y</mi><mo stretchy="false">(</mo><msub><mi>N</mi><mi>i</mi></msub><mn>,</mn><mi>Z</mi><mo stretchy="false">)</mo></mrow></mstyle></mrow><mrow><mstyle displaystyle="true"><msub><mo>∑</mo><mi>i</mi></msub><mrow><mi>Y</mi><mo stretchy="false">(</mo><msub><mi>N</mi><mi>i</mi></msub><mn>,</mn><mi>Z</mi><mo stretchy="false">)</mo></mrow></mstyle></mrow></mfrac><mo>/</mo><mi>Z</mi><mn>.</mn></mrow></math> where Y ( N i,Z ) is the yield of the fragment with a neutron number N i and charge number Z . The summation is for all isotopes at a given Z . As an example, Fig. 3 Fig. 3 present the <math id="M22"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> ratio of primary fragments (red lines) as a function of E c.m. / A for fragments with Z =3, 4, 5, 6, 7, 8. The left six panels are for 124 Sn+ 124 Sn and the right six panels are for the reaction system 112 Sn+ 112 Sn. In general, the slopes of the <math id="M23"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> ratio as a function of E c.m. / A are negative, because the Coulomb energy accelerates the neutron-poor fragments to a higher kinetic energy region. These behavior occurs for all elements of the IMFs. Furthermore, the curves of the <math id="M24"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> ratios as a function of kinetic energy E c.m. / A depend on the stiffness of the symmetry energy, and the results obtained for γ =0.5 and 2.0 cross over at a specific energy, E cross . The values of E cross increase from 4 MeV to approximately 10 MeV as Z increases from 3 to 8, owing to the Coulomb interaction. Below the cross energy, the values of <math id="M26"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> obtained for γ =0.5 are smaller than those obtained for γ =2.0, and turn over above E cross . The predicted magnitude of the slope of the <math id="M27"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> ratio as a function of E c.m. / A obtained for γ =0.5 is smaller compared to that for γ =2.0. This behavior can be understood based on the cluster formation mechanism. For the Z ≤4 fragments with high kinetic energy, and the isospin effects were governed by the dynamics. Thus, the behavior of <math id="M28"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> of the fragments in the high kinetic energy region is similar to the n/p ratios of the emitted nucleons. Consequently, one can expect the values of <math id="M29"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> obtained for γ =0.5 to be larger than those obtained for γ =2.0. Fragments with low kinetic energy or heavier fragments are formed at a later stage of the reaction, and their isospin content is determined by their free energy. The free energy of the primary fragments, f ( A, δ, T ), can be approximated as <math display="block" id="M30"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>A</mi><mn>,</mn><mi>δ</mi><mn>,</mn><mi>T</mi><mo stretchy="false">)</mo><mo>=</mo></mrow></mtd><mtd columnalign="left"><mrow><mo>−</mo><msubsup><mi>a</mi><mrow><mtext>sym</mtext></mrow><mi>v</mi></msubsup><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>T</mi><mo stretchy="false">)</mo><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mi>Z</mi><msup><mo stretchy="false">)</mo><mn>2</mn></msup><mo>/</mo><mi>A</mi><mo>+</mo><msubsup><mi>a</mi><mrow><mtext>sym</mtext></mrow><mrow><mtext>surf</mtext></mrow></msubsup><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>T</mi><mo stretchy="false">)</mo><msup><mi>δ</mi><mn>2</mn></msup><msup><mi>A</mi><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>−</mo><msub><mi>a</mi><mtext>c</mtext></msub><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>T</mi><mo stretchy="false">)</mo><mi>Z</mi><mo stretchy="false">(</mo><mi>Z</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo><mo>/</mo><msup><mi>A</mi><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>+</mo><msub><mi>a</mi><mtext>v</mtext></msub><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>T</mi><mo stretchy="false">)</mo><mi>A</mi><mo>−</mo><msub><mi>a</mi><mtext>s</mtext></msub><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>T</mi><mo stretchy="false">)</mo><msup><mi>A</mi><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup><mn>.</mn></mrow></mtd></mtr></mtable></mrow></math> where <math id="M31"><mrow><msubsup><mi>a</mi><mrow><mtext>sym</mtext></mrow><mtext>v</mtext></msubsup></mrow></math> and <math id="M32"><mrow><msubsup><mi>a</mi><mrow><mtext>sym</mtext></mrow><mrow><mtext>surf</mtext></mrow></msubsup></mrow></math> are the bulk-and surface-related symmetry energy coefficients, respectively. a c is a coefficient related to the Coulomb energy, and a v and a s are the volume and surface energy coefficients, respectively. If we assume that the fragments are at thermal equilibrium and stable, the isospin asymmetry of the primary fragment is determined by the minimum free energy of the cluster at temperature T , i.e., <math id="M33"><mrow><mfrac><mrow><mo>∂</mo><mi>f</mi><mo stretchy="false">(</mo><mi>A</mi><mn>,</mn><mi>δ</mi><mn>,</mn><mi>T</mi><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><mi>δ</mi></mrow></mfrac><mo>=</mo><mn>0</mn></mrow></math> . Therefore,the isospin asymmetry of the primary fragments should be proportional to <math display="block" id="M34"><mrow><msub><mi>δ</mi><mrow><mtext>frag</mtext></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>a</mi><mtext>c</mtext></msub><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>T</mi><mo stretchy="false">)</mo><msup><mi>A</mi><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup><mrow><mo>(</mo><mrow><mi>A</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow><mrow><msub><mi>a</mi><mtext>c</mtext></msub><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>T</mi><mo stretchy="false">)</mo><msup><mi>A</mi><mrow><mn>5</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>+</mo><mn>4</mn><msubsup><mi>a</mi><mrow><mtext>sym</mtext></mrow><mtext>v</mtext></msubsup><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>T</mi><mo stretchy="false">)</mo><mi>A</mi><mo>−</mo><mn>4</mn><msubsup><mi>a</mi><mrow><mtext>sym</mtext></mrow><mrow><mtext>surf</mtext></mrow></msubsup><msup><mi>A</mi><mrow><mn>2</mn><mo>/</mo><mn>3</mn></mrow></msup></mrow></mfrac><mn>.</mn></mrow></math> The Coulomb interaction increases the isospin asymmetry of the primary fragments and the symmetry energy term decreases the isospin asymmetry of the fragments. Since the fragments are formed during the expansion phase during which the system is at a subsaturation density, the values of the symmetry energy <math id="M35"><mrow><msubsup><mi>a</mi><mrow><mtext>sym</mtext></mrow><mi>v</mi></msubsup><mo stretchy="false">(</mo><mi>ρ</mi><mn>,</mn><mi>T</mi><mo stretchy="false">)</mo></mrow></math> obtained for γ =0.5 are larger than those obtained for γ =2.0, and it is expected that the <math id="M36"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> values of the fragments obtained for γ =2.0 will be larger than those for γ =0.5, according to Eq. (6). To investigate the secondary sequential decay effects, we present <math id="M37"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> as a function of E c.m. / A as blue lines in Fig. 3 Fig. 3 . For fragments with Z ≥5, the secondary decay reduces the sensitivity of <math id="M38"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> to the stiffness of the symmetry energy because of abundant neutron evaporation. Specifically, the decay weakly influences the isospin contents of the fragments obtained for γ =0.5, and strongly suppresses the <math id="M39"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> obtained for γ =2.0. This is because the E^*/A obtained for γ =0.5 is smaller than that obtained for γ =2.0. For γ =0.5, the values of E^*/A varied from 1.6 to 1.8 MeV. For γ =2.0, the values of E^*/A varied from 1.8 to 2.4 MeV. The hardly understanding point is that the sensitivity of <math id="M40"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> to the stiffness of the symmetry energy is the same or is enhanced for Z =3 and 4 after the decay. This could be related to the “ 3 He- 4 He puzzle”[ 60 60 - 62 62 60-62 ] and the decay of 8 Be[ 50 50 ]. Thus, the cluster formation mechanism in transport models and the sequential decay model require further refinement. To avoid the cluster formation flaw in the transport model, we calculated <math id="M41"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> for fragments with Z from 3 to 8. <math id="M42"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi></mrow></math> and <math id="M43"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> are defined as follows: <math display="block" id="M44"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>i</mi><mn>,</mn><mi>j</mi></mrow></munder></mstyle><msub><mi>N</mi><mi>i</mi></msub><mo>×</mo><mi>Y</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mi>j</mi></msub><mn>,</mn><msub><mi>N</mi><mi>i</mi></msub><mo stretchy="false">)</mo><mo>;</mo></mrow></math> <math display="block" id="M45"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>i</mi><mn>,</mn><mi>j</mi></mrow></munder></mstyle><msub><mi>Z</mi><mi>j</mi></msub><mo>×</mo><mi>Y</mi><mo stretchy="false">(</mo><msub><mi>Z</mi><mi>j</mi></msub><mn>,</mn><msub><mi>N</mi><mi>i</mi></msub><mo stretchy="false">)</mo><mn>.</mn></mrow></math> The summation is over all IMFs for a given kinetic energy bin (1.5 MeV/nucleon) for all events. This is similar to the method proposed in Ref. [ 6 6 ] for studying the isospin effects and symmetry energy in HICs. The data on the isospin observable <math id="M46"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> for the central collision of Sn+Sn were reported in Ref.[ 50 50 ]. It should be noted that the experimental central collisions are defined by the high multiplicity of charged particles, and the corresponding real impact parameter distribution is Gaussian [ 63 63 - 65 65 63-65 ]. For the central collision events selected in Refs.[ 50 50 ], the reduced impact parameter ranges from 0 to approximately 0.4 and peakes at approximately 0.2. Thus, we performed simulations at a representative value of the impact parameter, such as the mean of the impact parameter distribution, b =2 fm. We also selected fragments emitted at center-of-mass angles of 70° ≤ θ c.m. ≤ 110°, which is the same as that in Ref.[ 50 50 ]. The calculated results are shown in Fig. 4 (a)-(d) Fig. 4 (a)-(d) . Panels (a) and (c) show 124 Sn+ 124 Sn, and panels (b) and (d) show 112 Sn+ 112 Sn. The calculated results for <math id="M47"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> for the primary fragments (top panels) strongly depend on the density dependence of the symmetry, which is consistent with the results in Ref.[ 6 6 ]. However, <math id="M48"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> as a function of E c.m. / A obtained in the transverse direction is weakly sensitive to the isospin asymmetry of the system. This is because the IMFs preferentially form in the forward and backward regions owing to the IMFs formation mechanism. In panels (e) and (f), we present <math id="M49"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> as a function of E c.m. / A in the non-transverse direction region, i.e., θ c.m. < 70° and θ c.m. > 110°, and <math id="M50"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> in the non-transverse direction not only depend on the stiffness of the symmetry energy, but also on the isospin asymmetry of the systems. The results of <math id="M52"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> of the IMFs after decay are plotted in panels (c) and (d) for the transverse direction and in panels (g) and (h) for the non-transverse direction. The secondary sequential decay mainly influences the <math id="M53"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> ratios in the lower kinetic energy region because the fragments with lower kinetic energy mainly consist of neutron-rich isotopes. However, the <math id="M54"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> ratios still maintain a sensitivity to the density dependence of the symmetry energy after sequential decay. In the transverse direction, the calculations indicate that the sensitivity of <math id="M55"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> increases as E c.m. / A increases. We compared the calculated results in the transverse direction with the data (solid points in panels (c) and (d))[ 50 50 ], which includes the 8 Be contribution, to understand the model prediction. It is evident that the calculations with γ =0.5 are favored, and this conclusion is consistent with the constraints on γ obtained with the n/p ratio, isospin diffusion, and isospin transport ratio as a function of rapidity[ 5 5 , 54 54 ]. For <math id="M56"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> as a function of E c.m. / A in the non-transverse direction, our calculations show that the values of <math id="M57"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> are not only sensitive to the symmetry energy, but also the isospin asymmetry of the systems. However, the sensitivity of <math id="M58"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> in the non-transverse direction is weaker compared to that in the transverse direction. Finally, we analyze the double ratios of <math id="M59"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> for the two systems A= 124 Sn+ 124 Sn and B= 112 Sn+ 112 Sn, i.e., <math display="block" id="M60"><mrow><mi>D</mi><mi>R</mi><mo stretchy="false">(</mo><mi>N</mi><mo>/</mo><mi>Z</mi><mo stretchy="false">)</mo><mo>=</mo><mo stretchy="false">(</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi><msub><mo stretchy="false">)</mo><mi>A</mi></msub><mo>/</mo><mo stretchy="false">(</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi><msub><mo stretchy="false">)</mo><mi>B</mi></msub><mn>.</mn></mrow></math> The idea of this double ratio is similar to the double neutron-to-proton ratio, which has been used to isolate the isospin effects and constrain the density dependence of the symmetry energy [ 54 54 , 6 6 , 45 45 ]. To obtain good statistics of the results of DR ( N/Z ) in the transverse direction, we performed simulations with 100,000 events. In Fig. 5 (a) Fig. 5 (a) , we present the DR ( N/Z ) ratios for the cold fragments in the transverse direction for an experimental case such that 70°≤ θ c.m. ≤ 110°. The shaded region with a solid pattern (dashed pattern) is the result obtained for γ =0.5 ( γ =2.0). The DR ( N/Z ) ratios obtained for γ =2.0 are larger than those obtained for the γ =0.5 case. This can be observed from Fig. 3 Fig. 3 . These trends reflect the fact that the fragments originate from the composite system that survives the pre-equilibrium emission, as indicated in Ref.[ 6 6 ]. Compared to the data, the calculations for the case of γ =0.5 are favored. We also analyzed DR ( N/Z ) in the non-transverse direction, and the results are shown in panel (b). The calculations show that the sensitivity of DR(N/Z) to the stiffness of the symmetry energy is maintained in the non-transverse direction, and the DR(N/Z) in the non-transverse direction has good statistics because the IMFs are distributed in the forward and backward regions. Summary and prospective 1 In summary, we investigated the isospin properties of intermediate mass fragments for the reaction 112,124 Sn+ 112,124 Sn at a beam energy of 50 MeV per nucleon. Our results indicated that the average kinetic energy per nucleon of the fragments <math id="M61"><mrow><mo>〈</mo><msub><mi>E</mi><mrow><mtext>c</mtext><mtext>.m</mtext><mtext>.</mtext></mrow></msub><mo>/</mo><mi>A</mi><mo>〉</mo></mrow></math> as a function of Z is sensitive to the density dependence of the symmetry energy, and this sensitivity does not change significantly by secondary sequential decay. Furthermore, the <math id="M62"><mrow><mo>〈</mo><mi>N</mi><mo>〉</mo><mo>/</mo><mi>Z</mi></mrow></math> ratio of the IMFs as a function of the kinetic energy per nucleon is also sensitive to the density dependence of the symmetry energy. To circumvent the difficulties associated with cluster formation in transport model calculations, we also investigated the <math id="M63"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> and DR ( N/Z ) ratios as a function of kinetic energy. The calculations revealed that both observables are sensitive to the density dependence of the symmetry energy, and the sensitivity was maintained even after secondary sequential decay. By comparing the calculations to the obtained data, the symmetry energy with γ =0.5 is favored. The results are consistent with the conclusions obtained by simultaneously describing the n/p ratios, isospin diffusion, and isospin transport ratio as a function of rapidity. Nevertheless, the tight constraints on the symmetry energy still need more work. For example, the effects of the momentum-dependent symmetry potential should be well understood because the isospin-dependent momentum-dependent symmetry potential influences the cluster formation and the n/p ratio of the emitted nucleons [ 66 66 , 67 67 ]. Thus, the curves of the average kinetic energy of the IMFs as a function of Z are modified. The <math id="M64"><mrow><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>N</mi><mo>/</mo><mstyle displaystyle="true"><mo>∑</mo></mstyle><mi>Z</mi></mrow></math> ratios and DR ( N/Z ) may also be influenced by the different signs associated with effective mass splitting. These points should be carefully investigated in the future. Our analysis and experimental studies[ 50 50 ] also suggest that a sequential decay model should be developed. References Bao-An Li , Lie-Wen Chen , Che Ming Ko , Recent progress and new challenges in isospin physics with heavy-ion reactions . Phys. Rep. 464 , 113 ( 2008 ). doi: 10.1016/j.physrep.2008.04.005 http://doi.org/10.1016/j.physrep.2008.04.005 Recent progress and new challenges in isospin physics with heavy-ion reactions M. B. Tsang , J. R. Stone , F. Camera , et.al. , Constraints on the symmetry energy and neutron skins from experiments and theory . Phys. Rev. C 86 , 015803 ( 2012 ). doi: 10.1103/PhysRevC.86.015803 http://doi.org/10.1103/PhysRevC.86.015803 Constraints on the symmetry energy and neutron skins from experiments and theory P. Danielewicz , Flow and the equation of state of nuclear matter . Nucl. Phys. A 685 , 368c ( 2001 ). doi: 10.1016/S0375-9474(01)00554-1 http://doi.org/10.1016/S0375-9474(01)00554-1 Flow and the equation of state of nuclear matter Lie-Wen Chen , Che Ming Ko , and Bao-An Li , High-energy behavior of the nuclear symmetry potential in asymmetric nuclear matter . Phys. Rev. Lett. 94 , 032701 ( 2005 ). doi: 10.1103/PhysRevC.72.064606 http://doi.org/10.1103/PhysRevC.72.064606 High-energy behavior of the nuclear symmetry potential in asymmetric nuclear matter M. B. Tsang , Yingxun Zhang , P. Danielewicz , et. al. , Constraints on the Density Dependence of the Symmetry Energy . Phys. Rev. Lett. 102 , 122701 ( 2009 ). doi: 10.1103/PhysRevLett.102.122701 http://doi.org/10.1103/PhysRevLett.102.122701 Constraints on the Density Dependence of the Symmetry Energy M. Colonna , V. Banran , M. Di Toro , et. al. , Isospin distillation with radial flow: A test of the nuclear symmetry energy . Phys. Rev. C 78 , 064618 ( 2008 ). doi: 10.1103/PhysRevC.78.064618 http://doi.org/10.1103/PhysRevC.78.064618 Isospin distillation with radial flow: A test of the nuclear symmetry energy Z.G. Xiao , R.J. Hu , H.Y. Wu , et. al. , System dependence of the correlation function of IMFs in 36 Ar + 112,124 Sn at 35 MeV/u . Phys. Lett. B 639 , 436 ( 2006 ). doi: 10.1016/j.physletb.2006.06.076 http://doi.org/10.1016/j.physletb.2006.06.076 System dependence of the correlation function of IMFs in 36Ar +112,124Sn at 35 MeV/u E.W. Cornell , T.M. Hamilton , D. Fox , et.al. , Investigating the Evolution of Multifragmenting Systems with Fragment Emission Order . Phys. Rev. Lett. 77 , 4508 ( 1996 ). doi: 10.1103/PhysRevLett.77.4508 http://doi.org/10.1103/PhysRevLett.77.4508 Investigating the Evolution of Multifragmenting Systems with Fragment Emission Order E. Galichet , M. F. Rivet , B. Borderie , et. al. , Isospin diffusion in 58 Ni-induced reactions at intermediate energies. I. Experimental results Phys. Rev. C 79 , 064616 ( 2009 ). doi: 10.1103/PhysRevC.79.064614 http://doi.org/10.1103/PhysRevC.79.064614 Isospin diffusion in 58Ni-induced reactions at intermediate energies. I. Experimental results T. X. Liu , W. G. Lynch , M. B. Tsang , et. al. , Isospin diffusion observables in heavy-ion reactions . Phys. Rev. C 76 , 034603 ( 2007 ). doi: 10.1103/PhysRevC.76.034603 http://doi.org/10.1103/PhysRevC.76.034603 Isospin diffusion observables in heavy-ion reactions Yan Zhang , Junlong Tian , Wenjing Cheng , et. al. , Long-time drift of the isospin degree of freedom in heavy ion collisions . Phys. Rev. C 95 , 387 041602(R) ( 2017 ). doi: 10.1103/PhysRevC.95.041602 http://doi.org/10.1103/PhysRevC.95.041602 Long-time drift of the isospin degree of freedom in heavy ion collisions Yangyang Liu , Yongjia Wang , Ying Cui , et. al. , Insights into the pion production mechanism and the symmetry energy at high density . Phys. Rev. C 103 , 014616 ( 2021 ). doi: 10.1103/PhysRevC.103.014616 http://doi.org/10.1103/PhysRevC.103.014616 Insights into the pion production mechanism and the symmetry energy at high density J. Estee et al. (S π RIT Collaboration), Probing the Symmetry Energy with the Spectral Pion Ratio . Phys. Rev. Lett. 126 , 162701 . doi: 10.1103/PhysRevLett.126.162701 http://doi.org/10.1103/PhysRevLett.126.162701 Probing the Symmetry Energy with the Spectral Pion Ratio Yingxun Zhang , Min Liu , Cheng-Jun Xia , et. al. , Constraints on the symmetry energy and its associated parameters from nuclei to neutron stars . Phys. Rev. C 101 , 034303 ( 2020 ). doi: 10.1103/PhysRevC.101.034303 http://doi.org/10.1103/PhysRevC.101.034303 Constraints on the symmetry energy and its associated parameters from nuclei to neutron stars Bing Li , Na Tang , Feng-Shou Zhang , Isospin effects of projectile fragmentation in a Boltzmann-Langevin approach . Chin. Phys. C Vol. 45 , 394 No. 8 084103 ( 2021 ). doi: 10.1088/1674-1137/ac009a http://doi.org/10.1088/1674-1137/ac009a Isospin effects of projectile fragmentation in a Boltzmann-Langevin approach Jun Su , Constraining symmetry energy at subnormal density by isovector giant dipole resonances of spherical nuclei . Chin. Phys. C Vol. 43 , No. 6 , 064109 ( 2019 ). doi: 10.1088/1674-1137/43/6/064109 http://doi.org/10.1088/1674-1137/43/6/064109 Constraining symmetry energy at subnormal density by isovector giant dipole resonances of spherical nuclei Tong-Gang Yue , Lie-Wen Chen , Zhen Zhang , Ying Zhou , Constraints on the Symmetry Energy from PREX-II in the Multimessenger Era . arXiv:2102.05267 Constraints on the Symmetry Energy from PREX-II in the Multimessenger Era Zhao-Qing Feng . Nuclear dynamics and particle production near threshold energies in heavy-ion collisions . Nuclear Science and Techniques 29 ( 3 ): 40 ( 2018 ). doi: 10.1007/s41365-018-0379-z http://doi.org/10.1007/s41365-018-0379-z . Nuclear dynamics and particle production near threshold energies in heavy-ion collisions Hao Yu , De-Qing Fang , Yu-Gang Ma . Investigation of the symmetry energy of nuclear matter using isospin-dependent quantum molecular dynamics . Nuclear Science and Techniques 31 ( 6 ): 61 ( 2020 ) doi: 10.1007/s41365-020-00766-x http://doi.org/10.1007/s41365-020-00766-x . Investigation of the symmetry energy of nuclear matter using isospin-dependent quantum molecular dynamics Jie Liu , Chao Gao , Niu Wan , et al. Basic quantities of the Equation of State in isospin asymmetric nuclear matter . [J]. Nuclear Science and Techniques 32 ( 11 ): 117 ( 2021 ) doi: 10.1007/s41365-021-00955-2 http://doi.org/10.1007/s41365-021-00955-2 . Basic quantities of the Equation of State in isospin asymmetric nuclear matter H. Zheng , S. Burrello , M. Colonna , V. Baran , Pre-equilibrium effects in charge-asymmetric low-energy reactions . Phys. Lett. B 769 , 424 ( 2017 ). doi: 10.1016/j.physletb.2017.04.002 http://doi.org/10.1016/j.physletb.2017.04.002 Pre-equilibrium effects in charge-asymmetric low-energy reactions Jun Xu , Constraining Isovector Nuclear Interactions with Giant Dipole Resonance and Neutron Skin in 208Pb from a Bayesian Approach , Chin.Phys.Lett. 38 04210 ( 2021 ) doi: 10.1088/0256-307X/38/4/042101 http://doi.org/10.1088/0256-307X/38/4/042101 Constraining Isovector Nuclear Interactions with Giant Dipole Resonance and Neutron Skin in 208Pb from a Bayesian Approach B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral . Phys. Rev. Lett. 119 , 161101 ( 2017 ). doi: 10.1103/PhysRevLett.119.161101 http://doi.org/10.1103/PhysRevLett.119.161101 GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), GW170817: Measurements of Neutron Star Radii and Equation of State . Phys. Rev. Lett. 121 , 161101 ( 2018 ). doi: 10.1103/PhysRevLett.121.161101 http://doi.org/10.1103/PhysRevLett.121.161101 GW170817: Measurements of Neutron Star Radii and Equation of State E. Annala , T. Gorda , A. Kurkela , et. al. , Gravitational-Wave Constraints on the Neutron-Star-Matter Equation of StateGravitational-Wave Constraints on the Neutron-Star-Matter Equation of State Phys. Rev. Lett. 120 , 172703 ( 2018 ). doi: 10.1103/PhysRevLett.120.172703 http://doi.org/10.1103/PhysRevLett.120.172703 Gravitational-Wave Constraints on the Neutron-Star-Matter Equation of StateGravitational-Wave Constraints on the Neutron-Star-Matter Equation of State F. J. Fattoyev , J. Piekarewicz , and C. J. Horowitz , Neutron Skins and Neutron Stars in the Multimessenger Era , Phys. Rev. Lett. 120 , 172702 ( 2018 ). doi: 10.1103/PhysRevLett.120.172702 http://doi.org/10.1103/PhysRevLett.120.172702 Neutron Skins and Neutron Stars in the Multimessenger Era B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Properties of the Binary Neutron Star Merger GW170817 . Phys. Rev. X 9 , 011001 ( 2019 ). doi: 10.1103/PhysRevX.9.011001 http://doi.org/10.1103/PhysRevX.9.011001 Properties of the Binary Neutron Star Merger GW170817 T. Malik , B. K. Agrawal , J. N. De , et. al. , Tides in merging neutron stars: Consistency of the GW170817 event with experimental data on finite nuclei . Phys. Rev. C 99 , 052801(R) ( 2019 ). doi: 10.1103/PhysRevC.99.052801 http://doi.org/10.1103/PhysRevC.99.052801 Tides in merging neutron stars: Consistency of the GW170817 event with experimental data on finite nuclei N. B. Zhang and B. A. Li , Delineating effects of nuclear symmetry energy on the radii and tidal polarizabilities of neutron stars . J. Phys. G: Nucl. Part. Phys. 46 , 014002 ( 2019 ). doi: 10.1088/1361-6471/aaef54 http://doi.org/10.1088/1361-6471/aaef54 Delineating effects of nuclear symmetry energy on the radii and tidal polarizabilities of neutron stars C. Y. Tsang , M. B. Tsang , P. Danielewicz , et. al. , Insights on Skyrme parameters from GW170817 . Phys. Lett. B 796 , 1 ( 2019 ). doi: 10.1016/j.physletb.2019.05.055 http://doi.org/10.1016/j.physletb.2019.05.055 Insights on Skyrme parameters from GW170817 M. B. Tsang , W. G. Lynch , P. Danielewicz , et. al. , Symmetry energy constraints from GW170817 and laboratory experiments . Phys. Lett. B 795 , 533 ( 2019 ). doi: 10.1016/j.physletb.2019.06.059 http://doi.org/10.1016/j.physletb.2019.06.059 Symmetry energy constraints from GW170817 and laboratory experiments Chong-Ji Jiang , Yu Qiang , Da-Wei Guan , Qing-Zhen Chai , Chun-Yuan Qiao , and Jun-Chen Pei , Chin.Phys.Lett. 38 , 052101 ( 2021 ). http://cpl.iphy.ac.cn/10.1088/0256-307X/38/5/052101 http://cpl.iphy.ac.cn/10.1088/0256-307X/38/5/052101 Chin.Phys.Lett. Jing Zhang , Dehua Wen , Yuxi Li , Constraint on nuclear symmetry energy imposed by f-mode oscillation of neutron stars , Commun. Theor. Phys. 73 ( 2021 ) 115302 . doi: 10.1088/1572-9494/ac1669 http://doi.org/10.1088/1572-9494/ac1669 Constraint on nuclear symmetry energy imposed by f-mode oscillation of neutron stars Bao-Jun Cai , Lie-Wen Chen . Constraints on the skewness coefficient of symmetric nuclear matter within the nonlinear relativistic mean field model . [J]. Nuclear Science and Techniques 28 ( 12 ): 185 ( 2017 ) doi: 10.1007/s41365-017-0329-1 http://doi.org/10.1007/s41365-017-0329-1 . Constraints on the skewness coefficient of symmetric nuclear matter within the nonlinear relativistic mean field model Bao-An Li , Bao-Jun Cai , Lie-Wen Chen , Wen-Jie Xie , Jun Xu , Nai-Bo Zhang , A theoretical overview of isospin and EOS effects in heavy-ion reactions at intermediate energies . arXiv:2201.09133 A theoretical overview of isospin and EOS effects in heavy-ion reactions at intermediate energies B. T. Reed , F. J. Fattoyev , C. J. Horowitz and J. Piekarewicz , Implications of PREX-2 on the Equation of State of Neutron-Rich Matter . Phys. Rev. Lett. 126 , 172503 ( 2021 ). doi: 10.1103/PhysRevLett.126.172503 http://doi.org/10.1103/PhysRevLett.126.172503 Implications of PREX-2 on the Equation of State of Neutron-Rich Matter S. V. Pineda , K. Konig , D. M. Rossi , B. A. Brown , A. Incorvati , J. Lantis , K. Minamisono , W. Nortershauser , J. Piekarewicz and R. Powel , et al. , Charge Radius of Neutron-Deficient 54 Ni and Symmetry Energy Constraints Using the Difference in Mirror Pair Charge Radii , Phys. Rev. Lett. 127 , 182503 ( 2021 ). doi: 10.1103/PhysRevLett.127.182503 http://doi.org/10.1103/PhysRevLett.127.182503 Charge Radius of Neutron-Deficient 54Ni and Symmetry Energy Constraints Using the Difference in Mirror Pair Charge Radii R. Essick , I. Tews , P. Landry and A. Schwenk , Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of 208 Pb with Minimal Modeling Assumptions . Phys. Rev. Lett. 127 , 192701 ( 2021 ). doi: 10.1103/PhysRevLett.127.192701 http://doi.org/10.1103/PhysRevLett.127.192701 Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of 208Pb with Minimal Modeling Assumptions Jun Xu , Lie-Wen Chen , ManYee Betty Tsang , et. al. , Understanding transport simulations of heavy-ion collisions at 100A and 400A MeV: Comparison of heavy-ion transport codes under controlled conditions . Phys. Rev. C 93 044609 ( 2016 ). doi: 10.1103/PhysRevC.93.044609 http://doi.org/10.1103/PhysRevC.93.044609 Understanding transport simulations of heavy-ion collisions at 100A and 400A MeV: Comparison of heavy-ion transport codes under controlled conditions Y. Zhang , Y. Wang , M. Colonna , et. al. , Comparison of heavy-ion transport simulations: Collision integral in a box . Phys. Rev. C 97 034625 ( 2018 ). doi: 10.1103/PhysRevC.97.034625 http://doi.org/10.1103/PhysRevC.97.034625 Comparison of heavy-ion transport simulations: Collision integral in a box A. Ono , J. Xu , M. Colonna , et. al. , Comparison of heavy-ion transport simulations: Collision integral with pions and Δ resonances in a box Phys. Rev. C 100 044617 ( 2019 ). doi: 10.1103/PhysRevC.100.044617 http://doi.org/10.1103/PhysRevC.100.044617 Comparison of heavy-ion transport simulations: Collision integral with pions and Δ resonances in a box M. Colonna , Y. Zhang , Y. Wang et.al. , Comparison of heavy-ion transport simulations: Mean-field dynamics in a box . Phys. Rev. C 104 024603 ( 2021 ). doi: 10.1103/PhysRevC.104.024603 http://doi.org/10.1103/PhysRevC.104.024603 Comparison of heavy-ion transport simulations: Mean-field dynamics in a box Hermann Wolter , Maria Colonna , Dan Cozma , et al. , submitted to Prog. Part. Nucl. Phys. Prog. Part. Nucl. Phys. M. B. Tsang , W. A. Friedman , C. K. Gelbke , et.al. , Isotopic Scaling in Nuclear Reactions . Phys. Rev. Lett. 86 , 5023 ( 2001 ). doi: 10.1103/PhysRevLett.86.5023 http://doi.org/10.1103/PhysRevLett.86.5023 Isotopic Scaling in Nuclear Reactions M. A. Famiano , T. Liu , W. G. Lynch , et. al. , Neutron and Proton Transverse Emission Ratio Measurements and the Density Dependence of the Asymmetry Term of the Nuclear Equation of State Phys. Rev. Lett. 97 , 052701 ( 2006 ). doi: 10.1103/PhysRevLett.97.052701 http://doi.org/10.1103/PhysRevLett.97.052701 Neutron and Proton Transverse Emission Ratio Measurements and the Density Dependence of the Asymmetry Term of the Nuclear Equation of State M. B. Tsang , T. X. Liu , L. Shi , et. al. , Isospin Diffusion and the Nuclear Symmetry Energy in Heavy Ion Reactions . Phys. Rev. Lett. 92 , 6 ( 2004 ). doi: 10.1103/PhysRevLett.92.062701 http://doi.org/10.1103/PhysRevLett.92.062701 Isospin Diffusion and the Nuclear Symmetry Energy in Heavy Ion Reactions Ying-Xun Zhang , Cheng-Shuang Zhou , Ji-Xian Chen , et. al. , Correlation between the fragmentation modes and light charged particles emission in heavy ion collisions Sci. China-Phys. Mech. Astron. Vol. 58 No. 11 ( 2015 ). doi: 10.1007/s11433-015-5723-2 http://doi.org/10.1007/s11433-015-5723-2 Correlation between the fragmentation modes and light charged particles emission in heavy ion collisions Chun-Wang Ma , Yu-Gang Ma , Shannon Information Entropy in Heavy-ion Collisions . Prog. Part. Nucl. Phys. 99 , 120 ( 2018 ). doi: 10.1016/j.ppnp.2018.01.002 http://doi.org/10.1016/j.ppnp.2018.01.002 Shannon Information Entropy in Heavy-ion Collisions Chun-Wang Ma , Hui-Ling Wei , Xing-Quang Liu , et. al. , Nuclear fragments in projectile fragmentation reactions . Prog. Part. Nucl. Phys. 121 , 10391 ( 2021 ). doi: 10.1016/j.ppnp.2021.103911 http://doi.org/10.1016/j.ppnp.2021.103911 Nuclear fragments in projectile fragmentation reactions T. X. Liu , W. G. Lynch , R. H. Showalter , et.al. , Isospin observables from fragment energy spectra . Phys. Rev. C 86 024605 ( 2012 ). doi: 10.1103/PhysRevC.86.024605 http://doi.org/10.1103/PhysRevC.86.024605 Isospin observables from fragment energy spectra Yingxun Zhang , D. D. S. Coupland , P. Danielewicz , et. al. , Influence of in-medium NN cross sections, symmetry potential, and impact parameter on isospin observables . Phys. Rev. C 85 , 024602 ( 2012 ). doi: 10.1103/PhysRevC.85.024602 http://doi.org/10.1103/PhysRevC.85.024602 Influence of in-medium NN cross sections, symmetry potential, and impact parameter on isospin observables Yingxun Zhang and Zhuxia Li , Probing the density dependence of the symmetry potential with peripheral heavy-ion collisions . Phys. Rev. C 71 , 024604 ( 2005 ). doi: 10.1103/PhysRevC.71.024604 http://doi.org/10.1103/PhysRevC.71.024604 Probing the density dependence of the symmetry potential with peripheral heavy-ion collisions Yingxun Zhang and Zhuxia Li , Elliptic flow and system size dependence of transition energies at intermediate energies . Phys. Rev. C 74 , 014602 ( 2006 ). doi: 10.1103/PhysRevC.74.014602 http://doi.org/10.1103/PhysRevC.74.014602 Elliptic flow and system size dependence of transition energies at intermediate energies Yingxun Zhang , P. Danielewicz , M. Famiano , et. al. , The influence of cluster emission and the symmetry energy on neutron-proton spectral double ratios . Phys. Lett. B 664 , 145 ( 2008 ). doi: 10.1016/j.physletb.2008.03.075 http://doi.org/10.1016/j.physletb.2008.03.075 The influence of cluster emission and the symmetry energy on neutron-proton spectral double ratios Ying-Xun Zhang , Ning Wang , Qing-Feng Li , et.al. , Progress of quantum molecular dynamics model and its applications in heavy ion collisions Front . Phys. 15 ( 5 ), 54301 ( 2020 ). doi: 10.1007/s11467-020-0961-9 http://doi.org/10.1007/s11467-020-0961-9 Progress of quantum molecular dynamics model and its applications in heavy ion collisions Front R.J. Charity , M.A. McMahan , G.J. Wozniak , et. al. , systematics of complex fragment emission in niobium-induced reactions , Nucl. Phys. A 483 , 371 ( 1988 ). systematics of complex fragment emission in niobium-induced reactions R. J. Charity , in Joint ICTP-AIEA Advanced Workshop on Model Codes for Spallation Reactions ( IAEA , Vienna , 2008 ), Report INDC(NDC)-0530 . Joint ICTP-AIEA Advanced Workshop on Model Codes for Spallation Reactions Jun Su , Long Zhu , Ching-Yuan Huang , Wen-Jie Xie , et. al. , Correlation between symmetry energy and effective k-mass splitting with an improved isospin- and momentum-dependent interaction . Phys. Rev. C 94 , 034619 ( 2016 ). doi: 10.1103/PhysRevC.94.034619 http://doi.org/10.1103/PhysRevC.94.034619 Correlation between symmetry energy and effective k-mass splitting with an improved isospin- and momentum-dependent interaction Wen-Jie Xie , Jun Su , Long Zhu , et. al. , Neutron-proton effective mass splitting in a Boltzmann-Langevin approach . Phys. Rev. C 88 , 061601(R) ( 2013 ). doi: 10.1103/PhysRevC.88.061601 http://doi.org/10.1103/PhysRevC.88.061601 Neutron-proton effective mass splitting in a Boltzmann-Langevin approach H. H. Gutbrod , A. Sandoval , P. J. Johansen , et.al. , Final-State Interactions in the Production of Hydrogen and Helium Isotopes by Relativistic Heavy Ions on Uranium . Phys. Rev. Lett. 37 , 667 ( 1976 ). doi: 10.1103/PhysRevLett.37.667 http://doi.org/10.1103/PhysRevLett.37.667 Final-State Interactions in the Production of Hydrogen and Helium Isotopes by Relativistic Heavy Ions on Uranium E. Renshaw , S. J. Yennello , K. Kwiatkowski , et. al. , Analyzing powers and isotope ratios for the nat Ag(p → intermediate-mass fragment) reaction at 200 MeV Phys. Rev. C 44 , 2618 ( 1991 ). doi: 10.1103/PhysRevC.44.2618 http://doi.org/10.1103/PhysRevC.44.2618 Analyzing powers and isotope ratios for the natAg(p→ intermediate-mass fragment) reaction at 200 MeV W. Neubert and A. S. Botvina , What is the physics behind the 3 He- 4 He anomaly? Eur. Phys. J. A 7 , 101 ( 2000 ). doi: 10.1007/s100500050016 http://doi.org/10.1007/s100500050016 What is the physics behind the 3He-4He anomaly? Li Li , Yingxun Zhang , Zhuxia Li , et. al. , Impact parameter smearing effects on isospin sensitive observables in heavy ion collisions . Phys. Rev. C 97 , 044606 ( 2018 ). doi: 10.1103/PhysRevC.97.044606 http://doi.org/10.1103/PhysRevC.97.044606 Impact parameter smearing effects on isospin sensitive observables in heavy ion collisions E. Plagnol J. Łukasik , G. Auger , et. al. , Onset of midvelocity emissions in symmetric heavy ion reactions . Phys. Rev. C 61 , 014606 ( 1999 ). doi: 10.1103/PhysRevC.61.014606 http://doi.org/10.1103/PhysRevC.61.014606 Onset of midvelocity emissions in symmetric heavy ion reactions L. Phair , D.R. Bowman , C.K. Gelbke , et. al. , Impact-parameter filters for 36 Ar+ 197 Au collisions at E/A = 50, 80 and 110 MeV . Nuclear Physics A 548 , 489 ( 1992 ). doi: 10.1016/0375-9474(92)90697-I http://doi.org/10.1016/0375-9474(92)90697-I Impact-parameter filters for 36Ar+197Au collisions at E/A = 50, 80 and 110 MeV D. D. S. Coupland , W. G. Lynch , M. B. Tsang , et. al. , Influence of transport variables on isospin transport ratios . Phys. Rev. C 84 , 054603 ( 2011 ). doi: 10.1103/PhysRevC.84.054603 http://doi.org/10.1103/PhysRevC.84.054603 Influence of transport variables on isospin transport ratios Yingxun Zhang , M.B. Tsang , Zhuxia Lia , et. al. , Constraints on nucleon effective mass splitting with heavy ion collisions . Phys. Lett. B 732 , 186 - 190 ( 2014 ). doi: 10.1016/j.physletb.2014.03.030 http://doi.org/10.1016/j.physletb.2014.03.030 Constraints on nucleon effective mass splitting with heavy ion collisions All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Li Li, Fang-Yuan Wang and Ying-Xun Zhang. The first draft of the manuscript was written by Ying-Xun Zhang and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript. Y. Zhang is thankful for the helpful discussions with Prof. M. B. Tsang. 10.1007/s41365-022-01050-w This work was supported by the National Natural Science Foundation of China (Nos. 11875323, 11705163, 11790320, 11790323, and 11961141003), the National Key R&D Program of China (No. 2018 YFA0404404), the Continuous Basic Scientific Research Project (Nos. WDJC-2019-13, BJ20002501), funding from the China Institute of Atomic Energy and the Leading Innovation Project of CNNC (Nos. LC192209000701 and LC202309000201) . 2022-03-29 2022-04-23 2022-04-25 2022-06-15 2022-07-26 © 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. 2022 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-022-01050-w http://dx.doi.org/10.1007/s41365-022-01050-w . 2022 Nuclear Science and Techniques 05 33