skip to main content

Title: Electromagnetic transients and gravitational waves from white dwarf disruptions by stellar black holes in triple systems
ABSTRACT Mergers of binaries comprising compact objects can give rise to explosive transient events, heralding the birth of exotic objects that cannot be formed through single-star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effects of gravitational wave (GW) emission and Lidov–Kozai oscillations imposed by the tidal field of an outer tertiary companion orbiting the inner BH–WD binary. We explore the sensitivity of our results to the distributions of natal kick velocities imparted to the BH and WD upon formation, adiabatic mass loss, semimajor axes and eccentricities of the triples, and stellar-mass ratios. We find rates of WD–tidal disruption events (TDEs) in the range 1.2 × 10−3 − 1.4 Gpc−3 yr−1 for z ≤ 0.1, rarer than stellar TDEs in triples by a factor of ∼3–30. The uncertainty in the TDE rates may be greatly reduced in the future using GW observations of Galactic binaries and triples with LISA. WD–TDEs may give rise to high-energy X-ray or gamma-ray transients of duration similar to long gamma-ray bursts but lacking the signatures of a core-collapse supernova, while being accompanied by a supernova-like optical transient that lasts for only days. WD–BH and WD–NS binaries will also emit GWs in the LISA band before the TDE. The discovery and identification of triple-induced WD–TDE events by future time domain surveys and/or GWs could enable the study of the demographics of BHs in nearby galaxies.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
1061 to 1072
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    Many astrophysical environments, from star clusters and globular clusters to the discs of active galactic nuclei, are characterized by frequent interactions between stars and the compact objects that they leave behind. Here, using a suite of 3D hydrodynamics simulations, we explore the outcome of close interactions between $1\, \mathrm{M}_{\odot }$ stars and binary black holes (BBHs) in the gravitational wave regime, resulting in a tidal disruption event (TDE) or a pure scattering, focusing on the accretion rates, the back reaction on the BH binary orbital parameters, and the increase in the binary BH effective spin. We find that TDEs can make a significant impact on the binary orbit, which is often different from that of a pure scattering. Binaries experiencing a prograde (retrograde) TDE tend to be widened (hardened) by up to $\simeq 20{{\ \rm per\ cent}}$. Initially circular binaries become more eccentric by $\lesssim 10{{\ \rm per\ cent}}$ by a prograde or retrograde TDE, whereas the eccentricity of initially eccentric binaries increases (decreases) by a retrograde (prograde) TDE by $\lesssim 5{{\ \rm per\ cent}}$. Overall, a single TDE can generally result in changes of the gravitational-wave-driven merger time-scale by order unity. The accretion rates of both black holes are very highly super-Eddington, showing modulations (preferentially for retrograde TDEs) on a time-scale of the orbital period, which can be a characteristic feature of BBH-driven TDEs. Prograde TDEs result in the effective spin parameter χ to vary by ≲0.02, while χ ≳ −0.005 for retrograde TDEs.

    more » « less
  2. Aims. The modelling of spectroscopic observations of tidal disruption events (TDEs) to date suggests that the newly formed accretion disks are mostly quasi-circular. In this work we study the transient event AT 2020zso, hosted by an active galactic nucleus (AGN; as inferred from narrow emission line diagnostics), with the aim of characterising the properties of its newly formed accretion flow. Methods. We classify AT 2020zso as a TDE based on the blackbody evolution inferred from UV/optical photometric observations and spectral line content and evolution. We identify transient, double-peaked Bowen (N  III ), He  I , He  II, and H α emission lines. We model medium-resolution optical spectroscopy of the He  II (after careful de-blending of the N  III contribution) and H α lines during the rise, peak, and early decline of the light curve using relativistic, elliptical accretion disk models. Results. We find that the spectral evolution before the peak can be explained by optical depth effects consistent with an outflowing, optically thick Eddington envelope. Around the peak, the envelope reaches its maximum extent (approximately 10 15 cm, or ∼3000–6000 gravitational radii for an inferred black hole mass of 5−10 × 10 5 M ⊙ ) and becomes optically thin. The H α and He  II emission lines at and after the peak can be reproduced with a highly inclined ( i  = 85 ± 5 degrees), highly elliptical ( e  = 0.97 ± 0.01), and relatively compact ( R in = several 100 R g and R out = several 1000 R g ) accretion disk. Conclusions. Overall, the line profiles suggest a highly elliptical geometry for the new accretion flow, consistent with theoretical expectations of newly formed TDE disks. We quantitatively confirm, for the first time, the high inclination nature of a Bowen (and X-ray dim) TDE, consistent with the unification picture of TDEs, where the inclination largely determines the observational appearance. Rapid line profile variations rule out the binary supermassive black hole hypothesis as the origin of the eccentricity; these results thus provide a direct link between a TDE in an AGN and the eccentric accretion disk. We illustrate for the first time how optical spectroscopy can be used to constrain the black hole spin, through (the lack of) disk precession signatures (changes in inferred inclination). We constrain the disk alignment timescale to > 15 days in AT2020zso, which rules out high black hole spin values ( a  < 0.8) for M BH  ∼ 10 6 M ⊙ and disk viscosity α  ≳ 0.1. 
    more » « less
  3. Abstract

    We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility over 3 yr (from 2018 October to 2021 September). We infer the black hole (BH) mass (MBH) with host galaxy scaling relations, showing that the sampleMBHranges from 105.1Mto 108.2M. We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frameg-band luminosity function can be well described by a broken power law ofϕ(Lg)Lg/Lbk0.3+Lg/Lbk2.61, withLbk= 1043.1erg s−1. In the BH mass regime of 105.3≲ (MBH/M) ≲ 107.3, the TDE mass function followsϕ(MBH)MBH0.25, which favors a flat local BH mass function (dnBH/dlogMBHconstant). We confirm the significant rate suppression at the high-mass end (MBH≳ 107.5M), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass ofMgal∼ 1010M, the average optical TDE rate is ≈3.2 × 10−5galaxy−1yr−1. We constrain the optical TDE rate to be [3.7, 7.4, and 1.6] × 10−5galaxy−1yr−1in galaxies with red, green, and blue colors.

    more » « less

    The launch of space-based gravitational-wave (GW) detectors (e.g. Laser Interferometry Space Antenna; LISA) and current and upcoming Pulsar Timing Arrays will extend the GW window to low frequencies, opening new investigations into dynamical processes involving massive black hole binaries (MBHBs) and their mergers across cosmic time. MBHBs are expected to be among the primary sources for the upcoming low-frequency (10−4–10−1 Hz) window probed by LISA. It is important to investigate the expected supermassive BH merger rates and associated signals, to determine how potential LISA events are affected by physics included in current models. To study this, we post-process the large population of MBHBs in the Illustris simulation to account for dynamical friction time delays associated with BH infall/inspiral. We show that merger delays associated with binary evolution have the potential to decrease the expected merger rates, with $M_{\rm {BH}}\ \gt\ 10^6\ \mathrm{M}_\odot$ MBHBs (the lowest mass in Illustris) decreasing from ∼3 to ∼0.1 yr−1, and shifting the merger peak from z ∼2 to ∼1.25. During this time, we estimate that accretion grows the total merging mass by as much as 7x the original mass. Importantly, however, dynamical friction-associated delays (which shift the mergers toward lower redshift and higher masses) lead to a stronger signal/strain for the emitted GWs in the LISA band, increasing mean frequency from 10−3.1 to 10−3.4–10−4.0 Hz, and mean strain from 10−17.2 to 10−16.3–10−15.3. Finally, we show that after including a merger delay and associated MBH growth, mergers still tend to lie on the typical MBH–M* relation, but with an increased likelihood of an undermassive BH.

    more » « less

    When a star passes close to a supermassive black hole (BH), the BH’s tidal forces rip it apart into a thin stream, leading to a tidal disruption event (TDE). In this work, we study the post-disruption phase of TDEs in general relativistic hydrodynamics (GRHD) using our GPU-accelerated code h-amr. We carry out the first grid-based simulation of a deep-penetration TDE (β = 7) with realistic system parameters: a black hole-to-star mass ratio of 106, a parabolic stellar trajectory, and a non-zero BH spin. We also carry out a simulation of a tilted TDE whose stellar orbit is inclined relative to the BH midplane. We show that for our aligned TDE, an accretion disc forms due to the dissipation of orbital energy with ∼20 per cent of the infalling material reaching the BH. The dissipation is initially dominated by violent self-intersections and later by stream–disc interactions near the pericentre. The self-intersections completely disrupt the incoming stream, resulting in five distinct self-intersection events separated by approximately 12 h and a flaring in the accretion rate. We also find that the disc is eccentric with mean eccentricity e ≈ 0.88. For our tilted TDE, we find only partial self-intersections due to nodal precession near pericentre. Although these partial intersections eject gas out of the orbital plane, an accretion disc still forms with a similar accreted fraction of the material to the aligned case. These results have important implications for disc formation in realistic tidal disruptions. For instance, the periodicity in accretion rate induced by the complete stream disruption may explain the flaring events from Swift J1644+57.

    more » « less