Showing 1–50 of 105 results for author: Aloy, M Á

Force-free wave interaction in magnetar magnetospheres: Computational modeling in axisymmetry

Authors: Jens F. Mahlmann, Miguel Á. Aloy, Xinyu Li

Abstract: Crustal quakes of highly magnetized neutron stars can disrupt their magnetospheres, triggering energetic phenomena like X-ray and fast radio bursts (FRBs). Understanding plasma wave dynamics in these extreme environments is vital for predicting energy transport across scales to the radiation length. This study models relativistic plasma wave interaction in magnetar magnetospheres with force-free e… ▽ More Crustal quakes of highly magnetized neutron stars can disrupt their magnetospheres, triggering energetic phenomena like X-ray and fast radio bursts (FRBs). Understanding plasma wave dynamics in these extreme environments is vital for predicting energy transport across scales to the radiation length. This study models relativistic plasma wave interaction in magnetar magnetospheres with force-free electrodynamics simulations. For propagation along curved magnetic field lines, we observe the continuous conversion of Alfvén waves (AWs) to fast magnetosonic (FMS) waves. The conversion efficiency can be up to three times higher when counter-propagating AWs interact in the equatorial region. Alfvén waves generate FMS waves of twice their frequency during their first crossing of the magnetosphere. After the initial transient burst of FMS waves, AWs convert to FMS waves periodically, generating variations on timescales of the magnetospheric AW crossing time. This decaying FMS wave tail carries a significant portion (half) of the total energy emitted. Plastic damping of `bouncing' AWs by the magnetar crust has minimal impact on the FMS efficiency. We discuss the implications of the identified wave phenomena for magnetar observations. Outgoing FMS waves can develop electric zones, potential sources of coherent radiation. Long wavelength FMS waves could generate FRBs through reconnection beyond the light cylinder. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 18 pages, 11 figures, submitted to ApJ

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.

Self-consistent MHD simulation of jet launching in a neutron star - white dwarf merger

Authors: J. Moran-Fraile, F. K. Roepke, R. Pakmor, M. A. Aloy, S. T. Ohlmann, F. R. N. Schneider, G. Leidi

Abstract: The merger of a white dwarf (WD) and a neutron star (NS) is a relatively common event that will produce an observable electromagnetic signal. Furthermore, the compactness of these stellar objects makes them an interesting candidate for gravitational wave (GW) astronomy, potentially being in the frequency range of LISA and other missions. To date, three-dimensional simulations of these mergers have… ▽ More The merger of a white dwarf (WD) and a neutron star (NS) is a relatively common event that will produce an observable electromagnetic signal. Furthermore, the compactness of these stellar objects makes them an interesting candidate for gravitational wave (GW) astronomy, potentially being in the frequency range of LISA and other missions. To date, three-dimensional simulations of these mergers have not fully modelled the WD disruption, or have used lower resolutions and have not included magnetic fields even though they potentially shape the evolution of the merger remnant. In this work, we simulate the merger of a 1.4$M_\odot$ NS with a 1$M_\odot$ carbon oxygen WD in the magnetohydrodynamic moving mesh code \AREPO. We find that the disruption of the WD forms an accretion disk around the NS, and the subsequent accretion by the NS powers the launch of strongly magnetized, mildly relativistic jets perpendicular to the orbital plane. Although the exact properties of the jets could be altered by unresolved physics around the NS, the event could result in a transient with a larger luminosity than kilonovae. We discuss possible connections to fast blue optical transients (FBOTs) and long-duration gamma-ray bursts. We find that the frequency of GWs released during the merger is too high to be detectable by the LISA mission, but suitable for deci-hertz observatories such as LGWA, BBO or DECIGO. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: Accepted for publication in A&A. 13 pages, 11 figures

Journal ref: A&A 681, A41 (2024)

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

The magneto-rotational instability in massive stars

Authors: A. Griffiths, G. Meynet, P. Eggenberger, F. Moyano, M. A. Aloy

Abstract: The interaction between weak magnetic fields and rotation can lead to instabilities that transport angular momentum (AM) and chemical elements affecting the evolution of massive stars. We explore the effects of the magneto-rotational instability (MRI) in massive stars to determine its impact on stellar evolution. We use the GENEC code to simulate the evolution of a 15 $M_\odot$ model at solar meta… ▽ More The interaction between weak magnetic fields and rotation can lead to instabilities that transport angular momentum (AM) and chemical elements affecting the evolution of massive stars. We explore the effects of the magneto-rotational instability (MRI) in massive stars to determine its impact on stellar evolution. We use the GENEC code to simulate the evolution of a 15 $M_\odot$ model at solar metallicity up to the end of oxygen burning. The MRI is computed with different trigger conditions, (depending on the weighting of chemical gradients through an arbitrary but commonly used factor), and with different treatments of meridional circulation as either advective or diffusive. We also compare the MRI with the Tayler-Spruit (TS) dynamo, in models that included both instabilities interacting linearly. The MRI triggers throughout stellar evolution. Its activation is highly sensitive to the treatment of meridional circulation and the existence of chemical gradients. The MRI is very efficient at transporting both matter and AM, leading to noticeable differences in rotation rates and chemical structure, which may be observable in young main sequence stars. While the TS dynamo is the dominant mechanism for transferring AM, the MRI remains relevant in models where both instabilities are included. Extrapolation of our results suggests that models including the MRI tend to develop more compact cores, which likely produce failed explosions and black holes, than models where only the TS dynamo is included (where explosions an neutron stars may be more frequent. The MRI is an important factor in massive star evolution but is very sensitive to the implementation ofother processes in the model. The transport of AM and chemical elements due to the MRI alters the rotation rates and the chemical make-up of the star from the core to the surface, and may change the explodability properties of massive stars. △ Less

Submitted 28 July, 2022; v1 submitted 31 March, 2022; originally announced April 2022.

Comments: 20 pages, 10 Figures, submitted to A&A

Journal ref: A&A 665, A147 (2022)

Diffusivity in force-free simulations of global magnetospheres

Authors: J. F. Mahlmann, M. A. Aloy

Abstract: Abstract: Assuming that the numerical diffusivity triggered by violations of the force-free electrodynamics constraints is a proxy for the physical resistivity, we examine its impact on the overall dynamics of force-free aligned pulsar magnetospheres endowed with an equatorial current sheet. We assess the constraint violations as a diffusivity source. The effects of modifications on electric field… ▽ More Abstract: Assuming that the numerical diffusivity triggered by violations of the force-free electrodynamics constraints is a proxy for the physical resistivity, we examine its impact on the overall dynamics of force-free aligned pulsar magnetospheres endowed with an equatorial current sheet. We assess the constraint violations as a diffusivity source. The effects of modifications on electric fields used to restore force-free conditions are not confined to the equatorial current sheet, but modify the magnetospheric dynamics on timescales shorter than the pulsar rotational period. These corrections propagate especially via a channel that was unexplored, namely, changes induced to the electric charge density, $ρ$. We quantify the global consequences of diffusivity by comparing different techniques to model $ρ$. By default, we combine a conservative $ρ$-evolution with hyperbolic/parabolic cleaning of inaccuracies in the Maxwell equations. As an alternative, we enforce a constrained evolution where $ρ$ is directly computed as the electric field divergence. The conservative approach reduces the Poynting flux dissipated in the equatorial current sheet by an order of magnitude, along with an increase of the pulsar luminosity driven by a shift of the Y-point location. The luminosity changes according to $L_{\rm Y}\propto α^{0.11}$, where $α$ is the ratio of diffusion to advection timescales, controlling the amount of (numerical) diffusivity. Our models suggest interpreting the luminosity dependence on the Y-point location as differences in resistivities encountered at the equatorial current sheet. Alternatively, they could be interpreted in terms of the pair formation multiplicity, $κ$, smaller diffusion being consistent with $κ\gg 1$. △ Less

Submitted 28 September, 2021; originally announced September 2021.

Comments: (17 pages, 9 figures, Accepted by 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)

Computational General Relativistic Force-Free Electrodynamics: II. Characterization of Numerical Diffusivity

Authors: J. F. Mahlmann, M. A. Aloy, V. Mewes, P. Cerdá-Durán

Abstract: Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe's most extreme flows of energy. The discretization of Maxwell's equations - needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling - introduces numerical diffusion. It acts as a source of dissipation independe… ▽ More Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe's most extreme flows of energy. The discretization of Maxwell's equations - needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling - introduces numerical diffusion. It acts as a source of dissipation independent of the system's physical constituents. Understanding the numerical diffusion of scientific codes is the key to classify their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics, by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field ($\mathbf{j}_\parallel$), in combination with the break-down of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g. the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of $\mathbf{j}_\parallel$. △ Less

Submitted 8 March, 2021; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: 15 pages, 9 figures, Accepted by A&A

Journal ref: A&A 647, A58 (2021)

Computational General Relativistic Force-Free Electrodynamics: I. Multi-Coordinate Implementation and Testing

Authors: J. F. Mahlmann, M. A. Aloy, V. Mewes, P. Cerdá-Durán

Abstract: General relativistic force-free electrodynamics is one possible plasma-limit employed to analyze energetic outflows in which strong magnetic fields are dominant over all inertial phenomena. The amazing images of black hole shadows from the galactic center and the M87 galaxy provide a first direct glimpse into the physics of accretion flows in the most extreme environments of the universe. The effi… ▽ More General relativistic force-free electrodynamics is one possible plasma-limit employed to analyze energetic outflows in which strong magnetic fields are dominant over all inertial phenomena. The amazing images of black hole shadows from the galactic center and the M87 galaxy provide a first direct glimpse into the physics of accretion flows in the most extreme environments of the universe. The efficient extraction of energy in the form of collimated outflows or jets from a rotating BH is directly linked to the topology of the surrounding magnetic field. We aim at providing a tool to numerically model the dynamics of such fields in magnetospheres around compact objects, such as black holes and neutron stars. By this, we probe their role in the formation of high energy phenomena such as magnetar flares and the highly variable teraelectronvolt emission of some active galactic nuclei. In this work, we present numerical strategies capable of modeling fully dynamical force-free magnetospheres of compact astrophysical objects. We provide implementation details and extensive testing of our implementation of general relativistic force-free electrodynamics in Cartesian and spherical coordinates using the infrastructure of the Einstein Toolkit. The employed hyperbolic/parabolic cleaning of numerical errors with full general relativistic compatibility allows for fast advection of numerical errors in dynamical spacetimes. Such fast advection of divergence errors significantly improves the stability of the general relativistic force-free electrodynamics modeling of black hole magnetospheres. △ Less

Submitted 8 March, 2021; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: 22 pages, 19 figures, Accepted by A&A

Journal ref: A&A 647, A57 (2021)

Striped Blandford/Znajek jets from advection of small scale magnetic field

Authors: J. F. Mahlmann, A. Levinson, M. A. Aloy

Abstract: Black hole - accretion disc systems are the central engines of relativistic jets from stellar to galactic scales. We numerically quantify the unsteady outgoing Poynting flux through the horizon of a rapidly spinning black hole endowed with a rotating accretion disc. The disc supports small-scale, concentric, flux tubes with zero net magnetic flux. Our General Relativistic force-free electrodynamic… ▽ More Black hole - accretion disc systems are the central engines of relativistic jets from stellar to galactic scales. We numerically quantify the unsteady outgoing Poynting flux through the horizon of a rapidly spinning black hole endowed with a rotating accretion disc. The disc supports small-scale, concentric, flux tubes with zero net magnetic flux. Our General Relativistic force-free electrodynamics simulations follow the accretion onto the black hole over several hundred dynamical timescales in 3D. For the case of counter-rotating accretion discs, the average process efficiency reaches up to $\left\langleε\right\rangle\approx 0.43$, compared to a stationary energy extraction by the Blandford/Znajek process. The process efficiency depends on the cross-sectional area of the loops, i.e. on the product $l\times h$, where $l$ is the radial loop thickness and $h$ its vertical scale height. We identify a strong correlation between efficient electromagnetic energy extraction and the quasi-stationary setting of ideal conditions for the operation of the Blandford/Znajek process (e.g. optimal field line angular velocity and fulfillment of the so-called Znajek condition). Remarkably, the energy extraction operates intermittently (alternating episodes of high and low efficiency) without imposing any large-scale magnetic field embedding the central object. Scaling our results to supermassive black holes, we estimate that the typical variability timescale of the system is of the order of days to months. Such timescales may account for the longest variability scales of TeV emission observed, e.g. in M87. △ Less

Submitted 13 July, 2020; v1 submitted 9 January, 2020; originally announced January 2020.

Comments: 25 pages, 17 figures of which 17 are in color. 1 animated figure. Accepted by MNRAS

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)

Instability of twisted magnetar magnetospheres

Authors: J. F. Mahlmann, T. Akgün, J. A. Pons, M. A. Aloy, P. Cerdá-Durán

Abstract: We present three-dimensional force-free electrodynamics simulations of magnetar magnetospheres that demonstrate the instability of certain degenerate, high energy equilibrium solutions of the Grad-Shafranov equation. This result indicates the existence of an unstable branch of twisted magnetospheric solutions and allows to formulate an instability criterion. The rearrangement of magnetic field lin… ▽ More We present three-dimensional force-free electrodynamics simulations of magnetar magnetospheres that demonstrate the instability of certain degenerate, high energy equilibrium solutions of the Grad-Shafranov equation. This result indicates the existence of an unstable branch of twisted magnetospheric solutions and allows to formulate an instability criterion. The rearrangement of magnetic field lines as a consequence of this instability triggers the dissipation of up to 30% of the magnetospheric energy on a thin layer above the magnetar surface. During this process, we predict an increase of the mechanical stresses onto the stellar crust, which can potentially result in a global mechanical failure of a significant fraction of it. We find that the estimated energy release and the emission properties are compatible with the observed giant flare events. The newly identified instability is a candidate for recurrent energy dissipation, which could explain part of the phenomenology observed in magnetars. △ Less

Submitted 24 November, 2019; v1 submitted 31 July, 2019; originally announced August 2019.

Comments: 21 pages, 11 figures of which 11 are in color. 1 animated figure. Accepted by MNRAS

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)

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)

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)

3D MHD modeling of the expanding remnant of SN 1987A. Role of magnetic field and non-thermal radio emission

Authors: S. Orlando, M. Miceli, O. Petruk, M. Ono, S. Nagataki, M. A. Aloy, P. Mimica, S. -H. Lee, F. Bocchino, G. Peres, M. Guarrasi

Abstract: Aims. We investigate the role played by a pre-supernova (SN) ambient magnetic field on the dynamics of the expanding remnant of SN 1987A and the origin and evolution of the radio emission from the remnant, in particular, during the interaction of the blast wave with the nebula surrounding the SN. Methods. We model the evolution of SN 1987A from the breakout of the shock wave at the stellar surface… ▽ More Aims. We investigate the role played by a pre-supernova (SN) ambient magnetic field on the dynamics of the expanding remnant of SN 1987A and the origin and evolution of the radio emission from the remnant, in particular, during the interaction of the blast wave with the nebula surrounding the SN. Methods. We model the evolution of SN 1987A from the breakout of the shock wave at the stellar surface to the expansion of its remnant through the surrounding nebula by 3D MHD simulations. The model considers the radiative cooling, the deviations from equilibrium of ionization, the deviation from temperature-equilibration between electrons and ions, and a plausible configuration of the pre-SN ambient magnetic field. From the simulations, we synthesize both thermal X-ray and non-thermal radio emission and compare the model results with observations. Results. The presence of an ambient magnetic field reduces the erosion and fragmentation of the dense equatorial ring after the impact of the SN blast wave. As a result, the ring survives the passage of the blast, at least, during the time covered by the simulations (40 years). Our model is able to reproduce the morphology and lightcurves of SN 1987A in both X-ray and radio bands. The model reproduces the observed radio emission if the flux originating from the reverse shock is heavily suppressed. In this case, the radio emission originates mostly from the forward shock traveling through the H II region. Possible mechanisms for the suppression of emission from the reverse shock are investigated. We find that synchrotron self-absorption and free-free absorption have negligible effects on the emission during the interaction with the nebula. We suggest that the emission from the reverse shock at radio frequencies might be limited by highly magnetized ejecta. △ Less

Submitted 30 November, 2018; originally announced December 2018.

Comments: 15 pages, 9 Figures, accepted for publication on A&A, high quality figures will be available on the published paper

Journal ref: A&A 622, A73 (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)

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)

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

Authors: 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, P. A. Altin, A. Amato , et al. (1113 additional authors not shown)

Abstract: We present the first Advanced LIGO and Advanced Virgo search for ultracompact binary systems with component masses between 0.2 $M_\odot$ - 1.0 $M_\odot$ using data taken between September 12, 2015 and January 19, 2016. We find no viable gravitational wave candidates. Our null result constrains the coalescence rate of monochromatic (delta function) distributions of non-spinning (0.2 $M_\odot$, 0.2… ▽ More We present the first Advanced LIGO and Advanced Virgo search for ultracompact binary systems with component masses between 0.2 $M_\odot$ - 1.0 $M_\odot$ using data taken between September 12, 2015 and January 19, 2016. We find no viable gravitational wave candidates. Our null result constrains the coalescence rate of monochromatic (delta function) distributions of non-spinning (0.2 $M_\odot$, 0.2 $M_\odot$) ultracompact binaries to be less than $1.0 \times 10^6 \text{Gpc}^{-3} \text{yr}^{-1}$ and the coalescence rate of a similar distribution of (1.0 $M_\odot$, 1.0 $M_\odot$) ultracompact binaries to be less than $1.9 \times 10^4 \text{Gpc}^{-3} \text{yr}^{-1}$ (at 90 percent confidence). Neither black holes nor neutron stars are expected to form below ~ 1 solar mass through conventional stellar evolution, though it has been proposed that similarly low mass black holes could be formed primordially through density fluctuations in the early universe. Under a particular primordial black hole binary formation scenario, we constrain monochromatic primordial black hole populations of 0.2 $M_\odot$ to be less than $33\%$ of the total dark matter density and monochromatic populations of 1.0 $M_\odot$ to be less than $5\%$ of the dark matter density. The latter strengthens the presently placed bounds from micro-lensing surveys of MAssive Compact Halo Objects (MACHOs) provided by the MACHO and EROS collaborations. △ Less

Submitted 15 August, 2018; v1 submitted 14 August, 2018; originally announced August 2018.

Comments: 12 pages, 3 figures

Report number: LIGO-DCC-P1800158-v12

Journal ref: Phys. Rev. Lett. 121, 231103 (2018)

A standard siren measurement of the Hubble constant from GW170817 without the electromagnetic counterpart

Authors: M. Fishbach, R. Gray, I. Magaña Hernandez, H. Qi, A. Sur, F. Acernese, L. Aiello, A. Allocca, M. A. Aloy, A. Amato, S. Antier, M. Arène, N. Arnaud, S. Ascenzi, P. Astone, F. Aubin, S. Babak, P. Bacon, F. Badaracco, M. K. M. Bader, F. Baldaccini, G. Ballardin, F. Barone, M. Barsuglia, D. Barta , et al. (294 additional authors not shown)

Abstract: We perform a statistical standard siren analysis of GW170817. Our analysis does not utilize knowledge of NGC 4993 as the unique host galaxy of the optical counterpart to GW170817. Instead, we consider each galaxy within the GW170817 localization region as a potential host; combining the redshift from each galaxy with the distance estimate from GW170817 provides an estimate of the Hubble constant,… ▽ More We perform a statistical standard siren analysis of GW170817. Our analysis does not utilize knowledge of NGC 4993 as the unique host galaxy of the optical counterpart to GW170817. Instead, we consider each galaxy within the GW170817 localization region as a potential host; combining the redshift from each galaxy with the distance estimate from GW170817 provides an estimate of the Hubble constant, $H_0$. We then combine the $H_0$ values from all the galaxies to provide a final measurement of $H_0$. We explore the dependence of our results on the thresholds by which galaxies are included in our sample, as well as the impact of weighting the galaxies by stellar mass and star-formation rate. Considering all galaxies brighter than $0.01 L^\star_B$ as equally likely to host a BNS merger, we find $H_0= 76^{+48}_{-23}$ km s$^{-1}$ Mpc$^{-1}$ (maximum a posteriori and 68.3% highest density posterior interval; assuming a flat $H_0$ prior in the range $\left[ 10, 220 \right]$ km s$^{-1}$ Mpc$^{-1}$). Restricting only to galaxies brighter than $0.626 L^\star_B$ tightens the measurement to $H_0= 77^{+37}_{-18}$ km s$^{-1}$ Mpc$^{-1}$. We show that weighting the host galaxies by stellar mass or star-formation rate provides entirely consistent results with potentially tighter constraints. While these statistical estimates are inferior to the value from the counterpart standard siren measurement utilizing NGC 4993 as the unique host, $H_0=76^{+19}_{-13}$ km s$^{-1}$ Mpc$^{-1}$ (determined from the same publicly available data), our analysis is a proof-of-principle demonstration of the statistical approach first proposed by Bernard Schutz over 30 years ago. △ Less

Submitted 5 January, 2019; v1 submitted 15 July, 2018; originally announced July 2018.

Comments: 14 pages, 7 figures. matches ApJL version (in press). addition of Figs. 1 and 7 to address uniqueness of GW170817 for statistical H0 measurements

Report number: LIGO-P1800192

Calibration of Advanced Virgo and Reconstruction of the Gravitational Wave Signal h(t) during the Observing Run O2

Authors: Virgo Collaboration, F. Acernese, T. Adams, K. Agatsuma, L. Aiello, A. Allocca, M. A. Aloy, A. Amato, S. Antier, M. Arène, N. Arnaud, S. Ascenzi, P. Astone, F. Aubin, S. Babak, P. Bacon, F. Badaracco, M. K. M. Bader, F. Baldaccini, G. Ballardin, F. Barone, M. Barsuglia, D. Barta, A. Basti, M. Bawaj , et al. (259 additional authors not shown)

Abstract: In August 2017, Advanced Virgo joined Advanced LIGO for the end of the O2 run, leading to the first gravitational waves detections with the three-detector network. This paper describes the Advanced Virgo calibration and the gravitational wave strain h(t) reconstruction during O2. The methods are the same as the ones developed for the initial Virgo detector and have already been described in previo… ▽ More In August 2017, Advanced Virgo joined Advanced LIGO for the end of the O2 run, leading to the first gravitational waves detections with the three-detector network. This paper describes the Advanced Virgo calibration and the gravitational wave strain h(t) reconstruction during O2. The methods are the same as the ones developed for the initial Virgo detector and have already been described in previous publications, this paper summarizes the differences and emphasis is put on estimating systematic uncertainties. Three versions of the h(t) signal have been computed for the Virgo O2 run, an online version and two post-run reprocessed versions with improved detector calibration and reconstruction algorithm. A photon calibrator has been used to establish the sign of h(t) and to make an independent partial cross-check of the systematic uncertainties. The uncertainties reached for the latest h(t) version are 5.1% in amplitude, 40 mrad in phase and 20 microseconds in timing. △ Less

Submitted 3 September, 2018; v1 submitted 9 July, 2018; originally announced July 2018.

Comments: 20 pages, 6 figures. Submitted to Classical and Quantum Gravity. Version updated after referees' comments

A dust-enshrouded tidal disruption event with a resolved radio jet in a galaxy merger

Authors: S. Mattila, M. Pérez-Torres, A. Efstathiou, P. Mimica, M. Fraser, E. Kankare, A. Alberdi, M. Á. Aloy, T. Heikkilä, P. G. Jonker, P. Lundqvist, I. Martí-Vidal, W. P. S. Meikle, C. Romero-Cañizales, S. J. Smartt, S. Tsygankov, E. Varenius, A. Alonso-Herrero, M. Bondi, C. Fransson, R. Herrero-Illana, T. Kangas, R. Kotak, N. Ramírez-Olivencia, P. Väisänen , et al. (11 additional authors not shown)

Abstract: Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its… ▽ More Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its emission re-radiated at infrared wavelengths by dust. Efficient reprocessing by dense gas and dust may explain the difference between theoretical predictions and observed luminosities of TDEs. The radio observations resolve an expanding and decelerating jet, probing the jet formation and evolution around a SMBH. △ Less

Submitted 14 June, 2018; originally announced June 2018.

Comments: Published online by Science on 14 Jun 2018

Black holes, gravitational waves and fundamental physics: a roadmap

Authors: Leor Barack, Vitor Cardoso, Samaya Nissanke, Thomas P. Sotiriou, Abbas Askar, Krzysztof Belczynski, Gianfranco Bertone, Edi Bon, Diego Blas, Richard Brito, Tomasz Bulik, Clare Burrage, Christian T. Byrnes, Chiara Caprini, Masha Chernyakova, Piotr Chrusciel, Monica Colpi, Valeria Ferrari, Daniele Gaggero, Jonathan Gair, Juan Garcia-Bellido, S. F. Hassan, Lavinia Heisenberg, Martin Hendry, Ik Siong Heng , et al. (181 additional authors not shown)

Abstract: The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horiz… ▽ More The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress. △ Less

Submitted 1 February, 2019; v1 submitted 13 June, 2018; originally announced June 2018.

Comments: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

Neutron star collapse and gravitational waves with a non-convex equation of state

Authors: M. A. Aloy, J. M. Ibáñez, N. Sanchis-Gual, M. Obergaulinger, J. A. Font, S. Serna, A. Marquina

Abstract: The thermodynamical properties of the equation of state (EoS) of high-density matter (above nuclear saturation density) and the possible existence of exotic states such as phase transitions from nuclear/hadronic matter into quark-gluon plasma, or the appearance of hyperons, may critically influence the stability and dynamics of compact relativistic stars. From a theoretical point of view, establis… ▽ More The thermodynamical properties of the equation of state (EoS) of high-density matter (above nuclear saturation density) and the possible existence of exotic states such as phase transitions from nuclear/hadronic matter into quark-gluon plasma, or the appearance of hyperons, may critically influence the stability and dynamics of compact relativistic stars. From a theoretical point of view, establishing the existence of those states requires the analysis of the `convexity' of the EoS. We show indications of the existence of regions in the dense-matter EoS where the thermodynamics may be non-convex as a result of a non-monotonic dependence of the sound speed with the rest-mass density. When this happens, non-conventional dynamics may develop. In this paper we investigate the effects of a phenomenological, non-convex EoS on the equilibrium structure of stable compact stars and on the dynamics of unstable neutron stars that collapse gravitationally to black holes, both for spherically symmetric and uniformly-rotating configurations. We show how the dynamics of the collapse with a non-convex EoS departs from the convex case, leaving distinctive imprints on the gravitational waveforms. The astrophysical significance of these results for microphysical EoSs is discussed. △ Less

