Identification of habitable planets beyond our solar system is a key goal of current and future space missions. Yet habitability depends not only on the stellar irradiance, but equally on constituent parts of the planetary atmosphere. Here we show, for the first time, that radiatively active mineral dust will have a significant impact on the habitability of Earth-like exoplanets. On tidally-locked planets, dust cools the day-side and warms the night-side, significantly widening the habitable zone. Independent of orbital configuration, we suggest that airborne dust can postpone planetary water loss at the inner edge of the habitable zone, through a feedback involving decreasing ocean coverage and increased dust loading. The inclusion of dust significantly obscures key biomarker gases (e.g. ozone, methane) in simulated transmission spectra, implying an important influence on the interpretation of observations. We demonstrate that future observational and theoretical studies of terrestrial exoplanets must consider the effect of dust.

The majority of potentially habitable exoplanets detected orbit stars cooler than the Sun, and therefore are irradiated by a stellar spectrum peaking at longer wavelengths than that incident on Earth. Here, we present results from a set of simulations of tidally-locked terrestrial planets orbiting three different host stars to isolate the effect of the stellar spectra on the simulated climate. Specifically, we perform simulations based on TRAPPIST-1e, adopting an Earth-like atmosphere and using the UK Met Office Unified Model in an idealised 'aqua-planet' configuration. Whilst holding the planetary parameters constant, including the total stellar flux (900 W/m$^2$) and orbital period (6.10 Earth days), we compare results between simulations where the stellar spectrum is that of a quiescent TRAPPIST-1, Proxima Centauri and the Sun. The simulations with cooler host stars had an increased proportion of incident stellar radiation absorbed directly by the troposphere compared to the surface. This, in turn, led to an increase in the stability against convection, a reduction in overall cloud coverage on the dayside (reducing scattering), leading to warmer surface temperatures. The increased direct heating of the troposphere also led to more efficient heat transport from the dayside to the nightside and, therefore, a reduced day-night temperature contrast. We inferred that planets with an Earth-like atmosphere orbiting cooler stars had lower dayside cloud coverage, potentially allowing habitable conditions at increased orbital radii, compared to similar planets orbiting hotter stars for a given planetary rotation rate.

We present four new secondary eclipse observations for the ultrahot Jupiter WASP-121b acquired using the Hubble Space Telescope Wide Field Camera 3. The eclipse depth is measured to a median precision of 60ppm across 28 spectroscopic channels spanning the 1.12-1.64 micron wavelength range. This is a considerable improvement to the 90ppm precision we achieved previously for a single eclipse observation using the same observing setup. Combining these data with those reported at other wavelengths, a blackbody spectrum for WASP-121b is ruled out at >6-sigma confidence and we confirm the interpretation of previous retrieval analyses that found the data is best explained by a dayside thermal inversion. The updated spectrum clearly resolves the water emission band at 1.3-1.6 micron, with higher signal-to-noise than before. It also fails to reproduce a bump in the spectrum at 1.25 micron derived from the first eclipse observation, which had tentatively been attributed to VO emission. We conclude the latter was either a statistical fluctuation or a systematic artefact specific to the first eclipse dataset.

We present results from a set of simulations using a fully coupled three-dimensional (3D) chemistry-radiation-hydrodynamics model and investigate the effect of transport of chemical species by the large-scale atmospheric flow in hot Jupiter atmospheres. We couple a flexible chemical kinetics scheme to the Met Office Unified Model which enables the study of the interaction of chemistry, radiative transfer and fluid dynamics. We use a newly-released "reduced" chemical network comprising 30 chemical species that has been specifically developed for application in 3D atmosphere models. We simulate the atmospheres of the well-studied hot Jupiters HD~209458b and HD~189733b which both have dayside--nightside temperature contrasts of several hundred Kelvin and superrotating equatorial jets. We find qualitatively quite different chemical structures between the two planets, particularly for methane (CH$_4$), when advection of chemical species is included. Our results show that consideration of 3D chemical transport is vital in understanding the chemical composition of hot Jupiter atmospheres. 3D mixing leads to significant changes in the abundances of absorbing gas-phase species compared with what would be expected by assuming local chemical equilibrium, or from models including 1D - and even 2D - chemical mixing. We find that CH$_4$, carbon dioxide (CO$_2$) and ammonia (NH$_3$) are particularly interesting as 3D mixing of these species leads to prominent signatures of out-of-equilibrium chemistry in the transmission and emission spectra, detectable with near-future instruments.

Using a 3D general circulation model (GCM), we investigate the sensitivity of the climate of tidally-locked Earth-like exoplanets, Trappist-1e and Proxima Centauri b, to the choice of a convection parameterization. Compared to a mass-flux convection parameterization, a simplified convection adjustment parameterization leads to a $>$60% decrease of the cloud albedo, increasing the mean day-side temperature by $\approx$10 K. The representation of convection also affects the atmospheric conditions of the night side, via a change in planetary-scale wave patterns. As a result, using the convection adjustment scheme makes the night-side cold traps warmer by 17-36 K for the planets in our simulations. The day-night thermal contrast is sensitive to the representation of convection in 3D GCM simulations, so caution should be taken when interpreting emission phase curves. The choice of convection treatment, however, does not alter the simulated climate enough to result in a departure from habitable conditions, at least for the atmospheric composition and planetary parameters used in our study. The near-surface conditions both in the Trappist-1e and Proxima b cases remain temperate, allowing for an active water cycle. We further advance our analysis using high-resolution model experiments, in which atmospheric convection is simulated explicitly. Our results suggest that in a hypothetical global convection-permitting simulation the surface temperature contrast would be higher than in the coarse-resolution simulations with parameterized convection. In other words, models with parameterized convection may overestimate the inter-hemispheric heat redistribution efficiency.

Upcoming telescopes such as the James Webb Space Telescope (JWST), or the Extremely Large Telescope (ELTs), may soon be able to characterize, through transmission, emission or reflection spectroscopy, the atmospheres of rocky exoplanets orbiting nearby M dwarfs. One of the most promising candidates is the late M dwarf system TRAPPIST-1 which has seven known transiting planets for which Transit Timing Variation (TTV) measurements suggest that they are terrestrial in nature, with a possible enrichment in volatiles. Among these seven planets, TRAPPIST-1e seems to be the most promising candidate to have habitable surface conditions, receiving ~66 % of the Earth's incident radiation, and thus needing only modest greenhouse gas inventories to raise surface temperatures to allow surface liquid water to exist. TRAPPIST-1e is therefore one of the prime targets for JWST atmospheric characterization. In this context, the modeling of its potential atmosphere is an essential step prior to observation. Global Climate Models (GCMs) offer the most detailed way to simulate planetary atmospheres. However, intrinsic differences exist between GCMs which can lead to different climate prediction and thus observability of gas and/or cloud features in transmission and thermal emission spectra. Such differences should preferably be known prior to observations. In this paper we present a protocol to inter-compare planetary GCMs. Four testing cases are considered for TRAPPIST-1e but the methodology is applicable to other rocky exoplanets in the Habitable Zone. The four test cases included two land planets composed with a modern Earth and pure CO2 atmospheres, respectively, and two aqua planets with the same atmospheric compositions. Currently, there are four participating models (LMDG, ROCKE-3D, ExoCAM, UM), however this protocol is intended to let other teams participate as well.

