Showing 1–41 of 41 results for author: Tews, I

Towards accelerated nuclear-physics parameter estimation from binary neutron star mergers: Emulators for the Tolman-Oppenheimer-Volkoff equations

Authors: Brendan T. Reed, Rahul Somasundaram, Soumi De, Cassandra L. Armstrong, Pablo Giuliani, Collin Capano, Duncan A. Brown, Ingo Tews

Abstract: Gravitational-wave observations of binary neutron-star (BNS) mergers have the potential to revolutionize our understanding of the nuclear equation of state (EOS) and the fundamental interactions that determine its properties. However, Bayesian parameter estimation frameworks do not typically sample over microscopic nuclear-physics parameters that determine the EOS. One of the major hurdles in doin… ▽ More Gravitational-wave observations of binary neutron-star (BNS) mergers have the potential to revolutionize our understanding of the nuclear equation of state (EOS) and the fundamental interactions that determine its properties. However, Bayesian parameter estimation frameworks do not typically sample over microscopic nuclear-physics parameters that determine the EOS. One of the major hurdles in doing so is the computational cost involved in solving the neutron-star structure equations, known as the Tolman-Oppenheimer-Volkoff (TOV) equations. In this paper, we explore approaches to emulating solutions for the TOV equations: Multilayer Perceptrons (MLP), Gaussian Processes (GP), and a data-driven variant of the reduced basis method (RBM). We implement these emulators for three different parameterizations of the nuclear EOS, each with a different degree of complexity represented by the number of model parameters. We find that our MLP-based emulators are generally more accurate than the other two algorithms whereas the RBM results in the largest speedup with respect to the full, high-fidelity TOV solver. We employ these emulators for a simple parameter inference using a potentially loud BNS observation, and show that the posteriors predicted by our emulators are in excellent agreement with those obtained from the full TOV solver. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 13 pages and 9 figures. Comments Welcome

Report number: LA-UR-24-25009

An overview of existing and new nuclear and astrophysical constraints on the equation of state of neutron-rich dense matter

Authors: Hauke Koehn, Henrik Rose, Peter T. H. Pang, Rahul Somasundaram, Brendan T. Reed, Ingo Tews, Adrian Abac, Oleg Komoltsev, Nina Kunert, Aleksi Kurkela, Michael W. Coughlin, Brian F. Healy, Tim Dietrich

Abstract: Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. Here, we apply these different pieces of data individually to a broad set of physics-agnostic candidate EOSs and analyze the resulting constraints. Specifically, we make use of… ▽ More Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. Here, we apply these different pieces of data individually to a broad set of physics-agnostic candidate EOSs and analyze the resulting constraints. Specifically, we make use of information from chiral effective field theory, perturbative quantum chromodynamics, as well as data from heavy-ion collisions and the PREX-II and CREX experiments. We also investigate the impact of current mass and radius measurements of neutron stars, such as radio timing measurements of heavy pulsars, NICER data, and other X-ray observations. We augment these by reanalyses of the gravitational-wave (GW) signal GW170817, its associated kilonova AT2017gfo and gamma-ray burst afterglow, the GW signal GW190425, and the GRB211211A afterglow, where we use improved models for the tidal waveform and kilonova light curves. Additionally, we consider the postmerger fate of GW170817 and its consequences for the EOS. This large and diverse set of constraints is eventually combined in numerous ways to explore limits on quantities such as the typical neutron-star radius, the maximum neutron-star mass, the nuclear symmetry-energy parameters, and the speed of sound. Based on the priors from our EOS candidate set, we find the radius of the canonical 1.4 M$_\odot$ neutron star to be $R_{1.4}= 12.27_{-0.94}^{+0.83}$ km and the TOV mass $M_{\rm TOV}= 2.26_{-0.22}^{+0.45}$ M$_\odot$ at 95% credibility, when including those constraints where systematic uncertainties are deemed small. A less conservative approach, combining all the presented constraints, similarly yields $R_{1.4}= 12.20_{-0.50}^{+0.53}$ km and $M_{\rm TOV}= 2.31_{-0.20}^{+0.08}$ M$_\odot$. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: 49 pages, 32 figures, webinterface for custom constraint combinations in https://enlil.gw.physik.uni-potsdam.de/eos_constraints

Report number: LA-UR-24-20420

Equation of state at neutron-star densities and beyond from perturbative QCD

Authors: Oleg Komoltsev, Rahul Somasundaram, Tyler Gorda, Aleksi Kurkela, Jérôme Margueron, Ingo Tews

Abstract: We explore the consequences of imposing robust thermodynamic constraints arising from perturbative Quantum Chromodynamics (QCD) when inferring the dense-matter equation-of-state (EOS). We find that the termination density, up to which the EOS modeling is performed in an inference setup, strongly affects the constraining power of the QCD input. This sensitivity in the constraining power arises from… ▽ More We explore the consequences of imposing robust thermodynamic constraints arising from perturbative Quantum Chromodynamics (QCD) when inferring the dense-matter equation-of-state (EOS). We find that the termination density, up to which the EOS modeling is performed in an inference setup, strongly affects the constraining power of the QCD input. This sensitivity in the constraining power arises from EOSs that have a specific form, with drastic softening immediately above the termination density followed by a strong stiffening. We also perform explicit modeling of the EOS down from perturbative-QCD densities to construct a new QCD likelihood function that incorporates additional perturbative-QCD calculations of the sound speed and is insensitive to the termination density, which we make publicly available. △ Less

Submitted 1 February, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: 16 pages, 11 figures. Version 2: added Appendix A, corrected minor error in the marginalization code that introduced small noise in results; updated figures accordingly. The link to the publicly available Marginalized QCD likelihood function https://doi.org/10.5281/zenodo.10592568

Report number: LA-UR-23-32625

What can we learn about the unstable equation-of-state branch from neutron-star mergers?

Authors: Maximiliano Ujevic, Rahul Somasundaram, Tim Dietrich, Jerome Margueron, Ingo Tews

Abstract: The Equation of State (EOS) of dense strongly-interacting matter can be probed by astrophysical observations of neutron stars (NS), such as X-ray detections of pulsars or the measurement of the tidal deformability of NSs during the inspiral stage of NS mergers. These observations constrain the EOS at most up to the density of the maximum-mass configuration, $n_\textrm{TOV}$, which is the highest d… ▽ More The Equation of State (EOS) of dense strongly-interacting matter can be probed by astrophysical observations of neutron stars (NS), such as X-ray detections of pulsars or the measurement of the tidal deformability of NSs during the inspiral stage of NS mergers. These observations constrain the EOS at most up to the density of the maximum-mass configuration, $n_\textrm{TOV}$, which is the highest density that can be explored by stable NSs for a given EOS. However, under the right circumstances, binary neutron star (BNS) mergers can create a postmerger remnant that explores densities above $n_\textrm{TOV}$. In this work, we explore whether the EOS above $n_\textrm{TOV}$ can be measured from gravitational-wave or electromagnetic observations of the postmerger remnant. We perform a total of twenty-five numerical-relativity simulations of BNS mergers for a range of EOSs and find no case in which different descriptions of the matter above $n_{\rm TOV}$ have a detectable impact on postmerger observables. Hence, we conclude that the EOS above $n_\textrm{TOV}$ can likely not be probed through BNS merger observations for the current and next generation of detectors. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 8 pages, 4 figures

Report number: LA-UR-23-31683

