We present a 3D hydrodynamics study of gravitational instabilities (GIs) in a 0.14 M$_{\odot}$ circumbinary protoplanetary disc orbiting a 1 M$_{\odot}$ star and a 0.02 M$_{\odot}$ brown dwarf companion. We examine the thermodynamical state of the disc and determine the strengths of GI-induced density waves, nonaxisymmetric density structures, mass inflow and outflow, and gravitational torques. Results are compared with a parallel simulation of a protoplanetary disc without the brown dwarf binary companion. Simulations are performed using CHYMERA, a radiative 3D hydrodynamics code. The onset of GIs in the circumbinary disc is much more violent due to the stimulation of a strong one-armed density wave by the brown dwarf. Despite this early difference, detailed analyses show that both discs relax to a very similar quasi-steady phase by 2,500 years after the beginning of the simulations. Similarities include the thermodynamics of the quasi-steady phase, the final surface density distribution, radial mass influx, and nonaxisymmetric power and torques for spiral arm multiplicities of two or more. Effects of binarity in the disc are evident in gravitational torque profiles, temperature profiles in the inner discs, and radial mass transport. After 3,800 years, the semimajor axis of the binary decreases by about one percentage and the eccentricity roughly doubles. The mass transport in the outer circumbinary disc associated with the one-armed wave may influence planet formation.

Gravitational instabilities (GIs) are most likely a fundamental process during the early stages of protoplanetary disc formation. Recently, there have been detections of spiral features in young, embedded objects that appear consistent with GI-driven structure. It is crucial to perform hydrodynamic and radiative transfer simulations of gravitationally unstable discs in order to assess the validity of GIs in such objects, and constrain optimal targets for future observations. We utilise the radiative transfer code LIME to produce continuum emission maps of a $0.17\,\mathrm{M}_{\odot}$ self-gravitating protosolar-like disc. We note the limitations of using LIME as is and explore methods to improve upon the default gridding. We use CASA to produce synthetic observations of 270 continuum emission maps generated across different frequencies, inclinations and dust opacities. We find that the spiral structure of our protosolar-like disc model is distinguishable across the majority of our parameter space after 1 hour of observation, and is especially prominent at 230$\,$GHz due to the favourable combination of angular resolution and sensitivity. Disc mass derived from the observations is sensitive to the assumed dust opacities and temperatures, and therefore can be underestimated by a factor of at least 30 at 850$\,$GHz and 2.5 at 90$\,$GHz. As a result, this effect could retrospectively validate GIs in discs previously thought not massive enough to be gravitationally unstable, which could have a significant impact on the understanding of the formation and evolution of protoplanetary discs.

To date, most simulations of the chemistry in protoplanetary discs have used 1+1D or 2D axisymmetric $\alpha$-disc models to determine chemical compositions within young systems. This assumption is inappropriate for non-axisymmetric, gravitationally unstable discs, which may be a significant stage in early protoplanetary disc evolution. Using 3D radiative hydrodynamics, we have modelled the physical and chemical evolution of a 0.17 M$_{\odot}$ self-gravitating disc over a period of 2000 yr. The 0.8 M$_{\odot}$ central protostar is likely to evolve into a solar-like star, and hence this Class 0 or early Class I young stellar object may be analogous to our early Solar System. Shocks driven by gravitational instabilities enhance the desorption rates, which dominate the changes in gas-phase fractional abundances for most species. We find that at the end of the simulation, a number of species distinctly trace the spiral structure of our relatively low-mass disc, particularly CN. We compare our simulation to that of a more massive disc, and conclude that mass differences between gravitationally unstable discs may not have a strong impact on the chemical composition. We find that over the duration of our simulation, successive shock heating has a permanent effect on the abundances of HNO, CN and NH$_3$, which may have significant implications for both simulations and observations. We also find that HCO$^+$ may be a useful tracer of disc mass. We conclude that gravitational instabilities induced in lower mass discs can significantly, and permanently, affect the chemical evolution, and that observations with high-resolution instruments such as ALMA offer a promising means of characterising gravitational instabilities in protosolar discs.

