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.

We perform state-of-the-art, 3D, time-dependent simulations of magnetized disk winds, carried out to simulation scales of 60 Astronomical Units, in order to confront optical HST observations of protostellar jets. We ``observe'' the optical forbidden line emission produced by shocks within our simulated jets and compare these with actual observations. Our simulations reproduce the rich structure of time varying jets, including jet rotation far from the source, an inner (up to 400 km/s) and outer (less than 100 km/s) component of the jet, and jet widths of up to 20 Astronomical Units in agreement with observed jets. These simulations when compared with the data are able to constrain disk wind models. In particular, models featuring a disk magnetic field with a modest radial spatial variation across the disk are favored.

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.

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.

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.

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.