Article

Multifilter Transit Observations of WASP-39b and WASP-43b with Three San Pedro Mártir Telescopes

, , , , , , , , and

© 2015. The Astronomical Society of the Pacific. All rights reserved. Printed in U.S.A.
, , Citation D. Ricci et al 2015 PASP 127 143 DOI 10.1086/680233

1538-3873/127/948/143

Abstract

Three optical telescopes located at the San Pedro Mártir National Observatory were used for the first time to obtain multifilter defocused photometry of the transiting extrasolar planets WASP-39b and WASP-43b. We observed WASP-39b with the 2.12 m telescope in the U filter for the first time, and additional observations were carried out in the R and I filters using the 0.84 m telescope. WASP-43b was observed in VRI with the same instrument, and in the i filter with the robotic 1.50 m telescope. We reduced the data using different pipelines and performed aperture photometry with the help of custom routines in order to obtain the light curves. The fit of the light curves (1.5–2.5 mmag rms), and of the period analysis, allowed a revision of the orbital and physical parameters, revealing for WASP-39b a period (4.0552947 ± 9.65 × 10-7 days) which is 3.084 ± 0.774 seconds larger than previously reported. Moreover, we find for WASP-43b a planet/star radius (0.1738 ± 0.0033) which is 0.01637 ± 0.00371 larger in the i filter with respect to previous works, and that should be confirmed with additional observations. Finally, we confirm no evidence of constant period variations in WASP-43b.

Export citation and abstract BibTeX RIS

1. Introduction

Since the first detection of an exoplanet orbiting around a pulsar (Wolszczan & Frail 1992; Wolszczan 1994) and around a main-sequence star (Mayor & Queloz 1995) using radial velocity measurements, several other techniques, such as microlensing, coronagraphy, and polarimetric observations have been developed for the search of other worlds. Transit observations have contributed to a large fraction of the about 1800 confirmed planets. 1 Several systematic surveys have been implemented for the detection of exoplanet candidates by transit observations.

Important examples of these surveys include space-based missions such as COROT (Baglin et al. 2009) and Kepler (Borucki et al. 2010) satellites. The high photometric precision of their observations boosted the number of candidate and confirmed extrasolar planets detected in the last 5 years. The next generation of space-based missions, such as the orbiting Gaia (Baglin & Catala 2009) and the upcoming James Webb Space Telescope (JWST; Clampin 2008; Batalha et al. 2014), Planetary Transits and Oscillations of Stars (PLATO; Rauer et al. 2013), and Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2010) are expected to increase the detection rate and broaden the limits of the currently observed distribution of exoplanetary parameters, despite that some of these missions do not have the detection of exoplanets as a main goal.

On the other hand, ground-based photometric surveys such as SuperWASP (Pollacco et al. 2006) and HAT (Bakos et al. 2007) are focused on exoplanets with large radii and orbiting close to their host star, as these objects are particularly suited for these instruments. In order to improve the photometric precision of ground-based follow-up observations of exoplanetary transits, a technique consisting of slightly defocusing the image has been introduced and widely used in recent works (Southworth et al. 2009a, 2009b, 2010, 2012b, 2013, 2014). The achieved results often compete with those of space-based surveys, and increase the opportunity for small telescopes to produce light curves with a remarkable photometric precision.

Following the results obtained with the defocused photometry technique, we started an observational campaign of known exoplanetary transits, already detected by the SuperWASP and HAT surveys. Photometric follow-up of these objects in different bandpasses is useful to study a possible dependence of the planetary radius as a function of the wavelength.

This project makes use of all three telescopes operating at the San Pedro Mártir National Astronomical Observatory (see § 2 for details), which we employ for the first time in the field of transit observations. In this article, we focus on two particular targets: WASP-39b and WASP-43b, which are briefly presented in the following sections.

1.1. WASP-39b

