Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Planetary nebula wikipedia , lookup
Standard solar model wikipedia , lookup
Astrophysical X-ray source wikipedia , lookup
White dwarf wikipedia , lookup
Astronomical spectroscopy wikipedia , lookup
Kerr metric wikipedia , lookup
Hawking radiation wikipedia , lookup
Hayashi track wikipedia , lookup
First observation of gravitational waves wikipedia , lookup
Main sequence wikipedia , lookup
Stellar evolution wikipedia , lookup
MNRAS 434, 2948–2960 (2013) doi:10.1093/mnras/stt1209 Advance Access publication 2013 August 8 Roche accretion of stars close to massive black holes Lixin Dai‹ and Roger Blandford‹ Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, Menlo Park, CA 94025, USA Accepted 2013 July 1. Received 2013 May 21; in original form 2012 September 10 ABSTRACT In this paper, we consider Roche accretion in an extreme mass-ratio inspiral binary system formed by a star orbiting a massive black hole. The ultimate goal is to constrain the mass and spin of the black hole and provide a test of general relativity in the strong-field regime from the resultant quasi-periodic signals. Before accretion starts, the stellar orbit is presumed to be circular and equatorial, and shrinks due to gravitational radiation. New fitting formulae are presented for the inspiral time and the radiation reaction torque in the relativistic regime. If the inspiralling star fills its Roche lobe outside the innermost stable circular orbit of the hole, gas will flow through the inner Lagrange point (L1) to the hole. We give a new, accurate interpolation formula of the volume enclosed by the relativistic Roche lobe. If this mass transfer happens on a time-scale faster than the thermal time-scale but slower than the dynamical timescale, the star will evolve adiabatically, and, in most cases, will recede from the hole while filling its Roche lobe. We calculate how the stellar orbital period and mass transfer rate change through the ‘Roche evolution’ for various types of stars in the relativistic regime. We envisage that the mass stream eventually hits the accretion disc, where it forms a hotspot orbiting the hole and may ultimately modulate the luminosity with the stellar orbital frequency. The observability of such a modulation is discussed along with a possible interpretation of an intermittent 1 h period in the X-ray emission of RE J1034+396. Key words: accretion, accretion discs – radiative transfer – stars: winds, outflows – galaxies: active – galaxies: individual: RE J1034+396 – X-rays: binaries. 1 I N T RO D U C T I O N The measurement of black hole parameters – mass and spin – has been an important astrophysics research topic over the past decade. Such measurements can lead to tests of general relativity in the strong-field regime, and reveal information on the formation and evolution of galaxies. There has been great progress on this topic using various approaches. The mass of a black hole is more easily measured since it is a Newtonian parameter. However, the spin, a relativistic parameter, is much harder to be measured. One method is to use radiation from the inner region of the accretion disc: either by fitting the shifts and shapes of the iron lines (e.g. Reynolds & Nowak 2003) or by using a multitemperature blackbody model (Yuan, Quataert & Narayan 2004; Li, Narayan & McClintock 2009; Shcherbakov, Penna & McKinney 2010; Gou et al. 2011). The problem with the iron line model is the uncertainty of whether the iron lines indeed come from near the innermost stable circular orbit (ISCO). The continuum method can be employed to obtain the spin of stellar mass black holes effectively. However, for supermassive black holes (SMBHs), the radiation from the inner disc region is in the UV range and can be easily absorbed. Jets are commonly E-mail: [email protected] (LD); [email protected] (RB) associated with black hole spins (e.g. Blandford & Znajek 1977; Tchekhovskoy, Narayan & McKinney 2011). However, these models are more complicated. Here, we propose another way of constraining SMBH spins using tidal effects induced by star–SMBH close passages. There are two extreme cases of tidal interactions between an SMBH and a star. The most violent scenario is a tidal disruption event, which is usually associated with a star approaching an SMBH in a parabolic or very eccentric orbit. If the stellar radius at pericentre is less than the tidal radius, the star will be torn apart and about half of the stellar debris will eventually be swallowed by the hole (e.g. Rees 1988). One possible instance of a tidal disruption is the recent luminous Swift J1644+57 observation (e.g. Bloom et al. 2011; Zauderer et al. 2011). Optical and UV flares have also been associated with tidal disruptions (e.g. van Velzen et al. 2011; Gezari et al. 2012). These events are very luminous but hard to model. One the other side, if a star approaches an SMBH in equatorial, quasi-circular orbits which shrink due to gravitational radiation, material will be tidally stripped from the star when the star fills the Roche lobe. Such tidal stripping is a stable process lasting for ∼10 000–100 000 years. This steady mass transfer is not as luminous as a tidal disruption event, but lasts for much longer, and is easier to model with high precision for the mass and spin of the hole. C 2013 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society Roche accretion of stars close to massive black holes Mass transfer in binary systems has been well studied for a variety of stellar components and mass ratios. This investigation is mostly focused on binary systems composed of a stellar-mass compact object and a companion star (e.g. Webbink 1985). In particular, Hjellming & Webbink (1987) and Soberman, Phinney & van den Heuvel (1997) explored the stability of adiabatic mass transfers based on polytropic stellar models, and showed that such mass transfer processes could be dynamically unstable. However, these studies are based on Newtonian calculations, since the companion star reaches the Roche limit far away from the ISCO of the stellarmass compact object due to their comparable masses. However, in extreme mass-ratio inspiral (EMRI) systems, stars can fill the Roche lobe near ISCO of the massive black hole. Therefore, the signals from such binary systems can in principle be used to constrain the mass and spin of the SMBH, and general relativity needs to be involved when calculating the stellar dynamics. In this paper, we present an accurate general relativistic framework for studying such mass transfer. We shall be especially concerned when the mass transfer is ‘adiabatic’ – in other words it happens on a timescale short compared with the thermal (Kevin–Helmholtz) timescale of the star yet long compared with the dynamical time-scale. This problem has previously been considered by King & Done (1993) and Hameury et al. (1994) in connection with the (erroneous) identification of periodicity in NGC6814 (Madejski et al. 1993). King & Done (1993) mainly calculated the stellar orbital evolution of a Roche lobe filling convective star, using pseudo-Newtonian treatments. Hameury et al. (1994) further calculated the evolution using relativistic expressions (without considering the SMBH spin) and a full binary stellar evolution code. This paper builds upon King & Done (1993) and Hameury et al. (1994), and presents a more comprehensive treatment of stars and planets. We also show for various classes of stars, how general relativity (including the SMBH spin) and the stellar rotational velocity around its own axis affect the Roche volume calculation, the stellar orbital evolution and the accretion rate on to the SMBH. We will also conclude that such accretion events are stable over dynamical time-scales, in contrast to the suggestions of Hameury et al. (1994). The differences in results are mainly caused by more accurate calculation of the relativistic Roche potential (and the resultant Roche volume), as well as more accurate gravitational radiation power calculation. In an EMRI system, the stellar orbit shrinks due to the loss of energy and angular momentum through gravitational radiation and other torques. This gravitational radiation is a prime target of space-borne gravitational radiation observatories intended to test the general theory of relativity in the strong field regime. An EMRI system can also emit electromagnetic radiation (principally X-rays and optical), and its observation can also provide the basis of a good test of general relativity, if measurements are carried out accurately. We restrict our discussion of EMRI systems to stars in circular, equatorial orbits. If and when such a star overflows its Roche lobe close to the ISCO, gas will flow through the inner Lagrange point (L1), and the stellar orbit will usually expand – an extreme mass-ratio outspiral or EMRO! We shall assume initially that the stellar evolution is adiabatic, and discuss later the limitations of this assumption. Under these conditions, the entropy variation with mass within the star is frozen and the evolution is as described in our previous paper (Dai, Blandford & Eggleton 2011). In the Newtonian regime it is the sign of (∂ρ̄ /∂m )s , where m is the remaining mass of the star, and ρ̄ is the stellar mean density, that determines whether or not an EMRI becomes an EMRO. 2949 The gas flowing through L1 will generally have enough angular momentum to form a torus orbiting the hole. A hotspot will be formed where the stream hits the disc. In addition, it is quite possible that one or two spiral arm shock waves will be created in the disc. The hotspot and shocks will orbit the hole with the orbital frequency of the star, and an inclined observer should see a modulation of the radiation from the disc, most likely in the X-ray range. Clearly, the modulation will be the greatest when the observer is in the equatorial plane, observing maximal Doppler boosting. Here, we assume that the disc is formed by the stripped stellar material. However, an accretion disc could have existed before the star comes in. In such a scenario, the interaction between the disc and the star can alter the stellar evolution as discussed in Syer, Clarke & Rees (1991). We also have not considered how the external radiation pressure affects the star and the Roche evolution. Such irradiation effects are important for large mass-ratio binaries, as they can change the Roche geometry and may stop the Roche lobe overflow through L1 (in which case overflow will happen at L2 and an ‘excretion’ disc forms instead of an accretion disc) (see Drechsel et al. 1995; Phillips & Podsiadlowski 2002). One observation that may possibly be explained by this model is the quasi-periodic XMM–Newton X-ray signal from a nearby active type 1 Seyfert galaxy RE J1034+396 observed by Gierliński et al. (2008). As shown in fig. 1 in their paper, the rest-frame X-ray emission has a period of ∼ 1 h. The X-ray flux has a modulation of ∼0.1 associated with the nucleus of the galaxy, and this modulation could be well fitted by a cosine function. We will show in this paper that the Roche accretion of a lower main-sequence star or a brown dwarf on to a ∼106 –107 M SMBH can produce emissions with ∼1 h period, though the mass accretion rate under adiabatic condition and pure gravitational radiation torque might be too low to produce the observed X-ray flux. In Section 2, we compute the stellar evolution in the Newtonian limit. We show how the star spirals in, when it reaches the Roche limit, and how the stellar orbits and mass transfer rate evolve taking a lower main-sequence star as a representative. We also discuss the final fate of stars after they have undergone Roche evolution. In Section 3, we incorporate general relativity into the calculations, and compare results with those in Hameury et al. (1994). In Section 4, we apply our model to analyse the RE J1034+396 data, and obtain constraints on the stellar mass and SMBH parameters. Discussion and future research directions are laid out in the last section. 2 N E W T O N I A N T R E AT M E N T A Newtonian treatment brings out the salient features of the problem. Consider a fully convective lower main-sequence star with mass m ≡ mM as a representative, and assume that it orbits a hole with mass M ≡ 106 M6 M , in a circular, equatorial orbit with radius r. (The stellar orbit circularizes quickly either through gas dynamical dragging or gravitational radiation). We take G = c = 1 henceforth. From Kepler’s third law, the orbital radius [measured in units of the gravitational radius Rg (M6 ) = 1.47 × 1011 M6 cm] is r̃ = r 2 −2 = 23.9 P 3 M6 3 , Rg where the orbital period P is measured in hours. The orbital angular momentum √ 2/3 L = m Mr = 4.9 mM6 P 1/3 (1) (2) measured in units of M c Rg (M6 = 1) = 8.8 × 1054 g cm2 s−1 . 2950 is L. Dai and R. Blandford Likewise, the orbital energy, in units of M c2 = 1.79 × 1054 erg, E = −0.5 Mm r −1 = −0.021 × mM6 P −2/3 . 2/3 (3) 2.1 Pre-Roche evolution We assume that the star follows a circular orbit which shrinks slowly due to gravitational radiation. In the Newtonian limit, the torque is 32 4/3 2 2π 7/3 = 1.7 × 1044 m2 M6 4/3 P −7/3 g cm2 s−2 . M m τ= 5 P (4) The time it takes for the star to plunge into the hole via gravitational radiation is (assuming the star starts to inspiral from radius r rISCO , which is the radius of the ISCO, see Section 3.1) tinspiral = r4 5 −2/3 8 ∼ 987 m−1 M6 P 3 yr 256 M 2 m (5) (e.g. Misner, Thorne & Wheeler 1973). Typically, if the orbit is at least mildly relativistic, the inspiral time is less than a Kelvin–Helmholtz time. 2.2 Roche limit In order to obtain the Roche volume of the star, it is convenient to transform from an inertial frame to a reference frame rotating with the binary system. The Roche potential (r) for an arbitrary point with position vector r in the rotating frame is given by (r)N = − M 1 m − − |ω × r|2 |r − r | |r − r BH | 2 (6) (Frank, King & Raine 2002), where r and r BH are the position vectors of the centres of the star and the hole and ω is the rotational angular velocity of the star in the inertial frame. ω satisfies 0 ≤ ω = kω , where is the stellar orbiting angular velocity and 0 ≤ kω = ω/ 1 is the dimensionless rotational parameter. The first term in this expression of the Roche potential is the gravitational potential of the star itself, making the approximation that the stellar mass is centrally concentrated and only retaining the monopolar term. The second term is the gravitational potential of the hole. The third term is the potential in the rotating frame due to the centrifugal force. We erect a locally orthonormal coordinate system {x, y, z} centred at the star and directed along the {r, θ , φ} coordinate axes. If we expand the gravitational potential near the star, and note that the linear term is cancelled by a linear term in the expansion of the centrifugal potential, the Newtonian tidal tensor contributes quadratic terms − 21 (φ,xx x 2 + φ,yy y 2 + φ,zz z2 ), where φ , xx = −2φ , yy = −2φ , zz = 2M/r3 . (Note that the contribution to the Laplacian vanishes as it should.) The expansion of the centrifugal term also contributes a quadratic form − 12 ω2 (x 2 + z2 ) or − 12 kω 2 (x 2 + z2 )M/r 3 . A corotating star will expand to fill an equipotential surface. However, it is rather unlikely that an inspiralling star will be able to maintain corotation. If it rotates with angular frequency ω < , the quadratic centrifugal term is approximately − 12 ω2 (x 2 + z2 ), which includes the effects of Coriolis force. (There may be some internal circulation which will vitiate adiabaticity, but we do not consider this possibility further in this paper.) The gradient of the Roche potential vanishes at the first and second Lagrange point. Using this we can calculate that the first and second Lagrange points (L1 and L2) are located at 1/3 m xR = ±0.794 X (kω ) r (7) M from the star, where X (kω ) decreases from 1 to 0.87 as kω goes from 0 to 1. The transverse extents of the equipotential Roche surface are 1/3 m r (8) yR = ±0.497 Y(kω ) M and m zR = ±0.497 Z(kω ) M 1/3 r, (9) where Y(kω ) goes from 1 to 0.89 and Z(kω ) goes from 1 to 0.93 as kω goes from 0 to 1. Furthermore, the volume contained by the Roche surface is numerically calculated to be −1 kω m 3 (10) VN 0.683 × +1 r . 2.78 M Combining this with equation (1), we obtain the relationship between the orbital period and the mean density of the star when the star reaches the Roche limit: 1/2 kω +1 ρ̄ −1/2 (11) P0 8.18 2.78 independent of the mass of the hole, where P0 is the stellar orbital period in hours and ρ̄ the mean density of the star in cgs unit. For a dwarf star with mass ∼0.3 M , the Roche condition is satisfied when the orbital period decreases to a value of 1.96–2.27 h depending on kω . If we suppose that the star plunges when its orbital radius reaches 6Rg , the maximum black hole mass for there to be mass transfer from a ∼0.3 M dwarf star is 1.6–1.9 × 107 M (kω = 0 − 1). 2.3 Roche evolution To set the stage for the relativistic treatment below, let us consider the classical Roche model for a dwarf star that is non-rotating. As discussed in Dai et al. (2011), its mean stellar density decreases monotonically as mass is tidally stripped, leading to expanding orbits. Assuming that the star has initial mass m0 ∼ 0.3 M , we found that m 2 ρ̄ (12) ρ̄0 m0 (Dai et al. 2011), where ρ̄0 is the initial mean density of the star. Combining equations (11) and (12), the orbital period P of the dwarf star, in terms of the initial period P0 , varies with its remaining mass as m −1 P . (13) P0 m0 We assume that the orbit of the star will change in a way that the angular momentum of the system is conserved so that dL = 0, (14) dt where we continue to assume that gravitational radiation provides the only torque as in equation (4). τ+ Roche accretion of stars close to massive black holes Figure 1. Shows the period–mass relationship of a non-rotating lower mainsequence star with initial mass ∼0.3 M accreting on to a Schwarzschild SMBH with mass 107 M . The vertical dashed line represents the evolution prior to the mass transfer, when the mass does not change, but the orbit shrinks. The black curve describes the subsequent Roche evolution in the Newtonian limit, when the stellar orbit expands. The blue curve shows how stellar orbital period changes with mass in the relativistic limit. Using equations (2), (4), (13) and (14), we see that the time elapsed since the onset of mass transfer goes as 11 P 3 8 − 23 3 − 1 yr. (15) tR ∼ 1.46 × 10 M6 P0 3 P0 The associated rate of mass transfer of Roche accretion, also obtained assuming conservation of angular momentum, is ṁ ∼ −1.18 × 1022 M6 m2 P −8/3 2/3 ∼ −2.55 × 1023 P0 −8/3 M6 m14/3 g s−1 . 2/3 (16) We show in Figs 1 and 2 the Newtonian evolution of the orbital period and the mass transfer rate for this lower main-sequence star when it is adiabatically stripped by a Schwarzschild SMBH with mass 107 M . We see that the stellar orbital period can increase to around five times the Roche period as the star is stripped. However, the majority of the Roche evolution happens on a time-scale 2951 ∼105 yr, and the stellar mass transfer rate decreases rapidly with time. As discussed in Dai et al. (2011), the angular momentum of the binary system as in equation (2) follows L ∼ m P 1/3 ∼ m ρ̄−1/6 . If the mass accretion is stable on dynamic time-scales, it means that L should decrease as mass is accreted on to the hole. If we define a parameter η = d(ln ρ̄ )/d(ln m ), then η < 6 is the criterion for Newtonian dynamical stability. A more intuitive method of calculating the dynamical stability, as Hameury et al. (1994) suggested, is to ensure that the star will not exceed the Roche surface under the assumption of no angular momentum loss through gravitational radiation. In other words, dv /dm < dv N /dm , where v = m /ρ̄ is the stellar volume and v N is the Roche volume in equation (10). This leads to the same stability criterion η < 6 in the Newtonian limit. Continuing to consider the ∼0.3 M lower main-sequence star, we have η = 2 (equation 12), ensuring that the star recedes from the hole stably under Roche accretion. For many purposes, this Newtonian treatment applied to a broader range of stars will be adequate. However, as we want to exhibit relativistic effects, we will turn to a relativistic treatment, again building upon the discussion in Hameury et al. (1994). 2.4 The final fate of stars under Roche evolution After mass transfer begins, most stars will recede from the black hole so that they will continue to fill their Roche lobes as their densities drop. The star loses mass less and less rapidly. Eventually, the star will cool and its structure will evolve when the thermal time-scale of the star is comparable with the mass-loss time-scale. The motion of the star under the action of gravitational radiation will then cause the orbit to shrink again. Additional excursions are possible depending upon the structure of the star. The star might cool down adiabatically or on a thermal time-scale. Therefore, it might or might not fill its Roche lobe on its return. The equilibrium properties of these stars will be even more complicated than the ones we calculated in this paper. Ultimately, the star crosses the ISCO and plunges into the hole. When more than one star is in a ‘parking orbit’, there is a possibility of resonant dynamical interaction analogous to that which occurs in planetary systems. In addition interaction with the disc can alter the evolution. We illustrate the scenario in Fig. 3. Around the hole, there will be stars migrating inward and mostly outward under Roche evolution, stars that have gone at least one Roche evolution ‘parked’ in a belt or cloud (if the stars approached the hole in inclined orbits) waiting to move in again, new stars coming in under gravitational radiation and still burning and stars coming in again from the ‘parking cloud’ and probably not burning. 3 R E L AT I V I S T I C T R E AT M E N T 3.1 Gravitational radiation emission Figure 2. Shows the maximum accretion rate when a non-rotating lower main-sequence star with initial mass ∼0.3 M is tidally stripped by a Schwarzschild SMBH with mass 107 M . The black curve is the Newtonian mass accretion rate. The blue dashed curve is the relativistic mass accretion rate calculated without considering the angular momentum loss at the ISCO. The blue solid curve is the relativistic mass accretion rate calculated with the angular momentum loss at the ISCO assuming a steady accretion disc has formed. A star (with a mass m ) orbiting a massive hole (with mass M and spin parameter a) will lose angular momentum and energy due to gravitational radiation. The orbit is circularized quickly, as we have shown in Dai et al. (2011). However, the rate of alignment of the spin and orbital angular momentum is slower than the inspiral rate (e.g. Hughes 2001). For the moment, we restrict attention to circular equatorial orbits. ˜ = M, We define a dimensionless orbital angular velocity where is the regular orbital angular velocity. The relativistic 2952 L. Dai and R. Blandford Figure 3. For an ∼0.3 M star orbiting a Schwarzschild SMBH with mass 107 M as discussed, the star reaches the Roche limit when the orbital radius is ∼12Rg , and then spiral out. When the evolution time-scale is comparable with Kevin–Helmholtz time-scale, the star stops receding at the ‘parking belt’, which is located at around ∼25Rg from the hole. Here, we calculated the Kevin–Helmholtz time-scale of the star using the initial stellar internal energy and luminosity. Kepler’s law is ˜ = 1 , r̃ 3/2 ± a (17) where ± refers to prograde and retrograde orbits. The star plunges into the hole when the orbit shrinks to the size of the ISCO. The ISCO of a black hole has a radius rISCO that satisfies the equation as given in Bardeen, Press & Teukolsky (1972): rISCO 2 rISCO rISCO + 8a − 3a 2 = 0. −6 (18) M M M This radius equals Rg , 6Rg and 9Rg for a = 1, 0 and −1, respectively. The realistic upper limit for a has been argued to be ∼0.998 (Thorne 1974) for which rISCO ∼ 1.237Rg . Let us now consider how to modify the treatment of gravitational radiation. Finn & Thorne (2000) computed a relativistic inspiral time-scale for a star with orbital radius r: tinspiral = M2 5 T ˜ 8/3 256 m (19) in the Boyer–Lindquist coordinate, where T is the relativistic correction term. They presented their results of T in tabular form. We have used these results to give a convenient interpolation formula for T : T (x, a) = (x − 1)2 + b (x − 1) , x2 + d x + f where x = r/rISCO . Clearly, T (x = 1, a) = 0 and T (x, a) → 1 as x → +∞ (since at infinity the relativistic effects are negligible). This formula holds true with a maximum error of ∼2 per cent for x ≥ 2. A contour plot of T (x, a) is given in Fig. 4. The orbital evolutions of an ∼0.3 M lower main-sequence star spiraling in around a 107 M SMBH of different spins are shown in Fig. 5. It is clear that the standard ‘Landau–Lifshitz’ form equation is adequate for r > 100Rg . The gravitational radiation reaction torque is readily computed from ∂L ∂tplunge , (21) τ (r, a)gr = − ∂r a ∂r a where the relativistic orbital angular momentum is L = m M (1 − 2r̃ −1 + a ∗ r̃ −3/2 ) r 1/2 (1 − 3r̃ −1 + 2a ∗ r̃ −3/2 )1/2 (22) (e.g. Bardeen et al. 1972). 3.2 Roche limit We use the same approach as the Newtonian treatment. The local tidal tensor is evaluated in Boyer–Lindquist coordinates and then transformed into an orthonormal coordinate basis freely falling with the star. This ‘local’ coordinate system has spatial coordinates directed along the {r, θ , φ} axes. The only non-zero spatial elements of the tidal potential can be expressed in the form φ,xx = − b = 1.06 × (1 − a)−0.11 − 0.66, 1+k , r3 φ,yy = k , r3 φ,zz = 1 , r3 (23) where d = 0.57 × (1 − a)0.23 − 2.22, k= f = [−1.08 × (1 − a)1.08 − 3.18 × (1 − a)−0.018 ] × e−0.23×(1−a) + 4.91, Figure 4. The relativistic correction term T to the plunge time as a function of the distance to the hole and the spin parameter of the hole. We see that as the radius goes to infinity, T becomes 1 no matter what the spin is. Also T = 0 when r = RISCO . (20) r̃ 2 − 4a r̃ 1/2 + 3a 2 . r̃ 2 − 3r̃ + 2a r̃ 1/2 (24) This satisfies Laplace’s equation as it must. k is a parameter expressing the importance of relativity. In the Newtonian limit, Roche accretion of stars close to massive black holes Figure 5. The total inspiral time in years as a function of r/Rg for an ∼0.3 M lower main-sequence star orbiting a 107 M SMBH. The three solid curves are accurate inspiral times obtained using equation (19) with different colours indicating different spin parameters of the hole: −0.99 (green), 0 (black), and 0.999 (blue). The red dotted curve is the Newtonian inspiral time using equation (5). All relativistic evolution curves approach the Newtonian curve when the star is far away from the hole. k → 1 as r → ∞ as must also be true. When the motion of the star is extremely relativistic, or r → rISCO , note the remarkable property that k → 2 independent of a. Combining the stellar gravity and rotation with the tidal contribution, the Roche potential becomes (x, y, z) = −m (x 2 + y 2 + z2 )−1/2 + M ky 2 + z2 − (k + 1)x 2 r3 2 M kω (x 2 + y 2 ) . (25) r3 2 Here, kω is the dimensionless stellar rotational parameter in the inertial frame. Generalizing the Newtonian treatment, we can locate the first and second Lagrange points and the transverse extents of the Roche equipotential surface, and calculate the volume contained by this surface. When k = 2, the Roche volume at the ISCO is calculated to be −1 kω m 3 +1 r . (26) VISCO = 0.456 × 4.09 M − 2953 Figure 6. Contour plot of Roche volume (in units of m /Mr3 ) as a function of k and kω for a = 0.5. k = 1 is the Newtonian limit. k = 2 when the star is at the ISCO. kω = 0 when the star is non-rotating, and kω = 1 when the star is corotating. This expression of this relativistic Roche volume seems complicated but has similar terms as its Newtonian companion equation (10). It is basically m /Mr3 times a function which depends on k and kω . As an example, Fig. 6 is the contour plot of V(r, a, kω ) in units of m /Mr3 around a black hole with a = 0.5 as a function of k and kω . This dimensionless volume decreases from 0.683 to 0.456 from the Newtonian limit to the most extreme relativistic limit (r = rISCO ) for a non-rotating star. For a corotating star, the variation is from 0.503 (Newtonian) to 0.368 (at the ISCO). These limits of the Roche volume at four corners of the plot are independent of the spin. Comparing with Fishbone (1973)’s Roche volume used by Hameury et al. (1994), our model is more accurate in the relativistic limit, including black hole spin and the stellar rotation. We plot the ratio of the two Roche volumes in Fig. 7. Using equations (17) and (27), we obtain a relationship linking the orbital radius or period to the stellar mean density when the star And more generally, an adequate approximation to the Roche volume is VN (r, kω ) − VISCO (r, kω ) V (r, a, kω ) = VN (r, kω ) − 0.5 , r r + F (a, k ) − 1 ω rISCO rISCO (27) where 13.9 14.8 ×(1 − a)0.02 + 23.8 − F (a, kω ) = −23.3 + kω + 2.8 2.8 + kω 0.4 + 0.9 − × (1 − a)−0.16 . 2.6 + kω Figure 7. The ratio of our Roche volume to Fishbone (1973)’s Roche volume, as a function of r/Rg . The four solid curves are for non-rotating stars orbiting a black hole with a = −1 (red), 0 (blue), 0.9 (black) and 0.998 (green). The lower four dashed ones are for corotating stars correspondingly. 2954 L. Dai and R. Blandford fills its Roche lobe: −1 kω −1 kω 0.683 × 2.78 +1 − 0.456 × 4.09 +1 M = − 0.5 ρ̄ r 3 r r + F (a, kω ) rISCO −1 rISCO + 0.683 × −1 kω . +1 2.78 (28) This equation describes how the evolution of the stellar properties (mainly the mean density) governs the orbital evolution during the Roche accretion phase. 3.3 Roche evolution As a star approaches an SMBH in a circular, equatorial orbit and spirals in due to gravitational radiation, if the Roche radius is larger than the ISCO size, the star can start to transfer mass to the hole through Roche accretion. Otherwise, the star directly plunges into the hole. After the Roche accretion starts, whether the star continues to spiral in or spiral out is mainly governed by the change of the stellar mean density. The detailed picture of the orbital evolution depends on the mass and spin of the black hole, the stellar type and parameters, and the rotational parameter kω of the star. Assuming that the stellar mass-loss time-scale is faster than the thermal timescale, we have calculated in Dai et al. (2011) how the stellar mean density responds to the adiabatic mass-loss for a representative set of stars and planets. Then using equation (28), we can check where the star reaches the Roche limit, and how the stellar orbit evolves through the Roche accretion. As the stellar orbital evolution scenarios are different for different classes of stars, and when the SMBH has different masses or spins, we illustrate the results in Fig. 8 (SMBH has mass M = 108 M ), Fig. 9 (M = 107 M ), Fig. 10 (M = 106 M ) and Fig. 11 (M = 105 M ). We also show one example of the Roche accretion of planets in Fig. 12. The representative stars and planets used in this paper are the same as in Dai et al. (2011), except that here we use a white dwarf of initial mass 0.8 M instead of 1 M for better illustration. The initial conditions of the stars and planets are: a Sun-like star, a lower main-sequence star with initial mass ∼0.3 M and radius ∼0.3 R , an upper main-sequence star with initial mass ∼7.9 M and radius ∼3.3 R , a red giant with initial mass ∼7.1 M and radius ∼183 R , a white dwarf with initial mass 0.8 M and radius ∼0.01 R , a brown dwarf with initial mass 0.05 M and radius ∼0.09 R , a Jupiter-like planet and an Earth-like planet. As we discussed in Dai et al. (2011), stars with convective structures such as a lower main-sequence star, a brown dwarf or a white dwarf will become less dense as they adiabatically lose mass, so they will spiral out throughout the Roche accretion phase. A Sun-like star, an upper main-sequence star or a red giant will spiral in first and eventually spiral out. Gas planets will spiral out, and solid planets will stay in almost the same orbit. Since we are mostly interested in how general relativity affects the Roche accretion, we only plot the region between the ISCO and 50Rg from the hole. There are four different curves for each star in each plot: using general relativistic or Newtonian calculations, and using a non-rotating or corotating star. The overall behaviours of the four curves are similar, but quantitatively distinct. An SMBH with mass ∼108 M has a density too low for Roche accretion to happen for most types of stars (Hills 1975). Stars just plunge into the hole, unless the hole is rapidly spinning so the ISCO size is small. The extended envelope of a red giant will be stripped, but its compact core will still plunge in. If the hole is spinning very fast (prograde a ≥ 0.9), main-sequence stars can start Roche accretion very close to the ISCO, where general relativity is crucial. It is possible for a Sun-like star or an upper main-sequence star to plunge during the mass transfer since they continue to spiral in for some time after reaching the Roche limit close to the ISCO. For an SMBH with mass ∼107 M , main-sequence stars usually go through Roche accretion till the mass-loss rate is too low, afterwards the star will park within 50Rg of the hole. Most of the Roche accretion processes happen in the regime where general relativity is important. Very compact stars like white dwarfs and the core of red giants still directly plunge in. A brown dwarf can be tidally stripped if the black hole is spinning fast. For an SMBH with mass ∼106 M , main-sequence stars and brown dwarfs will reach the Roche limit outside 15Rg of the hole, where general relativistic effects are not prominent. White dwarfs usually plunge in, unless the white dwarf has a small initial mass (thus being puffier) and the hole is spinning with an almost maximal spin. Red giants can reach the Roche limit in certain cases. For a massive black hole with mass ∼105 M , white dwarfs can go through relativistic Roche accretion process. All other types of stars reach the Roche limit far away from the hole. We also plot the period–mass relationship for planets as in Fig. 12. A solid planet like Earth has the same period as mass decreases due to unchanging density. A gas giant like Jupiter evolves so that the period will increase. From the conservation of angular momentum, the stellar angular momentum L satisfies L = LGR + LISCO + Lwind + Lrad . (29) Here, L is the angular momentum of the star as in equation (22). LGR is caused by the gravitational radiation torque τ gr calculated from our new model (equation 21). LISCO represents the loss of angular momentum of materials falling into the ISCO after a steady accretion disc forms. Lwind is the angular momentum of materials lost from an accretion disc. In this inner region of an SMBH accretion disc, there are mainly two kinds of wind: radiation pressure driven and magnetic driven. The radiation pressure induced wind carries its own specific angular momentum as on the disc. The magnetic induced wind carries away very little mass but very large angular momentum which can be 10 times as large as its specific angular momentum on the disc. The last term is the disc radiation torque due to beaming effect, which is relatively small compared with the others. In order to understand the details of Lwind , three-dimensional dynamical simulations need to be employed. In this paper, for simplicity, we consider only the first two L’s, the mass accretion rate on to the hole will be the same as the stellar mass deposition rate. We then have dL + τgr + ṁ jms = 0. dt (30) jms is the specific angular momentum at the ISCO, which equals 0 in the Newtonian limit. Using this equation (30) and equations (22) and (21), we can calculate the relativistic Roche mass accretion rate. For the same lower main-sequence star with initial mass ∼0.3 M , we show its P and ṁ change in the relativistic limit in comparison with the Newtonian results in Figs 1 and 2. With the consideration of general relativity, the star does not recede from the hole as far as in the Newtonian limit, and the mass accretion rate decreases by at least half comparing with the Newtonian rate. After a steady accretion disc forms, the accretion rate drops due to the torque exerted on the star due to the angular momentum loss at Roche accretion of stars close to massive black holes 2955 Figure 8. The orbital evolutions of different stars during the Roche evolution close to an SMBH with M = 108 M , and a = −1, 0, 0.9 and 0.998. In each panel of this figure and the Figs 9–11, the x-axis is the stellar mass ranging from 0.005 to 15 M , and the y-axis is the stellar orbital period from PISCO to the orbital period at 50Rg . We use different colours to denote different types of stars in this figure and the following three figures: a Sun-like star (red, Sun), a lower main-sequence star with initial mass ∼0.3 M (purple, LM), an upper main-sequence star with initial mass ∼7.9 M (blue, HM), a red giant with initial mass ∼7.1 M (orange, RG), a white dwarf with initial mass 0.8 M (cyan, WD) and a brown dwarf with initial mass 0.05 M (brown, BD). We plot four curves for each star for comparison: the relativistic calculation for a corotating star (coloured solid curve), the relativistic calculation for a non-rotating star (coloured dashed curve), the Newtonian calculation for a corotating star (black solid curve) and the Newtonian calculation for a non-rotating star (black dashed curve). The vertical blue dashed lines represent the phase of inspiral due to gravitational radiation. Note that in panel (c) and (d), whether the Sun or a lower main-sequence star can go through Roche evolution as indicated by the model may depend on its stellar rotational parameter and whether including relativity in the calculations. the ISCO. For other convective stars like a white dwarf or a brown dwarf, since their density evolution basically follows equation (12), their Roche accretion rates can be calculated using equation (16) with relativistic corrections. For an upper main-sequence star, its density does not monotonically decrease during the Roche mass transfer, and an analytical formula for the mass accretion rate is hard to obtain. Therefore, we illustrate in Fig. 13 the Newtonian mass transfer rate 2956 L. Dai and R. Blandford Figure 9. The orbital evolutions of different stars during the Roche evolution close to an SMBH with M = 107 M , and a = −1, 0, 0.9 and 0.998. The axes, the plotting ranges and the stars are the same as in Fig. 8. for an ∼7.9 M upper main-sequence star when it is stripped by an SMBH with mass M = 106 , 107 or 108 M . The relativistic accretion rate will be in the same order as the Newtonian rate, if the star reaches the Roche limit outside the ISCO. After the mass transfer starts, the accretion rate increases first, as the upper main-sequence star continues to spiral in. Eventually, the accretion rate drops when the star starts to spiral out. If a star is stripped by an SMBH with larger mass, the accretion rate will be larger. We can apply the same method as in Section 2.3 to check if the stellar volume will exceed the Roche volume in the evolution. In contrast to Hameury et al. (1994), we find that the Roche evolution of a fully convective star is always stable on dynamical time-scales in the relativistic limit, even when the orbital separation is less than ∼3Rg (in which case the SMBH needs to have a value 0.8). 4 A P P L I C AT I O N T O R E J 1 0 3 4 +3 9 6 We illustrate the use of our calculations by applying them to the observations of RE J1034+396. The ∼1 h period of the X-ray signal can be associated with the Roche accretion of a lower mainsequence star or a brown dwarf. The lower main-sequence star needs to have a mass m 0.13 M if the star is non-rotating or m 0.11 M if the star is corotating. The SMBH has mass ∼106 M – 107 M (the latter with a rapid prograde spin). Using equation (16), Roche accretion of stars close to massive black holes 2957 Figure 10. The orbital evolutions of different stars during the Roche evolution close to an SMBH with M = 106 M , and a = −1, 0, 0.9 and 0.998. The axes, the plotting ranges and the stars are the same as in Fig. 8. the maximum Roche accretion power of a 0.13 M star is 1.9–8.8 × 1041 erg s−1 (for M = 106 –107 M ). Assuming a 10 per cent SMBH radiative efficiency, the X-ray radiation power LX associated with the Roche accretion is less than 8.8 × 1040 erg s−1 . The X-ray luminosity of RE J1034+396 is 4 × 1043 erg s−1 with ∼9.2 per cent variation. Therefore, the Roche accretion power from a lower mainsequence star is still one to two orders below the observed modulation of RE J1034+396. However, with a very large observer inclination angle (almost along the equatorial plane), Doppler shift will enhance the modulation and can boost the observed LX by a factor of ∼10 if the star is rotating close to the ISCO. 16 cycles of signals are insufficient to be confident of a stable period. However, if more signals can be observed, information on spin could be extracted. In principle, the closer the star is to the hole, or the larger the spin of the hole, the stronger the emission would be. This could introduce an important selection effect as rapidly spinning holes would be more likely to be observed. Fig. 14 shows simulated periodic signals from a hotspot orbiting around the hole circularly using our numerical package. The modulation of intensity comes from beaming effect – as we have discussed in Dai, Fuerst & Blandford (2010), the bolometric intensity varies as the fourth power of the Doppler-gravitational redshift factor 2958 L. Dai and R. Blandford Figure 11. The orbital evolutions of different stars during the Roche evolution close to a massive black hole with M = 105 M , and a = −1, 0, 0.9 and 0.998. The axes, the plotting ranges and the stars are the same as in Fig. 8. g = E/E0 (Cunningham 1975). Here, E is the energy of the photon in the observer’s frame and E0 in the hotspot’s inertial rest frame. The shape of the curves and the relative intensity are different for different spins and observer inclination angles. Also a major difference between producing a hotspot on the disc using the Roche accretion model and other QPO models (e.g. Wagoner, Silbergleit & Ortega-Rodrı́guez 2001) is that the period of the modulation predicted by the Roche model will follow that of the stellar orbit although there is likely to be phase noise. Since the stellar orbit evolves on an adiabatic time-scale, the period of the hotspot will also change very slowly but in a consistent way. For example, the orbital period of a ∼0.3 M lower main-sequence star can change by as much as ∼0.1 s in one year during Roche evolution, in which time we could observe ∼8000 cycles of modulation. A P − Ṗ analysis similar to the double pulsar’s could in principle be done to extract information of the stellar type and black hole mass and spin. 5 DISCUSSION AND FUTURE RESEARCH DIRECTIONS Observation of quasi-periodic X-ray emissions from stars orbiting black holes in active galactic nuclei (AGN) is a potential probe of general relativity. There is an intriguing possibility that the Roche accretion of stars close to massive black holes 2959 Figure 14. Simulated X-ray period emission from a hotspot orbiting a black hole is shown for various black hole spin parameters a and observer inclination angles i. In each figure, the hotspot orbits the hole in equatorial circular orbits with the same radius. The units of the axes are arbitrary but consistent. Figure 12. Roche evolutions of a perfect solid planet with the mass of Earth (green) and a perfect gas planet with the mass of Jupiter (pink) near a maximally spinning SMBH with M = 108 M and a = 0.998. the x-axis is the stellar mass ranging from 10−8 to 10−3 M , and the y-axis is the stellar orbital period from PISCO to the orbital period at 50Rg . We plotted four curves for each planets for comparison: the relativistic calculation for a corotating star (coloured solid curve), the relativistic calculation for a non-rotating star (coloured dashed curve), the Newtonian calculation for a corotating star (black solid curve) and the Newtonian calculation for a nonrotating star (black dashed curve). The vertical blue dashed lines represent the phase of inspiral due to gravitational radiation. Figure 13. Shows the maximum Newtonian accretion rate when a nonrotating upper main-sequence star with initial mass ∼7.9 M is tidally stripped by an SMBH of mass M = 106 (red), M = 107 (blue) or M = 108 (green). X-ray modulation of RE J1034+396 does come from an orbiting star filling its Roche lobe and being swallowed by an SMBH, in a similar fashion to the ‘black widow’ pulsar systems (e.g. Thorsett & Chakrabarty 1999) like PSR1957+20 (e.g. Kulkarni & Hester 1988; Kulkarni et al. 1992; Arons & Tavani 1993). It is reasonable to search the existing AGN X-ray data bases for periodic signals with relevant periods within the range of ∼1–20 h. Although we are mostly interested in this paper in observable tidally disrupted or stripped stars, it is quite reasonable to suppose that most dormant AGN black holes are orbited by a stable system of stars with a common orbital plane in much the same way as the sun is orbited by the planets. Other orbiting stars can also have dynamical interactions with the accreting star and make its evolution more complicated. The probability of observing such phenomena relates to the stellar distribution in accretion discs and the mechanisms of bringing stars close to the ISCO. When materials accrete on to the black holes, stars can be borne in the accretion disc, which could go away eventually leaving only the stars. Such stars would orbit the hole in bound orbits in the equatorial plane with low eccentricity. On the other hand, if an accretion disc still exists when the star gets close to the ISCO, the star will have an interaction with the disc as well, or might be travelling in the disc already producing a gap and perturbation patterns on the disc. The results so far suggest more follow-up investigations. If the stellar orbit has a small eccentricity, tidal heating would change the evolution of the star, and could even detonate or destroy the star. (The star is very fragile, as the ratio of stellar self-binding energy 2 2 /vorb ∼ 10−4 .) If the star can survive to the orbital energy ∼vesc orbiting the black hole in this slightly eccentric or inclined orbit, perturbations will be added to the sinusoidal behaviour of the Roche accretion flux. A main-sequence star expanding to become a red giant can increase the mass accretion rate. A jet can also be formed due to Roche accretion process or tidal disruption (as observed in the Swift J1644+57 event), and the jet power can be used to probe the black hole spin. Between a total tidal disruption and a steady Roche mass transfer, intermediate phases exist – maybe the star deposits part of its envelope at pericentre of each orbit. Threedimensional general relativistic hydrodynamical simulation can be done to understand the evolution of the star, the formation of the disc, and the interaction of the tidal stream and the disc. The most exciting and immediate prospect, however, is to open up a new line of investigation using existing X-ray data bases, to seek evidence for stars in bound orbits around massive black holes. Discovering such systems would elucidate the evolution of AGN and also open the possibility of a new measurement of black hole spin and ultimately of concluding that the Kerr metric correctly describes spinning black holes. 2960 L. Dai and R. Blandford AC K N OW L E D G E M E N T S This work was supported by the US Department of Energy contract to SLAC no. DE-AC02-76SF00515. We would like to give special thanks to our collaborator P. Eggleton for providing the stellar models, R. Wagoner for valuable comments, and D. Chakrabarty, J. Faulkner, L. Gou, J. Grindlay, S. Hughes, A. Loeb, H. Marshall, J. McClintock, M. Nowak, S. Phinney, M. Rees, R. Remillard, P. Schechter, R. Simcoe, A. Tchekhovskoy and N. Weinberg for useful discussions. REFERENCES Arons J., Tavani M., 1993, ApJ, 403, 249 Bardeen J. M., Press W. H., Teukolsky S. A., 1972, ApJ, 178, 347 Blandford R. D., Znajek R. L., 1977, MNRAS, 179, 433 Bloom J. S. et al., 2011, Sci, 333, 203 Cunningham C. T., 1975, ApJ, 202, 788 Dai L. J., Fuerst S. V., Blandford R., 2010, MNRAS, 402, 1614 Dai L., Blandford R. D., Eggleton P. P., 2011, preprint (arXiv:e-prints) Drechsel H., Haas S., Lorenz R., Gayler S., 1995, A&A, 294, 723 Finn L. S., Thorne K. S., 2000, Phys. Rev. D, 62, 124021 Fishbone L. G., 1973, ApJ, 185, 43 Frank J., King A., Raine D. J., 2002, Accretion Power in Astrophysics, 3rd edn. Cambridge Univ. Press, Cambridge, p. 398 Gezari S. et al., 2012, Nat, 485, 217 Gierliński M., Middleton M., Ward M., Done C., 2008, Nat, 455, 369 Gou L. et al., 2011, preprint (arXiv:e-prints) Hameury J., King A. R., Lasota J., Auvergne M., 1994, A&A, 292, 404 Hills J. G., 1975, Nat, 254, 295 Hjellming M. S., Webbink R. F., 1987, ApJ, 318, 794 Hughes S. A., 2001, Phys. Rev. D, 64, 064004 King A. R., Done C., 1993, MNRAS, 264, 388 Kulkarni S. R., Hester J. J., 1988, Nat, 335, 801 Kulkarni S. R., Phinney E. S., Evans C. R., Hasinger G., 1992, Nat, 359, 300 Li L.-X., Narayan R., McClintock J. E., 2009, ApJ, 691, 847 Madejski G. M., Done C., Turner T. J., Mushotzky R. F., Serlemitsos P., Fiore F., Sikora M., Begelman M. C., 1993, Nat, 365, 626 Misner C. W., Thorne K. S., Wheeler J. A., 1973, Gravitation. Freeman & Co., San Francisco Phillips S. N., Podsiadlowski P., 2002, MNRAS, 337, 431 Rees M. J., 1988, Nat, 333, 523 Reynolds C. S., Nowak M. A., 2003, Phys. Rep., 377, 389 Shcherbakov R. V., Penna R. F., McKinney J. C., 2010, preprint (arXiv:eprints) Soberman G. E., Phinney E. S., van den Heuvel E. P. J., 1997, A&A, 327, 620 Syer D., Clarke C. J., Rees M. J., 1991, MNRAS, 250, 505 Tchekhovskoy A., Narayan R., McKinney J. C., 2011, preprint (arXiv:eprints) Thorne K. S., 1974, ApJ, 191, 507 Thorsett S. E., Chakrabarty D., 1999, ApJ, 512, 288 van Velzen S. et al., 2011, ApJ, 741, 73 Wagoner R. V., Silbergleit A. S., Ortega-Rodrı́guez M., 2001, ApJ, 559, L25 Webbink R. F., 1985, Stellar Evolution and Binaries. Cambridge Univ. Press, Cambridge, p. 39 Yuan F., Quataert E., Narayan R., 2004, ApJ, 606, 894 Zauderer B. A. et al., 2011, Nat, 476, 425 This paper has been typeset from a TEX/LATEX file prepared by the author.