Arxiv:0705.0126V1 [Astro-Ph] 1 May 2007 Ahntnd,54 Ra Rnhr.N,Wsigo DC, Washington NW, Rd

Home , HATNet Project

arXiv:0705.0126v1 [astro-ph] 1 May 2007 ahntnD,54 ra rnhR.N,Wsigo DC, Washington NW, Rd. Branch 20015-1305 Broad USA 5241 DC, Washington h aionaIsiueo ehooy ektm a eni been has time NASA. Keck by Technology. Califor granted of of Institute University California the by the operated is which Observatory, 219 Box O. P. rni hi aetsasadtu il ausfrtheir for values yield radius thus and and mass stars parent their transit A970 USA 94720, CA .Mazeh T. -58Bdps,Hungary. Budapest, H-1518 e tet abig,M 23,UA [email protected]. USA; 02138, MA Cambridge, Street, den as(aoe l 05,ipyn hti a large a has it that its for implying radius 2005), ( small al. core a heavy et has there (Sato 149026b although Marley, HD mass therein), Fortney, references exceptions: (e.g. and are theory planets 2006, from giant Barnes expects & gas one irradiated what for approximately fit jority oi 0 02,Pn oiee Italy Torinese, Pino 10025, 20, torio axis (semi-major close-in some are . there that expected with consistent orbits motion. have Keplerian TEPs exo- All circular d transiting P=4.46 2007). with al. HAT-P-1b density et is (Bakos lowest far suggesting so detected and perhaps (TEP) planet period masses, heating. longest internal their extra The of source for unknown presently radii some unexpect- large have WASP-1) edly HAT-P-1b, 209458b, (HD eral osdrbylre hnayo h 4tastn planets transiting 14 the of any than larger considerably nvriy a rnic,C 43,USA 94132, CA Francisco, San University, gary 1.3 10 12 † 9 8 5 6 3 2 7 4 odt 4etaoa lnt aebe on which found been have planets extrasolar 14 date To rmeitn ailvlct R)dt,i ih be might it data, (RV) velocity radial existing From 1 rf eso a ,20 ESO umte oApJ to Submitted VERSION L – using 2007 typeset 1, Preprint May version Draft 0 M eateto ersra ants,Crei nttt o Institute Carnegie Magnetism, Terrestrial of Department eateto srnm,Uiest fClfri,Berkel California, of University Astronomy, of Department NF-Osraoi srnmc iTrn,Srd Osserva- Strada Torino, di Astronomico Osservatorio - INAF 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 ieOsraoy e vvUiest,TlAi,Ire 699 Israel Aviv, Tel University, Aviv Tel Observatory, Wise avr-mtsna etrfrAtohsc CA,6 Gar 60 (CfA), Astrophysics for Center Harvard-Smithsonian A--b UE-ASV LNTI NECNRCOBTTRAN ORBIT ECCENTRIC AN IN PLANET SUPER-MASSIVE A HAT-P-2b: . ae npr nosrain banda h .M Keck M. W. the at obtained observations on part in Based ugra srnmclAscain 41Budapest, 1461 Association, Astronomical Hungarian xrslrPaesEccoei:http://exoplanet.eu Encyclopedia: Planets Extrasolar G. 7A,or AU, 07 J n ai rm0.7 from radii and , .Bakos A. ´ rgt( bright rmtetastlgtcrew eemn htterdu ft of radius the that determine we curve light transit the From aisaei codwt hoiso tutr fmsiegi insolati massive of variation of headings: 9-fold structure Subject a of causes theories eccentricity with high The accord in are radius A--bhsams bu ie h vrg aso previous of mass average the times of 9 density about a mass a has HAT-P-2b 4 erpr h icvr fHTP2,amsie( massive a HAT-P-2b, of discovery the report We .Shporer A. , ∼ 70 V P 12 M 1,2 8 = assrnefo 0.3 from range Masses . 1. . ⊕ .Kov G. , ρ aglne l 05,adsev- and 2005), al. et Laughlin ; A p INTRODUCTION . T 0d)gatpaeswt masses with planets giant ) d 10 )F trH 456 iha ria eido .3dy n an and days 5.63 of period orbital an with 147506, HD star F8 7) E 6 = 4 tl mltajv 12/14/05 v. emulateapj style X .P Butler P. R. , R . tr:idvda:H 456–paeaysses individu systems: planetary – 147506 HD individual: stars: cm g 6 J acs ´ .W Marcy W. G. oaot1.4 about to 3 .Torres G. , − rf eso a ,20 ESO umte oApJ to Submitted VERSION – 2007 1, May version Draft 3 iia ota frcypaeslk h at.Neverthele Earth. the like planets rocky of that to similar , 9 .P Stefanik P. R. , R 8 M .Sip B. , J 1 h ma- The . J .A Fischer A. D. , oabout to i and nia part n ocz ˝ ABSTRACT edu. 32, ey, 78 7,1 - f .L J. , 1 .M Fern M. J. , 6 ehrti eutdi ih uv ihexceptional with of curve number unprecedented light an a days), in (570 time-coverage resulted this gether contributing HAT- WHAT, chance of and and (G149) By FLWO) field at joint (HAT-7 another data-sets. Net yet in these is de- candidate with the independently confirmed was and transit tected The respectively. points, bevtr,Ire;Sprre l 06,fra extended yielded an for that Wise- 2006), (WHAT, al. period et telescope Shporer HAT Wise Israel; Observatory, Kea, the Mauna by atop and site Hawaii, (SMA) telescope Array HAT-9 Submillimeter HATNet’s that the by (G192) at observed field jointly another been with has corner was overlapping star the the Fortunately in marginal. was of cadence) consisting 5.5-min curve light the a The of of (SAO). detection (FLWO) Observatory Observatory Astrophysical Smithsonian Whipple the Lawrence Fred by the obtained at data the in in candidate G193 HATNet’s transit field a labeled as internally identified our initially was It 1997). n ipro aalx7 parallax Hipparcos and sas h ogs eidTP n h rtoet ex- to one first the and orbit. TEP, eccentric period such This highly longest hibit first radius. the and the mass also of its detection is of the determination our here par- and their report TEP, “super-massive” We transiting Thus found stars. ent be 10%. sometimes about should is planets star its transiting ido . n eimjrai f009A.A this At AU. 0.069 of pe- axis orbital b semi-major the and and HIP14810 distance mass, d orbital is similar 6.7 example with of Another 2007) riod al. 1997). et (Wright al. et (Butler etdfo Vvrain,adhsamnmmms of mass minimum a has is and M variations, 2000), RV from al. tected et Chabrier e.g. et eo h etru unn hehl ( threshold burning Deuterium ob- the only below considering example, jects well-known A known. now .W Latham W. D. , M az ´ 13 D170 sa 8sa ihvsa antd 8.7 magnitude visual with star F8 an is 147506 HD p p ar ´ no h lntbtenpr-adapastron. and peri- between planet the of on 2.1. sin http://www.hatnet.hu ∼ 8 = 5 20adtoa aapit,rsetvl.Alto- respectively. points, data additional 2200 .Papp I. , i n lnt opsdo ueHadHe. and H pure of composed planets ant =3.9 eeto ftetasti h ANtdata HATNet the in transit the of Detection . andez ´ 17 13 epae is planet he 2. ykontastn xpaes and exoplanets, transiting ly-known ± Bkse l 02 04 A- telescope HAT-6 2004) 2002, al. et (Bakos M BEVTOSADANALYSIS AND OBSERVATIONS ∼ 0 5 1 J . .S P. , .Sozzetti A. , mgtastwt .3dpro in period d 5.63 a with transit mmag 5 72 1 n risol .4 Ufo t star its from AU 0.046 only orbits and .W Noyes W. R. , M priori a ∼ ari ´ Jup 70and 6700 R IIGABIH STAR BRIGHT A SITING 5 . lnttastn the transiting planet ) l HAT-P-2b al: 39 p cetiiyof eccentricity 1 = ± rbblt fsc planet a such of probability ∼ 1,5 0 τ 00dt-ons(iha (with data-points 7000 . . .P A. , 18 8ms(ermne al. et (Perryman mas 88 1 ∼ s t asand mass its ss, .D Sasselov D. D. , o ,wihwsde- was which b, Boo 90adtoa data- additional 3900 ± 0 al ´ . 16 7,1 e 0 = .Johnson J. , R Jup . 5. . ∼ ∼ 13 1 , 6200 † M J 8 , , 2 Bakos et al. 8.05 8.06 8.07 8.08

I mag 8.09 8.1 HATNet, all 8.11 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 Phase 8.07 8.075 8.08 I mag 8.085 HATNet, all HATNet, binned 8.09 -0.1 -0.05 0 0.05 0.1 Phase 8.07

8.075

8.08

z mag (inst) 8.085 FLWO 1.2m Fit 8.09 12.65 12.7 12.75 12.8 12.85 12.9 12.95 13 13.05 HJD - 2454200.0 Fig. 1.— The upper panel shows the unbinned HATNet and WHAT joint light curve with 26400 data-points, phased with the P = 5.63341 d period. The 5 mmag deep transit is detected with a signal-to-noise of 26. The middle panel shows the same HATNet and WHAT data zooming in on the transit and binned with a φ = 0.0005 bin-size. The lower panel displays the Sloan z-band photometry taken with the FLWO 1.2 m telescope. Overplotted is our best fit obtained with the Mandel & Agol (2002) formalism. data-points (26461 measurements at 5.5-min cadence), The Sloan z-band light curve is shown in the lower panel and an rms of 5 mmag. It is noteworthy that the net- of Fig. 1. From the combined HATNet and KeplerCam work coverage by WHAT (longitude 35 E), HATNet at photometry, spanning a baseline of 839 days, we derive a FLWO (111 W) and HATNet at Hawaii (155 W) played period of 5.63341±0.00013 d and an epoch of mid-transit an important role in detecting such a long period and of Tmid = 2,454,212.8561 ± 0.0023 d (HJD). From the shallow transit. Data were reduced with a highly fine- FLWO 1.2 m data alone (and the analytic light curve fit tuned aperture photometry. We applied our external as described later), the length of transit is 0.177±0.002 d, parameter decorrelation (EPD) technique on the light and the length of ingress is 0.012 ± 0.002 d. curves, whereby deviations from the median were cross- correlated with a number of “external parameters”, such as the X and Y sub-pixel position, hour-angle, and zenith 2.2. Early spectroscopy follow-up distance. We have also applied the Trend Filtering Algo- Initial follow-up observations were made with the CfA rithm (TFA; Kovács, Bakos, & Noyes 2005) along with Digital Speedometer (DS; Latham 1992) in order to char- the Box Least Squares (BLS; Kovács, Zucker, & Mazeh acterize the host star and to reject obvious astrophysi- 2002) transit-search algorithm in our analysis (these two cal false-positive scenarios that mimic planetary tran- being combined in signal-reconstruction mode). The de- sits. These observations yielded values of Teff = 6250 K, tection of this relatively shallow transit is a good demon- log g = 4.0 and v sin i = 22 kms−1, corresponding to a stration of the strengths of TFA. The upper panel of moderately-rotating main sequence F star. The radial Fig. 1 shows the unbinned light curve with all 26400 velocity (RV) measurements showed an rms residual of data points, whereas the middle panel displays the tran- ∼0.82 kms−1, slightly larger than the nominal DS pre- sit binned to 1/2000 of the period (4 minutes). We note cision for a star with this rotation, and suggested that that due to the large amount of data, the binned light they may be variable. With a few dozen additional DS curve is of similar precision as a single-transit observation observations, it was found that the RV appeared peri- by a 1m-class telescope. odic with P ≈ 5.63 d, semi-amplitude ∼1 kms−1, and After several failed attempts (due to bad weather and phasing in agreement with predictions from the HAT- instrumental failure) to carry out high-precision photo- Net+WHAT light curve. This gave strong evidence that metric follow-up observations from FLWO, Wise Ob- there really was an RV signal resulting from Keplerian servatory, Konkoly Observatory, and the Clay Center motion, although the precision was insufficient to estab- (Boston), we finally succeeded in observing a full tran- lish the orbit with confidence. Altogether we collected sit using the KeplerCam detector on the FLWO 1.2 m 53 individual spectra spanning a time-base of more than telescope (see Holman et al. 2007) on UT 2007 April 22. a year (Table 1). A Super-Massive Transiting Planet 3

TABLE 1 TABLE 2 Radial Velocities for HD 147506. Summary of stellar parameters for HD 147506.

BJD −2,400,000 RVa Uncert. Observatoryb Parameter Value Source (days) (m s−1) (ms−1) Teff (K) 6290 ± 110 SME log g 4.22 ± 0.14 SME 53981.7785 -540.6 8.6 Keck 1 53982.8728 -851.5 8.8 Keck v sin i(km s− ) 19.8 ± 1.6 SME 53983.8160 -47.8 9.1 Keck [Fe/H](dex) +0.12 ± 0.08 SME −1 53984.8962 294.4 8.8 Keck Vmac (km s ) 4.51 ± 0.4 SME −1 54023.6940 169.2 10.2 Keck Vmic (km s ) 0.85 SME 54186.9962 132.3 5.5 Keck Distance (pc) 135 ± 16 HIP 54187.1021 114.8 5.5 Keck L∗ (L⊙) 4.25 ± 0.66 — 54187.1578 139.8 5.4 Keck Mass (M⊙) 1.35 ± 0.10 See text 54188.0148 177.9 5.2 Keck Radius (R⊙) 1.80 ± 0.25 See text 54188.1575 206.8 5.5 Keck ±1.4 Age (Gyr) 2.7 0.6 Y2 isochrones, MV fixed 54189.0083 82.2 5.6 Keck 54189.0868 79.9 6.0 Keck 54189.1556 33.7 6.0 Keck 2006), the SME temperature, and a bolometric correc- 54168.9666 -518.7 24.7 Lick tion of BCV = −0.011 ± 0.011 mag (Flower 1996), appli- 54169.9505 127.0 23.9 Lick cation of the Stefan-Boltzmann law yields a stellar radius 54170.8605 162.9 34.1 Lick 54171.0350 231.4 35.1 Lick of R∗ =1.84 ± 0.24 R⊙. As an independent method to estimate the radius we a The RVs include the barycentric correction. used the surface-brightness relations from Kervella et al. b Only the Keck and Lick data-points are shown here. Con- (2004) that give the (limb-darkened) angular diameter sult the electronic edition for a full data-set that includes the of a star as a function of its V − K color index and CfA DS measurements. K apparent magnitude. We neglected reddening, since 2.3. High-precision spectroscopy follow-up HD 147506 is a nearby star at a moderately high galac- tic latitude. After conversion of the Ks-band magnitude In order to confirm or refute the planetary nature of from 2MASS (Skrutskie et al. 2006) to the homogenized the transiting object, we pursued follow-up observations Bessell & Brett system (Carpenter 2001) for this applica- with the HIRES instrument (Vogt et al. 1994) on the tion, we obtained the angular diameter which, combined W. M. Keck telescope and with the Hamilton Echelle with the Hipparcos parallax, yields a stellar radius of spectrograph at the Lick Observatory (Vogt 1987). The R∗ = 1.75 ± 0.25R⊙. The error budget here is domi- ′′ spectrometer slit used at Keck is 0. 86, yielding a re- nated by the uncertainty in the parallax. solving power of about 55,000 with a spectral coverage A more sophisticated approach to determine the stel- between about 3200 and 8800 A.˚ The Hamilton Echelle lar parameters uses stellar evolution models along with spectrograph at Lick has a similar resolution of about the observational constraints from spectroscopy. For this 50,000. These spectra were used to i) more fully char- we used the models by Yi et al. (2001) and Demarque et acterize the stellar properties of the system, ii) to ob- al. (2004), and explored a wide range of ages to find all tain a radial velocity orbit, and to iii) check for spec- models consistent with Teff , log g, and [Fe/H] within the tral line bisector variations that may be indicative of a observational errors. The stellar mass and radius we ob- blend. We gathered 13 spectra at Keck (plus an iodine- 0.17 0.36 tained are M∗ =1.29 ±0.12 M⊙ and R∗ =1.46 ±0.27 R⊙. free template) spanning 207 days, and 4 spectra at Lick The surface gravity is a sensitive measure of the degree of (plus template) taken on 4 consecutive nights. The ra- evolution of the star, as is luminosity, and therefore has dial velocities measured from these spectra are shown in a very strong influence on the radius. However, log g is Table 1, along with those from the CfA DS. a notoriously difficult quantity to measure spectroscop- 3. ically and is often strongly correlated with other spec- STELLAR PARAMETERS troscopic parameters. As a test, we repeated the com- A spectral synthesis modeling of the iodine-free Keck parison with the models using the absolute visual mag- template spectrum was carried out using the SME soft- nitude (MV ) instead of log g as a proxy for luminosity, ware (Valenti & Piskunov 1996), with the wavelength adopting the value resulting from V and the parallax ranges and atomic line data described by Valenti & Fis- (MV = 3.05 ± 0.26 mag). This resulted in a somewhat 0.08 cher (2005). Results are shown in Table 2. The val- larger mass and radius for the star of M∗ =1.42±0.12M⊙ ues obtained for the effective temperature (T ), surface 0.31 +1.4 eff and R∗ =1.85±0.28 R⊙, and a best-fit age of 2.7−0.6 Gyr. gravity (log g), and projected rotational velocity (v sin i) Given the better agreement in the radius with the inde- are consistent with those found from the CfA DS spectra. pendent estimate from the surface brightness relation, As a check on Teff , we collected all available photometry and our concerns about log g, we have chosen to rely in for HD 147506 in the Johnson, Cousins, 2MASS, and Ty- the following on the results based on MV . cho systems, and applied a number of color-temperature As an additional estimate of the stellar mass, we in- calibrations (Ram´ırez & Meléndez 2005; Masana et al. terpolated the star’s Hipparcos parallax–based absolute 2006; Casagrande et al. 2006) using 7 different color magnitude and B − V color, and SME–derived [Fe/H] indices. These resulted in an average temperature of onto the Padova14 stellar model grids (Girardi et al. ∼6400±100 K, somewhat higher than the spectroscopic 2002). We used a Bayesian methodology similar to Pont value but consistent within the errors. & Eyer (2004) to generate a complete posterior probabil- Based on the Hipparcos parallax (π =7.39±0.88mas), the apparent magnitude V = 8.71 ± 0.01 (Droege et al. 14 http://pleiadi.pd.astro.it/ 4 Bakos et al. ity distribution. Taking the median of this distribution +0.09 we derive a M∗ = 1.37−0.12M⊙, where the upper and lower bounds correspond to the 95% confidence intervals. We also used a different set of stellar evolution codes (Sienkiewicz, Paczyński, & Ratcliff 1988; Cody & Sas- selov 2005) and repeated the exercise, finding R∗ = 1.9±0.20R⊙, and M∗ =1.32±0.10M⊙, which are similar to the values from the other models. As a compromise we adopt for HD 147506 the values M∗ =1.35 ± 0.10M⊙ and R∗ =1.80 ± 0.25R⊙ (see Table 2). 3.1. Stellar jitter Stars with significant rotation are known to exhibit excess scatter (“jitter”) in their radial velocities (e.g., Wright 2005, and references therein), due to enhanced chromospheric activity and the associated surface inho- mogeneities (spottedness). This jitter is in addition to the internal errors in the measured velocities, and could potentially be significant in our case. We note that after pre-whitening the light curve with the transit com- Fig. 2.— The upper panel shows the RV measurements phased ponent, we found no significant sinusoidal signal above with the period of P = 5.63341 d. The zero-point in phase cor- 0.3 mmag amplitude. From this we conclude that there responds to the epoch of mid-transit. Large filled circles indicate is no very significant spot activity on the star (in the ob- Keck and Lick points. Small open circles denote CfA DS data (not used for the fit). Overlaid is the fit that was based only on the served 500 day window). In order to estimate the level Keck and Lick data assuming 40 m s−1 stellar jitter. The mid- of chromospheric activity in the star, we have derived dle panel shows the residual of the fit. The lower panel exhibits an activity index from the Ca I H and K lines in our the line bisector spans on the same scale as the upper panel. No Keck spectra of log R′ = −4.72 ± 0.05. For this value variation in the line bisectors is seen concomitant with that in the HK RVs, essentially confirming the planetary nature of the transiting the calibration by Wright (2005) predicts velocity jitter object. ranging from 8 to 16 ms−1. An earlier calibration by Saar, Butler, & Marcy (1998), parametrized in terms of and transit epoch fixed at the photometric values given the projected rotational velocity, predicts a jitter level earlier. The parameters adjusted are the velocity semi- of up to 50 ms−1 for our measured v sin i of 20 kms−1. amplitude K, the eccentricity e, the longitude of perias- A different calibration by the same authors in terms of tron ω, the center-of-mass velocity for the Keck relative ′ −1 RHK gives 20 ms . An additional way to estimate the velocities γ, and offsets ∆vKL between Keck and Lick jitter of is to compare HD 147506 to stars of the Lick pro- and ∆vKC between Keck and CfA DS. The fitted pa- gram that have similar properties (0.4 < B − V < 0.5, rameters were found to be fairly insensitive to the level v sin i> 15kms−1). There are four such stars (J. John- of jitter assumed. However, only for a jitter of ∼40ms−1 son, private communication), and their average jitter is (or ∼50 ms−1 when the CfA DS data are included) did 45ms−1. A more direct measure for the particular case the χ2 approach values expected from the number of de- of HD 147506 may be obtained from the multiple ex- grees of freedom. There are thus two possible conclu- posures we collected during a 3-night Keck run in 2007 sions: if we accept that HD 147506 has stellar jitter at March. Ignoring the small velocity variations due to or- the 40 ms−1 level, then a single-planet solution such as bital motion during any given night, the overall scat- ours adequately describes our observations. If, on the ter of these 8 exposures relative to the nightly means is other hand, the true jitter is much smaller (. 20 ms−1), ∼20 ms−1. This may be taken as an estimate of the jit- then the extra scatter requires further explanation (see ter on short timescales, although it could be somewhat below). Our adopted orbital parameters for the simplest larger over the entire span of our observations. Alto- single-planet Keplerian solution are based only on the gether, it is reasonable to expect the jitter to be at least more precise Keck and Lick data, and assume the jitter 10 ms−1, and possibly around 30–50ms−1 for this star. is40ms−1 (Table 3). The orbital fit is shown graphically 4. in the upper panel of Fig. 2. In this figure, the zero-point SPECTROSCOPIC ORBITAL SOLUTION of phase is chosen to occur at the epoch of mid-transit, We have three velocity data sets available for analysis: Tmid = 2,454,212.8561 (HJD). The most significant re- 13 relative radial velocity measurements from Keck, 4 sults are the large eccentricity (e = 0.507 ± 0.012), and from Lick, and 53 measurements from the CfA DS, which the large velocity semi-amplitude (K = 884 ± 59kms−1, are nominally on an absolute scale (Table 1). Given the indicating a very massive companion). As we show in potential effect of stellar jitter, we performed weighted the next section (§ 5), the companion is a planet, thus Keplerian orbital solutions for a range of jitter values we denote it hereafter as HAT-P-2b. of 5, 10, 20, 30, 40, and 50 ms−1. These jitter val- As a consistency check we also fitted the orbits by fix- ues were added in quadrature to all individual inter- ing only the period, and leaving the transit epoch as nal errors. We performed separate fits for the star or- a free parameter. We found that for most values of the bited by a single planet, both with and without the CfA stellar jitter the predicted time of transit as derived from DS measurements, since these have errors (∼600 ms−1) the RV fit was consistent with the photometric ephemeris significantly larger than Keck (5–9 m s−1) or Lick (25– within the uncertainties. However, we found that in these 35 ms−1). In all of these solutions we held the period fits the orbital parameters were much more sensitive to A Super-Massive Transiting Planet 5 the level of jitter and to whether or not the CfA DS data 3 were included. The eccentricity values ranged from 0.37 to 0.49, though they were always significantly non-zero 2.5 (at the 10-sigma level). Finally, we note that the first of our Keck observations

] 2

