Source author record

Qingfeng Li

Qingfeng Li 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

41works
8topics
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

41 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.

preprint2023arXiv

Machine learning in nuclear physics at low and intermediate energies

Machine learning is becoming a new paradigm for scientific research in various research fields due to its exciting and powerful capability of modeling tools used for big-data processing task. In this mini-review, we first briefly introduce different methodologies of the machine learning algorithms and techniques. As a snapshot of many applications by machine learning, some selected applications are presented especially for low and intermediate energy nuclear physics, which include topics on theoretical applications in nuclear structure, nuclear reactions, properties of nuclear matter as well as experimental applications in event identification/reconstruction, complex system control and firmware performance. Finally, we also give a brief summary and outlook on the possible directions of using machine learning in low-intermediate energy nuclear physics and possible improvements in ML algorithms.

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

An Experiment to Test the Mechanical Losses of Different Bonding Techniques in Fused Silica

High-purity glasses are used for their low optical and mechanical loss, which makes them an excellent material for oscillators in optical systems, such as inertial sensors. Complex geometries often require the assembly of multiple pieces of glass and their permanent bonding. One common method is hydroxide catalysis bonding, which leaves an enclosed medium layer. This layer has different mechanical properties to the bulk glass around it. The higher mechanical loss of this layer makes it more susceptible to displacement noise originating from the conversion of energy from oscillation to heat and vice versa. Therefore, other methods are needed to bond together glass assemblies. To investigate this, we have set up an experiment to measure the mechanical losses of several different types of bond commonly used in fused silica manufacturing, namely; plasma activated direct bonding, hydroxide catalysis bonding, laser welding, and adhesive bonding. In this paper we present the experimental design and show initial results of the first test sample.

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.

preprint2020arXiv

Progress of Quantum Molecular Dynamics model and its applications in Heavy Ion Collisions

In this review article, we first briefly introduce the transport theory and quantum molecular dynamics model applied in the study of the heavy ion collisions from low to intermediate energies. The developments of improved quantum molecular dynamics model (ImQMD) and ultra-relativistic quantum molecular dynamics model (UrQMD), are reviewed. The reaction mechanism and phenomena related to the fusion, multinucleon transfer, fragmentation, collective flow and particle production are reviewed and discussed within the framework of the two models. The constraints on the isospin asymmetric nuclear equation of state and in-medium nucleon-nucleon cross sections by comparing the heavy ion collision data with transport models calculations in last decades are also discussed, and the uncertainties of these constraints are analyzed as well. Finally, we discuss the future direction of the development of the transport models for improving the understanding of the reaction mechanism, the descriptions of various observables, the constraint on the nuclear equation of state, as well as for the constraint on in-medium nucleon-nucleon cross sections.

preprint2019arXiv

DS-GCNs: Connectome Classification Using Dynamic Spectral Graph Convolution Networks with Assistant Task Training

Functional Connectivity (FC) matrices measure the regional interactions in the brain and have been widely used in neurological brain disease classification. However, a FC matrix is neither a natural image which contains shape and texture information, nor a vector of independent features, which renders the extracting of efficient features from matrices as a challenging problem. A brain network, also named as connectome, could forma a graph structure naturally, the nodes of which are brain regions and the edges are interregional connectivity. Thus, in this study, we proposed novel graph convolutional networks (GCNs) to extract efficient disease-related features from FC matrices. Considering the time-dependent nature of brain activity, we computed dynamic FC matrices with sliding-windows and implemented a graph convolution based LSTM (long short term memory) layer to process dynamic graphs. Moreover, the demographics of patients were also used to guide the classification. However, unlike in conventional methods where personal information, i.e., gender and age were added as extra inputs, we argue that this kind of approach may not actually improve the classification performance, for such personal information given in dataset was usually balanced distributed. In this paper, we proposed to utilize the demographic information as extra outputs and to share parameters among three networks predicting subject status, gender and age, which serve as assistant tasks. We tested the performance of the proposed architecture in ADNI II dataset to classify Alzheimer's disease patients from normal controls. The classification accuracy, sensitivity and specificity reach 0.90, 0.92 and 0.89 on ADNI II dataset.

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

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.

preprint2012arXiv

Examination of scaling of Hanbury-Brown--Twiss radii with charged particle multiplicity

In the light of the recent LHC data on proton-proton and lead-lead collisions we examine the question of the multiplicity scaling of HBT radii in relativistic nuclei and particle interactions. Within the UrQMD transport approach we study a large variety of system sizes at different beam energies and extract the HBT radii. In the calculation, we find a good scaling of the radii as a function of charged particle multiplicity, if the change in the multiplicity is caused by a change of centrality at the same energy. However, the scaling is only approximate when the energy, $\sqrt{s}$, is changed and breaks down when comparing pp to AA reactions.

