Molecular Interstellar Medium Research Articles

Graph-based interpretation of the molecular interstellar medium segmentation

We present a generalization of the Giant Molecular Cloud (GMC) identification problem based on cluster analysis. The method we designed, SCIMES (Spectral Clustering for Interstellar Molecular Emission Segmentation) considers the dendrogram of emission in the broader framework of graph theory and utilizes spectral clustering to find discrete regions with similar emission properties. For Galactic molecular cloud structures, we show that the characteristic volume and/or integrated CO luminosity are useful criteria to define the clustering, yielding emission structures that closely reproduce "by-eye" identification results. SCIMES performs best on well-resolved, high-resolution data, making it complementary to other available algorithms. Using 12CO(1-0) data for the Orion-Monoceros complex, we demonstrate that SCIMES provides robust results against changes of the dendrogram-construction parameters, noise realizations and degraded resolution. By comparing SCIMES with other cloud decomposition approaches, we show that our method is able to identify all canonical clouds of the Orion-Monoceros region, avoiding the over-division within high resolution survey data that represents a common limitation of several decomposition algorithms. The Orion-Monoceros objects exhibit hierarchies and size-line width relationships typical to the turbulent gas in molecular clouds, although "the Scissors" region deviates from this common description. SCIMES represents a significant step forward in moving away from pixel-based cloud segmentation towards a more physical-oriented approach, where virtually all properties of the ISM can be used for the segmentation of discrete objects.

SWAN: NGC 253's Nucleated Star Bursting Environment

AbstractWe present the first results from SWAN: “Survey of Water and Ammonia in Nearby galaxies”. Nearby galaxies are conveniently located to probe molecular gas properties on scales of 10 to 200 pc, which are appropriate for the study of Giant Molecular Clouds (GMCs). The resolution of the Very Large Array in D and C configurations corresponds to a few 10s of parsecs in these galaxies. To advance studies of galaxy evolution it is paramount to understand how processes in the molecular Interstellar Medium(ISM) and star formation are linked on these scales. We have observed the metastable transitions of ammonia and the 22GHz water maser line in four nearby galaxies: NGC 253, IC 342, NGC 2146, and NGC 6946 using the VLA. These galaxies were chosen to span an order of magnitude in star formation rate, and a range of galactic ecosystems. We use the ammonia transitions to derive kinetic temperatures, which exposes the heating and cooling balance of the ISM. We then aim to relate these conditions to energetic feedback from star formation as indicated by water masers.Currently, our analysis is focused on NGC 253. NGC 253 is a barred spiral starburst galaxy with a nucleated star formation rate (SFR) of ≈ 3M⊙ per year. We use a distance of 3.9Mpc for analysis. We have observed ammonia transitions (1,1) to (5,5) and the 22GHz water maser line with a resolution of ≈63pc. We have identified nine regions across the nucleated starburst for study. The ammonia (3,3) line appears to be masing in the centermost 200pc. We have identified two regions of water maser emission. The first region is a minor axis extension, about the center of the galaxy, and is very close to the outflow. By means of these measurements, we gain an understanding of the molecular ISM associated with the nucleated starburst environment in NGC 253.

Atomic and Molecular Phases of the Interstellar Medium

AbstractThis review covers four current questions in the behavior of the atomic and molecular interstellar medium. These include whether the atomic gas originates primarily in cold streams or hot flows onto galaxies; what the filling factor of cold gas actually is in galactic regions observationally determined to be completely molecular; whether molecular hydrogen determines or merely traces star formation; and whether gravity or turbulence drives the dynamical motions observed in interstellar clouds, with implications on their star formation properties.

The linewidth-size scaling law of molecular gas revisited