was obtained about 20 minutes before the “first contact” Jup [R start of a transit, according to the transit ephemeris dis- P = 0.4 ρ = 0.7 = 5.5 R 1.5 ρ ρ cussed above. If our derived photometric period and ρ = 1.0 transit epoch are in error by more than 2 sigma there HAT-P-2b ρ = 1.33 is a chance that the transit could have started 20 min- 1 utes earlier, in which case the measured velocity may be [ρ]: g cm-3 affected by the Rossiter-McLaughlin effect, which occurs 0.5 when a transiting body passes in front of a rotating star. 0.1 1 10 100 MP [MJUP] 4.1. Solutions involving two planets Fig. 3.— The mass–radius diagram of known TEPs (from www.exoplanet.eu and references therein), Jupiter and Saturn If we assume that the true stellar jitter is small, then (large filled circles), and low-mass stars from Beatty et al. (2007). the excess scatter in the RV fit could be explained by a HAT-P-2b is an intermediate mass object that is still in the plan- third body in the system, i.e., a hypothetical HAT-P-2c. etary regime (well below 13MJ). Overlaid are equidensity lines (labeled), Baraffe et al. (1998) (stellar) and Baraffe et al. (2003) In addition, such a body could provide a natural dynam- (zero insolation planetary) isochrones for ages of 0.5 Gyr (upper, ical explanation for the large eccentricity of HAT-P-2b at dotted line) and 5 Gyr (lower dashed-dotted line), respectively. this relatively short period orbit. Preliminary two-planet broadening of the star as determined from the SME anal- orbital fits using all the data yielded solutions only sig- ysis. As shown in Fig. 2, while the measured velocities nificant at the 2-sigma level, not compelling enough to exhibit significant variation as a function of phase (up- consider as evidence for such a configuration. Additional per panel), the bisector spans are essentially constant RV measurements are needed to firmly establish or refute within the errors (lower panel). Therefore, this analysis the existence of HAT-P-2c. rules out a blend scenario, and confirms that the orbiting We also exploited the fact that the HATNet light curve body is indeed a planet. has a unique time coverage and precision, and searched 6. for signs of a second transit that might be due to an- PLANETARY PARAMETERS other orbiting body around the host star. Successive For a precise determination of the physical properties box-prewhitening based on the BLS spectrum and as- of HAT-P-2b we have modeled the FLWO 1.2 m Sloan z- suming trapezoidal-shape transits revealed no secondary band photometric data shown in Fig. 1 using the formal- transit deeper than the 0.1% level and period . 10 days. ism of Mandel & Agol (2002). Values for the stellar ra- 5. dius and mass (and their uncertainties) were adopted as EXCLUDING BLEND SCENARIOS determined above, along with quadratic limb-darkening As an initial test to explore the possibility that the coefficients taken from Claret (2004). The fitted parame- photometric signal we detect is a false positive (blend) ters were the planet-to-star radius ratio (Rp/R∗), the im- due to contamination from an unresolved eclipsing bi- pact parameter in units of the stellar radius (h/R∗), and nary, we modeled the light curve assuming there are three the ratio arel/R∗, where arel is the semi-major axis of the coeval stars in the system, as described by Torres et al. relative orbit. Strictly speaking this modeling assumes (2004). We were indeed able to reproduce the observed a circular orbit, whereas HAT-P-2b shows a substantial light curve with a configuration in which the brighter ob- eccentricity. However, only the arel/R∗ parameter is sig- ject is accompanied by a slightly smaller F star which is nificantly affected. This quantity has no relevance for our in turn being eclipsed by a late-type M dwarf. However, analysis (we interpret it here only as a free parameter for the predicted relative brightness of the two brighter ob- optimizing the fit), so the difference is of no consequence. jects at optical wavelengths would be ∼ 0.58, and this On the other hand the ratio of the radii, which is the pa- would have been easily detected in our spectra. This rameter in the fit that determines the size of the planet, is configuration can thus be ruled out. unaffected by the assumption of orbit circularity because The reality of the velocity variations was tested by it depends essentially only on the depth of the transit. carefully examining the spectral line bisectors of the Similarly, the impact parameter is little affected because star in our more numerous Keck spectra. If the veloc- the orbit is nearly edge-on. ity changes measured are due only to distortions in the The result for the radius ratio is Rp/R∗ = 0.0679 ± line profiles arising from contamination of the spectrum 0.0001 corresponding to a planetary radius of Rp = by the presence of a binary with a period of 5.63 days, 1.18 ± 0.16 RJ, where the above statistical errors from we would expect the bisector spans (which measure line the fit and systematic errors from the uncertainties of asymmetry) to vary with this period and with an am- the stellar radius have been added in quadrature. Our plitude similar to the velocities (see, e.g., Queloz et al. fits show that the orbital inclination is consistent with a ◦ 2001; Torres et al. 2005). The bisector spans were com- central transit, i.e., ip = 90.0±1 . The planet mass from puted from the cross-correlation function averaged over spectroscopy is then 8.17±0.72 MJ. These properties are 15 spectral orders blueward of 5000 A˚ and unaffected by summarized in Table 3. the iodine lines, which is representative of the average 7. spectral line profile of the star. The cross correlations DISCUSSION were performed against a synthetic spectrum matching In comparison with the other 14 previously known the effective temperature, surface gravity, and rotational transiting exoplanets, HAT-P-2b is quite remarkable 6 Bakos et al.

TABLE 3 HAT-P-2b Orbital fit and planetary parameters for the 100 HAT-P-2 system.

Parameter Value ] -2 Period (d)a 5.63341 ± 0.00013

