Of the seven known double neutron stars (DNS) with precisely measure masses in the Milky Way that will merge within a Hubble time, all but one has a mass ratio, $q$, close to unity. Recently, precise measurements of three post-Keplerian parameters in the DNS J1913$+$1102 constrain this system to have a significantly non-unity mass ratio of 0.78$\pm$0.03. One may be tempted to conclude that approximately one out of seven (14\%) DNS mergers detected by gravitational wave observatories will have mass ratios significantly different from unity. However J1913$+$1102 has a relatively long merger time of 470 Myr. We show that when merger times and observational biases are taken into account, the population of Galactic DNSs imply that $\simeq98\%$ of all merging DNSs will have $q>$0.9. We then apply two separate fitting formulas informed by 3D hydrodynamic simulations of DNS mergers to our results on Galactic DNS masses, finding that either $\simeq$0.004 ${M_{\odot}}$ or $\simeq$0.010 ${M_{\odot}}$ of material will be ejected at merger, depending on which formula is used. These ejecta masses have implications for both the peak bolometric luminosities of electromagnetic counterparts (which we find to be $\sim$10$^{41}$ erg s$^{-1}$) as well as the $r$-process enrichment of the Milky Way.

We investigate a suspected very massive star in one of the most metal-poor dwarf galaxies, PHL~293B. Excitingly, we find the sudden disappearance of the stellar signatures from our 2019 spectra, in particular the broad H lines with P~Cygni profiles that have been associated with a massive luminous blue variable (LBV) star. Such features are absent from our spectra obtained in 2019 with the ESPRESSO and X-shooter instruments of the ESO's VLT. We compute radiative transfer models using CMFGEN that fit the observed spectrum of the LBV and are consistent with ground-based and archival Hubble Space Telescope photometry. Our models show that during 2001--2011 the LBV had a luminosity $L_* = 2.5-3.5 \times 10^6 ~L_{\odot}$, a mass-loss rate $\dot{M} = 0.005-0.020 ~M_{\odot}$~yr$^{-1}$, a wind velocity of 1000~km~s$^{-1}$, and effective and stellar temperatures of $T_\mathrm{eff} = 6000-6800$~K and $T_\mathrm{*}=9500-15000$~K. These stellar properties indicate an eruptive state. We consider two main hypotheses for the absence of the broad emission components from the spectra obtained since 2011. One possibility is that we are seeing the end of an LBV eruption of a surviving star, with a mild drop in luminosity, a shift to hotter effective temperatures, and some dust obscuration. Alternatively, the LBV could have collapsed to a massive black hole without the production of a bright supernova.

We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$\alpha$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries.

We present the discovery and high-cadence follow-up observations of SN 2018ivc, an unusual Type II supernova that exploded in NGC 1068 (D=10.1 Mpc). The light curve of SN 2018ivc declines piecewise-linearly, changing slope frequently, with four clear slope changes in the first 30 days of evolution. This rapidly changing light curve indicates that interaction between the circumstellar material and ejecta plays a significant role in the evolution. Circumstellar interaction is further supported by a strong X-ray detection. The spectra are rapidly evolving and dominated by hydrogen, helium, and calcium emission lines. We identify a rare high-velocity emission-line feature blueshifted at ~7800 km/s (in Ha, Hb, Pb, Pg, HeI, CaII), which is visible from day 18 until at least day 78 and could be evidence of an asymmetric progenitor or explosion. From the overall similarity between SN 2018ivc and SN 1996al, the \Ha{} equivalent width of its parent HII region, and constraints from pre-explosion archival Hubble Space Telescope images, we find that the progenitor of SN 2018ivc could be as massive as 52 Msun but is more likely <12 Msun. SN 2018ivc demonstrates the importance of the early discovery and rapid follow-up observations of nearby supernovae to study the physics and progenitors of these cosmic explosions.

The luminous Type IIn SN 2010jl shows strong signs of interaction between the SN ejecta and dense circumstellar material. Dust may be present in the unshocked ejecta, the cool, dense shell between the shocks in the interaction region, or in the circumstellar medium. We present and model new optical and infrared photometry and spectroscopy of SN 2010jl from 82 to 1367 days since explosion. We evaluate the photometric and spectroscopic evolution using the radiative transfer codes MOCASSIN and DAMOCLES, respectively. We propose an interaction scenario and investigate the resulting dust formation scenarios and dust masses. We find that SN 2010jl has been continuously forming dust based on the evolution of its infrared emission and optical spectra. There is evidence for pre-existing dust in the circumstellar medium as well as new dust formation in the cool, dense shell and/or ejecta. We estimate that 0.005--0.01M$_{\odot}$ of predominantly carbon dust grains has formed in SN 2010jl by ~1400d post-outburst.

The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

We present multi-wavelength photometric and spectroscopic observations of SN 2019ein, a High Velocity Type Ia supernova discovered in the nearby galaxy NGC 5353 with a two day non-detection limit. SN 2019ein exhibited some of the highest measured expansion velocities of any Type Ia supernovae, with a Si II absorption minimum blueshifted by 24,000 km s$^{-1}$ at 14 days before peak brightness. More unusually, we observed the emission components of the P Cygni profiles to be blueshifted upwards of 10,000 km s$^{-1}$ with respect to the host galaxy redshift before B-band maximum light. This blueshift, among the highest in a sample of other Type Ia supernovae, is greatest at our earliest spectroscopic epoch and subsequently decreases toward maximum light. We discuss models of progenitor systems and explosion mechanisms that could explain these extreme absorption and emission velocities. Radio observations beginning 14 days before B-band maximum light yield non-detections at the position of SN 2019ein, which rules out symbiotic progenitor systems, most models of white dwarfs powering fast optically-thick accretion winds, and optically thin shells of mass $\lesssim 10^{-6}$ M$_\odot$ at radii $< 100$ AU. Comparing our spectra to models and observations of other High Velocity Type Ia supernovae, we find that the observational properties of SN 2019ein match those of a delayed detonation explosion. We propose that the high velocity ejecta and blueshifted emission may be the result of abundance enhancements due to ejecta mixing in an asymmetric explosion, or optical depth effects in the photosphere of the ejecta at early times. These findings may provide evidence for common explosion mechanisms and ejecta geometries among High Velocity Type Ia supernovae.

We present the final sample of 98 detached double white dwarf (WD) binaries found in the Extremely Low Mass (ELM) Survey, a spectroscopic survey targeting <0.3 Msun He-core WDs completed in the Sloan Digital Sky Survey footprint. Over the course of the survey we observed ancillary low mass WD candidates like GD278, which we show is a P=0.19 d double WD binary, as well as candidates that turn out to be field blue straggler/subdwarf A-type stars with luminosities too large to be WDs given their Gaia parallaxes. Here, we define a clean sample of ELM WDs that is complete within our target selection and magnitude range 15<g_0<20 mag. The measurements are consistent with 100% of ELM WDs being 0.0089 < P < 1.5 d double WD binaries, 35% of which belong to the Galactic halo. We infer these are mostly He+CO WD binaries given the measurement constraints. The merger rate of the observed He+CO WD binaries exceeds the formation rate of stable mass transfer AM CVn binaries by a factor of 25, and so the majority of He+CO WD binaries must experience unstable mass transfer and merge. The shortest-period systems like J0651+2844 are signature {\it LISA} verification binaries that can be studied with gravitational waves and light.

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.

The formation and evolution of binary stars is a critical component of several fields in astronomy. The most numerous sources for gravitational wave observatories are inspiraling and/or merging compact binaries, while binary stars are present in nearly every electromagnetic survey regardless of the target population. Simulations of large binary populations serve to both predict and inform observations of electromagnetic and gravitational wave sources. Binary population synthesis is a tool that balances physical modeling with simulation speed to produce large binary populations on timescales of days. We present a community-developed binary population synthesis suite: COSMIC which is designed to simulate compact-object binary populations and their progenitors. As a proof of concept, we simulate the Galactic population of compact binaries and their gravitational wave signal observable by the Laser Interferometer Space Antenna (LISA). We find that $\sim10^8$ compact binaries reside in the Milky Way today, while $\sim10^4$ of them may be resolvable by LISA.

The composition of comets in the Solar System come in multiple groups thought to encode information about their formation in different regions of the outer protosolar disk. The recent discovery of the second interstellar object, 2I/Borisov, allows for spectroscopic investigations into its gas content and a preliminary classification of it within the Solar System comet taxonomies to test the applicability of planetesimal formation models to other stellar systems. We present spectroscopic and imaging observations from 2019 September 20th to October 26th at the Bok, MMT, and LBT telescopes. We identify CN in the comet's spectrum and set precise upper limits on the abundance of C2 on all dates. We use a Haser model to convert our integrated fluxes to production rates and find Q(CN) = 5.0 +/- 2.0 * 10^24 mol/s on September 20th and Q(CN) = 1.1 - 1.9 * 10^24 mol/s on later dates, both consistent with contemporaneous observations. We set our lowest upper limit on a C2 production rate, Q(C2) < 1.6 * 10^23 mol/s, on October 10th. The measured ratio upper limit for that date, Q(C2)/Q(CN) < 0.095 indicates that 2I/Borisov is strongly in the (carbon chain) 'depleted' taxonomic group. The only comparable Solar System comets have detected ratios near this limit, making 2I/Borisov statistically likely to be more depleted than any known comet. Most 'depleted' comets are Jupiter Family Comets, perhaps indicating a similiarity in formation conditions between the most depleted of the JFCs and 2I/Borisov. More work is needed to understand the applicability of our knowledge of Solar System comet taxonomies onto interstellar objects, and we discuss future work that could help clarify the usefulness of the approach.

Using a Milky Way double neutron star (DNS) merger rate of 210 Myr$^{-1}$, as derived by the Laser Interferometer Gravitational-Wave Observatory (LIGO), we demonstrate that the Laser Interferometer Space Antenna (LISA) will detect on average 240 (330) DNSs within the Milky Way for a 4-year (8-year) mission with a signal-to-noise ratio greater than 7. Even adopting a more pessimistic rate of 42 Myr$^{-1}$, as derived by the population of Galactic DNSs, we find a significant detection of 46 (65) Milky Way DNSs. These DNSs can be leveraged to constrain formation scenarios. In particular, traditional NS-discovery methods using radio telescopes are unable to detect DNSs with $P_{\rm orb}$ $\lesssim$1 hour (merger times $\lesssim$10 Myr). If a fast-merging channel exists that forms DNSs at these short orbital periods, LISA affords, perhaps, the only opportunity to observationally characterize these systems; we show that toy models for possible formation scenarios leave unique imprints on DNS orbital eccentricities, which may be measured by LISA for values as small as $\sim$10$^{-2}$.

In astrometric binaries, the presence of a dark, unseen star can be inferred from the gravitational pull it induces on its luminous binary companion. While the orbit of such binaries can be characterized with precise astrometric measurements, constraints made from astrometry alone are not enough to measure the component masses. In this work, we determine the precision with which Gaia can astrometrically measure the orbits and -- with additional observations -- the component masses, for luminous stars hosting hidden companions. Using realistic mock Gaia observations, we find that Gaia can precisely measure the orbits of binaries hosting hidden brown-dwarfs out to tens of pc and hidden white dwarf and neutron star companions at distances as far as several hundred pc. Heavier black hole companions may be measured out to 1 kpc or farther. We further determine how orbital period affects this precision, finding that Gaia can characterize orbits with periods as short as 10 days and as long as a few 10$^3$ days, with the best measured orbits having periods just short of Gaia's mission lifetime. Extending Gaia's nominal five-year mission lifetime by an additional five years not only allows for the measurement of longer period orbits, but those longer period binaries can be seen at even greater distances.

If the progenitor of GW170817 harbored a pulsar, then a Poynting flux dominated bow-shock cavity would have been expected to form around the traveling binary. The characteristic size of this evacuated region depends strongly on the spin-down evolution of the pulsar companion, which in turn depends on the merging timescale of the system. If this evacuated region is able to grow to a sufficiently large scale, then the deceleration of the jet, and thus the onset of the afterglow, would be noticeably delayed. The first detection of afterglow emission, which was uncovered 9.2 days after the $\gamma$-ray burst trigger, can thus be used to constrain the size of a pre-existing pulsar-wind cavity. We use this information, together with a model of the jet to place limits on the presence of a pulsar in GW170817 and discuss the derived constraints in the context of the observed double neutron star binary population. We find that the majority of Galactic systems that are close enough to merge within a Hubble time would have carved a discernibly large pulsar-wind cavity, inconsistent with the onset timescale of the X-ray afterglow of GW170817. Conversely, the recently detected system J1913+1102, which host a low-luminosity pulsar, provides a congruous Milky Way analog of GW170817's progenitor model. This study highlights the potential of the proposed observational test for gaining insight into the origin of double neutron star binaries, in particular if the properties of Galactic systems are representative of the overall merging population.

In the past five years, the number of known double neutron stars (DNS) in the Milky Way has roughly doubled. We argue that the observed sample can be split into three distinct sub-populations based on their orbital characteristics: (i) short-period, low-eccentricity binaries; (ii) wide binaries; and (iii) short-period, high-eccentricity binaries. These sub-populations also exhibit distinct spin period and spindown rate properties. We focus on sub-population (iii), which contains the Hulse-Taylor binary. Contrary to previous analysis, we demonstrate that, if they are the product of primordial binary evolution, the $P_{\rm orb}$ and $e$ distribution of these systems requires that the second-born NSs must have been formed with small natal kicks ($\lesssim$25 km s$^{-1}$) and have pre-SN masses narrowly distributed around 3.2 M$_{\odot}$. These constraints challenge binary evolution theory and further predict closely aligned spin and orbital axes, inconsistent with the Hulse-Taylor binary's measured spin-orbit misalignment angle of $\approx$20$^{\circ}$. Motivated by the similarity of these DNSs to B2127+11C, a DNS residing in the globular cluster M15, we argue that this sub-population is consistent with being formed in, and then ejected from, globular clusters. This scenario provides a pathway for the formation and merger of DNSs in stellar environments without recent star formation, as observed in the host galaxy population of short gamma ray bursts and the recent detection by LIGO of a merging DNS in an old stellar population.

