CO Rotational Lines Research Articles

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.

ALMA detection of CO rotational line emission in red supergiant stars of the massive young star cluster RSGC1

Context. The fate of stars largely depends on the amount of mass lost during the end stages of evolution. For single stars with an initial mass between ∼8–30 M⊙, most mass is lost during the red supergiant (RSG) phase, when stellar winds deplete the H-rich envelope. However, the RSG mass-loss rate (Ṁ) is poorly understood theoretically, and so stellar evolution models rely on empirically derived mass-loss rate prescriptions. However, it has been shown that these empirical relations differ largely, with differences up to 2 orders of magnitude. Aims. We aim to derive a new mass-loss rate prescription for RSGs that is not afflicted with some uncertainties inherent in preceding studies. Methods. We have observed CO rotational line emission towards a sample of RSGs in the open cluster RSGC1 that all are of a similar initial mass. The ALMA CO(2–1) line detections allowed us to retrieve the gas mass-loss rates (ṀCO). In contrast to mass-loss rates derived from the analysis of dust spectral features (ṀSED), the data allowed us a direct determination of the wind velocity and no uncertain dust-to-gas correction factor was needed. Results. Five RSGs in RSGC1 have been detected in CO(2–1). The retrieved ṀCO values are systematically lower than ṀSED. Although only five RSGs in RSGC1 have been detected, the data allow us to propose a new mass-loss rate relation for M-type red supergiants with effective temperatures between ∼3200 and 3800 K that is dependent on the luminosity and initial mass, and that is valid during the phase where nuclear burning determines the evolution along the RSG branch. The new mass-loss rate relation is based on the new ṀCO values for the RSGs in RSGC1 and on prior ṀSED values for RSGs in four clusters, including RSGC1. The new Ṁ-prescription yields a good prediction for the mass-loss rate of some well-known Galactic RSGs that are observed in multiple CO rotational lines, including α Ori, μ Cep and VX Sgr. Moreover, there are indications that a stronger, potentially eruptive, mass-loss process is occurring during some fraction of the RSG lifetime, suggesting that RSGs might experience a phase change in mass loss leading to the wind mass-loss rate dominating the RSG evolution at that stage. Conclusions. Implementing a lower mass-loss rate in evolution codes for massive stars has important consequences as to the nature of their end-state. A reduction of the RSG mass-loss rate implies that quiescent RSG mass loss is not enough to strip a single star’s hydrogen-rich envelope. Upon core collapse such single stars would explode as RSGs. Mass-loss rates of order ∼6 times higher would be needed to strip the H-rich envelope and produce a Wolf-Rayet star while evolving back to the blue side of the Hertzsprung–Russell diagram. Future observations of a larger sample of RSGs in open clusters should allow a more stringent determination of the ṀCO–luminosity relation and a sharper diagnostic as to when the phase change in mass loss is occurring.

SOFIA FEEDBACK Survey: The Pillars of Creation in [C ii] and Molecular Lines

We investigate the physical structure and conditions of photodissociation regions (PDRs) and molecular gas within the Pillars of Creation in the Eagle Nebula using SOFIA FEEDBACK observations of the [C ii] 158 μm line. These observations are velocity resolved to 0.5 km s−1 and are analyzed alongside a collection of complimentary data with similar spatial and spectral resolution: the [O i] 63 μm line, also observed with SOFIA, and rotational lines of CO, HCN, HCO+, CS, and N2H+. Using the superb spectral resolution of SOFIA, APEX, CARMA, and BIMA, we reveal the relationships between the warm PDR and cool molecular gas layers in context of the Pillars’ kinematic structure. We assemble a geometric picture of the Pillars and their surroundings informed by illumination patterns and kinematic relationships and derive physical conditions in the PDRs associated with the Pillars. We estimate an average molecular gas density cm−3 and an average atomic gas density n H ∼ 1.8 × 104 cm−3 and infer that the ionized, atomic, and molecular phases are in pressure equilibrium if the atomic gas is magnetically supported. We find pillar masses of 103, 78, 103, and 18 M ⊙ for P1a, P1b, P2, and P3, respectively, and evaporation times of ∼1–2 Myr. The dense clumps at the tops of the pillars are currently supported by the magnetic field. Our analysis suggests that ambipolar diffusion is rapid and these clumps are likely to collapse within their photoevaporation timescales.

