Blue straggler stars (BSS) have been identified in star clusters and in field populations in our own Milky Way galaxy and in its satellite galaxies. They manifest as stars bluer and more luminous than the dominant old population, and usually have a spatial distribution that follows the old population. Their progenitors are likely to have been close binaries. We investigate trends of the BSS population in dwarf spheroidal galaxies (dSph) and in the bulge of the Milky Way and find an anti-correlation between the relative frequency of BSS and the metallicity of the parent population. The rate of occurrence of BSS in the metal-poor dwarf galaxies is approximately twice that found in the solar-metallicity bulge population. This trend of decreasing relative population of BSS with increasing metallicity mirrors that found for the close-binary fraction in the field population of the Milky Way. We argue that the dominant mode of BSS formation in low-density environments is likely to be mass transfer in close-binary systems. It then follows that the similarity between the trends for BSS in the dSph and field stars in our Galaxy supports the proposal that the small-scale fragmentation during star formation is driven by the same dominant physical process, despite the diversity in environments, plausibly gravitational instability of proto-stellar discs.

MCA-1B (also called UIT003) is a luminous hot star in the western outskirts of M33, classified over 20yr ago with a spectral type of Ofpe/WN9 and identified then as a candidate luminous blue variable (LBV). Palomar Transient Factory data reveal that this star brightened in 2010, with a light curve resembling that of the classic LBV star AFAnd in M31. Other Ofpe/WN9 stars have erupted as LBVs, but MCA-1B was unusual because it remained hot. It showed a WN-type spectrum throughout its eruption, whereas LBVs usually get much cooler. MCA-1B showed an almost four-fold increase in bolometric luminosity and a doubling of its radius, but its temperature stayed around 29kK. As it faded, it shifted to even hotter temperatures, exhibiting a WN7/WN8-type spectrum, and doubling its wind speed. MCA-1B is reminiscent of some supernova impostors, and its location resembles the isolated environment of SN 2009ip. It is most similar to HD5980 (in the SMC) and GR 290 (also in M33). Whereas these two LBVs exhibited B-type spectra in eruption, MCA-1B is the first clear case where a Wolf-Rayet (WR) spectrum persisted at all times. Together, MCA-1B, HD 5980, and GR 290 constitute a class of WN-type LBVs, distinct from S Doradus LBVs. They are most interesting in the context of LBVs at low metallicity, a possible post-LBV/WR transition in binaries, and as likely Type~Ibn supernova progenitors.

Close binaries suppress the formation of circumstellar (S-type) planets and therefore significantly bias the inferred planet occurrence rates and statistical trends. After compiling various radial velocity and high-resolution imaging surveys, we determine that binaries with a < 1 au fully suppress S-type planets, binaries with a = 10 au host close planets at 15% the occurrence rate of single stars, and wide binaries with a > 200 au have a negligible effect on planet formation. We show that F = 43% +/- 6% of solar-type primaries do not host close planets due to suppression by close stellar companions. By removing spectroscopic binaries from their samples, radial velocity surveys for giant planets boost their detection rates by a factor of 1/(1-F) = 1.8 +/- 0.2 compared to transiting surveys. This selection bias fully accounts for the discrepancy in hot Jupiter occurrence rates inferred from these two detection methods. Correcting for both planet suppression by close binaries and transit dilution by wide binaries, the occurrence rate of small planets orbiting single G-dwarfs is 2.1 +/- 0.3 times larger than the rate inferred from all G-dwarfs in the Kepler survey. About half (but not all) of the observed increase in small, short-period planets toward low-mass hosts can be explained by the corresponding decrease in the binary fraction. Finally, we demonstrate that the apparent enhancement of wide stellar companions to hot Jupiter hosts is due to multiple selection effects. Although very close binaries with secondary masses M_2 > 10 M_J preferentially have wide tertiary companions, genuine hot Jupiters with M_p = 0.2-4 M_J that formed via core accretion do not exhibit a statistically significant excess of wide stellar companions.

Joint statistics of periods and mass ratios of close binaries and its dependence on primary mass can be explained by assuming that seed binary companions are formed by disc fragmentation at random intervals during assemblage of stellar mass and migrate inwards as they accrete from the circumbinary disk. A toy model based on simple prescriptions for the companion growth and migration reproduces such aspects of close solar-mass binaries as the distribution of binary periods P, the brown dwarf desert at short P, the nearly uniform distribution of mass ratios, and a population of equal-mass binaries (twins) that decreases linearly in frequency with logP. For massive stars, the model predicts a large fraction of early mergers, a distribution of logP with a negative slope, and a mass-ratio distribution that is also uniform but with a substantially reduced twin fraction. By treating disc fragmentation as a stochastic process, we also reproduce the observed properties of compact triples. Success of our toy model suggests that most close binaries and compact triples indeed formed by disc fragmentation followed by accretion-driven inward migration.

We use a homogeneous catalog of 42,000 main-sequence wide binaries identified by Gaia to measure the mass ratio distribution, p(q), of binaries with primary masses $0.1<M_1/M_{\odot}<2.5$, mass ratios $0.1 \lesssim q<1$, and separations $50<s/{\rm AU}<50,000$. A well-understood selection function allows us to constrain p(q) in 35 independent bins of primary mass and separation, with hundreds to thousands of binaries in each bin. Our investigation reveals a sharp excess of equal-mass "twin" binaries that is statistically significant out to separations of 1,000 to 10,000 AU, depending on primary mass. The excess is narrow: a steep increase in p(q) at $0.95 \lesssim q<1$, with no significant excess at $q\lesssim 0.95$. A range of tests confirm the signal is real, not a data artifact or selection effect. Combining the Gaia constraints with those from close binaries, we show that the twin excess decreases with increasing separation, but its width ($q\gtrsim 0.95$) is constant over $0.01<a/{\rm AU}<10,000$. The wide twin population would be difficult to explain if the components of all wide binaries formed via core fragmentation, which is not expected to produce strongly correlated component masses. We conjecture that wide twins formed at closer separations ($a \lesssim 100$ AU), likely via accretion from circumbinary disks, and were subsequently widened by dynamical interactions in their birth environments. The separation-dependence of the twin excess then constrains the efficiency of dynamical widening and disruption of binaries in young clusters. We also constrain p(q) across $0.1 \lesssim q<1$. Besides changes in the twin fraction, p(q) is independent of separation at fixed primary mass over $100 \lesssim s/{\rm AU} < 50,000$. It is flatter than expected for random pairings from the IMF but more bottom-heavy for wide binaries than for binaries with $a\lesssim$100 AU.

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches.

