Source author record

David Charbonneau

David Charbonneau appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

122works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

122 published item(s)

preprint2022arXiv

A Possible Alignment Between the Orbits of Planetary Systems and their Visual Binary Companions

Astronomers do not have a complete picture of the effects of wide-binary companions (semimajor axes greater than 100 AU) on the formation and evolution of exoplanets. We investigate these effects using new data from Gaia EDR3 and the TESS mission to characterize wide-binary systems with transiting exoplanets. We identify a sample of 67 systems of transiting exoplanet candidates (with well-determined, edge-on orbital inclinations) that reside in wide visual binary systems. We derive limits on orbital parameters for the wide-binary systems and measure the minimum difference in orbital inclination between the binary and planet orbits. We determine that there is statistically significant difference in the inclination distribution of wide-binary systems with transiting planets compared to a control sample, with the probability that the two distributions are the same being 0.0037. This implies that there is an overabundance of planets in binary systems whose orbits are aligned with those of the binary. The overabundance of aligned systems appears to primarily have semimajor axes less than 700 AU. We investigate some effects that could cause the alignment and conclude that a torque caused by a misaligned binary companion on the protoplanetary disk is the most promising explanation.

preprint2022arXiv

Constraints on the Spindown of Fully-Convective M Dwarfs Using Wide Field Binaries

M dwarfs remain active over longer timescales than their Sunlike counterparts, with potentially devastating implications for the atmospheres of their planets. However, the age at which fully-convective M dwarfs transition from active and rapidly rotating to quiescent and slowly rotating is poorly understood, as these stars remain rapidly rotating in the oldest clusters that are near enough for a large sample of low-mass M dwarfs to be studied. To constrain the spindown of these low-mass stars, we measure photometric rotation periods for field M dwarfs in wide binary systems, primarily using TESS and MEarth. Our analysis includes M-M pairs, which are coeval but of unknown age, as well as M dwarfs with white dwarf or Sunlike primaries, for which we can estimate ages using techniques like white dwarf cooling curves, gyrochronology, and lithium abundance. We find that the epoch of spindown is strongly dependent on mass. Fully-convective M dwarfs initially spin down slowly, with the population of 0.2--0.3M$_\odot$ rapid rotators evolving from $P_{\rm rot} < 2$ days at 600 Myr to $2 < P_{\rm rot} < 10$ days at 1--3 Gyr before rapidly spinning down to long rotation periods at older ages. However, we also identify some variability in the spindown of fully-convective M dwarfs, with a small number of stars having substantially spun down by 600 Myr. These observations are consistent with models of magnetic morphology-driven spindown, where angular momentum loss is initially inefficient until changes in the magnetic field allow spindown to progress rapidly.

preprint2022arXiv

Galactic Kinematics and Observed Flare Rates of a Volume-Complete Sample of Mid-to-Late M-dwarfs: Constraints on the History of the Stellar Radiation Environment of Planets Orbiting Low-mass Stars

We present a study of the relationship between galactic kinematics, flare rates, chromospheric activity, and rotation periods for a volume-complete, nearly all-sky sample of 219 single stars within 15 parsecs and with masses between 0.1$-$0.3 M$_{odot}l$ observed during the primary mission of TESS. We find all stars are consistent with a common value of $α$=1.984 $\pm$ 0.019 for the exponent of the flare frequency distribution. Using our measured stellar radial velocities and Gaia astrometry, we determine galactic UVW space motions. We find 64% of stars are members of the Galactic thin disk, 5% belong to the thick disk, and for the remaining 31%, we cannot confidently assign membership to either component. If we assume star formation has been constant in the thin disk for the past 8 Gyr, then based on the fraction that we observe to be active, we estimate the average age at which these stars transition from the saturated to the unsaturated flaring regime to be 2.4 $\pm$ 0.3 Gyr. This is consistent with the ages that we assign from galactic kinematics: We find that stars with Prot $<$ 10 days have an age of 2.0 $\pm$ 1.2 Gyr, stars with 10 $<$ Prot $\leq$ 90 days have an age of 5.6 $\pm$ 2.7 Gyr, and stars with Prot $>$ 90 days have an age of 12.9 $\pm$ 3.5 Gyr. We find that the average age of stars with Prot $<$ 10 days increases with decreasing stellar mass from 0.6 $\pm$ 0.3 Gyr (0.2 - 0.3 M$_{odot}$) to 2.3 $\pm$ 1.3 Gyr (0.1--0.2 M$_{odot}l$).

preprint2022arXiv

Independent validation of the temperate Super-Earth HD79211 b using HARPS-N

We present high-precision radial velocities (RVs) from the HARPS-N spectrograph for HD79210 and HD79211, two M0V members of a gravitationally-bound binary system. We detect a planet candidate with a period of $24.421^{+0.016}_{-0.017}$ days around HD79211 in these HARPS-N RVs, validating the planet candidate originally identified in CARMENES RV data alone. Using HARPS-N, CARMENES and HIRES RVs spanning a total of 25 years, we further refine the planet candidate parameters to $P=24.422\pm0.014$ days, $K=3.19\pm0.27$ m/s, $M$ sin $i = 10.6 \pm 1.2 M_\oplus$, and $a = 0.142 \pm0.005$ au. We do not find any additional planet candidate signals in the data of HD79211 nor do we find any planet candidate signals in HD79210. This system adds to the number of exoplanets detected in binaries with M dwarf members, and serves as a case study for planet formation in stellar binaries.

preprint2022arXiv

The LHS 1678 System: Two Earth-Sized Transiting Planets and an Astrometric Companion Orbiting an M Dwarf Near the Convective Boundary at 20 pc

We present the TESS discovery of the LHS 1678 (TOI-696) exoplanet system, comprised of two approximately Earth-sized transiting planets and a likely astrometric brown dwarf orbiting a bright ($V_J$=12.5, $K_s$=8.3) M2 dwarf at 19.9 pc. The two TESS-detected planets are of radius 0.70$\pm$0.04 $R_\oplus$ and 0.98$\pm$0.06 $R_\oplus$ in 0.86-day and 3.69-day orbits, respectively. Both planets are validated and characterized via ground-based follow-up observations. HARPS RV monitoring yields 97.7 percentile mass upper limits of 0.35 $M_\oplus$ and 1.4 $M_\oplus$ for planets b and c, respectively. The astrometric companion detected by the CTIO/SMARTS 0.9m has an orbital period on the order of decades and is undetected by other means. Additional ground-based observations constrain the companion to being a high-mass brown dwarf or smaller. Each planet is of unique interest; the inner planet has an ultra-short period, and the outer planet is in the Venus zone. Both are promising targets for atmospheric characterization with the JWST and mass measurements via extreme-precision radial velocity. A third planet candidate of radius 0.9$\pm$0.1 $R_\oplus$ in a 4.97-day orbit is also identified in multi-Cycle TESS data for validation in future work. The host star is associated with an observed gap in the lower main sequence of the Hertzsprung-Russell diagram. This gap is tied to the transition from partially- to fully-convective interiors in M dwarfs, and the effect of the associated stellar astrophysics on exoplanet evolution is currently unknown. The culmination of these system properties makes LHS 1678 a unique, compelling playground for comparative exoplanet science and understanding the formation and evolution of small, short-period exoplanets orbiting low-mass stars.

preprint2022arXiv

TOI-1075 b: A Dense, Massive, Ultra-Short Period Hot Super-Earth Straddling the Radius Gap

Populating the exoplanet mass-radius diagram in order to identify the underlying relationship that governs planet composition is driving an interdisciplinary effort within the exoplanet community. The discovery of hot super-Earths - a high temperature, short-period subset of the super-Earth planet population - has presented many unresolved questions concerning the formation, evolution, and composition of rocky planets. We report the discovery of a transiting, ultra-short period hot super-Earth orbiting TOI-1075 (TIC 351601843), a nearby ($d$ = 61.4 pc) late K-/early M-dwarf star, using data from the Transiting Exoplanet Survey Satellite (TESS). The newly discovered planet has a radius of $1.791^{+0.116}_{-0.081}$ $R_{\oplus}$, and an orbital period of 0.605 days (14.5 hours). We precisely measure the planet mass to be $9.95^{+1.36}_{-1.30}$ $M_{\oplus}$ using radial velocity measurements obtained with the Planet Finder Spectrograph (PFS), mounted on the Magellan II telescope. Our radial velocity data also show a long-term trend, suggesting an additional planet in the system. While TOI-1075 b is expected to have a substantial H/He atmosphere given its size relative to the radius gap, its high density ($9.32^{+2.05}_{-1.85}$ $\rm{g/cm^3}$) is likely inconsistent with this possibility. We explore TOI-1075 b's location relative to the M-dwarf radius valley, evaluate the planet's prospects for atmospheric characterization, and discuss potential planet formation mechanisms. Studying the TOI-1075 system in the broader context of ultra-short period planetary systems is necessary for testing planet formation and evolution theories, density enhancing mechanisms, and for future atmospheric and surface characterization studies via emission spectroscopy with JWST.

preprint2022arXiv

TOI-1696: a nearby M4 dwarf with a $3R_\oplus$ planet in the Neptunian desert

We present the discovery and validation of a temperate sub-Neptune around the nearby mid-M dwarf TIC 470381900 (TOI-1696), with a radius of $3.09 \pm 0.11 \,R_\oplus$ and an orbital period of $2.5 \,\rm{days}$, using a combination of TESS and follow-up observations using ground-based telescopes. Joint analysis of multi-band photometry from TESS, MuSCAT, MuSCAT3, Sinistro, and KeplerCam confirmed the transit signal to be achromatic as well as refined the orbital ephemeris. High-resolution imaging with Gemini/'Alopeke and high-resolution spectroscopy with the Subaru/IRD confirmed that there are no stellar companions or background sources to the star. The spectroscopic observations with IRD and IRTF/SpeX were used to determine the stellar parameters, and found the host star is an M4 dwarf with an effective temperature of $T_{eff} = 3185 \pm 76\,\rm{K}$ and a metallicity of [Fe/H] $=0.336 \pm 0.060 \,\rm{dex}$. The radial velocities measured from IRD set a $2$-$σ$ upper limit on the planetary mass to be $48.8 \,M_\oplus$. The large radius ratio ($R_p/R_\star \sim 0.1$) and the relatively bright NIR magnitude ($J=12.2 \,\rm{mag}$) make this planet an attractive target for further followup observations. TOI-1696b is one of the planets belonging to the Neptunian desert with the highest transmission spectroscopy metric discovered to date, making it an interesting candidate for atmospheric characterizations with JWST.

preprint2022arXiv

Validation of 13 Hot and Potentially Terrestrial TESS Planets

The James Webb Space Telescope (JWST) will be able to probe the atmospheres and surface properties of hot, terrestrial planets via emission spectroscopy. We identify 18 potentially terrestrial planet candidates detected by the Transiting Exoplanet Survey Satellite (TESS) that would make ideal targets for these observations. These planet candidates cover a broad range of planet radii ($R_{\rm p} \sim 0.6 - 2.0 R_\oplus$) and orbit stars of various magnitudes ($K_s = 5.78 - 10.78$, $V = 8.4 - 15.69$) and effective temperatures ($T_{\rm eff }\sim 3000 - 6000$ K). We use ground-based observations collected through the TESS Follow-up Observing Program (TFOP) and two vetting tools -- DAVE and TRICERATOPS -- to assess the reliabilities of these candidates as planets. We validate 13 planets: TOI-206 b, TOI-500 b, TOI-544 b, TOI-833 b, TOI-1075 b, TOI-1411 b, TOI-1442 b, TOI-1693 b, TOI-1860 b, TOI-2260 b, TOI-2411 b, TOI-2427 b, and TOI-2445 b. Seven of these planets (TOI-206 b, TOI-500 b, TOI-1075 b, TOI-1442 b, TOI-2260 b, TOI-2411 b, and TOI-2445 b) are ultra-short-period planets. TOI-1860 is the youngest ($133 \pm 26$ Myr) solar twin with a known planet to date. TOI-2260 is a young ($321 \pm 96$ Myr) G dwarf that is among the most metal-rich ([Fe/H] = $0.22 \pm 0.06$ dex) stars to host an ultra-short-period planet. With an estimated equilibrium temperature of $\sim 2600$ K, TOI-2260 b is also the fourth hottest known planet with $R_{\rm p} < 2 \, R_\oplus$.

preprint2022arXiv

Variability Timescales of H$α$ on Active Mid-to-Late M dwarfs

We present a study of the variation timescales of the chromospheric activity indicator H$α$ on a sample of 13 fully-convective, active mid-to-late M stars with masses between 0.1--0.3 solar masses. Our goal was to determine the dominant variability timescale and, by inference, a possible mechanism responsible for the variation. We gathered 10 or more high-resolution spectra each of 10 stars using the TRES spectrograph at times chosen to span all phases of stellar rotation, as determined from photometric data from the MEarth Observatories. All stars varied in their H$α$ emission. For 9 of these stars, we found no correlation between H$α$ and rotational phase, indicating that constant emission from fixed magnetic structures, such as starspots and plage, are unlikely to be the dominant source of H$α$ emission variability. In contrast, one star, G 7-34, shows a clear relationship between H$α$ and stellar rotational phase. Intriguingly, we found that this star is a member of the AB Doradus moving group and hence has the young age of 149 Myr. High-cadence spectroscopic observations of three additional stars revealed that they are variable on timescales ranging from 20--45 minutes, which we posit may be due to flaring behavior. For one star, GJ 1111, simultaneous TESS photometry and spectroscopic monitoring show an increase in H$α$ emission with increased photometric brightness. We conclude that low-energy flares are able to produce variation in H$α$ on the timescales we observe and thus may be the dominant source of H$α$ variability on active fully-convective M dwarfs.

preprint2021arXiv

K2-79b and K2-222b: Mass measurements of two small exoplanets with periods beyond 10 days that overlap with periodic magnetic activity signals

We present mass and radius measurements of K2-79b and K2-222b, two transiting exoplanets orbiting active G-type stars. Their respective 10.99d and 15.39d orbital periods fall near periods of signals induced by stellar magnetic activity. The two signals might therefore interfere and lead to an inaccurate estimate of exoplanet mass. We present a method to mitigate these effects when radial velocity and activity indicator observations are available over multiple observing seasons and the orbital period of the exoplanet is known. We perform correlation and periodogram analyses on sub-sets composed of each target's two observing seasons, in addition to the full data sets. For both targets, these analyses reveal an optimal season with little to no interference at the orbital period of the known exoplanet. We make a confident mass detection of each exoplanet by confirming agreement between fits to the full radial velocity set and the optimal season. For K2-79b, we measure a mass of 11.8 $\pm$ 3.6 $M_{Earth}$ and a radius of 4.09 $\pm$ 0.17 $R_{Earth}$. For K2-222b, we measure a mass of 8.0 $\pm$ 1.8 $M_{Earth}$ and a radius of 2.35 $\pm$ 0.08 $R_{Earth}$. According to model predictions, K2-79b is a highly irradiated Uranus-analog and K2-222b hosts significant amounts of water ice. We also present an RV solution for a candidate second companion orbiting K2-222 at 147.5d.

preprint2021arXiv

Precise transit and radial-velocity characterization of a resonant pair: a warm Jupiter TOI-216c and eccentric warm Neptune TOI-216b

TOI-216 hosts a pair of warm, large exoplanets discovered by the TESS Mission. These planets were found to be in or near the 2:1 resonance, and both of them exhibit transit timing variations (TTVs). Precise characterization of the planets' masses and radii, orbital properties, and resonant behavior can test theories for the origins of planets orbiting close to their stars. Previous characterization of the system using the first six sectors of TESS data suffered from a degeneracy between planet mass and orbital eccentricity. Radial velocity measurements using HARPS, FEROS, and PFS break that degeneracy, and an expanded TTV baseline from TESS and an ongoing ground-based transit observing campaign increase the precision of the mass and eccentricity measurements. We determine that TOI-216c is a warm Jupiter, TOI-216b is an eccentric warm Neptune, and that they librate in the 2:1 resonance with a moderate libration amplitude of 60 +/- 2 degrees; small but significant free eccentricity of 0.0222 +0.0005/-0.0003 for TOI-216b; and small but significant mutual inclination of 1.2-3.9 degrees (95% confidence interval). The libration amplitude, free eccentricity, and mutual inclination imply a disturbance of TOI-216b before or after resonance capture, perhaps by an undetected third planet.

preprint2021arXiv

TESS discovery of a super-Earth and three sub-Neptunes hosted by the bright, Sun-like star HD 108236

We report the discovery and validation of four extrasolar planets hosted by the nearby, bright, Sun-like (G3V) star HD~108236 using data from the Transiting Exoplanet Survey Satellite (TESS). We present transit photometry, reconnaissance and precise Doppler spectroscopy as well as high-resolution imaging, to validate the planetary nature of the objects transiting HD~108236, also known as the TESS Object of Interest (TOI) 1233. The innermost planet is a possibly-rocky super-Earth with a period of $3.79523_{-0.00044}^{+0.00047}$ days and has a radius of $1.586\pm0.098$ $R_\oplus$. The outer planets are sub-Neptunes, with potential gaseous envelopes, having radii of $2.068_{-0.091}^{+0.10}$ $R_\oplus$, $2.72\pm0.11$ $R_\oplus$, and $3.12_{-0.12}^{+0.13}$ $R_\oplus$ and periods of $6.20370_{-0.00052}^{+0.00064}$ days, $14.17555_{-0.0011}^{+0.00099}$ days, and $19.5917_{-0.0020}^{+0.0022}$ days, respectively. With V and K$_{\rm s}$ magnitudes of 9.2 and 7.6, respectively, the bright host star makes the transiting planets favorable targets for mass measurements and, potentially, for atmospheric characterization via transmission spectroscopy. HD~108236 is the brightest Sun-like star in the visual (V) band known to host four or more transiting exoplanets. The discovered planets span a broad range of planetary radii and equilibrium temperatures, and share a common history of insolation from a Sun-like star ($R_\star = 0.888 \pm 0.017$ R$_\odot$, $T_{\rm eff} = 5730 \pm 50$ K), making HD 108236 an exciting, opportune cosmic laboratory for testing models of planet formation and evolution.

preprint2021arXiv

The Featureless HST/WFC3 Transmission Spectrum of the Rocky Exoplanet GJ 1132b: No Evidence For A Cloud-Free Primordial Atmosphere and Constraints on Starspot Contamination

Orbiting a M dwarf 12 pc away, the transiting exoplanet GJ 1132b is a prime target for transmission spectroscopy. With a mass of 1.7 Earth masses and radius of 1.1 Earth radii, GJ 1132b's bulk density indicates that this planet is rocky. Yet with an equilibrium temperature of 580 K, GJ 1132b may still retain some semblance of an atmosphere. Understanding whether this atmosphere exists and its composition will be vital for understanding how the atmospheres of terrestrial planets orbiting M dwarfs evolve. We observe five transits of GJ 1132b with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST). We find a featureless transmission spectrum from 1.1--1.7 microns, ruling out cloud-free atmospheres with metallicities <300x Solar with >4.8$σ$ confidence. We combine our WFC3 results with transit depths from TESS and archival broadband and spectroscopic observations to find a featureless spectrum from 0.7--4.5 microns. GJ 1132b has a high mean molecular weight atmosphere, possesses a high-altitude aerosol layer, or has effectively no atmosphere. Higher precision observations are required to differentiate between these possibilities. We explore the impact of hot and cold starspots on the observed transmission spectrum GJ 1132b, quantifying the amplitude of spot-induced transit depth features. Using a simple Poisson model we estimate spot temperature contrasts, spot covering fractions, and spot sizes for GJ 1132. These limits, and the modeling framework, may be useful for future observations of GJ 1132b or other planets transiting similarly inactive M dwarfs.

preprint2021arXiv

TOI-712: a system of adolescent mini-Neptunes extending to the habitable zone

As an all-sky survey, NASA's $TESS$ mission is able to detect the brightest and rarest types of transiting planetary systems, including young planets that enable study of the evolutionary processes that occur within the first billion years. Here, we report the discovery of a young, multi-planet system orbiting the bright K4.5V star, TOI-712 ($V = 10.838$, $M_\star = 0.733_{-0.025}^{+0.026} M_\odot$, $R_\star = 0.674\pm0.016 R_\odot$, $T_{\rm eff} = 4622_{-60}^{+61}$ K). From the $TESS$ light curve, we measure a rotation period of 12.48 days, and derive an age between about $500$ Myr and 1.1 Gyr. The photometric observations reveal three transiting mini-Neptunes ($R_b = 2.049^{+0.12}_{-0.080} R_\oplus$, $R_c = 2.701^{+0.092}_{-0.082} R_\oplus$, $R_d = 2.474^{+0.090}_{-0.082} R_\oplus $), with orbital periods of $P_b = 9.531$ days, $P_c = 51.699$ days, and $P_d = 84.839$ days. After modeling the three-planet system, an additional Earth-sized candidate is identified, TOI-712.05 ($P = 4.32$ days, $R_P = 0.81 \pm 0.11 R_\oplus$). We calculate that the habitable zone falls between 0.339 and 0.844 au (82.7 and 325.3 days), placing TOI-712 d near its inner edge. Among planetary systems harboring temperate planets, TOI-712 ($T = 9.9$) stands out as a relatively young star bright enough to motivate further characterization.

preprint2021arXiv