Probing Quarkyonic Matter in Neutron Stars with the Bayesian Nuclear-Physics Multi-Messenger Astrophysics Framework

Authors: Peter T. H. Pang, Lars Sivertsen, Rahul Somasundaram, Tim Dietrich, Srimoyee Sen, Ingo Tews, Michael Coughlin, Chris Van Den Broeck

Abstract: The interior of neutron stars contains matter at the highest densities realized in our Universe. Interestingly, theoretical studies of dense matter, in combination with the existence of two solar mass neutron stars, indicate that the speed of sound $c_s$ has to increase to values well above the conformal limit ($c_s^2\sim 1/3$) before decreasing again at higher densities. The decrease could be exp… ▽ More The interior of neutron stars contains matter at the highest densities realized in our Universe. Interestingly, theoretical studies of dense matter, in combination with the existence of two solar mass neutron stars, indicate that the speed of sound $c_s$ has to increase to values well above the conformal limit ($c_s^2\sim 1/3$) before decreasing again at higher densities. The decrease could be explained by either a strong first-order phase transition or a cross-over transition from hadronic to quark matter. The latter scenario leads to a pronounced peak in the speed of sound reaching values above the conformal limit, naturally explaining the inferred behavior. In this work, we use the Nuclear-Physics Multi-Messenger Astrophysics framework \textsc{NMMA} to compare predictions of the quarkyonic matter model with astrophysical observations of neutron stars, with the goal of constraining model parameters. Assuming quarkyonic matter to be realized within neutron stars, we find that there can be a significant amount of quarks inside the core of neutron stars with masses in the two solar mass range, amounting to up to $\sim 0.13M_\odot$, contributing $\sim 5.9\%$ of the total mass. Furthermore, for the quarkyonic matter model investigated here, the radius of a $1.4M_\odot$ neutron star would be $13.44^{+1.69}_{-1.54} (13.54^{+1.02}_{-1.04})$ km, at $95\%$ credibility, without (with) the inclusion of AT2017gfo. △ Less

Submitted 28 February, 2024; v1 submitted 29 August, 2023; originally announced August 2023.

Comments: 14 pages, 10 figures

Report number: LA-UR-23-29118, NP3M-P2300017

Journal ref: Phys. Rev. C 109, 025807 (2024)

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

Authors: Mark G. Alford, Liam Brodie, Alexander Haber, Ingo Tews

Abstract: 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… ▽ More 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. △ Less

Submitted 17 November, 2023; v1 submitted 16 April, 2023; originally announced April 2023.

Comments: 8 pages, 3 figures v2: version published in Phys. Scr. 98, added finite T chiral EFT comparison

Report number: LA-UR-23-23837

Journal ref: Phys. Scr. 98 125302 (2023)

Reproducing the results for NICER observation of PSR J0030+0451

Authors: Chaitanya Afle, Patrick R. Miles, Silvina Caino-Lores, Collin D. Capano, Ingo Tews, Karan Vahi, Ewa Deelman, Michela Taufer, Duncan A. Brown

Abstract: NASA's Neutron Star Interior Composition Explorer (NICER) observed X-ray emission from the pulsar PSR J0030+0451 in 2018. Riley et al. reported Bayesian parameter measurements of the mass and the star's radius using pulse-profile modeling of the X-ray data. This paper reproduces their result using the open-source software X-PSI and publicly available data within expected statistical errors. We not… ▽ More NASA's Neutron Star Interior Composition Explorer (NICER) observed X-ray emission from the pulsar PSR J0030+0451 in 2018. Riley et al. reported Bayesian parameter measurements of the mass and the star's radius using pulse-profile modeling of the X-ray data. This paper reproduces their result using the open-source software X-PSI and publicly available data within expected statistical errors. We note the challenges we faced in reproducing the results and demonstrate that the analysis can be reproduced and reused in future works by changing the prior distribution for the radius and the sampler configuration. We find no significant change in the measurement of the mass and radius, demonstrating that the original result is robust to these changes. Finally, we provide a containerized working environment that facilitates third-party reproduction of the measurements of mass and radius of PSR J0030+0451 using the NICER observations. △ Less

Submitted 31 January, 2024; v1 submitted 3 April, 2023; originally announced April 2023.

Comments: 13 pages, 4 figures, 2 tables. Final version accepted for publication in Computing in Science & Engineering

Report number: LA-UR-22-29359

Revealing the strength of three-nucleon interactions with the Einstein Telescope

Authors: Henrik Rose, Nina Kunert, Tim Dietrich, Peter T. H. Pang, Rory Smith, Chris Van Den Broeck, Stefano Gandolfi, Ingo Tews

Abstract: Three-nucleon forces are crucial for the accurate description of nuclear systems, including dense matter probed in neutron stars. We explore nuclear Hamiltonians that reproduce two-nucleon scattering data and properties of light nuclei, but differ in the three-nucleon interactions among neutrons. While no significantly improved constraints can be obtained from current astrophysical data, we show t… ▽ More Three-nucleon forces are crucial for the accurate description of nuclear systems, including dense matter probed in neutron stars. We explore nuclear Hamiltonians that reproduce two-nucleon scattering data and properties of light nuclei, but differ in the three-nucleon interactions among neutrons. While no significantly improved constraints can be obtained from current astrophysical data, we show that observations of neutron star mergers by next-generation detectors like the proposed Einstein Telescope could provide strong evidence to distinguish between these Hamiltonians. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Report number: LA-UR-22-21458

Bayesian model selection for GRB 211211A through multi-wavelength analyses

Authors: Nina Kunert, Sarah Antier, Vsevolod Nedora, Mattia Bulla, Peter T. H. Pang, Shreya Anand, Michael Coughlin, Ingo Tews, Jennifer Barnes, Thomas Hussenot-Desenonges, Brian Healy, Theophile Jegou du Laz, Meili Pilloix, Weizmann Kiendrebeogo, Tim Dietrich

Abstract: Although GRB 211211A is one of the closest gamma-ray bursts (GRBs), its classification is challenging because of its partially inconclusive electromagnetic signatures. In this paper, we investigate four different astrophysical scenarios as possible progenitors for GRB~211211A: a binary neutron-star merger, a black-hole--neutron-star merger, a core-collapse supernova, and an r-process enriched core… ▽ More Although GRB 211211A is one of the closest gamma-ray bursts (GRBs), its classification is challenging because of its partially inconclusive electromagnetic signatures. In this paper, we investigate four different astrophysical scenarios as possible progenitors for GRB~211211A: a binary neutron-star merger, a black-hole--neutron-star merger, a core-collapse supernova, and an r-process enriched core collapse of a rapidly rotating massive star (a collapsar). We perform a large set of Bayesian multi-wavelength analyses based on different models describing these scenarios and priors to investigate which astrophysical scenarios and processes might be related to GRB~211211A. Our analysis supports previous studies in which the presence of an additional component, likely related to $r$-process nucleosynthesis, is required to explain the observed light curves of GRB~211211A, as it can not solely be explained as a GRB afterglow. Fixing the distance to about $350~\rm Mpc$, namely the distance of the possible host galaxy SDSS J140910.47+275320.8, we find a statistical preference for a binary neutron-star merger scenario. △ Less

Submitted 15 December, 2023; v1 submitted 5 January, 2023; originally announced January 2023.

Comments: 14 pages, 6 figures

Report number: LA-UR-22-32364

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 2, January 2024, Pages 3900-3911

