We summarize some of the compelling new scientific opportunities for understanding stars and stellar systems that can be enabled by sub-milliarcsec (sub-mas) angular resolution, UV-Optical spectral imaging observations, which can reveal the details of the many dynamic processes (e.g., evolving magnetic fields, accretion, convection, shocks, pulsations, winds, and jets) that affect stellar formation, structure, and evolution. These observations can only be provided by long-baseline interferometers or sparse aperture telescopes in space, since the aperture diameters required are in excess of 500 m (a regime in which monolithic or segmented designs are not and will not be feasible) and since they require observations at wavelengths (UV) not accessible from the ground. Such observational capabilities would enable tremendous gains in our understanding of the individual stars and stellar systems that are the building blocks of our Universe and which serve as the hosts for life throughout the Cosmos.

In response to the need for the Astro2020 Decadal Survey to explicitly engage early career astronomers, the National Academies of Sciences, Engineering, and Medicine hosted the Early Career Astronomer and Astrophysicist Focus Session (ECFS) on October 8-9, 2018 under the auspices of Committee of Astronomy and Astrophysics. The meeting was attended by fifty six pre-tenure faculty, research scientists, postdoctoral scholars, and senior graduate students, as well as eight former decadal survey committee members, who acted as facilitators. The event was designed to educate early career astronomers about the decadal survey process, to solicit their feedback on the role that early career astronomers should play in Astro2020, and to provide a forum for the discussion of a wide range of topics regarding the astrophysics career path. This white paper presents highlights and themes that emerged during two days of discussion. In Section 1, we discuss concerns that emerged regarding the coming decade and the astrophysics career path, as well as specific recommendations from participants regarding how to address them. We have organized these concerns and suggestions into five broad themes. These include (sequentially): (1) adequately training astronomers in the statistical and computational techniques necessary in an era of "big data", (2) responses to the growth of collaborations and telescopes, (3) concerns about the adequacy of graduate and postdoctoral training, (4) the need for improvements in equity and inclusion in astronomy, and (5) smoothing and facilitating transitions between early career stages. Section 2 is focused on ideas regarding the decadal survey itself, including: incorporating early career voices, ensuring diverse input from a variety of stakeholders, and successfully and broadly disseminating the results of the survey.

Asteroseismology is the only observational tool in astronomy that can probe the interiors of stars, and is a benchmark method for deriving fundamental properties of stars and exoplanets. Over the coming decade, space-based and ground-based observations will provide a several order of magnitude increase of solar-like oscillators, as well as a dramatic increase in the number and quality of classical pulsator observations, providing unprecedented possibilities to study stellar physics and galactic stellar populations. In this white paper, we describe key science questions and necessary facilities to continue the asteroseismology revolution into the 2020's.

Direct determination of fundamental stellar parameters has many profound and wide-ranging impacts throughout astrophysics. These determinations are rooted in high angular resolution observations. In particular, as long-baseline optical interferometry has matured over the past decade, increasingly large survey samples are serving to empirically ground the basic parameters of these building blocks of the universe. True imaging and improved parametric fitting are becoming routinely available, an essential component of fully characterizing stars, stellar environments, and planets these stars may host.

Evolved stars dominate galactic spectra, enrich the galactic medium, expand to change their planetary systems, eject winds of a complex nature, produce spectacular nebulae and illuminate them, and transfer material between binary companions. While doing this, they fill the HR diagram with diagnostic loops that write the story of late stellar evolution. Evolved stars sometimes release unfathomable amounts of energy in neutrinos, light, kinetic flow, and gravitational waves. During these late-life times, stars evolve complexly, with expansion, convection, mixing, pulsation, mass loss. Some processes have virtually no spatial symmetries, and are poorly addressed with low-resolution measurements and analysis. Even a "simple" question as how to model mass loss resists solution. However, new methods offer increasingly diagnostic tools. Astrometry reveals populations and groupings. Pulsations/oscillations support study of stellar interiors. Optical/radio interferometry enable 2-3d imagery of atmospheres and shells. Bright stars with rich molecular spectra and velocity fields are a ripe opportunity for imaging with high spatial and spectral resolution, giving insight into the physics and modeling of later stellar evolution.

