Lie-Wen Chen

Constraints on the Symmetry Energy from PREX-II in the Multimessenger Era

The neutron skin thickness $Δr_{\rm{np}}$ of heavy nuclei is essentially determined by the symmetry energy density slope $L({ρ})$ at $ρ_c = 0.11~{\rm {fm}^{-3}}\approx 2/3ρ_0$ ( $ρ_0$ is nuclear saturation density), roughly corresponding to the average density of finite nuclei. The PREX collaboration recently reported a model-independent extraction of $Δr^{208}_{\rm{np}} = 0.283 \pm 0.071$ fm for the $Δr_{\rm{np}}$ of $^{208}$ Pb, suggesting a rather stiff symmetry energy $E_{\rm{sym}}({ρ})$ with $L({ρ_c }) \ge 52$ MeV. We show that the $E_{\rm{sym}}({ρ})$ cannot be too stiff and $L({ρ_c }) \le 73$ MeV is necessary to be compatible with (1) the ground-state properties and giant monopole resonances of finite nuclei, (2) the constraints on the equation of state of symmetric nuclear matter at suprasaturation densities from flow data in heavy-ion collisions, (3) the largest neutron star (NS) mass reported so far for PSR J0740+6620, (4) the NS tidal deformability extracted from gravitational wave signal GW170817 and (5) the mass-radius of PSR J0030+045 measured simultaneously by NICER. This allows us to obtain $52 \le L({ρ_c }) \le 73$ MeV and $0.212 \le Δr^{208}_{\rm{np}} \le 0.271$ fm, and further $E_{\rm{sym}}({ρ_0 }) = 34.3 \pm 1.7$ MeV, $L({ρ_0 }) = 83.1 \pm 24.7$ MeV, and $E_{\rm{sym}}({2ρ_0 }) = 62.8 \pm 15.9$ MeV. A number of critical implications on nuclear physics and astrophysics are discussed.

Effects of isoscalar- and isovector-scalar meson mixing on neutron star structure

Based on the accurately calibrated interaction FSUGold, we show that including isovector scalar $δ$ meson and its coupling to isoscalar scalar $σ$ meson in the relativistic mean field (RMF) model can soften the symmetry energy $E_{\rm{sym}}(n)$ at intermediate densities while stiffen the $E_{\rm{sym}}(n)$ at high densities. We find this new RMF model can be simultaneously compatible with (1) the constraints on the equation of state of symmetric nuclear matter at suprasaturation densities from flow data in heavy-ion collisions, (2) the neutron skin thickness of $^{208}$ Pb from the PREX-II experiment, (3) the largest mass of neutron star (NS) reported so far from PSR J0740+6620, (4) the limit of $Λ_{1.4}\leq580$ for the dimensionless tidal deformability of the canonical 1.4 $M_{\odot}$ NS from the gravitational wave signal GW170817, (5) the mass-radius relation of PSR J0030+0451 and PSR J0740+6620 measured by NICER, and thus remove the tension between PREX-II and GW170817 observed in the conventional RMF model.

Measuring neutron skin by grazing isobaric collisions

Neutron skin thickness ( $Δr_{\rm np}$ ) of nuclei and the inferred nuclear symmetry energy are of critical importance to nuclear physics and astrophysics. It is traditionally measured by nuclear processes with significant theoretical uncertainties. We recently proposed an indirect measurement of the $Δr_{\rm np}$ by charged hadron multiplicities in central isobaric collisions at relativistic energies, which are sensitive to nuclear densities. In this Letter, we propose a direct measurement of the $Δr_{\rm np}$ by using net-charge multiplicities in ultra-peripheral (grazing) collisions of those isobars, under the assumption that they are simple superimposition of nucleon-nucleon interactions. We illustrate this novel approach by the TRENTO and URQMD models.

Spinodal Enhancement of Light Nuclei Yield Ratio in Relativistic Heavy Ion Collisions

Using a relativistic transport model to describe the evolution of the quantum chromodynamic matter produced in Au+Au collisions at $\sqrt{s_{NN}}=3-200$ GeV, we study the effect of a first-order phase transition in the equation of state of this matter on the yield ratio $N_tN_p/ N_d^2$ ( $tp/d^2$ ) of produced proton ( $p$ ), deuteron ( $d$ ), and triton ( $t$ ). We find that the large density inhomogeneities generated by the spinodal instability during the first-order phase transition can survive the fast expansion of the subsequent hadronic matter and lead to an enhanced $tp/d^2$ in central collisions at $\sqrt{s_{NN}}=3-5$ GeV as seen in the experiments by the STAR Collaboration and the E864 Collaboration. However, this enhancement subsides with increasing collision centrality, and the resulting almost flat centrality dependence of $tp/d^2$ at $\sqrt{s_{NN}}=3$ GeV can also be used as a signal for the first-order phase transition.

Transport Model Comparison Studies of Intermediate-Energy Heavy-Ion Collisions

Transport models are the main method to obtain physics information from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. Calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, and more realistic calculations of heavy-ion collisions. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. We review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known as BUU and QMD type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. We often indicate improved strategies in performing transport simulations and thus provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review.

preprint2021arXiv

Supernova Preshock Neutronization Burst as a Probe of Non-Standard Neutrino Interactions

Core-collapse supernova (CCSN) provides a unique astrophysical site for studying neutrino-matter interactions. Prior to the shock-breakout neutrino burst during the collapse of the iron core, a preshock $ν_e$ burst arises from the electron capture of nuclei and it is sensitive to the low-energy coherent elastic neutrino-nucleus scattering (CE $ν$ NS) which dominates the neutrino opacity. Since the CE $ν$ NS depends strongly on the non-standard neutrino interactions (NSI) which are completely beyond the standard model and yet to be determined, the detection of the preshock burst thus provides a clean way to extract the NSI information. Within the spherically symmetric general-relativistic hydrodynamic simulation for the CCSN, we investigate the NSI effects on the preshock burst. We find that the NSI can maximally enhance the peak luminosity of the preshock burst almost by a factor of three, reaching a value comparable with that of the shock-breakout burst. Future detection of the preshock burst will have critical implications on astrophysics, neutrino physics and physics beyond the standard model.

preprint2020arXiv

Constraining the in-medium nucleon-nucleon cross section from the width of nuclear giant dipole resonance

We develop a new lattice Hamiltonian method for solving the Boltzmann-Uehling-Uhlenbeck (BUU) equation. Adopting the stochastic approach to treat the collision term and using the GPU parallel computing to carry out the calculations allows for a rather high accuracy in evaluating the collision term, especially its Pauli blocking, leading thus to a new level of precision in solving the BUU equation. Applying this lattice BUU method to study the width of giant dipole resonance (GDR) in nuclei, where the accurate treatment of the collision term is crucial, we find that the obtained GDR width of $^{208}{\rm Pb}$ shows a strong dependence on the in-medium nucleon-nucleon cross section $σ_{\rm NN}^*$ . A very large medium reduction of $σ_{\rm NN}^*$ is needed to reproduce the measured value of the GDR width of $^{208}{\rm Pb}$ at the Research Center for Nuclear Physics in Osaka, Japan.

preprint2019arXiv

Ruling out the supersoft high-density symmetry energy from the discovery of PSR J0740+6620 with mass $2.14^{+0.10}_{-0.09}M_\odot$

Using the very recently reported mass $2.14^{+0.10}_{-0.09}M_\odot$ of PSR J0740+6620 together with the data of finite nuclei and the constraints on the equation of state of symmetric nuclear matter at suprasaturation densities from flow data in heavy-ion collisions, we show that the symmetry energy $E_{\rm sym}(n)$ cannot be supersoft so that it becomes negative at suprasaturation densities in neutron stars (NSs) and thus may make the NS have a pure neutron matter core. This is in contrast to the fact that using the mass $2.01 \pm 0.04 M_\odot$ of PSR J0348+0432 as the NS maximum mass cannot rule out the supersoft high-density $E_{\rm sym}(n)$ . Furthermore, we find the stiffer high-density $E_{\rm sym}(n)$ based on the existence of $2.14M_\odot$ NSs leads to a strong constraint of $Λ_{1.4} \ge 348^{+88}_{-51}$ for the dimensionless tidal deformability of the canonical $1.4M_\odot$ NS.

Antimatter $^4_Λ$ H Hypernucleus Production and the $^3_Λ$ H/ $^3$ He Puzzle in Relativistic Heavy-Ion Collisions

We show that the measured yield ratio $^3_Λ$ H/ $^3$ He ( $^3_{\overlineΛ}\overline{\text{H}}$ / $^3\overline{\text{He}}$ ) in Au+Au collisions at $\sqrt{s_{NN}}=200$ GeV and in Pb+Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV can be understood within a covariant coalescence model if (anti-) $Λ$ particles freeze out earlier than (anti-)nucleons but their relative freezeout time is closer at $\sqrt{s_{NN}}=2.76$ TeV than at $\sqrt{s_{NN}}=200$ GeV. The earlier (anti-) $Λ$ freezeout can significantly enhance the yield of (anti)hypernucleus $^4_Λ$ H ( $^4_{\overlineΛ}\overline{\text{H}}$ ), leading to that $^4_{\overlineΛ}\overline{\text{H}}$ has a comparable abundance with $^4\overline{\text{He}}$ and thus provides an easily measured antimatter candidate heavier than $^4\overline{\text{He}}$ . The future measurement on $^4_Λ$ H ( $^4_{\overlineΛ}\overline{\text{H}}$ ) would be very useful to understand the (anti-) $Λ$ freezeout dynamics and the production mechanism of (anti)hypernuclei in relativistic heavy-ion collisions.

