Source author record

Eric Gaidos

Eric Gaidos 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

53works
8topics
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

53 published item(s)

preprint2022arXiv

A close-in puffy Neptune with hidden friends: The enigma of TOI 620

We present the validation of a transiting low-density exoplanet orbiting the M2.5 dwarf TOI 620 discovered by the NASA TESS mission. We utilize photometric data from both TESS and ground-based follow-up observations to validate the ephemerides of the 5.09-day transiting signal and vet false positive scenarios. High-contrast imaging data are used to resolve the stellar host and exclude stellar companions at separations $\gtrsim 0.2''$. We obtain follow-up spectroscopy and corresponding precise radial velocities (RVs) with multiple PRV spectrographs to confirm the planetary nature of the transiting exoplanet. We calculate a 5$σ$ upper limit of $M_P < 7.1$ M$_\oplus$ and $ρ_P < 0.74$ g cm$^{-3}$, and we identify a non-transiting 17.7-day candidate. We also find evidence for a substellar (1-20 M$_{\rm J}$) companion with a projected separation $\lesssim 20$ au from a combined analysis of Gaia, AO imaging, and RVs. With the discovery of this outer companion, we carry out a detailed exploration of the possibilities that TOI 620 b might instead be a circum-secondary planet or a pair of eclipsing binary stars orbiting the host in a hierarchical triple system. We find, under scrutiny, that we can exclude both of these scenarios from the multi-wavelength transit photometry, thus validating TOI 620 b as a low-density exoplanet transiting the central star in this system. The low density of TOI 620 b makes it one of the most amenable exoplanets for atmospheric characterization, such as with JWST and Ariel, validated or confirmed by the TESS mission to date.

preprint2022arXiv

AU Microscopii in the FUV: Observations in Quiescence, During Flares, and Implications for AU Mic b and c

High energy X-ray and ultraviolet (UV) radiation from young stars impacts planetary atmospheric chemistry and mass loss. The active $\sim 22$ Myr M dwarf AU Mic hosts two exoplanets orbiting interior to its debris disk. Therefore, this system provides a unique opportunity to quantify the effects of stellar XUV irradiation on planetary atmospheres as a function of both age and orbital separation. In this paper we present over 5 hours of Far-UV (FUV) observations of AU Mic taken with the Cosmic Origins Spectrograph (COS; 1070-1360 Angstrom) on the Hubble Space Telescope (HST). We provide an itemization of $120$ emission features in the HST/COS FUV spectrum and quantify the flux contributions from formation temperatures ranging from $10^4-10^7$ K. We detect 13flares in the FUV white-light curve with energies ranging from $10^{29} - 10^{31}$ ergs. The majority of the energy in each of these flares is released from the transition region between the chromosphere and the corona. There is a 100$\times$ increase in flux at continuum wavelengths $λ< 1100$ Angstrom in each flare which may be caused by thermal Bremsstrahlung emission. We calculate that the baseline atmospheric mass-loss rate for AU Mic b is $\sim 10^8$ g s$^{-1}$, although this rate can be as high as $\sim 10^{14}$ g s$^{-1}$ during flares with $L_\textrm{flare} \simeq 10^{33}$ erg s$^{-1}$. Finally, we model the transmission spectra for AU Mic b and c with a new panchromatic spectrum of AU Mic c and motivate future JWST observations of these planets.

preprint2022arXiv

Color Dependence of the Transit Detectability for Young Active M-dwarfs

We investigate the planetary transit detectability in the presence of stellar rotational activity from light curves for young M-dwarfs and estimate improvements of the detection at near-infrared (NIR) wavelengths. Making maps of the transit signal detection efficiency over the orbital period and planetary radius with light curves of members of four clusters, Hyades, Praesepe, Pleiades, and Upper Scorpius observed by the K2 mission, we evaluate the detectability for the rotation period and modulation semi-amplitude. We find that the detection efficiency remarkably decreases to about 20% for rapidly rotators with P_{rot} <= 1 d and the lack of planets in Pleiades is likely due to the high fraction of rapidly rotating M-dwarfs. We also evaluate the improvements of the planet detection with NIR photometry via tests using mock light curves assuming that the signal amplitude of stellar rotation decreases at NIR wavelengths. Our results suggest that NIR photometric monitoring would double relative detection efficiency for transiting planetary candidates with P_{rot} <= 1 d and find planets around M-dwarfs with approximately 100 Myr missing in the past transit surveys from the space.

preprint2022arXiv

Is [Y/Mg] a Reliable Age Diagnostic for FGK Stars?

Current spectroscopic surveys are producing large catalogs of chemical abundances for stars of all types. The yttrium to magnesium ratio, [Y/Mg], has emerged as a candidate age indicator for solar twins in the local stellar neighborhood. However, it is unclear whether it is a viable age diagnostic for more diverse stellar types, so we investigate [Y/Mg] as an age indicator for the FGK-type planet host stars observed by $Kepler$. We find that the [Y/Mg] "Clock" is most precise for solar twins, with a [Y/Mg]/Age slope of $m$ = $-$0.0370 $\pm$ 0.0071 dex/Gyr and $σ_{\mathrm{Age}}$ = 2.6 Gyr. We attribute the lower precision compared to literature results to non-solar twins contaminating our solar twin sample and recommend a 1.5 Gyr systematic uncertainty for stellar ages derived with any [Y/Mg]-Age relation. We also analyzed the [Y/Mg] Clock as a function of $T_{\mathrm{eff}}$, $\log g$, and metallicity individually and find no strong trends, but compute statistically significant [Y/Mg]-Age relations for subsamples defined by ranges in $T_{\mathrm{eff}}$, $\log g$, and metallicity. Finally, we compare [Y/Mg] and rotation ages and find statistically similar trends as for isochrone ages, although we find that rotation ages perform better for GK dwarfs while isochrones perform better for FG subgiants. We conclude that the [Y/Mg] Clock is most precise for solar twins and analogs but is also a useful age diagnostic for FGK stars.

preprint2022arXiv

Transit Timing Variations for AU Microscopii b & c

We explore the transit timing variations (TTVs) of the young (22 Myr) nearby AU Mic planetary system. For AU Mic b, we introduce three Spitzer (4.5 $μ$m) transits, five TESS transits, 11 LCO transits, one PEST transit, one Brierfield transit, and two transit timing measurements from Rossiter-McLaughlin observations; for AU Mic c, we introduce three TESS transits. We present two independent TTV analyses. First, we use EXOFASTv2 to jointly model the Spitzer and ground-based transits and to obtain the midpoint transit times. We then construct an O--C diagram and model the TTVs with Exo-Striker. Second, we reproduce our results with an independent photodynamical analysis. We recover a TTV mass for AU Mic c of 10.8$^{+2.3}_{-2.2}$ M$_{E}$. We compare the TTV-derived constraints to a recent radial-velocity (RV) mass determination. We also observe excess TTVs that do not appear to be consistent with the dynamical interactions of b and c alone, and do not appear to be due to spots or flares. Thus, we present a hypothetical non-transiting "middle-d" candidate exoplanet that is consistent with the observed TTVs, the candidate RV signal, and would establish the AU Mic system as a compact resonant multi-planet chain in a 4:6:9 period commensurability. These results demonstrate that the AU Mic planetary system is dynamically interacting producing detectable TTVs, and the implied orbital dynamics may inform the formation mechanisms for this young system. We recommend future RV and TTV observations of AU Mic b and c to further constrain the masses and to confirm the existence of possible additional planet(s).

preprint2021arXiv

Flares, Rotation, and Planets of the AU Mic System from TESS Observations

AU Mic is a young ($\sim$24 Myr), pre-Main Sequence M~dwarf star that was observed in the first month of science observations of the Transiting Exoplanet Survey Satellite (TESS) and re-observed two years later. This target has photometric variability from a variety of sources that is readily apparent in the TESS light curves; spots induce modulation in the light curve, flares are present throughout (manifesting as sharp rises with slow exponential decay phases), and transits of AU Mic b may be seen by eye as dips in the light curve. We present a combined analysis of both TESS Sector 1 and Sector 27 AU Mic light curves including the new 20-second cadence data from TESS Year 3. We compare flare rates between both observations and analyze the spot evolution, showing that the activity levels increase slightly from Sector 1 to Sector 27. Furthermore, the 20-second data collection allows us to detect more flares, smaller flares, and better resolve flare morphology in white light as compared to the 2-minute data collection mode. We also refine the parameters for AU Mic b by fitting three additional transits of AU Mic b from Sector 27 using a model that includes stellar activity. We show that the transits exhibit clear transit timing variations (TTVs) with an amplitude of $\sim$80 seconds. We also detect three transits of a 2.8 $R_\oplus$ planet, AU Mic c, which has a period of 18.86 days.

preprint2020arXiv

Evidence for Spin-orbit Alignment in the TRAPPIST-1 System