Cool, evolved stars are the main source of chemical enrichment of the interstellar medium (ISM), and understanding their mass loss and structure offers a unique opportunity to study the cycle of matter in the Universe. Pulsation, convection, and other dynamic processes in cool evolved stars create an atmosphere where molecules and dust can form, including those necessary to the formation of life (e.g.~Carbon-bearing molecules). Understanding the structure and composition of these stars is thus vital to several aspects of stellar astrophysics, ranging from ISM studies to modeling young galaxies and exoplanet research. Recent modeling efforts and increasingly precise observations now reveal that our understanding of cool stars photospheric, chromospheric, and atmospheric structures is limited by inadequate knowledge of the dynamic and chemical processes at work. Here we outline promising scientific opportunities for the next decade. We identify and discuss the following main opportunities: (1) identify and model the physical processes that must be included in current 1D and 3D atmosphere models of cool, evolved stars; (2) refine our understanding of photospheric, chromospheric, and outer atmospheric regions of cool evolved stars, their properties and parameters, through high-resolution spectroscopic observations, and interferometric observations at high angular resolution; (3) include the neglected role of chromospheric activity in the mass loss process of red giant branch and red super giant stars, and understand the role played by their magnetic fields; (4) identify the important shaping mechanisms for planetary nebulae and their relation with the parent asymptotic giant branch stars.

We are now in an era where we can image details on the surfaces of stars. When resolving stellar surfaces, we see that every surface is uniquely complicated. Each imaged star provides insight into not only the stellar surface structures, but also the stellar interiors suggesting constraints on evolution and dynamo models. As more resources become operational in the coming years, imaging stellar surfaces should become commonplace for revealing the true nature of stars. Here, we discuss the main types of stars for which imaging surface features is currently useful and what improved observing techniques would provide for imaging stellar surface features.

We aim at detecting and characterizing the main-sequence companions of a sample of known and suspected Galactic binary Cepheids. We used the multi-telescope interferometric combiners MIRC and PIONIER to detect and measure the astrometric positions of the high-contrast companions orbiting 16 bright Galactic Cepheids. We made use of the CANDID algorithm to search for the companions and set detection limits. We also present new high-precision RVs which were used to fit radial pulsation and orbital velocities. We report the detection of the companions orbiting the Cepheids U Aql, BP Cir, and S Mus, possible detections for FF Aql, Y Car, BG Cru, X Sgr, V350 Sgr, and V636 Sco, while no component is detected around U Car, YZ Car, T Mon, R Mus, S Nor, W Sgr and AH Vel. For U Aql and S Mus, we performed a preliminary orbital fit combining astrometry with new high-precision RVs, providing the full set of orbital and pulsation parameters. Assuming the distance, we estimated preliminary masses of M(U Aql) = 4.97+/-0.62Msun and M(S Mus) = 4.63+/-0.99Msun. For YZ Car, W Sgr, V350 Sgr, and V636 Sco, we revised the spectroscopic orbits using new high-precision RVs, while we updated the pulsation parameters for BP Cir, BG Cru, S Nor and AH Vel. Our interferometric observations also provide measurements of the angular diameters, that can be used in a Baade-Wesselink type analysis. We have now several astrometric detections of Cepheid companions. When RVs of the companions will be available, such systems will provide accurate and independent masses and distances. Orbital parallaxes with an accuracy <5% will be particularly useful for a better calibration of the P-L relation. The final Gaia parallaxes will be also particularly helpful for single-line spectroscopic systems, where mass and distance are degenerate. Mass measurements are necessary for a better understanding of the age and evolution of Cepheids

Starspots and flares are indicators of stellar magnetic activity that can both be studied in greater detail by utilizing the long-term, space-based archive of the Kepler satellite. Here, we aim to investigate a subset of the Kepler archive to reveal a connection between the starspots and the stellar flares, in order to provide insight into the overall stellar magnetic field. We use the flare-finding algorithm FLATW'RM in conjunction with a new suite of algorithms that aim to locate the local minima caused by starspot groups. We compare the phase difference between the time of maximum flux of a flare and the time of minimum stellar flux due to a starspot group. The strongest flares do not appear to be correlated to the largest starspot group present, but are also not uniformly distributed in phase with respect to the starspot group. The weaker flares, however, do show an increased occurrence close to the starspot groups.

