CO Column Density Research Articles

Low CI/CO abundance ratio revealed by HST UV spectroscopy of CO-rich debris discs

ABSTRACT The origin and evolution of CO gas in debris discs has been debated since its initial detection. The gas could have a primordial origin, as a remnant of the protoplanetary disc or a secondary exocometary origin. This paper investigates the origin of gas in two debris discs, HD110058 and HD131488, using Hubble Space Telescope (HST) observations of CI and CO, which play critical roles in the gas evolution. We fitted several electronic transitions of CI and CO rovibronic bands to derive column densities and temperatures for each system, revealing high CO column densities (∼3–4 orders of magnitude higher than β Pictoris), and low CI/CO ratios in both. Using the exogas model, we simulated the radial evolution of the gas in the debris disc assuming a secondary gas origin. We explored a wide range of CO exocometary release rates and α viscosities, which are the key parameters of the model. Additionally, we incorporated photodissociation due to stellar ultraviolet (UV) to the exogas model and found that it is negligible for typical CO-rich discs and host stars, even at a few au due to the high radial optical depths in the Extreme ultraviolet radiation (EUV). We find that the current steady-state secondary release model cannot simultaneously reproduce the CO and CI HST-derived column densities, as it predicts larger CI/CO ratios than observed. Our direct UV measurement of low CI/CO ratios agrees with results derived from recent Atacama Large Millimetre/submillimetre Array findings and may point to vertical layering of CI, additional CI removal, CO shielding processes, or different gas origin scenarios.

PRODIGE – planet-forming disks in Taurus with NOEMA

Context. The physics and chemistry of planet-forming disks are far from being fully understood. To make further progress, both broad line surveys and observations of individual tracers in a statistically significant number of disks are required. Aims. Our aim is to perform a line survey of eight planet-forming Class II disks in Taurus with the IRAM NOrthern Extended Millimeter Array (NOEMA), as a part of the MPG-IRAM Observatory Program PRODIGE (PROtostars and DIsks: Global Evolution; PIs: P. Caselli and Th. Henning). Methods. Compact and extended disks around T Tauri stars CI, CY, DG, DL, DM, DN, IQ Tau, and UZ Tau E are observed in ~80 lines from > 20 C-, O,- N-, and S-bearing species. The observations in four spectral settings at 210–280 GHz with a 1σ rms sensitivity of ~8–12 mJy beam−1 at a 0.9″ and 0.3 km s−1 resolution will be completed in 2024. The uv visibilities are fitted with the DiskFit model to obtain key stellar and disk properties. Results. In this first paper, the combined 12CO, 13CO, and C18O J = 2–1 data are presented. We find that the CO fluxes and disk masses inferred from dust continuum tentatively correlate with the CO emission sizes. We constrained dynamical stellar masses, geometries, temperatures, the CO column densities, and gas masses for each disk. The best-fit temperatures at 100 au are ~ 17–37 K, and decrease radially with the power-law exponent q ~ 0.05–0.76. The inferred CO column densities decrease radially with the power-law exponent p ~ 0.2–3.1. The gas masses estimated from 13CO (2–1) are ~0.001–0.2 M⊙. Conclusions. Using NOEMA, we confirm the presence of temperature gradients in our disk sample. The best-fit CO column densities point to severe CO freeze-out in these disks. The DL Tau disk is an outlier, and has either stronger CO depletion or lower gas mass than the rest of the sample. The CO isotopologue ratios are roughly consistent with the observed values in disks and the low-mass star-forming regions. The high 13CO/C18O ratio of ~23 in DM Tau could be indicative of strong selective photodissociation of C18O in this disk.

JWST detection of extremely excited outflowing CO and H2O in VV 114 E SW: A possible rapidly accreting IMBH