In an effort to measure the Rossiter-McLaughlin effect for the TRAPPIST-1 system, we performed high-resolution spectroscopy during transits of planets e, f, and b. The spectra were obtained with the InfraRed Doppler spectrograph on the Subaru 8.2-m telescope, and were supplemented with simultaneous photometry obtained with a 1-m telescope of the Las Cumbres Observatory Global Telescope. By analyzing the anomalous radial velocities, we found the projected stellar obliquity to be $λ=1\pm 28$ degrees under the assumption that the three planets have coplanar orbits, although we caution that the radial-velocity data show correlated noise of unknown origin. We also sought evidence for the expected deformations of the stellar absorption lines, and thereby detected the "Doppler shadow" of planet b with a false alarm probability of $1.7\,\%$. The joint analysis of the observed residual cross-correlation map including the three transits gave $λ=19_{-15}^{+13}$ degrees. These results indicate that the the TRAPPIST-1 star is not strongly misaligned with the common orbital plane of the planets, although further observations are encouraged to verify this conclusion.

preprint2020arXiv

Limits on the Spin-Orbit Angle and Atmospheric Escape for the 22 Myr-old Planet AU Mic b

We obtained spectra of the pre-main sequence star AU Microscopii during a transit of its Neptune-sized planet to investigate its orbit and atmosphere. We used the high-dispersion near-infrared spectrograph IRD on the Subaru telescope to detect the Doppler "shadow" from the planet and constrain the projected stellar obliquity. Modeling of the observed planetary Doppler shadow suggests a spin-orbit alignment of the system ($λ=-4.7_{-6.4}^{+6.8}$ degrees), but additional observations are needed to confirm this finding. We use both the IRD data and spectra obtained with NIRSPEC on Keck-II to search for absorption in the 1083 nm line of metastable triplet He I by the planet's atmosphere and place an upper limit for the equivalent width of 3.7 mÅat 99 $\%$ confidence. With this limit and a Parker wind model we constrain the escape rate from the atmosphere to $<0.15-0.45\, M_{\oplus}$ Gyr$^{-1}$, comparable to the rates predicted by an XUV energy-limited escape calculation and hydrodynamic models, but refinement of the planet mass is needed for rigorous tests.

preprint2020arXiv

Precision radial velocity measurements by the forward-modeling technique in the near-infrared

Precision radial velocity (RV) measurements in the near-infrared are a powerful tool to detect and characterize exoplanets around low-mass stars or young stars with higher magnetic activity. However, the presence of strong telluric absorption lines and emission lines in the near infrared that significantly vary in time can prevent extraction of RV information from these spectra by classical techniques, which ignore or mask the telluric lines. We present a methodology and pipeline to derive precision RVs from near-infrared spectra using a forward-modeling technique. We applied this to spectra with a wide wavelength coverage (Y, J, and H bands, simultaneously), taken by the InfraRed Doppler (IRD) spectrograph on the Subaru 8.2-m telescope. Our pipeline extracts the instantaneous instrumental profile of the spectrograph for each spectral segment, based on a reference spectrum of the laser-frequency comb that is injected into the spectrograph simultaneously with the stellar light. These profiles are used to derive the intrinsic stellar template spectrum, which is free from instrumental broadening and telluric features, as well as model and fit individual observed spectra in the RV analysis. Implementing a series of numerical simulations using theoretical spectra that mimic IRD data, we test the pipeline and show that IRD can achieve $<2$ m s$^{-1}$ precision for slowly rotating mid-to-late M dwarfs with a signal-to-noise ratio $> 100$ per pixel at 1000 nm. Dependences of RV precision on various stellar parameters (e.g., $T_{\rm eff}$, $v\sin i$, [Fe/H]) and the impact of telluric-line blendings on the RV accuracy are discussed through the mock spectra analyses. We also apply the RV-analysis pipeline to the observed spectra of GJ 699 and TRAPPIST-1, demonstrating that the spectrograph and the pipeline are capable of an RV accuracy of $<3$ m s$^{-1}$ at least on a time scale of a few months.

preprint2020arXiv

The Gaia-Kepler Stellar Properties Catalog. I. Homogeneous Fundamental Properties for 186,301 Kepler Stars

An accurate and precise Kepler Stellar Properties Catalog is essential for the interpretation of the Kepler exoplanet survey results. Previous Kepler Stellar Properties Catalogs have focused on reporting the best-available parameters for each star, but this has required combining data from a variety of heterogeneous sources. We present the Gaia-Kepler Stellar Properties Catalog, a set of stellar properties of 186,301 Kepler stars, homogeneously derived from isochrones and broadband photometry, Gaia Data Release 2 parallaxes, and spectroscopic metallicities, where available. Our photometric effective temperatures, derived from $g-K_s$ colors, are calibrated on stars with interferometric angular diameters. Median catalog uncertainties are 112 K for $T_{\mathrm{eff}}$, 0.05 dex for $\log g$, 4% for $R_\star$, 7% for $M_\star$, 13% for $ρ_\star$, 10% for $L_\star$, and 56% for stellar age. These precise constraints on stellar properties for this sample of stars will allow unprecedented investigations into trends in stellar and exoplanet properties as a function of stellar mass and age. In addition, our homogeneous parameter determinations will permit more accurate calculations of planet occurrence and trends with stellar properties.

preprint2020arXiv

The Gaia-Kepler Stellar Properties Catalog. II. Planet Radius Demographics as a Function of Stellar Mass and Age

Studies of exoplanet demographics require large samples and precise constraints on exoplanet host stars. Using the homogeneous Kepler stellar properties derived using Gaia Data Release 2 by Berger et al. (2020), we re-compute Kepler planet radii and incident fluxes and investigate their distributions with stellar mass and age. We measure the stellar mass dependence of the planet radius valley to be $d \log R_{\mathrm{p}}$/$d \log M_\star = 0.26^{+0.21}_{-0.16}$, consistent with the slope predicted by a planet mass dependence on stellar mass ($0.24-0.35$) and core-powered mass-loss (0.33). We also find first evidence of a stellar age dependence of the planet populations straddling the radius valley. Specifically, we determine that the fraction of super-Earths ($1-1.8 \mathrm{R_\oplus}$) to sub-Neptunes ($1.8-3.5 \mathrm{R_\oplus}$) increases from $0.61 \pm 0.09$ at young ages (< 1 Gyr) to $1.00 \pm 0.10$ at old ages (> 1 Gyr), consistent with the prediction by core-powered mass-loss that the mechanism shaping the radius valley operates over Gyr timescales. Additionally, we find a tentative decrease in the radii of relatively cool ($F_{\mathrm{p}} < 150 \mathrm{F_\oplus}$) sub-Neptunes over Gyr timescales, which suggests that these planets may possess H/He envelopes instead of higher mean molecular weight atmospheres. We confirm the existence of planets within the hot sub-Neptunian "desert" ($2.2 < R_{\mathrm{p}} < 3.8 \mathrm{R_\oplus}$, $F_{\mathrm{p}} > 650 \mathrm{F_\oplus}$) and show that these planets are preferentially orbiting more evolved stars compared to other planets at similar incident fluxes. In addition, we identify candidates for cool ($F_{\mathrm{p}} < 20 \mathrm{F_\oplus}$) inflated Jupiters, present a revised list of habitable zone candidates, and find that the ages of single- and multiple-transiting planet systems are statistically indistinguishable.

preprint2020arXiv

Very regular high-frequency pulsation modes in young intermediate-mass stars

Asteroseismology is a powerful tool for probing the internal structures of stars by using their natural pulsation frequencies. It relies on identifying sequences of pulsation modes that can be compared with theoretical models, which has been done successfully for many classes of pulsators, including low-mass solar-type stars, red giants, high-mass stars and white dwarfs. However, a large group of pulsating stars of intermediate mass--the so-called delta Scuti stars--have rich pulsation spectra for which systematic mode identification has not hitherto been possible. This arises because only a seemingly random subset of possible modes are excited, and because rapid rotation tends to spoil the regular patterns. Here we report the detection of remarkably regular sequences of high-frequency pulsation modes in 60 intermediate-mass main-sequence stars, allowing definitive mode identification. Some of these stars have space motions that indicate they are members of known associations of young stars, and modelling of their pulsation spectra confirms that these stars are indeed young.

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

Giant planet occurrence within 0.2 AU of low-luminosity red giant branch stars with K2

Every Sun-like star will eventually evolve into a red giant, a transition which can profoundly affect the evolution of a surrounding planetary system. The timescale of dynamical planet evolution and orbital decay has important implications for planetary habitability, as well as post-main sequence star and planet interaction, evolution and internal structure. Here, we investigate these effects by estimating planet occurrence around 2476 low-luminosity red giant branch (LLRGB) stars observed by the NASA K2 mission. We measure stellar masses and radii using asteroseismology, with median random uncertainties of 3.7% in mass and 2.2% in radius. We compare this planet population to the known population of planets around dwarf Sun-like stars, accounting for detection efficiency differences between the stellar populations. We find that 0.51% +/- 0.29% of LLRGB stars host planets larger than Jupiter with orbital periods less than 10 days, tentatively higher than main sequence stars hosting similar planets (0.15% +/- 0.06%). Our results suggest that the effects of stellar evolution on the occurrence of close-in planets larger than Jupiter is not significant until stars have begun ascending substantially up the red giant branch (>~ 5-6 Rsun).

