Search

We present our sixth set of results from our mid-infrared imaging survey of Milky Way Giant HII regions with our detailed analysis of NGC 3603, the most luminous GHII region in the Galaxy. We used imaging data from the FORCAST instrument on the Stratospheric Observatory For Infrared Astronomy (SOFIA) at 20 and 37 microns which mapped the central ~8.5'x8.5' infrared-emitting area of NGC 3603 at a spatial resolution of <~3". Utilizing these SOFIA data in conjunction with multi-wavelength observations from the near-infrared to radio, including Spitzer-IRAC and Herschel-PACS archival data, we investigate the physical nature of individual infrared sources and sub-components within NGC 3603. For individual compact sources we used the multi-wavelength photometry data to construct spectral energy distributions (SEDs) and fit them with massive young stellar object (MYSO) SED models, and find 14 sources that are likely to be MYSOs. We also detect dust emission from the 3 massive proplyd candidates, as well as from the disk and outflow of the evolved blue supergiant, Sher 25. Utilizing multi-wavelength data, we derived luminosity-to-mass ratio and virial parameters for the star-forming clumps within NGC 3603, estimating their relative ages and finding that NGC 3603 is an older GHII region overall, compared to our previously studied GHII regions. We discuss how NGC 3603, which we categorize as a 'cavity-type' GHII region, exhibits a more modest number of MYSOs and molecular clumps when compared to the 'distributed-type' GHII regions that share similar Lyman continuum photon rates., Comment: 32 pages, 15 figures, accepted for publication in ApJ

We present our fifth set of results from our mid-infrared imaging survey of Milky Way Giant HII (GHII) regions with our detailed analysis of DR7 and K3-50. We obtained 20/25 and 37um imaging maps of both regions using the FORCAST instrument on the Stratospheric Observatory For Infrared Astronomy (SOFIA). We investigate the multi-scale properties of DR7 and K3-50 using our data in conjunction with previous multi-wavelength observations. Near to far-infrared spectral energy distributions of individual compact infrared sources were constructed and fitted with massive young stellar object (MYSO) models. We find eight out of the ten (80%) compact sources in K3-50 and three out of the four (75%) sources in DR7 are likely to be MYSOs. We derived luminosity-to-mass ratios of the extended radio sub-regions of DR7 and K3-50 to estimate their relative ages. The large spread in evolutionary state for the sub-regions in K3-50 likely indicates that the star-forming complex has undergone multiple star-forming events separated more widely in time, whereas the smaller spread in DR7 likely indicates the star formation sub-regions are more co-eval. DR7 and K3-50 have Lyman continuum photon rates just above the formal threshold criterion for being categorized as a GHII region (10^50 photons/s) but with large enough errors that this classification is uncertain. By measuring other observational characteristics in the infrared, we find that K3-50 has properties more akin to previous bona fide GHII regions we have studied, whereas DR7 has values more like those of the non-GHII regions we have previously studied., Comment: 29 pages, 10 figures, accepted for publication in ApJ

At a distance of 50 kpc, Supernova 1987A is an ideal target to study how a young supernova (SN) evolves in time. Its equatorial ring, filled with material expelled from the progenitor star about 20,000 years ago, has been engulfed with SN blast waves. Shocks heat dust grains in the ring, emitting their energy at mid-infrared (IR) wavelengths We present ground-based 10--18$\mu$m monitoring of the ring of SN 1987A from day 6067 to 12814 at a resolution of 0.5", together with SOFIA photometry at 10-30 $\mu$m. The IR images in the 2000's (day 6067-7242) showed that the shocks first began brightening the east side of the ring. Later, our mid-IR images from 2017 to 2022 (day 10952-12714) show that dust emission is now fading in the east, while it has brightened on the west side of the ring. Because dust grains are heated in the shocked plasma, which can emit X-rays, the IR and X-ray brightness ratio represent shock diagnostics. Until 2007 the IR to X-ray brightness ratio remained constant over time, and during this time shocks seemed to be largely influencing the east side of the ring. However, since then, the IR to X-ray ratio has been declining, due to increased X-ray brightness. Whether the declining IR brightness is because of dust grains being destroyed or being cooled in the post-shock regions will require more detailed modelling., Comment: MNRAS accepted

We present a detailed study of the massive star-forming region G35.2-0.74N with ALMA 1.3 mm multi-configuration observations. At 0.2" (440 au) resolution, the continuum emission reveals several dense cores along a filamentary structure, consistent with previous ALMA 0.85 mm observations. At 0.03" (66 au) resolution, we detect 22 compact sources, most of which are associated with the filament. Four of the sources are associated with compact centimeter continuum emission, and two of these are associated with H30{\alpha} recombination line emission. The H30{\alpha} line kinematics show ordered motion of the ionized gas, consistent with disk rotation and/or outflow expansion. We construct models of photoionized regions to simultaneously fit the multi-wavelength free-free fluxes and the H30{\alpha} total fluxes. The derived properties suggest the presence of at least three massive young stars with nascent hypercompact Hii regions. Two of these ionized regions are surrounded by a large rotating structure that feeds two individual disks, revealed by dense gas tracers, such as SO2, H2CO, and CH3OH. In particular, the SO2 emission highlights two spiral structures in one of the disks and probes the faster-rotating inner disks. The 12CO emission from the general region reveals a complex outflow structure, with at least four outflows identified. The remaining 18 compact sources are expected to be associated with lower-mass protostars forming in the vicinity of the massive stars. We find potential evidence for disk disruption due to dynamical interactions in the inner region of this protocluster. The spatial distribution of the sources suggests a smooth overall radial density gradient without subclustering, but with tentative evidence of primordial mass segregation., Comment: 46 pages, 21 figures, 5 tables. Accepted to ApJ

We present $\sim10-40\,\mu$m SOFIA-FORCAST images of 11 isolated protostars as part of the SOFIA Massive (SOMA) Star Formation Survey, with this morphological classification based on 37 $\mu$m imaging. We develop an automated method to define source aperture size using the gradient of its background-subtracted enclosed flux and apply this to build spectral energy distributions (SEDs). We fit the SEDs with radiative transfer models, developed within the framework of turbulent core accretion (TCA) theory, to estimate key protostellar properties. Here, we release the sedcreator python package that carries out these methods. The SEDs are generally well fitted by the TCA models, from which we infer initial core masses $M_c$ ranging from $20-430\:M_\odot$, clump mass surface densities $\Sigma_{\rm cl}\sim0.3-1.7\:{\rm{g\:cm}}^{-2}$ and current protostellar masses $m_*\sim3-50\:M_\odot$. From a uniform analysis of the 40 sources in the full SOMA survey to date, we find that massive protostars form across a wide range of clump mass surface density environments, placing constraints on theories that predict a minimum threshold $\Sigma_{\rm cl}$ for massive star formation. However, the upper end of the $m_*-\Sigma_{\rm cl}$ distribution follows trends predicted by models of internal protostellar feedback that find greater star formation efficiency in higher $\Sigma_{\rm cl}$ conditions. We also investigate protostellar far-IR variability by comparison with IRAS data, finding no significant variation over an $\sim$40 year baseline., Comment: Accepted for publication in ApJ. Published version including the proofs

