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]

Stars do not form in isolation but together with other stars, and often in a clustered environment. Depending on the initial conditions in these environments, such as initial density and substructure, the distances of encounters between stars will differ. These encounters can also affect just-formed exoplanetary systems. Using N -body simulations, we show the effect of a single fly-by on a common type of exoplanetary system: close-in super-Earths/sub-Neptunes with or without a distant giant planet. Even a single encounter can significantly modify the architecture of these exoplanetary systems over their long lifetimes. We test fly-bys with different characteristics, such as distance and mass, and show how they perturb the inner planets long after the encounter, leading to collisions and mutual inclination excitation, which can significantly modify the observed architecture of these systems in transit. We find that our initially four-planet inner systems reduce to three or two inner planets depending on their initial separation and that the mutual inclinations of these remaining planets can be high enough to reduce the number of observable, transiting planets. In our 500 Myr simulations, we show that this reduction in the number of transiting planets due to stellar fly-bys can contribute to the observed excess of single-transit systems. [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]

Catastrophically evaporating rocky planets provide a unique opportunity to study the composition of small planets. The surface composition of these planets can be constrained via modelling their comet-like tails of dust. In this work, we present a new self-consistent model of the dusty tails: we physically model the trajectory of the dust grains after they have left the gaseous outflow, including an on-the-fly calculation of the dust cloud's optical depth. We model two catastrophically evaporating planets: KIC 1255 b and K2-22 b. For both planets, we find the dust is likely composed of magnesium–iron silicates (olivine and pyroxene), consistent with an Earth-like composition. We constrain the initial dust grain sizes to be ∼ 1.25–1.75 μm and the average (dusty) planetary mass-loss rate to be ∼ 3 |$\, M_{\oplus } \mathrm{Gyr^{-1}}$|⁠. Our model shows that the origin of the leading tail of dust of K2-22 b is likely a combination of the geometry of the outflow and a low radiation pressure force to stellar gravitational force ratio. We find the optical depth of the dust cloud to be a factor of a few in the vicinity of the planet. Our composition constraint supports the recently suggested idea that the dusty outflows of these planets go through a greenhouse effect–nuclear winter cycle, which gives origin to the observed transit depth time variability. Magnesium–iron silicates have the necessary visible-to-infrared opacity ratio to give origin to this cycle in the high mass-loss state. [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]

Disc accretion rate onto low mass protostar FU Ori suddenly increased hundreds of times 85 yr ago and remains elevated to this day. We show that the sum of historic and recent observations challenges existing FU Ori models. We build a theory of a new process, Extreme Evaporation (EE) of young gas giant planets in discs with midplane temperatures of ≳ 30 000 K. Such temperatures are reached in the inner 0.1 AU during thermal instability bursts. In our 1D time-dependent code the disc and an embedded planet interact through gravity, heat, and mass exchange. We use disc viscosity constrained by simulations and observations of dwarf novae instabilities, and we constrain planet properties with a stellar evolution code. We show that dusty gas giants born in the outer self-gravitating disc reach the innermost disc in a ∼ O (104) yr with radius of ∼10 R J. We show that their EE rates are |$\gtrsim 10^{-5} {\rm {\rm M}_{\odot }}$|  yr−1; if this exceeds the background disc accretion activity then the system enters a planet-sourced mode. Like a stellar secondary in mass-transferring binaries, the planet becomes the dominant source of matter for the star, albeit for ∼ O (100) yr. We find that a ∼6 Jupiter mass planet evaporating in a disc fed at a time-averaged rate of |$\sim 10^{-6} {\rm {\rm M}_{\odot }}$|  yr−1 appears to explain all that we currently know about FU Ori accretion outburst. More massive planets and/or planets in older less massive discs do not experience EE process. Future FUOR modelling may constrain planet internal structure and evolution of the earliest discs. [ABSTRACT FROM AUTHOR]

The Galactic Centre is known to have undergone a recent star formation episode a few Myr ago, which likely produced many T Tauri stars hosting circumstellar discs. It has been suggested that these discs may be the compact and dusty ionized sources identified as 'G-clouds'. Given the Galactic Centre's hostile environment, we study the possible evolutionary pathways these discs experience. We compute new external photoevaporation models applicable to discs in the Galactic Centre that account for the subsonic launching of the wind and absorption of UV photons by dust. Using evolutionary disc calculations, we find that photoevaporation's rapid truncation of the disc causes them to accrete onto the central star rapidly. Ultimately, an accreting circumstellar disc has a lifetime ≲ 1 Myr, which would fail to live long enough to explain the G-clouds. However, we identify a new evolutionary pathway for circumstellar discs in the Galactic Centre. Removal of disc material by photoevaporation prevents the young star from spinning down due to magnetic braking, ultimately causing the rapidly spinning young star to torque the disc into a 'decretion disc' state which prevents accretion. At the same time, any planetary companion in the disc will trap dust outside its orbit, shutting down photoevaporation. The disc can survive for up to ∼10 Myr in this state. Encounters with other stars are likely to remove the planet on Myr time-scales, causing photoevaporation to restart, giving rise to a G-cloud signature. A giant planet fraction of |$\sim 10{{\ \rm per\ cent}}$| can explain the number of observed G-clouds. [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]

