Over the last several years, spectroscopic observations of transiting exoplanets have begun to uncover information about their atmospheres, including atmospheric composition and indications of the presence of clouds and hazes. Spectral retrieval is the leading technique for interpretation of transmission spectra and is employed by several teams using a variety of forward models and parameter estimation algorithms. However, different model suites have mostly been used in isolation and so it is unknown whether the results from each are comparable. As we approach the launch of the James Webb Space Telescope we anticipate advances in wavelength coverage, precision, and resolution of transit spectroscopic data, so it is important that the tools that will be used to interpret these information rich spectra are validated. To this end, we present an inter-model comparison of three retrieval suites: TauREx, NEMESIS and CHIMERA. We demonstrate that the forward model spectra are in good agreement (residual deviations on the order of 20 - 40 ppm), and discuss the results of cross retrievals between the three tools. Generally, the constraints from the cross retrievals are consistent with each other and with input values to within 1 sigma However, for high precision scenarios with error envelopes of order 30 ppm, subtle differences in the simulated spectra result in discrepancies between the different retrieval suites, and inaccuracies in retrieved values of several sigma. This can be considered analogous to substantial systematic/astrophysical noise in a real observation, or errors/omissions in a forward model such as molecular linelist incompleteness or missing absorbers.

We present a novel retrieval technique that attempts to model phase curve observations of exoplanets more realistically and reliably, which we call the 2.5-dimension (2.5-D) approach. In our 2.5-D approach we retrieve the vertical temperature profile and mean gaseous abundance of a planet at all longitudes and latitudes \textbf{simultaneously}, assuming that the temperature or composition, $x$, at a particular longitude and latitude $(\Lambda,\Phi)$ is given by $x(\Lambda,\Phi) = \bar{x} + (x(\Lambda,0) - \bar{x})\cos^n\Phi$, where $\bar{x}$ is the mean of the morning and evening terminator values of $x(\Lambda,0)$, and $n$ is an assumed coefficient. We compare our new 2.5-D scheme with the more traditional 1-D approach, which assumes the same temperature profile and gaseous abundances at all points on the visible disc of a planet for each individual phase observation, using a set of synthetic phase curves generated from a GCM-based simulation. We find that our 2.5-D model fits these data more realistically than the 1-D approach, confining the hotter regions of the planet more closely to the dayside. We then apply both models to WASP-43b phase curve observations of HST/WFC3 and Spitzer/IRAC. We find that the dayside of WASP-43b is apparently much hotter than the nightside and show that this could be explained by the presence of a thick cloud on the nightside with a cloud top at pressure $< 0.2$ bar. We further show that while the mole fraction of water vapour is reasonably well constrained to $(1-10)\times10^{-4}$, the abundance of CO is very difficult to constrain with these data since it is degenerate with temperature and prone to possible systematic radiometric differences between the HST/WFC3 and Spitzer/IRAC observations. Hence, it is difficult to reliably constrain C/O.

The current sparse wavelength range coverage of exoplanet direct imaging observations, and the fact that models are defined using a finite wavelength range, lead both to uncertainties on effective temperature determination.We study these effects using black-bodies and atmospheric models and we detail how to infer this parameter. Through highlighting the key wavelength coverage that allows for a more accurate representation of the effective temperature, our analysis can be used to mitigate or manage extra uncertainties being added in the analysis from the models. We find that the wavelength range coverage will soon no longer be a problem. An effective temperature computed by integrating the spectroscopic observations of the James Webb Space Telescope (JWST) will give uncertainties similar to, or better than, the current state-of-the-art, which is to fit models to data. Accurately calculating the effective temperature will help to improve current modelling approaches. Obtaining an independent and precise estimation of this crucial parameter will help the benchmarking process to identify the best practice to model exoplanet atmospheres.

Current observational data of exoplanets are providing increasing detail of their 3D atmospheric structures. As characterisation efforts expand in scope, the need to develop consistent 3D radiative-transfer methods becomes more pertinent as the complex atmospheric properties of exoplanets are required to be modelled together consistently. We aim to compare the transmission and emission spectra results of a 3D Monte Carlo Radiative Transfer (MCRT) model to contemporary radiative-transfer suites. We perform several benchmarking tests of a MCRT code, Cloudy Monte Carlo Radiative Transfer (CMCRT), to transmission and emission spectra model output. We add flexibility to the model through the use of k-distribution tables as input opacities. We present a hybrid MCRT and ray tracing methodology for the calculation of transmission spectra with a multiple scattering component. CMCRT compares well to the transmission spectra benchmarks at the 10s of ppm level. Emission spectra benchmarks are consistent to within 10% of the 1D models. We suggest that differences in the benchmark results are likely caused by geometric effects between plane-parallel and spherical models. In a practical application, we post-process a cloudy 3DHD 189733b GCM model and compare to available observational data. Our results suggest the core methodology and algorithms of CMCRT produce consistent results to contemporary radiative transfer suites. 3D MCRT methods are highly suitable for detailed post-processing of cloudy and non-cloudy 1D and 3D exoplanet atmosphere simulations in instances where atmospheric inhomogeneities, significant limb effects/geometry or multiple scattering components are important considerations.

