We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

One-dimensional stellar evolution models have been successful at representing the structure and evolution of stars in diverse astrophysical contexts, but complications have been noted in the context of young, magnetically active stars, as well as close binary stars with significant tidal interactions. Numerous puzzles are associated with pre-main sequence and active main-sequence stars, relating to their radii, their colors, certain elemental abundances, and the coevality of young clusters, among others. A promising explanation for these puzzles is the distorting effects of magnetic activity and starspots on the structure of active stars. To assist the community in evaluating this hypothesis, we present the Stellar Parameters Of Tracks with Starspots (SPOTS) models, a grid of solar-metallicity stellar evolutionary tracks and isochrones which include a treatment of the structural effects of starspots. The models range from 0.1-1.3\msun and from spot-less to a surface covering fraction of 85\%, and are evolved from the pre-main sequence to the red giant branch (or 15~Gyr). We also produce two-temperature synthetic colors for our models using empirically-calibrated color tables. We describe the physical ingredients included in the SPOTS models and compare their predictions to other modern evolution codes. Finally, we apply these models to several open questions in the field of active stars, including the radii of young eclipsing binaries, the color scale of pre-main sequence stars, and the existence of sub-subgiants, demonstrating that our models can explain many peculiar features of active stars.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

The presence of certain elements within a star, and by extension its planet, strongly impacts the formation and evolution of the planetary system. The positive correlation between a host star's iron-content and the presence of an orbiting giant exoplanet has been confirmed; however, the importance of other elements in predicting giant planet occurrence is less certain despite their central role in shaping internal planetary structure. We designed and applied a machine learning algorithm to the Hypatia Catalog (Hinkel et a. 2014) to analyze the stellar abundance patterns of known host stars to determine those elements important in identifying potential giant exoplanet host stars. We analyzed a variety of different elements ensembles, namely volatiles, lithophiles, siderophiles, and Fe. We show that the relative abundances of oxygen, carbon, and sodium, in addition to iron, are influential indicators of the presence of a giant planet. We demonstrate the predictive power of our algorithm by analyzing stars with known giant planets and found that they had median 75% prediction score. We present a list of ~350 stars with no currently discovered planets that have a $\geq$90% prediction probability likelihood of hosting a giant exoplanet. We investigated archival HARPS data and found significant trends that HIP62345, HIP71803, and HIP10278 host long-period giant planet companions with estimated minimum $M_p\sin(i)$ values of 3.7, 6.8, and 8.5 M$_{J}$, respectively. We anticipate that our findings will revolutionize future target selection, the role that elements play in giant planet formation, and the determination of giant planet interior structure models.

Radius inflation continues to be explored as a peculiar occurrence among magnetically active, low-mass stars. Recently Somers et. al. 2017 showed that radius inflation among low-mass stars in the young open cluster M45 (Pleiades Cluster) is correlated to the rotation rate: faster rotators are more inflated. Here we extend that work to a sample of 68 stars of the older open Hyades Cluster. We derive the stars' spectral energy distributions to measure their bolometric fluxes. With spectroscopically defined $T_{\rm eff}$ and {\it Gaia\/} distances we calculate stellar radii using the Stefan-Boltzmann relation. We find numerous stars that exhibit significant (3-4$\sigma$) radius inflation relative to a nominal cluster isochrone. We compare these results to that of the younger Pleiades and consider radius inflation as a function of open cluster evolution. We find that unlike the Pleiades, there is not a statistically significant correlation between radius inflation and stellar rotation period. However, we do find that most inflated stars have (rapid) rotational Rossby numbers of 0.1--0.2, such that the correlation of radius inflation with Rossby number is statistically significant at 99.98\% confidence. Since the canonical rotation-activity relation of low-mass stars is understood to result from the connection between magnetic activity and surface convection, our results imply that magnetic activity within the convective layers of low-mass stars is what preferentially drives radius inflation

