BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Bao-An Li

Bao-An Li contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
nucl-exnucl-thastro-ph.HEastro-ph.SRhep-phgr-qchep-th

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 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