Source author record

Yongjia Wang

Yongjia Wang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2023arXiv

Impacts of momentum dependent interaction, symmetry energy and near-threshold $NN\to NΔ$ cross sections on isospin sensitive flow and pion observables

Based on the ultra-relativistic quantum molecular dynamics (UrQMD) model, the impacts of momentum dependent interaction, symmetry energy and near-threshold $NN\to NΔ$ cross sections on isospin sensitive collective flow and pion observables are investigated. Our results confirm that the elliptic flow of neutrons and charged particles, i.e. $v_2^n$ and $v_2^{ch}$, are sensitive to the strength of momentum dependence interaction and the elliptic flow ratio, i.e., $v_2^n/v_2^{ch}$, is sensitive to the stiffness of symmetry energy. For describing the pion multiplicity near the threshold energy, accurate $NN\to NΔ$ cross sections are crucial. With the updated momentum dependent interaction and $NN\to NΔ$ cross sections in UrQMD model, seven observables, such as directed flow and elliptic flow of neutrons and charged particles, the elliptic flow ratio of neutrons to charged particles, charged pion multiplicity and its ratio $π^-/π^+$, can be well described by the parameter sets with the slope of symmetry energy from 5 MeV to 70 MeV. To describe the constraints of symmetry energy at the densities probed by the collective flow and pion observables, the named characteristic density is investigated and used. Our analysis found that the flow characteristic density is around 1.2$ρ_0$ and pion characteristic density is around 1.5$ρ_0$, and we got the constrains of symmetry energy at characteristic densities are $S(1.2ρ_0)=34\pm 4$ MeV and $S(1.5ρ_0)=36\pm 8$ MeV. These results are consistent with previous analysis by using pion and flow observable with different transport models, and demonstrate a reasonable description of symmetry energy constraint should be presented at the characteristic density of isospin sensitive observables.

preprint2022arXiv

Accessing the in-medium effects on nucleon-nucleon elastic cross section with collective flows and nuclear stopping

A systematic study of the in-medium correction factor ($F$) on nucleon-nucleon elastic cross section is performed within the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model. The effects of the beam energy dependence of $F$ on the directed, elliptic flow and nuclear stopping in $^{197}$Au+$^{197}$Au collisions with energy ranging from $0.09$ to $0.8A$ GeV are explored. It is found that the directed, elliptic flow and nuclear stopping at relatively low energies are very sensitive to $F$, and the sensitivity gradually weakens with increasing beam energy. The beam energy dependent in-medium correction factor $F$ is deduced from the comparison of the excitation functions of the directed, elliptic flow and nuclear stopping between the calculated results and the FOPI experimental data.

preprint2022arXiv

Decoding the nuclear symmetry energy event-by-event in heavy-ion collisions with machine learning

Inferences of the nuclear symmetry energy from heavy-ion collisions are currently based on the comparison of measured observables and transport model simulations. Only the expectation values of observables over all considered events are used in these approaches, however, observables can be obtained event-by-event both in experiments and transport model simulations. By using the light gradient boosting machine (LightGBM), a modern machine-learning algorithm, we present a framework for inferring the density-dependent nuclear symmetry energy from observables in heavy-ion collisions on the event-by-event analysis. The ultrarelativistic quantum molecular dynamics (UrQMD) model simulations are used as training data. The symmetry energy slope parameter extracted with LightGBM event-by-event from test data also by UrQMD has an average spread of approximately 30~MeV from the truth, and is found to be robust against variations in model parameters. In addition, LightGBM can identify features that have the greatest effect on the physics of interest, thereby offering valuable insights. Our study suggests that the present framework can be a powerful tool and may offer a new paradigm to study the underlying physics in heavy-ion collisions.

preprint2022arXiv

Elliptic flow in heavy-ion collisions at intermediate energy: the role of impact parameter, mean field potential, and collision term

