Radii and Temperatures of Exoplanet Host Stars

Directly Determined Linear Radii and Effective Temperatures
of Exoplanet Host Stars

Gerard T. van Belle11affiliation: European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching; gerard.van.belle@eso.org Kaspar von Braun22affiliation: NASA Exoplanet Science Institute, California Institute of Technology, MC 100-22, Pasadena, CA 91125; kaspar@ipac.caltech.edu

We present interferometric angular sizes for 12 stars with known planetary companions, for comparison with 28 additional main-sequence stars not known to host planets. For all objects we estimate bolometric fluxes and reddenings through spectral energy distribution fits, and in conjunction with the angular sizes, measurements of effective temperature. The angular sizes of these stars are sufficiently small that the fundamental resolution limits of our primary instrument, the Palomar Testbed Interferometer, are investigated at the sub-milliarcsecond level and empirically established based upon known performance limits. We demonstrate that the effective temperature scale as a function of dereddened color is statistically identical for stars with and without planets. A useful byproduct of this investigation is a direct calibration of the scale for solar-like stars, as a function of both spectral type and color, with an precision of K over the range and K for the range F6V – G5V. Additionally, we provide in an appendix spectral energy distribution fits for the 166 stars with known planets which have sufficient photometry available in the literature for such fits; this derived “XO-Rad” database includes homogenous estimates of bolometric flux, reddening, and angular size.

Subject headings:
infrared: stars, stars: fundamental parameters, techniques: interferometric
slugcomment: Submitted to ApJ on 15.10.2008, accepted on 27.12.2008

1. Introduction

The formation, evolution, and environment of extrasolar planets are heavily influenced by their respective parent stars, including the location and extent of the habitable zone. To provide constraints on the characterization of these planets, it is therefore of significant scientific value to directly determine the astrophysical parameters of the host stars. Of particular interest are stellar radius and effective surface temperature since these two parameters help uniquely characterize our knowledge of extrasolar planet enviroments. In the case of radius, planetary radii are frequently not directly measured but established through observations of transit events as a ratio of planet to stellar radius. Measurements of planetary temperature are directly linked to the spectral characteristics of the star irradiating the planet.

For extrasolar planet hosting stars (EHSs) that can be resolved with interferometers, their angular sizes are directly measured. Since the Stefan-Boltzman Law (Stefan, 1879; Boltzmann, 1884) can be rewritten as , where is the reddening-corrected bolometric flux, the effective temperature can be directly measured for these stars. We obtained data with the Palomar Testbed Interferometer (PTI) for 9 nearby EHSs with the aim of directly measuring their angular diameters, and computed estimates of their through spectral energy distribution (SED) fitting to their available literature photometry. Additional EHS angular diameters from the Center for High Angular Resolution Astronomy (CHARA) Array are also folded into this investigation (Baines et al., 2007, 2008a). Further interferometric work relevant to the diameters of EHSs can be found in Mozurkewich et al. (1991, 2003).

Observational biases cause a large fraction of known EHSs to be nearby, enabling the use of Hipparcos parallaxes for a direct determination of their distances (Perryman et al., 1997; Perryman & ESA, 1997); in combination with angular size measurements, their linear radii can be determined.111It is misleading, however, to indicate that interferometric angular size measurements independently lead to characterizations of stellar luminosity. A common mistake is to assume that radius and temperature measurements derived from a single interferometric angular size can be combined through use of the Stefan-Boltzman law to ‘measure’ luminosity. A cursory examination of the relationship between angular size and radius and temperature will demonstrate the new information contained in an angular size measurement is discarded when calculating : only bolometric flux and distance information affect measures of .

The aim of this publication is to provide directly determined and astrophysical parameters of these 12 EHSs along with equivalently derived parameters for a control group of 28 main sequence stars not currently known to host extrasolar planets. In addition, we present estimates of astrophysical parameters for all currently known EHSs with sufficient literature photometry (166 of the 230 known)222As of Feb. 1, 2008.. The literature photometry and aforementioned SED fitting provides, for the sample of 166 EHS stars, estimates of and , done in the same way as done by van Belle et al. (2008), if a is assumed to be associated with the particular SED template being used to fit the stellar photometry. These estimates of and are presented in the “XO-Rad” database at the end of this paper.

We describe the observations and data reduction of these stars in §4; supporting data and spectral energy distribution fits are described in §3; derived effective temperatures and radii are presented in §5, along with comparisons of our values to previous investigations (where available); finally, a detailed statical comparison of the EHS stars versus our control group is seen in §5.3.

2. Description of the Datasets

We present interferometric results on two different datasets:

  1. Known EHSs for which we were able to obtain PTI data (§4) and calculate angular radii (§5). Knowledge of angular radii imposes an independent constraint on the SED fitting (§3) and allows to be a determined directly. We were also able to augment our PTI data with 7 stars published from the CHARA Array by Baines et al. (2008a). This dataset comprises 12 stars, 4 of which have data from both CHARA and PTI. Together this sample of EHSs with angular sizes is our ‘EHSA’ sample.

  2. A number of main-sequence stars for which it was deemed possible to resolve angular radii using PTI. This dataset comprises 28 stars and will be referred to as our ‘control sample’. These stars are not currently known to host extrasolar planets and thus serve as a comparison group for the EHSs with respect to astrophysical parameters.

