TRAPPIST-1 is a fantastic nearby (~39.14 light years) planetary system made of at least seven transiting terrestrial-size, terrestrial-mass planets all receiving a moderate amount of irradiation. To date, this is the most observationally favourable system of potentially habitable planets. Since the announcement of the discovery of TRAPPIST-1 planets in 2016, a growing number of techniques and approaches have been used and proposed to reveal its true nature. Here we have compiled a state-of-the-art overview of all the observational and theoretical constraints that have been obtained so far using these techniques and approaches. The goal is to get a better understanding of whether or not TRAPPIST-1 planets can have atmospheres, and if so, what they are made of. For this, we surveyed the literature on TRAPPIST-1 about topics as broad as irradiation environment, orbital architecture, transit observations, density measurements, stellar contamination, and numerical climate and escape models. Each of these topics adds a brick to our understanding of the likely atmospheres of the seven planets. We show that (i) HST transit observations, (ii) density measurements, (iii) atmospheric escape modelling, and (iv) gas accretion modelling altogether offer solid evidence against the presence of H2-dominated atmospheres around TRAPPIST-1 planets. This means they likely have either (i) a high molecular weight atmosphere or (ii) no atmosphere at all. There are several key challenges ahead to characterize the bulk compositions of the atmospheres (if present) of TRAPPIST-1 planets. The main one so far is characterizing and correcting for the effects of stellar contamination. Fortunately, a new wave of observations with the James Webb Space Telescope and near-infrared high-resolution ground-based spectrographs on very large telescopes will bring significant advances in the coming decade.

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.

"Orion fingers" are a system of dozens of bowshocks, with the wings of shocks pointing to a common system of origin, which is centered on a dynamically disintegrating system of several massive stars. The shock heads propagate with velocities of up to 300-400 km/s, 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 bowshocks ("fingers"). In the first scenario, bullets are self-gravitating, Jupiter-mass objects which were formed rapidly and then ejected during the strong dynamical interactions of massive stars and their disks. This scenario naturally explains the similar timescales 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.

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.

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.

The upcoming launch of the James Webb Space Telescope (JWST) combined with the unique features of the TRAPPIST-1 planetary system should enable the young field of exoplanetology to enter into the realm of temperate Earth-sized worlds. Indeed, the proximity of the system (12pc) and the small size (0.12 Rsun) and luminosity (0.05 Lsun) of its host star should make the comparative atmospheric characterization of its seven transiting planets within reach of an ambitious JWST program. Given the limited lifetime of JWST, the ecliptic location of the star that limits its visibility to 100d per year, the large number of observational time required by this study, and the numerous observational and theoretical challenges awaiting it, its full success will critically depend on a large level of coordination between the involved teams and on the support of a large community. In this context, we present here a community initiative aiming to develop a well-defined sequential structure for the study of the system with JWST and to coordinate on every aspect of its preparation and implementation, both on the observational (e.g. study of the instrumental limitations, data analysis techniques, complementary space-based and ground-based observations) and theoretical levels (e.g. model developments and comparison, retrieval techniques, inferences). Depending on the outcome of the first phase of JWST observations of the planets, this initiative could become the seed of a major JWST Legacy Program devoted to the study of TRAPPIST-1.

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.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

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 multi-transiting 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 & 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.

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.

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.

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.

Exoplanets with substantial Hydrogen/Helium atmospheres have been discovered in abundance, many residing extremely close to their parent stars. The extreme irradiation levels these atmospheres experience causes 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 lower-mass super-Earths/mini-Neptunes. Detailed hydrodynamic calculations that incorporate radiative transfer and ionization chemistry are now common in one-dimensional models, and multi-dimensional calculations that incorporate magnetic-fields and interactions with the interstellar environment are cutting edge. However, there remains very limited comparison between simulations and observations. 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 we observe today.

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.

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.

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.

Various theoretical models treating the effect of stellar irradiation on planetary envelopes predict the presence of a radius valley: i.e. a bimodal distribution of planet radii, with super-Earths and sub-Neptune planets separated by a valley at around $\approx 2~R_\oplus$. Such a valley was observed recently, owing to an improvement in the precision of stellar, and therefore planetary radii. Here we investigate the presence, location and shape of such a valley using a small sample with highly accurate stellar parameters determined from asteroseismology, which includes 117 planets with a median uncertainty on the radius of 3.3%. We detect a clear bimodal distribution, with super-Earths ($\approx 1.5~R_\oplus$) and sub-Neptunes ($\approx 2.5~R_\oplus$) separated by a deficiency around $2~R_\oplus$. We furthermore characterize the slope of the valley as a power law $R \propto P^\gamma$ with $\gamma = {-0.09^{+0.02}_{-0.04}}$. A negative slope is consistent with models of photo-evaporation, but not with the late formation of rocky planets in a gas-poor environment, which would lead to a slope of opposite sign. The exact location of the gap further points to planet cores consisting of a significant fraction of rocky material.

