This Letter reports on the photometric detection of transits of the Neptune-mass planet orbiting the nearby M-dwarf star GJ 436. It is by far the closest, smallest, and least massive transiting ...planet detected so far. Its mass is slightly larger than Neptune's at $M = 22.6 \pm 1.9~M_\oplus $. The shape and depth of the transit lightcurves show that it is crossing the host star disc near its limb (impact parameter $0.84 \pm 0.03$) and that the planet size is comparable to that of Uranus and Neptune, $R = 25\,200 \pm 2200$ km = $3.95 \pm 0.35~R_\oplus$. Its main constituant is therefore very likely to be water ice. If the current planet structure models are correct, an outer layer of H/He constituting up to ten percent in mass is probably needed on top of the ice to account for the observed radius.
Context.
With more than 1000 h of observation from Feb. 2016 to Oct. 2019, the
Spitzer
Exploration Program Red Worlds (ID: 13067, 13175 and 14223) exclusively targeted TRAPPIST-1, a nearby (12 pc) ...ultracool dwarf star, finding that it is orbited by seven transiting Earth-sized planets. At least three of these planets orbit within the classical habitable zone of the star, and all of them are well-suited for a detailed atmospheric characterization with the upcoming JWST.
Aims.
The main goals of the
Spitzer
Red Worlds program were (1) to explore the system for new transiting planets, (2) to intensively monitor the planets’ transits to yield the strongest possible constraints on their masses, sizes, compositions, and dynamics, and (3) to assess the infrared variability of the host star. In this paper, we present the global results of the project.
Methods.
We analyzed 88 new transits and combined them with 100 previously analyzed transits, for a total of 188 transits observed at 3.6 or 4.5
μ
m. For a comprehensive study, we analyzed all light curves both individually and globally. We also analyzed 29 occultations (secondary eclipses) of planet b and eight occultations of planet c observed at 4.5
μ
m to constrain the brightness temperatures of their daysides.
Results.
We identify several orphan transit-like structures in our
Spitzer
photometry, but all of them are of low significance. We do not confirm any new transiting planets. We do not detect any significant variation of the transit depths of the planets throughout the different campaigns. Comparing our individual and global analyses of the transits, we estimate for TRAPPIST-1 transit depth measurements mean noise floors of ~35 and 25 ppm in channels 1 and 2 of
Spitzer
/IRAC, respectively. We estimate that most of this noise floor is of instrumental origins and due to the large inter-pixel inhomogeneity of IRAC InSb arrays, and that the much better interpixel homogeneity of JWST instruments should result in noise floors as low as 10 ppm, which is low enough to enable the atmospheric characterization of the planets by transit transmission spectroscopy. Our analysis reveals a few outlier transits, but we cannot conclude whether or not they correspond to spot or faculae crossing events. We construct updated broadband transmission spectra for all seven planets which show consistent transit depths between the two
Spitzer
channels. Although we are limited by instrumental precision, the combined transmission spectrum of planet b to g tells us that their atmospheres seem unlikely to be CH
4
-dominated. We identify and model five distinct high energy flares in the whole dataset, and discuss our results in the context of habitability. Finally, we fail to detect occultation signals of planets b and c at 4.5
μ
m, and can only set 3-
σ
upper limits on their dayside brightness temperatures (611 K for b 586 K for c).
Abstract
The TRAPPIST-1 planetary system is a favourable target for the atmospheric characterization of temperate earth-sized exoplanets by means of transmission spectroscopy with the forthcoming ...James Webb Space Telescope (JWST). A possible obstacle to this technique could come from the photospheric heterogeneity of the host star that could affect planetary signatures in the transit transmission spectra. To constrain further this possibility, we gathered an extensive photometric data set of 25 TRAPPIST-1 transits observed in the near-IR J band (1.2 μm) with the UKIRT and the AAT, and in the NB2090 band (2.1 μm) with the VLT during the period 2015–18. In our analysis of these data, we used a special strategy aiming to ensure uniformity in our measurements and robustness in our conclusions. We reach a photometric precision of 0.003 (RMS of the residuals), and we detect no significant temporal variations of transit depths of TRAPPIST-1 b, c, e, and g over the period of 3 yr. The few transit depths measured for planets d and f hint towards some level of variability, but more measurements will be required for confirmation. Our depth measurements for planets b and c disagree with the stellar contamination spectra originating from the possible existence of bright spots of temperature 4500 K. We report updated transmission spectra for the six inner planets of the system which are globally flat for planets b and g and some structures are seen for planets c, d, e, and f.
Transit timing variations (TTVs) can be a very efficient way of constraining masses and eccentricities of multi-planet systems. Recent measurements of the TTVs of TRAPPIST-1 have led to an estimate ...of the masses of the planets, enabling an estimate of their densities and their water content. A recent TTV analysis using data obtained in the past two years yields a 34 and 13% increase in mass for TRAPPIST-1b and c, respectively. In most studies to date, a Newtonian
N
-body model is used to fit the masses of the planets, while sometimes general relativity is accounted for. Using the Posidonius
N
-body code, in this paper we show that in the case of the TRAPPIST-1 system, non-Newtonian effects might also be relevant to correctly model the dynamics of the system and the resulting TTVs. In particular, using standard values of the tidal Love number
k
2
(accounting for the tidal deformation) and the fluid Love number
k
2
f
(accounting for the rotational flattening) leads to differences in the TTVs of TRAPPIST-1b and c that are similar to the differences caused by general relativity. We also show that relaxing the values of tidal Love number
k
2
and the fluid Love number
k
2
f
can lead to TTVs which differ by as much as a few 10 s on a 3−4-yr timescale, which is a potentially observable level. The high values of the Love numbers needed to reach observable levels for the TTVs could be achieved for planets with a liquid ocean, which if detected might then be interpreted as a sign that TRAPPIST-1b and TRAPPIST-1c could have a liquid magma ocean. For TRAPPIST-1 and similar systems the models to fit the TTVs should potentially account for general relativity, for the tidal deformation of the planets, for the rotational deformation of the planets, and to a lesser extent for the rotational deformation of the star, which would add up to 7 × 2 + 1 = 15 additional free parameters in the case of TRAPPIST-1.
High-precision eclipse spectrophotometry of transiting terrestrial exoplanets represents a promising path for the first atmospheric characterizations of habitable worlds and the search for life ...outside our solar system. The detection of terrestrial planets transiting nearby late-type M-dwarfs could make this approach applicable within the next decade, with soon-to-come general facilities. In this context, we previously identified GJ 1214 as a high-priority target for a transit search, as the transit probability of a habitable planet orbiting this nearby M4.5 dwarf would be significantly enhanced by the transiting nature of GJ 1214 b, the super-Earth already known to orbit the star. Based on this observation, we have set up an ambitious high-precision photometric monitoring of GJ 1214 with the Spitzer Space Telescope to probe the inner part of its habitable zone in search of a transiting planet as small as Mars. We present here the results of this transit search. Unfortunately, we did not detect any other transiting planets. Assuming that GJ 1214 hosts a habitable planet larger than Mars that has an orbital period smaller than 20.9 days, our global analysis of the whole Spitzer dataset leads to an a posteriori no-transit probability of ~98%. Our analysis allows us to significantly improve the characterization of GJ 1214 b, to measure its occultation depth to be 70 ± 35 ppm at 4.5 μm, and to constrain it to be smaller than 205 ppm (3σ upper limit) at 3.6 μm. In agreement with the many transmission measurements published so far for GJ 1214 b, these emission measurements are consistent with both a metal-rich and a cloudy hydrogen-rich atmosphere.
We report the detection of the secondary eclipse of the hot Jupiter HD 209458 b in optical/visible light using the CHEOPS space telescope. Our measurement of 20.4
−3.3
+3.2
parts per million ...translates into a geometric albedo of
A
g
= 0.096 ± 0.016. The previously estimated dayside temperature of about 1500 K implies that our geometric albedo measurement consists predominantly of reflected starlight and is largely uncontaminated by thermal emission. This makes the present result one of the most robust measurements of
A
g
for any exoplanet. Our calculations of the bandpass-integrated geometric albedo demonstrate that the measured value of
A
g
is consistent with a cloud-free atmosphere, where starlight is reflected via Rayleigh scattering by hydrogen molecules, and the water and sodium abundances are consistent with stellar metallicity. We predict that the bandpass-integrated TESS geometric albedo is too faint to detect and that a phase curve of HD 209458 b observed by CHEOPS would have a distinct shape associated with Rayleigh scattering if the atmosphere is indeed cloud free.
Context.
In the search for small exoplanets orbiting cool stars whose spectral energy distributions peak in the near infrared, the strong absorption of radiation in this region due to water vapour in ...the atmosphere is a particularly adverse effect for ground-based observations of cool stars.
Aims.
To achieve the photometric precision required to detect exoplanets in the near infrared, it is necessary to mitigate the effect of variable precipitable water vapour (PWV) on radial-velocity and photometric measurements. The aim is to enable global PWV correction by monitoring the amount of precipitable water vapour at zenith and along the line of sight of any visible target.
Methods.
We developed an open-source Python package that uses imagery data obtained with the Geostationary Operational Environmental Satellites (GOES), which provide temperature and relative humidity at different pressure levels to compute the near real-time PWV above any ground-based observatory that is covered by GOES every 5 min or 10 min, depending on the location.
Results.
We computed PWV values on selected days above Cerro Paranal (Chile) and San Pedro Mártir (Mexico) to benchmark the procedure. We also simulated different pointing at test targets as observed from the sites to compute the PWV along the line of sight. To asses the accuracy of our method, we compared our results with the on-site radiometer measurements obtained from Cerro Paranal.
Conclusions.
Our results show that our publicly available code proves to be a good supporting tool for measuring the local PWV for any ground-based facility within the GOES coverage, which will help to reduce correlated noise contributions in near-infrared ground-based observations that do not benefit from on-site PWV measurements.
We measured the radii of 7 low and very low-mass stars using long baseline interferometry with the VLTI interferometer and its VINCI and AMBER near-infrared recombiners. We use these new data, ...together with literature measurements, to examine the luminosity-radius and mass-radius relations for K and M dwarfs. The precision of the new interferometric radii now competes with what can be obtained for double-lined eclipsing binaries. Interferometry provides access to much less active stars, as well as to stars with much better measured distances and luminosities, and therefore complements the information obtained from eclipsing systems. The radii of magnetically quiet late-K to M dwarfs match the predictions of stellar evolution models very well, providing direct confirmation that magnetic activity explains the discrepancy that was recently found for magnetically active eclipsing systems. The radii of the early K dwarfs are reproduced well for a mixing length parameter that approaches the solar value, as qualitatively expected.
Context.
Measurements of the occultation of an exoplanet at visible wavelengths allow us to determine the reflective properties of a planetary atmosphere. The observed occultation depth can be ...translated into a geometric albedo. This in turn aids in characterising the structure and composition of an atmosphere by providing additional information on the wavelength-dependent reflective qualities of the aerosols in the atmosphere.
Aims.
Our aim is to provide a precise measurement of the geometric albedo of the gas giant HD 189733b by measuring the occultation depth in the broad optical bandpass of CHEOPS (350–1100 nm).
Methods.
We analysed 13 observations of the occultation of HD 189733b performed by CHEOPS utilising the Python package PyCHEOPS. The resulting occultation depth is then used to infer the geometric albedo accounting for the contribution of thermal emission from the planet. We also aid the analysis by refining the transit parameters combining observations made by the TESS and CHEOPS space telescopes.
Results.
We report the detection of an 24.7 ± 4.5 ppm occultation in the CHEOPS observations. This occultation depth corresponds to a geometric albedo of 0.076 ± 0.016. Our measurement is consistent with models assuming the atmosphere of the planet to be cloud-free at the scattering level and absorption in the CHEOPS band to be dominated by the resonant Na doublet. Taking into account previous optical-light occultation observations obtained with the
Hubble
Space Telescope, both measurements combined are consistent with a super-stellar Na elemental abundance in the dayside atmosphere of HD 189733b. We further constrain the planetary Bond albedo to between 0.013 and 0.42 at 3
σ
confidence.
Conclusions.
We find that the reflective properties of the HD 189733b dayside atmosphere are consistent with a cloud-free atmosphere having a super-stellar metal content. When compared to an analogous CHEOPS measurement for HD 209458b, our data hint at a slightly lower geometric albedo for HD 189733b (0.076 ± 0.016) than for HD 209458b (0.096 ± 0.016), or a higher atmospheric Na content in the same modelling framework. While our constraint on the Bond albedo is consistent with previously published values, we note that the higher-end values of ~0.4, as derived previously from infrared phase curves, would also require peculiarly high reflectance in the infrared, which again would make it more difficult to disentangle reflected and emitted light in the total observed flux, and therefore to correctly account for reflected light in the interpretation of those phase curves. Lower reported values for the Bond albedos are less affected by this ambiguity.