No Arabic abstract
We describe 2D hydrodynamic simulations of the migration of low-mass planets ($leq 30 M_{oplus}$) in nearly laminar disks (viscosity parameter $alpha < 10^{-3}$) over timescales of several thousand orbit periods. We consider disk masses of 1, 2, and 5 times the minimum mass solar nebula, disk thickness parameters of $H/r = 0.035$ and 0.05, and a variety of $alpha$ values and planet masses. Disk self-gravity is fully included. Previous analytic work has suggested that Type I planet migration can be halted in disks of sufficiently low turbulent viscosity, for $alpha sim 10^{-4}$. The halting is due to a feedback effect of breaking density waves that results in a slight mass redistribution and consequently an increased outward torque contribution. The simulations confirm the existence of a critical mass ($M_{cr} sim 10 M_{oplus}$) beyond which migration halts in nearly laminar disks. For $alpha ga 10^{-3}$, density feedback effects are washed out and Type I migration persists. The critical masses are in good agreement with the analytic model of Rafikov (2002). In addition, for $alpha la 10^{-4}$ steep density gradients produce a vortex instability, resulting in a small time-varying eccentricity in the planets orbit and a slight outward migration. Migration in nearly laminar disks may be sufficiently slow to reconcile the timescales of migration theory with those of giant planet formation in the core accretion model.
Numerical simulations of planet-disk interactions are usually performed with hydro-codes that -- because they consider only an annulus of the disk, over a 2D grid -- can not take into account the global evolution of the disk. However, the latter governs planetary migration of type II, so that the accuracy of the planetary evolution can be questioned. To develop an algorithm that models the local planet-disk interactions together with the global viscous evolution of the disk, we surround the usual 2D grid with a 1D grid ranging over the real extension of the disk. The 1D and 2D grids are coupled at their common boundaries via ghost rings, paying particular attention to the fluxes at the interface, especially the flux of angular momentum carried by waves. The computation is done in the frame centered on the center of mass to ensure angular momentum conservation. The global evolution of the disk and the local planet-disk interactions are both well described and the feedback of one on the other can be studied with this algorithm, for a negligible additional computing cost with respect to usual algorithms.
We investigate the unsaturated horseshoe drag exerted on a low-mass planet by an isothermal gaseous disk. In the globally isothermal case, we use a formal- ism, based on the use of a Bernoulli invariant, that takes into account pressure effects, and that extends the torque estimate to a region wider than the horse- shoe region. We find a result that is strictly identical to the standard horseshoe drag. This shows that the horseshoe drag accounts for the torque of the whole corotation region, and not only of the horseshoe region, thereby deserving to be called corotation torque. We find that evanescent waves launched downstream of the horseshoe U-turns by the perturbations of vortensity exert a feed-back on the upstream region, that render the horseshoe region asymmetric. This asymmetry scales with the vortensity gradient and with the disks aspect ratio. It does not depend on the planetary mass, and it does not have any impact on the horseshoe drag. Since the horseshoe drag has a steep dependence on the width of the horseshoe region, we provide an adequate definition of the width that needs to be used in horseshoe drag estimates. We then consider the case of locally isothermal disks, in which the tempera- ture is constant in time but depends on the distance to the star. The horseshoe drag appears to be different from the case of a globally isothermal disk. The difference, which is due to the driving of vortensity in the vicinity of the planet, is intimately linked to the topology of the flow. We provide a descriptive inter- pretation of these effects, as well as a crude estimate of the dependency of the excess on the temperature gradient.
Planet migration originally refers to protoplanetary disks, which are more massive and dense than typical accretion disks in binary systems. We study planet migration in an accretion disk in a binary system consisting of a solar-like star hosting a planet and a red giant donor star. The accretion disk is fed by a stellar wind. %, disk self-gravity is neglected. We use the $alpha$-disk model and consider that the stellar wind is time-dependent. Assuming the disk is quasi-stationary we calculate its temperature and surface density profiles. In addition to the standard disk model, when matter is captured by the disk at its outer edge, we study the situation when the stellar wind delivers matter on the whole disc surface inside the accretion radius with the rate depending on distance from the central star. Implying that a planet experiences classical type I/II migration we calculate migration time for a planet on a circular orbit coplanar with the disk. Potentially, rapid inward planet migration can result in a planet-star merger which can be accompanied by an optical or/and UV/X-ray transient. We calculate timescale of migration for different parameters of planets and binaries. Our results demonstrate that planets can fall on their host stars within the lifetime of the late-type donor for realistic sets of parameters.
Earth-mass planets embedded in gaseous protoplanetary disks undergo Type I orbital migration. In radiative disks an additional component of the corotation torque scaling with the entropy gradient across the horseshoe region can counteract the general inward migration, Type I migration can then be directed either inward or outward depending on the local disk properties. Thus, special locations exist in the disk toward which planets migrate in a convergent way. Here we present N-body simulations of the convergent migration of systems of low-mass (M=1-10 m_earth) planets. We show that planets do not actually converge in convergence zones. Rather, they become trapped in chains of mean motion resonances (MMRs). This causes the planets eccentricities to increase to high enough values to affect the structure of the horseshoe region and weaken the positive corotation torque. The zero-torque equilibrium point of the resonant chain of planets is determined by the sum of the attenuated corotation torques and unattenuated differential Lindblad torques acting on each planet. The effective convergence zone is shifted inward. Systems with several planets can experience stochastic migration as a whole due to continuous perturbations from planets entering and leaving resonances.
The increasing number of newly detected exoplanets at short orbital periods raises questions about their formation and migration histories. A particular puzzle that requires explanation arises from one of the key results of the Kepler mission, namely the increase in the planetary occurrence rate with orbital period up to 10 days for F, G, K and M stars. We investigate the conditions for planet formation and migration near the dust sublimation front in protostellar disks around young Sun-like stars. For this analysis we use iterative 2D radiation hydrostatic disk models which include irradiation by the star, and dust sublimation and deposition depending on the local temperature and vapor pressure. We perform a parameter study by varying the magnetized turbulence onset temperature, the accretion stress, the dust mass fraction, and the mass accretion rate. Our models feature a gas-only inner disk, a silicate sublimation front and dust rim starting at around 0.08 au, an ionization transition zone with a corresponding density jump, and a pressure maximum which acts as a pebble trap at around 0.12 au. Migration torque maps show Earth- and super-Earth-mass planets halt in our model disks at orbital periods ranging from 10 to 22 days. Such periods are in good agreement with both the inferred location of the innermost planets in multiplanetary systems, and the break in planet occurrence rates from the Kepler sample at 10 days. In particular, models with small grains depleted produce a trap located at a 10-day orbital period, while models with a higher abundance of small grains present a trap at around a 17-day orbital period. The snow line lies at 1.6 au, near where the occurrence rate of the giant planets peaks. We conclude that the dust sublimation zone is crucial for forming close-in planets, especially when considering tightly packed super-Earth systems.