The mass-period or radius-period distribution of close-in exoplanets shows a paucity of intermediate mass/size (sub-Jovian) planets with periods ~< 3 days. 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 result from the fact that photoevaporation is more effective closer to the host star, and that in order for a gas giant to tidally circularise 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.

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.

Seven temperate Earth-sized exoplanets readily amenable for atmospheric studies transit the nearby ultracool dwarf star TRAPPIST-1 (refs 1,2). Their atmospheric regime is unknown and could range from extended primordial hydrogen-dominated to depleted atmospheres (refs 3-6). Hydrogen in particular is a powerful greenhouse gas that may prevent the habitability of inner planets while enabling the habitability of outer ones (refs 6-8). An atmosphere largely dominated by hydrogen, if cloud-free, should yield prominent spectroscopic signatures in the near-infrared detectable during transits. Observations of the innermost planets have ruled out such signatures (ref 9). However, the outermost planets are more likely to have sustained such a Neptune-like atmosphere (refs 10,11). Here, we report observations for the four planets within or near the system's habitable zone, the circumstellar region where liquid water could exist on a planetary surface (refs 12-14). These planets do not exhibit prominent spectroscopic signatures at near-infrared wavelengths either, which rules out cloud-free hydrogen-dominated atmospheres for TRAPPIST-1 d, e and f, with significance of 8, 6 and 4 sigma, respectively. Such an atmosphere is instead not excluded for planet g. As high-altitude clouds and hazes are not expected in hydrogen-dominated atmospheres around planets with such insolation (refs 15,16), these observations further support their terrestrial and potentially habitable nature.

We present high resolution ($0.11"\times0.06"$) 3mm ALMA observations of the highly inclined transition disk around the star T Cha. Our continuum image reveals multiple dust structures: an inner disk, a spatially resolved dust gap, and an outer ring. When fitting sky-brightness models to the real component of the 3mm visibilities, we infer that the inner emission is compact ($\le1$au in radius), the gap width is between 18-28 au, and the emission from the outer ring peaks at $\sim36$ au. We compare our ALMA image with previously published 1.6$\mu$m VLT/SPHERE imagery. This comparison reveals that the location of the outer ring is wavelength dependent. More specifically, the peak emission of the 3mm ring is at a larger radial distance than that of the 1.6$\mu$m ring, suggesting that millimeter-sized grains in the outer disk are located further away from the central star than micron-sized grains. We discuss different scenarios to explain our findings, including dead zones, star-driven photoevaporation, and planet-disk interactions. We find that the most likely origin of the dust gap is from an embedded planet, and estimate --- for a single planet scenario --- that T Cha's gap is carved by a $1.2M_{jup}$ planet.

Circumstellar discs with large dust depleted cavities and vigorous accretion onto the central star are often considered signposts for (multiple) giant planet formation. In this letter we show that X-ray photoevaporation operating in discs with modest (factors 3-10) gas-phase depletion of Carbon and Oxygen at large radii (> 15 AU) yield the inner radius and accretion rates for most of the observed discs, without the need to invoke giant planet formation. We present one-dimensional viscous evolution models of discs affected by X-ray photoevaporation assuming moderate gas-phase depletion of Carbon and Oxygen, well within the range reported by recent observations. Our models use a simplified prescription for scaling the X-ray photoevaporation rates and profiles at different metallicity, and our quantitative result depends on this scaling. While more rigorous hydrodynamical modelling of mass loss profiles at low metallicities is required to constrain the observational parameter space that can be explained by our models, the general conclusion that metal sequestering at large radii may be responsible for the observed diversity of transition discs is shown to be robust. Gap opening by giant planet formation may still be responsible for a number of observed transition discs with large cavities and very high accretion rate.