LIMpy: A Semianalytic Approach to Simulating Multiline Intensity Maps at Millimeter Wavelengths

Mapping of multiple lines such as the fine-structure emission from [C ii] (157.7 μm), [O iii] (52 and 88.4 μm), and rotational emission lines from CO are of particular interest for upcoming line intensity mapping (LIM) experiments at millimeter wavelengths, due to their brightness features. Several upcoming experiments aim to cover a broad range of scientific goals, from detecting signatures of the epoch of reionization to the physics of star formation and its role in galaxy evolution. In this paper, we develop a semianalytic approach to modeling line strengths as functions of the star formation rate (SFR) or infrared luminosity based on observations of local and high-z galaxies. This package, LIMpy (Line Intensity Mapping in Python), estimates the intensity and power spectra of [C ii], [O iii], and CO rotational transition lines up to the J levels (1–0) to (13–12) based both on analytic formalism and on simulations. We develop a relation among halo mass, SFR, and multiline intensities that permits us to construct a generic formula for the evolution of several line strengths up to z ∼ 10. We implement a variety of star formation models and multiline luminosity relations to estimate the astrophysical uncertainties on the intensity power spectrum of these lines. As a demonstration, we predict the signal-to-noise ratio of [C ii] detection for an EoR-Spec-like instrument on the Fred Young Submillimeter Telescope. Furthermore, the ability to use any halo catalog allows the LIMpy code to be easily integrated into existing simulation pipelines, providing a flexible tool to study intensity mapping in the context of complex galaxy formation physics.

High‐resolution stimulated Raman spectroscopy of low‐frequency rotational transitions: Obtention of collisional broadening coefficients in the rotational S branch of CO perturbed by H2

AbstractWe demonstrate the application of high‐resolution quasi‐continuous (q‐cw) stimulated Raman spectroscopy (SRS) to the study of low‐frequency, pure rotational transitions of molecules in diluted samples in order to obtain collisional information. The main problems encountered in the few pioneering q‐cw SRS rotational works, namely weak signals because of the small population differences between the upper and lower states of the transitions and the presence of alternating current (AC) Stark effect, have been addressed: technical improvements in the optical elements that perform the beam mixing, separation, and filtering have allowed us to use a collinear beam geometry and a triple pass configuration to enhance the Raman signals. The distortions in the line profiles caused by AC Stark effect have been measured and included in the numerical treatment of the spectra. With these extended capabilities, we have been able to register stimulated Raman transitions with wavenumbers as low as ∼ 11.5 cm−1, corresponding to the S(0) pure rotational transition of the molecule of CO. We have applied the setup to the study of collisional broadening of CO rotational lines perturbed by H2 and have obtained broadening coefficients for lines from S(0) to S(20) at 77, 180, and 298 K. A comparison between these anisotropic coefficients and those obtained from IR rovibrational experiments shows excellent agreement. This is, to our knowledge, the first demonstration of the use of rotational SRS for the obtention of collisional information in diluted gas samples.

PHANGS–JWST First Results: A Global and Moderately Resolved View of Mid-infrared and CO Line Emission from Galaxies at the Start of the JWST Era

