We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $\mu$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $\mu$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($<1.6$ kpc), the degree of polarization is found to be lower than for far-infrared sources in the Milky Way, and has minima at the locations of most intense star formation near the outer ends of the inner-bar. Inside the starburst ring, the field direction deviates from the model, becoming more radial along the leading edges of the inner-bar. The polarized flux and dust temperature peak $\sim 3-6$" NE of the AGN at the location of a bow shock between the AGN outflow and the surrounding interstellar medium, but the AGN itself is weakly polarized ($< 1$%) at both 53 and 89 \um.

We report on polarimetric maps made with HAWC+/SOFIA toward Rho Oph A, the densest portion of the Rho Ophiuchi molecular complex. We employed HAWC+ bands C (89 $\mu$m) and D (154 $\mu$m). The slope of the polarization spectrum was investigated by defining the quantity R_DC = p_D/p_C, where p_C and p_D represent polarization degrees in bands C and D, respectively. We find a clear correlation between R_DC and the molecular hydrogen column density across the cloud. A positive slope (R_DC > 1) dominates the lower density and well illuminated portions of the cloud, that are heated by the high mass star Oph S1, whereas a transition to a negative slope (R_DC < 1) is observed toward the denser and less evenly illuminated cloud core. We interpret the trends as due to a combination of: (1) Warm grains at the cloud outskirts, which are efficiently aligned by the abundant exposure to radiation from Oph S1, as proposed in the radiative torques theory; and (2) Cold grains deep in the cloud core, which are poorly aligned due to shielding from external radiation. To assess this interpretation, we developed a very simple toy model using a spherically symmetric cloud core based on Herschel data, and verified that the predicted variation of R_DC is consistent with the observations. This result introduces a new method that can be used to probe the grain alignment efficiency in molecular clouds, based on the analysis of trends in the far-infrared polarization spectrum.

We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions.

We present Far-Infrared polarimetry observations of M82 at 53 and $154~\mu \rm{m}$ and NGC 253 at $89~\mu \rm{m}$, which were taken with HAWC+ in polarimetry mode on the Stratospheric Observatory for Infrared Astronomy (SOFIA). The polarization of M82 at $53~\mu \rm{m}$ clearly shows a magnetic field geometry perpendicular to the disk in the hot dust emission. For M82 the polarization at $154~\mu \rm{m}$ shows a combination of field geometry perpendicular to the disk in the nuclear region, but closer to parallel to the disk away from the nucleus. The fractional polarization at $53~\mu \rm{m}$ $(154~\mu \rm{m})$ ranges from 7% (3%) off nucleus to 0.5% (0.3%) near the nucleus. A simple interpretation of the observations of M82 invokes a massive polar outflow, dragging the field along, from a region $\sim 700$~pc in diameter that has entrained some of the gas and dust, creating a vertical field geometry seen mostly in the hotter $(53~\mu \rm{m})$ dust emission. This outflow sits within a larger disk with a more typical planar geometry that more strongly contributes to the cooler $(154~\mu \rm{m})$ dust emission. For NGC 253, the polarization at $89~\mu \rm{m}$ is dominated by a planar geometry in the tilted disk, with weak indication of a vertical geometry above and below the plane from the nucleus. The polarization observations of NGC 253 at $53~\mu \rm{m}$ were of insufficient S/N for detailed analysis.