A new piece of evidence supporting the photoevaporation-driven evolution model for low-mass, close-in exoplanets was recently presented by the California-Kepler-Survey. The radius distribution of the Kepler planets is shown to be bimodal, with a ``valley' separating two peaks at 1.3 and 2.6 Rearth. Such an ``evaporation-valley' had been predicted by numerical models previously. Here, we develop a minimal model to demonstrate that this valley results from the following fact: the timescale for envelope erosion is the longest for those planets with hydrogen/helium-rich envelopes that, while only a few percent in weight, double its radius. The timescale falls for envelopes lighter than this because the planet's radius remains largely constant for tenuous envelopes. The timescale also drops for heavier envelopes because the planet swells up faster than the addition of envelope mass. Photoevaporation, therefore, herds planets into either bare cores ~1.3 Rearth, or those with double the core's radius (~2.6 Rearth). This process mostly occurs during the first 100 Myrs when the stars' high energy flux are high and nearly constant. The observed radius distribution further requires that the Kepler planets are clustered around 3 Mearth in mass, are born with H/He envelopes more than a few percent in mass, and that their cores are similar to the Earth in composition. Such envelopes must have been accreted before the dispersal of the gas disks, while the core composition indicates formation inside the ice-line. Lastly, the photoevaporation model fails to account for bare planets beyond ~30-60 days, if these planets are abundant, they may point to a significant second channel for planet formation, resembling the Solar-System terrestrial planets.

(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.

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.

The point at which a protoplanetary disk becomes a debris disk is difficult to identify. To better understand this, here we study the $\sim$40~AU separation binary T~54 in the Chamaeleon I cloud. We derive a K5 spectral type for T~54~A (which dominates the emission of the system) and an age of $\sim$2~Myr. However, the dust disk properties of T~54 are consistent with those of debris disks seen around older and earlier-type stars. At the same time, T~54 has evidence of gas remaining in the disk as indicated by [Ne II], [Ne III], and [O I] line detections. We model the spectral energy distribution of T~54 and estimate that $\sim$3$\times$10$^{-3}$ Earth-masses of small dust grains ($<$0.25~$\mu$m) are present in an optically thin circumbinary disk along with at least $\sim$3$\times$10$^{-7}$ Earth-masses of larger ($>$10~$\mu$m) grains within a circumprimary disk. Assuming a solar-like mixture, we use Ne line luminosities to place a minimum limit on the gas mass of the disk ($\sim$3$\times$10$^{-4}$ Earth-masses) and derive a gas-to-dust mass ratio of $\sim$0.1. We do not detect substantial accretion, but we do see H$\alpha$ in emission in one epoch, suggestive that there may be intermittent dumping of small amounts of matter onto the star. Considering the low dust mass, the presence of gas, and young age of T~54, we conclude that this system is on the bridge between the protoplanetary and debris disk stages.

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.

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.

Accretion onto Classical T Tauri stars is thought to take place through the action of magnetospheric processes, with gas in the inner disk being channeled onto the star's surface by the stellar magnetic field lines. Young stars are known to accrete material in a time-variable manner and the source of this variability remains an open problem, particularly on the shortest (~ day) timescales. Using one-dimensional time-dependent numerical simulations that follow the field line geometry, we find that for plausibly realistic young stars, steady-state transonic accretion occurs naturally in the absence of any other source of variability. However, we show that if the density in the inner disk varies smoothly in time with ~ day long time-scales (e.g., due to turbulence) this complication can lead to the development of shocks in the accretion column. These shocks propagate along the accretion column and ultimately hit the star, leading to rapid, large amplitude changes in the accretion rate. We argue that when these shocks hit the star the observed time-dependence will be a rapid increase in accretion luminosity followed by a slower decline and could be an explanation for some of the short period variability observed in accreting young stars. Our one-dimensional approach bridges previous analytic work to more complicated, multi-dimensional simulations, and observations.

Ultracool dwarfs (UCD; $T_{\rm eff}<\sim3000~$K) cool to settle on the main sequence after $\sim$1 Gyr. For brown dwarfs, this cooling never stops. Their habitable zone (HZ) thus sweeps inward at least during the first Gyr of their lives. Assuming they possess water, planets found in the HZ of UCDs have experienced a runaway greenhouse phase too hot for liquid water prior to entering the HZ. It has been proposed that such planets are desiccated by this hot early phase and enter the HZ as dry worlds. Here we model the water loss during this pre-HZ hot phase taking into account recent upper limits on the XUV emission of UCDs and using 1D radiation-hydrodynamic simulations. We address the whole range of UCDs but also focus on the planets recently found around the $0.08~M_\odot$ dwarf TRAPPIST-1. Despite assumptions maximizing the FUV-photolysis of water and the XUV-driven escape of hydrogen, we find that planets can retain significant amounts of water in the HZ of UCDs, with a sweet spot in the $0.04$-$0.06~M_\odot$ range. We also studied the TRAPPIST-1 system using observed constraints on the XUV-flux. We find that TRAPPIST-1b and c may have lost as much as 15 Earth Oceans and planet d -- which might be inside the HZ -- may have lost less than 1 Earth Ocean. Depending on their initial water contents, they could have enough water to remain habitable. TRAPPIST-1 planets are key targets for atmospheric characterization and could provide strong constraints on the water erosion around UCDs.

Here we present three-dimensional high resolution simulations of Type Ia supernova in the presence of a non-degenerate companion. We find that the presence of a nearby companion leaves a long-lived hole in the supernova ejecta. In particular, we aim to study the long term evolution of this hole as the supernova ejecta interacts with the surrounding interstellar medium. Using estimates for the x-ray emission, we find that the hole generated by the companion remains for many centuries after the interaction between the ejecta and the interstellar medium. We also show that the hole is discernible over a wide range of viewing angles and companion masses.

We discuss a generalization of the classic Keplerian disk test problem allowing for both pressure and rotational support, as a method of testing astrophysical codes incorporating both gravitation and hydrodynamics. We argue for the inclusion of pressure in rotating disk simulations on the grounds that realistic, astrophysical disks exhibit non-negligible pressure support. We then apply this test problem to examine the performance of various smoothed particle hydrodynamics (SPH) methods incorporating a number of improvements proposed over the years to help SPH better address problems noted in modeling the classical gravitation only Keplerian disk. We also apply this test to a newly developed extension of SPH based on reproducing kernels called CRKSPH. Counterintuitively, we find that pressure support worsens the performance of traditional SPH on this problem, causing unphysical collapse away from the steady-state disk solution even more rapidly than the purely gravitational problem, whereas CRKSPH greatly reduces this error.

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 6300A. 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.

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.

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.

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.

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.

Creating spherical initial conditions in smoothed particle hydrodynamics simulations that are spherically conformal is a difficult task. Here, we describe two algorithmic methods for evenly distributing points on surfaces, that when paired can be used to build 3D spherical objects with optimal equipartition of volume between particles, commensurate with an arbitrary, radial density function. We demonstrate the efficacy of our method against stretched lattice arrangements on the metrics of hydrodynamic stability, spherical conformity, and the harmonic power distribution of gravitational settling oscillations. We further demonstrate how our method is highly optimized for simulating multi-material spheres, such as planets with core-mantle boundaries.

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.

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.

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.

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.

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).

