Showing 1–50 of 86 results for author: Obergaulinger, M

Nucleosynthesis in magnetorotational supernovae: impact of the magnetic field configuration

Authors: M. Reichert, M. Bugli, J. Guilet, M. Obergaulinger, M. Á. Aloy, A. Arcones

Abstract: The production of heavy elements is one of the main by-products of the explosive end of massive stars. A long sought goal is finding differentiated patterns in the nucleosynthesis yields, which could permit identifying a number of properties of the explosive core. Among them, the traces of the magnetic field topology are particularly important for \emph{extreme} supernova explosions, most likely h… ▽ More The production of heavy elements is one of the main by-products of the explosive end of massive stars. A long sought goal is finding differentiated patterns in the nucleosynthesis yields, which could permit identifying a number of properties of the explosive core. Among them, the traces of the magnetic field topology are particularly important for \emph{extreme} supernova explosions, most likely hosted by magnetorotational effects. We investigate the nucleosynthesis of five state-of-the-art magnetohydrodynamic models with fast rotation that have been previously calculated in full 3D and that involve an accurate neutrino transport (M1). One of the models does not contain any magnetic field and synthesizes elements around the iron group, in agreement with other CC-SNe models in literature. All other models host a strong magnetic field of the same intensity, but with different topology. For the first time, we investigate the nucleosynthesis of MR-SNe models with a quadrupolar magnetic field and a 90 degree tilted dipole. We obtain a large variety of ejecta compositions reaching from iron nuclei to nuclei up to the third r-process peak. We assess the robustness of our results by considering the impact of different nuclear physics uncertainties such as different nuclear masses, $β^{-}$-decays and $β^{-}$-delayed neutron emission probabilities, neutrino reactions, fission, and a feedback of nuclear energy on the temperature. We find that the qualitative results do not change with different nuclear physics input. The properties of the explosion dynamics and the magnetic field configuration are the dominant factors determining the ejecta composition. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Determining the core-collapse supernova explosion mechanism with current and future gravitational-wave observatories

Authors: Jade Powell, Alberto Iess, Miquel Llorens-Monteagudo, Martin Obergaulinger, Bernhard Müller, Alejandro Torres-Forné, Elena Cuoco, José A. Font

Abstract: Gravitational waves are emitted from deep within a core-collapse supernova, which may enable us to determine the mechanism of the explosion from a gravitational-wave detection. Previous studies suggested that it is possible to determine if the explosion mechanism is neutrino-driven or magneto-rotationally powered from the gravitational-wave signal. However, long duration magneto-rotational wavefor… ▽ More Gravitational waves are emitted from deep within a core-collapse supernova, which may enable us to determine the mechanism of the explosion from a gravitational-wave detection. Previous studies suggested that it is possible to determine if the explosion mechanism is neutrino-driven or magneto-rotationally powered from the gravitational-wave signal. However, long duration magneto-rotational waveforms, that cover the full explosion phase, were not available during the time of previous studies, and explosions were just assumed to be magneto-rotationally driven if the model was rapidly rotating. Therefore, we perform an updated study using new 3D long-duration magneto-rotational core-collapse supernova waveforms that cover the full explosion phase, injected into noise for the Advanced LIGO, Einstein Telescope and NEMO gravitational-wave detectors. We also include a category for failed explosions in our signal classification results. We then determine the explosion mechanism of the signals using three different methods: Bayesian model selection, dictionary learning, and convolutional neural networks. The three different methods are able to distinguish between neutrino-driven explosions and magneto-rotational explosions, even if the neutrino-driven explosion model is rapidly rotating. However they can only distinguish between the non-exploding and neutrino-driven explosions for signals with a high signal to noise ratio. △ Less

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

Comments: Accepted in Physical Review D

Assessment of a new sub-grid model for magnetohydrodynamical turbulence. II. Kelvin-Helmholtz instability

Authors: Miquel Miravet-Tenés, Pablo Cerdá-Durán, Martin Obergaulinger, José A. Font

Abstract: The modelling of astrophysical systems such as binary neutron star mergers or the formation of magnetars from the collapse of massive stars involves the numerical evolution of magnetised fluids at extremely large Reynolds numbers. This is a major challenge for (unresolved) direct numerical simulations which may struggle to resolve highly dynamical features as, e.g. turbulence, magnetic field ampli… ▽ More The modelling of astrophysical systems such as binary neutron star mergers or the formation of magnetars from the collapse of massive stars involves the numerical evolution of magnetised fluids at extremely large Reynolds numbers. This is a major challenge for (unresolved) direct numerical simulations which may struggle to resolve highly dynamical features as, e.g. turbulence, magnetic field amplification, or the transport of angular momentum. Sub-grid models offer a means to overcome those difficulties. In a recent paper we presented MInIT, an MHD-instability-induced-turbulence mean-field, sub-grid model based on the modelling of the turbulent (Maxwell, Reynolds, and Faraday) stress tensors. While in our previous work MInIT was assessed within the framework of the magnetorotational instability, in this paper we further evaluate the model in the context of the Kelvin-Helmholtz instability (KHI). The main difference with other sub-grid models (as e.g. the alpha-viscosity model or the gradient model) is that in MInIT we track independently the turbulent energy density at sub-grid scales, which is used, via a simple closure relation, to compute the different turbulent stresses relevant for the dynamics. The free coefficients of the model are calibrated using well resolved box simulations of magnetic turbulence generated by the KHI. We test the model against these simulations and show that it yields order-of-magnitude accurate predictions for the evolution of the turbulent Reynolds and Maxwell stresses. △ Less

Submitted 19 October, 2023; v1 submitted 11 August, 2023; originally announced August 2023.

Comments: 12 pages, 13 figures Accepted for publication in MNRAS

An Optically Targeted Search for Gravitational Waves emitted by Core-Collapse Supernovae during the Third Observing Run of Advanced LIGO and Advanced Virgo

Authors: Marek J. Szczepańczyk, Yanyan Zheng, Javier M. Antelis, Michael Benjamin, Marie-Anne Bizouard, Alejandro Casallas-Lagos, Pablo Cerdá-Durán, Derek Davis, Dorota Gondek-Rosińska, Sergey Klimenko, Claudia Moreno, Martin Obergaulinger, Jade Powell, Dymetris Ramirez, Brad Ratto, Colter Richarson, Abhinav Rijal, Amber L. Stuver, Paweł Szewczyk, Gabriele Vedovato, Michele Zanolin, Imre Bartos, Shubhagata Bhaumik, Tomasz Bulik, Marco Drago , et al. (13 additional authors not shown)

Abstract: We present the results from a search for gravitational-wave transients associated with core-collapse supernovae observed optically within 30 Mpc during the third observing run of Advanced LIGO and Advanced Virgo. No gravitational wave associated with a core-collapse supernova has been identified. We then report the detection efficiency for a variety of possible gravitational-wave emissions. For ne… ▽ More We present the results from a search for gravitational-wave transients associated with core-collapse supernovae observed optically within 30 Mpc during the third observing run of Advanced LIGO and Advanced Virgo. No gravitational wave associated with a core-collapse supernova has been identified. We then report the detection efficiency for a variety of possible gravitational-wave emissions. For neutrino-driven explosions, the distance at which we reach 50% detection efficiency is up to 8.9 kpc, while more energetic magnetorotationally-driven explosions are detectable at larger distances. The distance reaches for selected models of the black hole formation, and quantum chromodynamics phase transition are also provided. We then constrain the core-collapse supernova engine across a wide frequency range from 50 Hz to 2 kHz. The upper limits on gravitational-wave energy and luminosity emission are at low frequencies down to $10^{-4}\,M_\odot c^2$ and $6 \times 10^{-4}\,M_\odot c^2$/s, respectively. The upper limits on the proto-neutron star ellipticity are down to 3 at high frequencies. Finally, by combining the results obtained with the data from the first and second observing runs of LIGO and Virgo, we improve the constraints of the parameter spaces of the extreme emission models. Specifically, the proto-neutron star ellipticities for the long-lasting bar mode model are down to 1 for long emission (1 s) at high frequency. △ Less

Submitted 3 July, 2024; v1 submitted 25 May, 2023; originally announced May 2023.

Comments: 20 pages, 11 figures; https://dcc.ligo.org/LIGO-P2200361

Nucleosynthesis in jet-driven and jet-associated supernovae

Authors: M. Obergaulinger, M. Reichert

Abstract: In contrast to regular core-collapse supernovae, explosions of rapidly rotating massive stars can develop jets, fast collimated outflows directed along the rotational axis. Depending on the rate of rotation and the magnetic field strength before collapse as well as on possible mechanisms amplifying the magnetic field, such a core can explode magnetorotationally rather than via the standard superno… ▽ More In contrast to regular core-collapse supernovae, explosions of rapidly rotating massive stars can develop jets, fast collimated outflows directed along the rotational axis. Depending on the rate of rotation and the magnetic field strength before collapse as well as on possible mechanisms amplifying the magnetic field, such a core can explode magnetorotationally rather than via the standard supernova mechanism based on neutrino heating. This scenario can explain the highest kinetic energies observed in the class of hypernovae. On longer time scales, rotation and magnetic fields can play an important role in the engine of long gamma-ray burst powered by proto-magnetars or hyperaccreting black holes in collapsars. Both classes of events are characterized by relativistic jets and winds driven by neutrinos or magnetic spin-down of the central objects. The nucleosynthesis in these events includes the production of Fe group elements, including a possibly enhanced synthesis of radioactive 56Ni leading to high peak luminosities. Additionally, these events are, out of all stellar core-collapse events the ones most likely to allow for the formation of the heaviest nuclei via rapid neutron captures. Increasingly sophisticated numerical simulations indicate that at least a limited r-process is possible, though it remains open how robust this result is against variations in the numerical methods and the initial conditions. If so, supernovae with jets could contribute to the observed galactic chemical enrichment, in particular at early times before neutron-star mergers might be able to set in. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: 41 pages, 7 figures. Chapter to appear in Handbook of Nuclear Physics (Springer)

Minimally implicit methods for the numerical integration of the neutrino transport equations

Authors: Samuel Santos-Pérez, Martin Obergaulinger, Isabel Cordero-Carrión

