Search

The vast majority of Milky Way stellar halo stars were likely accreted from a small number ($\lesssim$3) of relatively large dwarf galaxy accretion events. However, the timing of these events is poorly constrained, relying predominantly on indirect dynamical mixing arguments or imprecise age measurements of stars associated with debris structures. Here, we aim to infer robust stellar ages for stars associated with galactic substructures to more directly constrain the merger history of the Galaxy. By combining kinematic, asteroseismic, and spectroscopic data where available, we infer stellar ages for a sample of 10 red giant stars that were kinematically selected to be associated with the stellar halo, a subset of which are associated with the Gaia-Enceladus-Sausage halo substructure, and compare their ages to 3 red giant stars in the Galactic disk. Despite systematic differences in both absolute and relative ages determined by this work, age rankings of stars in this sample are robust. Passing the same observable inputs to multiple stellar age determination packages, we measure a weighted average age for the Gaia-Enceladus-Sausage stars in our sample of 8 $\pm$ 3 (stat.) $\pm$ 1 (sys.) Gyr. We also determine hierarchical ages for the populations of Gaia-Enceladus-Sausage, in situ halo and disk stars, finding a Gaia-Enceladus-Sausage population age of 8.0$^{+3.2}_{-2.3}$ Gyr. Although we cannot distinguish hierarchical population ages of halo or disk structures with our limited data and sample of stars, this framework should allow distinct characterization of Galactic substructures using larger stellar samples and additional data available in the near future., 28 pages, 21 figures. Published by AAS Journals

The light curve of KIC 8462852, aka Boyajian’s Star, undergoes deep dips the origin of which remains unclear. A faint star ≈2″ to the east was discovered in Keck/NIRC2 imaging in Boyajian et al., but its status as a binary, and possible contribution to the observed variability, was unclear. Here, we use three epochs of Keck/NIRC2 imaging, spanning 5 yr, in JHK near-infrared bands to obtain 1 mas precision astrometry. We show that the two objects exhibit common proper motion, measure a relative velocity of μ = 0.14 ± 0.44 mas yr−1 (μ = 0.30 ± 0.93 km s−1) and conclude that they are a binary pair at 880 ± 10 au projected separation. There is marginal detection of possible orbital motion, but our astrometry is insufficient to characterize the orbit. We show that two other point sources are not associated with KIC 8462852. We recommend that attempts to model KIC 8462852 A’s light curve should revisit the possibility that the bound stellar companion may play a role in causing the irregular brightness variations, for example, through disruption of the orbits of bodies around the primary due to long-term orbital evolution of the binary orbit.

High precision spectrographs can enable not only the discovery of exoplanets, but can also provide a fundamental measurement in Galactic dynamics. Over about ten year baselines, the expected change in the line-of-sight velocity due to the Galaxy's gravitational field for stars at $\sim$ kpc scale distances above the Galactic mid-plane is $\sim$ few - 10 cm/s, and may be detectable by the current generation of high precision spectrographs. Here, we provide theoretical expectations for this measurement based on both static models of the Milky Way and isolated Milky Way simulations, as well from controlled dynamical simulations of the Milky Way interacting with dwarf galaxies. We simulate a population synthesis model to analyze the contribution of planets and binaries to the Galactic acceleration signal. We find that while low-mass, long-period planetary companions are a contaminant to the Galactic acceleration signal, their contribution is very small. Our analysis of $\sim$ ten years of data from the LCES HIRES/Keck precision radial velocity (RV) survey shows that slopes of the RV curves of standard RV stars agree with expectations of the local Galactic acceleration near the Sun within the errors, and that the error in the slope scales inversely as the square root of the number of observations. Thus, we demonstrate that a survey of stars with low intrinsic stellar jitter at kpc distances above the Galactic mid-plane for realistic sample sizes can enable a direct determination of the dark matter density., accepted to ApJ Letters

