Search

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., Comment: 16 pages, 1 table, 14 figures; accepted for publication in MNRAS

Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets. While JWST Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a better understanding of their host star. Here, we propose a roadmap to characterize the TRAPPIST-1 system -- and others like it -- in an efficient and robust manner. We notably recommend that -- although more challenging to schedule -- multi-transit windows be prioritized to mitigate the effects of stellar activity and gather up to twice more transits per JWST hour spent. We conclude that, for such systems, planets cannot be studied in isolation by small programs, but rather need large-scale, jointly space- and ground-based initiatives to fully exploit the capabilities of JWST for the exploration of terrestrial planets.

Disc accretion rate onto low mass protostar FU Ori suddenly increased hundreds of times 85 years 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 exceeding 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 $\sim$ 10,000 years with radius of $\sim 10 R_J$. We show that their EE rates are $\sim 10^{-5}$ Msun/yr; 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 $\sim$ O(100) years. We find that a $\sim$ 6 Jupiter mass planet evaporating in a disc fed at a time-averaged rate of $\sim 10^{-6}$ Msun/yr 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., Comment: Accepted to MNRAS

Planetary magnetic fields can affect the predicted mass loss rate for close-in planets that experience large amounts of UV 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 10830 Angstrom line of helium. Motivated by previous work on hydrodynamic and magneto-hydrodynamic photoevaporation, we suggest that planets with magnetic fields that are too weak to control the outflow's topology lead to blue-shifted transits due to day-to-night-side 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 to test this concept, computing synthetic transit profiles in helium. As expected, we find that hydrodynamically dominated outflows lead to blue-shifted transits on the order of the sound speed of the gas. Strong surface magnetic fields lead to unshifted or slightly red-shifted transit profiles. High-resolution observations can distinguish between these profiles; however, eccentricity uncertainties generally mean that we cannot conclusively say velocity shifts are due to the outflow for individual planets. The majority of helium observations are blue-shifted, which could be a tentative indication that close-in planets generally have surface dipole magnetic field strengths $\lesssim 0.1$ gauss. More 3D hydrodynamic and magneto-hydrodynamic are needed to confirm this conclusion robustly., Comment: Published in MNRAS

The spin evolution of main sequence stars has long been of interest for basic stellar evolution, stellar aging, 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 which 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., Comment: 8 pages, 4 figures; Accepted for publication in MNRAS

The Galactic Centre is known to have undergone a recent star formation episode a few Myrs 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 sub-sonic 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 $\lesssim1~$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 $\sim$10 Myr in this state. Encounters with other stars are likely to remove the planet on Myr timescales, causing photoevaporation to restart, giving rise to a G-cloud signature. A giant planet fraction of $\sim10\%$ can explain the number of observed G-clouds., Comment: Accepted for publication in MNRAS

We used the Atacama Pathfinder Experiment (APEX) telescope to search for CI 1-0 (492.16GHz) emission towards 8 proplyds in NGC 1977, which is an FUV radiation environment two orders of magnitude weaker than that irradiating the Orion Nebular Cluster (ONC) proplyds. CI is expected to enable us to probe the wind launching region of externally photoevaporating discs. Of the 8 targets observed, no 3$\sigma$ detections of the CI line were made despite reaching sensitivities deeper than the anticipated requirement for detection from prior APEX CI observations of nearby discs and models of external photoevaporation of quite massive discs. By comparing both the proplyd mass loss rates and CI flux constraints with a large grid of external photoevaporation simulations, we determine that the non-detections are in fact fully consistent with the models if the proplyd discs are very low mass. Deeper observations in CI and probes of the disc mass with other tracers (e.g. in the continuum and CO) can test this. If such a test finds higher masses, this would imply carbon depletion in the outer disc, as has been proposed for other discs with surprisingly low CI fluxes, though more massive discs would also be incompatible with models that can explain the observed mass loss rates and CI non-detections. The expected remaining lifetimes of the proplyds are estimated to be similar to those of proplyds in the ONC at 0.1Myr. Rapid destruction of discs is therefore also a feature of common, intermediate UV environments., Comment: 10 pages. Accepted for publication in MNRAS

The system of two transiting Neptune-sized planets around the bright, young M-dwarf AU Mic provides a unique opportunity to test models of planet formation, early evolution, and star-planet interaction. However, the intense magnetic activity of the host star makes measuring the masses of the planets via the radial velocity (RV) method very challenging. We report on a 1-year, intensive monitoring campaign of the system using 91 observations with the HARPS spectrograph, allowing for detailed modelling of the $\sim 600\, \rm{m\,s^{-1}}$ peak-to-peak activity-induced RV variations. We used a multidimensional Gaussian Process framework to model these and the planetary signals simultaneously. We detect the latter with semi-amplitudes of $\rm{K_b} = 5.8 \pm 2.5\, \rm{m\,s^{-1}}$ and $\rm{K_c} = 8.5 \pm 2.5\, \rm{m\, s^{-1}}$, respectively. The resulting mass estimates, $\rm{M_b} = 11.7 \pm 5.0\, \rm{M_{\rm \oplus}}$ and $\rm{M_c} = 22.2 \pm 6.7\, \rm{M_{\rm \oplus}}$, suggest that planet b might be less dense, and planet c considerably denser than previously thought. These results are in tension with the current standard models of core-accretion. They suggest that both planets accreted a H/He envelope that is smaller than expected, and the trend between the two planets' envelope fractions is the opposite of what is predicted by theory.