Cepheid stars play a considerable role as extragalactic distances indicators, thanks to the simple empirical relation between their pulsation period and their luminosity. They overlap with that of secondary distance indicators, such as Type Ia supernovae, whose distance scale is tied to Cepheid luminosities. However, the Period-Luminosity (P-L) relation still lacks a calibration to better than 5%. Using an original combination of interferometric astrometry with optical and ultraviolet spectroscopy, we measured the geometrical distance d = 720.35+/-7.84 pc of the 3.33 d period Cepheid V1334 Cyg with an unprecedented accuracy of +/-1 %, providing the most accurate distance for a Cepheid. Placing this star in the P-L diagram provides an independent test of existing period-luminosity relations. We show that the secondary star has a significant impact on the integrated magnitude, particularly at visible wavelengths. Binarity in future high precision calibrations of the P-L relations is not negligible, at least in the short-period regime. Subtracting the companion flux leaves V1334 Cyg in marginal agreement with existing photometric-based P-L relations, indicating either an overall calibration bias or a significant intrinsic dispersion at a few percent level. Our work also enabled us to determine the dynamical masses of both components, M1 = 4.288 +/- 0.133 Msun (Cepheid) and M2 = 4.040 +/- 0.048 Msun (companion), providing the most accurate masses for a Galactic binary Cepheid system.

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.

Archives of long photometric surveys, like the Kepler database, are a gold mine for studying flares. However, identifying them is a complex task; while in the case of single-target observations it can be easily done manually by visual inspection, this is nearly impossible for years-long time series for several thousand targets. Although there exist automated methods for this task, several problems are difficult (or impossible) to overcome with traditional fitting and analysis approaches. We introduce a code for identifying and analyzing flares based on machine learning methods, which are intrinsically adept at handling such data sets. We used the RANSAC (RANdom SAmple Consensus) algorithm to model light curves, as it yields robust fits even in case of several outliers, like flares. The light curve is divided into search windows, approximately in the order of the stellar rotation period. This search window is shifted over the data set, and a voting system is used to keep false positives to a minimum: only those flare candidate points are kept that were identified in several windows as a flare. The code was tested on the K2 observations of the TRAPPIST-1, and on the long cadence data of KIC 1722506. The detected flare events and flare energies are consistent with earlier results from manual inspections.

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 $<10$ $\mu$-arcsecond precision on differential position with 10-minute observations. The high precision of our astrometric orbit suggests that exoplanet detection capabilities are within reach of MIRC at CHARA. We compute exoplanet detection limits around $\delta$ Del, and conclude that if this precision is extended to wider systems we should be able to detect most exoplanets $>2$ M$_{J}$ on orbits $>0.75$ AU around individual components of hot binary stars via differential astrometry.

The signatures of planets hosted by M dwarfs are more readily detected with transit photometry and radial velocity methods than those of planets around larger stars. Recently, transit photometry was used to discover seven planets orbiting the late-M dwarf TRAPPIST-1. Three of TRAPPIST-1's planets fall in the Habitable Zone, a region where liquid water could exist on the planetary surface given appropriate planetary conditions. We aim to investigate the habitability of the TRAPPIST-1 planets by studying the star's activity and its effect on the planets. We analyze previously-published space- and ground-based light curves and show the photometrically-determined rotation period of TRAPPIST-1 appears to vary over time due to complicated, evolving surface activity. The dramatic changes of the surface of TRAPPIST-1 suggest that rotation periods determined photometrically may not be reliable for this and similarly active stars. While the activity of the star is low, we use the premise of the "cosmic shoreline" to provide evidence that the TRAPPIST-1 environment has potentially led to the erosion of possible planetary atmospheres by extreme ultraviolet stellar emission.

Nearby, active stars with relatively rapid rotation and large starspot structures offer the opportunity to compare interferometric, spectroscopic, and photometric imaging techniques. In this paper, we image a spotted star with three different methods for the first time. The giant primary star of the RS Canum Venaticorum binary $\sigma$ Geminorum ($\sigma$ Gem) was imaged for two epochs of interferometric, high-resolution spectroscopic, and photometric observations. The light curves from the reconstructions show good agreement with the observed light curves, supported by the longitudinally-consistent spot features on the different maps. However, there is strong disagreement in the spot latitudes across the methods.

