Search

[Abridged] The aim of this work is to search for an absorption signal from exospheric sodium (Na) and singly ionized calcium (Ca$^+$) in the optical transmission spectrum of the hot rocky super-Earth 55 Cancri e. Although the current best-fitting models to the planet mass and radius require a possible atmospheric component, uncertainties in the radius exist, making it possible that 55 Cancri e could be a hot rocky planet without an atmosphere. High resolution (R$\sim$110000) time-series spectra of five transits of 55 Cancri e, obtained with three different telescopes (UVES/VLT, HARPS/ESO 3.6m & HARPS-N/TNG) were analysed. Targeting the sodium D lines and the calcium H and K lines, the potential planet exospheric signal was filtered out from the much stronger stellar and telluric signals, making use of the change of the radial component of the orbital velocity of the planet over the transit from -57 to +57 km/sec. Combining all five transit data sets, we detect a signal potentially associated with sodium in the planet exosphere at a statistical significance level of 3$\sigma$. Combining the four HARPS transits that cover the calcium H and K lines, we also find a potential signal from ionized calcium (4.1 $\sigma$). Interestingly, this latter signal originates from just one of the transit measurements - with a 4.9$\sigma$ detection at this epoch. Unfortunately, due to the low significance of the measured sodium signal and the potentially variable Ca$^+$ signal, we estimate the p-values of these signals to be too high (corresponding to <4$\sigma$) to claim unambiguous exospheric detections. By comparing the observed signals with artificial signals injected early in the analysis, the absorption by Na and Ca$^+$ are estimated to be at a level of approximately 2.3$\times 10^{-3}$ and 7.0$\times 10^{-2}$ respectively, relative to the stellar spectrum., Comment: 15 pages, 8 figures, submission updated after English language editing, submission updated to correct a mistaken cross-reference noticed in A&A proof

Giant exoplanets orbiting very close to their parent star (hot Jupiters) are subject to tidal forces expected to synchronize their rotational and orbital periods on short timescales (tidal locking). However, spin rotation has never been measured directly for hot Jupiters. Furthermore, their atmospheres can show equatorial super-rotation via strong eastward jet streams, and/or high-altitude winds flowing from the day- to the night-side hemisphere. Planet rotation and atmospheric circulation broaden and distort the planet spectral lines to an extent that is detectable with measurements at high spectral resolution. We observed a transit of the hot Jupiter HD 189733 b around 2.3 {\mu}m and at a spectral resolution of R~10$^5$ with CRIRES at the ESO Very Large Telescope. After correcting for the stellar absorption lines and their distortion during transit (the Rossiter-McLaughlin effect), we detect the absorption of carbon monoxide and water vapor in the planet transmission spectrum by cross-correlating with model spectra. The signal is maximized (7.6{\sigma}) for a planet rotational velocity of $(3.4^{+1.3}_{-2.1})$ km/s, corresponding to a rotational period of $(1.7^{+2.9}_{-0.4})$ days. This is consistent with the planet orbital period of 2.2 days and therefore with tidal locking. We find that the rotation of HD 189733 b is longer than 1 day (3{\sigma}). The data only marginally (1.5{\sigma}) prefer models with rotation versus models without rotation. We measure a small day- to night-side wind speed of $(-1.7^{+1.1}_{-1.2})$ km/s. Compared to the recent detection of sodium blue-shifted by (8$\pm$2) km/s, this likely implies a strong vertical wind shear between the pressures probed by near-infrared and optical transmission spectroscopy., Comment: 24 pages, 8 figures, accepted for publication in ApJ

