Inclination Evolution of Protoplanetary Disks Around Eccentric Binaries

It is usually thought that viscous torque works to align a circumbinary disk with the binary's orbital plane. However, recent numerical simulations suggest that the disk may evolve to a configuration perpendicular to the binary orbit ("polar alignment") if the binary is eccentric and the initial disk-binary inclination is sufficiently large. We carry out a theoretical study on the long-term evolution of inclined disks around eccentric binaries, calculating the disk warp profile and dissipative torque acting on the disk. For disks with aspect ratio $H/r$ larger than the viscosity parameter $\alpha$, bending wave propagation effectively makes the disk precess as a quasi-rigid body, while viscosity acts on the disk warp and twist to drive secular evolution of the disk-binary inclination. We derive a simple analytic criterion (in terms of the binary eccentricity and initial disk orientation) for the disk to evolve toward polar alignment with the eccentric binary. When the disk has a non-negligible angular momentum compared to the binary, the final"polar alignment"inclination angle is reduced from $90^\circ$. For typical protoplanetary disk parameters, the timescale of the inclination evolution is shorter than the disk lifetime, suggesting that highly-inclined disks and planets may exist orbiting eccentric binaries.


INTRODUCTION
To date, 11 transiting circumbinary planets have been detected around 9 binary star systems (Doyle et al. 2011;Kostov et al. 2013Kostov et al. , 2014Kostov et al. , 2016Orosz et al. 2012a,b;Schwamb et al. 2013;Welsh et al. 2012Welsh et al. , 2015. All planets detected have orbital planes very well aligned with their binary orbital planes, with mutual binary-planet inclinations not exceeding 3 • . The circumbinary planet detectability is a very sensitive function of the binary-planet inclination (Martin & Triaud 2015;Li, Holman, & Tao 2016). If the mutual inclination is always small ( 5 • ), then the occurance rate of circumbinary planets is comperable to that of planets around single stars, but if modest inclinations ( 5 • ) are common, the circumbinary planet occurance rate may be much larger (Armstrong et al. 2014). For these reasons, it is important to understand if and how a binary aligns with its circumbinary disk from which these planets form.
Observations show that most circumbinary disks tend to be aligned with their host binary orbital planes. The gas rich circumbinary disks HD 98800 B (Andrews et al. 2010), AK Sco (Czekala et al. 2015), DQ Tau (Czekala et Email: jjz54@cornell.edu al. 2016), and the debris circumbiniary disks α CrB and β Tri (Kennedy et al. 2012b) all have mutual disk-binary inclinations not exceeding 3 • . However, there are some notable exceptions. The circumbinary disk around KH 15D is mildly misaligned with the binary orbital plane by ∼ 10 • -20 • (Winn et al. 2004;Chiang & Murray-Clay 2004;Capelo et al. 2012). Shadows (Marino et al. 2015) and gas kinematics  of the disks in HD 142527 are consistent with a misalignment of ∼ 70 • between the outer circumbinary disk and binary orbital plane (Lacour et al. 2016). The disks (circumbinary and two circumstellar) in the binary protostar IRS 43 are misaligned with each other and with the binary (Brinch et al. 2016). Most intriguingly, the debris disk around the eccentric (e b = 0.77) binary 99 Herculis may be highly inclined: By modeling the resolved images from Hershel, Kennedy et al. (2012a) strongly favor a disk orientation where the disk angular momentum vector is inclined to the binary orbital angular momentum vector by 90 • (polar alignment). Kennedy et al. (2012a) also produced a model with a disk-binary inclination of 30 • which fits the observations, but this configuration is unlikely, since differential precession of dust due to the gravitational influence of the binary would rapidly destroy the disk.
Since star/binary formation takes place in turbulent arXiv:1706.07823v2 [astro-ph.EP] 10 Oct 2017 molecular clouds (McKee & Ostriker 2007), the gas that falls onto the central protostellar core/binary and assembles onto the disk at different times may rotate in different directions (e.g. Bate, Bonnell, & Bromm 2003, see also Bate, Lodato, & Pringle 2010;Fielding et al. 2015). In this scenario, it is reasonable to expect a newly formed binary to be surrounded by a highly misaligned circumbinary disk which forms as a result of continued gas accretion (Foucart & Lai 2013). The observed orientations of circumbinary disks then depend on the long-term inclination evolution driven by binary-disk interactions. Foucart & Lai (2013, 2014 studied the warping and the dissipative torque driving the inclination evolution of a circumbinary disk, assuming a circular binary. Foucart & Lai (2013) considered an infinite disk and included the effect of accretion onto the binary, while Foucart & Lai (2014) considered a more realistic disk of finite size and angular momentum, which can precess coherently around the binary. It was shown that under typical protoplanetary conditions, both viscous torque associated with disk warping and accretion torque tend to damp the mutual disk-binary inclination on timescale much shorter than the disk lifetime (a few Myr). By contrast, a circumstellar disc inside a binary can maintain large misalignment with respect to the binary orbital plane over its entire lifetime (Lubow & Ogilvie 2000;Foucart & Lai 2014). This is consistent with the observations that most circumbinary disks are nearly coplanar with their host binaries. On the other hand, the observed circumbinary disk misalignment (such as in KH 15D and IRS 43) can provide useful constraints on the uncertain aspects of the disc warp theory, such as non-linear effects (Ogilvie 2006) and parametric instabilities due to disk warping (Gammie, Goodman, & Ogilvie 2000;Ogilvie & Latter 2013) However, several recent numerical studies using smoothed particle hydrodynamics (SPH) suggest that other outcomes may be possible for disks around eccentric binaries. Aly et al. (2015) showed that disks around binary black holes (which typically lie in the "viscous regime" of disk warps, with the viscosity parameter α larger than the disk aspect ratio H/r; Papaloizou & Pringle 1983;Ogilvie 1999;see Sec. 3) around eccentric binaries may be driven into polar alignment. Martin & Lubow (2017) found numerically that a circumbinary protoplanetary disk (typically in the bending wave regime, with α H/r Papaloizou & Lin 1995;Lubow & Ogilvie 2000), inclined to an eccentric (e b = 0.5) binary by 60 • will evolve to a polar configuration. They suggested that this dynamical outcome arises from the combined influence of the gravitational torque on the disk from the binary and viscous torques from disk warping. They also proposed that 99 Herculis (with e b = 0.77) followed such an evolution to end in the orbital configuration (polar alignment) observed today.
In this paper, we provide a theoretical anaysis to the above numerical results. In particular, we generalize the study of Foucart & Lai (2014) to apply to circumbinary disks with arbitrary disk-binary inclinations and binary eccentricities. We derive the critical condition and calculate the timescale for the disk to evolve toward polar alignment with the binary. In Section 2, we review the secular dynamics of a test particle around an eccentric binary. In Section 3, we calculate the disk warp profile and dissipative disk torques acting on the disk, and derive the requirements for the disk to evolve into polar alignment with the binary. Section 4 considers the situation when the circumbinary disk has a non-negligible angular momentum compared to the inner binary. In Section 5, we examine the back-reaction torque from the disk on the binary and the effect of gas accretion. We discuss our results in Section 6, and summarize in Section 7.

TEST PARTICLE DYNAMICS
In preparation for later sections, we review the secular dynamics of a test particle surrounding an eccentric binary (Farago & Laskar 2010;Li, Zhou, & Zhang 2014;Naoz et al. 2017). Consider a circular test mass with semimajor axis r and orbital angular momentum unit vectorl, surrounding an eccentric binary with orbital angular momentum vectorl b , eccentricity vector e b , semimajor axis a b , total mass M b = M1 + M2 (where M1, M2 are individual masses), and reduced mass µ b = M1M2/M b . The orbit-averaged torque per unit mass on the test particle is (e.g. Liu, Muñoz, & Lai 2015;Petrovich 2015) where n GM b /r 3 is the test particle orbital frequency (mean-motion), and characterizes the precession frequency of the test particle around the binary. The torque T b in Eq. (1) is evaluated to the lowest order in a b /r. The time evolution of the test particle's orbital angular momentum vector is given by Equation (3) can be solved analytically (Landau & Lifshitz 1969;Farago & Laskar 2010;Li, Zhou, & Zhang 2014), but the dynamics may be easily understood by analyzing the energy curves. Equation (3) has an integral of motion To plot the energy curves, we set up the Cartesian coordinate system (x, y, z), wherel b =ẑ and e b = e bx . We may writel = (sin I sin Ω, − sin I cos Ω, cos I) = (cos Ie, sin Ie sin Ωe, sin Ie cos Ωe), where I is the angle betweenl andl b , Ω is the test particle's longitude of the ascending node (measured in the xy plane from the x-axis); similarly Ie is the angle betweenl and e b , and Ωe measures the longitude of the node in the yz plane (see Fig. 1). In terms of I and Ω, we have   In Figure 2, we plot the test particle trajectories in the I − Ω (left panel) and Ie − Ωe (right panel) planes for the binary eccentricity e b = 0.3. The critical separatrix Λ = 0 is displayed in black in both plots. When Λ > 0,l precesses aroundl b with I ∼ constant and Ω circulating the full range (0 − 360 • ), while Ωe librates around 0 • . When Λ < 0,l precesses around e b with Ie ∼ constant and Ωe circulating the full range (0 − 360 • ), while Ω librates around 90 • (see Fig. 1). Thus, the test particle angular momentum axisl transitions from precession aroundl b for Λ > 0 to precession around e b for Λ < 0. Because the Λ = 0 separatrix has Ω ∈ [0 • , 360 • ] (Fig. 2), a necessary condition forl to precess Figure 2 clearly reveals the stable fixed points of the system. In terms of the variables (Ω, I), the stable fixed points (where dI/dt = dΩ/dt = 0) are I = π/2 and Ω = π/2, 3π/2, corresponding tol = ±e b /e b . In terms of the variables (Ωe, Ie), the fixed points are Ie = π/2 and Ωe = 0, π, corresponding tol = ±l b . We will see in Section 3 that in the presence of dissipation, the disk may be driven toward one of these fixed points.

CIRCUMBINARY DISK DYNAMICS
We now consider a binary (with the same parameters as in Section 2, see also Fig. 3) surrounded by a circular circumbinary disk with inner truncation radius rin, outer truncation radius rout, with unit angular momentum vectorl =l(r, t), and surface density Σ = Σ(r). For concreteness, we adopt the surface density profile Σ(r) = Σin rin r .
We assume rin rout throughout this work. We could assume a more general surface density profile Σ ∝ r −p , with p observationally constrained to lie in the range 0.5 − 1. The binary has orbital angular momentum while the disk has mass Σrdr 2πΣinrinrout (11) and angular momentum (assuming a small disk warp; see below) where n(r) L d for typical circumbinary disk parameters, in this section we assumel b and e b are fixed in time, neglecting the back-reaction torque on the binary from the disk. We discuss the system's dynamics when L d is nonnegligible compared to L b in Section 4, and the effects of the back-reaction torque on the binary from the disk in Section 5.

Qualitative Discussion
Assuming the disk to be nearly flat, the time evolution of the disk unit angular momentum vector is given by where T b is given in Eq.
(1),l d (t) is a suitably averaged unit angular momentum of the disk [see Eq. (24)], and . . . implies a proper average over r [see Eq. (27)]. When the disk is flat, the time evolution ofl d is identical to that of a test particle (see discussion at the end of Sec. 3.2). When α H/r (H is the disk scaleheight, α is the viscosity parameter), the main internal torque enforcing disk rigidity and coherent precession comes from bending wave propigation (Papaloizou & Lin 1995;Lubow & Ogilvie 2000). As bending waves travel at 1/2 the sound speed, the wave crossing time is of order t bw = 2r/cs. When t bw is longer than the characteristic precession time (2)], strong disk warps can be induced. In the extreme nonlinear regime, disk breaking may be possible in circumbinary disks (Larwood & Papaloizou 1997;Facchini, Lodato, & Price 2013;Nixon, King, & Price 2013). To compare t bw with ω −1 b , we adopt the disk sound speed profile where h = H/r. We find Thus, we expect that the small warp approximation should be valid everywhere in the disk except the inner-most region. Throughout this paper, we scale our results to h = 0. Although the disk is flat to a good approximation, the interplay between disk twist/warp and viscous dissipation may modify the disk's dynamics over timescales much longer than ω −1 b . When the external torque T b is applied to the disk in the bending wave regime, the disk's viscosity causes the disk to develop a small twist of order while the precession of bending waves from a non-Keplarian epicyclic frequency κ causes the disk to develop a small warp, of order The viscous twist [Eq. (17)] interacts with the external torque, effecting the evolution ofl over the viscous timescale.
To an order of magnitude, we have In the above estimate, we have assumed the relevant misalignment angles (betweenl d andl b , or betweenl d and e b ) is of order unity.

Formalism
The torque per unit mass on the disk from the inner binary is given by Eq.
(1), with T b = T b (r, t). In addition, the gravitational potential from the binary induces a non-Keplarian angular frequency (Miranda & Lai 2015), with where When the Shakura-Sunaev α-viscosity parameter satisfies α H/r, the disk lies in the bending wave regime (Papaloizou & Lin 1995;Lubow & Ogilvie 2000). Any warp induced by an external torque is smoothed by bending waves passing through the disk. Protoplanetary disks typically lie in the bending wave regime. The time evolution ofl(r, t) is governed by the equations (Lubow & Ogilvie 2000; see also Lubow, Ogilvie, & Pringle 2002) Σr 2 n ∂l ∂t where G is the internal torque. From equation (16), we see that t bw < ω −1 b for standard circumbinary disk parameters, so the disk should be only mildly warped. We may therefore expand l(r, t) =l d (t) + l1(r, t) + . . . , wherel d is the unit vector along the total angular momentum of the disk, |l1| |l d | = 1 [see Eqs. (31)-(32) below]. As we will see, the internal torque G0(r, t) maintains the  (38)] as a function of radius. We take r in = 2 AU and rout = 100 AU. rigid body dynamical evolution ofl d , while G1(r, t) maintains the warp profile l1. Perturbative expansions to study warped disk structure and time evolution have been taken by Lubow & Ogilvie (2000, 2001 and Foucart & Lai (2014). Inserting (24) into Eq. (22), integrating over rdr, and using the zero torque boundary condition G0(rin, t) = G0(rout, t) = 0, we find the leading order time evolution ofl is given by Here, is the characteristic precession frequency of the rigid disk. Equation (27) is equivalent to Equation (3) if one replaces ω b with ω b , and the disk dynamics reduce to those of a test particle withl =l d when cs → ∞.

Disk Warp Profile
Withl d determined with boundary condition (26), we may solve for G0(r, t): With the leading order terms forl and G, we may solve for l1. We impose the normalization condition r out r in so thatl d is the unit vector along the total angular momentum of the disk, or Inserting Eq. (29) into Eq. (23) and integrating, we obtain l1(r, t) = (l1)twist + (l1)warp where and Here, and The third term in Eq. (33) arises from the fact thatl d ·e b and l d ·l b are not constant in time, and is dynamically unimportant. Although it is strait-forward to compute the integrals in Eqs. (36)-(38), this calculation is tedious and unilluminating. Instead, we notice that over most of the region in the integrals, the internal torque radial function g b (r) is of order Evaluating the warp functions and using the fact that rin rout, we obtain the approximate expressions In Figure 4, we plot the rescaled warp
We choose to normalize γ b by α = 0.01; real protoplanetary disks may have α in the range 10 −1 − 10 −5 (Rafikov 2017). Since the viscous dissipation from disk twisting effects the evolution ofl d according to where Λ is given by Eq. (4), except we replacel byl d : Equation (50) is the main result of our technical calculation. We see d dt (52)-(53) show the system has two different end-states depending on the initial value for Λ: The final state ofl d is alignment (or anti-alignment) with e b . Figure 5 shows several examples of the results for the evolution of disk orientation, obtained by integrating the time evolution ofl d , including gravitational [Eq. (27)] and viscous [Eq. (50)] torques. On the left panels, we plot the disk inclination I with time, for the binary eccentricities indicated. We choose the initial Ω(0) = 90 • for all cases, so that I < Icrit (I > Icrit) corresponds exactly to Λ > 0 (Λ < 0) (see Eqs. (4) and (8)). Thus we expect I → 0 • when I < Icrit, I → 90 • when Icrit < I < 180 • − Icrit, and I → 180 • when I > 180 • − Icrit. On the right panels of Figure 5, we plot the disk trajectories on the I − Ω plane [Eq. (6) withl →l d ]. Again, we see when I < Icrit (Λ > 0), l d aligns withl b , while when I > Icrit (Λ < 0),l d aligns with e b , as expected.

SECULAR DYNAMICS WITH MASSIVE INCLINED OUTER BODY
Sections 2-3 neglected the circumbinary disk's angular momentum, a valid assumption as long as L d L b [Eq. (13)]. When L d L b , the non-zero disk angular momentum will change the locations of the fixed points of the system, and hence may effect its dynamical evolution over viscous timescales.
Consider the setup of Section 2, except we now include the outer body's mass m and angular momentum L = m √ GM b rl. The evolution equations forl, j b = 1 − e 2 bl b , and e b are [Liu, Muñoz, & Lai 2015;Eqs. (17) where For a given K, one may solve Eq. (60) to get e 2 b = e 2 b (I). Assuming 0 ≤ I ≤ π/2 and requiring 0 ≤ e b < 1, we obtain Equation (59) then gives Ψ = Ψ(I, Ω). When J ∼ K −1 1, Eq. (61) reduces to while when J 1, Eq. (61) becomes The fixed points of the system in the I − Ω plane are determined by The condition ∂Ψ/∂Ω = 0 gives Ω = π/2 and Ω = 3π/2, as before (see Sec. 2). For arbitrary J, one must numerically solve ∂Ψ/∂I| Ω=π/2,3π/2 = 0 to calculate the fixed points I = I fp > 0 (I = 0 is always a fixed point of the system). However, when J 1, one may show analytically that (as found in Sec. 2) while when J 1, where e 2 b = e 2 b (0). Notice I fp is the Lidov-Kozai critical inclination when J 1 and e b (0) = 0 (Lidov 1962;Kozai 1962). Figure 6 plots trajectories of the system in the I − Ω and e b − Ω planes. When J 1, the system's dynamics reduce to that discussed in Section 2, with I fp 90 • (black x's), e b e b (0), and trajectories above and below I = 90 • are symmetric. As J increases in magnitude, I fp decreases, e b begins to oscillate, and the inclination symmetry above and below I = 90 • is lost. Although different trajectories may cross in the I − Ω plane, each still has a unique Ψ value [Eq. (59)], since the binary's e b value differs from Eq. (61) when I > π/2. When J 1, the system's dynamics approaches the classic Lidov-Kozai regime (Lidov 1962;Kozai 1962). The fixed point I fp of the system approaches Eq. (66), with e b reaching large values when I(0) > I fp , and with trajectories symmetric above and below I = 90 • . Figure 7 plots I fp as a function of J −1 , computed for Ω = π/2 with the e b = e b (0) values as indicated. The two limiting cases given by Eqs. (65) and (66) are achieved when J 1 and J 1, respectively, and I fp generally varies nonmonitonically with increasing J. Since e b should evolve in time under the influence of viscous dissipation from disk warping, one cannot determine the final value of I fp the system may evolve into starting from initial I(0) and e b (0) values without a detailed calculation similar to Section 3. Nevertheless, Figures 6 and 7 show there exist highly inclined fixed points for any value of J. For J 0.1, the system may evolve into near polar alignment, with I fp somewhat less than 90 • .

TORQUE ON BINARY AND EFFECT OF ACCRETION
In the previous sections, we have studied the evolution of the disk around a binary with fixedl b and e b . Here we study In addition, accretion onto the binary from the disk adds angular momentum to the binary's orbit: Here,Ṁ is the mass accretion rate onto the binary, λ ∼ 1 (e.g. Miranda, Muñoz, & Lai 2017), and we have assumed l(rin, t) l d (t) (see below). The torques (67) and (68) are equivalent to those considered in Foucart & Lai (2013), except we give different power-law prescriptions for Σ = Σ(r) and H = H(r). For disks in steady state, we havė M 3παh 2 Σinr 2 in n(rin), Using Eqs. (46) (with e b = 0), (67) and (68), we obtain the net disk-binary alignment timescale for small angle between l d andl b : where measures the strength of the accretion torque to the viscous torque on the binary (η/λ = f −1 , λ = g in the notation of Foucart & Lai 2013). Sincel(rin, t) =l(rout, t), the disk angular momentum loss through accretion causesl d to change with time: Because the magnitude of the tilt ofl(rin, t) froml d is of order we find Detailed calculation shows that the accretion torque (74) is always much less than the viscous torque (19) on the disk. We relegate the calculation and discussion of the accretion torque (74) to the Appendix. For eccentric binaries, the back-reaction toque from the disk is dL b /dt = −dL d /dt [Eq. (46)]. But this is not sufficient for determining the evolution of e b andl b . In addition, how accretion affects the binary eccentricity is also uncertain (e.g. Rafikov 2016;Miranda, Muñoz, & Lai 2017). Nevertheless, as long as L b L d , the timescale for the disk-binary inclination evolution should be of order γ −1 b , with an estimate given by Eq. (48).

Theoretical Uncertainties
Our theoretical analysis of disks around binaries assumes a linear disk warp. However, we find that at the inner disk region, |∂l/∂ ln r| reaches ∼ 0.1 for a wide range of binary and disk parameters. Inclusion of weakly non-linear warps in Equations (22)-(23) may introduce new features in the disk warp profile (Ogilvie 2006). In addition, disk warps of this magnitude may interact resonantly with inertial waves in the disk, leading to a parametric instability which may excite turbulence in the disk (Gammie, Goodman, & Ogilvie 2000;Ogilvie & Latter 2013). An investigation of these effects is outside the scope of this paper, but their inclusion is unlikely to change the direction of disk-binary inclination evolution (alignment vs polar alignment).

Observational Implications
In Section 3.4, we showed that the viscous torque associated with disk twist/warp tends to drive the circumbinary disk axisl d toward ±l b (alignment or anti-alignment) when Λ > 0, and toward ±e b (polar alignment) when Λ < 0. Note that Icrit < I < 180 • − Icrit is a necissary, but not sufficient condition for polar alignment of the disk [Eq. (8)]. An extreme example is when Ω = 0 • , since Λ ≥ 0 for all inclinations I. Because the circumbinary disk probably formed in a turbulent molecular cloud, the disk is unlikely to have a preferred Ω when it forms. The condition for polar alignment (Λ < 0) requires Ω to satisfy Assuming a uniform distribution of Ω-values from 0 to 2π, the probability of the disk to polar align is (for given I, e b ) P polar (I, e b ) = 1 − 2Ωmin/π. (76) where Ωmin(I, e b ) = We define the inclination I polar through P polar (I polar , e b ) = 0.5. Solving for I polar , we obtain In Figure 8, we plot contours of constant P polar in the I − e b space. The P polar = 0 curve (black) traces out Icrit [Eq. (8)], while the P polar = 0.5 curve (red) traces out I polar [Eq. (78)]. When I < Icrit, alignment ofl withl b is inevitable. When I > I polar , alignment ofl with e b is probable. Table 1 lists a number of circumbinary systems with highly eccentric binaries. With the exception of 99 Herculis, all the binaries listed have disks coplanar with the binary orbital plane within a few degrees. We also list Icrit [Eq. (8)] and I polar [Eq. (78)] for these systems. We do not list the binaries KH 15D (Winn et al. 2004;Chiang & Murray-Clay 2004;Capelo et al. 2012) and HD 142527B (Marino et al. 2015;Casassus et al. 2015;Lacour et al. 2016) since the orbital elements of these binaries are not well constrained. However, both binaries appear to have significant eccentricities (Chiang & Murray-Clay 2004;Lacour et al. 2016).
Since planets form in gaseous circumbinary disks, planets may form with orbital planes perpendicular to the binary orbital plane if the binary is sufficiently eccentric. Such planets may be detectable in transit surveys of eclipsing binaries due to nodal precession of the planet's orbit.
The twist and warp calculated in Section 3.3 is nonnegligible. Further observations of (gaseous) circumbinary disks may be able to detect such warps (Juhász & Facchini 2017), further constraining the orientation and dynamics of circumbinary disk systems.

SUMMARY
Using semi-analytic theory, we have studied the warp and long-term evolution of circumbinary disks around eccentric binaries. Our main results and conclusions are listed below.
(i) For protoplanetary disks with dimensionless thickness H/r larger than the viscosity parameter α, bending wave propagation effectively couples different regions of the disk, making it precess as a quasi-rigid body. Without viscous dissipation from disk warping, the dynamics of such a disk is similar to that of a test particle around an eccentric binary (Secs. 2 and 3.2).
(ii) When the binary is eccentric and the disk is significantly inclined, the disk warp profile exhibits new features not seen in previous works. The disk twist [Eq. (34)] and warp [Eq. (35)] have additional contributions due to additional torques on the disk when the binary is eccentric.
(iii) Including the dissipative torque from warping, the disk may evolve to one of two states, depending on the initial sign of Λ [Eq. (4)] (Sec. 3.4). When Λ is initially positive, the disk angular momentum vector aligns (or anti-aligns) with the binary orbital angular momentum vector. When Λ is initially negative, the disk angular momentum vector aligns with the binary eccentricity vector (polar alignment). Note that Λ depends on both I (the disk-binary inclination) and Ω (the longitude of ascending node of the disk). Thus for a given e b , the direction of inclination evolution depends not only on the initial I(0), but also on the initial Ω(0).
(iv) When the disk has a non-negligible angular momentum compared to the binary, the systems fixed points are modified (Sec. 4). The disk may then evolve to a state of near polar alignment, with the inclination somewhat less than 90 • .
(v) The timescale of evolution of the disk-binary inclination angle [see Eqs. (52)-(53)] depends on various disk parameters [see Eq. (48)], but is in general less than a few Myrs. This suggests that highly inclined disks and planets may exist around eccentric binaries. on the I − Ω plane, for the binary eccentricities indicated. All disk trajectories evolve toward the Λ ≈ 0 seperatrix.
The relative strength of the viscous to the accretion torques from disk warping is given by the ratio As long as |γ b | |γaW bb (rin)|, the viscous torque dominates, andl d aligns with eitherl b or e b , depending on the sign of Λ (Sec. 3.4). When |γ b | |γaW bb (rin)|, the accretion torques may dominate, andl d may be driven to the seperatrix Λ ≈ 0. |W bb (rin)γa|, the viscous torque dominates the disk's dynamics. As a result, Figure 10 is almost indistinguishable from Figure 5. Only for unrealistically hot protoplanetary disks with h 0.5 may accretion torques significantly effect the disk evolution over viscous timescales.