BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Bao-An Li

Bao-An Li appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
nucl-thnucl-exastro-ph.SRastro-ph.HEhep-phgr-qcastro-phastro-ph.COastro-ph.GAhep-th

Catalog footprint

What is connected

106works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

106 published item(s)

preprint2026arXiv

A Universal Upper Bound on the Pressure-to-Energy Density Ratio in Neutron Stars

The equation-of-state (EOS) parameter $ϕ\equiv P/\varepsilon$, defined as the ratio of pressure to energy density, encapsulates the fundamental response of matter under extreme compression. Its value at the center of the most massive neutron star (NS), $\x \equiv ϕ_{\rm c} = P_{\rm c}/\varepsilon_{\rm c}$, sets a universal upper bound on the maximum denseness attainable by any form of visible matter anywhere in the Universe. Remarkably, owing to the intrinsically nonlinear structure of the EOS in General Relativity (GR), this bound is forced to lie far below the naive Special Relativity (SR) limit of unity. In this work, we refine the theoretical upper bound on $\x$ in a self-consistent manner by incorporating, in addition to the causality constraint from SR, the mass-sphere stability condition associated with the mass evolution pattern in the vicinity of the NS center. This condition is formulated within the intrinsic-and-perturbative analysis of the dimensionless Tolman--Oppenheimer--Volkoff equations (IPAD-TOV) framework. The combined constraints yield an improved bound, $\x \lesssim 0.385$, which is slightly above but fully consistent with the previously derived causal-only limit, $\x \lesssim 0.374$. We further derive an improved scaling relation for NS compactness and verify its universality across a broad set of 284 realistic EOSs, including models with first-order phase transitions, exotic degrees of freedom, continuous crossover behavior, and deconfined quark cores. The resulting bound on $\x$ thus provides a new, EOS-independent window into the microphysics of cold superdense matter compressed by strong-field gravity in GR.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Is The Trace Anomaly at its Minimum Value at Neutron Star Centers?

While the equation of state (EOS) $P(\varepsilon)$ of neutron star (NS) matter has been extensively studied, the EOS-parameter $ϕ= P/\varepsilon$ or equivalently the dimensionless trace anomaly $Δ= 1/3 - ϕ$, which quantifies the balance between pressure $P$ and energy density $\varepsilon$, remains far less explored, especially in NS cores. Its bounds and density profile carry crucial information about the nature of superdense matter. Physically, the EOS-parameter $ϕ$ represents the mean stiffness of matter accumulated from the stellar surface up to a given density. Based on the intrinsic structure of the Tolman--Oppenheimer--Volkoff equations, we show that $ϕ$ decreases monotonically outward from the NS center, independent of any specific input NS EOS model. Furthermore, observational evidence of a peak in the speed-of-sound squared (SSS) density-profile near the center effectively rules out a valley and a subsequent peak in the radial profile of $ϕ$ at similar densities, reinforcing its monotonic decrease. These model-independent relations impose strong constraints on the near-center behavior of the EOS-parameter $ϕ$, particularly demonstrating that the mean stiffness (or equivalently $Δ$) reaches a local maximum (minimum) at the center.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Equation of State of Neutron-Rich Matter in $d$-Dimensions

Nuclear systems under constraints, with high degrees of symmetries and/or collectivities may be considered as moving effectively in spaces with reduced spatial dimensions. We first derive analytical expressions for the nucleon specific energy $E_0(ρ)$, pressure $P_0(ρ)$, incompressibility coefficient $K_0(ρ)$ and skewness coefficient $J_0(ρ)$ of symmetric nucleonic matter (SNM), the quadratic symmetry energy $E_{\rm{sym}}(ρ)$, its slope parameter $L(ρ)$ and curvature coefficient $K_{\rm{sym}}(ρ)$ as well as the fourth-order symmetry energy $E_{\rm{sym,4}}(ρ)$ of neutron-rich matter in general $d$ spatial dimensions (abbreviated as "$d$D") in terms of the isoscalar and isovector parts of the isospin-dependent single-nucleon potential according to the generalized Hugenholtz-Van Hove (HVH) theorem. The equation of state (EOS) of nuclear matter in $d$D can be linked to that in the conventional 3-dimensional (3D) space by the $ε$-expansion which is a perturbative approach successfully used previously in treating second-order phase transitions and related critical phenomena and more recently in studying the EOS of cold atoms. The $ε$-expansion of nuclear EOS in $d$D based on a reference dimension $d_{\rm{f}}=d-ε$ is shown to be effective with $-1\lesssimε\lesssim1$ starting from $1\lesssim d_{\rm{f}}\lesssim3$ in comparison with the exact expressions derived using the HVH theorem. Moreover, the EOS of SNM (with/without considering its potential part) is found to be reduced (enhanced) in lower (higher) dimensions, indicating in particular that the many-nucleon system tends to be deeper bounded but saturate at higher densities in spaces with lower dimensions. The links between the EOSs in 3D and $d$D spaces from the $ε$-expansion provide new perspectives to the EOS of neutron-rich matter.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Investigating Effects of Relativistic Kinematics, Dimensionality, Interactions, and Short-Range Correlations on the Ratio of Quartic over Quadratic Nuclear Symmetry Energies

While ample evidence for the so-called empirical parabolic law of the Equation of State (EOS) of isospin asymmetric nuclear matter (ANM) has been obtained in many studies within both non-relativistic and relativistic nuclear many-body theories using various interactions, it has been unclear if there is any fundamental physics reason for the small quartic symmetry energy compared to the quadratic one even as the ANM approaches pure neutron matter. Within both relativistic and non-relativistic Free Fermi Gas (FFG) models in coordinate spaces of arbitrary dimension $d$ with and without considering Short-Range Correlations (SRC) as well as non-linear Relativistic Mean Field (RMF) models, we study effects of relativistic kinematics, dimensionality, interactions and SRC on the ratio $Ψ(ρ)$ of quartic over quadratic symmetry energies in ANM EOSs. We found that the ratio $Ψ(ρ)$ in the FFG model depends strongly on the dimension $d$. While it is very small already in the normal 3D space, it could be even smaller in spaces with reduced dimensions for sub-systems of particles in heavy-ion reactions and/or whole neutron stars due to constraints, collectivities and/or symmetries. We also found that the ratio $Ψ(ρ)$ could theoretically become very large only at the ultra-relativistic limit far above the density reachable in neutron stars. On the other hand, nuclear interaction directly and/or indirectly through SRC-induced high-momentum nucleons affect significantly the density dependence of $Ψ(ρ)$ compared to the relativistic FFG model prediction. The SRC affects significantly not only the kinetic energy of symmetric nuclear matter but also the ratio $Ψ(ρ)$ while the relativistic corrections are found negligible. The results may help better understand the EOS of dense neutron-rich matter.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Probing high-density nuclear symmetry energy with $Ξ^{-}/Ξ^{0}$ ratio in heavy-ion collisions at $\sqrt{s_{NN}} \sim 3$ GeV

Recent beam energy scan (BES) experiments at RHIC by the STAR Collaboration (PLB {\bf 827},137003 (2022) and PRL {\bf 128}, 202303 (2022)) found that hadronic interactions dominate the collective flow and the proton cumulant ratios are driven by baryon number conservation in a region of high baryon density in $\sqrt{s_{NN}}$ = 3 GeV Au+Au reactions, indicating the dense medium formed in such collisions is likely hadronic matter. Within an updated ART (A Relativistic Transport) model with momentum dependent isoscalar and isovector single-nucleon mean-field potentials corresponding to different symmetry energies at suprasaturation densities, the $n/p$, $π^{-}/π^{+}$, $K_{s}^{0}/K^{+}$, $Σ^{-}/Σ^{+}$ and $Ξ^{-}/Ξ^{0}$ ratios are studied for central Au+Au collisions at $\sqrt{s_{NN}}$ = 3 GeV where the maximum central density reaches about $(3.6\sim 4.0)ρ_0$. The doubly strange $Ξ^{-}/Ξ^{0}$ ratio is found to have the strongest sensitivity to the variation of high-density nuclear symmetry energy. Thus, the $Ξ^{-}/Ξ^{0}$ ratio in relativistic heavy-ion reactions at $\sqrt{s_{NN}} \sim 3$ GeV may help probe sensitively the poorly known symmetry energy of dense neutron-rich matter critically important for understanding various properties of neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Bayesian inference of dense matter EOS encapsulating a first-order hadron-quark phase transition from observables of canonical neutron stars

[Purpose:] We infer the posterior probability distribution functions (PDFs) and correlations of nine parameters characterizing the EOS of dense neutron-rich matter encapsulating a first-order hadron-quark phase transition from the radius data of canonical NSs reported by LIGO/VIRGO, NICER and Chandra Collaborations. We also infer the quark matter (QM) mass fraction and its radius in a 1.4 M$_{\odot}$ NS and predict their values in more massive NSs. [Method:] Meta-modelings are used to generate both hadronic and QM EOSs in the Markov-Chain Monte Carlo sampling process within the Bayesian statistical framework. An explicitly isospin-dependent parametric EOS for the $npeμ$ matter in NSs at $β$ equilibrium is connected through the Maxwell construction to the QM EOS described by the constant speed of sound (CSS) model of Alford, Han and Prakash. [Results:] (1) The most probable values of the hadron-quark transition density $ρ_t/ρ_0$ and the relative energy density jump there $\De\ep/\ep_t$ are $ρ_t/ρ_0=1.6^{+1.2}_{-0.4}$ and $\De\ep/\ep_t=0.4^{+0.20}_{-0.15}$ at 68\% confidence level, respectively. The corresponding probability distribution of QM fraction in a 1.4 M$_{\odot}$ NS peaks around 0.9 in a 10 km sphere. Strongly correlated to the PDFs of $ρ_t$ and $\De\ep/\ep_t$, the PDF of the QM speed of sound squared $\cQMsq/c^2$ peaks at $0.95^{+0.05}_{-0.35}$, and the total probability of being less than 1/3 is very small. (2) The correlations between PDFs of hadronic and QM EOS parameters are very weak. [Conclusions:] The available astrophysical data considered together with all known EOS constraints from theories and terrestrial nuclear experiments prefer the formation of a large volume of QM even in canonical NSs.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Tasting Nuclear Pasta Made with Classical Molecular Dynamics Simulations