Internal mixing on the giant branch is an important process which affects the evolution of stars and the chemical evolution of the galaxy. While several mechanisms have been proposed to explain this mixing, better empirical constraints are necessary. Here, we use [C/N] abundances in 26097 evolved stars from the SDSS-IV/APOGEE-2 Data Release 14 to trace mixing and extra mixing in old field giants with -1.7< [Fe/H] < 0.1. We show that the APOGEE [C/N] ratios before any dredge-up occurs are metallicity dependent, but that the change in [C/N] at the first dredge-up is metallicity independent for stars above [Fe/H] ~ -1. We identify the position of the red giant branch (RGB) bump as a function of metallicity, note that a metallicity-dependent extra mixing episode takes place for low-metallicity stars ([Fe/H] <-0.4) 0.14 dex in log g above the bump, and confirm that this extra mixing is stronger at low metallicity, reaching $\Delta$ [C/N] = 0.58 dex at [Fe/H] = -1.4. We show evidence for further extra mixing on the upper giant branch, well above the bump, among the stars with [Fe/H] < -1.0. This upper giant branch mixing is stronger in the more metal-poor stars, reaching 0.38 dex in [C/N] for each 1.0 dex in log g. The APOGEE [C/N] ratios for red clump (RC) stars are significantly higher than for stars at the tip of the RGB, suggesting additional mixing processes occur during the helium flash or that unknown abundance zero points for C and N may exist among the red clump RC sample. Finally, because of extra mixing, we note that current empirical calibrations between [C/N] ratios and ages cannot be naively extrapolated for use in low-metallicity stars specifically for those above the bump in the luminosity function.

We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We present an analysis of a slightly eccentric ($e=0.05$), partially eclipsing long-period ($P = 69.73$ d) main sequence binary system (WOCS 12009, Sanders 1247) in the benchmark old open cluster M67. Using Kepler K2 and ground-based photometry along with a large set of new and reanalyzed spectra, we derived highly precise masses ($1.111\pm0.015$ and $0.748\pm0.005 M_\odot$) and radii ($1.071\pm0.008\pm0.003$ and $0.713\pm0.019\pm0.026 R_\odot$, with statistical and systematic error estimates) for the stars. The radius of the secondary star is in agreement with theory. The primary, however, is approximately $15\%$ smaller than reasonable isochrones for the cluster predict. Our best explanation is that the primary star was produced from the merger of two stars, as this can also account for the non-detection of photospheric lithium and its higher temperature relative to other cluster main sequence stars at the same $V$ magnitude. To understand the dynamical characteristics (low measured rotational line broadening of the primary star and the low eccentricity of the current binary orbit), we believe that the most probable (but not the only) explanation is the tidal evolution of a close binary within a primordial triple system (possibly after a period of Kozai-Lidov oscillations), leading to merger approximately 1Gyr ago. This star appears to be a future blue straggler that is being revealed as the cluster ages and the most massive main sequence stars die out.

We present optical continuum lags for two Seyfert 1 galaxies, MCG+08-11-011 and NGC 2617, using monitoring data from a reverberation mapping campaign carried out in 2014. Our light curves span the ugriz filters over four months, with median cadences of 1.0 and 0.6 days for MCG+08-11-011 and NGC\,2617, respectively, combined with roughly daily X-ray and near-UV data from Swift for NGC 2617. We find lags consistent with geometrically thin accretion-disk models that predict a lag-wavelength relation of $\tau \propto \lambda^{4/3}$. However, the observed lags are larger than predictions based on standard thin-disk theory by factors of 3.3 for MCG+08-11-011 and 2.3 for NGC\,2617. These differences can be explained if the mass accretion rates are larger than inferred from the optical luminosity by a factor of 4.3 in MCG+08-11-011 and a factor of 1.3 in NGC\,2617, although uncertainty in the SMBH masses determines the significance of this result. While the X-ray variability in NGC\,2617 precedes the UV/optical variability, the long 2.6 day lag is problematic for coronal reprocessing models.

We present the first results from an optical reverberation mapping campaign executed in 2014, targeting the active galactic nuclei (AGN) MCG+08-11-011, NGC 2617, NGC 4051, 3C 382, and Mrk 374. Our targets have diverse and interesting observational properties, including a "changing look" AGN and a broad-line radio galaxy. Based on continuum-H$\beta$ lags, we measure black hole masses for all five targets. We also obtain H$\gamma$ and He{\sc ii}\,$\lambda 4686$ lags for all objects except 3C 382. The He{\sc ii}\,$\lambda 4686$ lags indicate radial stratification of the BLR, and the masses derived from different emission lines are in general agreement. The relative responsivities of these lines are also in qualitative agreement with photoionization models. These spectra have extremely high signal-to-noise ratios (100--300 per pixel) and there are excellent prospects for obtaining velocity-resolved reverberation signatures.