We present observations of NGC 1068 covering the $19.7-53.0$ $\mu$m wavelength range using FORCAST and HAWC+ onboard SOFIA. Using these observations, high-angular resolution infrared (IR) and sub-mm observations, we find an observational turn-over of the torus emission in the $30-40$ $\mu$m wavelength range with a characteristic temperature of $70-100$ K. This component is clearly different from the diffuse extended emission in the narrow line and star formation regions at 10-100 $\mu$m within the central 700 pc. We compute $2.2-432$ $\mu$m 2D images using the best inferred \textsc{clumpy} torus model based on several nuclear spectral energy distribution (SED) coverages. We find that when $1-20$ $\mu$m SED is used, the inferred result gives a small torus size ($<4$ pc radius) and a steep radial dust distribution. The computed torus using the $1-432$ $\mu$m SED provides comparable torus sizes, $5.1^{+0.4}_{-0.4}$ pc radius, and morphology to the recently resolved 432 $\mu$m ALMA observations. This result indicates that the $1-20$ $\mu$m wavelength range is not able to probe the full extent of the torus. The characterization of the turn-over emission of the torus using the $30-60$ $\mu$m wavelength range is sensitive to the detection of cold dust in the torus. The morphology of the dust emission in our 2D image at 432 $\mu$m is spatially coincident with the cloud distribution, while the morphology of the emission in the $1-20$ $\mu$m wavelength range shows an elongated morphology perpendicular to the cloud distribution. We find that our 2D \textsc{clumpy} torus image at 12 $\mu$m can produce comparable results to those observed using IR interferometry.

The benchmark exoplanet GJ 1214b is one of the best studied transiting planets in the transition zone between rocky Earth-sized planets and gas or ice giants. This class of super-Earth/mini-Neptune planets is unknown in our Solar System, yet is one of the most frequently detected classes of exoplanets. Understanding the transition from rocky to gaseous planets is a crucial step in the exploration of extrasolar planetary systems, in particular with regard to the potential habitability of this class of planets. GJ 1214b has already been studied in detail from various platforms at many different wavelengths. Our airborne observations with SOFIA add information in the Paschen-alpha cont. 1.9 micron infrared wavelength band, which is not accessible by any other current ground- or space-based instrument due to telluric absorption or limited spectral coverage. We used FLIPO and FPI+ on SOFIA to comprehensively analyse the transmission signal of the possible water-world GJ 1214b through photometric observations during transit in three optical and one infrared channels. We present four simultaneous light curves and corresponding transit depths in three optical and one infrared channel, which we compare to previous observations and state-of-the-art synthetic atmospheric models of GJ 1214b. The final precision in transit depth is between 1.5 and 2.5 times the theoretical photon noise limit, not sensitive enough to constrain the theoretical models any better than previous observations. This is the first exoplanet observation with SOFIA that uses its full set of instruments available to exoplanet spectrophotometry. Therefore we use these results to evaluate SOFIA's potential in this field and suggest future improvements.

We use a modified version of MOOG to generate large grids of synthetic spectra in an attempt to derive quantitative abundances for three CVs (GK Per, RU Peg, and SS Cyg) by comparing the models to moderate resolution (R $\sim$ 25,000) $K$-band spectra obtained with NIRSPEC on Keck. For each of the three systems we find solar, or slightly sub-solar values for [Fe/H], but significant deficits of carbon: for SS Cyg we find [C/Fe] = $-$0.50, for RU Peg [C/Fe] = $-$0.75, and for GK Per [C/Fe] = $-$1.00. We show that it is possible to use lower resolution (R $\sim$ 2,000) spectra to quantify carbon deficits. We examine realistic veiling scenarios and find that emission from H I or CO cannot reproduce the observations.

Here we report on the first successful exoplanet transit observation with the Stratospheric Observatory for Infrared Astronomy (SOFIA). We observed a single transit of the hot Jupiter HD 189733 b, obtaining two simultaneous primary transit lightcurves in the B and z' bands as a demonstration of SOFIA's capability to perform absolute transit photometry. We present a detailed description of our data reduction, in particular the correlation of photometric systematics with various in-flight parameters unique to the airborne observing environment. The derived transit depths at B and z' wavelengths confirm a previously reported slope in the optical transmission spectrum of HD 189733 b. Our results give new insights to the current discussion about the source of this Rayleigh scattering in the upper atmosphere and the question of fixed limb darkening coefficients in fitting routines.

