Interruption of Tidal Disruption Flares By SMBHBs

Interruption of Tidal Disruption Flares By Supermassive Black Hole Binaries

F. K. Liu11affiliation: Astronomy Department, Peking University, 100871 Beijing, China; fkliu@bac.pku.edu.cn, lis@bac.pku.edu.cn, chenx@bac.pku.edu.cn 22affiliation: Kavli Institute for Astronomy and Astrophysics, Peking University, 100871 Beijing, China , S. Li11affiliation: Astronomy Department, Peking University, 100871 Beijing, China; fkliu@bac.pku.edu.cn, lis@bac.pku.edu.cn, chenx@bac.pku.edu.cn , and Xian Chen11affiliation: Astronomy Department, Peking University, 100871 Beijing, China; fkliu@bac.pku.edu.cn, lis@bac.pku.edu.cn, chenx@bac.pku.edu.cn 33affiliation: Department of Astronomy, University of California, Santa Cruz, CA 96064, USA
Abstract

Supermassive black hole binaries (SMBHBs) are products of galaxy mergers, and are important in testing cold dark matter cosmology and locating gravitational-wave-radiation sources. A unique electromagnetic signature of SMBHBs in galactic nuclei is essential in identifying the binaries in observations from the IR band through optical to X-ray. Recently, the flares in optical, UV, and X-ray caused by supermassive black holes (SMBHs) tidally disrupting nearby stars have been successfully used to observationally probe single SMBHs in normal galaxies. In this Letter, we investigate the accretion of the gaseous debris of a tidally disrupted star by a SMBHB. Using both stability analysis of three-body systems and numerical scattering experiments, we show that the accretion of stellar debris gas, which initially decays with time , would stop at a time . Here, and is the orbital period of the SMBHB. After a period of interruption, the accretion recurs discretely at time , where . Both and sensitively depend on the orbital parameters of the tidally disrupted star at the tidal radius and the orbit eccentricity of SMBHB. The interrupted accretion of the stellar debris gas gives rise to an interrupted tidal flare, which could be used to identify SMBHBs in non-active galaxies in the upcoming transient surveys.

accretion, accretion disks – black hole physics – galaxies: active – galaxies: evolution – galaxies: nuclei – gravitational waves

1 Introduction

Supermassive black hole binaries (SMBHBs) are predicted by the hierarchical galaxy formation model in cold dark matter (CDM) cosmology (Begelman et al., 1980; Volonteri et al., 2003). After a SMBHB at the center of merging systems become hard, it would stall at the hard radius for a timescale even longer than the Hubble time if spherical two-body relaxation dominates (Begelman et al., 1980). However, recent investigations suggested that the hardening rates of SMBHBs can be boosted and SMBHBs may coalesce within a Hubble time either due to various stellar dynamical processes other than the spherical two-body relaxation (Yu, 2002; Chatterjee et al., 2003; Merritt & Poon, 2004; Berczik et al., 2006; Sesana et al., 2008), or due to gas dynamics (Gould & Rix, 2000; Liu et al., 2003; Colpi & Dotti, 2009, and references therein).

The strong gravitational wave (GW) radiation generated by coalescing supermassive black holes (SMBHs) is the main target of the GW detector the Laser Interferometer Space Antenna (LISA) and of the Pulsar Timing Array (PTA) program. Because of the poor accuracy of both LISA and PTA in locating GW radiation sources, it is of key importance to detect electromagnetic counterparts (EMCs) of GW radiation sources. Identifying SMBHBs by their EMCs is also essential to constraining the poorly understood galaxy-merger history. Several EMCs have been suggested in the literature to probe SMBHBs and their coalescence: (1) precession of jet orientation and its acceleration in radio galaxies during the in-spiraling of SMBHBs (Begelman et al., 1980; Liu & Chen, 2007), (2) optical periodic outbursts in blazars due to the interaction between SMBHB and accretion disk (Sillanpaa et al., 1988; Liu et al., 1995, 2006; Liu & Wu, 2002; Valtonen et al., 2008; Haiman et al., 2009), (3) jet reorientation in X-shaped radio galaxies due to the exchange of angular momentum between SMBHB and accretion disk (Liu, 2004), (4) two systems of broad emission lines (BELs) in quasars (Boroson & Lauer, 2009), (5) intermittent activity in double-double radio galaxies at binary coalescence (Liu et al., 2003), (6) X-ray, UV, optical, and IR afterglow following binary coalescence (Milosavljević & Phinney, 2005; Shields & Bonning, 2008; Lippai et al., 2008; Schnittman & Krolik, 2008), and (7) systematically shifted BELs relative to narrow emission lines (Merritt et al., 2006; Komossa et al., 2008) and off-center active galactic nuclei (AGNs; Madau & Quataert, 2004; Loeb, 2007) because of SMBH GW radiation recoil.

