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

We present a method for detecting starspots on cool stars using the cross-correlation function (CCF) of high resolution molecular spectral templates applied to archival high-resolution spectra of G and K stars observed with HARPS/HARPS-N. We report non-detections of starspots on the Sun even when the Sun was spotted, the solar twin 18 Scorpii, and the very spotted Sun-like star HAT-P-11, suggesting that Sun-like starspot distributions will be invisible to the CCF technique, and should not produce molecular absorption signals which might be confused for signatures of exoplanet atmospheres. We detect strong TiO absorption in the T Tauri K-dwarfs LkCa 4 and AA Tau, consistent with significant coverage by cool regions. We show that despite the non-detections, the technique is sensitive to relatively small spot coverages on M dwarfs and large starspot areas on Sun-like stars.

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.

We assess broadband color filters for the two fast cameras on the PLAnetary Transits and Oscillations (PLATO) of stars space mission with respect to exoplanetary atmospheric characterization. We focus on Ultra Hot Jupiters and Hot Jupiters placed 25pc and 100pc away from the Earth and low mass low density planets placed 10pc and 25pc away. Our analysis takes as input literature values for the difference in transit depth between the broadband lower (500 to 675nm) wavelength interval (hereafter referred to as blue) and the upper (675-1125nm) broadband wavelength interval (hereafter referred to as red) for transmission, occultation and phase curve analyses. Planets orbiting main sequence central stars with stellar classes F, G, K and M are investigated. We calculate the signal-to-noise ratio with respect to photon and instrument noise for detecting the difference in transit depth between the two spectral intervals. Results suggest that bulk atmospheric composition and planetary geometric albedos could be detected for (Ultra) Hot Jupiters up to about 100pc (about 25pc) with strong (moderate) Rayleigh extinction. Phase curve information could be extracted for Ultra Hot Jupiters orbiting K and G dwarf stars up to 25pc away. For low mass low density planets, basic atmospheric types (primary and water-dominated) and the presence of sub-micron hazes in the upper atmosphere could be distinguished for up to a handful of cases up to about 10pc.

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.

A terrestrial planet is molten during formation and may remain so if subject to intense insolation or tidal forces. Observations continue to favour the detection and characterisation of hot planets, potentially with large outgassed atmospheres. We aim to determine the radius of hot Earth-like planets with large outgassed atmospheres and explore differences between molten and solid silicate planets and their influence on the mass-radius relationship and transmission and emission spectra. An interior-atmosphere model, combined with static structure calculations, tracks the evolving radius of a rocky mantle that is outgassing CO$_2$ and H$_2$O. Synthetic emission and transmission spectra are generated for CO$_2$ and H$_2$O dominated atmospheres. Atmospheres dominated by CO$_2$ suppress the outgassing of H$_2$O to a greater extent than previously realised, as previous studies have applied an erroneous relationship between volatile mass and partial pressure. We therefore predict more H$_2$O can be retained by the interior during the later stages of magma ocean crystallisation. Furthermore, formation of a lid at the surface can tie outgassing of H$_2$O to the efficiency of heat transport through the lid, rather than the atmosphere's radiative timescale. Contraction of the mantle as it solidifies gives $\sim5\%$ radius decrease, which can partly be offset by addition of a relatively light species to the atmosphere. We conclude that a molten silicate mantle can increase the radius of a terrestrial planet by around $5\%$ compared to its solid counterpart, or equivalently account for a $13\%$ decrease in bulk density. An outgassing atmosphere can perturb the total radius according to its speciation. Atmospheres of terrestrial planets around M-stars that are dominated by CO$_2$ or H$_2$O can be distinguished by observing facilities with extended wavelength coverage (e.g., JWST).