TOI-811b and TOI-852b: New transiting brown dwarfs with similar masses and very different radii and ages from the TESS mission

We report the discovery of two transiting brown dwarfs (BDs), TOI-811b and TOI-852b, from NASA's Transiting Exoplanet Survey Satellite mission. These two transiting BDs have similar masses, but very different radii and ages. Their host stars have similar masses, effective temperatures, and metallicities. The younger and larger transiting BD is TOI-811b at a mass of $M_b = 55.3 \pm 3.2{\rm M_J}$ and radius of $R_b = 1.35 \pm 0.09{\rm R_J}$ and it orbits its host star in a period of $P = 25.16551 \pm 0.00004$ days. Its age of $93^{+61}_{-29}$ Myr, which we derive from an application of gyrochronology to its host star, is why this BD's radius is relatively large, not heating from its host star since this BD orbits at a longer orbital period than most known transiting BDs. This constraint on the youth of TOI-811b allows us to test substellar mass-radius isochrones where the radius of BDs changes rapidly with age. TOI-852b is a much older (4.0 Gyr from stellar isochrone models of the host star) and smaller transiting BD at a mass of $M_b = 53.7 \pm 1.3{\rm M_J}$, a radius of $R_b = 0.75 \pm 0.03{\rm R_J}$, and an orbital period of $P = 4.94561 \pm 0.00008$ days. TOI-852b joins the likes of other old transiting BDs that trace out the oldest substellar mass-radius isochrones where contraction of the BD's radius asymptotically slows. Both host stars have a mass of $M_\star = 1.32{\rm M_\odot}\pm0.05$ and differ in their radii, $T_{\rm eff}$, and [Fe/H] with TOI-811 having $R_\star=1.27\pm0.09{\rm R_\odot}$, $T_{\rm eff} = 6107 \pm 77$K, and $\rm [Fe/H] = +0.40 \pm 0.09$ and TOI-852 having $R_\star=1.71\pm0.04{\rm R_\odot}$, $T_{\rm eff} = 5768 \pm 84$K, and $\rm [Fe/H] = +0.33 \pm 0.09$. We take this opportunity to examine how TOI-811b and TOI-852b serve as test points for young and old substellar isochrones, respectively.

preprint2020arXiv

A pair of TESS planets spanning the radius valley around the nearby mid-M dwarf LTT 3780

We present the confirmation of two new planets transiting the nearby mid-M dwarf LTT 3780 (TIC 36724087, TOI-732, $V=13.07$, $K_s=8.204$, $R_s$=0.374 R$_{\odot}$, $M_s$=0.401 M$_{\odot}$, d=22 pc). The two planet candidates are identified in a single TESS sector and are validated with reconnaissance spectroscopy, ground-based photometric follow-up, and high-resolution imaging. With measured orbital periods of $P_b=0.77$ days, $P_c=12.25$ days and sizes $r_{p,b}=1.33\pm 0.07$ R$_{\oplus}$, $r_{p,c}=2.30\pm 0.16$ R$_{\oplus}$, the two planets span the radius valley in period-radius space around low mass stars thus making the system a laboratory to test competing theories of the emergence of the radius valley in that stellar mass regime. By combining 63 precise radial-velocity measurements from HARPS and HARPS-N, we measure planet masses of $m_{p,b}=2.62^{+0.48}_{-0.46}$ M$_{\oplus}$ and $m_{p,c}=8.6^{+1.6}_{-1.3}$ M$_{\oplus}$, which indicates that LTT 3780b has a bulk composition consistent with being Earth-like, while LTT 3780c likely hosts an extended H/He envelope. We show that the recovered planetary masses are consistent with predictions from both photoevaporation and from core-powered mass loss models. The brightness and small size of LTT 3780, along with the measured planetary parameters, render LTT 3780b and c as accessible targets for atmospheric characterization of planets within the same planetary system and spanning the radius valley.

preprint2020arXiv

Flare Rates, Rotation Periods. and Spectroscopic Activity Indicators of a Volume-Complete Sample of Mid-to-Late M dwarfs within 15 Parsecs

We present a study of flare rates, rotation periods, and spectroscopic activity indicators of 125 single stars within 15 parsecs and with masses between 0.1$-$0.3 $M_\odot$ observed during the first year of the TESS mission, with the goal of elucidating the relationship between these various magnetically connected phenomena. We gathered multi-epoch high resolution spectra of each target and we measured equivalent widths of the activity indicators Helium I D$_3$, $Hα$, and the Calcium infrared triplet line at 8542.09 angstroms. We present 18 new rotation periods from MEarth photometry and 19 new rotation periods from TESS photometry. We present a catalog of 1392 flares. After correcting for sensitivity, we find the slope of the flare frequency distribution for all stars to have a standard value of $α$ = 1.98 $\pm$ 0.02. We determine R$_{31.5}$, the rate of flares per day with energies above E = 3.16$\times$10$^{31}$ ergs in the TESS bandpass. We find that below a critical value of $Hα$ EW = -0.71 angstroms, log R$_{31.5}$ increases linearly with increasing $Hα$ emission; above this value, log R$_{31.5}$ declines rapidly. The stars divide into two groups: 26% have $Hα$ in emission, high flare rates with typical values of log R$_{31.5}$ = -1.30 $\pm$ 0.08, and have Rossby numbers $<$ 0.50. The remaining 74% show little to no $Hα$ in emission and exhibit log R$_{31.5}$ $<$ -3.86, with the majority of these stars not showing a single flare during the TESS observations.

preprint2020arXiv

GJ 1252 b: A 1.2 $R_{\oplus}$ planet transiting an M3-dwarf at 20.4 pc

We report the discovery of GJ 1252 b, a planet with a radius of 1.193 $\pm$ 0.074 $R_{\oplus}$ and an orbital period of 0.52 days around an M3-type star (0.381 $\pm$ 0.019 $M_{\odot}$, 0.391 $\pm$ 0.020 $R_{\odot}$) located 20.385 $\pm$ 0.019 pc away. We use TESS data, ground-based photometry and spectroscopy, Gaia astrometry, and high angular resolution imaging to show that the transit signal seen in the TESS data must originate from a transiting planet. We do so by ruling out all false positive scenarios that attempt to explain the transit signal as originating from an eclipsing stellar binary. Precise Doppler monitoring also leads to a tentative mass measurement of 2.09 $\pm$ 0.56 $M_{\oplus}$. The host star proximity, brightness ($V$ = 12.19 mag, $K$ = 7.92 mag), low stellar activity, and the system's short orbital period make this planet an attractive target for detailed characterization, including precise mass measurement, looking for other objects in the system, and planet atmosphere characterization.

preprint2020arXiv

Simultaneous Optical Transmission Spectroscopy of a Terrestrial, Habitable-Zone Exoplanet with Two Ground-Based Multi-Object Spectrographs

Investigating the atmospheres of rocky exoplanets is key to performing comparative planetology between these worlds and the terrestrial planets that reside in the inner solar system. Terrestrial exoplanet atmospheres exhibit weak signals, and attempting to detect them pushes at the boundaries of what is possible for current instrumentation. We focus on the habitable-zone terrestrial exoplanet LHS 1140b. Given its 25-day orbital period and 2 hr transit duration, capturing transits of LHS 1140b is challenging. We observed two transits of this object, approximately 1 yr apart, which yielded four data sets thanks to our simultaneous use of the IMACS and LDSS3C multiobject spectrographs mounted on the twin Magellan telescopes at Las Campanas Observatory. We present a jointly fit white light curve, as well as jointly fit 20 nm wavelength-binned light curves from which we construct a transmission spectrum. Binning the joint white light-curve residuals to 3-minute time bins gives an rms of 145 ppm; binning down to 10-minute time bins gives an rms of 77 ppm. Our median uncertainty in Rp^2/Rs^2 in the 20 nm wavelength bins is 260 ppm, and we achieve an average precision of 1.3x the photon noise when fitting the wavelength-binned light curves with a Gaussian process regression. Our precision on Rp^2/Rs^2 is a factor of four larger than the feature amplitudes of a clear, hydrogen-dominated atmosphere, meaning that we are not able to test realistic models of LHS 1140b's atmosphere. The techniques and caveats presented here are applicable to the growing sample of terrestrial worlds in the Transiting Exoplanet Survey Satellite era, as well as to the upcoming generation of ground-based giant segmented mirror telescopes.

preprint2020arXiv

Spectroscopic Orbits of Eleven Nearby, Mid-to-Late M Dwarf Binaries

We present the spectroscopic orbits of eleven nearby, mid-to-late M dwarf binary systems in a variety of configurations: two single-lined binaries (SB1s), seven double-lined binaries (SB2s), one double-lined triple (ST2), and one triple-lined triple (ST3). Eight of these orbits are the first published for these systems, while five are newly identified multiples. We obtained multi-epoch, high-resolution spectra with the TRES instrument on the 1.5m Tillinghast Reflector at the Fred Lawrence Whipple Observatory located on Mt. Hopkins in AZ. Using the TiO molecular bands at 7065 -- 7165 Angstroms, we calculated radial velocities for these systems, from which we derived their orbits. We find LHS 1817 to have in a 7-hour period a companion that is likely a white dwarf, due to the ellipsoidal modulation we see in our MEarth-North light curve data. We find G 123-45 and LTT 11586 to host companions with minimum masses of 41 M_Jup and 44 M_Jup with orbital periods of 35 and 15 days, respectively. We find 2MA 0930+0227 to have a rapidly rotating stellar companion in a 917-day orbital period. GJ 268, GJ 1029, LP 734-34, GJ 1182, G 258-17, and LTT 7077 are SB2s with stellar companions with orbital periods of 10, 96, 34, 154, 5, and 84 days; LP 655-43 is an ST3 with one companion in an 18-day orbital period and an outer component in a longer undetermined period. In addition, we present radial velocities for both components of L 870-44AB and for the outer components of LTT 11586 and LP 655-43.

preprint2020arXiv

The First Habitable Zone Earth-Sized Planet From TESS II: $Spitzer$ Confirms TOI-700 d

We present $Spitzer$ 4.5$μ$m observations of the transit of TOI-700 d, a habitable zone Earth-sized planet in a multiplanet system transiting a nearby M-dwarf star (TIC 150428135, 2MASS J06282325-6534456). TOI-700 d has a radius of $1.144^{+0.062}_{-0.061}R_\oplus$ and orbits within its host star's conservative habitable zone with a period of 37.42 days ($T_\mathrm{eq} \sim 269$K). TOI-700 also hosts two small inner planets (R$_b$=$1.037^{+0.065}_{-0.064}R_\oplus$ & R$_c$=$2.65^{+0.16}_{-0.15}R_\oplus$) with periods of 9.98 and 16.05 days, respectively. Our $Spitzer$ observations confirm the TESS detection of TOI-700 d and remove any remaining doubt that it is a genuine planet. We analyze the $Spitzer$ light curve combined with the 11 sectors of TESS observations and a transit of TOI-700 c from the LCOGT network to determine the full system parameters. Although studying the atmosphere of TOI-700 d is not likely feasible with upcoming facilities, it may be possible to measure the mass of TOI-700 d using state-of-the-art radial velocity instruments (expected RV semi-amplitude of $\sim$70 cm/s).

preprint2020arXiv

The Volume-Complete Sample of M Dwarfs with Masses 0.1 <= M/M_sol <= 0.3 within 15 Parsecs

M dwarfs with masses 0.1 <= M/M_sol <= 0.3 are under increasing scrutiny because these fully convective stars pose interesting astrophysical questions regarding their magnetic activity and angular momentum history. They also afford the most accessible near-future opportunity to study the atmospheres of terrestrial planets. Because they are intrinsically low in luminosity, the identification of the nearest examples of these M dwarfs is essential for progress. We present the volume-complete, all-sky list of 512 M dwarfs with masses 0.1 <= M/M_sol <= 0.3 and with trigonometric distances placing them within 15 pc (parallax >= 66.67 mas) from which we have created a sample of 413 M dwarfs for spectroscopic study. We present the mass function for these 512 M dwarfs, which increases with decreasing stellar mass in linear mass space, but is flat in logarithmic mass space. As part of this sample, we present new VRI photometry for 17 targets, measured as a result of the RECONS group's long-term work at the CTIO/SMARTS 0.9m telescope. We also note the details of targets that are known to be members of multiple systems and find a preliminary multiplicity rate of 21 +/- 2% for the primary M dwarfs in our sample, when considering known stellar and brown dwarf companions at all separations from their primaries. We further find that 43 +/- 2% of all M dwarfs with masses 0.1 <= M/M_sol <= 0.3 are found in multiple systems with primary stars of all masses within 15 pc.

preprint2020arXiv

TIC 278956474: Two close binaries in one young quadruple system, identified by \textit{TESS}

We have identified a quadruple system with two close eclipsing binaries in TESS data. The object is unresolved in Gaia and appears as a single source at parallax 1.08~$\pm$0.01 mas. Both binaries have observable primary and secondary eclipses and were monitored throughout TESS Cycle 1 (sectors 1-13), falling within the TESS Continuous Viewing Zone. In one eclipsing binary (P = 5.488 d), the smaller star is completely occluded by the larger star during the secondary eclipse; in the other (P = 5.674 d) both eclipses are grazing. Using these data, spectroscopy, speckle photometry, SED analysis and evolutionary stellar tracks, we have constrained the masses and radii of the four stars in the two eclipsing binaries. The Li I EW indicates an age of 10-50 Myr and, with an outer period of $858^{+7}_{-5}$ days, our analysis indicates this is one of the most compact young 2+2 quadruple systems known.

preprint2020arXiv

TOI-1235 b: a keystone super-Earth for testing radius valley emergence models around early M dwarfs

Small planets on close-in orbits tend to exhibit envelope mass fractions of either effectively zero or up to a few percent depending on their size and orbital period. Models of thermally-driven atmospheric mass loss and of terrestrial planet formation in a gas-poor environment make distinct predictions regarding the location of this rocky/non-rocky transition in period-radius space. Here we present the confirmation of TOI-1235 b ($P=3.44$ days, $r_p=1.738^{+0.087}_{-0.076}$ R$_{\oplus}$), a planet whose size and period are intermediate between the competing model predictions thus making the system an important test case for emergence models of the rocky/non-rocky transition around early M dwarfs ($R_s=0.630\pm 0.015$ R$_{\odot}$, $M_s=0.640\pm 0.016$ M$_{\odot}$). We confirm the TESS planet discovery using reconnaissance spectroscopy, ground-based photometry, high-resolution imaging, and a set of 38 precise radial-velocities from HARPS-N and HIRES. We measure a planet mass of $6.91^{+0.75}_{-0.85}$ M$_{\oplus}$, which implies an iron core mass fraction of $20^{+15}_{-12}$% in the absence of a gaseous envelope. The bulk composition of TOI-1235 b is therefore consistent with being Earth-like and we constrain a H/He envelope mass fraction to be $<0.5$% at 90% confidence. Our results are consistent with model predictions from thermally-driven atmospheric mass loss but not with gas-poor formation, suggesting that the former class of processes remain efficient at sculpting close-in planets around early M dwarfs. Our RV analysis also reveals a strong periodicity close to the first harmonic of the photometrically-determined stellar rotation period that we treat as stellar activity, despite other lines of evidence favoring a planetary origin ($P=21.8^{+0.9}_{-0.8}$ days, $m_p\sin{i}=13.0^{+3.8}_{-5.3}$ M$_{\oplus}$) that cannot be firmly ruled out by our data.

preprint2020arXiv

Zodiacal Exoplanets in Time (ZEIT) IX: a flat transmission spectrum and a highly eccentric orbit for the young Neptune K2-25b as revealed by Spitzer

Transiting planets in nearby young clusters offer the opportunity to study the atmospheres and dynamics of planets during their formative years. To this end, we focused on K2-25b -- a close-in ($P$=3.48 days), Neptune-sized exoplanet orbiting a M4.5 dwarf in the 650Myr Hyades cluster. We combined photometric observations of K2-25 covering a total of 44 transits and spanning >2 yr, drawn from a mix of space-based telescopes (Spitzer Space Telescope and K2) and ground-based facilities (Las Cumbres Observatory Global Telescope network and MEarth). The transit photometry spanned 0.6--4.5$μ$m, which enabled our study of K2-25b's transmission spectrum. We combined and fit each dataset at a common wavelength within a Markov Chain Monte Carlo framework, yielding consistent planet parameters. The resulting transit depths ruled out a solar-composition atmosphere for K2-25b for the range of expected planetary masses and equilibrium temperature at a $>4σ$ confidence level, and are consistent with a flat transmission spectrum. Mass constraints and transit observations at a finer grid of wavelengths (e.g., from the Hubble Space Telescope) are needed to make more definitive statements about the presence of clouds or an atmosphere of high mean molecular weight. Our precise measurements of K2-25b's transit duration also enabled new constraints on the eccentricity of K2-25's orbit. We find K2-25b's orbit to be eccentric ($e>0.20$) for all reasonable stellar densities and independent of the observation wavelength or instrument. The high eccentricity is suggestive of a complex dynamical history and motivates future searches for additional planets or stellar companions.

preprint2019arXiv

The young planetary system K2-25: constraints on companions and starspots

The abundance of planets with orbital periods of a few to tens of days suggests that exoplanets experience complex dynamical histories. Planets in young stellar clusters or associations have well-constrained ages and therefore provide an opportunity to explore the dynamical evolution of exoplanets. K2-25b is a Neptune-sized planet in an eccentric, 3.48 day orbit around an M4.5 dwarf star in the Hyades cluster (650 Myr). In order to investigate its non-zero eccentricity and tight orbit, we analyze transit timing variations (TTVs) which could reveal clues to the migration processes that may have acted on the planet. We obtain 12 non-consecutive transits using the MEarth Observatories and long-term photometric monitoring, which we combine with 10 transits from the Spitzer Space Telescope and 20 transits from K2. Tables of MEarth photometry accompany this work. We fit each transit lightcurve independently. We first investigate whether inhomogeneities on the stellar surface (such as spots or plages) are differentially affecting our transit observations. The measured transit depth does not vary significantly between transits, though we see some deviations from the fiducial transit model. We then looked for TTVs as evidence of a non-transiting perturber in the system. We find no evidence for > 1 $M_\oplus$ companions within a 2:1 period ratio, or for > 5 $M_\oplus$ planets within a 7:2 period ratio.

preprint2016arXiv

A 1.9 Earth radius rocky planet and the discovery of a non-transiting planet in the Kepler-20 system

Kepler-20 is a solar-type star (V = 12.5) hosting a compact system of five transiting planets, all packed within the orbital distance of Mercury in our own Solar System. A transition from rocky to gaseous planets with a planetary transition radius of ~1.6 REarth has recently been proposed by several publications in the literature (Rogers 2015;Weiss & Marcy 2014). Kepler-20b (Rp ~ 1.9 REarth) has a size beyond this transition radius, however previous mass measurements were not sufficiently precise to allow definite conclusions to be drawn regarding its composition. We present new mass measurements of three of the planets in the Kepler-20 system facilitated by 104 radial velocity measurements from the HARPS-N spectrograph and 30 archival Keck/HIRES observations, as well as an updated photometric analysis of the Kepler data and an asteroseismic analysis of the host star (MStar = 0.948+-0.051 Msun and Rstar = 0.964+-0.018 Rsun). Kepler-20b is a 1.868+0.066-0.034 REarth planet in a 3.7 day period with a mass of 9.70+1.41-1.44 MEarth resulting in a mean density of 8.2+1.5-1.3 g/cc indicating a rocky composition with an iron to silicate ratio consistent with that of the Earth. This makes Kepler-20b the most massive planet with a rocky composition found to date. Furthermore, we report the discovery of an additional non-transiting planet with a minimum mass of 19.96+3.08-3.61 MEarth and an orbital period of ~34 days in the gap between Kepler-20f (P ~ 11 days) and Kepler-20d (P ~ 78 days).

preprint2016arXiv

Kepler-21b: A rocky planet around a V = 8.25 magnitude star

HD 179070, aka Kepler-21, is a V = 8.25 F6IV star and the brightest exoplanet host discovered by Kepler. An early detailed analysis by Howell et al. (2012) of the first thirteen months (Q0 - Q5) of Kepler light curves revealed transits of a planetary companion, Kepler-21b, with a radius of about 1.60 +/- 0.04 R_earth and an orbital period of about 2.7857 days. However, they could not determine the mass of the planet from the initial radial velocity observations with Keck-HIRES, and were only able to impose a 2-sigma upper limit of 10 M_earth. Here we present results from the analysis of 82 new radial velocity observations of this system obtained with HARPS-N, together with the existing 14 HIRES data points. We detect the Doppler signal of Kepler-21b with a radial velocity semi-amplitude K = 2.00 +/- 0.65 m/s, which corresponds to a planetary mass of 5.1 +/- 1.7 M_earth. We also measure an improved radius for the planet of 1.639 (+0.019, -0.015) R_earth, in agreement with the radius reported by Howell et al. (2012). We conclude that Kepler-21b, with a density of 6.4 +/- 2.1 g/cm^3, belongs to the population of terrestrial planets with iron, magnesium silicate interiors, which have lost the majority of their envelope volatiles via stellar winds or gravitational escape. The radial velocity analysis presented in this paper serves as example of the type of analysis that will be necessary to confirm the masses of TESS small planet candidates.

preprint2016arXiv

Results from the 2014 November 15th multi-chord stellar occultation by the TNO (229762) 2007 UK$_{126}$