Abstract: The evolution of many astrophysical systems is dominated by the interaction between matter and radiation such as photons or neutrinos. The dynamics can be described by the evolution equations of radiation hydrodynamics in which reactions between matter particles and radiation quanta couples the hydrodynamic equations to those of radiative transfer, see Munier & Weaver (1986a), Munier & Weaver (198… ▽ More The evolution of many astrophysical systems is dominated by the interaction between matter and radiation such as photons or neutrinos. The dynamics can be described by the evolution equations of radiation hydrodynamics in which reactions between matter particles and radiation quanta couples the hydrodynamic equations to those of radiative transfer, see Munier & Weaver (1986a), Munier & Weaver (1986b).. The numerical treatment has to account for their potential stiffness (e.g., in optically thick environments). In this article, we will present a new method to numerically integrate these equations in a stable way by using minimally implicit Runge-Kutta methods. With these methods, the inversion of the implicit operator can be done analytically. We also take into account the physical behavior of the evolved variables in the limit of the stiff regime. We will show the results of applying this method to the reactions between neutrinos and matter in core-collapse supernovae simulations. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages. 5 figures. Submitted article to Astrophysical Journal Series

Inference of proto-neutron star properties in core-collapse supernovae from a gravitational-wave detector network

Authors: T. Bruel, M-A. Bizouard, M. Obergaulinger, P. Maturana-Russel, A. Torres-Forné, P. Cerdá-Durán, N. Christensen, J. A. Font, R. Meyer

Abstract: The next Galactic core-collapse supernova (CCSN) will be a unique opportunity to study within a fully multi-messenger approach the explosion mechanism responsible for the formation of neutron stars and stellar-mass black holes. State-of-the-art numerical simulations of those events reveal the complexity of the gravitational-wave emission which is highly stochastic. This challenges the possibility… ▽ More The next Galactic core-collapse supernova (CCSN) will be a unique opportunity to study within a fully multi-messenger approach the explosion mechanism responsible for the formation of neutron stars and stellar-mass black holes. State-of-the-art numerical simulations of those events reveal the complexity of the gravitational-wave emission which is highly stochastic. This challenges the possibility to infer the properties of the compact remnant and of its progenitor using the information encoded in the waveforms. In this paper we take further steps in a program we recently initiated to overcome those difficulties. In particular we show how oscillation modes of the proto-neutron star, highly visible in the gravitational-wave signal, can be used to reconstruct the time evolution of their physical properties. Extending our previous work where only the information from a single detector was used we here describe a new data-analysis pipeline that coherently combines gravitational-wave detectors' data and infers the time evolution of a combination of the mass and radius of the compact remnant. The performance of the method is estimated employing waveforms from 2D and 3D CCSN simulations covering a progenitor mass range between 11$\mathrm{M_{\odot}}$\, and 40$\mathrm{M_{\odot}}$\, and different equations of state for both a network of up to five second-generation detectors and the proposed third-generation detectors Einstein Telescope and Cosmic Explorer. Our study shows that it will be possible to infer PNS properties for CCSN events occurring in the vicinity of the Milky Way, up to the Large Magellanic Cloud, with the current generation of gravitational-wave detectors. △ Less

Submitted 6 July, 2023; v1 submitted 24 January, 2023; originally announced January 2023.

Journal ref: Phys. Rev. D 107, 083029 (2023)

Core-collapse supernova simulations with reduced nucleosynthesis networks

Authors: Gerard Navó, Moritz Reichert, Martin Obergaulinger, Almudena Arcones

Abstract: We present core-collapse supernova simulations including nuclear reaction networks that impact explosion dynamics and nucleosynthesis. The different composition treatment can lead to changes in the neutrino heating in the vicinity of the shock by modifying the number of nucleons and thus the neutrino-opacity of the region. This reduces the ram pressure outside the shock and allows an easier expans… ▽ More We present core-collapse supernova simulations including nuclear reaction networks that impact explosion dynamics and nucleosynthesis. The different composition treatment can lead to changes in the neutrino heating in the vicinity of the shock by modifying the number of nucleons and thus the neutrino-opacity of the region. This reduces the ram pressure outside the shock and allows an easier expansion. The energy released by the nuclear reactions during collapse also slows down the accretion and aids the shock expansion. In addition, nuclear energy generation in the postshocked matter produces up to $20\%$ more energetic explosions. Nucleosynthesis is affected due to the different dynamic evolution of the explosion. Our results indicate that the energy generation from nuclear reactions helps to sustain late outflows from the vicinity of the proto-neutron star, synthesizing more neutron-rich species. Furthermore, we show that there are systematic discrepancies between the ejecta calculated with in-situ and ex-situ reaction networks. These differences stem from the intrinsic characteristics of evolving the composition in hydrodynamic simulations or calculating it with Lagrangian tracer particles. The mass fractions of some Ca, Ti, Cr, and Fe isotopes are consistently underproduced in postprocessing calculations, leading to different nucleosynthesis paths. Our results suggest that large in-situ nuclear reaction networks are important for a realistic feedback of the energy generation, the neutrino heating, and a more accurate ejecta composition. △ Less

Submitted 13 June, 2023; v1 submitted 21 October, 2022; originally announced October 2022.

Comments: Submitted to ApJ. Received 2022 October 20; revised 2023 May 13; accepted 2023 May 15

Three-dimensional core-collapse supernovae with complex magnetic structures: II. Rotational instabilities and multi-messenger signatures

Authors: Matteo Bugli, Jérôme Guilet, Thierry Foglizzo, Martin Obergaulinger

Abstract: The gravitational collapse of rapidly rotating massive stars can lead to the onset of the low $T/\|W\|$ instability within the central proto-neutron star (PNS), which leaves strong signatures in both the gravitational wave (GW) and neutrino emission. Strong large-scale magnetic fields are usually invoked to explain outstanding stellar explosions of rapidly rotating progenitors, but their impact on… ▽ More The gravitational collapse of rapidly rotating massive stars can lead to the onset of the low $T/\|W\|$ instability within the central proto-neutron star (PNS), which leaves strong signatures in both the gravitational wave (GW) and neutrino emission. Strong large-scale magnetic fields are usually invoked to explain outstanding stellar explosions of rapidly rotating progenitors, but their impact on the growth of such instability has not yet been cleared. We analyze a series of three-dimensional magnetohydrodynamic models to characterize the effects of different magnetic configurations on the development of the low $T/\|W\|$ and the related multi-messenger features. In the absence of magnetic fields, we observe the growth on dynamical time scales of the low $T/\|W\|$, associated with a strong burst of GW and a correlated modulation of the neutrino emission. However, models with a strong magnetic field show a quenching of the low $T/\|W\|$, due to a flattening of the rotation profile in the first $\sim100$ ms after shock formation caused by the magnetic transport of angular momentum. The associated GW emission is weakened by an order of magnitude, exhibits a broader spectral shape, and has no dominant feature associated with the PNS large-scale oscillation modes. Neutrino luminosities are damped along the equatorial plane due to a more oblate PNS, and the only clear modulation in the signal is due to SASI activity. Finally, magnetized models produce lower luminosities for $ν_e$ than for $\barν_e$, which is connected to a higher concentration of neutron-rich material in the PNS surroundings. △ Less

Submitted 9 February, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 14 pages, 17 figures; accepted for publication in MNRAS

Assessment of a new sub-grid model for magneto-hydrodynamical turbulence. I. Magnetorotational instability

Authors: Miquel Miravet-Tenés, Pablo Cerdá-Durán, Martin Obergaulinger, José A. Font

Abstract: Insufficient numerical resolution of grid-based, direct numerical simulations (DNS) hampers the development of instabilitydriven turbulence at small (unresolved) scales. As an alternative to DNS, sub-grid models can potentially reproduce the effects of turbulence at small scales in terms of the resolved scales, and hence can capture physical effects with less computational resources. We present a… ▽ More Insufficient numerical resolution of grid-based, direct numerical simulations (DNS) hampers the development of instabilitydriven turbulence at small (unresolved) scales. As an alternative to DNS, sub-grid models can potentially reproduce the effects of turbulence at small scales in terms of the resolved scales, and hence can capture physical effects with less computational resources. We present a new sub-grid model, the MHD-instability-induced-turbulence (MInIT) mean-field model. MInIT is a physically motivated model based on the evolution of the turbulent (Maxwell, Reynolds, and Faraday) stress tensors and their relation with the turbulent energy densities of the magneto-rotational (MRI) and parasitic instabilities, modeled with two partial differential evolution equations with stiff source terms. Their solution allows obtaining the turbulent stress tensors through the constant coefficients that link them to the energy densities. The model is assessed using data from MRI in-box DNS and applying a filtering operation to compare the filtered data with that from the model. Using the $L_2$-norm as the metric for the comparison, we find less than one order-of-magnitude difference between the two sets of data. No dependence on filter size or length scale of unresolved scales is found, as opposed to results using the gradient model (which we also use to contrast our model) in which the $L_2$-norm of some of the stresses increases with filter size. We conclude that MInIT can help DNS by properly capturing small-scale turbulent stresses which has potential implications on the dynamics of highly-magnetized rotating compact objects, such as those formed during binary neutron star mergers. △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: 20 pages, 15 figures

Magnetorotational supernovae: A nucleosynthetic analysis of sophisticated 3D models

Authors: Moritz Reichert, Martin Obergaulinger, Miguel-Angel Aloy, Michael Gabler, Almudena Arcones, Friedrich-Karl Thielemann

Abstract: Magnetorotational supernovae are a rare type of core-collapse supernovae where the magnetic field and rotation play a central role in the dynamics of the explosion. We present the post-processed nucleosynthesis of state-of-the-art neutrino-MHD supernova models that follow the post explosion evolution for few seconds. We find three different dynamical mechanisms to produce heavy r-process elements:… ▽ More Magnetorotational supernovae are a rare type of core-collapse supernovae where the magnetic field and rotation play a central role in the dynamics of the explosion. We present the post-processed nucleosynthesis of state-of-the-art neutrino-MHD supernova models that follow the post explosion evolution for few seconds. We find three different dynamical mechanisms to produce heavy r-process elements: i) a prompt ejection of matter right after core bounce, ii) neutron-rich matter that is ejected at late times due to a reconfiguration of the protoneutronstar shape, iii) small amount of mass ejected with high entropies in the center of the jet. We investigate total ejecta yields, including the ones of unstable nuclei such as $^{26}$Al, $^{44}$Ti, $^{56}$Ni, and $^{60}$Fe. The obtained $^{56}$Ni masses vary between $0.01 - 1\,\mathrm{M_\odot}$. The latter maximum is compatible with hypernova observations. Furthermore, all of our models synthesize Zn masses in agreement with observations of old metal-poor stars. We calculate simplified light curves to investigate whether our models can be candidates for superluminous supernovae. The peak luminosities obtained from taking into account only nuclear heating reach up to a few $\sim 10^{43} \,\mathrm{erg\,s^{-1}}$. Under certain conditions, we find a significant impact of the $^{66}$Ni decay chain that can raise the peak luminosity up to $\sim 38\%$ compared to models including only the $^{56}$Ni decay chain. This work reinforces the theoretical evidence on the critical role of magnetorotational supernovae to understand the occurrence of hypernovae, superluminous supernovae, and the synthesis of heavy elements. △ Less

Submitted 16 November, 2022; v1 submitted 23 June, 2022; originally announced June 2022.

R-process viable outflows are suppressed in global alpha-viscosity models of collapsar disks

Authors: Oliver Just, Miguel A. Aloy, Martin Obergaulinger, Shigehiro Nagataki

Abstract: Collapsar disks have been proposed to be rich factories of heavy elements, but the major question of whether their outflows are neutron-rich, and could therefore represent significant sites of the rapid neutron-capture (r-) process, or dominated by iron-group elements remains unresolved. We present the first global models of collapsars that start from a stellar progenitor and self-consistently des… ▽ More Collapsar disks have been proposed to be rich factories of heavy elements, but the major question of whether their outflows are neutron-rich, and could therefore represent significant sites of the rapid neutron-capture (r-) process, or dominated by iron-group elements remains unresolved. We present the first global models of collapsars that start from a stellar progenitor and self-consistently describe the evolution of the disk, its composition, and its outflows in response to the imploding stellar mantle, using energy-dependent M1 neutrino transport and an alpha-viscosity to approximate turbulent angular-momentum transport. We find that a neutron-rich, neutrino-dominated accretion flow (NDAF) is established only marginally--either for short times or relatively low viscosities--because the disk tends to disintegrate into an advective disk (ADAF) already at relatively high mass-accretion rates, launching powerful outflows but preventing it from developing a hot, dense, and therefore neutron-rich core. Viscous outflows disrupt the star within ~100s with explosion energies close to that of hypernovae. If viscosity is neglected, a stable NDAF with disk mass of about 1Msun is formed but is unable to release neutron-rich ejecta, while it produces a relatively mild explosion powered by a neutrino-driven wind blown off its surface. With ejecta electron fractions close to 0.5, all models presumably produce large amounts of Ni56. Our results suggest that collapsar models based on the alpha-viscosity are inefficient r-process sites and that genuinely magnetohydrodynamic effects may be required to generate neutron-rich outflows. A relatively weak effective viscosity generated by magnetohydrodynamic turbulence would improve the prospects for obtaining neutron-rich ejecta. △ Less

Submitted 4 August, 2022; v1 submitted 27 May, 2022; originally announced May 2022.

Comments: 7 pages, 4 figures, 1 table, slightly revised discussion, main results unchanged compared to v1, accepted to ApJL

Report number: RIKEN-iTHEMS-Report-22

A Comparison of 2D Magnetohydrodynamic Supernova Simulations with the CoCoNuT-FMT and Aenus-Alcar Codes

Authors: Vishnu Varma, Bernhard Mueller, Martin Obergaulinger