Extended Skyrme interactions for nuclear matter, finite nuclei and neutron stars

Recent progress in theory, experiment and observation challenges the mean field models using the conventional Skyrme interaction, suggesting that the extension of the conventional Skyrme interaction is necessary. In this work, by fitting the experimental data of a number of finite nuclei together with a few additional constraints on nuclear matter using the simulated annealing method, we construct three Skyrme interaction parameter sets, namely, eMSL07, eMSL08 and eMSL09, based on an extended Skyrme interaction which includes additional momentum and density dependent two-body forces to effectively simulate the momentum dependence of the three-body force. The three new interactions can reasonably describe the ground-state properties and the isoscalar giant monopole resonance energies of various spherical nuclei used in the fit as well as the ground-state properties of many other spherical nuclei, nicely conform to the current knowledge on the equation of state of asymmetric nuclear matter, eliminate the notorious unphysical instabilities of symmetric nuclear matter and pure neutron matter up to a very high density of $1.2$ fm $^{-3}$ , and simultaneously support heavier neutron stars with mass larger than two times solar mass. One important difference of the three new interactions is about the prediction of the symmetry energy at supra-saturation densities, and these new interactions are thus potentially useful for the determination of the largely uncertain high-density symmetry energy in future. In addition, a comparison is made for the predictions of nuclear matter, finite nuclei and neutron stars with the three new interactions versus those with three typical interactions BSk22, BSk24 and BSk26 from Brussels group.

Hidden-charm Pentaquark States in Heavy Ion Collisions at the Large Hadron Collider

In the framework of the quark combination, we derive the yield formulas and study the yield ratios of the hidden-charm pentaquark states in ultra-relativistic heavy ion collisions. We propose some interesting yield ratios which clearly exhibit the production relationships between different hidden-charm pentaquark states. We show how to employ a specific quark combination model to evaluate the yields of exotic $P_c^+(4380)$ , $P_c^+(4450)$ and their partners on the basis of reproducing the yields of normal identified hadrons, and execute the calculations in central Pb+Pb collisions at $\sqrt{s_{NN}}= 2.76$ TeV as an example.

Isospin dependence of nucleon effective masses in neutron-rich matter

In this talk, we first briefly review the isospin dependence of the total nucleon effective mass $M^{\ast}_{J}$ inferred from analyzing nucleon-nucleus scattering data within an isospin dependent non-relativistic optical potential model, and the isospin dependence of the nucleon E-mass $M^{\ast,\rm{E}}_{J}$ obtained from applying the Migdal-Luttinger theorem to a phenomenological single-nucleon momentum distribution in nuclei constrained by recent electron-nucleus scattering experiments. Combining information about the isospin dependence of both the nucleon total effective mass and E-mass, we then infer the isospin dependence of nucleon k-mass using the well-known relation $M^{\ast}_{J}=M^{\ast,\rm{E}}_{J}\cdot M^{\ast,\rm{k}}_{J}$ . Implications of the results on the nucleon mean free path (MFP) in neutron-rich matter are discussed.

Isospin properties of quark matter from a 3-flavor NJL model

We have studied the properties of hot and dense quark matter based on the 3-flavor Nambu-Jona-Lasinio (NJL) model as well as its Polyakov-loop extension (pNJL) with scalar-isovector and vector-isovector couplings. Provided a considerable large isospin asymmetry or isospin chemical potential, isospin splittings of constituent mass, chiral phase transition boundary, and critical point for $u$ and $d$ quarks can be observed for positive isovector coupling constants but are suppressed for negative ones. The quark matter symmetry energy decreases with the increasing isovector coupling constant, and is mostly enhanced in the pNJL model than in the NJL model. A positive scalar-isovector coupling constant is more likely to lead to an unstable isospin asymmetric quark matter. The isovector coupling has been further found to affect particle fractions as well as the equation of state in hybrid stars. Possible effects on the isospin properties of quark matter have also been discussed if the strangeness sector is further broken among the flavor symmetry.

Isospin splitting of nucleon effective mass from giant resonances in $^{208}$ Pb

Based on mean field calculations with Skyrme interactions, we extract a constraint on the isovector effective mass in nuclear matter at saturation density $ρ_0$ , i.e., $m_{v}^{\ast}(ρ_0)=(0.77\pm0.03) m$ by combining the experimental data of the centroid energy of the isovector giant dipole resonance (IVGDR) and the electric dipole polarizability $α_{\mathrm{D}}$ in $^{208}$ Pb. Meanwhile, the isoscalar effective mass at $ρ_0$ is determined to be $m_{s}^{\ast}(ρ_0)=(0.91\pm0.05) m$ by analyzing the measured excitation energy of the isoscalar giant quadrupole resonance (ISGQR) in $^{208}$ Pb. From the constrained $m_{s}^{\ast}(ρ_0)$ and $m_{v}^{\ast}(ρ_0)$ , we obtain the isospin splitting of nucleon effective mass in asymmetric nuclear matter of isospin asymmetry $δ$ at $ρ_0$ as $[m_n^{\ast}(ρ_0,δ)-m_p^{\ast}(ρ_0,δ)]/m = Δm^*_1(ρ_0) δ+ O(δ^3)$ with the linear isospin splitting coefficient $Δm^*_1(ρ_0) = 0.33\pm0.16$ . We notice that using the recently corrected data on the $α_{\mathrm{D}}$ in $^{208}$ Pb with the contribution of the quasideuteron effect subtracted slightly enhances the isovector effective mass to $m_{v}^{\ast}(ρ_0)=(0.80\pm0.03) m$ and reduces the linear isospin splitting coefficient to $Δm^*_1(ρ_0) = 0.27\pm0.15$ . Furthermore, the constraints on $m_{v}^{\ast}(ρ)$ , $m_{s}^{\ast}(ρ)$ and $Δm^*_1(ρ)$ at other densities are obtained from the similar analyses and we find that the $Δm^*_1(ρ)$ increases with the density.

Production of $ΛΛ$ and $\bar{Λ\text{n}}$ in central Pb+Pb collisions at $\sqrt{s_{NN}}$ =2.76 TeV within a covariant coalescence model

We study the production of $ΛΛ$ and $\overline{Λ\text{n}}$ exotic states in central Pb+Pb collisions at $\sqrt{s_{NN}}=2.76$ TeV at LHC via both hadron and quark coalescence within a covariant coalescence model with a blast-wave-like parametrization for the phase-space configurations of constituent particles at freezeout. In the hadron coalescence, the two states are considered as molecular states while they are considered as six-quark states in the quark coalescence.For $\overline{Λ\text{n}}$ , we find that the yields of both molecular and six-quark states are much larger than the experimental upper-limits. For $ΛΛ$ , while the molecule-state yield is much larger than the experimental upper-limits, the six-quark-state yield could be lower than the upper-limits. The higher molecule-state yields are mainly due to the large contribution of short-lived strong resonance decays into (anti-)nucleons and (anti-) $Λ$ which can significantly enhance the molecule-state yields of $ΛΛ$ and $\overline{Λ\text{n}}$ via hadron coalescence. Our results suggest that the current experimental measurement at LHC cannot exclude the existence of the $ΛΛ$ as an exotic six-quark state, and if $ΛΛ$ is a six-quark state, it is then on the brink of being discovered.

Proton-skins in momentum space and neutron-skins in coordinate space in heavy nuclei

Neutron-skins in coordinate space and proton-skins in momentum space are predicted to coexist in heavy nuclei and their correlation is governed by Liouville's theorem and Heisenberg's uncertainty principle. An analysis of their correlation within a further extended Thomas-Fermi (ETF $^+$ ) approximation incorporating effects of nucleon short-range correlations reveals generally that protons move faster than neutrons in neutron-skins of heavy nuclei.

Strange Quark Stars as Probe of Dark Matter

We demonstrate that the observation of old strange quark stars (SQSs) can set important limits on the scattering cross sections $σ_q$ between the light quarks and the non-interacting scalar dark matter (DM). By analyzing a set of 1403 of solitary pulsarlike compact stars in the Milky Way, we find the old solitary pulsar PSR J1801-0857D can set the most stringent upper limits on $σ_q$ or the DM-proton scattering cross sections $σ_p$ . By converting $σ_q$ into $σ_p$ based on effective operator analyses, we show the resulting $σ_p$ limit by assuming PSR J1801-0857D to be a SQS could be comparable with that of the current direct detection experiments but much weaker (by several orders of magnitude) than that obtained by assuming PSR J1801-0857D to be a neutron star (NS), which requires an extremely small $σ_p$ far beyond the limits of direct detection experiments. Our findings imply that the old pulsars are favored to be SQSs rather than NSs if the scalar DM were observed by future terrestrial experiments.

