HATS-9b and HATS-10b

HATS-9b and HATS-10b: Two compact hot Jupiters in field 7 of the K2 mission1


We report the discovery of two transiting extrasolar planets by the HATSouth survey. HATS-9b orbits an old ( Gyr) V=13.3 G dwarf star, with a period  d. The host star has a mass of 1.03 , radius of 1.503  and effective temperature  K. The planetary companion has a mass of 0.837 , and radius of 1.065  yielding a mean density of 0.85 . HATS-10b orbits a V=13.1 G dwarf star, with a period  d. The host star has a mass of 1.1 , radius of 1.11   and effective temperature  K. The planetary companion has a mass of 0.53 , and radius of 0.97  yielding a mean density of 0.7 . Both planets are compact in comparison with planets receiving similar irradiation from their host stars, and lie in the nominal coordinates of Field 7 of K2 but only HATS-9b falls on working silicon. Future characterization of HATS-9b with the exquisite photometric precision of the Kepler telescope may provide measurements of its reflected light signature.

Subject headings:
planetary systems — stars: individual ( HATS-9, GSC 6305-02502, HATS-10, GSC 6311-00085 ) techniques: spectroscopic, photometric

1. Introduction

Our current understanding of the structure and orbital evolution of extrasolar giant planets has been, to a large degree, informed by the characterization of transiting planetary systems. Besides the determination of the planet radius, true mass, and bulk density, follow-up studies of transiting extrasolar planets (TEPs) allow the extraction of valuable information, like the spin-orbit angle and the properties and composition of the planetary atmospheres, that are not be easily recovered unless the orbital plane is favorably oriented such that the planet eclipses its host star.

Detections of giant TEPs, mostly driven by transiting ground based surveys like SuperWASP (Pollacco et al., 2006) and HATNet (Bakos et al., 2004), have revealed a large number of systems in the region of parameter space with 0.8, 0.4, P5 d and FGK-type host stars. The measured properties of these systems, coupled with subsequent follow-up studies, have been fundamental for testing formation and interior models of these giant planets, which are known as hot Jupiters.

New ground-based transiting surveys like HATSouth (Bakos et al., 2013) have been designed with the goal of expanding the parameter space of well characterized TEPs by detecting planets with smaller radii () and/or longer periods ( days). In the process of searching for these kinds of planets, new hot Jupiters are detected which contribute to enlarging the sample of known systems. Even though many hot Jupiters are already known, more are still needed to make headway into understanding their physical properties, e.g. a firm understanding of the mechanism that causes some hot Jupiters to have inflated radii (e.g. HAT-P-32b and HAT-P-33b, Hartman et al., 2011).

New planet discoveries around bright stars, accessible by follow-up facilities are especially valuable given the wealth of detailed studies that they can be subject to. Indeed, some of the most analyzed and characterized giant TEPs are three planets (TrES-2b, HAT-P-7b and HAT-P-11b) that were detected by ground-based surveys (Pál et al., 2008; Bakos et al., 2010; O’Donovan et al., 2006) and later observed by NASA’s Kepler mission (Borucki et al., 2010). Even though the primary goal of Kepler satellite was the detection of planets near the habitable zone for estimating their frequency and distribution in our galaxy, the high photometric precision of Kepler allowed very detailed studies of the small population of giant planets on close orbits around moderately bright stars () that fell in its field of view.

Kepler was able to detect secondary transits and phase variations on TrES-2b and HAT-P-7b (Esteves et al., 2013) which were useful in the study of their atmospherical properties, such as the determination of the geometric albedos and planetary phase curve offsets. Doppler beaming and ellipsoidal variations measured with Kepler also constrained the mass of those planets. In the case of HAT-P-11b, Kepler observations were useful in characterizing the activity of the K-type host star; and the analysis of crossing stellar spots allowed the determination of the spin-orbit misalignment of this system (Deming et al., 2011; Sanchis-Ojeda & Winn, 2011). Simultaneous observations of the transits of HAT-P-11b by Kepler and Spitzer allowed also the detection of water vapor in the atmosphere of this Neptune-size planet (Fraine et al., 2014). Estimation of the planetary physical parameters depend strongly on the estimated stellar properties. In this regard, Kepler was also able to measure model independent stellar properties by the use of asteroseismology on the three mentioned systems (Christensen-Dalsgaard et al., 2010).

After the failure of two of its reaction wheels, the Kepler satellite is still working, but with a new observation strategy and a photometric precision within a factor of 2 of the nominal Kepler mission performance (e.g., Vanderburg & Johnson, 2014; Aigrain et al., 2015; Foreman-Mackey et al., 2015; Crossfield et al., 2015). This new mission concept, called K2 (Howell et al., 2014), will observe 10 fields, each for a period of approximately 70 days, and some of these fields lie in the Southern hemisphere. One of the limitations of K2 is that the number of stars that can be monitored in each field is substantially lower than for the original Kepler mission. For this reason, the pre-selection of targets, based on ground-based observations of K2 fields is especially important for an efficient use of the satellite.

In this work we present the discovery of HATS-9b and HATS-10b, two hot Jupiters discovered by the HATSouth survey which are located in the nominal coordinates of Field 7 of K2 mission. In Section 2 we summarize the observations that allowed the discovery and confirmation of these planets. In Section 3 we show the global analysis of the spectroscopic and photometric data that confirmed the planetary nature of the transiting candidates and also rejected blend scenarios that can mimic the photometric and radial velocity signals. Our findings are discussed in Section 4.

2. Observations

2.1. Photometric detection

HATS-9 and HATS-10 were identified as transiting planetary host candidates after obtaining 10000 images of the same field with three stations on the three HATSouth observing sites. The number of photometric observations that were taken for each star on each of the HATSouth stations is indicated in Table 1, where it can be seen that in both cases 45 of the observations came from the HATSouth station located at Las Campanas Observatory (LCO).