The spectral signature of an exoplanet can be separated from the spectrum of its host star using high-resolution spectroscopy. During such observations, the radial component of the planet's orbital velocity changes, resulting in a significant Doppler shift which allows its spectral features to be extracted. Aims: In this work, we aim to detect TiO in the optical transmission spectrum of HD 209458b. Gaseous TiO has been suggested as the cause of the thermal inversion layer invoked to explain the dayside spectrum of this planet. Method: We used archival data from the 8.2m Subaru Telescope taken with the High Dispersion Spectrograph of a transit of HD209458b in 2002. We created model transmission spectra which include absorption by TiO, and cross-correlated them with the residual spectral data after removal of the dominating stellar absorption features. We subsequently co-added the correlation signal in time, taking into account the change in Doppler shift due to the orbit of the planet. Results: We detect no significant cross-correlation signal due to TiO, though artificial injection of our template spectra into the data indicates a sensitivity down to a volume mixing ratio of ~10E-10. However, cross-correlating the template spectra with a HARPS spectrum of Barnard's star yields only a weak wavelength-dependent correlation, even though Barnard's star is an M4V dwarf which exhibits clear TiO absorption. We infer that the TiO line list poorly match the real positions of TiO lines at spectral resolutions of ~100,000. Similar line lists are also used in the PHOENIX and Kurucz stellar atmosphere suites and we show that their synthetic M-dwarf spectra also correlate poorly with the HARPS spectra of Barnard's star and five other M-dwarfs. We conclude that the lack of an accurate TiO line list is currently critically hampering this high-resolution retrieval technique., Comment: 11 pages, 10 figures. Accepted for publication in A&A

We present an analysis of the VIMS solar occultations dataset, which allows us to extract vertically resolved information on the characteristics of Titan's atmosphere between 100-700 km with a characteristic vertical resolution of 10 km. After a series of data treatment procedures, 4 occultations out of 10 are retained. This sample covers different seasons and latitudes of Titan. The transmittances show clearly the evolution of the haze and detect the detached layer at 310 km in Sept. 2011 at mid-northern latitudes. Through the inversion of the transmission spectra with a line-by-line radiative transfer code we retrieve the vertical distribution of CH$_4$ and CO mixing ratio. The two methane bands at 1.4 and 1.7 {\mu}m are always in good agreement and yield an average stratospheric abundance of $1.28\pm0.08$%. This is significantly less than the value of 1.48% obtained by the GCMS/Huygens instrument. The analysis of the residual spectra after the inversion shows that there are additional absorptions which affect a great part of the VIMS wavelength range. We attribute many of these additional bands to gaseous ethane, whose near-infrared spectrum is not well modeled yet. Ethane contributes significantly to the strong absorption between 3.2-3.5 {\mu}m that was previously attributed only to C-H stretching bands from aerosols. Ethane bands may affect the surface windows too, especially at 2.7 {\mu}m. Other residual bands are generated by stretching modes of C-H, C-C and C-N bonds. In addition to the C-H stretch from aliphatic hydrocarbons at 3.4 {\mu}m, we detect a strong and narrow absorption at 3.28 {\mu}m which we tentatively attribute to the presence of PAHs in the stratosphere. C-C and C-N stretching bands are possibly present between 4.3-4.5 {\mu}m. Finally, we obtain the CO mixing ratio between 70-170 km. The average result of $46\pm16$ ppm is in good agreement with previous studies., Comment: 51 pages, 28 figures

(Abridged) In recent years, ground-based high-resolution spectroscopy has become a powerful tool for investigating exoplanet atmospheres. It allows the robust identification of molecular species, and it can be applied to both transiting and non-transiting planets. Radial-velocity measurements of the star HD 179949 indicate the presence of a giant planet companion in a close-in orbit. Here we present the analysis of spectra of the system at 2.3 micron, obtained at a resolution of R~100,000, during three nights of observations with CRIRES at the VLT. We targeted the system while the exoplanet was near superior conjunction, aiming to detect the planet's thermal spectrum and the radial component of its orbital velocity. We detect molecular absorption from carbon monoxide and water vapor with a combined S/N of 6.3, at a projected planet orbital velocity of K_P = (142.8 +- 3.4) km/s, which translates into a planet mass of M_P = (0.98 +- 0.04) Jupiter masses, and an orbital inclination of i = (67.7 +- 4.3) degrees, using the known stellar radial velocity and stellar mass. The detection of absorption features rather than emission means that, despite being highly irradiated, HD 179949 b does not have an atmospheric temperature inversion in the probed range of pressures and temperatures. Since the host star is active (R_HK > -4.9), this is in line with the hypothesis that stellar activity damps the onset of thermal inversion layers owing to UV flux photo-dissociating high-altitude, optical absorbers. Finally, our analysis favors an oxygen-rich atmosphere for HD 179949 b, although a carbon-rich planet cannot be statistically ruled out based on these data alone., Comment: 10 pages, 9 figures. Accepted for publication in Astronomy and Astrophysics