Submitted 24 January, 2019; v1 submitted 8 June, 2018; originally announced June 2018.

Comments: 29 pages, 22 figures, Accepted by MNRAS on January 24tth 2019. The author order has changed with respect to the previous arXiv version

Core collapse with magnetic fields and rotation

Authors: M. Obergaulinger, O. Just, M. Á. Aloy

Abstract: We study the effects of magnetic fields and rotation on the core collapse of a star of an initial mass of M = 20 solar masses using axisymmetric simulations coupling special relativistic magnetohydrodynamics, an approximately relativistic gravitational potential, and spectral neutrino transport. We compare models of the same core with different, artificially added profiles of rotation and magnetic… ▽ More We study the effects of magnetic fields and rotation on the core collapse of a star of an initial mass of M = 20 solar masses using axisymmetric simulations coupling special relativistic magnetohydrodynamics, an approximately relativistic gravitational potential, and spectral neutrino transport. We compare models of the same core with different, artificially added profiles of rotation and magnetic field. A model with weak field and slow rotation does not produce an explosion, while stronger fields and fast rotation open the possibility of explosions. Whereas the neutrino luminosities of the exploding models are the same as or even less than those of the non-exploding model, magnetic fields locally in equipartition with the gas pressure provide a strong contribution to the shock revival and the acceleration of bipolar outflows. Among the amplification processes generating such strong fields, we find the magneto-rotational instability. However, our limited grid resolution allows us to find it only in limited regions of the model with the strongest pre-collapse field ($10^{11}$ G) and fastest rotation. △ Less

