CO Line Research Articles

Filamentary Molecular Cloud Formation via Collision-induced Magnetic Reconnection in a Cold Neutral Medium

We have investigated the possibility of molecular cloud formation via the collision-induced magnetic reconnection (CMR) mechanism of the cold neutral medium (CNM). Two atomic gas clouds with conditions typical of the CNM were set to collide at the interface of reverse magnetic fields. The cloud–cloud collision triggered magnetic reconnection and produced a giant 20 pc filamentary structure that was not seen in the control models without CMR. The cloud, with rich fiber-like substructures, developed a fully molecular spine at 5 Myr. Radiative transfer modeling of dust emission at far-infrared wavelengths showed that the middle part of the filament contained dense cores over a span of 5 pc. Some of the cores were actively forming stars and typically exhibited both connecting fibers in dust emission and high-velocity gas in CO line emission, indicative of active accretion through streamers. Supersonic turbulence was present in and around the CMR filament due to inflowing gas moving at supersonic velocities in the collision midplane. The shocked gas was condensed and transported to the main filament piece by piece by reconnected fields, making the filament and star formation a bottom-up process. Instead of forming a gravitationally bounded cloud that then fragments hierarchically (top-down) and forms stars, the CMR process creates dense gas pieces and magnetically transports them to the central axis to constitute the filament. Since no turbulence is manually driven, our results suggest that CMR is capable of self-generating turbulence. Finally, the resulting helical field should show field reversal on both sides of the filament from most viewing angles.

Millimetre observations of the S-type AGB star chi Cygni: Variability of the emission of the inner envelope

New observations are presented of millimetre line emissions of the circumstellar envelope (CSE) of the asymptotic giant branch (AGB) star chi Cygni, using the recently upgraded NOEMA array. chi Cygni is an S-type Mira variable, at the border between oxygen-rich and carbon-rich stars. It has been observed for over 40 years to display features that suggest evidence for the strong role played by pulsation-associated shock waves in the generation of its wind. These new observations provide evidence of a bright H12CN(3-2) line emission confined to the very close neighbourhood of the star; however, this emission appears significantly more extended in 2024 than in 2023. The interpretation of such variability in terms of maser emission has been considered and found to raise significant unanswered questions. Moreover, other unexpected features are observed in the very close neighbourhood of the star, including low Si16O(6-5)/Si17O(6-5), 28SiO(5-4)/29SiO(5-4), and 12CO(2-1)/13CO(2-1) line emission ratios. We discuss several features, which possibly confirm the important role played by shocks: a measurement of the SiO(5-4)/SiO(6-5)emission ratio; the observation of a recent mass ejection, particularly enhanced in the north-western red-shifted octant, which has left a depression in its wake; patterns of enhanced CO(2-1) line emission, which suggest an interpretation in terms of episodic outflows, on a time scale of a few decades, enhanced over solid angles associated with the surface of convective cells. Unravelling the mechanisms underlying such newly observed features is very challenging. Thus, a confirmation of the reported observations with improved sensitivity and angular resolution would be highly welcome. The observation of SiO maser emission in the (nu =1,$J$=6-5) transition is reported for the first time.

IPA: Class 0 Protostars Viewed in CO Emission Using JWST

We investigate the bright CO fundamental emission in the central regions of five protostars in their primary mass assembly phase using new observations from JWST’s Near-Infrared Spectrograph and Mid-Infrared Instrument. CO line emission images and fluxes are extracted for a forest of ∼150 rovibrational transitions from two vibrational bands, v = 1−0 and v = 2−1. However, 13CO is undetected, indicating that 12CO emission is optically thin. We use H2 emission lines to correct fluxes for extinction and then construct rotation diagrams for the CO lines with the highest spectral resolution and sensitivity to estimate rotational temperatures and numbers of CO molecules. Two distinct rotational temperature components are required for v = 1 (∼600 to 1000 K and 2000 to ∼104 K), while one hotter component is required for v = 2 (≳3500 K). 13CO is depleted compared to the abundances found in the interstellar medium, indicating selective UV photodissociation of 13CO; therefore, UV radiative pumping may explain the higher rotational temperatures in v = 2. The average vibrational temperature is ∼1000 K for our sources and is similar to the lowest rotational temperature components. Using the measured rotational and vibrational temperatures to infer a total number of CO molecules, we find that the total gas masses range from lower limits of ∼1022 g for the lowest mass protostars to ∼1026 g for the highest mass protostars. Our gas mass lower limits are compatible with those in more evolved systems, which suggest the lowest rotational temperature component comes from the inner disk, scattered into our line of sight, but we also cannot exclude the contribution to the CO emission from disk winds for higher mass targets.

