The black hole low-mass X-ray binary V404 Cygni is part of a wide triple

Evidence suggests that, when compact objects such as black holes and neutron stars form, they may receive a 'natal kick', during which the stellar remnant gains momentum. Observational evidence for neutron star kicks is substantial1,2, yet is limited for black hole natal kicks, and some proposed black hole formation scenarios result in very small kicks3-5. Here we report that the canonical black hole low-mass X-ray binary (LMXB) V404 Cygni is part of a wide hierarchical triple with a tertiary companion at least 3,500 astronomical units (AU) away from the inner binary. Given the orbital configuration, the black hole probably received a sub-5 km s-1 kick to have avoided unbinding the tertiary. This discovery lends support to the idea that at least some black holes form with nearly no natal kick. Furthermore, the tertiary in this system lends credence to evolutionary models of LMXBs involving a hierarchical triple structure6. Remarkably, the tertiary is evolved, indicating that the system formed 3-5 billion years ago and that the black hole has removed at least half a solar mass of matter from its evolved secondary companion. During the event in which the black hole formed, it is required that at least half of the mass of the black hole progenitor collapsed into the black hole; it may even have undergone a complete implosion, enabling the tertiary to remain loosely bound.

Probing the Pulsar Explanation of the Galactic-Center GeV Excess Using Continuous Gravitational-Wave Searches

Over 10 years ago, Fermi observed an excess of GeV gamma rays from the Galactic Center whose origin is still under debate. One explanation for this excess involves annihilating dark matter, another requires an unresolved population of millisecond pulsars concentrated at the Galactic Center. In this work, we use the results from LIGO and Virgo's most recent all-sky search for quasimonochromatic, persistent gravitational-wave signals from isolated neutron stars, which is estimated to be about 20%-50% of the population, to determine whether unresolved millisecond pulsars could actually explain this excess. First, we choose a luminosity function that determines the number of millisecond pulsars required to explain the observed excess. Then, we consider two models for deformations on millisecond pulsars to determine their ellipticity distributions, which are directly related to their gravitational-wave radiation. Lastly, based on null results from the O3 frequency-Hough all-sky search for continuous gravitational waves, we find that a large set of the parameter space in the pulsar luminosity function can be excluded. We also evaluate how these exclusion regions may change with respect to various model choices. Our results are the first of their kind and represent a bridge between gamma-ray astrophysics, gravitational-wave astronomy, and dark-matter physics.

Impact of Dynamical Tides on the Reconstruction of the Neutron Star Equation of State

Gravitational waves (GWs) from inspiraling neutron stars afford us a unique opportunity to infer the as-of-yet unknown equation of state of cold hadronic matter at supranuclear densities. During the inspiral, the dominant matter effects arise due to the star's response to their companion's tidal field, leaving a characteristic imprint in the emitted GW signal. This unique signature allows us to constrain the cold neutron star equation of state. At GW frequencies above ≳800  Hz, however, subdominant tidal effects known as dynamical tides become important. In this Letter, we demonstrate that neglecting dynamical tidal effects associated with the fundamental (f) mode leads to large systematic biases in the measured tidal deformability of the stars and hence in the inferred neutron star equation of state. Importantly, we find that f-mode dynamical tides will already be relevant for Advanced LIGO's and Virgo's fifth observing run (∼2025)-neglecting dynamical tides can lead to errors on the neutron radius of O(1  km), with dramatic implications for the measurement of the equation of state. Our results demonstrate that the accurate modeling of subdominant tidal effects beyond the adiabatic limit will be crucial to perform accurate measurements of the neutron star equation of state in upcoming GW observations.

A transient radio source consistent with a merger-triggered core collapse supernova

[Figure: see text].

Was GW190814 a Black Hole-Strange Quark Star System?

We investigate the possibility that the low mass companion of the black hole in the source of GW190814 was a strange quark star. This possibility is viable within the so-called two-families scenario in which neutron stars and strange quark stars coexist. Strange quark stars can reach the mass range indicated by GW190814, M∼(2.5-2.67)  M_{⊙} due to a large value of the adiabatic index, without the need for a velocity of sound close to the causal limit. Neutron stars (actually hyperonic stars in the two-families scenario) can instead fulfill the presently available astrophysical and nuclear physics constraints which require a softer equation of state. In this scheme it is possible to satisfy both the request of very large stellar masses and of small radii while using totally realistic and physically motivated equations of state. Moreover it is possible to get a radius for a 1.4  M_{⊙} star of the order or less than 11 km, which is impossible if only one family of compact stars exists.

