We present high spatial resolution LBTI/NOMIC $9-12$ $\mu m$ images of VY CMa and its massive outflow feature, the Southwest (SW) Clump. Combined with high-resolution imaging from HST ($0.4-1$ $\mu m$) and LBT/LMIRCam ($1-5$ $\mu m$), we isolate the spectral energy distribution (SED) of the clump from the star itself. Using radiative-transfer code DUSTY, we model both the scattered light from VY CMa and the thermal emission from the dust in the clump to estimate the optical depth, mass, and temperature of the SW Clump. The SW Clump is optically thick at 8.9 $\mu m$ with a brightness temperature of $\sim$200 K. With a dust chemistry of equal parts silicates and metallic iron, as well as assumptions on grain size distribution, we estimate a dust mass of $5.4\times10^{-5}\,M_\odot$. For a gas--to--dust ratio of 100, this implies a total mass of $5.4\times10^{-3}\,M_\odot$. Compared to the typical mass-loss rate of VY CMa, the SW Clump represents an extreme, localized mass-loss event from $\lesssim300$ years ago.

The obscuring structure surrounding active galactic nuclei (AGN) can be explained as a dust and gas flow cycle that fundamentally connects the AGN with their host galaxies. This structure is believed to be associated with dusty winds driven by radiation pressure. However, the role of magnetic fields, which are invoked in almost all models for accretion onto a supermassive black hole and outflows, is not thoroughly studied. Here we report the first detection of polarized thermal emission by means of magnetically aligned dust grains in the dusty torus of NGC 1068 using ALMA Cycle 4 polarimetric dust continuum observations ($0.07"$, $4.2$ pc; 348.5 GHz, $860$ $\mu$m). The polarized torus has an asymmetric variation across the equatorial axis with a peak polarization of $3.7\pm0.5$\% and position angle of $109\pm2^{\circ}$ (B-vector) at $\sim8$ pc east from the core. We compute synthetic polarimetric observations of magnetically aligned dust grains assuming a toroidal magnetic field and homogeneous grain alignment. We conclude that the measured 860 $\mu$m continuum polarization arises from magnetically aligned dust grains in an optically thin region of the torus. The asymmetric polarization across the equatorial axis of the torus arises from 1) an inhomogeneous optical depth, and 2) a variation of the velocity dispersion, i.e. variation of the magnetic field turbulence at sub-pc scales, from the eastern to the western region of the torus. These observations and modeling constrain the torus properties beyond spectral energy distribution results. This study strongly supports that magnetic fields up to a few pc contribute to the accretion flow onto the active nuclei.

HST/STIS spectra of the small clumps and filaments closest to the central star in VY CMa reveal that the very strong K I emission and TiO and VO molecular emission, long thought to form in a dusty circumstellar shell, actually originate in a few small clumps 100's of AU from the star. The K I lines are 10 to 20 times stronger in these nearest ejecta than on the star. The observations also confirm VO as a circumstellar molecule. In this letter we discuss the spectra of the features, their motions and ages, and the identification of the molecular emission. The strength of the atomic and molecular features in the small clumps present an astrophysical problem for the excitation process. We show that the clumps must have a nearly clear line of sight to the star's radiation.

We present the results of spectroscopy and multi-wavelength photometry of luminous and variable star candidates in the nearby spiral galaxies NGC 2403 and M81. We discuss specific classes of stars, the Luminous Blue Variables (LBVs), B[e] supergiants (sgB[e]), and the high luminosity yellow hypergiants. We identify two new LBV candidates, and three sgB[e] stars in M81. We also find that some stars previously considered LBV candidates are actually field stars. The confirmed and candidate LBVs and sgB[e] stars together with the other confirmed members are shown on the HR Diagrams for their respective galaxies. We also present the HR Diagrams for the two "SN impostors", V37 (SN2002kg) and V12(SN1954J) in NGC 2403 and the stars in their immediate environments.

The Stratospheric Observatory for Infrared Astronomy (SOFIA) is a Boeing 747SP aircraft modified to accommodate a 2.7 meter gyro-stabilized telescope, which is mainly focused to studying the Universe at infrared wavelengths. As part of the Strategic Director's Discretionary Time (S-DDT) program, SOFIA performs observations of relevant science cases and immediately offers science-ready data products to the astronomical community. We present the first data release of the S-DDT program on far-infrared imaging polarimetric observations of 30 Doradus using the High-resolution Airborne Wideband Camera-Plus (HAWC+) at 53, 89, 154, and 214 micron. We present the status and quality of the observations, an overview of the SOFIA data products, and examples of working with HAWC+ polarimetric data that will enhance the scientific analysis of this, and future, data sets. These observations illustrate the potential influence of magnetic fields and turbulence in a star-forming region within the Tarantula Nebula.

New MMT/MIRAC (9-11 {\mu}m), SOFIA/FORCAST (11-37 {\mu}m), and Herschel/PACS (70 and 160 {\mu}m) infrared (IR) imaging and photometry is presented for three famous OH/IR red supergiants (NML Cyg, VX Sgr, and S Per) and two normal red supergiants (RS Per and T Per). We model the observed spectral energy distributions (SEDs) using radiative transfer code DUSTY. Azimuthal average profiles from the SOFIA/FORCAST imaging, in addition to dust mass distribution profiles from DUSTY, constrain the mass-loss histories of these supergiants. For all of our observed supergiants, the DUSTY models suggest that constant mass-loss rates do not produce enough dust to explain the observed infrared emission in the stars' SEDs. Combining our results with Shenoy et al. (2016) (Paper I) we find mixed results with some red supergiants showing evidence for variable and high mass-loss events while others have constant mass loss over the past few thousand years.

We describe new results on two supernova impostors in NGC 2403, SN 1954J(V12) and SN 2002kg(V37). For the famous object SN 1954J we combine four critical observations: its current SED, its Halpha emission line profile, the Ca II triplet in absorption in its red spectrum, and the brightness compared to its pre-event state. Together these strongly suggest that the survivor is now a hot supergiant with T ~ 20000 K, a dense wind, substantial circumstellar extinction, and a G-type supergiant companion. The hot star progenitor of V12's giant eruption was likely in the post-red supergiant stage and had already shed a lot of mass. V37 is a classical LBV/S Dor variable. Our photometry and spectra observed during and after its eruption show that its outburst was an apparent transit on the HR Diagram due to enhanced mass loss and the formation of a cooler, dense wind. V37 is an evolved hot supergiant at ~10^6 Lsun with a probable initial mass of 60 -80 Msun.

In this series of papers we have presented the results of a spectroscopic survey of luminous and variable stars in the nearby spirals M31 and M33. In this paper, we present spectroscopy of 132 additional luminous stars, variables, and emission line objects. Most of the stars have emission line spectra, including LBVs and candidate LBVs, Fe II emission line stars and the B[e] supergiants, and the warm hypergiants. Many of these objects are spectroscopically similar and are often confused with each other. With this large spectroscopic data set including various types of emission line stars, we examine their similarities and differences and propose the following criteria that can be used to help distinguish these stars in future work: 1. The B[e] supergiants have emission lines of [O I] and [Fe II] in their spectra. Most of the spectroscopically confirmed sgB[e] stars also have warm circumstellar dust in their SEDs. 2. Confirmed LBVs do not have the [O I] emission lines in their spectra. Some LBVs have [Fe II] emission lines, but not all. Their SEDS shows free-free emission in the near-infrared but no evidence for warm dust. Their most important and defining characteristic is the S Dor-type variability. 3. The warm hypergiants spectroscopically resemble both the LBVs in their eruption or dense wind state and the B[e] supergiants. However, they are very dusty. Some have [Fe II] and [O I] emission in their spectra like the sgB[e] stars, but can be distinguished by their absorption line spectra characteristic of A and F-type supergiants. In contrast, the B[e] supergiant spectra have strong continua and few if any apparent absorption lines.

We present images and spectral energy distributions (SEDs) of massive young stellar objects (YSOs) in three star-forming H II regions of the Large Magellanic Cloud: N159A, N159 Papillon, and N160. We use photometry from SOFIA/FORCAST at 25.3--37.1 um to constrain model fits to the SEDs and determine luminosities, ages, and dust content of the embedded YSOs and their local environments. By placing these sources on mid-infrared color-magnitude and color-color diagrams, we analyze their dust properties and consider their evolutionary status. Since each object in the FORCAST images has an obvious bright near-infrared counterpart in Spitzer Space Telescope images, we do not find any evidence for new, very cool, previously-undiscovered Class 0 YSOs. Additionally, based on its mid-infrared colors and model parameters, N159A is younger than N160 and the Papillon. The nature of the first extragalactic protostars in N159, P1 and P2, is also discussed.

We present high-angular (0.17$-$0.35 arcsec) resolution imaging polarimetric observations of Mrk 231 in the 3.1 $\mu$m filter using MMT-Pol on the 6.5-m MMT, and in the 8.7 $\mu$m, 10.3 $\mu$m, and 11.6 $\mu$m filters using CanariCam on the 10.4-m Gran Telescopio CANARIAS. In combination with already published observations, we compile the 1$-$12 $\mu$m total and polarized nuclear spectral energy distribution (SED). The total flux SED in the central 400 pc is explained as the combination of 1) a hot (731 $\pm$ 4 K) dusty structure, directly irradiated by the central engine, which is at 1.6 $\pm$ 0.1 pc away and attributed to be in the pc-scale polar region, 2) an optically-thick, smooth and disk-like dusty structure (`torus') with an inclination of 48 $\pm$ 23$^{\circ}$ surrounding the central engine, and 3) an extinguished (A$_{\mbox{V}} =$ 36 $\pm$ 5 mag) starburst component. The polarized SED decreases from 0.77 $\pm$ 0.14 per cent at 1.2 $\mu$m to 0.31 $\pm$ 0.15 per cent at 11.6 $\mu$m and follows a power-law function, $\lambda^{\sim0.57}$. The polarization angle remains constant ($\sim$108$^{\circ}$) in the 1$-$12 $\mu$m wavelength range. The dominant polarization mechanism is explained as scattering off hot dust grains in the pc-scale polar regions.

The supernova impostor PSN J09132750+7627410 in NGC 2748 reached a maximum luminosity of approximately -14 mag. It was quickly realized that its was not a true supernova, but another example of a non-terminal giant eruption. PSN J09132750+7627410 is distinguished by multiple P Cygni absorption minima in the Balmer emission lines that correspond to outflow velocities of -400, -1100, and -1600 km/s. Multiple outflows have been observed in only a few other objects. In this paper we describe the evolution of the spectrum and the P Cygni profiles for three months past maximum, the post-maximum formation of a cool, dense wind, and the identification of a possible progenitor. One of the possible progenitors is an infrared source. Its pre-eruption spectral energy distribution suggests a bolometric luminosity of -8.3 mag and a dust temperature of 780 degrees K. If it is the progenitor it is above the AGB limit unlike the intermediate luminosity red transients. The three P Cygni profiles could be due to ejecta from the current eruption, the wind of the progenitor, or previous mass loss events. We suggest that they were all formed as part of the same high mass loss event and are due to material ejected at different velocities or energies. We also suggest that multiple outflows during giant eruptions may be more common than reported.

In a recent paper, Smith and Tombleson (2015) state that the Luminous Blue Variables (LBVs) in the Milky Way and the Magellanic Clouds are isolated; that they are not spatially associated with young O-type stars. They propose a novel explanation that would overturn the standard view of LBVs. In this paper we test their hypothesis for the LBVs in M31 and M33 as well as the LMC and SMC. In M31 and M33, the LBVs are associated with luminous young stars and supergiants appropriate to their luminosities and positions on the HR Diagram. Moreover, in the Smith and Tombleson scenario most of the LBVs should be runaway stars, but the stars' velocities are consistent with their positions in the respective galaxies. In the Magellanic Clouds, those authors' sample was a mixed population. We reassess their analysis, removing seven stars that have no clear relation to LBVs. When we separate the more massive classical and the less luminous LBVs, the classical LBVs have a distribution similar to the late O-type stars, while the less luminous LBVs have a distribution like the red supergiants. None of the confirmed LBVs have high velocities or are candidate runaway stars. These results support the accepted description of LBVs as evolved massive stars that have shed a lot of mass, and are now close to their Eddington limit.

Recent supernova and transient surveys have revealed an increasing number of non-terminal stellar eruptions. Though the progenitor class of these eruptions includes the most luminous stars, little is known of the pre-supernova mechanics of massive stars in their most evolved state, thus motivating a census of possible progenitors. From surveys of evolved and unstable luminous star populations in nearby galaxies, we select a sample of yellow and red supergiant candidates in M31 and M33 for review of their spectral characteristics and spectral energy distributions. Since the position of intermediate and late-type supergiants on the color-magnitude diagram can be heavily contaminated by foreground dwarfs, we employ spectral classification and multi-band photometry from optical and near-infrared surveys to confirm membership. Based on spectroscopic evidence for mass loss and the presence of circumstellar dust in their SEDs, we find that $30-40\%$ of the yellow supergiants are likely in a post-red supergiant state. Comparison with evolutionary tracks shows that these mass-losing, post-RSGs have initial masses between $20-40\,M_{\odot}$. More than half of the observed red supergiants in M31 and M33 are producing dusty circumstellar ejecta. We also identify two new warm hypergiants in M31, J004621.05+421308.06 and J004051.59+403303.00, both of which are likely in a post-RSG state.

We report the fifth confirmed Luminous Blue Variable/S Doradus variable in M31. In 2006, J004526.62+415006.3 had the spectrum of hot Fe II emission line star with strong P Cygni profiles in the Balmer lines. In 2010, its absorption line spectrum resembled an early A-type supergiant with H and Fe II emission lines with strong P Cygni profiles, and in 2013 the spectrum had fully transitioned to an F-type supergiant due to the formation of the optically thick, cool wind which characterizes LBVs at maximum light. The photometric record supports the LBV/S Dor nature of the variability. Its bolometric luminosity ~ -9.65 mag places it on the HR Diagram near the known LBVs, AE And, Var C in M33 and S Dor.

We discuss the spectrum of Var C in M33 obtained just before the onset of its current brightening and recent spectra during its present "eruption" or optically thick wind stage. These spectra illustrate the typical LBV transition in apparent spectral type or temperature that characterizes the classical LBV or S Dor-type variability. LBVs are known to have slow, dense winds during their maximum phase. Interestingly, Var C had a slow wind even during its hot, quiescent stage in comparison with the normal hot supergiants with similar temperatures. Its outflow or wind speeds also show very little change between these two states.

We present constraints on cosmological and astrophysical parameters from high-resolution microwave background maps at 148 GHz and 218 GHz made by the Atacama Cosmology Telescope (ACT) in three seasons of observations from 2008 to 2010. A model of primary cosmological and secondary foreground parameters is fit to the map power spectra and lensing deflection power spectrum, including contributions from both the thermal Sunyaev-Zeldovich (tSZ) effect and the kinematic Sunyaev-Zeldovich (kSZ) effect, Poisson and correlated anisotropy from unresolved infrared sources, radio sources, and the correlation between the tSZ effect and infrared sources. The power ell^2 C_ell/2pi of the thermal SZ power spectrum at 148 GHz is measured to be 3.4 +\- 1.4 muK^2 at ell=3000, while the corresponding amplitude of the kinematic SZ power spectrum has a 95% confidence level upper limit of 8.6 muK^2. Combining ACT power spectra with the WMAP 7-year temperature and polarization power spectra, we find excellent consistency with the LCDM model. We constrain the number of effective relativistic degrees of freedom in the early universe to be Neff=2.79 +\- 0.56, in agreement with the canonical value of Neff=3.046 for three massless neutrinos. We constrain the sum of the neutrino masses to be Sigma m_nu < 0.39 eV at 95% confidence when combining ACT and WMAP 7-year data with BAO and Hubble constant measurements. We constrain the amount of primordial helium to be Yp = 0.225 +\- 0.034, and measure no variation in the fine structure constant alpha since recombination, with alpha/alpha0 = 1.004 +/- 0.005. We also find no evidence for any running of the scalar spectral index, dns/dlnk = -0.004 +\- 0.012.

The Keck Array (SPUD) is a set of microwave polarimeters that observes from the South Pole at degree angular scales in search of a signature of Inflation imprinted as B-mode polarization in the Cosmic Microwave Background (CMB). The first three Keck Array receivers were deployed during the 2010-2011 Austral summer, followed by two new receivers in the 2011-2012 summer season, completing the full five-receiver array. All five receivers are currently observing at 150 GHz. The Keck Array employs the field-proven BICEP/BICEP2 strategy of using small, cold, on-axis refractive optics, providing excellent control of systematics while maintaining a large field of view. This design allows for full characterization of far-field optical performance using microwave sources on the ground. We describe our efforts to characterize the main beam shape and beam shape mismatch between co-located orthogonally-polarized detector pairs, and discuss the implications of measured differential beam parameters on temperature to polarization leakage in CMB analysis.

The BICEP2 and Keck Array experiments are designed to measure the polarization of the cosmic microwave background (CMB) on angular scales of 2-4 degrees (l=50-100). This is the region in which the B-mode signal, a signature prediction of cosmic inflation, is expected to peak. BICEP2 was deployed to the South Pole at the end of 2009 and is in the middle of its third year of observing with 500 polarization-sensitive detectors at 150 GHz. The Keck Array was deployed to the South Pole at the end of 2010, initially with three receivers--each similar to BICEP2. An additional two receivers have been added during the 2011-12 summer. We give an overview of the two experiments, report on substantial gains in the sensitivity of the two experiments after post-deployment optimization, and show preliminary maps of CMB polarization from BICEP2.

We present photometry for all bright red giant branch (RGB), horizontal branch (HB), and asymptotic giant branch (AGB) stars within 10' of the center of M13. We find support for the idea that the population of HB stars redder than the primary group are noticeably evolved, which resolves a disagreement between distance moduli derived from the tip of the RGB and from stars near the instability strip. The sharp cut at the red end of the HB provides strong evidence that stars from the dominant HB group must still be undergoing blue loops, implying that diffusion is being inhibited. We argue that M13's HB is a somewhat pathological case - the dominant HB population occurs very near the "knee" in optical CMDs, and evolved stars exclusively appear redward of that peak, leading to the incorrect appearance of a continuation of the unevolved HB. M13 has a distinct group of HB stars previously identified with the second U jump, which may be examples of early hot flashers that ignite core helium fusion shortly after leaving the RGB. However, there is not convincing evidence that a large fraction of stars leave the RGB before helium flash. We revisited the helium-sensitive R ratio, and find that M13's ratio is in agreement with theoretical values for primordial helium abundance Y_P = 0.245 and inconsistent with a helium enhancement DY = 0.04. The brightness of the HB (both in comparison to the end of the canonical HB and to the tip of the RGB) also appears to rule out the idea that the envelopes of the reddest HB stars have been significantly enriched in helium. The absolute colors of the turnoffs of M3 and M13 may potentially be used to look for differences in their mean helium abundances.(ABRIDGED)

We conduct a search for radio stars by combining radio and optical data from the Faint Images of the Radio Sky at Twenty cm survey (FIRST) and the Sloan Digital Sky Survey (SDSS). The faint limit of SDSS makes possible a homogeneous search for radio emission from stars of low optical luminosity. We select a sample of 112 candidate radio stars in the magnitude range $15<i\lesssim19.1$ and with radio flux $S_{20}\geq1.25$mJy, from about 7000 deg$^2$ of sky. The selection criteria are positional coincidence within $1\arcsec$, radio and optical point source morphology, and an SDSS spectrum classified as stellar. The sample contamination is estimated by random matching to be $108\pm13$, suggesting that at most a small fraction of the selected candidates are genuine radio stars. Therefore, we rule out a very rare population of extremely radio-loud stars: no more than 1.2 of every million stars in the magnitude range $15<i<19.1$ stars has radio flux $S_{20}\geq1.25$ mJy. We investigate the optical and radio colors of the sample to find candidates that show the largest likelihood of being real radio stars. The significant outliers from the stellar locus, as well as the magnetically active stars, are the best candidates for follow-up radio observations. We conclude that, while the present wide-area radio surveys are not sensitive enough to provide homogeneous samples of the extremely rare radio stars, upcoming surveys which exploit the great sensitivity of current and planned telescopes do have sufficient sensitivity and will allow the properties of this class of object to be investigated in detail.

We describe briefly the properties of the recently completed Southern African Large Telescope (SALT), along with its first light imager SALTICAM. Using this instrument, we present 4.3 hr of high speed unfiltered photometric observations of the eclipsing polar SDSSJ015543.40+002807.2 with time resolution as short as 112 ms, the highest quality observations of this kind of any polar to date. The system was observed during its high luminosity state. Two accreting poles are clearly seen in the eclipse light curve. The binary system parameters have been constrained: the white dwarf mass is at the low end of the range expected for cataclysmic variables. Correlations between the positions of the accretion regions on or near the surface of the white dwarf and the binary system parameters were established. The sizes of the accretion regions and their relative movement from eclipse to eclipse were estimated: they are typically 4-7 deg depending on the mass of the white dwarf. The potential of these observations will only fully be realised when low state data of the same kind are obtained and the contact phases of the eclipse of the white dwarf are measured.