The excitation, propagation and dissipation of waves in accretion discs: the non-linear axisymmetric case

We analyse the non-linear propagation and dissipation of axisymmetric waves in accretion discs using the ZEUS-2D hydrodynamics code. The waves are numerically resolved in the vertical and radial directions. Both vertically isothermal and thermally stratified accretion discs are considered. The waves are generated by means of resonant forcing, and several forms of forcing are considered. Compressional motions are taken to be locally adiabatic  ( γ =5/3)  . Prior to non-linear dissipation, the numerical results are in excellent agreement with the linear theory of wave channelling in predicting the types of modes that are excited, the energy flux by carried by each mode, and the vertical wave energy distribution as a function of radius. In all cases, waves are excited that propagate on both sides of the resonance (inwards and outwards). For vertically isothermal discs, non-linear dissipation occurs primarily through shocks that result from the classical steepening of acoustic waves. For discs that are substantially thermally stratified, wave channelling is the primary mechanism for shock generation. Wave channelling boosts the Mach number of the wave by vertically confining the wave to a small cool region at the base of the disc atmosphere. In general, outwardly propagating waves with Mach numbers near resonance  ℳ_r≳0.01  undergo shocks within a distance of order the resonance radius.     

Grains accretion processes in a protoplanetary nebula: Conducting and insulating materials

In a previous work the authors estimated cross-sections for constructive and destructive collisions (‘constructive’ and ‘destructive’ cross-sections) related to silica grains embedded in a protoplanetary nebula. The outstanding point of the problem has induced the authors to examine more closely the physical conditions in which the accretion processes occur, and their dependence on the materials, by which the particles are probably made. In this paper the interaction processes among conducting grains (iron, graphite) and among the insulating (silica) ones have been considered. The following results have been obtained: (1) insulating grains have smaller constructive crossections than conducting; and (2) conducting grains show very large cross-sections for masses up to about 10^?9 g.     

Grains accretion processes in a protoplanetary nebula: Effects of the presence of liquid drops

The temporal behaviour of the mass distribution function of iron solid and/or liquid grains is studied in the inner regions of the solar nebula during its gravitational collapse. At distances from the centre of the nebula equal to 0.2 AU, 0.5 AU and 0.7 AU a mass limit for the accretionary process equal, respectively, to 17.6 g, 5.2 g and 3×10^–3 g is obtained in a reasonable time interval. The importance of a drop phase during the growth of the particles and of the sedimentation of the grains during the disk phase is underlined.     

Theory of protoplanetary disks

The evolution of protoplanetary disks around young stars is briefly reviewed. The most important physical mechanisms that drive the mass accretion are gravitational, magnetic, and thermal convective instabilities. These mechanisms are dominant in different regions of the disk and at different evolutionary epochs.     

Magnetic fields in protoplanetary disks

Magnetic fields likely play a key role in the dynamics and evolution of protoplanetary disks. They have the potential to efficiently transport angular momentum by MHD turbulence or via the magnetocentrifugal acceleration of outflows from the disk surface. Magnetically-driven mixing has implications for disk chemistry and evolution of the grain population, and the effective viscous response of the disk determines whether planets migrate inwards or outwards. However, the weak ionisation of protoplanetary disks means that magnetic fields may not be able to effectively couple to the matter. I examine the magnetic diffusivity in a minimum solar nebula model and present calculations of the ionisation equilibrium and magnetic diffusivity as a function of height from the disk midplane at radii of 1 and 5 AU. Dust grains tend to suppress magnetic coupling by soaking up electrons and ions from the gas phase and reducing the conductivity of the gas by many orders of magnitude. However, once grains have grown to a few microns in size their effect starts to wane and magnetic fields can begin to couple to the gas even at the disk midplane. Because ions are generally decoupled from the magnetic field by neutral collisions while electrons are not, the Hall effect tends to dominate the diffusion of the magnetic field when it is able to partially couple to the gas, except at the disk surfaces where the low density of neutrals permits the ions to remain attached to the field lines. For a standard population of 0.1 μm grains the active surface layers have a combined column Σ_active≈2 g cm⁻² at 1 AU; by the time grains have aggregated to 3 μm, Σ_active≈80 g cm⁻². Ionisation in the active layers is dominated by stellar X-rays. In the absence of grains, X-rays maintain magnetic coupling to 10% of the disk material at 1 AU (i.e. Σ_active≈150 g cm⁻²). At 5 AU the Σ_active≈Σ_total once grains have aggregated to 1 μm in size.     

Magnetically-aided evolution of protoplanetary disks

