Search

[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., Comment: Accepted on 7 June 2019 for publication in Astronomy and Astrophysics (A&A). 17 pages, 11 figures. Submission updated after language editing by A&A

We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize their atmospheric structure, composition, and circulation, from gas giants to rocky planets. However, exoplanet atmospheric models capable of interpreting the upcoming observations are often limited by insufficiencies in the laboratory and theoretical data that serve as critical inputs to atmospheric physical and chemical tools. Here we provide an up-to-date and condensed description of areas where laboratory and/or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry, building off a larger 2016 white paper, and endorsed by the NAS Exoplanet Science Strategy report. Now is the ideal time for progress in these areas, but this progress requires better access to, understanding of, and training in the production of spectroscopic data as well as a better insight into chemical reaction kinetics both thermal and radiation-induced at a broad range of temperatures. Given that most published efforts have emphasized relatively Earth-like conditions, we can expect significant and enlightening discoveries as emphasis moves to the exotic atmospheres of exoplanets., Comment: Submitted as an Astro2020 Science White Paper

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., Comment: Accepted by MNRAS (36 pages with Appendix)

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., Comment: 17 pages, 7 figures, published in A&A

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., Comment: published in A&A

We present broad-band photometry of eleven planetary transits of the hot Jupiter WASP-74b, using three medium-class telescopes and employing the telescope-defocussing technique. Most of the transits were monitored through I filters and one was simultaneously observed in five optical (U, g', r', i', z') and three near infrared (J, H, K) passbands, for a total of 18 light curves. We also obtained new high-resolution spectra of the host star. We used these new data to review the orbital and physical properties of the WASP-74 planetary system. We were able to better constrain the main system characteristics, measuring smaller radius and mass for both the hot Jupiter and its host star than previously reported in the literature. Joining our optical data with those taken with the HST in the near infrared, we built up an observational transmission spectrum of the planet, which suggests the presence of strong optical absorbers, as TiO and VO gases, in its atmosphere., Comment: 12 pages, 9 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

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., Comment: Language corrected version and improved arrangements of figures, online data at: http://www.space.unibe.ch/research/research_groups/planets_in_time/numerical_data/index_eng.html

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., Comment: 27 pages, 12 figures, accepted in ApJ

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., Comment: 22 pages, 12 figures, updated version, accepted for publication in A & A

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., Comment: 33 pages, 25 figures. Accepted in Astronomy and Astrophysics

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., Comment: 23 pages, accepted by A&A

Aims: We aim to characterize the orbital and atmospheric properties of PDS 70 b, which was first identified on May 2015 in the course of the SHINE survey with SPHERE, the extreme adaptive-optics instrument at the VLT. Methods: We obtained new deep SPHERE/IRDIS imaging and SPHERE/IFS spectroscopic observations of PDS 70 b. The astrometric baseline now covers 6 years which allows us to perform an orbital analysis. For the first time, we present spectrophotometry of the young planet which covers almost the entire near-infrared range (0.96 to 3.8 micrometer). We use different atmospheric models covering a large parameter space in temperature, log(g), chemical composition, and cloud properties to characterize the properties of the atmosphere of PDS 70 b. Results: PDS 70 b is most likely orbiting the star on a circular and disk coplanar orbit at ~22 au inside the gap of the disk. We find a range of models that can describe the spectrophotometric data reasonably well in the temperature range between 1000-1600 K and log(g) no larger than 3.5 dex. The planet radius covers a relatively large range between 1.4 and 3.7 R_jupiter with the larger radii being higher than expected from planet evolution models for the age of the planet of 5.4 Myr. Conclusions: This study provides a comprehensive dataset on the orbital motion of PDS 70 b, indicating a circular orbit and a motion coplanar with the disk. The first detailed spectral energy distribution of PDS 70 b indicates a temperature typical for young giant planets. The detailed atmospheric analysis indicates that a circumplanetary disk may contribute to the total planet flux., Comment: accepted for publication in A&A (10 pages, 8 figures)

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., Comment: 17 pages, 7 figures and 8 tables. Accepted for publication in MNRAS on 2017 November 29

Given the forthcoming launch of the James Webb Space Telescope (JWST) which will allow observ- ing exoplanet atmospheres with unprecedented signal-over-noise ratio, spectral coverage and spatial resolution, the uncertainties in the atmosphere modelling used to interpret the data need to be as- sessed. As the first step, we compare three independent 1D radiative-convective models: ATMO, Exo-REM and petitCODE. We identify differences in physical and chemical processes taken into ac- count thanks to a benchmark protocol we developed. We study the impact of these differences on the analysis of observable spectra. We show the importance of selecting carefully relevant molecular linelists to compute the atmospheric opacity. Indeed, differences between spectra calculated with Hitran and ExoMol exceed the expected uncertainties of future JWST observations. We also show the limitation in the precision of the models due to uncertainties on alkali and molecule lineshape, which induce spectral effects also larger than the expected JWST uncertainties. We compare two chemical models, Exo-REM and Venot Chemical Code, which do not lead to significant differences in the emission or transmission spectra. We discuss the observational consequences of using equilibrium or out-of- equilibrium chemistry and the major impact of phosphine, detectable with the JWST.Each of the models has benefited from the benchmarking activity and has been updated. The protocol developed in this paper and the online results can constitute a test case for other models., Comment: 36 pages, 30 figures, accepted by ApJ

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., Comment: 44 pages, 26 figures (journal format). A&A in print. Language correction only relative to V1

The SHINE program is a large high-contrast near-infrared survey of 600 young, nearby stars. It is aimed at searching for and characterizing new planetary systems using VLT/SPHERE's unprecedented high-contrast and high-angular resolution imaging capabilities. It also intends at placing statistical constraints on the occurrence and orbital properties of the giant planet population at large orbits as a function of the stellar host mass and age to test planet formation theories. We use the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images and spectra of the young A2 star HIP65426. It is a member of the ~17 Myr old Lower Centaurus-Crux association. At a separation of 830 mas (92 au projected) from the star, we detect a faint red companion. Multi-epoch observations confirm that it shares common proper motion with HIP65426. Spectro-photometric measurements extracted with IFS and IRDIS between 0.95 and 2.2um indicate a warm, dusty atmosphere characteristic of young low surface-gravity L5-L7 dwarfs. Hot-start evolutionary models predict a luminosity consistent with a 6-12 MJup, Teff=1300-1600 K and R=1.5 RJup giant planet. Finally, the comparison with Exo-REM and PHOENIX BT-Settl synthetic atmosphere models gives consistent effective temperatures but with slightly higher surface gravity solutions of log(g)=4.0-5.0 with smaller radii (1.0-1.3 RJup). Given its physical and spectral properties, HIP65426b occupies a rather unique placement in terms of age, mass and spectral-type among the currently known imaged planets. It represents a particularly interesting case to study the presence of clouds as a function of particle size, composition, and location in the atmosphere, to search for signatures of non-equilibrium chemistry, and finally to test the theory of planet formation and evolution., Comment: 10 pages, 8 figures, 4 tables, accepted for publication in A&A, file updated

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., Comment: 18 pages, 9 figures, accepted for publication in A&A

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., Comment: 29 pages, 31 figures, accepted for publication in A&A

Detecting the atmospheres of low-mass low-temperature exoplanets is a high-priority goal on the path to ultimately detect biosignatures in the atmospheres of habitable exoplanets. High-precision HST observations of several super-Earths with equilibrium temperatures below 1000K have to date all resulted in featureless transmission spectra, which have been suggested to be due to high-altitude clouds. We report the detection of an atmospheric feature in the atmosphere of a 1.6 Mearth transiting exoplanet, GJ 1132b, with an equilibrium temperature of ~600K and orbiting a nearby M dwarf. We present observations of nine transits of the planet obtained simultaneously in the griz and JHK passbands. We find an average radius of 1.43 +/- 0.16 Rearth for the planet, averaged over all the passbands, and a radius of 0.255 +/- 0.023 Rsun for the star, both of which are significantly greater than previously found. The planet radius can be decomposed into a "surface radius" at ~1.375 Rearth overlaid by atmospheric features which increase the observed radius in the z and K bands. The z-band radius is 4sigma higher than the continuum, suggesting a strong detection of an atmosphere. We deploy a suite of tests to verify the reliability of the transmission spectrum, which are greatly helped by the existence of repeat observations. The large z-band transit depth indicates strong opacity from H2O and/or CH4 or a hitherto unconsidered opacity. A surface radius of 1.375 +/- 0.16 Rearth allows for a wide range of interior compositions ranging from a nearly Earth-like rocky interior, with ~70% silicate and ~30% Fe, to a substantially H2O-rich water world., Comment: 13 pages, 10 colour figures, 5 tables. Submitted on 2016/12/05. Version 2 is the version accepted for publication by AJ on 2017/02/28

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., Comment: 24 pages, 7 figures, accepted for publication in A&A

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., Comment: 32 pages, 18 figures. Accepted for publication in ApJ

We report 13 high-precision light curves of eight transits of the exoplanet WASP-52b, obtained by using four medium-class telescopes, through different filters, and adopting the defocussing technique. One transit was recorded simultaneously from two different observatories and another one from the same site but with two different instruments, including a multi-band camera. Anomalies were clearly detected in five light curves and modelled as starspots occulted by the planet during the transit events. We fitted the clean light curves with the jktebop code, and those with the anomalies with the prism+gemc codes in order to simultaneously model the photometric parameters of the transits and the position, size and contrast of each starspot. We used these new light curves and some from the literature to revise the physical properties of the WASP-52 system. Starspots with similar characteristics were detected in four transits over a period of 43 days. In the hypothesis that we are dealing with the same starspot, periodically occulted by the transiting planet, we estimated the projected orbital obliquity of WASP-52b to be lambda = 3.8 \pm 8.4 degree. We also determined the true orbital obliquity, psi = 20 \pm 50 degree, which is, although very uncertain, the first measurement of psi purely from starspot crossings. We finally assembled an optical transmission spectrum of the planet and searched for variations of its radius as a function of wavelength. Our analysis suggests a flat transmission spectrum within the experimental uncertainties., Comment: 16 pages, 12 figures, to appear in Monthly Notices of the Royal Astronomical Society

The WASP-98 planetary system represents a rare case of a hot Jupiter hosted by a metal-poor main-sequence star. We present a follow-up study of this system based on multi-band photometry and high-resolution spectroscopy. Two new transit events of WASP-98b were simultaneously observed in four passbands (g,r,i,z), using the telescope-defocussing technique, yielding eight high-precision light curves with point-to-point scatters of less than 1 mmag. We also collected three spectra of the parent star with a high-resolution spectrograph, which we used to remeasure its spectral characteristics, in particular its metallicity. We found this to be very low, Fe/H]=-0.49, but larger than was previously reported, [Fe/H]=-0.60. We used these new photometric and spectroscopic data to refine the orbital and physical properties of this planetary system, finding that the stellar and planetary mass measurements are significantly larger than those in the discovery paper. In addition, the multi-band light curves were used to construct an optical transmission spectrum of WASP-98b and probe the characteristics of its atmosphere at the terminator. We measured a lower radius at z compared with the other three passbands. The maximum variation is between the r and z bands, has a confidence level of roughly 6 sigma and equates to 5.5 pressure scale heights. We compared this spectrum to theoretical models, investigating several possible types of atmospheres, including hazy, cloudy, cloud-free, and clear atmospheres with titanium and vanadium oxide opacities. We could not find a good fit to the observations, except in the extreme case of a clear atmosphere with TiO and VO opacities, in which the condensation of Ti and V was suppressed. As this case is unrealistic, our results suggest the presence of an additional optical-absorbing species in the atmosphere of WASP-98b, of unknown chemical nature., Comment: 10 pages, 4 figures, to appear in MNRAS

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., Comment: 49 pages, 12 figures, accepted for publication in Space Science Reviews. Chapter in International Space Science Institute (ISSI) Book on "The Disk in Relation to the Formation of Planets and their Proto-atmospheres" to be published in Space Science Reviews by Springer

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., Comment: 11 pages, 9 figures, accepted for publication in MNRAS

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., Comment: 30 pages, 20 figures. Accepted for publication in ApJ

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., Comment: 30 pages, 16 figures. Accepted for publication in the International Journal of Astrobiology (Cambridge University Press)

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., Comment: Accepted to A&A. Identical as v1 except for corrected typos. 22 pages, 15 figures