The HATSouth observations consist of four-minute Sloan -band exposures obtained with 24 Takahashi E180 astrographs (18cm aperture) coupled to Apogee 4Kx4K U16M ALTA CCDs. Readout times are of the order of one minute which results in a cadence of about 5 minutes. Detailed descriptions of the image processing steps and the candidate identification procedures of HATSouth data can be found in Bakos et al. (2013) and Penev et al. (2013). Briefly, after applying aperture photometry on the images, the light curves generated are detrended using external parameter decorrelation (EPD) and trend filtering algorithm (TFA Kovács et al., 2005). Periodic transits on the detrended light curves are then searched using Box-fitted Least Squares (BLS) algorithm (Kovács et al., 2002).

Figure 1 shows the phase-folded detection light curves of HATS-9b and HATS-10b, where a clear 10 mmag flat-bottom transit can be observed in both cases.

Figure 1.— Phase-folded unbinned HATSouth light curves for HATS-9 (left) and HATS-10 (right). In each case we show two panels. The top panel shows the full light curve, while the bottom panel shows the light curve zoomed-in on the transit. The solid lines show the model fits to the light curves. The dark filled circles in the bottom panels show the light curves binned in phase with a bin size of 0.002.
Instrument/Field33 Date(s) # Images Cadence34 Filter Precision35
(sec) (mmag)
    HS-1/G579 2010 Mar–2011 Aug 4317 300  band 6.9
    HS-3/G579 2010 Mar–2011 Aug 2138 303  band 7.6
    HS-5/G579 2010 Sep–2011 Aug 2784 303  band 6.9
    FTS 2013 Apr 11 134 80  band 1.4
    PEST 2013 May 31 186 130  band 3.4
    HS-1/G579 2009 Sep–2011 Aug 4389 301  band 7.3
    HS-3/G579 2010 Mar–2011 Aug 2596 303  band 7.2
    HS-5/G579 2011 Mar–2011 Aug 3297 303  band 7.8
    CTIO 0.9m 2012 Aug 29 69 213  band 2.3
    FTS 2013 Apr 05 142 63  band 4.3
    GROND 2013 Jun 14 92 156  band 0.8
    GROND 2013 Jun 14 88 156  band 1.3
    GROND 2013 Jun 14 94 156  band 0.7
    GROND 2013 Jun 14 89 156  band 0.8
    PEST 2013 Jun 27 145 130  band 4.6
Table 1 Summary of photometric observations

2.2. Spectroscopic Observations

Transit-like light curves can be produced by different configurations of stellar binaries. Spectroscopic observations are required to reject false positives and to obtain the orbital parameters and masses of the true planets. Due to the great number of HATSouth candidates and the limited available observing time on spectroscopic facilities, this follow-up is performed in a two-step procedure as we now describe. All spectroscopic observations are summarized in Table 2.

First, initial spectra are acquired (with either low resolution, or low S/N) to make a rough estimation of the stellar parameters, identifying spectra composed of more than one star, and measuring RV variations produced by stellar mass companions. HATS-9 was observed with WIFeS (Dopita et al., 2007) on the ANU 2.3m telescope, obtaining =5821 300 K, =3.9 0.3 =0.5 0.5; and with ARCES on the APO 3.5m obtaining =5692 50 K, =4.14 0.1 =0.5 0.08. Both estimates of stellar parameters were consistent with a G-type dwarf, but the sub-solar surface gravity value points towards a slightly evolved system. Details on the observing strategy, reduction methods and the processing of the spectra for WIFeS can be found in Bayliss et al. (2013). The ARCES observation was carried out using the slit yielding an echelle spectrum with 107 orders covering the wavelength range 3200–10000Å at a resolution of . A single ThAr lamp spectrum was obtained immediately following the science exposure with the telescope still pointed toward HATS-9. The science observation was reduced to a wavelength-calibrated spectrum using the standard IRAF echelle package36, and analyzed using the Spectral Parameter Classification (SPC) program (Buchhave et al., 2012) to determine the radial velocity and stellar atmospheric parameters.

Reconnaissance spectroscopy was performed for HATS-10 using the echelle spectrograph mounted on the du Pont 2.5m telescope at Las Campanas Observatory. One observation using the slit () was enough to confirm that HATS-10 has a single lined spectrum with the following stellar parameters: =6100 100 K, =4.6 0.5 =0.0 0.5, =5.0 2.0 km/s. This spectrum was reduced and analysed with an automated pipeline developed to deal with data coming from a host of different echelle spectrographs (Brahm et al. in prep.). The pipeline for DuPont is very similar to the ones we have previously detailed for Coralie and FEROS data in Jordán et al. (2014).

Once both candidates were identified as single-lined late-type dwarfs, spectra from high precision instruments were required to measure RV variations with high precision ( m/s) in order to measure the mass of the substellar companions and obtain the orbital parameters. HATS-9 and HATS-10 were observed several times with Coralie (Queloz et al., 2001) on the 1.2m Euler telescope, FEROS (Kaufer & Pasquini, 1998) on the 2.2m MPG telescope, and HDS on the 8m Subaru telescope (Noguchi et al., 2002). Coralie and FEROS data were processed with the pipeline described in Jordán et al. (2014), where RV values are obtained using the cross correlation technique against a binary mask and bisector span (BS) measurements are computed from the cross-correlation peak following Queloz et al. (2001). HDS RVs were measured using the procedure detailed in Sato et al. (2002, 2012) which are in turn based on the method of Butler et al. (1996) while BS values were obtained following Bakos et al. (2007).

Phased high-precision RV and BS measurements are shown for each system in Figure 2 and the data are listed in Table 3. Both candidates show RV variations in phase with photometric ephemeris, however for HATS-10 the residuals are higher than expected. This deviation can be partly explained by moonlight contamination in 5 spectra acquired with Coralie in August 2013 which are marked with crosses in Figure 2. There are no significant correlations between RV and BS variations and thus we conclude the RV variations are not produced by stellar activity. The 95% confidence interval for the Pearson correlation coefficient between RV and BS was computed for both candidates using a bootstrap procedure. The confidence intervals are [-0.57,0.07] and [-0.43, 0.37] for HATS-9 and HATS-10, respectively. The individual FEROS spectra were median combined for both candidates to perform a precise estimation of the stellar parameters.