We explore the relationship between mid-infrared (mid-IR) and CO rotational line emission from massive star-forming galaxies, which is one of the tightest scalings in the local universe. We assemble a large set of unresolved and moderately (∼1 kpc) spatially resolved measurements of CO (1–0) and CO (2–1) intensity, I CO, and mid-IR intensity, I MIR, at 8, 12, 22, and 24 μm. The I CO versus I MIR relationship is reasonably described by a power law with slopes 0.7–1.2 and normalization I CO ∼ 1 K km s−1 at I MIR ∼ 1 MJy sr−1. Both the slopes and intercepts vary systematically with choice of line and band. The comparison between the relations measured for CO (1–0) and CO (2–1) allow us to infer that , in good agreement with other work. The 8 μm and 12 μm bands, with strong polycyclic aromatic hydrocarbon (PAH) features, show steeper CO versus mid-IR slopes than the 22 and 24 μm, consistent with PAH emission arising not just from CO-bright gas but also from atomic or CO-dark gas. The CO-to-mid-IR ratio correlates with global galaxy stellar mass (M ⋆) and anticorrelates with star formation rate/M ⋆. At ∼1 kpc resolution, the first four PHANGS–JWST targets show CO-to-mid-IR relationships that are quantitatively similar to our larger literature sample, including showing the steep CO-to-mid-IR slopes for the JWST PAH-tracing bands, although we caution that these initial data have a small sample size and span a limited range of intensities.

Self-generated ultraviolet radiation in molecular shock waves

Context. The energetics and physical conditions of the interstellar medium and feedback processes remain challenging to probe. Aims. Shocks, modelled over a broad range of parameters, are used to construct a new tool to deduce the mechanical energy and physical conditions from observed atomic or molecular emission lines. Methods. We compute magnetised, molecular shock models with velocities Vs = 5–80 km s−1, pre-shock proton densities nH = 102–106 cm−3, weak or moderate magnetic field strengths, and in the absence or presence of an external UV radiation field. These parameters represent the broadest published range of physical conditions for molecular shocks. As a key shock tracer, we focus on the production of CH+ and post-process the radiative transfer of its rovibrational lines. We develop a simple emission model of an ensemble of shocks for connecting any observed emission lines to the mechanical energy and physical conditions of the system. Results. For this range of parameters, we find the full diversity (C-, C*-, CJ-, and J-type) of magnetohydrodynamic shocks. H2 and H are dominant coolants, with up to 30% of the shock kinetic flux escaping in Lyα photons. The reformation of molecules in the cooling tail means H2 is even a good tracer of dissociative shocks and shocks that were initially fully atomic. The known shock tracer CH+ can also be a significant coolant, reprocessing up to 1% of the kinetic flux. Its production and excitation is intimately linked to the presence of H2 and C+. For each shock model we provide integrated intensities of rovibrational lines of H2, CO, and CH+, and atomic H lines, and atomic fine-structure and metastable lines. We demonstrate how to use these shock models to deduce the mechanical energy and physical conditions of extragalactic environments. As a template example, we interpret the CH+(1−0) emission from the Eyelash starburst galaxy. A mechanical energy injection rate of at least 1011 L⊙ into molecular shocks is required to reproduce the observed line. We find that shocks with velocities as low as 5 km s−1 irradiated by a strong UV field are compatible with the available energy budget. The low-velocity, externally irradiated shocks are at least an order magnitude more efficient than the most efficient shocks with no external irradiation in terms of the total mechanical energy required. We predict differences of more than two orders of magnitude in the intensities of the pure rotational lines of CO, Lyα, and the metastable lines of O, S+, and N between representative models of low-velocity (Vs ~ 10 km s−1) externally irradiated shocks and higher-velocity shocks (Vs ≥ 50 km s−1) with no external irradiation. Conclusions. Shock modelling over an extensive range of physical conditions allows for the interpretation of challenging observations of broad line emission from distant galaxies. Our new method opens up a promising avenue to quantitatively probe the physical conditions and mechanical energy of galaxy-scale gas flows.

Precision tests of CO and [CII] power spectra models against simulated intensity maps

