We present results from a 2-epoch HST H$\alpha$ emission line survey of the Andromeda Galaxy that overlaps the footprint of the Panchromatic Hubble Andromeda Treasury (PHAT) survey. We find 552 (542) classical Be stars and 8429 (8556) normal B-type stars in epoch # 1 (epoch # 2), yielding an overall fractional Be content of 6.15% $\pm$0.26% (5.96% $\pm$0.25%). The fractional Be content decreased with spectral sub-type from $\sim$23.6% $\pm$2.0% ($\sim$23.9% $\pm$2.0%) for B0-type stars to $\sim$3.1% $\pm$0.34% ($\sim$3.4% $\pm$0.35%) for B8-type stars in epoch # 1 (epoch # 2). We observe a clear population of cluster Be stars at early fractional main sequence lifetimes, indicating that a subset of Be stars emerge onto the ZAMS as rapid rotators. Be stars are 2.8x rarer in M31 for the earliest sub-types compared to the SMC, confirming that the fractional Be content decreases in significantly more metal rich environments (like the Milky Way and M31). However, M31 does not follow a clear trend of Be fraction decreasing with metallicity compared to the Milky Way, which may reflect that the Be phenomenon is enhanced with evolutionary age. The rate of disk-loss or disk-regeneration episodes we observe, 22% $\pm$ 2% yr$^{-1}$, is similar to that observed for seven other Galactic clusters reported in the literature, assuming these latter transient fractions scale by a linear rate. The similar number of disk-loss events (57) as disk-renewal events (43) was unexpected since disk dissipation time-scales can be $\sim$2x the typical time-scales for disk build-up phases.

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.

We present coronagraphic long slit spectra of AU Mic's debris disk taken with the STIS instrument aboard the Hubble Space Telescope (HST). Our spectra are the first spatially resolved, scattered light spectra of the system's disk, which we detect at projected distances between approximately 10 and 45 AU. Our spectra cover a wavelength range between 5200 and 10200 angstroms. We find that the color of AU Mic's debris disk is bluest at small (12-35 AU) projected separations. These results both confirm and quantify the findings qualitatively noted by Krist et al. (2005), and are different than IR observations that suggested a uniform blue or gray color as a function of projected separation in this region of the disk. Unlike previous literature that reported the color of AU Mic's disk became increasingly more blue as a function of projected separation beyond approximately 30 AU, we find the disk's optical color between 35-45 AU to be uniformly blue on the southeast side of the disk and decreasingly blue on the northwest side. We note that this apparent change in disk color at larger projected separations coincides with several fast, outward moving "features" that are passing through this region of the southeast side of the disk. We speculate that these phenomenon might be related, and that the fast moving features could be changing the localized distribution of sub-micron sized grains as they pass by, thereby reducing the blue color of the disk in the process. We encourage follow-up optical spectroscopic observations of the AU Mic to both confirm this result, and search for further modifications of the disk color caused by additional fast moving features propagating through the disk.

We analyze 45 spectropolarimetric observations of the eclipsing, interacting binary star V356 Sgr, obtained over a period of 21 years, to characterize the geometry of the system's circumstellar material. After removing interstellar polarization from these data, we find the system exhibits a large intrinsic polarization signature arising from electron scattering. In addition, the lack of repeatable eclipses in the polarization phase curves indicates the presence of a substantial pool of scatterers not occulted by either star. We suggest that these scatterers form either a circumbinary disk coplanar with the gainer's accretion disk or an elongated structure perpendicular to the orbital plane of V356 Sgr, possibly formed by bipolar outflows. We also observe small-scale, cycle-to-cycle variations in the magnitude of intrinsic polarization at individual phases, which we interpret as evidence of variability in the amount of scattering material present within and around the system. This may indicate a mass transfer or mass loss rate that varies on the time-scale of the system's orbital period. Finally, we compare the basic polarimetric properties of V356 Sgr with those of the well studied beta Lyr system; the significant differences observed between the two systems suggests diversity in the basic circumstellar geometry of Roche-lobe overflow systems