Reverse phase transitions in binary neutron-star systems with exotic-matter cores

Authors: Maximiliano Ujevic, Henrique Gieg, Federico Schianchi, Swami Vivekanandji Chaurasia, Ingo Tews, Tim Dietrich

Abstract: Multi-messenger observations of binary neutron star mergers provide a unique opportunity to constrain the dense-matter equation of state. Although it is known from quantum chromodynamics that hadronic matter will undergo a phase transition to exotic forms of matter, e.g., quark matter, the onset density of such a phase transition cannot be computed from first principles. Hence, it remains an open… ▽ More Multi-messenger observations of binary neutron star mergers provide a unique opportunity to constrain the dense-matter equation of state. Although it is known from quantum chromodynamics that hadronic matter will undergo a phase transition to exotic forms of matter, e.g., quark matter, the onset density of such a phase transition cannot be computed from first principles. Hence, it remains an open question if such phase transitions occur inside isolated neutron stars or during binary neutron star mergers, or if they appear at even higher densities that are not realized in the Cosmos. In this article, we perform numerical-relativity simulations of neutron-star mergers and investigate scenarios in which the onset density of such a phase transition is exceeded in at least one inspiralling binary component. Our simulations reveal that shortly before the merger it is possible that such stars undergo a "reverse phase transition", i.e., densities decrease and the quark core inside the star disappears leaving a purely hadronic star at merger. After the merger, when densities increase once more, the phase transition occurs again and leads, in the cases considered in this work, to a rapid formation of a black hole. We compute the gravitational-wave signal and the mass ejection for our simulations of such scenarios and find clear signatures that are related to the post-merger phase transition, e.g., smaller ejecta masses due to the softening of the equation of state through the quark core formation. Unfortunately, we do not find measurable imprints of the reverse phase transition. △ Less

Submitted 5 January, 2023; v1 submitted 8 November, 2022; originally announced November 2022.

Comments: 12 pages, 9 figures

Report number: Report-no: LA-UR-22-31740

Long Range Plan: Dense matter theory for heavy-ion collisions and neutron stars

Authors: Alessandro Lovato, Travis Dore, Robert D. Pisarski, Bjoern Schenke, Katerina Chatziioannou, Jocelyn S. Read, Philippe Landry, Pawel Danielewicz, Dean Lee, Scott Pratt, Fabian Rennecke, Hannah Elfner, Veronica Dexheimer, Rajesh Kumar, Michael Strickland, Johannes Jahan, Claudia Ratti, Volodymyr Vovchenko, Mikhail Stephanov, Dekrayat Almaalol, Gordon Baym, Mauricio Hippert, Jacquelyn Noronha-Hostler, Jorge Noronha, Enrico Speranza , et al. (39 additional authors not shown)

Abstract: Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theo… ▽ More Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theory of dense baryonic matter that connects low- and high-energy nuclear physics, astrophysics, gravitational waves physics, and data science △ Less

Submitted 7 November, 2022; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 70 pages, 3 figures, White Paper for the Long Range Plan for Nuclear Science

Report number: LA-UR-22-31648

Modeling Solids in Nuclear Astrophysics with Smoothed Particle Hydrodynamics

Authors: Irina Sagert, Oleg Korobkin, Ingo Tews, Bing-Jyun Tsao, Hyun Lim, Michael J. Falato, Julien Loiseau

Abstract: 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 neut… ▽ More 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. △ Less

Submitted 17 November, 2023; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 28 pages, 34 figures

Report number: LA-UR-22-30036

Journal ref: ApJS 267 47 (2023)

Snowmass 2021 Cosmic Frontier White Paper: The Dense Matter Equation of State and QCD Phase Transitions

Authors: Slavko Bogdanov, Emmanuel Fonseca, Rahul Kashyap, Aleksi Kurkela, James M. Lattimer, Jocelyn S. Read, Bangalore S. Sathyaprakash, H. Thankful Cromartie, Tim Dietrich, Arnab Dhani, Timothy Dolch, Tyler Gorda, Sebastien Guillot, Wynn C. G. Ho, Rachael Huxford, Frederick K. Lamb, Philippe Landry, Bradley W. Meyers, M. Coleman Miller, Joonas Nättilä, Risto Paatelainen, Chanda Prescod-Weinstein, Saga Säppi, Ingrid H. Stairs, Nikolaos Stergioulas , et al. (4 additional authors not shown)

Abstract: Our limited understanding of the physical properties of matter at ultra-high density, high proton/neutron number asymmetry, and low temperature is presently one of the major outstanding problems in physics. As matter in this extreme state is known to only exist stably in the cores of neutron stars (NSs), complementary measurements from electromagnetic and gravitational wave astrophysical observati… ▽ More Our limited understanding of the physical properties of matter at ultra-high density, high proton/neutron number asymmetry, and low temperature is presently one of the major outstanding problems in physics. As matter in this extreme state is known to only exist stably in the cores of neutron stars (NSs), complementary measurements from electromagnetic and gravitational wave astrophysical observations of NSs, combined with terrestrial laboratory constraints and further theoretical investigations, hold the promise to provide important insight into the properties of matter in a region of the quantum chromodynamics phase space that is otherwise inaccessible. This multidisciplinary endeavor imposes the following requirements for facilities and resources in the upcoming decade and beyond: * A next generation of gravitational wave detectors to uncover more double NS and neutron star-black hole mergers; * Sensitive radio telescopes to find the most massive and fastest spinning NSs; * Large-area, high-time-resolution and/or high angular resolution X-ray telescopes to constrain the NS mass-radius relation; * Suitable laboratory facilities for nuclear physics experiments to constrain the dense matter equation of state; * Funding resources for theoretical studies of matter in this regime; * The availability of modern large-scale high performance computing infrastructure. The same facilities and resources would also enable significant advances in other high-profile fields of inquiry in modern physics such as the nature of dark matter, alternative theories of gravity, nucleon superfluidity and superconductivity, as well as an array of astrophysics, including but not limited to stellar evolution, nucleosynthesis, and primordial black holes. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021) under Cosmic Frontier (CF07: Cosmic probes of fundamental physics); 30 pages, 8 figures

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

Authors: Mark G. Alford, Liam Brodie, Alexander Haber, Ingo Tews

Abstract: 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 whic… ▽ More 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. △ Less

Submitted 21 November, 2022; v1 submitted 20 May, 2022; originally announced May 2022.

Comments: V2: Minor corrections, version published in PRC. 12 pages and 5 figures

Report number: LA-UR-22-23972

Journal ref: Phys.Rev.C 106 (2022) 5, 055804

An updated nuclear-physics and multi-messenger astrophysics framework for binary neutron star mergers

Authors: Peter T. H. Pang, Tim Dietrich, Michael W. Coughlin, Mattia Bulla, Ingo Tews, Mouza Almualla, Tyler Barna, Weizmann Kiendrebeogo, Nina Kunert, Gargi Mansingh, Brandon Reed, Niharika Sravan, Andrew Toivonen, Sarah Antier, Robert O. VandenBerg, Jack Heinzel, Vsevolod Nedora, Pouyan Salehi, Ritwik Sharma, Rahul Somasundaram, Chris Van Den Broeck

Abstract: The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic em… ▽ More The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a $1.4M_\odot$ neutron star to be $R=11.98^{+0.35}_{-0.40}$km. △ Less