Mid-infrared (mid-IR) gas-phase molecular bands are powerful diagnostics of the warm interstellar medium. We report the James Webb Space Telescope detection of the CO v = 1 − 0 (4.4 − 5.0 μm) and H2O ν2 = 1 − 0 (5.0 − 7.8 μm) ro-vibrational bands, both in absorption, toward the “s2” core in the southwest nucleus of the merging galaxy VV 114 E. All ro-vibrational CO lines up to Jlow = 33 (Elow ≈ 3000 K) are detected, as well as a forest of H2O lines up to 130, 13 (Elow ≈ 2600 K). The highest-excitation lines are blueshifted by ∼180 km s−1 relative to the extended molecular cloud, which is traced by the rotational CO (J = 3 − 2) 346 GHz line observed with the Atacama Large Millimeter/submillimeter Array. The bands also show absorption in a low-velocity component (blueshifted by ≈30 km s−1) with lower excitation. The analysis shows that the bands are observed against a continuum with an effective temperature of Tbck ∼ 550 K extinguished with τ6 μmext ∼ 2.5−3 (Ak ∼ 6.9 − 8.3 mag). The high-excitation CO and H2O lines are consistent with v = 0 thermalization with Trot ≈ 450 K and column densities of NCO ≈ (1.7 − 3.5)×1019 cm−2 and NH2O ≈ (1.5 − 3.0)×1019 cm−2. Thermalization of the v = 0 levels of H2O requires either an extreme density of nH2 ≳ 109 cm−3, or radiative excitation by the mid-IR field in a very compact (< 1 pc) optically thick source emitting ∼1010 L⊙. The latter alternative is favored, implying that the observed absorption probes the very early stages of a fully enshrouded active black hole (BH). On the basis of a simple model for BH growth and applying a lifetime constraint to the s2 core, an intermediate-mass BH (IMBH, MBH ∼ 4.5 × 104 M⊙) accreting at super-Eddington rates is suggested, where the observed feedback has not yet been able to break through the natal cocoon.

Vertical Structure of Gas and Dust in Four Debris Disks

We present high-spectral-resolution M-band spectra from iSHELL on NASA’s Infrared Telescope Facility along the line of sight to the debris disk host star HD 32297. We also present a Gemini Planet Imager H-band polarimetric image of the HD 131488 debris disk. We search for fundamental CO absorption lines in the iSHELL spectra of HD 32297, but do not detect any. We place an upper limit on the CO column density of ∼6 × 1015 cm−2. By combining the column density upper limit, the CO mass measured with the Atacama Large Millimeter/submillimeter Array (ALMA), and the geometrical properties of the disk, we estimate the scale height of the CO to be ≲2 au across the radial extent of the disk (∼80–120 au). We use the same method to estimate the CO scale height of three other edge-on, CO-rich debris disks that all have CO observed in absorption with the Hubble Space Telescope as well as in emission with ALMA: β Pictoris, HD 110058, and HD 131488. We compare our estimated CO scale heights of these four systems to the millimeter dust scale heights and find that, under the assumption of hydrostatic equilibrium, there is a potential correlation between the CO and millimeter dust scale heights. There are multiple factors that affect the gas vertical structure such as turbulence, photodissociation with weak vertical mixing, as well as where the gas originates. One possible explanation for the potential correlation could be that the gas and dust are of a similar secondary origin in these four systems.

JWST Reveals Widespread CO Ice and Gas Absorption in the Galactic Center Cloud G0.253+0.016

We report JWST NIRCam observations of G0.253+0.016, the molecular cloud in the Central Molecular Zone known as “The Brick,” with the F182M, F187N, F212N, F410M, F405N, and F466N filters. We catalog 56,146 stars detected in all six filters using the crowdsource package. Stars within and behind The Brick exhibit prodigious absorption in the F466N filter that is produced by a combination of CO ice and gas. In support of this conclusion, and as a general resource, we present models of CO gas and ice and CO2 ice in the F466N, F470N, and F410M filters. Both CO gas and ice contribute to the observed stellar colors. We show, however, that CO gas does not absorb the Pfβ and Huϵ lines in F466N, but that these lines show excess absorption, indicating that CO ice is present and contributes to observed F466N absorption. The most strongly absorbed stars in F466N are extincted by ∼2 mag, corresponding to >80% flux loss. This high observed absorption requires very high column densities of CO, and thus a total CO column that is in tension with standard CO abundance and/or gas-to-dust ratios. This result suggests the CO/H2 ratio and dust-to-gas ratio are greater in the Galactic Center than in the Galactic disk. Ice and/or gas absorption is observed even in the cloud outskirts, implying that additional caution is needed when interpreting stellar photometry in filters that overlap with ice bands throughout galactic centers.

