Richard P. Nelson

Tidally Induced Gap Formation in Protostellar Disks: Gap Clearing and Suppression of Protoplanetary Growth

We present the results of numerical simulations of protostellar accretion disks that are perturbed by a protoplanetary companion that has a much smaller mass than the central object. We consider the limiting cases where the companion is in a coplanar, circular orbit and is initially embedded in the disk. Three independent numerical schemes are employed, and generic features of the flow are found in each case. In the first series of idealized models, the secondary companion is modeled as a massless, orbiting sink hole able to absorb all matter incident upon it without exerting any gravitational torque on the disk. In these simulations, accretion onto the companion induces surface density depression and gap formation centered on its orbital radius. After an initial transitory adjustment, the accretion rate onto the sink hole becomes regulated by the rate at which viscous evolution of the disk can cause matter to diffuse into the vicinity of the sink hole orbit, and thus the sink hole grows on a disk viscous timescale. In the second series of comprehensive models, the companions gravity is included. When the tidal torque exerted by the companion on the disk becomes important, the angular momentum exchange rate between the companion and the disk induces the protoplanetary accretion to reduce markedly below that in the idealized sink hole models. Whether this process is effective in inhibiting protoplanetary accretion depends on the equation of state and disk model parameters. For polytropic or isothermal equations of state, we find, in basic agreement with earlier work, that when the mean Roche lobe radius of the companion exceeds the disk thickness and when the mass ratio, q, between the companion and the central object exceeds ~40/, where is the effective Reynolds number, a clean deep gap forms. Although precise estimation is rendered difficult due to the limitation of numerical schemes in dealing with large density contrasts, the generic results of three series of simulations indicate that accretion onto sufficiently massive protoplanets can become ineffective over the expected disk lifetimes in their neighborhood. We suggest that such a process operates during planetary formation and is important in determining the final mass of giant planets.

A terrestrial planet candidate in a temperate orbit around Proxima Centauri

At a distance of 1.295 parsecs, the red dwarf Proxima Centauri (α Centauri C, GL 551, HIP 70890 or simply Proxima) is the Sun’s closest stellar neighbour and one of the best-studied low-mass stars. It has an effective temperature of only around 3,050 kelvin, a luminosity of 0.15 per cent of that of the Sun, a measured radius of 14 per cent of the radius of the Sun and a mass of about 12 per cent of the mass of the Sun. Although Proxima is considered a moderately active star, its rotation period is about 83 days (ref. 3) and its quiescent activity levels and X-ray luminosity are comparable to those of the Sun. Here we report observations that reveal the presence of a small planet with a minimum mass of about 1.3 Earth masses orbiting Proxima with a period of approximately 11.2 days at a semi-major-axis distance of around 0.05 astronomical units. Its equilibrium temperature is within the range where water could be liquid on its surface.

The migration and growth of protoplanets in protostellar discs

ABSTRA C T We investigate the gravitational interaction of a Jovian-mass protoplanet with a gaseous disc with aspect ratio and kinematic viscosity expected for the protoplanetary disc from which it formed. Different disc surface density distributions are investigated. We focus on the tidal interaction with the disc with the consequent gap formation and orbital migration of the protoplanet. Non-linear two-dimensional hydrodynamic simulations are employed using three independent numerical codes. A principal result is that the direction of the orbital migration is always inwards and such that the protoplanet reaches the central star in a near-circular orbit after a characteristic viscous time-scale of ,10 4 initial orbital periods. This is found to be independent of whether the protoplanet is allowed to accrete mass or not. Inward migration is helped by the disappearance of the inner disc, and therefore the positive torque it would exert, because of accretion on to the central star. Maximally accreting protoplanets reach about 4 Jovian masses on reaching the neighbourhood of the central star. Our results indicate that a realistic upper limit for the masses of closely orbiting giant planets is ,5 Jupiter masses, if they originate in protoplanetary discs similar to the minimum-mass solar nebula. This is because of the reduced accretion rates obtained for planets of increasing mass. Assuming that some process such as termination of the inner disc through a magnetospheric cavity stops the migration, the range of masses estimated for a number of close orbiting giant planets as well as their inward orbital migration can be accounted for by consideration of disc‐protoplanet interactions during the late stages of giant planet formation.