Submitted 8 January, 2024; v1 submitted 17 May, 2022; originally announced May 2022.

Comments: code available at https://github.com/nuclear-multimessenger-astronomy

Report number: LA-UR-22-23872, LIGO-P2200150

Journal ref: Nature Commun. 14 (2023) 1, 8352

Horizons: Nuclear Astrophysics in the 2020s and Beyond

Authors: H. Schatz, A. D. Becerril Reyes, A. Best, E. F. Brown, K. Chatziioannou, K. A. Chipps, C. M. Deibel, R. Ezzeddine, D. K. Galloway, C. J. Hansen, F. Herwig, A. P. Ji, M. Lugaro, Z. Meisel, D. Norman, J. S. Read, L. F. Roberts, A. Spyrou, I. Tews, F. X. Timmes, C. Travaglio, N. Vassh, C. Abia, P. Adsley, S. Agarwal , et al. (140 additional authors not shown)

Abstract: Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit… ▽ More Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 96 pages. Submitted to Journal of Physics G

Report number: LA-UR-22-23997

Quantifying modeling uncertainties when combining multiple gravitational-wave detections from binary neutron star sources

Authors: Nina Kunert, Peter T. H. Pang, Ingo Tews, Michael W. Coughlin, Tim Dietrich

Abstract: With the increasing sensitivity of gravitational-wave detectors, we expect to observe multiple binary neutron-star systems through gravitational waves in the near future. The combined analysis of these gravitational-wave signals offers the possibility to constrain the neutron-star radius and the equation of state of dense nuclear matter with unprecedented accuracy. However, it is crucial to ensure… ▽ More With the increasing sensitivity of gravitational-wave detectors, we expect to observe multiple binary neutron-star systems through gravitational waves in the near future. The combined analysis of these gravitational-wave signals offers the possibility to constrain the neutron-star radius and the equation of state of dense nuclear matter with unprecedented accuracy. However, it is crucial to ensure that uncertainties inherent in the gravitational-wave models will not lead to systematic biases when information from multiple detections are combined. To quantify waveform systematics, we perform an extensive simulation campaign of binary neutron-star sources and analyse them with a set of four different waveform models. Based on our analysis with about 38 simulations, we find that statistical uncertainties in the neutron-star radius decrease to $\pm 250\rm m$ ($2\%$ at $90\%$ credible interval) but that systematic differences between currently employed waveform models can be twice as large. Hence, it will be essential to ensure that systematic biases will not become dominant in inferences of the neutron-star equation of state when capitalizing on future developments. △ Less

Submitted 10 March, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: 9 pages, 5 figures

Report number: LA-UR-21-27560

Journal ref: Physical Review D (2022), Vol. 105, L061301

Modeling Neutron Star Oscillations in a Fixed General Relativistic Background Including Solid Crust Dynamics

Authors: Bing-Jyun Tsao, Irina Sagert, Oleg Korobkin, Ingo Tews, Hyun Lim, Gary Dilts, Julien Loiseau

Abstract: Measurements of the gravitational-wave signals from neutron star mergers allow scientists to learn about the interior of neutron stars and the properties of dense nuclear matter. The study of neutron star mergers is usually performed with computational fluid dynamics codes, mostly in Eulerian but also in Lagrangian formulation such as smoothed particle hydrodynamics (SPH). Codes include our best k… ▽ More Measurements of the gravitational-wave signals from neutron star mergers allow scientists to learn about the interior of neutron stars and the properties of dense nuclear matter. The study of neutron star mergers is usually performed with computational fluid dynamics codes, mostly in Eulerian but also in Lagrangian formulation such as smoothed particle hydrodynamics (SPH). Codes include our best knowledge of nuclear matter in the form of an equation of state as well as effects of general relativity (GR). However, one important aspect of neutron stars is usually ignored: the solid nature of their crust. The solid matter in the crust is the strongest material known in nature which could lead to a multitude of possible observational effects that have not been studied in dynamical simulations yet. The crust could change the way a neutron star deforms during a merger, leaving an imprint in the gravitational wave signal. It could even shatter during the inspiral, producing a potentially observable electromagnetic signal. Here, we present a first study of the dynamical behavior of neutron stars with a solid crust and fixed GR background with FleCSPH. FleCSPH is a general-purpose SPH code, developed at Los Alamos National Laboratory. It features an efficient algorithm for gravitational interactions via the Fast Multipole Method, which, together with the implemented nuclear equation of state, makes it appropriate for astrophysical applications. The solid material dynamics is described via the elastic-perfectly plastic model with maximum-strain breaking. Despite its simplicity, the model reproduces the stress-strain curve of crustal material as extracted from microphysical simulations very well. We present first tests of our implementation via simulations of neutron star oscillations and give an outlook on our study of the dynamical behavior of the solid crust in neutron star merger events. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 9 pages, 9 figures, proceedings of 2021 international SPHERIC workshop (virtual), June, 8-12 2021

Report number: LA-UR-21-24191

Constraining Neutron-Star Matter with Microscopic and Macroscopic Collisions

Authors: S. Huth, P. T. H. Pang, I. Tews, T. Dietrich, A. Le Fèvre, A. Schwenk, W. Trautmann, K. Agarwal, M. Bulla, M. W. Coughlin, C. Van Den Broeck

Abstract: Interpreting high-energy, astrophysical phenomena, such as supernova explosions or neutron-star collisions, requires a robust understanding of matter at supranuclear densities. However, our knowledge about dense matter explored in the cores of neutron stars remains limited. Fortunately, dense matter is not only probed in astrophysical observations, but also in terrestrial heavy-ion collision exper… ▽ More Interpreting high-energy, astrophysical phenomena, such as supernova explosions or neutron-star collisions, requires a robust understanding of matter at supranuclear densities. However, our knowledge about dense matter explored in the cores of neutron stars remains limited. Fortunately, dense matter is not only probed in astrophysical observations, but also in terrestrial heavy-ion collision experiments. In this work, we use Bayesian inference to combine data from astrophysical multi-messenger observations of neutron stars and from heavy-ion collisions of gold nuclei at relativistic energies with microscopic nuclear theory calculations to improve our understanding of dense matter. We find that the inclusion of heavy-ion collision data indicates an increase in the pressure in dense matter relative to previous analyses, shifting neutron-star radii towards larger values, consistent with recent NICER observations. Our findings show that constraints from heavy-ion collision experiments show a remarkable consistency with multi-messenger observations and provide complementary information on nuclear matter at intermediate densities. This work combines nuclear theory, nuclear experiment, and astrophysical observations, and shows how joint analyses can shed light on the properties of neutron-rich supranuclear matter over the density range probed in neutron stars. △ Less

Submitted 24 June, 2022; v1 submitted 13 July, 2021; originally announced July 2021.

Comments: 7 pages, 2 figures, Supplemental Material, minor changes

Report number: LA-UR-21-22072

Journal ref: Nature 606, 276-280 (2022)

A Detailed Examination of Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of $^{208}$Pb with Minimal Modeling Assumptions

Authors: Reed Essick, Philippe Landry, Achim Schwenk, Ingo Tews