Nuclear clusters or voids in the inner crust of neutron stars were predicted to have various shapes collectively nicknamed nuclear pasta. The recent review in Ref. \cite{Lopez1} by López, Dorso and Frank summarized their systematic investigations into properties especially the morphological and thermodynamical phase transitions of the nuclear pasta within a Classical Molecular Dynamics model, providing further stimuli to find more observational evidences of the predicted nuclear pasta in neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bayesian inference of the incompressibility, skewness and kurtosis of nuclear matter from empirical pressures in relativistic heavy-ion collisions

Within the Bayesian statistical framework we infer the incompressibility $K_0$, skewness $J_0$ and kurtosis $Z_0$ parameters of symmetric nuclear matter (SNM) at its saturation density $ρ_0$ using the constraining bands on the pressure in cold SNM in the density range of 1.3$ρ_0$ to 4.5$ρ_0$ from transport model analyses of kaon production and nuclear collective flow in relativistic heavy-ion collisions. As the default option assuming the $K_0$, $J_0$ and $Z_0$ have Gaussian prior probability distribution functions (PDFs) with the means and variances of $235\pm 30$, $-200\pm 200$ and $-146\pm 1728$ MeV, their posterior most probable values are narrowed down to 192$^{+12}_{-16}$ MeV, -180$^{+100}_{-110}$ MeV and 200$^{+250}_{-250}$ at 68\% confidence level, respectively. The results are largely independent of the prior PDFs of $J_0$ and $Z_0$ used. However, if one adopts the strong belief that the incompressibility $K_0$ has a uniform prior PDF within its absolute boundary of 220-260 MeV as one can find easily in the literature, the posterior most probable values of $K_0$, $J_0$ and $Z_0$ shift to $K_0=220^{+6}_{-0}$ MeV, $J_0=-390^{+60}_{-70}$ MeV and $Z_0=600^{+200}_{-200}$ MeV, respectively. While the posterior PDFs of the SNM EOS parameters depend somewhat on the prior PDF of $K_0$ used, the results from using different prior PDFs are qualitatively consistent. The uncertainties of all three parameters are significantly reduced especially for the $J_0$ and $Z_0$ parameters compared to their current values.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bayesian Inference of the Symmetry Energy of Super-Dense Neutron-Rich Matter from Future Radius Measurements of Massive Neutron Stars

Using an explicitly isospin-dependent parametric Equation of State (EOS) for the core of neutron stars (NSs) within the Bayesian statistical approach, we infer the EOS parameters of super-dense neutron-rich nuclear matter from three sets of imagined mass-radius correlation data representing typical predictions by various nuclear many-body theories, i.e, the radius stays the same, decreases or increases with increasing NS mass within $\pm 15\%$ between 1.4 M$_{\odot}$ and 2.0 M$_{\odot}$. The corresponding average density increases quickly, slowly or slightly decreases as the NS mass increases from 1.4 M$_{\odot}$ to 2.0 M$_{\odot}$. Using the posterior probability distribution functions (PDFs) of EOS parameters inferred from GW170817 and NICER radius data for canonical NSs as references, we investigate how future radius measurements of massive NS will improve our knowledge about the EOS of super-dense neutron-rich nuclear matter, especially its symmetry energy term, compared to what people have already learned from analyzing the GW170817 and NICER data. While the EOS of symmetric nuclear matter (SNM) inferred from the three data sets are approximately the same, the corresponding high-density symmetry energies at densities above about $2ρ_0$ are very different, indicating that the radii of massive NSs carry reliable information about the high-density behavior of nuclear symmetry energy with little influence from the remaining uncertainties of the SNM EOS.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraints on the muon fraction and density profile in neutron stars

Muons in neutron stars (NSs) play especially important roles in addressing several interesting new physics questions associated with detecting as well as understanding interactions and astrophysical effects of muonphilic dark matter particles. The key model inputs for studying the latter are the total muon mass $M_μ$, the muon mass fraction $M_μ/M_{\rm NS}$ over the NS mass $M_{\rm NS}$ and the muon radial density profile $ρ_μ(r)$ in NSs of varying masses. We investigate these quantities within a minimum model for the core of NSs consisting of neutrons, protons, electrons, and muons using an explicitly isospin-dependent parametric Equation of State (EOS) constrained by available nuclear laboratory experiments and the latest astrophysical observations of NS masses, radii and tidal deformabilities. We found that the absolutely maximum muon mass $M_μ$ and its mass fraction $M_μ/M_{\rm NS}$ in the most massive NSs allowed by causality are about 0.025 $M_\odot$ and 1.1\%, respectively. For the most massive NS of mass 2.14 $M_\odot$ observed so far, they reduce to about 0.020 $M_\odot$ and 0.9\%, respectively. We also study respective effects of individual parameters describing the EOS of high-density neutron-rich nucleonic matter on the muon contents in NSs with varying masses. We found that the most important but uncertain nuclear physics ingredient for determining the muon contents in NSs is the high-density nuclear symmetry energy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Symmetry energy of super-dense neutron-rich matter from integrating barotropic pressures in neutron stars and heavy-ion reactions

Within the minimum model of neutron stars (NS) consisting of neutrons, protons and electrons, a new approach is proposed for inferring the symmetry energy of super-dense neutron-rich nucleonic matter above twice the saturation density $ρ_0$ of nuclear matter directly from integrating iteratively barotropic pressures in both neutron stars and heavy-ion reactions. Simultaneously, the proton fraction of NSs at $β$ equilibrium is extracted as a function of baryon density from the same procedure. An application of this approach using the NS pressure from GW170817 and the pressure in cold symmetric nuclear matter (SNM) extracted earlier by analyzing nuclear collective flow data in relativistic heavy-ion collisions provides a useful constraining band for the symmetry energy above $2ρ_0$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Equations of motion of test particles for solving the spin-dependent Boltzmann-Vlasov equation

A consistent derivation of the equations of motion (EOMs) of test particles for solving the spin-dependent Boltzmann-Vlasov equation is presented. The resulting EOMs in phase space are similar to the canonical equations in Hamiltonian dynamics, and the EOM of spin is the same as that in the Heisenburg picture of quantum mechanics. Considering further the quantum nature of spin and choosing the direction of total angular momentum in heavy-ion reactions as a reference of measuring nucleon spin, the EOMs of spin-up and spin-down nucleons are given separately. The key elements affecting the spin dynamics in heavy-ion collisions are identified. The resulting EOMs provide a solid foundation for using the test-particle approach in studying spin dynamics in heavy-ion collisions at intermediate energies. Future comparisons of model simulations with experimental data will help constrain the poorly known in-medium nucleon spin-orbit coupling relevant for understanding properties of rare isotopes and their astrophysical impacts.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nuclear constraints on gravitational waves from deformed pulsars

The recent direct detection of gravitational waves (GWs) from binary black hole mergers (2016, Phys. Rev. Lett. 116, no. 6, 061102; no. 24, 241103) opens up an entirely new non-electromagnetic window into the Universe making it possible to probe physics that has been hidden or dark to electromagnetic observations. In addition to cataclysmic events involving black holes, GWs can be triggered by physical processes and systems involving neutron stars. Properties of neutron stars are largely determined by the equation of state (EOS) of neutron-rich matter, which is the major ingredient in calculating the stellar structure and properties of related phenomena, such as gravitational wave emission from elliptically deformed pulsars and neutron star binaries. Although the EOS of neutron-rich matter is still rather uncertain mainly due to the poorly known density dependence of nuclear symmetry energy at high densities, significant progress has been made recently in constraining the symmetry energy using data from terrestrial nuclear laboratories. These constraints could provide useful information on the limits of GWs expected from neutron stars. Here after briefly reviewing our previous work on constraining gravitational radiation from elliptically deformed pulsars with terrestrial nuclear laboratory data in light of the recent gravitational wave detection, we estimate the maximum gravitational wave strain amplitude, using an optimistic value for the breaking strain of the neutron star crust, for 15 pulsars at distances 0.16 kpc to 0.91 kpc from Earth, and find it to be in the range of $\sim[0.2-31.1]\times 10^{-24}$, depending on the details of the EOS used to compute the neutron star properties. Implications are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nucleon Effective E-Mass in Neutron-Rich Matter from the Migdal-Luttinger Jump

