It is likely that most protostellar systems undergo a brief phase where the protostellar disc is self-gravitating. If these discs are prone to fragmentation, then they are able to rapidly form objects that are initially of several Jupiter masses and larger. The fate of these disc fragments (and the fate of planetary bodies formed afterwards via core accretion) depends sensitively not only on the fragment's interaction with the disc, but with its neighbouring fragments. We return to and revise our population synthesis model of self-gravitating disc fragmentation and tidal downsizing. Amongst other improvements, the model now directly incorporates fragment-fragment interactions while the disc is still present. We find that fragment-fragment scattering dominates the orbital evolution, even when we enforce rapid migration and inefficient gap formation. Compared to our previous model, we see a small increase in the number of terrestrial-type objects being formed, although their survival under tidal evolution is at best unclear. We also see evidence for disrupted fragments with evolved grain populations - this is circumstantial evidence for the formation of planetesimal belts, a phenomenon not seen in runs where fragment-fragment interactions are ignored. In spite of intense dynamical evolution, our population is dominated by massive giant planets and brown dwarfs at large semimajor axis, which direct imaging surveys should, but only rarely, detect. Finally, disc fragmentation is shown to be an efficient manufacturer of free floating planetary mass objects, and the typical multiplicity of systems formed via gravitational instability will be low.

In this paper, we model the chemical evolution of a 0.25 M$_{\odot}$ protoplanetary disc surrounding a 1 M$_{\odot}$ star that undergoes fragmentation due to self-gravity. We use Smoothed Particle Hydrodynamics including a radiative transfer scheme, along with time-dependent chemical evolution code to follow the composition of the disc and resulting fragments over approximately 4000 yrs. Initially, four quasi-stable fragments are formed, of which two are eventually disrupted by tidal torques in the disc. From the results of our chemical modelling, we identify species that are abundant in the fragments (e.g. H$_{\rm 2}$O, H$_{\rm 2}$S, HNO, N$_{\rm 2}$, NH$_{\rm 3}$, OCS, SO), species that are abundant in the spiral shocks within the disc (e.g. CO, CH$_{\rm 4}$, CN, CS, H$_{\rm 2}$CO), and species which are abundant in the circumfragmentary material (e.g. HCO$^{\rm +}$). Our models suggest that in some fragments it is plausible for grains to sediment to the core before releasing their volatiles into the planetary envelope, leading to changes in, e.g., the C/O ratio of the gas and ice components. We would therefore predict that the atmospheric composition of planets generated by gravitational instability should not necessarily follow the bulk chemical composition of the local disc material.

The mechanism through which meter-sized boulders grow to km-sized planetesimals in protoplanetary discs is a subject of active research, since it is critical for planet formation. To avoid spiralling into the protostar due to aerodynamic drag, objects must rapidly grow from cm-sized pebbles, which are tightly coupled to the gas, to large boulders of 1-100m in diameter. It is already well known that over-densities in the gaseous component of the disc provide potential sites for the collection of solids, and that significant density structures in the gaseous component of the disc (e.g., spiral density waves) can trap solids efficiently enough for the solid component of the disc to undergo further gravitational collapse due to their own self-gravity. In this work, we employ the PENCIL CODE to conduct local shearing sheet simulations of massive self-gravitating protoplanetary discs, to study the effect of anticyclonic transient vortices, or eddies, on the evolution of solids in these discs. We find that these types of structures are extremely efficient at concentrating small and intermediate-sized dust particles with friction times comparable to, or less than, the local orbital period of the disc. This can lead to significant over-densities in the solid component of the disc, with density enhancements comparable to, and even higher, than those within spiral density waves; increasing the rate of gravitational collapse of solids into bound structures.

We study particle dynamics in self-gravitating gaseous discs with a simple cooling law prescription via two-dimensional simulations in the shearing sheet approximation. It is well known that structures arising in the gaseous component of the disc due to a gravitational instability can have a significant effect on the evolution of dust particles. Previous results have shown that spiral density waves can be highly efficient at collecting dust particles, creating significant local over-densities of particles. The degree of such concentrations has been shown to be dependent on two parameters: the size of the dust particles and the rate of gas cooling. We expand on these findings, including the self-gravity of dust particles, to see how these particle over-densities evolve. We use the PENCIL CODE to solve the local shearing sheet equations for gas on a fixed grid together with the equations of motion for solids coupled to the gas through an aerodynamic drag force. We find that the enhancements in the surface density of particles in spiral density wave crests can reach levels high enough to allow the solid component of the disc to collapse under its own self-gravity. This produces many gravitationally bound collections of particles within the spiral structure. The total mass contained in bound structures appears nearly independent of the cooling time, suggesting that the formation of planetesimals through dust particle trapping by self-gravitating density waves may be possible at a larger range of radii within a disc than previously thought. So, density waves due to gravitational instabilities in the early stages of star formation may provide excellent sites for the rapid formation of many large, planetesimal-sized objects.