Additionally, SED fits are provided for all the well-characterized EHSs (status 1 Feb 2008, according to the Exoplanet Encyclopedia333http://exoplanet.eu/). The source dataset comprises approximately 230 stars (including the ones for which we obtained PTI data), although most of the fainter (10) stars are excluded due to a lack of available photometry; they are presented in the “XO-Rad” database in Appendix A. SED fitting for these stars is performed based on literature photometry and spectral templates with associated estimates of effective temperatures.

3. Supporting Data and Spectral Energy Distribution Fitting

For all of the sources considered in this investigation, spectral energy distribution (SED) fits were performed. Each fit, accomplished using available photometry and an appropriate template spectrum, produces estimates for the bolometric flux (), the angular diameter () and the reddening (); effective temperature during the SED fit is fixed for each of the template spectra. In the absence of direct measurement of the angular diameter (i.e. calibrators and stars listed in the XO-Rad database), SED fitting is used to estimate the angular size. When the angular diameter is available from interferometrically measurements, SED fitting is used to determine the bolometric flux and the reddening; effective temperature as well as dereddened colors can then be derived.

These SED fits are accomplished using photometry available in the literature as the input values, with template spectra from the Pickles (1998) library appropriate for the spectral types indicated for the stars in question. Spectral types used in the SED fitting for all EHS stars are those values found in the Exoplanet Encyclopedia, which is in turn based upon the respective source discovery papers cataloged therein. The control sample stars as defined in §2 had their spectral types established from those values found in Hipparcos catalog (Perryman et al., 1997).

The template spectra are adjusted by the fitting routine to account for overall flux level, wavelength-dependent reddening, and expected angular size. Reddening corrections are based upon the empirical reddening determination described by Cardelli et al. (1989), which differs little from van de Hulst’s theoretical reddening curve number 15 (Johnson, 1968; Dyck et al., 1996). Both narrowband and wideband photometry in the 0.3 m to 30 m are used as available, including Johnson (see, for example, Eggen, 1963; Moreno, 1971) Stromgren (Piirola, 1976), 2MASS (Cutri et al., 2003), Geneva (Rufener, 1976), and Vilnius (Zdanavicius et al., 1972); flux calibrations are based upon the values given in Fukugita et al. (1995) and Cox (2000). The results of the fitting for the calibrator stars is given in Table 1; for the EHSA and control sample stars, Table 6, and for the “XO-Rad” database, Table 7.

4. Observations and Data Reduction

4.1. Visibility and Angular Sizes

The calibration of the target star visibility data is performed by estimating the interferometer system visibility () using the calibration sources with model angular diameters and then normalizing the raw target star visibility by to estimate the measured by an ideal interferometer at that epoch (Mozurkewich et al., 1991; Boden et al., 1998; van Belle & van Belle, 2005). Uncertainties in the system visibility and the calibrated target visibility are inferred from internal scatter among the data in an observation using standard error-propagation calculations (Boden et al., 1999). Calibrating our point-like calibration objects against each other produced no evidence of systematics, with all objects delivering reduced .

Visibility and uniform disk angular size are related using the first Bessel function : , where spatial frequency . We may establish uniform disk angular sizes for the target stars observed by the interferometer since the accompanying parameters (projected telescope-to-telescope separation, or baseline, and wavelength of observation ) are well-characterized during the observation. The uniform disk angular size can (and should) be connected to a more physical limb darkened angular size ; however, this is a minor effect since is small in the near-infrared (; see, for example, Scholz & Takeda, 1987; Tuthill, 1994; Dyck et al., 1996, 1998; Davis et al., 2000).

Strictly speaking, limb darkened angular size is utilized here as a reasonable proxy for the Rosseland angular diameter, which corresponds to the surface where the Rosseland mean optical depth equals unity, as advocated by Scholz & Takeda (1987) as the most appropriate surface for computing an effective temperature. The dense, compact atmospheres of the stars considered in this investigation are well characterized by a uniform disk fit, and the small correction factors tabulated in Davis et al. (2000) will be used to convert our sizes into the appropriate limb darkened numbers. The number of visibility points , derived sizes, associated goodness-of-fit and residuals , Davis et al. (2000) correction factors and resultant sizes are found in the first columns of Table 6.

4.2. PTI Observations

PTI is an 85 to 110 m - and -band 1.6 m and 2.2 m) interferometer located at Palomar Observatory in San Diego County, California, and is described in detail in Colavita (1999). It has three 40-cm apertures used in pairwise combination for detection of stellar fringe visibility on sources that range in angular size up to 5.0 milliarcseconds (mas), being able to resolve individual sources with angular diameter greater than 0.60 mas in size. PTI has been in nightly operation since 1997, with minimum downtime throughout the intervening years. The data from PTI considered herein cover the range from the beginning of 1998 (when the standardized data collection and pipeline reduction went into place) until the beginning of 2008 (when the analysis of this manuscript was begun). In addition to the target stars discussed herein, appropriate calibration sources were observed as well and can be found en masse in van Belle et al. (2008). Additional calibration sources of minimal angular size, as discussed in §4.3, were also selected and are listed in Table 1.