We present results derived from the first multi-chord stellar occultation by the trans-Neptunian object (229762) 2007 UK$_{126}$, observed on 2014 November 15. The event was observed by the Research and Education Collaborative Occultation Network (RECON) project and International Occultation Timing Association (IOTA) collaborators throughout the United States. Use of two different data analysis methods obtain a satisfactory fit to seven chords, yelding an elliptical fit to the chords with an equatorial radius of $R=338_{-10} ^{+15}$ km and equivalent radius of $R_{eq}=319_{-7} ^{+14}$ km. A circular fit also gives a radius of $R=324_{-23} ^{+30}$ km. Assuming that the object is a Maclaurin spheroid with indeterminate aspect angle, and using two published absolute magnitudes for the body, we derive possible ranges for geometric albedo between $p_{V}=0.159_{-0.013} ^{+0.007}$ and $p_{R}=0.189_{-0.015}^{+0.009}$, and for the body oblateness between $ε=0.105_{-0.040} ^{+0.050}$ and $ε=0.118_{-0.048} ^{+0.055}$. For a nominal rotational period of 11.05 h, an upper limit for density of $ρ=1740$ kg~m$^{-3}$ is estimated for the body.

preprint2016arXiv

The Impact of Stellar Rotation on the Detectability of Habitable Planets Around M Dwarfs

Stellar activity and rotation frustrate the detection of exoplanets through the radial velocity technique. This effect is particularly of concern for M dwarfs, which can remain magnetically active for billions of years. We compile rotation periods for late-type stars and for the M dwarf planet-host sample in order to investigate the rotation periods of older field stars across the main sequence. We show that for stars with masses between 0.25 and 0.5 solar masses (M4V to M1V), the stellar rotation period typical of field stars coincides with the orbital periods of planets in the habitable zone. This will pose a fundamental challenge to the discovery and characterization of potentially habitable planets around early M dwarfs. Due to the longer rotation periods reached by mid M dwarfs and the shorter orbital period at which the planetary habitable zone is found, stars with masses between 0.1 and 0.25 solar masses (M6V to M4V) offer better opportunities for the detection of habitable planets via radial velocities.

preprint2016arXiv

The Need for Laboratory Work to Aid in The Understanding of Exoplanetary Atmospheres

Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr)

preprint2016arXiv

The Orbit and Mass of the Third Planet in the Kepler-56 System

While the vast majority of multiple-planet systems have their orbital angular momentum axes aligned with the spin axis of their host star, Kepler-56 is an exception: its two transiting planets are coplanar yet misaligned by at least 40 degrees with respect to their host star. Additional follow-up observations of Kepler-56 suggest the presence of a massive, non-transiting companion that may help explain this misalignment. We model the transit data along with Keck/HIRES and HARPS-N radial velocity data to update the masses of the two transiting planets and infer the physical properties of the third, non-transiting planet. We employ a Markov Chain Monte Carlo sampler to calculate the best-fitting orbital parameters and their uncertainties for each planet. We find the outer planet has a period of 1002 $\pm$ 5 days and minimum mass of 5.61 $\pm$ 0.38 Jupiter masses. We also place a 95% upper limit of 0.80 m/s/yr on long-term trends caused by additional, more distant companions.

preprint2016arXiv

The rotation and Galactic kinematics of mid M dwarfs in the Solar Neighborhood

Rotation is a directly-observable stellar property, and drives magnetic field generation and activity through a magnetic dynamo. Main sequence stars with masses below approximately 0.35Msun (mid-to-late M dwarfs) are fully-convective, and are expected to have a different type of dynamo mechanism than solar-type stars. Measurements of their rotation rates provide insights into these mechanisms, but few rotation periods are available for these stars at field ages. Using photometry from the MEarth transit survey, we measure rotation periods for 387 nearby, mid-to-late M dwarfs in the Northern hemisphere, finding periods from 0.1 to 140 days. The typical detected rotator has stable, sinusoidal photometric modulations at a semi-amplitude of 0.5 to 1%. We find no period-amplitude relation for stars below 0.25Msun and an anti-correlation between period and amplitude for higher-mass M dwarfs. We highlight the existence of older, slowly-rotating stars without Hα emission that nevertheless have strong photometric variability. The Galactic kinematics of our sample is consistent with the local population of G and K dwarfs, and rotators have metallicities characteristic of the Solar Neighborhood. We use the W space velocities and established age-velocity relations to estimate that stars with P<10 days are on average <2 Gyrs, and that those with P>70 days are about 5 Gyrs. The period distribution is mass dependent: as the mass decreases, the slowest rotators at a given mass have longer periods, and the fastest rotators have shorter periods. We find a lack of stars with intermediate rotation periods. [Abridged]

preprint2016arXiv

Transiting Exoplanet Studies and Community Targets for JWST's Early Release Science Program

The James Webb Space Telescope will revolutionize transiting exoplanet atmospheric science due to its capability for continuous, long-duration observations and its larger collecting area, spectral coverage, and spectral resolution compared to existing space-based facilities. However, it is unclear precisely how well JWST will perform and which of its myriad instruments and observing modes will be best suited for transiting exoplanet studies. In this article, we describe a prefatory JWST Early Release Science (ERS) program that focuses on testing specific observing modes to quickly give the community the data and experience it needs to plan more efficient and successful future transiting exoplanet characterization programs. We propose a multi-pronged approach wherein one aspect of the program focuses on observing transits of a single target with all of the recommended observing modes to identify and understand potential systematics, compare transmission spectra at overlapping and neighboring wavelength regions, confirm throughputs, and determine overall performances. In our search for transiting exoplanets that are well suited to achieving these goals, we identify 12 objects (dubbed "community targets") that meet our defined criteria. Currently, the most favorable target is WASP-62b because of its large predicted signal size, relatively bright host star, and location in JWST's continuous viewing zone. Since most of the community targets do not have well-characterized atmospheres, we recommend initiating preparatory observing programs to determine the presence of obscuring clouds/hazes within their atmospheres. Measurable spectroscopic features are needed to establish the optimal resolution and wavelength regions for exoplanet characterization. Other initiatives from our proposed ERS program include testing the instrument brightness limits and performing phase-curve observations.(Abridged)

preprint2016arXiv

Zodiacal Exoplanets in Time (ZEIT) III: A short-period planet orbiting a pre-main-sequence star in the Upper Scorpius OB Association

We confirm and characterize a close-in ($P_{\rm{orb}}$ = 5.425 days), super-Neptune sized ($5.04^{+0.34}_{-0.37}$ Earth radii) planet transiting K2-33 (2MASS J16101473-1919095), a late-type (M3) pre-main sequence (11 Myr-old) star in the Upper Scorpius subgroup of the Scorpius-Centaurus OB association. The host star has the kinematics of a member of the Upper Scorpius OB association, and its spectrum contains lithium absorption, an unambiguous sign of youth (<20 Myr) in late-type dwarfs. We combine photometry from K2 and the ground-based MEarth project to refine the planet's properties and constrain the host star's density. We determine \name's bolometric flux and effective temperature from moderate resolution spectra. By utilizing isochrones that include the effects of magnetic fields, we derive a precise radius (6-7%) and mass (16%) for the host star, and a stellar age consistent with the established value for Upper Scorpius. Follow-up high-resolution imaging and Doppler spectroscopy confirm that the transiting object is not a stellar companion or a background eclipsing binary blended with the target. The shape of the transit, the constancy of the transit depth and periodicity over 1.5 years, and the independence with wavelength rules out stellar variability, or a dust cloud or debris disk partially occulting the star as the source of the signal; we conclude it must instead be planetary in origin. The existence of K2-33b suggests close-in planets can form in situ or migrate within $\sim 10$ Myr, e.g., via interactions with a disk, and that long-timescale dynamical migration such as by Lidov-Kozai or planet-planet scattering is not responsible for all short-period planets.

preprint2015arXiv

A disintegrating minor planet transiting a white dwarf

White dwarfs are the end state of most stars, including the Sun, after they exhaust their nuclear fuel. Between 1/4 and 1/2 of white dwarfs have elements heavier than helium in their atmospheres, even though these elements should rapidly settle into the stellar interiors unless they are occasionally replenished. The abundance ratios of heavy elements in white dwarf atmospheres are similar to rocky bodies in the Solar system. This and the existence of warm dusty debris disks around about 4% of white dwarfs suggest that rocky debris from white dwarf progenitors' planetary systems occasionally pollute the stars' atmospheres. The total accreted mass can be comparable to that of large asteroids in the solar system. However, the process of disrupting planetary material has not yet been observed. Here, we report observations of a white dwarf being transited by at least one and likely multiple disintegrating planetesimals with periods ranging from 4.5 hours to 4.9 hours. The strongest transit signals occur every 4.5 hours and exhibit varying depths up to 40% and asymmetric profiles, indicative of a small object with a cometary tail of dusty effluent material. The star hosts a dusty debris disk and the star's spectrum shows prominent lines from heavy elements like magnesium, aluminium, silicon, calcium, iron, and nickel. This system provides evidence that heavy element pollution of white dwarfs can originate from disrupted rocky bodies such as asteroids and minor planets.

preprint2015arXiv

A rocky planet transiting a nearby low-mass star

M-dwarf stars -- hydrogen-burning stars that are smaller than 60 per cent of the size of the Sun -- are the most common class of star in our Galaxy and outnumber Sun-like stars by a ratio of 12:1. Recent results have shown that M dwarfs host Earth-sized planets in great numbers: the average number of M-dwarf planets that are between 0.5 to 1.5 times the size of Earth is at least 1.4 per star. The nearest such planets known to transit their star are 39 parsecs away, too distant for detailed follow-up observations to measure the planetary masses or to study their atmospheres. Here we report observations of GJ 1132b, a planet with a size of 1.2 Earth radii that is transiting a small star 12 parsecs away. Our Doppler mass measurement of GJ 1132b yields a density consistent with an Earth-like bulk composition, similar to the compositions of the six known exoplanets with masses less than six times that of the Earth and precisely measured densities. Receiving 19 times more stellar radiation than the Earth, the planet is too hot to be habitable but is cool enough to support a substantial atmosphere, one that has probably been considerably depleted of hydrogen. Because the host star is nearby and only 21 per cent the radius of the Sun, existing and upcoming telescopes will be able to observe the composition and dynamics of the planetary atmosphere.

preprint2015arXiv

An Activity-Rotation Relationship and Kinematic Analysis of Nearby Mid-to-Late-type M Dwarfs

Using spectroscopic observations and photometric light curves of 238 nearby M dwarfs from the MEarth exoplanet transit survey, we examine the relationships between magnetic activity (quantified by H-alpha emission), rotation period, and stellar age. Previous attempts to investigate the relationship between magnetic activity and rotation in these stars were hampered by the limited number of M dwarfs with measured rotation periods (and the fact that vsini measurements probe only rapid rotation). However, the photometric data from MEarth allows us to probe a wide range of rotation periods for hundreds of M dwarf stars (from shorter than than one to longer than 100 days). Over all M spectral types that we probe, we find that the presence of magnetic activity is tied to rotation, including for late-type, fully convective M dwarfs. We also find evidence that the fraction of late-type M dwarfs that are active may be higher at longer rotation periods compared to their early-type counterparts, with several active, late-type, slowly rotating stars present in our sample. Additionally, we find that all M dwarfs with rotation periods shorter than 26 days (early-type; M1-M4) and 86 days (late-type; M5-M8) are magnetically active. This potential mismatch suggests that the physical mechanisms that connect stellar rotation to chromospheric heating may be different in fully convective stars. A kinematic analysis suggests that the magnetically active, rapidly rotating stars are consistent with a kinematically young population, while slow-rotators are less active or inactive and appear to belong to an older, dynamically heated stellar population.

preprint2015arXiv

An empirical calibration to estimate cool dwarf fundamental parameters from H-band spectra

Interferometric radius measurements provide a direct probe of the fundamental parameters of M dwarfs, but is within reach for only a limited sample of nearby, bright stars. We use interferometrically-measured radii, bolometric luminosities, and effective temperatures to develop new empirical calibrations based on low-resolution, near-infrared spectra. We use H-band Mg and Al features to derive calibrations for effective temperature, radius and log luminosity; the standard deviations in the residuals of our best fits are, respectively, 73K, 0.027Rsun, and 0.049 dex (11% error on luminosity). These relationships are valid for mid K to mid M dwarf stars, roughly corresponding to temperatures between 3100 and 4800K. We apply our calibrations to M dwarfs targeted by the MEarth transiting planet survey and to the cool Kepler Objects of Interest (KOIs). We independently validate our calibrations by demonstrating a clear relationship between our inferred parameters and the absolute K magnitudes of the MEarth stars, and we identify objects with magnitudes too bright for their estimated luminosities as candidate multiple systems. We also use our inferred luminosities to address the applicability of near-infrared metallicity calibrations to mid and late M dwarfs. The temperatures we infer for the KOIs agree remarkably well with those from the literature; however, our stellar radii are systematically larger than those presented in previous works that derive radii from model isochrones. This results in a mean planet radius that is 15% larger than one would infer using the stellar properties from recent catalogs. Our results confirm those of previous in-depth studies of Kepler-42, Kepler-45, and Kepler-186.

preprint2015arXiv

Calibration of the MEarth Photometric System: Optical Magnitudes and Photometric Metallicity Estimates for 1802 Nearby M-dwarfs

The MEarth Project is a photometric survey systematically searching the smallest stars nearest to the Sun for transiting rocky planets. Since 2008, MEarth has taken approximately two million images of 1844 stars suspected to be mid-to-late M dwarfs. We have augmented this survey by taking nightly exposures of photometric standard stars and have utilized this data to photometrically calibrate the $MEarth$ system, identify photometric nights, and obtain an optical magnitude with $1.5\%$ precision for each M dwarf system. Each optical magnitude is an average over many years of data, and therefore should be largely immune to stellar variability and flaring. We combine this with trigonometric distance measurements, spectroscopic metallicity measurements, and 2MASS infrared magnitude measurements in order to derive a color-magnitude-metallicity relation across the mid-to-late M dwarf spectral sequence that can reproduce spectroscopic metallicity determinations to a precision of 0.1 dex. We release optical magnitudes and metallicity estimates for 1567 M dwarfs, many of which did not have an accurate determination of either prior to this work. For an additional 277 stars without a trigonometric parallax, we provide an estimate of the distance assuming solar neighborhood metallicity. We find that the median metallicity for a volume limited sample of stars within 20 parsecs of the Sun is [Fe/H] = $-0.03 \pm 0.008$, and that 29 / 565 of these stars have a metallicity of [Fe/H] = $-0.5$ or lower, similar to the low-metallicity distribution of nearby G-dwarfs. When combined with the results of ongoing and future planet surveys targeting these objects, the metallicity estimates presented here will be important in assessing the significance of any putative planet-metallicity correlation.

preprint2015arXiv

Dynamical Masses of Young M Dwarfs: Masses and Orbital Parameters of GJ 3305 AB, the Wide Binary Companion to the Imaged Exoplanet Host 51 Eri

We combine new high resolution imaging and spectroscopy from Keck/NIRC2, Discovery Channel Telescope/DSSI, and Keck/HIRES with published astrometry and radial velocities to measure individual masses and orbital elements of the GJ 3305 AB system, a young (~20 Myr) M+M binary (unresolved spectral type M0) member of the beta Pictoris moving group comoving with the imaged exoplanet host 51 Eri. We measure a total system mass of 1.11 \pm 0.04 M_sun, a period of 29.03 \pm 0.50$ yr, a semimajor axis of 9.78 \pm 0.14 AU, and an eccentricity of 0.19 \pm 0.02. The primary component has a dynamical mass of 0.67 \pm 0.05 M_sun and the secondary has a mass of 0.44 \pm 0.05 M_sun. The recently updated BHAC15 models are consistent with the masses of both stars to within 1.5 sigma. Given the observed masses the models predict an age of the GJ 3305 AB system of 37 \pm 9 Myr. Based on the the observed system architecture and our dynamical mass measurement, it is unlikely that the orbit of 51 Eri b has been significantly altered by the Kozai-Lidov mechanism.

preprint2015arXiv

HARPS-N observes the Sun as a star

Radial velocity perturbations induced by stellar surface inhomogeneities including spots, plages and granules currently limit the detection of Earth-twins using Doppler spectroscopy. Such stellar noise is poorly understood for stars other than the Sun because their surface is unresolved. In particular, the effects of stellar surface inhomogeneities on observed stellar radial velocities are extremely difficult to characterize, and thus developing optimal correction techniques to extract true stellar radial velocities is extremely challenging. In this paper, we present preliminary results of a solar telescope built to feed full-disk sunlight into the HARPS-N spectrograph, which is in turn calibrated with an astro-comb. This setup enables long-term observation of the Sun as a star with state-of-the-art sensitivity to radial velocity changes. Over seven days of observing in 2014, we show an average 50\cms radial velocity rms over a few hours of observation. After correcting observed radial velocities for spot and plage perturbations using full-disk photometry of the Sun, we lower by a factor of two the weekly radial velocity rms to 60\cms. The solar telescope is now entering routine operation, and will observe the Sun every clear day for several hours. We will use these radial velocities combined with data from solar satellites to improve our understanding of stellar noise and develop optimal correction methods. If successful, these new methods should enable the detection of Venus over the next two to three years, thus demonstrating the possibility of detecting Earth-twins around other solar-like stars using the radial velocity technique.

preprint2015arXiv

Low False-Positive Rate of Kepler Candidates Estimated From A Combination Of Spitzer And Follow-Up Observations

(Abridged) NASA's Kepler mission has provided several thousand transiting planet candidates, yet only a small subset have been confirmed as true planets. Therefore, the most fundamental question about these candidates is the fraction of bona fide planets. Estimating the rate of false positives of the overall Kepler sample is necessary to derive the planet occurrence rate. We present the results from two large observational campaigns that were conducted with the Spitzer telescope during the the Kepler mission. These observations are dedicated to estimating the false positive rate (FPR) amongst the Kepler candidates. We select a sub-sample of 51 candidates, spanning wide ranges in stellar, orbital and planetary parameter space, and we observe their transits with Spitzer at 4.5 microns. We use these observations to measures the candidate's transit depths and infrared magnitudes. A bandpass-dependent depth alerts us to the potential presence of a blending star that could be the source of the observed eclipse: a false-positive scenario. For most of the candidates (85%), the transit depths measured with Kepler are consistent with the depths measured with Spitzer as expected for planetary objects, while we find that the most discrepant measurements are due to the presence of unresolved stars that dilute the photometry. The Spitzer constraints on their own yield FPRs between 5-40%, depending on the KOIs. By considering the population of the Kepler field stars, and by combining follow-up observations (imaging) when available, we find that the overall FPR of our sample is low. The measured upper limit on the FPR of our sample is 8.8% at a confidence level of 3 sigma. This observational result, which uses the achromatic property of planetary transit signals that is not investigated by the Kepler observations, provides an independent indication that Kepler's false positive rate is low.

preprint2015arXiv

The Kepler-454 System: A Small, Not-rocky Inner Planet, a Jovian World, and a Distant Companion

Kepler-454 (KOI-273) is a relatively bright (V = 11.69 mag), Sun-like starthat hosts a transiting planet candidate in a 10.6 d orbit. From spectroscopy, we estimate the stellar temperature to be 5687 +/- 50 K, its metallicity to be [m/H] = 0.32 +/- 0.08, and the projected rotational velocity to be v sin i <2.4 km s-1. We combine these values with a study of the asteroseismic frequencies from short cadence Kepler data to estimate the stellar mass to be 1.028+0:04-0:03 M_Sun, the radius to be 1.066 +/- 0.012 R_Sun and the age to be 5.25+1:41-1:39 Gyr. We estimate the radius of the 10.6 d planet as 2.37 +/- 0.13 R_Earth. Using 63 radial velocity observations obtained with the HARPS-N spectrograph on the Telescopio Nazionale Galileo and 36 observations made with the HIRES spectrograph at Keck Observatory, we measure the mass of this planet to be 6.8 +/- 1.4M_Earth. We also detect two additional non-transiting companions, a planet with a minimum mass of 4.46 +/- 0.12 M_J in a nearly circular 524 d orbit and a massive companion with a period >10 years and mass >12.1M_J . The twelve exoplanets with radii <2.7 R_Earth and precise mass measurements appear to fall into two populations, with those <1.6 R_Earth following an Earth-like composition curve and larger planets requiring a significant fraction of volatiles. With a density of 2.76 +/- 0.73 g cm-3, Kepler-454b lies near the mass transition between these two populations and requires the presence of volatiles and/or H/He gas.

preprint2015arXiv

The Occurrence of Potentially Habitable Planets Orbiting M Dwarfs Estimated from the Full Kepler Dataset and an Empirical Measurement of the Detection Sensitivity