Sunspots are cool areas caused by strong surface magnetic fields inhibiting convection. Moreover, strong magnetic fields can alter the average atmospheric structure, degrading our ability to measure stellar masses and ages. Stars more active than the Sun have more and stronger dark spots than in the solar case, including on the rotational pole itself. Doppler imaging, which has so far produced the most detailed images of surface structures on other stars than the Sun, cannot always distinguish the hemisphere in which the starspots are located, especially in the equatorial region and if the data quality is not optimal. This leads to problems in investigating the north-south distribution of starspot active latitudes (those latitudes with more spot activity), which are crucial constraints of dynamo theory. Polar spots, inferred only from Doppler tomography, could plausibly be observational artifacts, casting some doubt on their very existence. Here we report imaging of the old, magnetically-active star $\zeta$ Andromedae using long-baseline infrared interferometry. In our data, a dark polar spot is seen in each of two epochs, while lower-latitude spot structures in both hemispheres do not persist between observations revealing global starspot asymmetries. The north-south symmetry of active latitudes observed on the Sun is absent on $\zeta$ And, which hosts global spot patterns that cannot be produced by solar-type dynamos.

Magnetic fields are essential ingredients of many physical processes in the interiors and envelopes of cool stars. Yet their direct detection and characterisation is notoriously difficult, requiring high-quality observations and advanced analysis techniques. Significant progress has been recently achieved by several types of direct magnetic field studies on the surfaces of cool active stars. In particular, complementary techniques of the field topology mapping with polarisation data and total magnetic flux measurements from intensity spectra have been systematically applied to different classes of active stars leading to interesting and occasionally controversial results. In this paper we summarise the current status of direct magnetic field studies of cool stars, and investigations of surface inhomogeneities caused by the field, based on the material presented at the Cool Stars 19 splinter session.

Using the high-precision photometry from the Kepler space telescope, thousands of stars with stellar and planetary companions have been observed. The characterization of stars with companions is not always straightforward and can be contaminated by systematic and stellar influences on the light curves. Here, through a detailed analysis of starspots and eclipses, we identify KOI-1003 as a new, active RS CVn star---the first identified with data from Kepler. The Kepler light curve of this close binary system exhibits the system's primary transit, secondary eclipse, and starspot evolution of two persistent active longitudes. The near equality of the system's orbital and rotation periods indicates the orbit and primary star's rotation are nearly synchronized ($P_\mathrm{orb} = 8.360613\pm0.000003$ days; $P_\mathrm{rot} \sim 8.23$ days). By assuming the secondary star is on the main sequence, we suggest the system consists of a $1.45^{+0.11}_{-0.19} \ M_\odot$ subgiant primary and a $0.59^{+0.03}_{-0.04} \ M_\odot$ main-sequence companion. Our work gives a distance of $4400 \pm 600$ pc and an age of $t = 3.0^{-0.5}_{+2.0}$ Gyr, parameters which are discrepant with previous studies that included the star as a member of the open cluster NGC 6791.

We report new CHARA/MIRC interferometric observations of the Cepheid archetype $\delta$ Cep, which aimed at detecting the newly discovered spectroscopic companion. We reached a maximum dynamic range $\Delta H $ = 6.4, 5.8, and 5.2 mag, respectively within the relative distance to the Cepheid $r < 25$ mas, $25 < r < 50$ mas and $50 < r < 100$ mas. Our observations did not show strong evidence of a companion. We have a marginal detection at $3\sigma$ with a flux ratio of 0.21\%, but nothing convincing as we found other possible probable locations. We ruled out the presence of companion with a spectral type earlier than F0V, A1V and B9V, respectively for the previously cited ranges $r$. From our estimated sensitivity limits and the Cepheid light curve, we derived lower-limit magnitudes in the $H$ band for this possible companion to be $H_\mathrm{comp} > 9.15, 8.31$ and 7.77 mag, respectively for $r < 25$ mas, $25 < r < 50$ mas and $50 < r < 100$ mas. We also found that to be consistent with the predicted orbital period, the companion has to be located at a projected separation $< 24$ mas with a spectral type later than a F0V star.

