Observations suggest an abundance of water and paucity of methane in the majority of observed exoplanetary atmospheres. We isolate the effect of atmospheric processes to investigate possible causes. Previously, we studied the effect of effective temperature, surface gravity, metallicity, carbon-to-oxygen ratio, and stellar type assuming cloud-free thermochemical equilibrium and disequilibrium chemistry. However, under these assumptions, methane remains a persisting spectral feature in the transmission spectra of exoplanets over a certain parameter space, the Methane Valley. In this work we investigate the role of clouds on this domain and we find that clouds change the spectral appearance of methane in two direct ways: 1) by heating-up the photosphere of colder planets, and 2) by obscuring molecular features. The presence of clouds also affects methane features indirectly: 1) cloud heating results in more evaporation of condensates and hence releases additional oxygen, causing water dominated spectra of colder carbon-poor exoplanets, and 2) HCN/CO production results in a suppression of depleted methane features by these molecules. The presence of HCN/CO and a lack of methane could be an indication of cloud formation on hot exoplanets. Cloud heating can also deplete ammonia. Therefore, a simultaneous depletion of methane and ammonia is not unique to photochemical processes. We propose that the best targets for methane detection are likely to be massive but smaller planets with a temperature around 1450 K orbiting colder stars. We also construct Spitzer synthetic color-maps and find that clouds can explain some of the high contrast observations by IRAC's channel 1 and 2.

Temperature inversion layers are predicted to be present in ultra-hot giant planet atmospheres. Although such inversion layers have been observed in several ultra-hot Jupiters recently, the chemical species responsible for creating the inversion remain unidentified. Here, we present observation of the thermal emission spectrum of an ultra-hot Jupiter, WASP-189b, at high spectral resolution using the HARPS-N spectrograph. Using the cross-correlation technique, we detect a strong Fe I signal. The detected Fe I spectral lines are found in emission, which is a direct evidence of a temperature inversion in the planetary atmosphere. We further performed a retrieval on the observed spectrum using a forward model with an MCMC approach. When assuming a solar metallicity, the best-fit result returns a temperature of $4320_{-100}^{+120}$ K at the top of the inversion, which is significantly hotter than the planetary equilibrium temperature (2641 K). The temperature at the bottom of the inversion is determined as $2200_{-800}^{+1000}$ K. Such a strong temperature inversion is probably created by the absorption of atomic species like Fe I.

Knowing the bulk elemental abundances of an exoplanetary atmosphere is not an easy task, but it is crucial to understand the formation history of planets. The purpose of this work is to show that the observability of TiO features at optical wavelengths in the transmission spectra of hot Jupiter atmospheres is sensitive to the bulk chemical properties of the atmosphere. For that, we run a grid of chemical models which include TiO formation and destruction, for the ultra-hot Jupiter WASP-19b and an ultra-hot version of HD~209458b. We take into account non-equilibrium chemistry and changes in the temperature and pressure structure of these atmospheres caused by different C/O ratios. We calculate synthetic transmission spectra for these models, and study the relative strengths of TiO and \ce{H2O} features quantitatively. To compare with observations, we use a model independent metric for molecular abundances, $\Delta Z_{\rm TiO-H_2O}/H_{\rm eq}$ that has been previously used in observational studies of exoplanetary atmospheres. We find that with this metric we can differentiate between different chemical models and place constraints on the atmosphere's bulk carbon and oxygen abundance. From chemical considerations we expected that the TiO abundance would depend on the bulk nitrogen, however found that any change to N/H did not result in changes to the resulting TiO. We apply our method to a set of known exoplanets that have been observed in the relevant optical wavelengths and find good agreement between low resolution observations and our model for WASP-121b, marginally good agreement with WASP-79b, WASP-76b, and WASP-19b, and poorer agreement with HD 209458b. Our method can be particularly helpful for indirect studies of the bulk abundances of carbon and oxygen.

We present WASP-43b climate simulations with deep wind jets (down to 700~bar) that are linked to retrograde (westward) flow at the equatorial day side for $p<0.1$~bar. Retrograde flow inhibits efficient eastward heat transport and naturally explains the small hotspot shift and large day-night-side gradient of WASP-43b ($P_{\text{orb}}=P_{\text{rot}}=0.8135$~days) observed with Spitzer. We find that deep wind jets are mainly associated with very fast rotations ($P_{\text{rot}}=P_{\text{orb}}\leq 1.5$~days) which correspond to the Rhines length smaller than $2$ planetary radii. We also diagnose wave activity that likely gives rise to deviations from superrotation. Further, we show that we can achieve full steady state in our climate simulations by imposing a deep forcing regime for $p>10$~bar: convergence time scale $\tau_{\text{conv}}=10^6-10^8$~s to a common adiabat, as well as linear drag at depth ($p\geq 200$~bar), which mimics to first order magnetic drag. Lower boundary stability and the deep forcing assumptions were also tested with climate simulations for HD~209458b ($P_{\text{orb}}=P_{\text{rot}}=3.5$~days). HD~209458b simulations always show shallow wind jets (never deeper than 100~bar) and unperturbed superrotation. If we impose a fast rotation ($P_{\text{orb}}=P_{\text{rot}}=0.8135$~days), also the HD~209458b-like simulation shows equatorial retrograde flow at the day side. We conclude that the placement of the lower boundary at $p=200$~bar is justified for slow rotators like HD~209458b, but we suggest that it has to be placed deeper for fast-rotating, dense hot Jupiters ($P_{\text{orb}}\leq 1.5$~days) like WASP-43b. Our study highlights that the deep atmosphere may have a strong influence on the observable atmospheric flow in some hot Jupiters.

Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our radiative transfer code petitRADTRANS for generating spectra, which we couple to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that ${\rm C/O}=0.60_{-0.08}^{+0.07}$. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: cloudy, or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, while not constraining C/O. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO$_2$ or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.

