Search

We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $\mu$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation and sample selection of the program. Here, we focus on the general trends in the emissive polarization fraction. Far-infrared polarimetric observations trace the thermal polarized emission of magnetically aligned dust grains across the galaxy disks with polarization fractions of $P=0-15$% in the cold, $T_{\rm d} = [19,48]$ K, and dense, $\log_{10}(N_{\rm HI+H_{2}}) = [19.96,22.91]$, interstellar medium. The spiral galaxies show a median $\langle P_{154\mu m} \rangle = 3.3\pm0.9 $% across the disks. We report the first polarized spectrum of starburst galaxies showing a minimum within $89-154$ $\mu$m. The falling $53-154$ $\mu$m polarized spectrum may be due to a decrease in the dust grain alignment efficiency produced by variations in dust temperatures along the line-of-sight in the galactic outflow. We find that the starburst galaxies and the star-forming regions within normal galaxies have the lowest polarization fractions. We find that 50% (7 out of 14) of the galaxies require a broken power-law in the $P-N_{HI+H_{2}}$ and $P-T_{d}$ relations with three different trends. Group 1 has a relative increase of anisotropic random B-fields produced by compression or shear of B-fields in the galactic outflows, starburst rings, and inner-bar of galaxies; and Groups 2 and 3 have a relative increase of isotropic random B-fields driven by star-forming regions in the spiral arms, and/or an increase of dust grain alignment efficiency caused by shock-driven regions or evolutionary stages of a galaxy., Comment: 43 pages, 28 figures. Accepted for publication to ApJ

We describe the data processing of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program). This first data release presents 33% (51.34h out of 155.7h, including overheads) of the total awarded time taken from January 2020 to December 2021. Our observations were performed using the newly implemented on-the-fly mapping (OTFMAP) technique in the polarimetric mode. We present the pipeline steps to obtain homogeneously reduced high-level data products of polarimetric maps of galaxies for use in scientific analysis. Our approach has a general design and can be applied to sources smaller than the field-of-view of the HAWC+ array in any given band. We estimate that the OTFMAP polarimetric mode offers a reduction of observing overheads by a factor 2.34, and an improvement in sensitivity by a factor 1.80 when compared to previously obtained polarimetric observations using the chopping and nodding mode. The OTFMAP is a significant optimization of the polarimetric mode of HAWC+ as it ultimately reduces the cost of operations of SOFIA/HAWC+ by increasing the science collected per hour of observation up to an overall factor of 2.49. The OTFMAP polarimetric mode is the standard observing strategy of SALSA. The results and quantitative analysis of this first data release are presented in Papers IV and V of the series., Comment: 26 pages, 16 figures, 5 tables. Submitted to ApJ (Comments welcome)

Applying dendrogram analysis to the CARMA-NRO C$^{18}$O ($J$=1--0) data having an angular resolution of $\sim$ 8", we identified 692 dense cores in the Orion Nebula Cluster (ONC) region. Using this core sample, we compare the core and initial stellar mass functions in the same area to quantify the step from cores to stars. About 22 \% of the identified cores are gravitationally bound. The derived core mass function (CMF) for starless cores has a slope similar to Salpeter's stellar initial mass function (IMF) for the mass range above 1 $M_\odot$, consistent with previous studies. Our CMF has a peak at a subsolar mass of $\sim$ 0.1 $M_\odot$, which is comparable to the peak mass of the IMF derived in the same area. We also find that the current star formation rate is consistent with the picture in which stars are born only from self-gravitating starless cores. However, the cores must gain additional gas from the surroundings to reproduce the current IMF (e.g., its slope and peak mass), because the core mass cannot be accreted onto the star with a 100\% efficiency. Thus, the mass accretion from the surroundings may play a crucial role in determining the final stellar masses of stars., Comment: 13 pages, 5 figures, accepted by ApJL

SOFIA HAWC+ polarimetry at $154~\micron$ is reported for the face-on galaxy M51 and the edge-on galaxy NGC 891. For M51, the polarization vectors generally follow the spiral pattern defined by the molecular gas distribution, the far-infrared (FIR) intensity contours, and other tracers of star formation. The fractional polarization is much lower in the FIR-bright central regions than in the outer regions, and we rule out loss of grain alignment and variations in magnetic field strength as causes. When compared with existing synchrotron observations, which sample different regions with different weighting, we find the net position angles are strongly correlated, the fractional polarizations are moderately correlated, but the polarized intensities are uncorrelated. We argue that the low fractional polarization in the central regions must be due to significant numbers of highly turbulent segments across the beam and along lines of sight in the beam in the central 3 kpc of M51. For NGC 891, the FIR polarization vectors within an intensity contour of 1500 $\rm{MJy~sr^{-1}}$ are oriented very close to the plane of the galaxy. The FIR polarimetry is probably sampling the magnetic field geometry in NGC 891 much deeper into the disk than is possible with NIR polarimetry and radio synchrotron measurements. In some locations in NGC 891 the FIR polarization is very low, suggesting we are preferentially viewing the magnetic field mostly along the line of sight, down the length of embedded spiral arms. There is tentative evidence for a vertical field in the polarized emission off the plane of the disk., Comment: AJ, Accepted, In Press