In this work, we investigate the potassium excess absorption around 7699A of the exoplanets HD189733b and HD209458b. For this purpose, we used high spectral resolution transit observations acquired with the 2 x 8.4m Large Binocular Telescope (LBT) and the Potsdam Echelle Polarimetric and Spectroscopic Instrument (PEPSI). For a bandwidth of 0.8A, we present a detection > 7-sigma with an absorption level of 0.18% for HD189733b. Applying the same analysis to HD209458b, we can set 3-sigma upper limit of 0.09%, even though we do not detect a K- excess absorption. The investigation suggests that the K- feature is less present in the atmosphere of HD209458b than in the one of HD189733b. This comparison confirms previous claims that the atmospheres of these two planets must have fundamentally different properties.

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.

Exoplanet science is one of the most thriving fields of modern astrophysics. A major goal is the atmospheric characterization of dozens of small, terrestrial exoplanets in order to search for signatures in their atmospheres that indicate biological activity, assess their ability to provide conditions for life as we know it, and investigate their expected atmospheric diversity. None of the currently adopted projects or missions, from ground or in space, can address these goals. In this White Paper we argue that a large space-based mission designed to detect and investigate thermal emission spectra of terrestrial exoplanets in the MIR wavelength range provides unique scientific potential to address these goals and surpasses the capabilities of other approaches. While NASA might be focusing on large missions that aim to detect terrestrial planets in reflected light, ESA has the opportunity to take leadership and spearhead the development of a large MIR exoplanet mission within the scope of the "Voyage 2050" long-term plan establishing Europe at the forefront of exoplanet science for decades to come. Given the ambitious science goals of such a mission, additional international partners might be interested in participating and contributing to a roadmap that, in the long run, leads to a successful implementation. A new, dedicated development program funded by ESA to help reduce development and implementation cost and further push some of the required key technologies would be a first important step in this direction. Ultimately, a large MIR exoplanet imaging mission will be needed to help answer one of mankind's most fundamental questions: "How unique is our Earth?"

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.

KELT-9b is a recently discovered exoplanet with a 1.49 d orbit around a B9.5/A0-type star. The unparalleled levels of UV irradiation it receives from its host star put KELT-9b in its own unique class of ultra-hot Jupiters, with an equilibrium temperature > 4000 K. The high quantities of dissociated hydrogen and atomic metals present in the dayside atmosphere of KELT-9b bear more resemblance to a K-type star than a gas giant. We present a single observation of KELT-9b during its secondary eclipse, taken with the Wide Field Camera on the Isaac Newton Telescope (INT). This observation was taken in the U-band, a window particularly sensitive to Rayleigh scattering. We do not detect a secondary eclipse signal, but our 3$\sigma$ upper limit of 181 ppm on the depth allows us to constrain the dayside temperature of KELT-9b at pressures of ~30 mbar to 4995 K (3$\sigma$). Although we can place an observational constraint of $A_g<$ 0.14, our models suggest that the actual value is considerably lower than this due to H$^-$ opacity. This places KELT-9b squarely in the albedo regime populated by its cooler cousins, almost all of which reflect very small components of the light incident on their daysides. This work demonstrates the ability of ground-based 2m-class telescopes like the INT to perform secondary eclipse studies in the NUV, which have previously only been conducted from space-based facilities.

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 long-term goals of exoplanet science is the (atmospheric) characterization of a large sample (>100) of terrestrial planets to assess their potential habitability and overall diversity. Hence, it is crucial to quantitatively evaluate and compare the scientific return of various mission concepts. Here we discuss the exoplanet yield of a space-based mid-infrared (MIR) nulling interferometer. We use Monte-Carlo simulations, based on the observed planet population statistics from the Kepler mission, to quantify the number and properties of detectable exoplanets (incl. potentially habitable planets) and we compare the results to those for a large aperture optical/NIR space telescope. We investigate how changes in the underlying technical assumptions (sensitivity and spatial resolution) impact the results and discuss scientific aspects that influence the choice for the wavelength coverage and spectral resolution. Finally, we discuss the advantages of detecting exoplanets at MIR wavelengths, summarize the current status of some key technologies, and describe what is needed in terms of further technology development to pave the road for a space-based MIR nulling interferometer for exoplanet science.

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.

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.

