Observations Of Neutron Stars Research Articles

Towards mitigation of apparent tension between nuclear physics and astrophysical observations by improved modeling of neutron star matter

Observations of neutron stars (NSs) by the LIGO-Virgo and NICER collaborations have provided reasonably precise measurements of their various macroscopic properties. In this paper, we employ a Bayesian framework to combine them and place improved joint constraints on the properties of the NS equation of state (EOS). We use a hybrid EOS formulation that employs a parabolic expansion-based nuclear empirical parametrization around the nuclear saturation density augmented by a generic 3-segment piecewise polytrope model at higher densities. Within the 90% credible level this parametrization predicts ${R}_{1.4}={12.57}_{\ensuremath{-}0.92}^{+0.73}\text{ }\text{ }\mathrm{km}$ and ${\mathrm{\ensuremath{\Lambda}}}_{1.4}={550}_{\ensuremath{-}225}^{+223}$ for the radius and dimensionless tidal deformability, respectively, of a $1.4\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ NS. Finally, we show how the construction of the full NS EOS based solely on the nuclear empirical parameters at saturation density leads to certain tension with the astrophysical data, and how the hybrid approach provides a resolution to it.

Extracting Hypernuclear Properties from the (e, e′K+) Cross Section

Experimental studies of hypernuclear dynamics, besides being essential for the understanding of strong interactions in the strange sector, have important astrophysical implications. The observation of neutron stars with masses exceeding two solar masses poses a serious challenge to the models of hyperon dynamics in dense nuclear matter, many of which predict a maximum mass incompatible with the data. In this paper, it is argued that valuable new insight can be gained from the forthcoming extension of the experimental studies of kaon electro production from nuclei to include the 208Pb(e,e′K+)Λ208Tl process. A comprehensive framework for the description of kaon electro production, based on factorization of the nuclear cross section and the formalism of the nuclear many-body theory, is outlined. This approach highlights the connection between the kaon production and proton knockout reactions, which will allow us to exploit the available 208Pb(e,e′p)207Tl data to achieve a largely model-independent analysis of the measured cross section.

Astrophysical and Theoretical Physics Implications from Multimessenger Neutron Star Observations.

The Neutron Star Interior Composition Explorer (NICER) recently measured the mass and equatorial radius of the isolated neutron star PSR J0030+0451. We use these measurements to infer the moment of inertia, the quadrupole moment, and the surface eccentricity of an isolated neutron star for the first time, using relations between these quantities that are insensitive to the unknown equation of state of supranuclear matter. We also use these results to forecast the moment of inertia of neutron star A in the double pulsar binary J0737-3039, a quantity anticipated to be directly measured in the coming decade with radio observations. Combining this information with the measurement of the tidal Love number with LIGO/Virgo observations, we propose and implement the first theory-agnostic and equation-of-state-insensitive test of general relativity. Specializing these constraints to a particular modified theory, we find that consistency with general relativity places the most stringent constraint on gravitational parity violation to date, surpassing all other previously reported bounds by 7 orders of magnitude and opens the path for a future test of general relativity with multimessenger neutron star observations.

Constraining nuclear matter parameters from correlation systematics: a mean-field perspective

The nuclear matter parameters define the nuclear equation of state (EoS), they appear as coefficients of expansion around the saturation density of symmetric and asymmetric nuclear matter. We review their correlations with several properties of finite nuclei and of neutron stars within mean-field frameworks. The lower order nuclear matter parameters such as the binding energy per nucleon, incompressibility and the symmetry energy coefficients are found to be constrained in narrow limits through their strong ties with selective properties of finite nuclei. From the correlations of nuclear matter parameters with neutron star observables, we further review how precision knowledge of the radii and tidal deformability of neutron stars in the mass range $1 - 2 M_\odot$ may help cast them in narrower bounds. The higher order parameters such as the density slope and the curvature of the symmetry energy or the skewness of the symmetric nuclear matter EoS are, however, plagued with larger uncertainty. From inter-correlation of these higher order nuclear matter parameters with lower order ones, we explore how they can be brought to more harmonious bounds.

Limiting masses and radii of neutron stars and their implications

