Search NASA⌕ Search

Engineering topics

Coughlin, Michael W.

Publications and source records attributed to Coughlin, Michael W..

Probing quarkyonic matter in neutron stars with the Bayesian nuclear-physics multimessenger astrophysics framework

The interiors of neutron stars contain 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$ = 1/3) before decreasing again at higher densities. Further, the decrease could be explained by either a strong first-order phase transition or a crossover 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 multimessenger astrophysics (NMMA) framework 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 cores of neutron stars with masses in the two-solar-mass range, amounting to up to ≈0.13$M$ ⊙ , contributing ≈ 5.9% of the total mass. Furthermore, for the quarkyonic matter model investigated here, the radius of a 1.4$M$ ⊙ neutron star would be $13.44_{–1.54}^{+1.69}(13. 54_{–1.04}^{+1.02})$ km, at 95% credibility, without (with) the inclusion of AT2017gfo.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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

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 ⊙ neutron star to be $R$ $=$ $11.98^{+0.35}_{-0.04}$ km.

79 ASTRONOMY AND ASTROPHYSICS↗

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

Here, 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 is combined. To quantify waveform systematics, we perform an extensive simulation campaign of binary neutron-star sources and analyze them with a set of four different waveform models. For our analysis with 38 simulations, we find that statistical uncertainties in the neutron-star radius decrease to ± 250 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.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Challenges Ahead for Multimessenger Analyses of Gravitational Waves and Kilonova: A Case Study on GW190425

In recent years, there have been significant advances in multimessenger astronomy due to the discovery of the first, and so far only confirmed, gravitational wave event with a simultaneous electromagnetic (EM) counterpart, as well as improvements in numerical simulations, gravitational wave (GW) detectors, and transient astronomy. This has led to the exciting possibility of performing joint analyses of the GW and EM data, providing additional constraints on fundamental properties of the binary progenitor and merger remnant. Here, in this work, we present a new Bayesian framework that allows inference of these properties, while taking into account the systematic modeling uncertainties that arise when mapping from GW binary progenitor properties to photometric light curves. We extend the relative binning method presented in Zackay et al. to include extrinsic GW parameters for fast analysis of the GW signal. The focus of our EM framework is on light curves arising from r-process nucleosynthesis in the ejected material during and after merger, the so-called kilonova, and particularly on black hole–neutron star systems. As a case study, we examine the recent detection of GW190425, where the primary object is consistent with being either a black hole or a neutron star. We show quantitatively how improved mapping between binary progenitor and outflow properties, and/or an increase in EM data quantity and quality are required in order to break degeneracies in the fundamental source parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

The large-scale environment of thermonuclear and core-collapse supernovae

ABSTRACT The new generation of wide-field time-domain surveys has made it feasible to study the clustering of supernova (SN) host galaxies in the large-scale structure (LSS) for the first time. We investigate the LSS environment of SN populations, using 106 dark matter density realisations with a resolution of ∼3.8 Mpc, constrained by the 2M+ + galaxy survey. We limit our analysis to redshift z < 0.036, using samples of 498 thermonuclear and 782 core-collapse SNe from the Zwicky Transient Facility’s Bright Transient Survey and Census of the Local Universe catalogues. We detect clustering of SNe with high significance; the observed clustering of the two SNe populations is consistent with each other. Further, the clustering of SN hosts is consistent with that of the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey DR12 spectroscopic galaxy sample in the same redshift range. Using a tidal shear classifier, we classify the LSS into voids, sheets, filaments, and knots. We find that both SNe and SDSS galaxies are predominantly found in sheets and filaments. SNe are significantly under-represented in voids and over-represented in knots compared to the volume fraction in these structures. This work opens the potential for using forthcoming wide-field deep SN surveys as a complementary LSS probe.

Astronomy & Astrophysics↗

Nuclear Physics Multimessenger Astrophysics Constraints on the Neutron Star Equation of State: Adding NICER’s PSR J0740+6620 Measurement

In the past few years, new observations of neutron stars (NSs) and NS mergers have provided a wealth of data that allow one to constrain the equation of state (EOS) of nuclear matter at densities above nuclear saturation density. However, most observations were based on NSs with masses of about 1.4 M ⊙ , probing densities up to ~three to four times the nuclear saturation density. Even higher densities are probed inside massive NSs 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 multimessenger astrophysics constraints and derive updated constraints on the EOS describing the NS 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 AT 2017gfo, and the gamma-ray burst GRB 170817A, we are able to estimate the radius of a typical 1.4 M ⊙ NS to be ${11.94}_{-0.87}^{+0.76}\,\mathrm{km}$ at 90% confidence. Finally, our analysis allows us to revisit the upper bound on the maximum mass of NSs and disfavors the presence of a strong first-order phase transition from nuclear matter to exotic forms of matter, such as quark matter, inside NSs.

60 APPLIED LIFE SCIENCES↗

2900 Square Degree Search for the Optical Counterpart of Short Gamma-Ray Burst GRB 180523B with the Zwicky Transient Facility

There is significant interest in the models for production of short gamma-ray bursts (GRBs). Until now, the number of known short GRBs with multi-wavelength afterglows has been small. While the Fermi GRB Monitor detects many GRBs relative to the Neil Gehrels Swift Observatory, the large localization regions makes the search for counterparts difficult. With the Zwicky Transient Facility (ZTF - part of Palomar Observatory) recently achieving first light, it is now fruitful to use its combination of depth (m (sub AB) approximating 20.6), field of view (approximately 47 square degrees), and survey cadence (every approximately 3 days) to perform Target of Opportunity observations. We demonstrate this capability on GRB 180523B, which was recently announced by the Fermi GRB Monitor as a short GRB. ZTF imaged 2900 square degrees of the localization region, resulting in the coverage of 61.6 percent of the enclosed probability over two nights to a depth of m (sub AB) approximating 20.5. We characterized 14 previously unidentified transients, and none were found to be consistent with a short GRB counterpart. This search with the ZTF shows it is an efficient camera for searching for coarsely localized short GRB and gravitational-wave counterparts, allowing for a sensitive search with minimal interruption to its nominal cadence.

Coughlin, Michael W.↗