For the calculation of complex neutral/ionized gas phase chemical equilibria, we present a semi-analytical versatile and efficient computer program, called FastChem. The applied method is based on the solution of a system of coupled nonlinear (and linear) algebraic equations, namely the law of mass action and the element conservation equations including charge balance, in many variables. Specifically, the system of equations is decomposed into a set of coupled nonlinear equations in one variable each, which are solved analytically whenever feasible to reduce computation time. Notably, the electron density is determined by using the method of Nelder and Mead at low temperatures. The program is written in object-oriented C++ which makes it easy to couple the code with other programs, although a stand-alone version is provided. FastChem can be used in parallel or sequentially and is available under the GNU General Public License version 3 at https://github.com/exoclime/FastChem together with several sample applications. The code has been successfully validated against previous studies and its convergence behavior has been tested even for extreme physical parameter ranges down to 100 K and up to 1000 bar. FastChem converges stable and robust in even most demanding chemical situations, which posed sometimes extreme challenges for previous algorithms.

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

The benchmark exoplanet GJ 1214b is one of the best studied transiting planets in the transition zone between rocky Earth-sized planets and gas or ice giants. This class of super-Earth/mini-Neptune planets is unknown in our Solar System, yet is one of the most frequently detected classes of exoplanets. Understanding the transition from rocky to gaseous planets is a crucial step in the exploration of extrasolar planetary systems, in particular with regard to the potential habitability of this class of planets. GJ 1214b has already been studied in detail from various platforms at many different wavelengths. Our airborne observations with SOFIA add information in the Paschen-alpha cont. 1.9 micron infrared wavelength band, which is not accessible by any other current ground- or space-based instrument due to telluric absorption or limited spectral coverage. We used FLIPO and FPI+ on SOFIA to comprehensively analyse the transmission signal of the possible water-world GJ 1214b through photometric observations during transit in three optical and one infrared channels. We present four simultaneous light curves and corresponding transit depths in three optical and one infrared channel, which we compare to previous observations and state-of-the-art synthetic atmospheric models of GJ 1214b. The final precision in transit depth is between 1.5 and 2.5 times the theoretical photon noise limit, not sensitive enough to constrain the theoretical models any better than previous observations. This is the first exoplanet observation with SOFIA that uses its full set of instruments available to exoplanet spectrophotometry. Therefore we use these results to evaluate SOFIA's potential in this field and suggest future improvements.

Clouds have a strong impact on the climate of planetary atmospheres. The potential scattering greenhouse effect of CO2 ice clouds in the atmospheres of terrestrial extrasolar planets is of particular interest because it might influence the position and thus the extension of the outer boundary of the classic habitable zone around main sequence stars. Here, the impact of CO2 ice clouds on the surface temperatures of terrestrial planets with CO2 dominated atmospheres, orbiting different types of stars is studied. Additionally, their corresponding effect on the position of the outer habitable zone boundary is evaluated. For this study, a radiative-convective atmospheric model is used the calculate the surface temperatures influenced by CO2 ice particles. The clouds are included using a parametrised cloud model. The atmospheric model includes a general discrete ordinate radiative transfer that can describe the anisotropic scattering by the cloud particles accurately. A net scattering greenhouse effect caused by CO2 clouds is only obtained in a rather limited parameter range which also strongly depends on the stellar effective temperature. For cool M-stars, CO2 clouds only provide about 6 K of additional greenhouse heating in the best case scenario. On the other hand, the surface temperature for a planet around an F-type star can be increased by 30 K if carbon dioxide clouds are present. Accordingly, the extension of the habitable zone due to clouds is quite small for late-type stars. Higher stellar effective temperatures, on the other hand, can lead to outer HZ boundaries about 0.5 au farther out than the corresponding clear-sky values.

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]

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.

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

