ABSTRACT
We measured the angular diameter of the exoplanet host star ι Dra with Georgia State University's Center for High Angular Resolution Astronomy Array interferometer and, using the star's parallax and photometry from the literature, calculated its physical radius and effective temperature. We then combined our results with stellar oscillation frequencies from Zechmeister et al. and orbital elements from Kane et al. to determine the masses for the star and exoplanet. Our value for the central star's mass is 1.82 ± 0.23 M☉, which means the exoplanet's minimum mass is 12.6 ± 1.1 MJupiter. Using our new effective temperature, we recalculated the habitable zone for the system, though it is well outside the star–planet separation.
Export citation and abstract BibTeX RIS
1. INTRODUCTION
Frink et al. (2002) announced the discovery of a substellar companion to ι Draconis (K2 III, HD 137759) with a period of 536 days and a minimum mass for the companion of 8.9 MJupiter. They used a stellar mass of 1.05 M☉ from Allende Prieto & Lambert (1999), who compared the absolute visual magnitude and (B − V) color from Hipparcos data with theoretical isochrones from Bertelli et al. (1994). However, Frink et al. acknowledge that evolutionary tracks for a range of masses are close together on the H-R diagram, so any slight change in the evolutionary model can have a large impact on the derived mass.
Zechmeister et al. (2008, hereafter Z08) observed ι Dra in search of stellar oscillations using three separate instruments over almost eight years in order to refine the orbital parameters of the planet and determine the mass of the central star. They found low-amplitude, solar-like oscillations with a frequency of 3.8 day−1 in two of the data sets and derived a stellar mass of 2.2 M☉ using the equations
and
where νosc is the oscillation velocity amplitude, fmax is the frequency of the strongest mode, and Teff is the effective temperature. They used a luminosity of 64.2 ± 2.1 L☉ from the Hipparcos catalog and Teff = 4490 K from McWilliam (1990). Z08 then compared their 2.2 M☉ value to those derived using Teff, surface gravities (log g), and metallicities ([Fe/H]) from the literature and the PARAM stellar model by Girardi et al. (2000) and da Silva et al. (2006).4 The masses ranged from 1.05 ± 0.36 M☉ based on values from Allende Prieto & Lambert (1999) to 1.71 ± 0.38 M☉ based on values from Santos et al. (2004). Because Z08's mass was significantly higher than those derived using the model, they chose a mass of 1.4 M☉ to calculate the minimum mass of the companion, which they list as 10.3 MJupiter.
A more accurate way to estimate the star's mass would be to investigate the frequency splitting (Δf0) using the equation
combined with an interferometrically measured radius, but unfortunately Z08's data set was not suitable for measuring Δf0.
The advantage interferometry brings is the ability to directly measure the angular diameter of the star. Then, the physical radius can be determined using the distance from the parallax, and Teff can be calculated. We combine our results with those from stellar oscillation frequencies to more completely understand the system through a description of the central star's and exoplanet's masses and the extent of the habitable zone. Section 2 details our observing procedure, Section 3 discusses how ι Dra's angular diameter and Teff were determined, and Section 4 explores the physical implications of the new measurements.
2. INTERFEROMETRIC OBSERVATIONS
Observations were obtained using the Center for High Angular Resolution Astronomy (CHARA) Array, a six-element Y-shaped optical–infrared interferometer located on Mount Wilson, California (ten Brummelaar et al. 2005). We used the CHARA Classic and CLIMB beam combiners in the K' band (2.13 μm) while visible wavelengths (470–800 nm) were used for tracking and tip/tilt corrections. The observing procedure and data reduction process employed here are described in McAlister et al. (2005). We observed ι Dra over four nights spanning four years with two baselines: in 2007 and 2008, we used the longest telescope pair S1–E1 with a maximum baseline of 331 m, and in 2011, we used a shorter telescope pair W1–W2 with a maximum baseline of 108 m.5
Choosing proper calibrator stars is vital because they act as the standard against which the scientific target is measured. We selected four calibrators (HD 128998, HD 139778, HD 141472, and HD 145454) because they are single stars with expected visibility amplitudes >80% so they were very nearly unresolved on the baseline used. This meant uncertainties in the calibrators' diameters did not affect the target's diameter calculation as much as if the calibrator stars had a significant angular size on the sky. We interleaved calibrator and target star observations so that every target was flanked by calibrator observations made as close in time as possible, which allowed us to convert instrumental target and calibrator visibilities to calibrated visibilities for the target.
To check for possible unseen close companions that would contaminate our observations, we created spectral energy distribution (SED) fits based on published U BV RI JHK photometric values obtained from the literature for each calibrator to establish diameter estimates. We combined the photometry with Kurucz model atmospheres6 based on Teff and log g values to calculate angular diameters for the calibrators. The stellar models were fit to observed photometry after converting magnitudes to fluxes using Colina et al. (1996, UBVRI) and Cohen et al. (2003, JHK). The photometry, Teff and log g values, and resulting angular diameters for the calibrators are listed in Table 1. There were no hints of excess emission associated with a low-mass stellar companion or circumstellar disk in the calibrators' SED fits (see Figure 1).
Table 1. Calibrator Information
Parameter | HD 128998 | HD 139778 | HD 141472 | HD 145454 |
---|---|---|---|---|
U magnitude | 5.84 | N/A | N/A | 5.35 |
B magnitude | 5.84 | 6.94 | 7.31 | 5.42 |
V magnitude | 5.85 | 5.87 | 5.92 | 5.44 |
R magnitude | 5.82 | 5.21 | 5.09 | 5.46 |
I magnitude | 5.84 | 4.65 | 4.23 | 5.50 |
J magnitude | 5.76 | 4.15 | 3.55 | 5.37 |
H magnitude | 5.80 | 3.56 | 2.85 | 5.43 |
K magnitude | 5.76 | 3.30 | 2.49 | 5.43 |
Extinction AV | N/A | 0.13 | 0.13 | N/A |
Teff (K) | 9395 | 4525 | 4200 | 9772 |
log g (cm s−2) | 4.14 | 1.94 | 1.94 | 4.13 |
θUD (mas) | 0.23 ± 0.01 | 1.08 ± 0.08 | 1.48 ± 0.08 | 0.27 ± 0.02 |
Notes. The photometric values are from the following sources: (U BV) Mermilliod 1991; (RI) Monet et al. 2003; (JHK) Cutri et al. 2003. The AV values are from Famaey et al. (2005). All Teff and log g values are from Cox (2000) and are based on the star's spectral type, except for HD 145454, which is from Allende Prieto & Lambert (1999). The uniform-disk angular diameters (θUD) are the result of the SED-fitting procedure described in Section 2.
Download table as: ASCIITypeset image
3. DETERMINATION OF ANGULAR DIAMETER AND Teff
The observed quantity of an interferometer is defined as the visibility (V), which is fit with a model of a uniformly illuminated disk (UD) that represents the observed face of the star. Diameter fits to V were based upon the UD approximation given by V = 2J1(x)/x, where J1 is the first-order Bessel function and x = πBθUDλ−1, where B is the projected baseline at the star's position, θUD is the apparent UD angular diameter of the star, and λ is the effective wavelength of the observation (Shao & Colavita 1992). Table 2 lists the Modified Julian Date (MJD), projected baseline at the time of observation (B), projected baseline position angle (Θ), calibrated visibility (V), and error in V (σV) for ι Dra.
Table 2. ι Dra's Calibrated Visibilities
Calib | MJD | B | Θ | V | σV |
---|---|---|---|---|---|
HD | (m) | (deg) | |||
145454 | 54358.188 | 306.28 | 115.2 | 0.041 | 0.003 |
54358.198 | 304.57 | 117.5 | 0.039 | 0.003 | |
54358.211 | 301.74 | 120.7 | 0.033 | 0.002 | |
145454 | 54643.239 | 311.24 | 253.6 | 0.038 | 0.005 |
54643.252 | 312.42 | 256.8 | 0.043 | 0.004 | |
54643.265 | 313.31 | 256.0 | 0.040 | 0.004 | |
54643.276 | 313.91 | 262.9 | 0.047 | 0.008 | |
54643.287 | 314.27 | 265.5 | 0.053 | 0.010 | |
54643.298 | 314.47 | 268.2 | 0.041 | 0.007 | |
128998 | 54648.261 | 313.81 | 262.3 | 0.044 | 0.007 |
54648.273 | 314.23 | 265.1 | 0.046 | 0.008 | |
54648.283 | 314.44 | 267.7 | 0.044 | 0.008 | |
145454 | 54648.354 | 311.66 | 105.3 | 0.051 | 0.009 |
54648.366 | 310.44 | 108.2 | 0.042 | 0.005 | |
54648.378 | 308.91 | 111.1 | 0.038 | 0.006 | |
54648.389 | 307.30 | 113.7 | 0.042 | 0.005 | |
54648.400 | 305.43 | 116.4 | 0.038 | 0.005 | |
54648.410 | 303.34 | 118.9 | 0.035 | 0.005 | |
54648.421 | 300.95 | 121.5 | 0.030 | 0.004 | |
139778 | 55775.159 | 106.94 | 184.9 | 0.372 | 0.048 |
55775.169 | 106.39 | 188.2 | 0.377 | 0.045 | |
55775.180 | 105.70 | 191.6 | 0.384 | 0.032 | |
55775.192 | 104.84 | 195.3 | 0.356 | 0.051 | |
141472 | 55775.169 | 106.39 | 188.2 | 0.381 | 0.058 |
55775.180 | 105.70 | 191.6 | 0.360 | 0.055 | |
55775.192 | 104.84 | 195.3 | 0.343 | 0.056 |
Note. The projected baseline position angle (Θ) is measured eastward from north.
Download table as: ASCIITypeset image
Figure 2 shows the UD fit to the observed visibilities and it is clear that this is not a sufficient model. A more realistic model of a star's disk involves limb darkening (LD), particularly in this case because most of the observations are on the third lobe of the visibility curve where secondary effects such as LD play a more important role than in the curve between B = 0 and the first null. Lacour et al. (2008) analyzed multiple LD prescriptions for the K1.5 III star Arcturus and found the power-law model to be a sufficient approximation. We chose to use this model because Arcturus' spectral type is very close to that of ι Dra. The model was based on Hestroffer (1997):
where I(μ) is the brightness of a point source at wavelength λ, I(1) is the brightness at the center, , r is the angular distance from the star's center, and θLD is the limb-darkened angular diameter. In terms of the star's visibilities, the power-law prescription becomes
where vr is the radial spatial frequency, Γ is the Euler function (Γ(k+1) = k!), α = 0.258 ± 0.003 (Lacour et al. 2008), and k is the number of terms used. We tried a range of k values to check its effect on θLD and found that after 10 terms, θLD remained steady. Therefore we used k = 11.
Download figure:
Standard image High-resolution imageThe resulting UD and LD angular diameters and their errors (<1%) are listed in Table 3. Additionally, the combination of the interferometric measurement of the star's angular diameter plus the Hipparcos parallax (van Leeuwen 2007) allowed us to determine the star's physical radius, which is also listed in Table 3.
Table 3. ι Dra Stellar Parameters
Parameter | Value | Reference |
---|---|---|
[Fe/H] | +0.14 | Averaged from Soubiran et al. (2010) |
V magnitude | 3.29 ± 0.02 | Mermilliod (1991) |
K magnitude | 0.72 ± 0.04 | Neugebauer & Leighton (1969) |
AV | 0.03 | Famaey et al. (2005) |
BC | 0.42 ± 0.10 | Alonso et al. (1999) |
Luminosity (L☉) | 55.3 ± 5.3 | Calculated here |
FBOL (10−8 erg s−1 cm−2) | 183.8 ± 17.7 | Calculated here |
θUD (mas) | 3.434 ± 0.012 (0.3%) | Measured here |
θLD (mas) | 3.596 ± 0.015 (0.4%) | Measured here |
Rlinear (R☉) | 11.99 ± 0.06 (0.5%) | Calculated here |
Teff (K) | 4545 ± 110 (2%) | Calculated here |
Download table as: ASCIITypeset image
For the θLD fit, the errors were derived via the reduced χ2 minimization method (Wall & Jenkins 2003; Press et al. 1992): the diameter fit with the lowest χ2 was found and the corresponding diameter was the final θLD for the star. The errors were calculated by finding the diameter at χ2 + 1 on either side of the minimum χ2 and determining the difference between the χ2 diameter and χ2 + 1 diameter. Figure 3 shows the LD diameter fit for ι Dra.
Download figure:
Standard image High-resolution imageOnce θLD was determined interferometrically, the Teff was calculated using the relation
where FBOL is the bolometric flux and σ is the Stefan–Boltzmann constant. FBOL was determined in the following way: ι Dra's V and K magnitudes were dereddened using the extinction curve described in Cardelli et al. (1989) and its interstellar absorption (AV) value was from Famaey et al. (2005). The intrinsic broadband color (V − K) was calculated and the bolometric correction (BC) was determined by interpolating between the [Fe/H] = 0.0 and +0.2 tables from Alonso et al. (1999). They point out that in the range of 6000 K ⩾Teff ⩾ 4000 K, their BC calibration is symmetrically distributed around a ±0.10 mag band when compared to other calibrations, so we assigned the BC an error of 0.10. The FBOL was then determined by applying the BC and the Teff was calculated. The star's luminosity (L) was also calculated using the absolute V magnitude and the BC. See Table 3 for a summary of these parameters.
4. RESULTS AND DISCUSSION
We estimated the limb-darkened angular diameter for ι Dra using two additional methods as a check for our measurement. First, we created an SED fit for the star as described in Section 2, where UBV photometry was from Mermilliod (1991), RI photometry was from Monet et al. (2003), and JHK photometry was from Cutri et al. (2003). Figure 4 shows the resulting fit. Second, we used the relationship described in Blackwell & Lynas-Gray (1994) between the (V − K) color, Teff, and θ. Our measured θLD is 3.596 ± 0.015 mas, the SED fit estimates 3.81 ± 0.23 mas, and the color–temperature–diameter relationship produces 3.63 ± 0.53 mas. Because ι Dra is so bright in the K band (0.7 mag) and because Two Micron All Sky Survey (2MASS) measurements saturate at magnitudes brighter than ∼3.5 even when using the shortest exposure time,7 we used the K magnitude from the Two-Micron Sky Survey (Neugebauer & Leighton 1969) for the color–diameter determination.
Download figure:
Standard image High-resolution imageThe main sources of errors for the three methods are uncertainties in visibilities for the interferometric measurement, uncertainties in the comparison between observed fluxes and the model fluxes for a given Teff and log g for the SED estimate, and uncertainties in the parameters of the relation and the spread of stars around that relation for the color–temperature–diameter determination. All three diameters agree within their errors but our interferometric measurements provide an error approximately 15 and 35 times smaller than the other methods, respectively.
With our newly calculated Teff, we were able to estimate the mass of the central star using Equations (1) and (2) and obtained a mass of 1.82 ± 0.23 M☉. We then calculated the exoplanet's minimum mass using the orbital parameters presented in Kane et al. (2010) and the equation
where the period P was 510.72 ± 0.07 days, the amplitude K was 306.0 ± 3.8 m s−1, and the eccentricity e was 0.713 ± 0.008. Our calculation produced a minimum mass of 12.6 ± 1.1 MJupiter, which converged in two iterations.
Our Teff of 4545 ± 110 K is within the range listed in Z08, which spans 4466 ± 100 K (Allende Prieto & Lambert 1999) to 4775 ± 113 K (Santos et al. 2004). On the other hand, the stellar mass derived here is lower than that calculated in Z08 (2.2 M☉) and slightly higher than the range presented in their paper (1.2–1.7 M☉), though it overlaps within errors with the mass derived from the Teff and log g in Santos et al. (2004). Hopefully, future observations of this star will determine the frequency splitting (see Equation (3)), which will allow for the direct measurement of the star's mass when combined with the interferometrically measured radius.
Using the following equations from Jones et al. (2006), we were also able to calculate the size of the system's habitable zone:
and
where Sb, i(Teff) and Sb, o(Teff) are the critical fluxes at the inner and outer boundaries in units of the solar constant. The inner and outer physical boundaries ri, o in AU were then calculated using
We obtained habitable zone boundaries of 6.8 AU and 13.5 AU. ι Dra's planet has a semimajor axis of 1.34 AU (Z08), so there is no chance that the planet orbits anywhere near the habitable zone.
The CHARA Array is funded by the National Science Foundation through NSF grant AST-0606958 and by Georgia State University through the College of Arts and Sciences, and by the W.M. Keck Foundation. S.T.R. acknowledges partial support by NASA grant NNH09AK731. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. 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.
Footnotes
- 4
Available online at http://stev.oapd.inaf.it/cgi-bin/param_1.0.
- 5
The three arms of the CHARA Array are denoted by their cardinal directions: "S," "E," and "W" are south, east, and west, respectively. Each arm bears two telescopes, numbered "1" for the telescope farthest from the beam combining laboratory and "2" for the telescope closer to the lab. The "baseline" is the distance between the telescopes.
- 6
Available to download at http://kurucz.cfa.harvard.edu.
- 7
Explanatory Supplement to the 2MASS All Sky Data Release and Extended Mission Products, http://www.ipac.caltech.edu/2mass/releases/allsky/doc/.