skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Close encounters of stars with stellar-mass black hole binaries
ABSTRACT 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
Award ID(s):
2006839
PAR ID:
10371036
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
516
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
p. 2204-2217
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT 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
  2. ABSTRACT We investigate the effect of the cutting-edge circumbinary disc (CBD) evolution models on massive black hole binary (MBHB) populations and the gravitational wave background (GWB). We show that CBD-driven evolution leaves a tell-tale signature in MBHB populations, by driving binaries towards an equilibrium eccentricity that depends on the binary mass ratio. We find high orbital eccentricities ($$e_{\rm b} \sim 0.5$$) as MBHBs enter multimessenger observable frequency bands. The CBD-induced eccentricity distribution of MBHB populations in observable bands is independent of the initial eccentricity distribution at binary formation, erasing any memory of eccentricities induced in the large-scale dynamics of merging galaxies. Our results suggest that eccentric MBHBs are the rule rather than the exception in upcoming transient surveys, provided that CBDs regularly form in MBHB systems. We show that the GWB amplitude is sensitive to CBD-driven preferential accretion onto the secondary, resulting in an increase in GWB amplitude $$A_{\rm yr^{-1}}$$ by over 100 per cent with just 10 per cent Eddington accretion. As we self-consistently allow for binary hardening and softening, we show that CBD-driven orbital expansion does not diminish the GWB amplitude, and instead increases the amplitude by a small amount. We further present detection rates and population statistics of MBHBs with $$M_{\rm b} \gtrsim 10^6 \, {\rm M}_{\odot }$$ in Laser Interferometer Space Antenna, showing that most binaries have equal mass ratios and can retain residual eccentricities up to $$e_{\rm b} \sim 10^{-3}$$ due to CBD-driven evolution. 
    more » « less
  3. ABSTRACT We assess the possibility of detecting both eccentricity and gas effects (migration and accretion) in the gravitational wave (GW) signal from LISA massive black hole binaries at redshift $z=1$. Gas induces a phase correction to the GW signal with an effective amplitude ($$C_{\rm g}$$) and a semimajor axis dependence (assumed to follow a power-law with slope $$n_{\rm g}$$). We use a complete model of the LISA response and employ a gas-corrected post-Newtonian inspiral-only waveform model TaylorF2Ecc. By using the Fisher formalism and Bayesian inference, we constrain $$C_{\rm g}$$ together with the initial eccentricity $$e_0$$, the total redshifted mass $$M_z$$, the primary-to-secondary mass ratio q, the dimensionless spins $$\chi _{1,2}$$ of both component BHs, and the time of coalescence $$t_c$$. We find that simultaneously constraining $$C_{\rm g}$$ and $$e_0$$ leads to worse constraints on both parameters with respect to when considered individually. For a standard thin viscous accretion disc around $$M_z=10^5~{\rm M}_{\odot }$$, $q=8$, $$\chi _{1,2}=0.9$$, and $$t_c=4$$ years MBHB, we can confidently measure (with a relative error of $$\lt 50$$ per cent) an Eddington ratio $${\rm f}_{\rm Edd}\sim 0.1$$ for a circular binary and $${\rm f}_{\rm Edd}\sim 1$$ for an eccentric system assuming $$\mathcal {O}(10)$$ stronger gas torque near-merger than at the currently explored much-wider binary separations. The minimum measurable eccentricity is $$e_0\gtrsim 10^{-2.75}$$ in vacuum and $$e_0\gtrsim 10^{-2}$$ in gas. A weak environmental perturbation ($${\rm f}_{\rm Edd}\lesssim 1$$) to a circular binary can be mimicked by an orbital eccentricity during inspiral, implying that an electromagnetic counterpart would be required to confirm the presence of an accretion disc. 
    more » « less
  4. null (Ed.)
    ABSTRACT In this work, we combine spectroscopic information from the SkyMapper survey for Extremely Metal-Poor stars and astrometry from Gaia DR2 to investigate the kinematics of a sample of 475 stars with a metallicity range of $$-6.5 \le \rm [Fe/H] \le -2.05$$ dex. Exploiting the action map, we identify 16 and 40 stars dynamically consistent with the Gaia Sausage and Gaia Sequoia accretion events, respectively. The most metal poor of these candidates have metallicities of $$\rm [Fe/H]=-3.31\, \mathrm{ and }\, -3.74$$, respectively, helping to define the low-metallicity tail of the progenitors involved in the accretion events. We also find, consistent with other studies, that ∼21 per cent of the sample have orbits that remain confined to within 3 kpc of the Galactic plane, that is, |Zmax| ≤ 3 kpc. Of particular interest is a subsample (∼11 per cent of the total) of low |Zmax| stars with low eccentricities and prograde motions. The lowest metallicity of these stars has [Fe/H] = –4.30 and the subsample is best interpreted as the very low-metallicity tail of the metal-weak thick disc population. The low |Zmax|, low eccentricity stars with retrograde orbits are likely accreted, while the low |Zmax|, high eccentricity pro- and retrograde stars are plausibly associated with the Gaia Sausage system. We find that a small fraction of our sample (∼4 per cent of the total) is likely escaping from the Galaxy, and postulate that these stars have gained energy from gravitational interactions that occur when infalling dwarf galaxies are tidally disrupted. 
    more » « less
  5. ABSTRACT At least $$70\, {\rm per\, cent}$$ of massive OBA-type stars reside in binary or higher order systems. The dynamical evolution of these systems can lend insight into the origins of extreme phenomena such as X-ray binaries and gravitational wave sources. In one such dynamical process, the Eccentric Kozai–Lidov (EKL) mechanism, a third companion star alters the secular evolution of a binary system. For dynamical stability, these triple systems must have a hierarchical configuration. We explore the effects of a distant third companion’s gravitational perturbations on a massive binary’s orbital configuration before significant stellar evolution has taken place (≤10 Myr). We include tidal dissipation and general relativistic precession. With large (38 000 total) Monte Carlo realizations of massive hierarchical triples, we characterize imprints of the birth conditions on the final orbital distributions. Specifically, we find that the final eccentricity distribution over the range of 0.1–0.7 is an excellent indicator of its birth distribution. Furthermore, we find that the period distributions have a similar mapping for wide orbits. Finally, we demonstrate that the observed period distribution for approximately 10-Myr-old massive stars is consistent with EKL evolution. 
    more » « less