Line intensity mapping (LIM) is an emerging technique with a unique potential to probe a wide range of scales and redshifts. Realizing the full potential of LIM, however, relies on accurate modeling of the signal. We introduce an extended halo model for the power spectrum of intensity fluctuations of CO rotational lines and [CII] fine transition line in real space, modeling nonlinearities in matter fluctuations and biasing relation between the line intensity fluctuations and the underlying dark matter distribution. We also compute the stochastic contributions beyond the Poisson approximation using the halo model framework. To establish the accuracy of the model, we create the first cosmological-scale simulations of CO and [CII] intensity maps, MithraLIMSims, at redshifts 0.5 ≤ z≤6, using halo catalogs from Hidden-Valley simulations, and painting halos according to mass-redshift-luminosity relations for each line. We show that at z=1 on scales kmax≲ 0.8 Mpc-11h, the model predictions of clustering power (with only two free parameters) are in agreement with the measured power spectrum at better than 5%. At higher redshift of z=4.5, this remarkable agreement extends to smaller scale of kmax≲ 2 Mpc-11h. Furthermore, we show that on large scales, the stochastic contributions to CO and CII power spectra are non-Poissonian, with amplitudes reproduced reasonably well by the halo model prescription. Lastly, we assess the performance of the theoretical model of the baryon acoustic oscillations (BAO) and show that hypothetical LIM surveys probing CO lines at z=1, that can be deployed within this decade, will be able to make a high significance measurement of the BAO. On a longer time scale, a space-based mission probing [CII] line can uniquely measure the BAO on a wide range of redshifts at an unprecedented precision.

Neutrino Properties with Ground-based Millimeter-wavelength Line Intensity Mapping

Abstract Line intensity mapping (LIM) is emerging as a powerful technique to map the cosmic large-scale structure and to probe cosmology over a wide range of redshifts and spatial scales. We perform Fisher forecasts to determine the optimal design of wide-field ground-based millimeter-wavelength LIM surveys for constraining properties of neutrinos and light relics. We consider measuring the auto-power spectra of several CO rotational lines (from J = 2–1 to J = 6–5) and the [C ii] fine-structure line in the redshift range of 0.25 < z < 12. We study the constraints with and without interloper lines as a source of noise in our analysis, and for several one-parameter and multiparameter extensions of ΛCDM. We show that LIM surveys deployable this decade, in combination with existing cosmic microwave background (CMB; primary) data, could achieve order-of-magnitude improvements over Planck constraints on N eff and M ν . Compared to next-generation CMB and galaxy surveys, a LIM experiment of this scale could achieve bounds that are a factor of ∼3 better than those forecasted for surveys such as EUCLID (galaxy clustering), and potentially exceed the constraining power of CMB-S4 by a factor of ∼1.5 and ∼3 for N eff and M ν , respectively. We show that the forecasted constraints are not substantially affected when enlarging the parameter space, and additionally demonstrate that such a survey could also be used to measure ΛCDM parameters and the dark energy equation of state exquisitely well.

Observability of the vertical shear instability in protoplanetary disk CO kinematics

Context. Dynamical and turbulent motions of gas in a protoplanetary disk are crucial for their evolution and are thought to affect planet formation. Recent (sub-)millimeter observations show evidence of weak turbulence in the disk’s outer regions. However, the detailed physical mechanism of turbulence in these outer regions remains uncertain. The vertical shear instability (VSI) is a promising candidate mechanism to produce turbulence in the outer parts of the disk. Aims. Our objective is to study the observability of the gas velocity structure produced by the VSI via CO kinematics with the Atacama Large Millimetre/submillimetre Array (ALMA). Methods. We performed global 3D hydrodynamical simulations of an inviscid and locally isothermal VSI-unstable disk. We post-processed the simulation results with radiative transfer calculations and produced synthetic predictions of CO rotational emission lines. Next, we computed the line of sight velocity map and its deviations from a sub-Keplerian equilibrium solution. We explored the detectability of the VSI by identifying kinematic signatures using realistic simulated observations using the CASA package. Results. Our 3D hydrodynamical simulations of the VSI show the steady state dynamics of the gas in great detail. From the velocity structure, we infer a turbulent stress value of αrϕ = 1.4 × 10−4. On large scales, we observe clear velocity deviations of the order of 50 m s−1 as axisymmetric rings with radially interspersed signs. By comparing synthetic observations at different inclinations we find optimal conditions at i ≲ 20° to trace for the kinematic structures of the VSI. We found that current diagnostics to constrain gas turbulence from nonthermal broadening of the molecular line emission are not applicable to anisotropic VSI turbulence. Conclusions. We conclude that the detection of kinematic signatures produced by the VSI is possible with ALMA’s current capabilities. Observations including an extended antenna configuration are required to resolve the structure (beam sizes below ~10 au). The highest spectral resolution available is needed (~0.05 km s−1 with ALMA Band 6) for a robust detection. The characterization of the large-scale velocity perturbations is required to constrain the turbulence level produced by the VSI from gas observations.