Abstract: The symmetry energy and its density dependence are pivotal for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of $0.283\pm0.071$ fm for the neutron-skin thickness of $^{208}$Pb, $R_{\rm skin}^{^{208}\text{Pb}}$, implying a sy… ▽ More The symmetry energy and its density dependence are pivotal for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of $0.283\pm0.071$ fm for the neutron-skin thickness of $^{208}$Pb, $R_{\rm skin}^{^{208}\text{Pb}}$, implying a symmetry-energy slope parameter $L$ of $106\pm37$ MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy $S_0$, $L$, and $R_{\rm skin}^{^{208}\text{Pb}}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER favor smaller values of the neutron skin and $L$, as well as negative symmetry incompressibilities. Combining astrophysical data with chiral effective field theory ($χ$EFT) and PREX-II constraints yields $S_0 = 33.0^{+2.0}_{-1.8}$ MeV, $L=53^{+14}_{-15}$ MeV, and $R_{\rm skin}^{^{208}\text{Pb}} = 0.17^{+0.04}_{-0.04}$ fm. We also examine the consistency of several individual $χ$EFT calculations with astrophysical observations and terrestrial experiments. We find that there is only mild tension between $χ$EFT, astrophysical data, and PREX-II's $R_\mathrm{skin}^{^{208}\mathrm{Pb}}$ measurement ($p$-value $= 12.3\%$) and that there is excellent agreement between $χ$EFT, astrophysical data, and other nuclear experiments. △ Less

Submitted 15 December, 2021; v1 submitted 12 July, 2021; originally announced July 2021.

Comments: 18 pages, 12 figures, 1 table

Report number: LA-UR-21-21363

Journal ref: Phys. Rev. C 104, 065804 (2021)

Nuclear-Physics Multi-Messenger Astrophysics Constraints on the Neutron-Star Equation of State: Adding NICER's PSR J0740+6620 Measurement

Authors: Peter T. H. Pang, Ingo Tews, Michael W. Coughlin, Mattia Bulla, Chris Van Den Broeck, Tim Dietrich

Abstract: In the past few years, new observations of neutron stars and neutron-star mergers have provided a wealth of data that allow one to constrain the equation of state of nuclear matter at densities above nuclear saturation density. However, most observations were based on neutron stars with masses of about 1.4 solar masses, probing densities up to $\sim$ 3-4 times the nuclear saturation density. Even… ▽ More In the past few years, new observations of neutron stars and neutron-star mergers have provided a wealth of data that allow one to constrain the equation of state of nuclear matter at densities above nuclear saturation density. However, most observations were based on neutron stars with masses of about 1.4 solar masses, probing densities up to $\sim$ 3-4 times the nuclear saturation density. Even higher densities are probed inside massive neutron stars such as PSR J0740+6620. Very recently, new radio observations provided an update to the mass estimate for PSR J0740+6620 and X-ray observations by the NICER and XMM telescopes constrained its radius. Based on these new measurements, we revisit our previous nuclear-physics multi-messenger astrophysics constraints and derive updated constraints on the equation of state describing the neutron-star interior. By combining astrophysical observations of two radio pulsars, two NICER measurements, the two gravitational-wave detections GW170817 and GW190425, detailed modeling of the kilonova AT2017gfo, as well as the gamma-ray burst GRB170817A, we are able to estimate the radius of a typical 1.4-solar mass neutron star to be $11.94^{+0.76}_{-0.87} \rm{km}$ at $90\%$ confidence. Our analysis allows us to revisit the upper bound on the maximum mass of neutron stars and disfavors the presence of a strong first-order phase transition from nuclear matter to exotic forms of matter, such as quark matter, inside neutron stars. △ Less

Submitted 25 October, 2021; v1 submitted 18 May, 2021; originally announced May 2021.

Comments: 13 pages, 5 figures

Report number: LA-UR-21-20534

Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of $^{208}$Pb with Minimal Modeling Assumptions

Authors: Reed Essick, Ingo Tews, Philippe Landry, Achim Schwenk

Abstract: The symmetry energy and its density dependence are crucial inputs for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of $0.283 \pm 0.071$ fm for the neutron-skin thickness of $^{208}$Pb, implying a slope parameter… ▽ More The symmetry energy and its density dependence are crucial inputs for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of $0.283 \pm 0.071$ fm for the neutron-skin thickness of $^{208}$Pb, implying a slope parameter $L = 106 \pm 37$ MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy $S_0$, $L$, and $R_\mathrm{skin}^{^{208}\mathrm{Pb}}$ directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER clearly favor smaller values of the neutron skin and $L$, as well as negative symmetry incompressibilities. Combining astrophysical data with PREX-II and chiral effective field theory constraints yields $S_0 = 33.0^{+2.0}_{-1.8}$ MeV, $L=53^{+14}_{-15}$ MeV, and $R_\mathrm{skin}^{^{208}\mathrm{Pb}}=0.17^{+0.04}_{-0.04}$ fm. △ Less

Submitted 3 November, 2021; v1 submitted 19 February, 2021; originally announced February 2021.

Comments: 11 pages, 6 figures, 2 tables

Report number: LA-UR-21-20527

Combining Electromagnetic and Gravitational-Wave Constraints on Neutron-Star Masses and Radii

Authors: Mohammad Al-Mamun, Andrew W. Steiner, Joonas Nättilä, Jacob Lange, Richard O'Shaughnessy, Ingo Tews, Stefano Gandolfi, Craig Heinke, Sophia Han

Abstract: We perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass X-ray binaries (QLMXBs), photospheric radius expansion X-ray bursts (PREs), and X-ray timing observations of J0030+0451. With this data set, the form of the prior distribution still has an impact on the posterior mass-radius (MR) curves and equation of state (EOS)… ▽ More We perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass X-ray binaries (QLMXBs), photospheric radius expansion X-ray bursts (PREs), and X-ray timing observations of J0030+0451. With this data set, the form of the prior distribution still has an impact on the posterior mass-radius (MR) curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an "intrinsic scattering" contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS. △ Less

Submitted 17 November, 2020; v1 submitted 28 August, 2020; originally announced August 2020.

Comments: 6 pages, 4 figures and supplemental material

Report number: LA-UR-20-26423

Journal ref: Phys. Rev. Lett. 126, 061101 (2021)

On the nature of GW190814 and its impact on the understanding of supranuclear matter

Authors: Ingo Tews, Peter T. H. Pang, Tim Dietrich, Michael W. Coughlin, Sarah Antier, Mattia Bulla, Jack Heinzel, Lina Issa

Abstract: The observation of a compact object with a mass of $2.50-2.67M_{\odot}$ on August 14, 2019, 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… ▽ More The observation of a compact object with a mass of $2.50-2.67M_{\odot}$ on August 14, 2019, 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 $\sim 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 $c_s\geq \sqrt{0.6}$ times the speed of light, while the binary black hole scenario does not offer any new insight. △ Less

Submitted 5 February, 2021; v1 submitted 12 July, 2020; originally announced July 2020.

Comments: 9 pages, 3 figures, published version

Report number: LA-UR-20-24959

Journal ref: Astrophys. J. Lett. 908, L1 (2021)

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

Authors: Peter T. H. Pang, Tim Dietrich, Ingo Tews, Chris Van Den Broeck

Abstract: 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 c… ▽ More 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. △ Less

Submitted 29 September, 2020; v1 submitted 26 June, 2020; originally announced June 2020.

Comments: 17 pages, 11 figures

Report number: LA-UR-20-24019

Journal ref: Phys. Rev. Research 2, 033514 (2020)

Direct Astrophysical Tests of Chiral Effective Field Theory at Supranuclear Densities

