Warp Propagation Research Articles

On the conditions for warping and breaking protoplanetary discs

ABSTRACT Recent observations demonstrate that misalignments and other out-of-plane structures are common in protoplanetary discs. Many of these have been linked to a central host binary with an orbit that is inclined with respect to the disc. We present simulations of misaligned circumbinary discs with a range of parameters to gain a better understanding of the link between those parameters and the disc morphology in the wave-like regime of warp propagation that is appropriate to protoplanetary discs. The simulations confirm that disc tearing is possible in protoplanetary discs as long as the mass ratio, μ, and disc–binary inclination angle, i, are not too small. For the simulations presented here, this corresponds to μ > 0.1 and i ≳ 40°. For highly eccentric binaries, tearing can occur for discs with smaller misalignment. Existing theoretical predictions provide an estimate of the radial extent of the disc in which we can expect breaking to occur. However, there does not seem to be a simple relationship between the disc properties and the radius within the circumbinary disc at which the breaks appear, and furthermore the radius at which the disc breaks can change as a function of time in each case. We discuss the implications of our results for interpreting observations and suggest some considerations for modelling misaligned discs in the future.

Dynamical dust traps in misaligned circumbinary discs: analytical theory and numerical simulations

ABSTRACT Recent observations have shown that circumbinary discs can be misaligned with respect to the binary orbital plane. The lack of spherical symmetry, together with the non-planar geometry of these systems, causes differential precession which might induce the propagation of warps. While gas dynamics in such environments is well understood, little is known about dusty discs. In this work, we analytically study the problem of dust traps formation in misaligned circumbinary discs. We find that pile-ups may be induced not by pressure maxima, as the usual dust traps, but by a difference in precession rates between the gas and dust. Indeed, this difference makes the radial drift inefficient in two locations, leading to the formation of two dust rings whose position depends on the system parameters. This phenomenon is likely to occur to marginally coupled dust particles (St ≳ 1) as both the effect of gravitational and drag force are considerable. We then perform a suite of 3D Smoothed Particle Hydrodynamics (SPH) numerical simulations to compare the results with our theoretical predictions. We explore the parameter space, varying stellar mass ratio, disc thickness, radial extension, and we find a general agreement with the analytical expectations. Such dust pile-up prevents radial drift, fosters dust growth and may thus promote the planet formation in circumbinary discs.

Kozai–Lidov oscillations triggered by a tilt instability of detached circumplanetary discs

ABSTRACT Circumplanetary discs can be linearly unstable to the growth of disc tilt in the tidal potential of the star–planet system. We use 3D hydrodynamical simulations to characterize the disc conditions needed for instability, together with its long-term evolution. Tilt growth occurs for disc aspect ratios, evaluated near the disc outer edge, of H/r ≳ 0.05, with a weak dependence on viscosity in the wave-like regime of warp propagation. Lower mass giant planets are more likely to have circumplanetary discs that satisfy the conditions for instability. We show that the tilt instability can excite the inclination to above the threshold where the circumplanetary disc becomes unstable to Kozai–Lidov (KL) oscillations. Dissipation in the KL unstable regime caps further tilt growth, but the disc experiences large oscillations in both inclination and eccentricity. Planetary accretion occurs in episodic accretion events. We discuss implications of the joint tilt–KL instability for the detectability of circumplanetary discs, for the obliquity evolution of forming giant planets, and for the formation of satellite systems.

Sound Speed Dependence of Alignment in Accretion Disks Subjected to Lense–Thirring Torques

