The Galactic magnetar SGR1935+2154 has been reported to produce the first known example of a bright millisecond duration radio burst (FRB 200428) similar to the cosmological population of fast radio bursts (FRBs), bolstering the association of FRBs to active magnetars. The detection of a coincident bright X-ray burst has revealed the first observed multi-wavelength counterpart of a FRB. However, the search for similar emission at optical wavelengths has been hampered by the high inferred extinction on the line of sight. Here, we present results from the first search for second-timescale emission from the source at near-infrared wavelengths using the Palomar Gattini-IR observing system in J-band, made possible by a recently implemented detector read-out mode that allowed for short exposure times of 0.84 s with 99.9% observing efficiency. With a total observing time of 12 hours (47728 images) on source, we place median $3\,\sigma$ limits on the second-timescale emission of $< 20$ mJy (13.1 AB mag). We present non-detection limits from epochs of four simultaneous X-ray bursts detected by the Insight-{\it HXMT} and {\it NuSTAR} telescopes during our observing campaign. The limits translate to an extinction corrected fluence limit of $< 125$ Jy ms for an estimated extinction of $A_J = 2.0$ mag. These limits provide the most stringent constraints to date on the fluence of flares at frequencies of $\sim 10^{14}$ Hz, and constrain the ratio of the near-infrared (NIR) fluence to that of coincident X-ray bursts to $R_{\rm NIR} < 2.5 \times 10^{-2}$. Our observations were sensitive enough to easily detect a near-infrared counterpart of FRB 200428 if the NIR emission falls on the same power law as that observed across its radio to X-ray spectrum. The non-detection of NIR emission around the coincident X-ray bursts constrains the fluence index of the brightest burst to be steeper than $0.35$.

Helium is expected to be present in the ejecta of some hydrogen-poor superluminous supernovae (SLSN-I). However, so far only one event, PTF10hgi has been identified with He features in its photospheric spectra (Quimby et al. 2018). At $z=0.0866$, ZTF19aawfbtg (SN2019hge) has more than 10 optical spectra at phases from $-41$ to $+103$ days. Most of these spectra match well with that of PTF10hgi, suggesting SN2019hge as the second case of a He-rich SLSN-I. Confirmation comes from its near-IR spectrum taken at $+34$ days, revealing He I features with P-Cygni profiles at 1.083 and 2.058$\mu$m. Using the optical spectra of PTF10hgi and SN2019hge as templates, we examine 70 SLSN-I discovered by ZTF in the first two years of operation and found additional five SLSN-I with distinct He-features, yielding a lower limit to the He-rich fraction of $\sim9$\%. It is suggested that the excitation of He I atoms requires non-thermal radiation. These He-rich events can not be explained by traditional $^{56}$Ni mixing model because of their blue spectra, high peak luminosity and slow rising time scales. Magnetar models offer a possible solution since pulsar winds naturally generate high energy particles as sources of non-thermal excitation. However, how the magnetar energy is transported to the outer layers is unclear. An alternative model is ejecta interaction with H-poor CSM which may be supported by the observed light curve undulations. These six SLSN-Ib appear to have relatively low-peak luminosities, around $-20.06\pm0.16$ mag ($g$-band, rest-frame), in comparison to the full range of $-20$ to $-23$ mag covered by SLSN-I.

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

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).

(Abridged) Using the Zwicky Transient Facility alert stream, we are conducting a large campaign to spectroscopically classify all transients occurring in galaxies in the Census of the Local Universe (CLU) catalog. The aim of the experiment is to construct a spectroscopically complete, volume-limited sample of transients coincident within 100" of CLU galaxies out to 200 Mpc, and to a depth of 20 mag. We describe the survey design and spectroscopic completeness from the first 16 months of operations. We present results from a systematic search for Calcium rich gap transients in the sample of 22 low luminosity (peak absolute magnitude $M > -17$), hydrogen poor events found in the experiment (out of 754 spectroscopically classified SNe). We report the detection of eight Calcium rich gap transients, and constrain their volumetric rate to be at least $\approx 15\pm5$% of the SN Ia rate. Combining this sample with ten events from the literature, we find a likely continuum of spectroscopic properties ranging from events with SN Ia-like features (Ca-Ia objects) to SN Ib/c-like features (Ca-Ib/c objects) at peak light. Within the Ca-Ib/c events, we find two populations of events distinguished by their red ($g - r \approx 1.5$ mag) or green ($g - r \approx 0.5$ mag) spectral colors at $r$-band peak, wherein redder events show strong line blanketing signatures, slower light curves, weaker He lines and lower [Ca II]/[O I] in the nebular phase. Together, we find that the spectroscopic continuum, volumetric rates and striking old environments are consistent with the explosive burning of He shells on low mass white dwarfs. We posit that Ca-Ia and red Ca-Ib/c objects are consistent with the double detonation of He shells with high He burning efficiency, while green Ca-Ib/c objects could arise from less efficient He burning scenarios such as detonations in low density He shells or He shell deflagrations.

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.

We present a survey of the early evolution of 12 Type IIn supernovae (SNe IIn) in the Ultra-Violet (UV) and visible light. We use this survey to constrain the geometry of the circumstellar material (CSM) surrounding SN IIn explosions, which may shed light on their progenitor diversity. In order to distinguish between aspherical and spherical circumstellar material (CSM), we estimate the blackbody radius temporal evolution of the SNe IIn of our sample, following the method introduced by Soumagnac et al. We find that higher luminosity objects tend to show evidence for aspherical CSM. Depending on whether this correlation is due to physical reasons or to some selection bias, we derive a lower limit between 35% and 66% on the fraction of SNe IIn showing evidence for aspherical CSM. This result suggests that asphericity of the CSM surrounding SNe IIn is common - consistent with data from resolved images of stars undergoing considerable mass loss. It should be taken into account for more realistic modelling of these events.

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

(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.

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

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.

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 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.

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.

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 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.

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.

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.