Planet-Disk Interaction and Orbital Evolution

As planets form and grow within gaseous protoplanetary disks, the mutual gravitational interaction between the disk and planet leads to the exchange of angular momentum and migration of the planet. We review current understanding of disk-planet interactions, focusing in particular on physical processes that determine the speed and direction of migration. We describe the evolution of low-mass planets embedded in protoplanetary disks and examine the influence of Lindblad and corotation torques as a function of the disk properties. The role of the disk in causing the evolution of eccentricities and inclinations is also discussed. We describe the rapid migration of intermediate-mass planets that may occur as a runaway process and examine the transition to gap formation and slower migration driven by the viscous evolution of the disk for massive planets. The roles and influence of disk self-gravity and magnetohydrodynamic turbulence are discussed in detail, as a function of the planet mass, as is the evolution...

A comparative study of disc–planet interaction

We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analysed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5 per cent for the Eulerian simulations. The SPH results predict the correct shape of the gap although have less resolution in the low-density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10 per cent. The spread is larger in the tidal torques acting on the planet which agree within a factor of 2 at the end of the simulation. In the Neptune case, the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.

The interaction of giant planets with a disc with MHD turbulence – IV. Migration rates of embedded protoplanets

We present the results of global cylindrical disc simulations and local shearing box simulations of protoplanets interacting with a disc undergoing magnetohydrodynamic (MHD) turbulence. The specific emphasis of this paper is to examine and quantify the magnitude of the torque exerted by the disc on the embedded protoplanets as a function of the protoplanet mass, and thus to make a first study of the induced orbital migration of protoplanets resulting from their interaction with magnetic, turbulent discs. This issue is of crucial importance in understanding the formation of gas giant planets through the so-called core instability model, and the subsequent orbital evolution post-formation prior to the dispersal of the protostellar disc. Current estimates of the migration time of protoplanetary cores in the 3–30 Earth mass range in standard disc models are τmig≃ 104–105 yr, which is much shorter than the estimated gas accretion time-scale of Jupiter-type planets. The global simulations were carried out for a disc with constant aspect ratio H/r= 0.07 and protoplanet masses of Mp= 3, 10, 30 Earth masses, and 3 Jupiter masses. The local shearing box simulations were carried out for values of the dimensionless parameter (Mp/M*)/(H/R)3= 0.1, 0.3, 1.0 and 2.0, with M*, R and H being the central mass, the orbital radius and the local disc semithickness, respectively. These allow both embedded and gap-forming protoplanets for which the disc response is non-linear to be investigated. In all cases the instantaneous net torque experienced by a protoplanet showed strong fluctuations on an orbital time-scale, and in the low-mass embedded cases oscillated between negative and positive values. Consequently, in contrast to the laminar disc type I migration scenario, orbital migration would occur as a random walk. Running time averages for embedded protoplanets over typical run times of 20–25 orbital periods indicated that the averaged torques from the inner and outer disc took on values characteristic of type I migration. However, large fluctuations occurring on longer than orbital time-scales remained, preventing convergence of the average torque to well-defined values or even to a well-defined sign for these lower mass cases. Fluctuations became relatively smaller for larger masses indicating better convergence properties, to the extent that in the 30-M⊕ simulation consistently inward, albeit noisy, migration was indicated. Both the global and local simulations showed this trend with increasing protoplanet mass, which is due to its perturbation on the disc increasing to become comparable to and then dominate the turbulence in its neighbourhood. The turbulence then becomes unable to produce very large long-term fluctuations in the torques acting on the protoplanet. Eventually gap formation occurs and there is a transition to the usual type II migration at a rate determined by the angular momentum transport in the distant parts of the disc. The existence of significant fluctuations occurring in the turbulent discs on long time-scales is an important unexplored issue for the lower mass embedded protoplanets, which are unable to modify the turbulence in their neighbourhood, and which have been studied here. If significant fluctuations occur on the longest disc evolutionary time-scales, convergence of torque running averages for practical purposes will not occur and the migration behaviour of low-mass protoplanets considered as an ensemble will be very different from predictions of type I migration theory for laminar discs. The fact that noise levels were relatively smaller in the local simulations may indicate the presence of long-term global fluctuations, but the issue remains an important one for future investigation.