One clear observational prediction of the single degenerate progenitor scenario as the origin of type Ia supernovae (SNe) is the presence of relatively narrow ($\approx$1000 km s$^{-1}$) H$\alpha$ emission at nebular phases, although this feature is rarely seen. We present a compilation of nebular phase H$\alpha$ limits for SN Ia in the literature and demonstrate that this heterogenous sample has been biased towards SN Ia with relatively high luminosities and slow decline rates, as parameterized by $\Delta$m$_{15}(B)$, the difference in $B$-band magnitude between maximum light and fifteen days afterward. Motivated by the need to explore the full parameter space of SN~Ia and their subtypes, we present two new and six previously published nebular spectra of SN Ia with $\Delta$m$_{15}(B)$$ > $1.3 mag (including members of the transitional and SN1991bg-like subclasses) and measure nondetection limits of $L_{H\alpha}$$~<~$0.85--9.9$\times$10$^{36}$ ergs s$^{-1}$, which we confirmed by implanting simulated H$\alpha$ emission into our data. Based on the lastest models of swept-up material stripped from a nondegenerate companion star, these $L_{H\alpha}$ values correspond to hydrogen mass limits of $M_H$$~\lesssim~$1-3$\times$10$^{-4}$ $M_{\odot}$, roughly three orders of magnitude below that expected for the systems modeled, although we note that no simulations of H$\alpha$ nebular emission in such weak explosions have yet been performed. Despite the recent detection of strong H$\alpha$ in ASASSN-18tb (SN 2018fhw; $\Delta$m$_{15}(B)$ = 2.0 mag), we see no evidence that fast declining systems are more likely to have late time H$\alpha$ emission, although a larger sample is needed to confirm this result.

Multiplicity is a fundamental property that is set early during stellar lifetimes, and it is a stringent probe of the physics of star formation. The distribution of close companions around young stars is still poorly constrained by observations. We present an analysis of stellar multiplicity derived from APOGEE-2 spectra obtained in targeted observations of nearby star-forming regions. This is the largest homogeneously observed sample of high-resolution spectra of young stars. We developed an autonomous method to identify double lined spectroscopic binaries (SB2s). Out of 5007 sources spanning the mass range of $\sim$0.05--1.5 \msun, we find 399 binaries, including both RV variables and SB2s. The mass ratio distribution of SB2s is consistent with a uniform for $q<0.95$ with an excess of twins with $q>0.95$. The period distribution is consistent with what has been observed in close binaries ($<10$ AU) in the evolved populations. Three systems are found to have $q\sim$0.1, with a companion located within the brown dwarf desert. There are not any strong trends in the multiplicity fraction (MF) as a function of cluster age from 1 to 100 Myr. There is a weak dependence on stellar density, with companions being most numerous at $\Sigma_*\sim30$ stars/pc$^{-2}$, and decreasing in more diffuse regions. Finally, disk-bearing sources are deficient in SB2s (but not RV variables) by a factor of $\sim$2; this deficit is recovered by the systems without disks. This may indicate a quick dispersal of disk material in short-period equal mass systems that is less effective in binaries with lower $q$.

Our uncertainties about binary star systems (and triples and so on) limit our capabilities in literally every single one of the Thematic Areas identified for Astro2020. We need to understand the population statistics of stellar multiplicity and their variations with stellar type, chemistry, and dynamical environment: Correct interpretation of any exoplanet experiment depends on proper treatment of resolved and unresolved binaries; stellar multiplicity is a direct outcome of star and companion formation; the most precise constraints on stellar structure come from well-characterized binary systems; stellar populations heavily rely on stellar and binary evolution modeling; high-redshift galaxy radiation and reionization is controlled by binary-dependent stellar physics; compact objects are the outcomes of binary evolution; the interpretation of multi-messenger astronomy from gravitational waves, light, and neutrinos relies on understanding the products of binary star evolution; near-Universe constraints on the Hubble constant with Type Ia supernovae and gravitational-wave mergers are subject to systematics related to their binary star progenitors; local measures of dark-matter substructure masses are distorted by binary populations. In order to realize the scientific goals in each of these themes over the next decade, we therefore need to understand how binary stars and stellar multiplets are formed and distributed in the space of masses, composition, age, and orbital properties, and how the distribution evolves with time. This white paper emphasizes the interdisciplinary importance of binary-star science and advocates that coordinated investment from all astrophysical communities will benefit almost all branches of astrophysics.