Star RA DE SpType Model
HD4058 00 43 28.09 +47 01 28.7 75 A5V A5V 1.07
HD10205 01 40 34.80 +40 34 37.6 55 B8III B5III 1.83
HD10874 01 47 48.00 +46 13 47.6 21 F6V F6V 3.81
HD11529 01 56 00.00 +68 41 07.0 39 B8III B5III 1.06
HD13476 02 13 41.61 +58 33 38.1 78 A3Iab A2I 4.02
HD14212 02 19 16.85 +47 22 48.0 37 A1V A0V 1.15
HD15138 02 27 51.75 +50 34 12.7 26 F4V F2V 2.80
HD16399 02 38 00.70 +07 41 43.4 66 F6IV F5IV 0.27
HD16582 02 39 28.95 +00 19 42.7 79 B2IV B2IV 2.13
HD17163 02 45 20.87 +04 42 42.2 66 F0III: F0III 0.58
HD18331 02 56 37.45 -03 42 44.0 244 A3Vn A3V 3.53
HD20418 03 19 07.62 +50 05 42.1 33 B5V B57V 2.43
HD23005 03 46 00.82 +67 12 06.8 28 F0IV F02IV 0.24
HD23363 03 44 30.51 -01 09 47.1 48 B7V B57V 0.91
HD24479 03 57 25.44 +63 04 20.1 50 B9.5V B9V 1.79
HD35039 05 21 45.75 +00 22 56.9 104 B2IV-V B2IV 3.59
HD36777 05 34 16.79 +03 46 01.0 62 A2V A2V 4.54
HD37077 05 35 39.49 -04 51 21.9 51 F0III F0III 0.94
HD41040 06 03 27.36 +19 41 26.2 63 B8V B8V 3.28
HD42618 06 12 00.45 +06 47 01.3 60 G4V G2V 2.30
HD46300 06 32 54.23 +07 19 58.7 114 A0Ib A0I 3.85
HD86360 09 58 13.39 +12 26 41.4 43 B9IV B6IV 3.43
HD89389 10 20 14.88 +53 46 45.4 36 F9V F8V 0.50
HD91480 10 35 09.62 +57 04 57.2 98 F1V F0V 0.43
HD93702 10 49 15.43 +10 32 42.9 54 A2V A2V 4.08
HD96738 11 08 49.08 +24 39 30.4 33 A3IV A0IV 1.71
HD97334 11 12 32.53 +35 48 52.0 61 G0V G0V 0.24
HD97486 11 14 04.63 +62 16 55.7 15 G5III G5III 1.33
HD102634 11 49 01.40 +00 19 07.2 70 F7V F6V 0.71
HD103578 11 55 40.53 +15 38 48.5 61 A3V A3V 3.09
HD104181 11 59 56.92 +03 39 18.8 60 A1V A0V 1.72
HD106661 12 16 00.23 +14 53 56.9 68 A3V A3V 1.35
HD110392 12 41 26.98 +40 34 45.7 15 K0III K0III 5.20
HD111604 12 50 10.81 +37 31 00.8 50 A3V A3V 2.62
HD113771 13 05 40.89 +26 35 08.5 11 K0III K0III 3.69
HD114762 13 12 19.743 +17 31 01.6 100 F9V F8V 1.88
HD119288 13 42 12.98 +08 23 19.0 46 F3Vp F2V 2.29
HD119550 13 43 35.700 +14 21 56.1 48 G2V G2V 0.22
HD119550 13 43 35.700 +14 21 56.1 48 G2V G2V 0.22
HD119550 13 43 35.89 +14 21 56.3 48 G2V G2V 0.22
HD121560 13 55 49.994 +14 03 23.4 36 F6V F6V 0.79
HD125161 14 16 10.07 +51 22 01.3 37 A9V A7V 1.41
HD128332 14 34 15.70 +57 03 57.0 29 F7V F6V 1.19
HD140775 15 45 23.47 +05 26 50.4 101 A1V A0V 0.92
HD141187 15 47 17.35 +14 06 55.0 37 A3V A3V 5.94
HD142908 15 55 47.587 +37 56 49.0 123 F0IV F02IV 0.24
HD144579 16 04 56.793 +39 09 23.4 70 G8V G8V 2.09
HD144874 16 07 37.55 +09 53 30.3 46 A7V A7V 0.91
HD150557 16 41 42.54 +01 10 52.0 62 F2.7III-IV F2III 2.64
HD151900 16 50 22.25 -02 39 15.3 58 F1III-IV F0III 4.87
HD158352 17 28 49.69 +00 19 50.1 42 A8V A7V 1.18
HD164353 18 00 38.72 +02 55 53.7 110 B5Ib B5I 1.53
HD164613 17 55 11.14 +72 00 18.5 28 F2.5II-III F2III 3.31
HD169702 18 24 13.80 +39 30 26.1 25 A3IVn A0IV 1.64
HD173649 18 44 48.19 +37 35 40.4 76 F0IVvar F02IV 0.90
HD180482 19 16 31.02 +04 50 05.4 46 A3IV A0IV 2.21
HD180777 19 09 09.75 +76 33 38.9 67 A9V A7V 1.12
HD182564 19 20 40.07 +65 42 51.9 38 A2IIIs A0III 1.81
HD184663 19 35 25.13 +02 54 48.5 36 F6IV F5IV 1.14
HD186568 19 43 51.452 +34 09 45.8 58 B8III B5III 3.94
HD192640 20 14 32.033 +36 48 22.6 86 A2V A2V 2.44
HD198478 20 48 56.29 +46 06 50.9 79 B3Ia B3I 0.62
HD199081 20 53 14.75 +44 23 14.2 71 B5V B57V 2.76
HD200723 21 03 52.14 +41 37 41.9 21 F3IV F02IV 0.64
HD202240 21 13 26.43 +36 37 59.7 53 F0III F0III 4.96
HD210264 22 08 50.40 +22 08 19.6 15 G5III G5III 1.33
HD214734 22 38 39.05 +63 35 04.3 34 A3IV A0IV 2.21
HD217813 23 03 04.977 +20 55 06.8 40 G5V G5V 5.09
HD218261 23 06 31.71 +19 54 39.0 42 F7V F6V 1.96
HD218261 23 06 31.885 +19 54 39.0 42 F7V F6V 1.96
HD218396 23 07 28.715 +21 08 03.3 82 A5V A5V 3.21
HD218687 23 09 57.17 +14 25 36.3 30 G0V G0V 0.71
HD220102 23 20 20.82 +60 16 29.2 40 F5II F2II 1.07
HD220102 23 20 20.82 +60 16 29.2 40 F5II F2II 1.07
HD223346 23 48 49.36 +02 12 52.2 64 F5III-IV F5III 1.33

Note. – is the number of photometric data points available for the bolometric flux fitting; SpType is the spectral type as reported by SIMBAD; Model is the spectral template chosen from Pickles (1998) for the fitting; is the reduced chi-squared value of the fit, and is the estimated angular size from the fit.

Table 1New calibration sources used in this investigation, discussed in §4.2.

4.3. Limits of PTI Calibration

As discussed by Boden et al. (1998, 1999), PTI has an empirically established fundamental limiting visibility measurement error of %. The source of this limiting night-to-night measurement error is most likely a combination of effects: uncharacterized atmospheric seeing (in particular, scintillation), detector noise, and other instrumental effects.

This night-to-night repeatability limit restricts the ultimate resolution of the instrument. This is at odds with the desire to measure stellar diameters which, for a given brightness, are quite small in an angular sense relative to PTI’s resolution. Main sequence stars are squarely in this regime for PTI, with only a few examples - those considered in this investigation - that creep out of the nether regions of point-like obscurity into the realm of resolvability. Attempting to resolve stars at the edge of PTI’s performance envelope requires careful consideration of the demonstrated limits of the instrument, using the techniques described in van Belle & van Belle (2005, henceforth Paper VB2).

For PTI, operating at the -band with its 109-m N-S baseline, a target of 0.60 milliarcseconds (mas) in size should have a normalized visibility of (as introduced in §4.1). As discussed in in VB2, there is a strong motivation towards using calibration sources that are as point-like as possible - generally speaking, one wishes to have calibration sources that are significantly smaller than the targets being observed. For this investigation, to reach the regime of 0.60 mas targets, we restricted our use of calibrators to those that are, on average, 0.35 mas or less in size. These two size limits are selected to have sufficient numbers of sufficiently bright targets and calibrators, respectively.

