# Photometric follow-up of the transiting planetary system TrES-3:
transit timing variation and long-term stability of the system
^{†}^{†}thanks: Partly based on observations made at the Centro Astronómico
Hispano Alemán (CAHA), operated jointly by the Max-Planck Institut für Astronomie
and the Instituto de Astrofísica de Andalucía (CSIC).

###### Abstract

We present new observations of the transiting system TrES-3 obtained from 2009 to 2011 at several observatories. The orbital parameters of the system were re- determined and a new linear ephemeris was calculated. The best quality light curve was used for light curve analysis, and other datasets were used to determine mid-transit times () and study transit time variation (TTV). For planet parameter determination we used two independent codes and finally, we concluded that our parameters are in agreement with previous studies. Based on our observations, we determined 14 mid-transit times. Together with published we found that the timing residuals showed no significant deviation from the linear ephemeris. We concluded that a periodic TTV signal with an amplitude greater than 1 minute over a 4-year time span seems to be unlikely. Our analysis of an upper mass limit allows us to exclude an additional Earth-mass planet close to inner 3:1, 2:1, and 5:3 and outer 3:5, 1:2, and 1:3 mean-motion resonances. Using the long-term integration and applying the method of maximum eccentricity, the region from 0.015au to 0.05au was found unstable and the region beyond the 0.05au was found to have a chaotic behaviour and its depletion increases with increasing values of the initial eccentricity as well as inclination.

###### keywords:

planets and satellites: individual: TrES-3b â stars: individual: TrES-3^{†}

^{†}pagerange: Photometric follow-up of the transiting planetary system TrES-3: transit timing variation and long-term stability of the system

^{†}

^{†}thanks: Partly based on observations made at the Centro Astronómico Hispano Alemán (CAHA), operated jointly by the Max-Planck Institut für Astronomie and the Instituto de Astrofísica de Andalucía (CSIC). –References

^{†}

^{†}pubyear: 2012

## 1 Introduction

Since the first discovery of a transiting planet around HD 209458
(Charbonneau et al., 2000), 235 transiting extrasolar systems have already been
confirmed up to December 4th, 2012.^{1}^{1}1http://exoplanet.eu
Whilst transiting exoplanets offer unique scientific
possibilities, their study involves several complications.
In general, it is impossible to measure the mass
and radius of a planet based on a dataset obtained
with one observational technique. Transit light curves allow us to
determine just the relative
size of a star and planet, the orbital inclination and the stellar
limb-darkening coefficients. By combining this with radial-velocity measurements,
the observations offer the opportunity to measure the precise
stellar and planetary parameters. In order to obtain such parameters,
some constraints are needed, and are usually provided by forcing the
properties of the host stars to match theoretical expectations
(Southworth, 2010).
Significant uncertainties remain in the stellar mass and radius
determinations of many systems. In some cases, this is due to poorly
determined photospheric properties (i.e. effective temperature and metallicity),
and in other cases due to a lack of an accurate luminosity estimate
(Sozzetti et al., 2009). In addition, the different methods used for
these determinations as well as different approaches toward systematic
errors are leading to rather inhomogeneous set of planet properties.
Because of such inhomogenities, recently a few papers were published
where authors re-analyzed a large subset of known transiting planets,
and applied a uniform methodology
to all systems (e.g. Torres et al., 2008; Southworth, 2010, 2012).