Early Planet Formation in Embedded Disks. XI. A High-resolution View Toward the BHR 71 Class 0 Protostellar Wide Binary

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the binary Class 0 protostellar system BHR 71 IRS1 and IRS2 as part of the Early Planet Formation in Embedded Disks (eDisk) ALMA Large Program. We describe the 12CO (J = 2–1), 13CO (J = 2–1), C18O (J = 2–1), H2CO (J = 32,1–22,0), and SiO (J = 5–4) molecular lines along with the 1.3 mm continuum at high spatial resolution (∼0.″08 or ∼5 au). Dust continuum emission is detected toward BHR 71 IRS1 and IRS2, with a central compact component and extended continuum emission. The compact components are smooth and show no sign of substructures such as spirals, rings, or gaps. However, there is a brightness asymmetry along the minor axis of the presumed disk in IRS1, possibly indicative of an inclined geometrically and optically thick disk-like component. Using a position–velocity diagram analysis of the C18O line, clear Keplerian motions were not detected toward either source. If Keplerian rotationally supported disks are present, they are likely deeply embedded in their envelope. However, we can set upper limits of the central protostellar mass of 0.46 M ⊙ and 0.26 M ⊙ for BHR 71 IRS1 and BHR 71 IRS2, respectively. Outflows traced by 12CO and SiO are detected in both sources. The outflows can be divided into two components, a wide-angle outflow and a jet. In IRS1, the jet exhibits a double helical structure, reflecting the removal of angular momentum from the system. In IRS2, the jet is very collimated and shows a chain of knots, suggesting episodic accretion events.

Probing the physics of star formation (ProPStar)

Context. Turbulence is a key component of molecular cloud structure. It is usually described by a cascade of energy down to the dissipation scale. The power spectrum for subsonic incompressible turbulence is ∝k−5/3, while for supersonic turbulence it is ∝k−2. Aims. We determine the power spectrum in an actively star-forming molecular cloud, from parsec scales down to the expected magnetohydrodynamic (MHD) wave cutoff (dissipation scale). Methods. We analyzed observations of the nearby NGC 1333 star-forming region in three different tracers to cover the different scales from ∼10 pc down to 20 mpc. The largest scales are covered with the low-density gas tracer 13CO (1–0) obtained with a single dish, the intermediate scales are covered with single-dish observations of the C18O (3–2) line, while the smallest scales are covered in H13CO+ (1–0) and HNC (1–0) with a combination of NOEMA interferometer and IRAM 30m single-dish observations. The complementarity of these observations enables us to generate a combined power spectrum covering more than two orders of magnitude in spatial scale. Results. We derive the power spectrum in an active star-forming region spanning more than 2 decades of spatial scales. The power spectrum of the intensity maps shows a single power-law behavior, with an exponent of 2.9 ± 0.1 and no evidence of dissipation. Moreover, there is evidence that the power spectrum of the ions to have more power at smaller scales than the neutrals, which is opposite to the theoretical expectations. Conclusions. We show new possibilities for studying the dissipation of energy at small scales in star-forming regions provided by interferometric observations.

Molecular clouds as hubs in spiral galaxies: gas inflow and evolutionary sequence

