Planet formation processes or evolution mechanisms are surmised to be at the origin of the hot Neptune desert. Studying exoplanets currently living within or at the edge of this desert could allow disentangling the respective roles of formation and evolution. We present the HARPS transmission spectrum of the bloated super-Neptune WASP-166b, located at the outer rim of the Neptune desert. Neutral sodium is detected at the 3.4 $\sigma$ level ($0.455 \pm 0.135 \%$), with a tentative indication of line broadening, which could be caused by winds blowing sodium farther into space, a possible manifestation of the bloated character of these highly irradiated worlds. We put this detection into context with previous work claiming a non-detection of sodium in the same observations and show that the high noise in the trace of the discarded stellar sodium lines was responsible for the non-detection. We highlight the impact of this low signal-to-noise remnant on detections for exoplanets similar to WASP-166b.

We carry out a phase curve analysis of the KELT-9 system using photometric observations from NASA's Transiting Exoplanet Survey Satellite. The measured secondary eclipse depth and peak-to-peak atmospheric brightness modulation are $650^{+14}_{-15}$ ppm and $566\pm16$ ppm, respectively. The planet's brightness variation reaches maximum $31\pm5$ minutes before the midpoint of secondary eclipse, indicating a $5\overset{\circ}{.}2\pm0\overset{\circ}{.}9$ eastward shift in the dayside hotspot from the substellar point. We also detect stellar pulsations on KELT-9 with a period of $7.58695\pm0.00091$ hours. The dayside emission of KELT-9b in the TESS bandpass is consistent with a blackbody brightness temperature of $4600\pm100$ K. The corresponding nightside brightness temperature is $3040\pm100$ K, comparable to the dayside temperatures of the hottest known exoplanets. In addition, we detect a significant phase curve signal at the first harmonic of the orbital frequency and a marginal signal at the second harmonic. While the amplitude of the first harmonic component is consistent with the predicted ellipsoidal distortion modulation assuming equilibrium tides, the phase of this photometric variation is shifted relative to the expectation. Placing KELT-9b in the context of other exoplanets with phase curve observations, we find that the elevated nightside temperature and relatively low day-night temperature contrast agree with the predictions of atmospheric models that include H$_{2}$ dissociation and recombination. The nightside temperature of KELT-9b implies an atmospheric composition containing about 50% molecular and 50% atomic hydrogen at 0.1 bar, a nightside emission spectrum that deviates significantly from a blackbody, and a 0.5-2.0 $\mu$m transmission spectrum that is featureless at low resolution.

Aims: We survey the transmission spectrum of WASP-121 b for line-absorption by metals and molecules at high spectral resolution, and elaborate on existing interpretations of the optical transmission spectrum observed with HST/STIS and WFC3. Methods: We use the cross-correlation technique and direct differential spectroscopy to search for sodium and other neutral and ionised atoms, TiO, VO and SH in high-resolution transit spectra obtained with the HARPS spectrograph. We inject models assuming chemical and hydrostatic equilibrium with varying temperature and composition to enable model comparison, and employ two bootstrap methods to test the robustness of our detections. Results: We detect neutral Mg, Na, Ca, Cr, Fe, Ni and V, which we predict exists in equilibrium with a significant quantity of VO, supporting earlier observations by HST/WFC3. Non-detections of Ti and TiO support the hypothesis that Ti is depleted via a cold-trap mechanism as has been proposed in the literature. Atomic line depths are under-predicted by hydrostatic models by a factor of 1.5 to 8, confirming recent findings that the atmosphere is extended. We predict the existence of significant concentrations of gas-phase TiO$_2$, VO$_2$ and TiS, which could be important absorbers at optical and NIR wavelengths in hot Jupiter atmospheres, but for which accurate line-list data is currently not available. We find no evidence for absorption by SH, and find that inflated atomic lines can plausibly explain the slope of the transmission spectrum observed in the NUV with HST/STIS. The Na D lines are significantly broadened and show a difference in their respective depths of 15 scale heights, which is not expected from isothermal hydrostatic theory.

The lifetime of a planetary disc which orbits a white dwarf represents a crucial input parameter into evolutionary models of that system. Here we apply a purely analytical formalism to estimate lifetimes of the debris phase of these discs, before they are ground down into dust or are subject to sublimation from the white dwarf. We compute maximum lifetimes for three different types of white dwarf discs, formed from (i) radiative YORP breakup of exo-asteroids along the giant branch phases at 2-100 au, (ii) radiation-less spin-up disruption of these minor planets at $\sim 1.5-4.5R_{\odot}$, and (iii) tidal disruption of minor or major planets within about $1.3R_{\odot}$. We display these maximum lifetimes as a function of disc mass and extent, constituent planetesimal properties, and representative orbital excitations of eccentricity and inclination. We find that YORP discs with masses up to $10^{24}$ kg live long enough to provide a reservoir of surviving cm-sized pebbles and m- to km-sized boulders that can be perturbed intact to white dwarfs with cooling ages of up to 10 Gyr. Debris discs formed from the spin or tidal disruption of these minor planets or major planets can survive in a steady state for up to respectively 1 Myr or 0.01 Myr, although most tidal discs would leave a steady state within about 1 yr. Our results illustrate that dust-less planetesimal transit detections are plausible, and would provide particularly robust evolutionary constraints. Our formalism can easily be adapted to individual systems and future discoveries.

THOR is the first open-source general circulation model (GCM) developed from scratch to study the atmospheres and climates of exoplanets, free from Earth- or Solar System-centric tunings. It solves the general non-hydrostatic Euler equations (instead of the primitive equations) on a sphere using the icosahedral grid. In the current study, we report major upgrades to THOR, building upon the work of Mendon\c{c}a et al. (2016). First, while the Horizontally Explicit Vertically Implicit (HEVI) integration scheme is the same as that described in Mendon\c{c}a et al. (2016), we provide a clearer description of the scheme and improved its implementation in the code. The differences in implementation between the hydrostatic shallow (HSS), quasi-hydrostatic deep (QHD) and non-hydrostatic deep (NHD) treatments are fully detailed. Second, standard physics modules are added: two-stream, double-gray radiative transfer and dry convective adjustment. Third, THOR is tested on additional benchmarks: tidally-locked Earth, deep hot Jupiter, acoustic wave, and gravity wave. Fourth, we report that differences between the hydrostatic and non-hydrostatic simulations are negligible in the Earth case, but pronounced in the hot Jupiter case. Finally, the effects of the so-called "sponge layer", a form of drag implemented in most GCMs to provide numerical stability, are examined. Overall, these upgrades have improved the flexibility, user-friendliness, and stability of THOR.

The determination of the fundamental properties (mass, separation, age, gravity and atmospheric properties) of brown dwarf companions allows us to infer crucial informations on their formation and evolution mechanisms. Spectroscopy of substellar companions is available to date only for a limited number of objects (and mostly at very low resolution, R<50) because of technical limitations, i.e., contrast and angular resolution. We present medium resolution (R=350), coronagraphic long-slit spectroscopic observations with SPHERE of two substellar companions, HD 1160 B and HD 19467 B. We found that HD 1160 B has a peculiar spectrum that cannot be fitted by spectra in current spectral libraries. A good fit is possible only considering separately the Y+J and the H spectral band. The spectral type is between M5 and M7. We also estimated a T_eff of 2800-2900 K and a log(g) of 3.5-4.0 dex. The low surface gravity seems to favour young age (10-20 Myr) and low mass (~20 M Jup ) for this object. HD 19467 B is instead a fully evolved object with a T_eff of ~1000 K and log g of ~5.0 dex. Its spectral type is T6+/-1.

Warm Neptunes offer a rich opportunity for understanding exo-atmospheric chemistry. With the upcoming James Webb Space Telescope (JWST), there is a need to elucidate the balance between investments in telescope time versus scientific yield. We use the supervised machine learning method of the random forest to perform an information content analysis on a 11-parameter model of transmission spectra from the various NIRSpec modes. The three bluest medium-resolution NIRSpec modes (0.7 - 1.27 microns, 0.97 - 1.84 microns, 1.66 - 3.07 microns) are insensitive to the presence of CO. The reddest medium-resolution mode (2.87 - 5.10 microns) is sensitive to all of the molecules assumed in our model: CO, CO2, CH4, C2H2, H2O, HCN and NH3. It competes effectively with the three bluest modes on the information encoded on cloud abundance and particle size. It is also competitive with the low-resolution prism mode (0.6 - 5.3 microns) on the inference of every parameter except for the temperature and ammonia abundance. We recommend astronomers to use the reddest medium-resolution NIRSpec mode for studying the atmospheric chemistry of 800-1200 K warm Neptunes; its corresponding high-resolution counterpart offers diminishing returns. We compare our findings to previous JWST information content analyses that favor the blue orders, and suggest that the reliance on chemical equilibrium could lead to biased outcomes if this assumption does not apply. A simple, pressure-independent diagnostic for identifying chemical disequilibrium is proposed based on measuring the abundances of H2O, CO and CO2.

We report detections of atomic species in the atmosphere of MASCARA-2 b, using the first transit observations obtained with the newly commissioned EXPRES spectrograph. EXPRES is a highly stabilised optical echelle spectrograph, designed to detect stellar reflex motions with amplitudes down to 30 cm/s, and was recently deployed at the Lowell Discovery Telescope. By analysing the transmission spectrum of the ultra-hot Jupiter MASCARA-2 b using the cross-correlation method, we confirm previous detections of Fe I, Fe II and Na I, which likely originate in the upper regions of the inflated atmosphere. In addition, we report significant detections of Mg I and Cr II. The absorption strengths change slightly with time, possibly indicating different temperatures and chemistry in the day-side and night-side terminators. Using the effective stellar line-shape variation induced by the transiting planet, we constrain the projected spin-orbit misalignment of the system to $1.6\pm3.1$ degrees, consistent with an aligned orbit. We demonstrate that EXPRES joins a suite of instruments capable of phase-resolved spectroscopy of exoplanet atmospheres.

Over large timescales, a terrestrial planet may be driven towards spin-orbit synchronous rotation by tidal forces. In this particular configuration, the planet exhibits permanent dayside and nightside, which may induce strong day-night temperature gradients. The nightside temperature depends on the efficiency of the day-night heat redistribution and determines the stability of the atmosphere against collapse. To better constrain the atmospheric stability, climate, and surface conditions of rocky planets located in the habitable zone of their host star, it is thus crucial to understand the complex mechanism of heat redistribution. Building on early works and assuming dry thermodynamics, we developed a hierarchy of analytic models taking into account the coupling between radiative transfer, dayside convection, and large-scale atmospheric circulation in the case of slowly rotating planets. There are two types of these models: a zero-dimensional two-layer approach and a two-column radiative-convective-subsiding-upwelling (RCSU) model. They yield analytical solutions and scaling laws characterising the dependence of the collapse pressure on physical features, which are compared to the results obtained by early works using 3D global climate models (GCMs). The analytical theory captures (i) the dependence of temperatures on atmospheric opacities and scattering in the shortwave and in the longwave, (ii) the behaviour of the collapse pressure observed in GCM simulations at low stellar fluxes that are due to the non-linear dependence of the atmospheric opacity on the longwave optical depth at the planet's surface, (iii) the increase of stability generated by dayside sensible heating, and (iv) the decrease of stability induced by the increase of the planet size.

Spectral retrieval has long been a powerful tool for interpreting planetary remote sensing observations. Flexible, parameterised, agnostic models are coupled with inversion algorithms in order to infer atmospheric properties directly from observations, with minimal reliance on physical assumptions. This approach, originally developed for application to Earth satellite data and subsequently observations of other Solar System planets, has been recently and successfully applied to transit, eclipse and phase curve spectra of transiting exoplanets. In this review, we present the current state-of-the-art in terms of our ability to accurately retrieve information about atmospheric chemistry, temperature, clouds and spatial variability; we discuss the limitations of this, both in the available data and modelling strategies used; and we recommend approaches for future improvement.

We present a new method for performing atmospheric retrieval on ground-based, high-resolution data of exoplanets. Our method combines cross-correlation functions with a random forest, a supervised machine learning technique, to overcome challenges associated with high-resolution data. A series of cross-correlation functions are concatenated to give a "CCF-sequence" for each model atmosphere, which reduces the dimensionality by a factor of ~100. The random forest, trained on our grid of ~65,000 models, provides a likelihood-free method of retrieval. The pre-computed grid spans 31 values of both temperature and metallicity, and incorporates a realistic noise model. We apply our method to HARPS-N observations of the ultra-hot Jupiter KELT-9b, and obtain a metallicity consistent with solar (logM = $-0.2\pm0.2$). Our retrieved transit chord temperature (T = $6000^{+0}_{-200}$K) is unreliable as the ion cross-correlations lie outside of the training set, which we interpret as being indicative of missing physics in our atmospheric model. We compare our method to traditional nested-sampling, as well as other machine learning techniques, such as Bayesian neural networks. We demonstrate that the likelihood-free aspect of the random forest makes it more robust than nested-sampling to different error distributions, and that the Bayesian neural network we tested is unable to reproduce complex posteriors. We also address the claim in Cobb et al. (2019) that our random forest retrieval technique can be over-confident but incorrect. We show that this is an artefact of the training set, rather than the machine learning method, and that the posteriors agree with those obtained using nested-sampling.

We present an improved, hybrid CPU-GPU atmospheric retrieval code, Helios-r2, which is applicable to medium-resolution emission spectra of brown dwarfs, in preparation for precision atmospheric spectroscopy in the era of the James Webb Space Telescope. The model is available as open-source code on the Exoclimes Simulation Platform. We subject Helios-r2 to a battery of tests of varying difficulty. The simplest test involves a mock retrieval on a forward model generated using the same radiative transfer technique, the same implementation of opacities, and the same chemistry model. The least trivial test involves a mock retrieval on synthetic spectra from the Sonora model grid, which uses a different radiative transfer technique, a different implementation of opacities, and a different chemistry model. A calibration factor, which is included to capture uncertainties in the brown dwarf radius, distance to the brown dwarf and flux calibration of the spectrum, may compensate, sometimes erroneously, for discrepancies in modeling choices and implementation. We analyze spectra of the benchmark brown dwarf GJ 570 D and the binary brown dwarf companions in the Epsilon Indi system. The retrieved surface gravities are consistent with previous studies and/or values inferred from dynamical masses (for Epsilon Indi Ba and Bb only). There remains no clear criterion on how to reject unphysical values of the retrieved brown dwarf radii. The inferred radii and corresponding masses should be taken with great caution. The retrieved carbon-to-oxygen ratios and metallicity depend on whether chemical equilibrium is assumed.

Ultra-hot Jupiters offer interesting prospects for expanding our theories on dynamical evolution and the properties of extremely irradiated atmospheres. In this context, we present the analysis of new optical spectroscopy for the transiting ultra-hot Jupiter WASP-121b. We first refine the orbital properties of WASP-121b, which is on a nearly polar (obliquity $\psi^{\rm North}$=88.1$\pm$0.25$^{\circ}$ or $\psi^{\rm South}$=91.11$\pm$0.20$^{\circ}$) orbit, and exclude a high differential rotation for its fast-rotating (P$<$1.13 days), highly inclined ($i_\mathrm{\star}^{\rm North}$=8.1$\stackrel{+3.0}{_{-2.6}}^{\circ}$ or $i_\mathrm{\star}^{\rm South}$=171.9$\stackrel{+2.5}{_{-3.4}}^{\circ}$) star. We then present a new method that exploits the reloaded Rossiter-McLaughlin technique to separate the contribution of the planetary atmosphere and of the spectrum of the stellar surface along the transit chord. Its application to HARPS transit spectroscopy of WASP-121b reveals the absorption signature from metals, likely atomic iron, in the planet atmospheric limb. The width of the signal (14.3$\pm$1.2 km/s) can be explained by the rotation of the tidally locked planet. Its blueshift (-5.2$\pm$0.5 km/s) could trace strong winds from the dayside to the nightside, or the anisotropic expansion of the planetary thermosphere.

Self-consistent model grids of brown dwarfs involve complex physics and chemistry, and are often computed using proprietary computer codes, making it challenging to identify the reasons for discrepancies between model and data as well as between the models produced by different research groups. In the current study, we demonstrate a novel method for analyzing brown dwarf spectra, which combines the use of the Sonora, AMES-Cond and HELIOS model grids with the supervised machine learning method of the random forest. Besides performing atmospheric retrieval, the random forest enables information content analysis of the three model grids as a natural outcome of the method, both individually on each grid and by comparing the grids against one another, via computing large suites of mock retrievals. Our analysis reveals that the different choices made in modelling the alkali line shapes hinder the use of the alkali lines as gravity indicators. Nevertheless, the spectrum longward of 1.2 micron encodes enough information on the surface gravity to allow its inference from retrieval. Temperature may be accurately and precisely inferred independent of the choice of model grid, but not the surface gravity. We apply random forest retrieval to three objects: the benchmark T7.5 brown dwarf GJ 570D; and Epsilon Indi Ba (T1.5 brown dwarf) and Bb (T6 brown dwarf), which are part of a binary system and have measured dynamical masses. For GJ 570D, the inferred effective temperature and surface gravity are consistent with previous studies. For Epsilon Indi Ba and Bb, the inferred surface gravities are broadly consistent with the values informed by the dynamical masses.