We analyze the wavelength dependence of the far-infrared polarization fraction toward the OMC-1 star forming region using observations from HAWC+/SOFIA at 53, 89, 154, and 214 $\mu$m. We find that the shape of the far-infrared polarization spectrum is variable across the cloud and that there is evidence of a correlation between the slope of the polarization spectrum and the average line-of-sight temperature. The slope of the polarization spectrum tends to be negative (falling toward longer wavelengths) in cooler regions and positive or flat in warmer regions. This is very similar to what was discovered in $\rho$ Oph A via SOFIA polarimetry at 89 and 154 $\mu$m. Like the authors of this earlier work, we argue that the most natural explanation for our falling spectra is line-of-sight superposition of differing grain populations, with polarized emission from the warmer regions and less-polarized emission from the cooler ones. In contrast with the earlier work on $\rho$ Oph A, we do not find a clear correlation of polarization spectrum slope with column density. This suggests that falling spectra are attributable to variations in grain alignment efficiency in a heterogeneous cloud consistent with radiative torques theory. Alternative explanations in which variations in grain alignment efficiency are caused by varying gas density rather than by varying radiation intensity are disfavored., Comment: 21 pages, 10 figures. Accepted to ApJ. Major changes from first draft: 1. Additional discussion/introduction about ETAC/HCE, 2. More in-depth analysis on robustness of polarization slope/dust temperature correlation, 3. Updated notation on the spatial polarization ratio analysis, 4. Updated many figures for clarity and consistency

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., Comment: Accepted by ApJ (July 18 2019)

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., Comment: 6 Figures, 20 pages, ApJL accepted

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., Comment: 28 pages, 14 figures, ApJ, accepted

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., Comment: 12 pages, 7 Figures, Accepted for publication in ApJ

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., Comment: 10 pages, no figures, White Paper for Astro 2010 Survey Committee on State of the Profession

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., Comment: 29 pages; 5 figures; re-submitted to MNRAS after replying to referee's report

The Stratospheric Observatory for Infrared Astronomy (SOFIA) is an airborne observatory, carrying a 2.5 m telescope onboard a heavily modified Boeing 747SP aircraft. SOFIA is optimized for operation at infrared wavelengths, much of which is obscured for ground-based observatories by atmospheric water vapor. The SOFIA science instrument complement consists of seven instruments: FORCAST (Faint Object InfraRed CAmera for the SOFIA Telescope), GREAT (German Receiver for Astronomy at Terahertz Frequencies), HIPO (High-speed Imaging Photometer for Occultations), FLITECAM (First Light Infrared Test Experiment CAMera), FIFI-LS (Far-Infrared Field-Imaging Line Spectrometer), EXES (Echelon-Cross-Echelle Spectrograph), and HAWC (High-resolution Airborne Wideband Camera). FORCAST is a 540 m imager with grism spectroscopy, developed at Cornell University. GREAT is a heterodyne spectrometer providing high-resolution spectroscopy in several bands from 60240 m, developed at the Max Planck Institute for Radio Astronomy. HIPO is a 0.31.1 m imager, developed at Lowell Observatory. FLITECAM is a 15 m wide-field imager with grism spectroscopy, developed at UCLA. FIFI-LS is a 42210 m integral field imaging grating spectrometer, developed at the University of Stuttgart. EXES is a 528 m high-resolution spectrograph, developed at UC Davis and NASA ARC. HAWC is a 50240 m imager, developed at the University of Chicago, and undergoing an upgrade at JPL to add polarimetry capability and substantially larger GSFC detectors. We describe the capabilities, performance, and status of each instrument, highlighting science results obtained using FORCAST, GREAT, and HIPO during SOFIA Early Science observations conducted in 2011.

LFS was an educational outreach adventure which brought the excitement of astronomical exploration on NASA's Kuiper Airborne Observatory (KAO) to a nationwide audience of children, parents and children through live, interactive television, broadcast from the KAO at an altitude of 41,000 feet during an actual scientific observing mission. The project encompassed three KAO flights during the fall of 1995, including a short practice mission, a daytime observing flight between Moffett Field, California to Houston, Texas, and a nighttime mission from Houston back to Moffett Field. The University of Chicago infrared research team participated in planning the program, developing auxiliary materials including background information and lesson plans, developing software which allowed students on the ground to control the telescope and on-board cameras via the Internet from the Adler Planetarium in Chicago, and acting as on-camera correspondents to explain and answer questions about the scientific research conducted during the flights.

The scientific objectives, engineering analysis and design, results of technology development, and focal-plane instrumentation for a two-meter balloon-borne telescope for far-infrared and submillimeter astronomy are presented. The unique capabilities of balloon-borne observations are discussed. A program summary emphasizes the development of the two-meter design. The relationship of the Large Deployable Reflector (LDR) is also discussed. Detailed treatment is given to scientific objectives, gondola design, the mirror development program, experiment accommodations, ground support equipment requirements, NSBF design drivers and payload support requirements, the implementation phase summary development plan, and a comparison of three-meter and two-meter gondola concepts.