In this paper we focus on the transiting system TrES-3. The system consists of a nearby G-type dwarf and a massive hot Jupiter with an orbital period of 1.3 days. It was discovered by O’Donovan et al. (2007) and also detected by the SuperWASP survey (Collier Cameron et al., 2007). Later, Sozzetti et al. (2009) presented new spectroscopic and photometric observations of the host star. A detailed abundance analysis based on high-resolution spectra yields [Fe/H] = 0.19 0.08, = 5650 75 K and log = 4.0 0.1. The spectroscopic orbital solution was improved with new radial velocity measurements obtained by Sozzetti et al. (2009). Moreover, these authors redetermined the stellar parameters (i.e. = 0.928 and = 0.829 ) and finally, the new values of the planetary mass and radius were determined (see Tab 3). They also studied the transit timing variations (TTVs) of TrES-3 and noted significant outliers from a constant period. In the same year, Gibson et al. (2009) presented the follow-up transit photometry. It consisted of nine transits of TrES-3, taken as part of transit timing program using the RISE instrument on the Liverpool Telescope. These transits, together with eight transit times published before (Sozzetti et al., 2009), were used to place upper mass limit as a function of the period ratio of a potential perturbing planet and transiting planet. It was shown that timing residuals are sufficiently sensitive to probe sub-Earth mass planet in both interior and exterior 2:1 resonances, assuming that the additional planet is in an initially circular orbit. Christiansen et al. (2011) has observed TrES-3 as a part of the NASA EPOXI Mission of Opportunity. They detected a long-term variability in the TrES-3 light curve, which may be due to star spots. They also confirmed that the planetary atmosphere does not have a temperature inversion. Later, Turner et al. (2013) observed nine primary transits of the hot Jupiter TrES-3b in several optical and near-UV photometric bands from June 2009 to April 2012 in an attempt to detect its magnetic field. Authors determined an upper limit of TrES-3b’s magnetic field strength between 0.013 and 1.3 G using a timing difference of 138 s derived from the Nyquist–Shannon sampling theorem. They also presented a refinement of the physical parameters of TrES-3b, an updated ephemeris and its first published near-UV light curve. The near-UV planetary radius of = 1.386 was also determined. This value is consistent with the planet’s optical radius. Recently, Kundurthy et al. (2013) observed eleven transits of TrES-3b over a two year period in order to constrain system parameters and look for transit timing and depth variations. They also estimated the system parameters for TrES-3b and found consistency with previous estimates. Their analysis of the transit timing data show no evidence for transit timing variations and timing measurements are able to rule out Super-Earth and Gas Giant companions in low order mean motion resonance with TrES-3b.

The main aims of this study can be summarized in the following items: (i) determination of the system parameters for TrES-3b (two independent codes will be used) and comparison with previous studies (i.e. O’Donovan et al., 2007; Sozzetti et al., 2009; Gibson et al., 2009; Colón et al., 2010; Southworth, 2010, 2011; Lee et al., 2011; Christiansen et al., 2011; Sada et al., 2012; Turner et al., 2013; Kundurthy et al., 2013). (ii) based on the obtained transits, we will determine the mid-transit times () and with following analysis of transit time variation (TTV) we will discuss possible presence of a hypothetical additional planet (perturber). We will try to estimate its upper-mass limit as a function of orbital periods ratio of transiting planet and the hypothetical perturber. (iii) Finally, using the long-term integration and applying the method of maximum eccentricity we will search for stability of regions inside the TrES-3b planet in context of additional planet(s).

The remainder of this paper is organized as follows. In the Section 2, we describe observations and data reduction pipelines used to produce the light curves. Section 3 presents the methods for analysis of transit light curves as well as discussion and comparison of the parameters of TrES-3 system. Section 4 and 5 are devoted to TTV and long-term stability of the system, respectively. Finally, in Section 6 we summarize and discuss our results.

## 2 Observations and data reduction

We obtained our data using several telescopes with different instruments. This allowed us to obtain many light curves since this strategy can effectively cope with the weather problems. On the other hand, this approach results in rather heterogeneous data. We used most of the data in average quality for the TTV analysis which is not very demanding on homogeneity of the data. Only the best quality light curve was used for the planet parameter determination.

Observations used in this paper were carried out at the several observatories in Slovakia (Stará Lesná Observatory; 49°09’ 10”N, 20°17’ 28”E), Poland (Piwnice Observatory; 53°05’ 43”N, 18°13’ 46”E), Germany (Grossschwabhausen Observatory; 50°55’ 44”N, 11°29’ 03”E; Volkssternwarte Kirchheim Observatory, 50°55’44”N, 11°29’ 03”E and Michael Adrian Observatory, 49°55’27”N, 08°24’ 33”E) and Spain (Calar Alto Observatory; 37°13’25”N, 02°32’ 46”E). We collected 14 transit light curves obtained between May 2009 and September 2011. The transits on May 12, 2009 and August 20, 2010 were observed simultaneously at two different observatories. The telescope diameters of 0.5 to 2.2 m allowed us to obtain photometry with 1.2 – 7.8 mmag precision, depending on observing conditions. Observations generally started 1 hour before the expected beginning of a transit and ended 1 hour after the event. Unfortunately, weather conditions and schedule constraints meant that we were not able to fit this scheme in all cases.