Photo-evaporation and planet formation have both been proposed as mechanisms responsible for the creation of a transition disc. We have studied their combined effect through a suite of 2d simulations of protoplanetary discs undergoing X-ray photoevaporation with an embedded giant planet. In a previous work we explored how the formation of a giant planet triggers the dispersal of the inner disc by photo-evaporation at earlier times than what would have happened otherwise. This is particularly relevant for the observed transition discs with large holes and high mass accretion rates that cannot be explained by photo-evaporation alone. In this work we significantly expand the parameter space investigated by previous simulations. In addition, the updated model includes thermal sweeping, needed for studying the complete dispersal of the disc. After the removal of the inner disc the disc is a non accreting transition disc, an object that is rarely seen in observations. We assess the relative length of this phase, to understand if it is long lived enough to be found observationally. Depending on the parameters, especially on the X-ray luminosity of the star, we find that the fraction of time spent as a non-accretor greatly varies. We build a population synthesis model to compare with observations and find that in general thermal sweeping is not effective enough to destroy the outer disc, leaving many transition discs in a relatively long lived phase with a gas free hole, at odds with observations. We discuss the implications for transition disc evolution. In particular, we highlight the current lack of explanation for the missing non-accreting transition discs with large holes, which is a serious issue in the planet hypothesis.