All the above observational signatures require gas accretion disks around SMBHBs. The SMBHBs in gas-poor galactic nuclei are difficult to detect, because the SMBHs are dormant. However, a dormant SMBH could be temporarily activated by tidally disrupting a star passing by and accreting the disrupted stellar debris (Hills, 1975). A tidal flare decays typically as a power law (Rees, 1988; Phinney, 1989), which has been observed in several non-active galaxies (Komossa & Bade, 1999; Komossa, 2002; Halpern et al., 2004; Esquej et al., 2008; Gezari et al., 2008, 2009). Recently, Chen et al. (2008) and Chen et al. (2009) calculated the tidal disruption rate in SMBHB systems at different evolutionary stages, and found that it is significantly different from the typical rate for single SMBHs by orders of magnitude. This great difference of the flaring rates enables one to statistically constrain the SMBHB population in normal galaxies (Chen et al., 2008), but identifying SMBHBs individually is still difficult at present. In this Letter, we investigate the influence of a SMBHB on the accretion of tidally disrupted stellar plasma and on the tidal flare. We show that the accretion of the stellar debris is interrupted by the SMBHB and this interruption can be taken as a key distinguishable observational signature for SMBHBs in gas-poor galactic nuclei.

2 Tidal Disruption in a SMBHB system

A star with mass and radius will be tidally disrupted by a SMBH with mass if it approaches the black hole (BH) within the tidal radius

(1)

(Hills, 1975; Rees, 1988; Phinney, 1989), where is a dimensionless parameter of order unity. We focus on tidal disruptions of solar type stars with solar radius and solar mass by SMBHs with , which are detected in the tidal-flare surveys (Komossa & Bade, 1999; Komossa, 2002; Esquej et al., 2008; Gezari et al., 2008, 2009). For more massive BHs, the tidal radii become smaller than the event horizon (Ivanov & Chernyakova, 2006). After tidal disruption at the tidal radius, the specific energy across the star ranges from to , where is the orbital binding energy of the star, is the spread in specific energy across the stellar radius, and if the tidal spin-up of the star is negligible or if the star is spun-up to the break-up angular velocity (Rees, 1988; Lacy et al., 1982; Li et al., 2002). Since is times larger than the internal binding energy of the star and the interaction between the stellar debris is negligible after the tidal disruption (Evans & Kochanek, 1989), each fluid element moves ballistically.

If we neglect the influence of the companion SMBH, the bound debris with should move in Keplerian orbits with eccentricity of about unity and fall back to the tidal radius after one Keplerian period, . The debris continuously returning to the tidal radius generates a falling-back rate . For Keplerian orbits, . If we assume following a constant distribution in the energy range to (see e.g. Rees, 1988; Evans & Kochanek, 1989; Lodato et al., 2009), the falling-back rate evolves as

(2)

for , where is the disruption time and

(3)

is the returning time for the most bound debris (typically ). It is assumed that once the bound material comes back to the tidal disruption radius, it loses kinetic energy due to strong shocks because of the interaction between fluid elements and is circularized on a timescale much shorter than to form an orbiting torus at a radius about around the SMBH (e.g. Rees, 1988; Phinney, 1989; Ulmer, 1999; Li et al., 2002). Both the radiative dissipation of the shock energy and the accretion of the gas torus onto the BH will give rise to an X-ray/UV flare, decaying with .