We combine equation of state of dense matter up to twice nuclear saturation density ($n_{\rm sat}=0.16\, \text{fm}^{-3}$) obtained using chiral effective field theory ($\chi$EFT), and recent observations of neutron stars to gain insights about the high-density matter encountered in their cores. A key element in our study is the recent Bayesian analysis of correlated EFT truncation errors based on order-by-order calculations up to next-to-next-to-next-to-leading order in the $\chi$EFT expansion. We refine the bounds on the maximum mass imposed by causality at high densities, and provide stringent limits on the maximum and minimum radii of $\sim1.4\,{\rm M}_{\odot}$ and $\sim2.0\,{\rm M}_{\odot}$ stars. Including $\chi$EFT predictions from $n_{\rm sat}$ to $2\,n_{\rm sat}$ reduces the permitted ranges of the radius of a $1.4\,{\rm M}_{\odot}$ star, $R_{1.4}$, by $\sim3.5\, \text{km}$. If observations indicate $R_{1.4}<11.2\, \text{km}$, our study implies that either the squared speed of sound $c^2_{s}>1/2$ for densities above $2\,n_{\rm sat}$, or that $\chi$EFT breaks down below $2\,n_{\rm sat}$. We also comment on the nature of the secondary compact object in GW190814 with mass $\simeq 2.6\,{\rm M}_{\odot}$, and discuss the implications of massive neutron stars $>2.1 \,{\rm M}_{\odot}\,(2.6\,{\rm M}_{\odot})$ in future radio and gravitational-wave searches. Some form of strongly interacting matter with $c^2_{s}>0.35\, (0.55)$ must be realized in the cores of such massive neutron stars. In the absence of phase transitions below $2\,n_{\rm sat}$, the small tidal deformability inferred from GW170817 lends support for the relatively small pressure predicted by $\chi$EFT for the baryon density $n_{\rm B}$ in the range $1-2\,n_{\rm sat}$. Together they imply that the rapid stiffening required to support a high maximum mass should occur only when $n_{\rm B} \gtrsim 1.5-1.8\,n_{\rm sat}$.

Constraints on the Time Variation of the Gravitational Constant Using Gravitational Wave Observations of Binary Neutron Stars.

We propose a method to constrain the variation of the gravitational constant G with cosmic time using gravitational wave (GW) observations of merging binary neutron stars. The method essentially relies on the fact that the maximum and minimum allowed masses of neutron stars at a particular cosmic epoch have a simple dependence on the value of G at that epoch.GWs carry an imprint of the value of G at the time of the merger. Thus, if the value of G at merger is significantly different from its current value, the masses of the neutron stars inferred from the GW observations will be inconsistent with the theoretically allowed range. This enables us to place bounds on the variation of G between the merger epoch and the present epoch. Using the observation of the binary neutron star system GW170817, we constrain the fractional difference in G between the merger and the current epoch to be in the range -1≲ΔG/G≲8. Assuming a monotonic variation in G, this corresponds to a bound on the average rate of change of -7×10^{-9} yr^{-1}≤G[over ˙]/G≤5×10^{-8} yr^{-1} between these epochs. Future observations will put tight constraints on the deviation of G over vast cosmological epochs not probed by other observations.

Constraints on Nuclear Saturation Properties from Terrestrial Experiments and Astrophysical Observations of Neutron Stars

Abstract Taking into account the terrestrial experiments and the recent astrophysical observations of neutron stars and gravitational-wave signals, we impose restrictions on the equation of state (EoS) for isospin-asymmetric nuclear matter. Using the relativistic mean-field model with SU(3) flavor symmetry, we investigate the impacts of effective nucleon mass, nuclear incompressibility, and slope parameter of nuclear symmetry energy on the nuclear and neutron star properties. It is found that the astrophysical information of massive neutron stars and tidal deformabilities, as well as the nuclear experimental data, plays an important role to restrict the EoS for neutron stars. In particular, the softness of the nuclear EoS due to the existence of hyperons in the core gives stringent constraints on those physical quantities. Furthermore, it is possible to put limits on the curvature parameter of nuclear symmetry energy by means of nuclear and astrophysical calculations.

A Modern View of the Equation of State in Nuclear and Neutron Star Matter