Multimessenger Binary Mergers Containing Neutron Stars: Gravitational Waves, Jets, and γ-Ray Bursts

Neutron stars (NSs) are extraordinary not only because they are the densest form of matter in the visible Universe but also because they can generate magnetic fields ten orders of magnitude larger than those currently constructed on earth. The combination of extreme gravity with the enormous electromagnetic (EM) fields gives rise to spectacular phenomena like those observed on August 2017 with the merger of a binary neutron star system, an event that generated a gravitational wave (GW) signal, a short γ -ray burst (sGRB), and a kilonova. This event serves as the highlight so far of the era of multimessenger astronomy. In this review, we present the current state of our theoretical understanding of compact binary mergers containing NSs as gleaned from the latest general relativistic magnetohydrodynamic simulations. Such mergers can lead to events like the one on August 2017, GW170817, and its EM counterparts, GRB 170817 and AT 2017gfo. In addition to exploring the GW emission from binary black hole-neutron star and neutron star-neutron star mergers, we also focus on their counterpart EM signals. In particular, we are interested in identifying the conditions under which a relativistic jet can be launched following these mergers. Such a jet is an essential feature of most sGRB models and provides the main conduit of energy from the central object to the outer radiation regions. Jet properties, including their lifetimes and Poynting luminosities, the effects of the initial magnetic field geometries and spins of the coalescing NSs, as well as their governing equation of state, are discussed. Lastly, we present our current understanding of how the Blandford-Znajek mechanism arises from merger remnants as the trigger for launching jets, if, when and how a horizon is necessary for this mechanism, and the possibility that it can turn on in magnetized neutron ergostars, which contain ergoregions, but no horizons.

Black hole-neutron star coalescence: Effects of the neutron star spin on jet launching and dynamical ejecta mass

Black hole-neutron star (BHNS) mergers are thought to be sources of gravitational waves (GWs) with coincident electromagnetic (EM) counterparts. To further probe whether these systems are viable progenitors of short gamma-ray bursts (SGRBs) and kilonovas, and how one may use (the lack of) EM counterparts associated with LIGO/Virgo candidate BHNS GW events to sharpen parameter estimation, we study the impact of neutron star spin in BHNS mergers. Using dynamical spacetime magnetohydrodynamic simulations of BHNSs initially on a quasicircular orbit, we survey configurations that differ in the BH spin (a BH/M BH = 0 and 0.75), the NS spin (a NS/M NS = -0.17, 0, 0.23, and 0.33), and the binary mass ratio (q = M BH:M NS = 3:1 and 5:1). The general trend we find is that increasing the NS prograde spin increases both the rest mass of the accretion disk onto the remnant black hole, and the rest mass of dynamically ejected matter. By a time Δt ~ 3500-5500M ~ 88-138(M NS/1.4 M ⊙) ms after the peak gravitational-wave amplitude, a magnetically driven jet is launched only for q = 3:1 regardless of the initial NS spin. The lifetime of the jets [Δt ~ 0.5-0.8(M NS/1.4 M ⊙) s] and their outgoing Poynting luminosity [L Poyn ~ 1051.5±0.5 erg/s] are consistent with typical SGRBs' luminosities and expectations from the Blandford-Znajek mechanism. By the time we terminate our simulations, we do not observe either an outflow or a large-scale magnetic-field collimation for the other systems we consider. The mass range of dynamically ejected matter is 10-4.5-10-2(M NS/1.4 M ⊙) M ⊙, which can power kilonovas with peak bolometric luminosities L knova ~ 1040-1041.4 erg/s with rise times ≲6.5 h and potentially detectable by the LSST.

Asymmetric mass ratios for bright double neutron-star mergers

