Search for: All records

A planetary system can undergo multiple episodes of intense dynamical activities throughout its life, resulting in the production of star-grazing planetesimals (or exocomets) and pollution of the host star. Such activity is especially pronounced when giant planets interact with other small bodies during the system’s evolution. However, due to the chaotic nature of the dynamics, it is difficult to determine the properties of the perturbing planet(s) from the observed planetesimal-disruption activities. In this study, we examine the outcomes of planetesimal-planet scatterings in a general setting. We focus on one-planet systems, and determine the likelihood and time-scale of planetesimal disruption by the host star as a function of the planet properties. We obtain a new analytical expression for the minimum distance a scattering body can reach, extending previous results by considering finite planet eccentricity and non-zero planetesimal mass. Through N-body simulations, we derive the distribution of minimum distances and the likelihood and time-scales of three possible outcomes of planetesimal-planet scatterings: collision with the planet, ejection, and disruption by the star. For planetesimals with negligible mass, we identify four defining dimensionless parameters (the planet eccentricity, planet-to-star mass ratio, planet radius to semimajor axis ratio, and the stellar disruption radius to planet semimajor axis ratio) that enable us to scale the problem and generalize our findings to a wide range of orbital configurations. Using these results, we explore three applications: falling evaporating bodies in the β Pictoris system, white dwarf pollution due to planetesimal disruption and planet engulfment by main-sequence stars.

Closely packed multiplanet systems are known to experience dynamical instability if the spacings between the planets are too small. Such instability can be tempered by the frictional forces acting on the planets from gaseous discs. A similar situation applies to stellar-mass black holes embedded in active galactic nuclei discs around supermassive black holes. We use N-body integrations to evaluate how the frictional damping of orbital eccentricity affects the growth of dynamical instability for a wide range of K (the difference in the planetary semimajor axes in units of the mutual Hill radius) and (unequal) planet masses. We find that, in general, the stable region (large K) and unstable region (small K) are separated by a “grey zone”, where the (in)stability is not guaranteed. We report the numerical values of the critical spacing for stability Kcrit and the “grey zone” range in different systems, and provide fitting formulae for arbitrary frictional forcing strength. We show that the stability of a system depends on the damping time-scale τ relative to the zero-friction instability growth time-scale tinst: two-planet systems are stable if tinst ≳ τ; three-planet systems require tinst ≳ 10τ−100τ. When K is sufficiently small, tinst can be less than the synodic period between the planets, which makes frictional stabilization unlikely to occur. As K increases, tinst tends to grow exponentially, but can also fluctuate by a few orders of magnitude. We also devise a linear map to analyse the dynamical instability of the “planet + test mass” system, and find qualitative agreement with N-body simulations.

In recent years, a number of eccentric debris belts have been observed in extrasolar systems. The most common explanation for their shape is the presence of a nearby eccentric planetary companion. The gravitational perturbation from such a companion would induce periodic eccentricity variations on the planetesimals in the belt, with a range of precession frequencies. The overall expected shape is an eccentric belt with a finite minimum width. However, several observed eccentric debris discs have been found to exhibit a narrower width than the theoretical expectation. In this paper, we study two mechanisms that can produce this small width: (i) the protoplanetary disc can interact with the planet and/or the planetesimals, slowly driving the eccentricity of the former and damping the eccentricities of the latter; and (ii) the companion planet could have gained its eccentricity stochastically, through planet–planet scatterings. We show that under appropriate conditions, both of these scenarios offer a plausible way to reduce the minimum width of an eccentric belt exterior to a perturbing planet. However, the effects of protoplanetary discs are diminished at large separations (a > 10 au) due to the scarcity of gas and the limited disc lifetime. These findings suggest that one can use the shape and width of debris discs to shed light on the evolution of extrasolar systems, constraining the protoplanetary disc properties and the prevalence of planet–planet scatterings. Further observations of debris-harbouring systems could confirm whether thin debris belts are a common occurrence, or the results of rare initial conditions or evolutionary processes.

Planets migrating in their natal discs can be captured into mean-motion resonance (MMR), in which the planets’ periods are related by integer ratios. Recent observations indicate that planets in MMR can be either apsidally aligned or anti-aligned. How these different configurations arise is unclear. In this paper, we study the MMR capture process of migrating planets, focusing on the property of the apsidal angles of the captured planets. We show that the standard picture of MMR capture, in which the planets undergo convergent migration and experience eccentricity damping due to planet–disc interactions, always leads to apsidal anti-alignment of the captured planets. However, when the planets experience eccentricity driving from the disc, apsidally aligned configuration in MMR can be produced. In this configuration, both planets’ resonance angles circulate, but a ‘mixed’ resonance angle librates and traps the planets near the nominal resonance location. The MMR capture process in the presence of disc eccentricity driving is generally complex and irregular, and can lead to various outcomes, including apsidal alignment and anti-alignment, as well as the disruption of the resonance. We suggest that the two resonant planets in the K2-19 system, with their moderate eccentricities and aligned apsides, have experienced eccentricity driving from their natal disc in the past.

We study the long-term evolution of two or more stellar black holes (BHs) on initially separated but unstable circular orbits around a supermassive BH (SMBH). Such a close-packed orbital configuration can naturally arise from BH migrations in the AGN disk. Dynamical instability of the orbits leads to recurring close encounters between two BHs, during which the BH separationr_pbecomes less than the Hill radiusR_H. In rare very close encounters, a tight merging BH binary can form with the help of gravitational wave emission. We useN-body simulations to study the time evolution of close encounters of various degrees ofcloseness. For a typical “SMBH+2BH” system, the averaged cumulative number of close encounters (withr_p≲R_H) scales approximately as ∝t^0.5. The minimum encounter separationr_pfollows a cumulative distributionP(<r_p) ∝r_pforr_p≪R_H. We obtain a semi-analytical expression for the averaged rate of binary captures that lead to BH mergers. Our results suggest that close-packed BHs in AGN disks may take a long time (≳10⁷orbits around the SMBH) to experience a sufficiently close encounter and form a bound binary. This time can be shorter if the initial BH orbits are highly aligned. The BH binary mergers produced in this scenario have high eccentricities when entering the LIGO band and broad distribution of orbital inclinations relative to the original AGN disk. We explore the effects of the gas disk and find that simple gas drags on the BHs do not necessarily lead to an enhanced BH binary capture rate.