For a SMBHB system, the stellar-disrupting BH could be either the primary (with mass ) or the secondary (with mass ). Since the probability of stellar disruption by the secondary is relatively low if the two BHs are very unequal (Chen et al., 2008, 2009), here we assume the primary BH to be the stellar-disrupting one. Because the interaction between the fluid elements is negligible before they come back to the tidal disruption radius, the SMBHB and each bound fluid element constitute a restricted three-body system. When the orbit of a bound fluid element is inside the SMBHB orbit, the system is called S-type. An S-type three-body system could be stable for a long time only when the system is hierarchical, namely, the system consists of an inner binary (the bound fluid element and the primary BH) on a nearly Keplerian orbit with semi-major axis and eccentricity of near unit, and an outer binary with eccentricity and semi-major axis in which the secondary BH orbits the mass center of the inner binary. The orbital change of a stable S-type system is negligible on the fluid-element dynamical timescale, so the falling-back stellar debris may interact with one another at the tidal disruption radius and finally be accreted with accretion rate given in Equation (2). However, for those less-bound fluid elements with larger than a critical radius , the orbit of the triple system becomes chaotic. For a triple system with an angle between the angular momenta of the inner and outer binaries, can be semiempirically given by

(4)

(Mardling & Aarseth, 2001), where is the BH mass ratio.

Because of the nonlinear overlap of the multiple resonances (Mardling, 2007) in the chaotic triple systems, the fluid elements significantly exchange angular momentum and energy with the SMBHB and change their orbit dramatically on the dynamical timescale of the triple system, therefore would not return to the tidal radius to fuel the accreting torus and to form continuous accretion. Although a fraction of the fluid elements with chaotic orbits may return to the tidal radius on a timescale much longer than its Keplerian period, others would escape the system during the three-body interactions (see our numerical simulations in Section 3). Therefore, Equation (4) implies that the continuous accretion of tidally disrupted stellar plasma stops with the accretion of the last fluid element with , this occurring at a time

(5)

where is the SMBHB orbital period. The eccentricity of hard SMBHBs is moderate (Milosavljević & Merritt, 2001) and minor mergers with are the most common in the hierarchical galaxy formation model (Volonteri et al., 2003). We re-write Equation (5) as , where is in the range for and , and independent of . For a SMBHB with semi-major axis residing in a stellar cluster with velocity dispersion , where

(6)

is the hard radius and is the hardness, the accretion is interrupted at

(7)

where , and .

3 Numerical simulations and results

Because the fluid elements of a tidally disrupted star move like test particles before falling back to the tidal radius, their evolution could be correctly simulated with scattering experiments. We use particles to resolve the stellar debris of a disrupted star and integrate their trajectories in restricted three-body systems. The particles are logarithmically sampled in the binding-energy range , where is a binding energy corresponding to an orbital period of and is the spread in specific energy for a non-spinning star. For a typical tidal-disruption event with , the th particle has a binding energy , pericenter velocity , and mass , where is the size of the th energy bin.

The evolution of the three-body systems is computed in a frame centered on the mass center of the SMBHB with plane aligned with the SMBHB orbital plane. We assume in this work for simplicity, and that the secondary BH initially lies on the positive axis and moves in the direction of positive . In each experiment, a test particle starts with velocity from the pericenter at tidal radius about the primary BH, and the orbit is determined by three initial parameters (determined by the orbital parameters of the disrupted star at tidal radius) (1) the inclination angle between the orbital planes of the particle and the SMBHB, (2) the longitude of ascending node , and (3) the argument of pericenter . With the initial parameters, we integrate the equations of motion in the pseudo-Newtonian potential (Paczyński & Wiita, 1980), where is the Schwarzschild radius of BHs, using an explicit Runge-Kutta method of order 8 (Hairer et al., 1987). Stellar debris which falls back to within from the primary BH is assumed to be circularized via shocks and be accreted instantaneously. Therefore, if a particle with id passes by the primary BH within a distance 111We test our simulations with , , and and the results are nearly independent of ., the integration is stopped and the time is recorded. Otherwise, the integration continues until time . The accretion rate of the stellar debris is calculated using the recorded and .