Two theoretical quandaries involving transmission spectra of gas-giant exoplanets are elucidated. When computing the transit radius as a function of wavelength, one needs to specify a reference transit radius corresponding to a reference pressure. Mathematically, the reference transit radius is a constant of integration that originates from evaluating an integral for the transit depth. Physically, its interpretation has been debated in the literature. Jordan & Espinoza (2018) suggested that the optical depth is discontinuous across, and infinite below, the reference transit radius. Betremieux & Swain (2017, 2018) interpreted the spherical surface located at the reference transit radius to represent the boundary associated with an opaque cloud deck. It is demonstrated that continuous functions for the optical depth may be found. The optical depth below and at the reference transit radius need not take on special or divergent values. In the limit of a spatially uniform grey cloud with constant opacity, the transit chord with optical depth on the order of unity mimics the presence of a "cloud top". While the surface located at the reference pressure may mimic the presence of grey clouds, it is more natural to include the effects of these clouds as part of the opacity function because the cloud opacity may be computed from first principles. It is unclear how this mimicry extends to non-grey clouds comprising small particles.

We present the analysis of TESS optical photometry of WASP-121b, which reveal the phase curve of this transiting ultra-hot Jupiter. Its hotspot is located at the substellar point, showing inefficient heat transport from the dayside (2870 K) to the nightside ($<$ 2200 K) at the altitudes probed by TESS. The TESS eclipse depth, measured at the shortest wavelength to date for WASP-121b, confirms the strong deviation from blackbody planetary emission. Our atmospheric retrieval on the complete emission spectrum supports the presence of a temperature inversion, which can be explained by the presence of VO and possibly TiO and FeH. The strong planetary emission at short wavelengths could arise from an H$^{-}$ continuum.

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.

Motivated by both ground- and space-based detections of the sodium doublet in the transmission spectra of exoplanetary atmospheres, we revisit the theory and interpretation of sodium lines in non-local thermodynamic equilibrium (NLTE), where collisions are not efficient enough to maintain a Boltzmann distribution for the excited and ground states of the sodium atom. We consider non-Boltzmann distributions that account for the ineffectiveness of collisions. We analyze the sodium doublet in transmission spectra measured at low (HAT-P-1b, HAT-P-12b, HD 189733b, WASP-6b, WASP-17b and WASP-39b) and high (WASP-49b) spectral resolutions. Nested-sampling retrievals performed on low-resolution optical/visible transmission spectra are unable to break the normalization degeneracy if the spectral continuum is associated with Rayleigh scattering by small cloud particles. Using mock retrievals, we demonstrate that un-normalized ground-based, high-resolution spectra centered on the sodium doublet alone are unable to precisely inform us about the pressure levels probed by the transit chord and hence to identify the region (i.e., thermosphere, exosphere) of the atmosphere being probed. Retrievals performed on the HARPS transmission spectrum of WASP-49b support this conclusion. Generally, we are unable to distinguish between LTE versus NLTE interpretations of the sodium doublet based on the computed Bayesian evidence with the implication that LTE interpretations tend to under-estimate the temperature probed by the transit chord. With the current low-resolution data, the sodium line shapes are consistent with Voigt profiles without the need for sub-Lorentzian wings. The retrieved sodium abundances are consistent with being sub-solar to solar.

Context: KELT-9 b exemplifies a newly emerging class of short-period gaseous exoplanets that tend to orbit hot, early type stars - termed ultra-hot Jupiters. The severe stellar irradiation heats their atmospheres to temperatures of $\sim 4,000$ K, similar to the photospheres of dwarf stars. Due to the absence of aerosols and complex molecular chemistry at such temperatures, these planets offer the potential of detailed chemical characterisation through transit and day-side spectroscopy. Studies of their chemical inventories may provide crucial constraints on their formation process and evolution history. Aims: To search the optical transmission spectrum of KELT-9 b for absorption lines by metals using the cross-correlation technique. Methods: We analyse 2 transits observed with the HARPS-N spectrograph. We use an isothermal equilibrium chemistry model to predict the transmission spectrum for each of the neutral and singly-ionized atoms with atomic numbers between 3 and 78. Of these, we identify the elements that are expected to have spectral lines in the visible wavelength range and use those as cross-correlation templates. Results: We detect absorption of Na I, Cr II, Sc II and Y II, and confirm previous detections of Mg I, Fe I, Fe II and Ti II. In addition, we find evidence of Ca I, Cr I, Co I, and Sr II that will require further observations to verify. The detected absorption lines are significantly deeper than model predictions, suggesting that material is transported to higher altitudes where the density is enhanced compared to a hydrostatic profile. There appears to be no significant blue-shift of the absorption spectrum due to a net day-to-night side wind. In particular, the strong Fe II feature is shifted by $0.18 \pm 0.27$ km~s$^{-1}$, consistent with zero. Using the orbital velocity of the planet we revise the steller and planetary masses and radii.

We present new methodological features and physical ingredients included in the 1D radiative transfer code HELIOS, improving the hemispheric two-stream formalism. We conduct a thorough intercomparison survey with several established forward models, including COOLTLUSTY, PHOENIX, and find satisfactory consistency with their results. Then, we explore the impact of (i) different groups of opacity sources, (ii) a stellar path length adjustment, and (iii) a scattering correction on self-consistently calculated atmospheric temperatures and planetary emission spectra. First, we observe that temperature-pressure (T-P) profiles are very sensitive to the opacities included, with metal oxides, hydrides, the alkali atoms (and ionized hydrogen) playing an important role for the absorption of shortwave radiation (in very hot surroundings). Moreover, if these species are sufficiently abundant, they are likely to induce non-monotonic T-P profiles. Second, without the stellar path length adjustment, the incoming stellar flux is significantly underestimated for zenith angles above 80{\deg}, which somewhat affects the upper atmospheric temperatures and the planetary emission. Third, the scattering correction improves the accuracy of the computation of the reflected stellar light by ~10%. We use HELIOS to calculate a grid of cloud-free atmospheres in radiative-convective equilibrium for self-luminous planets for a range of effective temperatures, surface gravities, metallicities, and C/O ratios, to be used by planetary evolution studies. Furthermore, we calculate dayside temperatures and secondary eclipse spectra for a sample of exoplanets for varying chemistry and heat redistribution. These results may be used to make predictions on the feasibility of atmospheric characterizations with future observations.

High-resolution optical spectroscopy is a powerful tool to characterise exoplanetary atmospheres from the ground. The sodium D lines, with their large cross sections, are especially suited to study the upper layers of atmospheres in this context. We report on the results from HEARTS, a spectroscopic survey of exoplanet atmospheres, performing a comparative study of hot gas giants to determine the effects of stellar irradiation. In this second installation of the series, we highlight the detection of neutral sodium on the ultra-hot giant WASP-76b. We observed three transits of the planet using the HARPS high-resolution spectrograph at the ESO 3.6m telescope and collected 175 spectra of WASP-76. We repeatedly detect the absorption signature of neutral sodium in the planet atmosphere ($0.371\pm0.034\%$; $10.75 \sigma$ in a $0.75$ \r{A} passband). The sodium lines have a Gaussian profile with full width at half maximum (FWHM) of $27.6\pm2.8$ km s$^{-1}$. This is significantly broader than the line spread function of HARPS ($2.7$ km s$^{-1}$). We surmise that the observed broadening could trace the super-rotation in the upper atmosphere of this ultra-hot gas giant.

Spectral features in the observed spectra of exoplanets depend on the composition of their atmospheres. A good knowledge of the main atmospheric processes that drive the chemical distribution is therefore essential to interpret exoplanetary spectra. An atmosphere reaches chemical equilibrium if the rates of the forward and backward chemical reactions converge to the same value. However, there are atmospheric processes, such as atmospheric transport, that destabilize this equilibrium. In this work we study the changes in composition driven by a 3D wind field in WASP-43b using our Global Circulation Model, THOR. Our model uses validated temperature- and pressure-dependent chemical timescales that allow us to explore the disequilibrium chemistry of CO, CO$_2$, H$_2$O and CH$_4$. In WASP-43b the formation of the equatorial jet has an important impact in the chemical distribution of the different species across the atmosphere. At low latitudes the chemistry is longitudinally quenched, except for CO$_2$ at solar abundances. The polar vortexes have a distinct chemical distribution since these are regions with lower temperature and atmospheric mixing. Vertical and latitudinal mixing have a secondary impact in the chemical transport. We determine graphically the effect of disequilibrium on observed emission spectra. Our results do not show any significant differences in the emission spectra between the equilibrium and disequilibrium solutions for C/O = 0.5. However, if C/O is increased to 2.0, differences in the spectra due to the disequilibrium chemistry of CH$_4$ become non-negligible. In some spectral ranges the emission spectra can have more than 15$\%$ departures from the equilibrium solution.

A comprehensive analysis of 38 previously published Wide Field Camera 3 (WFC3) transmission spectra is performed using a hierarchy of nested-sampling retrievals: with versus without clouds, grey versus non-grey clouds, isothermal versus non-isothermal transit chords and with water, hydrogen cyanide and/or ammonia. We revisit the "normalisation degeneracy": the relative abundances of molecules are degenerate at the order-of-magnitude level with the absolute normalisation of the transmission spectrum. Using a suite of mock retrievals, we demonstrate that the normalisation degeneracy may be partially broken using WFC3 data alone, even in the absence of optical/visible data and without appealing to the presence of patchy clouds, although lower limits to the mixing ratios may be prior-dominated depending on the measurement uncertainties. With James Webb Space Telescope-like spectral resolutions, the normalisation degeneracy may be completely broken from infrared spectra alone. We find no trend in the retrieved water abundances across nearly two orders of magnitude in exoplanet mass and a factor of 5 in retrieved temperature (about 500 to 2500 K). We further show that there is a general lack of strong Bayesian evidence to support interpretations of non-grey over grey clouds (only for WASP-69b and WASP-76b) and non-isothermal over isothermal atmospheres (no objects). 35 out of 38 WFC3 transmission spectra are well-fitted by an isothermal transit chord with grey clouds and water only, while 8 are adequately explained by flat lines. Generally, the cloud composition is unconstrained.

A key legacy of the recently launched TESS mission will be to provide the astronomical community with many of the best transiting exoplanet targets for atmospheric characterization. However, time is of the essence to take full advantage of this opportunity. JWST, although delayed, will still complete its nominal five year mission on a timeline that motivates rapid identification, confirmation, and mass measurement of the top atmospheric characterization targets from TESS. Beyond JWST, future dedicated missions for atmospheric studies such as ARIEL require the discovery and confirmation of several hundred additional sub-Jovian size planets (R_p < 10 R_Earth) orbiting bright stars, beyond those known today, to ensure a successful statistical census of exoplanet atmospheres. Ground-based ELTs will also contribute to surveying the atmospheres of the transiting planets discovered by TESS. Here we present a set of two straightforward analytic metrics, quantifying the expected signal-to-noise in transmission and thermal emission spectroscopy for a given planet, that will allow the top atmospheric characterization targets to be readily identified among the TESS planet candidates. Targets that meet our proposed threshold values for these metrics would be encouraged for rapid follow-up and confirmation via radial velocity mass measurements. Based on the catalog of simulated TESS detections by Sullivan et al. (2015), we determine appropriate cutoff values of the metrics, such that the TESS mission will ultimately yield a sample of $\sim300$ high-quality atmospheric characterization targets across a range of planet size bins, extending down to Earth-size, potentially habitable worlds.

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 observations of WASP-63b by the Hubble Space Telescope (HST) as part of "A Preparatory Program to Identify the Single Best Transiting Exoplanet for JWST Early Release Science". WASP-63b is one of the community targets under consideration for the James Webb Space Telescope (JWST) Early Release Science (ERS) program. We present a spectrum derived from a single observation by HST Wide Field Camera 3 in the near infrared. We engaged groups across the transiting exoplanet community to participate in the analysis of the data and present results from each. There is general agreement amongst all results that we find an H2O absorption feature with 3.5-4.0 sigma significance. However, the feature is muted in comparison to a clear atmosphere at solar composition. Although the detection of the water feature is robust, the reasons for the muting of this feature are ambiguous due to a degeneracy between clouds and composition. The data does not yield robust detections of any molecular species other than H2O. The group was motivated to perform an additional set of retrieval exercises to investigate an apparent bump in the spectrum at ~ 1.55 um. We explore possible disequilibrium chemistry and find this feature is consistent with super-solar HCN abundance but it is questionable if the required mixing ratio of HCN is chemically and physically plausible. The ultimate goal of this study is to vet WASP-63b as a potential community target to best demonstrate the capabilities and systematics of JWST instruments for transiting exoplanet science. In the case of WASP-63b, the presence of a detectable water feature indicates that WASP-63b remains a plausible target for ERS observations.

The chemical composition of an exoplanet is a key ingredient in constraining its formation history. Iron is the most abundant transition metal, but has never been directly detected in an exoplanet due to its highly refractory nature. KELT-9b (HD 195689b) is the archetype of the class of ultra-hot Jupiters that straddle the transition between stars and gas-giant exoplanets and serve as distinctive laboratories for studying atmospheric chemistry, because of its high equilibrium temperature of 4050 +/- 180 K. These properties imply that its atmosphere is a tightly constrained chemical system that is expected to be nearly in chemical equilibrium and cloud-free. It was previously predicted that the spectral lines of iron will be detectable in the visible range of wavelengths. At these high temperatures, iron and several other transition metals are not sequestered in molecules or cloud particles and exist solely in their atomic forms. Here, we report the direct detection of atomic neutral and singly-ionized iron (Fe and Fe+), and singly-ionized titanium (Ti+) in KELT-9b via the cross-correlation technique applied to high-resolution spectra obtained during the primary transit of the exoplanet.

One of the desired outcomes of studying exoplanetary atmospheres is to set decisive constraints on exoplanet formation theories. Astronomers often speak of the "metallicity" in broad terms. When modeling the bulk metallicity, workers refer to the elemental iron abundance (Fe/H). By contrast, when exo-atmospheric astronomers speak of the "metallicity" derived from analysing low-resolution Hubble and Spitzer spectrophotometry, they are referring to the elemental abundances of oxygen (O/H), carbon (C/H) and nitrogen (N/H): in decreasing order of importance, since spectra from the Hubble Wide Field Camera 3 are primarily sensitive to water with secondary contributions from hydrogen cyanide and ammonia, while Spitzer photometry is sensitive to methane and carbon monoxide. From retrieving for the water abundances, workers such as Kreidberg et al. (2014), Wakeford et al. (2017, 2018), Arcangeli et al. (2018) and Mansfield et al. (2018) have published figures that show the "metallicity" versus the masses of the exoplanets. Even if we assume that bulk and atmospheric elemental abundances are equal, the motivation behind this research note is to demonstrate that the conversion between the water abundance and O/H is not straightforward. "Water at solar metallicity" is a temperature- and pressure-dependent statement, because there is no single value of the water volume mixing ratio corresponding to solar values of O/H and C/H. Overall, the conversion factor between the water volume mixing ratio and O/H is not unity and depends on temperature, pressure, O/H, C/H and C/O, even without invoking disequilibrium chemistry (photochemistry and/or atmospheric mixing). It remains unproven that disequilibrium chemistry or condensation will always return a conversion factor of unity. It is thus a fair statement to say that the "metallicity" estimates inferred from atmospheric retrievals are model-dependent.

The atmospheric temperatures of the ultra-hot Jupiter KELT-9b straddle the transition between gas giants and stars, and therefore between two traditionally distinct regimes of atmospheric chemistry. Previous theoretical studies assume the atmosphere of KELT-9b to be in chemical equilibrium. Despite the high ultraviolet flux from KELT-9, we show using photochemical kinetics calculations that the observable atmosphere of KELT-9b is predicted to be close to chemical equilibrium, which greatly simplifies any theoretical interpretation of its spectra. It also makes the atmosphere of KELT-9b, which is expected to be cloudfree, a tightly constrained chemical system that lends itself to a clean set of theoretical predictions. Due to the lower pressures probed in transmission (compared to emission) spectroscopy, we predict the abundance of water to vary by several orders of magnitude across the atmospheric limb depending on temperature, which makes water a sensitive thermometer. Carbon monoxide is predicted to be the dominant molecule under a wide range of scenarios, rendering it a robust diagnostic of the metallicity when analyzed in tandem with water. All of the other usual suspects (acetylene, ammonia, carbon dioxide, hydrogen cyanide, methane) are predicted to be subdominant at solar metallicity, while atomic oxygen, iron and magnesium are predicted to have relative abundances as high as 1 part in 10,000. Neutral atomic iron is predicted to be seen through a forest of optical and near-infrared lines, which makes KELT-9b suitable for high-resolution ground-based spectroscopy with HARPS-N or CARMENES. We summarize future observational prospects of characterizing the atmosphere of KELT-9b.

Two-stream radiative transfer is a workhorse in the Earth, planetary and exoplanetary sciences due to its simplicity and ease of implementation. However, a longstanding limitation of the two-stream approximation is its inaccuracy in the presence of medium-sized or large aerosols. This limitation was lifted by the discovery of the improved two-stream technique, where the accuracy of the scattering greenhouse effect is matched to that of multi-stream calculations by construction. In the present study, we derive the full solutions for improved two-stream radiative transfer, following its introduction by Heng & Kitzmann (2017), and include contributions from the direct stellar beam. The generalization of the original two-stream flux solutions comes in the form of a correction factor, traditionally set to unity, which is the ratio of a pair of first Eddington coefficients. We derive an analytical expression for this correction factor and also provide a simple fitting function for its ease of use by other workers. We prove that the direct stellar beam is associated with a second Eddington coefficient that is on the order of unity. Setting this second Eddington coefficient to 2/3 and $1/\sqrt{3}$ reproduces the Eddington and quadrature closures, respectively, associated with the direct beam. We use our improved two-stream solutions for the fluxes to derive two-stream source function solutions for both the intensity and fluxes.

