# Radiation-Driven Warping of Circumbinary Disks Around Eccentric Young Star Binaries

###### Abstract

We study a warping instability of a geometrically thin, non-self-gravitating, circumbinary disk around young binary stars on an eccentric orbit. Such a disk is subject to both the tidal torques due to a time-dependent binary potential and the radiative torques due to radiation emitted from each star. The tilt angle between the circumbinary disk plane and the binary orbital plane is assumed to be very small. We find that there is a radius within/beyond which the circumbinary disk is unstable to radiation-driven warping, depending on the disk density and temperature gradient indices. This marginally stable warping radius is very sensitive to viscosity parameters, a fiducial disk radius and the temperature measured there, the stellar luminosity, and the disk surface density at a radius where the disk changes from the optically thick to thin for the irradiation from the central stars. On the other hand, it is insensitive to the orbital eccentricity and binary irradiation parameter, which is a function of the binary mass ratio and luminosity of each star. Since the tidal torques can suppress the warping in the inner part of the circumbinary disk, the disk starts to be warped in the outer part. While the circumbinary disks are most likely to be subject to the radiation-driven warping on a AU to kilo-AU scale for binaries with young massive stars more luminous than , the radiation driven warping does not work for those around young binaries with the luminosity comparable to the solar luminosity.

^{†}

^{†}slugcomment:

## 1 Introduction

About 60 of main sequence stars are considered to be born as binary or multiple systems (Duquennoy & Mayor, 1991). Indeed, the direct imaging of the circumbinary disk was successfully done for a few young binary star systems with the binary mass comparable to a solar mass, including GG Tau (Dutrey et al., 1994) and UY Aurigae (Duvert et al., 1998) by the Plateau de Bure Interferometer (PdBI). It is confirmed by numerical simulations that the circumbinary disk is formed around such young low mass binary stars embedded in the dense molecular gas (Artymowicz & Lubow, 1994; Bate & Bonnell, 1997; Günther & Kley, 2002, 2004).

The existence of a significant number of OB eclipsing binaries (Hilditch et al., 2005) suggests that most known massive stars are also born as binary or multiple systems (Sana et al., 2012). Such binaries are expected to be accompanied by circumbinary disks at the early stage of the massive star formation. Contrary to the case of T Tauri stars, however, the massive star formation is still poorly understood because the observation of massive star forming regions is challenging. This is because they are obscured by a dusty environment and the chance to observe massive young stellar objects is small because of their short lives. Araya et al. (2010) suggested that observed periodic maser emissions are originated from periodic accretion onto the binary stars on a significantly eccentric orbit from the circumbinary disk (Artymowicz & Lubow, 1996), although there are two other scenarios: the colliding wind binary scenario (van der Walt, 2011) and the scenario of pulsation of protostars growing via high mass accretion rate (Inayoshi et al., 2013).

It is natural that the circumstellar and circumbinary disks in the star forming regions have a warped structure. Krist et al. (2005) suggested that there is observational evidence for the disk warping in GG Tau system. The tidal effects of close encounters between stars and a disk surrounding it can make the disk outer parts misaligned or warped (Moeckel&Bally, 2006). If the circumbinary disk is originally misaligned with the binary orbital plane, it should be warped by the tidal alignment, although the disk tilt secularly decays by the disk viscosity (Facchini et al., 2013; Lodato&Facchini, 2013). The origin of the misalignment and disk warping is, however, still a matter of debate.

Hayasaki et al. (2014a) (hereafter, Paper I) derived the condition that a circumbinary disk subject to both the tidal torque and the torque due to the radiation emitted from two accretion disks around the individual supermassive black holes is unstable to the radiation-driven warping in the context of observed warped maser disks in active galactic nuclei (e.g., Kuo et al. 2011; Kormendy&Ho 2013). They assumed that the binary is on a circular orbit and the irradiation luminosity is proportional to the mass accretion rate. There is, however, observational evidence that young binaries have a significant orbital eccentricity (Mathieu et al., 1990). In addition, the circumbinary disk is mainly irradiated by two young stars in a binary in the current problem. In order to apply our previous model to those binaries, therefore, we need to relax the above assumptions adopted in Paper I and make the binary mass scale down to the stellar mass from supermassive black hole mass.