Background: We analyze several constraints on the nuclear equation of state (EOS) currently available from neutron star (NS) observations and laboratory experiments and study the existence of possible correlations among properties of nuclear matter at saturation density with NS observables. Methods: We use a set of different models that include several phenomenological EOSs based on Skyrme and relativistic mean field models as well as microscopic calculations based on different many-body approaches, i.e., the (Dirac–)Brueckner–Hartree–Fock theories, Quantum Monte Carlo techniques, and the variational method. Results: We find that almost all the models considered are compatible with the laboratory constraints of the nuclear matter properties as well as with the largest NS mass observed up to now, 2.14−0.09+0.10M⊙ for the object PSR J0740+6620, and with the upper limit of the maximum mass of about 2.3–2.5M⊙ deduced from the analysis of the GW170817 NS merger event. Conclusion: Our study shows that whereas no correlation exists between the tidal deformability and the value of the nuclear symmetry energy at saturation for any value of the NS mass, very weak correlations seem to exist with the derivative of the nuclear symmetry energy and with the nuclear incompressibility.

Modelling of dark energy stars with tolman IV gravitational potential

The present study explores the properties of dark energy star with high density under the dark energy equation of state (EOS) (pr=ωρ) with anisotropic fluid distribution, where ω describes the dark energy parameter. The relativistic dark-energy stellar structures are studied by adopting Tolman IV-type gravitational potential. In this relation we examine many physical properties through graphical plots and mathematical calculation, such as energy condition, compactness, mass-radius ratio, anisotropic nature and surface redshift. It is observed that the strong energy condition is violated whereas all other energy conditions are satisfied. We have also fitted radius corresponding to observed mass value for observed neutron stars, namely, PSRJ1416−2230,VelaX−1,4U1608−52,HerX−1 and PSRJ1903+327, which are compatible with the predicted values.

Impact of dynamical dark energy on the neutron star equilibrium

We study the density distribution of the minimally-coupled scalar field dark energy inside a neutron star. The dark energy is considered in the hydrodynamical representation as a perfect fluid with three parameters (background density, equation of state, and effective sound speed). The neutron star matter is modeled with three unified equations of state, developed by the Brussels-Montreal group. With the calculated density distribution of the dark energy inside a neutron star (and its dependence on the dark energy parameters) we investigate how its presence impacts the macroscopic characteristics and the value of the mass limit for neutron stars. From this impact we derive the possible constraints on the effective speed of sound of dark energy with the help of maximal known masses of observed neutron stars. In this approach, we have found, that the squared effective speed of sound can not be smaller than ∼ 10−2 in units of squared speed of light.

Resonance suppression of the r -mode instability in superfluid neutron stars: Accounting for muons and entrainment

We calculate the finite-temperature r-mode spectrum of a superfluid neutron star accounting for both muons in the core and the entrainment between neutrons and protons. We show that the standard perturbation scheme, considering the rotation rate as an expansion parameter, breaks down in this case. We develop an original perturbation scheme which circumvents this problem by treating both the perturbations due to rotation and (weak) entrainment simultaneously. Applying this scheme, we propose a simple method for calculating the superfluid r-mode eigenfrequency in the limit of vanishing rotation rate. We also calculate the r-mode spectrum at finite rotation rate for realistic microphysics input (adopting, however, the Newtonian framework and Cowling approximation when considering perturbed oscillation equations) and show that the normal r-mode exhibits resonances with superfluid r-modes at certain values of temperatures and rotation frequencies in the parameter range relevant to neutron stars in low-mass X-ray binaries (LMXBs). This turns the recently suggested phenomenological model of resonance r-mode stabilization into a quantitative theory, capable of explaining observations. A strong dependence of resonance rotation rates and temperatures on the neutron superfluidity model allows us to constrain the latter by confronting our calculations with the observations of neutron stars in LMXBs.

Thermal properties of hot and dense matter: Influence of rapid rotation on protoneutron stars, hot neutron stars, and neutron star merger remnants

The knowledge of the equation of state is a key ingredient for many dynamical phenomena that depend sensitively on the hot and dense nuclear matter, such as the formation of protoneutron stars and hot neutron stars. In order to accurately describe them, we construct equations of state at FInite temperature and entropy per baryon for matter with varying proton fractions. This procedure is based on the momentum dependent interaction model and state-of-the-art microscopic data. In addition, we investigate the role of thermal and rotation effects on microscopic and macroscopic properties of neutron stars, including the mass and radius, the frequency, the Kerr parameter, the central baryon density, etc. The latter is also connected to the hot and rapidly rotating remnant after neutron star merger. The interplay between these quantities and data from late observations of neutron stars, both isolated and in matter of merging, could provide useful insight and robust constraints on the equation of state of nuclear matter.

Nuclear Symmetry Energies at Supra-saturation Densities Extracted from Observations of Neutron Stars and Gravitational Waves