We present medium resolution near-infrared (NIR) spectra, covering 1.1 to 3.4 microns, of the normal Type Ia supernova (SN Ia) SN 2014J in M82 obtained with the FLITECAM instrument aboard SOFIA approximately 17-25 days after maximum B light. Our 2.8-3.4 micron spectra may be the first ~3 micron spectra of a SN Ia ever published. The spectra spanning the 1.5-2.7 micron range are characterized by a strong emission feature at ~1.77 microns with a full width at half maximum of ~11,000-13,000 km/s. We compare the observed FLITECAM spectra to the recent non-LTE delayed detonation models of Dessart et al. (2014) and find that the models agree with the spectra remarkably well in the 1.5-2.7 micron wavelength range. Based on this comparison we identify the ~1.77 micron emission peak as a blend of permitted lines of Co II. Other features seen in the 2.0 - 2.5 micron spectra are also identified as emission from permitted transitions of Co II. However, the models are not as successful at reproducing the spectra in the 1.1 - 1.4 micron range or between 2.8 and 3.4 microns. These observations demonstrate the promise of SOFIA by allowing access to wavelength regions inaccessible from the ground, and serve to draw attention to the usefulness of the regions between the standard ground-based NIR passbands for constraining SN models.

We have used the PACS instrument on the Herschel Space Observatory to observe eight cataclysmic variables at 70 and 160 microns. Of these eight objects, only AM Her was detected. We have combined the Herschel results with ground-based, Spitzer, and WISE observations to construct spectral energy distributions for all of the targets. For the two dwarf novae in the sample, SS Cyg and U Gem, we find that their infrared luminosities are completely dominated by their secondary stars. For the two highly magnetic "polars" in our survey, AM Her and EF Eri, we find that their mid-infrared excesses, previously attributed to circumbinary dust emission, can be fully explained by cyclotron emission. The WISE light curves for both sources show large, orbitally modulated variations that are identically phased to their near-IR light curves. We propose that significant emission from the lowest cyclotron harmonics (n </= 3) is present in EF Eri and AM Her. Previously, such emission would have been presumed to be optically thick, and not provide significant orbitally modulated flux. This suggests that the accretion onto polars is more complicated than assumed in the simple models developed for these two sources. We develop a model for the near-/mid-IR light curves for WZ Sge with an L2 donor star that shows that the ellipsoidal variations from its secondary star are detected. We conclude that none of the targets surveyed have dusty circumbinary disks.

We present K-band spectroscopy of short period, "sub-gap" cataclysmic variable (CV) systems obtained using ISAAC on the VLT. We show the infrared spectra (IR) for nine systems below the 2-3 hour period gap: V2051 Oph, V436 Cen, EX Hya, VW Hyi, Z Cha, WX Hyi, V893 Sco, RZ Leo, and TY PsA. We are able to clearly detect the secondary star in all but WX Hyi, V893 Sco, and TY PsA. We present the first direct detection of the secondary stars of V2051 Oph, V436 Cen, and determine new spectral classifications for EX Hya, VW Hyi, Z Cha, and RZ Leo. We find that the CO band strengths of all but Z Cha appear normal for their spectral types, in contrast to their longer period cousins above the period gap. This brings the total number of CVs and pre-CVs with moderate resolution (R >~ 1500) IR spectroscopy to sixty-one systems: nineteen pre-CVs, thirty-one non-magnetic systems, and eleven magnetic or partially magnetic systems. We discuss the trends seen in the IR abundance patterns thus far, and highlight a potential link between anomalous abundances seen in the IR with the C IV/N V anomaly seen in the ultraviolet. We present a compilation of all systems with sufficient resolution IR observations to assess the CO band strengths, and, by proxy, obtain an estimate on the C abundance on the secondary star.