Probing Cosmic Reionization and Molecular Gas Growth with TIME

Abstract Line intensity mapping (LIM) provides a unique and powerful means to probe cosmic structures by measuring the aggregate line emission from all galaxies across redshift. The method is complementary to conventional galaxy redshift surveys that are object based and demand exquisite point-source sensitivity. The Tomographic Ionized-carbon Mapping Experiment (TIME) will measure the star formation rate during cosmic reionization by observing the redshifted [C ii] 158 μm line (6 ≲ z ≲ 9) in the LIM regime. TIME will simultaneously study the abundance of molecular gas during the era of peak star formation by observing the rotational CO lines emitted by galaxies at 0.5 ≲ z ≲ 2. We present the modeling framework that predicts the constraining power of TIME on a number of observables, including the line luminosity function and the auto- and cross-correlation power spectra, including synergies with external galaxy tracers. Based on an optimized survey strategy and fiducial model parameters informed by existing observations, we forecast constraints on physical quantities relevant to reionization and galaxy evolution, such as the escape fraction of ionizing photons during reionization, the faint-end slope of the galaxy luminosity function at high redshift, and the cosmic molecular gas density at cosmic noon. We discuss how these constraints can advance our understanding of cosmological galaxy evolution at the two distinct cosmic epochs for TIME, starting in 2021, and how they could be improved in future phases of the experiment.

Resolving the Dust-to-Metals Ratio and CO-to-H2 Conversion Factor in the Nearby Universe

Abstract We investigate the relationship between the dust-to-metals ratio (D/M) and the local interstellar medium environment at ∼2 kpc resolution in five nearby galaxies: IC 342, M31, M33, M101, and NGC 628. A modified blackbody model with a broken power-law emissivity is used to model the dust emission from 100 to 500 μm observed by Herschel. We utilize the metallicity gradient derived from auroral line measurements in H i regions whenever possible. Both archival and new CO rotational line and H i 21 cm maps are adopted to calculate gas surface density, including new wide-field CO and H i maps for IC 342 from IRAM and the VLA, respectively. We experiment with several prescriptions of the CO-to-H2 conversion factor and compare the resulting D/M–metallicity and D/M–density correlations, both of which are expected to be nonnegative from depletion studies. The D/M is sensitive to the choice of the conversion factor. The conversion factor prescriptions based on metallicity only yield too much molecular gas in the center of IC 342 to obtain the expected correlations. Among the prescriptions tested, the one that yields the expected correlations depends on both metallicity and surface density. The 1σ range of the derived D/M spans 0.40–0.58. Compared to chemical evolution models, our measurements suggest that the dust growth timescale is much shorter than the dust destruction timescale. The measured D/M is consistent with the D/M in galaxy-integrated studies derived from infrared dust emission. Meanwhile, the measured D/M is systematically higher than the D/M derived from absorption, which likely indicates a systematic offset between the two methods.

Cosmic Expansion History from Line-Intensity Mapping

