Recent Neutron Star Observations Research Articles

Constraints onBose–Einsteincondensate stars as neutron stars models from new observational data

AbstractWe evaluate the feasibility of Bose–Einstein condensate stars (BECS) as models for the interior of neutron stars (NSs). BECS are compact objects composed of bosons, formed through the spin‐parallel pairing of neutrons. Here, we utilize the astronomical data from GW170817, XMMU J173203.3‐344518 (the lightest NS known), and a novel lower limit on NS core heat capacity to scrutinize the compatibility of BECS with these recent observations of NSs. Our specific focus is to constrain the values of the scattering length , parameter determining the strength of particle interactions in the model. Our analysis suggests that if the stars involved in GW170817 were BECSs, the scattering length of their constituent bosons should fall within the to fm range. Additionally, at a scattering length of fm, stars with mass and radius characteristics akin to XMMU J173203.3‐344518 are identified. Moreover, we find that the heat capacity depends on the mass and temperature of BECS, and surpasses the established lower bound for NS cores when fm. In summary, our results endorse BECS models with fm, providing NS observables in agreement with diverse observations and contributing to the understanding of NS interiors.

Chiral condensates for neutron stars in hadron-quark crossover; from a parity doublet nucleon model to an NJL quark model

In this contribution, we summarize our recent studies on the chiral invariant mass and the chiral condensates in neutron star matter. We construct a unified equations of state assuming the crossover phase transition from hadronic matter described by a parity doublet model to quark matter by an Nambu--Jona-Lasinio type quark model. We first show that the chiral invariant mass is constrained to be $600\,\mbox{MeV}\, \lesssim\, m_0 \,\lesssim \, 900\,\mbox{MeV}$ from recent observations of neutron stars. We then determine the density dependence of the chiral condensate in the crossover description, and show that the chiral condensates are actually smoothly connected from the hadronic matter where the change is driven by the positive chiral scalar charge in a nucleon, to the quark matter where the change is by the modification of the quark Dirac sea, reflecting the hadron-quark crossover.

Extremal instability for topological black holes

The initial idea underlying the Weak Gravity Conjecture is that extremal black holes must always be “unstable”, in the sense that they should slowly decay by emitting either particles or smaller black holes. Here we show that, when this idea is applied to the planar asymptotically AdS black holes which play a central role in applications of holography, the result, via gauge-gravity duality, is a prediction that there should exist a lower bound on the possible densities of cold strongly coupled matter. Recent observations of neutron stars suggest that, in many cases, even the extreme densities in their cores may not be sufficient to generate quark matter, showing that there is indeed a (very high) lower bound on the possible density of cold quark matter.

Translating Neutron Star Observations to Nuclear Symmetry Energy via Deep Neural Networks

One of the most significant challenges involved in efforts to understand the equation of state of dense neutron-rich matter is the uncertain density dependence of the nuclear symmetry energy. In particular, the nuclear symmetry energy is still rather poorly constrained, especially at high densities. On the other hand, detailed knowledge of the equation of state is critical for our understanding of many important phenomena in the nuclear terrestrial laboratories and the cosmos. Because of its broad impact, pinning down the density dependence of the nuclear symmetry energy has been a long-standing goal of both nuclear physics and astrophysics. Recent observations of neutron stars, in both electromagnetic and gravitational-wave spectra, have already constrained significantly the nuclear symmetry energy at high densities. The next generation of telescopes and gravitational-wave observatories will provide an unprecedented wealth of detailed observations of neutron stars, which will improve further our knowledge of the density dependence of nuclear symmetry energy, and the underlying equation of state of dense neutron-rich matter. Training deep neural networks to learn a computationally efficient representation of the mapping between astrophysical observables of neutron stars, such as masses, radii, and tidal deformabilities, and the nuclear symmetry energy allows its density dependence to be determined reliably and accurately. In this work, we use a deep learning approach to determine the nuclear symmetry energy as a function of density directly from observational neutron star data. We show, for the first time, that artificial neural networks can precisely reconstruct the nuclear symmetry energy from a set of available neutron star observables, such as masses and radii as measured by, e.g., the NICER mission, or masses and tidal deformabilities as measured by the LIGO/VIRGO/KAGRA gravitational-wave detectors. These results demonstrate the potential of artificial neural networks to reconstruct the symmetry energy and the equation of state directly from neutron star observational data, and emphasize the importance of the deep learning approach in the era of multi-messenger astrophysics.

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

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.

Analytic I-Love-C relations for realistic neutron stars