By means of radiative transfer simulation we study the evolution of the far-infrared colours of protoplanetary discs undergoing inside-out dispersal, often referred to as transition discs. We show that a brightening of the mid and far-infrared emission from these objects is a natural consequence of the removal of the inner disc. Our results can fully explain recent observations of transition discs in the Chamaleon and Lupus star forming regions from the Herschel Gould Belt Survey, which show a higher median for the 70?um (Herschel PACS 1) band of known transition objects compared with primordial discs. Our theoretical results hence support the suggestion that the 70?um band may be a powerful diagnostic for the identification of transition discs from photometry data, provided that the inner hole is larger than tens of AU, depending on spectral type. Furthermore we show that a comparison of photometry in the K , 12?um and 7u0?m bands to model tracks can provide a rough, but quick estimate of the inner hole size of these objects, provided their inclination is below ?85 degrees and the inner hole size is again larger than tens of AU.

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 timescale for the wind to reach its sonic radius. We also estimate the momentum flux from time-dependent one-dimensional hydrodynamical simulations. We find that only in a very narrow regime in planet radius, mass and stellar radiation flux is a wind capable of exerting a significant torque on the planet's orbit. Similar to the Yarkovsky effect, the wind causes the planet to drift outward 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 percent of its semi-major axis. While this change is small, it places constraints on the evolution of resonant pairs such as Kepler 36 b and c.

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.

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.

Hot Jupiters can experience mass loss driven by heating from UV radiation from their host stars, and this flow is often controlled by magnetic fields. More specifically, near the planetry surface, the magnetic pressure dominates the ram pressure of the outflow by several orders of magnitude. After leaving the vicinity of the planet, the flow must connect onto the background environment provided by the stellar wind and the stellar magnetic field. This contribution considers magnetically controlled planetary outflows and extends previous work by comparing two different geometries for the background magnetic field provided by the star. In the first case, stellar field is assumed to retain the form of a dipole, which is anti-aligned with the dipole field of the planet. In the second case, the stellar outflow opens up the stellar magnetic field structure so that the background field at the location of the planet is perpendicular to the planetary dipole. Using numerical simulations, we consider the launch of the planetary wind with these field configurations.

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.

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.

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.

We investigate under what circumstances an embedded planet in a protoplanetary disc may sculpt the dust distribution such that it observationally presents as a `transition' disc. We concern ourselves with `transition' discs that have large holes ($\gtrsim 10$ AU) and high accretion rates ($\sim 10^{-9}-10^{-8}$ M$_\odot$ yr$^{-1}$). Particularly, those discs which photoevaporative models struggle to explain. Assuming the standard picture for how massive planets sculpt their parent discs, along with the observed accretion rates in `transition' discs, we find that the accretion luminosity from the forming planet is significant, and can dominate over the stellar luminosity at the gap edge. This planetary accretion luminosity can apply a significant radiation pressure to small ($s\lesssim 1\mu$m) dust particles provided they are suitably decoupled from the gas. Secular evolution calculations that account for the evolution of the gas and dust components in a disc with an embedded, accreting planet, show that only with the addition of the radiation pressure can we explain the full observed characteristics of a `transition' disc (NIR dip in the SED, mm cavity and high accretion rate). At suitably high planet masses ($\gtrsim 3-4$ M$_J$), radiation pressure from the accreting planet is able to hold back the small dust particles, producing a heavily dust-depleted inner disc that is optically thin (vertically and radially) to Infra-Red radiation. We use our models to calculate synthetic observations and present a observational evolutionary scenario for a forming planet, sculpting its parent disc. The planet-disc system will present as a `transition' disc with a dip in the SED, only when the planet mass and planetary accretion rate is high enough. At other times it will present as a disc with a primordial SED, but with a cavity in the mm, as observed in a handful of protoplanetary discs.

We analyse current measurements of accretion rates onto pre-main sequence stars as a function of stellar mass, and conclude that the steep dependance of accretion rates on stellar mass is real and not driven by selection/detection threshold, as has been previously feared. These conclusions are reached by means of statistical tests including a survival analysis which can account for upper limits. The power-law slope of the Mdot-M* relation is found to be in the range of 1.6-1.9 for young stars with masses lower than 1 Msun. The measured slopes and distributions can be easily reproduced by means of a simple disc model which includes viscous accretion and X-ray photoevaporation. We conclude that the Mdot-M* relation in pre-main sequence stars bears the signature of disc dispersal by X-ray photoevaporation, suggesting that the relation is a straight- forward consequence of disc physics rather than an imprint of initial conditions.

