Search for: All records

Planet–planetesimal interactions cause a planet to migrate, manifesting as a random walk in semimajor axis. In models for Neptune’s migration involving a gravitational upheaval, this planetesimal-driven migration is a side-effect of the dynamical friction required to damp Neptune’s orbital eccentricity. This migration is noisy, potentially causing Trans-Neptunian Objects (TNOs) in mean motion resonance to be lost. With N-body simulations, we validate a previously derived analytic model for resonance retention and determine unknown coefficients. We identify the impact of random-walk (noisy) migration on resonance retention for resonances up to fourth order lying between 39  and 75 au. Using a population estimate for the weak 7:3 resonance from the well-characterized Outer Solar System Origins Survey (OSSOS), we rule out two cases: (1) a planetesimal disc distributed between 13.3 and 39.9 au with ≳ 30 Earth masses in today’s size distribution and Tmig ≳ 40 Myr and (2) a top-heavy size distribution with ≳2000 Pluto-sized TNOs and Tmig ≳10 Myr, where Tmig is Neptune’s migration time-scale. We find that low-eccentricity TNOs in the heavily populated 5:2 resonance are easily lost due to noisy migration. Improved observations of the low-eccentricity region of the 5:2 resonance and of weak mean motion resonances with Rubin Observatory’s Legacy Survey of Space and Time will provide better population estimates, allowing for comparison with our model’s retention fractions and providing strong evidence for or against Neptune’s random interactions with planetesimals.

 

We present a semi-analytic model for the growth, drift, desorption, and fragmentation of millimeter- to meter-sized particles in protoplanetary disks. Fragmentation occurs where particle collision velocities exceed critical fragmentation velocities. Using this criterion, we produce fragmentation regions in disk orbital radius–particle size phase space for particles with a range of material properties, structures, and compositions (including SiO₂, Mg₂SiO₄, H₂O, CO₂, and CO). For reasonable disk conditions, compact aggregate H₂O, CO₂, and CO ice particles do not reach destructive relative velocities and are thus not likely to undergo collisional fragmentation. Uncoated silicate particles are more susceptible to collisional destruction and are expected to fragment in the inner disk, consistent with previous work. We then calculate the growth, drift, and sublimation of small particles, initially located in the outer disk. We find that ice-coated particles can avoid fragmentation as they grow and drift inward under a substantial range of disk conditions, as long as the particles are aggregates composed of 0.1μm-sized monomers. Such particles may undergo runaway growth in disk regions abundant in H₂O or CO₂ice, depending on the assumed disk temperature structure. These results indicate that icy collisional growth to planetesimally relevant sizes may happen efficiently throughout a disk’s lifetime, and is particularly robust at early times when the disk’s dust-to-gas ratio is comparable to that of the interstellar medium.

 

We present a direct imaging study of V892 Tau, a young Herbig Ae/Be star with a close-in stellar companion and circumbinary disk. Our observations consist of images acquired via Keck II/NIRC2 with nonredundant masking and the pyramid wavefront sensor at$K\prime$band (2.12μm) and$L\prime$band (3.78μm). Sensitivity to low-mass accreting companions and cool disk material is high at$L\prime$band, while complimentary observations at$K\prime$band probe hotter material with higher angular resolution. These multiwavelength, multiepoch data allow us to differentiate the secondary stellar emission from disk emission and deeply probe the structure of the circumbinary disk at small angular separations. We constrain architectural properties of the system by fitting geometric disk and companion models to the$K\prime$- and$L\prime$-band data. From these models, we constrain the astrometric and photometric properties of the stellar binary and update the orbit, placing the tightest estimates to date on the V892 Tau orbital parameters. We also constrain the geometric structure of the circumbinary disk, and resolve a circumprimary disk for the first time.

 