In 2017, the California-Kepler Survey (CKS) published its first data release (DR1) of high-resolution optical spectra of 1305 planet hosts. Refined CKS planet radii revealed that small planets are bifurcated into two distinct populations: super-Earths (smaller than 1.5 $R_E$) and sub-Neptunes (between 2.0 and 4.0 $R_E$), with few planets in between (the "Radius Gap.") Several theoretical models of the Radius Gap predict variation with stellar mass, but testing these predictions are challenging with CKS DR1 due to its limited $M_\star$ range of 0.8-1.4 $M_\odot$. Here, we present CKS DR2 with 411 additional spectra and derived properties focusing on stars of 0.5-0.8 $M_\odot$. We found the Radius Gap follows $R_p \propto P^m$ with $m = -0.10 \pm 0.03$, consistent with predictions of XUV- and core-powered mass-loss mechanisms. We found no evidence that $m$ varies with $M_\star$. We observed a correlation between the average sub-Neptune size and $M_\star$. Over 0.5 to 1.4 $M_\odot$, the average sub-Neptune grows from 2.1 to 2.6 $R_E$, following $R_p \propto M_\star^\alpha$ with $\alpha = 0.25 \pm 0.03$. In contrast, there is no detectable change for super-Earths. These $M_\star$-$R_p$ trends suggests that protoplanetary disks can efficiently produce cores up to a threshold mass of $M_c$, which grows linearly with stellar mass according to $M_c \approx 10 M_E~(M_\star / M_\odot)$. There is no significant correlation between sub-Neptune size and stellar metallicity (over $-$0.5 to $+$0.5 dex) suggesting a weak relationship between planet envelope opacity and stellar metallicity. Finally, there is no significant variation in sub-Neptune size with stellar age (over 1 to 10 Gyr), which suggests that the majority of envelope contraction concludes after $\sim$1 Gyr., Comment: 30 pages, 16 figures, 4 tables; accepted for publication in The Astronomical Journal; machine readable versions of Tables 1 and 2 uploaded

Mature super-Earths and sub-Neptunes are predicted to be $\simeq$Jovian radius when younger than 10 Myr. Thus, we expect to find 5-15$R_\oplus$ planets around young stars even if their older counterparts harbor none. We report the discovery and validation of TOI 1227 b, a $0.85\pm0.05R_J$ (9.5$R_\oplus$) planet transiting a very low-mass star ($0.170\pm0.015M_\odot$) every 27.4 days. TOI~1227's kinematics and strong lithium absorption confirm it is a member of a previously discovered sub-group in the Lower Centaurus Crux OB association, which we designate the Musca group. We derive an age of 11$\pm$2 Myr for Musca, based on lithium, rotation, and the color-magnitude diagram of Musca members. The TESS data and ground-based follow-up show a deep (2.5\%) transit. We use multiwavelength transit observations and radial velocities from the IGRINS spectrograph to validate the signal as planetary in nature, and we obtain an upper limit on the planet mass of $\simeq0.5 M_J$. Because such large planets are exceptionally rare around mature low-mass stars, we suggest that TOI 1227 b is still contracting and will eventually turn into one of the more common $<5R_\oplus$ planets., Comment: Accepted to the Astronomical Journal. Minor updates during referee process and proofs

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 UV 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 evolves in their outer regions and the observable imprints of the external photoevaporation process., Comment: Accepted for publication in MNRAS

Close-in super-Earths are the most abundant exoplanets known. It has been hypothesized that they form in the inner regions of protoplanetary discs, out of the dust that may accumulate at the boundary between the inner region susceptible to the magneto-rotational instability (MRI) and an MRI-dead zone further out. In Paper I we presented a model for the viscous inner disc which includes heating due to both irradiation and MRI-driven accretion; thermal and non-thermal ionization; dust opacities; and dust effects on ionization. Here we examine how the inner disc structure varies with stellar, disc and dust parameters. For high accretion rates and small dust grains, we find that: (1) the main sources of ionization are thermal ionization and thermionic and ion emission; (2) the disc features a hot, high-viscosity inner region, and a local gas pressure maximum at the outer edge of this region (in line with previous studies); and (3) an increase in the dust-to-gas ratio pushes the pressure maximum outwards. Consequently, dust can accumulate in such inner discs without suppressing the MRI, with the amount of accumulation depending on the viscosity in the MRI-dead regions. Conversely, for low accretion rates and large dust grains, there appears to be an additional steady-state solution in which: (1) stellar X-rays become the main source of ionization; (2) MRI-viscosity is high throughout the disc; and (3) the pressure maximum ceases to exist. Hence, if planets form in the inner disc, larger accretion rates (and thus younger disks) are favoured., Comment: accepted to MNRAS

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 magneto-rotational 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$_\odot$yr$^{-1}$) and small dust grains, we find that the inner disc is optically thick, and the accretion heat is primarily released near the midplane. As a result, both the disc midplane 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 onto 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., Comment: accepted for publication in MNRAS

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 $\dot{M}$. 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)., Comment: 18 pages, 11 figures