We present a transmission spectrum for the Neptune-size exoplanet HD 106315 c from optical to infrared wavelengths based on transit observations from the Hubble Space Telescope/Wide Field Camera 3, K2, and Spitzer. The spectrum shows tentative evidence for a water absorption feature in the $1.1 - 1.7\mu$m wavelength range with a small amplitude of 30 ppm (corresponding to just $0.8 \pm 0.04$ atmospheric scale heights). Based on an atmospheric retrieval analysis, the presence of water vapor is tentatively favored with a Bayes factor of 1.7 - 2.6 (depending on prior assumptions). The spectrum is most consistent with either enhanced metallicity, high altitude condensates, or both. Cloud-free solar composition atmospheres are ruled out at $>5\sigma$ confidence. We compare the spectrum to grids of cloudy and hazy forward models and find that the spectrum is fit well by models with moderate cloud lofting or haze formation efficiency, over a wide range of metallicities ($1 - 100\times$ solar). We combine the constraints on the envelope composition with an interior structure model and estimate that the core mass fraction is $\gtrsim0.3$. With a bulk composition reminiscent of that of Neptune and an orbital distance of 0.15 AU, HD 106315 c hints that planets may form out of broadly similar material and arrive at vastly different orbits later in their evolution.

{We investigate the atmospheric composition of the long period ($P_{\rm orb}=$~10 days), eccentric exo-Saturn WASP-117b. WASP-117b could be in atmospheric temperature and chemistry similar to WASP-107b. In mass and radius WASP-117b is similar to WASP-39b, which allows a comparative study of these planets.} % methods heading (mandatory) { We analyze a near-infrared transmission spectrum of WASP-117b taken with Hubble Space Telescope/WFC3 G141, which was reduced with two independent pipelines to ensure a robust detection of the $3 \sigma$-water spectrum. High resolution measurements were taken with VLT/ESPRESSO in the optical.} % results heading (mandatory) {Using a 1D atmosphere model with isothermal temperature, uniform cloud deck and equilibrium chemistry, the Bayesian evidence of a retrieval analysis of the transmission spectrum indicates a preference for a high atmospheric metallicity ${\rm [Fe/H]}=2.58^{+0.26}_{-0.37}$ and clear skies. The data are also consistent with a lower-metallicity composition ${\rm [Fe/H]}<1.75$ and a cloud deck between $10^{-2.2} - 10^{-5.1}$~bar, but with weaker Bayesian preference. We retrieve a low CH$_4$ abundance of $<10^{-4}$ volume fraction within $1 \sigma$. We cannot constrain the temperature between theoretically imposed limits of 700 and 1000~K. Further observations are needed to confirm quenching of CH$_4$ with $K_{zz}\geq 10^8$~cm$^2$/s. We report indications of Na and K in the VLT/ESPRESSO high resolution spectrum with substantial Bayesian evidence in combination with HST data.}

Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studies. Aims. We aim to further characterize the orbital, spectral, and physical properties of the HD 19467 system. Methods. We present new high-contrast imaging data with the SPHERE and NaCo instruments. We also analyze archival data from HARPS, NaCo, HIRES, UVES, and ASAS. We also use proper motion data of the star from Hipparcos and Gaia. Results. We refine the properties of the host star and derive an age of 8.0$^{+2.0}_{-1.0}$ Gyr based on isochrones, gyrochronology, and chemical and kinematic arguments. This estimate is slightly younger than previous estimates of ~9-11 Gyr. No orbital curvature is seen in the current imaging, RV, and astrometric data. From a joint fit of the data, we refine the orbital parameters for HD 19467B: period 398$^{+95}_{-93}$ yr, inclination 129.8$^{+8.1}_{-5.1}$ deg, eccentricity 0.56$\pm$0.09, longitude of the ascending node 134.8$\pm$4.5 deg, and argument of the periastron 64.2$^{+5.5}_{-6.3}$ deg. We assess a dynamical mass of 74$^{+12}_{-9}$ MJ. The fit with atmospheric models of the spectrophotometric data of HD 19467B indicates an atmosphere without clouds or with very thin clouds, an effective temperature of 1042$^{+77}_{-71}$ K, and a large surface gravity of 5.34$^{+0.08}_{-0.09}$ dex. The comparison to model predictions of the bolometric luminosity and dynamical mass of HD 19467B, assuming our system age estimate, indicates a better agreement with the Burrows et al. models; whereas the other evolutionary models used tend to underestimate its cooling rate.