Abstract: Code comparisons are a valuable tool for the verification of supernova simulation codes and the quantification of model uncertainties. Here we present a first comparison of axisymmetric magnetohydrodynamic (MHD) supernova simulations with the CoCoNuT-FMT and Aenus-Alcar codes, which use distinct methods for treating the MHD induction equation and the neutrino transport. We run two sets of simulati… ▽ More Code comparisons are a valuable tool for the verification of supernova simulation codes and the quantification of model uncertainties. Here we present a first comparison of axisymmetric magnetohydrodynamic (MHD) supernova simulations with the CoCoNuT-FMT and Aenus-Alcar codes, which use distinct methods for treating the MHD induction equation and the neutrino transport. We run two sets of simulations of a rapidly rotating 35M gamma-ray burst progenitor model with different choices for the initial field strength, namely 10^12 G for the maximum poloidal and toroidal field in the strong-field case and 10^10 G in the weak-field case. We also investigate the influence of the Riemann solver and the resolution in CoCoNuT-FMT. The dynamics is qualitatively similar for both codes and robust with respect to these numerical details, with a rapid magnetorotational explosion in the strong-field case and a delayed neutrino-driven explosion in the weak-field case. Despite relatively similar shock trajectories, we find sizeable differences in many other global metrics of the dynamics, like the explosion energy and the magnetic energy of the proto-neutron star. Further differences emerge upon closer inspection, for example, the disk-like surface structure of the proto-neutron star proves highly sensitives to numerical details. The electron fraction distribution in the ejecta as a crucial determinant for the nucleosynthesis is qualitatively robust, but the extent of neutron- or proton-rich tails is sensitive to numerical details. Due to the complexity of the dynamics, the ultimate cause of model differences can rarely be uniquely identified, but our comparison helps gauge uncertainties inherent in current MHD supernova simulations. △ Less

Submitted 15 October, 2021; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 16 pages, 20 figures, accepted for publication in MNRAS

Magnetorotational core collapse of possible GRB progenitors. IV. A wider range of progenitors

Authors: M. Obergaulinger, M. Á. Aloy

Abstract: The final collapse of the cores of massive stars can lead to a wide variety of outcomes in terms of electromagnetic and kinetic energies, nucleosynthesis, and remnants. The connection of this wide spectrum of explosion and remnant types to the properties of the progenitors remains an open issue. Rotation and magnetic fields in Wolf-Rayet stars of subsolar metallicity may explain extreme events suc… ▽ More The final collapse of the cores of massive stars can lead to a wide variety of outcomes in terms of electromagnetic and kinetic energies, nucleosynthesis, and remnants. The connection of this wide spectrum of explosion and remnant types to the properties of the progenitors remains an open issue. Rotation and magnetic fields in Wolf-Rayet stars of subsolar metallicity may explain extreme events such as superluminous supernovae and gamma-ray bursts powered by proto-magnetars or collapsars. Continuing numerical studies of magnetorotational core collapse including detailed neutrino physics, we focus on progenitors with zero-age main-sequence masses in the range between 5 and 39 solar masses. The pre-collapse stars are one dimensional models employing prescriptions for the effects of rotation and magnetic fields. Eight of the ten stars we consider being the results of chemically homogeneous evolution due to enhanced rotational mixing (Aguilera-Dena et al. 2018). All but one of them produce explosions driven by neutrino heating (more likely for low mass progenitors up to 8 solar masses) and non-spherical flows or by magnetorotational stresses (more frequent above 26 solar masses). In most of them and for the one non-exploding model, ongoing accretion leads to black-hole formation. Rapid rotation makes a subsequent collapsar activity plausible. Models not forming black holes show proto-magnetar driven jets. Conditions for the formation of nickel are more favourable in magneto-rotationally driven models, though our rough estimates fall short of the requirements for extremely bright events if these are powered by radioactive decay. However, approximate light curves of our models suggest that a proto-magnetar or black hole spin-down may fuel luminous transients (with peak luminosities ~ 10^{43...44} erg). △ Less

Submitted 4 April, 2022; v1 submitted 31 August, 2021; originally announced August 2021.

Comments: 21 pages, 14 figures

Journal ref: MNRAS, 512,2, 2489-2507 (2022)

Post-explosion evolution of core-collapse supernovae

Authors: M. Witt, A. Psaltis, H. Yasin, C. Horn, M. Reichert, T. Kuroda, M. Obergaulinger, S. M. Couch, A. Arcones

Abstract: We investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 successful explosions during several seconds. We present a broad study based on three progenitors (11.2 $M_\odot$, 15 $M_\odot$, and 27 $M_\odot$), different neutrino-heating efficie… ▽ More We investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 successful explosions during several seconds. We present a broad study based on three progenitors (11.2 $M_\odot$, 15 $M_\odot$, and 27 $M_\odot$), different neutrino-heating efficiencies, and various rotation rates. We show that the first seconds after shock revival determine the final explosion energy, remnant mass, and properties of ejected matter. Our results suggest that a continued mass accretion increases the explosion energy even at late times. We link the late-time mass accretion to initial conditions such as rotation strength and shock deformation at explosion time. Only some of our simulations develop a neutrino-driven wind that survives for several seconds. This indicates that neutrino-driven winds are not a standard feature expected after every successful explosion. Even if our neutrino treatment is simple, we estimate the nucleosynthesis of the exploding models for the 15 $M_\odot$ progenitor after correcting the neutrino energies and luminosities to get a more realistic electron fraction. △ Less

Submitted 1 July, 2021; originally announced July 2021.

Comments: 16 pages, 15 figures, submitted to ApJ

Three-dimensional core-collapse supernovae with complex magnetic structures: I. Explosion dynamics

Authors: Matteo Bugli, Jérôme Guilet, Martin Obergaulinger

Abstract: Magnetic fields can play a major role in the dynamics of outstanding explosions associated to violent events such as GRBs and hypernovae, since they provide a natural mechanism to harness the rotational energy of the central proto-neutron star and power relativistic jets through the stellar progenitor. As the structure of such fields is quite uncertain, most numerical models of MHD-driven core-col… ▽ More Magnetic fields can play a major role in the dynamics of outstanding explosions associated to violent events such as GRBs and hypernovae, since they provide a natural mechanism to harness the rotational energy of the central proto-neutron star and power relativistic jets through the stellar progenitor. As the structure of such fields is quite uncertain, most numerical models of MHD-driven core-collapse supernovae consider an aligned dipole as initial magnetic field, while the field's morphology can actually be much more complex. We present three-dimensional simulations of core-collapse supernovae with more realistic magnetic structures, such as quadrupolar fields and, for the first time, an equatorial dipolar field. Configurations other than an aligned dipole produce weaker explosions and less collimated outflows, but can at the same time be more efficient in extracting the rotational energy from the PNS. This energy is then stored in the surroundings of the PNS, rather than powering the polar jets. A significant axial dipolar component is also produced by models starting with a quadrupolar field, pointing to an effective dynamo mechanism operating in proximity of the PNS surface. △ Less

Submitted 23 July, 2021; v1 submitted 3 May, 2021; originally announced May 2021.

Comments: 13 pages, 12 figures; accepted for publication in MNRAS

Inference of proto-neutron star properties from gravitational-wave data in core-collapse supernovae

Authors: Marie-Anne Bizouard, Patricio Maturana-Russel, Alejandro Torres-Forné, Martin Obergaulinger, Pablo Cerdá-Durán, Nelson Christensen, José A. Font, Renate Meyer

Abstract: The eventual detection of gravitational waves from core-collapse supernovae (CCSN) will help improve our current understanding of the explosion mechanism of massive stars. The stochastic nature of the late post-bounce gravitational wave signal due to the non-linear dynamics of the matter involved and the large number of degrees of freedom of the phenomenon make the source parameter inference probl… ▽ More The eventual detection of gravitational waves from core-collapse supernovae (CCSN) will help improve our current understanding of the explosion mechanism of massive stars. The stochastic nature of the late post-bounce gravitational wave signal due to the non-linear dynamics of the matter involved and the large number of degrees of freedom of the phenomenon make the source parameter inference problem very challenging. In this paper we take a step towards that goal and present a parameter estimation approach which is based on the gravitational waves associated with oscillations of proto-neutron stars (PNS). Numerical simulations of CCSN have shown that buoyancy-driven g-modes are responsible for a significant fraction of the gravitational wave signal and their time-frequency evolution is linked to the physical properties of the compact remnant through universal relations, as demonstrated in [1]. We use a set of 1D CCSN simulations to build a model that relates the evolution of the PNS properties with the frequency of the dominant g-mode, which is extracted from the gravitational-wave data using a new algorithm we have developed for our study. The model is used to infer the time evolution of a combination of the mass and the radius of the PNS. The performance of the method is estimated employing simulations of 2D CCSN waveforms covering a progenitor mass range between 11 and 40 solar masses and different equations of state. Considering signals embedded in Gaussian gravitational wave detector noise, we show that it is possible to infer PNS properties for a galactic source using Advanced LIGO and Advanced Virgo data at design sensitivities. Third generation detectors such as Einstein Telescope and Cosmic Explorer will allow to test distances of ${\cal O}(100\, {\rm kpc})$. △ Less

Submitted 15 November, 2021; v1 submitted 1 December, 2020; originally announced December 2020.

Journal ref: Phys. Rev. D 103, 063006 (2021)

The chemical signature of jet-driven hypernovae

Authors: J. J. Grimmett, Bernhard Müller, Alexander Heger, Projjwal Banerjee, Martin Obergaulinger

Abstract: Hypernovae powered by magnetic jets launched from the surface of rapidly rotating millisecond magnetars are one of the leading models to explain broad-lined Type Ic supernovae (SNe Ic-BL), and have been implicated as an important source of metal enrichment in the early Universe. We investigate the nucleosynthesis in such jet-driven hypernovae using a parameterised, but physically motivated, approa… ▽ More Hypernovae powered by magnetic jets launched from the surface of rapidly rotating millisecond magnetars are one of the leading models to explain broad-lined Type Ic supernovae (SNe Ic-BL), and have been implicated as an important source of metal enrichment in the early Universe. We investigate the nucleosynthesis in such jet-driven hypernovae using a parameterised, but physically motivated, approach that analytically relates an artificially injected jet energy flux to the power available from the energy in differential rotation in the proto-neutron star. We find ejected $^{56}\mathrm{Ni}$ masses of $0.05\,\mathrm{M}_\odot - 0.45\,\mathrm{M}_\odot$ in our most energetic models with explosion energy $>10^{52}\,\mathrm{erg}$. This is in good agreement with the range of observationally inferred values for SNe Ic-BL. The $^{56}\mathrm{Ni}$ is mostly synthesised in the shocked stellar envelope, and is therefore only moderately sensitive to the jet composition. Jets with a high electron fraction $Y_\mathrm{e}=0.5$ eject more $^{56}\mathrm{Ni}$ by a factor of 2 than neutron-rich jets. We can obtain chemical abundance profiles in good agreement with the average chemical signature observed in extremely metal-poor (EMP) stars presumably polluted by hypernova ejecta. Notably, $\mathrm{[Zn/Fe]} \gtrsim 0.5$ is consistently produced in our models. For neutron-rich jets, there is a significant r-process component, and agreement with EMP star abundances in fact requires either a limited contribution from neutron-rich jets or a stronger dilution of r-process material in the interstellar medium than for the slow SN ejecta outside the jet. The high $\mathrm{[C/Fe]}\gtrsim 0.7$ observed in many EMP stars cannot be consistently achieved due to the large mass of iron in the ejecta, however, and remains a challenge for jet-driven hypernovae based on the magneto-rotational mechanism. △ Less

Submitted 6 December, 2020; v1 submitted 13 October, 2020; originally announced October 2020.

Nucleosynthesis in magneto-rotational supernovae

Authors: Moritz Reichert, Martin Obergaulinger, Marius Eichler, Miguel-Ángel Aloy, Almudena Arcones