All instruments are equipped with CCD cameras with the Johnson-Cousins
()
standard filter system.
The information from individual observatories and instruments as well as the
summary of observing runs are given in Table 1 and Table
2. The standard correction procedure
(bias, dark and flat field correction) and subsequently aperture photometry was
performed by IRAF^{2}^{2}2IRAF is distributed by the National
Optical Astronomy Observatories, which are operated by the Association of Universities for
Research in Astronomy, Inc., under cooperative agreement with the National Science
Foundation. and task chphot (Raetz et al., 2009) (GSH and VK),
C-munipack package^{3}^{3}3http://c-munipack.sourceforge.net/ (G1) and
Mira_Pro_7^{4}^{4}4http://www.mirametrics.com/mira_pro.htm (MA). Data from
remaining telescopes (P and CA) were reduced with the software pipeline
developed for the Semi–Automatic Variability Search sky survey (Niedzielski, Maciejewski &
Czart, 2003).
To generate an artificial comparison star, at least 20–30 per cent of stars
with the lowest light–curve scatter were selected iteratively from the field
stars brighter than 2.5–3 mag below the saturation level
(e.g. Broeg et al., 2005; Raetz et al., 2009). To measure
instrumental magnitudes, various aperture radii were used. The aperture which
was found to produce light curve with the smallest overall scatter was applied to
generate final light curve. The linear trend in the out-of-transit
parts was also removed.

Obs. | Telescope | Detector | N |
---|---|---|---|

CCD size | FoV [arcmin] | ||

G1 | Newton | SBIG ST10-MXE | 5 |

508/2500 | 2184 1472, 6.8 m | 20.4 13.8 | |

MA | Cassegrain | SBIG STL-6303E | 2 |

1200/9600 | 3072 2048, 9 m | 10 7 | |

GSH | CTK | SITe TK1024 | 1 |

250/2250 | 1024 1024, 24 m | 37.7 37.7 | |

STK | E2V CCD42-10 | 1 | |

600/1758 | 2048 2048, 13.5 m | 52.8 52.8 | |

VK | RCT | STL-6303E | 2 |

600/1800 | 3072 2048, 9 m | 71 52 | |

P | Cassegrain | SBIG STL-1001 | 1 |

600/13500 | 1024 1024, 24m | 11.8 11.8 | |

CA | RCF | SITe CCD | 2 |

2200/17037 | 2048 2048, 24m | 18.1 18.1 |

Obs. | Date | Filter | (s) | |
---|---|---|---|---|

G1 | 2009 May 12 | 319 | 40 | |

2009 Aug 01 | 345 | 45 | ||

2010 Apr 27 | 180 | 35 | ||

2010 Jun 30 | 238 | 40 | ||

2010 Aug 07 | 168 | 35 | ||

MA | 2010 July 13 | 349 | 25 | |

2010 Aug 20 | 296 | 20 | ||

GSH (CTK) | 2009 May 25 | 138 | 80 | |

(STK) | 2011 Mar 22 | 131 | 60 | |

VK | 2009 Aug 14 | 103 | 120 | |

2010 Aug 20 | 320 | 30 | ||

P | 2009 May 12 | 120 | 55 | |

CA | 2010 Sep 06 | 158 | 30-35 | |

2011 Sep 12 | 128 | 45 |

## 3 Light curve analysis

Our photometric observations of the TrES-3 system consist of data from different instruments and are of different photometric quality. For the purpose of radius determination, we decided to analyse only a light curve with the lowest scatter and the best rms = 1.2 mmag. We choose data obtained at Calar Alto observatory on September 06, 2010 (see Figs. 1 and 2). We first refined the light-curve system parameters and subsequently determined the individual times of transit as described below in order to improve ephemeris.

For the calculation of synthetic light curves we used two independent approaches: the first one is based on the routines from (Mandel & Agol, 2002) and the Monte Carlo simulations (described in Section 3.1), the second one uses JKTEBOP code (Southworth, Maxted & Smalley, 2004) (see Section 3.2).

### 3.1 Solution 1

First we used the downhill simplex minimization procedure (implemented in routine AMOEBA, Press et al., 1992) to determine system parameters (planet to star radius ratio), (inclination), (mid-transit time) and (star radius to semi-major axis ratio). The model light curve itself was computed via the analytic expressions from Mandel & Agol (2002). The quadratic limb darkening law was assumed and corresponding limb darkening coefficients , were linearly interpolated from Claret (2000) assuming the stellar parameters from Sozzetti et al. (2009): = 5 650 K, (g) = 4.4 and [Fe/H] = 0.19. As a goodness of the fit estimator we used the function:

(1) |

where is the model value and is the measured value of the flux, is the uncertainty of the measurement and the sum is taken over all measurements. The orbital period and the limb darkening coefficients were fixed through the minimization procedure. The transit duration was determined assuming the semi-major axis = 0.02282au (Sozzetti et al., 2009). The parameters and the corresponding light curve for which we found the minimum value of function, are in Table 3 and Figure 2, respectively (SOLUTION 1).