Combining visual and spectroscopic orbits of binary stars leads to a determination of the full 3D orbit, individual masses, and distance to the system. We present a full analysis of the evolved binary system $\delta$ Delphini using astrometric data from the MIRC and PAVO instruments on the CHARA long-baseline interferometer, 97 new spectra from the Fairborn Observatory, and 87 unpublished spectra from Lick Observatory. We determine the full set of orbital elements for $\delta$ Del, along with masses of $1.78 \pm 0.07$ $M_{\odot}$ and $1.62 \pm 0.07$ $M_{\odot}$ for each component, and a distance of $63.61 \pm 0.89$ pc. These results are important in two contexts: for testing stellar evolution models and defining the detection capabilities for future planet searches. We find that the evolutionary state of this system is puzzling, as our measured flux ratios, radii, and masses imply a $\sim$ 200 Myr age difference between the components using standard stellar evolution models. Possible explanations for this age discrepancy include mass transfer scenarios with a now ejected tertiary companion. For individual measurements taken over a span of 2 years we achieve $2$ M$_{J}$ on orbits $>0.75$ AU around individual components of hot binary stars via differential astrometry., Comment: 28 pages, 14 figures, accepted for publication in ApJ

We present new photometry of Procyon, obtained by MOST during a 38 day run in 2007, and frequency analyses of those data. The long time coverage and low point-to-point scatter of the light curve yield an average noise amplitude of about 1.5-2.0 ppm in the frequency range 500-1500 μHz. This is half the noise level obtained from each of the previous two Procyon campaigns by MOST in 2004 and 2005. The 2007 MOST amplitude spectrum shows some evidence for p-mode signal: excess power centered near 1000 μHz and an autocorrelation signal near 55 μHz (suggestive of a mode spacing around that frequency), both consistent with p-mode model predictions. However, we do not see regularly spaced frequencies aligned in common l-valued ridges in echelle diagrams of the most significant peaks in the spectrum unless we select modes from the spectrum using a priori assumptions. The most significant peaks in the spectrum are scattered by more than ±5 μHz about the predicted l-valued ridges, a value that is consistent with the scatter among individually identified frequencies obtained from ground-based radial velocity (RV) observations. We argue that the observed scatter is intrinsic to the star, due to short lifetimes of the modes and the dynamic structure of Procyon's thin convection zone. We compare the MOST Procyon amplitude and power density spectra with preliminary results of three-dimensional numerical models of convection by the Yale group. These models show that, unlike in the Sun, Procyon's granulation signal in luminosity has a peak coinciding with the expected frequency region for p-modes near 1000 μHz.

The Microvariability and Oscillations of stars (MOST) satellite has detected low-amplitude light variations ($\Delta m\sim$1 mmag) in the Be star $\beta$ CMi (B8Ve). The observations lasted 41 days and the variations have typical periods $\sim 0.3$ days. We demonstrate that the dominant frequencies are consistent with prograde high-order g-modes of $m=-1$ excited by the Fe-bump of opacity in an intermediate-mass ($\approx 3.5 M_\odot$) star with a nearly critical rotation period of 0.38 days. This is the first detection of nonradial g-mode pulsations in a Be star later than B6 leading to the possibility that pulsations are excited in all classical Be stars.

We used the Robo-AO laser adaptive optics system to image 99 main sequence and subgiant stars that have Kepler-detected asteroseismic signals. Robo-AO allows us to resolve blended secondary sources at separations as close as 0.15" that may contribute to the measured Kepler light curves and affect asteroseismic analysis and interpretation. We report 8 new secondary sources within 4.0" of these Kepler asteroseismic stars. We used Subaru and Keck adaptive optics to measure differential infrared photometry for these candidate companion systems. Two of the secondary sources are likely foreground objects and at least 6 of the secondaries are background sources; however we cannot exclude the possibility that three of the objects may be physically associated. We measured a range of i'-band amplitude dilutions for the candidate companion systems from 0.43% to 15.4%. We find that the measured amplitude dilutions are insufficient to explain the previously identified excess scatter in the relationship between asteroseismic oscillation amplitude and the frequency of maximum power., 16 pages, 8 figures, 5 tables. Updated version, accepted to AJ

