CRyogenic High-resolution InfraRed Echelle Spectrograph Research Articles

The ESO SupJup Survey

Context. It has been proposed that the distinct formation and evolutionary pathways of exoplanets and brown dwarfs may affect the chemical and isotopic content of their atmospheres. Recent work has indeed shown differences in the 12C/13C isotope ratio, which have provisionally been attributed to the top-down formation of brown dwarfs and the core accretion pathway of super-Jupiters. Aims. The ESO SupJup Survey is aimed at disentangling the formation pathways of isolated brown dwarfs and planetary-mass companions using chemical and isotopic tracers. The survey utilises high-resolution spectroscopy with the recently upgraded CRyogenic high-resolution InfraRed Echelle Spectrograph (CRIRES+) at the Very Large Telescope, covering a total of 49 targets. Here, we present the first results of this survey: an atmospheric characterisation of DENIS J0255-4700, an isolated brown dwarf near the L-T transition. Methods. We analysed its observed CRIRES+ K-band spectrum using an atmospheric retrieval framework in which the radiative transfer code petitRADTRANS was coupled with the PyMultiNest sampling algorithm. Gaussian processes were employed to model inter-pixel correlations. In addition, we adopted an updated parameterisation of the pressure-temperature profile. Results. Abundances of CO, H2O, CH4, and NH3 were retrieved for this fast-rotating L-dwarf. The ExoMol H2O line list provides a significantly better fit than that of HITEMP. A free-chemistry retrieval is strongly favoured over equilibrium chemistry, caused by an under-abundance of CH4. The free-chemistry retrieval constrains a super-solar C/O-ratio of ~0.68 and a solar metallicity. We find tentative evidence (~3σ) for the presence of 13CO, with a constraint on the isotopologue ratio of 12CO/13CO = 184−40+61 and a lower limit of ≳97, which suggests a depletion of 13C compared to the local interstellar medium (12C/13C ~ 68). Conclusions. High-resolution, high signal-to-noise K-band spectra provide an excellent means of constraining the chemistry and isotopic content of sub-stellar objects, which is the main objective of the ESO SupJup Survey.

β Pictoris b through the eyes of the upgraded CRIRES+

Context. High-resolution spectrographs fed by adaptive optics (AO) provide a unique opportunity to characterize directly imaged exoplanets. Observations with such instruments allow us to probe the atmospheric composition, spin rotation, and radial velocity of the planet, thereby helping to reveal information on its formation and migration history. The recent upgrade of the Cryogenic High-Resolution Infrared Echelle Spectrograph (CRIRES+) at the VLT makes it a highly suitable instrument for characterizing directly imaged exoplanets. Aims. In this work, we report on observations of β Pictoris b with CRIRES+ and use them to constrain the planets atmospheric properties and update the estimation of its spin rotation. Methods. The data were reduced using the open-source pycrires package. We subsequently forward-modeled the stellar, planetary, and systematic contribution to the data to detect molecules in the planet’s atmosphere. We also used atmospheric retrievals to provide new constraints on its atmosphere. Results. We confidently detected water and carbon monoxide in the atmosphere of β Pictoris b and retrieved a slightly sub-solar carbon-to-oxygen ratio, which is in agreement with previous results. The interpretation is hampered by our limited knowledge of the C/O ratio of the host star. We also obtained a much improved constraint on its spin rotation of 19.9 ± 1.0 kms−1, which gives a rotation period of 8.7 ± 0.8 h, assuming no obliquity. We find that there is a degeneracy between the metallicity and clouds, but this has minimal impact on the retrieved C/O, υ sin i, and radial velocity. Our results show that CRIRES+ is performing well and stands as a highly useful instrument for characterizing directly imaged planets.

A high-resolution study of near-IR diffuse interstellar bands, search for small-scale structure, time variability, and stellar features