We present an improved estimate of the occurrence rate of small planets orbiting small stars by searching the full four-year Kepler data set for transiting planets using our own planet detection pipeline and conducting transit injection and recovery simulations to empirically measure the search completeness of our pipeline. We identified 156 planet candidates, including one object that was not previously identified as a Kepler Object of Interest. We inspected all publicly available follow-up images, observing notes, and centroid analyses, and corrected for the likelihood of false positives. We evaluated the sensitivity of our detection pipeline on a star-by-star basis by injecting 2000 transit signals into the light curve of each target star. For periods shorter than 50 days, we find 0.56 (+0.06/-0.05) Earth-size planets (1-1.5 Earth radii) and 0.46 (+0.07/-0.05) super-Earths (1.5-2 Earth radii) per M dwarf. In total, we estimate a cumulative planet occurrence rate of $2.5\pm0.2$ planets per M dwarf with radii 1-4 Earth radii and periods shorter than 200 days. Within a conservatively defined habitable zone based on the moist greenhouse inner limit and maximum greenhouse outer limit, we estimate an occurrence rate of 0.16 (+0.17/-0.07) Earth-size planets and 0.12 (+0.10/-0.05) super-Earths per M dwarf habitable zone. Adopting the broader insolation boundaries of the recent Venus and early Mars limits yields a higher estimate of 0.24 (+0.18/-0.08) Earth-size planets and 0.21 (+0.11/-0.06) super-Earths per M dwarf habitable zone. This suggests that the nearest potentially habitable non-transiting and transiting Earth-size planets are $2.6\pm0.4$ pc and 10.6 (+1.6/-1.8) pc away, respectively. If we include super-Earths, these distances diminish to $2.1\pm0.2$ pc and 8.6 (+0.7/-0.8) pc.

preprint2014arXiv

A Precise Water Abundance Measurement for the Hot Jupiter WASP-43b

The water abundance in a planetary atmosphere provides a key constraint on the planet's primordial origins because water ice is expected to play an important role in the core accretion model of planet formation. However, the water content of the Solar System giant planets is not well known because water is sequestered in clouds deep in their atmospheres. By contrast, short-period exoplanets have such high temperatures that their atmospheres have water in the gas phase, making it possible to measure the water abundance for these objects. We present a precise determination of the water abundance in the atmosphere of the 2 $M_\mathrm{Jup}$ short-period exoplanet WASP-43b based on thermal emission and transmission spectroscopy measurements obtained with the Hubble Space Telescope. We find the water content is consistent with the value expected in a solar composition gas at planetary temperatures (0.4-3.5x solar at 1 $σ$ confidence). The metallicity of WASP-43b's atmosphere suggested by this result extends the trend observed in the Solar System of lower metal enrichment for higher planet masses.

preprint2014arXiv

Atmospheric Characterization of 5 Hot Jupiters with Wide Field Camera 3 on the Hubble Space Telescope

We probe the structure and composition of the atmospheres of 5 hot Jupiter exoplanets using the Hubble Space Telescope Wide Field Camera 3 (WFC3) instrument. We use the G141 grism (1.1-1.7 $μ$m) to study TrES-2b, TrES-4b, and CoRoT-1b in transit, TrES-3b in secondary eclipse, and WASP-4b in both. This wavelength region includes a predicted absorption feature from water at 1.4 $μ$m, which we expect to be nondegenerate with the other molecules that are likely to be abundant for hydrocarbon-poor (e.g. solar composition) hot Jupiter atmospheres. We divide our wavelength regions into 10 bins. For each bin we produce a spectrophotometric light curve spanning the time of transit and/or eclipse. We correct these light curves for instrumental systematics without reference to an instrument model. For our transmission spectra, our mean $1-σ$ precision per bin corresponds to variations of 2.1, 2.8, and 3.0 atmospheric scale heights for TrES-2b, TrES-4b, and CoRoT-1b, respectively. We find featureless spectra for these three planets. We are unable to extract a robust transmission spectrum for WASP-4b. For our dayside emission spectra, our mean $1-σ$ precision per bin corresponds to a planet-to-star flux ratio of $1.5\times10^{-4}$ and $2.1\times10^{-4}$ for WASP-4b and TrES-3b, respectively. We combine these estimates with previous broadband measurements and conclude that for both planets isothermal atmospheres are disfavored. We find no signs of features due to water. We confirm that WFC3 is suitable for studies of transiting exoplanets, but in staring mode multi-visit campaigns are necessary to place strong constraints on water abundance.

preprint2014arXiv

Characterizing K2 Planet Discoveries: A super-Earth transiting the bright K-dwarf HIP 116454

We report the first planet discovery from the two-wheeled Kepler (K2) mission: HIP 116454 b. The host star HIP 116454 is a bright (V = 10.1, K = 8.0) K1-dwarf with high proper motion, and a parallax-based distance of 55.2 +/- 5.4 pc. Based on high-resolution optical spectroscopy, we find that the host star is metal-poor with [Fe/H] = -.16 +/- .18, and has a radius R = 0.716 +/- .0024 R_sun and mass M = .775 +/- .027 Msun. The star was observed by the Kepler spacecraft during its Two-Wheeled Concept Engineering Test in February 2014. During the 9 days of observations, K2 observed a single transit event. Using a new K2 photometric analysis technique we are able to correct small telescope drifts and recover the observed transit at high confidence, corresponding to a planetary radius of Rp = 2.53 +/- 0.18 Rearth. Radial velocity observations with the HARPS-N spectrograph reveal a 11.82 +/- 1.33 Mearth planet in a 9.1 day orbit, consistent with the transit depth, duration, and ephemeris. Follow-up photometric measurements from the MOST satellite confirm the transit observed in the K2 photometry and provide a refined ephemeris, making HIP 116454 b amenable for future follow-up observations of this latest addition to the growing population of transiting super-Earths around nearby, bright stars.

preprint2014arXiv

Kepler-93b: A Terrestrial World Measured to within 120 km, and a Test Case for a New Spitzer Observing Mode

We present the characterization of the Kepler-93 exoplanetary system, based on three years of photometry gathered by the Kepler spacecraft. The duration and cadence of the Kepler observations, in tandem with the brightness of the star, enable unusually precise constraints on both the planet and its host. We conduct an asteroseismic analysis of the Kepler photometry and conclude that the star has an average density of 1.652+/-0.006 g/cm^3. Its mass of 0.911+/-0.033 M_Sun renders it one of the lowest-mass subjects of asteroseismic study. An analysis of the transit signature produced by the planet Kepler-93b, which appears with a period of 4.72673978+/-9.7x10^-7 days, returns a consistent but less precise measurement of the stellar density, 1.72+0.02-0.28 g/cm^3. The agreement of these two values lends credence to the planetary interpretation of the transit signal. The achromatic transit depth, as compared between Kepler and the Spitzer Space Telescope, supports the same conclusion. We observed seven transits of Kepler-93b with Spitzer, three of which we conducted in a new observing mode. The pointing strategy we employed to gather this subset of observations halved our uncertainty on the transit radius ratio R_p/R_star. We find, after folding together the stellar radius measurement of 0.919+/-0.011 R_Sun with the transit depth, a best-fit value for the planetary radius of 1.481+/-0.019 R_Earth. The uncertainty of 120 km on our measurement of the planet's size currently renders it one of the most precisely measured planetary radii outside of the Solar System. Together with the radius, the planetary mass of 3.8+/-1.5 M_Earth corresponds to a rocky density of 6.3+/-2.6 g/cm^3. After applying a prior on the plausible maximum densities of similarly-sized worlds between 1--1.5 R_Earth, we find that Kepler-93b possesses an average density within this group.

preprint2014arXiv

Masses, Radii, and Orbits of Small Kepler Planets: The Transition from Gaseous to Rocky Planets

We report on the masses, sizes, and orbits of the planets orbiting 22 Kepler stars. There are 49 planet candidates around these stars, including 42 detected through transits and 7 revealed by precise Doppler measurements of the host stars. Based on an analysis of the Kepler brightness measurements, along with high-resolution imaging and spectroscopy, Doppler spectroscopy, and (for 11 stars) asteroseismology, we establish low false-positive probabilities for all of the transiting planets (41 of 42 have a false-positive probability under 1%), and we constrain their sizes and masses. Most of the transiting planets are smaller than 3X the size of Earth. For 16 planets, the Doppler signal was securely detected, providing a direct measurement of the planet's mass. For the other 26 planets we provide either marginal mass measurements or upper limits to their masses and densities; in many cases we can rule out a rocky composition. We identify 6 planets with densities above 5 g/cc, suggesting a mostly rocky interior for them. Indeed, the only planets that are compatible with a purely rocky composition are smaller than ~2 R_earth. Larger planets evidently contain a larger fraction of low-density material (H, He, and H2O).

preprint2014arXiv

The 4.5 $μ$m full-orbit phase curve of the hot Jupiter HD 209458b

The hot Jupiter HD 209458b is particularly amenable to detailed study as it is among the brightest transiting exoplanet systems currently known (V-mag = 7.65; K-mag = 6.308) and has a large planet-to-star contrast ratio. HD 209458b is predicted to be in synchronous rotation about its host star with a hot spot that is shifted eastward of the substellar point by superrotating equatorial winds. Here we present the first full-orbit observations of HD 209458b, in which its 4.5 $μ$m emission was recorded with $Spitzer$/IRAC. Our study revises the previous 4.5 $μ$m measurement of HD 209458b's secondary eclipse emission downward by $\sim$35% to $0.1391%^{+0.0072%}_{-0.0069%}$, changing our interpretation of the properties of its dayside atmosphere. We find that the hot spot on the planet's dayside is shifted eastward of the substellar point by $40.9^{\circ}\pm{6.0^{\circ}}$, in agreement with circulation models predicting equatorial superrotation. HD 209458b's dayside (T$_{bright}$ = 1499 $\pm$ 15 K) and nightside (T$_{bright}$ = 972 $\pm$ 44 K) emission indicates a day-to-night brightness temperature contrast smaller than that observed for more highly irradiated exoplanets, suggesting that the day-to-night temperature contrast may be partially a function of the incident stellar radiation. The observed phase curve shape deviates modestly from global circulation model predictions potentially due to disequilibrium chemistry or deficiencies in the current hot CH$_{4}$ line lists used in these models. Observations of the phase curve at additional wavelengths are needed in order to determine the possible presence and spatial extent of a dayside temperature inversion, as well as to improve our overall understanding of this planet's atmospheric circulation.

preprint2014arXiv

The Kepler-10 planetary system revisited by HARPS-N: A hot rocky world and a solid Neptune-mass planet

Kepler-10b was the first rocky planet detected by the Kepler satellite and con- firmed with radial velocity follow-up observations from Keck-HIRES. The mass of the planet was measured with a precision of around 30%, which was insufficient to constrain models of its internal structure and composition in detail. In addition to Kepler-10b, a second planet transiting the same star with a period of 45 days was sta- tistically validated, but the radial velocities were only good enough to set an upper limit of 20 Mearth for the mass of Kepler-10c. To improve the precision on the mass for planet b, the HARPS-N Collaboration decided to observe Kepler-10 intensively with the HARPS-N spectrograph on the Telescopio Nazionale Galileo on La Palma. In to- tal, 148 high-quality radial-velocity measurements were obtained over two observing seasons. These new data allow us to improve the precision of the mass determina- tion for Kepler-10b to 15%. With a mass of 3.33 +/- 0.49 Mearth and an updated radius of 1.47 +0.03 -0.02 Rearth, Kepler-10b has a density of 5.8 +/- 0.8 g cm-3, very close to the value -0.02 predicted by models with the same internal structure and composition as the Earth. We were also able to determine a mass for the 45-day period planet Kepler-10c, with an even better precision of 11%. With a mass of 17.2 +/- 1.9 Mearth and radius of 2.35 +0.09 -0.04 Rearth, -0.04 Kepler-10c has a density of 7.1 +/- 1.0 g cm-3. Kepler-10c appears to be the first strong evidence of a class of more massive solid planets with longer orbital periods.

preprint2014arXiv

The Mass of Kepler-93b and The Composition of Terrestrial Planets

Kepler-93b is a 1.478 +/- 0.019 Earth radius planet with a 4.7 day period around a bright (V=10.2), astroseismically-characterized host star with a mass of 0.911+/-0.033 solar masses and a radius of 0.919+/-0.011 solar radii. Based on 86 radial velocity observations obtained with the HARPS-N spectrograph on the Telescopio Nazionale Galileo and 32 archival Keck/HIRES observations, we present a precise mass estimate of 4.02+/-0.68 Earth masses. The corresponding high density of 6.88+/-1.18 g/cc is consistent with a rocky composition of primarily iron and magnesium silicate. We compare Kepler-93b to other dense planets with well-constrained parameters and find that between 1-6 Earth masses, all dense planets including the Earth and Venus are well-described by the same fixed ratio of iron to magnesium silicate. There are as of yet no examples of such planets with masses > 6 Earth masses: All known planets in this mass regime have lower densities requiring significant fractions of volatiles or H/He gas. We also constrain the mass and period of the outer companion in the Kepler-93 system from the long-term radial velocity trend and archival adaptive optics images. As the sample of dense planets with well-constrained masses and radii continues to grow, we will be able to test whether the fixed compositional model found for the seven dense planets considered in this paper extends to the full population of 1-6 Earth mass planets.

preprint2014arXiv

The MEarth-North and MEarth-South transit surveys: searching for habitable super-Earth exoplanets around nearby M-dwarfs

Detection and characterization of potentially habitable Earth-size extrasolar planets is one of the major goals of contemporary astronomy. By applying the transit method to very low-mass M-dwarfs, it is possible to find these planets from the ground with present-day instrumentation and observational techniques. The MEarth project is one such survey with stations in both hemispheres: MEarth-North at the Fred Lawrence Whipple Observatory, Mount Hopkins, Arizona, and MEarth-South at Cerro Tololo Inter-American Observatory, Chile. We present an update on recent results of this survey, for planet occurrence rates, and interesting stellar astrophysics, for which our sample of 3000 nearby mid-to-late M-dwarfs has been very fruitful. All light curves gathered during the survey are made publicly available after one year, and we describe how to access and use these data.

preprint2014arXiv

The Transiting Exoplanet Survey Satellite

The Transiting Exoplanet Survey Satellite (TESS) will search for planets transiting bright and nearby stars. TESS has been selected by NASA for launch in 2017 as an Astrophysics Explorer mission. The spacecraft will be placed into a highly elliptical 13.7-day orbit around the Earth. During its two-year mission, TESS will employ four wide-field optical CCD cameras to monitor at least 200,000 main-sequence dwarf stars with I = 4-13 for temporary drops in brightness caused by planetary transits. Each star will be observed for an interval ranging from one month to one year, depending mainly on the star's ecliptic latitude. The longest observing intervals will be for stars near the ecliptic poles, which are the optimal locations for follow-up observations with the James Webb Space Telescope. Brightness measurements of preselected target stars will be recorded every 2 min, and full frame images will be recorded every 30 min. TESS stars will be 10-100 times brighter than those surveyed by the pioneering Kepler mission. This will make TESS planets easier to characterize with follow-up observations. TESS is expected to find more than a thousand planets smaller than Neptune, including dozens that are comparable in size to the Earth. Public data releases will occur every four months, inviting immediate community-wide efforts to study the new planets. The TESS legacy will be a catalog of the nearest and brightest stars hosting transiting planets, which will endure as highly favorable targets for detailed investigations.

preprint2014arXiv

Thermal structure of an exoplanet atmosphere from phase-resolved emission spectroscopy

Exoplanets that orbit close to their host stars are much more highly irradiated than their Solar System counterparts. Understanding the thermal structures and appearances of these planets requires investigating how their atmospheres respond to such extreme stellar forcing. We present spectroscopic thermal emission measurements as a function of orbital phase ("phase-curve observations") for the highly-irradiated exoplanet WASP-43b spanning three full planet rotations using the Hubble Space Telescope. With these data, we construct a map of the planet's atmospheric thermal structure, from which we find large day-night temperature variations at all measured altitudes and a monotonically decreasing temperature with pressure at all longitudes. We also derive a Bond albedo of 0.18 +0.07,-0.12 and an altitude dependence in the hot-spot offset relative to the substellar point.

preprint2013arXiv

A sub-Mercury-sized exoplanet

Since the discovery of the first exoplanet we have known that other planetary systems can look quite unlike our own. However, until recently we have only been able to probe the upper range of the planet size distribution. The high precision of the Kepler space telescope has allowed us to detect planets that are the size of Earth and somewhat smaller, but no previous planets have been found that are smaller than those we see in our own Solar System. Here we report the discovery of a planet significantly smaller than Mercury. This tiny planet is the innermost of three planets that orbit the Sun-like host star, which we have designated Kepler-37. Owing to its extremely small size, similar to that of Earth's Moon, and highly irradiated surface, Kepler-37b is probably a rocky planet with no atmosphere or water, similar to Mercury.

preprint2013arXiv

An Earth-sized planet with an Earth-like density

Kepler-78 (KIC 8435766) was identified by Sanchis-Ojeda et al. (2013) as harbouring a transiting planet of 1.16 times the size of the Earth and an orbital period of only 8.5 hours. While the exquisite Kepler photometry was able to determine its radius and period, the mass of the planet (and thus its mean density) remained unconstrained in the absence of precise radial-velocity measurements. Here we present an accurate mass measurement of Kepler-78b using the HARPS-N spectrograph, recently installed on the Telescopio Nazionale Galileo (INAF) at the Roque de los Muchachos Observatory, La Palma, Spain. These new data yield a mass of 1.86 Earth masses. The resulting mean density of the planet is 5.57 grams per cubic centimetre, which is similar to that of the Earth and implies a composition of iron and rock. Kepler-78b, which orbits a Sun-like star called Kepler 78 located in the Cygnus constellation at a distance of about 400 light years from us, is now the smallest exoplanet for which both the mass and radius are known accurately.

preprint2013arXiv

Constraints on Planet Occurrence around Nearby Mid-to-Late M Dwarfs from the MEarth Project

The MEarth Project is a ground-based photometric survey to find planets transiting the closest and smallest main-sequence stars. In its first four years, MEarth discovered one transiting exoplanet, the 2.7 Earth radius planet GJ1214b. Here, we answer an outstanding question: in light of the bounty of small planets transiting small stars uncovered by the Kepler mission, should MEarth have found more than just one planet so far? We estimate MEarth's ensemble sensitivity to exoplanets by performing end-to-end simulations of 1.25 million observations of 988 nearby mid-to-late M dwarfs, gathered by MEarth between October 2008 and June 2012. For 2-4 Earth radius planets, we compare this sensitivity to results from Kepler and find that MEarth should have found planets at a rate of 0.05 - 0.36 planets/year in its first four years. As part of this analysis, we provide new analytic fits to the Kepler early M dwarf planet occurrence distribution. When extrapolating between Kepler's early M dwarfs and MEarth's mid-to-late M dwarfs, we find that assuming the planet occurrence distribution stays fixed with respect to planetary equilibrium temperature provides a good match to our detection of a planet with GJ1214b's observed properties. For larger planets, we find that the warm (600-700K), Neptune-sized (4 Earth radius) exoplanets that transit early M dwarfs like Gl436 and GJ3470 occur at a rate of <0.15/star (at 95% confidence) around MEarth's later M dwarf targets. We describe a strategy with which MEarth can increase its expected planet yield by 2.5X without new telescopes, by shifting its sensitivity toward the smaller and cooler exoplanets that Kepler has demonstrated to be abundant.

preprint2013arXiv

Exoplanet Characterization by Proxy: a Transiting 2.15 R_Earth Planet Near the Habitable Zone of the Late K dwarf Kepler-61

We present the validation and characterization of Kepler-61b: a 2.15 R_Earth planet orbiting near the inner edge of the habitable zone of a low-mass star. Our characterization of the host star Kepler-61 is based upon a comparison with the set of spectroscopically similar stars with directly-measured radii and temperatures. We apply a stellar prior drawn from the weighted mean of these properties, in tandem with the Kepler photometry, to infer a planetary radius for Kepler-61b of 2.15+/-0.13 R_Earth and an equilibrium temperature of 273+/-13 K (given its period of 59.87756+/-0.00020 days and assuming a planetary albedo of 0.3). The technique of leveraging the physical properties of nearby "proxy" stars allows for an independent check on stellar characterization via the traditional measurements with stellar spectra and evolutionary models. In this case, such a check had implications for the putative habitability of Kepler-61b: the planet is 10% warmer and larger than inferred from K-band spectral characterization. From the Kepler photometry, we estimate a stellar rotation period of 36 days, which implies a stellar age of >1 Gyr. We summarize the evidence for the planetary nature of the Kepler-61 transit signal, which we conclude is 30,000 times more likely to be due to a planet than a blend scenario. Finally, we discuss possible compositions for Kepler-61b with a comparison to theoretical models as well as to known exoplanets with similar radii and dynamically measured masses.

preprint2013arXiv

Infrared Transmission Spectroscopy of the Exoplanets HD209458b and XO-1b Using the Wide Field Camera-3 on the Hubble Space Telescope

Exoplanetary transmission spectroscopy in the near-infrared using Hubble/NICMOS is currently ambiguous because different observational groups claim different results from the same data, depending on their analysis methodologies. Spatial scanning with Hubble/WFC3 provides an opportunity to resolve this ambiguity. We here report WFC3 spectroscopy of the giant planets HD209458b and XO-1b in transit, using spatial scanning mode for maximum photon-collecting efficiency. We introduce an analysis technique that derives the exoplanetary transmission spectrum without the necessity of explicitly decorrelating instrumental effects, and achieves nearly photon-limited precision even at the high flux levels collected in spatial scan mode. Our errors are within 6-percent (XO-1) and 26-percent (HD209458b) of the photon-limit at a spectral resolving power of 70, and are better than 0.01-percent per spectral channel. Both planets exhibit water absorption of approximately 200 ppm at the water peak near 1.38 microns. Our result for XO-1b contradicts the much larger absorption derived from NICMOS spectroscopy. The weak water absorption we measure for HD209458b is reminiscent of the weakness of sodium absorption in the first transmission spectroscopy of an exoplanet atmosphere by Charbonneau et al. (2002). Model atmospheres having uniformly-distributed extra opacity of 0.012 cm^2 per gram account approximately for both our water measurement and the sodium absorption in this planet. Our results for HD209458b support the picture advocated by Pont et al. (2013) in which weak molecular absorptions are superposed on a transmission spectrum that is dominated by continuous opacity due to haze and/or dust. However, the extra opacity needed for HD209458b is grayer than for HD189733b, with a weaker Rayleigh component.