A Rotating Accretion Disk around MWC 297, a Young B1.5Ve Star

High-resolution spectra with iSHELL on IRTF in the K and M bands of the young, heavily accreting B1.5Ve star MWC 297 show numerous double-peaked CO lines. These CO lines originate in an inclined gaseous disk in Keplerian rotation. MWC 297 is the only early B star known to show a Keplerian disk in CO. Analysis of the spectra shows that 12CO 1–0 is optically thick for the low excitation lines. Even the 13CO 1–0 and 12CO 2–1 have somewhat optically thick lines at low J levels. We find that the CO emission in the disk can be fitted with CO being in a narrow ring at a radius of 12 au, with a temperature of 1500 K and a CO column density of 1.6 × 1018 cm−2. This model underestimates the line strength of high-J lines, indicating that they are excited by fluorescence. The CO overtone lines have a similar temperature. The 13CO lines are much brighter than expected from interstellar isotope ratios. The 13CO lines are wider than the 12CO ones, suggesting different excitation conditions. The same is true for 12CO 2–1. We see strong absorption in 12CO and 13CO 1–0 at low J levels, which is due to two cold foreground clouds. These clouds, one with a temperature of 8.3 K and a column density of 6.7 1017 cm−2 and the other one colder and with lower column density, can fully account for the observed extinction toward MWC 297.

Spatiotemporal distribution of air pollutants during a heat wave-induced forest fire event in Uttarakhand.

Prevailing dry conditions and rainfall deficit during the spring season in North India led to heat wave conditions which resulted in widespread and intense forest fire events in the Himalayan state of Uttarakhand during April 16-30, 2022. A total of 7589 active fires were detected by VIIRS during the second half of April 2022 compared to 1558 during the first half. The TROPOMI observed total column values of CO and NO2 increased by 4.4% and 11.7%, respectively during April 16-30, 2022 with respect to April 1-15, 2022. A noticeable increase in surface level concentration of trace gases was also observed at Dehradun. In situ measurements of CO, NOx, and O3 during April 16-30, 2022 show an increase of 133, 35, and 6%compared to previous year observations during the same period. Weather Research and Forecasting model with chemistry (WRF-Chem) is utilized to quantitatively estimate the contribution of this event on the distribution of air pollutants over this state. The model results were evaluated against ERA5 reanalysis, upper air soundings, and TROPOMI-retrieved total column density (TCD) of CO, NO2, and O3. Two simulations with (Fire) and without (NoFire) biomassburning emissions input were performed to quantify the contribution of forest fires to the concentration of trace gases and particulates. The CO, NO2, and O3 emitted/produced from forest fire over Uttarakhand during April 2022 contributed approximately 39.95, 35.73, and 9.97% to the surface concentration of respective gas. In the case of aerosols, it was around 71.20, 71.44, and 33.62% for PM2.5, PM10, and BCrespectively. The vertical profile analysis of pollutants revealed that extreme forest fire events can perturb the distribution of air pollutants from the surface up to 450 hPa.

The CO(6−5) Protostellar Jet Driven by HH 212 as Imaged by ALMA-Band 9

During the earliest phases of the formation of a Sun-like star, jets play a key role in removing angular momentum from the accretion disk, allowing accretion onto the central star. Here we report, for the first time, the bipolar fast (∼150 km s−1) jet driven by the archetypical HH 212, as imaged on scales of ∼100 au by CO(6–5). Adopting Local Thermodynamic Equilibrium conditions, and optically thin emission, we derive the total CO column densities in the jet lobes: ∼1017 cm−2. Both jet lobes are 1200 au long. The observations have been performed using the ALMA Early Science operations, showing how the ALMA interferometer can be efficiently used in its Band 9 (>600 GHz) to investigate of protostellar jet.

ALMACAL – X. Constraints on molecular gas in the low-redshift circumgalactic medium