Understanding transport simulations of heavy-ion collisions at 100 and 400 AMeV: Comparison of heavy ion transport codes under controlled conditions

Transport simulations are very valuable for extracting physics information from heavy-ion collision experiments. With the emergence of many different transport codes in recent years, it becomes important to estimate their robustness in extracting physics information from experiments. We report on the results of a transport code comparison project. 18 commonly used transport codes were included in this comparison: 9 Boltzmann-Uehling-Uhlenbeck-type codes and 9 Quantum-Molecular-Dynamics-type codes. These codes have been required to simulate Au+Au collisions using the same physics input for mean fields and for in-medium nucleon-nucleon cross sections, as well as the same initialization set-up, the impact parameter, and other calculational parameters at 100 and 400 AMeV incident energy. Among the codes we compare one-body observables such as rapidity and transverse flow distributions. We also monitor non-observables such as the initialization of the internal states of colliding nuclei and their stability, the collision rates and the Pauli blocking. We find that not completely identical initializations constitute partly for different evolutions. Different strategies to determine the collision probabilities, and to enforce the Pauli blocking, also produce considerably different results. There is a substantial spread in the predictions for the observables, which is much smaller at the higher incident energy. We quantify the uncertainties in the collective flow resulting from the simulation alone as about $30\%$ at 100 AMeV and $13\%$ at 400 AMeV, respectively. We propose further steps within the code comparison project to test the different aspects of transport simulations in a box calculation of infinite nuclear matter. This should, in particular, improve the robustness of transport model predictions at lower incident energies where abundant amounts of data are available.

Electric Dipole Polarizability in $^{208}$ Pb as a Probe of the Symmetry Energy and Neutron Matter around $ρ_0/3$

It is currently a big challenge to accurately determine the symmetry energy $E_{\text{sym}}(ρ)$ and the pure neutron matter equation of state $E_{\text{PNM}}(ρ)$ , even their values around saturation density $ρ_0$ . We find that the electric dipole polarizability $α_ {\text{D}}$ in $^{208}$ Pb can be determined uniquely by the magnitude of the $E_{\text{sym}}(ρ)$ or almost equivalently the $E_{\text{PNM}}(ρ)$ at subsaturation densities around $ρ_0/3$ , shedding a light upon the genuine correlation between the $α_ {\text{D}}$ and the $E_{\text{sym}}(ρ)$ . By analyzing the experimental data of the $α_ {\text{D}}$ in $^{208}$ Pb from RCNP using a number of non-relativistic and relativistic mean-field models, we obtain very stringent constraints on $E_{\text{sym}}(ρ)$ and $E_{\text{PNM}}(ρ)$ around $ρ_0/3$ . The obtained constraints are found to be in good agreement with the results extracted in other analyses. In particular, our results provide for the first time the experimental constraints on $E_{\text{PNM}}(ρ)$ around $ρ_0/3$ , which are in harmony with the recent determination of $E_{\text{PNM}}(ρ)$ from microscopic theoretical studies and potentially useful in constraining the largely uncertain many-nucleon interactions in microscopic calculations of neutron matter.

Influence of neutron-skin thickness on $π^{-}/π^{+}$ ratio in Pb+Pb collisions

Within an isospin- and momentum-dependent transport model IBUU11 using as an input nucleon density profiles from Hartree-Fock calculations based on a modified Skyrme-like (MSL) model, we study the influence of the uncertainty of the neutron skin thickness on the $π^{-}/π^{+}$ ratio in both central and peripheral Pb+Pb collisions at beam energies of 400 MeV/nucleon and 1000 MeV/nucleon. Within the current experimental uncertainty range of neutron skin in $^{208}$ Pb, while the neutron skin effect on the \rpi ratio is negligible in central reactions at both energies, it increases gradually with increasing impact parameter and becomes comparable with or even larger than the symmetry energy effect in peripheral collisions especially at 400 MeV/nucleon. Moreover, we found that while the \rpi ratio is larger with a softer \esym in central collisions, above certain impact parameters depending on the size of the neutron skin, a stiffer \esym can lead to a larger \rpi ratio as most of the pions are produced at densities below the saturation density in these peripheral reactions. Thus, a clear impact parameter selection is important to extract reliable information about the \esym at suprasaturation densities (size of neutron skin) from the $π^-/π^+$ ratio in central (peripheral) heavy-ion collisions.

Neutron-proton effective mass splitting in neutron-rich matter and its impacts on nuclear reactions

The neutron-proton effective mass splitting in neutron-rich nucleonic matter reflects the space-time nonlocality of the isovector nuclear interaction. It affects the neutron/proton ratio during the earlier evolution of the Universe, cooling of protoneutron stars, structure of rare isotopes and dynamics of heavy-ion collisions. While there is still no consensus on whether the neutron-proton effective mass splitting is negative, zero or positive and how it depends on the density as well as the isospin-asymmetry of the medium, significant progress has been made in recent yeas in addressing these issues. We first recall the connections among the neutron-proton effective mass splitting, the momentum dependence of the isovector potential and the density dependence of the symmetry energy. We then make a few observations about the progress in calculating the neutron-proton effective mass splitting using various nuclear many-body theories and its effects on the isospin-dependence of in-medium nucleon-nucleon cross sections. Perhaps, our most reliable knowledge so far about the neutron-proton effective mass splitting at saturation density of nuclear matter comes from optical model analyses of huge sets of nucleon-nucleus scattering data accumulated over the last five decades. The momentum dependence of the symmetry potential from these analyses provide a useful boundary condition at saturation density for calibrating nuclear many-body calculations. Several observables in heavy-ion collisions have been identified as sensitive probes of the neutron-proton effective mass splitting in dense neutron-rich matter based on transport model simulations. We review these observables and comment on the latest experimental findings.

Neutron-proton effective mass splitting in neutron-rich matter at normal density from analyzing nucleon-nucleus scattering data within an isospin dependent optical model

The neutron-proton effective $k$ -mass splitting in asymmetric nucleonic matter of isospin asymmetry $δ$ and normal density is found to be $m^{*}_{n-p}\equiv(m^{*}_{n}-m^{*}_{p})/m=(0.41 \pm0.15)δ$ from analyzing globally 1088 sets of reaction and angular differential cross sections of proton elastic scattering on 130 targets with beam energies from 0.783 MeV to 200 MeV, and 1161 sets of data of neutron elastic scattering on 104 targets with beam energies from 0.05 MeV to 200 MeV within an isospin dependent non-relativistic optical potential model. It sets a useful reference for testing model predictions on the momentum dependence of the nucleon isovector potential necessary for understanding novel structures and reactions of rare isotopes.

Positioning the neutron drip line and the r-process paths in the nuclear landscape

Exploring nucleon drip lines and astrophysical rapid neutron capture process (r-process) paths in the nuclear landscape is extremely challenging in nuclear physics and astrophysics. While various models predict similar proton drip line, their predictions for neutron drip line and the r-process paths involving heavy neutron-rich nuclei exhibit a significant variation which hampers our accurate understanding of the r-process nucleosynthesis mechanism. Using microscopic density functional theory with a representative set of non-relativistic and relativistic interactions, we demonstrate for the first time that this variation is mainly due to the uncertainty of nuclear matter symmetry energy $E_{\rm{sym}}(ρ_{\rm{sc}})$ at the subsaturation cross density $ρ_{\rm{sc}}=0.11/0.16\timesρ_0$ ( $ρ_0$ is saturation density), which reflects the symmetry energy of heavy nuclei. Using the recent accurate constraint on $E_{\rm{sym}}(ρ_{\rm{sc}})$ from the binding energy difference of heavy isotope pairs, we obtain quite precise predictions for the location of the neutron drip line, the r-process paths and the number of bound nuclei in the nuclear landscape. Our results have important implications on extrapolating the properties of unknown neutron-rich rare isotopes from the data on known nuclei.

Production of antimatter $^{5,6}$ Li nuclei in central Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV

Combining the covariant coalescence model and a blast-wave-like analytical parametrization for (anti-)nucleon phase-space freezeout configuration, we explore light (anti-)nucleus production in central Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV. Using the nucleon freezeout configuration (denoted by FO1) determined from the measured spectra of protons (p), deutrons (d) and $^{3}$ He, we find the predicted yield of $^{4}$ He is significantly smaller than the experimental data. We show this disagreement can be removed by using a nucleon freezeout configuration (denoted by FO2) in which the nucleons are assumed to freeze out earlier than those in FO1 to effectively consider the effect of large binding energy value of $^{4}$ He. Assuming the binding energy effect also exists for the production of $^5\text{Li}$ , $^5\overline{\text{Li}}$ , $^6\text{Li}$ and $^6\overline{\text{Li}}$ due to their similar binding energy values as $^{4}$ He, we find the yields of these heavier (anti-)nuclei can be enhanced by a factor of about one order, implying that although the stable (anti-) $^6$ Li nucleus is unlikely to be observed, the unstable (anti-) $^5$ Li nucleus could be produced in observable abundance in Au+Au collisions at $\sqrt{s_{NN}} = 200$ GeV where it may be identified through the p- $^4\text{He}$ ( $\overline{\text{p}}$ - $^4\overline{\text{He}}$ ) invariant mass spectrum. The future experimental measurement on (anti-) $^5\text{Li}$ would be very useful to understand the production mechanism of heavier antimatter.