The Kilodegree Extremely Little Telescope (KELT) has been surveying more than $70\%$ of the celestial sphere for nearly a decade. While the primary science goal of the survey is the discovery of transiting, large-radii planets around bright host stars, the survey has collected more than $10^6$ images, with a typical cadence between $10-30$ minutes, for more than $4$ million sources with apparent visual magnitudes in the approximate range $7<V<13$. Here we provide a catalog of 52,741 objects showing significant large-amplitude fluctuations likely caused by stellar variability and 62,229 objects identified with likely stellar rotation periods. The detected variability ranges in $rms$-amplitude from 3 mmag to 2.3 mag, and the detected periods range from $\sim$0.1 days to $\gtrsim$2000 days. We provide variability upper limits for all other $\sim$4 million sources. These upper limits are principally a function of stellar brightness, but we achieve typical 1$\sigma$ sensitivity on 30-minute timescales down to $\sim5$ mmag at $V\sim 8$, and down to $\sim43$ mmag at $V\sim 13$. We have matched our catalog to the $TESS$ Input catalog and the AAVSO Variable Star Index to precipitate the follow up and classification of each source. The catalog is maintained as a living database on the Filtergraph visualization portal at the URL https://filtergraph.com/kelt$\_$vars

We study the occurrence of spectroscopic binaries in young star-forming regions using the INfrared Spectroscopy of Young Nebulous Clusters(IN-SYNC) survey, carried out in SDSS-III with the APOGEE spectrograph. Multi-epoch observations of thousands of low-mass stars in Orion A, NGC 2264, NGC 1333, IC 348, and the Pleiades have been carried out, yielding H-band spectra with a nominal resolution of R=22,500 for sources with H $\le$ 12 mag. Radial velocity precisions of $\sim$0.3 $km\:s^{-1}$ were achieved, which we use to identify radial velocity variations indicative of undetected companions. We use Monte Carlo simulations to assess the types of spectroscopic binaries to which we are sensitive, finding sensitivity to binaries with orbital periods $< 10^{4}$ d, for stars with $2500 {\rm K} \le T_\mathrm{eff} \le 6000 {\rm K}$ and $\it{v} \sin \it{i}$ $\le$ 100 $km\:s^{-1}$. Using Bayesian inference, we find evidence for a decline in the spectroscopic binary fraction, by a factor of 3-4 from the age of our pre-main-sequence sample to the Pleiades age . The significance of this decline is weakened if spot-induced radial-velocity jitter is strong in the sample, and is only marginally significant when comparing any one of the pre-main-sequence clusters against the Pleiades. However, the same decline in both sense and magnitude is found for each of the five pre-main-sequence clusters, and the decline reaches statistical significance of greater than 95% confidence when considering the pre-main-sequence clusters jointly. Our results suggest that dynamical processes disrupt the widest spectroscopic binaries ($P_{\rm orb} \approx 10^3 - 10^4$ d) as clusters age, indicating that this occurs early in the stars' evolution, while they still reside within their nascent clusters.

Recent observations of the low-mass rotation distributions of the Pleiades and Praesepe clusters have revealed a ubiquitous correlation between mass and rotation, such that late M dwarfs rotate an order-of-magnitude faster than early M dwarfs. In this paper, we demonstrate that this mass-rotation correlation is present in the 10 Myr Upper Scorpius association, as revealed by new ${\it K2}$ rotation measurements. Using rotational evolution models we show that the low-mass ($0.1-0.6 M_{\odot}$) rotation distribution of the 125 Myr Pleiades cluster can only be produced if it hosted an equally strong mass-rotation correlation at 10 Myr. This suggests that physical processes important in the early pre-main sequence (star formation, accretion, disk-locking) are primarily responsible for the M dwarf rotation morphology, and not quirks of later angular momentum evolution. Such early mass trends must be taken into account when constructing initial conditions for future studies of stellar rotation. Finally, we show that the average M star loses $\sim 25-40$% of its angular momentum between 10 and 125 Myr, a figure accurately and generically predicted by modern solar-calibrated wind models. Their success rules out a lossless pre-main sequence, and validates the extrapolation of magnetic wind laws designed for solar-type stars to the low-mass regime at early times.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