Over forty years of research suggests that the common envelope phase, in which an evolved star engulfs its companion upon expansion, is the critical evolutionary stage forming short-period, compact-object binary systems, such as coalescing double compact objects, X-ray binaries, and cataclysmic variables. In this work, we adapt the one-dimensional hydrodynamic stellar evolution code, MESA, to model the inspiral of a 1.4M$_{\odot}$ neutron star (NS) inside the envelope of a 12$M_{\odot}$ red supergiant star. We self-consistently calculate the drag force experienced by the NS as well as the back-reaction onto the expanding envelope as the NS spirals in. Nearly all of the hydrogen envelope escapes, expanding to large radii ($\sim$10$^2$ AU) where it forms an optically thick envelope with temperatures low enough that dust formation occurs. We simulate the NS orbit until only 0.8M$_{\odot}$ of the hydrogen envelope remains around the giant star's core. Our results suggest that the inspiral will continue until another $\approx$0.3M$_{\odot}$ are removed, at which point the remaining envelope will retract. Upon separation, a phase of dynamically stable mass transfer onto the NS accretor is likely to ensue, which may be observable as an ultraluminous X-ray source. The resulting binary, comprised of a detached 2.6M$_{\odot}$ helium-star and a NS with a separation of 3.3-5.7R$_{\odot}$, is expected to evolve into a merging double neutron-star, analogous to those recently detected by LIGO/Virgo. For our chosen combination of binary parameters, our estimated final separation (including the phase of stable mass transfer) suggests a very high $\alpha_{\rm CE}$-equivalent efficiency of $\approx$5.

A new all-sky visible and Near-InfraRed (NIR) space astrometry mission with a wavelength cutoff in the K-band is not just focused on a single or small number of key science cases. Instead, it is extremely broad, answering key science questions in nearly every branch of astronomy while also providing a dense and accurate visible-NIR reference frame needed for future astronomy facilities. For almost 2 billion common stars the combination of Gaia and a new all-sky NIR astrometry mission would provide much improved proper motions, answering key science questions -- from the solar system and stellar systems, including exoplanet systems, to compact galaxies, quasars, neutron stars, binaries and dark matter substructures. The addition of NIR will result in up to 8 billion newly measured stars in some of the most obscured parts of our Galaxy, and crucially reveal the very heart of the Galactic bulge region. In this white paper we argue that rather than improving on the accuracy, a greater overall science return can be achieved by going deeper than Gaia and by expanding the wavelength range to the NIR.

Type IIP supernovae (SNe IIP), which represent the most common class of core-collapse (CC) SNe, show a rapid increase in continuum polarization just after entering the tail phase. This feature can be explained by a highly asymmetric helium core, which is exposed when the hydrogen envelope becomes transparent. Here we report the case of a SN IIP (SN~2017gmr) that shows an unusually early rise of the polarization, $\gtrsim 30$ days before the start of the tail phase. This implies that SN~2017gmr is an SN IIP that has very extended asphericity. The asymmetries are not confined to the helium core, but reach out to a significant part of the outer hydrogen envelope, hence clearly indicating a marked intrinsic diversity in the aspherical structure of CC explosions. These observations provide new constraints on the explosion mechanism, where viable models must be able to produce such extended deviations from spherical symmetry, and account for the observed geometrical diversity.

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.

The optical and ultra-violet broadband photometric and spectroscopic observations of the Type II supernova (SN) 2016gfy are presented. The $V$-band light curve (LC) shows a distinct plateau phase with a slope, $s_2$ $\sim$ 0.12 mag / 100 d and a duration of 90 $\pm$ 5 d. Detailed analysis of SN 2016gfy provided a mean $\rm^{56}Ni$ mass of 0.033 $\pm$ 0.003 $\rm M_{\odot}$, a progenitor radius of $\sim$ 350-700 $\rm R_{\odot}$, a progenitor mass of $\sim$ 12-15 $\rm M_{\odot}$ and an explosion energy of 0.9-1.4$\rm \times 10^{51}\ erg\ s^{-1}$. The P-Cygni profile of H$\rm \alpha$ in the early phase spectra ($\sim$ 11-21 d) shows a boxy emission. Assuming that this profile arises from the interaction of the SN ejecta with the pre-existing circumstellar material (CSM), it is inferred that the progenitor underwent a recent episode (30-80 years prior to the explosion) of enhanced mass loss. Numerical modeling suggests that the early LC peak is reproduced better with an existing CSM of 0.15 $\rm M_{\odot}$ spread out to $\sim$ 70 AU. A late-plateau bump is seen in the $VRI$ LCs during $\sim$ 50-95 d. This bump is explained as a result of the CSM interaction and/or partial mixing of radioactive $\rm ^{56}Ni$ in the SN ejecta. Using strong-line diagnostics, a sub-solar oxygen abundance is estimated for the supernova H II region (12 + log(O/H) = 8.50 $\pm$ 0.11), indicating an average metallicity for the host of a Type II SN. A star formation rate of $\sim$ 8.5 $\rm M_{\odot}\ yr^{-1}$ is estimated for NGC 2276 using the archival $GALEX$ FUV data.

We present high-cadence ultraviolet (UV), optical, and near-infrared (NIR) data on the luminous Type II-P supernova SN 2017gmr from hours after discovery through the first 180 days. SN 2017gmr does not show signs of narrow, high-ionization emission lines in the early optical spectra, yet the optical lightcurve evolution suggests that an extra energy source from circumstellar medium (CSM) interaction must be present for at least 2 days after explosion. Modeling of the early lightcurve indicates a ~500R$_{\odot}$ progenitor radius, consistent with a rather compact red supergiant, and late-time luminosities indicate up to 0.130 $\pm$ 0.026 M$_{\odot}$ of $^{56}$Ni are present, if the lightcurve is solely powered by radioactive decay, although the $^{56}$Ni mass may be lower if CSM interaction contributes to the post-plateau luminosity. Prominent multi-peaked emission lines of H$\alpha$ and [O I] emerge after day 154, as a result of either an asymmetric explosion or asymmetries in the CSM. The lack of narrow lines within the first two days of explosion in the likely presence of CSM interaction may be an example of close, dense, asymmetric CSM that is quickly enveloped by the spherical supernova ejecta.

We show that the recently discovered binary 2M05215658+4359220 (2M0521), comprised of a giant star (GS) orbiting a suspected black hole (BH) in a ~80 day orbit, may be instrumental in shedding light on uncertain BH-formation physics and can be a test case for studying wind accretion models. Using binary population synthesis with a realistic prescription for the star formation history and metallicity evolution of the Milky Way, we analyze the evolution of binaries containing compact objects (COs) in orbit around GSs with properties similar to 2M0521. We find ~100-1000 CO-GS binaries in the Milky Way observable by Gaia, and 0-12 BH-GS and 0-1 neutron star-GS binaries in the Milky Way with properties similar to 2M0521. We find that all CO-GSs with Porb<5 yr, including 2M0521, go through a common envelope (CE) and hence form a class of higher mass analogs to white dwarf post-CE binaries. We further show how the component masses of 2M0521-like binaries depend strongly on the supernova-engine model we adopt. Thus, an improved measurement of the orbit of 2M0521, imminent with Gaia's third data release, will strongly constrain its component masses and as a result inform supernova-engine models widely used in binary population synthesis studies. These results have widespread implications for the origins and properties of CO binaries, especially those detectable by LIGO and LISA. Finally, we show that the reported X-ray non-detection of 2M0521 is a challenge for wind accretion theory, making 2M0521-like CO-GS binaries a prime target for further study with accretion models.

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.

We present extensive optical photometric and spectroscopic observations, from 4 to 482 days after explosion, of the Type II-plateau (II-P) supernova (SN) 2017eaw in NGC 6946. SN 2017eaw is a normal SN II-P intermediate in properties between, for example, SN 1999em and SN 2012aw and the more luminous SN 2004et, also in NGC 6946. We have determined that the extinction to SN 2017eaw is primarily due to the Galactic foreground and that the SN site metallicity is likely subsolar. We have also independently confirmed a tip-of-the-red-giant-branch (TRGB) distance to NGC 6946 of 7.73+/-0.78 Mpc. The distances to the SN that we have also estimated via both the standardized candle method and expanding photosphere method corroborate the TRGB distance. We confirm the SN progenitor identity in pre-explosion archival Hubble Space Telescope (HST) and Spitzer Space Telescope images, via imaging of the SN through our HST Target of Opportunity program. Detailed modeling of the progenitor's spectral energy distribution indicates that the star was a dusty, luminous red supergiant consistent with an initial mass of ~15 Msuns.

The formation and evolution of double neutron stars (DNS) have traditionally been studied using binary population synthesis. In this work, we take an alternative approach by focusing only on the second supernova (SN) forming the DNS and the subsequent orbital decay and merger due to gravitational wave radiation. Using analytic and numerical methods, we explore how different NS natal kick velocity distributions, pre-SN orbital separations, and progenitor He-star masses affect the post-SN orbital periods, eccentricities, merger times, systemic velocities, and distances traveled by the system before merging. Comparison with the set of 17 known DNSs in the Milky Way shows that DNSs have pre-SN orbital separations ranging between 1 and 44 $R_{\odot}$. Those DNSs with pre-SN separations $\sim$1 $R_{\odot}$ have merger time distributions that peak $\sim$10-100 Myr after formation, regardless of the kick velocity received by the NS. These DNSs are typically formed with systemic velocities $\sim$10$^2$ km s$^{-1}$ and may travel $\sim$1-10 kpc before merging. Depending on progenitor mass of the second-born NS, the short merger time can account for the $r$-process enrichment observed in compact stellar systems such as ultra-faint dwarf galaxies. For Milky Way-mass galaxies only DNSs with the tightest pre-SN orbits and large kick velocities ($\gtrsim$10$^2$ km s$^{-1}$) can escape. However, those DNSs that do escape may travel as far as $\sim$Mpc before merging, which as previous studies have pointed out has implications for identifying the host galaxies to short gamma ray bursts and gravitational wave events.

Supernovae are among the most powerful and influential explosions in the universe. They are also ideal multi-messenger laboratories to study extreme astrophysics. However, many fundamental properties of supernovae related to their diverse progenitor systems and explosion mechanisms remain poorly constrained. Here we outline how late-time spectroscopic observations obtained during the nebular phase (several months to years after explosion), made possible with the next generation of Extremely Large Telescopes, will facilitate transformational science opportunities and rapidly accelerate the community towards our goal of achieving a complete understanding of supernova explosions. We highlight specific examples of how complementary GMT and TMT instrumentation will enable high fidelity spectroscopy from which the line profiles and luminosities of elements tracing mass loss and ejecta can be used to extract kinematic and chemical information with unprecedented detail, for hundreds of objects. This will provide uniquely powerful constraints on the evolutionary phases stars may experience approaching a supernova explosion; the subsequent explosion dynamics; their nucleosynthesis yields; and the formation of compact objects that may act as central engines.

In order to understand the contribution of core-collapse supernovae to the dust budget of the early universe, it is important to understand not only the mass of dust that can form in core-collapse supernovae but also the location and rate of dust formation. SN 2005ip is of particular interest since dust has been inferred to have formed in both the ejecta and the post-shock region behind the radiative reverse shock. We have collated eight optical archival spectra that span the lifetime of SN 2005ip and we additionally present a new X-shooter optical-near-IR spectrum of SN 2005ip at 4075d post-discovery. Using the Monte Carlo line transfer code DAMOCLES, we have modelled the blueshifted broad and intermediate width H$\alpha$, H$\beta$ and He I lines from 48d to 4075d post-discovery using an ejecta dust model. We find that dust in the ejecta can account for the asymmetries observed in the broad and intermediate width H$\alpha$, H$\beta$ and He I line profiles at all epochs and that it is not necessary to invoke post-shock dust formation to explain the blueshifting observed in the intermediate width post-shock lines. Using a Bayesian approach, we have determined the evolution of the ejecta dust mass in SN 2005ip over 10 years presuming an ejecta dust model, with an increasing dust mass from ~10$^{-8}$ M$_{\odot}$ at 48d to a current dust mass of $\sim$0.1 M$_{\odot}$.

We present the star cluster catalogs for 17 dwarf and irregular galaxies in the $HST$ Treasury Program "Legacy ExtraGalactic UV Survey" (LEGUS). Cluster identification and photometry in this subsample are similar to that of the entire LEGUS sample, but special methods were developed to provide robust catalogs with accurate fluxes due to low cluster statistics. The colors and ages are largely consistent for two widely used aperture corrections, but a significant fraction of the clusters are more compact than the average training cluster. However, the ensemble luminosity, mass, and age distributions are consistent suggesting that the systematics between the two methods are less than the random errors. When compared with the clusters from previous dwarf galaxy samples, we find that the LEGUS catalogs are more complete and provide more accurate total fluxes. Combining all clusters into a composite dwarf galaxy, we find that the luminosity and mass functions can be described by a power law with the canonical index of $-2$ independent of age and global SFR binning. The age distribution declines as a power law, with an index of $\approx-0.80\pm0.15$, independent of cluster mass and global SFR binning. This decline of clusters is dominated by cluster disruption since the combined star formation histories and integrated-light SFRs are both approximately constant over the last few hundred Myr. Finally, we find little evidence for an upper-mass cutoff ($<2\sigma$) in the composite cluster mass function, and can rule out a truncation mass below $\approx10^{4.5}$M$_{\odot}$ but cannot rule out the existence of a truncation at higher masses.

