We present a dust spectral energy distribution (SED) and binary stellar population analysis revisiting the dust production rates (DPRs) in the winds of carbon-rich Wolf-Rayet (WC) binaries and their impact on galactic dust budgets. DustEM SED models of 19 Galactic WC ``dustars" reveal DPRs of $\dot{M}_d\sim10^{-10}-10^{-6}$ M$_\odot$ yr$^{-1}$ and carbon dust condensation fractions, $\chi_C$, between $0.002 - 40\%$. A large ($0.1 - 1.0$ $\mu$m) dust grain size composition is favored for efficient dustars where $\chi_C\gtrsim1\%$. Results for dustars with known orbital periods verify a power-law relation between $\chi_C$, orbital period, WC mass-loss rate, and wind velocity consistent with predictions from theoretical models of dust formation in colliding-wind binaries. We incorporated dust production into Binary Population and Spectral Synthesis (BPASS) models to analyze dust production rates from WC dustars, asymptotic giant branch stars (AGBs), red supergiants (RSGs), and core-collapse supernovae (SNe). BPASS models assuming constant star formation (SF) and a co-eval $10^6$ M$_\odot$ stellar population were performed at low, Large Magellanic Cloud (LMC)-like, and solar metallicities (Z = 0.001, 0.008, and 0.020). Both constant SF and co-eval models show that SNe are net dust destroyers at all metallicities. Constant SF models at LMC-like metallicities show that AGB stars slightly outproduce WC binaries and RSGs by factors of $2-3$, whereas at solar metallicites WC binaries are the dominant source of dust for $\sim60$ Myr until the onset of AGBs, which match the dust input of WC binaries. Co-eval population models show that for "bursty" SF, AGB stars dominate dust production at late times ($t\gtrsim 70$ Myr).

The Galactic Center contains some of the most extreme conditions for star formation in our Galaxy as well as many other phenomena that are unique to this region. Given our relative proximity to the Galactic Center, we are able to study details of physical processes to a level that is simply not yet possible for more distant galaxies, yielding an otherwise inaccessible view of the nuclear region of a galaxy. We recently carried out a targeted imaging survey of mid-infrared bright portions of the Galactic Center at 25 and 37 $\mu$m using the FORCAST instrument on SOFIA. This survey was one of the inaugural Legacy Programs from SOFIA cycle 7, observing a total area of 403 arcmin$^2$ (2180 pc$^2$), including the Sgr A, B, and C complexes. Here we present an overview of the survey strategy, observations, and data reduction as an accompaniment to the initial public release of the survey data. We discuss interesting regions and features within the data including extended features near the circumnuclear disk, structures in the Arched Filaments and Sickle H II regions, and signs of embedded star formation in Sgr B2 and Sgr C. We also feature a handful of less well studied mid-infrared sources located between Sgr A and Sgr C that could be sites of relatively isolated star formation activity. Last, we discuss plans for subsequent publications and future data releases from the survey.

In 2010 May, an intermediate luminosity optical transient was discovered in the nearby galaxy NGC 300 by a South African amateur astronomer. In the decade since its discovery, multi-wavelength observations of the misnamed ``SN 2010da'' have continually re-shaped our understanding of this high mass X-ray binary system. In this review, we present an overview of the multi-wavelength observations and attempts to understand the 2010 transient event and, later, the re-classification of this system as NGC~300 ULX-1: a red supergiant + neutron star ultraluminous X-ray source.

We present a systematic study of the most luminous ($M_{\mathrm{IR}}$ [Vega magnitudes] brighter than $-14$) infrared (IR) transients discovered by the SPitzer InfraRed Intensive Transients Survey (SPIRITS) between 2014 and 2018 in nearby galaxies ($D < 35$ Mpc). The sample consists of nine events that span peak IR luminosities of $M_{[4.5],\mathrm{peak}}$ between $-14$ and $-18.2$, show IR colors between $0.2 < ([3.6]{-}[4.5]) < 3.0$, and fade on timescales between $55$ days $< t_{\mathrm{fade}} < 480$ days. The two reddest events ($A_V > 12$) show multiple, luminous IR outbursts over several years and have directly detected, massive progenitors in archival imaging. With analyses of extensive, multiwavelength follow-up, we suggest the following possible classifications: five obscured core-collapse supernovae (CCSNe), two erupting massive stars, one luminous red nova, and one intermediate-luminosity red transient. We define a control sample of all optically discovered transients recovered in SPIRITS galaxies and satisfying the same selection criteria. The control sample consists of eight CCSNe and one Type Iax SN. We find that 7 of the 13 CCSNe in the SPIRITS sample have lower bounds on their extinction of $2 < A_V < 8$. We estimate a nominal fraction of CCSNe in nearby galaxies that are missed by optical surveys as high as $38.5^{+26.0}_{-21.9}$% (90% confidence). This study suggests that a significant fraction of CCSNe may be heavily obscured by dust and therefore undercounted in the census of nearby CCSNe from optical searches.

(Abridged) Palomar Gattini-IR is a new wide-field, near-infrared robotic time domain survey operating at Palomar Observatory. Using a 30 cm telescope mounted with a H2RG detector, Gattini-IR achieves a field of view of 25 sq. deg. with a pixel scale of 8.7" in J-band. Here, we describe the system design, survey operations, data processing system and on-sky performance of Palomar Gattini-IR. As a part of the nominal survey, Gattini-IR scans $\approx 7500$ square degrees of the sky every night to a median 5$\sigma$ depth of $15.7$ AB mag outside the Galactic plane. The survey covers $\approx 15000$ square degrees of the sky visible from Palomar with a median cadence of 2 days. A real-time data processing system produces stacked science images from dithered raw images taken on sky, together with PSF-fit source catalogs and transient candidates identified from subtractions within a median delay of $\approx 4$ hours from the time of observation. The calibrated data products achieve an astrometric accuracy (RMS) of $\approx 0.7$" with respect to Gaia DR2 for sources with S/N $> 10$, and better than $\approx 0.35$" for sources brighter than $\approx 12$ Vega mag. The photometric accuracy (RMS) achieved in the PSF-fit source catalogs is better than $\approx 3$% for sources brighter than $\approx 12$ Vega mag, as calibrated against the 2MASS catalog. With a field of view $\approx 40\times$ larger than any other existing near infrared imaging instrument, Gattini-IR is probing the reddest and dustiest transients in the local universe such as dust obscured supernovae in nearby galaxies, novae behind large columns of extinction within the galaxy, reddened micro-lensing events in the Galactic plane and variability from cool and dust obscured stars. We present results from transients and variables identified since the start of the commissioning period.