To estimate the uncertainties of the calculated transit parameters, we employed the Monte Carlo simulation method (Press et al., 1992). We produced 10 000 synthetic data sets with the same probability distribution as the residuals of the fit in Figure 2. From each synthetic data set obtained by this way we estimated the synthetic transit parameters. The minimum value corresponding to each set of Monte Carlo parameters was calculated as:

(2) |

where is the original best fit model value and is the ’Monte–Carlo simulated’ value. Figure 3 shows the dependence of the parameters and on the reduced . This quantity is defined as:

(3) |

where is the number of data points and is the number of fitted parameters. To fully understand the errors of the system parameters we constructed confidence intervals in 2-dimensional space. Figure 4 depicts the confidence regions for 2 parameters ( vs. ) as a projection of the original 4-dimensional region. The gray-colored data points stand for the , and region with corresponding value of , and , respectively ( is the of the original best fit model value). From the shape of the dependence it could be seen that these two parameters correlate (see Section 3.3).

Finally, we took into account the uncertainty of the stellar mass and semi-major axis according to the simple error propagation rule. The results from the first analysis are shown in the Table 3 as SOLUTION 1.

### 3.2 SOLUTION 2 (JKTEBOP code)

We have modelled the light curve using the JKTEBOP^{5}^{5}5JKTEBOP
is written in FORTRAN77 and the source code is available at
http://www.astro.keele.ac.uk/jkt/codes/jktebop.html code
(Southworth, 2011) as well.
JKTEBOP grew out of the original EBOP program written for eclipsing binary
systems (Etzel, 1981; Popper & Etzel, 1981) via implementing the NDE (Nelson & Davis, 1972)
model. JKTEBOP uses biaxial spheroids to model the component stars (or star and planet)
and performs a numerical integration in concentric annuli over the surface of each
body to obtain the flux coming from the system. This feature of the code allows
us to avoid
small and spherical planet approximations which are used in analytic light-curve
generators based on Mandel & Agol (2002), and hence to derive planetâs oblateness.
A model is fitted to the data by the Levenberg-Marquardt least-square procedure.
The code converges rapidly toward a reliable solution and diminishes the correlation
between fitted parameters (Southworth, 2008).

The main parameters of a JKTEBOP fit are the orbital inclination , and the fractional radii of the host star and planet, and . The fractional radii are defined as:

(4) |

where and are the stellar and planetary radii and is the orbital semi-major axis. Parameters and correspond to radii of spheres of the same volume as the biaxial spheroids. In JKTEBOP the fractional radii are reparametrized as their sum and ratio:

(5) |

because these are only weakly correlated with each other (Southworth, 2008). The directly fitted orbital inclination, , allows the transit impact parameter = cos to be calculated. The initial values of parameters listed above were taken from Sozzetti et al. (2009). Because of different quality of data, the synthetic light curve was calculated only for the best light curve obtained at Calar Alto on September 06, 2010 (the same like in SOLUTION 1). A value of = 0.02282 au (Sozzetti et al., 2009) was used in subsequent calculations. The best-fit model was used as a template and fitted to other light curves for which only the mid-transit time was allowed to vary. The resulting mid-transit times together with obtained from literature (Sozzetti et al., 2009; Gibson et al., 2009) are analysed in details in the Section 4. The determined parameters obtained by JKTEBOP code are presented in Table 3 as SOLUTION 2.

The errors of derived parameters were determined in two ways for each combination of data set and adopted LD law (we have used the same coefficients like in the case of SOLUTION 1). Firstly, we ran 1000 Monte Carlo (MC) simulations, a spread range of a given parameter within 68.3% was taken as its error estimate. Secondly, the prayer-bead method (e.g. Désert et al., 2011; Winn et al., 2009) was used to check whether red noise was present in our data. MC errors were found to be 2 – 3 times smaller than the values returned by the prayer bead method. This indicates that the light curve is affected not only by Poisson noise but also by additional correlated noise. Therefore, our prayer-bead error estimates were taken as our final errors.

Source | |||||
---|---|---|---|---|---|

[min] | [days] | ||||

O’Donovan et al. (2007) | 0.1650 0.0027 | 0.1660 0.0024 | 82.15 0.21 | – | 1.30619(1) |

Sozzetti et al. (2009) | 0.1687 | 0.1655 0.0020 | 81.85 0.16 | – | 1.30618581(51) |

Gibson et al. (2009) | – | 0.1664 | 81.73 | 79.92 | 1.3061864(5) |