TRANSITIONAL AND PRE-TRANSITIONAL DISKS: GAP OPENING BY MULTIPLE PLANETS?

We use two-dimensional hydrodynamic simulations of viscous disks to examine whether dynamically interacting multiple giant planets can explain the large gaps (spanning over one order of magnitude in radius) inferred for the transitional and pre-transitional disks around T Tauri stars. In the absence of inner disk dust depletion, we find that it requires three to four giant planets to open up large enough gaps to be consistent with inferences from spectral energy distributions, because the gap width is limited by the tendency of the planets to be driven together into 2:1 resonances. With very strong tidal torques and/or rapid planetary accretion, fewer planets can also generate a large cavity interior to the locally formed gap(s) by preventing outer disk material from moving in. In these cases, however, the reduction of surface density produces a corresponding reduction in the inner disk accretion rate onto the star; this makes it difficult to explain the observed accretion rates of the pre-transitional/transitional disks. We find that even with four planets in disks, additional substantial dust depletion is required to explain observed disk gaps/holes. Substantial dust settling and growth, with consequent significant reductions in optical depths, is inferred for typical T Tauri disks in any case, and an earlier history of dust growth is consistent with the hypothesis that pre-transitional/transitional disks are explained by the presence of giant planets. We conclude that the depths and widths of gaps and disk accretion rates in pre-transitional/transitional disks cannot be reproduced by a planet-induced gap opening scenario alone. Significant dust depletion is also required within the gaps/holes. Order-of-magnitude estimates suggest that the mass of small dust particles (1 μm) relative to the gas must be depleted to 10 −5 to 10 −2 of the interstellar medium value, implying a very efficient mechanism of small dust removal or dust growth.

The tidally induced warping, precession and truncation of accretion discs in binary systems: three-dimensional simulations

We present the results of non linear, hydrodynamic simulations, in three dimensions, of the tidal perturbation of accretion discs in binary systems where the orbit is circular and not necessarily coplanar with the disc mid–plane. The accretion discs are assumed to be geometrically thin, and of low mass relative to the stellar mass so that they are governed by thermal pressure and viscosity, but not self–gravity. The parameters that we consider in our models are the ratio of the orbital distance to the disc radius, D/R, the binary mass ratio, Ms/Mp, the initial inclination angle between the orbit and disc planes, �, and the Mach number in the outer parts of the unperturbed disc, M. Since we consider non self-gravitating discs, these calculations are relevant to protostellar binaries with separations below a few hundred astronomical units. For binary mass ratios of around unity and D/R in the range 3 to 4, we find that the global evolution of the discs is governed primarily by the value of M. For relatively low Mach numbers (i.e. M = 10 to 20) we find that the discs develop a mildly warped structure, are tidally truncated, and undergo a near rigid body precession at a rate which is in close agreement with analytical arguments. For higher Mach numbers (M ≈ 30), the evolution is towards a considerably more warped structure, but the disc none the less maintains itself as a long–lived, coherent entity. A further increase in Mach number to M = 50 leads to a dramatic disruption of the disc due to differential precession, since the sound speed is too low to allow efficient communication between the disc’s constituent parts. Additionally, it is found that the inclination angle between the disc and orbital angular momentum vectors evolves on a longer timescale which is probably the viscous evolution timescale of the disc. The calculations are relevant to a number of observed astrophysical phenomena, including the precession of jets associated with young stars, the high spectral index of some T Tauri stars, and the light curves of X–ray binaries such as Hercules X-1 which suggest the presence of precessing accretion discs.

Orbital eccentricity growth through disc-companion tidal interaction

We investigate the driving of orbital eccentricity of giant protoplanets and brown dwarfs through disc-companion tidal interactions by means of two dimensional numerical simulations. We consider disc models that are thought to be typical of protostellar discs during the planet forming epoch, with characteristic surface densities similar to standard minimum mass solar nebula models. We consider companions, ranging in mass between 1 and 30 Jupiter masses

On the orbital evolution of low mass protoplanets in turbulent, magnetised disks

M_{\rm J},