Context. The diffuse interstellar bands (DIBs) are a set of hundreds of unidentified absorption features that appear almost ubiquitously throughout the interstellar medium. Most DIBs appear at optical wavelengths, but some are in the near-infrared. Aims. We aim to characterise near-infrared DIBs at high spectral resolving power towards multiple targets. Methods. We observed 76 early-type stars at a resolving power of 50 000 (velocity resolution ~6 km s−1) and signal-to-noise ratios of several hundreds using the CRyogenic high-resolution InfraRed Echelle Spectrograph (CRIRES). These data allow us to investigate the DIBs around 1318.1, 1527.4, 1561.1, 1565.1, 1567.0, 1574.4, and 1624.2 nm. We detected a total of six DIB features and 17 likely stellar features through comparisons with a model spectrum computed with CMFGEN. Additionally, we measured equivalent widths of the DIBs at 1318.1 and 1527.4 nm using observations with X-shooter towards ten very highly reddened (3.2 < E(B–V) < 6.5) Cepheid variable stars and towards four stars observed at low values of precipitable water vapour as well as by using other archive data. Results. We measured correlations (correlation coefficient r ~ 0.73–0.96, depending on the subsample used) between DIB equivalent width and reddening for the DIBs at 1318.1, 1561.1, 1565.1, and 1567.0 nm. Comparing the near-infrared DIBs with 50 of the strongest optical DIBs, we find correlations r > ~0.8 between the 1318, 1527, 1561, 1565, and 1567 nm and the optical DIBs 5705, 5780, 6203, 6283, and 6269 Å. The 5797 Å DIB is less well correlated with the near-infrared DIBs. The DIB at 9632.1 Å, which is likely C60+, is not well correlated with the 1318.1 nm DIB. Partial correlation coefficients using E(B-V) as the covariate were also determined. For stars earlier than B2, the 1318.1 nm DIB is affected by an emission line on its blue wing that is likely stellar in nature, although we cannot rule out an interstellar or circumstellar origin for this line caused by, for example, a DIB in emission. The 1318.1 nm DIB also has an extended red wing. The line is reasonably well fitted by two Gaussian components, although neither the component equivalent width (EW) ratios nor the separation between components are obviously correlated with such indicators as λλ5780/5797 and reddening. The EW at 1318 nm correlates with H I with EW(1318 nm)/E(B – V) decreasing with f(H2). Five pairs of stars within one arcmin of each other show very similar 1318.1 nm DIB profiles. Possible variation in the 1318.1 nm feature is seen between HD 145501 and HD 145502 (separated by 41 arcsec, equivalent to 7200 au) and HD 168607 and HD 168625 (separated by 67 arcsec, equivalent to 0.52 pc on the plane of the sky). Seventeen sightlines have repeat CRIRES observations separated by six to 14 months, and two sightlines have repeat X-shooter observations separated by 9.9 yr. No time variability was detected in the 1318.1 nm DIB in the CRIRES data nor in the 5780.5 Å, 5797.1 Å, 1318.1 nm, and 1527.4 nm DIBs. Tentative time variation is observed in the C60+ DIBs at 9577 and 9632 Å towards HD 183143, although it is very close to the noise level and requires confirmation. Conclusions. The Near Infrared (NIR) DIBs observed occur more in more UV-irradiated regions than the 5797 Å DIB allowing the study of heavily reddened sightlines. Future searches for time variability in DIBs will require either higher quality data, larger intervals between epochs, or both.

CRIRES high-resolution near-infrared spectroscopy of diffuse interstellar band profiles

Aims. A high spectral resolution investigation of diffuse interstellar bands (DIBs) in the near-infrared (YJ band) is conducted to test new methods, to confirm and improve existing parameters, and to search for new DIBs. Methods. The CRyogenic high-resolution InfraRed Echelle Spectrograph (CRIRES) on the European Southern Observatory's Very Large Telescope was employed to obtain spectra of four reddened background supergiant stars (HD 183143, HD 165784, HD 92207, HD 111613) and an unreddened comparison star (HD 87737) at the highest resolution of R ≈ 100 000 currently achievable at near-infrared wavelengths, more than twice as high as accomplished in previous near-infrared DIB studies. The correction for telluric absorption was performed by a modelling approach. Non-local thermodynamic equilibrium spectral modelling of available optical and the new near-infrared stellar spectra facilitated a comprehensive characterisation of the atmospheric properties of the background stars. As a consequence, a more precise and accurate determination of the reddening and the reddening law along the respective sight lines could be achieved than feasible before by comparison of the observed and model spectral energy distributions. For DIBs that overlap with stellar lines the DIB profile shapes could be recovered. Results. Seventeen known near-infrared DIBs were confirmed, and 12 previously unknown and generally weaker DIBs were identified in the YJ band. Three DIBs that show uniform profiles along all sight lines were identified, possibly connected to transitions from a common lower state of the same carrier. The divergent extinction curve towards the frequently discussed DIB standard star HD 183143 could be reproduced for the first time, requiring extra absorption by ~3.5 mag due to polycyclic aromatic hydrocarbons (PAHs) to match the ultraviolet extinction bump. This extra absorption probably stems from a circumstellar bubble lying in front of the star which is intersected tangentially by the line of sight, making this particular sight line more peculiar than standard.