Axisymmetric dust rings containing tens to hundreds of Earth masses of solids have been observed in protoplanetary discs with (sub-)millimetre imaging. Here, we investigate the growth of a planetary embryo in a massive (150 M⊕) axisymmetric dust trap through dust and gas hydrodynamics simulations. When accounting for the accretion luminosity of the planetary embryo from pebble accretion, the thermal feedback on the surrounding gas leads to the formation of an anticyclonic vortex. Since the vortex forms at the location of the planet, this has significant consequences for the planet's growth: as dust drifts towards the pressure maximum at the centre of the vortex, which is initially co-located with the planet, a rapid accretion rate is achieved, in a distinct phase of 'vortex-assisted' pebble accretion. Once the vortex separates from the planet due to interactions with the disc, it accumulates dust, shutting off accretion on to the planet. We find that this rapid accretion, mediated by the vortex, results in a planet containing ≈100 M⊕ of solids. We follow the evolution of the vortex, as well as the efficiency with which dust grains accumulate at its pressure maximum as a function of their size, and investigate the consequences this has for the growth of the planet as well as the morphology of the protoplanetary disc. We speculate that this extreme formation scenario may be the origin of giant planets that are identified to be significantly enhanced in heavy elements. [ABSTRACT FROM AUTHOR]

The photoevaporation of protoplanetary discs by nearby massive stars present in their birth cluster plays a vital role in their evolution. Previous modelling assumes that the disc behaves like a classical Keplerian accretion disc out to a radius where the photoevaporative outflow is launched. There is then an abrupt change in the angular velocity profile, and the outflow is modelled by forcing the fluid parcels to conserve their specific angular momenta. Instead, we model externally photoevaporating discs using the slim disc formalism. The slim disc approach self-consistently includes the advection of radial and angular momentum as well as angular momentum redistribution by internal viscous torques. Our resulting models produce a smooth transition from a rotationally supported Keplerian disc to a photoevaporative driven outflow, where this transition typically occurs over ∼4–5 scale heights. The penetration of ultraviolet photons predominately sets the radius of the transition and the viscosity's strength plays a minor role. By studying the entrainment of dust particles in the outflow, we find a rapid change in the dust size and surface density distribution in the transition region due to the steep gas density gradients present. This rapid change in the dust properties leaves a potentially observable signature in the continuum spectral index of the disc at mm wavelengths. Using the slim disc formalism in future evolutionary calculations will reveal how both the gas and dust evolve in their outer regions and the observable imprints of the external photoevaporation process. [ABSTRACT FROM AUTHOR]

Short-period super-Earth-sized planets are common. Explaining how they form near their present orbits requires understanding the structure of the inner regions of protoplanetary discs. Previous studies have argued that the hot inner protoplanetary disc is unstable to the magnetorotational instability (MRI) due to thermal ionization of potassium, and that a local gas pressure maximum forms at the outer edge of this MRI-active zone. Here we present a steady-state model for inner discs accreting viscously, primarily due to the MRI. The structure and MRI-viscosity of the inner disc are fully coupled in our model; moreover, we account for many processes omitted in previous such models, including disc heating by both accretion and stellar irradiation, vertical energy transport, realistic dust opacities, dust effects on disc ionization, and non-thermal sources of ionization. For a disc around a solar-mass star with a standard gas accretion rate (⁠|$\dot{M}\, \sim \, 10^{-8}$|  M⊙ yr−1) and small dust grains, we find that the inner disc is optically thick, and the accretion heat is primarily released near the mid-plane. As a result, both the disc mid-plane temperature and the location of the pressure maximum are only marginally affected by stellar irradiation, and the inner disc is also convectively unstable. As previously suggested, the inner disc is primarily ionized through thermionic and potassium ion emission from dust grains, which, at high temperatures, counteract adsorption of free charges on to grains. Our results show that the location of the pressure maximum is determined by the threshold temperature above which thermionic and ion emission become efficient. [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]

Recently, K2 and TESS have discovered transiting planets with radii between ∼5 and 10  R ⊕ around stars with ages <100 Myr. These young planets are likely to be the progenitors of the ubiquitous super-Earths/sub-Neptunes, which are well studied around stars with ages ≳1 Gyr. The formation and early evolution of super-Earths/sub-Neptunes are poorly understood. Various planetary origin scenarios predict a wide range of possible formation entropies. We show how the formation entropies of young (∼20–60 Myr), highly irradiated planets can be constrained if their mass, radius, and age are measured. This method works by determining how low-mass an H/He envelope a planet can retain against mass-loss, this lower bound on the H/He envelope mass can then be converted into an upper bound on the entropy. If planet mass measurements with errors ≲20 per cent can be achieved for the discovered young planets around DS Tuc A and V1298 Tau, then insights into their origins can be obtained. For these planets, higher measured planet masses would be consistent with the standard core-accretion theory. In contrast, lower planet masses (≲6–7 M⊕) would require a 'boil-off' phase during protoplanetary disc dispersal to explain. [ABSTRACT FROM AUTHOR]

The Sun shows a ∼10 per cent depletion in refractory elements relative to nearby solar twins. It has been suggested that this depletion is a signpost of planet formation. The exoplanet statistics are now good enough to show that the origin of this depletion does not arise from the sequestration of refractory material inside the planets themselves. This conclusion arises because most sun-like stars host close-in planetary systems that are on average more massive than the Sun's. Using evolutionary models for the protoplanetary discs that surrounded the young Sun and solar twins, we demonstrate that the origin of the depletion likely arises due to the trapping of dust exterior to the orbit of a forming giant planet. In this scenario, a forming giant planet opens a gap in the gas disc, creating a pressure trap. If the planet forms early enough, while the disc is still massive, the planet can trap ≳100 M⊕ of dust exterior to its orbit, preventing the dust from accreting on to the star in contrast to the gas. Forming giant planets can create refractory depletions of |$\sim 5{-}15{{\ \rm per\ cent}}$|⁠ , with the larger values occurring for initial conditions that favour giant planet formation (e.g. more massive discs that live longer). The incidence of solar twins that show refractory depletion matches both the occurrence of giant planets discovered in exoplanet surveys and 'transition' discs that show similar depletion patterns in the material that is accreting on to the star. [ABSTRACT FROM AUTHOR]

