The following article is Open access

HATS-36b and 24 Other Transiting/Eclipsing Systems from the HATSouth-K2 Campaign 7 Program

, , , , , , , , , , , , , , , , , , , , , , and

Published 2018 February 16 © 2018. The American Astronomical Society. All rights reserved.
, , Citation D. Bayliss et al 2018 AJ 155 119 DOI 10.3847/1538-3881/aaa8e6

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

1538-3881/155/3/119

Abstract

We report on the result of a campaign to monitor 25 HATSouth candidates using the Kepler space telescope during Campaign 7 of the K2 mission. We discover HATS-36b (EPIC 215969174b, K2-145b), an eccentric ($e=0.105\pm 0.028$) hot Jupiter with a mass of $3.216\pm 0.062$ ${M}_{{\rm{J}}}$ and a radius of $1.235\pm 0.043$ ${R}_{{\rm{J}}}$, which transits a solar-type G0V star (V = 14.386) in a $4.1752$-day period. We also refine the properties of three previously discovered HATSouth transiting planets (HATS-9b, HATS-11b, and HATS-12b) and search the K2 data for TTVs and additional transiting planets in these systems. In addition, we also report on a further three systems that remain as Jupiter-radius transiting exoplanet candidates. These candidates do not have determined masses, however pass all of our other vetting observations. Finally, we report on the 18 candidates that we are now able to classify as eclipsing binary or blended eclipsing binary systems based on a combination of the HATSouth data, the K2 data, and follow-up ground-based photometry and spectroscopy. These range in periods from 0.7 day to 16.7 days, and down to 1.5 mmag in eclipse depths. Our results show the power of combining ground-based imaging and spectroscopy with higher precision space-based photometry, and serve as an illustration as to what will be possible when combining ground-based observations with TESS data.

Export citation and abstract BibTeX RIS

Original content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

1. Introduction

Transiting planet systems are the most valuable systems for exoplanet studies due to the wide range of characterization observations that can been made for them, both during transit and also in secondary eclipse. However, as the geometric probability of transit is low, large numbers of stars must be monitored in order to discover such systems. This task began over a decade ago with lenses-based wide-field ground-based surveys, most successfully HATNet (Bakos et al. 2004), WASP (Pollacco et al. 2006), and KELT (Pepper et al. 2007). These were followed by the space-based surveys CoRoT (Barge et al. 2008; Moutou et al. 2013) and Kepler (Borucki et al. 2010). In this paper, we combine data from two ongoing transit surveys: HATSouth (Bakos et al. 2013, hereafter HS) and K2 (Howell et al. 2014). While the HATSouth survey has been monitoring selected fields in the southern Hemisphere over the last six years, the K2 survey has been monitoring selected ecliptic plane fields since 2014. K2 Campaign 7 (C7) monitored a field centered at $\alpha ={19}^{{\rm{h}}}{11}^{{\rm{m}}}{19}^{{\rm{s}}}$, $\delta =-23^\circ 21^{\prime} 36^{\prime\prime} $ between 2015 October 4 to 2015 December 26. By chance, approximately 25% of this field had already been monitored by the HATSouth survey five years earlier between 2010 March to 2011 August. It was therefore possible for us to propose for K2 targets based on the HATSouth light curves that showed transit-like signals, and this was done as part of the K2 Guest Observer program (GO7067; PI: Bakos). In addition, we were also able to select as K2 targets the HATSouth transiting planets that had already been confirmed as bona fide planets, namely HATS-9b (Brahm et al. 2015), HATS-11b, and HATS-12b (Rabus et al. 2016); this formed the K2 Guest Observer program GO7066 (PI:Bakos). In total, 25 stars were monitored by K2 as part of these programs, which we collectively name "HS-K2C7" targets. Details for each of the HS-K2C7 candidates are set out in Table 1.

Table 1.  Basic Data for the HS-K2C7 Transiting Exoplanet Candidates

