Research Proposal Characterization of Kepler's Transiting

Document Sample
Research Proposal Characterization of Kepler's Transiting Powered By Docstoc
					                      Research Proposal:
 Characterization of Kepler ’s Transiting Extrasolar Giant Planets
                                        Jason W. Barnes

1     Statement of Problem

I propose to help characterize the transiting extrasolar planets detected by NASA’s Kepler
mission by refining and applying transit lightcurve fitting algorithms.
    Unlike HD209458, the parent star of the best-known transiting planet, the stars that
Kepler will monitor do not have accurately measured parallaxes. Therefore the radii of
these stars is unknown. This uncertainty in stellar radius directly translates to uncertainty
in the radius of any transiting planets, which is the most important measurement that Kepler
will make. Hence it is critical that Kepler employ a reliable method for estimating stellar
    Previous investigations into transit lightcurves (Seager and Mall´n-Ornelas, 2003i.e.)
have noted the advantages of minimizing stellar limb darkening when analyzing low-precision
photometric data. I suggest that in circumstances where the photometric precision is high,
stellar limb darkening can be used to an transit observer’s advantage to potentially constrain
the stellar radius without any previous knowledge of the star’s mass or radius. Since the
technique only requires the transit lightcurve, no additional observations are needed. Essen-
tially, if the nature of the stellar limb darkening itself is regular, or predictable, between stars
in Kepler ’s spectral filter, then it may be used as a constraint that will set the transit impact
parameter, constraining the relationship between M∗ , the stellar mass, and R∗ , the stellar
radius. By then applying stellar models similar to those used by Cody and Sasselov (2002),
the M∗ - R∗ degeneracy can be broken, establishing both the star’s radius and its mass. This
procedure will be particularly useful for planets in the middle of Kepler ’s size detectability
range: those whose ingress and egress are not temporally resolved but for which the limb
darkening is photometrically resolved. I would like to investigate this idea quantitatively
during my NRC Postdoctoral Research Associateship.
   Transits can also reveal moons (Sartoretti and Schneider, 1999). Brown et al. (2001)
detected no moons orbiting HD209458b, consistent with theoretical expectation (Barnes
and O’Brien, 2002). Kepler is capable of detecting large moons of extrasolar planets. I will
extend the work of Sartoretti and Schneider (1999), incorporating modern lightcurve-fitting
techniques to the direct transit of extrasolar moons. The resulting algorithm will be used to
search for moons around planets detected by Kepler.
   Subtle (1 part in 105 for zero obliquity, up to 1 in 104 for nonzero obliquity) deviations
during a planet’s transit ingress and egress can reveal the planet’s oblateness (Barnes and
Fortney, 2003), and from that the planet’s rotation rate. For Kepler this precision level

will require the coaddition of many separate transit events in order to achieve the necessary
signal-to-noise ratio, and will still be difficult. However, for low-density, rapidly rotating
planets at intermediate semimajor axis (0.2 < a < 1.6) it may be possible to measure
rotation from the Kepler transit lightcurve.
    In Barnes and Fortney (2004) I showed that large, Saturn-like ring systems, in the right
viewing geometry, can show whopping transit lightcurve detectabilities on the order of 1 ×
10−3 . The potentially intricate complexity of such ring systems, like that of Saturn, make
fitting a lightcurve for ring sizes and optical depths computationally intensive. In order to
fit potentially anomalous Kepler transit lightcurves, I propose to develop and refine ringed
planet transit fitting algorithms until they are capable of fitting Kepler data on a reasonable
timescale. I will then apply the fitting scheme to the Kepler transiting planets.

2    Background and Relevance