Within the ultrarelativistic quantum molecular dynamics (UrQMD) model, by reverse tracing nucleons that are finally emitted at mid-rapidity (|$y_0$| < 0.1) in the entire reaction process, the time evolution of elliptic flow ($v_2$) of these traced nucleons produced in Au+Au collisions at beam energy of 0.4 GeV$/$nucleon with different impact parameters ($b$) is studied. The initial value of $v_2$ is positive and increases with $b$, then it decreases as time passes and tends to saturate at a negative value. It is found that nucleon-nucleon collisions always depress the value of $v_2$ (enhance the out-of-plane emission), while the nuclear mean field potential may slightly raise the value of $v_2$ during the expansion stage in peripheral reactions. The related density mostly probed by $v_2$ of nucleons at mid-rapidity is found to be $\sim$ 60% of the maximum density reached during the collisions.

preprint2022arXiv

How Do Constraints of Nuclear Symmetry Energy Reconcile with Different Models?

By simultaneously describing the data of isospin sensitive nucleonic flow and pion observables, such as $v_2^n/v_2^{ch}$ and $π^-/π^+$, with ultra-relativistic quantum molecular dynamics (UrQMD) model, we got the symmetry energy at flow and pion characteristic densities which are $S(1.2ρ_0)=34\pm 4$ MeV and $S(1.5ρ_0)=36\pm 8$ MeV. Within the uncertainties, the constraints of symmetry energy at characteristic densities are consistent with the previous constraints by using other transport models. The consistency suggests that the reliable constraints on symmetry energy should be presented at the characteristic density of isospin sensitive observables. By using the constraints of symmetry energy at two different characteristic densities, the extrapolated value of $L$ is provided. Within $2σ$ uncertainty, the extrapolated value of $L$ is in $5-70$ MeV which is consistent with the recent combination analysis from PREX-II and astrophyiscs data. Further, the calculations with the constrained parameter sets can describe the data of charged pion multiplicities from S$π$RIT collaboration.

preprint2021arXiv

Proton correlations and apparent intermittency in the UrQMD model with hadronic potentials

It is shown that the inclusion of hadronic interactions, and in particular nuclear potentials, in simulations of heavy ion collisions at the SPS energy range can lead to obvious correlations of protons. These correlations contribute significantly to an intermittency analysis as performed at the NA61 experiment. The beam energy and system size dependence is studied by comparing the resulting intermittency index for heavy ion collisions of different nuclei at beam energies of $40A$, $80A$ and $150A$ GeV. The resulting intermittency index from our simulations is similar to the reported values of the NA61 collaboration, if nuclear interactions are included. The observed apparent intermittency signal is the result of the correlated proton pairs with small relative transverse momentum $Δp_{t}$, which would be enhanced by hadronic potentials, and this correlation between the protons is slightly influenced by the coalescence parameters and the relative invariant four-momentum $q_{inv}$ cut.

preprint2020arXiv

Application of microscopic transport model in the study of nuclear equation of state from heavy ion collisions at intermediate energies

The equation of state (EOS) of nuclear matter, i.e., the thermodynamic relationship between the binding energy per nucleon, temperature, density, as well as the isospin asymmetry, has been a hot topic in nuclear physics and astrophysics for a long time. The knowledge of the nuclear EOS is essential for studying the properties of nuclei, the structure of neutron stars, the dynamics of heavy ion collision (HIC), as well as neutron star mergers. HIC offers a unique way to create nuclear matter with high density and isospin asymmetry in terrestrial laboratory, but the formed dense nuclear matter exists only for a very short period, one cannot measure the nuclear EOS directly in experiments. Practically, transport models which often incorporate phenomenological potentials as an input are utilized to deduce the EOS from the comparison with the observables measured in laboratory. The ultrarelativistic quantum molecular dynamics (UrQMD) model has been widely employed for investigating HIC from the Fermi energy (40 MeV per nucleon) up to the CERN Large Hadron Collider energies (TeV). With further improvement in the nuclear mean-field potential term, the collision term, and the cluster recognition term of the UrQMD model, the newly measured collective flow and nuclear stopping data of light charged particles by the FOPI Collaboration can be reproduced. In this article we highlight our recent results on the studies of the nuclear EOS and the nuclear symmetry energy with the UrQMD model. New opportunities and challenges in the extraction of the nuclear EOS from transport models and HIC experiments are discussed.

preprint2020arXiv

Beam energy dependence of cumulants of the net-baryon, net-charge and deuteron multiplicity distributions in Au+Au collisions at $\sqrt{s_{NN}}=3.0-5.0$ GeV