WASP-39b is a Saturn-mass planet recently discovered by the SuperWASP survey, and confirmed by Faedi et al. (2011) with a photometric and spectroscopic follow-up observations. The photometry was carried on with the Pan-STARRS-z filter on the 2 m Faulkes Telescope North, and in the Gunn-i filter on the 1.2 m Euler Telescope. The spectroscopy was carried out with the SOPHIE (Perruchot et al. 2008) and the CORALIE (Baranne et al. 1996; Queloz et al. 2000) spectrographs, installed on the 1.93 m telescope at the Haute-Provence Observatory (OHP) and on the 1.2 m Euler Telescope, respectively.

Faedi et al. (2011) obtain a period of P = 4.055259 days, a T0 = 2455342.9688 in Heliocentric Julian Date (HJD), a transit duration of 168.192 minutes, and the following physical parameters: a mass of 0.28 MJ , a radius of 1.27 RJ , and thus a resulting mean density of 0.14ρJ . This is one of the lowest values for the currently known extrasolar planets. The reported semimajor axis is 0.0486 AU, and evidences of a highly inflated radius are also found.

This work presents the first photometric follow-up of WASP-39b since the discovery.

1.2. WASP-43b

WASP-43b is a very short-period, hot Jupiter-like planet discovered by Hellier et al. (2011) using observations from the WASP-South camera array, complemented with photometric observations with the 0.6 m TRAnsiting Planets and PlanetesImals Small Telescope (TRAPPIST) in the I + z band and the Euler Telescope in the Gunn r band. Radial velocity measurements were also obtained with CORALIE (Echelle Spectrograph on the 1.2-m Leonard Euler Swiss Telescope). Hellier et al. (2011) find a period of P = 0.813475 days, an initial epoch of T0 = 2455528.86774 HJD, and a duration of 69.552 min. The authors derive a mass of 1.8 MJ , a radius of 0.9 RJ leading to a density of 2.21ρJ . The reported semimajor axis is 0.0142 AU.

In order to improve the characterization of WASP-43b and its host star, for which Husnoo et al. (2012) find evidences of excess rotation, Gillon et al. (2012) observed over 30 eclipses between transits and occultations, using mainly TRAPPIST in the I + z and z bands, respectively. The survey was completed with transit observations obtained with the Euler Telescope in the Gunn r' band and with the VLT/HAWK-I instrument installed on the 8.2 m UT4 telescope in Paranal. The work allowed significant improvements in the determination of the stellar density and in the refinement of the parameters of the planetary system. In particular, Gillon et al. (2012) find a mass of 2.034 MJ and a radius of 1.036 RJ . Occultations also detected the thermal emission of the planet.

Thermal emission is also the subject of Spitzer observations, carried out by Blecic et al. (2012, 2014), which highlighted the possibility of thermal inversion in the atmosphere.

Line et al. (2014) uniformly analyzed occultations of nine exoplanets including WASP-43b, with the goal to characterize their C to O ratio, which is found to be >1 in the case of WASP-43b, and which is attributed to the upper limit value of CH4.

Other observations using the GROND (Gamma-Ray Burst Optical/Near-Infrared Detector) instrument on the MPG/ESO 2.2 m telescope in the g', r', i', z', J, H, and K bands (Chen et al. 2014) allowed to refine the period to a value of P = 0.81347437. This value agrees with spectroscopic observations carried on by Murgas et al. (2014). These observations also confirmed planetary dayside thermal emission with poor heat redistribution, as observed by Wang et al. (2013a, 2013b).

Finally, an X-ray study by Czesla et al. (2013) on the host star of WASP-43b using XMM-Newton observations found that the planet is exposed to a high-energy radiation field, and water abundance measurements were recently reported by Kreidberg et al. (2014).

In § 2, we describe the used telescopes and our observations of WASP-39b and WASP-43b, while in §§ 3 and 4, we focus on data reduction techniques and on a discussion about the period of the two objects. Results are presented in § 5, while conclusions are described in § 6.

2. Observations

Observations were carried on using all three equatorial Ritchey–Chrétien telescopes installed at the San Pedro Mártir National Astronomical Observatory (SPM-OAN). These facilities are located at an altitude of 2800 m in the middle of the peninsula of Baja California in northwestern Mexico (N 31°2'39'', W 11°527'49''). The 0.84 m and the 2.12 m are operated on-site, and the 1.50 m is a robotic telescope monitored and operated remotely.