Volatile species in protoplanetary discs can undergo a phase change from vapour to solid. These "snow-lines" can play vital roles in planet formation at all scales, from dust coagulation to planetary migration. In the outer regions of protoplanetary discs, the temperature profile is set by the absorption of reprocessed stellar light by the solids. Further, the temperature profile sets the distribution of solids through sublimation and condensation at various snow-lines. Hence the snow-line position depends on the temperature profile and vice-versa. We show that this coupling can be thermally unstable, such that a patch of the disc at a snow-line will produce either run-away sublimation or condensation. This thermal instability arises at moderate optical depths, where heating by absorption of reprocessed stellar light from the disc's atmosphere is optically thick, yet cooling is optically thin. Since volatiles in the solid phase drift much faster than volatiles in the vapour phase, this thermal instability results in a limit-cycle. The snow-line progressively moves in, condensing volatiles, before receding, as the volatiles sublimate. Using numerical simulations, we study the evolution of the CO snow-line. We find the CO snow-line is thermally unstable under typical disc conditions and evolves inwards from $\sim50$ to $\sim30$~AU on timescales from 1,000-10,000 years. The CO snow-line spends between $\sim 10\%-50\%$ of its time at smaller separations, where the exact value is sensitive to the total opacity and turbulent viscosity. The evolving snow-line also creates ring-like structures in the solid distribution interior to the snow-line. Multiple ring-like structures created by moving snow-lines could potentially explain the sub-structures seen in many {\it ALMA} images., Comment: Accepted for publication in MNRAS

The Sun shows a $\sim 10$% 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 $\gtrsim 100 $M$_\oplus$ of dust exterior to its orbit, preventing the dust from accreting onto the star in contrast to the gas. Forming giant planets can create refractory depletions of $\sim 5-15$%, 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 onto the star., Comment: 11 pages, 5 figures. Accepted for publication in MNRAS

We use a suite of SPH simulations to investigate the susceptibility of protoplanetary discs to the effects of self-gravity as a function of star-disc properties. We also include passive irradiation from the host star using different models for the stellar luminosities. The critical disc-to-star mass ratio for axisymmetry (for which we produce criteria) increases significantly for low-mass stars. This could have important consequences for increasing the potential mass reservoir in a proto Trappist-1 system, since even the efficient Ormel et al. (2017) formation model will be influenced by processes like external photoevaporation, which can rapidly and dramatically deplete the dust reservoir. The aforementioned scaling of the critical $M_d/M_*$ for axisymmetry occurs in part because the Toomre $Q$ parameter has a linear dependence on surface density (which promotes instability) and only an $M_*^{1/2}$ dependence on shear (which reduces instability), but also occurs because, for a given $M_d/M_*$, the thermal evolution depends on the host star mass. The early phase stellar irradiation of the disc (for which the luminosity is much higher than at the zero age main sequence, particularly at low stellar masses) can also play a key role in significantly reducing the role of self-gravity, meaning that even Solar mass stars could support axisymmetric discs a factor two higher in mass than usually considered possible. We apply our criteria to the DSHARP discs with spirals, finding that self-gravity can explain the observed spirals so long as the discs are optically thick to the host star irradiation., Comment: 21 pages, Accepted for publication MNRAS

We demonstrate the utility of CI 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, CI traces unsubtle kinematic signatures of a wind in intermediate UV environments ($\sim1000$G$_0$) and can be used to place constraints on the kinematics and temperature of the wind. In CI position-velocity diagrams external photoevaporation results in velocities that are faster than those from Keplerian rotation alone, as well as emission from quadrants of position-velocity 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 ($\sim 30-60^\circ$) 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). CI is readily observable with ALMA, which means that it is now possible to identify and characterise the effect of external photoevaporation on planet-forming discs in intermediate UV environments., Comment: 12 pages, accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS

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 ~1e5 years, while the photoevaporation created cavity still resides at 10s 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., Comment: 14 pages, accpeted for publication in MNRAS