The formation and earliest stages of protoplanetary discs remain poorly constrained by observations. ALMA will soon revolutionise this field. Therefore, it is important to provide predictions which will be valuable for the interpretation of future high sensitivity and high angular resolution observations. Here we present simulated ALMA observations based on radiative transfer modelling of a relatively massive (0.39 M_solar) self-gravitating disc embedded in a 10 M_solar dense core, with structure similar to the pre-stellar core L1544. We focus on simple species and conclude that C17O 3-2, HCO+ 3-2, OCS 26-25 and H2CO 404-303 lines can be used to probe the disc structure and kinematics at all scales.

Until now, axisymmetric, alpha-disc models have been adopted for calculations of the chemical composition of protoplanetary discs. While this approach is reasonable for many discs, it is not appropriate when self-gravity is important. In this case, spiral waves and shocks cause temperature and density variations that affect the chemistry. We have adopted a dynamical model of a solar-mass star surrounded by a massive (0.39 Msun), self-gravitating disc, similar to those that may be found around Class 0 and early Class I protostars, in a study of disc chemistry. We find that for each of a number of species, e.g. H2O, adsorption and desorption dominate the changes in the gas-phase fractional abundance; because the desorption rates are very sensitive to temperature, maps of the emissions from such species should reveal the locations of shocks of varying strengths. The gas-phase fractional abundances of some other species, e.g. CS, are also affected by gas-phase reactions, particularly in warm shocked regions. We conclude that the dynamics of massive discs have a strong impact on how they appear when imaged in the emission lines of various molecular species.

Characterization of migration in gravitationally unstable disks is necessary to understand the fate of protoplanets formed by disk instability. As part of a larger study, we are using a 3D radiative hydrodynamics code to investigate how an embedded gas giant planet interacts with a gas disk that undergoes gravitational instabilities (GIs). This Letter presents results from simulations with a Jupiter-mass planet placed in orbit at 25 AU within a 0.14 $M_{\odot}$ disk. The disk spans 5 to 40 AU around a 1 $M_{\odot}$ star and is initially marginally unstable. In one simulation, the planet is inserted prior to the eruption of GIs; in another, it is inserted only after the disk has settled into a quasi-steady GI-active state, where heating by GIs roughly balances radiative cooling. When the planet is present from the beginning, its own wake stimulates growth of a particular global mode with which it strongly interacts, and the planet plunges inward six AU in about 10$^3$ years. In both cases with embedded planets, there are times when the planet's radial motion is slow and varies in direction. At other times, when the planet appears to be interacting with strong spiral modes, migration both inward and outward can be relatively rapid, covering several AUs over hundreds of years. Migration in both cases appears to stall near the inner Lindblad resonance of a dominant low-order mode. Planet orbit eccentricities fluctuate rapidly between about 0.02 to 0.1 throughout the GI-active phases of the simulations.

We investigate the coupling between rock-size solids and gas during the formation of gas giant planets by disk fragmentation in the outer regions of massive disks. In this study, we use three-dimensional radiative hydrodynamics simulations and model solids as a spatial distribution of particles. We assume that half of the total solid fraction is in small grains and half in large solids. The former are perfectly entrained with the gas and set the opacity in the disk, while the latter are allowed to respond to gas drag forces, with the back reaction on the gas taken into account. To explore the maximum effects of gas-solid interactions, we first consider 10cm-size particles. We then compare these results to a simulation with 1 km-size particles, which explores the low-drag regime. We show that (1) disk instability planets have the potential to form large cores due to aerodynamic capturing of rock-size solids in spiral arms before fragmentation; (2) that temporary clumps can concentrate tens of $M_{\oplus}$ of solids in very localized regions before clump disruption; (3) that the formation of permanent clumps, even in the outer disk, is dependent on the grain-size distribution, i.e., the opacity; (4) that nonaxisymmetric structure in the disk can create disk regions that have a solids-to-gas ratio greater than unity; (5) that the solid distribution may affect the fragmentation process; (6) that proto-gas giants and proto-brown dwarfs can start as differentiated objects prior to the H$_2$ collapse phase; (7) that spiral arms in a gravitationally unstable disk are able to stop the inward drift of rock-size solids, even redistributing them to larger radii; and, (8) that large solids can form spiral arms that are offset from the gaseous spiral arms. We conclude that planet embryo formation can be strongly affected by the growth of solids during the earliest stages of disk accretion.