For such calibrators, observed by PTI, the visibility calibration limit is (from VB2, Equation 7), which contributes an angular size error due to calibration of roughly 0.012 mas. The night-to-night limiting measurement error of , however, contributes an angular size error of 0.086 mas. This is significant in that the measurement error dominates any possible calibration bias, which is particularly important when considering smaller targets. If we were instead to have selected calibrators closer to mas in size - more typical of PTI investigations that observe larger targets that are mas in size - then the calibration angular size error be mas, and would start to compete with the measurement error in dominating the error budget. This would put our results at substantial risk of directly reporting any measurement bias inherent in the process we used to estimate the angular sizes of our calibration sources. Since our goal is direct measurement of the target angular sizes, we have taken great care to ensure that this is not the case.

A second aspect of this consideration of PTI limiting performance is the reported angular sizes of our target stars. For stars that, after calibration, report formal errors that are sufficiently small to be in violation of PTI’s known night-to-night repeatability, we increased their reported angular size errors to the level consistent with that repeatability. As a function of target angular size, we show the limits of angular size accuracy possible with PTI’s repeatability limit in Table 2. The first column shows various target angular sizes, followed by the corresponding visibilities. A calibrator of 0.35 mas, as noted above, contributes a the limit on knowledge of visibility of ; the associated limit in angular size knowledge is then listed in column three. The next two columns list the night-to-night repeatability limit of , and the associated angular size error. The final column combines the calibration limit and the night-to-night limit in quadrature.

Calibration Calibration Night-to-night Night-to-night
Target Target floor
(mas) (mas) (mas) (mas)
0.600 0.94893 0.00186 0.012 0.01500 0.085 0.086
0.650 0.94028 0.00186 0.010 0.01500 0.079 0.080
0.700 0.93103 0.00186 0.010 0.01500 0.075 0.076
0.750 0.92116 0.00186 0.010 0.01500 0.071 0.072
0.800 0.91072 0.00186 0.009 0.01500 0.068 0.069
0.850 0.89971 0.00186 0.008 0.01500 0.064 0.065
0.900 0.88815 0.00186 0.008 0.01500 0.062 0.063
0.950 0.87607 0.00186 0.008 0.01500 0.060 0.060
1.000 0.86348 0.00186 0.008 0.01500 0.058 0.058
Table 2Calibration floor by target angular size as discussed in §4.3.

4.4. CHARA EHS Data

Additional angular diameters of EHSs were obtained with the Georgia State University CHARA Array (Baines et al., 2008a) with an intent of detecting possible face-on binarity masquerading as planetary companionship (Baines et al., 2008b). The CHARA Array is a optical/near-infrared interferometer similar to PTI (ten Brummelaar et al., 2005), but with longer baselines (up to 330m), allowing for resolution of smaller objects. For inclusion of the appropriate CHARA data into our dataset, we will apply observation criteria similar to the PTI data: First, the calibration sources must be sufficiently unresolved, which we set for CHARA to be 0.50mas or less. Second, the ratio of angular sizes of science targets and their calibrators must be greater than 1.5. In applying these two criteria, we are confident that the resulting measured angular sizes are sufficiently independent of the calibrator angular sizes predicted by SED fitting.

The resulting dataset for inclusion in this analysis consists of seven EHS angular sizes from the CHARA investigation, of which four stars are common to both the PTI and CHARA samples (as noted §2). The ratios of the CHARA to PTI UD angular sizes for those four stars (HD3651, HD75732, HD143761, HD217014) are , respectively, with an overall weighted average ratio of , indicating possibly a slight tendency for the PTI sizes to be too small (or the CHARA sizes to be too large), but this is a weak 1- result.

As a further check on the consistency of the CHARA results and our techniques, we modeled the predicted SED sizes of the calibrators found in Baines et al. (2008a). These results are seen in Table 3; on average, our calibrator predictions are within , and no individual results are more than away from Baines et al. (2008a). Overall, we find that the CHARA and PTI results are excellent agreement with each other, despite independently developed methodologies.

Target Calibrator Calibrator Size CHARA Difference
HD HD Est. (mas) Est. (mas) (mas)
3651 4568 -0.016 1.6
11964 13456 -0.027 1.9
19994 19411 0.001 0.0
75732 72779 -0.002 0.1
143761 136849 0.019 -0.5
189733 190993 0.001 0.0
217014 218261 -0.003 0.2
Table 3Comparison of spectral energy distribution fits for calibrators from Baines et al. (2008a) as discussed in §4.4.

5. Stellar Parameters

For both the EHSA and our control sample stars, the basic astrophysical parameters of effective temperature and linear radius are computed from the angular size data and ancillary supporting data. These parameters are then compared between the two samples as a function of color and, in the case of temperature, spectral type; the results of sections §5.1 and §5.2 are found in Table 4.

5.1. Effective Temperatures

Stellar effective temperature, , is defined in terms of the star’s luminosity and radius by . As noted in §1, rewriting this equation in terms of angular diameter and bolometric flux , can be expressed as , where is in ergs cm s and is in mas (van Belle et al., 1999). The derived temperature values for the resolved stars of this study are found in Table 4, along with color. These temperatures are plotted versus in Figure 1, and to explore any potential difference between the EHSA stars and the control sample, a fit of the versus trend is performed.

For the control sample, the initial fit reveals HD87901 as a significant outlier. This is most likely due to two factors: (1) HD87901 is bluest and hottest star, at and K, and (2) HD87901 is a rapid rotator with km/s (Abt et al., 2002), and will show departures from sphericity that induce gravity darkening which render individual determinations meaningless (Aufdenberg et al., 2006). Omitting HD87901 from the fit, the best fit for the control sample stars is


with =1.72, with the fitting and error ellipses following the techniques described in Press et al. (1992). (Inclusion of HD87901 in this fit returns =4.98.)

If we include the EHSA stars in the fit, we find the CHARA data point for 55 Cnc (HD75732) a significant outlier as well, which we will discuss further in §5.4.1. Omitting 55 Cnc from the unified fit, we find a single fit gives:


with =1.82. This fit line is plotted in Figure 1. These fits indicate there is no statistically significant difference between the two populations (noting that the EHSA fit is poorly constrained with a small number of data points over a small range of , preventing a fit to those data alone). We revisit the question of population similarity in further detail in §5.3.

For the fit in Equation 2, the median value of the differences between the values predicted by this fit and the measured values is K. Since the median value of the errors in the individual measurements is K, we believe the limit of precision in the line fit is not due to any intrinsic astrophysical scatter in the versus relationship, but rather the limits of the current measurements.

Alternatively, a fit may be made for a cubic relationship between and , (see, for example, the corresponding equation in Levesque et al., 2005) but this produces no significant improvement:


with only =1.68, in spite of the extra degree of freedom.

For those spectral types for which we have more than one stellar angular size measurement, we can compare the resultant weighted mean values to the ‘canonical’ values cited in Cox (2000), which can be traced back to the investigation by de Jager & Nieuwenhuijzen (1987). This comparison is seen in Table 5. It is interesting to note that our values of all track increasingly lower between types F8V to G2V in comparison to the de Jager & Nieuwenhuijzen (1987) values, before returning to agreement with those values at G5V and cooler.

Finally, given the large number of individual samples of our data set between types F6V and G5V, we present an empirical calibration of versus spectral type for this full range, also in Table 5. Spectral types that have no measurements (e.g., F7V) have values interpolated from the adjoining spectral types. The average error by spectral type is K. This table and Equation 2 represent a direct calibration of the scale for solar-like main sequence stars for the spectral type range F6V-G5V and color range . No attempt was made for calibration for the later types due to the sparseness of the data, although our data at K1V, K7V and M2V represent calibration for those specific spectral types.

Star ID
(mag) (K) (pc)
Control Sample: Stars not known to host planets:
EHSA Sample: Known planet hosting stars (PTI):
EHSA Sample: Known planet hosting stars (CHARA):
Table 4Dereddened colors, effective temperatures and radii for luminosity class V stars, discussed in §5.
Figure 1.— Effective temperature versus color for control sample and EHSA stars. Also shown is a fit to the luminosity class V stars (solid line, discussed in §5.1), the relationship for giants found in van Belle et al. (1999) (dashed line) and for a blackbody radiator (dotted line). The median deviation of the stellar data points from the solid line fit is K.

5.2. Linear Radii

From the parallax values found in Table 6 from Hipparcos (Perryman et al., 1997), linear radii are derived for the resolved stars of this investigation and are found in Table 4. A cubic relationship fit to the combined EHSA and control samples is:


with a =15.1. Clearly this metric indicates a poor fit, which is consistent with some of the stars beginning to evolve well off of the zero-age-main-sequence (ZAMS) line. This effect is seen in a plot of the data in Figure 2, with the presumably older stars being situated to the right of the line fit. As such, Equation 4 should be regarded as only a rough indication of stellar radius, and not applicable in any general sense to determining linear radii of random field stars.

Spectral N
Type (K) (K)
F6V 6 6515aaNo specific value given in Cox (2000), interpolated from neighboring data points.
F7V 6385aaNo specific value given in Cox (2000), interpolated from neighboring data points.
F8V 4 6250
F9V 6095aaNo specific value given in Cox (2000), interpolated from neighboring data points.
G0V 7 5940
G1V 5865aaNo specific value given in Cox (2000), interpolated from neighboring data points.
G2V 2 5790
G3V 5715aaNo specific value given in Cox (2000), interpolated from neighboring data points.
G4V 5635aaNo specific value given in Cox (2000), interpolated from neighboring data points.
G5V 4 5560
K1V 4 4990aaNo specific value given in Cox (2000), interpolated from neighboring data points.
K7V 3 4125aaNo specific value given in Cox (2000), interpolated from neighboring data points.
M2V 3 3520

Note. – See discussion at end of §5.1. Data after G5V were sufficiently sparse to not merit empirical calibration of the full range. is the number of angular size measurements per spectral type; rows with no value for are interpolated values. Columns 3 and 4 are from this work and Cox (2000), respectively.

Table 5Effective temperature versus spectral type, with an empirical calibration of effective temperature versus spectral type for types F6V through G5V.
Figure 2.— Linear radius versus color for control sample and EHSA stars as discussed in §5.2. Also shown is a fit to the control sample and EHSA stars (solid line). One of our EHSA stars, HD28305, is a giant star with and is off the scale of this plot.

5.3. Kolmogorov-Smirnov Comparison Between Exoplanet Hosting Stars and Control Stars

As detailed in Press et al. (1992), the Kolmogorov-Smirnov (KS) test can be executed to compare two arrays of data values, and examine the probability that the two arrays are drawn from the same distribution. The KS test returns two values: the KS statistic , which specifies the maximum deviation between the cumulative distribution of the two sample of data, and probability , giving the significance of the KS statistic. Small values of show that the two distributions differ significantly.

Examining the versus data of the EHSA stars versus the control sample stars, we find that with – strong indication that two data sets are indeed statistically indistinguishable. The astrophysical implication is that, within the limits of our measurements, the effective temperature scale of stars with known planets does not differ from those without known planets.

The corresponding versus KS test, however, reports and , which seems to indicate the two samples are inconsistent with each other. However, the significance of this result is simple: our control sample is specifically selected to be main sequence stars, whereas the EHSA sample includes a number of evolved sources, as clearly seen in Figure 2. One corollary implication of these two KS tests is that stars on main sequence and those evolving off of it do not differ significantly in their versus relationships.

5.4. Comparison with Previous Studies

There is a variety of data available for the known EHSs in the literature, derived from different methods by different authors. Thus, discrepancies, though sometimes small, exist. In order to be as consistent as possible, we chose the following two catalogs as data sources for astrophysical parameters:

  • Mass, age, and [Fe/H] from Valenti & Fischer (2005)

  • Linear radius from Takeda et al. (2007)

A comparison of the values measured in this investigation can be directly contrasted against those found in Valenti & Fischer (2005). Combining our EHSA and control star samples, we find:


with . As illustrated in Figure 3, there is no significant difference between the values obtained with interferometry and spectroscopy.

A marginal offset is found between our and the radii of Takeda et al. (2007):


with – roughly a 2- offset between the line slope and intercept values for from theory versus those determined interferometrically. The general trend is for the larger stars to have a larger theoretical, rather than interferometric, linear size. These values and the general trend can be seen in Figure 4.