The distributions of the initial main-sequence binary parameters are one of the key ingredients in obtaining evolutionary predictions for compact binary (BH-BH / BH-NS / NS-NS) merger rates. Until now, such calculations were done under the assumption that initial binary parameter distributions were independent. Here, we implement empirically derived inter-correlated distributions of initial binary parameters primary mass (M1), mass ratio (q), orbital period (P), and eccentricity (e). Unexpectedly, the introduction of inter-correlated initial binary parameters leads to only a small decrease in the predicted merger rates by a factor of 2 $-$ 3 relative to the previously used non-correlated initial distributions. The formation of compact object mergers in the isolated classical binary evolution favors initial binaries with stars of comparable masses (q = 0.5 $-$ 1) at intermediate orbital periods (log P (days) = 2 $-$ 4). New distributions slightly shift the mass ratios towards smaller values with respect to the previously used flat q distribution, which is the dominant effect decreasing the rates. New orbital periods only negligibly increase the number of progenitors. Additionally, we discuss the uncertainty of merger rate predictions associated with possible variations of the massive-star initial mass function (IMF). We argue that evolutionary calculations should be normalized to a star formation rate (SFR) that is obtained from the observed amount of UV light at wavelength 1500{\AA} (SFR indicator). In this case, contrary to recent reports, the uncertainty of the IMF does not affect the rates by more than a factor of 2. Any change to the IMF slope for massive stars requires a change of SFR in a way that counteracts the impact of IMF variations on the merger rates. In contrast, we suggest that the uncertainty in cosmic SFR at low metallicity can be a significant factor at play.

We use GAIA DR2 proper motions of the RIOTS4 field OB stars in the Small Magellanic Cloud (SMC) to study the kinematics of runaway stars. The data reveal that the SMC Wing has a systemic peculiar motion relative to the SMC Bar of (v_RA, v_Dec) = (62 +/-7, -18+/-5) km/s and relative radial velocity +4.5 +/- 5.0 km/s. This unambiguously demonstrates that these two regions are kinematically distinct: the Wing is moving away from the Bar, and towards the Large Magellanic Cloud with a 3-D velocity of 64 +/- 10 km/s. This is consistent with models for a recent, direct collision between the Clouds. We present transverse velocity distributions for our field OB stars, confirming that unbound runaways comprise on the order of half our sample, possibly more. Using eclipsing binaries and double-lined spectroscopic binaries as tracers of dynamically ejected runaways, and high-mass X-ray binaries (HMXBs) as tracers of runaways accelerated by supernova kicks, we find significant contributions from both populations. The data suggest that HMXBs have lower velocity dispersion relative to dynamically ejected binaries, consistent with the former corresponding to less energetic supernova kicks that failed to unbind the components. Evidence suggests that our fast runaways are dominated by dynamical, rather than supernova, ejections.

There is now strong evidence that the close binary fraction (P < 10$^4$ days; a < 10 AU) of solar-type stars ($M_1$ = 0.6-1.5M$_{\odot}$) decreases significantly with metallicity. Although early surveys showed that the observed spectroscopic binary (SB) fractions in the galactic disk and halo are similar (e.g., Carney-Latham sample), these studies did not correct for incompleteness. In this study, we examine five different surveys and thoroughly account for their underlying selection biases to measure the intrinsic occurrence rate of close solar-type binaries. We re-analyze: (1) a volume-limited sample of solar-type stars, (2) an SB survey of high-proper-motion stars, (3) various SB samples of metal-poor giants, (4) the APOGEE survey of radial velocity (RV) variables, and (5) Kepler eclipsing binaries (EBs). The observed APOGEE RV variability fraction and Kepler EB fraction both decrease by a factor of $\approx$4 across $-$1.0 < [Fe/H] < 0.5 at the 22$\sigma$ and 9$\sigma$ confidence levels, respectively. After correcting for incompleteness, all five samples exhibit a quantitatively consistent anti-correlation between the intrinsic close binary fraction (a < 10 AU) and metallicity: $F_{\rm close}$ = 53%$\pm$12%, 40%$\pm$6%, 24%$\pm$4%, and 10%$\pm$3% at [Fe/H] = $-$3.0, $-$1.0, $-$0.2 (mean field metallicity), and +0.5, respectively. We present fragmentation models that explain why the close binary fraction of solar-type stars strongly decreases with metallicity while the wide binary fraction, close binary fraction of OB stars, and initial mass function are all constant across $-$1.5 < [Fe/H] < 0.5. The majority of solar-type stars with [Fe/H] < $-$1.0 will interact with a stellar companion, which has profound implications for binary evolution in old and metal-poor environments such as the galactic halo, bulge, thick disk, globular clusters, dwarf galaxies, and high-redshift universe.

High-resolution spectroscopy and speckle interferometry reveal the young star HD 86588 as a quadruple system with a 3-tier hierarchy. The 0.3" resolved binary A,B with an estimated period around 300 years contains the 8-year pair Aa,Abc (also potentially resolvable), where Ab,Ac is a double-lined binary with equal components, for which we compute the spectroscopic orbit. Despite the short period of 2.4058 day, the orbit of Ab,Ac is eccentric (e=0.086+-0.003). It has a large inclination, but there are no eclipses; only a 4.4 mmag light modulation apparently caused by star spots on the components of this binary is detected with Evryscope. Assuming a moderate extinction of A_V = 0.5 mag and a parallax of 5.2 mas, we find that the stars are on or close to the main sequence (age >10 Myr) and their masses are from 1 to 1.3 solar. We measure the strength of the Lithium line in the visual secondary B which, together with rotation, suggests that the system is younger than 150 Myr. This object is located behind the extension of the Chamaeleon I dark cloud (which explains extinction and interstellar Sodium absorption), but apparently does not belong to it. We propose a scenario where the inner orbit has recently acquired its high eccentricity through dynamical interaction with the outer two components; it is now undergoing rapid tidal circularization on a time scale of ~1 Myr. Alternatively, the eccentricity could be excited quasi-stationary by the outer component Aa.