We present a series of simulations in both pure hydrodynamics (HD) and magnetohydrodynamics (MHD) exploring the degree to which alignment of disks subjected to external precessional torques (e.g., as in the `Bardeen-Petterson' effect) is dependent upon the disk sound speed c_s. Across the range of sound speeds examined, we find that the influence of the sound speed can be encapsulated in a simple "lumped-parameter" model proposed by Sorathia et al. (2013a). In this model, alignment fronts propagate outward at a speed ~0.2 rOmega_precess(r), where Omega_precess is the local test-particle precession frequency. Meanwhile, transonic radial motions transport angular momentum both inward and outward at a rate that may be described roughly in terms of an orientation diffusion model with diffusion coefficient ~2c_s^2/Omega, for local orbital frequency Omega. The competition between the two leads, in isothermal disks, to a stationary position for the alignment front at a radius proportional to c_s^(-4/5). For alignment to happen at all, the disk must either be turbulent due to the magnetorotational instability in MHD, or, in HD, it must be cool enough for the bending waves driven by disk warp to be nonlinear at their launch point. Contrary to long-standing predictions, warp propagation in MHD disks is diffusive independent of the parameter c_s/(alpha v_orb$, for orbital speed v_orb and ratio of stress to pressure of alpha. In purely HD disks, i.e., those with no internal stresses other than bulk viscosity, warmer disks align weakly or not at all; cooler disks align qualitatively similarly to MHD disks.

A local model of warped magnetized accretion discs

We derive expressions for the local ideal magnetohydrodynamic (MHD) equations for a warped astrophysical disc using a warped shearing box formalism. A perturbation expansion of these equations to first order in the warping amplitude leads to a linear theory for the internal local structure of magnetised warped discs in the absence of MRI turbulence. In the special case of an external magnetic field oriented normal to the disc surface these equations are solved semi-analytically via a spectral method. The relatively rapid warp propagation of low-viscosity Keplerian hydrodynamic warped discs is diminished by the presence of a magnetic field. The magnetic tension adds a stiffness to the epicyclic oscillations, detuning the natural frequency from the orbital frequency and thereby removing the resonant forcing of epicyclic modes characteristic of hydrodynamic warped discs. In contrast to a single hydrodynamic resonance we find a series of Alfv\'{e}nic-epicyclic modes which may be resonantly forced by the warped geometry at critical values of the orbital shear rate $q$ and magnetic field strength. At these critical points large internal torques are generated and anomalously rapid warp propagation occurs. As our treatment omits MRI turbulence, these results are of greatest applicability to strongly magnetised discs.

Wave-like warp propagation in circumbinary discs – II. Application to KH 15D

KH 15D is a protostellar binary system that shows a peculiar light curve. In order to model it, a narrow circumbinary precessing disc has been invoked, but a proper dynamical model has never been developed. In this paper, we analytically address the issue of whether such a disc can rigidly precess around KH 15D, and we relate the precessional period to the main parameters of the system. Then, we simulate the disc's dynamics by using a 1D model developed in a companion paper, such that the warp propagates into the disc as a bending wave, which is expected to be the case for protostellar discs. The validity of such an approach has been confirmed by comparing its results with full 3D SPH simulations on extended discs. In the present case, we use this 1D code to model the propagation of the warp in a narrow disc. If the inner truncation radius of the disc is set by the binary tidal torques at {\sim} 1 AU, we find that the disc should extend out to 6-10 AU (depending on the models), and is therefore wider than previously suggested. Our simulations show that such a disc does reach an almost steady state, and then precesses as a rigid body. The disc displays a very small warp, with a tilt inclination that increases with radius in order to keep the disc in equilibrium against the binary torque. However, for such wider discs, the presence of viscosity leads to a secular decay of the tilt on a timescale of {\approx} 3000 ({\alpha}/0.05)^(-1) years, where {\alpha} is the disc viscosity parameter. The presence of a third body (such as a planet), orbiting at roughly 10 AU might simultaneously explain the outer truncation of the disc and the maintenance of the tilt for a prolonged time.

Wave-like warp propagation in circumbinary discs – I. Analytic theory and numerical simulations

In this paper we analyse the propagation of warps in protostellar circumbinary discs. We use these systems as a test environment in which to study warp propagation in the bending-wave regime, with the addition of an external torque due to the binary gravitational potential. In particular, we want to test the linear regime, for which an analytic theory has been developed. In order to do so, we first compute analytically the steady state shape of an inviscid disc subject to the binary torques. The steady state tilt is a monotonically increasing function of radius. In the absence of viscosity, the disc does not present any twist. Then, we compare the time-dependent evolution of the warped disc calculated via the known linearised equations both with the analytic solutions and with full 3D numerical simulations, which have been performed with the PHANTOM SPH code using 2 million particles. We find a good agreement both in the tilt and in the phase evolution for small inclinations, even at very low viscosities. Moreover, we have verified that the linearised equations are able to reproduce the diffusive behaviour when {\alpha} > H/R, where {\alpha} is the disc viscosity parameter. Finally, we have used the 3D simulations to explore the non-linear regime. We observe a strongly non-linear behaviour, which leads to the breaking of the disc. Then, the inner disc starts precessing with its own precessional frequency. This behaviour has already been observed with numerical simulations in accretion discs around spinning black holes. The evolution of circumstellar accretion discs strongly depends on the warp evolution. Therefore the issue explored in this paper could be of fundamental importance in order to understand the evolution of accretion discs in crowded environments, when the gravitational interaction with other stars is highly likely, and in multiple systems.

Black hole mergers: do gas discs lead to spin alignment?

Abstract In this Letter, we revisit arguments suggesting that the Bardeen–Petterson effect can coalign the spins of a central supermassive black hole binary accreting from a circumbinary (or circumnuclear) gas disc. We improve on previous estimates by adding the dependence on system parameters and noting that the non-linear nature of warp propagation in a thin viscous disc affects alignment. This reduces the disc's ability to communicate the warp, and can severely reduce the effectiveness of disc-assisted spin alignment. We test our predictions with a Monte Carlo realization of random misalignments and accretion rates, and we find that the outcome depends strongly on the spin magnitude. We estimate a generous upper limit to the probability of alignment by making assumptions which favour it throughout. Even with these assumptions, about 40 per cent of black holes with a ≳ 0.5 do not have time to align with the disc. If the residual misalignment is not small and it is maintained down to the final coalescence phase, this can give a powerful recoil velocity to the merged hole. Highly spinning black holes are thus more likely being subject to strong recoils, the occurrence of which is currently debated.

Broken discs: warp propagation in accretion discs

We simulate the viscous evolution of an accretion disc around a spinning black hole. In general any such disc is misaligned, and warped by the Lense-Thirring effect. Unlike previous studies we use effective viscosities constrained to be consistent with the internal fluid dynamics of the disc. We find that nonlinear fluid effects, which reduce the effective viscosities in warped regions, can promote the breaking of the disc into two distinct planes. This occurs when the Shakura & Sunyaev dimensionless viscosity parameter alpha is <~ 0.3 and the initial angle of misalignment between the disc and hole is >~ 45 degrees. The break can be a long-lived feature, propagating outwards in the disc on the usual alignment timescale, after which the disc is fully co- or counter-aligned with the hole. Such a break in the disc may be significant in systems where we know the inclination of the outer accretion disc to the line of sight, such as some X-ray binaries: the inner disc, and so any jets, may be noticeably misaligned with respect to the orbital plane.

On the diffusive propagation of warps in thin accretion discs

In this paper we revisit the issue of the propagation of warps in thin and viscous accretion discs. In this regime warps are know to propagate diffusively, with a diffusion coefficient approximately inversely proportional to the disc viscosity. Previous numerical investigations of this problem (Lodato & Pringle 2007) did not find a good agreement between the numerical results and the predictions of the analytic theories of warp propagation, both in the linear and in the non-linear case. Here, we take advantage of a new, low-memory and highly efficient SPH code to run a large set of very high resolution simulations (up to 20 million SPH particles) of warp propagation, implementing an isotropic disc viscosity in different ways, to investigate the origin of the discrepancy between the theory and the numerical results. Our new and improved analysis now shows a remarkable agreement with the analytic theory both in the linear and in the non-linear regime, in terms of warp diffusion coefficient and precession rate. It is worth noting that the resulting diffusion coefficient is inversely proportional to the disc viscosity only for small amplitude warps and small values of the disc $\alpha$ coefficient ($\alpha < 0.1$). For non-linear warps, the diffusion coefficient is a function of both radius and time, and is significantly smaller than the standard value. Warped accretion discs are present in many contexts, from protostellar discs to accretion discs around supermassive black holes. In all such cases, the exact value of the warp diffusion coefficient may strongly affect the evolution of the system and therefore its careful evaluation is critical in order to correctly estimate the system dynamics (abridged).

Warp and eccentricity propagation in discs around black holes

We consider the inward propagation of warping and eccentric disturbances in discs around black holes under a wide variety of conditions. In our calculations, we use secular theories of warped and eccentric discs and assume the deformations to be stationary and propagating in a disc model similar to regions (a) and (b) of Shakura & Sunyaev discs. We find that the propagation of deformations to the innermost regions of the disc is facilitated for low viscous damping and high accretion rate. We relate our results to the possible excitation of trapped inertial modes, and to the observations of high-frequency quasi-periodic oscillations (QPOs) in black hole systems in the very high spectral state.

Hydrodynamic simulations of the Bardeen-Petterson effect

We present SPH simulations of accretion discs in orbit about rotating compact objects, such as black holes and neutron stars, and study the structure of warped discs produced by the Bardeen-Petterson effect. We calculate the transition radius out to which the disc specific angular momentum vector is aligned with that of the black hole. We focus on the parameter regime where the warp dynamics are controlled by bending wave propagation, but also consider models in which warps are subject to diffusion rather than wave transport, and are able to consider the fully non-linear regime. Because of hydrodynamic or pressure effects, for the parameter range investigated, the transition radius is always found to be much smaller than that obtained by Bardeen & Petterson. For discs with mid-plane Mach numbers of ∼10, the transition occurs between 10 and 16R+ (gravitational radii), whereas for a Mach number of ∼30 it occurs at around 30R+. A thicker disc with a Mach number of 5 is found to produce no discernible warped structure. The rate of black hole-disc alignment is found to be consistent with the ideas of Rees, with the alignment torque behaving as if it arises from the accreted material transferring its misaligned component of angular momentum at the larger transition radius of Bardeen & Petterson. The inclusion of Einstein precession in the calculations modified both the warped disc structure and, consistent with linear analysis, produced an increased alignment rate by up to a factor of 4 because of the effect that the non-Keplerian potential has on the propagation of warps.

Non-Linear Excitation of Warps by Spiral Waves in Galaxies

AbstractWe present a mechanism to explain the permanent excitation of warps in spiral or barred galaxies. We show that non-linear mode coupling with the spiral wave can be very efficient to drive two warp waves. We derive an expression of the coupling coefficient in the WKB approximation. We show that the coupling is too weak in the stellar disk, except at the Outer Lindblad Resonance where the spiral slows down and is efficiently coupled to warp waves. There, the spiral can almost completely transfer its energy and angular momentum to a ”transmitted“ warp, which we can observe in HI, and a ”reflected“ one, which we can observe as a corrugation. Our mechanism reproduces the observed amplitudes of the warp and of the corrugation.We first present side results concerning the propagation of warps in moderately thick disks.