A search for TiO in the optical high-resolution transmission spectrum of HD 209458b: Hindrance due to inaccuracies in the line database
Key Words.:Line: Identification – Molecular data – Planets and satellites: Atmospheres – Planets and satellites: HD 209458b – Methods: Observational – Methods: Spectroscopic
Context: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.
Methods: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 . 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 . 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.
1.1 High Dispersion Spectroscopy
Recently, carbon monoxide was identified in the transmission spectrum of the exoplanet HD209458b (Snellen et al. 2010). During transit, the radial velocity of the hot-Jupiter changes by a few tens of , resulting in a changing Doppler shift which is measurable at a resolution of . This time-variation of the planet’s spectral signature provides a powerful tool to discriminate between the planet and the overwhelmingly bright host star, because the absorption features in the stellar spectrum are quasi-invariant on these time scales (see Figure 1). Futhermore, because the orbital period of the planet is well known, its radial velocity is known at all times and the Doppler shift of the planets’ absorption features can be targeted specifically.
This technique has also been successfully applied during other parts of of the orbit, where direct thermal emission of the planet is probed. Absorption in the dayside spectra of both transiting and non-transiting hot-Jupiters has been observed: CO and absorption in the atmospheres of Boötes b, HD 189733 b, 51 Peg b and HD 179733 b (Birkby et al. 2013; Brogi et al. 2012, 2013, 2014; de Kok et al. 2013; Lockwood et al. 2014; Rodler et al. 2013). In case of the non-transiting planets, this has allowed for their previously unknown orbital inclinations to be determined: for Boötes b (Brogi et al. 2012), consistent with the measurement of Rodler et al. (2012), for 51 Peg b (Brogi et al. 2013), and for HD 179733 b (Brogi et al. 2014), which in turn lead to a determination of their masses.
1.2 Inversion layer in HD 209458b
HD 209458b is a well-studied exoplanet. It orbits a solar-type star once every , located from the Sun (see Table 1 for the physical properties of the system). Using the Spitzer Space Telescope Knutson et al. (2008) found evidence for emission in the IRAC bands at and , indicative of a thermal inversion layer in the atmosphere of the planet. A temperature inversion is caused by the absorption of incident starlight in a high-altitude layer of the atmosphere. On Earth, solar UV radiation is absorbed in the oxygen-ozone cycle, effectively absorbing all sunlight short of , heating the atmosphere and forming the stratosphere between and altitude (see e.g. Portmann & Solomon 2007, and references therein). In the atmospheres of hot-Jupiters, similar processes may occur, albeit in different chemical environments.
A number of compounds have been proposed to be capable of causing the inversion layer in the atmosphere of HD 209458b, one of which is TiO (Hubeny et al. 2003; Burrows et al. 2007), a diatomic molecule which is known to cause major absorption features throughout the optical and near-infrared spectra of M-dwarfs. An atmosphere containing a strong optical absorber like TiO would furthermore be consistent with HD 209458b’s low albedo (Rowe et al. 2006). Désert et al. (2008) reported a tentative detection of TiO absorption in the transmission spectrum of HD 209458b, though more recent broadband studies seem to support a lack of TiO absorption in the transmission spectra of other hot-Jupiters (see e.g. Huitson et al. 2013; Sing et al. 2013; Gibson et al. 2013; Bento et al. 2014). Also, models of hot-Jupiter atmospheres suggest that maintaining a significant gaseous TiO concentration at high altitudes is difficult due to condensation of TiO either at depth or at the night-side (Spiegel et al. 2009; Burrows & Orton 2009; Fortney et al. 2010; Parmentier et al. 2013). More recent observations of the dayside emission spectrum of HD 209458b, call into question the existence of the inversion layer altogether (Zellem et al. 2014; Diamond-Lowe et al. 2014). Interestingly, Stevenson et al. (2014) observed molecular absorption features in the near-infrared dayside spectrum of Wasp-12b, possibly attributable to TiO. This contrasts with the optical transmission spectra of the same planet, put forth by Sing et al. (2013). Thus, it is as yet unclear whether TiO is an important chemical component of hot-Jupiter atmospheres.
In this paper we apply a cross-correlation-based retrieval method to search for TiO in the transmission spectrum of HD 209458b at high spectral resolution. Section 2 describes our dataset and the application of the method. Section 3 discusses our results, followed by our conclusions in Section 4.
|Effective temperature (K)|
|Semi-major axis (AU)|
|Equilibrium temperature (K)|
Adopted from Mazeh et al. (2000).
Adopted from Knutson et al. (2007).
Adopted from Crossfield et al. (2012).
Adopted from the Exoplanet Orbit Database (Wright et al. 2011).
2 Observations and data reduction
2.1 Subaru data of HD 209458
The data used in this analysis were taken on October 25, 2002 using the High Dispersion Spectrograph (HDS) on the Subaru 8.2m telescope. The observations were performed by Narita et al. (2005) and were originally aimed at detecting absorption of sodium and a number of other atomic species in the transmission spectrum of HD 209458b. These were later used by Snellen et al. (2008) who did indeed detect sodium. The data consist of 31 exposures with a mean exposure time of which cover one 3h transit, plus 1.5h and 0.5h of baseline before and after, respectively. The spectrum is observed over 20 echelle orders, covering a wavelength range between and at a spectral resolution of ( resolution). Each spectral order is sampled at per pixel over 4100 pixels. The basic data reduction (such as bias subtraction, flat fielding and 1D spectral extraction) were performed by Snellen et al. (2008) and are described therein in detail.
2.2 HARPS data of M-dwarfs
Our analysis relies on cross-correlating the transmission spectra of HD 209458b with high-resolution template spectra of TiO-bearing models of HD 209458b’s absorption spectrum which rely on line lists published by Freedman et al. (2008, see sections 2.3 and 2.4). The optical spectra of M-dwarfs are dominated by TiO absorption features and are therefore natural benchmarks against which we tested the accuracy of these TiO templates. Many nearby M-dwarfs have been observed extensively at high spectral resolution in order to detect the radial velocity signature induced by potential orbiting companions. We downloaded high resolution spectra of 5 different M-dwarfs from the ESO data archive. These spectra were obtained by the HARPS spectrograph at ESO’s 3.6m La Silla telescope, between and at resolutions of . Information about the 5 stars is shown in Table 2. Barnard’s star is used throughout this work as the primary example, and is shown in the top panel of Figure 2 before the removal of all broadband features using a high-pass filter (the spectra of the other M-dwarfs were filtered in the same way), and after correcting for its radial velocity of (Nidever et al. 2002).
|Name||Spectral type||PID (ESO)||PI|
2.3 Post-processing and cross-correlation
The analysis of the high dispersion transmission spectra of HD 209458b differs at several points from that of infrared data used in previous high dispersion spectroscopy studies by our group (see Birkby et al. 2013; Brogi et al. 2012, 2013, 2014; de Kok et al. 2013), which are dominated by time-dependent telluric absorption features, mainly due to variations in airmass. In the optical, the atmosphere presents no significant absorption features except in some specific wavelength ranges. Rather, the spectrum is dominated by stellar absorption lines. Assuming that the stellar spectrum is stable during the transit, it can be removed by dividing each of the individual spectra in the time-series with the time-averaged stellar spectrum. In principle, this operation does not affect the planet signal because it significantly shifts in wavelength during the observations due to the change in the radial component of the orbital velocity during transit (from to ).
The full sequence of processing steps is summarized below. Intermediate data products are shown in Figure 3.
Broadband removal and normalization: To remove strong broadband components such as the stellar continuum and the spectrograph’s blaze function, the brightest exposure of each spectral order was divided into 50 wavelength regions. A 7-th order polynomial was fit through the maximum of each part to obtain the continuum level. This fit was removed from all other exposures through division, removing the strongest broadband features. Subsequently, the continuum of each of the exposures was obtained by dividing each exposure into the same 50 regions as above and again finding the maximum of each part. These continua were interpolated, and removed through division. This acts to remove residual broad-band variations not removed by the 7th-order polynomial fit, and normalizes the continuum of each exposure to unity. At this stage there were still significant residual broad-band variations, which were removed at a later stage (see step 7).
Alignment of spectra: The spectra slowly drift in wavelength due to the change in radial velocity of the observatory during the observations, and instrumental instability. These misalignments are of the order of a pixel and would adversely affect the removal of the time-average stellar spectrum if left untreated. In each spectral order a dozen strong stellar absorption lines were identified by eye. To these lines in each of the 31 exposures, Gaussian profiles were fitted to obtain the position of each line center. A linear fit through the misalignments of the different absorption lines was used to align the spectra to a common reference frame (a higher order fit was not warranted because the scatter in the determination of the line centroids was too high to identify higher order components).
Wavelength solution: As HD 209458 is a close solar analogue, we matched our time-averaged spectrum (see next step) of the HD 209458 system to a model solar spectrum (adopted from the Kurucz stellar atmosphere atlas) in a step-wise fashion. The strongest stellar absorption lines were identified and matched by eye. This produced a crude wavelength solution which was subsequently used to identify 30 - 40 less prominent lines in each order. We subsequently fitted Gaussian profiles to the cores of these lines in both the data and the solar template, producing a wavelength solution with a standard deviation of , corresponding to or of the spectral resolution. Because the solar spectrum is synthesized at rest, the wavelength solution is calibrated to the rest frame of HD 209458 and the mean radial velocity with respect to the observatory on Earth is thus automatically corrected for. Any time-dependent Doppler shifts had already been corrected in step 2. Therefore, each exposure is now in the rest frame of the host star.
Removal of stellar spectrum: Each pixel in each spectral order was time-averaged by taking the median of the aligned spectra, which corresponds to taking the median of each of the columns of the array in panel 2 of Figure 3. This average contains all time-constant stellar features, but only a small part of the time-varying planet signal, since it is shifted to a different column in each exposure. It is removed from each exposure through division. The width of planetary absorption lines is limited to , which is the resolution of the instrument. Between exposures, the planet signal is shifted by . Therefore, the planet signal is present in no more than 4 exposures per pixel column. These 4 exposures are part of a median with 27 other exposures which don’t contain the planet signal. Therefore, we expect that up to of the planet signal is removed along with the median stellar spectrum. This has no importance for the rest of the analysis, because we measure the significance of any TiO detection by our ability to retrieve injected model spectra, which suffer in the same way from this degradation (see section 3).
Removal of bad pixels and cosmics: The resulting residuals contain bad pixels and cosmic rays. All pixels deviating by more than from the mean were reset to the mean, removing all significant bad pixels and cosmic rays. 0.21% of all pixel values were affected.
Normalization by signal to noise: The noise in the residual spectrum varies with pixel position due to the presence of deep stellar absorption lines up until step 4. We divide each pixel in the residual spectrum by the variance of that pixel value (variance of each column of panel 3, Figure 2.3) to weigh down such pixels. Low SNR regions would otherwise add a significant amount of noise to the cross-correlation function.
Suppression of broadband residuals: The residuals feature broadband variations due to imperfect removal of instrumental effects and broad stellar lines. Latent broadband variations were characterized by smoothing each exposure with a Gaussian with a width of 40px () and were divided out of the data. The last panel of Figure 3 shows the residual of this final processing step, on which the cross-correlation is performed.
Cross correlation: The planet’s potential absorption features were extracted by cross-correlating the residual spectra with template spectra of TiO (see section 2.4) across the full wavelength range of the data i.e. between and . The templates were Doppler-shifted to velocities ranging between relative to the star, in steps of (which roughly corresponds to the sampling rate of the data, per pixel). At the Doppler shift equal to the radial velocity of the planet, the cross-correlation effectively co-adds all of the individual absorption lines of the molecule.
2.4 Model spectra
For the purpose of applying the cross-correlation method to the case of TiO in the atmosphere of HD 209458b (see below), we generated several model transmission spectra of atmospheres containing various concentrations of TiO, at a sampling resolution of . The line database stems from Freedman et al. (2008), who presented a modified line list based on the one published by Schwenke (1998). The model atmosphere has a temperature-pressure profile as shown in Figure 4, and includes scattering and molecular absorption due to TiO with volume mixing ratios (VMR) of , , , and . The solar photospheric abundance of titanium is (Asplund et al. 2006), so the transmission spectrum of TiO is modelled between and times the solar abundance of Ti. These mixing ratios were assumed to be uniform throughout the atmosphere. As such, no chemistry, phase-changes or dynamics were taken into account. The atmosphere was assumed to have a temperature inversion as displayed in the right panel of Figure 4 (adopted from Burrows et al. 2007), but because the background continuum originates from a star with an effective temperature of , the assumption of a particular T-P-profile does not significantly influence the template spectrum. To confirm the robustness of the template to a temperature mismatch, we also modelled a TiO-bearing atmosphere at a constant temperature of and cross-correlated it with the template shown in Figure 4. The result is shown in Figure 6, and discussed in more detail in section 3.1.
For benchmark purposes (see section 2.2), we also cross-correlated with synthetic M-dwarf spectra originating from the Kurucz and PHOENIX stellar atmosphere model suites which necessarily include models of the TiO molecule. Both models use a line list calculated by Plez (1998). The Kurucz model is a template spectrum of the M2V dwarf GL411 at a resolution of (Castelli & Kurucz 2003) and the PHOENIX model, a template of a generic M4V dwarf (, , , obtained from the PHOENIX online database as presented by Husser et al. (2013)) at a resolution of . They are shown together with the HARPS spectrum of Barnard’s star in Figure 2.
3 Results and discussion
The left panels of Figure 5 show the correlations between the HDS spectra of HD 209458b and our TiO template at five modeled VMRs (see section 2.4). The correlation coefficients are plotted as a function of the radial velocity to which the template is shifted (horizontal axis) for every individual exposure (i.e. time, vertical axis).
To assess the ability of the cross-correlation procedure to retrieve the template spectra from real data, we inject the templates into the data prior to cross-correlation using the known orbital velocity of HD 209458b, between steps 3 and 4 in section 2.3, similar to Snellen et al. (2010). Prior to injection, we convolve the template with a Gaussian to match the spectral resolution of the data , and convolve again with a box-function to account for the widening of the planetary absorption lines due to the changing Doppler-shift during the exposures. This gives rise to the middle panels in Figure 5. Injection of a TiO-bearing template at a VMR greater than clearly produces a slanted feature, the gradient of which is a measure of the orbital velocity of the planet at which the template was injected.
Given that the radial velocity of the planet is known at all times, each exposure (i.e. row in Figure 5) can be shifted to the rest-frame of the planet. This allows for co-addition of all exposures, the result of which is shown in the right panels in Figure 5. Here, the solid lines represents the correlation of the templates with the residuals of the data, and the dashed line represents the correlation after injection of the template. The signal-to-noise of the observed correlation peak is indicative of the significance at which the template would be retrieved, had it been identically present in the residuals. Evidently, the peak correlation strength decreases with decreasing VMR, hence the VMR at which the correlation peaks at the level, we consider to be the theoretical limiting VMR at which TiO would be retrievable by applying the current analysis on this data.
It is clear from the left panels of Figure 5 that the residuals do not correlate with our TiO templates. This non-detection is significant down to VMRs of , as evidenced by the successful retrieval of our templates after having artificially injected them into the data. This sensitivity limit corresponds to times the solar abundance of Ti (Asplund et al. 2006). The models of the atmosphere of HD209458b are greatly simplified approximates to the real physical environment. Especially the assumption of a uniform VMR across the entire atmosphere is questionable, as the atmosphere is not expected to be homogeneous. However, by extending the absorbing TiO layer uniformly throughout the atmosphere, the TiO absorption at that VMR is maximized. The limiting VMR of is therefore a true lower limit: If TiO would only be present at a VMR of in parts of the atmosphere, it would not be retrieved in this analysis.
Before invoking a plethora of reasons for this non-detection of TiO such as hazes, obscuring clouds, screening by other species, cold-traps or destructive chemistry, it should be noted that the sensitivity of the cross-correlation procedure critically depends on the accuracy of the line positions in the template spectrum. If the wavelength positions of the lines are inaccurate, the cross-correlation procedure may not add their contributions constructively. As shown below, this is a real concern.
3.1 Evaluation of template spectra
Because TiO absorption visibly dominates M-dwarf optical spectra, a high-resolution optical spectrum of Barnard’s star can be used to test the accuracy of our TiO template. Their mutual correlation should be very clear.
Figure 6 shows the cross-correlation between our TiO template at different temperatures, the stellar atmosphere models and the HARPS-spectrum of Barnard’s star. The TiO-bearing model of HD 209458b (VMR=, T-P profile as shown in Figure 4) shows a correlation peak of 0.35 with the stellar atmosphere models. The cross-correlation between our TiO templates at and peaks at 0.9. M-dwarf transmission spectra feature a host of other optical absorbers, which decrease the correlation between an M-dwarf spectrum and a pure-TiO template. We regard this as the main reason for the non-prefect correlation between our TiO template and the stellar atmosphere models, rather than a mismatch of the T-P profile associated with the template.
In contrast with the stellar atmosphere models, a marginal correlation peak of only is observed between the TiO template and the HARPS spectrum. From this we may conclude that our TiO-bearing model spectra do not accurately reproduce real-world TiO absorption at high resolution.
To test the extent at which stellar atmosphere models and real M-dwarf spectra match each other at high resolution, we cross-correlate Barnard’s star’s spectrum with the spectra of the other 5 M-dwarfs, and with the PHOENIX M4V model. This is shown in Figure 7. Regardless of differences in effective temperature, metallicity and signal-to-noise, real M-dwarf spectra correlate strongly with each other; with a peak of ¿0.75 in all cases. However, as is the case with the TiO template, the PHOENIX model correlates at less than 0.15 with all of the M-dwarf spectra. This strengthens our conclusion that the modelled spectrum of TiO is not representative of real TiO. The high correlation between the TiO template and the PHOENIX and Kurucz models, suggests that the TiO line lists have a common source, and the same uncertainties in the line positions seem to have persisted in the works of both Plez (1998) and Schwenke (1998). A visual comparison of the models and the spectrum of Barnard’s star, shows great similarity between all three models at high resolution, but a severe mismatch with the HARPS data (Figure 8).
Although the TiO models do not match the HARPS data well, our TiO template is able to partially retrieve the presence of TiO in Barnard’s star spectrum. The extent of the inaccuracies of the models as shown in Figure 8 near , raises the suspicion that the marginal correlation with Barnard’s star is not due to the haphazard alignment of model absorption lines. Instead, it is possible that certain absorption bands are + more accurately than others. To test the hypothesis that inaccuracies are associated with certain absorption bands only, we divide the available range of our HARPS data into smaller sections and again cross-correlate our model of HD 209458b with each section separately. The result is shown in Figure 9.
Over most of the wavelength range, our model does not correlate with the TiO spectrum. However above , the TiO template appears to better match the data, increasing in correlation towards longer wavelengths (a trend which may continue further towards the near infra-red). The four sections shown in in Figure 9 correlate maximally at slightly different radial velocity shifts (few km/s), meaning that ensembles of lines are collectively offset. This specifically indicates inaccuracies in the calculation of the energy levels of the TiO molecule.
3.2 Alternative template spectra
As shown above, the cross-correlation procedure is very sensitive to species with rich optical absorption spectra in the HDS residuals. As for TiO, the sensitivity reaches down to VMRs below the ppb level. Although we have also shown that the line-lists which lie at the basis of our TiO (and other’s) model spectra display extensive inaccuracies, the template is able to partially retrieve TiO features from Barnard’s star’s optical spectrum at wavelengths longwards of and is thus to a certain extent useful to probe the atmosphere of HD 209458b for the presence of TiO. We therefore cross-correlate the residuals with the TiO template at a VMR of , using only wavelengths between , the result of which is shown in the top row of figure 10.
Injecting the template into the residuals like before, reveals a theoretical sensitivity of (Top right panel of Figure 10). However, because the line list is only partially accurate at these wavelengths, this is only an upper limit: I.e. if line list was perfect, the presence of TiO in the atmosphere of HD 209458b at a volume-mixing ratio of would have been retrieved with a significance of . Though since the line-list suffers from uncharacterised inaccuracies, the true limiting VMR is not determined.
Earlier it has been shown (Figure 6) that if the temperature-pressure environment of the TiO template does not match that of the data, the correlation is not severely affected. Therefore, although the temperature of the photosphere of Barnard’s star is some higher than the atmosphere of HD 209458b at the terminator region, the TiO absorption present in the spectra of M-dwarfs may also be useful to identify TiO absorption in exoplanet atmosphere at high resolution. As shown before (Figure 7), the presence of noise in a cross-correlation template does not destroy the cross-correlation at high-resolution, and thus we propose to use the HARPS spectrum of Barnard’s star as a template for cross-correlation with the HDS data of HD 209458b.
Injecting the absorption spectrum of Barnard’s star into the residuals of the HDS data does not provide a physically meaningful sensitivity. Although our template at correlates significantly with the stellar atmosphere models, the absorption spectrum of an M-dwarf is not a good model for the absorption spectrum of the upper atmosphere of HD 209458b during transit, because the underlying physical parameters are very different: In addition to the different temperature-pressure environment, the presence of other species in an M-dwarf’s photosphere and the different continuum level (HD 209458 is a solar type star) make it difficult to adjust Barnard’s star’s absorption spectrum to form a template from which parameters such as the limiting VMR can be derived as before. Extracting a physically meaningful TiO absorption spectrum from real-world absorption spectra of M-dwarfs would require more work and is beyond the scope of this paper. For now, we only use Barnard’s star’s spectrum to qualitatively probe for the presence of TiO absorption in the HDS dataset. The resulting cross-correlation is shown in Figure 10. Again, no significant correlation is detected.
We have probed the optical transmission spectrum of HD 209458b for TiO absorption using the time-dependent radial velocity of the planet, which can be spectrally resolved at high resolution. At first glance our null-result seems to be in contrast with the possible detection of TiO by Désert et al. (2008), and in line with more recent observations of hot-Jupiters and models. Injection of model TiO spectra into the data shows that our cross-correlation retrieval method is sensitive to volume mixing ratios down to , and is thus theoretically capable of ruling out the presence of TiO at a stringent level.
However, this method relies on the accuracy of template spectra of TiO at resolutions in the order of . Although such templates are available in the literature, we observe that the line databases used to synthesise these spectra are based on model calculations of the TiO molecule which display widespread inaccuracies. We infer that the energy levels of the TiO molecule are not determined well enough to accurately synthesise absorption spectra at these resolutions. Nonetheless, we have found that the TiO line list used to generate our template spectrum is accurate to some extent at wavelengths between and , and that the spectra of M-dwarfs can also serve as cross-correlation templates. However, using a high-resolution spectrum of Barnard’s star and restricting our template to the aforementioned wavelength range, do not result in a significant correlation. It would therefore seem that the presence of a large concentration of TiO in the upper atmosphere of HD 209458b at the terminator region is unlikely. However, due to the uncertainties regarding the physical meaning of these template models, such a conclusion is tentative.
The lack of an accurate TiO line list is thus a critical hindrance in the application of this retrieval technique to search for the presence of TiO in the optical transmission spectrum of HD 209458b. Also, since the TiO molecule is among the best characterized molecules in astrophysics, we suspect that the absorption spectra of less common and more complicated species are likely to be even less accurately characterized. Furthermore, inaccuracies in the calculations of line positions do not only severely hinder the application of high-resolution retrieval techniques like the one employed in this work, but they also pervade widely used stellar atmosphere models. Modern-day stellar atmosphere codes employ highly advanced modelling techniques to poor line databases, which in the case of TiO is over a decade old.
Therefore we conclude that an increasing community of observational and theoretical astronomers would benefit from an immediate and extensive effort in improving the quality of absorption spectra of molecules like TiO, either through more sophisticated theoretical calculations or by means of direct measurements in the laboratory.
Acknowledgements.This work is part of the research programmes PEPSci and VICI 639.043.107, which are financed by the Netherlands Organisation for Scientific Research (NWO), and was also supported by the Leiden Observatory Huygens Fellowship. This work was performed in part under contract with the California Institute of Technology (Caltech)/Jet Propulsion Laboratory (JPL) funded by NASA through the Sagan Fellowship Program executed by the NASA Exoplanet Science Institute. Support for this work was provided in part by NASA through Hubble Fellowship grant HST-HF2-51336 awarded by the Space Telescope Science Institute. Also, this research has made use of the Exoplanet Orbit Database and the Exoplanet Data Explorer at exoplanets.org. Finally, we would like to thank Leigh Fetcher, Jaemin Lee and Patrick Irwin for providing the TiO line-list necessary to perform this analysis.
- Asplund, M., Grevesse, N., & Jacques Sauval, A. 2006, Nuclear Physics A, 777, 1
- Bento, J., Wheatley, P. J., Copperwheat, C. M., et al. 2014, MNRAS, 437, 1511
- Birkby, J. L., de Kok, R. J., Brogi, M., et al. 2013, MNRAS, 436, L35
- Brogi, M., de Kok, R. J., Birkby, J. L., Schwarz, H., & Snellen, I. A. G. 2014, A&A, 565, A124
- Brogi, M., Snellen, I. A. G., de Kok, R. J., et al. 2012, Nature, 486, 502
- Brogi, M., Snellen, I. A. G., de Kok, R. J., et al. 2013, ApJ, 767, 27
- Burrows, A., Hubeny, I., Budaj, J., Knutson, H. A., & Charbonneau, D. 2007, ApJ, 668, L171
- Burrows, A. & Orton, G. 2009, ArXiv e-prints
- Castelli, F. & Kurucz, R. L. 2003, in IAU Symposium, Vol. 210, Modelling of Stellar Atmospheres, ed. N. Piskunov, W. W. Weiss, & D. F. Gray, 20P
- Crossfield, I. J. M., Knutson, H., Fortney, J., et al. 2012, ApJ, 752, 81
- de Kok, R. J., Brogi, M., Snellen, I. A. G., et al. 2013, A&A, 554, A82
- Désert, J.-M., Vidal-Madjar, A., Lecavelier Des Etangs, A., et al. 2008, A&A, 492, 585
- Diamond-Lowe, H., Stevenson, K. B., Bean, J. L., Line, M. R., & Fortney, J. J. 2014, ArXiv e-prints
- Fortney, J. J., Shabram, M., Showman, A. P., et al. 2010, ApJ, 709, 1396
- Freedman, R. S., Marley, M. S., & Lodders, K. 2008, ApJS, 174, 504
- Gibson, N. P., Aigrain, S., Barstow, J. K., et al. 2013, MNRAS, 436, 2974
- Hubeny, I., Burrows, A., & Sudarsky, D. 2003, ApJ, 594, 1011
- Huitson, C. M., Sing, D. K., Pont, F., et al. 2013, MNRAS, 434, 3252
- Husser, T.-O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6
- Knutson, H. A., Charbonneau, D., Allen, L. E., Burrows, A., & Megeath, S. T. 2008, ApJ, 673, 526
- Knutson, H. A., Charbonneau, D., Noyes, R. W., Brown, T. M., & Gilliland, R. L. 2007, ApJ, 655, 564
- Lockwood, A. C., Johnson, J. A., Bender, C. F., et al. 2014, ApJ, 783, L29
- Mazeh, T., Naef, D., Torres, G., et al. 2000, ApJ, 532, L55
- Narita, N., Suto, Y., Winn, J. N., et al. 2005, PASJ, 57, 471
- Nidever, D. L., Marcy, G. W., Butler, R. P., Fischer, D. A., & Vogt, S. S. 2002, ApJS, 141, 503
- Parmentier, V., Showman, A. P., & Lian, Y. 2013, A&A, 558, A91
- Plez, B. 1998, A&A, 337, 495
- Portmann, R. W. & Solomon, S. 2007, Geochim. Res. Lett., 34, 2813
- Rodler, F., Kürster, M., & Barnes, J. R. 2013, MNRAS, 432, 1980
- Rodler, F., Lopez-Morales, M., & Ribas, I. 2012, ApJ, 753, L25
- Rowe, J. F., Matthews, J. M., Seager, S., et al. 2006, ApJ, 646, 1241
- Schwenke, D. W. 1998, Faraday Discussions, 109, 321
- Sing, D. K., Lecavelier des Etangs, A., Fortney, J. J., et al. 2013, MNRAS, 436, 2956
- Snellen, I. A. G., Albrecht, S., de Mooij, E. J. W., & Le Poole, R. S. 2008, A&A, 487, 357
- Snellen, I. A. G., de Kok, R. J., de Mooij, E. J. W., & Albrecht, S. 2010, Nature, 465, 1049
- Spiegel, D. S., Silverio, K., & Burrows, A. 2009, in AAS/Division for Planetary Sciences Meeting Abstracts, Vol. 41, AAS/Division for Planetary Sciences Meeting Abstracts 41, 68.03
- Stevenson, K. B., Bean, J. L., Seifahrt, A., et al. 2014, AJ, 147, 161
- Torres, G., Winn, J. N., & Holman, M. J. 2008, ApJ, 677, 1324
- Wright, J. T., Fakhouri, O., Marcy, G. W., et al. 2011, PASP, 123, 412
- Zellem, R. T., Lewis, N. K., Knutson, H. A., et al. 2014, ApJ, 790, 53