AbstractThe origin of the linewidth-size (LWS) scaling law, first noticed by Larson three decades ago and ascribed to turbulence, is still a highly debated issue. Not unexpectedly, its properties depend on the environment and on the line tracer used. When the optically thick 12CO (J=1-0) line is used, a specific medium is sampled: the translucent molecular gas of moderate density that builds up the bulk of the molecular interstellar medium in galaxies like the Milky Way. The sensitivity of the 12CO line to this gas is such that the LWS is found to hold over almost five orders of magnitude in lengthscale, although with a considerable scatter (± 0.5 dex). It also appears to split into two regimes, depending on the gas mass surface density: below a given threshold that is proposed to be linked to the galactic structure, it bears the signature of a turbulent cascade, while above it, the scaling law is ascribed to virial balance. Large deviations from the LWS scaling law are observed at small scales where signatures of turbulent intermittency appear. The mass-size scaling law built with the 12CO (J=1-0) line also splits into two regimes. The mass surface density is uniform (also with a large scatter) above lengthscales ~ 10pc and increases with size at smaller scales, following turbulence predictions. The two thresholds define an average gas density nH ~ 300 cm−3.

A SURVEY OF THE MOLECULAR ISM PROPERTIES OF NEARBY GALAXIES USING THEHERSCHELFTS

The 12CO J=4-3 to J=13-12 lines of the interstellar medium from nearby galaxies, newly observable with the Herschel SPIRE Fourier Transform Spectrometer (FTS), offer an opportunity to study warmer, more luminous molecular gas than that traced by 12CO J=1-0. Here we present a survey of 17 nearby infrared-luminous galaxy systems (21 pointings). In addition to photometric modeling of dust, we modeled full 12CO spectral line energy distributions from J=1-0 to J=13-12 with two components of warm and cool CO gas, and included LTE analysis of [CI], [CII], [NII] and H2 lines. CO is emitted from a low-pressure/high-mass component traced by the low-J lines and a high-pressure/low-mass component which dominates the luminosity. We found that, on average, the ratios of the warm/cool pressure, mass, and 12CO luminosity are 60 +/- 30, 0.11 +/- 0.02, and 15.6 +/- 2.7. The gas-to-dust-mass ratios are < 120 throughout the sample. The 12CO luminosity is dominated by the high-J lines and is 4 $\times 10^{-4}$ LFIR on average. We discuss systematic effects of single-component and multi-component CO modeling (e.g., single-component J < 3 models overestimate gas pressure by ~ 0.5 dex), as well as compare to Galactic star-forming regions. With this comparison, we show the molecular interstellar medium of starburst galaxies is not simply an ensemble of Galactic-type GMCs. The warm gas emission is likely dominated by regions resembling the warm extended cloud of Sgr B2.

Multimolecular studies of Galactic star-forming regions

Molecular emission-line observations of isolated Galactic star-forming regions are used to model the physical properties of the molecular interstellar medium in these systems. Observed line ratios ...

The cycling of carbon into and out of dust.

Observational evidence seems to indicate that the depletion of interstellar carbon into dust shows rather wide variations and that carbon undergoes rather rapid recycling in the interstellar medium (ISM). Small hydrocarbon grains are processed in photo-dissociation regions by UV photons, by ion and electron collisions in interstellar shock waves and by cosmic rays. A significant fraction of hydrocarbon dust must therefore be re-formed by accretion in the dense, molecular ISM. A new dust model (Jones et al., Astron. Astrophys., 2013, 558, A62) shows that variations in the dust observables in the diffuse interstellar medium (n(H) < or = 10(3) cm(-3)), can be explained by systematic and environmentally-driven changes in the small hydrocarbon grain population. Here we explore the consequences of gas-phase carbon accretion onto the surfaces of grains in the transition regions between the diffuse ISM and molecular clouds (e.g., Jones, Astron. Astrophys., 2013, 555, A39). We find that significant carbonaceous dust re-processing and/or mantle accretion can occur in the outer regions of molecular clouds and that this dust will have significantly different optical properties from the dust in the adjacent diffuse ISM. We conclude that the (re-)processing and cycling of carbon into and out of dust is perhaps the key to advancing our understanding of dust evolution in the ISM.

ISM chemistry in metal-rich environments: molecular tracers of metallicity