Gas giant planets play a fundamental role in shaping the orbital architecture of planetary systems and in affecting the delivery of volatile materials to terrestrial planets in the habitable zones. Current theories of gas giant planet formation rely on either of two mechanisms: the core accretion model and the disk instability model. In this chapter, we describe the essential principles upon which these models are built and discuss the successes and limitations of each model in explaining observational data of giant planets orbiting the Sun and other stars.

There has been disagreement currently about whether cooling in protoplanetary disks can be sufficiently fast to induce the formation of gas giant protoplanets via gravitational instabilities. Simulations by our own group and others indicate that this method of planet formation does not work for disks around young, low- mass stars inside several tens of AU, while simulations by other groups show fragmentation into protoplanetary clumps in this region. To allow direct comparison in hopes of isolating the cause of the differences, we here present a high resolution three-dimensional hydrodynamics simulation of a protoplanetary disk, where the disk model, initial perturbation, and simulation conditions are essentially identical to those used in a set of simulations by Boss. As in earlier papers by the same author, Boss (2007, hereafter B07) purports to show that cooling is fast enough to produce protoplanetary clumps. Here, we evolve the same B07 disk using an improved version of one of our own radiative schemes and find that the disk does not fragment in our code but instead quickly settles into a state with only low amplitude nonaxisymmetric structure, which persists for at least several outer disk rotations. We see no rapid radiative or convective cooling. We conclude that the differences in results are due to different treatments of regions at and above the disk photosphere, and we explain at least one way in which the scheme in B07 may lead to artificially fast cooling.

We test the effect of assumptions about stellar motion on the behavior of gravitational instabilities in protoplanetary disks around solar-type stars by performing two simulations that are identical in all respects except the treatment of the star. In one simulation, the star is assumed to remain fixed at the center of the inertial reference frame. In the other, stellar motion is handled properly by including an indirect potential in the hydrodynamic equations to model the star's reference frame as one which is accelerated by star/disk interactions. The disks in both simulations orbit a solar mass star, initially extend from 2.3 to 40 AU with a r^-1/2 surface density profile, and have a total mass of 0.14 M_sun. The gamma = 5/3 ideal gas is assumed to cool everywhere with a constant cooling time of two outer rotation periods. The overall behavior of the disk evolution is similar, except for weakening in various measures of GI activity by about at most tens of percent for the indirect potential case. Overall conclusions about disk evolution in earlier papers by our group, where the star was always assumed to be fixed in an inertial frame, remain valid. There is no evidence for independent one-armed instabilities, like SLING, in either simulation. On the other hand, the stellar motion about the system center of mass (COM) in the simulation with the indirect potential is substantial, up to 0.25 AU during the burst phase, as GIs initiate, and averaging about 0.9 AU during the asymptotic phase, when the GIs reach an overall balance of heating and cooling. These motions appear to be a stellar response to nonlinear interactions between discrete global spiral modes in both the burst and asymptotic phases of the evolution, and the star's orbital motion about the COM reflects the orbit periods of disk material near the corotation radii of the dominant spiral waves. This motion is, in principle, large enough to be observable and could be confused with stellar wobble due to the presence of one or more super-Jupiter mass protoplanets orbiting at 10's AU. We discuss why the excursions in our simulation are so much larger than those seen in simulations by Rice et al. 2003a.

We explore the initial conditions for fragments in the extended regions ($r\gtrsim 50$ AU) of gravitationally unstable disks. We combine analytic estimates for the fragmentation of spiral arms with 3D SPH simulations to show that initial fragment masses are in the gas giant regime. These initial fragments will have substantial angular momentum, and should form disks with radii of a few AU. We show that clumps will survive for multiple orbits before they undergo a second, rapid collapse due to H$_2$ dissociation and that it is possible to destroy bound clumps by transporting them into the inner disk. The consequences of disrupted clumps for planet formation, dust processing, and disk evolution are discussed. We argue that it is possible to produce Earth-mass cores in the outer disk during the earliest phases of disk evolution.