a [m s mid ± P

T (HJD) 2,454,212.8561 0.0023 g Transit duration (day) 0.177 ± 0.002 Ingress duration (day) 0.012 ± 0.002 Stellar jitter (m s−1)b 40 10 γ (m s−1)c −237 ± 18 K (m s−1) 884 ± 59 ω (deg) 184.6 ± 3.1 1 2 3 4 5 6 e 0.507 ± 0.012 P [days] Tperi (HJD) 2,454,213.444 ± 0.044 −1 ∆vKL (m s ) −16 ± 29 Fig. 4.— The surface gravity of TEPs as a function of orbital −1 d period. Data taken from Southworth, Wheatley, & Sams (2007) ∆vKC (km s ) 19.861 ± 0.087 −9 with the exception of HAT-P-2b. This object is clearly not obeying f(M) (M⊙) (258 ± 56) × 10 2/3 the suspected correlation of gp and P for “Jupiter-mass” objects. Mp sin i (MJ) 6.66 ± 0.48([M∗ + Mp]/M⊙]) 6 a∗ sin i (km) (0.0590 ± 0.0042) × 10 arel (AU) 0.0685 ± 0.0017 HAT-P-2b is also unique in having an orbit with remark- ip(deg) 90.0 ± 1 ably high eccentricity. Aside from the issue of how a high ± Mp(MJ) 8.17 0.72 eccentricity might be maintained over 2 to 3 Gyr, there Rp(RJ) 1.18 ± 0.16 −3 is the question of how such an eccentricity was created ρp(g cm ) 6.61 ± 2.71 −2 gp(ms ) 149 ± 13 in the first place. One possible explanation could be a that the planet was scattered inward from a larger or- Fixed in the orbital fit. b Adopted (see text). bit, acquiring a high eccentricity in the process (Ford & c The γ velocity is not in an absolute reference frame. Rasio 2007). If so, then the scattering event might have d The offset between Keck and CfA DS is given for reference caused its new orbital plane to inclined relative to the from a fit that includes all data sets, but does not affect our plane of the original disk, and hence out of the equatorial solution. plane of the parent star. This angle between these two (Fig. 3, Fig. 4). Its mass of 8.17 ± 0.72MJ is ∼5 times planes should be readily measurable from the Rossiter- greater than any of these 14 other exoplanets. Its mean McLaughlin effect (Winn et al. 2005). Indeed, the star density ρ =6.6±2.7gcm−3 is ∼5 times that of the dens- HD 147506 is an ideal subject for studying this effect, be- est known exoplanet (OGLE-TR-113b, ρ =1.35 gcm−3) cause its rapid rotation should lead to a relatively large and indeed greater than that of Earth (ρ = 5.5 gcm−3) Rossiter signal. or other solar system rocky planets. Its surface gravity There are a number of other interesting issues related of 149 ± 13 ms−2 is 5 times that of any of the previously to the very high eccentricity of HAT-P-2b. During its known TEPs, and 20 times that of HAT-P-1b. 5.63-day orbit, the insolation reaching the planet’s sur- We may compare the mass and radius for HAT-P-2b face varies by a factor of 9. Assuming an albedo of 0.1 with evolutionary models, including irradiation, as re- (Rowe et al. 2006) and complete redistribution of insola- cently presented by Fortney, Marley, & Barnes (2006). tion energy over the surface of the planet, the equilibrium Given the inferred stellar luminosity (Table 2), and the temperature varies from about 2150 K at periastron to time-integral of the insolation over an entire period (tak- 1240 K at apastron. This would have a major influence ing into account the orbital parameters, notably e and on atmospheric dynamics and photochemistry. arel), the equivalent semi-major axis aeq for the same Because of its high orbital eccentricity, the planet does amount of irradiation if the central star were solar is not transit across the face of the star at a uniform rate: 0.031 AU. At that separation, Fortney, Marley, & Barnes given the value of ω = 184.6◦, the projected motion will (2006) find for a pure hydrogen/helium planet of mass change from start to end of each transit, and this will al- 8MJ and age of 4.5 Gyr a planetary radius about 1.17 RJ, ter the symmetry of the light curve. This effect may be that is, very close to the radius we infer. The radius is measurable by very high precision photometry. Similarly, about 10% larger at age 500 Myr, and about 10% smaller the Rossiter-McLaughlin effect should also show anoma- for a 4.5 Gyr old and 8MJ mass planet with a 10M⊕ core; lous behavior engendered by the planet’s non-uniform both of these models would be acceptable given the un- projected motion across the disk. certainty in our derived radius. A 4.5 Gyr model with a It is interesting to compare the properties of the HAT- 100M⊕ core would have a radius of only about 1.05RJ; P-2 system with the τ Boo system, which – as already such a model cannot be ruled out because it is still within noted – harbors a close-in planet with minimum mass 1 sigma of our radius determination (Table 3). Mp sin i = 3.9 MJ. Similarities of the two parent stars Fig. 3 also shows a theoretical mass-radius relation for include the nearly identical masses, effective temper- objects ranging from gas giant planets to stars (Baraffe ature, and the rapid rotation, although τ Boo, with et al. 1998, 2003). Note that HAT-P-2b falls nicely on [Fe/H] = +0.28, is somewhat more metal rich than the relation connecting giant planets to brown dwarfs to HD 147506, with [Fe/H] = +0.12. A striking differ- stars. It thus appears to be intermediate in its properties ence is that, while the orbital eccentricity of HAT-P-2b between Jupiter-like planets and more massive objects is 0.5, the eccentricity of τ Boo b is not measurably dif- like brown dwarfs or even low mass stars. ferent from zero. However, τ Boo b’s orbital period, 3.3 Compared to the other 14 known transiting planets, days, is almost half that of HAT-P-2b. A very large A Super-Massive Transiting Planet 7 fraction of close-in planets with 5