We investigate the global evolution of a turbulent protoplanetary disk incorporating the effects of Maxwell stress due to a large-scale magnetic field permeating the disk. A magnetic field is produced continuously by an dynamo and the resultant Maxwell stress assists the viscous stress in p roviding the means for disk evolution. The most striking feature of magnetized disk evolution is the presence of the surface density bulge located in the magnetic gap, the region of the disk where the degree of ionization is too low to allow for coupli ng between the magnetic field and the gas. The bulge persists for a time of the order of 10⁵–10⁶ yr. The presence and persistence of the surface density bulge may have important implications for the process of planet formation and the overall characteristics of resultant planetary systems.Operated by USRA under contract No. NASW-4574 with NASA.     

SPH simulations of structures in protoplanetary disks

Using the GADGET-2 code modified by us, we have computed hydrodynamic models of a protoplanetary disk perturbed by a low-mass companion. We have considered the cases of circular and eccentric orbits coplanar with the disk and inclined relative to its midplane. During our simulations we computed the column density of test particles on the line of sight between the central star and observer. On this basis we computed the column density of circumstellar dust by assuming the dust and gas to be well mixed with a mass ratio of 1: 100. To study the influence of the disk orientation relative to the observer on the interstellar extinction, we performed our computations for four inclinations of the line of sight to the disk plane and eight azimuthal directions. The column densities in the circumstellar disk of the central star and the circumbinary disk were computed separately. Our computations have shown that periodic column density oscillations can arise in both inner and circumbinary disks. The amplitude and shape of these oscillations depend on the system’s parameters (the orbital eccentricity and inclination, the component mass ratio) and its orientation in space. The results of our simulations can be used to explain the cyclic brightness variations of young UX Ori stars.     

Secular ring instability in the protoplanetary accretion disk

The basic parameters describing the angular momentum distribution within the Uranus system and of its tidal evolution have been estimated. The nine satellites orbiting under the synchronous zone of Uranus is the maximum number in the solar system and it makes the Uranus system different compared with any other in the Solar system, however the satellites in question are relatively small and their contribution of the tidal dynamics of the system is small compared with that due to UI and UV. The time for existence of the nine satellites as integrated bodies can be estimated as 1.4 × 10⁹ y (UVI) and more. The total tidal decrease in the Uranus angular velocity of rotation is estimated as 7 × 10^–9s^–1.     

Wave propagation in protoplanetary disks: Formation of twin planets by ‘disk-brown dwarf’ collisions?

In this paper we derive a special linear non-vortical wave propagation solution in the shearing sheet, a model of a compressible two-dimensional fluid system with constant density, constant shear and constant Coriolis force, but without self-gravity. The linear analysis of the shearing sheet leads to a single differential equation for the azimuthal velocity perturbation. A detailed derivation of a special solution with a prescribed azimuthal wavenumber k is presented. More general wave solutions, eventually excited by large local ‘impacts’, can be derived by superimposing all k-modes. The special wave functions so obtained describe the formation of two independent spiral wave arms originating out of a ring-shaped structure. The motivation for this investigation lies in the fact that similar wave propagators can be excited by the transit of a solid or ‘clumpy’ object through a protoplanetary disk. We speculate that a disk-brown dwarf collision can produce in the disk a pair of two spiral density wave fragments triggering the rapid accretion of two giant planets by a gravitational shear instability simultaneously (Hypothesis of a mechanism for the production of giant planets in pairs).     

Episodic accretion in magnetically layered protoplanetary discs

