-
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.
-
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.
-
Numerical relativity simulations of prompt collapse mergers: threshold mass and phenomenological constraints on neutron star properties after GW170817
Authors:
Rahul Kashyap,
Abhishek Das,
David Radice,
Surendra Padamata,
Aviral Prakash,
Domenico Logoteta,
Albino Perego,
Daniel A. Godzieba,
Sebastiano Bernuzzi,
Ignazio Bombaci,
Farrukh J. Fattoyev,
Brendan T. Reed,
André da Silva Schneider
Abstract:
We determine the threshold mass for prompt (no bounce) black hole formation in equal-mass neutron star (NS) mergers using a new set of 227 numerical relativity simulations. We consider 23 phenomenological and microphysical finite temperature equations of state (EOS), including models with hyperons and first-order phase transitions to deconfined quarks. We confirm the existence of EOS-insensitive r…
▽ More
We determine the threshold mass for prompt (no bounce) black hole formation in equal-mass neutron star (NS) mergers using a new set of 227 numerical relativity simulations. We consider 23 phenomenological and microphysical finite temperature equations of state (EOS), including models with hyperons and first-order phase transitions to deconfined quarks. We confirm the existence of EOS-insensitive relations between the threshold mass, the binary tidal parameter at the threshold ($Λ_{th}$), the maximum mass of nonrotating NSs, and the radii of reference mass NSs. We correct the systematic errors in previously reported fitting coefficients that were obtained with approximate general-relativity simulations. We combine the EOS-insensitive relations, phenomenological constraints on NS properties and observational data from GW170817 to derive an improved lower limit on radii of maximum mass and 1.6 M$_\odot$ NS of 9.81 km and 10.90 km, respectively. We also constrain the radius and quadrupolar tidal deformability ($Λ$) of a 1.4 $M_\odot$ NS to be larger than 10.74 km and 172, respectively. We consider uncertainties in all independent parameters -- fitting coefficients as well as GW170817 masses while reporting the range of radii constraints. We introduce new methods to constrain the upper as well as lower limit of NS maximum mass using future BNS detections and their identification as prompt or delayed collapse. With future observations it will be possible to derive even tighter constraints on the properties of matter at and above nuclear density using the method proposed in this work.
△ Less
Submitted 8 April, 2022; v1 submitted 9 November, 2021;
originally announced November 2021.
-
Modeling the Galactic Neutron Star Population for Use in Continuous Gravitational Wave Searches
Authors:
Brendan T. Reed,
Alex Deibel,
C. J. Horowitz
Abstract:
Searches for continuous gravitational waves from \textit{unknown} Galactic neutron stars provide limits on the shapes of neutron stars. A rotating neutron star will produce gravitational waves if asymmetric deformations exist in its structure that are characterized by the star's ellipticity. In this study, we use a simple model of the spatial and spin distribution of Galactic neutron stars to esti…
▽ More
Searches for continuous gravitational waves from \textit{unknown} Galactic neutron stars provide limits on the shapes of neutron stars. A rotating neutron star will produce gravitational waves if asymmetric deformations exist in its structure that are characterized by the star's ellipticity. In this study, we use a simple model of the spatial and spin distribution of Galactic neutron stars to estimate the total number of neutron stars probed, using gravitational waves, to a given upper limit on the ellipticity. This may help optimize future searches with improved sensitivity. The improved sensitivity of third-generation gravitational wave detectors may increase the number of neutron stars probed, to a given ellipticity, by factors of 100 to 1000.
△ Less
Submitted 9 August, 2021; v1 submitted 1 April, 2021;
originally announced April 2021.
-
Implications of PREX-II on the equation of state of neutron-rich matter
Authors:
Brendan T. Reed,
F. J. Fattoyev,
C. J. Horowitz,
J. Piekarewicz
Abstract:
Laboratory experiments sensitive to the equation of state of neutron rich matter in the vicinity of nuclear saturation density provide the first rung in a "density ladder" that connects terrestrial experiments to astronomical observations. In this context, the neutron skin thickness of 208Pb (Rskin) provides a stringent laboratory constraint on the density dependence of the symmetry energy. In tur…
▽ More
Laboratory experiments sensitive to the equation of state of neutron rich matter in the vicinity of nuclear saturation density provide the first rung in a "density ladder" that connects terrestrial experiments to astronomical observations. In this context, the neutron skin thickness of 208Pb (Rskin) provides a stringent laboratory constraint on the density dependence of the symmetry energy. In turn, an improved value of Rskin has been reported recently by the PREX collaboration. Exploiting the strong correlation between Rskin and the slope of the symmetry energy L within a specific class of relativistic energy density functionals, we report a value of L=(106 +/- 37)MeV -- that systematically overestimates current limits based on both theoretical approaches and experimental measurements. The impact of such a stiff symmetry energy on some critical neutron-star observables is also examined.
△ Less
Submitted 31 March, 2021; v1 submitted 8 January, 2021;
originally announced January 2021.