We principally used the 0.84 m telescope, an f/15 which provides an 8.4' × 19.0' field of view on the Mexman instrument. We decided to window the Esopo CCD to avoid vignetting and to reduce the readout time. The instrument comes with a set of standard UBVRI Johnson filters, among others. In this article, we present observations carried out using the VRI filters.

Observations were also obtained with the robotic 1.50 m telescope and its multichannel imager RATIR (Reionization and Transients InfraRed) (Rapchun et al. 2011; Watson et al. 2012; Butler et al. 2012; Farah et al. 2012). This instrument, designed for gamma-ray burst observations and follow-up, is equipped with three dichroics to transmit the light to four detectors named C0, C1, C2, and C3, equipped with optical and infrared filters.

The CCDs work in 2 × 2 binning mode by default. In this article, we present observations carried out using the Gunn-i only, corresponding to the C1 detector which provides a 5.6' × 5.6' field of view.

Taking advantage from free observing time during maintenance operations, we had the possibility to use the 2.12 m telescope, improved with active optics on the primary mirror, which is normally dedicated to spectroscopy. We used this telescope in direct-imaging mode. The observation presented in this article was carried out using the U filter on the Marconi CCD providing a 6.2' × 6.2' field of view.

Additional details about telescopes and detectors are listed in Table 1, while Table 2 shows the log of the observations.

2.1. WASP-39b

We performed three observations of WASP-39b, covering two subsequent transits. Its parent star is a G-type dwarf with a reported magnitude of V = 12.1 (Faedi et al. 2011).

For the transit of 2014 March 17, we had simultaneous access to the 0.84 m and to the 2.12 m telescopes. Figure 1 shows the target and its nearby stars as seen in the fields of the two telescopes, some of which (A, B, and C) were tested as comparison stars, as seen by the fields of the two telescopes. We observed in the I band with the smaller telescope, obtaining 103 frames with exposure times of 120 s and 180 s, and in the U band with the larger one, in order to take advantage of its sensitivity, obtaining 57 frames with exposure times of 240 s and 300 s. As the object was particularly faint in this filter, we decided to use a 2 × 2 binning and a small defocus.

Fig. 1. 

Fig. 1.  Superposed fields of WASP-39b obtained with the 0.84 m and the 2.12 m telescopes, in two random CCD images, and using an inverse color scale. The target and field stars A, B, and C are also shown. Star A was chosen as a reference for photometry for both the telescopes.

The transit of 2014 March 21 was observed with the 0.84 m telescope only. We obtained 142 frames in the R band with an exposure time of 120 s, slightly defocusing the telescope to avoid saturation.

The night conditions during all observations suffered from very thin high haze and a nearly full Moon, with a seeing of about 3–4''. Depending on the instrument, the band, and seeing conditions, different defocusing levels were used. Moreover, rapid changes in the atmospheric conditions required a variation of the exposure time in order to avoid saturation.

2.2. WASP-43b

Due to its short period with respect to WASP-39b, we had the opportunity to carry out more observations of WASP-43b . We obtained a total of seven transits spanning six nights. Figure 2 shows two fields containing the target and nearby stars as seen with the 0.84 m and 1.50 m telescopes.

Fig. 2. 

Fig. 2.  Superposed fields of WASP-43b as obtained with the 0.84 m and the 1.50 m telescopes, in two random CCD images, and using an inverse color scale. The defocusing is particularly appreciable in the 1.50 m telescope images. The target field stars A, B, and C are also shown. Star A was chosen as a reference for photometry for both the telescopes.

Four transits were observed with the 0.84 m telescope only: on 2014 February 13, 2014 March 20, 2014 March 21, and 2014 May 12 in the R, I, V, and again I filters, respectively, with exposure times spanning from 60 s to 120 s.

Another transit was obtained with the 1.50 m telescope only (2014 March 29), strongly defocused, in the Gunn-i filter, and with a 90 s exposure time. This transit was unfortunately affected by clouds from the first contact until almost the midtime.

Finally, we observed with both the 0.84 m and 1.50 m telescopes on 2014 March 7, using the R and Gunn-i filters, respectively, and with exposure times of 60 s and 140 s.