ABSTRACT We decomposed the molecular gas in the spiral galaxy NGC 628 (M74) into multiscale hub-filament structures using the CO (2$-$1) line by the dendrogram algorithm. All leaf structures as potential hubs were classified into three categories, i.e. leaf-HFs-A, leaf-HFs-B and leaf-HFs-C. Leaf-HFs-A exhibit the best hub-filament morphology, which also have the highest density contrast, the largest mass and the lowest virial ratio. We employed the filfinder algorithm to identify and characterize filaments within 185 leaf-HFs-A structures, and fitted the velocity gradients around the intensity peaks. Measurements of velocity gradients provide evidence for gas inflow within these structures, which can serve as a kinematic evidence that these structures are hub-filament structures. The numbers of the associated 21 μm and H α structures and the peak intensities of 7.7 μm, 21 μm, and H α emissions decrease from leaf-HFs-A to leaf-HFs-C. The spatial separations between the intensity peaks of CO and 21 μm structures of leaf-HFs-A are larger than those of leaf-HFs-C. These evidence indicate that leaf-HFs-A are more evolved than leaf-HFs-C. There may be an evolutionary sequence from leaf-HFs-C to leaf-HFs-A. Currently, leaf-HFs-C lack a distinct gravitational collapse process that would result in a significant density contrast. The density contrast can effectively measure the extent of the gravitational collapse and the depth of the gravitational potential of the structure which, in turn, shapes the hub-filament morphology. Combined with the kinematic analysis presented in previous studies, a picture emerges that molecular gas in spiral galaxies is organized into network structures through the gravitational coupling of multiscale hub-filament structures. Molecular clouds, acting as knots within these networks, serve as hubs, which are local gravitational centres and the main sites of star formation.

H_3^+ absorption and emission in local (U)LIRGs with JWST NIRSpec: Evidence for high H_2 ionization rates

We study the $3.4-4.4$ fundamental rovibrational band of H$_3^+$, a key tracer of the ionization of the molecular interstellar medium (ISM), in a sample of 12 local ($d$$<$400\,Mpc) (ultra)luminous infrared galaxies ((U)LIRGs) observed with JWST NIRSpec. The P, Q, and R branches of the band are detected in 13 out of 20 analyzed regions within these (U)LIRGs, which increases the number of extragalactic H$_3^+$ detections by a factor of 6. For the first time in the ISM, the H$_3^+$ band is observed in emission; we detect this emission in three regions. In the remaining ten regions, the band is seen in absorption. The absorptions are produced toward the $3.4-4.4$ hot dust continuum rather than toward the stellar continuum, indicating that they likely originate in clouds associated with the dust continuum source. The H$_3^+$ band is undetected in Seyfert-like (U)LIRGs where the mildly obscured X-ray radiation from the active galactic nuclei might limit the abundance of this molecule. For the detections, the H$_3^+$ abundances, $N($H$_3^+$) H (0.5--5.5)times 10$^ imply relatively high ionization rates, $ H_2 $, of between 3times 10$^ $ and $>$4times 10$^ $, which is lower than the molecular outflow velocities measured using other tracers such as OH\,119 or rotational CO lines. This suggests that H$_3^+$ traces gas close to the outflow-launching sites before it has been fully accelerated. We used nonlocal thermodynamic equilibrium models to investigate the physical conditions of these clouds. In seven out of ten objects, the H$_3^+$ excitation is consistent with inelastic collisions with H$_2$ in warm translucent molecular clouds kin $sim 250--500\,K and H_2 )$sim 10$^ $). In three objects, dominant infrared pumping excitation is required to explain the absorptions from the (3,0) and (2,1) levels of H$_3^+$ detected for the first time in the ISM.

Comparison of NH3 and 12CO, 13CO, C18O Molecular Lines in the Aquila Rift Cloud Complex