We report the discovery of a warm Neptune and a hot sub-Neptune transiting TOI-421 (BD-14 1137, TIC 94986319), a bright (V=9.9) G9 dwarf star in a visual binary system observed by the TESS space mission in Sectors 5 and 6. We performed ground-based follow-up observations -- comprised of LCOGT transit photometry, NIRC2 adaptive optics imaging, and FIES, CORALIE, HARPS, HIRES, and PFS high-precision Doppler measurements -- and confirmed the planetary nature of the 16-day transiting candidate announced by the TESS team. We discovered an additional radial velocity signal with a period of 5 days induced by the presence of a second planet in the system, which we also found to transit its host star. We found that the outer warm Neptune, TOI-421 b, has an orbital period of Pb=16.06815+0.00034-0.00035 days, a mass of Mb=16.23+1.14-1.08 MEarth, a radius of Rb=5.17+-0.13 REarth and a density of rho_b=0.645+0.069-0.062 g/cm3, whereas the inner mini-Neptune, TOI-421 c, has a period of Pc=5.19676+0.00049-0.00048 days, a mass of Mc=7.05+0.71-0.70 MEarth and a radius of Rc=2.72+0.19-0.18 REarth. With its characteristics the latter planet (rho_c=1.93+0.49-0.39 g/cm3) is placed in the intriguing class of the super-puffy mini-Neptunes. TOI-421 b and c are found to be well suitable for atmospheric characterization. Our atmospheric simulations predict significant Ly-alpha transit absorption, due to strong hydrogen escape in both planets, and the presence of detectable CH4 in the atmosphere of TOI-421 b if equilibrium chemistry is assumed.

KELT-9 b, the hottest known exoplanet with $T\sim4400$ K, is the archetype of a new planet class known as ultra-hot Jupiters. These exoplanets are presumed to have an atmosphere dominated by neutral and ionized atomic species. In particular, H$\alpha$ and H$\beta$ Balmer lines have been detected in the KELT-9 b upper atmosphere, suggesting that hydrogen is filling the planetary Roche lobe and escaping from the planet. In this work, we detected $\delta$ Scuti-type stellar pulsation (with a period $P=7.54\pm0.12$ h) and studied the Rossiter-McLaughlin effect (finding a spin-orbit angle $\lambda=-85.01^{\deg}\pm0.23^{\deg}$) prior to focussing on the Balmer lines (H$\alpha$ to H$\zeta$) in the optical transmission spectrum of KELT-9 b. Our HARPS-N data show significant absorption for H$\alpha$ to H$\delta$. The precise line shapes of the H$\alpha$, H$\beta$, and H$\gamma$ absorptions allow us to put constraints on the thermospheric temperature. Moreover, the mass loss rate, and the excited hydrogen population of KELT-9 b are also constrained, thanks to a retrieval analysis performed with a new atmospheric model. We retrieved a thermospheric temperature of $T=13200^{+800}_{-720}$ K and a mass loss rate of $\dot{M}=10^{12.8\pm0.3}$ g s$^{-1}$ when the atmosphere was assumed to be in hydrodynamical expansion and in local thermodynamic equilibrium (LTE). Since the thermospheres of hot Jupiters are not expected to be in LTE, we explored atmospheric structures with non-Boltzmann equilibrium for the population of the excited hydrogen. We do not find strong statistical evidence in favor of a departure from LTE. However, our non-LTE scenario suggests that a departure from the Boltzmann equilibrium may not be sufficient to explain the retrieved low number densities of the excited hydrogen. In non-LTE, Saha equilibrium departure via photo-ionization, is also likely to be necessary to explain the data.

We explore the chemistry and observability of nitrogen dominated atmospheres for ultra-short-period super-Earths. We base the assumption, that super-Earths could have nitrogen filled atmospheres, on observations of 55 Cnc e that favour a scenario with a high-mean-molecular-weight atmosphere. We take Titan's elemental budget as our starting point and using chemical kinetics compute a large range of possible compositions for a hot super-Earth. We use analytical temperature profiles and explore a parameter space spanning orders of magnitude in C/O \& N/O ratios, while always keeping nitrogen the dominant component. We generate synthetic transmission and emission spectra and assess their potential observability with the future James Webb Space Telescope and ARIEL. Our results suggest that HCN is a strong indicator of a high C/O ratio, which is similar to what is found for H-dominated atmospheres. We find that these worlds are likely to possess C/O > 1.0, and that HCN, CN, CO should be the primary molecules to be searched for in thermal emission. For lower temperatures (T < 1500 K), we additionally find NH3 in high N/O ratio cases, and C2H4, CH4 in low N/O ratio cases to be strong absorbers. Depletion of hydrogen in such atmospheres would make CN, CO and NO exceptionally prominent molecules to look for in the 0.6 - 5.0 $\mu$m range. Our models show that the upcoming JWST and ARIEL missions will be able to distinguish atmospheric compositions of ultra-short period super-Earths with unprecedented confidence.

Directly imaged planets and substellar companions are key targets for the characterization of self-luminous atmospheres. Their photometric appearance at 4-5 $\mu$m is sensitive to the chemical composition and cloud content of their atmosphere. We aim at systematically characterizing the atmospheres of directly imaged low-mass companions at 4-5 $\mu$m. We want to homogeneously process the data and compile a photometric library at thermal wavelengths of these mostly young, low-gravity objects. In this way, we want to find trends related to their spectral type and surface gravity by comparing with isolated brown dwarfs and predictions from atmospheric models. We have used the high-resolution, high-contrast capabilities of NACO at the VLT to directly image the companions of HIP 65426, PZ Tel, and HD 206893 in the NB4.05 and/or $M'$ filters. For the same targets, and additionally $\beta$ Pic, we have also analyzed six archival VLT/NACO datasets which were taken with the NB3.74, $L'$, NB4.05, and $M'$ filters. The $L'$-NB4.05 and $L'$-$M'$ colors of the studied sample are all red while the NB4.05-$M'$ color is blue for $\beta$ Pic b, gray for PZ Tel B, and red for HIP 65426 b and HD 206893 B (although typically with low significance). The absolute NB4.05 and $M'$ fluxes of our sample are all larger than those of field dwarfs with similar spectral types. Finally, the surface gravity of $\beta$ Pic b has been constrained to $\log{g} = 4.17_{-0.13}^{+0.10}$ dex from its photometry and dynamical mass. A red color at 3-4 $\mu$m and a blue color at 4-5 $\mu$m might be (partially) caused by H$_2$O and CO absorption, respectively. The red characteristics of $\beta$ Pic b, HIP 65426 b, and HD 206893 B at 3-5$\mu$m, as well as their higher fluxes in NB4.05 and $M'$ compared to field dwarfs, indicate that cloud densities are enhanced close to the photosphere as a result of their low surface gravity.

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $\beta$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $\beta$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $\beta$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