This is the fourth paper exploring the infrared properties of giant HII regions with the FORCAST instrument on the Stratospheric Observatory For Infrared Astronomy (SOFIA). Our survey utilizes the census of 56 Milky Way giant HII regions identified by Conti & Crowther (2004), and in this paper we present the 20 and 37 micron imaging data we have obtained from SOFIA for sources Sgr D and W42. Based upon the SOFIA data and other multi-wavelength data, we derive and discuss the detailed physical properties of the individual compact sources and sub-regions as well as the large scale properties of Sgr D and W42. However, improved measurements have revealed much closer distances to both regions than previously believed, and consequently both sources are not powerful enough to be considered giant HII regions any longer. Motivated by this, we revisit the census of giant HII regions, performing a search through the last two decades of literature to update each source with the most recent and/or most accurate distance measurements. Based on these new distance estimates, we determine that 14 sources in total (25%) are at sufficiently reliable and closer distances that they are not powerful enough to be considered giant HII regions. We briefly discuss the observational and physical characteristics specific to Sgr D and W42 and show that they have properties distinct from the giant HII regions previously studied as a part of this survey., Comment: 31 pages, 8 figures, 7 tables; accepted for publication in ApJ

We present our third set of results from our mid-infrared imaging survey of Milky Way Giant HII (GHII) regions with our detailed analysis of W49A, one of the most distant, yet most luminous, GHII regions in the Galaxy. We used the FORCAST instrument on the Stratospheric Observatory For Infrared Astronomy (SOFIA) to obtain 20 and 37$\mu$m images of the entire ~5.0' x 3.5' infrared-emitting area of W49A at a spatial resolution of ~3". Utilizing these SOFIA data in conjunction with previous multi-wavelength observations from the near-infrared to radio, including Spitzer-IRAC and Herschel-PACS archival data, we investigate the physical nature of individual infrared sources and sub-components within W49A. For individual compact sources we used the multi-wavelength photometry data to construct spectral energy distributions (SEDs) and fit them with massive young stellar object (MYSO) SED models, and find 22 sources that are likely to be MYSOs. Ten new sources are identified for the first time in this work. Even at 37$\mu$m we are unable to detect infrared emission from the sources on the western side of the extremely extinguished ring of compact radio emission sources known as the Welch Ring. Utilizing multi-wavelength data, we derived luminosity-to-mass ratio and virial parameters of the extended radio sub-regions of W49A to estimate their relative ages and find that overall the sub-components of W49A have a very small spread in evolutionary state compared to our previously studied GHII regions., Comment: 35 pages, 15 figures, 10 tables, accepted for publication in ApJ

We present our third set of results from our mid-infrared imaging survey of Milky Way Giant HII (GHII) regions with our detailed analysis of W49A, one of the most distant, yet most luminous, GHII regions in the Galaxy. We used the FORCAST instrument on the Stratospheric Observatory For Infrared Astronomy (SOFIA) to obtain 20 and 37$\mu$m images of the entire ~5.0' x 3.5' infrared-emitting area of W49A at a spatial resolution of ~3". Utilizing these SOFIA data in conjunction with previous multi-wavelength observations from the near-infrared to radio, including Spitzer-IRAC and Herschel-PACS archival data, we investigate the physical nature of individual infrared sources and sub-components within W49A. For individual compact sources we used the multi-wavelength photometry data to construct spectral energy distributions (SEDs) and fit them with massive young stellar object (MYSO) SED models, and find 22 sources that are likely to be MYSOs. Ten new sources are identified for the first time in this work. Even at 37$\mu$m we are unable to detect infrared emission from the sources on the western side of the extremely extinguished ring of compact radio emission sources known as the Welch Ring. Utilizing multi-wavelength data, we derived luminosity-to-mass ratio and virial parameters of the extended radio sub-regions of W49A to estimate their relative ages and find that overall the sub-components of W49A have a very small spread in evolutionary state compared to our previously studied GHII regions., Comment: 35 pages, 15 figures, 10 tables, accepted for publication in ApJ

We present $\sim10-40\,\mu$m SOFIA-FORCAST images of 14 intermediate-mass protostar candidates as part of the SOFIA Massive (SOMA) Star Formation Survey. We build spectral energy distributions (SEDs), also utilizing archival Spitzer, Herschel and IRAS data. We then fit the SEDs with radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, to estimate key protostellar properties. With the addition of these intermediate-mass sources, SOMA protostars span luminosities from $\sim10^{2}-10^{6}\:L_{\odot}$, current protostellar masses from $\sim0.5-30\:M_{\odot}$ and ambient clump mass surface densities, $\Sigma_{\rm cl}$ from $0.1-3\:{\rm{g\:cm}^{-2}}$. A wide range of evolutionary states of the individual protostars and of the protocluster environments are also probed. We have also considered about 50 protostars identified in Infrared Dark Clouds and expected to be at the earliest stages of their evolution. With this global sample, most of the evolutionary stages of high- and intermediate-mass protostars are probed. From the best fitting models, there is no evidence of a threshold value of protocluster clump mass surface density being needed to form protostars up to $\sim25\:M_\odot$. However, to form more massive protostars, there is tentative evidence that $\Sigma_{\rm{cl}}$ needs to be $\gtrsim1\:{\rm{g\,cm}}^{-2}$. We discuss how this is consistent with expectations from core accretion models that include internal feedback from the forming massive star., Comment: 40 pages, 21 figures, 4 tables, accepted to ApJ

We present $\sim10-40\,\mu$m SOFIA-FORCAST images of 14 intermediate-mass protostar candidates as part of the SOFIA Massive (SOMA) Star Formation Survey. We build spectral energy distributions (SEDs), also utilizing archival Spitzer, Herschel and IRAS data. We then fit the SEDs with radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, to estimate key protostellar properties. With the addition of these intermediate-mass sources, SOMA protostars span luminosities from $\sim10^{2}-10^{6}\:L_{\odot}$, current protostellar masses from $\sim0.5-30\:M_{\odot}$ and ambient clump mass surface densities, $\Sigma_{\rm cl}$ from $0.1-3\:{\rm{g\:cm}^{-2}}$. A wide range of evolutionary states of the individual protostars and of the protocluster environments are also probed. We have also considered about 50 protostars identified in Infrared Dark Clouds and expected to be at the earliest stages of their evolution. With this global sample, most of the evolutionary stages of high- and intermediate-mass protostars are probed. From the best fitting models, there is no evidence of a threshold value of protocluster clump mass surface density being needed to form protostars up to $\sim25\:M_\odot$. However, to form more massive protostars, there is tentative evidence that $\Sigma_{\rm{cl}}$ needs to be $\gtrsim1\:{\rm{g\,cm}}^{-2}$. We discuss how this is consistent with expectations from core accretion models that include internal feedback from the forming massive star., Comment: 40 pages, 21 figures, 4 tables, accepted to ApJ

We present our second set of results from our mid-infrared imaging survey of Milky Way Giant HII regions. We used the FORCAST instrument on the Stratospheric Observatory For Infrared Astronomy to obtain 20 and 37$\mu$m images of the central ~10'X10' area of M17. We investigate the small- and large-scale properties of M17 using our data in conjunction with previous multi-wavelength observations. The spectral energy distributions of individual compact sources were constructed with Spitzer-IRAC, SOFIA-FORCAST, and Herschel-PACS photometry data and fitted with massive young stellar object (MYSO) models. Seven sources were found to match the criteria for being MYSO candidates, four of which are identified here for the first time, and the stellar mass of the most massive object, UC1, is determined to be 64 solar mass. We resolve the extended mid-infrared emission from the KW Object, and suggest that the angle of this extended emission is influenced by outflow. IRS5 is shown to decrease in brightness as a function of wavelength from the mid- to far-infrared, and has several other indicators that point to it being an intermediate mass Class II object and not a MYSO. We find that the large-scale appearance of emission in M17 at 20$\mu$m is significantly affected by contamination from the [SIII] emission line from the ionized gas of the Giant HII region. Finally, a number of potential evolutionary tracers yield a consistent picture suggesting that the southern bar of M17 is likely younger than the northern bar., Comment: 25 pages, 16 figures, 9 tables, accepted for publication in ApJ