preprint2016arXiv

Atmosphere-interior exchange on hot rocky exoplanets

We provide estimates of atmospheric pressure and surface composition on short-period rocky exoplanets with dayside magma pools and silicate vapor atmospheres. Atmospheric pressure tends toward vapor-pressure equilibrium with surface magma, and magma-surface composition is set by the competing effects of fractional vaporization and surface-interior exchange. We use basic models to show how surface-interior exchange is controlled by the planet's temperature, mass, and initial composition. We assume that mantle rock undergoes bulk melting to form the magma pool, and that winds flow radially away from the substellar point. With these assumptions, we find that: (1) atmosphere-interior exchange is fast when the planet's bulk-silicate FeO concentration is low, and slow when FeO concentration is high; (2) magma pools are compositionally well-mixed for substellar temperatures $\lesssim$ 2400 K, but compositionally variegated and rapidly variable for substellar temperatures $\gtrsim$ 2400 K; (3) currents within the magma pool tend to cool the top of the solid mantle ("tectonic refrigeration"); (4) contrary to earlier work, many magma planets have time-variable surface compositions.

preprint2016arXiv

The K2-ESPRINT Project III: A Close-in Super-Earth around a Metal-rich Mid-M Dwarf

We validate a $R_p=2.32\pm 0.24R_\oplus$ planet on a close-in orbit ($P=2.260455\pm 0.000041$ days) around K2-28 (EPIC 206318379), a metal-rich M4-type dwarf in the Campaign 3 field of the K2 mission. Our follow-up observations included multi-band transit observations from the optical to the near infrared, low-resolution spectroscopy, and high-resolution adaptive-optics (AO) imaging. We perform a global fit to all the observed transits using a Gaussian process-based method and show that the transit depths in all passbands adopted for the ground-based transit follow-ups ($r'_2, z_\mathrm{s,2}, J, H, K_\mathrm{s}$) are within $\sim 2σ$ of the K2 value. Based on a model of the background stellar population and the absence of nearby sources in our AO imaging, we estimate the probability that a background eclipsing binary could cause a false positive to be $< 2\times 10^{-5}$. We also show that K2-28 cannot have a physically associated companion of stellar type later than M4, based on the measurement of almost identical transit depths in multiple passbands. There is a low probability for a M4 dwarf companion ($\approx 0.072_{-0.04}^{+0.02}$), but even if this were the case, the size of K2-28b falls within the planetary regime. K2-28b has the same radius (within $1σ$) and experiences a similar irradiation from its host star as the well-studied GJ~1214b. Given the relative brightness of K2-28 in the near infrared ($m_\mathrm{Kep}=14.85$ mag and $m_H=11.03$ mag) and relatively deep transit ($0.6-0.7\%$), a comparison between the atmospheric properties of these two planets with future observations would be especially interesting.

preprint2016arXiv

Zodiacal Exoplanets In Time (ZEIT) I: A Neptune-sized planet orbiting an M4.5 dwarf in the Hyades Star Cluster

Studying the properties of young planetary systems can shed light on how the dynamics and structure of planets evolve during their most formative years. Recent K2 observations of nearby young clusters (10-800 Myr) have enabled the discovery of such planetary systems. Here we report the discovery of a Neptune-sized planet transiting an M4.5 dwarf (K2-25) in the Hyades cluster (650-800 Myr). The lightcurve shows a strong periodic signal at 1.88 days, which we attribute to spot coverage and rotation. We confirm the planet host is a member of the Hyades by measuring the radial velocity of the system with the high-resolution near-infrared spectrograph IGRINS. This enables us to calculate a distance based on EPIC 210490365's kinematics and membership to the Hyades, which in turn provides a stellar radius and mass to 5-10%, better than what is currently possible for most Kepler M dwarfs (12-20%). We use the derived stellar density as a prior on fitting the K2 transit photometry, which provides weak constraints on eccentricity. Utilizing a combination of adaptive optics imaging and high-resolution spectra we rule out the possibility that the signal is due to a bound or background eclipsing binary, confirming the transits' planetary origin. EPIC 210490365b has a radius ($3.43^{+0.95}_{-0.31}$R$_{E}$) much larger than older Kepler planets with similar orbital periods (3.484 days) and host-star masses (0.29$M_{\odot}$). This suggests that close-in planets lose some of their atmospheres past the first few hundred Myr. Additional transiting planets around the Hyades, Pleiades, and Praesepe clusters from K2 will help confirm if this planet is atypical or representative of other close-in planets of similar age.

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

How to Constrain Your M Dwarf: measuring effective temperature, bolometric luminosity, mass, and radius

Precise and accurate parameters for late-type (late K and M) dwarf stars are important for characterization of any orbiting planets, but such determinations have been hampered by these stars' complex spectra and dissimilarity to the Sun. We exploit an empirically calibrated method to estimate spectroscopic effective temperature ($T_{\rm{eff}}$) and the Stefan-Boltzmann law to determine radii of 183 nearby K7-M7 single stars with a precision of 2-5%. Our improved stellar parameters enable us to develop model-independent relations between $T_{\rm{eff}}$ or absolute magnitude and radius, as well as between color and $T_{\rm{eff}}$. The derived $T_{\rm{eff}}$-radius relation depends strongly on [Fe/H], as predicted by theory. The relation between absolute $K_S$ magnitude and radius can predict radii accurate to $\simeq$3%. We derive bolometric corrections to the $VR_CI_CgrizJHK_S$ and Gaia passbands as a function of color, accurate to 1-3%. We confront the reliability of predictions from Dartmouth stellar evolution models using a Markov Chain Monte Carlo to find the values of unobservable model parameters (mass, age) that best reproduce the observed effective temperature and bolometric flux while satisfying constraints on distance and metallicity as Bayesian priors. With the inferred masses we derive a semi-empirical mass-absolute magnitude relation with a scatter of 2% in mass. The best-agreement models over-predict stellar $T_{\rm{eff}}$s by an average of 2.2% and under-predict stellar radii by 4.6%, similar to differences with values from low-mass eclipsing binaries. These differences are not correlated with metallicity, mass, or indicators of activity, suggesting issues with the underlying model assumptions e.g., opacities or convective mixing length.

preprint2015arXiv

Kepler-445, Kepler-446 and the Occurrence of Compact Multiples Orbiting Mid-M Dwarf Stars

We confirm and characterize the exoplanetary systems Kepler-445 and Kepler-446: two mid-M dwarf stars, each with multiple, small, short-period transiting planets. Kepler-445 is a metal-rich ([Fe/H]=+0.25 $\pm$ 0.10) M4 dwarf with three transiting planets, and Kepler-446 is a metal-poor ([Fe/H]=-0.30 $\pm$ 0.10) M4 dwarf also with three transiting planets. Kepler-445c is similar to GJ 1214b: both in planetary radius and the properties of the host star. The Kepler-446 system is similar to the Kepler-42 system: both are metal-poor with large galactic space velocities and three short-period, likely-rocky transiting planets that were initially assigned erroneously large planet-to-star radius ratios. We independently determined stellar parameters from spectroscopy and searched for and fitted the transit light curves for the planets, imposing a strict prior on stellar density in order to remove correlations between the fitted impact parameter and planet-to-star radius ratio for short-duration transits. Combining Kepler-445, Kepler-446 and Kepler-42, and isolating all mid-M dwarf stars observed by Kepler with the precision necessary to detect similar systems, we calculate that 21 $^{+7}_{-5}$ % of mid-M dwarf stars host compact multiples (multiple planets with periods of less than 10 days) for a wide range of metallicities. We suggest that the inferred planet masses for these systems support highly efficient accretion of protoplanetary disk metals by mid-M dwarf protoplanets.

preprint2015arXiv

The Enigmatic and Ephemeral M Dwarf System KOI 6705: Cheshire Cat or Wild Goose?

We confirm a 0.995 d periodic planetary transit-like signal, KOI 6705.01, in the Kepler lightcurve of the star KIC 6423922. Optical and infrared spectra show that this star is a mid M-type dwarf with an effective temperature $= 3327 \pm 60$K, metallicity [Fe/H] $= -0.08 \pm 0.10$, radius $= 0.31 \pm 0.03 R_{\odot}$, and mass $= 0.28 \pm 0.05 M_{\odot}$. The star is $\approx 70$ pc away and its space motion, rotation period, and lack of H$α$ emission indicate it is an older member of the "thin disk" population. On the other hand, the star exhibits excess infrared emission suggesting a dust disk more typical of a very young star. If the KOI 6705.01 signal is produced by a planet, the transit depth of 60 ppm means its radius is only $0.26^{+0.034}_{-0.029} R_{\oplus}$, or about the size of the Moon. However, the duration ($\gtrsim 3$~hr) and time variation of KOI 6705.01 are anomalous: the signal was undetected in the first two years of the mission and increased through the latter two years. These characteristics require implausible orbits and material properties for any planet and rule out such an explanation, although a dust cloud is possible. We excluded several false positive scenarios including background stars, scattered light from stars that are nearby on the sky, and electronic cross-talk between detector readout channels. We find the most likely explanation to be that KOI 6705.01 is a false positive created by charge transfer inefficiency in a detector column on which KIC 6423922 and a 1.99 d eclipsing binary both happened to fall.

preprint2015arXiv

What are little worlds made of? Stellar abundances and the building blocks of planets

If the photospheres of solar-type stars represent the composition of circumstellar disks from which any planets formed, spectroscopic determinations of stellar elemental abundances offer information on the composition of those planets, including smaller, rocky planets. In particular, the C/O ratio is proposed to be a key determinant of the composition of solids that condense from disk gas and are incorporated into planets. Also, planets may leave chemical signatures on the photospheres of their host stars by sequestering heavy elements, or by being accreted by the stars. The presence, absence, and composition of planets could be revealed by small differences in the relative abundances between stars. I critically examine these scenarios and show that (i) a model of Galactic chemical evolution predicts that the C/O ratio is expected to be close to the solar value and vary little between dwarf stars in the solar neighborhood; (ii) spectroscopic surveys of M dwarf stars limit the occurrence of stars with C/O $\gtrsim 1$ to $<10^{-3}$; and (iii) planetesimal chemistry will be controlled by the composition of oxygen-rich dust inherited from the molecular cloud and processed in a dust-rich environment, not a gas with the stellar composition. A second generation of more reduced planetesimals could be produced by re-equilibration of some material with dust-depleted gas. Finally, I discuss how minor differences in relative abundances between stars that correlate with condensation temperature can be explained by dust-gas segregation, perhaps in circumstellar disks, rather than planet formation.

preprint2014arXiv

A Statistical Reconstruction of the Planet Population Around Kepler Solar-Type Stars

Using the cumulative catalog of planets detected by the NASA Kepler mission, we reconstruct the intrinsic occurrence of Earth- to Neptune-size (1 - 4$R_{\oplus}$) planets and their distributions with radius and orbital period. We analyze 76,711 solar-type ($0.8<R_*/R_{\odot}<1.2 $) stars with 430 planets on 20-200~d orbits, excluding close-in planets that may have been affected by the proximity to the host star. Our analysis considers errors in planet radii and includes an "iterative simulation" technique that does not bin the data. We find a radius distribution that peaks at 2-2.8 Earth radii, with lower numbers of smaller and larger planets. These planets are uniformly distributed with logarithmic period, and the mean number of such planets per star is $0.46 \pm 0.03$. The occurrence is $\sim 0.66$ if planets interior to 20~d are included. We estimate the occurrence of Earth-size planets in the "habitable zone" (defined as 1-2$R_{\oplus}$, 0.99-1.7 AU for solar-like stars) as $6.4^{+3.4}_{-1.1} \%$. Our results largely agree with those of Petigura et al. (2013), although we find a higher occurrence of 2.8-4 Earth-radii planets. The reasons for this excess are the inclusion of errors in planet radius, updated Huber et al. (2014) stellar parameters, and also the exclusion of planets which may have been affected by proximity to the host star.

preprint2014arXiv

M Dwarf Metallicities and Giant Planet Occurrence: Ironing Out Uncertainties and Systematics

Comparisons between the planet populations around solar-type stars and those orbiting M dwarfs shed light on the possible dependence of planet formation and evolution on stellar mass. However, such analyses must control for other factors, i.e. metallicity, a stellar parameter which strongly influences the occurrence of gas giant planets. We obtained infrared spectra of 121 M dwarfs stars monitored by the California Planet Search (CPS) and determined metallicities with an accuracy of 0.08 dex. The mean and standard deviation of the sample is -0.05 and 0.20 dex, respectively. We parameterized the metallicity dependence of the occurrence of giant planets on orbits with period less than 2 yr around solar-type stars and applied this to our M dwarf sample to estimate the expected number of giant planets. The number of detected planets (3) is lower than the predicted number (6.4) but the difference is not very significant (12% probability of finding as many or fewer planets). The three M dwarf planet hosts are not especially metal rich and the most likely value of the power-law index relating planet occurrence to metallicity is 1.06 dex per dex for M dwarfs compared to 1.80 for solar-type stars; this difference, however, is comparable to uncertainties. Giant planet occurrence around both types of stars allows, but does not necessarily require, mass dependence of $\sim 1$ dex per dex. The actual planet-mass-metallicity relation may be complex and elucidating it will require larger surveys like those to be conducted by ground-based infrared spectrographs and the Gaia space astrometry mission.

preprint2014arXiv

Prospecting in late-type dwarfs: a calibration of infrared and visible spectroscopic metallicities of late-K and M dwarfs spanning 1.5 dex

Knowledge of late K and M dwarf metallicities can be used to guide planet searches and constrain planet formation models. However, the determination of metallicities of late-type stars is difficult because visible wavelength spectra of their cool atmospheres contain many overlapping absorption lines, preventing the measurement of equivalent widths. We present new methods, and improved calibrations of existing methods, to determine metallicities of late-K and M dwarfs from moderate resolution (1300 < R < 2000) visible and infrared spectra. We select a sample of 112 wide binary systems that contain a late-type companion to a solar-type primary star. Our sample includes 62 primary stars with previously published metallicities, as well as 50 stars with metallicities determined from our own observations. We use our sample to empirically determine which features in the spectrum of the companion are best correlated with the metallicity of the primary. We derive metallicity calibrations for different wavelength ranges, and show that it is possible to get metallicities reliable to < 0.10 dex using either visible, J, H, or K band spectra. Our calibrations are applicable to dwarfs with metallicities of -1.04 < [Fe/H]< +0.56 and spectral types from K7 to M5. Lastly, we use our sample of wide binaries to test and refine existing calibrations to determine M dwarf metallicities. We find that the zeta parameter, which measures the ratio of TiO can CaH bands, is correlated with [Fe/H] for super-solar metallicities, and zeta does not always correctly identify metal-poor M dwarfs. We also find that existing calibrations in the K and H band are quite reliable for stars with [Fe/H] > -0.5, but are less useful for more metal-poor stars.

preprint2014arXiv

Prospecting in Ultracool Dwarfs: Measuring the Metallicities of Mid- and Late-M Dwarfs

Metallicity is a fundamental parameter that contributes to the physical characteristics of a star. However, the low temperatures and complex molecules present in M dwarf atmospheres make it difficult to measure their metallicities using techniques that have been commonly used for Sun-like stars. Although there has been significant progress in developing empirical methods to measure M dwarf metallicities over the last few years, these techniques have been developed primarily for early- to mid-M dwarfs. We present a method to measure the metallicity of mid- to late-M dwarfs from moderate resolution (R~2000) K-band (2.2 microns) spectra. We calibrate our formula using 44 wide binaries containing an F, G, K, or early M primary of known metallicity and a mid- to late-M dwarf companion. We show that similar features and techniques used for early M dwarfs are still effective for late-M dwarfs. Our revised calibration is accurate to 0.07 dex for M4.5-M9.5 dwarfs with -0.58<[Fe/H]<+0.56 and shows no systematic trends with spectral type, metallicity, or the method used to determine the primary star metallicity. We show that our method gives consistent metallicities for the components of M+M wide binaries. We verify that our new formula works for unresolved binaries by combining spectra of single stars. Lastly, we show that our calibration gives consistent metallicities with the Mann et al. (2013a) study for overlapping (M4-M5) stars, establishing that the two calibrations can be used in combination to determine metallicities across the entire M dwarf sequence.

preprint2014arXiv

Revised Stellar Properties of Kepler Targets for the Quarter 1-16 Transit Detection Run

We present revised properties for 196,468 stars observed by the NASA Kepler Mission and used in the analysis of Quarter 1-16 (Q1-Q16) data to detect and characterize transiting exoplanets. The catalog is based on a compilation of literature values for atmospheric properties (temperature, surface gravity, and metallicity) derived from different observational techniques (photometry, spectroscopy, asteroseismology, and exoplanet transits), which were then homogeneously fitted to a grid of Dartmouth stellar isochrones. We use broadband photometry and asteroseismology to characterize 11,532 Kepler targets which were previously unclassified in the Kepler Input Catalog (KIC). We report the detection of oscillations in 2,762 of these targets, classifying them as giant stars and increasing the number of known oscillating giant stars observed by Kepler by ~20% to a total of ~15,500 stars. Typical uncertainties in derived radii and masses are ~40% and ~20%, respectively, for stars with photometric constraints only, and 5-15% and ~10% for stars based on spectroscopy and/or asteroseismology, although these uncertainties vary strongly with spectral type and luminosity class. A comparison with the Q1-Q12 catalog shows a systematic decrease in radii for M dwarfs, while radii for K dwarfs decrease or increase depending on the Q1-Q12 provenance (KIC or Yonsei-Yale isochrones). Radii of F-G dwarfs are on average unchanged, with the exception of newly identified giants. The Q1-Q16 star properties catalog is a first step towards an improved characterization of all Kepler targets to support planet occurrence studies.

preprint2014arXiv

Spectro-Thermometry of M dwarfs and their candidate planets: too hot, too cool, or just right?

We use moderate-resolution spectra of nearby late K and M dwarf stars with parallaxes and interferometrically determined radii to refine their effective temperatures, luminosities, and metallicities. We use these revised values to calibrate spectroscopic techniques to infer the fundamental parameters of more distant late-type dwarf stars. We demonstrate that, after masking out poorly modeled regions, the newest version of the PHOENIX atmosphere models accurately reproduce temperatures derived bolometrically. We apply methods to late-type hosts of transiting planet candidates in the Kepler field, and calculate effective temperature, radius, mass, and luminosity with typical errors of 57 K, 7%, 11%, and 13%, respectively. We find systematic offsets between our values and those from previous analyses of the same stars, which we attribute to differences in atmospheric models utilized for each study. We investigate which of the planets in this sample are likely to orbit in the circumstellar habitable zone. We determine that four candidate planets (KOI 854.01, 1298.02, 1686.01, and 2992.01) are inside of or within 1-sigma of a conservative definition of the habitable zone, but that several planets identified by previous analyses are not (e.g. KOI 1422.02 and KOI 2626.01). Only one of the four habitable-zone planets is Earth sized, suggesting a downward revision in the occurrence of such planets around M dwarfs. These findings highlight the importance of measuring accurate stellar parameters when deriving parameters of their orbiting planets.

preprint2014arXiv

The Near-Ultraviolet Luminosity Function of Young, Early M-Type Dwarf Stars

Planets orbiting within the close-in habitable zones of M dwarf stars will be exposed to elevated high-energy radiation driven by strong magneto-hydrodynamic dynamos during stellar youth. Near-ultraviolet (NUV) irradiation can erode and alter the chemistry of planetary atmospheres, and a quantitative description of the evolution of NUV emission from M dwarfs is needed when modeling these effects. We investigated the NUV luminosity evolution of early M-type dwarfs by cross-correlating the Lepine & Gaidos (2011) catalog of bright M dwarfs with the GALEX catalog of NUV (1771-2831A) sources. Of the 4805 sources with GALEX counterparts, 797 have NUV emission significantly (> 2.5 sigma) in excess of an empirical basal level. We inspected these candidate active stars using visible-wavelength spectra, high-resolution adaptive optics imaging, time-series photometry, and literature searches to identify cases where the elevated NUV emission is due to unresolved background sources or stellar companions; we estimated the overall occurrence of these "false positives" as ~16%. We constructed a NUV luminosity function that accounted for false positives, detection biases of the source catalogs, and GALEX upper limits. We found the NUV luminosity function to be inconsistent with predictions from a constant star-formation rate and simplified age-activity relation defined by a two-parameter power law.

preprint2014arXiv

Touchstone Stars: Highlights from the Cool Stars 18 Splinter Session

We present a summary of the splinter session on "touchstone stars" -- stars with directly measured parameters -- that was organized as part of the Cool Stars 18 conference. We discuss several methods to precisely determine cool star properties such as masses and radii from eclipsing binaries, and radii and effective temperatures from interferometry. We highlight recent results in identifying and measuring parameters for touchstone stars, and ongoing efforts to use touchstone stars to determine parameters for other stars. We conclude by comparing the results of touchstone stars with cool star models, noting some unusual patterns in the differences.

preprint2013arXiv

A Spectroscopic Catalog of the Brightest (J<9) M Dwarfs in the Northern Sky

We present a spectroscopic catalog of the 1,564 brightest (J<9) M dwarf candidates in the northern sky, as selected from the SUPERBLINK proper motion catalog. Observations confirm 1,408 of the candidates to be late-K and M dwarfs with spectral subtypes K7-M6. From the low μ>40 mas/yr proper motion limit and high level of completeness of the SUPERBLINK catalog in that magnitude range, we estimate that our spectroscopic census most likely includes >90% of all existing, northern-sky M dwarfs with apparent magnitude J<9. Only 682 stars in our sample are listed in the Third Catalog of Nearby Stars (CNS3); most others are relative unknowns and have spectroscopic data presented here for the first time. Spectral subtypes are assigned based on spectral index measurements of CaH and TiO molecular bands. A comparison of spectra from the same stars obtained at different observatories however reveals that spectral band index measurements are dependent on spectral resolution, spectrophotometric calibration, and other instrumental factors. After systematic corrections and a recalibration of the subtype-index relationships, we find that we can consistently and reliably classify all our stars to a half-subtype precision. The use of corrected spectral indices further requires us to recalibrate the ζparameter, from which we estimate metallicities +/-0.5dex accuracy in dwarfs of subtypes M2 and earlier, and +/-0.2dex accuracy for later (M3-M5) subtypes. Existing geometric parallax measurements are extracted from the literature for 624 stars, and are used to determine spectroscopic and photometric distances for all the other stars. Our catalog will be most useful to guide the selection of the best M dwarf targets for exoplanet searches, in particular those using high-precision radial velocity measurements.

preprint2013arXiv

An Understanding of the Shoulder of Giants: Jovian Planets around Late K Dwarf Stars and the Trend with Stellar Mass

Analyses of exoplanet statistics suggest a trend of giant planet occurrence with host star mass, a clue to how planets like Jupiter form. One missing piece of the puzzle is the occurrence around late K dwarf stars (masses of 0.5-0.75Msun and effective temperatures of 3900-4800K). We analyzed four years of Doppler radial velocities data of 110 late K dwarfs, one of which hosts two previously reported giant planets. We estimate that 4.0+/-2.3% of these stars have Saturn-mass or larger planets with orbital periods <245d, depending on the planet mass distribution and RV variability of stars without giant planets. We also estimate that 0.7+/-0.5% of similar stars observed by Kepler have giant planets. This Kepler rate is significantly (99% confidence) lower than that derived from our Doppler survey, but the difference vanishes if only the single Doppler system (HIP 57274) with completely resolved orbits is considered. The difference could also be explained by the exclusion of close binaries (without giant planets) from the Doppler but not Kepler surveys, the effect of long-period companions and stellar noise on the Doppler data, or an intrinsic difference between the two populations. Our estimates for late K dwarfs bridge those for solar-type stars and M dwarfs and support a positive trend with stellar mass. Small sample size precludes statements about finer structure, e.g. a "shoulder" in the distribution of giant planets with stellar mass. Future surveys such as the Next Generation Transit Survey and the Transiting Exoplanet Satellite Survey will ameliorate this deficiency.

preprint2013arXiv

Candidate Planets in the Habitable Zones of Kepler Stars

A key goal of the Kepler mission is the discovery of Earth-size transiting planets in "habitable zones" where stellar irradiance maintains a temperate climate on an Earth-like planet. Robust estimates of planet radius and irradiance require accurate stellar parameters, but most Kepler systems are faint, making spectroscopy difficult and prioritization of targets desirable. The parameters of 2035 host stars were estimated by Bayesian analysis and the probabilities p(HZ) that 2738 candidate or confirmed planets orbit in the habitable zone were calculated. Dartmouth Stellar Evolution Program models were compared to photometry from the Kepler Input Catalog, priors for stellar mass, age, metallicity and distance, and planet transit duration. The analysis yielded probability density functions for calculating confidence intervals of planet radius and stellar irradiance, as well as p(HZ). Sixty-two planets have p(HZ) > 0.5 and a most probable stellar irradiance within habitable zone limits. Fourteen of these have radii less than twice the Earth; the objects most resembling Earth in terms of radius and irradiance are KOIs 2626.01 and 3010.01, which orbit late K/M-type dwarf stars. The fraction of Kepler dwarf stars with Earth-size planets in the habitable zone (eta-Earth) is 0.46, with a 95% confidence interval of 0.31-0.64. Parallaxes from the Gaia mission will reduce uncertainties by more than a factor of five and permit definitive assignments of transiting planets to the habitable zones of Kepler stars.

preprint2013arXiv

Narrow-K-Band Observations of the GJ 1214 System

GJ 1214 is a nearby M dwarf star that hosts a transiting super-Earth-size planet, making this system an excellent target for atmospheric studies. Most studies find that the transmission spectrum of GJ 1214b is flat, which favors either a high mean molecular weight or cloudy/hazy hydrogen (H) rich atmosphere model. Photometry at short wavelengths (< 0.7 micron) and in the K-band can discriminate the most between these different atmosphere models for GJ 1214b, but current observations do not have sufficiently high precision. We present photometry of seven transits of GJ 1214b through a narrow K-band (2.141 micron) filter with the Wide Field Camera on the 3.8 m United Kingdom Infrared Telescope. Our photometric precision is typically 1.7x10^-3 (for a single transit), comparable with other ground-based observations of GJ 1214b. We measure a planet-star radius ratio of 0.1158+/-0.0013, which, along with other studies, also supports a flat transmission spectrum for GJ 1214b. Since this does not exclude a scenario where GJ 1214b has a H-rich envelope with heavy elements that are sequestered below a cloud/haze layer, we compare K-band observations with models of H_2 collision-induced absorption in an atmosphere for a range of temperatures. While we find no evidence for deviation from a flat spectrum (slope s = 0.0016+/-0.0038), an H_2 dominated upper atmosphere (< 60 mbar) cannot be excluded. More precise observations at < 0.7 micron and in the K-band as well as a uniform analysis of all published data would be useful for establishing more robust limits on atmosphere models for GJ 1214b.

preprint2013arXiv

Testing the Metal of Late-Type Kepler Planet Hosts with Iron-Clad Methods

It has been shown that F, G, and early K dwarf hosts of Neptune-sized planets are not preferentially metal-rich. However, it is less clear whether the same holds for late K and M dwarf planet hosts. We report metallicities of Kepler targets and candidate transiting planet hosts with effective temperatures below 4500 K. We use new metallicity calibrations to determine [Fe/H] from visible and near-infrared spectra. We find that the metallicity distribution of late K and M dwarfs monitored by Kepler is consistent with that of the solar neighborhood. Further, we show that hosts of Earth- to Neptune-sized planets have metallicities consistent with those lacking detected planets and rule out a previously claimed 0.2 dex offset between the two distributions at 6sigma confidence. We also demonstrate that the metallicities of late K and M dwarfs hosting multiple detected planets are consistent with those lacking detected planets. Our results indicate that multiple terrestrial and Neptune-sized planets can form around late K and M dwarfs with metallicities as low as 0.25 of the solar value. The presence of Neptune-sized planets orbiting such low-metallicity M dwarfs suggests that accreting planets collect most or all of the solids from the disk and that the potential cores of giant planets can readily form around M dwarfs. The paucity of giant planets around M dwarfs compared to solar-type stars must be due to relatively rapid disk evaporation or a slower rate of core accretion, rather than insufficient solids to form a core.

preprint2012arXiv

Galactic Chemical Evolution and the Oxygen Isotopic Composition of the Solar System

We review current observational and theoretical constraints on the Galactic chemical evolution (GCE) of oxygen isotopes in order to explore whether GCE plays a role in explaining the lower 17O/18O ratio of the Sun, relative to the present-day interstellar medium, or the existence of distinct 16O-rich and 16O-poor reservoirs in the Solar System. Although the production of both 17O and 18O are related to the metallicity of progenitor stars, 17O is most likely produced in stars that evolve on longer timescales than those that produce 18O. Therefore the 17O/18O ratio need not have remained constant over time, contrary to preconceptions and the simplest models of GCE. An apparent linear, slope-one correlation between delta17O and delta18O in the ISM need not necessarily reflect an O isotopic gradient, and any slope-one galactocentric gradient need not correspond to evolution in time. Instead, increasing 17O/18O is consistent both with observational data from molecular clouds and with modeling of the compositions of presolar grains. Models in which the rate of star formation has decelerated over the past few Gyr or in which an enhanced period of star formation occurred shortly before solar birth ("starburst") can explain the solar-ISM O-isotopic difference without requiring a local input of supernova ejecta into the protosolar cloud. "Cosmic chemical memory" models in which interstellar dust is on average older than interstellar gas predict that primordial Solar System solids should be 16O-rich, relative to the Sun, in conflict with observations. However, scenarios in which the 16O-rich contribution of very massive stars could lead to 16O-poor solids and a 16O-rich bulk Sun, if the Solar System formed shortly after a starburst, independent of the popular scenario of photochemical self-shielding of CO.

preprint2012arXiv

Objects in Kepler's Mirror May be Larger Than They Appear: Bias and Selection Effects in Transiting Planet Surveys

Statistical analyses of large surveys for transiting planets such as the Kepler mission must account for systematic errors and biases. Transit detection depends not only on the planet's radius and orbital period, but also on host star properties. Thus, a sample of stars with transiting planets may not accurately represent the target population. Moreover, targets are selected using criteria such as a limiting apparent magnitude. These selection effects, combined with uncertainties in stellar radius, lead to biases in the properties of transiting planets and their host stars. We quantify possible biases in the Kepler survey. First, Eddington bias produced by a steep planet radius distribution and uncertainties in stellar radius results in a 15-20% overestimate of planet occurrence. Second, the magnitude limit of the Kepler target catalog induces Malmquist bias towards large, more luminous stars and underestimation of the radii of about one third of candidate planets, especially those larger than Neptune. Third, because metal-poor stars are smaller, stars with detected planets will be very slightly (<0.02 dex) more metal-poor than the target average. Fourth, uncertainties in stellar radii produce correlated errors in planet radius and stellar irradiation. A previous finding, that highly-irradiated giant are more likely to have "inflated" radii, remains significant, even accounting for this effect. In contrast, transit depth is negatively correlated with stellar metallicity even in the absence of any intrinsic correlation, and a previous claim of a negative correlation between giant planet transit depth and stellar metallicity is probably an artifact.

preprint2012arXiv

They might be giants: luminosity class, planet frequency, and planet-metallicity relation of the coolest Kepler target stars

We estimate the stellar parameters of late K and early M type Kepler target stars. We obtain medium resolution visible spectra of 382 stars with Kp-J>2 (~K5 and later spectral type). We determine luminosity class by comparing the strength of gravity-sensitive indices (CaH, K I, Ca II, and Na I) to their strength in a sample of stars of known luminosity class. We find that giants constitute 96+-% of the bright (Kp<14) Kepler target stars, and 7+-3% of dim (Kp>14) stars, significantly higher than fractions based on the stellar parameters quoted in the Kepler Input Catalog (KIC). The KIC effective temperatures are systematically (110 +15 -35} K) higher than temperatures we determine from fitting our spectra to PHOENIX stellar models. Through Monte Carlo simulations of the Kepler exoplanet candidate population, we find a planet occurrence of 0.36+-0.08 when giant stars are properly removed, somewhat higher than when a KIC log(g)>4 criterion is used (0.27+-0.05). Lastly, we show that there is no significant difference in g-r color (a probe of metallicity) between late-type Kepler stars with transiting Earth-to-Neptune sized exoplanet candidates and dwarf stars with no detected transits. We show that a previous claimed offset between these two populations is most likely an artifact of including a large number of misidentified giants.