The high-density behavior of nuclear symmetry energy is the most uncertain part of the Equation of State (EOS) of dense neutron-rich nucleonic matter. It has significant ramifications in understanding properties of nuclear reactions induced by rare isotopes, neutron stars and gravitational waves from various sources. To infer the underlying EOS of dense neutron-rich nuclear matter from observables of neutron stars is the longstanding neutron star inverse-structure problem. To make progress in solving this problem, we advanced two inversion approaches: (1) the direct inversion in the three-dimensional high-density EOS parameter space and (2) Bayesian inference of six EOS parameters statistically. Neutron star observational data obtained so far, especially their radii from analyzing gravitational waves (e.g., via the tidal deformability) from GW170817 and X-rays from quiescent low-mass X-ray binaries can already help constrain significantly the nuclear EOS up to about twice the saturation density of nuclear matter. In particular, the nuclear symmetry energy at twice the saturation density of nuclear matter E sym(2ρ 0) is found to be MeV at 68% confidence level.

Speed of sound and quark confinement inside neutron stars

Several observations of high-mass neutron stars (NSs), as well as the first historic detection of the binary neutron star merger GW170817, have delivered stringent constraints on the equation of state (EoS) of cold and dense matter. Recent studies suggest that, in order to simultaneously accommodate a 2M⊙ NS and the upper limit on the compactness, the pressure has to swiftly increase with density and the corresponding speed of sound likely exceeds the conformal limit. In this work, we employ a unified description of hadron-quark matter, the hybrid quark-meson-nucleon (QMN) model, to investigate the EoS under NS conditions. We show that the dynamical confining mechanism of the model plays an important role in explaining the observed properties of NSs.

Direct astrophysical tests of chiral effective field theory at supranuclear densities

Recent observations of neutron stars with gravitational waves and X-ray timing provide unprecedented access to the equation of state (EoS) of cold dense matter at densities difficult to realize in terrestrial experiments. At the same time, predictions for the EoS with reliable uncertainty estimates from chiral effective field theory ($\chi$EFT) bound our theoretical ignorance. In this work, we analyze astrophysical data using a nonparametric representation of the neutron-star EoS conditioned on $\chi$EFT to directly constrain the underlying physical properties of the compact objects. We discuss how the data alone constrain the EoS at high densities when we condition on $\chi$EFT at low densities. We also demonstrate how to exploit astrophysical data to directly test the predictions of $\chi$EFT for the EoS up to twice nuclear saturation density, and estimate the density at which these predictions might break down. We find that the existence of massive pulsars, gravitational waves from GW170817, and NICER observations of PSR J0030+0451 favor $\chi$EFT predictions for the EoS up to nuclear saturation density over a more agnostic analysis by as much as a factor of 7 for the quantum Monte Carlo (QMC) calculations used in this work. While $\chi$EFT predictions using QMC are fully consistent with gravitational-wave data up to twice nuclear saturation density, NICER observations suggest that the EoS stiffens relative to these predictions at nuclear saturation density. Additionally, we marginalize over the uncertainty in the density at which $\chi$EFT begins to break down, constraining the radius of a $1.4\,M_\odot$ neutron star to $R_{1.4}=11.40^{+1.38}_{-1.04}$ ($12.54^{+0.71}_{-0.63}$) km and the pressure at twice nuclear saturation density to $p(2n_\mathrm{sat})=14.2^{+18.1}_{-8.4}$ ($28.7^{+15.3}_{-15.0}$) MeV/fm$^3$ with massive pulsar and gravitational-wave (and NICER) data.

Looking for Dark Matter in Neutron Star Light

Radio observations of neutron stars demonstrate a way to search for axions through their expected conversion to electromagnetic waves in a star’s magnetic field.

Curvature-slope correlation of nuclear symmetry energy and its imprints on the crust-core transition, radius, and tidal deformability of canonical neutron stars

