We report updated results for the first 7 years of our program to monitor 27 galaxies within 10 Mpc using the Large Binocular Telescope to search for failed supernovae (SNe) -- core-collapses of massive stars that form black holes without luminous supernovae. In the new data, we identify no new compelling candidates and confirm the existing candidate. Given the 6 successful core-collapse SNe in the sample and one likely failed SN, the implied fraction of core-collapses that result in failed SNe is $f = 0.14^{+0.33}_{-0.10}$ at 90% confidence. If the current candidate is a failed SN, the fraction of failed SN naturally explains the missing high-mass red supergiant SN progenitors and the black hole mass function. If the current candidate is ultimately rejected, the data implies a 90% confidence upper limit on the failed SN fraction of $f < 0.35$.

We present Hubble Space Telescope imaging confirming the optical disappearance of the failed supernova (SN) candidate identified by Gerke et al. (2015). This $\sim 25~M_{\odot}$ red supergiant experienced a weak $\sim 10^{6}~L_{\odot}$ optical outburst in 2009 and is now at least 5 magnitudes fainter than the progenitor in the optical. The mid-IR flux has slowly decreased to the lowest levels since the first measurements in 2004. There is faint ($2000-3000~L_{\odot}$) near-IR emission likely associated with the source. We find the late-time evolution of the source to be inconsistent with obscuration from an ejected, dusty shell. Models of the spectral energy distribution indicate that the remaining bolometric luminosity is $>6$ times fainter than that of the progenitor and is decreasing as $\sim t^{-4/3}$. We conclude that the transient is unlikely to be a SN impostor or stellar merger. The event is consistent with the ejection of the envelope of a red supergiant in a failed SN and the late-time emission could be powered by fallback accretion onto a newly-formed black hole. Future IR and X-ray observations are needed to confirm this interpretation of the fate for the star.

We identify and phase a sample of 81 Cepheids in the maser-host galaxy NGC 4258 using the Large Binocular Telescope (LBT), and obtain calibrated mean magnitudes in up to 4 filters for a subset of 43 Cepheids using archival HST data. We employ 3 models to study the systematic effects of extinction, the assumed extinction law, and metallicity on the Cepheid distance to NGC 4258. We find a correction to the Cepheid colors consistent with a grayer extinction law in NGC 4258 compared to the Milky Way ($R_V =4.9$), although we believe this is indicative of other systematic effects. If we combine our Cepheid sample with previously known Cepheids, we find a significant metallicity adjustment to the distance modulus of $\gamma_1 = -0.61 \pm 0.21$ mag/dex, for the Zaritsky et al. (1994) metallicity scale, as well as a weak trend of Cepheid colors with metallicity. Conclusions about the absolute effect of metallicity on Cepheid mean magnitudes appear to be limited by the available data on the metallicity gradient in NGC 4258, but our Cepheid data require at least some metallicity adjustment to make the Cepheid distance consistent with independent distances to the LMC and NGC 4258. From our ensemble of models and the geometric maser distance of NGC 4258 ($\mu_{N4258} = 29.40 \pm 0.06$ mag), we estimate $\mu_{LMC} = 18.57 \pm 0.14$ mag ($51.82 \pm 3.23$ kpc).

We are monitoring 27 galaxies within 10 Mpc using the Large Binocular Telescope to search for failed supernovae (SNe), massive stars that collapse to form a black hole without a SN explosion. We present the results from the first 4 years of survey data, during which these galaxies were observed to produce 3 successful core-collapse SNe. We search for stars that have "vanished" over the course of our survey, by examining all stars showing a decrease in luminosity of $\Delta \nu L_{\nu} \ge 10^4L_{\odot}$ from the first to the last observation. We also search for the low luminosity, long duration transients predicted by \citet{lovegrove2013} for failed explosions of red supergiants. After analyzing the first 4 years of data in this first direct search for failed SNe, we are left with one candidate requiring further study. This candidate has an estimated mass of 18-25$M_{\odot}$, a mass range likely associated with failed SNe and, if real, implies that failed SN represent a median fraction of $f \simeq 0.30$ of core-collapses, with symmetric 90% confidence limits of $0.07 \le f \le 0.62$. If follow up data eliminate this candidate, we find an upper limit on the fraction of core collapses leading to a failed SN of $f<0.40$ at 90% confidence. As the duration of the survey continues to increase, it will begin to constrain the $f \simeq$ 10-30% failure rates needed to explain the deficit of massive SN progenitors and the observed black hole mass function.