We present multi-wavelength images observed with SOFIA-FORCAST from $\sim$10 to 40 $\mu$m of seven high luminosity massive protostars, as part of the SOFIA Massive (SOMA) Star Formation Survey. Source morphologies at these wavelengths appear to be influenced by outflow cavities and extinction from dense gas surrounding the protostars. Using these images, we build spectral energy distributions (SEDs) of the protostars, also including archival data from Spitzer, Herschel and other facilities. Radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, are then fit to the SEDs to estimate key properties of the protostars. Considering the best five models fit to each source, the protostars have masses $m_{*} \sim 12-64 \: M_{\odot}$ accreting at rates of $\dot{m}_{*} \sim 10^{-4}-10^{-3} \: M_{\odot} \: \rm yr^{-1}$ inside cores of initial masses $M_{c} \sim 100-500 \: M_{\odot}$ embedded in clumps with mass surface densities $\Sigma_{\rm cl} \sim 0.1-3 \: \rm g \: cm^{-2}$ and span a luminosity range of $10^{4} -10^{6} \: L_{\odot}$. Compared with the first eight protostars in Paper I, the sources analyzed here are more luminous, and thus likely to be more massive protostars. They are often in a clustered environment or have a companion protostar relatively nearby. From the range of parameter space of the models, we do not see any evidence that $\Sigma_{\rm cl}$ needs to be high to form these massive stars. For most sources the RT models provide reasonable fits to the SEDs, though the cold clump material often influences the long wavelength fitting. However, for sources in very clustered environments, the model SEDs may not be such a good description of the data, indicating potential limitations of the models for these regions., Comment: 30 pages, 19 figures, Accepted for publication in ApJ

Almost all massive stars have bound stellar companions, existing in binaries or higher-order multiples. While binarity is theorized to be an essential feature of how massive stars form, essentially all information about such properties is derived from observations of already formed stars, whose orbital properties may have evolved since birth. Little is known about binarity during formation stages. Here we report high angular resolution observations of 1.3 mm continuum and H30alpha recombination line emission, which reveal a massive protobinary with apparent separation of 180 au at the center of the massive star-forming region IRAS07299-1651. From the line-of-sight velocity difference of 9.5 km/s of the two protostars, the binary is estimated to have a minimum total mass of 18 solar masses, consistent with several other metrics, and maximum period of 570 years, assuming a circular orbit. The H30alpha line from the primary protostar shows kinematics consistent with rotation along a ring of radius of 12 au. The observations indicate that disk fragmentation at several hundred au may have formed the binary, and much smaller disks are feeding the individual protostars., Comment: Published in Nature Astronomy. This is author's version. Full article is available here (https://rdcu.be/brENk). 47 pages, 10 figures, including methods and supplementary information

We discuss the first results from our mid-infrared imaging survey of Milky Way Giant HII regions with our detailed analysis of W51A, which is one of the largest GHII regions in our Galaxy. We used the FORCAST instrument on SOFIA to obtain 20 and 37$\mu$m images of the central $10' \times 20'$ area, which encompasses both of the G49.5-0.4 and G49.4-0.3 sub-regions. Based on these new data, and in conjunction with previous multi-wavelength observations, we conjecture on the physical nature of several individual sources and sub-components within W51A. We find that extinction seems to play an important role in the observed structures we see in the near- to mid-infrared, both globally and locally. We used the SOFIA photometry combined with Spitzer-IRAC and Herschel-PACS photometry data to construct spectral energy distributions (SEDs) of sub-components and point sources detected in the SOFIA images. We fit those SEDs with young stellar object models, and found 41 sources that are likely to be massive young stellar objects, many of which are identified as such in this work for the first time. Close to half of the massive young stellar objects do not have detectable radio continuum emission at cm wavelengths, implying a very young state of formation. We derived luminosity-to-mass ratio and virial parameters of the extended radio sub-regions of W51A to estimate their relative ages., Comment: 35 pages, 15 figures, 6 tables, Accepted by ApJ

Almost all massive stars have bound stellar companions, existing in binaries or higher-order multiples. While binarity is theorized to be an essential feature of how massive stars form, essentially all information about such properties is derived from observations of already formed stars, whose orbital properties may have evolved since birth. Little is known about binarity during formation stages. Here we report high angular resolution observations of 1.3 mm continuum and H30alpha recombination line emission, which reveal a massive protobinary with apparent separation of 180 au at the center of the massive star-forming region IRAS07299-1651. From the line-of-sight velocity difference of 9.5 km/s of the two protostars, the binary is estimated to have a minimum total mass of 18 solar masses, consistent with several other metrics, and maximum period of 570 years, assuming a circular orbit. The H30alpha line from the primary protostar shows kinematics consistent with rotation along a ring of radius of 12 au. The observations indicate that disk fragmentation at several hundred au may have formed the binary, and much smaller disks are feeding the individual protostars., Comment: Published in Nature Astronomy. This is author's version. Full article is available here (https://rdcu.be/brENk). 47 pages, 10 figures, including methods and supplementary information

We present multi-wavelength images observed with SOFIA-FORCAST from $\sim$10 to 40 $\mu$m of seven high luminosity massive protostars, as part of the SOFIA Massive (SOMA) Star Formation Survey. Source morphologies at these wavelengths appear to be influenced by outflow cavities and extinction from dense gas surrounding the protostars. Using these images, we build spectral energy distributions (SEDs) of the protostars, also including archival data from Spitzer, Herschel and other facilities. Radiative transfer (RT) models of Zhang & Tan (2018), based on Turbulent Core Accretion theory, are then fit to the SEDs to estimate key properties of the protostars. Considering the best five models fit to each source, the protostars have masses $m_{*} \sim 12-64 \: M_{\odot}$ accreting at rates of $\dot{m}_{*} \sim 10^{-4}-10^{-3} \: M_{\odot} \: \rm yr^{-1}$ inside cores of initial masses $M_{c} \sim 100-500 \: M_{\odot}$ embedded in clumps with mass surface densities $\Sigma_{\rm cl} \sim 0.1-3 \: \rm g \: cm^{-2}$ and span a luminosity range of $10^{4} -10^{6} \: L_{\odot}$. Compared with the first eight protostars in Paper I, the sources analyzed here are more luminous, and thus likely to be more massive protostars. They are often in a clustered environment or have a companion protostar relatively nearby. From the range of parameter space of the models, we do not see any evidence that $\Sigma_{\rm cl}$ needs to be high to form these massive stars. For most sources the RT models provide reasonable fits to the SEDs, though the cold clump material often influences the long wavelength fitting. However, for sources in very clustered environments, the model SEDs may not be such a good description of the data, indicating potential limitations of the models for these regions., Comment: 30 pages, 19 figures, Accepted for publication in ApJ

We discuss the first results from our mid-infrared imaging survey of Milky Way Giant HII regions with our detailed analysis of W51A, which is one of the largest GHII regions in our Galaxy. We used the FORCAST instrument on SOFIA to obtain 20 and 37$\mu$m images of the central $10' \times 20'$ area, which encompasses both of the G49.5-0.4 and G49.4-0.3 sub-regions. Based on these new data, and in conjunction with previous multi-wavelength observations, we conjecture on the physical nature of several individual sources and sub-components within W51A. We find that extinction seems to play an important role in the observed structures we see in the near- to mid-infrared, both globally and locally. We used the SOFIA photometry combined with Spitzer-IRAC and Herschel-PACS photometry data to construct spectral energy distributions (SEDs) of sub-components and point sources detected in the SOFIA images. We fit those SEDs with young stellar object models, and found 41 sources that are likely to be massive young stellar objects, many of which are identified as such in this work for the first time. Close to half of the massive young stellar objects do not have detectable radio continuum emission at cm wavelengths, implying a very young state of formation. We derived luminosity-to-mass ratio and virial parameters of the extended radio sub-regions of W51A to estimate their relative ages., Comment: 35 pages, 15 figures, 6 tables, Accepted by ApJ