Solar-type binaries with short orbital periods ($P_{\rm close}$ $\equiv$ 1 - 10 days; $a$ $\lesssim$ 0.1 AU) cannot form directly via fragmentation of molecular clouds or protostellar disks, yet their component masses are highly correlated, suggesting interaction during the pre-main-sequence (pre-MS) phase. Moreover, the close binary fraction of pre-MS stars is consistent with that of their MS counterparts in the field ($F_{\rm close}$ = 2.1%). Thus we can infer that some migration mechanism operates during the early pre-MS phase ($\tau$ $\lesssim$ 5 Myr) that reshapes the primordial separation distribution. We test the feasibility of this hypothesis by carrying out a population synthesis calculation which accounts for two formation channels: Kozai-Lidov (KL) oscillations and dynamical instability in triple systems. Our models incorporate (1) more realistic initial conditions compared to previous studies, (2) octupole-level effects in the secular evolution, (3) tidal energy dissipation via weak-friction equilibrium tides at small eccentricities and via non-radial dynamical oscillations at large eccentricities, and (4) the larger tidal radius of a pre-MS primary. Given a 15% triple star fraction, we simulate a close binary fraction from KL oscillations alone of $F_{\rm close}$ $\approx$ 0.4% after $\tau$ = 5 Myr, which increases to $F_{\rm close}$ $\approx$ 0.8% by $\tau$ = 5 Gyr. Dynamical ejections and disruptions of unstable coplanar triples in the disk produce solitary binaries with slightly longer periods $P$ $\approx$ 10 - 100 days. The remaining $\approx$60% of close binaries with outer tertiaries, particularly those in compact coplanar configurations with log $P_{\rm out}$ (days) $\approx$ 2 - 5 ($a_{\rm out}$ $<$ 50 AU), can be explained only with substantial extra energy dissipation due to interactions with primordial gas.

We use the multi-epoch radial velocities acquired by the APOGEE survey to perform a large scale statistical study of stellar multiplicity for field stars in the Milky Way, spanning the evolutionary phases between the main sequence and the red clump. We show that the distribution of maximum radial velocity shifts (\drvm) for APOGEE targets is a strong function of \logg, with main sequence stars showing \drvm\ as high as $\sim$300 \kms, and steadily dropping down to $\sim$30 \kms\ for \logg$\sim$0, as stars climb up the Red Giant Branch (RGB). Red clump stars show a distribution of \drvm\ values comparable to that of stars at the tip of the RGB, implying they have similar multiplicity characteristics. The observed attrition of high \drvm\ systems in the RGB is consistent with a lognormal period distribution in the main sequence and a multiplicity fraction of 0.35, which is truncated at an increasing period as stars become physically larger and undergo mass transfer after Roche Lobe Overflow during H shell burning. The \drvm\ distributions also show that the multiplicity characteristics of field stars are metallicity dependent, with metal-poor ([Fe/H]$\lesssim-0.5$) stars having a multiplicity fraction a factor 2-3 higher than metal-rich ([Fe/H]$\gtrsim0.0$) stars. This has profound implications for the formation rates of interacting binaries observed by astronomical transient surveys and gravitational wave detectors, as well as the habitability of circumbinary planets.

The orbital parameters of binaries at intermediate periods ($10^2$--$10^3$ d) are difficult to measure with conventional methods and are very incomplete. We have undertaken a new survey, applying our pulsation timing method to {\it Kepler} light curves of 2224 main-sequence A/F stars and found 341 non-eclipsing binaries. We calculate the orbital parameters for 317 PB1 systems (single-pulsator binaries) and 24 PB2s (double-pulsators). The method reaches down to small mass ratios $q=0.02$ and yields a highly homogeneous sample. We parametrize the mass-ratio distribution using both inversion and MCMC forward-modelling techniques, and find it to be skewed towards low-mass companions, peaking at $q=0.2$. While solar-type primaries show a brown dwarf desert across intermediate periods, we find a small but significant (2.6$\sigma$) population of extreme-mass-ratio companions ($q<0.1$) to our intermediate-mass primaries. Across periods of 100--1500 d and at $q>0.1$, we measure the binary fraction of current A/F primaries to be 15.4%$\pm$1.4%, though we find that a large fraction of the companions (21%$\pm$6%) are white dwarfs in post-mass-transfer systems with primaries that are now blue stragglers, some of which are the progenitors of Type Ia supernovae, barium stars, symbiotics, and related phenomena. Excluding these white dwarfs, we find the binary fraction of original A/F primaries to be 13.9%$\pm$2.1% over the same parameter space. Combining literature measurements with ours, we find the binary fraction across these periods is a constant 5% for primaries $M_1<0.8$ M$_{\odot}$, but then increases linearly with $\log M_1$, demonstrating that natal discs around more massive protostars $M_1>1$ M$_{\odot}$ become increasingly more prone to fragmentation. Finally, we find the eccentricity distribution of the main-sequence pairs to be much less eccentric than the thermal distribution.

We compile observations of early-type binaries identified via spectroscopy, eclipses, long-baseline interferometry, adaptive optics, common proper motion, etc. Each observational technique is sensitive to companions across a narrow parameter space of orbital periods P and mass ratios q = M_comp/M_1. After combining the samples from the various surveys and correcting for their respective selection effects, we find the properties of companions to O-type and B-type main-sequence (MS) stars differ among three regimes. First, at short orbital periods P < 20 days (separations a < 0.4 AU), the binaries have small eccentricities e < 0.4, favor modest mass ratios <q> = 0.5, and exhibit a small excess of twins q > 0.95. Second, the companion frequency peaks at intermediate periods log P (days) = 3.5 (a = 10 AU), where the binaries have mass ratios weighted toward small values q = 0.2-0.3 and follow a Maxwellian "thermal" eccentricity distribution. Finally, companions with long orbital periods log P (days) = 5.5-7.5 (a = 200-5,000 AU) are outer tertiary components in hierarchical triples, and have a mass ratio distribution across q = 0.1-1.0 that is nearly consistent with random pairings drawn from the initial mass function. We discuss these companion distributions and properties in the context of binary star formation and evolution. We also reanalyze the binary statistics of solar-type MS primaries, taking into account that (30+/-10)% of single-lined spectroscopic binaries likely contain white dwarf companions instead of low-mass stellar secondaries. The mean frequency of stellar companions with q > 0.1 and log P (days) < 8.0 per primary increases from 0.50+/-0.04 for solar-type MS primaries to 2.1+/-0.3 for O-type MS primaries. We fit joint probability density functions f(M_1,q,P,e) to the corrected distributions, which can be incorporated into binary population synthesis studies.