One of the most well-studied young stellar associations, Taurus-Auriga, will be observed by the extended Kepler mission, K2, in the spring of 2017. K2 Campaign 13 (C13) will be a unique opportunity to study many stars in this young association at high photometric precision and cadence. Using observations from the Kilodegree Extremely Little Telescope (KELT) survey, we identify "dippers", stochastic variables, and periodic variables among K2 C13 target stars. This release of the KELT data (lightcurve data in e-tables) provides the community with long-time baseline observations to assist in the understanding of the more exotic variables in the association. Transient-like phenomena on timescales of months to years are known characteristics in the light curves of young stellar objects, making contextual pre- and post-K2 observations critical to understanding their underlying processes. We are providing a comprehensive set of the KELT light curves for known Taurus-Auriga stars in K2 C13. The combined data sets from K2 and KELT should permit a broad array of investigations related to star formation, stellar variability, and protoplanetary environments.

Almost every star in our Galaxy is likely to harbor a terrestrial planet, but accurate measurements of an exoplanet's mass and radius demands accurate knowledge of the properties of its host star. The imminent TESS and CHEOPS missions are slated to discover thousands of new exoplanets. Along with WFIRST, which will directly image nearby planets, these surveys make urgent the need to better characterize stars in the nearby solar neighborhood (< 30 pc). We have compiled the CATalog of Stellar Unified Properties (CATSUP) for 951 stars, including such data as: Gaia astrometry; multiplicity within stellar systems; stellar elemental abundance measurements; standardized spectral types; Ca II H and K stellar activity indices; GALEX NUV and FUV photometry; and X-ray fluxes and luminosities from ROSAT, XMM, and Chandra. We use this data-rich catalog to find correlations, especially between stellar emission indices, colors, and galactic velocity. Additionally, we demonstrate that thick-disk stars in the sample are generally older, have lower activity, and have higher velocities normal to the galactic plane. We anticipate CATSUP will be useful in discerning other trends among stars within the nearby solar neighborhood, for comparing thin-disk vs. thick-disk stars, for comparing stars with and without planets, and for finding correlations between chemical and kinematic properties.

High-precision photometry from space-based missions such as K2 and TESS enables detailed studies of young star variability. However, because space-based observing campaigns are often short (e.g., 80 days for K2), complementary long-baseline photometric surveys are critical for obtaining a complete understanding of young star variability, which can change on timescales of minutes to years. We therefore present and analyze light curves of members of the Upper Scorpius association made over 5.5 years by the ground-based Kilodegree Extremely Little Telescope (KELT), which complement the high-precision observations of this region taken by K2 during its Campaigns~2 and 15. We show that KELT data accurately identify the periodic signals found with high-precision K2 photometry, demonstrating the power of ground-based surveys in deriving stellar rotation periods of young stars. We also use KELT data to identify sources exhibiting variability that is likely related to circumstellar material and/or stellar activity cycles; these signatures are often unseen in the short-term K2 data, illustrating the importance of long-term monitoring surveys for studying the full range of young star variability. We provide the KELT light curves as electronic tables in an ongoing effort to establish legacy time-series datasets for young stellar clusters.

During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

The Kepler light curve of DF Cyg is unparalleled in precision and cadence for any RV Tau star to date spanning a baseline of $\sim$4 years and clearly displaying the signature pulsating behavior of alternating deep and shallow minima as well as the long-term trend indicative of an RVb-type variable. We measured DF Cyg's formal period (the time interval between two successive deep minima) to be $~49.84 \pm 0.02~$days. The trend in the arrival times emulates that of the long-term period. There appears to be precisely 16 deep+shallow minima cycles in one long-term cycle, suggesting a long-term cycle period of $\approx$ $795 \pm 5~$days. We argue that binarity may naturally explain the long-term periodicity in DF Cyg. The spectral energy distribution of DF Cyg features an infrared excess indicative of a disk possibly linked to a binary companion. From a recent Gaia parallax measurement for DF Cyg, we calculated that it has a distance of $990 \pm 372$pc and a physical radius of $R_\star = 10.3 \pm 3.8$ R$_\odot$. From kinematics and geometric arguments, we argue that the most likely interpretation for the decrease in flux from the long-period maximum to the long-period minimum, as well as the reduction of short-term pulsation amplitude, is caused by an occulting body such as a disk surrounding DF Cyg and its binary companion.