EPIC ID/K2 ID Mag R.A. Decl. Period ${R}_{p}/{R}_{\star }$ Tc
(HS-ID) $r$ band (J2000) (J2000) (days)   (BJD)
EPIC 217671466/K2-142 (HATS-9) $13.072$ ${19}^{{\rm{h}}}{23}^{{\rm{m}}}14\buildrel{\rm{s}}\over{.} 28$ $-20^\circ 09^{\prime} 58\buildrel{\prime\prime}\over{.} 7$ $1.9153073\pm 0.0000052$ $0.0725\pm 0.0041$ $2457380.70247\pm 0.00004$
EPIC 216414930/K2-143 (HATS-11) $13.865$ ${19}^{{\rm{h}}}{17}^{{\rm{m}}}36\buildrel{\rm{s}}\over{.} 24$ $-22^\circ 23^{\prime} 23\buildrel{\prime\prime}\over{.} 7$ $3.6191613\pm 0.0000099$ $0.1076\pm 0.0028$ $2457378.41910\pm 0.00007$
EPIC 218131080/K2-144 (HATS-12) $12.654$ ${19}^{{\rm{h}}}{16}^{{\rm{m}}}48\buildrel{\rm{s}}\over{.} 72$ $-19^\circ 21^{\prime} 21\buildrel{\prime\prime}\over{.} 2$ $3.1428330\pm 0.000011$ $0.0630\pm 0.0022$ $2457364.66541\pm 0.00012$
EPIC 215969174/K2-145 (HATS-36) $14.2$ ${19}^{{\rm{h}}}{25}^{{\rm{m}}}54\buildrel{\rm{s}}\over{.} 84$ $-23^\circ 12^{\prime} 10\buildrel{\prime\prime}\over{.} 0$ $4.1752387\pm 0.0000022$ $0.10966\pm 0.00067$ $2457360.390330\pm 0.000085$
EPIC 218210199 (HATS578-002) $13.872$ ${19}^{{\rm{h}}}{12}^{{\rm{m}}}51\buildrel{\rm{s}}\over{.} 48$ $-{19}^{^\circ }{12}^{{\prime} }47\buildrel{\prime\prime}\over{.} 3$ $1.9858275\pm 0.0000019$ $0.11261\pm 0.00155$ $2457170.51503\pm 0.00047$
EPIC 217231249 (HATS578-003) $13.881$ ${19}^{{\rm{h}}}{12}^{{\rm{m}}}55\buildrel{\rm{s}}\over{.} 08$ $-{20}^{^\circ }{56}^{{\prime} }22\buildrel{\prime\prime}\over{.} 0$ $4.8332040\pm 0.0000091$ $0.13627\pm 0.00202$ $2456589.65837\pm 0.00149$
EPIC 216579956 (HATS578-004) $15.217$ ${19}^{{\rm{h}}}{13}^{{\rm{m}}}20\buildrel{\rm{s}}\over{.} 64$ $-{22}^{^\circ }{05}^{{\prime} }40\buildrel{\prime\prime}\over{.} 4$ $0.7057527\pm 0.0000063$ $0.04817\pm 0.00197$ $2457332.86183\pm 0.00106$
EPIC 214652580 (HATS579-001) $14.312$ ${19}^{{\rm{h}}}{21}^{{\rm{m}}}10\buildrel{\rm{s}}\over{.} 44$ $-{25}^{^\circ }{57}^{{\prime} }59\buildrel{\prime\prime}\over{.} 8$ $8.9120143\pm 0.0000230$ $0.14571\pm 0.00100$ $2457320.69926\pm 0.00047$
EPIC 215626177 (HATS579-007) $12.971$ ${19}^{{\rm{h}}}{17}^{{\rm{m}}}44\buildrel{\rm{s}}\over{.} 16$ $-{23}^{^\circ }{50}^{{\prime} }51\buildrel{\prime\prime}\over{.} 2$ $2.0772332\pm 0.0000157$ $0.06156\pm 0.00357$ $2455943.10546\pm 0.00243$
EPIC 215716837 (HATS579-008) $15.027$ ${19}^{{\rm{h}}}{20}^{{\rm{m}}}29\buildrel{\rm{s}}\over{.} 04$ $-{23}^{^\circ }{40}^{{\prime} }28\buildrel{\prime\prime}\over{.} 6$ $8.6829946\pm 0.0000362$ $0.14714\pm 0.00160$ $2457285.32239\pm 0.00121$
EPIC 215101303 (HATS579-009) $14.945$ ${19}^{{\rm{h}}}{17}^{{\rm{m}}}38\buildrel{\rm{s}}\over{.} 76$ $-{24}^{^\circ }{56}^{{\prime} }09\buildrel{\prime\prime}\over{.} 0$ $15.2073500\pm 0.0000307$ $0.16459\pm 0.00179$ $2457047.10940\pm 0.00117$
EPIC 214912104 (HATS579-010) $14.210$ ${19}^{{\rm{h}}}{19}^{{\rm{m}}}18\buildrel{\rm{s}}\over{.} 12$ $-{25}^{^\circ }{21}^{{\prime} }21\buildrel{\prime\prime}\over{.} 2$ $13.3370855\pm 0.0000020$ $0.1900\pm 0.0010$ $2400010.45400\pm 0.00000$
EPIC 216442060 (HATS579-014) $12.747$ ${19}^{{\rm{h}}}{24}^{{\rm{m}}}10\buildrel{\rm{s}}\over{.} 08$ $-{22}^{^\circ }{20}^{{\prime} }27\buildrel{\prime\prime}\over{.} 2$ $5.2027430\pm 0.0000044$ $0.14502\pm 0.00086$ $2457265.46459\pm 0.00028$
EPIC 217149884 (HATS579-015) $14.204$ ${19}^{{\rm{h}}}{22}^{{\rm{m}}}38\buildrel{\rm{s}}\over{.} 28$ $-{21}^{^\circ }{05}^{{\prime} }01\buildrel{\prime\prime}\over{.} 9$ $16.6924091\pm 0.0000246$ $0.19491\pm 0.00131$ $2457248.29441\pm 0.00053$
EPIC 215358983 (HATS579-036) $13.788$ ${19}^{{\rm{h}}}{37}^{{\rm{m}}}55\buildrel{\rm{s}}\over{.} 20$ $-{24}^{^\circ }{23}^{{\prime} }05\buildrel{\prime\prime}\over{.} 9$ $6.4218981\pm 0.0000077$ $0.15515\pm 0.00084$ $2457001.60205\pm 0.00079$
EPIC 215234145 (HATS579-037) $13.140$ ${19}^{{\rm{h}}}{44}^{{\rm{m}}}53\buildrel{\rm{s}}\over{.} 16$ $-{24}^{^\circ }{38}^{{\prime} }57\buildrel{\prime\prime}\over{.} 4$ $1.2539910\pm 0.0000020$ $0.09970\pm 0.00226$ $2455879.30405\pm 0.00069$
EPIC 215353525 (HATS579-039) $14.508$ ${19}^{{\rm{h}}}{36}^{{\rm{m}}}18\buildrel{\rm{s}}\over{.} 72$ $-{24}^{^\circ }{23}^{{\prime} }47\buildrel{\prime\prime}\over{.} 0$ $0.9085948\pm 0.0000042$ $0.10881\pm 0.00469$ $2455844.28763\pm 0.00128$
EPIC 216231580 (HATS579-040) $14.825$ ${19}^{{\rm{h}}}{40}^{{\rm{m}}}37\buildrel{\rm{s}}\over{.} 56$ $-{22}^{^\circ }{43}^{{\prime} }18\buildrel{\prime\prime}\over{.} 0$ $3.9052839\pm 0.0000085$ $0.13867\pm 0.00141$ $2457247.84650\pm 0.00081$
EPIC 215714765 (HATS579-041) $14.349$ ${19}^{{\rm{h}}}{35}^{{\rm{m}}}41\buildrel{\rm{s}}\over{.} 64$ $-{23}^{^\circ }{40}^{{\prime} }42\buildrel{\prime\prime}\over{.} 5$ $6.6890947\pm 0.0000092$ $0.16474\pm 0.00313$ $2457086.78510\pm 0.00081$
EPIC 216562832 (HATS579-043) $15.714$ ${19}^{{\rm{h}}}{25}^{{\rm{m}}}51\buildrel{\rm{s}}\over{.} 60$ $-{22}^{^\circ }{07}^{{\prime} }30\buildrel{\prime\prime}\over{.} 3$ $0.6633600\pm 0.0000010$ $0.10784\pm 0.00111$ $2457303.61090\pm 0.00036$
EPIC 217393088 (HATS579-044) $15.363$ ${19}^{{\rm{h}}}{17}^{{\rm{m}}}45\buildrel{\rm{s}}\over{.} 24$ $-{20}^{^\circ }{39}^{{\prime} }15\buildrel{\prime\prime}\over{.} 6$ $1.3194747\pm 0.0000021$ $0.10742\pm 0.00135$ $2457235.14088\pm 0.00065$
EPIC 215816368 (HATS579-048) $15.642$ ${19}^{{\rm{h}}}{43}^{{\rm{m}}}39\buildrel{\rm{s}}\over{.} 36$ $-{23}^{^\circ }{29}^{{\prime} }13\buildrel{\prime\prime}\over{.} 2$ $10.1460176\pm 0.0000464$ $0.16817\pm 0.00220$ $2457263.96758\pm 0.00162$
EPIC 215474548 (HATS579-050) $15.881$ ${19}^{{\rm{h}}}{15}^{{\rm{m}}}34\buildrel{\rm{s}}\over{.} 92$ $-{24}^{^\circ }{08}^{{\prime} }34\buildrel{\prime\prime}\over{.} 9$ $1.2085467\pm 0.0000024$ $0.12716\pm 0.00114$ $2457322.57557\pm 0.00035$
EPIC 214512594 (HATS624-002) $15.774$ ${19}^{{\rm{h}}}{14}^{{\rm{m}}}33\buildrel{\rm{s}}\over{.} 72$ $-{26}^{^\circ }{18}^{{\prime} }32\buildrel{\prime\prime}\over{.} 6$ $1.8769829\pm 0.0000021$ $0.25501\pm 0.00292$ $2456676.27256\pm 0.00087$
EPIC 214439239 (HATS624-003) $15.453$ ${19}^{{\rm{h}}}{17}^{{\rm{m}}}11\buildrel{\rm{s}}\over{.} 76$ $-{26}^{^\circ }{29}^{{\prime} }21\buildrel{\prime\prime}\over{.} 4$ $0.6425813\pm 0.0000014$ $0.11122\pm 0.00094$ $2457316.03141\pm 0.00044$

Note. Period, ${R}_{p}/{R}_{\star }$, and Tc from combined HATSouth + K2 data. Mag $r$ band from APASS (Henden et al. 2009).

Download table as:  ASCIITypeset image

This paper is organized as follows. In Section 2, we describe the observations made for each HS-K2C7 target, including the initial HATSouth discovery photometry, reconnaissance spectroscopy and photometry, radial velocity measurements, and the K2 photometry. In Section 3, we analyze each HS-K2C7 candidate, including the discovery of a new transiting exoplanet (HATS-36b), refinement of the parameters of three known HATSouth transiting planets (HATS-9b, HATS-11b, HATS-12b), identification of three Jupiter-radii candidates, and classification of 18 HS-K2C7 candidates as eclipsing binaries or blended eclipsing binaries. Finally in Section 4, we discuss the results and implications of this joint ground/space-based photometry project.

2. Observations

2.1. HATSouth Photometry

The HATSouth global telescope network consists of six HS4 units, two each at Siding Spring Observatory (SSO), Las Campanas Observatory (LCO), and the High Energy Spectroscopic Survey (H.E.S.S.) site. Each HS4 unit holds four Takahashi astrographs (f/2.8, 18 cm apertures), which are each coupled to an Apogee U16M 4K × 4K CCD camera. Imaging is performed using Sloan r-band filters and with exposure times of 240 s. A full description of the HATSouth telescope network hardware and operations can be found in Bakos et al. (2013).