Ultra-hot Jupiters are emerging as a new class of exoplanets. Studying their chemical compositions and temperature structures will improve the understanding of their mass loss rate as well as their formation and evolution. We present the detection of ionized calcium in the two hottest giant exoplanets - KELT-9b and WASP-33b. By utilizing transit datasets from CARMENES and HARPS-N observations, we achieved high confidence level detections of Ca II using the cross-correlation method. We further obtain the transmission spectra around the individual lines of the Ca II H&K doublet and the near-infrared triplet, and measure their line profiles. The Ca II H&K lines have an average line depth of 2.02 $\pm$ 0.17 % (effective radius of 1.56 Rp) for WASP-33b and an average line depth of 0.78 $\pm$ 0.04 % (effective radius of 1.47 Rp) for KELT-9b, which indicates that the absorptions are from very high upper atmosphere layers close to the planetary Roche lobes. The observed Ca II lines are significantly deeper than the predicted values from the hydrostatic models. Such a discrepancy is probably a result of hydrodynamic outflow that transports a significant amount of Ca II into the upper atmosphere. The prominent Ca II detection with the lack of significant Ca I detection implies that calcium is mostly ionized in the upper atmospheres of the two planets.

Almost all planetary atmospheres are affected by disequilibrium chemical processes. In this paper we introduce our recently developed Chemical Kinetic Model (\texttt{ChemKM}). We show that the results of our HD189733b model are in good agreement with previously published results, except at $\mu$bar regime, where molecular diffusion and photochemistry are the dominant processes. We thus recommend careful consideration of these processes when abundances at the top of the atmosphere are desired. We also propose a new metric for a quantitative measure of quenching levels. By applying this metric, we find that quenching pressure decreases with the effective temperature of planets, but it also varies significantly with other atmospheric parameters such as [Fe/H], log(g), and C/O. In addition, we find that the "Methane Valley", a region between 800 and 1500K where above a certain C/O threshold value a greater chance of CH$_4$ detection is expected, still exists after including the vertical mixing. The first robust CH$_4$ detection on an irradiated planet (HD102195b) places this object within this region; supporting our prediction. We also investigate the detectability of disequilibrium spectral fingerprints by JWST, and suggest focusing on the targets with T$_{eff}$ between 1000 and 1800K, orbiting around M-dwarfs, having low surface gravity but high metallicity and a C/O ratio value around unity. Finally, constructing Spitzer color-maps suggests that the main two color-populations are largely insensitive to the vertical mixing. Therefore any deviation of observational points from these populations are likely due to the presence of clouds and not disequilibrium processes. However, some cold planets (T$_{eff}<$900K) with very low C/O ratios ($<$0.25) show significant deviations; making these planets interesting cases for further investigation.

We present the easy-to-use, publicly available, Python package petitRADTRANS, built for the spectral characterization of exoplanet atmospheres. The code is fast, accurate, and versatile; it can calculate both transmission and emission spectra within a few seconds at low resolution ($\lambda/\Delta\lambda$ = 1000; correlated-k method) and high resolution ($\lambda/\Delta\lambda = 10^6$; line-by-line method), using only a few lines of input instruction. The somewhat slower correlated-k method is used at low resolution because it is more accurate than methods such as opacity sampling. Clouds can be included and treated using wavelength-dependent power law opacities, or by using optical constants of real condensates, specifying either the cloud particle size, or the atmospheric mixing and particle settling strength. Opacities of amorphous or crystalline, spherical or irregularly-shaped cloud particles are available. The line opacity database spans temperatures between 80 and 3000 K, allowing to model fluxes of objects such as terrestrial planets, super-Earths, Neptunes, or hot Jupiters, if their atmospheres are hydrogen-dominated. Higher temperature points and species will be added in the future, allowing to also model the class of ultra hot-Jupiters, with equilibrium temperatures $T_{\rm eq} \gtrsim 2000$ K. Radiative transfer results were tested by cross-verifying the low- and high-resolution implementation of petitRADTRANS, and benchmarked with the petitCODE, which itself is also benchmarked to the ATMO and Exo-REM codes. We successfully carried out test retrievals of synthetic JWST emission and transmission spectra (for the hot Jupiter TrES-4b, which has a $T_{\rm eq}$ of $\sim$ 1800 K). The code is publicly available at http://gitlab.com/mauricemolli/petitRADTRANS, and its documentation can be found at https://petitradtrans.readthedocs.io.