Photoevaporation of planet forming discs by high energy radiation from the central star is potentially a crucial mechanism for disc evolution and it may play an important role in the formation and evolution of planetary system. We present here a new generation of X-ray photoevaporation models for solar-type stars, based on a new set of hydrodynamical simulations, which account for stellar irradiation via a new, significantly improved, parameterisation of gas temperatures, based on detailed photoionisation and radiation transfer calculations. This is the first of a series of papers aiming at providing a library of models which cover the observed parameter space in stellar and disc mass, metallicity and stellar X-ray properties. We focus here on solar-type stars (0.7 M$_\odot$) with relatively low-mass discs (1% of the stellar mass) and explore the dependence of the wind mass loss rates on stellar X-ray luminosity. We model primordial discs as well as transition discs at various stages of evolution. Our 2D hydrodynamical models are then used to derive simple recipes for the mass loss rates that are suitable for inclusion in one-dimensional disc evolution and/or planet formation models typically employed for population synthesis studies. Line profiles from typical wind diagnostics ([OI] 6300 $\overset{\lower.5em\circ}{\mathrm{A}}$ and [NeII] 12.8 $\mu$m) are also calculated for our models and found to be roughly in agreement with previous studies. Finally, we perform a population study of transition discs by means of one-dimensional viscous evolution models including our new photoevaporation prescription and find that roughly a half of observed transition discs cavities and accretion rates could be reproduced by our models., Comment: 12 pages, 14 figures, published in MNRAS

It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current instrumentation. In line with theCommittee on an Exoplanet Science Strategy recognition that holistic observational approaches are needed, we point out here that a full understanding of exoplanet atmospheres, their evolution and determination of habitability requires a powerful high-resolution X-ray imaging and spectroscopic observatory. This is the only capability that can: (1) characterize by proxy the crucial, difficult to observe, EUV stellar flux, its history and its variations for planet hosting stars; (2) observe the stellar wind; (3) detect the subtle Doppler signatures of coronal mass ejections., Comment: Astro2020 Decadal Survey Science White Paper

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 which 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 onto 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 onto 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., Comment: Accepted for publication in MNRAS

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-6Rearth) 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 days. 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 sigma 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., Comment: To appear in MNRAS

The large population of Earth to super-Earth sized planets found very close to their host stars has motivated consideration of $in$ $situ$ formation models. In particular, Inside-Out Planet Formation is a scenario in which planets coalesce sequentially in the disk, at the local gas pressure maximum near the inner boundary of the dead zone. The pressure maximum arises from a decline in viscosity, going from the active innermost disk (where thermal ionization of alkalis yields high viscosities via the magneto-rotational instability (MRI)) to the adjacent dead zone (where the MRI is quenched). Previous studies of the pressure maximum, based on $\alpha$-disk models, have assumed ad hoc values for the viscosity parameter $\alpha$ in the active zone, ignoring the detailed physics of the MRI. Here we explicitly couple the MRI criteria to the $\alpha$-disk equations, to find steady-state (constant accretion rate) solutions for the disk structure. We consider the effects of both Ohmic and ambipolar resistivities, and find solutions for a range of disk accretion rates ($\dot{M}$ = $10^{-10}$ - $10^{-8}$ ${\rm M}_{\odot}$/yr), stellar masses ($M_{\ast}$ = 0.1 - 1 ${\rm M}_{\odot}$), and fiducial values of the $non$-MRI $\alpha$-viscosity in the dead zone ($\alpha_{\rm {DZ}} = 10^{-5}$ - $10^{-3}$). We find that: (1) A midplane pressure maximum forms radially $outside$ the inner boundary of the dead zone; (2) Hall resistivity dominates near the midplane in the inner disk, which may explain why close-in planets do $not$ form in $\sim$50% of systems; (3) X-ray ionization can be competitive with thermal ionization in the inner disk, because of the low surface density there in steady-state; and (4) our inner disk solutions are viscously unstable to surface density perturbations., Comment: 34 pages, 28 figures, 3 appendices. Accepted by the Astrophysical Journal

Many protostellar gapped and binary discs show misalignments between their inner and outer discs; in some cases, $\sim70$ degree misalignments have been observed. Here we show that these misalignments can be generated through a "secular precession resonance" between the nodal precession of the inner disc and the precession of the gap-opening (stellar or massive planetary) companion. An evolving protostellar system may naturally cross this resonance during its lifetime due to disc dissipation and/or companion migration. If resonance crossing occurs on the right timescale, of order a few Myrs, characteristic for young protostellar systems, the inner and outer discs can become highly misaligned ($\gtrsim 60$ degrees). When the primary star has a mass of order a solar mass, generating a significant misalignment typically requires the companion to have a mass of $\sim 0.01-0.1$ M$_\odot$ and an orbital separation of tens of AU. The recently observed companion in the cavity of the gapped, highly misaligned system HD 142527 satisfies these requirements, indicating that a previous resonance crossing event misaligned the inner and outer discs. Our scenario for HD 142527's misaligned discs predicts that the companion's orbital plane is aligned with the outer disc's; this prediction should be testable with future observations as the companion's orbit is mapped out. Misalignments observed in several other gapped disc systems could be generated by the same secular resonance mechanism., Comment: 12 pages, 8 figures, submitted to MNRAS