Abriged version for astroph: The young late-type star V1118 Orionis was in outburst from 2005 to 2006. We followed the outburst with optical and near-infrared photometry; the X-ray emission was further probed with observations taken with XMM-Newton and Chandra during and after the outburst. In addition, we obtained mid-infrared photometry and spectroscopy with Spitzer at the peak of the outburst and in the post-outburst phase. The spectral energy distribution of V1118 Ori varied significantly over the course of the outburst. The optical flux showed the largest variations, most likely due to enhanced emission by a hot spot. The latter dominated the optical and near-infrared emission at the peak of the outburst, while the disk emission dominated in the mid-infrared. The X-ray flux correlated with the optical and infrared fluxes, indicating that accretion affected the magnetically active corona and the stellar magnetosphere. The thermal structure of the corona was variable with some indication of a cooling of the coronal temperature in the early phase of the outburst with a gradual return to normal values. Color-color diagrams in the optical and infrared showed variations during the outburst, with no obvious signature of reddening due to circumstellar matter. Using MC realizations of star+disk+hotspot models to fit the SED in ``quiescence'' and at the peak of the outburst, we determined that the mass accretion rate varied from about 2.5E-7 Msun/yr to 1E-6 Msun/yr; in addition the fractional area of the hotspot increased significantly as well. The multi-wavelength study of the V1118 Ori outburst helped us to understand the variations in spectral energy distributions and demonstrated the interplay between the disk and the stellar magnetosphere in a young, strongly accreting star.

The prototype nova-like variable VY Sculptoris was observed by the IUE during four different optical brightness states of the system. The FUV flux level from the highest state to the lowest state declines by a factor of 28. We have carried out model accretion disk and white dwarf atmosphere fitting to the spectra. The corresponding accretion rates range from $\dot{M}= 8 \times 10^{-9}$M$_{\sun}$ yr$^{-1}$ at the highest FUV flux level down to $\dot{M}= 1.9\times 10^{-10}$M$_{\sun}$ yr$^{-1}$ at the lowest flux level. We report tentative evidence for the detection of the underlying white dwarf with $T_{\rm eff} = 45,000$K in the spectrum with the lowest flux level.

We present a synthetic spectral analysis of nearly the entire far ultraviolet International Ultraviolet Explorer (IUE) archive of spectra of dwarf novae in or near outburst. The study includes 46 systems of all dwarf nova subtypes both above and below the period gap. The spectra were uniformly analyzed using synthetic spectral codes for optically thick accretion disks and stellar photospheres along with the best-available distance measurements or estimates. We present newly estimated accretion rates and discuss the implications of our study for disk accretion physics and CV evolution.

We present results from a multi-wavelength campaign to monitor the 2005 outburst of the low-mass young star V1118 Ori. Although our campaign covers the X-ray, optical, infrared, and radio regimes, we focus in this Letter on the properties of the X-ray emission in V1118 Ori during the first few months after the optical outburst. Chandra and XMM-Newton detected V1118 Ori at three epochs in early 2005. The X-ray flux and luminosity stayed similar within a factor of two, and at the same level as in a pre-outburst observation in 2002. The hydrogen column density showed no evidence for variation from its modest pre-outburst value of $N_\mathrm{H} \sim 3 \times 10^{21}$ cm$^{-2}$. However, a spectral change occurred from a dominant hot plasma ($\sim 25$ MK) in 2002 and in January 2005 to a cooler plasma ($\sim 8$ MK) in February 2005 and in March 2005. We hypothesize that the hot magnetic loops high in the corona were disrupted by the closing in of the accretion disk due to the increased accretion rate during the outburst, whereas the lower cooler loops were probably less affected and became the dominant coronal component.

Using newly determined parallaxes for dwarf novae, we have derived outburst accretion rates for VY Aqr, RU Peg and T Leo and for T Leo during quiescence. The two short-period dwarf novae, VY Aqr and T Leo, show good agreement with optically thick steady-state accretion disks in outburst, whereas RU Peg shows a significant departure from a steady-state disk. We have determined that the white dwarf in T Leo has T$_{eff} = 16,000 \pm 1000$ K, a value consistent with long term compressional heating when gravitational wave emission drives mass transfer. The white dwarf in T Leo has a temperature in the same narrow range as other WZ Sge-like dwarf novae.