Submitted 1 June, 2018; originally announced June 2018.

Comments: 18 pages, 9 figures. J. Phys. G: Nucl. Part. Phys. in press

Journal ref: Journal of Physics G: Nuclear and Particle Physics, Volume 45, Number 8. Focus on microphysics in core-collapse supernovae: 30 years since SN1987A

GW170817: Measurements of Neutron Star Radii and Equation of State

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. (1127 additional authors not shown)

Abstract: On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of th… ▽ More On 17 August 2017, the LIGO and Virgo observatories made the first direct detection of gravitational waves from the coalescence of a neutron star binary system. The detection of this gravitational-wave signal, GW170817, offers a novel opportunity to directly probe the properties of matter at the extreme conditions found in the interior of these stars. The initial, minimal-assumption analysis of the LIGO and Virgo data placed constraints on the tidal effects of the coalescing bodies, which were then translated to constraints on neutron star radii. Here, we expand upon previous analyses by working under the hypothesis that both bodies were neutron stars that are described by the same equation of state and have spins within the range observed in Galactic binary neutron stars. Our analysis employs two methods: the use of equation-of-state-insensitive relations between various macroscopic properties of the neutron stars and the use of an efficient parametrization of the defining function $p(ρ)$ of the equation of state itself. From the LIGO and Virgo data alone and the first method, we measure the two neutron star radii as $R_1=10.8^{+2.0}_{-1.7}$ km for the heavier star and $R_2= 10.7^{+2.1}_{-1.5}$ km for the lighter star at the 90% credible level. If we additionally require that the equation of state supports neutron stars with masses larger than $1.97 \,M_\odot$ as required from electromagnetic observations and employ the equation-of-state parametrization, we further constrain $R_1= 11.9^{+1.4}_{-1.4}$ km and $R_2= 11.9^{+1.4}_{-1.4}$ km at the 90% credible level. Finally, we obtain constraints on $p(ρ)$ at supranuclear densities, with pressure at twice nuclear saturation density measured at $3.5^{+2.7}_{-1.7}\times 10^{34} \,\mathrm{dyn}/\mathrm{cm}^{2}$ at the 90% level. △ Less