(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.

(Abridged) Inspired by the Kepler 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. We construct a theoretical population of planets with varying core masses, envelope masses, orbital separations, and stellar spectral types, and compare these against 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 bimodal distribution in planet sizes with a deficit of planets around 2R_E. 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 distribution is naturally explained by the evaporation model, where, depending on their X-ray exposure, close-in planets can either hold on to hydrogen envelopes 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 overlaid with hydrogen envelopes, but few of them should have initial masses above 20 M_E, and the majority of them should have core masses of a few Earth masses.

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.

We assume a scenario in which transition discs (i.e. discs around young stars that have signatures of cool dust but lack significant near infra-red emission from warm dust) are associated with the presence of planets (or brown dwarfs). These are assumed to filter the dust content of any gas flow within the planetary orbit and produce an inner `opacity hole'. In order to match the properties of transition discs with the largest (~50 A.U. scale) holes, we place such `planets' at large radii in massive discs and then follow the evolution of the tidally coupled disc-planet system, comparing the system's evolution in the plane of mm flux against hole radius with the properties of observed transition discs. We find that, on account of the high disc masses in these systems, all but the most massive `planets' (100 Jupiter masses) are conveyed to small radii by Type II migration without significant fading at millimetre wavelengths. Such behaviour would contradict the observed lack of mm bright transition discs with small (<10 A.U.) holes. On the other hand, imaging surveys clearly rule out the presence of such massive companions in transition discs. We conclude that this is a serious problem for models that seek to explain transition discs in terms of planetary companions unless some mechanism can be found to halt inward migration and/or suppress mm flux production. We suggest that the dynamical effects of substantial accretion on to the planet/through the gap may offer the best prospect for halting such migration but that further long term simulations are required to clarify this issue.

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.

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.

We examine the distribution of transition discs as a function of 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 we 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 sub-sample have systematically higher accretion rates and inner hole radii than those in the faint mm sub-sample, along with being systematically weighted to earlier spectral types.

We consider the evaporation of close in planets by the star's intrinsic EUV and X-ray radiation. We calculate evaporation rates by solving the hydrodynamical problem for planetary evaporation including heating from both X-ray and EUV radiation. We show that most close-in planets (a<0.1 AU) are evaporating hydrodynamically, with the evaporation occurring in two distinct regimes: X-ray driven, in which the X-ray heated flow contains a sonic point, and EUV driven, in which the X-ray region is entirely sub-sonic. The mass-loss rates scale as L_X/a^2 for X-ray driven evaporation, and as Phi_*^{1/2}/a for EUV driven evaporation at early times, with mass-loss rates of order 10e10-10e14 g/s. No exact scaling exists for the mass-loss rate with planet mass and planet radius, however, in general evaporation proceeds more rapidly for planets with lower densities and higher masses. Furthermore, we find that in general the transition from X-ray driven to EUV driven evaporation occurs at lower X-ray luminosities for planets closer to their parent stars and for planets with lower densities. Coupling our evaporation models to the evolution of the high energy radiation - which falls with time - we are able to follow the evolution of evaporating planets. We find that most planets start off evaporating in the X-ray driven regime, but switch to EUV driven once the X-ray luminosity falls below a critical value. The evolution models suggest that while `hot Jupiters' are evaporating, they are not evaporating at a rate sufficient to remove the entire gaseous envelope on Gyr time-scales. However, we do find that close in Neptune mass planets are more susceptible to complete evaporation of their envelopes. Thus we conclude that planetary evaporation is more important for lower mass planets, particularly those in the `hot Neptune'/`super Earth' regime.

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.

We address the question of whether protoplanetary discs around low mass stars (e.g. M-dwarfs) may be longer lived than their solar-type counterparts. This question is particularly relevant to assess the planet-making potential of these lower mass discs. Given the uncertainties inherent to age-dating young stars, we propose an alternative approach that is to analyse the spatial distribution of disc-bearing low-mass stars and compare it to that of disc-bearing solar-type stars in the same cluster. A significant age difference between the two populations should be reflected in their average nearest neighbour distance (normalised to the number of sources), where the older population should appear more spread out. To this aim, we perform a Minimum Spanning Tree (MST) analysis on the spatial distribution of disc-bearing young stellar objects (YSOs) in six nearby low mass star forming regions. We find no evidence for significant age differences between the disc-bearing low-mass (later than M2) and 'solar-type' (earlier than M2) stars in these regions. We model our results by constructing and analysing synthetic fractal distributions that we evolve for typical values of the velocity dispersions. A comparison of simple models to our MST analysis suggests that the lifetime of discs around M-stars is similar to that of discs around solar-type stars. Furthermore, a model-independent spatial analysis of the observations robustly shows that any age differences between the two samples must be smaller than the average age difference between disc-bearing classical T-Tauri stars and disc-less Weak-Lined T-Tauri stars.

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.

The Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated, with special consideration for high velocity gas. We find the cloud to have a mass of 203 solar masses, and to be gravitationally bound, and that the kinetic energy of the outflowing gas is approximately seventy percent of the turbulent kinetic energy of the cloud. We identify compact outflows towards some of the submillimetre Class 0/I sources in the region

We present a calculation of the atomic and low-ionisation emission line spectra of photoevaporating protoplanetary discs. Line luminosities and profiles are obtained from detailed photoionisation 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 OI 6300A 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 OI 6300A line and other forbidden lines observed by Hartigan et al., which covered a large sample of T-Tauri stars. We find that the OI 6300A and the NeII 12.8um lines are predominantly produced in the X-ray-driven wind and thus appear blue-shifted by a few km/s for some of the systems when observed at non-edge-on inclinations. We note however that blue-shifts 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.

(Abridged) We present the first radiation-hydrodynamic model of a protoplanetary disc irradiated with an X-EUV spectrum. In a model where the total ionizing luminosity is divided equally between X-ray and EUV luminosity, we find a photoevaporation rate of 1.4e-8 M_sun/yr, which is two orders of magnitude greater than the case of EUV photoevaporation alone. Thus it is clear that the X-rays are the dominant driving mechanism for photoevaporation. This can be understood inasmuch as X-rays are capable of penetrating much larger columns (~1e22 cm^-2) and can thus effect heating in denser regions and at larger radius than the EUV can. The radial extent of the launching region of the X-ray heated wind is 1-70AU compared with the pure EUV case where the launch region is concentrated around a few AU. When we couple our wind mass-loss rates with models for the disc's viscous evolution, we find that, as in the pure EUV case, there is a photoevaporative switch, such that an inner hole develops at ~ 1 AU at the point that the accretion rate in the disc drops below the wind mass loss rate. At this point, the remaining disc material is quickly removed in the final 15-20% of the disc's lifetime. This is consistent with the 1e5 yr transitional timescale estimated from observations of T-Tauri stars. We also caution that although our mass-loss rates are high compared to some accretion rates observed in young stars, our model has a rather large X-ray luminosity of 2e30 erg/s; further modeling is required in order to investigate the evolutionary implications of the large observed spread of X-ray luminosities in T-Tauri stars.

We describe a modified form of Smoothed Particle Hydrodynamics (SPH) in which the specific thermal energy equation is based on a compatibly differenced formalism, guaranteeing exact conservation of the total energy. We compare the errors and convergence rates of the standard and compatible SPH formalisms on analytic test problems involving shocks. We find that the new compatible formalism reliably achieves the expected first-order convergence in such tests, and in all cases improves the accuracy of the numerical solution over the standard formalism.

We have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including SPH and grid methods with fixed, deformable or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gas dynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of two. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to about 10%, and the ratio of specific kinetic to thermal energies of the gas agree to about 5%. Various factors contribute to the spread in calculated cluster properties, including differences in the internal timing of the simulations. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modelled.

The growth of structure from scale-free initial conditions is one of the most important tests of cosmological simulation methods, providing a realistically complex problem in which numerical results can be compared to rigorous analytic scaling laws. Previous studies of this problem have incorporated gravitational dynamics and adiabatic gas dynamics, but radiative cooling, an essential element of the physics of galaxy formation, normally violates the conditions necessary for self-similar evolution. We show that for any specified value of the initial power spectrum index n [where P(k)~k^n], there is a family of power-law cooling functions that preserves self-similarity by ensuring that the cooling time in an object of the characteristic mass M_* is a fixed fraction t_C of the Hubble time. We perform hydrodynamic numerical simulations with an Einstein-de Sitter cosmology, a baryon fraction of 5%, Gaussian initial conditions, two different power spectrum indices, and four values of t_C for each index, ranging from no cooling to strong cooling. We restrict the numerical simulations to two dimensions in order to allow exploration of a wide parameter space with adequate dynamic range. In all cases, the simulations are remarkably successful at reproducing the analytically predicted scalings of the mass function of dissipated objects and the gas temperature distributions and cooled gas fractions in collapsed systems. While similar success with 3-D simulations must still be demonstrated, our results have encouraging implications for numerical studies of galaxy formation, indicating that simulations with resolution comparable to that in many current studies can accurately follow the collapse and dissipation of baryons into the dense, cold systems where star formation is likely to occur.

We analyze hierarchical structure formation based on scale-free initial conditions in an Einstein-de Sitter universe, including a baryonic component. We present three independent, smoothed particle hydrodynamics (SPH) simulations, performed with two different SPH codes (TreeSPH and P3MSPH) at two resolutions. Each simulation is based upon identical initial conditions, which consist of Gaussian distributed initial density fluctuations that have an n=-1 power spectrum. The baryonic material is modeled as an ideal gas subject only to shock heating and adiabatic heating and cooling. The evolution is expected to be self-similar in time, and under certain restrictions we identify the expected scalings for many properties of the distribution of collapsed objects in all three realizations. The distributions of dark matter masses, baryon masses, and mass and emission weighted temperatures scale quite reliably. However, the density estimates in the central regions of these structures are determined by the degree of numerical resolution. As a result, mean gas densities and luminosities obey the expected scalings only when calculated within a limited dynamic range in density contrast. The temperatures and luminosities of the groups show tight correlations with the baryon masses, which can be well-represented by power-laws. The Press-Schechter (PS) approximation predicts the distribution of group dark matter and baryon masses fairly well, though it tends to overestimate the baryon masses. Combining the PS mass distribution with the measured relations for T(M) and L(M) predicts the temperature and luminosity distributions reasonably, though there are some discrepancies at high temperatures/luminosities. The three simulations agree well for the properties of groups that are resolved by 32 or more particles.

This paper presents an alternative formulation of the ASPH algorithm for evolving anisotropic smoothing kernels, in which the geometric approach of Shapiro et al. (1996; Paper I) is replaced by an approach involving a local transformation of coordinates to those in which the underlying anisotropic volume changes appear to be isotropic. The ASPH method is presented in 2D and 3D, including a number of details not previously included in Paper I, some of which represent either advances or different choices with respect to Paper I. Among the advances included here are an asynchronous time-integration scheme with different time steps for different particles and the generalization of the ASPH method to 3D. The shock-tracking algorithm described in Paper I for locally adapting the artificial viscosity to restrict viscous heating just to particles encountering shocks, is not included here. Instead, we adopt a different interpolation kernel for use with the artificial viscosity, which has the effect of spatially localizing effects of the artificial viscosity. This version of the ASPH method in 2D and 3D is then applied to a series of 1D, 2D, and 3D test problems, and the results are compared to those of standard SPH applied to the same problems. These include the problem of cosmological pancake collapse, the Riemann shock tube, cylindrical and spherical Sedov blast waves, the collision of two strong shocks, and problems involving shearing disks intended to test the angular momentum conservation properties of the method. These results further support the idea that ASPH has significantly better resolving power than standard SPH for a wide range of problems, including that of cosmological structure formation. (Abridged)

We investigate the properties of hybrid gravitational/hydrodynamical simulations, examining both the numerics and the general physical properties of gravitationally driven, hierarchical collapse in a mixed baryonic/dark matter fluid. We demonstrate that, under certain restrictions, such simulations converge with increasing resolution to a consistent solution. The dark matter achieves convergence provided that the relevant scales dominating nonlinear collapse are resolved. If the gas has a minimum temperature (as expected when intergalactic gas is heated by photoionization due to the ultraviolet background) and the corresponding Jeans mass is resolved, then the baryons also converge. However, if there is no minimum baryonic collapse mass or if this scale is not resolved, then the baryon results err in a systematic fashion. In such a case, as resolution is increased the baryon distribution tends toward a higher density, more tightly bound state. We attribute this to the fact that under hierarchical structure formation on all scales there is always an earlier generation of smaller scale collapses, causing shocks which irreversibly alter the state of the baryon gas. In a simulation with finite resolution we miss such earlier generation collapses, unless a physical scale is introduced below which structure formation is suppressed in the baryons. We also find that the baryon/dark matter ratio follows a characteristic pattern, such that collapsed structures possess a baryon enriched core (enriched by factors of 2 or more over the universal average) which is embedded within a dark matter halo, even without accounting for radiative cooling of the gas. The dark matter is unaffected by changing the baryon distribution (at least in the dark matter dominated case investigated here).