All observations were taken in condition of a nearly full Moon.

3. Data Reduction

CCD images of both transits were debiased and flat-fielded using standard IRAF procedures from the ccdproc package. Master bias and master sky flat fields were obtained from images taken at the beginning and at the end of each observing night. Cosmic rays were also removed using the method implemented by the lakos package (van Dokkum 2001), then images were aligned using the imalign IRAF procedure.

We decided to use different aperture photometry routines to obtain the light curves: a custom IRAF routine and the defot routine (Southworth et al. 2010) written in the IDL language, that we modified to work with the FITS headers generated by the SPM-OAN telescopes. We found that the results of these two pipelines are in good agreement. In the case of images with a small amount of defocus, it was possible to calculate the full width half-maximum (FWHM) of the point-spread function by fitting an elliptical Gaussian. We used this information to dynamically adapt the aperture radius to a value of 2.5 times the calculated FWHM. We found that this technique increases the quality of the photometry with respect to the fixed-radius aperture. In the case of images with a consistent defocus, several fixed apertures were tested in order to find the optimal value.

Several field stars were tested in order to find a reference star to perform differential photometry, and a nonsaturated star with a color index and magnitude values as close as possible to those of the target star was chosen. The chosen reference stars for WASP-39b and WASP-43b are labeled as "A" in Figures 1 and 2, respectively.

We then used the target and the reference star to produce differential light curves.

The resulting light curves may present some residual trends, which were removed with a first-order airmass correction described by Ramón-Fox & Sada (2013).

The timestamps of the light curve were converted to the dynamical time-based system (BJDTDB, hereafter BJD) using the transformation given by Eastman et al. (2010).

We used the aperture photometry technique presented above for both WASP-39b, whose three final light curves are presented in Figure 3, and WASP-43b, whose seven final light curves are shown in Figure 4.

Fig. 3. 

Fig. 3.  Light curves of WASP-39b obtained with the 0.84 m telescope in R and I bands and with the 2.12 m telescope in the U band. For displaying purposes, the light curves were reported to the time of the last observation (2014 March 21), and an arbitrarily shift is applied. The fit of the light curves (solid lines) is also shown, as well as the residuals in the lower box.

Fig. 4. 

Fig. 4.  Light curves of WASP-43b obtained with the 0.84 m telescope in VRI bands and with the 1.50 m telescope in the Gunn-i band, reported to the time of the observation of 2014 March 29, and arbitrarily shifted. The fit and the residuals are also shown. The dashed line represent an interpolation of the model.

3.1. Fit

The light curves were fitted using the IDL software Transit Analysis Package (TAP) implemented by Gazak et al. (2011), which uses a Markov-chain Monte Carlo (MCMC) method to find the best-fit parameters for the Mandel & Agol (2002) model.

This code allows the fitting of multiple light curves simultaneously, which is particularly useful for fixing parameters such as the orbital inclination i and the scaled semimajor axis a/R to the same values for a group of curves, and thus obtain a global fit for these quantities.

Other parameters, such as the midtransit time Tmid and the scaled planetary radius Rp /R, are allowed to be fitted individually for each curve. This technique allows fitting of the same orbital model to a group of light curves of the same object obtained in different bands, and the finding of potential planetary radius differences as a function of the wavelength.

All light curves of each system were fitted simultaneously. We used a set of fixed values in the MCMC analysis for several parameters: the period P, taken from Faedi et al. (2011) for WASP-39b and from Chen et al. (2014) for WASP-43b; the eccentricity e and the argument of periastron ω, both set = 0. TAP allows the use of linear and quadratic models for the stellar limb darkening. A quadratic model is assumed in this analysis, and the two terms l1 (linear) and l2 (quadratic) are fixed in the MCMC analysis to the values corresponding to the filter used for each light curve. The values were obtained from the exofast (Eastman et al. 2012, 2013) online tool, which are interpolated from stellar atmosphere models of Claret & Bloemen (2011). 2

For each system, we then fitted the orbital inclination i, the scaled semimajor axis a/R, and a value of the scaled radius Rp /R for each one of the used filters. TAP was initialized with the most recent parameters reported in Faedi et al. (2011) for WASP-39b and Chen et al. (2014) for WASP-43b .