Colón et al. (2010) | – | 0.1662 | – | 83.77 | – |

Southworth (2010) | – | – | 82.07 0.17 | – | 1.3061864(5) |

Lee et al. (2011) | – | 0.1603 0.0042 | 81.77 0.14 | – | 1.30618700(15) |

Christiansen et al. (2011) | 0.1664 0.0204 | 0.1661 0.0343 | 81.99 0.30 | 81.9 1.1 | 1.30618608(38) |

Southworth (2011) | – | – | 81.93 0.13 | – | 1.30618700(72) |

Sada et al. (2012) | – | – | – | 77.9 1.9 | 1.3061865(2) |

Kundurthy et al. (2013) | |||||

Solution_1 | 0.1675 0.0008 | 0.1652 0.0009 | 81.95 0.06 | – | 1.3062132(2) |

Kundurthy et al. (2013) | |||||

Solution_2 | 0.1698 0.0014 | 0.1649 0.0015 | 81.51 0.14 | – | 1.3062128(2) |

Turner et al. (2013) | 0.1721 | 0.1693 | 81.35 | 81.30 0.23 | 1.306 1854(1) |

This work | |||||

Solution_1 | 0.1682 0.0032 | 0.1644 0.0047 | 81.86 0.28 | 79.20 1.38 | 1.306186 |

This work | |||||

Solution_2 | 0.1696 | 0.1669 | 81.76 | 79.08 0.72 | 1.306186 |

### 3.3 Light curve analysis results

The resulting values of the parameters together with their uncertainties are given in Table 3. For comparison, the parameters from previous studies are added there as well. Figure 2 shows the resulting best fit obtained by SOLUTION 1. In order to compare both solutions, we also plotted differences from the fit obtained by routines resulting from SOLUTION 1 and SOLUTION 2.

The final parameters are in good agreement with already published values. The output from SOLUTION 1, in particular and , correspond a bit better to parameters of Sozzetti et al. (2009) and Christiansen et al. (2011). The radii determined by JKTEBOP are also within the range of errors. Based on these two parameters we also inferred the critical value of inclination for total transit of TrES-3b as:

(6) |

As listed in Table 3, all values of inclination including uncertainties are in agreement with critical inclination calculated above. This is an evidence of grazing transit of TrES-3b.

We used MC simulations to produce of 10000 synthetic data sets with the same probability distribution as the residuals of the fit. From each synthetic data set we estimated the synthetic transit parameters. Using the results of the simulation, the dependence of the inclination versus the planet to star radius ratio is plotted in Figure 4. The figure demonstrates that there is a correlation between the parameters and that solution for the system TrES-3 is not unique and can be located in a relatively wide range (degeneracy of the parameters). This fact is also caused by grazing transit and by subsequent sensitivity of the solution to the determination of LD coefficients.

For determination of our mid-transit times, the best-fit model was used as a template and fitted to other light curves for which only the mid-transit time was allowed to vary. In order to have homogeneous analysis of TTV, we collected the data points of transit light curves published in previous studies (see Figure 5). Transit light curve data of Sozzetti et al. (2009) and Colón et al. (2010) were available on-line and the remaining tabulated data were obtained from mentioned authors by e-mail. We re-analysed all transit light curves and derived mid-transit times by the same procedure (the code JKTEBOP) to get a homogeneous dataset for TTV analysis (see Section 4).

We also re-analysed all light curves under consideration with as a free parameter and tried to search for any variation. All determinations of were consistent within error bars with a mean value, so we did not detect any transit depth variation. We also saw no significant signal in periodograms.

## 4 Transit Timing Variation

Our new 14 mid-transit times and 42 redetermined literature values
were used to refine the transit ephemeris.
The mid-transit times were transformed from JD or HJD (based on UTC)
into BJD (based on Barycentric Dynamical Time – TDB) using the on-line
converter^{6}^{6}6http://astroutils.astronomy.ohio-state.edu/time/utc2bjd.html by Eastman, Siverd & Gaudi (2010).
As a result of fitting a linear function of the epoch, we obtained the mid-transit time for the initial
epoch BJD and the orbital period d.
The individual mid-transit errors were taken as weights. The linear fit yields reduced of that is
similar to a value of 2.3 reported by Gibson et al. (2009). These values, noticeable greater than 1 might suggest the
existence of an additional planet which perturbs the orbital motion of TrES-3 b (Sozzetti et al., 2009; Gibson et al., 2009).

Date | Obs. | Epoch | error | (d) | |
---|---|---|---|---|---|

2009 May 12 | P | 326 | 2454964.39885 | 0.00084 | |