$\theta$ Cygni is an F3 spectral-type main-sequence star with visual magnitude V=4.48. This star was the brightest star observed by the original Kepler spacecraft mission. Short-cadence (58.8 s) photometric data using a custom aperture were obtained during Quarter 6 (June-September 2010) and subsequently in Quarters 8 and 12-17. We present analyses of the solar-like oscillations based on Q6 and Q8 data, identifying angular degree $l$ = 0, 1, and 2 oscillations in the range 1000-2700 microHz, with a large frequency separation of 83.9 plus/minus 0.4 microHz, and frequency with maximum amplitude 1829 plus/minus 54 microHz. We also present analyses of new ground-based spectroscopic observations, which, when combined with angular diameter measurements from interferometry and Hipparcos parallax, give T_eff = 6697 plus/minus 78 K, radius 1.49 plus/minus 0.03 solar radii, [Fe/H] = -0.02 plus/minus 0.06 dex, and log g = 4.23 plus/minus 0.03. We calculate stellar models matching the constraints using several methods, including using the Yale Rotating Evolution Code and the Asteroseismic Modeling Portal. The best-fit models have masses 1.35-1.39 solar masses and ages 1.0-1.6 Gyr. theta Cyg's T_eff and log g place it cooler than the red edge of the gamma Doradus instability region established from pre-Kepler ground-based observations, but just at the red edge derived from pulsation modeling. The pulsation models show gamma Dor gravity-mode pulsations driven by the convective-blocking mechanism, with frequencies of 1 to 3 cycles/day (11 to 33 microHz). However, gravity modes were not detected in the Kepler data, one signal at 1.776 cycles/day (20.56 microHz) may be attributable to a faint, possibly background, binary. Asteroseismic studies of theta Cyg and other A-F stars observed by Kepler and CoRoT, will help to improve stellar model physics and to test pulsation driving mechanisms., Comment: Accepted for Publication in The Astrophysical Journal, July 1, 2016

Asteroseismology has proven to be an excellent tool to determine not only the global stellar properties with a good precision but also to infer stellar structure, dynamics, and evolution for a large sample of Kepler stars. Prior to the launch of the mission the properties of Kepler targets were inferred from broadband photometry, leading to the Input Catalog (KIC Brown et al. 2011). The KIC was later revised in the Kepler Star Properties Catalog (Huber et al. 2014), based on literature values and an asteroseismic analysis of stars which were unclassified in the KIC. Here we present an asteroseismic analysis of 45,400 stars which were classified as dwarfs in the Kepler Star Properties Catalog. We found that around 2% of the sample shows acoustic modes in the typical frequency range that put them in the red-giant category rather than cool dwarfs. We analyse the asteroseismic properties of these stars, derive their surface gravities, masses, and radii and present updated effective temperatures and distances. We show that the sample is significantly fainter than the previously known oscillating giants in the Kepler field, with the faintest stars reaching down to a Kepler magnitude, Kp~16. We demonstrate that 404 stars are at distances beyond 5 kpc and that the stars are significantly less massive than for the original Kepler red-giant sample, consistent with a population of distant halo giants. A comparison with a galactic population model shows that up to 40 stars might be genuine halo giants, which would increase the number of known asteroseismic halo stars by a factor of 4. The detections presented here will provide a valuable sample for galactic archeology studies., Comment: 14 pages, 13 figures. Accepted for publication in ApJ

We present the first results from our Keck program investigating the orbital architectures of planet-hosting multiple star systems. Kepler-444 is a metal-poor triple star system that hosts five sub-Earth-sized planets orbiting the primary star (Kepler-444A), as well as a spatially unresolved pair of M dwarfs (Kepler-444BC) at a projected distance of 1.8" (66 AU). We combine our Keck/NIRC2 adaptive optics astrometry with multi-epoch Keck/HIRES RVs of all three stars to determine a precise orbit for the BC pair around A, given their empirically constrained masses. We measure minimal astrometric motion ($1.0\pm0.6$ mas yr$^{-1}$, or $0.17\pm0.10$ km s$^{-1}$), but our RVs reveal significant orbital velocity ($1.7\pm0.2$ km s$^{-1}$) and acceleration ($7.8\pm0.5$ m s$^{-1}$ yr$^{-1}$). We determine a highly eccentric stellar orbit ($e=0.864\pm0.023$) that brings the tight M dwarf pair within $5.0^{+0.9}_{-1.0}$ AU of the planetary system. We validate that the system is dynamically stable in its present configuration via n-body simulations. We find that the A$-$BC orbit and planetary orbits are likely aligned (98%) given that they both have edge-on orbits and misalignment induces precession of the planets out of transit. We conclude that the stars were likely on their current orbits during the epoch of planet formation, truncating the protoplanetary disk at $\approx$2 AU. This truncated disk would have been severely depleted of solid material from which to form the total $\approx$1.5 $M_{\rm Earth}$ of planets. We thereby strongly constrain the efficiency of the conversion of dust into planets and suggest that the Kepler-444 system is consistent with models that explain the formation of more typical close-in Kepler planets in normal, not truncated, disks., accepted to ApJ