The observations of the Aquila Rift cloud complex at 23.708 and 115.271 GHz made using the Nanshan 26 m radio telescope and the 13.7 m millimeter-wavelength telescope are presented. We find that the CO(1 − 0) gas distribution is similar to the NH3 gas distribution in the Aquila Rift cloud complex. In some diffusion regions characterized by CO, we identified several dense clumps based on the distribution of detected ammonia molecular emission. Through the comparison of spectral line parameters for NH3, 13CO, and C18O, our study reveals that the line center velocities of the NH3, 13CO, and C18O lines are comparable and positively correlated, indicating that they originate from the same emission region. No significant correlation was identified for other parameters, including integrated intensity, line widths, main beam brightness temperature, as well as the column densities of NH3, 13CO, and C18O. The absolute difference in line-center velocities between the 13CO and NH3 lines is less than both the average line width of NH3 and that of 13CO. This suggests that there are no significant movements of NH3 clumps in relation to their envelopes. The velocity deviation is likely due to turbulent activity within the clumps.

The Formation of Milky Way “Bones”: Ubiquitous HI Narrow Self-absorption Associated with CO Emission

Long and skinny molecular filaments running along Galactic spiral arms are known as “bones,” since they make up the skeleton of the Milky Way. However, their origin is still an open question. Here, we compare spectral images of HI taken by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) with archival CO and Herschel dust emission to investigate the conversion from HI to H2 in two typical Galactic bones, CFG028.68-0.28 and CFG047.06+0.26. Sensitive FAST HI images and an improved methodology enabled us to extract HI narrow self-absorption (HINSA) features associated with CO line emission on and off the filaments, revealing the ubiquity of HINSA toward distant clouds for the first time. The derived cold HI abundances, [HI]/[H2], of the two bones range from ∼(0.5 to 44.7) × 10−3, which reveal different degrees of HI–H2 conversion, and are similar to those of nearby, low-mass star-forming clouds, Planck Galactic cold clumps, and a nearby active high-mass star-forming region G176.51+00.20. The HI–H2 conversion has been ongoing for 2.2–13.2 Myr in the bones, a timescale comparable to that of massive star formation therein. Therefore, we are witnessing young giant molecular clouds (GMCs) with rapid massive star formation. Our study paves the way of using HINSA to study cloud formation in Galactic bones and, more generally, in distant GMCs in the FAST era.

Coarsely scannable comb-locked ECDL with fixed local oscillator for CRDS measurement of CO line strength

Coarsely scannable comb-locked ECDL with fixed local oscillator for CRDS measurement of CO line strength

Intrinsic line profiles for X-ray fluorescent lines in SKIRT

X-ray microcalorimeter instruments are expected to spectrally resolve the intrinsic line shapes of the strongest fluorescent lines. X-ray models should therefore incorporate these intrinsic line profiles to obtain meaningful constraints from observational data. We included the intrinsic line profiles of the strongest fluorescent lines in the X-ray radiative transfer code SKIRT to model the cold-gas structure and kinematics based on high-resolution line observations from XRISM/Resolve and Athena /X-IFU. The intrinsic line profiles of the $ K and $ K lines of Cr, Mn, Fe, Co, Ni, and Cu were implemented based on a multi-Lorentzian parameterisation. Line energies are sampled from these Lorentzian components during the radiative transfer routine. In the optically thin regime, the SKIRT results match the intrinsic line profiles as measured in the laboratory. With a more complex 3D model that also includes kinematics, we find that the intrinsic line profiles are broadened and shifted to an extent that will be detectable with XRISM/Resolve; this model also demonstrates the importance of the intrinsic line shapes for constraining kinematics. We find that observed line profiles directly trace the cold-gas kinematics, without any additional radiative transfer effects. With the advent of the first XRISM/Resolve data, this update to the X-ray radiative transfer framework of SKIRT is timely and provides a unique tool for constraining the velocity structure of cold gas from X-ray microcalorimeter spectra.

An Hi-absorption-selected Cold Rotating Disk Galaxy at z ≈ 2.193