We present new analysis of multi-epoch, H-band, scattered light images of the AB Aur system. We used a Monte Carlo, radiative transfer code to simultaneously model the system's SED and H-band polarized intensity imagery. We find that a disk-dominated model, as opposed to one that is envelope dominated, can plausibly reproduce AB Aur's SED and near-IR imagery. This is consistent with previous modeling attempts presented in the literature and supports the idea that at least a subset of AB Aur's spirals originate within the disk. In light of this, we also analyzed the movement of spiral structures in multi-epoch H-band total light and polarized intensity imagery of the disk. We detect no significant rotation or change in spatial location of the spiral structures in these data, which span a 5.8 year baseline. If such structures are caused by disk-planet interactions, the lack of observed rotation constrains the location of the orbit of planetary perturbers to be >47 AU.

Because most massive stars have been or will be affected by a companion during the course of their evolution, we cannot afford to neglect binaries when discussing the progenitors of supernovae and GRBs. Analyzing linear polarization in the emission lines of close binary systems allows us to probe the structures of these systems' winds and mass flows, making it possible to map the complex morphologies of the mass loss and mass transfer structures that shape their subsequent evolution. In Wolf-Rayet (WR) binaries, line polarization variations with orbital phase distinguish polarimetric signatures arising from lines that scatter near the stars from those that scatter far from the orbital plane. These far-scattering lines may form the basis for a "binary line-effect method" of identifying rapidly rotating WR stars (and hence GRB progenitor candidates) in binary systems.

The SDSS III APOGEE survey recently identified two new $\sigma$ Ori E type candidates, HD 345439 and HD 23478, which are a rare subset of rapidly rotating massive stars whose large (kGauss) magnetic fields confine circumstellar material around these systems. Our analysis of multi-epoch photometric observations of HD 345439 from the KELT, SuperWASP, and ASAS surveys reveals the presence of a $\sim$0.7701 day period in each dataset, suggesting the system is amongst the faster known $\sigma$ Ori E analogs. We also see clear evidence that the strength of H-alpha, H I Brackett series lines, and He I lines also vary on a $\sim$0.7701 day period from our analysis of multi-epoch, multi-wavelength spectroscopic monitoring of the system from the APO 3.5m telescope. We trace the evolution of select emission line profiles in the system, and observe coherent line profile variability in both optical and infrared H I lines, as expected for rigidly rotating magnetosphere stars. We also analyze the evolution of the H I Br-11 line strength and line profile in multi-epoch observations of HD 23478 from the SDSS-III APOGEE instrument. The observed periodic behavior is consistent with that recently reported by Sikora and collaborators in optical spectra.

We present an overview of four phase-constrained Chandra HETGS X-ray observations of Delta Ori A. Delta Ori A is actually a triple system which includes the nearest massive eclipsing spectroscopic binary, Delta Ori Aa, the only such object which can be observed with little phase-smearing with the Chandra gratings. Since the fainter star, Delta Ori Aa2, has a much lower X-ray luminosity than the brighter primary, Delta Ori A provides a unique system with which to test the spatial distribution of the X-ray emitting gas around Delta Ori Aa1 via occultation by the photosphere of and wind cavity around the X-ray dark secondary. Here we discuss the X-ray spectrum and X-ray line profiles for the combined observation, having an exposure time of nearly 500 ksec and covering nearly the entire binary orbit. Companion papers discuss the X-ray variability seen in the Chandra spectra, present new space-based photometry and ground-based radial velocities simultaneous with the X-ray data to better constrain the system parameters, and model the effects of X-rays on the optical and UV spectrum. We find that the X-ray emission is dominated by embedded wind shock emission from star Aa1, with little contribution from the tertiary star Ab or the shocked gas produced by the collision of the wind of Aa1 against the surface of Aa2. We find a similar temperature distribution to previous X-ray spectrum analyses. We also show that the line half-widths are about $0.3-0.5\times$ the terminal velocity of the wind of star Aa1. We find a strong anti-correlation between line widths and the line excitation energy, which suggests that longer-wavelength, lower-temperature lines form farther out in the wind. Our analysis also indicates that the ratio of the intensities of the strong and weak lines of \ion{Fe}{17} and \ion{Ne}{10} are inconsistent with model predictions, which may be an effect of resonance scattering

