Arxiv:Astro-Ph/0609369 V2 15 Sep 2006 Ahntnd,54 Ra Rnhr.N,Wsigo DC, Washington NW, Rd

Home , HATNet Project

arXiv:astro-ph/0609369 v2 15 Sep 2006 ahntnD,54 ra rnhR.N,Wsigo DC, Washington NW, Rd. Branch 20015-1305 Broad USA 5241 DC, Washington n h aionaIsiueo ehooy ektm a be has time Califor Keck of Technology. University of NASA. the Institute by by granted California operated the and is which Observatory, 8hAeu,Clmu,O 43210 OH Columbus, Avenue, 18th hc soeae yteNtoa srnmclObservatory Astronomical National the by Japan. operated is which A970 USA 94720, CA bevtr,Kmgt,Aauh,Oaaa7903,Japan 719-0232, Okayama Asakuchi, Kamogata, Observatory, -58Bdps,Hungary. Budapest, H-1518 .O o 219 Box O. P. e tet abig,M 23,UA [email protected]. USA; 02138, MA Cambridge, Street, den nvriy a rnic,C 43,UA fischer@stars.sf USA; 94132, CA Francisco, San University, hs lnt r l ahrfit aigeetv fol- effective making faint, rather all are planets these gary pcnrainppr o hs bet,seteEx- the see objects, 2002a,b,c, follow- Encyclopedia these al. the Planets et for to Udalski trasolar papers references For confirmation OGLE-TR-10b, 132b; up are respectively). and 2004, (these 113b, 2003, project 111b, discov- OGLE 56b, candidates Five the of spectro- by follow-up through ered radii. found photometric and were and planets scopic masses transiting their unam- these of yield on transits exo- information known photometric 200 biguous their nearly the because among planets place special a occupy A--b AG-AIS O-EST XPAE TRANSITI EXOPLANET LOW-DENSITY LARGE-RADIUS, A HAT-P-1b: 10 13 † ⋆ 9 8 7 5 6 3 2 4 1 h 0cretykontastn xrslrplanets extrasolar transiting known currently 10 The rf eso etme 8 06–VRINAcpe o public for Accepted VERSION L – using 2006 typeset 18, Preprint September version Draft eateto ersra ants,Crei nttt o Institute Carnegie Magnetism, Terrestrial of Department eateto srnm,Uiest fClfri,Berkel California, of University Astronomy, of Department kym srpyia bevtr,Ntoa Astronomica National Observatory, Astrophysical Okayama ugra srnmclAscain 41Budapest, 1461 Association, Astronomical Hungarian eateto hsc srnm,SnFacsoState Francisco San Astronomy, & Physics of Department okl bevtr,Bdps,PO o 7 -15 Hun- H-1125, 67, Box P.O. Budapest, Observatory, Konkoly Fellow. Hubble eateto srnm,Eötvös Astronomy, Pf. Loránd of University, Department avr-mtsna etrfrAtohsc CA,6 Gar 60 (CfA), Astrophysics for Center Harvard-Smithsonian eateto srnm,Oi tt nvriy 4 West 140 University, State Ohio Astronomy, of Department ae npr nosrain banda h .M Keck M. W. the at obtained observations on part in Based ae npr ndt olce tteSbr Telescope, Subaru the at collected data on part in Based http://exoplanet.eu/ G. .P Stefanik P. R. .Bakos A. ´ lnt eintdHTP1,apast ea es slarge as headings: least Subject planet. at previously-known be any to than appears HAT-P-1b, designated planet, pnigtreyas peeie o h rni etrar center transit the for ephemerides years, three spanning 85 htmtyadsetocp,w n htteojc sa“hot a is object the that find M we spectroscopy, and photometry ansqec tr ihaeaot3Gra itneof distance a at 16 ADS Gyr system 3 star about double age the ∼ with of stars member main-sequence one transiting object 50A.Tetastsga a eido .62 asadde and days 4.46529 of period a has signal transit The AU. 1550 sn ml uoae eecpsi rzn n aai the Hawaii, and Arizona in telescopes automated small Using ◦ J . 1. ,tu rniigtesa tipc aaee .4o h st the of 0.74 parameter impact at star the transiting thus 9, n radius and NRDCINADSUMMARY AND INTRODUCTION 1,2 .W Noyes W. R. , A 1 T .Sato B. , E tl mltajv 12/14/05 v. emulateapj style X ∼ rf eso etme 8 06–VRINAcpe o publ for Accepted VERSION – 2006 18, September version Draft 1.36 binaries tr:idvda:AS142A lntr ytm:indivi systems: planetary – AB 16402 ADS individual: stars: .Z Stanek Z. K. R 13 7 J .A Johnson A. J. , h aetsasof stars parent The . 1 rvln na ri ihsmmjrai .5 Uadinclin and AU 0.055 axis semimajor with orbit an in traveling .Kov G. , 10 .L J. , acs ´ 3 .W Latham W. D. , az 8 ´ .P A. , su.edu ar ´ ABSTRACT edu. BINARY 5 nia 32, ey, .Papp I. , en of al ´ - f l 4,1 .W Marcy W. G. , Teohrtrewr rtdtce yrda velocity radial by detected nearby first were around measurements). surveys. three planets such by other transiting detected (The were known XO-1b, and Indeed, 5 TrES-1 2006b). stars, the al. et of Charbonneau are two telescopes, e.g. small a (see, using underway reason, surveys, al. this wide-field For et of special facilitated. as Sato number – are of well follow-up as – of are XO-1b determination types planets other (2005), parameter 149026b These accurate al. as HD et (2006). interest, Bouchy al. (2004), – et TrES- McCullough 189733b al. – (2000), et HD Alonso 209458b al. (2005), et HD – Henry (2000); 1 stars: al. bright et Charbonneau nearby, orbit exoplanets ftastn lnt rudbih tr.I currently 8 each discovery It monitors telescopes, the which automated stars. for wide-field of bright small aims 6 around that comprises planets survey transiting wide-field of a also is five other The difficult. somewhat lowup uigterfren fti paper. this of refereeing the during magni- large with 10th telescopes. follow-up space a for well-suited and orbits is ground-based HAT- thus HAT-P-1b, by and which as star, detection planet, tude This label first planet. hereafter the extrasolar we on transiting a (2002, report of al. we Net et Here instrumenta- Bakos see for flow, 2004). but data be and paper, observations, will forthcoming HATNet tion, a of in aspects Array Technical described Submillimeter Hawaii. the at at in (FLWO) (SMA) telescopes Observatory two Whipple and L. Arizona, F. the at fou spread telescopes with are network, instruments longitude-distributed two-station, The a they stars. in dips parent their photometric on 1% induce typically planetary the of detection through permit transits to enough bright stars 5000 † 5 15 14 h ANtproject, HATNet The h xsec fHTP1 a rtifre rma from inferred first was HAT-P-1b of existence The .S P. , ⋆ 1 .D Sasselov D. D. , to nApJ in ation icvr fTE- ODnvne l 06 a announced was 2006) al. et (O’Donovan TrES-2 of Discovery www.hatnet.hu e: ari ´ T nrdu,adsalri endensity, mean in smaller and radius, in C ∼ uie”pae ihms bu 0.53 about mass with planet Jupiter” 5 3 cadpoetdsprto of separation projected and pc 139 2453984 = .Sip B. & 0.Ti ytmi aro G0 of pair a is system This 402. la ais ae nadt set data a on Based radius. ellar 8 t f005mg rmfollow-up From mag. 0.015 of pth ANtpoethsdtce an detected has project HATNet .P Butler P. R. , GOEMME FASTELLAR A OF MEMBER ONE NG cto nApJ in ication 1 ◦ ocz ˝ .Torres G. , . × 9 + 397 15 4,1 8 ◦ ntae n20 yG. by 2003 in initiated fsy yial containing typically sky, of ul A--b– HAT-P-1b dual: N 9 .A Esquerdo A. G. , tr ∗ 1 .A Fischer A. D. , 4 . to about ation 62.The 46529. 14 transiting 1 , 6 .B, A. , ´ r 2 Bakos et al.

Fig. 1.— The left panel shows a 70′ × 70′ (300 pixels wide) stamp of a typical image taken by HAT-5. The stamp width is about 1/7th of the original frame. ADS 16402 is in the middle (the two components are merged into a bright star). Peaks of saturated stars are masked, and shown in white. The right panel is a stamp from an image taken with the FLWO 1.2 m telescope. It is 6.′8 × 6.′8 (600 pixels wide), and corresponds to the small box on the left panel. ADS 16402B(marked with the cross) is well separated from ADS 16402A. 0.6% (6 mmag) deep transit-like dip seen in the I-band of Fig. 1. Exposure times were 5 min at a cadence of light curve of the stellar system ADS 16402 (Aitken Dou- 5.5 min. Light curves were derived by aperture pho- ble Star Catalogue; Aitken & Doolittle 1932) obtained tometry for the 6400 stars in G205 bright enough to by the combination of the HAT-5 instrument at FLWO yield photometric precision of better than 2% (reach- and HAT-8 at the SMA. The members of this binary ing 0.3% in some cases). In deriving the light curves, stellar system are a pair of nearly identical G0V stars, we made use of the Trend Filtering Algorithm (TFA; ADS 16402A (I ≈ 9) and ADS 16402B (I ≈ 9.6) sep- Kovács, Bakos, & Noyes 2005) to correct for spurious arated by 11′′. 2. The system appears as a single, elon- trends in the data. This was particularly important for gated source in the HAT images because this separa- the light curve for ADS 16402, for which variable blend- tion is less than the 14′′pixel−1 image resolution. Ini- ing distorts the combined shape of the two close-lying tial follow-up spectroscopy at the FLWO 1.5 m telescope unresolved components; indeed the shallow transit sig- showed the absence of radial velocity variations of either nal in the combined light of the two stars is only barely star at a level of 1 kms−1 or larger, demonstrating that visible in the raw light curve, while it is readily seen after the source of the photometric dips could not be due to a application of TFA. common false positive, namely an M dwarf star transiting We then searched all the light curves for characteristic in front of one of the stars. Subsequent spectroscopy at transit signals, using the Box-fitting Least Squares (BLS; the Subaru telescope using N2K project observing time Kovács, Zucker, & Mazeh 2002) algorithm, which seeks (Fischer et al. 2004; Sato et al. 2005), and later at the for box-shaped dips in the parameter space of frequency, W. M. Keck telescope revealed sinusoidal velocity varia- transit duration, and phase of ingress. Candidate transit tions of ADS 16402B with an amplitude of 60m s−1, thus signals with the highest detection significance were then suggesting the existence of a 0.5MJ planet transiting that examined individually, to isolate those with the best com- star. After excluding blend scenarios, the combination of bination of stellar type (preferably main-sequence stars the photometric and spectroscopic data leads to the con- of spectral type mid-F or later), and depth, shape, and clusion that this object, HAT-P-1b, is indeed a planet, duration of transit. One of these was ADS 16402, for and has a radius comparable to that of HD 209458b, and which we identified a prominent boxcar-like periodicity a somewhat lower mean density. In the remainder of this with period 4.4656 days and transit depth in the com- paper we discuss all of the above points in more detail. bined light of the two stars of 6mmag. Portions of five transits were observed during the observing interval. The 2. OBSERVATIONS upper panel of Fig. 2 shows the combined data for these 2.1. HATNet Discovery Photometry five transits, with different symbols for observations from HAT-5 (Arizona) and HAT-8 (Hawaii). ADS 16402 (also BD+37 4734, α = 22h57m46s.8, δ = +38◦40′28′′, J2000.0) is a visual double star system that also bears the name HJ1832, as it was discovered 2.2. Rejection-mode Spectroscopic Observations by John Herschel (Herschel 1831). The system lies in Initial follow-up observations were made with the CfA HATNet survey field G205 (centered at α = 22h55m, Digital Speedometer (DS; Latham 1992) in order to char- δ = +37◦30′). 2059 observations were made of this field acterize the two individual stars comprising ADS 16402, by the HAT-5 telescope and 901 by HAT-8 between 2003 as well as to search for evidence that the transit signal October 13 and 2004 January 30. A small stamp taken could be due to the transit of a late M dwarf in front from a typical HATNet image is shown on the left panel of either of the stars, as has been found for many oth- Large Radius Transiting Planet in a Stellar Binary 3 erwise plausible planet transit candidates in the HAT- event associated only with the star ADS 16402B. Net database. These observations indicated that both Recent much higher quality observations using Kepler- sources indeed have values of Teff and log g consistent Cam on the FLWO 1.2 m telescope have confirmed this, with main-sequence late F or early G dwarfs. The DS revealing the egress from the transit on 2006 July 23, data revealed no sign of radial velocity variation in ei- and later, on 2006 Aug 31. A small stamp from an ther of the two components, with upper limit of detection FLWO 1.2 m image that clearly shows the two compo- about 1.7kms−1. For a period of ∼4.5d this limit corre- nents of ADS 16402AB is exhibited in the right panel sponds to a secondary mass of 15MJ, thus ruling out the of Fig. 1. The FLWO 1.2 m data were reduced using possibility that the transit signal is due to a low-mass our own aperture photometry software (fihat/fiphot, stellar companion transiting one of the two stars. Pál et al., in preparation). The photometry on individual frames was then transformed to a selected reference 2.3. Follow-up Photometry Observations frame using some one hundred stars as comparison in the 23′ × 23′ field provided by KeplerCam (we refer to this

8.515 step as smooth magnitude transformation). The result- HAT-5 data HAT-8 data 8.52 ing light curve for ADS 16402B already showed a smooth

8.525 transit curve with very high S/N, but low amplitude sys- inst I 8.53 tematic variations were present. As a ﬁnal correction, we

8.535 exploited the fact that ADS 16402B has an ideal com-

8.54 parison star, namely ADS 16402A, with almost identical -0.1 -0.05 0 0.05 0.1 Phase brightness, spectral type, and spatial position. The data 9.52 shown in the second and third panels of Fig. 2 thus cor- 9.53 respond to the diﬀerential magnitude of ADS 16402B to 9.54 ADS 16402A after the smooth magnitude transformation inst z 9.55 procedure (shifted back to an arbitrary z ≈ 9.53 out-of-

9.56 transit level). FLWO 1.2m 2006-Jul-23 data Fit 9.57 Unfortunately both FLWO 1.2 m transit events were 0.7 0.75 0.8 0.85 0.9 HJD - 2453939 partial, thus we seeked for further opportunities to ob- 9.52 serve a full transit. Such an event was observed re- 9.53 cently in I-band, on 2006 Sep 5, using the 60/90/180cm 9.54 Schmidt telescope at the Piszk´estet˝omountain station inst z 9.55 of Konkoly Observatory. Data are shown on the bottom 9.56 FLWO 1.2m 2006-Aug-31 data panel of Fig. 2, and were reduced in a similar manner to Fit 9.57 40.9 40.95 41 41.05 41.1 the FLWO 1.2 m data. HJD - 2453939 The HATNet, TopHAT, FLWO 1.2 m and Konkoly 9.51 Schmidt data together span 1058 days, or 237 periods. 9.52 Combining all the transit timings in these data we ﬁnd 9.53

inst improved ephemerides: T (HJD)= E + N · P , where z C 0 tr 9.54 E0 = 2453984.397 ± 0.009, and P = 4.46529 ± 0.00009 9.55 Konkoly/Piszkesteto 0.6m data Fit (days), and Ntr refers to the number of transits since the 9.56 45.3 45.35 45.4 45.45 45.5 45.55 45.6 recent Konkoly Schmidt observation that is a full transit HJD - 2453939 observation, hence it has the best timing measurement. Fig. 2.— The upper panel shows the HAT-5 (open circles) and HAT-8 (filled circles) TFA corrected light curve from FLWO, Ari- 2.4. High-precision Spectroscopic Observations zona and SMA, Hawaii, respectively. Altogether five transit events were observed successfully between the two sites. The data are In parallel to obtaining photometry follow-up data phased at period P=4.46529 days. Zero phase corresponds to the with the FLWO 1.2 m telescope, we obtained precise time of mid-transit using the ephemeris discussed in the text. The Doppler measurements of both ADS 16402A and B with second and third panels show the follow-up observations taken with KeplerCam on the FLWO 1.2 m telescope in z-band in 2006 July the High Dispersion Spectrograph (HDS) on the 8.2 m 23 and 2006 Aug 31, respectively. Unfortunately only the central Subaru telescope (Noguchi 2002) to confirm the plane- and egress portions of the transit were observable. The bottom tary nature of the HATNet photometric signal. panel shows the I-band data of a full-transit taken by the Konkoly The HDS spectral format spans 3500 − 6100A˚ over a Observatory 0.6 m Schmidt telescope at 2006 Sep 5. Superimposed ′′ on the lower three panels are our best fits (§ 6) to the transit light mosaic of two CCDs, and the 0. 8 slit yielded a resolu- curve. tion (λ/∆λ) of 55000. With only four observations we determined that the brighter component, ADS 16402A, After identification of the ADS 16402 system as a can- had low rms velocity scatter while component B exhib- didate for harboring a transiting extrasolar planet, we ited velocity variations that were well represented by a used the TopHAT telescope to search for evidence of a sine curve with an amplitude of about 60m s−1, a pe- transit. TopHAT is a 0.26 m aperture automated tele- riod equal to the photometric period, and a time of con- scope at FLWO, part of the HATNet, designed for follow- junction consistent with that predicted from the tran- up of transit candidates (Bakos et al. 2006). With its sit ephemeris. To improve phase coverage, nine obser- pixel scale of 2.2′′pixel−1 it is able to resolve the two stars vations of ADS 16402B were subsequently obtained at in the system. However, data were taken during the com- W. M. Keck telescope over eight nights in 2006 July using missioning phase of TopHAT with strongly variable fo- HIRES (Vogt et al. 1994). The spectrometer slit at Keck cus, and furthermore, conditions were non-photometric. is 0′′. 86 yielding a similar spectral resolution of about Nevertheless, we identified in the data a transit ingress 55000 with a spectral format from ∼ 3200 − 8800A.˚ The 4 Bakos et al.

2MASS Point Source Catalogue (Skrutskie et al. 2006), TABLE 1 and by deriving a weighted linear regression between J, Summary of stellar parameters for ADS 16402AB. H, K and I. The linear relation based on ∼ 400 stars is very well defined, and has an rms of 0.06mag. The I mag- Parameter ADS 16402A ADS 16402B nitude difference, ∆I = 0.53 ± 0.03, is more accurately determined than the individual magnitudes, because the From SME analysis (§ 3.1) Teff (K) 6047 ± 56 5975 ± 45 stars are nearly identical in spectral type and brightness log g 4.13 ± 0.10 4.45 ± 0.06 so that systematic errors are significantly canceled. v sin i(km s−1) 7.1 ± 0.3 2.2 ± 0.2 [Fe/H](dex) +0.12 ± 0.05 +0.13 ± 0.02 Simultaneous evolutionary track fitting(§ 3.2) Mass (M⊙) 1.16 ± 0.11 1.12 ± 0.09 0.14 0.10 Radius (R⊙) 1.23±0.10 1.15±0.07 0.17 log L∗(L⊙) 0.26 ± 0.15 0.18±0.14 MI (mag) 3.4 ± 0.3 3.7 ± 0.3 Age (Gyr) 3.6 Z 0.025 Other parameters I (mag) 9.035 ± 0.05 9.563 ± 0.05 ∆I (mag) 0.53 ± 0.03 22 Mean Dist. (pc) 139±19 250 Proj. Sep.(AU) 1550±210

Fig. 3.— Stellar evolution tracks (solid curved lines) for Subaru and Keck data were used to characterize more ADS 16402A and ADS 16402B corresponding to the best mutual fully the stellar properties of the system (§ 3), as well as fit (filled circles with error bars). The two stars are physically to obtain a definitive radial velocity orbit (§ 4). bound, are assumed to be coeval, and to have the same abundance, as determined by the SME analysis. Initial parameters from 3. STELLAR PROPERTIES OF THE ADS 16402 SYSTEM SME are shown with open circles with their respective (overlap- ping) error-boxes. The two tracks stop at 3.6 Gyr, so that the line 3.1. SME analysis of the individual stars joining their endpoints defines the local isochrone. The straight dashed lines plot the ZAMS for metallicity Z=0.020 and Z=0.025. We have carried out spectral synthesis modeling of the For clarity the ZAMS for Z=0.030 is not shown in the figure; it iodine-free Keck template spectra for both ADS 16402A is displaced by the same amount above the ZAMS for Z=0.025 as and B using the SME software (Valenti & Piskunov the latter is above the ZAMS for Z=0.020. Dotted lines illustrate lines of equal radius. The position of the Sun and its track to its 1996), plus the wavelength ranges and atomic line data current age are also shown. described by Valenti & Fischer (2005). Results for Teff , log g, v sin i and [Fe/H] for both stars are shown in the By exploiting the physical association of the two stars, first four lines of Table 1. Given the derived values of Teff their common abundances, age, and their known lumi- the stars appear to fit the MK classification of G0, rather nosity ratio, we carried out a simultaneous evolutionary than the F8 classification reported in the literature. track fitting. This in turn leads to a much better determination of the parameters of ADS 16402B than would 3.2. Simultaneous evolutionary track fitting be the case if it were a single star. The two stars are very likely to comprise a physical Initial parameters (Fig. 3, open circles) and their error pair, for the following reasons: First, they have com- boxes (solid line rectangles in Fig. 3) were based on the mon proper motion (Halbwachs 1986). Second, their atmospheric parameters obtained with the SME analysis relative apparent magnitudes compared to the abso- of the Subaru spectra. While the errors in Teff are taken lute magnitude expected for main-sequence stars with directly from Table 1, the absolute luminosities rely on Teff and log g given in Table 1 indicate that they are the estimated gravity and iteration with the mass de- at a common distance. This also implies that their termination, and this procedure is imprecise. However, true tangential velocities are very similar, and also that the ratio of the luminosities of the two stars is well con- their real separation is small (projected separation be- strained, and is used in initializing the isochrone fitting ing ∼1550 AU). Third, based on the six DS veloci- to both stars. ties for each star acquired over a time interval of 234 The value of [Fe/H] found from spectral synthesis days, they have common radial velocity within the er- corresponds to the range of metal abundance by mass −1 rors: Vrad,A = −3.43 ± 0.32 (rms) ± 0.14(sys)km s , Z=0.020–0.030. We have calculated stellar evolution −1 Vrad,B = −2.94 ± 0.56 (rms) ± 0.25(sys)km s . Finally, tracks of ADS 16402A and ADS 16402B using Z in this they have identical metallicity within uncertainties. The range (Fig. 3). The straight line that would connect the physical companionship implies that the two stars are centers of the two error boxes is steeper than the theo- coeval, and provides the opportunity to perform simul- retical isochrones (the straight line that would connect taneous evolutionary track fitting. the end tips of the two evolutionary tracks) for any Z in Furthermore, because of their common distance, the the range Z=0.020–0.030. All possible solutions are some difference in their apparent magnitudes yields an accu- sort of compromise within the observational error boxes, rate measurement of the ratio of their luminosities. The imposed by the requirement for equal age, and with the I magnitudes of ADS 16402A and ADS 16402B respec- obvious additional constraint that each track must begin tively are 9.03 and 9.56. These values were determined by on the ZAMS appropriate to the chosen value of Z. For cross-correlating the I magnitudes for the Landolt (1992) different Z’s, different mass ranges sweep the two error standard stars with the J, H and K magnitudes of the boxes – smaller ones for Z=0.020 and progressively larger Large Radius Transiting Planet in a Stellar Binary 5

Subaru and Keck were typically 500 s, yielding an S/N TABLE 2 of about 150. Radial Velocities for ADS 16402B. Emission in the core of the Ca H&K lines was mea- HJD - 2453000 RV Uncert. Observatory sured in each of the Keck spectra to assess the chromo- (days) (m s−1) (ms−1) spheric activity of ADS 16402A and B. With five spectra for ADS 16402A, we measure an average value for ′ − 897.11285 -2.89 5.00 Subaru SHK =0.16 corresponding to log RHK = 4.923 for the 899.12678 52.53 4.77 Subaru A component and an average value of nine spectra yields ′ 900.12008 -34.51 4.65 Subaru SHK = 0.14 and log RHK = −5.03 for the B compo- 901.10260 -37.51 4.88 Subaru nent. These values indicate low chromospheric activity. 927.06848 -37.76 3.84 Keck Our estimate for intrinsic astrophysical velocity jitter is 927.96558 -28.27 3.90 Keck −1 931.03690 -12.10 3.83 Keck 3.7ms for both stars (Wright 2005). 931.94061 -46.59 3.92 Keck Radial velocity measurements from Subaru and Keck 932.03590 -42.95 3.83 Keck were used together to derive an orbital fit. The reference 932.99985 3.53 3.85 Keck 933.92455 66.16 4.04 Keck velocities were those from Keck, and an offset (∆v) be- 934.90368 49.20 4.92 Keck tween Subaru and Keck was also fitted for in the orbital 934.90692 49.21 4.00 Keck solution. Velocity measurements are given in Table 2 along with observation times and errors. The velocity −1 for higher Z, with a largely normal distribution. The offset of ∆v = 14.5ms has been already applied to observational error box of ADS 16402B from the SME the Subaru data shown in the table, and the expected analysis (Fig. 3) overlaps with the ZAMS for any Z, and velocity jitter of 3.7ms−1 has already been added in naturally, the star cannot lie below the ZAMS for the quadrature to the uncertainties listed in the table. The accepted abundance, so this constrains the possible so- photometric period of 4.46529 days, and the mid-transit lution from the lower bounds of log g and log L. The time from the photometric ephemerides were fixed in this distribution of possible stellar radii is asymmetric due to model. We also assumed zero eccentricity based on the these constraints and this leads later to asymmetric error theoretical expectation of circular orbits for periods as bars of the planetary radius. short as this. In summary, the fitted parameters were the Our best fit theoretical isochrone is for Z=0.025, and γ velocity of the Keck data points, the K semi-amplitude, roughly bisects the line connecting the centers of the and the ∆v velocity difference between Keck and Sub- SME error boxes. This gives masses of 1.12 ± 0.09 aru. To determine the uncertainties, 100 Monte Carlo 0.10 M⊙ and 1.16 ± 0.11 M⊙ and radii of 1.15±0.07R⊙ and fitting trials were run. In each trial, the best fit Kep- 0.14 1.23±0.10R⊙, for ADS 16402B and ADS 16402A respec- lerian model was subtracted from the velocities and the tively, and an age of 3.6 Gyrs for the system (see Table residuals were shuffled and added back with replacement 1), and we use these in the further analysis. to the theoretical velocities. From these trials, we find −1 The absolute I magnitude of ADS 16402B is MI = a velocity semi-amplitude, K = 60.3 ± 2.1ms , an rms 3.74±0.3 mag (with bolometric correction of +0.55 mag). of 5.11ms−1 and a reduced χ2 = 1.54. The final fit is From its apparent I magnitude of 9.56 ± 0.06mag (see over-plotted on the RV data in the upper panel of Fig. 4. above), and by neglecting reddening (being a near-by The lower panel shows the residuals. ◦ 24 source at b = 20 ) we derive a distance of 146±21pc. The Adopting a stellar mass of 1.12 ± 0.09M⊙, we derive 21 same calculation for ADS 16402A yields 131±17pc. We arel =0.0551 ± 0.0015 AU and m sin i =0.53 ± 0.04MJ. ◦ ◦ estimate the distance of the ADS 16402AB system as the Using an inclination iP = 85.9 ± 0.8 (§ 6 below), this ′′ average of the two: 139 ± 20pc. The 11. 2 separation of yields a planetary mass MP =0.53±0.04MJ. The orbital the two components then implies a projected separation elements and uncertainties are summarized in Table 3. of 1550 ± 250 AU. There is a small but non-negligible systematic trend All derived stellar parameters, along with absolute and in the residuals of the orbital fit (lower panel of Fig. 4). apparent I magnitudes, magnitude difference, mean dis- We note that the fit is improved if the eccentricity and tance, and projected separation of the stellar components ω are allowed to float, and we get e = 0.09 ± 0.02, of ADS 16402 are given in Table 1. and ω = 80◦.7 ± 7◦.6 with rms of 3.2ms−1 and reduced χ2 = 0.53. The non-zero eccentricity is only suggestive, 4. SPECTROSCOPIC ORBIT OF ADS 16402B however, because of the i) small number of data-points, The Subaru HDS Doppler observations were obtained and ii) the unrealistically small reduced χ2. Nevertheless, through an Iodine absorption cell (in front of the spec- implications of this finding are potentially very impor- trometer slit for this observing run) to provide a fidu- tant, and we elaborate further on them in the discussion cial wavelength scale and to model the instrumental (§ 7). point spread function (Kambe et al. 2002). We measured 5. the star’s radial velocity from each high-resolution spec- EXCLUDING BLEND SCENARIOS trum, using a synthetic stellar template that was cre- We investigated the possibility that the measured ated by modifying a close-match spectrum from a li- radial-velocity variations are due not to a planetary com- brary of high resolution template spectra obtained at panion but rather to distortions of spectral line profiles Keck (Johnson et al. 2006). Doppler velocities at Keck arising from contamination of the spectrum by the pres- are also modeled from observations taken through an Io- ence of an unresolved binary companion to ADS 16402B, dine cell, however, a standard template observation (ob- with an orbital period of ∼4.4653 days and an ampli- tained without the Iodine cell) was used for the reference tude of several tens of kms−1. However, this would give spectrum (Butler et al. 1996). Exposure times at both rise to time-varying asymmetry in the spectral line bi- 6 Bakos et al.