Within the ultra-relativistic quantum molecular dynamics (UrQMD) model, in which the Lorentz-covariant treatment of nuclear mean-field potential is considered, the fluctuations of net-baryon, net-charge and deuterons multiplicity distributions in Au+Au head-on collisions at $\sqrt{s_{NN}}=3.0-5.0$ GeV are calculated. The results show that the nuclear mean-field potential can significantly enhance the magnitude of baryon number fluctuations in narrow rapidity windows, and this enhancement rapidly weakens with increasing beam energy. However, for proton and net-charge number fluctuations, the mean-field effects are less noticeable than that for baryon number. In addition, for net-charge number fluctuations, the negative binomial distribution agrees well with the calculated results at mid-pseudorapidity window. Finally, the event-by-event fluctuations of deuteron number in the coalescence production picture are calculated as well, it is found that its cumulant ratios decrease linearly with increasing the average multiplicity of deuterons per event, i.e., increase with increasing beam energy.

preprint2020arXiv

Elliptic flow splitting between protons and antiprotons from hadronic potentials

The difference in elliptic flow $v_{2}$ between protons and antiprotons, produced in $^{197}\text{Au}+^{197}\text{Au}$ collisions at center-of-mass energies $\sqrt{s_{NN}}=5-12~\text{GeV}$, is studied within a modified version of the ultrarelativistic quantum molecular dynamics (UrQMD) model. Two different model scenarios are compared: the cascade mode and the mean field mode which includes potential interactions for both formed and pre-formed hadrons. The model results for the elliptic flow of protons and the relative $v_{2}$ difference between protons and antiprotons obtained from the mean field mode agree with the available experimental data, while the $v_{2}$ difference is near zero for the cascade mode. Our results show that the elliptic flow splitting, observed for particles and antiparticles, can be explained by the inclusion of proper hadronic interactions. In addition, the difference in $v_{2}$ between protons and antiprotons depends on the centrality and the rapidity window. With smaller centrality and/or rapidity acceptance, the observed elliptic flow splitting is more sensitive to the beam energy, indicating a strong net baryon density dependence of the effect. We propose to confirm this splitting at the upcoming experiments from Beam Energy Scan (BES) Phase-\Rmnum{2} at Relativistic Heavy Ion Collider (RHIC), the Compressed Baryonic Matter (CBM) at Facility for Antiproton and Ion Research (FAIR), High Intensity heavy ion Accelerator Facility (HIAF) and Nuclotron-based Ion Collider fAcility (NICA).

preprint2020arXiv

Insights on pion production mechanism and symmetry energy at high density

The $NΔ\to NN$ cross sections, which take into account the $Δ$-mass dependence of M-matrix and momentum $p_{NΔ}$, are applied on the calculation of pion production within the framework of the UrQMD model. Our study shows that UrQMD calculations with the $Δ$-mass dependent $NΔ\to NN$ cross sections enhance the pion multiplicities and decrease the $π^-/π^+$ ratios. By analyzing the time evolution of the pion production rate and the density in the overlapped region for Au+Au at the beam energy of 0.4A GeV, we find that the pion multiplicity probes the symmetry energy in the region of 1-2 times normal density. The process of pion production in the reaction is tracked including the loops of $NN\leftrightarrow NΔ$ and $Δ\leftrightarrow Nπ$, our calculations show that the sensitivity of $π^-/π^+$ to symmetry energy is weakened after 4-5 N-$Δ$-$π$ loops in the pion production path, while the $π^{-}/π^{+}$ ratio in reactions at near threshold energies remains its sensitivity to the symmetry energy. By comparing the calculations to the FOPI data, we obtain a model dependent conclusion on the symmetry energy and the symmetry energy at two times normal density is $S(2ρ_0)$=38-73 MeV within $1σ$ uncertainties. Under the constraints of tidal deformability and maximum mass of neutron star, the symmetry energy at two times normal density is reduced to $48-58$ MeV and slope of symmetry energy $L=54-81$ MeV, and it is consistent with the constraints from ASY-EOS flow data.

preprint2016arXiv

An explanation of the elliptic flow difference between proton and anti-proton from the UrQMD model with hadron potentials