We discuss moderate resolution spectra, multicolor photometry, and light curves of thirty-one of the most luminous stars and variables in the giant spiral M101. The majority are intermediate A to F-type supergiants. We present new photometry and light curves for three known "irregular blue variables" V2, V4 and V9) and identify a new candidate. Their spectra and variability confirm that they are LBV candidates and V9 may be in an LBV-like maximum light state or eruption.

Understanding the late-stage evolution of the most massive stars such as $\eta$ Carinae is challenging because no true analogs of $\eta$ Car have been clearly identified in the Milky Way or other galaxies. In Khan et. al. (2013), we utilized Spitzer IRAC images of $7$ nearby ($\lesssim4$ Mpc) galaxies to search for such analogs, and found $34$ candidates with flat or red mid-IR spectral energy distributions. Here, in Paper II, we present our characterization of these candidates using multi-wavelength data from the optical through the far-IR. Our search detected no true analogs of $\eta$ Car, which implies an eruption rate that is a fraction $0.01\lesssim F \lesssim 0.19$ of the ccSN rate. This is roughly consistent with each $M_{ZAMS} \gtrsim 70M_\odot$ star undergoing $1$ or $2$ outbursts in its lifetime. However, we do identify a significant population of $18$ lower luminosity $\left(\log(L/L_\odot)\simeq5.5-6.0\right)$ dusty stars. Stars enter this phase at a rate that is fraction $0.09 \lesssim F \lesssim 0.55$ of the ccSN rate, and this is consistent with all $25 < M_{ZAMS} < 60M_\odot$ stars undergoing an obscured phase at most lasting a few thousand years once or twice. These phases constitute a negligible fraction of post-main sequence lifetimes of massive stars, which implies that these events are likely to be associated with special periods in the evolution of the stars. The mass of the obscuring material is of order $\sim M_\odot$, and we simply do not find enough heavily obscured stars for theses phases to represent more than a modest fraction ($\sim 10\%$ not $\sim 50\%$) of the total mass lost by these stars. In the long term, the sources that we identified will be prime candidates for detailed physical analysis with JWST.

We measure the integrated contributions of dusty AGB stars and other luminous red mid-IR sources to the mid-IR luminosities of 6 galaxies (M81, NGC 2403, NGC 300, M33 and the Magellanic Clouds). We find the dusty AGB stars whose mid-IR fluxes are dominated by dust rather than photospheric emission contribute from 0.6% (M81) to 5.6% (SMC) of the 3.6 micron flux and 1.0% (M81) to 10.1% (SMC) of the 4.5 micron flux. We find a trend of decreasing AGB contribution with increasing galaxy metallicity, luminosity and mass and decreasing SSFR. However, these galaxy properties are strongly correlated in our sample and the simplest explanation of the trend is galaxy metallicity. Bright, red sources other than dusty AGB stars represent a smaller fraction of the luminosity, ~1.2% at 3.6 microns, however their dust is likely cooler and their contributions are likely larger at longer wavelengths. Excluding the SMC, the contribution from these red sources correlates with the specific star formation rate as we would expect for massive stars. In total, after correcting for dust emission at other wavelengths, the dust around AGB stars radiates 0.1-0.8% of the bolometric luminosities of the galaxies. Thus, hot dust emission from AGB and other luminous dusty stars represent a small fraction of the total luminosities of the galaxies but a significant fraction of their mid-IR emissions.

We identify and phase a sample of 107 Cepheids with 10<P/days<100 in M81 using the LBT and calibrate their BVI mean magnitudes with archival HST data. The use of a ground-based telescope to identify and phase the Cepheids and HST only for the final calibration reduces the demand on HST by nearly an order of magnitude and yields Period-Luminosity (PL) relations with dispersions comparable to the best LMC samples. We fit the sample using the OGLE-II LMC PL relations and are unable to find a self-consistent distance for different band combinations or radial locations within M81. We can do so after adding a radial dependence to the PL zero point that corresponds to a luminosity dependence on metallicity of g_mu=-0.56+/-0.36 mag/dex. We find marginal evidence for a shift in color as a function of metallicity, distinguishable from the effects of extinction, of g_2=+0.07+/-0.03 mag/dex. We find a distance modulus for M81, relative to the LMC, of mu(M81-LMC)=9.39+/-0.14 mag, including uncertainties due to the metallicity corrections. This corresponds to a distance to M81 of 3.6+/-0.2 Mpc, assuming a LMC distance modulus of 18.41 mag. We carry out a joint analysis of M81 and NGC4258 Cepheids and simultaneously solve for the distance of M81 relative to NGC4258 and the metallicity corrections. Given the current data, the uncertainties of such joint fits are dominated by the relative metallicities and the abundance gradients rather than by measurement errors of the Cepheid magnitudes or colors. We find mu(M81-LMC)=9.40 (-0.11/+0.15) mag, mu(N4258-LMC)=11.08 (-0.17/+0.21) mag and mu(N4258-M81)=1.68+/-0.08 mag and joint metallicity corrections of g_mu=-0.62 (-0.35/+0.31) mag/dex and g_2=0.01+/-0.01 mag/dex. Quantitative analyses of Cepheid distances must take into account both the metallicity dependencies of the Cepheids and the uncertainties in the abundance estimates. (ABRIDGED)

