We develop a simplified model for studying the long-term evolution of giant planets in protoplanetary discs. The model accounts for the eccentricity evolution of the planets and the dynamics of eccentric discs under the influences of secular planet–disc interactions and internal disc pressure, self-gravity, and viscosity. Adopting the ansatz that the disc precesses coherently with aligned apsides, the eccentricity evolution equations of the planet–disc system reduce to a set of linearized ordinary differential equations, which allows for fast computation of the evolution of planet–disc eccentricities over long time-scales. Applying our model to ‘giant planet + external disc’ systems, we are able to reproduce and explain the secular behaviours found in previously published hydrodynamical simulations. We re-examine the possibility of eccentricity excitation (due to secular resonance) of multiple planets embedded in a dispersing disc, and find that taking into account the dynamics of eccentric discs can significantly affect the evolution of the planets’ eccentricities.
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.
more » « less- NSF-PAR ID:
- 10372429
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 516
- Issue:
- 4
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 5544-5554
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
ABSTRACT 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.
-
more » « less
ABSTRACT ϵ Eridani is the closest star to our Sun known to host a debris disc. Prior observations in the (sub-)millimetre regime have potentially detected clumpy structure in the disc and attributed this to interactions with an (as yet) undetected planet. However, the prior observations were unable to distinguish between structure in the disc and background confusion. Here, we present the first ALMA image of the entire disc, which has a resolution of 1.6 × 1.2 arcsec2. We clearly detect the star, the main belt, and two-point sources. The resolution and sensitivity of this data allow us to clearly distinguish background galaxies (that show up as point sources) from the disc emission. We show that the two-point sources are consistent with background galaxies. After taking account of these, we find that resolved residuals are still present in the main belt, including two clumps with a >3σ significance – one to the east of the star and the other to the north-west. We perform N-body simulations to demonstrate that a migrating planet can form structures similar to those observed by trapping planetesimals in resonances. We find that the observed features can be reproduced by a migrating planet trapping planetesimals in the 2:1 mean motion resonance and the symmetry of the most prominent clumps means that the planet should have a position angle of either ∼10° or ∼190°. Observations over multiple epochs are necessary to test whether the observed features rotate around the star.
-
more » « less
ABSTRACT Spirals in protoplanetary discs have been used to locate the potential planet in discs. Since only the spiral shape from a circularly orbiting perturber is known, most previous works assume that the planet is in a circular orbit. We develop a simple semi-analytical method to calculate the shape of the spirals launched by an eccentric planet. We assume that the planet emits wavelets during its orbit, and the wave fronts of these propagating wavelets form the spirals. The resulting spiral shape from this simple method agrees with numerical simulations exceptionally well. The spirals excited by an eccentric planet can detach from the planet, bifurcate, or even cross each other, which are all reproduced by this simple method. The spiral’s bifurcation point corresponds to the wavelet that is emitted when the planet’s radial speed reaches the disc’s sound speed. Multiple spirals can be excited by an eccentric planet (more than five spirals when e ≳ 0.2). The pitch angle and pattern speed are different between different spirals and can vary significantly across one spiral. The spiral wakes launched by high-mass eccentric planets steepen to spiral shocks and the crossing of spiral shocks leads to distorted or broken spirals. With the same mass, a more eccentric planet launches weaker spirals and induces a shallower gap over a long period of time. The observed unusually large/small pitch angles of some spirals, the irregular multiple spirals, and the different pattern speeds between different spirals may suggest the existence of eccentric perturbers in protoplanetary discs.
-
null (Ed.)ABSTRACT The maximum temperature and radial temperature profile in a protoplanetary disc are important for the condensation of different elements in the disc. We simulate the evolution of a set of protoplanetary discs from the collapse of their progenitor molecular cloud cores as well as the dust decoupling within the discs as they evolve. We show how the initial properties of the cloud cores affect the thermal history of the protoplanetary discs using a simple viscous disc model. Our results show that the maximum mid-plane temperature in the disc occurs within 0.5 au. It increases with the initial cloud temperature and decreases with its angular velocity and the viscosity of the disc. From the observed properties of the molecular cloud cores, we find the median value of the maximum temperature is around 1250 K, with roughly 90 per cent of them being less than 1500 K – a value that is lower than the 50 per cent condensation temperatures of most refractory elements. Therefore, only cloud cores with high initial temperatures or low-angular velocities and/or low viscosities within the planet-forming discs will result in refractory-rich planetesimals. To reproduce the volatile depletion pattern of CM, CO, and CV chondrites and the terrestrial planets in Solar system, one must either have rare properties of the initial molecular cloud cores like high core temperature, or other sources of energy to heat the disc to sufficiently high temperatures. Alternatively, the volatile depletion observed in these chondrites may be inherited from the progenitor molecular cloud.more » « less