2009 May 12 | G1 | 326 | 2454964.39859 | 0.00081 | |

2009 May 25 | GSH (CTK) | 336 | 2454977.4593 | 0.0016 | |

2009 Aug 01 | G1 | 388 | 2455045.38091 | 0.00045 | |

2009 Aug 14 | VK | 398 | 2455058.44418 | 0.00077 | |

2010 Apr 27 | G1 | 594 | 2455314.45510 | 0.00062 | |

2010 Jun 30 | G1 | 643 | 2455378.45937 | 0.00031 | |

2010 July 13 | MA | 653 | 2455391.52067 | 0.00017 | |

2010 Aug 07 | G1 | 672 | 2455416.3397 | 0.0011 | |

2010 Aug 20 | MA | 682 | 2455429.40118 | 0.00033 | |

2010 Aug 20 | VK | 682 | 2455429.3997 | 0.0004 | |

2010 Sep 06 | CA | 695 | 2455446.38075 | 0.00017 | |

2011 Mar 22 | GSH (STK) | 846 | 2455643.6145 | 0.0005 | |

2011 Sep 12 | CA | 979 | 2455817.3372 | 0.0014 |

Results for new mid-transit times are shown in Table 4. The observed minus calculated (O–C) diagram, plotted in Fig. 5, shows no significant deviation from the linear ephemeris. All our data points deviated by less than . We also searched for a periodicity that could be a sign of an additional body in the system. We generated a Lomb–Scargle periodogram (Lomb, 1976; Scargle, 1982) for the residuals in a frequency domain limited by the Nyquist frequency and found the highest peak at 30 d and peak-to-peak amplitude of s. This period could coincide with a stellar rotation period which is roughly estimated to be 28 d. Examples of such TTV signals induced by the stellar activity are observed in the Kepler data (e.g. Mazeh et al., 2013). However, the false alarm probability (calculated empirically by a bootstrap resampling method with trials) of the putative signal is disqualifying with a value of 18.2%. In addition, the amplitude is close to the mean 1- timing error of our observations. These findings allow us to conclude that a strictly periodic TTV signal with the amplitude greater than 1 minute over a 4-year time span seems to be unlikely.

Following Gibson et al. (2009), we put upper constraints on a mass of a potential perturbing planet in the system with refined assumptions. We adopted a value of 1 min for the maximal amplitude of the possible TTV signal. We also increased the sampling resolution to probe resonances other than inner and outer 1:2 commensurabilities studied in detail by Gibson et al. (2009). We simplified the three-body problem by assuming that the planetary system is coplanar and initial orbits of both planets are circular. The masses of both the star and the transiting planet, as well as its semi-major axis were taken from Sozzetti et al. (2009). The orbital period of the hypothetical perturber was varied in a range between 0.2 and 5 orbital periods of TrES-3 b. We produced 2800 synthetic O–C diagrams based on calculations done with the Mercury code (Chambers, 1999). We applied the Bulirsch–Stoer algorithm to integrate the equations of motion.

The most important feature of the Bulirsch-Stoer algorithm for N-body simulations is that it is capable of keeping an upper limit on the local errors introduced due to taking finite time-steps by adaptively reducing the step size when interactions between the particles increase in strength. Calculations covered 1500 days, i.e. a total time span of available transit observations. The results of simulations are presented in Fig. 6. Our analysis allows us to exclude an additional Earth-mass planet close to inner 1:3, 1:2, and 3:5 and outer 5:3, 2:1, and 3:1 mean-motion resonances (MMRs).

## 5 Long-term stability of the system

In this section, we investigated the long-term gravitational influence of TrES-3b planet on a potential second planet in the system. Thus, we performed numerical simulation for studying the stability of orbits and checking their chaotic behavior using the method of maximum eccentricity (e.g. Dvorak et al., 2003).

We used long-term integration of small-mass (Earth-mass) planet orbits for inspecting the stability regions in the TrES-3b system. The time span of the integration was around 140000 revolutions of the planet around the star. For the integration of this system we again applied an efficient variable-time-step algorithm (Bulirsch-Stoer integration method). The parameter which controls the accuracy of the integration was set to , in this case.