The use of machine learning is becoming ubiquitous in astronomy, but remains rare in the study of the atmospheres of exoplanets. Given the spectrum of an exoplanetary atmosphere, a multi-parameter space is swept through in real time to find the best-fit model. Known as atmospheric retrieval, it is a technique that originates from the Earth and planetary sciences. Such methods are very time-consuming and by necessity there is a compromise between physical and chemical realism versus computational feasibility. Machine learning has previously been used to determine which molecules to include in the model, but the retrieval itself was still performed using standard methods. Here, we report an adaptation of the random forest method of supervised machine learning, trained on a pre-computed grid of atmospheric models, which retrieves full posterior distributions of the abundances of molecules and the cloud opacity. The use of a pre-computed grid allows a large part of the computational burden to be shifted offline. We demonstrate our technique on a transmission spectrum of the hot gas-giant exoplanet WASP-12b using a five-parameter model (temperature, a constant cloud opacity and the volume mixing ratios or relative abundance by number of water, ammonia and hydrogen cyanide). We obtain results consistent with the standard nested-sampling retrieval method. Additionally, we can estimate the sensitivity of the measured spectrum to constraining the model parameters and we can quantify the information content of the spectrum. Our method can be straightforwardly applied using more sophisticated atmospheric models and also to interpreting an ensemble of spectra without having to retrain the random forest.

Motivated by the work of Cooper & Showman, we revisit the chemical relaxation method, which seeks to enhance the computational efficiency of chemical-kinetics calculations by replacing the chemical network with a handful of independent source/sink terms. Chemical relaxation solves the evolution of the system and can treat disequilibrium chemistry, as the source/sink terms are driven towards chemical equilibrium on a prescribed chemical timescale, but it has surprisingly never been validated. First, we generalize the treatment by forgoing the use of a single chemical timescale, instead developing a pathway analysis tool that allows us to identify the rate-limiting reaction as a function of temperature and pressure. For the interconversion between methane and carbon monoxide and between ammonia, and molecular nitrogen, we identify the key rate-limiting reactions for conditions relevant to currently characterizable exo-atmospheres (500-3000 K, 0.1 mbar to 1 kbar). Second, we extend chemical relaxation to include carbon dioxide and water. Third, we examine the role of metallicity and carbon-to-oxygen ratio in chemical relaxation. Fourth, we apply our pathway analysis tool to diagnose the differences between our chemical network and that of Moses and Venot. Finally, we validate the chemical relaxation method against full chemical kinetics calculations in one dimension. For WASP-18b-, HD 189733b- and GJ 1214-b-like atmospheres, we show that chemical relaxation is mostly accurate to within an order of magnitude, a factor of 2 and $\sim 10\%$, respectively. The level of accuracy attained allows for the chemical relaxation method to be included in three-dimensional general circulation models.

The study of high contrast imaged brown dwarfs and exoplanets depends strongly on evolutionary models. To estimate the mass of a directly imaged substellar object, its extracted photometry or spectrum is used and adjusted with model spectra together with the estimated age of the system. These models still need to be properly tested and constrained. HD 4747B is a brown dwarf close to the H burning mass limit, orbiting a nearby, solar-type star and has been observed with the radial velocity method over almost two decades now. Its companion was also recently detected by direct imaging, allowing a complete study of this particular object. We aim to fully characterize HD 4747B by combining a well constrained dynamical mass and a study of its observed spectral features in order to test evolutionary models for substellar objects and characterize its atmosphere. We combine the radial velocity measurements of HIRES and CORALIE taken over two decades and high contrast imaging of several epochs from NACO, NIRC2 and SPHERE to obtain a dynamical mass. From the SPHERE data we obtain a low resolution spectrum of the companion from Y to H band, as well as two narrow band-width photometric measurements in the K band. A study of the primary star allows in addition to constrain the age of the system as well as its distance. Thanks to the new SPHERE epoch and NACO archival data combined with previous imaging data and high precision radial velocity measurements, we have been able to derive a well constrained orbit. We derive a dynamical mass of mB=70.0$\pm$1.6 MJup which is higher than a previous study, but in better agreement with the models. By comparing the object with known brown dwarfs spectra, we derive a spectral type of L9 and an effective temperature of 1350$\pm$50 K. With a retrieval analysis we constrain the oxygen and carbon abundances and compare them with the ones from the HR 8799 planets.

Space-borne low-to medium-resolution (R~10^2-10^3) transmission spectroscopy of atmospheres detect the broadest spectral features (alkali doublets, molecular bands, scattering), while high-resolution (R~10^5), ground-based observations probe the sharpest features (cores of the alkali lines, molecular lines).The two techniques differ by:(1) The LSF of ground-based observations is 10^3 times narrower than for space-borne observations;(2)Space-borne transmission spectra probe up to the base of thermosphere, while ground-based observations can reach pressures down to 10^(-11);(3)Space-borne observations directly yield the transit depth of the planet, while ground-based observations measure differences in the radius of the planet at different wavelengths.It is challenging to combine both techniques.We develop a method to compare theoretical models with observations at different resolutions.We introduce PyETA, a line-by-line 1D radiative transfer code to compute transmission spectra at R~10^6 (0.01 A) over a broad wavelength range.An hybrid forward modeling/retrieval optimization scheme is devised to deal with the large computational resources required by modeling a broad wavelength range (0.3-2 $\mu$m) at high resolution.We apply our technique to HD189733b.Here, HST observations reveal a flattened spectrum due to scattering by aerosols, while high-resolution ground-based HARPS observations reveal the sharp cores of sodium lines.We reconcile these results by building models that reproduce simultaneously both data sets, from the troposphere to the thermosphere. We confirm:(1)the presence of scattering by tropospheric aerosols;(2)that the sodium core feature is of thermospheric origin.Accounting for aerosols, the sodium cores indicate T up to 10000K in the thermosphere.The precise value of the thermospheric temperature is degenerate with the abundance of sodium and altitude of the aerosol deck.

We apply the picket fence treatment to model the effects brought about by spectral lines on the thermal structure of irradiated atmospheres. The lines may be due to purely absorption processes, purely coherent scattering processes or some combination of absorption and scattering. If the lines arise as a pure absorption process, the surface layers of the atmosphere are cooler whereas this surface cooling is completely absent if the lines are due to pure coherent isotropic scattering. The lines also lead to a warming of the deeper atmosphere. The warming of the deeper layers is, however, independent of the nature of line formation. Accounting for coherent isotropic scattering in the shortwave and longwave continuum results in anti-greenhouse cooling and greenhouse warming on an atmosphere-wide scale. The effects of coherent isotropic scattering in the line and continuum operate in tandem to determine the resulting thermal structure of the irradiated atmosphere.

Recently acquired Hubble and Spitzer phase curves of the short-period hot Jupiter WASP-43b make it an ideal target for confronting theory with data. On the observational front, we re-analyze the 3.6 and 4.5 $\mu$m Spitzer phase curves and demonstrate that our improved analysis better removes residual red noise due to intra-pixel sensitivity, which leads to greater fluxes emanating from the nightside of WASP-43b, thus reducing the tension between theory and data. On the theoretical front, we construct cloudfree and cloudy atmospheres of WASP-43b using our Global Circulation Model (GCM), THOR, which solves the non-hydrostatic Euler equations (compared to GCMs that typically solve the hydrostatic primitive equations). The cloudfree atmosphere produces a reasonable fit to the dayside emission spectrum. The multi-phase emission spectra constrain the cloud deck to be confined to the nightside and have a finite cloud-top pressure. The multi-wavelength phase curves are naturally consistent with our cloudy atmospheres, except for the 4.5 $\mu$m phase curve, which requires the presence of enhanced carbon dioxide in the atmosphere of WASP-43b. Multi-phase emission spectra at higher spectral resolution, as may be obtained using the James Webb Space Telescope, and a reflected-light phase curve at visible wavelengths would further constrain the properties of clouds in WASP-43b.

Context. The TRAPPIST-1 system hosts seven Earth-sized, temperate exoplanets orbiting an ultra-cool dwarf star. As such, it represents a remarkable setting to study the formation and evolution of terrestrial planets that formed in the same protoplanetary disk. While the sizes of the TRAPPIST-1 planets are all known to better than 5% precision, their densities have significant uncertainties (between 28% and 95%) because of poor constraints on the planet's masses. Aims.The goal of this paper is to improve our knowledge of the TRAPPIST-1 planetary masses and densities using transit-timing variations (TTV). The complexity of the TTV inversion problem is known to be particularly acute in multi-planetary systems (convergence issues, degeneracies and size of the parameter space), especially for resonant chain systems such as TRAPPIST-1. Methods. To overcome these challenges, we have used a novel method that employs a genetic algorithm coupled to a full N-body integrator that we applied to a set of 284 individual transit timings. This approach enables us to efficiently explore the parameter space and to derive reliable masses and densities from TTVs for all seven planets. Results. Our new masses result in a five- to eight-fold improvement on the planetary density uncertainties, with precisions ranging from 5% to 12%. These updated values provide new insights into the bulk structure of the TRAPPIST-1 planets. We find that TRAPPIST-1\,c and e likely have largely rocky interiors, while planets b, d, f, g, and h require envelopes of volatiles in the form of thick atmospheres, oceans, or ice, in most cases with water mass fractions less than 5%.

We propose to search for biosignatures in the spectra of reflected light from about 100 Earth-sized planets that are already known to be orbiting in their habitable zones (HZ). For a sample of G and K type hosts, most of these planets will be between 25 and 50 milli-arcsec (mas) from their host star and 1 billion to 10 billion times fainter. To separate the planet's image from that of its host star at the wavelength (763nm) of the oxygen biosignature we need a telescope with an aperture of 16 metres. Furthermore, the intensity of the light from the host star at the position in the image of the exoplanet must be suppressed otherwise the exoplanet will be lost in the glare. This presents huge technical challenges. The Earth's atmosphere is turbulent which makes it impossible to achieve the required contrast from the ground at 763nm. The telescope therefore needs to be in space and to fit the telescope in the rocket fairing it must be a factor of 4 or more times smaller when folded than when operational. To obtain spectroscopy of the planet's biosignature at 763nm we need to use an integral field spectrometer (IFS) with a field of view (FOV) of 1000 x 1000 milli-arcsec (mas) and a spectral resolution of 100. This is a device that simultaneously takes many pictures of the exoplanet each at a slightly different wavelength which are then recorded as a data cube with two spatial dimensions and one wavelength dimension. In every data cube wavelength slice, the background light from the host star at the location of the planet image must be minimised. This is achieved via a coronagraph which blocks the light from the host star and active/adaptive optics techniques which continuously maintain very high accuracy optical alignment to make the images as sharp as possible. These are the technical challenges to be addressed in a design study.

The angle between the spin of a star and its planets' orbital planes traces the history of the planetary system. Exoplanets orbiting close to cool stars are expected to be on circular, aligned orbits because of strong tidal interactions with the stellar convective envelope. Spin-orbit alignment can be measured when the planet transits its star, but such ground-based spectroscopic measurements are challenging for cool, slowly-rotating stars. Here we report the characterization of a planet three-dimensional trajectory around an M dwarf star, derived by mapping the spectrum of the stellar photosphere along the chord transited by the planet. We find that the eccentric orbit of the Neptune-mass exoplanet GJ 436b is nearly perpendicular to the stellar equator. Both eccentricity and misalignment, surprising around a cool star, can result from dynamical interactions (via Kozai migration) with a yet-undetected outer companion. This inward migration of GJ 436b could have triggered the atmospheric escape that now sustains its giant exosphere. Eccentric, misaligned exoplanets orbiting close to cool stars might thus hint at the presence of unseen perturbers and illustrate the diversity of orbital architectures seen in exoplanetary systems.