preprint2013arXiv

Kepler-68: Three Planets, One With a Density Between That of Earth and Ice Giants

NASA's Kepler Mission has revealed two transiting planets orbiting Kepler-68. Follow-up Doppler measurements have established the mass of the innermost planet and revealed a third jovian-mass planet orbiting beyond the two transiting planets. Kepler-68b, in a 5.4 day orbit has mass 8.3 +/- 2.3 Earth, radius 2.31 +/- 0.07 Earth radii, and a density of 3.32 +/- 0.92 (cgs), giving Kepler-68b a density intermediate between that of the ice giants and Earth. Kepler-68c is Earth-sized with a radius of 0.953 Earth and transits on a 9.6 day orbit; validation of Kepler-68c posed unique challenges. Kepler-68d has an orbital period of 580 +/- 15 days and minimum mass of Msin(i) = 0.947 Jupiter. Power spectra of the Kepler photometry at 1-minute cadence exhibit a rich and strong set of asteroseismic pulsation modes enabling detailed analysis of the stellar interior. Spectroscopy of the star coupled with asteroseismic modeling of the multiple pulsation modes yield precise measurements of stellar properties, notably Teff = 5793 +/- 74 K, M = 1.079 +/- 0.051 Msun, R = 1.243 +/- 0.019 Rsun, and density 0.7903 +/- 0.0054 (cgs), all measured with fractional uncertainties of only a few percent. Models of Kepler-68b suggest it is likely composed of rock and water, or has a H and He envelope to yield its density of about 3 (cgs).

preprint2013arXiv

Near-infrared Metallicities, Radial Velocities and Spectral Types for 447 Nearby M Dwarfs

We present metallicities, radial velocities and near-infrared spectral types for 447 M dwarfs determined from moderate resolution (R~2000) near-infrared (NIR) spectra obtained with IRTF/SpeX. These M dwarfs are targets of the MEarth Survey, a transiting planet survey searching for super Earths around mid-to-late M dwarfs within 33pc. We present NIR spectral types and new IRTF spectral templates in the Z, J, H and K-bands, created using M dwarfs with near-solar metallicities. We developed two spectroscopic distance calibrations that use NIR spectral type or an index based on the curvature of the K-band continuum. Our distance calibration has a scatter of 14%. We searched 27 NIR spectral lines and 10 spectral indices for metallicity sensitive features, taking into account correlated noise in our estimates of the errors on these parameters. We calibrated our relation using 36 M dwarfs in common proper pairs with an F, G or K-type star of known metallicity. We validated the physical association of these pairs using proper motions, radial velocities and spectroscopic distance estimates. Our resulting metallicity calibration uses the sodium doublet at 2.2um as the sole indicator for metallicity. It has an accuracy of 0.12dex inferred from the scatter between the metallicities of the primaries and the estimated metallicities of the secondaries. Our relation is valid for NIR spectral types from M1V to M5V and for -1.0<[Fe/H]<+0.35dex. We present a new metallicity relation using J-H and J-K colors that directly relates two observables: the distance from the M dwarf main sequence and equivalent width of the sodium line at 2.2um. We measured radial velocities by modeling telluric features to determine the absolute wavelength calibration of our spectra, and used M dwarf binaries, observations at different epochs, and comparison to precisely measured radial velocities to demonstrate 4.4km/s accuracy.

preprint2013arXiv

Orbital Phase Variations of the Eccentric Giant Planet HAT-P-2b

We present the first secondary eclipse and phase curve observations for the highly eccentric hot Jupiter HAT-P-2b in the 3.6, 4.5, 5.8, and 8.0 μm bands of the Spitzer Space Telescope. The 3.6 and 4.5 μm data sets span an entire orbital period of HAT-P-2b, making them the longest continuous phase curve observations obtained to date and the first full-orbit observations of a planet with an eccentricity exceeding 0.2. We present an improved non-parametric method for removing the intrapixel sensitivity variations in Spitzer data at 3.6 and 4.5 μm that robustly maps position-dependent flux variations. We find that the peak in planetary flux occurs at 4.39+/-0.28, 5.84+/-0.39, and 4.68+/-0.37 hours after periapse passage with corresponding maxima in the planet/star flux ratio of 0.1138%+/-0.0089%, 0.1162%+/-0.0080%, and 0.1888%+/-0.0072% in the 3.6, 4.5, and 8.0 μm bands respectively. We compare our measured secondary eclipse depths to the predictions from a one-dimensional radiative transfer model, which suggests the possible presence of a transient day side inversion in HAT-P-2b's atmosphere near periapse. We also derive improved estimates for the system parameters, including its mass, radius, and orbital ephemeris. Our simultaneous fit to the transit, secondary eclipse, and radial velocity data allows us to determine the eccentricity and argument of periapse of HAT-P-2b's orbit with a greater precision than has been achieved for any other eccentric extrasolar planet. We also find evidence for a long-term linear trend in the radial velocity data. This trend suggests the presence of another substellar companion in the HAT-P-2 system, which could have caused HAT-P-2b to migrate inward to its present-day orbit via the Kozai mechanism.

preprint2013arXiv

Secondary Eclipse Photometry of the Exoplanet WASP-5b with Warm Spitzer

We present photometry of the extrasolar planet WASP-5b in the 3.6 and 4.5 micron bands taken with the Spitzer Space Telescope's Infrared Array Camera as part of the extended warm mission. By examining the depth of the planet's secondary eclipse at these two wavelengths, we can place joint constraints on the planet's atmospheric pressure-temperature profile and chemistry. We measure secondary eclipse depths of 0.197% +/- 0.028% and 0.227% +/- 0.025% in the 3.6 micron and 4.5 micron bands, respectively. Our observations are best matched by models showing a hot dayside and, depending on our choice of model, a weak thermal inversion or no inversion at all. We measure a mean offset from the predicted center of eclipse of 0.078 +/- 0.032 hours, translating to ecos(omega) = 0.0031 +/- 0.0013 and consistent with a circular orbit. We see no evidence for any eclipse timing variations comparable to those reported in a previous transit study.

preprint2013arXiv

The false positive rate of Kepler and the occurrence of planets

The Kepler Mission is uniquely suited to study the frequencies of extrasolar planets. This goal requires knowledge of the incidence of false positives such as eclipsing binaries in the background of the targets, or physically bound to them, which can mimic the photometric signal of a transiting planet. We perform numerical simulations to predict the occurrence of astrophysical false positives detectable by the Mission, and quantify the fraction of them that would pass the Kepler candidate vetting procedure. By comparing their distribution with that of the Kepler Objects of Interest detected during the first six quarters of operation of the spacecraft, we infer the false positive rate of Kepler and study its dependence on spectral type, candidate planet size, and orbital period. We find that the global false positive rate of Kepler is 9.4 %, peaking for giant planets (6-22 Earth radii) at 17.7 %, reaching a low of 6.7 % for small Neptunes (2-4 Earth radii), and increasing again for Earth-size planets (0.8-1.25 Earth radii) to 12.3 %. Most importantly, we also quantify and characterize the distribution and rate of occurrence of planets down to Earth size with no prior assumptions on their frequency, by subtracting from the population of actual Kepler candidates our simulated population of astrophysical false positives. We find that 16.5 +/- 3.6 % of main-sequence FGK stars have at least one planet between 0.8 and 1.25 Earth radii with orbital periods up to 85 days. There is no significant dependence of the rates of planet occurrence between 0.8 and 4 Earth radii with spectral type. In the process, we derive also a prescription for the signal recovery rate of Kepler that enables a good match to both the KOI size and orbital period distribution, as well as their signal-to-noise distribution.

preprint2013arXiv

The Occurrence Rate of Small Planets around Small Stars

We use the optical and near-infrared photometry from the Kepler Input Catalog to provide improved estimates of the stellar characteristics of the smallest stars in the Kepler target list. We find 3897 dwarfs with temperatures below 4000K, including 64 planet candidate host stars orbited by 95 transiting planet candidates. We refit the transit events in the Kepler light curves for these planet candidates and combine the revised planet/star radius ratios with our improved stellar radii to revise the radii of the planet candidates orbiting the cool target stars. We then compare the number of observed planet candidates to the number of stars around which such planets could have been detected in order to estimate the planet occurrence rate around cool stars. We find that the occurrence rate of 0.5-4 Earth radius planets with orbital periods shorter than 50 days is 0.90 (+0.04/-0.03) planets per star. The occurrence rate of Earth-size (0.5-1.4 Earth radius) planets is constant across the temperature range of our sample at 0.51 (+0.06/-0.05) Earth-size planets per star, but the occurrence of 1.4-4 Earth radius planets decreases significantly at cooler temperatures. Our sample includes 2 Earth-size planet candidates in the habitable zone, allowing us to estimate that the mean number of Earth-size planets in the habitable zone is 0.15 (+0.13/-0.06) planets per cool star. Our 95% confidence lower limit on the occurrence rate of Earth-size planets in the habitable zones of cool stars is 0.04 planets per star. With 95% confidence, the nearest transiting Earth-size planet in the habitable zone of a cool star is within 21 pc. Moreover, the nearest non-transiting planet in the habitable zone is within 5 pc with 95% confidence.

preprint2013arXiv

Trigonometric Parallaxes for 1,507 Nearby Mid-to-Late M-dwarfs

The MEarth survey is a search for small rocky planets around the smallest, nearest stars to the Sun as identified by high proper motion with red colors. We augmented our planetary search time series with lower cadence astrometric imaging and obtained two million images of approximately 1800 stars suspected to be mid-to-late M dwarfs. We fit an astrometric model to MEarth's images for 1507 stars and obtained trigonometric distance measurements to each star with an average precision of 5 milliarcseconds. Our measurements, combined with the 2MASS photometry, allowed us to obtain an absolute K_s magnitude for each star. In turn, this allows us to better estimate the stellar parameters than those obtained with photometric estimates alone and to better prioritize the targets chosen to monitor at high cadence for planetary transits. The MEarth sample is mostly complete out to a distance of 25 parsecs for stars of type M5.5V and earlier, and mostly complete for later type stars out to 20 parsecs. We find eight stars that are within ten parsecs of the Sun for which there did not exist a published trigonometric parallax distance estimate. We release with this work a catalog of the trigonometric parallax measurements for 1,507 mid-to-late M-dwarfs, as well as new estimates of their masses and radii.

preprint2013arXiv

Warm Spitzer Photometry of Three Hot Jupiters: HAT-P-3b, HAT-P-4b and HAT-P-12b

We present Warm Spitzer/IRAC secondary eclipse time series photometry of three short-period transiting exoplanets, HAT-P-3b, HAT-P-4b and HAT-P-12b, in both the available 3.6 and 4.5 micron bands. HAT-P-3b and HAT-P-4b are Jupiter-mass, objects orbiting an early K and an early G dwarf stars, respectively. For HAT-P-3b we find eclipse depths of 0.112%+0.015%-0.030% (3.6 micron) and 0.094%+0.016%-0.009% (4.5 micron). The HAT-P-4b values are 0.142%+0.014%-0.016% (3.6 micron) and 0.122%+0.012%-0.014% (4.5micron). The two planets' photometry is consistent with inefficient heat redistribution from their day to night sides (and low albedos), but it is inconclusive about possible temperature inversions in their atmospheres. HAT-P-12b is a Saturn-mass planet and is one of the coolest planets ever observed during secondary eclipse, along with hot Neptune GJ 436b and hot Saturn WASP-29b. We are able to place 3$σ$ upper limits on the secondary eclipse depth of HAT-P-12b in both wavelengths: < 0.042% (3.6 micron) and <0.085% (4.5 micron). We discuss these results in the context of the {\it Spitzer} secondary eclipse measurements of GJ 436b and WASP-29b. It is possible that we do not detect the eclipses of HAT-P-12b due to high eccentricity, but find that weak planetary emission in these wavelengths is a more likely explanation. We place 3 sigma upper limits on the |e cos(omega)| quantity (where e is eccentricity and omega is the argument of periapsis) for HAT-P-3b (<0.0081) and HAT-P-4b (<0.0042), based on the secondary eclipse timings.

preprint2012arXiv

3.6 and 4.5 Micron Phase Curves and Evidence for Non-Equilibrium Chemistry in the Atmosphere of Extrasolar Planet HD 189733b

We present new, full-orbit observations of the infrared phase variations of the canonical hot Jupiter HD 189733b obtained in the 3.6 and 4.5 micron bands using the Spitzer Space Telescope. When combined with previous phase curve observations at 8.0 and 24 micron, these data allow us to characterize the exoplanet's emission spectrum as a function of planetary longitude. We utilize improved methods for removing the effects of intrapixel sensitivity variations and accounting for the presence of time-correlated noise in our data. We measure a phase curve amplitude of 0.1242% +/- 0.0061% in the 3.6 micron band and 0.0982% +/- 0.0089% in the 4.5 micron band. We find that the times of minimum and maximum flux occur several hours earlier than predicted for an atmosphere in radiative equilibrium, consistent with the eastward advection of gas by an equatorial super-rotating jet. The locations of the flux minima in our new data differ from our previous observations at 8 micron, and we present new evidence indicating that the flux minimum observed in the 8 micron is likely caused by an over-shooting effect in the 8 micron array. We obtain improved estimates for HD 189733b's dayside planet-star flux ratio of 0.1466% +/- 0.0040% at 3.6 micron and 0.1787% +/- 0.0038% at 4.5 micron; these are the most accurate secondary eclipse depths obtained to date for an extrasolar planet. We compare our new dayside and nightside spectra for HD 189733b to the predictions of models from Burrows et al. (2008) and Showman et al. (2009). We find that HD 189733b's 4.5 micron nightside flux is 3.3 sigma smaller than predicted by the Showman et al. models, which assume that the chemistry is in local thermal equilibrium. We conclude that this discrepancy is best-explained by vertical mixing, which should lead to an excess of CO and correspondingly enhanced 4.5 micron absorption in this region. [abridged]

preprint2012arXiv

Kepler-20: A Sun-like Star with Three Sub-Neptune Exoplanets and Two Earth-size Candidates

We present the discovery of the Kepler-20 planetary system, which we initially identified through the detection of five distinct periodic transit signals in the Kepler light curve of the host star 2MASSJ19104752+4220194. We find a stellar effective temperature Teff=5455+-100K, a metallicity of [Fe/H]=0.01+-0.04, and a surface gravity of log(g)=4.4+-0.1. Combined with an estimate of the stellar density from the transit light curves we deduce a stellar mass of Mstar=0.912+-0.034 Msun and a stellar radius of Rstar=0.944^{+0.060}_{-0.095} Rsun. For three of the transit signals, our results strongly disfavor the possibility that these result from astrophysical false positives. We conclude that the planetary scenario is more likely than that of an astrophysical false positive by a factor of 2e5 (Kepler-20b), 1e5 (Kepler-20c), and 1.1e3 (Kepler-20d), sufficient to validate these objects as planetary companions. For Kepler-20c and Kepler-20d, the blend scenario is independently disfavored by the achromaticity of the transit: From Spitzer data gathered at 4.5um, we infer a ratio of the planetary to stellar radii of 0.075+-0.015 (Kepler-20c) and 0.065+-0.011 (Kepler-20d), consistent with each of the depths measured in the Kepler optical bandpass. We determine the orbital periods and physical radii of the three confirmed planets to be 3.70d and 1.91^{+0.12}_{-0.21} Rearth for Kepler-20b, 10.85 d and 3.07^{+0.20}_{-0.31} Rearth for Kepelr-20c, and 77.61 d and 2.75^{+0.17}_{-0.30} Rearth for Kepler-20d. From multi-epoch radial velocities, we determine the masses of Kepler-20b and Kepler-20c to be 8.7\+-2.2 Mearth and 16.1+-3.5 Mearth, respectively, and we place an upper limit on the mass of Kepler-20d of 20.1 Mearth (2 sigma).

preprint2012arXiv

Planetary Candidates Observed by Kepler, III: Analysis of the First 16 Months of Data

New transiting planet candidates are identified in sixteen months (May 2009 - September 2010) of data from the Kepler spacecraft. Nearly five thousand periodic transit-like signals are vetted against astrophysical and instrumental false positives yielding 1,091 viable new planet candidates, bringing the total count up to over 2,300. Improved vetting metrics are employed, contributing to higher catalog reliability. Most notable is the noise-weighted robust averaging of multi-quarter photo-center offsets derived from difference image analysis which identifies likely background eclipsing binaries. Twenty-two months of photometry are used for the purpose of characterizing each of the new candidates. Ephemerides (transit epoch, T_0, and orbital period, P) are tabulated as well as the products of light curve modeling: reduced radius (Rp/R*), reduced semi-major axis (d/R*), and impact parameter (b). The largest fractional increases are seen for the smallest planet candidates (197% for candidates smaller than 2Re compared to 52% for candidates larger than 2Re) and those at longer orbital periods (123% for candidates outside of 50-day orbits versus 85% for candidates inside of 50-day orbits). The gains are larger than expected from increasing the observing window from thirteen months (Quarter 1-- Quarter 5) to sixteen months (Quarter 1 -- Quarter 6). This demonstrates the benefit of continued development of pipeline analysis software. The fraction of all host stars with multiple candidates has grown from 17% to 20%, and the paucity of short-period giant planets in multiple systems is still evident. The progression toward smaller planets at longer orbital periods with each new catalog release suggests that Earth-size planets in the Habitable Zone are forthcoming if, indeed, such planets are abundant.

preprint2012arXiv

Transit Detection in the MEarth Survey of Nearby M Dwarfs: Bridging the Clean-First, Search-Later Divide

In the effort to characterize the masses, radii, and atmospheres of potentially habitable exoplanets, there is an urgent need to find examples of such planets transiting nearby M dwarfs. The MEarth Project is an ongoing effort to do so, as a ground-based photometric survey designed to detect exoplanets as small as 2 Earth radii transiting mid-to-late M dwarfs within 33 pc of the Sun. Unfortunately, identifying transits of such planets in photometric monitoring is complicated both by the intrinsic stellar variability that is common among these stars and by the nocturnal cadence, atmospheric variations, and instrumental systematics that often plague Earth-bound observatories. Here we summarize the properties of MEarth data gathered so far, and we present a new framework to detect shallow exoplanet transits in wiggly and irregularly-spaced light curves. In contrast to previous methods that clean trends from light curves before searching for transits, this framework assesses the significance of individual transits simultaneously while modeling variability, systematics, and the photometric quality of individual nights. Our Method for Including Starspots and Systematics in the Marginalized Probability of a Lone Eclipse (MISS MarPLE) uses a computationally efficient semi-Bayesian approach to explore the vast probability space spanned by the many parameters of this model, naturally incorporating the uncertainties in these parameters into its evaluation of candidate events. We show how to combine individual transits processed by MISS MarPLE into periodic transiting planet candidates and compare our results to the popular Box-fitting Least Squares (BLS) method with simulations. By applying MISS MarPLE to observations from the MEarth Project, we demonstrate the utility of this framework for robustly assessing the false alarm probability of transit signals in real data. [slightly abridged]

preprint2012arXiv

Transit Timing Observations from Kepler: II. Confirmation of Two Multiplanet Systems via a Non-parametric Correlation Analysis

We present a new method for confirming transiting planets based on the combination of transit timingn variations (TTVs) and dynamical stability. Correlated TTVs provide evidence that the pair of bodies are in the same physical system. Orbital stability provides upper limits for the masses of the transiting companions that are in the planetary regime. This paper describes a non-parametric technique for quantifying the statistical significance of TTVs based on the correlation of two TTV data sets. We apply this method to an analysis of the transit timing variations of two stars with multiple transiting planet candidates identified by Kepler. We confirm four transiting planets in two multiple planet systems based on their TTVs and the constraints imposed by dynamical stability. An additional three candidates in these same systems are not confirmed as planets, but are likely to be validated as real planets once further observations and analyses are possible. If all were confirmed, these systems would be near 4:6:9 and 2:4:6:9 period commensurabilities. Our results demonstrate that TTVs provide a powerful tool for confirming transiting planets, including low-mass planets and planets around faint stars for which Doppler follow-up is not practical with existing facilities. Continued Kepler observations will dramatically improve the constraints on the planet masses and orbits and provide sensitivity for detecting additional non-transiting planets. If Kepler observations were extended to eight years, then a similar analysis could likely confirm systems with multiple closely spaced, small transiting planets in or near the habitable zone of solar-type stars.

preprint2012arXiv

Transit Timing Observations from Kepler: III. Confirmation of 4 Multiple Planet Systems by a Fourier-Domain Study of Anti-correlated Transit Timing Variations