Authors: Reed Essick, Ingo Tews, Philippe Landry, Sanjay Reddy, Daniel E. Holz

Abstract: 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 ($χ$EFT) bound our theoretical ignorance. In this work, we analyz… ▽ More 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 ($χ$EFT) bound our theoretical ignorance. In this work, we analyze astrophysical data using a nonparametric representation of the neutron-star EoS conditioned on $χ$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 $χ$EFT at low densities. We also demonstrate how to exploit astrophysical data to directly test the predictions of $χ$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 $χ$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 $χ$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 $χ$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. △ Less

Submitted 8 November, 2020; v1 submitted 16 April, 2020; originally announced April 2020.

Report number: LA-UR-20-22615

Journal ref: Phys. Rev. C 102, 055803 (2020)

Multi-messenger constraints on the neutron-star equation of state and the Hubble constant

Authors: Tim Dietrich, Michael W. Coughlin, Peter T. H. Pang, Mattia Bulla, Jack Heinzel, Lina Issa, Ingo Tews, Sarah Antier

Abstract: Observations of neutron-star mergers based on distinct messengers, including gravitational waves and electromagnetic signals, can be used to study the behavior of matter denser than an atomic nucleus, and to measure the expansion rate of the Universe described by the Hubble constant. We perform a joint analysis of the gravitational-wave signal GW170817 with its electromagnetic counterparts AT2017g… ▽ More Observations of neutron-star mergers based on distinct messengers, including gravitational waves and electromagnetic signals, can be used to study the behavior of matter denser than an atomic nucleus, and to measure the expansion rate of the Universe described by the Hubble constant. We perform a joint analysis of the gravitational-wave signal GW170817 with its electromagnetic counterparts AT2017gfo and GRB170817A, and the gravitational-wave signal GW190425, both originating from neutron-star mergers. We combine these with previous measurements of pulsars using X-ray and radio observations, and nuclear-theory computations using chiral effective field theory to constrain the neutron-star equation of state. We find that the radius of a $1.4$ solar mass neutron star is $11.75^{+0.86}_{-0.81}\ \rm km$ at $90\%$ confidence and the Hubble constant is $66.2^{+4.4}_{-4.2}\ \rm km \,Mpc^{-1}\, s^{-1}$ at $1σ$ uncertainty. △ Less

Submitted 21 December, 2020; v1 submitted 26 February, 2020; originally announced February 2020.

Report number: LA-UR-20-21470

Journal ref: Science (2020), Vol. 370, Issue 6523, pp. 1450-1453

Nuclear and neutron-star matter from local chiral interactions

Authors: D. Lonardoni, I. Tews, S. Gandolfi, J. Carlson

Abstract: We report a quantum Monte Carlo calculation of the equation of state of symmetric nuclear matter using local interactions derived from chiral effective field theory up to next-to-next-to-leading order fit to few-body observables only. The empirical saturation density and energy are well reproduced within statistical and systematic uncertainties. We have also derived the symmetry energy as a functi… ▽ More We report a quantum Monte Carlo calculation of the equation of state of symmetric nuclear matter using local interactions derived from chiral effective field theory up to next-to-next-to-leading order fit to few-body observables only. The empirical saturation density and energy are well reproduced within statistical and systematic uncertainties. We have also derived the symmetry energy as a function of the density, finding good agreement with available experimentally derived constraints at saturation and twice saturation density. We find that the corresponding pressure is also in excellent agreement with recent constraints extracted from gravitational waves of the neutron-star merger GW170817 by the LIGO-Virgo detection. △ Less

Submitted 7 May, 2020; v1 submitted 19 December, 2019; originally announced December 2019.

Comments: 7 pages, 3 figures, 3 tables

Report number: LA-UR-19-32538

Journal ref: Phys. Rev. Research 2, 022033 (2020)

Stringent constraints on neutron-star radii from multimessenger observations and nuclear theory

Authors: Collin D. Capano, Ingo Tews, Stephanie M. Brown, Ben Margalit, Soumi De, Sumit Kumar, Duncan A. Brown, Badri Krishnan, Sanjay Reddy

Abstract: The properties of neutron stars are determined by the nature of the matter that they contain. These properties can be constrained by measurements of the star's size. We obtain stringent constraints on neutron-star radii by combining multimessenger observations of the binary neutron-star merger GW170817 with nuclear theory that best accounts for density-dependent uncertainties in the equation of st… ▽ More The properties of neutron stars are determined by the nature of the matter that they contain. These properties can be constrained by measurements of the star's size. We obtain stringent constraints on neutron-star radii by combining multimessenger observations of the binary neutron-star merger GW170817 with nuclear theory that best accounts for density-dependent uncertainties in the equation of state. We construct equations of state constrained by chiral effective field theory and marginalize over these using the gravitational-wave observations. Combining this with the electromagnetic observations of the merger remnant that imply the presence of a short-lived hyper-massive neutron star, we find that the radius of a $1.4\,\rm{M}_\odot$ neutron star is $R_{1.4\,\mathrm{M}_\odot} = 11.0^{+0.9}_{-0.6}~{\rm km}$ (90% credible interval). Using this constraint, we show that neutron stars are unlikely to be disrupted in neutron-star black-hole mergers; subsequently, such events will not produce observable electromagnetic emission. △ Less

Submitted 24 March, 2020; v1 submitted 27 August, 2019; originally announced August 2019.

Comments: 39 pages, 4 figures

Report number: INT-PUB-19-037, LA-UR-19-28442

Journal ref: Nature Astronomy, 2020

Spin-polarized neutron matter, the maximum mass of neutron stars, and GW170817

Authors: I. Tews, A. Schwenk

Abstract: We investigate how a phase transition from neutron-star matter to spin-polarized neutron matter affects the equation of state and mass-radius relation of neutron stars. While general extension schemes for the equation of state allow for high pressures inside neutron stars, we find that a phase transition to spin-polarized neutron matter excludes extreme regimes. Hence, such a transition limits the… ▽ More We investigate how a phase transition from neutron-star matter to spin-polarized neutron matter affects the equation of state and mass-radius relation of neutron stars. While general extension schemes for the equation of state allow for high pressures inside neutron stars, we find that a phase transition to spin-polarized neutron matter excludes extreme regimes. Hence, such a transition limits the maximum mass of neutron stars to lie below 2.6$-$$2.9 \, M_{\odot}$, depending on the microscopic nuclear forces used, while significantly larger masses could be reached without these constraints. These limits are in good agreement with recent constraints extracted from the neutron-star merger GW170817 and its electromagnetic counterpart. Assuming the description in terms of spin-polarized neutron matter to be valid in the center of neutron stars, we find that stars with a large spin-polarized domain in their core are ruled out by GW170817. △ Less

Submitted 24 March, 2020; v1 submitted 7 August, 2019; originally announced August 2019.

Comments: 10 pages, 3 figures, published version

Report number: LA-UR-19-27912

Journal ref: Astrophys. J. 892, 14 (2020)

Constraining the properties of dense matter and neutron stars by combining nuclear physics and gravitational waves from GW170817

Authors: I. Tews, J. Margueron, S. Reddy