[Abridged] Aims. We searched for circumplanetary sodium and ionized calcium gas around the disintegrating rocky exoplanet K2-22 b to constrain its gas-loss and sublimation processes. Methods. We observed four transits of K2-22 b with X-shooter on ESO's Very Large Telescope to obtain time-series of intermediate-resolution (R $\sim$ 11400) spectra. Our analysis focused on the two sodium D lines (588.995 nm and 589.592 nm) and the Ca$^{+}$ triplet (849.802 nm, 854.209 nm and 866.214 nm). Planet-related absorption is searched for in the velocity rest frame of the planet, which changes from $\pm$66 kms$^{-1}$ during the transit. Results. Since K2-22 b exhibits highly variable transit depths, we analyzed the individual nights and their average. By injecting signals we reached 5$\sigma$ upper-limits on the individual nights that ranged from 11% - 13% and 1.7% - 2.0% for the tail's sodium and ionized calcium absorption, respectively. Night 1 was contaminated by its companion star so we considered weighted averages with and without Night 1 and quote conservative 5$\sigma$ limits without Night 1 of 9% and 1.4%, respectively. Assuming their mass fractions to be similar to those in the Earth's crust, these limits correspond to scenarios in which 0.04% and 35% of the transiting dust is sublimated and observed as absorbing gas. However, this assumes the gas to be co-moving with the planet. We show that for the high irradiation environment of K2-22 b, sodium and ionized calcium could be quickly accelerated to 100s of km s$^{-1}$ due to radiation pressure and entrainment by the stellar wind, making them much more difficult to detect. No evidence for such possibly broad and blue-shifted signals are seen in our data. Conclusions. Future observations aimed at observing circumplanetary gas should take into account the possible broad and blue-shifted velocity field of atomic and ionized species.

To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,\mu$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.

Future instruments like NIRCam and MIRI on JWST or METIS at the ELT will be able to image exoplanets that are too faint for current direct imaging instruments. Evolutionary models predicting the planetary intrinsic luminosity as a function of time have traditionally concentrated on gas-dominated giant planets. We extend these cooling curves to Saturnian and Neptunian planets. We simulate the cooling of isolated core-dominated and gas giant planets with masses of 5 Earthmasses to 2 Jupitermasses. The luminosity includes the contribution from the cooling and contraction of the core and of the H/He envelope, as well as radiogenic decay. For the atmosphere we use grey, AMES-Cond, petitCODE, and HELIOS models. We consider solar and non-solar metallicities as well as cloud-free and cloudy atmospheres. The most important initial conditions, namely the core-to-envelope ratio and the initial luminosity are taken from planet formation simulations based on the core accretion paradigm. We first compare our cooling curves for Uranus, Neptune, Jupiter, Saturn, GJ 436b, and a 5 Earthmass-planet with a 1% H/He envelope with other evolutionary models. We then present the temporal evolution of planets with masses between 5 Earthmasses and 2 Jupitermasses in terms of their luminosity, effective temperature, radius, and entropy. We discuss the impact of different post formation entropies. For the different atmosphere types and initial conditions magnitudes in various filter bands between 0.9 and 30 micrometer wavelength are provided. Using black body fluxes and non-grey spectra, we estimate the detectability of such planets with JWST. It is found that a 20 (100) Earthmass-planet can be detected with JWST in the background limit up to an age of about 10 (100) Myr with NIRCam and MIRI, respectively.

A carbon-to-oxygen ratio (C/O) of around unity is believed to act as a natural separator of water- and methane-dominated spectra when characterizing exoplanet atmospheres. In this paper we quantify the C/O ratios at which this separation occurs by calculating a large self-consistent grid of cloud-free atmospheric models in chemical equilibrium, using the latest version of petitCODE. Our study covers a broad range of parameter space: 400 K$<$T$_{eff}<$2600 K, 2.0$<$log(g)$<$5.0, -1.0$<$[Fe/H]$<$2.0, 0.25$<$C/O$<$1.25, and stellar types from M to F. We make the synthetic transmission and emission spectra, as well as the temperature structures publicly available. We find that the transition C/O ratio depends on many parameters such as effective temperature, surface gravity, metallicity and spectral type of the host star, and could have values less, equal, or higher than unity. By mapping all the transition C/O ratios we propose a "four-class" classification scheme for irradiated planets in this temperature range. We find a parameter space where methane always remains the cause of dominant spectral features. CH$_4$ detection in this region, or the lack of it, provides a diagnostic tool to identify the prevalence of cloud formation and non-equilibrium chemistry. As another diagnostic tool, we construct synthetic Spitzer IRAC color-diagrams showing two distinguishable populations of planets. Since most of the exoplanet atmospheres appear cloudy when studied in transmission, we regard this study as a starting point of how such a C/O-sensitive observation-based classification scheme should be constructed. This preparatory work will have to be refined by future cloudy and non-equilibrium modeling, to further investigate the existence and exact location of the classes, as well as the color-diagram analysis.

Cross-correlation is a well-tested method for exoplanet characterization. A new, potentially powerful application is the measurement of atmospheric isotope ratios. In particular D/H can give unique insights into a planet's formation and evolution. Here we aim to study the detectability of isotopologues in the high-dispersion spectra of exoplanets, to identify the optimal wavelengths ranges, and to predict the required observational efforts with current and future ground-based instruments. High-dispersion (R=10$^5$) thermal emission (and sometimes reflection) spectra were simulated by self-consistently modeling exoplanet atmospheres over a wide range of temperatures. These were synthetically observed with telescopes equivalent to the VLT or ELT, and analyzed with cross-correlation, resulting in S/N predictions for the detection of $^{13}$CO, HDO, and CH$_3$D. For the best observable exoplanets, $^{13}$CO is in range of current telescopes. It will be most favorably detected at 2.4 microns, just longward of the spectral range probed by several high-dispersion observations in the literature. CH$_3$D can best be seen at 4.7 microns, using 40m-class telescopes for planets with $T_{\rm equ}$ below 600 K. In this case, sky emission is often dominating the noise. HDO can be targeted at 3.7 microns, where sky emission is smaller. 40m-class telescopes may detect it in planets with $T_{\rm equ}$ below 900~K, potentially even 8m-class telescopes in the case of methane quenching. If Proxima Cen b is water-rich, HDO could be detected with the ELT in 1 night in reflected light. Isotopologues will soon belong to the exoplanet characterisation tools. Measuring D/H, and ratios of other isotopes, could be a prime science case for the METIS instrument on the ELT, especially for nearby rocky and ice giant planets. This can give unique insights in their history of ice enrichment and atmospheric evaporation.