Using analytic arguments and numerical simulations, we examine whether chondrule formation and the FU Orionis phenomenon can be caused by the burst-like onset of gravitational instabilities (GIs) in dead zones. At least two scenarios for bursting dead zones can work, in principle. If the disk is on the verge of fragmention, GI activation near $r\sim4$ to 5 AU can produce chondrule-forming shocks, at least under extreme conditions. Mass fluxes are also high enough during the onset of GIs to suggest that the outburst is related to an FU Orionis phenomenon. This situation is demonstrated by numerical simulations. In contrast, as supported by analytic arguments, if the burst takes place close to $r\sim1$ AU, then even low pitch angle spiral waves can create chondrule-producing shocks and outbursts. We also study the stability of the massive disks in our simulations against fragmentation and find that although disk evolution is sensitive to changes in opacity, the disks we study do not fragment, even at high resolution and even for extreme assumptions.

It is generally thought that protoplanetary disks embedded in envelopes are more massive and thus more susceptible to gravitational instabilities (GIs) than exposed disks. We present three-dimensional radiative hydrodynamics simulations of protoplanetary disks with the presence of envelope irradiation. For a disk with a radius of 40 AU and a mass of 0.07 Msun around a young star of 0.5 Msun, envelope irradiation tends to weaken and even suppress GIs as the irradiating flux is increased. The global mass transport induced by GIs is dominated by lower-order modes, and irradiation preferentially suppresses higher-order modes. As a result, gravitational torques and mass inflow rates are actually increased by mild irradiation. None of the simulations produce dense clumps or rapid cooling by convection, arguing against direct formation of giant planets by disk instability, at least in irradiated disks. However, dense gas rings and radial mass concentrations are produced, and these might be conducive to accelerated planetary core formation. Preliminary results from a simulation of a massive embedded disk with physical characteristics similar to one of the disks in the embedded source L1551 IRS5 indicate a long radiative cooling time and no fragmentation. The GIs in this disk are dominated by global two and three-armed modes.

Numerical hydrodynamics simulations have established that disks which are evolved under the condition of local isothermality will fragment into small dense clumps due to gravitational instabilities when the Toomre stability parameter $Q$ is sufficiently low. Because fragmentation through disk instability has been suggested as a gas giant planet formation mechanism, it is important to understand the physics underlying this process as thoroughly as possible. In this paper, we offer analytic arguments for why, at low $Q$, fragments are most likely to form first at the corotation radii of growing spiral modes, and we support these arguments with results from 3D hydrodynamics simulations.

Recent three-dimensional radiative hydrodynamics simulations of protoplanetary disks report disparate disk behaviors, and these differences involve the importance of convection to disk cooling, the dependence of disk cooling on metallicity, and the stability of disks against fragmentation and clump formation. To guarantee trustworthy results, a radiative physics algorithm must demonstrate the capability to handle both the high and low optical depth regimes. We develop a test suite that can be used to demonstrate an algorithm's ability to relax to known analytic flux and temperature distributions, to follow a contracting slab, and to inhibit or permit convection appropriately. We then show that the radiative algorithm employed by Meji\'a (2004) and Boley et al. (2006) and the algorithm employed by Cai et al. (2006) and Cai et al. (2007, in prep.) pass these tests with reasonable accuracy. In addition, we discuss a new algorithm that couples flux-limited diffusion with vertical rays, we apply the test suite, and we discuss the results of evolving the Boley et al. (2006) disk with this new routine. Although the outcome is significantly different in detail with the new algorithm, we obtain the same qualitative answers. Our disk does not cool fast due to convection, and it is stable to fragmentation. We find an effective $\alpha\approx 10^{-2}$. In addition, transport is dominated by low-order modes.

The gas equation of state may be one of the critical factors for the disk instability theory of gas giant planet formation. This letter addresses the treatment of H$_2$ in hydrodynamical simulations of gravitationally unstable disks. In our discussion, we point out possible consequences of erroneous specific internal energy relations, approximate specific internal energy relations with discontinuities, and assumptions of constant $\Gamma_1$. In addition, we consider whether the ortho/para ratio for H$_2$ in protoplanetary disks should be treated dynamically as if the species are in equilibrium. Preliminary simulations indicate that the correct treatment is particularly critical for the study of gravitational instability when $T = 30$-50 K.

This paper presents a fully three-dimensional radiative hydrodymanics simulation with realistic opacities for a gravitationally unstable 0.07 Msun disk around a 0.5 Msun star. We address the following aspects of disk evolution: the strength of gravitational instabilities under realistic cooling, mass transport in the disk that arises from GIs, comparisons between the gravitational and Reynolds stresses measured in the disk and those expected in an alpha-disk, and comparisons between the SED derived for the disk and SEDs derived from observationally determined parameters. The mass transport in this disk is dominated by global modes, and the cooling times are too long to permit fragmentation for all radii. Moreover, our results suggest a plausible explanation for the FU Ori outburst phenomenon.