The midtransit time Tmid for each light curve was also fitted (in the case of multiple light curves of the same object obtained on the same night using different telescopes, the same Tmid was fitted for both the light curves). For the MCMC analysis, the Tmid is initialized with a value estimated by TAP from the input light curve.

The best-fit models and corresponding residuals are shown in Figures 3 and 4, and the values of the parameters are shown in Table 3.

4. Period And Period Variation

Using the fitted values for the midtransit time Tmid, we retrieved the period P of WASP-39b and WASP-43b by fitting the following linear law:

where T0 is the initial epoch and N is the number of periods since T0. The terms T0 and P are best-fit parameters.

Further analysis of the ephemerids was carried out in order to investigate a long-term variation of the period with respect to the time. A constant decrease of the period for close-orbiting planets may be indicative of processes that remove orbital energy such as tidal dissipation (Adams et al. 2010a, 2010b; Sasselov 2003; Levrard et al. 2009). A simple model assuming a constant variation of the period is given by the following quadratic equation

proposed by Adams et al. (2010a), where .

For the present analysis, a Levenberg–Marquardt least-squares fitting algorithm implemented in the code of Markwardt (2009) was initially used to find the best-fit parameters of equations (1) and (2). However, this algorithm can be trapped in a local minimum of χ2. For this reason, a Monte Carlo code was set-up to search for the χ2 minimum and sample the parameter space. Both approaches were tested on the OGLE-TR-113b data provided by Adams et al. (2010b), finding good agreement with their results.

We find that both approaches are consistent and give similar best-fit values. Comforted by these results, we decided to apply this method to our data, and we describe the results in the following subsections.

4.1. WASP-39b

Concerning WASP-39b, the Tmid reported by Faedi et al. (2011) and the one calculated from our observations are fitted using Eq. 1. This gives the following values:

  • 1.  
    T0 = 2455342.9687998 ± 0.0001999 (BJD),
  • 2.  
    P = 4.0552989 ± 9.963 × 10-7 days
  • 3.  
    χ2 = 0.541

where χ2 stands for the "reduced χ2," i.e., χ2/Df , where Df are the degrees of freedom (number of data–number of parameters).

However, due to the significantly small number of observations, we decided to increase the number of transits considered in our analysis by introducing additional data from the Exoplanet Transit Database (ETD) Web page (Poddaný et al. 2010), in order to improve the robustness of our results. We selected a total of three light curves. According to this database, the timings of the light curves are reported in HJD. They were transformed to BJD using the software provided by Eastman et al. (2010). We used TAP to obtain only Tmid for these light curves by running the MCMC analysis with all other parameters to be fixed to the values provided by Faedi et al. (2011). By adding these new data to previous ones, we obtain a total of seven observations. The linear fit of equation (1) was repeated with these additional transits, and the following results were obtained:

  • 1.  
    T0 = 2455342.9695727 ± 0.000199 (BJD)
  • 2.  
    P = 4.0552947 ± 9.651 × 10-7 days
  • 3.  
    χ2 = 11.82

These indicate a period larger than that reported by Faedi et al. (2011). The residuals for this fit are shown in Figure 5. As the last two Tresca points show large differences with respect to the fit, we tested the fit process without their contribution, obtaining the following results:

  • 1.  
    T0 = 2455342.9695554 ± 0.0001994 (BJD),
  • 2.  
    P = 4.0552965 ± 9.956 × 10-7 days
  • 3.  
    χ2 = 1.299

In this case, too, we report a larger period with respect to Faedi et al. (2011) results. In both cases, the results are consistent, but we suggest to improve their robustness with additional observations.

Fig. 5. 

Fig. 5.  Observed–calculated Tmid and linear fit of the WASP-39b period, including data from Faedi et al. (2011), from the ETD-TRESCA database, and from our observations.

An attempt to fit equation (2) was made, but the small amount of data do not allow a reasonable result for the fit. More data are needed as well as precise timings, in order to assess if there is a long-term constant variation of the period, or to detect the presence of other transit timing variations.

4.2. WASP-43b

