HATP27b: A Hot Jupiter Transiting a G Star
on a 3 day orbit^{1}
Abstract
We report the discovery of HATP27b, an exoplanet transiting the moderately bright G8 dwarf star GSC 033300351 (). The orbital period is d, the reference epoch of transit is (BJD), and the transit duration is d. The host star with its effective temperature K is somewhat cooler than the Sun, and is more metalrich with a metallicity of . Its mass is and radius is . For the planetary companion we determine a mass of and radius of . For the 30 known transiting exoplanets between 0.3 and 0.8 , a negative correlation between host star metallicity and planetary radius, and an additional dependence of planetary radius on equilibrium temperature are confirmed at a high level of statistical significance.
Subject headings:
planetary systems — stars: individual (HATP27, GSC 033300351) — techniques: spectroscopic, photometric1. Introduction
Studying exoplanets is vital for understanding our own Solar
System, particularly its formation. The sample of more than 500
confirmed exoplanets
Among these planets, transiting ones are of special significance, because the transit parameters yield planetary mass and radius estimates. They also provide a means to determine some of the stellar parameters more accurately than is possible with spectroscopy alone, such as the stellar surface gravity. The more than 100 transiting exoplanets confirmed to date provide a sample large enough to draw meaningful conclusions about the planetary parameters that could not be determined by radial velocity (RV) data alone; for example, the correlation between stellar chromospheric activity and planetary surface gravity (Hartman, 2010), or the correlation of planetary parameters with host star metallicity and planetary equilibrium temperature, as described in Section 4.
The Hungarianmade Automated Telescope Network (HATNet; Bakos et al., 2011) is a system of fully automated widefield small telescopes designed to detect the small photometric dips when exoplanets transit their host stars. Since 2006, HATNet has announced and published 26 planetary systems with 28 planets in total, 26 of which transit their host stars.
Here we report the detection of our twentyseventh transiting exoplanet, named HATP27b, around the relatively bright G8 dwarf known as GSC 033300351. This planet is a textbook example of a transiting exoplanet with its radius and orbital period d being close to the median values for currently known transiting exoplanets, and with its mass of being not much less than the median mass of transiting exoplanets.
In Section 2 we present the photometric detection of the transit, along with photometric and spectroscopic followup observations of the host star HATP27. In Section 3, we describe the analysis of the data, first ruling out false positive scenarios, then determining parameters of the host star, and finally performing a global fit for all observational data. We conclude the paper by discussing HATP27b in the context of other known transiting exoplanets and investigating correlations of planetary parameters with host star metallicity and equilibrium temperature in Section 4.
2. Observations
2.1. Photometric detection
Transits of HATP27b were detected in two HATNet fields containing its host star GSC 033300351, also known as 2MASS 14510418+0556505; , , J2000, V=12.214 (Droege et al., 2006); hereafter HATP27. These fields were observed in Sloan band on a nightly basis, weather conditions permitting, from 2009 January to August, with the HAT6 and HAT10 instruments on Mount Hopkins, and with the HAT9 instrument on Mauna Kea. In total, we took 10 600 science frames with 5 minute exposure time and 5.5 minute cadence. For approximately 1200 of the images, photometric measurements of individual stars had significant error, therefore these frames were rejected. Each image contains approximately 20 000 stars down to . For the brightest stars, the perpoint photometric precision is approximately 4.5 mmag.
Calibration, astrometry, aperture photometry, External Parameter Decorrelation (EPD), the Trend Filtering Algorithm (TFA) and the Boxfitting Least Squares method were applied to the data as described in Bakos et al. (2010a). We detected a transit signature in the light curve of HATP27, with a signal depth of mmag and period of days. This presumed transit has relative firsttolastcontact duration , corresponding to total duration days ( hours). The folded light curve is presented in Figure 1.
2.2. Reconnaissance Spectroscopy
Reconnaissance spectra (Latham et al., 2009a) of HATP27 were taken using three facilities: the Tillinghast Reflector Echelle Spectrograph (TRES; Fűrész, 2008) on the 1.5 m Tillinghast Reflector at FLWO, the echelle spectrograph on the 2.5 m du Pont telescope at Las Campanas Observatory (LCO) in Chile, and the echelle spectrograph on the Australian National University (ANU) 2.3 m telescope at Siding Spring Observatory (SSO) in Australia. We gathered two spectra of HATP27 with TRES in 2009 July and 2010 February, two spectra with the du Pont telescope in 2009 July, and fourteen spectra with the ANU 2.3 m telescope in 2009 July. The exact dates and the results of these observations are summarized in Table 1.
Instrument  Date  Number of 