The well-known Migdal-Luttinger theorem states that the jump of the single-nucleon momentum distribution at the Fermi surface is equal to the inverse of the nucleon effective E-mass. Recent experiments studying short-range correlations (SRC) in nuclei using electron-nucleus scatterings at the Jefferson National Laboratory (JLAB) together with model calculations constrained significantly the Migdal-Luttinger jump at saturation density of nuclear matter. We show that the corresponding nucleon effective E-mass is consequently constrained to $M_0^{\ast,\rm{E}}/M\approx2.22\pm0.35$ in symmetric nuclear matter (SNM) and the E-mass of neutrons is smaller than that of protons in neutron-rich matter. Moreover, the average depletion of the nucleon Fermi sea increases (decreases) approximately linearly with the isospin asymmetry $δ$ according to $κ_{\rm{p}/\rm{n}}\approx 0.21\pm0.06 \pm (0.19\pm0.08)δ$ for protons (neutrons). These results will help improve our knowledge about the space-time non-locality of the single-nucleon potential in neutron-rich nucleonic matter useful in both nuclear physics and astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spin transport in intermediate-energy heavy-ion collisions as a probe of in-medium spin-orbit interactions

The spin up-down splitting of collective flows in intermediate-energy heavy-ion collisions as a result of the nuclear spin-orbit interaction is investigated within a spin- and isospin-dependent Boltzmann-Uehing-Uhlenbeck transport model SIBUU12. Using a Skyrme-type spin-orbit coupling quadratic in momentum, we found that the spin splittings of the directed flow and elliptic flow are largest in peripheral Au+Au collisions at beam energies of about 100-200 MeV/nucleon, and the effect is considerable even in smaller systems especially for nucleons with high transverse momenta. The collective flows of light clusters of different spin states are also investigated using an improved dynamical coalescence model with spin. Our study can be important in understanding the properties of in-medium nuclear spin-orbit interactions once the spin-dependent observables proposed in this work can be measured.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Symmetry energy of cold nucleonic matter within a relativistic mean field model encapsulating effects of high momentum nucleons induced by short-range correlations

Significant progress has been made recently in constraining the isospin-dependent parameters characterizing the SRC (short-range correlation)-modified single-nucleon momentum distribution in neutron-rich nucleonic matter using both experimental data and microscopic model calculations. Using the constrained single-nucleon momentum distribution in a nonlinear relativistic mean field (RMF) model, we study the equation of state (EOS) of asymmetric nucleonic matter (ANM), especially the density dependence of nuclear symmetry energy $E_{\rm{sym}}(ρ)$. Firstly, as a test of the model, the average nucleon kinetic energy extracted recently from electron-nucleus scattering experiments using a neutron-proton dominance model is well reproduced by the RMF model incorporating effects of the SRC-induced high momentum nucleons, while it is significantly under predicted by the RMF model using a step function for the single-nucleon momentum distribution as in free Fermi gas (FFG) models. Secondly, the kinetic symmetry energy of quasi-nucleons is found to be $E^{\rm{kin}}_{\rm{sym}}(ρ_0)=-16.94\pm13.66\,\rm{MeV}$ which is dramatically different from the prediction of $E^{\rm{kin}}_{\rm{sym}}(ρ_0)\approx 12.5$ MeV by FFG models at nuclear matter saturation density $ρ_0=0.16\,\rm{fm}^{-3}$. Thirdly, comparing the RMF calculations with and without the high momentum nucleons using two sets of model parameters both reproducing identically all empirically constraints on the EOS of symmetric nuclear matter (SNM) and the symmetry energy of ANM at $ρ_0$, the SRC-modified single-nucleon momentum distribution is found to make the $E_{\rm{sym}}(ρ)$ more concave around $ρ_0$ by softening it significantly at both sub-saturation and supra-saturation densities, leading to an isospin-dependent incompressibility of ANM in better agreement with existing experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraints on the symmetry energy from observational probes of the neutron star crust

A number of observed phenomena associated with individual neutron star systems or neutron star populations find explanations in models in which the neutron star crust plays an important role. We review recent work examining the sensitivity to the slope of the symmetry energy $L$ of such models, and constraints extracted on $L$ from confronting them with observations. We focus on six sets of observations and proposed explanations: (i) The cooling rate of the neutron star in Cassiopeia A, confronting cooling models which include enhanced cooling in the nuclear pasta regions of the inner crust, (ii) the upper limit of the observed periods of young X-ray pulsars, confronting models of magnetic field decay in the crust caused by the high resistivity of the nuclear pasta layer, (iii) glitches from the Vela pulsar, confronting the paradigm that they arise due to a sudden re-coupling of the crustal neutron superfluid to the crustal lattice after a period during which they were decoupled due to vortex pinning, (iv) The frequencies of quasi-periodic oscillations in the X-ray tail of light curves from giant flares from soft gamma-ray repeaters, confronting models of torsional crust oscillations, (v) the upper limit on the frequency to which millisecond pulsars can be spun-up due to accretion from a binary companion, confronting models of the r-mode instability arising above a threshold frequency determined in part by the viscous dissipation timescale at the crust-core boundary, and (vi) the observations of precursor electromagnetic flares a few seconds before short gamma-ray bursts, confronting a model of crust shattering caused by resonant excitation of a crustal oscillation mode by the tidal gravitational field of a companion neutron star just before merger.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Critical Density and Impact of $Δ(1232)$ Resonance Formation in Neutron Stars

The critical densities and impact of forming \D resonances in neutron stars are investigated within an extended nonlinear relativistic mean-field (RMF) model. The critical densities for the formation of four different charge states of \D are found to depend differently on the separate kinetic and potential parts of nuclear symmetry energy, the first example of a microphysical property of neutron stars to do so. Moreover, they are sensitive to the in-medium Delta mass $m_Δ$ and the completely unknown $Δ$-$ρ$ coupling strength $g_{ρΔ}$. In the universal baryon-meson coupling scheme where the respective $Δ$-meson and nucleon-meson coupling constants are assumed to be the same, the critical density for the first $Δ^-(1232)$ to appear is found to be \rc=$(2.08\pm0.02)ρ_0$ using RMF model parameters consistent with current constraints on all seven macroscopic parameters usually used to characterize the equation of state (EoS) of isospin-asymmetric nuclear matter (ANM) at saturation density $ρ_0$. Moreover, the composition and the mass-radius relation of neutron stars are found to depend significantly on the values of the $g_{ρΔ}$ and $m_Δ$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamical effects of spin-dependent interactions in low- and intermediate-energy heavy-ion reactions

It is well known that non-central nuclear forces, such as the spin-orbital coupling and the tensor force, play important roles in understanding many interesting features of nuclear structures. However, their dynamical effects in nuclear reactions are poorly known since only the spin-averaged observables are normally studied both experimentally and theoretically. Realizing that spin-sensitive observables in nuclear reactions may carry useful information about the in-medium properties of non-central nuclear interactions, besides earlier studies using the time-dependent Hartree-Fock approach to understand effects of spin-orbital coupling on the threshold energy and spin polarization in fusion reactions, some efforts have been made recently to explore dynamical effects of non-central nuclear forces in intermediate-energy heavy-ion collisions using transport models. The focuses of these studies have been on investigating signatures of the density and isospin dependence of the form factor in the spin-dependent single-nucleon potential. Interestingly, some useful probes were identified in the model studies while so far there is still no data to compare with. In this brief review, we summarize the main physics motivations as well as the recent progress in understanding the spin dynamics and identifying spin-sensitive observables in heavy-ion reactions at intermediate energies. We hope the interesting, important, and new physics potentials identified in the spin dynamics of heavy-ion collisions will stimulate more experimental work in this direction.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effects of short-range correlation reduced kinetic symmetry energy in heavy-ion collisions at intermediate energies

Besides earlier predictions based on both phenomenological models and modern microscopic many-body theories, circumstantial evidence was recently found for a reduced kinetic symmetry energy of isospin-asymmetric nucleonic matter compared to the free Fermi gas model prediction due to the short-range correlation of high-momentum neutron-proton pairs. While keeping the total symmetry energy near the saturation density of nuclear matter consistent with existing experimental constraints, we examine the correspondingly enhanced role of the isospin degree of freedom in heavy-ion collisions at intermediate energies due to the reduced (enhanced) kinetic (potential) symmetry energy. Important observable consequences are investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Efficacy of crustal superfluid neutrons in pulsar glitch models

In order to assess the ability of purely crust-driven glitch models to match the observed glitch activity in the Vela pulsar, we conduct a systematic analysis of the dependence of the fractional moment of inertia of the inner crustal neutrons on the stiffness of the nuclear symmetry energy at saturation density $L$. We take into account both crustal entrainment and the fact that only a fraction $Y_{\rm g}$ of the core neutrons may couple to the crust on the glitch-rise timescale. We use a set of consistently-generated crust and core compositions and equations-of-state which are fit to results of low-density pure neutron matter calculations. When entrainment is included at the level suggested by recent microscopic calculations and the core is fully coupled to the crust, the model is only able to account for the Vela glitch activity for a 1.4$M_{\odot}$ star if the equation of state is particularly stiff $L>100$ MeV. However, an uncertainty of about 10\% in the crust-core transition density and pressure allows for the Vela glitch activity to be marginally accounted for in the range $L\approx30-60$MeV consistent with a range of experimental results. Alternatively, only a small amount of core neutrons need be involved. If less than 50\% of the core neutrons are coupled to the crust during the glitch, we can also account for the Vela glitch activity using crustal neutrons alone for EOSs consistent with the inferred range of $L$. We also explore the possibility of Vela being a high-mass neutron star, and of crustal entrainment being reduced or enhanced relative to its currently predicted values.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Shedding Light on the EOS-Gravity Degeneracy and Constraining the Nuclear Symmetry Energy from the Gravitational Binding Energy of Neutron Stars

