Source author record

Keivan G. Stassun

Keivan G. Stassun 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

191works
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

191 published item(s)

preprint2026arXiv

Observable Metal Pollution in Main-Sequence Stars: Simulations of Rocky Planets Engulfed by Stars in the $0.5$ to $1.4$ M$_\odot$ Range

The engulfment of planets by their host stars is an expected outcome of various dynamical processes and has been invoked to explain a variety of observed stellar properties, such as rapid rotation, chemical abundance abnormalities, and other transient phenomena. Recent observations support engulfment as the cause of such signatures; however, many engulfment process details remain uncertain. Here, we present a model for determining the chemical signatures produced due to the pollution of main sequence stars by rocky planets, a common engulfment scenario due to the high frequency of observed short-period rocky exoplanets. A key novel element of our model is that we calculate the gradual evaporation of the planet due to drag interactions with the stellar envelope, which can lead to observable pollution on the stellar surface even if the bulk of the planet is only destroyed below the star's outer convective zone. Our results indicate that rocky planet pollution is most easily measurable for stars in the $1.0$ to $1.4$ M$_\odot$ range and that elements such as aluminium, calcium, and vanadium, in addition to lithium, are most suited to detect pollution. We predict that it is also possible to differentiate between the engulfment of one large planet versus several small planets, for the same total pollution mass, for some stellar hosts. We find that rocky planet engulfment events generally take years to decades for most stars. Our results can guide future observational campaigns that may search for sites of past or current engulfment events.

preprint2026arXiv

The Advanced X-ray Imaging Satellite (AXIS) Community Science Book

The AXIS Community Science Book represents the collective effort of 592 scientists worldwide to define the transformative science enabled by the Advanced X-ray Imaging Satellite (AXIS), a next-generation X-ray mission selected by NASA's Astrophysics Probe Program for Phase A study. AXIS will advance the legacy of high-angular-resolution X-ray astronomy with ~1.5'' imaging over a wide 24' field of view and an order of magnitude greater collecting area than Chandra in the 0.3-12 keV band. Combining sharp imaging, high throughput, and rapid response capabilities, AXIS will open new windows on virtually every aspect of modern astrophysics, exploring the birth and growth of supermassive black holes, the feedback processes that shape galaxies, the life cycles of stars and exoplanet environments, and the nature of compact stellar remnants, supernova remnants, and explosive transients. This book compiles 138 community-contributed science cases developed by five Science Working Groups focused on AGN and supermassive black holes, galaxy evolution and feedback, compact objects and supernova remnants, stellar physics and exoplanets, and time-domain and multi-messenger astrophysics. Together, these studies establish the scientific foundation for next-generation X-ray exploration in the 2030s and highlight strong synergies with facilities of the 2030s, such as JWST, Roman, Rubin/LSST, SKA, ALMA, ngVLA, and next-generation gravitational-wave and neutrino networks.

preprint2026arXiv

TOI-3862 b: A dense super-Neptune deep in the hot Neptune desert

The structure and evolution of close-in exoplanets are shaped by atmospheric loss and migration processes, which give rise to key population features such as the hot Neptune desert, ridge, and savanna - regions of the period-radius space whose boundaries offer critical insights into planetary formation and survival. As part of the KESPRINT collaboration, we selected the TESS transiting planet candidate TOI-3862.01 for radial velocity follow-up to confirm its planetary nature and characterize its mass and bulk properties. This planet candidate is of particular interest due to its position in the middle of the hot Neptune desert, making it a valuable probe for testing theories of planet migration and atmospheric loss. We confirmed the planetary nature and determined the mass of TOI-3862.01 (hereinafter TOI-3862b) by performing a joint fit with both transit and radial velocity data, precisely characterizing the bulk properties of this planet. TOI-3862b is a super-Neptune on a 1.56-day orbit around a Sun-like star with an effective temperature of 5300$\pm$50K. It has a mass of 53.7$_{-2.9}^{+2.8}$ M$_{\oplus}$ and a radius of 5.53$\pm$0.18 R$_{\oplus}$, corresponding to a density of 1.7$\pm$0.2 g/cm^3. This places it among the rare population of hot and dense super-Neptune desert planets. TOI-3862b, residing deep in the hot Neptune desert, represents a rare occurrence in an otherwise sparsely populated region, offering a valuable opportunity to probe the processes that may allow planets to survive in such environments.

preprint2026arXiv

Two temperate Earth- and Neptune-sized planets orbiting fully convective M dwarfs

As the diversity of exoplanets continues to grow, it is important to revisit assumptions about habitability and classical HZ definitions. In this work, we introduce an expanded 'temperate' zone, defined by instellation fluxes between $0.1<S/\mathrm{S}_\oplus<5$, thus encompassing a broader range of potentially habitable worlds. We also introduce the TEMPOS survey, which aims to produce a catalogue of precise radii for temperate planets orbiting M dwarfs with $T_\mathrm{eff}\leq3400\,$ K. This work reports the discovery and characterisation of two planets in this temperate regime orbiting mid-type M dwarfs: TOI-6716\,b, a $0.98\pm0.07\,\mathrm{R}_\oplus$ planet orbiting its M4 host star ($R_\star=0.231\,\pm0.015\mathrm{R}_\odot$, $M_\star=0.223\pm0.011\,\mathrm{M}_\odot$, $T_\mathrm{eff}=3110\pm80\,\mathrm{K}$) with a period $P=4.7185898^{+0.0000054}_{-0.0000041}\,\mathrm{d}$, and TOI-7384 b, a $3.56\pm0.21\,\mathrm{R}_\oplus$ planet orbiting an M4 ($R_\star=0.319\,\pm0.018\mathrm{R}_\odot$, $M_\star=0.318\pm0.016\,\mathrm{M}_\odot$, $T_\mathrm{eff}=3185\pm75\,\mathrm{K}$) star every $P=6.2340258^{+0.0000034}_{-0.0000036}\,\mathrm{d}$. The radii of TOI-6716 b and TOI-7384 b have precisions of $6.8\%$ and $5.9\%$ respectively. We validate these planets with multi-band ground-based photometric observations, high-resolution imaging and statistical analyses. We find these planets to have instellation fluxes close to the inner (hotter) edge of the temperate zone, with $4.4\pm1.1\,\mathrm{S}_\oplus$ and $4.9\pm1.1\,\mathrm{S}_\oplus$ for TOI-6716 b and TOI-7384 b respectively. Also, with a predicted TSM similar to the TRAPPIST-1 planets, TOI-6716 b is likely to be a good rocky-world JWST target, should it have retained its atmosphere.

preprint2024arXiv

ABYSS II: Identification of young stars in optical SDSS spectra and their properties

We developed a tool that measures equivalent widths of various lines in low resolution optical spectra, and it was applied to stellar spectra obtained as part of SDSS-V and LAMOST programs. These lines, such as Li I which directly indicates stellar youth, or optical H I and Ca II which in emission indicate activity associated with stellar youth, are commonly seen in YSOs. We observe several notable differences in the properties of these lines between YSOs and the field stars. Using these data, we devise a set of criteria through which it is possible to confirm the youth of stars that have been observed by the ABYSS program, as well as to identify likely young stars that have serendipitously been observed by other programs. We examine the decrement of H lines seen in emission in CTTSs, and estimate the properties of the accretion stream that is responsible for the production of these lines. Finally, we examine the evolution of Li I as a function of age, and characterize the scatter in its abundance that appears to be intrinsic in young M dwarfs.

preprint2024arXiv

Spectacular nucleosynthesis from early massive stars

Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys.

preprint2024arXiv

Weakened Magnetic Braking in the Exoplanet Host Star 51 Peg

The consistently low activity level of the old solar analog 51 Peg not only facilitated the discovery of the first hot Jupiter, but also led to the suggestion that the star could be experiencing a magnetic grand minimum. However, the 50 year time series showing minimal chromospheric variability could also be associated with the onset of weakened magnetic braking (WMB), where sufficiently slow rotation disrupts cycling activity and the production of large-scale magnetic fields by the stellar dynamo, thereby shrinking the Alfven radius and inhibiting the efficient loss of angular momentum to magnetized stellar winds. In this Letter, we evaluate the magnetic evolutionary state of 51 Peg by estimating its wind braking torque. We use new spectropolarimetric measurements from the Large Binocular Telescope to reconstruct the large-scale magnetic morphology, we reanalyze archival X-ray measurements to estimate the mass-loss rate, and we detect solar-like oscillations in photometry from the Transiting Exoplanet Survey Satellite, yielding precise stellar properties from asteroseismology. Our estimate of the wind braking torque for 51 Peg clearly places it in the WMB regime, driven by changes in the mass-loss rate and the magnetic field strength and morphology that substantially exceed theoretical expectations. Although our revised stellar properties have minimal consequences for the characterization of the exoplanet, they have interesting implications for the current space weather environment of the system.

preprint2023arXiv

Science with the Ultraviolet Explorer (UVEX)

UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program.

preprint2022arXiv

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

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

preprint2022arXiv

APOGEE Net: An expanded spectral model of both low mass and high mass stars

We train a convolutional neural network, APOGEE Net, to predict $T_\mathrm{eff}$, $\log g$, and, for some stars, [Fe/H], based on the APOGEE spectra. This is the first pipeline adapted for these data that is capable of estimating these parameters in a self-consistent manner not only for low mass stars, (such as main sequence dwarfs, pre-main sequence stars, and red giants), but also high mass stars with $T_\mathrm{eff}$ in excess of 50,000 K, including hot dwarfs and blue supergiants. The catalog of ~650,000 stars presented in this paper allows for a detailed investigation of the star forming history of not just the Milky Way, but also of the Magellanic clouds, as different type of objects tracing different parts of these galaxies can be more cleanly selected through their distinct placement in $T_\mathrm{eff}$-$\log g$ parameter space than in previous APOGEE catalogs produced through different pipelines.

preprint2022arXiv

Detailed Chemical Abundances for a Benchmark Sample of M Dwarfs from the APOGEE Survey

Individual chemical abundances for fourteen elements (C, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Ni) are derived for a sample of M-dwarfs using high-resolution near-infrared $H$-band spectra from the SDSS-IV/APOGEE survey. The quantitative analysis included synthetic spectra computed with 1-D LTE plane-parallel MARCS models using the APOGEE DR17 line list to determine chemical abundances. The sample consists of eleven M-dwarfs in binary systems with warmer FGK-dwarf primaries and ten measured interferometric angular diameters. To minimize atomic diffusion effects, [X/Fe] ratios are used to compare M-dwarfs in binary systems and literature results for their warmer primary stars, indicating good agreement ($<$0.08 dex) for all studied elements. The mean abundance differences in Primaries-this work M-dwarfs is -0.05$\pm$0.03 dex. It indicates that M-dwarfs in binary systems are a reliable way to calibrate empirical relationships. A comparison with abundance, effective temperature, and surface gravity results from the ASPCAP pipeline (DR16) finds a systematic offset of [M/H], $T_{\rm eff}$, log$g$ = +0.21 dex, -50 K, and 0.30 dex, respectively, although ASPCAP [X/Fe] ratios are generally consistent with this study. The metallicities of the M dwarfs cover the range of [Fe/H] = -0.9 to +0.4 and are used to investigate Galactic chemical evolution via trends of [X/Fe] as a function of [Fe/H]. The behavior of the various elemental abundances [X/Fe] versus [Fe/H] agrees well with the corresponding trends derived from warmer FGK-dwarfs, demonstrating that the APOGEE spectra can be used to Galactic chemical evolution using large samples of selected M-dwarfs.

preprint2022arXiv

Dynamical star forming history of Per OB2

We analyze the internal dynamics of young stars towards Perseus using Gaia EDR3 data, including Per OB2 and California Cloud. Interpreting the current dynamics, we speculate that Per OB2 may have formed from two separate clouds that have begun forming stars in a close proximity of each other. IC 348 is caught in the middle between the two of them, inheriting kinematics of both, and it stands out as a possible site of cloud cloud interaction. We also consider a possibility of a past supernova in Per OB2 - while one had likely occurred, it did not appear to have caused any obvious triggered star formation, but it has created a shock which has swept away the molecular gas away from IC 348. Finally, we examine a recently proposed shell between Taurus and Perseus. While its origin is unknown, we find no support for an expanding bubble in stellar kinematics, nor can we identify a likely progenitor for a supernova that may have caused it, disfavoring this scenario in the formation of this apparent shell.

preprint2022arXiv

HD 28109 hosts a trio of transiting Neptunian planets including a near-resonant pair, confirmed by ASTEP from Antarctica

We report on the discovery and characterisation of three planets orbiting the F8 star HD~28109, which sits comfortably in \tess's continuous viewing zone. The two outer planets have periods of $\rm 56.0067 \pm 0.0003~days$ and $\rm 84.2597_{-0.0008}^{+0.0010}~days$, which implies a period ratio very close to that of the first-order 3:2 mean motion resonance, exciting transit timing variations (TTVs) of up to $\rm 60\,mins$. These two planets were first identified by \tess, and we identified a third planet in the \textcolor{black}{\tess photometry} with a period of $\rm 22.8911 \pm 0.0004~days$. We confirm the planetary nature of all three planetary candidates using ground-based photometry from Hazelwood, ASTEP and LCO, including a full detection of the $\rm \sim9\,h$ transit of HD~28109 c from Antarctica. The radii of the three planets are \textcolor{black}{$\rm R_b=2.199_{-0.10}^{+0.098} ~R_{\oplus}$, $\rm R_c=4.23\pm0.11~ R_{\oplus}$ and $\rm R_d=3.25\pm0.11 ~R_{\oplus}$}; we characterise their masses using TTVs and precise radial velocities from ESPRESSO and HARPS, and find them to be $\rm M_b=18.5_{-7.6}^{+9.1}~M_{\oplus}$, $\rm M_c=7.9_{-3.0}^{+4.2}~M_{\oplus}$ and $\rm M_d=5.7_{-2.1}^{+2.7}~M_{\oplus}$, making planet b a dense, massive planet while c and d are both under-dense. We also demonstrate that the two outer planets are ripe for atmospheric characterisation using transmission spectroscopy, especially given their position in the CVZ of JWST. The data obtained to date are consistent with resonant (librating) and non-resonant (circulating) solutions; additional observations will show whether the pair is actually locked in resonance or just near-resonant.

preprint2022arXiv

HD 56414 b: A Warm Neptune Transiting an A-type Star

We report the discovery in TESS data and validation of HD 56414 b (a.k.a. TOI-1228 b), a Neptune-size ($R_{\rm p} = 3.71 \pm 0.20\, R_\oplus$) planet with a 29-day orbital period transiting a young (Age = $420 \pm 140$ Myr) A-type star in the TESS southern continuous viewing zone. HD 56414 is one of the hottest stars ($T_{\rm eff} = 8500 \pm 150 \, {\rm K}$) to host a known sub-Jovian planet. HD 56414 b lies on the boundary of the hot Neptune desert in planet radius -- bolometric insolation flux space, suggesting that the planet may be experiencing mass loss. To explore this, we apply a photoevaporation model that incorporates the high near ultraviolet continuum emission of A-type stars. We find that the planet can retain most of its atmosphere over the typical 1-Gyr main sequence lifetime of an A-type star if its mass is $\ge 8 \, M_\oplus$. Our model also predicts that close-in Neptune-size planets with masses $< 14 \, M_\oplus$ are susceptible to total atmospheric stripping over 1 Gyr, hinting that the hot Neptune desert, which has been previously observed around FGKM-type stars, likely extends to A-type stars.

preprint2022arXiv

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

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

preprint2022arXiv

Pre-main Sequence Brackett Emitters in the APOGEE DR17 Catalog: Line Strengths and Physical Properties of Accretion Columns

Very young (t $\lesssim$ 10 Myrs) stars possess strong magnetic fields that channel ionized gas from the interiors of their circumstellar discs to the surface of the star. Upon impacting the stellar surface, the shocked gas recombines and emits hydrogen spectral lines. To characterize the density and temperature of the gas within these accretion streams, we measure equivalent widths of Brackett (Br) 11-20 emission lines detected in 1101 APOGEE spectra of 326 likely pre-main sequence accretors. For sources with multiple observations, we measure median epoch-to-epoch line strength variations of 10% in Br11 and 20% in Br20. We also fit the measured line ratios to predictions of radiative transfer models by Kwan & Fischer. We find characteristic best-fit electron densities of $n_e$ = 10$^{11} - 10^{12}$ cm$^{-3}$, and excitation temperatures that are inversely correlated with electron density (from T$\sim$5000 K for $n_e \sim 10^{12}$ cm$^{-3}$, to T$\sim$12500 K at $n_e \sim 10^{11}$ cm$^{-3}$). These physical parameters are in good agreement with predictions from modelling of accretion streams that account for the hydrodynamics and radiative transfer within the accretion stream. We also present a supplementary catalog of line measurements from 9733 spectra of 4255 Brackett emission line sources in the APOGEE DR17 dataset.

preprint2022arXiv

Resolving the Multiplicity of Exoplanet Host Stars in Gemini/NIRI Data

The majority of stars have one or more stellar companions. As exoplanets continue to be discovered, it is crucial to examine planetary systems to identify their stellar companions. By observing a change in proper motion, companions can be detected by the acceleration they induce on their host stars. We selected 701 stars from the Hipparcos-Gaia Catalog of Accelerations (HGCA) that have existing adaptive optics imaging data gathered with Gemini/NIRI. Of these, we examined 21 stars known to host planet candidates and reduced their archival NIRI data with Gemini's DRAGONS software. We assessed these systems for companions using the NIRI images as well as Renormalized Unit Weight Error values in Gaia and accelerations in the HGCA. We detected three known visible companions and found two more systems with no visible companions but astrometric measurements indicating likely unresolved companions.

preprint2022arXiv

Rubin Observatory LSST Transients and Variable Stars Roadmap

The Vera C. Rubin Legacy Survey of Space and Time holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the Transient and Variable Stars Science Collaboration, one of the eight Rubin LSST Science Collaborations, has identified research areas of interest and requirements, and paths to enable them. While our roadmap is ever-evolving, this document represents a snapshot of our plans and preparatory work in the final years and months leading up to the survey's first light.

preprint2022arXiv

Solar-like oscillations and ellipsoidal variations in TESS observations of the binary 12 Boötis

Binary stars in which oscillations can be studied in either or both components can provide powerful constraints on our understanding of stellar physics. The bright binary 12 Boötis (12 Boo) is a particularly promising system because the primary is roughly 60 per cent brighter than the secondary despite being only a few per cent more massive. Both stars have substantial surface convection zones and are therefore, presumably, solar-like oscillators. We report here the first detection of solar-like oscillations and ellipsoidal variations in the TESS light curve of 12 Boo. Though the solar-like oscillations are not clear enough to unambiguously measure individual mode frequencies, we combine global asteroseismic parameters and a precise fit to the spectral energy distribution (SED) to provide new constraints on the properties of the system that are several times more precise than values in the literature. The SED fit alone provides new effective temperatures, luminosities and radii of $6115\pm45\,\mathrm{K}$, $7.531\pm0.110\,\mathrm{L}_\odot$ and $2.450\pm0.045\,\mathrm{R}_\odot$ for 12 Boo A and $6200\pm60\,\mathrm{K}$, $4.692\pm0.095\,\mathrm{L}_\odot$ and $1.901\pm0.045\,\mathrm{R}_\odot$ for 12 Boo B. When combined with our asteroseismic constraints on 12 Boo A, we obtain an age of $2.67^{+0.12}_{-0.16}\,\mathrm{Gyr}$, which is consistent with that of 12 Boo B.

preprint2022arXiv

Stellar multiplicity and stellar rotation: Insights from APOGEE

We measure rotational broadening in spectra taken by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey to characterise the relationship between stellar multiplicity and rotation. We create a sample of 2786 giants and 24 496 dwarfs with stellar parameters and multiple radial velocities from the APOGEE pipeline, projected rotation speeds \vsini\ determined from our own pipeline, and distances, masses, and ages measured by Sanders \& Das. We use the statistical distribution of the maximum shift in the radial velocities, \drvm, as a proxy for the close binary fraction to explore the interplay between stellar evolution, rotation, and multiplicity. Assuming that the minimum orbital period allowed is the critical period for Roche Lobe overflow and rotational synchronization, we calculate theoretical upper limits on expected \vsini\ and \drvm\ values. These expectations agree with the positive correlation between the maximum \drvm\ and \vsini\ values observed in our sample as a function of \logg. We find that the fast rotators in our sample have a high occurrence of short-period ($\log(P/\text{d})\lesssim 4$) companions. We also find that old, rapidly-rotating main sequence stars have larger completeness-corrected close binary fractions than their younger peers. Furthermore, rapidly-rotating stars with large \drvm\ consistently show differences of 1-10 Gyr between the predicted gyrochronological and measured isochronal ages. These results point towards a link between rapid rotation and close binarity through tidal interactions. We conclude that stellar rotation is strongly correlated with stellar multiplicity in the field, and caution should be taken in the application of gyrochronology relations to cool stars.

preprint2022arXiv

Stellar Properties for a Comprehensive Collection of Star Forming Regions in the SDSS APOGEE-2 Survey

The Sloan Digital Sky Survey IV (SDSS-IV) APOGEE-2 primary science goal was to observe red giant stars throughout the Galaxy to study its dynamics, morphology, and chemical evolution. The APOGEE instrument, a high-resolution 300 fiber H-band (1.55-1.71 micron) spectrograph, is also ideal to study other stellar populations in the Galaxy, among which are a number of star forming regions and young open clusters. We present the results of the determination of six stellar properties ($T_{eff}$, $\log{g}$, [Fe/H], $L/L_\odot$, $M/M_\odot$, and ages) for a sample that is composed of 3360 young stars, of sub-solar to super-solar types, in sixteen Galactic star formation and young open cluster regions. Those sources were selected by using a clustering method that removes most of the field contamination. Samples were also refined by removing targets affected by various systematic effects of the parameter determination. The final samples are presented in a comprehensive catalog that includes all six estimated parameters. This overview study also includes parameter spatial distribution maps for all regions and Hertzprung-Russell ($L/L_\odot$ vs. $T_{eff}$) diagrams. This study serves as a guide for detailed studies on individual regions, and paves the way for the future studies on the global properties of stars in the pre-main sequence phase of stellar evolution using more robust samples.

preprint2022arXiv

Systematic KMTNet Planetary Anomaly Search. V. Complete Sample of 2018 Prime-Field

We complete the analysis of all 2018 prime-field microlensing planets identified by the KMTNet AnomalyFinder. Among the 10 previously unpublished events with clear planetary solutions, 8 are either unambiguously planetary or are very likely to be planetary in nature: OGLE-2018-BLG-1126, KMT-2018-BLG-2004, OGLE-2018-BLG-1647, OGLE-2018-BLG-1367, OGLE-2018-BLG-1544, OGLE-2018-BLG-0932, OGLE-2018-BLG-1212, and KMT-2018-BLG-2718. Combined with the 4 previously published new AnomalyFinder events and 12 previously published (or in preparation) planets that were discovered by eye, thismakes a total of 24 2018 prime-field planets discovered or recovered by AnomalyFinder. Together with a paper in preparation on 2018 sub-prime planets, this work lays the basis for the first statistical analysis of the planet mass-ratio function based on planets identified in KMTNet data. By systematically applying the heuristic analysis of Hwang et al. (2022) to each event, we identify the small modification in their formalism that is needed to unify the so-called close/wide and inner/outer degeneracies, as conjectured by

preprint2022arXiv

The Discovery of a Planetary Companion Interior to Hot Jupiter WASP-132 b

Hot Jupiters are generally observed to lack close planetary companions, a trend that has been interpreted as evidence for high-eccentricity migration. We present the discovery and validation of WASP-132 c (TOI-822.02), a 1.85 $\pm$ 0.10 $R_{\oplus}$ planet on a 1.01 day orbit interior to the hot Jupiter WASP-132 b. Transiting Exoplanet Survey Satellite (TESS) and ground-based follow-up observations, in conjunction with vetting and validation analysis, enable us to rule out common astrophysical false positives and validate the observed transit signal produced by WASP-132 c as a planet. Running the validation tools \texttt{vespa} and \texttt{triceratops} on this signal yield false positive probabilities of $9.02 \times 10^{-5}$ and 0.0107, respectively. Analysis of archival CORALIE radial velocity data leads to a 3$σ$ upper limit of 28.23 ms$^{-1}$ on the amplitude of any 1.01-day signal, corresponding to a 3$σ$ upper mass limit of 37.35 $M_{\oplus}$. Dynamical simulations reveal that the system is stable within the 3$σ$ uncertainties on planetary and orbital parameters for timescales of $\sim$100 Myr. The existence of a planetary companion near the hot Jupiter WASP-132 b makes the giant planet's formation and evolution via high-eccentricity migration highly unlikely. Being one of just a handful of nearby planetary companions to hot Jupiters, WASP-132 c carries with it significant implications for the formation of the system and hot Jupiters as a population.

preprint2022arXiv

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

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

preprint2022arXiv

The Origin of Weakened Magnetic Braking in Old Solar Analogs

The rotation rates of main-sequence stars slow over time as they gradually lose angular momentum to their magnetized stellar winds. The rate of angular momentum loss depends on the strength and morphology of the magnetic field, the mass-loss rate, and the stellar rotation period, mass, and radius. Previous observations suggested a shift in magnetic morphology between two F-type stars with similar rotation rates but very different ages (88 Leo and rho CrB). In this Letter, we identify a comparable transition in an evolutionary sequence of solar analogs with ages between 2-7 Gyr. We present new spectropolarimetry of 18 Sco and 16 Cyg A & B from the Large Binocular Telescope, and we reanalyze previously published Zeeman Doppler images of HD 76151 and 18 Sco, providing additional constraints on the nature and timing of this transition. We combine archival X-ray observations with updated distances from Gaia to estimate mass-loss rates, and we adopt precise stellar properties from asteroseismology and other sources. We then calculate the wind braking torque for each star in the evolutionary sequence, demonstrating that the rate of angular momentum loss drops by more than an order of magnitude between the ages of HD 76151 and 18 Sco (2.6-3.7 Gyr) and continues to decrease modestly to the age of 16 Cyg A & B (7 Gyr). We suggest that this magnetic transition may represent a disruption of the global dynamo arising from weaker differential rotation, and we outline plans to probe this phenomenon in additional stars spanning a wide range of spectral types.

preprint2022arXiv

The TESS-Keck Survey. XI. Mass Measurements for Four Transiting sub-Neptunes orbiting K dwarf TOI-1246

Multi-planet systems are valuable arenas for investigating exoplanet architectures and comparing planetary siblings. TOI-1246 is one such system, with a moderately bright K dwarf ($\rm{V=11.6,~K=9.9}$) and four transiting sub-Neptunes identified by TESS with orbital periods of $4.31~\rm{d},~5.90~\rm{d},~18.66~\rm{d}$, and $~37.92~\rm{d}$. We collected 130 radial velocity observations with Keck/HIRES and TNG/HARPS-N to measure planet masses. We refit the 14 sectors of TESS photometry to refine planet radii ($\rm{2.97 \pm 0.06~R_\oplus},\rm{2.47 \pm 0.08~R_\oplus}, \rm{3.46 \pm 0.09~R_\oplus}$, $\rm{3.72 \pm 0.16~R_\oplus}$), and confirm the four planets. We find that TOI-1246 e is substantially more massive than the three inner planets ($\rm{8.1 \pm 1.1 M_\oplus}$, $\rm{8.8 \pm 1.2 M_\oplus}$, $\rm{5.3 \pm 1.7 M_\oplus}$, $\rm{14.8 \pm 2.3 M_\oplus}$). The two outer planets, TOI-1246 d and TOI-1246 e, lie near to the 2:1 resonance ($\rm{P_{e}/P_{d}=2.03}$) and exhibit transit timing variations. TOI-1246 is one of the brightest four-planet systems, making it amenable for continued observations. It is one of only six systems with measured masses and radii for all four transiting planets. The planet densities range from $\rm{0.70 \pm 0.24}$ to $3.21 \pm 0.44 \rm{g/cm^3}$, implying a range of bulk and atmospheric compositions. We also report a fifth planet candidate found in the RV data with a minimum mass of 25.6 $\pm$ 3.6 $\rm{M_\oplus}$. This planet candidate is exterior to TOI-1246 e with a candidate period of 93.8 d, and we discuss the implications if it is confirmed to be planetary in nature.

preprint2022arXiv

TKS V. Twin sub-Neptunes Transiting the Nearby G Star HD 63935