We report a 4.8 sigma detection of water absorption features in the day side spectrum of the hot Jupiter HD 189733 b. We used high-resolution (R~100,000) spectra taken at 3.2 microns with CRIRES on the VLT to trace the radial-velocity shift of the water features in the planet's day side atmosphere during 5 h of its 2.2 d orbit as it approached secondary eclipse. Despite considerable telluric contamination in this wavelength regime, we detect the signal within our uncertainties at the expected combination of systemic velocity (Vsys=-3 +5-6 km/s) and planet orbital velocity (Kp=154 +14-10 km/s), and determine a H2O line contrast ratio of (1.3+/-0.2)x10^-3 with respect to the stellar continuum. We find no evidence of significant absorption or emission from other carbon-bearing molecules, such as methane, although we do note a marginal increase in the significance of our detection to 5.1 sigma with the inclusion of carbon dioxide in our template spectrum. This result demonstrates that ground-based, high-resolution spectroscopy is suited to finding not just simple molecules like CO, but also to more complex molecules like H2O even in highly telluric contaminated regions of the Earth's transmission spectrum. It is a powerful tool that can be used for conducting an immediate census of the carbon- and oxygen-bearing molecules in the atmospheres of giant planets, and will potentially allow the formation and migration history of these planets to be constrained by the measurement of their atmospheric C/O ratios., Comment: 5 pages, 4 figures, accepted for publication in MNRAS Letters

In this paper we present ground-based high-resolution spectroscopy of 51 Pegasi using CRIRES at the Very Large Telescope. The system was observed for 3x5 hours at 2.3 {\mu}m at a spectral resolution of R = 100,000, targeting potential signatures from carbon monoxide, water vapour and methane in the planet's dayside spectrum. In the first 2x5 hours of data, we find a combined signal from carbon monoxide and water in absorption at a formal 5.9{\sigma} confidence level, indicating a non-inverted atmosphere. We derive a planet mass of M_P = (0.46 +- 0.02) M_Jup and an orbital inclination i between 79.6 and 82.2 degrees, with the upper limit set by the non-detection of the planet transit in previous photometric monitoring. However, there is no trace of the signal in the final 5 hours of data. A statistical analysis indicates that the signal from the first two nights is robust, but we find no compelling explanation for its absence in the final night. The latter suffers from stronger noise residuals and greater instrumental instability than the first two nights, but these cannot fully account for the missing signal. It is possible that the integrated dayside emission from 51 Peg b is instead strongly affected by weather. However, more data are required before we can claim any time variability in the planet's atmosphere., Comment: 16 pages, 6 figures. Accepted for publication in ApJ

In recent years, day-side emission from about a dozen hot Jupiters has been detected through ground-based secondary eclipse observations in the near-infrared. These near-infrared observations are vital for determining the energy budgets of hot Jupiters, since they probe the planet's spectral energy distribution near its peak. The aim of this work is to measure the Ks-band secondary eclipse depth of WASP-33b, the first planet discovered to transit an A-type star. This planet receives the highest level of irradiation of all transiting planets discovered to date. Furthermore, its host-star shows pulsations and is classified as a low-amplitude delta-Scuti. As part of our GROUnd-based Secondary Eclipse (GROUSE) project we have obtained observations of two separate secondary eclipses of WASP-33b in the Ks-band using the LIRIS instrument on the William Herschel Telescope (WHT). The telescope was significantly defocused to avoid saturation of the detector for this bright star (K~7.5). To increase the stability and the cadence of the observations, they were performed in staring mode. We collected a total of 5100 and 6900 frames for the first and the second night respectively, both with an average cadence of 3.3 seconds. On the second night the eclipse is detected at the 12-sigma level, with a measured eclipse depth of 0.244+0.027-0.020 %. This eclipse depth corresponds to a brightness temperature of 3270+115-160 K. The measured brightness temperature on the second night is consistent with the expected equilibrium temperature for a planet with a very low albedo and a rapid re-radiation of the absorbed stellar light. For the other night the short out-of-eclipse baseline prevents good corrections for the stellar pulsations and systematic effects, which makes this dataset unreliable for eclipse depth measurements. This demonstrates the need of getting a sufficient out-of-eclipse baseline., Comment: 12 pages, 10 figures. Accepted for publication in Astronomy and Astrophysics

