skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on November 1, 2025

Title: Quantifying the informativity of emission lines to infer physical conditions in giant molecular clouds: I. Application to model predictions
Context.Observations of ionic, atomic, or molecular lines are performed to improve our understanding of the interstellar medium (ISM). However, the potential of a line to constrain the physical conditions of the ISM is difficult to assess quantitatively, because of the complexity of the ISM physics. The situation is even more complex when trying to assess which combinations of lines are the most useful. Therefore, observation campaigns usually try to observe as many lines as possible for as much time as possible. Aims.We have searched for a quantitative statistical criterion to evaluate the full constraining power of a (combination of) tracer(s) with respect to physical conditions. Our goal with such a criterion is twofold. First, we want to improve our understanding of the statistical relationships between ISM tracers and physical conditions. Secondly, by exploiting this criterion, we aim to propose a method that helps observers to make their observation proposals; for example, by choosing to observe the lines with the highest constraining power given limited resources and time. Methods.We propose an approach based on information theory, in particular the concepts of conditional differential entropy and mutual information. The best (combination of) tracer(s) is obtained by comparing the mutual information between a physical parameter and different sets of lines. The presented analysis is independent of the choice of the estimation algorithm (e.g., neural network orχ2minimization). We applied this method to simulations of radio molecular lines emitted by a photodissociation region similar to the Horsehead Nebula. In this simulated data, we considered the noise properties of a state-of-the-art single dish telescope such as the IRAM 30m telescope. We searched for the best lines to constrain the visual extinction,AVtot, or the ultraviolet illumination field,G0. We ran this search for different gas regimes, namely translucent gas, filamentary gas, and dense cores. Results.The most informative lines change with the physical regime (e.g., cloud extinction). However, the determination of the optimal (combination of) line(s) to constrain a physical parameter such as the visual extinction depends not only on the radiative transfer of the lines and chemistry of the associated species, but also on the achieved mean signal-to-noise ratio. The short integration time of the CO isotopologueJ= 1 − 0 lines already yields much information on the total column density for a large range of (AVtot,G0) space. The best set of lines to constrain the visual extinction does not necessarily combine the most informative individual lines. Precise constraints on the radiation field are more difficult to achieve with molecular lines. They require spectral lines emitted at the cloud surface (e.g., [CII] and [CI] lines). Conclusions.This approach allows one to better explore the knowledge provided by ISM codes, and to guide future observation campaigns.  more » « less
Award ID(s):
2142300
PAR ID:
10569719
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; « less
Publisher / Repository:
Astronomy and Astrophysics
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
691
ISSN:
0004-6361
Page Range / eLocation ID:
A109
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The CO(1–0) line has been carefully calibrated as a tracer of molecular gas mass. However, recent studies often favor higherJtransitions of the CO molecule, which are brighter and accessible for redshift ranges where CO(1–0) is not. These lines are not perfect analogs for CO(1–0), owing to their more stringent excitation conditions, and must be calibrated for use as molecular gas tracers. Here, we introduce the Arizona Molecular ISM Survey with the SMT, a multi-CO line survey ofz∼ 0 galaxies conducted to calibrate the CO(2–1) and CO(3–2) lines. The final survey includes CO(2–1) spectra of 176 galaxies and CO(3–2) spectra for a subset of 45. We supplement these with archival CO(1–0) spectra from xCOLD GASS for all sources and additional CO(1–0) observations with the Kitt Peak 12 m Telescope. Targets were selected to be representative of the 109M≤M*≤ 1011.5Mgalaxy population. Our project emphasized careful characterization of statistical and systematic uncertainties to enable studies of trends in CO line ratios. We show that optical and CO disk sizes are on average equal, for both the CO(1–0) and CO(2–1) line. We measure the distribution of CO line luminosity ratios, finding medians (16th–84th percentile) of 0.71 (0.51–0.96) for the CO(2–1)-to-CO(1–0) ratio, 0.39 (0.24–0.53) for the CO(3–2)-to-CO(1–0) ratio, and 0.53 (0.41–0.74) for the CO(3–2)-to-CO(2–1) ratio. A companion paper presents our study of CO(2–1)'s applicability as a molecular gas mass tracer and search for trends in the CO(2–1)-to-CO(1–0) ratio. Our catalog of CO line luminosities is publicly available. 
    more » « less
  2. Context. Carbon monoxide (CO) is a poor tracer of H2in the diffuse interstellar medium (ISM), where most of the carbon is not incorporated into CO molecules, unlike the situation at higher extinctions. Aims. We present a novel, indirect method for constraining H2column densities (NH2) without employing CO observations. We show that previously recognized nonlinearities in the relation between the extinction,AV(H2), derived from dust emission and the H Icolumn density (NH I) are due to the presence of molecular gas. Methods. We employed archival (NH2) data, obtained from the UV spectra of stars, and calculatedAV(H2) toward these sight lines using 3D extinction maps. The following relation fits the data: logNH2= 1.38742 (logAV(H2))3− 0.05359 (logAV(H2))2+ 0.25722 logAV(H2) + 20.67191. This relation is useful for constrainingNH2in the diffuse ISM as it requires onlyNH Iand dust extinction data, which are both easily accessible. In 95% of the cases, the estimates produced by the fitted equation have deviations of less than a factor of 3.5. We constructed aNH2map of our Galaxy and compared it to the CO integrated intensity (WCO) distribution. Results. We find that the average ratio (XCO) betweenNH2andWCOis approximately equal to 2 × 1020cm−2(K km s−1)−1, consistent with previous estimates. However, we find that theXCOfactor varies by orders of magnitude on arcminute scales between the outer and the central portions of molecular clouds. For regions withNH2≳ 1020cm−2, we estimate that the average H2fractional abundance,fH2= 2NH2/(2NH2+NH I), is 0.25. Multiple (distinct) largely atomic clouds are likely found along high-extinction sightlines (AV≥ 1 mag), hence limitingfH2in these directions. Conclusions. More than 50% of the lines of sight withNH2≥ 1020cm−2are untraceable by CO with aJ= 1−0 sensitivity limitWCO= 1 K km s−1
    more » « less
  3. We present Atacama Large Millimeter/submillimeter Array (ALMA) sub-kiloparsec- to kiloparsec-scale resolution observations of the [C II], CO (9–8), and OH+(11–01) lines along with their dust continuum emission toward the far-infrared (FIR) luminous quasar SDSS J231038.88+185519.7 atz = 6.0031, to study the interstellar medium distribution, the gas kinematics, and the quasar-host system dynamics. We decompose the intensity maps of the [C II] and CO (9–8) lines and the dust continuum with two-dimensional elliptical Sérsic models. The [C II] brightness follows a flat distribution with a Sérsic index of 0.59. The CO (9–8) line and the dust continuum can be fit with an unresolved nuclear component and an extended Sérsic component with a Sérsic index of ∼1, which may correspond to the emission from an active galactic nucleus dusty molecular torus and a quasar host galaxy, respectively. The different [C II] spatial distribution may be due to the effect of the high dust opacity, which increases the FIR background radiation on the [C II] line, especially in the galaxy center, significantly suppressing the [C II] emission profile. The dust temperature drops with distance from the center. The effective radius of the dust continuum is smaller than that of the line emission and the dust mass surface density, but is consistent with that of the star formation rate surface density. This may indicate that the dust emission is a less robust tracer of the dust and gas distribution but is a decent tracer of the obscured star formation activity. The OH+(11–01) line shows a P-Cygni profile with an absorption at ∼–400 km s−1, which may indicate an outflow with a neutral gas mass of (6.2 ± 1.2)×108 Malong the line of sight. We employed a three-dimensional tilted ring model to fit the [C II] and CO (9–8) data cubes. The two lines are both rotation dominated and trace identical disk geometries and gas motions. This suggest that the [C II] and CO (9–8) gas are coplanar and corotating in this quasar host galaxy. The consistent circular velocities measured with [C II] and CO (9–8) lines indicate that these two lines trace a similar gravitational potential. We decompose the circular rotation curve measured from the kinematic model fit to the [C II] line into four matter components (black hole, stars, gas, and dark matter). The quasar-starburst system is dominated by baryonic matter inside the central few kiloparsecs. We constrain the black hole mass to be 2.97+0.51-0.77 × 109M; this is the first time that the dynamical mass of a black hole has been measured atz ∼ 6. This mass is consistent with that determined using the scaling relations from quasar emission lines. A massive stellar component (on the order of 109 M) may have already existed when the Universe was only ∼0.93 Gyr old. The relations between the black hole mass and the baryonic mass of this quasar indicate that the central supermassive black hole may have formed before its host galaxy. 
    more » « less
  4. Abstract We present a new suite of numerical simulations of the star-forming interstellar medium (ISM) in galactic disks using the TIGRESS-NCR framework. Distinctive aspects of our simulation suite are (1) sophisticated and comprehensive numerical treatments of essential physical processes including magnetohydrodynamics, self-gravity, and galactic differential rotation, as well as photochemistry, cooling, and heating coupled with direct ray-tracing UV radiation transfer and resolved supernova feedback and (2) wide parameter coverage including the variation in metallicity over Z Z / Z 0.1 - 3 , gas surface density Σgas∼ 5–150Mpc−2, and stellar surface density Σstar∼ 1–50Mpc−2. The range of emergent star formation rate surface density is ΣSFR∼ 10−4–0.5Mkpc−2yr−1, and ISM total midplane pressure isPtot/kB= 103–106cm−3K, withPtotequal to the ISM weight W . For given Σgasand Σstar, we find Σ SFR Z 0.3 . We provide an interpretation based on the pressure-regulated feedback-modulated (PRFM) star formation theory. The total midplane pressure consists of thermal, turbulent, and magnetic stresses. We characterize feedback modulation in terms of the yield ϒ, defined as the ratio of each stress to ΣSFR. The thermal feedback yield varies sensitively with both weight and metallicity as ϒ th W 0.46 Z 0.53 , while the combined turbulent and magnetic feedback yield shows weaker dependence ϒ turb + mag W 0.22 Z 0.18 . The reduction in ΣSFRat low metallicity is due mainly to enhanced thermal feedback yield, resulting from reduced attenuation of UV radiation. With the metallicity-dependent calibrations we provide, PRFM theory can be used for a new subgrid star formation prescription in cosmological simulations where the ISM is unresolved. 
    more » « less
  5. Understanding the nature of high-redshift dusty galaxies requires a comprehensive view of their interstellar medium (ISM) and molecular complexity. However, the molecular ISM at high redshifts is commonly studied using only a few species beyond12C16O, limiting our understanding. In this paper, we present the results of deep 3 mm spectral line surveys using the NOrthern Extended Millimeter Array (NOEMA) targeting two strongly lensed dusty galaxies observed when the Universe was less than 1.8 Gyr old: APM 08279+5255, a quasar at redshiftz= 3.911, and NCv1.143 (H-ATLAS J125632.7+233625), az= 3.565 starburst galaxy. The spectral line surveys cover rest-frame frequencies from about 330 to 550 GHz for both galaxies. We report the detection of 38 and 25 emission lines in APM 08279+5255 and NCv1.143, respectively. These lines originate from 17 species, namely CO,13CO, C18O, CN, CCH, HCN, HCO+, HNC, CS, C34S, H2O, H3O+, NO, N2H+, CH, c-C3H2, and the vibrationally excited HCN and neutral carbon. The spectra reveal the chemical richness and the complexity of the physical properties of the ISM. By comparing the spectra of the two sources and combining the analysis of the molecular gas excitation, we find that the physical properties and the chemical imprints of the ISM are different: the molecular gas is more excited in APM 08279+5255, which exhibits higher molecular gas temperatures and densities compared to NCv1.143; the molecular abundances in APM 08279+5255 are akin to the values of local active galactic nuclei (AGN), showing boosted relative abundances of the dense gas tracers that might be related to high-temperature chemistry and/or the X-ray-dominated regions, while NCv1.143 more closely resembles local starburst galaxies. The most significant differences between the two sources are found in H2O: the 448 GHz ortho-H2O(423 − 330) line is significantly brighter in APM 08279+5255, which is likely linked to the intense far-infrared radiation from the dust powered by AGN. Our astrochemical model suggests that, at such high column densities, far-ultraviolet radiation is less important in regulating the ISM, while cosmic rays (and/or X-rays and shocks) are the key players in shaping the molecular abundances and the initial conditions of star formation. Both our observed CO isotopologs line ratios and the derived extreme ISM conditions (high gas temperatures, densities, and cosmic-ray ionization rates) suggest the presence of a top-heavy stellar initial mass function. From the ∼330–550 GHz continuum, we also find evidence of nonthermal millimeter flux excess in APM 08279+5255 that might be related to the central supermassive black hole. Such deep spectral line surveys open a new window into the physics and chemistry of the ISM and the radiation field of galaxies in the early Universe. 
    more » « less