Recently it has been suggested that the fragmentation boundary in Smoothed Particle Hydrodynamic (SPH) and FARGO simulations of self-gravitating accretion discs with beta-cooling do not converge as resolution is increased. Furthermore, this recent work suggests that by carefully optimising the artificial viscosity parameters in these codes it can be shown that fragmentation may occur for much longer cooling times than earlier work suggests. If correct, this result is intriguing as it suggests that gas giant planets could form, via direct gravitational collapse, reasonably close to their parent stars. This result is, however, slightly surprising and there have been a number of recent studies suggesting that the result is likely an indication of a numerical problem with the simulations. One suggestion, in particular, is that the SPH results are influenced by the manner in which the cooling is implemented. We extend this work here and show that if the cooling is implemented in a manner that removes a known numerical artefact in the shock regions, the fragmentation boundary converges to a value consistent with earlier work and that fragmentation is unlikely for the long cooling times suggested by this recent work. We also investigate the optimisation of the artificial viscosity parameters and show that the values that appear optimal are likely introducing numerical problems in both the SPH and FARGO simulations. We therefore conclude that earlier predictions for the cooling times required for fragmentation are likely correct and that, as suggested by this earlier work, fragmentation cannot occur in the inner parts (r < 50 au) of typical protostellar discs.

It is well accepted that 'hot Jupiters' did not form in situ, as the temperature in the protoplanetary disc at the radius at which they now orbit would have been too high for planet formation to have occurred. These planets, instead, form at larger radii and then move into the region in which they now orbit. The exact process that leads to the formation of these close-in planets is, however, unclear and it seems that there may be more than one mechanism that can produce these short-period systems. Dynamical interactions in multiple-planet systems can scatter planets into highly eccentric orbits which, if the pericentre is sufficiently close to the parent star, can be tidally circularised by tidal interactions between the planet and star. Furthermore, systems with distant planetary or stellar companions can undergo Kozai cycles which can result in a planet orbiting very close to its parent star. However, the most developed model for the origin of short period planets is one in which the planet exchanges angular momentum with the surrounding protoplanetary disc and spirals in towards the central star. In the case of 'hot Jupiters', the planet is expected to open a gap in the disc and migrate through Type II .migration. If this is the dominant mechanism for producing `hot Jupiters' then we would expect the currect properties of observed close-in giant planets to be consistent with an initial population resulting from Type II migration followed by evolution due to tidal interactions with the central star. We consider initial distributions that are consistent with Type II migration and find that after tidal evolution, the final distributions can be consistent with that observed. Our results suggest that a modest initial pile-up at a ~ 0.05 au is required and that the initial eccentricity distribution must peak at e \sim 0.

We study particle dynamics in local two-dimensional simulations of self-gravitating accretion discs with a simple cooling law. It is well known that the structure which arises in the gaseous component of the disc due to a gravitational instability can have a significant effect on the evolution of dust particles. Previous results using global simulations indicate that spiral density waves are highly efficient at collecting dust particles, creating significant local over-densities which may be able to undergo gravitational collapse. We expand on these findings, using a range of cooling times to mimic the conditions at a large range of radii within the disc. Here we use the Pencil Code to solve the 2D local shearing sheet equations for gas on a fixed grid together with the equations of motion for solids coupled to the gas solely through aerodynamic drag force. We find that spiral density waves can create significant enhancements in the surface density of solids, equivalent to 1-10cm sized particles in a disc following the profiles of Clarke (2009) around a solar mass star, causing it to reach concentrations several orders of magnitude larger than the particles mean surface density. We also study the velocity dispersion of the particles, finding that the spiral structure can result in the particle velocities becoming highly ordered, having a narrow velocity dispersion. This implies low relative velocities between particles, which in turn suggests that collisions are typically low energy, lessening the likelihood of grain destruction. Both these findings suggest that the density waves that arise due to gravitational instabilities in the early stages of star formation provide excellent sites for the formation of large, planetesimal-sized objects.

Recent simulations of self-gravitating accretion discs, carried out using a three-dimensional Smoothed Particle Hydrodynamics (SPH) code by Meru and Bate, have been interpreted as implying that three-dimensional global discs fragment much more easily than would be expected from a two-dimensional local model. Subsequently, global and local two-dimensional models have been shown to display similar fragmentation properties, leaving it unclear whether the three-dimensional results reflect a physical effect or a numerical problem associated with the treatment of cooling or artificial viscosity in SPH. Here, we study how fragmentation of self-gravitating disc flows in SPH depends upon the implementation of cooling. We run disc simulations that compare a simple cooling scheme, in which each particle loses energy based upon its internal energy per unit mass, with a method in which the cooling is derived from a smoothed internal energy density field. For the simple per particle cooling scheme, we find a significant increase in the minimum cooling time scale for fragmentation with increasing resolution, matching previous results. Switching to smoothed cooling, however, results in lower critical cooling time scales, and tentative evidence for convergence at the highest spatial resolution tested. We conclude that precision studies of fragmentation using SPH require careful consideration of how cooling (and, probably, artificial viscosity) is implemented, and that the apparent non-convergence of the fragmentation boundary seen in prior simulations is likely a numerical effect. In real discs, where cooling is physically smoothed by radiative transfer effects, the fragmentation boundary is probably displaced from the two-dimensional value by a factor that is only of the order of unity.

Self-gravity becomes competitive as an angular momentum transport process in accretion discs at large radii, where the temperature is low enough that external irradiation likely contributes to the thermal balance. Irradiation is known to weaken the strength of disc self-gravity, and can suppress it entirely if the disc is maintained above the threshold for linear instability. However, its impact on the susceptibility of the disc to fragmentation is less clear. We use two-dimensional numerical simulations to investigate the evolution of self-gravitating discs as a function of the local cooling time and strength of irradiation. In the regime where the disc does not fragment, we show that local thermal equilibrium continues to determine the stress - which can be represented as an effective viscous alpha - out to very long cooling times (at least 240 dynamical times). In this regime, the power spectrum of the perturbations is uniquely set by the effective viscous alpha and not by the cooling rate. Fragmentation occurs for cooling times tau < beta_crit / Omega, where beta_crit is a weak function of the level of irradiation. We find that beta_crit declines by approximately a factor of two, as irradiation is increased from zero up to the level where instability is almost quenched. The numerical results imply that irradiation cannot generally avert fragmentation of self-gravitating discs at large radii; if other angular momentum transport sources are weak mass will build up until self-gravity sets in, and fragmentation will ensue.

Motivated by the recent results of \citet{Lesur_Ogilvie10} on the transport properties of incompressible convection in protoplanetary discs, in this paper we study the role of compressibility and hence of another basic mode -- spiral density waves -- in convective instability in discs. We analyse the linear dynamics of non-axisymmetric convection and spiral density waves in a Keplerian disc with superadiabatic vertical stratification using the local shearing box approach. It is demonstrated that the shear associated with Keplerian differential rotation introduces a novel phenomenon, it causes these two perturbation modes to become coupled: during evolution the convective mode generates (trailing) spiral density waves and can therefore be regarded as a new source of spiral density waves in discs. The wave generation process studied here owes its existence solely to shear of the disc's differential rotation, and is a special manifestation of a more general linear mode coupling phenomena universally taking place in flows with an inhomogeneous velocity profile. We quantify the efficiency of spiral density wave generation by convection as a function of azimuthal and vertical wavenumbers of these modes and find that it is maximal and most powerful when both these length-scales are comparable to the disc scale height. We also show that unlike the convective mode, which tends to transport angular momentum inwards in the linear regime, the spiral density waves transport angular momentum outwards. Based on these findings, we suggest that in the non-linear regime spiral density waves generated by convection may play a role in enhancing the transport of angular momentum due the convective mode alone, which is actually being changed to outward by non-linearity, as indicated by above-mentioned recent developments.

It has recently been noted that many discs around T Tauri stars appear to comprise only a few Jupiter-masses of gas and dust. Using millimetre surveys of discs within six local star-formation regions, we confirm this result, and find that only a few percent of young stars have enough circumstellar material to build gas giant planets, in standard core accretion models. Since the frequency of observed exo-planets is greater than this, there is a `missing mass' problem. As alternatives to simply adjusting the conversion of dust-flux to disc mass, we investigate three other classes of solution. Migration of planets could hypothetically sweep up the disc mass reservoir more efficiently, but trends in multi-planet systems do not support such a model, and theoretical models suggest that the gas accretion timescale is too short for migration to sweep the disc. Enhanced inner-disc mass reservoirs are possible, agreeing with predictions of disc evolution through self-gravity, but not adding to millimetre dust-flux as the inner disc is optically thick. Finally, the incidence of massive discs is shown to be higher at the {\it proto}stellar stages, Classes 0 and I, where discs substantial enough to form planets via core accretion are abundant enough to match the frequency of exo-planets. Gravitational instability may also operate in the Class 0 epoch, where half the objects have potentially unstable discs of $\ga$30 % of the stellar mass. However, recent calculations indicate that forming gas giants inside 50 AU by instability is unlikely, even in such massive discs. Overall, the results presented suggest that the canonically 'proto-planetary' discs of Class II T Tauri stars {\bf have globally low masses in dust observable at millimetre wavelengths, and conversion to larger bodies (anywhere from small rocks up to planetary cores) must already have occurred.}