Context. The Kepler object KIC 12557548 b is peculiar. It exhibits transit-like features every 15.7 hours that vary in depth between 0.2% and 1.2%. Rappaport et al. (2012) explain the observations in terms of a disintegrating, rocky planet that has a trailing cloud of dust created and constantly replenished by thermal surface erosion. The variability of the transit depth is then a consequence of changes in the cloud optical depth. Aims. We aim to validate the disintegrating-planet scenario by modeling the detailed shape of the observed light curve, and thereby constrain the cloud particle properties to better understand the nature of this intriguing object. Methods. We analysed the six publicly-available quarters of raw Kepler data, phase-folded the light curve and fitted it to a model for the trailing dust cloud. Constraints on the particle properties were investigated with a light-scattering code. Results. The light curve exhibits clear signatures of light scattering and absorption by dust, including a brightening in flux just before ingress correlated with the transit depth and explained by forward scattering, and an asymmetry in the transit light curve shape, which is easily reproduced by an exponentially decaying distribution of optically thin dust, with a typical grain size of 0.1 micron. Conclusions. Our quantitative analysis supports the hypothesis that the transit signal of KIC 12557548 b is due to a variable cloud of dust, most likely originating from a disintegrating object., Comment: 5 pages, 4 figures. Accepted for publication in Astronomy and Astrophysics

[Abridged] The transmission of light through a planetary atmosphere can be studied as a function of altitude and wavelength using stellar or solar occultations, giving often unique constraints on the atmospheric composition. For exoplanets, a transit yields a limb-integrated, wavelength-dependent transmission spectrum of an atmosphere. When scattering haze and/or cloud particles are present in the planetary atmosphere, the amount of transmitted flux not only depends on the total optical thickness of the slant light path that is probed, but also on the amount of forward-scattering by the scattering particles. Here, we present results of calculations with a three-dimensional Monte Carlo code that simulates the transmitted flux during occultations or transits. For isotropically scattering particles, like gas molecules, the transmitted flux appears to be well-described by the total atmospheric optical thickness. Strongly forward-scattering particles, however, such as commonly found in atmospheres of Solar System planets, can increase the transmitted flux significantly. For exoplanets, such added flux can decrease the apparent radius of the planet by several scale heights, which is comparable to predicted and measured features in exoplanet transit spectra. We performed detailed calculations for Titan's atmosphere between 2.0 and 2.8 micron and show that haze and gas abundances will be underestimated by about 8% if forward-scattering is ignored in the retrievals. At shorter wavelengths, errors in the gas and haze abundances and in the spectral slope of the haze particles can be several tens of percent, also for other Solar System planetary atmospheres. We also find that the contribution of forward-scattering can be fairly well described by modelling the atmosphere as a plane-parallel slab., Comment: Icarus, accepted for publication