Instrument UT Date(s) # Spec. Res. S/N Range37 38 RV Precision39
//1000 () ()
APO 3.5 m/ARCES 2012 Aug 25 1 31.5 27 -11.5 500
ANU 2.3 m/WiFeS 2012 Sep 8 1 3 140
Euler 1.2 m/Coralie 2012 Nov 6–10 4 60 14–20 -10.634 37
MPG 2.2 m/FEROS 2012 Aug–2013 May 9 48 32–76 -10.653 32
Subaru 8 m/HDS 2012 Sep 19 3 60 100–114
Subaru 8 m/HDS+I 2012 Sep 20–22 9 60 60–100 11
du Pont 2.5 m/Echelle 2013 Aug 21 1 40 48 -29.2 500
Euler 1.2 m/Coralie 2012 Aug–2013 Aug 12 60 17–23 -28.131 68
MPG 2.2 m/FEROS 2013 Mar–Jul 5 48 29–85 -28.044 50
Subaru 8 m/HDS 2012 Sep 22 3 60 74–94
Subaru 8 m/HDS+I 2012 Sep 19–21 9 60 41–99 14
Table 2 Summary of spectroscopy observations

Figure 2.— Phased high-precision RV measurements for HATS-9 (left), and HATS-10 (right) from HDS (filled circles), FEROS (open triangles), and Coralie (filled triangles). In each case we show three panels. The top panel shows the phased measurements together with our best-fit model (see Table 6) for each system. Zero-phase corresponds to the time of mid-transit. The center-of-mass velocity has been subtracted. The second panel shows the velocity residuals from the best fit. The error bars include the jitter terms listed in Table 6 added in quadrature to the formal errors for each instrument. The third panel shows the bisector spans (BS), with the mean value subtracted. Note the different vertical scales of the panels. RV measurements highly contaminated with moonlight are marked with crosses.
BJD RV40 41 BS Phase Instrument
(2,456,000) () () () ()

Note. – Note that for the iodine-free template exposures we do not measure the RV but do measure the BS and S index. Such template exposures can be distinguished by the missing RV value. The Subaru/HDS observations of HATS-10 without BS measurements have too low S/N in the I-free blue spectral region to pass our quality threshold for calculating accurate BS values.

Table 3 Relative radial velocities and bisector spans for HATS-9 and HATS-10.

2.3. Photometric follow-up observations

Figure 3.— Left: Unbinned transit light curves for HATS-9. The light curves have been corrected for quadratic trends in time fitted simultaneously with the transit model. The dates of the events, filters and instruments used are indicated. The second curve is displaced vertically for clarity. Our best fit from the global modeling described in Section 3.3 is shown by the solid lines. Right: residuals from the fits are displayed in the same order as the left curves. The error bars represent the photon and background shot noise, plus the readout noise.
Figure 4.— Similar to Figure 3; here we show the follow-up light curves for HATS-10.

In order to confirm the occurrence of the transits and to better constrain the orbital and physical parameters of the companions, higher precision light curves for both candidates were acquired using several telescopes around the globe. Table 1 summarizes the key aspects of this photometric follow-up, including the dates of the observations, the cadence and the filter.

Two partial transits of HATS-9 were detected using the 0.3m Perth Exoplanet Survey Telescope (PEST) and the Spectral camera on the 2m Faulkes Telescope South (FTS), part of Las Cumbres Observatory Global Telescope (LCOGT). Results of these observations are presented in Table 4 and shown in Figure 3. Two partial transits of HATS-10 were observed with FTS and the CTIO 0.9m telescope. Another two full transits were measured with PEST and the GROND instrument on the MPG 2.2 m. These HATS-10 light curves are shown in Figure 4. All the facilities used for high precision photometric follow-up have been previously used by HATSouth; the instrument specifications, observation strategies and reduction procedures adopted can be found in Bayliss et al. (2013), Zhou et al. (2014b), Hartman et al. (2014) and Mohler-Fischer et al. (2013) for FTS, PEST, CTIO 0.9m, and GROND, respectively.

Object43 BJD44 Mag45 Mag(orig)46 Filter Instrument

Note. – This table is available in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.

Table 4 Light curve data for HATS-9 and HATS-10.

3. Analysis

3.1. Properties of the parent star

We determine precise stellar parameters for HATS-9 and HATS-10 using a new code called ZASPE (Zonal Atmospherical Stellar Parameter Estimator) on median combined FEROS spectra. The detailed structure and performance of ZASPE will be presented elsewhere (Brahm et al. in preparation), but in summary ZASPE is a Python-based code that computes the between the observed spectra and the PHOENIX grid of synthetic spectra (Husser et al., 2013) only on the most sensitive spectral zones to each stellar parameter. The optimal set of stellar parameters (,, and ) is found iteratively and the sensitive zones are determined in each iteration. One of the most novel features of ZASPE is that the errors on the stellar parameters are computed from the data itself and include the systematic mismatches between the observations and the best fitted model. We have validated the results of ZASPE against a set of stars with interferometrically determined stellar parameters (Boyajian et al., 2012) and which have publicly available FEROS spectra. Results of this comparison are shown in Figure 5. The resulting parameters for HATS-9 are: =5363 90 K, =3.97 0.2, =0.33 0.09, =4.67 0.5 km/s; while for HATS-10 we get: =5974 110 K, =4.44 0.13, =0.19 0.07, =5.66 0.5 km/s.

Figure 5.— (Left): Comparison between T values obtained with ZASPE and the ones derived from interferometric information. (Right): Comparison between log(g) values obtained with ZASPE and the ones derived from interferometric information.

These sets of stellar parameters were refined using the information contained in the transit light-curves. The stellar mean density () can be computed directly from one of the light-curve model parameters () and the period and eccentricity of the orbit by using Kepler’s third law with only a slight dependence on the stellar parameters through the limb-darkening coefficients (Sozzetti et al., 2007). The spectroscopically determined and were coupled with and Yonsei-Yale stellar evolution models (Y2; Yi et al., 2001) to determine the stellar physical parameters (, and the age of the star), which were used to compute a new and more precise estimation of for HATS-9 (=4.12 0.04) and HATS-10 (=4.38 0.03). A new set of , , was determined using ZASPE with fixed to the precise values obtained by modeling the light curves, followed by a new estimation of and a new modeling of stellar isochrones. The new set of stellar parameters fixing , which are the ones we adopted for further analysis, were consistent with the initial values quoted in the previous paragraph and are listed in Table 5, where distances are determined by comparing the measured broad-band photometry listed in that table to the predicted magnitudes in each filter from the isochrones. We assume a extinction law from Cardelli et al. (1989) to determine the extinction. The 1 and 2 confidence ellipsoids in and are plotted in Figure 6 for both planet hosts, along with the Y2 isochrones for the ZASPE determined . We find that HATS-9 is a , , quite evolved ( Gyr) star, while HATS-10 is a , main-sequence star.