We perform linear analysis of axisymmetric vertical normal modes in stratified compressible self-gravitating polytropic discs in the shearing box approximation. We study specific dynamics for subadiabatic, adiabatic and superadiabatic vertical stratifications. In the absence of self-gravity, four well-known principal modes can be identified in a stratified disc: acoustic p-, surface gravity f-, buoyancy g- and inertial r-modes. After characterizing modes in the non-self-gravitating case, we include self-gravity and investigate how it modifies the properties of these modes. We find that self-gravity, to a certain degree, reduces their frequencies and changes the structure of the dispersion curves and eigenfunctions at radial wavelengths comparable to the disc height. Its influence on the basic branch of the r-mode, in the case of subadiabatic and adiabatic stratifications, and on the basic branch of the g-mode, in the case of superadiabatic stratification (which in addition exhibits convective instability), does appear to be strongest. Reducing the three-dimensional Toomre's parameter Q_{3D} results in the latter modes becoming unstable due to self-gravity, so that they determine the onset criterion and nature of gravitational instability of a stratified disc. By contrast, the p-, f- and convectively stable g-modes, although their corresponding \omega^2 are reduced by self-gravity, never become unstable however small the value of Q_{3D}. This is a consequence of the three-dimensionality of the disc. The eigenfunctions corresponding to the gravitationally unstable modes are intrinsically 3D. We also contrast the more exact instability criterion based on our 3D model with that of density waves in 2D (razor-thin) discs. Based on these findings, we comment on the origin of surface distortions seen in numerical simulations of self-gravitating discs.

