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.

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 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.

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.

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.