We present the discovery of two nearly identically-sized sub-Neptune transiting planets orbiting HD 63935, a bright ($V=8.6$ mag), sun-like ($T_{eff}=5560K$) star at 49 pc. TESS identified the first planet, HD 63935 b (TOI-509.01), in Sectors 7 and 34. We identified the second signal (HD 63935 c) in Keck HIRES and Lick APF radial velocity data as part of our followup campaign. It was subsequently confirmed with TESS photometry in Sector 34 as TOI-509.02. Our analysis of the photometric and radial velocity data yields a robust detection of both planets with periods of $9.0600 \pm 0.007$ and $21.40 \pm 0.0019$ days, radii of $2.99 \pm 0.14$ and $2.90 \pm 0.13$ $R_\oplus$, and masses of $10.8 \pm 1.8$ and $11.1 \pm 2.4$ $M_\oplus$. We calculate densities for planets b and c consistent with a few percent of the planet mass in hydrogen/helium envelopes. We also describe our survey's efforts to choose the best targets for JWST atmospheric followup. These efforts suggest that HD 63935 b will have the most clearly visible atmosphere of its class. It is the best target for transmission spectroscopy (ranked by Transmission Spectroscopy Metric, a proxy for atmospheric observability) in the so-far uncharacterized parameter space comprising sub-Neptune-sized (2.6 $R_\oplus$ $<$ $R_p$ $<$ 4 $R_\oplus$), moderately-irradiated (100 $F_\oplus$ $<$ $F_p$ $<$ 1000 $F_\oplus$) planets around G-stars. Planet c is also a viable target for transmission spectroscopy, and given the indistinguishable masses and radii of the two planets, the system serves as a natural laboratory for examining the processes that shape the evolution of sub-Neptune planets.

preprint2022arXiv

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

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

preprint2022arXiv

TOI-1452 b: SPIRou and TESS reveal a super-Earth in a temperate orbit transiting an M4 dwarf

Exploring the properties of exoplanets near or inside the radius valley provides insights on the transition from the rocky super-Earths to the larger, hydrogen-rich atmosphere mini-Neptunes. Here, we report the discovery of TOI-1452 b, a transiting super-Earth ($R_{\rm p} = 1.67 \pm 0.07$ R$_{\oplus}$) in an 11.1--day temperate orbit ($T_{\rm eq} = 326 \pm 7$ K) around the primary member ($H = 10.0$, $T_{\rm eff} = 3185 \pm 50$ K) of a nearby visual binary M dwarf. The transits were first detected by TESS, then successfully isolated between the two $3.2^{\prime\prime}$ companions with ground-based photometry from OMM and MuSCAT3. The planetary nature of TOI-1452 b was established through high-precision velocimetry with the near-infrared SPIRou spectropolarimeter as part of the ongoing SPIRou Legacy Survey. The measured planetary mass ($4.8 \pm 1.3$ M$_{\oplus}$) and inferred bulk density ($5.6^{+1.8}_{-1.6}$ g/cm$^3$) is suggestive of a rocky core surrounded by a volatile-rich envelope. More quantitatively, the mass and radius of TOI-1452 b, combined with the stellar abundance of refractory elements (Fe, Mg and Si) measured by SPIRou, is consistent with a core mass fraction of $18\pm6$ % and a water mass fraction of $22^{+21}_{-13}$%. The water world candidate TOI-1452 b is a prime target for future atmospheric characterization with JWST, featuring a Transmission Spectroscopy Metric similar to other well-known temperate small planets such as LHS 1140 b and K2-18 b. The system is located near Webb's northern Continuous Viewing Zone, implying that is can be followed at almost any moment of the year.

preprint2022arXiv

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

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

preprint2022arXiv

TOI-2119: A transiting brown dwarf orbiting an active M-dwarf from NASA's TESS mission

We report the discovery of TOI-2119b, a transiting brown dwarf (BD) that orbits and is completely eclipsed by an active M-dwarf star. Using light curve data from the Transiting Exoplanet Survey Satellite mission and follow-up high-resolution Doppler spectroscopic observations, we find the BD has a radius of $R_b = 1.08 \pm 0.03{\rm R_J}$, a mass of $M_b = 64.4 \pm 2.3{\rm M_J}$, an orbital period of $P = 7.200865 \pm 0.00002$ days, and an eccentricity of $e=0.337\pm 0.002$. The host star has a mass of $M_\star = 0.53 \pm 0.02{\rm M_\odot}$, a radius of $R_\star= 0.50 \pm 0.01{\rm R_\odot}$, an effective temperature of $T_{\rm eff} = 3621 \pm 48$K, and a metallicity of $\rm [Fe/H]=+0.06\pm 0.08$. TOI-2119b joins an emerging population of transiting BDs around M-dwarf host stars, with TOI-2119 being the ninth such system. These M-dwarf--brown dwarf systems typically occupy mass ratios near $q = M_b/M_\star \approx 0.1-0.2$, which separates them from the typical mass ratios for systems with transiting substellar objects and giant exoplanets that orbit more massive stars. The nature of the secondary eclipse of the BD by the star enables us to estimate the effective temperature of the substellar object to be $2030\pm 84$K, which is consistent with predictions by substellar evolutionary models.

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).

preprint2022arXiv

Twinkle -- a small satellite spectroscopy mission for the next phase of exoplanet science

With a focus on off-the-shelf components, Twinkle is the first in a series of cost competitive small satellites managed and financed by Blue Skies Space Ltd. The satellite is based on a high-heritage Airbus platform that will carry a 0.45 m telescope and a spectrometer which will provide simultaneous wavelength coverage from 0.5-4.5 $\rm{μm}$. The spacecraft prime is Airbus Stevenage while the telescope is being developed by Airbus Toulouse and the spectrometer by ABB Canada. Scheduled to begin scientific operations in 2025, Twinkle will sit in a thermally-stable, sun-synchronous, low-Earth orbit. The mission has a designed operation lifetime of at least seven years and, during the first three years of operation, will conduct two large-scale survey programmes: one focused on Solar System objects and the other dedicated to extrasolar targets. Here we present an overview of the architecture of the mission, refinements in the design approach, and some of the key science themes of the extrasolar survey.

preprint2022arXiv

White dwarf binaries Across the H-R Diagram

We created the APOGEE-GALEX-\emph{Gaia} catalog to study white dwarfs binaries. This database aims to create a minimally biased sample of WD binary systems identified from a combination of GALEX, {\it Gaia}, and APOGEE data to increase the number of WD binaries with orbital parameters and chemical compositions. We identify 3,414 sources as WD binary candidates, with nondegenerate companions of spectral types between F and M, including main sequence, main sequence binaries, subgiants, sub-subgiants, red giants, and red clump stars. Among our findings are (a) a total of 1,806 systems having inferred WD radii $R < 25$ R$_{\Earth}$, which constitute a more reliable group of WD binary candidates within the main sample; (b) a difference in the metallicity distribution function between WD binary candidates and the control sample of most luminous giants ($M_H < -3.0$); (c) the existence of a population of sub-subgiants with WD companions; (d) evidence for shorter periods in binaries that contain WDs compared to those that do not, as shown by the cumulative distributions of APOGEE radial velocity shifts; (e) evidence for systemic orbital evolution in a sample of 252 WD binaries with orbital periods, based on differences in the period distribution between systems with red clump, main sequence binary, and sub-subgiant companions and systems with main sequence or red giant companions; and (f) evidence for chemical enrichment during common envelope (CE) evolution, shown by lower metallicities in wide WD binary candidates ($P > 100$ days) compared to post-CE ($P < 100$ days) WD binary candidates.

preprint2021arXiv

A 20-Second Cadence View of Solar-Type Stars and Their Planets with TESS: Asteroseismology of Solar Analogs and a Re-characterization of pi Men c

We present an analysis of the first 20-second cadence light curves obtained by the TESS space telescope during its extended mission. We find a precision improvement of 20-second data compared to 2-minute data for bright stars when binned to the same cadence (~10-25% better for T<~8 mag, reaching equal precision at T~13 mag), consistent with pre-flight expectations based on differences in cosmic ray mitigation algorithms. We present two results enabled by this improvement. First, we use 20-second data to detect oscillations in three solar analogs (gamma Pav, zeta Tuc and pi Men) and use asteroseismology to measure their radii, masses, densities and ages to ~1%, ~3%, ~1% and ~20% respectively, including systematic errors. Combining our asteroseismic ages with chromospheric activity measurements we find evidence that the spread in the activity-age relation is linked to stellar mass and thus convection-zone depth. Second, we combine 20-second data and published radial velocities to re-characterize pi Men c, which is now the closest transiting exoplanet for which detailed asteroseismology of the host star is possible. We show that pi Men c is located at the upper edge of the planet radius valley for its orbital period, confirming that it has likely retained a volatile atmosphere and that the "asteroseismic radius valley" remains devoid of planets. Our analysis favors a low eccentricity for pi Men c (<0.1 at 68% confidence), suggesting efficient tidal dissipation (Q/k <~ 2400) if it formed via high-eccentricity migration. Combined, these early results demonstrate the strong potential of TESS 20-second cadence data for stellar astrophysics and exoplanet science.

preprint2021arXiv

A hot mini-Neptune in the radius valley orbiting solar analogue HD 110113

We report the discovery of HD 110113 b (TOI-755.01), a transiting mini-Neptune exoplanet on a 2.5-day orbit around the solar-analogue HD 110113 (Teff = 5730K). Using TESS photometry and HARPS radial velocities gathered by the NCORES program, we find HD 110113 b has a radius of $2.05\pm0.12$ $R_\oplus$ and a mass of $4.55\pm0.62$ $M_\oplus$. The resulting density of $2.90^{+0.75}_{-0.59}$ g cm^{-3} is significantly lower than would be expected from a pure-rock world; therefore, HD 110113 b must be a mini-Neptune with a significant volatile atmosphere. The high incident flux places it within the so-called radius valley; however, HD 110113 b was able to hold onto a substantial (0.1-1\%) H-He atmosphere over its $\sim4$ Gyr lifetime. Through a novel simultaneous gaussian process fit to multiple activity indicators, we were also able to fit for the strong stellar rotation signal with period $20.8\pm1.2$ d from the RVs and confirm an additional non-transiting planet with a mass of $10.5\pm1.2$ $M_\oplus$ and a period of $6.744^{+0.008}_{-0.009}$ d.

preprint2021arXiv

Chemical Cartography with APOGEE: Mapping Disk Populations with a Two-Process Model and Residual Abundances

We apply a novel statistical analysis to measurements of 16 elemental abundances in 34,410 Milky Way disk stars from the final data release (DR17) of APOGEE-2. Building on recent work, we fit median abundance ratio trends [X/Mg] vs. [Mg/H] with a 2-process model, which decomposes abundance patterns into a "prompt" component tracing core collapse supernovae and a "delayed" component tracing Type Ia supernovae. For each sample star, we fit the amplitudes of these two components, then compute the residuals Δ[X/H] from this two-parameter fit. The rms residuals range from ~0.01-0.03 dex for the most precisely measured APOGEE abundances to ~0.1 dex for Na, V, and Ce. The correlations of residuals reveal a complex underlying structure, including a correlated element group comprised of Ca, Na, Al, K, Cr, and Ce and a separate group comprised of Ni, V, Mn, and Co. Selecting stars poorly fit by the 2-process model reveals a rich variety of physical outliers and sometimes subtle measurement errors. Residual abundances allow comparison of populations controlled for differences in metallicity and [α/Fe]. Relative to the main disk (R=3-13 kpc, |Z|<2 kpc), we find nearly identical abundance patterns in the outer disk (R=15-17 kpc), 0.05-0.2 dex depressions of multiple elements in LMC and Gaia Sausage/Enceladus stars, and wild deviations (0.4-1 dex) of multiple elements in ωCen. Residual abundance analysis opens new opportunities for discovering chemically distinctive stars and stellar populations, for empirically constraining nucleosynthetic yields, and for testing chemical evolution models that include stochasticity in the production and redistribution of elements.

preprint2021arXiv

Final Targeting Strategy for the SDSS-IV APOGEE-2N Survey

APOGEE-2 is a dual-hemisphere, near-infrared (NIR), spectroscopic survey with the goal of producing a chemo-dynamical mapping of the Milky Way Galaxy. The targeting for APOGEE-2 is complex and has evolved with time. In this paper, we present the updates and additions to the initial targeting strategy for APOGEE-2N presented in Zasowski et al. (2017). These modifications come in two implementation modes: (i) "Ancillary Science Programs" competitively awarded to SDSS-IV PIs through proposal calls in 2015 and 2017 for the pursuit of new scientific avenues outside the main survey, and (ii) an effective 1.5-year expansion of the survey, known as the Bright Time Extension, made possible through accrued efficiency gains over the first years of the APOGEE-2N project. For the 23 distinct ancillary programs, we provide descriptions of the scientific aims, target selection, and how to identify these targets within the APOGEE-2 sample. The Bright Time Extension permitted changes to the main survey strategy, the inclusion of new programs in response to scientific discoveries or to exploit major new datasets not available at the outset of the survey design, and expansions of existing programs to enhance their scientific success and reach. After describing the motivations, implementation, and assessment of these programs, we also leave a summary of lessons learned from nearly a decade of APOGEE-1 and APOGEE-2 survey operations. A companion paper, Santana et al. (submitted), provides a complementary presentation of targeting modifications relevant to APOGEE-2 operations in the Southern Hemisphere.

preprint2021arXiv

Following up TESS Single Transits With Archival Photometry and Radial Velocities

NASA's Transiting Exoplanet Survey Satellite (TESS) mission is expected to discover hundreds of planets via single transits first identified in their light curves. Determining the orbital period of these single transit candidates typically requires a significant amount of follow-up work to observe a second transit or measure a radial velocity orbit. In Yao et al. (2019), we developed simulations that demonstrated the ability to use archival photometric data in combination with TESS to "precover" the orbital period for these candidates with a precision of several minutes, assuming circular orbits. In this work, we incorporate updated models for TESS single transits, allowing for eccentric orbits, along with an updated methodology to improve the reliability of the results. Additionally, we explore how radial velocity (RV) observations can be used to follow up single transit events, using strategies distinct from those employed when the orbital period is known. We find that the use of an estimated period based on a circular orbit to schedule reconnaissance RV observations can efficiently distinguish eclipsing binaries from planets. For candidates that pass reconnaissance RV observations, we simulate RV monitoring campaigns that enable one to obtain an approximate orbital solution. We find this method can regularly determine the orbital periods for planets more massive than 0.5 M_J with orbital periods as long as 100 days.

preprint2021arXiv

Gaia EDR3 confirms that Westerlund 1 is closer and older than previously thought

Using {\it Gaia} Early Data Release 3 (EDR3) parallaxes and Bayesian inference, we infer a parallax of the Westerlund 1 (Wd1) cluster. We find a parallax of $0.34\pm{0.05}$ mas corresponding to a distance of $2.8^{+0.7}_{-0.6}$ kpc. The new {\it Gaia} EDR3 distance is consistent with our previous result using {\it Gaia} DR2 parallaxes. This confirms that Wd1 is less massive and older than previously assumed. Compared to DR2, the EDR3 individual parallax uncertainties for each star decreased by 30\%. However, the aggregate parallax uncertainty for the cluster remained the same. This suggests that the uncertainty is dominated by systematics, which is possibly due to crowding, motions within the cluster, or motions due to binary orbits.

preprint2021arXiv

Parallax Systematics and Photocenter Motions of Benchmark Eclipsing Binaries in Gaia EDR3

Previous analyses of various standard candles observed by the Gaia satellite have reported statistically significant systematics in the parallaxes that have improved from $\sim$250 $μ$as in the first data release (DR1) to 50--80 $μ$as in the second data release (DR2). Here we examine the parallaxes newly reported in the Gaia early third data release (EDR3) using the same sample of benchmark eclipsing binaries (EBs) we used to assess the DR1 and DR2 parallaxes. We find a mean offset of $-37\pm20$ $μ$as (Gaia$-$EB), which decreases to $-15\pm18$ $μ$as after applying the corrections recommended by the Gaia Mission team; global systematics in the Gaia parallaxes have clearly improved and are no longer statistically significant for the EB sample, which spans $5\lesssim{G}\lesssim12$ in brightness and 0.03--3 kpc in distance. We also find that the RUWE goodness-of-fit statistic reported in Gaia EDR3 is highly sensitive to unresolved companions (tertiaries in the case of our EB sample) as well as to photocenter motion of the binaries themselves. RUWE is nearly perfectly correlated ($r^2=0.82$) with photocenter motions down to $\lesssim$0.1 mas, and surprisingly this correlation exists entirely within the nominal "good" RUWE range of 1.0--1.4. This suggests that RUWE values even slightly greater than 1.0 may signify unresolved binaries in Gaia, and that the RUWE value can serve as a quantitative predictor of the photocenter motion.

preprint2021arXiv

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

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

preprint2021arXiv

TESS Eclipsing Binary Stars. I. Short cadence observations of 4584 eclipsing binaries in Sectors 1-26

In this paper we present a catalog of 4584 eclipsing binaries observed during the first two years (26 sectors) of the TESS survey. We discuss selection criteria for eclipsing binary candidates, detection of hither-to unknown eclipsing systems, determination of the ephemerides, the validation and triage process, and the derivation of heuristic estimates for the ephemerides. Instead of keeping to the widely used discrete classes, we propose a binary star morphology classification based on a dimensionality reduction algorithm. Finally, we present statistical properties of the sample, we qualitatively estimate completeness, and discuss the results. The work presented here is organized and performed within the TESS Eclipsing Binary Working Group, an open group of professional and citizen scientists; we conclude by describing ongoing work and future goals for the group. The catalog is available from http://tessEBs.villanova.edu and from MAST.

preprint2021arXiv

The Impact of Observing Strategy on Reliable Classification of Standard Candle Stars: Detection of Amplitude, Period, and Phase Modulation (Blazhko Effect) of RR Lyrae Stars with LSST

The Vera C. Rubin Observatory will carry out its Legacy Survey of Space and Time (LSST) with a single-exposure depth of $r{\sim}24.7$ and an anticipated baseline of 10 years, allowing to access the Milky Way's old halo not only deeper, but also with a longer baseline and better cadence than e.g. PS1 3$π$ (Chambers et al. 2016). This will make LSST ideal to study populations of variable stars such as RR Lyrae stars (RRL). Here, we address the question of observing strategy optimization of LSST, as survey footprint definition, single visit exposure time as well as the cadence of repeat visits in different filters are yet to be finalized. We present metrics used to assess the impact of different observing strategies on the reliable detectability and classification of standard-candle variable stars, including detection of amplitude period, phase modulation effects of RR Lyrae stars, the so-called Blazhko effect (Blazhko 1907, Kollath et al. 2011), by evaluating metrics for simulated potential survey designs. So far, due to depth and cadence of typical all-sky surveys, it was nearly impossible to study this effect on a larger sample. All-sky surveys with relatively few observations over a moderately long baseline allow only for fitting phase-folded RRL light curves, thus integrating over the complete survey length and hiding any information regarding possible period or phase modulation during the survey. On the other hand, surveys with a cadence to detect slightly changing light curves usually have a relatively small footprint. LSST's survey strategy, however, will allow for studying variable stars in a way that makes population studies possible.

preprint2021arXiv

Thermal Phase Curves of XO-3b: an Eccentric Hot Jupiter at the Deuterium Burning Limit

We report \textit{Spitzer} full-orbit phase observations of the eccentric hot Jupiter XO-3b at 3.6 and 4.5 $μ$m. Our new eclipse depth measurements of $1770 \pm 180$ ppm at 3.6 $μ$m and $1610 \pm 70$ ppm at 4.5 $μ$m show no evidence of the previously reported dayside temperature inversion. We also empirically derive the mass and radius of XO-3b and its host star using Gaia DR3's parallax measurement and find a planetary mass $M_p=11.79 \pm 0.98 ~M_{\rm{Jup}}$ and radius $R_p=1.295 \pm 0.066 ~R_{\rm{Jup}}$. We compare our \textit{Spitzer} observations with multiple atmospheric models to constrain the radiative and advective properties of XO-3b. While the decorrelated 4.5 $μ$m observations are pristine, the 3.6 $μ$m phase curve remains polluted with detector systematics due to larger amplitude intrapixel sensitivity variations in this channel. We focus our analysis on the more reliable 4.5 $μ$m phase curve and fit an energy balance model with solid body rotation to estimate the zonal wind speed and the pressure of the bottom of the mixed layer. Our energy balance model fit suggests an eastward equatorial wind speed of $3.13 ^{+0.26} _{-0.83}$ km/s, an atmospheric mixed layer down to $2.40 ^{+0.92} _{-0.16}$ bar, and Bond albedo of $0.106 ^{+0.008} _{-0.106}$. We assume that the wind speed and mixed layer depth are constant throughout the orbit. We compare our observations with a 1D planet-averaged model predictions at apoapse and periapse and 3D general circulation model (GCM) predictions for XO-3b. We also investigate the inflated radius of XO-3b and find that it would require an unusually large amount of internal heating to explain the observed planetary radius.

preprint2021arXiv

TOI-1842b: A Transiting Warm Saturn Undergoing Re-Inflation around an Evolving Subgiant

The imminent launch of space telescopes designed to probe the atmospheres of exoplanets has prompted new efforts to prioritise the thousands of transiting planet candidates for follow-up characterisation. We report the detection and confirmation of TOI-1842b, a warm Saturn identified by TESS and confirmed with ground-based observations from Minerva-Australis, NRES, and the Las Cumbres Observatory Global Telescope. This planet has a radius of $1.04^{+0.06}_{-0.05}\,R_{Jup}$, a mass of $0.214^{+0.040}_{-0.038}\,M_{Jup}$, an orbital period of $9.5739^{+0.0002}_{-0.0001}$ days, and an extremely low density ($ρ$=0.252$\pm$0.091 g cm$^{-3}$). TOI-1842b has among the best known combinations of large atmospheric scale height (893 km) and host-star brightness ($J=8.747$ mag), making it an attractive target for atmospheric characterisation. As the host star is beginning to evolve off the main sequence, TOI-1842b presents an excellent opportunity to test models of gas giant re-inflation. The primary transit duration of only 4.3 hours also makes TOI-1842b an easily-schedulable target for further ground-based atmospheric characterisation.

preprint2021arXiv

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

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

preprint2021arXiv

TOI-954 b and K2-329 b: Short-Period Saturn-Mass Planets that Test whether Irradiation Leads to Inflation

We report the discovery of two short-period Saturn-mass planets, one transiting the G subgiant TOI-954 (TIC 44792534, $ V = 10.343 $, $ T = 9.78 $) observed in TESS sectors 4 and 5, and one transiting the G dwarf K2-329 (EPIC 246193072, $ V = 12.70 $, $ K = 10.67 $) observed in K2 campaigns 12 and 19. We confirm and characterize these two planets with a variety of ground-based archival and follow-up observations, including photometry, reconnaissance spectroscopy, precise radial velocity, and high-resolution imaging. Combining all available data, we find that TOI-954 b has a radius of $0.852_{-0.062}^{+0.053} \, R_{\mathrm{J}}$ and a mass of $0.174_{-0.017}^{+0.018} \, M_{\mathrm{J}}$ and is in a 3.68 day orbit, while K2-329 b has a radius of $0.774_{-0.024}^{+0.026} \, R_{\mathrm{J}}$ and a mass of $0.260_{-0.022}^{+0.020} \, M_{\mathrm{J}}$ and is in a 12.46 day orbit. As TOI-954 b is 30 times more irradiated than K2-329 b but more or less the same size, these two planets provide an opportunity to test whether irradiation leads to inflation of Saturn-mass planets and contribute to future comparative studies that explore Saturn-mass planets at contrasting points in their lifetimes.

preprint2021arXiv

Weighing stars from birth to death: mass determination methods across the HRD

The mass of a star is the most fundamental parameter for its structure, evolution, and final fate. It is particularly important for any kind of stellar archaeology and characterization of exoplanets. There exists a variety of methods in astronomy to estimate or determine it. In this review we present a significant number of such methods, beginning with the most direct and model-independent approach using detached eclipsing binaries. We then move to more indirect and model-dependent methods, such as the quite commonly used isochrone or stellar track fitting. The arrival of quantitative asteroseismology has opened a completely new approach to determine stellar masses and to complement and improve the accuracy of other methods. We include methods for different evolutionary stages, from the pre-main sequence to evolved (super)giants and final remnants. For all methods uncertainties and restrictions will be discussed. We provide lists of altogether more than 200 benchmark stars with relative mass accuracies between $[0.3,2]\%$ for the covered mass range of $M\in [0.1,16]\,\msun$, $75\%$ of which are stars burning hydrogen in their core and the other $25\%$ covering all other evolved stages. We close with a recommendation how to combine various methods to arrive at a "mass-ladder" for stars.

preprint2020arXiv

A Giant Planet Candidate Transiting a White Dwarf

Astronomers have discovered thousands of planets outside the solar system, most of which orbit stars that will eventually evolve into red giants and then into white dwarfs. During the red giant phase, any close-orbiting planets will be engulfed by the star, but more distant planets can survive this phase and remain in orbit around the white dwarf. Some white dwarfs show evidence for rocky material floating in their atmospheres, in warm debris disks, or orbiting very closely, which has been interpreted as the debris of rocky planets that were scattered inward and tidally disrupted. Recently, the discovery of a gaseous debris disk with a composition similar to ice giant planets demonstrated that massive planets might also find their way into tight orbits around white dwarfs, but it is unclear whether the planets can survive the journey. So far, the detection of intact planets in close orbits around white dwarfs has remained elusive. Here, we report the discovery of a giant planet candidate transiting the white dwarf WD 1856+534 (TIC 267574918) every 1.4 days. The planet candidate is roughly the same size as Jupiter and is no more than 14 times as massive (with 95% confidence). Other cases of white dwarfs with close brown dwarf or stellar companions are explained as the consequence of common-envelope evolution, wherein the original orbit is enveloped during the red-giant phase and shrinks due to friction. In this case, though, the low mass and relatively long orbital period of the planet candidate make common-envelope evolution less likely. Instead, the WD 1856+534 system seems to demonstrate that giant planets can be scattered into tight orbits without being tidally disrupted, and motivates searches for smaller transiting planets around white dwarfs.

preprint2020arXiv

A KELT-TESS Eclipsing Binary in a Young Triple System Associated with a "Stellar String" Theia 301

HD 54236 is a nearby, wide common-proper-motion visual pair that has been previously identified as likely being very young by virtue of strong X-ray emission and lithium absorption. Here we report the discovery that the brighter member of the wide pair, HD~54236A, is itself an eclipsing binary (EB), comprising two near-equal solar-mass stars on a 2.4 d orbit. It represents a potentially valuable opportunity to expand the number of benchmark-grade EBs at young stellar ages. Using new observations of Ca2H&K emission and lithium absorption in the wide K-dwarf companion, HD 54236B, we obtain a robust age estimate of 225 +/- 50 Myr for the system. This age estimate and Gaia proper motions show HD 54236 is associated with Theia~301, a newly discovered local "stellar string", which itself may be related to the AB Dor moving group through shared stellar members. Applying this age estimate to AB~Dor itself alleviates reported tension between observation and theory that arises for the luminosity of the 90M_Jup star/brown dwarf AB Dor C when younger age estimates are used.

preprint2020arXiv

A remnant planetary core in the hot-Neptune desert

The interiors of giant planets remain poorly understood. Even for the planets in the Solar System, difficulties in observation lead to large uncertainties in the properties of planetary cores. Exoplanets that have undergone rare evolutionary processes provide a route to understanding planetary interiors. Planets found in and near the typically barren hot-Neptune 'desert' (a region in mass-radius space that contains few planets) have proved to be particularly valuable in this regard. These planets include HD149026b, which is thought to have an unusually massive core, and recent discoveries such as LTT9779b and NGTS-4b, on which photoevaporation has removed a substantial part of their outer atmospheres. Here we report observations of the planet TOI-849b, which has a radius smaller than Neptune's but an anomalously large mass of $39.1^{+2.7}_{-2.6}$ Earth masses and a density of $5.2^{+0.7}_{-0.8}$ grams per cubic centimetre, similar to Earth's. Interior structure models suggest that any gaseous envelope of pure hydrogen and helium consists of no more than $3.9^{+0.8}_{-0.9}$ per cent of the total planetary mass. The planet could have been a gas giant before undergoing extreme mass loss via thermal self-disruption or giant planet collisions, or it could have avoided substantial gas accretion, perhaps through gap opening or late formation. Although photoevaporation rates cannot account for the mass loss required to reduce a Jupiter-like gas giant, they can remove a small (a few Earth masses) hydrogen and helium envelope on timescales of several billion years, implying that any remaining atmosphere on TOI-849b is likely to be enriched by water or other volatiles from the planetary interior. We conclude that TOI-849b is the remnant core of a giant planet.

preprint2020arXiv

Age dating of an early Milky Way merger via asteroseismology of the naked-eye star $ν$ Indi