Background: The nuclear symmetry energy $E_{sym}(\rho)$ encodes information about the energy necessary to make nuclear systems more neutron-rich. While its slope parameter L at the saturation density $\rho_0$ of nuclear matter has been relatively well constrained by recent astrophysical observations and terrestrial nuclear experiments, its curvature $K_{\rm{sym}}$ characterizing the $E_{sym}(\rho)$ around $2\rho_0$ remains largely unconstrained. Over 520 calculations for $E_{sym}(\rho)$ using various nuclear theories and interactions in the literature have predicted several significantly different $K_{\rm{sym}}-L$ correlations. Purpose: If a unique $K_{\rm{sym}}-L$ correlation of $E_{sym}(\rho)$ can be firmly established, it will enable us to progressively better constrain the high-density behavior of $E_{sym}(\rho)$ using the available constraints on its slope parameter L. We investigate if and by how much the different $K_{\rm{sym}}-L$ correlations may affect neutron star observables. Method: A meta-model of nuclear Equation of States (EOSs) with three representative $K_{\rm{sym}}-L$ correlation functions is used to generate multiple EOSs for neutron stars. We then examine effects of the $K_{\rm{sym}}-L$ correlation on the crust-core transition density and pressure as well as the radius and tidal deformation of canonical neutron stars. Results:The $K_{\rm{sym}}-L$ correlation affects significantly both the crust-core transition density and pressure. It also has strong imprints on the radius and tidal deformability of canonical neutron stars especially at small L values. The available data from LIGO/VIRGO and NICER set some useful limits for the slope L but can not distinguish the three representative $K_{\rm{sym}}-L$ correlations considered.

Quantum nucleation of up-down quark matter and astrophysical implications

Quark matter with only $u$ and $d$ quarks ($ud$QM) might be the ground state of baryonic matter at large baryon number $A>A_{\rm min}$. With $A_{\rm min}\gtrsim 300$, this has no direct conflict with the stability of ordinary nuclei. An intriguing test of this scenario is to look for quantum nucleation of $ud$QM inside neutron stars due to their large baryon densities. In this paper, we study the transition rate of cold neutron stars to $ud$ quark stars ($ud$QSs) and the astrophysical implications, considering the relevant theoretical uncertainties and observational constraints. It turns out that a large portion of parameter space predicts an instantaneous transition, and so the observed neutron stars are mostly $ud$QSs. We find this possibility still viable under the recent gravitational wave and pulsar observations, although there are debates on its compatibility with some observations that involve complicated structure of quark matter. The tension could be partially relieved in the two-families scenario, where the high-mass stars ($M\gtrsim2 M_{\odot}$) are all $ud$QSs and the low-mass ones ($M\sim1.4\, M_{\odot}$) are mostly hadronic stars. In this case, the slow transition of the low-mass hadronic stars points to a very specific class of hadronic models with moderately stiff EOSs, and $ud$QM properties are also strongly constrained.

Dark matter capture in celestial objects: light mediators, self-interactions, and complementarity with direct detection

We generalize the formalism for DM capture in celestial bodies to account for arbitrary mediator mass, and update the existing and projected astrophysical constraints on DM-nucleon scattering cross section from observations of neutron stars. We show that the astrophysical constraints on the DM-nucleon interaction strength, that were thought to be the most stringent, drastically weaken for light mediators and can be completely voided. For asymmetric DM, existing astrophysical constraints are completely washed out for mediators lighter than 5 MeV, and for annihilating DM the projected constraints are washed out for mediators lighter than 0.25 MeV. Related terrestrial direct detection bounds also weaken, but in a complementary fashion; they supersede the astrophysical capture bounds for small or large DM mass, respectively for asymmetric or annihilating DM. Repulsive self-interactions of DM have an insignificant impact on the total capture rate, but a significant impact on the black hole formation criterion. This further weakens the constraints on DM-nucleon interaction strength for asymmetric self-repelling DM, whereas constraints remain unaltered for annihilating self-repelling DM. We use the correct Hawking evaporation rate of the newly formed black hole, that was approximated as a blackbody in previous studies, and show that, despite a more extensive alleviation of collapse as a result, the observation of a neutron star collapse can probe a wide range of DM self-interaction strengths.

Equation of State Constraints from Nuclear Physics, Neutron Star Masses, and Future Moment of Inertia Measurements

Abstract We explore constraints on the equation of state (EOS) of neutron-rich matter based on microscopic calculations up to nuclear densities and observations of neutron stars. In a previous work we showed that predictions based on modern nuclear interactions derived within chiral effective field theory and the observation of two-solar-mass neutron stars result in a robust uncertainty range for neutron star radii and the EOS over a wide range of densities. In this work we extend this study, employing both the piecewise polytrope extension from Hebeler et al. as well as the speed of sound model of Greif et al., and show that moment of inertia measurements of neutron stars can significantly improve the constraints on the EOS and neutron star radii.