Explore the vast range of titles available.

The density dependence of nuclear symmetry energy is among the most uncertain parts of the Equation of State (EOS) of dense neutron-rich nuclear matter. It is currently poorly known especially at suprasaturation densities partially because of our poor knowledge about isovector nuclear interactions at short distances. Because of its broad impacts on many interesting issues, pinning down the density dependence of nuclear symmetry energy has been a longstanding and shared goal of both astrophysics and nuclear physics. New observational data of neutron stars including their masses, radii, and tidal deformations since GW170817 have helped improve our knowledge about nuclear symmetry energy, especially at high densities. Based on various model analyses of these new data by many people in the nuclear astrophysics community, while our brief review might be incomplete and biased unintentionally, we learned in particular the following: (1) The slope parameter L of nuclear symmetry energy at saturation density ρ0 of nuclear matter from 24 new analyses of neutron star observables was about L≈57.7±19 MeV at a 68% confidence level, consistent with its fiducial value from surveys of over 50 earlier analyses of both terrestrial and astrophysical data within error bars. (2) The curvature Ksym of nuclear symmetry energy at ρ0 from 16 new analyses of neutron star observables was about Ksym≈−107±88 MeV at a 68% confidence level, in very good agreement with the systematics of earlier analyses. (3) The magnitude of nuclear symmetry energy at 2ρ0, i.e., Esym(2ρ0)≈51±13 MeV at a 68% confidence level, was extracted from nine new analyses of neutron star observables, consistent with the results from earlier analyses of heavy-ion reactions and the latest predictions of the state-of-the-art nuclear many-body theories. (4) While the available data from canonical neutron stars did not provide tight constraints on nuclear symmetry energy at densities above about 2ρ0, the lower radius boundary R2.01=12.2 km from NICER’s very recent observation of PSR J0740+6620 of mass 2.08±0.07M⊙ and radius R=12.2–16.3 km at a 68% confidence level set a tight lower limit for nuclear symmetry energy at densities above 2ρ0. (5) Bayesian inferences of nuclear symmetry energy using models encapsulating a first-order hadron–quark phase transition from observables of canonical neutron stars indicated that the phase transition shifted appreciably both L and Ksym to higher values, but with larger uncertainties compared to analyses assuming no such phase transition. (6) The high-density behavior of nuclear symmetry energy significantly affected the minimum frequency necessary to rotationally support GW190814’s secondary component of mass (2.50–2.67) M⊙ as the fastest and most massive pulsar discovered so far. Overall, thanks to the hard work of many people in the astrophysics and nuclear physics community, new data of neutron star observations since the discovery of GW170817 have significantly enriched our knowledge about the symmetry energy of dense neutron-rich nuclear matter.

We investigate the possibility that GW170817 was not the merger of two conventional neutron stars (NS), but involved at least one if not two hybrid stars with a quark matter core that might even belong to a third family of compact stars. To this end, we develop a Bayesian analysis method for selecting the most probable equation of state (EoS) under a set of constraints from compact star physics, which now also include the tidal deformability from GW170817 and the first result for the mass and radius determination for PSR J0030+0451 by the NICER Collaboration. We apply this method for the first time to a two-parameter family of hybrid EoS based on the DD2 model with nucleonic excluded volume for hadronic matter and the color superconducting generalized nlNJL model for quark matter. The model has a variable onset density for deconfinement and can mimic the effects of pasta phases with the possibility of producing a third family of hybrid stars in the mass-radius diagram. The main findings of this study are that: (1) the presence of multiple configurations for a given mass (twins or even triples) corresponds to a set of disconnected lines in the Λ 1 – Λ 2 diagram of tidal deformabilities for binary mergers, so that merger events from the same mass range may result in a probability landscape with different peak positions; (2) the Bayesian analysis with the above observational constraints favors an early onset of the deconfinement transition, at masses of M onset ≤ 0.8 M ⊙ with an M–R relationship that in the range of observed neutron star masses is almost indistinguishable from that of a soft hadronic Akmal, Pandharipande, and Ravenhall (APR) EoS; (3) a few, yet fictitious measurements of the NICER experiment two times more accurate than the present value and a different mass and radius that would change the posterior likelihood so that hybrid EoS with a phase transition onset in the range M onset = 1.1–1.6 M ⊙ would be favored.