Over the course of its history, the Milky Way has ingested multiple smaller satellite galaxies. While these accreted stellar populations can be forensically identified as kinematically distinct structures within the Galaxy, it is difficult in general to precisely date the age at which any one merger occurred. Recent results have revealed a population of stars that were accreted via the collision of a dwarf galaxy, called \textit{Gaia}-Enceladus, leading to a substantial pollution of the chemical and dynamical properties of the Milky Way. Here, we identify the very bright, naked-eye star $ν$\,Indi as a probe of the age of the early in situ population of the Galaxy. We combine asteroseismic, spectroscopic, astrometric, and kinematic observations to show that this metal-poor, alpha-element-rich star was an indigenous member of the halo, and we measure its age to be $11.0 \pm 0.7$ (stat) $\pm 0.8$ (sys)$\,\rm Gyr$. The star bears hallmarks consistent with it having been kinematically heated by the \textit{Gaia}-Enceladus collision. Its age implies that the earliest the merger could have begun was 11.6 and 13.2 Gyr ago at 68 and 95% confidence, respectively. Input from computations based on hierarchical cosmological models tightens (i.e. reduces) slightly the above limits.

preprint2020arXiv

Analysis of Previously Classified White Dwarf-Main Sequence Binaries Using Data from the APOGEE Survey

We present analyses of near-infrared, spectroscopic data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey for 45 previously confirmed or candidate white dwarf-main sequence (WDMS) binaries identified by the optical SDSS and LAMOST surveys. Among these 45 systems, we classify three as having red giant primaries in the LAMOST sample and fourteen to be young stellar object contaminants in the photometrically identified SDSS sample. From among the subsample of 28 systems that we confirm to have MS primaries, we derive and place limits on orbital periods and velocity amplitudes for fourteen. Seven systems have significant velocity variations that warrant a post-common-envelope (PCE) binary classification -- four of which are newly classified, three of which are newly confirmed, and five for which we can derive full orbital parameters. If confirmed, one of these newly discovered systems (2M14544500+4626456) will have the second longest orbital period reported for a typical, compact PCE WDMS binary ($P=15.1$ days). In addition to the seven above, we also recover and characterize with APOGEE data the well known PCE WDMS systems EG UMa and HZ 9. We also investigate the overall metallicity distribution of the WDMS sample, which is a parameter space not often explored for these systems. Of note, we find one system (2M14244053+4929580) to be extremely metal-poor (${\rm [Fe/H]}=-1.42$) relative to the rest of the near-solar sample. Additionally, the PCE systems in our sample are found to be, on average, higher in metallicity than their wide-binary counterparts, though we caution that with this small number of systems, the sample may not be representative of the overall distribution of WDMS systems.

preprint2020arXiv

APOGEE Net: Improving the derived spectral parameters for young stars through deep learning

Machine learning allows efficient extraction of physical properties from stellar spectra that have been obtained by large surveys. The viability of ML approaches has been demonstrated for spectra covering a variety of wavelengths and spectral resolutions, but most often for main sequence or evolved stars, where reliable synthetic spectra provide labels and data for training. Spectral models of young stellar objects (YSOs) and low mass main sequence (MS) stars are less well-matched to their empirical counterparts, however, posing barriers to previous approaches to classify spectra of such stars. In this work we generate labels for YSOs and low mass MS stars through their photometry. We then use these labels to train a deep convolutional neural network to predict log g, Teff, and Fe/H for stars with APOGEE spectra in the DR14 dataset. This "APOGEE Net" has produced reliable predictions of log g for YSOs, with uncertainties of within 0.1 dex and a good agreement with the structure indicated by pre-main sequence evolutionary tracks, and correlate well with independently derived stellar radii. These values will be useful for studying pre-main sequence stellar populations to accurately diagnose membership and ages.

preprint2020arXiv

Close Binary Companions to APOGEE DR16 Stars: 20,000 Binary-star Systems Across the Color-Magnitude Diagram

Many problems in contemporary astrophysics---from understanding the formation of black holes to untangling the chemical evolution of galaxies---rely on knowledge about binary stars. This, in turn, depends on discovery and characterization of binary companions for large numbers of different kinds of stars in different chemical and dynamical environments. Current stellar spectroscopic surveys observe hundreds of thousands to millions of stars with (typically) few observational epochs, which allows binary discovery but makes orbital characterization challenging. We use a custom Monte Carlo sampler (The Joker) to perform discovery and characterization of binary systems through radial-velocities, in the regime of sparse, noisy, and poorly sampled multi-epoch data. We use it to generate posterior samplings in Keplerian parameters for 232,531 sources released in APOGEE Data Release 16. Our final catalog contains 19,635 high-confidence close-binary (P < few years, a < few AU) systems that show interesting relationships between binary occurrence rate and location in the color-magnitude diagram. We find notable faint companions at high masses (black-hole candidates), at low masses (substellar candidates), and at very close separations (mass-transfer candidates). We also use the posterior samplings in a (toy) hierarchical inference to measure the long-period binary-star eccentricity distribution. We release the full set of posterior samplings for the entire parent sample of 232,531 stars. This set of samplings involves no heuristic "discovery" threshold and therefore can be used for myriad statistical purposes, including hierarchical inferences about binary-star populations and sub-threshold searches.

preprint2020arXiv

Evidence for H$_{2}$ Dissociation and Recombination Heat Transport in the Atmosphere of KELT-9b

Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $μ$m. KELT-9b is the hottest known transiting planet, with a 4.5-$μ$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5σ$ confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.

preprint2020arXiv

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

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

preprint2020arXiv

HD 191939: Three Sub-Neptunes Transiting a Sun-like Star Only 54 pc Away

We present the discovery of three sub-Neptune-sized planets transiting the nearby and bright Sun-like star HD 191939 (TIC 269701147, TOI 1339), a $K_{s}=7.18$ magnitude G8 V dwarf at a distance of only 54 parsecs. We validate the planetary nature of the transit signals by combining five months of data from the Transiting Exoplanet Survey Satellite with follow-up ground-based photometry, archival optical images, radial velocities, and high angular resolution observations. The three sub-Neptunes have similar radii ($R_{b} = 3.42^{+0.11}_{-0.11}\,R_{\oplus}$, $R_{c}=3.23_{-0.11}^{+0.11}\,R_{\oplus}$, and $R_{d}=3.16_{-0.11}^{+0.11}\,R_{\oplus}$) and their orbits are consistent with a stable, circular, and co-planar architecture near mean motion resonances of 1:3 and 3:4 ($P_{b}=8.88$ days, $P_{c}=28.58$ days, and $P_{d}=38.35$ days). The HD~191939 system is an excellent candidate for precise mass determinations of the planets with high-resolution spectroscopy due to the host star's brightness and low chromospheric activity. Moreover, the system's compact and near-resonant nature can provide an independent way to measure planetary masses via transit timing variations while also enabling dynamical and evolutionary studies. Finally, as a promising target for multi-wavelength transmission spectroscopy of all three planets' atmospheres, HD 191939 can offer valuable insight into multiple sub-Neptunes born from a proto-planetary disk that may have resembled that of the early Sun.

preprint2020arXiv

Inferring the parallax of Westerlund 1 from Gaia DR2

Westerlund 1 (Wd1) is potentially the largest star cluster in the Galaxy. That designation critically depends upon the distance to the cluster, yet the cluster is highly obscured, making luminosity-based distance estimates difficult. Using {\it Gaia} Data Release 2 (DR2) parallaxes and Bayesian inference, we infer a parallax of $0.35^{+0.07}_{-0.06}$ mas corresponding to a distance of $2.6^{+0.6}_{-0.4}$ kpc. To leverage the combined statistics of all stars in the direction of Wd1, we derive the Bayesian model for a cluster of stars hidden among Galactic field stars; this model includes the parallax zero-point. Previous estimates for the distance to Wd1 ranged from 1.0 to 5.5 kpc, although values around 5 kpc have usually been adopted. The {\it Gaia} DR2 parallaxes reduce the uncertainty from a factor of 3 to 18\% and rules out the most often quoted value of 5 kpc with 99\% confidence. This new distance allows for more accurate mass and age determinations for the stars in Wd1. For example, the previously inferred initial mass at the main-sequence turn-off was around 40 M$_{\odot}$; the new {\it Gaia} DR2 distance shifts this down to about 22 M$_{\odot}$. This has important implications for our understanding of the late stages of stellar evolution, including the initial mass of the magnetar and the LBV in Wd1. Similarly, the new distance suggests that the total cluster mass is about four times lower than previously calculated.

preprint2020arXiv

LHS 1815b: The First Thick-Disk Planet Detected By TESS

We report the first discovery of a thick-disk planet, LHS 1815b (TOI-704b, TIC 260004324), detected in the TESS survey. LHS 1815b transits a bright (V = 12.19 mag, K = 7.99 mag) and quiet M dwarf located $ 29.87\pm0.02 pc$ away with a mass of $0.502\pm0.015 M_{\odot}$ and a radius of $0.501\pm0.030 R_{\odot}$. We validate the planet by combining space and ground-based photometry, spectroscopy, and imaging. The planet has a radius of $1.088\pm 0.064 R_{\oplus}$ with a $3 σ$ mass upper-limit of $8.7 M_{\oplus}$. We analyze the galactic kinematics and orbit of the host star LHS1815 and find that it has a large probability ($P_{thick}/P_{thin} = 6482$) to be in the thick disk with a much higher expected maximal height ($Z_{max} = 1.8 kpc$) above the Galactic plane compared with other TESS planet host stars. Future studies of the interior structure and atmospheric properties of planets in such systems using for example the upcoming James Webb Space Telescope (JWST), can investigate the differences in formation efficiency and evolution for planetary systems between different Galactic components (thick and thin disks, and halo).

preprint2020arXiv

Long-Period High-Amplitude Red Variables in the KELT Survey

We present a sample of 4,132 Mira-like variables (red variables with long periods and high amplitudes) in the Kilodegree Extremely Little Telescope (KELT) survey. Of these, 814 are new detections. We used 2MASS colors to identify candidate asymptotic giant branch (AGB) stars. We tested for photometric variability among the sample and used Lomb-Scargle to determine the periodicity of the variable sample. We selected variables with high amplitudes and strong periodic behavior using a Random Forest classifier. Of the sample of 4,132 Mira-like variables, we estimate that 70% are Miras, and 30% are semi-regular (SR) variables. We also adopt the method of using (W_{RP} - W_{K_s}) vs. (J - K_s) colors (Lebzelter et al. 2018) in distinguishing between O-rich and C-rich Miras and find it to be an improvement over 2MASS colors.

preprint2020arXiv

Orbital Refinement and Stellar Properties for the HD 9446, HD 43691, and HD 179079 Planetary Systems

The Transit Ephemeris Refinement and Monitoring Survey (TERMS) is a project which aims to detect transits of intermediate-long period planets by refining orbital parameters of the known radial velocity planets using additional data from ground based telescopes, calculating a revised transit ephemeris for the planet, then monitoring the planet host star during the predicted transit window. Here we present the results from three systems that had high probabilities of transiting planets: HD 9446 b & c, HD 43691 b, & HD 179079 b. We provide new radial velocity (RV) measurements that are then used to improve the orbital solution for the known planets. We search the RV data for indications of additional planets in orbit and find that HD 9446 shows a strong linear trend of 4.8$σ$. Using the newly refined planet orbital solutions, which include a new best-fit solution for the orbital period of HD 9446 c, and an improved transit ephemerides, we found no evidence of transiting planets in the photometry for each system. Transits of HD 9446 b can be ruled out completely and transits HD 9446 c & HD 43691 b can be ruled out for impact parameters up to b = 0.5778 and b = 0.898 respectively due to gaps in the photometry. A transit of HD 179079 b cannot be ruled out however due to the relatively small size of this planet compared to the large star and thus low signal to noise. We determine properties of the three host stars through spectroscopic analysis and find through photometric analysis that HD 9446 exhibits periodic variability.

preprint2020arXiv

Rotationally Driven Ultraviolet Emission of Red Giant Stars

Main sequence stars exhibit a clear rotation-activity relationship, in which rapidly rotating stars drive strong chromospheric/coronal ultraviolet and X-ray emission. While the vast majority of red giant stars are inactive, a few percent exhibit strong ultraviolet emission. Here we use a sample of 133 red giant stars observed by SDSS APOGEE and GALEX to demonstrate an empirical relationship between NUV excess and rotational velocity (vsini). Beyond this simple relationship, we find that NUV excess also correlates with rotation period and with Rossby number in a manner that shares broadly similar trends to those found in M dwarfs, including activity saturation among rapid rotators. Our data also suggest that the most extremely rapidly rotating giants may exhibit so-called "super-saturation", which could be caused by centrifugal stripping of these stars rotating at a high fraction of breakup speed. As an example application of our empirical rotation-activity relation, we demonstrate that the NUV emission observed from a recently reported system comprising a red giant with a black hole companion is fully consistent with arising from the rapidly rotating red giant in that system. Most fundamentally, our findings suggest a common origin of chromospheric activity in rotation and convection for cool stars from main sequence to red giant stages of evolution.

preprint2020arXiv

Stellar Characterization of M-dwarfs from the APOGEE Survey: A Calibrator Sample for the M-dwarf Metallicities

We present spectroscopic determinations of the effective temperatures, surface gravities and metallicities for 21 M-dwarfs observed at high-resolution (R $\sim$ 22,500) in the \textit{H}-band as part of the SDSS-IV APOGEE survey. The atmospheric parameters and metallicities are derived from spectral syntheses with 1-D LTE plane parallel MARCS models and the APOGEE atomic/molecular line list, together with up-to-date H$_{2}$O and FeH molecular line lists. Our sample range in $T_{\rm eff}$ from $\sim$ 3200 to 3800K, where eleven stars are in binary systems with a warmer (FGK) primary, while the other 10 M-dwarfs have interferometric radii in the literature. We define an $M_{K_{S}}$--Radius calibration based on our M-dwarf radii derived from the detailed analysis of APOGEE spectra and Gaia DR2 distances, as well as a mass-radius relation using the spectroscopically-derived surface gravities. A comparison of the derived radii with interferometric values from the literature finds that the spectroscopic radii are slightly offset towards smaller values, with $Δ$ = -0.01 $\pm$ 0.02 $R{\star}$/$R_{\odot}$. In addition, the derived M-dwarf masses based upon the radii and surface gravities tend to be slightly smaller (by $\sim$5-10\%) than masses derived for M-dwarf members of eclipsing binary systems for a given stellar radius. The metallicities derived for the 11 M-dwarfs in binary systems, compared to metallicities obtained for their hotter FGK main-sequence primary stars from the literature, shows excellent agreement, with a mean difference of [Fe/H](M-dwarf - FGK primary) = +0.04 $\pm$ 0.18 dex, confirming the APOGEE metallicity scale derived here for M-dwarfs.

preprint2020arXiv

Stellar Flares from the First Tess Data Release: Exploring a New Sample of M-dwarfs

We perform a study of stellar flares for the 24,809 stars observed with 2 minute cadence during the first two months of the TESS mission. Flares may erode exoplanets' atmospheres and impact their habitability, but might also trigger the genesis of life around small stars. TESS provides a new sample of bright dwarf stars in our galactic neighborhood, collecting data for thousands of M-dwarfs that might host habitable exoplanets. Here, we use an automated search for flares accompanied by visual inspection. Then, our public allesfitter code robustly selects the appropriate model for potentially complex flares via Bayesian evidence. We identify 1228 flaring stars, 673 of which are M-dwarfs. Among 8695 flares in total, the largest superflare increased the stellar brightness by a factor of 16.1. Bolometric flare energies range from 10^31.0 to 10^36.9 erg, with a median of 10^33.1 erg. Furthermore, we study the flare rate and energy as a function of stellar type and rotation period. We solidify past findings that fast rotating M-dwarfs are the most likely to flare, and that their flare amplitude is independent of the rotation period. Finally, we link our results to criteria for prebiotic chemistry, atmospheric loss through coronal mass ejections, and ozone sterilization. Four of our flaring M dwarfs host exoplanet candidates alerted on by TESS, for which we discuss how these effects can impact life. With upcoming TESS data releases, our flare analysis can be expanded to almost all bright small stars, aiding in defining criteria for exoplanet habitability.

preprint2020arXiv

Stellar parameters for the First Release of the MaStar Library: An Empirical Approach

We report the stellar atmospheric parameters for 7503 spectra contained in the first release of the MaNGA stellar library (MaStar) in SDSS DR15. The first release of MaStar contains 8646 spectra measured from 3321 unique stars, each covering the wavelength range 3622 Å to 10354 Å with a resolving power of $R \sim$ 1800. In this work, we first determined the basic stellar parameters: effective temperature ($\rm T_{eff}$), surface gravity ($\log g$), and metallicity ($\rm[Fe/H]$), which best fit the data using an empirical interpolator based on the Medium-resolution Isaac Newton Telescope library of empirical spectra (MILES), as implemented by the University of Lyon Spectroscopic analysis Software (Koleva et al. 2008, ULySS) package. While we analyzed all 8646 spectra from the first release of MaStar, since MaStar has a wider parameter-space coverage than MILES, not all of these fits are robust. In addition, not all parameter regions covered by MILES yield robust results, likely due to the non-uniform coverage of the parameter space by MILES. We tested the robustness of the method using the MILES spectra itself and identified a proxy based on the local density of the training set. With this proxy, we identified 7503 MaStar spectra with robust fitting results. They cover the range from 3179K to 20,517K in effective temperature ($\rm T_{eff}$), from 0.40 to 5.0 in surface gravity ($\log g$), and from $-$2.49 to $+$0.73 in metallicity ($\rm[Fe/H]$).

preprint2020arXiv

TESS Asteroseismic Analysis of the Known Exoplanet Host Star HD 222076

The Transiting Exoplanet Survey Satellite (TESS) is an all-sky survey mission aiming to search for exoplanets that transit bright stars. The high-quality photometric data of TESS are excellent for the asteroseismic study of solar-like stars. In this work, we present an asteroseismic analysis of the red-giant star HD~222076 hosting a long-period (2.4 yr) giant planet discovered through radial velocities. Solar-like oscillations of HD~222076 are detected around $203 \, μ$Hz by TESS for the first time. Asteroseismic modeling, using global asteroseismic parameters as input, yields a determination of the stellar mass ($M_\star = 1.12 \pm 0.12\, M_\odot$), radius ($R_\star = 4.34 \pm 0.21\,R_\odot$), and age ($7.4 \pm 2.7\,$Gyr), with precisions greatly improved from previous studies. The period spacing of the dipolar mixed modes extracted from the observed power spectrum reveals that the star is on the red-giant branch burning hydrogen in a shell surrounding the core. We find that the planet will not escape the tidal pull of the star and be engulfed into it within about $800\,$Myr, before the tip of the red-giant branch is reached.

preprint2020arXiv

TESS spots a hot Jupiter with an inner transiting Neptune

Hot Jupiters are rarely accompanied by other planets within a factor of a few in orbital distance. Previously, only two such systems have been found. Here, we report the discovery of a third system using data from the Transiting Exoplanet Survey Satellite (TESS). The host star, TOI-1130, is an 11th magnitude K-dwarf in the Gaia G band. It has two transiting planets: a Neptune-sized planet ($3.65\pm 0.10$ $R_E$) with a 4.1-day period, and a hot Jupiter ($1.50^{+0.27}_{-0.22}$ $R_J$) with an 8.4-day period. Precise radial-velocity observations show that the mass of the hot Jupiter is $0.974^{+0.043}_{-0.044}$ $M_J$. For the inner Neptune, the data provide only an upper limit on the mass of 0.17 $M_J$ (3$σ$). Nevertheless, we are confident the inner planet is real, based on follow-up ground-based photometry and adaptive optics imaging that rule out other plausible sources of the TESS transit signal. The unusual planetary architecture of and the brightness of the host star make TOI-1130 a good test case for planet formation theories, and an attractive target for future spectroscopic observations.

preprint2020arXiv

The Evolution of Rotation and Magnetic Activity in 94 Aqr Aa from Asteroseismology with TESS

Most previous efforts to calibrate how rotation and magnetic activity depend on stellar age and mass have relied on observations of clusters, where isochrones from stellar evolution models are used to determine the properties of the ensemble. Asteroseismology employs similar models to measure the properties of an individual star by matching its normal modes of oscillation, yielding the stellar age and mass with high precision. We use 27 days of photometry from the Transiting Exoplanet Survey Satellite to characterize solar-like oscillations in the G8 subgiant of the 94 Aqr triple system. The resulting stellar properties, when combined with a reanalysis of 35 yr of activity measurements from the Mount Wilson HK project, allow us to probe the evolution of rotation and magnetic activity in the system. The asteroseismic age of the subgiant agrees with a stellar isochrone fit, but the rotation period is much shorter than expected from standard models of angular momentum evolution. We conclude that weakened magnetic braking may be needed to reproduce the stellar properties, and that evolved subgiants in the hydrogen shell-burning phase can reinvigorate large-scale dynamo action and briefly sustain magnetic activity cycles before ascending the red giant branch.

preprint2020arXiv

The First Habitable Zone Earth-sized Planet from TESS. I: Validation of the TOI-700 System

We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

preprint2020arXiv

The Sixteenth Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

preprint2020arXiv

The TESS-Keck Survey I: A Warm Sub-Saturn-mass Planet and a Caution about Stray Light in TESS Cameras

We report the detection of a Saturn-size exoplanet orbiting HD 332231 (TOI 1456) in light curves from the Transiting Exoplanet Survey Satellite (TESS). HD 332231, an F8 dwarf star with a V-band magnitude of 8.56, was observed by TESS in Sectors 14 and 15. We detect a single-transit event in the Sector 15 presearch data conditioning (PDC) light curve. We obtain spectroscopic follow-up observations of HD 332231 with the Automated Planet Finder, Keck I, and SONG telescopes. The orbital period we infer from the radial velocity (RV) observations leads to the discovery of another transit in Sector 14 that was masked by PDC due to scattered light contamination. A joint analysis of the transit and RV data confirms the planetary nature of HD 332231 b, a Saturn-size ($0.867^{+0.027}_{-0.025} \; R_{\rm J}$), sub-Saturn-mass ($0.244\pm0.021 \; M_{\rm J}$) exoplanet on a 18.71 day circular orbit. The low surface gravity of HD 332231 b and the relatively low stellar flux it receives make it a compelling target for transmission spectroscopy. Also, the stellar obliquity is likely measurable via the Rossiter-McLaughlin effect, an exciting prospect given the 0.14 au orbital separation of HD 332231 b. The spectroscopic observations do not provide substantial evidence for any additional planets in the HD 332231 system, but continued RV monitoring is needed to further characterize this system. We also predict that the frequency and duration of masked data in the PDC light curves for TESS Sectors 14-16 could hide transits of some exoplanets with orbital periods between 10.5 and 17.5 days.

preprint2020arXiv

The TESS-Keck Survey II: An Ultra-Short Period Rocky Planet and its Siblings Transiting the Galactic Thick-Disk Star TOI-561

We report the discovery of TOI-561, a multi-planet system in the galactic thick disk that contains a rocky, ultra-short period planet (USP). This bright ($V=10.2$) star hosts three small transiting planets identified in photometry from the NASA TESS mission: TOI-561 b (TOI-561.02, P=0.44 days, $R_b = 1.45\pm0.11\,R_\oplus$), c (TOI-561.01, P=10.8 days, $R_c=2.90\pm0.13\,R_\oplus$), and d (TOI-561.03, P=16.3 days, $R_d=2.32\pm0.16\,R_\oplus$). The star is chemically ([Fe/H]$=-0.41\pm0.05$, [$α$/H]$=+0.23\pm0.05$) and kinematically consistent with the galactic thick disk population, making TOI-561 one of the oldest ($10\pm3\,$Gyr) and most metal-poor planetary systems discovered yet. We dynamically confirm planets b and c with radial velocities from the W. M. Keck Observatory High Resolution Echelle Spectrometer. Planet b has a mass and density of $3.2\pm0.8\,M_\oplus$ and $5.5^{+2.0}_{-1.6}\,$g$\,$cm$^{-3}$, consistent with a rocky composition. Its lower-than-average density is consistent with an iron-poor composition, although an Earth-like iron-to-silicates ratio is not ruled out. Planet c is $7.0\pm2.3\,M_\oplus$ and $1.6\pm0.6\,$g$\,$cm$^{-3}$, consistent with an interior rocky core overlaid with a low-mass volatile envelope. Several attributes of the photometry for planet d (which we did not detect dynamically) complicate the analysis, but we vet the planet with high-contrast imaging, ground-based photometric follow-up and radial velocities. TOI-561 b is the first rocky world around a galactic thick-disk star confirmed with radial velocities and one of the best rocky planets for thermal emission studies.

preprint2020arXiv

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

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

preprint2020arXiv

TOI 694 b and TIC 220568520 b: Two Low-Mass Companions Near the Hydrogen Burning Mass Limit Orbiting Sun-like Stars

We report the discovery of TOI 694 b and TIC 220568520 b, two low-mass stellar companions in eccentric orbits around metal-rich Sun-like stars, first detected by the Transiting Exoplanet Survey Satellite (TESS). TOI 694 b has an orbital period of 48.05131$\pm$0.00019 days and eccentricity of 0.51946$\pm$0.00081, and we derive a mass of 89.0$\pm$5.3 $M_J$ (0.0849$\pm$0.0051 $M_\odot$) and radius of 1.111$\pm$0.017 $R_J$ (0.1142$\pm$0.0017 $R_\odot$). TIC 220568520 b has an orbital period of 18.55769$\pm$0.00039 days and eccentricity of 0.0964$\pm$0.0032, and we derive a mass of 107.2$\pm$5.2 $M_J$ (0.1023$\pm$0.0050 $M_\odot$) and radius of 1.248$\pm$0.018 $R_J$ (0.1282$\pm$0.0019 $R_\odot$). Both binary companions lie close to and above the Hydrogen burning mass threshold that separates brown dwarfs and the lowest mass stars, with TOI 694 b being 2-$σ$ above the canonical mass threshold of 0.075 $M_\odot$. The relatively long periods of the systems mean that the magnetic fields of the low-mass companions are not expected to inhibit convection and inflate the radius, which according to one leading theory is common in similar objects residing in short-period tidally-synchronized binary systems. Indeed we do not find radius inflation for these two objects when compared to theoretical isochrones. These two new objects add to the short but growing list of low-mass stars with well-measured masses and radii, and highlight the potential of the TESS mission for detecting such rare objects orbiting bright stars.

preprint2020arXiv

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

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

preprint2020arXiv

TOI-503: The first known brown dwarf-Am star binary from the TESS mission

We report the discovery of an intermediate-mass transiting brown dwarf, TOI-503b, from the TESS mission. TOI-503b is the first brown dwarf discovered by TESS and orbits a metallic-line A-type star with a period of $P=3.6772 \pm 0.0001$ days. The light curve from TESS indicates that TOI-503b transits its host star in a grazing manner, which limits the precision with which we measure the brown dwarf's radius ($R_b = 1.34^{+0.26}_{-0.15} R_J$). We obtained high-resolution spectroscopic observations with the FIES, Ondřejov, PARAS, Tautenburg, and TRES spectrographs and measured the mass of TOI-503b to be $M_b = 53.7 \pm 1.2 M_J$. The host star has a mass of $M_\star = 1.80 \pm 0.06 M_\odot$, a radius of $R_\star = 1.70 \pm 0.05 R_\odot$, an effective temperature of $T_{\rm eff} = 7650 \pm 160$K, and a relatively high metallicity of $0.61\pm 0.07$ dex. We used stellar isochrones to derive the age of the system to be $\sim$180 Myr, which places its age between that of RIK 72b (a $\sim$10 Myr old brown dwarf in the Upper Scorpius stellar association) and AD 3116b (a $\sim$600 Myr old brown dwarf in the Praesepe cluster). We argue that this brown dwarf formed in-situ, based on the young age of the system and the long circularization timescale for this brown dwarf around its host star. TOI-503b joins a growing number of known short-period, intermediate-mass brown dwarfs orbiting main sequence stars, and is the second such brown dwarf known to transit an A star, after HATS-70b. With the growth in the population in this regime, the driest region in the brown dwarf desert ($35-55 M_J \sin{i}$) is reforesting and its mass range shrinking.

preprint2020arXiv

TOI-824 b: A New Planet on the Lower Edge of the Hot Neptune Desert

We report the detection of a transiting hot Neptune exoplanet orbiting TOI-824 (SCR J1448-5735), a nearby (d = 64 pc) K4V star, using data from the \textit{Transiting Exoplanet Survey Satellite} (TESS). The newly discovered planet has a radius, $R_{\rm{p}}$ = 2.93 $\pm$ 0.20 R$_{\oplus}$, and an orbital period of 1.393 days. Radial velocity measurements using the Planet Finder Spectrograph (PFS) and the High Accuracy Radial velocity Planet Searcher (HARPS) spectrograph confirm the existence of the planet and we estimate its mass to be $M_{\rm{p}}$ = 18.47 $\pm$ 1.84 M$_{\oplus}$. The planet's mean density is $ρ_{\rm{p}}$ = 4.03$^{+0.98}_{-0.78}$ g cm$^{-3}$ making it more than twice as dense as Neptune. TOI-824 b's high equilibrium temperature makes the planet likely to have a cloud free atmosphere, and thus an excellent candidate for follow up atmospheric studies. The detectability of TOI-824 b's atmosphere from both ground and space is promising and could lead to the detailed characterization of the most irradiated, small planet at the edge of the hot Neptune desert that has retained its atmosphere to date.