We attempted to measure the Lithium absorption line at 6707.8 Å for testing the age estimation of HATS-9 but the quality of our spectra was only enough to rule out a strong absorption feature.

Figure 6.— Model isochrones from Yi et al. (2001) for the measured metallicities of HATS-9 (left) and HATS-10 (right). In each case we show models for ages of 0.2 Gyr and 1.0 to 14.0 Gyr in 1.0 Gyr increments (ages increasing from left to right). The adopted values of and  are shown together with their 1 and 2 confidence ellipsoids. The initial values of  and  from the first ZASPE and light curve analyzes are represented with a triangle.
        Parameter Value Value Source

Astrometric properties and cross-identifications
    2MASS-ID 2MASS 19231442-2009587 2MASS 19371363-2212161
    GSC-ID GSC 6305-02502 GSC 6311-00085
    R.A. (J2000) 2MASS
    Dec. (J2000) 2MASS
     () UCAC4
     () UCAC4
Spectroscopic properties
     (K) ZASPE47
     () ZASPE
     () 4.6 3.8 Assumed48
     () 1.0 1.0 Assumed49
     () Coralie,FEROS
Photometric properties
     (mag) APASS50
     (mag) APASS51
     (mag) APASS52
     (mag) APASS53
     (mag) APASS54
     (mag) 2MASS
     (mag) 2MASS
     (mag) 2MASS
Derived properties
     () YY++ZASPE 55
     () YY++ZASPE
     (cgs) YY++ZASPE
     () YY++ZASPE 56
     () YY++ZASPE
     (mag) YY++ZASPE
     (mag,ESO) YY++ZASPE
    Age (Gyr) YY++ZASPE
     (mag) YY++ZASPE
    Distance (pc) YY++ZASPE
Table 5 Stellar parameters for HATS-9 and HATS-10

3.2. Excluding blend scenarios

In order to exclude blend scenarios we carried out a blend analysis of the observations following Hartman et al. (2012). For HATS-9 we find that scenarios involving blends between a stellar eclipsing binary and a foreground or background star can be ruled out with greater than 5 confidence based on the photometric data alone. The primary constraint in this case is the lack of out-of-transit variations seen in the HATSouth light curve. Due to the short orbital period, the best-fit blend model which reproduces the shape of the transit has a  mmag amplitude ellipsoidal variation, and a  mmag deep secondary eclipse, neither of which are detected in the HATSouth observations. Moreover the Subaru/HDS observations of HATS-9 show no significant bisector span variation (the RMS scatter of the BS measurements is 12 ) providing further evidence that the system is not a blended eclipsing binary. For HATS-10 the photometric observations can be fit by a G+M star eclipsing binary blended with another G star that is slightly brighter than the primary in the eclipsing system. Based on the difference in , this model is indistinguishable from a single G star with a transiting planet. We simulated spectra for blend models that could plausibly fit the photometric observations, finding that it all cases the blended systems would have easily been detected as having composite spectra. They also would produce RV and BS variations of several , whereas the observed RV variation is , and the Subaru/HDS BS scatter is only 18 . We conclude that neither HATS-9 nor HATS-10 is a blended eclipsing binary system. As is often the case, however, we are not able to rule out the possibility that either transiting planet system has a fainter stellar-mass companion. For both systems a stellar companion of any mass, up to the mass of the planet-hosting star, is possible. If a massive stellar companion is present in a given system, the true planet radius would be up to % larger than inferred here. The planet mass would also be larger. High resolution adaptive optics imaging, and/or long-term RV observations are needed to determine whether either system has a stellar companion (e.g. Howell et al., 2011; Horch et al., 2014; Everett et al., 2015).

3.3. Global modeling of the data

We modeled the HATSouth photometry, the follow-up photometry, and the high-precision RV measurements following Pál et al. (2008); Bakos et al. (2010); Hartman et al. (2012). We fit Mandel & Agol (2002) transit models to the light curves, allowing for a dilution of the HATSouth transit depth as a result of blending from neighboring stars and over-correction by the trend-filtering method. For the follow-up light curves we include a quadratic trend in time in our model for each event to correct for systematic errors in the photometry. We fit Keplerian orbits to the RV curves allowing the zero-point for each instrument to vary independently in the fit, and allowing for RV jitter which we we also vary as a free parameter for each instrument.

We used a Differential Evolution Markov Chain Monte Carlo procedure (ter Braak, 2006; Eastman et al., 2013) to explore the fitness landscape and to determine the posterior distribution of the parameters.

The resulting parameters for each system are listed in Table 6. HATS-9b has a radius of and a mass of , while HATS-10b has a radius of and a mass of . Both planets have bulk densities slightly lower than the one of Jupiter ( g cm and g cm, respectively)

HATS-9b HATS-10b
               Parameter Value Value

Light curve parameters
    () 57
    (days) 58
    (days) 59
Limb-darkening coefficients 61
    (linear term)
    (quadratic term)
RV parameters
   RV jitter HDS () 63
   RV jitter FEROS ()
   RV jitter Coralie ()
Planetary parameters
    (cgs) 66
Table 6Orbital and planetary parameters for HATS-9b and HATS-10b

4. Discussion

We have presented the discovery of two new transiting planets which are shown on mass-radius and equilibrium temperature versus radius diagrams in Figure 7. From the mass-radius diagram, HATS-9b and HATS-10b can be classified as typical non-inflated hot Jupiters. HATS-9b is slightly less massive than Jupiter (0.84 ) and has almost the same radius. Its orbital period of days is rather short compared to the period distribution of known hot Jupiters. HATS-10b has a mass in the range between Saturn and Jupiter (0.53 ), a radius consistent with that of Jupiter and a period of days, which is close to the mean period of known hot Jupiters.

