Search

Ly  |$\alpha$| transits provide an opportunity to test models of atmospheric escape directly. However, translating observations into constraints on the properties of the escaping atmosphere is challenging. The major reason for this is that the observable parts of the outflow often comes from material outside the planet's Hill sphere, where the interaction between the planetary outflow and circumstellar environment is important. As a result, 3D models are required to match observations. Whilst 3D hydrodynamic simulations are able to match observational features qualitatively, they are too computationally expensive to perform a statistical retrieval of properties of the outflow. Here, we develop a model that determines the trajectory, ionization state, and 3D geometry of the outflow as a function of its properties and system parameters. We then couple this model to a ray tracing routine in order to produce synthetic transits. We demonstrate the validity of this approach, reproducing the trajectory of the outflows seen in 3D simulations. We illustrate the use of this model by performing a retrieval on the transit spectrum of GJ 436 b. The bound on planetary outflow velocity and mass-loss rates are consistent with a photoevaporative wind. [ABSTRACT FROM AUTHOR]

The atmospheres of small, close-in exoplanets are vulnerable to rapid mass loss during protoplanetary disc dispersal via a process referred to as 'boil-off', in which confining pressure from the local gas disc reduces, inducing atmospheric loss and contraction. We construct self-consistent models of planet evolution during gaseous core accretion and boil-off. As the surrounding disc gas dissipates, we find that planets lose mass via subsonic breeze outflows which allow causal contact to exist between disc and planet. Planets initially accrete of order |$\sim 10~{{\ \rm per\ cent}}$| in atmospheric mass, however, boil-off can remove |$\gtrsim 90~{{\ \rm per\ cent}}$| of this mass during disc dispersal. We show that a planet's final atmospheric mass fraction is strongly dictated by the ratio of cooling time-scale to disc dispersal time-scale, as well as the planet's core mass and equilibrium temperature. With contributions from core cooling and radioactivity, we show that core luminosity eventually leads to the transition from boil-off to core-powered mass loss. We find that smaller mass planets closest to their host star may have their atmospheres completely stripped through a combination of boil-off and core-powered mass loss during disc dispersal, implying the existence of a population-level radius gap emerging as the disc disperses. We additionally consider the transition from boil-off/core-powered mass loss to X-ray and extreme ultraviolet (XUV) photoevaporation by considering the penetration of stellar XUV photons below the planet's sonic surface. Finally, we show that planets may open gaps in their protoplanetary discs during the late stages of boil-off, which may enhance mass-loss rates. [ABSTRACT FROM AUTHOR]

Understanding atmospheric escape in close-in exoplanets is critical to interpreting their evolution. We map out the parameter space over which photoevaporation and core-powered mass-loss dominate atmospheric escape. Generally, the transition between the two regimes is determined by the location of the Bondi radius (i.e. the sonic point of core-powered outflow) relative to the penetration depth of extreme ultra-violet (XUV) photons. Photoevaporation dominates the loss when the XUV penetration depth lies inside the Bondi radius (R XUV < R B) and core-powered mass-loss when XUV radiation is absorbed higher up in the flow (R B < R XUV). The transition between the two regimes occurs at a roughly constant ratio of the planet's radius to its Bondi radius, with the exact value depending logarithmically on planetary and stellar properties. In general, core-powered mass-loss dominates for lower gravity planets with higher equilibrium temperatures, and photoevaporation dominates for higher gravity planets with lower equilibrium temperatures. However, planets can transition between these two mass-loss regimes during their evolution, and core-powered mass-loss can 'enhance' photoevaporation over a significant region of parameter space. Interestingly, a planet that is ultimately stripped by core-powered mass-loss has likely only ever experienced core-powered mass-loss. In contrast, a planet that is ultimately stripped by photoevaporation could have experienced an early phase of core-powered mass-loss. Applying our results to the observed super-Earth population suggests that it contains significant fractions of planets where each mechanism controlled the final removal of the H/He envelope, although photoevaporation appears to be responsible for the final carving of the exoplanet radius valley. [ABSTRACT FROM AUTHOR]