Observational studies show that the probability of finding gas giant planets around a star increases with the star's metallicity. Our latest simulations of disks undergoing gravitational instabilities (GIs) with realistic radiative cooling indicate that protoplanetary disks with lower metallicity generally cool faster and thus show stronger overall GI-activity. More importantly, the global cooling times in our simulations are too long for disk fragmentation to occur, and the disks do not fragment into dense protoplanetary clumps. Our results suggest that direct gas giant planet formation via disk instabilities is unlikely to be the mechanism that produced most observed planets. Nevertheless, GIs may still play an important role in a hybrid scenario, compatible with the observed metallicity trend, where structure created by GIs accelerates planet formation by core accretion.

In this paper, we describe the nonlinear outcome of spiral shocks in protoplanetary disks. Spiral shocks, for most protoplanetary disk conditions, create a loss of vertical force balance in the post-shock region and result in rapid expansion of the gas perpendicular to the disk midplane. This expansion has characteristics similar to hydraulic jumps, which occur in incompressible fluids. We present a theory to describe the behavior of these hybrids between shocks and hydraulic jumps (shock bores) and then compare the theory to three-dimensional hydrodynamics simulations. We discuss the fully three-dimensional shock structures that shock bores produce and discuss possible consequences for disk mixing, turbulence, and evolution of solids.

