We present SOFIA/FORCAST imaging of the circumstellar dust shells surrounding the luminous blue variable (LBV) candidates MN 90 and HD 168625 to quantify the mineral abundances of the dust and to constrain the evolutionary state of these objects. Our image at 37.1 micron of MN 90 shows a limb-brightened, spherical dust shell. A least-squares fit to the spectral energy distribution of MN 90 yields a dust temperature of 59 \pm 10 K, with the peak of the emission at 42.7 micron. Using 2-Dust radiative transfer code, we estimate for MN 90 that mass-loss occurred at a rate of (7.3 \pm 0.4)x10^-7 M_sun/yr x (v_exp/50 km/s) to create a dust shell with a dust mass of (3.2 \pm 0.1)x10^-2 M_sun. Our images between 7.7 - 37.1 micron of HD 168625 complement previously obtained mid-IR imaging of its bipolar nebulae. The SOFIA/FORCAST imaging of HD 168625 shows evidence for the limb-brightened peaks of an equatorial torus. We estimate a dust temperature of 170 \pm 40 K for the equatorial dust surrounding HD 168625, with the peak of the emission at 18.3 micron. Our 2-Dust model for HD 168625 estimates that mass-loss occurred at a rate of (3.2 \pm 0.2)x10^-7 M_sun/yr to create a dust torus/shell with a dust mass of (2.5 \pm 0.1)x10^-3 M_sun.

We present 5-28 micron SOFIA FORECAST spectroscopy complemented by panchromatic X-ray through infrared observations of the CO nova V5668 Sgr documenting the formation and destruction of dust during 500 days following outburst. Dust condensation commenced by 82 days after outburst at a temperature of 1090 K. The condensation temperature indicates that the condensate was amorphous carbon. There was a gradual decrease of the grain size and dust mass during the recovery phase. Absolute parameter values given here are for an assumed distance of 1.2 kpc. We conclude that the maximum mass of dust produced was 1.2 x 10-7 solar masses if the dust was amorphous carbon. The average grain radius grew to a maximum of 2.9 microns at a temperature of 720 K around day 113 when the shell visual optical depth was Tau = 5.4. Maximum grain growth was followed by followed by a period of grain destruction. X-rays were detected with Swift from day 95 to beyond day 500. The Swift X-ray count rate due to the hot white dwarf peaked around day 220, when its spectrum was that of a kT = 35 eV blackbody. The temperature, together with the super-soft X-ray turn-on and turn-off times, suggests a WD mass of 1.1 solar masses. We show that the X-ray fluence was sufficient to destroy the dust. Our data show that the post-dust event X-ray brightening is not due to dust destruction, which certainly occurred, as the dust is optically thin to X-rays.

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 present a SOFIA FORCAST grism spectroscopic survey to examine the mineralogy of the circumstellar dust in a sample of post-asymptotic giant branch yellow supergiants that are believed to be the precursors of planetary nebulae. Our mineralogical model of each star indicates the presence of both carbon rich and oxygen rich dust species-contrary to simple dredge-up models-with a majority of the dust in the form of amorphous carbon and graphite. The oxygen rich dust is primarily in the form of amorphous silicates. The spectra do not exhibit any prominent crystalline silicate emission features. For most of the systems, our analysis suggests that the grains are relatively large and have undergone significant processing, supporting the hypothesis that the dust is confined to a Keplerian disk and that we are viewing the heavily processed, central regions of the disk from a nearly face-on orientation. These results help to determine the physical properties of the post-AGB circumstellar environment and to constrain models of post-AGB mass loss and planetary nebula formation.