Strange stars ought to exist in the universe according to the strange quark matter hypothesis, which states that matter made of roughly equal numbers of up, down, and strange quarks could be the true ground state of baryonic matter rather than ordinary atomic nuclei. Theoretical models of strange quark matter, such as the standard MIT bag model, the density-dependent quark mass model, or the quasi-particle model, however, appear to be unable to reproduce some of the properties (masses, radii, and tidal deformabilities) of recently observed compact stars. This is different if alternative gravity theory (e.g., non-Newtonian gravity) or dark matter (e.g., mirror dark matter) are considered, which resolve these issues. The possible existence of strange stars could thus provide a clue to new physics, as discussed in this review.

The presence of dark matter (DM) within neutron stars (NSs) can be introduced by different accumulation scenarios in which DM and baryonic matter (BM) may interact only through the gravitational force. In this work, we consider asymmetric self-interacting bosonic DM, which can reside as a dense core inside the NS or form an extended halo around it. It is seen that depending on the boson mass (mχ), self-coupling constant (λ) and DM fraction (Fχ), the maximum mass, radius and tidal deformability of NSs with DM admixture will be altered significantly. The impact of DM causes some modifications in the observable features induced solely by the BM component. Here, we focus on the widely used nuclear matter equation of state (EoS) called DD2 for describing NS matter. We show that by involving DM in NSs, the corresponding observational parameters will be changed to be consistent with the latest multi-messenger observations of NSs. It is seen that for mχ≳200 MeV and λ≲2π, DM-admixed NSs with 4%≲Fχ≲20% are consistent with the maximum mass and tidal deformability constraints.

On 14 August 2019, the LIGO/Virgo collaboration observed a compact object with mass ∼2.59−0.09+0.08M⊙, as a component of a system where the main companion was a black hole with mass ∼23M⊙. A scientific debate initiated concerning the identification of the low mass component, as it falls into the neutron star–black hole mass gap. The understanding of the nature of GW190814 event will offer rich information concerning open issues, the speed of sound and the possible phase transition into other degrees of freedom. In the present work, we made an effort to probe the nuclear equation of state along with the GW190814 event. Firstly, we examine possible constraints on the nuclear equation of state inferred from the consideration that the low mass companion is a slow or rapidly rotating neutron star. In this case, the role of the upper bounds on the speed of sound is revealed, in connection with the dense nuclear matter properties. Secondly, we systematically study the tidal deformability of a possible high mass candidate existing as an individual star or as a component one in a binary neutron star system. As the tidal deformability and radius are quantities very sensitive on the neutron star equation of state, they are excellent counters on dense matter properties. We conjecture that similar isolated neutron stars or systems may exist in the universe and their possible future observation will shed light on the maximum neutron star mass problem.

Neutron stars (NSs), among the densest objects in the universe, are exceptional laboratories for investigating the properties of dark matter (DM). Recent theoretical and observational developments have heightened interest in exploring the impact of DM on NS structure, giving rise to the concept of dark matter admixed neutron stars (DANSs). This review examines how NSs can accumulate DM over time, potentially altering their fundamental properties. We explore the leading models describing DM behavior within NSs, focusing on the effects of both bosonic and fermionic candidates on key features such as mass, radius, and tidal deformability. Additionally, we review how DM can modify the cooling and heating processes, trigger the formation of a black hole, and impact gravitational wave (GW) emissions from binary systems. By synthesizing recent research, this work highlights how DANSs might produce observable signatures, offering new opportunities to probe DM’s properties through astrophysical phenomena.