We present time-resolved and phase-resolved variability studies of an extensive X-ray high-resolution spectral dataset of the $\delta$ Orionis Aa binary system. The four observations, obtained with Chandra ACIS HETGS, have a total exposure time of ~479 ks and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range 5-25 $\AA$ is confirmed, with maximum amplitude of about +/-15% within a single ~125 ks observation. Periods of 4.76d and 2.04d are found in the total X-ray flux, as well as an apparent overall increase in flux level throughout the 9-day observational campaign. Using 40 ks contiguous spectra derived from the original observations, we investigate variability of emission line parameters and ratios. Several emission lines are shown to be variable, including S XV, Si XIII, and Ne IX. For the first time, variations of the X-ray emission line widths as a function of the binary phase are found in a binary system, with the smallest widths at phase=0.0 when the secondary $\delta$ Orionis Aa2 is at inferior conjunction. Using 3D hydrodynamic modeling of the interacting winds, we relate the emission line width variability to the presence of a wind cavity created by a wind-wind collision, which is effectively void of embedded wind shocks and is carved out of the X-ray-producing primary wind, thus producing phase-locked X-ray variability.

We report on both high-precision photometry from the MOST space telescope and ground-based spectroscopy of the triple system delta Ori A consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7d, and a more distant tertiary (O9 IV P > 400 yrs). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for 3 weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy we have a well constrained radial velocity curve of Aa1. While we are unable to recover radial velocity variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries which show evidence for tidally induced pulsations.

We present results from a study of the eclipsing, colliding-wind binary V444 Cyg that uses a combination of X-ray and optical spectropolarimetric methods to describe the 3-D nature of the shock and wind structure within the system. We have created the most complete X-ray light curve of V444 Cyg to date using 40 ksec of new data from Swift, and 200 ksec of new and archived XMM-Newton observations. In addition, we have characterized the intrinsic, polarimetric phase-dependent behavior of the strongest optical emission lines using data obtained with the University of Wisconsin's Half-Wave Spectropolarimeter. We have detected evidence of the Coriolis distortion of the wind-wind collision in the X-ray regime, which manifests itself through asymmetric behavior around the eclipses in the system's X-ray light curves. The large opening angle of the X-ray emitting region, as well as its location (i.e. the WN wind does not collide with the O star, but rather its wind) are evidence of radiative braking/inhibition occurring within the system. Additionally, the polarimetric results show evidence of the cavity the wind-wind collision region carves out of the Wolf-Rayet star's wind.

We present results from six years of recalibrated and new spectropolarimetric data taken with the University of Wisconsin's Half-Wave Spectropolarimeter (HPOL) and six years of new data taken with the photoelastic modulating polarimeter (PEMP) at the Flower and Cook Observatory. Combining these data with polarimetric data from the literature allows us to characterize the intrinsic BVRI polarized light curves. A repeatable discrepancy of 0.245 days (approximately 6 hours) between the secondary minima in the total light curve and the polarization curve in the V band, with similar behavior in the other bands, may represent the first direct evidence for an accretion hot spot on the disk edge.

Beta Lyrae is an eclipsing, semi-detached binary system whose state of active mass transfer can reveal details of the nonconservative evolution of binary stars. Roche lobe overflow has caused the system to evolve to a complex state. A thick accretion disk almost completely obscures the secondary, mass-gaining star while the rapid mass transfer likely drives mass loss through the system's bipolar outflows. Polarimetry can provide important information about the physical structure of complex systems; in fact, the discovery of bipolar outflows in beta Lyrae was confirmed through polarimetry. Here we present results from 6 years of new and recalibrated spectropolarimetric data taken with the University of Wisconsin's Half-Wave Spectropolarimeter (HPOL). We discuss their implications for our current understanding of the system's disk-jet geometry. Using both broadband and line polarization analysis techniques, we derive new information about the structure of the disk, the presence and location of a hot spot, and the distribution of hot line-emitting gas.