We have generated mass-less particles for representing small planets in this system. We assumed that semi-major axis ranges from 0.00625au to 0.1au. The inner border of the generated system of small planets is around 1.5 times star radius. This value of inner border is approximately at the location of Roche limit for a Earth-like planet and includes also the region in the vicinity of 3:1 MMR (0.01 au for this case) in our analysis. The step-size in semi-major axis was = 0.00025au, in eccentricity = 0.02 and = 5 in inclination. The upper limit of the grid in eccentricity was 0.5 and 50 in inclination. We integrated the orbits of the small planets for the time-span of 500 years (about 140 000 revolutions of TrES-3b around the parent star). We obtained orbital evolution of each small planet in the system and also were able to study the stability regions at the end of the integration using the method of maximum eccentricity, mentioned above. The maximum eccentricity that the potential perturber of TrES-3b reached over the time of the integration is plotted on the () stability map in the Fig. 7 (Fig. 8).

Figure 7 shows the stability maps in plane for selected values of inclinations: , , and (from top left to bottom right). Also the MMRs with TrES-3b planet are marked in these plots. We can see a stable region inside the 2:1 MMR for each value of inclination. Outside the 2:1 MMR the gravitational influence of the TrES-3b planet is very strong and leads to depletion of these regions. Only the MMRs 2:1, 3:2 and 1:1 are moderately populated but the population decreases with the increase of the inclination. For completeness, we note that 5:3 ( 0.016 au) and 3:5 ( 0.032 au) MMRs are depleted and thus not stable for additional Earth-mass planet in their vicinity. In Fig. 8 we present the stability maps in plane for several values of eccentricities: , , and . One can see the stable region inside the 2:1 MMR with TrES-3b planet, which we found stable also in the plane. Considering the region beyond the 1:3 MMR, the depletion of the planet population is not so strong than in the region between 2:1 and 1:3 MMR. This feature is in a good agreement with the weaker gravitational influence of the TrES-3b planet.

Based on our results, we showed that the region inside the TrES-3b planet orbit (especially inside the 2:1 up to 3:1 MMR) can be stable on longer timescales (hundreds of years or hundred thousands of revolutions of TrES-3b around the star) from the dynamical point of view. The region from 0.015au (near 2:1 MMR) to 0.05au (near 1:3 MMR) was found unstable apart from moderately populated MMRs located in this area (see figures Fig. 7 and Fig. 8). The relatively small increase of population in the mentioned MMRs depends on the initial values of semi-major axis, eccentricity and inclination. The region beyond the 0.05au was found to have a chaotic behavior and the depletion of the planet population increases with increasing values of initial eccentricity as well as inclination.

## 6 Discussion and Conclusion

Based on the transit light curves obtained at several observatories between May 2009 and September 2011 we redetermined orbital parameters and the radius of the transiting planet TrES-3b. The best light curve (obtained at Calar Alto, Sep 2010) was used for light curve analysis, and the data from the other observatories were used for determination and TTV investigation. We used two independent solutions for parameters determination and finally, we concluded that our values are consistent with previous results of Sozzetti et al. (2009) and Christiansen et al. (2011).

The aim of this present paper was also to discuss possible presence of a second planet in the extrasolar system TrES-3. For this purpose we used our new 14 mid-transit times and the individual determinations from Sozzetti et al. (2009) and Gibson et al. (2009). The resulting diagram showed no significant deviation of data points from the linear ephemeris. In addition, we tried to search for a periodicity that could be caused by an additional body in the system. We can conclude that a strictly periodic TTV signal with the amplitude greater than 1 minute over a 4-year time span seems to be unlikely. This result, together with refined assumptions of Gibson et al. (2009) allow us to put upper constraints on the mass of a potential perturbing planet. The additional Earth-mass planet close to inner 3:1, 2:1, and 5:3 and outer 3:5, 1:2, and 1:3 MMRs can be excluded.

Finally, we used the long-term integration of the theoretical set of massless particles generated in TrES-3 system for studying the dynamical stability of potential second planet in the system (influenced by TrES-3b gravitation). From our analysis we found that the region inside the TrES-3b planet orbit (especially inside the 2:1 MMR) up to 3:1 MMR can be stable on longer timescales (hundreds of years or hundred thousands of revolutions of TrES-3b around the star) from the dynamical point of view. The region from 0.015au (near 2:1 MMR) to 0.05au (near 1:3 MMR) was found to be unstable apart from moderately populated MMRs located in this area. The relatively small increase of population in these MMRs depends on the initial values of semi-major axis, eccentricity and inclination. The region beyond 0.05au was found to have a chaotic behavior and depletion of the planet population increases with increasing values of initial eccentricity as well as inclination.

## Acknowledgements