The recent aLIGO/aVirgo discovery of gravitational waves from the neutron star merger (NSM) GW170817 and the follow up kilonova observations have shown that NSMs produce copious amount of r-process material. However, it is difficult to reconcile the large natal kicks and long average merging times of Double Neutron Stars (DNSs), with the levels of r-process enrichment seen in ultra faint dwarf (UFD) galaxies such as Reticulum II and Tucana III. Assuming that such dwarf systems have lost a significant fraction of their stellar mass through tidal stripping, we conclude that contrary to most current models, it is the DNSs with rather large natal kicks but very short merging timescales that can enrich UFD-type galaxies. These binaries are either on highly eccentric orbits, or form with very short separations due to an additional mass-transfer between the first-born neutron star and a naked helium star, progenitor of the second-born neutron star. These DNSs are born with a frequency that agrees with the statistics of the r-process UFDs, and merge well within the virial radius of their host halos, therefore contributing significantly to their r-process enrichment.

Stars of a common origin are thought to have similar, if not nearly identical, chemistry. Chemical tagging seeks to exploit this fact to identify Milky Way subpopulations through their unique chemical fingerprints. In this work, we compare the chemical abundances of dwarf stars in wide binaries to test the abundance consistency of stars of a common origin. Our sample of 31 wide binaries is identified from a catalog produced by cross-matching APOGEE stars with UCAC5 astrometry, and we confirm the fidelity of this sample with precision parallaxes from Gaia DR2. For as many as 14 separate elements, we compare the abundances between components of our wide binaries, finding they have very similar chemistry (typically within 0.1 dex). This level of consistency is more similar than can be expected from stars with different origins (which show typical abundance differences of 0.3-0.4 dex within our sample). For the best measured elements, Fe, Si, K, Ca, Mn, and Ni, these differences are reduced to 0.05-0.08 dex when selecting pairs of dwarf stars with similar temperatures. Our results suggest that APOGEE dwarf stars may currently be used for chemical tagging at the level of $\sim$0.1 dex or at the level of $\sim$0.05 dex when restricting for the best-measured elements in stars of similar temperatures. Larger wide binary catalogs may provide calibration sets, in complement to open cluster samples, for on-going spectroscopic surveys.

We present an exquisite, 30-min cadence Kepler (K2) light curve of the Type Ia supernova (SN Ia) 2018oh (ASASSN-18bt), starting weeks before explosion, covering the moment of explosion and the subsequent rise, and continuing past peak brightness. These data are supplemented by multi-color Pan-STARRS1 and CTIO 4-m DECam observations obtained within hours of explosion. The K2 light curve has an unusual two-component shape, where the flux rises with a steep linear gradient for the first few days, followed by a quadratic rise as seen for typical SNe Ia. This "flux excess" relative to canonical SN Ia behavior is confirmed in our $i$-band light curve, and furthermore, SN 2018oh is especially blue during the early epochs. The flux excess peaks 2.14$\pm0.04$ days after explosion, has a FWHM of 3.12$\pm0.04$ days, a blackbody temperature of $T=17,500^{+11,500}_{-9,000}$ K, a peak luminosity of $4.3\pm0.2\times10^{37}\,{\rm erg\,s^{-1}}$, and a total integrated energy of $1.27\pm0.01\times10^{43}\,{\rm erg}$. We compare SN 2018oh to several models that may provide additional heating at early times, including collision with a companion and a shallow concentration of radioactive nickel. While all of these models generally reproduce the early K2 light curve shape, we slightly favor a companion interaction, at a distance of $\sim$$2\times10^{12}\,{\rm cm}$ based on our early color measurements, although the exact distance depends on the uncertain viewing angle. Additional confirmation of a companion interaction in future modeling and observations of SN 2018oh would provide strong support for a single-degenerate progenitor system.

Supernova (SN) 2018oh (ASASSN-18bt) is the first spectroscopically-confirmed type Ia supernova (SN Ia) observed in the $Kepler$ field. The $Kepler$ data revealed an excess emission in its early light curve, allowing to place interesting constraints on its progenitor system (Dimitriadis et al. 2018, Shappee et al. 2018b). Here, we present extensive optical, ultraviolet, and near-infrared photometry, as well as dense sampling of optical spectra, for this object. SN 2018oh is relatively normal in its photometric evolution, with a rise time of 18.3$\pm$0.3 days and $\Delta$m$_{15}(B)=0.96\pm$0.03 mag, but it seems to have bluer $B - V$ colors. We construct the "uvoir" bolometric light curve having peak luminosity as 1.49$\times$10$^{43}$erg s$^{-1}$, from which we derive a nickel mass as 0.55$\pm$0.04M$_{\odot}$ by fitting radiation diffusion models powered by centrally located $^{56}$Ni. Note that the moment when nickel-powered luminosity starts to emerge is +3.85 days after the first light in the Kepler data, suggesting other origins of the early-time emission, e.g., mixing of $^{56}$Ni to outer layers of the ejecta or interaction between the ejecta and nearby circumstellar material or a non-degenerate companion star. The spectral evolution of SN 2018oh is similar to that of a normal SN Ia, but is characterized by prominent and persistent carbon absorption features. The C II features can be detected from the early phases to about 3 weeks after the maximum light, representing the latest detection of carbon ever recorded in a SN Ia. This indicates that a considerable amount of unburned carbon exists in the ejecta of SN 2018oh and may mix into deeper layers.

We present multi-epoch photometry and spectroscopy of a light echo from eta Carinae's 19th century Great Eruption. This echo shows a steady decline over a decade, sampling the 1850s plateau of the eruption. Spectra show the bulk outflow speed increasing from 150 km/s at early times, up to 600 km/s in the plateau. Later phases also develop remarkably broad emission wings indicating mass accelerated to more than 10,000 km/s. Together with other clues, this provides direct evidence for an explosive ejection. This is accompanied by a transition from narrow absorption lines to emission lines, often with broad P Cygni profiles. These changes imply that the pre-1845 luminosity spikes are distinct from the 1850s plateau. The key reason for this change may be that shock interaction dominates the plateau. The spectral evolution of eta Car closely resembles that of UGC2773-OT, which had clear signatures of shock interaction. We propose a 2-stage scenario for eta Car's eruption: (1) a slow outflow in the decades before the eruption, driven by binary interaction that produced a dense equatorial outflow, followed by (2) explosive energy injection that drove CSM interaction, powering the plateau and sweeping slower CSM into a fast shell that became the Homunculus. We discuss how this sequence could arise from a stellar merger in a triple system, leaving the eccentric binary seen today. This gives a self-consistent scenario that may explain interacting transients across a wide range of initial mass.

In our ongoing study of eta Carinae's light echoes, there is a relatively bright echo that has been fading slowly, reflecting the 1845-1858 plateau of the eruption. A separate paper discusses its detailed evolution, but here we highlight one important result: the H-alpha line shows extremely broad emission wings that reach -10,000km/s to the blue and +20,000km/s to the red. The line profile shape is inconsistent with electron scattering wings, indicating high-velocity outflowing material. These are the fastest outflow speeds ever seen in a non-terminal massive star eruption. The broad wings are absent in early phases of the eruption, but strengthen in the 1850s. These speeds are two orders of magnitude faster than the escape speed from a warm supergiant, and 5-10 times faster than winds from O-type or Wolf-Rayet stars. Instead, they are reminiscent of fast supernova ejecta or outflows from accreting compact objects, profoundly impacting our understanding of eta Car and related transients. This echo views eta Car from latitudes near the equator, so the high speed does not trace a collimated polar jet aligned with the Homunculus. Combined with fast material in the Outer Ejecta, it indicates a wide-angle explosive outflow. The fast material may constitute a small fraction of the total outflowing mass, most of which expands at 600 km/s. This is reminiscent of fast material revealed by broad absorption during the presupernova eruptions of SN2009ip.

Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this "dynamically driven double-degenerate double-detonation" (D^6) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its > 1000 km/s pre-SN orbital velocity. We perform a search for these hypervelocity runaway WDs using Gaia's second data release. In this paper, we discuss seven candidates followed up with ground-based instruments. Three sources are likely to be some of the fastest known stars in the Milky Way, with total Galactocentric velocities between 1000 and 3000 km/s, and are consistent with having previously been companion WDs in pre-SN Ia systems. However, although the radial velocity of one of the stars is > 1000 km/s, the radial velocities of the other two stars are puzzlingly consistent with 0. The combined five-parameter astrometric solutions from Gaia and radial velocities from follow-up spectra yield tentative 6D confirmation of the D^6 scenario. The past position of one of these stars places it within a faint, old SN remnant, further strengthening the interpretation of these candidates as hypervelocity runaways from binary systems that underwent SNe Ia.

We report a luminous Type II supernova, ASASSN-15nx, with a peak luminosity of M_V=-20 mag, that is between typical core-collapse supernovae and super-luminous supernovae. The post-peak optical light curves show a long, linear decline with a steep slope of 2.5 mag/100 days (i.e., an exponential decline in flux), through the end of observations at phase ~260 days. In contrast, the light curves of hydrogen rich supernovae (SNe II-P/L) always show breaks in their light curves at phase ~100 days, before settling onto Co56 radioactive decay tails with a decline rate of about 1 mag/100 days. The spectra of ASASSN-15nx do not exhibit the narrow emission-line features characteristic of Type IIn SNe, which can have a wide variety of light-curve shapes usually attributed to strong interactions with a dense circumstellar medium (CSM). ASASSN-15nx has a number of spectroscopic peculiarities, including a relatively weak and triangularly-shaped H-alpha emission profile with no absorption component. The physical origin of these peculiarities is unclear, but the long and linear post-peak light curve without a break suggests a single dominant powering mechanism. Decay of a large amount of Ni56 (M_Ni56 = 1.6 +/- 0.2 M_sun) can power the light curve of ASASSN-15nx, and the steep light-curve slope requires substantial gamma-ray escape from the ejecta, which is possible given a low-mass hydrogen envelope for the progenitor. Another possibility is strong CSM interactions powering the light curve, but the CSM needs to be sculpted to produce the unique light-curve shape and to avoid producing SN IIn-like narrow emission lines.

We report the discovery of a previously unidentified pulsar as part of a radio campaign to identify neutron star companions to low-mass white dwarfs (LMWDs) using the Robert C.\ Byrd Green Bank Telescope (GBT). PSR J0802-0955, which is coincident with the position of a WD with a mass of 0.2 solar masses, has a pulse period of 571 ms. Because of its relatively long pulse period, the lack of radial velocity (RV) variations in the radio data, and GBT's large beam size at the observing frequency of 340 MHz, we conclude that PSR J0802-0955 is unassociated with the LMWD at roughly the same position and distance.

By employing Monte Carlo random sampling, traditional binary population synthesis (BPS) offers a substantial improvement in efficiency over brute force, grid-based studies. Even so, BPS models typically require a large number of simulation realizations, a computationally expensive endeavor, to generate statistically robust results. Recent advances in statistical methods have led us to revisit the traditional approach to BPS. In this work we describe our publicly available code dart_board which combines rapid binary evolution codes, typically used in traditional BPS, with modern Markov chain Monte Carlo methods. dart_board takes a novel approach that treats the initial binary parameters and the supernova kick vector as model parameters. This formulation has several advantages, including the ability to model either populations of systems or individual binaries, the natural inclusion of observational uncertainties, and the flexible addition of new constraints which are problematic to include using traditional BPS. After testing our code with mock systems, we demonstrate the flexibility of dart_board by applying it to three examples: (i) a generic population of high mass X-ray binaries (HMXBs), (ii) the population of HMXBs in the Large Magellanic Cloud (LMC) in which the spatially resolved star formation history is used as a prior, and (iii) one particular HMXB in the LMC, Swift J0513.4-6547, in which we include observations of the system's component masses and orbital period. Although this work focuses on HMXBs, dart_board can be applied to a variety of stellar binaries including the recent detections by gravitational wave observatories of merging compact object binaries.

Featuring greatly improved astrometry and precise radial velocities (RVs), the second data release (DR2) from Gaia affords us a unique opportunity to test the validity of samples of stellar wide binaries. In a previous work, we presented a set of wide binaries identified in the joint Tycho-Gaia Astrometric Solution catalog. Here, we use DR2 to confirm the low contamination rate of this sample, thereby verifying the effectiveness of our algorithm for application to DR2 data.

We have identified a progenitor candidate in archival Hubble Space Telescope (HST) images for the Type Ic SN 2017ein in NGC 3938, pinpointing the candidate's location via HST Target-of-Opportunity imaging of the SN itself. This would be the first identification of a stellar-like object as a progenitor candidate for any Type Ic supernova to date. We also present observations of SN 2017ein during the first ~49 days since explosion. We find that SN 2017ein most resembles the well-studied Type Ic SN 2007gr. We infer that SN 2017ein experienced a total visual extinction of A_V~1.0--1.9 mag, predominantly because of dust within the host galaxy. Although the distance is not well known, if this object is the progenitor, it was likely of high initial mass, ~47--48 M_sun if a single star, or ~60--80 M_sun if in a binary system. However, we also find that the progenitor candidate could be a very blue and young compact cluster, further implying a very massive (>65 M_sun) progenitor. Furthermore, the actual progenitor might not be associated with the candidate at all and could be far less massive. From the immediate stellar environment, we find possible evidence for three different populations; if the SN progenitor was a member of the youngest population, this would be consistent with an initial mass of ~57 M_sun. After it has faded, the SN should be reobserved at high spatial resolution and sensitivity, to determine whether the candidate is indeed the progenitor.