preprint2011arXiv

An All-Sky Catalog of Bright M Dwarfs

We present an all-sky catalog of M dwarf stars with apparent infrared magnitude J<10. The 8,889 stars are selected from the ongoing SUPERBLINK survey of stars with proper motion >40 mas/yr, supplemented on the bright end with the TYCHO-2 catalog. Completeness tests which account for kinematic (proper motion) bias suggest that our catalog represents ~75% of the estimated ~11,900 M dwarfs with J<10 expected to populate the entire sky. Our catalog is, however, significantly more complete for the Northern sky (~90%) than it is for the South (~60%). Stars are identified as cool, red M dwarfs from a combination of optical and infrared color cuts, and are distinguished from background M giants and highly-reddened stars using either existing parallax measurements or, if such measurements are lacking, on their location in an optical-to-infrared reduced proper motion diagram. These bright M dwarfs are all prime targets for exoplanet surveys using the Doppler radial velocity or transit methods; the combination of low-mass and bright apparent magnitude should make possible the detection of Earth-size planets on short-period orbits using currently available techniques. Parallax measurements, when available, and photometric distance estimates are provided for all stars, and these place most systems within 60 parsecs of the Sun. Spectral type estimated from V-J color shows that most of the stars range from K7 to M4, with only a few late M dwarfs, all within 20 pc. Proximity to the Sun also makes these stars good targets for high-resolution exoplanet imaging searches, especially if younger objects can be identified on the basis of X-ray or UV excess. For that purpose, we include X-ray flux from ROSAT and FUV/NUV ultraviolet magnitudes from GALEX for all stars for which a counterpart can be identified in those catalogs.