Spectra  (K)  (cgs)  ()  ()  
TRES  2009 Jul 05  1  
du Pont  2009 Jul 10  1  
du Pont  2009 Jul 11  1  
ANU 2.3 m  2009 Jul 16  5  
ANU 2.3 m  2009 Jul 17  6  
ANU 2.3 m  2009 Jul 18  3  
TRES  2010 Feb 13  1 
Following Quinn et al. (2010) and Buchhave et al. (2010), we calculated the mean radial velocity, effective temperature, surface gravity, and projected rotational velocity of the host star, based on spectra taken by TRES. The inferred radial velocity RMS residual of is consistent with no detectable RV variation within the precision of the measurements. We established the following atmospheric parameters for HATP27: effective temperature K, surface gravity (cgs), and projected rotational velocity , indicating a G8 dwarf. The mean heliocentric RV of HATP27 after subtracting the gravitational redshift of the Sun is .
Because this is the first time we used the du Pont 2.5 m and ANU 2.3 m telescopes for reconnaissance spectroscopy of HATNet targets, we briefly describe the instruments and our data reduction procedure.
The spectrograph on the du Pont telescope was used with a
long and wide slit. The obtained spectra have
wavelength coverage 3700–7000 Å at a resolution of
. During the first observation
the seeing ranged between 2–3″ and we used an exposure time of
1200 s, which provided electrons per resolution element at
the wavelength of 5187 Å. The seeing during the second observation
was and we used an exposure time of 150 s to obtain a lower
S/N spectrum sufficient to detect a velocity variation of several
. We obtained a ThAr lamp spectrum before and after each
observation to use in determining the wavelength solution. We used the
CCDPROC package of IRAF
The extracted du Pont spectra were then crosscorrelated against a library of synthetic stellar spectra to estimate the effective temperature, surface gravity, projected rotational velocity, and radial velocity of the star. We followed a procedure similar to that described by Torres et al. (2002), using the same synthetic templates, but broadened to the resolution of the du Pont echelle. These templates, which were generated for the CfA Digital Speedometer (Latham, 1992), only cover a wavelength range of 5150–5360 Å, so we restricted our analysis to a single order of the spectrum covering a similar range.
Spectra were also taken using the echelle spectrograph on the ANU 2.3 m telescope. The echelle was used in a standard configuration with a wide slit and crossdisperser setting of , which delivered 27 full orders between 3900–6720 Å with a nominal spectral resolution of . The CCD is a 2K2K format with pixels. The gain was two electrons per ADU resulting in a read noise of approximately 2.3 ADU for each pixel. The spectra were binned by two along the spatial direction. A total of fourteen 1200 s exposures were taken. The seeing ranged from to . The signaltonoise on a single pixel was between 5 and 10 for each individual exposure. ThAr lamp calibration exposures were taken every hour for wavelength calibration. A high signaltonoise exposure was also taken of the bright radial velocity standard star HD 223311.
Spectra were reduced using tasks in the IRAF packages CCDPROC and ECHELLE. The spectra were crosscorrelated against the radial velocity standard star HD 223311 using the IRAF task FXCOR in the RV package. We used at least 20, typically 25 of the 27 orders for the crosscorrelation, rejecting the bluest orders for many exposures where the signaltonoise was too low. Each spectral order was treated separately and the mean of the velocities from the individual orders was calculated. Their standard deviations were less than for each exposure.
Each night, the exposures were taken within a two hour interval, much shorter than the orbital period of the presumed companion. For detecting large radial velocity variations to rule out the possibility of an eclipsing binary, we consider the mean radial velocities per night. The standard deviations between exposures were less than for each night. Stellar parameters could have been estimated only with large uncertainty based on data with such low signaltonoise. Therefore these parameters are not calculated from the ANU 2.3 m observations.
The results of the observations taken with these three telescopes are listed in Table 1. Note that for each telescope, the radial velocity measurement uncertainty is much higher than the radial velocity variations of the Sun due to Solar System bodies. Therefore we only calculated heliocentric radial velocities of the target. For the more accurate measurements described in the next section, we will use barycentric radial velocities instead. This accuracy, however, is enough to rule out the case of an eclipsing binary star with great certainty. The largest RV variation within an instrument was only 2 , much less than the orbital speed in a typical binary system. Note that the zeropoint shift between instruments is as large as , due to the different methods used for data reduction. Also, all eighteen spectra were singlelined and spectral lines were symmetric, providing no evidence for additional stars in the system up to the precision of the measurements.
2.3. High resolution, high S/N spectroscopy
We acquired highresolution, highS/N spectra of HATP27 using the HIRES instrument (Vogt et al., 1994) on the Keck I telescope located on Mauna Kea, Hawaii, between 2009 December and 2010 June. The spectrometer was configured with the wide slit, yielding a resolving power of over the wavelength range of 3800–8000 Å.
Nine exposures were taken using an gas cell (see Marcy & Butler, 1992), and a single template exposure was obtained without the absorption cell. We followed Butler et al. (1996) to establish RVs in the Solar System barycentric frame. We also calculated the index for each spectrum (following Isaacson & Fischer, 2010). The resulting values and their uncertainties are listed in Table 2. They are plotted periodfolded in Figure 2, together with the fit established in Section 3.
The effective temperature established later in Section 3.2 compared to Figure 4 of Valenti & Fischer (2005) implies for the star. This can in turn be used in the formula given in Noyes et al. (1984) together with the median of 0.231 to conclude . This activity value is consistent with the RV jitter value established in Section 3.3, according to the observations given by Wright (2005). Based on Figure 9 in Isaacson & Fischer (2010), this value qualifies HATP27 as chromospherically active relative to California Planet Search targets of the same spectral class. The activity index does not correlate significantly with orbital phase.
BJD (UTC)  RV 