Although over 150 extrasolar planets have been discovered (see,
characterization has been possible only for those 9 that transit. Radial velocity monitoring
alone measures planets’ minimum masses and orbital characteristics (Marcy et al., 2000,
e.g.). Precision transit photometry reveals a planet’s orbital inclination and radius. To-
gether, the two techniques are able to nail down the planet’s actual mass and therefore its
    This analysis was first done for transiting planet HD209458b (Charbonneau et al., 2000;
Henry et al., 2000). The edge-on inclination of HD209458b’s orbit removed the sin i degen-
eracy from the radial-velocity results, pinning the planet’s mass at 0.69 Jupiter masses (M J ).
Subsequent Hubble Space Telescope photometry of the transit showed no evidence for moons
or rings (Brown et al., 2001), consistent with theoretical constraints placed on extrasolar
moons (Barnes and O’Brien, 2002). Though its parallax was measured by Hipparcos, the
leading source of error in the determination of HD209458b’s radius remained the uncertainty
in the radius of its parent star, HD209458. The current best estimate for HD209458b’s ra-
dius, 1.42+0.10 RJ , was calculated by Cody and Sasselov (2002) by coupling the HST transit
lightcurve fit to a stellar structure model. However, this sophisticated treatment requires
knowledge of the star’s luminosity, which cannot be calculated without knowledge of the
star’s distance.
   The transits of HD209458b allow additional follow-up observations that have made HD209458b
the best-characterized extrasolar planet to date. Charbonneau et al. (2002) detected sodium
in HD209458b’s outer atmosphere from HST spectrophotometry. Hydrogen, carbon, and
oxygen were all detected in the process of hydrodynamic escape from HD209458b’s atmo-
sphere (Vidal-Madjar et al., 2003; ?). The Spitzer Space Telescope detected infrared emission
from HD209458b’s dayside during its secondary eclipse, directly measuring its day temper-
ature to be 1100K (?). Very recently, ?) used transit observations to measure the relative

angle between the planet’s orbital plane and the stellar equatorial plane to be 4.4◦ .
   Since 1999, 8 additional transiting planets have been discovered. Like HD209458b, two
others were first found using radial velocity techniques but were later found to transit. Both
were discovered in the last 6 months. HD149026b (Mp = 0.36MJ , Rp = 0.73RJ ) (?) and
HD189733b (Mp = 1.15MJ , Rp = 1.26RJ ) (?) are both close-in hot Jupiters. HD149026b in
particular is important because although its mass is similar to that of Saturn, its radius is
lower despite high insolation, implying that the planet has a large core (?). Both of these
planets are around stars bright enough for significant follow-up observations.
    Six extrasolar planets have been discovered using the transit technique so far. TrES-1 (?),
an 0.75 MJ , 1.08 RJ planet in a 3 day orbit, was found from its transits of a relatively bright
(V=11.79) star based on a wide-field transit search. As such, TrES-1 is bright enough for
follow-up observations which have included the detection of dayside infrared emissions using
the Spitzer Space Telescope (Charbonneau et al., 2005). The Optical Gravitational Lensing
Experiment (OGLE) team has spent one month per year on a narrow-field, deep search
for planetary transits (Udalski et al., 2002a,b). Weeding out false positives from possible
transits is the hardest part of the process – less than 10% of prospective transits have yielded
planets. There are currently 5 generally-accepted OGLE planets: OGLE-TR-56b (Konacki
et al., 2003), OGLE-TR-113b, OGLE-TR-132b (Bouchy et al., 2004), OGLE-TR-111b (?),
and OGLE-TR-10b (?). The OGLE planets are all far away and their parent stars are faint.
Ground-based follow-up for them is challenging.
    The prospects for characterizing extrasolar giant planets and their environs are very
favorable with the 2008 launch of NASA’s Kepler mission(?). Kepler is a 0.95-m aperture
photometer designed to continuously monitor 100,000 stars for 4 years searching for transiting
planets. Although designed to find Earth-sized worlds, Kepler is also expected to discover
30 transiting giant planets in large (¿ 1.6 AU) orbits, 35 in moderate (0.2 ¡ a ¡ 1.6 AU)
orbits, and 100 in close-in orbits. Any Neptune-mass objects that Kepler finds will also be
amenable to characterization. Furthermore, Kepler ’s high photometric precision (∼ 1 × 10 −4
or better for 15 minute data points, depending on the target) will result in lightcurves of
comparable or superior quality to Hubble Space Telescope observations of HD209458 (Brown
et al., 2001) without the need for targeted follow-up.

3     Methodology, Procedures, and Techniques

Transit lightcurve fitting algorithms measure just 4 parameters: b, the transit impact pa-
rameter, related to the planet’s orbital inclination; R∗ , the stellar radius; Rp , the planetary
radius; and c1, a composite limb-darkening coefficient (after Brown et al. (2001)). Roughly
speaking, the depth (fraction of the star’s light blocked in mid-transit) of the transit sets the
ratio R∗ . The curvature of the transit bottom constrains c1. The duration of ingress and
egress establishes b, and the total transit duration sets the length of the transit chord which

then yields R∗ . These constraints rely on previous knowledge of the star’s mass, M∗ .
    Seager and Mall´n-Ornelas (2003) showed that, in cases where limb darkening is negligi-