We conduct a spectroscopic search of quasars observed by the Sloan Digital Sky Survey (SDSS) with broad absorption line (BAL) troughs due to Mg II and troughs due to Fe II that simultaneously exhibit strong Balmer narrow emission lines (NELs). We find that in a redshift range of 0.4 less than or equal to z less than or equal to 0.9 approximately 23 of the 70 Mg II BALs and 4 of a subset of 15 Fe II BALs exhibit strong Balmer emission. We also find significant fractions of Mg II BALs (approximately 23%) and those Mg II BALs with Fe II troughs (approximately 27%) have strong continuum reddening, E(B - V) greater than or equal to 0.1. From measurements of the Balmer decrement in three objects, we find similarly significant reddening of the NEL region in three of the four objects; the NELs in the fourth object are not measurable. We also include one object in this study not taken from the SDSS sample that shows Fe II absorption and strong narrow emission, but due to measurement uncertainty and low continuum reddening the comparison is consistent but inconclusive. We find a trend in both the Mg II and Fe II BAL samples between the NEL reddening and continuum reddening. Because the narrow line reddening is consistent with the continuum reddening in every object in the two SDSS samples, it suggests that the reddening sources in these objects likely exist at larger radial distances than the narrow line regions from the central nucleus.

We analyze 221 eclipsing binaries (EBs) in the Large Magellanic Cloud with B-type main-sequence (MS) primaries ($M_1$ $\approx$ 4 - 14 M$_{\odot}$) and orbital periods $P$ = 20 - 50 days that were photometrically monitored by the Optical Gravitational Lensing Experiment. We utilize our three-stage automated pipeline to (1) classify all 221 EBs, (2) fit physical models to the light curves of 130 detached well-defined EBs from which unique parameters can be determined, and (3) recover the intrinsic binary statistics by correcting for selection effects. We uncover two statistically significant trends with age. First, younger EBs tend to reside in dustier environments with larger photometric extinctions, an empirical relation that can be implemented when modeling stellar populations. Second, younger EBs generally have large eccentricities. This demonstrates that massive binaries at moderate orbital periods are born with a Maxwellian "thermal" orbital velocity distribution, which indicates they formed via dynamical interactions. In addition, the age-eccentricity anticorrelation provides a direct constraint for tidal evolution in highly eccentric binaries containing hot MS stars with radiative envelopes. The intrinsic fraction of B-type MS stars with stellar companions $q$ = $M_2$/$M_1$ $>$ 0.2 and orbital periods $P$ = 20 - 50 days is (7 $\pm$ 2)%. We find early-type binaries at $P$ = 20 - 50 days are weighted significantly toward small mass ratios $q$ $\approx$ 0.2 - 0.3, which is different than the results from previous observations of closer binaries with $P$ $<$ 20 days. This indicates that early-type binaries at slightly wider orbital separations have experienced substantially less coevolution and competitive accretion during their formation in the circumbinary disk.

Early B-type main-sequence (MS) stars (M$_1$ = 5-16 M$_{\odot}$) with closely orbiting low-mass stellar companions (q = M$_2$/M$_1$ < 0.25) can evolve to produce Type Ia supernovae, low-mass X-ray binaries, and millisecond pulsars. However, the formation mechanism and intrinsic frequency of such close extreme mass-ratio binaries have been debated, especially considering none have hitherto been detected. Utilizing observations of the Large Magellanic Cloud galaxy conducted by the Optical Gravitational Lensing Experiment, we have discovered a new class of eclipsing binaries in which a luminous B-type MS star irradiates a closely orbiting low-mass pre-MS companion that has not yet fully formed. The primordial pre-MS companions have large radii and discernibly reflect much of the light they intercept from the B-type MS primaries ($\Delta$I$_{\rm refl}$ = 0.02-0.14 mag). For the 18 definitive MS + pre-MS eclipsing binaries in our sample with good model fits to the observed light curves, we measure short orbital periods P = 3.0-8.5 days, young ages $\tau$ = 0.6-8 Myr, and small secondary masses M$_2$ = 0.8-2.4 M$_{\odot}$ (q = 0.07-0.36). The majority of these nascent eclipsing binaries are still associated with stellar nurseries, e.g. the system with the deepest eclipse $\Delta$I$_1$ = 2.8 mag and youngest age $\tau$ = 0.6$\pm$0.4 Myr is embedded in the bright H II region 30 Doradus. After correcting for selection effects, we find that (2.0$\pm$0.6)% of B-type MS stars have companions with short orbital periods P = 3.0-8.5 days and extreme mass ratios q = 0.06-0.25. This is $\approx$10 times greater than that observed for solar-type MS primaries. We discuss how these new eclipsing binaries provide invaluable insights, diagnostics, and challenges for the formation and evolution of stars, binaries, and H II regions.

In order to understand the rates and properties of Type Ia and Type Ib/c supernovae, X-ray binaries, gravitational wave sources, and gamma ray bursts as a function of galactic environment and cosmic age, it is imperative that we measure how the close binary properties of O and B-type stars vary with metallicity. We have studied eclipsing binaries with early-B main-sequence primaries in three galaxies with different metallicities: the Large and Small Magellanic Clouds (LMC and SMC, respectively) as well as the Milky Way (MW). The observed fractions of early-B stars which exhibit deep eclipses 0.25 < Delta(m) (mag) < 0.65 and orbital periods 2 < P (days) < 20 in the MW, LMC, and SMC span a narrow range of (0.7-1.0)%, which is a model independent result. After correcting for geometrical selection effects and incompleteness toward low-mass companions, we find for early-B stars in all three environments: (1) a close binary fraction of (22+/-5)% across orbital periods 2 < P (days) < 20 and mass ratios q = M_2/M_1 > 0.1, (2) an intrinsic orbital period distribution slightly skewed toward shorter periods relative to a distribution that is uniform in log P, (3) a mass-ratio distribution weighted toward low-mass companions, and (4) a small, nearly negligible excess fraction of twins with q > 0.9. Our fitted parameters derived for the MW eclipsing binaries match the properties inferred from nearby, early-type spectroscopic binaries, which further validates our results. There are no statistically significant trends with metallicity, demonstrating that the close binary properties of massive stars do not vary across metallicities -0.7 < log(Z/Z_sun) < 0.0 beyond the measured uncertainties.