We use time-dependent, one-dimensional disc models to investigate the evolution of protostellar discs that form through the collapse of molecular cloud cores and in which the primary transport mechanism is self-gravity. We assume that these discs settle into a state of thermal equilibrium with Q = 2 and that the strength of the angular momentum transport is set by the cooling rate of the disc. The results suggest that these discs will attain a quasi-steady state that persists for a number of free-fall times and in which most of the mass within 100 au is located inside 10-20 au. This pile-up of mass in the inner disc could result in temperatures that are high enough for the growth of MHD turbulence which could rapidly drain the inner disc and lead to FU Orionis-like outbursts. In all our simulations, the inner regions of the discs (r < 40 au) were stable against fragmentation, while fragmentation was possible in the outer regions (r > 40 au) of discs that formed from cores that had enough initial angular momentum to deposit sufficient mass in these outer regions. The large amounts of mass in these outer regions, however, suggests that fragmentation will lead to the formation of sub-stellar and stellar mass companions, rather than planetary mass objects. Although mass accretion rates were largely consistent with observations, the large disc masses suggest that an additional transport mechanism (such as MRI occuring in the upper layers of the disc) must operate in order to drain the remaining disc material within observed disc lifetimes.

Angular momentum transport within young massive protoplanetary discs may be dominated by self-gravity at radii where the disk is too weakly ionized to allow the development of the magneto-rotational instability. We use time-dependent one-dimensional disc models, based on a local cooling time calculation of the efficiency of transport, to study the radial structure and stability (against fragmentation) of protoplanetary discs in which self-gravity is the sole transport mechanism. We find that self-gravitating discs rapidly attain a quasi-steady state in which the surface density in the inner disc is high and the strength of turbulence very low (alpha ~ 10^{-3} or less inside 5 au). Temperatures high enough to form crystalline silicates may extend out to several au at early times within these discs. None of our discs spontaneously develop regions that would be unambiguously unstable to fragmentation into substellar objects, though the outer regions (beyond 20 au) of the most massive discs are close enough to the threshold that fragmentation cannot be ruled out. We discuss how the mass accretion rates through such discs may vary with disc mass and with mass of the central star, and note that a determination of the \dot{M}-M_* relation for very young systems may allow a test of the model.

(abridged) Vortices are believed to play a role in the formation of km-sized planetesimals. However, vortex dynamics is commonly studied in non-self-gravitating discs. The main goal here is to examine the effects of disc self-gravity on vortex dynamics. For this purpose, we employ the 2D self-gravitating shearing sheet approximation. A simple cooling law with a constant cooling time is adopted, such that the disc settles down into a quasi-steady gravitoturbulent state. In this state, vortices appear as transient structures undergoing recurring phases of formation, growth to sizes comparable to a local Jeans scale and eventual shearing and destruction due to the combined effects of self-gravity and background Keplerian shear. Each phase typically lasts about 2 orbital periods or less. As a result, in self-gravitating discs the overall dynamical picture of vortex evolution is irregular consisting of many transient vortices at different evolutionary stages and, therefore, with various sizes up to the local Jeans scale. Vortices generate density waves during evolution, which turn into shocks. Therefore, the dynamics of density waves and vortices are coupled implying that, in general, one should consider both vortex and spiral density wave modes in order to get a proper understanding of self-gravitating disc dynamics. Our results suggest that given such an irregular and rapidly varying character of vortex evolution in self-gravitating discs, it may be difficult for such vortices to effectively trap dust particles. Further study of the behaviour of dust particles embedded in a self-gravitating gaseous disc is, however, required to strengthen this conclusion.