Supernova (SN) 2014C is a unique explosion where a seemingly typical hydrogen-poor stripped envelope SN started to interact with a dense, hydrogen-rich circumstellar medium (CSM) a few months after the explosion. The delayed interaction suggests a detached CSM shell, unlike in a typical SN IIn where the CSM is much closer and the interaction commences earlier post-explosion; indicating a different mass loss history. We present near- to mid-infrared observations of SN 2014C from 1-5 years after the explosion, including uncommon 9.7 $\mu$m imaging with COMICS on the Subaru telescope. Spectroscopy shows that the interaction is still ongoing, with the intermediate-width He I 1.083 $\mu$m emission present out to our latest epoch 1639 days post-explosion. The last Spitzer/IRAC photometry at 1920 days post-explosion further confirms ongoing CSM interaction. The 1-10 $\mu$m spectral energy distributions (SEDs) can be explained by a dust model with a mixture of 69% carbonaceous and 31% silicate dust, pointing to a chemically inhomogeneous CSM. The inference of silicate dust is the first among interacting SNe. An SED model with purely carbonaceous CSM dust is possible, but would require more than 0.22 $M_{\odot}$ of dust, which is an order of magnitude larger than what observed in any other SNe, measured in the same way, at this epoch. The light curve beyond 500 days is well fit by an interaction model with a wind-driven CSM and a mass loss rate of $\sim 10^{-3} \, M_{\odot}\,\rm yr^{-1}$, which presents an additional CSM density component exterior to the constant density shell reported previously in the literature. SN 2014C could originate in a binary system, similar to RY Scuti, which would explain the observed chemical and density profile inhomogeneity in the CSM.

SN2010da/NGC 300 ULX-1 was first detected as a supernova impostor in May 2010 and was recently discovered to be a pulsating ultraluminous X-ray source. In this letter, we present VLT/X-shooter spectra of this source obtained in October 2018, covering the wavelength range 350-2300 nm. The $J$- and $H$-bands clearly show the presence of a red supergiant donor star that is best matched by a MARCS stellar atmosphere with $T_{\rm eff} = 3650 - 3900$ K and $\log(L_{\rm bol}/L_{\odot}) = 4.25\pm0.10$, which yields a stellar radius $R = 310 \pm 70 R_{\odot}$. To fit the full spectrum, two additional components are required: a blue excess that can be fitted either by a hot blackbody (T $\gtrsim 20,000$ K) or a power law (spectral index $\alpha \approx 4$) and is likely due to X-ray emission reprocessed in the outer accretion disk or the donor star; and a red excess that is well fitted by a blackbody with a temperature of $\sim 1100$ K, and is likely due to warm dust in the vicinity of SN2010da. The presence of a red supergiant in this system implies an orbital period of at least 0.8-2.1 years, assuming Roche lobe overflow. Given the large donor-to-compact object mass ratio, orbital modulations of the radial velocity of the red supergiant are likely undetectable. However, the radial velocity amplitude of the neutron star is large enough (up to 40-60 km s$^{-1}$) to potentially be measured in the future, unless the system is viewed at a very unfavorable inclination.

We present the discovery of an optical transient (OT) in Messier 51, designated M51 OT2019-1 (also ZTF19aadyppr, AT 2019abn, ATLAS19bzl), by the Zwicky Transient Facility (ZTF). The OT rose over 15 days to an observed luminosity of $M_r=-13$ (${\nu}L_{\nu}=9\times10^6~L_{\odot}$), in the luminosity gap between novae and typical supernovae (SNe). Spectra during the outburst show a red continuum, Balmer emission with a velocity width of $\approx400$ km s$^{-1}$, Ca II and [Ca II] emission, and absorption features characteristic of an F-type supergiant. The spectra and multiband light curves are similar to the so-called "SN impostors" and intermediate-luminosity red transients (ILRTs). We directly identify the likely progenitor in archival Spitzer Space Telescope imaging with a $4.5~\mu$m luminosity of $M_{[4.5]}\approx-12.2$ and a $[3.6]-[4.5]$ color redder than 0.74 mag, similar to those of the prototype ILRTs SN 2008S and NGC 300 OT2008-1. Intensive monitoring of M51 with Spitzer further reveals evidence for variability of the progenitor candidate at [4.5] in the years before the OT. The progenitor is not detected in pre-outburst Hubble Space Telescope optical and near-IR images. The optical colors during outburst combined with spectroscopic temperature constraints imply a higher reddening of $E(B-V)\approx0.7$ mag and higher intrinsic luminosity of $M_r\approx-14.9$ (${\nu}L_{\nu}=5.3\times10^7~L_{\odot}$) near peak than seen in previous ILRT candidates. Moreover, the extinction estimate is higher on the rise than on the plateau, suggestive of an extended phase of circumstellar dust destruction. These results, enabled by the early discovery of M51 OT2019-1 and extensive pre-outburst archival coverage, offer new clues about the debated origins of ILRTs and may challenge the hypothesis that they arise from the electron-capture induced collapse of extreme asymptotic giant branch stars.

We present a catalog of 417 luminous infrared variable stars with periods exceeding 250 days. These were identified in 20 nearby galaxies by the ongoing SPIRITS survey with the Spitzer Space Telescope. Of these, 359 variables have $M_{[4.5]}$ (phase-weighted mean magnitudes) fainter than $-12$ and periods and luminosities consistent with previously reported variables in the Large Magellanic Cloud. However, 58 variables are more luminous than $M_{[4.5]} = -12$, including 11 that are brighter than $M_{[4.5]} = -13$ with the brightest having $M_{[4.5]} = -15.51$. Most of these bright variable sources have quasi-periods longer than 1000 days, including four over 2000 days. We suggest that the fundamental period-luminosity relationship, previously measured for the Large Magellanic Cloud, extends to much higher luminosities and longer periods in this large galaxy sample. We posit that these variables include massive AGB stars (possibly super-AGB stars), red supergiants experiencing exceptionally high mass-loss rates, and interacting binaries. We also present 3.6, 4.5, 5.8 and 8.0 $\mu$m photometric catalogs for all sources in these 20 galaxies.