The time evolution of both proton and anti-proton $v_2$ flows from Au+Au collisions at $\sqrt{s_{NN}}$=7.7 GeV are examined by using both pure cascade and mean-field potential versions of the UrQMD model. Due to a stronger repulsion at the early stage introduced by the repulsive potentials and hence much less annihilation probabilities, anti-protons are frozen out earlier with smaller $v_2$ values. Therefore, the experimental data of anti-proton $v_2$ as well as the flow difference between proton and anti-proton can be reasonably described with the potential version of UrQMD.

preprint2016arXiv

Helium-3 production from Pb+Pb collisions at SPS energies with the UrQMD model and the traditional coalescence afterburner

A potential version of the UrQMD (UrQMD/M) transport model and a traditional coalescence model are combined to calculate the production of $^3$He fragments in central Pb+Pb collisions at SPS energies 20-80 GeV/nucleon. It is found that the Lorentz transformation in the afterburner influences visibly the $^3$He yield and should be considered in calculations. The rapidity distribution of $^3$He multiplicities (including the concave shape) can be described well with UrQMD/M when it stops during t$_{\rm cut}$=100$\pm$25 fm$/c$ and the coalescence afterburner with one parameter set of ($R_0$,$P_0$)=(3.8 fm, 0.3 GeV$/$c) is taken into use afterwards.

preprint2016arXiv

Influence of coalescence parameters on the production of protons and Helium-3 fragments

The time evolution of protons and $^3$He fragments from Au+Au/Pb+Pb reactions at 0.25, 2, and 20 GeV$/$nucleon is investigated with the potential version of the Ultrarelativistic Quantum Molecular Dynamics (UrQMD) model combined with the traditional coalescence afterburner. In the coalescence process, the relative distance $R_0$ and relative momentum $P_0$ are surveyed in the range of 3-4 fm and 0.25-0.35 GeV$/$c, respectively. For both clusters, a strong reversed correlation between $R_0$ and $P_0$ is seen and it is time-dependent as well. For protons, the accepted ($R_0$, $P_0$) bands lie in the time interval 30-60 fm$/$c, while for $^3$He, a longer time evolution (at about 60-90 fm$/$c) is needed. Otherwise, much smaller $R_0$ and $P_0$ values should be chosen. If we further look at the rapidity distributions from both central and semi-central collisions, it is found that the accepted [$t_{\rm cut}, (R_0, P_0$)] assemble can provide consistent results for proton yield and collective flows especially at mid-rapdities, while for $^3$He, the consistency is destroyed at both middle and projectile-target rapidities.

preprint2016arXiv

Influence of differential elastic nucleon-nucleon cross section on stopping and collective flows in heavy-ion collisions at intermediate energies

By considering three different Nucleon-Nucleon (NN) elastic differential cross sections: the Cugnon \emph{et al.} parameterized differential cross section [Nucl. Instrum. Methods Phys. Res., Sect. \textbf{B111}, 215 (1996)], and the differential cross section derived from the collision term of the self-consistent relativistic Boltzmann-Uehling-Uhlenbeck equation proposed by Mao \emph{et al.} [Z.\ Phys.\ A {\bf 347}, 173 (1994)], as well as the isotropic differential cross section, within the newly updated version of the ultrarelativistic quantum molecular dynamics (UrQMD) model, the influence of the differential elastic NN cross section on various observables (e.g., nuclear stopping, both the rapidity and transverse-velocity dependence of the directed and elliptic flows) in Au+Au collisions at beam energies 150, 250, 400, and 800 MeV$/$nucleon is investigated. By comparing calculations with those three differential cross sections, it is found that the nuclear stopping power, the directed and elliptic flows are affected to some extent by the differential cross sections, and the impact of differential cross section on those observables becomes more visible as the beam energy increases. The effect on the elliptic flow difference $v_{2}^{n}$-$v_{2}^{H}$ and ratio $v_{2}^{n}$/$v_{2}^{H}$ of neutrons versus hydrogen isotopes ($Z=1$), which have been used as sensitive observables for probing nuclear symmetry energy at high densities, is weak.

preprint2016arXiv

Rapidity distribution of protons from the potential version of UrQMD model and the traditional coalescence afterburner