We present a study of the dust-to-gas ratios in five nearby galaxies NGC 628 (M74), NGC 6503, NGC 7793, UGC 5139 (Holmberg I), and UGC 4305 (Holmberg II). Using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury program LEGUS (Legacy ExtraGalactic UV Survey) combined with archival HST/ACS data, we correct thousands of individual stars for extinction across these five galaxies using an isochrone-matching (reddening-free Q) method. We generate extinction maps for each galaxy from the individual stellar extinctions using both adaptive and fixed resolution techniques, and correlate these maps with neutral HI and CO gas maps from literature, including The HI Nearby Galaxy Survey (THINGS) and the HERA CO-Line Extragalactic Survey (HERACLES). We calculate dust-to-gas ratios and investigate variations in the dust-to-gas ratio with galaxy metallicity. We find a power law relationship between dust-to-gas ratio and metallicity, consistent with other studies of dust-to-gas ratio compared to metallicity. We find a change in the relation when H$_2$ is not included. This implies that underestimation of $N_{H_2}$ in low-metallicity dwarfs from a too-low CO-to-H$_2$ conversion factor $X_{CO}$ could have produced too low a slope in the derived relationship between dust-to-gas ratio and metallicity. We also compare our extinctions to those derived from fitting the spectral energy distribution (SED) using the Bayesian Extinction and Stellar Tool (BEAST) for NGC 7793 and find systematically lower extinctions from SED-fitting as compared to isochrone matching.

Various binary black hole formation channels have been proposed since the first gravitational event GW150914 was discovered by the Advanced Laser Interferometer Gravitational-Wave Observatory (AdLIGO). For all evolutionary channels based on the evolution of isolated binaries, the immediate progenitor of the binary black hole is a close binary system composed of a black hole and a helium star. We perform detailed binary evolution and systematically explore the parameter space of initial binary properties, including initial black hole and helium star masses, initial rotation of the helium star as well as metallicity. We argue that the spin of the first-born black hole at its birth is negligible ($\lesssim 0.1$), hence the second-born black hole's spin dominates the measured effective spin, $\chi_{\rm eff}$, from gravitational wave events of double black hole mergers. We find that tides can be important only when orbital periods are shorter than 2 days. Upon core collapse, the helium star produces a black hole (the second-born black hole in the system) with a spin that can span the entire range from zero to maximally spinning. We show that the bimodal distribution of the spin of the second-born black hole obtained in recent papers is mainly due to oversimplifying assumptions. We find an anti-correlation between the merging timescale of the two black holes, T$_{\rm merger}$, and the effective spin $\chi_{\rm eff}$. Finally, we provide new prescriptions for the tidal coefficient E$_2$ for both H-rich and the helium-rich stars. We predict that, with future improvements to AdLIGO's sensitivity, the sample of merging binary black hole systems will show an overdensity of sources with positive but small $\chi_{\rm eff}$. originating from lower mass black hole mergers born at low redshift.

The Legacy ExtraGalactic UV Survey (LEGUS) is a multiwavelength Cycle 21 Treasury program on the Hubble Space Telescope. It studied 50 nearby star-forming galaxies in five bands from the near UV to the I-band, combining new Wide Field Camera 3 observations with archival Advanced Camera for Surveys data. LEGUS was designed to investigate how star formation occurs and develops on both small and large scales, and how it relates to the galactic environments. In this paper we present the photometric catalogs for all the apparently single stars identified in the 50 LEGUS galaxies. Photometric catalogs and mosaicked images for all filters are available for download. We present optical and near UV color-magnitude diagrams for all the galaxies. For each galaxy we derived the distance from the tip of the red giant branch. We then used the NUV color-magnitude diagrams to identify stars more massive than 14 Mo, and compared their number with the number of massive stars expected from the GALEX FUV luminosity. Our analysis shows that the fraction of massive stars forming in star clusters and stellar associations is about constant with the star formation rate. This lack of a relation suggests that the time scale for evaporation of unbound structures is comparable or longer than 10 Myr. At low star formation rates this translates to an excess of mass in clustered environments as compared to model predictions of cluster evolution, suggesting that a significant fraction of stars form in unbound systems.

We present a moderate-resolution spectrum of the peculiar Type II supernova iPTF14hls taken on day 1153 after discovery. This spectrum reveals the clear signature of shock interaction with dense circumstellar material (CSM). We suggest that this CSM interaction may be an important clue for understanding the extremely unusual photometric and spectroscopic evolution seen over the first 600 days of iPTF14hls. The late-time spectrum shows a double-peaked intermediate-width H-alpha line indicative of expansion speeds around 1000 km/s, with the double-peaked shape hinting at a disc-like geometry in the CSM. If the CSM was highly asymmetric, perhaps in a disc or torus that was ejected from the star 3-6 years prior to explosion, then the CSM interaction could have been overrun and hidden below the SN ejecta photosphere from a wide range of viewing angles. In that case, CSM interaction luminosity would have been thermalized well below the photosphere, possibly sustaining the high luminosity without exhibiting the traditional observational signatures of strong CSM interaction (narrow H-alpha emission and X-rays). Variations in density structure of the CSM could account for the multiple rebrightenings of the lightcurve. We propose that enveloped CSM interaction as seen in some recent SNe, rather than an entirely new explosion mechanism, may be adequate to explain the peculiar evolution of iPTF14hls.

We present photometry and spectroscopy of SN2013fs and SN2013fr in the first 100 days post-explosion. Both objects showed transient, relatively narrow H$\alpha$ emission lines characteristic of SNeIIn, but later resembled normal SNeII-P or SNeII-L, indicative of fleeting interaction with circumstellar material (CSM). SN2013fs was discovered within 8hr of explosion. Its light curve exhibits a plateau, with spectra revealing strong CSM interaction at early times. It is a less luminous version of the transitional SNIIn PTF11iqb, further demonstrating a continuum of CSM interaction intensity between SNeII-P and IIn. It requires dense CSM within 6.5$\times$10$^{14}$~cm of the progenitor, from a phase of advanced pre-SN mass loss shortly before explosion. Spectropolarimetry of SN2013fs shows little continuum polarization, but noticeable line polarization during the plateau phase. SN2013fr morphed from a SNIIn at early times to a SNII-L. After the first epoch its narrow lines probably arose from host-galaxy emission, but the bright, narrow H$\alpha$ emission at early times may be intrinsic. As for SN2013fs, this would point to a short-lived phase of strong CSM interaction if proven to be intrinsic, suggesting a continuum between SNeIIn and II-L. It is a low-velocity SNII-L, like SN2009kr but more luminous. SN2013fr also developed an IR excess at later times, due to warm CSM dust that require a more sustained phase of strong pre-SN mass loss.

We investigate the relation between gas and star formation in sub-galactic regions, ~360 pc to ~1.5 kpc in size, within the nearby starburst dwarf NGC4449, in order to separate the underlying relation from the effects of sampling at varying spatial scales. Dust and gas mass surface densities are derived by combining new observations at 1.1 mm, obtained with the AzTEC instrument on the Large Millimeter Telescope, with archival infrared images in the range 8-500 micron from the Spitzer Space Telescope and the Herschel Space Observatory. We extend the dynamic range of our mm (and dust) maps at the faint end, using a correlation between the far-infrared/millimeter colors F(70)/F(1100) [and F(160)/F(1100)] and the mid-infrared color F(8)/F(24) that we establish for the first time for this and other galaxies. Supplementing our data with maps of the extinction-corrected star formation rate (SFR) surface density, we measure both the SFR-molecular gas and the SFR-total gas relations in NGC4449. We find that the SFR-molecular gas relation is described by a power law with exponent that decreases from ~1.5 to ~1.2 for increasing region size, while the exponent of the SFR-total gas relation remains constant with value ~1.5 independent of region size. We attribute the molecular law behavior to the increasingly better sampling of the molecular cloud mass function at larger region sizes; conversely, the total gas law behavior likely results from the balance between the atomic and molecular gas phases achieved in regions of active star formation. Our results indicate a non-linear relation between SFR and gas surface density in NGC4449, similar to what is observed for galaxy samples.

We present photometry, spectra, and spectropolarimetry of supernova (SN) 2012ab, mostly obtained over the course of $\sim 300$ days after discovery. SN 2012ab was a Type IIn (SN IIn) event discovered near the nucleus of spiral galaxy 2MASXJ12224762+0536247. While its light curve resembles that of SN 1998S, its spectral evolution does not. We see indications of CSM interaction in the strong intermediate-width emission features, the high luminosity (peak at absolute magnitude $M=-19.5$), and the lack of broad absorption features in the spectrum. The H$\alpha$ emission undergoes a peculiar transition. At early times it shows a broad blue emission wing out to $-14{,}000$ km $\mathrm{s^{-1}}$ and a truncated red wing. Then at late times ($>$ 100$\,$days) it shows a truncated blue wing and a very broad red emission wing out to roughly $+20{,}000$ km $\mathrm{s^{-1}}$. This late-time broad red wing probably arises in the reverse shock. Spectra also show an asymmetric intermediate-width H$\alpha$ component with stronger emission on the red side at late times. The evolution of the asymmetric profiles requires a density structure in the distant CSM that is highly aspherical. Our spectropolarimetric data also suggest asphericity with a strong continuum polarization of $\sim 1-3$% and depolarization in the H$\alpha$ line, indicating asphericity in the CSM at a level comparable to that in other SNe IIn. We estimate a mass-loss rate of $\dot{M} = 0.050\, {\rm M}_{\odot}\,\mathrm{yr^{-1}}$ for $v_{\rm pre} = 100$$\,$km$\,$$\mathrm{s^{-1}}$ extending back at least 75$\,$yr prior to the SN. The strong departure from axisymmetry in the CSM of SN 2012ab may suggest that the progenitor was an eccentric binary system undergoing eruptive mass loss.

From our recent catalog based on the first Gaia data release (TGAS), we select wide binaries in which both stars have been observed by the Radial Velocity Experiment (RAVE) or the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST). Using RAVE and LAMOST metallicities and RAVE Mg, Al, Si, Ti, and Fe abundances, we find that the differences in the metallicities and elemental abundances of components of wide binaries are consistent with being due to observational uncertainties, in agreement with previous results for smaller and more restricted samples. The metallicity and elemental abundance consistency between wide binary components presented in this work confirms their common origin and bolsters the status of wide binaries as "mini-open clusters." Furthermore, this is evidence that wide binaries are effectively co-eval and co-chemical, supporting their use for e.g., constraining age-activity-rotation relations, the initial-final mass relation for white dwarfs, and M-dwarf metallicity indicators. Additionally, we demonstrate that the common proper motion, common parallax pairs in TGAS with the most extreme separations (s $\gtrsim$ 0.1 pc) typically have inconsistent metallicities, radial velocities or both and are therefore likely to be predominantly comprised of random alignments of unassociated stars with similar astrometry, in agreement with our previous results. Finally, we propose that wide binaries form an ideal data set with which to test chemical tagging as a method to identify stars of common origin, particularly because the stars in wide binaries span a wide range of metallicities, much wider than that spanned by nearby open clusters.

Recently acquired WFC3 UV (F275W and F336W) imaging mosaics under the Legacy Extragalactic UV Survey (LEGUS) combined with archival ACS data of M51 are used to study the young star cluster (YSC) population of this interacting system. Our newly extracted source catalogue contains 2834 cluster candidates, morphologically classified to be compact and uniform in colour, for which ages, masses and extinction are derived. In this first work we study the main properties of the YSC population of the whole galaxy, considering a mass-limited sample. Both luminosity and mass functions follow a power law shape with slope -2, but at high luminosities and masses a dearth of sources is observed. The analysis of the mass function suggests that it is best fitted by a Schechter function with slope -2 and a truncation mass at $1.00\pm0.12\times10^5\ M_\odot$. Through Monte Carlo simulations we confirm this result and link the shape of the luminosity function to the presence of a truncation in the mass function. A mass limited age function analysis, between 10 and 200 Myr, suggests that the cluster population is undergoing only moderate disruption. We observe little variation in the shape of the mass function at masses above $1\times10^4\ M_\odot$, over this age range. The fraction of star formation happening in the form of bound clusters in M51 is $\sim20\%$ in the age range 10 to 100 Myr and little variation is observed over the whole range from 1 to 200 Myr.

We present the results of an extensive observational campaign on the nearby Type Ibn SN 2015G, including data from radio through ultraviolet wavelengths. SN 2015G was asymmetric, showing late-time nebular lines redshifted by ~1000 km/s. It shared many features with the prototypical SN In 2006jc, including extremely strong He I emssion lines and a late-time blue pseudocontinuum. The young SN 2015G showed narrow P-Cygni profiles of He I, but never in its evolution did it show any signature of hydrogen - arguing for a dense, ionized, and hydrogen-free circumstellar medium moving outward with a velocity of ~1000 km/s and created by relatively recent mass loss from the progenitor star. Ultraviolet through infrared observations show that the fading SN 2015G (which was probably discovered some 20 days post-peak) had a spectral energy distribution that was well described by a simple, single-component blackbody. Archival HST images provide upper limits on the luminosity of SN 2015G's progenitor, while nondetections of any luminous radio afterglow and optical nondetections of outbursts over the past two decades provide constraints upon its mass-loss history.

A search for the progenitor of SN~2010jl, an unusually luminous core-collapse supernova of Type~IIn, using pre-explosion {\it Hubble}/WFPC2 and {\it Spitzer}/IRAC images of the region, yielded upper limits on the UV and near-infrared (IR) fluxes from any candidate star. These upper limits constrain the luminosity and effective temperature of the progenitor, the mass of any preexisting dust in its surrounding circumstellar medium (CSM), and dust proximity to the star. A {\it lower} limit on the CSM dust mass is required to hide a luminous progenitor from detection by {\it Hubble}. {\it Upper} limits on the CSM dust mass and constraints on its proximity to the star are set by requiring that the absorbed and reradiated IR emission not exceed the IRAC upper limits. Using the combined extinction-IR emission constraints we present viable $M_d-R_1$ combinations, where $M_d$ and $R_1$ are the CSM dust mass and its inner radius. These depend on the CSM outer radius, dust composition and grain size, and the properties of the progenitor. The results constrain the pre-supernova evolution of the progenitor, and the nature and origin of the observed post-explosion IR emission from SN~2010jl. In particular, an $\eta$~Car-type progenitor will require at least 4~mag of visual extinction to avoid detection by the {\it Hubble}. This can be achieved with dust masses $\gtrsim 10^{-3}$~\msun\ (less than the estimated 0.2-0.5~\msun\ around $\eta$~Car) which must be located at distances of $\gtrsim 10^{16}$~cm from the star to avoid detection by {\it Spitzer}.

