In simulations of viscously evolving accretion disks, the inner boundary condition is particularly important. If treated incorrectly, it induces incorrect behavior very quickly, because the viscous time is shortest near the inner boundary. Recent work has determined the correct inner boundary in Eulerian simulations. But in quasi-Lagrangian simulations (e.g., SPH, moving mesh, and mesh-less), where the inner boundary is modeled by removing mass within a finite zone, the inner density profile typically becomes anomalously depleted. Here we show how the boundary condition should be applied in such codes, via a simple modification of the usual approach: when one removes mass, one must speed up the remaining material so that the disk's angular momentum is unchanged. We show with both 1D and 2D moving-mesh (AREPO) simulations that this scheme works as desired in viscously evolving disks. It produces no spurious density depletions and is independent of the mass removal rate, provided that the disk is adequately resolved and that the mass removal rate is not so extreme as to trigger instabilities. This "torque-free" mass removal technique permits the use of quasi-Lagrangian codes to simulate viscously evolving disks, while including a variety of additional effects. As an example, we apply our scheme to a 2D simulation of an accretion disk perturbed by a very massive planet, in which the disk is evolved to viscous steady state.

We investigate how planets interact with viscous accretion disks, in the limit that the disk is sufficiently low mass that the planet migrates more slowly than the disk material. In that case, the disk's surface density profile is determined by the disk being in viscous steady state (VSS), while overflowing the planet's orbit. We compute the VSS profiles with 2D hydrodynamical simulations, and show that disk material piles up behind the planet, with the planet effectively acting as a leaky dam. Previous 2D hydrodynamical simulations missed the pileup effect because of incorrect boundary conditions, while previous 1D models greatly overpredicted the pileup due to the neglect of non-local deposition. Our simulations quantify the magnitude of the pileup for a variety of planet masses and disk viscosities. We also calculate theoretically the magnitude of the pileup for moderately deep gaps, showing good agreement with simulations. For very deep gaps, current theory is inadequate, and we show why and what must be understood better. The pileup is important for two reasons. First, it is observable in directly imaged protoplanetary disks, and hence can be used to diagnose the mass of a planet that causes it or the viscosity within the disk. And second, it determines the planet's migration rate. Our simulations determine a new Type-II migration rate (valid for low mass disks), and show how it connects continuously with the well-verified Type-I rate.

Accretion disks can be eccentric: they support $m=1$ modes that are global and slowly precessing. But whether the modes remain trapped in the disk---and hence are long-lived---depends on conditions at the outer edge of the disk. Here we show that in disks with realistic boundaries, in which the surface density drops rapidly beyond a given radius, eccentric modes are trapped and hence long-lived. We focus on pressure-only disks around a central mass, and show how this result can be understood with the help of a simple second-order WKB theory. We show that the longest lived mode is the zero-node mode in which all of the disk's elliptical streamlines are aligned, and that this mode decays coherently on the viscous timescale of the disk. Hence such a mode, once excited, will live for the lifetime of the disk. It may be responsible for asymmetries seen in recent images of protoplanetary disks.

Accretion disks around stars, or other central massive bodies, can support long-lived, slowly precessing $m=1$ disturbances in which the fluid motion is nearly Keplerian with non-zero eccentricity. We study such `slow modes' in disks that are subject to both pressure and self-gravity forces. We derive a second-order WKB dispersion relation that describes the dynamics quite accurately, and show that the apparently complicated nature of the various modes can be understood in a simple way with the help of a graphical method. We also solve the linearized fluid equations numerically, and show that the results agree with the theory. We find that when self-gravity is weak ($Q\gtrsim 1/h$, where $Q$ is Toomre's parameter, and $h$ is the disk aspect ratio) the modes are pressure dominated. But when self-gravity is strong ($1<Q\lesssim 1/h$), two kinds of gravity-dominated modes appear: one is an aligned elliptical pattern and the other is a one-armed spiral. In the context of protoplanetary disks, we suggest that if the radial eccentricity profile can be measured, it could be used to determine the total disk mass.

We run a suite of dissipative N-body simulations to determine which regions of phase space for smooth disk migration are consistent with the GJ876 system, an M-dwarf hosting three planets orbiting in a chaotic 4:2:1 Laplace resonance. We adopt adaptive mesh refinement (AMR) methods which are commonly used in hydrodynamical simulations to efficiently explore the parameter space defined by the semi-major axis and eccentricity damping timescales. We find that there is a large region of phase space which produces systems in the chaotic Laplace resonance and a smaller region consistent with the observed eccentricities and libration amplitudes for the resonant angles. Under the assumptions of Type I migration for the outer planet, we translate these damping timescales into constraints on the protoplanetary disk surface density and thickness. When we strongly (weakly) damp the eccentricities of the inner two Laplace planets, these timescales correspond to disk surface densities around ten thousand (a few hundred) grams per square centimeter and disk aspect ratios between 1-10%. Additionally, smooth migration produces systems with a range of chaotic timescales, from decades and centuries to upwards of thousands of years. In agreement with previous studies, the less chaotic regions of phase space coincide with the system being in a low energy double apsidal corotation resonance. Our detailed modeling of multi-planetary systems coupled with our AMR exploration method enhances our ability to map out the parameter space of planet formation models, and is well suited to study other resonant chain systems such as Trappist-1, Kepler-60, and others.

We study thermal convection in a rotating fluid in order to better understand the properties of convection zones in rotating stars and planets. We first derive mixing-length theory for rapidly-rotating convection, arriving at the results of Stevenson (1979) via simple physical arguments. The theory predicts the properties of convection as a function of the imposed heat flux and rotation rate, independent of microscopic diffusivities. In particular, it predicts the mean temperature gradient; the rms velocity and temperature fluctuations; and the size of the eddies that dominate heat transport. We test all of these predictions with high resolution three-dimensional hydrodynamical simulations of Boussinesq convection in a Cartesian box. The results agree remarkably well with the theory across more than two orders of magnitude in rotation rate. For example, the temperature gradient is predicted to scale as the rotation rate to the 4/5th power at fixed flux, and the simulations yield $0.75\pm 0.06$. We conclude that the mixing length theory is a solid foundation for understanding the properties of convection zones in rotating stars and planets.