Submitted 15 October, 2018; v1 submitted 29 May, 2018; originally announced May 2018.

Comments: 10 pages, 3 figures; v2 matches published version; data associated with the figures can be found at https://dcc.ligo.org/LIGO-P1800115/public

Report number: LIGO-P1800115

Journal ref: Phys. Rev. Lett. 121, 161101 (2018)

Properties 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. (1126 additional authors not shown)

Abstract: On August 17, 2017, the Advanced LIGO and Advanced Virgo gravitational-wave detectors observed a low-mass compact binary inspiral. The initial sky localization of the source of the gravitational-wave signal, GW170817, allowed electromagnetic observatories to identify NGC 4993 as the host galaxy. In this work, we improve initial estimates of the binary's properties, including component masses, spin… ▽ More On August 17, 2017, the Advanced LIGO and Advanced Virgo gravitational-wave detectors observed a low-mass compact binary inspiral. The initial sky localization of the source of the gravitational-wave signal, GW170817, allowed electromagnetic observatories to identify NGC 4993 as the host galaxy. In this work, we improve initial estimates of the binary's properties, including component masses, spins, and tidal parameters, using the known source location, improved modeling, and recalibrated Virgo data. We extend the range of gravitational-wave frequencies considered down to 23 Hz, compared to 30 Hz in the initial analysis. We also compare results inferred using several signal models, which are more accurate and incorporate additional physical effects as compared to the initial analysis. We improve the localization of the gravitational-wave source to a 90% credible region of $16~\mathrm{deg}^2$. We find tighter constraints on the masses, spins, and tidal parameters, and continue to find no evidence for nonzero component spins. The component masses are inferred to lie between 1.00 and 1.89 $M_\odot$ when allowing for large component spins, and to lie between 1.16 and 1.60 $M_\odot$ (with a total mass $2.73^{+0.04}_{-0.01} \, M_\odot$) when the spins are restricted to be within the range observed in Galactic binary neutron stars. Under minimal assumptions about the nature of the compact objects, our constraints for the tidal deformability parameter $\tilde Λ$ are $(0,630)$ when we allow for large component spins, and $300^{+420}_{-230}$ (using a 90% highest posterior density interval) when restricting the magnitude of the component spins, ruling out several equation-of-state models at the 90% credible level. Finally, with LIGO and GEO600 data, we use a Bayesian analysis to place upper limits on the amplitude and spectral energy density of a possible post-merger signal. (Abridged) △ Less