A thorough understanding of properties of neutron stars requires both a reliable knowledge of the equation of state (EOS) of super-dense nuclear matter and the strong-field gravity theories simultaneously. To provide information that may help break this EOS-gravity degeneracy, we investigate effects of nuclear symmetry energy on the gravitational binding energy of neutron stars within GR and the scalar-tensor subset of alternative gravity models. We focus on effects of the slope $L$ of nuclear symmetry energy at saturation density and the high-density behavior of nuclear symmetry energy. We find that the variation of either the density slope $L$ or the high-density behavior of nuclear symmetry energy leads to large changes in the binding energy of neutron stars. The difference in predictions using the GR and the scalar-tensor theory appears only for massive neutron stars, and even then is significantly smaller than the difference resulting from variations in the symmetry energy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Small radii of neutron stars as an indication of novel in-medium effects

At present, neutron star radii from both observations and model predictions remain very uncertain. Whereas different models can predict a wide range of neutron star radii, it is not possible for most models to predict radii that are smaller than about 10 km, thus if such small radii are established in the future they will be very difficult to reconcile with model estimates. By invoking a new term in the equation of state that enhances the energy density, but leaves the pressure unchanged we simulate the current uncertainty in the neutron star radii. This new term can be possibly due to the exchange of the weakly interacting light U-boson with appropriate in-medium parameters, which does not compromise the success of the conventional nuclear models. The validity of this new scheme will be tested eventually by more precise measurements of neutron star radii.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Symmetry Energy of Nucleonic Matter With Tensor Correlations

The nuclear symmetry energy (E_sym(\roh)) is a vital ingredient of our understanding of many processes, from heavy-ion collisions to neutron stars structure. While the total nuclear symmetry energy at nuclear saturation density (ρ_0) is relatively well determined, its value at supranuclear densities is not. The latter can be better constrained by separately examining its kinetic and potential terms and their density dependencies. The kinetic term of the symmetry energy, E_sym^kin(ρ_0), equals the difference in the per-nucleon kinetic energy between pure neutron matter (PNM) and symmetric nuclear matter (SNM), often calculated using a simple Fermi gas model. However, experiments show that tensor force induced short-range correlations (SRC) between proton-neutron pairs shift nucleons to high-momentum in SNM, but have almost no effect in PNM. We present an approximate analytical expression for E_sym^kin(ρ) of correlated nucleonic matter. In our model, E_sym^kin(ρ_0) = -10 MeV, which differs significantly from +12.5 MeV for the widely-used free Fermi gas model. This result is consistent with our analysis of recent data on the free proton-to-neutron ratios measured in intermediate energy nucleus-nucleus collisions as well as with microscopic many-body calculations, and previous phenomenological extractions. We then use our calculated E_sym^kin(ρ) in combination with the known total symmetry energy and its density dependence at saturation density to constrain the value and density dependence of the potential part and to extrapolate the total symmetry energy to supranuclear densities.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Symmetry potential of $Δ(1232)$ resonance and its effects on the $π^-/π^+$ ratio in heavy-ion collisions near the pion production threshold

Effects of the completely unknown symmetry (isovector) potential of the \D on the total and differential \rpi in heavy-ion collisions at beam energies from 100 to 1000 MeV/A are explored within an isospin-dependent transport model IBUU. The effects are found to be negligible at beam energies above the pion production threshold due to the very short lifetimes of less than 2 fm/c for $Δ$ resonances with masses around $m_Δ=1232$ MeV, leaving the $π^-/π^+$ ratios of especially the energetic pions still a reliable probe of the high-density behavior of nuclear symmetry energy $E_{sym}(ρ)$. However, as the beam energy becomes deeply sub-threshold for pion production, effects of the $Δ$ symmetry potential becomes appreciable especially on the \rpi of low-energy pions from the decays of low-mass $Δ$ resonances which have lived long enough to be affected by their mean-field potentials, providing a useful tool to study the symmetry potential and spectroscopy of $Δ$ resonances in neutron-rich nuclear matter. Interestingly though, even at the deeply sub-threshold beam energies, the differential \rpi of energetic pions remains sensitive to the \esym at supra saturation densities with little influence from the uncertain symmetry potential of the $Δ$ resonance.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The isospin quartic term in the kinetic energy of neutron-rich nucleonic matter

The energy of a free gas of neutrons and protons is well known to be approximately isospin parabolic with a negligibly small quartic term of only $0.45$ MeV at the saturation density of nuclear matter $ρ_0=0.16/\rm{fm}^3$. Using an isospin-dependent single-nucleon momentum distribution including a high (low) momentum tail (depletion) with its shape parameters constrained by recent high-energy electron scattering and medium-energy nuclear photodisintegration experiments as well as the state-of-the-art calculations of the deuteron wave function and the equation of state of pure neutron matter near the unitary limit within several modern microscopic many-body theories, we show for the first time that the kinetic energy of interacting nucleons in neutron-rich nucleonic matter has a significant quartic term of $7.18\pm2.52\,\rm{MeV}$. Such a large quartic term has significant ramifications in determining the equation of state of neutron-rich nucleonic matter using both terrestrial and astrophysical observables.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Impact of the equation-of-state -- gravity degeneracy on constraining the nuclear symmetry energy from astrophysical observables

There is a degeneracy between the equation of state (EOS) of superdense neutron-rich nuclear matter and the strong-field gravity in understanding properties of neutron stars. While the EOS is still poorly known, there are also longstanding ambiguities in choosing General Relativity or alternative gravity theories in the not-so-well tested strong-field regime. Besides possible appearance of hyperons and new phases, the most uncertain part of the nucleonic EOS is currently the density dependence of nuclear symmetry energy. To provide information that may help break the EOS-gravity degeneracy, we investigate effects of symmetry energy within its uncertain range determined by terrestrial nuclear laboratory experiments on the gravitational binding energy and spacetime curvature of neutron stars within GR and the scalar-tensor (ST) theory of gravity. In particular, we focus on effects of the following parameters characterizing the EOS of neutron-rich nucleonic matter: (1) the incompressibility $K_0$ of symmetric nuclear matter, (2) the slope $L$ of symmetry energy at saturation density and (3) the high-density behavior of symmetry energy. We find that the variation of either the density slope $L$ or the high-density behavior of symmetry energy leads to large changes in both the binding energy and curvature of neutron stars while effects of varying the more constrained $K_0$ are negligibly small. The difference in predictions using the GR and the ST theory appears only for massive neutron stars, and is significantly smaller than the differences resulting from variations in the symmetry energy. We conclude that within the ST theory of gravity, the EOS-gravity degeneracy has been broken by the recent relativistic pulsar measurements, and that measurements of neutron star properties sensitive to the compactness constrain mainly the density dependence of the symmetry energy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nuclear Theory and Science of the Facility for Rare Isotope Beams

The Facility for Rare Isotope Beams (FRIB) will be a world-leading laboratory for the study of nuclear structure, reactions and astrophysics. Experiments with intense beams of rare isotopes produced at FRIB will guide us toward a comprehensive description of nuclei, elucidate the origin of the elements in the cosmos, help provide an understanding of matter in neutron stars, and establish the scientific foundation for innovative applications of nuclear science to society. FRIB will be essential for gaining access to key regions of the nuclear chart, where the measured nuclear properties will challenge established concepts, and highlight shortcomings and needed modifications to current theory. Conversely, nuclear theory will play a critical role in providing the intellectual framework for the science at FRIB, and will provide invaluable guidance to FRIB's experimental programs. This article overviews the broad scope of the FRIB theory effort, which reaches beyond the traditional fields of nuclear structure and reactions, and nuclear astrophysics, to explore exciting interdisciplinary boundaries with other areas. \keywords{Nuclear Structure and Reactions. Nuclear Astrophysics. Fundamental Interactions. High Performance Computing. Rare Isotopes. Radioactive Beams.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantifying Correlations Between Isovector Observables and the Density Dependence of Nuclear Symmetry Energy away from Saturation Density

According to the Hugenholtz-Van Hove theorem, the nuclear symmetry energy $S(ρ)$ and its slope $L(ρ)$ at arbitrary densities can be decomposed in terms of the density and momentum dependence of the single-nucleon potentials in isospin-asymmetric nuclear matter which are potentially accessible to experiment. We quantify the correlations between several well-known isovector observables and $L(ρ)$ to locate the density range in which each isovector observable is most sensitive to the density dependence of the $S(ρ)$. We then study the correlation coefficients between those isovector observables and all the components of the $L(ρ)$. The neutron skin thickness of $^{208}$Pb is found to be strongly correlated with the $L(ρ)$ at a subsaturation density of $ρ= 0.59 ρ_0$ through the density dependence of the first-order symmetry potential. Neutron star radii are found to be strongly correlated with the $L(ρ)$ over a wide range of supra-saturation densities mainly through both the density and momentum dependence of the first-order symmetry potential. Finally, we find that although the crust-core transition pressure has a complex correlation with the $L(ρ)$, it is strongly correlated with the momentum derivative of the first-order symmetry potential, and the density dependence of the second-order symmetry potential.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spin-orbit coupling and the up-down differential transverse flow in intermediate-energy heavy-ion collisions