Revisit of the neutron/proton ratio puzzle in intermediate-energy heavy-ion collisions

Incorporating a newly improved isospin- and momentum-dependent interaction in the isospin-dependent Boltzmann-Uehling-Uhlenbeck transport model IBUU11, we have investigated relative effects of the density dependence of nuclear symmetry energy $E_{sym}(ρ)$ and the neutron-proton effective mass splitting $m^*_n-m^*_p$ on the neutron/proton ratio of free nucleons and those in light clusters. It is found that the $m^*_n-m^*_p$ has a relatively stronger effect than the $E_{sym}(ρ)$ and the assumption of $m^*_n\leq m^*_p$ leads to a higher neutron/proton ratio. Moreover, this finding is independent of the in-medium nucleon-nucleon cross sections used. However, results of our calculations using the $E_{sym}(ρ)$ and $m^*_n-m^*_p$ both within their current uncertainty ranges are all too low compared to the recent NSCL/MSU double neutron/proton ratio data from central $^{124}$ Sn+ $^{124}$ Sn and $^{112}$ Sn+ $^{112}$ Sn collisions at 50 and 120 MeV/u, thus calling for new mechanisms to explain the puzzlingly high neutron/proton ratio observed in the experiments.

Symmetry energy systematics and its high density behavior

We explore the systematics of the density dependence of nuclear matter symmetry energy in the ambit of microscopic calculations with various energy density functionals, and find that the symmetry energy from subsaturation density to supra-saturation density can be well determined by three characteristic parameters of the symmetry energy at saturation density $ρ_0$ , i.e., the magnitude $E_{\text{sym}}({ρ_0 })$ , the density slope $L$ and the density curvature $K_{\text{sym}}$ . This finding opens a new window to constrain the supra-saturation density behavior of the symmetry energy from its (sub-)saturation density behavior. In particular, we obtain $L=46.7 \pm 12.8$ MeV and $K_{\text{sym}}=-166.9 \pm 168.3$ MeV as well as $E_{\text{sym}}({2ρ_{0}}) \approx 40.2 \pm 12.8$ MeV and $L({2ρ_{0}}) \approx 8.9 \pm 108.7$ MeV based on the present knowledge of $E_{\text{sym}}({ρ_{0}}) = 32.5 \pm 0.5$ MeV, $E_{\text{sym}}({ρ_c}) = 26.65 \pm 0.2$ MeV and $L({ρ_c}) = 46.0 \pm 4.5$ MeV at $ρ_{\rm{c}}= 0.11$ fm $^{-3}$ extracted from nuclear mass and the neutron skin thickness of Sn isotopes. Our results indicate that the symmetry energy cannot be stiffer than a linear density dependence.In addition, we also discuss the quark matter symmetry energy since the deconfined quarks could be the right degree of freedom in dense matter at high baryon densities.

Thermal properties of asymmetric nuclear matter with an improved isospin- and momentum-dependent interaction

Thermal properties of asymmetric nuclear matter, including the temperature dependence of the symmetry energy, single-particle properties, and differential isospin fractionation, are investigated with different neutron-proton effective mass splittings using an improved isospin- and momentum-dependent interaction. In this improved interaction, the momentum-dependence of the isoscalar single-particle potential at saturation density is well fitted to that extracted from optical model analyses of proton-nucleus scattering data up to nucleon kinetic energy of 1 GeV, and the isovector properties, i.e., the slope of the nuclear symmetry energy, the momentum-dependence of the symmetry potential, and the symmetry energy at saturation density can be flexibly adjusted via three parameters $x$ , $y$ , and $z$ , respectively. Our results indicate that the nucleon phase-space distribution in equilibrium, the temperature dependence of the symmetry energy, and the differential isospin fractionation can be significantly affected by the isospin splitting of nucleon effective mass.

Breaking the EOS-Gravity Degeneracy with Masses and Pulsating Frequencies of Neutron Stars

A thorough understanding of many astrophysical phenomena associated with compact objects requires reliable knowledge about both the equation of state (EOS)of super-dense nuclear matter and the theory of strong-field gravity simultaneously because of the EOS-gravity degeneracy. Currently, variations of the neutron star (NS) mass-radius correlation from using alternative gravity theories are much larger than those from changing the NS matter EOS within known constraints. At least two independent observables are required to break the EOS-gravity degeneracy. Using model EOSs for hybrid stars and a Yukawa-type non-Newtonian gravity, we investigate both the mass-radius correlation and pulsating frequencies of NSs. While the maximum mass of NSs increases, the frequencies of the $f$ , $p_1$ , $p_2$ , and $w_I$ pulsating modes are found to decrease with the increasing strength of the Yukawa-type non-Newtonian gravity, providing a useful reference for future determination simultaneously of both the strong-field gravity and the supranuclear EOS by combining data of x-ray and gravitational wave emissions of neutron stars.

Constraining the density slope of nuclear symmetry energy at subsaturation densities using electric dipole polarizability in $^{208}$ Pb

Nuclear structure observables usually most effectively probe the properties of nuclear matter at subsaturation densities rather than at saturation density. We demonstrate that the electric dipole polarizibility $α_ {\text{D}}$ in $^{208}$ Pb is sensitive to both the magnitude $E_{\text{sym}}(ρ_{\text{c}})$ and density slope $L(ρ_{\text{c}})$ of the symmetry energy at a subsaturation cross density $ρ_{\text{c}} = 0.11$ fm $^{-3}$ . Using the experimental data of $α_ {\text{D}}$ in $^{208}$ Pb from RCNP and the recent accurate constraint of $E_{\text{sym}}(ρ_{\text{c}})$ from the binding energy difference of heavy isotope pairs, we extract a value of $L(ρ_{\text{c}}) = 47.3 \pm 7.8$ MeV. The implication of the present constraint of $L(ρ_{\text{c}})$ to the symmetry energy at saturation density, the neutron skin thickness of $^{208}$ Pb and the core-crust transition density in neutron stars is discussed.

Form Factor Effects in the Direct Detection of Isospin-Violating Dark Matter

Isospin-violating dark matter (IVDM) provides a possible mechanism to ameliorate the tension among recent direct detection experiments. For IVDM, we demonstrate that the results of direct detection experiments based on neutron-rich target nuclei may depend strongly on the density dependence of the symmetry energy which is presently largely unknown and controls the neutron skin thickness that reflects the relative difference of neutron and proton form factors in the neutron-rich nuclei. In particular, using the neutron and proton form factors obtained from Skyrme-Hartree-Fock calculations by varying the symmetry energy within the uncertainty region set by the latest model-independent measurement of the neutron skin thickness of $^{208}$ Pb from PREX experiment at JLab, we find that, for IVDM with neutron-to-proton coupling ratio fixed to $f_n/f_p=-0.7$ , the form factor effect may enhance the sensitivity of Xe-based detectors (e.g., XENON100 and LUX) to the DM-proton cross section by a factor of $3$ in the DM mass region constrained by CMDS-II(Si) and even by more than an order of magnitude for heavy DM with mass larger than $80$ GeV, compared with the results using the empirical Helm form factor. Our results further indicate that the form factor effect can significantly modify the recoil spectrum of Xe-based detectors for heavy IVDM with $f_n/f_p=-0.7$ .

Old neutron stars as probes of isospin-violating dark matter

Isospin-violating dark matter (IVDM), which couples differently with protons and neutrons, provides a promising mechanism to ameliorate the tension among recent direct detection experiments. Assuming DM is non-interacting bosonic asymmetric IVDM, we investigate how the existence of old neutron stars limits the DM-proton scattering cross-section $σ_{\rm p}$ , especially the effects of the isospin violating DM-nucleon interactions and the symmetry energy in the equation of state of isospin asymmetric nuclear matter. Our calculations are completely based on general relativity and especially the structure of neutron stars is obtained by solving the Tolman-Oppenheimer-Volkoff equations with nuclear matter equation of state constrained by terrestrial experiments. We find that, by considering the more realistic neutron star model rather than a simple uniform neutron sphere as usual, the $σ_{\rm p}$ bounds from old neutron stars can be varied by more than an order of magnitude depending on the specific values of the DM neutron-to-proton coupling ratio $f_{\rm n}/f_{\rm p}$ , and they can be further varied by more than a factor of two depending on the density dependence of the symmetry energy. In particular, we demonstrate that the observed nearby isolated old neutron star PSR B1257+12 can set a very strong limit on $σ_{\rm p}$ for low-mass DM particles ( $\le 20 \, {\rm GeV}$ ) that reaches a sensitivity beyond the current best limits from direct detection experiments and disfavors the DM interpretation of previously-reported positive experimental results, including the IVDM.

Probing isospin- and momentum-dependent nuclear effective interactions in neutron-rich matter