The discovery of a radioactively powered kilonova associated with the binary neutron-star merger GW170817 remains the only confirmed electromagnetic counterpart to a gravitational-wave event^1,2. Observations of the late-time electromagnetic emission, however, do not agree with the expectations from standard neutron-star merger models. Although the large measured ejecta mass^3,4 could be explained by a progenitor system that is asymmetric in terms of the stellar component masses (that is, with a mass ratio q of 0.7 to 0.8)⁵, the known Galactic population of merging double neutron-star systems (that is, those that will coalesce within billions of years or less) has until now consisted only of nearly equal-mass (q > 0.9) binaries⁶. The pulsar PSR J1913+1102 is a double system in a five-hour, low-eccentricity (0.09) orbit, with an orbital separation of 1.8 solar radii⁷, and the two neutron stars are predicted to coalesce in [Formula: see text] million years owing to gravitational-wave emission. Here we report that the masses of the pulsar and the companion neutron star, as measured by a dedicated pulsar timing campaign, are 1.62 ± 0.03 and 1.27 ± 0.03 solar masses, respectively. With a measured mass ratio of q = 0.78 ± 0.03, this is the most asymmetric merging system reported so far. On the basis of this detection, our population synthesis analysis implies that such asymmetric binaries represent between 2 and 30 per cent (90 per cent confidence) of the total population of merging binaries. The coalescence of a member of this population offers a possible explanation for the anomalous properties of GW170817, including the observed kilonova emission from that event.

On the origin of GW190425

The LIGO/Virgo collaborations recently announced the detection of a binary neutron star merger, GW190425. The mass of GW190425 is significantly larger than the masses of Galactic double neutron stars known through radio astronomy. We hypothesize that GW190425 formed differently from Galactic double neutron stars, via unstable ‘case BB’ mass transfer. According to this hypothesis, the progenitor of GW190425 was a binary consisting of a neutron star and a ∼4–$5\, {\mathrm{ M}_\odot }$ helium star, which underwent common-envelope evolution. Following the supernova of the helium star, an eccentric double neutron star was formed, which merged in ${\lesssim }10\, {\rm Myr}$. The helium star progenitor may explain the unusually large mass of GW190425, while the short time to merger may explain why similar systems are not observed in radio. To test this hypothesis, we measure the eccentricity of GW190425 using publicly available LIGO/Virgo data. We constrain the eccentricity at $10\, {\rm Hz}$ to be e ≤ 0.007 with $90{{\ \rm per\ cent}}$ confidence. This provides no evidence for or against the unstable mass transfer scenario, because the binary is likely to have circularized to e ≲ 10−4 by the time it was detected. Future detectors will help to reveal the formation channel of mergers similar to GW190425 using eccentricity measurements.

Lense-Thirring frame dragging induced by a fast-rotating white dwarf in a binary pulsar system

Radio pulsars in short-period eccentric binary orbits can be used to study both gravitational dynamics and binary evolution. The binary system containing PSR J1141-6545 includes a massive white dwarf (WD) companion that formed before the gravitationally bound young radio pulsar. We observed a temporal evolution of the orbital inclination of this pulsar that we infer is caused by a combination of a Newtonian quadrupole moment and Lense-Thirring (LT) precession of the orbit resulting from rapid rotation of the WD. LT precession, an effect of relativistic frame dragging, is a prediction of general relativity. This detection is consistent with an evolutionary scenario in which the WD accreted matter from the pulsar progenitor, spinning up the WD to a period of <200 seconds.

Observational diversity of magnetized neutron stars

Young and rotation-powered neutron stars (NSs) are commonly observed as rapidly-spinning pulsars. They dissipate their rotational energy by emitting pulsar wind with electromagnetic radiation and spin down at a steady rate, according to the simple steadily-rotating magnetic dipole model. In reality, however, multiwavelength observations of radiation from the NS surface and magnetosphere have revealed that the evolution and properties of NSs are highly diverse, often dubbed as 'NS zoo'. In particular, many of young and highly magnetized NSs show a high degree of activities, such as sporadic electromagnetic outbursts and irregular changes in pulse arrival times. Importantly, their magnetic field, which are the strongest in the universe, makes them ideal laboratories for fundamental physics. A class of highly-magnetized isolated NSs is empirically divided into several subclasses. In a broad classification, they are, in the order of the magnetic field strength (B) from the highest, 'magnetars' (historically recognized as soft gamma-ray repeaters and/or anomalous x-ray pulsars), 'high-B pulsars', and (nearby) x-ray isolated NSs. This article presents an introductory review for non-astrophysicists about the observational properties of highly-magnetized NSs, and their implications. The observed dynamic nature of NSs must be interpreted in conjunction with transient magnetic activities triggered during magnetic-energy dissipation process. In particular, we focus on how the five fundamental quantities of NSs, i.e. mass, radius, spin period, surface temperature, and magnetic fields, as observed with modern instruments, change with evolution of, and vary depending on the class of, the NSs. They are the foundation for a future unified theory of NSs.