Abstract: Gravitational waves from neutron-star mergers are expected to provide stringent constraints on the structure of neutron stars. At the same time, recent advances in nuclear theory have enabled reliable calculations of the low density equation of state using effective field theory based Hamiltonians and advanced techniques to solve the quantum many-body problem. In this paper, we address how the fir… ▽ More Gravitational waves from neutron-star mergers are expected to provide stringent constraints on the structure of neutron stars. At the same time, recent advances in nuclear theory have enabled reliable calculations of the low density equation of state using effective field theory based Hamiltonians and advanced techniques to solve the quantum many-body problem. In this paper, we address how the first observation of gravitational waves from GW170817 can be combined with modern calculations of the equation of state to extract useful insights about the equation of state of matter encountered inside neutron stars. We analyze the impact of various uncertainties and we show that the tidal deformability extracted from GW170817 is compatible, while less constraining, than modern nuclear physics knowledge. △ Less

Submitted 11 June, 2019; v1 submitted 23 May, 2019; originally announced May 2019.

Comments: arXiv admin note: substantial text overlap with arXiv:1901.09874

Report number: LA-UR-19-24461, INT-PUB-19-017

Journal ref: Xiamen-CUSTIPEN Workshop on the EOS of Dense Neutron-Rich Matter in the Era of Gravitational Wave Astronomy, AIP Conference Proceedings 2019, American Institute of Physics

From the microscopic to the macroscopic world: from nucleons to neutron stars

Authors: S. Gandolfi, J. Lippuner, A. W. Steiner, I. Tews, X. Du, M. Al-Mamun

Abstract: Recent observations of neutron-star properties, in particular the recent detection of gravitational waves emitted from binary neutron stars, GW 170817, open the way to put strong constraints on nuclear interactions. In this paper, we review the state of the art in calculating the equation of state of strongly interacting matter from first principle calculations starting from microscopic interactio… ▽ More Recent observations of neutron-star properties, in particular the recent detection of gravitational waves emitted from binary neutron stars, GW 170817, open the way to put strong constraints on nuclear interactions. In this paper, we review the state of the art in calculating the equation of state of strongly interacting matter from first principle calculations starting from microscopic interactions among nucleons. We then review selected properties of neutron stars that can be directly compared with present and future observations. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: 28 pages, 9 figures, Special Issue Article for Journal of Physics G

Report number: LA-UR-19-22359

Journal ref: J. Phys. G 46 (2019) 103001

Confronting gravitational-wave observations with modern nuclear physics constraints

Authors: I. Tews, J. Margueron, S. Reddy

Abstract: Multi-messenger observations of neutron star (NS) mergers have the potential to revolutionize nuclear astrophysics. They will improve our understanding of nucleosynthesis, provide insights about the equation of state (EOS) of strongly-interacting matter at high densities, and enable tests of the theory of gravity and of dark matter. Here, we focus on the EOS, where both gravitational waves (GWs) f… ▽ More Multi-messenger observations of neutron star (NS) mergers have the potential to revolutionize nuclear astrophysics. They will improve our understanding of nucleosynthesis, provide insights about the equation of state (EOS) of strongly-interacting matter at high densities, and enable tests of the theory of gravity and of dark matter. Here, we focus on the EOS, where both gravitational waves (GWs) from neutron-star mergers and X-ray observations from space-based detectors such as NICER will provide more stringent constraints on the structure of neutron stars. Furthermore, recent advances in nuclear theory have enabled reliable calculations of the EOS at low densities using effective field theory based Hamiltonians and advanced techniques to solve the quantum many-body problem. In this paper, we address how the first observation of GWs from GW170817 can be combined with modern calculations of the EOS to extract useful insights about the EOS of matter encountered inside neutron stars. We analyze the impact of various uncertainties, the role of phase transitions in the NS core, and discuss how future observations will improve our understanding of dense matter. △ Less

Submitted 21 June, 2019; v1 submitted 28 January, 2019; originally announced January 2019.

Comments: 18 pages, 10 figures, published version; Invited contribution prepared for the EPJ A Topical Issue "First joint gravitational wave and electromagnetic observations: Implications for nuclear and particle physics."

Report number: LA-UR-19-20198

Journal ref: Eur. Phys. J. A 55, 97 (2019)

Dense matter with eXTP

Authors: Anna L. Watts, Wenfei Yu, Juri Poutanen, Shu Zhang, Sudip Bhattacharyya, Slavko Bogdanov, Long Ji, Alessandro Patruno, Thomas E. Riley, Pavel Bakala, Altan Baykal, Federico Bernardini, Ignazio Bombaci, Edward Brown, Yuri Cavecchi, Deepto Chakrabarty, Jérôme Chenevez, Nathalie Degenaar, Melania Del Santo, Tiziana Di Salvo, Victor Doroshenko, Maurizio Falanga, Robert D. Ferdman, Marco Feroci, Angelo F. Gambino , et al. (51 additional authors not shown)

Abstract: In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, b… ▽ More In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

A critical examination of constraints on the equation of state of dense matter obtained from GW170817

Authors: I. Tews, J. Margueron, S. Reddy

Abstract: The correlation of the tidal polarizabilities $Λ_1$-$Λ_2$ for GW170817 is predicted by combining dense-matter equations of state (EOSs) that satisfy nuclear physics constraints with the chirp mass and mass asymmetry for this event. Our models are constrained by calculations of the neutron-matter EOS using chiral effective field theory Hamiltonians with reliable error estimates up to once or twice… ▽ More The correlation of the tidal polarizabilities $Λ_1$-$Λ_2$ for GW170817 is predicted by combining dense-matter equations of state (EOSs) that satisfy nuclear physics constraints with the chirp mass and mass asymmetry for this event. Our models are constrained by calculations of the neutron-matter EOS using chiral effective field theory Hamiltonians with reliable error estimates up to once or twice the nuclear saturation density. In the latter case, we find that GW170817 does not improve our understanding of the EOS. We contrast two distinct extrapolations to higher density: a minimal model (MM) which assumes that the EOS is a smooth function of density described by a Taylor expansion and a more general model parametrized by the speed of sound that admits phase transitions. This allows us to identify regions in the $Λ_1$-$Λ_2$ plots that could favor the existence of new phases of matter in neutron stars. We predict the combined tidal polarizability of the two neutron stars in GW170817 to be $80\le \tildeΛ\le 580$ ($280\le \tildeΛ\le 480$ for the MM), which is smaller than the range suggested by the LIGO-Virgo data analysis. Our analysis also shows that GW170817 requires a NS with $M=1.4M_\odot$ to have a radius $9.0<R_{1.4}<13.6$~km ($ 11.3 <R_{1.4}< 13.6$~km for the MM). △ Less

Submitted 24 October, 2018; v1 submitted 8 April, 2018; originally announced April 2018.

Comments: 6 pages, 4 figures, published version

Report number: INT-PUB-18-014

Journal ref: Phys. Rev. C 98, 045804 (2018)

Symmetry Parameter Constraints From A Lower Bound On The Neutron-Matter Energy

Authors: Ingo Tews, James M. Lattimer, Akira Ohnishi, Evgeni E. Kolomeitsev