We present a simple model for low-mass planet formation and subsequent evolution within "transition" discs. We demonstrate quantitatively that the predicted and observed structure of such discs are prime birthsites of planets. Planet formation is likely to proceed through pebble accretion, should a planetary embryo ($M\gtrsim 10^{-4}\,$M$_\oplus$) form. Efficient pebble accretion is likely to be unavoidable in transition disc dust traps, as the size of 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 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 $\sim 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 ~10AU separation should be a common outcome of this cycle., Comment: 15 pages, accepted for publication in MNRAS

Massive giant planets, such as the ones being discovered by direct imaging surveys, likely experience the majority of their growth through a circumplanetary disc. We argue that the entropy of accreted material is determined by boundary layer processes, unlike the "cold-" or "hot-start" hypotheses usually invoked in the core accretion and direct collapse scenarios. A simple planetary evolution model illustrates how a wide range of radius and luminosity tracks become possible, depending on details of the accretion process. Specifically, the proto-planet evolves towards "hot-start" tracks if the scale-height of the boundary layer is $\gtrsim0.24$, a value not much larger than the scale-height of the circumplanetary disc. Understanding the luminosity and radii of young giant planets will thus require detailed models of circumplanetary accretion., Comment: 7 pages, 4 figures, to appear in ApJL

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 characterising planets in the "Habitable Zone" (HZ). However, Kepler finds that terrestrial-mass exoplanets are often born with voluminous H/He envelopes, comprising mass-fractions ($M_{env}/M_{core}$) $\gtrsim 1$%. If these planets retain such envelopes over Gyr timescales, 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 $\lesssim 0.9$ M$_\oplus$, can lose enough H/He to become habitable if their initial envelope mass-fraction is $\sim$1%, (3) at the outer edge of the HZ, evaporation cannot remove a $\sim$1% H/He envelope even from cores down to 0.8 M$_\oplus$. Thus, if planets form with bulky H/He envelopes, only those with low-mass cores may eventually be habitable. Cores $\gtrsim$ 1 M$_\oplus$, with $\gtrsim$1% natal H/He envelopes, will not be habitable in the HZ of M dwarfs., Comment: 52 pages, 21 figures, accepted for publication in MNRAS

Transition discs are protoplanetary discs that show evidence for large holes or wide gaps (with widths comparable to their radii) in their dust component. These discs could be giving us clues about the disc destruction mechanism or hints about the location and time-scales for the formation of planets. However, at the moment there remain key gaps in our theoretical understanding. The vast majority of transition discs are accreting onto their central stars, indicating that - at least close to the star - dust has been depleted from the gas by a very large amount. In this review, we discuss evidence for two distinct populations of transition discs: mm-faint - those with low mm-fluxes, small holes ($\lesssim$10AU) and low accretion rates (~1e-10-1e-9 Msun/yr) and mm-bright - discs with large mm-fluxes, large holes ($\gtrsim$20 AU) and high accretion rates ~1e-8 Msun/yr. MM-faint transition discs are consistent with what would naively be expected from a disc undergoing dispersal; however, mm-bright discs are not, and are likely to be rare and long-lived objects. We discuss the two commonly proposed mechanisms for creating transition discs: photoevaporation and planet-disc interactions, with a particular emphasis on how they would evolve in these models, comparing these predictions to the observed population. More theoretical work on explaining the lack of optically thick, non-accreting transition discs is required in both the photoevaporation and planetary hypothesis, before we can start to use transition discs to constrain models of planet formation. Finally, we suggest that the few discs with primordial looking SEDs, but serendipitously imaged, showing large cavities in the mm (e.g. MWC758 & WSB 60) may represent a hidden population of associated objects. Characterising and understanding how these objects fit into the overall paradigm may allow us to unravel the mystery of transition discs., Comment: 24 pages. Accepted for publication in PASA as an invited, refereed, review as part of the special issue on "Disc dynamics and planet formation"

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 inward 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)., Comment: 16 pages, accepted for publication in MNRAS

The lack of observed transition discs with inner gas holes of radii greater than ~50AU implies that protoplanetary discs dispersed from the inside out must remove gas from the outer regions rapidly. We investigate the role of photoevaporation in the final clearing of gas from low mass discs with inner holes. In particular, we study the so-called "thermal sweeping" mechanism which results in rapid clearing of the disc. Thermal sweeping was originally thought to arise when the radial and vertical pressure scale lengths at the X-ray heated inner edge of the disc match. We demonstrate that this criterion is not fundamental. Rather, thermal sweeping occurs when the pressure maximum at the inner edge of the dust heated disc falls below the maximum possible pressure of X-ray heated gas (which depends on the local X-ray flux). We derive new critical peak volume and surface density estimates for rapid radiative clearing which, in general, result in rapid dispersal happening less readily than in previous estimates. This less efficient clearing of discs by X-ray driven thermal sweeping leaves open the issue of what mechanism can clear gas from the outer disc sufficiently quickly to explain the non-detection of cold gas around weak line T Tauri stars., Comment: 13 pages, Accepted for publication in MNRAS