We present a method to confirm the planetary nature of objects in systems with multiple transiting exoplanet candidates. This method involves a Fourier-Domain analysis of the deviations in the transit times from a constant period that result from dynamical interactions within the system. The combination of observed anti-correlations in the transit times and mass constraints from dynamical stability allow us to claim the discovery of four planetary systems Kepler-25, Kepler-26, Kepler-27, and Kepler-28, containing eight planets and one additional planet candidate.

preprint2011arXiv

A Closely-Packed System of Low-Mass, Low-Density Planets Transiting Kepler-11

When an extrasolar planet passes in front of its star (transits), its radius can be measured from the decrease in starlight and its orbital period from the time between transits. Multiple planets transiting the same star reveal more: period ratios determine stability and dynamics, mutual gravitational interactions reflect planet masses and orbital shapes, and the fraction of transiting planets observed as multiples has implications for the planarity of planetary systems. But few stars have more than one known transiting planet, and none has more than three. Here we report Kepler spacecraft observations of a single Sun-like star that reveal six transiting planets, five with orbital periods between 10 and 47 days plus a sixth one with a longer period. The five inner planets are among the smallest whose masses and sizes have both been measured, and these measurements imply substantial envelopes of light gases. The degree of coplanarity and proximity of the planetary orbits imply energy dissipation near the end of planet formation.

preprint2011arXiv

A Search for Additional Planets in Five of the Exoplanetary Systems Studied by the NASA EPOXI Mission

We present time series photometry and constraints on additional planets in five of the exoplanetary systems studied by the EPOCh (Extrasolar Planet Observation and Characterization) component of the NASA EPOXI mission: HAT-P-4, TrES-3, TrES-2, WASP-3, and HAT-P-7. We conduct a search of the high-precision time series for photometric transits of additional planets. We find no candidate transits with significance higher than our detection limit. From Monte Carlo tests of the time series using putative periods from 0.5 days to 7 days, we demonstrate the sensitivity to detect Neptune-sized companions around TrES-2, sub-Saturn-sized companions in the HAT-P-4, TrES-3, and WASP-3 systems, and Saturn-sized companions around HAT-P-7. We investigate in particular our sensitivity to additional transits in the dynamically favorable 3:2 and 2:1 exterior resonances with the known exoplanets: if we assume coplanar orbits with the known planets, then companions in these resonances with HAT-P-4b, WASP-3b, and HAT-P-7b would be expected to transit, and we can set lower limits on the radii of companions in these systems. In the nearly grazing exoplanetary systems TrES-3 and TrES-2, additional coplanar planets in these resonances are not expected to transit. However, we place lower limits on the radii of companions that would transit if the orbits were misaligned by 2.0 degrees and 1.4 degrees for TrES-3 and TrES-2, respectively.

preprint2011arXiv

A Spitzer Transmission Spectrum for the Exoplanet GJ 436b, Evidence for Stellar Variability, and Constraints on Dayside Flux Variations

In this paper we describe a uniform analysis of eight transits and eleven secondary eclipses of the extrasolar planet GJ 436b obtained in the 3.6, 4.5, and 8.0 micron bands using the IRAC instrument on the Spitzer Space Telescope between UT 2007 June 29 and UT 2009 Feb 4. We find that the best-fit transit depths for visits in the same bandpass can vary by as much as 8% of the total (4.7 sigma significance) from one epoch to the next. Although we cannot entirely rule out residual detector effects or a time-varying, high-altitude cloud layer in the planet's atmosphere as the cause of these variations, we consider the occultation of active regions on the star in a subset of the transit observations to be the most likely explanation. We reconcile the presence of magnetically active regions with the lack of significant visible or infrared flux variations from the star by proposing that the star's spin axis is tilted with respect to our line of sight, and that the planet's orbit is therefore likely to be misaligned. These observations serve to illustrate the challenges associated with transmission spectroscopy of planets orbiting late-type stars; we expect that other systems, such as GJ 1214, may display comparably variable transit depths. Our measured 8 micron secondary eclipse depths are consistent with a constant value, and we place a 1 sigma upper limit of 17% on changes in the planet's dayside flux in this band. Averaging over the eleven visits gives us an improved estimate of 0.0452% +/- 0.0027% for the secondary eclipse depth. We combine timing information from our observations with previously published data to produce a refined orbital ephemeris, and determine that the best-fit transit and eclipse times are consistent with a constant orbital period. [ABRIDGED]

preprint2011arXiv

Characteristics of planetary candidates observed by Kepler, II: Analysis of the first four months of data

On 1 February 2011 the Kepler Mission released data for 156,453 stars observed from the beginning of the science observations on 2 May through 16 September 2009. There are 1235 planetary candidates with transit like signatures detected in this period. These are associated with 997 host stars. Distributions of the characteristics of the planetary candidates are separated into five class-sizes; 68 candidates of approximately Earth-size (radius < 1.25 Earth radii), 288 super-Earth size (1.25 Earth radii < radius < 2 Earth radii), 662 Neptune-size (2 Earth radii < radius < 6 Earth radii), 165 Jupiter-size (6 Earth radii < radius < 15 Earth radii), and 19 up to twice the size of Jupiter (15 Earth radii < radius < 22 Earth radii). In the temperature range appropriate for the habitable zone, 54 candidates are found with sizes ranging from Earth-size to larger than that of Jupiter. Five are less than twice the size of the Earth. Over 74% of the planetary candidates are smaller than Neptune. The observed number versus size distribution of planetary candidates increases to a peak at two to three times Earth-size and then declines inversely proportional to area of the candidate. Our current best estimates of the intrinsic frequencies of planetary candidates, after correcting for geometric and sensitivity biases, are 6% for Earth-size candidates, 7% for super-Earth size candidates, 17% for Neptune-size candidates, and 4% for Jupiter-size candidates. Multi-candidate, transiting systems are frequent; 17% of the host stars have multi-candidate systems, and 33.9% of all the candidates are part of multi-candidate systems.

preprint2011arXiv

Discovery and Atmospheric Characterization of Giant Planet Kepler-12b: An Inflated Radius Outlier

We report the discovery of planet Kepler-12b (KOI-20), which at 1.695\pm0.030 RJ is among the handful of planets with super-inflated radii above 1.65 RJ. Orbiting its slightly evolved G0 host with a 4.438-day period, this 0.431\pm0.041 MJ planet is the least-irradiated within this largest-planet-radius group, which has important implications for planetary physics. The planet's inflated radius and low mass lead to a very low density of 0.111\pm0.010 g cm-3. We detect the occultation of the planet at a significance of 3.7σ in the Kepler bandpass. This yields a geometric albedo of 0.14\pm0.04; the planetary flux is due to a combination of scattered light and emitted thermal flux. We use multiple observations with Warm Spitzer to detect the occultation at 7σ and 4σ in the 3.6 and 4.5 μm bandpasses, respectively. The occultation photometry timing is consistent with a circular orbit, at e < 0.01 (1σ), and e < 0.09 (3σ). The occultation detections across the three bands favor an atmospheric model with no dayside temperature inversion. The Kepler occultation detection provides significant leverage, but conclusions regarding temperature structure are preliminary, given our ignorance of opacity sources at optical wavelengths in hot Jupiter atmospheres. If Kepler-12b and HD 209458b, which intercept similar incident stellar fluxes, have the same heavy element masses, the interior energy source needed to explain the large radius of Kepler-12b is three times larger than that of HD 209458b. This may suggest that more than one radius-inflation mechanism is at work for Kepler-12b, or that it is less heavy-element rich than other transiting planets.

preprint2011arXiv

Kepler 18-b, c, and d: A System Of Three Planets Confirmed by Transit Timing Variations, Lightcurve Validation, Spitzer Photometry and Radial Velocity Measurements

We report the detection of three transiting planets around a Sunlike star, which we designate Kepler-18. The transit signals were detected in photometric data from the Kepler satellite, and were confirmed to arise from planets using a combination of large transit-timing variations, radial-velocity variations, Warm-Spitzer observations, and statistical analysis of false-positive probabilities. The Kepler-18 star has a mass of 0.97M_sun, radius 1.1R_sun, effective temperature 5345K, and iron abundance [Fe/H]= +0.19. The planets have orbital periods of approximately 3.5, 7.6 and 14.9 days. The innermost planet "b" is a "super-Earth" with mass 6.9 \pm 3.4M_earth, radius 2.00 \pm 0.10R_earth, and mean density 4.9 \pm 2.4 g cm^-3. The two outer planets "c" and "d" are both low-density Neptune-mass planets. Kepler-18c has a mass of 17.3 \pm 1.9M_earth, radius 5.49 \pm 0.26R_earth, and mean density 0.59 \pm 0.07 g cm^-3, while Kepler-18d has a mass of 16.4 \pm 1.4M_earth, radius 6.98 \pm 0.33R_earth, and mean density 0.27 \pm 0.03 g cm^-3. Kepler-18c and Kepler-18d have orbital periods near a 2:1 mean-motion resonance, leading to large and readily detected transit timing variations.

preprint2011arXiv

Kepler-10c, a 2.2-Earth radius transiting planet in a multiple system

The Kepler Mission has recently announced the discovery of Kepler-10 b, the smallest exoplanet discovered to date and the first rocky planet found by the spacecraft. A second, 45-day period transit-like signal present in the photometry from the first eight months of data could not be confirmed as being caused by a planet at the time of that announcement. Here we apply the light-curve modeling technique known as BLENDER to explore the possibility that the signal might be due to an astrophysical false positive (blend). To aid in this analysis we report the observation of two transits with the Spitzer Space Telescope at 4.5 μm. When combined they yield a transit depth of 344 \pm 85 ppm that is consistent with the depth in the Kepler passband (376 \pm 9 ppm, ignoring limb darkening), which rules out blends with an eclipsing binary of a significantly different color than the target. Using these observations along with other constraints from high resolution imaging and spectroscopy we are able to exclude the vast majority of possible false positives. We assess the likelihood of the remaining blends, and arrive conservatively at a false alarm rate of 1.6 \times 10-5 that is small enough to validate the candidate as a planet (designated Kepler-10 c) with a very high level of confidence. The radius of this object is measured to be Rp = 2.227+0.052 -0.057 Earth radii. Kepler-10 c represents another example (with Kepler-9 d and Kepler-11 g) of statistical "validation" of a transiting exoplanet, as opposed to the usual "confirmation" that can take place when the Doppler signal is detected or transit timing variations are measured. It is anticipated that many of Kepler's smaller candidates will receive a similar treatment since dynamical confirmation may be difficult or impractical with the sensitivity of current instrumentation.

preprint2011arXiv

Kepler-22b: A 2.4 Earth-radius Planet in the Habitable Zone of a Sun-like Star

A search of the time-series photometry from NASA's Kepler spacecraft reveals a transiting planet candidate orbiting the 11th magnitude G5 dwarf KIC 10593626 with a period of 290 days. The characteristics of the host star are well constrained by high-resolution spectroscopy combined with an asteroseismic analysis of the Kepler photometry, leading to an estimated mass and radius of 0.970 +/- 0.060 MSun and 0.979 +/- 0.020 RSun. The depth of 492 +/- 10ppm for the three observed transits yields a radius of 2.38 +/- 0.13 REarth for the planet. The system passes a battery of tests for false positives, including reconnaissance spectroscopy, high-resolution imaging, and centroid motion. A full BLENDER analysis provides further validation of the planet interpretation by showing that contamination of the target by an eclipsing system would rarely mimic the observed shape of the transits. The final validation of the planet is provided by 16 radial velocities obtained with HIRES on Keck 1 over a one year span. Although the velocities do not lead to a reliable orbit and mass determination, they are able to constrain the mass to a 3σ upper limit of 124 MEarth, safely in the regime of planetary masses, thus earning the designation Kepler-22b. The radiative equilibrium temperature is 262K for a planet in Kepler-22b's orbit. Although there is no evidence that Kepler-22b is a rocky planet, it is the first confirmed planet with a measured radius to orbit in the Habitable Zone of any star other than the Sun.

preprint2011arXiv

LSPM J1112+7626: detection of a 41-day M-dwarf eclipsing binary from the MEarth transit survey

We report the detection of eclipses in LSPM J1112+7626, which we find to be a moderately bright (I_C = 12.14 +/- 0.05) very low-mass binary system with an orbital period of 41.03236 +/- 0.00002 days, and component masses M_1 = 0.395 +/- 0.002 Msol and M_2 = 0.275 +/- 0.001 Msol in an eccentric (e = 0.239 +/- 0.002) orbit. A 65 day out of eclipse modulation of approximately 2% peak-to-peak amplitude is seen in I-band, which is probably due to rotational modulation of photospheric spots on one of the binary components. This paper presents the discovery and characterization of the object, including radial velocities sufficient to determine both component masses to better than 1% precision, and a photometric solution. We find that the sum of the component radii, which is much better-determined than the individual radii, is inflated by 3.8 +0.9 -0.5 % compared to the theoretical model predictions, depending on the age and metallicity assumed. These results demonstrate that the difficulties in reproducing observed M-dwarf eclipsing binary radii with theoretical models are not confined to systems with very short orbital periods. This object promises to be a fruitful testing ground for the hypothesized link between inflated radii in M-dwarfs and activity.

preprint2011arXiv

Observational evidence for a metal rich atmosphere on the super-Earth GJ1214b

We report observations of two consecutive transits of the warm super-Earth exoplanet GJ1214b at 3.6 and 4.5 microns with the Infrared Array Camera instrument on-board the Spitzer Space Telescope. The two transit light curves allow for the determination of the transit parameters for this system. We find these paremeters to be consistent with the previously determined values and no evidence for transit timing variations. The main investigation consists of measuring the transit depths in each bandpass to constrain the planet's transmission spectrum. Fixing the system scale and impact parameters, we measure R_p/R_star=0.1176 (+0.0008/-0.0009) and 0.1163 (+0.0010/-0.0008) at 3.6 and 4.5 microns, respectively. Combining these data with the previously reported MEarth Observatory measurements in the red optical yields constraints on the GJ1214b's transmission spectrum and allows us to rule-out a cloud-free, solar composition (i.e., hydrogen-dominated) atmosphere at 4.5 sigma confidence. This independently confirms a recent finding that was based on a measurement of the planet's transmission spectrum using the VLT. The Spitzer, MEarth, and VLT observations together yield a remarkably flat transmission spectrum over the large wavelength domain spanned by the data. Consequently, cloud-free atmospheric models require more than 30% metals (assumed to be in the form of H2O by volume to be consistent with all the observations.

preprint2011arXiv

Planet Occurrence within 0.25 AU of Solar-type Stars from Kepler

We report the distribution of planets as a function of planet radius (R_p), orbital period (P), and stellar effective temperature (Teff) for P < 50 day orbits around GK stars. These results are based on the 1,235 planets (formally "planet candidates") from the Kepler mission that include a nearly complete set of detected planets as small as 2 Earth radii (Re). For each of the 156,000 target stars we assess the detectability of planets as a function of R_p and P. We also correct for the geometric probability of transit, R*/a. We consider first stars within the "solar subset" having Teff = 4100-6100 K, logg = 4.0-4.9, and Kepler magnitude Kp < 15 mag. We include only those stars having noise low enough to permit detection of planets down to 2 Re. We count planets in small domains of R_p and P and divide by the included target stars to calculate planet occurrence in each domain. Occurrence of planets varies by more than three orders of magnitude and increases substantially down to the smallest radius (2 Re) and out to the longest orbital period (50 days, ~0.25 AU) in our study. For P < 50 days, the radius distribution is given by a power law, df/dlogR= k R^α. This rapid increase in planet occurrence with decreasing planet size agrees with core-accretion, but disagrees with population synthesis models. We fit occurrence as a function of P to a power law model with an exponential cutoff below a critical period P_0. For smaller planets, P_0 has larger values, suggesting that the "parking distance" for migrating planets moves outward with decreasing planet size. We also measured planet occurrence over Teff = 3600-7100 K, spanning M0 to F2 dwarfs. The occurrence of 2-4 Re planets in the Kepler field increases with decreasing Teff, making these small planets seven times more abundant around cool stars than the hottest stars in our sample. [abridged]

preprint2011arXiv

Rotational Variability of Earth's Polar Regions: Implications for Detecting Snowball Planets

We have obtained the first time-resolved, disc-integrated observations of Earth's poles with the Deep Impact spacecraft as part of the EPOXI Mission of Opportunity. These data mimic what we will see when we point next-generation space telescopes at nearby exoplanets. We use principal component analysis (PCA) and rotational lightcurve inversion to characterize color inhomogeneities and map their spatial distribution from these unusual vantage points, as a complement to the equatorial views presented in Cowan et al. (2009). We also perform the same PCA on a suite of simulated rotational multi-band lightcurves from NASA's Virtual Planetary Laboratory 3D spectral Earth model. This numerical experiment allows us to understand what sorts of surface features PCA can robustly identify. We find that the EPOXI polar observations have similar broadband colors as the equatorial Earth, but with 20-30% greater apparent albedo. This is because the polar observations are most sensitive to mid-latitudes, which tend to be more cloudy than the equatorial latitudes emphasized by the original EPOXI Earth observations. The cloudiness of the mid-latitudes also manifests itself in the form of increased variability at short wavelengths in the polar observations, and as a dominant gray eigencolor in the south polar observation. We construct a simple reflectance model for a snowball Earth. By construction, our model has a higher Bond albedo than the modern Earth; its surface albedo is so high that Rayleigh scattering does not noticeably affect its spectrum. The rotational color variations occur at short wavelengths due to the large contrast between glacier ice and bare land in those wavebands. Thus we find that both the broadband colors and diurnal color variations of such a planet would be easily distinguishable from the modern-day Earth, regardless of viewing angle.

preprint2011arXiv

Spitzer Infrared Observations and Independent Validation of the Transiting Super-Earth CoRoT-7b

The detection and characterization of the first transiting super-Earth, CoRoT-7 b, has required an unprecedented effort in terms of telescope time and analysis. Although the star does display a radial velocity signal at the period of the planet, this has been difficult to disentangle from the intrinsic stellar variability, and pinning down the velocity amplitude has been very challenging. As a result, the precise value of the mass of the planet - and even the extent to which it can be considered to be confirmed - have been debated in the recent literature, with six mass measurements published so far based on the same spectroscopic observations, ranging from about 2 to 8 Earth masses. Here we report on an independent validation of the planet discovery, using one of the fundamental properties of a transit signal: its achromaticity. We observed four transits of CoRoT-7 b with Spitzer, in order to determine whether the depth of the transit signal in the near-infrared is consistent with that observed in the CoRoT bandpass, as expected for a planet. We detected the transit and found an average depth of 0.426 {\pm} 0.115 mmag at 4.5 μm, which is in good agreement with the depth of 0.350 {\pm} 0.011 mmag found by CoRoT. These observations place important constraints on the kinds of astrophysical false positives that could mimic the signal. Combining this with additional constraints reported earlier, we performed an exhaustive exploration of possible blends scenarios for CoRoT-7 b using the BLENDER technique. We are able to rule out the vast majority of false positives, and the remaining ones are found to be much less likely than a true transiting planet. We thus validate CoRoT-7 b as a bona-fide planet with a very high degree of confidence, independently of any radial-velocity information. Our Spitzer observations have additionally allowed us to significantly improve the ephemeris of the planet.

preprint2011arXiv

The atmospheres of the hot-Jupiters Kepler-5b and Kepler-6b observed during occultations with Warm-Spitzer and Kepler

This paper reports the detection and the measurements of occultations of the two transiting hot giant exoplanets Kepler-5b and Kepler-6b by their parent stars. The observations are obtained in the near infrared with Spitzer Space Telescope and at optical wavelengths by combining more than a year of Kepler photometry. The investigation consists of constraining the eccentricities of these systems and of obtaining broad band emergent spectra for individual planets. For both targets, the occultations are detected at 3 sigma level at each wavelength with mid-occultation times consistent with circular orbits. The brightness temperatures of these planets are deduced from the infrared observations and reach T=1930+/-100K and T=1660+/-120K for Kepler-5b and Kepler-6b respectively. We measure optical geometric albedos A_g in the Kepler bandpass and find A_g=0.12+/-0.04 for Kepler-5b and A_g=0.11+/-0.04 for Kepler-6b leading to an upper limit for the Bond albedo of A_B < 0.17 in both cases. The observations for both planets are best described by models for which most of the incident energy is redistributed on the dayside, with only less than 10% of the absorbed stellar flux redistributed to the night side of these planets. The data for Kepler-5b favor a model without a temperature inversion, whereas for Kepler-6b they do not allow distinguishing between models with and without inversion.

preprint2011arXiv

The Flat Transmission Spectrum of the Super-Earth GJ1214b from Wide Field Camera 3 on the Hubble Space Telescope

Capitalizing on the observational advantage offered by its tiny M dwarf host, we present HST/WFC3 grism measurements of the transmission spectrum of the super-Earth exoplanet GJ1214b. These are the first published WFC3 observations of a transiting exoplanet atmosphere. After correcting for a ramp-like instrumental systematic, we achieve nearly photon-limited precision in these observations, finding the transmission spectrum of GJ1214b to be flat between 1.1 and 1.7 microns. Inconsistent with a cloud-free solar composition atmosphere at 8.2 sigma, the measured achromatic transit depth most likely implies a large mean molecular weight for GJ1214b's outer envelope. A dense atmosphere rules out bulk compositions for GJ1214b that explain its large radius by the presence of a very low density gas layer surrounding the planet. High-altitude clouds can alternatively explain the flat transmission spectrum, but they would need to be optically thick up to 10 mbar or consist of particles with a range of sizes approaching 1 micron in diameter.