preprint2020arXiv

Two intermediate-mass transiting brown dwarfs from the TESS mission

We report the discovery of two intermediate-mass brown dwarfs (BDs), TOI-569b and TOI-1406b, from NASA's Transiting Exoplanet Survey Satellite mission. TOI-569b has an orbital period of $P = 6.55604 \pm 0.00016$ days, a mass of $M_b = 64.1 \pm 1.9 M_J$, and a radius of $R_b = 0.75 \pm 0.02 R_J$. Its host star, TOI-569, has a mass of $M_\star = 1.21 \pm 0.03 M_\odot$, a radius of $R_\star = 1.47 \pm 0.03 R_\odot$, $\rm [Fe/H] = +0.29 \pm 0.09$ dex, and an effective temperature of $T_{\rm eff} = 5768 \pm 110K$. TOI-1406b has an orbital period of $P = 10.57415 \pm 0.00063$ days, a mass of $M_b =46.0 \pm 2.7 M_J$, and a radius of $R_b = 0.86 \pm 0.03 R_J$. The host star for this BD has a mass of $M_\star =1 .18 \pm 0.09 M_\odot$, a radius of $R_\star = 1.35 \pm 0.03 R_\odot$, $ \rm [Fe/H] = -0.08 \pm 0.09$ dex and an effective temperature of $T_{\rm eff} = 6290 \pm 100K$. Both BDs are in circular orbits around their host stars and are older than 3 Gyr based on stellar isochrone models of the stars. TOI-569 is one of two slightly evolved stars known to host a transiting BD (the other being KOI-415). TOI-1406b is one of three known transiting BDs to occupy the mass range of $40-50 M_J$ and one of two to have a circular orbit at a period near 10 days (with the first being KOI-205b).Both BDs have reliable ages from stellar isochrones in addition to their well-constrained masses and radii, making them particularly valuable as tests for substellar isochrones in the BD mass-radius diagram.

preprint2019arXiv

TESS Reveals HD 118203 b to be a Transiting Planet

The exoplanet HD 118203 b, orbiting a bright (V = 8.05) host star, was discovered using the radial velocity method by da Silva et al. (2006), but was not previously known to transit. TESS photometry has revealed that this planet transits its host star. Five planetary transits were observed by TESS, allowing us to measure the radius of the planet to be $1.133 \pm 0.031 R_J$, and to calculate the planet mass to be $2.173 \pm 0.078 M_J$. The host star is slightly evolved with an effective temperature of $T_{\rm eff} = 5692 \pm 83$ K and a surface gravity of ${\rm log}(g) = 3.891 \pm 0.019$. With an orbital period of $6.134980 \pm 0.000038$ days and an eccentricity of $0.316 \pm 0.021$, the planet occupies a transitional regime between circularized hot Jupiters and more dynamically active planets at longer orbital periods. The host star is among the ten brightest known to have transiting giant planets, providing opportunities for both planetary atmospheric and asteroseismic studies.

preprint2019arXiv

TOI-132 b: A short-period planet in the Neptune desert transiting a $V=11.3$ G-type star

The Neptune desert is a feature seen in the radius-mass-period plane, whereby a notable dearth of short period, Neptune-like planets is found. Here we report the {\it TESS} discovery of a new short-period planet in the Neptune desert, orbiting the G-type dwarf TYC\,8003-1117-1 (TOI-132). {\it TESS} photometry shows transit-like dips at the level of $\sim$1400 ppm occurring every $\sim$2.11 days. High-precision radial velocity follow-up with HARPS confirmed the planetary nature of the transit signal and provided a semi-amplitude radial velocity variation of $\sim$11.5 m s$^{-1}$, which, when combined with the stellar mass of $0.97\pm0.06$ $M_{\odot}$, provides a planetary mass of 22.83$^{+1.81}_{-1.80}$ $M_{\oplus}$. Modeling the {\it TESS} high-quality light curve returns a planet radius of 3.43$^{+0.13}_{-0.14}$ $R_{\oplus}$, and therefore the planet bulk density is found to be 3.11$^{+0.44}_{-0.450}$ g cm$^{-3}$. Planet structure models suggest that the bulk of the planet mass is in the form of a rocky core, with an atmospheric mass fraction of 4.3$^{+1.2}_{-2.3}$\%. TOI-132 b is a {\it TESS} Level 1 Science Requirement candidate, and therefore priority follow-up will allow the search for additional planets in the system, whilst helping to constrain low-mass planet formation and evolution models, particularly valuable for better understanding the Neptune desert.

preprint2016arXiv

A statistical analysis of the accuracy of the digitized magnitudes of photometric plates on the time scale of decades with an application to the century-long light curve of KIC 8462852

We present a statistical analysis of the accuracy of the digitized magnitudes of photometric plates on the time scale of decades. In our examination of archival Johnson B photometry from the Harvard DASCH archive, we find a median RMS scatter of lightcurves of order 0.15mag over the range B~9-17 for all calibrations. Slight underlying systematics (trends or flux discontinuities) are on a level of ~<=0.2mag per century (1889-1990) for the majority of constant stars. These historic data can be unambiguously used for processes that happen on scales of magnitudes, and need to be carefully examined in cases approaching the noise floor. The characterization of these limits in photometric stability may guide future studies in their use of plate archives. We explain these limitations for the example case of KIC8462852, which has been claimed to dim by 0.16mag per century, and show that this trend cannot be considered as significant.

preprint2016arXiv

An Extreme Analogue of $ε$ Aurigae: An M-giant Eclipsed Every 69 Years by a Large Opaque Disk Surrounding a Small Hot Source

We present TYC 2505-672-1 as a newly discovered and remarkable eclipsing system comprising an M-type red giant that undergoes a ~3.45 year long, near-total eclipse (depth of ~4.5 mag) with a very long period of ~69.1 yr. TYC 2505-672-1 is now the longest-period eclipsing binary system yet discovered, more than twice as long as that of the currently longest-period system, $ε$ Aurigae. We show from analysis of the light curve including both our own data and historical data spanning more than 120 yr and from modeling of the spectral energy distribution, both before and during eclipse, that the red giant primary is orbited by a moderately hot source (T$_{eff}$~8000 K) that is itself surrounded by an extended, opaque circumstellar disk. From the measured ratio of luminosities, the radius of the hot companion must be in the range 0.1-0.5 Rsun (depending on the assumed radius of the red giant primary), which is an order of magnitude smaller than that for a main sequence A star and 1-2 orders of magnitude larger than that for a white dwarf. The companion is therefore most likely a "stripped red giant" subdwarf-B type star destined to become a He white dwarf. It is however somewhat cooler than most sdB stars, implying a very low mass for this "pre-He-WD" star. The opaque disk surrounding this hot source may be a remnant of the stripping of its former hydrogen envelope. However, it is puzzling how this object became stripped, given that it is at present so distant (orbital semi-major axis of ~24 AU) from the current red giant primary star. Extrapolating from our calculated ephemeris, the next eclipse should begin in early UT 2080 April and end in mid UT 2083 September (eclipse center UT 2081 December 24). This system is poised to become an exemplar of a very rare class of systems, even more extreme in several respects than the well studied archetype $ε$ Aurigae.

preprint2016arXiv

DM Ori: A Young Star Occulted by a Disturbance in its Protoplanetary Disk

In some planet formation theories, protoplanets grow gravitationally within a young star's protoplanetary disk, a signature of which may be a localized disturbance in the disk's radial and/or vertical structure. Using time-series photometric observations by the Kilodegree Extremely Little Telescope South (KELT-South) project and the All-Sky Automated Survey for SuperNovae (ASAS-SN), combined with archival observations, we present the discovery of two extended dimming events of the young star, DM Ori. This young system faded by $\sim$1.5 mag from 2000 March to 2002 August and then again in 2013 January until 2014 September (depth $\sim$1.7 mag). We constrain the duration of the 2000-2002 dimming to be $<$ 860 days, and the event in 2013-2014 to be $<$ 585 days, separated by $\sim$12.5 years. A model of the spectral energy distribution (SED) indicates a large infrared excess consistent with an extensive circumstellar disk. Using basic kinematic arguments, we propose that DM Ori is likely being periodically occulted by a feature (possibly a warp or perturbation) in its circumstellar disk. In this scenario, the occulting feature is located $>$6 AU from the host star, moving at $\sim$14.6 km/s, and is $\sim$4.9 AU in width. This localized structure may indicate a disturbance such as may be caused by a protoplanet early in its formation.

preprint2016arXiv

Eclipsing Binary Stars as Benchmarks for Trigonometric Parallaxes in the Gaia Era

We present fits to the broadband photometric spectral energy distributions (SEDs) of 158 eclipsing binaries (EBs) in the Tycho-2 catalog. These EBs were selected because they have highly precise stellar radii, effective temperatures, and in many cases metallicities previously determined in the literature, and thus have bolometric luminosities that are typically good to $\lesssim$ 10%. In most cases the available broadband photometry spans a wavelength range 0.4-10 $μ$m, and in many cases spans 0.15-22 $μ$m. The resulting SED fits, which have only extinction as a free parameter, provide a virtually model-independent measure of the bolometric flux at Earth. The SED fits are satisfactory for 156 of the EBs, for which we achieve typical precisions in the bolometric flux of $\approx$ 3%. Combined with the accurately known bolometric luminosity, the result for each EB is a predicted parallax that is typically precise to $\lesssim$ 5%. These predicted parallaxes---with typical uncertainties of 200 $μ$as---are 4-5 times more precise than those determined by Hipparcos for 99 of the EBs in our sample, with which we find excellent agreement. There is no evidence among this sample for significant systematics in the Hipparcos parallaxes of the sort that notoriously afflicted the Pleiades measurement. The EBs are distributed over the entire sky, span more than 10 mag in brightness, reach distances of more than 5 kpc, and in many cases our predicted parallaxes should also be more precise than those expected from the Gaia first data release. The EBs studied here can thus serve as empirical, independent benchmarks for these upcoming fundamental parallax measurements.

preprint2016arXiv

Evidence for a systematic offset of $-$0.25 mas in the Gaia DR1 parallaxes

We test the parallaxes reported in the Gaia first data release using the sample of eclipsing binaries with accurate, empirical distances from Stassun \& Torres (2016). We find an average offset of $-$0.25$\pm$0.05 mas in the sense of the Gaia parallaxes being too small (i.e., the distances too long). The offset does not depend strongly on obvious parameters such as color or brightness. However, we find with high confidence that the offset may depend on ecliptic latitude: the mean offset is $-$0.38$\pm$0.06 mas in the ecliptic north and $-$0.05$\pm$0.09 mas in the ecliptic south. The ecliptic latitude dependence may also be represented by the linear relation, $Δπ\approx -0.22(\pm0.05) -0.003(\pm0.001)\timesβ$ mas ($β$ in degrees). Finally, there is a possible dependence of the parallax offset on distance, with the offset becoming negligible for $π\lesssim 1$ mas; we discuss whether this could be caused by a systematic error in the eclipsing binary distance scale, and reject this interpretation as unlikely.

preprint2016arXiv

Examining the relationships between colour, $T_{\rm eff}$, and [M/H] for APOGEE K and M dwarfs

We present the effective temperatures ($T_{\rm eff}$), metallicities, and colours in SDSS, 2MASS, and WISE filters, of a sample of 3834 late-K and early-M dwarfs selected from the Sloan Digital Sky Survey APOGEE spectroscopic survey ASPCAP catalog. We confirm that ASPCAP $T_{\rm eff}$ values between 3550 K$<T_{\rm eff}<$4200 K are accurate to $\sim$100 K compared to interferometric $T_{\rm eff}$ values. In that same $T_{\rm eff}$ range, ASPCAP metallicities are accurate to 0.18 dex between $-1.0<$[M/H]$<0.2$. For these cool dwarfs, nearly every colour is sensitive to both $T_{\rm eff}$ and metallicity. Notably, we find that $g-r$ is not a good indicator of metallicity for near-solar metallicity early-M dwarfs. We confirm that $J-K_S$ colour is strongly dependent on metallicity, and find that $W1-W2$ colour is a promising metallicity indicator. Comparison of the late-K and early-M dwarf colours, metallicities, and $T_{\rm eff}$ to those from three different model grids shows reasonable agreement in $r-z$ and $J-K_S$ colours, but poor agreement in $u-g$, $g-r$, and $W1-W2$. Comparison of the metallicities of the KM dwarf sample to those from previous colour-metallicity relations reveals a lack of consensus in photometric metallicity indicators for late-K and early-M dwarfs. We also present empirical relations for $T_{\rm eff}$ as a function of $r-z$ colour combined with either [M/H] or $W1-W2$ colour, and for [M/H] as a function of $r-z$ and $W1-W2$ colour. These relations yield $T_{\rm eff}$ to $\sim$100 K and [M/H] to $\sim$0.18 dex precision with colours alone, for $T_{\rm eff}$ in the range of 3550-4200 K and [M/H] in the range of $-$0.5-0.2.

preprint2016arXiv

Integration of a $^6$LiInSe$_2$ Thermal Neutron Detector into a CubeSat Instrument

We present a preliminary design for a novel neutron detection system that is compact, lightweight, and low power consuming, utilizing the CubeSat platform making it suitable for space-based applications. This is made possible using the scintillating crystal lithium indium diselenide ($^6$LiInSe$_2$), the first crystal to include $^6$Li in the crystalline structure, and a silicon avalanche photodiode (Si-APD). The schematics of this instrument are presented as well as the response of the instrument to initial testing under alpha radiation. A principal aim of this work is to demonstrate the feasibility of such a neutron detection system within a CubeSat platform. The entire end-to-end system presented here is 10 cm x 10 cm x 15 cm, weighs 670 grams and requires 5 Volts direct current at 3 Watts.

preprint2016arXiv

KELT-10b: The First Transiting Exoplanet from the KELT-South Survey -- A Hot Sub-Jupiter Transiting a V = 10.7 Early G-Star

We report the discovery of KELT-10b, the first transiting exoplanet discovered using the KELT-South telescope. KELT-10b is a highly inflated sub-Jupiter mass planet transiting a relatively bright $V = 10.7$ star (TYC 8378-64-1), with T$_{eff}$ = $5948\pm74$ K, $\log{g}$ = $4.319_{-0.030}^{+0.020}$ and [Fe/H] = $0.09_{-0.10}^{+0.11}$, an inferred mass M$_{*}$ = $1.112_{-0.061}^{+0.055}$ M$_{\odot}$ and radius R$_{*}$ = $1.209_{-0.035}^{+0.047}$ R$_{\odot}$. The planet has a radius R$_{P}$ = $1.399_{-0.049}^{+0.069}$ R$_{J}$ and mass M$_{P}$ = $0.679_{-0.038}^{+0.039}$ M$_{J}$. The planet has an eccentricity consistent with zero and a semi-major axis $a$ = $0.05250_{-0.00097}^{+0.00086}$ AU. The best fitting linear ephemeris is $T_{0}$ = 2457066.72045$\pm$0.00027 BJD$_{TDB}$ and P = 4.1662739$\pm$0.0000063 days. This planet joins a group of highly inflated transiting exoplanets with a radius much larger and a mass much less than those of Jupiter. The planet, which boasts deep transits of 1.4%, has a relatively high equilibrium temperature of T$_{eq}$ = $1377_{-23}^{+28}$ K, assuming zero albedo and perfect heat redistribution. KELT-10b receives an estimated insolation of $0.817_{-0.054}^{+0.068}$ $\times$ 10$^9$ erg s$^{-1}$ cm$^{-2}$, which places it far above the insolation threshold above which hot Jupiters exhibit increasing amounts of radius inflation. Evolutionary analysis of the host star suggests that KELT-10b is unlikely to survive beyond the current subgiant phase, due to a concomitant in-spiral of the planet over the next $\sim$1 Gyr. The planet transits a relatively bright star and exhibits the third largest transit depth of all transiting exoplanets with V $<$ 11 in the southern hemisphere, making it a promising candidate for future atmospheric characterization studies.

preprint2016arXiv

KELT-14b and KELT-15b: An Independent Discovery of WASP-122b and a New Hot Jupiter

We report the discovery of KELT-14b and KELT-15b, two hot Jupiters from the KELT-South survey. KELT-14b, an independent discovery of the recently announced WASP-122b, is an inflated Jupiter mass planet that orbits a $\sim5.0^{+0.3}_{-0.7}$ Gyr, $V$ = 11.0, G2 star that is near the main sequence turnoff. The host star, KELT-14 (TYC 7638-981-1), has an inferred mass $M_{*}$=$1.18_{-0.07}^{+0.05}$$M_{\odot}$ and radius $R_{*}$=$1.37\pm{-0.08}$$R_{\odot}$, and has $T_{eff}$=$5802_{-92}^{+95}$K, $\log{g_*}$=$4.23_{-0.04}^{+0.05}$ and =$0.33\pm{0.09}$. The planet orbits with a period of $1.7100588 \pm 0.0000025$ days ($T_{0}$=2457091.02863$\pm$0.00047) and has a radius R$_{p}$=$1.52_{-0.11}^{+0.12}$$R_{J}$ and mass M$_{p}$=$1.196\pm0.072$$M_{J}$, and the eccentricity is consistent with zero. KELT-15b is another inflated Jupiter mass planet that orbits a $\sim$ $4.6^{+0.5}_{-0.4}$ Gyr, $V$ = 11.2, G0 star (TYC 8146-86-1) that is near the "blue hook" stage of evolution prior to the Hertzsprung gap, and has an inferred mass $M_{*}$=$1.181_{-0.050}^{+0.051}$$M_{\odot}$ and radius $R_{*}$=$1.48_{-0.04}^{+0.09}$$R_{\odot}$, and $T_{eff}$=$6003_{-52}^{+56}$K, $\log{g_*}$=$4.17_{-0.04}^{+0.02}$ and [Fe/H]=$0.05\pm0.03$. The planet orbits on a period of $3.329441 \pm 0.000016$ days ($T_{0}$ = 2457029.1663$\pm$0.0073) and has a radius R$_{p}$=$1.443_{-0.057}^{+0.11}$$R_{J}$ and mass M$_{p}$=$0.91_{-0.22}^{+0.21}$$M_{J}$ and an eccentricity consistent with zero. KELT-14b has the second largest expected emission signal in the K-band for known transiting planets brighter than $K<10.5$. Both KELT-14b and KELT-15b are predicted to have large enough emission signals that their secondary eclipses should be detectable using ground-based observatories.

preprint2016arXiv

KELT-17b: A hot-Jupiter transiting an A-star in a misaligned orbit detected with Doppler tomography

We present the discovery of a hot-Jupiter transiting the V=9.23 mag main-sequence A-star KELT-17 (BD+14 1881). KELT-17b is a 1.31 -0.29/+0.28 Mj, 1.525 -0.060/+0.065 Rj hot-Jupiter in a 3.08 day period orbit misaligned at -115.9 +/- 4.1 deg to the rotation axis of the star. The planet is confirmed via both the detection of the radial velocity orbit, and the Doppler tomographic detection of the shadow of the planet over two transits. The nature of the spin-orbit misaligned transit geometry allows us to place a constraint on the level of differential rotation in the host star; we find that KELT-17 is consistent with both rigid-body rotation and solar differential rotation rates (alpha < 0.30 at 2 sigma significance). KELT-17 is only the fourth A-star with a confirmed transiting planet, and with a mass of 1.635 -0.061/+0.066 Msun, effective temperature of 7454 +/- 49 K, and projected rotational velocity v sin I_* = 44.2 -1.3/+1.5 km/s; it is amongst the most massive, hottest, and most rapidly rotating of known planet hosts.

preprint2016arXiv

New Pleiades Eclipsing Binaries and a Hyades Transiting System Identified by K2

We present the discovery in Kepler's $K2$ mission observations and our follow-up radial velocity observations from Keck/HIRES for four eclipsing binary (EB) star systems in the young benchmark Pleiades cluster. Based on our modeling results, we announce two new low mass ($M_{tot} < 0.6 M_\odot$) EBs among Pleiades members (HCG 76 and MHO 9) and we report on two previously known Pleiades binaries that are also found to be EB systems (HII 2407 and HD 23642). We measured the masses of the binary HCG 76 to $\lesssim$2.5% precision, and the radii to $\lesssim$4.5% precision, which together with the precise effective temperatures yield an independent Pleiades distance of 132$\pm$5 pc. We discuss another EB towards the Pleiades that is a possible but unlikely Pleiades cluster member (AK II 465). The two new confirmed Pleiades systems extend the mass range of Pleiades EB components to 0.2-2 $M_\odot$. Our initial measurements of the fundamental stellar parameters for the Pleiades EBs are discussed in the context of the current stellar models and the nominal cluster isochrone, finding good agreement with the stellar models of Baraffe et al. (2015) at the nominal Pleiades age of 120 Myr. Finally, in the Hyades, we report a new low mass eclipsing system (vA 50) that was concurrently discovered and studied by Mann et al. (2016). We confirm that the eclipse is likely caused by a Neptune-sized transiting planet, and with the additional radial velocity constraints presented here we improve the constraint on the maximum mass of the planet to be $\lesssim$1.2 $M_\mathrm{Jup}$.

preprint2016arXiv

The K2-TESS Stellar Properties Catalog

We introduce a catalog of stellar properties for stars observed by the Kepler follow-on mission, K2. We base the catalog on a cross-match between the K2 Campaign target lists and the current working version of the NASA TESS target catalog. The resulting K2-TESS Stellar Properties Catalog includes value-added information from the TESS Target Catalog, including stellar colors, proper motions, effective temperatures, an estimated luminosity class (dwarf/subgiant versus giant) for each star based on reduced-proper-motion, and many other properties via cross-matching to other all-sky catalogs. Also included is the Guest Observer program identification number(s) associated with each K2 target. The K2-TESS Stellar Properties Catalog is available to the community as a freely accessible data portal on the Filtergraph system at: http://filtergraph.vanderbilt.edu/tess_k2campaigns .

preprint2016arXiv

The Stability of F-star Brightness on Century Timescales

The century-long photometric record of the DASCH project provides a unique window into the variability of stars normally considered to be photometrically inactive. In this paper, we look for long-term trends in the brightness of F stars, with particular attention to KIC 8462852,an F3 main sequence star that has been identified as significant short-term variability according to Kepler observations. Although a simple search for variability suggests long-term dimming of a number of F stars, we find that such trends are artifacts of the 'Menzel Gap' in the DASCH data. That includes the behavior of KIC 8462852, which we believe is consistent with constant flux over the full duration of observations. We do, however, present a selection of F stars that do have significant photometric trends, even after systematics are taken into account.

preprint2016arXiv

The Time-Domain Spectroscopic Survey: Understanding the Optically Variable Sky with SEQUELS in SDSS-III

The Time-Domain Spectroscopic Survey (TDSS) is an SDSS-IV eBOSS subproject primarily aimed at obtaining identification spectra of ~220,000 optically-variable objects systematically selected from SDSS/Pan-STARRS1 multi-epoch imaging. We present a preview of the science enabled by TDSS, based on TDSS spectra taken over ~320 deg^2 of sky as part of the SEQUELS survey in SDSS-III, which is in part a pilot survey for eBOSS in SDSS-IV. Using the 15,746 TDSS-selected single-epoch spectra of photometrically variable objects in SEQUELS, we determine the demographics of our variability-selected sample, and investigate the unique spectral characteristics inherent in samples selected by variability. We show that variability-based selection of quasars complements color-based selection by selecting additional redder quasars, and mitigates redshift biases to produce a smooth quasar redshift distribution over a wide range of redshifts. The resulting quasar sample contains systematically higher fractions of blazars and broad absorption line quasars than from color-selected samples. Similarly, we show that M-dwarfs in the TDSS-selected stellar sample have systematically higher chromospheric active fractions than the underlying M-dwarf population, based on their H-alpha emission. TDSS also contains a large number of RR Lyrae and eclipsing binary stars with main-sequence colors, including a few composite-spectrum binaries. Finally, our visual inspection of TDSS spectra uncovers a significant number of peculiar spectra, and we highlight a few cases of these interesting objects. With a factor of ~15 more spectra, the main TDSS survey in SDSS-IV will leverage the lessons learned from these early results for a variety of time-domain science applications.

preprint2016arXiv

VISION: A Six-Telescope Fiber-Fed Visible Light Beam Combiner for the Navy Precision Optical Interferometer

Visible-light long baseline interferometry holds the promise of advancing a number of important applications in fundamental astronomy, including the direct measurement of the angular diameters and oblateness of stars, and the direct measurement of the orbits of binary and multiple star systems. To advance, the field of visible-light interferometry requires development of instruments capable of combining light from 15 baselines (6 telescopes) simultaneously. The Visible Imaging System for Interferometric Observations at NPOI (VISION) is a new visible light beam combiner for the Navy Precision Optical Interferometer (NPOI) that uses single-mode fibers to coherently combine light from up to six telescopes simultaneously with an image-plane combination scheme. It features a photometric camera for calibrations and spatial filtering from single-mode fibers with two Andor Ixon electron multiplying CCDs. This paper presents the VISION system, results of laboratory tests, and results of commissioning on-sky observations. A new set of corrections have been determined for the power spectrum and bispectrum by taking into account non-Gaussian statistics and read noise present in electron-multipying CCDs to enable measurement of visibilities and closure phases in the VISION post-processing pipeline. The post-processing pipeline has been verified via new on-sky observations of the O-type supergiant binary $ζ$ Orionis A, obtaining a flux ratio of $2.18\pm0.13$ mag with a position angle of $223.9\pm1.0^{\circ}$ and separation $40.6\pm1.8$ mas over 570-750 nm, in good agreement with expectations from the previously published orbit.

preprint2015arXiv

A Disk-based Dynamical Mass Estimate for the Young Binary AK Sco

We present spatially and spectrally resolved Atacama Large Millimeter/submillimeter Array (ALMA) observations of gas and dust in the disk orbiting the pre-main sequence binary AK Sco. By forward-modeling the disk velocity field traced by CO J=2-1 line emission, we infer the mass of the central binary, $M_\ast = 2.49 \pm 0.10~M_\odot$, a new dynamical measurement that is independent of stellar evolutionary models. Assuming the disk and binary are co-planar within $\sim$2°, this disk-based binary mass measurement is in excellent agreement with constraints from radial velocity monitoring of the combined stellar spectra. These ALMA results are also compared with the standard approach of estimating masses from the location of the binary in the Hertzsprung-Russell diagram, using several common pre-main sequence model grids. These models predict stellar masses that are marginally consistent with our dynamical measurement (at $\sim 2\,σ$), but are systematically high (by $\sim$10%). These same models consistently predict an age of $18\pm1$ Myr for AK Sco, in line with its membership in the Upper Centaurus-Lupus association but surprisingly old for it to still host a gas-rich disk. As ALMA accumulates comparable data for large samples of pre-main sequence stars, the methodology employed here to extract a dynamical mass from the disk rotation curve should prove extraordinarily useful for efforts to characterize the fundamental parameters of early stellar evolution.

preprint2015arXiv

A Granulation "Flicker"-based Measure of Stellar Surface Gravity

In Bastien et al. (2013) we found that high quality light curves, such as those obtained by Kepler, may be used to measure stellar surface gravity via granulation-driven light curve "flicker". Here, we update and extend the relation originally presented in Bastien et al. (2013) after calibrating flicker against a more robust set of asteroseismically derived surface gravities. We describe in detail how we extract the flicker signal from the light curves, including how we treat phenomena, such as exoplanet transits and shot noise, that adversely affect the measurement of flicker. We examine the limitations of the technique, and, as a result, we now provide an updated treatment of the flicker-based logg error. We briefly highlight further applications of the technique, such as astrodensity profiling or its use in other types of stars with convective outer layers. We discuss potential uses in current and upcoming space-based photometric missions. Finally, we supply flicker-based logg values, and their uncertainties, for 27 628 Kepler stars not identified as transiting-planet hosts, with 4500<teff<7150 K, 2.5<logg<4.6, Kepler magnitude <13.5, and overall photometric amplitudes <10 parts per thousand.

preprint2015arXiv

Detailed Abundances of Planet-Hosting Wide Binaries. II. HD80606 + HD80607

We present a detailed chemical abundance analysis of 15 elements in the planet-hosting wide binary system HD80606 + HD80607 using Keck/HIRES spectra. As in our previous analysis of the planet-hosting wide binary HD20782 + HD20781, we presume that these two G5 dwarf stars formed together and therefore had identical primordial abundances. In this binary, HD80606 hosts an eccentric ($e\approx0.93$) giant planet at $\sim$0.5 AU, but HD80607 has no detected planets. If close-in giant planets on eccentric orbits are efficient at scattering rocky planetary material into their host stars, then HD80606 should show evidence of having accreted rocky material while HD80607 should not. Here we show that the trends of abundance versus element condensation temperature for HD80606 and HD80607 are statistically indistinguishable, corroborating the recent result of Saffe et al. This could suggest that both stars accreted similar amounts of rocky material; indeed, our model for the chemical signature of rocky planet accretion indicates that HD80606 could have accreted up to 2.5~$M_{\oplus}$ of rocky material---about half that contained in the Solar System and primordial asteroid belt---relative to HD80607 and still be consistent with the data. Since HD80607 has no known giant planets that might have pushed rocky planet material via migration onto that star, we consider it more likely that HD80606/07 experienced essentially no rocky planet accretion. This in turn suggests that the migration history of the HD80606 giant planet must have been such that it ejected any close-in planetary material that might have otherwise been shepherded onto the star.