The presence of young massive stars orbiting on eccentric rings within a few tenths of a parsec of the supermassive black hole in the Galactic centre is challenging for theories of star formation. The high tidal shear from the black hole should tear apart the molecular clouds that form stars elsewhere in the Galaxy, while transporting the stars to the Galactic centre also appears unlikely during their stellar lifetimes. We present numerical simulations of the infall of a giant molecular cloud that interacts with the black hole. The transfer of energy during closest approach allows part of the cloud to become bound to the black hole, forming an eccentric disc that quickly fragments to form stars. Compressional heating due to the black hole raises the temperature of the gas to 100-1000K, ensuring that the fragmentation produces relatively high stellar masses. These stars retain the eccentricity of the disc and, for a sufficiently massive initial cloud, produce an extremely top-heavy distribution of stellar masses. This potentially repetitive process can therefore explain the presence of multiple eccentric rings of young stars in the presence of a supermassive black hole.

Vortices are believed to greatly help the formation of km sized planetesimals by collecting dust particles in their centers. However, vortex dynamics is commonly studied in non-self-gravitating disks. The main goal here is to examine the effects of disk self-gravity on the vortex dynamics via numerical simulations. In the self-gravitating case, when quasi-steady gravitoturbulent state is reached, vortices appear as transient structures undergoing recurring phases of formation, growth to sizes comparable to a local Jeans scale, and eventual shearing and destruction due to gravitational instability. Each phase lasts over 2-3 orbital periods. Vortices and density waves appear to be coupled implying that, in general, one should consider both vortex and density wave modes for a proper understanding of self-gravitating disk dynamics. Our results imply that given such an irregular and rapidly changing, transient character of vortex evolution in self-gravitating disks it may be difficult for such vortices to effectively trap dust particles in their centers that is a necessary process towards planet formation.

We have imaged the disc of the young star HL Tau using the VLA at 1.3 cm, with 0.08" resolution (as small as the orbit of Jupiter). The disc is around half the stellar mass, assuming a canonical gas-mass conversion from the measured mass in large dust grains. A simulation shows that such discs are gravitationally unstable, and can fragment at radii of a few tens of AU to form planets. The VLA image shows a compact feature in the disc at 65 AU radius (confirming the `nebulosity' of Welch et al. 2004), which is interpreted as a localised surface density enhancement representing a candidate proto-planet in its earliest accretion phase. If correct, this is the first image of a low-mass companion object seen together with the parent disc material out of which it is forming. The object has an inferred gas plus dust mass of approximately 14 M(Jupiter), similar to the mass of a proto-planet formed in the simulation. The disc instability may have been enhanced by a stellar flyby: the proper motion of the nearby star XZ Tau shows it could have recently passed the HL Tau disc as close as ~600 AU.