ABSTRACT Despite its crucial role in galaxy evolution, the complex circumgalactic medium (CGM) remains underexplored. Although, it is known to be multiphase, the importance of the molecular gas phase to the total CGM mass budget is, to date, unconstrained. We present the first constraints on the molecular gas covering fraction in the CGM of low-redshift galaxies, using measurements of CO column densities along sightlines towards mm-bright background quasars with intervening galaxies. We do not detect molecular absorption against the background quasars. For the individual, low-redshift, ‘normal’ galaxy haloes probed here, we can therefore rule out the presence of an extremely molecular gas-rich CGM, as recently reported in high-redshift protoclusters and around luminous active galactic nuclei. We also set statistical limits on the volume filling factor of molecular material in the CGM as a whole, and as a function of radius. Interstellar medium-like molecular clouds of ∼30 pc in radius with column densities of N(CO) ≳ 1016 cm−2 have volume filling factors of <0.2 per cent. Large-scale smooth gas reservoirs are ruled out much more stringently. The development of this technique in the future will allow deeper constraining limits to be set on the importance (or unimportance) of molecular gas in the CGM.

AP-HTPB propellant combustion under strain conditions with laser absorption spectroscopy.

Understanding the combustion behaviors of solid propellant with different levels of strains is of practical interest. In this work, an experimental study of the effects of static and dynamic strains on the burning rate, temperature, CO, and C O 2 formation of aluminized ammonium perchlorate (AP)-hydroxyl terminated poly-butadiene (HTPB) propellant combustion was presented at initial pressures of 0.1MPa, 0.2MPa, and 0.5MPa. The strains were being applied onto solid propellant by exerting static and cyclic loadings. The propellant burning rate was acquired by a 4kHz high-speed photography system, and the combustion temperature, CO, and C O 2 column densities were measured at 10kHz through laser absorption spectroscopy (LAS). At atmospheric pressure, it was demonstrated that the propellant burning rate increased with tensile stress and decreased with compressive stress. The measured flame temperature showed a similar correlation with strains as compared to the propellant burning rate. At elevated pressures, the increase of the propellant burning rate due to tensile stress was more evident, while the difference in combustion temperatures was less significant. For the cyclic strain condition, the variations of the measured C O 2 and CO column densities were consistent with the static strain condition.

LIDA: The Leiden Ice Database for Astrochemistry

Context. High-quality vibrational spectra of solid-phase molecules in ice mixtures and for temperatures of astrophysical relevance are needed to interpret infrared observations toward protostars and background stars. Such data are collected worldwide by several laboratory groups in support of existing and upcoming astronomical observations. Over the last 25 yr, the Laboratory for Astrophysics at Leiden Observatory has provided more than 1100 (high-resolution) spectra of diverse ice samples. Aims. In time with the recent launch of the James Webb Space Telescope, we have fully upgraded the Leiden Ice Database for Astrochemistry (LIDA) adding recently measured spectra. The goal of this paper is to describe what options exist regarding accessing and working with a large collection of infrared (IR) spectra, and the ultraviolet-visible (UV/vis) to the mid-infrared refractive index of H2O ice. This also includes astronomy-oriented online tools to support the interpretation of IR ice observations. Methods. This ice database is based on open-source Python software, such as Flask and Bokeh, used to generate the web pages and graph visualization, respectively. Structured Query Language (SQL) is used for searching ice analogs within the database and Jmol allows for three-dimensional molecule visualization. The database provides the vibrational modes of molecules known and expected to exist as ice in space. These modes are characterized using density functional theory with the orca software. The IR data in the database are recorded via transmission spectroscopy of ice films condensed on cryogenic substrates. The real UV/vis refractive indices of H2O ice are derived from interference fringes created from the simultaneous use of a monochromatic HeNe laser beam and a broadband Xe-arc lamp, whereas the real and imaginary mid-IR values are theoretically calculated. LIDA not only provides information on fundamental ice properties, but it also offers online tools. The first tool, SPECFY, is directly linked to the data in the database to create a synthetic spectrum of ices towards protostars. The second tool allows the uploading of external files and the calculation of mid-infrared refractive index values. Results. LIDA provides an open-access and user-friendly platform to search, download, and visualize experimental data of astrophysically relevant molecules in the solid phase. It also provides the means to support astronomical observations; in particular, those that will be obtained with the James Webb Space Telescope. As an example, we analysed the Infrared Space Observatory spectrum of the protostar AFGL 989 using the resources available in LIDA and derived the column densities of H2O, CO and CO2 ices.