We have used the Atacama Large Millimeter/submillimeter Array to map CO(3–2) emission from a galaxy, DLA-B1228g, associated with the high-metallicity damped Lyα absorber at z ≈ 2.1929 toward the QSO PKS B1228–113. At an angular resolution of ≈0.″32 × 0.″24, DLA-B1228g shows extended CO(3–2) emission with a deconvolved size of ≈0.″78 × 0.″18, i.e., a spatial extent of ≈6.4 kpc. We detect extended stellar emission from DLA-B1228g in a Hubble Space Telescope Wide Field Camera 3 F160W image and find that Hα emission is detected in a Very Large Telescope SINFONI image from only one side of the galaxy. While the clumpy nature of the F160W emission and the offset between the kinematic and physical centers of the CO(3–2) emission are consistent with a merger scenario, this appears unlikely due to the lack of strong Hα emission, the symmetric double-peaked CO(3–2) line profile, the high molecular gas depletion timescale, and the similar velocity dispersions in the two halves of the CO(3–2) image. Kinematic modeling reveals that the CO(3–2) emission is consistent with arising from an axisymmetric rotating disk with an exponential profile, a rotation velocity of v rot = 328 ± 7 km s−1, and a velocity dispersion of σ v = 62 ± 7 km s−1. The high value of the ratio v rot/σ v , ≈5.3, implies that DLA-B1228g is a rotation-dominated cold disk galaxy, the second case of a high-z Hi-absorption-selected galaxy identified with a cold rotating disk. We obtain a dynamical mass of M dyn = (1.5 ± 0.1) × 1011 M ⊙, similar to the molecular gas mass of ≈1011 M ⊙ inferred from earlier CO(1–0) studies; this implies that the galaxy is baryon-dominated in its inner regions.

Multiline observations of hydrogen, helium, and carbon radio-recombination lines toward Orion A: A detailed dynamical study and direct determination of physical conditions