We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 au such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion on to the protostar occurs in ∼10⁴-yr bursts with ̇ ∼10⁻⁵ M_⊙ yr⁻¹, separated by quiescent intervals lasting ∼10⁵ yr where ̇ ≈10⁻⁸ M_⊙ yr⁻¹. Such bursts could lead to repeated episodes of strong mass outflow in young stellar objects. The transition to this episodic mode of accretion occurs at an early epoch ( t ≪1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum-mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at r ∼1 au. The predicted rate of low-mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1–3 au.     

Two-dimensional transport of solids in viscous protoplanetary disks

Large-scale radial transport of solids appears to be a fundamental consequence of protoplanetary disk evolution based on the presence of high temperature minerals in comets and the outer regions of protoplanetary disks around other stars. Further, inward transport of solids from the outer regions of the solar nebula has been postulated to be the manner in which short-lived radionuclides were introduced to the terrestrial planet region and the cause of the variations in oxygen isotope ratios in the primitive materials. Here, both outward and inward transport of solids are investigated in the context of a two-dimensional, viscously evolving protoplanetary disk. The dynamics of solids are investigated to determine how they depend on particle size and the particular stage of protoplanetary disk evolution, corresponding to different rates of mass transport. It is found that the outward flows that arise around the disk midplane of a protoplanetary disk aid in the outward transport of solids up to the size of CAIs s and can increase the crystallinity fraction of silicate dust at 10 AU around a solar mass star to as much as ∼40% in the case of rapidly evolving disks, decreasing as the accretion rate onto the star slows. High velocity, inward flows along the disk surface aid in the rapid transport of solids from the outer disk to the inner disk, particularly for small dust. Despite the diffusion that occurs throughout the disk, the large-scale, meridonal flows associated with mass transport prevent complete homogenization of the disk, allowing compositional gradients to develop that vary in intensity for a timescale of one million of years. The variations in the rates and the preferred direction of radial transport with height above the disk midplane thus have important implications for the dynamics and chemical evolution of primitive materials.     

On oligarchic growth of planets in protoplanetary disks

In this paper we present a new semianalytical model of oligarchic growth of planets considering a distribution of planetesimal sizes, fragmentation of planetesimals in mutual collisions, sublimation of ices through the snow line, random velocities out of equilibrium and merging of planetary embryos. We show that the presence of several planetary embryos growing simultaneously at different locations in the protoplanetary disk affects the whole accretion history, specially for the innermost planets. The results presented here clearly indicate the relevance of considering a distribution of planetesimal sizes. Fragmentation occurring during planetesimal-planetesimal collisions represent only a marginal effect in shaping the surface density of solid material in the protoplanetary disc.     

Resonant character of the hall instability in protoplanetary disks

We consider nonaxisymmetric magnetosonic oscillations of a radially stratified, weakly ionized protoplanetary disk with a vertical magnetic field. The combined effect of the Hall electric field and the density and magnetic field inhomogeneities present in the disk has been previously predicted to lead to an instability of its small azimuthal perturbations. We revise the previous results and take into account the effect of inhomogeneous ionization of the protoplanetary material related to the inhomogeneity of the disk medium. We show that the instability criterion is governed by three parameters: the magnetic field and ionization fraction gradients and the plasma β. We have found that at high values of β typical of protoplanetary disks, the instability does not manifest itself if the gradients are directed oppositely. In the case of codirectional gradients, the interaction of magnetosonic fluctuations with inhomogeneities of a fixed size is resonant in character, giving rise to an instability in a narrow range of wave numbers.     

Signatures of planets in protoplanetary and debris disks

We discuss selected possibilities to detect planets in circumstellar disks. We consider the search for characteristic signatures in these disks caused by the interaction of giant planets with the disk as the most promising approach. Numerical simulations show that these signatures are usually much larger in size than the planet itself and thus much easier to detect. The particular result of the planet–disk interaction depends on the evolutionary stage of the disk. Primary signatures of planets embedded in disks are gaps in the case of young disks and characteristic asymmetric density patterns in debris disks.We present simulations which demonstrate that high spatial resolution observations performed with instruments/telescopes that will become available in the near future will be able to trace the location and other properties of young and evolved planets. These observations will allow to directly investigate the formation and evolution of planets in protoplanetary and debris disks.     

A submillimeter view of protoplanetary dust disks

We present some results from our submillimeter single-dish and aperture synthesis imaging surveys of protoplanetary disks using the JCMT, CSO, and Submillimeter Array (SMA) on Mauna Kea, Hawaii. Employing a simple disk model, we simultaneously fit the spectral energy distributions and spatially resolved submillimeter continuum emission from our SMA survey to constrain disk structure properties, including surface density profiles and sizes. The typical disk structure we infer is consistent with a fiducial accretion disk model with a viscosity parameter α≈0.01. Combined with a large, multiwavelength single-dish survey of similar disks, we show how these observations provide evidence for significant grain growth and rapid evolution in the outer regions of disks, perhaps due to an internal photoevaporation process. In addition, we discuss SMA observations of the disks in the Orion Trapezium (proplyds) in the context of disk evolution in a more extreme environment.     

Chaotic behavior in accretion disks

The eccentric luminosity variations of quasars are still a mystery. Analytic results of this behavior ranged from multi-periodic behavior to a purely random process. Recently, we have used the non-linear time-series analysis techniques to analyze the optical light curve of 3C 273 and found that its eccentric behavior may be a low-dimensional chaos. This result induces us to seek some non-linear mechanism for the eccentric luminosity variation. In this paper, we propose a simple non-linear accretion disk model and find that under some circumstances, the chaos appears both in the disk and in the light curve. Then we compute the outburst energy ΔI, defined by the difference between the maximum luminosity and the minimum luminosity, and the mean luminosity 〈I〉 and find that when chaos appears, ΔI～〈I〉 ^α ～M ^0.5α , where M is the mass of black hole and α≈1. These results are confirmed by and/or compatible with the observational data analysis.     

Dissipation of thick accretion disks

A thick accretion disk which is isentropic cannot have simple laminar flow because fluid elements follow orbits which intersect the orbits of other fluid elements, leading to turbulence in astrophysical disks which have very large Reynolds numbers. The turbulence in such disks is estimated using molecular analogies for the behavior of the fluid elements. The usual empirical dissipation parameter ‘α’ is found to be equal to 0.25 under normal circumstances. Characteristic local disk parameters are calculated for a variety of conditions at different distances from a central star of one solar mass. Circumstances involving low midplane optical depths or external heating which can lead to large reductions in the turbulence are discussed.