In the updated APOGEE-Kepler catalog, we have asteroseismic and spectroscopic data for over 3000 first ascent red giants. Given the size and accuracy of this sample, these data offer an unprecedented test of the accuracy of stellar models on the post-main-sequence. When we compare these data to theoretical predictions, we find a metallicity dependent temperature offset with a slope of around 100 K per dex in metallicity. We find that this effect is present in all model grids tested and that theoretical uncertainties in the models, correlated spectroscopic errors, and shifts in the asteroseismic mass scale are insufficient to explain this effect. Stellar models can be brought into agreement with the data if a metallicity dependent convective mixing length is used, with $ \Delta\alpha_{\rm ML, YREC} \sim 0.2$ per dex in metallicity, a trend inconsistent with the predictions of three dimensional stellar convection simulations. If this effect is not taken into account, isochrone ages for red giants from the Gaia data will be off by as much as a factor of 2 even at modest deviations from solar metallicity ([Fe/H]=$-$0.5).

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548.

Ages and masses of young stars are often estimated by comparing their luminosities and effective temperatures to pre-main sequence stellar evolution tracks, but magnetic fields and starspots complicate both the observations and evolution. To understand their influence, we study the heavily-spotted weak-lined T-Tauri star LkCa 4 by searching for spectral signatures of radiation originating from the starspot or starspot groups. We introduce a new methodology for constraining both the starspot filling factor and the spot temperature by fitting two-temperature stellar atmosphere models constructed from Phoenix synthetic spectra to a high-resolution near-IR IGRINS spectrum. Clearly discernable spectral features arise from both a hot photospheric component $T_{\mathrm{hot}} \sim4100$ K and to a cool component $T_{\mathrm{cool}} \sim2700-3000$ K, which covers $\sim80\%$ of the visible surface. This mix of hot and cool emission is supported by analyses of the spectral energy distribution, rotational modulation of colors and of TiO band strengths, and features in low-resolution optical/near-IR spectroscopy. Although the revised effective temperature and luminosity make LkCa 4 appear much younger and lower mass than previous estimates from unspotted stellar evolution models, appropriate estimates will require the production and adoption of spotted evolutionary models. Biases from starspots likely afflict most fully convective young stars and contribute to uncertainties in ages and age spreads of open clusters. In some spectral regions starspots act as a featureless veiling continuum owing to high rotational broadening and heavy line-blanketing in cool star spectra. Some evidence is also found for an anti-correlation between the velocities of the warm and cool components.

Precise measurements of eclipsing binary parameters and statistical studies of young clusters have suggested that some magnetically active low-mass dwarfs possess radii inflated by $\sim$5-15% relative to theoretical expectations. If true, this effect should be pronounced in young open clusters, due to the rapid rotation and strong magnetic activity of their most extreme members. We explore this possibility by determining empirical radii for 83 members of the nearby Pleiades open cluster, using spectral energy distribution fitting to establish $\mathcal{F}_{bol}$ with a typical accuracy of $\approx$3\% together with color and spectro-photometric indices to determine $T_{\rm eff}$. We find several Pleiades members with radii inflated above radius-$T_{\rm eff}$ models from state-of-the-art calculations, and apparent dispersions in radii for the K-dwarfs of the cluster. Moreover, we demonstrate that this putative radius inflation correlates strongly with rotation rate, consistent with inflation of young stars by magnetic activity and/or starspots. We argue that this signal is not a consequence of starspot-induced color anomalies, binarity, or depth effects in the cluster, employing Gaia DR1 distances as a check. Finally, we consider the lithium abundances of these stars, demonstrating a triple correlation between rotation rate, radius inflation, and enhanced lithium abundance. Our result---already significant to $\sim$99.99% confidence---provides strong support for a magnetic origin of the inflated radii and lithium dispersion observed in young, low-mass stars.

