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.

We present the Gravitational Wave Treasure Map, a tool to coordinate, visualize, and assess the electromagnetic follow-up of gravitational wave (GW) events. With typical GW localization regions of hundreds to thousands of square degrees and dozens of active follow-up groups, the pursuit of electromagnetic (EM) counterparts is a challenging endeavor, but the scientific payoff for early discovery of any counterpart is clear. With this tool, we provide a website and API interface that allows users to easily see where other groups have searched and better inform their own follow-up search efforts. A strong community of Treasure Map users will increase the overall efficiency of EM counterpart searches and will play a fundamental role in the future of multi-messenger astronomy.

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

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

Hydrogen-rich, core-collapse supernovae are typically divided into four classes: IIP, IIL, IIn, and IIb. In general, interaction with circumstellar material is only considered for Type IIn supernovae. However, recent hydrodynamic modeling of IIP and IIL supernovae requires circumstellar material to reproduce their early light curves. In this scenario, IIL supernovae experience large amounts of mass loss before exploding. We test this hypothesis on ASASSN-15oz, a Type IIL supernova. With extensive follow-up in the X- ray, UV, optical, IR, and radio we present our search for signs of interaction, and the mass-loss history indicated by their detection. We find evidence of short-lived intense mass-loss just prior to explosion from light curve modeling, amounting in 1.5 M$_{\odot}$ of material within 1800 R$_{\odot}$ of the progenitor. We also detect the supernova in the radio, indicating mass-loss rates of $10^{-6}-10^{-7}$ M$_{\odot}$ yr$^{-1}$ prior to the extreme mass-loss period. Our failure to detect the supernova in the X-ray and the lack of narrow emission lines in the UV, optical, and NIR do not contradict this picture and place an upper limit on the mass-loss rate outside the extreme period of $<10^{-4}$ M$_{\odot}$ yr$^{-1}$. This paper highlights the importance gathering comprehensive data on more Type II supernovae to enable detailed modeling of the progenitor and supernova which can elucidate their mass-loss histories and envelope structures and thus inform stellar evolution models.

The Hubble Space Telescope photometric survey of Galactic globular clusters (GCs) has revealed a peculiar "chromosome map" for NGC6934. Besides a typical sequence, similar to that observed in Type I GCs, NGC 6934 displays additional stars on the red side, analogous to the anomalous, Type II GCs, as defined in our previous work. We present a chemical abundance analysis of four red giants in this GC. Two stars are located on the chromosome map sequence common to all GCs, and another two on the additional sequence. We find: (i) star-to-star Fe variations, with the two anomalous stars being enriched by ~0.2 dex. Due to our small-size sample, this difference is at the ~2.5 sigma level; (ii) no evidence for variations in the slow neutron-capture abundances over Fe, at odds with what is often observed in anomalous Type II GCs, e.g. M 22 and Omega Centauri; (iii) no large variations in light elements C, O and Na, compatible with the targets location on the lower part of the chromosome map where such variations are not expected. Since the analyzed stars are homogeneous in light elements, the only way to reproduce the photometric splits on the sub-giant (SGB) and the red-giant (RGB) branches is to assume that red-RGB/faint-SGB stars are enhanced in [Fe/H] by ~0.2. This fact corroborates the spectroscopic evidence of a metallicity variation in NGC6934. The observed chemical pattern resembles only partially the other Type II GCs, suggesting that NGC6934 might belong either to a third class of GCs, or be a link between normal Type I and anomalous Type II GCs.

The Cygnus OB2 Association is one of the nearest and largest collections of massive stars in the Galaxy. Situated at the heart of the "Cygnus X" complex of star-forming regions and molecular clouds, its distance has proven elusive owing to the ambiguous nature of kinematic distances along this $\ell\simeq80$ degree sightline and the heavy, patchy extinction. In an effort to refine the three-dimensional geometry of key Cygnus~X constituents, we have measured distances to four eclipsing double-lined OB-type spectroscopic binaries that are probable members of Cyg~OB2. We find distances of $1.33\pm0.17$, $1.32\pm0.07$, $1.44\pm0.18$, and $1.32\pm0.13$ kpc toward MT91~372, MT91~696, CPR2002~A36, and Schulte~3 respectively. We adopt a weighted average distance of 1.33$\pm$0.06~kpc. This agrees well with spectrophotometric estimates for the Association as a whole and with parallax measurements of protostellar masers in the surrounding interstellar clouds, thereby linking the ongoing star formation in these clouds with Cyg~OB2. We also identify Schulte 3C (O9.5V), a 4" visual companion to the 4.75 day binary Schulte~3(A+B), as a previously unrecognized Association member.