We present an ongoing, systematic search for extragalactic infrared transients, dubbed SPIRITS --- SPitzer InfraRed Intensive Transients Survey. In the first year, using Spitzer/IRAC, we searched 190 nearby galaxies with cadence baselines of one month and six months. We discovered over 1958 variables and 43 transients. Here, we describe the survey design and highlight 14 unusual infrared transients with no optical counterparts to deep limits, which we refer to as SPRITEs (eSPecially Red Intermediate Luminosity Transient Events). SPRITEs are in the infrared luminosity gap between novae and supernovae, with [4.5] absolute magnitudes between -11 and -14 (Vega-mag) and [3.6]-[4.5] colors between 0.3 mag and 1.6 mag. The photometric evolution of SPRITEs is diverse, ranging from < 0.1 mag/yr to > 7 mag/yr. SPRITEs occur in star-forming galaxies. We present an in-depth study of one of them, SPIRITS 14ajc in Messier 83, which shows shock-excited molecular hydrogen emission. This shock may have been triggered by the dynamic decay of a non-hierarchical system of massive stars that led to either the formation of a binary or a proto-stellar merger.

We present infrared spectroscopy of the classical nova V339 Delphini, obtained over a $\sim2$ year period. The infrared emission lines were initially symmetrical, with HWHM velocities of 525 km s$^{-1}$. In later ($t\gtrsim77$days, where $t$ is the time from outburst) spectra however, the lines displayed a distinct asymmetry, with a much stronger blue wing, possibly due to obscuration of the receding component by dust. Dust formation commenced at $\sim$ day 34.75 at a condensation temperature of $1480\pm20$K, consistent with graphitic carbon. Thereafter the dust temperature declined with time as $T_{\rm d}\propto{t}^{-0.346}$, also consistent with graphitic carbon. The mass of dust initally rose, as a result of an increase in grain size and/or number, peaked at $\sim$ day 100, and then declined precipitously. This decline was most likely caused by grain shattering due to electrostatic stress after the dust was exposed to X-radiation. An Appendix summarises Planck Means for carbon, and the determination of grain mass and radius for a carbon dust shell.

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.

We present mid- and far- IR imaging of four famous hypergiant stars: the red supergiants $\mu$ Cep and VY CMa, and the warm hypergiants IRC +10420 and $\rho$ Cas. Our 11 to 37 $\mu$m SOFIA/FORCAST imaging probes cool dust not detected in visual and near-IR imaging studies. Adaptive optics (AO) 8 - 12 $\mu$m imaging of $\mu$ Cep and IRC +10420 with MMT/MIRAC reveals extended envelopes that are the likely sources of these stars' strong silicate emission features. We find $\mu$ Cep's mass-loss rate to have declined by about a factor of 5 over a 13,000 history, ranging from 5 $\times$ 10$^{-6}$ down to $\sim$1 $\times$ 10$^{-6}$ $M_{\odot}$ yr$^{-1}$. The morphology of VY CMa indicates a cooler dust component coincident with the highly asymmetric reflection nebulae seen in the visual and near-IR. The lack of cold dust at greater distances around VY CMa indicates its mass-loss history is limited to the last $\sim$1200 years, with an average rate of 6 $\times$ 10$^{-4}$ $M_{\odot}$ yr$^{-1}$. We find two distinct periods in the mass-loss history of IRC +10420 with a high rate of 2 $\times$ 10$^{-3}$ $M_{\odot}$ yr$^{-1}$ until approximately 2000 yr ago, followed by an order of magnitude decrease in the recent past. We interpret this change as evidence of its evolution beyond the RSG stage. Our new infrared photometry of $\rho$ Cas is consistent with emission from the expanding dust shell ejected in its 1946 eruption, with no evidence of newer dust formation from its more recent events.

We describe recent panchromatic imaging and spectroscopic studies of the supergiant, mass-transferring, binary star RY Scuti, which is in a brief transitional phase to become a Wolf-Rayet star and a stripped-envelope supernova.