We present optical and NIR photometry and spectroscopy of SN 2013L for the first four years post-explosion. SN 2013L was a moderately luminous (M$_{r}$ = -19.0) Type IIn supernova (SN) that showed signs of strong shock interaction with the circumstellar medium (CSM). The CSM interaction was equal to or stronger to SN 1988Z for the first 200 days and is observed at all epochs after explosion. Optical spectra revealed multi-component hydrogen lines appearing by day 33 and persisting and slowly evolving over the next few years. By day 1509 the H$\alpha$ emission was still strong and exhibiting multiple peaks, hinting that the CSM was in a disc or torus around the SN. SN 2013L is part of a growing subset of SNe IIn that shows both strong CSM interaction signatures and the underlying broad lines from the SN ejecta photosphere. The presence of a blue H$\alpha$ emission bump and a lack of a red peak does not appear to be due to dust obscuration since an identical profile is seen in Pa$\beta$. Instead this suggests a high concentration of material on the near-side of the SN or a disc inclination of roughly edge-on and hints that SN 2013L was part of a massive interactive binary system. Narrow H$\alpha$ P-Cygni lines that persist through the entirety of the observations measure a progenitor outflow speed of 80--130 km s$^{-1}$, speeds normally associated with extreme red supergiants, yellow hypergiants, or luminous blue variable winds. This progenitor scenario is also consistent with an inferred progenitor mass-loss rate of 0.3 - 8.0 $\times$ 10$^{-3}$ M$_{sun}$ yr$^{-1}$.

We report the large effort which is producing comprehensive high-level young star cluster (YSC) catalogues for a significant fraction of galaxies observed with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We present the methodology developed to extract cluster positions, verify their genuine nature, produce multiband photometry (from NUV to NIR), and derive their physical properties via spectral energy distribution fitting analyses. We use the nearby spiral galaxy NGC628 as a test case for demonstrating the impact that LEGUS will have on our understanding of the formation and evolution of YSCs and compact stellar associations within their host galaxy. Our analysis of the cluster luminosity function from the UV to the NIR finds a steepening at the bright end and at all wavelengths suggesting a dearth of luminous clusters. The cluster mass function of NGC628 is consistent with a power-law distribution of slopes $\sim -2$ and a truncation of a few times $10^5$ M$_\odot$. After their formation YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer timeframe, confirming their being potentially bound systems. Associations disappear on time scales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find mass-independent cluster disruption in the inner region of NGC628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass ($\leq$ $10^4$ M$_\odot$) clusters suggesting that a mass-dependent component is necessary to fully describe the YSC disruption process in NGC628.

We mine the Tycho-{\it Gaia} astrometric solution (TGAS) catalog for wide stellar binaries by matching positions, proper motions, and astrometric parallaxes. We separate genuine binaries from unassociated stellar pairs through a Bayesian formulation that includes correlated uncertainties in the proper motions and parallaxes. Rather than relying on assumptions about the structure of the Galaxy, we calculate Bayesian priors and likelihoods based on the nature of Keplerian orbits and the TGAS catalog itself. We calibrate our method using radial velocity measurements and obtain 6196 high-confidence candidate wide binaries with projected separations $s\lesssim1$ pc. The normalization of this distribution suggests that at least 0.6\% of TGAS stars have an associated, distant TGAS companion in a wide binary. We demonstrate that {\it Gaia}'s astrometry is precise enough that it can detect projected orbital velocities in wide binaries with orbital periods as large as 10$^6$ yr. For pairs with $s\ \lesssim\ 4\times10^4$~AU, characterization of random alignments indicate our contamination to be $\approx$5\%. For $s \lesssim 5\times10^3$~AU, our distribution is consistent with \"{O}pik's Law. At larger separations, the distribution is steeper and consistent with a power-law $P(s)\propto s^{-1.6}$; there is no evidence in our data of any bimodality in this distribution for $s \lesssim$ 1 pc. Using radial velocities, we demonstrate that at large separations, i.e., of order $s \sim$ 1 pc and beyond, any potential sample of genuine wide binaries in TGAS cannot be easily distinguished from ionized former wide binaries, moving groups, or contamination from randomly aligned stars.

SN2005ip was a TypeIIn event notable for its sustained strong interaction with circumstellar material (CSM), coronal emission lines, and IR excess, interpreted as shock interaction with the very dense and clumpy wind of an extreme red supergiant. We present a series of late-time spectra of SN2005ip and a first radio detection of this SN, plus late-time X-rays, all of which indicate that its CSM interaction is still strong a decade post-explosion. We also present and discuss new spectra of geriatric SNe with continued CSM interaction: SN1988Z, SN1993J, and SN1998S. From 3-10 yr post-explosion, SN2005ip's H-alpha luminosity and other observed characteristics were nearly identical to those of the radio-luminous SN1988Z, and much more luminous than SNe1993J and 1998S. At 10 yr after explosion, SN2005ip showed a drop in H$\alpha$ luminosity, followed by a quick resurgence over several months. We interpret this variability as ejecta crashing into a dense shell located at around 0.05 pc from the star, which may be the same shell that caused the IR echo at earlier epochs. The extreme H-alpha luminosities in SN2005ip and SN1988Z are still dominated by the forward shock at 10 yr post-explosion, whereas SN1993J and SN1998S are dominated by the reverse shock at a similar age. Continuous strong CSM interaction in SNe~2005ip and 1988Z is indicative of enhanced mass loss for about 1e3 yr before core collapse, longer than Ne, O, or Si burning phases. Instead, the episodic mass loss must extend back through C burning and perhaps even part of He burning.

Filaments in Herschel molecular cloud images are found to exhibit a "characteristic width". This finding is in tension with spatial power spectra of the data, which show no indication of this characteristic scale. We demonstrate that this discrepancy is a result of the methodology adopted for measuring filament widths. First, we perform the previously used analysis technique on artificial scale-free data, and obtain a peaked width distribution of filament-like structures. Next, we repeat the analysis on three Herschel maps and reproduce the narrow distribution of widths found in previous studies $-$ when considering the average width of each filament. However, the distribution of widths measured at all points along a filament spine is broader than the distribution of mean filament widths, indicating that the narrow spread (interpreted as a "characteristic" width) results from averaging. Furthermore, the width is found to vary significantly from one end of a filament to the other. Therefore, the previously identified peak at 0.1 pc cannot be understood as representing the typical width of filaments. We find an alternative explanation by modelling the observed width distribution as a truncated power-law distribution, sampled with uncertainties. The position of the peak is connected to the lower truncation scale and is likely set by the choice of parameters used in measuring filament widths. We conclude that a "characteristic" width of filaments is not supported by the available data.

The nature of the progenitor star (or system) for the Type IIn supernova (SN) subclass remains uncertain. While there are direct imaging constraints on the progenitors of at least four Type IIn supernovae, one of them being SN 2010jl, ambiguities remain in the interpretation of the unstable progenitors and the explosive events themselves. A blue source in pre-explosion HST/WFPC2 images falls within the 5 sigma astrometric error circle derived from post-explosion ground-based imaging of SN 2010jl. At the time the ground-based astrometry was published, however, the SN had not faded sufficiently for post-explosion HST follow-up observations to determine a more precise astrometric solution and/or confirm if the pre-explosion source had disappeared, both of which are necessary to ultimately disentangle the possible progenitor scenarios. Here we present HST/WFC3 imaging of the SN 2010jl field obtained in 2014 and 2015, when the SN had faded sufficiently to allow for new constraints on the progenitor. The SN, which is still detected in the new images, is offset by 0."099 +/- 0."008 (24 +/- 2 pc) from the underlying and extended source of emission that contributes at least partially, if not entirely, to the blue source previously suggested as the candidate progenitor in the WFPC2 data. This point alone rules out the possibility that the blue source in the pre-explosion images is the exploding star, but may instead suggest an association with a young (<5-6 Myr) cluster and still argues for a massive (>30 solar masses) progenitor. We obtain new upper limits on the flux from a single star at the SN position in the pre-explosion WFPC2 and Spitzer/IRAC images that may ultimately be used to constrain the progenitor properties.

SPIRITS---SPitzer InfraRed Intensive Transients Survey---is an ongoing survey of nearby galaxies searching for infrared (IR) transients with Spitzer/IRAC. We present the discovery and follow-up observations of one of our most luminous ($M_{[4.5]} = -17.1\pm0.4$ mag, Vega) and red ($[3.6] - [4.5] = 3.0 \pm 0.2$ mag) transients, SPIRITS 15c. The transient was detected in a dusty spiral arm of IC 2163 ($D\approx35.5$ Mpc). Pre-discovery ground-based imaging revealed an associated, shorter-duration transient in the optical and near-IR (NIR). NIR spectroscopy showed a broad ($\approx 8400$ km s$^{-1}$), double-peaked emission line of He I at $1.083 \mu$m, indicating an explosive origin. The NIR spectrum of SPIRITS 15c is similar to that of the Type IIb SN 2011dh at a phase of $\approx 200$ days. Assuming $A_V = 2.2$ mag of extinction in SPIRITS 15c provides a good match between their optical light curves. The IR light curves and the extreme $[3.6]-[4.5]$ color cannot be explained using only a standard extinction law. Another luminous ($M_{4.5} = -16.1\pm0.4$ mag) event, SPIRITS 14buu, was serendipitously discovered in the same galaxy. The source displays an optical plateau lasting $\gtrsim 80$ days, and we suggest a scenario similar to the low-luminosity Type IIP SN 2005cs obscured by $A_V \approx 1.5$ mag. Other classes of IR-luminous transients can likely be ruled out in both cases. If both events are indeed SNe, this may suggest $\gtrsim 18\%$ of nearby core-collapse SNe are missed by currently operating optical surveys.

We identify 709 arc-shaped mid-infrared nebula in 24 micron Spitzer Space Telescope or 22 micron Wide Field Infrared Explorer surveys of the Galactic Plane as probable dusty interstellar bowshocks powered by early-type stars. About 20% are visible at 8 microns or shorter mid-infrared wavelengths as well. The vast majority (660) have no previous identification in the literature. These extended infrared sources are strongly concentrated near Galactic mid-Plane with an angular scale height of ~0.6 degrees. All host a symmetrically placed star implicated as the source of a stellar wind sweeping up interstellar material. These are candidate "runaway" stars potentially having high velocities in the reference frame of the local medium. Among the 286 objects with measured proper motions, we find an unambiguous excess having velocity vectors aligned with the infrared morphology --- kinematic evidence that many of these are "runaway" stars with large peculiar motions responsible for the bowshock signature. We discuss a population of "in-situ" bowshocks (103 objects) that face giant HII regions where the relative motions between the star and ISM may be caused by bulk outflows from an overpressured bubble. We also identify 58 objects that face 8 micron bright-rimmed clouds and apparently constitute a sub-class of in-situ bowshocks where the stellar wind interacts with a photo-evaporative flow from an eroding molecular cloud interface (i.e., "PEF bowshocks"). Orientations of the arcuate nebulae exhibit a correlation over small angular scales, indicating that external influences such as HII regions are responsible for producing some bowshock nebulae. However, the vast majority of this sample appear to be isolated (499 objects) from obvious external influences.

We present late-time HST images of the site of supernova (SN) 2009ip taken almost 3 yr after its bright 2012 luminosity peak. SN 2009ip is now slightly fainter in broad filters than the progenitor candidate detected by HST in 1999. The current source continues to be dominated by ongoing late-time CSM interaction that produces strong H-alpha emission and a weak pseudo-continuum, as found previously for 1-2 yr after explosion. The intent of these observations was to search for evidence of recent star formation in the local (1kpc; 10 arcsec) environment around SN 2009ip, in the remote outskirts of its host spiral galaxy NGC 7259. We can rule out the presence of any massive star-forming complexes like 30 Dor or the Carina Nebula at the SN site or within a few kpc. If the progenitor of SN 2009ip was really a 50-80 Msun star as archival HST images suggested, then it is strange that there is no sign of this type of massive star formation anywhere in the vicinity. A possible explanation is that the progenitor was the product of a merger or binary mass transfer, rejuvenated after a lifetime that was much longer than 4-5 Myr, allowing its natal H II region to have faded. A smaller region like the Orion Nebula would be an unresolved but easily detected point source. This is ruled out within 1.5 kpc around SN 2009ip, but a small H II region could be hiding in the glare of SN 2009ip itself. Later images after a few more years have passed are needed to confirm that the progenitor candidate is truly gone and to test for the presence of a small H II region or cluster at the SN position.

For sufficiently wide orbital separations {\it a}, the two members of a stellar binary evolve independently. This implies that in a wide double white dwarf (DWD), the more massive WD should always be produced first, when its more massive progenitor ends its main-sequence life, and should therefore be older and cooler than its companion. The bound, wide DWD HS 2220$+$2146 ($a\approx500$ AU) does not conform to this picture: the more massive WD is the younger, hotter of the pair. We show that this discrepancy is unlikely to be due to past mass-transfer phases or to the presence of an unresolved companion. Instead, we propose that HS 2220$+$2146 formed through a new wide DWD evolutionary channel involving the merger of the inner binary in a hierarchical triple system. The resulting blue straggler and its wide companion then evolved independently, forming the WD pair seen today. Although we cannot rule out other scenarios, the most likely formation channel has the inner binary merging while both stars are still on the main sequence. This provides us with the tantalizing possibility that Kozai-Lidov oscillations may have played a role in the inner binary's merger. {\it Gaia} may uncover hundreds more wide DWDs, leading to the identification of other systems like HS 2220$+$2146. There are already indications that other WD systems may have been formed through different, but related, hierarchical triple evolutionary scenarios. Characterizing these populations may allow for thorough testing of the efficiency with which KL oscillations induce stellar mergers.