The prevalence of clouds in currently observable exoplanetary atmospheres motivates the compilation and calculation of their optical properties. First, we present a new open-source Mie scattering code known as LX-MIE, which is able to consider large size parameters ($\sim 10^7$) using a single computational treatment. We validate LX-MIE against the classical MIEV0 code as well as previous studies. Second, we embark on an expanded survey of the published literature for both the real and imaginary components of the refractive indices of 32 condensate species. As much as possible, we rely on experimental measurements of the refractive indices and resort to obtaining the real from the imaginary component (or vice versa), via the Kramers-Kronig relation, only in the absence of data. We use these refractive indices as input for LX-MIE to compute the absorption, scattering and extinction efficiencies of all 32 condensate species. Finally, we use a three-parameter function to provide convenient fits to the shape of the extinction efficiency curve. We show that the errors associated with these simple fits in the Wide Field Camera 3 (WFC3), J, H and K wavebands are $\sim 10\%$. These fits allow for the extinction cross section or opacity of the condensate species to be easily included in retrieval analyses of transmission spectra. We discuss prospects for future experimental work. The compilation of the optical constants and LX-MIE are publicly available as part of the open-source Exoclime Simulation Platform (http://www.exoclime.org).

We analyze the interiors of HD~219134~b and c, which are among the coolest super Earths detected thus far. Without using spectroscopic measurements, we aim at constraining if the possible atmospheres are hydrogen-rich or hydrogen-poor. In a first step, we employ a full probabilistic Bayesian inference analysis in order to rigorously quantify the degeneracy of interior parameters given the data of mass, radius, refractory element abundances, semi-major axes, and stellar irradiation. We obtain constraints on structure and composition for core, mantle, ice layer, and atmosphere. In a second step, we aim to draw conclusions on the nature of possible atmospheres by considering atmospheric escape. Specifically, we compare the actual possible atmospheres to a threshold thickness above which a primordial (H$_2$-dominated) atmosphere can be retained against evaporation over the planet's lifetime. The best constrained parameters are the individual layer thicknesses. The maximum radius fraction of possible atmospheres are 0.18 and 0.13 $R$ (radius), for planets b and c, respectively. These values are significantly smaller than the threshold thicknesses of primordial atmospheres: 0.28 and 0.19 $R$, respectively. Thus, the possible atmospheres of planets b and c are unlikely to be H$_2$-dominated. However, whether possible volatile layers are made of gas or liquid/solid water cannot be uniquely determined. Our main conclusions are: (1) the possible atmospheres for planets b and c are enriched and thus possibly secondary in nature, and (2) both planets may contain a gas layer, whereas the layer of HD 219134 b must be larger. HD 219134 c can be rocky.

The M dwarf GJ 436 hosts a transiting warm Neptune known to experience atmospheric escape. Previous observations revealed the presence of a giant hydrogen exosphere transiting the star for more than 5 h, and absorbing up to 56% of the flux in the blue wing of the stellar Lyman-{\alpha} line of neutral hydrogen (H i Ly{\alpha}). The unexpected size of this comet-like exosphere prevented observing the full transit of its tail. In this Letter, we present new Ly{\alpha} observations of GJ 436 obtained with the Space Telescope Imaging Spectrograph (STIS) instrument onboard the Hubble Space Telescope. The stability of the Ly{\alpha} line over six years allowed us to combine these new observations with archival data sets, substantially expanding the coverage of the exospheric transit. Hydrogen atoms in the tail of the exospheric cloud keep occulting the star for 10-25 h after the transit of the planet, remarkably confirming a previous prediction based on 3D numerical simulations with the EVaporating Exoplanet code (EVE). This result strengthens the interpretation that the exosphere of GJ 436b is shaped by both radiative braking and charge exchanges with the stellar wind. We further report flux decreases of 15 +/- 2% and 47 +/- 10% in the red wing of the Ly{\alpha} line and in the line of ionised silicon (Si iii). Despite some temporal variability possibly linked with stellar activity, these two signals occur during the exospheric transit and could be of planetary origin. Follow-up observations will be required to assess the possibility that the redshifted Ly{\alpha} and Si iii absorption signatures arise from interactions between the exospheric flow and the magnetic field of the star.

We analyze the emission spectrum of the hot Jupiter WASP-12b using our HELIOS-R retrieval code and HELIOS-K opacity calculator. When interpreting Hubble and Spitzer data, the retrieval outcomes are found to be prior-dominated. When the prior distributions of the molecular abundances are assumed to be log-uniform, the volume mixing ratio of HCN is found to be implausibly high. A VULCAN chemical kinetics model of WASP-12b suggests that chemical equilibrium is a reasonable assumption even when atmospheric mixing is implausibly rigorous. Guided by (exo)planet formation theory, we set Gaussian priors on the elemental abundances of carbon, oxygen and nitrogen with the Gaussian peaks being centered on the measured C/H, O/H and N/H values of the star. By enforcing chemical equilibrium, we find substellar O/H and stellar to slightly superstellar C/H for the dayside atmosphere of WASP-12b. The superstellar carbon-to-oxygen ratio is just above unity, regardless of whether clouds are included in the retrieval analysis, consistent with Madhusudhan et al. (2011). Furthermore, whether a temperature inversion exists in the atmosphere depends on one's assumption for the Gaussian width of the priors. Our retrieved posterior distributions are consistent with the formation of WASP-12b in a solar-composition protoplanetary disk, beyond the water iceline, via gravitational instability or pebble accretion (without core erosion) and migration inwards to its present orbital location via a disk-free mechanism, and are inconsistent with both in-situ formation and core accretion with disk migration, as predicted by Madhusudhan et al. (2017). We predict that the interpretation of James Webb Space Telescope WASP-12b data will not be prior-dominated.

We present a novel generalization of the two-stream method of radiative transfer, which allows for the accurate treatment of radiative transfer in the presence of strong infrared scattering by aerosols. We prove that this generalization involves only a simple modification of the coupling coefficients and transmission functions in the hemispheric two-stream method. This modification originates from allowing the ratio of the first Eddington coefficients to depart from unity. At the heart of the method is the fact that this ratio may be computed once and for all over the entire range of values of the single-scattering albedo and scattering asymmetry factor. We benchmark our improved two-stream method by calculating the fraction of flux reflected by a single atmospheric layer (the reflectivity) and comparing these calculations to those performed using a 32-stream discrete-ordinates method. We further compare our improved two-stream method to the two-stream source function (16 streams) and delta-Eddington methods, demonstrating that it is often more accurate at the order-of-magnitude level. Finally, we illustrate its accuracy using a toy model of the early Martian atmosphere hosting a cloud layer composed of carbon-dioxide ice particles. The simplicity of implementation and accuracy of our improved two-stream method renders it suitable for implementation in three-dimensional general circulation models. In other words, our improved two-stream method has the ease of implementation of a standard two-stream method, but the accuracy of a 32-stream method.

We present wide-field, spatially and highly resolved spectroscopic observations of Balmer filaments in the northeastern rim of Tycho's supernova remnant in order to investigate the signal of cosmic-ray (CR) acceleration. The spectra of Balmer-dominated shocks (BDSs) have characteristic narrow (FWHM $\sim$ 10 kms$^{-1}$) and broad (FWHM $\sim$ 1000 kms$^{-1}$) H$\alpha$ components. CRs affect the H$\alpha$-line parameters: heating the cold neutrals in the interstellar medium results in broadening of the narrow H$\alpha$-line width beyond 20 kms$^{-1}$, but also in reduction of the broad H$\alpha$-line width due to energy being removed from the protons in the post-shock region. For the first time we show that the width of the narrow H$\alpha$ line, much larger than 20 kms$^{-1}$, is not a resolution or geometric effect nor a spurious result of a neglected intermediate (FWHM $\sim$ 100 kms$^{-1}$) component resulting from hydrogen atoms undergoing charge exchange with warm protons in the broad-neutral precursor. Moreover, we show that a narrow line width $\gg$ 20 kms$^{-1}$ extends across the entire NE rim, implying CR acceleration is ubiquitous, and making it possible to relate its strength to locally varying shock conditions. Finally, we find several locations along the rim, where spectra are significantly better explained (based on Bayesian evidence) by inclusion of the intermediate component, with a width of 180 kms$^{-1}$ on average.

We present H$\alpha$ spectroscopic observations and detailed modelling of the Balmer filaments in the supernova remnant Tycho. We used Galaxy H$\alpha$ Fabry-P\'erot Spectrometer on the William Herschel Telescope with a 3.4'$\times$3.4' field-of-view, 0.2" pixel scale and $\sigma_\rm{instr}=8.1$ km/s resolution at 1" seeing for $\sim10$ hours, resulting in 82 spatial-spectral bins that resolve the narrow H$\alpha$ line in the entire Tycho's northeastern rim. For the first time, we can mitigate artificial line broadening from unresolved differential motion, and probe H$\alpha$ emission parameters in varying shock and ambient medium conditions. Broad H$\alpha$ line remains unresolved within spectral coverage of 392 km/s. We employed Bayesian inference to obtain reliable parameter confidence intervals, and quantify the evidence for models with multiple line components. The median H$\alpha$ narrow-line full-width at half-maximum of all bins and models is $W_\rm{NL}=(54.8\pm1.8)$ km/s at the $95\%$ confidence level, varying within [35, 72] km/s between bins and clearly broadened compared to the intrinsic (thermal) $\approx20$ km/s. Possible line splits are accounted for, significant in $\approx18\%$ of the filament, and presumably due to remaining projection effects. We also find wide-spread evidence for intermediate-line emission of a broad-neutral precursor, with median $W_\rm{IL}=(180\pm14)$ km/s ($95\%$ confidence). Finally, we present a measurement of the remnant's systemic velocity, $V_\rm{LSR}=-34$ km/s, and map differential line-of-sight motions. Our results confirm the existence and interplay of shock precursors in Tycho's remnant. In particular, we show that suprathermal narrow-line emission is near-universal in Tycho and that, in absence of an alternative explanation, collisionless supernova remnant shocks constitute a viable acceleration source for Galactic TeV Cosmic-Ray protons.

We present an open-source retrieval code named HELIOS-Retrieval (hereafter HELIOS-R), designed to obtain chemical abundances and temperature-pressure profiles from inverting the measured spectra of exoplanetary atmospheres. In the current implementation, we use an exact solution of the radiative transfer equation, in the pure absorption limit, in our forward model, which allows us to analytically integrate over all of the outgoing rays (instead of performing Gaussian quadrature). Two chemistry models are considered: unconstrained chemistry (where the mixing ratios are treated as free parameters) and equilibrium chemistry (enforced via analytical formulae, where only the elemental abundances are free parameters). The nested sampling algorithm allows us to formally implement Occam's Razor based on a comparison of the Bayesian evidence between models. We perform a retrieval analysis on the measured spectra of the HR 8799b, c, d and e directly imaged exoplanets. Chemical equilibrium is disfavored by the Bayesian evidence for HR 8799b, c and d. We find supersolar C/O, C/H and O/H values for the outer HR 8799b and c exoplanets, while the inner HR 8799d and e exoplanets have substellar C/O, substellar C/H and superstellar O/H values. If these retrieved properties are representative of the bulk compositions of the exoplanets, then they are inconsistent with formation via gravitational instability (without late-time accretion) and consistent with a core accretion scenario in which late-time accretion of ices occurred differently for the inner and outer exoplanets. For HR 8799e, we find that spectroscopy in the K band is crucial for constraining C/O and C/H. HELIOS-R is publicly available as part of the Exoclimes Simulation Platform (ESP; www.exoclime.org).

Diagnostic potential of the spectral signatures of Raman scattering imprinted in planetary albedo spectra at short optical wavelengths has been demonstrated in research on Solar System planets and has recently been proposed as a probe of exoplanet atmospheres, complementary to albedo studies at longer wavelengths. Spectral features caused by Raman scattering offer insight into the properties of planetary atmospheres, such as the atmospheric depth, composition, and temperature, as well as the possibility of detecting and spectroscopically identifying spectrally inactive species, such as H$_2$ and N$_2$, in the visible wavelength range. Raman albedo features, however, depend on both the properties of the atmosphere and the shape of the incident stellar spectrum. Identical planetary atmospheres can produce very different albedo spectra depending on the spectral properties of the host star. Here we present a set of geometric albedo spectra calculated for atmospheres with H$_2$/He, N$_2$, and CO$_2$ composition, irradiated by different stellar types ranging from late A to late K stars. Prominent albedo features caused by Raman scattering appear at different wavelengths for different types of host stars. We investigate how absorption due to alkali elements, sodium and potassium, may affect the intensity of Raman features, and discuss the preferred observing strategies for detecting Raman features in future observations.

We have found many Earth-sized worlds but we have no way of determining if their surfaces are Earth-like. This makes it impossible to quantitatively compare habitability, and pretending we can risks damaging the field.

Terrestrial planets in the habitable zones (HZs) of low-mass stars and cool dwarfs have received significant scrutiny recently because their shorter orbital periods increase their chances of detection and characterization compared to planets around G-dwarfs. As these planets are likely tidal-locked, improved 3D numerical simulations of such planetary atmospheres are needed to guide target selection. Here we use a 3-D climate system model, updated with new water-vapor absorption coefficients derived from the HITRAN 2012 database, to study ocean covered planets at the inner edge of the HZ around late-M to mid-K stars ($2600$ K <= Teff <= 4500K). Our results indicate that these updated water-vapor coefficients result in significant warming compared to previous studies, so the inner HZ around M-dwarfs is not as close as suggested by earlier work. Assuming synchronously rotating planets with background 1 bar N2 atmospheres, we find that planets at the inner HZ of stars with Teff > 3000K undergo the classical "moist-greenhouse" (H2O mixing ratio > 10-3 in the stratosphere) at significantly lower surface temperature (~ 280K) in our 3-D model compared with 1-D climate models (~ 340K). This implies that some planets around low mass stars can simultaneously undergo water-loss and remain habitable. However, for star with Teff <= 3000K, planets at the inner HZ may directly transition to a runaway state, while bypassing the moist greenhouse water-loss entirely. We analyze transmission spectra of planets in a moist green-house regime, and find that there are several prominent H2O features within JWST instruments range. Thus, relying only upon standard Earth-analog spectra with 24-hour rotation period around M-dwarfs for habitability studies will miss the strong H2O features that one would expect to see on synchronously rotating planets around M-dwarf stars, with JWST.

Remote sensing of the atmospheres of distant worlds motivates a firm understanding of radiative transfer. In this review, we provide a pedagogical cookbook that describes the principal ingredients needed to perform a radiative transfer calculation and predict the spectrum of an exoplanet atmosphere, including solving the radiative transfer equation, calculating opacities (and chemistry), iterating for radiative equilibrium (or not), and adapting the output of the calculations to the astronomical observations. A review of the state of the art is performed, focusing on selected milestone papers. Outstanding issues, including the need to understand aerosols or clouds and elucidating the assumptions and caveats behind inversion methods, are discussed. A checklist is provided to assist referees/reviewers in their scrutiny of works involving radiative transfer. A table summarizing the methodology employed by past studies is provided.

The computation of transmission spectra is a central ingredient in the study of exoplanetary atmospheres. First, we revisit the theory of transmission spectra, unifying ideas from several workers in the literature. Transmission spectra lack an absolute normalization due to the a priori unknown value of a reference transit radius, which is tied to an unknown reference pressure. We show that there is a degeneracy between the uncertainty in the transit radius, the assumed value of the reference pressure (typically set to 10 bar) and the inferred value of the water abundance when interpreting a WFC3 transmission spectrum. Second, we show that the transmission spectra of isothermal atmospheres are nearly isobaric. We validate the isothermal, isobaric analytical formula for the transmission spectrum against full numerical calculations and show that the typical errors are ~0.1% (~10 ppm) within the WFC3 range of wavelengths for temperatures of 1500 K (or higher). Third, we generalize the previous expression for the transit radius to include a small temperature gradient. Finally, we analyze the measured WFC3 transmission spectrum of WASP-12b and demonstrate that we obtain consistent results with the retrieval approach of Kreidberg et al. (2015) if the reference transit radius and reference pressure are fixed to assumed values. The unknown functional relationship between the reference transit radius and reference pressure implies that it is the product of the water abundance and reference pressure that is being retrieved from the data, and not just the water abundance alone. This degeneracy leads to a limitation on how accurately we may extract molecular abundances from transmission spectra using WFC3 data alone. Finally, we compare our study to that of Griffith (2014) and discuss why the degeneracy was missed in previous retrieval studies. [abridged]

The TRAPPIST-1 system is the first transiting planet system found orbiting an ultra-cool dwarf star. At least seven planets similar to Earth in radius and in mass were previously found to transit this host star. Subsequently, TRAPPIST-1 was observed as part of the K2 mission and, with these new data, we report the measurement of an 18.77 d orbital period for the outermost planet, TRAPPIST-1h, which was unconstrained until now. This value matches our theoretical expectations based on Laplace relations and places TRAPPIST-1h as the seventh member of a complex chain, with three-body resonances linking every member. We find that TRAPPIST-1h has a radius of 0.727 Earth radii and an equilibrium temperature of 173 K. We have also measured the rotational period of the star at 3.3 d and detected a number of flares consistent with a low-activity, middle-aged, late M dwarf.

Space-based photometric surveys have discovered large numbers of planets transiting other stars, but these observe in a single band-pass and yield only the planet radius, orbital period, and transit duration. Information on the masses, compositions, and any atmospheres of these planets requires additional observations from the ground or space. The Transiting Exoplanet Survey Satellite (TESS) will yield thousands of planets around bright stars suitable for such follow-up. In the absence of spectroscopy or spectrophotometry from space, observations through the different pass-bands of multiple space telescopes provide some spectral information useful for identifying false positive signals, distinguishing between reflected light and thermal emission from hot Jupiters, and detecting Rayleigh scattering by planetary atmospheres. We calculated the expected difference in transit depths measured by the TESS and Characterizing Exoplanets Satellites (CHEOPS) missions, which will be more sensitive to redder and bluer optical wavelengths, respectively. The difference due to companion or background stars is small (<3% for main sequence companions) and likely to be negligible and undetectable. For only a few "hot" Jupiters, can combined photometry disambiguate between the reflected and thermal signals from planets. However, Rayleigh scattering by hazy atmospheres with particles sizes near 0.04 $\mu$m and at pressure altitudes above ~1 mbar can be detected for ~100 transiting planets, assuming every planet has such an atmosphere. Hazes with this characteristic particle size do not obscure observations at longer (near-infrared) wavelengths; CHEOPS follow-up of TESS-detected planets could thus identify candidates suitable for further study with the James Webb Space Telescope.

High-resolution optical spectroscopy during the transit of HD 189733b, a prototypical hot Jupiter, allowed the resolution of the Na I D sodium lines in the planet, giving access to the extreme conditions of the planet upper atmosphere. We have undertaken HEARTS, a spectroscopic survey of exoplanet upper atmospheres, to perform a comparative study of hot gas giants and determine how stellar irradiation affect them. Here, we report on the first HEARTS observations of the hot Saturn-mass planet WASP-49b. We observed the planet with the HARPS high-resolution spectrograph at ESO 3.6m telescope. We collected 126 spectra of WASP-49, covering three transits of WASP-49b. We analyzed and modeled the planet transit spectrum, while paying particular attention to the treatment of potentially spurious signals of stellar origin. We spectrally resolve the Na I D lines in the planet atmosphere and show that these signatures are unlikely to arise from stellar contamination. The large contrasts of $2.0\pm0.5\%$ (D$_2$) and $1.8\pm0.7\%$ (D$_1$) require the presence of hot neutral sodium ($2,950^{+400}_{-500}$ K) at high altitudes ($\sim$1.5 planet radius or $\sim$45,000 km). From estimating the cloudiness index of WASP-49b, we determine its atmosphere to be cloud free at the altitudes probed by the sodium lines. WASP-49b is close to the border of the evaporation desert and exhibits an enhanced thermospheric signature with respect to a farther-away planet such as HD 189733b.

We present an open-source and validated chemical kinetics code for studying hot exoplanetary atmospheres, which we name VULCAN. It is constructed for gaseous chemistry from 500 to 2500 K using a reduced C-H-O chemical network with about 300 reactions. It uses eddy diffusion to mimic atmospheric dynamics and excludes photochemistry. We have provided a full description of the rate coefficients and thermodynamic data used. We validate VULCAN by reproducing chemical equilibrium and by comparing its output versus the disequilibrium-chemistry calculations of Moses et al. and Rimmer & Helling. It reproduces the models of HD 189733b and HD 209458b by Moses et al., which employ a network with nearly 1600 reactions. We also use VULCAN to examine the theoretical trends produced when the temperature-pressure profile and carbon-to-oxygen ratio are varied. Assisted by a sensitivity test designed to identify the key reactions responsible for producing a specific molecule, we revisit the quenching approximation and find that it is accurate for methane but breaks down for acetylene, because the disequilibrium abundance of acetylene is not directly determined by transport-induced quenching, but is rather indirectly controlled by the disequilibrium abundance of methane. Therefore, we suggest that the quenching approximation should be used with caution and must always be checked against a chemical kinetics calculation. A one-dimensional model atmosphere with 100 layers, computed using VULCAN, typically takes several minutes to complete. VULCAN is part of the Exoclimes Simulation Platform (ESP; http://www.exoclime.net/) and publicly available at http://github.com/exoclime/VULCAN .

We aim to present a generalized Bayesian inference method for constraining interiors of super Earths and sub-Neptunes. Our methodology succeeds in quantifying the degeneracy and correlation of structural parameters for high dimensional parameter spaces. Specifically, we identify what constraints can be placed on composition and thickness of core, mantle, ice, ocean, and atmospheric layers given observations of mass, radius, and bulk refractory abundance constraints (Fe, Mg, Si) from observations of the host star's photospheric composition. We employed a full probabilistic Bayesian inference analysis that formally accounts for observational and model uncertainties. Using a Markov chain Monte Carlo technique, we computed joint and marginal posterior probability distributions for all structural parameters of interest. We included state-of-the-art structural models based on self-consistent thermodynamics of core, mantle, high-pressure ice, and liquid water. Furthermore, we tested and compared two different atmospheric models that are tailored for modeling thick and thin atmospheres, respectively. First, we validate our method against Neptune. Second, we apply it to synthetic exoplanets of fixed mass and determine the effect on interior structure and composition when (1) radius, (2) atmospheric model, (3) data uncertainties, (4) semi-major axes, (5) atmospheric composition (i.e., a priori assumption of enriched envelopes versus pure H/He envelopes), and (6) prior distributions are varied. Our main conclusions are: [...]

We present the open-source radiative transfer code named HELIOS, which is constructed for studying exoplanetary atmospheres. In its initial version, the model atmospheres of HELIOS are one-dimensional and plane-parallel, and the equation of radiative transfer is solved in the two-stream approximation with non-isotropic scattering. A small set of the main infrared absorbers is employed, computed with the opacity calculator HELIOS-K and combined using a correlated-$k$ approximation. The molecular abundances originate from validated analytical formulae for equilibrium chemistry. We compare HELIOS with the work of Miller-Ricci & Fortney using a model of GJ 1214b, and perform several tests, where we find: model atmospheres with single-temperature layers struggle to converge to radiative equilibrium; $k$-distribution tables constructed with $\gtrsim 0.01$ cm$^{-1}$ resolution in the opacity function ($ \lesssim 10^3$ points per wavenumber bin) may result in errors $\gtrsim 1$-10 % in the synthetic spectra; and a diffusivity factor of 2 approximates well the exact radiative transfer solution in the limit of pure absorption. We construct "null-hypothesis" models (chemical equilibrium, radiative equilibrium and solar element abundances) for 6 hot Jupiters. We find that the dayside emission spectra of HD 189733b and WASP-43b are consistent with the null hypothesis, while it consistently under-predicts the observed fluxes of WASP-8b, WASP-12b, WASP-14b and WASP-33b. We demonstrate that our results are somewhat insensitive to the choice of stellar models (blackbody, Kurucz or PHOENIX) and metallicity, but are strongly affected by higher carbon-to-oxygen ratios. The code is publicly available as part of the Exoclimes Simulation Platform (ESP; exoclime.net).

Due to its proximity, SN 1987A offers a unique opportunity to directly observe the geometry of a stellar explosion as it unfolds. Here we present spectral and imaging observations of SN 1987A obtained ~10,000 days after the explosion with HST/STIS and VLT/SINFONI at optical and near-infrared wavelengths. These observations allow us to produce the most detailed 3D map of H-alpha to date, the first 3D maps for [Ca II] \lambda \lambda 7292, 7324, [O I] \lambda \lambda 6300, 6364 and Mg II \lambda \lambda 9218, 9244, as well as new maps for [Si I]+[Fe II] 1.644 \mu m and He I 2.058 \mu m. A comparison with previous observations shows that the [Si I]+[Fe II] flux and morphology have not changed significantly during the past ten years, providing evidence that it is powered by 44Ti. The time-evolution of H-alpha shows that it is predominantly powered by X-rays from the ring, in agreement with previous findings. All lines that have sufficient signal show a similar large-scale 3D structure, with a north-south asymmetry that resembles a broken dipole. This structure correlates with early observations of asymmetries, showing that there is a global asymmetry that extends from the inner core to the outer envelope. On smaller scales, the two brightest lines, H-alpha and [Si I]+[Fe II] 1.644 \mu m, show substructures at the level of ~ 200 - 1000 km/s and clear differences in their 3D geometries. We discuss these results in the context of explosion models and the properties of dust in the ejecta.

An important source of opacity in exoplanet atmospheres at short visible and near-UV wavelengths is Rayleigh scattering of light on molecules. It is accompanied by a related, albeit weaker process -- Raman scattering. We analyze the signatures of Raman scattering imprinted in the reflected light and the geometric albedo of exoplanets, which could provide information about atmospheric properties. Raman scattering affects the geometric albedo spectra of planets in following ways. Firstly, it causes filling-in of strong absorption lines in the incident radiation, thus producing sharp peaks in the albedo. Secondly, it shifts the wavelengths of spectral features in the reflected light causing the so-called Raman ghost lines. Raman scattering can also cause a broadband reduction of the albedo due to wavelength shifting of a stellar spectrum with red spectral index. Observing the Raman peaks in the albedo could be used to measure the column density of gas, thus providing constrains on the presence of clouds in the atmosphere. Observing the Raman ghost lines could be used to spectroscopically identify the main scatterer in the atmosphere, even molecules like H$_2$ or N$_2$ that do not have prominent spectral signatures in the optical wavelength range. If detected, ghost lines could also provide information about the temperature of the atmosphere. In this paper we investigate the effects of Raman scattering in hydrogen- and nitrogen-dominated atmospheres. We analyze the feasibility of detecting the signatures of Raman scattering with the existing and future observational facilities, and of using these signatures as probes of exoplanetary atmospheres.

We have designed and developed, from scratch, a global circulation model named THOR that solves the three-dimensional non-hydrostatic Euler equations. Our general approach lifts the commonly used assumptions of a shallow atmosphere and hydrostatic equilibrium. We solve the "pole problem" (where converging meridians on a sphere lead to increasingly smaller time steps near the poles) by implementing an icosahedral grid. Irregularities in the grid, which lead to grid imprinting, are smoothed using the "spring dynamics" technique. We validate our implementation of spring dynamics by examining calculations of the divergence and gradient of test functions. To prevent the computational time step from being bottlenecked by having to resolve sound waves, we implement a split-explicit method together with a horizontally explicit and vertically implicit integration. We validate our global circulation model by reproducing the Earth and also the hot Jupiter-like benchmark tests. THOR was designed to run on Graphics Processing Units (GPUs), which allows for physics modules (radiative transfer, clouds, chemistry) to be added in the future, and is part of the open-source Exoclimes Simulation Platform (ESP; www.exoclime.org).

We present novel, analytical, equilibrium-chemistry formulae for the abundances of molecules in hot exoplanetary atmospheres that include the carbon, oxygen and nitrogen networks. Our hydrogen-dominated solutions involve acetylene (C$_2$H$_2$), ammonia (NH$_3$), carbon dioxide (CO$_2$), carbon monoxide (CO), ethylene (C$_2$H$_4$), hydrogen cyanide (HCN), methane (CH$_4$), molecular nitrogen (N$_2$) and water (H$_2$O). By considering only the gas phase, we prove that the mixing ratio of carbon monoxide is governed by a decic equation (polynomial equation of degree 10). We validate our solutions against numerical calculations of equilibrium chemistry that perform Gibbs free energy minimization and demonstrate that they are accurate at the $\sim 1\%$ level for temperatures from 500 to 3000 K. In hydrogen-dominated atmospheres, the ratio of abundances of HCN to CH$_4$ is nearly constant across a wide range of carbon-to-oxygen ratios, which makes it a robust diagnostic of the metallicity in the gas phase. Our validated formulae allow for the convenient benchmarking of chemical kinetics codes and provide an efficient way of enforcing chemical equilibrium in atmospheric retrieval calculations.

We present a dimensionless index that quantifies the degree of cloudiness of the atmosphere of a transiting exoplanet. Our cloudiness index is based on measuring the transit radii associated with the line center and wing of the sodium or potassium line. In deriving this index, we revisited the algebraic formulae for inferring the isothermal pressure scale height from transit measurements. We demonstrate that the formulae of Lecavelier et al. and Benneke & Seager are identical: the former is inferring the temperature while assuming a value for the mean molecular mass and the latter is inferring the mean molecular mass while assuming a value for the temperature. More importantly, these formulae cannot be used to distinguish between cloudy and cloudfree atmospheres. We derive values of our cloudiness index for a small sample of 7 hot Saturns/Jupiters taken from Sing et al. We show that WASP-17b, WASP-31b and HAT-P-1b are nearly cloudfree at visible wavelengths. We find the tentative trend that more irradiated atmospheres tend to have less clouds consisting of sub-micron-sized particles. We also derive absolute sodium and/or potassium abundances $\sim 10^2$ cm$^{-3}$ for WASP-17b, WASP-31b and HAT-P-1b (and upper limits for the other objects). Higher-resolution measurements of both the sodium and potassium lines, for a larger sample of exoplanetary atmospheres, are needed to confirm or refute this trend.

The James Webb Space Telescope will revolutionize transiting exoplanet atmospheric science due to its capability for continuous, long-duration observations and its larger collecting area, spectral coverage, and spectral resolution compared to existing space-based facilities. However, it is unclear precisely how well JWST will perform and which of its myriad instruments and observing modes will be best suited for transiting exoplanet studies. In this article, we describe a prefatory JWST Early Release Science (ERS) program that focuses on testing specific observing modes to quickly give the community the data and experience it needs to plan more efficient and successful future transiting exoplanet characterization programs. We propose a multi-pronged approach wherein one aspect of the program focuses on observing transits of a single target with all of the recommended observing modes to identify and understand potential systematics, compare transmission spectra at overlapping and neighboring wavelength regions, confirm throughputs, and determine overall performances. In our search for transiting exoplanets that are well suited to achieving these goals, we identify 12 objects (dubbed "community targets") that meet our defined criteria. Currently, the most favorable target is WASP-62b because of its large predicted signal size, relatively bright host star, and location in JWST's continuous viewing zone. Since most of the community targets do not have well-characterized atmospheres, we recommend initiating preparatory observing programs to determine the presence of obscuring clouds/hazes within their atmospheres. Measurable spectroscopic features are needed to establish the optimal resolution and wavelength regions for exoplanet characterization. Other initiatives from our proposed ERS program include testing the instrument brightness limits and performing phase-curve observations.(Abridged)

Over the past decade, observations of giant exoplanets (Jupiter-size) have provided key insights into their atmospheres, but the properties of lower-mass exoplanets (sub-Neptune) remain largely unconstrained because of the challenges of observing small planets. Numerous efforts to observe the spectra of super-Earths (exoplanets with masses of one to ten times that of Earth) have so far revealed only featureless spectra. Here we report a longitudinal thermal brightness map of the nearby transiting super-Earth 55 Cancri e revealing highly asymmetric dayside thermal emission and a strong day-night temperature contrast. Dedicated space-based monitoring of the planet in the infrared revealed a modulation of the thermal flux as 55 Cancri e revolves around its star in a tidally locked configuration. These observations reveal a hot spot that is located 41 +- 12 degrees east of the substellar point (the point at which incident light from the star is perpendicular to the surface of the planet). From the orbital phase curve, we also constrain the nightside brightness temperature of the planet to 1380 +- 400 kelvin and the temperature of the warmest hemisphere (centred on the hot spot) to be about 1300 kelvin hotter (2700 +- 270 kelvin) at a wavelength of 4.5 microns, which indicates inefficient heat redistribution from the dayside to the nightside. Our observations are consistent with either an optically thick atmosphere with heat recirculation confined to the planetary dayside, or a planet devoid of atmosphere with low-viscosity magma flows at the surface.

Hunting for Alien Habitats Without Defining Life?

General circulation models of the atmosphere of hot Jupiters have shown the existence of a supersonic eastward equatorial jet. In this paper, we investigate the effects of compressibility on the atmospheric dynamics by solving the standard Euler equations. This is done by means of a series of simulations performed in the framework of the equatorial beta-plane approximation using the finite volume shock-capturing code RAMSES. At low resolution, we recover the classical results described in the literature: we find a strong and steady supersonic equatorial jet of a few km/s that displays no signature of shocks. We next show that the jet zonal velocity depends significantly on the grid meridional resolution. When that resolution is fine enough to properly resolve the jet, the latter is subject to a Kelvin-Helmholtz instability. The jet zonal mean velocity displays regular oscillations with a typical timescale of few days and a significant amplitude of about 15% of the jet velocity. We also find compelling evidence for the development of a vertical shear instability at pressure levels of a few bars. It seems to be responsible for an increased downward kinetic energy flux, significantly affecting the temperature of the deep atmosphere, and appears to act as a form of drag on the equatorial jet. This instability also creates velocity fluctuations that propagate upward and steepen into weak shocks at pressure levels of a few mbars. We conclude that hot Jupiter equatorial jets are potentially unstable to both a barotropic Kelvin-Helmholtz instability and a vertical shear instability. Upon confirmation using more realistic models, both instabilities could result in significant time variability of the atmospheric winds, may provide a small scale dissipation mechanism in the flow, and might have consequences for the internal evolution of hot Jupiters.

Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr)

Over the duration of the Kepler mission, KIC8462852 was observed to undergo irregularly shaped, aperiodic dips in flux of up to $\sim 20$\%. The dipping activity can last for between 5 and 80 days. We characterize the object with high-resolution spectroscopy, spectral energy distribution fitting, radial velocity measurements, high-resolution imaging, and Fourier analyses of the Kepler light curve. We determine that KIC8462852 is a typical main-sequence F3 V star that exhibits no significant IR excess, and has no very close interacting companions. In this paper, we describe various scenarios to explain the dipping events observed in the Kepler light curve. We confirm that the dipping signals in the data are not caused by any instrumental or data processing artifact, and thus are astrophysical in origin. We construct scenario-independent constraints on the size and location of a body in the system that is needed to reproduce the observations. We deliberate over several assorted stellar and circumstellar astrophysical scenarios, most of which have problems explaining the data in hand. By considering the observational constraints on dust clumps in orbit around a normal main-sequence star, we conclude that the scenario most consistent with the data in hand is the passage of a family of exocomet or planetesimal fragments, all of which are associated with a single previous break-up event, possibly caused by tidal disruption or thermal processing. The minimum total mass associated with these fragments likely exceeds $10^{-6}$~\mearth, corresponding to an original rocky body of $>100$~km in diameter. We discuss the necessity of future observations to help interpret the system.

Optical phase curves have become one of the common probes of exoplanetary atmospheres, but the information they encode has not been fully elucidated. Building on a diverse body of work, we upgrade the Flexible Modeling System (FMS) to include scattering in the two-stream, dual-band approximation and generate plausible, three-dimensional structures of irradiated atmospheres to study the radiative effects of aerosols or condensates. In the optical, we treat the scattering of starlight using a generalisation of Beer's law that allows for a finite Bond albedo to be prescribed. In the infrared, we implement the two-stream solutions and include scattering via an infrared scattering parameter. We present a suite of four-parameter general circulation models for Kepler-7b and demonstrate that its climatology is expected to be robust to variations in optical and infrared scattering. The westward and eastward shifts of the optical and infrared phase curves, respectively, are shown to be robust outcomes of the simulations. Assuming micron-sized particles and a simplified treatment of local brightness, we further show that the peak offset of the optical phase curve is sensitive to the composition of the aerosols or condensates. However, to within the measurement uncertainties, we cannot distinguish between aerosols made of silicates (enstatite or forsterite), iron, corundum or titanium oxide, based on a comparison to the measured peak offset ($41^\circ \pm 12^\circ$) of the optical phase curve of Kepler-7b. Measuring high-precision optical phase curves will provide important constraints on the atmospheres of cloudy exoplanets and reduce degeneracies in interpreting their infrared spectra.

We present a comprehensive study of the abundance of carbon dioxide in exoplanetary atmospheres in hot, hydrogen-dominated atmospheres. We construct novel analytical models of systems in chemical equilibrium that include carbon monoxide, carbon dioxide, water, methane and acetylene and relate the equilibrium constants of the chemical reactions to temperature and pressure via the tabulated Gibbs free energies. We prove that such chemical systems may be described by a quintic equation for the mixing ratio of methane. By examining the abundances of these molecules across a broad range of temperatures (spanning equilibrium temperatures from 600 to 2500 K), pressures (via temperature-pressure profiles that explore albedo and opacity variations) and carbon-to-oxygen ratios, we conclude that carbon dioxide is subdominant compared to carbon monoxide and water. Atmospheric mixing does not alter this conclusion if carbon dioxide is subdominant everywhere in the atmosphere. Carbon dioxide and carbon monoxide may attain comparable abundances if the metallicity is greatly enhanced, but this property is negated by temperatures above 1000 K. For hydrogen-dominated atmospheres, our generic result has the implication that retrieval studies may wish to set the subdominance of carbon dioxide as a prior of the calculation and not let its abundance completely roam free as a fitting parameter, because it directly affects the inferred value of the carbon-to-oxygen ratio and may produce unphysical conclusions. We discuss the relevance of these implications for the hot Jupiter WASP-12b and suggest that some of the previous results are chemically impossible. The relative abundance of carbon dioxide to acetylene is potentially a sensitive diagnostic of the carbon-to-oxygen ratio.

We present a self-consistent formalism for computing and understanding the atmospheric chemistry of exoplanets from the viewpoint of an astrophysicist. Starting from the first law of thermodynamics, we demonstrate that the van't Hoff equation (which describes the equilibrium constant), Arrhenius equation (which describes the rate coefficients) and procedures associated with the Gibbs free energy (minimisation, rescaling) have a common physical and mathematical origin. We address an ambiguity associated with the equilibrium constant, which is used to relate the forward and reverse rate coefficients, and restate its two definitions. By necessity, one of the equilibrium constants must be dimensionless and equate to an exponential function involving the Gibbs free energy, while the other is a ratio of rate coefficients and must therefore possess physical units. We demonstrate that the Arrhenius equation takes on a functional form that is more general than previously stated without recourse to tagging on ad hoc functional forms. Finally, we derive analytical models of chemical systems, in equilibrium, with carbon, hydrogen and oxygen. We include acetylene and are able to reproduce several key trends, versus temperature and carbon-to-oxygen ratio, published in the literature. The rich variety of behavior that mixing ratios exhibit as a function of the carbon-to-oxygen ratio is merely the outcome of stoichiometric book-keeping and not the direct consequence of temperature or pressure variations.

Ocean planets are volatile rich planets, not present in our Solar System, which are thought to be dominated by deep, global oceans. This results in the formation of high-pressure water ice, separating the planetary crust from the liquid ocean and, thus, also from the atmosphere. Therefore, instead of a carbonate-silicate cycle like on the Earth, the atmospheric carbon dioxide concentration is governed by the capability of the ocean to dissolve carbon dioxide (CO2). In our study, we focus on the CO2 cycle between the atmosphere and the ocean which determines the atmospheric CO2 content. The atmospheric amount of CO2 is a fundamental quantity for assessing the potential habitability of the planet's surface because of its strong greenhouse effect, which determines the planetary surface temperature to a large degree. In contrast to the stabilising carbonate-silicate cycle regulating the long-term CO2 inventory of the Earth atmosphere, we find that the CO2 cycle feedback on ocean planets is negative and has strong destabilising effects on the planetary climate. By using a chemistry model for oceanic CO2 dissolution and an atmospheric model for exoplanets, we show that the CO2 feedback cycle can severely limit the extension of the habitable zone for ocean planets.

What strange new worlds will our next-generation telescopes find?

We present an ultrafast opacity calculator that we name HELIOS-K. It takes a line list as an input, computes the shape of each spectral line and provides an option for grouping an enormous number of lines into a manageable number of bins. We implement a combination of Algorithm 916 and Gauss-Hermite quadrature to compute the Voigt profile, write the code in CUDA and optimise the computation for graphics processing units (GPUs). We restate the theory of the k-distribution method and use it to reduce $\sim 10^5$ to $10^8$ lines to $\sim 10$ to $10^4$ wavenumber bins, which may then be used for radiative transfer, atmospheric retrieval and general circulation models. The choice of line-wing cutoff for the Voigt profile is a significant source of error and affects the value of the computed flux by $\sim 10\%$. This is an outstanding physical (rather than computational) problem, due to our incomplete knowledge of pressure broadening of spectral lines in the far line wings. We emphasize that this problem remains regardless of whether one performs line-by-line calculations or uses the k-distribution method and affects all calculations of exoplanetary atmospheres requiring the use of wavelength-dependent opacities. We elucidate the correlated-k approximation and demonstrate that it applies equally to inhomogeneous atmospheres with a single atomic/molecular species or homogeneous atmospheres with multiple species. Using a NVIDIA K20 GPU, HELIOS-K is capable of computing an opacity function with $\sim 10^5$ spectral lines in $\sim 1$ second and is publicly available as part of the Exoclimes Simulation Platform (ESP; www.exoclime.org).

We present imaging and spectroscopic observations with HST and VLT of the ring of SN 1987A from 1994 to 2014. After an almost exponential increase of the shocked emission from the hotspots up to day ~8,000 (~2009), both this and the unshocked emission are now fading. From the radial positions of the hotspots we see an acceleration of these up to 500-1000 km/s, consistent with the highest spectroscopic shock velocities from the radiative shocks. In the most recent observations (2013 and 2014), we find several new hotspots outside the inner ring, excited by either X-rays from the shocks or by direct shock interaction. All of these observations indicate that the interaction with the supernova ejecta is now gradually dissolving the hotspots. We predict, based on the observed decay, that the inner ring will be destroyed by ~2025.

We present an inversion method based on Bayesian analysis to constrain the interior structure of terrestrial exoplanets, in the form of chemical composition of the mantle and core size. Specifically, we identify what parts of the interior structure of terrestrial exoplanets can be determined from observations of mass, radius, and stellar elemental abundances. We perform a full probabilistic inverse analysis to formally account for observational and model uncertainties and obtain confidence regions of interior structure models. This enables us to characterize how model variability depends on data and associated uncertainties. We test our method on terrestrial solar system planets and find that our model predictions are consistent with independent estimates. Furthermore, we apply our method to synthetic exoplanets up to 10 Earth masses and up to 1.7 Earth radii as well as to exoplanet Kepler-36b. Importantly, the inversion strategy proposed here provides a framework for understanding the level of precision required to characterize the interior of exoplanets. Our main conclusions are: (1) observations of mass and radius are sufficient to constrain core size; (2) stellar elemental abundances (Fe, Si, Mg) are key constraints to reduce degeneracy in interior structure models and to constrain mantle composition; (3) the inherent degeneracy in determining interior structure from mass and radius observations does not only depend on measurement accuracies but also on the actual size and density of the exoplanet. We argue that precise observations of stellar elemental abundances are central in order to place constraints on planetary bulk composition and to reduce model degeneracy. [...]

WASP-80b is a missing link in the study of exo-atmospheres. It falls between the warm Neptunes and the hot Jupiters and is amenable for characterisation, thanks to its host star's properties. We observed the planet through transit and during occultation with Warm Spitzer. Combining our mid-infrared transits with optical time series, we find that the planet presents a transmission spectrum indistinguishable from a horizontal line. In emission, WASP-80b is the intrinsically faintest planet whose dayside flux has been detected in both the 3.6 and 4.5 $\mu$m Spitzer channels. The depths of the occultations reveal that WASP-80b is as bright and as red as a T4 dwarf, but that its temperature is cooler. If planets go through the equivalent of an L-T transition, our results would imply this happens at cooler temperatures than for brown dwarfs. Placing WASP-80b's dayside into a colour-magnitude diagram, it falls exactly at the junction between a blackbody model and the T-dwarf sequence; we cannot discern which of those two interpretations is the more likely. Flux measurements on other planets with similar equilibrium temperatures are required to establish whether irradiated gas giants, like brown dwarfs, transition between two spectral classes. An eventual detection of methane absorption in transmission would also help lift that degeneracy. We obtained a second series of high-resolution spectra during transit, using HARPS. We reanalyse the Rossiter-McLaughlin effect. The data now favour an aligned orbital solution and a stellar rotation nearly three times slower than stellar line broadening implies. A contribution to stellar line broadening, maybe macroturbulence, is likely to have been underestimated for cool stars, whose rotations have therefore been systematically overestimated. [abridged]

We present a theory for interpreting the sodium lines detected in transmission spectra of exoplanetary atmospheres. Previous analyses employed the isothermal approximation and dealt only with the transit radius. By recognising the absorption depth and the transit radius as being independent observables, we develop a theory for jointly interpreting both quantities, which allows us to infer the temperatures and number densities associated with the sodium lines. We are able to treat a non-isothermal situation with a constant temperature gradient. Our novel diagnostics take the form of simple-to-use algebraic formulae and require measurements of the transit radii (and their corresponding absorption depths) at line center and in the line wing for both sodium lines. We apply our diagnostics to the HARPS data of HD 189733b, confirm the upper atmospheric heating reported by Huitson et al. (2012), derive a temperature gradient of $0.4376 \pm 0.0154$ K km$^{-1}$ and find densities $\sim 1$ to $10^4$ cm$^{-3}$.

Exoplanet science is now in its full expansion, particularly after the CoRoT and Kepler space missions that led us to the discovery of thousands of extra-solar planets. The last decade has taught us that UV observations play a major role in advancing our understanding of planets and of their host stars, but the necessary UV observations can be carried out only by HST, and this is going to be the case for many years to come. It is therefore crucial to build a treasury data archive of UV exoplanet observations formed by a dozen "golden systems" for which observations will be available from the UV to the infrared. Only in this way we will be able to fully exploit JWST observations for exoplanet science, one of the key JWST science case.

We present new {\it Hubble Space Telescope} images of high-velocity H-$\alpha$ and Lyman-$\alpha$ emission in the outer debris of SN~1987A. The H-$\alpha$ images are dominated by emission from hydrogen atoms crossing the reverse shock. For the first time we observe emission from the reverse shock surface well above and below the equatorial ring, suggesting a bipolar or conical structure perpendicular to the ring plane. Using the H$\alpha$ imaging, we measure the mass flux of hydrogen atoms crossing the reverse shock front, in the velocity intervals ($-$7,500~$<$~$V_{obs}$~$<$~$-$2,800 km s$^{-1}$) and (1,000~$<$~$V_{obs}$~$<$~7,500 km s$^{-1}$), $\dot{M_{H}}$ = 1.2~$\times$~10$^{-3}$ M$_{\odot}$ yr$^{-1}$. We also present the first Lyman-$\alpha$ imaging of the whole remnant and new $Chandra$ X-ray observations. Comparing the spatial distribution of the Lyman-$\alpha$ and X-ray emission, we observe that the majority of the high-velocity Lyman-$\alpha$ emission originates interior to the equatorial ring. The observed Lyman-$\alpha$/H-$\alpha$ photon ratio, $\langle$$R(L\alpha / H\alpha)$$\rangle$ $\approx$~17, is significantly higher than the theoretically predicted ratio of $\approx$ 5 for neutral atoms crossing the reverse shock front. We attribute this excess to Lyman-$\alpha$ emission produced by X-ray heating of the outer debris. The spatial orientation of the Lyman-$\alpha$ and X-ray emission suggests that X-ray heating of the outer debris is the dominant Lyman-$\alpha$ production mechanism in SN 1987A at this phase in its evolution.

Supernova remnants are among the most spectacular examples of astrophysical pistons in our cosmic neighborhood. The gas expelled by the supernova explosion is launched with velocities ~ 1000 km/s into the ambient, tenuous interstellar medium, producing shocks that excite hydrogen lines. We have used an optical integral-field spectrograph to obtain high-resolution spatial-spectral maps that allow us to study in detail the shocks in the northwestern rim of supernova 1006. The two-component H\alpha\ line is detected at 133 sky locations. Variations in the broad line widths and the broad-to-narrow line intensity ratios across tens of atomic mean free paths suggest the presence of suprathermal protons, the potential seed particles for generating high-energy cosmic-rays.

We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior) and solutions for the temperature-pressure profiles. Generally, the problem is mathematically under-determined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat and the properties of scattering both in optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing and incoming fluxes in the convective regime.