Direct imaging has paved the way for atmospheric characterization of young and self-luminous gas giants. Scattering in a horizontally-inhomogeneous atmosphere causes the disk-integrated polarization of the thermal radiation to be linearly polarized, possibly detectable with the newest generation of high-contrast imaging instruments. We aim to investigate the effect of latitudinal and longitudinal cloud variations, circumplanetary disks, atmospheric oblateness, and cloud particle properties on the integrated degree and direction of polarization in the near-infrared. We have developed a three-dimensional Monte Carlo radiative transfer code (ARTES) for scattered light simulations in (exo)planetary atmospheres. The code is applicable to calculations of reflected light and thermal radiation in a spherical grid with a parameterized distribution of gas, clouds, hazes, and circumplanetary material. The disk-integrated degree of polarization of a horizontally-inhomogeneous atmosphere is maximal when the planet is flattened, the optical thickness of the equatorial clouds is large compared to the polar clouds, and the clouds are located at high altitude. For a flattened planet, the integrated polarization can both increase or decrease with respect to a spherical planet which depends on the horizontal distribution and optical thickness of the clouds. The direction of polarization can be either parallel or perpendicular to the projected direction of the rotation axis when clouds are zonally distributed. Rayleigh scattering by submicron-sized cloud particles will maximize the polarimetric signal whereas the integrated degree of polarization is significantly reduced with micron-sized cloud particles as a result of forward scattering. The presence of a cold or hot circumplanetary disk may also produce a detectable degree of polarization ($\lesssim$1%) even with a uniform cloud layer in the atmosphere.

Young circumstellar disks are of prime interest to understand the physical and chemical conditions under which planet formation takes place. Only very few detections of planet candidates within these disks exist, and most of them are currently suspected to be disk features. In this context, the transition disk around the young star PDS 70 is of particular interest, due to its large gap identified in previous observations, indicative of ongoing planet formation. We aim to search for the presence of planets and search for disk structures indicative for disk-planet interactions and other evolutionary processes. We analyse new and archival near-infrared (NIR) images of the transition disk PDS 70 obtained with the VLT/SPHERE, VLT/NaCo and Gemini/NICI instruments in polarimetric differential imaging (PDI) and angular differential imaging (ADI) modes. We detect a point source within the gap of the disk at about 195 mas (about 22 au) projected separation. The detection is confirmed at five different epochs, in three filter bands and using different instruments. The astrometry results in an object of bound nature, with high significance. The comparison of the measured magnitudes and colours to evolutionary tracks suggests that the detection is a companion of planetary mass. We confirm the detection of a large gap of about 54 au in size within the disk in our scattered light images, and detect a signal from an inner disk component. We find that its spatial extent is very likely smaller than about 17 au in radius. The images of the outer disk show evidence of a complex azimuthal brightness distribution which may in part be explained by Rayleigh scattering from very small grains. Future observations of this system at different wavelengths and continuing astrometry will allow us to test theoretical predictions regarding planet-disk interactions, planetary atmospheres and evolutionary models.

The G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model.

We present high-precision photometry of eight separate transit events in the HAT-P-32 planetary system. One transit event was observed simultaneously by two telescopes of which one obtained a simultaneous multi-band light curve in three optical bands, giving a total of 11 transit light curves. Due to the filter selection and in conjunction with using the defocussed photometry technique we were able to obtain an extremely high precision, ground-based transit in the \textit{u}-band (350\,nm), with an rms scatter of $\approx 1$\,mmag. All 11 transits were modelled using \textsc{prism} and \textsc{gemc}, and the physical properties of the system calculated. We find the mass and radius of the host star to be $1.182\pm 0.041\Msun$ and $1.225\pm0.015\Rsun$, respectively. For the planet we find a mass of $0.80\pm 0.14\Mjup$, a radius of $1.807\pm0.022\Rjup$ and a density of $0.126\pm0.023\pjup$. These values are consistent with those found in the literature. We also obtain a new orbital ephemeris for the system $ T_0 = {\rm BJD/TDB} \,\, 2\,454\,420.447187 (96) \, + \,2.15000800 (10) \times E $. We measured the transmission spectrum of HAT-P-32\,A\,b and compared it to theoretical transmission spectra. Our results indicate a bimodal cloud particle distribution consisting of Rayleigh--like haze and grey absorbing cloud particles within the atmosphere of HAT-P-32\,A\,b.