preprint2015arXiv

First Results from the Disk Eclipse Search with KELT (DESK) Survey

Using time-series photometry from the Kilodegree Extremely Little Telescope (KELT) exoplanet survey, we are looking for eclipses of stars by their protoplanetary disks, specifically in young stellar associations. To date, we have discovered two previously unknown, large dimming events around the young stars RW Aurigae and V409 Tau. We attribute the dimming of RW Aurigae to an occultation by its tidally disrupted disk, with the disruption perhaps resulting from a recent flyby of its binary companion. Even with the dynamical environment of RW Aurigae, the distorted disk material remains very compact and presumably capable of forming planets. This system also shows that strong binary interactions with disks can also influence planet and core composition by stirring up and mixing materials during planet formation. We interpret the dimming of V409 Tau to be due to a feature, possibly a warp or perturbation, lying at least 10 AU from the host star in its nearly edge-on circumstellar disk.

preprint2015arXiv

HII 2407: A Low-Mass Eclipsing Binary Revealed by K2 Observations of the Pleiades

The star HII 2407 is a member of the relatively young Pleiades star cluster and was previously discovered to be a single-lined spectroscopic binary. It is newly identified here within $Kepler$/$K2$ photometric time series data as an eclipsing binary system. Mutual fitting of the radial velocity and photometric data leads to an orbital solution and constraints on fundamental stellar parameters. While the primary has arrived on the main sequence, the secondary is still pre-main-sequence and we compare our results for the $M/M_\odot$ and $R/R_\odot$ values with stellar evolutionary models. We also demonstrate that the system is likely to be tidally synchronized. Follow-up infrared spectroscopy is likely to reveal the lines of the secondary, allowing for dynamically measured masses and elevating the system to benchmark eclipsing binary status.

preprint2015arXiv

Impact of Declining Proposal Success Rates on Scientific Productivity

Over the last decade proposal success rates in the fundamental sciences have dropped significantly. Astronomy and related fields funded by NASA and NSF are no exception. Data across agencies show that this is not principally the result of a decline in proposal merit (the proportion of proposals receiving high rankings is largely unchanged), nor of a shift in proposer demographics (seniority, gender, and institutional affiliation have all remained unchanged), nor of an increase (beyond inflation) in the average requested funding per proposal, nor of an increase in the number of proposals per investigator in any one year. Rather, the statistics are consistent with a scenario in which agency budgets for competed research are flat or decreasing in inflation-adjusted dollars, the overall population of investigators has grown, and a larger proportion of these investigators are resubmitting meritorious but unfunded proposals. This White Paper presents statistics which support this conclusion, as well as recent research on the time cost of proposal writing versus that of producing publishable results. We conclude that an aspirational proposal success rate of 30-35% would still provide a healthily competitive environment for researchers, would more fully utilize the scientific capacity of the community's facilities and missions, and provide relief to the funding agencies who face the logistics of ever-increasing volumes of proposals.

preprint2015arXiv

IN-SYNC IV - The Young Stellar Population in the Orion A Molecular Cloud

We present the results of the SDSS APOGEE INfrared Spectroscopy of Young Nebulous Clusters program (IN-SYNC) survey of the Orion A molecular cloud. This survey obtained high resolution near infrared (NIR) spectroscopy of about 2700 young pre-main sequence stars throughout the region, acquired across five distinct fields spanning 6deg field of view (FOV). With these spectra, we have measured accurate stellar parameters (T_eff, log g, v sin i) and extinctions, and placed the sources in the Hertzsprung-Russel Diagram (HRD). We have also extracted radial velocities for the kinematic characterization of the population. We compare our measurements with literature results for a sub-sample of targets in order to assess the performances and accuracy of the survey. Source extinction shows evidence for dust grains that are larger than those in the diffuse interstellar medium (ISM): we estimate an average R_V=5.5 in the region. Importantly, we find a clear correlation between HRD inferred ages and spectroscopic surface-gravity inferred ages. This clearly indicates a real spread of stellar radii at fixed temperature, and together with additional correlations with extinction and with disk presence, strongly suggests a real spread of ages large than a few Myr. Focussing on the young population around NGC1980 iota Ori, which has previously been suggested to be a separate, foreground, older cluster, we confirm its older (5Myr) age and low A_V, but considering that its radial velocity distribution is indistinguishable from the Orion A's population, we suggest that NGC1980 is part of Orion A's star formation activity. Based on their stellar parameters and kinematic properties, we identify 383 new candidate members of Orion A, most of which are diskless sources in areas of the region poorly studied by previous works.

preprint2015arXiv

KELT-4Ab: An inflated Hot Jupiter transiting the bright (V~10) component of a hierarchical triple

We report the discovery of KELT-4Ab, an inflated, transiting Hot Jupiter orbiting the brightest component of a hierarchical triple stellar system. The host star is an F star with $T_{\rm eff}=6206\pm75$ K, $\log g=4.108\pm0.014$, $\left[{\rm Fe}/{\rm H}\right]=-0.116_{-0.069}^{+0.065}$, ${\rm M_*}=1.201_{-0.061}^{+0.067} \ {\rm M}_{\odot}$, and ${\rm R_*}=1.610_{-0.068}^{+0.078} \ {\rm R}_{\odot}$. The best-fit linear ephemeris is $\rm {BJD_{TDB}} = 2456193.29157 \pm 0.00021 + E\left(2.9895936 \pm 0.0000048\right)$. With a magnitude of $V\sim10$, a planetary radius of $1.699_{-0.045}^{+0.046} \ {\rm R_J}$, and a mass of $0.902_{-0.059}^{+0.060} \ {\rm M_J}$, it is the brightest host among the population of inflated Hot Jupiters ($R_P > 1.5R_J$), making it a valuable discovery for probing the nature of inflated planets. In addition, its existence within a hierarchical triple and its proximity to Earth ($210$ pc) provides a unique opportunity for dynamical studies with continued monitoring with high resolution imaging and precision radial velocities. In particular, the motion of the binary stars around each other and of both stars around the primary star relative to the measured epoch in this work should be detectable when it rises in October 2015.

preprint2015arXiv

KELT-7b: A hot Jupiter transiting a bright V=8.54 rapidly rotating F-star

We report the discovery of KELT-7b, a transiting hot Jupiter with a mass of $1.28 \pm 0.18$ MJ, radius of $1.53_{-0.047}^{+0.046}$ RJ, and an orbital period of $2.7347749 \pm 0.0000039$ days. The bright host star (HD33643; KELT-7) is an F-star with $V=8.54$, Teff $=6789_{-49}^{+50}$ K, [Fe/H] $=0.139_{-0.081}^{+0.075}$, and $\log{g}=4.149 \pm 0.019$. It has a mass of $1.535_{-0.054}^{+0.066}$ Msun, a radius of $1.732_{-0.045}^{+0.043}$ Rsun, and is the fifth most massive, fifth hottest, and the ninth brightest star known to host a transiting planet. It is also the brightest star around which KELT has discovered a transiting planet. Thus, KELT-7b is an ideal target for detailed characterization given its relatively low surface gravity, high equilibrium temperature, and bright host star. The rapid rotation of the star ($73 \pm 0.5$ km/s) results in a Rossiter-McLaughlin effect with an unusually large amplitude of several hundred m/s. We find that the orbit normal of the planet is likely to be well-aligned with the stellar spin axis, with a projected spin-orbit alignment of $λ=9.7 \pm 5.2$ degrees. This is currently the second most rapidly rotating star to have a reflex signal (and thus mass determination) due to a planetary companion measured.

preprint2015arXiv

KELT-8b: A highly inflated transiting hot Jupiter and a new technique for extracting high-precision radial velocities from noisy spectra

We announce the discovery of a highly inflated transiting hot Jupiter discovered by the KELT-North survey. A global analysis including constraints from isochrones indicates that the V = 10.8 host star (HD 343246) is a mildly evolved, G dwarf with $T_{\rm eff} = 5754_{-55}^{+54}$ K, $\log{g} = 4.078_{-0.054}^{+0.049}$, $[Fe/H] = 0.272\pm0.038$, an inferred mass $M_{*}=1.211_{-0.066}^{+0.078}$ M$_{\odot}$, and radius $R_{*}=1.67_{-0.12}^{+0.14}$ R$_{\odot}$. The planetary companion has mass $M_P = 0.867_{-0.061}^{+0.065}$ $M_{J}$, radius $R_P = 1.86_{-0.16}^{+0.18}$ $R_{J}$, surface gravity $\log{g_{P}} = 2.793_{-0.075}^{+0.072}$, and density $ρ_P = 0.167_{-0.038}^{+0.047}$ g cm$^{-3}$. The planet is on a roughly circular orbit with semimajor axis $a = 0.04571_{-0.00084}^{+0.00096}$ AU and eccentricity $e = 0.035_{-0.025}^{+0.050}$. The best-fit linear ephemeris is $T_0 = 2456883.4803 \pm 0.0007$ BJD$_{\rm TDB}$ and $P = 3.24406 \pm 0.00016$ days. This planet is one of the most inflated of all known transiting exoplanets, making it one of the few members of a class of extremely low density, highly-irradiated gas giants. The low stellar $\log{g}$ and large implied radius are supported by stellar density constraints from follow-up light curves, plus an evolutionary and space motion analysis. We also develop a new technique to extract high precision radial velocities from noisy spectra that reduces the observing time needed to confirm transiting planet candidates. This planet boasts deep transits of a bright star, a large inferred atmospheric scale height, and a high equilibrium temperature of $T_{eq}=1675^{+61}_{-55}$ K, assuming zero albedo and perfect heat redistribution, making it one of the best targets for future atmospheric characterization studies.

preprint2015arXiv

Kepler Eclipsing Binary Stars. VI. Identification of Eclipsing Binaries in the K2 Campaign 0 Data-set

The original {\it Kepler} mission observed and characterized over 2400 eclipsing binaries in addition to its prolific exoplanet detections. Despite the mechanical malfunction and subsequent non-recovery of two reaction wheels used to stabilize the instrument, the {\it Kepler} satellite continues collecting data in its repurposed {\it K2} mission surveying a series of fields along the ecliptic plane. Here we present an analysis of the first full baseline {\it K2} data release: the Campaign 0 data-set. In the 7761 light curves, we have identified a total of 207 eclipsing binaries. Of these, 97 are new discoveries that were not previously identified. Our pixel-level analysis of these objects has also resulted in identification of several false positives (observed targets contaminated by neighboring eclipsing binaries), as well as the serendipitous discovery of two short period exoplanet candidates. We provide catalog cross-matched source identifications, orbital periods, morphologies and ephemerides for these eclipsing systems. We also describe the incorporation of the K2 sample into the Kepler Eclipsing Binary Catalog\footnote{\url{keplerebs.villanova.edu/k2}}, present spectroscopic follow-up observations for a limited selection of nine systems, and discuss prospects for upcoming {\it K2} campaigns.

preprint2015arXiv

Metrics for Optimization of Large Synoptic Survey Telescope Observations of Stellar Variables and Transients

The Large Synoptic Survey Telescope (LSST) will be the largest time-domain photometric survey ever. In order to maximize the LSST science yield for a broad array of transient stellar phenomena, it is necessary to optimize the survey cadence, coverage, and depth via quantitative metrics that are specifically designed to characterize the time-domain behavior of various types of stellar transients. In this paper we present three such metrics built on the LSST Metric Analysis Framework (MAF) model (Jones et al. 2014). Two of the metrics quantify the ability of LSST to detect non-periodic and/or non-recurring transient events, and the ability of LSST to reliably measure periodic signals of various timescales. The third metric provides a way to quantify the range of stellar parameters in the stellar populations that LSST will probe. We provide example uses of these metrics and discuss some implications based on these metrics for optimization of the LSST survey for observations of stellar variables and transients.

preprint2015arXiv

Recurring Occultations of RW Aurigae by Coagulated Dust in the Tidally Disrupted Circumstellar Disk

We present photometric observations of RW Aurigae, a Classical T Tauri system, that reveal two remarkable dimming events. These events are similar to that which we observed in 2010-2011, which was the first such deep dimming observed in RW Aur in a century's worth of photometric monitoring. We suggested the 2010-2011 dimming was the result of an occultation of the star by its tidally disrupted circumstellar disk. In 2012-2013, the RW Aur system dimmed by ~0.7 mag for ~40 days and in 2014/2015 the system dimmed by ~2 mag for >250 days. The ingress/egress duration measurements of the more recent events agree well with those from the 2010-2011 event, providing strong evidence that the new dimmings are kinematically associated with the same occulting source. Therefore, we suggest that both the 2012-2013 and 2014-2015 dimming events, measured using data from the Kilodegree Extremely Little Telescope and the Kutztown University Observatory, are also occultations of RW Aur A by tidally disrupted circumstellar material. Recent hydrodynamical simulations of the eccentric fly-by of RW Aur B suggest the occulting body to be a bridge of material connecting RW Aur A and B. These simulations suggest the possibility of additional occultations, supported by the observations presented in this work. The color evolution of the dimmings suggest that the tidally stripped disk material includes dust grains ranging in size from small grains at the leading edge, typical of star forming regions, to large grains, ices or pebbles producing grey or nearly grey extinction deeper within the occulting material. It is not known whether this material represents arrested planet building prior to the tidal disruption event, or perhaps accelerated planet building as a result of the disruption event, but in any case the evidence suggests the presence of advanced planet building material in the space between RW Aur A and B.

preprint2015arXiv

SLoWPoKES-II: 100,000 Wide Binaries Identified in SDSS without Proper Motions

We present the SLoWPoKES-II catalog of low-mass visual binaries identified from the Sloan Digital Sky Survey by matching photometric distances. The candidate pairs are vetted by comparing the stellar density at their respective Galactic positions to Monte Carlo realizations of a simulated Milky Way. In this way, we are able to identify large numbers of bona fide wide binaries without the need of proper motions. 105,537 visual binaries with angular separations of $\sim$1-20", are identified, each with a probability of chance alignment of $\lesssim$5%. This is the largest catalog of bona fide wide binaries to date, and it contains a diversity of systems---in mass, mass ratios, binary separations, metallicity, and evolutionary states---that should facilitate follow-up studies to characterize the properties of M dwarfs and white dwarfs. There is a subtle but definitive suggestion of multiple populations in the physical separation distribution, supporting earlier findings. We suggest that wide binaries are comprised of multiple populations, most likely representing different formation modes. There are 141 M7 or later wide binary candidates, representing a 7-fold increase in the number currently known. These binaries are too wide to have been formed via the ejection mechanism. Finally, we find that ~6% of spectroscopically confirmed M dwarfs are not included in the SDSS STAR catalog; they are misclassified as extended sources due to the presence of a nearby or partially resolved companion. The SLoWPoKES-II catalog is publicly available to the entire community on the world wide web via the Filtergraph data visualization portal.

preprint2015arXiv

Target Selection for the SDSS-III MARVELS Survey

We present the target selection process for the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the Sloan Digital Sky Survey (SDSS) III. MARVELS is a medium-resolution ($R \sim 11000$) multi-fiber spectrograph capable of obtaining radial velocities for 60 objects at a time in order to find brown dwarfs and giant planets. The survey was configured to target dwarf stars with effective temperatures approximately between $4500$ and $6250 \, \mbox{K}$. For the first 2 years MARVELS relied on low-resolution spectroscopic pre-observations to estimate the effective temperature and $\log(g)$ for candidate stars and then selected suitable dwarf stars from this pool. Ultimately, the pre-observation spectra proved ineffective at filtering out giant stars; many giants were incorrectly classified as dwarfs, resulting in a giant contamination rate of $\sim$30\% for the first phase of the MARVELS survey. Thereafter, the survey instead applied a reduced proper motion cut to eliminate giants and used the Infrared Flux Method to estimate effective temperatures, using only extant photometric and proper-motion catalog information. The target selection method introduced here may be useful for other surveys that need to rely on extant catalog data for selection of specific stellar populations.

preprint2015arXiv

The APOGEE Spectroscopic Survey of Kepler Planet Hosts: Feasibility, Efficiency, and First Results

The Kepler mission has yielded a large number of planet candidates from among the Kepler Objects of Interest (KOIs), but spectroscopic follow-up of these relatively faint stars is a serious bottleneck in confirming and characterizing these systems. We present motivation and survey design for an ongoing project with the SDSS-III multiplexed APOGEE near-infrared spectrograph to monitor hundreds of KOI host stars. We report some of our first results using representative targets from our sample, which include current planet candidates that we find to be false positives, as well as candidates listed as false positives that we do not find to be spectroscopic binaries. With this survey, KOI hosts are observed over ~20 epochs at a radial velocity precision of 100-200 m/s. These observations can easily identify a majority of false positives caused by physically-associated stellar or substellar binaries, and in many cases, fully characterize their orbits. We demonstrate that APOGEE is capable of achieving RV precision at the 100-200 m/s level over long time baselines, and that APOGEE's multiplexing capability makes it substantially more efficient at identifying false positives due to binaries than other single-object spectrographs working to confirm KOIs as planets. These APOGEE RVs enable ancillary science projects, such as studies of fundamental stellar astrophysics or intrinsically rare substellar companions. The coadded APOGEE spectra can be used to derive stellar properties (T_eff, log(g)) and chemical abundances of over a dozen elements to probe correlations of planet properties with individual elemental abundances.

preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

preprint2015arXiv

The Puzzling Li-rich Red Giant Associated with NGC 6819

A Li-rich red giant star (2M19411367+4003382) recently discovered in the direction of NGC 6819 belongs to the rare subset of Li-rich stars that have not yet evolved to the luminosity bump, an evolutionary stage where models predict Li can be replenished. The currently favored model to explain Li enhancement in first-ascent red giants like 2M19411367+4003382 requires deep mixing into the stellar interior. Testing this model requires a measurement of 12C/13C, which is possible to obtain from APOGEE spectra. However, the Li-rich star also has abnormal asteroseismic properties that call into question its membership in the cluster, even though its radial velocity and location on color-magnitude diagrams are consistent with membership. To address these puzzles, we have measured a wide array of abundances in the Li-rich star and three comparison stars using spectra taken as part of the APOGEE survey to determine the degree of stellar mixing, address the question of membership, and measure the surface gravity. We confirm that the Li-rich star is a red giant with the same overall chemistry as the other cluster giants. However, its log g is significantly lower, consistent with the asteroseismology results and suggestive of a very low mass if the star is indeed a cluster member. Regardless of the cluster membership, the 12C/13C and C/N ratios of the Li-rich star are consistent with standard first dredge-up, indicating that Li dilution has already occurred, and inconsistent with internal Li enrichment scenarios that require deep mixing.

preprint2015arXiv

Transiting Planets with LSST II. Period Detection of Planets Orbiting 1 Solar Mass Hosts

The Large Synoptic Survey Telescope (LSST) will photometrically monitor approximately 1 billion stars for ten years. The resulting light curves can be used to detect transiting exoplanets. In particular, as demonstrated by Lund et al. (2015), LSST will probe stellar populations currently undersampled in most exoplanet transit surveys, including out to extragalactic distances. In this paper we test the efficiency of the box-fitting least-squares (BLS) algorithm for accurately recovering the periods of transiting exoplanets using simulated LSST data. We model planets with a range of radii orbiting a solar-mass star at a distance of 7 kpc, with orbital periods ranging from 0.5 to 20 days. We find that standard-cadence LSST observations will be able to reliably recover the periods of Hot Jupiters with periods shorter than approximately 3 days, however it will remain a challenge to confidently distinguish these transiting planets from false positives. At the same time, we find that the LSST deep drilling cadence is extremely powerful: the BLS algorithm successfully recovers at least 30% of sub-Saturn-size exoplanets with orbital periods as long as 20 days, and a simple BLS power criterion robustly distinguishes approximately 98% of these from photometric (i.e. statistical) false positives.

preprint2015arXiv

V409 Tau As Another AA Tau: Photometric Observations of Stellar Occultations by the Circumstellar Disk

AA Tau is a well studied young stellar object that presents many of the photometric characteristics of a Classical T Tauri star (CTTS), including short-timescale stochastic variability attributed to spots and/or accretion as well as long duration dimming events attributed to occultations by vertical features (e.g., warps) in its circumstellar disk. We present new photometric observations of AA Tau from the Kilodegree Extremely Little Telescope North (KELT-North) which reveal a deep, extended dimming event in 2011, which we show supports the interpretation by Bouvier et al. (2013) of an occultation by a high-density feature in the circumstellar disk located >8 AU from the star. We also present KELT-North observations of V409 Tau, a relatively unstudied young stellar object also in Taurus-Auriga, showing short timescale erratic variability, along with two separate long and deep dimming events, one from January 2009 through late October 2010, and the other from March 2012 until at least September 2013. We interpret both dimming events to have lasted more than 600 days, each with a depth of ~1.4 mag. From a spectral energy distribution analysis, we propose that V409 Tau is most likely surrounded by a circumstellar disk viewed nearly edge-on, and using Keplerian timescale arguments we interpret the deep dimmings of V409 Tau as occultations from one or more features within this disk >10 AU from the star. In both AA Tau and V409 Tau, the usual CTTS short-timescale variations associated with accretion processes close to the stars continue during the occultations, further supporting the distant occulting material interpretation. Like AA Tau, V409 Tau serves as a laboratory for studying the detailed structure of the protoplanetary environments of T Tauri disks, specifically disk structures that may be signposts of planet formation at many AU out in the disk.

preprint2014arXiv

A Strict Test of Stellar Evolution Models: The Absolute Dimensions of Massive Benchmark Eclipsing Binary V578 Mon

We determine the absolute dimensions of the eclipsing binary V578 Mon, a detached system of two early B-type stars (B0V + B1V, P$=$2.40848 d) in the star-forming region NGC 2244 of the Rosette Nebula. From the light curve analysis of 40 yr of photometry and the analysis of HERMES spectra, we find radii of $5.41\pm0.04$ Rsun and $4.29\pm 0.05$ Rsun, and temperatures of $ 30000\pm 500$~K and $ 25750\pm 435$ K respectively. We find that our disentangled component spectra for V578 Mon agree well previous spectral disentangling from the literature. We also reconfirm the previous spectroscopic orbit of V578 Mon finding that masses of $ 14.54\pm 0.08$ Msun and $ 10.29\pm 0.06$ Msun are fully compatible with the new analysis. We compare the absolute dimensions to the rotating models of the Geneva and Utrecht groups and the models of Granada group. We find all three sets of models marginally reproduce the absolute dimensions of both stars with a common age within uncertainty for gravity-effective temperature isochrones. However - there are some apparent age discrepancies for the corresponding mass-radius isochrones. Models with larger convective overshoot $>0.35$ worked best. Combined with our previously determined apsidal motion of $0.07089^{+0.00021}_{-0.00013}$ deg cycle$^{-1}$, we compute the internal structure constants (tidal Love number) for the newtonian and general relativistic contribution to the apsidal motion, $\log{k_2}=-1.975\pm0.017$ and $\log{k_2}=-3.412\pm0.018$ respectively. We find the relativistic contribution to the apsidal motion of be small $<4\%$. We find that the prediction of $\log{k_{\rm 2,theo}}=-2.005\pm0.025$ of the Granada models fully agrees with our observed $\log{k_2}$.

preprint2014arXiv

Accurate Atmospheric Parameters at Moderate Resolution Using Spectral Indices: Preliminary Application to the MARVELS Survey

Studies of Galactic chemical and dynamical evolution in the solar neighborhood depend on the availability of precise atmospheric parameters (Teff, [Fe/H] and log g) for solar-type stars. Many large-scale spectroscopic surveys operate at low to moderate spectral resolution for efficiency in observing large samples, which makes the stellar characterization difficult due to the high degree of blending of spectral features. While most surveys use spectral synthesis, in this work we employ an alternative method based on spectral indices to determine the atmospheric parameters of a sample of nearby FGK dwarfs and subgiants observed by the MARVELS survey at moderate resolving power (R~12,000). We have developed three codes to automatically normalize the observed spectra, measure the equivalent widths of the indices and, through the comparison of those with values calculated with pre-determined calibrations, derive the atmospheric parameters of the stars. The calibrations were built using a sample of 309 stars with precise stellar parameters obtained from the analysis of high-resolution FEROS spectra. A validation test of the method was conducted with a sample of 30 MARVELS targets that also have reliable atmospheric parameters from high-resolution spectroscopic analysis. Our approach was able to recover the parameters within 80 K for Teff, 0.05 dex for [Fe/H] and 0.15 dex for log g, values that are lower or equal to the typical external uncertainties found between different high-resolution analyzes. An additional test was performed with a subsample of 138 stars from the ELODIE stellar library and the literature atmospheric parameters were recovered within 125 K for Teff, 0.10 dex for [Fe/H] and 0.29 dex for log g. These results show that the spectral indices are a competitive tool to characterize stars with the intermediate resolution spectra.

preprint2014arXiv

An Improved Determination of the Lithium Depletion Boundary Age of Blanco 1 and a First Look on the Effects of Magnetic Activity

The Lithium Depletion Boundary (LDB) is a robust method for accurately determining the ages of young clusters, but most pre-main-sequence models used to derive LDB ages do not include the effects of magnetic activity on stellar properties. In light of this, we present results from our spectroscopic study of the very low-mass members of the southern open cluster Blanco 1 using the Gemini-North telescope, program IDs: GN-2009B-Q-53 and GN-2010B-Q-96. We obtained GMOS spectra at intermediate resolution for cluster candidate members with $I$$\approx$13--20 mag. From our sample of 43 spectra, we find 14 probable cluster members by considering proximity to the cluster sequence in an $I/I-K\rm_s$ color-magnitude diagram, agreement with the cluster's systemic radial velocity, and magnetic activity as a youth indicator. We systematically analyze the H$α$ and Li features and update the LDB age of Blanco 1 to be $126^{+13}_{-14}$ Myr. Our new LDB age for Blanco 1 shows remarkable coevality with the benchmark Pleiades open cluster. Using available empirical activity corrections, we investigate the effects of magnetic activity on the LDB age of Blanco 1. Accounting for activity, we infer a corrected LDB age of $114^{+9}_{-10}$ Myr. This work demonstrates the importance of accounting for magnetic activity on LDB inferred stellar ages, suggesting the need to re-investigate previous LDB age determinations.

preprint2014arXiv

Detailed Abundances of Planet-Hosting Wide Binaries. I. Did Planet Formation Imprint Chemical Signatures in the Atmospheres of HD 20782/81?

Using high-resolution echelle spectra obtained with Magellan/MIKE, we present a chemical abundance analysis of both stars in the planet-hosting wide binary system HD20782 + HD20781. Both stars are G dwarfs, and presumably coeval, forming in the same molecular cloud. Therefore we expect that they should possess the same bulk metallicities. Furthermore, both stars also host giant planets on eccentric orbits with pericenters $\lesssim 0.2\,$ AU. We investigate if planets with such orbits could lead to the host stars ingesting material, which in turn may leave similar chemical imprints in their atmospheric abundances. We derived abundances of 15 elements spanning a range of condensation temperatures ($T_{C}\approx 40-1660\,$ K). The two stars are found to have a mean element-to-element abundance difference of $0.04\pm0.07\,$ dex, which is consistent with both stars having identical bulk metallicities. In addition, for both stars, the refractory elements ($T_{C} > 900\,$ K) exhibit a positive correlation between abundance (relative to solar) and $T_{C}$, with similar slopes of $\approx$ $1\times10^{-4}\,$ dex K$^{-1}$. The measured positive correlations are not perfect; both stars exhibit a scatter of $\approx$ $5\times10^{-5}\,$ dex K$^{-1}$ about the mean trend, and certain elements (Na, Al, Sc) are similarly deviant in both stars. These findings are discussed in the context of models for giant planet migration that predict the accretion of H-depleted rocky material by the host star. We show that a simple simulation of a solar-type star accreting material with Earth-like composition predicts a positive---but imperfect---correlation between refractory elemental abundances and $T_{C}$. Our measured slopes for HD 20782/81 are consistent with what is predicted for the ingestion of 10--20 Earths by both stars.

preprint2014arXiv

Discovery of a Transiting Planet Near the Snow-Line

In most theories of planet formation, the snow-line represents a boundary between the emergence of the interior rocky planets and the exterior ice giants. The wide separation of the snow-line makes the discovery of transiting worlds challenging, yet transits would allow for detailed subsequent characterization. We present the discovery of Kepler-421b, a Uranus-sized exoplanet transiting a G9/K0 dwarf once every 704.2 days in a near-circular orbit. Using public Kepler photometry, we demonstrate that the two observed transits can be uniquely attributed to the 704.2 day period. Detailed light curve analysis with BLENDER validates the planetary nature of Kepler-421b to >4 sigmas confidence. Kepler-421b receives the same insolation as a body at ~2AU in the Solar System and for a Uranian albedo would have an effective temperature of ~180K. Using a time-dependent model for the protoplanetary disk, we estimate that Kepler-421b's present semi-major axis was beyond the snow-line after ~3Myr, indicating that Kepler-421b may have formed at its observed location.