A sufficient amount of WASP-43b transit timings is available in the literature spanning several years of observations (Hellier et al. 2011; Gillon et al. 2012; Chen et al. 2014; Murgas et al. 2014).

Some previous studies of possible long-term variations have been carried out by Chen et al. (2014) and Murgas et al. (2014), suggesting that the period of WASP-43b is slowly decreasing with time. In this work, we repeat the analysis combining our data with the timings available in the literature. A best-fit model was obtained using the Marquardt–Levenberg algorithm previously described. For the linear fit, best-fit values are:

  • 1.  
    T0 = 2455528.86839966 ± 0.000036904078 (BJD)
  • 2.  
    P = 0.81347409 ± 7.3 × 10-8 days
  • 3.  
    χ2 = 3.784

For the quadratic fit, best-fit values are

  • 1.  
    T0 = 2455528.86848347 ± 0.0000369043 (BJD)
  • 2.  
    P = 0.81347374 ± 7.285 × 10-8 days
  • 3.  
    δP = 6.7 ± 6.9 × 10-10
  • 4.  
  • 5.  
    χ2 = 3.819

The small difference in χ2 between the linear and quadratic fits shows that, statistically, there is no significant improvement by using the latter, as shown in Figure 6.

Fig. 6. 

Fig. 6.  Observed–calculated Tmid, linear and quadratic fit of the WASP-43b period including data from Hellier et al. (2011); Gillon et al. (2012); Chen et al. (2014), and from our observations. See § 4.2 for curve details.

Using the values provided by Blecic et al. (2014), we obtain the following:

  • 1.  
    P = 0.81347530 ± 3.9 × 10-7 days
  • 2.  
  • 3.  
    χ2 = 21.25

The parameters reported by Chen et al. (2014) are as follows:

  • 1.  
    P = 0.81347399 ± 2.2 × 10-7 days
  • 2.  
  • 3.  
    χ2 = 18.26

The fit was also tested with the Monte Carlo approach, which found a minimum χ2 corresponding to P = 0.81347382 and , with a χ2 = 4.046. This analysis showed that the χ2 minimum is very close to a value of .

These fits are visualized in Figure 6. The model of Blecic et al. (2014) fits most of the points reasonably well, except for those corresponding to our observations. The model of Chen et al. (2014) shows more deviations with respect to recent observations. For comparison, we also plot a third model using the same period reported by Blecic et al. (2014) but with δP = -1.0 × 10-9, which is in the order of magnitude of previously reported values, giving a χ2 = 46.96. This corresponds to the same orbital period of Blecic et al. (2014) but with a smaller value of . It can be seen that it fits reasonably points from earlier observations. However, a detailed analysis of transit time variations has not been considered in this work, which could explain some of the observed variations (Chen et al. 2014). The present analysis then confirms the results of Chen et al. (2014), also showing that a quadratic model does not improve significantly the fit of the transit timings for WASP-43b.

5. Results And Discussion

We report for the first time a light curve of WASP-39b in the U band obtained with the 2.12 m telescope at SPM-OAN. Observations in this band are important to study the presence of high-hazes in the atmosphere, as Rayleigh scattering becomes important at smaller wavelengths, which can produce larger observed radii in this band (Southworth et al. 2012a; Copperwheat et al. 2013). The resulting light curve is shown in Figure 3. Although there is more dispersion in the data compared to the I and R bands, a depth flux of approximately 3% is observed. The best-fit values obtained in the previous section give a relative radius of Rp /R = 0.1462 ± 0.0116.

No significant variations of the planetary radius between the U, R, and I bands is found. It is also noted that the U and I band light curves were obtained in the same night, and no evident asymmetries are observed that could suggest the presence of tails or other features around the planet. Additional observations in the U band can help to reduce the uncertainty for this highly inflated extrasolar planet.

The analysis of the ephemerids in the previous section gives a period for WASP-39b approximately 3.084 ± 0.774 s larger than that of Faedi et al. (2011).