This paper continues a series in which we predict the main observable characteristics of exoplanets based on their formation. In Paper I we described our global planet formation and evolution model. In Paper II we studied the planetary mass-radius relationship. Here we present an extensive study of the statistics of planetary luminosities during both formation and evolution. Our results can be compared with individual directly imaged (proto)planets as well as statistical results from surveys. We calculated three synthetic planet populations assuming different efficiencies of the accretional heating by gas and planetesimals. We describe the temporal evolution of the planetary mass-luminosity relation. We study the shock and internal luminosity during formation. We predict a statistical version of the post-formation mass versus entropy "tuning fork" diagram. We find high nominal post-formation luminosities for hot and cold gas accretion. Individual formation histories can still lead to a factor of a few spread in the post-formation luminosity at a given mass. However, if the gas and planetesimal accretional heating is unknown, the post-formation luminosity may exhibit a spread of as much as 2-3 orders of magnitude at a fixed mass covering cold, warm, and hot states. As a key result we predict a flat log-luminosity distribution for giant planets, and a steep increase towards lower luminosities due to the higher occurrence rate of low-mass planets. Future surveys may detect this upturn. During formation an estimate of the planet mass may be possible for cold gas accretion if the gas accretion rate can be estimated. Due to the "core-mass effect" planets that underwent cold gas accretion can still have high post-formation entropies. Once the number of directly imaged exoplanets with known ages and luminosities increases, the observed distributions may be compared with our predictions.

51 Eridani b is an exoplanet around a young (20 Myr) nearby (29.4 pc) F0-type star, recently discovered by direct imaging. Being only 0.5" away from its host star it is well suited for spectroscopic analysis using integral field spectrographs. We aim to refine the atmospheric properties of this and to further constrain the architecture of the system by searching for additional companions. Using the SPHERE instrument at the VLT we extend the spectral coverage of the planet to the complete Y- to H-band range and provide photometry in the K12-bands (2.11, 2.25 micron). The object is compared to other cool and peculiar dwarfs. Furthermore, the posterior probability distributions of cloudy and clear atmospheric models are explored using MCMC. We verified our methods by determining atmospheric parameters for the two benchmark brown dwarfs Gl 570D and HD 3651B. For probing the innermost region for additional companions, archival VLT-NACO (L') SAM data is used. We present the first spectrophotometric measurements in the Y- and K-bands for the planet and revise its J-band flux to values 40% fainter than previous measurements. Cloudy models with uniform cloud coverage provide a good match to the data. We derive the temperature, radius, surface gravity, metallicity and cloud sedimentation parameter f_sed. We find that the atmosphere is highly super-solar (Fe/H~1.0) with an extended, thick cloud cover of small particles. The model radius and surface gravity suggest planetary masses of about 9 M_jup. The evolutionary model only provides a lower mass limit of >2 M_jup (for pure hot-start). The cold-start model cannot explain the planet's luminosity. The SPHERE and NACO/SAM detection limits probe the 51 Eri system at Solar System scales and exclude brown-dwarf companions more massive than 20 M_jup beyond separations of ~2.5 au and giant planets more massive than 2 M_jup beyond 9 au.

The James Webb Space Telescope will enable astronomers to obtain exoplanet spectra of unprecedented precision. Especially the MIRI instrument may shed light on the nature of the cloud particles obscuring planetary transmission spectra in the optical and near-infrared. We provide self-consistent atmospheric models and synthetic JWST observations for prime exoplanet targets in order to identify spectral regions of interest and estimate the number of transits needed to distinguish between model setups. We select targets which span a wide range in planetary temperature and surface gravity, ranging from super-Earths to giant planets, and have a high expected SNR. For all targets we vary the enrichment, C/O ratio, presence of optical absorbers (TiO/VO) and cloud treatment. We calculate atmospheric structures and emission and transmission spectra for all targets and use a radiometric model to obtain simulated observations. We analyze JWST's ability to distinguish between various scenarios. We find that in very cloudy planets such as GJ 1214b less than 10 transits with NIRSpec may be enough to reveal molecular features. Further, the presence of small silicate grains in atmospheres of hot Jupiters may be detectable with a single JWST MIRI transit. For a more detailed characterization of such particles less than 10 transits are necessary. Finally, we find that some of the hottest hot Jupiters are well fitted by models which neglect the redistribution of the insolation and harbor inversions, and that 1-4 eclipse measurements with NIRSpec are needed to distinguish between the inversion models. Wet thus demonstrate the capabilities of JWST for solving some of the most intriguing puzzles in current exoplanet atmospheric research. Further, by publishing all models calculated for this study we enable the community to carry out similar or retrieval analyses for all planets included in our target list.

The composition of a planet's atmosphere is determined by its formation, evolution, and present-day insolation. A planet's spectrum therefore may hold clues on its origins. We present a "chain" of models, linking the formation of a planet to its observable present-day spectrum. The chain links include (1) the planet's formation and migration, (2) its long-term thermodynamic evolution, (3) a variety of disk chemistry models, (4) a non-gray atmospheric model, and (5) a radiometric model to obtain simulated spectroscopic observations with JWST and ARIEL. In our standard chemistry model the inner disk is depleted in refractory carbon as in the Solar System and in white dwarfs polluted by extrasolar planetesimals. Our main findings are: (1) Envelope enrichment by planetesimal impacts during formation dominates the final planetary atmospheric composition of hot Jupiters. We investigate two, under this finding, prototypical formation pathways: a formation inside or outside the water iceline, called "dry" and "wet" planets, respectively. (2) Both the "dry" and "wet" planets are oxygen-rich (C/O<1) due to the oxygen-rich nature of the solid building blocks. The "dry" planet's C/O ratio is <0.2 for standard carbon depletion, while the "wet" planet has typical C/O values between 0.1 and 0.5 depending mainly on the clathrate formation efficiency. Only non-standard disk chemistries without carbon depletion lead to carbon-rich C/O ratios >1 for the "dry" planet. (3) While we consistently find C/O ratios <1, they still vary significantly. To link a formation history to a specific C/O, a better understanding of the disk chemistry is thus needed.