BS 

Phase  

()  ()  ()  ()  ()  
Note. – For the iodinefree template exposures there is no RV measurement, but the BS and index can still be determined.
2.4. Photometric followup observations
A highprecision photometric followup of a complete transit was carried out, permitting refined estimates of the light curve parameters and thus orbital and planetary properties: The transit of HATP27 was observed on the night of 2010 March 2 with the KeplerCam CCD camera on the FLWO 1.2 m telescope. We acquired 168 science frames in Sloan band with 60 second exposure time, 73 second cadence.
Following the procedure described by Bakos et al. (2010a), these images were first calibrated, then astrometry and aperture photometry was performed to arrive at light curves, which were finally cleaned of trends using EPD and TFA, carried out simultaneously with the global modeling described in Section 3.3. The result is shown in Figure 3, along with the bestfit transit light curve; the individual measurements are reported in Table 3.
BJD (UTC)  Mag 
Mag(orig) 
Filter  

()  
Note. – This table is available in a machinereadable form in the online journal. A portion is shown here for guidance regarding its form and content.
3. Analysis
3.1. Excluding blend scenarios
To further exclude possible blends, we perform bisector analysis the same way as in §5 of Bakos et al. (2007). A significant scatter is found, strongly correlated to the presence of moonlight, which we account for using the method described by Hartman et al. (2011). The bisector spans, corrected for the effect of the moonlight, are shown in the third panel of Figure 2. They do not exhibit significant correlation with the RV values, and the RMS scatter of the bisector spans (4.6 ) is significantly smaller than the RV amplitude. These findings rule out a blend scenario with high certainty, implying that the measured photometric and spectroscopic features are due to a planet orbiting HATP27.
3.2. Properties of the parent star
We first determine spectroscopic parameters of HATP27, which will allow us to calculate stellar mass and radius. The Spectroscopy Made Easy analysis package (SME; Valenti & Piskunov, 1996) is used to establish the effective temperature, metallicity and stellar surface gravity based on the Keck/HIRES template spectrum of HATP27, using atomic line data from the database of Valenti & Fischer (2005). After an initial estimate for these parameters, we perform a Monte Carlo calculation, relying also on the normalized semimajor axis inferred from transit light curves, for the reasons described by Bakos et al. (2010b). The final values adopted after two iterations are K, , and , also listed in Table 4.
Based on the final spectroscopic parameters the model isochrones yield a stellar mass = and radius = for HATP27, along with other properties listed at the bottom of Table 4. These values classify the star as a G8 dwarf, and suggest an age of Gyr. The model isochrones of Yi et al. (2001) for a metallicity of are plotted in Figure 4, along with the best estimate of and of HATP27, and their and confidence ellipsoids. For comparison, the initial SME result, corresponding to a somewhat younger state, is also indicated.
The intrinsic absolute magnitude predictions of this model
(given in the ESO photometric system)
can be compared to observations to calculate the distance of HATP27.
For this we use the nearinfrared brightness measurements from
the 2MASS Catalogue (Skrutskie et al., 2006):
,
and
.
These values are converted to ESO
(Carpenter, 2001), then
compared to the absolute magnitude estimates to calculate the
distance. We account for interstellar dust extinction in the line of
sight using from the dust map by
Schlegel et al. (1998)
Parameter  Value  Source 

Spectroscopic properties  
(K)  SME 

SME  
()  SME  
()  SME  
()  SME  
()  TRES  
Photometric properties  
(mag)  12.214  TASS 
(mag)  TASS  
(mag)  2MASS  
(mag)  2MASS  
(mag)  2MASS  
Derived properties  
()  YY++SME 