Protoplanetary gas disks are likely to experience gravitational instabilites (GI's) during some phase of their evolution. Density perturbations in an unstable disk grow on a dynamic time scale into spiral arms that produce efficient outward transfer of angular momentum and inward transfer of mass through gravitational torques. In a cool disk with rapid enough cooling, the spiral arms in an unstable disk form self-gravitating clumps. Whether gas giant protoplanets can form by such a disk instability process is the primary question addressed by this review. We discuss the wide range of calculations undertaken by ourselves and others using various numerical techniques, and we report preliminary results from a large multi-code collaboration. Additional topics include -- triggering mechanisms for GI's, disk heating and cooling, orbital survival of dense clumps, interactions of solids with GI-driven waves and shocks, and hybrid scenarios where GI's facilitate core accretion. The review ends with a discussion of how well disk instability and core accretion fare in meeting observational constraints.

We investigate the Mach number dependence of the radiative overstability of shocks. We find that the stability of radiative shocks increases with decreasing Mach number, with the result that M=2 shocks require cooling exponents < -1.2 to be overstable. The frequency of oscillation of the fundamental mode has a strong Mach number dependence. We find that feedback between the cooling region and the cold dense layer (CDL) further downstream is a function of Mach number, with stronger feedback and oscillation of the boundary between the CDL and the cooling region occuring at lower Mach numbers. This feedback can be quantified in terms of the reflection coefficient of sound waves. An interesting finding is that the stability properties of low Mach number shocks can be dramatically altered if the shocked gas is able to cool to temperatures less than the pre-shock value, and it is probable that low Mach number astrophysical shocks will be overstable in a variety of situations. The results of this work are relevant to astrophysical shocks with low Mach numbers, such as supernova remnants (SNRs) immersed in a hot interstellar medium (e.g., within a starburst region), and shocks in molecular clouds, where time-dependent chemistry can lead to overstability. (Abridged)

Chondrules are mm-sized spherules found throughout primitive, chondritic meteorites. Flash heating by a shock front is the leading explanation of their formation. However, identifying a mechanism for creating shock fronts inside the solar nebula has been difficult. In a gaseous disk capable of forming Jupiter, the disk must have been marginally gravitationally unstable at and beyond Jupiter's orbit. We show that this instability can drive inward spiral shock fronts with shock speeds of up to about 10 km/s at asteroidal orbits, sufficient to account for chondrule formation. Mixing and transport of solids in such a disk, combined with the planet-forming tendencies of gravitational instabilities, results in a unified scenario linking chondrite production with gas giant planet formation.

The core-accretion mechanism for gas giant formation may be too slow to create all observed gas giant planets during reasonable gas disk lifetimes, but it has yet to be firmly established that the disk instability model can produce permanent bound gaseous protoplanets under realistic conditions. Based on our recent simulations of gravitational instabilities in disks around young stars, we suggest that, even if instabilities due to disk self-gravity do not produce gaseous protoplanets directly, they may create persistent dense rings that are conducive to accelerated growth of gas giants through core accretion. The rings occur at and near the boundary between stable and unstable regions of the disk and appear to be produced by resonances with discrete spiral modes on the unstable side.

In order to investigate mass transport and planet formation by gravitational instabilities (GIs), we have extended our 3-D hydrodynamic simulations of protoplanetary disks from a previous paper. Our goal is to determine the asymptotic behavior of GIs and how it is affected by different constant cooling times. Initially, Rdisk = 40 AU, Mdisk = 0.07 Mo, M* = 0.5 Mo, and Qmin = 1.8. Sustained cooling, with tcool = 2 orps (outer rotation periods, 1 orp ~ 250 yrs), drives the disk to instability in ~ 4 orps. This calculation is followed for 23.5 orps. After 12 orps, the disk settles into a quasi-steady state with sustained nonlinear instabilities, an average Q = 1.44 over the outer disk, a well-defined power-law Sigma(r), and a roughly steady Mdot ~ 5(-7) Mo/yr. The transport is driven by global low-order spiral modes. We restart the calculation at 11.2 orps with tcool = 1 and 1/4 orp. The latter case is also run at high azimuthal resolution. We find that shorter cooling times lead to increased Mdots, denser and thinner spiral structures, and more violent dynamic behavior. The asymptotic total internal energy and the azimuthally averaged Q(r) are insensitive to tcool. Fragmentation occurs only in the high-resolution tcool = 1/4 orp case; however, none of the fragments survive for even a quarter of an orbit. Ring-like density enhancements appear and grow near the boundary between GI active and inactive regions. We discuss the possible implications of these rings for gas giant planet formation.

Previous simulations of self-gravitating protostellar disks have shown that, once developed, gravitational instabilities are enhanced by cooling the disk constantly during its evolution (Pickett et al. 2002). These earlier calculations included a very simple form of volumetric cooling, with a constant cooling time throughout the disk, which acted against the stabilizing effects of shock heating. The present work incorporates more realistic treatments of energy transport. The initial disk model extends from 2.3 to 40 AU, has a mass of 0.07 Msun and orbits a 0.5 Msun star. The models evolve for a period of over 2500 years, during which extensive spiral arms form. The disks structure is profoundly altered, transient clumps form in one case, but no permanent bound companion objects develop.

We present the results of an intensive coordinated monitoring campaign in the optical and X-ray wavelength ranges of the low-mass, pre-main sequence star V410Tau carried out with the aim to study the relation between various indicators for magnetic activity that probe emission from different atmospheric layers: optical photometric star spot (rotation) cycle, chromospheric Halpha emission, and coronal X-rays. Two X-ray pointings were carried out with the Chandra satellite simultaneously with the optical observations, and centered near the maximum and minimum levels of the optical lightcurve. A relation of their different count levels to the rotation period of the dominating spot is not confirmed by a third Chandra observation carried out some months later, during another minimum of the 1.87d cycle. Similarly we find no indications for a correlation of the Halpha emission with the spots' rotational phase. The extraordinary stability of the largest spot is confirmed by long-term photometric and radial velocity measurements. Joining our optical photometry with previous data we provide a new estimate for the dominant periodicity of V410Tau. This updated value removes systematic offsets of the time of minimum observed in data taken over the last decade. Furthermore, the combination of the new data with published measurements taken during the last decade allows us to examine long-term changes in the mean light level of the photometry of V410Tau. A variation on the timescale of 5.4yr is suggested. Assuming that this behavior is truely cyclic V410Tau is the first pre-main sequence star on which an activity cycle is detected.

In November 2001 we undertook a coordinated observing campaign to study the connection between X-ray and optical variability in the weak-line T Tauri star V410Tau. The observing plan included three 15 ksec observations with Chandra using the Advanced CCD Imaging Spectrometer for Spectroscopy scheduled for different phases of the known 1.87 d starspot cycle. Photometric and spectroscopic monitoring of V410Tau involving telescopes on three different continents was scheduled simultaneously with the Chandra exposures.