To explore the strength, isospin dependence, and density dependence of the in-medium spin-orbit coupling that is poorly known but relevant for understanding the structure of rare isotopes, the nucleon spin up-down differential transverse flow is studied systematically by varying the beam energy and impact parameter of Au+Au collisions within a spin-isospin dependent transport model recently developed for investigating spin-related phenomena in intermediate-energy heavy-ion collisions. The optimal reaction conditions for delineating and disentangling effects of different terms in the spin-orbit coupling are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Constraining the neutron-proton effective mass splitting using empirical constraints on the density dependence of nuclear symmetry energy around normal density

According to the Hugenholtz-Van Hove theorem, nuclear symmetry energy \esym and its slope \lr at an arbitrary density $ρ$ are determined by the nucleon isovector (symmetry) potential \usym and its momentum dependence $\frac{\partial U_{sym}}{\partial k}$. The latter determines uniquely the neutron-proton effective k-mass splitting $m^*_{n-p}(ρ,δ)\equiv (m_{\rm n}^*-m_{\rm p}^*)/m$ in neutron-rich nucleonic matter of isospin asymmetry $δ$. Using currently available constraints on the \es0 and \l0 at normal density $ρ_0$ of nuclear matter from 28 recent analyses of various terrestrial nuclear laboratory experiments and astrophysical observations, we try to infer the corresponding neutron-proton effective k-mass splitting $m^*_{n-p}(ρ_0,δ)$. While the mean values of the $m^*_{n-p}(ρ_0,δ)$ obtained from most of the studies are remarkably consistent with each other and scatter very closely around an empirical value of \emass$=0.27\cdotδ$, it is currently not possible to scientifically state surely that the \emass is positive within the present knowledge of the uncertainties. Quantifying, better understanding and then further reducing the uncertainties using modern statistical and computational techniques in extracting the \es0 and \l0 from analyzing the experimental data are much needed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effects of nuclear symmetry energy on eta meson production and its rare decay to the dark U-boson in heavy-ion reactions

Using a relativistic transport model ART1.0, we explore effects of nuclear symmetry energy on eta meson production and its rare decay to the dark U-boson in heavy-ion reactions from 0.2 to 10 GeV/nucleon available at several current and future facilities. The yield of eta mesons at sub-threshold energies is found to be very sensitive to the density dependence of nuclear symmetry energy. Above a beam energy of about 5 GeV/nucleon in Au+Au reactions, the sensitivity to symmetry energy disappears. Using the branching ratio of the rare eta decay (eta->gamma U) available in the literature, we estimate the maximum cross section for the U-boson production in the energy range considered, providing a useful reference for future U-boson search using heavy-ion reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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}$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Probing in-medium spin-orbit potential with intermediate-energy heavy-ion collisions

Incorporating for the first time both the spin and isospin degrees of freedom explicitly in transport model simulations of intermediate-energy heavy-ion collisions, we observe that a local spin polarization appears during collision process. Most interestingly, it is found that the nucleon spin up-down differential transverse flow is a sensitive probe of the spin-orbit interaction, providing a novel approach to probe both the density and isospin dependence of the in-medium spin-orbit coupling that is important for understanding the structure of rare isotopes and synthesis of superheavy elements.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Probing the high-density behavior of symmetry energy with gravitational waves

Gravitational wave (GW) astronomy opens up an entirely new window on the Universe to probe the equations of state (EOS) of neutron-rich matter. With the advent of next generation GW detectors, measuring the gravitational radiation from coalescing binary neutron star systems, mountains on rotating neutron stars, and stellar oscillation modes may become possible in the near future. Using a set of model EOSs satisfying the latest constraints from terrestrial nuclear experiments, state of the art nuclear many-body calculations of the pure neutron matter EOS, and astrophysical observations consistently, we study various GW signatures of the high-density behavior of the nuclear symmetry energy, which is considered among the most uncertain properties of dense neutron-rich nucleonic matter. In particular, we find the tidal polarizability of neutron stars, potentially measurable in binary systems just prior to merger, is more sensitive to the high density component of the nuclear symmetry energy than the symmetry energy at nuclear saturation density. We also find that the upper limit on the GW strain amplitude from elliptically deformed stars is very sensitive to the density dependence of the symmetry energy. This suggests that future developments in modeling of the neutron star crust, and direct gravitational wave signals from accreting binaries will provide a wealth of information on the EOS of neutron-rich matter. We also review the sensitivity of the $r$-mode instability window to the density dependence of the symmetry energy. Whereas models with larger values of the density slope of the symmetry energy at saturation seem to be disfavored by the current observational data, within a simple $r$-mode model, we point out that a subsequent softer behavior of the symmetry energy at high densities (hinted at by recent observational interpretations) could rule them in.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The cooling of the Cassiopeia A neutron star as a probe of the nuclear symmetry energy and nuclear pasta

X-ray observations of the neutron star in the Cas A supernova remnant over the past decade suggest the star is undergoing a rapid drop in surface temperature of $\approx$ $2-5.5\%$. One explanation suggests the rapid cooling is triggered by the onset of neutron superfluidity in the core of the star, causing enhanced neutrino emission from neutron Cooper pair breaking and formation (PBF). Using consistent neutron star crust and core equations of state (EOSs) and compositions, we explore the sensitivity of this interpretation to the density dependence of the symmetry energy $L$ of the EOS used, and to the presence of enhanced neutrino cooling in the bubble phases of crustal "nuclear pasta". Modeling cooling over a conservative range of neutron star masses and envelope compositions, we find $L\lesssim70$ MeV, competitive with terrestrial experimental constraints and other astrophysical observations. For masses near the most likely mass of $M\gtrsim 1.65 M_{\odot}$, the constraint becomes more restrictive $35\lesssim L\lesssim 55$ MeV. The inclusion of the bubble cooling processes decreases the cooling rate of the star during the PBF phase, matching the observed rate only when $L\lesssim45$ MeV, taking all masses into consideration, corresponding to neutron star radii $\lesssim 11$km.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A survey of the parameter space of the compressible liquid drop model as applied to the neutron star inner crust

We present a systematic survey the range of predictions of the neutron star inner crust composition, crust-core transition densities and pressures, and density range of the nuclear `pasta' phases at the bottom of the crust provided by the compressible liquid drop model in the light of current experimental and theoretical constraints on model parameters. Using a Skyrme-like model for nuclear matter, we construct baseline sequences of crust models by consistently varying the density dependence of the bulk symmetry energy at nuclear saturation density, $L$, under two conditions: (i) that the magnitude of the symmetry energy at saturation density $J$ is held constant, and (ii) $J$ correlates with $L$ under the constraint that the pure neutron matter (PNM) EoS satisfies the results of ab-initio calculations at low densities. Such baseline crust models facilitate consistent exploration of the $L$ dependence of crustal properties. The remaining surface energy and symmetric nuclear matter parameters are systematically varied around the baseline, and different functional forms of the PNM EoS at sub-saturation densities implemented, to estimate theoretical `error bars' for the baseline predictions. Inner crust composition and transition densities are shown to be most sensitive to the surface energy at very low proton fractions and to the behavior of the sub-saturation PNM EoS. Recent calculations of the energies of neutron drops suggest that the low-proton-fraction surface energy might be higher than predicted in Skyrme-like models, which our study suggests may result in a greatly reduced volume of pasta in the crust than conventionally predicted.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraining the High-Density Behavior of Nuclear Symmetry Energy with the Tidal Polarizability of Neutron Stars

Using a set of model equations of state satisfying the latest constraints from both terrestrial nuclear experiments and astrophysical observations as well as state-of-the-art nuclear many-body calculations of the pure neutron matter equation of state, the tidal polarizability of canonical neutron stars in coalescing binaries is found to be a very sensitive probe of the high-density behavior of nuclear symmetry energy which is among the most uncertain properties of dense neutron-rich nucleonic matter. Moreover, it changes less than $\pm 10%$ by varying various properties of symmetric nuclear matter and symmetry energy around the saturation density within their respective ranges of remaining uncertainty.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraints on the symmetry energy from neutron star observations

The modeling of many neutron star observables incorporates the microphysics of both the stellar crust and core, which is tied intimately to the properties of the nuclear matter equation of state (EoS). We explore the predictions of such models over the range of experimentally constrained nuclear matter parameters, focusing on the slope of the symmetry energy at nuclear saturation density $L$. We use a consistent model of the composition and EoS of neutron star crust and core matter to model the binding energy of pulsar B of the double pulsar system J0737-3039, the frequencies of torsional oscillations of the neutron star crust and the instability region for r-modes in the neutron star core damped by electron-electron viscosity at the crust-core interface. By confronting these models with observations, we illustrate the potential of astrophysical observables to offer constraints on poorly known nuclear matter parameters complementary to terrestrial experiments, and demonstrate that our models consistently predict $L<70$ MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Delineating effects of tensor force on the density dependence of nuclear symmetry energy

In this talk, we report results of our recent studies to delineate effects of the tensor force on the density dependence of nuclear symmetry energy within phenomenological models. The tensor force active in the isosinglet neutron-proton interaction channel leads to appreciable depletion/population of nucleons below/above the Fermi surface in the single-nucleon momentum distribution in cold symmetric nuclear matter (SNM). We found that as a consequence of the high momentum tail in SNM the kinetic part of the symmetry energy $E^{kin}_{sym}(ρ)$ is significantly below the well-known Fermi gas model prediction of approximately $12.5 (ρ/ρ_0)^{2/3}$. With about 15% nucleons in the high momentum tail as indicated by the recent experiments at J-Lab by the CLAS Collaboration, the $E^{kin}_{sym}(ρ)$ is negligibly small. It even becomes negative when more nucleons are in the high momentum tail in SNM. These features have recently been confirmed by three independent studies based on the state-of-the-art microscopic nuclear many-body theories. In addition, we also estimate the second-order tensor force contribution to the potential part of the symmetry energy. Implications of these findings in extracting information about nuclear symmetry energy from nuclear reactions are discussed briefly.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Disentangling effects of collision geometry and symmetry energy in U+U collisions