We report on our discovery and observations of the Pan-STARRS1 supernova (SN) PS1-12sk, a transient with properties that indicate atypical star formation in its host galaxy cluster or pose a challenge to popular progenitor system models for this class of explosion. The optical spectra of PS1-12sk classify it as a Type Ibn SN (c.f. SN 2006jc), dominated by intermediate-width (3x10^3 km/s) and time variable He I emission. Our multi-wavelength monitoring establishes the rise time dt = 9-23 days and shows an NUV-NIR SED with temperature > 17x10^3 K and a peak rise magnitude of Mz = -18.9 mag. SN Ibn spectroscopic properties are commonly interpreted as the signature of a massive star (17 - 100 M_sun) explosion within a He-enriched circumstellar medium. However, unlike previous Type Ibn supernovae, PS1-12sk is associated with an elliptical brightest cluster galaxy, CGCG 208-042 (z = 0.054) in cluster RXC J0844.9+4258. The expected probability of an event like PS1-12sk in such environments is low given the measured infrequency of core-collapse SNe in red sequence galaxies compounded by the low volumetric rate of SN Ibn. Furthermore, we find no evidence of star formation at the explosion site to sensitive limits (Sigma Halpha < 2x10^-3 M_sun/yr/kpc^2). We therefore discuss white dwarf binary systems as a possible progenitor channel for SNe Ibn. We conclude that PS1-12sk represents either a fortuitous and statistically unlikely discovery, evidence for a top-heavy IMF in galaxy cluster cooling flow filaments, or the first clue suggesting an alternate progenitor channel for Type Ibn SNe.

In an attempt to determine how many planetary nebulae derive from binary interactions, we have started a project to measure their unbiased binary fraction. This number, when compared to the binary fraction of the presumed parent population can give a first handle on the origin of planetary nebulae. By detecting 27 bona fide central stars in the I band we have found that 30% of our sample have an I band excess between one and a few sigmas, possibly denoting companions brighter than M3-4V and with separations smaller than approximately 1000 AU. By accounting for the undetectable companions, we determine a de-biased binary fraction of 67-78% for all companions at all separations. We compare this number to a main sequence binary fraction of (50+/-4)% determined for spectral types F6V-G2V, appropriate if the progenitors of today's PN central star population is indeed the F6V-G2V stars. The error on our estimate could be between 10 and 30%. We conclude that the central star binary fraction may be larger than expected from the putative parent population. Using the more sensitive J band of a subset of 11 central stars, the binary fraction is 54% for companions brighter than approximately M5-6V and with separations smaller than about 900 AU. De-biassing this number we obtain a binary fraction of 100-107%. The two numbers should be the same and the discrepancy is likely due to small number statistics. We also present an accurately vetted compilation of observed main sequence star magnitudes, colours and masses, which can serve as a reference for future studies. We also present synthetic colours of hot stars as a function of temperature (20-170kK) and gravity (log g= 6-8) for Solar and PG1159 compositions.

We present photometric and spectroscopic observations of a peculiar variable (designated DASCH J075731.1+201735 or J0757) discovered from our DASCH project using the digitized Harvard College Observatory archival photographic plates. It brightened by about 1.5 magnitudes in B within a year starting in 1942, and then slowly faded back to its pre-outburst brightness from 1943 to the 1950s. The mean brightness level was stable before and after the outburst, and ellipsoidal variations with a period of $P=119.18\pm0.07$ days are seen, suggesting that the star is tidally distorted. Radial-velocity measurements indicate that the orbit is nearly circular ($e=0.02\pm0.01$) with a spectroscopic period that is the same as the photometric period. The binary consists of a $1.1\pm0.3 M_\odot$ M0III star, and a $0.6\pm0.2 M_\odot$ companion, very likely a white dwarf (WD). Unlike other symbiotic binaries, there is no sign of emission lines or a stellar wind in the spectra. With an outburst timescale of ~10 years and estimated B band peak luminosity M_B~0.7, J0757 is different from any other known classic or symbiotic novae. The most probable explanation of the outburst is Hydrogen shell-burning on the WD, although an accretion-powered flare cannot be ruled out.

We report unique EVLA observations of SN 2011fe representing the most sensitive radio study of a Type Ia supernova to date. Our data place direct constraints on the density of the surrounding medium at radii ~10^15-10^16 cm, implying an upper limit on the mass loss rate from the progenitor system of Mdot <~ 6 x 10^-10 Msol/yr (assuming a wind speed of 100 km/s), or expansion into a uniform medium with density n_CSM <~ 6 cm^-3. Drawing from the observed properties of non-conservative mass transfer among accreting white dwarfs, we use these limits on the density of the immediate environs to exclude a phase space of possible progenitors systems for SN 2011fe. We rule out a symbiotic progenitor system and also a system characterized by high accretion rate onto the white dwarf that is expected to give rise to optically-thick accretion winds. Assuming that a small fraction, 1%, of the mass accreted is lost from the progenitor system, we also eliminate much of the potential progenitor parameter space for white dwarfs hosting recurrent novae or undergoing stable nuclear burning. Therefore, we rule out the most popular single degenerate progenitor models for SN 2011fe, leaving a limited phase space inhabited by some double degenerate systems and exotic progenitor scenarios.