Recent observations of neutron stars (NS) with radio, X-rays and gravitational waves have begun to constrain the equation of state (EoS) for nuclear matter beyond the nuclear saturation density. To one's surprise, there exist approximate universal relations connecting certain bulk properties of NSs that are insensitive to the underlying EoS and having important applications on probing fundamental physics including nuclear and gravitational physics. To date, analytic works on universal relations for realistic NSs are lacking, which may lead to a better understanding of the origin of the universality. Here, we focus on the universal relations between the compactness(C), moment of inertia (I), and tidal deformability(related to the Love number), and derive analytic, approximate I-Love-C relations. To achieve this, we construct slowly-rotating/ tidally-deformed NS solutions analytically starting from an extended Tolman VII model that accurately describes non-rotating realistic NSs, allowing us to extract the moment of inertia and the tidal deformability on top of the compactness. We solve the field equations analytically by expanding them about the Newtonian limit and keeping up to 6th order in the stellar compactness. Based on these analytic solutions, we can mathematically demonstrate the O(10%) EoS variation in the I-C and Love-C relations and the O(1%) variation in the I-Love relation that have previously been found numerically. Our new analytic relations agree more accurately with numerical results for realistic NSs(especially the I-C and Love-C ones) than the analytic relations for constant density stars derived in previous work. Based on these analytic findings, we attribute a possible origin of the universality for the I-C and Love-C relations to the fact that the energy density of realistic NSs can be approximated as a quadratic function, as is the case for the Tolman VII solution.

Effects of Symmetry Energy on the Equation of State for Simulations of Core-collapse Supernovae and Neutron-star Mergers

Abstract We construct a new equation of state (EOS) for numerical simulations of core-collapse supernovae and neutron-star mergers based on an extended relativistic mean-field model with a small symmetry energy slope L, which is compatible with both experimental nuclear data and recent observations of neutron stars. The new EOS table (EOS4) based on the extended TM1 (TM1e) model with L = 40 MeV is designed in the same tabular form and compared with the commonly used Shen EOS (EOS2) based on the original TM1 model with . This is convenient and useful for performing numerical simulations and examining the influences of symmetry energy and its density dependence on astrophysical phenomena. In comparison with the TM1 model used in EOS2, the TM1e model provides a similar maximum neutron-star mass but smaller radius and tidal deformability for a 1.4M ⊙ neutron star, which is more consistent with current constraints. By comparing the phase diagram and thermodynamic quantities between EOS4 and EOS2, it is found that the TM1e model predicts a relatively larger region of nonuniform matter and softer EOS for neutron-rich matter. Significant differences between EOS4 and EOS2 are observed in the case with low proton fraction, while the properties of symmetric matter remain unchanged.

Bounds on the speed of sound in dense matter, and neutron star structure

The accurate determination of the maximum mass of the neutron stars is one of the most important tasks in astrophysics. It is directly related to the identification of the black holes in the universe, the production of neutron stars from the supernovae explosion, and the equation of state (EoS) of dense matter. However, not only the EoS is directly connected with neutron star masses, but also the speed of sound in dense matter is a crucial quantity which characterizes the stiffness of the EoS. The upper bound of the speed of sound imposes strong constraints on the maximum mass of neutron stars. However, this upper bound remains still an open issue. Recent observations, of binary neutron star systems, offer the possibility of measuring with high accuracy both the mass and the tidal polarizability of the stars. We study possible effects of the upper bound of the speed of sound on the upper bound of the mass and the tidal polarizability. We conclude that these kinds of measurements, combined with recent observations of neutron stars with masses close to $2 M_{\odot}$, will provide robust constraints on the equation of state of hadronic matter at high densities.

Density dependence of nuclear symmetry energy

High density behavior of nuclear symmetry energy is studied on the basis of the stiffest density dependence of asymmetric contribution to energy per nucleon in charge neutral n + p + e + [Formula: see text] matter under beta equilibrium. The density dependence of nuclear symmetry energy obtained in this way is neither very stiff nor soft at high densities and is found to be in conformity with recent observations of neutron stars.

A NEW CLASS OF NEUTRON STAR BINARIES AND ITS IMPLICATIONS

Recent discovery of 2M ⊙ neutron stars in white dwarf-neutron star binaries, PSR J1614-2230 and PSR J0348+0432, has given strong constraints on the maximum mass of neutron stars. On the other hand, all well-measured neutron star masses in double neutron star binaries are still less than 1.5M ⊙ . These observations suggest that the neutron star masses in binaries may depend on the evolution process of neutron star binaries. In addition, recent works on LMXB (low-mass X-ray binaries) provides us the possibility of estimating the masses and radii of accreting neutron stars in LMXBs. In this talk, we discuss the implications of recent neutron star observations to the neutron star equation of states and the related astrophysical problems. For the evolution of neutron star binaries, we also discuss the possibilities of super-Eddington accretion onto the primary neutron stars.