The characterization of exoplanetary atmospheres has come of age in the last decade, as astronomical techniques now allow for albedos, chemical abundances, temperature profiles and maps, rotation periods and even wind speeds to be measured. Atmospheric dynamics sets the background state of density, temperature and velocity that determines or influences the spectral and temporal appearance of an exoplanetary atmosphere. Hot exoplanets are most amenable to these characterization techniques; in the present review, we focus on highly-irradiated, large exoplanets (the "hot Jupiters"), as astronomical data begin to confront theoretical questions. We summarize the basic atmospheric quantities inferred from the astronomical observations. We review the state of the art by addressing a series of current questions and look towards the future by considering a separate set of exploratory questions. Attaining the next level of understanding will require a concerted effort of constructing multi-faceted, multi-wavelength datasets for benchmark objects. Understanding clouds presents a formidable obstacle, as they introduce degeneracies into the interpretation of spectra, yet their properties and existence are directly influenced by atmospheric dynamics. Confronting general circulation models with these multi-faceted, multi-wavelength datasets will help us understand these and other degeneracies. The coming decade will witness a decisive confrontation of theory and simulation by the next generation of astronomical data.

We present high resolution transmission spectra of giant planet atmospheres from a coupled 3-D atmospheric dynamics and transmission spectrum model that includes Doppler shifts which arise from winds and planetary motion. We model jovian planets covering more than two orders of magnitude in incident flux, corresponding to planets with 0.9 to 55 day orbital periods around solar-type stars. The results of our 3-D dynamical models reveal certain aspects of high resolution transmission spectra that are not present in simple 1-D models. We find that the hottest planets experience strong substellar to anti-stellar (SSAS) winds, resulting in transmission spectra with net blue shifts of up to 3 km s$^{-1}$, whereas less irradiated planets show almost no net Doppler shifts. Compared to 1-D models, peak line strengths are significantly reduced for the hottest atmospheres owing to Doppler broadening from a combination of rotation (which is faster for close-in planets under the assumption of tidal locking) and atmospheric winds. Finally, high resolution transmission spectra may be useful in studying the atmospheres of exoplanets with optically thick clouds since line cores for very strong transitions should remain optically thick to very high altitude. High resolution transmission spectra are an excellent observational test for the validity of 3-D atmospheric dynamics models, because they provide a direct probe of wind structures and heat circulation. Ground-based exoplanet spectroscopy is currently on the verge of being able to verify some of our modeling predictions, most notably the dependence of SSAS winds on insolation. We caution that interpretation of high resolution transmission spectra based on 1-D atmospheric models may be inadequate, as 3-D atmospheric motions can produce a noticeable effect on the absorption signatures.