The Kepler planetary system consists of two exoplanets at similar separations (0.115 & 0.128 AU), which have dramatically different densities. The inner planet has a density consistent with an Earth-like composition, while the outer planet is extremely low-density, such that it must contain a voluminous H/He envelope. Such a density difference would pose a problem for any formation mechanism if their current densities were representative of their composition at formation. However, both planets are at close enough separations to have undergone significant evaporation in the past. We constrain the core-mass, core composition, initial envelope-mass, and initial cooling-time of each planet using evaporation models conditioned on their present-day masses and radii, as inferred from Kepler photometry and transit timing analysis. The inner planet is consistent with being an evaporatively stripped core, while the outer planet has retained some of its initial envelope due to its higher core-mass. Therefore, both planets could have had a similar formation pathway, with the inner planet having an initial envelope-mass fraction of $\lesssim 10\%$ and core-mass of $\sim4.4$ M$_\oplus$, while the outer had an initial envelope-mass fraction of order 15-30\% and core-mass $\sim7.3$ M$_\oplus$. Finally, our results indicate that the outer planet had a long ($\gtrsim30$ Myr) initial cooling-time, much longer than would naively be predicted from simple timescale arguments. The long initial cooling-time could be evidence for a dramatic early cooling episode such as the recently proposed "boil-off" process., Comment: Accepted version in ApJL, 6 pages, 3 figures

Late-type main sequence stars exhibit an x-ray to bolometric flux ratio that depends on ${\tilde Ro}$, the ratio of rotation period to convective turnover time, as ${\tilde Ro}^{-\zeta}$ with $2\le \zeta \le 3$ for ${\tilde Ro} >0.13$, but saturates with $|\zeta| <0.2$ for ${\tilde Ro} < 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 time-evolution 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 x-ray 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 toward a unified explanation of all of the aforementioned observational trends., Comment: 12 pages, 5 figures; published in MNRAS; this version updated to correct some typographical errors identified in a post-publication erratum

Chondrule formation remains one of the most elusive early Solar System events. Here, we take the novel approach of employing numerical simulations to investigate chondrule origin beyond purely cosmochemical methods. We model the transport of generically-produced chondrules and dust in a 1D viscous protoplanetary disk model, in order to constrain the chondrule formation events. For a single formation event we are able to match analytical predictions of the memory chondrule and dust populations retain of each other (complementarity), finding that a large mass accretion rate ($\gtrsim 10^{-7}$~M$_\odot$~yr$^{-1}$) allows for delays on the order of the disk's viscous timescale between chondrule formation and chondrite accretion. Further, we find older disks to be severely diminished of chondrules, with accretion rates $\lesssim 10^{-9}$~M$_\odot$~yr$^{-1}$ for nominal parameters. We then characterize the distribution of chondrule origins in both space and time, as functions of disk parameters and chondrule formation rates, in runs with continuous chondrule formation and both static and evolving disks. Our data suggest that these can account for the observed diversity between distinct chondrite classes, if some diversity in accretion time is allowed for., Comment: Accepted to MNRAS. 17 pages, 15 figures, 1 table

We show that, for a low-mass planet that orbits its host star within a few tenths of an AU (like the majority of the {\it Kepler} planets), the atmosphere it was able to accumulate while embedded in the proto-planetary disk may not survive unscathed after the disk disperses. This gas envelope, if more massive than a few percent of the core (with a mass below $10 M_\oplus$), has a cooling time that is much longer than the time-scale on which the planet exits the disk. As such, it could not have contracted significantly from its original size, of order the Bondi radius. So a newly exposed proto-planet would be losing mass via a Parker wind that is catalyzed by the stellar continuum radiation. This represents an intermediate stage of mass-loss, occurring soon after the disc has dispersed, but before the EUV/X-ray driven photoevaporation becomes relevant. The surface mass-loss induces a mass movement within the envelope that advects internal heat outward. As a result, the planet atmosphere rapidly cools down and contracts, until it has reached a radius of order $0.1$ Bondi radius, at which time the mass-loss effectively shuts down. Within a million years after the disk disperses, we find a planet that has only about ten percent of its original envelope, and a Kelvin-Helmholtz time that is much longer than its actual age. We suggest that this "boil-off" process may be partially responsible for the lack of planets above a radius of $2.5 R_\oplus$ in the {\it Kepler} data, provided planet formation results in initial envelope masses of tens of percent., Comment: Published in ApJ