We use the Met Office Unified Model to explore the potential of a tidally locked M dwarf planet, nominally Proxima Centauri b irradiated by a quiescent version of its host star, to sustain an atmospheric ozone layer. We assume a slab ocean surface layer, and an Earth-like atmosphere of nitrogen and oxygen with trace amounts of ozone and water vapour. We describe ozone chemistry using the Chapman mechanism and the hydrogen oxide (HO$_x$, describing the sum of OH and HO$_2$) catalytic cycle. We find that Proxima Centauri radiates with sufficient UV energy to initialize the Chapman mechanism. The result is a thin but stable ozone layer that peaks at 0.75 parts per million at 25 km. The quasi-stationary distribution of atmospheric ozone is determined by photolysis driven by incoming stellar radiation and by atmospheric transport. Ozone mole fractions are smallest in the lowest 15 km of the atmosphere at the sub-stellar point and largest in the nightside gyres. Above 15 km the ozone distribution is dominated by an equatorial jet stream that circumnavigates the planet. The nightside ozone distribution is dominated by two cyclonic Rossby gyres that result in localized ozone hotspots. On the dayside the atmospheric lifetime is determined by the HO$_x$ catalytic cycle and deposition to the surface, with nightside lifetimes due to chemistry much longer than timescales associated with atmospheric transport. Surface UV values peak at the substellar point with values of 0.01 W/m$^2$, shielded by the overlying atmospheric ozone layer but more importantly by water vapour clouds.

Context: The anomalously large radii of hot Jupiters has long been a mystery. However, by combining both theoretical arguments and 2D models, a recent study has suggested that the vertical advection of potential temperature leads to an adiabatic temperature profile in the deep atmosphere hotter than the profile obtained with standard 1D models. Aims: In order to confirm the viability of that scenario, we extend this investigation to three dimensional, time-dependent, models. Methods: We use a 3D GCM, DYNAMICO to perform a series of calculations designed to explore the formation and structure of the driving atmospheric circulations, and detail how it responds to changes in both upper and deep atmospheric forcing. Results: In agreement with the previous, 2D, study, we find that a hot adiabat is the natural outcome of the long-term evolution of the deep atmosphere. Integration times of order $1500$ years are needed for that adiabat to emerge from an isothermal atmosphere, explaining why it has not been found in previous hot Jupiter studies. Models initialised from a hotter deep atmosphere tend to evolve faster toward the same final state. We also find that the deep adiabat is stable against low-levels of deep heating and cooling, as long as the Newtonian cooling time-scale is longer than $\sim 3000$ years at $200$ bar. Conclusions: We conclude that the steady-state vertical advection of potential temperature by deep atmospheric circulations constitutes a robust mechanism to explain hot Jupiter inflated radii. We suggest that future studies of hot Jupiters are evolved for a longer time than currently done, and, when possible, include models initialised with a hot deep adiabat. We stress that this mechanism stems from the advection of entropy by irradiation induced mass flows and does not require (finely tuned) dissipative process, in contrast with most previously suggested scenarios.

Context. Atmospheric superrotating flows at the equator are an almost ubiquitous result of simulations of hot Jupiters, and a theory explaining how this zonally coherent flow reaches an equilibrium has been developed in the literature. However, this understanding relies on the existence of either an initial superrotating or a sheared flow, coupled with a slow evolution such that a linear steady state can be reached. Aims. A consistent physical understanding of superrotation is needed for arbitrary drag and radiative timescales, and the relevance of considering linear steady states needs to be assessed. Methods. We obtain an analytical expression for the structure, frequency and decay rate of propagating waves in hot Jupiter atmospheres around a state at rest in the 2D shallow water beta plane limit. We solve this expression numerically and confirm the robustness of our results with a 3D linear wave algorithm. We then compare with 3D simulations of hot Jupiter atmospheres and study the non linear momentum fluxes. Results. We show that under strong day night heating the dynamics does not transit through a linear steady state when starting from an initial atmosphere in solid body rotation. We further show that non linear effects favour the initial spin up of superrotation and that the acceleration due to the vertical component of the eddy momentum flux is critical to the initial development of superrotation. Conclusions. Overall, we describe the initial phases of the acceleration of superrotation, including consideration of differing radiative and drag timescales, and conclude that eddy-momentum driven superrotating equatorial jets are robust, physical phenomena in simulations of hot Jupiter atmospheres.

Context. The study of linear waves and instabilities is necessary to understand the physical evolution of an atmosphere, and can provide physical interpretation of the complex flows found in simulations performed using Global Circulation Models (GCM). In particular, the acceleration of superrotating flow at the equator of hot Jupiters has mostly been studied under several simplifying assumptions, the relaxing of which may impact final results. Aims. We develop and benchmark a publicly available algorithm to identify the eigenmodes of an atmosphere around any initial steady state. We also solve for linear steady states. Methods. We linearise the hydrodynamical equations of a planetary atmosphere in a steady state with arbitrary velocities and thermal profile. We then discretise the linearised equations on an appropriate staggered grid, and solve for eigenvectors and linear steady solutions with the use of a parallel library for linear algebra: ScaLAPACK. We also implement a posteriori calculation of an energy equation in order to obtain more information on the underlying physics of the mode. Results. Our code is benchmarked against classical wave and instability test cases in multiple geometries. The steady linear circulation calculations also reproduce expected results for the atmosphere of hot Jupiters. We finally show the robustness of our energy equation, and its power to obtain physical insight into the modes. Conclusions. We have developed and benchmarked a code for the study of linear processes in planetary atmospheres, with an arbitrary steady state. The calculation of an a posteriori energy equation provides both increased robustness and physical meaning to the obtained eigenmodes. This code can be applied to various problems, and notably to further study the initial spin up of superrotation of GCM simulations of hot Jupiter.

In this white paper, we recommend the European Space Agency plays a proactive role in developing a global collaborative effort to construct a large high-contrast imaging space telescope, e.g. as currently under study by NASA. Such a mission will be needed to characterize a sizable sample of temperate Earth-like planets in the habitable zones of nearby Sun-like stars and to search for extraterrestrial biological activity. We provide an overview of relevant European expertise, and advocate ESA to start a technology development program towards detecting life outside the Solar system.

We present results from 3D radiative-hydrodynamical simulations of HD 209458b with a fully coupled treatment of clouds using the EddySed code, critically, including cloud radiative feedback via absorption and scattering. We demonstrate that the thermal and optical structure of the simulated atmosphere is markedly different, for the majority of our simulations, when including cloud radiative effects, suggesting this important mechanism can not be neglected. Additionally, we further demonstrate that the cloud structure is sensitive to not only the cloud sedimentation efficiency (termed $f_{\textrm{sed}}$ in EddySed), but also the temperature-pressure profile of the deeper atmosphere. We briefly discuss the large difference between the resolved cloud structures of this work, adopting a phase-equilibrium and parameterised cloud model, and our previous work incorporating a cloud microphysical model, although a fairer comparison where, for example, the same list of constituent condensates is included in both treatments, is reserved for a future work. Our results underline the importance of further study into the potential condensate size distributions and vertical structures, as both strongly influence the radiative impact of clouds on the atmosphere. Finally, we present synthetic observations from our simulations reporting an improved match, over our previous cloud-free simulations, to the observed transmission, HST WFC3 emission and 4.5 $\mu$m Spitzer phase curve of HD 209458b. Additionally, we find all our cloudy simulations have an apparent albedo consistent with observations.

