Supermassive black hole binaries may form as a consequence of galaxy mergers. Both prograde and retrograde orbits have been proposed. We study a binary of a small mass ratio, q, in a retrograde orbit immersed in and interacting with a gaseous accreti
on disk in order to estimate time scales for inward migration leading to coalescence and the accretion rate to the secondary component. We employ both semi-analytic methods and two dimensional numerical simulations, focusing on the case where the binary mass ratio is small but large enough to significantly perturb the disk. We develop the theory of type I migration for this case and determine conditions for gap formation finding that then inward migration occurs on a time scale equal to the time required for one half of the secondary mass to be accreted through the unperturbed disk, with accretion onto the secondary playing only a minor role. The semi-analytic and fully numerical approaches are in good agreement, the former being applicable over long time scales. Inward migration induced by interaction with the disk alleviates the final parsec problem. Accretion onto the secondary does not significantly affect the orbital evolution, but may have observational consequences for high accretion efficiency. The binary may then appear as two sources of radiation rotating around each other. This study should be extended to consider orbits with significant eccentricity and the effects of gravitational radiation at small length scales. Note too that torques acting between a circumbinary disk and a retrograde binary orbit may cause the mutual inclination to increase on a timescale that can be similar to, or smaller than that for orbital evolution, depending on detailed parameters. This is also an aspect for future study (abridged).
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 wo
rk 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.
We consider the radial migration of vortices in two-dimensional isothermal gaseous disks. We find that a vortex core, orbiting at the local gas velocity, induces velocity perturbations that propagate away from the vortex as density waves. The resulti
ng spiral wave pattern is reminiscent of an embedded planet. There are two main causes for asymmetries in these wakes: geometrical effects tend to favor the outer wave, while a radial vortensity gradient leads to an asymmetric vortex core, which favors the wave at the side that has the lowest density. In the case of asymmetric waves, which we always find except for a disk of constant pressure, there is a net exchange of angular momentum between the vortex and the surrounding disk, which leads to orbital migration of the vortex. Numerical hydrodynamical simulations show that this migration can be very rapid, on a time scale of a few thousand orbits, for vortices with a size comparable to the scale height of the disk. We discuss the possible effects of vortex migration on planet formation scenarios.
We study the effects of diffusion on the non-linear corotation torque, or horseshoe drag, in the two-dimensional limit, focusing on low-mass planets for which the width of the horseshoe region is much smaller than the scale height of the disc. In the
absence of diffusion, the non-linear corotation torque saturates, leaving only the Lindblad torque. Diffusion of heat and momentum can act to sustain the corotation torque. In the limit of very strong diffusion, the linear corotation torque is recovered. For the case of thermal diffusion, this limit corresponds to having a locally isothermal equation of state. We present some simple models that are able to capture the dependence of the torque on diffusive processes to within 20% of the numerical simulations.
Outward migration of low-mass planets has recently been shown to be a possibility in non-barotropic disks. We examine the consequences of this result in evolutionary models of protoplanetary disks. Planet migration occurs towards equilibrium radii wi
th zero torque. These radii themselves migrate inwards because of viscous accretion and photoevaporation. We show that as the surface density and temperature fall, the planet orbital migration and disk depletion timescales eventually become comparable, with the precise timing depending on the mass of the planet. When this occurs, the planet decouples from the equilibrium radius. At this time, however, the gas surface density is already too low to drive substantial further migration. A higher mass planet, of 10 Earth masses, can open a gap during the late evolution of the disk, and stops migrating. Low mass planets, with 1 or 0.1 Earth masses, released beyond 1 AU in our models, avoid migrating into the star. Our results provide support for the reduced migration rates adopted in recent planet population synthesis models.
We study the collisional evolution of km-sized planetesimals in tight binary star systems to investigate whether accretion towards protoplanets can proceed despite the strong gravitational perturbations from the secondary star. The orbits of planetes
imals are numerically integrated in two dimensions under the influence of the two stars and gas drag. The masses and orbits of the planetesimals are allowed to evolve due to collisions with other planetesimals and accretion of collisional debris. In addition, the mass in debris can evolve due to planetesimal-planetesimal collisions and the creation of new planetesimals. We show that it is possible in principle for km-sized planetesimals to grow by two orders of magnitude in size if the efficiency of planetesimal formation is relatively low. We discuss the limitations of our two-dimensional approach.
One class of protoplanetary disc models, the X-wind model, predicts strongly subkeplerian orbital gas velocities, a configuration that can be sustained by magnetic tension. We investigate disc-planet interactions in these subkeplerian discs, focusing
on orbital migration for low-mass planets and gap formation for high-mass planets. We use linear calculations and nonlinear hydrodynamical simulations to measure the torque and look at gap formation. In both cases, the subkeplerian nature of the disc is treated as a fixed external constraint. We show that, depending on the degree to which the disc is subkeplerian, the torque on low-mass planets varies between the usual Type I torque and the one-sided outer Lindblad torque, which is also negative but an order of magnitude stronger. In strongly subkeplerian discs, corotation effects can be ignored, making migration fast and inward. Gap formation near the planets orbit is more difficult in such discs, since there are no resonances close to the planet accommodating angular momentum transport. In stead, the location of the gap is shifted inwards with respect to the planet, leaving the planet on the outside of a surface density depression. Depending on the degree to which a protoplanetary disc is subkeplerian, disc-planet interactions can be very different from the usual Keplerian picture, making these discs in general more hazardous for young planets.
We study the disc planet interactions of low-mass protoplanets embedded in a circumstellar disc. We extend the standard theory of planet migration from the usual locally isothermal assumption to include non-barotropic effects, focusing on the validit
y of linear theory. We compared solutions of the linear equations with results from non-linear hydrodynamic simulations, where in both cases we adopted a background entropy gradient and solved the energy equation. We show that the migration behavior of embedded planets depends critically on the background radial entropy gradient in the disc. The presence of such a gradient not only changes the corotation torque on the planet, but also always guarantees a departure from linear behavior, which gives a singular density response at corotation, in the absence of thermal or viscous diffusion. A negative entropy gradient tends to give rise to positive, non-linear corotation torques apparently produced as material executes horseshoe turns at approximately constant entropy. These torques have no counterpart in linear theory, but can be strong enough to push the planet outwards until saturation starts to occur after a horseshoe libration period. Increased thermal diffusion acts to reduce these non-linear torques, but, at the same time, it can help to prevent their saturation. In combination with a small kinematic viscosity that is able to maintain a smooth density profile the positive torque could be sustained.
Observations of extrasolar planets reveal that planets can be found in close binary systems, where the semi-major axis of the binary orbit is less than 20 AU. The existence of these planets challenges planet formation theory, because the strong gravi
tational perturbations due to the companion increase encounter velocities between planetesimals and make it difficult for them to grow through accreting collisions. We study planetesimal encounter velocities in binary systems, where the planetesimals are embedded in a circumprimary gas disc that is allowed to evolve under influence of the gravitational perturbations of the companion star. We find that the encounter velocities between planetesimals of different size strongly depend on the gas disc eccentricity. In all cases studied, inclusion of the full gas dynamics increases the encounter velocity compared to the case of a static, circular gas disc. Full numerical parameter exploration is still impossible, but we derive analytical formulae to estimate encounter velocities between bodies of different sizes given the gas disc eccentricity. The gas dynamical evolution of a protoplanetary disc in a binary system tends to make planetesimal accretion even more difficult than in a static, axisymmetric gas disc.
We study the interaction of a low-mass planet with a protoplanetary disk with a realistic treatment of the energy balance by doing radiation-hydrodynamical simulations. We look at accretion and migration rates and compare them to isothermal studies.
We used a three-dimensional version of the hydrodynamical method RODEO, together with radiative transport in the flux-limited diffusion approach. The accretion rate, as well as the torque on the planet, depend critically on the ability of the disk to cool efficiently. For densities appropriate to 5 AU in the solar nebula, the accretion rate drops by more than an order of magnitude compared to isothermal models, while at the same time the torque on the planet is positive, indicating outward migration. It is necessary to lower the density by a factor of 2 to recover inward migration and more than 2 orders of magnitude to recover the usual Type I migration. The torque appears to be proportional to the radial entropy gradient in the unperturbed disk. These findings are critical for the survival of protoplanets, and they should ultimately find their way into population synthesis models.