WASP-43b shows a slight dependence of the Rp /R as a function of the used filter (see Fig. 7), which should be investigated with additional observations. In particular, we find for filter i a value higher by 0.01637 ± 0.00371 with respect to the value reported by Chen et al. (2014) in the same filter (0.1738 ± 0.0033 against 0.15847 ± 0.00041). Although, as our two i filter curves are not complete and show signatures of red noise, we suggest to carry on more observations to confirm the significancy of this result.

Fig. 7. 

Fig. 7.  WASP-43b planet radius in the observing filters (see Table 3). The value reported by Chen et al. (2014) in the i filter is also shown.

Nevertheless, we find a tendency for the Rp /R similar to that of Chen et al. (2014) in the sense that the planetary radius is significantly higher at the i band with respect to the value in other bands. The semimajor axis is slightly smaller, but comparable to the value obtained by Gillon et al. (2012) and Chen et al. (2014).

The analysis of the ephemerids gives a period consistent with previous results, and we confirm that there is no improvement while using a quadratic model for fitting the ephemerids. Further analysis and observations would be required to ascertain if a constant period variation exists for WASP-43b . Given the proximity of this planet to its host star, this is particularly relevant to explore potential mechanisms that remove orbital energy such as tidal dissipation.

The present work is part of an ongoing survey started in 2014. A total of 23 extrasolar planets were observed mainly with the 0.84 m and the 1.50 m telescopes, giving a total of 40 good-quality light curves in several filters. Upgrades in the instrumentation, refinement of the observing technique, and the possibility to benefit from dark and photometric nights will increase the quality of the results in the near future, and will complement other telescope data for multisite studies and analysis of this class of objects. Moreover, three new 1.30 m telescopes will be completed by 2016 at the San Pedro Mártir observatory, which are part of the TAOS II project (Lehner et al. 2012a; Geary et al. 2012; Lehner et al. 2013).

The characteristics of this survey (Lehner et al. 2012b) allow specific pipelines to be set up for the early detection of new extrasolar planets, for example with a binning of the light curves as already tested for TAOS data (Ricci et al. 2014). Also, a project for the installation of a 6.5 m telescope at the San Pedro Mártir observatory has recently reached the phase of preliminary design.

Current and in fieri SPM-OAN telescopes then represent a unique opportunity for the extrasolar planet community to complement the northern hemisphere observatories involved in this research area.

6. Conclusions

Multifilter exoplanet transit observations of WASP-39b and WASP-43b, carried on for the first time with all three San Pedro Mártir telescopes, have shown the capability of these facilities of 0.84 m, 1.50 m, and 2.12 m for this kind of investigation. The fit of the two objects shows a scatter of 1.5–2.5 mmag rms in condition of full or nearly full Moon. This makes these instruments suitable for future ground-based observing campaigns, in particular for the follow-up of alerts triggered by the upcoming projects currently in development.

The analysis shows most of the parameters to be in good agreement with previous works, in particular for what concerns the first observation WASP-39b in the U filter. The period for WASP-39b is found to be significantly larger (>3 s) with a 4σ accuracy with respect to previous works. The value of the planet/star radius of WASP-43b is larger in the i filter with respect to previous works. A tendency of a higher planet radius in the i band is obtained, consistent with the tendency reported in recent works. Additional observations of this object are needed to confirm a slight dependence of this parameter on the observing filter. We plan additional observations in order to confirm the accuracy of these results with respect to the literature.

Research carried out thanks to the support of UNAM-DGAPA-PAPIIT project IN115413. FGRF acknowledges the support of a CONACYT scholarship for postgraduate studies in Mexico. We also acknowledge Dr. Wolfgang Steffen Burg for the discussion about future exoplanetary-related projects, John Southworth for his useful comments about the form and the content of this manuscript, and the OAN technical and service staff for daily support. We finally acknowledge the anonymous referee for the useful remarks.

Facilities: San Pedro Mártir 0.84 m (MEXMAN), San Pedro Mártir 1.50 m (RATIR), San Pedro Mártir 2.12 m (direct imaging mode).

Footnotes

  • http://exoplanetarchive.ipac.caltech.edu/cgi-bin/ExoTables/nph-exotbls?dataset=planets

  • http://astroutils.astronomy.ohio-state.edu/exofast/limbdark.shtml

Please wait… references are loading.