preprint2014arXiv

Empirical Tests of Pre-Main-Sequence Stellar Evolution Models with Eclipsing Binaries

We examine the performance of standard PMS stellar evolution models against the accurately measured properties of a benchmark sample of 26 PMS stars in 13 EB systems. We provide a definitive compilation of all fundamental properties for the EBs. We also provide a definitive compilation of the various PMS model sets. In the H-R diagram, the masses inferred for the individual stars by the models are accurate to better than 10% above 1 Msun, but below 1 Msun they are discrepant by 50-100%. We find evidence that the failure of the models to match the data is linked to the triples in the EB sample; at least half of the EBs possess tertiary companions. Excluding the triples, the models reproduce the stellar masses to better than ~10% in the H-R diagram, down to 0.5 Msun, below which the current sample is fully contaminated by tertiaries. We consider several mechanisms by which a tertiary might cause changes in the EB properties and thus corrupt the agreement with stellar model predictions. We show that the energies of the tertiary orbits are comparable to that needed to potentially explain the scatter in the EB properties through injection of heat, perhaps involving tidal interaction. It seems from the evidence at hand that this mechanism, however it operates in detail, has more influence on the surface properties of the stars than on their internal structure, as the lithium abundances are broadly in good agreement with model predictions. The EBs that are members of young clusters appear individually coeval to within 20%, but collectively show an apparent age spread of ~50%, suggesting true age spreads in young clusters. However, this apparent spread in the EB ages may also be the result of scatter in the EB properties induced by tertiaries. [Abridged]

preprint2014arXiv

Flicker as a tool for characterizing planets through Asterodensity Profiling

Variability in the time series brightness of a star on a timescale of 8 hours, known as 'flicker', has been previously demonstrated to serve as a proxy for the surface gravity of a star by Bastien et al. (2013). Although surface gravity is crucial for stellar classification, it is the mean stellar density which is most useful when studying transiting exoplanets, due to its direct impact on the transit light curve shape. Indeed, an accurate and independent measure of the stellar density can be leveraged to infer subtle properties of a transiting system, such as the companion's orbital eccentricity via asterodensity profiling. We here calibrate flicker to the mean stellar density of 439 Kepler targets with asteroseismology, allowing us to derive a new empirical relation given by $\log_{10}(ρ_{\star}\,[\mathrm{kg}\,\mathrm{m}^{-3}]) = 5.413 - 1.850 \log_{10}(F_8\,[\mathrm{ppm}])$. The calibration is valid for stars with $4500$K$<T_{\mathrm{eff}}<6500$K, $K_P<14$ and flicker estimates corresponding to stars with $3.25<\log g_{\star}<4.43$. Our relation has a model error in the stellar density of 31.7% and so has $\sim8$ times lower precision than that from asteroseismology but is applicable to a sample $\sim40$ times greater. Flicker therefore provides an empirical method to enable asterodensity profiling on hundreds of planetary candidates from present and future missions.

preprint2014arXiv

IN-SYNC I: Homogeneous Stellar Parameters from High Resolution APOGEE Spectra for Thousands of Pre-main Sequence Star

Over two years 8,859 high-resolution H-band spectra of 3,493 young (1 - 10 Myr) stars were gathered by the multi-object spectrograph of the APOGEE project as part of the IN-SYNC ancillary program of that SDSS-III survey. Here we present the forward modeling approach used to derive effective temperatures, surface gravities, radial velocities, rotational velocities, and H-band veiling from these near-infrared spectra. We discuss in detail the statistical and systematic uncertainties in these stellar parameters. In addition we present accurate extinctions by measuring the E(J-H) of these young stars with respect to the single-star photometric locus in the Pleiades. Finally we identify an intrinsic stellar radius spread of about 25% for late-type stars in IC 348 using three (nearly) independent measures of stellar radius, namely the extinction-corrected J-band magnitude, the surface gravity and the $R \sin i$ from the rotational velocities and literature rotation periods. We exclude that this spread is caused by uncertainties in the stellar parameters by showing that the three estimators of stellar radius are correlated, so that brighter stars tend to have lower surface gravities and larger $R \sin i$ than fainter stars at the same effective temperature. Tables providing the spectral and photometric parameters for the Pleiades and IC 348 have been provided online.

preprint2014arXiv

IN-SYNC II: Virial Stars from Sub-Virial Cores -- The Velocity Dispersion of Embedded Pre-Main-Sequence Stars in NGC 1333

The initial velocity dispersion of newborn stars is a major unconstrained aspect of star formation theory. Using near-infrared spectra obtained with the APOGEE spectrograph, we show that the velocity dispersion of young (1-2 Myr) stars in NGC 1333 is 0.92+/-0.12 km/s after correcting for measurement uncertainties and the effect of binaries. This velocity dispersion is consistent with the virial velocity of the region and the diffuse gas velocity dispersion, but significantly larger than the velocity dispersion of the dense, star-forming cores, which have a sub-virial velocity dispersion of 0.5 km/s. Since the NGC 1333 cluster is dynamically young and deeply embedded, this measurement provides a strong constraint on the initial velocity dispersion of newly-formed stars. We propose that the difference in velocity dispersion between stars and dense cores may be due to the influence of a 70 micro-Gauss magnetic field acting on the dense cores, or be the signature of a cluster with initial sub-structure undergoing global collapse.

preprint2014arXiv

KELT-6b: A P~7.9 d Hot Saturn Transiting a Metal-Poor Star with a Long-Period Companion

We report the discovery of KELT-6b, a mildly-inflated Saturn-mass planet transiting a metal-poor host. The initial transit signal was identified in KELT-North survey data, and the planetary nature of the occulter was established using a combination of follow-up photometry, high-resolution imaging, high-resolution spectroscopy, and precise radial velocity measurements. The fiducial model from a global analysis including constraints from isochrones indicates that the V=10.38 host star (BD+31 2447) is a mildly evolved, late-F star with T_eff=6102 \pm 43 K, log(g_*)=4.07_{-0.07}^{+0.04} and [Fe/H]=-0.28 \pm 0.04, with an inferred mass M_*=1.09 \pm 0.04 M_sun and radius R_star=1.58_{-0.09}^{+0.16} R_sun. The planetary companion has mass M_P=0.43 \pm 0.05 M_J, radius R_P=1.19_{-0.08}^{+0.13} R_J, surface gravity log(g_P)=2.86_{-0.08}^{+0.06}, and density rho_P=0.31_{-0.08}^{+0.07} g~cm^{-3}. The planet is on an orbit with semimajor axis a=0.079 \pm 0.001 AU and eccentricity e=0.22_{-0.10}^{+0.12}, which is roughly consistent with circular, and has ephemeris of T_c(BJD_TDB)=2456347.79679 \pm 0.00036 and P=7.845631 \pm 0.000046 d. Equally plausible fits that employ empirical constraints on the host star parameters rather than isochrones yield a larger planet mass and radius by ~4-7%. KELT-6b has surface gravity and incident flux similar to HD209458b, but orbits a host that is more metal poor than HD209458 by ~0.3 dex. Thus, the KELT-6 system offers an opportunity to perform a comparative measurement of two similar planets in similar environments around stars of very different metallicities. The precise radial velocity data also reveal an acceleration indicative of a longer-period third body in the system, although the companion is not detected in Keck adaptive optics images.

preprint2014arXiv

Kepler Eclipsing Binary Stars. V. Identification of 31 Eclipsing Binaries in the K2 Engineering Data-set

Over 2500 eclipsing binaries were identified and characterized from the ultra-precise photometric data provided by the Kepler space telescope. Kepler is now beginning its second mission, K2, which is proving to again provide ultra-precise photometry for a large sample of eclipsing binary stars. In the 1951 light curves covering 12 days in the K2 engineering data-set, we have identified and determined the ephemerides for 31 eclipsing binaries that demonstrate the capabilities for eclipsing binary science in the upcoming campaigns in K2. Of those, 20 are new discoveries. We describe both manual and automated approaches to harvesting the complete set of eclipsing binaries in the K2 data, provide identifications and details for the full set of eclipsing binaries present in the engineering data-set, and discuss the prospects for application of eclipsing binary searches in the K2 mission.

preprint2014arXiv

Larger Planet Radii Inferred from Stellar "Flicker" Brightness Variations of Bright Planet Host Stars

Most extrasolar planets have been detected by their influence on their parent star, typically either gravitationally (the Doppler method) or by the small dip in brightness as the planet blocks a portion of the star (the transit method). Therefore, the accuracy with which we know the masses and radii of extrasolar planets depends directly on how well we know those of the stars, the latter usually determined from the measured stellar surface gravity, logg. Recent work has demonstrated that the short-timescale brightness variations ("flicker") of stars can be used to measure logg to a high accuracy of ~0.1-0.2 dex (Bastien et al. 2013). Here, we use flicker measurements of 289 bright (Kepmag<13) candidate planet-hosting stars with Teff=4500-6650 K to re-assess the stellar parameters and determine the resulting impact on derived planet properties. This re-assessment reveals that for the brightest planet-host stars, an astrophysical bias exists that contaminates the stellar sample with evolved stars: nearly 50% of the bright planet-host stars are subgiants. As a result, the stellar radii, and hence the radii of the planets orbiting these stars, are on average 20-30% larger than previous measurements had suggested.

preprint2014arXiv

Multiplicity in Early Stellar Evolution

Observations from optical to centimeter wavelengths have demonstrated that multiple systems of two or more bodies is the norm at all stellar evolutionary stages. Multiple systems are widely agreed to result from the collapse and fragmentation of cloud cores, despite the inhibiting influence of magnetic fields. Surveys of Class 0 protostars with mm interferometers have revealed a very high multiplicity frequency of about 2/3, even though there are observational difficulties in resolving close protobinaries, thus supporting the possibility that all stars could be born in multiple systems. Near-infrared adaptive optics observations of Class I protostars show a lower binary frequency relative to the Class 0 phase, a declining trend that continues through the Class II/III stages to the field population. This loss of companions is a natural consequence of dynamical interplay in small multiple systems, leading to ejection of members. We discuss observational consequences of this dynamical evolution, and its influence on circumstellar disks, and we review the evolution of circumbinary disks and their role in defining binary mass ratios. Special attention is paid to eclipsing PMS binaries, which allow for observational tests of evolutionary models of early stellar evolution. Many stars are born in clusters and small groups, and we discuss how interactions in dense stellar environments can significantly alter the distribution of binary separations through dissolution of wider binaries. The binaries and multiples we find in the field are the survivors of these internal and external destructive processes, and we provide a detailed overview of the multiplicity statistics of the field, which form a boundary condition for all models of binary evolution. Finally we discuss various formation mechanisms for massive binaries, and the properties of massive trapezia.

preprint2014arXiv

New BVIc Photometry of Low-mass Pleiades Stars: Exploring the Effects of Rotation on Broadband Colors

We present new BVIc photometry for 350 Pleiades proper-motion members with 9 < V < 17. Importantly, our new catalog includes a large number of K and early-M type stars, roughly doubling the number of low-mass stars with well calibrated Johnson/Cousins photometry in this benchmark cluster. We combine our new photometry with existing photometry from the literature to define a purely empirical isochrone at Pleiades age (~100 Myr) extending from V=9 to 17. We use the empirical isochrone to identify 48 new probable binaries and 14 likely non-members.The photometrically identified single stars are compared against their expected positions in the color-magnitude diagram (CMD). At 100 Myr, the mid K and early M stars are predicted to lie above the zero-age main sequence (ZAMS) having not yet reached the ZAMS. We find in the B-V vs. V CMD that mid K and early M dwarfs are instead displaced below (or blueward of) the ZAMS. Using the stars' previously reported rotation periods, we find a highly statistically significant correlation between rotation period and CMD displacement, in the sense that the more rapidly rotating stars have the largest displacements in the B- V CMD.

preprint2014arXiv

New Red Jewels in Coma Berenices

We have used Sloan Digital Sky Survey-III (SDSS-III) Apache Point Observatory Galactic Evolution Experiment (APOGEE) radial velocity observations in the near-infrared $H$-band to explore the membership of the nearby ($86.7 \pm 0.9$ pc) open cluster Coma Berenices (Melotte 111), concentrating on the poorly-populated low-mass end of the main sequence. Using SDSS-III APOGEE radial velocity measurements, we confirm the membership of eight K/M dwarf members, providing the first confirmed low-mass members of the Coma Berenices cluster. Using $R\sim2000$ spectra from IRTF-SpeX, we confirm the independently luminosity classes of these targets, and find their metallicities to be consistent with the known solar mean metallicity of Coma Berenices and of M dwarfs in the Solar neighborhood. In addition, the APOGEE spectra have enabled measurement of $v\sin i$ for each target and detection for the first time of the low-mass secondary components of the known binary systems Melotte 111 102 and Melotte 111 120, as well as identification of the previously unknown binary system 2MASS J12214070+2707510. Finally, we use Kilodegree Extremely Little Telescope (KELT) photometry to measure photometric variability and rotation periods for a subset of the Coma Berenices members.

preprint2014arXiv

Observations of the M82 SN 2014J with the Kilodegree Extremely Little Telescope

We report observations of the bright M82 supernova 2014J serendipitously obtained with the Kilodegree Extremely Little Telescope (KELT). The SN was observed at high cadence for over 100 days, from pre-explosion, to early rise and peak times, through the secondary bump. The high cadence KELT data with high S/N is completely unique for SN 2014J and for any other SNIa, with the exception of the (yet) unpublished Kepler data. Here, we report determinations of the SN explosion time and peak time. We also report measures of the "smoothness" of the light curve on timescales of minutes/hours never before probed, and we use this to place limits on energy produced from short-lived isotopes or inhomogeneities in the explosion or the circumstellar medium. From the non-observation of significant perturbations of the light curves, we derive a 3sigma upper-limit corresponding to 8.7 x 10^36 erg/s for any such extra sources of luminosity at optical wavelengths.

preprint2014arXiv

Sodium and Oxygen Abundances in the Open Cluster NGC 6791 from APOGEE H-Band Spectroscopy

The open cluster NGC 6791 is among the oldest, most massive and metal-rich open clusters in the Galaxy. High-resolution $H$-band spectra from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) of 11 red giants in NGC 6791 are analyzed for their chemical abundances of iron, oxygen, and sodium. The abundances of these three elements are found to be homogeneous (with abundance dispersions at the level of $\sim$ 0.05 - 0.07 dex) in these cluster red giants, which span much of the red-giant branch (T$_{\rm eff}$ $\sim$ 3500K - 4600K), and include two red-clump giants. From the infrared spectra, this cluster is confirmed to be among the most metal-rich clusters in the Galaxy ($<$[Fe/H]$>$ = 0.34 $\pm$ 0.06), and is found to have a roughly solar value of [O/Fe] and slightly enhanced [Na/Fe]. Non-LTE calculations for the studied Na I lines in the APOGEE spectral region ($λ$16373.86Å and $λ$16388.85Å) indicate only small departures from LTE ($\leq$ 0.04 dex) for the parameter range and metallicity of the studied stars. The previously reported double population of cluster members with different Na abundances is not found among the studied sample.

preprint2014arXiv

Spitzer and z' Secondary Eclipse Observations of the Highly Irradiated Transiting Brown Dwarf KELT-1b

We present secondary eclipse observations of the highly irradiated transiting brown dwarf KELT-1b. These observations represent the first constraints on the atmospheric dynamics of a highly irradiated brown dwarf, and the atmospheres of irradiated giant planets at high surface gravity. Using the Spitzer Space Telescope, we measure secondary eclipse depths of 0.195+/-0.010% at 3.6um and 0.200+/-0.012% at 4.5um. We also find tentative evidence for the secondary eclipse in the z' band with a depth of 0.049+/-0.023%. These measured eclipse depths are most consistent with an atmosphere model in which there is a strong substellar hotspot, implying that heat redistribution in the atmosphere of KELT-1b is low. While models with a more mild hotspot or even with dayside heat redistribution are only marginally disfavored, models with complete heat redistribution are strongly ruled out. The eclipse depths also prefer an atmosphere with no TiO inversion layer, although a model with TiO inversion is permitted in the dayside heat redistribution case, and we consider the possibility of a day-night TiO cold trap in this object. For the first time, we compare the IRAC colors of brown dwarfs and hot Jupiters as a function of effective temperature. Importantly, our measurements reveal that KELT-1b has a [3.6]-[4.5] color of 0.07+/-0.11, identical to that of isolated brown dwarfs of similarly high temperature. In contrast, hot Jupiters generally show redder [3.6]-[4.5] colors of ~0.4, with a very large range from ~0 to ~1. Evidently, despite being more similar to hot Jupiters than to isolated brown dwarfs in terms of external forcing of the atmosphere by stellar insolation, KELT-1b has an atmosphere most like that of other brown dwarfs. This suggests that surface gravity is very important in controlling the atmospheric systems of substellar mass bodies.

preprint2014arXiv

Testing the Asteroseismic Mass Scale Using Metal-Poor Stars Characterized with APOGEE and Kepler

Fundamental stellar properties, such as mass, radius, and age, can be inferred using asteroseismology. Cool stars with convective envelopes have turbulent motions that can stochastically drive and damp pulsations. The properties of the oscillation frequency power spectrum can be tied to mass and radius through solar-scaled asteroseismic relations. Stellar properties derived using these scaling relations need verification over a range of metallicities. Because the age and mass of halo stars are well-constrained by astrophysical priors, they provide an independent, empirical check on asteroseismic mass estimates in the low-metallicity regime. We identify nine metal-poor red giants (including six stars that are kinematically associated with the halo) from a sample observed by both the Kepler space telescope and the Sloan Digital Sky Survey-III APOGEE spectroscopic survey. We compare masses inferred using asteroseismology to those expected for halo and thick-disk stars. Although our sample is small, standard scaling relations, combined with asteroseismic parameters from the APOKASC Catalog, produce masses that are systematically higher (<ΔM>=0.17+/-0.05 Msun) than astrophysical expectations. The magnitude of the mass discrepancy is reduced by known theoretical corrections to the measured large frequency separation scaling relationship. Using alternative methods for measuring asteroseismic parameters induces systematic shifts at the 0.04 Msun level. We also compare published asteroseismic analyses with scaling relationship masses to examine the impact of using the frequency of maximum power as a constraint. Upcoming APOKASC observations will provide a larger sample of ~100 metal-poor stars, important for detailed asteroseismic characterization of Galactic stellar populations.

preprint2014arXiv

The Impact of Chromospheric Activity on Observed Initial Mass Functions

Using recently established empirical calibrations for the impact of chromospheric activity on the radii, effective temperatures, and estimated masses of active low-mass stars and brown dwarfs, we reassess the shape of the initial mass function (IMF) across the stellar/substellar boundary in the Upper Sco star-forming region (age 5-10 Myr). We adjust the observed effective temperatures to warmer values using the observed strength of the chromospheric H$α$ emission, and redetermine the estimated masses of objects using pre--main-sequence evolutionary tracks in the H-R diagram. The effect of the activity-adjusted temperatures is to shift the objects to higher masses by 3-100%. While the slope of the resulting IMF at substellar masses is not strongly changed, the peak of the IMF does shift from ~0.06 to ~0.11 Msun. Moreover, for objects with masses <~0.2 Msun, the ratio of brown dwarfs to stars changes from ~80% to ~33%. These results suggest that activity corrections are essential for studies of the substellar mass function, if the masses are estimated from spectral types or from effective temperatures.

preprint2014arXiv

The Tenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Apache Point Observatory Galactic Evolution Experiment

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

preprint2014arXiv

The Transiting Exoplanet Survey Satellite

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

preprint2014arXiv

Transiting Planets with LSST I: Potential for LSST Exoplanet Detection

The Large Synoptic Survey Telescope (LSST) has been designed in order to satisfy several different scientific objectives that can be addressed by a ten-year synoptic sky survey. However, LSST will also provide a large amount of data that can then be exploited for additional science beyond its primary goals. We demonstrate the potential of using LSST data to search for transiting exoplanets, and in particular to find planets orbiting host stars that are members of stellar populations that have been less thoroughly probed by current exoplanet surveys. We find that existing algorithms can detect in simulated LSST light curves the transits of Hot Jupiters around solar-type stars, Hot Neptunes around K dwarfs, and planets orbiting stars in the Large Magellanic Cloud. We also show that LSST would have the sensitivity to potentially detect Super-Earths orbiting red dwarfs, including those in habitable zone orbits, if they are present in some fields that LSST will observe. From these results, we make the case that LSST has the ability to provide a valuable contribution to exoplanet science.

preprint2013arXiv

A Cautionary Tale: MARVELS Brown Dwarf Candidate Reveals Itself To Be A Very Long Period, Highly Eccentric Spectroscopic Stellar Binary

We report the discovery of a highly eccentric, double-lined spectroscopic binary star system (TYC 3010-1494-1), comprising two solar-type stars that we had initially identified as a single star with a brown dwarf companion. At the moderate resolving power of the MARVELS spectrograph and the spectrographs used for subsequent radial-velocity (RV) measurements (R ~ <30,000), this particular stellar binary mimics a single-lined binary with an RV signal that would be induced by a brown dwarf companion (Msin(i)~50 M_Jup) to a solar-type primary. At least three properties of this system allow it to masquerade as a single star with a very low-mass companion: its large eccentricity (e~0.8), its relatively long period (P~238 days), and the approximately perpendicular orientation of the semi-major axis with respect to the line of sight (omega~189 degrees). As a result of these properties, for ~95% of the orbit the two sets of stellar spectral lines are completely blended, and the RV measurements based on centroiding on the apparently single-lined spectrum is very well fit by an orbit solution indicative of a brown dwarf companion on a more circular orbit (e~0.3). Only during the ~5% of the orbit near periastron passage does the true, double-lined nature and large RV amplitude of ~15 km/s reveal itself. The discovery of this binary system is an important lesson for RV surveys searching for substellar companions; at a given resolution and observing cadence, a survey will be susceptible to these kinds of astrophysical false positives for a range of orbital parameters. Finally, for surveys like MARVELS that lack the resolution for a useful line bisector analysis, it is imperative to monitor the peak of the cross-correlation function for suspicious changes in width or shape, so that such false positives can be flagged during the candidate vetting process.

preprint2013arXiv

An observational correlation between stellar brightness variations and surface gravity

Surface gravity is one of a star's basic properties, but it is difficult to measure accurately, with typical uncertainties of 25-50 per cent if measured spectroscopically and 90-150 per cent photometrically. Asteroseismology measures gravity with an uncertainty of about two per cent but is restricted to relatively small samples of bright stars, most of which are giants. The availability of high-precision measurements of brightness variations for >150,000 stars provides an opportunity to investigate whether the variations can be used to determine surface gravities. The Fourier power of granulation on a star's surface correlates physically with surface gravity; if brightness variations on timescales of hours arise from granulation, then such variations should correlate with surface gravity. Here we report an analysis of archival data that reveals an observational correlation between surface gravity and the root-mean-square brightness variations on timescales of less than eight hours for stars with temperatures of 4500-6750K, log of surface gravities of 2.5-4.5 (cgs units), and having overall brightness variations <3 parts per thousand. A straightforward observation of optical brightness variations therefore allows a determination of the surface gravity with a precision of <25 percent for inactive Sun-like stars at main-sequence to giant stages of evolution.

preprint2013arXiv

Constraining Explosion Type of Young Supernova Remnants Using 24 Micron Emission Morphology

Determination of the explosion type of supernova remnants (SNRs) can be challenging, as SNRs are hundreds to thousands of years old and supernovae (SNe) are classified based on spectral properties days after explosion. Previous studies of thermal X-ray emission from Milky Way and Large Magellanic Cloud (LMC) SNRs have shown that Type Ia and core-collapse (CC) SNRs have statistically different symmetries, and thus these sources can be typed based on their X-ray morphologies. In this paper, we extend the same technique, a multipole expansion technique using power ratios, to infrared (IR) images of SNRs to test whether they can be typed using the symmetry of their warm dust emission as well. We analyzed archival Spitzer Space Telescope Multiband Imaging Photometer (MIPS) 24 micron observations of the previously used X-ray sample, and we find that the two classes of SNRs separate according to their IR morphologies. The Type Ia SNRs are statistically more circular and mirror symmetric than the CC SNRs, likely due to the different circumstellar environments and explosion geometries of the progenitors. Broadly, our work indicates that the IR emission retains information of the explosive origins of the SNR and offers a new method to type SNRs based on IR morphology.

preprint2013arXiv

Filtergraph: An Interactive Web Application for Visualization of Astronomy Datasets

Filtergraph is a web application being developed and maintained by the Vanderbilt Initiative in Data-intensive Astrophysics (VIDA) to flexibly and rapidly visualize a large variety of astronomy datasets of various formats and sizes. The user loads a flat-file dataset into Filtergraph which automatically generates an interactive data portal that can be easily shared with others. From this portal, the user can immediately generate scatter plots of up to 5 dimensions as well as histograms and tables based on the dataset. Key features of the portal include intuitive controls with auto-completed variable names, the ability to filter the data in real time through user-specified criteria, the ability to select data by dragging on the screen, and the ability to perform arithmetic operations on the data in real time. To enable seamless data visualization and exploration, changes are quickly rendered on screen and visualizations can be exported as high quality graphics files. The application is optimized for speed in the context of large datasets: for instance, a plot generated from a stellar database of 3.1 million entries renders in less than 2 seconds on a standard web server platform. This web application has been created using the Web2py web framework based on the Python programming language. Filtergraph is free to use at http://filtergraph.vanderbilt.edu/.

preprint2013arXiv

KELT-3b: A Hot Jupiter Transiting a V=9.8 Late-F Star

We report the discovery of KELT-3b, a moderately inflated transiting hot Jupiter with a mass of 1.477 (-0.067, +0.066) M_J, and radius of 1.345 +/- 0.072 R_J, with an orbital period of 2.7033904 +/- 0.000010 days. The host star, KELT-3, is a V=9.8 late F star with M_* = 1.278 (-0.061, +0.063) M_sun, R_* = 1.472 (-0.067, +0.065) R_sun, T_eff = 6306 (-49, +50) K, log(g) = 4.209 (-0.031, +0.033), and [Fe/H] = 0.044 (-0.082, +0.080), and has a likely proper motion companion. KELT-3b is the third transiting exoplanet discovered by the KELT survey, and is orbiting one of the 20 brightest known transiting planet host stars, making it a promising candidate for detailed characterization studies. Although we infer that KELT-3 is significantly evolved, a preliminary analysis of the stellar and orbital evolution of the system suggests that the planet has likely always received a level of incident flux above the empirically-identified threshold for radius inflation suggested by Demory & Seager (2011).

preprint2013arXiv

Kepler Eclipsing Binary Stars. IV. Precise Eclipse Times for Close Binaries and Identification of Candidate Three-Body Systems

We present a catalog of precise eclipse times and analysis of third body signals among 1279 close binaries in the latest Kepler Eclipsing Binary Catalog. For these short period binaries, Kepler's 30 minute exposure time causes significant smearing of light curves. In addition, common astrophysical phenomena such as chromospheric activity, as well as imperfections in the light curve detrending process, can create systematic artifacts that may produce fictitious signals in the eclipse timings. We present a method to measure precise eclipse times in the presence of distorted light curves, such as in contact and near-contact binaries which exhibit continuously changing light levels in and out of eclipse. 236 systems for which we find a timing variation signal compatible with the presence of a third body are identified. These are modeled for the light time travel effect and the basic properties of the third body are derived. This study complements Orosz et al. (2013; in prep), which focuses on eclipse timing variations of longer period binaries with flat out-of-eclipse regions. Together, these two papers provide comprehensive eclipse timings for all binaries in the Kepler Eclipsing Binary Catalog, as an ongoing resource freely accessible online to the community.

preprint2013arXiv

MARVELS-1: A face-on double-lined binary star masquerading as a resonant planetary system; and consideration of rare false positives in radial velocity planet searches

We have analyzed new and previously published radial velocity observations of MARVELS-1, known to have an ostensibly substellar companion in a ~6- day orbit. We find significant (~100 m/s) residuals to the best-fit model for the companion, and these residuals are naively consistent with an interior giant planet with a P = 1.965d in a nearly perfect 3:1 period commensuribility (|Pb/Pc - 3| < 10^{-4}). We have performed several tests for the reality of such a companion, including a dynamical analysis, a search for photometric variability, and a hunt for contaminating stellar spectra. We find many reasons to be critical of a planetary interpretation, including the fact that most of the three-body dynamical solutions are unstable. We find no evidence for transits, and no evidence of stellar photometric variability. We have discovered two apparent companions to MARVELS-1 with adaptive optics imaging at Keck; both are M dwarfs, one is likely bound, and the other is likely a foreground object. We explore false-alarm scenarios inspired by various curiosities in the data. Ultimately, a line profile and bisector analysis lead us to conclude that the ~100 m/s residuals are an artifact of spectral contamination from a stellar companion contributing ~15-30% of the optical light in the system. We conclude that origin of this contamination is the previously detected radial velocity companion to MARVELS-1, which is not, as previously reported, a brown dwarf, but in fact a G dwarf in a face-on orbit.