We are continuing a J, K and narrow-band imaging survey of 300 square degrees of the plane of the Galaxy, searching for new Wolf-Rayet stars. Our survey spans 150 degrees in Galactic longitude and reaches 1 degree above and below the Galactic plane. The survey has a useful limiting magnitude of K = 15 over most of the observed Galactic plane, and K = 14 (due to severe crowding) within a few degrees of the Galactic center. Thousands of emission line candidates have been detected. In spectrographic follow-ups of 146 relatively bright WR star candidates we have re-examined 11 previously known WC and WN stars and discovered 71 new WR stars, 17 of type WN and 54 of type WC. Our latest image analysis pipeline now picks out WR stars with a 57% success rate. Star subtype assignments have been confirmed with K band spectra, and distances approximated using the method of spectroscopic parallax. Some of the new WR stars are amongst the most distant known in our Galaxy. The distribution of these new WR stars is beginning to trace the locations of massive stars along the distant spiral arms of the Milky Way.

We show that the candidate progenitor of the core-collapse SN 2011dh in M51 (8 Mpc away) was fading by 0.039 +- 0.006 mag/year during the three years prior to the supernova, and that this level of variability is moderately unusual for other similar stars in M 51. While there are uncertainties about whether the true progenitor was a blue companion to this candidate, the result illustrates that there are no technical challenges to obtaining fairly high precision light curves of supernova progenitors using ground based observations of nearby (<10 Mpc) galaxies with wide field cameras on 8m-class telescopes. While other sources of variability may dominate, it is even possible to reach into the range of evolution rates required by the quasi-static evolution of the stellar envelope. For M 81, where we have many more epochs and a slightly longer time baseline, our formal 3 sigma sensitivity to slow changes is presently 3 millimag/year for a M_V ~= -8 mag star. In short, there is no observational barrier to determining whether the variability properties of stars in their last phases of evolution (post Carbon ignition) are different from earlier phases.

The discovery of new Wolf-Rayet (WR) stars in our Galaxy via large-scale narrowband optical surveys has been severely limited by dust extinction. Recent improvements in infrared technology have made narrowband-broadband imaging surveys viable again. We report a new J, K and narrow-band imaging survey of 300 square degrees of the plane of the Galaxy, spanning 150 degrees in Galactic longitude and reaching 1 degree above and below the Galactic plane. The survey has a useful limiting magnitude of K = 15 over most of the observed Galactic plane, and K = 14 within a few degrees of the Galactic center. Thousands of emission line candidates have been detected. In spectrographic follow-ups of 173 WR star candidates we have discovered 41 new WR stars, 15 of type WN and 26 of type WC. Star subtype assignments have been confirmed with K band spectra, and distances approximated using the method of spectroscopic parallax. A few of the new WR stars are amongst the most distant known in our Galaxy. The distribution of these new WR stars is seen to follow that of previously known WR stars along the spiral arms of the Galaxy. Tentative radial velocities were also measured for most of the new WR stars.

Long term optical and near-infrared photometric and blue spectroscopic observations were obtained for QS Tel, V834 Cen, and BL Hyi. The optical light curves of all three polars displayed large magnitude changes during our observations. These same high/low state transitions were also apparent in near-infrared JHK photometry, though with decreased amplitude. The color of the polar with respect to its state was examined and found not to be a good indicator of the instantaneous state. During low to high state transitions, a nearly constant magnitude difference was observed in all three polars. This $\Delta$m value was found to be consistent with the level expected to occur if accretion onto the white dwarf reached the Eddington luminosity during the high state. The high state Balmer decrement was measured for each star and used to estimate that the temperature of the emission line forming region was ~12,000K with N$_H$ near 12.8 dex. No relationship between the Balmer emission line strength and the white dwarf magnetic field strength was seen, in contrast to a good correlation between these two parameters observed for UV emission lines.