Search

Polycyclic aromatic hydrocarbons are an important component of the interstellar medium of galaxies and photochemistry plays a key role in the evolution of these species in space. Here, we explore the photofragmentation behaviour of the coronene cation (C24H12+) using time of flight mass spectrometry. The experiments show photodissociation fragmentation channels including the formation of bare carbon clusters (Cn+) and hydrocarbon chains (CnHx+). The mass spectrum of coronene is dominated by peaks from C11+ and C7H+. Density functional theory was used to calculate relative energies, potential dissociation pathways, and possible structures for relevant species. We identify 6-6 to 5-7 ring isomerisation as a key step in the formation of both the bare carbon clusters and the hydrocarbon chains observed in this study. We present the dissociation mechanism outlined here as a potential formation route for C60 and other astrochemically relevant species., Comment: 17 pages, 12 figure. Version submitted after referee report to PCCP. Supplementary information available at https://www.rsc.org/suppdata/d4/cp/d4cp01301h/d4cp01301h1.pdf

The aromatic infrared Bands (AIBs) dominate the mid-infrared spectra of many galactic and extragalactic sources. These AIBs are generally attributed to fluorescent emission from aromatic molecules. Unified efforts from experimentalists and theoreticians to assign these AIB features have recently gotten additional impetus with the launch of the James Webb Space Telescope (JWST) as the Mid-InfraRed Instrument (MIRI) delivers mid-IR spectrum with greatly increased sensitivity and spectral resolution. PAHs in space can exist in either neutral or ionic forms, absorb UV photons and undergo fragmentation, becoming a rich source of small hydrocarbons. This top-down mechanism of larger PAHs fragmenting into smaller species is of utmost importance in photo-dissociation regions (PDR) in space. In this work, we experimentally and theoretically investigate the photo-fragmentation pathways of two astronomically significant PAH cations - corannulene (C20H10) and sumanene (C21H12) that are structural motifs of fullerene C60, to understand their sequential fragmentation pathways. The photo-fragmentation experiments exhibit channels that are much different from planar PAHs. The breakdown of carbon skeleton is found to have different pathways for C20H10 and C21H12 because of the number and positioning of pentagon rings; yet the most abundant low mass cations produced by these two species are found to be similar. The low mass cations showcased in this work could be of interest for their astronomical detections. For completeness, the qualitative photo fragmentation behaviour of the dicationic corannulene and sumanene have also been experimented, but the potential energy surface of these dications are beyond the scope of this paper., Comment: 17 pages, 12 figures. Version submitted after referee report to PCCP. Supplementary information available at https://www.rsc.org/suppdata/d4/cp/d4cp01247j/d4cp01247j1.pdf

The gas-phase abundance of carbon, x_C = C/H, and its depletion factors are essential parameters for understanding the gas and solid compositions that are ultimately incorporated into planets. The majority of protoplanetary disks are born in clusters and, as a result, are exposed to external FUV radiation. These FUV photons potentially affect the disk's evolution, chemical composition, and line excitation. We present the first detection of the [CI]609um fine-structure line of neutral carbon (CI), achieved with ALMA, toward one of these disks, d203-506, in the Orion Nebula Cluster. We also report the detection of CI forbidden and permitted lines (from electronically excited states up to 10 eV) observed with JWST in the IR. These lines trace the irradiated outer disk and photo-evaporative wind. Contrary to the common belief that these IR lines are C+ recombination lines, we find that they are dominated by FUV-pumping of CI followed by fluorescence cascades. They trace the transition from atomic to molecular gas, and their intensities scale with G0. The lack of outstanding IR OI fluorescent emission, however, implies a sharper attenuation of external FUV radiation with E > 12 eV (~Lyman-beta). This is related to a lower effective FUV dust absorption cross section compared to that of interstellar grains, implying a more prominent role for FUV shielding by the CI photoionization continuum. The [CI]609um intensity is proportional to N(CI) and can be used to infer x_C. We derive x_C ~ 1.4E-4. This implies that there is no major depletion of volatile carbon compared to x_C measured in the natal cloud, hinting at a young disk. We also show that external FUV radiation impacts the outer disk and wind by vertically shifting the water freeze-out depth, which results in less efficient grain growth and settling. This shift leads to nearly solar gas-phase C/O abundance ratios in these irradiated layers., Comment: Accepted for publication in A&A Letters. 14 pages including Appendices

Ro-vibrational absorption lines of H$_2$O in the 5-8 $\mu$m wavelength range selectively probe gas against the mid-infrared continuum emitting background of the inner regions of YSOs and AGN and deliver important information about these warm, dust-obscured environments. JWST/MIRI detects these lines in many lines of sight at a moderate spectral resolving power of $R\sim3500$ (FWHM of 85 km/s). Based on our analysis of high-resolution SOFIA/EXES observations, we find that the interpretation of JWST/MIRI absorption spectra can be severely hampered by the blending of individual transitions and the lost information on the intrinsic line width or the partial coverage of the background continuum source. In this paper, we point out problems such as degeneracy that arise in deriving physical properties from an insufficiently resolved spectrum. This can lead to differences in the column density by two orders of magnitude. We emphasize the importance of weighting optically thin and weak lines in spectral analyses and provide recipes for breaking down the coupled parameters. We also provide an online tool to generate the H$_2$O absorption line spectra that can be compared to observations., Comment: Accepted for publication in ApJS. 26 pages, 23 figures. Comments are more than welcome!

We present a study of hydrogen, helium, and carbon millimeter-wave radio-recombination lines (RRLs) toward ten representative positions throughout the Orion Nebula complex, using the Yebes 40m telescope in the Q band (31.3 GHz to 50.6 GHz) at an angular resolution of about $45\arcsec$ ($\sim$0.09\,pc). The observed positions include the Orion Nebula (M42) with the Orion Molecular Core 1, M43, and the Orion Molecular Core 3 bordering on NGC 1973, 1975, and 1977. While hydrogen and helium RRLs arise in the ionized gas surrounding the massive stars in the Orion Nebula complex, carbon RRLs stem from the neutral gas of the adjacent photo-dissociation regions (PDRs). The high velocity resolution ($0.3\,\mathrm{km\,s^{-1}}$) enables us to discern the detailed dynamics of the RRL emitting neutral and ionized gas. We compare the carbon RRLs with SOFIA/upGREAT observations of the [CII] $158\,\mu\mathrm{m}$ line and IRAM 30m observations of the $^{13}$CO (J=2-1) line. Using the [CII] and [$^{13}$CII] intensities with the carbon RRL intensities, we can infer physical conditions (electron temperature and electron density) in the PDR gas using non-LTE excitation models. Our observations are sensitive enough to detect faint lines toward two positions in OMC1, that may be attributed to RRLs of C$^+$ or O$^+$. In general, the RRL line widths of both the ionized and neutral gas, as well as the [CII] and $^{13}$CO line widths, are broader than thermal, indicating significant turbulence in the interstellar medium, which transitions from super-Alfv\'enic and subsonic in the ionized gas to sub-Alfv\'enic and supersonic in the molecular gas. At the scales probed by our observations, the turbulent pressure dominates the pressure balance in the neutral and molecular gas, while in the ionized gas the turbulent pressure is much smaller than the thermal pressure.