The growing body of observational data on extrasolar planets and protoplanetary disks has stimulated intense research on planet formation and evolution in the past few years. The extremely diverse, sometimes unexpected physical and orbital characteristics of exoplanets lead to frequent updates on the mainstream scenarios for planet formation and evolution, but also to the exploration of alternative avenues. The aim of this review is to bring together classical pictures and new ideas on the formation, orbital and internal evolutions of planets, highlighting the key role of the protoplanetary disk in the various parts of the theory. We begin by briefly reviewing the conventional mechanism of core accretion by the growth of planetesimals, and discuss a relatively recent model of core growth through the accretion of pebbles. We review the basic physics of planet-disk interactions, recent progress in this area, and discuss their role in observed planetary systems. We address the most important effects of planets internal evolution, like cooling and contraction, the mass-luminosity relation, and the bulk composition expressed in the mass-radius and mass-mean density relations.

We present broad-band photometry of five transits in the planetary system WASP-36, totaling 17 high-precision light curves. Four of the transits were simultaneously observed in four passbands (g, r, i, z), using the telescope-defocussing technique, and achieving scatters of less than 1 mmag per observation. We used these data to improve the measured orbital and physical properties of the system, and obtain an optical transmission spectrum of the planet. We measured a decreasing radius from bluer to redder passbands with a confidence level of more than 5 sigma. The radius variation is roughly 11 pressure scale heights between the g and the z bands. This is too strong to be Rayleigh scattering in the planetary atmosphere, and implies the presence of a species which absorbs strongly at bluer wavelengths.

Many parameters constraining the spectral appearance of exoplanets are still poorly understood. We therefore study the properties of irradiated exoplanet atmospheres over a wide parameter range including metallicity, C/O ratio and host spectral type. We calculate a grid of 1-d radiative-convective atmospheres and emission spectra. We perform the calculations with our new Pressure-Temperature Iterator and Spectral Emission Calculator for Planetary Atmospheres (PETIT) code, assuming chemical equilibrium. The atmospheric structures and spectra are made available online. We find that atmospheres of planets with C/O ratios $\sim$ 1 and $T_{\rm eff}$ $\gtrsim$ 1500 K can exhibit inversions due to heating by the alkalis because the main coolants CH$_4$, H$_2$O and HCN are depleted. Therefore, temperature inversions possibly occur without the presence of additional absorbers like TiO and VO. At low temperatures we find that the pressure level of the photosphere strongly influences whether the atmospheric opacity is dominated by either water (for low C/O) or methane (for high C/O), or both (regardless of the C/O). For hot, carbon-rich objects this pressure level governs whether the atmosphere is dominated by methane or HCN. Further we find that host stars of late spectral type lead to planetary atmospheres which have shallower, more isothermal temperature profiles. In agreement with prior work we find that for planets with $T_{\rm eff}$ $<$ 1750 K the transition between water or methane dominated spectra occurs at C/O $\sim$ 0.7, instead of $\sim$ 1, because condensation preferentially removes oxygen.

Despite the increase in observational data on exoplanets, the processes that lead to the formation of planets are still not well understood. But thanks to the high number of known exoplanets, it is now possible to look at them as a population that puts statistical constraints on theoretical models. A method that uses these constraints is planetary population synthesis. Its key element is a global model of planet formation and evolution that directly predicts observable planetary properties based on properties of the natal protoplanetary disk. To do so, global models build on many specialized models that address one specific physical process. We thoroughly review the physics of the sub-models included in global formation models. The sub-models can be classified as models describing the protoplanetary disk (gas and solids), the (proto)planet (solid core, gaseous envelope, and atmosphere), and finally the interactions (migration and N-body interaction). We compare the approaches in different global models and identify physical processes that require improved descriptions in future. We then address important results of population synthesis like the planetary mass function or the mass-radius relation. In these results, the global effects of physical mechanisms occurring during planet formation and evolution become apparent, and specialized models describing them can be put to the observational test. Due to their nature as meta models, global models depend on the development of the field of planet formation theory as a whole. Because there are important uncertainties in this theory, it is likely that global models will in future undergo significant modifications. Despite this, they can already now yield many testable predictions. With future global models addressing the geophysical characteristics, it should eventually become possible to make predictions about the habitability of planets.

Aims. Our aim is to study deuterium burning in objects forming according to the core accretion scenario in the hot and cold start assumption and what minimum deuterium burning mass limit is found for these objects. We also study how the burning process influences the structure and luminosity of the objects. Furthermore we want to test and verify our results by comparing them to already existing hot start simulations which did not consider, however, the formation process. Methods. We present a new method to calculate deuterium burning of objects in a self-consistently coupled model of planet formation and evolution. We discuss which theory is used to describe the process of deuterium burning and how it was implemented. Results. We find that the objects forming according to a hot start scenario behave approximately in the same way as found in previous works of evolutionary calculations, which did not consider the formation. However, for cold start objects one finds that the objects expand during deuterium burning instead of being partially stabilized against contraction. In both cases, hot and cold start, the mass of the solid core has an influence on the minimum mass limit of deuterium burning. The general position of the mass limit, 13 MJ, stays however approximately the same. None of the investigated parameters was able to change this mass limit by more than 0.8 MJ. Due to deuterium burning, the luminosity of hot and cold start objects becomes comparable after ~ 200 Myrs.