REFERENCES Bakos, G. A.,´ Lázár, J., Papp, I., Sári, P., & Green, E. M. 2002, Mandel, K., & Agol, E. 2002, ApJ, 580, L171 PASP, 114, 974 Masana, E., Jordi, C., Ribas, I. 2006, A&A, 450, 735, astro- Bakos, G. A.,´ Noyes, R. W., Kovács, G., Stanek, K. Z., Sasselov, ph/0601049 D. D., & Domsa, I. 2004, PASP, 116, 266 Perryman, M. A. C., et al. 1997, A&A, 323, L49 Bakos, G. A.,´ et al. 2007, ApJ, 656, 552 Pont, F., & Eyer, L. 2004, MNRAS, 351, 487 Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, Queloz, D., Henry, G. W., Sivan, J. P., Baliunas, S. L., Beuzit, J. A&A, 337, 403 L., Donahue, R. A., Mayor, M., Naef, D., Perrier, C., & Udry, S. Baraffe, I., Chabrier, G., Barman, T. S., Allard, F., & Hauschildt, 2001, A&A, 379, 279 P. H. 2003, A&A, 402, 701 Ram´ırez, I.,& Meléndez, J. 2005, ApJ, 626, 465 Beatty, T. G. et al. 2007, ArXiv e-prints, 704, arXiv:0704.0059 Rowe, J. F., et al. 2006, ApJ, 646, 1241 Bodenheimer, P., Laughlin, G., & Lin, D. N. C. 2003, ApJ, 592, Saar, S. H., & Donahue, R. A. 1997, ApJ, 485, 319 555 Saar, S. H., Butler, R. P., & Marcy, G. W. 1998, ApJ, 498, L153 Butler, R. P., Marcy, G. W., Williams, E., Hauser, H., & Shirts, Sato, B., et al. 2005, ApJ, 633, 465 P. 1997, ApJ, 474, L115 Shporer, A., Mazeh, T., Moran, A., Bakos, G. A.,´ Kovacs, G., & Carpenter, J. M. 2001, AJ, 121, 2851 Mashal, E. 2006, Tenth Anniversary of 51 Peg-b: Status of and Casagrande, L., Portinari, L., & Flynn, C. 2006, MNRAS, 373, 13 prospects for hot Jupiter studies, 196 Chabrier, G., Baraffe, I., Allard, F., & Hauschildt, P.2000, ApJ, Shporer, A., Tamuz, O., Zucker, S., & Mazeh, T. 2007, MNRAS, 542, L119 136, 1296 Claret, A. 2004, A&A, 428, 1001 Sienkiewicz, R., Paczyński, B., & Ratcliff, S. J. 1988, ApJ, 326, 392 Cody, A. M., & Sasselov, D. D. 2005, ApJ, 622, 704 Skrutskie, M. F., et al. 2006, AJ, 131, 1163 Demarque, P., Woo, J.-H., Kim, Y.-C., & Yi, S. K. 2004, ApJS, Southworth, J., Wheatley, P. J., & Sams, G. 2007, ArXiv e-prints, 155, 667 704, arXiv:0704.1570 Droege T. F., Richmond M. W., & Sallman M. 2006, PASP, 118, Terquem, C., Papaloizou, J. C. B., Nelson, R. P., & Lin, 1666 D. N. C. 1998, ApJ, 502, 788 Ford, E. B., & Rasio, F. A. 2007, ArXiv Astrophysics e-prints, Torres, G., Konacki, M., Sasselov, D. D., & Jha, S. 2004, ApJ, 609, arXiv:astro-ph/0703163 1071 Flower, P. J. 1996, ApJ, 469, 355 Torres, G., Konacki, M., Sasselov, D. D., & Jha, S. 2005, ApJ, 619, Fortney, J. J., Marley, M. S., & Barnes, J. W. 2006, ArXiv 558 Astrophysics e-prints, arXiv:astro-ph/0612671 Valenti, J. A., & Piskunov, N. 1996, A&AS, 118, 595 Girardi, L. et al. 2002, A&A, 391, 195 Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141 Holman, M. J. et al. 2007, ApJ, in press (arXiv:0704.2907) Vogt, S. S. 1987, PASP, 99, 1214 Kervella, P., Thévenin, Di Folco, E., & Ségransan, D. 2004, A&A, Vogt, S. S., et al. 1994, Proc. SPIE, 2198, 362 426, 297 Winn, J. N., et al. 2005, ApJ, 631, 1215 Kovács, G., Zucker, S., & Mazeh, T. 2002, A&A, 391, 369 Wright, J. T. 2005, PASP, 117, 657 Kovács, G., Bakos, G. A.,´ & Noyes, R. W. 2005, MNRAS, 356, 557 Wright, J. T., et al. 2007, ApJ, 657, 533 Latham, D. W. 1992, ASP Conf. Ser. 32: IAU Colloq. 135: Yi, S. K., Demarque, P., Kim, Y.-C., Lee, Y.-W., Ree, C. H., Complementary Approaches to Double and Multiple Star Lejeune, T., & Barnes, S. 2001, ApJS, 136, 417 Research, 32, 110 Laughlin, G., Wolf, A., Vanmunster, T., Bodenheimer, P., Fischer, D., Marcy, G., Butler, P., & Vogt, S.2005, ApJ, 621, 1072