The single-particle potentials for nucleons and hyperons in neutron-rich matter generally depends on the density and isospin asymmetry of the medium as well as the momentum and isospin of the particle. It further depends on the temperature of the matter if the latter is in thermal equilibrium. We review here the extension of a Gogny-type isospin- and momentum-dependent interaction in several aspects made in recent years and their applications in studying intermediate-energy heavy ion collisions, thermal properties of asymmetric nuclear matter and properties of neutron stars. The importance of the isospin- and momentum-dependence of the single-particle potential, especially the momentum dependence of the isovector potential, is clearly revealed throughout these studies.

Quark magnetar in three-flavor Nambu--Jona-Lasinio model with vector interaction and magnetized gluon potential

We investigate properties of strange quark matter in the framework of SU(3) Nambu--Jona-Lasinio(NJL) model with vector interaction under strong magnetic fields. The effects of vector-isoscalar and vector-isovector interaction on the equation of state of strange quark matter are investigated, and it is found that the equation of state is not sensitive to the vector-isovector interaction, however, a repulsive interaction in the vector-isoscalar channel gives a stiffer equation of state for cold dense quark matter. In the presence of magnetic field, gluons will be magnetized via quark loops, and the contribution from magnetized gluons to the equation of state is also estimated. The sound velocity square is a quantity to measure the hardness or softness of dense quark matter, and in the NJL model without vector interaction at zero magnetic field the sound velocity square is always less than 1/3. It is found that a repulsive vector-isoscalar interaction and a positive pressure contribution from magnetized gluons can enhance the sound velocity square, which can even reach 1. To construct quark magnetars under strong magnetic fields, we consider anisotropic pressures and use a density-dependent magnetic field profile to mimic the magnetic field distribution in a quark star. We also analyze the parameter region for the magnitude of vector-isoscalar interaction and the contribution from magnetized gluons in order to produce 2 solar mass quark magnetars.

Quark stars under strong magnetic fields

Within the confined isospin- and density-dependent mass model, we study the properties of strange quark matter (SQM) and quark stars (QSs) under strong magnetic fields. The equation of state of SQM under a constant magnetic field is obtained self-consistently and the pressure perpendicular to the magnetic field is shown to be larger than that parallel to the magnetic field, implying that the properties of magnetized QSs generally depend on both the strength and the orientation of the magnetic fields distributed inside the stars. Using a density-dependent magnetic field profile which is introduced to mimic the magnetic field strength distribution in a star, we study the properties of static spherical QSs by assuming two extreme cases for the magnetic field orientation in the stars, i.e., the radial orientation in which the local magnetic fields are along the radial direction and the transverse orientation in which the local magnetic fields are randomly oriented but perpendicular to the radial direction. Our results indicate that including the magnetic fields with radial (transverse) orientation can significantly decrease (increase) the maximum mass of QSs, demonstrating the importance of the magnetic field orientation inside the magnetized compact stars.

Constraining the symmetry energy at subsaturation densities using isotope binding energy difference and neutron skin thickness

We show that the neutron skin thickness $Δr_{np}$ of heavy nuclei is uniquely fixed by the symmetry energy density slope $L(ρ)$ at a subsaturation cross density $ρ_c \approx 0.11$ fm $^{-3}$ rather than at saturation density $ρ_0$ , while the binding energy difference $ΔE$ between a heavy isotope pair is essentially determined by the magnitude of the symmetry energy $E_{\text{sym}}(ρ)$ at the same $ρ_c$ . Furthermore, we find a value of $L({ρ_c})$ leads to a negative $E_{\text{sym}}({ρ_{0}})$ - $L({ρ_{0}})$ correlation while a value of $E_{\text{sym}}({ρ_{c})}$ leads to a positive one. Using data on $Δr_{np}$ of Sn isotopes and $ΔE$ of a number of heavy isotope pairs, we obtain simultaneously $E_{\text{sym}}({ρ_{c})}=26.65\pm0.20$ MeV and $L({ρ_c})= 46.0\pm4.5$ MeV at 95% confidence level, whose extrapolation gives $E_{\text{sym}}({ρ_{0}})=32.3\pm1.0$ MeV and $L({ρ_{0}})=45.2\pm10.0$ MeV. The implication of these new constraints on the $Δr_{np}$ of $^{208}$ Pb and the core-crust transition density in neutron stars is discussed.

Energy dependence of pion in-medium effects on π^-/π^+ ratio in heavy-ion collisions

Within the framework of a thermal model with its parameters fitted to the results from an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model, we study the pion in-medium effect on the charged-pion ratio in heavy-ion collisions at various energies. We find that due to the cancellation between the effects from pion-nucleon s-wave and p-wave interactions in nuclear medium, the π^-/π^+ ratio generally decreases after including the pion in-medium effect. The effect is larger at lower collision energies as a result of narrower pion spectral functions at lower temperatures.

Extracting the nuclear symmetry potential and energy from neutron-nucleus scattering data

According to the Hugenholtz-Van Hove theorem, the symmetry energy and its density slope can be decomposed uniquely in terms of the single-nucleon potential in asymmetric nuclear matter which, at normal density, can be constrained by the nucleon optical model potential extracted from analyzing the nucleon-nucleus scattering data. To more accurately extract information on the symmetry energy and its density slope at normal density from neutron-nucleus scattering data, going beyond the well-known Lane potential, we include consistently the second order terms in isospin asymmetry in both the optical model potential and the symmetry energy decomposition. We find that the strength of the second-order symmetry potential $U_{sym,2}$ in asymmetric nuclear matter is significant compared to the first-order one U_{sym,1} $, especially at high nucleon momentum. While the$ U_{sym,1} $at normal density decreases with nucleon momentum, the$ U_{sym,2} $is found to have the opposite behavior. Moreover, we discuss effects of the$ U_{sym,1} $and$ U_{sym,2} $on determining the density dependence of the symmetry energy, and we find that the available neutron-nucleus scattering data favor a softer density dependence of the symmetry energy at normal density with a negligible contribution from the$ U_{sym,2}$.

Nuclear constraints on non-Newtonian gravity at femtometer scale

Effects of the non-Newtonian gravity on properties of finite nuclei are studied by consistently incorporating both the direct and exchange contribution of the Yukawa potential in the Hartree-Fock approach using a well-tested Skyrme force for the strong interaction. It is shown for the first time that the strength of the Yukawa term in the non-Newtonian gravity is limited to $\log(|α|)<1.75/[λ(\rm fm)]^{0.54} + 33.6$ within the length scale of $λ=1-10$ fm in order for the calculated properties of finite nuclei not to be in conflict with accurate experimental data available.

Probing nuclear symmetry energy at high densities using pion, kaon, eta and photon productions in heavy-ion collisions

The high-density behavior of nuclear symmetry energy is among the most uncertain properties of dense neutron-rich matter. Its accurate determination has significant ramifications in understanding not only the reaction dynamics of heavy-ion reactions especially those induced by radioactive beams but also many interesting phenomena in astrophysics, such as the explosion mechanism of supernova and the properties of neutron stars. The heavy-ion physics community has devoted much effort during the last few years to constrain the high-density symmetry using various probes. In particular, the pion-/pion+ ratio has been most extensively studied both theoretically and experimentally. All models have consistently predicted qualitatively that the pion-/pion+ ratio is a sensitive probe of the high-density symmetry energy especially with beam energies near the pion production threshold. However, the predicted values of the pion-/pion+ ratio are still quite model dependent mostly because of the complexity of modeling pion production and reabsorption dynamics in heavy-ion collisions, leading to currently still controversial conclusions regarding the high-density behavior of nuclear symmetry energy from comparing various model calculations with available experimental data. As more pion-/pion+ data become available and a deeper understanding about the pion dynamics in heavy-ion reactions is obtained, more penetrating probes, such as the kaon+/kaon0 ratio, eta meson and high energy photons are also being investigated or planned at several facilities. Here, we review some of our recent contributions to the community effort of constraining the high-density behavior of nuclear symmetry energy in heavy-ion collisions. In addition, the status of some worldwide experiments for studying the high-density symmetry energy, including the HIRFL-CSR external target experiment (CEE) are briefly introduced.

Quark matter symmetry energy and quark stars

We extend the confined-density-dependent-mass (CDDM) model to include isospin dependence of the equivalent quark mass. Within the confined-isospin-density-dependent-mass (CIDDM) model, we study the quark matter symmetry energy, the stability of strange quark matter, and the properties of quark stars. We find that including isospin dependence of the equivalent quark mass can significantly influence the quark matter symmetry energy as well as the properties of strange quark matter and quark stars. While the recently discovered large mass pulsars PSR J1614-2230 and PSR J0348+0432 with masses around $2M_{\odot}$ cannot be quark stars within the CDDM model, they can be well described by quark stars in the CIDDM model. In particular, our results indicate that the two-flavor $u$ - $d$ quark matter symmetry energy should be at least about twice that of a free quark gas or normal quark matter within conventional Nambu-Jona-Lasinio model in order to describe the PSR J1614-2230 and PSR J0348+0432 as quark stars.

Relationship between the symmetry energy and the single-nucleon potential in isospin-asymmetric nucleonic matter

