Structure Of Dense Gas Research Articles

Understanding the star formation efficiency in dense gas: Initial results from the CAFFEINE survey with ArTéMiS

Context. Despite recent progress, the question of what regulates the star formation efficiency (SFE) in galaxies remains one of the most debated problems in astrophysics. According to the dominant picture, star formation (SF) is regulated by turbulence and feedback, and the SFE is ~1–2% or less per local free-fall time on all scales from Galactic clouds to high-redshift galaxies. In an alternate scenario, the star formation rate (SFR) in galactic disks is linearly proportional to the mass of dense gas above some critical density threshold ~104 cm–3. Aims. We aim to discriminate between these two pictures thanks to high-resolution submillimeter and mid-infrared imaging observations, which trace both dense gas and young stellar objects (YSOs) for a comprehensive sample of 49 nearby massive SF complexes out to a distance of d ~ 3 kpc in the Galactic disk. Methods. We used data from CAFFEINE, a complete 350/450 µm survey with APEX/ArTéMiS of the densest portions of all southern molecular clouds at d ≲ 3 kpc, in combination with Herschel data to produce column density maps at a factor of ~4 higher resolution (8") than standard Herschel column density maps (36″). Our maps are free of any saturation effect around luminous high-mass pro-tostellar objects and resolve the structure of dense gas and the typical ~0.1 pc width of molecular filaments out to 3 kpc, which is the most important asset of the present study and is impossible to achieve with Herschel data alone. Coupled with SFR estimates derived from Spitzer mid-infrared observations of the YSO content of the same clouds, this allowed us to study the dependence of the SFE on density in the CAFFEINE clouds. We also combine our findings with existing SF efficiency measurements in nearby clouds to extend our analysis down to lower column densities. Results. Our results suggest that the SFE does not increase with density above the critical threshold and support a scenario in which the SFE in dense gas is approximately constant (independent of free-fall time). However, the SF efficiency measurements traced by Class I YSOs in nearby clouds are more inconclusive, since they are consistent with both the presence of a density threshold and a dependence on density above the threshold. Overall, we suggest that the SF efficiency in dense gas is primarily governed by the physics of filament fragmentation into protostellar cores.

Gas inflows from cloud to core scales in G332.83-0.55: Hierarchical hub-filament structures and tide-regulated gravitational collapse

The massive star-forming region G332.83-0.55 contains at least two levels of hub-filament structures. The hub-filament structures may form through the “gravitational focusing” process. High-resolution LAsMA and ALMA observations can directly trace the gas inflows from cloud to core scales. We investigated the effects of shear and tides from the protocluster on the surrounding local dense gas structures. Our results seem to deny the importance of shear and tides from the protocluster. However, for a gas structure, it bears the tidal interactions from all external material, not only the protocluster. To fully consider the tidal interactions, we derived the tide field according to the surface density distribution. Then, we used the average strength of the external tidal field of a structure to measure the total tidal interactions that are exerted on it. For comparison, we also adopted an original pixel-by-pixel computation to estimate the average tidal strength for each structure. Both methods give comparable results. After considering the total tidal interactions, for the scaling relation between the velocity dispersion σ, the effective radius R, and the column density N of all the structures, the slope of the σ − N * R relation changes from 0.20 ± 0.04 to 0.52 ± 0.03, close to 0.5 of the pure free-fall gravitational collapse, and the correlation also becomes stronger. Thus, the deformation due to the external tides can effectively slow down the pure free-fall gravitational collapse of gas structures. The external tide tries to tear up the structure, but the external pressure on the structure prevents this process. The counterbalance between the external tide and external pressure hinders the free-fall gravitational collapse of the structure, which can also cause the pure free-fall gravitational collapse to be slowed down. These mechanisms can be called “tide-regulated gravitational collapse”.

Feedback from protoclusters does not significantly change the kinematic properties of the embedded dense gas structures

A total of 64 ATOMS sources at different evolutionary stages were selected to investigate the kinematics and dynamics of gas structures under feedback. We identified dense gas structures based on the integrated intensity map of H13CO+ J = 1−0 emission, and then extracted the average spectra of all the structures to investigate their velocity components and gas kinematics. For the scaling relations between the velocity dispersion, σ, the effective radius, R, and the column density, N, of all the structures, σ − N * R always has a stronger correlation compared to σ − N and σ − R. There are significant correlations between velocity dispersion and column density, which may imply that the velocity dispersion originates in gravitational collapse, also revealed by the velocity gradients. The measured velocity gradients for dense gas structures in early-stage sources and late-stage sources are comparable, indicating gravitational collapse through all evolutionary stages. Late-stage sources do not have large-scale hub-filament structures, but the embedded dense gas structures in late-stage sources show similar kinematic modes to those in early- and middle-stage sources. These results may be explained by the multi-scale hub-filament structures in the clouds. We quantitatively estimated the velocity dispersion generated by the outflows, inflows, ionized gas pressure, and radiation pressure, and found that the ionized gas feedback is stronger than other feedback mechanisms. However, although feedback from HII regions is the strongest, it does not significantly affect the physical properties of the embedded dense gas structures. Combined with the conclusions in our previous work on cloud-clump scales, we suggest that although feedback from cloud to core scales will break up the original cloud complex, the substructures of the original complex can be reorganized into new gravitationally governed configurations around new gravitational centers. This process is accompanied by structural destruction and generation, and changes in gravitational centers, but gravitational collapse is always ongoing.

High-resolution APEX/LAsMA 12CO and 13CO (3–2) observation of the G333 giant molecular cloud complex

Context. Feedback from young massive stars has an important impact on the star formation potential of their parental molecular clouds. Aims. We investigate the physical properties of gas structures under feedback in the G333 complex using data of the 13CO J = 3–2 line observed with the LAsMA heterodyne camera on the APEX telescope. Methods. We used the Dendrogram algorithm to identify molecular gas structures based on the integrated intensity map of the 13CO (3–2) emission, and extracted the average spectra of all structures to investigate their velocity components and gas kinematics. Results. We derive the column density ratios between different transitions of the 13CO emission pixel by pixel, and find the peak values N2−1/N1−0 ≈ 0.5, N3−2/N1−0 ≈ 0.3, and N3−2/N2−1 ≈ 0.5. These ratios can also be roughly predicted by the nonlocal thermodynamic equilibrium (NLTE) molecular radiative transfer code RADEX for an average H2 volume density of ~4.2 × 103 cm−3. A classical virial analysis does not reflect the true physical state of the identified structures, and we find that external pressure from the ambient cloud plays an important role in confining the observed gas structures. For high-column-density structures, velocity dispersion and density show a clear correlation that is not seen for low-column-density structures, indicating the contribution of gravitational collapse to the velocity dispersion. Branch structures show a more significant correlation between 8 μm surface brightness and velocity dispersion than leaf structures, implying that feedback has a greater impact on large-scale structures. For both leaf and branch structures, σ − N * R always has a stronger correlation compared to σ − N and σ − R. The scaling relations are stronger, and have steeper slopes when considering only self-gravitating structures, which are the structures most closely associated with the Heyer relation. Conclusions. Although the feedback disrupting the molecular clouds will break up the original cloud complex, the substructures of the original complex can be reorganized into new gravitationally governed configurations around new gravitational centers. This process is accompanied by structural destruction and generation, and changes in gravitational centers, but gravitational collapse is always ongoing.

Dense gas properties around the centre of the Circinus galaxy

ABSTRACT Studies of spatial distribution, kinematics, and physical structure of dense molecular gas in the centres of galaxies is crucial in understanding the nature of physical processes, especially those related to active galactic nuclei (AGNs) and starbursts, in the vicinity of central massive black holes. However, our current understanding of dense gas in the circumnuclear discs (CND) remains limited, particularly those around AGNs. Here-, we report high-resolution deep observations of the CS (3−2), HC3N (15−14), H2CO (2−1), CH3CCN (8−7), CO (1−0), and (3−2) lines, toward the central region (R ∼ 0.4 kpc) of the nearest Seyfert-2 galaxy, the Circinus galaxy, using the Atacama Large Millimeter/submillimeter Array. The dense molecular gas traced by CS (3−2) and HC3N (15−14) exhibit an asymmetric and cloud structure, which is in sharp contrast against the symmetric low-density structure traced by CO (1−0). Four molecular cloud structures surrounding the AGN, which have a typical size of ∼ 20 pc, are detected with the optically thin HC3N (15−14) line. HC3N (15−14)/CO (1−0) line ratio in the CND is found to be higher than that in the nuclear region, indicating that the bulk of dense gas accumulates in the CND, instead of the nucleus. HC3N (15−14) line and optically thick tracers in four molecular clouds, provide a good comparison between density layers and show evidence for a density stratification of the CND.

Fragmentation and dynamics of dense gas structures in the proximity of massive young stellar object W42-MME

ABSTRACT We present an analysis of the dense gas structures in the immediate surroundings of the massive young stellar object (MYSO) W42-MME, using the high-resolution (0″.31 × 0″.25) Atacama Large Millimetre/submillimetre Array dust continuum and molecular line data. We performed a dendrogram analysis of H13CO+ (4–3) line data to study multiscale structures and their spatio–kinematic properties, and analysed the fragmentation and dynamics of dense structures down to ∼2000 au scale. Our results reveal 19 dense gas structures, out of which 12 are leaves and 7 are branches in dendrogram terminology. These structures exhibit transonic–supersonic gas motions (1$\lt \mathcal {M}\lt 5$) with overvirial states (αvir ≥ 2). The non-thermal velocity dispersion–size relation (σnt–L) of dendrogram structures shows a weak negative correlation, while the velocity dispersion across the sky ($\delta \mathit {V_{\rm lsr}}$) correlates positively with structure size (L). Velocity structure function (S2(l)1/2) analysis of H13CO+ data reveals strong power-law dependences with lag (l) up to a scale length of ≲6000 au. The mass–size (M–R) relation of dendrogram structures shows a positive correlation with power-law index of 1.73 ± 0.23, and the leaf L17 hosting W42-MME meets the mass–size conditions for massive star formation. Blue asymmetry is observed in the H12CO+ (4–3) line profiles of most of the leaves, indicating infall. Overall, our results observationally support the hierarchical and chaotic collapse scenario in the proximity of the MYSO W42-MME.

The role of magnetic fields in the stability and fragmentation of filamentary molecular clouds: two case studies at OMC-3 and OMC-4

ABSTRACT We present the stability analysis of two regions, OMC-3 and OMC-4, in the massive and long molecular cloud complex of Orion A. We obtained 214 $\mu$m HAWC + /SOFIA polarization data, and we make use of archival data for the column density and C18O (1–0) emission line. We find clear depolarization in both observed regions and that the polarization fraction is anticorrelated with the column density and the polarization-angle dispersion function. We find that the filamentary cloud and dense clumps in OMC-3 are magnetically supercritical and strongly subvirial. This region should be in the gravitational collapse phase and is consistent with many young stellar objects (YSOs) forming in the region. Our histogram of relative orientation (HRO) analysis shows that the magnetic field is dynamically sub-dominant in the dense gas structures of OMC-3. We present the first polarization map of OMC-4. We find that the observed region is generally magnetically subcritical except for an elongated dense core, which could be a result of projection effect of a filamentary structure aligned close to the line of sight. The relative large velocity dispersion and the unusual positive shape parameters at high column densities in the HROs analysis suggest that our viewing angle may be close to axes of filamentary substructures in OMC-4. The dominating strong magnetic field in OMC-4 is unfavourable for star formation and is consistent with much fewer YSOs than in OMC-3.

Chemically Fresh Gas Inflows Detected in a Nearby High-mass Star-forming Region

Abstract We report the detection of a chemically fresh inflow that is feeding high-mass young-stellar-object (HMYSO) growth in the nearby high-mass star-forming region G352.63 made with both the Atacama Large Millimeter/submillimeter Array (ALMA) and the Submillimeter Array (SMA). High-quality images of the dust and molecular lines from both ALMA and SMA have consistently revealed a gravitationally controlled cold (∼10 K) gas inflow of chemically fresh molecules (e.g., CCH and HC3N) toward the central HMYSO and its surrounding dense gas structure, which has a possible torus- or disk-like morphology. The HMYSO is also observed to have an outflow, which is nearly perpendicular to the torus and its parental filament, and thus can be clearly separated from the inflows. These kinematic features provide observational evidence to support the conjecture that the infalling streamers in high-mass star-forming regions could proceed in a similar process to that observed in low-mass counterparts. The chemically fresh infalling streamers could also be involved in the disk or torus configuration, fragmentation, and accretion bursts that occur in both simulations and observations.

Probing the structure of a massive filament: ArTéMiS 350 and 450 μm mapping of the integral-shaped filament in Orion A

Context. The Orion molecular cloud is the closest region of high-mass star formation. It is an ideal target for investigating the detailed structure of massive star-forming filaments at high resolution and the relevance of the filament paradigm for the earliest stages of intermediate- to high-mass star formation. Aims. Within the Orion A molecular cloud, the integral-shaped filament (ISF) is a prominent, degree-long structure of dense gas and dust with clear signs of recent and ongoing high-mass star formation. Our aim is to characterise the structure of this massive filament at moderately high angular resolution (8′′ or ~0.016 pc) in order to measure the intrinsic width of the main filament, down to scales well below 0.1 pc, which has been identified as the characteristic width of filaments. Methods. We used the ArTéMiS bolometer camera at APEX to map a ~0.6 × 0.2 deg2 region covering OMC-1, OMC-2, and OMC-3 at 350 and 450 μm. We combined these data with Herschel-SPIRE maps to recover extended emission. The combined Herschel-ArTéMiS maps provide details on the distribution of dense cold material, with a high spatial dynamic range, from our 8′′ resolution up to the transverse angular size of the map, ~10–15′. By combining Herschel and ArTéMiS data at 160, 250, 350, and 450 μm, we constructed high-resolution temperature and H2 column density maps. We extracted radial intensity profiles from the column density map in several representative portions of the ISF, which we fitted with Gaussian and Plummer models to derive their intrinsic widths. We also compared the distribution of material traced by ArTéMiS with that seen in the higher-density tracer N2H+(1–0) that was recently observed with the ALMA interferometer. Results. All the radial profiles that we extracted show a clear deviation from a Gaussian, with evidence for an inner plateau that had not previously been seen clearly using Herschel-only data. We measure intrinsic half-power widths in the range 0.06–0.11 pc. This is significantly larger than the Gaussian widths measured for fibres seen in N2H+, which probably only traces the dense innermost regions of the large-scale filament. These half-power widths are within a factor of two of the value of ~0.1 pc found for a large sample of nearby filaments in various low-mass star-forming regions, which tends to indicate that the physical conditions governing the fragmentation of pre-stellar cores within transcritical or supercritical filaments are the same over a large range of masses per unit length.

Kinematic study of the molecular gas associated with two cometary globules in Sh2−236

Aims. Cometary globules, dense molecular gas structures exposed to UV radiation, are found inside H II regions. Understanding the nature and origin of these structures through a kinematic study of the molecular gas could be useful to advance in our knowledge of the interplay between radiation and molecular gas. Methods. Using the Atacama Submillimeter Telescope Experiment (Chile), we carried out molecular observations toward two cometary globules (Sim129 and Sim130) in the H II region Sh2−236. We mapped two regions of about 1′ × 1′ with the 12CO J = 3−2 and HCO+ J = 4−3 lines. Additionally, we carried out two single pointings with the C2H N = 4–3, HNC, and HCN J = 4−3 transitions. The angular resolution was about 22′′. We combined our molecular observations with public infrared and optical data to analyze the distribution and kinematics of the molecular gas. Results. We find kinematic signatures of infalling gas in the 12CO J = 3−2 and C2H N = 4−3 spectra toward Sim 129. We detect HCO+, HCN, and HNC J = 4−3 only toward Sim 130. The HCN/HNC integrated ratio of about three found in Sim 130 suggests that the possible star-formation activity inside this globule has not yet ionized the gas. The location of the NVSS source 052255+33315, which peaks toward the brightest border of the globule, supports this scenario. The non-detection of these molecules toward Sim 129 could be due to the radiation field arising from the star-formation activity inside this globule. The ubiquitous presence of the C2H molecule toward Sim 129 and Sim 130 evidences the action of the nearby O-B stars irradiating the external layer of both globules. Based on the mid-infrared 5.8 μm emission, we identify two new structures: (1) a region of diffuse emission (R1) located, in projection, in front of the head of Sim 129 and (2) a pillar-like feature (P1) placed besides Sim 130. Based on the 12CO J = 3−2 transition, we find molecular gas associated with Sim 129, Sim 130, R1, and P1 at radial velocities of −1.5, −11, +10, and +4 km s−1, respectively. Therefore, while Sim 129 and P1 are located at the far side of the shell, Sim 130 is placed at the near side, consistent with earlier results. Finally, the molecular gas related to R1 exhibits a radial velocity that differs in more than 11 km s−1 with the radial velocity of S129, which suggests that while S129 is located at the far side of the expanding shell, R1 would be placed well beyond.

Relative Alignment between the Magnetic Field and Molecular Gas Structure in the Vela C Giant Molecular Cloud Using Low- and High-density Tracers

We compare the magnetic field orientation for the young giant molecular cloud Vela C inferred from 500-$\mu$m polarization maps made with the BLASTPol balloon-borne polarimeter to the orientation of structures in the integrated line emission maps from Mopra observations. Averaging over the entire cloud we find that elongated structures in integrated line-intensity, or zeroth-moment maps, for low density tracers such as $^{12}$CO and $^{13}$CO $J$ $\rightarrow$ 1 - 0 are statistically more likely to align parallel to the magnetic field, while intermediate or high density tracers show (on average) a tendency for alignment perpendicular to the magnetic field. This observation agrees with previous studies of the change in relative orientation with column density in Vela C, and supports a model where the magnetic field is strong enough to have influenced the formation of dense gas structures within Vela C. The transition from parallel to no preferred/perpendicular orientation appears to happen between the densities traced by $^{13}$CO and by C$^{18}$O $J$ $\rightarrow$ 1 - 0. Using RADEX radiative transfer models to estimate the characteristic number density traced by each molecular line we find that the transition occurs at a molecular hydrogen number density of approximately $10^3$ cm$^{-3}$. We also see that the Centre-Ridge (the highest column density and most active star-forming region within Vela C) appears to have a transition at a lower number density, suggesting that this may depend on the evolutionary state of the cloud.

Fire in the Heart: A Characterization of the High Kinetic Temperatures and Heating Sources in the Nucleus of NGC 253

Abstract The nuclear starburst within the central ∼15″ (∼250 pc; 1″ ≃ 17 pc) of NGC 253 has been extensively studied as a prototype for the starburst phase in galactic evolution. Atacama Large Millimeter/submillimeter Array (ALMA) imaging within receiver Bands 6 and 7 has been used to investigate the dense gas structure, kinetic temperature, and heating processes that drive the NGC 253 starburst. A total of 29 transitions from 15 molecular species/isotopologues have been identified and imaged at 1.″5–0.″4 resolution, allowing for the identification of five of the previously studied giant molecular clouds within the central molecular zone (CMZ) of NGC 253. Ten transitions from the formaldehyde (H2CO) molecule have been used to derive the kinetic temperature within the ∼0.″5–5″ dense gas structures imaged. On ∼5″ scales we measure T K ≳ 50 K, while on size scales ≲1″ we measure T K ≳ 300 K. These kinetic temperature measurements further delineate the association between potential sources of dense gas heating. We have investigated potential heating sources by comparing our measurements to models that predict the physical conditions associated with dense molecular clouds that possess a variety of heating mechanisms. This comparison has been supplemented with tracers of recently formed massive stars (Brγ) and shocks ([Fe ii]). Derived molecular column densities point to a radially decreasing abundance of molecules with sensitivity to cosmic-ray and mechanical heating within the NGC 253 CMZ. These measurements are consistent with radio spectral index calculations that suggest a higher concentration of cosmic-ray-producing supernova remnants within the central 10 pc of NGC 253.

Investigating Fragmentation of Gas Structures in OB Cluster-forming Molecular Clump G33.92+0.11 with 1000 au Resolution Observations of ALMA

Abstract We report new ∼1000 au spatial resolution observations of 225 GHz dust continuum emission toward the OB cluster-forming molecular clump G33.92+0.11. On parsec scales, this molecular clump presents a morphology with several arm-like dense gas structures surrounding the two central massive (≳100 M ⊙) cores. From the new higher resolution observations, we identified 28 localized, spatially compact dust continuum emission sources, which may be candidates of young stellar objects (YSOs). Only one of them is not embedded within known arm-like (or elongated) dense gas structures. The spatial separations of these compact sources can be very well explained by Jeans lengths. We found that G33.92+0.11 may be consistently described by a marginally centrifugally supported, Toomre unstable accretion flow that is approximately in a face-on projection. The arm-like overdensities are a natural consequence of the Toomre instability, which can fragment to form YSOs in shorter timescales than the timescale of the global clump contraction. On our resolved spatial scales, there is no evidence yet that the fragmentation is halted by turbulence, the magnetic field, or stellar feedback.

The impact of radiation feedback on the assembly of star clusters in a galactic context

Massive star clusters are observed in a broad range of galaxy luminosity and types, and are assumed to form in dense gas-rich environments. Using a parsec-resolution hydrodynamical simulation of an isolated gas-rich low mass galaxy, we discuss here the non-linear effects of stellar feedback on the properties of star clusters with a focus on the progenitors of nuclear clusters. Our simulation shows two categories of star clusters: those for which feedback expels gas leftovers associated with their formation sites, and those, in a denser environment around which feedback fails at totally clearing the gas. We confirm that radiation feedback (photo-ionization and radiative pressure) plays a more important role than type-II supernovae in destroying dense gas structures, and altering or quenching the subsequent cluster formation. It also disturbs the cluster mass growth, by increasing the internal energy of the gas component to the point when radiation pressure overcomes the cluster gravity. We discuss how these effects may depend on the local properties of the interstellar medium, and also on the details of the subgrid recipes, which can affect the available cluster gas reservoirs, the evolution of potential nuclear clusters progenitors, and the overall galaxy morphology.

Core Emergence in a Massive Infrared Dark Cloud: A Comparison between Mid-IR Extinction and 1.3 mm Emission

Abstract Stars are born from dense cores in molecular clouds. Observationally, it is crucial to capture the formation of cores in order to understand the necessary conditions and rate of the star formation process. The Atacama Large Millimeter/submillimeter Array (ALMA) is extremely powerful for identifying dense gas structures, including cores, at millimeter wavelengths via their dust continuum emission. Here, we use ALMA to carry out a survey of dense gas and cores in the central region of the massive (∼105 M ⊙) infrared dark cloud (IRDC) G28.37+0.07. The observation consists of a mosaic of 86 pointings of the 12 m array and produces an unprecedented view of the densest structures of this IRDC. In this first Letter about this data set, we focus on a comparison between the 1.3 mm continuum emission and a mid-infrared (MIR) extinction map of the IRDC. This allows estimation of the “dense gas” detection probability function (DPF), i.e., as a function of the local mass surface density, Σ, for various choices of thresholds of millimeter continuum emission to define “dense gas.” We then estimate the dense gas mass fraction, f dg, in the central region of the IRDC and, via extrapolation with the DPF and the known Σ probability distribution function, to the larger-scale surrounding regions, finding values of about 5% to 15% for the fiducial choice of threshold. We argue that this observed dense gas is a good tracer of the protostellar core population and, in this context, estimate a star formation efficiency per free-fall time in the central IRDC region of ϵ ff ∼ 10%, with approximately a factor of two systematic uncertainties.

Giant molecular cloud collisions as triggers of star formation. VI. Collision-induced turbulence

Abstract We investigate collisions between giant molecular clouds (GMCs) as potential generators of their internal turbulence. Using magnetohydrodynamic (MHD) simulations of self-gravitating, magnetized, turbulent GMCs, we compare kinematic and dynamic properties of dense gas structures formed when such clouds collide compared to those that form in non-colliding clouds as self-gravity overwhelms decaying turbulence. We explore the nature of turbulence in these structures via distribution functions of density, velocity dispersions, virial parameters, and momentum injection. We find that the dense clumps formed from GMC collisions have higher effective Mach number, greater overall velocity dispersions, sustain near-virial equilibrium states for longer times, and are the conduit for the injection of turbulent momentum into high density gas at high rates.

The Green Bank Ammonia Survey: Observations of Hierarchical Dense Gas Structures in Cepheus-L1251

Abstract We use Green Bank Ammonia Survey observations of NH3 (1, 1) and (2, 2) emission with 32″ FWHM resolution from a ∼10 pc2 portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH3 data results in 22 top-level structures, which reside within 13 lower-level parent structures. The structures are compact and are spatially correlated with the highest H2 column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.″2 FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH3 column density, derived from detailed modeling of the NH3 data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median T dust and T K measurements of 11.7 ± 1.1 K and 10.3 ± 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS (20 − 10) and HC5N emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH3 (1, 1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores.

Low-redshift Lyman limit systems as diagnostics of cosmological inflows and outflows

We use cosmological hydrodynamic simulations with stellar feedback from the FIRE project to study the physical nature of Lyman limit systems (LLSs) at z<1. At these low redshifts, LLSs are closely associated with dense gas structures surrounding galaxies, such as galactic winds, dwarf satellites, and cool inflows from the intergalactic medium. Our analysis is based on 14 zoom-in simulations covering the halo mass range M_h~10^9-10^13 Msun at z=0, which we convolve with the dark matter halo mass function to produce cosmological statistics. We find that the majority of cosmologically-selected LLSs are associated with halos in the mass range 10^10 < M_h < 10^12 Msun. The incidence and HI column density distribution of simulated absorbers with columns 10^16.2 < N_HI < 2x10^20 cm^-2 are consistent with observations. High-velocity outflows (with radial velocity exceeding the halo circular velocity by a factor >~2) tend to have higher metallicities ([X/H] ~ -0.5) while very low metallicity ([X/H] < -2) LLSs are typically associated with gas infalling from the intergalactic medium. However, most LLSs occupy an intermediate region in metallicity-radial velocity space, for which there is no clear trend between metallicity and radial kinematics. Metal-enriched inflows arise in the FIRE simulations as a result of galactic winds that fall back onto galaxies at low redshift. The overall simulated LLS metallicity distribution has a mean (standard deviation) [X/H] = -0.9 (0.4) and does not show significant evidence for bimodality, in contrast to recent observational studies but consistent with LLSs arising from halos with a broad range of masses and metallicities.

Interpreting observations of edge-on gravitationally unstable accretion flows

Gravitational collapse of molecular cloud or cloud core/clump may lead to the formation of geometrically flattened, rotating accretion flow surrounding the new born star or star cluster. Gravitational instability may occur in such accretion flow when the gas to stellar mass ratio is high (e.g. over $\sim$10\%). This paper takes the OB cluster-forming region G10.6-0.4 as an example. We introduce the enclosed gas mass around its central ultra compact (UC) H\textsc{ii} region, addresses the gravitational stability of the accreting gas, and outline the observed potential signatures of gravitational instability. The position-velocity (PV) diagrams of various molecular gas tracers on G10.6-0.4 consistently show asymmetry in the spatial and the velocity domain. We deduce the morphology of the dense gas accretion flow by modeling velocity distribution of the azimuthally asymmetric gas structures, and by directly de-projecting the PV diagrams. We found that within the 0.3 pc radius, an infall velocity of 1-2 km\,s$^{-1}$ may be required to explain the observed PV diagrams. In addition, the velocity distribution traced in the PV diagrams can be interpreted by spiral arm-like structures, which may be connected with exterior infalling gas filaments. The morphology of dense gas structures we propose appears very similar to the spatially resolved gas structures around the OB cluster-forming region G33.92+0.11 with similar gas mass and size, which however is likely to be approximately in a face-on projection. The dense gas accretion flow around G10.6-0.4 appears to be Toomre unstable, which is consistent with the existence of large-scale spiral arm-like structures, and the formation of localize gas condensations.

Densitometry and Thermometry of Starburst Galaxies

With a goal toward deriving the physical conditions in external galaxies, we have conducted a survey and subsequent high spatial resolution imaging of formaldehyde (H2CO) and ammonia (NH3) emission and absorption in a sample of starburst galaxies. In this article we present the results from a subset of this survey which focuses on high spatial resolution measurements of volume density-and kinetic temperature-sensitive transitions of the H2CO molecule. The volume density structure toward the nuclear region of NGC 253 has been derived from ? ≠4 arcsec NRAO Very Large Array (VLA) measurements of the 110-111 and 211-212 K-doublet transitions of H2CO. The kinetic temperature structure toward NGC 253 and NGC 4945 has been derived from ? ≠0.5-1.0 arcsec measurements of the H2CO 3K-1K+1-2K-1K+1 (near 218 GHz) and 5K-1K+1-4K-1K+1 (near 365 GHz) transitions acquired using the Atacama Large Millimeter/submillimeter Array (ALMA). These measurements have allowed us to characterize the dense gas and kinetic temperature structure within these star forming galaxies, which is a first step toward associating dense star-forming gas and the heating processes at work within galaxies.