Planetary magnetic fields can affect the predicted mass-loss rate for close-in planets that experience large amounts of ultraviolet irradiation. In this work, we present a method to detect the magnetic fields of close-in exoplanets undergoing atmospheric escape using transit spectroscopy at the 10 830 Å line of helium. Motivated by previous work on hydrodynamic and magnetohydrodynamic photoevaporation, we suggest that planets with magnetic fields that are too weak to control the outflow's topology lead to blueshifted transits due to dayside-to-nightside flows. In contrast, strong magnetic fields prevent this day-to-night flow, as the gas is forced to follow the magnetic field's roughly dipolar topology. We post-process existing 2D photoevaporation simulations, computing synthetic transit profiles in helium to test this concept. As expected, we find that hydrodynamically dominated outflows lead to blueshifted transits of the order of the sound speed of the gas. Strong surface magnetic fields lead to unshifted or slightly redshifted transit profiles. High-resolution observations can distinguish between these profiles; however, eccentricity uncertainties generally mean that we cannot conclusively say that velocity shifts are due to the outflow for individual planets. The majority of helium observations are blueshifted, which could be a tentative indication that close-in planets generally have surface dipole magnetic field strengths |$\lesssim \!\! 0.3$| G. More 3D hydrodynamic and magnetohydrodynamic simulations are needed to confirm this conclusion robustly. [ABSTRACT FROM AUTHOR]

Atmospheric mass-loss is known to play a leading role in sculpting the demographics of small, close-in exoplanets. Knowledge of how such planets evolve allows one to 'rewind the clock' to infer the conditions in which they formed. Here, we explore the relationship between a planet's core mass and its atmospheric mass after protoplanetary disc dispersal by exploiting XUV photoevaporation as an evolutionary process. Historically, this inference problem would be computationally infeasible due to the large number of planet models required; however, we use a novel atmospheric evolution emulator which utilizes neural networks to provide three orders of magnitude in speedup. First, we provide a proof of concept for this emulator on a real problem by inferring the initial atmospheric conditions of the TOI-270 multi-planet system. Using the emulator, we find near-indistinguishable results when compared to the original model. We then apply the emulator to the more complex inference problem, which aims to find the initial conditions for a sample of Kepler, K2 , and TESS planets with well-constrained masses and radii. We demonstrate that there is a relationship between core masses and the atmospheric mass they retain after disc dispersal. This trend is consistent with the 'boil-off' scenario, in which close-in planets undergo dramatic atmospheric escape during disc dispersal. Thus, it appears that the exoplanet population is consistent with the idea that close-in exoplanets initially acquired large massive atmospheres, the majority of which is lost during disc dispersal, before the final population is sculpted by atmospheric loss over 100 Myr to Gyr time-scales. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Ultrashort period planets offer a window into the poorly understood interior composition of exoplanets through material evaporated from their rocky interiors. Among these objects are a class of disintegrating planets, observed when their dusty tails transit in front of their host stars. These dusty tails are thought to originate from dust condensation in thermally driven winds emanating from the sublimating surfaces of these planets. Existing models of these winds have been unable to explain their highly variable nature and have not explicitly modelled how dust forms in the wind. Here, we present new radiation-hydrodynamic simulations of the winds from these planets, including a minimal model for the formation and destruction of dust, assuming that nucleation can readily take place. We find that dust forms readily in the winds, a consequence of large dust grains obtaining lower temperatures than the planet's surface. As hyphothesized previously, we find that the coupling of the planet's surface temperature to the outflow properties via the dust's opacity can drive time-variable flows when dust condensation is sufficiently fast. In agreement with previous work, our models suggest that these dusty tails are a signature of catastrophically evaporating planets that are close to the end of their lives. Finally, we discuss the implications of our results for the dust's composition. More detailed hydrodynamic models that self-consistently compute the nucleation and composition of the dust and gas are warranted in order to use these models to study the planet's interior composition. [ABSTRACT FROM AUTHOR]