Dramatic deficiencies of Li in the mid-F dwarf stars of the Hyades cluster were discovered by Boesgaard & Tripicco. Boesgaard & King discovered corresponding, but smaller, deficiencies in Be in the same narrow temperature region in the Hyades. With the Space Telescope Imaging Spectrograph on the Hubble Space Telescope we investigate B abundances in the Hyades F stars to look for a potential B dip using the B I resonance line at 2496.8 A. The light elements, Li, Be, and B, are destroyed inside stars at increasingly hotter temperatures: 2.5, 3.5, and 5x10^6 K respectively. Consequently, these elements survive to increasingly greater depths in a star and their surface abundances indicate the depth and thoroughness of mixing in the star. We have (re)determined Li abundances/upper limits for 79 Hyades dwarfs, Be for 43 stars, and B in five stars. We find evidence for a small drop in the B abundance across the Li-Be dip. The B abundances for the four stars in the temperature range 6100-6730 K fit the B-Be correlation found previously by Boesgaard et al. Models of rotational mixing produce good agreement with the relative depletions of Be and B in the dip region. We have compared our nLTE B abundances for the three high B stars on either side of the Li-Be dip with those found by Duncan et al. for the two Hyades giants. This confirms the factor of ~10 decline in the B abundance in the Hyades giants as predicted by dilution due to the deepening of the surface convection zone.

We use high quality K2 light curves for hundreds of stars in the Pleiades to understand better the angular momentum evolution and magnetic dynamos of young, low mass stars. The K2 light curves provide not only rotational periods but also detailed information from the shape of the phased light curve not available in previous studies. A slowly rotating sequence begins at $(V-K_{\rm s})_0\sim$1.1 (spectral type F5) and ends at $(V-K_{\rm s})_0\sim$ 3.7 (spectral type K8), with periods rising from $\sim$2 to $\sim$11 days in that interval. Fifty-two percent of the Pleiades members in that color interval have periods within 30\% of a curve defining the slow sequence; the slowly rotating fraction decreases significantly redward of $(V-K_{\rm s})_0$=2.6. Nearly all of the slow-sequence stars show light curves that evolve significantly on timescales less than the K2 campaign duration. The majority of the FGK Pleiades members identified as photometric binaries are relatively rapidly rotating, perhaps because binarity inhibits star-disk angular momentum loss mechanisms during pre-main sequence evolution. The fully convective, late M dwarf Pleiades members (5.0 $<(V-K_{\rm s})_0<$ 6.0) nearly always show stable light curves, with little spot evolution or evidence of differential rotation. During PMS evolution from $\sim$3 Myr (NGC2264 age) to $\sim$125 Myr (Pleiades age), stars of 0.3 $M_{\odot}$ shed about half their angular momentum, with the fractional change in period between 3 and 125 Myr being nearly independent of mass for fully convective stars. Our data also suggest that very low mass binaries form with rotation periods more similar to each other and faster than would be true if drawn at random from the parent population of single stars.

Rotational mixing is a prime candidate for explaining the gradual depletion of lithium from the photospheres of cool stars during the main sequence. However, previous mixing calculations have relied primarily on treatments of angular momentum transport in stellar interiors incompatible with solar and stellar data, in the sense that they overestimate internal differential rotation. Instead, recent studies suggest that stars are strongly differentially rotating at young ages, but approach solid body rotation during their lifetimes. We modify our rotating stellar evolution code to include an additional source of angular momentum transport, a necessary ingredient for explaining the open cluster rotation pattern, and examine the consequences for mixing. We confirm that core-envelope recoupling with a $\sim$20 Myr timescale is required to explain the evolution of the mean solar-mass rotation pattern along the main sequence, and demonstrate that it also provides a more accurate description of the Li depletion pattern seen in open clusters. Recoupling produces a characteristic pattern of efficient mixing at early ages and little mixing at late ages, thus predicting a flattening of Li depletion at a few Gyr, in agreement with the observed late-time evolution. Using Li abundances, we argue that the timescale for core-envelope recoupling during the main sequence decreases sharply with increasing mass. We discuss implications of this finding for stellar physics, including the viability of gravity waves and magnetic fields as agents of angular momentum transport. We also raise the possibility of intrinsic differences in initial conditions in star clusters, using M67 as an example.