Different Degrees of Nitrogen and Carbon Depletion in the Warm Molecular Layers of Protoplanetary Disks

Observations have revealed that the elemental abundances of carbon and oxygen in the warm molecular layers of some protoplanetary disks are depleted compared to those in the interstellar medium by a factor of ∼10–100. Meanwhile, little is known about nitrogen. To investigate the time evolution of nitrogen, carbon, and oxygen elemental abundances in disks, we develop a one-dimensional plane-parallel model that incorporates dust settling, turbulent diffusion of dust and ices, as well as gas-ice chemistry including the chemistry driven by stellar UV/X-rays and galactic cosmic rays. We find that gaseous CO in the warm molecular layer is converted to CO2 ice and locked up near the midplane via the combination of turbulent mixing (i.e., the vertical cold finger effect) and ice chemistry driven by stellar UV photons. On the other hand, gaseous N2, the main nitrogen reservoir in the warm molecular layer, is less processed by ice chemistry and exists as it is. Then, nitrogen depletion occurs solely through the vertical cold finger effect of N2. As the binding energy of N2 is lower than that of CO and CO2, the degree of nitrogen depletion is smaller than that of carbon and oxygen depletion, leading to higher elemental abundance of nitrogen than that of carbon and oxygen. This evolution occurs within 1 Myr and proceeds further, when the α parameter for the diffusion coefficient is ≳10−3. Consequently, the N2H+/CO column density ratio increases with time. How the vertical transport affects the midplane ice composition is briefly discussed.

The Star Formation–Gas Density Relation in Four Galactic GMCs: Effects of Stellar Feedback

We present maps of four Galactic giant molecular clouds (GMCs) in the J = 2-1 emission of both CO and 13CO. We use an LTE analysis to derive maps of the CO excitation temperature and column density and the distribution of total molecular gas column density, Σgas. The depletion of CO by freeze-out onto cold dust grains is accounted for by an approximation to the results of Lewis et al., which were derived from far-IR observations with Herschel. The surface density of young stellar objects (YSOs) is obtained from published catalogs. The mean YSO surface density exhibits a power-law dependence on Σgas, with exponents in the range 0.9–1.9. Gas column density probability distribution functions show power-law tails extending to high column densities. The distributions of sonic Mach number, M S , are sharply peaked at MS ∼ 5–8 for 3 GMCs; a fourth has a broad distribution up to M S = 30, possibly a result of feedback effects from multiple OB stars. An analysis following the methodology of Pokhrel et al. finds that our sample of GMCs shows power-law relations that are somewhat shallower than those found by Pokhrel et al. for the star formation rate versus 〈Σgas〉 and versus 〈Σgas〉/t ff in a different sample of clouds. We discuss possible differences in the two samples of star-forming clouds and the effects of stellar feedback on the relation between gas density and star formation rate.

Influence of Magnetic Fields on Precipitating Solar Wind Hydrogen at Mars.

Solar wind protons can interact directly with the hydrogen corona of Mars through charge exchange, resulting in energetic neutral atoms (ENAs) able to penetrate deep into the upper atmosphere of Mars. ENAs can undergo multiple charge changing interactions, leading to an observable beam of penetrating protons in the upper atmosphere. We seek to characterize the behavior of these protons in the presence of magnetic fields using data collected by the Mars Atmosphere and Volatile EvolutioN spacecraft. We find that backscattered penetrating proton flux is enhanced in regions where the magnetic field strength is greater than 200 nT. We also find a strong correlation at CO2 column densities less than 5.5 × 1014 cm−2 between magnetic field strength and the observed backscattered and downward flux. We do not see significant changes in penetrating proton flux with magnetic field strengths on the order of 10 nT.

