Constraining the speed of sound inside neutron stars with chiral effective field theory interactions and observations

The dense matter equation of state (EOS) determines neutron star (NS) structure but can be calculated reliably only up to one to two times the nuclear saturation density, using accurate many-body methods that employ nuclear interactions from chiral effective field theory constrained by scattering data. In this work, we use physically motivated ansatzes for the speed of sound $c_S$ at high density to extend microscopic calculations of neutron-rich matter to the highest densities encountered in stable NS cores. We show how existing and expected astrophysical constraints on NS masses and radii from X-ray observations can constrain the speed of sound in the NS core. We confirm earlier expectations that $c_S$ is likely to violate the conformal limit of $c_S^2\leq c^2/3$, possibly reaching values closer to the speed of light $c$ at a few times the nuclear saturation density, independent of the nuclear Hamiltonian. If QCD obeys the conformal limit, we conclude that the rapid increase of $c_S$ required to accommodate a $2$ M$_\odot$ NS suggests a form of strongly interacting matter where a description in terms of nucleons will be unwieldy, even between one and two times the nuclear saturation density. For typical NSs with masses in the range $1.2-1.4~$ M$_\odot$, we find radii between $10$ and $14$ km, and the smallest possible radius of a $1.4$ M$_{\odot}$ NS consistent with constraints from nuclear physics and observations is $8.4$ km. We also discuss how future observations could constrain the EOS and guide theoretical developments in nuclear physics.Comment: 24 pages, 14 figures, published versio

Tabulated Equations of State From Models Informed by Chiral Effective Field Theory

We construct four equation of state (EoS) tables, tabulated over a range of temperatures, densities, and charge fractions, relevant for neutron star applications such as simulations of neutron star mergers. The EoS are computed from a relativistic mean-field theory constrained by the pure neutron matter EoS from chiral effective field theory, inferred properties of isospin-symmetric nuclear matter, and astrophysical observations of neutron star structure. To model nuclear matter at low densities, we attach an EoS that models inhomogeneous nuclear matter at arbitrary temperatures and charge fractions. The four EoS tables we develop are available from the CompOSE EoS repository compose.obspm.fr/eos/297 and gitlab.com/ahaber/qmc-rmf-tables.Comment: 8 pages, 3 figures v2: version published in Phys. Scr. 98, added finite T chiral EFT compariso

Relativistic mean-field theories for neutron-star physics based on chiral effective field theory

We describe and implement a procedure for determining the couplings of a Relativistic Mean-Field Theory (RMFT) that is optimized for application to neutron star phenomenology. In the standard RMFT approach, the couplings are constrained by comparing the theory's predictions for symmetric matter at saturation density with measured nuclear properties. The theory is then applied to neutron stars which consist of neutron-rich matter at densities ranging up to several times saturation density, which allows for additional astrophysical constraints. In our approach, rather than using the RMFT to extrapolate from symmetric to neutron-rich matter and from finite-sized nuclei to uniform matter, we fit the RMFT to properties of uniform pure neutron matter obtained from chiral effective field theory. Chiral effective field theory incorporates the experimental data for nuclei in the framework of a controlled expansion for nuclear forces valid at nuclear densities and enables us to account for theoretical uncertainties when fitting the RMFT. We construct four simple RMFTs that span the uncertainties provided by chiral effective field theory for neutron matter, and are consistent with current astrophysical constraints on the equation of state. Our RMFTs can be used to model the properties of neutron-rich matter across the vast range of densities and temperatures encountered in neutron stars and their mergers.Comment: V2: Minor corrections, version published in PRC. 12 pages and 5 figure

Maximally local two-nucleon interactions at fourth order in delta-less chiral effective field theory