Young (125 Myr), populous ($>$1000 members), and relatively nearby, the Pleiades has provided an anchor for stellar angular momentum models for both younger and older stars. We used K2 to explore the distribution of rotation periods in the Pleiades. With more than 500 new periods for Pleiades members, we are vastly expanding the number of Pleiads with periods, particularly at the low mass end. About 92\% of the members in our sample have at least one measured spot-modulated rotation period. For the $\sim$8\% of the members without periods, non-astrophysical effects often dominate (saturation, etc.), such that periodic signals might have been detectable, all other things being equal. We now have an unusually complete view of the rotation distribution in the Pleiades. The relationship between $P$ and $(V-K_{\rm s})_0$ follows the overall trends found in other Pleiades studies. There is a slowly rotating sequence for $1.1\lesssim(V-K_{\rm s})_0\lesssim 3.7$, and a primarily rapidly rotating population for $(V-K_{\rm s})_0\gtrsim 5.0$. There is a region in which there seems to be a disorganized relationship between $P$ and $(V-K_{\rm s})_0$ for $3.7 \lesssim(V-K_{\rm s})_0\lesssim 5.0$. Paper II continues the discussion, focusing on multi-period structures, and Paper III speculates about the origin and evolution of the period distribution in the Pleiades.

We use K2 to continue the exploration of the distribution of rotation periods in Pleiades that we began in Paper I. We have discovered complicated multi-period behavior in Pleiades stars using these K2 data, and we have grouped them into categories, which are the focal part of this paper. About 24% of the sample has multiple, real frequencies in the periodogram, sometimes manifesting as obvious beating in the light curves. Those having complex and/or structured periodogram peaks, unresolved multiple periods, and resolved close multiple periods are likely due to spot/spot group evolution and/or latitudinal differential rotation; these largely compose the slowly rotating sequence in $P$ vs.~$(V-K_{\rm s})_0$ identified in Paper I. The fast sequence in $P$ vs.~$(V-K_{\rm s})_0$ is dominated by single-period stars; these are likely to be rotating as solid bodies. Paper III continues the discussion, speculating about the origin and evolution of the period distribution in the Pleiades.

We investigate the impact of starspots on the evolution of late-type stars during the pre-main sequence (pre-MS). We find that heavy spot coverage increases the radii of stars by 4-10%, consistent with inflation factors in eclipsing binary systems, and suppresses the rate of pre-MS lithium depletion, leading to a dispersion in zero-age MS Li abundance (comparable to observed spreads) if a range of spot properties exist within clusters from 3-10 Myr. This concordance with data implies that spots induce a range of radii at fixed mass during the pre-MS. These spots decrease the luminosity and $T_{\rm eff}$ of stars, leading to a displacement on the HR diagram. This displacement causes isochrone derived masses and ages to be systematically under-estimated, and can lead to the spurious appearance of an age spread in a co-eval population.

We assess the impact of starspots on the evolution of late-type stars during the pre-main sequence (pre-MS) using a modified stellar evolution code. We find that heavily spotted models of mass 0.1-1.2\msun\ are inflated by up to $10$% during the pre-MS, and up to 4% and 9% for fully- and partially-convective stars at the zero-age MS, consistent with measurements from active eclipsing binary systems. Spots similarly decrease stellar luminosity and $T_{\rm eff}$, causing isochrone-derived masses to be under-estimated by up to a factor of $2 \times$, and ages to be under-estimated by a factor of 2-10$\times$, at 3 Myr. Consequently, pre-MS clusters and their active stars are systematically older and more massive than often reported. Cluster ages derived with the lithium depletion boundary technique are erroneously young by $\sim 15$% and $10$% at $30$ and $100$ Myr respectively, if 50% spotted stars are interpreted with un-spotted models. Finally, lithium depletion is suppressed in spotted stars with radiative cores, leading to a fixed-temperature lithium dispersion on the MS if a range of spot properties are present on the pre-MS. Such dispersions are large enough to explain Li abundance spreads seen in young open clusters, and imply a range of radii at fixed mass and age during the pre-MS Li burning epoch. By extension, this implies that mass, composition, and age do not uniquely specify the HR diagram location of pre-MS stars.