The radius distribution of small, close-in exoplanets has recently been shown to be bimodal. The photoevaporation model predicted this bimodality. In the photoevaporation scenario, some planets are completely stripped of their primordial H/He atmospheres, whereas others retain them. Comparisons between the photoevaporation model and observed planetary populations have the power to unveil details of the planet population inaccessible by standard observations, such as the core mass distribution and core composition. In this work, we present a hierarchical inference analysis on the distribution of close-in exoplanets using forward models of photoevaporation evolution. We use this model to constrain the planetary distributions for core composition, core mass, and initial atmospheric mass fraction. We find that the core-mass distribution is peaked, with a peak-mass of ∼4 M ⊕. The bulk core-composition is consistent with a rock/iron mixture that is ice-poor and 'Earth-like'; the spread in core-composition is found to be narrow (⁠|$\lesssim 16{{\ \rm per\ cent}}$| variation in iron-mass fraction at the 2 σ level) and consistent with zero. This result favours core formation in a water/ice poor environment. We find the majority of planets accreted a H/He envelope with a typical mass fraction of |$\sim 4{{\ \rm per\ cent}}$|⁠ ; only a small fraction did not accrete large amounts of H/He and were 'born-rocky'. We find four times as many super-Earths were formed through photoevaporation, as formed without a large H/He atmosphere. Finally, we find core-accretion theory overpredicts the amount of H/He cores would have accreted by a factor of ∼5, pointing to additional mass-loss mechanisms (e.g. 'boil-off') or modifications to core-accretion theory. [ABSTRACT FROM AUTHOR]

The photoevaporation model is one of the leading explanations for the evolution of small, close-in planets and the origin of the radius-valley. However, without planet mass measurements, it is challenging to test the photoevaporation scenario. Even if masses are available for individual planets, the host star's unknown EUV/X-ray history makes it difficult to assess the role of photoevaporation. We show that systems with multiple transiting planets are the best in which to rigorously test the photoevaporation model. By scaling one planet to another in a multitransiting system, the host star's uncertain EUV/X-ray history can be negated. By focusing on systems that contain planets that straddle the radius-valley, one can estimate the minimum masses of planets above the radius-valley (and thus are assumed to have retained a voluminous hydrogen/helium envelope). This minimum mass is estimated by assuming that the planet below the radius-valley entirely lost its initial hydrogen/helium envelope, then calculating how massive any planet above the valley needs to be to retain its envelope. We apply this method to 104 planets above the radius gap in 73 systems for which precise enough radii measurements are available. We find excellent agreement with the photoevaporation model. Only two planets (Kepler -100c and 142c) appear to be inconsistent, suggesting they had a different formation history or followed a different evolutionary pathway to the bulk of the population. Our method can be used to identify TESS systems that warrant radial-velocity follow-up to further test the photoevaporation model. The software to estimate minimum planet masses is publicly available at https://github.com/jo276/EvapMass. [ABSTRACT FROM AUTHOR]

The ubiquity of M dwarfs, combined with the relative ease of detecting terrestrial-mass planets around them, has made them prime targets for finding and characterizing planets in the 'habitable zone' (HZ). However, Kepler finds that terrestrial-mass exoplanets are often born with voluminous H/He envelopes, comprising mass-fractions (Menv/Mcore) ≳1 per cent. If these planets retain such envelopes over Gyr time-scales, they will not be 'habitable' even within the HZ. Given the strong X-ray/UV fluxes of M dwarfs, we study whether sufficient envelope mass can be photoevaporated away for these planets to become habitable. We improve upon previous work by using hydrodynamic models that account for radiative cooling as well as the transition from hydrodynamic to ballistic escape. Adopting a template active M dwarf XUV spectrum, including stellar evolution, and considering both evaporation and thermal evolution, we show that: (1) the mass-loss is (considerably) lower than previous estimates that use an 'energy-limited' formalism and ignore the transition to Jeans escape; (2) at the inner edge of the HZ, planets with core mass ≲ 0.9M⊕ can lose enough H/He to become habitable if their initial envelope mass-fraction is ~1 per cent; (3) at the outer edge of the HZ, evaporation cannot remove a ~1 per cent H/He envelope even from cores down to 0.8M⊕. Thus, if planets form with bulky H/He envelopes, only those with low-mass cores may eventually be habitable. Cores ≳1 M⊕, with ≳1 per cent natal H/He envelopes, will not be habitable in the HZ of M dwarfs. [ABSTRACT FROM AUTHOR]