In this paper we use observations of molecular tracers in metal rich and alpha-enhanced galaxies to study the effect of abundance changes on molecular chemistry. We selected a sample of metal rich spiral and star bursting objects from the literature, and present here new data for a sample of early-type galaxies (ETGs). We conducted the first survey of CS and methanol emission in ETGs, detecting 7 objects in CS, and 5 in methanol emission. We find evidence to support the hypothesis that CS is a better tracer of dense star-forming gas than HCN. We suggest that the methanol emission in these sources is driven by dust mantle destruction due to ionisation from high mass star formation, but cannot rule out shocks dominating in some sources. The derived source averaged CS/methanol column densities and rotation temperatures are similar to those found in normal spiral and starburst galaxies, suggesting dense clouds are little affected by the differences between galaxy types. Finally we used the total column density ratios for our galaxy samples to show for the first time that some molecular tracers do seem to show systematic variations that appear to correlate with metallicity, and that these variations roughly match those predicted by chemical models. Using this fact, the chemical models of Bayet et al. (2012b), and assumptions about the optical depth we are able to roughly predict the metallicity of our spiral and ETG sample, with a scatter of ~0.3 dex. We provide the community with linear approximations to the relationship between the HCN and CS column density ratio and metallicity. Further study will clearly be required to determine if this, or any, molecular tracer can be used to robustly determine gas-phase metallically, but that a relationship exists at all suggests that in the future it may be possible to calibrate a metallicity indicator for the molecular interstellar medium (abridged).

Principal component analysis of spectral line data: analytic formulation

Principal component analysis is a powerful statistical system to investigate the structure and dynamics of the molecular interstellar medium, with particular emphasis on the study of turbulence, as revealed by spectroscopic imaging of molecular line emission. To date, the method to retrieve the power-law index of the velocity structure function or power spectrum has relied on an empirical calibration and testing with model turbulent velocity fields, while lacking a firm theoretical basis. In this paper, we present an analytic formulation that reveals the detailed mechanics of the method and confirms previous empirical calibrations of its recovery of the scale dependence of turbulent velocity fluctuations.

The Supershell–Molecular Cloud Connection: Large-Scale Stellar Feedback and the Formation of the Molecular ISM

AbstractThe accumulation, compression, and cooling of the ambient interstellar medium (ISM) in large-scale flows powered by OB cluster feedback can drive the production of dense molecular clouds. We review the current state of the field, with a strong focus on the explicit modelling and observation of the neutral ISM. Magnetohydrodynamic simulations of colliding ISM flows provide a strong theoretical framework in which to view feedback-driven cloud formation, as do models of the gravitational fragmentation of expanding shells. Rapid theoretical developments are accompanied by a growing body of observational work that provides good evidence for the formation of molecular gas via stellar feedback—both in the Milky Way and the Large Magellanic Cloud. The importance of stellar feedback compared with other major astrophysical drivers of dense gas formation remains to be investigated further, and will be an important target for future work.

The APEX-CHAMP+view of the Orion Molecular Cloud 1 core

A high density portion of the Orion Molecular Cloud 1 (OMC-1) contains the prominent, warm Kleinmann-Low (KL) nebula that is internally powered by an energetic event plus a farther region in which intermediate to high mass stars are forming. Its outside is affected by ultraviolet radiation from the neighboring Orion Nebula Cluster and forms the archetypical photon-dominated region (PDR) with the prominent bar feature. Its nearness makes the OMC-1 core region a touchstone for research on the dense molecular interstellar medium and PDRs. Using the Atacama Pathfinder Experiment telescope (APEX), we have imaged the line emission from the multiple transitions of several carbon monoxide (CO) isotopologues over the OMC-1 core region. Our observations employed the 2x7 pixel submillimeter CHAMP+ array to produce maps (~ 300 arcsec x 350 arcsec) of 12CO, 13CO, and C18O from mid-J transitions (J=6-5 to 8-7). We also obtained the 13CO and C18O J=3-2 images toward this region. The 12CO line emission shows a well-defined structure which is shaped and excited by a variety of phenomena, including the energetic photons from hot, massive stars in the nearby Orion Nebula's central Trapezium cluster, active high- and intermediate-mass star formation, and a past energetic event that excites the KL nebula. Radiative transfer modeling of the various isotopologic CO lines implies typical H2 densities in the OMC-1 core region of ~10^4-10^6 cm^-3 and generally elevated temperatures (~ 50-250 K). We estimate a warm gas mass in the OMC-1 core region of 86-285 solar masses.