Mineral snowflakes on exoplanets and brown dwarfs

Context. Exoplanet atmosphere characterisation has become an important tool in understanding exoplanet formation, evolution, and it also is a window into potential habitability. However, clouds remain a key challenge for characterisation: upcoming space telescopes (e.g. the James Webb Space Telescope, JWST, and the Atmospheric Remote-sensing Infrared Exoplanet Large-survey) and ground-based high-resolution spectrographs (e.g. the next-generation CRyogenic high-resolution InfraRed Echelle Spectrograph) will produce data requiring detailed understanding of cloud formation and cloud effects for a variety of exoplanets and brown dwarfs. Aims. We aim to understand how the micro-porosity of cloud particles affects the cloud structure, particle size, and material composition on exoplanets and brown dwarfs. We further examine the spectroscopic effects of micro-porous particles, the particle size distribution, and non-spherical cloud particles. Methods. We expanded our kinetic non-equilibrium cloud formation model to study the effect of micro-porosity on the cloud structure using prescribed 1D (Tgas–pgas) profiles from the DRIFT-PHOENIX model atmosphere grid. We applied the effective medium theory and the Mie theory to model the spectroscopic properties of cloud particles with micro-porosity and a derived particle size distribution. In addition, we used a statistical distribution of hollow spheres to represent the effects of non-spherical cloud particles. Results. Highly micro-porous cloud particles (90% vacuum) have a larger surface area, enabling efficient bulk growth higher in the atmosphere than for compact particles. Increases in single scattering albedo and cross-sectional area for these mineral snowflakes cause the cloud deck to become optically thin only at a wavelength of ~100 μm instead of at the ~20 μm for compact cloud particles. A significant enhancement in albedo is also seen when cloud particles occur with a locally changing Gaussian size distribution. Non-spherical particles increase the opacity of silicate spectral features, which further increases the wavelength at which the clouds become optically thin. Conclusions. Retrievals of cloud properties, particularly particle size and mass of clouds, are biased by the assumption of compact spherical particles. The JWST mid-infrared instrument will be sensitive to signatures of micro-porous and non-spherical cloud particles based on the wavelength at which clouds are optically thin. Details of spectral features are also dependent on particle shape, and greater care must be taken in modelling clouds as observational data improves.

The full infrared spectrum of molecular hydrogen

Context. The high spectral resolution R ∼ 45 000 provided by IGRINS (Immersion Grating INfrared Spectrometer) at MacDonald Observatory and R ∼ 100 000 achieved by CRIRES (CRyogenic high-resolution InfraRed Echelle Spectrograph) at VLT (Very Large Telescope) challenges the present knowledge of infrared spectra. Aims. We aim to predict the full infrared spectrum of molecular hydrogen at a comparable accuracy. Methods. We take advantage of the recent theoretical ab initio studies on molecular hydrogen to compute both the electric quadrupole and magnetic dipole transitions taking place within the ground electronic molecular state of hydrogen. Results. We computed the full infrared spectrum of molecular hydrogen at an unprecedented accuracy and derive for the first time the emission probabilities including both electric quadrupole (ΔJ = 0, ±2) and magnetic dipole transitions (ΔJ = 0) as well as the total radiative lifetime of each rovibrational state. Inclusion of magnetic dipole transitions increases the emission probabilities by factors of a few for highly excited rotational levels, which occur in the 3–20 μ range.

Evidence for Multiple Molecular Species in the Hot Jupiter HD 209458b