preprint2011arXiv

Climate instability on tidally locked exoplanets

Feedbacks that can destabilize the climates of synchronously-rotating rocky planets may arise on planets with strong day-night surface temperature contrasts. Earth-like habitable-zone (HZ) planets maintain stable surface liquid water over geological time. This requires equilibrium between the temperature-dependent rate of greenhouse-gas consumption by weathering,and greenhouse-gas resupply by other processes. Detected small-radius exoplanets, and anticipated M-dwarf HZ rocky planets, are expected to be tidally locked. We investigate two feedbacks that can destabilize climate on tidally-locked planets. (1) If small changes in pressure alter the temperature distribution across a planet's surface such that the weathering rate increases when the pressure decreases, a positive feedback occurs involving increasing weathering rate near the substellar point, decreasing pressure, and increasing substellar surface temperature. (2) When decreases in pressure increase the surface area above the melting point (through reduced advective cooling of the substellar point), and the corresponding increase in volume of liquid causes net dissolution of the atmosphere, a further decrease in pressure occurs. We use an idealized energy balance model to map out the conditions under which these instabilities may occur. The weathering runaway can shrink the habitable zone, and cause geologically rapid 10^3-fold pressure shifts within the HZ. Mars may have undergone a weathering runaway in the past. Substellar dissolution is usually a negative feedback or weak positive feedback on changes in pressure. Both instabilities are suppressed if the atmosphere has a high radiative efficiency. Our results are most relevant for atmospheres that are thin and have low greenhouse-gas radiative efficiency. These results identify a new pathway by which HZ planets can undergo rapid climate shifts and become uninhabitable.

