We study the dynamical evolution of Jupiter and Saturn embedded in a gaseous, solar-nebula-type disc by means of hydrodynamics simulations with the FARGO2D1D code. We study the evolution for different initial separations of the planets' orbits, Delta a_SJ , to investigate whether they become captured in mean motion resonance (MMR) and the direction of the subsequent migration of the planet (inward or outward). We also provide an assessment of the planet's orbital dynamics at different epochs of Saturn's growth. We find that the evolution of initially compact orbital configurations is dependent on the value of Delta a_SJ . This implies that an evolution as that proposed in the Grand Tack model depends on the precise initial orbits of Jupiter and Saturn and on the timescales for their formation. Capture in the 1:2 MMR and inward or (nearly) stalled migration are highly favoured. Within its limits, our work suggests that the reversed migration, associated with the resonance capture of Jupiter and Saturn, may be a low probability evolutionary scenario, so that other planetary systems with giant planets are not expected to have experienced a Grand Tack-like evolutionary path.

An extended gravitational object embedded in a parent system comprised of gas and collisionless particles may undergo both dynamical friction (DF) and mass loss by tidal forces. If the object is compact enough, it can increase its mass through accretion of material from the surrounding medium. We extend the classical linear analysis of DF on a constant-mass body in a gaseous medium to the case where its mass changes with time. We show that the structure of the wake may differ significantly from the constant-mass case. For instance, the front-back symmetry of density about subsonic constant-mass perturbers is broken down for variable-mass perturbers. The density wake keeps a memory of the past mass history of the perturber. For dissolving perturbers, the density field is more dense than expected using the instantaneous mass of the perturber in the classical formula. As a consequence, the instantaneous-mass approximation underestimates the drag force for mass-losing perturbers and overestimates it for mass-gaining perturbers. We present cases in which the percentage error in the drag force using the instantaneous-mass approximation is greater than 50%.

We evaluate the torque acting on a gravitational perturber on a retrograde circular orbit in the midplane of a gaseous disk. We assume that the mass of this satellite is so low it weakly disturbs the disk (type I migration). The perturber may represent the companion of a binary system with a small mass ratio. We compare the results of hydrodynamical simulations with analytic predictions. Our two-dimensional (2D) simulations indicate that the torque acting on a perturber with softening radius $R_{\rm soft}$ can be accounted for by a scattering approach if $R_{\rm soft}<0.3H$, where $H$ is defined as the ratio between the sound speed and the angular velocity at the orbital radius of the perturber. For $R_{\rm soft}>0.3H$, the torque may present large and persistent oscillations, but the resultant time-averaged torque decreases rapidly with increasing $R_{\rm soft}/H$, in agreement with previous analytical studies. We then focus on the torque acting on small-size perturbers embedded in full three-dimensional (3D) disks and argue that the density waves propagating at distances $\lesssim H$ from the perturber contribute significantly to the torque because they transport angular momentum. We find a good agreement between the torque found in 3D simulations and analytical estimates based on ballistic orbits. We compare the radial migration timescales of prograde versus retrograde perturbers. For a certain range of the perturber's mass and aspect ratio of the disk, the radial migration timescale in the retrograde case may be appreciably shorter than in the prograde case. We also provide the smoothing length required in 2D simulations in order to account for 3D effects.

We study gap formation in gaseous protoplanetary discs by a Jupiter mass planet. The planet's orbit is circular and inclined relative to the midplane of the disc. We use the impulse approximation to estimate the gravitational tidal torque between the planet and the disc, and infer the gap profile. For low-mass discs, we provide a criterion for gap opening when the orbital inclination is $\leq 30^{\circ}$. Using the FARGO3D code, we simulate the disc response to an inclined massive planet. The dependence of the depth and width of the gap obtained in the simulations on the inclination of the planet is broadly consistent with the scaling laws derived in the impulse approximation. Although we mainly focus on planets kept on fixed orbits, the formalism permits to infer the temporal evolution of the gap profile in cases where the inclination of the planet changes with time. This study may be useful to understand the migration of massive planets on inclined orbit, because the strength of the interaction with the disc depends on whether a gap is opened or not.

Using time-dependent linear theory, we investigate the morphology of the gravitational wake induced by a binary, whose center of mass moves at velocity Vcm against a uniform background of gas. For simplicity, we assume that the binary's components are on circular orbits about their common center of mass. The consequences of dynamical friction is twofold. First, gas dynamical friction may drag the binary's center of mass and cause the binary to migrate. Second, drag forces also induce a braking torque, which causes the orbits of the binary components to shrink. We compute the drag forces acting on one component of the binary due to the gravitational interaction with its own wake. We show that the dynamical friction force responsible to decelerate the binary's center of mass is smaller than it is in the point-mass case because of the loss of gravitational focusing. We show that the braking internal torque depends on the Mach numbers of each binary component about their center of mass, and also on the Mach number of the center of mass of the binary. In general, the internal torque decreases with increasing the velocity of the binary relative to the ambient gas cloud. However, this is not always the case. We also mention the relevance of our results on the period distribution of binaries.