We demonstrate the utility of C  i as a tracer of photoevaporative winds that are being driven from discs by their ambient UV environment. Commonly observed CO lines only trace these winds in relatively weak UV environments and are otherwise dissociated in the wind at the intermediate to high UV fields that most young stars experience. However, C  i traces unsubtle kinematic signatures of a wind in intermediate UV environments (∼1000 G0) and can be used to place constraints on the kinematics and temperature of the wind. In C  i position–velocity (PV) diagrams external photoevaporation results in velocities that are faster than those from Keplerian rotation alone, as well as emission from quadrants of PV space in which there would be no Keplerian emission. This is independent of viewing angle because the wind has components that are perpendicular to the azimuthal rotation of the disc. At intermediate viewing angles (∼30–60°) moment 1 maps also exhibit a twisted morphology over large scales (unlike other processes that result in twists, which are typically towards the inner disc). C  i is readily observable with ALMA, which means that it is now possible to identify and characterize the effect of external photoevaporation on planet-forming discs in intermediate UV environments. [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 observed distribution of radii for exoplanets shows a bimodal form that can be explained by mass-loss from planetary atmospheres due to high-energy radiation emitted by their host stars. The location of the minimum of this radius distribution depends on the mass–radius relation, which in turn depends on the composition of planetary cores. Current studies suggest that super-Earth and mini-Neptune planets have iron-rich and hence largely Earth-like composition cores. This paper explores how non-zero planetary magnetic fields can decrease the expected mass-loss rates from these planets. These lower mass-loss rates, in turn, affect the location of the minimum of the radius distribution and the inferred chemical composition of the planetary cores. [ABSTRACT FROM AUTHOR]

Theoretical models of protoplanetary disc dispersal predict a phase where photoevaporation has truncated the disc at several au, creating a pressure trap which is dust-rich. Previous models predicted this phase could be long-lived (∼Myr), contrary to the observational constraints. We show that dust in the pressure trap can be removed from the disc by radiation pressure exerting a significant acceleration, and hence radial velocity, on small dust particles that reside in the surface layers of the disc. The dust in the pressure trap is not subject to radial drift so it can grow to reach sizes large enough to fragment. Hence small particles removed from the surface layers are replaced by the fragments of larger particles. This link means radiation pressure can deplete the dust at all particle sizes. Through a combination of 1D and 2D models, along with secular models that follow the disc's long-term evolution, we show that radiation pressure can deplete dust from pressure traps created by photoevaporation in ∼105 yr, while the photoevaporation created cavity still resides at 10 s of au. After this phase of radiation pressure removal of dust, the disc is gas-rich and dust depleted and radially optically thin to stellar light, having observational signatures similar to a gas-rich, young debris disc. Indeed many of the young stars (≲10 Myr old) classified as hosting a debris disc may rather be discs that have undergone this process. [ABSTRACT FROM AUTHOR]

Exoplanets with substantial hydrogen/helium atmospheres have been discovered in abundance, many residing extremely close to their parent stars. The extreme irradiation levels that these atmospheres experience cause them to undergo hydrodynamic atmospheric escape. Ongoing atmospheric escape has been observed to be occurring in a few nearby exoplanet systems through transit spectroscopy both for hot Jupiters and for lower-mass super-Earths and mini-Neptunes. Detailed hydrodynamic calculations that incorporate radiative transfer and ionization chemistry are now common in one-dimensional models, and multidimensional calculations that incorporate magnetic fields and interactions with the interstellar environment are cutting edge. However, comparison between simulations and observations remains very limited. While hot Jupiters experience atmospheric escape, the mass-loss rates are not high enough to affect their evolution. However, for lower-mass planets, atmospheric escape drives and controls their evolution, sculpting the exoplanet population that we observe today. ▪ Observations of some exoplanets have detected atmospheric escape driven by hydrodynamic outflows, causing the exoplanets to lose mass over time. ▪ Hydrodynamic simulations of atmospheric escape are approaching the sophistication required to compare them directly to observations. ▪ Atmospheric escape sculpts sharp features into the exoplanet population that we can observe today; these features have recently been detected. [ABSTRACT FROM AUTHOR]

We explore in situ formation and subsequent evolution of close-in super-Earths and mini-Neptunes. We adopt a steady-state inner protoplanetary gas disc structure that arises from viscous accretion due to the magneto-rotational instability (MRI). We consider the evolution of dust in the inner disc, including growth, radial drift, and fragmentation, and find that dust particles that radially drift into the inner disc fragment severely due to the MRI-induced turbulence. This result has two consequences: (1) radial drift of grains within the inner disc is quenched, leading to an enhancement of dust in the inner regions that scales as dust-to-gas-mass-flux-ratio at ∼1 au; (2) however, despite this enhancement, planetesimal formation is impeded by the small grain size. Nevertheless, assuming that planetary cores are present in the inner disc, we then investigate the accretion of atmospheres on to cores and their subsequent photoevaporation. We then compare our results to the observed exoplanet mass–radius relationship. We find that (1) the low gas surface densities and high temperatures in the inner disc reduce gas accretion on to cores compared to the minimum mass solar nebula, preventing the cores from growing into hot Jupiters, in agreement with the data; (2) however, our predicted envelope masses are still typically larger than observed ones. Finally, we sketch a qualitative picture of how grains may grow and planetesimals form in the inner disc if grain effects on the ionization levels and the MRI and the back reaction of the dust on the gas (both neglected in our calculations) are accounted for. [ABSTRACT FROM AUTHOR]

We present a new two-dimensional (2D) axisymmetric code, cuDisc , for studying protoplanetary discs, focusing on the self-consistent calculation of dust dynamics, grain-size distribution and disc temperature. Self-consistently studying these physical processes is essential for many disc problems, such as structure formation and dust removal, given that the processes heavily depend on one another. To follow the evolution over substantial fractions of the disc lifetime, cuDisc uses the cuda language and libraries to speed up the code through GPU acceleration. cuDisc employs a second-order finite-volume Godonuv solver for dust dynamics, solves the Smoluchowski equation for dust growth, and calculates radiative transfer using a multifrequency hybrid ray-tracing/flux-limited-diffusion method. We benchmark our code against current state-of-the-art codes. Through studying steady-state problems, we find that including 2D structure reveals that when collisions are important, the dust vertical structure appears to reach a diffusion-settling-coagulation equilibrium that can differ substantially from standard models that ignore coagulation. For low fragmentation velocities, we find an enhancement of intermediate-sized dust grains at heights of ∼1 gas scale height due to the variation in collision rates with height, and for large fragmentation velocities, we find an enhancement of small grains around the disc mid-plane due to collisional 'sweeping' of small grains by large grains. These results could be important for the analysis of disc spectral energy distributions or scattered light images, given these observables are sensitive to the vertical grain distribution. [ABSTRACT FROM AUTHOR]

Using data from the California-Kepler Survey (CKS), we study trends in planetary properties with host star metallicity for close-in planets. By incorporating knowledge of the properties of the planetary radius gap identified by the CKS survey, we are able to investigate the properties of planetary cores and their gaseous envelopes separately. Our primary findings are that the solid core masses of planets are higher around higher metallicity stars and that these more massive cores were able to accrete larger gas envelopes. Furthermore, investigating the recently reported result that planets with radii in the range ( 2–6 R⊕) are more common at short periods around higher metallicity stars in detail, we find that the average host star metallicity of H/He atmosphere-hosting planets increases smoothly inside an orbital period of ∼20 d. We interpret the location of the metallicity increase within the context of atmospheric photoevaporation: higher metallicity stars are likely to host planets with higher atmospheric metallicity, which increases the cooling in the photoevaporative outflow, lowering the mass-loss rates. Therefore, planets with higher metallicity atmospheres are able to resist photoevaporation at shorter orbital periods. Finally, we find evidence at 2.8 σ that planets that do not host H/He atmospheres at long periods are more commonly found around lower metallicity stars. Such planets are difficult to explain by photoevaporative stripping of planets which originally accreted H/He atmospheres. Alternatively, this population of planets could be representative of planets that formed in a terrestrial-like fashion, after the gas disc dispersed. [ABSTRACT FROM AUTHOR]

The mass–period or radius–period distribution of close-in exoplanets shows a paucity of intermediate mass/size (sub-Jovian) planets with periods ≲3 d. We show that this sub-Jovian desert can be explained by the photoevaporation of highly irradiated sub-Neptunes and the tidal disruption barrier for gas giants undergoing high-eccentricity migration. The distinctive triangular shape of the sub-Jovain desert results from the fact that photoevaporation is more effective closer to the host star, and that in order for a gas giant to tidally circularize closer to the star without tidal disruption it needs to be more massive. Our work indicates that super-Earths/mini-Neptunes and hot-Jupiters had distinctly separate formation channels and arrived at their present locations at different times. [ABSTRACT FROM AUTHOR]

In this work, we develop a rocky planet interior model and use it to investigate the evolution of catastrophically evaporating rocky exoplanets. These planets, detected through the dust tails produced by evaporative outflows from their molten surfaces, can be entirely destroyed in a fraction of their host star's lifetime. To allow for the major decrease in mass, our interior model can simultaneously calculate the evolution of the pressure and density structure of a planet alongside its thermal evolution, which includes the effects of conduction, convection and partial melting. We first use this model to show that the underlying planets are likely to be almost entirely solid. This means that the dusty tails are made up of material sampled only from a thin dayside lava pool. If one wishes to infer the bulk compositions of rocky exoplanets from their dust tails, it is important to take the localized origin of this material into account. Secondly, by considering how frequently one should be able to detect mass loss from these systems, we investigate the occurrence of sub-Earth mass exoplanets, which is difficult with conventional planet detection surveys. We predict that, depending on model assumptions, the number of progenitors of the catastrophically evaporating planets is either in line with, or higher than, the observed population of close-in (substellar temperatures around 2200 K) terrestrial exoplanets. [ABSTRACT FROM AUTHOR]

We present a simple model for low-mass planet formation and subsequent evolution within 'transition' discs. We demonstrate quantitatively that the predicted and observed structures of such discs are prime birthsites of planets. Planet formation is likely to proceed through pebble accretion, should a planetary embryo (M ≳ 10-4 M⊕) form. Efficient pebble accretion is likely to be unavoidable in transition disc dust traps, as the dust particles required for pebble accretion are those which are most efficiently trapped in the transition disc dust trap. Rapid pebble accretion within the dust trap gives rise not only to low-mass planets, but also to a large accretion luminosity. This accretion luminosity is sufficient to heat the disc outside the gravitational influence of the planet and makes the disc locally baroclinic, and a source of vorticity. Using numerical simulations, we demonstrate that this source of vorticity can lead to the growth of a single large-scale vortex in ~100 orbits, which is capable of trapping particles. Finally, we suggest an evolutionary cycle: Planet formation proceeds through pebble accretion, followed by vortex formation and particle trapping in the vortex quenching the planetary accretion and thus removing the vorticity source. After the vortex is destroyed, the process can begin anew. This means transition discs should present with large-scale vortices for a significant fraction of their lifetimes, and remnant planets at large 10 au radii should be a common outcome of this cycle. [ABSTRACT FROM AUTHOR]

Photoevaporation of protoplanetary discs by high-energy radiation from the central young stellar object is currently the favourite model to explain the sudden dispersal of discs from the inside out. While several theoretical works have provided a detailed pictured of this process, the direct observational validation is still lacking. Emission lines produced in these slow-moving protoplanetary disc winds may bear the imprint of the wind structure and thus provide a potential diagnostic of the underlying dispersal process. In this paper, we primarily focus on the collisionally excited neutral oxygen line at 6300 Å. We compare our models predictions to observational data and demonstrate a thermal origin for the observed blueshifted low-velocity component of this line from protoplanetary discs. Furthermore, our models show that while this line is a clear tell-tale sign of a warm, quasi-neutral disc wind, typical of X-ray photoevaporation, its strong temperature dependence makes it unsuitable to measure detailed wind quantities like mass-loss rate. [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]

Late-typemain-sequence stars exhibit anX-ray to bolometric flux ratio that depends on R̃o, the ratio of rotation period to convective turnover time, as R̃o -ζ with 2 ≤ ζ ≤ 3 for R̃o > 0.13, but saturates with ∣ζ ∣ < 0.2 for R̃o < 0.13. Saturated stars are younger than unsaturated stars and show a broader spread of rotation rates and X-ray activity. The unsaturated stars have magnetic fields and rotation speeds that scale roughly with the square root of their age, though possibly flattening for stars older than the Sun. The connection between faster rotators, stronger fields, and higher activity has been established observationally, but a theory for the unified timeevolution of X-ray luminosity, rotation, magnetic field and mass loss that captures the above trends has been lacking. Here we derive a minimalist holistic framework for the time evolution of these quantities built from combining a Parker wind with new ingredients: (1) explicit sourcing of both the thermal energy launching the wind and the X-ray luminosity via dynamo produced magnetic fields; (2) explicit coupling of X-ray activity and mass-loss saturation to dynamo saturation (via magnetic helicity build-up and convection eddy shredding); (3) use of coronal equilibrium to determine how magnetic energy is divided into wind and Xray contributions. For solar-type stars younger than the Sun, we infer conduction to be a subdominant power loss compared to X-rays and wind. For older stars, conduction is more important, possibly quenching the wind and reducing angular momentum loss. We focus on the time evolution for stars younger than the Sun, highlighting what is possible for further generalizations. Overall, the approach shows promise towards a unified explanation of all of the aforementioned observational trends. [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]

The D/H ratio of water in protostellar systems is a result of both inheritance from the parent molecular cloud and isotopic exchange in the disc. A possibly widespread feature of disc evolution, ignored in previous studies, is accretion bursts (or FU Orionis outbursts), which may thermally process a large fraction of the water. One proposed underlying mechanism for FU Orionis outbursts relies on the presence of a magnetically dead zone. Here we examine the evolution of (D/H)water in 1D simulations of a disc's evolution that includes dead zones and infall from an envelope with given D/H ratio in the infalling water (~10-3), and compare the results with similar calculations without dead zones. We find that the accretion bursts result in a significantly lower (D/H)water ratio and a more extended region (radius up to ~1-3 au) where water is equilibrated with hydrogen gas (D/H = 2 × 10-5), when compared to burst-free models. Solar system constraints suggest that our solar nebula either experienced no accretion bursts and had a Schmidt number ≲ 0.2 or had a Schmidt number closer to 'nominal' values (~1) and experienced several accretion bursts. Finally, future observations of (D/H)water in protoplanetary discs will allow inferences about angular momentum properties of the disc during disc building and the role of accretion bursts. [ABSTRACT FROM AUTHOR]

We show that if young low-mass stars are subject to vigorous X-ray driven disc winds, then such winds may be rendered detectable in cluster environments through their interaction with ionizing radiation from massive stars. In particular we argue that in the Orion nebula cluster (ONC) one expects to see of order tens of ‘X-ray proplyds’ (i.e. objects with offset ionization fronts detectable through optical imaging) in the range 0.3–0.6 pc from θ1C Ori (the dominant O star in the ONC). Objects at this distance lie outside the central ‘FUV zone’ in the ONC where proplyd structures are instead well explained by neutral winds driven by external far-ultraviolet (FUV) emission from θ1C Ori. We show that the predicted numbers and sizes of X-ray proplyds in this region are compatible with the numbers of proplyds observed and that this may also provide an explanation for at least some of the far flung proplyds observed in the Carina nebula. We compare the sizes of observed proplyds outside the FUV region of the ONC with model predictions based on the current observed X-ray luminosities of these sources (bearing in mind that the current size is actually set by the X-ray luminosity a few hundred years previously, corresponding to the flow time to the ionization front). We discuss whether variability on this time-scale can plausibly explain the proplyd size data on a case-by-case basis. We also calculate the predicted radio free–free emission signature of X-ray proplyds and show that this is readily detectable. Monitoring is, however, required in order to distinguish such emission from non-thermal radio emission from active coronae. We also predict that it is only at distances more than a parsec from θ1C Ori that the free–free emission signature of such offset ionized structures would be clearly distinguishable from an externally driven ionized disc wind. We argue that the fortuitous proximity of massive stars in the ONC can be used as a beacon to light up internally driven X-ray winds and that this represents a promising avenue for observational tests of the X-ray photoevaporation scenario. [ABSTRACT FROM AUTHOR]

We consider the role of thermal diffusion due to turbulence and radiation on accretion bursts that occur in protoplanetary discs which contain dead-zones. Using 1D viscous disc models, we show that diffusive radial transport of heat is important during the gravomagnetic limit cycle, and can strongly modify the duration and frequency of accretion outbursts. When the Prandtl number is large – such that turbulent diffusion of heat is unimportant – radial radiative diffusion reduces the burst duration compared to models with no diffusive transport of heat. When the Prandtl number is small ≲25, we find that turbulent diffusion dominates the radial transport of heat, reducing the burst duration to ≲103 yr as well as increasing the burst frequency. Furthermore, inclusion of radial transport of heat extends the range of infall rates under which the disc undergoes accretion bursts from 10−8 to 10−6 M⊙ yr−1 with no diffusion, to 10−8 to ≳10−4 M⊙ yr−1 with radiative and strong turbulent diffusion. The relative roles of radiative and turbulent thermal diffusion are likely to vary during an accretion burst, but simple estimates suggest that the expected Prandtl numbers are of the order of 10 in protoplanetary discs, and hence that turbulent diffusion is likely to be an important process during accretion outbursts due to the gravomagnetic limit cycle. [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]

Inspired by the Kepler mission's planet discoveries, we consider the thermal contraction of planets close to their parent star, under the influence of evaporation. The mass-loss rates are based on hydrodynamic models of evaporation that include both X-ray and EUV irradiation. We find that only low mass planets with hydrogen envelopes are significantly affected by evaporation, with evaporation being able to remove massive hydrogen envelopes inward of ∼0.1 AU for Neptune-mass objects, while evaporation is negligible for Jupiter-mass objects. Moreover, most of the evaporation occurs in the first 100 Myr of stars' lives when they are more chromospherically active. We construct a theoretical population of planets with varying core masses, envelope masses, orbital separations, and stellar spectral types, and compare this population with the sizes and densities measured for low-mass planets, both in the Kepler mission and from radial velocity surveys. This exercise leads us to conclude that evaporation is the driving force of evolution for close-in Kepler planets. In fact, some 50% of the Kepler planet candidates may have been significantly eroded. Evaporation explains two striking correlations observed in these objects: a lack of large radius/low density planets close to the stars and a possible bimodal distribution in planet sizes with a deficit of planets around 2 R⊕. Planets that have experienced high X-ray exposures are generally smaller than this size, and those with lower X-ray exposures are typically larger. A bimodal planet size distribution is naturally predicted by the evaporation model, where, depending on their X-ray exposure, close-in planets can either hold on to hydrogen envelopes ∼0.5%-1% in mass or be stripped entirely. To quantitatively reproduce the observed features, we argue that not only do low-mass Kepler planets need to be made of rocky cores surrounded with hydrogen envelopes, but few of them should have initial masses above 20 M⊕ and the majority of them should have core masses of a few Earth masses. [ABSTRACT FROM AUTHOR]

The spin evolution of main-sequence stars has long been of interest for basic stellar evolution, stellar ageing, stellar activity, and consequent influence on companion planets. Observations of older-than-solar late-type main-sequence stars have been interpreted to imply that a change from a dipole-dominated magnetic field to one with more prominent higher multipoles might be necessary to account for the data. The spin-down models that lead to this inference are essentially tuned to the Sun. Here, we take a different approach that considers individual stars as fixed points rather than just the Sun. We use a time-dependent theoretical model to solve for the spin evolution of low-mass main-sequence stars that includes a Parker-type wind and a time-evolving magnetic field coupled to the spin. Because the wind is exponentially sensitive to the stellar mass over radius and the coronal base temperature, the use of each observed star as a separate fixed point is more appropriate and, in turn, produces a set of solution curves that produces a solution envelope rather than a simple line. This envelope of solution curves, unlike a single line fit, is consistent with the data and does not unambiguously require a modal transition in the magnetic field to explain it. [ABSTRACT FROM AUTHOR]

ABSTRACT We examine the distribution of transition discs as a function of millimetre (mm) flux. We confirm that as expected in any model in which most primordial discs turn into transition discs and in which mm flux declines with time, transition discs have lower mm fluxes on average than primordial discs. However, we find that the incidence of transition discs does not, as expected, fall monotonically towards large mm fluxes and also investigate the hypothesis that these mm bright transition discs may have a distinct physical origin. We find that mm bright transition discs occupy a separate region of parameter space. Transition discs in the bright mm subsample have systematically higher accretion rates than those in the faint mm subsample, along with being systematically weighted to earlier spectral types. [ABSTRACT FROM AUTHOR]

ABSTRACT We discuss a hydrodynamical model for the dispersal of protoplanetary discs around young, low-mass (<1.5 M⊙) stars by photoevaporation from the central object's energetic radiation, which considers the far-ultraviolet (FUV) as well as the X-ray component of the radiation field. We present analytical scaling relations and derive estimates for the total mass-loss rates, as well as discussing the existence of similarity solutions for flows from primordial discs and discs with inner holes. Furthermore, we perform numerical calculations which span a wide range of parameter space and allow us to provide accurate scalings of the mass-loss rates with the physical parameters of the systems (X-ray and FUV luminosity, stellar mass, disc mass, disc temperature and inner hole radius). The model suggests that the X-ray component dominates the photoevaporative mass-loss rates from the inner disc. The mass-loss rates have values in the range from 10−11 to 10−7 M⊙ yr−1 and scale linearly with X-ray luminosity, with only a weak dependence on the other parameters explored. However, in the case of high FUV-to-X-ray luminosity ratios ( LFUV/ LX > 100), the FUV constricts the X-ray flow and may dominate the mass loss. Simulations of low-mass discs with inner holes demonstrate a further stage of disc clearing, which we call 'thermal sweeping'. This process occurs when the mid-plane pressure drops to sufficiently low values. At this stage, a bound, warm, X-ray heated region becomes sufficiently large and unstable, such that the remaining disc material is cleared on approximately dynamical time-scales. This process significantly reduces the time taken to clear the outer regions of the disc, resulting in an expected transition disc population that will be dominated by accreting objects, as indicated by recent observations. [ABSTRACT FROM AUTHOR]

We explore the role of X-ray photoevaporation in the evolution and dispersal of viscously evolving T Tauri discs. We show that the X-ray photoevaporation wind rates scale linearly with X-ray luminosity, such that the observed range of X-ray luminosities for solar-type T Tauri stars (10-10 erg s) gives rise to vigorous disc winds with rates of the order of 10 to 10 M yr. These mass-loss rates are comparable to typically observed T Tauri accretion rates, immediately demonstrating the relevance of X-ray photoevaporation to disc evolution. We use the wind solutions from radiation-hydrodynamic models, coupled to a viscous evolution model, to construct a population synthesis model so that we may study the physical properties of evolving discs and so-called 'transition discs'. Current observations of disc lifetimes and accretion rates can be matched by our model assuming a viscosity parameter . Our models confirm that X-rays play a dominant role in the evolution and dispersal of protoplanetary discs giving rise to the observed diverse population of inner-hole 'transition' sources which include those with massive outer discs, those with gas in their inner holes and those with detectable accretion signatures. To help understand the nature of observed transition discs we present a diagnostic diagram based on accretion rates versus inner-hole sizes that demonstrate that, contrary to recent claims, many of the observed accreting and non-accreting transition discs can easily be explained by X-ray photoevaporation. However, we draw attention to a smaller but still significant population of strongly accreting (10 M yr) transition discs with large inner holes (>20 au) that lie outside the predicted X-ray photoevaporation region, suggesting a different origin for their inner holes. Finally, we confirm the conjecture of that accretion is suppressed by the X-rays through 'photoevaporation-starved accretion' and predict that this effect can give rise to a negative correlation between X-ray luminosity and accretion rate, as reported in the Orion data. We also demonstrate that our model can replicate the observed difference in X-ray properties between accreting and non-accreting T Tauri stars. [ABSTRACT FROM AUTHOR]

We have performed hydrodynamic and radiative transfer calculations of a photoevaporating disc around a Herbig Ae/Be star to determine the evolution and observational impact of dust entrained in the wind. We find that the wind selectively entrains grains of different sizes at different radii resulting in a dust population that varies spatially and increases with height above the disc at radii >10 au. This variable grain population results in a 'wingnut' morphology to the dust density distribution. We calculate images of this dust distribution at near-infrared wavelengths that also show a 'wingnut' morphology at all wavelengths considered. We have also considered the contribution that crystalline dust grains will have in the wind and show that photoevaporative wind can result in a significant crystallinity fraction at all radii, when the disc is edge-on. However, when the disc's photosphere is unobscured, the photoevaporative wind makes no contribution to the observable crystallinity fraction in the disc. Finally, we conclude that the analysis of extended emission around edge-on discs could provide a new and independent method of testing photoevaporation models. [ABSTRACT FROM AUTHOR]

We present a calculation of the atomic and low-ionization emission line spectra of photoevaporating protoplanetary discs. Line luminosities and profiles are obtained from detailed photoionization calculations of the disc and wind structures surrounding young active solar-type stars. The disc and wind density and velocity fields were obtained from the recently developed radiation-hydrodynamic models of Owen et al. that include stellar X-ray and EUV irradiation of protoplanetary discs at various stages of clearing, from primordial sources to inner hole sources of various hole sizes. Our models compare favourably with currently available observations, lending support to an X-ray driven photoevaporation model for disc dispersal. In particular, we find that X-rays drive a warm, predominantly neutral flow where the O i 6300 Å line can be produced by neutral hydrogen collisional excitation. Our models can, for the first time, provide a very good match to both luminosities and profiles of the low-velocity component of the O i 6300 Å line and other forbidden lines observed by Hartigan et al., which covered a large sample of T-Tauri stars. We find that the O i 6300 Å and the Ne ii 12.8 μm lines are predominantly produced in the X-ray-driven wind and thus appear blueshifted by a few km s−1 for some of the systems when observed at non-edge-on inclinations. We note, however, that blueshifts are only produced under certain conditions: X-ray luminosity, spectral shape and inner hole size all affect the location of the emitting region and the physical conditions in the wind. We caution therefore that while a blueshifted line is a tell-tale sign of an outflow, the lack of a blueshift should not be necessarily interpreted as a lack of outflow. Comparison of spectrally resolved observations of multiple emission lines with detailed model sets, like the one presented here, should provide useful diagnostics of the clearing of gaseous discs. [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]

Hot Jupiters present a unique opportunity for measuring how planet formation history shapes present-day atmospheric composition. However, due to the myriad pathways influencing composition, a well-constructed sample of planets is needed to determine whether formation history can be accurately traced back from atmospheric composition. To this end, the BOWIE-ALIGN survey (A spectral Light Investigation into hot gas Giant origiNs by the collaboration of Bristol, Oxford, Warwick, Imperial, Exeter,  +) will compare the compositions of eight hot Jupiters around F stars, four with orbits aligned with the stellar rotation axis, and four misaligned. Using the alignment as an indicator for planets that underwent disc migration or high-eccentricity migration, one can determine whether migration history produces notable differences in composition between the two samples of planets. This paper describes the planet formation model that motivates our observing programme. Our model traces the accretion of chemical components from the gas and dust in the disc over a broad parameter space to create a full, unbiased model sample from which we can estimate the range of final atmospheric compositions. For high metallicity atmospheres (⁠|$\mathrm{ O}\mathrm{ /H}\ge 10 \times$|  solar), the C/O ratios of aligned and misaligned planets diverge, with aligned planets having lower C/O (⁠|$\lt 0.25$|⁠) due to the accretion of oxygen-rich silicates from the inner disc. However, silicates may rain out instead of releasing their oxygen into the atmosphere. This would significantly increase the C/O of aligned planets (C/O |$\gt 0.6$|⁠), inverting the trend between the aligned and misaligned planets. Nevertheless, by comparing statistically significant samples of aligned and misaligned planets, we expect atmospheric composition to constrain how planets form. [ABSTRACT FROM AUTHOR]

The radius valley, i.e. a dearth of planets with radii between 1.5 and 2 Earth radii, provides insights into planetary formation and evolution. Using homogenously revised planetary parameters from Kepler 1-min short cadence light curves, we remodel transits of 72 small planets mostly orbiting low-mass stars, improving the precision and accuracy of planet parameters. By combining this sample with a similar sample of planets around higher mass stars, we determine the depth of the radius valley as a function of stellar mass. We find that the radius valley is shallower for low-mass stars compared to their higher mass counterparts. Upon comparison, we find that theoretical models of photoevaporation underpredict the number of planets observed inside the radius valley for low-mass stars: with decreasing stellar mass, the predicted fraction of planets inside the valley remains approximately constant whereas the observed fraction increases. We argue that this provides evidence for the presence of icy planets around low-mass stars. Alternatively, planets orbiting low-mass stars undergo more frequent collisions and scatter in the stars' high-energy output may also cause planets to fill the valley. We predict that more precise mass measurements for planets orbiting low-mass stars would be able to distinguish between these scenarios. [ABSTRACT FROM AUTHOR]

This document is a correction to a previously published paper titled "Why the observed spin evolution of older-than-solar-like stars might not require a dynamo mode change." The correction addresses minor typographical errors that do not impact the results of the paper but could cause confusion in reproducing the calculations. Specifically, the correction involves changing the notation for stellar mass and radius in subsection 2.3. The corrected notation uses m⊙* and r⊙* to represent normalized values to solar mass and solar radius. The authors of the paper are Ketevan Kotorashvili, Eric G Blackman, and James E Owen. [Extracted from the article]

Photometric observations of accreting, low-mass, pre-main-sequence stars (i.e., Classical T Tauri stars; CTTS) have revealed different categories of variability. Several of these classifications have been linked to changes in. To test how accretion variability conditions lead to different light-curve morphologies, we used 1D hydrodynamic simulations of accretion along a magnetic field line coupled with radiative transfer models and a simple treatment of rotation to generate synthetic light curves. We adopted previously developed metrics in order to classify observations to facilitate comparisons between observations and our models. We found that stellar mass, magnetic field geometry, corotation radius, inclination, and turbulence all play roles in producing the observed light curves and that no single parameter is entirely dominant in controlling the observed variability. While the periodic behavior of the light curve is most strongly affected by the inclination, it is also a function of the magnetic field geometry and inner disk turbulence. Objects with either pure dipole fields, strong aligned octupole components, or high turbulence in the inner disk all tend to display accretion bursts. Objects with anti-aligned octupole components or aligned, weaker octupole components tend to show light curves with slightly fewer bursts. We did not find clear monotonic trends between the stellar mass and empirical classification. This work establishes the groundwork for more detailed characterization of well-studied targets as more light curves of CTTS become available through missions such as the Transiting Exoplanet Survey Satellite (TESS). [ABSTRACT FROM AUTHOR]

'Orion fingers' are a system of dozens of bow shocks, with the wings of shocks pointing to a common system of origin, which is centred on a dynamically disintegrating system of several massive stars. The shock heads propagate with velocities of up to 300–400 km s−1, but the formation and physical properties of the 'bullets' leading the shocks are not known. Here, we summarize two possible scenarios for the formation of the 'bullets' and the resulting bow shocks ('fingers'). In the first scenario, bullets are self-gravitating, Jupiter-mass objects that were formed rapidly and then ejected during the strong dynamical interactions of massive stars and their discs. This scenario naturally explains the similar time-scales for the outflow of bullets and for the dynamical interaction of the massive stars, but has some difficulty explaining the observed high velocities of the bullets. In the second scenario, bullets are formed via hydrodynamic instabilities in a massive, infrared-driven wind, naturally explaining the high velocities and the morphology of outflow, but the bullets are not required to be self-gravitating. The processes that created the Orion fingers are likely not unique to this particular star-forming region and may result in free-floating, high-velocity, core-less planets. [ABSTRACT FROM AUTHOR]