Fig. 5.— Comparison between the relative radial velocities Fig. 4.— Radial Velocities from Subaru and Keck, folded with and bisector spans of ADS 16402A (left panel) and ADS 16402B the period of 4.46529 days. Phase 0 corresponds to the predicted (right panel), based on the Subaru spectra. ADS 16402A shows no time of mid-transit from the orbital ephemeris given in text. Open variation in either the velocities or the bisectors spans (left panels), 1 circles represent Subaru data (after the 14.5m s− shift described while ADS 16402B shows obvious velocity variations but does not in the text, whereas filled circles show Keck measurements. show the variation in the bisector spans that would be expected if the velocity signal were the result of a blend with an eclipsing binary. TABLE 3 nary, with the requirement that the brightness of the Parameters of the HAT-P-1b secondary be no more than ∼10% that of the primary planetary system. (to have escaped detection in the DS spectra) and that the light curve of the eclipsing binary, attenuated by the Parameter Value much brighter primary, should be indistinguishable from the observed light curve in the FLWO 1.2 m data. How- Period (d) 4.46529 ± 0.00009a a ever, we have found it impossible to construct an accept- Tmid (HJD) 2453984.397 ± 0.009 ω (deg) 0a able light curve in this fashion. ecc 0a Based on the above, we conclude that ADS 16402B is 1 K1 (m s− ) 60.3 ± 2.1 orbited by a transiting object that gives rise to the ob- arel (AU) 0.0551 ± 0.0015 served photometric and radial velocity variations. From m sin i(MJ) 0.53 ± 0.04 the 60ms−1 amplitude of the radial velocity orbit (§ 4), iP (deg) 85.9 ± 0.8 MP (MJ) 0.53 ± 0.04 we conclude that the object must have mass about 0.53 0.11 RP (RJ) 1.36±0.09 MJ and hence must be an extrasolar planet. a Fixed in the orbital fit. 6. THE TRANSITING PLANET HAT-P-1b AND ITS sectors (see, e.g., Santos et al. 2002; Torres et al. 2005). PROPERTIES We have searched for such line bisector variations in our In order to determine more precisely the physical char- Subaru spectra and find no evidence for them. Compar- acteristics of HAT-P-1b and its orbital parameters, we ison between the relative radial velocities and bisector have carried out model fitting of the different sets of pho- spans of ADS 16402A and B is shown in Fig. 5. The tometric data shown in the four panels of Fig. 2. A simple results are based on the Subaru spectra only, because goodness-of-fit was determined against the predicted flux they were obtained on the same nights for both stars. values evaluated with the algorithm of Mandel & Agol Line bisectors for each spectrum were computed from (2002). We used the values for the stellar radius and the cross-correlation function averaged over all orders, mass (and their uncertainties) determined above, and co- which is representative of the average spectral line pro- efficients for stellar limb darkening for the Sloan z-filter file of the star. The cross-correlations were performed by Claret (2004). First we fitted the original discovery against a synthetic spectrum matching the properties of data from HATNet with the appropriate correction for the stars. The bisector spans were then computed as the flux of the binary star (unresolved on the HATNet ◦ the velocity difference between points chosen near the images): the RP = 1.42RJ and iP = 85.4, have large top and bottom of the bisectors (see Torres et al. 2005). uncertainties, as anticipated, and might be subject to ADS 16402A shows no variation in either the velocities large systematic errors. We then fitted to the recent or the bisectors spans. ADS 16402B, on the other hand, Konkoly Schmidt observations of a complete HAT-P-1b shows obvious velocity variations but does not show the transit in I-band: the values are RP =1.31±0.07RJ and ◦ ◦ concomitant change in the bisector spans that would be iP = 86.0±1.0, and a precise estimate of transit duration expected if it were the result of a blend with an eclips- of d =0.11671 days and center of transit Tc = 3984.3967. ing binary. If the light and velocity variations were the The latter value allows us to estimate accurately the Tc result of a blend with an eclipsing binary, the bisector of the transit we observed just an orbital period earlier spans would be expected to vary with a similar ampli- (Aug 31) at FLWO, where the beginning of the ingress is tude as the velocities (see, e.g., Queloz et al. 2001). That missing (see Fig. 2, third panel from the top). Combin- scenario is therefore ruled out by these observations. ing both z-band transit curves from FLWO with a known We have also modeled the light curve as a blend fol- Tc allows an independent fit: we get RP =1.36 ± 0.05RJ ◦ ◦ lowing Torres et al. (2004b), subject to all available ob- and iP = 85.9 ± 0.8. The fit is better because the qual- servational constraints. Specifically, we calculated model ity and amount of photometric data is higher, and the light curves for a three-star system consisting of a slightly limb darkening in z-band is less. The fits to all data evolved early G star primary plus a later-type secondary sets give us the same planet parameters within the error; and a tertiary which together comprise an eclipsing bi- henceforth we use our best fit with the FLWO z-band Large Radius Transiting Planet in a Stellar Binary 7 transits. 1.6 The result for the planetary radius is RP = ρ = 0.4 ρ = 0.7 0.11 ± 1.4 HAT-P-1 1.36 0.09RJ, where the above statistical errors from the 209458b ρ = 1.0

fit and systematic errors from the uncertainties of the TrES-2 OGLE 56b . 0 09 XO-1b 189733b stellar radius and mass (± . ), have been added in 1.2 OGLE 10b 0 07 ] OGLE 132b ρ = 1.33 quadrature. Due to the asymmetric error bars on the Jup OGLE 111b TrES-1 OGLE 113b [R stellar radius (see above) we have asymmetric error bars pl 1 ◦ ◦ R Jupiter on RP too. The orbital inclination is iP = 85.9 ± 0. 8. Saturn The radius, orbital inclination, and derived mass of the 0.8 149026b planet are included in Table 3. [ρ]: g cm-3 An even better light curve with e.g. the Hubble Space 0.6 Telescope, and a more detailed analysis of the parameters of the stellar pair could reduce that error considerably. 0 0.5 1 1.5 2 The results of this on-going work will be reported in a Mpl [MJUP] forthcoming paper, as soon as the multi-band observa- Fig. 6.— Mass-radius diagram depicting known transiting exoplanets, plus Saturn and Jupiter (for comparison). tions enable us to refine the binary system parameters. Sources for the mass and radius values (in this order) are listed after the name of the planet; HAT-P-1b: this work, 7. DISCUSSION HD 209458b: Laughlin et al. (2005), Knutson et al. (2006), HD 189733b: Bakos et al. (2006), OGLE-TR-111b: Pont et al. Fig. 6 shows the location of HAT-P-1b in a plot of (2004), OGLE-TR-10b: Santos et al. (2006), OGLE-TR-132b: planetary radius versus planetary mass for the now 11 Moutou et al. (2004), OGLE-TR-56b: Torres et al. (2004a), known transiting exoplanets, plus Jupiter and Saturn. OGLE-TR-113b: Bouchy et al. (2004); Konacki et al. (2004), It is immediately seen that HAT-P-1b is in an extreme TrES-1: Sozzetti et al. (2004), HD 149026: Sato et al. (2005), Charbonneau et al. (2006a), XO-1b: Holman et al. (2006), position. It has a radius apparently as large as or larger TrES-2: O’Donovan et al. (2006). than the largest known exoplanet, HD 209458b, and at nisms requiring unusual circumstances seem less likely. the same time it has a mass that is significantly smaller. Moreover, for HD 209458b, explanation (b) appears Thus, it has a lower mean density than any other known ruled out because Doppler data (Laughlin et al. 2005b; planet. Until now, HD 209458b was the only significantly Winn et al. 2005) and transit timing variations data anomalous transiting planet. Thus, Guillot (2005) noted (Miller-Ricci et al. 2005) have excluded companion plan- that all of the transiting planets known at that time, with ets capable of enacting it. the exception of HD 209458b, fall within a mass-radius In the case of HAT-P-1b, it remains possible that the relation derived from modeling the interiors of strongly system has a second planet capable of forcing its eccen- irradiated planets. All of the transiting planets detected tricity and leading to internal tidal heating. Our current subsequently, except HAT-P-1b, also fit the same rela- orbital solution (§ 4) suggests e = 0.09 ± 0.02. If con- tion. Among them, OGLE-TR-10b had been thought to firmed, such eccentricity would be a strong indication be similarly anomalous (Konacki et al. 2005), but recent for the presence of a third body (second planet), because evidence (Holman et al. 2005; Santos et al. 2006, M. Hol- otherwise tidal dissipation would have circularized HAT- man, private communication) points to a smaller plan- Q P-1b’s orbit in ∼ 225( 6 ) Myr. The dissipation factor, etary radius and normal density. Similarly, the XO-1b 10 discovery paper (McCullough et al. 2006) derives a larger Q, is poorly known for “hot Jupiters”, but is expected to lie between the values for Jupiter and solar-like stars radius than the value later refined by (Holman et al. 5 6 2006). (10 to 10 ). With forced orbital eccentricity, e = 0.09, HD 209458b has a radius some 10 to 20 percent larger HAT-P-1b will be subject to internal tidal heating due than expected, given its age, mass, and expected temper- to an energy dissipation rate ature including effects of irradiation. Possible solutions 2 5 2 dE 63 e n Rp GM∗ include: (a) a fraction of a percent of the stellar flux irra- = dt 4 Q a a diating HD 209458b could be transported to deep layers of order 100 bar, and thereby keep the planet’s radius where G is the gravitational constant, M⋆ is the mass puffed up (Showman & Guillot 2002); (b) there is a small of the star, a is the semimajor axis, n is the mean eccentricity of its orbit, forced by the presence of a so-far motion, and Q is the specific dissipation factor. With undetected planetary companion and leading to internal e = 0.09, the rate of energy dissipation in HAT-P-1b tidal heating (Bodenheimer, Laughlin, & Lin 2003), and would be 3.0 × 1026erg s−1, corresponding to energy dis- 26 −1 −1 (c) HD 209458b might be in a Cassini State—a dynami- sipation per unit Jupiter mass of 5.7 × 10 erg s MJ . cal state where its rotation axis lies close to the plane of This rate is almost identical to the dissipation rate 26 −1 −1 its orbit (Winn & Holman 2005), leading to large obliq- per unit Jupiter mass of 5.8 × 10 erg s MJ re- uity tides that cause heating at significant depth leading quired to achieve a 20% increase in RP for HD 209458b to the planet’s large radius. A natural question about (no core model) and explain its anomalous low den- the first of these mechanisms is why all close-in “hot sity (Bodenheimer, Laughlin, & Lin 2003). Therefore we Jupiters” are not similarly bloated by that mechanism. conclude that if the orbital eccentricity is confirmed, the A problem with the last two mechanisms is that both anomalous RP and low density of HAT-P-1b could be are ad hoc; they require unusual circumstances for which easily explained by tidal heating. The perturbing sec- there is no independent evidence. ond planet would be detectable by both transit timing With the discovery of another planet, HAT-P-1b, variations and Doppler velocities. which has a radius similar to and probably even larger For now, however, all that can be said is that bloated than HD 209458b, and an even smaller density, mecha- planets are not that unusual, and, since a second perturb- 8 Bakos et al. ing planet is lacking in HD 209458b’s case, perhaps some and K. Peek for help collecting data at Subaru and Keck. other explanation should be sought that can explain such The Keck Observatory was made possible by the gener- planets more generally. ous financial support of the W. M. Keck Foundation. D. A. F is a Cottrell Science Scholar of Research Cor- poration. We acknowledge support from NASA grant Operation of the HATNet project is funded in part NNG05G164G to DAF. We would like to thank Carl Ak- by NASA grant NNG04GN74G. Support for program ´ erlof and the de-commissioned ROTSE-I project the gen- number HST-HF-01170.01-A to G. A. B. was provided erous loan of some of the lenses and CCDs that we use by NASA through a Hubble Fellowship grant from the for operating HATNet. Space Telescope Science Institute, which is operated by We owe special thanks to Emilio Falco, Dan Fabricant, the Association of Universities for Research in Astron- James Moran and Antony Schinckel for their help in es- omy, Incorporated, under NASA contract NAS5-26555. tablishing and operating the HATNet stations at FLWO G. K. wishes to thank support from Hungarian Scientific and SMA. G. A.´ B. wishes to thank the support given Research Foundation (OTKA) grant K-60750. We ac- by telescope operators Mike Calkins and Perry Berlind knowledge partial support from the Kepler Mission under in the operation of the FLWO HATNet station. We also NASA Cooperative Agreement NCC2-1390 (D. W. L., thank Gergely Gálfi for useful discussions. PI). G. T. acknowledges partial support from NASA Ori- This publication made use of the VizieR interac- gins grant NNG04LG89G. We thank Akito Tajitsu for his tive catalogue (Ochsenbein, Bauer, & Marcout 2000) at expertise and support of the Subaru HDS observations. CDS, Strasbourg, and the 2MASS catalogue. We thank graduate students M. Ohmiya, S. Robinson

REFERENCES Aitken, R. G., & Doolittle, E. 1932, “New general catalogue of Miller-Ricci, E., et al. 2005, Bulletin of the American Astronomical double stars within 120 of the North pole”, Carnegie Institution Society, 37, 1339 of Washington. Moutou, C., Pont, F., Bouchy, F., & Mayor, M. 2004, A&A, 424, Alonso, R., et al. 2004, ApJ, 613, L153 L31 Bakos, G. A.,´ Lázár, J., Papp, I., Sári, P., & Green, E. M. 2002, Noguchi, K. et al. 2002, PASJ, 54, 855 PASP, 114, 974 Ochsenbein, F., Bauer, P., & Marcout, J. 2000, A&AS, 143, 23 Bakos, G., Noyes, R. W., Kovács, G., Stanek, K. Z., Sasselov, D. D., O’Donovan, F. T. et al. 2006, accepted for publication in ApJL, & Domsa, I. 2004, PASP, 116, 266 astro-ph/0609335 Bakos, G. A.,´ Pál, A., Latham, D. W., Noyes, R. W., & Stefanik, Pont, F., Bouchy, F., Queloz, D., Santos, N. C., Melo, C., Mayor, R. P. 2006, ApJ, 641, L57 M., & Udry, S. 2004, A&A, 426, L15 Bodenheimer, P., Laughlin, G., & Lin, D. N. C. 2003, ApJ, 592, Queloz, D., Henry, G. W., Sivan, J. P., Baliunas, S. L., Beuzit, J. 555 L., Donahue, R. A., Mayor, M., Naef, D., Perrier, C., & Udry, S. Bouchy, F., Pont, F., Santos, N. C., Melo, C., Mayor, M., Queloz, 2001, A&A, 379, 279 D., & Udry, S. 2004, A&A, 421, L13 Santos, N. C., Mayor, M., Naef, D., Pepe, F., Queloz, D., Udry, Bouchy, F., et al. 2005, A&A, 444, L15 S.,Burnet, M., Clausen, J. V., Helt, B. E., Olsen, E. H., & Butler, R. P., Marcy, G. W., Williams, E., McCarthy, C., Dosanjh, Pritchard, J. D. 2002, A&A, 392, 215 P., Vogt, S. S. 1996, PASP, 108, 500 Santos, N. C., et al. 2006, A&A, 450, 825 Charbonneau, D., Brown, T. M., Latham, D. W., & Mayor, Sato, B., et al. 2005, ApJ, 633, 465 M. 2000, ApJ, 529, L45 Showman, A. P., & Guillot, T. 2002, A&A, 385, 166 Charbonneau, D., et al. 2006, ApJ, 636, 445 Skrutskie, M. F., et al. 2006, AJ, 131, 1163 Charbonneau, D., Brown, T. M., Burrows, A., & Laughlin, G. 2006, Sozzetti, A., et al. 2004, ApJ, 616, L167 Review Chapter for Protostars and Planets V, astro-ph/0603376 Torres, G., Konacki, M., Sasselov, D. D., & Jha, S. 2004, ApJ, 609, Fischer, D. A., Laughlin, G, et al. 2004, ApJ, 620, 481 1071 Guillot, T. 2005, Annual Review of Earth and Planetary Sciences, Torres, G., Konacki, M., Sasselov, D. D., & Jha, S. 2004, ApJ, 614, 33, 493 979 Halbwachs, J. L. 1986, A&AS, 66, 131 Torres, G., Konacki, M., Sasselov, D. D., & Jha, S. 2005, ApJ, 619, Herschel, J. F. W. 1831, MmRAS, 4, 331 558 Holman, M. et al. 2005, submitted to ApJ, astro-ph/0506569 Udalski, A., et al. 2002, Acta Astronomica, 52, 1 Holman, M. et al. 2006, accepted for publication in ApJ, Udalski, A., Zebrun, K., Szymanski, M., Kubiak, M., Soszynski, astro-ph/0607571 I., Szewczyk, O., Wyrzykowski, L., & Pietrzynski, G. 2002, Acta Konacki, M., et al. 2004, ApJ, 609, L37 Astronomica, 52, 115 Konacki, M., Torres, G., Sasselov, D. D., & Jha, S. 2005, ApJ, 624, Udalski, A., Szewczyk, O., Zebrun, K., Pietrzynski, G., Szymanski, 372 M., Kubiak, M., Soszynski, I., & Wyrzykowski, L. 2002, Acta Kovács, G., Zucker, S., & Mazeh, T. 2002, A&A, 391, 369 Astronomica, 52, 317 Kovács, G., Bakos, G., & Noyes, R. W. 2005, MNRAS, 356, 557 Udalski, A., Pietrzynski, G., Szymanski, M., Kubiak, M., Zebrun, Knutson, H. et al. 2006, submitted to ApJ, astro-ph/0603542 K., Soszynski, I., Szewczyk, O., & Wyrzykowski, L. 2003, Acta Henry, G. W., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2000, Astronomica, 53, 1 ApJ, 529, L41 Udalski, A., Szymanski, M. K., Kubiak, M., Pietrzynski, G., Kambe, E. et al. 2002, PASJ, 54, 873 Soszynski, I., Zebrun, K., Szewczyk, O., & Wyrzykowski, L. 2004, Johnson, J. et al. 2006, ApJin press, astro-ph/0604348 Acta Astronomica, 54, 313 Landolt, A. U. 1992, AJ, 104, 340 Valenti, J. A., & Piskunov, N. 1996, A&AS, 118, 595 Latham, D. W. 1992, ASP Conf. Ser. 32: IAU Colloq. 135: Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141 Complementary Approaches to Double and Multiple Star Vogt, S. S., et al. 1994, Proc. SPIE, 2198, 362 Research, 32, 110 Winn, J. N., & Holman, M. J. 2005, ApJ, 628, L159 Laughlin, G. et al. 2005, ApJ, 621, 1072 Winn, J. N., et al. 2005, ApJ, 631, 1215 Laughlin, G., Marcy, G. W., Vogt, S. S., Fischer, D. A., & Butler, Wright, J. T. 2005, PASP, 117, 657 R. P. 2005, ApJ, 629, L121 Mandel, K., & Agol, E. 2002, ApJ, 580, L171 McCullough, P. R., et al. 2006, ApJ, accepted, astro-ph/0605414