In this paper, we study the warping instability of a circumbinary disk around young binary stars on an eccentric orbit. In section 2, we describe the external torques acting on a circumbinary disk. We consider both the tidal torques originating from a time-dependent binary potential and the radiative torques from the binary stars. In section 3, we examine the evolution of a slightly tilted circumbinary disk subject to those two torques, and derive the warping condition and timescale of the local precession of the linear warping mode. Finally, section 4 is devoted to summary and discussion of our scenario.

## 2 External Torques acting on the circumbinary disk

Let us consider the torques from the binary potential acting on the circumbinary disk surrounding two stars in a binary on an eccentric orbit. Figure 1 illustrates a schematic picture of the setting of our model; binary stars orbiting each other are surrounded by a misaligned circumbinary disk. The binary is put on the - plane with its center of mass being at the origin in the Cartesian coordinate. The masses of the primary and secondary stars are represented by and , respectively, and . We put a circumbinary disk around the origin. The unit vector of specific angular momentum of the circumbinary disk is expressed by (e.g. Pringle 1996)

(1) |

where is the tilt angle between the circumbinary disk plane and the binary orbital plane, and is the azimuth of tilt. Here, , , and are unit vectors in the , , and , respectively. The position vector of the circumbinary disk can be expressed by

(2) |

where the azimuthal angle is measured from the descending node. The only difference from the Paper I is the position vector of each star, which is given by

(3) |

where is the true anomaly and is written as

(4) |

with and . Here, is the binary orbital eccentricity, is the binary mass ratio, and is the binary semi-major axis with and . These and other model parameters are listed in Table 1.

Definition | Symbol |
---|---|

Total stellar mass | |

Primary stellar mass | |

Secondary stellar mass | |

Binary mass ratio | |

Mass ratio parameters | , |

Binary semi-major axis | |

Orbital eccentricity | |

Orbital frequency | |

Orbital period | |

True anomaly | |

Tilt angle | |

Azimuth of tilt | |

Azimuthal angle | |

Shakura-Sunyaev viscosity parameter | |

Horizontal shear viscosity | |

Vertical shear viscosity | |

Ratio of vertical to horizontal shear viscosities | |

Luminosities emitted from two stars in a binary | , |

Total luminosity | |

Binary irradiation parameter |

### 2.1 Gravitational Torques

The gravitational force on the unit mass at position on the circumbinary disk can be written by

(5) |

The corresponding torque is given by

(6) |

We consider the tidal warping/precession with timescales much longer than local rotation period of the circumbinary disk. This allows us to use the torque averaged in the azimuthal direction and over the orbital period:

(7) | |||||

where is the angular frequency of mean binary motion. Here, we used for the integration the following relationship:

(8) |

and the approximations:

(9) |

Note that equation (7) is derived in the same manner as in equation (7) of Hayasaki et al. (2014b).

For a small tilt angle , equation (7) is reduced to

(10) |

where and can be written from equation (1) as and . The magnitude of the specific tidal torque with is then given by

(11) |

The tidal torque tends to align the tilted circumbinary disk with the orbital plane (c.f. Bate et al. 2000). The tidal alignment timescale for is given by

(12) |

where with the disk angular frequency and are the specific angular momentum and binary orbital period, respectively. Equation (12) is reduced to equation (8) of Paper I for . Note that the tidal alignment timescale depends on the azimuth of tilt for . Since the inner edge of the circumbinary disk is estimated to be (Artymowicz & Lubow, 1994), the tidal alignment timescale is longer than the binary orbital period.

### 2.2 Radiative Torques

The circumbinary disk around young binary stars can be mainly illuminated by light emitted from each star. The re-radiation from the circumbinary-disk surface, which absorbs photons emitted from these stars, causes a reaction force. This is the origin of the radiative torques. Here, the central two stars are can be regarded as point irradiation sources, because each star is much smaller than the size of circumbinary disk. Circumstellar disks in binary systems have been observed by many authors in the past (e.g. Mayama et al. 2010). These circumstellar disks could be other irradiation sources. We estimate their luminosities as follows: the maximum luminosity of each disk is given by with , where the Eddington luminosity is given by . Since it is substantially smaller than the stellar luminosity, the effect of the circumstellar-disk irradiation is negligible.