Recent observations of the potentially habitable planets TRAPPIST-1 e, f, and g suggest that they possess large water mass fractions of possibly several tens of wt% of water, even though the host star’s activity should drive rapid atmospheric escape. These processes can photolyze water, generating free oxygen and possibly desiccating the planet. After the planets formed, their mantles were likely completely molten with volatiles dissolving and exsolving from the melt. In order to understand these planets and prepare for future observations, the magma ocean phase of these worlds must be understood. To simulate these planets, we have combined existing models of stellar evolution, atmospheric escape, tidal heating, radiogenic heating, magma ocean cooling, planetary radiation, and water-oxygen-iron geochemistry. We present Magm Oc, a versatile magma ocean evolution model, validated against the rocky Super-Earth GJ 1132b and early Earth. We simulate the coupled magma ocean-atmospheric evolution of TRAPPIST-1 e, f, and g for a range of tidal and radiogenic heating rates, as well as initial water contents between 1 and 100 Earth oceans. We also reanalyze the structures of these planets and find they have water mass fractions of 0–0.23, 0.01–0.21, and 0.11–0.24 for planets e, f, and g, respectively. Our model does not make a strong prediction about the water and oxygen content of the atmosphere of TRAPPIST-1 e at the time of mantle solidification. In contrast, the model predicts that TRAPPIST-1 f and g would have a thick steam atmosphere with a small amount of oxygen at that stage. For all planets that we investigated, we find that only 3 ́5% of the initial water will be locked in the mantle after the magma ocean solidified.