Issue |
A&A
Volume 589, May 2016
|
|
---|---|---|
Article Number | A62 | |
Number of page(s) | 5 | |
Section | Planets and planetary systems | |
DOI | https://doi.org/10.1051/0004-6361/201527881 | |
Published online | 13 April 2016 |
The GTC exoplanet transit spectroscopy survey
III. No asymmetries in the transit of CoRoT-29b⋆
1 Instituto de Astrofísica de Canarias (IAC), Vía Láctea s/n, 38205 La Laguna, Spain
e-mail: epalle@iac.es
2 Departamento de Astrofísica, Universidad de La Laguna, 38205 La Laguna, Tenerife, Spain
3 Key Laboratory of Planetary Sciences, Purple Mountain Observatory, Chinese Academy of Sciences, 210008 Nanjing, PR China
4 Institute of Planetary Research, German Aerospace Center, Rutherfordstrasse 2, 12489 Berlin, Germany
5 Univ. Grenoble Alpes, IPAG, 38000 Grenoble, France
6 CNRS, IPAG, 38000 Grenoble, France
7 Sub-department of Astrophysics, Department of Physics, University of Oxford, Oxford, OX1 3RH, UK
8 Institut für Astrophysik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
9 Theoretical Meteorology group, Klimacampus, University of Hamburg, Grindelberg 5, 20144 Hamburg, Germany
Received: 2 December 2015
Accepted: 30 January 2016
Context. The launch of the exoplanet space missions obtaining exquisite photometry from space has resulted in the discovery of thousands of planetary systems with very different physical properties and architectures. Among them, the exoplanet CoRoT-29b was identified in the light curves the mission obtained in summer 2011, and presented an asymmetric transit light curve, which was tentatively explained via the effects of gravity darkening.
Aims. Transits of CoRoT-29b are measured with precision photometry, to characterize the reported asymmetry in their transit shape.
Methods. Using the OSIRIS spectrograph at the 10-m GTC telescope, we perform spectro-photometric differential observations, which allow us to both calculate a high-accuracy photometric light curve, and a study of the color-dependence of the transit.
Results. After careful data analysis, we find that the previously reported asymmetry is not present in either of two transits, observed in July 2014 and July 2015 with high photometric precisions of 300 ppm over 5 min. Due to the relative faintness of the star, we do not reach the precision necessary to perform transmission spectroscopy of its atmosphere, but we see no signs of color-dependency of the transit depth or duration.
Conclusions. We conclude that the previously reported asymmetry may have been a time-dependent phenomenon, which did not occur in more recent epochs. Alternatively, instrumental effects in the discovery data may need to be reconsidered.
Key words: planets and satellites: atmospheres / planet-star interactions / extraterrestrial intelligence / astrobiology / eclipses / techniques: spectroscopic
Light curves are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/589/A62
© ESO, 2016
1. Introduction
The launches of the exoplanet space missions CoRoT (Baglin 2003) and Kepler (Borucki et al. 2003), returning exquisite photometry from space of thousands of stars, represented a tremendous milestone in exoplanet research. These missions have discovered hundreds of confirmed exoplanets, from hot Jupiters to rocky planets (Léger et al. 2009; Borucki et al. 2011), some of them possibly habitable (Borucki et al. 2013; Jenkins et al. 2015).
Among the possibilities that such exquisite photometry allow is the exploration of the interaction between stars and planets. Gravity darkening is produced by the oblateness of a star due to fast rotation or the gravitational pull of a close companion. Indeed, gravity darkening, which takes into account the relationship between surface brightness and local effective gravity, is necessary to explain the light curves of many binary systems (Espinosa Lara & Rieutord 2012).
In the context of exoplanets, two transiting hot Jupiter systems monitored with the Kepler spacecraft, Kepler-13Ab (Mazeh et al. 2012) and HAT-P-7b (Pál et al. 2008) have shown evidences of gravity darkening in their asymmetric transit light curves (Barnes et al. 2011; Morris et al. 2013; Masuda 2015).
More recently, Cabrera et al. (2015) reported the discovery of CoRoT-29b, a Jupiter-sized planet (M = 0.85 ± 0.20 MJup; R = 0.90 ± 0.16 RJup) orbiting an oblate star with an orbital period of ≈2.85 days. Because of the asymmetric shape of its transit light curve, observed both in the CoRoT data and subsequent ground-based followup, transit fitting models had to include the possibility that the stellar surface was not strictly spherical. Thus gravity darkening effects were called upon for a possible explanation, which the authors themselves did not consider definitive.
Still, CoRoT-29 is a relatively faint star (mV = 15.56) and some doubts remained about the significance of the asymmetry detection. In such cases, the use of large collecting area ground-based telescopes may reach the sensitivity that space facilities cannot. Here, we have used the 10-m GTC telescope to investigate the transit light curve of CoRoT-29b in search for gravity darkening signatures.
![]() |
Fig. 1 Slit image (left) and its corresponding dispersed 2D spectra (right). The left of the two bright stars is CoRoT-29. The yellow lines mark the best aperture size that minimizes the light-curve standard deviation. The aperture of the reference star includes the flux of the faint close-companion stars. |
Observing logs of the two transit datasets, including date, times, observing setup and atmospheric conditions.
2. Observations and data reduction
Two full transits of CoRoT-29b were observed with OSIRIS (Optical System for Imaging and low-Intermediate-Resolution Integrated Spectroscopy; Sánchez et al. 2012) spectrograph at the 10-m Gran Telescopio Canarias (GTC). The observed transits epochs were 31 July 2014 and 8 July 2015.
Observations were taken in the 2 × 2 binning mode, a readout speed of 200 kHz, and a gain of 0.95 e-/ADU, a combination that resulted in a readout noise of 4.5 e-. The R1000R grating was used and a custom-made 12′′-wide slit was used. This long-slit spectroscopy configuration covers the spectral range of 520−1040 nm, although wavelengths redder than 910 nm are affected by fringing (Murgas et al. 2014) and were not used in this work.
As illustrated in Fig. 1, both CoRoT-29 and the comparison star (2MASS-18353616+062854), separated in the sky 9.7 arcsec, were located on OSIRIS CCD1, close to the telescope optical center, which reduces considerably all telescope systematics (Chen et al. 2016; Nortmann et al. 2016). The position angle of the reference star with respect to the target (measured East of North) was −33.98 degrees. A detailed log of the observations can be found in Table 1.
Data were reduced using standard procedures, and include bias substraction, flat fielding, and corrections for spatial distortion and wavelength shifts to bring all the spectra to a reference. This methodology is described in detail in Murgas et al. (2014) and (Chen et al. 2016). Spectra extraction was performed using IRAF1. While CoRoT-29b has no close companion, our comparison star has two nearby faint companions which were completely included in the extraction aperture. The optimal aperture diameters were 22 and 21 pixels for the first and second night, respectively, which were retrieved using the built-in optimal extraction procedure in IRAF. Figure 2 shows the extracted spectrum of both CoRoT-29 and the slightly redder reference star.
![]() |
Fig. 2 Extracted spectra of CoRoT-29 (in red), and the reference star (in blue). The spectra were averaged over all out-of-transit observations. The color shaded areas indicate the pass-bands that are used to generate the spectro-photometric light curves. Note that the spectral type of the comparison is a bit redder than CoRoT-29, but has similar spectral features. Some spectral lines, common in stellar and planetary atmospheres, are also marked. The spectra are not corrected for instrumental response. |
For each image we calculate the drifts in spatial and also in dispersion direction. We monitored the drift in spatial direction by fitting a Gaussian function to the stellar profile tracing the position of the peaks. The full width at half maximum (FWHM) of the fitted Gaussian profile in the spatial direction is used to monitor the seeing variations. The drift in the dispersion direction was determined by fitting a Gaussian function to the absorption lines (e.g. Hα).
CoRoT-29b refined planet parameters.
3. Transit white and color light curves
Transit light curves were obtained by integrating the fluxes of each spectrum of the time series and then taking the ratio between the integrated flux of the target and the reference star. In the case of the white light curve, the flux was integrated between 515 and 915 nm.
To study the possible color signature of CoRoT-29b, we defined eight broad filters, whose spectral coverage is illustrated in Fig. 2. Redwards of 915 nm we do not use the data as they might be affected by fringing. The region around the telluric oxygen band at 760 nm is also discarded due to its poor signal-to-noise ratio (Parviainen et al. 2016).
3.1. Light curve fitting
To derive the transit parameters, we fit the flux ratio between the target and comparison star light curves with an analytic transit model, simultaneously multiplied by a parametric decorrelation baseline model, of the form:
where the analytic transit model
is taken from Mandel & Agol (2002), and is described by the planet’s position relative to the star, assuming an opaque, dark sphere eclipsing a spherical star. The best parametric decorrelation baseline model, B(cj) = c0 + c1sy + c2t + c3t2, was selected among several others using the BIC evaluation criteria (Schwarz 1978), and consists of a polynomial combination of various state vectors, including the FWHM in the spatial direction sy and time t. To determine the best-fitting model parameters and their associated uncertainties, we use a modified version of the Transit Analysis Package (Gazak et al. 2012) to perform a Markov chain Monte Carlo (MCMC) analysis.
The free parameters in our light-curve model were the mid-transit time Tmid, the orbital inclination i, the scaled semi-major axis a/R⋆, the planet-to-star radius ratio Rp/R⋆, two limb-darkening coefficients (u1, u2), and four baseline coefficients (c0, c1, c2, c3). The orbital period P was first fixed to the literature value and then refined with the newly fitted mid-transit times. The orbital eccentricity e and argument of periastron w were fixed to zero. The transit model uses a quadratic limb-darkening law, with the derived coefficients treated as Gaussian priors in the fitting of light curves. The full methodology employed in the model fitting is extensively detailed in (Chen et al. 2016).
In Table 2, the refined planet parameters for CoRoT-29b are given. The orbital period P and the zero epoch T0 are calculated by fitting a linear function to our two transit epochs together with literature epochs. The rest of parameters are calculated by individually fitting the light curve of each epoch, but also by jointly fitting both transits together. These values are compared to those in the discovery paper by Cabrera et al. (2015).
![]() |
Fig. 3 White light curve of CoRoT-29b. It is a summation of all the color-shaded passbands as shown in Fig. 2. From top to bottom: (1) normalized flux ratio between CoRoT-29 and its reference star, together with the best-fitting model; (2) white light curve after removing the modeled systematics, for display purpose; (3, 4) best-fitting light-curve residuals. |
![]() |
Fig. 4 Color light curves of CoRoT-29b. Left and right panels: data from July 31st, 2014 and July 8th 2015, respectively. The light curves are shown after removing the modeled systematics, and overplotted are the best-fitting light-curve models. |
![]() |
Fig. 5 Transmission spectrum of CoRoT-29b. The dashed line indicates the best-fitting Rp/R⋆ derived from the white light curve. Blue and red colors indicate the two observing epochs. The gray-shaded area marks the ± 5H scale height. Within 1σ error bars the spectrum is flat. |
4. Discussion and conclusions
Figure 3 presents the white light curve for the two datasets. In the top panel the raw flux ratios between CoRoT-29 and the comparison star are shown, where effects of airmass and seeing are still present. In the middle panel of the same figure, the light curves have been normalized and all systematic effects, except the analytic transit model, have been removed. Overplotted to the data are our best fits to each of the transits. It is remarkable that even at the faint magnitude of the CoRoT-29 (mv = 15.56), using GTC we are able to reach precisions down to ≈300 ppm in 5 min exposure time; this is a factor ≈2.2 above photon noise. The residuals of the transit fit do not show any clearly identifiable features attributable to red noise, or an unaccounted for planetary signal.
To be sure that the decorrelation process is not removing real asymmetries in the transit curve, we have generated mock light curves introducing an asymmetry in the form of a stellar spot (not shown). In that case using the same fitting methodology and a symmetric transit, the asymmetry is recovered, and only when using an asymetric transit profile, taking into accout the stellar spot effects, are the residuals of the fit similar to those in Fig. 3.
To further explore possible transit anomalies, we plot in Fig. 4 the color light curves for each of our spectral bins for the two observed epochs. It can be seen that the symmetry of the transit is maintained over all spectral bins and in both epochs. The only exception would be on the bluest light curves from July 8th 2015, where deviations from the transit model are larger. However, these anomalies are not only of the opposite sign to those reported in the original discovery paper by Cabrera et al. (2015), but are more likely due to low signal to noise ratios in this wavelength bin.
While for such a faint star the signal-to-noise ratio is far too low to perform transmission spectroscopy to explore the atmosphere of CoRoT-29b, it is still useful to look at the transit parameters as a function of wavelength. In Fig. 5 the retrieved transit depth as a function of wavelength is plotted, for each of the two epochs. All data points are consistent within 1−σ error bars, indicating that there is no change in transit depth with color. Overall, Figs. 4 and 5 lead us to the conclusion that there are no asymmetries in the observed transits of CoRoT-29b, at any wavelength interval, that can be associated to gravity darkening or other effects.
In conclusion, we observe two transits of CoRoT-29b, with a precision of around 300 ppm in the white light curve and about 400−700 ppm in several color passbands covering the optical range from 515−915 nm. We find no signs of an anomalous or asymmetric light curve as reported in previous studies. We note that the transits reported by Cabrera et al. (2015) were taken by the CoRoT mission in summer 2011. Cabrera et al. (2015) report also two ground-based light-curves, one taken with the 2-m Faulkes Telescope North (FTN) on 28 May 2012, and one taken on 14 July 2014 with the IAC80 telescope. Both CoRoT and the FTN data showed transits with clear and consistent asymmetries. In the IAC80 data, – taken only 2 weeks before our first GTC run – the relatively large noise permits however only a marginal identification of the asymmetry, if it was present at all. While it is not impossible that both the CoRoT and the FTN data suffered from unfortunate instrumental effects causing a similar asymmetry, another explanation is that the poorly understood source of this asymmetry disappeared sometimes between 2012 and 2014.
The Kepler mission has recently shown examples of uncommon light curves, with the detection of a possible planetary debris disk around a white dwarf (Vanderburg et al. 2015), or the puzzling behavior of the light curve of KIC 8462852 (Boyajian et al. 2016). In the near future, with the launch of missionw such as TESS and PLATO, which will cover very large fractions of the sky, the number of exoplanet discoveries around very bright stars will increase. As the total number of planets increases, the number of planets which may present challenging or unexpected light curves will also increase. Some might be first confirmations of theoretically predicted phenomena, but some may turn out to be systematic noises due to the low signal to noise of the observations or remain difficult to explain. In all cases, we demonstrated here that, using spectro-photometric techniques and having adequate instrumentation, large aperture ground-based telescopes can greatly contribute to the study of these phenomena.
Acknowledgments
Based on observations made with the Gran Telescopio Canarias (GTC), installed in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias, in the island of La Palma. This work is partly financed by the Spanish Ministry of Economics and Competitiveness through projects ESP2013-48391-C4-2-R, ESP2014-57495-C2-1-R, and AYA2012-39346-C02-02. G.C. also acknowledges the support by the Natural Science Foundation of Jiangsu Province (Grant No. BK20151051) and the National Natural Science Foundation of China (Grant No. 11503088). R.A. acknowledges the Spanish Ministry of Economy and Competitiveness (MINECO) for the financial support under the Ramn y Cajal program RYC-2010-06519. F.M. acknowledges the support of the French Agence Nationale de la Recherche (ANR), under the program ANR-12-BS05-0012 Exo-atmos. S.H. acknowledges financial support from the Spanish Ministry of Economy and Competitiveness (MINECO) under the 2011 Severo Ochoa Program MINECO SEV-2011-0187.
References
- Baglin, A. 2003, Adv. Space Res., 31, 345 [NASA ADS] [CrossRef] [Google Scholar]
- Barnes, J. W., Linscott, E., & Shporer, A. 2011, ApJS, 197, 10 [NASA ADS] [CrossRef] [Google Scholar]
- Borucki, W. J., Koch, D. G., Lissauer, J. J., et al. 2003, in Future EUV/UV and Visible Space Astrophysics Missions and Instrumentation, eds. J. C. Blades, & O. H. W. Siegmund, SPIE Conf. Ser., 4854, 129 [Google Scholar]
- Borucki, W. J., Koch, D. G., Basri, G., et al. 2011, ApJ, 728, 117 [NASA ADS] [CrossRef] [Google Scholar]
- Borucki, W. J., Agol, E., Fressin, F., et al. 2013, Science, 340, 587 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Boyajian, T. S., LaCourse, D. M., Rappaport, S. A., et al. 2016, MNRAS, 457, 3988 [NASA ADS] [CrossRef] [Google Scholar]
- Cabrera, J., Csizmadia, S., Montagnier, G., et al. 2015, A&A, 579, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Chen, G., Palle, E., Murgas, F., et al. 2016, A&A, submitted [Google Scholar]
- Espinosa Lara, F., & Rieutord, M. 2012, A&A, 547, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gazak, J. Z., Johnson, J. A., Tonry, J., et al. 2012, Adv. Astron., 2012, 30 [Google Scholar]
- Jenkins, J. M., Twicken, J. D., Batalha, N. M., et al. 2015, AJ, 150, 56 [NASA ADS] [CrossRef] [Google Scholar]
- Léger, A., Rouan, D., Schneider, J., et al. 2009, A&A, 506, 287 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
- Mandel, K., & Agol, E. 2002, ApJ, 580, L171 [NASA ADS] [CrossRef] [Google Scholar]
- Masuda, K. 2015, ApJ, 805, 28 [NASA ADS] [CrossRef] [Google Scholar]
- Mazeh, T., Nachmani, G., Sokol, G., Faigler, S., & Zucker, S. 2012, A&A, 541, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Morris, B. M., Mandell, A. M., & Deming, D. 2013, ApJ, 764, L22 [NASA ADS] [CrossRef] [Google Scholar]
- Murgas, F., Pallé, E., Zapatero Osorio, M. R., et al. 2014, A&A, 563, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nortmann, L., Pallé, E., Murgas, F., et al. 2016, A&A, submitted [Google Scholar]
- Pál, A., Bakos, G. Á., Torres, G., et al. 2008, ApJ, 680, 1450 [NASA ADS] [CrossRef] [Google Scholar]
- Parviainen, H., Pallé, E., Nortmann, L., et al. 2016, A&A, 585, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sánchez, B., Aguiar-González, M., Barreto, R., et al. 2012, in SPIE Conf. Ser., 8446, 4 [Google Scholar]
- Schwarz, U. J. 1978, A&A, 65, 345 [NASA ADS] [Google Scholar]
- Vanderburg, A., Johnson, J. A., Rappaport, S., et al. 2015, Nature, 526, 546 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Observing logs of the two transit datasets, including date, times, observing setup and atmospheric conditions.
All Figures
![]() |
Fig. 1 Slit image (left) and its corresponding dispersed 2D spectra (right). The left of the two bright stars is CoRoT-29. The yellow lines mark the best aperture size that minimizes the light-curve standard deviation. The aperture of the reference star includes the flux of the faint close-companion stars. |
In the text |
![]() |
Fig. 2 Extracted spectra of CoRoT-29 (in red), and the reference star (in blue). The spectra were averaged over all out-of-transit observations. The color shaded areas indicate the pass-bands that are used to generate the spectro-photometric light curves. Note that the spectral type of the comparison is a bit redder than CoRoT-29, but has similar spectral features. Some spectral lines, common in stellar and planetary atmospheres, are also marked. The spectra are not corrected for instrumental response. |
In the text |
![]() |
Fig. 3 White light curve of CoRoT-29b. It is a summation of all the color-shaded passbands as shown in Fig. 2. From top to bottom: (1) normalized flux ratio between CoRoT-29 and its reference star, together with the best-fitting model; (2) white light curve after removing the modeled systematics, for display purpose; (3, 4) best-fitting light-curve residuals. |
In the text |
![]() |
Fig. 4 Color light curves of CoRoT-29b. Left and right panels: data from July 31st, 2014 and July 8th 2015, respectively. The light curves are shown after removing the modeled systematics, and overplotted are the best-fitting light-curve models. |
In the text |
![]() |
Fig. 5 Transmission spectrum of CoRoT-29b. The dashed line indicates the best-fitting Rp/R⋆ derived from the white light curve. Blue and red colors indicate the two observing epochs. The gray-shaded area marks the ± 5H scale height. Within 1σ error bars the spectrum is flat. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.