preprint2013arXiv

Occultation of the T Tauri Star RW Aurigae A by its Tidally Disrupted Disk

RW Aur A is a classical T Tauri star, believed to have undergone a reconfiguration of its circumstellar environment as a consequence of a recent fly-by of its stellar companion, RW Aur B. This interaction stripped away part of the circumstellar disk of RW Aur A, leaving a tidally disrupted arm and a short truncated circumstellar disk. We present photometric observations of the RW Aur system from the Kilodegree Extremely Little Telescope (KELT) survey showing a long and deep dimming that occurred from September 2010 until March 2011. The dimming has a depth of ~2 magnitudes, a duration of ~180 days and was confirmed by archival observations from American Association of Variable Star Observers (AAVSO). We suggest that this event is the result of a portion of the tidally disrupted disk occulting RW Aur A, specifically a fragment of the tidally disrupted arm. The calculated transverse linear velocity of the occulter is in excellent agreement with the measured relative radial velocity of the tidally disrupted arm. Using simple kinematic and geometric arguments, we show that the occulter cannot be a feature of the RW Aur A circumstellar disk, and we consider and discount other hypotheses. We also place constraints on the thickness and semi-major axis of the portion of the arm that occulted the star.

preprint2013arXiv

Radial Velocity Variations of Photometrically Quiet, Chromospherically Inactive Kepler Stars: A Link Between RV Jitter and Photometric Flicker

We compare stellar photometric variability, as measured from Kepler light curves by Basri et al. (2011), with measurements of radial velocity (RV) root-mean-square (RMS) variations of all California Planet Search overlap stars. We newly derive rotation periods from the Kepler light curves for all of the stars in our study sample. The RV variations reported herein range from less than 4 m/s to 135 m/s, yet the stars all have amplitudes of photometric variability less than 3 mmag, reflecting the preference of the RV program for chromospherically "quiet" stars. Despite the small size of our sample, we find with high statistical significance that the RV RMS manifests strongly in the Fourier power spectrum of the light curve: stars that are noisier in RV have a greater number of frequency components in the light curve. We also find that spot models of the observed light curves systematically underpredict the observed RV variations by factors of ~2--1000, likely because the low level photometric variations in our sample are driven by processes not included in simple spot models. The stars best fit by these models tend to have simpler light curves, dominated by a single relatively high amplitude component of variability. Finally, we demonstrate that the RV RMS behavior of our sample can be explained in the context of the photometric variability evolutionary diagram introduced by Bastien et al. (2013). We use this diagram to derive the surface gravities of the stars in our sample, revealing many of them to have moved off the main-sequence. More generally, we find that the stars with the largest RV RMS are those that have evolved onto the "flicker floor" sequence in that diagram, characterized by relatively low amplitude but highly complex photometric variations which grow as the stars evolve to become subgiants.

preprint2013arXiv

Reanalysis of the Radii of the Benchmark Eclipsing Binary V578 Mon

V578 Mon is an eclipsing binary system in which both stars have masses above 10 Msun determined with an accuracy better than 3%. It is one of only five such massive eclipsing binaries known that also possess eccentric orbits and measured apsidal motions, thus making it an important benchmark for theoretical stellar evolution models. However, recently reported determinations of the radii of V578 Mon differ significantly from previously reported values. We reanalyze the published data for V578 Mon and trace the discrepancy to the use of an incorrect formulation for the stellar potentials in the most recent analysis. Here we report corrected radii for this important benchmark eclipsing binary.

preprint2013arXiv

Stellar Granulation as the Source of High-Frequency Flicker in Kepler Light Curves

A large fraction of cool, low-mass stars exhibit brightness fluctuations that arise from a combination of convective granulation, acoustic oscillations, magnetic activity, and stellar rotation. Much of the short-timescale variability takes the form of stochastic noise, whose presence may limit the progress of extrasolar planet detection and characterization. In order to lay the groundwork for extracting useful information from these quasi-random signals, we focus on the origin of the granulation-driven component of the variability. We apply existing theoretical scaling relations to predict the star-integrated variability amplitudes for 508 stars with photometric light curves measured by the Kepler mission. We also derive an empirical correction factor that aims to account for the suppression of convection in F-dwarf stars with magnetic activity and shallow convection zones. So that we can make predictions of specific observational quantities, we performed Monte Carlo simulations of granulation light curves using a Lorentzian power spectrum. These simulations allowed us to reproduce the so-called "flicker floor" (i.e., a lower bound in the relationship between the full light-curve range and power in short-timescale fluctuations) that was found in the Kepler data. The Monte Carlo model also enabled us to convert the modeled fluctuation variance into a flicker amplitude directly comparable with observations. When the magnetic suppression factor described above is applied, the model reproduces the observed correlation between stellar surface gravity and flicker amplitude. Observationally validated models like these provide new and complementary evidence for a possible impact of magnetic activity on the properties of near-surface convection.

preprint2013arXiv

Very Low Mass Stellar and Substellar Companions to Solar-like Stars From MARVELS IV: A Candidate Brown Dwarf or Low-Mass Stellar Companion to HIP 67526

We report the discovery of a candidate brown dwarf or a very low mass stellar companion (MARVELS-5b) to the star HIP 67526 from the Multi-object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The radial velocity curve for this object contains 31 epochs spread over 2.5 years. Our Keplerian fit using a Markov Chain Monte Carlo approach, reveals that the companion has an orbital period of $90.2695^{+0.0188}_{-0.0187}$ days, an eccentricity of $0.4375 \pm 0.0040$ and a semi-amplitude of $2948.14^{+16.65}_{-16.55}$ m s$^{-1}$. Using additional high-resolution spectroscopy, we find the host star has an effective temperature $T_{\rm{eff}}=6004 \pm 34$ K, a surface gravity $\log g$ [cgs] $=4.55 \pm 0.17$ and a metallicity [Fe/H] $=+0.04 \pm 0.06$. The stellar mass and radius determined through the empirical relationship of Torres et al. (2010), yields 1.10$\pm$0.09 $M_{\sun}$ and 0.92$\pm$0.19 $R_{\sun}$. The minimum mass of MARVELS-5b is $65.0 \pm 2.9 M_{Jup}$, indicating that it is likely to be either a brown dwarf or a very low mass star, thus occupying a relatively sparsely-populated region of the mass function of companions to solar-type stars. The distance to this system is 101$\pm$10 pc from the astrometric measurements of Hipparcos. No stellar tertiary is detected in the high-contrast images taken by either FastCam lucky imaging or Keck adaptive optics imaging, ruling out any star with mass greater than 0.2$M_{\sun}$ at a separation larger than 40 AU.

preprint2013arXiv

Very Low Mass Stellar and Substellar Companions to Solar-Like Stars From MARVELS V: A Low Eccentricity Brown Dwarf from the Driest Part of the Desert, MARVELS-6b

We describe the discovery of a likely brown dwarf (BD) companion with a minimum mass of 31.7 +/- 2.0 M_Jup to GSC 03546-01452 from the MARVELS radial velocity survey, which we designate as MARVELS-6b. For reasonable priors, our analysis gives a probability of 72% that MARVELS-6b has a mass below the hydrogen-burning limit of 0.072 M_Sun, and thus it is a high-confidence BD companion. It has a moderately long orbital period of 47.8929 +0.0063/-0.0062 days with a low eccentricty of 0.1442 +0.0078/-0.0073, and a semi-amplitude of 1644 +12/-13 m/s. Moderate resolution spectroscopy of the host star has determined the following parameters: T_eff = 5598 +/- 63, log g = 4.44 +/- 0.17, and [Fe/H] = +0.40 +/- 0.09. Based upon these measurements, GSC 03546-01452 has a probable mass and radius of M_star = 1.11 +/- 0.11 M_Sun and R_star = 1.06 +/- 0.23 R_Sun with an age consistent with less than ~6 Gyr at a distance of 219 +/- 21 pc from the Sun. Although MARVELS-6b is not observed to transit, we cannot definitively rule out a transiting configuration based on our observations. There is a visual companion detected with Lucky Imaging at 7.7 arcsec from the host star, but our analysis shows that it is not bound to this system. The minimum mass of MARVELS-6b exists at the minimum of the mass functions for both stars and planets, making this a rare object even compared to other BDs.

preprint2012arXiv

An Empirical Correction for Activity Effects on the Temperatures, Radii, and Estimated Masses of Low-Mass Stars and Brown Dwarfs

We present empirical relations for determining the amount by which the effective temperatures and radii -- and therefore the estimated masses -- of low-mass stars and brown dwarfs are altered due to chromospheric activity. We base our relations on a large set of low-mass stars in the field with Halpha activity measurements, and on a set of low-mass eclipsing binaries with X-ray activity measurements from which we indirectly infer the Halpha activity. Both samples yield consistent relations linking the amount by which an active object's temperature is suppressed, and its radius inflated, to the strength of its Halpha emission. These relations are found to approximately preserve bolometric luminosity. We apply these relations to the peculiar brown-dwarf eclipsing binary 2M0535-05, in which the active, higher-mass brown dwarf has a cooler temperature than its inactive, lower-mass companion. The relations correctly reproduce the observed temperatures and radii of 2M0535-05 after accounting for the Halpha emission; 2M0535-05 would be in precise agreement with theoretical isochrones were it inactive. The relations that we present are applicable to brown dwarfs and low-mass stars with masses below 0.8 Msun and for which the activity, as measured by the fractional Halpha luminosity, is in the range -4.6 < log LHa/Lbol < -3.3. We expect these relations to be most useful for correcting radius and mass estimates of low-mass stars and brown dwarfs over their active lifetimes (few Gyr) and when the ages or distances (and therefore luminosities) are unknown. We also discuss the implications of this work for improved determinations of young cluster initial mass functions.

preprint2012arXiv

Characterizing the Cool KOIs III. KOI-961: A Small Star with Large Proper Motion and Three Small Planets

We present the characterization of the star KOI 961, an M dwarf with transit signals indicative of three short-period exoplanets, originally discovered by the Kepler Mission. We proceed by comparing KOI 961 to Barnard's Star, a nearby, well-characterized mid-M dwarf. By comparing colors, optical and near-infrared spectra, we find remarkable agreement between the two, implying similar effective temperatures and metallicities. Both are metal-poor compared to the Solar neighborhood, have low projected rotational velocity, high absolute radial velocity, large proper motion and no quiescent H-alpha emission--all of which is consistent with being old M dwarfs. We combine empirical measurements of Barnard's Star and expectations from evolutionary isochrones to estimate KOI 961's mass (0.13 +/- 0.05 Msun), radius (0.17 +/- 0.04 Rsun) and luminosity (2.40 x 10^(-3.0 +/- 0.3) Lsun). We calculate KOI 961's distance (38.7 +/- 6.3 pc) and space motions, which, like Barnard's Star, are consistent with a high scale-height population in the Milky Way. We perform an independent multi-transit fit to the public Kepler light curve and significantly revise the transit parameters for the three planets. We calculate the false-positive probability for each planet-candidate, and find a less than 1% chance that any one of the transiting signals is due to a background or hierarchical eclipsing binary, validating the planetary nature of the transits. The best-fitting radii for all three planets are less than 1 Rearth, with KOI 961.03 being Mars-sized (Rp = 0.57 +/- 0.18 Rearth), and they represent some of the smallest exoplanets detected to date.

preprint2012arXiv

Correcting for Activity Effects on the Temperatures, Radii, and Estimated Masses of Low-Mass Stars and Brown Dwarfs

We present empirical relations for determining the amount by which the effective temperatures and radii---and therefore the estimated masses---of low-mass stars and brown dwarfs are altered due to chromospheric activity. Accurate estimates of stellar radii are especially important in the context of searches for transiting exoplanets, which rely upon the assumed stellar radius/density to infer the planet radius/density. Our relations are based on a large set of well studied low-mass stars in the field and on a set of benchmark low-mass eclipsing binaries. The relations link the amount by which an active object's temperature is suppressed, and its radius inflated, to the strength of its Halpha emission. These relations are found to approximately preserve bolometric luminosity. We apply these relations to the peculiar brown-dwarf eclipsing binary 2M0535-05, in which the active, higher-mass brown dwarf has a cooler temperature than its inactive, lower-mass companion. The relations correctly reproduce the observed temperatures and radii of 2M0535-05 after accounting for the Halpha emission; 2M0535-05 would be in precise agreement with theoretical isochrones were it inactive. The relations that we present are applicable to brown dwarfs and low-mass stars with masses below 0.8 Msun and for which the activity, as measured by Halpha, is in the range -4.6 < log Lha/Lbol < -3.3. We expect these relations to be most useful for correcting radius and mass estimates of low-mass stars and brown dwarfs over their active lifetimes (few Gyr). We also discuss the implications of this work for determinations of young cluster IMFs.

preprint2012arXiv

Filtergraph: A Flexible Web Application for Instant Data Visualization of Astronomy Datasets

Filtergraph is a web application being developed by the Vanderbilt Initiative in Data-intensive Astrophysics (VIDA) to flexibly handle a large variety of astronomy datasets. While current datasets at Vanderbilt are being used to search for eclipsing binaries and extrasolar planets, this system can be easily reconfigured for a wide variety of data sources. The user loads a flat-file dataset into Filtergraph which instantly generates an interactive data portal that can be easily shared with others. From this portal, the user can immediately generate scatter plots, histograms, and tables based on the dataset. Key features of the portal include the ability to filter the data in real time through user-specified criteria, the ability to select data by dragging on the screen, and the ability to perform arithmetic operations on the data in real time. The application is being optimized for speed in the context of very large datasets: for instance, plot generated from a stellar database of 3.1 million entries render in less than 2 seconds on a standard web server platform. This web application has been created using the Web2py web framework based on the Python programming language. Filtergraph is freely available at http://filtergraph.vanderbilt.edu/.

preprint2012arXiv

High-Resolution Spectroscopy during Eclipse of the Young Substellar Eclipsing Binary 2MASS 0535-0546. II. Secondary Spectrum: No Evidence that Spots Cause the Temperature Reversal

We present high-resolution optical spectra of the young brown-dwarf eclipsing binary 2M0535-05, obtained during eclipse of the higher-mass (primary) brown dwarf. Combined with our previous spectrum of the primary alone (Paper I), the new observations yield the spectrum of the secondary alone. We investigate, through a differential analysis of the two binary components, whether cool surface spots are responsible for suppressing the temperature of the primary. In Paper I, we found a significant discrepancy between the empirical surface gravity of the primary and that inferred via fine analysis of its spectrum. Here we find precisely the same discrepancy in surface gravity, both qualitatively and quantitatively. While this may again be ascribed to either cool spots or model opacity errors, it implies that cool spots cannot be responsible for preferentially lowering the temperature of the primary: if they were, spot effects on the primary spectrum should be preferentially larger, and they are not. The Teff we infer for the primary and secondary, from the TiO-epsilon bands alone, show the same reversal, in the same ratio, as is empirically observed, bolstering the validity of our analysis. In turn, this implies that if suppression of convection by magnetic fields on the primary is the fundamental cause of the Teff reversal, then it cannot be a local suppression yielding spots mainly on the primary (though both components may be equally spotted), but a global suppression in the interior of the primary. We briefly discuss current theories of how this might work.

preprint2012arXiv

KELT-1b: A Strongly Irradiated, Highly Inflated, Short Period, 27 Jupiter-mass Companion Transiting a mid-F Star

We present the discovery of KELT-1b, the first transiting low-mass companion from the wide-field Kilodegree Extremely Little Telescope-North (KELT-North) survey. The V=10.7 primary is a mildly evolved, solar-metallicity, mid-F star. The companion is a low-mass brown dwarf or super-massive planet with mass of 27.23+/-0.50 MJ and radius of 1.110+0.037-0.024 RJ, on a very short period (P=1.21750007) circular orbit. KELT-1b receives a large amount of stellar insolation, with an equilibrium temperature assuming zero albedo and perfect redistribution of 2422 K. Upper limits on the secondary eclipse depth indicate that either the companion must have a non-zero albedo, or it must experience some energy redistribution. Comparison with standard evolutionary models for brown dwarfs suggests that the radius of KELT-1b is significantly inflated. Adaptive optics imaging reveals a candidate stellar companion to KELT-1, which is consistent with an M dwarf if bound. The projected spin-orbit alignment angle is consistent with zero stellar obliquity, and the vsini of the primary is consistent with tidal synchronization. Given the extreme parameters of the KELT-1 system, we expect it to provide an important testbed for theories of the emplacement and evolution of short-period companions, and theories of tidal dissipation and irradiated brown dwarf atmospheres.

preprint2012arXiv

KELT-2Ab: A Hot Jupiter Transiting the Bright (V=8.77) Primary Star of a Binary System

We report the discovery of KELT-2Ab, a hot Jupiter transiting the bright (V=8.77) primary star of the HD 42176 binary system. The host is a slightly evolved late F-star likely in the very short-lived "blue-hook" stage of evolution, with $\teff=6148\pm48{\rm K}$, $\log{g}=4.030_{-0.026}^{+0.015}$ and $\feh=0.034\pm0.78$. The inferred stellar mass is $M_*=1.314_{-0.060}^{+0.063}$\msun\ and the star has a relatively large radius of $R_*=1.836_{-0.046}^{+0.066}$\rsun. The planet is a typical hot Jupiter with period $4.11379\pm0.00001$ days and a mass of $M_P=1.524\pm0.088$\mj\ and radius of $R_P=1.290_{-0.050}^{+0.064}$\rj. This is mildly inflated as compared to models of irradiated giant planets at the $\sim$4 Gyr age of the system. KELT-2A is the third brightest star with a transiting planet identified by ground-based transit surveys, and the ninth brightest star overall with a transiting planet. KELT-2Ab's mass and radius are unique among the subset of planets with $V<9$ host stars, and therefore increases the diversity of bright benchmark systems. We also measure the relative motion of KELT-2A and -2B over a baseline of 38 years, robustly demonstrating for the first time that the stars are bound. This allows us to infer that KELT-2B is an early K-dwarf. We hypothesize that through the eccentric Kozai mechanism KELT-2B may have emplaced KELT-2Ab in its current orbit. This scenario is potentially testable with Rossiter-McLaughlin measurements, which should have an amplitude of $\sim$44 m s$^{-1}$.

preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

preprint2012arXiv

The SLoWPoKES Catalog of Low-mass Ultra-wide Binaries: A Cool Stars Resource for Testing Fundamental Properties and for Constraining Binary Formation Theory

We present results from the Sloan Low-mass Wide Pairs of Kinematically Equivalent Stars (SLoWPoKES) catalog of ultra-wide (10^3-10^5.5 AU), low-mass (K5-M7) common proper motion binaries. We constructed a Galactic model, based on empirical stellar number density and 3D velocity distributions, to select bona fide pairs with probability of chance alignment <5%, making SLoWPoKES an efficient sample for followup observations. Our initial catalog contains 1342 disk dwarf, subdwarf, and white dwarf-red dwarf systems and is the largest collection of low-mass, wide binaries ever assembled. The diversity---in mass, metallicity, age, and evolutionary states---of SLoWPoKES pairs makes it a valuable resource of coeval laboratories to examine and constrain the physical properties of low-mass stars. SLoWPoKES pairs show signatures of two (or more) formation modes in the distribution of the physical separation and higher-order multiplicity. Neither dynamical dissipation of primordial triples/quadruples or dynamical capture of ejected stars can explain the observed populations by itself. We use followup spectroscopic observations to recalibrate the metallicity-sensitive ζ_(TiO/CaH) index by assuming that both members of the binary system have the same composition. Our new formulation is a significantly better tracer of absolute metallicity, particularly for the early-type M dwarfs. The SLoWPoKES catalog is publicly available on a custom data visualization portal.

preprint2012arXiv

Very Low-Mass Stellar and Substellar Companions to Solar-Like Stars from MARVELS I: A Low Mass Ratio Stellar Companion to TYC 4110-01037-1 in a 79-day Orbit

TYC 4110-01037-1 has a low-mass stellar companion, whose small mass ratio and short orbital period are atypical amongst solar-like (Teff ~< 6000 K) binary systems. Our analysis of TYC 4110-01037-1 reveals it to be a moderately aged (~<5 Gyr) solar-like star having a mass of 1.07 +/- 0.08 MSun and radius of 0.99 +/- 0.18 RSun. We analyze 32 radial velocity measurements from the SDSS-III MARVELS survey as well as 6 supporting radial velocity measurements from the SARG spectrograph on the 3.6m TNG telescope obtained over a period of ~2 years. The best Keplerian orbital fit parameters were found to have a period of 78.994 +/- 0.012 days, an eccentricity of 0.1095 +/- 0.0023, and a semi-amplitude of 4199 +/- 11 m/s. We determine the minimum companion mass (if sin i = 1) to be 97.7 +/- 5.8 MJup. The system's companion to host star mass ratio, >0.087 +/- 0.003, places it at the lowest end of observed values for short period stellar companions to solar-like (Teff ~< 6000 K) stars. One possible way to create such a system would be if a triple-component stellar multiple broke up into a short period, low q binary during the cluster dispersal phase of its lifetime. A candidate tertiary body has been identified in the system via single-epoch, high contrast imagery. If this object is confirmed to be co-moving, we estimate it would be a dM4 star. We present these results in the context of our larger-scale effort to constrain the statistics of low mass stellar and brown dwarf companions to FGK-type stars via the MARVELS survey.

preprint2012arXiv

Very Low-mass Stellar and Substellar Companions to Solar-like Stars from Marvels III: A Short-Period Brown Dwarf Candidate Around An Active G0Iv Subgiant

We present an eccentric, short-period brown dwarf candidate orbiting the active, slightly evolved subgiant star TYC 2087-00255-1, which has effective temperature T_eff = 5903+/-42 K, surface gravity log (g) = 4.07+/-0.16 (cgs), and metallicity [Fe/H] = -0.23+/-0.07. This candidate was discovered using data from the first two years of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of Sloan Digital Sky Survey. From our 38 radial velocity measurements spread over a two-year time baseline, we derive a Keplerian orbital fit with semi-amplitude K=3.571+/-0.041 km/s, period P=9.0090+/-0.0004 days, and eccentricity e=0.226+/-0.011. Adopting a mass of 1.16+/-0.11 Msun for the subgiant host star, we infer that the companion has a minimum mass of 40.0+/-2.5 M_Jup. Assuming an edge-on orbit, the semimajor axis is 0.090+/-0.003 AU. The host star is photometrically variable at the \sim1% level with a period of \sim13.16+/-0.01 days, indicating that the host star spin and companion orbit are not synchronized. Through adaptive optics imaging we also found a point source 643+/-10 mas away from TYC 2087-00255-1, which would have a mass of 0.13 Msun if it is physically associated with TYC 2087-00255-1 and has the same age. Future proper motion observation should be able to resolve if this tertiary object is physically associated with TYC 2087-00255-1 and make TYC 2087-00255-1 a triple body system. Core Ca II H and K line emission indicate that the host is chromospherically active, at a level that is consistent with the inferred spin period and measured v_{rot}*sin i, but unusual for a subgiant of this T_eff. This activity could be explained by ongoing tidal spin-up of the host star by the companion.

preprint2011arXiv

A Chandra ACIS Study of the Young Star Cluster Trumpler 15 in Carina and Correlation with Near-infrared Sources