We have investigated the evaporation of close-in exoplanets irradiated by ionizing photons. We find that the properties of the flow are controlled by the ratio of the recombination time to the flow time-scale. When the recombination time-scale is short compared to the flow time-scale the the flow is in approximate local ionization equilibrium with a thin ionization front, where the photon mean free path is short compared to flow scale. In this "recombination limited" flow the mass-loss scales roughly with the square root of the incident flux. When the recombination time is long compared to the flow time-scale the ionization front becomes thick and encompasses the entire flow, with the mass-loss rate scaling linearly with flux. If the planet's potential is deep the flow is approximately "energy-limited"; however, if the planet's potential is shallow we identify a new limiting mass-loss regime, which we term "photon-limited". In this scenario the mass-loss rate is purely limited by the incoming flux of ionizing photons. We have developed a new numerical approach that takes into account the frequency dependence of the incoming ionizing spectrum and performed a large suite of 1D simulations to characterise UV driven mass-loss around low mass planets. We find the flow is "recombination-limited" at high fluxes but becomes "energy-limited" at low fluxes; however, the transition is broad occurring over several order of magnitude in flux. Finally, we point out the transitions between the different flow types does not occur at a single flux value, but depends on the planet's properties, with higher mass planets becoming "energy-limited" at lower fluxes., Comment: Published in ApJ

We show that if young low mass stars undergo vigorous X-ray driven disc winds, these may be detected in clusters through their interaction with ionising radiation from massive stars. We argue that in the ONC (Orion Nebula Cluster) one should see $\sim$ 10s of `X-ray proplyds' ( objects with optically imaged offset ionisation fronts) in the range $0.3-0.6$pc from $\theta_1$C Ori (the dominant O star in the ONC). Such objects lie outside the central `FUV zone' in the ONC where proplyds are instead well explained by neutral winds driven by external Far Ultraviolet (FUV) emission from $\theta_1$C. We show that the predicted numbers and sizes of X-ray proplyds are compatible with those observed and that this may also explain at least some of the far flung proplyds seen 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 ago, $\sim$ the flow time to the ionisation front). We discuss whether variability on this timescale can plausibly explain proplyd size data on a case by case basis. We also show that the predicted radio free-free emission signature of X-ray proplyds is readily detectable. Monitoring is however required to distinguish such emission from non-thermal emission from active coronae. We also predict that it is only at $> 1$ pc from $\theta_1$C that the free-free emission from such offset ionised structures would be clearly distinguishable from an externally driven ionised 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 is a promising avenue for observationally testing X-ray photoevaporation., Comment: 8 pages, 3 figures, MNRAS accepted

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)$_{\rm water}$ in 1D simulations of a disc's evolution that include dead zones and infall from an envelope with given D/H ratio in the infalling water ($\sim 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)$_{\rm water}$ ratio and a more extended region (radius up to $\sim 1-3$ AU) where water is equilibrated with hydrogen gas (D/H=$2\times 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 $\lesssim 0.2$ or had a Schmidt number closer to "nominal" values ($\sim 1$) and experienced several accretion bursts. Finally, future observations of (D/H)$_{\rm water}$ in protoplanetary discs will allow inferences about angular momentum properties of the disc during disc building and the role of accretion bursts., Comment: 14 pages, 13 figures, accepted for publication in MNRAS

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 gravo-magnetic 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 ($\lesssim 25$) we find that turbulent diffusion dominates the radial transport of heat, reducing the burst duration to $\lesssim 10^3$ years 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$_\odot$ yr$^{-1}$ with no diffusion, to $10^{-8}$ to $\gtrsim10^{-4}$ M$_\odot$ 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 gravo-magnetic limit cycle., Comment: 11 Pages, 11 Figures, Accepted for publication in MNRAS

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 $B_P\gtrsim 0.3$ gauss and fluxes $F_{UV}\sim10^{6}$ 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 steady-state, resulting in time-variable mass-loss., Comment: Accepted for publication in MNRAS, 20 pages, 13 figures

Sharp chemical discontinuities can occur in protoplanetary discs, particularly at `snow-lines' where a gas-phase species freezes out to form ice grains. Such sharp discontinuities will diffuse out due to the turbulence suspected to drive angular momentum transport in accretion discs. We demonstrate that the concentration gradient - in the vicinity of the snow-line - of a species present outside a snow-line but destroyed inside is strongly sensitive to the level of turbulent diffusion (provided the chemical and transport time-scales are decoupled) and provides a direct measurement of the radial `Schmidt number' (the ratio of the angular momentum transport to radial turbulent diffusion). Taking as an example the tracer species N$_2$H$^+$, which is expected to be destroyed inside the CO snow-line (as recently observed in TW Hya) we show that ALMA observations possess significant angular resolution to constrain the Schmidt number. Since different turbulent driving mechanisms predict different Schmidt numbers, a direct measurement of the Schmidt number in accretion discs would allow inferences about the nature of the turbulence to be made., Comment: 6 Pages, 2 Figures, Accepted for publication in ApJL

