(Abridged) We analyse the stability and evolution of power-law accretion disc models. These have midplane densities that follow radial power-laws, and have either temperature or entropy distributions that are power-law functions of cylindrical radius
. We employ two different hydrodynamic codes to perform 2D-axisymmetric and 3D simulations that examine the long-term evolution of the disc models as a function of the power-law indices of the temperature or entropy, the thermal relaxation time of the fluid, and the viscosity. We present a stability analysis of the problem that we use to interpret the simulation results. We find that disc models whose temperature or entropy profiles cause the equilibrium angular velocity to vary with height are unstable to the growth of modes with wavenumber ratios |k_R/k_Z| >> 1 when the thermodynamic response of the fluid is isothermal, or the thermal evolution time is comparable to or shorter than the local dynamical time scale. These discs are subject to the Goldreich-Schubert-Fricke (GSF) or `vertical shear linear instability. Development of the instability involves excitation of vertical breathing and corrugation modes in the disc, with the corrugation modes in particular being a feature of the nonlinear saturated state. Instability operates when the dimensionless disc kinematic viscosity nu < 10^{-6} (Reynolds numbers Re>H c_s/nu > 2500). In 3D the instability generates a quasi-turbulent flow, and the Reynolds stress produces a fluctuating effective viscosity coefficient whose mean value reaches alpha ~ 6 x 10^{-4} by the end of the simulation. The vertical shear instability in disc models which include realistic thermal physics has yet to be examined. Should it occur, however, our results suggest that it will have significant consequences for their internal dynamics, transport properties, and observational appearance.
Context. Abridged. Many stars are members of binary systems. During early phases when the stars are surrounded by discs, the binary orbit and disc midplane may be mutually inclined. The discs around T Tauri stars will become mildly warped and undergo
solid body precession around the angular momentum vector of the binary system. It is unclear how planetesimals in such a disc will evolve and affect planet formation. Aims. We investigate the dynamics of planetesimals embedded in discs that are perturbed by a binary companion on a circular, inclined orbit. We examine collisional velocities of the planetesimals to determine when they can grow through accretion. We vary the binary inclination, binary separation, D, disc mass, and planetesimal radius. Our standard model has D=60 AU, inclination=45 deg, and a disc mass equivalent to the MMSN. Methods. We use a 3D hydrodynamics code to model the disc. Planetesimals are test particles which experience gas drag, the gravitational force of the disc, the companion star gravity. Planetesimal orbit crossing events are detected and used to estimate collisional velocities. Results. For binary systems with modest inclination (25 deg), disc gravity prevents planetesimal orbits from undergoing strong differential nodal precession (which occurs in absence of the disc), and forces planetesimals to precess with the disc on average. For bodies of different size the orbit planes become modestly mutually inclined, leading to collisional velocities that inhibit growth. For larger inclinations (45 degrees), the Kozai effect operates, leading to destructively large relative velocities. Conclusions. Planet formation via planetesimal accretion is difficult in an inclined binary system with parameters similar to those considered in this paper. For systems in which the Kozai mechanism operates, the prospects for forming planets are very remote.
The gravitational interaction between a protoplanetary disc and planetary sized bodies that form within it leads to the exchange of angular momentum, resulting in migration of the planets and possible gap formation in the disc for more massive planet
s. In this article, we review the basic theory of disc-planet interactions, and discuss the results of recent numerical simulations of planets embedded in protoplanetary discs. We consider the migration of low mass planets and recent developments in our understanding of so-called type I migration when a fuller treatment of the disc thermodynamics is included. We discuss the runaway migration of intermediate mass planets (so-called type III migration), and the migration of giant planets (type II migration) and the associated gap formation in the disc. The availability of high performance computing facilities has enabled global simulations of magnetised, turbulent discs to be computed, and we discuss recent results for both low and high mass planets embedded in such discs.