We present the highest-angular-resolution infrared monitoring of LkCa 15, a young solar analog hosting a transition disk. This system has been the subject of a number of direct-imaging studies from the millimeter through the optical, which have revealed multiple protoplanetary disk rings as well as three orbiting protoplanet candidates detected in infrared continuum emission (one of which was simultaneously seen at Hα). We use high-angular-resolution infrared imaging from 2014 to 2020 to systematically monitor these infrared signals and determine their physical origin. We find that three self-luminous protoplanets cannot explain the positional evolution of the infrared sources since the longer time baseline images lack the coherent orbital motion that would be expected for companions. However, the data still strongly prefer a time-variable morphology that cannot be reproduced by static scattered-light disk models. The multiepoch observations suggest the presence of complex and dynamic substructures moving through the forward-scattering side of the disk at ∼20 au or quickly varying shadowing by closer-in material. We explore whether the previous Hαdetection of one candidate would be inconsistent with this scenario and in the process develop an analytical signal-to-noise penalty for Hαexcesses detected near forward-scattered light. Under these new noise considerations, the Hαdetection is not strongly inconsistent with forward scattering, making the dynamic LkCa 15 disk a natural explanation for both the infrared and Hαdata.

 

Lyman α transits have been detected from several nearby exoplanets and are one of our best insights into the atmospheric escape process. However, due to ISM absorption, we typically only observe the transit signature in the blue-wing, making them challenging to interpret. This challenge has been recently highlighted by non-detections from planets thought to be undergoing vigorous escape. Pioneering 3D simulations have shown that escaping hydrogen is shaped into a cometary tail receding from the planet. Motivated by this work, we develop a simple model to interpret Lyman α transits. Using this framework, we show that the Lyman α transit depth is primarily controlled by the properties of the stellar tidal field rather than details of the escape process. Instead, the transit duration provides a direct measurement of the velocity of the planetary outflow. This result arises because the underlying physics is the distance a neutral hydrogen atom can travel before it is photoionized in the outflow. Thus, higher irradiation levels, expected to drive more powerful outflows, produce weaker, shorter Lyman α transits because the outflowing gas is ionized more quickly. Our framework suggests that the generation of energetic neutral atoms may dominate the transit signature early, but the acceleration of planetary material produces long tails. Thus, Lyman α transits do not primarily probe the mass-loss rates. Instead, they inform us about the velocity at which the escape mechanism is ejecting material from the planet, providing a clean test of predictions from atmospheric escape models.

 

The role of charge exchange in shaping exoplanet photoevaporation remains a topic of contention. Exchange of electrons between stellar wind protons from the exoplanet’s host star and neutral hydrogen from the planet’s wind has been proposed as a mechanism to create ‘energetic neutral atoms’ (ENAs), which could explain the high absorption line velocities observed in systems where mass-loss is occurring. In this paper, we present results from three-dimensional hydrodynamic simulations of the mass-loss of a planet similar to HD 209458b. We self-consistently launch a planetary wind by calculating the ionization and heating resulting from incident high-energy radiation, inject a stellar wind into the simulation, and allow electron exchange between the stellar and planetary winds. We predict the potential production of ENAs by the wind–wind interaction analytically, and then present the results of our simulations, which confirm the analytic limits. Within the limits of our hydrodynamic simulation, we find that charge exchange with the stellar wind properties examined here is unable to explain the absorption observed at high Doppler velocities.

 

ABSTRACT The role of radiation pressure in shaping exoplanet photoevaporation remains a topic of contention. Radiation pressure from the exoplanet’s host star has been proposed as a mechanism to drive the escaping atmosphere into a ‘cometary’ tail and explain the high velocities observed in systems where mass-loss is occurring. In this paper, we present results from high-resolution 3D hydrodynamic simulations of a planet similar to HD 209458b. We self-consistently launch a wind flowing outwards from the planet by calculating the ionization and heating resulting from incident high-energy radiation, and account for radiation pressure. We first present a simplified calculation, setting a limit on the Lyman-α flux required to drive the photoevaporated planetary material to larger radii and line-of-sight velocities. We then present the results of our simulations, which confirm the limits determined by our analytic calculation. We thus demonstrate that, within the limits of our hydrodynamic simulation and for the Lyman-α fluxes expected for HD 209458, radiation pressure is unlikely to significantly affect photoevaporative winds or to explain the high velocities at which wind material is observed, though further possibilities remain to be investigated.