Rapidity distributions of both E895 proton data at AGS energies and NA49 net proton data at SPS energies can be described reasonably well with a potential version of the UrQMD in which mean-field potentials for both pre-formed hadrons and confined baryons are considered, with the help of a traditional coalescence afterburner in which one parameter set for both relative distance $R_0$ and relative momentum $P_0$, (3.8 fm, 0.3 GeV$/$c), is used. Because of the large cancellation between the expansion in $R_0$ and the shrinkage in $P_0$ through the Lorentz transformation, the relativistic effect in clusters has little effect on the rapidity distribution of free (net) protons. Using a Woods-Saxon-like function instead of a pure logarithmic function as seen by FOPI collaboration at SIS energies, one can fit well both the data at SIS energies and the UrQMD calculation results at AGS and SPS energies. Further, it is found that for central Au+Au or Pb+Pb collisions at top SIS, SPS and RHIC energies, the proton fractions in clusters are about 33$\%$, 10$\%$, and 0.7$\%$, respectively.

preprint2016arXiv

Results of the ASY-EOS experiment at GSI: The symmetry energy at suprasaturation density

Directed and elliptic flows of neutrons and light charged particles were measured for the reaction 197Au+197Au at 400 MeV/nucleon incident energy within the ASY-EOS experimental campaign at the GSI laboratory. The detection system consisted of the Large Area Neutron Detector LAND, combined with parts of the CHIMERA multidetector, of the ALADIN Time-of-flight Wall, and of the Washington-University Microball detector. The latter three arrays were used for the event characterization and reaction-plane reconstruction. In addition, an array of triple telescopes, KRATTA, was used for complementary measurements of the isotopic composition and flows of light charged particles. From the comparison of the elliptic flow ratio of neutrons with respect to charged particles with UrQMD predictions, a value γ= 0.72 \pm 0.19 is obtained for the power-law coefficient describing the density dependence of the potential part in the parametrization of the symmetry energy. It represents a new and more stringent constraint for the regime of supra-saturation density and confirms, with a considerably smaller uncertainty, the moderately soft to linear density dependence deduced from the earlier FOPI-LAND data. The densities probed are shown to reach beyond twice saturation.

preprint2016arXiv

The production of unknown neutron-rich isotopes in $^{238}$U+$^{238}$U collisions at near-barrier energy

The production cross sections for primary and residual fragments with charge number from $Z$=70 to 120 produced in the collision of $^{238}$U+$^{238}$U at 7.0 MeV/nucleon are calculated by the improved quantum molecular dynamics (ImQMD) model incorporated with the statistical evaporation model (HIVAP code). The calculation results predict that about sixty unknown neutron-rich isotopes from element Ra ($Z$=88) to Db ($Z$=105) can be produced with the production cross sections above the lower bound of $10^{-8}$ mb in this reaction. And almost all of unknown neutron-rich isotopes are emitted at the laboratory angles $θ_{lab}\leq$ 60$^\circ$. Two cases, i.e. the production of the unknown uranium isotopes with $A\geq$ 244 and that of rutherfordium with $A\geq$ 269 are investigated for understanding the production mechanism of unknown neutron-rich isotopes. It is found that for the former case the collision time between two uranium nuclei is shorter and the primary fragments producing the residues have smaller excitation energies of $\leq$ 30 MeV and the outgoing angles of those residues cover a range of 30$^\circ$-60$^\circ$. For the later case, the longer collision time is needed for a large number of nucleons being transferred and thus it results in the higher excitation energies and smaller outgoing angles of primary fragments, and eventually results in a very small production cross section for the residues of Rf with $A\geq$ 269 which have a small interval of outgoing angles of $θ_{lab}$=40$^\circ$-50$^\circ$.

preprint2015arXiv

Influence of clustering and hadron potentials on the rapidity distribution of protons from the UrQMD model

The Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model supplemented by potentials for both pre-formed hadrons and confined baryons (called UrQMD/M) are used to describe rapidity distributions of both the E895 proton data at AGS energies and the NA49 net proton data at SPS energies. With the help of a coalescence afterburner using only one parameter set of ($R_0$, $P_0$)=(3.8 fm, 0.3 GeV$/$c), both sets of experimental data can be described fairly well except for a small discrepancy seen for the net protons at mid-rapidity from heavy ion collisions (HICs) at high SPS energies. Furthermore, in contrast to the logarithmic dependence with beam energy at SIS energies there are still about 10$\%$ of protons in clusters from central HICs at the beam energy of $80$ GeV$/$nucleon.