Line-intensity mapping (LIM) of emission from star-forming galaxies can be used to measure the baryon acoustic oscillation (BAO) scale as far back as the epoch of reionization. This provides a standard cosmic ruler to constrain the expansion rate of the Universe at redshifts which cannot be directly probed otherwise. In light of growing tension between measurements of the current expansion rate using the local distance ladder and those inferred from the cosmic microwave background, extending the constraints on the expansion history to bridge between the late and early Universe is of paramount importance. Using a newly derived methodology to robustly extract cosmological information from LIM, which minimizes the inherent degeneracy with unknown astrophysics, we show that present and future experiments can gradually improve the measurement precision of the expansion rate history, ultimately reaching percent-level constraints on the BAO scale. Specifically, we provide detailed forecasts for the SPHEREx satellite, which will target the Hα and Lyman-α lines, for a near-future stage-2 experiment targeting CII, and for the ground-based COMAP instrument-as well as a future stage-3 experiment-that will target the CO rotational lines. Besides weighing in on the so-called Hubble tension, reliable LIM cosmic rulers can enable wide-ranging tests of dark matter, dark energy, and modified gravity.

Radiation thermo-chemical models of protoplanetary disks

Context. Disks around pre-main-sequence stars evolve over time by turbulent viscous spreading. The main contender to explain the strength of the turbulence is the magnetorotational instability model, whose efficiency depends on the disk ionization fraction. Aims. Our aim is to compute self-consistently the chemistry including polycyclic aromatic hydrocarbon (PAH) charge chemistry, the grain charging, and an estimate of an effective value of the turbulence α parameter in order to find observational signatures of disk turbulence. Methods. We introduced PAH and grain charging physics and their interplay with other gas-phase reactions in the physico-chemical code PRODIMO. Non-ideal magnetohydrodynamics effects such as ohmic and ambipolar diffusion are parametrized to derive an effective value for the turbulent parameter αeff. We explored the effects of turbulence heating and line broadening on CO isotopologue submillimeter lines. Results. The spatial distribution of αeff depends on various unconstrained disk parameters such as the magnetic parameter βmag or the cosmic ray density distribution inside the protoplanetary disk s. The inner disk midplane shows the presence of the so-called dead zone where the turbulence is almost inexistent. The disk is heated mostly by thermal accommodation on dust grains in the dead zone, by viscous heating outside the dead zone up to a few hundred astronomical units, and by chemical heating in the outer disk. The CO rotational lines probe the warm molecular disk layers where the turbulence is at its maximum. However, the effect of turbulence on the CO line profiles is minimal and difficult to distinguish from the thermal broadening. Conclusions. Viscous heating of the gas in the disk midplane outside the dead zone is efficient. The determination of α from CO rotational line observations alone is challenging.

An ALMA Multiline Survey of the Interstellar Medium of the Redshift 7.5 Quasar Host Galaxy J1342+0928

Abstract We use Atacama Large Millimeter Array observations of the host galaxy of the quasar ULAS J1342+0928 at z = 7.54, to study the dust continuum and far-infrared lines emitted from its interstellar medium (ISM). The Rayleigh–Jeans tail of the dust continuum is well sampled with eight different spectral setups, and from a modified blackbody fit we obtain an emissivity coefficient of β = 1.85 ± 0.3. Assuming a standard dust temperature of 47 K we derive a dust mass of M dust = 0.35 × 108 M ☉ and a star formation rate of . We have >4σ detections of the , , and atomic fine structure lines and limits on the , , and emission. We also report multiple limits of CO rotational lines with J up ≥ 7, as well as a tentative 3.3σ detection of the stack of four CO lines (J up = 11, 10, 8, and 7). We find line deficits that are in agreement with local ultra-luminous infrared galaxies. Comparison of the and lines indicates that the emission arises predominantly from the neutral medium, and we estimate that the photodisassociation regions in J1342+0928 have densities ≲5 × 104 cm−3. The data suggest that ∼16% of hydrogen is in ionized form and that the H ii regions have high electron densities of n e > 180 cm−3. Our observations favor a low gas-to-dust ratio of <100, and a metallicity of the ISM comparable to the solar value. All the measurements presented here suggest that the host galaxy of J1342+0928 is highly enriched in metal and dust, despite being observed just 680 Myr after the big bang.