Within the context of exoplanetary atmospheres, we present a comprehensive linear analysis of forced, damped, magnetized shallow water systems, exploring the effects of dimensionality, geometry (Cartesian, pseudo-spherical and spherical), rotation, magnetic tension and hydrodynamic and magnetic sources of friction. Across a broad range of conditions, we find that the key governing equation for atmospheres and quantum harmonic oscillators are identical, even when forcing (stellar irradiation), sources of friction (molecular viscosity, Rayleigh drag and magnetic drag) and magnetic tension are included. The global atmospheric structure is largely controlled by a single, key parameter that involves the Rossby and Prandtl numbers. This near-universality breaks down when either molecular viscosity or magnetic drag acts non-uniformly across latitude or a poloidal magnetic field is present, suggesting that these effects will introduce qualitative changes to the familiar chevron-shaped feature witnessed in simulations of atmospheric circulation. We also find that hydrodynamic and magnetic sources of friction have dissimilar phase signatures and affect the flow in fundamentally different ways, implying that using Rayleigh drag to mimic magnetic drag is inaccurate. We exhaustively lay down the theoretical formalism (dispersion relations, governing equations and time-dependent wave solutions) for a broad suite of models. In all situations, we derive the steady state of an atmosphere, which is relevant to interpreting infrared phase and eclipse maps of exoplanetary atmospheres. We elucidate a pinching effect that confines the atmospheric structure to be near the equator. Our suite of analytical models may be used to decisively develop physical intuition and as a reference point for three-dimensional, magnetohydrodynamic (MHD) simulations of atmospheric circulation.

A number of observations have shown that Rayleigh scattering by aerosols dominates the transmission spectrum of HD 189733b at wavelengths shortward of 1 $\mu$m. In this study, we retrieve a range of aerosol distributions consistent with transmission spectroscopy between 0.3-24 $\mu$m that were recently re-analyzed by Pont et al. (2013). To constrain the particle size and the optical depth of the aerosol layer, we investigate the degeneracies between aerosol composition, temperature, planetary radius, and molecular abundances that prevent unique solutions for transit spectroscopy. Assuming that the aerosol is composed of MgSiO$_3$, we suggest that a vertically uniform aerosol layer over all pressures with a monodisperse particle size smaller than about 0.1 $\mu$m and an optical depth in the range 0.002-0.02 at 1 $\mu$m provides statistically meaningful solutions for the day/night terminator regions of HD 189733b. Generally, we find that a uniform aerosol layer provide adequate fits to the data if the optical depth is less than 0.1 and the particle size is smaller than 0.1 $\mu$m, irrespective of the atmospheric temperature, planetary radius, aerosol composition, and gaseous molecules. Strong constraints on the aerosol properties are provided by spectra at wavelengths shortward of 1 $\mu$m as well as longward of 8 $\mu$m, if the aerosol material has absorption features in this region. We show that these are the optimal wavelengths for quantifying the effects of aerosols, which may guide the design of future space observations. The present investigation indicates that the current data offer sufficient information to constrain some of the aerosol properties of HD189733b, but the chemistry in the terminator regions remains uncertain.

Is numerical mimicry a third way of establishing truth?

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

We present new visible and infrared observations of the hot Jupiter Kepler-7b to determine its atmospheric properties. Our analysis allows us to 1) refine Kepler-7b's relatively large geometric albedo of Ag=0.35+-0.02, 2) place upper limits on Kepler-7b thermal emission that remains undetected in both Spitzer bandpasses and 3) report a westward shift in the Kepler optical phase curve. We argue that Kepler-7b's visible flux cannot be due to thermal emission or Rayleigh scattering from H2 molecules. We therefore conclude that high altitude, optically reflective clouds located west from the substellar point are present in its atmosphere. We find that a silicate-based cloud composition is a possible candidate. Kepler-7b exhibits several properties that may make it particularly amenable to cloud formation in its upper atmosphere. These include a hot deep atmosphere that avoids a cloud cold trap, very low surface gravity to suppress cloud sedimentation, and a planetary equilibrium temperature in a range that allows for silicate clouds to potentially form in the visible atmosphere probed by Kepler. Our analysis does not only present evidence of optically thick clouds on Kepler-7b but also yields the first map of clouds in an exoplanet atmosphere.

Unlike previously explored relationships between the properties of hot Jovian atmospheres, the geometric albedo and the incident stellar flux do not exhibit a clear correlation, as revealed by our re-analysis of Q0 to Q14 Kepler data. If the albedo is primarily associated with the presence of clouds in these irradiated atmospheres, a holistic modeling approach needs to relate the following properties: the strength of stellar irradiation (and hence the strength and depth of atmospheric circulation), the geometric albedo (which controls both the fraction of starlight absorbed and the pressure level at which it is predominantly absorbed) and the properties of the embedded cloud particles (which determine the albedo). The anticipated diversity in cloud properties renders any correlation between the geometric albedo and the stellar flux to be weak and characterized by considerable scatter. In the limit of vertically uniform populations of scatterers and absorbers, we use an analytical model and scaling relations to relate the temperature-pressure profile of an irradiated atmosphere and the photon deposition layer and to estimate if a cloud particle will be lofted by atmospheric circulation. We derive an analytical formula for computing the albedo spectrum in terms of the cloud properties, which we compare to the measured albedo spectrum of HD 189733b by Evans et al. (2013). Furthermore, we show that whether an optical phase curve is flat or sinusoidal depends on whether the particles are small or large as defined by the Knudsen number. This may be an explanation for why Kepler-7b exhibits evidence for the longitudinal variation in abundance of condensates, while Kepler-12b shows no evidence for the presence of condensates, despite the incident stellar flux being similar for both exoplanets.

Directly-imaged exoplanets are unexplored laboratories for the application of the spectral and temperature retrieval method, where the chemistry and composition of their atmospheres are inferred from inverse modeling of the available data. As a pilot study, we focus on the extrasolar gas giant HR 8799b for which more than 50 data points are available. We upgrade our non-linear optimal estimation retrieval method to include a phenomenological model of clouds that requires the cloud optical depth and monodisperse particle size to be specified. Previous studies have focused on forward models with assumed values of the exoplanetary properties; there is no consensus on the best-fit values of the radius, mass, surface gravity and effective temperature of HR 8799b. We show that cloudfree models produce reasonable fits to the data if the atmosphere is of super-solar metallicity and non-solar elemental abundances. Intermediately cloudy models with moderate values of the cloud optical depth and micron-sized particles provide an equally reasonable fit to the data and require a lower mean molecular weight. We report our best-fit values for the radius, mass, surface gravity and effective temperature of HR 8799b. The mean molecular weight is about 3.8, while the carbon-to-oxygen ratio is about unity due to the prevalence of carbon monoxide. Our study emphasizes the need for robust claims about the nature of an exoplanetary atmosphere to be based on analyses involving both photometry and spectroscopy and inferred from beyond a few photometric data points, such as are typically reported for hot Jupiters.

We present integral-field spectroscopic observations with the VIMOS-IFU at the VLT of fast (2000-3000 km/s) Balmer-dominated shocks surrounding the northwestern rim of the remnant of supernova 1006. The high spatial and spectral resolution of the instrument enable us to show that the physical characteristics of the shocks exhibit a strong spatial variation over few atomic scale lengths across 133 sky locations. Our results point to the presence of a population of non-thermal protons (10-100 keV) which might well be the seed particles for generating high-energy cosmic rays. We also present observations of Tycho's supernova remnant taken with the narrow-band tunable filter imager OSIRIS at the GTC and the Fabry-Perot interferometer GHaFaS at the WHT to resolve respectively the broad and narrow H\alpha\ lines across a large part of the remnant.