preprint2015arXiv

Residue cross sections of $^{50}$Ti-induced fusion reactions based on the two-step model

$^{50}$Ti-induced fusion reactions to synthesize superheavy elements are studied systematically with the two-step model developed recently, where fusion process is divided into approaching phase and formation phase. Furthermore, the residue cross sections for different neutron evaporation channels are evaluated with the statistical evaporation model. In general, the calculated cross sections are much smaller than that of $^{48}$Ca-induced fusion reactions, but the results are within the detection capability of experimental facilities nowadays. The maximum calculated residue cross section for producing superheavy element $Z=119$ is in the reaction $^{50}$Ti+$^{247}$Bk in $3n$ channels with $σ_{\rm res}(3n)=0.043$ pb at $E^{*}$ = 37.0 MeV.

preprint2014arXiv

$^3$H/$^3$He ratio as a probe of the nuclear symmetry energy at sub-saturation densities

Within the newly updated version of the Ultra-relativistic quantum molecular dynamics (UrQMD) model in which the Skyrme potential energy-density functional is introduced, the yield ratio between $^3$H and $^3$He clusters emitted from central $^{40}$Ca+$^{40}$Ca, $^{96}$Zr+$^{96}$Zr, $^{96}$Ru+$^{96}$Ru, and $^{197}$Au+$^{197}$Au collisions in the beam energy range from 0.12 to 1 GeV$/$nucleon is studied. The recent FOPI data for the $^3$H$/$$^3$He ratio are compared with UrQMD calculations using 13 Skyrme interactions (all exhibiting similar values of iso-scalar incompressibility but very different density dependences of the symmetry energy). It is found that the $^3$H$/$$^3$He ratio is sensitive to the nuclear symmetry energy at sub-saturation densities. Model calculations with moderately soft to linear symmetry energies are in agreement with the experimental FOPI data. This result is in line with both, the recent constraints on the low-density symmetry energy available in the literature and our previous results for the high-density symmetry energy obtained with the elliptic flow of free nucleons and hydrogen isotopes as a sensitive probe.

preprint2014arXiv

Collective flows of light particles in the Au+Au collision at intermediate energies

The Skyrme potential energy density functional is introduced into the Ultrarelativistic Quantum Molecular Dynamics (UrQMD) model and the updated version is applied to studying the directed and elliptic flows of light particles (protons, neutrons, deuterons, tritons, $^3$He and $^4$He) in $^{197}$Au+$^{197}$Au collisions at beam energies 150, 250 and 400 MeV/nucleon. The results are compared with the recent FOPI experimental data. It is found that the yields and collective flows of light particles can be described quite well. The influence of the equation of state (EoS), medium-modified nucleon-nucleon elastic cross sections (NNECS) and cluster recognition criteria on the directed and elliptic flows is studied in detail. It is found that the flows of light particles are sensitive to the medium-modified NNECS, but not sensitive to the isospin dependent cluster recognition criteria. It seems difficult, however, even with the new data and calculations, to obtain a more accurate constraint on the nuclear incompressibility $K_0$ than the interval 200-260 MeV.

preprint2014arXiv

Constraining the high-density nuclear symmetry energy with the transverse-momentum dependent elliptic flow

Within the newly updated version of the ultrarelativistic quantum molecular dynamics (UrQMD) model, the transverse-velocity dependence of the elliptic flow of free nucleons from $^{197}$Au+$^{197}$Au collisions at the incident energy 400 MeV$/$nucleon is studied within different windows of the normalized c.m. rapidity $y_0$. It is found that the elliptic flow difference $v_{2}^{n}$-$v_{2}^{p}$ and ratio $v_{2}^{n}$/$v_{2}^{p}$ of neutrons versus protons are sensitive to the density dependence of the symmetry energy, especially the ratio $v_{2}^{n}$/$v_{2}^{p}$ at small transverse velocity in the intermediate rapidity intervals $0.4<|y_0|<0.6$. By comparing either transverse-momentum dependent or integrated FOPI/LAND elliptic flow data of nucleons and hydrogen isotopes with calculations using various Skyrme interactions, all exhibiting similar values of isoscalar incompressibility but very different density dependences of the symmetry energy, a moderately soft to linear symmetry energy is extracted, in good agreement with previous UrQMD or Tübingen QMD model calculations but contrasting results obtained with $π^-/π^+$ yield ratios available in the literature.