preprint2011arXiv

Differentiation of Planetesimals and the Thermal Consequences of Melt Migration

We model the heating of a primordial planetesimal by decay of the short-lived radionuclides Al-26 and Fe-60 to determine (i) the timescale on which melting will occur; (ii) the minimum size of a body that will produce silicate melt and differentiate; (iii) the migration rate of molten material within the interior; and (iv) the thermal consequences of the transport of Al-26 in partial melt. Our models incorporate results from previous studies of planetary differentiation and are constrained by petrologic (i.e. grain size distributions), isotopic (e.g. Pb-Pb and Hf-W ages) and mineralogical properties of differentiated achondrites. We show that formation of a basaltic crust via melt percolation was limited by the formation time of the body, matrix grain size and viscosity of the melt. We show that low viscosity (< 1 Pa-s) silicate melt can buoyantly migrate on a timescale comparable to the mean life of Al-26. The equilibrium partitioning of Al into silicate partial melt and the migration of that melt acts to dampen internal temperatures. However, subsequent heating from the decay of Fe-60 generated melt fractions in excess of 50%, thus completing differentiation for bodies that accreted within 2 Myr of CAI formation (i.e. the onset of isotopic decay). Migration and concentration of Al-26 into a crust results in remelting of that crust for accretion times less than 2 Myr and for bodies >100 km in size. Differentiation would be most likely for planetesimals larger than 20 km in diameter that accreted within ~2.7 Myr of CAI formation.