Effects of the collision geometry on experimental observables that are known to be sensitive to the high-density behavior of nuclear symmetry energy are examined in U+U collisions at 0.52 GeV/nucleon using an isospin- and momentum-dependent interaction within the framework of IBUU transport model. It is found that the neutron-proton differential flow in tip-tip collisions and the difference of neutron and proton elliptic flow in body-body collisions are more sensitive to the symmetry energy at supra-saturation densities compared with collisions of spherical nuclei of the same masses. In addition, the n/p ratio of pre-equilibrium nucleons is found to be slightly more sensitive to the symmetry energy in tip-tip collisions, and the collision geometry affects the π^-/π^+ ratio significantly.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Generic Constraints on the Relativistic Mean-Field and Skyrme-Hartree-Fock Models from the Pure Neutron Matter Equation of State

We study the nuclear symmetry energy S(rho) and related quantities of nuclear physics and nuclear astrophysics predicted generically by relativistic mean-field (RMF) and Skyrme-Hartree-Fock (SHF) models. We establish a simple prescription for preparing equivalent RMF and SHF parametrizations starting from a minimal set of empirical constraints on symmetric nuclear matter, nuclear binding energy and charge radii, enforcing equivalence of their Lorenz effective masses, and then using the pure neutron matter (PNM) equation of state (EoS) obtained from ab-initio calculations to optimize the pure isovector parameters in the RMF and SHF models. We find the resulting RMF and SHF parametrizations give broadly consistent predictions of the symmetry energy J and its slope parameter L at saturation density within a tight range of <~2 MeV and <~6 MeV respectively, but that clear model dependence shows up in the predictions of higher-order symmetry energy parameters, leading to important differences in (a) the slope of the correlation between J and L from the confidence ellipse, (b) the isospin-dependent part of the incompressibility of nuclear matter K_tau, (c) the symmetry energy at supra-saturation densities, and (d) the predicted neutron star radii. The model dependence can lead to about 1-2 km difference in predictions of the neutron star radius given identical predicted values of J, L and symmetric nuclear matter (SNM) saturation properties. Allowing the full freedom in the effective masses in both models leads to constraints of 30<~J<~31.5 MeV, 35<~L<~60 MeV, -330<~K_tau<~-216 MeV for the RMF model as a whole and 30<~J<~33 MeV, 28<~L<~65 MeV, -420<~K_tau<~-325 MeV for the SHF model as a whole. Notably, given PNM constraints, these results place RMF and SHF models as a whole at odds with some constraints on K_tau inferred from giant monopole resonance and neutron skin experimental results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pure Neutron Matter Constraints and Nuclear Symmetry Energy

In this review, we will discuss the results of our recent work to study the general optimization of the pure isovector parameters of the popular relativistic mean-field (RMF) and Skyrme-Hartree-Fock (SHF) nuclear energy-density functionals (EDFs), using constraints on the pure neutron matter (PNM) equation of state (EoS) from recent {\sl ab initio} calculations. By using RMF and SHF parameterizations that give equivalent predictions for ground-state properties of doubly magic nuclei and properties of symmetric nuclear matter (SNM) and PNM, we found that such optimization leads to broadly consistent symmetry energy $J$ and its slope parameter $L$ at saturation density within a tight range of $σ(J) < 2$ MeV and $σ(L) < 6$ MeV. We demonstrate that a clear model dependence shows up (a) in the curvature parameter of the symmetry energy $K_{\rm sym}$, (b) the symmetry energy at supra-saturation densities, and (c) the radius of neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sensitivity of the neutron star r-mode instability window to the density dependence of the nuclear symmetry energy

Using a simple model of a neutron star with a perfectly rigid crust constructed and a set of crust and core equations of state that span the range of nuclear experimental uncertainty in the density dependence of the symmetry energy from 25 MeV (soft EOS) to 115 MeV (stiff EOS), we calculate the instability window for the onset of the Chandrasekhar-Friedmann-Schutz (CFS) instability in r-mode oscillations for canonical neutron stars (1.4 M_{\odot}) and massive neutron stars (2.0 M_{\odot}). In these models the crust-core transition density, and thus crustal thickness, is calculated consistently with the core equation of state (EOS). For the canonical neutron star, the lower bound of the r-mode instability window is reduced in frequency by \approx150 Hz from the softest to the stiffest symmetry energy used, independent of mass and temperature. The instability window also drops by \approx 100 Hz independent of EOS when the mass is raised from 1.4 M_{\odot} to 2.0 M_{\odot}. Where temperature estimates are available, the observed neutron stars in low mass X-ray binaries (LMXBs) have frequencies below the instability window for the 1.4 M_{\odot} models, while some LMXBs fall within the instability window for 2.0 M_{\odot} stars if the symmetry energy is relatively stiff, indicating that a softer symmetry energy is more consistent with observations within this model. Thus we conclude that smaller values of L help stabilize neutron stars against runaway r-mode oscillations. The critical temperature, below which no star can reach the instability window without exceeding its Kepler frequency, varies by nearly an order of magnitude from soft to stiff symmetry energies. When the crust thickness and core EOS are treated consistently, a thicker crust corresponds to a lower critical temperature, the opposite result to previous studies in which the transition density was independent of the core EOS.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magnetic effects in heavy-ion collisions at intermediate energies

The time-evolution and space-distribution of internal electromagnetic fields in heavy-ion reactions at beam energies between 200 and 2000 MeV/nucleon are studied within an Isospin-dependent Boltzmann-Uhling-Uhlenbeck transport model IBUU11. While the magnetic field can reach about $7\times 10^{16}$ G which is significantly higher than the estimated surface magnetic field ($\sim 10^{15}$ G) of magnetars, it has almost no effect on nucleon observables as the Lorentz force is normally much weaker than the nuclear force. Very interestingly, however, the magnetic field generated by the projectile-like (target-like) spectator has a strong focusing/diverging effect on positive/negative pions at forward (backward) rapidities. Consequently, the differential $π^-/π^+$ ratio as a function of rapidity is significantly altered by the magnetic field while the total multiplicities of both positive and negative pions remain about the same. At beam energies above about 1 GeV/nucleon, while the integrated ratio of total $π^-$ to $π^+$ multiplicities is not, the differential $π^-/π^+$ ratio is sensitive to the density dependence of nuclear symmetry energy $E_{\rm{sym}}(ρ)$. Our findings suggest that magnetic effects should be carefully considered in future studies of using the differential $π^-/π^+$ ratio as a probe of the $E_{\rm{sym}}(ρ)$ at supra-saturation densities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Short-range tensor interaction and high-density nuclear symmetry energy

Effects of the short-range tensor interaction on the density-dependence of nuclear symmetry energy are examined by applying an approximate expression for the second-order tensor contribution to the symmetry energy derived earlier by G.E. Brown and R. Machleidt. It is found that the uncertainty in the short-range tensor force leads directly to a divergent high-density behavior of the nuclear symmetry energy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Tensor force induced isospin-dependence of short-range nucleon-nucleon correlation and high-density behavior of nuclear symmetry energy

Quantitative information on the tensor force induced isospin-dependence of short-range nucleon-nucleon correlation (SRC) extracted from recent J-Lab experiments was used in constraining the high-momentum tails of single nucleon momentum distributions in both symmetric nuclear matter (SNM) and pure neutron matter (PNM). Its effects on the Equations of State of SNM and PNM as well as the nuclear symmetry energy are investigated. It is found that the tensor force induced isospin-dependence of SRC softens significantly the nuclear symmetry energy especially at supra-saturation densities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The nuclear symmetry energy, the inner crust, and global neutron star modeling

The structure and composition of the inner crust of neutron stars, as well as global stellar properties such as radius and moment of inertia, have been shown to correlate with parameters characterizing the symmetry energy of nuclear matter such as its magnitude J and density dependence L at saturation density. It is thus mutually beneficial to nuclear physicists and astrophysicists to examine the combined effects of such correlations on potential neutron star observables in the light of recent experimental and theoretical constraints on J, L, and relationships between them. We review some basic correlations between these nuclear and astrophysical observables, and illustrate the impact of recent progress in constraining the J-L parameter space on the composition of the inner crust, crust-core transition density and pressure, and extent of the hypothesized pasta region. We use a simple compressible liquid drop model in conjunction with a simple model of nuclear matter which allows for independent, smooth, variation of the J and L. We extend the model into the core using the same nuclear matter model to explore the effects on global crust and core properties, and on potential observables such as crust oscillation frequencies and mechanically supported crust deformation. Throughout we illustrate the importance of the relationship between J and L implicit in a particular model of nuclear matter to the predictions of neutron star properties.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Upper limits on the observational effects of nuclear pasta in neutron stars