In this contribution, we review the most important physics presented originally in our recent publications. Some new analyses, insights and perspectives are also provided. We showed recently that the symmetry energy $E_{sym}(ρ)$ and its density slope $L(ρ)$ at an arbitrary density $ρ$ can be expressed analytically in terms of the magnitude and momentum dependence of the single-nucleon potentials using the Hugenholtz-Van Hove (HVH) theorem. These relationships provide new insights about the fundamental physics governing the density dependence of nuclear symmetry energy. Using the isospin and momentum (k) dependent MDI interaction as an example, the contribution of different terms in the single-nucleon potential to the $E_{sym}(ρ)$ and $L(ρ)$ are analyzed in detail at different densities. It is shown that the behavior of $E_{sym}(ρ)$ is mainly determined by the first-order symmetry potential $U_{sym,1}(ρ,k)$ of the single-nucleon potential. The density slope $L(ρ)$ depends not only on the first-order symmetry potential $U_{sym,1}(ρ,k)$ but also the second-order one $U_{sym,2}(ρ,k)$ . Both the $U_{sym,1}(ρ,k)$ and $U_{sym,2}(ρ,k)$ at normal density $ρ_0$ are constrained by the isospin and momentum dependent nucleon optical potential extracted from the available nucleon-nucleus scattering data. The $U_{sym,2}(ρ,k)$ especially at high density and momentum affects significantly the $L(ρ)$ , but it is theoretically poorly understood and currently there is almost no experimental constraints known.

Shear viscosity of neutron-rich nucleonic matter near its liquid-gas phase transition

Within a relaxation time approach using free nucleon-nucleon cross sections modified by the in-medium nucleon masses that are determined from an isospin- and momentum-dependent effective nucleon-nucleon interaction, we investigate the specific shear viscosity ( $η/s$ ) of neutron-rich nucleonic matter near its liquid-gas phase transition. It is found that as the nucleonic matter is heated at fixed pressure or compressed at fixed temperature, its specific shear viscosity shows a valley shape in the temperature or density dependence, with the minimum located at the boundary of the phase transition. Moreover, the value of $η/s$ drops suddenly at the first-order liquid-gas phase transition temperature, reaching as low as $4\sim5$ times the KSS bound of $\hbar/4π$ . However, it varies smoothly for the second-order liquid-gas phase transition. Effects of the isospin degree of freedom and the nuclear symmetry energy on the value of $η/s$ are also discussed.

Contributions of hyperon-hyperon scattering to subthreshold cascade production in heavy ion collisions

Using a gauged flavor SU(3)-invariant hadronic Lagrangian, we calculate the cross sections for the strangeness-exchange reactions YY to NΞ(Y=Λ, Σ) in the Born approximation. These cross sections are then used in the Relativistic Vlasov-Uehling-Uhlenbeck (RVUU) transport model to study Ξproduction in Ar+KCl collisions at incident energy of 1.76A GeV and impact parameter b=3.5 fm. We find that including the contributions of hyperon-hyperon scattering channels strongly enhances the yield of Ξ, leading to the abundance ratio Ξ^{-}/(Λ+Σ^{0})=3.38E-3, which is essentially consistent with the recently measured value of $(5.6 \pm 1.2_{-1.7}^{+1.8})\times 10^{-3}$ by the HADES collaboration at GSI.

Correlations between the nuclear breathing mode energy and properties of asymmetric nuclear matter

Based on microscopic Hartree-Fock + random phase approximation calculations with Skyrme interactions, we study the correlations between the nuclear breathing mode energy $E_{ISGMR}$ and properties of asymmetric nuclear matter with a recently developed analysis method. Our results indicate that the $E_{ISGMR}$ of $^{208}$ Pb exhibits moderate correlations with the density slope $L$ of the symmetry energy and the isoscalar nucleon effective mass $m_{s,0}^{\ast}$ besides a strong dependence on the incompressibility $K_{0}$ of symmetric nuclear matter. Using the present empirical values of $L=60\pm 30$ MeV and $m_{s,0}^{\ast}=(0.8\pm 0.1)m$ , we obtain a theoretical uncertainty of about $\pm 16$ MeV for the extraction of $K_{0}$ from the $E_{ISGMR}$ of $^{208}$ Pb. Furthermore, we find the $E_{ISGMR}$ difference between $^{100}$ Sn and $^{132}$ Sn strongly correlates with $L$ and thus provides a potentially useful probe of the symmetry energy.

Effects of hadronic potentials on elliptic flows in relativistic heavy ion collisions

Within the framework of a multiphase transport (AMPT) model that includes both initial partonic and final hadronic interactions, we show that including mean-field potentials in the hadronic phase leads to a splitting of the elliptic flows of particles and their antiparticles, providing thus a plausible explanation of the different elliptic flows between $p$ and ${\bar p}$ , $K^+$ and $K^-$ , and $π^+$ and $π^-$ observed in recent Beam Energy Scan (BES) program at the Relativistic Heavy-Ion Collider (RHIC).

Effects of the U-boson on the inner edge of neutron star crusts

We explore effects of the light vector $U$ -boson, which is weakly coupled to nucleons, on the transition density $ρ_{t}$ and pressure $P_{t}$ at the inner edge separating the liquid core from the solid crust of neutron stars. Three methods, i.e., the thermodynamical approach, the curvature matrix approach and the Vlasov equation approach are used to determine the transition density $ρ_{t}$ with the Skyrme effective nucleon-nucleon interactions. We find that the $ρ_{t}$ and $P_{t}$ depend on not only the ratio of coupling strength to mass squared of the $U$ -boson $g^{2}/μ^{2}$ but also its mass $μ$ due to the finite range interaction from the $U$ -boson exchange. In particular, our results indicate that the $ρ_{t}$ and $P_{t}$ are sensitive to both $g^{2}/μ^{2}$ and $μ$ if the $U$ -boson mass $μ$ is larger than about 2 MeV. Furthermore, we show that both $g^{2}/μ^{2}$ and $μ$ can have significant influence on the mass-radius relation and the crustal fraction of total moment of inertia of neutron stars. In addition, we study the exchange term contribution of the $U$ -boson based on the density matrix expansion method, and demonstrate that the exchange term effects on the nuclear matter equation of state as well as the $ρ_{t}$ and $P_{t}$ are generally negligible.

Large-mass neutron stars with hyperonization

Within a density-dependent relativistic mean-field model using in-medium meson-hadron coupling constants and meson masses, we explore effects of in-medium hyperon interactions on properties of neutron stars. It is found that the hyperonic constituents in large-mass neutron stars can not be simply ruled out, while the recently measured mass of the millisecond pulsar J1614-2230 can constrain significantly the in-medium hyperon interactions. Moreover, effects of nuclear symmetry energy on hyperonization in neutron stars are also discussed.

Lorentz covariant nucleon self-energy decomposition of the nuclear symmetry energy

Using the Hugenholtz-Van Hove theorem, we derive analytical expressions for the nuclear symmetry energy $E_{sym}(ρ)$ and its density slope $L(ρ)$ in terms of the Lorentz covariant nucleon self-energies in isospin asymmetric nuclear matter. These general expressions are useful for determining the density dependence of the symmetry energy and understanding the Lorentz structure and the microscopic origin of the nuclear symmetry energy in relativistic covariant formulism. As an example, we analyze the Lorentz covariant nucleon self-energy decomposition of $E_{sym}(ρ)$ and $L(ρ)$ and derive the corresponding analytical expressions within the nonlinear $σ$ - $ω$ - $ρ$ - $δ$ relativistic mean field model.

Non-Newtonian gravity in finite nuclei

In this talk, we report our recent study of constraining the non-Newtonian gravity at femtometer scale. We incorporate the Yukawa-type non-Newtonian gravitational potential consistently to the Skyrme functional form using the exact treatment for the direct contribution and density-matrix expansion method for the exchange contribution. The effects from the non-Newtonian potential on finite nuclei properties are then studied together with a well-tested Skyrme force. Assuming that the framework without non-Newtonian gravity can explain the binding energies and charge radii of medium to heavy nuclei within 2% error, we set an upper limit for the strength of the non-Newtonian gravitational potential at femtometer scale.

Nuclear matter fourth-order symmetry energy in the relativistic mean field models

Within the nonlinear relativistic mean field model, we derive the analytical expression of the nuclear matter fourth-order symmetry energy $E_{sym,4}(ρ)$ . Based on two accurately calibrated interactions FSUGold and IU-FSU, our results show that the value of $E_{sym,4}(ρ)$ at normal nuclear matter density $ρ_{0}$ is generally less than 1 MeV, confirming the empirical parabolic approximation to the equation of state for asymmetric nuclear matter at $ρ_{0}$ . On the other hand, we find that the $E_{sym,4}(ρ)$ may become nonnegligible at high densities. Furthermore, the analytical form of the $E_{sym,4}(ρ)$ provides the possibility to study the higher-order effects on the isobaric incompressibility of asymmetric nuclear matter, i.e., $K_{sat}(δ)=K_{0}+K_{sat,2}δ^{2}+K_{sat,4}δ^{4}+\mathcal{O}(δ^{6})$ where $δ=(ρ_{n}-ρ_{p})/ρ$ is the isospin asymmetry, and we find that the value of $K_{sat,4}$ is generally small compared with that of the $K_{sat,2}$ . In addition, we study the effects of the $E_{sym,4}(ρ)$ on the proton fraction $x_{p}$ and the core-crust transition density $ρ_{t}$ and pressure $P_{t}$ in neutron stars. Interestingly, we find that, compared with the results from the empirical parabolic approximation, including the $E_{sym,4}(ρ)$ contribution can significantly enhance the $x_{p}$ at high densities and strongly reduce the $ρ_{t}$ and $P_{t}$ in neutron stars, demonstrating that the widely used empirical parabolic approximation may cause large errors in determining the $x_{p}$ at high densities as well as the $ρ_{t}$ and $P_{t}$ in neutron stars within the nonlinear relativistic mean field model, consistent with previous nonrelativistic calculations.