We present a generalized analytic formalism for the inverse Compton X-ray emission from hydrogen-poor supernovae and apply this framework to SN2011fe using Swift-XRT, UVOT and Chandra observations. We characterize the optical properties of SN2011fe in the Swift bands and find them to be broadly consistent with a "normal" SN Ia, however, no X-ray source is detected by either XRT or Chandra. We constrain the progenitor system mass loss rate to be lower than 2x10^-9 M_sun/yr (3sigma c.l.) for wind velocity v_w=100 km/s. Our result rules out symbiotic binary progenitors for SN2011fe and argues against Roche-lobe overflowing subgiants and main sequence secondary stars if >1% of the transferred mass is lost at the Lagrangian points. Regardless of the density profile, the X-ray non-detections are suggestive of a clean environment (particle density < 150 cm-3) for (2x10^15<R<5x10^16) cm around the progenitor site. This is either consistent with the bulk of material being confined within the binary system or with a significant delay between mass loss and supernova explosion. We furthermore combine X-ray and radio limits from Chomiuk et al. 2012 to constrain the post shock energy density in magnetic fields. Finally, we searched for the shock breakout pulse using gamma-ray observations from the Interplanetary Network and find no compelling evidence for a supernova-associated burst. Based on the compact radius of the progenitor star we estimate that the shock break out pulse was likely not detectable by current satellites.

The explosion of a Type Ia supernova, SN 2011ef, in the nearby Pinwheel galaxy (M101 at 6.4 Mpc) provides an opportunity to study pre-explosion images and search for the progenitor, which should consist of a white dwarf (WD), possibly surrounded by an accretion disk, in orbit with another star. We report on our use of deep Chandra observations to limit the luminosity and temperature of the pre-explosion white dwarf (WD). It is found that if the spectrum was a blackbody, then WDs of highest possible temperatures and luminosities are excluded but, even if the WD was emitting at the Eddington luminosity, values of kT less than roughly 60 eV are permitted. This allows the progenitor to be an accreting nuclear-burning WD with an expanded photosphere. Pre-SN HST observations were used to derive a lower limit of about 10 eV for the expanded photosphere. Li et al.\, (2011) have already ruled out the possibility of a giant donor. We consider the combined emission from the WD, disk, and donor, and find that even the combined emission from a bright subgiant, WD and disk would not likely have been observed prior to explosion, and neither would some local candidates for the nuclear-burning WD model.

We present absorption line analysis of the outflow in the quasar IRAS F04250-5718. Far-ultraviolet data from the Cosmic Origins Spectrograph onboard the Hubble Space Telescope reveal intrinsic narrow absorption lines from high ionization ions (e.g., C IV, N V, and O VI) as well as low ionization ions (e.g., C II and Si III). We identify three kinematic components with central velocities ranging from ~-50 to ~-230 km/s. Velocity dependent, non-black saturation is evident from the line profiles of the high ionization ions. From the non-detection of absorption from a metastable level of C II, we are able to determine that the electron number density in the main component of the outflow is < 30 per cubic cm. Photoionization analysis yields an ionization parameter log U ~ -1.6 +/- 0.2, which accounts for changes in the metallicity of the outflow and the shape of the incident spectrum. We also consider solutions with two ionization parameters. If the ionization structure of the outflow is due to photoionization by the active galactic nucleus, we determine that the distance to this component from the central source is > 3 kpc. Due to the large distance determined for the main kinematic component, we discuss the possibility that this outflow is part of a galactic wind.

