Dust settling in local simulations of turbulent protoplanetary disks

In this paper, we study the effect of MHD turbulence on the dynamics of dust particles in protoplanetary disks. We vary the size of the particles and relate the dust evolution to the turbulent velocity fluctuations. We performed numerical simulations using two Eulerian MHD codes, both based on finite difference techniques: ZEUS--3D and NIRVANA. These were local shearing box simulations incorporating vertical stratification. Both ideal and non ideal MHD simulations with midplane dead zones were carried out. The codes were extended to incorporate different models for the dust as an additional fluid component. Good agreement between results obtained using the different approaches was obtained. The simulations show that a thin layer of very small dust particles is diffusively spread over the full vertical extent of the disk. We show that a simple description obtained using the diffusion equation with a diffusion coefficient simply expressed in terms of the velocity correlations accurately matches the results. Dust settling starts to become apparent for particle sizes of the order of 1 to 10 centimeters for which the gas begins to decouple in a standard solar nebula model at 5.2 AU. However, for particles which are 10 centimeters in size, complete settling toward a very thin midplane layer is prevented by turbulent motions within the disk, even in the presence of a midplane dead zone of significant size. These results indicate that, when present, MHD turbulence affects dust dynamics in protoplanetary disks. We find that the evolution and settling of the dust can be accurately modelled using an advection diffusion equation that incorporates vertical settling. The value of the diffusion coefficient can be calculated from the turbulent velocity field when that is known for a time of several local orbits.Comment: 15 pages, 16 figures, accepted in Astronomy & Astrophysic

How to adapt broad-band gravitational-wave searches for r-modes

Up to now there has been no search for gravitational waves from the r-modes of neutron stars in spite of the theoretical interest in the subject. Several oddities of r-modes must be addressed to obtain an observational result: The gravitational radiation field is dominated by the mass current (gravitomagnetic) quadrupole rather than the usual mass quadrupole, and the consequent difference in polarization affects detection statistics and parameter estimation. To astrophysically interpret a detection or upper limit it is necessary to convert the wave amplitude to an r-mode amplitude. Also, it is helpful to know indirect limits on gravitational-wave emission to gauge the interest of various searches. Here I address these issues, thereby providing the ingredients to adapt broad-band searches for continuous gravitational waves to obtain r-mode results. I also show that searches of existing data can already have interesting sensitivities to r-modes.Comment: 8 pages, no figure

Corotation Resonance and Diskoseismology Modes of Black Hole Accretion Disks

We demonstrate that the corotation resonance affects only some non-axisymmetric g-mode oscillations of thin accretion disks, since it is located within their capture zones. Using a more general (weaker radial WKB approximation) formulation of the governing equations, such g-modes, treated as perfect fluid perturbations, are shown to formally diverge at the position of the corotation resonance. A small amount of viscosity adds a small imaginary part to the eigenfrequency which has been shown to induce a secular instability (mode growth) if it acts hydrodynamically. The g-mode corotation resonance divergence disappears, but the mode magnitude can remain largest at the place of the corotation resonance. For the known g-modes with moderate values of the radial mode number and axial mode number (and any vertical mode number), the corotation resonance lies well outside their trapping region (and inside the innermost stable circular orbit), so the observationally relevant modes are unaffected by the resonance. The axisymmetric g-mode has been seen by Reynolds & Miller in a recent inviscid hydrodynamic accretion disk global numerical simulation. We also point out that the g-mode eigenfrequencies are approximately proportional to m for axial mode numbers |m|>0.Comment: 16 pages, no figures. Submitted to The Astrophysical Journa

On the orbital evolution and growth of protoplanets embedded in a gaseous disc