We measure planet occurrence rates using the planet candidates discovered by the Q1-Q16 Kepler pipeline search. This study examines planet occurrence rates for the Kepler GK dwarf target sample for planet radii, 0.75, Comment: 19 Pages, 17 Figures, Submitted ApJ. Python source to support Kepler pipeline completeness estimates available at http://github.com/christopherburke/KeplerPORTs/

Kepler-91b is a rare example of a transiting hot Jupiter around a red giant star, providing the possibility to study the formation and composition of hot Jupiters under different conditions compared to main-sequence stars. However, the planetary nature of Kepler-91b, which was confirmed using phase-curve variations by Lillo-Box et al., was recently called into question based on a re-analysis of Kepler data. We have obtained ground-based radial velocity observations from the Hobby-Eberly Telescope and unambiguously confirm the planetary nature of Kepler-91b by simultaneously modeling the Kepler and radial velocity data. The star exhibits temporally correlated noise due to stellar granulation which we model as a Gaussian Process. We hypothesize that it is this noise component that led previous studies to suspect Kepler-91b to be a false positive. Our work confirms the conclusions presented by Lillo-Box et al. that Kepler-91b is a 0.73+/-0.13 Mjup planet orbiting a red giant star., Published in ApJ

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

We have searched the Kepler light curves of ~3900 M-star targets for evidence of periodicities that indicate, by means of the effects of starspots, rapid stellar rotation. Several analysis techniques, including Fourier transforms, inspection of folded light curves, 'sonograms', and phase tracking of individual modulation cycles, were applied in order to distinguish the periodicities due to rapid rotation from those due to stellar pulsations, eclipsing binaries, or transiting planets. We find 178 Kepler M-star targets with rotation periods, P_rot, of < 2 days, and 110 with P_rot < 1 day. Some 30 of the 178 systems exhibit two or more independent short periods within the same Kepler photometric aperture, while several have three or more short periods. Adaptive optics imaging and modeling of the Kepler pixel response function for a subset of our sample support the conclusion that the targets with multiple periods are highly likely to be relatively young physical binary, triple, and even quadruple M star systems. We explore in detail the one object with four incommensurate periods all less than 1.2 days, and show that two of the periods arise from one of a close pair of stars, while the other two arise from the second star, which itself is probably a visual binary. If most of these M-star systems with multiple periods turn out to be bound M stars, this could prove a valuable way of discovering young hierarchical M-star systems; the same approach may also be applicable to G and K stars. The ~5% occurrence rate of rapid rotation among the ~3900 M star targets is consistent with spin evolution models that include an initial contraction phase followed by magnetic braking, wherein a typical M star can spend several hundred Myr before spinning down to periods longer than 2 days., 17 pages, 12 figures, 2 tables; accepted for publication in The Astrophysical Journal

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}(\rho_{\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$, Comment: 6 pages, 3 figures, 1 table. Accepted to ApJ Letters. Code available at https://www.cfa.harvard.edu/~dkipping/flicker.html

We present a novel method for estimating lower-limit surface gravities log g of Kepler targets whose data do not allow the detection of solar-like oscillations. The method is tested using an ensemble of solar-type stars observed in the context of the Kepler Asteroseismic Science Consortium. We then proceed to estimate lower-limit log g for a cohort of Kepler solar-type planet-candidate host stars with no detected oscillations. Limits on fundamental stellar properties, as provided by this work, are likely to be useful in the characterization of the corresponding candidate planetary systems. Furthermore, an important byproduct of the current work is the confirmation that amplitudes of solar-like oscillations are suppressed in stars with increased levels of surface magnetic activity., Accepted for publication in ApJ; 35 pages, 10 figures, 5 tables

With the success of the Kepler and CoRoT missions, the number of stars with detected solar-like oscillations has increased by several orders of magnitude, for the first time we are able to perform large-scale ensemble asteroseismology of these stars. In preparation for this golden age of asteroseismology we have computed expected values of various asteroseismic observables from models of varying mass and metallicity. The relationships between these asteroseismic observables, such as the separations between mode frequencies, are able to significantly constrain estimates of the ages and masses of these stars. We investigate the scaling relation between the large frequency separation, Delta nu, and mean stellar density. Furthermore we present model evolutionary tracks for several asteroseismic diagrams. We have extended the so-called C-D diagram beyond the main sequence to the subgiants and the red-giant branch. We also consider another asteroseismic diagram, the epsilon diagram, which is more sensitive to variations in stellar properties at the subgiant stages and can aid in determining the correct mode identification. The recent discovery of gravity-mode period spacings in red giants forms the basis for a third asteroseismic diagram. We compare the evolutionary model tracks in these asteroseismic diagrams with results from pre-Kepler studies of solar-like oscillations, and early results from Kepler., Comment: 14 pages, 15 figures and 2 tables, accepted for publication in the Astrophysical Journal

We have analyzed solar-like oscillations in ~1700 stars observed by the Kepler Mission, spanning from the main-sequence to the red clump. Using evolutionary models, we test asteroseismic scaling relations for the frequency of maximum power (nu_max), the large frequency separation (Delta_nu) and oscillation amplitudes. We show that the difference of the Delta_nu-nu_max relation for unevolved and evolved stars can be explained by different distributions in effective temperature and stellar mass, in agreement with what is expected from scaling relations. For oscillation amplitudes, we show that neither (L/M)^s scaling nor the revised scaling relation by Kjeldsen & Bedding (2011) are accurate for red-giant stars, and demonstrate that a revised scaling relation with a separate luminosity-mass dependence can be used to calculate amplitudes from the main-sequence to red-giants to a precision of ~25%. The residuals show an offset particularly for unevolved stars, suggesting that an additional physical dependency is necessary to fully reproduce the observed amplitudes. We investigate correlations between amplitudes and stellar activity, and find evidence that the effect of amplitude suppression is most pronounced for subgiant stars. Finally, we test the location of the cool edge of the instability strip in the Hertzsprung-Russell diagram using solar-like oscillations and find the detections in the hottest stars compatible with a domain of hybrid stochastically excited and opacity driven pulsation., 11 pages, 9 figures, accepted for publication in ApJ

We have analysed oscillations of the red giant star HD 186355 observed by the NASA Kepler satellite. The data consist of the first five quarters of science operations of Kepler, which cover about 13 months. The high-precision time-series data allow us to accurately extract the oscillation frequencies from the power spectrum. We find the frequency of the maximum oscillation power, {\nu}_max, and the mean large frequency separation, {\Delta}{\nu}, are around 106 and 9.4 {\mu}Hz respectively. A regular pattern of radial and non-radial oscillation modes is identified by stacking the power spectra in an echelle diagram. We use the scaling relations of {\Delta}{\nu} and {\nu}_max to estimate the preliminary asteroseismic mass, which is confirmed with the modelling result (M = 1.45 \pm 0.05 M_sun) using the Yale Rotating stellar Evolution Code (YREC7). In addition, we constrain the effective temperature, luminosity and radius from comparisons between observational constraints and models. A number of mixed l = 1 modes are also detected and taken into account in our model comparisons. We find a mean observational period spacing for these mixed modes of about 58 s, suggesting that this red giant branch star is in the shell hydrogen-burning phase., Comment: 26 pages, 5 figures and 2 tables

We compare the simultaneous 2007 space-based MOST photometry and ground-based radial velocity observations of the F5 star Procyon. We identify slow variations in the MOST data that are similar to those reported in the radial velocity (RV) time series, and confirm by comparison with the Sun that these variations are likely the signature of stellar activity. The MOST power spectrum yields clear evidence for individual oscillation frequencies that match those found in the radial velocity data by Bedding et al. (2010). We identify the same ridges due to modes of different spherical degree in both datasets, but are not able to confirm a definite ridge identification using the MOST data. We measure the luminosity amplitude per radial mode A_{l=0, phot} = 9.1 +/- 0.5 ppm. Combined with the estimate for the RV data by Arentoft et al. (2008) this gives a mean amplitude ratio of A_{l=0, phot}/A_{l=0, RV} = 0.24 +/- 0.02 ppm cm^{-1} s, considerably higher than expected from scaling relations but in reasonable agreement with theoretical models by Houdek (2010). We also compare the amplitude ratio as a function of frequency, and find that the maximum of the oscillation envelope is shifted to higher frequencies in photometry than in velocity., Comment: 10 pages, 9 figures, accepted for publication in ApJ; v2 includes minor changes made in the proofs (updated references & corrected typos)