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.

From our Secret Lives of Cepheids program, the prototype Classical Cepheid, $\delta$ Cep, is found to be an X-ray source with periodic pulsation-modulated X-ray variations. This finding complements our earlier reported phase-dependent FUV--UV emissions of the star that increase $\sim$10--20 times with highest fluxes at $\sim0.90-0.95\phi$, just prior to maximum brightness. Previously, $\delta$ Cep was found as potentially X-ray variable, using XMM-Newton observations (Engle et al. 2014). Additional phase-constrained data were secured with Chandra near X-ray emission peak, to determine if the emission and variability were pulsation-phase-specific to $\delta$ Cep and not transient or due to a possible coronally-active, cool companion. The Chandra data were combined with prior XMM-Newton observations, and very closely match the previously observed X-ray behavior. From the combined dataset, a $\sim$4$\times$ increase in X-ray flux is measured, reaching a peak $L_X$ = 1.7 $\times$ 10$^{29}$ erg s$^{-1}$ near 0.45$\phi$. The precise X-ray flux phasing with the star's pulsation indicates that the emissions arise from the Cepheid and not a companion. However, it is puzzling that maximum X-ray flux occurs $\sim$0.5$\phi$ ($\sim$3 days) later than the FUV--UV maximum. There are several other potential Cepheid X-ray detections with properties similar to $\delta$ Cep, and comparable X-ray variability is indicated for two other Cepheids: $\beta$ Dor and V473 Lyr. X-ray generating mechanisms in $\delta$ Cep and other Cepheids are discussed. If additional Cepheids are confirmed to show phased X-ray variations, then $\delta$ Cep will be the prototype of new class of pulsation-induced X-ray variables.

The Stratospheric Observatory for Infrared Astronomy (SOFIA) is an airborne observatory consisting of a specially modified Boeing 747SP with a 2.7-m telescope, flying at altitudes as high as 13.7 km (45,000 ft). Designed to observe at wavelengths from 0.3 micron to 1.6 mm, SOFIA operates above 99.8 % of the water vapor that obscures much of the infrared and submillimeter. SOFIA has seven science instruments under development, including an occultation photometer, near-, mid-, and far-infrared cameras, infrared spectrometers, and heterodyne receivers. SOFIA, a joint project between NASA and the German Aerospace Center DLR, began initial science flights in 2010 December, and has conducted 30 science flights in the subsequent year. During this early science period three instruments have flown: the mid-infrared camera FORCAST, the heterodyne spectrometer GREAT, and the occultation photometer HIPO. This article provides an overview of the observatory and its early performance.

This white paper is submitted to the Astronomy and Astrophysics 2010 Decadal Survey (Astro2010)1 Committee on the State of the Profession to emphasize the potential of the Stratospheric Observatory for Infrared Astronomy (SOFIA) to contribute to the training of instrumentalists and observers, and to related technology developments. This potential goes beyond the primary mission of SOFIA, which is to carry out unique, high priority astronomical research. SOFIA is a Boeing 747SP aircraft with a 2.5 meter telescope. It will enable astronomical observations anywhere, any time, and at most wavelengths between 0.3 microns and 1.6 mm not accessible from ground-based observatories. These attributes, accruing from the mobility and flight altitude of SOFIA, guarantee a wealth of scientific return. Its instrument teams (nine in the first generation) and guest investigators will do suborbital astronomy in a shirt-sleeve environment. The project will invest $10M per year in science instrument development over a lifetime of 20 years. This, frequent flight opportunities, and operation that enables rapid changes of science instruments and hands-on in-flight access to the instruments, assure a unique and extensive potential - both for training young instrumentalists and for encouraging and deploying nascent technologies. Novel instruments covering optical, infrared, and submillimeter bands can be developed for and tested on SOFIA by their developers (including apprentices) for their own observations and for those of guest observers, to validate technologies and maximize observational effectiveness.

Multi-wavelength imaging polarimetry at far-infrared wavelengths has proven to be an excellent tool for studying the physical properties of dust, molecular clouds, and magnetic fields in the interstellar medium. Although these wavelengths are only observable from airborne or space-based platforms, no first-generation instrument for the Stratospheric Observatory for Infrared Astronomy (SOFIA) is presently designed with polarimetric capabilities. We study several options for upgrading the High-resolution Airborne Wideband Camera (HAWC) to a sensitive FIR polarimeter. HAWC is a 12 x 32 pixel bolometer camera designed to cover the 53 - 215 micron spectral range in 4 colors, all at diffraction-limited resolution (5 - 21 arcsec). Upgrade options include: (1) an external set of optics which modulates the polarization state of the incoming radiation before entering the cryostat window; (2) internal polarizing optics; and (3) a replacement of the current detector array with two state-of-the-art superconducting bolometer arrays, an upgrade of the HAWC camera as well as polarimeter. We discuss a range of science studies which will be possible with these upgrades including magnetic fields in star-forming regions and galaxies and the wavelength-dependence of polarization.

We have performed deep imaging of a diverse sample of 26 low surface brightness galaxies (LSBGs) in the optical and the near-infrared. Using stellar population synthesis models, we find that it is possible to place constraints on the ratio of young to old stars (which we parameterise in terms of the average age of the galaxy), as well as the metallicity of the galaxy, using optical and near-infrared colours. LSBGs have a wide range of morphologies and stellar populations, ranging from older, high metallicity earlier types to much younger and lower metallicity late type galaxies. Despite this wide range of star formation histories, we find that colour gradients are common in LSBGs. These are most naturally interpreted as gradients in mean stellar age, with the outer regions of LSBGs having younger ages than their inner regions. In an attempt to understand what drives the differences in LSBG stellar populations, we compare LSBG average ages and metallicities with their physical parameters. Strong correlations are seen between a LSBG's star formation history and its K band surface brightness, K band absolute magnitude and gas fraction. These correlations are consistent with a scenario in which the star formation history of a LSBG primarily correlates with its surface density and its metallicity correlates both with its mass and surface density.