We present results from one-dimensional atmospheric simulations investigating the effect of varying the carbon-to-oxygen (C/O) ratio on the thermal structure, chemical composition and transmission and emission spectra, for irradiated hydrogen-dominated atmospheres. We find that each of these properties of the atmosphere are strongly dependent on the individual abundances (relative to hydrogen) of carbon and oxygen. We confirm previous findings that different chemical equilibrium compositions result from different sets of element abundances but with the same C/O ratio. We investigate the effect of this difference in composition on the thermal structure and simulated spectra. We also simulate observations using the PandExo tool and show that these differences are observationally significant with current (i.e. Hubble Space Telescope) and future (i.e. James Webb Space Telescope) instruments. We conclude that it is important to consider the full set of individual element abundances, with respect to hydrogen, rather than the ratios of only two elements, such as the C/O ratio, particularly when comparing model predictions with observed transmission and emission spectra.

In the original manuscript arXiv:1710.10269, we presented a grid of forward model transmission spectra for hot Jupiter exoplanets. However, we recently identified an error in the treatment of rainout in our 1D atmosphere model ATMO. We explain the error, the correction, its validation and the changed conclusions in this erratum. The correction of this error led to changes in the equilibrium chemical abundances using rainout condensation and thereby the transmission spectra. We note that this error only affects the online library that includes rainout condensation, the library with local condensation (without rainout) is unaffected. The online library was updated with the correction in July 2018.

We present significant differences in the simulated atmospheric flow for warm, tidally-locked small Neptunes and super Earths (based on a nominal GJ 1214b) when solving the simplified, and commonly used, primitive dynamical equations or the full Navier-Stokes equations. The dominant prograde, superrotating zonal jet is markedly different between the simulations which are performed using practically identical numerical setups, within the same model. The differences arise due to the breakdown of the so-called `shallow-fluid' and traditional approximations, which worsens when rotation rates are slowed, and day-night temperature contrasts are increased. The changes in the zonal advection between simulations solving the full and simplified equations, give rise to significant differences in the atmospheric redistribution of heat, altering the position of the hottest part of the atmosphere and temperature contrast between the day and night sides. The implications for the atmospheric chemistry and, therefore, observations need to be studied with a model including a more detailed treatment of the radiative transfer and chemistry. Small Neptunes and super Earths are extremely abundant and important, potentially bridging the structural properties (mass, radius, composition) of terrestrial and gas giant planets. Our results indicate care is required when interpreting the output of models solving the primitive equations of motion for such planets.

Simulated exoplanet transmission spectra are critical for planning and interpretation of observations and to explore the sensitivity of spectral features to atmospheric thermochemical processes. We present a publicly available generic model grid of planetary transmission spectra, scalable to a wide range of H$_2$/He dominated atmospheres. The grid is computed using the 1D/2D atmosphere model ATMO for two different chemical scenarios, first considering local condensation only, secondly considering global condensation and removal of species from the atmospheric column (rainout). The entire grid consists of 56,320 model simulations across 22 equilibrium temperatures (400 - 2600 K), four planetary gravities (5 - 50 ms$^{-2}$), five atmospheric metallicities (1x - 200x), four C/O ratios (0.35 - 1.0), four scattering haze parameters, four uniform cloud parameters, and two chemical scenarios. We derive scaling equations which can be used with this grid, for a wide range of planet-star combinations. We validate this grid by comparing it with other model transmission spectra available in the literature. We highlight some of the important findings, such as the rise of SO$_2$ features at 100x solar metallicity, differences in spectral features at high C/O ratios between two condensation approaches, the importance of VO features without TiO to constrain the limb temperature and features of TiO/VO both, to constrain the condensation processes. Finally, this generic grid can be used to plan future observations using the HST, VLT, JWST and various other telescopes. The fine variation of parameters in the grid also allows it to be incorporated in a retrieval framework, with various machine learning techniques.

In this paper we present three-dimensional atmospheric simulations of the hot Jupiter HD~189733b under two different scenarios: local chemical equilibrium and including advection of the chemistry by the resolved wind. Our model consistently couples the treatment of dynamics, radiative transfer and chemistry, completing the feedback cycle between these three important processes. The effect of wind--driven advection on the chemical composition is qualitatively similar to our previous results for the warmer atmosphere of HD~209458b, found using the same model. However, we find more significant alterations to both the thermal and dynamical structure for the cooler atmosphere of HD~189733b, with changes in both the temperature and wind velocities reaching $\sim10\%$. We also present the contribution function, diagnosed from our simulations, and show that wind--driven chemistry has a significant impact on its three--dimensional structure, particularly for regions where methane is an important absorber. Finally, we present emission phase curves from our simulations and show the significant effect of wind--driven chemistry on the thermal emission, particularly within the 3.6 \textmu m Spitzer/IRAC channel.

The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission.

We present results of simulations of the climate of the newly discovered planet Proxima Centauri B, performed using the Met Office Unified Model (UM). We examine the responses of both an `Earth-like' atmosphere and simplified nitrogen and trace carbon dioxide atmosphere to the radiation likely received by Proxima Centauri B. Additionally, we explore the effects of orbital eccentricity on the planetary conditions using a range of eccentricities guided by the observational constraints. Overall, our results are in agreement with previous studies in suggesting Proxima Centauri B may well have surface temperatures conducive to the presence of liquid water. Moreover, we have expanded the parameter regime over which the planet may support liquid water to higher values of eccentricity (>= 0.1) and lower incident fluxes (881.7 Wm-2) than previous work. This increased parameter space arises because of the low sensitivity of the planet to changes in stellar flux, a consequence of the stellar spectrum and orbital configuration. However, we also find interesting differences from previous simulations, such as cooler mean surface temperatures for the tidally-locked case. Finally, we have produced high resolution planetary emission and reflectance spectra, and highlight signatures of gases vital to the evolution of complex life on Earth (oxygen, ozone and carbon dioxide).

We present high resolution transmission spectra, calculated directly from a 3D radiative-hydrodynamics simulation that includes kinetic cloud formation, for HD209458b. We find that the high opacity of our vertically extensive cloud deck, composed of a large number density of sub-micron particles, flattens the transmission spectrum and obscures spectral features identified in observed data. We use the PandExo simulator to explore features of our HD209458b spectrum which may be detectable with the James Webb Space Telescope (JWST). We determine that an 8 - 12 micron absorption feature attributed to the mixed-composition, predominantly silicate cloud particles is a viable marker for the presence of cloud. Further calculations explore, and trends are identified with, variations in cloud opacity, composition heterogeneity and artificially scaled gravitational settling on the transmission spectrum. Principally, by varying the upper extent of our cloud decks, rainout is identified to be a key process for the dynamical atmospheres of hot-Jupiters and shown to dramatically alter the resulting spectrum. Our synthetic transmission spectra, obtained from the most complete, forward atmosphere simulations to--date, allow us to explore the model's ability to conform with observations. Such comparisons can provide insight into the physical processes either missing, or requiring improvement.