Recent progress on the determination of the symmetry Energy

We summarize the current status on constraining the density dependence of the symmetry energy from terrestrial laboratory measurements and astrophysical observations. While the value $E_{sym}({ρ_{0}})$ and density slope $L$ of the symmetry energy at saturation density $ρ_{0}$ can vary largely depending on the data or methods, all the existing constraints are essentially consistent with $E_{sym}({ρ_{0}}) = 31 \pm 2$ MeV and $L = 50 \pm 20$ MeV. The determination of the supra-saturation density behavior of the symmetry energy remains a big challenge.

Single-nucleon potential decomposition of the nuclear symmetry energy

The nuclear symmetry energy $E_{sym}(ρ)$ and its density slope $L(ρ)$ can be decomposed analytically in terms of the single-nucleon potential in isospin asymmetric nuclear matter. Using three popular nuclear effective interaction models which have been extensively used in nuclear structure and reaction studies, namely, the isospin and momentum dependent MDI interaction model, the Skyrme Hartree-Fock approach and the Gogny Hartree-Fock approach, we analyze the contribution of different terms in the single-nucleon potential to the $E_{sym}(ρ)$ and $L(ρ)$ . Our results show that the observed different density behaviors of $E_{sym}(ρ)$ for different interactions are essentially due to the variation of the symmetry potential $U_{sym,1}(ρ,k)$ . Furthermore, we find that the contribution of the second-order symmetry potential $U_{sym,2}(ρ,k)$ to the $L(ρ)$ generally cannot be neglected. Moreover, our results demonstrate that the magnitude of the $U_{sym,2}(ρ,k)$ is generally comparable with that of $U_{sym,1}(ρ,k)$ , indicating the second-order symmetry potential $U_{sym,2}(ρ,k)$ may have significant corrections to the widely used Lane approximation to the single-nucleon potential in extremely neutron(proton)-rich nuclear matter.

Can the maximum mass of neutron stars rule out any equation of state of dense stellar matter before gravity is well understood?

Probably No! As an example, using soft EOSs consistent with existing terrestrial nuclear laboratory experiments for hybrid neutron stars containing a quark core described with MIT bag model using reasonable parameters, we show that the recently discovered new holder of neutron star maximum mass PSR J1614-2230 of $1.97\pm0.04M_{\odot}$ can be well described by incorporating a Yukawa gravitational correction that is consistent with existing constraints from neutron-proton and neutron-lead scatterings as well as the spectroscopy of antiproton atoms.

Higher order bulk characteristic parameters of asymmetric nuclear matter

The bulk parameters characterizing the energy of symmetric nuclear matter and the symmetry energy defined at normal nuclear density $ρ_0$ provide important information on the equation of state (EOS) of isospin asymmetric nuclear matter. While significant progress has been made in determining some lower order bulk characteristic parameters, such as the energy $E_0(ρ_0)$ and incompressibility $K_0$ of symmetric nuclear matter as well as the symmetry energy $E_{sym}(ρ_0)$ and its slope parameter $L$ , yet the higher order bulk characteristic parameters are still poorly known. Here, we analyze the correlations between the lower and higher order bulk characteristic parameters within the framework of Skyrme Hartree-Fock energy density functional and then estimate the values of some higher order bulk characteristic parameters. In particular, we obtain $J_0=-355 \pm 95$ MeV and $I_0=1473 \pm 680$ MeV for the third-order and fourth-order derivative parameters of symmetric nuclear matter at $ρ_0$ and $K_{sym} = -100 \pm 165$ MeV, $J_{sym} = 224 \pm 385$ MeV, $I_{sym} = -1309 \pm 2025$ MeV for the curvature parameter, third-order and fourth-order derivative parameters of the symmetry energy at $ρ_0$ , using the empirical constraints on $E_0(ρ_0)$ , $K_0$ , $E_{sym}(ρ_0)$ , $L$ , and the isoscalar and isovector nucleon effective masses. Furthermore, our results indicate that the three parameters $E_0(ρ_0)$ , $K_0$ , and $J_0$ can reasonably characterize the EOS of symmetric nuclear matter up to $2ρ_0$ while the symmetry energy up to $2ρ_0$ can be well described by $E_{sym}(ρ_0)$ , $L$ , and $K_{sym}$ .

Imprints of Nuclear Symmetry Energy on Properties of Neutron Stars

Significant progress has been made in recent years in constraining the density dependence of nuclear symmetry energy using terrestrial nuclear laboratory data. Around and below the nuclear matter saturation density, the experimental constraints start to merge in a relatively narrow region. At supra-saturation densities, there are, however, still large uncertainties. After summarizing the latest experimental constraints on the density dependence of nuclear symmetry energy, we highlight a few recent studies examining imprints of nuclear symmetry energy on the binding energy, energy release during hadron-quark phase transitions as well as the $w$ -mode frequency and damping time of gravitational wave emission of neutron stars.

Isospin dependent global neutron-nucleus optical model potential

In this paper, we construct a new phenomenological isospin dependent global neutron-nucleus optical model potential. Based on the existing experimental data of elastic scattering angular distributions for neutron as projectile, we obtain a set of the isospin dependent global neutron-nucleus optical model potential parameters, which can basically reproduce the experimental data for target nuclei from $^{24}$ Mg to $^{242}$ Pu with the energy region up to 200 MeV.

Nuclear matter symmetry energy and the symmetry energy coefficient in the mass formula

Within the Skyrme-Hartree-Fock (SHF) approach, we show that for a fixed mass number A, both the symmetry energy coefficient a_{sym}(A) in the semi-empirical mass formula and the nuclear matter symmetry energy E_{sym}(ρ_A) at a subsaturation reference density rho_A can be determined essentially by the symmetry energy E_{sym}(rho_0) and its density slope L at saturation density rho_0. Meanwhile, we find the dependence of a_{sym}(A) on E_{sym}(rho_0) or L is approximately linear and is very similar to the corresponding linear dependence displayed by E_{sym}(ρ_A), providing an explanation for the relation E_{sym}(ρ_A) \approx a_{sym}(A). Our results indicate that a value of E_{sym}(ρ_A) leads to a linear correlation between E_{sym}(rho_0) and L and thus can put important constraints on E_{sym}(rho_0) and L. Particularly, the values of E_{sym}(rho_0)= 30.5 +- 3 MeV and L= 52.5 +- 20 MeV are simultaneously obtained by combining the constraints from recently extracted E_{sym}(ρ_A=0.1 fm^{-3}) with those from recent analyses of neutron skin thickness of Sn isotopes in the same SHF approach.

Transport Model Studies of the Baryon-Rich Quark-Gluon Plasma formed in Heavy Ion Collisions

Heavy ion collisions in the low energy run at Relativistic Heavy Ion Collider (RHIC) and future Facility for Antiproton and Ion Research (FAIR) in Germany are expected to produce a quark-gluon plasma that has a finite baryon chemical potential, allowing thus the possibility to study the location of the critical endpoint in the QCD phase diagram. In this talk, using a multiphase transport model, that includes interactions in both initial partonic and final hadronic matters and the transition between these two phases of matter, we discuss the effects of partonic interactions on observables such as the elliptic flow that have played essential roles in studying the properties of the net baryon free quark-gluon plasma produced in heavy ion collisions at higher energies at RHIC. Also, we study the effect of density fluctuations due to a first-order transition between the quark-gluon plasma and hadronic matter on fluctuations of hadron mean transverse momentum and produced deuteron number as well as on two-pion correlations. We further discuss the possibility of studying the mechanism of charm energy loss in the baryon-rich quark-gluon plasma and the properties of phi mesons in hot-dense matter produced in these collisions.

Analytical relations between nuclear symmetry energy and single-nucleon potentials in isospin asymmetric nuclear matter

Using the Hugenholtz-Van Hove theorem, we derive general expressions for the quadratic and quartic symmetry energies in terms of single-nucleon potentials in isospin asymmetric nuclear matter. These analytical relations are useful for gaining deeper insights into the microscopic origins of the uncertainties in our knowledge on nuclear symmetry energies especially at supra-saturation densities. As examples, the formalism is applied to two model single-nucleon potentials that are widely used in transport model simulations of heavy-ion reactions.

Density slope of the nuclear symmetry energy from the neutron skin thickness of heavy nuclei