We present J' and K' imaging linear polarimetric adaptive optics observations of NGC 1068 using MMT-Pol on the 6.5-m MMT. These observations allow us to study the torus from a magnetohydrodynamical (MHD) framework. In a 0.5" (30 pc) aperture at K', we find that polarisation arising from the passage of radiation from the inner edge of the torus through magnetically aligned dust grains in the clumps is the dominant polarisation mechanism, with an intrinsic polarisation of 7.0%$\pm$2.2%. This result yields a torus magnetic field strength in the range of 4$-$82 mG through paramagnetic alignment, and 139$^{+11}_{-20}$ mG through the Chandrasekhar-Fermi method. The measured position angle (P.A.) of polarisation at K$'$ is found to be similar to the P.A. of the obscuring dusty component at few parsec scales using infrared interferometric techniques. We show that the constant component of the magnetic field is responsible for the alignment of the dust grains, and aligned with the torus axis onto the plane of the sky. Adopting this magnetic field configuration and the physical conditions of the clumps in the MHD outflow wind model, we estimate a mass outflow rate $\le$0.17 M$_{\odot}$ yr$^{-1}$ at 0.4 pc from the central engine for those clumps showing near-infrared dichroism. The models used were able to create the torus in a timescale of $\geq$10$^{5}$ yr with a rotational velocity of $\leq$1228 km s$^{-1}$ at 0.4 pc. We conclude that the evolution, morphology and kinematics of the torus in NGC 1068 can be explained within a MHD framework.

We present 2 - 5 micron adaptive optics (AO) imaging and polarimetry of the famous hypergiant stars IRC +10420 and VY Canis Majoris. The imaging polarimetry of IRC +10420 with MMT-Pol at 2.2 micron resolves nebular emission with intrinsic polarization of 30%, with a high surface brightness indicating optically thick scattering. The relatively uniform distribution of this polarized emission both radially and azimuthally around the star confirms previous studies that place the scattering dust largely in the plane of the sky. Using constraints on scattered light consistent with the polarimetry at 2.2 micron, extrapolation to wavelengths in the 3 - 5 micron band predicts a scattered light component significantly below the nebular flux that is observed in our LBT/LMIRCam 3 - 5 micron AO imaging. Under the assumption this excess emission is thermal, we find a color temperature of ~ 500 K is required, well in excess of the emissivity-modified equilibrium temperature for typical astrophysical dust. The nebular features of VY CMa are found to be highly polarized (up to 60%) at 1.3 micron, again with optically thick scattering required to reproduce the observed surface brightness. This star's peculiar nebular feature dubbed the "Southwest Clump" is clearly detected in the 3.1 micron polarimetry as well, which, unlike IRC+10420, is consistent with scattered light alone. The high intrinsic polarizations of both hypergiants' nebulae are compatible with optically thick scattering for typical dust around evolved dusty stars, where the depolarizing effect of multiple scatters is mitigated by the grains' low albedos.

We present adaptive optics images of the extreme red supergiant VY Canis Majoris in the Ks, L' and M bands (2.15 to 4.8 micron) made with LMIRCam on the Large Binocular Telescope (LBT). The peculiar "Southwest Clump" previously imaged from 1 to 2.2 micron appears prominently in all three filters. We find its brightness is due almost entirely to scattering, with the contribution of thermal emission limited to at most 25%. We model its brightness as optically thick scattering from silicate dust grains using typical size distributions. We find a lower limit mass for this single feature of 5E-03 Msun to 2.5E-02 Msun depending on the assumed gas-to-dust ratio. The presence of the Clump as a distinct feature with no apparent counterpart on the other side of the star is suggestive of an ejection event from a localized region of the star and is consistent with VY CMa's history of asymmetric high mass loss events.

We present Spitzer Space Telescope mid-infrared IRS spectra, supplemented by ground-based optical observations, of the classical novae V1974 Cyg, V382 Vel, and V1494 Aql more than 11, 8, and 4 years after outburst respectively. The spectra are dominated by forbidden emission from neon and oxygen, though in some cases, there are weak signatures of magnesium, sulfur, and argon. We investigate the geometry and distribution of the late time ejecta by examination of the emission line profiles. Using nebular analysis in the low density regime, we estimate lower limits on the abundances in these novae. In V1974 Cyg and V382 Vel, our observations confirm the abundance estimates presented by other authors and support the claims that these eruptions occurred on ONe white dwarfs. We report the first detection of neon emission in V1494 Aql and show that the system most likely contains a CO white dwarf.