Submitted 22 January, 2019; v1 submitted 29 May, 2018; originally announced May 2018.

Comments: 32 pages, 15 figures. Version 3 is the final published version; results are unchanged compared to previous versions. Data behind the figures, including posterior samples, are available at dcc.ligo.org/LIGO-P1800061/public. The gravitational wave strain data for this event are available at the LIGO Open Science Center, gw-openscience.org/events/GW170817

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

Numerically solving the relativistic Grad-Shafranov equation in Kerr spacetimes: Numerical techniques

Authors: J. F. Mahlmann, P. Cerdá-Durán, M. A. Aloy

Abstract: The study of the electrodynamics of static, axisymmetric and force-free Kerr magnetospheres relies vastly on solutions of the so called relativistic Grad-Shafranov equation (GSE). Different numerical approaches to the solution of the GSE have been introduced in the literature, but none of them has been fully assessed from the numerical point of view in terms of efficiency and quality of the soluti… ▽ More The study of the electrodynamics of static, axisymmetric and force-free Kerr magnetospheres relies vastly on solutions of the so called relativistic Grad-Shafranov equation (GSE). Different numerical approaches to the solution of the GSE have been introduced in the literature, but none of them has been fully assessed from the numerical point of view in terms of efficiency and quality of the solutions found. We present a generalization of these algorithms and give detailed background on the algorithmic implementation. We assess the numerical stability of the implemented algorithms and quantify the convergence of the presented methodology for the most established setups (split-monopole, paraboloidal, BH-disk, uniform). △ Less