We present a new computation of the linear tidal interaction of a protoplanetary core with a thin gaseous disc in which it is fully embedded. For the first time a discussion of the orbital evolution of cores with eccentricity (e) significantly larger than the gas-disc scale height to radius ratio (H/r) is given. We find that the direction of orbital migration reverses for e>1.1H/r. This occurs as a result of the orbital crossing of resonances in the disc that do not overlap the orbit when the eccentricity is very small. Simple expressions giving approximate fits to the eccentricity damping rate and the orbital migration rate are presented. We go on to calculate the rate of increase of the mean eccentricity for a system of protoplanetary cores due to dynamical relaxation. By equating the eccentricity damping time-scale with the dynamical relaxation time-scale we deduce that an equilibrium between eccentricity damping and excitation through scattering is attained on a 10^3 to 10^4 yr time-scale, at 1au. The equilibrium thickness of the protoplanet distribution is such that it is generally well confined within the gas disc. By use of a three dimensional N-body code we simulate the evolution of a system of protoplanetary cores, incorporating our eccentricity damping and migration rates. Assuming that collisions lead to agglomeration, we find that the vertical confinement of the protoplanet distribution permits cores to build up from 0.1 to 1 earth mass in only ~10^4 yr, within 1au. The time-scale required to achieve this is comparable to the migration time-scale. We deduce that it is not possible to build up a massive enough core to form a gas giant planet before orbital migration ultimately results in the preferential delivery of all such bodies to the neighbourhood of the central star. [Abridged]Comment: Latex in MNRAS style, 13 pages with 6 figures, also available from http://www.maths.qmw.ac.uk/~jdl

Discs and Planetary Formation

The formation, structure and evolution of protoplanetary discs is considered. The formation of giant planets within the environment of these models is also discussed.Comment: 22 pages, LaTeX (including 6 figures), uses paspconf.sty, epsf.sty and rotate.sty, to be published in Proceedings of the EC Summer School on 'Astrophysical Discs', eds J. A. Sellwood and J. Goodman, ASP Conf. Serie

Global m=1 modes and migration of protoplanetary cores in eccentric protoplanetary discs

We calculate global $m=1$ modes with low pattern speed corresponding to introducing a finite eccentricity into a protoplanetary disc. We consider disc models which are either isolated or contain one or two protoplanets orbiting in an inner cavity. Global modes that are strongly coupled to inner protoplanets are found to have disc orbits which tend to have apsidal lines antialigned with respect to those of the inner protoplanets. Other modes corresponding to free disc modes may be global over a large range of length scales and accordingly be long lived. We consider the motion of a protoplanet in the earth mass range embedded in an eccentric disc and determine the equilibrium orbits which maintain fixed apsidal alignment with respect to the disc gas orbits. Equilibrium eccentricities are found to be comparable or possibly exceed the disc eccentricity. We then approximately calculate the tidal interaction with the disc in order to estimate the orbital migration rate. Results are found to deviate from the case of axisymmetric disc with near circular protoplanet orbit once eccentricities of protoplanet and disc orbits become comparable to the disc aspect ratio in magnitude. Aligned protoplanet orbits with very similar eccentricity to that of the gas disc are found to undergo litle eccentricity change while undergoing inward migration in general. However, for significantly larger orbital eccentricities, migration may be significantly reduced or even reverse from inwards to outwards. Thus the existence of global non circular motions in discs with radial excursions comparable to the semi-thickness may have important consequences for the migration and survival of protoplanetary cores in the earth mass range.Comment: Accepted for publication by A &

A comparison of one and two dimensional models of transonic accretion discs around collapsed objects

We construct models of the inner part of a transonic adiabatic accretion disc assuming constant specific angular momentum taking the vertical structure fully into account.\\ For comparison purposes, we construct the corresponding one dimensional viscous disc models derived under vertical averaging assumptions. The conditions under which a unique location for the critical/sonic point is obtained, given an appropriate set of exterior boundary conditions for these models, is also discussed. This is not unique if the standard '$\alpha$' prescription with viscous stress proportional to the angular velocity gradient is used.\\ We use a simple model to discuss the possible limitations on the form of the viscous stress arising from the requirement that viscous information must travel at a finite speed. Contrary to results in the existing literature, the viscous stress tends to be {\it increased} rather than reduced for the type of flows we consider in which the angular momentum and angular velocity gradients have opposite signs. However, finite propagation effects may result in a unique location for the sonic point.\\ We found good agreement between the radial flow and specific angular momentum profiles in the inner regions of the one dimensional models and those in the equatorial plane for corresponding two dimensional models which may be matched for a range of $\alpha$ between 0.1 and 10$^{-4}$.Comment: (To appear in M.N.R.A.S.) 25 pages, plain LaTeX, (8 figures available upon request), QMW 93-