The discontinuous Galerkin finite element method (DG-FEM) is successfully applied to treat a broad variety of transport problems numerically. In this work, we use the full capacity of the DG-FEM to solve the radiative transfer equation in spherical symmetry. We present a discontinuous Galerkin method to directly solve the spherically-symmetric radiative transfer equation as a two-dimensional problem. The transport equation in spherical atmospheres is more complicated than in the plane-parallel case due to the appearance of an additional derivative with respect to the polar angle. The DG-FEM formalism allows for the exact integration of arbitrarily complex scattering phase functions, independent of the angular mesh resolution. We show that the discontinuous Galerkin method is able to describe accurately the radiative transfer in extended atmospheres and to capture discontinuities or complex scattering behaviour which might be present in the solution of certain radiative transfer tasks and can, therefore, cause severe numerical problems for other radiative transfer solution methods.

The habitable zone (HZ) describes the range of orbital distances around a star where the existence of liquid water on the surface of an Earth-like planet is in principle possible. While 3D climate studies can calculate the water vapor, ice albedo, and cloud feedback self-consistently and therefore allow for a deeper understanding and the identification of relevant climate processes, 1D model studies rely on fewer model assumptions and can be more easily applied to the large parameter space possible for exoplanets. We evaluate the applicability of 1D climate models to estimate the potential habitability of Earth-like exoplanets by comparing our 1D model results to those of 3D climate studies in the literature. We applied a cloud-free 1D radiative-convective climate model to calculate the climate of Earth-like planets around different types of main-sequence stars with varying surface albedo and relative humidity profile. These parameters depend on climate feedbacks that are not treated self-consistently in most 1D models. We compared the results to those of 3D model calculations in the literature and investigated to what extent the 1D model can approximate the surface temperatures calculated by the 3D models. The 1D parameter study results in a large range of climates possible for an Earth-sized planet with an Earth-like atmosphere and water reservoir at a certain stellar insolation. At some stellar insolations the full spectrum of climate states could be realized, i.e., uninhabitable conditions as well as habitable surface conditions, depending only on the relative humidity and surface albedo assumed. When treating the surface albedo and the relative humidity profile as parameters in 1D model studies and using the habitability constraints found by recent 3D modeling studies, the same conclusions about the potential habitability of a planet can be drawn as from 3D model calculations.

Carbon dioxide ice clouds are thought to play an important role for cold terrestrial planets with thick CO2 dominated atmospheres. Various previous studies showed that a scattering greenhouse effect by carbon dioxide ice clouds could result in a massive warming of the planetary surface. However, all of these studies only employed simplified two-stream radiative transfer schemes to describe the anisotropic scattering. Using accurate radiative transfer models with a general discrete ordinate method, this study revisits this important effect and shows that the positive climatic impact of carbon dioxide clouds was strongly overestimated in the past. The revised scattering greenhouse effect can have important implications for the early Mars, but also for planets like the early Earth or the position of the outer boundary of the habitable zone.

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.

The potential habitability of a terrestrial planet is usually defined by the possible existence of liquid water on its surface. The potential presence of liquid water depends on many factors such as, most importantly, surface temperatures. The properties of the planetary atmosphere and its interaction with the radiative energy provided by the planet's host star are thereby of decisive importance. In this study we investigate the influence of different main-sequence stars upon the climate of Earth-like extrasolar planets and their potential habitability by applying a 3D Earth climate model accounting for local and dynamical processes. The calculations have been performed for planets with Earth-like atmospheres at orbital distances where the total amount of energy received from the various host stars equals the solar constant. In contrast to previous 3D modeling studies, we include the effect of ozone radiative heating upon the vertical temperature structure of the atmospheres. The global orbital mean results obtained have been compared to those of a 1D radiative convective climate model. The different stellar spectral energy distributions lead to different surface temperatures and due to ozone heating to very different vertical temperature structures. As previous 1D studies we find higher surface temperatures for the Earth-like planet around the K-type star, and lower temperatures for the planet around the F-type star compared to an Earth-like planet around the Sun. However, this effect is more pronounced in the 3D model results than in the 1D model because the 3D model accounts for feedback processes such as the ice-albedo and the water vapor feedback. Whether the 1D model may approximate the global mean of the 3D model results strongly depends on the choice of the relative humidity profile in the 1D model, which is used to determine the water vapor profile.