We present our second set of results from our mid-infrared imaging survey of Milky Way Giant HII regions. We used the FORCAST instrument on the Stratospheric Observatory For Infrared Astronomy to obtain 20 and 37$\mu$m images of the central ~10'X10' area of M17. We investigate the small- and large-scale properties of M17 using our data in conjunction with previous multi-wavelength observations. The spectral energy distributions of individual compact sources were constructed with Spitzer-IRAC, SOFIA-FORCAST, and Herschel-PACS photometry data and fitted with massive young stellar object (MYSO) models. Seven sources were found to match the criteria for being MYSO candidates, four of which are identified here for the first time, and the stellar mass of the most massive object, UC1, is determined to be 64 solar mass. We resolve the extended mid-infrared emission from the KW Object, and suggest that the angle of this extended emission is influenced by outflow. IRS5 is shown to decrease in brightness as a function of wavelength from the mid- to far-infrared, and has several other indicators that point to it being an intermediate mass Class II object and not a MYSO. We find that the large-scale appearance of emission in M17 at 20$\mu$m is significantly affected by contamination from the [SIII] emission line from the ionized gas of the Giant HII region. Finally, a number of potential evolutionary tracers yield a consistent picture suggesting that the southern bar of M17 is likely younger than the northern bar., Comment: 25 pages, 16 figures, 9 tables, accepted for publication in ApJ

We report molecular line observations of the massive protostellar source G339.88-1.26 with the Atacama Large Millimeter/Submillimeter Array. The observations reveal a highly collimated SiO jet extending from the 1.3 mm continuum source, which connects to a slightly wider but still highly collimated CO outflow. Rotational features perpendicular to the outflow axis are detected in various molecular emissions, including SiO, SO2, H2S, CH3OH, and H2CO emissions. Based on their spatial distributions and kinematics, we find that they trace different parts of the envelope-disk system. The SiO emission traces the disk and inner envelope in addition to the jet. The CH3OH and H2CO emissions mostly trace the infalling-rotating envelope, and are enhanced around the transition region between envelope and disk, i.e., the centrifugal barrier. The SO2 and H2S emissions are enhanced around the centrifugal barrier, and also trace the outer part of the disk. Envelope kinematics are consistent with rotating-infalling motion, while those of the disk are consistent with Keplerian rotation. The radius and velocity of the centrifugal barrier are estimated to be about 530 au and 6 km/s, leading to a central mass of about 11 solar masses, consistent with estimates based on spectral energy distribution fitting. These results indicate that an ordered transition from an infalling-rotating envelope to a Keplerian disk through a centrifugal barrier, accompanied by changes of types of molecular line emissions, is a valid description of this massive protostellar source. This implies that at least some massive stars form in a similar way as low-mass stars via Core Accretion., Comment: Accepted to ApJ

The equatorial ring of Supernova (SN) 1987A has been exposed to forward shocks from the SN blast wave, and it has been suggested that these forward shocks have been causing on-going destruction of dust in the ring. We obtained SOFIA FORCAST 11.1, 19.7 and 31.5 micron photometry of SN\,1987A in 2016. Compared with Spitzer measurements 10 years earlier, the 31.5 micron flux has significantly increased. The excess at 31.5 micron appears to be related to the Herschel 70 micron excess, which was detected 5 years earlier. The dust mass needed to account for the the 31.5--70 micron excess is 3--7x10^-4 Msun, more than ten times larger than the ring dust mass (1x10^-5 Msun) estimate from the data 10-years earlier. We argue that dust grains are re-formed or grown in the post-shock regions in the ring after forward shocks have destroyed pre-existing dust grains in the ring and released refractory elements into gas. In the post-shock region, atoms can stick to surviving dust grains, and the dust mass may have increased (grain growth), or dust grains might have condensed directly from the gas. An alternative possibility is that the outer part of the expanding ejecta dust might have been heated by X-ray emission from the circumstellar ring. The future development of this excess could reveal whether grains are reformed in the post-shocked region of the ring or eject dust is heated by X-ray., Comment: 10 mages; MNRAS accented

We report molecular line observations of the massive protostellar source G339.88-1.26 with the Atacama Large Millimeter/Submillimeter Array. The observations reveal a highly collimated SiO jet extending from the 1.3 mm continuum source, which connects to a slightly wider but still highly collimated CO outflow. Rotational features perpendicular to the outflow axis are detected in various molecular emissions, including SiO, SO2, H2S, CH3OH, and H2CO emissions. Based on their spatial distributions and kinematics, we find that they trace different parts of the envelope-disk system. The SiO emission traces the disk and inner envelope in addition to the jet. The CH3OH and H2CO emissions mostly trace the infalling-rotating envelope, and are enhanced around the transition region between envelope and disk, i.e., the centrifugal barrier. The SO2 and H2S emissions are enhanced around the centrifugal barrier, and also trace the outer part of the disk. Envelope kinematics are consistent with rotating-infalling motion, while those of the disk are consistent with Keplerian rotation. The radius and velocity of the centrifugal barrier are estimated to be about 530 au and 6 km/s, leading to a central mass of about 11 solar masses, consistent with estimates based on spectral energy distribution fitting. These results indicate that an ordered transition from an infalling-rotating envelope to a Keplerian disk through a centrifugal barrier, accompanied by changes of types of molecular line emissions, is a valid description of this massive protostellar source. This implies that at least some massive stars form in a similar way as low-mass stars via Core Accretion., Comment: Accepted to ApJ

The equatorial ring of Supernova (SN) 1987A has been exposed to forward shocks from the SN blast wave, and it has been suggested that these forward shocks have been causing on-going destruction of dust in the ring. We obtained SOFIA FORCAST 11.1, 19.7 and 31.5 micron photometry of SN\,1987A in 2016. Compared with Spitzer measurements 10 years earlier, the 31.5 micron flux has significantly increased. The excess at 31.5 micron appears to be related to the Herschel 70 micron excess, which was detected 5 years earlier. The dust mass needed to account for the the 31.5--70 micron excess is 3--7x10^-4 Msun, more than ten times larger than the ring dust mass (1x10^-5 Msun) estimate from the data 10-years earlier. We argue that dust grains are re-formed or grown in the post-shock regions in the ring after forward shocks have destroyed pre-existing dust grains in the ring and released refractory elements into gas. In the post-shock region, atoms can stick to surviving dust grains, and the dust mass may have increased (grain growth), or dust grains might have condensed directly from the gas. An alternative possibility is that the outer part of the expanding ejecta dust might have been heated by X-ray emission from the circumstellar ring. The future development of this excess could reveal whether grains are reformed in the post-shocked region of the ring or eject dust is heated by X-ray., Comment: 10 mages; MNRAS accented

Debris disk morphology is wavelength dependent due to the wide range of particle sizes and size-dependent dynamics influenced by various forces. Resolved images of nearby debris disks reveal complex disk structures that are difficult to distinguish from their spectral energy distributions. Therefore, multi-wavelength resolved images of nearby debris systems provide an essential foundation to understand the intricate interplay between collisional, gravitational, and radiative forces that govern debris disk structures. We present the SOFIA 35 um resolved disk image of epsilon Eri, the closest debris disk around a star similar to the early Sun. Combining with the Spitzer resolved image at 24 um and 15-38 um excess spectrum, we examine two proposed origins of the inner debris in epsilon Eri: (1) in-situ planetesimal belt(s) and (2) dragged-in grains from the cold outer belt. We find that the presence of in-situ dust-producing planetesmial belt(s) is the most likely source of the excess emission in the inner 25 au region. Although a small amount of dragged-in grains from the cold belt could contribute to the excess emission in the inner region, the resolution of the SOFIA data is high enough to rule out the possibility that the entire inner warm excess results from dragged-in grains, but not enough to distinguish one broad inner disk from two narrow belts., Comment: AJ in press