The {\alpha}-formalism is a common way to parametrize the common envelope interaction between a giant star and a more compact companion. The {\alpha} parameter describes the fraction of orbital energy released by the companion that is available to eject the giant star's envelope. By using new, detailed stellar evolutionary calculations we derive a user-friendly prescription for the {\lambda} parameter and an improved approximation for the envelope binding energy, thus revising the {\alpha} equation. We then determine {\alpha} both from simulations and observations in a self consistent manner. By using our own stellar structure models as well as population considerations to reconstruct the primary's parameters at the time of the common envelope interaction, we gain a deeper understanding of the uncertainties. We find that systems with very low values of q (the ratio of the companion's mass to the mass of the primary at the time of the common envelope interaction) have higher values of {\alpha}. A fit to the data suggests that lower mass companions are left at comparable or larger orbital separations to more massive companions. We conjecture that lower mass companions take longer than a stellar dynamical time to spiral in to the giant's core, and that this is key to allowing the giant to use its own thermal energy to help unbind its envelope. As a result, although systems with light companions might not have enough orbital energy to unbind the common envelope, they might stimulate a stellar reaction that results in the common envelope ejection.

Using high resolution VLT spectra, we study the multi-component outflow systems of two quasars exhibiting intrinsic Fe II absorption (QSO 2359-1241 and SDSS J0318-0600). From the extracted ionic column densities and using photoionization modeling we determine the gas density, total column density, and ionization parameter for several of the components. For each object the largest column density component is also the densest, and all other components have densities of roughly 1/4 of that of the main component. We demonstrate that all the absorbers lie roughly at the same distance from the source. Further, we calculate the total kinetic luminosities and mass outflow rates of all components and show that these quantities are dominated by the main absorption component.

We present high spectral resolution VLT observations of the BAL quasar SDSS J0318-0600. This high quality data set allows us to extract accurate ionic column densities and determine an electron number density of n_e=10^3.3 +/- 0.2 cm^-3 for the main outflow absorption component. The heavily reddened spectrum of SDSS J0318-0600 requires purely silicate dust with a reddening curve characteristic of predominately large grains, from which we estimate the bolometric luminosity. We carry out photoionization modeling to determine the total column density, ionization parameter and distance of the gas and find that the photionization models suggest abundances greater than solar. Due to the uncertainty in the location of the dust extinction, we arrive at two viable distances for the main ouflow component from the central source, 6 and 18 kpc, where we consider the 6 kpc location as somewhat more physically plausable. Assuming the canonical global covering of 20% for the outflow and a distance of 6 kpc, our analysis yields a mass flux of 120 M_sun yr^-1 and a kinetic luminosity that is ~0.1% of the bolometric luminosity of the object. Should the dust be part of the outflow, then these values are ~4x larger. The large mass flux and kinetic luminosity make this outflow a significant contributor to AGN feedback processes.

We present a detailed analysis of the Astrophysical Research Consortium 3.5 m telescope spectrum of QSO SDSS J0838+2955. The object shows three broad absorption line (BAL) systems at 22,000, 13,000, and 4900 km s^-1 blueshifted from the systemic redshift of z=2.043. Of particular interest is the lowest velocity system that displays absorption from low-ionization species such as Mg II, Al II, Si II, Si II*, Fe II and Fe II*. Accurate column densities were measured for all transitions in this lowest velocity BAL using an inhomogeneous absorber model. The ratio of column densities of Si II* and Fe II* with respect to their ground states gave an electron number density of log n_e (cm^-3) = 3.75 +/- 0.22 for the outflow. Photoionization modeling with careful regards to chemical abundances and the incident spectral energy distribution predicts an ionization parameter of log U_H = -1.93 +/- 0.21 and a hydrogen column density of log N_H (cm^-2) = 20.80 +/- 0.28. This places the outflow at 3.3+1.5-1.0 kpc from the central AGN. Assuming that the fraction of solid angle subtended by the outflow is 0.2, these values yield a kinetic luminosity of (4.5+3.1-1.8) x 10^45 erg s^-1, which is (1.4+1.1-0.6)% the bolometric luminosity of the QSO itself. Such large kinetic luminosity suggests that QSO outflows are a major contributor to AGN feedback mechanisms.

We analyze the physical conditions of the outflow seen in QSO 2359-1241 (NVSS J235953-124148), based on high resolution spectroscopic VLT observations. This object was previously studied using Keck/HIRES data. The main improvement over the HIRES results is our ability to accurately determine the number density of the outflow. For the major absorption component, level population from five different Fe II excited level yields n_H=10^4.4 cm^-3 with less than 20% scatter. We find that the Fe ii absorption arises from a region with roughly constant conditions and temperature greater than 9000 K, before the ionization front where temperature and electron density drop. Further, we model the observed spectra and investigate the effects of varying gas metalicities and the spectral energy distribution of the incident ionizing radiation field. The accurately measured column densities allow us to determine the ionization parameter log(U) = -2.4 and total column density of the outflow (log(N_H) = 20.6 cm^-2). Combined with the number density finding, these are stepping stones towards determining the mass flux and kinetic luminosity of the outflow, and therefore its importance to AGN feedback processes.

We present high resolution spectroscopic VLT observations of the outflow seen in QSO 2359-1241. These data contain absorption troughs from five resonance Fe II lines with a resolution of ~7 km/s and signal-to-noise ratio per resolution element of order 100. We use this unprecedented high quality data set to investigate the physical distribution of the material in front of the source, and by that determine the column densities of the absorbed troughs. We find that the apparent optical depth model gives a very poor fit to the data and greatly underestimates the column density measurements. Power-law distributions and partial covering models give much better fits with some advantage to power-law models, while both models yield similar column density estimates. The better fit of the power-law model solves a long standing problem plaguing the partial covering model when applied to large distance scale outflow: How to obtain a velocity dependent covering factor for an outflow situated at distances thousands of time greater than the size of the AGN emission source. This problem does not affect power-law models. Therefore, based on the better fit and plausibility of the physical model, we conclude that in QSO 2359-1241, the outflow covers the full extent of the emission source but in a non-homogeneous way.

We present a population synthesis calculation to derive the total number of planetary nebulae (PN) in the Galaxy that descend from single stars and stars in binary systems. Using the most recent literature results on galactic and stellar formation as well as stellar evolution, we predict the total number of galactic PNe with radii <0.9 pc to be (46,000 +/- 13,000). We do not claim this to be the complete population, since there can be visible PNe with radii larger than this limit. However, by taking this limit, we make our predicted population inherently comparable to the observationally-based value of Peimbert, who determined (7200 +/- 1800) PNe should reside in the Galaxy today. Our prediction is discrepant with the observations at the 2.9-sigma level, a disagreement which we argue is meaningful in view of our specific treatment of the uncertainty. We conclude that it is likely that only a subset of the stars thought to be capable of making a visible PN, actually do. In the second paper in this series, an argument will be presented that the bulk of the galactic PN population might be better explained if only binaries produce PNe. The predicted PN formation rate density from single stars and binaries is (1.1 +/- 0.5) x 10^{-12} PN/yr per cubic pc in the local neighborhood. This number is lower than the most recent PN birthrate density estimates of 2.1 x 10^{-12} PN/yr per cubic pc, which are based on local PN counts and the PN distance scale, but more in line with the white dwarf birthrate densities determined by Liebert et al. ((1.0 +/- 0.25) x 10^{-12} WD/yr per cubic pc). The predicted PN birthrate density will be revised down, if we assume that only binaries make PNe. This revision will imply that the PN distance scale has to be revised to larger values.

The common envelope interaction is responsible for evolved close binaries. Among them are a minority of central stars of planetary nebula (PN). Recent observational results, however, point to most PN actually being in binary systems. We therefore ask the question if it is feasible that most, or even all Galactic PN derive from a common envelope interaction. Our recent calculation finds that if all single and binary primary stars with mass between ~1-8 Mo eject a PN, there would be many more PN in the galaxy than observed. On the other hand, the predicted number of post-common envelope PN is more in agreement with the total number of PN in the Galaxy. This is a new indication that binary interactions play a functional role in the creation of PN and an encouragement to intensify efforts to detect binary companions.