Using the highest-resolution X-ray observation of the Trumpler 15 star cluster taken by the Chandra X-ray Observatory, we estimate the total size of its stellar population by comparing the X-ray luminosity function of the detected sources to a calibrator cluster, and identify for the first time a significant fraction (~14%) of its individual members. The highest-resolution near-IR observation of Trumpler 15 (taken by the HAWK-I instrument on the VLT) was found to detect most of our X-ray selected sample of cluster members, with a K-excess disk frequency of 3.8+-0.7%. The near-IR data, X-ray luminosity function, and published spectral types of the brightest members support a cluster age estimate (5-10 Myr) that is older than those for the nearby Trumpler 14 and Trumpler 16 clusters, and suggest that high-mass members may have already exploded as supernovae. The morphology of the inner ~0.7 pc core of the cluster is found to be spherical. However, the outer regions (beyond 2 pc) are elongated, forming an `envelope' of stars that, in projection, appears to connect Trumpler 15 to Trumpler 14; this morphology supports the view that these clusters are physically associated. Clear evidence of mass segregation is seen. This study appears in a Special Issue of the ApJS devoted to the Chandra Carina Complex Project (CCCP), a 1.42 square degree Chandra X-ray survey of the Great Nebula in Carina.

preprint2011arXiv

A Pan-Carina YSO Catalog: Intermediate-Mass Young Stellar Objects in the Carina Nebula Identified Via Mid-Infrared Excess Emission

We present a catalog of 1439 young stellar objects (YSOs) spanning the 1.42 deg^2 field surveyed by the Chandra Carina Complex Project (CCCP), which includes the major ionizing clusters and the most active sites of ongoing star formation within the Great Nebula in Carina. Candidate YSOs were identified via infrared (IR) excess emission from dusty circumstellar disks and envelopes, using data from the Spitzer Space Telescope Vela--Carina survey and the Two-Micron All Sky Survey. We model the 1--24 /mu m IR spectral energy distributions of the YSOs to constrain physical properties. Our Pan-Carina YSO Catalog (PCYC) is dominated by intermediate-mass (2 Msun < m < 10 Msun) objects with disks, including Herbig Ae/Be stars and their less evolved progenitors. The PCYC provides a valuable complementary dataset to the CCCP X-ray source catalogs, identifying 1029 YSOs in Carina with no X-ray detection. We also catalog 410 YSOs with X-ray counterparts, including 62 candidate protostars. Candidate protostars with X-ray detections tend to be more evolved than those without. In most cases, X-ray emission apparently originating from intermediate-mass, disk-dominated YSOs is consistent with the presence of low-mass companions, but we also find that X-ray emission correlates with cooler stellar photospheres and higher disk masses. We suggest that intermediate-mass YSOs produce X-rays during their early pre-main sequence evolution, perhaps driven by magnetic dynamo activity during the convective atmosphere phase, but this emission dies off as the stars approach the main sequence. Extrapolating over the stellar initial mass function scaled to the PCYC population, we predict a total population of >2x10^4 YSOs and a present-day star formation rate (SFR) of >0.008 Msun/yr. The global SFR in the Carina Nebula, averaged over the past ~5 Myr, has been approximately constant.

preprint2011arXiv

Apsidal Motion of the Massive, Benchmark Eclipsing Binary V578 Mon

V578 Mon is a system of two early B-type stars in the Rosette Nebula star-forming region (NGC 2244), and is one of only nine eclipsing binaries with component masses greater than 10 M\odot whose physical parameters have been determined with an accuracy of better than 3%. It is therefore a benchmark system for evolutionary and stellar structure models of newly formed massive stars. Combining our multi-band light curves spanning 40 yr with previous light curve data from the literature, we fit a model light curve that for the first time includes the effects of apsidal motion of the system. We measure an apsidal period of 33.48+0.10-0.06 yr. As a consequence of incorporating the apsidal motion into the modeling of the system's orbital parameters, we determine an updated eccentricity of e = 0.07755+0.00022-0.00027, which differs significantly from the value previously reported in the literature. Evidently, the inclusion of apsidal motion in the light curve modeling significantly affects the eccentricity determination. Incorporating these key parameters into a comprehensive model of the system's physical parameters-including internal structure constraints- will bring V578 Mon to the next level of benchmark precision and utility.

preprint2011arXiv

Eclipsing Binary Science Via the Merging of Transit and Doppler Exoplanet Survey Data - A Case Study With the MARVELS Pilot Project and SuperWASP

Exoplanet transit and Doppler surveys discover many binary stars during their operation that can be used to conduct a variety of ancillary science. Specifically, eclipsing binary stars can be used to study the stellar mass-radius relationship and to test predictions of theoretical stellar evolution models. By cross-referencing 24 binary stars found in the MARVELS Pilot Project with SuperWASP photometry, we find two new eclipsing binaries, TYC 0272-00458-1 and TYC 1422-01328-1, which we use as case studies to develop a general approach to eclipsing binaries in survey data. TYC 0272-00458-1 is a single-lined spectroscopic binary for which we calculate a mass of the secondary and radii for both components using reasonable constraints on the primary mass through several different techniques. For a primary mass of M_1 = 0.92 +/- 0.1 M_solar, we find M_2 = 0.610 +/- 0.036 M_solar, R_1 = 0.932 +/- 0.076 R_solar and R_2 = 0.559 +/- 0.102 R_solar, and find that both stars have masses and radii consistent with model predictions. TYC 1422-01328-1 is a triple-component system for which we can directly measure the masses and radii of the eclipsing pair. We find that the eclipsing pair consists of an evolved primary star (M_1 = 1.163 +/- 0.034 M_solar, R_1 = 2.063 +/- 0.058 R_solar) and a G-type dwarf secondary (M_2 = 0.905 +/- 0.067 M_solar, R_2 = 0.887 +/- 0.037 R_solar). We provide the framework necessary to apply this analysis to much larger datasets.

preprint2011arXiv

Expected Large Synoptic Survey Telescope (LSST) Yield of Eclipsing Binary Stars

In this paper we estimate the Large Synoptic Survey Telescope (LSST) yield of eclipsing binary stars, which will survey ~20,000 square degrees of the southern sky during the period of 10 years in 6 photometric passbands to r ~ 24.5. We generate a set of 10,000 eclipsing binary light curves sampled to the LSST time cadence across the whole sky, with added noise as a function of apparent magnitude. This set is passed to the Analysis of Variance (AoV) period finder to assess the recoverability rate for the periods, and the successfully phased light curves are passed to the artificial intelligence-based pipeline EBAI to assess the recoverability rate in terms of the eclipsing binaries' physical and geometric parameters. We find that, out of ~24 million eclipsing binaries observed by LSST with S/N>10 in mission life-time, ~28% or 6.7 million can be fully characterized by the pipeline. Of those, ~25% or 1.7 million will be double-lined binaries, a true treasure trove for stellar astrophysics.

preprint2011arXiv

High Cadence Time-Series Photometry of V1647 Orionis

We present high cadence (1-10 hr^-1) time-series photometry of the eruptive young variable star V1647 Orionis during its 2003-2004 and 2008-2009 outbursts. The 2003 light curve was obtained mid-outburst at the phase of steepest luminosity increase of the system, during which time the accretion rate of the system was presumably continuing to increase toward its maximum rate. The 2009 light curve was obtained after the system luminosity had plateaued, presumably when the rate of accretion had also plateaued. We detect a 'flicker noise' signature in the power spectrum of the lightcurves, which may suggest that the stellar magnetosphere continued to interact with the accretion disk during each outburst event. Only the 2003 power spectrum, however, evinces a significant signal with a period of 0.13 d. While the 0.13 d period cannot be attributed to the stellar rotation period, we show that it may plausibly be due to short-lived radial oscillations of the star, possibly caused by the surge in the accretion rate.

preprint2011arXiv

Refined Metallicity Indices for M Dwarfs Using the SLoWPoKES Catalog of Wide, Low-mass Binaries

We report the results from spectroscopic observations of 113 ultra-wide, low-mass binary systems, composed largely of M0--M3 dwarfs, from the SLoWPoKES catalog of common proper motion pairs identified in the Sloan Digital Sky Survey. Radial velocities of each binary member were used to confirm that they are co-moving and, consequently, to further validate the high fidelity of the SLoWPoKES catalog. Ten stars appear to be spectroscopic binaries based on broad or split spectral features, supporting previous findings that wide binaries are likely to be hierarchical systems. We measured the Hα equivalent width of the stars in our sample and found that components of 81% of the observed pairs has similar Hα levels. The difference in Hα equivalent width amongst components with similar masses was smaller than the range of Hα variability for individual objects. We confirm that the Lepine et al. ζ(CaH2+CaH3, TiO5) index traces iso-metallicity loci for most of our sample of M dwarfs. However, we find a small systematic bias in ζ, especially in the early-type M dwarfs. We use our sample to recalibrate the definition of ζ. While representing a small change in the definition, the new ζ is a significantly better predictor of iso-metallicity for the higher mass M dwarfs.

preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

preprint2011arXiv

Spectral Energy Distributions of Young Stars in IC 348: The Role of Disks in Angular Momentum Evolution of Young, Low-Mass Stars

Theoretical work suggests that a young star's angular momentum and rotation rate may be strongly influenced by magnetic interactions with its circumstellar disk. A generic prediction of these 'disk-locking' (DL) theories is that a disk-locked star will be forced to co-rotate with the Keplerian angular velocity of the inner edge of the disk. These theories have also been interpreted to suggest a correlation between young stars' rotation periods and the structural properties of their disks, such that slowly rotating stars possess close-in disks that enforce the star's slow rotation, whereas rapidly rotating stars possess anemic or evacuated inner disks that are unable to brake the stars and they spin up as they contract. To test these expectations, we model the SEDs of 33 young stars in IC 348 with known rotation periods and infrared excesses indicating the presence of disks. For each star, we match the observed spectral energy distribution, typically sampling 0.6-8.0 \mum, to a grid of 200,000 pre-computed star+disk radiative transfer models, from which we infer the disk's inner-truncation radius (R_trunc). We then compare this R_trunc to the disk's co-rotation radius (R_co), calculated from the star's rotation period. We do not find obvious differences in the disk R_trunc of slow vs. rapid rotators. This holds true both at the level of whether close-in disk material is present at all, and in analyzing the precise location of the inner disk edge relative to the R_co amongst the subset of stars with close-in disk material. One interpretation is that DL is unimportant for the IC 348 stars in our sample. Alternatively, if DL does operate, then it must operate on both the slow and rapid rotators, potentially producing both spin-up and spin-down torques, and the transition from the disk-locked state to the disk-released state must occur more rapidly than the stellar contraction timescale.

preprint2011arXiv

The Chandra Carina Complex Project View of Trumpler 16

Trumpler 16 is a well--known rich star cluster containing the eruptive supergiant $η$ Carinæ and located in the Carina star-forming complex. In the context of the Chandra Carina Complex Project, we study Trumpler 16 using new and archival X-ray data. A revised X-ray source list of the Trumpler 16 region contains 1232 X-ray sources including 1187 likely Carina members. These are matched to 1047 near-infrared counterparts detected by the HAWK-I instrument at the VLT allowing for better selection of cluster members. The cluster is irregular in shape. Although it is roughly circular, there is a high degree of sub-clustering, no noticeable central concentration and an extension to the southeast. The high--mass stars show neither evidence of mass segregation nor evidence of strong differential extinction. The derived power-law slope of the X-ray luminosity function for Trumpler 16 reveals a much steeper function than the Orion Nebula Cluster implying different ratio of solar- to higher-mass stars. We estimate the total Trumpler 16 pre-main sequence population to be > 6500 Class II and Class III X-ray sources. An overall K-excess disk frequency of ~ 8.9% is derived using the X-ray selected sample, although there is some variation among the sub-clusters, especially in the Southeastern extension. X-ray emission is detected from 29 high--mass stars with spectral types between B2 and O3.

preprint2011arXiv

The Color-Period Diagram and Stellar Rotational Evolution - New Rotation Period Measurements in the Open Cluster M34

We present results from a 5-month photometric survey for stellar rotation periods combined with a 4-year radial-velocity survey for membership and binarity in the 220Myr open cluster M34. We report surface rotation periods for 120 stars, 83 of which are late-type cluster members. A comparison to previous work serves to illustrate the importance of high cadence long baseline photometric observations and membership information. The new M34 periods are less biased against slow rotation and cleaned for non-members. The rotation periods of the cluster members span more than an order of magnitude from 0.5 day up to 11.5 days, and trace two distinct rotational sequences - fast (C) and moderate-to-slow (I) - in the color-period diagram. The sequences represent two different states in the rotational evolution of the late-type cluster members. We use the color-period diagrams for M34 and for younger and older clusters to estimate the timescale for the transition from the C to the I sequence and find ~<150Myr, ~150-300Myr, and ~300-600Myr for G, early-mid K, and late K dwarfs, respectively. The small number of stars in the gap between C and I suggest a quick transition. We estimate a lower limit on the maximum spin-down rate (dP/dt) during this transition to be ~0.06 days/Myr and ~0.08 days/Myr for early and late K dwarfs, respectively. We compare the I sequence rotation periods in M34 and the Hyades for G and K dwarfs and find that K dwarfs spin down slower than the Skumanich rate. We determine a gyrochronology age of 240Myr for M34. We measure the effect of cluster age uncertainties on the gyrochronology age for M34 and find the resulting error to be consistent with the error estimate for the technique. We use the M34 I sequence to redetermine the coefficients in the expression for rotational dependence on color used in gyrochronology (abridged).

preprint2011arXiv

The Initial Mass Function of the Orion Nebula Cluster across the H-burning limit

We present a new census of the Orion Nebula Cluster (ONC) over a large field of view (>30'x30'), significantly increasing the known population of stellar and substellar cluster members with precisely determined properties. We develop and exploit a technique to determine stellar effective temperatures from optical colors, nearly doubling the previously available number of objects with effective temperature determinations in this benchmark cluster. Our technique utilizes colors from deep photometry in the I-band and in two medium-band filters at lambda~753 and 770nm, which accurately measure the depth of a molecular feature present in the spectra of cool stars. From these colors we can derive effective temperatures with a precision corresponding to better than one-half spectral subtype, and importantly this precision is independent of the extinction to the individual stars. Also, because this technique utilizes only photometry redward of 750nm, the results are only mildly sensitive to optical veiling produced by accretion. Completing our census with previously available data, we place some 1750 sources in the Hertzsprung-Russel diagram and assign masses and ages down to 0.02 solar masses. At faint luminosities, we detect a large population of background sources which is easily separated in our photometry from the bona fide cluster members. The resulting initial mass function of the cluster has good completeness well into the substellar mass range, and we find that it declines steeply with decreasing mass. This suggests a deficiency of newly formed brown dwarfs in the cluster compared to the Galactic disk population.

preprint2011arXiv

Time-Series Photometry of Stars in and around the Lagoon Nebula. I. Rotation Periods of 290 Low-Mass Pre-Main-Sequence Stars in NGC 6530

We have conducted a long-term, wide-field, high-cadence photometric monitoring survey of ~50,000 stars in the Lagoon Nebula \ion{H}{2} region. This first paper presents rotation periods for 290 low-mass stars in NGC 6530, the young cluster illuminating the nebula, and for which we assemble a catalog of infrared and spectroscopic disk indicators, estimated masses and ages, and X-ray luminosities. The distribution of rotation periods we measure is broadly uniform for 0.5 < P < 10 d; the short-period cutoff corresponds to breakup. We observe no obvious bimodality in the period distribution, but we do find that stars with disk signatures rotate more slowly on average. The stars' X-ray luminosities are roughly flat with rotation period, at the saturation level ($\log L_X / L_{\rm bol} \approx -3.3$). However, we find a significant positive correlation between $L_X / L_{\rm bol}$ and co-rotation radius, suggesting that the observed X-ray luminosities are regulated by centrifugal stripping of the stellar coronae. The period-mass relationship in NGC 6530 is broadly similar to that of the Orion Nebula Cluster (ONC), but the slope of the relationship among the slowest rotators differs from that in the ONC and other young clusters. We show that the slope of the period-mass relationship for the slowest rotators can be used as a proxy for the age of a young cluster, and we argue that NGC 6530 may be slightly younger than the ONC, making it a particularly important touchstone for models of angular momentum evolution in young, low-mass stars.

preprint2011arXiv

X-ray Star Clusters in the Carina Complex

The distribution of young stars found in the Chandra Carina Complex Project (CCCP) is examined for clustering structure. X-ray surveys are advantageous for identifying young stellar populations compared to optical and infrared surveys in suffering less contamination from nebular emission and Galactic field stars. The analysis is based on smoothed maps of a spatially complete subsample of about 3000 brighter X-ray sources classified as Carina members, and about 10,000 stars from the full CCCP sample. The principal known clusters are recovered, and some additional smaller groups are identified. No rich embedded clusters are present, although a number of sparse groups are found. The CCCP reveals considerable complexity in clustering properties. The Trumpler 14 and 15 clusters have rich stellar populations in unimodal, centrally concentrated structures several parsecs across. Non-spherical internal structure is seen, and large-scale low surface density distributions surround these rich clusters. Trumpler 16, in contrast, is comprised of several smaller clusters within a circular boundary. Collinder 228 is a third type of cluster which extends over tens of parsecs with many sparse compact groups likely arising from triggered star formation processes. A widely dispersed, but highly populous, distribution of X-ray stars across the about 50 pc CCCP mosaic supports a model of past generations of star formation in the region. Collinder 234, a group of massive stars without an associated cluster of pre-main sequence stars, may be part of this dispersed population.

preprint2010arXiv

A Multi-color Optical Survey of the Orion Nebula Cluster. II. the H-R diagram

We present a new analysis of the stellar population of the Orion Nebula Cluster (ONC) based on multi-band optical photometry and spectroscopy. We study the color-color diagrams in BVI, plus a narrow-band filter centered at 6200A, finding evidences that intrinsic color scales valid for main-sequence dwarfs are incompatible with the ONC, while a better agreement is found employing synthetic intrinsic colors obtained constraining the typical lower surface gravity of young stars. We refine these model colors even further, empirically, by comparison with a sample of ONC stars with no accretion and no extinction. We consider the stars with known spectral types from the literature, and add 65 newly classified stars from slit spectroscopy and 182 M-type from narrow-band photometry; in this way we isolate a sample of about 1000 stars with known spectral type. We introduce a new method to self-consistently derive reddening and accretion excess from the location of each star in the BVI color-color diagram. This enables us to accurately determine the extinction of the ONC members. We adopt a lower distance for the ONC than previously assumed, based on recent parallax measurements. With a careful choice also of the spectral type-temperature transformation, we produce the new H-R diagram of the ONC population, more populated than previous works. With respect to previous works, we find higher luminosity for late-type stars and a lower luminosity for early types. We determine the age distribution of the population, peaking at 2-3 Myr, a higher age than previously estimated. We study the distribution of the members in the mass-age plane, and find that taking into account selection effects due to incompleteness removes an apparent correlation between mass and age. We derive the IMF for low- and intermediate-mass members of the ONC, which turns out to be model-dependent, and shows a turn-over at ~<0.2Msun.

preprint2010arXiv

A Search for Star-Disk Interaction Among the Strongest X-ray Flaring Stars in the Orion Nebula Cluster

The Chandra Orion Ultradeep Project observed hundreds of young, low-mass stars undergoing highly energetic X-ray flare events. The 32 most powerful cases have been modeled with the result that the magnetic structures responsible for these flares can be many stellar radii in extent. In this paper, we model the observed spectral energy distributions of these 32 stars in order to determine, in detail for each star, whether there is circumstellar disk material situated in sufficient proximity to the stellar surface for interaction with the large magnetic loops inferred from the observed X-ray flares. Our spectral energy distributions span the wavelength range 0.3-8 um (plus 24 um for some stars), allowing us to constrain the presence of dusty circumstellar material out to >10 AU from the stellar surface in most cases. For 24 of the 32 stars in our sample the available data are sufficient to constrain the location of the inner edge of the dusty disks. Six of these (25%) have spectral energy distributions consistent with inner disks within reach of the observed magnetic loops. Another four stars may have gas disks interior to the dust disk and extending within reach of the magnetic loops, but we cannot confirm this with the available data. The remaining 14 stars (58%) appear to have no significant disk material within reach of the large flaring loops. Thus, up to ~40% of the sample stars exhibit energetic X-ray flares that possibly arise from a magnetic star-disk interaction, and the remainder are evidently associated with extremely large, free-standing magnetic loops anchored only to the stellar surface.

preprint2010arXiv

High Resolution Spectroscopy during Eclipse of the Young Substellar Eclipsing Binary 2MASS 0535-0546. I. Primary Spectrum: Cool Spots versus Opacity Uncertainties

We present high-resolution Keck optical spectra of the very young substellar eclipsing binary 2MASS J05352184-0546085, obtained during eclipse of the lower-mass (secondary) brown dwarf. The observations yield the spectrum of the higher-mass (primary) brown dwarf alone, with negligible (~1.6%) contamination by the secondary. We perform a simultaneous fine-analysis of the TiO-epsilon band and the red lobe of the KI doublet, using state-of-the-art PHOENIX Dusty and Cond synthetic spectra. Comparing the effective temperature and surface gravity derived from these fits to the {\it empirically} determined surface gravity of the primary (logg=3.5) then allows us to test the model spectra as well as probe the prevailing photospheric conditions. We find that: (1) fits to TiO-epsilon alone imply Teff=2500 \pm 50K; (2) at this Teff, fits to KI imply logg=3.0, 0.5 dex lower than the true value; and (3) at the true logg, KI fits yield Teff=2650 \pm 50K, ~150K higher than from TiO-epsilon alone. On the one hand, these are the trends expected in the presence of cool spots covering a large fraction of the primary's surface (as theorized previously to explain the observed Teff reversal between the primary and secondary). Specifically, our results can be reproduced by an unspotted stellar photosphere with Teff=2700K and (empirical) logg=3.5, coupled with axisymmetric cool spots that are 15% cooler (2300K), have an effective logg=3.0 (0.5 dex lower than photospheric), and cover 70% of the surface. On the other hand, the trends in our analysis can also be reproduced by model opacity errors: there are lacks in the synthetic TiO-epsilon opacities, at least for higher-gravity field dwarfs. Stringently discriminating between the two possibilities requires combining the present results with an equivalent analysis of the secondary (predicted to be relatively unspotted compared to the primary).

preprint2010arXiv

MARVELS-1b: A Short-Period, Brown Dwarf Desert Candidate from the SDSS-III MARVELS Planet Search

We present a new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup.

preprint2010arXiv

Resolved Spectroscopy of M Dwarf/L Dwarf Binaries. IV. Discovery of an M9 + L6 BInary Separated by Over 100 AU

We report the discovery of a faint L6 \pm 1 companion to the previously known M9 dwarf, 2MASS J01303563-4445411, based on our near-infrared imaging and spectroscopic observations with the 3m Infrared Telescope Facility SpeX imager/spectrometer. The visual binary is separated by 3. 28 \pm 0. 05 on the sky at a spectrophotometric distance of 40 \pm 14 pc. The projected physical separation is 130 \pm 50 AU, making it one of the widest VLM field multiples containing a brown dwarf companion. 2MASS J0130-4445 is only one of ten wide VLM pairs and only one of six in the field. The secondary is considerably fainter (ΔK ~ 2.35 mag) and redder (Δ (J - Ks) ~ 0.81 dex), consistent with component near-infrared types of M9.0 \pm 0.5 and L6 \pm 1 based on our resolved spectroscopy. The component types suggest a secondary mass well within the hydrogen-burning limit and an age-dependent mass ratio of 0.6-0.9. The system's space motion and spectroscopic indicators suggest an age of 2-4 Gyr while the model-dependent masses and binding energies suggest that this system is unlikely to have formed via dynamical ejection. The age, composition, and separation of the 2MASS J01303563-4445411 system make it useful for tests of VLM formation theories and of condensate cloud formation in L dwarfs.

preprint2010arXiv

Sloan Low-Mass Wide Pairs of Kinematically Equivalent Stars (SLoWPoKES): A Catalog of Very Wide, Low-mass Pairs

We present the Sloan Low-mass Wide Pairs of Kinematically Equivalent Stars (SLoWPoKES), a catalog of 1342 very-wide (projected separation >500 AU), low-mass (at least one mid-K--mid-M dwarf component) common proper motion pairs identified from astrometry, photometry, and proper motions in the Sloan Digital Sky Survey. A Monte Carlo based Galactic model is constructed to assess the probability of chance alignment for each pair; only pairs with a probability of chance alignment </= 0.05 are included in the catalog. The overall fidelity of the catalog is expected to be 98.35%. The selection algorithm is purposely exclusive to ensure that the resulting catalog is efficient for follow-up studies of low-mass pairs. The SLoWPoKES catalog is the largest sample of wide, low-mass pairs to date and is intended as an ongoing community resource for detailed study of bona fide systems. Here we summarize the general characteristics of the SLoWPoKES sample and present preliminary results describing the properties of wide, low-mass pairs. While the majority of the identified pairs are disk dwarfs, there are 70 halo subdwarf pairs and 21 white dwarf-disk dwarf pairs, as well as four triples. Most SLoWPoKES pairs violate the previously defined empirical limits for maximum angular separation or binding energies. However, they are well within the theoretical limits and should prove very useful in putting firm constraints on the maximum size of binary systems and on different formation scenarios. We find a lower limit to the wide binary frequency for the mid-K-mid-M spectral types that constitute our sample to be 1.1%. This frequency decreases as a function of Galactic height, indicating a time evolution of the wide binary frequency. [See text for full abstract.]

preprint2010arXiv

Solar Flares and Coronal Mass Ejections: A Statistically Determined Flare Flux-CME Mass Correlation

In an effort to examine the relationship between flare flux and corresponding CME mass, we temporally and spatially correlate all X-ray flares and CMEs in the LASCO and GOES archives from 1996 to 2006. We cross-reference 6,733 CMEs having well-measured masses against 12,050 X-ray flares having position information as determined from their optical counterparts. For a given flare, we search in time for CMEs which occur 10-80 minutes afterward, and we further require the flare and CME to occur within +/-45 degrees in position angle on the solar disk. There are 826 CME/flare pairs which fit these criteria. Comparing the flare fluxes with CME masses of these paired events, we find CME mass increases with flare flux, following an approximately log-linear, broken relationship: in the limit of lower flare fluxes, log(CME mass)~0.68*log(flare flux), and in the limit of higher flare fluxes, log(CME mass)~0.33*log(flare flux). We show that this broken power-law, and in particular the flatter slope at higher flare fluxes, may be due to an observational bias against CMEs associated with the most energetic flares: halo CMEs. Correcting for this bias yields a single power-law relationship of the form log(CME mass)~0.70*log(flare flux). This function describes the relationship between CME mass and flare flux over at least 3 dex in flare flux, from ~10^-7 to 10^-4 W m^-2.

preprint2010arXiv

Spitzer Space Telescope observations of the Carina Nebula: The steady march of feedback-driven star formation

We report the first results of imaging the Carina Nebula with Spitzer/IRAC, providing a catalog of point sources and YSOs based on SED fits. We discuss several aspects of the extended emission, including dust pillars that result when a clumpy molecular cloud is shredded by massive star feedback. There are few "extended green objects" (EGOs) normally taken as signposts of outflow activity, and none of the HH jets detected optically are seen as EGOs. A population of "extended red objects" tends to be found around OB stars, some with clear bow-shocks. These are dusty shocks where stellar winds collide with flows off nearby clouds. Finally, the relative distributions of O stars and subclusters of YSOs as compared to dust pillars shows that while some YSOs are located within pillars, many more stars and YSOs reside just outside pillar heads. We suggest that pillars are transient phenomena, part of a continuous outwardly propagating wave of star formation driven by massive star feedback. As pillars are destroyed, they leave newly formed stars in their wake, which are then subsumed into the young OB association. Altogether, the current generation of YSOs shows no strong deviation from a normal IMF. The number of YSOs suggests a roughly constant star-formation rate over the past 3Myr, implying that star formation in pillars constitutes an important mechanism to construct unbound OB associations. Accelerated pillars may give birth to O-type stars that, after several Myr, could appear to have formed in isolation.

preprint2010arXiv

The Exotic Eclipsing Nucleus of the Ring Planetary Nebula SuWt2

SuWt2 is a planetary nebula (PN) consisting of a bright ionized thin ring seen nearly edge-on. It has a bright (V=12) central star, too cool to ionize the PN, which we discovered to be an eclipsing binary. A spectrum from IUE did not reveal a UV source. We present extensive ground-based photometry and spectroscopy of the central binary collected over the ensuing two decades, resulting in the determination that the orbital period of the eclipsing pair is 4.9 d, and consists of two nearly identical A1 V stars, each of mass ~2.7 M_sun. The physical parameters of the A stars, combined with evolutionary tracks, show that both are in the short-lived "blue-hook" evolutionary phase that occurs between the main sequence and the Hertzsprung gap, and that the age of the system is about 520 Myr. One puzzle is that the stars' rotational velocities are different from each other, and considerably slower than synchronous with the orbital period. It is possible that the center-of-mass velocity of the eclipsing pair is varying with time, suggesting that there is an unseen third orbiting body in the system. We propose a scenario in which the system began as a hierarchical triple, consisting of a ~2.9 M_sun star orbiting the close pair of A stars. Upon reaching the AGB stage, the primary engulfed the pair into a common envelope, leading to a rapid contraction of the orbit and catastrophic ejection of the envelope into the orbital plane. In this picture, the exposed core of the initial primary is now a white dwarf of ~0.7 M_sun, orbiting the eclipsing pair, which has already cooled below the detectability possible by IUE at our derived distance of 2.3 kpc and a reddening of E(B-V)=0.40. The SuWt2 system may be destined to perish as a Type Ia supernova. (Abridged)

preprint2010arXiv

The M4 Transition: Toward a comprehensive understanding of the transition into the fully convective regime

The difference in stellar structure above and below spectral type ~M4 is expected to be a very important one, connected directly or indirectly to a variety of observational phenomena in cool stars---such as rotation, activity, magnetic field generation and topology, timescales for evolution of these, and even the basic mass-radius relationship. In this Cool Stars XVI Splinter Session, we aimed to use the M4 transition as an opportunity for discussion about the interiors of low-mass stars and the mechanisms which determine their fundamental properties. By the conclusion of the session, several key points were elucidated. Although M dwarfs exhibit significant changes across the fully convective boundary, this "M4 transition" is not observationally sharp or discrete. Instead, the properties of M dwarfs (i.e. radius, effective temperature, rotation, activity lifetime, magnetic field strength and topology) show smooth changes across M3--M6 spectral types. In addition, a wide range of stellar masses share similar spectral types around the fully convective transition. There appears to be a second transition at M6--M8 spectral types, below which there exists a clear dichotomy of magnetic field topologies. Finally, we used the information and ideas presented in the session to construct a framework for how the structure of an M dwarf star, born with specific mass and chemical composition, responds to the presence of its magnetic field, itself driven by a feedback process that links the star's rotation, interior structure, and field topology.

preprint2009arXiv

Circumstellar Environment and Effective Temperature of the Young Substellar Eclipsing Binary 2MASS J05352184-0546085

We present new Spitzer IRAC/PU/MIPS photometry from 3.6 to 24 um, and new Gemini GMOS photometry at 0.48 um, of the young brown dwarf eclipsing binary 2MASS J05352184-0546085, located in the Orion Nebula Cluster. No excess disk emission is detected: The measured fluxes at lambda < 8 um are within 1 sigma (< 0.1 mJy) of a bare photosphere, and the 3 sigma upper limit at 16 um is a mere 0.04 mJy above the bare photospheric level. Together with the known properties of the system, this implies the absence of optically thick disks around the individual components. It also implies that if any circumbinary disk is present, it must either be optically thin and extremely tenuous (10^{-10} Msun) if it extends in to within ~0.1 AU of the binary (the approximate tidal truncation radius), or it must be optically thick with a large inner hole, > 0.6-10 AU in radius depending on degree of flaring. The consequence in all cases is that disk accretion is likely to be negligible or absent. This supports the recent proposal that the strong Halpha emission in the primary (more massive) brown dwarf results from chromospheric activity, and thereby bolsters the hypothesis that the surprising Teff inversion observed between the components is due to strong magnetic fields on the primary. Our data also set constraints on the Teff of the components independent of spectral type, and thereby on models of the aforementioned magnetic field effects. We discuss the consequences for the derived fundamental properties of young brown dwarfs and very low-mass stars in general. Specifically, if very active isolated young brown dwarfs and very low-mass stars suffer the same activity/field related effects as the 2M0535-05 primary, the low-mass stellar/substellar IMF currently derived from standard evolutionary tracks may be substantially in error.

preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

preprint2008arXiv

Stellar Rotation in M35: Mass-Period Relations, Spin-Down Rates, and Gyrochronology

We present the results of a 5 month photometric time-series survey for stellar rotation over a 40'x40' field on the 150 Myr open cluster M35. We report rotation periods for 441 stars and determine their cluster membership and binarity based on a decade-long radial-velocity survey, proper-motion measurements, and multi-band photometric observations. We find that 310 of the stars with measured rotation periods are late-type members of M35. Their distribution of rotation periods span more than two orders of magnitude from ~0.1-15 days, not constrained by the sampling frequency and the time-span of the survey. With an age between the zero-age main-sequence and the Hyades, and with ~6 times more rotation periods than measured in the Pleiades, M35 permit detailed studies of early rotational evolution of late-type stars. Nearly 80% of the 310 rotators lie on two distinct sequences in the color-period plane, defining clear relations between stellar rotation period and color (mass). The M35 color-period diagram enables us to determine timescales for the transition between the two rotational states for G and K dwarfs, respectively. These timescales are inversely related to the mass of the convective envelope, and offer constraints on the rates of internal and external angular momentum transport and of the evolution of stellar dynamos. A comparison to the Hyades, confirm the Skumanich (1972) spindown-dependence for G dwarfs on one rotational state, but suggest that K dwarfs spin down more slowly. The locations of the rotational sequences in the M35 color-period diagram support the use of rotational isochrones to determine ages for coeval stellar populations. We use such gyrochronology to determine "gyro-ages" for M35. We use the M35 data to evaluate new color dependencies for the rotational isochrones.