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.

There is growing evidence that some Be stars were spun up through mass transfer in a close binary system, leaving the former mass donor star as a hot, stripped-down object. There are five known cases of Be stars with hot subdwarf (sdO) companions that were discovered through International Ultraviolet Explorer (IUE) spectroscopy. Here we expand the search for Be+sdO candidates using archival FUV spectra from IUE. We collected IUE spectra for 264 stars and formed cross-correlation functions (CCFs) with a model spectrum for a hot subdwarf. Twelve new candidate Be+sdO systems were found, and eight of these display radial velocity variations associated with orbital motion. The new plus known Be+sdO systems have Be stars with spectral subtypes of B0 to B3, and the lack of later-type systems is surprising given the large number of cooler B-stars in our sample. We discuss explanations for the observed number and spectral type distribution of the Be+sdO systems, and we argue that there are probably many Be systems with stripped companions that are too faint for detection through our analysis.

We used archival International Ultraviolet Explorer (IUE) high-dispersion, short wavelength spectra data to search for evidence of the spectra of hot subdwarf companions of six rapidly rotating Be stars in binary systems. We searched for the signature of a hot companion through an analysis of the cross-correlation functions of observed and model spectra that were separated into primary and secondary components using a Doppler tomography algorithm and adopted spectroscopic orbital solutions. A positive detection of the flux from a hot companion was made for the reconstructed secondary cross-correlation function of just one target, 60 Cygni (B1 Ve). We estimate that the companion of the Be star in 60 Cygni has $T_{\rm eff} = 42 \pm 4$ kK, mass ratio $M_2/M_1 = 0.15 \pm 0.02$, and monochromatic flux ratio ${f_2}/{f_1} = 0.034 \pm 0.002$ in the spectral region near 1525 \AA. If the companions of the other target Be stars are also hot, then they must be faint and contribute less than $\approx 1\%$ of the UV flux ($<0.6 \%$ in the case of $\gamma$ Cas). We also discuss in an appendix a shell episode of Pleione (28 Tau) recorded in the IUE spectra.

We present a spectroscopic investigation of the Be+sdO binary system HR 2142 that is based upon large sets of ultraviolet observations from the International Ultraviolet Explorer and ground-based H-alpha observations. We measured radial velocities for the Be star component from these spectra, and computed a revised orbit. In order to search for the spectral signature of the hot subdwarf, we cross-correlated the short wavelength end of each IUE spectrum with a model hot star spectrum, and then we used the predicted Doppler shifts of the subdwarf to shift-and-add all the cross-correlation functions to the frame of the subdwarf. This merged function shows the weak signal from the spectral lines of the hot star, and a best fit is obtained with a mass ratio M_2/M_1 = 0.07 +/- 0.02, companion temperature T_{eff} > 43 +/- 5 kK, projected rotational velocity V sin i < 30 km/s, and a monochromatic flux ratio near 1170 Angstroms of f_2/f_1 > 0.009 +/- 0.001. This hot subdwarf creates a one-armed spiral, tidal wake in the disk of the Be star, and we present a circumbinary disk model that can explain the occurrence of shell absorption lines by gas enhancements that occur where gas crossing the gap created by the subdwarf strikes the disk boundaries. The faint companion of HR 2142 may be representative of a significant fraction of Be stars with undetected former mass donor companion stars.

Many short-period binary stars have distant orbiting companions that have played a role in driving the binary components into close separation. Indirect detection of a tertiary star is possible by measuring apparent changes in eclipse times of eclipsing binaries as the binary orbits the common center of mass. Here we present an analysis of the eclipse timings of 41 eclipsing binaries observed throughout the NASA Kepler mission of long duration and precise photometry. This subset of binaries is characterized by relatively deep and frequent eclipses of both stellar components. We present preliminary orbital elements for seven probable triple stars among this sample, and we discuss apparent period changes in seven additional eclipsing binaries that may be related to motion about a tertiary in a long period orbit. The results will be used in ongoing investigations of the spectra and light curves of these binaries for further evidence of the presence of third stars.

We report on the detection of a hot subdwarf component in the Be binary system, 59 Cygni. The spectral signature is found in cross-correlation functions of photospheric model spectra with far-ultraviolet spectra obtained by the International Ultraviolet Explorer Satellite, and we used radial velocities from the cross-correlation functions to determine a double-lined spectroscopic orbit. The individual spectra of the binary components were extracted using a Doppler tomography algorithm. The flux of the system is dominated by the rapidly rotating Be star. However, the subdwarf contributes approximately 4% of the UV flux, and its spectrum bears a strong resemblance to that of the hot sdO star BD +75 325. Based upon the appearance of the UV spectrum and the orbital elements, we present estimates for the stellar masses, radii, and temperatures. The presence of the hot companion causes excess emission from the outer part of the Be disk facing the companion. We present a set of red spectra that show the orbital phase variations of the He I 6678 emission formed in the heated region of the disk, which probably occurs near the disk outer boundary. 59 Cygni, FY Canis Majoris, and phi Persei comprise the known set of Be binaries with detected hot evolved companions, which are the stripped down remains of mass transfer. Their properties demonstrate that some fraction of Be stars were spun up through angular momentum transfer by Roche lobe overflow.