preprint2012arXiv

Formation time dependence of femtoscopic $ππ$ correlations in p+p collisions at $\sqrt{s_{NN}}$=7 TeV

We investigate femtoscopic $ππ$ correlations using the UrQMD approach combined with a correlation afterburner. The dependence of $ππ$ correlations on the charged particle multiplicity and formation time in p+p collisions at $\sqrt{s_{NN}}$= 7 TeV is explored and compared to present ALICE data. The data allows to constrain the formation time in the string fragmentation to $τ_f \leq 0.8$ fm/c.

preprint2012arXiv

UrQMD calculations of two-pion HBT correlations in central Pb-Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV

Two-pion Hanbury-Brown-Twiss (HBT) correlations for central Pb-Pb collisions at the Large-Hadron-Collider (LHC) energy of $\sqrt{s_{NN}}=2.76$ TeV are investigated for the first time with the microscopic transport model UrQMD (Ultra-relativistic Quantum Molecular Dynamics). The transverse momentum dependence of the Pratt-Bertsch HBT radii is extracted from a three dimensional Gaussian fit to the correlator in the longitudinal co-moving system (LCMS). Qualitative agreement with the ALICE data is obtained, however $R_{out}$ is overpredicted by nearly 50%. The LHC results are also compared to data from the STAR experiment at RHIC. For both energies we find that the calculated $R_O/R_S$ ratio is always larger than data, indicating that the emission in the model is less explosive than observed in the data.

preprint2012arXiv

UrQMD calculations of two-pion HBT correlations in p+p and Pb+Pb collisions at LHC energies

Two-pion Hanbury-Brown-Twiss (HBT) correlations for p+p and central Pb+Pb collisions at the Large-Hadron-Collider (LHC) energies are investigated with the ultra-relativistic quantum molecular dynamics model combined with a correlation afterburner. The transverse momentum dependence of the Pratt-Bertsch HBT radii $R_{long}$, $R_{out}$, and $R_{side}$ is extracted from a three-dimensional Gaussian fit to the correlator in the longitudinal co-moving system. In the p+p case, the dependence of correlations on the charged particle multiplicity and formation time is explored and the data allows to constrain the formation time in the string fragmentation to $τ_f \leq 0.8$ fm/c. In the Pb+Pb case, it is found that $R_{out}$ is overpredicted by nearly 50%. The LHC results are also compared to data from the STAR experiment at RHIC. For both energies we find that the calculated $R_{out}/R_{side}$ ratio is always larger than data, indicating that the emission in the model is less explosive than observed in the data.

preprint2011arXiv

Fusion Hindrance in the Heavy Ion Reactions -- Border Between the Normal and Hindered Fusions

The fusion hindrance in heavy ion collisions is studied in the framework of the two-center liquid drop model. It appears that the neck and the radial degrees of freedom might both be hampered by an inner potential barrier on their path between the contact configuration to the compound nucleus. Heavy ion reactions with and without the two kinds of fusion hindrance are classified through systematic calculations. It is found that the number of reactions without radial fusion hindrance is much smaller than that without neck fusion hindrance, and for both kinds of fusion hindrance the number of reactions without fusion hindrance at small mass-asymmetry parameter $α$ is smaller than that at large $α$. In the formation of a given compound nucleus, if a reaction with $α_c$ is not hindered, then other reactions with $α> α_c$ are also not hindered as it is well known experimentally.

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.

preprint2010arXiv

Production and rescattering of strange baryons at SPS energies in a transport model with hadron potentials

A mean-field potential version of the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model is used to investigate the production of strange baryons, especially the $Λ$s and $\overlineΛ$s, from heavy ion collisions at SPS energies. It is found that, with the consideration of both formed and pre-formed hadron potentials in UrQMD, the transverse mass and longitudinal rapidity distributions of experimental data of both $Λ$s and $\overlineΛ$s can be quantitatively explained fairly well. Our investigation also shows that both the production mechanism and the rescattering process of hadrons play important roles in the final yield of strange baryons.

preprint2010arXiv

System-size dependence of the pion freeze-out volume as a potential signature for the phase transition to a Quark Gluon Plasma