More than star formation: High-J CO SLEDs of high-z galaxies

AbstractTheoretical work suggests that AGNs play an important role in quenching star formation in massive galaxies. In addition to molecular outflows observed in the local universe, emission from very high-J CO rotational transitions have been a key piece of evidence for AGN directly affecting the molecular gas reservoirs that fuel star formation. However, very few observations exist of CO rotational lines past the peak of the CO spectral line energy distribution (SLED) for galaxies in the early universe. Here we present new ALMA observations of high-J CO rotational lines (from CO(5–4) to CO(16–15)) in six z > 2 IR-bright systems, including several sources not known to contain a strong AGN for comparison. We detect significant amounts of high-excitation CO emission that suggests the presence of energy sources beyond UV-heating.

Dark Molecular Gas in Simulations of z ∼ 0 Disk Galaxies

Abstract The H2 mass of molecular clouds has traditionally been traced by the CO(J = 1−0) rotational transition line. This said, CO is relatively easily photodissociated and can also be destroyed by cosmic rays, thus rendering some fraction of molecular gas to be “CO-dark.” We investigate the amount and physical properties of CO-dark gas in two z ∼ 0 disk galaxies and develop predictions for the expected intensities of promising alternative tracers ([C i] 609 μm and [C ii] 158 μm emission). We do this by combining cosmological zoom simulations of disk galaxies with thermal-radiative-chemical equilibrium interstellar medium (ISM) calculations to model the predicted H i and H2 abundances and CO (J = 1−0), [C i] 609 μm, and [C ii] 158 μm emission properties. Our model treats the ISM as a collection of radially stratified clouds whose properties are dictated by their volume and column densities, the gas-phase metallicity, and the interstellar radiation field (ISRF) and CR ionization rates. Our main results follow. Adopting an observationally motivated definition of CO-dark gas, i.e., H2 gas with W CO < 0.1 K km s−1, we find that a significant amount (≳50%) of the total H2 mass lies in CO-dark gas, most of which is diffuse gas, poorly shielded due to low dust column density. The CO-dark molecular gas tends to be dominated by [C ii], though [C i] also serves as a bright tracer of the dark gas in many instances. At the same time, [C ii] also tends to trace neutral atomic gas. As a result, when we quantify the conversion factors for the three carbon-based tracers of molecular gas, we find that [C i] suffers the least contamination from diffuse atomic gas and is relatively insensitive to secondary parameters.

Planck’s dusty GEMS

We present an extensive CO emission-line survey of the Planck’s dusty Gravitationally Enhanced subMillimetre Sources, a small set of 11 strongly lensed dusty star-forming galaxies at z = 2–4 discovered with Planck and Herschel satellites, using EMIR on the IRAM 30-m telescope. We detected a total of 45 CO rotational lines from Jup = 3 to Jup = 11, and up to eight transitions per source, allowing a detailed analysis of the gas excitation and interstellar medium conditions within these extremely bright (μLFIR = 0.5 − 3.0 × 1014L⊙), vigorous starbursts. The peak of the CO spectral-line energy distributions (SLEDs) fall between Jup = 4 and Jup = 7 for nine out of 11 sources, in the same range as other lensed and unlensed submillimeter galaxies (SMGs) and the inner regions of local starbursts. We applied radiative transfer models using the large velocity gradient approach to infer the spatially-averaged molecular gas densities, nH2 ≃ 102.6 − 104.1 cm−3, and kinetic temperatures, Tk ≃ 30–1000 K. In five sources, we find evidence of two distinct gas phases with different properties and model their CO SLED with two excitation components. The warm (70–320 K) and dense gas reservoirs in these galaxies are highly excited, while the cooler (15–60 K) and more extended low-excitation components cover a range of gas densities. In two sources, the latter is associated with diffuse Milky Way-like gas phases of density nH2 ≃ 102.4 − 102.8 cm−3, which provides evidence that a significant fraction of the total gas masses of dusty starburst galaxies can be embedded in cool, low-density reservoirs. The delensed masses of the warm star-forming molecular gas range from 0.6to12 × 1010 M⊙. Finally, we show that the CO line luminosity ratios are consistent with those predicted by models of photon-dominated regions (PDRs) and disfavor scenarios of gas clouds irradiated by intense X-ray fields from active galactic nuclei. By combining CO, [C I] and [C II] line diagnostics, we obtain average PDR gas densities significantly higher than in normal star-forming galaxies at low-redshift, as well as far-ultraviolet radiation fields 102–104 times more intense than in the Milky Way. These spatially-averaged conditions are consistent with those in high-redshift SMGs and in a range of low-redshift environments, from the central regions of ultra-luminous infrared galaxies and bluer starbursts to Galactic giant molecular clouds.