We use numerical simulations to model the migration of massive planets at small radii and compare the results with the known properties of 'hot Jupiters' (extrasolar planets with semi-major axes a < 0.1 AU). For planet masses Mp sin i > 0.5 MJup, the evidence for any `pile-up' at small radii is weak (statistically insignificant), and although the mass function of hot Jupiters is deficient in high mass planets as compared to a reference sample located further out, the small sample size precludes definitive conclusions. We suggest that these properties are consistent with disc migration followed by entry into a magnetospheric cavity close to the star. Entry into the cavity results in a slowing of migration, accompanied by a growth in orbital eccentricity. For planet masses in excess of 1 Jupiter mass we find eccentricity growth timescales of a few x 10^5 years, suggesting that these planets may often be rapidly destroyed. Eccentricity growth appears to be faster for more massive planets which may explain changes in the planetary mass function at small radii and may also predict a pile-up of lower mass planets, the sample of which is still incomplete.

The quadruple young stellar system HD 98800 consists of two spectroscopic binary pairs with a circumbinary disk around the B component. Recent work by Boden and collaborators using infrared interferometry and radial velocity data resulted in a determination of the physical orbit for HD 98800 B. We use the resulting inclination of the binary and the measured extinction toward the B component stars to constrain the distribution of circumbinary material. Although a standard optically and geometrically thick disk model can reproduce the spectral energy distribution, it can not account for the observed extinction if the binary and the disk are co-planar. We next constructed a dynamical model to investigate the influence of the A component, which is not in the Ba-Bb orbital plane, on the B disk. We find that these interactions have a substantial impact on the inclination of the B circumbinary disk with respect to the Ba-Bb orbital plane. The resulting warp would be sufficient to place material into the line of sight and the non-coplanar disk orientation may also cause the upper layers of the disk to intersect the line of sight if the disk is geometrically thick. These simulations also support that the dynamics of the Ba-Bb orbit clear the inner region to a radius of ~3 AU. We then discuss whether the somewhat unusual properties of the HD 98800 B disk are consistent with material remnant from the star formation process or with more recent creation by collisions from larger bodies.

In this paper, we use high-resolution smoothed particle hydrodynamics (SPH) simulations to investigate the response of a marginally stable self-gravitating protostellar disc to a close parabolic encounter with a companion discless star. Our main aim is to test whether close brown dwarfs or massive planets can form out of the fragmentation of such discs. We follow the thermal evolution of the disc by including the effects of heating due to compression and shocks and a simple prescription for cooling and find results that contrast with previous isothermal simulations. In the present case we find that fragmentation is inhibited by the interaction, due to the strong effect of tidal heating, which results in a strong stabilization of the disc. A similar behaviour was also previously observed in other simulations involving discs in binary systems. As in the case of isolated discs, it appears that the condition for fragmentation ultimately depends on the cooling rate.

The spectral energy distributions (SEDs) of some T Tauri stars display a deficit of near-IR flux that could be a consequence of an embedded Jupiter-mass planet partially clearing an inner hole in the circumstellar disc. Here, we use two-dimensional numerical simulations of the planet-disc interaction, in concert with simple models for the dust dynamics, to quantify how a planet influences the dust at different radii within the disc. We show that pressure gradients at the outer edge of the gap cleared by the planet act as a filter - letting particles smaller than a critical size through to the inner disc while holding back larger particles in the outer disc. The critical particle size depends upon the disc properties, but is typically of the order of 10 microns. This filtration process will lead to discontinuous grain populations across the planet's orbital radius, with small grains in the inner disc and an outer population of larger grains. We show that this type of dust population is qualitatively consistent with SED modelling of systems that have optically thin inner holes in their circumstellar discs. This process can also produce a very large gas-to-dust ratio in the inner disc, potentially explaining those systems with optically thin inner cavities that still have relatively high accretion rates.

An unsolved issue in the standard core accretion model for gaseous planet formation is how kilometre-sized planetesimals form from, initially, micron-sized dust grains. Solid growth beyond metre sizes can be difficult both because the sticking efficiency becomes very small, and because these particles should rapidly migrate into the central star. We consider here how metre-sized particles evolve in self-gravitating accretion discs using simulations in which the gravitational influence of the solid particles is also included. Metre-sized particles become strongly concentrated in the spiral structures present in the disc and, if the solid to gas density ratio is sufficiently high, can fragment due to their own self-gravity to form planetesimals directly. This result suggests that planetesimal formation may occur very early in the star formation process while discs are still massive enough to be self-gravitating. The dependence of this process on the surface density of the solids is also consistent with the observation that extrasolar planets are preferentially found around high metallicity stars.

The issue of fragmentation in self-gravitating gaseous accretion discs has implications both for the formation of stars in discs in the nuclei of active galaxies, and for the formation of gaseous planets or brown dwarfs in circumstellar discs. It is now well established that fragmentation occurs if the disc is cooled on a timescale smaller than the local dynamical timescale, while for longer cooling times the disc reaches a quasi-steady state in thermal equilibrium, with the cooling rate balanced by the heating due to gravitational stresses. We investigate here how the fragmentation boundary depends on the assumed equation of state. We find that the cooling time required for fragmentation increases as the specific heat ratio, gamma, decreases, exceeding the local dynamical timescale for gamma = 7/5. This result can be easily interpreted as a consequence of there being a maximum stress (in units of the local disc pressure) that can be sustained by a self-gravitating disc in quasi-equilibrium. Fragmentation occurs if the cooling time is such that the stress required to reach thermal equilibrium exceeds this value, independent of gamma. This result suggest that a quasi-steady, self-gravitating disc can never produce a stress that results in the viscous alpha parameter exceeding ~ 0.06.

Gravitational torques between a planet and gas in the protoplanetary disk result in orbital migration of the planet, and are likely to play an important role in the formation and early evolution of planetary systems. For masses comparable to those of observed giant extrasolar planets, the interaction with the disk is strong enough to form a gap, leading to coupled evolution of the planet and disk on a viscous time scale (Type II migration). Both the existence of hot Jupiters, and the statistical distribution of observed orbital radii, are consistent with an important role for Type II migration in the history of currently observed systems. We discuss the possibility of improving constraints on migration by including information on the host stars' metallicity, and note that migration could also form a population of massive planets at large orbital radii that may be indirectly detected via their influence on debris disks. For lower mass planets with masses of the order of that of the Earth, surface density perturbations created by the planet are small, and migration in a laminar disk is driven by an intrinsic and apparently robust asymmetry between interior and exterior torques. Analytic and numerical calculations of this Type I migration are in reasonable accord, and predict rapid orbital decay during the final stages of the formation of giant planet cores. The difficulty of reconciling Type I migration with giant planet formation may signal basic errors in our understanding of protoplanetary disks, core accretion, or both. We discuss physical effects that might alter Type I behavior, in particular the possibility that for sufficiently low masses turbulent fluctuations in the gas surface density dominate the torque, leading to random walk migration of very low mass bodies.

We investigate how the statistical distribution of extrasolar planets may be combined with knowledge of the host stars' metallicity to yield constraints on the migration histories of gas giant planets. At any radius, planets that barely manage to form around the lowest metallicity stars accrete their envelopes just as the gas disk is being dissipated, so the lower envelope of planets in a plot of metallicity vs semi-major axis defines a sample of non-migratory planets that will have suffered less than average migration subsequent to gap opening. Under the assumption that metallicity largely controls the initial surface density of planetesimals, we use simplified core accretion models to calculate how the minimum metallicity needed for planet formation varies as a function of semi-major axis. Models that do not include core migration prior to gap opening (Type I migration) predict that the critical metallicity is largely flat between the snow line and a semimajor axis of about 6 AU, with a weak dependence on the initial surface density profile of planetesimals. When slow Type I migration is included, the critical metallicity is found to increase steadily from 1-10 AU. Large planet samples, that include planets at modestly greater orbital radii than present surveys, therefore have the potential to quantify the extent of migration in both Type I and Type II regimes.

In this paper, we extend our previous analysis (Lodato & Rice 2004) of the transport properties induced by gravitational instabilities in cooling, gaseous accretion discs to the case where the disc mass is comparable to the central object. In order to do so, we have performed global, three-dimensional smoothed particle hydrodynamics simulations of massive discs. These new simulations show a much more complex temporal evolution with respect to the less massive case. Whereas in the low disc mass case a self-regulated, marginally stable state (characterized by an approximately constant radial profile of the stability parameter $Q$) is easily established, in the high disc mass case we observe the development of an initial transient and subsequent settling down in a self-regulated state in some simulations, or a series or recurrent spiral episodes, with low azimuthal wave number $m$, in others. Accretion in this last case can therefore be a highly variable process. On the other hand, we find that the secular evolution of the disc is relatively slow. In fact, the time-average of the stress induced by self-gravity results in accretion time-scales much longer than the dynamical timescale, in contrast with previous isothermal simulations of massive accretion discs. We have also compared the resulting stress tensor with the expectations based on a local theory of transport, finding no significant evidence for global wave energy transport.

In this paper we consider the evolution of small planetesimals in marginally stable, self-gravitating protoplanetary discs. The drag force between the disc gas and the embedded planetesimals generally causes the planetesimals to drift inwards through the disc at a rate that depends on the particle size. In a marginally stable, self-gravitating disc, however, the planetesimals are significantly influenced by the non-axisymmetric spiral structures resulting from the growth of the gravitational instability. For small particles, that are strongly coupled to the disc gas, and for large particles, that have essentially decoupled from the disc gas, the effect is not particularly significant. Intermediate sized particles, which would generally have the largest radial drift rates, do, however, become significantly concentrated at the peaks of the spiral arms. These high density regions may persist for, of order, an orbital period and may attain densities comparable to that of the disc gas. Although at the end of the simulation only ~ 25 % of the planetesimal particles lie in regions of enhanced density, during the course of the simulation at least 75 % of the planetesimal particles have at some stage been in a such a region. We find that the concentration of particles in the spiral arms results in an increased collision rate, an effect that could significantly accelerate planetesimal growth. The density enhancements may also be sufficient for the growth of planetesimals through direct gravitational collapse. The interaction between small planetesimals and self-gravitating spiral structures may therefore play an important role in the formation of large planetesimals that will ultimately coagulate to form terrestrial planets or the cores of gas/ice giant planets.

In this paper we examine the issue of characterising the transport associated with gravitational instabilities in relatively cold discs, discussing in particular the conditions under which it can be described within a local, viscous framework. We present the results of global, three-dimensional, SPH simulations of self-gravitating accretion discs, in which the disc is cooled using a simple parametrisation for the cooling function. Our simulations show that the disc settles in a ``self-regulated'' state, where the axisymmetric stability parameter $Q\approx 1$ and where transport and energy dissipation are dominated by self-gravity. We have computed the gravitational stress tensor and compared our results with expectations based on a local theory of transport. We find that, as long as the disc mass is smaller than $0.25M_{\star}$ and the aspect ratio $H/R\lesssim 0.1$, transport is determined locally, thus allowing for a viscous treatment of the disc evolution.

Self-gravitating protostellar discs are unstable to fragmentation if the gas can cool on a time scale that is short compared to the orbital period. We use a combination of hydrodynamic simulations and N-body orbit integrations to study the long term evolution of a fragmenting disc with an initial mass ratio to the star of M_disc/M_star = 0.1. For a disc which is initially unstable across a range of radii, a combination of collapse and subsequent accretion yields substellar objects with a spectrum of masses extending (for a Solar mass star) up to ~0.01 M_sun. Subsequent gravitational evolution ejects most of the lower mass objects within a few million years, leaving a small number of very massive planets or brown dwarfs in eccentric orbits at moderately small radii. Based on these results, systems such as HD 168443 -- in which the companions are close to or beyond the deuterium burning limit -- appear to be the best candidates to have formed via gravitational instability. If massive substellar companions originate from disc fragmentation, while lower-mass planetary companions originate from core accretion, the metallicity distribution of stars which host massive substellar companions at radii of ~1 au should differ from that of stars with lower mass planetary companions.

Numerical simulations show that the migration of growing planetary cores may be dominated by turbulent fluctuations in the protoplanetary disk, rather than by any mean property of the flow. We quantify the impact of this stochastic core migration on the formation time scale and core mass of giant planets at the onset of runaway gas accretion. For standard Solar Nebula conditions, the formation of Jupiter can be accelerated by almost an order of magnitude if the growing core executes a random walk with an amplitude of a few tenths of an au. A modestly reduced surface density of planetesimals allows Jupiter to form within 10 Myr, with an initial core mass below 10 Earth masses, in better agreement with observational constraints. For extrasolar planetary systems, the results suggest that core accretion could form massive planets in disks with lower metallicities, and shorter lifetimes, than the Solar Nebula.

In this paper we examine the issue of characterizing the transport associated with gravitational instabilities in relatively cold discs, discussing in particular under which condition it can be described within a local, viscous framework. We present the results of global, three-dimensional, SPH simulations of self-gravitating accretion discs, in which the disc is cooled using a simple parameterization for the cooling function. Our simulations show that the disc settles in a ``self-regulated'' state, where the axisymmetric stability parameter $Q\approx 1$ and where transport and energy dissipation are dominated by self-gravity. We have computed the gravitational stress tensor and compared our results with expectations based on a local theory of transport. We find that, for disc masses smaller than $0.25M_{\star}$ and aspect ratio $H/r\lesssim 0.1$, transport is determined locally, thus allowing for a viscous treatment of the disc evolution.

It has already been shown, using a local model, that accretion discs with cooling times t_cool <= 3 Omega^-1 fragment into gravitationally bound objects, while those with cooling times t_cool > 3 Omega^-1 evolve into a quasi-steady state. We present results of three-dimensional simulations that test if the local result still holds globally. We find that the fragmentation boundary is close to that determined using the local model, but that fragmentation may occur for longer cooling times when the disc is more massive or when the mass is distributed in such a way as to make a particular region of the disc more susceptible to the growth of the gravitational instability. These results have significant implications for the formation of of gaseous planets in protoplanetary discs and also for the redistribution of angular momentum which could be driven by the presence of relatively massive, bound objects within the disc.

The unusual spectral energy distribution (SED) of the classical T Tauri star GM Aurigae provides evidence for the presence of an inner disc hole extending to several au. Using a combination of hydrodynamical simulations and Monte Carlo radiative transport, we investigate whether the observed SED is consistent with the inner hole being created and maintained by an orbiting planet. We show that an ~ 2 M_Jupiter planet, orbiting at 2.5 au in a disc with mass 0.047 M_sun and radius 300 au, provides a good match both to the SED and to CO observations which constrain the velocity field in the disc. A range of planet masses is allowed by current data, but could in principle be distinguished with further observations between 3 and ~ 20 microns. Future high precision astrometric instruments should also be able to detect the motion of the central star due to an orbiting Jupiter mass planet. We argue that the small number of T Tauri stars with SEDs resembling that of GM Aur is broadly consistent with the expected statistics of embedded migrating planets.

Using a local model Gammie (2001) has shown that accretion discs with cooling times t_cool <= 3/Omega fragment into gravitationally bound objects, while those with cooling times t_cool > 3/Omega evolve into a quasi-steady state. We use three-dimensional smoothed particle hydrodynamic simulations of protoplanetary accretion discs to test if the local results hold globally. We find that for disc masses appropriate for T Tauri discs, the fragmentation boundary still occurs at a cooling time close to t_cool = 3/Omega. For more massive discs, which are likely to be present at an earlier stage of the star formation process, fragmentation occurs for longer cooling times, but still within a factor of two of that predicted using a local model. These results have implications not only for planet formation in protoplanetary discs and star formation in AGN discs, but also for the redistribution of angular momentum which could be driven by the presence of relatively massive objects within the accretion disc.

We use high resolution numerical simulations to study whether gravitational instabilities within circumstellar discs can produce astrometrically detectable motion of the central star. For discs with masses of M_disc = 0.1 M_star, which are permanantly stable against fragmentation, we find that the magnitude of the astrometric signal depends upon the efficiency of disc cooling. Short cooling times produce prominent filamentary spiral structures in the disc, and lead to stellar motions that are potentially observable with future high precision astrometric experiments. For a disc that is marginally unstable within radii of \~10 au, we estimate astrometric displacements of 10-100 microarcsec on decade timescales for a star at a distance of 100 pc. The predicted displacement is suppressed by a factor of several in more stable discs in which the cooling time exceeds the local dynamical time by an order of magnitude. We find that the largest contribution comes from material in the outer regions of the disc and hence, in the most pessimistic scenario, the stellar motions caused by the disc could confuse astrometric searches for low mass planets orbiting at large radii. They are, however, unlikely to present any complications in searches for embedded planets orbiting at small radii, relative to the disc size, or Jupiter mass planets or greater orbiting at large radii.