To measure the stellar and orbital properties of the metal-poor RS CVn binary o Draconis (o Dra), we directly detect the companion using interferometric observations obtained with the Michigan InfraRed Combiner at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array. The H-band flux ratio between the primary and secondary stars is the highest confirmed flux ratio (370 +/- 40) observed with long-baseline optical interferometry. These detections are combined with radial velocity data of both the primary and secondary stars, including new data obtained with the Tillinghast Reflector Echelle Spectrograph on the Tillinghast Reflector at the Fred Lawrence Whipple Observatory and the 2-m Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We determine an orbit from which we find model-independent masses and ages of the components (M_A = 1.35 +\- 0.05 M_Sun, M_B = 0.99 +\- 0.02 M_Sun, system age = 3.0 -\+ 0.5 Gyr). An average of a 23-year light curve of o Dra from the Tennessee State University Automated Photometric Telescope folded over the orbital period newly reveals eclipses and the quasi-sinusoidal signature of ellipsoidal variations. The modeled light curve for our system's stellar and orbital parameters confirm these ellipsoidal variations due to the primary star partially filling its Roche lobe potential, suggesting most of the photometric variations are not due to stellar activity (starspots). Measuring gravity darkening from the average light curve gives a best-fit of beta = 0.07 +\- 0.03, a value consistent with conventional theory for convective envelope stars. The primary star also exhibits an anomalously short rotation period, which, when taken with other system parameters, suggests the star likely engulfed a low-mass companion that had recently spun-up the star.

Long-baseline interferometry is an important technique to spatially resolve binary or multiple systems in close orbits. By combining several telescopes together and spectrally dispersing the light, it is possible to detect faint components around bright stars. Aims. We provide a rigorous and detailed method to search for high-contrast companions around stars, determine the detection level, and estimate the dynamic range from interferometric observations. We developed the code CANDID (Companion Analysis and Non-Detection in Interferometric Data), a set of Python tools that allows us to search systematically for point-source, high-contrast companions and estimate the detection limit. The search pro- cedure is made on a N x N grid of fit, whose minimum needed resolution is estimated a posteriori. It includes a tool to estimate the detection level of the companion in the number of sigmas. The code CANDID also incorporates a robust method to set a 3{\sigma} detection limit on the flux ratio, which is based on an analytical injection of a fake companion at each point in the grid. We used CANDID to search for the companions around the binary Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul. First, we showed that our previous discoveries of the components orbiting V1334 Cyg and AX Cir were detected at > 13 sigmas. The companion around AW Per is detected at more than 15 sigmas with a flux ratio of f = 1.22 +/- 0.30 %. We made a possible detection of the companion orbiting RT Aur with f = 0.22 +/- 0.11 %. It was detected at 3.8{\sigma} using the closure phases only, and so more observations are needed to confirm the detection. We also set the detection limit for possible undetected companions. We found that there is no companion with a spectral type earlier than B7V, A5V, F0V, B9V, A0V, and B9V orbiting V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul, respectively.

A classical nova occurs when material accreting onto the surface of a white dwarf in a close binary system ignites in a thermonuclear runaway. Complex structures observed in the ejecta at late stages could result from interactions with the companion during the common envelope phase. Alternatively, the explosion could be intrinsically bipolar, resulting from a localized ignition on the surface of the white dwarf or as a consequence of rotational distortion. Studying the structure of novae during the earliest phases is challenging because of the high spatial resolution needed to measure their small sizes. Here we report near-infrared interferometric measurements of the angular size of Nova Delphini 2013, starting from one day after the explosion and continuing with extensive time coverage during the first 43 days. Changes in the apparent expansion rate can be explained by an explosion model consisting of an optically thick core surrounded by a diffuse envelope. The optical depth of the ejected material changes as it expands. We detect an ellipticity in the light distribution, suggesting a prolate or bipolar structure that develops as early as the second day. Combining the angular expansion rate with radial velocity measurements, we derive a geometric distance to the nova of 4.54 +/- 0.59 kpc from the Sun.