Questions Related to the Equation of State of High-Density Matter

Astronomical data about neutron stars can be combined with laboratory nuclear data to give us a strong base from which to infer the equation of state of cold catalyzed matter beyond nuclear density. However, the nuclear and astrophysical communities are largely distinct; each has their own methods, which means that there is often imperfect communication between the communities regarding caveats about claimed measurements and constraints. Here we present a brief summary from one astronomer’s perspective of relevant observations of neutron stars, with warnings as appropriate, followed by a set of questions that are intended to help enhance the dialog between nuclear physicists and astrophysicists.

GW170817: Measurements of Neutron Star Radii and Equation of State

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_{-1.7}^{+2.0}  km for the heavier star and R_{2}=10.7_{-1.5}^{+2.1}  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_{⊙} 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_{-1.7}^{+2.7}×10^{34}  dyn cm^{-2} at the 90% level.

A hot and fast ultra-stripped supernova that likely formed a compact neutron star binary

Compact neutron star binary systems are produced from binary massive stars through stellar evolution involving up to two supernova explosions. The final stages in the formation of these systems have not been directly observed. We report the discovery of iPTF 14gqr (SN 2014ft), a type Ic supernova with a fast-evolving light curve indicating an extremely low ejecta mass (≈0.2 solar masses) and low kinetic energy (≈2 × 10⁵⁰ ergs). Early photometry and spectroscopy reveal evidence of shock cooling of an extended helium-rich envelope, likely ejected in an intense pre-explosion mass-loss episode of the progenitor. Taken together, we interpret iPTF 14gqr as evidence for ultra-stripped supernovae that form neutron stars in compact binary systems.

Tidal Deformabilities and Radii of Neutron Stars from the Observation of GW170817

We use gravitational-wave observations of the binary neutron star merger GW170817 to explore the tidal deformabilities and radii of neutron stars. We perform a Bayesian parameter estimation with the source location and distance informed by electromagnetic observations. We also assume that the two stars have the same equation of state; we demonstrate that, for stars with masses comparable to the component masses of GW170817, this is effectively implemented by assuming that the stars' dimensionless tidal deformabilities are determined by the binary's mass ratio q by Λ_{1}/Λ_{2}=q^{6}. We investigate different choices of prior on the component masses of the neutron stars. We find that the tidal deformability and 90% credible interval is Λ[over ˜]=222_{-138}^{+420} for a uniform component mass prior, Λ[over ˜]=245_{-151}^{+453} for a component mass prior informed by radio observations of Galactic double neutron stars, and Λ[over ˜]=233_{-144}^{+448} for a component mass prior informed by radio pulsars. We find a robust measurement of the common areal radius of the neutron stars across all mass priors of 8.9≤R[over ^]≤13.2  km, with a mean value of ⟨R[over ^]⟩=10.8  km. Our results are the first measurement of tidal deformability with a physical constraint on the star's equation of state and place the first lower bounds on the deformability and areal radii of neutron stars using gravitational waves.

GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral

On August 17, 2017 at 12∶41:04 UTC the Advanced LIGO and Advanced Virgo gravitational-wave detectors made their first observation of a binary neutron star inspiral. The signal, GW170817, was detected with a combined signal-to-noise ratio of 32.4 and a false-alarm-rate estimate of less than one per 8.0×10^{4}  years. We infer the component masses of the binary to be between 0.86 and 2.26  M_{⊙}, in agreement with masses of known neutron stars. Restricting the component spins to the range inferred in binary neutron stars, we find the component masses to be in the range 1.17-1.60  M_{⊙}, with the total mass of the system 2.74_{-0.01}^{+0.04}M_{⊙}. The source was localized within a sky region of 28  deg^{2} (90% probability) and had a luminosity distance of 40_{-14}^{+8}  Mpc, the closest and most precisely localized gravitational-wave signal yet. The association with the γ-ray burst GRB 170817A, detected by Fermi-GBM 1.7 s after the coalescence, corroborates the hypothesis of a neutron star merger and provides the first direct evidence of a link between these mergers and short γ-ray bursts. Subsequent identification of transient counterparts across the electromagnetic spectrum in the same location further supports the interpretation of this event as a neutron star merger. This unprecedented joint gravitational and electromagnetic observation provides insight into astrophysics, dense matter, gravitation, and cosmology.