An official website of the United States government
ABSTRACT We study the evolution of eccentricity and inclination of massive planets in low-density cavities of protoplanetary discs using three-dimensional (3D) simulations. When the planet’s orbit is aligned with the equatorial plane of the disc, the eccentricity increases to high values of 0.7–0.9 due to the resonant interaction with the inner parts of the disc. For planets on inclined orbits, the eccentricity increases due to the Kozai–Lidov mechanism, where the disc acts as an external massive body, which perturbs the planet’s orbit. At small inclination angles, $${\lesssim}30^\circ$$, the resonant interaction with the inner disc strongly contributes to the eccentricity growth, while at larger angles, eccentricity growth is mainly due to the Kozai–Lidov mechanism. We conclude that planets inside low-density cavities tend to acquire high eccentricity if favourable conditions give sufficient time for growth. The final value of the planet’s eccentricity after the disc dispersal depends on the planet’s mass and the properties of the cavity and protoplanetary disc.
more »
« less
Eccentricity growth of massive planets inside cavities of protoplanetary discs
ABSTRACT We carry out hydrodynamical simulations to study the eccentricity growth of a 1–30 Jupiter mass planet located inside the fixed cavity of a protoplanetary disc. The planet exchanges energy and angular momentum with the disc at resonant locations, and its eccentricity grows due to Lindblad resonances. We observe several phases of eccentricity growth where different eccentric Lindblad resonances dominate from 1:3 up to 3:5. The maximum values of eccentricity reached in our simulations are 0.65–0.75. We calculate the eccentricity growth rate for different planet masses and disc parameters and derive analytical dependencies on these parameters. We observe that the growth rate is proportional to both the planet’s mass and the characteristic disc mass for a wide range of parameters. In a separate set of simulations, we derived the width of the 1:3 Lindblad resonance.
more »
« less
- Award ID(s):
- 2009820
- PAR ID:
- 10421342
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 523
- Issue:
- 2
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 2832-2849
- Size(s):
- p. 2832-2849
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT We developed convolutional neural networks (CNNs) to rapidly and directly infer the planet mass from radio dust continuum images. Substructures induced by young planets in protoplanetary discs can be used to infer the potential young planets’ properties. Hydrodynamical simulations have been used to study the relationships between the planet’s properties and these disc features. However, these attempts either fine-tuned numerical simulations to fit one protoplanetary disc at a time, which was time consuming, or azimuthally averaged simulation results to derive some linear relationships between the gap width/depth and the planet mass, which lost information on asymmetric features in discs. To cope with these disadvantages, we developed Planet Gap neural Networks (PGNets) to infer the planet mass from two-dimensional images. We first fit the gridded data in Zhang et al. as a classification problem. Then, we quadrupled the data set by running additional simulations with near-randomly sampled parameters, and derived the planet mass and disc viscosity together as a regression problem. The classification approach can reach an accuracy of 92 per cent, whereas the regression approach can reach 1σ as 0.16 dex for planet mass and 0.23 dex for disc viscosity. We can reproduce the degeneracy scaling α ∝ $$M_\mathrm{ p}^3$$ found in the linear fitting method, which means that the CNN method can even be used to find degeneracy relationship. The gradient-weighted class activation mapping effectively confirms that PGNets use proper disc features to constrain the planet mass. We provide programs for PGNets and the traditional fitting method from Zhang et al., and discuss each method’s advantages and disadvantages.more » « less
-
ABSTRACT The hunt is on for dozens of protoplanets hypothesized to reside in protoplanetary discs with imaged gaps. How bright these planets are, and what they will grow to become, depend on their accretion rates, which may be in the runaway regime. Using 3D global simulations, we calculate maximum gas accretion rates for planet masses Mp from 1$$\, \mathrm{ M}_{{\oplus }}$$ to $$10\, \mathrm{ M}_{\rm J}$$. When the planet is small enough that its sphere of influence is fully embedded in the disc, with a Bondi radius rBondi smaller than the disc’s scale height Hp – such planets have thermal mass parameters qth ≡ (Mp/M⋆)/(Hp/Rp)3 ≲ 0.3, for host stellar mass M⋆ and orbital radius Rp – the maximum accretion rate follows a Bondi scaling, with $$\max \dot{M}_{\rm p} \propto \rho _{\rm g}M_{\rm p}^2 / (H_{\rm p}/R_{\rm p})^3$$ for ambient disc density ρg. For more massive planets with 0.3 ≲ qth ≲ 10, the Hill sphere replaces the Bondi sphere as the gravitational sphere of influence, and $$\max \dot{M}_{\rm p} \propto \rho _{\rm g}M_{\rm p}^1$$, with no dependence on Hp/Rp. In the strongly superthermal limit when qth ≳ 10, the Hill sphere pops well out of the disc, and $$\max \dot{M}_{\rm p} \propto \rho _{\rm g}M_{\rm p}^{2/3} (H_{\rm p}/R_{\rm p})^1$$. Applied to the two confirmed protoplanets PDS 70b and c, our numerically calibrated maximum accretion rates imply that their Jupiter-like masses may increase by up to a factor of ∼2 before their parent disc dissipates.more » « less
-
ABSTRACT We use the general relativistic radiation magnetohydrodynamics code KORAL to simulate the accretion disc formation resulting from the tidal disruption of a solar mass star around a supermassive black hole (BH) of mass 106 M⊙. We simulate the disruption of artificially more bound stars with orbital eccentricity e ≤ 0.99 (compared to the more realistic case of parabolic orbits with e = 1) on close orbits with impact parameter β ≥ 3. We use a novel method of injecting the tidal stream into the domain, and we begin the stream injection at the peak fallback rate in this study. For two simulations, we choose e = 0.99 and inject mass at a rate that is similar to parabolic TDEs. We find that the disc only becomes mildly circularized with eccentricity e ≈ 0.6 within the 3.5 d that we simulate. The rate of circularization is faster for pericenter radii that come closer to the BH. The emitted radiation is mildly super-Eddington with $$L_{\rm {bol}}\approx 3{-}5\, L_{\rm {Edd}}$$ and the photosphere is highly asymmetric with the photosphere being significantly closer to the inner accretion disc for viewing angles near pericenter. We find that soft X-ray radiation with Trad ≈ 3–5 × 105 K may be visible for chance viewing angles. Our simulations suggest that TDEs should be radiatively inefficient with η ≈ 0.009–0.014.more » « less
-
ABSTRACT Mean motion resonances are important in the analysis and understanding of the dynamics of planetary systems. While perturbative approaches have been dominant in many previous studies, recent non-perturbative approaches have revealed novel properties in the low-eccentricity regime for interior mean motion resonances of Jupiter in the fundamental model of the circular planar restricted three-body model. Here, we extend the non-perturbative investigation to exterior mean motion resonances in the low-eccentricity regime (up to about 0.1) and for perturber mass in the range of ∼5 × 10−5 to 1 × 10−3 (in units of the central mass). Our results demonstrate that first-order exterior resonances have two branches at low eccentricity as well as low-eccentricity bridges connecting neighbouring first-order resonances. With increasing perturber mass, higher order resonances dissolve into chaos, whereas low-order resonances persist with larger widths in their radial extent but smaller azimuthal widths. For low-order resonances, we also detect secondary resonances arising from small-integer commensurabilities between resonant librations and the synodic frequency. These secondary resonances contribute significantly to generating the chaotic sea that typically occurs near mean motion resonances of higher mass perturbers.more » « less