Submitted 3 May, 2018; v1 submitted 2 February, 2018; originally announced February 2018.

Comments: 19 pages, 12 figure of which 9 are in color. Accepted for publication in MNRAS

Anomalous dynamics triggered by a non-convex equation of state in relativistic flows

Authors: J. M. Ibáñez, A. Marquina, S. Serna, M. A. Aloy

Abstract: The non-monotonicity of the local speed of sound in dense matter at baryon number densities much higher than the nuclear saturation density ($n_0 \approx 0.16\,$fm$^{-3}$) suggests the possible existence of a non-convex thermodynamics which will lead to a non-convex dynamics. Here, we explore the rich and complex dynamics that an equation of state (EoS) with non-convex regions in the pressure-dens… ▽ More The non-monotonicity of the local speed of sound in dense matter at baryon number densities much higher than the nuclear saturation density ($n_0 \approx 0.16\,$fm$^{-3}$) suggests the possible existence of a non-convex thermodynamics which will lead to a non-convex dynamics. Here, we explore the rich and complex dynamics that an equation of state (EoS) with non-convex regions in the pressure-density plane may develop as a result of genuinely relativistic effects, without a classical counterpart. To this end, we have introduced a phenomenological EoS, whose parameters can be restricted heeding to causality and thermodynamic stability constraints. This EoS shall be regarded as a toy-model with which we may mimic realistic (and far more complex) EoS of practical use in the realm of Relativistic Hydrodynamics. △ Less