We analyze orbital solutions for 48 massive multiple-star systems in the Cygnus OB2 Association, 23 of which are newly presented here, to find that the observed distribution of orbital periods is approximately uniform in log P for P<45 d, but it is not scale-free. Inflections in the cumulative distribution near 6 d, 14, d, and 45 d, suggest key physical scales of ~0.2, ~0.4, and ~1 A.U. where yet-to-be-identified phenomena create distinct features. No single power law provides a statistically compelling prescription, but if features are ignored, a power law with exponent beta = -0.22 provides a crude approximation over P=1.4 -- 2000 d, as does a piece-wise linear function with a break near 45 d. The cumulative period distribution flattens at P > 45 d, even after correction for completeness, indicating either a lower binary fraction or a shift toward low-mass companions. A high degree of similarity (91% likelihood) between the Cyg OB2 period distribution and that of other surveys suggests that the binary properties at P<25 d are determined by local physics of disk/clump fragmentation and are relatively insensitive to environmental and evolutionary factors. Fully 30% of the unbiased parent sample is a binary with period P < 45 d. Completeness corrections imply a binary fraction near 55% for P < 5000 d. The observed distribution of mass ratios 0.2<q<1 is consistent with uniform, while the observed distribution of eccentricities 0.1<e<0.6 is consistent with uniform plus an excess of e ~ 0 systems. We identify six stars, all supergiants, that exhibit aperiodic velocity variations of ~30 km/s attributed to atmospheric fluctuations.

We present an all-sky sample of 984 candidate intermediate-mass Galactic star-forming regions color-selected from the Infrared Astronomical Satellite (IRAS) Point Source Catalog and morphologically classify each object using mid-infrared Wide-field Infrared Survey Explorer (WISE) images. Of the 984 candidates, 616 are probable star-forming regions (62.6%), 128 are filamentary structures (13.0%), 39 are point-like objects of unknown nature (4.0%), and 201 are galaxies (20.4%). We conduct a study of four of these regions, IRAS 00259+5625, IRAS 00420+5530, IRAS 01080+5717, and IRAS 05380+2020, at Galactic latitudes |b| > 5 degrees using optical spectroscopy from the Wyoming Infrared Observatory along with near-infrared photometry from the Two-Micron All Sky Survey to investigate their stellar content. New optical spectra, color-magnitude diagrams, and color-color diagrams reveal their extinctions, spectrophotometric distances, and the presence of small stellar clusters containing 20-78 solar masses of stars. These low-mass diffuse star clusters contain 65-250 stars for a typical initial mass function, including one or more mid-B stars as their most massive constituents. Using infrared spectral energy distributions we identify young stellar objects near each region and assign probable masses and evolutionary stages to the protostars. The total infrared luminosity lies in the range 190 to 960 solar luminosities, consistent with the sum of the luminosities of the individually identified young stellar objects.

Massive binary stars may constitute a substantial fraction of progenitors to supernovae and gamma-ray bursts, and the distribution of their orbital characteristics holds clues to the formation process of massive stars. As a contribution to securing statistics on OB-type binaries, we report the discovery and orbital parameters for five new systems as part of the Cygnus OB2 Radial Velocity Survey. Four of the new systems (MT070, MT174, MT267, and MT734 (a.k.a. VI Cygni #11) are single-lined spectroscopic binaries while one (MT103) is a double-lined system (B1V+B2V). MT070 is noteworthy as the longest period system yet measured in Cyg OB2, with P=6.2 yr. The other four systems have periods ranging between 4 and 73 days. MT174 is noteworthy for having a probable mass ratio q<0.1, making it a candidate progenitor to a low-mass X-ray binary. These measurements bring the total number of massive binaries in Cyg OB2 to 25, the most currently known in any single cluster or association.

We report the discovery and orbital solutions for two new OB binaries in the Cygnus OB2 Association, MT311 (B2V+B3V) and MT605 (B0.5V+B2.5:V). We also identify the system MT429 as a probable triple system consisting of a tight eclipsing 2.97 day B3V+B6V pair and a B0V at a projected separation of 138 AU. We further provide the first spectroscopic orbital solutions to the eclipsing, double-lined, O-star binary MT696 (O9.5V+B1:V), the double-lined, early B binary MT720 (B0-1V+B1-2V), and the double-lined, O-star binary MT771 (O7V+O9V). These systems exhibit orbital periods between 1.5 days and 12.3 days, with the majority having P<6 days. The two new binary discoveries and six spectroscopic solutions bring the total number of known massive binaries in the central region of the Cygnus OB2 Association to 20, with all but two having full orbital solutions.

Mid-infrared images from the Spitzer Space Telescope Galactic Legacy Infrared MidPlane Survey Extraordinaire program reveal that the infrared source IRAS 03063+5735 is a bowshock nebula produced by an early B star, 2MASS 03101044+5747035. We present new optical spectra of this star, classify it as a B1.5 V, and determine a probable association with a molecular cloud complex at V_LSR=-38 -- -42 km/s in the outer Galaxy near l=140.59 degr, b=-0.250 degr. On the basis of spectroscopic parallax, we estimate a distance of 4.0 +/-1 kpc to both the bowshock nebula and the molecular complex. One plausible scenario is that this a high-velocity runaway star impinging upon a molecular cloud. We identify the HII region and stellar cluster associated with IRAS 03064+5638 at a projected distance of 64 pc as one plausible birth site. The spectrophotometric distance and linkage to a molecular feature provides another piece of data helping to secure the ill-determined rotation curve in the outer Galaxy. As a by-product of spectral typing this star, we present empirical spectral diagnostic diagrams suitable for approximate spectral classification of O and B stars using He lines in the little-used yellow-red portion of the optical spectrum.