preprint2014arXiv

Effect of the spin-orbit interaction on flows in heavy-ion collisions at intermediate energies

The effect of the spin-orbit coupling in heavy ion collisions is investigated based on an updated version of the ultra-relativistic quantum molecular dynamics (UrQMD) model, in which the Skyrme potential energy density functional is employed. And in special, the spin-orbit coupling effects on the directed and elliptic flows of free nucleons emitted from $^{197}$Au+$^{197}$Au collisions as functions of both the beam energy and the impact parameter are studied. Our results show that the net contribution of the spin-orbit term to flows of nucleons is negligible, whereas a directed flow splitting between spin-up and spin-down nucleons is visible especially at large impact parameters and a peak of the splitting is found at the beam energy around 150 MeV$/$nucleon. We also found that the directed flow splitting between spin-up and spin-down neutrons is comparable with the neutron directed flow difference calculated by a soft and a stiff symmetry energy, indicating that the directed flow of neutrons cannot be used to pin down the stiffness of symmetry energy any more without considering the spin degree of freedom in models in case of spin polarization.

preprint2014arXiv

Normal or abnormal isospin-fractionation as a qualitative probe of nuclear symmetry energy at supradensities

Within two different frameworks of isospin-dependent transport model, effect of nuclear symmetry energy at supradensities on the isospin-fractionation (IsoF) was investigated. With positive/negative symmetry potential at supradensities (i.e., values of symmetry energy increase/decrease with density above saturation density), for energetic nucleons, the value of neutron to proton ratio of free nucleons is larger/smaller than that of bound nucleon fragments. Compared with extensively studied quantitative observables of nuclear symmetry energy, the normal or abnormal isospin-fractionation of energetic nucleons can be a qualitative probe of nuclear symmetry energy at supradensities.

preprint2013arXiv

Model dependence of isospin sensitive observables at high densities

Within two different frameworks of isospin-dependent transport model, i.e., Boltzmann-Uehling-Uhlenbeck (IBUU04) and Ultrarelativistic Quantum Molecular Dynamics (UrQMD) transport models, sensitive probes of nuclear symmetry energy are simulated and compared. It is shown that neutron to proton ratio of free nucleons, pi-/pi+ ratio as well as isospin-sensitive transverse and elliptic flows given by the two transport models with their "best settings", all have obvious differences. Discrepancy of numerical value of isospin-sensitive n/p ratio of free nucleon from the two models mainly originates from different symmetry potentials used and discrepancies of numerical value of charged pi-/pi+ ratio and isospin-sensitive flows mainly originate from different isospin-dependent nucleon-nucleon cross sections. These demonstrations call for more detailed studies on the model inputs (i.e., the density- and momentum-dependent symmetry potential and the isospin-dependent nucleon-nucleon cross section in medium) of isospin-dependent transport model used. The studies of model dependence of isospin sensitive observables can help nuclear physicists to pin down the density dependence of nuclear symmetry energy through comparison between experiments and theoretical simulations scientifically.

preprint2011arXiv

Non-equilibrium dynamics in heavy ion collisions at low SIS energies

The Ultrarelativistic Quantum Molecular Dynamics (UrQMD) model, a microscopic transport model, is used to study the directed and elliptic collective flows and the nuclear stopping in Au+Au collisions at incident energies covered by INDRA and lower-energy FOPI experiments. It is seen clearly that these observables are sensitive to both, the potential terms (including iso-scalar and iso-vector parts as well as the momentum dependent term) in the equation of state (EoS) and the collision term (including the Pauli-blocking and the medium-modified nucleon-nucleon elastic cross section (NNECS)). The momentum modifications of both, the mean-field potentials and the density dependent NNECS, are found to be sensitive to the collectivity of heavy-ion collisions. At INDRA energies ($\le 150$~MeV/nucleon), the dynamic transport with a soft EoS with momentum dependence and with the momentum-modified density-dependent NNECS describes the directed flow exhibited by hydrogen isotopes ($Z=1$) emitted at mid-rapidity fairly well.