ble, the length of ingress and egress can be directly measured, and an analytical solution to
the mappings discussed in the previous paragraph exists. However, Kepler operates in visi-
ble wavelengths where limb darkening is not negligible. For planets with Rp ∼ RJ Kepler ’s
excellent photometric precision will still resolve the end of ingress and beginning of egress
sufficiently for a model-based fitting routine to measure the 4 basic transit parameters.
    However, as the duration of ingress and egress depends on the radius of the transiting
planet, for planets with Rp < 0.3RJ or so the ingress and egress will not be temporally
resolved. In this case, the stellar limb darkening may be capable of constraining the transit
impact parameter, thus taking the place of the ingress and egress duration in model fits.
To test this hypothesis, I will alter my transit-fitting algorithm developed for Barnes and
Fortney (2003) and Barnes and Fortney (2004) to accomodate the new procedure. I will then
generate synthetic lightcurves simulating Kepler transit photometry and fit those data using
both fixed and variable limb darkening parameters. The process may require prior knowledge
of the precise value of c1, or it may be that the mere knowledge that the star’s limb darkening
structure can be wholly determined by the single parameter c1 may be sufficient.
    Either way, this technique will only establish a relationship between M∗ and R∗ , and will
not by itself unambiguously determine the stellar radius. The application of stellar modelling,
of a type similar to that pioneered for use with planetary transits by Cody and Sasselov
(2002), may resolve the degeneracy. Along with ground-based astrometric measurements of
stellar parallax, application of the star’s spectral type should be able to limit the allowable
range of M∗ and R∗ as set by the lightcurve fit. I do not yet know what level of uncertainty
in R∗ will remain after this type of analysis, but determination of the expected residual
uncertainty will be one of the research goals.
    A possible back-door method for determining M∗ involves moons. If a moon were detected
in transit, its period, along with multiple determinations of its distance from the parent
planet, would independently and unambiguously determine Mp . The amplitude and period
of later-measured planet-induced radial velocity variations of the star could then set the
stellar mass, and thereby allow easy determination of the stellar radius from the planet’s
transit lightcurve. In Barnes and O’Brien (2002) I showed that large, detectable moons
should not exist around close-in giant planets. However, moons larger than 1R⊕ can exist
around giant planets with semimajor axes larger than ∼ 0.2 AU, and can serve to measure
the planet’s mass.
    Direct measurement of the semimajor axis of an extrasolar moon may be challenging. I
will first develop the capability to fit the transit lightcurve of an extrasolar moon incorpo-
rating orbital motion on a variably inclined orbital plane, and then apply the technique to
synthetic Kepler lightcurves to evaluate its feasibility and limitations. The resulting program
will be used to fit anomalous lightcurves from Kepler after its launch in 2008.

   The algorithm for fitting planetary oblateness has already been written (Barnes and
Fortney, 2003), and can be used on Kepler lightcurves with little modification.
    However, no algorithm exists for fitting the transits of ringed planets. The reduced
symmetry of the problem makes a brute-force fitting approach glacially slow – fits would
take many months. I will make optimizing improvements in the fitting algorithm in order to
fit transiting ringed planets on a reasonable timescale. The technique will assume a single
ring of uniform optical depth (to be fit) and an interior and exterior radius to be fit, along
with the planetary obliquity.

4    Expected Results and Significance

The Kepler mission will establish the types and distribution of terrestrial planets and plan-
etary systems in the Sun’s vicinity. It will also discover and characterize hundreds of giant
    My proposed contribution to the project will be to develop a sophisticated transit lightcurve
fitting algorithm for eventual application to Kepler data. In the meantime, I will publish
my studies of the usefulness and limitations of the new techniques.
    Planetary radius is the fundamental measurable quantity that Kepler will provide, and
its measurement of that value will depend on the stellar radius. My algorithm and analysis
will help to constrain stellar radii throughout the Kepler project.
   Moons of extrasolar planets are exciting in their own right and as possible habitable
worlds. Moons also might be used to independently measure the mass of an extrasolar
planet. The types and distribution of extrasolar moons will help to constrain the formation
and evolution of planets and planetary systems.
   I also expect to be able to measure the rotation rate of perhaps a handful of extrasolar
giant planets based on a determination of oblateness from a transit lightcurve (Barnes and
Fortney, 2003). This analysis can also constrain the planet’s obliquity (Seager and Hui,
2002). These measurements also constrain formation and evolution histories, and may also
help to reveal the yet-unknown tidal dissipation mechanism within giant planets (Barnes
and Fortney, 2003).
    In our own solar system, the origin, age, and longevity of Saturn’s spectacular ring