The effects of the existence of exotic nuclear shapes at the bottom of the neutron star inner crust - nuclear `pasta' - on observational phenomena are estimated by comparing the limiting cases that those phases have a vanishing shear modulus and that they have the shear modulus of a crystalline solid . We estimate the effect on torsional crustal vibrations and on the maximum quadrupole ellipticity sustainable by the crust. The crust composition and transition densities are calculated consistently with the global properties, using a liquid drop model with a bulk nuclear equation of state (EoS) which allows a systematic variation of the nuclear symmetry energy. The symmetry energy J and its density dependence L at nuclear saturation density are the dominant nuclear inputs which determine the thickness of the crust, the range of densities at which pasta might appear, as well as global properties such as the radius and moment of inertia. We show the importance of calculating the global neutron star properties on the same footing as the crust EoS, and demonstrate that in the range of experimentally acceptable values of L, the pasta phase can alter the crust frequencies by up to a factor of three, exceeding the effects of superfluidity on the crust modes, and decrease the maximum quadrupole ellipticity sustainable by the crust by up to an order of magnitude. The signature of the pasta phases and the density dependence of the symmetry energy on the potential observables highlights the possibility of constraining the EoS of dense, neutron-rich matter and the properties of the pasta phases using astrophysical observations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An improved single particle potential for transport model simulations of nuclear reactions induced by rare isotope beams

Taking into account more accurately the isospin dependence of nucleon-nucleon interactions in the in-medium many-body force term of the Gogny effective interaction, new expressions for the single nucleon potential and the symmetry energy are derived. Effects of both the spin(isospin) and the density dependence of nuclear effective interactions on the symmetry potential and the symmetry energy are examined. It is shown that they both play a crucial role in determining the symmetry potential and the symmetry energy at supra-saturation densities. The improved single nucleon potential will be useful for simulating more accurately nuclear reactions induced by rare isotope beams within transport models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Charged hadrons and nuclear parton distributions in p(d)A collisions

Nuclear gluon modifications are the least constrained component of current global fits to nuclear parton distributions, due to the inadequate constraining power of presently available experimental data from nuclear deep inelastic scattering and nuclear Drell-Yan lepton-pair production. A recent advance is the use of observables from relativistic nucleus-nucleus collisions to supplement the data pool for global fits. It is thus of interest to investigate the sensitivity of various experimental observables to different strengths of nuclear gluon modifications from large to small Bjorken $x$. In this work we utilize three recent global fits with different gluon strengths to investigate the sensitivity of three observables: nuclear modification factor, pseudorapidity asymmetry, and charge ratio. We observe that both nuclear modification factor and pseudorapidity asymmetry are quite sensitive to the strength of gluon modifications in a wide pseudorapidity interval. The sensitivity is greatly enhanced at LHC (Large Hadron Collider) energies relative to that at RHIC (Relativistic Heavy Ion Collider). The charge ratio is mildly sensitive only at large Bjorken x. Thus measurement of these observables in proton-lead collisions at the LHC affords the potential to further constrain gluon modifications in global fits.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraining the gravitational binding energy of PSR J0737-3039B using terrestrial nuclear data

We show that the gravitational binding energy of a neutron star of a given mass is correlated with the slope of the nuclear symmetry energy at 1-2 times nuclear saturation density for equations of state without significant softening (i.e., those that predict maximum masses $M_{\rm max} > 1.44M_{\odot}$ in line with the largest accurately measured neutron star mass). Applying recent laboratory constraints on the slope of the symmetry energy to this correlation we extract a constraint on the baryon mass of the lower mass member of the double pulsar binary system, PSR J0737-3039B. We compare with independent constraints derived from modeling the progenitor star of J0737-3039B up to and through its collapse under the assumption that it formed in an electron capture supernova. The two sets of constraints are consistent only if $L \lesssim$ 70 MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Energy release from hadron-quark phase transition in neutron stars and the axial $w$-mode of gravitational waves

Describing the hyperonic and quark phases of neutron stars with an isospin- and momentum-dependent effective interaction for the baryon octet and the MIT bag model, respectively, and using the Gibbs conditions to construct the mixed phase, we study the energy release due to the hadron-quark phase transition. Moreover, the frequency and damping time of the first axial $w$-mode of gravitational waves are studied for both hyperonic and hybrid stars. We find that the energy release is much more sensitive to the bag constant than the density dependence of the nuclear symmetry energy. Also, the frequency of the $w$-mode is found to be significantly different with or without the hadron-quark phase transition and depends strongly on the value of the bag constant. Effects of the density dependence of the nuclear symmetry energy become, however, important for large values of the bag constant that lead to higher hadron-quark transition densities.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Impact parameter profiles from nuclear shadowing ratio and applications to deuteron-gold collisions

While current nuclear parton distribution functions (nPDFs) from global fits to experimental data are spatially homogeneous, many experimental observables in nucleus-nucleus collisions are presented in terms of centrality cuts. These cuts can be related to impact parameter using the Glauber theory and it is thus usual in the description of such observables to convolute an assumed impact parameter distribution with the homogeneous nPDFs. In this study we use the Gribov theory of nuclear shadowing supplemented with information from diffraction to model the impact parameter distributions of nuclear shadowing ratio in the small-$x$ region. The modeled distributions are applied to the description of the centrality dependence of observables in deuteron-gold (d+Au) collisions at $\sqrt{s_{NN}} = 200$ AGeV.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nuclear limits on properties of pulsars and gravitational waves

Pulsars are among the most mysterious astrophysical objects in the Universe and are believed to be rotating neutron stars formed in supernova explosions. They are unique testing grounds of dense matter theories and gravitational physics and also provide links among nuclear physics, particle physics and General Relativity. Neutron stars may exhibit some of the most extreme and exotic characteristics that could not be found elsewhere in the Universe. Their properties are largely determined by the equation of state (EOS) of neutron-rich matter, which is the chief ingredient in calculating neutron star structure and properties of related phenomena, such as gravitational wave emission from deformed pulsars. Presently, the EOS of neutron-rich matter is still very uncertain mainly due to the poorly known density dependence of the nuclear symmetry energy especially at supra-saturation densities. Nevertheless, significant progress has been made recently in constraining the density dependence of the nuclear symmetry energy mostly at sub-saturation densities using terrestrial nuclear reactions. While there are still some uncertainties especially at supra-saturation densities, these constraints could provide useful information on the limits of the global properties of pulsars and the gravitational waves to be expected from them. Here we review our recent work on constraining properties of pulsars and gravitational radiation with data from terrestrial nuclear laboratories.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Why is the nuclear symmetry energy so uncertain at supra-saturation densities?

Within the interacting Fermi gas model for isospin asymmetric nuclear matter, effects of the in-medium three-body interaction and the two-body short-range tensor force due to the $ρ$ meson exchange as well as the short-range nucleon correlation on the high-density behavior of the nuclear symmetry energy are demonstrated respectively in a transparent way. Possible physics origins of the extremely uncertain nuclear symmetry energy at supra-saturation densities are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Effects of medium-induced $ρ-ω$ meson mixing on the equation of state in isospin-asymmetric nuclear matter

We reexamine effects of the $ρ-ω$ meson mixing mediated by nucleon polarizations on the symmetry energy in isospin-asymmetric nuclear matter. Taking into account the rearrangement term neglected in previous studies by others, we evaluate the $ρ-ω$ mixing angle in a novel way within the Relativistic Mean-Field Models with and without chiral limits. It is found that the symmetry energy is significantly softened at high densities contrary to the finding in earlier studies. As the first step of going beyond the lowest-order calculations, we also solve the RPA equation for the $ρ-ω$ mixing. In this case, it is found that the symmetry energy is not only significantly softened by the $ρ-ω$ mixing at supra-saturation densities, similar to the lowest-order $ρ-ω$ mixing, but interestingly also softened at subsaturation densities. In addition, the softening of the symmetry energy at subsaturation densities can be partly suppressed by the nonlinear self-interaction of the $σ$ meson.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Triton-3He relative and differential flows and the high density behavior of nuclear symmetry energy

Using a transport model coupled with a phase-space coalescence after-burner we study the triton-3He relative and differential transverse flows in semi-central 132Sn+124Sn reactions at a beam energy of 400 MeV/nucleon. We find that the triton-3He pairs carry interesting information about the density dependence of the nuclear symmetry energy. The t-3He relative flow can be used as a particularly powerful probe of the high-density behavior of the nuclear symmetry energy.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Nuclear limits on gravitational waves from elliptically deformed pulsars

Gravitational radiation is a fundamental prediction of General Relativity. Elliptically deformed pulsars are among the possible sources emitting gravitational waves (GWs) with a strain-amplitude dependent upon the star's quadrupole moment, rotational frequency, and distance from the detector. We show that the gravitational wave strain amplitude $h_0$ depends strongly on the equation of state of neutron-rich stellar matter. Applying an equation of state with symmetry energy constrained by recent nuclear laboratory data, we set an upper limit on the strain-amplitude of GWs produced by elliptically deformed pulsars. Depending on details of the EOS, for several millisecond pulsars at distances $0.18kpc$ to $0.35kpc$ from Earth, the {\it maximal} $h_0$ is found to be in the range of $\sim[0.4-1.5]\times 10^{-24}$. This prediction serves as the first {\it direct} nuclear constraint on the gravitational radiation. Its implications are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Recent Progress and New Challenges in Isospin Physics with Heavy-Ion Reactions

The ultimate goal of studying isospin physics via heavy-ion reactions with neutron-rich, stable and/or radioactive nuclei is to explore the isospin dependence of in-medium nuclear effective interactions and the equation of state of neutron-rich nuclear matter, particularly the isospin-dependent term in the equation of state, i.e., the density dependence of the symmetry energy. Because of its great importance for understanding many phenomena in both nuclear physics and astrophysics, the study of the density dependence of the nuclear symmetry energy has been the main focus of the intermediate-energy heavy-ion physics community during the last decade, and significant progress has been achieved both experimentally and theoretically. In particular, a number of phenomena or observables have been identified as sensitive probes to the density dependence of the nuclear symmetry energy. Experimental studies have confirmed some of these interesting isospin-dependent effects and allowed us to constrain relatively stringently the symmetry energy at sub-saturation densities. The impacts of this constrained density dependence of the symmetry energy on the properties of neutron stars have also been studied, and they were found to be very useful for the astrophysical community. With new opportunities provided by the various radioactive beam facilities being constructed around the world, the study of isospin physics is expected to remain one of the forefront research areas in nuclear physics. In this report, we review the major progress achieved during the last decade in isospin physics with heavy ion reactions and discuss future challenges to the most important issues in this field.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

A Transport Model for Nuclear Reactions Induced by Radioactive Beams

Major ingredients of an isospin and momentum dependent transport model for nuclear reactions induced by radioactive beams are outlined. Within the IBUU04 version of this model we study several experimental probes of the equation of state of neutron-rich matter, especially the density dependence of the nuclear symmetry energy. Comparing with the recent experimental data from NSCL/MSU on isospin diffusion, we found a nuclear symmetry energy of $% E_{sym}(ρ)\approx 31.6(ρ/ρ_{0})^{1.05}$ at subnormal densities. Predictions on several observables sensitive to the density dependence of the symmetry energy at supranormal densities accessible at GSI and the planned Rare Isotope Accelerator (RIA) are also made.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Constraining the neutron-proton effective mass splitting in neutron-rich matter

Within Bombaci's phenomenological single-nucleon potential model we study the neutron-proton effective mass splitting $m_n^*- m_p^*$ in neutron-rich matter. It is shown that an effective mass splitting of $m_n^*<m_p^*$ leads to a symmetry potential that is inconsistent with the energy dependence of the Lane potential constrained by the nucleon-nucleus scattering experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Light cluster production in intermediate energy heavy-ion collisions induced by neutron-rich nuclei

The coalescence model based on nucleon distribution functions from an isospin-dependent transport model is used to study the production of light clusters such as deuteron, triton, and $^{3}$He from heavy-ion collisions induced by neutron-rich nuclei at intermediate energies. It is found that the emission time of light clusters depends on their masses. For clusters with the same momentum per nucleon, heavier ones are emitted earlier. Both the yield and energy spectrum of light clusters are sensitive to the density dependence of nuclear symmetry energy, with more light clusters produced in the case of a stiff symmetry energy. On the other hand, effects due to the stiffness of the isoscalar part of nuclear equation of state and the medium dependence of nucleon-nucleon cross sections on light cluster production are unimportant. We have also studied the correlation functions of clusters, and they are affected by the density dependence of nuclear symmetry energy as well, with the stiff symmetry energy giving a stronger anti-correlation of light clusters, particularly for those with large kinetic energies. Dependence of light cluster production on the centrality and incident energy of heavy ion collisions as well as the mass of the reaction system is also investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Probing the EOS of dense neutron-rich matter with high-energy radioactive beams

Nuclear reactions induced by high energy radioactive beams create a transient state of nuclear matter with high density and appreciable neutron to proton asymmetry. This will provide a unique opportunity to explore novel properties of dense neutron-rich matter and the isospin-dependence of the nuclear equation of state (EOS). Here we study the $π^-/π^+$ ratio as a probe of the EOS of dense neutron-rich matter.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

High Density Behaviour of Nuclear Symmetry Energy and High Energy Heavy-Ion Collisions

High energy heavy-ion collisions are proposed as a novel means to obtain information about the high density ({\rm HD}) behaviour of nuclear symmetry energy. Within an isospin-dependent hadronic transport model using phenomenological equations of state ({\rm EOS}) for dense neutron-rich matter, it is shown that the isospin asymmetry of the HD nuclear matter formed in high energy heavy-ion collisions is determined mainly by the HD behaviour of nuclear symmetry energy. Experimental signatures in several sensitive probes, i.e., $π^-$ to $π^+$ ratio, transverse collective flow and its excitation function as well as neutron-proton differential flow, are investigated. A precursor of the possible isospin separation instability in dense neutron-rich matter is predicted to appear as the local minima in the excitation functions of the transverse flow parameter for both neutrons and protons above the pion production threshold. Because of its {\it qualitative} nature unlike other {\it quantitative} observables, this precursor can be used as a unique signature of the isospin dependence of the nuclear {\rm EOS}. Measurements of these observables will provide the first terrestrial data to constrain stringently the HD behaviour of nuclear symmetry energy and thus also the {\rm EOS} of dense neutron-rich matter. Implications of our findings to neutron star studies are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Probing the High Density Behaviour of Nuclear Symmetry Energy with High Energy Heavy-Ion Collisions

High energy heavy-ion collisions are proposed as a novel means to constrain stringently the high density (HD) behaviour of nuclear symmetry energy. Within an isospin-dependent hadronic transport model, it is shown for the first time that the isospin asymmetry of the HD nuclear matter formed in high energy heavy-ion collisions is uniquely determined by the HD behaviour of nuclear symmetry energy. Experimental signatures in two sensitive probes, i.e., $π^-$ to $π^+$ ratio and neutron-proton differential collective flow, are also investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Charged particle rapidity distributions at RHIC

Using a multiphase transport model (AMPT), which includes both initial partonic and final hadronic interactions, we study the rapidity distributions of charged particles such as protons, antiprotons, pions, and kaons in heavy ion collisions at RHIC. The theoretical results for the total charged particle multiplicity at midrapidity are consistent with those measured by the PHOBOS collaboration in central Au+Au collisions at $\sqrt s$ = 56 and 130 AGeV. We find that these hadronic observables are much more sensitive to the hadronic interactions than the partonic interactions.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Equation of state of asymmetric nuclear matter and collisions of neutron-rich nuclei

The ratio of pre-equilibrium neutrons to protons from collisions of neutron-rich nuclei is studied as a function of their kinetic energies. This ratio is found to be sensitive to the density dependence of the nuclear symmetry energy, but is independent of the compressibility of symmetric nuclear matter and the in-medium nucleon-nucleon cross sections. The experimental measurement of this ratio thus provides a novel means for determining the nuclear equation of state of asymmetric nuclear matter.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Isospin Physics in Heavy-Ion Collisions at Intermediate Energies

In nuclear collisions induced by stable or radioactive neutron-rich nuclei a transient state of nuclear matter with an appreciable isospin asymmetry as well as thermal and compressional excitation can be created. This offers the possibility to study the properties of nuclear matter in the region between symmetric nuclear matter and pure neutron matter. In this review, we discuss recent theoretical studies of the equation of state of isospin-asymmetric nuclear matter and its relations to the properties of neutron stars and radioactive nuclei. Chemical and mechanical instabilities as well as the liquid-gas phase transition in asymmetric nuclear matter are investigated. The in-medium nucleon-nucleon cross sections at different isospin states are reviewed as they affect significantly the dynamics of heavy ion collisions induced by radioactive beams. We then discuss an isospin-dependent transport model, which includes different mean-field potentials and cross sections for the proton and neutron, and its application to these reactions. Furthermore, we review the comparisons between theoretical predictions and available experimental data. In particular, we discuss the study of nuclear stopping in terms of isospin equilibration, the dependence of nuclear collective flow and balance energy on the isospin-dependent nuclear equation of state and cross sections, the isospin dependence of total nuclear reaction cross sections, and the role of isospin in preequilibrium nucleon emissions and subthreshold pion production.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Excitation functions in central Au+Au collisions from SIS/GSI to AGS/Brookhaven

Using the relativistic transport model (ART), we predict the energy dependence of the stopping power, maximum baryon and energy densities, the population of resonance matter as well as the strength of the transverse and radial flow for central Au+Au reactions at beam momentum from 2 to 12 GeV/c available at Brookhaven's AGS. The maximum baryon and energy densities are further compared to the predictions of relativistic hydrodynamics assuming the formation of shock waves. We also discuss the Fermi-Landau scaling of the pion multiplicity in these reactions.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Pion Transparency in 500 MeV C(pi,pi') Reactions?

The question whether there is a pion transparency in 500 MeV $C(π, π^{'})$ scatterings is studied using a semiclassical, hadronic transport model. The double differential cross sections of this reaction measured at LAMPF can be largely accounted for, if one uses energy-dependent, anisotropic angular distributions which are fitted to pion-nucleon scattering data for the decay of $Δ(1232)$ and $N^*(1440)$ resonances. The remaining discrepancy between the data and the calculation sets a limit on effects of more exotic processes.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Mean field effects in heavy-ion collisions at AGS energies

The question whether mean field effects exist in heavy-ion collisions at AGS energies is studied in the framework of A Relativistic Transport (ART) model. It is found that in central collisions of Au+Au at $P_{beam}/A=$11.6 GeV/c a simple, Skyrme-type nuclear mean field satisfying the causality requirement reduces the maximum baryon and energy densities reached in the cascade model by about 30\% and 40\%, respectively. The mean field affects the inclusive, single particle observables of various hadrons by at most 20\%. We show, however, that the mean field causes a factor of 2.5 increase in the strength of the baryon transverse collective flow.

0 citations0 reviewsNo code linkedOpen access