Debris disk morphology is wavelength dependent due to the wide range of particle sizes and size-dependent dynamics influenced by various forces. Resolved images of nearby debris disks reveal complex disk structures that are difficult to distinguish from their spectral energy distributions. Therefore, multi-wavelength resolved images of nearby debris systems provide an essential foundation to understand the intricate interplay between collisional, gravitational, and radiative forces that govern debris disk structures. We present the SOFIA 35 um resolved disk image of epsilon Eri, the closest debris disk around a star similar to the early Sun. Combining with the Spitzer resolved image at 24 um and 15-38 um excess spectrum, we examine two proposed origins of the inner debris in epsilon Eri: (1) in-situ planetesimal belt(s) and (2) dragged-in grains from the cold outer belt. We find that the presence of in-situ dust-producing planetesmial belt(s) is the most likely source of the excess emission in the inner 25 au region. Although a small amount of dragged-in grains from the cold belt could contribute to the excess emission in the inner region, the resolution of the SOFIA data is high enough to rule out the possibility that the entire inner warm excess results from dragged-in grains, but not enough to distinguish one broad inner disk from two narrow belts., Comment: AJ in press

We present an overview and first results of the Stratospheric Observatory For Infrared Astronomy Massive (SOMA) Star Formation Survey, which is using the FORCAST instrument to image massive protostars from $\sim10$--$40\:\rm{\mu}\rm{m}$. These wavelengths trace thermal emission from warm dust, which in Core Accretion models mainly emerges from the inner regions of protostellar outflow cavities. Dust in dense core envelopes also imprints characteristic extinction patterns at these wavelengths, causing intensity peaks to shift along the outflow axis and profiles to become more symmetric at longer wavelengths. We present observational results for the first eight protostars in the survey, i.e., multiwavelength images, including some ancillary ground-based MIR observations and archival {\it{Spitzer}} and {\it{Herschel}} data. These images generally show extended MIR/FIR emission along directions consistent with those of known outflows and with shorter wavelength peak flux positions displaced from the protostar along the blueshifted, near-facing sides, thus confirming qualitative predictions of Core Accretion models. We then compile spectral energy distributions and use these to derive protostellar properties by fitting theoretical radiative transfer models. Zhang and Tan models, based on the Turbulent Core Model of McKee and Tan, imply the sources have protostellar masses $m_*\sim10$--50$\:M_\odot$ accreting at $\sim10^{-4}$--$10^{-3}\:M_\odot\:{\rm{yr}}^{-1}$ inside cores of initial masses $M_c\sim30$--500$\:M_\odot$ embedded in clumps with mass surface densities $\Sigma_{\rm{cl}}\sim0.1$--3$\:{\rm{g\:cm}^{-2}}$. Fitting Robitaille et al. models typically leads to slightly higher protostellar masses, but with disk accretion rates $\sim100\times$ smaller. We discuss reasons for these differences and overall implications of these first survey results for massive star formation theories., Comment: Accepted to ApJ, 32 pages

We present an overview and first results of the Stratospheric Observatory For Infrared Astronomy Massive (SOMA) Star Formation Survey, which is using the FORCAST instrument to image massive protostars from $\sim10$--$40\:\rm{\mu}\rm{m}$. These wavelengths trace thermal emission from warm dust, which in Core Accretion models mainly emerges from the inner regions of protostellar outflow cavities. Dust in dense core envelopes also imprints characteristic extinction patterns at these wavelengths, causing intensity peaks to shift along the outflow axis and profiles to become more symmetric at longer wavelengths. We present observational results for the first eight protostars in the survey, i.e., multiwavelength images, including some ancillary ground-based MIR observations and archival {\it{Spitzer}} and {\it{Herschel}} data. These images generally show extended MIR/FIR emission along directions consistent with those of known outflows and with shorter wavelength peak flux positions displaced from the protostar along the blueshifted, near-facing sides, thus confirming qualitative predictions of Core Accretion models. We then compile spectral energy distributions and use these to derive protostellar properties by fitting theoretical radiative transfer models. Zhang and Tan models, based on the Turbulent Core Model of McKee and Tan, imply the sources have protostellar masses $m_*\sim10$--50$\:M_\odot$ accreting at $\sim10^{-4}$--$10^{-3}\:M_\odot\:{\rm{yr}}^{-1}$ inside cores of initial masses $M_c\sim30$--500$\:M_\odot$ embedded in clumps with mass surface densities $\Sigma_{\rm{cl}}\sim0.1$--3$\:{\rm{g\:cm}^{-2}}$. Fitting Robitaille et al. models typically leads to slightly higher protostellar masses, but with disk accretion rates $\sim100\times$ smaller. We discuss reasons for these differences and overall implications of these first survey results for massive star formation theories., Comment: Accepted to ApJ, 32 pages

We present 30 and 40 micron imaging of the massive protostar G35.20-0.74 with SOFIA-FORCAST. The high surface density of the natal core around the protostar leads to high extinction, even at these relatively long wavelengths, causing the observed flux to be dominated by that emerging from the near-facing outflow cavity. However, emission from the far-facing cavity is still clearly detected. We combine these results with fluxes from the near-infrared to mm to construct a spectral energy distribution (SED). For isotropic emission the bolometric luminosity would be 3.3x10^4 Lsun. We perform radiative transfer modeling of a protostar forming by ordered, symmetric collapse from a massive core bounded by a clump with high mass surface density, Sigma_cl. To fit the SED requires protostellar masses ~20-34 Msun depending on the outflow cavity opening angle (35 - 50 degrees), and Sigma_cl ~ 0.4-1 g cm-2. After accounting for the foreground extinction and the flashlight effect, the true bolometric luminosity is ~ (0.7-2.2)x10^5 Lsun. One of these models also has excellent agreement with the observed intensity profiles along the outflow axis at 10, 18, 31 and 37 microns. Overall our results support a model of massive star formation involving the relatively ordered, symmetric collapse of a massive, dense core and the launching bipolar outflows that clear low density cavities. Thus a unified model may apply for the formation of both low and high mass stars., Comment: 6 pages, 4 figures, 1 table, accepted to ApJ

The Stratospheric Observatory For Infrared Astronomy, or SOFIA, is the largest flying observatory ever built,consisting of a 2.7-meter diameter telescope embedded in a modified Boeing 747-SP aircraft. SOFIA is a joint project between NASA and the German Aerospace Center Deutsches Zentrum fur Luft und-Raumfahrt (DLR). By flying at altitudes up to 45000 feet, the observatory gets above 99.9 percent of the infrared-absorbing water vapor in the Earth's atmosphere. This opens up an almost uninterrupted wavelength range from 0.3-1600 microns that is in large part obscured from ground based observatories. Since its 'Initial Science Flight' in December 2010, SOFIA has flown several dozen science flights, and has observed a wide array of objects from Solar System bodies, to stellar nurseries, to distant galaxies. This paper reviews a few of the exciting new science results from these first flights which were made by three instruments: the mid-infrared camera FORCAST, the far-infrared heterodyne spectrometer GREAT, and the optical occultation photometer HIPO., Comment: 6 pages; From the Proceedings of the Advanced Maui Optical and Space Surveillance Technologies (AMOS) Conference 2012