We explore the dynamics of magnetically controlled outflows from hot Jupiters, where these flows are driven by UV heating from the central star. In these systems, some of the open field lines do not allow the flow to pass smoothly through the sonic point, so that steady-state solutions do not exist in general. This paper focuses on this type of magnetic field configuration, where the resulting flow becomes manifestly time-dependent. We consider the case of both steady heating and time-variable heating, and find the time-scales for the corresponding time variations of the outflow. Because the flow cannot pass through the sonic transition, it remains subsonic and leads to so-called breeze solutions. One manifestation of the time variability is that the flow samples a collection of different breeze solutions over time, and the mass outflow rate varies in quasi-periodic fashion. Because the flow is subsonic, information can propagate inwards from the outer boundary, which determines, in part, the time-scale of the flow variability. This work finds the relationship between the outer boundary scale and the time-scale of flow variations. In practice, the location of the outer boundary is set by the extent of the sphere of influence of the planet. The measured time variability can be used, in principle, to constrain the parameters of the system (e.g. the strengths of the surface magnetic fields). [ABSTRACT FROM AUTHOR]

Winds from short-period Earth and Neptune mass exoplanets, driven by high-energy radiation from a young star, may evaporate a significant fraction of a planet's mass. If the momentum flux from the evaporative wind is not aligned with the planet/star axis, then it can exert a torque on the planet's orbit. Using steady-state one-dimensional evaporative wind models, we estimate this torque using a lag angle that depends on the product of the speed of the planet's upper atmosphere and a flow time-scale for the wind to reach its sonic radius. We estimate the regime of planet radius, mass and stellar radiation flux in which a wind is capable of exerting a significant torque on the planet's orbit, and we find that it could be important for some of the observed planets. We also estimate the momentum flux from time-dependent onedimensional hydrodynamical simulations. Similar to the Yarkovsky effect, the wind causes the planet to drift outwards if atmospheric circulation is prograde (super-rotating) and in the opposite direction if the circulation is retrograde. A close-in super-Earth mass planet that loses a large fraction of its mass in a wind could drift a few per cent of its semimajor axis. While this change is small, it places constraints on the evolution of resonant pairs such as Kepler 36b and c. [ABSTRACT FROM AUTHOR]

We consider the role magnetic fields play in guiding and controlling mass-loss via evaporative outflows from exoplanets that experience UV irradiation. First, we present analytic results that account for planetary and stellar magnetic fields, along with mass-loss from both the star and planet. We then conduct series of numerical simulations for gas giant planets, and vary the planetary field strength, background stellar field strength, UV heating flux, and planet mass. These simulations show that the flow is magnetically controlled for moderate field strengths and even the highest UV fluxes, i.e. planetary surface fields BP ≳ 0.3 G and fluxes FUV ∼ 106 erg s−1. We thus conclude that outflows from all hot Jupiters with moderate surface fields are magnetically controlled. The inclusion of magnetic fields highly suppresses outflow from the night side of the planet. Only the magnetic field lines near the pole are open and allow outflow to occur. The fraction of open field lines depends sensitively on the strength (and geometry) of the background magnetic field from the star, along with the UV heating rate. The net effect of the magnetic field is to suppress the mass-loss rate by (approximately) an order of magnitude. Finally, some open field lines do not allow the flow to pass smoothly through the sonic point; flow along these streamlines does not reach a steady state, resulting in time-variable mass-loss. [ABSTRACT FROM AUTHOR]