Since the surface element on the circumbinary disk is given in the polar coordinates by

(13) |

the radiative flux at is given by

(14) |

where is the sum of luminosities of the radiation emitted from the primary star, , and that from the secondary star, . Here, we assume that the surface element is not shadowed by other interior parts of the circumbinary disk. If we ignore limb darkening, the force acting on the disk surface by the radiation reaction has the magnitude of and is antiparallel to the local disk normal (Pringle, 1996). The total radiative force on can then be written by

(15) |

Consequently, the total radiative torque acting on a ring of radial width is given by

(16) |

where holds for and , and the first term, which we call , of the right-hand side of equation (16) corresponds to equation (2.15) of Pringle (1996):

(17) |

and the second term, which we call , is originated from the orbital motion of the binary.

Here, we consider the radiation-driven warping/precession with timescales much longer than the orbital period, as in the case of tidally driven warping/precession. The orbit-average of the torque is then given by

(18) | |||||

where is a binary irradiation parameter ( by definition), and is reduced to for .

From equation (18), the specific radiative torque averaged over azimuthal angle and orbital phase is given by

(19) | |||||

where and are the growth speed of a warping mode induced by the radiative torque and disk surface density, respectively. The growth timescale of the warping mode for an optically thick gas disk can be estimated to be

(20) |

Since it is clear that the growth timescale for is much longer that the orbital period, our assumption for the orbit-averaged radiative torque is ensured.

## 3 Tilt angle evolution of circumbinary disks

In this section, we examine the response of the circumbinary disk to the external forces for case. The evolution equation for disk tilt is given by(Pringle, 1996)

(21) |

where and are respectively the horizontal and vertical shear viscosities, the latter of which tends to reduce disk tilt, and is the term to which the external torques contribute. For a geometrically thin disk, is approximately given by with the Shakura-Sunyaev viscosity parameter (Shakura & Sunyaev, 1973), the isothermal sound speed , the temperature , the gas constant , and the molecular weight . The primes indicate differentiation with respect to . We adopt for the circumbinary disk structure the following assumptions that ,

(22) | |||||

(23) |

where and are a constant, and , , and are the fiducial radius, fiducial surface density, and fiducial temperature, respectively. Equation (21) can be then reduced to

(24) |

where is used. Here, with the ratio of vertical to horizontal viscosities: and is written as

(25) | |||||

We look for solutions of equation (24) of the form , with . Replacing with , with , and with , we have the following set of linearized equations:

(26) |

where

(27) |

and, represents the magnitude of the local precession frequency of the liner warping mode for the case of the standard disk model and (see equation (35) of Paper I):

(28) |

The determinant of the coefficient matrix on the left hand side of equation (26) must vanish because of . The local dispersion relation is then obtained as

(29) |

where

(30) |

Here, we adopt the approximation that the terms proportional to or are negligible in comparison with the other terms. The dispersion relation is then rewritten in the following simple form:

(31) |

The imaginary part of corresponds to the excitation or damping of oscillation, whereas the real part provides the local precession frequency due to the external torques.

In order for the perturbation to grow, must be negative. The growth condition is given by

(32) |

In terms of , the growth timescale of the warping mode induced by the radiative torques in the binary system is given by

(33) |

Figure 2 shows the dependence of on the orbital eccentricity and in a circumbinary disk with , , , , and . Each growth timescale is normalized by . The growth timescale is longer than that of circular binary case in the range of . This is because the averaged incident flux normal to the disk surface per binary orbit is lower than that of the circular binary case. Note that is approximately equal to for . We therefore treat as the growth timescale of the binary star case in what follows.

We focus our attention on a perturbation with , where is the radial wavelength of the perturbation. The condition that the circumbinary disk is unstable to the warping mode can be, then, rewritten as

(34) |

Substituting equations (22) and (23) into the above equation, we obtain

(35) |

where shows the marginally stable warping radius for a single star case:

(36) |

where in deriving the second equation, we have used the relationship between and : for (Ogilvie, 1999). In the case of , no unstable solution exists except for very special combination of parameters.