We present optical/near-infrared spectroscopy and photometry of the supernova (SN) PS15si. This object was originally identified as a Type IIn SN, but here we argue that it should be reclassified as a Type Ia SN with narrow hydrogen lines originating from interaction with circumstellar matter (CSM; i.e., SN Ia/IIn or SN Ia-CSM). Based on deep nondetections 27 days before discovery, we infer that this SN was discovered around or slightly before optical maximum, and we estimate the approximate time that it reached R-band maximum based on comparison with other SNe Ia/IIn. In terms of spectral morphology, we find that PS15si can be matched to a range of SN Ia spectral types, although SN 1991T-like SNe Ia provides the most self-consistent match. While this observation agrees with analysis of most other SNe Ia/IIn, we find that the implied CSM luminosity is too low to account for the overall luminosity of the SN at a time when the CSM should outshine the underlying SN by a few magnitudes. We infer that the similarity between PS15si and the hot, overluminous, high-ionisation spectrum of SN 1991T is a consequence of a spectrum that originates in ejecta layers that are heated by ultraviolet/X-ray radiation from CSM interaction. In addition, PS15si may have rebrightened over a short timescale in the B and V bands around 85 days after discovery, perhaps indicating that the SN ejecta are interacting with a local enhancement in CSM produced by clumps or a shell at large radii.

The Galactic transient V1309 Sco was the result of a merger in a low-mass star system, while V838 Mon was thought to be a similar merger event from a more massive B-type progenitor. In this paper we study an optical/IR transient discovered in the nearby galaxy NGC4490, which appeared similar to these merger events (unobscured progenitor, irregular multi-peaked light curve, increasingly red color, similar optical spectrum, IR excess at late times), but which had a higher peak luminosity and longer duration in outburst. NGC4490-OT has less in common with the class of SN~2008S-like transients. A progenitor detected in pre-eruption HST images, combined with upper limits in the IR, requires a luminous and blue progenitor that has faded in late-time HST images. The same source was detected by Spitzer and ground-based data as a luminous IR transient, indicating a transition to a self-obscured state qualitatively similar to the evolution seen in other stellar mergers and in LBVs. The post-outburst dust-obscured source is too luminous and too warm at late times to be explained with an IR echo, suggesting that the object survived the event. The luminosity of the enshrouded IR source is similar to that of the progenitor. Compared to proposed merger events, the more massive progenitor of NGC4490-OT seems to extend a correlation between stellar mass and peak luminosity, and may suggest that both of these correlate with duration. We show that spectra of NGC4490-OT and V838 Mon also resemble light-echo spectra of eta Car, prompting us to speculate that eta Car may be an extreme extension of this phenomenon.

While SN impostors resemble the Great Eruption of eta Car in the sense that their spectra show narrow H lines and they have typical peak absolute magnitudes of -13 to -14 mag, most extragalactic events observed so far are quite different from eta Car in duration. Their bright phases typically last for 100~d or less, rather than persisting for several years. The transient object UGC2773-OT had a similar peak absolute magnitude to other SN impostors, but with a gradual 5-yr prediscovery rise. In the 6 yr since discovery, it has faded very slowly (0.26 mag/yr). Overall, we suggest that its decade-long eruption is so far the best known analog of eta Car's 19th century eruption. We discuss extensive spectroscopy of the ongoing eruption. The spectra show interesting changes in velocity and line shape that we discuss in detail, including an asymmetric Halpha emission line that we show is consistent with the ejection of a bipolar nebula that could be very much like the Homunculus of eta Car. Moreover, changes in the line width, line profile, blue excess emission resembling that of Type IIn supernovae, and the intensity of Halpha suggest the presence of strong circumstellar interaction in the eruption at late times. This supports the hypothesis that the extended plateau of eta Car's eruption may have been powered by shock interaction as well. One interesting difference compared to eta Car, however, is that UGC2773-OT so far does not exhibit the repeated brief spikes in luminosity that have been associated with binary periastron events.

Supernovae Type Iax (SNe Iax) are less energetic and less luminous than typical thermonuclear explosions. A suggested explanation for the observed characteristics of this subclass is a binary progenitor system consisting of a CO white dwarf primary accreting from a helium star companion. A single-degenerate explosion channel might be expected to result in a dense circumstellar medium (CSM), although no evidence for such a CSM has yet been observed for this subclass. Here we present recent Spitzer observations of the SN Iax 2014dt obtained by the SPIRITS program nearly one year post-explosion that reveal a strong mid-IR excess over the expected fluxes of more normal SNe Ia. This excess is consistent with 1E-5 M_solar of newly formed dust, which would be the first time that newly formed dust has been observed to form in a normal Type Ia. The excess, however, is also consistent with a dusty CSM that was likely formed in pre-explosion mass-loss, thereby suggesting a single degenerate progenitor system. Compared to other SNe Ia that show significant shock interaction (SNe Ia-CSM) and interacting core-collapse events (SNe IIn), this dust shell in SN 2014dt is less massive. We consider the implications that such a pre-existing dust shell has for the progenitor system, including a binary system with a mass donor that is a red giant, a red supergiant, and an asymptotic giant branch star.

The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly extended, with trace hydrocarbons, a global haze layer, and a surface pressure near 10 microbars. Pluto's diverse surface geology and long-term activity raise fundamental questions about how small planets remain active many billions of years after formation. Pluto's large moon Charon displays tectonics and evidence for a heterogeneous crustal composition, its north pole displays puzzling dark terrain. Small satellites Hydra and Nix have higher albedos than expected.

We present observational constraints on the initial-final mass relation (IFMR) using wide double white dwarfs (DWDs). We identify 65 new candidate wide DWDs within the Sloan Digital Sky Survey, bringing the number of candidate wide DWDs to 142. We then engage in a spectroscopic follow-up campaign and collect existing spectra for these objects; using these spectra, we derive masses and cooling ages for 54 hydrogen (DA) WDs in DWDs. We also identify one new DA/DB pair, four candidate DA/DC pairs, four candidate DA/DAH pairs, and one new candidate triple degenerate system. Because wide DWDs are co-eval and evolve independently, the difference in the pre-WD lifetimes should equal the difference in the WD cooling ages. We use this to develop a Bayesian hierarchical framework and construct a likelihood function to determine the probability that any particular IFMR fits a sample of wide DWDs. We then define a parametric model for the IFMR and find the best parameters indicated by our sample of DWDs. We place robust constraints on the IFMR for initial masses of 2--4 \Msun. The WD masses produced by our model for stars within this mass range differ from those predicted by semi-empirical fits to open cluster WDs. Within this mass range, where there are few constraining open cluster WDs and disagreements in the cluster ages, wide DWDs may provide more reliable constraints on the IFMR. Expanding this method to the many wide DWDs expected to be discovered by Gaia may transform our understanding of the IFMR.

SN 2006gy was the most luminous SN ever observed at the time of its discovery and the first of the newly defined class of superluminous supernovae (SLSNe). The extraordinary energetics of SN 2006gy and all SLSNe (>10^51 erg) require either atypically large explosion energies (e.g., pair-instability explosion) or the efficient conversion of kinetic into radiative energy (e.g., shock interaction). The mass-loss characteristics can therefore offer important clues regarding the progenitor system. For the case of SN 2006gy, both a scattered and thermal light echo from circumstellar material (CSM) have been reported at later epochs (day ~800), ruling out the likelihood of a pair-instability event and leading to constraints on the characteristics of the CSM. Owing to the proximity of the SN to the bright host-galaxy nucleus, continued monitoring of the light echo has not been trivial, requiring the high resolution offered by the Hubble Space Telescope (HST) or ground-based adaptive optics (AO). Here we report detections of SN 2006gy using HST and Keck AO at ~3000 days post-explosion and consider the emission mechanism for the very late-time light curve. While the optical light curve and optical spectral energy distribution are consistent with a continued scattered-light echo, a thermal echo is insufficient to power the K'-band emission by day 3000. Instead, we present evidence for late-time infrared emission from dust that is radiatively heated by CSM interaction within an extremely dense dust shell, and we consider the implications on the CSM characteristics and progenitor system.

SN 2011ja was a bright (I = -18.3) Type II supernova occurring in the nearby edge on spiral galaxy NGC 4945. Flat-topped and multi-peaked H-alpha and H-beta spectral emission lines appear between 64 - 84 days post-explosion, indicating interaction with a disc-like circumstellar medium inclined 30-45 degrees from edge-on. After day 84 an increase in the H- and K-band flux along with heavy attenuation of the red wing of the emission lines are strong indications of early dust formation, likely located in the cool dense shell created between the forward shock of the SN ejecta and the reverse shock created as the ejecta plows into the existing CSM. Radiative transfer modeling reveals both ~1.5 x 10^-4 Msun of pre-existing dust located ~ 10^16.7 cm away and ~ 5 x 10^-5 Msun of newly formed dust. Spectral observations after 1.5 years reveal the possibility that the fading SN is located within a young (3-6 Myr) massive stellar cluster, which when combined with tentative 56Ni mass estimates of 0.2 Msun may indicate a massive (> 25 Msun) progenitor for SN 2011ja.