We present a mid-infrared (IR) sample study of nearby ultraluminous X-ray sources (ULXs) using multi-epoch observations with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. Spitzer/IRAC observations taken after 2014 were obtained as part of the Spitzer Infrared Intensive Transients Survey (SPIRITS). Our sample includes 96 ULXs located within 10 Mpc. Of the 96~ULXs, 12 have candidate counterparts consistent with absolute mid-IR magnitudes of supergiants, and 16 counterparts exceeded the mid-IR brightness of single supergiants and are thus more consistent with star clusters or non-ULX background active galactic nuclei (AGN). The supergiant candidate counterparts exhibit a bi-modal color distribution in a Spitzer/IRAC color-magnitude diagram, where "red" and "blue" ULXs fall in IRAC colors $[3.6] - [4.5]\sim0.7$ and $[3.6] - [4.5]\sim0.0$, respectively. The mid-IR colors and absolute magnitudes of 4 "red" and 5 "blue" ULXs are consistent with that of supergiant B[e] (sgB[e]) and red supergiant (RSG) stars, respectively. While "blue", RSG-like mid-IR ULX counterparts likely host RSG mass donors, we propose the "red" counterparts are ULXs exhibiting the "B[e] phenomenon'' rather than hosts of sgB[e] mass donors. We show that the mid-IR excess from the "red" ULXs is likely due to thermal emission from circumstellar or circumbinary dust. Using dust as a probe for total mass, we estimate mass-loss rates of $\dot{M}\sim1\times10^{-4}$ M$_\odot$ yr$^{-1}$ in dust-forming outflows of red ULXs. Based on the transient mid-IR behavior and its relatively flat spectral index, $\alpha=-0.19\pm0.1$, we suggest that the mid-IR emission from Holmberg IX X-1 originates from a variable jet.

Opening up the dynamic infrared sky for systematic time-domain exploration would yield many scientific advances. Multi-messenger pursuits such as localizing gravitational waves from neutron star mergers and quantifying the nucleosynthetic yields require the infrared. Another multi-messenger endeavor that needs infrared surveyors is the study of the much-awaited supernova in our own Milky Way. Understanding shocks in novae, true rates of supernovae and stellar mergers are some other examples of stellar evolution and high energy physics wherein the answers are buried in the infrared. We discuss some of the challenges in the infrared and pathfinders to overcome them. We conclude with recommendations on both infrared discovery engines and infrared follow-up machines that would enable this field to flourish in the next decade.

We present infrared (IR) photometry and spectroscopy of the Type II-P SN 2017eaw and its progenitor in the nearby galaxy NGC 6946. Progenitor observations in the Ks band in 4 epochs from 1 year to 1 day before the explosion reveal no significant variability in the progenitor star greater than 6% that last longer than 200 days. SN 2017eaw is a typical SN II-P with near-IR and mid-IR photometric evolution similar to those of SNe 2002hh and 2004et, other normal SNe II-P in the same galaxy. Spectroscopic monitoring between 389 and 480 days post explosion reveals strong CO first overtone emission at 389 d, with a line profile matching that of SN 1987A from the same epoch, indicating $\sim 10^{-3} \, M_{\odot}$ of CO at 1,800 K. From the 389 d epoch until the most recent observation at 566 d, the first overtone feature fades while the 4.5 $\mu$m excess, likely from the CO fundamental band, remains. This behavior indicates that the CO has not been destroyed, but that the gas has cooled enough that the levels responsible for first overtone emissions are no longer populated. Finally, the evolution of Spitzer 3.6 $\mu$m photometry shows evidence for dust formation in SN 2017eaw, with a dust mass of $10^{-6}$ or $10^{-4}\,M_{\odot}$ assuming carbonaceous or silicate grains respectively.

We report our Spitzer Space Telescope observations and detections of the binary neutron star merger GW170817. At 4.5um, GW170817 is detected at 21.9 mag AB at +43 days and 23.9 mag AB at +74 days after merger. At 3.6um, GW170817 is not detected to a limit of 23.2 mag AB at +43 days and 23.1 mag AB at +74 days. Our detections constitute the latest and reddest constraints on the kilonova/macronova emission and composition of heavy elements. The 4.5um luminosity at this late phase cannot be explained by elements exclusively from the first abundance peak of the r-process. Moreover, the steep decline in the Spitzer band, with a power-law index of 3.4 +/- 0.2, can be explained by a few of the heaviest isotopes in the third abundance peak with half-life around 14 days dominating the luminosity (e.g. 140Ba, 143Pr, 147Nd, 156Eu, 191Os, 223Ra, 225Ra, 233Pa, 234Th) or a model with lower deposition efficiency. This data offers evidence that the heaviest elements in the second and third r-process abundance peak were indeed synthesized. Our conclusion is verified by both analytics and network simulations and robust despite intricacies and uncertainties in the nuclear physics. Future observations with Spitzer and James Webb Space Telescope will further illuminate the relative abundance of the synthesized heavy elements.

We present the first spatially resolved mid-infrared (37.1 $\mu$m) image of the Fomalhaut debris disk. We use PSF fitting and subtraction to distinctly measure the flux from the unresolved component and the debris disk. We measure an infrared excess in the point source of $0.9 \pm 0.2$ Jy, consistent with emission from warm dust in an inner disk structure (Su et al. 2016), and inconsistent with a stellar wind origin. We cannot confirm or rule out the presence of a pileup ring (Su et al. 2016) near the star. In the cold region, the 37 $\mu$m imaging is sensitive to emission from small, blowout grains, which is an excellent probe of the dust production rate from planetesimal collisions. Under the assumptions that the dust grains are icy aggregates and the debris disk is in steady state, this result is consistent with the dust production rates predicted by Kenyon & Bromley (2008) from theoretical models of icy planet formation. We find a dust luminosity of $(7.9 \pm 0.8) \times 10^{-4}$ L$_\odot$ and a dust mass of 8 -- 16 lunar masses, depending on grain porosity, with $\sim 1$ lunar mass in grains with radius 1 $\mu$m -- 1 mm. If the grains are icy and highly porous, meter-sized objects must be invoked to explain the far-IR, submm, and mm emission. If the grains are composed of astronomical silicates, there is a dearth of blowout grains (Pawellek et al. 2014) and the mass loss rate is well below the predicted dust production values.