AU Mon is a long-period (11.113 d) Algol-type binary system with a persistent accretion disk that is apparent as double-peaked H-alpha emission. We present previously unpublished optical spectra of AU Mon which were obtained over several years with dense orbital phase coverage. We utilize these data, along with archival UV spectra, to model the temperature and structure of the accretion disk and the gas stream. Synthetic spectral profiles for lines including H-alpha, H-beta, and the Al III and Si IV doublets were computed with the Shellspec program. The best match between the model spectra and the observations is obtained for an accretion disk of inner/outer radius 5.1/23 R_sun, thickness of 5.2 R_sun, density of 1.0e-13 g/cm^3, and maximum temperature of 14000 K, along with a gas stream at a temperature of ~8000 K transferring ~2.4e-9 M_sun/yr. We show H-alpha Doppler tomograms of the velocity structure of the gas, constructed from difference profiles calculated through sequentially subtracting contributions from the stars and accretion structures. The tomograms provide independent support for the Shellspec modeling, while also illustrating that residual emission at sub-Keplerian velocities persists even after subtracting the disk and stream emission. Spectral variability in the H-alpha profile beyond that expected from either the orbital or the long-period cycle is present on both multi-week and multi-year timescales, and may reflect quasi-random changes in the mass transfer rate or the disk structure. Finally, a transient UV spectral absorption feature may be modeled as an occasional outflow launched from the vicinity of the disk-stream interaction region.

V393 Scorpii is a Double Periodic Variable characterized by a relatively stable non-orbital photometric cycle of 253 days. Mennickent et al. argue for the presence of a massive optically thick disc around the more massive B-type component and describe the evolutionary stage of the system. In this paper we analyze the behavior of the main spectroscopic optical lines during the long non-orbital photometric cycle. We study the radial velocity of the donor determining their orbital elements and find a small but significant orbital eccentricity (e = 0.04). The donor spectral features are modeled and removed from the spectrum at every observing epoch using the light-curve model given by Mennickent et al. We find that the line emission is larger during eclipses and mostly comes from a bipolar wind. We find that the long cycle is explained in terms of a modulation of the wind strength; the wind has a larger line and continuum emissivity on the high state. We report the discovery of highly variable chromospheric emission in the donor, as revealed by Doppler maps of the emission lines MgII 4481 and CI 6588. We discuss notable and some novel spectroscopic features like discrete absorption components, especially visible at blue-depressed OI 7773 absorption wings during the second half-cycle, Balmer double emission with V/R-curves showing "Z-type" and "S-type" excursions around secondary and main eclipse, respectively, and H_beta emission wings extending up to +- 2000 km/s. We discuss possible causes for these phenomena and for their modulations with the long cycle.

We present the first results of a Kepler survey of 41 eclipsing binaries that we undertook to search for third star companions. Such tertiaries will periodically alter the eclipse timings through light travel time and dynamical effects. We discuss the prevalence of starspots and pulsation among these binaries and how these phenomena influence the eclipse times. There is no evidence of short period companions (P < 700 d) among this sample, but we do find evidence for long term timing variations in 14 targets (34%). We argue that this finding is consistent with the presence of tertiary companions among a significant fraction of the targets, especially if many have orbits measured in decades. This result supports the idea that the formation of close binaries involves the deposition of angular momentum into the orbital motion of a third star.

We summarize some of the compelling new scientific opportunities for understanding stars and stellar systems that can be enabled by sub-mas angular resolution, UV/Optical spectral imaging observations, which can reveal the details of the many dynamic processes (e.g., variable magnetic fields, accretion, convection, shocks, pulsations, winds, and jets) that affect their 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. Two mission concepts which could provide these invaluable observations are NASA's Stellar Imager (SI; http://hires.gsfc.nasa.gov/si/) interferometer and ESA's Luciola sparse aperture hypertelescope, which each could resolve hundreds of stars and stellar systems. These observatories will also open an immense new discovery space for astrophysical research in general and, in particular, for Active Galactic Nuclei (Kraemer et al. Decadal Survey Science Whitepaper). The technology developments needed for these missions are challenging, but eminently feasible (Carpenter et al. Decadal Survey Technology Whitepaper) with a reasonable investment over the next decade to enable flight in the 2025+ timeframe. That investment 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.

We report the results of our study of magnetic fields in a sample of Be stars using spectropolarimetric data obtained at the European Southern Observatory with the multi-mode instrument FORS1 installed at the 8m Kueyen telescope. The detected magnetic fields are rather weak, not stronger than ~150G. A few classical Be stars display cyclic variability of the magnetic field with periods of tens of minutes.