The equilibrium temperature versus radius diagram shows that both planets tend to depart from the known correlation between the planet radius and its degree of irradiation. This correlation, first proposed in Guillot (2005), indicates that the inflated radius of some hot Jupiters can be at least partially explained by the enhanced insolation from their parent star. HATS-9b has a moderately high equilibrium temperature (=K) due to the small star-planet separation coupled to the large stellar radius, while HATS-10b has a more typical equilibrium temperature for a hot Jupiter (=K). According to the empirical relations proposed in Enoch et al. (2012), which give the radius of a giant planet from its equilibrium temperature and semi-major axis, HATS-9b and HATS-10b should have radii of 1.36 and 1.22, respectively. The observed radii are and below these values, which indicates that these planets are very compact given their irradiation levels and that thus additional variables must be responsible of setting the radii of short period giant planets.

One possible explanation is that HATS-9b and HATS-10b may have significant amounts of heavy elements in their cores. According to the interior models of Fortney et al. (2007), both planets will require a core mass of 60 M to explain their radii based on their masses, stellar host masses and orbital periods for an age of 4.5 Gyr. This explanation can be further motivated by the relatively high metallicity of their parent stars ( dex and dex, respectively). Several works (Guillot et al., 2006; Burrows et al., 2007; Enoch et al., 2011, 2012) have proposed a correlation between the inferred core mass of giant planets and the metallicity of the parent star. The principal idea behind the proposed correlation is that a more metal rich proto-planetary disk will be more efficient in creating massive cores following the core-accretion scenario of planetary formation. Even though this process is expected to occur in the formation and migration steps, the final relation between the stellar metallicity and the radius of giant planets is not at all clear and other phenomena can act in the opposite direction. As shown by Burrows et al. (2007), the presence of heavy elements in the atmosphere of young giant planets will increase its opacity, slowing the contraction and making the planetary radius more inflated than expected. Moreover, the validity of the proposed correlation has been put into question by the analysis Zhou et al. (2014a) who find no significant correlation between and [Fe/H] for the complete sample of detected giant TEPs.

The age of the system may be another important variable, since the radius of giant planets should undergo Kelvin-Helmholtz contraction as they age, controlled by their upper radiative atmosphere (Hubbard, 1977). Figure 8 presents the mass-radius diagram of transiting hot Jupiters having similar insolation levels to HATS-9b (1750K 1900K). This figure shows that in general the bloating of the atmosphere of strongly irradiated planets is prevented for more massive hot Jupiters. This correlation presents some outliers, with HATS-9b the most extreme one. A peculiarity of HATS-9b is the advanced age of the system (11 Gyr) contrasted with the ages of the rest of the planets in Figure 8 (5 Gyr). Among the complete sample of well characterized hot Jupiters, HATS-9b and CoRoT-17 b (10.7 1.0 Gyr) are the oldest systems known to have an age uncertainty better than 20%. Figure 9 shows the radius as function of age for hot Jupiters with 0.52 and orbital period P 10 days having age uncertainties smaller than 40%. Systems older than 3 Gyr exhibit the expected contraction of the envelope through time but most of them are systematically more inflated than expected from theoretical models of structure and evolution. By fitting a straight line through the planets with ages higher than 3 Gyr we obtain an empirical contraction function for hot Jupiters: , where is the age of the system in Gyr. The difference between the theoretical function and the empirical relation decreases with the age of the system and for the case of HATS-9b both functions are consistent with the observed values. The proposed empirical relation between the age of the system and the radius of the planet shown in Figure 9 supports the study of Burrows et al. (2007) where for young giant planets the higher opacity produced by heavy elements delays the contraction, while at later ages the higher mean molecular weight dominates and leads to smaller radii. However, in order to perform a precise study of the evolution of the the radii of giant extrasolar planets, particular models with the properties of each system should be constructed.

A possible confusing factor in Figures 7, 8 and 9 is the assumption of zero albedo and complete heat redistribution. The measurement of secondary transits on these systems in different wavelengths will be informative for explaining the departure of HATS-9b from the correlation. A more precise determination of the radius of HATS-9b is also required. The somewhat larger uncertainty in the radius is a result of the incomplete photometric follow-up for this system. The errors in the planet radius are governed at this point by the light-curve data, but future precise measurements of the transit of HATS-9b will be able to lower this uncertainty until it becomes dominated by the uncertainties on the stellar parameters.

Future precise RV measurements of HATS-10b are required to determine a more precise mass of the planet and to explain the high jitter measured with FEROS and Coralie with respect to Subaru/HDS. One possible explanation may be the presence of another planetary companion. Subaru/HDS observations, which don’t seem to show enhanced jitter, were performed in three continuous days, while Coralie and FEROS observations were separated by months and in this case the influence of a second more distant companion should be stronger. The jitter values quoted in Table 6 refer to RV uncertainties for each instrument that have to be added in quadrature to the formal RV errors in order for them to be consistent with the RV signal computed with the orbital parameters of the system.

4.1. K2 possibilities

Even thought HATS-9b and HATS-10b are located in the nominal coordinates of field 7 of K2, only HATS-9b falls on working silicon. A proposal to observe this star in short cadence was recently submitted. The high photometric precision of K2 will allow us to estimate a much more precise radius for HATS-9b, which will help us in determining if this planet is a true outlier in the correlation between planet radius, equilibrium temperature and planet mass. The high insolation of this planet makes it a very good target for measuring secondary transits and phase curve variations with K2, which will allow us to estimate the albedo and provide a more reliable estimate of its equilibrium temperature. Figure 10 shows a measure of the reflected light signature , , for hot Jupiters observed by Kepler as a function of planetary radius. From this Figure we can see that the potential of detecting reflected light signatures of HATS-9b is high and its amplitude should be similar to the one of the giant planets observed by Kepler so far. Other subtle photometric effects, like ellipsoidal variations, Doppler beaming and the measurement of asteroseismological frequencies, if present, will also be very valuable for the detailed characterization of this particular planet.