Figure 3.— Effective temperature as determined by this study, versus those values found spectroscopically by Valenti & Fischer (2005) for EHSA stars (red triangles) and our control sample stars (blue diamonds), as discussed in §5.3. The solid line is the 1:1 line, with the dotted line the fit to the versus data points.
Figure 4.— Linear radii as determined by this study, versus those values found spectroscopically by Takeda et al. (2007) for EHS stars (red triangles) and our control group (blue diamonds). The solid line is the 1:1 line, with the dotted line the fit to the values. A trend is seen with the larger stars being larger in the Takeda et al. (2007) study.

5.4.1 Discussion of 55 Cnc (HD 75732)

Inclusion of the PTI and CHARA data points for 55 Cnc44455 Cnc’s distance is pc (Perryman et al., 1997). It is K0IV–V star (Gray et al., 2003) with (Bessell, 2000). It has a mass of , an age of Gyr, K, and [Fe/H]= (Valenti & Fischer, 2005). Its radius is 0.91 in Pasinetti Fracassini et al. (2001) and when using the equations in Lang (1980). The values from this investigation are: , . in the fit of Equation 2 pushes the from 1.82 up to 2.91, with the CHARA data points remaining as outliers; inclusion of just the PTI points results in =1.88. As such, we decided to omit both the CHARA and PTI data points for 55 Cnc from the fit. There are two possible reasons for 55 Cnc turning up as “too cool” to fall onto the versus fit of Equation 2.

First, the CHARA data points could be in error: including just the PTI data for 55 Cnc does not significantly alter the resulting value. However, the angular size and values for 55 Cnc from PTI and CHARA are in direct agreement with each other, although the PTI size data point has a larger error, indicative of its lesser resolving power for this mas star. To ‘force’ the 55 Cnc data onto the versus fit line, its angular size would need to be reduced to mas. Calibrator size error does not appear to be the source of the problem: the size of calibrator HD72779 quoted in Baines et al. (2008a) is mas – confirmed independently in this investigation with a value of mas – and would have to be mas to push the 55 Cnc visibility data to deliver the larger angular size. Alternatively, the calculation for 55 Cnc could be too low, but require an increase from erg cm s to erg cm s, which is far outside the allowable bounds of SED fitting, regardless of the template selected.

The second possible reason is that the visibility data could be contaminated by the presence of a secondary stellar companion. Such a companion would reduce the observed visibility, resulting in an apparent increase in angular size, which in turn would effect an apparent decrease in derived temperature - as seen with the 55 Cnc data. Examination of the plots associated with the CHARA dates and configuration cited in Baines et al. (2008a) indicate a small amount () of baseline rotation, with nearly zero change in baseline length, which would have lead to a null result in detection in Baines et al. (2008b) for a secondary stellar companion - even in some cases where one is present. There is a known companion to 55 Cnc at a distance of 1000 pc, or 95 on the sky; however, with (based on a spectral type of M4), in the worst case we would see a visibility change of of only , which would only lower the apparent size from 0.85 to 0.82 mas. Additionally, our naïve expectation is that the intensive spectroscopic studies of 55 Cnc that have turned up no less than 5 planets (Fischer et al., 2008) would have uncovered such a companion, so we are at a loss as to how to reconcile interferometric data with the spectroscopic discoveries. For the moment we will be content to simply remove it from the effective temperature scale calibrations presented in §5.1.

6. Summary and Conclusion

We present directly determined stellar radii and effective temperatures for 12 exoplanet host stars, along with the same estimates for 28 main-sequence control stars not known to host planets. In the process, we demonstrate the empirical limit of PTI’s stellar angular resolution and the implications for angular sizes measured near that limit. While our results show consistency between the direct measurements of effective temperature and indirectly determined literature values, a small difference exists between our radii measurements and theoretical estimates in the sense that for larger stars, the theoretical estimate falls slightly above the direct measurement. From our effective temperature measurements, an empirical calibration of effective temperature versus color and spectral type is presented, with a spread of K over the range and K for F6V – G5V. No such calibration is possible for linear radius versus color, due to the large spread in radius values for any given color (presumably due to stellar evolution effects). Among the stars considered, 55 Cnc is found to be problematic in terms of its interferometrically determined effective temperature, for reasons that are unclear. Finally, the spectral energy distribution fitting tools employed in this investigation also enable indirect estimates of stellar angular size to be attempted for the full ensemble of stars known to host extrasolar planets, and this database of 166 stars is presented in the “XO-Rad” appendix.

We would like to acknowledge constructive input and the occasional snide comment from David Ciardi. This investigation has made extensive use of the sedFit code, graciously provided by perl guru Andrew F. Boden. The preparation of this manuscript was greatly helped by the use of the Extrasolar Planet Encyclopedia555Available at http://exoplanet.eu.. This research made use of the NASA/IPAC/NExScI Star and Exoplanet Database (NStED)666Available at http://nsted.ipac.caltech.edu., which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. The Palomar Testbed Interferometer is operated by the NASA Exoplanet Science Institute/Michelson Science Center on and the PTI collaboration and was constructed with funds from the Jet Propulsion Laboratory, Caltech as provided by the National Aeronautics and Space Administration. This work has made use of services produced by the NASA Exoplanet Science Institute at the California Institute of Technology.

Appendix A The XO-Rad Database