Abstract Molecular species in planetary atmospheres provide key insights into their atmospheric processes and formation conditions. In recent years, high-resolution Doppler spectroscopy in the near-infrared has allowed detections of H2O and CO in the atmospheres of several hot Jupiters. This method involves monitoring the spectral lines of the planetary thermal emission Doppler-shifted due to the radial velocity of the planet over its orbit. However, aside from CO and H2O, which are the primary oxygen- and carbon-bearing species in hot H2-rich atmospheres, little else is known about the molecular compositions of hot Jupiters. Several recent studies have suggested the importance and detectability of nitrogen-bearing species in such atmospheres. In this Letter, we confirm potential detections of CO and H2O in the hot Jupiter HD 209458b using high-resolution spectroscopy. We also report a cross-correlation peak with a signal-to-noise ratio of 4.7 from a search for HCN. The results are obtained using high-resolution phase-resolved spectroscopy with the Very Large Telescope CRyogenic high-resolution InfraRed Echelle Spectrograph (VLT CRIRES) and standard analysis methods reported in the literature. A more robust treatment of telluric contamination and other residuals would improve confidence and enable unambiguous molecular detections. The presence of HCN could provide constraints on the C/O ratio of HD 209458b and its potential origins.

A proposed new diagnostic for Herbig disc geometry

Aims. CO ro-vibrational lines observed from Herbig group II discs are seen to be often broad while the same lines observed from group I discs are often narrow. This difference is not well understood. In this paper we explore the underlying cause for this difference and provide a pathway for a better understanding of the geometry and structure of the inner discs around Herbig Ae/Be stars. Methods. High-spectral-resolution infrared spectra of CO ro-vibrational emission from six Herbig Ae/Be candidate stars were taken with CRIRES (CRyogenic high-resolution InfraRed Echelle Spectrograph) at the VLT (Very Large Telescope). From these spectra, we produce individual and co-added CO ro-vibrational line profiles. We investigate line profile shape differences, and we explore the FWHM (Full Width Half Maximum) variations with J quantum number in the context of disc geometry. Furthermore, we put our new sources into the context of earlier observed sources to study a large sample. For comparison, we also investigate the FWHM variations with J of modelled CO ro-vibrational lines from two typical disc geometries produced with the thermo-chemical disc modelling code ProDiMo. Results. - abbreviated - Comparing dust and gas inner disc geometries (inferred by SED and CO ro-vibrational emission) for the expanded sample of observed Herbig discs, we find no clear correspondence between the SED groups of the sources and their inner CO radius. Conclusions. The FWHM versus J is a potential new gas diagnostic for the inner disc with e.g. a constant FWHM versus J indicating the presence of a large gas hole or gap. Both models and observations indicate the potential of this new diagnostic. - abbreviated -

CRIRES-VLT high-resolution spectro-astrometry as a tool in the search for disks inside the cores of planetary nebulae

The onset of the asymmetry in planetary nebulae (PNe) occurs during the short transition between the end of the asymptotic giant branch (AGB) phase and the beginning of the PN phase. Sources in this transition phase are compact and emit intensely in infrared wavelengths, making high spatial resolution observations in the infrared mandatory to investigate the shaping process of PNe. Interferometric VLTI IR observations have revealed compelling evidence of disks at the cores of PNe, but the limited sensitivity, strong observational constraints, and limited spatial coverage place severe limits on the universal use of this technique. Inspired by the successful detection of proto-planetary disks using spectro-astrometric observations, we apply here for the first time this technique to search for sub-arcsecond structures in PNe. Our exploratory study using CRIRES (CRyogenic high-resolution Infra-Red Echelle Spectrograph) commissioning data of the proto-PN IRAS 17516-2525 and the young PN SwSt 1 has revealed small-sized structures after the spectro-astrometric analysis of the two sources. In IRAS 17516-2525, the spectro-astrometric signal has a size of only 12 mas, as detected in the Brackett-gamma line, whereas the structures found in SwSt 1 have sizes of 230 mas in the [Fe III] line and 130 mas in the Brackett-gamma line. The spectroscopic observations required to perform spectro-astrometry of sources in the transition towards the PN phase are less time consuming and much more sensitive than VLTI IR observations. The results presented here open a new window in the search of the small-sized collimating agents that shape the complex morphologies of extremely axisymmetric PNe.

CO ro-vibrational lines in HD 100546