We are building a next-generation laser adaptive optics system, Robo-AO-2, for the UH 2.2-m telescope that will deliver robotic, diffraction-limited observations at visible and near-infrared wavelengths in unprecedented numbers. The superior Maunakea observing site, expanded spectral range and rapid response to high-priority events represent a significant advance over the prototype. Robo-AO-2 will include a new reconfigurable natural guide star sensor for exquisite wavefront correction on bright targets and the demonstration of potentially transformative hybrid AO techniques that promise to extend the faintness limit on current and future exoplanet adaptive optics systems.

X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nuclear reaction sequences down to a depth in the inner crust where the mass density is 2E12 g/cm^3 using a full nuclear reaction network for a range of initial compositions. The reaction sequences differ substantially from previous work. We find a robust reduction of crust impurity at the transition to the inner crust regardless of initial composition, though shell effects can delay the formation of a pure crust somewhat to densities beyond 2E12 g/cm^3. This naturally explains the small inner crust impurity inferred from observations of a broad range of systems. The exception are initial compositions with A >= 102 nuclei, where the inner crust remains impure with an impurity parameter of Qimp~20 due to the N = 82 shell closure. In agreement with previous work we find that nuclear heating is relatively robust and independent of initial composition, while cooling via nuclear Urca cycles in the outer crust depends strongly on initial composition. This work forms a basis for future studies of the sensitivity of crust models to nuclear physics and provides profiles of composition for realistic crust models.

Based on phase-resolved broadband spectroscopy using $XMM$-$Newton$ and $NuSTAR$, we report on a potential cyclotron resonant scattering feature at $E \sim 13$ keV in the pulsed spectrum of the recently discoverd ULX pulsar NGC 300 ULX1. If this interpretation is correct, the implied magnetic field of the central neutron star is $B \sim 10^{12}$ G (assuming scattering off electrons), similar to that estimated from the observed spin-up of the star, and also similar to known Galactic X-ray pulsars. We discuss the implications of this result for the connection between NGC 300 ULX1 and the other known ULX pulsars, particularly in light of the recent discovery of a likely proton Cyclotron line in another ULX, M51 ULX-8.

We present the discovery by the SPitzer InfraRed Intensive Transients Survey (SPIRITS) of a likely supernova (SN) in NGC 3556 at only 8.8 Mpc, which, despite its proximity, was not detected by optical searches. A luminous infrared (IR) transient at $M_{[4.5]} = -16.7$ mag (Vega), SPIRITS 16tn is coincident with a dust lane in the inclined, star-forming disk of the host. Using IR, optical, and radio observations, we attempt to determine the nature of this event. We estimate $A_V \approx$ 8 - 9 mag of extinction, placing it among the three most highly obscured IR-discovered SNe to date. The [4.5] light curve declined at a rate of 0.013 mag day$^{-1}$, and the $[3.6] - [4.5]$ color grew redder from 0.7 to $\gtrsim$ 1.0 mag by 184.7 days post discovery. Optical/IR spectroscopy shows a red continuum, but no clearly discernible features, preventing a definitive spectroscopic classification. Deep radio observations constrain the radio luminosity of SPIRITS 16tn to $L_{\nu} \lesssim 10^{24}$ erg s$^{-1}$ Hz$^{-1}$ between 3 - 15 GHz, excluding many varieties of radio core-collapse SNe. A type Ia SN is ruled out by the observed red IR color, and lack of features normally attributed to Fe-peak elements in the optical and IR spectra. SPIRITS 16tn was fainter at [4.5] than typical stripped-envelope SNe by $\approx$ 1 mag. Comparison of the spectral energy distribution to SNe II suggests SPIRITS 16tn was both highly obscured, and intrinsically dim, possibly akin to the low-luminosity SN 2005cs. We infer the presence of an IR dust echo powered by a peak luminosity of the transient of $5 \times 10^{40}$ erg s$^{-1} < L_{\mathrm{peak}} < 4\times10^{43}$ erg s$^{-1}$, consistent with the observed range for SNe II. This discovery illustrates the power of IR surveys to overcome the compounding effects of visible extinction and optically sub-luminous events in completing the inventory of nearby SNe.

We present an optical and infrared (IR) study of IC 10 X-2, a high-mass X-ray binary in the galaxy IC 10. Previous optical and X-ray studies suggest X-2 is a Supergiant Fast X-ray Transient: a large-amplitude (factor of $\sim$ 100), short-duration (hours to weeks) X-ray outburst on 2010 May 21. We analyze R- and g-band light curves of X-2 from the intermediate Palomar Transient Factory taken between 2013 July 15 and 2017 Feb 14 show high-amplitude ($\gtrsim$ 1 mag), short-duration ($\lesssim8$ d) flares and dips ($\gtrsim$ 0.5 mag). Near-IR spectroscopy of X-2 from Palomar/TripleSpec show He I, Paschen-$\gamma$, and Paschen-$\beta$ emission lines with similar shapes and amplitudes as those of luminous blue variables (LBVs) and LBV candidates (LBVc). Mid-IR colors and magnitudes from Spitzer/IRAC photometry of X-2 resemble those of known LBV/LBVcs. We suggest that the stellar companion in X-2 is an LBV/LBVc and discuss possible origins of the optical flares. Dips in the optical light curve are indicative of eclipses from optically thick clumps formed in the winds of the stellar counterpart. Given the constraints on the flare duration ($0.02 - 0.8$ d) and the time between flares ($15.1\pm7.8$ d), we estimate the clump volume filling factor in the stellar winds, $f_V$, to be $0.01 < f_V < 0.71$, which overlaps with values measured from massive star winds. In X-2, we interpret the origin of the optical flares as the accretion of clumps formed in the winds of an LBV/LBVc onto the compact object.