To measure the properties of both components of the RS CVn binary sigma Geminorum (sigma Gem), we directly detect the faint companion, measure the orbit, obtain model-independent masses and evolutionary histories, detect ellipsoidal variations of the primary caused by the gravity of the companion, and measure gravity darkening. We detect the companion with interferometric observations obtained with the Michigan InfraRed Combiner (MIRC) at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array with a primary-to-secondary H-band flux ratio of 270+/-70. A radial velocity curve of the companion was obtained with spectra from the Tillinghast Reflector Echelle Spectrograph (TRES) on the 1.5-m Tillinghast Reflector at Fred Lawrence Whipple Observatory (FLWO). We additionally use new observations from the Tennessee State University Automated Spectroscopic and Photometric Telescopes (AST and APT, respectively). From our orbit, we determine model-independent masses of the components (M_1 = 1.28 +/- 0.07 M_Sun, M_2 = 0.73 +/- 0.03 M_Sun), and estimate a system age of 5 -/+ 1 Gyr. An average of the 27-year APT light curve of sigma Gem folded over the orbital period (P = 19.6027 +/- 0.0005 days) reveals a quasi-sinusoidal signature, which has previously been attributed to active longitudes 180 deg apart on the surface of sigma Gem. With the component masses, diameters, and orbit, we find that the predicted light curve for ellipsoidal variations due to the primary star partially filling its Roche lobe potential matches well with the observed average light curve, offering a compelling alternative explanation to the active longitudes hypothesis. Measuring gravity darkening from the light curve gives beta < 0.1, a value slightly lower than that expected from recent theory.

We present the results of long-baseline optical interferometry observations using the Precision Astronomical Visual Observations (PAVO) beam combiner at the Center for High Angular Resolution Astronomy (CHARA) Array to measure the angular sizes of three bright Kepler stars: {\theta} Cygni, and both components of the binary system 16 Cygni. Supporting infrared observations were made with the Michigan Infrared Combiner (MIRC) and Classic beam combiner, also at the CHARA Array. We find limb-darkened angular diameters of 0.753+/-0.009 mas for {\theta} Cyg, 0.539+/-0.007 mas for 16 Cyg A and 0.490+/-0.006 mas for 16 Cyg B. The Kepler Mission has observed these stars with outstanding photometric precision, revealing the presence of solar-like oscillations. Due to the brightness of these stars the oscillations have exceptional signal-to-noise, allowing for detailed study through asteroseismology, and are well constrained by other observations. We have combined our interferometric diameters with Hipparcos parallaxes, spectrophotometric bolometric fluxes and the asteroseismic large frequency separation to measure linear radii ({\theta} Cyg: 1.48+/-0.02 Rsun, 16 Cyg A: 1.22+/-0.02 Rsun, 16 Cyg B: 1.12+/-0.02 Rsun), effective temperatures ({\theta} Cyg: 6749+/-44 K, 16 Cyg A: 5839+/-42 K, 16 Cyg B: 5809+/-39 K), and masses ({\theta} Cyg: 1.37+/-0.04 Msun, 16 Cyg A: 1.07+/-0.05 Msun, 16 Cyg B: 1.05+/-0.04 Msun) for each star with very little model dependence. The measurements presented here will provide strong constraints for future stellar modelling efforts.

The Kepler target KIC 5110407, a K-type star, shows strong quasi-periodic light curve fluctuations likely arising from the formation and decay of spots on the stellar surface rotating with a period of 3.4693 days. Using an established light-curve inversion algorithm, we study the evolution of the surface features based on Kepler space telescope light curves over a period of two years (with a gap of .25 years). At virtually all epochs, we detect at least one large spot group on the surface causing a 1-10% flux modulation in the Kepler passband. By identifying and tracking spot groups over a range of inferred latitudes, we measured the surface differential rotation to be much smaller than that found for the Sun. We also searched for a correlation between the seventeen stellar flares that occurred during our observations and the orientation of the dominant surface spot at the time of each flare. No statistically-significant correlation was found except perhaps for the very brightest flares, suggesting most flares are associated with regions devoid of spots or spots too small to be clearly discerned using our reconstruction technique. While we may see hints of long-term changes in the spot characteristics and flare statistics within our current dataset, a longer baseline of observation will be needed to detect the existence of a magnetic cycle in KIC 5110407.

We have discovered a dwarf nova (DN) of type SU UMa in Kepler data which is 7.0 arcsec from the G-type exoplanet survey target KIC 4378554. The DN appears as a background source in the pixel aperture of the foreground G star. We extracted only the pixels where the DN is present and observed the source to undergo five outbursts -- one a superoutburst -- over a timespan of 22 months. The superoutburst was triggered by a normal outburst, a feature that has been seen in all DNe superoutburst observed by Kepler. Superhumps during the super outburst had a period of 1.842+/-0.004 h and we see a transition from disc-dominated superhump signal to a mix of disc and accretion stream impact. Predictions of the number of DNe present in Kepler data based on previously published space densities vary from 0.3 to 258. An investigation of the background pixels targets would lead to firmer constraints on the space density of DN.