system remain unresolved. If they exist, Kepler will find large, Saturn-like ring systems
around extrasolar giant planets. The distribution of these extrasolar ring systems will shed
light on the prevalence of rings in general, as well as provide fuel for a new understanding
of the workings of rings.

Barnes, J. W. and J. J. Fortney, 2003: Measuring the Oblateness and Rotation of Transiting
  Extrasolar Giant Planets. ApJ , 588, 545–556.

—, 2004: Transit Detectability of Ring Systems around Extrasolar Giant Planets. ApJ , 616,

Barnes, J. W. and D. P. O’Brien, 2002: Stability of Satellites around Close-in Extrasolar
  Giant Planets. ApJ , 575, 1087–1093.

Bouchy, F., F. Pont, N. C. Santos, C. Melo, M. Mayor, D. Queloz, and S. Udry, 2004: Two
  new “very hot Jupiters” among the OGLE transiting candidates. A&A, 421, L13–L16.

Brown, T. M., D. Charbonneau, R. L. Gilliland, R. W. Noyes, and A. Burrows, 2001:
  Hubble Space Telescope Time-Series Photometry of the Transiting Planet of HD 209458.
  ApJ , 552, 699–709.

Charbonneau, D., L. E. Allen, S. T. Megeath, F. Torres, T. M. Brown, R. L. Gilliland, D. W.
  Latham, G. Mandushev, F. T. O’Donovan, R. Alonso, and A. Sozzetti, 2005: Detectyion
  of Thermal Emission from an Extrasolar Planet. ApJ , in press.

Charbonneau, D., T. M. Brown, D. W. Latham, and M. Mayor, 2000: Detection of Planetary
  Transits Across a Sun-like Star. ApJ , 529, L45–L48.

Charbonneau, D., T. M. Brown, R. W. Noyes, and R. L. Gilliland, 2002: Detection of an
  Extrasolar Planet Atmosphere. ApJ , 568, 377–384.

Cody, A. M. and D. D. Sasselov, 2002: HD 209458: Physical Parameters of the Parent Star
  and the Transiting Planet. ApJ , 569, 451–458.

Henry, G. W., G. W. Marcy, R. P. Butler, and S. S. Vogt, 2000: A Transiting “51 Peg-like”
  Planet. ApJ , 529, L41–L44.

Konacki, M., G. Torres, S. Jha, and D. Sasselov, 2003: A New Transiting Extrasolar Giant
 Planet. Nature.

Marcy, G. W., W. D. Cochran, and M. Mayor, 2000: Extrasolar Planets around Main-
 Sequence Stars. In Protostars and Planets IV. (Tucson: University of Arizona Press), p.

Sartoretti, P. and J. Schneider, 1999: On the detection of satellites of extrasolar planets with
  the method of transits. A&AS , 134, 553–560.

Seager, S. and L. Hui, 2002: Constraining the Rotation Rate of Transiting Extrasolar Planets
  by Oblateness Measurements. ApJ , 574, 1004–1010.

Seager, S. and G. Mall´n-Ornelas, 2003: A Unique Solution of Planet and Star Parameters
  from an Extrasolar Planet Transit Light Curve. ApJ , 585, 1038–1055.

Udalski, A., B. Paczynski, K. Zebrun, M. Szymaski, M. Kubiak, I. Soszynski, O. Szewczyk,
 L. Wyrzykowski, and G. Pietrzynski, 2002a: The Optical Gravitational Lensing Exper-
 iment. Search for Planetary and Low-Luminosity Object Transits in the Galactic Disk.
 Results of 2001 Campaign. Acta Astronomica, 52, 1–37.

Udalski, A., K. Zebrun, M. Szymanski, Kubiak, I. Soszynski, O. Szewczyk, L. Wyrzykowski,
 and G. Pietrzynski, 2002b: The Optical Gravitational Lensing Experiment. Search for
 Planetary and Low- Luminosity Object Transits in the Galactic Disk. Results of 2001
 Campaign – Supplement. Acta Astronomica, 52, 115–128.

Vidal-Madjar, A., A. Lecavelier des Etangs, J.-M. D´sert, G. E. Ballester, R. Ferlet,
  G. H´brard, and M. Mayor, 2003: An extended upper atmosphere around the extraso-
  lar planet HD209458b. Nature, 422, 143–146.