The ice-albedo feedback is a potentially important de-stabilizing effect for the climate of terrestrial planets. It is based on the positive feedback between decreasing surface temperatures, an increase of snow and ice cover and an associated increase in planetary albedo, which then further decreases surface temperature. A recent study shows that for M stars, the strength of the ice-albedo feedback is reduced This study investigates the influence of the atmosphere (in terms of surface pressure and atmospheric composition) for this feedback. A plane-parallel radiative transfer model is used for the calculation of planetary albedos. We varied CO2 partial pressures as well as the H2O, CH4, and O3 content in the atmosphere for planets orbiting Sun-like and M-type stars. Results suggest that for planets around M stars, the ice-albedo effect is significantly reduced, compared to planets around Sun-like stars. Including the effects of an atmosphere further suppresses the sensitivity to the ice-albedo effect. Atmospheric key properties such as surface pressure, but also the abundance of radiative trace gases can considerably change the strength of the ice-albedo feedback. For dense CO2 atmospheres of the order of a few to tens of bar, atmospheric rather than surface properties begin to dominate the planetary radiation budget. At high CO2 pressures, the ice-albedo feedback is strongly reduced for planets around M stars. The presence of trace amounts of H2O and CH4 in the atmosphere also weakens the ice-albedo effect for both stellar types considered. For planets around Sun-like stars, O3 could also lead to a very strong decrease of the ice-albedo feedback at high CO2 pressures. (abridged)

Owing to their wavelengths dependent absorption and scattering properties, clouds have a strong impact on the climate of planetary atmospheres. Especially, the potential greenhouse effect of CO2 ice clouds in the atmospheres of terrestrial extrasolar planets is of particular interest because it might influence the position and thus the extension of the outer boundary of the classic habitable zone around main sequence stars. We study the radiative effects of CO2 ice particles obtained by different numerical treatments to solve the radiative transfer equation. The comparison between the results of a high-order discrete ordinate method and simpler two-stream approaches reveals large deviations in terms of a potential scattering efficiency of the greenhouse effect. The two-stream methods overestimate the transmitted and reflected radiation, thereby yielding a higher scattering greenhouse effect. For the particular case of a cool M-type dwarf the CO2 ice particles show no strong effective scattering greenhouse effect by using the high-order discrete ordinate method, whereas a positive net greenhouse effect was found in case of the two-stream radiative transfer schemes. As a result, previous studies on the effects of CO2 ice clouds using two-stream approximations overrated the atmospheric warming caused by the scattering greenhouse effect. Consequently, the scattering greenhouse effect of CO2 ice particles seems to be less effective than previously estimated. In general, higher order radiative transfer methods are necessary to describe the effects of CO2 ice clouds accurately as indicated by our numerical radiative transfer studies.

We report here that the equation for H2O Rayleigh scattering was incorrectly stated in the original paper [arXiv:1009.5814]. Instead of a quadratic dependence on refractivity r, we accidentally quoted an r^4 dependence. Since the correct form of the equation was implemented into the model, scientific results are not affected.