Molecular Gas Properties and CO-to-H2 Conversion Factors in the Central Kiloparsec of NGC 3351

Abstract The CO-to-H2 conversion factor (α CO) is critical to studying molecular gas and star formation in galaxies. The value of α CO has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on ~kiloparsec scales revealed low values of α CO in the centers of some barred spiral galaxies, including NGC 3351. We present new Atacama Large Millimeter/submillimeter Array Band 3, 6, and 7 observations of 12CO, 13CO, and C18O lines on 100 pc scales in the inner ∼2 kpc of NGC 3351. Using multiline radiative transfer modeling and a Bayesian likelihood analysis, we infer the H2 density, kinetic temperature, CO column density per line width, and CO isotopologue abundances on a pixel-by-pixel basis. Our modeling implies the existence of a dominant gas component with a density of 2–3 × 103 cm−3 in the central ∼1 kpc and a high temperature of 30–60 K near the nucleus and near the contact points that connect to the bar-driven inflows. Assuming a CO/H2 abundance of 3 × 10−4, our analysis yields α CO ∼ 0.5–2.0 M ⊙ (K km s−1 pc2)−1 with a decreasing trend with galactocentric radius in the central ∼1 kpc. The inflows show a substantially lower α CO ≲ 0.1 M ⊙ (K km s−1 pc2)−1, likely due to lower optical depths caused by turbulence or shear in the inflows. Over the whole region, this gives an intensity-weighted α CO of ∼1.5 M ⊙ (K km s−1 pc2)−1, which is similar to previous dust-modeling-based results at kiloparsec scales. This suggests that low α CO on kiloparsec scales in the centers of some barred galaxies may be due to the contribution of low-optical-depth CO emission in bar-driven inflows.

Synthetic CO emission and the XCO factor of young molecular clouds: a convergence study

ABSTRACT The properties of synthetic CO emission from 3D simulations of forming molecular clouds are studied within the SILCC-Zoom project. Since the time-scales of cloud evolution and molecule formation are comparable, the simulations include a live chemical network. Two sets of simulations with an increasing spatial resolution (dx = 3.9 pc to dx = 0.06 pc) are used to investigate the convergence of the synthetic CO emission, which is computed by post-processing the simulation data with the radmc-3d radiative transfer code. To determine the excitation conditions, it is necessary to include atomic hydrogen and helium alongside H2, which increases the resulting CO emission by ∼7–26 per cent. Combining the brightness temperature of 12CO and 13CO, we compare different methods to estimate the excitation temperature, the optical depth of the CO line and hence, the CO column density. An intensity-weighted average excitation temperature results in the most accurate estimate of the total CO mass. When the pixel-based excitation temperature is used to calculate the CO mass, it is over-/underestimated at low/high CO column densities where the assumption that 12CO is optically thick while 13CO is optically thin is not valid. Further, in order to obtain a converged total CO luminosity and hence 〈XCO〉 factor, the 3D simulation must have dx ≲ 0.1 pc. The 〈XCO〉 evolves over time and differs for the two clouds; yet pronounced differences with numerical resolution are found. Since high column density regions with a visual extinction larger than 3 mag are not resolved for dx ≳ 1 pc, in this case the H2 mass and CO luminosity both differ significantly from the higher resolution results and the local XCO is subject to strong noise. Our calculations suggest that synthetic CO emission maps are only converged for simulations with dx ≲ 0.1 pc.

Atmospheric CO[formula omitted] retrieval from satellite spectral measurements by a two-step machine learning approach

In the present study, a two-step machine learning approach is proposed and tested based on the Greenhouse gases Observing SATllite (GOSAT) spectral measurements to efficiently retrieve the atmospheric CO2 column density. A simple one-dimensional line-by-line forward radiation model is developed to simulate the GOSAT observed spectra, which is used to generate the training data for the two-step machine learning model. By using the two-step machine learning model, the atmospheric spectral optical thickness can be retrieved first, and then the CO2 column density is retrieved from the optical thickness spectrum later. As a proof-of-concept study, the model has been applied to retrieve CO2 column density in Australia for clear sky conditions from the GOSAT observed 1.6 μm spectra. Results have shown that the accuracy of the proposed method for CO2 column density retrieval is about 3 ppm but the efficiency is found to be excellent. Further improvements for improving the retrieval accuracy are discussed.