The dynamics of the ionized and molecular interstellar medium in powerful obscured quasars at z≥ 3.5★

We present an analysis of the kinematics and excitation of the warm ionized gas in two obscured, powerful quasars at z>=3.5 from the SWIRE survey, SWIRE J022513.90-043419.9 and SWIRE J022550.67-042142, based on imaging spectroscopy on the VLT. Line ratios in both targets are consistent with luminous narrow-line regions of AGN. SWIRE J022550.67-042142 has very broad (FWHM=2000 km/s), spatially compact [OIII] line emission. SWIRE J022513.90-043419.9 is spatially resolved, has complex line profiles of H-beta and [OIII], including broad wings with blueshifts of up to -1500 km/s relative to the narrow [OIII]5007 component, and widths of up to FWHM=5000 km/s. Estimating the systemic redshift from the narrow H-beta line, as is standard for AGN host galaxies, implies that a significant fraction of the molecular gas is blueshifted by up to ~ -1000 km/s relative to the systemic velocity. Thus the molecular gas could be participating in the outflow. Significant fractions of the ionized and molecular gas reach velocities greater than the escape velocity. We compare empirical and modeling constraints for different energy injection mechanisms, such as merging, star formation, and momentum-driven AGN winds. We argue that the radio source is the most likely culprit, in spite of the sources rather modest radio power of 10^25 W/Hz. Such a radio power is not uncommon for intense starburst galaxies at z~2. We discuss these results in light of the co-evolution of AGN and their host galaxy.

H3O+line emission from starbursts and AGNs

We use the H3O+ molecule to investigate the impact of starburst and AGN activity on the chemistry of the molecular interstellar medium. Using the JCMT, we have observed the 3+_2 - 2-_2 364 GHz line of p-H3O+ towards the centers of seven active galaxies. We have detected p-H3O+ towards IC342, NGC253, NGC1068, NGC4418, and NGC6240. Upper limits were obtained for IRAS15250 and Arp299. We find large H3O+ abundances (N(H3O+)/N(H2)>10^{-8}) in all detected galaxies apart from in IC342 where it is about one order of magnitude lower. We note, however, that uncertainties in N(H3O+) may be significant due to lack of definite information on source size and excitation. We furthermore compare the derived N(H3O+) with N(HCO+) and find that the H3O+ to HCO+ column density ratio is large in NGC1068 (24), moderate in NGC4418 and NGC253 (4-5), slightly less than unity in NGC6240 (0.7) and lowest in IC342 (0.2-0.6). We compare our results with models of X-ray and photon dominated regions (XDRs and PDRs). For IC342 we find that a starburst PDR chemistry can explain the observed H3O+ abundance. For the other galaxies, the large H3O+ columns are generally consistent with XDR models. In particular for NGC1068 the elevated N(H3O+)/N(HCO+) ratio suggests a low column density XDR. For NGC4418 however, large HC3N abundances are inconsistent with the XDR interpretation. An alternative possibility is that H3O+ forms through H2O evaporating off dust grains and reacting with HCO+ in warm, dense gas. This scenario could also potentially fit the results for NGC253. Further studies of the excitation and distribution of H3O+ - as well as Herschel observations of water abundances - will help to further constrain the models.

Star formation in extreme environments: the effects of cosmic rays and mechanical heating

Context. The molecular interstellar medium in extreme environments, such as Arp 220, but also NGC 253 appears to have extremely high cosmic ray (CR) rates (103−104 × Milky Way) and substantial mechanical heating from supernova driven turbulence.

A CENSUS OF THE HIGH-DENSITY MOLECULAR GAS IN M82