The rapid rotation of Be stars may be caused in some cases by past mass and angular momentum accretion in an interacting binary in which the mass donor is currently viewed as a small, hot subdwarf stripped of its outer envelope. Here we report on the spectroscopic detection of such a subdwarf in the Be binary system FY Canis Majoris from the analysis of data acquired by the IUE spacecraft and KPNO Coude Feed Telescope over the course of 16 and 21 years, respectively. We present a double-lined spectroscopic orbit for the binary based upon radial velocities from the IUE spectra and use the orbital solutions with a Doppler tomography algorithm to reconstruct the components' UV spectra. The subdwarf is hot (T_eff = 45+/-5 kK) and has a mass of about 1.3 M_sun and a radius of about 0.6 R_sun. It contributes about 4% as much flux as the Be star does in the FUV. We also present observations of the H-alpha and He I 6678 emission features that are formed in the circumstellar disk of the Be star. Orbital flux and velocity variations in the He I 6678 profile indicate that much of the emission forms along the disk rim facing the hot subdwarf where the disk is probably heated by the incident radiation from the subdwarf. A study of the FUV infall shell lines discovered in the 1980s confirms their episodic presence but reveals that they tend to be found around both quadrature phases, unlike the pattern in Algol binaries. Phase-dependent variations in the UV N V doublet suggest the presence of a N-enhanced wind from the subdwarf and a possible shock-interaction region between the stars where the subdwarf's wind collides with the disk of the Be star.

We review the recent discoveries of magnetic fields in different types of massive stars and briefly discuss strategies for spectropolarimetric observations to be carried out in the future.

We report the results of our study of magnetic fields in a sample of 15 Be stars using spectropolarimetric data obtained at the European Southern Observatory with the multi-mode instrument FORS1 installed at the 8m Kueyen telescope. We detect weak photospheric magnetic fields in four stars, HD56014, HD148184, HD155806, and HD181615. We note that for HD181615 the evolutionary status is not obvious due to the fact that it is a binary system currently observed in the initial rapid phase of mass exchange between the two components. Further, we notify the possible presence of distinct circular polarisation features in the circumstellar components of Ca II H&K in three stars, HD58011, HD117357, and HD181615, hinting at a probable presence of magnetic fields in the circumstellar mass loss disks of these stars. We emphasize the need for future spectropolarimetric observations of Be stars with detected magnetic fields to study the temporal evolution of their magnetic fields and the correlation of magnetic field properties with dynamical phenomena taking place in the gaseous circumstellar disks of these stars.

We have developed a systematic procedure to study the disks in Algol-type binaries using spectroscopic analysis, synthetic spectra, and tomography. We analyzed 119 H-alpha spectra of TT Hya, an Algol-type eclipsing interacting binary, collected from 1985-2001. The new radial velocities enabled us to derive reliable orbital elements, including a small non-zero eccentricity, and to improve the accuracy of the absolute dimensions of the system. High resolution IUE spectra were also analyzed to study the formation of the ultraviolet lines and continuum. Synthetic spectra of the iron curtain using our new shellspec program enabled us to derive a characteristic disk temperature of 7000K. We have demonstrated that the UV emission lines seen during total primary eclipse cannot originate from the accretion disk, but most likely arise from a hotter disk-stream interaction region. The synthetic spectra of the stars, disk, and stream allowed us to derive a lower limit to the mass transfer rate of 2e-10 solar masses per year. Doppler tomography of the observed H-alpha profiles revealed a distinct accretion disk. The difference spectra produced by subtracting the synthetic spectra of the stars resulted in an image of the disk, which virtually disappeared once the composite synthetic spectra of the stars and disk were used to calculate the difference spectra. An intensity enhancement of the resulting tomogram revealed images of the gas stream and an emission arc. We successfully modeled the gas stream using shellspec and associated the emission arc with an asymmetry in the accretion disk.

We present an investigation of H-alpha emission line variations observed in the massive Algol binary, RY Per. We give new radial velocity data for the secondary based upon our optical spectra and for the primary based upon high dispersion UV spectra. We present revised orbital elements and an estimate of the primary's projected rotational velocity (which indicates that the primary is rotating 7 times faster than synchronous). We use a Doppler tomography algorithm to reconstruct the individual primary and secondary spectra in the region of H-alpha, and we subtract the latter from each of our observations to obtain profiles of the primary and its disk alone. Our H-alpha observations of RY Per show that the mass gaining primary is surrounded by a persistent but time variable accretion disk. The profile that is observed outside-of-eclipse has weak, double-peaked emission flanking a deep central absorption, and we find that these properties can be reproduced by a disk model that includes the absorption of photospheric light by the band of the disk seen in projection against the face of the star. We developed a new method to reconstruct the disk surface density distribution from the ensemble of H-alpha profiles observed around the orbit, and this method accounts for the effects of disk occultation by the stellar components, the obscuration of the primary by the disk, and flux contributions from optically thick disk elements. The resulting surface density distribution is elongated along the axis joining the stars, in the same way as seen in hydrodynamical simulations of gas flows that strike the mass gainer near trailing edge of the star. This type of gas stream configuration is optimal for the transfer of angular momentum, and we show that rapid rotation is found in other Algols that have passed through a similar stage.