This work has been supported by the VEGA grants No. 2/0078/10, 2/0094/11, 2/0011/10. MV, MJ, JB, TP and ŠP would like to thank also the project APVV-0158-11. TK thanks the Student Project Grant at MU MUNI/A/0968/2009 and the National scholarship programme of the Slovak Republic. GM acknowledges Iuventus Plus grants IP2010 023070 and IP2011 031971. SR thanks the German National Science Foundation (DFG) for support in project NE 515 / 33-1.

## References

- Broeg et al. (2005) Broeg C., Fernández M., Neuhäuser R. 2005, AN, 326, 134
- Chambers (1999) Chambers J.E. 1999, MNRAS, 304, 793
- Charbonneau et al. (2000) Charbonneau D., Brown T. M., Latham D. W., Mayor M. 2000, ApJ, 529, L45
- Christiansen et al. (2011) Christiansen J. L., Ballard S., Charbonneau D. 2011, ApJ, 726, 94
- Claret (2000) Claret A. 2000, A&A, 363, 1081
- Collier Cameron et al. (2007) Collier Cameron A., Wilson D. M., West R. G., et al. 2007, MNRAS, 380, 1230
- Colón et al. (2010) Colón K. D., Ford E. B., Lee B., et al. 2010, MNRAS, 408, 1494
- Désert et al. (2011) Désert J. M., Sing D., Vidal-Madjar A., et al. 2011, A&A, 526, 12
- Dvorak et al. (2003) Dvorak, R., Pilat-Lohinger, E., Funk, B., Freistetter, F. 2003, A&A, 398, L1
- Eastman et al. (2010) Eastman J., Siverd R., Gaudi B.S. 2010, PASP, 122, 935
- Etzel (1981) Etzel P. B. 1981, in Carling E. B., Kopal Z., eds, NATO ASI Ser. C., 69, Photometric and Spectroscopic Binary Systems, Kluwer, Dordrecht, 111
- Gibson et al. (2009) Gibson N. P., Pollacco D., Simpson E. K., et al. 2009, ApJ, 700, 1078
- Kundurthy et al. (2013) Kundurthy P., Becker A. C., Agol E., et al. 2013, ApJ, 764, 8
- Lee et al. (2011) Lee J. W., Youn J. H., Kim S. L., et al. 2011, PASJ, 63, 301
- Lomb (1976) Lomb N.R. 1976, Ap&SS, 39, 447
- Mazeh et al. (2013) Mazeh T., Nachmani G., Holczer T., et al. 2013, arXiv:1301.5499v1
- Mandel & Agol (2002) Mandel K., Agol E. 2002, ApJ, 580, L171
- Mugrauer (2009) Mugrauer M. 2009, AN, 330, 419
- Mugrauer & Berthold (2010) Mugrauer M., Berthold T. 2010, AN, 331, 449
- Nelson & Davis (1972) Nelson B., Davis W. D. 1972, ApJ, 174, 617
- Niedzielski, Maciejewski & Czart (2003) Niedzielski A., Maciejewski G., Czart K. 2003, Acta Astronomica, 53, 281
- O’Donovan et al. (2007) O’Donovan F. T., Charbonneau D., Bakos G. Á., et al. 2007, ApJ, 663, L37
- Popper & Etzel (1981) Popper D. M., Etzel P. B. 1981, AJ, 86, 102
- Press et al. (1992) Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B. P. 1992, Numerical recipes in FORTRAN, Cambridge: University Press, London
- Raetz et al. (2009) Raetz S., Mugrauer M., Schmidt T. O. B., et al. 2009, AN, 330, 475
- Sada et al. (2012) Sada P. V., Deming, D. Jennings D. E., et al. 2012, PASP, 124, 212
- Scargle (1982) Scargle J.D. 1982, ApJ, 263, 835
- Southworth, Maxted & Smalley (2004) Southworth J., Maxted P. F. L., Smalley B. 2004, MNRAS, 349, 547
- Southworth (2008) Southworth J. 2008, MNRAS, 386, 1644
- Southworth (2010) Southworth J. 2010, MNRAS, 408, 1689
- Southworth (2011) Southworth J. 2011, MNRAS, 417, 2166
- Southworth (2012) Southworth J. 2012, MNRAS, 426, 1291
- Sozzetti et al. (2009) Sozzetti A., Torres G., Charbonneau D., et al. 2009, ApJ, 691, 1145
- Torres et al. (2008) Torres G., Winn J. N., Holman M. 2008, ApJ, 677, 1324
- Turner et al. (2013) Turner J. D., Smart B. M., Hardegree-Ullman K. K. et al. 2013, MNRAS, 428, 678
- Winn et al. (2009) Winn J. N., Holman M. J., Henry G. W., et al. 2009, ApJ, 693, 794