Submitted 7 January, 2018; v1 submitted 8 December, 2017; originally announced December 2017.

Comments: 22 pages, 8 figures. Accepted by MNRAS

An HLLC Riemann Solver for Resistive Relativistic Magnetohydrodynamics

Authors: S. Miranda-Aranguren, M. A. Aloy, T. Rembiasz

Abstract: We present a new approximate Riemann solver for the augmented system of equations of resistive relativistic magnetohydrodynamics (RRMHD) that belongs to the family of Harten-Lax-van Leer contact wave (HLLC) solvers. In HLLC solvers, the solution is approximated by two constant states flanked by two shocks separated by a contact wave. The accuracy of the new approximate solver is calibrated through… ▽ More We present a new approximate Riemann solver for the augmented system of equations of resistive relativistic magnetohydrodynamics (RRMHD) that belongs to the family of Harten-Lax-van Leer contact wave (HLLC) solvers. In HLLC solvers, the solution is approximated by two constant states flanked by two shocks separated by a contact wave. The accuracy of the new approximate solver is calibrated through one- and two-dimensional test problems. △ Less

Submitted 30 May, 2018; v1 submitted 17 November, 2017; originally announced November 2017.

Comments: 26 pages, 20 figures, 3 tables, submitted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 476, Issue 3, 21 May 2018, Pages 3837-3860