The giant planet orbiting tau Bootis was among the first extrasolar planets to be discovered through the reflex motion of its host star. It is one of the brightest known and most nearby planets with an orbital period of just a few days. Over the course of more than a decade, measurements of its orbital inclination have been announced and refuted, and have subsequently remained elusive until now. Here we report on the detection of carbon monoxide absorption in the thermal day-side spectrum of tau Bootis b. At a spectral resolution of R~100,000, we trace the change in the radial velocity of the planet over a large range in phase, determining an orbital inclination of i=44.5+-1.5 degrees and a true planet mass of 5.95+-0.28 MJup. This result extends atmospheric characterisation to non-transiting planets. The strong absorption signal points to an atmosphere with a temperature that is decreasing towards higher altitudes. This is a stark contrast to the temperature inversion invoked for other highly irradiated planets, and supports models in which the absorbing compounds believed to cause such atmospheric inversions are destroyed by the ultraviolet emission from the active host star., Comment: To appear in the June 28 issue of Nature: main article + supplementary information

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

GJ1214b is thought to be either a mini-Neptune with a thick, hydrogen-rich atmosphere, or a planet with a composition dominated by water. In the case of a hydrogen-rich atmosphere, molecular absorption and scattering processes may result in detectable radius variations as a function of wavelength. The aim of this paper is to measure these variations. We have obtained observations of the transit of GJ1214b in the r- and I-band with the INT, in the g, r, i and z bands with the 2.2 meter MPI/ESO telescope, in the Ks-band with the NOT, and in the Kc-band with the WHT. By comparing the transit depth between the the different bands, which is a measure for the planet-to-star size ratio, the atmosphere is investigated. We do not detect clearly significant variations in the planet-to-star size ratio as function of wavelength. Although the ratio at the shortest measured wavelength, in g-band, is 2sigma larger than in the other bands. The uncertainties in the Ks and Kc bands are large, due to systematic features in the light curves. The tentative increase in the planet-to-star size ratio at the shortest wavelength could be a sign of an increase in the effective planet-size due to Rayleigh scattering, which would require GJ1214b to have a hydrogen-rich atmosphere. If true, then the atmosphere has to have both clouds, to suppress planet-size variations at red optical wavelengths, as well as a sub-solar metallicity, to suppress strong molecular features in the near- and mid-infrared. However, star spots, which are known to be present on the hoststar's surface, can (partly) cancel out the expected variations in planet-to-star size ratio, due to the lower surface temperature of the spots . A hypothetical spot-fraction of 10% would be able to raise the infrared points sufficiently with respect to the optical measurements to be inconsistent with a water-dominated atmosphere. [abridged], Comment: 13 pages, 8 figures. Accepted for publication in A&A

Planets can emit polarized thermal radiation, just like brown dwarfs. We present calculated thermal polarization signals from hot exoplanets, using an advanced radiative transfer code that fully includes all orders of scattering by gaseous molecules and cloud particles. The code spatially resolves the disk of the planet, allowing simulations for horizontally inhomogeneous planets. Our results show that the degree of linear polarization, P, of an exoplanet's thermal radiation is expected to be highest near the planet's limb and that this P depends on the temperature and its gradient, the scattering properties and the distribution of the cloud particles. Integrated over the disk of a spherically symmetric planet, P of the thermal radiation equals zero. However, for planets that appear spherically asymmetric, e.g. due to flattening, cloud bands or spots in their atmosphere, differences in their day and night sides, and/or obscuring rings, P is often larger than 0.1 %, in favorable cases even reaching several percent at near-infrared wavelengths. Detection of thermal polarization signals can give access to planetary parameters that are otherwise hard to obtain: it immediately confirms the presence of clouds, and P can then constrain atmospheric inhomogeneities and the flattening due to the planet's rotation rate. For zonally symmetric planets, the angle of polarization will yield the components of the planet's spin axis normal to the line-of-sight. Finally, our simulations show that P is generally more sensitive to variability in a cloudy planet's atmosphere than the thermal flux is, and could hence better reveal certain dynamical processes., Comment: 9 pages, 5 figures, accepted for publication in ApJ