For the full list of stars found at the Extrasolar Planet Encyclopedia (as of 1 Feb 2008), we collected photometry and performed SED fits as described in the main manuscript in §3, and in detail in van Belle et al. (2008). 64 of the stars have insufficient photometry and were dropped from the SED fitting. The resultant 166 fits provide estimates of bolometric flux , -band reddening , angular size , and linear radius . Effective temperatures are constrained to be those associated with the best fitting Pickles (1998) empirical template. Spectral types used in the SED fitting for EHS stars are those values found in the Exoplanet Encyclopedia, which is in turn based upon the respective source discovery papers cataloged therein. The non-planet-hosting main sequence stars have their spectral types established from those values found in Hipparcos catalog (Perryman et al., 1997). Linear radius is computed by combining the angular size estimates with the Hipparcos data found in van Leeuwen (2007). For a few of the stars, the linear radius is too large to be consistent with the main sequence spectral types indicated in the literature; for these objects, a second iteration on the SED fit is performed with a subgiant (luminosity class IV) template, resulting in a more appropriate set of fit parameters . The full XO-Rad dataset of exoplanet radii is seen in Table 7.


  • Abt et al. (2002) Abt, H. A., Levato, H., & Grosso, M. 2002, ApJ, 573, 359
  • Aufdenberg et al. (2006) Aufdenberg, J. P., Mérand, A., Coudé du Foresto, V., Absil, O., Di Folco, E., Kervella, P., Ridgway, S. T., Berger, D. H., ten Brummelaar, T. A., McAlister, H. A., Sturmann, J., Sturmann, L., & Turner, N. H. 2006, ApJ, 645, 664
  • Baines et al. (2008a) Baines, E. K., McAlister, H. A., ten Brummelaar, T. A., Turner, N. H., Sturmann, J., Sturmann, L., Goldfinger, P. J., & Ridgway, S. T. 2008a, ApJ, 680, 728
  • Baines et al. (2008b) Baines, E. K., McAlister, H. A., ten Brummelaar, T. A., Turner, N. H., Sturmann, J., Sturmann, L., & Ridgway, S. T. 2008b, ApJ, 682, 577
  • Baines et al. (2007) Baines, E. K., van Belle, G. T., ten Brummelaar, T. A., McAlister, H. A., Swain, M., Turner, N. H., Sturmann, L., & Sturmann, J. 2007, ApJ, 661, L195
  • Bessell (2000) Bessell, M. S. 2000, PASP, 112, 961
  • Boden et al. (1998) Boden, A. F., Colavita, M. M., van Belle, G. T., & Shao, M. 1998, in Presented at the Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 3350, Proc. SPIE Vol. 3350, p. 872-880, Astronomical Interferometry, Robert D. Reasenberg; Ed., ed. R. D. Reasenberg, 872–880
  • Boden et al. (1999) Boden, A. F., Koresko, C. D., van Belle, G. T., Colavita, M. M., Dumont, P. J., Gubler, J., Kulkarni, S. R., Lane, B. F., Mobley, D., Shao, M., Wallace, J. K., The PTI Collaboration, & Henry, G. W. 1999, ApJ, 515, 356
  • Boltzmann (1884) Boltzmann, L. 1884, Annalen der Physik und Chemie, 22, 291
  • Cardelli et al. (1989) Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245
  • Colavita (1999) Colavita, M. M. 1999, PASP, 111, 111
  • Cox (2000) Cox, A. N. 2000, Allen’s astrophysical quantities (Allen’s astrophysical quantities, 4th ed. Publisher: New York: AIP Press; Springer, 2000. Editedy by Arthur N. Cox.  ISBN: 0387987460)
  • Cutri et al. (2003) Cutri, R. M., Skrutskie, M. F., van Dyk, S., Beichman, C. A., Carpenter, J. M., Chester, T., Cambresy, L., Evans, T., Fowler, J., Gizis, J., Howard, E., Huchra, J., Jarrett, T., Kopan, E. L., Kirkpatrick, J. D., Light, R. M., Marsh, K. A., McCallon, H., Schneider, S., Stiening, R., Sykes, M., Weinberg, M., Wheaton, W. A., Wheelock, S., & Zacharias, N. 2003, 2MASS All Sky Catalog of point sources. (The IRSA 2MASS All-Sky Point Source Catalog, NASA/IPAC Infrared Science Archive. http://irsa.ipac.caltech.edu/applications/Gator/)
  • Davis et al. (2000) Davis, J., Tango, W. J., & Booth, A. J. 2000, MNRAS, 318, 387
  • de Jager & Nieuwenhuijzen (1987) de Jager, C. & Nieuwenhuijzen, H. 1987, A&A, 177, 217
  • Dyck et al. (1996) Dyck, H. M., Benson, J. A., van Belle, G. T., & Ridgway, S. T. 1996, AJ, 111, 1705
  • Dyck et al. (1998) Dyck, H. M., van Belle, G. T., & Thompson, R. R. 1998, AJ, 116, 981
  • Eggen (1963) Eggen, O. J. 1963, AJ, 68, 483
  • Fischer et al. (2008) Fischer, D. A., Marcy, G. W., Butler, R. P., Vogt, S. S., Laughlin, G., Henry, G. W., Abouav, D., Peek, K. M. G., Wright, J. T., Johnson, J. A., McCarthy, C., & Isaacson, H. 2008, ApJ, 675, 790
  • Fukugita et al. (1995) Fukugita, M., Shimasaku, K., & Ichikawa, T. 1995, PASP, 107, 945
  • Gray et al. (2003) Gray, R. O., Corbally, C. J., Garrison, R. F., McFadden, M. T., & Robinson, P. E. 2003, AJ, 126, 2048
  • Johnson (1968) Johnson, H. L. Interstellar Extinction, ed. B. M. Middlehurst (University of Chicago Press, Chicago, IL USA), 167
  • Lang (1980) Lang, K. R. 1980, Astrophysical Formulae (A Compendium for the Physicist and Astrophysicist, XXIX, 783 pp. 46 figs., 69 tabs..  Springer-Verlag Berlin Heidelberg New York.  Also Springer Study Edition)
  • Levesque et al. (2005) Levesque, E. M., Massey, P., Olsen, K. A. G., Plez, B., Josselin, E., Maeder, A., & Meynet, G. 2005, ApJ, 628, 973
  • Moreno (1971) Moreno, H. 1971, A&A, 12, 442
  • Mozurkewich et al. (2003) Mozurkewich, D., Armstrong, J. T., Hindsley, R. B., Quirrenbach, A., Hummel, C. A., Hutter, D. J., Johnston, K. J., Hajian, A. R., Elias, II, N. M., Buscher, D. F., & Simon, R. S. 2003, AJ, 126, 2502
  • Mozurkewich et al. (1991) Mozurkewich, D., Johnston, K. J., Simon, R. S., Bowers, P. F., Gaume, R., Hutter, D. J., Colavita, M. M., Shao, M., & Pan, X. P. 1991, AJ, 101, 2207
  • Pasinetti Fracassini et al. (2001) Pasinetti Fracassini, L. E., Pastori, L., Covino, S., & Pozzi, A. 2001, A&A, 367, 521
  • Perryman & ESA (1997) Perryman, M. A. C. & ESA, eds. 1997, ESA Special Publication, Vol. 1200, The HIPPARCOS and TYCHO catalogues. Astrometric and photometric star catalogues derived from the ESA HIPPARCOS Space Astrometry Mission
  • Perryman et al. (1997) Perryman, M. A. C., Lindegren, L., Kovalevsky, J., Hoeg, E., Bastian, U., Bernacca, P. L., Crézé, M., Donati, F., Grenon, M., van Leeuwen, F., van der Marel, H., Mignard, F., Murray, C. A., Le Poole, R. S., Schrijver, H., Turon, C., Arenou, F., Froeschlé, M., & Petersen, C. S. 1997, A&A, 323, L49
  • Pickles (1998) Pickles, A. J. 1998, PASP, 110, 863
  • Piirola (1976) Piirola, V. 1976, Observatory and Astrophysics Laboratory University of Helsinki Report, 1, 0
  • Press et al. (1992) Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in C. The art of scientific computing (Cambridge: University Press, —c1992, 2nd ed.)
  • Rufener (1976) Rufener, F. 1976, A&AS, 26, 275
  • Scholz & Takeda (1987) Scholz, M. & Takeda, Y. 1987, A&A, 186, 200
  • Stefan (1879) Stefan, J. 1879, Sitzungsberichte der mathematisch-naturwissenschaftlichen Classe der kaiserlichen Akademie der Wissenschaften, 79, 391
  • Takeda et al. (2007) Takeda, G., Ford, E. B., Sills, A., Rasio, F. A., Fischer, D. A., & Valenti, J. A. 2007, ApJS, 168, 297
  • ten Brummelaar et al. (2005) ten Brummelaar, T. A., McAlister, H. A., Ridgway, S. T., Bagnuolo, Jr., W. G., Turner, N. H., Sturmann, L., Sturmann, J., Berger, D. H., Ogden, C. E., Cadman, R., Hartkopf, W. I., Hopper, C. H., & Shure, M. A. 2005, ApJ, 628, 453
  • Tuthill (1994) Tuthill, P. G. 1994, PhD thesis, , Univ. of Cambridge, (1994)
  • Valenti & Fischer (2005) Valenti, J. A. & Fischer, D. A. 2005, ApJS, 159, 141
  • van Belle et al. (1999) van Belle, G. T., Lane, B. F., Thompson, R. R., Boden, A. F., Colavita, M. M., Dumont, P. J., Mobley, D. W., Palmer, D., Shao, M., Vasisht, G. X., Wallace, J. K., Creech-Eakman, M. J., Koresko, C. D., Kulkarni, S. R., Pan, X. P., & Gubler, J. 1999, AJ, 117, 521
  • van Belle & van Belle (2005) van Belle, G. T. & van Belle, G. 2005, PASP, 117, 1263
  • van Belle et al. (2008) van Belle, G. T., van Belle, G., Creech-Eakman, M. J., Coyne, J., Boden, A. F., Akeson, R. L., Ciardi, D. R., Rykoski, K. M., Thompson, R. R., Lane, B. F., & Collaboration, T. P. 2008, ApJS, 176, 276
  • van Leeuwen (2007) van Leeuwen, F. 2007, A&A, 474, 653
  • Zdanavicius et al. (1972) Zdanavicius, K., Nikonov, V. B., Sudzius, J., Straizys, V., Sviderskiene, Z., Kalytis, R., Jodinskiene, E., Meistas, E., Kavaliauskaite, G., Jasevicius, V., Kakaras, G., Bartkevicius, A., Gurklyte, A., Bartkus, R., Azusienis, A., Sperauskas, J., Kazlauskas, A., & Zitkevicius, V. 1972, Vilnius Astronomijos Observatorijos Biuletenis, 34, 3
Star ID Spectral
Points (mas) (mas) (mag) 10 erg cm s Type (mag) (mag)
Control Sample: Stars not known to host planets:
HD1326 216 -
HD4628 98 K1V
HD16160 42 K3V
HD16895 118 F8
HD19373 14 F9.5V
HD20630 2 G5V
HD22484 8 F8V
HD30652 38 F6V
HD39587 84 G0IV-V
HD87901 262 B8IVn
HD88230 64 K8V
HD95735 80 Mb
HD97603 126 A5 IV(n)
HD102647 66 A3Va
HD109358 166 G0V
HD114710 28 G0
HD119850 142 K2
HD126660 134 F8
HD141004 6 G0IV-V
HD142860 58 F5
HD149661 18 K1V
HD157881 26 M1V
HD185144 6 G9V
HD201091 50 K5V
HD201092 16 K7V
HD210027 172 F5
HD215648 248 F5
HD222368 128 F8
EHSA Sample: Known planet hosting stars (PTI):
HD3651 222 K0V
HD9826 540 G0
HD28305 32 K0III
HD75732 16 K0IV-V
HD95128 48 G0
HD117176 192 G0
HD120136 264 F5
HD143761 354 G0V
HD217014 454 G3V
EHSA Sample: Known planet hosting stars (CHARA):
HD3651 K0V
HD11964 G5
HD19994 F8V
HD75732 G8V
HD143761 G0V
HD189733 K1V
HD217014 G3V
Table 6Observed and derived supporting parameters for luminosity class V stars.
HD Template Template
Number (K) ( cm s) (mag) (mas) ()
142 F6V 0.38 78
1237 G8V 3.24 23
2039 G0IV 1.08 23
2638 G8V 1.91 8
3651 K0V 2.52 135
4113 G5V 0.55 22
4208 G5V 1.44 46
4308 G5V 1.12 66
5319 G5IV 5.23 64
6434 G2V 1.65 41
8574 F8V 0.06 23
9826 F8IV 0.64 134
10647 F8V 0.57 61
10697 G5IV 0.28 47
11506 G0V 0.66 30
11964 G5IV 0.47 52
11977 G8III 0.80 36
12661 K0V 1.11 26
13189 K3III 0.98 7
13445 K0V 0.52 82
16141 G5IV 0.53 59
16175 F8IV 0.94 15
17051 F8V 0.24 73
17092 K0III 8.12 5
17156 F8IV 3.13 5