Expressing explicitly the parameters of the standard Skyrme interaction in terms of the macroscopic properties of asymmetric nuclear matter, we show in the Skyrme-Hartree-Fock approach that unambiguous correlations exist between observables of finite nuclei and nuclear matter properties. We find that existing data on neutron skin thickness $Δr_{np}$ of Sn isotopes give an important constraint on the symmetry energy $E_{sym}({ρ_{0}})$ and its density slope $L$ at saturation density ${ρ_{0}}$ . Combining these constraints with those from recent analyses of isospin diffusion and double neutron/proton ratio in heavy-ion collisions at intermediate energies leads to a more stringent limit on $L$ approximately independent of $E_{sym}({ρ_{0}})$ . The implication of these new constraints on the $Δr_{np}$ of $^{208}$ Pb as well as the core-crust transition density and pressure in neutron stars is discussed.

Isospin- and momentum-dependent effective interactions for the baryon octet and the properties of hybrid stars

The isospin- and momentum-dependent MDI interaction, which has been extensively used in intermediate-energy heavy-ion reactions to study the properties of asymmetric nuclear matter, is extended to include the nucleon-hyperon and hyperon-hyperon interactions by assuming same density, momentum and isospin dependence as for the nucleon-nucleon interaction. The parameters in these interactions are determined from the empirical hyperon single-particle potentials in symmetric nuclear matter at saturation density. The extended MDI interaction is then used to study in the mean-field approximation the equation of state of hypernuclear matter and also the properties of hybrid stars by including the phase transition from the hypernuclear matter to the quark matter at high densities. In particular, the effects of attractive and repulsive $Σ$ N interactions and different values of symmetry energies on the hybrid star properties are investigated.

Nuclear symmetry energy and its density slope at normal density extracted from global nucleon optical potentials

Based on the Hugenholtz-Van Hove theorem, it is shown that both the symmetry energy E $_{sym}(ρ)$ and its density slope $L(ρ)$ at normal density $ρ_0$ are completely determined by the global nucleon optical potentials that can be extracted directly from nucleon-nucleus scatterings, (p,n) charge exchange reactions and single-particle energy levels of bound states. Adopting a value of $m^*/m=0.7$ for the nucleon effective k-mass in symmetric nuclear matter at $ρ_0$ and averaging all phenomenological isovector nucleon potentials constrained by world data available in the literature since 1969, the best estimates of $E_{sym}(ρ_0)=31.3$ MeV and $L(ρ_0)=52.7$ MeV are simultaneously obtained. Uncertainties involved in the estimates are discussed.

Transition density and pressure in hot neutron stars

Using the momentum-dependent MDI effective interaction for nucleons, we have studied the transition density and pressure at the boundary between the inner crust and liquid core of hot neutron stars. We find that their values are larger in neutrino-trapped neutron stars than in neutrino-free neutron stars. Furthermore, both are found to decrease with increasing temperature of a neutron star as well as increasing slope parameter of the nuclear symmetry energy, except that the transition pressure in neutrino-trapped neutron stars for the case of small symmetry energy slope parameter first increases and then decreases with increasing temperature. We have also studied the effect of the nuclear symmetry energy on the critical temperature above which the inner crust in a hot neutron star disappears and found that with increasing value of the symmetry energy slope parameter, the critical temperature decreases slightly in neutrino-trapped neutron stars but first decreases and then increases in neutrino-free neutron stars.

Circumstantial evidence for a soft nuclear symmetry energy at supra-saturation densities

Within an isospin- and momentum-dependent hadronic transport model it is shown that the recent FOPI data on the $π^-/π^+$ ratio in central heavy-ion collisions at SIS/GSI energies (Willy Reisdorf {\it et al.}, NPA {\bf 781}, 459 (2007)) provide circumstantial evidence suggesting a rather soft nuclear symmetry energy \esym at $ρ\geq 2ρ_0$ compared to the Akmal-Pandharipande-Ravenhall prediction. Some astrophysical implications and the need for further experimental confirmations are discussed.

Incompressibility of asymmetric nuclear matter

The incompressibility $K_sat(δ)$ of isospin asymmetric nuclear matter at its saturation density. Our results show that in the expansion of $K_sat(δ)$ in powers of isospin asymmetry $δ$ , i.e., $K_sat(δ)$ =K_{0}+K_{sat,2}δ^{2}+K_{sat,4}δ^{4}+O(δ^{6})$, the magnitude of the 4th-order K_{sat,4} parameter is generally small. The 2nd-order K_{sat,2} parameter thus essentially characterizes the isospin dependence of the incompressibility of asymmetric nuclear matter at saturation density. Furthermore, the K_{sat,2} can be expressed as K_{sat,2}=K_{sym}-6L-J_{0}/{K_{0}L in terms of the slope parameter $L $and the curvature parameter$ K_{\mathrm{sym}} $of the symmetry energy and the third-order derivative parameter$ J_0 $of the energy of symmetric nuclear matter at saturation density, and we find the higher order$ J_0 $contribution to K_{sat,2} generally cannot be neglected. Also, we have found a linear correlation between K_{sym} and$ L $as well as between$ J_{0}/K_{0} $and$ K_{0} $. Using these correlations together with the empirical constraints on$ K_{0} $and$ L $, the nuclear symmetry energy$ E_sym(ρ_{0})$ at normal nuclear density, and the nucleon effective mass, we have obtained an estimated value of K_{sat,2}=-370 +- 120 MeV for the 2nd-order parameter in the isospin asymmetry expansion of the incompressibility of asymmetric nuclear matter at its saturation density.

Nuclear constraints on properties of neutron star crusts

The transition density $ρ_{t}$ and pressure $P_{t}$ at the inner edge separating the liquid core from the solid crust of neutron stars are systematically studied using a modified Gogny (MDI) and 47 popular Skyrme interactions within well established dynamical and thermodynamical methods. It is shown that the widely used parabolic approximation to the full Equation of State (EOS) of isospin asymmetric nuclear matter may lead to huge errors in estimating the ρ_{t} and P_{t}, especially for stiffer symmetry energy functionals $E_{sym}(ρ)$ . The ρ_{t} and P_{t} decrease roughly linearly with the increasing slope parameter $L$ of the $E_{sym}(ρ)$ using the full EOS within both methods. It is also shown that the thickness, fractional mass and moment of inertia of neutron star crust are all very sensitive to the parameter $L$ through the $ρ_{t}$ . Moreover, it is shown that the $E_{sym}(ρ)$ constrained in the same sub-saturation density range as the neutron star crust by the isospin diffusion data in heavy-ion collisions at intermediate energies limits the transition density and pressure to 0.040$ fm^-3}< ρ_{t} < 0.065$ fm^-3 and 0.01 MeV/fm^3 < P_{t} < 0.26 $MeV/fm^3, respectively. These constrained values for the transition density and pressure are significantly lower than their fiducial values currently used in the literature. Furthermore, the mass-radius relation and several other properties closely related to the neutron star crust are studied by using the MDI interaction. It is found that the newly constrained$ ρ_t $and$ P_t $together with the earlier estimate of$ ΔI/I>0.014 $for the crustal fraction of the moment of inertia of the Vela pulsar impose a stringent constraint of R>= 4.7+4.0M/M_sun km for the radius$ R $and mass$ M$ of neutron stars.

Nuclear constraints on the inner edge of neutron star crusts

We show that the widely used parabolic approximation to the Equation of State (EOS) of asymmetric nuclear matter leads systematically to significantly higher core-crust transition densities and pressures. Using an EOS for neutron-rich nuclear matter constrained by the isospin diffusion data from heavy-ion reactions in the same sub-saturation density range as the neutron star crust, the density and pressure at the inner edge separating the liquid core from the solid crust of neutron stars are determined to be 0.040 fm $^{-3}$ $\leq ρ_{t}\leq 0.065$ fm $^{-3}$ and 0.01 MeV/fm $^{3}$ $\leq P_{t}\leq 0.26$ MeV/fm $^{3}$ , respectively. Implications of these constraints on the Vela pulsar are discussed.

Super-soft symmetry energy encountering non-Newtonian gravity in neutron stars

Considering the non-Newtonian gravity proposed in the grand unification theories, we show that the stability and observed global properties of neutron stars can not rule out the super-soft nuclear symmetry energies at supra-saturation densities. The degree of possible violation of the Inverse-Square-Law of gravity in neutron stars is estimated using an Equation of State (EOS) of neutron-rich nuclear matter consistent with the available terrestrial laboratory data.

Transition Density and Pressure at the Inner Edge of Neutron Star Crusts

Using the nuclear symmetry energy that has been recently constrained by the isospin diffusion data in intermediate-energy heavy ion collisions, we have studied the transition density and pressure at the inner edge of neutron star crusts, and they are found to be 0.040 fm $^{-3}$ $\leq ρ_{t}\leq 0.065$ fm $^{-3}$ and 0.01 MeV/fm $^{3}$ $\leq P_{t}\leq 0.26$ MeV/fm $^{3}$ , respectively, in both the dynamical and thermodynamical approaches. We have also found that the widely used parabolic approximation to the equation of state of asymmetric nuclear matter gives significantly higher values of core-crust transition density and pressure, especially for stiff symmetry energies. With these newly determined transition density and pressure, we have obtained an improved relation between the mass and radius of neutron stars.