Broad absorption signatures from alkali metals, such as the sodium (Na I) and potassium (K I) resonance doublets, have long been predicted in the optical atmospheric spectra of cloud-free irradiated gas-giant exoplanets1,2,3. However, observations have only revealed the narrow cores of these features rather than the full pressure-broadened profiles4-6. Cloud and haze opacity at the day-night planetary terminator are considered responsible for obscuring the absorption-line wings, which hinders constraints on absolute atmospheric abundances7-9. Here we present an optical transmission spectrum for the 'hot-Saturn' WASP-96b obtained with the Very Large Telescope, which exhibits the complete pressure-broadened profile of the sodium absorption feature. The spectrum is in excellent agreement with cloud-free, solar-abundance models assuming chemical equilibrium. We are able to measure a precise, absolute sodium abundance of log\epsilon_Na=6.9+0.6-0.4, and use it as a proxy to the planet's atmospheric metallicity relative to the solar value (Z_p/Z_\star=2.3+8.9/--1.7). This result is consistent with the mass-metallicity trend observed for solar-system planets and exoplanets10-12.

To understand and compare the 3D atmospheric structure of HD 209458 b and HD 189733 b, focusing on the formation and distribution of cloud particles, as well as their feedback on the dynamics and thermal profile. We couple the 3D Met Office Unified Model (UM), including detailed treatments of atmospheric radiative transfer and dynamics, to a kinetic cloud formation scheme. The resulting model self--consistently solves for the formation of condensation seeds, surface growth and evaporation, gravitational settling and advection, cloud radiative feedback via absorption and, crucially, scattering. Fluxes directly obtained from the UM are used to produce synthetic SEDs and phase curves. Our simulations show extensive cloud formation in both planets, however, cooler temperatures in the HD 189733 b result in higher cloud particle number densities. Sub-micron particles are suspended by vertical flows leading to extensive upper-atmosphere cloud cover. A combination of meridional advection and efficient cloud formation in cooler high latitude regions, result in enhanced cloud coverage for latitudes > 30 degrees and leads to a zonally banded structure for all our simulations. The cloud bands extend around the entire planet(s), as the temperatures, even on the day side, remain below the condensation temperature of silicates and oxides. Therefore, our simulated optical phase curve for HD 209458 b shows no `offset', in contrast to observations. Efficient scattering by cloud results an atmospheric cooling of up to 250K, and an advection-driven fluctuating cloud opacity causes temporal variability in the thermal emission. The inclusion of this fundamental cloud-atmosphere radiative feedback leads to significant differences with approaches neglecting these physical elements and suggests both a note of caution of interpretations neglecting such cloud feedback and scattering, and merits further study.

We present a study of the effect of wind-driven advection on the chemical composition of hot Jupiter atmospheres using a fully-consistent 3D hydrodynamics, chemistry and radiative transfer code, the Met Office Unified Model (UM). Chemical modelling of exoplanet atmospheres has primarily been restricted to 1D models that cannot account for 3D dynamical processes. In this work we couple a chemical relaxation scheme to the UM to account for the chemical interconversion of methane and carbon monoxide. This is done consistently with the radiative transfer meaning that departures from chemical equilibrium are included in the heating rates (and emission) and hence complete the feedback between the dynamics, thermal structure and chemical composition. In this letter we simulate the well studied atmosphere of HD~209458b. We find that the combined effect of horizontal and vertical advection leads to an increase in the methane abundance by several orders of magnitude; directly opposite to the trend found in previous works. Our results demonstrate the need to include 3D effects when considering the chemistry of hot Jupiter atmospheres. We calculate transmission and emission spectra, as well as the emission phase curve, from our simulations. We conclude that gas-phase non-equilibrium chemistry is unlikely to explain the model-observation discrepancy in the 4.5\,{\textmu m} {\it Spitzer}/IRAC channel. However, we highlight other spectral regions, observable with the James Webb Space Telescope, where signatures of wind-driven chemistry are more prominant.

Previous studies have demonstrated that continental carbon-silicate weathering is important to the continued habitability of a terrestrial planet. Despite this, few studies have considered the influence of land on the climate of a tidally-locked planet. In this work we use the Met Office Unified Model, coupled to a land surface model, to investigate the climate effects of a continent located at the sub-stellar point. We choose to use the orbital and planetary parameters of Proxima Centauri B as a template, to allow comparison with the work of others. A region of the surface where $T_{\text{s}} > 273.15\,\text{K}$ is always retained, and previous conclusions on the habitability of Proxima Centauri B remain intact. We find that sub-stellar land causes global cooling, and increases day-night temperature contrasts by limiting heat redistribution. Furthermore, we find that sub-stellar land is able to introduce a regime change in the atmospheric circulation. Specifically, when a continent offset to the east of the sub-stellar point is introduced, we observe the formation of two mid-latitude counterrotating jets, and a substantially weakened equatorial superrotating jet.

In this work we have performed a series of simulations of the atmosphere of GJ~1214b assuming different metallicities using the Met Office Unified Model (UM). The UM is a general circulation model (GCM) that solves the deep, non-hydrostatic equations of motion and uses a flexible and accurate radiative transfer scheme, based on the two-stream and correlated-$k$ approximations, to calculate the heating rates. In this work we consistently couple a well-tested Gibbs energy minimisation scheme to solve for the chemical equilibrium abundances locally in each grid cell for a general set of elemental abundances, further improving the flexibility and accuracy of the model. As the metallicity of the atmosphere is increased we find significant changes in the dynamical and thermal structure, with subsequent implications for the simulated phase curve. The trends that we find are qualitatively consistent with previous works, though with quantitative differences. We investigate in detail the effect of increasing the metallicity by splitting the mechanism into constituents, involving the mean molecular weight, the heat capacity and the opacities. We find the opacity effect to be the dominant mechanism in altering the circulation and thermal structure. This result highlights the importance of accurately computing the opacities and radiative transfer in 3D GCMs.

The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a 2D steady state atmospheric circulation model, the validity of which is assessed by comparison to 3D calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occuring in the Earth's atmosphere or oceans. At depth, the vanishing heating flux forces the atmospheric structure to converge to a hotter adiabat than the one obtained with 1D calculations, implying a larger radius for the planet. Not only do the calculations reproduce the observed radius of HD209458b, but also the observed correlation between radius inflation and irradiation for transiting planets. Vertical advection of potential temperature induced by non uniform atmospheric heating thus provides a robust mechanism explaining the inflated radii of irradiated hot Jupiters.

We present highlights from a large set of simulations of a hot Jupiter atmosphere, nominally based on HD 209458b, aimed at exploring both the evolution of the deep atmosphere, and the acceleration of the zonal flow or jet. We find the occurrence of a super-rotating equatorial jet is robust to changes in various parameters, and over long timescales, even in the absence of strong inner or bottom boundary drag. This jet is diminished in one simulation only, where we strongly force the deep atmosphere equator-to-pole temperature gradient over long timescales. Finally, although the eddy momentum fluxes in our atmosphere show similarities with the proposed mechanism for accelerating jets on tidally-locked planets, the picture appears more complex. We present tentative evidence for a jet driven by a combination of eddy momentum transport and mean flow.

The correlated-k method is frequently used to speed up radiation calculations in both one-dimensional and three-dimensional atmosphere models. An inherent difficulty with this method is how to treat overlapping absorption, i.e. absorption by more than one gas in a given spectral region. We have evaluated the applicability of three different methods in hot Jupiter and brown dwarf atmosphere models, all of which have been previously applied within models in the literature: (i) Random overlap, both with and without resorting and rebinning, (ii) equivalent extinction and (iii) pre-mixing of opacities, where (i) and (ii) combine k-coefficients for different gases to obtain k-coefficients for a mixture of gases, while (iii) calculates k-coefficients for a given mixture from the corresponding mixed line-by-line opacities. We find that the random overlap method is the most accurate and flexible of these treatments, and is fast enough to be used in one-dimensional models with resorting and rebinning. In three-dimensional models such as GCMs it is too slow, however, and equivalent extinction can provide a speed-up of at least a factor of three with only a minor loss of accuracy while at the same time retaining the flexibility gained by combining k-coefficients computed for each gas individually. Pre-mixed opacities are significantly less flexible, and we also find that particular care must be taken when using this method in order to to adequately resolve steep variations in composition at important chemical equilibrium boundaries. We use the random overlap method with resorting and rebinning in our one-dimensional atmosphere model and equivalent extinction in our GCM, which allows us to e.g. consistently treat the feedback of non-equilibrium chemistry on the total opacity and therefore the calculated P-T profiles in our models.

To study the complexity of hot Jupiter atmospheres revealed by observations of increasing quality, we have adapted the UK Met Office global circulation model (GCM), the Unified Model (UM), to these exoplanets. The UM solves the full 3D Navier-Stokes equations with a height-varying gravity, avoiding the simplifications used in most GCMs currently applied to exoplanets. In this work we present the coupling of the UM dynamical core to an accurate radiation scheme based on the two-stream approximation and correlated-k method with state-of-the-art opacities from ExoMol. Our first application of this model is devoted to the extensively studied hot Jupiter HD 209458b. We have derived synthetic emission spectra and phase curves, and compare them to both previous models also based on state-of-the-art radiative transfer, and to observations. We find a reasonable a agreement between observations and both our day side emission and hot spot offset, however, our night side emission is too large. Overall our results are qualitatively similar to those found by Showman et al. (2009) with the SPARC/MITgcm, however, we note several quantitative differences: Our simulations show significant variation in the position of the hottest part of the atmosphere with pressure, as expected from simple timescale arguments, and in contrast to the "vertical coherency" found by Showman et al. (2009). We also see significant quantitative differences in calculated synthetic observations. Our comparisons strengthen the need for detailed intercomparisons of dynamical cores, radiation schemes and post-processing tools to understand these differences. This effort is necessary in order to make robust conclusions about these atmospheres based on GCM results.

In this work we investigate the impact of calculating non-equilibrium chemical abundances consistently with the temperature structure for the atmospheres of highly-irradiated, close-in gas giant exoplanets. Chemical kinetics models have been widely used in the literature to investigate the chemical compositions of hot Jupiter atmospheres which are expected to be driven away from chemical equilibrium via processes such as vertical mixing and photochemistry. All of these models have so far used pressure--temperature (P-T) profiles as fixed model input. This results in a decoupling of the chemistry from the radiative and thermal properties of the atmosphere, despite the fact that in nature they are intricately linked. We use a one-dimensional radiative-convective equilibrium model, ATMO, which includes a sophisticated chemistry scheme to calculate P-T profiles which are fully consistent with non-equilibrium chemical abundances, including vertical mixing and photochemistry. Our primary conclusion is that, in cases of strong chemical disequilibrium, consistent calculations can lead to differences in the P-T profile of up to 100 K compared to the P-T profile derived assuming chemical equilibrium. This temperature change can, in turn, have important consequences for the chemical abundances themselves as well as for the simulated emission spectra. In particular, we find that performing the chemical kinetics calculation consistently can reduce the overall impact of non-equilibrium chemistry on the observable emission spectrum of hot Jupiters. Simulated observations derived from non-consistent models could thus yield the wrong interpretation. We show that this behaviour is due to the non-consistent models violating the energy budget balance of the atmosphere.

3D atmosphere model results are used to comparatively study the kinetic, non-equilibrium cloud formation in the atmospheres of two example planets guided by the giant gas planets HD209458b and HD189733b. Rather independently of hydrodynamic model differences, our cloud modelling suggests that both planets are covered in mineral clouds throughout the entire modelling domain. Both planets harbour chemically complex clouds that are made of mineral particles that have a height-dependent material composition and size. The remaining gas-phase element abundances strongly effects the molecular abundances of the atmosphere in the cloud forming regions. Hydrocarbon and cyanopolyyne molecules can be rather abundant in the inner, dense part of the atmospheres of HD189733b and HD209458b. No one value for metallicity and the C/O ratio can be used to describe an extrasolar planet. Our results concerning the presence and location of water in relation to the clouds explain some of the observed differences between the two planets. In HD189733b, strong water features have been reported while such features appear less strong for HD209458b. By considering the location of the clouds in the two atmospheres, we see that obscuring clouds exist high in the atmosphere of HD209458b, but much deeper in HD189733b. We further conclude that the (self-imposed) degeneracy of cloud parameters in retrieval methods can only be lifted if the cloud formation processes are accurately modelled in contrast to prescribing them by independent parameters.

We present an isochrone server for semi-empirical pre-main-sequence model isochrones in the following systems: Johnson-Cousins, Sloan Digital Sky Survey, Two-Micron All-Sky Survey, Isaac Newton Telescope (INT) Wide-Field Camera, and INT Photometric H$\alpha$ Survey (IPHAS)/UV-Excess Survey (UVEX). The server can be accessed via the Cluster Collaboration webpage {http://www.astro.ex.ac.uk/people/timn/isochrones/}. To achieve this we have used the observed colours of member stars in young clusters with well-established age, distance and reddening to create fiducial loci in the colour-magnitude diagram. These empirical sequences have been used to quantify the discrepancy between the models and data arising from uncertainties in both the interior and atmospheric models, resulting in tables of semi-empirical bolometric corrections (BCs) in the various photometric systems. The model isochrones made available through the server are based on existing stellar interior models coupled with our newly derived semi-empirical BCs. As part of this analysis we also present new cluster parameters for both the Pleiades and Praesepe, yielding ages of $135^{+20}_{-11}$ and $665^{+14}_{-7}\,\rm{Myr}$ as well as distances of $132 \pm 2$ and $184 \pm 2\,\rm{pc}$ respectively (statistical uncertainty only).

The treatment of radiation transport in global circulation models (GCMs) is crucial to correctly describe Earth and exoplanet atmospheric dynamics processes. The two-stream approximation and correlated-k method are currently state-of-the-art approximations applied in both Earth and hot Jupiter GCM radiation schemes to facilitate rapid calculation of fluxes and heating rates. Their accuracy have been tested extensively for Earth-like conditions, but verification of the methods' applicability to hot Jupiter-like conditions is lacking in the literature. We are adapting the UK Met Office GCM, the Unified Model (UM), for the study of hot Jupiters, and present in this work the adaptation of the Edwards-Slingo radiation scheme based on the two-stream approximation and the correlated-k method. We discuss the calculation of absorption coefficients from high temperature line lists and highlight the large uncertainty in the pressure-broadened line widths. We compare fluxes and heating rates obtained with our adapted scheme to more accurate discrete ordinate (DO) line-by-line (LbL) calculations ignoring scattering effects. We find that, in most cases, errors stay below 10 % for both heating rates and fluxes using ~ 10 k-coefficients in each band and a diffusivity factor D = 1.66. The two-stream approximation and the correlated-k method both contribute non-negligibly to the total error. We also find that using band-averaged absorption coefficients, which have previously been used in radiative-hydrodynamical simulations of a hot Jupiter, may yield errors of ~ 100 %, and should thus be used with caution.

We are adapting the Global Circulation Model (GCM) of the UK Met Office, the so--called Unified Model (UM), for the study of hot Jupiters. In this work we demonstrate the successful adaptation of the most sophisticated dynamical core, the component of the GCM which solves the equations of motion for the atmosphere, available within the UM, ENDGame (Even Newer Dynamics for General atmospheric modelling of the environment). Within the same numerical scheme ENDGame supports solution to the dynamical equations under varying degrees of simplification. We present results from a simple, shallow (in atmospheric domain) hot Jupiter model (SHJ), and a more realistic (with a deeper atmosphere) HD 209458b test case. For both test cases we find that the large--scale, time--averaged (over the 1200 days prescribed test period), dynamical state of the atmosphere is relatively insensitive to the level of simplification of the dynamical equations. However, problems exist when attempting to reproduce the results for these test cases derived from other models. For the SHJ case the lower (and upper) boundary intersects the dominant dynamical features of the atmosphere meaning the results are heavily dependent on the boundary conditions. For the HD 209458b test case, when using the more complete dynamical models, the atmosphere is still clearly evolving after 1200 days, and in a transient state. Solving the complete (deep atmosphere and non--hydrostatic) dynamical equations allows exchange between the vertical and horizontal momentum of the atmosphere, via Coriolis and metric terms. Subsequently, interaction between the upper atmosphere and the deeper more slowly evolving (radiatively inactive) atmosphere significantly alters the results, and acts over timescales longer than 1200 days.

We demonstrate that both the current (New Dynamics), and next generation (ENDGame) dynamical cores of the UK Met Office global circulation model, the UM, reproduce consistently, the long-term, large-scale flows found in several published idealised tests. The cases presented are the Held-Suarez test, a simplified model of Earth (including a stratosphere), and a hypothetical tidally locked Earth. Furthermore, we show that using simplifications to the dynamical equations, which are expected to be justified for the physical domains and flow regimes we have studied, and which are supported by the ENDGame dynamical core, also produces matching long-term, large-scale flows. Finally, we present evidence for differences in the detail of the planetary flows and circulations resulting from improvements in the ENDGame formulation over New Dynamics.

Spectropolarimetric observations combined with tomographic imaging techniques have revealed that all pre-main sequence (PMS) stars host multipolar magnetic fields, ranging from strong and globally axisymmetric with ~>kilo-Gauss dipole components, to complex and non-axisymmetric with weak dipole components (<~0.1 kG). Many host dominantly octupolar large-scale fields. We argue that the large-scale magnetic properties of a PMS star are related to its location in the Hertzsprung-Russell diagram. This conference paper is a synopsis of Gregory et al. (2012), updated to include the latest results from magnetic mapping studies of PMS stars.

We present a deep Cousins RI photometric survey of the open cluster NGC 1960, complete to R_C \simeq 22, I_C \simeq 21, that is used to select a sample of very low-mass cluster candidates. Gemini spectroscopy of a subset of these is used to confirm membership and locate the age-dependent "lithium depletion boundary" (LDB) --the luminosity at which lithium remains unburned in its low-mass stars. The LDB implies a cluster age of 22 +/-4 Myr and is quite insensitive to choice of evolutionary model. NGC 1960 is the youngest cluster for which a LDB age has been estimated and possesses a well populated upper main sequence and a rich low-mass pre-main sequence. The LDB age determined here agrees well with precise age estimates made for the same cluster based on isochrone fits to its high- and low-mass populations. The concordance between these three age estimation techniques, that rely on different facets of stellar astrophysics at very different masses, is an important step towards calibrating the absolute ages of young open clusters and lends confidence to ages determined using any one of them.

We have derived ages for 13 young (<30 Myr) star-forming regions and find they are up to a factor two older than the ages typically adopted in the literature. This result has wide-ranging implications, including that circumstellar discs survive longer (~10-12 Myr) and that the average Class I lifetime is greater (~1 Myr) than currently believed. For each star-forming region we derived two ages from colour-magnitude diagrams. First we fitted models of the evolution between the zero-age main-sequence and terminal-age main-sequence to derive a homogeneous set of main-sequence ages, distances and reddenings with statistically meaningful uncertainties. Our second age for each star-forming region was derived by fitting pre-main-sequence stars to new semi-empirical model isochrones. For the first time (for a set of clusters younger than 50 Myr) we find broad agreement between these two ages, and since these are derived from two distinct mass regimes that rely on different aspects of stellar physics, it gives us confidence in the new age scale. This agreement is largely due to our adoption of empirical colour-Teff relations and bolometric corrections for pre-main-sequence stars cooler than 4000 K. The revised ages for the star-forming regions in our sample are: ~2 Myr for NGC 6611 (Eagle Nebula; M 16), IC 5146 (Cocoon Nebula), NGC 6530 (Lagoon Nebula; M 8), and NGC 2244 (Rosette Nebula); ~6 Myr for {\sigma} Ori, Cep OB3b, and IC 348; ~10 Myr for {\lambda} Ori (Collinder 69); ~11 Myr for NGC 2169; ~12 Myr for NGC 2362; ~13 Myr for NGC 7160; ~14 Myr for {\chi} Per (NGC 884); and ~20 Myr for NGC 1960 (M 36).

We have used medium resolution spectra to search for evidence that proto-stellar objects accrete at high rates during their early 'assembly phase'. Models predict that depleted lithium and reduced luminosity in T-Tauri stars are key signatures of 'cold' high-rate accretion occurring early in a star's evolution. We found no evidence in 168 stars in NGC 2264 and the Orion Nebula Cluster for strong lithium depletion through analysis of veiling corrected 6708 angstrom lithium spectral line strengths. This suggests that 'cold' accretion at high rates (M_dot > 5 x 10-4 M_sol yr-1) occurs in the assembly phase of fewer than 0.5 per cent of 0.3 < M < 1.9 M_sol stars. We also find that the dispersion in the strength of the 6708 angstrom lithium line might imply an age spread that is similar in magnitude to the apparent age spread implied by the luminosity dispersion seen in colour magnitude diagrams. Evidence for weak lithium depletion (< 10 per cent in equivalent width) that is correlated with luminosity is also apparent, but we are unable to determine whether age spreads or accretion at rates less than 5 x 10-4 M_sol yr-1 are responsible.

ZDI studies have shown that the magnetic fields of T Tauri stars can be significantly more complex than a simple dipole and can vary markedly between sources. We collect and summarize the magnetic field topology information obtained to date and present Hertzsprung-Russell (HR) diagrams for the stars in the sample. Intriguingly, the large scale field topology of a given pre-main sequence (PMS) star is strongly dependent upon the stellar internal structure, with the strength of the dipole component of its multipolar magnetic field decaying rapidly with the development of a radiative core. Using the observational data as a basis, we argue that the general characteristics of the global magnetic field of a PMS star can be determined from its position in the HR diagram. Moving from hotter and more luminous to cooler and less luminous stars across the PMS of the HR diagram, we present evidence for four distinct magnetic topology regimes. Stars with large radiative cores, empirically estimated to be those with a core mass in excess of ~40 per cent of the stellar mass, host highly complex and dominantly non-axisymmetric magnetic fields, while those with smaller radiative cores host axisymmetric fields with field modes of higher order than the dipole dominant (typically, but not always, the octupole). Fully convective stars stars above ~0.5 MSun appear to host dominantly axisymmetric fields with strong (kilo-Gauss) dipole components. Based on similarities between the magnetic properties of PMS stars and main sequence M-dwarfs with similar internal structures, we speculate that a bistable dynamo process operates for lower mass stars (<~0.5 MSun at an age of a few Myr) and that they will be found to host a variety of magnetic field topologies. If the magnetic topology trends across the HR diagram are confirmed they may provide a new method of constraining PMS stellar evolution models.

We present a critical assessment of commonly used pre-main-sequence isochrones by comparing their predictions to a set of well-calibrated colour-magnitude diagrams of the Pleiades in the wavelength range 0.4 to 2.5 microns. Our analysis shows that for temperatures less than 4000 K the models systematically overestimate the flux by a factor two at 0.5 microns, though this decreases with wavelength, becoming negligible at 2.2 microns. In optical colours this will result in the ages for stars younger than 10 Myr being underestimated by factors between two and three. We show that using observations of standard stars to transform the data into a standard system can introduce significant errors in the positioning of pre-main-sequences in colour-magnitude diagrams. Therefore we have compared the models to the data in the natural photometric system in which the observations were taken. Thus we have constructed and tested a model of the system responses for the Wide-Field Camera on the Isaac Newton Telescope. As a benchmark test for the development of pre-main-sequence models we provide both our system responses and the Pleiades sequence.

We present derived stellar and disc parameters for a sample of Taurus brown dwarfs both with and without evidence of an associated disc. These parameters have been derived using an online fitting tool (http://bd-server.astro.ex.ac.uk/), which includes a statistically robust derivation of uncertainties, an indication of pa- rameter degeneracies, and a complete treatment of the input photometric and spectroscopic observations. The observations of the Taurus members with indications of disc presence have been fitted using a grid of theoretical models including detailed treatments of physical processes accepted for higher mass stars, such as dust sublimation, and a simple treatment of the accretion flux. This grid of models has been designed to test the validity of the adopted physical mechanisms, but we have also constructed models using parameterisation, for example semi-empirical dust sublimation radii, for users solely interested in parameter derivation and the quality of the fit. The parameters derived for the naked and disc brown dwarf systems are largely consistent with literature observations. However, our inner disc edge locations are consistently closer to the star than previous results and we also derive elevated accretion rates over non-SED based accretion rate derivations. For inner edge locations we attribute these differences to the detailed modelling we have performed of the disc structure, particularly at the crucial inner edge where departures in geometry from the often adopted vertical wall due to dust sublimation (and therefore accretion flux) can compensate for temperature (and therefore distance) changes to the inner edge of the dust disc. In the case of the elevated derived accretion rates, in some cases, this may be caused by the intrinsic stellar luminosities of the targets exceeding that predicted by the isochrones we have adopted.

The wide luminosity dispersion seen for stars at a given effective temperature in the H-R diagrams of young clusters and star forming regions is often interpreted as due to significant (~10 Myr) spreads in stellar contraction age. In the scenario where most stars are born with circumstellar discs, and that disc signatures decay monotonically (on average) over timescales of only a few Myr, then any such age spread should lead to clear differences in the age distributions of stars with and without discs. We have investigated large samples of stars in the Orion Nebula Cluster (ONC) using three methods to diagnose disc presence from infrared measurements. We find no significant difference in the mean ages or age distributions of stars with and without discs, consistent with expectations for a coeval population. Using a simple quantitative model we show that any real age spread must be smaller than the median disc lifetime. For a log-normal age distribution, there is an upper limit of <0.14 dex (at 99% confidence) to any real age dispersion, compared to the ~=0.4 dex implied by the H-R diagram. If the mean age of the ONC is 2.5 Myr, this would mean at least 95% of its low-mass stellar population has ages between 1.3--4.8 Myr. We suggest that the observed luminosity dispersion is caused by a combination of observational uncertainties and physical mechanisms that disorder the conventional relationship between luminosity and age for pre main-sequence stars. This means that individual stellar ages from the H-R diagram are unreliable and cannot be used to directly infer a star formation history. Irrespective of what causes the wide luminosity dispersion, the finding that any real age dispersion is less than the median disc lifetime argues strongly against star formation scenarios for the ONC lasting longer than a few Myr.

We present an analysis of the rotation of young stars in the associations Cepheus OB3b, NGC 2264, NGC 2362 and the Orion Nebula Cluster (ONC). We discover a correlation between rotation rate and position in a colour-magnitude diagram (CMD) such that stars which lie above an empirically determined median pre-main sequence rotate more rapidly than stars which lie below this sequence. The same correlation is seen, with a high degree of statistical significance, in each association studied here. If position within the CMD is interpreted as being due to genuine age spreads within a cluster, then the stars above the median pre-main sequence would be the youngest stars. This would in turn imply that the most rapidly rotating stars in an association are the youngest, and hence those with the largest moments of inertia and highest likelihood of ongoing accretion. Such a result does not fit naturally into the existing picture of angular momentum evolution in young stars, where the stars are braked effectively by their accretion discs until the disc disperses. Instead, we argue that, for a given association of young stars, position within the CMD is not primarily a function of age, but of accretion history. We show that this hypothesis could explain the correlation we observe between rotation rate and position within the CMD.

We present a grid of models of accreting brown dwarf systems with circumstellar discs. The calculations involve a self-consistent solution of both vertical hydrostatic and radiative equilibrium along with a sophisticated treatment of dust sublimation. We have simulated observations of the spectral energy distributions and several broadband photometric systems. Analysis of the disc structures and simulated observations reveal a natural dichotomy in accretion rates, with \logmdot $>-$9 and $\leq -$9 classed as extreme and typical accretors respectively. Derivation of ages and masses from our simulated photometry using isochrones is demonstrated to be unreliable even for typical accretors. Although current brown dwarf disc candidate selection criteria have been shown to be largely reliable when applied to our model grid we suggest improved selection criteria in several colour indices. We show that as accretion rates increase brown dwarf disc systems are less likely to be correctly identified. This suggests that, within our grid, systems with higher accretion rates would be preferentially lost during brown dwarf target selection. We suggest that observations used to assert a $\dot{M}\propto M_*^2$ relationship may contain an intrinsic selection bias.

We present a discussion of the similarities and key differences between the transition onto (at the turn-on) and away from (at the turn-off) the main sequence, the latter termed the Hertzsprung gap. Using a set of model isochrones and adopting an initial mass function leads us to predict a dearth of G-type stars for any star forming region. This is caused by the (relatively) constant spectral type at which the transition from a fully convective star to a star with a radiative core begins. We also present analysis of the details of this transition in the ONC. In particular we show that a gap in the photometric and spectral type distributions is centred on, and a change in the fractional X-ray luminosity and rotation rate distribution occurs approximately at, the position of a peak in radiative core size as a function of mass. Whilst photometric signatures of this transition are lost at ages over ~20 Myrs, we show that changes in fractional X-ray luminosity and magnetic field configuration persist to older ages. Analysis of literature data show that the mass at which the change in fractional X-ray luminosity is observed decreases with age.

We present a photometric study of I-band variability in the young association Cepheus OB3b. The study is sensitive to periodic variability on timescales of less than a day, to more than 20 days. After rejection of contaminating objects using V, I, R and narrowband H-alpha photometry, we find 475 objects with measured rotation periods, which are very likely pre-main-sequence members of the Cep OB3b star forming region. We revise the distance and age to Cep OB3b, putting it on the self-consistent age and distance ladder of Mayne & Naylor (2008). This yields a distance modulus of 8.8 +/- 0.2 mags, corresponding to a distance of 580 +/- 60 pc, and an age of 4-5Myrs. The rotation period distribution confirms the general picture of rotational evolution in young stars, exhibiting both the correlation between accretion (determined in this case through narrowband H-alpha photometry) and rotation expected from disc locking, and the dependence of rotation upon mass that is seen in other star forming regions. However, this mass dependence is much weaker in our data than found in other studies. Comparison to the similarly aged NGC 2362 shows that the low-mass stars in Cep OB3b are rotating much more slowly. This points to a possible link between star forming environment and rotation properties. Such a link would call into question models of stellar angular momentum evolution, which assume that the rotational period distributions of young clusters and associations can be assembled into an evolutionary sequence, thus ignoring environmental effects.

We present a consistent age ordering for young clusters and groups determined using the contraction of stars through their pre-main-sequence phase. We compare these with ages derived from the evolution of the upper main-sequence stars, and find the upper MS ages are older by a factor 1.5 to 2. We show that increasing the binary fraction and number of equal-mass binaries amongst the O-stars compared to the rest of the MS cannot remove this discrepancy.

We present three new methods for determining the age of groups of pre-main-sequence stars. The first, creating empirical isochrones allows us to create a robust age ordering, but not to derive actual ages. The second, using the width of the gap in colour-magnitude space between the pre-main-sequence and main-sequence (the radiative convective gap) has promise as a distance and extinction independent measure of age, but is as yet uncalibrated. Finally we discuss tau-squared fitting of the main sequence as the stars approach the terminus of the main sequence. This method suggests that there is a factor two difference between these "nuclear" ages, and more conventional pre-main-sequence contraction ages.

We use I band imaging to perform a variability survey of the 13 Myr-old cluster h Per. We find a significant fraction of the cluster members to be variable. Most importantly, we find that variable members lie almost entirely on the convective side of the gap in the cluster sequence between fully convective stars and those which have a radiative core. This result is consistent with a scenario in which the magnetic field changes topology when the star changes from being fully convective, to one containing a radiative core. When the star is convective the magnetic field appears dominated by large-scale structures, resulting in global-size spots that drive the observed variability. For those stars with radiative cores we observe a marked absence of variability due to spots, which suggests a switch to a magnetic field dominated by smaller-scale structures, resulting in many smaller spots and thus less apparent variability. This implies that wide field variability surveys may only be sensitive to fully convective stars. On the one hand this reduces the chances of picking out young groups (since the convective stars are the lower mass and therefore fainter objects), but conversely the absolute magnitude of the head of the convective sequence provides a straightforward measure of age for those groups which are discovered.

We use several main-sequence models to derive distances (and extinctions), with statistically meaningful uncertainties for 11 star-forming-regions and young clusters. The model dependency is shown to be small, allowing us to adopt the distances derived using one model. Using these distances we have revised the age order for some of the clusters of Mayne et al (2007). The new (approximate) nominal ages are: 2 Myrs for NGC6530 and the ONC, 3 Myrs for Lambda Orionis, NGC2264 and Sigma Orionis, 4-5 Myrs for NGC2362, 13 Myrs for h and chi Per, 20 Myrs for NGC1960 and 40 Myrs for NGC2547. In cases of significantly variable extinction we have derived individual extinctions using a revised Q-method (Johnson and Morgan, 1953). These new data show that the largest remaining uncertainty in deriving an age ordering (and necessarily ages) is metallicity. We also discuss the use of a feature we term the R-C gap overlap to provide a diagnostic of isochronal age spreads or varying accretion histories within a given star-formation-region. Finally, recent derivations of the distance to the ONC lie in two groups. Our new more precise distance of 391(+12,-9) pc allows us to decisively reject the further distance, we adopt 400 pc as a convenient value.

We have selected pre-main-sequence stars in 12 groups of notional ages ranging from 1 Myr to 35 Myrs, using heterogeneous membership criteria. Using these members we have constructed empirical isochrones in V, V-I colour magnitude diagrams (CMDs). This allows us to identify clearly the gap between the radiative main sequence and the convective pre-main-sequence (the R-C gap). We follow the evolution of this gap with age and show that it can be a useful age indicator for groups less than ~15 Myrs old. We also observe a reduction in absolute spreads about the sequences with age. Finally the empirical isochrones allow us to place the groups in order of age, independently of theory. The youngest groups can be collated into three sets of similar ages. The youngest set is the ONC, NGC6530 and IC5146 (nominally 1 Myrs); next Cep OB3b, NGC2362, Lambda Ori and NGC2264 (nominally 3 Myrs); and finally Sigma Ori and IC348 (nominally 4-5 Myrs). This suggests Cep OB3b is younger than previously thought, and IC348 older. For IC348 the stellar rotation rate distribution and fraction of stars with discs imply a younger age than we derive. We suggest this is because of the absence of O-stars in this cluster, whose winds and/or ionising radiation may be an important factor in the removal of discs in other clusters.

We have used R_C I_C CCD photometry from the Isaac Newton telescope and intermediate resolution spectroscopy from the Gemini North telescope to identify and characterise low-mass (0.15<M/M_sun<1.3) pre-main sequence stars in the young open cluster NGC 2169. Isochrone fitting to the high- and low-mass populations yields an intrinsic distance modulus of 10.13{+0.06}{-0.09} mag and a model-dependent age of 9+/-2 Myr. Compared with the nearby, kinematically defined groups of a similar age, NGC 2169 has a large low-mass population which potentially offers a more precise statistical investigation of several aspects of star formation and early stellar evolution. By modelling the distribution of low-mass stars in the I_C versus R_C-I_C diagram we find that any age spread among cluster members has a Gaussian full width at half maximum <=2.5 Myr. A young age and small age spread (<10 Myr) are supported by the lack of significant lithium depletion in the vast majority of cluster members. There is no clear evidence for accretion or warm circumstellar dust in the low-mass members of NGC 2169, bolstering the idea that strong accretion has ceased and inner discs have dispersed in almost all low-mass stars by ages of 10 Myr.