Figure 1: Accretion rates of gaseous debris in units of solar mass per year vs. time. Time starts at tidal disruption and is in unit of the binary orbital period . The solid lines are the simulation results and the dashed lines give the standard in single BH systems. marks the interruption time estimated with Equation (5). The upper, middle, and bottom panels are, respectively, for , , and .

Figure 1 shows the accretion rate of stellar debris in our fiducial simulations with , , , , , , , and . Our results suggest that the interruption of accretion occurs at , depending on . In Figure 1, a better agreement between the numerical value and the analytical estimate for and is by coincidence, because not only depends on but also on and . To illustrate this, we did 100 numerical simulations with random , , and . The results show that ranges from to with a mean value , consistent with the analytical mean value 0.27 given by Equation (5). Figure 1 also shows that after being interrupted for about , accretion recurs and “accretion islands” emerge discretely. The accretion rate at the islands is variable and can be larger than the corresponding value for single BH. The duration of interruption and the periods of the accretion islands decrease with , the shortest occurring at .

Figure 2: Accretion rates of gaseous debris vs. time for different at . The red and green solid lines correspond to and , respectively. The other parameters are the same as those in Figure 1.

To investigate the effect of different on the results, we ran a scattering experiment for and , keeping the other parameters as in Figure 1. The result is given in Figure 2. Figure 2 shows that the interruption occurs at a slightly earlier time as hardness decreases, which seems inconsistent with the analytical prediction given by Equation (5). However, we should note that Equation (5) is obtained statistically by averaging over different and . When we did the simulations with different and , the averaged interruption time for different resides in the predicted range. Accordingly, in the scattering experiments with different BH mass () and mass ratios (), we find that weakly depends on and . Our numerical results are consistent with the analytical estimate given by Equation (5), implying that the dependence of on , , and is much weaker than on , , and .

4 Discussions

We investigated the accretion of tidally disrupted stellar debris in SMBHB systems. For simplicity, we assumed for the tidal debris gas: (1) a constant distribution of mass in binding energy at the tidal radius, (2) ballistic motion and negligible interaction of the fluid elements after tidal disruption, and (3) instant circularization of the debris gas probably due to shocks because of interaction between the fluid elements returned to within a radius two times the tidal radius. The first two assumptions have been justified by numerical hydrodynamic simulations (e.g. Evans & Kochanek, 1989; Lodato et al., 2009), but the third one needs to be verified by hydrodynamic simulations capable of capturing strong shocks. With these three assumptions, the fluid elements and SMBHB compose restricted three-body systems, so we can investigate the accretion of the gaseous stellar debris by analyzing the stability of the three-body systems using the resonance overlap stability criterion (Mardling, 2007). Because of the chaotic nature induced by the nonlinear overlap of several orbital mean motion resonances, the fluid elements inside the chaotic regions significantly change their orbits on a dynamical timescale and do not return to the tidal radius to fuel the BH, leading to the interruption of accretion. We also investigated the evolution of the stellar debris using three-body scattering experiments. Our results obtained both analytically and numerically show that the accretion rate of the debris gas decreases with a power law until a critical time . At time , the accretion pauses until about one SMBHB orbital period . relates to with , where is typically and in the range , depending on the initial orbital parameters () of the fluid elements at tidal radius and on the eccentricity of the SMBHB, but being nearly independent of the SMBHB semimajor , BH masses, and mass ratio . This suspension of accretion would result in an interruption of the tidal flare, although the residual accretion disk may still radiate weak optical, UV, and X-ray emission during the interruption.

Our numerical results indicate that the accretion of the gaseous stellar debris restarts at a time , leading to a flicker of flare. The exact recurring time, , depends on the initial orbital parameters , , and , but our numerical results suggest that is of order unity and . The interruption timescale of the tidal flare in a SMBHB system is

(8)