We present a secondary eclipse observation for the hot Jupiter HD189733b across the wavelength range 290-570nm made using the Space Telescope Imaging Spectrograph on the Hubble Space Telescope. We measure geometric albedos of Ag = 0.40 \pm 0.12 across 290-450nm and Ag < 0.12 across 450-570nm at 1-sigma confidence. The albedo decrease toward longer wavelengths is also apparent when using six wavelength bins over the same wavelength range. This can be interpreted as evidence for optically thick reflective clouds on the dayside hemisphere with sodium absorption suppressing the scattered light signal beyond ~450nm. Our best-fit albedo values imply that HD189733b would appear a deep blue color at visible wavelengths.

The ubiquity of worlds beyond our Solar System confounds us.

Motivated by the reported dearth of debris discs around M stars, we use survival models to study the occurrence of planetesimal discs around them. These survival models describe a planetesimal disc with a small number of parameters, determine if it may survive a series of dynamical processes and compute the associated infrared excess. For the WISE satellite, we demonstrate that the dearth of debris discs around M stars may be attributed to the small semi-major axes generally probed if either: 1. the dust grains behave like blackbodies emitting at a peak wavelength coincident with the observed one; 2. or the grains are hotter than predicted by their blackbody temperatures and emit at peak wavelengths that are shorter than the observed one. At these small distances from the M star, planetesimals are unlikely to survive or persist for time scales of 300 Myr or longer if the disc is too massive. Conversely, our survival models allow for the existence of a large population of low-mass debris discs that are too faint to be detected with current instruments. However, our interpretation becomes less clear and large infrared excesses are allowed if only one of these scenarios holds: 3. the dust grains are hotter than blackbody and predominantly emit at the observed wavelength; 4. or are blackbody in nature and emit at peak wavelengths longer than the observed one. Both scenarios imply that the parent planetesimals reside at larger distances from the star than inferred if the dust grains behaved like blackbodies. In all scenarios, we show that the infrared excesses detected at 22 and 70 microns from AU Mic are easily reconciled with its young age. We elucidate the conditions under which stellar wind drag may be neglected when considering dust populations around M stars. The WISE satellite should be capable of detecting debris discs around young M stars with ages on the order of 10 Myr.

Supersonic flows are expected to exist in the atmospheres of irradiated exoplanets, but the question of whether shocks develop lingers. Specifically, it reduces to whether continuous flow in a closed loop may become supersonic and if some portions of the supersonic flow steepen into shocks. We first demonstrate that continuous, supersonic flow may exist in two flavors: isentropic and non-isentropic, with shocks being included in the latter class of solutions. Supersonic flow is a necessary but insufficient condition for shocks to develop. The development of a shock requires the characteristics of neighboring points in a flow to intersect. We demonstrate that the intersection of characteristics may be quantified via knowledge of the Mach number. Finally, we examine 3D simulations of hot Jovian atmospheres and demonstrate that shock formation is expected to occur mostly on the dayside hemisphere, upstream of the substellar point, because the enhanced temperatures near the substellar point provide a natural pressure barrier for the returning flow. Understanding the role of shocks in irradiated exoplanetary atmospheres is relevant to correctly modeling observables such as the peak offsets of infrared phase curves.

Characterizing atmospheres beyond the Solar System is an endeavor no longer confined to the realm of science fiction.

We investigate the stability of super Earth atmospheres around M stars using a 7-parameter, analytical framework. We construct stability diagrams in the parameter space of exoplanetary radius versus semi-major axis and elucidate the regions in which the atmospheres are stable against the condensation of their major constituents, out of the gas phase, on their permanent nightside hemispheres. We find that super Earth atmospheres which are nitrogen-dominated ("Earth-like") occupy a smaller region of allowed parameter space, compared to hydrogen-dominated atmospheres, because of the dual effects of diminished advection and enhanced radiative cooling. Furthermore, some super Earths which reside within the habitable zones of M stars may not possess stable atmospheres, depending on the mean molecular weight and infrared photospheric pressure of their atmospheres. We apply our stability diagrams to GJ 436b and GJ 1214b, and demonstrate that atmospheric compositions with high mean molecular weights are disfavoured if these exoplanets possess solid surfaces and shallow atmospheres. Finally, we construct stability diagrams tailored to the Kepler dataset, for G and K stars, and predict that about half of the exoplanet candidates are expected to habour stable atmospheres if Earth-like conditions are assumed. We include 55 Cancri e and CoRoT-7b in our stability diagram for G stars.

Hot Jupiters, due to the proximity to their parent stars, are subjected to a strong irradiating flux which governs their radiative and dynamical properties. We compute a suite of 3D circulation models with dual-band radiative transfer, exploring a relevant range of irradiation temperatures, both with and without temperature inversions. We find that, for irradiation temperatures T \lesssim 2000 K, heat redistribution is very efficient, producing comparable day- and night-side fluxes. For Tirr \approx 2200-2400 K, the redistribution starts to break down, resulting in a high day-night flux contrast. Our simulations indicate that the efficiency of redistribution is primarily governed by the ratio of advective to radiative timescales. Models with temperature inversions display a higher day-night contrast due to the deposition of starlight at higher altitudes, but we find this opacity-driven effect to be secondary compared to the effects of irradiation. The hotspot offset from the substellar point is large when insolation is weak and redistribution is efficient, and decreases as redistribution breaks down. The atmospheric flow can be potentially subjected to the Kelvin-Helmholtz instability (as indicated by the Richardson number) only in the uppermost layers, with a depth that penetrates down to pressures of a few millibars at most. Shocks penetrate deeper, down to several bars in the hottest model. Ohmic dissipation generally occurs down to deeper levels than shock dissipation (to tens of bars), but the penetration depth varies with the atmospheric opacity. The total dissipated Ohmic power increases steeply with the strength of the irradiating flux and the dissipation depth recedes into the atmosphere, favoring radius inflation in the most irradiated objects.

In astrophysical regimes where the collisional excitation of hydrogen atoms is relevant, the cross sections for the interactions of hydrogen atoms with electrons and protons are necessary for calculating line profiles and intensities. In particular, at relative velocities exceeding ~1000 km/s, collisional excitation by protons dominates over that by electrons. Surprisingly, the hydrogen-proton cross sections at these velocities do not exist for atomic levels of n >= 4, forcing researchers to utilize extrapolation via inaccurate scaling laws. In this study, we present a faster and improved algorithm for computing cross sections for the hydrogen-proton collisional system, including excitation and charge transfer to the n >= 2 levels of the hydrogen atom. We develop a code named BDSCx which directly solves the Schrodinger equation with variable (but non-adaptive) resolution and utilizes a hybrid spatial-Fourier grid. Our novel hybrid grid reduces the number of grid points needed from ~4000 n^6 (for a "brute force", Cartesian grid) to ~2000 n^4 and speeds up the computation by a factor ~50 for calculations going up to n = 4 . We present (l,m)-resolved results for charge-transfer and excitation final states for n = 2--4 and for projectile energies of 5--80 keV, as well as fitting functions for the cross sections. The ability to accurately compute proton-hydrogen cross sections to n = 4 allows us to calculate the Balmer decrement, the ratio of Balmer alpha to Balmer beta line intensities. We find that the Balmer decrement starts to increase beyond its largely constant value of 2--3 below 10 keV, reaching values of 4--5 at 5 keV, thus complicating its use as a diagnostic of dust extinction when fast (~1000$ km/s) shocks are impinging upon the ambient interstellar medium.

Using semi-analytical, one-dimensional models, we elucidate the influence of scattering and absorption on the degree of Ohmic dissipation in hot Jovian atmospheres. With the assumption of Saha equilibrium, the variation in temperature is the main driver of the variations in the electrical conductivity, induced current and Ohmic power dissipated. Atmospheres possessing temperature inversions tend to dissipate most of the Ohmic power superficially, at high altitudes, whereas those without temperature inversions are capable of greater dissipation deeper down. Scattering in the optical range of wavelengths tends to cool the lower atmosphere, thus reducing the degree of dissipation at depth. Purely absorbing cloud decks (in the infrared), of a finite extent in height, allow for localized reductions in dissipation and may reverse a temperature inversion if they are dense and thick enough, thus greatly enhancing the dissipation at depth. If Ohmic dissipation is the mechanism for inflating hot Jupiters, then variations in the atmospheric opacity (which may be interpreted as arising from variations in metallicity and cloud/haze properties) and magnetic field strength naturally produce a scatter in the measured radii at a given strength of irradiation. Future work will determine if these effects are dominant over evolutionary effects, which also contribute a scatter to the measured radii.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

We present the most sensitive ultraviolet observations of Supernova 1987A to date. Imaging spectroscopy from the Hubble Space Telescope-Cosmic Origins Spectrograph shows many narrow (dv \sim 300 km/s) emission lines from the circumstellar ring, broad (dv \sim 10 -- 20 x 10^3 km/s) emission lines from the reverse shock, and ultraviolet continuum emission. The high signal-to-noise (> 40 per resolution element) broad LyA emission is excited by soft X-ray and EUV heating of mostly neutral gas in the circumstellar ring and outer supernova debris. The ultraviolet continuum at \lambda > 1350A can be explained by HI 2-photon emission from the same region. We confirm our earlier, tentative detection of NV \lambda 1240 emission from the reverse shock and we present the first detections of broad HeII \lambda1640, CIV \lambda1550, and NIV] \lambda1486 emission lines from the reverse shock. The helium abundance in the high-velocity material is He/H = 0.14 +/- 0.06. The NV/H-alpha line ratio requires partial ion-electron equilibration (T_{e}/T_{p} \approx 0.14 - 0.35). We find that the N/C abundance ratio in the gas crossing the reverse shock is significantly higher than that in the circumstellar ring, a result that may be attributed to chemical stratification in the outer envelope of the supernova progenitor. The N/C abundance ratio may have been stratified prior to the ring expulsion, or this result may indicate continued CNO processing in the progenitor subsequent to the expulsion of the circumstellar ring.

Improving upon our purely dynamical work, we present three-dimensional simulations of the atmospheric circulation on Earth-like (exo)planets and hot Jupiters using the GFDL-Princeton Flexible Modeling System (FMS). As the first steps away from the dynamical benchmarks of Heng, Menou & Phillipps (2011), we add dual-band radiative transfer and dry convective adjustment schemes to our computational setup. Our treatment of radiative transfer assumes stellar irradiation to peak at a wavelength shorter than and distinct from that at which the exoplanet re-emits radiation ("shortwave" versus "longwave"), and also uses a two-stream approximation. Convection is mimicked by adjusting unstable lapse rates to the dry adiabat. The bottom of the atmosphere is bounded by a uniform slab with a finite thermal inertia. For our models of hot Jupiters, we include an analytical formalism for calculating temperature-pressure profiles, in radiative equilibrium, which accounts for the effect of collision-induced absorption via a single parameter. We discuss our results within the context of: the predicted temperature-pressure profiles and the absence/presence of a temperature inversion; the possible maintenance, via atmospheric circulation, of the putative high-altitude, shortwave absorber expected to produce these inversions; the angular/temporal offset of the hot spot from the substellar point, its robustness to our ignorance of hyperviscosity and hence its utility in distinguishing between different hot Jovian atmospheres; and various zonal-mean flow quantities. Our work bridges the gap between three-dimensional simulations which are purely dynamical and those which incorporate multi-band radiative transfer, thus contributing to the construction of a required hierarchy of three-dimensional theoretical models.

Motivated by the work of Guillot (2010), we present a semi-analytical formalism for calculating the temperature-pressure profiles in hot Jovian atmospheres which includes the effects of clouds/hazes and collision-induced absorption. Using the dual-band approximation, we assume that stellar irradiation and thermal emission from the hot Jupiter occur at distinct wavelengths ("shortwave" versus "longwave"). For a purely absorbing cloud/haze, we demonstrate its dual effect of cooling and warming the upper and lower atmosphere, respectively, which modifies, in a non-trivial manner, the condition for whether a temperature inversion is present in the upper atmosphere. The warming effect becomes more pronounced as the cloud/haze deck resides at greater depths. If it sits below the shortwave photosphere, the warming effect becomes either more subdued or ceases altogether. If shortwave scattering is present, its dual effect is to warm and cool the upper and lower atmosphere, respectively, thus counteracting the effects of enhanced longwave absorption by the cloud/haze. We make a tentative comparison of a 4-parameter model to the temperature-pressure data points inferred from the observations of HD 189733b and estimate that its Bond albedo is approximately 10%. Besides their utility in developing physical intuition, our semi-analytical models are a guide for the parameter space exploration of hot Jovian atmospheres via three-dimensional simulations of atmospheric circulation.

We present the late-time optical light curve of the ejecta of SN 1987A measured from HST imaging observations spanning the past 17 years. We find that the flux from the ejecta declined up to around year 2001, powered by the radioactive decay of 44Ti. Then the flux started to increase, more than doubling by the end of 2009. We show that the increase is the result of energy deposited by X-rays produced in the interaction with the circumstellar medium. We suggest that the change of the dominant energy input to the ejecta, from internal to external, marks the transition from supernova to supernova remnant. The details of the observations and the modelling are described in the accompanying supplementary information.

We present a method to estimate the mass of the debris disc in the HD 69830 system, which also hosts three exoplanets with Neptune-like minimum masses. By considering the range of published stellar ages, we interpret the infrared emission from the debris disc as originating from a steady state, collisional cascade of dust grains. Using dynamical survival models subjected to observational constraints, we estimate the allowed range of disc masses. If the disc has an age of about 1 Gyr, then its mass is several times that of our asteroid belt. The maximum allowed age for the disc and the number of planetesimals it contains are determined by the assumed value for the binding energy of the planetesimals. If one insists on interpreting the disc as being transient, then this mass estimate becomes an upper limit.

Employing realistic and consistent atmosphere boundary conditions, we have generated evolutionary models for brown dwarfs and very-low-mass stars (VLMs) for different metallicities ([Fe/H]), with and without clouds. We find that the spread in radius at a given mass and age can be as large as $\sim$10% to $\sim$25%, with higher-metallicity, higher-cloud-thickness atmospheres resulting quite naturally in larger radii. For each 0.1 dex increase in [Fe/H], radii increase by $\sim$1% to $\sim$2.5%, depending upon age and mass. We also find that, while for smaller masses and older ages brown dwarf radii decrease with increasing helium fraction ($Y$) (as expected), for more massive brown dwarfs and a wide range of ages they increase with helium fraction. The increase in radius in going from $Y=0.25$ to $Y=0.28$ can be as large as $\sim$0.025 \rj\ ($\sim$2.5%). Furthermore, we find that for VLMs an increase in atmospheric metallicity from 0.0 to 0.5 dex, increases radii by $\sim$4%, and from -0.5 to 0.5 dex by $\sim$10%. Therefore, we suggest that opacity due to higher metallicity might naturally account for the apparent radius anomalies in some eclipsing VLM systems. Ten to twenty-five percent variations in radius exceed errors stemming from uncertainities in the equation of state alone. This serves to emphasize that transit and eclipse measurements of brown dwarf radii constrain numerous effects collectively, importantly including the atmosphere and condensate cloud models, and not just the equation of state. At all times, one is testing a multi-parameter theory, and not a universal radius$-$mass relation.

We use three-dimensional simulations to study the atmospheric circulation on the first Earth-sized exoplanet discovered in the habitable zone of an M star. We treat Gliese 581g as a scaled-up version of Earth by considering increased values for the exoplanetary radius and surface gravity, while retaining terrestrial values for parameters which are unconstrained by current observations. We examine the long-term, global temperature and wind maps near the surface of the exoplanet --- the climate. The specific locations for habitability on Gliese 581g depend on whether the exoplanet is tidally-locked and how fast radiative cooling occurs on a global scale. Independent of whether the existence of Gliese 581g is confirmed, our study highlights the use of general circulation models to quantify the atmospheric circulation on potentially habitable, Earth-sized exoplanets, which will be the prime targets of exoplanet discovery and characterization campaigns in the next decade.

The complexity of atmospheric modelling and its inherent non-linearity, together with the limited amount of data of exoplanets available, motivate model intercomparisons and benchmark tests. In the geophysical community, the Held-Suarez test is a standard benchmark for comparing dynamical core simulations of the Earth's atmosphere with different solvers, based on statistically-averaged flow quantities. In the present study, we perform analogues of the Held-Suarez test for tidally-locked exoplanets with the GFDL-Princeton Flexible Modeling System (FMS) by subjecting both the spectral and finite difference dynamical cores to a suite of tests, including the standard benchmark for Earth, a hypothetical tidally-locked Earth, a "shallow" hot Jupiter model and a "deep" model of HD 209458b. We find qualitative and quantitative agreement between the solvers for the Earth, tidally-locked Earth and shallow hot Jupiter benchmarks, but the agreement is less than satisfactory for the deep model of HD 209458b. Further investigation reveals that closer agreement may be attained by arbitrarily adjusting the values of the horizontal dissipation parameters in the two solvers, but it remains the case that the magnitude of the horizontal dissipation is not easily specified from first principles. Irrespective of radiative transfer or chemical composition considerations, our study points to limitations in our ability to accurately model hot Jupiter atmospheres with meteorological solvers at the level of ten percent for the temperature field and several tens of percent for the velocity field. Direct wind measurements should thus be particularly constraining for the models. Our suite of benchmark tests also provides a reference point for researchers wishing to adapt their codes to study the atmospheric circulation regimes of tidally-locked Earths/Neptunes/Jupiters.