We present a three-pointing study of the molecular gas in the starburst nucleus of M82 based on 190 - 307 GHz spectra obtained with Z-Spec at the Caltech Submillimeter Observatory. We present intensity measurements, detections and upper limits, for 20 transitions, including several new detections of CS, HNC, C2H, H2CO, and CH3CCH lines. We combine our measurements with previously-published measurements at other frequencies for HCN, HNC, CS, C34S, and HCO+ in a multi-species likelihood analysis constraining gas mass, density and temperature, and the species' relative abundances. We find some 1.7 - 2.7 x 10^8 M_sun of gas with n_H2 between 1 - 6 x 10^4 cm^-3 and T > 50 K. While the mass and temperature are comparable to values inferred from mid-J CO transitions, the thermal pressure is a factor of 10 - 20 greater. The molecular interstellar medium is largely fragmented and is subject to ultraviolet irradiation from the star clusters. It is also likely subject to cosmic rays and mechanical energy input from the supernovae, and is warmer on average than the molecular gas in the massive star formation regions in the Milky Way. The typical conditions in the dense gas in M82's central kpc appear unfavorable for further star formation; if any appreciable stellar populations are currently forming, they are likely biased against low mass stars, producing a top-heavy initial mass function.

The molecular interstellar medium of the Local Group dwarf NGC 6822

Do molecular clouds collapse to form stars at the same rate in all environments? In large spiral galaxies, the rate of transformation of H2 into stars (hereafter SFE) varies little. However, the SFE in distant objects (z~1) is much higher than in the large spiral disks that dominate the local universe. Some small local group galaxies share at least some of the characteristics of intermediate-redshift objects, such as size or color. Recent work has suggested that the Star Formation Efficiency (SFE, defined as the SFRate per unit H2) in local Dwarf galaxies may be as high as in the distant objects. A fundamental difficulty in these studies is the independent measure of the H2 mass in metal-deficient environments. At 490 kpc, NGC6822 is an excellent choice for this study; it has been mapped in the CO(2-1) line using the multibeam receiver HERA on the 30 meter IRAM telescope, yielding the largest sample of giant molecular clouds (GMCs) in this galaxy. Despite the much lower metallicity, we find no clear difference in the properties of the GMCs in NGC 6822 and those in the Milky Way except lower CO luminosities for a given mass. Several independent methods indicate that the total H2 mass in NGC 6822 is about 5 x 10^6 Msun in the area we mapped and less than 10^7 Msun in the whole galaxy. This corresponds to a NH2/ICO ~ 4 x 10^{21} cm^-2 /(Kkm/s) over large scales, such as would be observed in distant objects, and half that in individual GMCs. No evidence was found for H2 without CO emission. Our simulations of the radiative transfer in clouds are entirely compatible with these NH2/ICO values. The SFE implied is a factor 5 - 10 higher than what is observed in large local universe spirals.

BONDI–HOYLE–LYTTLETON ACCRETION ONTO A PROTOPLANETARY DISK