()  YY++SME  
(cgs)  YY++SME  
()  YY++SME  
(mag)  YY++SME  
(mag, ESO)  YY++SME  
Age (Gyr)  YY++SME  
Distance (pc)  YY++SME 
3.3. Global modeling of the data
We fit the combined model described by Bakos et al. (2010a) to the HATNet photometry, followup photometry, and radial velocity measurements simultaneously. The eight main parameters describing the model are the time of the first observed transit center, the time of the followup transit center, the normalized planetary radius , the square of the impact parameter , the reciprocal of the half duration of the transit , the RV semiamplitude , and the Lagrangian elements and (where is the longitude of periastron).
Instrumental parameters of the model include the HATNet outoftransit magnitude and the relative zeropoint of the Keck RVs. The joint fit provides the full a posteriori probability distributions of all variables (including ), which are used to update the limbdarkening coefficients for another iteration of the joint fit. This leads to estimated distributions for the stellar, light curve, and RV parameters, which are combined to calculate values for planetary parameters. These final values are summarized in Table 5.
The orbital eccentricity is consistent with zero (using the method of Lucy & Sweeney, 1971 we find that there is a 25% conditional probability of detecting an eccentricity of at least 0.078 given a circular orbit and an uncertainty of 0.047). Nevertheless, we stress that throughout the global modeling, the orbit was allowed to be eccentric, and all system parameters and their respective errors inherently contain the uncertainty arising from the floating and values.
Parameter  Value 

Light curve parameters  
(days)  


(days)


(days)


(degree)  
Limbdarkening coefficients 

(linear term)  
(quadratic term)  
RV parameters  
()  




(degree)  
RV jitter ()


Secondary eclipse parameters  
(BJD)  
Planetary parameters  
()  
()  


()  
(cgs)  
(AU)  
(K)  


( )


( )


( )

4. Discussion
4.1. Properties of HATP27b
Figure 5 presents the currently known transiting exoplanets and Solar System gas planets on a mass—radius diagram, with HATP27b highlighted. Also shown are the planetary isochrones of Fortney et al. (2007) interpolated to the insolation of HATP27 at the orbit of HATP27b. Taking into consideration the age established in Section 3.2, the planetary parameters are consistent with a hot Jupiter with a 10 core in a 4 Gyr old system.
HATP27b can be seen to lie inside the large accumulation of planets with similar masses and radii. To further compare it to other Hot Jupiters, in Figure 6 we present histograms of mass, radius and period for the 112 transiting exoplanets confirmed to date.
When comparing these parameters, we note that there is only one transiting exoplanet known that is more massive and has a smaller radius and a smaller period than HATP27b: this is HATP20b with 7.246 , 0.867 on a 2.875 d orbit (Bakos et al., 2010b). This means that HATP27b is less inflated than other planets of similar mass and orbital period, possibly due to a larger than average core.
Regarding eccentricity, there are 31 transiting exoplanets known under 8 with an orbital period within 0.5 d of that of HATP27b, out of which 8 – more than a quarter of them – are thought to be eccentric. This hints that there is a possibility for the orbit of HATP27b to be eccentric as well, justifying our choice not to fix eccentricity to zero in Section 3. Future observations of radial velocity or occultation timing would be required to determine whether the orbit is indeed eccentric.
The impact parameter of HATP27b is unusually large. As Ribas et al. (2008) pointed out, such a near grazing transit has the advantage of its depth and duration being more sensitive to the presense of further planetary companions on inclined orbits. This makes HATP27 a promising target for transit timing variation and transit duration variation studies.
Knutson et al. (2010) found a strong negative correlation between chromospheric activity of the host star and temperature inversion in the planetary atmosphere. However, since early type stars dominate magnitude limited surveys, cool, that is, active planetary hosts are rare. The bottom right panel of Figure 6 shows that HATP27 is relatively active compared to known planetary hosts for which has been reported, making it an exciting target for Spitzer Space Telescope to test this correlation.
4.2. Correlation of planetary parameters with host star metallicity
The relation between host star metallicity (denoted as for clarity, not to be confused with the assumed metal content of the planet) and planetary composition was studied by Guillot et al. (2006). A positive correlation, with Pearson correlation coefficient , was found between the inferred mass of the planetary core and stellar metallicity for the seven transiting exoplanets known at that time with positive inferred core mass. The idea is that planets have formed from the same cloud as their host stars, and therefore their metal content should correlate. However, it is not clear how stellar metallicity is connected to planetary metallicity, especially because a larger rocky core is likely to accrete more gas during the planet’s formation.
Burrows et al. (2007) also investigated this relation, based on 12 transiting exoplanets known at the time. They used an atmospheric opacity dependent core mass model to explain radius anomalies. They also found a strong correlation between host star metallicity and inferred core mass, but the correlation coefficient was not reported.
Enoch et al. (2010) found that there is a strong negative correlation with between and for the 18 known transiting exoplanets below 0.6 , whereas this correlation is negligible for more massive planets. This can be explained by noticing that the theoretical planet models of Fortney et al. (2007), Bodenheimer et al. (2003), and Baraffe et al. (2008) all suggest that the radius of a planet is more sensitive to its composition for low mass planets than it is for more massive ones.
name  reference  

