Global modeling of radiatively driven accretion of metals from compact debris disks onto the white dwarfs.
Recent infrared observations have revealed presence of compact (radii ) debris disks around more than a dozen of metal-rich white-dwarfs (WD), likely produced by tidal disruption of asteroids. Accretion of high-Z material from these disks may account for the metal contamination of these WDs. It was previously shown using local calculations that the Poynting-Robertson (PR) drag acting on the dense, optically thick disk naturally drives metal accretion onto the WD at the typical rate g s. Here we extend this local analysis by exploring global evolution of the debris disk under the action of the PR drag for a variety of assumptions about the disk properties. We find that massive disks (mass g), which are optically thick to incident stellar radiation inevitably give rise to metal accretion at rates . The magnitude of and its time evolution are determined predominantly by the initial pattern of the radial distribution of the debris (i.e. ring-like vs. disk-like) but not by the total mass of the disk. The latter determines only the disk lifetime, which can be several Myr or longer. Evolution of an optically thick disk generically results in the development of a sharp outer edge of the disk. We also find that the low mass ( g), optically thin disks exhibit and evolve on characteristic timescale yr, independent of their total mass.
Subject headings:White dwarfs — Accretion, accretion disks — Protoplanetary disks
Recent ground based and Spitzer observations (Zuckerman & Becklin 1987; Graham et al. 1990; Farihi et al. 2010) have revealed the presence of infrared excesses in spectra of a number of metal-rich white dwarfs (WD). These excesses have been interpreted (Jura 2003) as arising from reprocessing of stellar emission by a compact disk of refractory debris orbiting the WD. Detailed spectral modeling (Jura 2003; Farihi et al. 2010) suggests that these disks are geometrically thin and essentially flat, although in some cases there is certain evidence for flaring or warping (Jura et al. 2007b, 2009). Spectral fitting also finds the disks to be predominantly optically thick (Jura 2003; Farihi et al. 2010), even though this statement may not hold in some parts of the disk (Jura et al. 2007b, 2009). These properties make circum-WD debris disks look very similar to the ring of Saturn (Cuzzi 2010).
Outer radii of these disks are always found to be R (Jura et al. 2007; Farihi et al. 2010). This is close to the Roche radius , within which any object of normal density ( g cm) would be disrupted by the tidal forces of the WD. This observation motivated the idea (Jura 2003) that compact debris disks owe their origin to the tidal destruction of massive asteroid-like bodies, which were scattered into the low-periastron orbits by planets that have survived the AGB stage of stellar evolution.
Inner radii of disks are usually small, , where is the WD radius, although the exact number depends on the details of spectral modeling (Jura et al. 2007). In some cases is close to the distance at which a refractory body heated by the direct starlight would have a temperature equal to the sublimation temperature , suggesting that the inner cutoff of the disk is caused by particle sublimation. However, in some cases such interpretation implies very high value of (Jura et al. 2007; Brinkworth et al. 2009; Dufour et al. 2010). Nevertheless, it is always safe to assume that the ratio of the outer to inner disk radii is rather modest, . This will turn out being important later on, see §4.1.1.
A significant reservoir of refractory material in the immediate vicinity of the WD can naturally be responsible for contaminating the atmospheres of these WDs by metals (Jura 2003). This idea is strongly supported by the fact that the infrared excesses due to compact debris disks have been observed only around metal-rich WDs with rather high metal accretion rates (Farihi et al. 2009, 2010). Alternative model of the WD metal pollution by interstellar accretion (Dupuis et al. 1993a) runs into problem trying to account for the abundance ratios of different elements in the WD atmospheres, in particular, greatly suppressed hydrogen abundance (Dupuis et al. 1993b).
The new model of metal accretion from a circumstellar reservoir of high-Z material raises a natural question: how do metals get transported to the WD surface from the debris disk, which has an inner gap with the radius ? It is reasonable to assume that at disk particles get converted into metal gas by the intense heating due to the WD radiation, which then viscously accretes onto the WD. However, the issue of whether the debris disk can sustain the high metal accretion rates up to several g s inferred in some systems has not been clear for a long time.
Recently Rafikov (2011a; hereafter R11) showed that coupling of the disk to stellar radiation via the Poynting-Robertson (PR) effect (Burns et al. 1979) can result in metal accretion rates g s. Later Rafikov (2011b) has also demonstrated that under certain circumstances even higher values of can arise from the interaction between the debris disk and metal gas that is produced by the sublimation at .
The model of radiatively driven accretion of R11 did not cover the large scale evolution of the debris disk resulting from the PR drag, and made certain assumptions (e.g. high optical depths of the disk) which were not rigorously verified. The goal of this work is to extend the analysis of R11 and to develop a detailed global model of the compact debris disk evolution caused by the PR drag.
The paper is organized as follows. In §2 we outline the basic picture of the PR-driven debris accretion and derive master equation (17) that describes global evolution of the disk. We then explore in §3 both analytically and numerically the evolution of a low mass, optically thin disk of debris. In §4 we study global evolution of massive, optically thick debris disks starting with different initial spatial distributions of debris around the WD (ring-like, §4.1.1, 4.1.2, or disk-like, §4.1.3) to see the effect on the global disk evolution. We discuss our results and their observational implications in §5.
2. Description of the model.
In the following we consider an axisymmetric disk of particles extending from to in radius. The inner radius may coincide with the sublimation radius , at which the effective temperature of particles equals the sublimation temperature :
where is the WD radius, K for silicate grains, and K) is the normalized stellar temperature . Taking one finds R, in agreement with observationally inferred inner radii of compact debris disks (Jura et al. 2007, 2009a).
When particles reach the sublimation radius they produce metallic gas, which joins the gaseous disk extending down to the WD surface. Metal accretion onto the WD surface proceeds through this disk. Although this metal gas also spreads outward from the sublimation radius (Melis et al. 2010) and under certain circumstances can substantially affect debris disk evolution (Rafikov 2011b), in this work, following R11, we concentrate only on effects associated with the PR drag. Thus, here we neglect presence of the gas exterior of and its interaction with the debris disk. According to Rafikov (2011b) this is a valid approximation as long as the viscous timescale in the gaseous disk is shorter than the time on which this disk can be replenished by the sublimation of debris. Then the gas does not accumulate at and its density is always low enough for the gas drag to not affect the debris disk appreciably.
We characterize the disk at each point by its surface density and optical depth
where is the bulk density of particles and is their characteristic size. The particle size is not well constrained by the existing observations, but its actual value becomes important only for low mass disks with small , see §3.
Following Friedjung (1985) and R11 we represent irradiation of the disk by the WD using a single incidence angle at each radius from the disk (the so-called “lamp-post” illumination model described in R11):
where is the radius of the star. We can then introduce optical depth to incoming starlight according to the following definition:
This variable is very important for our subsequent analysis. In the following we will call the disk optically thick or thin based on whether and not is greater or smaller than unity.
2.1. Accretion through the disk.
R11 demonstrated that the mass accretion rate through the disk of solids driven by the PR drag may be written as
Here is the WD luminosity, is the speed of light and function
gives the fraction of incoming starlight that is absorbed by the disk.
It is instructive to examine different limits of the expression (5). If the disk is optically thick to incident stellar radiation, (which cannot be the case everywhere at all times, as we will see below), then the accretion rate is independent of and is given simply by
In particular, at the sublimation radius given by equation (1) this accretion rate is equal to
R11 has shown this numerical estimate to be consistent with the lower envelope of the values inferred for the metal-rich WDs exhibiting detectable IR emission associated with their debris disks.
In the opposite limit of a disk, which is optically thin to incident stellar radiation, , one finds
It is easy to show that this expression coincides with that one would calculate by simply assuming disk particles to be directly illuminated by starlight (i.e. unobscured by other particles) and considering each of them as independent from others. Azimuthal PR drag force acting on a single, perfectly absorbing particle of mass is given by
where is Keplerian angular velocity. This force gives rise to radial drift velocity
and results in mass accretion rate , which is easily shown to reduce to equation (10).
2.2. Evolution equations
Evolution of the debris disk is described by the continuity equation
with given by equation (5).
We introduce new dimensionless time and space variables
is the characteristic timescale of the problem. If coincides with the sublimation radius defined in (1) then
where we took g cm.
A notable fact about equation (17) is that it does not contain any free parameters.
Once is obtained by solving this equation, one can easily infer mass accretion rate at any point in the disk since
where is defined by equation (8). Because the viscous timescale of the gas produced by debris sublimation is very short compared to the disk evolution time (R11), the metal accretion rate onto the WD surface is given simply by
In the following we will use and interchangingly.
It has an implicit solution
where is an arbitrary function, which is set by the initial surface density distribution of the debris. Unfortunately, inferring from is not a trivial task, which makes analysis of solution (21) very difficult. Nevertheless, we can still deduce useful information about the behavior of this solution in some specific limits, see §3 and Appendix A.
We now proceed to investigate the details of the global disk evolution with different initial conditions.
3. Low mass disks.
We start by considering the case of a tenuous disk, in which the optical depth for incident stellar radiation . In that case according to equation (18) and the nonlinear term can be neglected in equation (20) reducing it to
This linear equation has an explicit solution
where is the initial spatial distribution of variable set by the initial surface density distribution .
If initially the debris was concentrated in a structure (a disk or a ring) with characteristic dimensionless scale (physical scale ), then it follows from equation (22) that the characteristic time on which evolves is
This time is independent of the disk mass but is sensitive to its size and particle properties, see equation (16). This is not surprising since when particles interact with the WD radiation independently of each other and drift towards the WD on a PR timescale of an individual particle. The latter coincides with up to constant factors and is indeed a function of , and the distance that needs to be traveled.
Disks which are optically thin to incident stellar radiation must satisfy a certain mass constraint. The mass of a disk extending out to can be written as
where the characteristic disk mass is defined as
and the numerical estimate assumes g cm, and . This numerical value of corresponds to the disruption of an asteroid with the diameter of 40 km.
The disk which is optically thin everywhere () should have mass since is of order unity for WD debris disks.
3.1. Numerical results.
We numerically integrated equation (17) for an initial ring-like density profile given by
with . In Figure 1 we show the results obtained for the case , , , i.e. a rather narrow, tenuous ring, with the radius large compared to . We also plot analytical solution for
which follows from (23). One can see that analytical solution reproduces the numerical results quite well. As time goes by the ring of material drifts towards the WD and broadens, simply because the speed of inward migration due to the PR drag increases as decreases. This eventually results in substantial surface density reaching the inner radius and giving rise to sublimation there. Metal gas produced at then directly accretes onto the WD.
Metal accretion rate in the optically thin case, simply because the initial value of everywhere in the disk and equation (23) then predicts that for any . Using solution (23) we can derive an explicit formula for for the initial debris density distribution (27):
In Figure 2 we compare this expression with the numerically derived and again the agreement is reasonably good. Mass accretion rate peaks at when the expression inside the exponential in (29) goes to zero. The maximum accretion rate achieved at this moment is about , in agreement the numerical results shown in Figure 2.
4. High mass disks.
Now we consider evolution of massive debris disks with mass exceeding defined by (26). Such disks have at least at some radii. This condition cannot be true throughout the whole disk as there should always be regions in which , e.g. near the disk edges. Even if they do not exist initially such optically thin regions naturally appear as a result of disk evolution as we show below.
Whenever one can neglect the term proportional to in equation (17) and obtain the following simple solution:
or, in dimensional units,
This solution implies that in parts of the disk, which are optically thick to incoming stellar radiation the surface density steadily decreases in time at a constant rate set only by the distance from the WD. This behavior is a direct consequence of the saturation at the level independent of in the case , see equation (7).
At a glance this kind of evolution looks very strange as the debris does not seem to move anywhere — simply goes down at all radii. This interesting behavior can be understood only by accounting for the optically thin regions that should inevitably be present in the disk and their interplay with the optically thick part of the disk. We demonstrate this in §4.1.
Solution (30) also implies that the characteristic timescale of the disk evolution is now
or in dimensionless units. Clearly, is (up to factors of order unity) just the time, on which the characteristic disk mass gets exhausted by accretion at the rate . Unlike defined by equation (24) is independent of particle properties — and , but scales linearly with the disk mass.
4.1. Numerical results.
We numerically integrated evolution of several initially optically thick density distributions. The morphologies we consider range from a narrow ring to an extended disk of debris. They represent different initial spatial distributions of the debris that might possibly result from the tidal disruption of an asteroid. We choose all disks to have the same total mass initially and concentrate on exploring the differences of their evolutionary routes caused purely by the morphology.
Large and narrow ring.
We start by exploring the initial distribution in the form (27) with , but now we take so that the disk mass . Resulting evolution of is shown in Figure 3. Apart from the theoretically expected steady decrease of with time in the bulk of the disk where one immediately notices two interesting features of the solutions: (1) a sharp drop of the disk surface density at the outer edge of the disk, and (2) an optically thin () tail of debris extending from the optically thick region all the way to , which remains stable throughout the evolution of the disk. We now discuss the nature of these features in more detail.
Sharp outer edge appears because the evolution timescale of a particular disk region , see equation (32). Initial density distribution in the form (27) has increasing inwards in the outer part of the ring. This causes the outermost extremity of the ring (which is optically thin and evolves on timescale ) to migrate towards the WD faster than the parts of the ring closer to the WD do. As a result, the outermost disk material catches up with the debris at smaller having higher , leading to a debris pileup there and formation of a sharp outer edge. This edge appears on a timescale since this is the characteristic time, on which the material from the optically thin outer disk region reaches where most of the ring mass is concentrated. We provide a more rigorous justification for the outer edge appearance in Appendix A.
Origin of the quasi-stationary tail with at small can be understood as follows. Even if initially the disk has everywhere such tail would rapidly form since, according to equation (30), at any distance from the WD should switch into the optically thin regime within finite time. This happens fastest at the inner edge of the disk. After the tail has formed it evolves on a timescale given by equation (24). Inner parts of the tail evolve more rapidly than the outer ones and a quasi-steady state is attained in the tail with almost independent of . This point is illustrated in Figure 4, where we show the run of for a calculation displayed in Figure 3 at different moments of time. One can easily see that after some evolution has taken place indeed becomes independent of in the tail.
Thus, the behavior of in the tail can be described by the equation (cf. equation (18))
where is the mass accretion rate through the tail (and is the dimensionless analog of this quantity), which in general is a function of time. Equation (33) describes the behavior of in Figure 3 as a function of quite well.
At some radius the tail becomes optically thick, i.e. . A this point the tail merges with the optically thick part of the disk and the evolution timescale, which is now given by , rapidly increases, see equation (32). Accretion rate through the tail is set almost exclusively by the value of since by definition independent of , so that
at all times. Thus, it is the dependence of on that determines how the mass accretion rate through the tail and at changes with time.
To find we need to set the left hand side of equation (30) to unity and solve the resultant equation for as a function of . Since here we consider the case of the relation between and is provided with sufficient accuracy by solving the equation
for a given initial density distribution in the disk.
For the Gaussian initial distribution in the form (27) one finds that
Whenever and , which is a reasonable approximation for the case shown in Figure 3, equation (36) predicts that and that depends on time only weakly (logarithmically). This explains the time invariance of the optically thin tail in Figure 3 and only weak time dependence of in the tail in Figure 4: the mass flux through the tail
is almost independent of if the initial mass distribution has the form of a large and narrow ring (i.e. ). As a result, the flux of metals onto the WD also varies with time only weakly for ring-like initial distribution of the debris, as can be seen in Figure 5. For debris disks around WDs produced by tidal disruption of asteroids and so equation (37) guarantees that stays above or so, see Figure 5.
Equation (33) also predicts that does not go significantly below unity even at the inner edge of the disk: depends only on , see equation (33), and does not change appreciably with time for a narrow ring.
It is clear that the general evolutionary picture presented in Figure 3 would still hold if we assume a ring density profile other than Gaussian, as long as the width of the ring is small and its radius is large compared to .
We now consider evolution of initial distribution in the Gaussian form (27) but with , , and truncated at (i.e. for ), while keeping the total disk mass the same, . In this case the width of the ring is no longer small compared to its radius. Snapshots of the profile at different time are presented in Figure 6.
Beyond this point a marginally optically thin tail develops near but it never reaches a steady state. This can be understood by inspection of the transcendental equation (36) determining the crossover radius (and the through the tail): whenever the width of the ring is not small compared to its radius this equation does not have a time-independent solution. As a result, noticeably varies in time, see Figure 5.
Finally we explore the case of an extended initial density distribution in the form of a disk described by a “top-hat” initial condition
where and . Evolution of this distribution is shown in Figure 7 and one can again clearly see the development of an optically thin tail at the inner edge of the disk and the sharp outer edge.
The tail appears on a timescale . Initially this tail attaches directly to because it takes time (see equation (35)) for to go down to unity. During this period of time stays virtually constant. These simple analytical conclusions are in excellent agreement with the numerical results presented in Figure 5.
This analytical prediction is shown by dotted line in Figure 5 and clearly agrees well with the numerical result for the initial top-hat density distribution.
In this paper we focused on studying how the debris disk evolves due to the PR drag alone. Among other radiative effects potentially affecting such disks we can mention the direct transfer of the radiative angular momentum to the disk that occurs if the WD is rapidly rotating. Previously Walker & Meszaros (1989) suggested this process to be important for the evolution of accretion disks around neutron stars experiencing Type I X-ray bursts, which often rotate quite rapidly. In the case of WDs there is no good evidence for rapid rotation (Kawaler 2004) with several disk-hosting systems showing km s (Gänsicke et al. 2006, 2007). However, one can easily demonstrate that even if WD were to rotate at breakup speed the radiative angular momentum transport would still not produce WD disk evolution comparable to that caused by the PR drag.
Indeed, the angular momentum carried by photons emitted by the WD rotating with equatorial speed can be estimated as . Disk absorbs some fraction of this angular momentum which exerts radiative torque on the disk. The ratio of this torque to the torque due to the PR drag is easily shown to be
where is the Keplerian angular frequency. Even if the WD rotates at breakup and one still finds , which is less than unity because the debris disk lies quite far from the WD surface, i.e. . Moreover, measured rotational speeds of the aforementioned WDs are of their breakup speed, so the direct radiative transfer of angular momentum to debris has negligible effect on disks in these systems.
One of the key results of this work is that the low metal accretion rates onto the WD are possible only if the debris disk hosted by the WD has small mass, below given by equation (26), and is everywhere optically thin to incident stellar radiation. This is demonstrated in §3 where we find (see equation (23)) throughout the whole evolution of the low mass disk.
On the contrary, if the disk is massive enough to contain optically thick regions, which typically requires g, then is guaranteed to not fall below , as we have shown in §4. This is in agreement with the statement in R11 that should not deviate significantly from for optically thick disks.
This lower bound on for massive disks is very robust, in particular it does not depend on the disk mass. This point is illustrated in Figure 8 where we show the time evolution of for three disks with the same initial Gaussian density distribution (27) but with different initial masses set to be in the ratio . One can easily see that the only thing, which is different between these three cases is the disk lifetime (which is always proportional to the disk mass), while is essentially the same.
The time evolution of in the optically thick case is found to depend predominantly on the shape of the initial density distribution. Debris distributions in the form of a narrow ring tend to produce only weakly varying with time. Extended initial distributions of the debris with radial width comparable to mean radius generally lead to varying in time and decaying by a factor of several by the end of the disk lifetime, see Figure 5.
Based on these results one should expect WDs with strong IR excesses, signifying massive, optically thick debris disk around them, to exhibit . This is exactly what the current data seem to suggest: no WD with a debris disk detected via its IR signature is inferred to have below (R11).
However, if the systems exhibiting IR excesses were found to show , the debris disks around them should be optically thin. It is then natural to expect a positive correlation between the strength of the IR excess and the in such systems, since both are linearly proportional to in the optically thin regime. Indeed, according to (23) while the amount of the WD luminosity that the optically thin disk intercepts and reradiates in the near-IR is
which also scales positively with .
The timescale on which debris disks get depleted can be quite different in the optically thick and thin cases. Characteristic lifetime of an optically thin disk composed of cm particles is several yr, as equations (16) and (24) demonstrate. This timescale depends only on the debris particle properties — size and density, being longer for bigger and denser particles — and the outer extent of the disk (for one finds disk lifetime approaching a Myr).
In the case of an optically thick disk the lifetime can be considerably longer since scales linearly with the disk mass, see equation (32). Thus, massive disks require long time to get exhausted by the PR drag alone. For example, an optically thick disk created by the tidal disruption of an asteroid with the diameter of 300 km would have a mass of g (assuming bulk density of 3 g cm). Accretion at the steady rate g s ( given by the numerical estimate in equation (8)) would deplete this disk only within 40 Myr. Such massive disks, however, are likely to be affected by the interaction with the gas that is produced by sublimation of solid debris at the inner edge of the disk. They may then evolve in a runaway fashion as described in Rafikov (2011b) and get exhausted much earlier.
An interesting finding of this work is the universal appearance of the sharp outer edge in optically thick disks, caused by the faster inward migration of the disk annuli having smaller optical depth. This outer edge steadily moves inwards which, combined with the outward expansion of the optically thin tail at small radii, leads to the gradual reduction of the radial width of the optically thick part of the debris disk. Thus, even an extended disk can turn into a narrow ring with time as can be easily seen in Figure 7.
Observational evidence for sharp outer edges would strongly support the picture of disk evolution due to the PR drag outlined in this work. Interestingly, Jura et al. (2007b, 2009) have previously claimed the need for existence of the outer optically thin regions in debris disks around some WDs. However, inferring the sharpness of the outer disk edge and the precise value of its optical depth is a nontrivial task since any knowledge about the radial density distribution in these disks is based on modeling their spectral energy distributions. Inverting disk spectrum to determine the fine details of the radial distribution of debris is a rather ambiguous procedure, additionally complicated by the possibility of disk warping or flaring (Jura et al. 2007b, 2009).
Appendix A Formation of sharp outer edge.
Equation (20) is similar to Hopf (or simple wave) equation , which is known to result in appearance of multivalued solutions usually interpreted as evidence for the shock formation. We may expect something similar in our case as well because the nonlinearity is clearly present in equation (20). The equation for characteristic, crossing the point (), is
Since the left hand side of (A1) is negative for all , characteristic does not exist at times
From that we can derive an upper limit on the time at which formation of a sharp edge occurs:
Thus, we expect a discontinuity in surface density to form during disk evolution and our numerical results in §4.1 support this expectation.
- affiliation: Departmen of General and Applied Physics, Moscow Institute of Physics and Technology, Dolgoprudny, 141700, Russia; firstname.lastname@example.org
- affiliation: Department of Astrophysical Sciences, Princeton University, Ivy Lane, Princeton, NJ 08540; email@example.com
- affiliation: Sloan Fellow
- Brinkworth, C. S., Gänsicke, B. T., Marsh, T. R., Hoard, D. W., & Tappert, C. 2009, ApJ, 696, 1402
- Burns, J. A., Lamy, P. L., & Soter, S. 1979, Icarus, 40, 1
- Cuzzi, J. N. et al. 2010, Science, 327, 1470
- Dufour, P., Kilic, M., Fontaine, G., Bergeron, P., Lachapelle, F.-R., Kleinman, S. J., & Leggett, S. K. 2010, ApJ, 719, 803
- Dupuis, J., Fontaine, G., Pelletier, C., & Wesemael, F. 1993a, ApJS, 84, 73
- Dupuis, J., Fontaine, G., & Wesemael, F. 1993b, ApJS, 87, 345
- Ehrenreich, D. et al. 2011, A&A, 525, 85
- Farihi, J., Zuckerman, B., & Becklin, E. E. 2008, ApJ, 674, 431
- Farihi, J., Jura, M., & Zuckerman, B. 2009, ApJ, 694, 805
- Farihi, J., Jura, M., Lee, J.-E., & Zuckerman, B. 2010, ApJ, 714, 1386
- Friedjung, M. 1985, A&A, 146, 366
- Gänsicke, B. T., Marsh, T. R., Southworth, J., & Rebassa-Mansergas, A. 2006, Science, 314, 1908
- Gänsicke, B. T., Marsh, T. R., & Southworth, J. 2007, MNRAS, 380, L35
- Graham, J. R., Matthews, K., Neugebauer, G., & Soifer, B. T. 1990, ApJ, 357, 216
- Jura, M. 2003, ApJL, 584, L91
- Jura, M., Farihi, J., & Zuckerman, B. 2007a, ApJ, 663, 1285
- Jura, M., Farihi, J., Zuckerman, B., & Becklin, E. E. 2007b, AJ, 133, 1927
- Jura, M., Farihi, J., & Zuckerman, B. 2009a, AJ, 137, 3191
- Jura, M., Muno, M. P., Farihi, J., & Zuckerman, B. 2009b, ApJ, 699, 1473
- Kawaler, S. D. 2004, in Stellar Rotation, Proc. IAU Symp., 215, ed. A. Maeder & P. Eenes, 561
- Melis, C., Jura, M., Albert, L., Klein, B., & Zuckerman, B. 2010, ApJ, 722, 1078
- Rafikov, R. R. 2011a, ApJ, 732, L3 (R11)
- Rafikov, R. R. 2011b, arXiv:1102.4343
- Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337
- Zuckerman, B. & Becklin, E. E. 1987, Nature, 330, 138
- Zuckerman, B., Koester, D., Melis, C., Hansen, B. M. S., & Jura, M. 2007, ApJ, 671, 872
- Walker, M. A. & Meszaros, P. 1989, ApJ, 346, 844