(abridged) We calculate near-infrared thermal emission spectra using a doubling-adding radiative transfer code, which includes scattering by clouds and haze. Initial temperature profiles and cloud optical depths are taken from the drift-phoenix brown dwarf model. As is well known, cloud particles change the spectrum compared to when clouds are ignored. The clouds reduce fluxes in the near-infrared spectrum and make it redder than for the clear sky case. We also confirm that not including scattering in the spectral calculations can result in errors on the spectra of many tens of percent, both in magnitude and in variations with wavelength. This is especially apparent for particles that are larger than the wavelength and only have little iron in them. Scattering particles will show deeper absorption features than absorbing (e.g. iron) particles and particle size will also affect the calculated infrared colours. Large particles also tend to be strongly forward-scattering, and we show that assuming isotropic scattering in this case also leads to very large errors in the spectrum. Thus, care must be taken in the choice of radiative transfer method for heat balance or spectral calculations when clouds are present in the atmosphere. Besides the choice of radiative transfer method, the type of particles that are predicted by models will change conclusions about e.g. infrared colours and trace gas abundances. As a result, knowledge of the scattering properties of the clouds is essential when deriving temperature profiles or gas abundances from direct infrared observations of exoplanets or brown dwarfs and from secondary eclipse measurements of transiting exoplanets, since scattering clouds will change the depth of gas absorption features, among other things. Thus, ignoring the presence of clouds can yield retrieved properties that differ significantly from the real atmospheric properties., Comment: Accepted for publication in Astronomy and Astrophysics. The abstract and a part of the introduction have been re-worded compared to the accepted version to avoid misinterpretation of the paper as much as possible

Context: Only recently it has become possible to measure the thermal emission from hot-Jupiters at near-Infrared wavelengths using ground-based telescopes, by secondary eclipse observations. This allows the planet flux to be probed around the peak of its spectral energy distribution, which is vital for the understanding of its energy budget. Aims: The aim of the reported work is to measure the eclipse depth of the planet HAT-P-1b at 2.2micron. This planet is an interesting case, since the amount of stellar irradiation it receives falls in between that of the two best studied systems (HD209458 and HD189733), and it has been suggested to have a weak thermal inversion layer. Methods: We have used the LIRIS instrument on the William Herschel Telescope (WHT) to observe the secondary eclipse of HATP-1b in the Ks-band, as part of our Ground-based secondary eclipse (GROUSE) project. The observations were done in staring mode, while significantly defocusing the telescope to avoid saturation on the K=8.4 star. With an average cadence of 2.5 seconds, we collected 6520 frames during one night. Results: The eclipse is detected at the 4sigma level, the measured depth being 0.109+/-0.025%. The uncertainties are dominated by residual systematic effects, as estimated from different reduction/analysis procedures. The measured depth corresponds to a brightness temperature of 2136+150-170K. This brightness temperature is significantly higher than those derived from longer wavelengths, making it difficult to fit all available data points with a plausible atmospheric model. However, it may be that we underestimate the true uncertainties of our measurements, since it is notoriously difficult to assign precise statistical significance to a result when systematic effects are important., Comment: 7 pages, 10 figures, Accepted for publication in A&A

A condensate cloud on Titan identified by its 220 cm−1 far-infrared signature continues to undergo seasonal changes at both the north and south poles. In the north, the cloud, which extends from 55 N to the pole, has been gradually decreasing in emission intensity since the beginning of the Cassini mission with a half-life of 3.8 years. The cloud in the south did not appear until 2012 but its intensity has increased rapidly, doubling every year. The shape of the cloud at the south pole is very different from that in the north. Mapping in 2013 December showed that the condensate emission was confined to a ring with a maximum at 80 S. The ring was centered 4deg from Titanʼs pole. The pattern of emission from stratospheric trace gases like nitriles and complex hydrocarbons (mapped in 2014 January) was also offset by 4deg, but had a central peak at the pole and a secondary maximum in a ring at about 70 S with a minimum at 80 S. The shape of the gas emission distribution can be explained by abundances that are high at the atmospheric pole and diminish toward the equator, combined with correspondingly increasing temperatures. We discuss possible causes for the condensate ring. The present rapid build up of the condensate cloud at the south pole is likely to transition to a gradual decline from 2015 to 2016. Key words: molecular processes - planets and satellites: atmospheres - planets and satellites: composition - planets and satellites: individual (Titan) - radiation mechanisms: thermal