The equation of state of neutron matter, symmetry energy and neutron star structure

We review the calculation of the equation of state of pure neutron matter using quantum Monte Carlo (QMC) methods. QMC algorithms permit the study of many-body nuclear systems using realistic two- and three-body forces in a non-perturbative framework. We present the results for the equation of state of neutron matter, and focus on the role of three-neutron forces at supranuclear density. We discuss the correlation between the symmetry energy, the neutron star radius and the symmetry energy. We also combine QMC and theoretical models of the three-nucleon interactions, and recent neutron star observations to constrain the value of the symmetry energy and its density dependence.

How neutron stars constrain the nuclear equation of state

Recent neutron star observations set new constraints for the equation of state of baryonic matter. A chiral effective field theory approach is used for the description of neutron-dominated nuclear matter present in the outer core of neutron stars. Possible hybrid stars with quark matter in the inner core are discussed using a three-flavor Nambu–Jona-Lasinio model.

CORE-COLLAPSE SUPERNOVA EQUATIONS OF STATE BASED ON NEUTRON STAR OBSERVATIONS

Many of the currently available equations of state for core-collapse supernova simulations give large neutron star radii and do not provide large enough neutron star masses, both of which are inconsistent with some recent neutron star observations. In addition, one of the critical uncertainties in the nucleon-nucleon interaction, the nuclear symmetry energy, is not fully explored by the currently available equations of state. In this article, we construct two new equations of state which match recent neutron star observations and provide more flexibility in studying the dependence on nuclear matter properties. The equations of state are also provided in tabular form, covering a wide range in density, temperature and asymmetry, suitable for astrophysical simulations. These new equations of state are implemented into our spherically symmetric core-collapse supernova model, which is based on general relativistic radiation hydrodynamics with three-flavor Boltzmann neutrino transport. The results are compared with commonly used equations of state in supernova simulations of 15 and 40 solar mass progenitors. We do not find any simple correlations between individual nuclear matter properties at saturation and the outcome of these simulations. However, the new equations of state lead to the most compact neutron stars among the relativistic mean-field models which we considered. The new models also obey the previously observed correlation between the time to black hole formation and the maximum mass of an s=4 neutron star.

Hyperons and massive neutron stars: Vector repulsion and strangeness

Finding the composition of cold dense matter is still a serious challenge. Recent neutron star observations, especially the observed mass of pulsar PSR J1614-2230 with M=(1.97±0.04)M⊙, constrain the equation of state (EoS) of hadronic matter and question the existence of hyperons in the core of neutron stars. As in previous work, we investigate the conditions for the presence of hyperons in the neutron star core by exploring the vector baryon–meson coupling assuming flavor SU(3) symmetry. However, as an extension to our work, we investigate the full range of SU(3) parameters which had previously been kept fixed.

The equation of state of neutron star matter and the symmetry energy

We present an overview of microscopic calculations of the Equation of State (EOS) of neutron matter performed using Quantum Monte Carlo techniques. We focus on the role of the model of the three-neutron force in the high-density part of the EOS up to a few times the saturation density. We also discuss the interplay between the symmetry energy and the neutron star mass-radius relation.The combination of theoretical models of the EOS with recent neutron star observations permits us to constrain the value of the symmetry energy and its slope. We show that astrophysical observations are starting to provide important insights into the properties of neutron star matter.

NEUTRON STARS WITH SMALL RADII—THE ROLE OF Δ RESONANCES

Recent neutron star observations suggest that the masses and radii of neutron stars may be smaller than previously considered, which would disfavor a purely nucleonic equation of state. In our model, we use a the flavor SU(3) sigma model that includes delta resonances and hyperons in the equation of state. We find that if the coupling of the delta resonances to the vector mesons is slightly smaller than that of the nucleons, we can reproduce both the measured mass-radius relationship and the extrapolated equation of state.

Recent progress constraining the nuclear equation of state from astrophysics and heavy-ion reactions

The quest for the nuclear equation of state (EoS) at high densities and/or extreme isospin is one of the longstanding problems of nuclear physics. Ab initio calculations for the nuclear many-body problem make predictions for the density and isospin dependence of the EoS far away from the saturation point of nuclear matter. On the other hand, in recent years substantial progress has been made to constrain the EoS, both from the astrophysical side and from accelerator-based experiments. Heavy-ion experiments support a soft EoS at moderate densities, while recent neutron star observations require a ‘stiff’ high density behavior. Both constraints are discussed and shown to be in agreement with the predictions from the many-body theory.