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

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.

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.

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.

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.

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.

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.

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.

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 exploit deep combined observations with Spitzer and Chandra of the SWIRE survey in the ELAIS-N1 region, to investigate the nature of the faint X-ray and IR sources in common, to identify AGN/starburst diagnostics, and to study the sources of the X-ray and IR cosmic backgrounds. In the 17'x17' area of the Chandra ACIS-I image there are 3400 SWIRE near-IR sources with 4 sigma detections in at least 2 IRAC bands and 988 sources detected at 24micron with MIPS brighter than 0.1 mJy. Of these, 102 IRAC and 59 MIPS sources have Chandra counterparts, out of a total of 122 X-ray sources present in the area with S(0.5-8 kev)>10^(-15) erg/cm^2/s. We have constructed SEDs for each source using data from the 4 IRAC wavebands, Chandra fluxes, and optical follow-up data in the wavebands U, g', r', i', Z, and H. We fit a number of spectral templates to the SEDs at optical and infrared wavelengths to determine photometric redshifts and spectral categories, and also make use of diagnostics based on the X-ray luminosities, hardness ratios, X-ray to infrared spectral slopes and optical morphologies. Although we have spectroscopic redshifts for only a minority of the Chandra sources, the available SEDs constrain the redshifts for most of the sample sources, which turn out to be typically at 0.5<z<2. We find that 39% of the Chandra sources are dominated by type-1 AGN emission, 23% display optical/IR spectra typical of type-2 AGNs, while the remaining 38% fraction show starburst-like or even normal galaxy spectra. Since we prove that all these galaxies are dominated by AGN emission in X-rays this brings the fraction of type-1 AGNs to be 80% of the type-2: even assuming that all the Chandra sources undetected by Spitzer are type-2 AGNs, the type-1 fraction would exceed 1/3 of the total population (abridged).

We present a cross-correlation of the European Large Area ISO survey (ELAIS) with the ELAIS Deep X-ray Survey of the N1 and N2 fields. There are 7 Chandra point sources with matches in the ELAIS Final Analysis 15um catalogue, out of a total of 28 extragalactic ISO sources present in the Chandra fields. Five of these are consistent with AGN giving an AGN fraction of ~19 per cent in the 15um flux range 0.8-6 mJy. We have co-added the hard X-ray fluxes of the individually-undetected ISO sources and find a low significance detection consistent with star formation in the remaining population. We combine our point source cross-correlation fraction with the XMM-Newton observations of the Lockman Hole and Chandra observations of the Hubble Deep Field North to constrain source count models of the mid-infrared galaxy population. The low dust-enshrouded AGN fraction in ELAIS implied by the number of cross-identifications between the ELAIS mid-infrared sample and the Chandra point sources is encouraging for the use of mid-infrared surveys to constrain the cosmic star formation history, provided there are not further large undetected populations of Compton-thick AGN.

Hard X-ray surveys have uncovered a large population of heavily obscured AGN. They also reveal a population of quasars with moderate obscuration at both visible and X-ray wavelengths. We use Chandra selected samples of quasars from the ELAIS Deep X-ray Survey (EDXS) and CDF-N to investigate the obscuration towards the nuclei of moderately obscured AGN. We find an inverse correlation between the optical to X-ray flux ratio and the X-ray hardness ratio which can be interpreted as due to obscuration at visible and X-ray wavelengths. We present detailed optical and near-infrared data for a sample of optically-faint (R>23) quasars from the EDXS. These are used to constrain the amount of rest-frame UV/optical reddening towards these quasars. It is found that optically-faint quasars are mostly faint due to obscuration, not because they are intrinsically weak. After correcting for reddening, the optical magnitudes of most of these quasars are similar to the brighter quasars at these X-ray fluxes. Combining with gas column densities inferred from the X-ray observations we consider the gas-to-dust ratios of the obscuring matter. We find that the quasars generally have higher gas-to-dust absorption than that seen in the Milky Way - similar to what has been found for nearby Seyfert galaxies. We consider the possible existence of a large population of X-ray sources which have optical properties of Type 1 (unobscured) quasars, but X-ray properties of Type 2 (obscured) quasars. We show that the observed distribution of optical-to-X-ray flux ratios of quasars at z>1 is skewed to low values compared to the intrinsic distribution...

We present an analysis of two deep (75 ks) Chandra observations of the European Large Area ISO Survey (ELAIS) fields N1 and N2 as the first results from the ELAIS deep X-ray survey. This survey is being conducted in well studied regions with extensive multi-wavelength coverage. Here we present the Chandra source catalogues along with an analysis of source counts, hardness ratios and optical classifications. A total of 233 X-ray point sources are detected in addition to 2 soft extended sources, which are found to be associated with galaxy clusters. An over-density of sources is found in N1 with 30% more sources than N2, which we attribute to large-scale structure. A similar variance is seen between other deep Chandra surveys. The source count statistics reveal an increasing fraction of hard sources at fainter fluxes. The number of galaxy-like counterparts also increases dramatically towards fainter fluxes, consistent with the emergence of a large population of obscured sources.