preprint2011arXiv

Ground-Based Submillimagnitude CCD Photometry of Bright Stars Using Snapshot Observations

We demonstrate ground-based submillimagnitude (<10^-3) photometry of widely separated bright stars using snapshot CCD imaging. We routinely achieved this photometric precision by (1) choosing nearby comparison stars of a similar magnitude and spectral type, (2) defocusing the telescope to allow high signal (>10^7 electrons) to be acquired in a single integration, (3) pointing the telescope so that all stellar images fall on the same detector pixels, and (4) using a region of the CCD detector that is free of nonlinear or aberrant pixels. We describe semiautomated observations with the Supernova Integrated Field Spectrograph (SNIFS) on the University of Hawaii 2.2m telescope on Mauna Kea, with which we achieved photometric precision as good as 5.2x10^-4 (0.56mmag) with a 5 minute cadence over a 2hr interval. In one experiment, we monitored eight stars, each separated by several degrees, and achieved submillimagnitude precision with a cadence (per star) of ~17 minutes. Our snapshot technique is suitable for automated searches for planetary transits among multiple bright stars.

preprint2011arXiv

Heterogeneous distribution of Al-26 at the birth of the Solar System

It is believed that Al-26, a short-lived (t1/2 = 0.73 Ma) and now extinct radionuclide, was uniformly distributed in the nascent Solar System with the initial Al-26/Al-27 ratio of ~5.2\times10-5, suggesting its external stellar origin. However, the stellar source of Al-26 and the manner in which it was injected into the solar system remain controversial: the Al-26 could have been produced by an asymptotic giant branch star, a supernova, or a Wolf-Rayet star and injected either into the protosolar molecular cloud or protoplanetary disk. Corundum (Al2O3) is thermodynamically predicted to be the first condensate from a cooling gas of solar composition. Here we show that micron-sized corundum condensates from O-16-rich gas (Big Delta O-17 ~ -25%) of solar composition recorded heterogeneous distribution of Al-26 at the birth of the solar system: the inferred initial Al-26/Al-27 ratio ranges from ~6.5x10-5 to <2x10-6; ~50% of the corundum grains measured are Al-26-poor. Other Al-26-poor, O-16-rich refractory objects include grossite (CaAl4O7)- and hibonite(CaAl12O19)-rich calcium-aluminum-rich inclusions (CAIs) in CH chondrites, platy hibonite crystals in CM chondrites, and FUN (fractionation and unidentified nuclear isotopic anomalies) CAIs in CV, CO, and CR chondrites. Considering the apparently early and short duration (<0.3 Ma) of condensation of refractory O-16-rich solids in the solar system, we infer that Al-26 was injected into the collapsing protosolar molecular cloud and later homogenized in the protoplanetary disk. The apparent lack of correlation between Al-26 abundance and O-isotope compositions of corundum grains put important constraints on the stellar source of Al-26 in the solar system.

preprint2011arXiv

Hydrogen Greenhouse Planets Beyond the Habitable Zone

We show that collision-induced absorption allows molecular hydrogen to act as an incondensible greenhouse gas, and that bars or tens of bars of primordial H2-He mixtures can maintain surface temperatures above the freezing point of water well beyond the "classical" habitable zone defined for CO2 greenhouse atmospheres. Using a 1-D radiative-convective model we find that 40 bars of pure H2 on a 3 Earth-mass planet can maintain a surface temperature of 280K out to 1.5AU from an early-type M dwarf star and 10 AU from a G-type star. Neglecting the effects of clouds and of gaseous absorbers besides H2, the flux at the surface would be sufficient for photosynthesis by cyanobacteria (in the G star case) or anoxygenic phototrophs (in the M star case). We argue that primordial atmospheres of one to several hundred bars of H2-He are possible, and use a model of hydrogen escape to show that such atmospheres are likely to persist further than 1.5 AU from M stars, and 2 AU from G stars, assuming these planets have protecting magnetic fields. We predict that the microlensing planet OGLE-05-390L could have retained a H2-He atmosphere and be habitable at ~2.6 AU from its host M star.

preprint2011arXiv

M2K: II. A Triple-Planet System Orbiting HIP 57274

Doppler observations from Keck Observatory have revealed a triple planet system orbiting the nearby mid-type K dwarf, HIP 57274. The inner planet, HIP 57274b, is a super-Earth with \msini\ = 11.6 \mearth (0.036 \mjup), an orbital period of 8.135 $\pm$ 0.004 d, and slightly eccentric orbit $e=0.19 \pm 0.1$. We calculate a transit probability of 6.5% for the inner planet. The second planet has \msini\ = 0.4 \mjup\ with an orbital period of 32.0 $\pm 0.02$ d in a nearly circular orbit, and $e = 0.05 \pm 0.03$. The third planet has \msini\ = 0.53 \mjup\ with an orbital period of 432 $\pm 8$ d (1.18 years) and an eccentricity $e = 0.23 \pm 0.03$. This discovery adds to the number of super Earth mass planets with $\msini < 12 \mearth\ $ that have been detected with Doppler surveys. We find that 56 $\pm 18$% super-Earths are members of multi-planet systems. This is certainly a lower limit because of observational detectability limits, yet significantly higher than the fraction of Jupiter mass exoplanets, $20 \pm 8$%, that are members of Doppler-detected, multi-planet systems.

preprint2011arXiv

Mantle Convection, Plate Tectonics, and Volcanism on Hot Exo-Earths

Recently discovered exoplanets on close-in orbits should have surface temperatures of 100's to 1000's of K. They are likely tidally locked and synchronously rotating around their parent stars and, if an atmosphere is absent, have surface temperature contrasts of many 100's to 1000's K between permanent day and night sides. We investigated the effect of elevated surface temperature and strong surface temperature contrasts for Earth-mass planets on the (i) pattern of mantle convection, (ii) tectonic regime, and (iii) rate and distribution of partial melting, using numerical simulations of mantle convection with a composite viscous/pseudo-plastic rheology. Our simulations indicate that, if a close-in rocky exoplanet lacks an atmosphere to redistribute heat, a >~ 400 K surface temperature contrast can maintain an asymmetric degree 1 pattern of mantle convection in which the surface of the planet moves preferentially toward subduction zones on the cold night side. The planetary surface features a hemispheric dichotomy, with plate-like tectonics on the night side and a continuously evolving mobile lid day side with diffuse surface deformation and vigorous volcanism. If volcanic outgassing establishes an atmosphere and redistributes heat, plate tectonics is globally replaced by diffuse surface deformation and volcanism accelerates and becomes distributed more uniformly across the planetary surface.