Observations with the Hubble Space Telescope (HST), conducted since 1990, now offer an unprecedented glimpse into fast astrophysical shocks in the young remnant of supernova 1987A. Comparing observations taken in 2010 using the refurbished instruments on HST with data taken in 2004, just before the Space Telescope Imaging Spectrograph failed, we find that the Ly-a and H-a lines from shock emission continue to brighten, while their maximum velocities continue to decrease. We observe broad blueshifted Ly-a, which we attribute to resonant scattering of photons emitted from hotspots on the equatorial ring. We also detect NV~\lambda\lambda 1239,1243 A line emission, but only to the red of Ly-A. The profiles of the NV lines differ markedly from that of H-a, suggesting that the N^{4+} ions are scattered and accelerated by turbulent electromagnetic fields that isotropize the ions in the collisionless shock.

Turbulent, two-dimensional, hydrodynamic flows are characterized by the emergence of coherent, long-lived vortices without a need to invoke special initial conditions. Vortices have the ability to sequester particles, with typical radii from about 1 mm to 10 cm, that are slightly decoupled from the gas. A generic feature of discs with surface density and effective temperature profiles that are decreasing, power-law functions of radial distance is that four vortex zones exist for a fixed particle size. In particular, two of the zones form an annulus at intermediate radial distances within which small particles reside. Particle capture by vortices occurs on a dynamical time scale near and at the boundaries of this annulus. As the disc ages and the particles grow via coagulation, the size of the annulus shrinks. Older discs prefer to capture smaller particles because the gas surface density decreases with time, a phenomenon we term "vortex aging". More viscous, more dust-opaque and/or less massive discs can have vortices that age faster and trap a broader range of particle sizes throughout the lifetime of the disc. Thus, how efficiently a disc retains its mass in solids depends on the relative time scales between coagulation and vortex aging. If vortices form in protoplanetary discs, they are important in discs with typical masses and for particles that are likely to condense out of the protostellar nebula. Particle capture also occurs at distances relevant to planet formation. Future infrared, submillimetre and centimetre observations of grain opacity as a function of radial distance will test the hypothesis that vortices serve as nurseries for particle growth in protoplanetary discs.

We investigate the survival of planetesimal discs over Gyr timescales, using a unified approach that is applicable to all Keplerian discs of solid bodies -- dust grains, asteroids, planets, etc. Planetesimal discs can be characterized locally by four parameters: surface density, semi-major axis, planetesimal size and planetesimal radial velocity dispersion. Any planetesimal disc must have survived all dynamical processes, including gravitational instability, dynamical chaos, gravitational scattering, physical collisions, and radiation forces, that would lead to significant evolution over its lifetime. These processes lead to a rich set of constraints that strongly restrict the possible properties of long-lived discs. Within this framework, we also discuss the detection of planetesimal discs using radial velocity measurements, transits, microlensing, and the infrared emission from the planetesimals themselves or from dust generated by planetesimal collisions.

Motivated by debris disk studies, we investigate the gravitational microlensing of background starlight by a planetesimal disk around a foreground star. We use dynamical survival models to construct a plausible example of a planetesimal disk and study its microlensing properties using established ideas of microlensing by small bodies. When a solar-type source star passes behind a planetesimal disk, the microlensing light curve may exhibit short-term, low-amplitude residuals caused by planetesimals several orders of magnitude below Earth mass. The minimum planetesimal mass probed depends on the photometric sensitivity and the size of the source star, and is lower when the planetesimal lens is located closer to us. Planetesimal lenses may be found more nearby than stellar lenses because the steepness of the planetesimal mass distribution changes how the microlensing signal depends on the lens/source distance ratio. Microlensing searches for planetesimals require essentially continuous monitoring programs that are already feasible and can potentially set constraints on models of debris disks, the supposed extrasolar analogues of Kuiper belts.

A concise and critical review of Balmer-dominated shocks (BDSs) is presented, summarizing the state of theory and observations, including models with/without shock precursors and their synergy with atomic physics. Observations of BDSs in supernova remnants are reviewed on an object-by-object basis. The relevance of BDSs towards understanding the acceleration of cosmic rays in shocks is emphasized. Probable and possible detections of BDSs in astrophysical objects other than supernova remnants, including pulsar wind nebulae and high-redshift galaxies, are described. The case for the continued future of studying BDSs in astrophysics is made, including their relevance towards understanding electron-ion temperature equilibration in collisionless shocks.

We present a linear analysis of inviscid, incompressible, magnetohydrodynamic (MHD) shallow water systems. In spherical geometry, a generic property of such systems is the existence of five wave modes. Three of them (two magneto-Poincare modes and one magneto-Rossby mode) are previously known. The other two wave modes are strongly influenced by the magnetic field and rotation, and have substantially lower angular frequencies; as such, we term them "magnetostrophic modes". We obtain analytical functions for the velocity, height and magnetic field perturbations in the limit that the magnitude of the MHD analogue of Lamb's parameter is large. On a sphere, the magnetostrophic modes reside near the poles, while the other modes are equatorially confined. Magnetostrophic modes may be an ingredient in explaining the frequency drifts observed in Type I X-ray bursts from neutron stars.

We present theoretical calculations of the X-ray scattering properties of porous grain aggregates with olivine monomers. The small and large angle scattering properties of these aggregates are governed by the global structure and substructure of the grain, respectively. We construct two diagnostics, R_X and T_X, based on the optical and X-ray properties of the aggregates, and apply them to a Chandra measurement of the dust halo around the Galactic binary GX13+1. Grain aggregates with porosities higher than about 0.55 are ruled out. Future high-precision observations of X-ray dust haloes together with detailed modeling of the X-ray scattering properties of porous grain mixtures will further constrain the presence of porous grain aggregates in a given dust population.

In order to attract and retain excellent researchers and diverse individuals in astrophysics, we recommend action be taken in several key areas impacting young scientists: (1) Maintain balance between large collaborations and individual projects through distribution of funding; encourage public releases of observational and simulation data for use by a broader community. (2) Improve the involvement of women, particularly at leading institutions. (3) Address the critical shortage of child care options and design reasonable profession-wide parental leave policies. (4) Streamline the job application and hiring process. We summarize our reasons for bringing these areas to the attention of the committee, and we suggest several practical steps that can be taken to address them.

Balmer-dominated shocks in supernova remnants (SNRs) produce strong hydrogen lines with a two-component profile composed of a narrow contribution from cold upstream hydrogen atoms, and a broad contribution from hydrogen atoms that have undergone charge transfer reactions with hot protons. Observations of emission lines from edge-wise shocks in SNRs can constrain the gas velocity and collisionless electron heating at the shock front. Downstream hydrogen atoms engage in charge transfer, excitation and ionization reactions, defining an interaction region called the shock transition zone. The properties of hot hydrogen atoms produced by charge transfers (called broad neutrals) are critical for accurately calculating the structure and radiation from the shock transition zone. This paper is the third in a series describing the kinetic, fluid and emission properties of Balmer-dominated shocks, and is the first to properly treat the effect of broad neutral kinetics on shock transition zone structure. We use our models to extract shock parameters from observations of Balmer-dominated SNRs. We find that inferred shock velocities and electron temperatures are lower than those of previous calculations by <10% for v_s<1500 km/s, and by 10-30% for v_s>1500 km/s. This effect is primarily due to the fact that excitation by proton collisions and charge transfer to excited levels favor the high speed part of the neutral hydrogen velocity distribution. Our results have a strong dependence on the ratio of electron to proton temperatures, \beta=T_e/T_p, which allows us to construct a relation \beta(v_s) between the temperature ratio and shock velocity. We compare our calculations to previous results by Ghavamian et al. (2007).

We present a proof-of-concept study that dust extinction curves can be extracted from the infrared (IR), optical, ultraviolet (UV) and X-ray afterglow observations of GRBs without assuming known extinction laws. We focus on GRB 050525A (z = 0.606), for which we also present IR observations from the Spitzer Space Telescope at about 2.3 days post-burst. We construct the spectral energy distribution (SED) of the afterglow and use it to derive the dust extinction curve of the host galaxy in 7 optical/UV wavebands. By comparing our derived extinction curve to known templates, we see that the Galactic or Milky Way extinction laws are disfavored versus the Small and Large Magellanic Cloud (SMC and LMC) ones, but that we cannot rule out the presence of a LMC-like 2175 angstrom bump in our extinction curve. The dust-to-gas ratio present within the host galaxy of GRB 050525A is similar to that found in the LMC, while about 10 to 40% more dust is required if the SMC template is assumed. Our method is useful to observatories that are capable of simultaneously observing GRB afterglows in multiple wavebands from the IR to the X-ray.

Charge transfer (or exchange) reactions between hydrogen atoms and protons in collisionless shocks of supernova remnants (SNRs) are a natural way of producing broad Balmer, Lyman and other lines of hydrogen. We wish to quantify the importance of shock-induced, non-thermal hydrogen emission from SNRs in young galaxies. We present a method to estimate the luminosity of broad ($\sim 1000$ km s$^{-1}$) Ly$\alpha$, Ly$\beta$, Ly$\gamma$, H$\beta$ and P$\alpha$ lines, as well as the broad and narrow luminosities of the two-photon (2$\gamma$) continuum, from existing measurements of the H$\alpha$ flux. The expected broad Ly$\alpha$ luminosity per object is at most $\sim 10^{36}$ erg s$^{-1}$. In principle, broad, ``non-radiative'' Ly$\alpha$ from SNRs in young galaxies can be directly observed in the optical range of wavelengths. However, by taking into consideration the different rates between core collapse and thermonuclear supernovae, as well as the duration we expect to observe such Ly$\alpha$ emission from SNRs, we expect their contribution to the total Ly$\alpha$ luminosity from $z \sim 3$ to 5 galaxies to be negligibly small ($\sim 0.001 %$), compared to the radiative shock mechanism described by Shull & Silk (1979). Though broad, non-thermal Ly$\alpha$ emission has never been observed, these photons are produced in SNRs and hence the non-radiative Ly$\alpha$ luminosity is a part of the intrinsic Ly$\alpha$ spectrum of young galaxies.

We report on the latest (2007 Jan) observations of supernova remnant (SNR) 1987A from the XMM-Newton mission. Since the 2003 May observations of Haberl et al. (2006), 11 emission lines have experienced increases in flux by factors ~ 3 to 10, with the 775 eV line of O VIII showing the greatest increase; we have observed 6 lines of Fe XVII and Fe XVIII previously unreported by XMM-Newton. A two-shock model representing plasmas in non-equilibrium ionization is fitted to the EPIC-pn spectra, yielding temperatures of ~ 0.4 and ~ 3 keV, as well as elemental abundances for N, O, Ne, Mg, Si, S and Fe. We demonstrate that the abundance ratio of N and O can be constrained to less than ~20% accuracy. Within the same confidence interval, the same analysis suggests that the C+N+O abundance varies from ~ 1.1 to 1.4 X 10^-4. Normalizing our obtained abundances by the Large Magellanic Cloud (LMC) values of Hughes, Hayashi & Koyama (1998), we find that O, Ne, Mg and Fe are under-abundant, while Si and S are over-abundant, consistent with the findings of Aschenbach (2007). Such a result has implications for both the single-star and binary accretion/merger models for the progenitor of SNR 1987A. In the context of the binary merger scenario proposed by Morris & Podsiadlowski (2006, 2007), material forming the inner, equatorial ring was expelled after the merger, implying that either our derived Fe abundance is inconsistent with typical LMC values or that iron is under-abundant at the site of the progenitor star of SNR 1987A.

We examine the structure of the post-shock region in supernova remnants (SNRs). The ``shock transition zone'' is set up by charge transfer and ionization events between atoms and ions, and has a width $\sim 10^{15}$ cm$^{-2}$ $n^{-1}_0$, where $n_0$ is the total pre-shock density (including both atoms and ions). For Balmer-dominated SNRs with shock velocity $v_s \gtrsim 1000$ km s$^{-1}$, the Rankine-Hugoniot conditions for ion velocity and temperature are obeyed instantly, leaving the full width at half-maximum (FWHM) of the broad H$\alpha$ line versus $v_s$ relation intact. However, the spatial variation in the post-shock densities is relevant to the problem of Ly$\alpha$ resonant scattering in young, core-collapse SNRs. Both two- (pre-shock atoms and ions) and three-component (pre-shock atoms, broad neutrals and ions) models are considered. We compute the spatial emissivities of the broad ($\xi_b$) and narrow ($\xi_n$) H$\alpha$ lines; a calculation of these emissivities in SN 1006 is in general agreement with the computed ones of Raymond et al. (2007). The (dimensionless) spatial shift, $\Theta_{\rm{shift}}$, between the centroids of $\xi_b$ and $\xi_n$ is unique for a given shock velocity and $f_{\rm{ion}}$, the pre-shock ion fraction. Measurements of $\Theta_{\rm{shift}}$ can be used to constrain $n_0$.

The reverse shock of supernova remnant (SNR) 1987A emits in H-alpha and Ly-alpha, and comes in two flavors: surface and interior. The former is due to direct, impact excitation of hydrogen atoms crossing the shock, while the latter is the result of charge transfer reactions between these atoms and slower, post-shock ions. Interior and surface emission are analogous to the broad- and narrow-line components observed in Balmer-dominated SNRs. I summarize a formalism to derive line intensities and ratios in these SNRs, as well as a study of the transition zone in supernova shocks; I include an appendix where I derive in detail the ratio of broad to narrow H-alpha emission. Further study of the reverse shock emission from SNR 1987A will allow us to predict when it will vanish and further investigate the origins of the interior emission.

Long gamma-ray bursts (GRBs) are likely associated with the collapse of massive stars, which produce dust and are born in dusty environments. Absorption and scattering of ultraviolet/X-ray photons from the prompt, optical flash and afterglow emission of the GRB produce dust echoes. We perform time-dependent calculations of these echoes, accounting for the evolution of the dust grain distribution due to selective grain destruction by the GRB radiation, and for off-axis beaming. We explore cloud configurations of differing density and size -- the echo light curve and spectrum depend on the cloud radius, with larger clouds peaking at longer wavelengths. For a region ~ 3 pc in size with hydrogen density ~ 1000 per cubic centimeter, the echo spectrum peaks at ~ 3.6 microns and ~ 8.8 eV for thermal and scattered components, respectively. Dust echoes should be detectable with the Very Large Telescope up to z ~ 0.1, IRAC onboard the Spitzer Space Telescope up to z ~ 0.2, and NICMOS onboard the Hubble Space Telescope up to z ~ 0.3. Furthermore, the shape of the echo light curve allows one to infer: the jet opening angle; the inclination of the jet axis with respect to the line of sight; the size of the dust-emitting region. For sources with symmetric, bipolar jets, dust echoes exhibit two bumps in the light curve, making them easily distinguishable from the rebrightening due to an underlying supernova.

We present a new formalism to describe the ratios and profiles of emission lines from hydrogen in Balmer-dominated shocks. We use this model to interpret the measured widths and ratios of broad and narrow H-alpha, H-beta and Ly-alpha emission lines in supernova remnants (SNRs). Our model results agree fairly well with those obtained previously by Chevalier, Kirshner & Raymond (1980) and are consistent with observations of several SNRs. The same model fails to account for the ratio of broad to narrow line emission from the reverse shock in SNR 1987A as observed by Heng et al. (2006). We suggest that this discrepancy between theory and observation results from a faulty assumption that Balmer-dominated shocks can be treated as sharp discontinuities. If the spatial structure of the shock transition zone is taken into account, the predicted ratios of broad to narrow line emission in most SNRs will change by modest factors, but the ratio in SNR 1987A will increase substantially. Significantly greater shock velocities will be required to account for the observed full widths at half-maximum of the broad emission lines in most SNRs.

We present new (2004 July) G750L and G140L Space Telescope Imaging Spectrograph (STIS) data of the H-alpha and Ly-alpha emission from supernova remnant (SNR) 1987A. With the aid of earlier data, from Oct 1997 to Oct 2002, we track the local evolution of Ly-alpha emission and both the local and global evolution of H-alpha emission. In addition to emission which we can clearly attribute to the surface of the reverse shock, we also measure comparable emission, in both H-alpha and Ly-alpha, which appears to emerge from supernova debris interior to the surface. New observations taken through slits positioned slightly eastward and westward of a central slit show a departure from cylindrical symmetry in the H-alpha surface emission. Using a combination of old and new observations, we construct a light curve of the total H-alpha flux, F, from the reverse shock, which has increased by a factor ~ 4 over about 8 years. However, due to large systematic uncertainties, we are unable to discern between the two limiting behaviours of the flux - F ~ t (self-similar expansion) and F ~ t^5 (halting of the reverse shock). Such a determination is relevant to the question of whether the reverse shock emission will vanish in less than about 7 years (Smith et al. 2005). Future deep, low- or moderate-resolution spectra are essential for accomplishing this task.

We use low-dispersion spectra obtained at the Magellan Observatory to study the broad H-alpha emission from the reverse shock of the infant supernova remnant SNR1987A. These spectra demonstrate that the spatio-kinematic structure of the reverse shock can be distinguished from that of the circumstellar ring and hotspots, even at ground-based spatial resolution. We measure a total dereddened H-alpha flux of 1.99(pm0.22)e-13 ergs/s/cm2 at an epoch 18.00 years after outburst. At 50 kpc, the total reverse shock luminosity in H-alpha is roughly 15 Lsun, which implies a total flux of neutral hydrogen atoms across the reverse shock of 8.9e46/s, or roughly 2.3e-3 Msun/yr. This represents an increase by a factor ~4 since 1997. Lyman continuum radiation from gas shocked by the forward blast wave can ionize neutral hydrogen atoms in the supernova debris before they reach the reverse shock. If the inward flux of ionizing photons exceeds the flux of hydrogen atoms approaching the reverse shock, this pre-ionization will shut off the broad Ly-alpha and H-alpha emission. The observed X-ray emission of SNR1987A implies that the ratio of ionizing flux to hydrogen atom flux across the reverse shock is presently at least 0.04. The X-ray emission is increasing much faster than the flux of atoms, and if these trends continue, we estimate that the broad Ly-alpha and H-alpha emission will vanish within about 7 years.