Abstract: We propose the existence of a lower bound on the energy of pure neutron matter (PNM) on the basis of unitary-gas considerations. We discuss its justification from experimental studies of cold atoms as well as from theoretical studies of neutron matter. We demonstrate that this bound results in limits to the density-dependent symmetry energy, which is the difference between the energies of symmetri… ▽ More We propose the existence of a lower bound on the energy of pure neutron matter (PNM) on the basis of unitary-gas considerations. We discuss its justification from experimental studies of cold atoms as well as from theoretical studies of neutron matter. We demonstrate that this bound results in limits to the density-dependent symmetry energy, which is the difference between the energies of symmetric nuclear matter and PNM. In particular, this bound leads to a lower limit to the volume symmetry energy parameter $S_0$. In addition, for assumed values of $S_0$ above this minimum, this bound implies both upper and lower limits to the symmetry energy slope parameter $L$, which describes the lowest-order density dependence of the symmetry energy. A lower bound on the neutron-matter incompressibility is also obtained. These bounds are found to be consistent with both recent calculations of the energies of PNM and constraints from nuclear experiments. Our results are significant because several equations of state that are currently used in astrophysical simulations of supernovae and neutron star mergers, as well as in nuclear physics simulations of heavy-ion collisions, have symmetry energy parameters that violate these bounds. Furthermore, below the nuclear saturation density, the bound on neutron-matter energies leads to a lower limit to the density-dependent symmetry energy, which leads to upper limits to the nuclear surface symmetry parameter and the neutron-star crust-core boundary. We also obtain a lower limit to the neutron-skin thicknesses of neutron-rich nuclei. Above the nuclear saturation density, the bound on neutron-matter energies also leads to an upper limit to the symmetry energy, with implications for neutron-star cooling via the direct Urca process. △ Less

Submitted 25 October, 2017; v1 submitted 21 November, 2016; originally announced November 2016.

Comments: 27 pages, 13 figures, published version

Report number: YITP-16-121, INT-PUB-16-045

Journal ref: Astrophys.J. 848, 105 (2017)

Spectrum of shear modes in the neutron-star crust: Estimating the nuclear-physics uncertainties

Authors: Ingo Tews

Abstract: I construct a model of the inner crust of neutron stars using interactions from chiral effective field theory (EFT) in order to calculate its equation of state (EOS), shear properties, and the spectrum of crustal shear modes. I systematically study uncertainties associated with the nuclear physics input, the crust composition, and neutron entrainment, and estimate their impact on crustal shear pro… ▽ More I construct a model of the inner crust of neutron stars using interactions from chiral effective field theory (EFT) in order to calculate its equation of state (EOS), shear properties, and the spectrum of crustal shear modes. I systematically study uncertainties associated with the nuclear physics input, the crust composition, and neutron entrainment, and estimate their impact on crustal shear properties and the shear-mode spectrum. I find that the uncertainties originate mainly in two sources: The neutron-matter EOS and neutron entrainment. I compare the spectrum of crustal shear modes to observed frequencies of quasi-periodic oscillations in the afterglow of giant $γ$-ray bursts and find that all of these frequencies could be described within uncertainties, which are, however, at present too sizable to infer neutron-star properties from observations. △ Less

Submitted 18 January, 2017; v1 submitted 24 July, 2016; originally announced July 2016.

Comments: 13 pages, 10 figures, published version

Report number: INT-PUB-16-020

Journal ref: Phys. Rev. C 95, 015803 (2017)

Neutron matter with chiral EFT interactions: Perturbative and first QMC calculations

Authors: I. Tews, T. Krüger, A. Gezerlis, K. Hebeler, A. Schwenk

Abstract: Neutron matter presents a unique system in chiral effective field theory (EFT), because all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We discuss perturbative and first Quantum Monte Carlo (QMC) calculations of neutron matter with chiral EFT interactions and their astrophysical impact for the equation of state and neutron stars. Neutron matter presents a unique system in chiral effective field theory (EFT), because all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We discuss perturbative and first Quantum Monte Carlo (QMC) calculations of neutron matter with chiral EFT interactions and their astrophysical impact for the equation of state and neutron stars. △ Less

Submitted 14 October, 2013; originally announced October 2013.

Comments: 10 pages, 7 figures, proceedings of the International Conference on Nuclear Theory in the Supercomputing Era 2013, Iowa State University, May 13-17, 2013

The chiral condensate in neutron matter

Authors: T. Krüger, I. Tews, B. Friman, K. Hebeler, A. Schwenk

Abstract: We calculate the chiral condensate in neutron matter at zero temperature based on nuclear forces derived within chiral effective field theory. Two-, three- and four-nucleon interactions are included consistently to next-to-next-to-next-to-leading order (N3LO) of the chiral expansion. We find that the interaction contributions lead to a modest increase of the condensate, thus impeding the restorati… ▽ More We calculate the chiral condensate in neutron matter at zero temperature based on nuclear forces derived within chiral effective field theory. Two-, three- and four-nucleon interactions are included consistently to next-to-next-to-next-to-leading order (N3LO) of the chiral expansion. We find that the interaction contributions lead to a modest increase of the condensate, thus impeding the restoration of chiral symmetry in dense matter and making a chiral phase transition in neutron-rich matter unlikely for densities that are not significantly higher than nuclear saturation density. △ Less

Submitted 30 September, 2013; v1 submitted 8 July, 2013; originally announced July 2013.

Comments: published version, 6 pages, 4 figures

Journal ref: Phys. Lett. B 726, 412 (2013)

Neutron matter from chiral effective field theory interactions

Authors: T. Krüger, I. Tews, K. Hebeler, A. Schwenk

Abstract: The neutron-matter equation of state constrains the properties of many physical systems over a wide density range and can be studied systematically using chiral effective field theory (EFT). In chiral EFT, all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We present details and additional results of the first complete N3LO calculation of the neutron… ▽ More The neutron-matter equation of state constrains the properties of many physical systems over a wide density range and can be studied systematically using chiral effective field theory (EFT). In chiral EFT, all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We present details and additional results of the first complete N3LO calculation of the neutron-matter energy, which includes the subleading three-nucleon as well as the leading four-nucleon forces, and provides theoretical uncertainties. In addition, we discuss the impact of our results for astrophysics: for the supernova equation of state, the symmetry energy and its density derivative, and for the structure of neutron stars. Finally, we give a first estimate for the size of the N3LO many-body contributions to the energy of symmetric nuclear matter, which shows that their inclusion will be important in nuclear structure calculations. △ Less

Submitted 12 August, 2013; v1 submitted 8 April, 2013; originally announced April 2013.

Comments: published version; 21 pages, 11 figures, 5 tables

Journal ref: Phys. Rev. C 88, 025802 (2013)

Neutron matter at next-to-next-to-next-to-leading order in chiral effective field theory

Authors: I. Tews, T. Krueger, K. Hebeler, A. Schwenk

Abstract: Neutron matter presents a unique system for chiral effective field theory (EFT), because all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We present the first complete N3LO calculation of the neutron matter energy. This includes the subleading three-nucleon (3N) forces for the first time and all leading four-nucleon (4N) forces. We find relatively… ▽ More Neutron matter presents a unique system for chiral effective field theory (EFT), because all many-body forces among neutrons are predicted to next-to-next-to-next-to-leading order (N3LO). We present the first complete N3LO calculation of the neutron matter energy. This includes the subleading three-nucleon (3N) forces for the first time and all leading four-nucleon (4N) forces. We find relatively large contributions from N3LO 3N forces. Our results provide constraints for neutron-rich matter in astrophysics with controlled theoretical uncertainties. △ Less

Submitted 4 February, 2013; v1 submitted 31 May, 2012; originally announced June 2012.

Comments: 5 pages, 4 figures; improved version, 3N ring and 2pi-contact contributions corrected, conclusions unchanged; v3: minor changes, published version

Journal ref: Phys. Rev. Lett. 110, 032504 (2013)