Conditions for the occurrence of mean-motion resonances in a low mass planetary system

The dynamical interactions that occur in newly formed planetary systems may reflect the conditions occurring in the protoplanetary disk out of which they formed. With this in mind, we explore the attainment and maintenance of orbital resonances by migrating planets in the terrestrial mass range. Migration time scales varying between millions of years and thousands of years are considered. In the former case, for which the migration time is comparable to the lifetime of the protoplanetary gas disk, a 2:1 resonance may be formed. In the latter, relatively rapid migration regime commensurabilities of high degree such as 8:7 or 11:10 may be formed. However, in any one large-scale migration several different commensurabilities may be formed sequentially, each being associated with significant orbital evolution. We also use a simple analytic theory to develop conditions for first order commensurabilities to be formed. These depend on the degree of the commensurability, the imposed migration and circularization rates, and the planet mass ratios. These conditions are found to be consistent with the results of our simulations.Comment: 11 pages with 4 figures, pdflatex, to appear in the proceedings of the conference "Extra-solar Planets in Multi-body Systems: Theory and Observations"; eds. K. Gozdziewski, A. Niedzielski and J. Schneider, EAS Publication Serie

On the evolution of mean motion resonances through stochastic forcing: Fast and slow libration modes and the origin of HD128311

Aims. We clarify the response of extrasolar planetary systems in a 2:1 mean motion commensurability with masses ranging from the super Jovian range to the terrestrial range to stochastic forcing that could result from protoplanetary disk turbulence. The behaviour of the different libration modes for a wide range of system parameters and stochastic forcing magnitudes is investigated. The growth of libration amplitudes is parameterized as a function of the relevant physical parameters. The results are applied to provide an explanation of the configuration of the HD128311 system. Methods. We first develop an analytic model from first principles without making the assumption that both eccentricities are small. We also perform numerical N-body simulations with additional stochastic forcing terms to represent the effects of putative disk turbulence. Results. Systems are quickly destabilized by large magnitudes of stochastic forcing but some stability is imparted should systems undergo a net orbital migration. The slow mode, which mostly corresponds to motion of the angle between the apsidal lines of the two planets, is converted to circulation more readily than the fast mode which is associated with oscillations of the semi-major axes. This mode is also vulnerable to the attainment of small eccentricities which causes oscillations between periods of libration and circulation. Conclusions. Stochastic forcing due to disk turbulence may have played a role in shaping the configurations of observed systems in mean motion resonance. It naturally provides a mechanism for accounting for the HD128311 system.Comment: 15 pages, 8 figures, added discussion in h and k coordinates, recommended for publicatio

On the tilting of protostellar disks by resonant tidal effects

We consider the dynamics of a protostellar disk surrounding a star in a circular-orbit binary system. Our aim is to determine whether, if the disk is initially tilted with respect to the plane of the binary orbit, the inclination of the system will increase or decrease with time. The problem is formulated in the binary frame in which the tidal potential of the companion star is static. We consider a steady, flat disk that is aligned with the binary plane and investigate its linear stability with respect to tilting or warping perturbations. The dynamics is controlled by the competing effects of the m=0 and m=2 azimuthal Fourier components of the tidal potential. In the presence of dissipation, the m=0 component causes alignment of the system, while the m=2 component has the opposite tendency. We find that disks that are sufficiently large, in particular those that extend to their tidal truncation radii, are generally stable and will therefore tend to alignment with the binary plane on a time-scale comparable to that found in previous studies. However, the effect of the m=2 component is enhanced in the vicinity of resonances where the outer radius of the disk is such that the natural frequency of a global bending mode of the disk is equal to twice the binary orbital frequency. Under such circumstances, the disk can be unstable to tilting and acquire a warped shape, even in the absence of dissipation. The outer radius corresponding to the primary resonance is always smaller than the tidal truncation radius. For disks smaller than the primary resonance, the m=2 component may be able to cause a very slow growth of inclination through the effect of a near resonance that occurs close to the disk center. We discuss these results in the light of recent observations of protostellar disks in binary systems.Comment: 21 pages, 7 figures, to be published in the Astrophysical Journa