Young stellar systems orbiting in the potential of their birth cluster can accrete from the dense molecular interstellar medium during the period between the star's birth and the dispersal of the cluster's gas. Over this time, which may span several Myr, the amount of material accreted can rival the amount in the initial protoplanetary disk; the potential importance of this `tail-end' accretion for planet formation was recently highlighted by Throop & Bally (2008). While accretion onto a point mass is successfully modeled by the classical Bondi-Hoyle-Lyttleton solutions, the more complicated case of accretion onto a star-disk system defies analytic solution. In this paper we investigate via direct hydrodynamic simulations the accretion of dense interstellar material onto a star with an associated gaseous protoplanetary disk. We discuss the changes to the structure of the accretion flow caused by the disk, and vice versa. We find that immersion in a dense accretion flow can redistribute disk material such that outer disk migrates inwards, increasing the inner disk surface density and reducing the outer radius. The accretion flow also triggers the development of spiral density features, and changes to the disk inclination. The mean accretion rate onto the star remains roughly the same with and without the presence of a disk. We discuss the potential impact of this process on planet formation, including the possibility of triggered gravitational instability; inclination differences between the disk and the star; and the appearance of spiral structure in a gravitationally stable system.

Turbulent driving scales in molecular clouds

<i>Context. <i/> Supersonic turbulence in molecular clouds is a dominant agent that strongly affects the clouds' evolution and star formation activity. Turbulence may be initiated and maintained by a number of processes, acting at a wide range of physical scales. By examining the dynamical state of molecular clouds, it is possible to assess the primary candidates for how the turbulent energy is injected. <i>Aims. <i/> The aim of this paper is to constrain the scales at which turbulence is driven in the molecular interstellar medium, by comparing simulated molecular spectral line observations of numerical magnetohydrodynamic models and molecular spectral line observations of real molecular clouds.<i>Methods. <i/> We use principal component analysis, applied to both models and observational data, to extract a quantitative measure of the driving scale of turbulence. <i>Results. <i/> We find that only models driven at large scales (comparable to, or exceeding, the size of the cloud) are consistent with observations. This result applies also to clouds with little or no internal star formation activity.<i>Conclusions. <i/> Astrophysical processes acting on large scales, including supernova-driven turbulence, magneto-rotational instability, or spiral shock forcing, are viable candidates for the generation and maintenance of molecular cloud turbulence. Small-scale driving by sources internal to molecular clouds, such as outflows, can be important on small scales, but cannot replicate the observed large-scale velocity fluctuations in the molecular interstellar medium.

Detection of the 13 CO [FORMULA][F]J=6→5[/F][/FORMULA] transition in the Starburst Galaxy NGC 253

We report the detection of 13CO J = 6→ 5 emission from the nucleus of the starburst galaxy NGC 253 with the redshift (z) and Early Universe Spectrometer (ZEUS), a new submillimeter grating spectrometer. This is the first extragalactic detection of the 13CO J = 6→ 5 transition, which traces warm, dense molecular gas. We employ a multiline LVG analysis and find ≈35%-60% of the molecular interstellar medium is both warm (T ∼ 110 K) and dense (nH2 ∼ 104 cm−3). We analyze the potential heat sources and conclude that ultraviolet and X-ray photons are unlikely to be energetically important. Instead, the molecular gas is most likely heated by an elevated density of cosmic rays or by the decay of supersonic turbulence through shocks. If the cosmic rays and turbulence are created by stellar feedback within the starburst, then our analysis suggests the starburst may be self-limiting.

Simulating High-Redshift Disk Galaxies: Applications to Long Duration Gamma-Ray Burst Hosts

AbstractThe efficiency of star formation governs many observable properties of the cosmological galaxy population, yet many current models of galaxy formation largely ignore the important physics of star formation and the interstellar medium (ISM). Using hydrodynamical simulations of disk galaxies that include a treatment of the molecular ISM and star formation in molecular clouds (Robertson &amp; Kravtsov 2008), we study the influence of star formation efficiency and molecular hydrogen abundance on the properties of high-redshift galaxy populations. In this work, we focus on a model of low-mass, star forming galaxies at 1 ≲ z ≲ 2 that may host long duration gamma-ray bursts (GRBs). Observations of GRB hosts have revealed a population of faint systems with star formation properties that often differ from Lyman-break galaxies (LBGs) and more luminous high-redshift field galaxies. Observed GRB sightlines are deficient in molecular hydrogen, but it is unclear to what degree this deficiency owes to intrinsic properties of the galaxy or the impact the GRB has on its environment. We find that hydrodynamical simulations of low-stellar mass systems at high-redshifts can reproduce the observed star formation rates and efficiencies of GRB host galaxies at redshifts 1 ≲ z ≲ 2. We show that the compact structure of low-mass high-redshift GRB hosts may lead to a molecular ISM fraction of a few tenths, well above that observed in individual GRB sightlines. However, the star formation rates of observed GRB host galaxies imply molecular gas masses of 108 – 109M⊙ similar to those produced in the simulations, and may therefore imply fairly large average H2 fractions in their ISM.