preprint2011arXiv

On the Nature of Small Planets around the Coolest Kepler Stars

We constrain the densities of Earth- to Neptune-size planets around very cool (Te =3660-4660K) Kepler stars by comparing 1202 Keck/HIRES radial velocity measurements of 150 nearby stars to a model based on Kepler candidate planet radii and a power-law mass-radius relation. Our analysis is based on the presumption that the planet populations around the two sets of stars are the same. The model can reproduce the observed distribution of radial velocity variation over a range of parameter values, but, for the expected level of Doppler systematic error, the highest Kolmogorov-Smirnov probabilities occur for a power-law index alpha ~ 4, indicating that rocky-metal planets dominate the planet population in this size range. A single population of gas-rich, low-density planets with alpha = 2 is ruled out unless our Doppler errors are >= 5m/s, i.e., much larger than expected based on observations and stellar chromospheric emission. If small planets are a mix of gamma rocky planets (alpha = 3.85) and 1-gamma gas-rich planets (alpha = 2), then gamma > 0.5 unless Doppler errors are >=4 m/s. Our comparison also suggests that Kepler's detection efficiency relative to ideal calculations is less than unity. One possible source of incompleteness is target stars that are misclassified subgiants or giants, for which the transits of small planets would be impossible to detect. Our results are robust to systematic effects, and plausible errors in the estimated radii of Kepler stars have only moderate impact.

preprint2010arXiv

M2K: I. A Jovian mass planet around the M3V star HIP79431

Doppler observations from Keck Observatory reveal the presence of a planet with Msini of 2.1 Mjup orbiting the M3V star HIP79431. This is the sixth giant planet to be detected in Doppler surveys of M dwarfs and it is one of the most massive planets discovered around an M dwarf star. The planet has an orbital period of 111.7 days and an orbital eccentricity of 0.29. The host star is metal rich, with an estimated [Fe/H] = +0.4. This is the first planet to emerge from our new survey of 1600 M-to-K dwarf stars.

preprint2010arXiv

Thermodynamic Limits on Magnetodynamos in Rocky Exoplanets

To ascertain whether magnetic dynamos operate in rocky exoplanets more massive or hotter than the Earth, we developed a parametric model of a differentiated rocky planet and its thermal evolution. Our model reproduces the established properties of Earth's interior and magnetic field at the present time. When applied to Venus, assuming that planet lacks plate tectonics and has a dehydrated mantle with an elevated viscosity, the model shows that the dynamo shuts down or never operated. Our model predicts that at a fixed planet mass, dynamo history is sensitive to core size, but not to the initial inventory of long-lived, heat-producing radionuclides. It predicts that rocky planets larger than 2.5 Earth masses will not develop inner cores because the temperature-pressure slope of the iron solidus becomes flatter than that of the core adiabat. Instead, iron "snow" will condense near or at the top of these cores, and the net transfer of latent heat upwards will suppress convection and a dynamo. More massive planets can have anemic dynamos due to core cooling, but only if they have mobile lids (plate tectonics). The lifetime of these dynamos is shorter with increasing planet mass but longer with higher surface temperature. Massive Venus-like planets with stagnant lids and more viscous mantles will lack dynamos altogether. We identify two alternative sources of magnetic fields on rocky planets: eddy currents induced in the hot or molten upper layers of planets on very short period orbits, and dynamos in the ionic conducting layers of "ocean" planets with ~10% mass in an upper mantle of water (ice).

preprint2009arXiv

Al-26 and the formation of the Solar System from a molecular cloud contaminated by Wolf-Rayet winds

In agreement with previous work, we show that the presence of the short-lived radionuclide Al-26 in the early Solar System was unlikely (<2% a priori probability) to be the result of direct introduction of supernova ejecta into the gaseous disk during the Class II stage of protosolar evolution. We also show that any Bondi-Hoyle accretion of contaminated residual gas from the natal star cluster made a negligible contribution to the primordial Al-26 inventory of the Solar System. These results are consistent with the absence of the oxygen isotopic signature expected with any late introduction of supernova ejecta into the protoplanetary disk. Instead, the presence of Al-26 in the oldest Solar System solids (calcium-aluminum-rich inclusions or CAIs) and its apparent uniform distribution with the inferred canonical Al-26/Al-27 ratio of (4.5-5)E-5 support the inheritance of Al-26 from the parent giant molecular cloud. We propose that this radionuclide originated in a prior generation of massive stars that formed in the same host molecular cloud as the Sun and contaminated that cloud by Wolf-Rayet winds. We calculated the Galactic distribution of Al-26/Al-27 ratios that arise from such contamination using the established embedded cluster mass and stellar initial mass functions, published nucleosynthetic yields from the winds of massive stars, and by assuming rapid and uniform mixing into the cloud. Although our model predicts that the majority of stellar systems contain no Al-26 from massive stars, and that the a priori probability that the Al-26/Al-27 ratio will reach or exceed the canonical Solar System value is only ~6%, the maximum in the distribution of non-zero values is close to the canonical ratio.

preprint2009arXiv

On the oxygen isotopic composition of the Solar System

The 18O/17O ratio of the Solar System is 5.2 while that of the interstellar medium (ISM) and young stellar objects is ~4. This difference cannot be explained by pollution of the Sun's natal molecular cloud by 18O-rich supernova ejecta because (1) the necessary B-star progenitors live longer than the duration of star formation in molecular clouds; (2) the delivery of ejecta gas is too inefficient and the amount of dust in supernova ejecta is too small compared to the required pollution (2% of total mass or ~20% of oxygen); and (3) the predicted amounts of concomitant short-lived radionuclides (SLRs) conflicts with the abundances of 26Al and 41Ca in the early Solar System. Proposals for the introduction of 18O-rich material must also be consistent with any explanation for the origin of the observed slope-one relationship between 17O/16O and 18O/16O in the high-temperature components of primitive meteorites. The difference in 18O/17O ratios can be explained by enrichment of the ISM by the 17O-rich winds of asymptotic giant branch (AGB) stars, the sequestration of comparatively 18O-rich gas from star-forming regions into long-lived, low-mass stars, and a monotonic decrease in the 18O/17O ratio of interstellar gas. At plausible rates of star formation and gas infall, Galactic chemical evolution does not follow a slope-one line in an three-isotope plot, but instead moves along a steeper trajectory towards an 17O-rich state. Evolution of the ISM and star-forming gas by AGB winds also explains the difference in the carbon isotope ratios of the Solar System and ISM.

preprint2008arXiv

The Effect of Lunar-like Satellites on the Orbital Infrared Light Curves of Earth-analog Planets

We investigate the influence of lunar-like satellites on the infrared orbital light curves of Earth-analog extra-solar planets. Such light curves will be obtained by NASA's Terrestrial Planet Finder (TPF) and ESA's Darwin missions as a consequence of repeat observations to confirm the companion status of a putative planet. We use an energy balance model to calculate disk-averaged infrared (bolometric) fluxes from planet-satellite systems over a full orbital period (one year). The satellites are assumed to lack an atmosphere, have a low thermal inertia like that of the Moon and span a range of plausible radii. The planets are assumed to have thermal and orbital properties that mimic those of the Earth while their obliquities and orbital longitudes of inferior conjunction remain free parameters. Even if the gross thermal properties of the planet can be independently constrained (e.g. via spectroscopy or visible-wavelength detection of specular glint from a surface ocean) only the largest (approximately Mars-size) lunar-like satellites can be detected by light curve data from a TPF-like instrument (i.e. one that achieves a photometric signal-to-noise of 10-20 at infrared wavelengths). Non-detection of a lunar-like satellite can obfuscate the interpretation of a given system's infrared light curve so that it may resemble a single planet with high obliquity, different orbital longitude of vernal equinox relative to inferior conjunction and in some cases drastically different thermal characteristics. If the thermal properties of the planet are not independently established then the presence of a lunar-like satellite cannot be inferred from infrared data, thus demonstrating that photometric light curves alone can only be used for preliminary study of extra-solar Earth-like planets.

preprint2003arXiv

A Survey of 10-Micron Silicate Emission from Dust around Young Sun-Like Stars

We obtained low resolution (R = 100) mid-infrared (8-13 micron wavelengths) spectra of 8 nearby young main sequence stars with the Keck 1 telescope and Long-Wavelength Spectrometer (LWS) to search for 10 micron silicate (Si-O stretch) emission from circumstellar dust. No stars exhibited readily apparent emission: Spectra were then analyzed by least-squares fitting of a template based on a spectrum of Comet Hale-Bopp. Using this technique, we were able to constrain the level of silicate emission to a threshold ten times below what was previously possible from space. We found one star, HD 17925, with a spectrum statistically different from its calibrator and consistent with a silicate emission peak of 7% of the photosphere at a wavelength of 10 microns. Excess emission at 60 microns from this star has already been reported.