Mid-infrared emission features probe the properties of ionized gas, and hot or warm molecular gas. The Orion Bar is a frequently studied photodissociation region (PDR) containing large amounts of gas under these conditions, and was observed with the MIRI IFU aboard JWST as part of the "PDRs4All" program. The resulting IR spectroscopic images of high angular resolution (0.2") reveal a rich observational inventory of mid-IR emission lines, and spatially resolve the substructure of the PDR, with a mosaic cutting perpendicularly across the ionization front and three dissociation fronts. We extracted five spectra that represent the ionized, atomic, and molecular gas layers, and measured the most prominent gas emission lines. An initial analysis summarizes the physical conditions of the gas and the potential of these data. We identified around 100 lines, report an additional 18 lines that remain unidentified, and measured the line intensities and central wavelengths. The H I recombination lines originating from the ionized gas layer bordering the PDR, have intensity ratios that are well matched by emissivity coefficients from H recombination theory, but deviate up to 10% due contamination by He I lines. We report the observed emission lines of various ionization stages of Ne, P, S, Cl, Ar, Fe, and Ni, and show how certain line ratios vary between the five regions. We observe the pure-rotational H$_2$ lines in the vibrational ground state from 0-0 S(1) to 0-0 S(8), and in the first vibrationally excited state from 1-1 S(5) to 1-1 S(9). We derive H$_2$ excitation diagrams, and approximate the excitation with one thermal (~700 K) component representative of an average gas temperature, and one non-thermal component (~2700 K) probing the effect of UV pumping. We compare these results to an existing model for the Orion Bar PDR and highlight the differences with the observations., Comment: 26 pages, 12 figures, 3 tables. Submitted to A&A, under review (1st revision)

Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far-ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photo-dissociation regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, impacting planet formation within the disks. We report JWST and Atacama Large Millimetere Array observations of a FUV-irradiated protoplanetary disk in the Orion Nebula. Emission lines are detected from the PDR; modelling their kinematics and excitation allows us to constrain the physical conditions within the gas. We quantify the mass-loss rate induced by the FUV irradiation, finding it is sufficient to remove gas from the disk in less than a million years. This is rapid enough to affect giant planet formation in the disk.

Anharmonicity strongly influences the absorption and emission spectra of polycyclic aromatic hydrocarbon (PAH) molecules. Here, IR-UV ion-dip spectroscopy experiments together with detailed anharmonic computations reveal the presence of fundamental, overtone, as well as 2- and 3-quanta combination band transitions in the far- and mid-infrared absorption spectrum of phenylacetylene and its singly deuterated isotopologue. Strong absorption features in the 400-900 cm$^{\rm -1}$ range originate from CH(D) in-plane and out-of-plane wags and bends, as well as bending motions including the C$\equiv$C and CH bonds of the acetylene substituent and the aromatic ring. For phenylacetylene, every absorption feature is assigned either directly or indirectly to a single or multiple vibrational mode(s). The measured spectrum is dense, broad, and structureless in many regions but well characterized by computations. Upon deuteration, large isotopic shifts are observed. At frequencies above 1500 cm$^{\rm -1}$ for d$_1$-phenylacetylene, a one-to-one match is seen when comparing computations and experiment with all features assigned to combination bands and overtones. The C$\equiv$C stretch observed in phenylacetylene is not observed in d$_1$-phenylacetylene due to a computed 40-fold drop in intensity. Overall, a careful treatment of anharmonicity that includes 2- and 3-quanta modes is found to be crucial to understand the rich details of the infrared spectrum of phenylacetylene. Based on these results, it can be expected that such an all-inclusive anharmonic treatment will also be key for unraveling the infrared spectra of PAHs in general., Comment: Accepted version in the Journal of Chemical Physics

Water is a key ingredient for the emergence of life as we know it. Yet, its destruction and reformation in space remains unprobed in warm gas. Here, we detect the hydroxyl radical (OH) emission from a planet-forming disk exposed to external far-ultraviolet (FUV) radiation with the James Webb Space Telescope. The observations are confronted with the results of quantum dynamical calculations. The highly excited OH infrared rotational lines are the tell-tale signs of H2O destruction by FUV. The OH infrared ro-vibrational lines are attributed to chemical excitation via the key reaction O+H=OH+H which seeds the formation of water in the gas-phase. We infer that the equivalent of the Earth ocean's worth of water is destroyed per month and replenished. These results show that under warm and irradiated conditions water is destroyed and efficiently reformed via gas-phase reactions. This process, assisted by diffusive transport, could reduce the HDO/H2O ratio in the warm regions of planet-forming disks., Comment: Version submitted to Nature Astronomy

The signatures of polycyclic aromatic hydrocarbons (PAHs) have been observed in protoplanetary discs, and their emission features obtained from spectral energy distributions (SED) have been used in the literature to characterise their size and determine their abundance. Two simple disc models (uniform PAH distribution against a PAH gap in the inner disc) are compared to investigate the difference of their SED and obtainable information. We used the radiative transfer code RADMC-3D to model the SED of two protoplanetary discs orbiting a typical Herbig star, one of which features a depletion of PAHs in the inner disc. We further created artificial images of the discs at face-on view to extract radial profiles of the PAH emission in the infrared. We find that the extracted PAH features from an SED provide limited information about the PAHs in protoplanetary disc environments, except for the ionisation state. The distribution of PAHs in a protoplanetary disc influences the total observed PAH luminosity in a non-linear fashion and alters the relative strength between the 3.3\,$\mu$m and 11.3\,$\mu$m features. Furthermore, we produced radial profiles at the 3\,$\mu$m, 6\,$\mu$m and, 11\,$\mu$m PAH emission features and find that they follow a double power-law profile where the slope reflects the radiative environment (single photon regime vs. multi-photon regime) in which the PAHs lie. Using spatially resolved techniques such as IFU or imaging in the era of the James Webb Space Telescope, we find that multi-wavelength radial emission profiles will not only provide information on the spatial distribution of the PAHs, but may also provide information on their size and underlying UV environment, which is crucial for photo-evaporative disc wind models., Comment: 9 pages, 5 figures, accepted for Astronomy and Astrophysics

(Abridged) We investigate the impact of radiative feedback from massive stars on their natal cloud and focus on the transition from the HII region to the atomic PDR (crossing the ionisation front (IF)), and the subsequent transition to the molecular PDR (crossing the dissociation front (DF)). We use high-resolution near-IR integral field spectroscopic data from NIRSpec on JWST to observe the Orion Bar PDR as part of the PDRs4All JWST Early Release Science Program. The NIRSpec data reveal a forest of lines including, but not limited to, HeI, HI, and CI recombination lines, ionic lines, OI and NI fluorescence lines, Aromatic Infrared Bands (AIBs including aromatic CH, aliphatic CH, and their CD counterparts), CO2 ice, pure rotational and ro-vibrational lines from H2, and ro-vibrational lines HD, CO, and CH+, most of them detected for the first time towards a PDR. Their spatial distribution resolves the H and He ionisation structure in the Huygens region, gives insight into the geometry of the Bar, and confirms the large-scale stratification of PDRs. We observe numerous smaller scale structures whose typical size decreases with distance from Ori C and IR lines from CI, if solely arising from radiative recombination and cascade, reveal very high gas temperatures consistent with the hot irradiated surface of small-scale dense clumps deep inside the PDR. The H2 lines reveal multiple, prominent filaments which exhibit different characteristics. This leaves the impression of a "terraced" transition from the predominantly atomic surface region to the CO-rich molecular zone deeper in. This study showcases the discovery space created by JWST to further our understanding of the impact radiation from young stars has on their natal molecular cloud and proto-planetary disk, which touches on star- and planet formation as well as galaxy evolution., Comment: 52 pages, 30 figures, submitted to A&A

We explore the potential of the Gaussian Mixture Model (GMM), an unsupervised machine learning method, to identify coherent physical structures in the ISM. The implementation we present can be used on any kind of spatially and spectrally resolved data set. We provide a step-by-step guide to use these models on different sources and data sets. Following the guide, we run the models on NGC 1977, RCW 120 and RCW 49 using the [CII] 158 $\mu$m mapping observations from the SOFIA telescope. We find that the models identified 6, 4 and 5 velocity coherent physical structures in NGC 1977, RCW 120 and RCW 49, respectively, which are validated by analysing the observed spectra towards these structures and by comparison to earlier findings. In this work we demonstrate that GMM is a powerful tool that can better automate the process of spatial and spectral analysis to interpret mapping observations., Comment: 19 pages, 14 figures

It has long been discussed whether stellar feedback in the form of winds and/or radiation can shred the nascent molecular cloud, thereby controlling the star formation rate. However, directly probing and quantifying the impact of stellar feedback on the neutral gas of the nascent clouds is challenging. We present an investigation doing exactly that toward the RCW 79 HII region using the ionized carbon line at 158 $\mu$m ([CII]) from the FEEDBACK Legacy Survey. We combine this data with information on the dozen ionizing O stars responsible for the evolution of the region, and observe in [CII] for the first time both blue- and red-shifted mostly neutral high-velocity gas which reaches velocities up to 25 km s$^{-1}$ relative to the bulk emission of the molecular cloud. This high-velocity gas mostly contains neutral gas and partly forms a fragmented shell, similar to recently found shells in a few Galactic HII regions. However, this shell does not account for all of the observed neutral high-velocity gas. We also find high-velocity gas streaming out of the nascent cloud through holes and obtain a range of dynamical timescales below 1.0 Myr for the high-velocity gas which is well below the 2.3$\pm$0.5 Myr age of the OB cluster. This suggests a different scenario for the evolution of RCW 79, where the high-velocity gas is not solely stemming from a spherical expanding bubble, but also from gas recently ablated at the edge of the turbulent molecular cloud into the surrounding interstellar medium through low-pressure holes or chimneys. The resulting mass ejection rate estimate for the cloud is 0.9-3.5$\times$10$^{-2}$ M$_{\odot}$~yr$^{-1}$, which leads to short erosion timescales, i.e. $<$5 Myr, for the nascent molecular cloud. This finding provides direct observational evidence of rapid molecular cloud dispersal., Comment: 9 pages, 7 figures, accepted in A&A Letters

We investigate the physical structure and conditions of photodissociation regions (PDRs) and molecular gas within the Pillars of Creation in the Eagle Nebula using SOFIA FEEDBACK observations of the [C II] 158 micron line. These observations are velocity resolved to 0.5 km s$^{-1}$ and are analyzed alongside a collection of complimentary data with similar spatial and spectral resolution: the [O I] 63 micron line, also observed with SOFIA, and rotational lines of CO, HCN, HCO$^{+}$, CS, and N$_2$H$^{+}$. Using the superb spectral resolution of SOFIA, APEX, CARMA, and BIMA, we reveal the relationships between the warm PDR and cool molecular gas layers in context of the Pillars' kinematic structure. We assemble a geometric picture of the Pillars and their surroundings informed by illumination patterns and kinematic relationships and derive physical conditions in the PDRs associated with the Pillars. We estimate an average molecular gas density $n_{{\rm H}_2} \sim 1.3 \times 10^5$ cm$^{-3}$ and an average atomic gas density $n_{\rm H} \sim 1.8 \times 10^4$ cm$^{-3}$ and infer that the ionized, atomic, and molecular phases are in pressure equilibrium if the atomic gas is magnetically supported. We find pillar masses of 103, 78, 103, and 18 solar masses for P1a, P1b, P2, and P3 respectively, and evaporation times of $\sim$1-2 Myr. The dense clumps at the tops of the pillars are currently supported by the magnetic field. Our analysis suggests that ambipolar diffusion is rapid and these clumps are likely to collapse within their photoevaporation timescales., Comment: 42 pages, 16 figures. Accepted for publication in The Astronomical Journal

The JWST has captured the most detailed and sharpest infrared images ever taken of the inner region of the Orion Nebula, the nearest massive star formation region, and a prototypical highly irradiated dense photo-dissociation region (PDR). We investigate the fundamental interaction of far-ultraviolet photons with molecular clouds. The transitions across the ionization front (IF), dissociation front (DF), and the molecular cloud are studied at high-angular resolution. These transitions are relevant to understanding the effects of radiative feedback from massive stars and the dominant physical and chemical processes that lead to the IR emission that JWST will detect in many Galactic and extragalactic environments. Due to the proximity of the Orion Nebula and the unprecedented angular resolution of JWST, these data reveal that the molecular cloud borders are hyper structured at small angular scales of 0.1-1" (0.0002-0.002 pc or 40-400 au at 414 pc). A diverse set of features are observed such as ridges, waves, globules and photoevaporated protoplanetary disks. At the PDR atomic to molecular transition, several bright features are detected that are associated with the highly irradiated surroundings of the dense molecular condensations and embedded young star. Toward the Orion Bar PDR, a highly sculpted interface is detected with sharp edges and density increases near the IF and DF. This was predicted by previous modeling studies, but the fronts were unresolved in most tracers. A complex, structured, and folded DF surface was traced by the H2 lines. This dataset was used to revisit the commonly adopted 2D PDR structure of the Orion Bar. JWST provides us with a complete view of the PDR, all the way from the PDR edge to the substructured dense region, and this allowed us to determine, in detail, where the emission of the atomic and molecular lines, aromatic bands, and dust originate.

(Abridged) Mid-infrared observations of photodissociation regions (PDRs) are dominated by strong emission features called aromatic infrared bands (AIBs). The most prominent AIBs are found at 3.3, 6.2, 7.7, 8.6, and 11.2 $\mu$m. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST. We provide an inventory of the AIBs found in the Orion Bar, along with mid-IR template spectra from five distinct regions in the Bar: the molecular PDR, the atomic PDR, and the HII region. We use JWST NIRSpec IFU and MIRI MRS observations of the Orion Bar from the JWST Early Release Science Program, PDRs4All (ID: 1288). We extract five template spectra to represent the morphology and environment of the Orion Bar PDR. The superb sensitivity and the spectral and spatial resolution of these JWST observations reveal many details of the AIB emission and enable an improved characterization of their detailed profile shapes and sub-components. While the spectra are dominated by the well-known AIBs at 3.3, 6.2, 7.7, 8.6, 11.2, and 12.7 $\mu$m, a wealth of weaker features and sub-components are present. We report trends in the widths and relative strengths of AIBs across the five template spectra. These trends yield valuable insight into the photochemical evolution of PAHs, such as the evolution responsible for the shift of 11.2 $\mu$m AIB emission from class B$_{11.2}$ in the molecular PDR to class A$_{11.2}$ in the PDR surface layers. This photochemical evolution is driven by the increased importance of FUV processing in the PDR surface layers, resulting in a "weeding out" of the weakest links of the PAH family in these layers. For now, these JWST observations are consistent with a model in which the underlying PAH family is composed of a few species: the so-called 'grandPAHs'., Comment: 25 pages, 10 figures, to appear in A&A

We present new $^{13}$CO(1-0), C$^{18}$O(1-0), HCO$^{+}$(1-0) and H$^{13}$CO$^{+}$(1-0) maps from the IRAM 30m telescope, and a spectrally-resolved [CII] 158 $\mu$m map observed with the SOFIA telescope towards the massive DR21 cloud. This traces the kinematics from low- to high-density gas in the cloud which allows to constrain the formation scenario of the high-mass star forming DR21 ridge. The molecular line data reveals that the sub-filaments are systematically redshifted relative to the dense ridge. We demonstrate that [CII] unveils the surrounding CO-poor gas of the dense filaments in the DR21 cloud. We also show that this surrounding gas is organized in a flattened cloud with curved redshifted dynamics perpendicular to the ridge. The sub-filaments thus form in this curved and flattened mass reservoir. A virial analysis of the different lines indicates that self-gravity should drive the evolution of the ridge and surrounding cloud. Combining all results we propose that bending of the magnetic field, due to the interaction with a mostly atomic colliding cloud, explains the velocity field and resulting mass accretion on the ridge. This is remarkably similar to what was found for at least two nearby low-mass filaments. We tentatively propose that this scenario might be a widespread mechanism to initiate star formation in the Milky Way. However, in contrast to low-mass clouds, gravitational collapse plays a role on the pc scale of the DR21 ridge because of the higher density. This allows more effective mass collection at the centers of collapse and should facilitate massive cluster formation., Comment: 32 pages, 28 figures, accepted in ApJ

We present a charge distribution-based emission model that calculates the infrared spectrum of fullerenes (C$_{60}$). Analysis of the modelled spectrum of C$_{60}$ in various charge states shows that the relative intensity of the features in the 5-10 $\mu$m versus 15-20 $\mu$m can be used to probe the C$_{60}$ charge state in interstellar spectra. We further used our model to simulate emission from polycyclic aromatic hydrocarbons (PAHs) and C$_{60}$ at five positions in the cavity of reflection nebula NGC~7023. Specifically, we modelled the 6.2/11.2 band ratio for circumcoronene and circumcircumcoronene and the 7.0/19.0 band ratio for C$_{60}$ as a function of the ionization parameter $\gamma$. A comparison of the model results with the observed band ratios shows that the $\gamma$ values in the cavity do not vary significantly, suggesting that the emission in the cavity does not originate from locations at the projected distances. Furthermore, we find that the C$_{60}$ derived $\gamma$ values are lower than the PAH-derived values by an order of magnitude. We discuss likely scenarios for this discrepancy. In one scenario, we attribute the differences in the derived $\gamma$ values to the uncertainties in the electron recombination rates of PAHs and C$_{60}$. In the other scenario, we suggest that PAHs and C$_{60}$ are not co-spatial resulting in different $\gamma$ values from their respective models. We highlight that experiments to determine necessary rates will be required in validating either one of the scenarios., Comment: 10 pages, 11 figures, Accepted for publication in MNRAS

Polycyclic aromatic hydrocarbons (PAHs) have been detected in numerous circumstellar discs. We propose the continuous processing of PAHs through clustering, adsorption on dust grains, and their reverse-processes as key mechanisms to reduce the emission-capable PAH abundance in protoplanetary discs. This cycle of processing is driven by vertical turbulence in the disc mixing PAHs between the disc midplane and the photosphere. We used a theoretical Monte Carlo model for photodesorption and a coagulation code in the disc midplane to estimate the relevance and timescale of these processes in a Herbig Ae/Be disc environment. By combining these components in a 1D vertical model, we calculated the gas-phase depletion of PAHs that stick as clusters on dust grains. Our results show that the clustering of gas-phase PAHs is very efficient, and that clusters with more than 100 monomers can grow for years before they are able to freeze out in the disc midplane. Once a PAH cluster is frozen on the dust grain surface, the large heat capacity of these clusters prevents them from evaporating off the grains in UV-rich environments such as the photosphere. Therefore, the clustering of PAHs followed by freeze-out can lead to a depletion of gas-phase PAHs in discs. Evaluated over the lifetime of protoplanetary discs, we find a depletion of PAHs by a factor that ranges between 50 and 1000 compared to the standard ISM abundance of PAHs in the inner disc through turbulent processing. Through these processes, we favour PAHs smaller than circumovalene as the major gas-phase emitters of the disc photosphere as larger PAH monomers cannot photodesorb from the grain surface. These gas-phase PAHs co-exist with large PAH clusters sticking on dust grains. We find a close relation between the amount of PAHs frozen out on dust grains and the dust population, as well as the strength of the vertical turbulence., Comment: Accepted for Astronomy and Astrophysics. 18 pages, 14 figures

Molecular hydrogen clouds are a key component of the interstellar medium because they are the birthplaces for stars. They are embedded in atomic gas that pervades the interstellar space. However, the details of how molecular clouds assemble from and interact with the atomic gas are still largely unknown. As a result of new observations of the 158~$\mu$m line of ionized carbon CII in the Cygnus region within the FEEDBACK program on SOFIA (Stratospheric Observatory for Infrared Astronomy), we present compelling evidence that CII unveils dynamic interactions between cloud ensembles. This process is neither a head-on collision of fully molecular clouds nor a gentle merging ofonly atomic clouds. Moreover, we demonstrate that the dense molecular clouds associated with the DR21 and W75N star-forming regions and a cloud at higher velocity are embedded in atomic gas and all components interact over a large range of velocities (20 km/s). The atomic gas has a density of 100 cm$^{-3}$ and a temperature of 100 K. We conclude that the CII 158 $\mu$m line is an excellent tracer to witness the processes involved in cloud interactions and anticipate further detections of this phenomenon in other regions, Comment: Nature Astronomy in press

To understand star formation rates, studying feedback mechanisms that regulate star formation is necessary. The radiation emitted by nascent massive stars play a significant role in feedback by photo-dissociating and ionizing their parental molecular clouds. To gain a detailed picture of the physical processes, we mapped the photo-dissociation region (PDR) M17-SW in several fine structure and high-J CO lines with FIFI-LS, the far-infrared imaging spectrometer aboard SOFIA. An analysis of the CO and [O I]146$\mu$m line intensities, combined with the far infrared intensity, allows us to create a density and UV intensity map using a one dimensional model. The density map reveals a sudden change in the gas density crossing the PDR. The strengths and limits of the model and the locations of the ionization and photo-dissociation front of the edge-on PDR are discussed., Comment: 18 pages, 17 figures, accepted by ApJ

The properties of dust change during the transition from diffuse to dense clouds as a result of ice formation and dust coagulation, but much is still unclear about this transformation. We present 2-20 micron spectra of 49 field stars behind the Perseus and Serpens Molecular Clouds and establish relationships between the near-infrared continuum extinction (AK) and the depths of the 9.7 micron silicate (tau97) and 3.0 micron H2O ice (tau30) absorption bands. The tau97/AK ratio varies from large, diffuse interstellar medium-like values (~0.55), to much lower ratios (~0.26). Above extinctions of AK~1.2 (AV~10; Perseus, Lupus, dense cores) and ~2.0 (AV~17; Serpens), the tau97/AK ratio is lowest. The tau97/AK reduction from diffuse to dense clouds is consistent with a moderate degree of grain growth (sizes up to ~0.5 micron), increasing the near-infrared color excess (and thus AK), but not affecting ice and silicate band profiles. This grain growth process seems to be related to the ice column densities and dense core formation thresholds, highlighting the importance of density. After correction for Serpens foreground extinction, the H2O ice formation threshold is in the range of AK=0.31-0.40 (AV=2.6-3.4) for all clouds, and thus grain growth takes place after the ices are formed. Finally, abundant CH3OH ice (~21% relative to H2O) is reported for 2MASSJ18285266+0028242 (Serpens), a factor of >4 larger than for the other targets., Comment: 29 pages+17 pages appendix. Published in The Astrophysical Journal, 930, 2, 2022

We quantified the effects of stellar feedback in RCW 49 by determining the physical conditions in different regions using the [CII] 158 $\mu$m and [OI] 63 $\mu$m observations from SOFIA, the $^{12}$CO (3-2) observations from APEX and the H$_2$ line observations from Spitzer telescopes. Large maps of RCW 49 were observed with the SOFIA and APEX telescopes, while the Spitzer observations were only available towards three small areas. From our qualitative analysis, we found that the H$_2$ 0-0 S(2) emission line probes denser gas compared to the H$_2$ 0-0 S(1) line. In four regions ("northern cloud", "pillar", "ridge", and "shell"), we compared our observations with the updated PDR Toolbox models and derived the integrated far-ultraviolet flux between 6-13.6 eV ($G_{\rm 0}$), H nucleus density ($n$), temperatures and pressures. We found the ridge to have the highest $G_{\rm 0}$ (2.4 $\times$ 10$^3$ Habing units), while the northern cloud has the lowest $G_{\rm 0}$ (5 $\times$ 10$^2$ Habing units). This is a direct consequence of the location of these regions with respect to the Wd2 cluster. The ridge also has a high density (6.4 $\times$ 10$^3$ cm$^{-3}$), which is consistent with its ongoing star formation. Among the Spitzer positions, we found the one closest to the Wd2 cluster to be the densest, suggesting an early phase of star formation. Furthermore, the Spitzer position that overlaps with the shell was found to have the highest $G_{\rm 0}$ and we expect this to be a result of its proximity to an O9V star., Comment: 18 pages, 12 figures

We present [CII] 158 $\mu$m and [OI] 63 $\mu$m observations of the bipolar HII region RCW 36 in the Vela C molecular cloud, obtained within the SOFIA legacy project FEEDBACK, which is complemented with APEX $^{12/13}$CO(3-2) and Chandra X-ray (0.5-7 keV) data. This shows that the molecular ring, forming the waist of the bipolar nebula, expands with a velocity of 1 - 1.9 km s$^{-1}$. We also observe an increased linewidth in the ring indicating that turbulence is driven by energy injection from the stellar feedback. The bipolar cavity hosts blue-shifted expanding [CII] shells at 5.2$\pm$0.5$\pm$0.5 km s$^{-1}$ (statistical and systematic uncertainty) which indicates that expansion out of the dense gas happens non-uniformly and that the observed bipolar phase might be relatively short ($\sim$0.2 Myr). The X-ray observations show diffuse emission that traces a hot plasma, created by stellar winds, in and around RCW 36. At least 50 \% of the stellar wind energy is missing in RCW 36. This is likely due to leakage which is clearing even larger cavities around the bipolar RCW 36 region. Lastly, the cavities host high-velocity wings in [CII] which indicates relatively high mass ejection rates ($\sim$5$\times$10$^{-4}$ M$_{\odot}$ yr$^{-1}$). This could be driven by stellar winds and/or radiation pressure, but remains difficult to constrain. This local mass ejection, which can remove all mass within 1 pc of RCW 36 in 1-2 Myr, and the large-scale clearing of ambient gas in the Vela C cloud indicates that stellar feedback plays a significant role in suppressing the star formation efficiency (SFE)., Comment: 38 pages, 27 figures, 8 tables, accepted in ApJ

We present in this paper the results of high spectral resolution ($R$=88,100) spectroscopy at 4.7 $\mu$m with iSHELL/IRTF of hot molecular gas close to the massive binary protostar W3 IRS5. The binary was spatially resolved and the spectra of the two sources (MIR1 and MIR2) were obtained simultaneously for the first time. Hundreds of $^{12}$CO $\nu$=0-1, $\nu$=1-2 lines, and $\nu$=0-1 transitions of the isotopes of $^{12}$CO were detected in absorption, and are blue-shifted compared to the cloud velocity $v_{LSR}=-$38 km/s. We decompose and identify kinematic components from the velocity profiles, and apply rotation diagram and curve of growth analyses to determine their physical properties. Temperatures and column densities of the identified components range from 30$-$700 K and 10$^{21}-$10$^{22}$ cm$^{-2}$, respectively. Our curve of growth analyses consider two scenarios. One assumes a foreground slab with a partial covering factor, which well reproduces the absorption of most of the components. The other assumes a circumstellar disk with an outward decreasing temperature in the vertical direction, and reproduces the absorption of all the hot components. We attribute the physical origins of the identified components to the foreground envelope ($<$100 K), post-J-shock regions (200$-$300 K), and clumpy structures on the circumstellar disks ($\sim$600 K). We propose that the components with a J-shock origin are akin to water maser spots in the same region, and are complementing the physical information of water masers along the direction of their movements., Comment: Accepted for publication in ApJ. 26 pages, 12 figures, and 9 tables. Comments are more than welcome!

We present results of a wide-field (approximately 60 x 90 pc) ALMA mosaic of CO(2-1) and $^{13}$CO(2-1) emission from the molecular cloud associated with the 30 Doradus star-forming region. Three main emission complexes, including two forming a bowtie-shaped structure extending northeast and southwest from the central R136 cluster, are resolved into complex filamentary networks. Consistent with previous studies, we find that the central region of the cloud has higher line widths at fixed size relative to the rest of the molecular cloud and to other LMC clouds, indicating an enhanced level of turbulent motions. However, there is no clear trend in gravitational boundedness (as measured by the virial parameter) with distance from R136. Structures observed in $^{13}$CO are spatially coincident with filaments and are close to a state of virial equilibrium. In contrast, CO structures vary greatly in virialization, with low CO surface brightness structures outside of the main filamentary network being predominantly unbound. The low surface brightness structures constitute ~10% of the measured CO luminosity; they may be shredded remnants of previously star-forming gas clumps, or alternatively the CO-emitting parts of more massive, CO-dark structures., Comment: 24 pages; published in ApJ on 15 Jun 2022; associated data can be found at https://doi.org/10.13012/B2IDB-1671495_V1

Astronomical observations and analysis of stardust isolated from meteorites have revealed a highly diverse interstellar and circumstellar grain inventory, including a wide range of amorphous materials and crystalline compounds (silicates and carbon). This diversity reflects the wide range of stellar sources injecting solids into the interstellar medium each with its own physical characteristics such as density, temperature and elemental composition and highlights the importance of kinetics rather than thermodynamics in the formation of these compounds. Based upon the extensive literature on soot formation in terrestrial settings, detailed kinetic pathways have been identified for the formation of carbon dust in C-rich stellar ejecta. These have been incorporated in astronomical models for these environments. In recent years, the chemical routes in the nucleation of oxides and silicates have been the focus of much astronomical research. These aspects of stardust formation will be reviewed and lessons for dust formation in planetary atmospheres will be drawn with the emphasis on the influence of kinetics on the characteristics and structure of dust in these environments.

Ionized gas probes the influence of massive stars on their environment. The Cygnus X region (d~1.5 kpc) is one of the most massive star forming complexes in our Galaxy, in which the Cyg OB2 association (age of 3-5 Myr and stellar mass $2 \times 10^{4}$ M$_{\odot}$) has a dominant influence. We observe the Cygnus X region at 148 MHz using the Low Frequency Array (LOFAR) and take into account short-spacing information during image deconvolution. Together with data from the Canadian Galactic Plane Survey, we investigate the morphology, distribution, and physical conditions of low-density ionized gas in a $4^{\circ} \times 4^{\circ}$ (100 pc $\times$ 100 pc) region at a resolution of 2' (0.9 pc). The Galactic radio emission in the region analyzed is almost entirely thermal (free-free) at 148 MHz, with emission measures of $10^3 < EM~{\rm[pc~cm^{-6}]} < 10^6$. As filamentary structure is a prominent feature of the emission, we use DisPerSE and FilChap to identify filamentary ridges and characterize their radial ($EM$) profiles. The distribution of radial profiles has a characteristic width of 4.3 pc and a power-law distribution ($\beta = -1.8 \pm 0.1$) in peak $EM$ down to our completeness limit of 4200 pc cm$^{-6}$. The electron densities of the filamentary structure range from $10 < n_e~{\rm[cm^{-3}]} < 400$ with a median value of 35 cm$^{-3}$, remarkably similar to [N II] surveys of ionized gas. Cyg OB2 may ionize at most two-thirds of the total ionized gas and the ionized gas in filaments. More than half of the filamentary structures are likely photoevaporating surfaces flowing into a surrounding diffuse (~5 cm$^{-3}$) medium. However, this is likely not the case for all ionized gas ridges. A characteristic width in the distribution of ionized gas points to the stellar winds of Cyg OB2 creating a fraction of the ionized filaments through swept-up ionized gas or dissipated turbulence., Comment: 19 pages, 14 figures, 1 table; accepted for publication in A&A

We present a charge distribution based model that computes the infrared spectrum of polycyclic aromatic hydrocarbon (PAH) molecules using recent measurements or quantum chemical calculations of specific PAHs. The model is applied to a sample of well-studied photodissociation regions (PDRs) with well-determined physical conditions (the radiation field strength, $G_{0}$, electron density $n_{e}$, and the gas temperature, $T_{\rm gas}$). Specifically, we modelled the emission of five PAHs ranging in size from 18 to 96 carbon atoms, over a range of physical conditions characterized by the ionization parameter $\gamma = G_{0}\times T_{\rm gas}^{1/2}/n_{e}$. The anions emerge as the dominant charge carriers in low $\gamma $ ($< 2\times 10^{2}$) environments, neutrals in the intermediate $\gamma$ ($10^{3} - 10^{4}$) environments, and cations in the high $\gamma$ ($ > 10^{5}$) environments. Furthermore, the PAH anions and cations exhibit similar spectral characteristics. The similarity in the cationic and anionic spectra translates into the interpretation of the 6.2/(11.0+11.2) band ratio, with high values of this ratio associated with large contributions from either cations or anions. The model's predicted values of 6.2/(11.0+11.2) and 3.3/6.2 compared well to the observations in the PDRs NGC 7023, NGC 2023, the horsehead nebula, the Orion bar, and the diffuse ISM, demonstrating that changes in the charge state can account for the variations in the observed PAH emission. We also reassess the diagnostic potential of the 6.2/(11.0+11.2) vs 3.3/(11.0+11.2) ratios and show that without any prior knowledge about $\gamma$, the 3.3/(11.0+11.2) can predict the PAH size, but the 6.2/(11.0+11.2) cannot predict the $\gamma$ of the astrophysical environment., Comment: 27 pages, 24 figures, Accepted for publication in MNRAS

Interest in stellar feedback has recently increased because new studies suggest that radiative and mechanical feedback from young massive stars regulate the physical and chemical composition of the interstellar medium (ISM) significantly. Recent SOFIA [CII] 158 micron observations of the Orion Veil revealed that the expanding bubble is powered by stellar winds and influenced by previously active molecular outflows of ionizing massive stars. We aim to investigate the mechanical feedback on the whole Veil shell by searching for jets/outflows interacting with the Veil shell and determining the origin/driving mechanisms of these collisions. In the light of these findings, as well as the momenta of the dents and their dynamical timescales, we propose that the dents are created by the interaction of collimated jets/outflows from protostars with luminosities ranging from 10$^3$ to 10$^4$ $L_\odot$ indicating B-type stars in the Orion star-forming cloud with the surrounding Veil shell. However, it is challenging to pinpoint the driving stars as they may have moved from the original ejection points of the jets/outflows. We conclude that the dynamics of the expanding Veil shell is influenced not just by the O-type stars in the Trapezium cluster, but also by less massive stars, especially B-type, in the Orion Nebula. Mechanical feedback from protostars with a range of masses appears to play an important role in determining the morphology of [HII] regions and injecting turbulence into the medium., Comment: 18 pages, 18 figures, Accepted by A&A on March 22, 2022

We present the Infrared Database of Extragalactic Observables from Spitzer (IDEOS), a homogeneous, publicly available, database of 77 fitted mid-infrared observables in the 5.4-36um range, comprising measurements for 3335 galaxies observed in the low-resolution staring mode of the Infrared Spectrometer onboard the Spitzer Space Telescope. Among the included observables are PAH fluxes and their equivalent widths, the strength of the 9.8um silicate feature, emission line fluxes, solid-state features, rest frame continuum fluxes, synthetic photometry, and a mid-infrared spectral classification. The IDEOS spectra were selected from the Cornell Atlas of Spitzer-IRS Sources. To our surprise we have detected at a >95% confidence level crystalline silicates in the spectra of 786 IDEOS galaxies. The detections range from single band detections to detections of all fitted crystalline bands (16, 19, 23, 28 and 33um). We find the strength of the crystalline silicate bands to correlate with the amorphous silicate strength, and the change from an emission to an absorption feature to occur at higher obscuration as the wavelength of the crystalline silicate band is longer. These observed characteristics are consistent with an origin for the amorphous and crystalline silicate features in a centrally heated dust geometry, either an edge-on disk or a cocoon. We find the 23 and 33um crystalline silicate bands to be well-suited to classify the obscuration level of galactic nuclei, even in the presence of strong circumnuclear star formation. Based on our detection statistics, we conclude that crystalline silicates are a common component of the interstellar medium of galactic nuclei., Comment: 46 pages, 34 figures, accepted for publication in ApJS. Associated spectra and data table can be downloaded from the IDEOS portal at http://ideos.astro.cornell.edu

The role of feedback in the self-regulation of star formation is a fundamental question in astrophysics. The Orion Nebula is the nearest site of ongoing and recent massive star formation. It is a unique laboratory for the study of stellar feedback. Recent SOFIA [CII] 158 $\mu$m observations revealed an expanding bubble, the Veil shell, being powered by stellar winds and ionization feedback. We have identified a protrusion-like substructure in the Northwest portion of the Orion Veil Shell that may indicate additional feedback mechanisms that are highly directional. Our goal is to investigate the origin of the protrusion by quantifying its possible driving mechanisms. We use the [CII] 158 $\mu$m map of the Orion Nebula obtained with the upGREAT instrument onboard SOFIA. The spectral and spatial resolution of the observations are 0.3 km/s and 16 arcsec, respectively. We consider three possible origins for this protrusion: Fossil outflow cavities created by jets/outflows during the protostellar accretion phase, pre-existing clumpiness in the OMC-1 core, and the stellar wind during the main sequence phase. Based on the energetics and the morphology, we conclude that the northwestern part of the pre-existing cloud was locally perturbed by outflows ejected from massive protostars in the Trapezium cluster. This suggests that the protrusion of the Veil is the result of mechanical rather than radiative feedback. Furthermore, we argue that the location of the protrusion is a suitable place to break the Orion Veil owing to the photo-ablation from the walls of the protrusion. We conclude that the outflows of massive protostars can influence the morphology of the future \hii\,region and even cause breakages in the ionization front. Specifically, the interaction of stellar winds of main-sequence stars with the molecular core pre-processed by the protostellar jet is important., Comment: Accepted by A&A on Feb 9, 2022 (20 pages, 19 figures)

Context: How do expanding HII regions interact with their environmental cloud? This is one of the central questions driving the SOFIA legacy program FEEDBACK. Here, we present a case study toward the prototypical H{\sc ii} region NGC7538. Methods: With SOFIA we mapped an area of ~210'^2 around NGC7538 in the [CII] line at 1.9THz. Complementary observed atomic carbon [CI] and high-J CO(8-7) data as well as archival NIR/FIR, cm continuum, CO(3-2) and HI data are folded into the analysis. Results: While the overall [CII] morphology follows the general ionized gas, the channel maps show multiple bubble-like structures with sizes on the order of ~80-100" (~1.0-1.28pc). While at least one of them may be an individual feedback bubble driven by the main exciting sources of the region, the other bubble-morphologies may also be due to the intrinsically porous structure of the HII region. An analysis of the expansion velocities around 10km s^{-1} indicates that thermal expansion is not sufficient but that wind-driving from the central O-stars is required. The most blue-shifted [CII] component has barely any molecular or atomic counterparts. At the interface to the molecular cloud, we find a typical photon-dominated region (PDR) with a bar-shape. Ionized, atomic and molecular carbon show a layered structure in this PDR. The carbon in the PDR is dominated by its ionized form with atomic and molecular masses of ~0.45+-0.1M_{\odot} and ~1.2+-0.1M_{\odot}, respectively, compared to the ionized carbon in the range of 3.6-9.7M_{\odot}. Conclusions: The NGC7538 HII region exhibits a diverse set of sub-structures that interact with each other as well as with the adjacent cloud. Compared to other recent [CII] observations of HII regions (e.g., Orion Veil, RCW120, RCW49), bubble-shape morphologies revealed in [CII] emission, indicative of expanding shells, are recurring structures of PDRs., Comment: 21 pages, 20 figures, accepted for Astronomy & Astrophysics, a high-resolution version can be found at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

Massive stars disrupt their natal molecular cloud material through radiative and mechanical feedback processes. These processes have profound effects on the evolution of interstellar matter in our Galaxy and throughout the Universe, from the era of vigorous star formation at redshifts of 1-3 to the present day. The dominant feedback processes can be probed by observations of the Photo-Dissociation Regions (PDRs) where the far-ultraviolet photons of massive stars create warm regions of gas and dust in the neutral atomic and molecular gas. PDR emission provides a unique tool to study in detail the physical and chemical processes that are relevant for most of the mass in inter- and circumstellar media including diffuse clouds, proto-planetary disks and molecular cloud surfaces, globules, planetary nebulae, and star-forming regions. PDR emission dominates the infrared (IR) spectra of star-forming galaxies. Most of the Galactic and extragalactic observations obtained with the James Webb Space Telescope (JWST) will therefore arise in PDR emission. In this paper we present an Early Release Science program using the MIRI, NIRSpec, and NIRCam instruments dedicated to the observations of an emblematic and nearby PDR: the Orion Bar. These early JWST observations will provide template datasets designed to identify key PDR characteristics in JWST observations. These data will serve to benchmark PDR models and extend them into the JWST era. We also present the Science-Enabling products that we will provide to the community. These template datasets and Science-Enabling products will guide the preparation of future proposals on star-forming regions in our Galaxy and beyond and will facilitate data analysis and interpretation of forthcoming JWST observations., Comment: Submitted to PASP

Revealing the 3D dynamics of HII regions and their associated molecular clouds is important for understanding the longstanding problem as to how stellar feedback affects the density structure and kinematics of the interstellar medium. We employed observations of the HII region RCW 120 in [CII], observed within the SOFIA legacy program FEEDBACK, and the $^{12}$CO and $^{13}$CO (3$\to$2) lines, obtained with APEX. In addition we used HI data from the Southern Galactic Plane Survey. Two radiative transfer models were used to fit the observed data. A line profile analysis with the 1D non-LTE radiative transfer code SimLine proves that the CO emission cannot stem from a spherically symmetric molecular cloud configuration. With a two-layer multicomponent model, we then quantified the amount of warm background and cold foreground gas. There is a deficit of CO emission along the line-of-sight toward the center of the HII region which indicates that the HII region is associated with a flattened molecular cloud. Self-absorption in the CO line may hide signatures of infalling and expanding molecular gas. The [CII] emission arises from an expanding [CII] bubble and from the PDRs. A significant part of [CII] emission is absorbed in a cool (~60-100 K), low-density (<500 cm$^{-3}$) atomic foreground layer with a thickness of a few parsec. We propose that the RCW 120 HII region formed in a flattened molecular cloud and is now bursting out of its parental cloud. The compressed surrounding molecular layer formed a torus around the spherically expanding HII bubble. This scenario can possibly be generalized for other HII bubbles and would explain the observed "flat" structure of molecular clouds associated with HII bubbles. We suggest that the [CII] absorption observed in many star-forming regions is at least partly caused by low-density, cool, HI-envelopes surrounding the molecular clouds.

The [CII] 158$\mu$m fine-structure line is the dominant cooling line of moderate-density photodissociation regions (PDRs) illuminated by moderately bright far-ultraviolet (FUV) radiation fields. We aim to understand the origin of [CII] emission and its relation to other tracers of gas and dust in PDRs. One focus is a study of the heating efficiency of interstellar gas as traced by the [CII] line to test models of the photoelectric heating of neutral gas by polycyclic aromatic hydrocarbon (PAH) molecules and very small grains. We make use of a one-square-degree map of velocity-resolved [CII] line emission toward the Orion Nebula complex, and split this out into the individual spatial components, the expanding Veil Shell, the surface of OMC4, and the PDRs associated with the compact HII region of M43 and the reflection nebula NGC 1977. We employed Herschel far-infrared photometric images to determine dust properties. Moreover, we compared with Spitzer mid-infrared photometry to trace hot dust and large molecules, and velocity-resolved IRAM 30m CO(2-1) observations of the molecular gas. The [CII] intensity is tightly correlated with PAH emission in the IRAC 8$\mu$m band and far-infrared emission from warm dust. The correlation between [CII] and CO(2-1) is very different in the four subregions and is very sensitive to the detailed geometry. Constant-density PDR models are able to reproduce the observed [CII], CO(2-1), and integrated far-infrared (FIR) intensities. We observe strong variations in the photoelectric heating efficiency in the Veil Shell behind the Orion Bar and these variations are seemingly not related to the spectral properties of the PAHs. The [CII] emission from the Orion Nebula complex stems mainly from moderately illuminated PDR surfaces. Future observations with the James Webb Space Telescope can shine light on the PAH properties that may be linked to these variations.

Chemistry along the star- and planet-formation sequence regulates how prebiotic building blocks -- carriers of the elements CHNOPS -- are incorporated into nascent planetesimals and planets. Spectral line observations across the electromagnetic spectrum are needed to fully characterize interstellar CHNOPS chemistry, yet to date there are only limited astrochemical constraints at THz frequencies. Here, we highlight advances to the study of CHNOPS astrochemistry that will be possible with the Orbiting Astronomical Satellite for Investigating Stellar Systems (OASIS). OASIS is a NASA mission concept for a space-based observatory that will utilize an inflatable 14-m reflector along with a heterodyne receiver system to observe at THz frequencies with unprecedented sensitivity and angular resolution. As part of a survey of H2O and HD towards ~100 protostellar and protoplanetary disk systems, OASIS will also obtain statistical constraints on the inventories of light hydrides including NH3 and H2S towards protoplanetary disks, as well as complex organics in protostellar hot corinos and envelopes. Line surveys of additional star-forming regions, including high-mass hot cores, protostellar outflow shocks, and prestellar cores, will also leverage the unique capabilities of OASIS to probe high-excitation organics and small hydrides, as is needed to fully understand the chemistry of these objects., Comment: Accepted to Frontiers in Astronomy and Space Sciences

We present 5--14~$\mu$m spectra at two different positions across the Orion Bar photodissociation region (PDR) obtained with the Infrared Spectrograph onboard the Spitzer Space Telescope and 3.3~$\mu$m PAH observations obtained with the Stratospheric Observatory for Infrared Astronomy (SOFIA). We aim to characterize emission from Polycyclic Aromatic Hydrocarbon (PAH), dust, atomic and molecular hydrogen, argon, sulfur, and neon as a function of distance from the primary illuminating source. We find that all the major PAH bands peak between the ionization front and the PDR front, as traced by H$_{2}$, while variations between these bands become more pronounced moving away from this peak into the face-on PDRs behind the PDR front and at the backside of the \HII\, region. While the relative PAH intensities are consistent with established PAH characteristics, we report unusual behaviours and attribute these to the PDR viewing angle and the strength of the FUV radiation field impinging on the PDRs. We determine the average PAH size which varies across the Orion Bar. We discuss subtle differences seen between the cationic PAH bands and highlight the photo-chemical evolution of carbonaceous species in this PDR environment. We find that PAHs are a good tracer of environmental properties such as the strength of the FUV radiation field and the PAH ionization parameter., Comment: 25 pages, 19 Figures. To be published in MNRAS

The infrared signature of polycyclic aromatic hydrocarbons (PAHs) are present in many protostellar disks and these speciesare thought to play an important role in heating of the gas in the photosphere. We aim to consider PAH cluster formation as one possible cause for non-detections of PAH features in protoplanetary disks. We test the necessary conditions for cluster formation and cluster dissociation by stellar optical and FUV photons in protoplanetarydisks using a Herbig Ae/Be and a T Tauri star disk model. We perform Monte-Carlo (MC) and statistical calculations to determine dissociation rates for coronene, circumcoronene and circumcoronene clusters with sizes between 2 and 200 cluster members. By applying general disk models to our Herbig Ae/Be and T Tauri star model, we estimate the formation rate of PAH dimers and compare these with the dissociation rates. We show that the formation of PAH dimers can take place in the inner 100 AU of protoplanetary disks in sub-photospheric layers. Dimer formation takes seconds to years allowing them to grow beyond dimer size in a short time. We further demonstrate that PAH cluster increase their stability while they grow if they are located beyond a critical distance that depends on stellar properties and PAH species. The comparison with the local vertical mixing time scale allows a determination of the minimum cluster size necessaryfor survival of PAH clusters. Considering the PAH cluster formation sites, cluster survival in the photosphere of the inner disk of Herbig stars isunlikely because of the high UV radiation. For the T Tauri stars, survival of coronene, circumcoronene and circumcircumcoronene clusters is possible and cluster formation should be considered as one possible explanation for low PAH detection rates in T Tauri star disks., Comment: Accepted for Astronomy and Astrophysics. 12 pages, 8 figures

The electron density ($n_{e^{-}}$) plays an important role in setting the chemistry and physics of the interstellar medium. However, measurements of $n_{e^{-}}$ in neutral clouds have been directly obtained only toward a few lines of sight or they rely on indirect determinations. We use carbon radio recombination lines and the far-infrared lines of C$^{+}$ to directly measure $n_{e^{-}}$ and the gas temperature in the envelope of the integral shaped filament (ISF) in the Orion A molecular cloud. We observed the C$102\alpha$ and C$109\alpha$ carbon radio recombination lines (CRRLs) using the Effelsberg 100m telescope at ~2' resolution toward five positions in OMC-2 and OMC-3. Since the CRRLs have similar line properties, we averaged them to increase the signal-to-noise ratio of the spectra. We compared the intensities of the averaged CRRLs, and the 158 {\mu}m-[CII] and [$^{13}$CII] lines to the predictions of a homogeneous model for the C$^{+}$/C interface in the envelope of a molecular cloud and from this comparison we determined the electron density, temperature and C$^{+}$ column density of the gas. We detect the CRRLs toward four positions, where their velocity and widths (FWHM 2.3 km s$^{-1}$) confirms that they trace the envelope of the ISF. Toward two positions we detect the CRRLs, and the [CII] and [$^{13}$CII] lines with a signal-to-noise ratio >5, and we find $n_{e^{-}}=0.65\pm0.12$ cm$^{-3}$ and $0.95\pm0.02$ cm$^{-3}$, which corresponds to a gas density $n_{H}\approx5\times10^{3}$ cm$^{-3}$ and a thermal pressure of $p_{th}\approx4\times10^{5}$ K cm$^{-3}$. We also constrained the ionization fraction in the denser portions of the molecular cloud using the HCN(1-0) and C$_{2}$H(1-0) lines to $x(e^{-})<3\times10^{-6}$. The derived electron densities and ionization fraction imply that $x(e^{-})$ drops by a factor >100 between the C$^{+}$ layer and the regions probed by HCN(1-0)., Comment: 14 pages, 6 figures, 8 tables, accepted for publication in A&A

Context. The [CII] 158micron far-infrared fine-structure line is one of the dominant cooling lines of the star-forming interstellar medium (ISM). Hence [CII] emission originates in and thus can be used to trace a range of ISM processes. Velocity-resolved large-scale mapping of [CII] in star-forming regions provides a unique perspective of the kinematics of these regions and their interactions with the exciting source of radiation. Aims. We explore the scientific applications of large-scale mapping of velocity-resolved [CII] observations. With the [CII] observations, we investigate the effect of stellar feedback on the ISM. We present the details of observation, calibration, and data reduction using a heterodyne array receiver mounted on an airborne observatory. Results. A square-degree [CII] map with a spectral resolution of 0.3 km/s is presented. The scientific potential of this data is summarized with discussion of mechanical and radiative stellar feedback, filament tracing using [CII], [CII] opacity effects, [CII] and carbon recombination lines, and [CII] interaction with the large molecular cloud. The data quality and calibration is discussed in detail, and new techniques are presented to mitigate the effects of unavoidable instrument deficiencies (e.g. baseline stability) and thus to improve the data quality. A comparison with a smaller [CII] map taken with the Herschel/Heterodyne Instrument for the Far-Infrared (HIFI) spectrometer is presented.

We present observations from the First Light Infrared TEst CAMera (FLITECAM) on board the Stratospheric Observatory for Infrared Astronomy (SOFIA), the Spitzer Infrared Array Camera (IRAC) and the Spitzer Infrared Spectrograph (IRS) SH mode in three well-known Photodissocation Regions (PDRs), the reflection nebulae (RNe) NGC 7023 and NGC 2023 and to the southeast of the Orion Bar, which are well suited to probe emission from Polycyclic Aromatic Hydrocarbon molecules (PAHs). We investigate the spatial behaviour of the FLITECAM 3.3 um filter as a proxy for the 3.3 um PAH band, the integrated 11.2 um PAH band, and the IRAC 8.0 um filter as a proxy for the sum of the 7.7 and 8.6 um PAH bands. The resulting ratios of 11.2/3.3 and IRAC 8.0/11.2 provide an approximate measure of the average PAH size and PAH ionization respectively. In both RNe, we find that the relative PAH ionization and the average PAH size increases with decreasing distance to the illuminating source. The average PAH sizes derived for NGC 2023 are greater than those found for NGC 7023 at all points. Both results indicate that PAH size is dependent on the radiation field intensity. These results provide additional evidence of a rich carbon-based chemistry driven by the photo-chemical evolution of the omnipresent PAH molecules within the interstellar medium. In contrast, we did not detect a significant variation in the average PAH size found in the region southeast of the Orion Bar and report a peculiar PAH ionization radial profile., Comment: 20 pages, 11 figures, accepted for publication in ApJ

The [CII] $158\,\mu\mathrm{m}$ fine-structure line is one of the dominant coolants of the neutral interstellar medium. It is hence one of the brightest far-infrared emission lines and can be observed not only in star-forming regions throughout the Galaxy, but also in the diffuse interstellar medium and in distant galaxies. [CII] line emission has been suggested to be a powerful tracer of star-formation. We aim to understand the origin of [CII] emission and its relation to other tracers of interstellar gas and dust. This includes a study of the heating efficiency of interstellar gas as traced by the [CII] line to test models of gas heating. We make use of a one-square-degree map of velocity-resolved [CII] line emission towards the Orion Nebula complex, including M43 and NGC 1977. The [CII] intensity is tightly correlated with PAH emission in the IRAC $8\,\mu\mathrm{m}$ band and far-infrared emission from warm dust. The correlation between [CII] and CO(2-1) is affected by the detailed geometry of the region. We find particularly low [CII]-over-FIR intensity ratios towards large columns of (warm and cold) dust, which suggest the interpretation of the "[CII] deficit" in terms of a "FIR excess". A slight decrease in the FIR line-over-continuum intensity ratio can be attributed to a decreased heating efficiency of the gas. We find that, at the mapped spatial scales, predictions of the star-formation rate from [CII] emission underestimate the star-formation rate calculated from YSO counts in the Orion Nebula complex by an order of magnitude. [CII] emission from the Orion Nebula complex arises dominantly in the cloud surfaces, many viewed in edge-on geometry. [CII] emission from extended faint cloud surfaces may contribute significantly to the total [CII] emission on galactic scales.

Radiative and mechanical feedback of massive stars regulates star formation and galaxy evolution. Positive feedback triggers the creation of new stars by collecting dense shells of gas, while negative feedback disrupts star formation by shredding molecular clouds. Although key to understanding star formation, their relative importance is unknown. Here, we report velocity-resolved observations from the SOFIA (Stratospheric Observatory for Infrared Astronomy) legacy program FEEDBACK of the massive star-forming region RCW 120 in the [CII] 1.9-THz fine-structure line, revealing a gas shell expanding at 15 km/s. Complementary APEX (Atacama Pathfinder Experiment) CO J=3-2 345-GHz observations exhibit a ring structure of molecular gas, fragmented into clumps that are actively forming stars. Our observations demonstrate that triggered star formation can occur on much shorter time scales than hitherto thought (<0.15 million years), suggesting that positive feedback operates on short time periods., Comment: 8 Pages, 4 Figures, 1 Table (+16 Page Supplementary Materials with 10 Figures). Published in Science Advances on April 9th, 2021, with free open access to the online article available at https://advances.sciencemag.org/content/7/15/eabe9511

We unveil the stellar wind driven shell of the luminous massive star-forming region of RCW 49 using SOFIA FEEDBACK observations of the [CII] 158 $\mu$m line. The complementary dataset of the $^{12}$CO and $^{13}$CO J = 3 - 2 transitions is observed by the APEX telescope and probes the dense gas toward RCW 49. Using the spatial and spectral resolution provided by the SOFIA and APEX telescopes, we disentangle the shell from a complex set of individual components of gas centered around RCW 49. We find that the shell of radius ~ 6 pc is expanding at a velocity of 13 km s$^{-1}$ toward the observer. Comparing our observed data with the ancillary data at X-Ray, infrared, sub-millimeter and radio wavelengths, we investigate the morphology of the region. The shell has a well defined eastern arc, while the western side is blown open and is venting plasma further into the west. Though the stellar cluster, which is ~ 2 Myr old gave rise to the shell, it only gained momentum relatively recently as we calculate the shell's expansion lifetime ~ 0.27 Myr, making the Wolf-Rayet star WR20a a likely candidate responsible for the shell's re-acceleration., Comment: 31 pages, 17 figures

We describe and discuss remarkable infrared spectra, covering key portions of the $2-5$ $\mu$m wavelength interval, of the probable OH/IR supergiant 2MASS J17470898$-$2829561 (2M1747), located in direction of the Sgr B molecular cloud complex within the Central Molecular Zone (CMZ) of the Galaxy. This star was originally singled out for examination based on its suitability for spectroscopy of lines of H$_3^+$ in the CMZ. Analysis of the spectra shows that 2M1747 is deeply embedded within Sgr B1, with A$_V$ $\gtrsim$ 100 mag, making it the only star within Sgr B for which infrared spectra have been obtained at present, and thereby a unique infrared probe of the dense interstellar medium within the CMZ. Despite the high extinction, spectra of 2M1747 reveal a veiled photosphere in the $K$ band and circumstellar gas in the $M$ band, giving clues as to its nature. Its $ 3.5-4.0$ $\mu$m spectrum contains the strongest absorption lines of H$_3^+$ observed toward any object to date. The $4.5-4.8$ $\mu$m spectrum has impressively deep and wide absorption lines of interstellar CO, most of which arise in dense gas within Sgr B1. The $3-5$ $\mu$m spectrum also contains several solid state absorption features, which are characteristic of both dense and diffuse clouds, and which raise questions about the identifications of some of these features. We discuss the nature of the star, the extinction to it, the extinction law for dust in the CMZ, and the identifications of the various solid-state features and where they are produced along this complex line of sight., Comment: 17 pages, 10 figures; accepted by ApJ 2021 March 10

We present the gas-phase infrared spectra of the phenyl cation, phenylium, in its perprotio C$_6$H$_5^+$ and perdeutero C$_6$D$_5^+$ forms, in the 260-1925 cm$^{-1}$ (5.2-38 $\mu$m) spectral range, and investigate the observed photofragmentation. The spectral and fragmentation data were obtained using Infrared Multiple Photon Dissociation (IRMPD) spectroscopy within a Fourier Transform Ion Cyclotron Resonance Mass Spectrometer (FTICR MS) located inside the cavity of the free electron laser FELICE (Free Electron Laser for Intra-Cavity Experiments). The $^1$A$_1$ singlet nature of the phenylium ion is ascertained by comparison of the observed IR spectrum with DFT calculations, using both harmonic and anharmonic frequency calculations. To investigate the observed loss of predominantly [2C,nH] (n=2-4) fragments, we explored the potential energy surface (PES) to unravel possible isomerization and fragmentation reaction pathways. The lowest energy pathways toward fragmentation include direct H elimination, and a combination of facile ring-opening mechanisms ($\leq2.4$ eV), followed by elimination of H or CCH$_2$. Energetically, all H-loss channels found are more easily accessible than CCH$_2$-loss. Calculations of the vibrational density of states for the various intermediates show that at high internal energies, ring opening is the thermodynamically the most advantageous, eliminating direct H-loss as a competing process. The observed loss of primarily [2C,2H] can be explained through entropy calculations that show favored loss of [2C,2H] at higher internal energies., Comment: 10 pages, 6 figures and 1 table. Accepted for publication on PCCP on 04 Feb 2021; supplementary Material available on PCCP website

FEEDBACK is a SOFIA legacy program dedicated to study the interaction of massive stars with their environment. It performs a survey of 11 galactic high mass star forming regions in the 158 $\mu$m (1.9 THz) line of CII and the 63 $\mu$m (4.7 THz) line of OI. We employ the 14 pixel LFA and 7 pixel HFA upGREAT instrument to spectrally resolve (0.24 MHz) these FIR structure lines. With an observing time of 96h, we will cover $\sim$6700 arcmin$^2$ at 14.1$''$ angular resolution for the CII line and 6.3$''$ for the OI line. The observations started in spring 2019 (Cycle 7). Our aim is to understand the dynamics in regions dominated by different feedback processes from massive stars such as stellar winds, thermal expansion, and radiation pressure, and to quantify the mechanical energy injection and radiative heating efficiency. The CII line provides the kinematics of the gas and is one of the dominant cooling lines of gas for low to moderate densities and UV fields. The OI line traces warm and high-density gas, excited in photodissociations regions with a strong UV field or by shocks. The source sample spans a broad range in stellar characteristics from single OB stars, to small groups of O stars, to rich young stellar clusters, to ministarburst complexes. It contains well-known targets such as Aquila, the Cygnus X region, M16, M17, NGC7538, NGC6334, Vela, and W43 as well as a selection of HII region bubbles, namely RCW49, RCW79, and RCW120. These CII maps, together with the less explored OI 63 $\mu$m line, provide an outstanding database for the community. They will be made publically available and will trigger further studies and follow-up observations.