The asymptotic giant branch (AGB) star R Sculptoris (R Scl) is one of the most extensively studied stars on the AGB. R Scl is a carbon star with a massive circumstellar shell ($M_{shell}\sim 7.3\times10^{-3}~M_{\odot}$) which is thought to have been produced during a thermal pulse event $\sim2200$ years ago. To study the thermal dust emission associated with its circumstellar material, observations were taken with the Faint Object InfraRed CAMera for the SOFIA Telescope (FORCAST) at 19.7, 25.2, 31.5, 34.8, and 37.1 $\mu$m. Maps of the infrared emission at these wavelengths were used to study the morphology and temperature structure of the spatially extended dust emission. Using the radiative transfer code DUSTY and fitting the spatial profile of the emission, we find that a geometrically thin dust shell cannot reproduce the observed spatially resolved emission. Instead, a second dust component in addition to the shell is needed to reproduce the observed emission. This component, which lies interior to the dust shell, traces the circumstellar envelope of R Scl. It is best fit by a density profile with $n \propto r^{\alpha}$ where $\alpha=0.75^{+0.45}_{-0.25}$ and dust mass of $M_d=9.0^{+2.3}_{-4.1}\times10^{-6}~M_{\odot}$. The strong departure from an $r^{-2}$ law indicates that the mass-loss rate of R Scl has not been constant. This result is consistent with a slow decline in the post-pulse mass-loss which has been inferred from observations of the molecular gas.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

We present constraints on Ks-band emission from one of the nearest short hard gamma-ray bursts, GRB 160821B, at z=0.16, at three epochs. We detect a reddened relativistic afterglow from the jetted emission in the first epoch but do not detect any excess kilonova emission in the second two epochs. We compare upper limits obtained with Keck I/MOSFIRE to multi-dimensional radiative transfer models of kilonovae, that employ composition-dependent nuclear heating and LTE opacities of heavy elements. We discuss eight models that combine toroidal dynamical ejecta and two types of wind and one model with dynamical ejecta only. We also discuss simple, empirical scaling laws of predicted emission as a function of ejecta mass and ejecta velocity. Our limits for GRB 160821B constrain the ejecta mass to be lower than 0.03 Msun for velocities greater than 0.1c. At the distance sensitivity range of advanced LIGO, similar ground-based observations would be sufficiently sensitive to the full range of predicted model emission including models with only dynamical ejecta. The color evolution of these models shows that I-K color spans 7--16 mag, which suggests that even relatively shallow infrared searches for kilonovae could be as constraining as optical searches.

We present mid-infrared (IR) light curves of the Ultraluminous X-ray Source (ULX) Holmberg II X-1 from observations taken between 2014 January 13 and 2017 January 5 with the \textit{Spitzer Space Telescope} at 3.6 and 4.5 $\mu$m in the \textit{Spitzer} Infrared Intensive Transients Survey (SPIRITS). The mid-IR light curves, which reveal the first detection of mid-IR variability from a ULX, is determined to arise primarily from dust emission rather than from a jet or an accretion disk outflow. We derived the evolution of the dust temperature ($T_\mathrm{d}\sim600 - 800$ K), IR luminosity ($L_\mathrm{IR}\sim3\times10^4$ $\mathrm{L}_\odot$), mass ($M_\mathrm{d}\sim1-3\times10^{-6}$ $\mathrm{M}_\odot$), and equilibrium temperature radius ($R_\mathrm{eq}\sim10-20$ AU). A comparison of X-1 with a sample spectroscopically identified massive stars in the Large Magellanic Cloud on a mid-IR color-magnitude diagram suggests that the mass donor in X-1 is a supergiant (sg) B[e]-star. The sgB[e]-interpretation is consistent with the derived dust properties and the presence of the [Fe II] ($\lambda=1.644$ $\mu$m) emission line revealed from previous near-IR studies of X-1. We attribute the mid-IR variability of X-1 to increased heating of dust located in a circumbinary torus. It is unclear what physical processes are responsible for the increased dust heating; however, it does not appear to be associated with the X-ray flux from the ULX given the constant X-ray luminosities provided by serendipitous, near-contemporaneous X-ray observations around the first mid-IR variability event in 2014. Our results highlight the importance of mid-IR observations of luminous X-ray sources traditionally studied at X-ray and radio wavelengths.

Massive stellar clusters provide radiation ($\mathrm{\sim 10^7-10^8~L_{\odot}}$) and winds ($\mathrm{\sim 1000~km/s}$) that act to heat dust and shape their surrounding environment. In this paper, the Arched Filaments in the Galactic center were studied to better understand the influence of the Arches cluster on its nearby interstellar medium (ISM). The Arched Filaments were observed with the Faint Object InfraRed CAMera for the SOFIA Telescope (FORCAST) at 19.7, 25.2, 31.5, and 37.1 $\mu$m. Color-temperature maps of the region created with the 25.2 and 37.1 $\mu$m data reveal relatively uniform dust temperatures (70-100 K) over the extent of the filaments ($\sim 25$ pc). Distances between the cluster and the filaments were calculated assuming equilibrium heating of standard size ISM dust grains ($\sim$0.1 $\mu$m). The distances inferred by this method are in conflict with the projected distance between the filaments and the cluster, although this inconsistency can be explained if the characteristic grain size in the filaments is smaller ($\sim$0.01 $\mu$m) than typical values. DustEM models of selected locations within the filaments show evidence of depleted abundances of polycyclic aromatic hydrocarbons (PAHs) by factors of $\sim$1.6-10 by mass compared to the diffuse ISM. The evidence for both PAH depletion and a smaller characteristic grain size points to processing of the ISM within the filaments. We argue that the eroding of dust grains within the filaments is not likely attributable to the radiation or winds from the Arches cluster, but may be related to the physical conditions in the Galactic center.

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 high spatial resolution mid-infrared images of the nebula around the late-type carbon-rich Wolf-Rayet (WC)-OB binary system WR~112 taken by the recently upgraded VLT spectrometer and imager for the mid-infrared (VISIR) with the PAH1, NeII\_2, and Q3 filters. The observations reveal a morphology resembling a series of arc-like filaments and broken shells. Dust temperatures and masses are derived for each of the identified filamentary structures, which exhibit temperatures ranging from $179_{-6}^{+8}$ K at the exterior W2 filament to $355_{-25}^{+37}$ K in the central 3". The total dust mass summed over the features is $2.6\pm0.4\times10^{-5}$ $\mathrm{M}_\odot$. A multi-epoch analysis of mid-IR photometry of WR~112 over the past $\sim20$ yr reveals no significant variability in the observed dust temperature and mass. The morphology of the mid-IR dust emission from WR~112 also exhibits no significant expansion from imaging data taken in 2001, 2007, and 2016, which disputes the current interpretation of the nebula as a high expansion velocity ($\sim1200$ km s$^{-1}$) "pinwheel"-shaped outflow driven by the central WC-OB colliding-wind binary. An upper limit of $\lesssim120$ km s$^{-1}$ is derived for the expansion velocity assuming a distance of $4.15$ kpc. The upper limit on the average total mass-loss rate from the central 3" of WR~112 is estimated to be $\lesssim8\times10^{-6}$ $\mathrm{M}_\odot$ yr$^{-1}$. We leave its true nature as an open question, but propose that the WR~112 nebula may have formed in the outflow during a previous red or yellow supergiant phase of the central Wolf-Rayet star.