WASP21b  Bouchy et al. (2010)  
HD 149026b  Ammlervon Eiff et al. (2009),  
Carter et al. (2009)  
Kepler7b  Latham et al. (2009b),  
Kipping & Bakos (2010)  
WASP13b  Skillen et al. (2009)  
Kepler8b  Jenkins et al. (2010),  
Kipping & Bakos (2010)  
CoRoT5b  Rauer et al. (2009)  
WASP31b  Anderson et al. (2010)  
WASP11/HATP10b  Bakos et al. (2009)  
WASP17b  Anderson et al. (2010)  
WASP6b  Gillon et al. (2009)  
HATP1b  Torres et al. (2008),  
Ammlervon Eiff et al. (2009)  
HATP17b  Howard et al. (2010)  
WASP15b  West et al. (2009)  
OGLETR111b  Santos et al. (2006),  
Torres et al. (2008),  
Adams et al. (2010)  
HATP4b  Kovács et al. (2007),  
Torres et al. (2008),  
Winn et al. (2010)  
WASP22b  Maxted et al. (2010)  
XO2b  Torres et al. (2008),  
Ammlervon Eiff et al. (2009)  
HATP25b  Quinn et al. (2010)  
WASP25b  Enoch et al. (2010)  
WASP34b  Smalley et al. (2010)  
HATP3b  Torres et al. (2007),  
Torres et al. (2008)  
HATP28b  Buchhave et al. (2011)  
HATP27b  this paper  
HATP24b  Kipping et al. (2010)  
HD 209458b  Laughlin et al. (2005),  
Torres et al. (2008)  
Kepler6b  Dunham et al. (2010),  
Kipping & Bakos (2010)  
OGLETR10b  Torres et al. (2008)  
CoRoT4b  Moutou et al. (2008)  
TrES1b  Torres et al. (2008)  
HATP9b  Shporer et al. (2009) 
In this subsection, we examine further the correlation between host star metallicity and planetary mass or planetary radius. We use a substantially expanded sample of 30 known transiting exoplanets with masses between 0.3 and 0.8 , see Table 6. The upper limit is selected to exclude more massive planets whose radius is expected to depend less on the composition, see above. We explain the role of the lower limit and the effect of the two lowest mass planets in Table 6 later.
The null hypothesis is that the host star metallicity and the selected planetar parameter are independent. The alternative hypothesis is that they are related by some underlying phenomenon. A false positive, also known as an error of the first kind, is rejection of the null hypothesis in spite of it being true. We implement three independent statistical methods to estimate the false positive probability: test, bootstrap technique and test. We denote the probability estimates by , and , respectively. This is the statistical significance of the correlation: the lower this probability is, the more confidently the null hypothesis (i.e., no correlation) can be rejected.
For the test, we assume that the investigated parameters have normal distribution. For each set of data pairs, we calculate the value from the correlation coefficient and sample size using the formula
The conditional distribution of this variable given the null hypothesis is Student’s distribution with degrees of freedom (Press et al., 1992, p. 640). Then the estimate for false positive probability is determined by looking up the twotailed probability of this distribution yielding larger absolute value that the one measured. For comparison, we also performed the test for the samples and parameters studied by Guillot et al. (2006) and Enoch et al. (2010). The resulting values are listed in Table 7.
For the sample set of Table 6, we also implement the bootstrap resampling technique (Efron & Tibshirani, 2003). This has the advantage that no assumption about the a priori distribution of the data is necessary. To perform bootstrap resampling, consider the data , where is the host star metallicity, and is the corresponding planetary mass or radius. Again, we would like to calculate an estimate of the probability that a sample of similar distribution, but independent parameters for each pair, has a correlation coefficient that exceeds that of our measurements in absolute value. For this, we build 10 000 000 sample sets of pairs by drawing and values independently with replacements from the set of measured and values, respectively. The percentile rank of the absolute value of the correlation coefficient for the measured data among the random samples gives our estimates , listed in Table 7.
Finally, we test these correlations with an additional method, the test (see e.g. Lupton, 1993, p. 100). This requires that the null hypothesis (no correlation) be nested in the tested hypothesis (linear correlation), which indeed is the case. Let RSS denote the residual sum of squares for the best fit of the null hypothesis, that is, the variance of about its mean, and RSS denote the residual sum of squares of the linear fit. The no correlation model has free parameters: the mean, whereas the linear fit has . The conditional distribution of
given the null hypothesis is the distribution with degrees of freedom. This enables us to calculate , the third estimate for the false positive probability.
The estimates , , given by the three methods are listed in Table 7. For each correlation, they coincide up to the uncertainty of the methods. The values reflect the significant and correlations found by Guillot et al. (2006) and Enoch et al. (2010).
As for the 30 planets listed in Table 6, it is important to note that the correlations depend strongly on the choice of the lower mass limit. The two least massive planets in the table are WASP21b with a mass of , very low host star metallicity of , and average radius of ; and the dense HD 149026b with a mass of , high host star metallicity of , and low radius of . Increasing the lower mass limit for our sample first excludes WASP21b, which would much support the positive — correlation with its low mass and low host star metallicity. Further increasing the limit then excludes HD 149026b, which would much weaken it with its low mass and high host star metallicity. To have an unbiased result, outliers cannot be excluded without a justified reason, therefore we need to compare the false positive probabilities of the three nested samples. They scatter between 15% and 58%, neither supporting, nor rejecting a — correlation.
Similarily, both WASP21b and HD 149026b have a strong effect on the negative — correlation, because of the extreme value of their host star metallicities. In this case, we see that the maximum of the false positive probabilities is , therefore this correlation is statistically significant for all our choices of lower mass limits. This is at least a fivefold improvement over the sample investigated by Enoch et al. (2010), due to the larger sample size.
Restriction 