We present a study of hydrogen, helium, and carbon millimeter-wave radio-recombination lines (RRLs) toward 10 representative positions throughout the Orion Nebula complex, using the Yebes 40 m telescope in the Q band (31.3 GHz to 50.6 GHz) at an angular resolution of about 45″ (~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 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 [C II] 158 µm line and IRAM 30 m observations of the 13CO (J = 2−1) line (the complete map is presented here for the first time). We observe small differences in peak velocities between the different tracers, which cannot always be attributed to geometry but potentially to shear motions. Using the far-infrared [C II] and [13C II] intensities with the carbon RRL intensities, we can infer physical conditions (electron temperature Te and electron density ne, converted to hydrogen nuclei density nH by dividing by the carbon gas-phase abundance 𝒜c ≃ 1.4 × 10−4) in the PDR gas using nonlocal thermal equilibrium excitation models. For positions in OMC1, we infer ne ≃ 20–40 cm−3 and Te ≃ 210–240 K. On the border between OMC1 and M43, we observe two gas components with ne ≃ 2 cm−3 and ne ≃ 8 cm−3, and Te ≃ 100 K and Te ≃ 150 K. In M43, we infer ne ≃ 2–3 cm−3 and Te ≃ 140 K. The Extended Orion Nebula southeast of OMC1 is characterized by ne ≃ 2 cm−3 and Te ≃ 180 K, while OMC3 has ne ≃ 1 cm−3 and Te ≃ 130 K. Our observations are sensitive enough to detect faint lines toward two positions in OMC1, in the BN/KL PDR and the PDR close to the Trapezium stars, 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 [C II] and 13CO line widths, are broader than thermal, indicating significant turbulence in the interstellar medium, which transitions from super-Alfvénic and subsonic in the ionized gas to sub-Alfvénic 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.

Probing the dynamical and kinematical structures of detached shells around AGB stars

Context. The chemical evolution of asymptotic giant branch (AGB) stars is driven by repeated thermal pulses (TPs). The duration of a TP is only a few hundred years, whereas an inter-pulse period lasts 104 − 105 yr. Direct observations of TPs are hence unlikely. However, the detached shells seen in CO line emission that are formed as a result of a TP provide indirect constraints on the changes experienced by the star during the pulse. Aims. We aim to resolve the spatial and kinematic sub-structures in five detached-shell sources to provide detailed constraints for hydrodynamic models that describe the formation and evolution of the shells. Methods. We used observations of the 12CO (1 − 0) emission towards five carbon-AGB stars with ALMA (Atacama Large Millimeter/submillimeter Array), including previously published observations of the carbon AGB star U Ant. The data have angular resolutions of 0″​​.3 to 1″ and a velocity resolution of 0.3 km s−1. This enabled us to quantify spatial and kinematic structures in the shells. Combining the ALMA data with single-dish observations of the 12CO (1 − 0) to 12CO (4 − 3) emission towards the sources, we used radiative transfer models to compare the observed structures with previous estimates of the shell masses and temperatures. Results. The observed emission is separated into two distinct components: a more coherent, bright outer shell and a more filamentary, fainter inner shell. The kinematic information shows that the inner sub-shells move at a higher velocity relative to the outer sub-shells. The observed sub-structures reveal a negative velocity gradient outwards across the detached shells, confirming the predictions from hydrodynamical models. However, the models do not predict a double-shell structure, and the CO emission likely only traces the inner and outer edges of the shell, implying a lack of CO in the middle layers of the detached shell. Previous estimates of the masses and temperatures are consistent with originating mainly from the brighter subshell, but the total shell masses are likely lower limits. Also, additional structures in the form of partial shells outside the detached shell around V644 Sco, arcs within the shell of R Scl, and a partially filled shell for DR Ser indicate a more complicated evolution of the shells and mass-loss process throughout the TP cycle than previously assumed. Conclusions. The observed spatial and kinematical splittings of the shells appear consistent with results from the hydrodynamical models, provided the CO emission does not trace the H2 density distribution in the shell but rather traces the edges of the shells. The hydrodynamical models predict very different density profiles depending on the evolution of the shells and the different physical processes involved in the wind-wind interaction (e.g. heating and cooling processes). It is therefore not possible to constrain the total shell mass based on the CO observations alone. Additional features outside and inside the shells complicate the interpretation of the data. Complementary observations of, for example, CI as a dissociation product of CO would be necessary to understand the distribution of CO compared to H2, in addition to new detailed hydrodynamical models of the pre-pulse, pulse, and post-pulse wind. Only a comprehensive combination of observations and models will allow us to constrain the evolution of the shells and the changes in the star during the thermal-pulse cycle.

Molecular gas scaling relations for local star-forming galaxies in the low-M* regime

We derived molecular gas fractions (fmol = Mmol/M*) and depletion times (τmol = Mmol/SFR) for 353 galaxies representative of the local star-forming population with 108.5 M⊙ < M* < 1010.5 M⊙ drawn from the ALLSMOG and xCOLDGASS surveys of CO(2−1) and CO(1−0) line emission. By adding constraints from low-mass galaxies and upper limits for CO non-detections, we find the median molecular gas fraction of the local star-forming population to be constant at log fmol = −0.99−0.19+0.22, challenging previous reports of increased molecular gas fractions in low-mass galaxies. Above M* ∼ 1010.5 M⊙, we find the fmol versus M* relation to be sensitive to the selection criteria for star-forming galaxies. We tested the robustness of our results against different prescriptions for the CO-to-H2 conversion factor and different selection criteria for star-forming galaxies. The depletion timescale τmol weakly depends on M*, following a power law with a best-fit slope of 0.16 ± 0.03. This suggests that small variations in specific star formation rate (sSFR = SFR/M*) across the local main sequence of star-forming galaxies with M* < 1010.5 M⊙ are mainly driven by differences in the efficiency of converting the available molecular gas into stars. We tested these results against a possible dependence of fmol and τmol on the surrounding (group) environment of the targets by splitting them into centrals, satellites, and isolated galaxies, and find no significant variation between these populations. We conclude that the group environment is unlikely to have a large systematic effect on the molecular gas content of star-forming galaxies in the local Universe.

Gas Kinematics and Dynamics of Carina Pillars: A Case Study of G287.76-0.87

We study the kinematics of a pillar, namely G287.76-0.87, using three rotational lines of 12CO(5-4), 12CO(8-7), 12CO(11-10), and a fine structure line of [O i] 63 μm in southern Carina observed by SOFIA/GREAT. This pillar is irradiated by the associated massive star cluster Trumpler 16, which includes η Carina. Our analysis shows that the relative velocity of the pillar with respect to this ionization source is small, ∼1 km s−1, and the gas motion in the tail is more turbulent than in the head. We also performed analytical calculations to estimate the gas column density in local thermal equilibrium (LTE) conditions, which yields N CO as (∼0.2–5) × 1017 cm−2. We further constrain the gas’s physical properties in non-LTE conditions using RADEX. The non-LTE estimations result in nH2≃105cm−3 and N CO ≃ 1016 cm−2. We found that the thermal pressure within the G287.76-0.87 pillar is sufficiently high to make it stable for the surrounding hot gas and radiation feedback if the winds are not active. While they are active, stellar winds from the clustered stars sculpt the surrounding molecular cloud into pillars within the giant bubble around η Carina.

Physical properties of the southwest outflow streamer in the starburst galaxy NGC 253 with ALCHEMI

Aims. The physical properties of galactic molecular outflows are important as they could constrain outflow formation mechanisms. In this work, we study the properties of the southwest (SW) outflow streamer including gas kinematics, optical depth, dense gas fraction, and shock strength through molecular emission in the central molecular zone of the starburst galaxy NGC 253. Methods. We imaged the molecular emission in NGC 253 at a spatial resolution of 1.6″(∼27 pc at D ∼ 3.5 Mpc) based on data from the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) large program. We traced the velocity and velocity dispersion of molecular gas with the CO(1–0) line and studied the molecular spectra in the region of the SW streamer, the brightest CO streamer in NGC 253. We constrained the optical depth of the CO emission with the CO/13CO(1–0) ratio, the dense gas fraction with the HCN/CO(1–0), H13CN/13CO(1–0) and N2H+/13CO(1–0) ratios, as well as the shock strength with the SiO(2–1)/13CO(1–0) and CH3OH(2k–1k)/13CO(1–0) ratios. Results. The CO/13CO(1–0) integrated intensity ratio is ∼21 in the SW streamer region, which approximates the C/13C isotopic abundance ratio. The higher integrated intensity ratio compared to the disk can be attributed to the optically thinner environment of CO(1–0) emission inside the SW streamer. The HCN/CO(1–0) and SiO(2–1)/13CO(1–0) integrated intensity ratios both approach ∼0.2 in three giant molecular clouds (GMCs) at the base of the outflow streamers, which implies a higher dense gas fraction and strength of fast shocks in those GMCs than in the disk, while the HCN/CO(1–0) integrated intensity ratio is moderate in the SW streamer region. The contours of those two integrated intensity ratios are extended in the directions of outflow streamers, which connect the enhanced dense gas fraction and shock strength with molecular outflow. Moreover, the molecular gas with an enhanced dense gas fraction and shock strength located at the base of the SW streamer shares the same velocity as the outflow. Conclusions. The enhanced dense gas fraction and shock strength at the base of the outflow streamers suggest that star formation inside the GMCs can trigger shocks and further drive the molecular outflow. The increased CO/13CO(1–0) integrated intensity ratio coupled with the moderate HCN/CO(1–0) integrated intensity ratio in the SW streamer region are consistent with the picture that the gas velocity gradient inside the streamer may decrease the optical depth of CO(1–0) emission, as well as the dense gas fraction in the extended streamer region.