The effective Gogny interactions of the D1 family were established by D. Gogny more than forty years ago with the aim to describe simultaneously the mean field and the pairing field corresponding to the nuclear interaction. The most popular Gogny parametrizations, namely D1S, D1N and D1M, describe accurately the ground-state properties of spherical and deformed finite nuclei all across the mass table obtained with Hartree–Fock–Bogoliubov (HFB) calculations. However, these forces produce a rather soft equation of state (EoS) in neutron matter, which leads to predict maximum masses of neutron stars well below the observed value of two solar masses. To remove this limitation, we built new Gogny parametrizations by modifying the density dependence of the symmetry energy predicted by the force in such a way that they can be applied to the neutron star domain and can also reproduce the properties of finite nuclei as good as their predecessors. These new parametrizations allow us to obtain stiffer EoS’s based on the Gogny interactions, which predict maximum masses of neutron stars around two solar masses. Moreover, other global properties of the star, such as the moment of inertia and the tidal deformability, are in harmony with those obtained with other well tested EoSs based on the SLy4 Skyrme force or the Barcelona–Catania–Paris–Madrid (BCPM) energy density functional. Properties of the core-crust transition predicted by these Gogny EoSs are also analyzed. Using these new Gogny forces, the EoS in the inner crust is obtained with the Wigner–Seitz approximation in the Variational Wigner–Kirkwood approach along with the Strutinsky integral method, which allows one to estimate in a perturbative way the proton shell and pairing corrections. For the outer crust, the EoS is determined basically by the nuclear masses, which are taken from the experiments, wherever they are available, or by HFB calculations performed with these new forces if the experimental masses are not known.

We develop a Bayesian analysis method for selecting the most probable equation of state under a set of constraints from compact star physics, which now include the tidal deformability from GW170817. We apply this method for the first time to a two-parameter family of hybrid equations of state that is based on realistic models for the hadronic phase (KVORcut02) and the quark matter phase (SFM α ) which produce a third family of hybrid stars in the mass–radius diagram. One parameter ( α ) characterizes the screening of the string tension in the string-flip model of quark matter while the other ( Δ P ) belongs to the mixed phase construction that mimics the thermodynamics of pasta phases and includes the Maxwell construction as a limiting case for Δ P = 0 . We present the corresponding results for compact star properties like mass, radius and tidal deformabilities and use empirical data for them in the newly developed Bayesian analysis method to obtain the probabilities for the model parameters within their considered range.

Isospin symmetry breaking effects on the mass-radius relation of a cold, non-accreting neutron star are studied on the basis of two Skyrme Energy Density Functionals (EDFs). One functional contains isospin symmetry breaking terms other than those typically included in Skyrme EDFs while its counterpart is of standard form. Both functionals are based on the same fitting protocol except for the observables and pseudo-observables sensitive to the isospin symmetry breaking channel. The quality of those functionals is similar in the description of terrestrial observables but choosing either of them has a non-negligible effect on the mass-radius relation and tidal deformability of a neutron star. Further investigations are needed to clarify the effects of isospin symmetry breaking on these and other observables of neutron stars that are, and will become, available.

We construct a new effective field theory approach to the equation of state (EoS), dubbed pseudo-confomal model “PCM”, for nuclear and compact star matter entirely in terms of effective hadron degrees of freedom. The possible transition at n ∼(2–4) n 0 (where n 0 is the normal nuclear matter density) from hadron degrees of freedom to strongly-coupled quark degrees of freedom, giving rise to a soft-to-hard changeover in the EoS that can accommodate the massive stars observed, is effectuated by the topology change at n 1/2 ≳ 2 n 0 from skyrmions to half-skyrmions without involving local order-parameter fields. The mechanism exploits possible emergence of hidden scale and local symmetries of QCD at high density, leading to a precocious “pseudo-conformal” sound velocity v s 2 = 1/3 (in unit of c = 1) for n ≳ 3 n 0 . The resulting prediction signals a drastic departure from standard nuclear many-body theory in the density regime involved in the massive stars. We suggest that the tidal deformability implemented in gravitational waves coming from coalescing neutron stars in LIGO/Virgo-type observations could pin down the location of the topology change density n 1/2 .