The nearby dwarf starburst galaxy NGC5253 hosts a number of young, massive star clusters, the two youngest of which are centrally concentrated and surrounded by thermal radio emission (the `radio nebula'). To investigate the role of these clusters in the starburst energetics, we combine new and archival Hubble Space Telescope images of NGC5253 with wavelength coverage from 1500 Ang to 1.9 micron in 13 filters. These include H-alpha, P-beta, and P-alpha, and the imaging from the Hubble Treasury Program LEGUS (Legacy Extragalactic UV Survey). The extraordinarily well-sampled spectral energy distributions enable modeling with unprecedented accuracy the ages, masses, and extinctions of the 9 optically brightest clusters (M_V < -8.8) and the two young radio nebula clusters. The clusters have ages ~1-15 Myr and masses ~1x10^4 - 2.5x10^5 M_sun. The clusters' spatial location and ages indicate that star formation has become more concentrated towards the radio nebula over the last ~15 Myr. The most massive cluster is in the radio nebula; with a mass 2.5x10^5 M_sun and an age ~1 Myr, it is 2-4 times less massive and younger than previously estimated. It is within a dust cloud with A_V~50 mag, and shows a clear nearIR excess, likely from hot dust. The second radio nebula cluster is also ~1 Myr old, confirming the extreme youth of the starburst region. These two clusters account for about half of the ionizing photon rate in the radio nebula, and will eventually supply about 2/3 of the mechanical energy in present-day shocks. Additional sources are required to supply the remaining ionizing radiation, and may include very massive stars.

Supernova (SN) 2002hh was unusual among core-collapse SNe because it was highly reddened, and displayed a bright infrared (IR) excess due to radiatively heated dust in its circumstellar medium (CSM). Estimates for the mass of dust responsible for the IR echo suggested the presence of a massive shell within 0.26 pc of the star. For a velocity of 5000 - 10000 km/s, this material should be hit by the SN blast wave at late times, starting at roughly 12 years post-explosion. We have obtained deep late-time spectra with the MMT Blue Channel spectrograph to search for any spectral signatures of ongoing shock interaction. Interaction with a strength comparable to SN 1987A's collision with the equatorial ring would be detected in our data. However, in the spectra reported here, we do not detect clear signs of strong CSM interaction, contrary to expectations based on the reported radii of the dust shell. We do, however, detect emission associated with the old SN, and we find that the broad lines in the spectrum indicate a continuation of an ongoing reflected light echo, which appears similar to the spectrum at peak luminosity for this Type II-P event.

We have discovered a luminous light echo around the normal Type II-Plateau Supernova (SN) 2012aw in Messier 95 (M95; NGC 3351), detected in images obtained approximately two years after explosion with the Wide Field Channel 3 on-board the Hubble Space Telescope (HST) by the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). The multi-band observations span from the near-ultraviolet through the optical (F275W, F336W, F438W, F555W, and F814W). The apparent brightness of the echo at the time was ~21--22 mag in all of these bands. The echo appears circular, although less obviously as a ring, with an inhomogeneous surface brightness, in particular, a prominent enhanced brightness to the southeast. The SN itself was still detectable, particularly in the redder bands. We are able to model the light echo as the time-integrated SN light scattered off of diffuse interstellar dust in the SN environment. We have assumed that this dust is analogous to that in the Milky Way with R_V=3.1. The SN light curves that we consider also include models of the unobserved early burst of light from the SN shock breakout. Our analysis of the echo suggests that the distance from the SN to the scattering dust elements along the echo is ~45 pc. The implied visual extinction for the echo-producing dust is consistent with estimates made previously from the SN itself. Finally, our estimate of the SN brightness in F814W is fainter than that measured for the red supergiant star at the precise SN location in pre-SN images, possibly indicating that the star has vanished and confirming it as the likely SN progenitor.

We present parallax observations and a detailed model atmosphere analysis of 54 cool and ultracool ($T_{\rm eff}$ < 4000 K) white dwarfs (WDs) in the solar neighbourhood. For the first time, a large number of cool and ultracool WDs have distance and tangential velocities measurements available. Our targets have distances ranging from 21 pc to >100 pc, and include five stars within 30 pc. Contrary to expectations, all but two of them have tangential velocities smaller than 150 km s$^{-1}$ thus suggesting Galactic disc membership. The oldest WDs in this sample have WD cooling ages of 10 Gyr, providing a firm lower limit to the age of the thick disc population. Many of our targets have uncharacteristically large radii, indicating that they are low mass WDs. It appears that we have detected the brighter population of cool and ultracool WDs near the Sun. The fainter population of ultracool CO-core WDs remain to be discovered in large numbers. The Large Synoptic Survey Telescope should find these elusive, more massive ultracool WDs in the solar neighbourhood.

We analyze binary population models of double-neutron stars and compare results to the accurately measured orbital periods and eccentricities of the eight known such systems in our Galaxy. In contrast to past similar studies, we especially focus on the dominant evolutionary channels (we identify three); for the first time, we use a detailed understanding of the evolutionary history of three double neutron stars as actual constraints on the population models. We find that the evolutionary constraints derived from the double pulsar are particularly tight, and less than half of the examined models survive the full set of constraints. The top-likelihood surviving models yield constraints on the key binary evolution parameters, but most interestingly reveal (i) the need for electron-capture supernovae from relatively low-mass degenerate, progenitor cores, and (ii) the most likely evolutionary paths for the rest of the known double neutron stars. In particular, we find that J1913+16 likely went through a phase of Case BB mass transfer, and J1906+0746 and J1756-2251 are consistent with having been formed in electron-capture supernovae.

PTF11iqb was initially classified as a TypeIIn event caught very early after explosion. It showed narrow Wolf-Rayet (WR) spectral features on day 2, but the narrow emission weakened quickly and the spectrum morphed to resemble those of Types II-L and II-P. At late times, Halpha emission exhibited a complex, multipeaked profile reminiscent of SN1998S. In terms of spectroscopic evolution, we find that PTF11iqb was a near twin of SN~1998S, although with weaker interaction with circumstellar material (CSM) at early times, and stronger CSM interaction at late times. We interpret the spectral changes as caused by early interaction with asymmetric CSM that is quickly (by day 20) enveloped by the expanding SN ejecta photosphere, but then revealed again after the end of the plateau when the photosphere recedes. The light curve can be matched with a simple model for weak CSM interaction added to the light curve of a normal SN~II-P. This plateau requires that the progenitor had an extended H envelope like a red supergiant, consistent with the slow progenitor wind speed indicated by narrow emission. The cool supergiant progenitor is significant because PTF11iqb showed WR features in its early spectrum --- meaning that the presence of such WR features in an early SN spectrum does not necessarily indicate a WR-like progenitor. [abridged] Overall, PTF11iqb bridges SNe~IIn with weaker pre-SN mass loss seen in SNe II-L and II-P, implying a continuum between these types.

Measuring the masses of companions to single-line spectroscopic binary stars is (in general) not possible because of the unknown orbital plane inclination. Even when the mass of the visible star can be measured, only a lower limit can be placed on the mass of the unseen companion. However, since these inclination angles should be isotropically distributed, for a large enough, unbiased sample, the companion mass distribution can be deconvolved from the distribution of observables. In this work, we construct a hierarchical probabilistic model to infer properties of unseen companion stars given observations of the orbital period and projected radial velocity of the primary star. We apply this model to three mock samples of low-mass white dwarfs (LMWDs, $M\lesssim0.45~\Msun$) and a sample of post-common-envelope binaries. We use a mixture of two Gaussians to model the WD and neutron star (NS) companion mass distributions. Our model successfully recovers the initial parameters of these test data sets. We then apply our model to 55 WDs in the extremely low-mass (ELM) WD Survey. Our maximum a posteriori model for the WD companion population has a mean mass $\mu_{\rm WD} = 0.74~\Msun$, with a standard deviation $\sigma_{\rm WD} = 0.24~\Msun$. Our model constrains the NS companion fraction $f_{\rm NS}$ to be $<$16\% at 68\% confidence. We make samples from the posterior distribution publicly available so that future observational efforts may compute the NS probability for newly discovered LMWDs.

The Legacy ExtraGalactic UV Survey (LEGUS) is a Cycle 21 Treasury program on the Hubble Space Telescope, aimed at the investigation of star formation and its relation with galactic environment in nearby galaxies, from the scales of individual stars to those of ~kpc-size clustered structures. Five-band imaging, from the near-ultraviolet to the I-band, with the Wide Field Camera 3, plus parallel optical imaging with the Advanced Camera for Surveys, is being collected for selected pointings of 50 galaxies within the local 12 Mpc. The filters used for the observations with the Wide Field Camera 3 are: F275W(2,704 A), F336W(3,355 A), F438W(4,325 A), F555W(5,308 A), and F814W(8,024 A); the parallel observations with the Advanced Camera for Surveys use the filters: F435W(4,328 A), F606W(5,921 A), and F814W(8,057 A). The multi-band images are yielding accurate recent (<~50 Myr) star formation histories from resolved massive stars and the extinction-corrected ages and masses of star clusters and associations. The extensive inventories of massive stars and clustered systems will be used to investigate the spatial and temporal evolution of star formation within galaxies. This will, in turn, inform theories of galaxy evolution and improve the understanding of the physical underpinning of the gas-star formation relation and the nature of star formation at high redshift. This paper describes the survey, its goals and observational strategy, and the initial science results. Because LEGUS will provide a reference survey and a foundation for future observations with JWST and with ALMA, a large number of data products are planned for delivery to the community.

In wide double white dwarf (DWD) binaries, in which the co-eval WDs evolve independently, the more massive, faster-evolving WD can be used to obtain a main-sequence lifetime for the less-massive WD. By converting this lifetime into an initial mass for the less-massive WD, and combining it with the spectroscopically derived mass for this WD, one can constrain the initial-final mass relation (IFMR). However, the number of known and well-characterized DWDs is small, severely limiting their usefulness for this analysis. To obtain new constraints on the IFMR, we search for wide DWDs in the Sloan Digital Sky Survey (SDSS) Data Release 9. We find 65 new candidate systems, thereby raising the number of known wide DWDs to 142. We then engage in a spectroscopic campaign to characterize these pairs, identifying 32 DA/DA pairs, two DA/DB pairs, four DA/DAH candidate pairs, a previously unidentified candidate triple WD system, and five DA/DC WDs. We present a reanalysis of the constraint on the IFMR placed by Finley & Koester (1997) using the DWD PG 0922+162, and finish by discussing how it could be expanded to a generic set of wide DWDs.

We have used multi-wavelength Hubble Space Telescope WFC3 data of the starbursting spiral galaxy M83 in order to measure variations in the upper end of the stellar initial mass function (uIMF) using the production rate of ionizing photons in unresolved clusters with ages $\leq$ 8 Myr. As in earlier papers on M51 and NGC 4214, the upper end of the stellar IMF in M83 is consistent with an universal IMF, and stochastic sampling of the stellar populations in the $\lessapprox$ 10$^{3}$ Msun clusters are responsible for any deviations in this universality. The ensemble cluster population, as well as individual clusters, also imply that the most massive star in a cluster does not depend on the cluster mass. In fact, we have found that these small clusters seem to have an over-abundance of ionizing photons when compared to an expected universal or truncated IMF. This also suggests that the presence of massive stars in these clusters does not affect the star formation in a destructive way.

Hierarchical structure in ultraviolet images of 12 late-type LEGUS galaxies is studied by determining the numbers and fluxes of nested regions as a function of size from ~1 to ~200 pc, and the number as a function of flux. Two starburst dwarfs, NGC 1705 and NGC 5253, have steeper number-size and flux-size distributions than the others, indicating high fractions of the projected areas filled with star formation. Nine subregions in 7 galaxies have similarly steep number-size slopes, even when the whole galaxies have shallower slopes. The results suggest that hierarchically structured star-forming regions several hundred parsecs or larger represent common unit structures. Small galaxies dominated by only a few of these units tend to be starbursts. The self-similarity of young stellar structures down to parsec scales suggests that star clusters form in the densest parts of a turbulent medium that also forms loose stellar groupings on larger scales. The presence of super star clusters in two of our starburst dwarfs would follow from the observed structure if cloud and stellar subregions more readily coalesce when self-gravity in the unit cell contributes more to the total gravitational potential.

We present an analysis of 98 galaxies of low-dust content, selected from the Spitzer Local Volume Legacy survey, aimed at examining the relation between the ultraviolet (UV) color, beta, and dust attenuation in normal star-forming galaxies. The IRX-beta diagram relates the total dust attenuation in a galaxy, traced by the far-IR (FIR) to UV ratio, to the observed UV color, indicated by beta. Previous research has indicated that while starburst galaxies exhibit a relatively tight IRX-beta relation, normal star-forming galaxies do not, and have a much larger spread in the total-IR to far-UV (FUV) luminosity for a fixed UV color. We examine the role that the age of the stellar population plays as the "second parameter" responsible for the observed deviation and spread of star-forming galaxies from the starburst relation. We model the FUV to FIR spectral energy distribution (SED) of each galaxy according to two broad bins of star formation history (SFH): constant and instantaneous burst. We find clear trends between stellar population mean age estimators (extinction-corrected FUV/NIR, U-B, and EW(H-alpha)) and the UV color, beta; the trends are mostly driven by the galaxies best-described by instantaneous burst populations. We also find a significant correlation between beta and the mean age directly determined from the best-fit instantaneous models. As already indicated by other authors, the UV attenuation in star-forming galaxies may not be recovered with the UV color alone and is highly influenced by the stellar population's mean age and SFH. Overall, the scatter in the IRX-beta diagram is better correlated with beta than with the perpendicular distance.

The production rate of ionizing photons in young < 8 Myr, unresolved stellar clusters in the nearby irregular galaxy NGC 4214 is probed using multi-wavelength Hubble Space Telescope WFC3 data. We normalize the ionizing photon rate by the cluster mass to investigate the upper end of the stellar initial mass function (IMF). We have found that within the uncertainties the upper end of the stellar IMF appears to be universal in this galaxy, and that deviations from a universal IMF can be attributed to stochastic sampling of stars in clusters with masses < 10^3 M_sun. Furthermore, we have found that there does not seem to be a dependence of the maximum stellar mass on the cluster mass. We have also found that for massive clusters, feedback may cause an underrepresentation in Ha luminosities, which needs to be taken into account when conducting this type of analysis.

Wide binaries made up of two white dwarfs (WDs) receive far less attention than their tight counterparts. However, our tests using the binary population synthesis code {\tt StarTrack} indicate that, for any set of reasonable initial conditions, there exists a significant observable population of double white dwarfs (WDWDs) with orbital separations of 10$^2$ to 10$^5$ AU. We adapt the technique of Dhital et al.\ to search for candidate common proper motion WD companions separated by $<10\amin$ around the $>$12,000 spectroscopically confirmed hydrogen-atmosphere WDs recently identified in the Sloan Digital Sky Survey. Using two techniques to separate random alignments from high-confidence pairs, we find nine new high-probability wide WDWDs and confirm three previously identified candidate wide WDWDs. This brings the number of known wide WDWDs to 45; our new pairs are a significant addition to the sample, especially at small proper motions ($<$200 mas/yr) and large angular separations ($>$10\asec). Spectroscopic follow-up and an extension of this method to a larger, photometrically selected set of SDSS WDs may eventually produce a large enough dataset for WDWDs to realize their full potential as testbeds for theories of stellar evolution.

We present nebular phase optical imaging and spectroscopy and near/mid-IR imaging of the Type II SN 2006bc. Observations reveal the central wavelength of the symmetric H$\alpha$ line profile to be red-shifted with respect to the host galaxy H$\alpha$ emission by day 325. Such an phenomenon has been argued to result from an asymmetric explosion in the iron-peak elements resulting in a larger mass of $^{56}$Ni and higher excitation of hydrogen on the far side of the SN explosion. We also observe a gradual blue-shifting of this H$\alpha$ peak which is indicative of dust formation in the ejecta. Although showing a normal peak brightness, V $\sim$ -17.2, for a core-collapse SN, 2006bc fades by $\sim$6 mag during the first 400 days suggesting either a relatively low $^{56}$Ni yield, an increase in extinction due to new dust, or both. A short duration flattening of the light curve is observed from day 416 to day 541 suggesting an optical light echo. Based on the narrow time window of this echo, we discuss implications on the location and geometry of the reflecting ISM. With our radiative transfer models, we find an upper limit of 2 x 10$^{-3}$ M$_{\odot}$ of dust around SN 2006bc. In the event that all of this dust were formed during the SN explosion, this quantity of dust is still several orders of magnitude lower than that needed to explain the large quantities of dust observed in the early universe.

We present a detailed model atmosphere analysis of two of the oldest stars known in the solar neighborhood, the high proper motion white dwarfs SDSS J110217.48+411315.4 (hereafter J1102) and WD 0346+246 (hereafter WD0346). We present trigonometric parallax observations of J1102, which places it at a distance of only 33.7 +- 2.0 pc. Based on the state of the art model atmospheres, optical, near-, mid-infrared photometry, and distances, we constrain the temperatures, atmospheric compositions, masses, and ages for both stars. J1102 is an 11 Gyr old (white dwarf plus main-sequence age), 0.62 Msol white dwarf with a pure H atmosphere and Teff = 3830 K. WD0346 is an 11.5 Gyr old, 0.77 Msol white dwarf with a mixed H/He atmosphere and Teff = 3650 K. Both stars display halo kinematics and their ages agree remarkably well with the ages of the nearest globular clusters, M4 and NGC 6397. J1102 and WD0346 are the closest examples of the oldest halo stars that we know of.

We report optical and mid-infrared photometry of SN 1980K between 2004 and 2010, which show slow monotonic fading consistent with previous spectroscopic and photometric observations made 8 to 17 years after outburst. The slow rate-of-change over two decades suggests that this evolution may result from scattered and thermal light echoes off of extended circumstellar material. We present a semi- analytic dust radiative-transfer model that uses an empirically corrected effective optical depth to provide a fast and robust alternative to full Monte-Carlo radiative transfer modeling for homogenous dust at low to intermediate optical depths. We find that unresolved echoes from a thin circumstellar shell 14-15 lt-yr from the progenitor, and containing about 0.02 Msun of carbon-rich dust, can explain the broadband spectral and temporal evolution. The size, mass and dust composition are in good agreement with the contact discontinuity observed in scattered echoes around SN 1987A. The origin of slowly-changing high-velocity [O I] and Halpha lines is also considered. We propose an origin in shocked high-velocity metal-rich clumps of ejecta, rather than arising in the impact of ejecta on slowly-moving circumstellar material, as is the case with hot spots in SN 1987A.

We present BVRIJHK band photometry of 6 core-collapse supernovae, SNe 1999bw, 2002hh, 2003gd, 2004et, 2005cs, and 2006bc measured at late epochs (>2 yrs) based on Hubble Space Telescope (HST), Gemini north, and WIYN telescopes. We also show the JHK lightcurves of a supernova impostor SN 2008S up to day 575. Of our 43 HST observations in total, 36 observations are successful in detecting the light from the SNe alone and measuring magnitudes of all the targets. HST observations show a resolved scattered light echo around SN 2003gd at day 1520 and around SN 2002hh at day 1717. Our Gemini and WIYN observations detected SNe 2002hh and 2004et, as well. Combining our data with previously published data, we show VRIJHK-band lightcurves and estimate decline magnitude rates at each band in 4 different phases. Our prior work on these lightcurves and other data indicate that dust is forming in our targets from day ~300-400, supporting SN dust formation theory. In this paper we focus on other physical properties derived from the late time light curves. We estimate 56Ni masses for our targets (0.5-14 x 10^{-2} Msun) from the bolometric lightcurve of each for days ~150-300 using SN 1987A as a standard (7.5 x 10^{-2} Msun). The flattening or sometimes increasing fluxes in the late time light curves of SNe 2002hh, 2003gd, 2004et and 2006bc indicate the presence of light echos. We estimate the circumstellar hydrogen density of the material causing the light echo and find that SN 2002hh is surrounded by relatively dense materials (n(H) >400 cm^{-3}) and SNe 2003gd and 2004et have densities more typical of the interstellar medium (~1 cm^{-3}). The 56Ni mass appears well correlated with progenitor mass with a slope of 0.31 x 10^{-2}, supporting the previous work by Maeda et al. (2010), who focus on more massive Type II SNe. The dust mass does not appear to be correlated with progenitor mass.

In 2007, R Coronae Borealis (R CrB) went into an historically deep and long decline. In this state, the dust acts like a natural coronagraph at visible wavelengths, allowing faint nebulosity around the star to be seen. Imaging has been obtained from 0.5 to 500 micron with Gemini/GMOS, HST/WFPC2, Spitzer/MIPS, and Herschel/SPIRE. Several of the structures around R CrB are cometary globules caused by wind from the star streaming past dense blobs. The estimated dust mass of the knots is consistent with their being responsible for the R CrB declines if they form along the line of sight to the star. In addition, there is a large diffuse shell extending up to 4 pc away from the star containing cool 25 K dust that is detected all the way out to 500 micron. The SED of R CrB can be well fit by a 150 AU disk surrounded by a very large diffuse envelope which corresponds to the size of the observed nebulosity. The total masses of the disk and envelope are 10^-4 and 2 M(Sun), respectively, assuming a gas-to-dust ratio of 100. The evidence pointing toward a white-dwarf merger or a final-helium-shell flash origin for R CrB is contradictory. The shell and the cometary knots are consistent with a fossil planetary nebula. Along with the fact that R CrB shows significant Lithium in its atmosphere, this supports the final-helium-shell flash. However, the relatively high inferred mass of R CrB and its high fluorine abundance support a white-dwarf merger.

We present an analysis of multi-epoch observations of the Type II-P supernova SN 2004et. New and archival optical spectra of SN 2004et are used to study the evolution of the Halpha and [O I] 6300A line profiles between days 259 and 646. Mid-infrared imaging was carried out between 2004 to 2010. We include Spitzer `warm' mission photometry at 3.6 and 4.5um obtained on days 1779, 1931 and 2151, along with ground-based and HST optical and near-infrared observations obtained between days 79 and 1803. Multi-wavelength light curves are presented, as well as optical-infrared spectral energy distributions (SEDs) for multiple epochs. Starting from about day 300, the optical light curves provide evidence for an increasing amount of circumstellar extinction attributable to newly formed dust, with the additional extinction reaching 0.8-1.5 magnitudes in the V-band by day 690. The overall SEDs were fitted with multiple blackbody components, in order to investigate the luminosity evolution of the supernova, and then with Monte Carlo radiative transfer models using smooth or clumpy dust distributions, in order to estimate how much new dust condensed in the ejecta. The luminosity evolution was consistent with the decay of 56Co in the ejecta up until about day 690, after which an additional emission source is required, in agreement with the findings of Kotak et al. (2009). Clumped dust density distributions consisting of 20% amorphous carbons and 80% silicates by mass were able to match the observed optical and infrared SEDs, with dust masses that increased from 8x10^{-5} Msun on day 300 to 1.5x10^{-3} Msun on day 690, still significantly lower than the values needed for core collapse supernovae to make a significant contribution to the dust enrichment of galaxies.

The Cosmic Origins Spectrograph (COS) is a moderate-resolution spectrograph with unprecedented sensitivity that was installed into the Hubble Space Telescope (HST) in May 2009, during HST Servicing Mission 4 (STS-125). We present the design philosophy and summarize the key characteristics of the instrument that will be of interest to potential observers. For faint targets, with flux F_lambda ~ 1.0E10-14 ergs/s/cm2/Angstrom, COS can achieve comparable signal to noise (when compared to STIS echelle modes) in 1-2% of the observing time. This has led to a significant increase in the total data volume and data quality available to the community. For example, in the first 20 months of science operation (September 2009 - June 2011) the cumulative redshift pathlength of extragalactic sight lines sampled by COS is 9 times that sampled at moderate resolution in 19 previous years of Hubble observations. COS programs have observed 214 distinct lines of sight suitable for study of the intergalactic medium as of June 2011. COS has measured, for the first time with high reliability, broad Lya absorbers and Ne VIII in the intergalactic medium, and observed the HeII reionization epoch along multiple sightlines. COS has detected the first CO emission and absorption in the UV spectra of low-mass circumstellar disks at the epoch of giant planet formation, and detected multiple ionization states of metals in extra-solar planetary atmospheres. In the coming years, COS will continue its census of intergalactic gas, probe galactic and cosmic structure, and explore physics in our solar system and Galaxy.

SN 2010jl was an extremely bright, Type IIn SNe which showed a significant IR excess no later than 90 days after explosion. We have obtained Spitzer 3.6 and 4.5 \mum and JHK observations of SN 2010jl \sim90 days post explosion. Little to no reddening in the host galaxy indicated that the circumstellar material lost from the progenitor must lie in a torus inclined out of the plane of the sky. The likely cause of the high mid-IR flux is the reprocessing of the initial flash of the SN by pre-existing circumstellar dust. Using a 3D Monte Carlo Radiative Transfer code, we have estimated that between 0.03-0.35 Msun of dust exists in a circumstellar torus around the SN located 6 \times 10 ^17 cm away from the SN and inclined between 60-80\cdot to the plane of the sky. On day 90, we are only seeing the illumination of approximately 5% of this torus, and expect to see an elevated IR flux from this material up until day \sim 450. It is likely this dust was created in an LBV-like mass loss event of more than 3 Msun, which is large but consistent with other LBV progenitors such as {\eta} Carinae.

We identify SDSS J010657.39-100003.3 (hereafter J0106-1000) as the shortest period detached binary white dwarf (WD) system currently known. We targeted J0106-1000 as part of our radial velocity program to search for companions around known extremely low-mass (ELM, ~ 0.2 Msol) WDs using the 6.5m MMT. We detect peak-to-peak radial velocity variations of 740 km/s with an orbital period of 39.1 min. The mass function and optical photometry rule out a main-sequence star companion. Follow-up high-speed photometric observations obtained at the McDonald 2.1m telescope reveal ellipsoidal variations from the distorted primary but no eclipses. This is the first example of a tidally distorted WD. Modeling the lightcurve, we constrain the inclination angle of the system to be 67 +- 13 deg. J0106-1000 contains a pair of WDs (0.17 Msol primary + 0.43 Msol invisible secondary) at a separation of 0.32 Rsol. The two WDs will merge in 37 Myr and most likely form a core He-burning single subdwarf star. J0106-1000 is the shortest timescale merger system currently known. The gravitational wave strain from J0106-1000 is at the detection limit of the Laser Interferometer Space Antenna (LISA). However, accurate ephemeris and orbital period measurements may enable LISA to detect J0106-1000 above the Galactic background noise.

SN 2007it is a bright, Type IIP supernova which shows indications of both pre-existing and newly formed dust. The visible photometry shows a bright late-time luminosity, powered by the 0.09 M$_{\sun}$ of $^{56}$Ni present in the ejecta. There is also a sudden drop in optical brightness after day 339, and a corresponding brightening in the IR due to new dust forming in the ejecta. CO and SiO emission, generally thought to be precursors to dust formation, may have been detected in the mid-IR photometry of SN 2007it. The optical spectra show stronger than average [O I] emission lines and weaker than average [Ca II] lines, which may indicate a 16 - 27 M$_{\sun}$ progenitor, on the higher end of expected Type IIP masses. Multi-component [O I] lines are also seen in the optical spectra, most likely caused by an asymmetric blob or a torus of oxygen core material being ejected during the SN explosion. Interaction with circumstellar material prior to day 540 may have created a cool dense shell between the forward and reverse shocks where new dust is condensing. At late times there is also a flattening of the visible lightcurve as the ejecta luminosity fades and a surrounding light echo becomes visible. Radiative transfer models of SN 2007it SEDs indicate that up to 10$^{-4}$ M$_{\sun}$ of new dust has formed in the ejecta, which is consistent with the amount of dust formed in other core collapse supernovae.

The Cosmic Origins Spectrograph (COS) was installed in the Hubble Space Telescope in May, 2009 as part of Servicing Mission 4 to provide high sensitivity, medium and low resolution spectroscopy at far- and near-ultraviolet wavelengths (FUV, NUV). COS is the most sensitive FUV/NUV spectrograph flown to date, spanning the wavelength range from 900{\AA} to 3200{\AA} with peak effective area approaching 3000 cm^2. This paper describes instrument design, the results of the Servicing Mission Orbital Verification (SMOV), and the ongoing performance monitoring program.

SN 2007od exhibits characteristics that have rarely been seen in a Type IIP supernova (SN). Optical V band photometry reveals a very steep brightness decline between the plateau and nebular phases of ~4.5 mag, likely due to SN 2007od containing a low mass of 56Ni. The optical spectra show an evolution from normal Type IIP with broad Halpha emission, to a complex, four component Halpha emission profile exhibiting asymmetries caused by dust extinction after day 232. This is similar to the spectral evolution of the Type IIn SN 1998S, although no early-time narrow (~200 km s-1) Halpha component was present in SN 2007od. In both SNe, the intermediate-width Halpha emission components are thought to arise in the interaction between the ejecta and its circumstellar medium (CSM). SN 2007od also shows a mid-IR excess due to new dust. The evolution of the Halpha profile and the presence of the mid-IR excess provide strong evidence that SN 2007od formed new dust before day 232. Late-time observations reveal a flattening of the visible lightcurve. This flattening is a strong indication of the presence of a light echo, which likely accounts for much of the broad, underlying Halpha component seen at late-times. We believe the multi-peaked Halpha emission is consistent with the interaction of the ejecta with a circumstellar ring or torus (for the inner components at \pm1500 km s-1), and a single blob or cloud of circumstellar material out of the plane of the CSM ring (for the outer component at -5000 km s-1). The most probable location for the formation of new dust is in the cool dense shell created by the interaction between the expanding ejecta and its CSM. Monte Carlo radiative transfer modeling of the dust emission from SN 2007od implies that up to 4x 10-4Msun of new dust has formed. This is similar to the amounts of dust formed in other CCSNe such as SNe 1999em, 2004et, and 2006jc.

SN 2008S erupted in early 2008 in the grand design spiral galaxy NGC 6946. The progenitor was detected by Prieto et al. in Spitzer Space Telescope images taken over the four years prior to the explosion, but was not detected in deep optical images, from which they inferred a self-obscured object with a mass of about 10 Msun. We obtained Spitzer observations of SN 2008S five days after its discovery, as well as coordinated Gemini and Spitzer optical and infrared observations six months after its outburst. We have constructed radiative transfer dust models for the object before and after the outburst, using the same r^-2 density distribution of pre-existing amorphous carbon grains for all epochs and taking light-travel time effects into account for the early post-outburst epoch. We rule out silicate grains as a significant component of the dust around SN 2008S. The inner radius of the dust shell moved outwards from its pre-outburst value of 85 AU to a post-outburst value of 1250 AU, attributable to grain vaporisation by the light flash from SN 2008S. Although this caused the circumstellar extinction to decrease from Av = 15 before the outburst to 0.8 after the outburst, we estimate that less than 2% of the overall circumstellar dust mass was destroyed. The total mass-loss rate from the progenitor star is estimated to have been (0.5-1.0)x10^-4 Msun yr^-1. The derived dust mass-loss rate of 5x10^-7 Msun yr^-1 implies a total dust injection into the ISM of up to 0.01 Msun over the suggested duration of the self-obscured phase. We consider the potential contribution of objects like SN 2008S to the dust enrichment of galaxies.

We present results of a population synthesis study aimed at examining the role of spin-kick alignment in producing a correlation between the spin period of the first-born neutron star and the orbital eccentricity of observed double neutron star binaries in the Galactic disk. We find spin-kick alignment to be compatible with the observed correlation, but not to alleviate the requirements for low kick velocities suggested in previous population synthesis studies. Our results furthermore suggest low- and high-eccentricity systems may form through two distinct formation channels distinguished by the presence or absence of a stable mass transfer phase before the formation of the second neutron star. The presence of highly eccentric systems in the observed sample of double neutron stars may furthermore support the notion that neutron stars accrete matter when moving through the envelope of a giant companion.