Shockingly Bright Warm Carbon Monoxide Molecular Features in the Supernova Remnant Cassiopeia A Revealed by JWST

We present JWST NIRCam (F356W and F444W filters) and MIRI (F770W) images and NIRSpec Integral Field Unit (IFU) spectroscopy of the young Galactic supernova remnant Cassiopeia A (Cas A) to probe the physical conditions for molecular CO formation and destruction in supernova ejecta. We obtained the data as part of a JWST survey of Cas A. The NIRCam and MIRI images map the spatial distributions of synchrotron radiation, Ar-rich ejecta, and CO on both large and small scales, revealing remarkably complex structures. The CO emission is stronger at the outer layers than the Ar ejecta, which indicates the re-formation of CO molecules behind the reverse shock. NIRSpec-IFU spectra (3–5.5 μm) were obtained toward two representative knots in the NE and S fields that show very different nucleosynthesis characteristics. Both regions are dominated by the bright fundamental rovibrational band of CO in the two R and P branches, with strong [Ar vi] and relatively weaker, variable strength ejecta lines of [Si ix], [Ca iv], [Ca v], and [Mg iv]. The NIRSpec-IFU data resolve individual ejecta knots and filaments spatially and in velocity space. The fundamental CO band in the JWST spectra reveals unique shapes of CO, showing a few tens of sinusoidal patterns of rovibrational lines with pseudocontinuum underneath, which is attributed to the high-velocity widths of CO lines. Our results with LTE modeling of CO emission indicate a temperature of ∼1080 K and provide unique insight into the correlations between dust, molecules, and highly ionized ejecta in supernovae and have strong ramifications for modeling dust formation that is led by CO cooling in the early Universe.