preprint2011arXiv

The GJ1214 Super-Earth System: Stellar Variability, New Transits, and a Search for Additional Planets

The super-Earth GJ1214b transits a nearby M dwarf that exhibits 1% intrinsic variability in the near-infrared. Here, we analyze new observations to refine the physical properties of both the star and planet. We present three years of out-of-transit photometric monitoring of the stellar host GJ1214 from the MEarth Observatory and find the rotation period to be long, mostly likely an integer multiple of 53 days, suggesting low levels of magnetic activity and an old age for the system. We show such variability will not pose significant problems to ongoing studies of the planet's atmosphere with transmission spectroscopy. We analyze 2 high-precision transit light curves from ESO's Very Large Telescope along with 7 others from the MEarth and FLWO 1.2 meter telescopes, finding physical parameters for the planet that are consistent with previous work. The VLT light curves show tentative evidence for spot occultations during transit. Using two years of MEarth light curves, we place limits on additional transiting planets around GJ1214 with periods out to the habitable zone of the system. We also improve upon the previous photographic V -band estimate for the star, finding V = 14.71 \pm 0.03.

preprint2011arXiv

The hot-Jupiter Kepler-17b: discovery, obliquity from stroboscopic starspots, and atmospheric characterization

This paper reports the discovery and characterization of the transiting hot giant exoplanet Kepler-17b. The planet has an orbital period of 1.486 days, and radial velocity measurements from the Hobby-Eberly Telescope (HET) show a Doppler signal of 420+/-15 m.s-1. From a transit-based estimate of the host star's mean density, combined with an estimate of the stellar effective temperature T_eff=5630+/-100 K from high-resolution spectra, we infer a stellar host mass of 1.061+/-0.067 M_sun and a stellar radius of 1.019+/-0.033 R_jup. We estimate the planet mass and radius to be Mp=2.450+/-0.114 M_jup and Rp=1.312+/-0.018 R_jup and a planet density near 1.35 g.cm-3. The host star is active, with dark spots that are frequently occulted by the planet. The continuous monitoring of the star reveals a stellar rotation period of 11.89 days, 8 times the the planet's orbital period; this period ratio produces stroboscopic effects on the occulted starspots. The temporal pattern of these spot-crossing events shows that the planet's orbit is prograde and the star's obliquity is smaller than 15 deg. We detected planetary occultations of Kepler-17b with both the Kepler and Spitzer Space Telescopes. We use these observations to constrain the eccentricity, e, and find that it is consistent with a circular orbit (e<0.0011). The brightness temperatures of the planet the infrared bandpasses are T_3.6um=1880+/-100 K and T4.5um=1770+/-150 K. We measure the optical geometric albedo A_g in the Kepler bandpass and find A_g = 0.10+/-0.02. The observations are best described by atmospheric models for which most of the incident energy is re-radiated away from the day side.

preprint2011arXiv

The Kepler-19 System: A Transiting 2.2 R_Earth Planet and a Second Planet Detected via Transit Timing Variations

We present the discovery of the Kepler-19 planetary system, which we first identified from a 9.3-day periodic transit signal in the Kepler photometry. From high-resolution spectroscopy of the star, we find a stellar effective temperature Teff=5541 \pm 60 K, a metallicity [Fe/H]=-0.13 \pm 0.06, and a surface gravity log(g)=4.59 \pm 0.10. We combine the estimate of Teff and [Fe/H] with an estimate of the stellar density derived from the photometric light curve to deduce a stellar mass of M_star = 0.936 \pm 0.040 M_Sun and a stellar radius of R_star = 0.850 \pm 0.018 R_Sun. We rule out the possibility that the transits result from an astrophysical false positive by first identifying the subset of stellar blends that reproduce the precise shape of the light curve. We conclude that the planetary scenario is more than three orders of magnitude more likely than a blend. The blend scenario is independently disfavored by the achromaticity of the transit: we measure a transit depth with Spitzer at 4.5 μm of 547+113-110 ppm, consistent with the depth measured in the Kepler optical bandpass of 567\pm6 ppm. We determine a physical radius of the planet Kepler-19b of R_p = 2.209 \pm 0.048 R_Earth. From radial-velocity observations of the star, we find an upper limit on the planet mass of 20.3 M_Earth, corresponding to a maximum density of 10.4 g cm^-3. We report a significant sinusoidal deviation of the transit times from a predicted linear ephemeris, which we conclude is due to an additional perturbing body in the system. We cannot uniquely determine the orbital parameters of the perturber, as various dynamical mechanisms match the amplitude, period, and shape of the transit timing signal and satisfy the host star's radial velocity limits. However, the perturber in these mechanisms has period <160 days and mass <6 M_Jup, confirming its planetary nature as Kepler-19c. [Abridged]

preprint2011arXiv

TrES-5: A Massive Jupiter-sized Planet Transiting A Cool G-dwarf

We report the discovery of TrES-5, a massive hot Jupiter that transits the star GSC 03949-00967 every 1.48 days. From spectroscopy of the star we estimate a stellar effective temperature of$5171 +/- 36 K, and from high-precision B, R and I photometry of the transit we constrain the ratio of the semi-major axis and the stellar radius to be 6.07 +/- 0.14. We compare these values to model stellar isochrones to obtain a stellar mass of 0.893 +/- 0.024 solar masses. Based on this estimate and the photometric time series, we constrain the stellar radius to be 0.866 +/- 0.013 solar radii, and the planet radius to be 1.209 +/- 0.021 Jupiter radii. We model our radial-velocity data assuming a circular orbit and find a planetary mass of 1.778 +/- 0.063 Jupiter masses. Our radial-velocity observations rule out line-bisector variations that would indicate a specious detection resulting from a blend of an eclipsing binary system. TrES-5 orbits one of the faintest stars with transiting planets found to date from the ground and demonstrates that precise photometry and followup spectroscopy are possible, albeit challenging, even for such faint stars.

preprint2011arXiv

Two Earth-sized planets orbiting Kepler-20

Since the discovery of the first extrasolar giant planets around Sun-like stars, evolving observational capabilities have brought us closer to the detection of true Earth analogues. The size of an exoplanet can be determined when it periodically passes in front of (transits) its parent star, causing a decrease in starlight proportional to its radius. The smallest exoplanet hitherto discovered has a radius 1.42 times that of the Earth's radius (R Earth), and hence has 2.9 times its volume. Here we report the discovery of two planets, one Earth-sized (1.03R Earth) and the other smaller than the Earth (0.87R Earth), orbiting the star Kepler-20, which is already known to host three other, larger, transiting planets. The gravitational pull of the new planets on the parent star is too small to measure with current instrumentation. We apply a statistical method to show that the likelihood of the planetary interpretation of the transit signals is more than three orders of magnitude larger than that of the alternative hypothesis that the signals result from an eclipsing binary star. Theoretical considerations imply that these planets are rocky, with a composition of iron and silicate. The outer planet could have developed a thick water vapour atmosphere.

preprint2011arXiv

Warm Spitzer Photometry of XO-4b, HAT-P-6b and HAT-P-8b

We have analyzed Warm Spitzer/IRAC observations of the secondary eclipses of three planets, XO-4b, HAT-P-6b and HAT-P-8b. We measure secondary eclipse amplitudes at 3.6μm and 4.5μm for each target. XO-4b exhibits a stronger eclipse depth at 4.5μm than at 3.6μm, which is consistent with the presence of a temperature inversion. HAT-P-8b shows a stronger eclipse amplitude at 3.6μm, and is best-described by models without a temperature inversion. The eclipse depths of HAT-P-6b can be fitted with models with a small or no temperature inversion. We consider our results in the context of a postulated relationship between stellar activity and temperature inversions and a relationship between irradiation level and planet dayside temperature, as discussed by Knutson et al. (2010) and Cowan & Agol (2011), respectively. Our results are consistent with these hypotheses, but do not significantly strengthen them. To measure accurate secondary eclipse central phases, we require accurate ephemerides. We obtain primary transit observations and supplement them with publicly available observations to update the orbital ephemerides of the three planets. Based on the secondary eclipse timing, we set upper boundaries for e cos(ω) for HAT-P-6b, HAT-P-8b and XO-4b and find that the values are consistent with circular orbits.

preprint2010arXiv

A Search for a Sub-Earth Sized Companion to GJ 436 and a Novel Method to Calibrate Warm Spitzer IRAC Observations

We discovered evidence for a possible additional 0.75 R_Earth transiting planet in the NASA EPOXI observations of the known M dwarf exoplanetary system GJ 436. Based on an ephemeris determined from the EPOXI data, we predicted a transit event in an extant Spitzer Space Telescope 8 micron data set of this star. Our subsequent analysis of those Spitzer data confirmed the signal of the predicted depth and at the predicted time, but we found that the transit depth was dependent on the aperture used to perform the photometry. Based on these suggestive findings, we gathered new Warm Spitzer Observations of GJ 436 at 4.5 microns spanning a time of transit predicted from the EPOXI and Spitzer 8 micron candidate events. The 4.5 micron data permit us to rule out a transit at high confidence, and we conclude that the earlier candidate transit signals resulted from correlated noise in the EPOXI and Spitzer 8 micron observations. In the course of this investigation, we developed a novel method for correcting the intrapixel sensitivity variations of the 3.6 and 4.5 micron channels of the Infrared Array Camera (IRAC) instrument. We demonstrate the sensitivity of Warm Spitzer observations of M dwarfs to confirm sub-Earth sized planets. Our analysis will inform similar work that will be undertaken to use Warm Spitzer observations to confirm rocky planets discovered by the Kepler mission.

preprint2010arXiv

A Search for Additional Planets in the NASA EPOXI Observations of the Exoplanet System GJ 436

We present time series photometry of the M dwarf transiting exoplanet system GJ 436 obtained with the the EPOCh (Extrasolar Planet Observation and Characterization) component of the NASA EPOXI mission. We conduct a search of the high-precision time series for additional planets around GJ 436, which could be revealed either directly through their photometric transits, or indirectly through the variations these second planets induce on the transits of the previously known planet. In the case of GJ 436, the presence of a second planet is perhaps indicated by the residual orbital eccentricity of the known hot Neptune companion. We find no candidate transits with significance higher than our detection limit. From Monte Carlo tests of the time series, we rule out transiting planets larger than 1.5 R_Earth interior to GJ 436b with 95% confidence, and larger than 1.25 R_Earth with 80% confidence. Assuming coplanarity of additional planets with the orbit of GJ 436b, we cannot expect that putative planets with orbital periods longer than about 3.4 days will transit. However, if such a planet were to transit, we rule out planets larger than 2.0 R_Earth with orbital periods less than 8.5 days with 95% confidence. We also place dynamical constraints on additional bodies in the GJ 436 system. Our analysis should serve as a useful guide for similar analyses for which radial velocity measurements are not available, such as those discovered by the Kepler mission. These dynamical constraints on additional planets with periods from 0.5 to 9 days rule out coplanar secular perturbers as small as 10 M_Earth and non-coplanar secular perturbers as small as 1 M_Earth in orbits close in to GJ 436b. We present refined estimates of the system parameters for GJ 436. We also report a sinusoidal modulation in the GJ 436 light curve that we attribute to star spots. [Abridged]

preprint2010arXiv

Modeling Kepler transit light curves as false positives: Rejection of blend scenarios for Kepler-9, and validation of Kepler-9d, a super-Earth-size planet in a multiple system

Light curves from the Kepler Mission contain valuable information on the nature of the phenomena producing the transit-like signals. To assist in exploring the possibility that they are due to an astrophysical false positive, we describe a procedure (BLENDER) to model the photometry in terms of a "blend" rather than a planet orbiting a star. A blend may consist of a background or foreground eclipsing binary (or star-planet pair) whose eclipses are attenuated by the light of the candidate and possibly other stars within the photometric aperture. We apply BLENDER to the case of Kepler-9, a target harboring two previously confirmed Saturn-size planets (Kepler-9b and Kepler-9c) showing transit timing variations, and an additional shallower signal with a 1.59-day period suggesting the presence of a super-Earth-size planet. Using BLENDER together with constraints from other follow-up observations we are able to rule out all blends for the two deeper signals, and provide independent validation of their planetary nature. For the shallower signal we rule out a large fraction of the false positives that might mimic the transits. The false alarm rate for remaining blends depends in part (and inversely) on the unknown frequency of small-size planets. Based on several realistic estimates of this frequency we conclude with very high confidence that this small signal is due to a super-Earth-size planet (Kepler-9d) in a multiple system, rather than a false positive. The radius is determined to be 1.64 (+0.19/-0.14) R(Earth), and current spectroscopic observations are as yet insufficient to establish its mass.

preprint2010arXiv

NLTT 41135: a field M-dwarf + brown dwarf eclipsing binary in a triple system, discovered by the MEarth observatory

We report the discovery of an eclipsing companion to NLTT 41135, a nearby M5 dwarf that was already known to have a wider, slightly more massive common proper motion companion, NLTT 41136, at 2.4 arcsec separation. Analysis of combined-light and radial velocity curves of the system indicates that NLTT 41135B is a 31-34 +/- 3 MJup brown dwarf (where the range depends on the unknown metallicity of the host star) on a circular orbit. The visual M-dwarf pair appears to be physically bound, so the system forms a hierarchical triple, with masses approximately in the ratio 8:6:1. The eclipses are grazing, preventing an unambiguous measurement of the secondary radius, but follow-up observations of the secondary eclipse (e.g. with the James Webb Space Telescope) could permit measurements of the surface brightness ratio between the two objects, and thus place constraints on models of brown dwarfs.

preprint2010arXiv

On the angular momentum evolution of fully-convective stars: rotation periods for field M-dwarfs from the MEarth transit survey

We present rotation period measurements for 41 field M-dwarfs, all of which have masses inferred (from their parallaxes and 2MASS K-band magnitudes) to be between the hydrogen burning limit and 0.35 Msol, and thus should remain fully-convective throughout their lifetimes. We measure a wide range of rotation periods, from 0.28 days to 154 days, with the latter commensurate with the typical sensitivity limit of our observations. Using kinematics as a proxy for age, we find that the majority of objects likely to be thick disk or halo members (and hence, on average, older) rotate very slowly, with a median period of 92 days, compared to 0.7 days for those likely to be thin disk members (on average, younger), although there are still some rapid rotators in the thick disk sample. When combined with literature measurements for M-dwarfs, these results indicate an increase in spin-down times with decreasing stellar mass, in agreement with previous work, and that the spin-down time becomes comparable to the age of the thick disk sample below the fully-convective boundary. We additionally infer that the spin-down must remove a substantial amount of angular momentum once it begins in order to produce the slow rotators we observe in the thick disk candidates, suggesting that fully-convective M-dwarfs may still experience strong winds.

preprint2010arXiv

Precise Estimates of the Physical Parameters for the Exoplanet System HD-17156 Enabled by HST FGS Transit and Asteroseismic Observations

We present observations of three distinct transits of HD 17156b obtained with the Fine Guidance Sensors (FGS) on board the Hubble Space Telescope} (HST). We analyzed both the transit photometry and previously published radial velocities to find the planet-star radius ratio R_p/R_s = 0.07454 +/- 0.00035, inclination i=86.49 +0.24/-0.20 deg, and scaled semi-major axis a/R = 23.19 +0.32/-0.27. This last value translates directly to a mean stellar density determination of 0.522 +0.021/-0.018 g cm^-3. Analysis of asteroseismology observations by the companion paper of Gilliland et al. (2009) provides a consistent but significantly refined measurement of the stellar mean density. We compare stellar isochrones to this density estimate and find M_s = 1.275 +/- 0.018 M_sun and a stellar age of $3.37 +0.20/-0.47 Gyr. Using this estimate of M_s and incorporating the density constraint from asteroseismology, we model both the photometry and published radial velocities to estimate the planet radius R_p= 1.0870 +/- 0.0066 Jupiter radii and the stellar radius R_s = 1.5007 +/- 0.0076 R_sun. The planet radius is larger than that found in previous studies and consistent with theoretical models of a solar-composition gas giant of the same mass and equilibrium temperature. For the three transits, we determine the times of mid-transit to a precision of 6.2 s, 7.6 s, and 6.9 s, and the transit times for HD 17156 do not show any significant departures from a constant period. The joint analysis of transit photometry and asteroseismology presages similar studies that will be enabled by the NASA Kepler Mission.

preprint2010arXiv

Secondary Eclipse Photometry of WASP-4b with Warm Spitzer

We present photometry of the giant extrasolar planet WASP-4b at 3.6 and 4.5 micron taken with the Infrared Array Camera on board the Spitzer Space Telescope as part of Spitzer's extended warm mission. We find secondary eclipse depths of 0.319+/-0.031% and 0.343+/-0.027% for the 3.6 and 4.5 micron bands, respectively and show model emission spectra and pressure-temperature profiles for the planetary atmosphere. These eclipse depths are well fit by model emission spectra with water and other molecules in absorption, similar to those used for TrES-3 and HD 189733b. Depending on our choice of model, these results indicate that this planet has either a weak dayside temperature inversion or no inversion at all. The absence of a strong thermal inversion on this highly irradiated planet is contrary to the idea that highly irradiated planets are expected to have inversions, perhaps due the presence of an unknown absorber in the upper atmosphere. This result might be explained by the modestly enhanced activity level of WASP-4b's G7V host star, which could increase the amount of UV flux received by the planet, therefore reducing the abundance of the unknown stratospheric absorber in the planetary atmosphere as suggested in Knutson et al. (2010). We also find no evidence for an offset in the timing of the secondary eclipse and place a 2 sigma upper limit on |ecos(omega)| of 0.0024, which constrains the range of tidal heating models that could explain this planet's inflated radius.

preprint2010arXiv

System parameters, transit times and secondary eclipse constraints of the exoplanet systems HAT-P-4, TrES-2, TrES-3 and WASP-3 from the NASA EPOXI Mission of Opportunity

As part of the NASA EPOXI Mission of Opportunity, we observed seven known transiting extrasolar planet systems in order to construct time series photometry of extremely high phase coverage and precision. Here we present the results for four "hot-Jupiter systems" with near-solar stars - HAT-P-4, TrES-3, TrES-2 and WASP-3. We observe ten transits of HAT-P-4, estimating the planet radius Rp = 1.332 \pm 0.052 RJup, the stellar radius R \star = 1.602 \pm 0.061 R \odot, the inclination i = 89.67 \pm 0.30 degrees and the transit duration from first to fourth contact T = 255.6 \pm 1.9 minutes. For TrES-3, we observe seven transits, and find Rp = 1.320 \pm 0.057 RJup, R\star = 0.817 \pm 0.022 R\odot, i = 81.99 \pm 0.30 degrees and T = 81.9 \pm 1.1 minutes. We also note a long term variability in the TrES-3 light curve, which may be due to star spots. We observe nine transits of TrES-2, and find Rp = 1.169 \pm 0.034 RJup, R\star = 0.940 \pm 0.026 R\odot, i = 84.15 \pm 0.16 degrees and T = 107.3 \pm 1.1 minutes. Finally we observe eight transits of WASP-3, finding Rp = 1.385 \pm 0.060 RJup, R\star = 1.354 \pm 0.056 R\odot, i = 84.22 \pm 0.81 degrees and T = 167.3 \pm 1.3 minutes. We present refined orbital periods and times of transit for each target. We state 95% confidence upper limits on the secondary eclipse depths in our broadband visible bandpass centered on 650 nm. These limits are 0.073% for HAT-P-4, 0.062% for TrES-3, 0.16% for TrES-2 and 0.11% for WASP-3. We combine the TrES-3 secondary eclipse information with the existing published data and confirm that the atmosphere likely does not have a temperature inversion.

preprint2010arXiv

The NIRSPEC Ultracool Dwarf Radial Velocity Survey

We report the results of an infrared Doppler survey designed to detect brown dwarf and giant planetary companions to a magnitude-limited sample of ultracool dwarfs. Using the NIRSPEC spectrograph on the Keck II telescope, we obtained approximately 600 radial velocity measurements over a period of six years for a sample of 59 late-M and L dwarfs spanning spectral types M8/L0 to L6. A subsample of 46 of our targets have been observed on three or more epochs. We rely on telluric CH4 absorption features in the Earth's atmosphere as a simultaneous wavelength reference and exploit the rich set of CO absorption features found in the K-band spectra of cool stars and brown dwarfs to measure radial velocities and projected rotational velocities. For a bright, slowly rotating M dwarf standard we demonstrate a radial velocity precision of 50 m/s, and for slowly rotating L dwarfs we achieve a typical radial velocity precision of approximately 200 m/s. This precision is sufficient for the detection of close-in giant planetary companions to mid-L dwarfs as well as more equal mass spectroscopic binary systems with small separations (a<2 AU). We present an orbital solution for the subdwarf binary LSR1610-0040 as well as an improved solution for the M/T binary 2M0320-04. We also combine our radial velocity measurements with distance estimates and proper motions from the literature to estimate the dispersion of the space velocities of the objects in our sample. Using a kinematic age estimate we conclude that our UCDs have an age of 5.0+0.7-0.6 Gyr, similar to that of nearby sun-like stars. We simulate the efficiency with which we detect spectroscopic binaries and find that the rate of tight (a<1 AU) binaries in our sample is 2.5+8.6-1.6%, consistent with recent estimates in the literature of a tight binary fraction of 3-4%. (abridged)

preprint2010arXiv