We present 30 and 40 micron imaging of the massive protostar G35.20-0.74 with SOFIA-FORCAST. The high surface density of the natal core around the protostar leads to high extinction, even at these relatively long wavelengths, causing the observed flux to be dominated by that emerging from the near-facing outflow cavity. However, emission from the far-facing cavity is still clearly detected. We combine these results with fluxes from the near-infrared to mm to construct a spectral energy distribution (SED). For isotropic emission the bolometric luminosity would be 3.3x10^4 Lsun. We perform radiative transfer modeling of a protostar forming by ordered, symmetric collapse from a massive core bounded by a clump with high mass surface density, Sigma_cl. To fit the SED requires protostellar masses ~20-34 Msun depending on the outflow cavity opening angle (35 - 50 degrees), and Sigma_cl ~ 0.4-1 g cm-2. After accounting for the foreground extinction and the flashlight effect, the true bolometric luminosity is ~ (0.7-2.2)x10^5 Lsun. One of these models also has excellent agreement with the observed intensity profiles along the outflow axis at 10, 18, 31 and 37 microns. Overall our results support a model of massive star formation involving the relatively ordered, symmetric collapse of a massive, dense core and the launching bipolar outflows that clear low density cavities. Thus a unified model may apply for the formation of both low and high mass stars., Comment: 6 pages, 4 figures, 1 table, accepted to ApJ

The Stratospheric Observatory For Infrared Astronomy, or SOFIA, is the largest flying observatory ever built,consisting of a 2.7-meter diameter telescope embedded in a modified Boeing 747-SP aircraft. SOFIA is a joint project between NASA and the German Aerospace Center Deutsches Zentrum fur Luft und-Raumfahrt (DLR). By flying at altitudes up to 45000 feet, the observatory gets above 99.9 percent of the infrared-absorbing water vapor in the Earth's atmosphere. This opens up an almost uninterrupted wavelength range from 0.3-1600 microns that is in large part obscured from ground based observatories. Since its 'Initial Science Flight' in December 2010, SOFIA has flown several dozen science flights, and has observed a wide array of objects from Solar System bodies, to stellar nurseries, to distant galaxies. This paper reviews a few of the exciting new science results from these first flights which were made by three instruments: the mid-infrared camera FORCAST, the far-infrared heterodyne spectrometer GREAT, and the optical occultation photometer HIPO., Comment: 6 pages; From the Proceedings of the Advanced Maui Optical and Space Surveillance Technologies (AMOS) Conference 2012

We present infrared observations of the ultra-compact H II region W3(OH) made by the FORCAST instrument aboard SOFIA and by Spitzer/IRAC. We contribute new wavelength data to the spectral energy distribution, which constrains the optical depth, grain size distribution, and temperature gradient of the dusty shell surrounding the H II region. We model the dust component as a spherical shell containing an inner cavity with radius ~ 600 AU, irradiated by a central star of type O9 and temperature ~ 31,000 K. The total luminosity of this system is 71,000 L_solar. An observed excess of 2.2 - 4.5 microns emission in the SED can be explained by our viewing a cavity opening or clumpiness in the shell structure whereby radiation from the warm interior of the shell can escape. We claim to detect the nearby water maser source W3 (H2O) at 31.4 and 37.1 microns using beam deconvolution of the FORCAST images. We constrain the flux densities of this object at 19.7 - 37.1 microns. Additionally, we present in situ observations of four young stellar and protostellar objects in the SOFIA field, presumably associated with the W3 molecular cloud. Results from the model SED fitting tool of Robitaille et al. (2006, 2007} suggest that two objects (2MASS J02270352+6152357 and 2MASS J02270824+6152281) are intermediate-luminosity (~ 236 - 432 L_solar) protostars; one object (2MASS J02270887+6152344) is either a high-mass protostar with luminosity 3000 L_solar or a less massive young star with a substantial circumstellar disk but depleted envelope; and one object (2MASS J02270743+6152281) is an intermediate-luminosity (~ 768 L_solar) protostar nearing the end of its envelope accretion phase or a young star surrounded by a circumstellar disk with no appreciable circumstellar envelope., Comment: 12 pages, 8 figures, 2 tables, accepted by ApJ

Milliarcsecond VLBI maps of regions containing 6.7 GHz methanol maser emission have lead to the recent discovery of ring-like distributions of maser spots and the plausible hypothesis that they may be tracing circumstellar disks around forming high mass stars. We aimed to test this hypothesis by imaging these regions in the near and mid-infrared at high spatial resolution and compare the observed emission to the expected infrared morphologies as inferred from the geometries of the maser rings. In the near infrared we used the Gemini North adaptive optics system of Altair/NIRI, while in the mid-infrared we used the combination of the Gemini South instrument T-ReCS and super-resolution techniques. Resultant images had a resolution of approximately 150 mas in both the near-infrared and mid-infrared. We discuss the expected distribution of circumstellar material around young and massive accreting (proto)stars and what infrared emission geometries would be expected for the different maser ring orientations under the assumption that the masers are coming from within circumstellar disks. Based upon the observed infrared emission geometries for the four targets in our sample and the results of SED modeling of the massive young stellar objects associated with the maser rings, we do not find compelling evidence in support of the hypothesis that methanol masers rings reside in circumstellar disks., Comment: 20 pages, 7 figures, accepted for publication by ApJ; article with full-resolution figures can be downloaded at http://www.jim-debuizer.net