The rapidly rotating cool dwarfs of the Pleiades are rich in lithium relative to their slowly rotating counterparts. Motivated by observations of inflated radii in young, active stars, and by calculations showing that radius inflation inhibits pre-main sequence (pre-MS) Li destruction, we test whether this pattern could arise from a connection between stellar rotation rate and radius inflation on the pre-MS. We demonstrate that pre-MS radius inflation can efficiently suppress lithium destruction by rotationally induced mixing in evolutionary models, and that the net effect of inflation and rotational mixing is a pattern where rotation correlates with lithium abundance for $M_{*} < {\rm M}_{\odot}$, and anti-correlates with lithium abundance for $M_{*} > {\rm M}_{\odot}$, similar to the empirical trend in the Pleiades. Next, we adopt different prescriptions for the dependence of inflation on rotation, and compare their predictions to the Pleiades lithium/rotation pattern. We find that if a connection between rotation and radius inflation exists, then the important qualitative features of this pattern naturally and generically emerge in our models. This is the first consistent physical model to date that explains the Li--rotation correlation in the Pleiades. We discuss plausible mechanisms for inducing this correlation and suggest an observational test using granulation.

We investigate lithium depletion in standard stellar models (SSMs) and main sequence (MS) open clusters, and explore the origin of the Li dispersion in young, cool stars of equal mass, age and composition. We first demonstrate that SSMs accurately predict the Li abundances of solar analogs at the zero-age main sequence (ZAMS) within theoretical uncertainties. We then measure the rate of MS Li depletion by removing the [Fe/H]-dependent ZAMS Li pattern from three well-studied clusters, and comparing the detrended data. MS depletion is found to be mass dependent, in the sense of more depletion at low mass. A dispersion in Li abundance at fixed $T_{\rm eff}$ is nearly universal, and sets in by $\sim$200 Myr. We discuss mass and age dispersion trends, and the pattern is mixed. We argue that metallicity impacts the ZAMS Li pattern, in agreement with theoretical expectations but contrary to the findings of some previous studies, and suggest Li as a test of cluster metallicity. Finally, we argue that a radius dispersion in stars of fixed mass and age, during the epoch of pre-MS Li destruction, is responsible for the spread in Li abundances and the correlation between rotation and Li in young cool stars, most well known in the Pleiades. We calculate stellar models, inflated to match observed radius anomalies in magnetically active systems, and the resulting range of Li abundances reproduces the observed patterns of young clusters. We discuss ramifications for pre-MS evolutionary tracks and age measurements of young clusters, and suggest an observational test.

We have used Chandra observations of eight bulge-less disk galaxies to identify new ultraluminous X-ray source (ULX) candidates, study their high mass X-ray binary (HMXB) population, and search for low-luminosity active galactic nuclei (AGN). We report the discovery of 16 new ULX candidates in our sample of galaxies. Eight of these are found in the star forming galaxy NGC 337, none of which are expected to be background contaminants. The HMXB luminosity function of NGC 337 implies a star formation rate (SFR) of 6.8$^{+4.4}_{-3.5}$ \msun\ yr$^{-1}$, consistent at 1.5$\sigma$ with a recent state of the art SFR determination. We also report the discovery of a bright ULX candidate (X-1) in ESO 501-23. X-1's spectrum is well fit by an absorbed power law with $\Gamma = 1.18^{+0.19}_{-0.11}$ and N$\rm{_H}$ = 1.13$^{+7.07}_{-1.13} \times 10^{20}$ cm$^{-2}$, implying a 0.3-8 keV flux of $1.08^{+0.05}_{-0.07} \times 10^{-12}$ \esc. Its X-ray luminosity (L$_X$) is poorly constrained due to uncertainties in the host galaxy's distance, but we argue that its spectrum implies L$_X > 10^{40}$ \es. An optical counterpart to this object may be present in HST imaging. We also identify ULX candidates in IC 1291, PGC 3853, NGC 5964 and NGC 2805. We find no evidence of nuclear activity in the galaxies in our sample, placing a flux upper limit of $4 \times 10^{-15}$ \esc\ on putative AGN. Additionally, the type II-P supernova SN 2011DQ in NGC 337, which exploded 2 months before our X-ray observation, is undetected.