The sources discovered in deep hard X-ray surveys with 2-8 keV fluxes of 10^-14 erg cm^-2 s^-1 make up the bulk of the X-ray background at these energies. We present here detailed multi-wavelength observations of three such sources from the ELAIS Deep X-ray Survey. The observations include sensitive near-infrared spectroscopy with the Subaru Telescope and X-ray spectral information from the Chandra X-ray Observatory. The sources observed all have optical-to-near-IR colours redder than an unobscured quasar and comprise a reddened quasar, a radio galaxy and an optically-obscured AGN. The reddened quasar is at a redshift z=2.61 and shows a very large X-ray absorbing column of N_H approx 3.10^23 cm^-2. This contrasts with the relatively small amount of dust reddening, implying a gas-to-dust ratio along the line-of-sight a hundred times greater than that of the Milky Way. The radio galaxy at z=1.57 shows only narrow emission lines, but has a surprisingly soft X-ray spectrum. The softness of this spectrum either indicates an unusually low gas-to-dust ratio for the absorbing medium or X-ray emission related to the young radio source. The host galaxy is extremely red (R-K=6.4) and its optical/near-IR spectrum is best fit by a strongly reddened (A_V~2) starburst. The third X-ray source discussed is also extremely red (R-K=6.1) and lies in a close grouping of three other R-K>6 galaxies. No emission or absorption lines were detected from this object, but its redshift (and that of one of the nearby galaxies) are constrained by SED-fitting to be just greater than z=1. The extremely red colours of these two galaxies can be accounted for by old stellar populations. These observations illustrate the diverse properties of hard X-ray selected AGN.

We explore the relationship between the hard X-ray and sub-mm populations using deep Chandra observations of a large, contiguous SCUBA survey. In agreement with other recent findings, we confirm that the direct overlap is small. Of the 17 sub-mm sources detected in this field at 850 microns, only one is coincident with a Chandra source. The resulting limits imply that the majority of SCUBA sources are not powered by AGN, unless the central engine is obscured by Compton-thick material with a low (<1 per cent) scattered component. Furthermore, since Chandra detects only ~5 per cent of SCUBA sources, the typical obscuration would need to be almost isotropic. The X-ray upper limits are so strong that in most cases we can also rule out a starburst SED at low redshift, suggesting that the majority of SCUBA sources lie at z>1 even if they are purely starburst galaxies. Despite the low detection rate, we find evidence for strong angular clustering between the X-ray and sub-mm populations. The implication is that AGN and SCUBA sources trace the same large-scale structure but do not generally coincide. If bright sub-mm sources represent massive elliptical galaxies in formation, we suggest that (for a given galaxy) the major episode of star-formation must be distinct from the period of observable quasar activity.

We present an analysis of X-ray variability in a sample of 156 radio quiet quasars taken from the ROSAT archive, covering a redshift range 0.1 < z < 4.1. A maximum likelihood method is used to constrain the amplitude of variability in these low signal to noise light curves. Through combining these data in ensembles we are able to identify trends in variability amplitude with luminosity and with redshift. The decline in variability amplitude with luminosity identified in local AGN (z < 0.1) is confirmed out to z=2. There is tentative evidence for an increase in QSO X-ray variability amplitude towards high redshifts (z > 2) in the sense that QSOs of the same X-ray luminosity are more variable at z > 2. We discuss possible explanations for this effect. The simplest explanation may be that high redshift QSOs are accreting at a larger fraction of the Eddington limit than local AGN.

We present initial follow-up results of the ELAIS Deep X-ray Survey which is being undertaken with the Chandra and XMM-Newton Observatories. 235 X-ray sources are detected in our two 75 ks ACIS-I observations in the well-studied ELAIS N1 and N2 areas. 90% of the X-ray sources are identified optically to R=26 with a median magnitude of R=24. We show that objects which are unresolved optically (i.e. quasars) follow a correlation between their optical and X-ray fluxes, whereas galaxies do not. We also find that the quasars with fainter optical counterparts have harder X-ray spectra, consistent with absorption at both wavebands. Initial spectroscopic follow-up has revealed a large fraction of high-luminosity Type 2 quasars. The prospects for studying the evolution of the host galaxies of X-ray selected Type 2 AGN are considered.

We present multiwaveband photometric and optical spectropolarimetric observations of the R=15.9 narrow emission line galaxy R117_A which lies on the edge of the error circle of the ROSAT X-ray source R117 (from McHardy et al 1998). The overall spectral energy distribution of the galaxy is well modelled by a combination of a normal spiral galaxy and a moderate-strength burst of star formation. The far infra-red and radio emission is extended along the major axis of the galaxy, indicating an extended starburst. On positional grounds, the galaxy is a good candidate for the identification of R117 and the observed X-ray flux is very close to what would be expected from a starburst of the observed far infra-red and radio fluxes. Although an obscured high redshift QSO cannot be entirely ruled out as contributing some fraction of the X-ray flux, we find no candidates to K=20.8 within the X-ray errorbox and so conclude that R117_A is responsible for a large fraction, if not all, of the X-ray emission from R117. Searches for indicators of an obscured AGN in R117_A have so far proven negative; deep spectropolarimetric observations show no signs of broad lines to a limit of one per cent and, for the observed far infra-red and radio emission, we would expect a ten times greater X-ray flux if the overall emission were powered by an AGN. We therefore conclude that the X-ray emission from R117 is dominated by starburst emission from the galaxy R117_A.