Cross-correlation Techniques to Mitigate the Interloper Contamination for Line Intensity Mapping Experiments

Line intensity mapping (LIM) serves as a potent probe in astrophysics, relying on the statistical analysis of integrated spectral line emissions originating from distant star-forming galaxies. While LIM observations hold the promise of achieving a broad spectrum of scientific objectives, a significant hurdle for future experiments lies in distinguishing the targeted spectral line emitted at a specific redshift from undesired line emissions originating at different redshifts. The presence of these interloping lines poses a challenge to the accuracy of cosmological analyses. In this study, we introduce a novel approach to quantify line–line cross-correlations (LIM-LLX), enabling us to investigate the target signal amid instrumental noise and interloping emissions. For example, at a redshift of z ∼ 3.7, we observed that the measured auto-power spectrum of C ii 158 exhibited substantial bias, from interloping line emission. However, cross-correlating C ii 158 with CO(6–5) lines using an FYST-like experiment yielded a promising result, with a signal-to-noise ratio of ∼10. This measurement is notably unbiased. Additionally, we explore the extensive capabilities of cross-correlation by leveraging various CO transitions to probe the tomographic Universe at lower redshifts through LIM-LLX. We further demonstrate that incorporating low-frequency channels, such as 90 and 150 GHz, into FYST’s EoR-Spec-like experiment can maximize the potential for cross-correlation studies, effectively reducing the bias introduced by instrumental noise and interlopers.

Hierarchical hub-filament structures and gas inflows on galaxy-cloud scales

Abstract We investigated the kinematics and dynamics of gas structures on galaxy-cloud scales in two spiral galaxies NGC5236 (M83) and NGC4321 (M100) using CO (2–1) line. We utilized the FILFINDER algorithm on integrated intensity maps for the identification of filaments in two galaxies. Clear fluctuations in velocity and density were observed along these filaments, enabling the fitting of velocity gradients around intensity peaks. The variations in velocity gradient across different scales suggest a gradual and consistent increase in velocity gradient from large to small scales, indicative of gravitational collapse, something also revealed by the correlation between velocity dispersion and column density of gas structures. Gas structures at different scales in the galaxy may be organized into hierarchical systems through gravitational coupling. All the features of gas kinematics on galaxy-cloud scale are very similar to that on cloud-clump and clump-core scales studied in previous works. Thus, the interstellar medium from galaxy to dense core scales presents multi-scale/hierarchical hub-filament structures. Like dense core as the hub in clump, clump as the hub in molecular cloud, now we verify that cloud or cloud complex can be the hub in spiral galaxies. Although the scaling relations and the measured velocity gradients support the gravitational collapse of gas structures on galaxy-cloud scales, the collapse is much slower than a pure free-fall gravitational collapse.