Figure 7.— (Left): Mass-radius diagram of giant TEPs. HATS-9b is marked with a filled square and HATS-10b with a filled triangle. Isodensity curves are plotted with dashed lines for gr cm and the 4.5 Gyr isochrones (Fortney et al., 2007) for core masses of 0 and 100 M with solid lines. (Right): Equilibrium temperature versus radius diagram for giant TEPs. Again, HATS-9b is marked with a filled square and HATS-10b with a filled triangle.
Figure 8.— Mass-radius diagram of giant TEPs having similar insolation levels to HATS-9b (1750K T1900K). HATS-9b is marked with a triangle. Filled symbols are colored according to the metallicity of the host star. HATS-9b does not follow the correlation formed by the other hot Jupiters with similar irradiation levels.
Figure 9.— Radius as function of the age of the system for hot Jupiters having , P10 days and an age estimation with a precision better than 40%. Green lines are the theoretical models of Fortney et al. (2007) for , AU and a core mass of 0 (dashed) and 50 (solid) times the mass of the earth. The red line is an empirical relation computed with these data points. HATS-9b is marked with a triangle. Hot Jupiters older than 3 Gyr follow the contraction of their radius over time but the observed contraction rate is steeper than the one predicted from the theoretical models. The theoretical radii for Hot Jupiters with ages greater than  10 Gyr (like HATS-9b) is consistent with the observations.
Figure 10.— Reflected light signature as function of the planetary radius for the hot Jupiters observed with Kepler. The symbols are colored according to the planetary mass. HATS-9b is marked with a triangle. Given that the photometric precision of K2 is similar to the one of the original Kepler mission, phase curve variations and secondary the secondary transit of HATS-9b should be measured by K2.
Development of the HATSouth project was funded by NSF MRI grant NSF/AST-0723074, operations have been supported by NASA grants NNX09AB29G / NNX12AH91H and internal Princeton funds. Follow-up observations receive partial support from grant NSF/AST-1108686. A.J. acknowledges support from FONDECYT project 1130857, BASAL CATA PFB-06, and project IC120009 “Millennium Institute of Astrophysics (MAS)” of the Millenium Science Initiative, Chilean Ministry of Economy. R.B. and N.E. are supported by CONICYT-PCHA/Doctorado Nacional. R.B. and N.E. acknowledge additional support from project IC120009 “Millenium Institute of Astrophysics (MAS)” of the Millennium Science Initiative, Chilean Ministry of Economy. V.S. acknowledges support form BASAL CATA PFB-06. K.P. acknowledges support from NASA grant NNX13AQ62G. This work is based on observations made with ESO Telescopes at the La Silla Observatory. This paper also uses observations obtained with facilities of the Las Cumbres Observatory Global Telescope. Work at the Australian National University is supported by ARC Laureate Fellowship Grant FL0992131. Operations at the MPG 2.2 m Telescope are jointly performed by the Max Planck Gesellschaft and the European Southern Observatory. The imaging system GROND has been built by the high-energy group of MPE in collaboration with the LSW Tautenburg and ESO. We thank Régis Lachaume for his technical assistance during the observations at the MPG 2.2 m Telescope. We thank Helmut Steinle and Jochen Greiner for supporting the GROND observations presented in this manuscript. We are grateful to P.Sackett for her help in the early phase of the HATSouth project. We thank Adam Burrows for his useful comments regarding the evolutionary models of hot Jupiters.


  1. affiliation: The HATSouth network is operated by a collaboration consisting of Princeton University (PU), the Max Planck Institute for Astronomy (MPIA), the Australian National University (ANU), and the Pontificia Universidad Católica de Chile (PUC). The station at Las Campanas Observatory (LCO) of the Carnegie Institute is operated by PU in conjunction with PUC, the station at the High Energy Spectroscopic Survey (H.E.S.S.) site is operated in conjunction with MPIA, and the station at Siding Spring Observatory (SSO) is operated jointly with ANU. Based in part on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. Based in part on observations made with the MPG 2.2 m Telescope at the ESO Observatory in La Silla. This paper uses observations obtained with facilities of the Las Cumbres Observatory Global Telescope. Based on observations obtained with the Apache Point Observatory 3.5-meter telescope, which is owned and operated by the Astrophysical Research Consortium.
  2. affiliation: Instituto de Astrofísica, Facultad de Física, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile; rbrahm@astro.puc.cl
  3. affiliation: Millennium Institute of Astrophysics, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile
  4. affiliation: Instituto de Astrofísica, Facultad de Física, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile; rbrahm@astro.puc.cl
  5. affiliation: Millennium Institute of Astrophysics, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile
  6. affiliation: Department of Astrophysical Sciences, Princeton University, NJ 08544, USA
  7. affiliation: Department of Astrophysical Sciences, Princeton University, NJ 08544, USA
  8. affiliation: Alfred P. Sloan Research Fellow
  9. affiliation: Packard Fellow
  10. affiliation: Observatoire Astronomique de l’Université de Genève, 51 ch. des Maillettes, 1290 Versoix, Switzerland
  11. affiliation: The Australian National University, Canberra, Australia
  12. affiliation: Department of Astrophysical Sciences, Princeton University, NJ 08544, USA
  13. affiliation: The Australian National University, Canberra, Australia
  14. affiliation: Max Planck Institute for Astronomy, Koenigstuhl 17, 69117, Heidelberg, Germany
  15. affiliation: Instituto de Astrofísica, Facultad de Física, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile; rbrahm@astro.puc.cl
  16. affiliation: Max Planck Institute for Astronomy, Koenigstuhl 17, 69117, Heidelberg, Germany
  17. affiliation: Instituto de Astrofísica, Facultad de Física, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile; rbrahm@astro.puc.cl
  18. affiliation: Millennium Institute of Astrophysics, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile
  19. affiliation: Max Planck Institute for Astronomy, Koenigstuhl 17, 69117, Heidelberg, Germany
  20. affiliation: Department of Astrophysical Sciences, Princeton University, NJ 08544, USA
  21. affiliation: Department of Astrophysical Sciences, Princeton University, NJ 08544, USA
  22. affiliation: Department of Earth and Planetary Sciences, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo 152- 8551
  23. affiliation: Perth Exoplanet Survey Telescope, Perth, Australia
  24. affiliation: Department of Astrophysical Sciences, Princeton University, NJ 08544, USA
  25. affiliation: Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138 USA
  26. affiliation: Max Planck Institute for Astronomy, Koenigstuhl 17, 69117, Heidelberg, Germany
  27. affiliation: The Australian National University, Canberra, Australia
  28. affiliation: Instituto de Astrofísica, Facultad de Física, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile; rbrahm@astro.puc.cl
  29. affiliation: Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138 USA
  30. affiliation: Hungarian Astronomical Association, Budapest, Hungary
  31. affiliation: Hungarian Astronomical Association, Budapest, Hungary
  32. affiliation: Hungarian Astronomical Association, Budapest, Hungary
  33. For HATSouth data we list the HATSouth unit and field name from which the observations are taken. HS-1 and -2 are located at Las Campanas Observatory in Chile, HS-3 and -4 are located at the H.E.S.S. site in Namibia, and HS-5 and -6 are located at Siding Spring Observatory in Australia. Each field corresponds to one of 838 fixed pointings used to cover the full 4 celestial sphere. All data from a given HATSouth field are reduced together, while detrending through External Parameter Decorrelation (EPD) is done independently for each unique field+unit combination.
  34. The median time between consecutive images rounded to the nearest second. Due to weather, the day–night cycle, guiding and focus corrections, and other factors, the cadence is only approximately uniform over short timescales.
  35. The RMS of the residuals from the best-fit model.
  36. IRAF 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.
  37. S/N per resolution element near 5180 Å.
  38. For Coralie and FEROS this is the systemic RV from fitting an orbit to the observations in Section 3.3. For ARCES and the du Pont Echelle it is the measured RV of the single observation. We do not provide this quantity for instruments for which only relative RVs are measured, or for WiFeS which was only used to measure stellar atmospheric parameters.
  39. For High-precision RV observations included in the orbit determination this is the RV residuals from the best-fit orbit, for other instruments used for reconnaissance spectroscopy this is an estimate of the precision.
  40. The zero-point of these velocities is arbitrary. An overall offset fitted independently to the velocities from each instrument has been subtracted.
  41. Internal errors excluding the component of astrophysical jitter considered in Section 3.3.
  42. footnotetext: Coralie observations acquired in August 2013 were contaminated with moonlight.
  43. Either HATS-9, or HATS-10.
  44. Barycentric Julian Date is computed directly from the UTC time without correction for leap seconds.
  45. The out-of-transit level has been subtracted. For observations made with the HATSouth instruments (identified by “HS” in the “Instrument” column) these magnitudes have been corrected for trends using the EPD and TFA procedures applied prior to fitting the transit model. This procedure may lead to an artificial dilution in the transit depths. For HATS-9 our fit is consistent with no dilution, for HATS-10 the HATSouth transit depth is % that of the true depth. For observations made with follow-up instruments (anything other than “HS” in the “Instrument” column), the magnitudes have been corrected for a quadratic trend in time fit simultaneously with the transit.
  46. Raw magnitude values without correction for the quadratic trend in time. These are only reported for the follow-up observations.
  47. ZASPE = **** routine for the analysis of high-resolution spectra (***cite***), applied to the FEROS spectra of HATS-9 and HATS-10. These parameters rely primarily on ZASPE, but have a small dependence also on the iterative analysis incorporating the isochrone search and global modeling of the data, as described in
  48. Computed following Valenti & Fischer (2005).
  49. Husser et al. (2013).
  50. From APASS DR6 for HATS-9, HATS-10 as listed in the UCAC 4 catalog (Zacharias et al., 2012)..
  51. From APASS DR6 for HATS-9, HATS-10 as listed in the UCAC 4 catalog (Zacharias et al., 2012)..
  52. From APASS DR6 for HATS-9, HATS-10 as listed in the UCAC 4 catalog (Zacharias et al., 2012)..
  53. From APASS DR6 for HATS-9, HATS-10 as listed in the UCAC 4 catalog (Zacharias et al., 2012)..
  54. From APASS DR6 for HATS-9, HATS-10 as listed in the UCAC 4 catalog (Zacharias et al., 2012)..
  55. YY++ZASPE = Based on the YY isochrones (Yi et al., 2001),  as a luminosity indicator, and the ZASPE results.
  56. In the case of the parameter is primarily determined from the global fit to the light curves and RV data. The value shown here also has a slight dependence on the stellar models and ZASPE parameters due to restricting the posterior distribution to combinations of ++ that match to a YY stellar model.
  57. Times are in Barycentric Julian Date calculated directly from UTC without correction for leap seconds. : Reference epoch of mid transit that minimizes the correlation with the orbital period. : total transit duration, time between first to last contact; : ingress/egress time, time between first and second, or third and fourth contact.
  58. Times are in Barycentric Julian Date calculated directly from UTC without correction for leap seconds. : Reference epoch of mid transit that minimizes the correlation with the orbital period. : total transit duration, time between first to last contact; : ingress/egress time, time between first and second, or third and fourth contact.
  59. Times are in Barycentric Julian Date calculated directly from UTC without correction for leap seconds. : Reference epoch of mid transit that minimizes the correlation with the orbital period. : total transit duration, time between first to last contact; : ingress/egress time, time between first and second, or third and fourth contact.
  60. Reciprocal of the half duration of the transit used as a jump parameter in our MCMC analysis in place of . It is related to by the expression (Bakos et al., 2010).
  61. Values for a quadratic law, adopted from the tabulations by Claret (2004) according to the spectroscopic (ZASPE) parameters listed in Table 5.
  62. As discussed in Section 3.3 the adopted parameters for all four systems are determined assuming circular orbits. We also list the 95% confidence upper limit on the eccentricity determined when and are allowed to vary in the fit.
  63. Term added in quadrature to the formal RV uncertainties for each instrument. This is treated as a free parameter in the fitting routine.
  64. Correlation coefficient between the planetary mass  and radius  estimated from the posterior parameter distribution.
  65. The Safronov number is given by (see Hansen & Barman, 2007).
  66. Incoming flux per unit surface area, averaged over the orbit.


  1. Aigrain, S., Hodgkin, S. T., Irwin, M. J., Lewis, J. R., & Roberts, S. J. 2015, MNRAS, 447, 2880
  2. Bakos, G., Noyes, R. W., Kovács, G., et al. 2004, PASP, 116, 266
  3. Bakos, G. Á., Kovács, G., Torres, G., et al. 2007, ApJ, 670, 826
  4. Bakos, G. Á., Torres, G., Pál, A., et al. 2010, ApJ, 710, 1724
  5. Bakos, G. Á., Csubry, Z., Penev, K., et al. 2013, PASP, 125, 154
  6. Bayliss, D., Zhou, G., Penev, K., et al. 2013, AJ, 146, 113
  7. Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977
  8. Boyajian, T. S., McAlister, H. A., van Belle, G., et al. 2012, ApJ, 746, 101
  9. Buchhave, L. A., Latham, D. W., Johansen, A., et al. 2012, Nature, 486, 375
  10. Burrows, A., Hubeny, I., Budaj, J., & Hubbard, W. B. 2007, ApJ, 661, 502
  11. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245
  12. Christensen-Dalsgaard, J., Kjeldsen, H., Brown, T. M., et al. 2010, ApJ, 713, L164
  13. Claret, A. 2004, A&A, 428, 1001
  14. Crossfield, I. J. M., Petigura, E., Schlieder, J. E., et al. 2015, ApJ, 804, 10
  15. Deming, D., Sada, P. V., Jackson, B., et al. 2011, ApJ, 740, 33
  16. Dopita, M., Hart, J., McGregor, P., et al. 2007, Ap&SS, 310, 255
  17. Eastman, J., Gaudi, B. S., & Agol, E. 2013, PASP, 125, 83
  18. Enoch, B., Collier Cameron, A., & Horne, K. 2012, A&A, 540, A99
  19. Enoch, B., Cameron, A. C., Anderson, D. R., et al. 2011, MNRAS, 410, 1631
  20. Esteves, L. J., De Mooij, E. J. W., & Jayawardhana, R. 2013, ApJ, 772, 51
  21. Everett, M. E., Barclay, T., Ciardi, D. R., et al. 2015, AJ, 149, 55
  22. Foreman-Mackey, D., Montet, B. T., Hogg, D. W., et al. 2015, ArXiv e-prints, 1502.04715
  23. Fortney, J. J., Marley, M. S., & Barnes, J. W. 2007, ApJ, 659, 1661
  24. Fraine, J., Deming, D., Benneke, B., et al. 2014, Nature, 513, 526
  25. Guillot, T. 2005, Annual Review of Earth and Planetary Sciences, 33, 493
  26. Guillot, T., Santos, N. C., Pont, F., et al. 2006, A&A, 453, L21
  27. Hansen, B. M. S., & Barman, T. 2007, ApJ, 671, 861
  28. Hartman, J. D., Bakos, G. Á., Torres, G., et al. 2011, ApJ, 742, 59
  29. Hartman, J. D., Bakos, G. Á., Béky, B., et al. 2012, AJ, 144, 139
  30. Hartman, J. D., Bayliss, D., Brahm, R., et al. 2014, ArXiv e-prints, 1408.1758
  31. Horch, E. P., Howell, S. B., Everett, M. E., & Ciardi, D. R. 2014, ApJ, 795, 60
  32. Howell, S. B., Everett, M. E., Sherry, W., Horch, E., & Ciardi, D. R. 2011, AJ, 142, 19
  33. Howell, S. B., Sobeck, C., Haas, M., et al. 2014, PASP, 126, 398
  34. Hubbard, W. B. 1977, icarus, 30, 305
  35. Husser, T.-O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6
  36. Jordán, A., Brahm, R., Bakos, G. Á., et al. 2014, AJ, 148, 29
  37. Kaufer, A., & Pasquini, L. 1998, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 3355, Optical Astronomical Instrumentation, ed. S. D’Odorico, 844–854
  38. Kovács, G., Bakos, G., & Noyes, R. W. 2005, MNRAS, 356, 557
  39. Kovács, G., Zucker, S., & Mazeh, T. 2002, A&A, 391, 369
  40. Mandel, K., & Agol, E. 2002, ApJ, 580, L171
  41. Mohler-Fischer, M., Mancini, L., Hartman, J. D., et al. 2013, A&A, 558, A55
  42. Noguchi, K., Aoki, W., Kawanomoto, S., et al. 2002, PASJ, 54, 855
  43. O’Donovan, F. T., Charbonneau, D., Mandushev, G., et al. 2006, ApJ, 651, L61
  44. Pál, A., Bakos, G. Á., Torres, G., et al. 2008, ApJ, 680, 1450
  45. Penev, K., Bakos, G. Á., Bayliss, D., et al. 2013, AJ, 145, 5
  46. Pollacco, D. L., Skillen, I., Collier Cameron, A., et al. 2006, PASP, 118, 1407
  47. Queloz, D., Mayor, M., Udry, S., et al. 2001, The Messenger, 105, 1
  48. Sanchis-Ojeda, R., & Winn, J. N. 2011, ApJ, 743, 61
  49. Sozzetti, A., Torres, G., Charbonneau, D., et al. 2007, ApJ, 664, 1190
  50. ter Braak, C. J. F. 2006, Statistics and Computing, 16, 239
  51. Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141
  52. Vanderburg, A., & Johnson, J. A. 2014, PASP, 126, 948
  53. Yi, S., Demarque, P., Kim, Y.-C., et al. 2001, ApJS, 136, 417
  54. Zacharias, N., Finch, C. T., Girard, T. M., et al. 2012, VizieR Online Data Catalog, 1322, 0
  55. Zhou, G., Bayliss, D., Penev, K., et al. 2014a, ArXiv e-prints, 1401.1582
  56. Zhou, G., Bayliss, D., Hartman, J. D., et al. 2014b, MNRAS, 437, 2831
This is a comment super asjknd jkasnjk adsnkj
The feedback cannot be empty
Comments 0
The feedback cannot be empty
Add comment

You’re adding your first comment!
How to quickly get a good reply:
  • Offer a constructive comment on the author work.
  • Add helpful links to code implementation or project page.