Our numerical simulations suggest that the accreted plasma during the discrete accretion consists of both the tidal debris gas falling back to the tidal radius after one Keplerian time and a fraction of those fluid elements with chaotic orbits falling back to the accreting torus on a timescale longer than its Keplerian time. When and is determined observationally, constraints on and could be made if we take .

In our simulations, we assume that the orbital binding energy of the tidally disrupted star is negligible compared to the spread in specific energy . For such tidal-disruption events, Equations (3) and (7) imply that the standard power-law decay and the interruption of tidal flares are detectable if

(9)

If the tidally disrupted star is initially very bound to one of the binary BHs so that is comparable to (Chen et al., 2009), the time would be much shorter and the interruption of tidal flare is detectable even in an ultra-hard SMBHB with . For a SMBHB with , the semi-major axis of the binary is about . The GW radiation emitted by such SMBHBs could be detected by PTA.

An interrupted tidal flare could be caught if is shorter than the mission duration of an transient survey , which corresponds to SMBHBs with . Here we use the relation (Tremaine et al., 2002). Because a hard SMBHB with spends most of its life time at (Yu, 2002; Sesana et al., 2008), we have a good chance to detect them with upcoming transient surveys. However, if one wants to catch an interrupted tidal flare from SMBHBs emitting strong GW radiation (), the time resolution of the survey should be . Because a SMBHB spends a small fraction of its lifetime at this stage, high-sensitivity and deep transient surveys are needed to accumulate many more tidal-disruption events.

We are grateful to X.-B. Wu, D.N.C. Lin, S. Aarseth, M. Valtonen and S. Komossa for useful discussions. Many thanks are due to the referee, R.A. Mardling, for valuable comments. Thanks are given to S. Justham for careful read of the manuscript. This work is supported by NSFC (10573001) and Chinese National 973 program (2007CB815405). XC thanks the China Scholarship Council for financial support. The numerical computation was carried on the SGI Altix 330 system at the Astronomy Department of PKU.