Abstract: We present the nucleosynthesis of magneto-rotational supernovae (MR-SNe) including neutrino-driven and magneto-rotational-driven ejecta based, for the first time, on 2D simulations with accurate neutrino transport. The models analysed here have different rotation and magnetic fields, allowing us to explore the impact of these two key ingredients. The accurate neutrino transport of the simulations… ▽ More We present the nucleosynthesis of magneto-rotational supernovae (MR-SNe) including neutrino-driven and magneto-rotational-driven ejecta based, for the first time, on 2D simulations with accurate neutrino transport. The models analysed here have different rotation and magnetic fields, allowing us to explore the impact of these two key ingredients. The accurate neutrino transport of the simulations is critical to analyse the slightly neutron-rich and proton-rich ejecta that are similar to the, also neutrino-driven, ejecta in standard supernovae. In the model with strong magnetic field, the r-process produces heavy elements up to the third r-process peak ($A\sim 195$), in agreement with previous works. This model presents a jet-like explosion with proton-rich jets surrounded by neutron-rich material where the r-process occurs. We have estimated a lower limit for $^{56}$Ni of $2.5\times10^{-2} M_\odot$, which is still well below the expected hypernova value. Longer simulations including the accretion disc evolution are required to get a final prediction. In addition, we have found that the late evolution is critical in a model with weak magnetic field in which late-ejected neutron-rich matter produces elements up to the second r-process peak. Even if we cannot yet provide conclusions for hypernova nucleosynthesis, our results agree with observations of old stars and radioactive isotopes in supernova remnants. This makes MR-SNe a good additional scenario to neutron star mergers for the synthesis of heavy elements and brings us closer to understand their origin and the role of MR-SNe in the early Galaxy nucleosynthesis. △ Less

Submitted 20 January, 2021; v1 submitted 5 October, 2020; originally announced October 2020.

Magnetorotational core collapse of possible GRB progenitors. III. Three-dimensional models

Authors: Martin Obergaulinger, Miguel-Ángel Aloy

Abstract: We explore the influence of non-axisymmetric modes on the dynamics of the collapsed core of rotating, magnetized high-mass stars in three-dimensional simulations of a rapidly rotating star with an initial mass of $M_{ZAMS}$ = 35 solar masses endowed with four different pre-collapse configurations of the magnetic field, ranging from moderate to very strong field strength and including the field pre… ▽ More We explore the influence of non-axisymmetric modes on the dynamics of the collapsed core of rotating, magnetized high-mass stars in three-dimensional simulations of a rapidly rotating star with an initial mass of $M_{ZAMS}$ = 35 solar masses endowed with four different pre-collapse configurations of the magnetic field, ranging from moderate to very strong field strength and including the field predicted by the stellar evolution model. The model with the weakest magnetic field achieves shock revival due to neutrino heating in a gain layer characterized by a large-scale, hydrodynamic m = 1 spiral mode. Later on, the growing magnetic field of the proto-neutron star launches weak outflows into the early ejecta. Their orientation follows the evolution of the rotational axis of the proto-neutron star, which starts to tilt from the original orientation due to the asymmetric accretion flows impinging on its surface. The models with stronger magnetization generate mildly relativistic, magnetically driven polar outflows propagating over a distance of $10^4$ km within a few 100 ms. These jets are stabilized against disruptive non-axisymmetric instabilities by their fast propagation and by the shear of their toroidal magnetic field. Within the simulation times of around 1 s, the explosions reach moderate energies and the growth of the proto-neutron star masses ceases at values substantially below the threshold for black hole formation, which, in combination with the high rotational energies, might suggest a possible later proto-magnetar activity. △ Less

Submitted 1 February, 2021; v1 submitted 17 August, 2020; originally announced August 2020.

Comments: 22 pages, 24 figures; accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 503, Issue 4, pp.4942-4963 (2021)

Magnetorotational core collapse of possible GRB progenitors. II. Formation of protomagnetars and collapsars

Authors: Miguel-Ángel Aloy, Martin Obergaulinger

Abstract: We assess the variance of the post-collapse evolution remnants of compact, massive, low-metallicity stars, under small changes in the degrees of rotation and magnetic field of selected pre-supernova cores. These stellar models are commonly considered progenitors of long gamma-ray bursts. The fate of the proto-neutron star (PNS) formed after collapse, whose mass may continuously grow due to accreti… ▽ More We assess the variance of the post-collapse evolution remnants of compact, massive, low-metallicity stars, under small changes in the degrees of rotation and magnetic field of selected pre-supernova cores. These stellar models are commonly considered progenitors of long gamma-ray bursts. The fate of the proto-neutron star (PNS) formed after collapse, whose mass may continuously grow due to accretion, critically depends on the poloidal magnetic field strength at bounce. Should the poloidal magnetic field be sufficiently weak, the PNS collapses to a black hole (BH) within a few seconds. Models on this evolutionary track contain promising collapsar engines. Poloidal magnetic fields smooth over large radial scales (e.g. dipolar fields) or slightly augmented with respect to the original pre-supernova core yield long-lasting PNSs. In these models, BH formation is avoided or staved off for a long time, hence, they may produce proto-magnetars (PMs). Some of our PM candidates have been run for $\lesssim 10\,$s after core bounce, but they have not entered the Kelvin-Helmholtz phase yet. Among these models, some display episodic events of spin-down during which we find properties broadly compatible with the theoretical expectations for PMs ($M_{PNS} \approx 1.85\,M_\odot - 2.5\,M_\odot$, $\bar{P}_{PNS} \approx 1.5 - 4\,$ms, and $b^{\rm surf}_{PNS} \lesssim 10^{15}\,$G) and their very collimated supernova ejecta has nearly reached the stellar surface with (still growing) explosion energies $\gtrsim 2\times 10^{51}\,$erg. △ Less

Submitted 16 October, 2020; v1 submitted 9 August, 2020; originally announced August 2020.

Comments: 33 pages, 23 figures, accepted by MNRAS

Numerical viscosity in simulations of the two-dimensional Kelvin-Helmholtz instability

Authors: Martin Obergaulinger, Miguel-Ángel Aloy

Abstract: The Kelvin-Helmholtz instability serves as a simple, well-defined setup for assessing the accuracy of different numerical methods for solving the equations of hydrodynamics. We use it to extend our previous analysis of the convergence and the numerical dissipation in models of the propagation of waves and in the tearing-mode instability in magnetohydrodynamic models. To this end, we perform two-di… ▽ More The Kelvin-Helmholtz instability serves as a simple, well-defined setup for assessing the accuracy of different numerical methods for solving the equations of hydrodynamics. We use it to extend our previous analysis of the convergence and the numerical dissipation in models of the propagation of waves and in the tearing-mode instability in magnetohydrodynamic models. To this end, we perform two-dimensional simulations with and without explicit physical viscosity at different resolutions. A comparison of the growth of the modes excited by our initial perturbations allows us to estimate the effective numerical viscosity of two spatial reconstruction schemes (fifth-order monotonicity preserving and second-order piecewise linear schemes). △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: 9 pages, 4 figures. Proceedings for the "ASTRONUM 2019" conference, July 2019, Paris, France

GW190425: Observation of a Compact Binary Coalescence with Total Mass $\sim 3.4 M_{\odot}$

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, S. Anand , et al. (1177 additional authors not shown)