We present new mid-infrared images of the central region of the Orion Nebula using the newly commissioned SOFIA airborne telescope and its 5 -- 40 micron camera FORCAST. The 37.1 micron images represent the highest resolution observations (<4") ever obtained of this region at these wavelengths. After BN/KL (which is described in a separate letter in this issue), the dominant source at all wavelengths except 37.1 micron is the Ney-Allen Nebula, a crescent-shaped extended source associated with theta 1D. The morphology of the Ney-Allen nebula in our images is consistent with the interpretation that it is ambient dust swept up by the stellar wind from theta 1D, as suggested by Smith et al. (2005). Our observations also reveal emission from two "proplyds" (proto-planetary disks), and a few embedded young stellar objects (YSOs; IRc9, and OMC1S IRS1, 2, and 10). The spectral energy distribution for IRc9 is presented and fitted with standard YSO models from Robitaille et al. (2007) to constrain the total luminosity, disk size, and envelope size. The diffuse, nebular emission we observe at all FORCAST wavelengths is most likely from the background photodissociation region (PDR) and shows structure that coincides roughly with H_alpha and [N II] emission. We conclude that the spatial variations in the diffuse emission are likely due to undulations in the surface of the background PDR., Comment: Accepted for publication in the ApJ Letters special SOFIA issue. 15 pages; 5 figures (color)

The BN/KL region of the Orion Nebula is the nearest region of high mass star formation in our galaxy. As such, it has been the subject of intense investigation at a variety of wavelengths, which have revealed it to be brightest in the infrared to sub-mm wavelength regime. Using the newly commissioned SOFIA airborne telescope and its 5-40 micron camera FORCAST, images of the entire BN/KL complex have been acquired. The 31.5 and 37.1 micron images represent the highest resolution observations (<=4") ever obtained of this region at these wavelengths. These observations reveal that the BN object is not the dominant brightness source in the complex at wavelengths >31.5 microns, and that this distinction goes instead to the source IRc4. It was determined from these images and derived dust color temperature maps that IRc4 is also likely to be self-luminous. A new source of emission has also been identified at wavelengths >31.5 microns that coincides with the northeastern outflow lobe from the protostellar disk associated with radio source I., Comment: 6 pages, 3 figures, accepted for publication in ApJL

We present 37\micron\ imaging of the S140 complex of infrared sources centered on IRS1 made with the FORCAST camera on SOFIA. These observations are the longest wavelength imaging to resolve clearly the three main sources seen at shorter wavelengths, IRS 1, 2 and 3, and are nearly at the diffraction limit of the 2.5-m telescope. We also obtained a small number of images at 11 and 31\micron\ that are useful for flux measurement. Our images cover the area of several strong sub-mm sources seen in the area -- SMM 1, 2, and 3 -- that are not coincident with any mid-infrared sources and are not visible in our longer wavelength imaging either. Our new observations confirm previous estimates of the relative dust optical depth and source luminosity for the components in this likely cluster of early B stars. We also investigate the use of super-resolution to go beyond the basic diffraction limit in imaging on SOFIA and find that the van Cittert algorithm, together with the "multi-resolution" technique, provides excellent results.

We examine eight young stellar objects in the OMC-2 star forming region based on observations from the SOFIA/FORCAST early science phase, the Spitzer Space Telescope, the Herschel Space Observatory, 2MASS, APEX, and other results in the literature. We show the spectral energy distributions of these objects from near-infrared to millimeter wavelengths, and compare the SEDs with those of sheet collapse models of protostars and circumstellar disks. Four of the objects can be modelled as protostars with infalling envelopes, two as young stars surrounded by disks, and the remaining two objects have double-peaked SEDs. We model the double-peaked sources as binaries containing a young star with a disk and a protostar. The six most luminous sources are found in a dense group within a 0.15 x 0.25 pc region; these sources have luminosities ranging from 300 L_sun to 20 L_sun. The most embedded source (OMC-2 FIR 4) can be fit by a class 0 protostar model having a luminosity of ~50 L_sun and mass infall rate of ~10^-4 solar masses per year., Comment: Accepted by ApJ Letters

We present infrared observations of the ultra-compact H II region W3(OH) made by the FORCAST instrument aboard SOFIA and by Spitzer/IRAC. We contribute new wavelength data to the spectral energy distribution, which constrains the optical depth, grain size distribution, and temperature gradient of the dusty shell surrounding the H II region. We model the dust component as a spherical shell containing an inner cavity with radius ~ 600 AU, irradiated by a central star of type O9 and temperature ~ 31,000 K. The total luminosity of this system is 71,000 L_solar. An observed excess of 2.2 - 4.5 microns emission in the SED can be explained by our viewing a cavity opening or clumpiness in the shell structure whereby radiation from the warm interior of the shell can escape. We claim to detect the nearby water maser source W3 (H2O) at 31.4 and 37.1 microns using beam deconvolution of the FORCAST images. We constrain the flux densities of this object at 19.7 - 37.1 microns. Additionally, we present in situ observations of four young stellar and protostellar objects in the SOFIA field, presumably associated with the W3 molecular cloud. Results from the model SED fitting tool of Robitaille et al. (2006, 2007} suggest that two objects (2MASS J02270352+6152357 and 2MASS J02270824+6152281) are intermediate-luminosity (~ 236 - 432 L_solar) protostars; one object (2MASS J02270887+6152344) is either a high-mass protostar with luminosity 3000 L_solar or a less massive young star with a substantial circumstellar disk but depleted envelope; and one object (2MASS J02270743+6152281) is an intermediate-luminosity (~ 768 L_solar) protostar nearing the end of its envelope accretion phase or a young star surrounded by a circumstellar disk with no appreciable circumstellar envelope., Comment: 12 pages, 8 figures, 2 tables, accepted by ApJ

Milliarcsecond VLBI maps of regions containing 6.7 GHz methanol maser emission have lead to the recent discovery of ring-like distributions of maser spots and the plausible hypothesis that they may be tracing circumstellar disks around forming high mass stars. We aimed to test this hypothesis by imaging these regions in the near and mid-infrared at high spatial resolution and compare the observed emission to the expected infrared morphologies as inferred from the geometries of the maser rings. In the near infrared we used the Gemini North adaptive optics system of Altair/NIRI, while in the mid-infrared we used the combination of the Gemini South instrument T-ReCS and super-resolution techniques. Resultant images had a resolution of approximately 150 mas in both the near-infrared and mid-infrared. We discuss the expected distribution of circumstellar material around young and massive accreting (proto)stars and what infrared emission geometries would be expected for the different maser ring orientations under the assumption that the masers are coming from within circumstellar disks. Based upon the observed infrared emission geometries for the four targets in our sample and the results of SED modeling of the massive young stellar objects associated with the maser rings, we do not find compelling evidence in support of the hypothesis that methanol masers rings reside in circumstellar disks., Comment: 20 pages, 7 figures, accepted for publication by ApJ; article with full-resolution figures can be downloaded at http://www.jim-debuizer.net

We present new mid-infrared images of the central region of the Orion Nebula using the newly commissioned SOFIA airborne telescope and its 5 -- 40 micron camera FORCAST. The 37.1 micron images represent the highest resolution observations (<4") ever obtained of this region at these wavelengths. After BN/KL (which is described in a separate letter in this issue), the dominant source at all wavelengths except 37.1 micron is the Ney-Allen Nebula, a crescent-shaped extended source associated with theta 1D. The morphology of the Ney-Allen nebula in our images is consistent with the interpretation that it is ambient dust swept up by the stellar wind from theta 1D, as suggested by Smith et al. (2005). Our observations also reveal emission from two "proplyds" (proto-planetary disks), and a few embedded young stellar objects (YSOs; IRc9, and OMC1S IRS1, 2, and 10). The spectral energy distribution for IRc9 is presented and fitted with standard YSO models from Robitaille et al. (2007) to constrain the total luminosity, disk size, and envelope size. The diffuse, nebular emission we observe at all FORCAST wavelengths is most likely from the background photodissociation region (PDR) and shows structure that coincides roughly with H_alpha and [N II] emission. We conclude that the spatial variations in the diffuse emission are likely due to undulations in the surface of the background PDR., Comment: Accepted for publication in the ApJ Letters special SOFIA issue. 15 pages; 5 figures (color)

We present 37\micron\ imaging of the S140 complex of infrared sources centered on IRS1 made with the FORCAST camera on SOFIA. These observations are the longest wavelength imaging to resolve clearly the three main sources seen at shorter wavelengths, IRS 1, 2 and 3, and are nearly at the diffraction limit of the 2.5-m telescope. We also obtained a small number of images at 11 and 31\micron\ that are useful for flux measurement. Our images cover the area of several strong sub-mm sources seen in the area -- SMM 1, 2, and 3 -- that are not coincident with any mid-infrared sources and are not visible in our longer wavelength imaging either. Our new observations confirm previous estimates of the relative dust optical depth and source luminosity for the components in this likely cluster of early B stars. We also investigate the use of super-resolution to go beyond the basic diffraction limit in imaging on SOFIA and find that the van Cittert algorithm, together with the "multi-resolution" technique, provides excellent results.

We examine eight young stellar objects in the OMC-2 star forming region based on observations from the SOFIA/FORCAST early science phase, the Spitzer Space Telescope, the Herschel Space Observatory, 2MASS, APEX, and other results in the literature. We show the spectral energy distributions of these objects from near-infrared to millimeter wavelengths, and compare the SEDs with those of sheet collapse models of protostars and circumstellar disks. Four of the objects can be modelled as protostars with infalling envelopes, two as young stars surrounded by disks, and the remaining two objects have double-peaked SEDs. We model the double-peaked sources as binaries containing a young star with a disk and a protostar. The six most luminous sources are found in a dense group within a 0.15 x 0.25 pc region; these sources have luminosities ranging from 300 L_sun to 20 L_sun. The most embedded source (OMC-2 FIR 4) can be fit by a class 0 protostar model having a luminosity of ~50 L_sun and mass infall rate of ~10^-4 solar masses per year., Comment: Accepted by ApJ Letters

The BN/KL region of the Orion Nebula is the nearest region of high mass star formation in our galaxy. As such, it has been the subject of intense investigation at a variety of wavelengths, which have revealed it to be brightest in the infrared to sub-mm wavelength regime. Using the newly commissioned SOFIA airborne telescope and its 5-40 micron camera FORCAST, images of the entire BN/KL complex have been acquired. The 31.5 and 37.1 micron images represent the highest resolution observations (<=4") ever obtained of this region at these wavelengths. These observations reveal that the BN object is not the dominant brightness source in the complex at wavelengths >31.5 microns, and that this distinction goes instead to the source IRc4. It was determined from these images and derived dust color temperature maps that IRc4 is also likely to be self-luminous. A new source of emission has also been identified at wavelengths >31.5 microns that coincides with the northeastern outflow lobe from the protostellar disk associated with radio source I., Comment: 6 pages, 3 figures, accepted for publication in ApJL

"Green fuzzies" or "extended green objects" were discovered in the recent Spitzer GLIMPSE survey data. These extended sources have enhanced emission in the 4.5um IRAC channel images (which are generally assigned to be green when making 3-color RGB images from Spitzer data). Green fuzzies are frequently found in the vicinities of massive young stellar objects, and it has been established that they are in some cases associated with outflows. Nevertheless, the spectral carrier(s) of this enhanced emission is still uncertain. Although it has been suggested that Br Alpha, H2, [Fe II], and/or broad CO emission may be contributing to and enhancing the 4.5um flux from these objects, to date there have been no direct observations of the 4-5um spectra of these objects. We report here on the first direct spectroscopic identification of the origin of the green fuzzy emission. We obtained spatially resolved L and M band spectra for two green fuzzy sources using NIRI on the Gemini North telescope. In the case of one source, G19.88-0.53, we detect three individual knots of green fuzzy emission around the source. The knots exhibit a pure molecular hydrogen line emission spectrum, with the 4.695um v=0-0 S(9) line dominating the emission in the 4-5um wavelength range, and no detected continuum component. Our data for G19.88-0.53 prove that green fuzzy emission can be due primarily to emission lines of molecular hydrogen within the bandpass of the IRAC 4.5um channel. However, the other target observed, G49.27-0.34, does not exhibit any line emission and appears to be an embedded massive young stellar object with a cometary UC HII region. We suggest that the effects of extinction in the 3-8um wavelength range and an exaggeration in the color stretch of the 4.5um filter in IRAC RGB images could lead to embedded sources such as this one falsely appearing "green"., Comment: 10 pages, 5 figures; Accepted for publication by AJ

We have used the Spitzer satellite to monitor the mid-IR evolution of SN 1987A over a 5 year period spanning the epochs between days 6000 and 8000 since the explosion. The supernova (SN) has evolved into a supernova remnant (SNR) and its radiative output is dominated by the interaction of the SN blast wave with the pre-existing equatorial ring (ER). The mid-IR spectrum is dominated by emission from ~180 K silicate dust, collisionally-heated by the hot X-ray emitting gas with a temperature and density of ~5x10^6 K and 3x10^4 cm-3, respectively. The mass of the radiating dust is ~1.2x10^(-6) Msun on day 7554, and scales linearly with IR flux. The infrared to soft-X-ray flux ratio is roughly constant with a value of 2.5. Gas-grain collisions therefore dominate the cooling of the shocked gas. The constancy of of this ratio suggests that very little grain processing or gas cooling have occurred throughout this epoch. The shape of the dust spectrum remained unchanged during the observations while the total flux increased with a time dependence of t^(0.87), t being the time since the first encounter between the blast wave and the ER. These observations are consistent with the transitioning of the blast wave from free expansion to a Sedov phase as it propagates into the main body of the ER., Comment: Accepted for publication in the ApJ, 11 pages, 11 figures

"Green fuzzies" or "extended green objects" were discovered in the recent Spitzer GLIMPSE survey data. These extended sources have enhanced emission in the 4.5um IRAC channel images (which are generally assigned to be green when making 3-color RGB images from Spitzer data). Green fuzzies are frequently found in the vicinities of massive young stellar objects, and it has been established that they are in some cases associated with outflows. Nevertheless, the spectral carrier(s) of this enhanced emission is still uncertain. Although it has been suggested that Br Alpha, H2, [Fe II], and/or broad CO emission may be contributing to and enhancing the 4.5um flux from these objects, to date there have been no direct observations of the 4-5um spectra of these objects. We report here on the first direct spectroscopic identification of the origin of the green fuzzy emission. We obtained spatially resolved L and M band spectra for two green fuzzy sources using NIRI on the Gemini North telescope. In the case of one source, G19.88-0.53, we detect three individual knots of green fuzzy emission around the source. The knots exhibit a pure molecular hydrogen line emission spectrum, with the 4.695um v=0-0 S(9) line dominating the emission in the 4-5um wavelength range, and no detected continuum component. Our data for G19.88-0.53 prove that green fuzzy emission can be due primarily to emission lines of molecular hydrogen within the bandpass of the IRAC 4.5um channel. However, the other target observed, G49.27-0.34, does not exhibit any line emission and appears to be an embedded massive young stellar object with a cometary UC HII region. We suggest that the effects of extinction in the 3-8um wavelength range and an exaggeration in the color stretch of the 4.5um filter in IRAC RGB images could lead to embedded sources such as this one falsely appearing "green"., Comment: 10 pages, 5 figures; Accepted for publication by AJ

We have used the Spitzer satellite to monitor the mid-IR evolution of SN 1987A over a 5 year period spanning the epochs between days 6000 and 8000 since the explosion. The supernova (SN) has evolved into a supernova remnant (SNR) and its radiative output is dominated by the interaction of the SN blast wave with the pre-existing equatorial ring (ER). The mid-IR spectrum is dominated by emission from ~180 K silicate dust, collisionally-heated by the hot X-ray emitting gas with a temperature and density of ~5x10^6 K and 3x10^4 cm-3, respectively. The mass of the radiating dust is ~1.2x10^(-6) Msun on day 7554, and scales linearly with IR flux. The infrared to soft-X-ray flux ratio is roughly constant with a value of 2.5. Gas-grain collisions therefore dominate the cooling of the shocked gas. The constancy of of this ratio suggests that very little grain processing or gas cooling have occurred throughout this epoch. The shape of the dust spectrum remained unchanged during the observations while the total flux increased with a time dependence of t^(0.87), t being the time since the first encounter between the blast wave and the ER. These observations are consistent with the transitioning of the blast wave from free expansion to a Sedov phase as it propagates into the main body of the ER., Comment: Accepted for publication in the ApJ, 11 pages, 11 figures

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the field of cosmology. Community-based White Paper to Astro2010 in response to the call for such papers., Comment: White Paper to the Cosmology and Fundamental Physics (GCT) Science Frontiers Panel of the Astro2010 Decadal Survey

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the fields of planetary systems and star formation. Community-based White Paper to Astro2010 in response to the call for such papers., Comment: White Paper to the Planetary Systems and Star Formation (PSF) Science Frontiers Panel of the Astro2010 Decadal Survey

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the field of Galactic science. Community-based White Paper to Astro2010 in response to the call for such papers., Comment: White Paper to the Galactic Neighborhood (GAN) Science Frontiers Panel of the Astro2010 Decadal Survey

Science opportunities and recommendations concerning optical/infrared polarimetry for the upcoming decade in the fields of stars and stellar evolution. Community-based White Paper to Astro2010 in response to the call for such papers., Comment: White Paper to the Stars and Stellar Evolution (SSE) Science Frontiers Panel of the Astro2010 Decadal Survey