(Abridged) We consider the properties of protoplanetary discs that are undergoing inside-out clearing by photoevaporation. In particular, we aim to characterise the conditions under which a protoplanetary disc may undergo `thermal sweeping', a rapid (< 1e4 years) disc destruction mechanism proposed to occur when a clearing disc reaches sufficiently low surface density at its inner edge and where the disc is unstable to runaway penetration by the X-rays. We use a large suite of 1D radiation-hydrodynamic simulations to probe the observable parameter space, which is unfeasible in higher dimensions. These models allow us to determine the surface density at which thermal sweeping will take over the disc's evolution and to evaluate this critical surface density as a function of X-ray luminosity, stellar mass and inner hole radius. We find that this critical surface density scales linearly with X-ray luminosity, increases with inner hole radius and decreases with stellar mass and we develop an analytic model that reproduces these results. This surface density criterion is then used to determine the evolutionary state of protoplanetary discs at the point that they become unstable to destruction by thermal sweeping. We find that transition discs created by photoevaporation will undergo thermal sweeping when their inner holes reach 20-40 AU, implying that transition discs with large holes and no accretion (which were previously a predicted outcome of the later stages of all flavours of photoevaporation model) will not form. We emphasise that the surface density criteria that we have developed apply to all situations where the disc develops an inner hole that is optically thin to X-rays. It thus applies not only to the case of holes originally created by photoevaporation but also to holes formed, for example, by the tidal influence of planets., Comment: 10 Pages, MNRAS accp

We consider continuum free-free radio emission from the upper atmosphere of protoplanetary discs as a probe of the ionized luminosity impinging upon the disc. Making use of previously computed hydrodynamic models of disc photoevaporation within the framework of EUV and X-ray irradiation, we use radiative transfer post-processing techniques to predict the expected free-free emission from protoplanetary discs. In general, the free-free luminosity scales roughly linearly with ionizing luminosity in both EUV and X-ray driven scenarios, where the emission dominates over the dust tail of the disc and is partial optically thin at cm wavelengths. We perform a test observation of GM Aur at 14-18 Ghz and detect an excess of radio emission above the dust tail to a very high level of confidence. The observed flux density and spectral index are consistent with free-free emission from the ionized disc in either the EUV or X-ray driven scenario. Finally, we suggest a possible route to testing the EUV and X-ray driven dispersal model of protoplanetary discs, by combining observed free-free flux densities with measurements of mass-accretion rates. On the point of disc dispersal one would expect to find a M_dot^2 scaling with free-free flux in the case of EUV driven disc dispersal or a M_dot scaling in the case of X-ray driven disc dispersal., Comment: Accepted MNRAS, 12 pages, 11 figures, (pdf generation fixed)

We assess the potential of planet formation instigating the early formation of a photoevaporation driven gap, up to radii larger than typical for photoevaporation alone. For our investigation we make use of hydrodynamics models of photoevaporating discs with a giant planet embedded. We find that, by reducing the mass accretion flow onto the star, discs that form giant planets will be dispersed at earlier times than discs without planets by X-ray photoevaporation. By clearing the portion of the disc inner of the planet orbital radius, planet formation induced photoevaporation (PIPE) is able to produce transition disc that for a given mass accretion rate have larger holes when compared to standard X-ray photoevaporation. This constitutes a possible route for the formation of the observed class of accreting transition discs with large holes, which are otherwise difficult to explain by planet formation or photoevaporation alone. Moreover, assuming that a planet is able to filter dust completely, PIPE produces a transition disc with a large hole and may provide a mechanism to quickly shut down accretion. This process appears to be too slow however to explain the observed desert in the population of transition disc with large holes and low mass accretion rates., Comment: 11 pages, 10 figures, accepted by MNRAS on 31/12/2012

We perform a comparison between two radiative transfer algorithms commonly employed in hydrodynamical calculations of star formation: grey flux limited diffusion and the hybrid scheme, in addition we compare these algorithms to results from the Monte-Carlo radiative transfer code MOCASSIN. In disc like density structures the hybrid scheme performs significantly better than the FLD method in the optically thin regions, with comparable results in optically thick regions. In the case of a forming high mass star we find the FLD method significantly underestimates the radiation pressure by a factor of ~100., Comment: 4 Pages; to appear in the proceedings of 'The Labyrinth of Star Formation', Crete, 18-22 June 2012

We discuss a hydrodynamical model for the dispersal of protoplanetary discs around young, low mass (<1.5 M_sun) stars by photoevaporation from the central object's energetic radiation, which considers the far-ultraviolet 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 suggest 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 10e-11 to 10e-7 M_sun/yr 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 (L_FUV/L_X>100) luminosity ratios, 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., Comment: Accepted MNRAS, 25 pages

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 NIR 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 a 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, a 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., Comment: 8 pages, 6 figures, accepted for publication in MNRAS

(Abridged) 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 (10e28-10e31 erg\s) gives rise to vigorous disc winds with rates of order 10e-10-10e-7 M_sun/yr. 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 alpha = 2.5e-3. 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. Finally, we confirm the conjecture of Drake et al. (2009), that accretion is suppressed by the X-rays through `photoevaporation starved accretion' and predict this effect can give rise to a negative correlation between X-ray luminosity and accretion rate, as reported in the Orion data., Comment: Figure 12 and 13 have been updated. In the original version the results from an unused model run were plotted by mistake