Planetary 






on planets  parameter  
7  0.  78 
3.  9%  
18  0.  53 
2.  4%  
0.3  30  0.  270  15.  0%  15.  0%  15.  0%  
0.35  29  0.  106  58.  4%  58.  2%  58.  4%  
0.4  28  0.  247  20.  4%  20.  4%  20.  4%  
0.3  30  0.  505  0.  44%  0.  43%  0.  44%  
0.35  29  0.  620  0.  03%  0.  03%  0.  034%  
0.4  28  0.  575  0.  14%  0.  14%  0.  14% 
4.3. Dependence on planetary equilibrium temperature
Other factors like insolation are likely to influence planetary radius as well, see e.g. Fortney et al. (2007), Kovács et al. (2010), Enoch et al. (2010), and Faedi et al. (2011). To further investigate this relation, we compare two models: for null hypothesis, we accept the linear planetary radius—host star metallicity relation of the previous section:
(1) 
The second model – alternative hypothesis – is similar to that of Enoch et al. (2010), accounting for the equilibrium temperature in addition to the host star metallicity:
(2) 
The equilibrium temperature of the planet is calculated from the timeaveraged stellar flux on its orbit, assuming gray body spectrum for the planets, and neglecting tidal and other heating mechanisms. For simplicity, we now include all 30 planets of Table 6 in our models. With the best fit parameters, the two models are
To quantify the statistical significance, we apply the test again. The resulting false positive probability is 0.18%. This means that once we accept the dependence of planetary radius on host star metallicity, then the probability of such a correlation with planetary equilibrium temperature if they were not physically related is only 0.18%. This strongly supports the threeparameter linear fit model. For reference, the residual sums of squares for the one, two and threeparameter fits of the — data are , and , respectively.
Figure 7 presents versus metallicity for the 30 planets. Equation (2) predicts this quantity to be , which is also plotted. HATP27b apparently follows the model’s prediction. For reference, the correlation coefficient between the displayed transformed variables is now , which has a larger absolute value than between and metallicity, as expected.
This analysis supports the statement that planetary radius depends on equilibrium temperature in addition to host star metallicity, as found by Enoch et al. (2010) and Faedi et al. (2011). However, this correlation does not imply that insolation itself would inflate planets: the underlying phenomenon could be related to anything correlated to equilibrium temperature, or equivalently, orbital radius. For instance, Batygin & Stevenson (2010) suggest that it is Ohmic dissipation in the interior of the planet that inflates hot Jupiters. This theory is further supported by Laughlin et al. (2011).
Altogether, HATP27b is an important addition to the growing sample of lowmass Jupiters. It orbits a metal rich star, and supports the suggested correlations between host star metallicity, planetary equilibrium temperature, and planetary radius. Also, HATP27 is chromospherically active, providing an excellent case for refining the confidence level of the hypothesized correlation between stellar activity and planetary temperature inversion.
Footnotes
 affiliation: Based in part on observations obtained at the W. M. Keck Observatory, which is operated by the University of California and the California Institute of Technology. Keck time has been granted by NOAO (A201Hr) and NASA (N018Hr, N167Hr).
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: email: bbeky@cfa.harvard.edu
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: Departamento de Astronomía y Astrofísica, Pontificia Universidad Católica de Chile, Santiago, Chile
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: Departamento de Astronomía y Astrofísica, Pontificia Universidad Católica de Chile, Santiago, Chile
 affiliation: Research School of Astronomy and Astrophysics, The Australian National University, Weston Creek, ACT, Australia
 affiliation: Konkoly Observatory, Budapest, Hungary
 affiliation: Sydney Institute for Astronomy, School of Physics, University of Sydney, Sydney, Australia
 affiliation: Konkoly Observatory, Budapest, Hungary
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: Department of Astronomy, University of California, Berkeley, CA
 affiliation: Department of Astronomy, University of California, Berkeley, CA
 affiliation: Department of Astronomy, Yale University, New Haven, CT
 affiliation: Astronomy Department, California Institute of Technology, Pasadena, CA
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: Niels Bohr Institute, Copenhagen University, Denmark
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: Harvard–Smithsonian Center for Astrophysics, Cambridge, MA
 affiliation: Department of Physics, University of South Dakota, Vermillion, SD
 affiliation: Hungarian Astronomical Association, Budapest, Hungary
 affiliation: Hungarian Astronomical Association, Budapest, Hungary
 affiliation: Hungarian Astronomical Association, Budapest, Hungary
 According to http://www.exoplanet.eu/catalogall.php.
 The mean heliocentric RV of the target. Systematic differences between the velocities from different instruments are consistent with the velocity zeropoint uncertainties.
 IRAF is distributed by the National Optical Astronomical Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation.
 The zeropoint of these velocities is arbitrary. An overall offset fitted to these velocities in Section 3.3 has not been subtracted.
 Internal errors excluding the component of astrophysical/instrumental jitter considered in Section 3.3.
 The bisector spans have been corrected for sky contamination following Hartman et al. (2011).
 Relative chromospheric activity index, calibrated to the scale of Vaughan, Preston & Wilson (1978).
 The outoftransit level has been subtracted. These magnitudes have been obtained by the EPD and TFA procedures, carried out simultaneously with the transit fit.
 Raw magnitude values without application of the EPD and TFA procedures.
 http://irsa.ipac.caltech.edu/applications/DUST
 SME = Spectroscopy Made Easy package for the analysis of highresolution spectra (Valenti & Piskunov, 1996). These parameters rely primarily on SME, but have a small dependence also on the iterative analysis incorporating the isochrone search and global modeling of the data, as described in the text.
 YY++SME = Based on the YY isochrones (Yi et al., 2001), as a luminosity indicator, and the SME results.
 Reference epoch of mid transit that minimizes the correlation with the orbital period.
 Total transit duration, time between first to last contact.
 Ingress/egress time, time between first and second, or third and fourth contact.
 Adopted from the tabulations by Claret (2004) according to the spectroscopic (SME) parameters listed in Table 4.
 Lagrangian orbital parameters derived from the global modeling, and primarily determined by the RV data.
 Lagrangian orbital parameters derived from the global modeling, and primarily determined by the RV data.
 The contribution of the intrinsic stellar jitter and possible instrumental errors that needs to be added in quadrature to the calculated RV uncertainties so that in the joint fit.
 Correlation coefficient between the planetary mass and radius.
 The Safronov number is given by Hansen & Barman (2007) as .
 Stellar irradiation flux per unit surface area at periastron, apastron and timeaveraged over the orbit, respectively.
 Stellar irradiation flux per unit surface area at periastron, apastron and timeaveraged over the orbit, respectively.
 Stellar irradiation flux per unit surface area at periastron, apastron and timeaveraged over the orbit, respectively.
 sample size
 correlation coefficient
 estimates for false positive probability given by test, bootstrap method and test, respectively
 estimates for false positive probability given by test, bootstrap method and test, respectively
 estimates for false positive probability given by test, bootstrap method and test, respectively
 reported by Guillot et al. (2006)
 reported by Enoch et al. (2010)