We report spatially resolved (FWHM$\sim3.8-4.6"$) mid-IR imaging observations of the planetary nebula (PN) NGC 7027 taken with the 2.5-m telescope aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). Images of NGC 7027 were acquired at 6.3, 6.6, 11.1, 19.7, 24.2, 33.6, and 37.1 $\mu\mathrm{m}$ using the Faint Object Infrared Camera for the SOFIA Telescope (FORCAST).The observations reveal emission from Polycyclic Aromatic Hydrocarbon (PAH) and warm dust ($T_D\sim90$ K) from the illuminated inner edge of the molecular envelope surrounding the ionized gas and central star. The DustEM code was used to fit the spectral energy distribution of fluxes obtained by FORCAST and the archival infrared spectrum of NGC 7027 acquired by the Short Wavelength Spectrometer (SWS) on the Infrared Space Observatory (ISO). Best-fit dust models provide a total dust mass of $5.8^{+2.3}_{-2.6}\times10^{-3}$ $\mathrm{M}_\odot$, where carbonaceous large ($a=1.5$ $\mu$m) and very small ($a \sim12\AA$) grains, and PAHs ($3.1\AA<a<12\AA$) compose 96.5, 2.2, and 1.3 $\%$ of the dust by mass, respectively. The 37 $\mu$m optical depth map shows minima in the dust column density at regions in the envelope that are coincident with a previously identified collimated outflow from the central star. The optical depth minima are also spatially coincident with enhancements in the 6.2 $\mu$m PAH feature, which is derived from the 6.3 and 6.6 $\mu$m maps. We interpret the spatial anti-correlation of the dust optical depth and PAH 6.2 $\mu$m feature strength and their alignment with the outflow from the central star as evidence of dust processing and rapid PAH formation via grain-grain collisions in the post-shock environment of the dense ($n_H\sim10^5\,\mathrm{cm}^{-3}$) photo-dissociation region (PDR) and molecular envelope.

While optical and radio transient surveys have enjoyed a renaissance over the past decade, the dynamic infrared sky remains virtually unexplored. The infrared is a powerful tool for probing transient events in dusty regions that have high optical extinction, and for detecting the coolest of stars that are bright only at these wavelengths. The fundamental roadblocks in studying the infrared time-domain have been the overwhelmingly bright sky background (250 times brighter than optical) and the narrow field-of-view of infrared cameras (largest is 0.6 sq deg). To begin to address these challenges and open a new observational window in the infrared, we present Palomar Gattini-IR: a 25 sq degree, 300mm aperture, infrared telescope at Palomar Observatory that surveys the entire accessible sky (20,000 sq deg) to a depth of 16.4 AB mag (J band, 1.25um) every night. Palomar Gattini-IR is wider in area than every existing infrared camera by more than a factor of 40 and is able to survey large areas of sky multiple times. We anticipate the potential for otherwise infeasible discoveries, including, for example, the elusive electromagnetic counterparts to gravitational wave detections. With dedicated hardware in hand, and a F/1.44 telescope available commercially and cost-effectively, Palomar Gattini-IR will be on-sky in early 2017 and will survey the entire accessible sky every night for two years. Palomar Gattini-IR will pave the way for a dual hemisphere, infrared-optimized, ultra-wide field high cadence machine called Turbo Gattini-IR. To take advantage of the low sky background at 2.5 um, two identical systems will be located at the polar sites of the South Pole, Antarctica and near Eureka on Ellesmere Island, Canada. Turbo Gattini-IR will survey 15,000 sq. degrees to a depth of 20AB, the same depth of the VISTA VHS survey, every 2 hours with a survey efficiency of 97%.

Since their discovery, the Quintuplet proper members (QPMs) have been somewhat mysterious in nature. Originally dubbed the "cocoon stars" due to their cool featureless spectra, high-resolution near-infrared imaging observations have shown that at least two of the objects exhibit "pinwheel" nebulae consistent with binary systems with a carbon-rich Wolf-Rayet star and O/B companion. In this paper, we present 19.7, 25.2, 31.5, and 37.1 $\mu$m observations of the QPMs (with an angular resolution of 3.2-3.8") taken with the Faint Object Infrared Camera for the SOFIA Telescope (FORCAST) in conjunction with high-resolution ($\sim$0.1-0.2") images at 8.8 and 11.7 $\mu$m from the Thermal-Region Camera Spectrograph (TReCS). DUSTY models of the thermal dust emission of two of the four detected QPMs, Q2 and Q3, are fitted by radial density profiles which are consistent with constant mass loss rates ($\rho_d \propto r^{-2}$). For the two remaining sources, Q1 and Q9, extended structures ($\sim$ 1") are detected around these objects in high-resolution imaging data. Based on the fitted dust masses, Q9 has an unusually large dust reservoir ($\mathrm{M_d}=1.3^{+0.8}_{-0.4}\times 10^{-3} \mathrm{M_{\odot}}$) compared to typical dusty Wolf-Rayet stars which suggests that it may have recently undergone an episode of enhanced mass loss.

We present multi-epoch mid-infrared (IR) photometry and the optical discovery observations of the "impostor" supernova (SN) 2010da in NGC 300 using new and archival Spitzer Space Telescope images and ground-based observatories. The mid-IR counterpart of SN 2010da was detected as SPIRITS 14bme in the SPitzer InfraRed Intensive Transient Survey (SPIRITS), an ongoing systematic search for IR transients. A sharp increase in the 3.6 $\mu$m flux followed by a rapid decrease measured ~150 d before and ~80 d after the initial outburst, respectively, reveal a mid-IR counterpart to the coincident optical and high luminosity X-ray outbursts. At late times after the outburst (~2000 d), the 3.6 and 4.5 $\mu$m emission increased to over a factor of 2 times the progenitor flux. We attribute the re-brightening mid-IR emission to continued dust production and increasing luminosity of the surviving system associated with SN 2010da. We analyze the evolution of the dust temperature, mass, luminosity, and equilibrium temperature radius in order to resolve the nature of SN 2010da. We address the leading interpretation of SN 2010da as an eruption from a luminous blue variable (LBV) high-mass X-ray binary (HMXB) system. We propose that SN 2010da is instead a supergiant (sg)B[e]-HMXB based on similar luminosities and dust masses exhibited by two other known sgB[e]-HMXB systems. Additionally, the SN 2010da progenitor occupies a similar region on a mid-IR color-magnitude diagram (CMD) with known sgB[e] stars in the Large Magellanic Cloud. The lower limit estimated for the orbital eccentricity of the sgB[e]-HMXB (e>0.82) from X-ray luminosity measurements is high compared to known sgHMXBs and supports the claim that SN 2010da may be associated with a newly formed HMXB system.

We present a systematic study of mid-infrared (mid-IR) emission from 141 nearby supernovae (SNe) observed with the InfraRed Array Camera (IRAC) on Spitzer.These SNe reside in one of the 190 galaxies within 20 Mpc drawn from the ongoing SPIRITS program. We detect 8 Type Ia SNe and 36 core-collapse SNe. All Type I SNe become undetectable within 3 years of explosion. About 22$\pm$11% of Type II SNe continue to be detected at late-times. Dust luminosity, temperature, and a lower liit on mass are obtained by fitting the SED using photometry with IRAC bands 1 and 2. The mass estimate does not distinguish between pre-existing and newly produced dust. We observe warm dust masses between $10^{-2}$ and $10^{-6}$ $\rm M_{\odot}$ and dust temperatures from 200 K to 1280 K.We present detailed case studies of two extreme Type II-P SNe: SN 2011ja and 2014bi. SN 2011ja was over-luminous ([4.5] = -15.6 mag) at 900 days post-explosion accompanied by the growing dust mass. This suggests either an episode of dust formation or an intensifying CSM interactions heating up pre-existing dust. SN 2014bi showed a factor of 10 decrease in dust mass over one month suggesting either an episode of dust destruction or a fading source of dust heating. A rebrightening of the Type Ib SN 2014C is observed and attributed to CSM interactions. This observation adds to a small number of stripped-envelope SNe that have mid-IR excess. The observations suggest that the CSM shell around SN 2014C is originated from an LBV-like eruption roughly 100 years before the explosion. We also report detections of SN1974E, 1979C, 1980K, 1986J, and 1993J detected more than 20 years post-explosion. The number of outlying SNe identified in this work demonstrates the power of late time mid-IR observations of a large sample of SNe to identify events with unusual evolution.

Massive, evolved stars play a crucial role in the metal-enrichment, dust budget, and energetics of the interstellar medium, however, the details of their evolution are uncertain because of their rarity and short lifetimes before exploding as supernovae. Discrepancies between theoretical predictions from single-star evolutionary models and observations of massive stars have evoked a shifting paradigm that implicates the importance of binary interaction. We present mid- to far-infrared observations from the Stratospheric Observatory for Infrared Astronomy (SOFIA) of a conical "helix" of warm dust ($\sim180$ K) that appears to extend from the Wolf-Rayet star WR102c. Our interpretation of the helix is a precessing, collimated outflow that emerged from WR102c during a previous evolutionary phase as a rapidly rotating luminous blue variable. We attribute the precession of WR102c to gravitational interactions with an unseen compact binary companion whose orbital period can be constrained to $800\,\mathrm{d}<P<1400$ d from the inferred precession period, $\tau_p\sim1.4\times10^4$ yr, and limits imposed on the stellar and orbital parameters of the system. Our results concur with the range of orbital periods ($P\lesssim1500$ d) where spin-up via mass exchange is expected to occur for massive binary systems.

We present mid-IR (19 - 37 microns) imaging observations of S106 from SOFIA/FORCAST, complemented with IR observations from Spitzer/IRAC (3.6 - 8.0 microns), IRTF/MIRLIN (11.3 and 12.5 microns), and Herschel/PACS (70 and 160 microns). We use these observations, observations in the literature, and radiation transfer modeling to study the heating and composition of the warm (~ 100 K) dust in the region. The dust is heated radiatively by the source S106 IR, with little contributions from grain-electron collisions and Ly-alpha radiation. The dust luminosity is >~ (9.02 +/- 1.01) x 10^4 L_sun, consistent with heating by a mid- to late-type O star. We find a temperature gradient (~ 75 - 107 K) in the lobes, which is consistent with a dusty equatorial geometry around S106 IR. Furthermore, the SOFIA observations resolve several cool (~ 65 - 70 K) lanes and pockets of warmer (~ 75 - 90 K) dust in the ionization shadow, indicating that the environment is fragmented. We model the dust mass as a composition of amorphous silicates, amorphous carbon, big grains, very small grains, and PAHs. We present the relative abundances of each grain component for several locations in S106.

Dust formation in supernova ejecta is currently the leading candidate to explain the large quantities of dust observed in the distant, early Universe. However, it is unclear whether the ejecta-formed dust can survive the hot interior of the supernova remnant (SNR). We present infrared observations of ~0.02 $M_\odot$ of warm (~100 K) dust seen near the center of the ~10,000 yr-old Sgr A East SNR at the Galactic center. Our findings signify the detection of dust within an older SNR that is expanding into a relatively dense surrounding medium ($n_e$ ~ 100 $\mathrm{cm}^{-3}$) and has survived the passage of the reverse shock. The results suggest that supernovae may indeed be the dominant dust production mechanism in the dense environment of early Universe galaxies.

We present imaging at 19, 25, 31, and 37 {\mu}m of the compact HII region complex G-0.02-0.07 located 6 pc in projection from the center of the Galaxy obtained with SOFIA using FORCAST. G-0.02-0.07 contains three compact HII regions (A, B, and C) and one ultra-compact HII region (D). Our observations reveal the presence of two faint, infrared sources located 23'' and 35'' to the east of region C (FIRS 1 and 2) and detect dust emission in two of the three "ridges" of ionized gas west of region A. The 19/37 color temperature and 37 {\mu}m optical depth maps of regions A - C are used to characterize the dust energetics and morphology. Regions A and B exhibit average 19/37 color temperatures of ~105 K, and regions C and D exhibit color temperatures of ~115 K and ~130 K, respectively. Using the DustEM code we model the SEDs of regions A - D and FIRS 1, all of which require populations of very small, transiently heated grains and large, equilibrium-heated grains. We also require the presence of polycyclic aromatic hydrocarbons (PAHs) in regions A - C in order to fit the 3.6, 4.5, 5.8, and 8.0 {\mu}m fluxes observed by Spitzer/IRAC. The location of the heating source for region A is determined by triangulation from distances and temperatures derived from DustEM models fit to SEDs of three different points around the region, and is found to be displaced to the northeast of the center of curvature near the color temperature peak. Based on total luminosity, expected 1.90 {\mu}m fluxes, and proximity to the mid-IR color temperature peaks we identify heating source candidates for regions A, B, and C. However, for region D, the observed fluxes at 1.87 and 1.90 {\mu}m of the previously proposed ionizing star are a factor of ~40 times too bright to be the heating source and hence is likely just a star lying along the line of sight towards region D.

Three Luminous Blue Variables (LBVs) are located in and near the Quintuplet Cluster at the Galactic Center: the Pistol star, G0.120-0.048, and qF362. We present imaging at 19, 25, 31, and 37 {\mu}m of the region containing these three LBVs, obtained with SOFIA using FORCAST. We argue that the Pistol and G0.120-0.048 are identical ``twins" that exhibit contrasting nebulae due to the external influence of their different environments. Our images reveal the asymmetric, compressed shell of hot dust surrounding the Pistol Star and provide the first detection of the thermal emission from the symmetric, hot dust envelope surrounding G0.120-0.048. Dust and gas composing the Pistol nebula are primarily heated and ionized by the nearby Quintuplet Cluster stars. The northern region of the Pistol nebula is decelerated due to the interaction with the high-velocity (2000 km/s) winds from adjacent Wolf-Rayet Carbon (WC) stars. With the DustEM code we determine that the Pistol nebula is composed of a distribution of very small, transiently-heated grains (10-~35 {\AA}) and that it exhibits a gradient of decreasing grain size from the south to the north due to differential sputtering by the winds from the WC stars. Dust in the G0.120-0.048 nebula is primarily heated by the central star; however, the nebular gas is ionized externally by the Arches Cluster. Unlike the Pistol nebula, the G0.120-0.048 nebula is freely expanding into the surrounding medium. Given independent dust and gas mass estimates we find that the Pistol and G0.120-0.048 nebulae exhibit similar gas-to-dust mass ratios of ~310 and ~290, respectively. Both nebulae share identical size scales (~ 0.7 pc) which suggests that they have similar dynamical timescales of ~10^5 yrs, assuming a shell expansion velocity of v_exp 60 km/s.

The temperature in the crust of an accreting neutron star, which comprises its outermost kilometer, is set by heating from nuclear reactions at large densities, neutrino cooling, and heat transport from the interior. The heated crust has been thought to affect observable phenomena at shallower depths, such as thermonuclear bursts in the accreted envelope. Here we report that cycles of electron capture and its inverse, $\beta^-$ decay, involving neutron-rich nuclei at a typical depth of about 150 m, cool the outer neutron star crust by emitting neutrinos while also thermally decoupling the surface layers from the deeper crust. This Urca mechanism has been studied in the context of white dwarfs and Type Ia supernovae, but hitherto was not considered in neutron stars, because previous models computed the crust reactions using a zero-temperature approximation and assumed that only a single nuclear species was present at any given depth. This thermal decoupling means that X-ray bursts and other surface phenomena are largely independent of the strength of deep crustal heating. The unexpectedly short recurrence times, of the order of years, observed for very energetic thermonuclear superbursts are therefore not an indicator of a hot crust, but may point instead to an unknown local heating mechanism near the neutron star surface.

We present 19.7, 31.5, and 37.1 {\mu}m images of the inner 6 pc of the Galactic Center of the Milky Way with a spatial resolution of 3.2 - 4.6'' taken by the Faint Object Infrared Camera on the Stratospheric Observatory for Infrared Astronomy (SOFIA). The images reveal in detail the "clumpy" structure of the Circumnuclear Ring (CNR)--the torus of hot gas and dust orbiting the supermassive black hole at the Galactic Center with an inner radius of 1.4 pc. The CNR exhibits features of a classic HII region: the dust emission at 19.7 {\mu}m closely traces the ionized gas emission observed in the radio while the 31.5 and 37.1 {\mu}m emission traces the photo-dissociation region beyond the ionized gas. The 19.7/37.1 color temperature map reveals a radial temperature gradient across the CNR with temperatures ranging from 65-85 K, consistent with the prevailing paradigm in which the dust is centrally heated by the inner cluster of hot, young stars. We produce a 37.1 {\mu}m intensity model of the CNR with the derived geometric properties and find that it is consistent with the observed 37.1 {\mu}m map of the CNR. Dense ($5 to 9 \times 10^{4} \, \mathrm{cm}^{-3}$) clumps with a FWHM of ~0.15 pc exist along the inner edge of the CNR and shadow the material deeper into the ring. The clumps are unlikely to be long-lived structures since they are not dense enough to be stable against tidal shear from the supermassive black hole.

The location of electron capture heat sources in the crust of accreting neutron stars depends on the masses of extremely neutron-rich nuclei. We present first results from a new implementation of the time-of-flight technique to measure nuclear masses of rare isotopes at the National Superconducting Cyclotron Laboratory. The masses of 16 neutron-rich nuclei in the scandium -- nickel range were determined simultaneously, improving the accuracy compared to previous data in 12 cases. The masses of $^{61}${V}, $^{63}${Cr}, $^{66}${Mn}, and $^{74}${Ni} were measured for the first time with mass excesses of $-30.510(890)$ MeV, $-35.280(650)$ MeV, $-36.900(790)$ MeV, and $-49.210(990)$ MeV, respectively. With the measurement of the $^{66}$Mn mass, the locations of the two dominant electron capture heat sources in the outer crust of accreting neutron stars that exhibit superbursts are now experimentally constrained. We find that the location of the $^{66}$Fe$\rightarrow^{66}$Mn electron capture transition occurs significantly closer to the surface than previously assumed because our new experimental Q-value is 2.1 MeV (2.6$\sigma$) smaller than predicted by the FRDM mass model.