Sophisticated atmospheric retrieval algorithms, such as Nested Sampling, explore large parameter spaces by iterating over millions of radiative transfer (RT) calculations. Probability distribution functions for retrieved parameters are highly sensitive to assumptions made within the RT forward model. One key difference between RT models is the computation of the gaseous absorption throughout the atmosphere. We compare two methods of calculating gaseous absorption, cross-sections and correlated-$k$, by examining their resulting spectra of a number of typical \ce{H2}-He dominated exoplanetary and brown dwarf atmospheres. We also consider the effects of including \ce{H2}-He pressure-broadening in some of these examples. We use NEMESIS to compute forward models. Our $k$-tables are verified by comparison to ExoMol cross-sections provided online and a line-by-line calculation. For test cases with typical resolutions ($\Delta \nu = 1$cm$^{-1}$), we show that the cross-section method overestimates the amount of absorption present in the atmosphere and should be used with caution. For mixed-gas atmospheres the morphology of the spectra changes, producing `ghost' features. The two methods produce differences in flux of up to a few orders of magnitude. The addition of pressure broadening of lines adds up to an additional order of magnitude change in flux. These effects are more pronounced for brown dwarfs and secondary eclipse geometries. We note that correlated-$k$ can produce similar results to very high-resolution cross-sections, but is much less computationally expensive. We conclude that inaccurate use of cross-sections and omission of pressure broadening can be key sources of error in the modelling of brown dwarf and exoplanet atmospheres.

Recent analysis of Gemini-North/NIFS H-band (1.45 - 1.8 $\mu$m) observations of Uranus, recorded in 2010, with recently updated line data has revealed the spectral signature of hydrogen sulphide (H$_2$S) in Uranus's atmosphere (Irwin et al., 2018). Here, we extend this analysis to Gemini-North/NIFS observations of Neptune recorded in 2009 and find a similar detection of H$_2$S spectral absorption features in the 1.57 - 1.58 $\mu$m range, albeit slightly less evident, and retrieve a mole fraction of $\sim1-3$ ppm at the cloud tops. We find a much clearer detection (and much higher retrieved column abundance above the clouds) at southern polar latitudes compared with equatorial latitudes, which suggests a higher relative humidity of H$_2$S here. We find our retrieved H$_2$S abundances are most consistent with atmospheric models that have reduced methane abundance near Neptune's south pole, consistent with HST/STIS determinations (Karkoschka and Tomasko, 2011). We also conducted a Principal Component Analysis (PCA) of the Neptune and Uranus data and found that in the 1.57 - 1.60 $\mu$m range, some of the Empirical Orthogonal Functions (EOFs) mapped closely to physically significant quantities, with one being strongly correlated with the modelled H$_2$S signal and clearly mapping the spatial dependence of its spectral detectability. Just as for Uranus, the detection of H$_2$S at the cloud tops constrains the deep bulk sulphur/nitrogen abundance to exceed unity (i.e. $ > 4.4 - 5.0$ times the solar value) in Neptune's bulk atmosphere, provided that ammonia is not sequestered at great depths, and places a lower limit on its mole fraction below the observed cloud of (0.4 - 1.3) $\times 10^{-5}$. The detection of gaseous H$_2$S at these pressure levels adds to the weight of evidence that the principal constituent of the 2.5 - 3.5-bar cloud is likely to be H$_2$S ice.

An analysis of currently available ammonia (NH$_3$) visible-to-near-infrared gas absorption data was recently undertaken by Irwin et al. (Icarus, 302 (2018) 426) to help interpret Very Large Telescope (VLT) MUSE observations of Jupiter from 0.48 - 0.93 $\mu$m, made in support of the NASA/Juno mission. Since this analysis a newly revised set of ammonia line data, covering the previously poorly constrained range 0.5 - 0.833 $\mu$m, has been released by the ExoMol project, "C2018" (Coles et al., JQSRT 219, 199 - 122, 2018), which demonstrates significant advantages over previously available data sets, and providing for the first time complete line data for the previously poorly constrained 5520- and 6475-\AA\ bands of NH$_3$. In this paper we compare spectra calculated using the ExoMol-C2018 data set (Coles et al., JQSRT 219, 199 - 122, 2018) with spectra calculated from previous sources to demonstrate its advantages. We conclude that at the present time the ExoMol-C2018 dataset provides the most reliable ammonia absorption source for analysing low- to medium-resolution spectra of Jupiter in the visible/near-IR spectral range, but note that the data are less able to model high-resolution spectra owing to small, but significant inaccuracies in the line wavenumber estimates. This work is of significance not only for solar system planetary physics, but for future proposed observations of Jupiter-like planets orbiting other stars, such as with NASA's planned Wide-Field Infrared Survey Telescope (WFIRST).

The James Webb Space Telescope (JWST) may be capable of finding biogenic gases in the atmospheres of habitable exoplanets around low mass stars. Considerable attention has been given to the detectability of biogenic oxygen, which could be found using an ozone proxy, but ozone detection with JWST will be extremely challenging, even for the most favorable targets. Here, we investigate the detectability of biosignatures in anoxic atmospheres analogous to those that likely existed on the early Earth. Arguably, such anoxic biosignatures could be more prevalent than oxygen biosignatures if life exists elsewhere. Specifically, we simulate JWST retrievals of TRAPPIST-1e to determine whether the methane plus carbon dioxide disequilibrium biosignature pair is detectable in transit transmission. We find that ~10 transits using the Near InfraRed Spectrograph (NIRSpec) prism instrument may be sufficient to detect carbon dioxide and constrain methane abundances sufficiently well to rule out known, non-biological CH$_{4}$ production scenarios to ~90% confidence. Furthermore, it might be possible to put an upper limit on carbon monoxide abundances that would help rule out non-biological methane-production scenarios, assuming the surface biosphere would efficiently drawdown atmospheric CO. Our results are relatively insensitive to high altitude clouds and instrument noise floor assumptions, although stellar heterogeneity and variability may present challenges.

Previous analyses of the spectra of OB-type stars in the Magellanic Clouds have identified targets with low projected rotational velocities and relatively high nitrogen abundances. The VFTS obtained spectroscopy for B-type 434 stars. We have estimated atmospheric parameters and nitrogen abundances using TLUSTY model atmospheres for 54 B-type targets that appear to be single, have projected rotational velocities <80 \kms and were not classified as supergiants. In addition, nitrogen abundances for 34 similar stars from a previous survey have been re-evaluated. Approximately 75-80% of the targets have nitrogen enhancements of less than 0.3 dex, consistent with them having experienced only small amounts of mixing. However, stars with low projected rotational velocities <40 \kms and significant nitrogen enrichments are found in both our samples and simulations imply that these cannot all be rapidly rotating objects observed near pole-on. Adopting an enhancement threshold of 0.6 dex, we observed five stars, yet stellar evolution models with rotation predict only 1.25$\pm$1.11 for random stellar viewing inclinations. The excess of such objects is estimated to be 20-30% of all stars with current rotational velocities of less than 40 \kms and 2-4% of the total single B-type sample. These estimates constitute lower limits for stars that appear inconsistent with current grids of stellar evolutionary models. Including targets with smaller nitrogen enhancements of greater than 0.2 dex implies larger percentages of targets that are inconsistent with current evolutionary models, viz. 70% of the stars with rotational velocities less than 40 \kms and 6-8% of the total single stellar population. We consider possible explanations of which the most promising would appear to be breaking due to magnetic fields or stellar mergers with subsequent magnetic braking.

A spectroscopic analysis has been undertaken for the B-type multiple systems (excluding those with supergiant primaries) in the VLT-FLAMES Tarantula Survey (VFTS). Projected rotational velocities, $v$sin$i$, for the primaries have been estimated using a Fourier Transform technique and confirmed by fitting rotationally broadened profiles. A subset of 33 systems with $v$sin$i$$\leq$ 80 km/s have been analysed using a TLUSTY grid of model atmospheres to estimate stellar parameters and surface abundances for the primaries. The effects of a potential flux contribution from an unseen secondary have also been considered. For 20 targets it was possible to reliably estimate their effective temperatures (T$_{eff}$) but for the other 13 objects it was only possible to provide a constraint of 20,000$\leq$T$_{eff}$$\leq$26,000 K -- the other parameters estimated for these targets will be consequently less reliable. The estimated stellar properties are compared with evolutionary models and are generally consistent with their membership of 30 Doradus, while the nature of the secondaries of 3 SB2 system is discussed. A comparison with a sample of single stars with $v$sin$i$$\leq$ 80 km/s obtained from the VFTS and analysed with the same techniques implies that the atmospheric parameters and nitrogen abundances of the two samples are similar. However, the binary sample may have a lack of primaries with significant nitrogen enhancements, which would be consistent with them having low rotational velocities and having effectively evolved as single stars without significant rotational mixing. This result, which may be actually a consequence of the limitations of the pathfinder investigation presented in this paper, should be considered as a motivation for spectroscopic abundance analysis of large samples of binary stars, with high quality observational data.