Warm CO in evolved stars from the THROES catalogue

In this work (Paper I), we analyse Herschel-PACS spectroscopy for a subsample of 23 O-rich and 3 S-type evolved stars, in different evolutionary stages from the asymptotic giant branch (AGB) to the planetary nebula (PN) phase, from the THROES catalogue. (C-rich targets are separately studied in Paper II). The broad spectral range covered by PACS (∼55–210 μm) includes a large number of high-J CO lines, from J = 14 − 13 to J = 45 − 44 (v = 0), that allow us to study the warm inner layers of the circumstellar envelopes (CSEs) of these objects, at typical distances from the star of ≈1014–1015 cm and ≈1016 cm for AGBs and post-AGB-PNe, respectively. We have generated CO rotational diagrams for each object to derive the rotational temperature, total mass within the CO-emitting region and average mass-loss rate during the ejection of these layers. We present first order estimations of these basic physical parameters using a large number of high-J CO rotational lines, with upper-level energies from Eup ∼ 580 to 5000 K, for a relatively big set of evolved low-to-intermediate mass stars in different AGB-to-PN evolutionary stages. We derive rotational temperatures ranging from Trot ∼ 200 to 700 K, with typical values around 500 K for AGBs and systematically lower, ∼200 K, for objects in more advanced evolutionary stages (post-AGBs and PNe). Our values of Trot are one order or magnitude higher than the temperatures of the outer CSE layers derived from low-J CO line studies. The total mass of the inner CSE regions where the PACS CO lines arise is found to range from MH2 ∼ 10−6 to ≈10−2 M⊙, which is expected to represent a small fraction of the total CSE mass. The mass-loss rates estimated are in the range Ṁ ∼ 10−7 − 10−4 M⊙yr−1, in agreement (within uncertainties) with values found in the literature. We find a clear anticorrelation between MH2 and Ṁ vs. Trot that probably results from a combination of most efficient line cooling and higher line opacities in high mass-loss rate objects. For some strong CO emitters in our sample, a double temperature (hot and warm) component is inferred. The temperatures of the warm and hot components are ∼400–500 K and ∼600–900 K, respectively. The mass of the warm component (∼10−5–8 × 10−2 M⊙) is always larger than that of the hot component, by a factor of between two and ten. The warm-to-hot MH2 and Trot ratios in our sample are correlated and are consistent with an average temperature radial profile of ∝ r−0.5 ± 0.1, that is, slightly shallower than in the outer envelope layers, in agreement with recent studies.

Constraining the expansion history and early dark energy with line intensity mapping

We consider the potential for line intensity mapping (IM) experiments to measure the baryon acoustic oscillations (BAO) from 3 < z < 6. This would constrain the expansion history in a redshift range that is currently unexplored. We calculate the map depths that future IM experiments targeting the CO(1-0) rotational transition line and [CII] ionized carbon fine-structure line would need to achieve in order to measure the BAO. We find that near-future IM experiments could constrain the BAO scale to 5% or better depending on CO/[CII] model amplitude. This measurement is at a precision that could make competitive constraints on models of early dark energy.