Abstract: On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from 1.12 to 2.52 $M_{\odot}$ (1.45 to 1.88 $M_{\odot}$ if w… ▽ More On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from 1.12 to 2.52 $M_{\odot}$ (1.45 to 1.88 $M_{\odot}$ if we restrict the dimensionless component spin magnitudes to be smaller than 0.05). These mass parameters are consistent with the individual binary components being neutron stars. However, both the source-frame chirp mass $1.44^{+0.02}_{-0.02} M_{\odot}$ and the total mass $3.4^{+0.3}_{-0.1}\,M_{\odot}$ of this system are significantly larger than those of any other known binary neutron star system. The possibility that one or both binary components of the system are black holes cannot be ruled out from gravitational-wave data. We discuss possible origins of the system based on its inconsistency with the known Galactic binary neutron star population. Under the assumption that the signal was produced by a binary neutron star coalescence, the local rate of neutron star mergers is updated to $250-2810 \text{Gpc}^{-3}\text{yr}^{-1}$. △ Less

Submitted 7 April, 2020; v1 submitted 6 January, 2020; originally announced January 2020.

Comments: 24 pages, 19 figures, published in ApJL

Report number: LIGO-P190425

Journal ref: Astrophysical Journal Letters 892 (2020) L3

The impact of non-dipolar magnetic fields in core-collapse supernovae

Authors: M. Bugli, J. Guilet, M. Obergaulinger, P. Cerdá-Durán, M. Á. Aloy

Abstract: The magnetic field is believed to play an important role in at least some core-collapse supernovae if its magnitude reaches $10^{15}\,\rm{G}$, which is a typical value for a magnetar. In the presence of fast rotation, such a strong magnetic field can drive powerful jet-like explosions if it has the large-scale coherence of a dipole. The topology of the magnetic field is, however, probably much mor… ▽ More The magnetic field is believed to play an important role in at least some core-collapse supernovae if its magnitude reaches $10^{15}\,\rm{G}$, which is a typical value for a magnetar. In the presence of fast rotation, such a strong magnetic field can drive powerful jet-like explosions if it has the large-scale coherence of a dipole. The topology of the magnetic field is, however, probably much more complex with strong multipolar and small-scale components and the consequences for the explosion are so far unclear. We investigate the effects of the magnetic field topology on the dynamics of core-collapse supernovae and the properties of forming proto-neutron star (PNS) by comparing pre-collapse fields of different multipolar orders and radial profiles. Using axisymmetric special relativistic MHD simulations and a two-moment neutrino transport, we find that higher multipolar magnetic configurations lead to generally less energetic explosions, slower expanding shocks and less collimated outflows. Models with a low order multipolar configuration tend to produce more oblate PNS, surrounded in some cases by a rotationally supported toroidal structure of neutron-rich material. Moreover, magnetic fields which are distributed on smaller angular scales produce more massive and faster rotating central PNS, suggesting that higher-order multipolar configurations tend to decrease the efficiency of the magnetorotational launching mechanism. Even if our dipolar models systematically display a far more efficient extraction of the rotational energy of the PNS, fields distributed on smaller angular scales are still capable of powering magnetorotational explosions and shape the evolution of the central compact object. △ Less

Submitted 16 January, 2020; v1 submitted 6 September, 2019; originally announced September 2019.

Comments: 15 pages, 15 figures; published in MNRAS

Journal ref: MNRAS, Volume 492, Issue 1, p.58-71, February 2020

Magnetorotational core collapse of possible GRB progenitors. I. Explosion mechanisms

Authors: M. Obergaulinger, M. Á. Aloy

Abstract: We investigate the explosion of stars with zero-age main-sequence masses between 20 and 35 solar masses and varying degrees of rotation and magnetic fields including ones commonly considered progenitors of gamma-ray bursts (GRBs). The simulations, combining special relativistic magnetohydrodynamics, a general relativistic approximate gravitational potential, and two-moment neutrino transport, demo… ▽ More We investigate the explosion of stars with zero-age main-sequence masses between 20 and 35 solar masses and varying degrees of rotation and magnetic fields including ones commonly considered progenitors of gamma-ray bursts (GRBs). The simulations, combining special relativistic magnetohydrodynamics, a general relativistic approximate gravitational potential, and two-moment neutrino transport, demonstrate the viability of different scenarios for the post-bounce evolution. Having formed a highly massive proto-neutron star (PNS), several models launch successful explosions, either by the standard supernova mechanism based on neutrino heating and hydrodynamic instabilities or by magnetorotational processes. It is, however, quite common for the PNS to collapse to a black hole (BH) within a few seconds. Others might produce proto-magnetar-driven explosions. We explore several ways to describe the different explosion mechanisms. The competition between the timescales for advection of gas through the gain layer and heating by neutrinos provides an approximate explanation for models with insignificant magnetic fields. The fidelity of this explosion criterion in the case of rapid rotation can be improved by accounting for the strong deviations from spherical symmetry and mixing between pole and equator. We furthermore study an alternative description including the ram pressure of the gas falling through the shock. Magnetically driven explosions tend to arise from a strongly magnetised region around the polar axis. In these cases, the onset of the explosion corresponds to the equality between the advection timescale and the timescale for the propagation of Alfvén waves through the gain layer. △ Less

Submitted 11 January, 2020; v1 submitted 3 September, 2019; originally announced September 2019.

Comments: 23 pages, 19 figures. Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 492, Issue 4, p.4613-4634 (2020)

A guide to LIGO-Virgo detector noise and extraction of transient gravitational-wave signals

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Alford, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, A. Ananyeva , et al. (1113 additional authors not shown)

Abstract: The LIGO Scientific Collaboration and the Virgo Collaboration have cataloged eleven confidently detected gravitational-wave events during the first two observing runs of the advanced detector era. All eleven events were consistent with being from well-modeled mergers between compact stellar-mass objects: black holes or neutron stars. The data around the time of each of these events have been made… ▽ More The LIGO Scientific Collaboration and the Virgo Collaboration have cataloged eleven confidently detected gravitational-wave events during the first two observing runs of the advanced detector era. All eleven events were consistent with being from well-modeled mergers between compact stellar-mass objects: black holes or neutron stars. The data around the time of each of these events have been made publicly available through the gravitational-wave open science center. The entirety of the gravitational-wave strain data from the first and second observing runs have also now been made publicly available. There is considerable interest among the broad scientific community in understanding the data and methods used in the analyses. In this paper, we provide an overview of the detector noise properties and the data analysis techniques used to detect gravitational-wave signals and infer the source properties. We describe some of the checks that are performed to validate the analyses and results from the observations of gravitational-wave events. We also address concerns that have been raised about various properties of LIGO-Virgo detector noise and the correctness of our analyses as applied to the resulting data. △ Less

Submitted 10 February, 2020; v1 submitted 29 August, 2019; originally announced August 2019.

Journal ref: B P Abbott et al 2020 Class. Quantum Grav. 37 055002

A gravitational-wave measurement of the Hubble constant following the second observing run of Advanced LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, S. Anand , et al. (1164 additional authors not shown)

Abstract: This paper presents the gravitational-wave measurement of the Hubble constant ($H_0$) using the detections from the first and second observing runs of the Advanced LIGO and Virgo detector network. The presence of the transient electromagnetic counterpart of the binary neutron star GW170817 led to the first standard-siren measurement of $H_0$. Here we additionally use binary black hole detections i… ▽ More This paper presents the gravitational-wave measurement of the Hubble constant ($H_0$) using the detections from the first and second observing runs of the Advanced LIGO and Virgo detector network. The presence of the transient electromagnetic counterpart of the binary neutron star GW170817 led to the first standard-siren measurement of $H_0$. Here we additionally use binary black hole detections in conjunction with galaxy catalogs and report a joint measurement. Our updated measurement is $H_0 = 68.7^{+17.0}_{-7.8}$ km/s/Mpc (68.3\% of the highest density posterior interval with a flat-in-log prior) which is an improvement by a factor of 1.04 (about 4\%) over the GW170817-only value of $68.7^{+17.5}_{-8.3}$ km/s/Mpc. A significant additional contribution currently comes from GW170814, a loud and well-localized detection from a part of the sky thoroughly covered by the Dark Energy Survey. With numerous detections anticipated over the upcoming years, an exhaustive understanding of other systematic effects are also going to become increasingly important. These results establish the path to cosmology using gravitational-wave observations with and without transient electromagnetic counterparts. △ Less

Submitted 8 November, 2021; v1 submitted 16 August, 2019; originally announced August 2019.

Comments: 21 pages, 8 figures; this version corrects Fig 2; there are minor changes to also Figs 3 & 4 and the final results

Report number: LIGO-P1900015

Journal ref: Astrophys J 909 Number 2 218 (2021)

An Optically Targeted Search for Gravitational Waves emitted by Core-Collapse Supernovae during the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, S. Anand, A. Ananyeva , et al. (1173 additional authors not shown)

Abstract: We present the results from a search for gravitational-wave transients associated with core-collapse supernovae observed within a source distance of approximately 20 Mpc during the first and second observing runs of Advanced LIGO and Advanced Virgo. No significant gravitational-wave candidate was detected. We report the detection efficiencies as a function of the distance for waveforms derived fro… ▽ More We present the results from a search for gravitational-wave transients associated with core-collapse supernovae observed within a source distance of approximately 20 Mpc during the first and second observing runs of Advanced LIGO and Advanced Virgo. No significant gravitational-wave candidate was detected. We report the detection efficiencies as a function of the distance for waveforms derived from multidimensional numerical simulations and phenomenological extreme emission models. For neutrino-driven explosions the distance at which we reach 50% detection efficiency is approaching 5 kpc, and for magnetorotationally-driven explosions is up to 54 kpc. However, waveforms for extreme emission models are detectable up to 28 Mpc. For the first time, the gravitational-wave data enabled us to exclude part of the parameter spaces of two extreme emission models with confidence up to 83%, limited by coincident data coverage. Besides, using ad hoc harmonic signals windowed with Gaussian envelopes we constrained the gravitational-wave energy emitted during core-collapse at the levels of $4.27\times 10^{-4}\,M_\odot c^2$ and $1.28\times 10^{-1}\,M_\odot c^2$ for emissions at 235 Hz and 1304 Hz respectively. These constraints are two orders of magnitude more stringent than previously derived in the corresponding analysis using initial LIGO, initial Virgo and GEO 600 data. △ Less

Submitted 20 August, 2019; v1 submitted 9 August, 2019; originally announced August 2019.

Comments: 13 pages, 5 figures

Report number: LIGO-P1700177

Journal ref: Phys. Rev. D 101, 084002 (2020)

Model comparison from LIGO-Virgo data on GW170817's binary components and consequences for the merger remnant

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, S. Anand, A. Ananyeva , et al. (1169 additional authors not shown)

Abstract: GW170817 is the very first observation of gravitational waves originating from the coalescence of two compact objects in the mass range of neutron stars, accompanied by electromagnetic counterparts, and offers an opportunity to directly probe the internal structure of neutron stars. We perform Bayesian model selection on a wide range of theoretical predictions for the neutron star equation of stat… ▽ More GW170817 is the very first observation of gravitational waves originating from the coalescence of two compact objects in the mass range of neutron stars, accompanied by electromagnetic counterparts, and offers an opportunity to directly probe the internal structure of neutron stars. We perform Bayesian model selection on a wide range of theoretical predictions for the neutron star equation of state. For the binary neutron star hypothesis, we find that we cannot rule out the majority of theoretical models considered. In addition, the gravitational-wave data alone does not rule out the possibility that one or both objects were low-mass black holes. We discuss the possible outcomes in the case of a binary neutron star merger, finding that all scenarios from prompt collapse to long-lived or even stable remnants are possible. For long-lived remnants, we place an upper limit of 1.9 kHz on the rotation rate. If a black hole was formed any time after merger and the coalescing stars were slowly rotating, then the maximum baryonic mass of non-rotating neutron stars is at most 3.05 $M_\odot$, and three equations of state considered here can be ruled out. We obtain a tighter limit of 2.67 $M_\odot$ for the case that the merger results in a hypermassive neutron star. △ Less

Submitted 6 March, 2020; v1 submitted 2 August, 2019; originally announced August 2019.

Comments: 35 pages, 4 figures

Report number: LIGO-P1800379

Journal ref: Classical and Quantum Gravity, Vol. 37, No 4, p 045006 (2020)

Search for Eccentric Binary Black Hole Mergers with Advanced LIGO and Advanced Virgo during their First and Second Observing Runs

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, S. Anand , et al. (1161 additional authors not shown)

Abstract: When formed through dynamical interactions, stellar-mass binary black holes may retain eccentric orbits ($e>0.1$ at 10 Hz) detectable by ground-based gravitational-wave detectors. Eccentricity can therefore be used to differentiate dynamically-formed binaries from isolated binary black hole mergers. Current template-based gravitational-wave searches do not use waveform models associated to eccentr… ▽ More When formed through dynamical interactions, stellar-mass binary black holes may retain eccentric orbits ($e>0.1$ at 10 Hz) detectable by ground-based gravitational-wave detectors. Eccentricity can therefore be used to differentiate dynamically-formed binaries from isolated binary black hole mergers. Current template-based gravitational-wave searches do not use waveform models associated to eccentric orbits, rendering the search less efficient to eccentric binary systems. Here we present results of a search for binary black hole mergers that inspiral in eccentric orbits using data from the first and second observing runs (O1 and O2) of Advanced LIGO and Advanced Virgo. The search uses minimal assumptions on the morphology of the transient gravitational waveform. We show that it is sensitive to binary mergers with a detection range that is weakly dependent on eccentricity for all bound systems. Our search did not identify any new binary merger candidates. We interpret these results in light of eccentric binary formation models. △ Less

Submitted 24 January, 2020; v1 submitted 22 July, 2019; originally announced July 2019.

Comments: 7 pages, 2 figures

Report number: LIGO Document P1900110

Search for gravitational-wave signals associated with gamma-ray bursts during the second observing run of Advanced LIGO and Advanced Virgo

Authors: B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, S. Anand, A. Ananyeva, S. B. Anderson , et al. (1174 additional authors not shown)

Abstract: We present the results of targeted searches for gravitational-wave transients associated with gamma-ray bursts during the second observing run of Advanced LIGO and Advanced Virgo, which took place from 2016 November to 2017 August. We have analyzed 98 gamma-ray bursts using an unmodeled search method that searches for generic transient gravitational waves and 42 with a modeled search method that t… ▽ More We present the results of targeted searches for gravitational-wave transients associated with gamma-ray bursts during the second observing run of Advanced LIGO and Advanced Virgo, which took place from 2016 November to 2017 August. We have analyzed 98 gamma-ray bursts using an unmodeled search method that searches for generic transient gravitational waves and 42 with a modeled search method that targets compact-binary mergers as progenitors of short gamma-ray bursts. Both methods clearly detect the previously reported binary merger signal GW170817, with p-values of $<9.38 \times 10^{-6}$ (modeled) and $3.1 \times 10^{-4}$ (unmodeled). We do not find any significant evidence for gravitational-wave signals associated with the other gamma-ray bursts analyzed, and therefore we report lower bounds on the distance to each of these, assuming various source types and signal morphologies. Using our final modeled search results, short gamma-ray burst observations, and assuming binary neutron star progenitors, we place bounds on the rate of short gamma-ray bursts as a function of redshift for $z \leq 1$. We estimate 0.07-1.80 joint detections with Fermi-GBM per year for the 2019-20 LIGO-Virgo observing run and 0.15-3.90 per year when current gravitational-wave detectors are operating at their design sensitivities. △ Less

Submitted 22 November, 2019; v1 submitted 2 July, 2019; originally announced July 2019.

Report number: LIGO-P1900034

Journal ref: Astrophys. J. 886, 75 (2019)

Search for intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato , et al. (1174 additional authors not shown)

Abstract: Gravitational wave astronomy has been firmly established with the detection of gravitational waves from the merger of ten stellar mass binary black holes and a neutron star binary. This paper reports on the all-sky search for gravitational waves from intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network. The search uses three independen… ▽ More Gravitational wave astronomy has been firmly established with the detection of gravitational waves from the merger of ten stellar mass binary black holes and a neutron star binary. This paper reports on the all-sky search for gravitational waves from intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network. The search uses three independent algorithms: two based on matched filtering of the data with waveform templates of gravitational wave signals from compact binaries, and a third, model-independent algorithm that employs no signal model for the incoming signal. No intermediate mass black hole binary event was detected in this search. Consequently, we place upper limits on the merger rate density for a family of intermediate mass black hole binaries. In particular, we choose sources with total masses $M=m_1+m_2\in[120,800]$M$_\odot$ and mass ratios $q = m_2/m_1 \in[0.1,1.0]$. For the first time, this calculation is done using numerical relativity waveforms (which include higher modes) as models of the real emitted signal. We place a most stringent upper limit of $0.20$~Gpc$^{-3}$yr$^{-1}$ (in co-moving units at the 90% confidence level) for equal-mass binaries with individual masses $m_{1,2}=100$M$_\odot$ and dimensionless spins $χ_{1,2}= 0.8$ aligned with the orbital angular momentum of the binary. This improves by a factor of $\sim 5$ that reported after Advanced LIGO's first observing run. △ Less

Submitted 24 January, 2020; v1 submitted 19 June, 2019; originally announced June 2019.

Comments: main paper: 14 pages, 2 figures and 1 table : total with appendices 19 pages, 2 figures and 2 tables

Report number: LIGO-P1900045

Journal ref: Phys. Rev. D 100, 064064 (2019)

All-sky search for short gravitational-wave bursts in the second Advanced LIGO and Advanced Virgo run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, S. Anand , et al. (1164 additional authors not shown)

Abstract: We present the results of a search for short-duration gravitational-wave transients in the data from the second observing run of Advanced LIGO and Advanced Virgo. We search for gravitational-wave transients with a duration of milliseconds to approximately one second in the 32-4096 Hz frequency band with minimal assumptions about the signal properties, thus targeting a wide variety of sources. We a… ▽ More We present the results of a search for short-duration gravitational-wave transients in the data from the second observing run of Advanced LIGO and Advanced Virgo. We search for gravitational-wave transients with a duration of milliseconds to approximately one second in the 32-4096 Hz frequency band with minimal assumptions about the signal properties, thus targeting a wide variety of sources. We also perform a matched-filter search for gravitational-wave transients from cosmic string cusps for which the waveform is well-modeled. The unmodeled search detected gravitational waves from several binary black hole mergers which have been identified by previous analyses. No other significant events have been found by either the unmodeled search or the cosmic string search. We thus present search sensitivity for a variety of signal waveforms and report upper limits on the source rate-density as function of the characteristic frequency of the signal. These upper limits are a factor of three lower than the first observing run, with a $50\%$ detection probability for gravitational-wave emissions with energies of $\sim10^{-9}M_\odot c^2$ at 153 Hz. For the search dedicated to cosmic string cusps we consider several loop distribution models, and present updated constraints from the same search done in the first observing run. △ Less

Submitted 9 May, 2019; originally announced May 2019.

Comments: 10 pages, 6 figures

Report number: LIGO-P1800308

Journal ref: Phys. Rev. D 100, 024017 (2019)

Search for sub-solar mass ultracompact binaries in Advanced LIGO's second observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, S. Anand , et al. (1165 additional authors not shown)

Abstract: We present an Advanced LIGO and Advanced Virgo search for sub-solar mass ultracompact objects in data obtained during Advanced LIGO's second observing run. In contrast to a previous search of Advanced LIGO data from the first observing run, this search includes the effects of component spin on the gravitational waveform. We identify no viable gravitational wave candidates consistent with sub-solar… ▽ More We present an Advanced LIGO and Advanced Virgo search for sub-solar mass ultracompact objects in data obtained during Advanced LIGO's second observing run. In contrast to a previous search of Advanced LIGO data from the first observing run, this search includes the effects of component spin on the gravitational waveform. We identify no viable gravitational wave candidates consistent with sub-solar mass ultracompact binaries with at least one component between 0.2 - 1.0 solar masses. We use the null result to constrain the binary merger rate of (0.2 solar mass, 0.2 solar mass) binaries to be less than 3.7 x 10^5 Gpc^-3 yr^-1 and the binary merger rate of (1.0 solar mass, 1.0 solar mass) binaries to be less than 5.2 x 10^3 Gpc^-3 yr^-1. Sub-solar mass ultracompact objects are not expected to form via known stellar evolution channels, though it has been suggested that primordial density fluctuations or particle dark matter with cooling mechanisms and/or nuclear interactions could form black holes with sub-solar masses. Assuming a particular primordial black hole formation model, we constrain a population of merging 0.2 solar mass black holes to account for less than 16% of the dark matter density and a population of merging 1.0 solar mass black holes to account for less than 2% of the dark matter density. We discuss how constraints on the merger rate and dark matter fraction may be extended to arbitrary black hole population models that predict sub-solar mass binaries. △ Less

Submitted 25 May, 2019; v1 submitted 18 April, 2019; originally announced April 2019.

Report number: LIGO-P1900037

Journal ref: Phys. Rev. Lett. 123, 161102 (2019)

All-sky search for continuous gravitational waves from isolated neutron stars using Advanced LIGO O2 data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, A. Ananyeva , et al. (1108 additional authors not shown)

Abstract: We present results of an all-sky search for continuous gravitational waves (CWs), which can be produced by fast-spinning neutron stars with an asymmetry around their rotation axis, using data from the second observing run of the Advanced LIGO detectors. We employ three different semi-coherent methods ($\textit{FrequencyHough}$, $\textit{SkyHough}$, and $\textit{Time-Domain $\mathcal{F}… ▽ More We present results of an all-sky search for continuous gravitational waves (CWs), which can be produced by fast-spinning neutron stars with an asymmetry around their rotation axis, using data from the second observing run of the Advanced LIGO detectors. We employ three different semi-coherent methods ($\textit{FrequencyHough}$, $\textit{SkyHough}$, and $\textit{Time-Domain $\mathcal{F}$-statistic}$) to search in a gravitational-wave frequency band from 20 to 1922 Hz and a first frequency derivative from $-1\times10^{-8}$ to $2\times10^{-9}$ Hz/s. None of these searches has found clear evidence for a CW signal, so we present upper limits on the gravitational-wave strain amplitude $h_0$ (the lowest upper limit on $h_0$ is $1.7\times10^{-25}$ in the 123-124 Hz region) and discuss the astrophysical implications of this result. This is the most sensitive search ever performed over the broad range of parameters explored in this study. △ Less

Submitted 24 May, 2019; v1 submitted 5 March, 2019; originally announced March 2019.

Report number: LIGO-P1900012-v7

Journal ref: Phys. Rev. D 100, 024004 (2019)

Universal relations for gravitational-wave asteroseismology of proto-neutron stars

Authors: Alejandro Torres-Forné, Pablo Cerdá-Durán, Martin Obergaulinger, Bernhard Müller, José A. Font

Abstract: State-of-the-art numerical simulations of core-collapse supernovae reveal that the main source of gravitational waves is the excitation of proto-neutron star modes during post-bounce evolution. In this work we derive universal relations that relate the frequencies of the most common oscillation modes observed, i.e. g-modes, p-modes and the f-mode, with fundamental properties of the system, such as… ▽ More State-of-the-art numerical simulations of core-collapse supernovae reveal that the main source of gravitational waves is the excitation of proto-neutron star modes during post-bounce evolution. In this work we derive universal relations that relate the frequencies of the most common oscillation modes observed, i.e. g-modes, p-modes and the f-mode, with fundamental properties of the system, such as the surface gravity of the proto-neutron star or the mean density in the region enclosed by the shock. These relations are independent of the equation of state, the neutrino treatment, and the progenitor mass and hence can be used to build methods to infer proto-neutron star properties from gravitational-wave observations alone. We outline how these measurements could be done and the constraints that could be placed on the proto-neutron star properties. △ Less

Submitted 8 November, 2021; v1 submitted 26 February, 2019; originally announced February 2019.

Journal ref: Phys. Rev. Lett. 123, 051102 (2019)

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015-2017 LIGO Data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, A. Ananyeva , et al. (1127 additional authors not shown)

Abstract: We present a search for gravitational waves from 221 pulsars with rotation frequencies $\gtrsim 10$ Hz. We use advanced LIGO data from its first and second observing runs spanning 2015-2017, which provides the highest-sensitivity gravitational-wave data so far obtained. In this search we target emission from both the $l = m = 2$ mass quadrupole mode, with a frequency at twice that of the pulsar's… ▽ More We present a search for gravitational waves from 221 pulsars with rotation frequencies $\gtrsim 10$ Hz. We use advanced LIGO data from its first and second observing runs spanning 2015-2017, which provides the highest-sensitivity gravitational-wave data so far obtained. In this search we target emission from both the $l = m = 2$ mass quadrupole mode, with a frequency at twice that of the pulsar's rotation, and from the $l = 2$, $m = 1$ mode, with a frequency at the pulsar rotation frequency. The search finds no evidence for gravitational-wave emission from any pulsar at either frequency. For the $l = m = 2$ mode search, we provide updated upper limits on the gravitational-wave amplitude, mass quadrupole moment, and fiducial ellipticity for 167 pulsars, and the first such limits for a further 55. For 20 young pulsars these results give limits that are below those inferred from the pulsars' spin-down. For the Crab and Vela pulsars our results constrain gravitational-wave emission to account for less than 0.017% and 0.18% of the spin-down luminosity, respectively. For the recycled millisecond pulsar J0711-6830 our limits are only a factor of 1.3 above the spin-down limit, assuming the canonical value of $10^{38}$ kg m$^2$ for the star's moment of inertia, and imply a gravitational-wave-derived upper limit on the star's ellipticity of $1.2\!\times\!10^{-8}$. We also place new limits on the emission amplitude at the rotation frequency of the pulsars. △ Less

Submitted 28 July, 2020; v1 submitted 22 February, 2019; originally announced February 2019.

Comments: Main paper: 39 pages, 9 figures, 2 tables. Accepted for publication in The Astrophysical Journal. Results from the tables are contained in the ancillary file Erratum: 11 pages, 4 figures, 1 table

Report number: LIGO-P1800344

Journal ref: The Astrophysical Journal 879 (2019) 10

Narrow-band search for gravitational waves from known pulsars using the second LIGO observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, A. Ananyeva , et al. (1113 additional authors not shown)

Abstract: Isolated spinning neutron stars, asymmetric with respect to their rotation axis, are expected to be sources of continuous gravitational waves. The most sensitive searches for these sources are based on accurate matched filtering techniques, that assume the continuous wave to be phase-locked with the pulsar beamed emission. While matched filtering maximizes the search sensitivity, a significant sig… ▽ More Isolated spinning neutron stars, asymmetric with respect to their rotation axis, are expected to be sources of continuous gravitational waves. The most sensitive searches for these sources are based on accurate matched filtering techniques, that assume the continuous wave to be phase-locked with the pulsar beamed emission. While matched filtering maximizes the search sensitivity, a significant signal-to-noise ratio loss will happen in case of a mismatch between the assumed and the true signal phase evolution. Narrow-band algorithms allow for a small mismatch in the frequency and spin-down values of the pulsar while integrating coherently the entire data set. In this paper we describe a narrow-band search using LIGO O2 data for the continuous wave emission of 33 pulsars. No evidence for a continuous wave signal has been found and upper-limits on the gravitational wave amplitude, over the analyzed frequency and spin-down volume, have been computed for each of the targets. In this search we have surpassed the spin-down limit for some of the pulsars already present in the O1 LIGO narrow-band search, such as J1400\textminus6325 J1813\textminus1246, J1833\textminus1034, J1952+3252, and for new targets such as J0940\textminus5428 and J1747\textminus2809. For J1400\textminus6325, J1833\textminus1034 and J1747\textminus2809 this is the first time the spin-down limit is surpassed. △ Less

Submitted 6 June, 2019; v1 submitted 22 February, 2019; originally announced February 2019.

Comments: 6 figures, 4 tables

Report number: LIGO-P1800391

Journal ref: Phys. Rev. D 99, 122002 (2019)

Titanium hidden in dust

Authors: A. Iyudin, E. Müller, M. Obergaulinger

Abstract: Cassiopeia A, one of the most intriguing galactic supernova remnants, has been a target of many observational efforts including most recent observations by ALMA, Hubble, Herschel, Spitzer, NuSTAR, Integral, and other observatories. We use recent gamma-ray lines observations of the radioactive products of Cas A supernova explosive nucleosynthesis as well as spectral energy densities derived for Cas… ▽ More Cassiopeia A, one of the most intriguing galactic supernova remnants, has been a target of many observational efforts including most recent observations by ALMA, Hubble, Herschel, Spitzer, NuSTAR, Integral, and other observatories. We use recent gamma-ray lines observations of the radioactive products of Cas A supernova explosive nucleosynthesis as well as spectral energy densities derived for Cas A at infrared wavelengths to speculate about the possibility of radioactive 44Ti being locked into large dust grains. This suggestion is also supported by the possible observation of a pre-supernova outburst about 80 years before the actual Cas A supernova explosion in 1671 AD by Italian astronomer G.D. Cassini. The plausibility of such a scenario is discussed also with reference to recent supernovae, and to the contribution of core-collapse supernovae to the overall dust production in the Galaxy. △ Less

Submitted 6 February, 2019; originally announced February 2019.

Comments: 8 pages, 5 figures; accepted for publication in MNRAS

Search for transient gravitational wave signals associated with magnetar bursts during Advanced LIGO's second observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, A. Ananyeva , et al. (1107 additional authors not shown)

Abstract: We present the results of a search for short and intermediate-duration gravitational-wave signals from four magnetar bursts in Advanced LIGO's second observing run. We find no evidence of a signal and set upper bounds on the root sum squared of the total dimensionless strain ($h_{\text{rss}}$) from incoming intermediate-duration gravitational waves ranging from $1.1 \times 10^{-22}$ at 150 Hz to… ▽ More We present the results of a search for short and intermediate-duration gravitational-wave signals from four magnetar bursts in Advanced LIGO's second observing run. We find no evidence of a signal and set upper bounds on the root sum squared of the total dimensionless strain ($h_{\text{rss}}$) from incoming intermediate-duration gravitational waves ranging from $1.1 \times 10^{-22}$ at 150 Hz to $4.4 \times 10^{-22}$ at 1550 Hz at 50% detection efficiency. From the known distance to the magnetar SGR 1806-20 (8.7 kpc) we can place upper bounds on the isotropic gravitational wave energy of $3.4 \times 10^{44} \text{erg}$ at 150 Hz assuming optimal orientation. This represents an improvement of about a factor of 10 in strain sensitivity from the previous search for such signals, conducted during initial LIGO's sixth science run. The short duration search yielded upper limits of $2.1 \times 10^{44} \,\text{erg}$ for short white noise bursts, and $2.3\times 10^{47}\,\text{erg}$ for $100 \, \text{ms}$ long ringdowns at 1500 Hz, both at 50% detection efficiency. △ Less

Submitted 9 February, 2019; v1 submitted 5 February, 2019; originally announced February 2019.

Report number: LIGO-P1800165

Low-Latency Gravitational Wave Alerts for Multi-Messenger Astronomy During the Second Advanced LIGO and Virgo Observing Run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, A. Ananyeva , et al. (1113 additional authors not shown)

Abstract: Advanced LIGO's second observing run (O2), conducted from November 30, 2016 to August 25, 2017, combined with Advanced Virgo's first observations in August 2017 witnessed the birth of gravitational-wave multi-messenger astronomy. The first ever gravitational-wave detection from the coalescence of two neutron stars, GW170817, and its gamma-ray counterpart, GRB 170817A, led to an electromagnetic fol… ▽ More Advanced LIGO's second observing run (O2), conducted from November 30, 2016 to August 25, 2017, combined with Advanced Virgo's first observations in August 2017 witnessed the birth of gravitational-wave multi-messenger astronomy. The first ever gravitational-wave detection from the coalescence of two neutron stars, GW170817, and its gamma-ray counterpart, GRB 170817A, led to an electromagnetic follow-up of the event at an unprecedented scale. Several teams from across the world searched for EM/neutrino counterparts to GW170817, paving the way for the discovery of optical, X-ray, and radio counterparts. In this article, we describe the online identification of gravitational-wave transients and the distribution of gravitational-wave alerts by the LIGO and Virgo collaborations during O2. We also describe the gravitational-wave observables which were sent in the alerts to enable searches for their counterparts. Finally, we give an overview of the online candidate alerts shared with observing partners during O2. Alerts were issued for 14 candidates, six of which have been confirmed as gravitational-wave events associated with the merger of black holes or neutron stars. Eight of the 14 alerts were issued less than an hour after data acquisition. △ Less

Submitted 11 March, 2019; v1 submitted 10 January, 2019; originally announced January 2019.

Comments: 27 pages, 5 figures, 3 tables, accepted for publication in the Astrophysical Journal

Report number: LIGO-P1800255

First measurement of the Hubble constant from a dark standard siren using the Dark Energy Survey galaxies and the LIGO/Virgo binary-black-hole merger GW170814

Authors: The DES Collaboration, the LIGO Scientific Collaboration, the Virgo Collaboration, M. Soares-Santos, A. Palmese, W. Hartley, J. Annis, J. Garcia-Bellido, O. Lahav, Z. Doctor, M. Fishbach, D. E. Holz, H. Lin, M. E. S. Pereira, A. Garcia, K. Herner, R. Kessler, H. V. Peiris, M. Sako, S. Allam, D. Brout, A. Carnero Rosell, H. Y. Chen, C. Conselice, J. deRose , et al. (1181 additional authors not shown)

Abstract: We present a multi-messenger measurement of the Hubble constant H_0 using the binary-black-hole merger GW170814 as a standard siren, combined with a photometric redshift catalog from the Dark Energy Survey (DES). The luminosity distance is obtained from the gravitational wave signal detected by the LIGO/Virgo Collaboration (LVC) on 2017 August 14, and the redshift information is provided by the DE… ▽ More We present a multi-messenger measurement of the Hubble constant H_0 using the binary-black-hole merger GW170814 as a standard siren, combined with a photometric redshift catalog from the Dark Energy Survey (DES). The luminosity distance is obtained from the gravitational wave signal detected by the LIGO/Virgo Collaboration (LVC) on 2017 August 14, and the redshift information is provided by the DES Year 3 data. Black-hole mergers such as GW170814 are expected to lack bright electromagnetic emission to uniquely identify their host galaxies and build an object-by-object Hubble diagram. However, they are suitable for a statistical measurement, provided that a galaxy catalog of adequate depth and redshift completion is available. Here we present the first Hubble parameter measurement using a black-hole merger. Our analysis results in $H_0 = 75.2^{+39.5}_{-32.4}~{\rm km~s^{-1}~Mpc^{-1}}$, which is consistent with both SN Ia and CMB measurements of the Hubble constant. The quoted 68% credible region comprises 60% of the uniform prior range [20,140] ${\rm km~s^{-1}~Mpc^{-1}}$, and it depends on the assumed prior range. If we take a broader prior of [10,220] ${\rm km~s^{-1}~Mpc^{-1}}$, we find $H_0 = 78^{+ 96}_{-24}~{\rm km~s^{-1}~Mpc^{-1}}$ ($57\%$ of the prior range). Although a weak constraint on the Hubble constant from a single event is expected using the dark siren method, a multifold increase in the LVC event rate is anticipated in the coming years and combinations of many sirens will lead to improved constraints on $H_0$. △ Less

Submitted 22 March, 2019; v1 submitted 6 January, 2019; originally announced January 2019.

Comments: 18 pages, 2 Figures, submitted to ApJL

Report number: FERMILAB-PUB-18-629-AE

Journal ref: The Astrophysical Journal Letters, Volume 876, 1, L7 (2019)

Searches for Continuous Gravitational Waves from Fifteen Supernova Remnants and Fomalhaut b with Advanced LIGO

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, A. Ananyeva , et al. (1107 additional authors not shown)

Abstract: We describe directed searches for continuous gravitational waves from sixteen well localized candidate neutron stars assuming none of the stars has a binary companion. The searches were directed toward fifteen supernova remnants and Fomalhaut~b, an extrasolar planet candidate which has been suggested to be a nearby old neutron star. Each search covered a broad band of frequencies and first and sec… ▽ More We describe directed searches for continuous gravitational waves from sixteen well localized candidate neutron stars assuming none of the stars has a binary companion. The searches were directed toward fifteen supernova remnants and Fomalhaut~b, an extrasolar planet candidate which has been suggested to be a nearby old neutron star. Each search covered a broad band of frequencies and first and second time derivatives. After coherently integrating spans of data from the first Advanced LIGO observing run of 3.5--53.7 days per search, applying data-based vetoes and discounting known instrumental artifacts, we found no astrophysical signals. We set upper limits on intrinsic gravitational wave strain as strict as $1\times10^{-25}$, on fiducial neutron star ellipticity as strict as $2\times10^{-9}$, and on fiducial $r$-mode amplitude as strict as $3\times10^{-8}$. △ Less

Submitted 12 August, 2021; v1 submitted 30 December, 2018; originally announced December 2018.

Comments: 21 pages, 22 figures, updated following erratum

Report number: LIGO-P1800333

Journal ref: Astrophysical Journal 875, 122 (2019)

Binary Black Hole Population Properties Inferred from the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, A. Ananyeva , et al. (1116 additional authors not shown)

Abstract: We present results on the mass, spin, and redshift distributions with phenomenological population models using the ten binary black hole mergers detected in the first and second observing runs completed by Advanced LIGO and Advanced Virgo. We constrain properties of the binary black hole (BBH) mass spectrum using models with a range of parameterizations of the BBH mass and spin distributions. We f… ▽ More We present results on the mass, spin, and redshift distributions with phenomenological population models using the ten binary black hole mergers detected in the first and second observing runs completed by Advanced LIGO and Advanced Virgo. We constrain properties of the binary black hole (BBH) mass spectrum using models with a range of parameterizations of the BBH mass and spin distributions. We find that the mass distribution of the more massive black hole in such binaries is well approximated by models with no more than 1% of black holes more massive than $45\,M_\odot$, and a power law index of $α= {1.3}^{+1.4}_{-1.7}$ (90% credibility). We also show that BBHs are unlikely to be composed of black holes with large spins aligned to the orbital angular momentum. Modelling the evolution of the BBH merger rate with redshift, we show that it is at or increasing with redshift with 93% probability. Marginalizing over uncertainties in the BBH population, we find robust estimates of the BBH merger rate density of $R = {53.2}^{+55.8}_{-28.2}$ Gpc$^{-3}$ yr$^{-1}$ (90% credibility). As the BBH catalog grows in future observing runs, we expect that uncertainties in the population model parameters will shrink, potentially providing insights into the formation of black holes via supernovae, binary interactions of massive stars, stellar cluster dynamics, and the formation history of black holes across cosmic time. △ Less

Submitted 21 September, 2020; v1 submitted 30 November, 2018; originally announced November 2018.

Comments: Accepted Manuscript, with content equivalent to the published version

Report number: LIGO-P1800324

Journal ref: Astrophys. J. Lett. 882, L24 (2019)

GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observed by LIGO and Virgo during the First and Second Observing Runs

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, M. A. Aloy, P. A. Altin, A. Amato, A. Ananyeva , et al. (1126 additional authors not shown)

Abstract: We present the results from three gravitational-wave searches for coalescing compact binaries with component masses above 1$\mathrm{M}_\odot$ during the first and second observing runs of the Advanced gravitational-wave detector network. During the first observing run (O1), from September $12^\mathrm{th}$, 2015 to January $19^\mathrm{th}$, 2016, gravitational waves from three binary black hole mer… ▽ More We present the results from three gravitational-wave searches for coalescing compact binaries with component masses above 1$\mathrm{M}_\odot$ during the first and second observing runs of the Advanced gravitational-wave detector network. During the first observing run (O1), from September $12^\mathrm{th}$, 2015 to January $19^\mathrm{th}$, 2016, gravitational waves from three binary black hole mergers were detected. The second observing run (O2), which ran from November $30^\mathrm{th}$, 2016 to August $25^\mathrm{th}$, 2017, saw the first detection of gravitational waves from a binary neutron star inspiral, in addition to the observation of gravitational waves from a total of seven binary black hole mergers, four of which we report here for the first time: GW170729, GW170809, GW170818 and GW170823. For all significant gravitational-wave events, we provide estimates of the source properties. The detected binary black holes have total masses between $18.6_{-0.7}^{+3.2}\mathrm{M}_\odot$, and $84.4_{-11.1}^{+15.8} \mathrm{M}_\odot$, and range in distance between $320_{-110}^{+120}$ Mpc and $2840_{-1360}^{+1400}$ Mpc. No neutron star - black hole mergers were detected. In addition to highly significant gravitational-wave events, we also provide a list of marginal event candidates with an estimated false alarm rate less than 1 per 30 days. From these results over the first two observing runs, which include approximately one gravitational-wave detection per 15 days of data searched, we infer merger rates at the 90% confidence intervals of $110\, -\, 3840$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$ for binary neutron stars and $9.7\, -\, 101$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$ for binary black holes assuming fixed population distributions, and determine a neutron star - black hole merger rate 90% upper limit of $610$ $\mathrm{Gpc}^{-3}\,\mathrm{y}^{-1}$. △ Less

Submitted 4 October, 2019; v1 submitted 30 November, 2018; originally announced November 2018.

Comments: Total: 49 pages and 17 figures. Strain data are available from the Gravitational Wave Open Science Center at https://doi.org/10.7935/82H3-HH23. Figures, data behind the figures, including posterior samples, are available from https://dcc.ligo.org/LIGO-P1800307/public

Report number: LIGO-P1800307

Journal ref: Phys. Rev. X 9, 031040 (2019)

Search for Multi-messenger Sources of Gravitational Waves and High-energy Neutrinos with Advanced LIGO during its first Observing Run, ANTARES and IceCube

Authors: ANTARES, IceCube, LIGO, Virgo Collaborations, :, A. Albert, M. Andre, M. Anghinolfi, M. Ardid, J. -J. Aubert, J. Aublin, T. Avgitas, B. Baret, J. Barrios-Marti, S. Basa, B. Belhorma, V. Bertin, S. Biagi, R. Bormuth, J. Boumaaza, S. Bourret, M. C. Bouwhuis, H. Brânzas, R. Bruijn, J. Brunner , et al. (1570 additional authors not shown)

Abstract: Astrophysical sources of gravitational waves, such as binary neutron star and black hole mergers or core-collapse supernovae, can drive relativistic outflows, giving rise to non-thermal high-energy emission. High-energy neutrinos are signatures of such outflows. The detection of gravitational waves and high-energy neutrinos from common sources could help establish the connection between the dynami… ▽ More Astrophysical sources of gravitational waves, such as binary neutron star and black hole mergers or core-collapse supernovae, can drive relativistic outflows, giving rise to non-thermal high-energy emission. High-energy neutrinos are signatures of such outflows. The detection of gravitational waves and high-energy neutrinos from common sources could help establish the connection between the dynamics of the progenitor and the properties of the outflow. We searched for associated emission of gravitational waves and high-energy neutrinos from astrophysical transients with minimal assumptions using data from Advanced LIGO from its first observing run O1, and data from the ANTARES and IceCube neutrino observatories from the same time period. We focused on candidate events whose astrophysical origin could not be determined from a single messenger. We found no significant coincident candidate, which we used to constrain the rate density of astrophysical sources dependent on their gravitational wave and neutrino emission processes. △ Less

Submitted 15 November, 2018; v1 submitted 24 October, 2018; originally announced October 2018.

Comments: 20 pages, 2 figures

Journal ref: Astrophys.J. 870 (2019) 134

A Fermi Gamma-ray Burst Monitor Search for Electromagnetic Signals Coincident with Gravitational-Wave Candidates in Advanced LIGO's First Observing Run

Authors: The Fermi Gamma-ray Burst Monitor Team, The LIGO Scientific Collaboration, the Virgo Collaboration, :, E. Burns, A. Goldstein, C. M. Hui, L. Blackburn, M. S. Briggs, V. Connaughton, R. Hamburg, D. Kocevski, P. Veres, C. A. Wilson-Hodge, E. Bissaldi, W. H. Cleveland, M. M. Giles, B. Mailyan, C. A. Meegan, W. A. Paciesas, S. Poolakkil, R. D. Preece, J. L. Racusin, O. J. Roberts, A. von Kienlin , et al. (1139 additional authors not shown)

Abstract: We present a search for prompt gamma-ray counterparts to compact binary coalescence gravitational wave (GW) candidates from Advanced LIGO's first observing run (O1). As demonstrated by the multimessenger observations of GW170817/GRB 170817A, electromagnetic and GW observations provide complementary information about the astrophysical source and, in the case of weaker candidates, may strengthen the… ▽ More We present a search for prompt gamma-ray counterparts to compact binary coalescence gravitational wave (GW) candidates from Advanced LIGO's first observing run (O1). As demonstrated by the multimessenger observations of GW170817/GRB 170817A, electromagnetic and GW observations provide complementary information about the astrophysical source and, in the case of weaker candidates, may strengthen the case for an astrophysical origin. Here we investigate low-significance GW candidates from the O1 compact-binary coalescence searches using the Fermi Gamma-ray Burst Monitor (GBM), leveraging its all-sky and broad energy coverage. Candidates are ranked and compared to background to measure significance. Those with false alarm rates of less than 10^-5 Hz (about one per day) are used as the search sample for gamma-ray follow-up. No GW candidates were found to be coincident with gamma-ray transients independently identified by blind searches of the GBM data. In addition, GW candidate event times were followed up by a separate targeted search of GBM data. Among the resulting GBM events, the two with lowest false alarm rates were the gamma-ray transient GW150914-GBM presented in Connaughton et al. (2016) and a solar flare in chance coincidence with a GW candidate. △ Less

Submitted 18 November, 2019; v1 submitted 5 October, 2018; originally announced October 2018.

Comments: Published in ApJ

Search for gravitational waves from a long-lived remnant of the binary neutron star merger GW170817

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca, M. A. Aloy , et al. (1116 additional authors not shown)

Abstract: One unanswered question about the binary neutron star coalescence GW170817 is the nature of its post-merger remnant. A previous search for post-merger gravitational waves targeted high-frequency signals from a possible neutron star remnant with a maximum signal duration of 500 s. Here we revisit the neutron star remnant scenario with a focus on longer signal durations up until the end of the Secon… ▽ More One unanswered question about the binary neutron star coalescence GW170817 is the nature of its post-merger remnant. A previous search for post-merger gravitational waves targeted high-frequency signals from a possible neutron star remnant with a maximum signal duration of 500 s. Here we revisit the neutron star remnant scenario with a focus on longer signal durations up until the end of the Second Advanced LIGO-Virgo Observing run, 8.5 days after the coalescence of GW170817. The main physical scenario for such emission is the power-law spindown of a massive magnetar-like remnant. We use four independent search algorithms with varying degrees of restrictiveness on the signal waveformand different ways of dealing with noise artefacts. In agreement with theoretical estimates, we find no significant signal candidates. Through simulated signals, we quantify that with the current detector sensitivity, nowhere in the studied parameter space are we sensitive to a signal from more than 1 Mpc away, compared to the actual distance of 40 Mpc. This study however serves as a prototype for post-merger analyses in future observing runs with expected higher sensitivity. △ Less

Submitted 4 October, 2019; v1 submitted 5 October, 2018; originally announced October 2018.

Comments: main paper: 9 pages and 3 figures; total with appendices: 24 pages and 9 figures. Full UL tables available as MRT ancillary files

Report number: LIGO-P1800195

Journal ref: The Astrophysical Journal 875:160 (2019)

Three-Dimensional Core-Collapse Supernova Simulations with Multi-Dimensional Neutrino Transport Compared to the Ray-by-Ray-plus Approximation

Authors: Robert Glas, Oliver Just, H. -Thomas Janka, Martin Obergaulinger

Abstract: Self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) are conducted with the Aenus-Alcar code, comparing, for the first time, calculations with fully multi-dimensional (FMD) neutrino transport and the ray-by-ray-plus (RbR+) approximation, both based on a two-moment solver with algebraic M1 closure. We find good agreement between 3D results with FMD and RbR+ tr… ▽ More Self-consistent, time-dependent supernova (SN) simulations in three spatial dimensions (3D) are conducted with the Aenus-Alcar code, comparing, for the first time, calculations with fully multi-dimensional (FMD) neutrino transport and the ray-by-ray-plus (RbR+) approximation, both based on a two-moment solver with algebraic M1 closure. We find good agreement between 3D results with FMD and RbR+ transport for both tested grid resolutions in the cases of a 20 solar-mass progenitor, which does not explode with the employed simplified set of neutrino opacities, and of an exploding 9 solar-mass model. This is in stark contrast to corresponding axisymmetric (2D) simulations, which confirm previous claims that the RbR+ approximation can foster explosions in 2D in particular in models with powerful axial sloshing of the stalled shock due to the standing accretion shock instability (SASI). However, while local and instantaneous variations of neutrino fluxes and heating rates can still be considerably higher with RbR+ transport in 3D, the time-averaged quantities are very similar to FMD results because of the absence of a fixed, artificial symmetry axis that channels the flow. Therefore, except for stochastic fluctuations, the neutrino signals and the post-bounce evolution of 3D simulations with FMD and RbR+ transport are also very similar, in particular for our calculations with the better grid resolution. Higher spatial resolution has clearly a more important impact than the differences by the two transport treatments. Our results back up the use of the RbR+ approximation for neutrino transport in 3D SN modeling. △ Less

Submitted 5 February, 2019; v1 submitted 26 September, 2018; originally announced September 2018.

Comments: 25 pages, 16 figures; referee comments included, new appendix added; accepted by ApJ

Constraining the p-mode--g-mode tidal instability with GW170817

Authors: The LIGO Scientific Collaboration, The Virgo Collaboration, B. P. Abbott, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, T. Adams, P. Addesso, R. X. Adhikari, V. B. Adya, C. Affeldt, B. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, B. Allen, G. Allen, A. Allocca, M. A. Aloy , et al. (1114 additional authors not shown)

Abstract: We analyze the impact of a proposed tidal instability coupling $p$-modes and $g$-modes within neutron stars on GW170817. This non-resonant instability transfers energy from the orbit of the binary to internal modes of the stars, accelerating the gravitational-wave driven inspiral. We model the impact of this instability on the phasing of the gravitational wave signal using three parameters per sta… ▽ More We analyze the impact of a proposed tidal instability coupling $p$-modes and $g$-modes within neutron stars on GW170817. This non-resonant instability transfers energy from the orbit of the binary to internal modes of the stars, accelerating the gravitational-wave driven inspiral. We model the impact of this instability on the phasing of the gravitational wave signal using three parameters per star: an overall amplitude, a saturation frequency, and a spectral index. Incorporating these additional parameters, we compute the Bayes Factor ($\ln B^{pg}_{!pg}$) comparing our $p$-$g$ model to a standard one. We find that the observed signal is consistent with waveform models that neglect $p$-$g$ effects, with $\ln B^{pg}_{!pg} = 0.03^{+0.70}_{-0.58}$ (maximum a posteriori and 90% credible region). By injecting simulated signals that do not include $p$-$g$ effects and recovering them with the $p$-$g$ model, we show that there is a $\simeq 50\%$ probability of obtaining similar $\ln B^{pg}_{!pg}$ even when $p$-$g$ effects are absent. We find that the $p$-$g$ amplitude for 1.4 $M_\odot$ neutron stars is constrained to $\lesssim \text{few}\times10^{-7}$, with maxima a posteriori near $\sim 10^{-7}$ and $p$-$g$ saturation frequency $\sim 70\, \mathrm{Hz}$. This suggests that there are less than a few hundred excited modes, assuming they all saturate by wave breaking. For comparison, theoretical upper bounds suggest a $p$-$g$ amplitude $\lesssim 10^{-6}$ and $\lesssim 10^{3}$ modes saturating by wave breaking. Thus, the measured constraints only rule out extreme values of the $p$-$g$ parameters. They also imply that the instability dissipates $\lesssim 10^{51}\, \mathrm{ergs}$ over the entire inspiral, i.e., less than a few percent of the energy radiated as gravitational waves. △ Less

Submitted 19 September, 2018; v1 submitted 26 August, 2018; originally announced August 2018.

Comments: 7 pages, 2 figures

Journal ref: Phys. Rev. Lett. 122, 061104 (2019)