References

  • Begelman et al. (1980) Begelman, M. C., Blandford, R. D., & Rees, M. J. 1980, Nature, 287, 307
  • Berczik et al. (2006) Berczik, P., Merritt, D., Spurzem, R., & Bischof, H.-P. 2006, ApJ, 642, L21
  • Boroson & Lauer (2009) Boroson, T. A., & Lauer, T. R. 2009, Nature, 458, 53
  • Chatterjee et al. (2003) Chatterjee, P., Hernquist, L., & Loeb, A. 2003, ApJ, 592, 32
  • Chen et al. (2008) Chen, X., Liu, F. K., & Magorrian, J. 2008, ApJ, 676, 54
  • Chen et al. (2009) Chen, X., Madau, P., Sesana, A., & Liu, F. K. 2009, ApJ, 697, L149
  • Colpi & Dotti (2009) Colpi, M., & Dotti, M. 2009, arXiv:0906.4339 [astro-ph]
  • Esquej et al. (2008) Esquej, P., et al. 2008, A&A, 489, 543
  • Evans & Kochanek (1989) Evans, C.R., & Kochanek, C. S. 1989, ApJ, 346, L13
  • Gezari et al. (2008) Gezari, S., et al. 2008, ApJ, 676, 944
  • Gezari et al. (2009) Gezari, S., et al. 2009, ApJ, 698, 1367
  • Gould & Rix (2000) Gould, A., & Rix, H.-W. 2000, ApJ, 532, L29
  • Haiman et al. (2009) Haiman, Z., Kocsis, B., & Menou, K. 2009, ApJ, 700, 1952
  • Hairer et al. (1987) Hairer, E., Norsett, S. P., & Wanner, G. 1987, Solving Ordinary Differential Equations I (1st ed.; Berlin: Springer)
  • Halpern et al. (2004) Halpern, J. P., Gezari, S., & Komossa, S. 2004, ApJ, 604, 572
  • Hills (1975) Hills, J. G. 1975, Nature, 254, 295
  • Ivanov & Chernyakova (2006) Ivanov, P. B., & Chernyakova, M. A. 2006, A&A, 448, 843
  • Komossa (2002) Komossa, S. 2002, Rev. Mod. Astron., 15, 27
  • Komossa & Bade (1999) Komossa, S., & Bade, N. 1999, A&A, 343, 775
  • Komossa et al. (2008) Komossa, S., Zhou, H., & Lu, H. 2008, ApJ, 678, L81
  • Lacy et al. (1982) Lacy, J. H., Townes, C. H., & Hollenbach, D. J. 1982, ApJ, 262, 120
  • Li et al. (2002) Li, L.-X., Narayan, R., & Menou, K. 2002, ApJ, 576, 753
  • Lippai et al. (2008) Lippai, Z., Frei, Z., & Haiman, Z. 2008, ApJ, 676, L5
  • Liu (2004) Liu, F. K. 2004, MNRAS, 347, 1357
  • Liu & Chen (2007) Liu, F. K., & Chen, X. 2007, ApJ, 671, 1272
  • Liu & Wu (2002) Liu, F. K., & Wu, X.-B. 2002, A&A, 388, L48
  • Liu et al. (2003) Liu, F. K., Wu, X.-B., & Cao, S. L. 2003, MNRAS, 340, 411
  • Liu et al. (1995) Liu, F. K., Xie, G. Z., & Bai, J. M. 1995, A&A, 295, 1
  • Liu et al. (2006) Liu, F. K., Zhao, G., & Wu X.-B. 2006, ApJ, 650, 749
  • Lodato et al. (2009) Lodato, G., King, A. R., & Pringle, J. E. 2009, MNRAS, 392, 332
  • Loeb (2007) Loeb, A. 2007, Phys. Rev. Lett., 99, 041103
  • Madau & Quataert (2004) Madau, P., & Quataert, E. 2004, ApJ, 606, L17
  • Mardling (2007) Mardling, R. A. 2007, in IAU Symp. 246, Dynamical Evolution of Dense Stellar Systems, (eds.) E. Vesperini, M. Giersz & A. Sills, (Cambridge Univ. Press: Cambridge), p.199
  • Mardling & Aarseth (2001) Mardling, R. A., & Aarseth, S. J. 2001, MNRAS, 321, 398
  • Merritt & Poon (2004) Merritt, D., & Poon, M. Y. 2004, ApJ, 606, 788
  • Merritt et al. (2006) Merritt, D., et al. 2006, MNRAS, 367, 1746
  • Milosavljević & Merritt (2001) Milosavljević, M., & Merritt, D., 2001, ApJ, 563, 34
  • Milosavljević & Phinney (2005) Milosavljević, M., & Phinney, E. S. 2005, ApJ, 622, L93
  • Paczyński & Wiita (1980) Paczyński, B., & Wiita, P. J. 1980, A&A, 88, 23
  • Phinney (1989) Phinney, E. S. 1989, in IAU Symp. 136, The Center of the Galaxy, (ed.) M. Morris (Dordrecht: Kluwer), p.543
  • Rees (1988) Rees, M. J. 1988, Nature, 333, 523
  • Schnittman & Krolik (2008) Schnittman, J. D., & Krolik, J. H. 2008, ApJ, 684, 835
  • Sesana et al. (2008) Sesana, A., Haardt, F., & Madau, P. 2008, ApJ, 686, 432
  • Shields & Bonning (2008) Shields, G. A., & Bonning, E. W. 2008, ApJ, 682, 758
  • Sillanpaa et al. (1988) Sillanpaa, A., Haarala, S., Valtonen, M. J., Sundelius, B., & Byrd, G. G. 1988, ApJ, 325, 628
  • Tremaine et al. (2002) Tremaine, S., et al. 2002, ApJ, 574, 740
  • Ulmer (1999) Ulmer, A. 1999, ApJ, 514, 180
  • Valtonen et al. (2008) Valtonen, M. J., et al. 2008, Nature, 452, 851
  • Volonteri et al. (2003) Volonteri, M., Haardt, F., & Madau, P. 2003, ApJ, 582, 559
  • Yu (2002) Yu, Q. 2002, MNRAS, 331, 935
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
236993
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description