No variations in transit times for Qatar-1 b
The transiting hot Jupiter planet Qatar-1 b was presented to exhibit variations in transit times that could be of perturbative nature. A hot Jupiter with a planetary companion on a nearby orbit would constitute an unprecedented planetary configuration, important for theories of formation and evolution of planetary systems. We performed a photometric follow-up campaign to confirm or refute transit timing variations. We extend the baseline of transit observations by acquiring 18 new transit light curves acquired with 0.6-2.0 m telescopes. These photometric time series, together with data available in the literature, were analyzed in a homogenous way to derive reliable transit parameters and their uncertainties. We show that the dataset of transit times is consistent with a linear ephemeris leaving no hint for any periodic variations with a range of 1 min. We find no compelling evidence for the existence of a close-in planetary companion to Qatar-1 b. This finding is in line with a paradigm that hot Jupiters are not components of compact multi-planetary systems. Based on dynamical simulations, we place tighter constraints on a mass of any fictitious nearby planet in the system. Furthermore, new transit light curves allowed us to redetermine system parameters with the precision better than that reported in previous studies. Our values generally agree with previous determinations.
The transiting planet Qatar-1 b (Alsubai et al. 2011) is the first planetary body discovered by the Qatar Exoplanet Survey (Alsubai et al. 2013). It is a typical hot Jupiter on a 1.42 d circular orbit. Alsubai et al. (2011) give a mass of (Jupiter’s mass), but new radial velocity (RV) measurements by Covino et al. (2013) show that the planet has the greater mass of . The published studies provide consistent determinations of a planetary radius which is greater by 16–18% than (Jupiter’s radius). The planet produces transits which are nearly 1.5 h long and 2.3% deep in flux. The host star, GSC 4240-470 ( mag), is a metal-rich dwarf of spectral type K3. It exhibits a moderate chromospheric activity, which could be partly induced by the close-in planet (Covino et al. 2013). The RV observations throughout a transit show a signature of the Rositter-McLaughlin effect that is consistent with a sky-projected obliquity close to zero (Covino et al. 2013).
|#||Telescope||Diameter (m)||CCD detector|
|Observatory||Size of matrix and field of view|
|1||Ritchey-Chrétien-Coudé Telescope||2.0||Roper Scientific VersArray 1300B||1|
|National Astronomical Observatory Rozhen, Bulgaria||pixels,|
|2||Ritchey-Chrétien Telescope||1.5||Roper Scientific VersArray 2048B||4|
|Sierra Nevada Observatory (OSN), Spain||pixels,|
|3||Trebur 1-meter Telescope||1.2||SBIG STL-6303||3|
|Michael Adrian Observatory, Trebur, Germany||pixels,|
|4||Schmidt Teleskop Kamera (Mugrauer & Berthold 2010)||0.9||E2V CCD-42-40||3|
|University Observatory Jena, Germany||pixels,|
|5||Cassegrain Telescope||0.6||SBIG STL-1001||6|
|Nicolaus Copernicus University, Poland||pixels,|
|6||Cassegrain Telescope||0.6||SBIG STL-6303||1|
|Volkssternwarte Kirchheim, Germany||pixels,|
|Date UT (epoch)||Telescope||Filter||Sky conditions|
|2011 Aug 01 (181)||Volk.Kirch. 0.6 m||none||0.87||3.68||Clear|
|2012 Dec 13 (533)||Torun 0.6 m||none||1.06||2.43||Clear|
|2013 Sep 03 (719)||Torun 0.6 m||none||1.51||2.51||Clear|
|2013 Sep 30 (738)||Torun 0.6 m||none||1.24||2.09||Photometric|
|2013 Nov 26 (778)||Trebur 1.2 m||none||1.39||1.18||Photometric|
|2013 Dec 20 (795)||Jena 0.9 m||1.12||2.67||Clear|
|2014 Mar 25 (862)||Trebur 1.2 m||none||0.40||4.36||Passing clouds|
|2014 Apr 02 (867)||Rozhen 2.0 m||1.22||1.47||Photometric|
|Torun 0.6 m||none||1.82||2.36||Clear|
|2014 May 16 (898)||OSN 1.5 m||1.09||2.21||Mostly clear, some high clouds|
|2014 May 26 (905)||OSN 1.5 m||1.33||1.39||Mostly photometric|
|2014 Jun 22 (924)||OSN 1.5 m||1.09||1.61||Mostly photometric, occasionally|
|passing high clouds before ingress|
|2014 Jul 01 (931)||Trebur 1.2 m||0.88||2.64||Photometric, first half of series|
|partly obscured by the dome slit|
|2014 Sep 10 (981)||OSN 1.5 m||1.33||1.29||Photometric|
|2014 Sep 18 (986)||Jena 0.9 m||1.16||2.25||Photometric|
|2014 Sep 27 (993)||Jena 0.9 m||1.14||2.02||Photometric|
|Torun 0.6 m||0.95||3.74||Clear, high humidity|
|2014 Oct 27 (1014)||Torun 0.6 m||0.86||3.45||High humidity, occasionally passing|
More recently, von Essen et al. (2013) have found indications for long-term variations in transit timing residuals from a linear ephemeris. The periodicity of this signal is postulated to be close to 190 or 390 d with a peak-to-peak amplitude of about 2 min. These variations were interpreted as a result of the exchange of energy and angular momentum via gravitational interaction with an unseen, close-in planetary companion (Holman & Murray 2005; Agol et al. 2005). Such planetary configuration would be unprecedented because hot Jupiters appear to be alone or accompanied by other planetary bodies on wide orbits (e.g., Steffen et al. 2012). On the other hand, studies of candidates for hot planets from a Kepler sample show that some fraction of these objects reveal periodic variations in transit timing that could be dynamically induced (Szabó et al. 2013). There may be only a small fraction of hot Jupiters which are a part of compact multiple-planet systems. If such systems are confirmed, their properties will place constraints on planet formation and migration theories. In some circumstances, migration models, based on tidal interactions between a giant planet and a circumstellar gas disk, predict the presence of neighboring low-mass planets, close to or trapped in interior or exterior mean-motion resonances (MMRs) with hot Jupiters (e.g., Podlewska & Szuszkiewicz 2009). Such resonant or near-resonant companions would induce orbital perturbations producing detectable variations in transit times. The analysis of multiple-planet candidate systems shows that most planetary systems are not resonant, but the distribution of planet period ratios shows clumping just wide of low-order resonances, such as 2:1 or 3:2 (Fabrycky et al. 2014). The origin of this pattern is unclear. The question whether hot Jupiters are accompanied by planets on nearby orbits or not, and what is the fraction of such planetary configurations still remains open.
In this context, Qatar-1 b with its preliminary detection of the transit time variation (TTV) signal could become a prototype of a new group of multiple-planet packed systems with hot Jupiters. As the system deserves further investigation, we organized an observing campaign to confirm or refute perturbations in the orbital motion of Qatar-1 b by the method of transit timing.
We acquired 18 new light curves for 16 transit events with six 0.6-2.0 m telescopes located across Europe (in Bulgaria, Poland, Germany, and Spain). The characteristics of individual instruments and detectors are given in Table 1. Two transits were observed simultaneously from two distinct sites. Most of the observations were performed in filter, in which the instrumental setups were found to be the most effective. We also acquired light curves in white light, i.e. without any filter, to increase the signal-to-noise ratio and obtain more precise data for timing purposes. Telescopes were moderately defocused to observe the Qatar-1 star with exposure times between 10 and 60 s, depending on mirror diameter, filter, and observing conditions. For a given light curve, the exposure time was kept fixed throughout a run. The timestamps were saved in coordinated universal time (UTC), provided by a GPS system or the network time protocol software. Details on individual runs are presented in Table 2.
Photometric time series were obtained with the differential aperture photometry applied to science images after bias and dark subtraction followed by flat-field correction with sky flats. The radius of the photometric aperture was optimized for each dataset to achieve the smallest scatter in the out-of-transit light curves. Magnitudes were determined against available comparison stars. Each light curve was de-trended by fitting a parabola along with a trial transit model, then converted to fluxes and normalized to unity outside of the transit event. The timestamps were converted to barycentric Julian dates in barycentric dynamical time (, Eastman et al. 2010).
3.1 System parameters
The Transit Analysis Package (TAP, Gazak et al. 2012) was used to model the set of 18 new transit light curves, supplemented with 31 complete or partial light curves available in the literature. Covino et al. (2013) provide 5 high-quality photometric time series acquired between May 2011 and September 2012 with the 1.82 m Asiago (Italy) and 1.23 m Calar Alto (Spain) telescopes. von Essen et al. (2013) make available 26 light curves. They were obtained from March 2011 to August 2012 with the 1.2 m Oskar-Lühning telescope at Hamburg Observatory (Germany) and the 0.6 m Planet Transit Study Telescope at the Mallorca Observatory (Spain). The TAP code uses the Markov Chain Monte Carlo (MCMC) method with the Metropolis-Hastings algorithm and a Gibbs sampler to find the best-fit transit model, based on the analytical approach of Mandel & Agol (2002). As the photometric time series may be afflicted by time-correlated noise, the wavelet-based algorithm of Carter & Winn (2009) is employed to robustly estimate parameter uncertainties. A quadratic limb darkening (LD) law (Kopal 1950) is implemented to approximate the flux distribution across the stellar disk. The values of linear and quadratic LD coefficients, and , respectively, were linearly interpolated from tables of Claret & Bloemen (2011) using a dedicated tool of the EXOFAST package (Eastman et al. 2013). The stellar parameters, such as effective temperature K, surface gravity (in cgs units), and metallicity , were taken from Covino et al. (2013). The efficient maximum of the energy distribution in white light was found to coincide with the band, so -band LD coefficients were adopted for white light curves.
|Parameter||This paper||Alsubai et al. (2011)
||Covino et al. (2013)||von Essen et al. (2013)|
|Orbital inclination, (deg)|
|Scaled semi-major axis,||—|
|Planet to star radii ratio,|
|Transit total duration, (min)||—||—|
|Semi-major axis, (AU)||—|
|Planetary radius, ()||—|
|Planetary density, ()||—|
|Stellar mass, ()||—|
|Stellar radius, ()||—|
|Stellar density, ()||—|
|Stellar luminosity, ()||—||—||—|
In a final iteration of the fitting procedure, the orbital inclination , the semimajor-axis scaled by stellar radius , and the planetary to stellar radii ratio were linked for all light curves and fitted simultaneously. The mid-transit times were determined independently for each epoch. The orbital period was kept fixed at a value refined in Sect. 3.2. The LD coefficients were allowed to vary around the theoretical values under the Gaussian penalty of to account for uncertainties in stellar parameters and possible systematic errors of the theoretical predictions (e.g., Müller et al. 2013). The orbit of the planet was assumed to be circular (Covino et al. 2013). Although the individual light curves are de-trended, the TAP code was allowed to take possible linear trends into account and include their uncertainties in the total error budget of the fit. We used ten MCMC chains, each containing steps. The first 10% of the results were discarded from each chain to minimize the influence of the initial values of the parameters. The best-fitting parameters, determined as the median values of marginalized posteriori probability distributions, and 1 uncertainty estimates, defined by a range of 68.3% values of the distributions, are given in Table 3. The new light curves are plotted in Fig. 1, together with the transit model and residuals.
The parameters derived from the transit model, combined with literature quantities determined from spectroscopic data, allowed us to redetermine the stellar mass , luminosity , and system age. The mean stellar density was calculated from transit parameters, independent of theoretical stellar models, using the formula
where is the universal gravitational constant. Then, the host star was placed at a modified Hertzsprung-Russel diagram (Fig. 2) displaying the versus plane, together with PARSEC isochrones in version 1.2S (Bressan et al. 2012). The stellar parameters and their uncertainties were derived by interpolating the sets of isochrones with ranging from 0.1 to 0.3.
The parameters and were checked for long timescale trends or periodic variability ( was not considered because it is correlated with ). The fitting procedure was repeated for a subset of 37 complete light curves, with and allowed to be determined independently for each dataset. The parameter was allowed to vary around a best-fit value under the Gaussian penalty, as given in Table 3. Individual determinations are plotted in Fig. 3. A model assuming a constant has a reduced chi-squared of 1.15 and a -value of 0.25. The parameter is also constant with that corresponds to a -value of 0.15. For both parameters, there is no reason to reject a null hypothesis assuming constancy. In addition, both parameters were searched for periodic variations with the Lomb-Scargle algorithm (Lomb 1976; Scargle 1982) and no significant signal was found. The lack of any variations justifies linking these parameters in the final model.
3.2 Transit timing
New mid-transit times, together with redetermined values from the literature data, allowed us to refine the transit ephemeris. As a result of a least-square linear fit, in which timing uncertainties were taken as weights, we derived
and the mid-transit time at cycle zero
The value of was found to be equal to 1.07, clearly showing that there is no hint for any deviation from the constant orbital period. The procedure was repeated for a limited sample of 17 mid-transit times with uncertainties smaller than 40 s, hence the most reliable. Again, the linear ephemeris was found to reproduce observed transits perfectly with . This shows that any possible TTV signal is not masked by a scatter of data points of lower quality and hence lower timing precision. Individual transit times and residuals are collected in Table 4, and residuals from the ephemeris are plotted in Fig. 4.
The Lomb-Scargle periodogram (Fig. 5) reveals no significant power at periods of and d, reported by von Essen et al. (2013). A bootstrap resampling method, based on trials of the randomly permuted timing residuals at the original observing epochs, was used to determine a false alarm probability (FAP) of 28% for the strongest peak.
|Epoch||(d)||O–C (d)||Data source|
3.3 Constraints on an additional planet
The TTV method is particularly sensitive to perturbers close to MMRs. The non-detection of periodic variations in transit timing allows such planetary configurations to be ruled out. In remaining configurations, planetary companions may be detected easier in the domain of precise Doppler measurements. Combining these two methods together places constraints on an upper mass of any hypothetical second planet in the system as a function of its semi-major axis.
We followed the method which we applied to WASP-3 and WASP-1 systems (Maciejewski et al. 2013, 2014). The Bulirsch–Stoer integrator, implemented in the Mercury 6 code (Chambers 1999), was used to predict deviations from a Keplerian motion of Qatar-1 b caused by a fictitious perturbing planet. The system was assumed to be coplanar, with both orbits initially circular. The mass of the perturber was set at 0.5, 1, 5, 10, 50, 100, and 500 (Earth masses), and the initial semi-major axis varied between 0.023 and 0.107 AU (corresponding to the orbital period ratio between 1 and 10) with a step of AU. Only configurations with outer perturber were considered because of the short orbital period of Qatar-1 b. The initial orbital longitude of the transiting planet was set at a value calculated for cycle zero, and the initial longitude of the fictitious perturber was shifted by . The calculations covered a time span of 1500 days in which transits of Qatar-1 b were observed. The residuals from a linear ephemeris for synthetic datasets were compared with the observed of transit timing, and a TTV upper mass of a fictitious planet was calculated.
Re-analysis of the RV dataset from Covino et al. (2013), performed with the Systemic 2.16 (Meschiari et al. 2009), gives a circular single-planet orbit with m s. The RV measurements were acquired with the High Accuracy Radial velocity Planet Searcher-North (HARPS-N) spectrograph installed at the 3.58 m Telescopio Nazionale Galileo (TNG) on La Palma (Spain). The RV uncertainties are between 4.4 and 20 m s. The dataset from Alsubai et al. (2011), acquired with the Tillinghast Reflector Echelle Spectrograph (TRES) coupled with the 1.5-m telescope at the Fred L. Whipple Observatory (USA), was omitted because of the much lower precision of the RV measurements ( m s). The value of was transformed into the RV mass limit of the hypothetical second planet on a circular orbit as a function of its semi-major axis.
Finally, the TTV and RV criteria were combined together to find planetary configurations with a hypothetical second planet below the detection threshold. Results are illustrated in Fig. 6. The TTV technique is more sensitive for configurations within the 2:1 resonance. We can eliminate perturbers with masses between 10 and 20 out of MMRs, and probe for a sub-Earth mass regime in 2:1 and 5:3 MMRs. For orbital periods longer than about 2.2 , the Doppler dataset provides tighter constraints except a 3:1 MMR, in which a perturber with a mass down to about 3 can be eliminated thanks to the TTV method.
4 Concluding discussion
Despite the conservative approach, our light curve modeling provides system parameters with the smallest uncertainties published so far. This is the result of the homogenous analysis of the rich set of 18 new and 31 literature light curves. The orbital inclination agrees perfectly with the value of reported by von Essen et al. (2013), and differs by 1.1 and 1.4 from results of Covino et al. (2013) and Alsubai et al. (2011), respectively. The scaled semi-major axis is between and brought by Covino et al. (2013) and von Essen et al. (2013), respectively. The planet to star radii ratio is consistent with derived by Alsubai et al. (2011), and differs by 1.5 and 3.5 from values reported by von Essen et al. (2013) and Covino et al. (2013), respectively. Planetary physical properties derived in this study confirm results of Covino et al. (2013), and deviate slightly from values reported by Alsubai et al. (2011) due to the planetary mass underestimated in that study. Stellar parameters were found to agree with values determined in previous studies. The system’s age is Gyr with the lower uncertainty of Gyr and the upper one limited by the age of the Universe. The host star is in the middle of its lifetime on the main sequence that is estimated to be Gy for a 0.8 star.
Spectral observations of the Ca II H and K lines show that the host star exhibits a moderate chromospheric activity with (Covino et al. 2013). However, we find no statistically significant variation in planet to star radii ratio that could be caused by the evolving distribution of dark spots or bright faculae outside the projected path of the planet on the stellar disk (e.g., Czesla et al. 2009). We also identified no features in the transit light curves which could be attributed to starspot occultations by the planetary disk.
Transit timing observations clearly show that the orbital motion of Qatar-1 b is not perturbed by any body which could produce periodic deviations with a range greater than 1 minute. We find no evidence which confirms the TTV signal with the periodicity and amplitude reported by von Essen et al. (2013). Combining the transit timing and RV datasets, we can rule out three proposed scenarios with perturbers in 2:1, 5:2, and 3:1 resonances even for circular orbits (in general, eccentric orbits make perturbations more pronounced). The lack of the TTV signal also makes a proposed massive perturber on a wide orbit unlikely. We note, however, that further RV measurements, acquired in a course of a few hundred days, would unequivocally reject this scenario.
We conclude that Qatar-1 b has no detectable planetary companions on nearby orbits. This finding is in line with results based on a sample of hot Jupiters observed with the space-borne facilities or ground based telescopes. As the formation of hot Jupiters is not yet fully understood (e.g., Steffen et al. 2012), the loneliness of Qatar-1 b speaks in favor of inward migration theories of massive outer planets through planet-planet scattering caused by mutual dynamical perturbations (Weidenschilling & Marzari 1996; Rasio & Ford 1996).
Acknowledgements.MK thanks the German national science foundation Deutsche Forschungsgemeinschaft (DFG) in projects NE515/34-1.34-2 and RE 882/12-2 for financial support. SR is currently a Research Fellow at ESA/ESTEC. SR would like to thank DFG for support in the Priority Programme SPP 1385 on the ”First Ten Million Years of the Solar System” in projects NE 515/33-1 and -2. RE thanks the Abbe-School of Photonics for support. AP acknowledges support from DFG in SFB TR7 project C2. JGS thanks DFG in SFB TR7 project B9 for support. We would like to acknowledge financial support from the Thuringian government (B 515-07010) for the STK CCD camera used in this project.
- thanks: Partly based on (1) data collected with telescopes at the Rozhen National Astronomical Observatory and (2) observations obtained with telescopes of the University Observatory Jena, which is operated by the Astrophysical Institute of the Friedrich-Schiller-University.
- email: email@example.com
- Diameter is the size of a telescope’s main mirror. is the total number of complete or partial transit light curves, acquired with the given instrument.
- Date UT is given for mid-transit time. Epoch is the transit number from the initial ephemeris given in Alsubai et al. (2011). and denote Cousins and Bessell -band filters, respectively. is the median number of exposures per minute. is the photometric scatter in millimag per minute as defined by Fulton et al. (2011). shows the airmass change during transit observations.
- Values for a solution with a circular orbit were taken.
- In calculations, the value of the RV semi-amplitude was taken from Covino et al. (2013).
- In calculations, the value of the RV semi-amplitude was taken from Covino et al. (2013).
- Based on parameters derived from transit photometry.
- Epoch is the transit number from the initial ephemeris given in Alsubai et al. (2011), is the mid-transit time in BJD minus 2450000 d. O–C is the value of timing residuals from the refined transit ephemeris. Data source: (1) von Essen et al. (2013); (2) Covino et al. (2013); (3) this paper.
- Agol, E., Steffen, J., Sari, R., & Clarkson, W. 2005, MNRAS, 359, 567
- Alsubai, K. A., Parley, N. R., Bramich, D. M., et al. 2013, Acta Astron., 63, 465
- Alsubai, K. A., Parley, N. R., Bramich, D. M., et al. 2011, MNRAS, 417, 709
- Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127
- Carter, J. A. & Winn, J. N. 2009, ApJ, 704, 51
- Chambers, J. E. 1999, MNRAS, 304, 793
- Claret, A. & Bloemen, S. 2011, A&A, 529, A75
- Covino, E., Esposito, M., Barbieri, M., et al. 2013, A&A, 554, A28
- Czesla, S., Huber, K. F., Wolter, U., Schröter, S., & Schmitt, J. H. M. M. 2009, A&A, 505, 1277
- Eastman, J., Gaudi, B. S., & Agol, E. 2013, PASP, 125, 83
- Eastman, J., Siverd, R., & Gaudi, B. S. 2010, PASP, 122, 935
- Fabrycky, D. C., Lissauer, J. J., Ragozzine, D., et al. 2014, ApJ, 790, 146
- Fulton, B. J., Shporer, A., Winn, J. N., et al. 2011, AJ, 142, 84
- Gazak, J. Z., Johnson, J. A., Tonry, J., et al. 2012, Advances in Astronomy, 2012, 30
- Holman, M. J. & Murray, N. W. 2005, Science, 307, 1288
- Kopal, Z. 1950, Harvard College Observatory Circular, 454, 1
- Lomb, N. R. 1976, Ap&SS, 39, 447
- Maciejewski, G., Niedzielski, A., Wolszczan, A., et al. 2013, AJ, 146, 147
- Maciejewski, G., Ohlert, J., Dimitrov, D., et al. 2014, Acta Astron., 64, 27
- Mandel, K. & Agol, E. 2002, ApJ, 580, L171
- Meschiari, S., Wolf, A. S., Rivera, E., et al. 2009, PASP, 121, 1016
- Mugrauer, M. & Berthold, T. 2010, Astronomische Nachrichten, 331, 449
- Müller, H. M., Huber, K. F., Czesla, S., Wolter, U., & Schmitt, J. H. M. M. 2013, A&A, 560, A112
- Podlewska, E. & Szuszkiewicz, E. 2009, MNRAS, 397, 1995
- Rasio, F. A. & Ford, E. B. 1996, Science, 274, 954
- Scargle, J. D. 1982, ApJ, 263, 835
- Steffen, J. H., Ragozzine, D., Fabrycky, D. C., et al. 2012, Proceedings of the National Academy of Science, 109, 7982
- Szabó, R., Szabó, G. M., Dálya, G., et al. 2013, A&A, 553, A17
- von Essen, C., Schröter, S., Agol, E., & Schmitt, J. H. M. M. 2013, A&A, 555, A92
- Weidenschilling, S. J. & Marzari, F. 1996, Nature, 384, 619