Hanburry-Brown-Twiss (HBT) correlation functions and radii of negatively charged pions from C+C, Si+Si, Cu+Cu, and In+In at lower RHIC/SPS energies are calculated with the UrQMD transport model and the CRAB analyzing program. We find a minimum in the excitation function of the pion freeze-out volume at low transverse momenta and around $E_{lab}\sim 20-30A$GeV which can be related to the transition from hadronic to string matter (which might be interpreted as a pre-cursor of the QGP). The existence of the minimum is explained by the competition of two mechanisms of the particle production, resonance decays and string formation/fragmentation.

preprint2010arXiv

Transport model study of nuclear stopping in heavy ion collisions over an energy range from 0.09A GeV to 160A GeV

Nuclear stopping in the heavy ion collisions over a beam energy range from SIS, AGS up to SPS is studied in the framework of the modified UrQMD transport model, in which mean field potentials of both formed and "pre-formed" hadrons (from string fragmentation) and medium modified nucleon-nucleon elastic cross sections are considered. It is found that the nuclear stopping is influenced by both the stiffness of the equation of state and the medium modifications of nucleon-nucleon cross sections at SIS energies. At the high SPS energies, the two-bump structure is shown in the experimental rapidity distribution of free protons, which can be understood with the consideration of the "pre-formed" hadron potentials.

preprint2009arXiv

Probing the momentum dependence of medium modifications of the nucleon-nucleon elastic cross sections

The momentum dependence of the medium modifications on nucleon-nucleon elastic cross sections is discussed with microscopic transport theories and numerically investigated with an updated UrQMD microscopic transport model. The semi-peripheral Au+Au reaction at beam energy $E_b=400A$ MeV is adopted as an example. It is found that the uncertainties of the momentum dependence on medium modifications of cross sections influence the yields of free nucleons and their collective flows as functions of their transverse momentum and rapidity. Among these observables, the elliptic flow is sensitively dependent on detailed forms of the momentum dependence and more attention should be paid. The elliptic flow is hardly influenced by the probable splitting effect of the neutron-neutron and proton-proton cross sections so that one might pin down the mass splitting effect of the mean-field level at high beam energies and high nuclear densities by exploring the elliptic flow of nucleons or light clusters.

preprint2005arXiv

Particle transfer and fusion cross-section for Super-heavy nuclei in dinuclear system

Within the dinuclear system (DNS) conception, instead of solving Fokker-Planck Equation (FPE) analytically, the Master equation is solved numerically to calculate the fusion probability of super-heavy nuclei, so that the harmonic oscillator approximation to the potential energy of the DNS is avoided. The relative motion concerning the energy, the angular momentum, and the fragment deformation relaxations is explicitly treated to couple with the diffusion process, so that the nucleon transition probabilities, which are derived microscopically, are time-dependent. Comparing with the analytical solution of FPE, our results preserve more dynamical effects. The calculated evaporation residue cross sections for one-neutron emission channel of Pb-based reactions are basically in agreement with the known experimental data within one order of magnitude.

preprint2005arXiv

Probing the density dependence of the symmetry potential at low and high densities

We investigate the sensitivity of several observables to the density dependence of the symmetry potential within the microscopic transport model UrQMD (ultrarelativistic quantum molecular dynamics model). The same systems are used to probe the symmetry potential at both low and high densities. The influence of the symmetry potentials on the yields of $π^{-}$, $π^{+}$, the $π^{-}/π^{+}$ ratio, the $n/p$ ratio of free nucleons and the $t/^3$He ratio are studied for neutron-rich heavy ion collisions ($^{208}Pb+^{208}Pb$, $^{132}Sn+^{124}Sn$, $^{96}Zr+^{96}Zr$) at $E_b=0.4A {\rm GeV}$. We find that these multiple probes provides comprehensive information on the density dependence of the symmetry potential.

preprint2001arXiv

Probing equilibration with respect to isospin degree of freedom in intermediate energy heavy ion collisions

We have studied equilibration with respect to isospin degree of freedom in four 96 mass systems $^{96}Ru+^{96}Ru$, $^{96}Ru+^{96}Zr$, $^{96}Zr+^{96}Ru$, $^{96}Zr+^{96}Zr$ at 100 AMeV and 400 AMeV with isospin dependent QMD. We propose that the neutron-proton differential rapidity distribution is a sensitive probe to the degree of equilibration with respect to the isospin degree of freedom. By analyzing the average $N/Z$ ratio of emitted nucleons, light charged particles (LCP) and intermediate mass fragments (IMF), it is found that there exist memory effect in multifragmentation process. The average $N/Z$ ratio of IMF reduces largely as beam energy increases from 100 AMeV to 400 AMeV which may result from the change of the behavior of the isotope distribution of IMF charges. The isotope distribution of IMF charges does also show certain memory effect at 100 AMeV case but not at 400 AMeV case.