Molecules with ALMA at Planet-forming Scales (MAPS). IV. Emission Surfaces and Vertical Distribution of Molecules

Abstract The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a unique opportunity to study the vertical distribution of gas, chemistry, and temperature in the protoplanetary disks around IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. By using the asymmetry of molecular line emission relative to the disk major axis, we infer the emission height (z) above the midplane as a function of radius (r). Using this method, we measure emitting surfaces for a suite of CO isotopologues, HCN, and C2H. We find that 12CO emission traces the most elevated regions with z / r > 0.3 , while emission from the less abundant 13CO and C18O probes deeper into the disk at altitudes of z / r ≲ 0.2 . C2H and HCN have lower opacities and signal-to-noise ratios, making surface fitting more difficult, and could only be reliably constrained in AS 209, HD 163296, and MWC 480, with z / r ≲ 0.1 , i.e., relatively close to the planet-forming midplanes. We determine peak brightness temperatures of the optically thick CO isotopologues and use these to trace 2D disk temperature structures. Several CO temperature profiles and emission surfaces show dips in temperature or vertical height, some of which are associated with gaps and rings in line and/or continuum emission. These substructures may be due to local changes in CO column density, gas surface density, or gas temperatures, and detailed thermochemical models are necessary to better constrain their origins and relate the chemical compositions of elevated disk layers with those of planet-forming material in disk midplanes. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.

Molecules with ALMA at Planet-forming Scales (MAPS). VIII. CO Gap in AS 209—Gas Depletion or Chemical Processing?

Abstract Emission substructures in gas and dust are common in protoplanetary disks. Such substructures can be linked to planet formation or planets themselves. We explore the observed gas substructures in AS 209 using thermochemical modeling with RAC2D and high-spatial-resolution data from the Molecules with ALMA at Planet-forming Scales (MAPS) program. The observations of C18O J = 2–1 emission exhibit a strong depression at 88 au overlapping with the positions of multiple gaps in millimeter dust continuum emission. We find that the observed CO column density is consistent with either gas surface-density perturbations or chemical processing, while C2H column density traces changes in the C/O ratio rather than the H2 gas surface density. However, the presence of a massive planet (>0.2 M Jup) would be required to account for this level of gas depression, which conflicts with constraints set by the dust emission and the pressure profile measured by gas kinematics. Based on our models, we infer that a local decrease of CO abundance is required to explain the observed structure in CO, dominating over a possible gap-carving planet present and its effect on the H2 surface density. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.

Molecules with ALMA at Planet-forming Scales (MAPS). VII. Substellar O/H and C/H and Superstellar C/O in Planet-feeding Gas

Abstract The elemental composition of the gas and dust in a protoplanetary disk influences the compositions of the planets that form in it. We use the Molecules with ALMA at Planet-forming Scales (MAPS) data to constrain the elemental composition of the gas at the locations of potentially forming planets. The elemental abundances are inferred by comparing source-specific gas-grain thermochemical models with variable C/O ratios and small-grain abundances from the DALI code with CO and C2H column densities derived from the high-resolution observations of the disks of AS 209, HD 163296, and MWC 480. Elevated C/O ratios (∼2.0), even within the CO ice line, are necessary to match the inferred C2H column densities over most of the pebble disk. Combined with constraints on the CO abundances in these systems, this implies that both the O/H and C/H ratios in the gas are substellar by a factor of 4–10, with the O/H depleted by a factor of 20–50, resulting in the high C/O ratios. This necessitates that even within the CO ice line, most of the volatile carbon and oxygen is still trapped on grains in the midplane. Planets accreting gas in the gaps of the AS 209, HD 163296, and MWC 480 disks will thus acquire very little carbon and oxygen after reaching the pebble isolation mass. In the absence of atmosphere-enriching events, these planets would thus have a strongly substellar O/H and C/H and superstellar C/O atmospheric composition. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.