As part of the HATSouth survey (http://hatsouth.org/), we monitored a 64 square degree field (centered at ${19}^{{\rm{h}}}{30}^{{\rm{m}}}{00}^{{\rm{s}}}$, $-22^\circ 30^{\prime} 00^{\prime\prime} $) between 2010 March to 2011 August. In total, we obtained 10137 images, with a cadence of approximately 300 s. In Table 2, we set out a summary of the HATSouth observations of this field, including the number of images from each site, the observation dates, and the cadence.

Table 2.  Summary of HS-K2C7 Photometric Observations

Facility Target(s) Date Number Cadencea Filter
    Range of Images (s)  
HS-1 (LCO) all 2010 Mar–2011 Aug 4293 300 r
HS-3 (HESS) all 2010 Mar–2011 Aug 2556 303 r
HS-5 (SSO) all 2010 Sep–2011 Aug 3287 303 r
K2 Long Cadence (Kepler) all 2015 Oct 4–2015 Dec 2 3754 1800 Kep
K2 Short Cadence (Kepler) HATS-9 2015 Oct 4–2015 Dec 2 117401 60.5 Kep
K2 Short Cadence (Kepler) HATS-11 2015 Oct 4–2015 Dec 2 117601 60.5 Kep
PEST HATS-36 2013 Jul 1 137 130 RC
GROND HATS-36 2014 Jul 24 75 139 $g,r,i,z$
GROND HATS-36 2014 Jul 28 80 139 $g,r,i$
ANU2.3 m HATS579-037 (EPIC 215234145) 2012 Sept 8 138 134 I
PEST HATS579-037 (EPIC 215234145) 2013 Jul 7 194 130 RC
PEST HATS579-037 (EPIC 215234145) 2015 Aug 11 168 131 I
SWOPE HATS579-039 (EPIC 215353525) 2013 Aug 20 67 211 I

Note.

aThe mode time difference rounded to the nearest second between consecutive points in each light curve. Includes exposure time and overheads.

Download table as:  ASCIITypeset image

HATSouth raw images are reduced to light curves using an automated aperture photometric pipeline detailed in Penev et al. (2013). The light curves are detrended using External Parameter Decorrelation (EPD; Bakos et al. 2010) and the Trend Filtering Algorithm (TFA; Kovács et al. 2005). These light curves are then combed for transit-like features using the box-fitting least squares algorithm (Kovács et al. 2002) with the methodology described in Bakos et al. (2004). We found 25 candidates with transit-like signals that were also were on-silicon for the K2 Campaign 7 (see Section 2.5). We designate these as "HS-K2C7 candidates," and summarize them in Table 1. Three of these have already been published as confirmed transiting planets HATS-9b (Brahm et al. 2015), HATS-11b (Rabus et al. 2016), and HATS-12b (Rabus et al. 2016). We therefore do not discuss these further in this Section.

HATSouth light curves for all stars overlapping with the K2 Campaign 7, including the HS-K2C7 candidates, are publicly available at https://hatsurveys.org/data.html. By way of example, we present the HATSouth light curve for HATS-36b in Figure 1, which shows the 18 mmag transit-like dip when phase-folded at $P=4.1752387$ days. The per-point rms of the HATSouth light curve for HATS-36 from the three HATSouth stations is approximately 10mmag. This is primarily shot noise from the sky background as the star is faint (V = $14.386\pm 0.020$) as well as a small component of systematic noise on the order of 1 mmag. All HATSouth photometric data used in this paper is set out in Table 3. We note that as a result of the EPD and TFA detrending, and also due to blending from neighbors, the apparent transit depth in the HATSouth light curves is somewhat shallower than that of the true depth in the Sloan r filter (the apparent depth is typically 85% that of the true depth).

Figure 1.

Figure 1. Top panel: unbinned instrumental $r$ band light curve of HATS-36 folded with at the period $P=4.1752387$ days. The solid magenta line shows the best-fit transit model (see Section 3). Middle panel: a zoom-in on the transit; the dark blue filled points here show the light curve binned in phase using a bin-size of 0.002. Lower panel: the residuals to the best-fit model.

Standard image High-resolution image

Table 3.  Differential Photometry of HS-K2C7 Candidates

HS-K2C7 BJD Mag ${\sigma }_{\mathrm{Mag}}$ Mag Filter Instrument
ID     (2 400 000+)         (orig)    
HATS-36 55800.93945 0.00446 0.01045 r HS
HATS-36 55788.41475 −0.00216 0.01768 r HS
HATS-36 55755.01420 −0.02741 0.02600 r HS
HATS-36 55679.86013 0.00595 0.01239 r HS
HATS-36 55800.94295 −0.00508 0.00957 r HS
HATS-36 55725.78922 −0.00775 0.01555 r HS
HATS-36 55788.41822 −0.00621 0.01639 r HS
HATS-36 55755.01796 −0.02324 0.03203 r HS
HATS-36 55679.86505 −0.00567 0.01186 r HS
HATS-36 55725.79268 −0.02899 0.01553 r HS

Note. The data are also available on the HATSouth website at http://www.hatsouth.org.

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

2.2. Reconnaissance Spectroscopy

As part of the usual HATSouth program to follow-up transiting planet candidates, follow-up spectroscopy was obtained for the HS-K2C7 candidates, primarily using WiFeS (Dopita et al. 2007) on the ANU 2.3 m telescope at SSO. Full details for this observing program are given in Bayliss et al. (2013). In summary, we obtain a flux calibrated, high S/N, $R\,=\lambda /{\rm{\Delta }}\lambda =3000$ spectrum for each candidate in order to determine the spectral type and class of the host star. These spectra are compared with a grid of synthetic templates from the MARCS model atmospheres (Gustafsson et al. 2008) in order to estimate ${T}_{\mathrm{eff}}$. The results for the HS-K2C7 candidates are set out in Table 4. For scheduling reasons, we obtained reconnaissance spectra for two candidates (HATS579-014 and HATS579-036) using FEROS (Kaufer et al. 1998) on the MPG 2.2 m telescope at the ESO observatory in La Silla, Chile (LSO). Spectral parameters were derived from these spectra using the CERES code (Brahm et al. 2017a) and these are also tabulated in Table 4. We note that the spectrum of HATS579-014 showed that the candidate was a spectroscopic binary. We were not able to obtain a spectrum for two candidates, HATS579-050 and HATS624-003, as they were too optically faint (V = $16.221$ and V = $15.537$ respectively). However, based on their 2MASS colors, apparent magnitudes, and (for HATS579-050) proper motion, HATS579-050 is most likely a late K-dwarf and HATS624-003 is most likely an early G-dwarf.

Table 4.  HS-K2C7 Candidate Classifications

EPIC ID HS Recon. HS Recon. Sec. Eclipse Class. Comment
(HS-ID) ${T}_{\mathrm{eff}}$ (K) K ($\mathrm{km}\,{{\rm{s}}}^{-1}$) in K2 LC    
EPIC 217671466 (HATS-9) ... ... NO TEP HATS-9b(Brahm et al. 2015)
EPIC 216414930 (HATS-11) ... ... NO TEP HATS-11b(Rabus et al. 2016)
EPIC 218131080 (HATS-12) ... ... NO TEP HATS-12b(Rabus et al. 2016)
EPIC 215969174 (HATS-36) 6000 ± 300 $\lt 2.0$ NO TEP HATS-36b (this work)
EPIC 218210199 (HATS578-002) $6014\pm 300$ $-16.64\pm 0.63$ YES EB ...
EPIC 217231249 (HATS578-003) $5928\pm 300$ $33.04\pm 1.41$ YES EB Shallow K2 Section eclipse
EPIC 216579956 (HATS578-004) $5728\pm 300$ $\lt 2.0$ YES BEB Blended in wide K2 apertures
EPIC 214652580 (HATS579-001) $5861\pm 300$ $44.22\pm 0.81$ YES EB ...
EPIC 215626177 (HATS579-007) $6245\pm 300$ $\lt 2.0$ YES BEB Blended in wide K2 apertures
EPIC 215716837 (HATS579-008) $5334\pm 300$ $\lt 2.0$ YES EB Shallow K2 Section eclipse
EPIC 215101303 (HATS579-009) $6260\pm 300$ $20.38\pm 2.71$ NO EB ...
EPIC 214912104 (HATS579-010) $5945\pm 300$ $\lt 2.0$ YES BEB Blend with nearby EB (P = 43 days)
EPIC 216442060 (HATS579-014) ${SB}2$ $...$ YES EB Shallow K2 Section eclipse and OOT
EPIC 217149884 (HATS579-015) $5751\pm 300$ 13.497 ± 0.011 YES EB Shallow K2 Section eclipse
EPIC 215358983(HATS579-036) $6700\pm 100$ $22.49\pm 0.01$ YES EB Shallow K2 Section eclipse
EPIC 215234145 (HATS579-037) $5524\pm 300$ $\lt 2.0$ NO BEB Color-dependent depth, K2 OOT at P = ×2
EPIC 215353525 (HATS579-039) $5229\pm 212$ $\lt 2.0$ YES EB K2 OOT
EPIC 216231580 (HATS579-040) $5304\pm 300$ $\lt 2.0$ NO CAND P = 3.9 days candidate (depth = 1.9%)
EPIC 215714765 (HATS579-041) $5645\pm 300$ $54.92\pm 1.14$ YES EB Shallow K2 Section eclipse and OOT
EPIC 216562832 (HATS579-043) $6300\pm 300$ $\lt 2.0$ YES EB K2 Section eclipse and OOT at P = ×2
EPIC 217393088 (HATS579-044) $5945\pm 300$ $\lt 2.0$ NO CAND P = 1.3 days candidate (depth = 1.1%)
EPIC 215816368 (HATS579-048) $5320\pm 300$ $7.24\pm 2.21$ NO CAND P = 10.1 days candidate (depth = 2.5%)
EPIC 215474548 (HATS579-050) $...$ $...$ YES EB ...
EPIC 214512594 (HATS624-002) $4765\pm 300$ $\lt 2.0$ YES EB ...
EPIC 214439239 (HATS624-003) $...$ $...$ YES BEB Blended in wide K2 apertures

Download table as:  ASCIITypeset image

We also use WiFeS on the ANU 2.3 m telescope to obtain multi-epoch medium resolution ($R=\lambda /{\rm{\Delta }}\lambda =7000$) spectra to check for large amplitude (K > 2 $\mathrm{km}\,{{\rm{s}}}^{-1}$) radial velocity variations in phase with the photometric signal. This allows us to identify candidates that are eclipsing binaries without the need of more resource-intensive high-precision radial velocity monitoring. For those targets with multi-epoch spectra, we list the measured semi-amplitudes in Table 4 assuming circular orbits. Again for the candidate HATS579-036, FEROS was used instead of WiFeS to measure the radial velocity semi-amplitude. These reconnaissance radial velocity data are available on the HATSouth public archive at http://www.hatsouth.org.

2.3. Reconnaissance Photometry

In order to further rule out eclipsing binaries, and refine our ephemerides, we obtained ground-based photometric follow-up for four HS-K2C7 candidates. In this section, we detail all of these observations.

  • 1.  
    HATS-36 (EPIC 215969174): we obtained initial photometric follow-up on the night of 2013 July 1 using the 0.3 m Perth Exoplanet Survey Telescope (PEST) in Perth, Australia. For a full description of the PEST facility see Zhou et al. (2014) and the PEST website (http://pestobservatory.com). Imaging was carried out in the ${R}_{C}$ band with exposure times of 120 s. In total, 137 exposures of HATS-36 were taken. Data was reduced via aperture photometry as described in Zhou et al. (2014). The resulting light curve is plotted in Figure 2 and the data is provided in Table 3. A full transit is clearly detected with a depth and duration consistent with the HATSouth discovery data. These data allowed us to refine the transit ephemeris. The following year, we observed two consecutive partial transits of HATS-36b with the multiband GROND imager (Greiner et al. 2008) on the MPG 2.2 m telescope at LSO in Chile. On the night of 2014 July 24, we observed a transit egress (75 × 110 s images in $g,r,i,z$ band), while on the night of 2014 July 28, we observed a transit ingress (80 × 110 s images in $g,r,i$ band). These data were reduced to light curves via aperture photometry following the method set out in Mohler-Fischer et al. (2013). The light curves are plotted in Figure 2 and the data is set out in Table 3. These high-precision light curves confirmed the transit depth was color independent and were consistent with a transiting planet. Both the PEST and GROND light curves for HATS-36 are used in the global fitting described in Section 3.
  • 2.  
    HATS579-037(EPIC 215234145): we obtained reconnaissance photometric follow-up from the ANU2.3 m imaging camera at SSO on 2012 September 8 in $I$ band, which showed a "V"-shaped transit with a depth of 20 mmag (compared with the 10 mmag transit observed in the HATSouth discovery data). This color-dependent depth difference was confirmed with observations from the PEST 0.3 m telescope, which observed a 10 mmag transit in ${R}_{C}$ band on 2013 July 7, and a 20 mmag transiting ${I}_{C}$ band on 2015 August 11. These depth differences are statistically significant at a $\gt 3\sigma $-level. The photometric data for this candidate is set out in Table 3, while the follow-up light curves are plotted in Figure 3. The color-dependent depths are a strong indication that this candidate is an eclipsing binary.
  • 3.  
    HATS579-007 (EPIC 215626177): we obtained reconnaissance photometric follow-up from the ANU2.3 m imaging camera at SSO on 2013 May 24, the imaging camera on the SWOPE 1 m at LCO on 2013 August 21, and the PEST 0.3 m PEST telescope on 2014 July 3 and 2015 August 27. None of these observations showed a transit feature. It is likely that uncertainties in the HATSouth ephemeris for this candidate were responsible for us missing the transit event for this candidate during these photometric follow-up observations.
  • 4.  
    HATS579-039 (EPIC 215353525): we obtained reconnaissance photometric follow-up in $I$ band from the SWOPE 1 m at LCO on 2013 August 20. A deep (25 mmag) V-shaped transit was observed, which was consistent with the $r$ band discovery data. This observation confirmed the transit feature and refined the ephemeris for this candidate. The data for this observation are set out in Table 3.

Figure 2.

Figure 2. Left: unbinned ground-based follow-up transit light curves of HATS-36. The dates, filters, and instruments used for each event are indicated. The light curves have been detrended using the EPD process. Curves after the first are shifted for clarity. Our best fit is shown by the solid red lines. Right: residuals from the fits in the same order as the curves at left.

Standard image High-resolution image
Figure 3.

Figure 3. Left: unbinned ground-based follow-up transit light curves of HATS579-037. Plots as for Figure 2. The over-plotted model is based on a global fit of all available data, and accounts for differences in limb darkening between the R- and i-bands, but assumes all other light curve parameters are the same. The poor fit of this model shows clear differences in the observed transit depths between the R-band and i-band observations, which allow us to classify this candidate as a eclipsing binary.

Standard image High-resolution image

A summary of all photometric observations are set out in Table 2. All of the follow-up photometric data are set out in Table 3.

2.4. High-resolution Spectroscopy—Radial Velocities

We obtain radial velocity measurements for HS-K2C7 candidates that remain after the reconnaissance spectroscopy and photometry set out in Sections 2.2 and 2.3 respectively. An additional magnitude constraint of V < 14.5 is placed on candidates at this stage, as the radial velocity monitoring of candidates fainter than V = 14.5 is beyond the reach of most telescopes/instruments. Only in exceptional cases such as an M-dwarf host (e.g., HATS-6b; Hartman et al. 2015) do we attempt high-precision radial velocity measurements for such faint candidates. By this criteria, just seven candidates remained: three of which have already been published (HATS-9b, HATS-11b, and HATS-12b) and four which are presented below:

  • 1.  
    HATS-36 (EPIC 215969174): We measured the radial velocity of HATS-36 using FEROS on the MPG 2.2 m at LSO between 2013 July 16 and 2014 July 24. In total, 16 measurements were made spread over the phase of the photometric period ($4.1752$ day). These data were reduced using the CERES FEROS echelle spectrograph pipeline described in Brahm et al. (2017a). We find a radial velocity variation in phase with the photometric period and with an amplitude of K = $356.9\pm 3.6$ ${\rm{m}}\,{{\rm{s}}}^{-1}$, indicating the transiting companion was of a planetary mass. We present these radial velocity data in Table 5 and plot the data along with a best-fit circular orbit in Figure 4. These data are used in Section 3 to model the global parameters of the system—primarily determining the mass of the transiting planet.
  • 2.  
    HATS579-007 (EPIC 215626177): We obtained multiple high-resolution spectra of this target from several different instruments; however, the data showed no radial velocity variation above 10 ${\rm{m}}\,{{\rm{s}}}^{-1}$. Coupled with the lack of a transit in the reconnaissance photometry (see Section 2.3), we put the monitoring of this candidate on hold until K2 data became available (see Section 2.5).
  • 3.  
    HATS579-015 (EPIC 217149884): We obtained four FEROS observations for this candidate between 2016 May 18 and 2017 May 31. The observations show an in-phase radial velocity variation with $K=13.497\,\pm 0.011\,\mathrm{km}\,{{\rm{s}}}^{-1}$, indicating the companion is of stellar mass and the candidate is therefore an eclipsing binary.
  • 4.  
    HATS579-037 (EPIC 215234145): We obtained multiple high-resolution spectra of this target from several different instruments; however, the data showed no radial velocity variation above 10 ${\rm{m}}\,{{\rm{s}}}^{-1}$. Due to this lack of variation, along with the color-dependent transit depths (see Section 2.3), we put this candidate on hold until K2 data became available (see Section 2.5).

Figure 4.

Figure 4. Top panel: high-precision RV measurements for HATS-36 from MPG 2.2 m/FEROS, together with our best-fit circular orbit model. Zero phase corresponds to the time of mid-transit. The center-of-mass velocity has been subtracted. Middle panel: velocity $O-C$ residuals from the best-fit model. The error bars for each instrument include the jitter, which is varied in the fit. Bottom panel: bisector spans (BS), with the mean value subtracted. Note the different vertical scales of the panels.

Standard image High-resolution image

Table 5.  Relative Radial Velocities and Bisector Span Measurements of HATS-36

BJD RVa ${\sigma }_{\mathrm{RV}}$ b BS ${\sigma }_{\mathrm{BS}}$ Phase Instrument
(2, 450, 000+) (${\rm{m}}\,{{\rm{s}}}^{-1}$) (${\rm{m}}\,{{\rm{s}}}^{-1}$) (${\rm{m}}\,{{\rm{s}}}^{-1}$)      
6489.79353 −113.43 15.60 −3.0 17.0 0.486 FEROS
6492.81654 −302.03 17.90 42.0 18.0 0.210 FEROS
6841.72454 328.57 15.00 37.0 16.0 0.776 FEROS
6842.59071 47.37 20.00 4.0 21.0 0.983 FEROS
6844.76223 −57.83 13.30 31.0 15.0 0.503 FEROS
6846.77951 29.47 13.30 13.0 15.0 0.987 FEROS
6847.58993 −226.73 15.30 2.0 15.0 0.181 FEROS
6852.66976 −315.43 18.40 109.0 16.0 0.397 FEROS
6852.82588 −232.63 19.60 86.0 16.0 0.435 FEROS
6853.77517 294.77 21.50 58.0 19.0 0.662 FEROS
6855.59317 −140.43 16.30 10.0 17.0 0.097 FEROS
6856.83423 −291.33 16.40 −76.0 17.0 0.395 FEROS
6857.73972 199.17 14.20 38.0 15.0 0.612 FEROS
6858.64854 327.37 13.90 4.0 15.0 0.829 FEROS
6859.67881 −121.03 16.60 13.0 18.0 0.076 FEROS
6862.65477 424.37 26.90 141.0 26.0 0.789 FEROS

Notes.

aThe zero-point of these velocities is arbitrary. An overall offset ${\gamma }_{\mathrm{rel}}$ fitted separately to the FEROS velocities in Section 3 has been subtracted. bInternal errors excluding the component of astrophysical/instrumental jitter considered in Section 3.

Download table as:  ASCIITypeset image

2.5. K2 Photometry

The K2 mission (Howell et al. 2014) uses the Kepler space telescope (Borucki et al. 2010) to monitor selected fields in the ecliptic plane for campaigns of approximately 80 days each. Between 2015 October 4 and 2015 December 2, K2 Campaign 7 monitored a field centered on ${19}^{{\rm{h}}}{11}^{{\rm{m}}}{19}^{{\rm{s}}}$, $-23^\circ 21^{\prime} 36^{\prime\prime} $. As this K2 field overlapped with the previously monitored HATSouth field described in 2.1, we were able to obtain K2 data for the 25 HS-K2C7 candidates listed in Table 2 via the K2 Guest Observer programs GO7066 and GO7067. A summary of the K2 imaging is set out in Table 1.

The K2 data was made available on 2016 April 20. We used three different versions of the K2 light curves downloaded from the Mikulski Archive for Space Telescopes (http://archive.stsci.edu). We used the K2 PDC light curves, the light curves produced by the "self-flat-fielding" technique described in Vanderburg & Johnson (2014), and the K2 light curves from the everest open-source pipeline fully described in Luger et al. (2016, 2017). In addition, we used the K2SC light curves (Aigrain et al. 2016) that were provided to us upon request (B. Pope 2016, private communication).

All four of these data products are derived from the same raw pixel data. However, the different apertures and detrending techniques result in light curves with sometimes quite marked differences. We therefore analyzed all four in order to help categorize our HS-K2C7 candidates. Each HS-K2C7 candidate was detrended individually to take into account the fact that each light curve potentially contained a combination of K2 systematics, variability due to stellar rotation, and ellipsoidal variability. We detrended the light curves by first masking the transit/eclipse signal, and then flattening the curve using a fifth order SavitzkyGolay filter (Savitzky & Golay 1964), with the window length selected to detrend variations due to K2 systematics and variability due to stellar rotation, but to avoid flattening any potential ellipsoidal variation or secondary eclipse. The resulting 25 light curves are presented in Figure 5.

Figure 5.

Figure 5.

Two of the 24 K2 light curves for the HS-K2C7 candidates. The full set is available in the online Figure set. Light curves are from Everest (EV) or the K2 PDC pipeline (PDC). Classifications are transiting exoplanet (TEP), candidate (CAND), eclipsing binary (EB) or blended eclipsing binary (BEB). (The complete figure set (24 images) is available.)

Standard image High-resolution image

In most instances, we utilized the Everest light curves; however, for some light curves the Everest algorithm removed real astrophysical variability, so in those cases we used the PDC light curves. Features seen in these light curves, such as secondary eclipses and out-of-transit variability, are noted in Table 4.

Due to the space-based environment, the large aperture (1 m), and near continuous 80 day coverage, all of the K2 light curves we used are of very high precision—for K2 Campaign 7, the median 6.5 hr combined differential photometric precision (CDPP) for a Kp = 12 mag dwarf star was 120 ppm. With this very high precision, we are able to see features not visible or ambiguous in the HATSouth discovery light curves. Most importantly, we can search for secondary eclipses, out-of-transit ellipsoidal variation, and odd/even transit depth differences. These features, in an optical light curve and at significant amplitudes, are characteristic of eclipsing binary systems rather than transiting exoplanets.

The timing was such that the K2 observations and data followed after we had already completed the photometric and spectroscopic follow-up of the HS-K2C7 candidates detailed in Sections 2.3, 2.2, and 2.4. In this respect, some candidates had already been robustly identified as either transiting exoplanets or eclipsing binaries before the K2 data was analyzed. However, in other cases the K2 data were critical to our classification of the candidate. Here, we detail the findings for each HS-K2C7 candidate.

  • 1.  
    HATS-36 (EPIC 215969174): The phase-folded K2 light curve for HATS-36 is presented in Figure 6, and the data is tabulated in Table 3. It shows a 15 mmag U-shaped transit consistent with the discovery and follow-up photometry presented in Sections 2.1 and 2.3, respectively. There is no secondary eclipse, odd/even depth difference, or out-of-transit variation present to the limits of the photometry. This light curve is used in our global analysis of this newly discovered transiting exoplanet system in Section 3.
  • 2.  
    Known HATSouth Planets: The K2 data for HATS-9, HATS-11, and HATS-12 are consistent with the previous published exoplanet discoveries—i.e., there was no evidence of any secondary eclipses, odd/even depth differences, or out-of-transit variations. We analyze these light curves for additional planets, TTVs and phase modulations in Section 3.
  • 3.  
    HATSouth Candidates: HATS579-040, HATS579-044, and HATS579-048 all have K2 light curves that are consistent with the HATSouth discovery data with no sign of secondary eclipses, odd/even depth differences, or out-of-transit variations. The light curves are set out in Figure 5.
  • 4.  
    Eclipsing Binaries: The 18 remaining HS-K2C7 candidates are eclipsing binaries. There are 16 candidates that show secondary eclipses in the K2 light curves indicating they are eclipsing binaries. The details for each candidate are set out in Table 4, and the light curves are set out in Figure 5. For HATS579-009, we do not detect a secondary eclipse; however, we detected a high-amplitude (K = 20 $\mathrm{km}\,{{\rm{s}}}^{-1}$) in-phase radial velocity variation (see Section 2.4) indicating it is an eclipsing binary. Likewise, HATS579-037 does not have a detectable secondary eclipse; however, it had color-dependent transit depths (see Figure 3) and the K2 light curve also shows out-of-transit variation when phase-folded at ×2 the discovery period. Therefore, we classify this candidate as a eclipsing binary.

Figure 6.

Figure 6. K2 ${Kep}$ band light curve (everest) of HATS-36 folded with the period $P=4.1752387$ days resulting from the global fit described in Section 3. The top panel shows the full phase-wrapped light curve. The middle panel shows a zoom-in on the transit. The lower panel shows a zoom around phase 0.5 with no detection of a secondary eclipse. Points are individual K2 measurements. The black solid line is the best-fit global model described in Section 3.

Standard image High-resolution image

3. Analysis

In this section, we analyze the newly discovered transiting exoplanet HATS-36b, the three known HATSouth with K2 data (HATS-9b, HATS-11b, and HATS-12b), and the three HS-K2C7 targets that remain as transiting exoplanet candidates.

3.1. HATS-36b: A High-mass, Eccentric, Transiting Hot Jupiter

To derive the physical properties of HATS-36, we obtain initial stellar parameters from the high-resolution spectra of HATS-36 from FEROS, together with ZASPE (Brahm et al. 2017b). This provides a first estimate of the temperature (${T}_{\mathrm{eff}\star }$), surface gravity ($\mathrm{log}{g}_{\star }$), metallicity ($[\mathrm{Fe}/{\rm{H}}]$), and projected equatorial rotation velocity ($v\sin i$) of HATS-36. The ${T}_{\mathrm{eff}\star }$ and $[\mathrm{Fe}/{\rm{H}}]$ values are then used with the stellar density ${\rho }_{\star }$, determined from the combined light-curve and radial velocity analysis, to determine a first estimate of the stellar physical parameters following the method described in Sozzetti et al. (2007). We use the Yonsei–Yale isochrones (Y2; Yi et al. 2001) to derive the stellar mass, radius, and age that best fit our estimated ${T}_{\mathrm{eff}\star }$, $[\mathrm{Fe}/{\rm{H}}]$ and ${\rho }_{\star }$ values. We then determine a revised value of $\mathrm{log}{g}_{\star }$ and perform a second iteration of ZASPE holding $\mathrm{log}{g}_{\star }$ fixed to this value while fitting for ${T}_{\mathrm{eff}\star }$, $[\mathrm{Fe}/{\rm{H}}]$, and $v\sin i$. We again compare this new value of ${\rho }_{\star }$ to the Y2 isochrones to produce our final adopted values for the physical stellar parameters. Figure 7 shows the final $\mathrm{log}{g}_{\star }$ value plotted against ${T}_{\mathrm{eff}\star }$, with 1σ and 2σ confidence ellipsoids and the appropriate Y2 isochrones for various stellar ages. The final parameters indicate that HATS-36 is a G0V dwarf star (${T}_{\mathrm{eff}\star }$ = $5970\pm 160$ K, $\mathrm{log}{g}_{\star }$ = $4.400\pm 0.023$) with a mass and radius of ${M}_{\star }$ = $1.155\pm 0.039$ ${M}_{\odot }$ and ${R}_{\star }$ = $1.222\pm 0.029$ ${R}_{\odot }$, respectively. The metallicity is slightly above solar with $[\mathrm{Fe}/{\rm{H}}]$ = $0.280\pm 0.037$. The rotational velocity ($v\sin i$ = $4.98\pm 0.29$ $\mathrm{km}\,{{\rm{s}}}^{-1}$) and age ($0.62\pm 0.55$ Gyr) are not atypical for a star of this type and class. The full list of final stellar parameters are set out in Table 6.

Figure 7.

Figure 7. Comparison between the measured values of ${T}_{\mathrm{eff}\star }$ and ${\rho }_{\star }$ (from ZASPE applied to the FEROS spectra, and from our modeling of the light curves and RV data, respectively), and the Y2 model isochrones from Yi et al. (2001). The best-fit value (dark blue filled circle), and approximate 1σ and 2σ confidence ellipsoids are shown. The value from our initial ZASPE iteration is shown with an open red triangle. Black solid lines show the Y2 isochrones for ages indicated (0.2 to 10.0 Gyr). The green dashed lines show evolutionary tracks for stars with masses indicated in solar mass units.

Standard image High-resolution image

Table 6.  Stellar Parameters for HATS-36

Parameter         Value Sourcea
Astrometric properties and cross-identifications
    2MASS-ID... 2MASS 19255488-2312100  
    K2-ID... EPIC 215969174 (K2-145)  
    R.A. (J2000)... ${19}^{{\rm{h}}}{25}^{{\rm{m}}}54\buildrel{\rm{s}}\over{.} 84$ 2MASS
    Decl. (J2000)... $-23^\circ 12^{\prime} 10\buildrel{\prime\prime}\over{.} 0$ 2MASS
    ${\mu }_{{\rm{R}}.{\rm{A}}.}$ ($\mathrm{mas}\,{\mathrm{yr}}^{-1}$) $-1.5\pm 3.0$ UCAC4
    ${\mu }_{\mathrm{Decl}.}$ ($\mathrm{mas}\,{\mathrm{yr}}^{-1}$) $-7.5\pm 3.0$ UCAC4
Spectroscopic properties
    ${T}_{\mathrm{eff}\star }$ (K)... $6149\pm 76$ ZASPE
    $[\mathrm{Fe}/{\rm{H}}]$... $0.280\pm 0.037$ ZASPE
    $v\sin i$ ($\mathrm{km}\,{{\rm{s}}}^{-1}$)... $4.98\pm 0.29$ ZASPE
    ${\gamma }_{\mathrm{RV}}$ ($\mathrm{km}\,{{\rm{s}}}^{-1}$)... $-24.392\pm 31$ FEROS
Photometric properties
    B (mag)... $15.060\pm 0.030$ APASS
    V (mag)... $14.386\pm 0.020$ APASS
    g (mag)... $14.675\pm 0.010$ APASS
    r (mag)... $14.231\pm 0.010$ APASS
    i (mag)... $14.146\pm 0.050$ APASS
    Kep (mag)... $14.300\pm 0.030$ EPIC
    G (mag)... $14.148\pm 0.013$ Gaia DR1
    J (mag)... $13.181\pm 0.026$ 2MASS
    H (mag)... $12.855\pm 0.025$ 2MASS
    Ks (mag)... $12.809\pm 0.026$ 2MASS
Derived properties
    ${M}_{\star }$ (${M}_{\odot }$)... $1.222\pm 0.029$ Y2+${\rho }_{\star }$+ZASPE
    ${R}_{\star }$ (${R}_{\odot }$)... $1.155\pm 0.039$ Y2+${\rho }_{\star }$+ZASPE
    ${\rho }_{\star }$ (${\rm{g}}\,{\mathrm{cm}}^{-3}$)b ... $1.26\pm 0.27$ Light curves
    ${\rho }_{\star }$ (${\rm{g}}\,{\mathrm{cm}}^{-3}$)b ... $1.118\pm 0.091$ Y2+Light curves+ZASPE
    $\mathrm{log}{g}_{\star }$ (cgs)... $4.400\pm 0.023$ Y2+${\rho }_{\star }$+ZASPE
    ${L}_{\star }$ (${L}_{\odot }$)... $1.64\pm 0.17$ Y2+${\rho }_{\star }$+ZASPE
    MV (mag)... $4.24\pm 0.12$ Y2+${\rho }_{\star }$+ZASPE
    MK (mag,ESO) $2.897\pm 0.083$ Y2+${\rho }_{\star }$+ZASPE
    Age (Gyr)... $0.62\pm 0.55$ Y2+${\rho }_{\star }$+ZASPE
    AV (mag)... $0.232\pm 0.060$ Y2+${\rho }_{\star }$+ZASPE
    Distance (pc)... $962\pm 37$ Y2+${\rho }_{\star }$+ZASPE

Notes.

a2MASS (Skrutskie et al. 2006), APASS (Henden et al. 2009), Gaia (Lindegren et al. 2016), ZASPE = Zonal Atmospherical Stellar Parameter Estimator routine for the analysis of high-resolution spectra (Brahm et al. 2017b), Y2 isochrones (Yi et al. 2001). bIn the case of ${\rho }_{\star }$, we list two values. The first value is determined from the global fit to the light curves and radial velocity data, without imposing a constraint that the parameters match the stellar evolution models. The second value results from restricting the posterior distribution to combinations of ${\rho }_{\star }$+${T}_{\mathrm{eff}\star }$+$[\mathrm{Fe}/{\rm{H}}]$ that match to a Y2 stellar model.

Download table as:  ASCIITypeset image

To exclude blended eclipsing binary scenarios for HATS-36b, we carried out an analysis of all the data following the methodology set out in Hartman et al. (2012). We model the photometric data, including the K2 observations, as an eclipsing binary system blended with a third star. The stars in the model are constrained using the Padova isochrones (Girardi et al. 2000), and also must have a blended spectrum consistent with the determined atmospheric parameters. We also simulate composite cross-correlation functions (CCFs) and use them to predict radial velocities and bisector spans for each blend scenario. All blend models tested can be rejected with $\gt 6\sigma $ confidence based on the photometry alone. Moreover, none of the blend models tested would produce RV variations or non-variable bisectors consistent with the observations.

To determine the properties of HATS-36b, we globally model the photometric and spectroscopic data following Pál et al. (2008), Bakos et al. (2010), and Hartman et al. (2012). We fit Mandel & Agol (2002) transit models to the light curves, and a Keplerian orbit is fit to the radial velocity measurements presented in Section 2.4, allowing for RV jitter. For the ground-based light curves, we fixed the quadratic limb-darkening coefficients to tabulated values based on the stellar atmospheric parameters, while for the K2 light curve, we allowed the quadratic limb-darkening coefficients to vary in the fit. For the long-cadence K2 observations, we made use of the everest light curves, with the transit model numerically integrated over the exposure times. This light curve showed large amplitude quasi-periodic variations, likely due to a combination of low-frequency systematic errors in the K2 photometry, and the rotational modulation of stellar active regions on the surface of HATS-36. We discuss the stellar activity signature later in this section. In order to model these variations in our analysis, we made use of a Morlet-type wavelet basis and a low order polynomial. This model has the form:

Equation (1)

where t is the time of observation, T0 is a fixed reference epoch, ${\tau }_{j}$ are a fixed evenly spaced set of times for centering the wavelets, σ is a fixed wavelet width which we set equal to ${\tau }_{j+1}-{\tau }_{j}$, ν is fixed to the dominant quasi-periodic frequency, and ci, ${a}_{j,k}$, and ${b}_{j,k}$ are linearly optimized parameters in the model. For our analysis of HATS-36b, we adopted ${N}_{\mathrm{poly}}=2$, ${N}_{\mathrm{Morlet}}=7$, and ${N}_{\mathrm{harm}}=4$. The resulting K2 light curve is shown in Figure 6, while the data are part of the photometry set presented in Table 3. The model is fit simultaneously with the transit model as part of the global analysis.

We use a Differential Evolution Markov Chain Monte Carlo procedure to determine the posterior distribution of the parameters. We find a statistically significant non-zero eccentricity of $e=0.105\pm 0.028$; hence, we do not assume a circular orbit for this system. The planetary parameters resulting from this global fit are set out in Table 7. HATS-36b is a high-mass (${M}_{p}$ = $3.216\pm 0.062$ ${M}_{{\rm{J}}}$), Jupiter-sized (${R}_{p}$ = $1.235\pm 0.043$ ${R}_{{\rm{J}}}$) transiting planet with an orbital period of $4.1752387\pm 0.0000022$ days. It has a bulk density of $2.12\pm 0.20$ ${\rm{g}}\,{\mathrm{cm}}^{-3}$.

Table 7.  Parameters for the Transiting Planet HATS-36b

Parameter         Value
Light curve parameters
   P (days) ... $4.1752387\pm 0.0000022$
   Tc ($\mathrm{BJD}$) ... $2457360.390330\pm 0.000085$
   T14 (days) ... $0.14428\pm 0.00078$
   ${T}_{12}={T}_{34}$ (days) ... $0.01753\pm 0.00076$
   $a/{R}_{\star }$ ... $10.10\pm 0.28$
   $\zeta /{R}_{\star }$ ... ${15.749}_{-0.065}^{+0.086}$
   ${R}_{p}$/${R}_{\star }$ ... $0.10966\pm 0.00067$
   b2 ... ${0.204}_{-0.034}^{+0.029}$
   $b\equiv a\cos i/{R}_{\star }$ ... ${0.452}_{-0.040}^{+0.031}$
   i (deg) ... $87.61\pm 0.23$
Limb-darkening coefficients
   ${c}_{1},i$ (linear term) ... $0.2234$
   ${c}_{2},i$ (quadratic term) ... $0.3621$
   ${c}_{1},r$ ... $0.3057$
   ${c}_{2},r$ ... $0.3623$
   ${c}_{1},\mathrm{kep}$ ... 0.335 ± 0.051
   ${c}_{2},\mathrm{kep}$ ... 0.291 ± 0.090
RV parameters
   K (${\rm{m}}\,{{\rm{s}}}^{-1}$) ... $356.9\pm 3.6$
   e ... $0.105\pm 0.028$
   RV jitter (${\rm{m}}\,{{\rm{s}}}^{-1}$) ... $15\pm 44$
Planetary parameters
   ${M}_{p}$ (${M}_{{\rm{J}}}$) ... $3.216\pm 0.062$
   ${R}_{p}$ (${R}_{{\rm{J}}}$) ... $1.235\pm 0.043$
   $C({M}_{p},{R}_{p})$ ... $0.61$
   ${\rho }_{p}$ (${\rm{g}}\,{\mathrm{cm}}^{-3}$) ... $2.12\pm 0.20$
   $\mathrm{log}{g}_{p}$ (cgs) ... $3.718\pm 0.026$
   a (AU) ... $0.05425\pm 0.00043$
   ${T}_{\mathrm{eq}}$ (K) ... $1356\pm 30$
   Θ ... $0.2307\pm 0.0082$
   $\langle F\rangle $ (${10}^{8}$ $\mathrm{erg}\,{{\rm{s}}}^{-1}\,{\mathrm{cm}}^{-2}$) ... $7.64\pm 0.69$

Note. For each parameter, we give the median value and 68.3% (1σ) confidence intervals from the posterior distribution. Detailed notes on parameters can be found in Table 6 of Espinoza et al. (2016).

Download table as:  ASCIITypeset image

From our radial velocity monitoring, we robustly measure a non-zero eccentricity ($e=0.105\pm 0.028$). This is not entirely surprising, as it has been noted that massive hot Jupiters (${M}_{p}\gt 2\,{M}_{{\rm{J}}}$) often show non-zero eccentricity (Turner et al. 2016). Indeed, HATS-36b is very similar in orbital period and eccentricity to HATS-22b (Bento et al. 2017) and HAT-P-14b (Torres et al. 2010), both of which are massive hot Jupiters at 2.7 and 2.2 ${M}_{{\rm{J}}}$, respectively. Given the young age of the HATS-36b system ($0.62\pm 0.55$ Gyr), it is entirely possible that the planet has not had time to tidally circularize as is thought to have occurred for other lower mass and older hot Jupiter systems.

We take the raw (pre-detrending) HATS-36 K2 light curve (Vanderburg & Johnson 2014), remove the transit events, and examine the longer timescale variability using a Lomb Scargle analysis. From this analysis, we detect a strong modulation at 14.4 days with a peak-to-peak amplitude of approximately $0.5 \% $. Such a rotational period is typical for a star of the spectral type and class of HATS-36 (McQuillan et al. 2014). We confirm this rotational signal is also present in the HATSouth discovery light curve, again showing modulation at P = 14 days. The Lomb–Scargle power spectrum is presented in Figure 8 for both the K2 data and the HATSouth data. Assuming ${R}_{\star }=1.155\pm 0.039$, this rotation period would result in an equatorial rotational of ${v}_{e}=4.17\,\mathrm{km}\,{{\rm{s}}}^{-1}$, in contrast the spectroscopic derived $v\sin i=4.98\pm 0.29\,\mathrm{km}\,{{\rm{s}}}^{-1}$. The fact that the derived rotational period is $2-\sigma $ below the spectroscopic $v\sin i$ may be a sign that we have picked up a harmonic of the rotational period rather than the true rotational period. Alternatively, it may be due to non-equatorial spots and solar-like differential rotation.

Figure 8.

Figure 8. The power spectrum from a Lomb–Scargle analysis of the HATSouth photometry (blue solid line) and the K2 photometry (red dashed line) for HATS-36. The peak at ∼14 days is the likely rotational period of the star and is prominent in both data sets.

Standard image High-resolution image

3.2. HATS-9b,11b,12b

High-precision K2 data allowed us to revisit the three transiting exoplanets already published by the HATSouth team, namely HATS-9b (Brahm et al. 2015) and HATS-11b and HATS-12b (Rabus et al. 2016). We augment our originally published photometric and spectroscopic data with the new K2 photometry and re-run the global modeling with the methodology described in those papers and in Section 3.1. Here, we make use of the (Espinoza et al. 2016) light curves rather than the everest light curves, and we apply filter low-frequency variations due to stellar activity and instrumental errors prior to fitting the transit model. For HATS-9b and HATS-11b, short cadence data was available, and we made use of these observations, rather than the long cadence, in our analysis.

The results allow us to improve the precision for the planetary parameters for these systems, and we list these in Table 8. Most of the revisions to the planetary parameters are relatively minor. The largest change is for the radius of HATS-9b, which is revised upwards by almost 10%. This is primarily due to the limited photometric follow-up that was available when the parameters of HATS-9b were calculated in the original analysis of Brahm et al. (2015).

Table 8.  Updated Planet Parameters

Parameter         Previous Value This Work
HATS-9b (K2-142b)
   P (days) ... $1.9153073\pm 0.0000052$ 1.91531100 ± 0.00000094
   Tc ($\mathrm{BJD}$) ... $2456124.25896\pm 0.00086$ 2457380.702470 ± 0.000036
   T14 (days) ... $0.1457\pm 0.0024$ 0.14618 ± 0.00017
   ${R}_{p}$/${R}_{\star }$ ... $0.0725\pm 0.0041$ 0.08316 ± 0.00014
   i (deg) ... ${86.5}_{-2.5}^{+1.6}$ 88.94 ± 0.44
   ${M}_{p}$ (${M}_{{\rm{J}}}$) ... $0.837\pm 0.029$ 0.816 ± 0.038
   ${R}_{p}$ (${R}_{{\rm{J}}}$) ... $1.065\pm 0.098$ 1.1724 ± 0.0098
   ${\rho }_{p}$ (${\rm{g}}\,{\mathrm{cm}}^{-3}$) ... $0.85\pm 0.19$ 0.626 ${}_{-0.022}^{0.029}$
HATS-11b (K2-143b)
   P (days) ... $3.6191613\pm 0.0000099$ 3.6191634 ± 0.0000031
   Tc ($\mathrm{BJD}$) ... $2456574.9657\pm 0.0013$ 2457378.419100 ± 0.000069
   T14 (days) ... $0.1819\pm 0.0039$ 0.18202 ± 0.00025
   ${R}_{p}$/${R}_{\star }$ ... $0.1076\pm 0.0028$ 0.10721 ± 0.00028
   i (deg) ... $88.31\pm 0.86$ 89.03 ${}_{-0.47}^{+0.36}$
   ${M}_{p}$ (${M}_{{\rm{J}}}$) ... $0.85\pm 0.12$ 0.83 ± 0.10
   ${R}_{p}$ (${R}_{{\rm{J}}}$) ... $1.510\pm 0.078$ 1.487 ± 0.031
   ${\rho }_{p}$ (${\rm{g}}\,{\mathrm{cm}}^{-3}$) ... ${0.299}_{-0.050}^{+0.071}$ 0.315 ± 0.037
HATS-12b (K2-144b)
   P (days) ... $3.1428330\pm 0.000011$ 3.1428347 ± 0.0000022
   Tc ($\mathrm{BJD}$) ... $2456798.9556\pm 0.0012$ 2457364.66541 ± 0.00012
   T14 (days) ... $0.1899\pm 0.0031$ 0.18849 ± 0.00064
   ${R}_{p}$/${R}_{\star }$ ... $0.0630\pm 0.0022$ 0.06316 ± 0.00063
   i (deg) ... $82.7\pm 1.9$ 82.27 ± 0.64
   ${M}_{p}$ (${M}_{{\rm{J}}}$) ... $2.38\pm 0.11$ 2.390 ± 0.087
   ${R}_{p}$ (${R}_{{\rm{J}}}$) ... $1.35\pm 0.17$ 1.384 ± 0.059
   ${\rho }_{p}$ (${\rm{g}}\,{\mathrm{cm}}^{-3}$) ... ${1.19}_{-0.32}^{+0.54}$ 1.12 ± 0.14

Download table as:  ASCIITypeset image

For the case of HATS-11b, our global modeling, assuming a circular orbit, also finds a secondary eclipse with a planet-to-star flux ratio of 0.0032 ± 0.0012.

In addition to improving the parameters for these planets, we utilize the K2 light curves to search for additional transiting planets. The discoveries of WASP-47d and WASP-47e (Becker et al. 2015) show that hot Jupiters can have nearby planetary neighbors. After removing the transit events from the hot Jupiters, we search the K2 light curves of HATS-9, HATS-11, and HATS-12, using the BLS search algorithm as discussed in Section 2.1. We do not find any significant evidence for additional transiting planets in any of these systems.

Using the original light curves and the K2 data, we check for any changes in the timing of the transits for these exoplanets. We do this by fitting the best-fit transit model (using the parameters set out in Table 8) to each individual transit event and measuring the difference between the best-fit central transit time and the expected transit time from a purely Keplerian orbit. For all three systems, we find no evidence of any transit timing variations.

Finally, we analyze the K2 light curves for evidence for of additional variability. HATS-11 and HATS-12 only show long-term drift in the K2 data that is likely caused by systematics from the spacecraft. HATS-9 shows an additional modulation with a period of approximately 8.8 days and a peak-to-peak amplitude of approximately $0.1 \% $. If this modulation is indeed due to stellar rotation, it would imply an equatorial rotation of ${v}_{e}=8.6\,\mathrm{km}\,{{\rm{s}}}^{-1}$. Given the spectroscopic $v\sin i$ of HATS-9 is $4.58\pm 0.90\,\mathrm{km}\,{{\rm{s}}}^{-1}$, this would mean the spin axis of the star is inclined with respect to our line of sight, and that the planet is on a misaligned orbit.

3.3. Candidates

Three of the HATSouth candidates from the HS-K2C7 campaign remain viable transiting exoplanet candidates, but have not been confirmed via radial velocity measurements due to the faintness of the host stars. In this Section, we discuss the details for each candidate.

  • 1.  
    HATS579-040(EPIC 216231580): The transiting candidate has a period of P = 3.905 days and shows a 17 mmag "U"-shaped transit. The host star is a faint (V = $14.963$) G9 dwarf (${T}_{\mathrm{eff}}$ = $5304\pm 300$ K). Assuming a host star mass of 0.9 ${M}_{\star }$, and with a radial velocity semi-amplitude of K $\lt 2.0$ $\mathrm{km}\,{{\rm{s}}}^{-1}$, we can constrain the companion mass to be less than 15 ${M}_{{\rm{J}}}$ for an unblended single-host/single companion system. Gaia DR1 (Lindegren et al. 2016) shows a single source ($G=14.702$) having no neighbors within 15'' down to the Gaia DR1 magnitude limit ($G\sim 20$) and separation limit (∼1'').
  • 2.  
    HATS579-044(EPIC 217393088): The transiting candidate has a period of P = 1.320 days and shows a 10 mmag "U"-shaped transit. The host star is faint (V = $15.575$) G0 dwarf (${T}_{\mathrm{eff}}$ = $5945\pm 300$ K). Assuming a host star with 1.1 ${M}_{\star }$, and with a radial velocity semi-amplitude of K $\lt 2.0$ $\mathrm{km}\,{{\rm{s}}}^{-1}$, we can constrain the companion mass to be less than 11.5 ${M}_{{\rm{J}}}$ for an unblended single-host/single companion system. Gaia DR1 shows the host is a single source ($G=15.288$), having no neighbors within 15''.
  • 3.  
    HATS579-048(EPIC 215816368): The transiting candidate has a relatively long period of P = 10.148 days and shows a 20 mmag "U"-shaped transit. The host star is a faint (V = $15.835$) G9 dwarf (${T}_{\mathrm{eff}}$ = $5320\pm 300$ K). Assuming a host star mass of 0.9 ${M}_{\star }$, with a radial velocity semi-amplitude of K $7.24\pm 2.21$ $\mathrm{km}\,{{\rm{s}}}^{-1}$, we would determine a companion mass of 70 ${M}_{{\rm{J}}}$-approximately at the lower limit for H-burning. If confirmed, this would join the small population of known transiting brown dwarfs, and would follow the trend of having a longer orbital period than typical hot Jupiters (Bayliss et al. 2017). From the Gaia DR1, we see a primary source ($G=15.610$), with a neighbor at 9farcs5 ($G=18.182$). However, by analyzing multiple pixel-aperture sizes from the SFF K2 light curves (Vanderburg & Johnson 2014), we determine the detected transit signal originates from the primary candidate star rather than the fainter neighbor.

4. Discussion

This is the first time we have vetted HATSouth candidates using the high-precision photometry afforded by the Kepler telescope, although it has been done for HATNet under K2 program GO0116 (PI Bakos) resulting in the discovery of HAT-P-56b (Huang et al. 2015). In the cases of both HATS-36b and HAT-P-56b, the radial velocity semi-amplitudes are high (∼300 ${\rm{m}}\,{{\rm{s}}}^{-1}$), but the stellar jitter is also high (∼100 ${\rm{m}}\,{{\rm{s}}}^{-1}$), and thus the K2 data is especially helpful in robustly confirming the nature of the systems.

4.1. HATS-36b

HATS-36b is a hot Jupiter with a typical orbital period (P = $4.1752387\pm 0.0000022$ days). The star is active, which we see manifest in the the variability of the K2 light curve and the rotational periodicity recovered from the HATSouth light curve. It has a moderate but well-measured eccentricity, which is consistent with other high-mass hot Jupiters with measured eccentricities. Due to its high mass, HATS-36b lies in a relatively sparsely populated region of the mass-density relationship for gas giant exoplanets (see Figure 9). However, its bulk density fits well on the mass-density sequence of gas giants.

Figure 9.

Figure 9. Mass-density relationship of all well-characterized (density uncertainty <20%) giant exoplanets. Blue circles are data from NASA Exoplanet Archive (2016 November 23) and the red square is HATS-36b.

Standard image High-resolution image

4.2. Candidates

Secondary eclipses in K2 allow us to robustly rule out 17 of the 25 HS-K2C7 candidates. However, we have three candidates that remain active. All three candidates require future radial velocity monitoring in order to determine if they are transiting exoplanets/brown dwarfs or (blended) eclipsing binaries. However, such a task is extremely difficult to due to faintness of the host stars. HATS579-048 (EPIC 215816368) is the most promising, as we have an indication of a velocity variation of K = $7.24\pm 2.21$ $\mathrm{km}\,{{\rm{s}}}^{-1}$. However, the host star is also the faintest of the three candidates at V = $15.835$.

4.3. Outlook

This program shows the benefit of using high-precision K2 space-based photometry to vet candidates identified from ground-based surveys. This concept will naturally extend to the TESS mission (Ricker et al. 2014). The primary differences will be twofold. First, TESS will monitor most stars for only 27 days, about one-third the duration of the K2 campaigns. For many of our HS-K2C7 candidates, such duration would result in only one or two transits to be observed. In these cases, the value of combining the TESS data with ground-based monitoring is greatly enhanced. Second, the spatial resolution of TESS is just $21\buildrel{\prime\prime}\over{.} 1\,$ pixel−1, meaning many blended systems that can be resolved in K2, such as the blended eclipsing binaries presented in this work, will not be readily identifiable from TESS data alone. In these cases, ground-based data such as that from HATSouth ($3\buildrel{\prime\prime}\over{.} 7\,$ pixel−1) will be highly beneficial in resolving the nature of the systems.

Development of the HATSouth project was funded by NSF MRI grant NSF/AST-0723074. Operations have been supported by NASA grants NNX09AB29G, NNX12AH91H, and NNX17AB61G, and follow-up observations receive partial support from grant NSF/AST-1108686. This work has been carried out within the framework of the National Centre for Competence in Research PlanetS supported by the Swiss National Science Foundation. D.B. acknowledges the financial support of the SNSF. J.H. acknowledges support from NASA grant NNX14AE87G. A.J. acknowledges support from FONDECYT project 1171208, BASAL CATA PFB-06, and project IC120009 "Millennium Institute of Astrophysics (MAS)" of the Millennium Science Initiative, Chilean Ministry of Economy. N.E. is supported by CONICYT-PCHA/Doctorado Nacional. R.B. and N.E. acknowledge support from project IC120009 "Millennium Institute of Astrophysics (MAS)" of the Millennium Science Initiative, Chilean Ministry of Economy. V.S. acknowledges support form BASAL CATA PFB-06. A.V. is supported by the NSF Graduate Research Fellowship, grant No. DGE 1144152. This paper includes data collected by the K2 mission. Funding for the K2 mission is provided by the NASA Science Mission directorate. The K2 observations presented here were obtained through the GO program, with analysis supported by NASA grant NNX16AE68G. 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. We acknowledge the use of the AAVSO Photometric All-Sky Survey (APASS), funded by the Robert Martin Ayers Sciences Fund, and the SIMBAD database, operated at CDS, Strasbourg, France. 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 the MPG 2.2 m telescope support team for their technical assistance during observations.

Please wait… references are loading.
10.3847/1538-3881/aaa8e6