We present new maximally-local two-nucleon interactions derived in delta-less chiral effective field theory up to next-to-next-to-next-to leading order that include all contact and pion-exchange contributions to the nuclear Hamiltonian up to this order. Our interactions are fit to nucleon-nucleon phase shifts using a Bayesian statistical approach, and explore a wide cutoff range from $0.6-0.9$ fm ($\sim$ 440 MeV to $\sim$ 660 MeV). These interactions can be straightforwardly employed in accurate quantum Monte Carlo methods, such as the auxiliary field diffusion Monte Carlo method. Together with local three-nucleon forces, calculations with these new interactions will provide improved benchmarks for the structure of atomic nuclei and serve as crucial input to analyses of exciting astrophysical phenomena involving neutron stars, such as binary neutron-star mergers.Comment: 19 pages, 8 figures. Comments welcom

Parameter estimation for strong phase transitions in supranuclear matter using gravitational-wave astronomy

At supranuclear densities, explored in the core of neutron stars, a strong phase transition from hadronic matter to more exotic forms of matter might be present. To test this hypothesis, binary neutron-star mergers offer a unique possibility to probe matter at densities that we can not create in any existing terrestrial experiment. In this work, we show that, if present, strong phase transitions can have a measurable imprint on the binary neutron-star coalescence and the emitted gravitational-wave signal. We construct a new parameterization of the supranuclear equation of state that allows us to test for the existence of a strong phase transition and extract its characteristic properties purely from the gravitational-wave signal of the inspiraling neutron stars. We test our approach using a Bayesian inference study simulating 600 signals with three different equations of state and find that for current gravitational-wave detector networks already twelve events might be sufficient to verify the presence of a strong phase transition. Finally, we use our methodology to analyze GW170817 and GW190425, but do not find any indication that a strong phase transition is present at densities probed during the inspiral.Comment: 17 pages, 11 figure

Modeling Solids in Nuclear Astrophysics with Smoothed Particle Hydrodynamics

Smoothed Particle Hydrodynamics (SPH) is a frequently applied tool in computational astrophysics to solve the fluid dynamics equations governing the systems under study. For some problems, for example when involving asteroids and asteroid impacts, the additional inclusion of material strength is necessary in order to accurately describe the dynamics. In compact stars, that is white dwarfs and neutron stars, solid components are also present. Neutron stars have a solid crust which is the strongest material known in nature. However, their dynamical evolution, when modeled via SPH or other computational fluid dynamics codes, is usually described as a purely fluid dynamics problem. Here, we present the first 3D simulations of neutron-star crustal toroidal oscillations including material strength with the Los Alamos National Laboratory SPH code FleCSPH. In the first half of the paper, we present the numerical implementation of solid material modeling together with standard tests. The second half is on the simulation of crustal oscillations in the fundamental toroidal mode. Here, we dedicate a large fraction of the paper to approaches which can suppress numerical noise in the solid. If not minimized, the latter can dominate the crustal motion in the simulations.Comment: 24 pages, 29 figure

On the Nature of GW190814 and Its Impact on the Understanding of Supranuclear Matter

The observation of a compact object with a mass of 2.50-2.67Me on 2019 August 14, by the LIGO Scientific and Virgo collaborations (LVC) has the potential to improve our understanding of the supranuclear equation of state. While the gravitational-wave analysis of the LVC suggests that GW190814 likely was a binary black hole system, the secondary component could also have been the heaviest neutron star observed to date. We use our previously derived nuclear-physics-multimessenger astrophysics framework to address the nature of this object. Based on our findings, we determine GW190814 to be a binary black hole merger with a probability of >99.9%. Even if we weaken previously employed constraints on the maximum mass of neutron stars, the probability of a binary black hole origin is still ∼81%. Furthermore, we study the impact that this observation has on our understanding of the nuclear equation of state by analyzing the allowed region in the mass-radius diagram of neutron stars for both a binary black hole or neutron star-black hole scenario. We find that the unlikely scenario in which the secondary object was a neutron star requires rather stiff equations of state with a maximum speed of sound cs ≥0.6 times the speed of light, while the binary black hole scenario does not offer any new insight