Linear Theory of the Rayleigh-Taylor Instability at a Discontinuous Surface of a Relativistic Flow

Authors: Jin Matsumoto, Miguel A. Aloy, Manel Perucho

Abstract: We address the linear stability of a discontinuous surface of a relativistic flow in the context of a jet that oscillates radially as it propagates. The restoring force of the oscillation is expected to drive a Rayleigh-Taylor instability (RTI) at the interface between the jet and its cocoon. We perform a linear analysis and numerical simulations of the growth of the RTI in the transverse plane to… ▽ More We address the linear stability of a discontinuous surface of a relativistic flow in the context of a jet that oscillates radially as it propagates. The restoring force of the oscillation is expected to drive a Rayleigh-Taylor instability (RTI) at the interface between the jet and its cocoon. We perform a linear analysis and numerical simulations of the growth of the RTI in the transverse plane to the jet flow with a uniform acceleration. In this system, an inertia force due to the uniform acceleration acts as the restoring force for the oscillation. We find that not only the difference in the inertia between the two fluids separated by the interface but also the pressure at the interface helps to drive the RTI because of a difference in the Lorenz factor across the discontinuous surface of the jet. The dispersion relation indicates that the linear growth rate of each mode becomes maximum when the Lorentz factor of the jet is much larger than that of the cocoon and the pressure at the jet interface is relativistic. By comparing the linear growth rates of the RTI in the analytical model and the numerical simulations, the validity of our analytically derived dispersion relation for the relativistic RTI is confirmed. △ Less

Submitted 2 August, 2017; v1 submitted 15 July, 2017; originally announced July 2017.

Comments: 11 pages, 6 figures. Accepted for publication in MNRAS