The recently discovered gamma-ray binary 1FGL J1018.6-5856 has a proposed optical/near-infrared (OIR) counterpart 2MASS 10185560-5856459. We present Stromgren photometry of this star to investigate its photometric variability and measure the reddening and distance to the system. We find that the gamma-ray binary has E(B-V) = 1.34 +/- 0.04 and d = 5.4^+4.6_-2.1 kpc. While E(B-V) is consistent with X-ray observations of the neutral hydrogen column density, the distance is somewhat closer than some previous authors have suggested.

We present the results of a study of differential rotation on the K2 IV primary of the RS CVn binary II Pegasi (HD 224085) performed by inverting light curves to produce images of the dark starspots on its surface. The data were obtained in the standard Johnson B and V filter passbands via the Tennessee State University T3 0.4-m Automated Photometric Telescope from JD 2447115.8086 to 2454136.6221 (1987 November 16 to 2007 February 5). The observations were subdivided into 68 data sets consisting of pairs of B and V light curves, which were then inverted using a constrained non-linear inversion algorithm that makes no a priori assumptions regarding the number of spots or their shapes. The resulting surface images were then assigned to 21 groups corresponding to time intervals over which we could observe the evolution of a given group of spots (except for three groups consisting of single data sets). Of these 21 groups, six showed convincing evidence of differential rotation over time intervals of several months. For the others, the spot configuration was such that differential rotation was neither exhibited nor contraindicated. The differential rotation we infer is in the same sense as that on the Sun: lower latitudes have shorter rotation periods. From plots of the range in longitude spanned by the spotted regions vs. time, we obtain estimates of the differential rotation coefficient k defined in earlier work by Henry et al., and show that our results for its value are consistent with the value obtained therein.

HD 15137 is an intriguing runaway O-type binary system that offers a rare opportunity to explore the mechanism by which it was ejected from the open cluster of its birth. Here we present recent blue optical spectra of HD 15137 and derive a new orbital solution for the spectroscopic binary and physical parameters of the O star primary. We also present the first XMM-Newton observations of the system. Fits of the EPIC spectra indicate soft, thermal X-ray emission consistent with an isolated O star. Upper limits on the undetected hard X-ray emission place limits on the emission from a proposed compact companion in the system, and we rule out a quiescent neutron star in the propellor regime or a weakly accreting neutron star. An unevolved secondary companion is also not detected in our optical spectra of the binary, and it is difficult to conclude that a gravitational interaction could have ejected this runaway binary with a low mass optical star. HD 15137 may contain an elusive neutron star in the ejector regime or a quiescent black hole with conditions unfavorable for accretion at the time of our observations.

Non-radial pulsations (NRPs) are a proposed mechanism for the formation of decretion disks around Be stars and are important tools to study the internal structure of stars. NGC 3766 has an unusually large fraction of transient Be stars, so it is an excellent location to study the formation mechanism of Be star disks. High resolution spectroscopy can reveal line profile variations from NRPs, allowing measurements of both the degree, l, and azimuthal order, m. However, spectroscopic studies require large amounts of time with large telescopes to achieve the necessary high S/N and time domain coverage. On the other hand, multi-color photometry can be performed more easily with small telescopes to measure l only. Here, we present representative light curves of Be stars and non-emitting B stars in NGC 3766 from the CTIO 0.9m telescope in an effort to study NRPs in this cluster.

LS I +61 303 and LS 5039 are two of only a handful of known high mass X-ray binaries (HMXBs) that exhibit very high energy emission in the MeV-TeV range, and these "gamma-ray binaries" are of renewed interest due to the recent launch of the Fermi Gamma-ray Space Telescope. Here we present new radial velocities of both systems based on recent red and blue optical spectra. Both systems have somewhat discrepant orbital solutions available in the literature, and our new measurements result in improved orbital elements and resolve the disagreements. The improved geometry of each orbit will aid in studies of the high energy emission region near each source.