The planetary system around the M star Gliese 581 contains at least three close-in potentially low-mass planets, GL 581 c, d, and e. In order to address the question of the habitability of GL 581 d, we performed detailed atmospheric modeling studies for several planetary scenarios. A 1D radiative-convective model was used to calculate temperature and pressure profiles of model atmospheres, assumed to be composed of molecular nitrogen, water, and carbon dioxide. The model allows for changing surface pressures caused by evaporation/condensation of water and carbon dioxide. Furthermore, the treatment of the energy transport has been improved in the model to account in particular for high CO2, high-pressure Super-Earth conditions. For four high-pressure scenarios of our study, the resulting surface temperatures were above 273 K, indicating a potential habitability of the planet. These scenarios include three CO2-dominated atmospheres (95% CO2 concentration with 5, 10, and 20 bar surface pressure) and a high-pressure CO2-enriched atmosphere (5% CO2 concentration with 20 bar surface pressure). For all other considered scenarios, the calculated GL 581 d surface temperatures were below the freezing point of water, suggesting that GL 581 d would not be habitable then. The results for our CO2-dominated scenarios confirm very recent model results by Wordsworth et al. (2010). However, our model calculations imply that also atmospheres that are not CO2-dominated (i.e., 5% vmr instead of 95% vmr) could result in habitable conditions for GL 581 d.

Clouds and hazes are commonplace in the atmospheres of solar system planets and are likely ubiquitous in the atmospheres of extrasolar planets as well. Clouds affect every aspect of a planetary atmosphere, from the transport of radiation, to atmospheric chemistry, to dynamics and they influence - if not control - aspects such as surface temperature and habitability. In this review we aim to provide an introduction to the role and properties of clouds in exoplanetary atmospheres. We consider the role clouds play in influencing the spectra of planets as well as their habitability and detectability. We briefly summarize how clouds are treated in terrestrial climate models and consider the far simpler approaches that have been taken so far to model exoplanet clouds, the evidence for which we also review. Since clouds play a major role in the atmospheres of certain classes of brown dwarfs we briefly discuss brown dwarf cloud modeling as well. We also review how the scattering and extinction efficiencies of cloud particles may be approximated in certain limiting cases of small and large particles in order to facilitate physical understanding. Since clouds play such important roles in planetary atmospheres, cloud modeling may well prove to be the limiting factor in our ability to interpret future observations of extrasolar planets.

We study the influence of low-level water and high-level ice clouds on low-resolution reflection spectra and planetary albedos of Earth-like planets orbiting different types of stars in both the visible and near infrared wavelength range. We use a one-dimensional radiative-convective steady-state atmospheric model coupled with a parametric cloud model, based on observations in the Earth's atmosphere to study the effect of both cloud types on the reflection spectra and albedos of Earth-like extrasolar planets at low resolution for various types of central stars. We find that the high scattering efficiency of clouds substantially causes both the amount of reflected light and the related depths of the absorption bands to be substantially larger than in comparison to the respective clear sky conditions. Low-level clouds have a stronger impact on the spectra than the high-level clouds because of their much larger scattering optical depth. The detectability of molecular features in near the UV - near IR wavelength range is strongly enhanced by the presence of clouds. However, the detectability of various chemical species in low-resolution reflection spectra depends strongly on the spectral energy distribution of the incident stellar radiation. In contrast to the reflection spectra the spectral planetary albedos enable molecular features to be detected without a direct influence of the spectral energy distribution of the stellar radiation. Here, clouds increase the contrast between the radiation fluxes of the planets and the respective central star by about one order of magnitude, but the resulting contrast values are still too low to be observable with the current generation of telescopes.

We study the impact of multi-layered clouds (low-level water and high-level ice clouds) on the thermal emission spectra of Earth-like planets orbiting different types of stars. Clouds have an important influence on such planetary emission spectra due to their wavelength dependent absorption and scattering properties. We also investigate the influence of clouds on the ability to derive information about planetary surface temperatures from low-resolution spectra.

The effects of multi-layered clouds in the atmospheres of Earth-like planets orbiting different types of stars are studied. The radiative effects of cloud particles are directly correlated with their wavelength-dependent optical properties. Therefore the incident stellar spectra may play an important role for the climatic effect of clouds. We discuss the influence of clouds with mean properties measured in the Earth's atmosphere on the surface temperatures and Bond albedos of Earth-like planets orbiting different types of main sequence dwarf stars.