Warm Spitzer Photometry of the Transiting Exoplanets CoRoT-1 and CoRoT-2 at Secondary Eclipse

We measure secondary eclipses of the hot giant exoplanets CoRoT-1 at 3.6 and 4.5 microns, and CoRoT-2 at 3.6 microns, both using Warm Spitzer. We find that the Warm Spitzer mission is working very well for exoplanet science. For consistency of our analysis we also re-analyze archival cryogenic Spitzer data for secondary eclipses of CoRoT-2 at 4.5 and 8 microns. We compare the total data for both planets, including optical eclipse measurements by the CoRoT mission, and ground-based eclipse measurements at 2 microns, to existing models. Both planets exhibit stronger eclipses at 4.5 than at 3.6 microns, which is often indicative of an atmospheric temperature inversion. The spectrum of CoRoT-1 is best reproduced by a 2460K blackbody, due either to a high altitude layer that strongly absorbs stellar irradiance, or an isothermal region in the planetary atmosphere. The spectrum of CoRoT-2 is unusual because the 8 micron contrast is anomalously low. Non-inverted atmospheres could potentially produce the CoRoT-2 spectrum if the planet exhibits line emission from CO at 4.5 microns, caused by tidal-induced mass loss. However, the viability of that hypothesis is questionable because the emitting region cannot be more than about 30-percent larger than the planet's transit radius, based on the ingress and egress times at eclipse. An alternative possibility to account for the spectrum of CoRoT-2 is an additional opacity source that acts strongly at wavelengths less than 5 microns, heating the upper atmosphere while allowing the deeper atmosphere seen at 8 microns to remain cooler. We obtain a similar result as Gillon et al. for the phase of the secondary eclipse of CoRoT-2, implying an eccentric orbit with e*cos(omega)=-0.0030 +/- 0.0004.

preprint2009arXiv

A super-Earth transiting a nearby low-mass star

A decade ago, the detection of the first transiting extrasolar planet provided a direct constraint on its composition and opened the door to spectroscopic investigations of extrasolar planetary atmospheres. As such characterization studies are feasible only for transiting systems that are both nearby and for which the planet-to-star radius ratio is relatively large, nearby small stars have been surveyed intensively. Doppler studies and microlensing have uncovered a population of planets with minimum masses of 1.9-10 times the Earth's mass (M_Earth), called super-Earths. The first constraint on the bulk composition of this novel class of planets was afforded by CoRoT-7b, but the distance and size of its star preclude atmospheric studies in the foreseeable future. Here we report observations of the transiting planet GJ 1214b, which has a mass of 6.55 M_Earth and a radius 2.68 times Earth's radius (R_Earth), indicating that it is intermediate in stature between Earth and the ice giants of the Solar System. We find that the planetary mass and radius are consistent with a composition of primarily water enshrouded by a hydrogen-helium envelope that is only 0.05% of the mass of the planet. The atmosphere is probably escaping hydrodynamically, indicating that it has undergone significant evolution during its history. As the star is small and only 13 parsecs away, the planetary atmosphere is amenable to study with current observatories.

preprint2009arXiv

Atmospheric circulation of hot Jupiters: Coupled radiative-dynamical general circulation model simulations of HD 189733b and HD 209458b

We present global, three-dimensional numerical simulations of HD 189733b and HD 209458b that couple the atmospheric dynamics to a realistic representation of non-gray cloud-free radiative transfer. The model, which we call the Substellar and Planetary Atmospheric Radiation and Circulation (SPARC) model, adopts the MITgcm for the dynamics and uses the radiative model of McKay, Marley, Fortney, and collaborators for the radiation. Like earlier work with simplified forcing, our simulations develop a broad eastward equatorial jet, mean westward flow at higher latitudes, and substantial flow over the poles at low pressure. For HD 189733b, our simulations without TiO and VO opacity can explain the broad features of the observed 8 and 24-micron light curves, including the modest day-night flux variation and the fact that the planet/star flux ratio peaks before the secondary eclipse. Our simulations also provide reasonable matches to the Spitzer secondary-eclipse depths at 4.5, 5.8, 8, 16, and 24 microns and the groundbased upper limit at 2.2 microns. However, we substantially underpredict the 3.6-micron secondary-eclipse depth, suggesting that our simulations are too cold in the 0.1-1 bar region. Predicted temporal variability in secondary-eclipse depths is ~1% at Spitzer bandpasses, consistent with recent observational upper limits at 8 microns. We also show that nonsynchronous rotation can significantly alter the jet structure. For HD 209458b, we include TiO and VO opacity; these simulations develop a hot (>2000 K) dayside stratosphere. Despite this stratosphere, we do not reproduce current Spitzer photometry of this planet. Light curves in Spitzer bandpasses show modest phase variation and satisfy the observational upper limit on day-night phase variation at 8 microns. (abridged)

preprint2009arXiv

Detection of Planetary Emission from the Exoplanet TrES-2 using Spitzer /IRAC

We present here the results of our observations of TrES-2 using the Infrared Array Camera on Spitzer. We monitored this transiting system during two secondary eclipses, when the planetary emission is blocked by the star. The resulting decrease in flux is 0.127%+-0.021%, 0.230%+-0.024%, 0.199%+-0.054%, and 0.359%+-0.060%, at 3.6 microns, 4.5 microns, 5.8 microns, and 8.0 microns, respectively. We show that three of these flux contrasts are well fit by a black body spectrum with T_{eff}=1500 K, as well as by a more detailed model spectrum of a planetary atmosphere. The observed planet-to-star flux ratios in all four IRAC channels can be explained by models with and without a thermal inversion in the atmosphere of TrES-2, although with different atmospheric chemistry. Based on the assumption of thermochemical equilibrium, the chemical composition of the inversion model seems more plausible, making it a more favorable scenario. TrES-2 also falls in the category of highly irradiated planets which have been theoretically predicted to exhibit thermal inversions. However, more observations at infrared and visible wavelengths would be needed to confirm a thermal inversion in this system. Furthermore, we find that the times of the secondary eclipses are consistent with previously published times of transit and the expectation from a circular orbit. This implies that TrES-2 most likely has a circular orbit, and thus does not obtain additional thermal energy from tidal dissipation of a non-zero orbital eccentricity, a proposed explanation for the large radius of this planet.

preprint2009arXiv

GJ 3236: a new bright, very low-mass eclipsing binary system discovered by the MEarth observatory

We report the detection of eclipses in GJ 3236, a bright (I = 11.6) very low mass binary system with an orbital period of 0.77 days. Analysis of light- and radial velocity curves of the system yielded component masses of 0.38 +/- 0.02 and 0.28 +/- 0.02 Msol. The central values for the stellar radii are larger than the theoretical models predict for these masses, in agreement with the results for existing eclipsing binaries, although the present 5% observational uncertainties limit the significance of the larger radii to approximately 1 sigma. Degeneracies in the light curve models resulting from the unknown configuration of surface spots on the components of GJ 3236 currently dominate the uncertainties in the radii, and could be reduced by obtaining precise, multi-band photometry covering the full orbital period. The system appears to be tidally synchronized and shows signs of high activity levels as expected for such a short orbital period, evidenced by strong Halpha emission lines in the spectra of both components. These observations probe an important region of mass-radius parameter space around the predicted transition to fully-convective stellar interiors, where there are a limited number of precise measurements available in the literature.

preprint2009arXiv

Studying the atmosphere of the exoplanet HAT-P-7b via secondary eclipse measurements with EPOXI, Spitzer and Kepler

The highly irradiated transiting exoplanet, HAT-P-7b, currently provides one of the best opportunities for studying planetary emission in the optical and infrared wavelengths. We observe six near-consecutive secondary eclipses of HAT-P-7b at optical wavelengths with the EPOXI spacecraft. We place an upper limit on the relative eclipse depth of 0.055% (95% confidence). We also analyze Spitzer observations of the same target in the infrared, obtaining secondary eclipse depths of 0.098+/-0.017%, 0.159+/-0.022%, 0.245+/-0.031% and 0.225+/-0.052% in the 3.6, 4.5, 5.8 and 8.0 micron IRAC bands respectively. We combine these measurements with the recently published Kepler secondary eclipse measurement, and generate atmospheric models for the day-side of the planet that are consistent with both the optical and infrared measurements. The data are best fit by models with a temperature inversion, as expected from the high incident flux. The models predict a low optical albedo of ~< 0.13, with subsolar abundances of Na, K, TiO and VO. We also find that the best fitting models predict that 10% of the absorbed stellar flux is redistributed to the night side of the planet, which is qualitatively consistent with the inefficient day-night redistribution apparent in the Kepler phase curve. Models without thermal inversions fit the data only at the 1.25 sigma level, and also require an overabundance of methane, which is not expected in the very hot atmosphere of HAT-P-7b. We also analyze the eight transits of HAT-P-7b present in the EPOXI dataset and improve the constraints on the system parameters, finding a period of P = 2.2047308+/-0.0000025 days, a stellar radius of R* = 1.824+/-0.089Rsun, a planetary radius of Rp = 1.342+/-0.068RJup and an inclination of i = 85.7+3.5-2.2 deg.

preprint2009arXiv

The 8 Micron Phase Variation of the Hot Saturn HD 149026b

We monitor the star HD 149026 and its Saturn-mass planet at 8.0 micron over slightly more than half an orbit using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We find an increase of 0.0227% +/- 0.0066% (3.4 sigma significance) in the combined planet-star flux during this interval. The minimum flux from the planet is 45% +/- 19% of the maximum planet flux, corresponding to a difference in brightness temperature of 480 +/- 140 K between the two hemispheres. We derive a new secondary eclipse depth of 0.0411% +/- 0.0076% in this band, corresponding to a dayside brightness temperature of 1440 +/- 150 K. Our new secondary eclipse depth is half that of a previous measurement (3.0 sigma difference) in this same bandpass by Harrington et al. (2007). We re-fit the Harrington et al. (2007) data and obtain a comparably good fit with a smaller eclipse depth that is consistent with our new value. In contrast to earlier claims, our new eclipse depth suggests that this planet's dayside emission spectrum is relatively cool, with an 8 micron brightness temperature that is less than the maximum planet-wide equilibrium temperature. We measure the interval between the transit and secondary eclipse and find that that the secondary eclipse occurs 20.9 +7.2 / -6.5 minutes earlier (2.9 sigma) than predicted for a circular orbit, a marginally significant result. This corresponds to e*cos(omega) = -0.0079 +0.0027 / -0.0025 where e is the planet's orbital eccentricity and omega is the argument of pericenter.

preprint2009arXiv

The Broadband Infrared Emission Spectrum of the Exoplanet TrES-3

We use the Spitzer Space Telescope to estimate the dayside thermal emission of the exoplanet TrES-3 integrated in the 3.6, 4.5, 5.8, and 8.0 micron bandpasses of the Infrared Array Camera (IRAC) instrument. We observe two secondary eclipses and find relative eclipse depths of 0.00346 +/- 0.00035, 0.00372 +/- 0.00054, 0.00449 +/- 0.00097, and 0.00475 +/- 0.00046, respectively in the 4 IRAC bandpasses. We combine our results with the earlier K band measurement of De Mooij et al. (2009), and compare them with models of the planetary emission. We find that the planet does not require the presence of an inversion layer in the high atmosphere. This is the first very strongly irradiated planet that does not have a temperature inversion, which indicates that stellar or planetary characteristics other than temperature have an important impact on temperature inversion. De Mooij & Snellen (2009) also detected a possible slight offset in the timing of the secondary eclipse in K band. However, based on our 4 Spitzer channels, we place a 3sigma upper limit of |ecos(w)| < 0.0056 where e is the planets orbital eccentricity and w is the longitude of the periastron. This result strongly indicates that the orbit is circular, as expected from tidal circularization theory.

preprint2008arXiv

Design Considerations for a Ground-based Transit Search for Habitable Planets Orbiting M dwarfs

By targeting nearby M dwarfs, a transit search using modest equipment is capable of discovering planets as small as 2 Earth radii in the habitable zones of their host stars. The MEarth Project, a future transit search, aims to employ a network of ground-based robotic telescopes to monitor M dwarfs in the northern hemisphere with sufficient precision and cadence to detect such planets. Here we investigate the design requirements for the MEarth Project. We evaluate the optimal bandpass, and the necessary field of view, telescope aperture, and telescope time allocation on a star-by-star basis, as is possible for the well-characterized nearby M dwarfs. Through these considerations, 1,976 late M dwarfs (R < 0.33 Rsun) emerge as favorable targets for transit monitoring. Based on an observational cadence and on total telescope time allocation tailored to recover 90% of transit signals from planets in habitable zone orbits, we find that a network of ten 30 cm telescopes could survey these 1,976 M dwarfs in less than 3 years. A null result from this survey would set an upper limit (at 99% confidence) of 17% for the rate of occurrence of planets larger than 2 Earth radii in the habitable zones of late M dwarfs, and even stronger constraints for planets lying closer than the habitable zone. If the true occurrence rate of habitable planets is 10%, the expected yield would be 2.6 planets.

preprint2008arXiv

Detection of a Temperature Inversion in the Broadband Infrared Emission Spectrum of TrES-4

We estimate the strength of the bandpass-integrated thermal emission from the extrasolar planet TrES-4 at 3.6, 4.5, 5.8, and 8.0 micron using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We find relative eclipse depths of 0.137 +/- 0.011%, 0.148 +/- 0.016%, 0.261 +/- 0.059%, and 0.318 +/- 0.044% in these four bandpasses, respectively. We also place a 2 sigma upper limit of 0.37% on the depth of the secondary eclipse in the 16 micron IRS peak-up array. These eclipse depths reveal that TrES-4 has an emission spectrum similar to that of HD 209458b, which requires the presence of water emission bands created by an thermal inversion layer high in the atmosphere in order to explain the observed features. TrES-4 receives more radiation from its star than HD 209458b and has a correspondingly higher effective temperature, therefore the presence of a temperature inversion in this planet's atmosphere lends support to the idea that inversions might be correlated with the irradiance received by the planet. We find no evidence for any offset in the timing of the secondary eclipse, and place a 3 sigma upper limit of |ecos(omega)|<0.0058 where e is the planet's orbital eccentricity and omega is the argument of pericenter. From this we conclude that tidal heating from ongoing orbital circulatization is unlikely to be the explanation for TrES-4's inflated radius.

preprint2007arXiv

A map of the day-night contrast of the extrasolar planet HD 189733b

"Hot Jupiter" extrasolar planets are expected to be tidally locked because they are close (<0.05 astronomical units, where 1 AU is the average Sun-Earth distance) to their parent stars, resulting in permanent daysides and nightsides. By observing systems where the planet and star periodically eclipse each other, several groups have been able to estimate the temperatures of the daysides of these planets. A key question is whether the atmosphere is able to transport the energy incident upon the dayside to the nightside, which will determine the temperature at different points on the planet's surface. Here we report observations of HD 189733, the closest of these eclipsing planetary systems, over half an orbital period, from which we can construct a 'map' of the distribution of temperatures. We detected the increase in brightness as the dayside of the planet rotated into view. We estimate a minimum brightness temperature of 973 +/- 33 K and a maximum brightness temperature of 1212 +/- 11 K at a wavelength of 8 microns, indicating that energy from the irradiated dayside is efficiently redistributed throughout the atmosphere, in contrast to a recent claim for another hot Jupiter. Our data indicate that the peak hemisphere-integrated brightness occurs 16$\pm$6 degrees before opposition, corresponding to a hot spot shifted east of the substellar point. The secondary eclipse (when the planet moves behind the star) occurs 120 +/- 24 s later than predicted, which may indicate a slightly eccentric orbit.

preprint2007arXiv

The 3.6-8.0 Micron Broadband Emission Spectrum of HD 209458b: Evidence for an Atmospheric Temperature Inversion

We estimate the strength of the bandpass-integrated thermal emission from the extrasolar planet HD 209458b at 3.6, 4.5, 5.8, and 8.0 microns using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We observe a single secondary eclipse simultaneously in all four bandpasses and find relative eclipse depths of 0.00094 +/- 0.00009, 0.00213 +/- 0.00015, 0.00301 +/- 0.00043, and 0.00240 +/- 0.00026, respectively. These eclipse depths reveal that the shape of the inferred emission spectrum for the planet differs significantly from the predictions of standard atmosphere models; instead the most plausible explanation would require the presence of an inversion layer high in the atmosphere leading to significant water emission in the 4.5 and 5.8 micron bandpasses. This is the first clear indication of such a temperature inversion in the atmosphere of a hot Jupiter, as previous observations of other planets appeared to be in reasonably good agreement with the predictions of models without such an inversion layer.

preprint2007arXiv

TrES-3: A Nearby, Massive, Transiting Hot Jupiter in a 31-Hour Orbit

We describe the discovery of a massive transiting hot Jupiter with a very short orbital period (1.30619 d), which we name TrES-3. From spectroscopy of the host star GSC 03089-00929, we measure T_eff = 5720 +- 150 K, logg=4.6 +- 0.3, and vsini < 2 km/s, and derive a stellar mass of 0.90 +- 0.15 M_sun. We estimate a planetary mass of 1.92 +- 0.23 M_Jup, based on the sinusoidal variation of our high-precision radial velocity measurements. This variation has a period and phase consistent with our transit photometry. Our spectra show no evidence of line bisector variations that would indicate a blended eclipsing binary star. From detailed modeling of our B and z photometry of the 2.5%-deep transits, we determine a stellar radius 0.802 +- 0.046 R_sun and a planetary radius 1.295 +- 0.081 R_Jup. TrES-3 has one of the shortest orbital periods of the known transiting exoplanets, facilitating studies of orbital decay and mass loss due to evaporation, and making it an excellent target for future studies of infrared emission and reflected starlight.

preprint2006arXiv

Precise Radius Estimates for the Exoplanets WASP-1b and WASP-2b

We present precise z-band photometric time series spanning times of transit of the two exoplanets recently discovered by the SuperWASP collaboration. We find planetary radii of 1.44 +/- 0.08 R_J and 1.04 +/- 0.06 R_J for WASP-1b and WASP-2b, respectively. These error estimates include both random errors in the photometry and also the uncertainty in the stellar masses. Our results are 5 times more precise than the values derived from the discovery data alone. Our measurement of the radius of WASP-2b agrees with previously published models of hot Jupiters that include both a 20-M_Earth core of solid material and the effects of stellar insolation. In contrast, we find that the models cannot account for the large size of WASP-1b, even if the planet has no core. Thus, we add WASP-1b to the growing list of hot Jupiters that are larger than expected. This suggests that ``inflated'' hot Jupiters are more common than previously thought, and that any purported explanations involving highly unusual circumstances are disfavored.

preprint2005arXiv

A Method For Eclipsing Component Identification In Large Photometric Datasets

We describe an automated method for assigning the most likely physical parameters to the components of an eclipsing binary (EB), using only its photometric light curve and combined color. In traditional methods (e.g. WD and EBOP) one attempts to optimize a multi-parameter model over many iterations, so as to minimize the chi-squared value. We suggest an alternative method, where one selects pairs of coeval stars from a set of theoretical stellar models, and compares their simulated light curves and combined colors with the observations. This approach greatly reduces the EB parameter-space over which one needs to search, and allows one to determine the components' masses, radii and absolute magnitudes, without spectroscopic data. We have implemented this method in an automated program using published theoretical isochrones and limb-darkening coefficients. Since it is easy to automate, this method lends itself to systematic analyses of datasets consisting of photometric time series of large numbers of stars, such as those produced by OGLE, MACHO, TrES, HAT, and many others surveys.

preprint2005arXiv

Measurement of Spin-Orbit Alignment in an Extrasolar Planetary System

We determine the stellar, planetary, and orbital properties of the transiting planetary system HD 209458, through a joint analysis of high-precision radial velocities, photometry, and timing of the secondary eclipse. Of primary interest is the strong detection of the Rossiter-McLaughlin effect, the alteration of photospheric line profiles that occurs because the planet occults part of the rotating surface of the star. We develop a new technique for modeling this effect, and use it to determine the inclination of the planetary orbit relative to the apparent stellar equator [lambda = (-4.4 +/- 1.4) degrees], and the line-of-sight rotation speed of the star [v*sin(I) = (4.70 +/- 0.16) km/s]. The uncertainty in these quantities has been reduced by an order of magnitude relative to the pioneering measurements by Queloz and collaborators. The small but nonzero misalignment is probably a relic of the planet formation epoch, because the expected timescale for tidal coplanarization is larger than the age of the star. Our determination of v*sin(I) is a rare case in which rotational line broadening has been isolated from other broadening mechanisms.

preprint2001arXiv

Detection of an Extrasolar Planet Atmosphere

We report high precision spectrophotometric observations of four planetary transits of HD 209458, in the region of the sodium resonance doublet at 589.3 nm. We find that the photometric dimming during transit in a bandpass centered on the sodium feature is deeper by (2.32 +/- 0.57) x 10^{-4} relative to simultaneous observations of the transit in adjacent bands. We interpret this additional dimming as absorption from sodium in the planetary atmosphere, as recently predicted from several theoretical modeling efforts. Our model for a cloudless planetary atmosphere with a solar abundance of sodium in atomic form predicts more sodium absorption than we observe. There are several possibilities that may account for this reduced amplitude, including reaction of atomic sodium into molecular gases and/or condensates, photoionization of sodium by the stellar flux, a low primordial abundance of sodium, or the presence of clouds high in the atmosphere.