Aims. We have studied the emission of CO ro-vibrational lines in the disc around the Herbig Be star HD 100546 to determine physical properties, disc asymmetries, the CO excitation mechanism, and the spatial extent of the emission, with the final goal of using the CO ro-vibrational lines as a diagnostic to understand inner disc structure in the context of planet formation. Methods. High-spectral-resolution infrared spectra of CO ro-vibrational emission at eight different position angles were taken with the CRyogenic high-resolution InfraRed Echelle Spectrograph (CRIRES) at the Very Large Telescope (VLT). From these spectra flux tables, line profiles for individual CO ro-vibrational transitions, co-added line profiles, and population diagrams were produced. We have investigated variations in the line profile shapes and line strengths as a function of slit position angle. We used the thermochemical disc modelling code ProDiMo based on the chemistry, radiation field, and temperature structure of a previously published model for HD 100546. We calculated line fluxes and profiles for the whole set of observed CO ro-vibrational transitions using a large CO model molecule that includes the lowest two electronic states, each with 7 vibrational levels and within them 60 rotational levels. Comparing observations and the model, we investigated the possibility of disc asymmetries, the excitation mechanism (UV fluorescence), the geometry, and physical conditions of the inner disc. Results. The observed CO ro-vibrational lines are largely emitted from the inner rim of the outer disc at 10–13 AU. The line shapes are similar for all v levels and line fluxes from all vibrational levels vary only within one order of magnitude. All line profile asymmetries and variations can be explained with a symmetric disc model to which a slit correction and pointing offset is applied. Because the angular size of the CO emitting region (10–13 AU) and the slit width are comparable the line profiles are very sensitive to the placing of the slit. The model reproduces the line shapes and the fluxes of the v = 1–0 lines as well as the spatial extent of the CO rovibrational emission. It does not reproduce the observed band ratios of 0.5–0.2 with higher vibrational bands. We find that lower gas volume densities at the surface of the inner rim of the outer disc can make the fluorescence pumping more efficient and reproduce the observed band ratios.

Identifying new opportunities for exoplanet characterisation at high spectral resolution

[Abridged] Recently, there have been a series of detections of molecules in the atmospheres of extrasolar planets using high spectral resolution (R~100,000) observations, mostly using the CRyogenic high-resolution InfraRed Echelle Spectrograph (CRIRES) on the Very Large Telescope. These measurements are able to resolve molecular bands into individual absorption lines. Observing many lines simultaneously as their Doppler shift changes with time allows the detection of specific molecules in the atmosphere of the exoplanet. We performed simulations of high-resolution CRIRES observations of a planet's thermal emission and transit between 1-5 micron and performed a cross-correlation analysis on these results to assess how well the planet signal can be extracted. We also simulated day-side and night-side spectra at high spectral resolution for planets with and without a day-side temperature inversion, based on the cases of HD 189733b and HD 209458b. Several small wavelength regions in the L-band promise to yield cross-correlation signals from the thermal emission of hot Jupiters that can exceed those of the current detections by up to a factor of 2-3 for the same integration time. For transit observations, the H-band is also attractive, with the H, K, and L-band giving cross-correlation signals of similar strength. High-resolution night-side spectra of hot Jupiters can give cross-correlation signals as high as the day-side, or even higher. We show that there are many new possibilities for high-resolution observations of exoplanet atmospheres that have expected planet signals at least as high as those already detected. Hence, high-resolution observations at well-chosen wavelengths and at different phases can improve our knowledge about hot Jupiter atmospheres significantly, already with currently available instrumentation.

Exoplanet atmospheres at high spectral resolution: A CRIRES survey of hot-Jupiters

AbstractRecently, we presented the detection of carbon monoxide in the transmission spectrum of extrasolar planet HD209458b, using CRIRES, the Cryogenic high-resolution Infrared Echelle Spectrograph at ESO's Very Large Telescope (VLT). The high spectral resolution observations (R=100,000) provide a wealth of information on the planet's orbit, mass, composition, and even on its atmospheric dynamics. The new observational strategy and data analysis techniques open up a whole world of opportunities. We therefore started an ESO large program using CRIRES to explore these, targeting both transiting and non-transiting planets in carbon monoxide, water vapour, and methane. Observations of the latter molecule will also serve as a test-bed for METIS, the proposed mid-infrared imager and spectrograph for the European Extremely Large Telescope.