-
Theoretical models of stars constitute a fundamental bedrock upon which much
of astrophysics is built, but large swaths of model parameter space remain
uncalibrated by observations. The best calibrators are eclipsing binaries in
clusters, allowing measurement of masses, radii, luminosities, and
temperatures, for stars of known metallicity and age. We present the discovery
and detailed characterization of PTFEB132.707+19.810, a P=6.0 day eclipsing
binary in the Praesepe cluster ($\tau$~600--800 Myr; [Fe/H]=0.14$\pm$0.04). The
system contains two late-type stars (SpT$_P$=M3.5$\pm$0.2;
SpT$_S$=M4.3$\pm$0.7) with precise masses ($M_p=0.3953\pm0.0020$~$M_{\odot}$;
$M_s=0.2098\pm0.0014$~$M_{\odot}$) and radii ($R_p=0.363\pm0.008$~$R_{\odot}$;
$R_s=0.272\pm0.012$~$R_{\odot}$). Neither star meets the predictions of stellar
evolutionary models. The primary has the expected radius, but is cooler and
less luminous, while the secondary has the expected luminosity, but is cooler
and substantially larger (by 20%). The system is not tidally locked or
circularized. Exploiting a fortuitous 4:5 commensurability between $P_{orb}$
and $P_{rot,prim}$, we demonstrate that fitting errors from the unknown spot
configuration only change the inferred radii by <1--2%. We also analyze subsets
of data to test the robustness of radius measurements; the radius sum is more
robust to systematic errors and preferable for model comparisons. We also test
plausible changes in limb darkening, and find corresponding uncertainties of
~1%. Finally, we validate our pipeline using extant data for GU Boo, finding
that our independent results match previous radii to within the mutual
uncertainties (2--3%). We therefore suggest that the substantial discrepancies
are astrophysical; since they are larger than for old field stars, they may be
tied to the intermediate age of PTFEB132.707+19.810.
-
When determining absolute ages of identifiably young stellar populations,
results strongly depend on which stars are studied. Cooler (K, M) stars
typically yield ages that are systematically younger than warmer (A, F, G)
stars by a factor of two. I explore the possibility that these age
discrepancies are the result of magnetic inhibition of convection in cool young
stars by using magnetic stellar evolution isochrones to determine the age of
the Upper Scorpius subgroup of the Scorpius-Centaurus OB Association. A median
age of 10 Myr consistent across spectral types A through M is found, except for
a subset of F-type stars that appear significantly older. Agreement is shown
for ages derived from the Hertzsprung-Russell diagram and from the empirical
mass-radius relationship defined by eclipsing multiple-star systems. Surface
magnetic field strengths required to produce agreement are of order 2.5 kG and
are predicted from a priori estimates of equipartition values. A region in the
HR diagram is identified that plausibly connects stars whose structures are
weakly influenced by the presence of magnetic fields with those whose
structures are strongly influenced by magnetic fields. The models suggest this
region is characterized by stars with rapidly thinning outer convective
envelopes where the radiative core mass is greater than 75% of the total
stellar mass. Furthermore, depletion of lithium predicted from magnetic models
appears in better agreement with observed lithium equivalent widths than
predictions from non-magnetic models. These results suggest that magnetic
inhibition of convection plays an important role in the early evolution of
low-mass stars and that it may be responsible for noted age discrepancies in
young stellar populations.
-
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.
-
Stellar evolution models are a cornerstone of young star astrophysics, which
necessitates that they yield accurate and reliable predictions of stellar
properties. Here, I review the current performance of stellar evolution models
against young astrophysical benchmarks and highlight recent progress
incorporating non-standard physics, such as magnetic field and starspots, to
explain observed deficiencies. While addition of these physical processes leads
to improved agreement between models and observations, there are several
fundamental limitations in our understanding about how these physical processes
operate. These limitations inhibit our ability to form a coherent picture of
the essential physics needed to accurately compute young stellar models, but
provide rich avenues for further exploration.
-
Stellar fundamental properties (masses, radii, effective temperatures) can be
extracted from observations of eclipsing binary systems with remarkable
precision, often better than 2%. Such precise measurements afford us the
opportunity to confront the validity of basic predictions of stellar evolution
theory, such as the mass-radius relationship. A brief historical overview of
confrontations between stellar models and data from eclipsing binaries is
given, highlighting key results and physical insight that have led directly to
our present understanding. The current paradigm that standard stellar evolution
theory is insufficient to describe the most basic relation, that of a star's
mass to its radius, along the main sequence is then described. Departures of
theoretical expectations from empirical data, however, provide a rich
opportunity to explore various physical solutions, improving our understanding
of important stellar astrophysical processes.
-
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.
-
We report accurate measurements of the physical properties (mass, radius,
temperature) of components of the G+M eclipsing binary V530 Ori. The M-type
secondary shows a larger radius and a cooler temperature than predicted by
standard stellar evolution models, as has been found for many other low-mass
stars and ascribed to the effects of magnetic activity and/or spots. We show
that models from the Dartmouth series that incorporate magnetic fields are able
to match the observations with plausible field strengths of 1-2 kG, consistent
with a rough estimate we derive for that star.
-
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.
-
We present direct radii measurements of the well-known transiting exoplanet
host stars HD 189733 and HD 209458 using the CHARA Array interferometer. We
find the limb-darkened angular diameters to be theta_LD = 0.3848 +/- 0.0055 and
0.2254 +/- 0.0072 milliarcsec for HD 189733 and HD 209458, respectively. HD
189733 and HD 209458 are currently the only two transiting exoplanet systems
where detection of the respective planetary companion's orbital motion from
high resolution spectroscopy has revealed absolute masses for both star and
planet. We use our new measurements together with the orbital information from
radial velocity and photometric time series data, Hipparcos distances, and
newly measured bolometric fluxes to determine the stellar effective
temperatures (T_eff = 4875 +/- 43, 6093 +/- 103 K), stellar linear radii (R_* =
0.805 +/- 0.016, 1.203 +/- 0.061 R_sun), mean stellar densities (rho_* = 1.62
+/- 0.11, 0.58 +/- 0.14 rho_sun), planetary radii (R_p = 1.216 +/- 0.024, 1.451
+/- 0.074 R_Jup), and mean planetary densities (rho_p = 0.605 +/- 0.029, 0.196
+/- 0.033 rho_Jup) for HD 189733 b and HD 209458 b, respectively. The stellar
parameters for HD 209458, a F9 dwarf, are consistent with indirect estimates
derived from spectroscopic and evolutionary modeling. However, we find that
models are unable to reproduce the observational results for the K2 dwarf, HD
189733. We show that, for stellar evolutionary models to match the observed
stellar properties of HD 189733, adjustments lowering the solar-calibrated
mixing length parameter from 1.83 to 1.34 need to be employed.
-
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.
-
We report an age revision for the low-mass detached eclipsing binary CM
Draconis and its common proper motion companion, WD 1633+572. An age of 8.5
$\pm$ 3.5 Gyr is found by combining an age estimate for the lifetime of WD
1633+572 and an estimate from galactic space motions. The revised age is
greater than a factor of two older than previous estimates. Our results provide
consistency between the white dwarf age and the system's galactic kinematics,
which reveal the system is a highly probable member of the galactic thick disk.
We find the probability that CM Draconis and WD 1633+572 are members of the
thick disk is 8500 times greater than the probability that they are members of
the thin disk and 170 times greater than the probability they are halo
interlopers. If CM Draconis is a member of the thick disk, it is likely
enriched in $\alpha$-elements compared to iron by at least 0.2 dex relative to
the Sun. This leads to the possibility that previous studies under-estimate the
[Fe/H] value, suggesting the system has a near-solar [Fe/H]. Implications for
the long-standing discrepancies between the radii of CM Draconis and
predictions from stellar evolution theory are discussed. We conclude that CM
Draconis is only inflated by about 2% compared to stellar evolution
predictions.
-
We present the current status of an effort to create an updated grid of
low-mass stellar evolution mass tracks and isochrones computed using the
Dartmouth stellar evolution code. Emphasis is placed on reliably extending the
present grid to the pre-main-sequence, where modeling uncertainties have the
greatest impact. Revisions to the original code release include: updated
surface boundary conditions, the introduction of deuterium burning, and
magnetic fields. The mass track grid contains models with a mass above 0.1
$M_{\odot}$ and metallicities in the range of -0.5 to +0.5 dex. Magnetic mass
tracks are calculated for surface magnetic field strengths between 0.1 kG and
4.0 kG using two different prescriptions for magneto-convection. Standard and
magnetic model isochrones are available for ages older than 1 Myr. Tabulated
quantities include the stellar fundamental properties, absolute photometric
magnitudes, magnetic field properties, convective turnover times, apsidal
motion constants, and lithium abundances. The complete grid of mass tracks and
isochrones will be made publicly available.
-
Based on high resolution optical spectra obtained with ESPaDOnS at CFHT, we
determine fundamental parameters (\Teff, R, \Lbol, \logg\ and metallicity) for
59 candidate members of nearby young kinematic groups. The candidates were
identified through the BANYAN Bayesian inference method of \citet{2013malo},
which takes into account the position, proper motion, magnitude, color, radial
velocity and parallax (when available) to establish a membership probability.
The derived parameters are compared to Dartmouth Magnetic evolutionary models
and to field stars with the goal to constrain the age of our candidates. We
find that, in general, low-mass stars in our sample are more luminous and have
inflated radii compared to older stars, a trend expected for pre-main sequence
stars. The Dartmouth Magnetic evolutionary models show a good fit to
observations of field K and M stars assuming a magnetic field strength of a few
kG, as typically observed for cool stars. Using the low-mass members of
$\beta$Pictoris moving group, we have re-examined the age inconsistency problem
between Lithium Depletion age and isochronal age (Hertzspring-Russell diagram).
We find that the inclusion of the magnetic field in evolutionary models
increase the isochronal age estimates for the K5V-M5V stars. Using these models
and field strengths, we derive an average isochronal age between 15 and 28 Myr
and we confirm a clear Lithium Depletion Boundary from which an age of
26$\pm$3~Myr is derived, consistent with previous age estimates based on this
method.
-
We present H- and K-band spectra for late-type Kepler Objects of Interest
(the "Cool KOIs"): low-mass stars with transiting-planet candidates discovered
by NASA's Kepler Mission that are listed on the NASA Exoplanet Archive. We
acquired spectra of 103 Cool KOIs and used the indices and calibrations of
Rojas-Ayala et al. to determine their spectral types, stellar effective
temperatures and metallicities, significantly augmenting previously published
values. We interpolate our measured effective temperatures and metallicities
onto evolutionary isochrones to determine stellar masses, radii, luminosities
and distances, assuming the stars have settled onto the main-sequence. As a
choice of isochrones, we use a new suite of Dartmouth predictions that reliably
include mid-to-late M dwarf stars. We identify five M4V stars: KOI-961
(confirmed as Kepler 42), KOI-2704, KOI-2842, KOI-4290, and the secondary
component to visual binary KOI-1725, which we call KOI-1725 B. We also identify
a peculiar star, KOI-3497, which has a Na and Ca lines consistent with a dwarf
star but CO lines consistent with a giant. Visible-wavelength adaptive optics
imaging reveals two objects within a 1 arc second diameter; however, the
objects' colors are peculiar. The spectra and properties presented in this
paper serve as a resource for prioritizing follow-up observations and planet
validation efforts for the Cool KOIs, and are all available for download online
using the "data behind the figure" feature.
-
We examine the hypothesis that magnetic fields are inflating the radii of
fully convective main sequence stars in detached eclipsing binaries (DEBs). The
magnetic Dartmouth stellar evolution code is used to analyze two systems in
particular: Kepler-16 and CM Draconis. Magneto-convection is treated assuming
stabilization of convection and also by assuming reductions in convective
efficiency due to a turbulent dynamo. We find that magnetic stellar models are
unable to reproduce the properties of inflated fully convective main sequence
stars, unless strong interior magnetic fields in excess of 10 MG are present.
Validation of the magnetic field hypothesis given the current generation of
magnetic stellar evolution models therefore depends critically on whether the
generation and maintenance of strong interior magnetic fields is physically
possible. An examination of this requirement is provided. Additionally, an
analysis of previous studies invoking the influence of star spots is presented
to assess the suggestion that star spots are inflating stars and biasing light
curve analyses toward larger radii. From our analysis, we find that there is
not yet sufficient evidence to definitively support the hypothesis that
magnetic fields are responsible for the observed inflation among fully
convective main sequence stars in DEBs.
-
Magnetic fields are hypothesized to inflate the radii of low-mass
stars---defined as less massive than 0.8$M_\odot$---in detached eclipsing
binaries (DEBs). We investigate this hypothesis using the recently introduced
magnetic Dartmouth stellar evolution code. In particular, we focus on stars
thought to have a radiative core and convective outer envelope by studying in
detail three individual DEBs: UV Psc, YY Gem, and CU Cnc. The results suggest
that the stabilization of thermal convection by a magnetic field is a plausible
explanation for the observed model-radius discrepancies. However, surface
magnetic field strengths required by the models are significantly stronger than
those estimated from the observed coronal X-ray emission. Agreement between
model predicted surface magnetic field strengths and those inferred from X-ray
observations can be found by assuming that the magnetic field sources its
energy from convection. This approach makes the transport of heat by convection
less efficient and is akin to reduced convective mixing length methods used in
other studies. Predictions for the metallicity and magnetic field strengths of
the aforementioned systems are reported. We also develop an expression relating
a reduction in the convective mixing length to a magnetic field strength in
units of the equipartition value. Our results are compared with those from
previous investigations to incorporate magnetic fields to explain the low-mass
DEB radius inflation. Finally, we explore how the effects of magnetic fields
might affect mass determinations using asteroseismic data and the implication
of magnetic fields on exoplanet studies.
-
We present results obtained using the magnetic Dartmouth stellar evolution
code that address the possibility that magnetic fields are inflating low-mass
stars in detached eclipsing binaries. While it seems plausible that magnetic
fields are inflating stars with radiative cores, the level of inflation
observed among fully convective stars appears too large to be explained by
magnetic fields. We provide an alternative explanation, stellar metallicity,
and propose observations that can help further constrain stellar models.
-
Magnetic fields have been hypothesized to inflate the radii of low-mass
stars---defined as less than 0.8 solar masses---in detached eclipsing binaries
(DEBs). We evaluate this hypothesis using the magnetic Dartmouth stellar
evolution code. Results suggest that magnetic suppression of thermal convection
can inflate low-mass stars that possess a radiative core and convective outer
envelope. A scaling relation between X-ray luminosity and surface magnetic flux
indicates that model surface magnetic field strength predictions are consistent
with observations. This supports the notion that magnetic fields may be
inflating these stars. However, magnetic models are unable to reproduce radii
of fully convective stars in DEBs. Instead, we propose that model discrepancies
below the fully convective boundary are related to metallicity.
-
We propose a novel method for determining the ages of low-mass, pre-main
sequence stellar systems using the apsidal motion of low-mass detached
eclipsing binaries. The apsidal motion of a binary system with an eccentric
orbit provides information regarding the interior structure constants of the
individual stars. These constants are related to the normalized stellar
interior density distribution and can be extracted from the predictions of
stellar evolution models. We demonstrate that low-mass, pre-main sequence stars
undergoing radiative core contraction display rapidly changing interior
structure constants (greater than 5% per 10 Myr) that, when combined with
observational determinations of the interior structure constants (with 5 -- 10%
precision), allow for a robust age estimate. This age estimate, unlike those
based on surface quantities, is largely insensitive to the surface layer where
effects of magnetic activity are likely to be most pronounced. On the main
sequence, where age sensitivity is minimal, the interior structure constants
provide a valuable test of the physics used in stellar structure models of
low-mass stars. There are currently no known systems where this technique is
applicable. Nevertheless, the emphasis on time domain astronomy with current
missions, such as Kepler, and future missions, such as LSST, has the potential
to discover systems where the proposed method will be observationally feasible.
-
We introduce a new one-dimensional stellar evolution code, based on the
existing Dartmouth code, that self-consistently accounts for the presence of a
globally pervasive magnetic field. The methods involved in perturbing the
equations of stellar structure, the equation of state, and the mixing-length
theory of convection are presented and discussed. As a first test of the code's
viability, stellar evolution models are computed for the components of a
solar-type, detached eclipsing binary (DEB) system, EF Aquarii, shown to
exhibit large disagreements with stellar models. The addition of the magnetic
perturbation corrects the radius and effective temperature discrepancies
observed in EF Aquarii. Furthermore, the required magnetic field strength at
the model photosphere is within a factor of two of the magnetic field strengths
estimated from the stellar X-ray luminosities measured by ROSAT and those
predicted from Ca II K line core emission. These models provide firm evidence
that the suppression of thermal convection arising from the presence of a
magnetic field is sufficient to significantly alter the structure of solar-type
stars, producing noticeably inflated radii and cooler effective temperatures.
The inclusion of magnetic effects within a stellar evolution model has a wide
range of applications, from DEBs and exoplanet host stars to the donor stars of
cataclysmic variables.
-
The CM Draconis system comprises two eclipsing mid-M dwarfs of nearly equal
mass in a 1.27-day orbit. This well-studied eclipsing binary has often been
used for benchmark tests of stellar models, since its components are amongst
the lowest mass stars with well-measured masses and radii (~ 1% relative
precision). However, as with many other low-mass stars, non-magnetic models
have been unable to match the observed radii and effective temperatures for CM
Dra at the 5-10% level. To date, the uncertain metallicity of the system has
complicated comparison of theoretical isochrones with observations. In this
Letter, we use data from the SpeX instrument on the NASA Infrared Telescope
Facility (IRTF) to measure the metallicity of the system during primary and
secondary eclipses, as well as out of eclipse, based on an empirical
metallicity calibration in the H and K near-infrared (NIR) bands. We derive a
[Fe/H] = -0.30 +- 0.12 that is consistent across all orbital phases. The
determination of [Fe/H] for this system constrains a key dimension of parameter
space when attempting to reconcile model isochrone predictions and
observations.
-
We examine the agreement between the observed and theoretical low-mass (< 0.8
solar masses) stellar main sequence mass-radius relationship by comparing
detached eclipsing binary (DEB) data with a new, large grid of stellar
evolution models. The new grid allows for a realistic variation in the age and
metallicity of the DEB population, characteristic of the local galactic
neighborhood. Overall, our models do a reasonable job of reproducing the
observational data. A large majority of the models match the observed stellar
radii to within 4%, with a mean absolute error of 2.3%. These results represent
a factor of two improvement compared to previous examinations of the low-mass
mass-radius relationship. The improved agreement between models and
observations brings the radius deviations within the limits imposed by
potential starspot-related uncertainties for 92% of the stars in our DEB
sample.
-
The recent discovery of an eclipsing hierarchical triple system with two
low-mass stars in a close orbit (KOI-126) by Carter et al. (2011) appeared to
reinforce the evidence that theoretical stellar evolution models are not able
to reproduce the observational mass-radius relation for low-mass stars. We
present a set of stellar models for the three stars in the KOI-126 system that
show excellent agreement with the observed radii. This agreement appears to be
due to the equation of state implemented by our code. A significant dispersion
in the observed mass-radius relation for fully convective stars is
demonstrated; indicative of the influence of physics currently not incorporated
in standard stellar evolution models. We also predict apsidal motion constants
for the two M-dwarf companions. These values should be observationally
determined to within 1% by the end of the Kepler mission.