References
 Adams, E. R., LópezMorales, M., Elliot, J. L., Seager, S., & Osip, D. J. 2010, arXiv:1003.0457
 Ammlervon Eiff, M., Santos, N. C., Sousa, S. G., Fernandes, J., Guillot, T., Israelian, G., Mayor, M., & Melo, C. 2009, A&A, 507, 523
 Anderson, D. R., et al. 2010, ApJ, 709, 159
 Anderson, D. R., et al. 2010, arXiv:1011.5882
 Bakos, G. Á., et al. 2007, ApJ, 670, 826
 Bakos, G. Á., et al. 2009, ApJ, 696, 1950
 Bakos, G. Á., et al. 2010a, ApJ, 710, 1724
 Bakos, G. Á., et al. 2010b, ApJ, submitted (arXiv:1008.3388B)
 Bakos, G. Á., et al. 2011, ApJ, submitted (arXiv:1101.0322B)
 Baraffe, I., Chabrier, G., & Barman, T. 2008, A&A, 482, 315
 Batygin, K. & Stevenson, D. J. 2010, ApJ, 714, 238
 Bodenheimer, P., Laughlin, G., Lin, D. N. C. 2003, ApJ, 592, 555
 Bonifacio, P., Monai, S., & Beers, T. C. 2000, AJ, 120, 2065
 Bouchy, F., et al. 2010, A&A, 519, 98
 Buchhave, L. A., et al. 2010, ApJ, 720, 1118
 Buchhave, L., et al. 2011, ApJ, submitted
 Burrows, A., Hubeny, I., Budaj, J., & Hubbard, W. B. 2007, ApJ, 661, 502
 Butler, R. P. et al. 1996, PASP, 108, 500
 Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, AJ, 345, 245
 Carpenter, J. M. 2001, AJ, 121, 2851
 Carter, J. A., Winn, J. N., Gilliland, R., & Holman, M. J. 2009, ApJ, 696, 241
 Claret, A. 2004, A&A, 428, 1001
 Droege, T. F., Richmond, M. W., & Sallman, M. 2006, PASP, 118, 1666
 Dunham, E. W., et al. 2010, arXiv:1001.0333
 Efron, B. & Tibshirani, R. J. 2003, An Introduction to the Bootstrap (Chapman & Hall: New York, NY)
 Enoch, B., et al. 2010, submitted to MNRAS (arXiv:1009.5917)
 Faedi, F., et al. 2011, submitted to A&A (arXiv:1102.1375)
 Fortney, J. J. 2010, BAAS, 42, 4701
 Fortney, J. J., Marley, M. S., & Barnes, J. W. 2007, ApJ, 659, 1661
 Fűrész, G. 2008 Ph.D. thesis, University of Szeged, Hungary
 Gillon, M., et al. 2009, A&A, 501, 785
 Guillot, T., Santos, N. C., Pont, F., Iro, N., Melo, C., & Ribas, I. 2006, A&A, 453, L21
 Hansen, B. M. S., & Barman, T. 2007, ApJ, 671, 861
 Hartman, J. D. 2010, ApJ, 717, 138
 Hartman, J. D., et al. 2011, ApJ, 726, 52
 Howard, A. W., et al. 2010, arXiv:1008.3898
 Ida, S., & Lin, D. N. C. 2004, BAAS, 35, 0105
 Ida, S., & Lin, D. N. C. 2005, ApJ, 626, 1045
 Ida, S., & Lin, D. N. C. 2008, ApJ, 673, 487
 Isaacson, H., & Fischer, D. 2010, ApJ, 725, 875
 Jenkins, J. M., et al. 2010, arXiv:1001.0416
 Kipping, D. M., & Bakos, G. Á. 2010, arXiv:1004.3538
 Kipping, D. M., et al. 2010, arXiv:1008.3389
 Knutson, Heather A., Howard, A. W., Isaacson, H. 2010, ApJ, 720, 1569
 Kovács, G., et al. 2007, ApJ, 670, L41
 Kovács, G., et al. 2010, ApJ, 724, 866
 Laughlin, G., Crismani, M., and Adams, F. C. 2011, submitted to ApJ (arXiv:1101.5827)
 Latham, D. W. 1992, in IAU Coll. 135, Complementary Approaches to Double and Multiple Star Research, ASP Conf. Ser. 32, eds. H. A. McAlister & W. I. Hartkopf (San Francisco: ASP), 110
 Latham, D. W., et al. 2009, ApJ, 704, 1107
 Latham, D. W., et al. 2009, arXiv:1001.0190
 Laughlin, G., Marcy, G. W., Vogt, S. S., Fischer, D. A., & Butler, R. P. 2005, ApJ, 629, L121
 Lucy, L. B., & Sweeney, M. A. 1971, AJ, 76, 544
 Lupton, R. 1993, Statistics in Theory and Practice (Princeton University Press: Princeton, NJ)
 Marcy, G. W., & Butler, R. P. 1992, PASP, 104, 270
 Mardling, R. A. 2007, MNRAS, 382, 1768
 Maxted, P. F. L., et al. 2010, arXiv:1004.1514
 Moutou, C., et al. 2008, A&A, 488, L47
 Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., & Vaughan, A. H. 1984, ApJ, 279, 763
 Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical Recipes in C: The Art of Scientific Computing (2nd ed.; New York: Cambridge University Press)
 Quinn, S. N., et al. 2010, ApJ, submitted (arXiv:1008.3565)
 Rauer, H., et al. 2009, A&A, 506, 281
 Ribas, I., FontRibera, A., & Beaulieu, J. 2008, ApJ, 677, 525
 Santos, N. C., et al. 2006, A&A, 450, 825
 Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525
 Shporer, A., et al. 2009, ApJ, 690, 1393
 Skillen, I., et al. 2009, A&A, 502, 391
 Skrutskie, M. F., et al. 2006, AJ, 131, 1163
 Smalley, B., et al. 2010, arXiv:1012.2278
 Torres, G., Neuhäuser, R., & Guenther, E. W. 2002, AJ, 123, 1701
 Torres, G., Winn, J. N., & Holman, M. J. 2008, ApJ, 677, 1324
 Torres, G., et al. 2007, ApJ, 666, L121
 Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141
 Valenti, J. A., & Piskunov, N. 1996, A&AS, 118, 595
 Vogt, S. S. et al. 1994, Proc. SPIE, 2198, 362
 Yi, S. K. et al. 2001, ApJS, 136, 417
 Vaughan, A. H., Preston, G. W., & Wilson, O. C. 1978, PASP, 90, 267
 West, R. G., et al. 2009, AJ, 137, 4834
 Winn, J. N., et al. 2010, arXiv:1010.1318
 Wright, J. T. 2005, PASP, 117, 657