The radiative torques work only for a region of optically thick to irradiation emitted from the central stars. In order for the circumbinary disk to be optically thick, the surface density must be higher than , where the electron scattering opacity is assumed for simplicity. This condition is rewritten as

(37) |

using equation (22). The equality is held at the radius where the circumbinary disk changes from optically thick to optically thin:

(38) |

In order for radiation-driven warping to be a possible mechanism for disk warping, the marginally stable warping radius must be less than the outer radius of the optically thick region, . This gives the upper (lower) limit of for as follows:

(39) |

where is the critical surface density given by

(40) |

Figure 3 shows the possible range in , for the circumbinary disk to be unstable to warping modes, as a function of and . The upper two panels give the range for case, whereas the lower two panels for case. From the figure, we can express as

(41) |

with a factor in the range for and for . Substituting equation (41) into equation (36) with equation (35), we can obtain the marginally stable warping radius as

(42) |

where can be rewritten as

(43) | |||||

Adopting for equation (43) , which corresponds to the power law index of the radial temperature profile of an optically thick region in the standard disk (Shakura & Sunyaev, 1973), the marginally stable warping radius is very sensitive to the values of , , , , and . It is clear from equations (42) and (43) that a circumbinary disk around young binary stars with and is stable for radiation driven warping, because the marginally stable warping radius is of the order of for and of even for . This shows that the circumbinary disks of a classical T Tauri star system such as GG Tau (Dutrey et al., 1994) and UY Aurigae (Duvert et al., 1998) are not warped by radiation-driven warping instability.

On the other hand, it is likely that a circumbinary disk around binary stars with much larger luminosities is unstable to radiation-driven warping. Figure 4 shows the dependence of the marginally stable warping radius in such a disk on the semi-major axis. Here, we adopt , , , , , , , , , and . While panels (a) and (c) are for , panels (b) and (d) are for . We also adopt and in panels (a) and (b), and and in panels (c) and (d). In these four panels, the black solid line and orange dashed line show and in units of , respectively. The red dotted line shows the radius where the growth timescale of the radiation-driven warping mode, , equals the timescale for the disk to align with the orbital plane by the tidal torque, . This tidal alignment radius is given by

(44) |

from equations (12), (20) and (22). The growth of a finite-amplitude warping mode induced by the radiative torque can be significantly suppressed by the tidal torque in the region inside the tidal alignment radius. The blue dot-dashed line show the inner radius of the circumbinary disk. Here, the inner radius is assumed to be equal to the tidal/resonant truncation radius, where the tidal/resonant torque is balanced with the horizontal viscous torque of the circumbinary disk (Artymowicz & Lubow, 1994). In the case of eccentric binaries with a nearly equal mass ratio, the tidal/resonant truncation radius is estimated to be . The shaded area shows the region where the circumbinary disk is unstable to the warping modes induced by the radiative torques.

From the figure, the marginally stable warping radius is on a sub-AU to kilo-AU scale for appropriate model parameters. In panes (a) and (b), the circumbinary disk is warped in the region outside all of the marginally stable warping radius, the tidally alignment radius, and the inner disk radius, and inside the outer radius of the optically thick part (in the region above the black line, below the orange dashed line, and on the left-side of both the red dotted line and the blue dot-dashed line). On the other hand, in panels (c) and (d), the circumbinary disk is warped in the region inside both the marginally stable warping radius and the outer radius of the optically thick region, and outside both the tidal alignment radius and the inner disk radius (in the region below the black line and the orange dashed line, and on the left-side of both the red dotted line and the blue dot-dashed line). Note that in any case the middle part of the circumbinary disk is warped because the warping mode in the inner part is suppressed by the tidal torques, and the radiative torques are effective only within the outer radius of the optically thick region.

The local precession frequency of the linear warping mode is obtained from the dispersion relation by

(45) |

where

(46) | |||||

(47) | |||||

(48) |

where , , are the local precession frequencies due to the tidal, viscous, and radiative torques, respectively. Note that has the same expression with the opposite sign to the precession frequency of the tidally induced wave (see equation (11) of Hayasaki & Okazaki 2009). This satisfies the condition for simultaneous resonant excitation of eccentric and tilt waves proposed by Kato (2014). Therefore, both waves are likely to coexist in the circumbinary disk.

Adopting from equation (32), we obtain