The evolution of protoplanetary disks is dominated by the conservation of angular momentum, where the accretion of material onto the central star is driven by viscous expansion of the outer disk or by disk winds extracting angular momentum without changing the disk size. Studying the time evolution of disk sizes allows us therefore to distinguish between viscous stresses or disk winds as the main mechanism of disk evolution. Observationally, estimates of the disk gaseous outer radius are based on the extent of the CO rotational emission, which, during the evolution, is also affected by the changing physical and chemical conditions in the disk. We use physical-chemical DALI models to study how the extent of the CO emission changes with time in a viscously expanding disk and investigate to what degree this observable gas outer radius is a suitable tracer of viscous spreading and whether current observations are consistent with viscous evolution. We find that the gas outer radius (R_co) measured from our models matches the expectations of a viscously spreading disk: R_co increases with time and for a given time R_co is larger for a disk with a higher viscosity alpha_visc. However, in the extreme case where the disk mass is low (less than 10^-4 Msun) and alpha_visc is high (larger than 10^-2), R_co will instead decrease with time as a result of CO photodissociation in the outer disk. For most disk ages R_co is up to 12x larger than the characteristic size R_c of the disk, and R_co/R_c is largest for the most massive disk. As a result of this difference, a simple conversion of R_co to alpha_visc will overestimate the true alpha_visc of the disk by up to an order of magnitude. We find that most observed gas outer radii in Lupus can be explained using a viscously evolving disk that starts out small (R_c = 10 AU) and has a low viscosity (alpha_visc = 10^-4 - 10^-3).

Water is a key volatile that provides insights into the initial stages of planet formation. The low water abundances inferred from water observations toward low-mass protostellar objects may point to a rapid locking of water as ice by large dust grains during star and planet formation. However, little is known about the water vapor abundance in newly formed planet-forming disks. We aim to determine the water abundance in embedded Keplerian disks through spatially-resolved observations of H$_2^{18}$O lines to understand the evolution of water during star and planet formation. We present H$_2^{18}$O line observations with ALMA and NOEMA millimeter interferometers toward five young stellar objects. NOEMA observed the 3$_{1,3}$ - $2_{2,0}$ line (E$_{\rm up}$ = 203.7 K) while ALMA targeted the $4_{1,4}$ - $3_{2,1}$ line (E$_{\rm up}$ = 322.0 K). Water column densities are derived considering optically thin and thermalized emission. Our observations are sensitive to the emission from the known Keplerian disks around three out of the five Class I objects in the sample. No H$_2^{18}$O emission is detected toward any of our five Class I disks. We report upper limits to the integrated line intensities. The inferred water column densities in Class I disks are N < 10$^{15}$ cm$^{-2}$ on 100 au scales which include both disk and envelope. The upper limits imply a disk-averaged water abundance of $\lesssim 10^{-6}$ with respect to H$_2$ for Class I objects. After taking into account the physical structure of the disk, the upper limit to the water abundance averaged over the inner warm disk with $T>$ 100 K is between 10$^{-7}$ up to 10$^{-5}$. Water vapor is not abundant in warm protostellar envelopes around Class I protostars. Upper limits to the water vapor column densities in Class I disks are at least two orders magnitude lower than values found in Class 0 disk-like structures.

The atmospheric composition of giant planets carries the information of their formation history. Superstellar C/H ratios are seen in atmospheres of Jupiter, Saturn, and various giant exoplanets. Also, giant exoplanets show a wide range of C/O ratio. To explain these ratios, one hypothesis is that protoplanets accrete carbon-enriched gas when a large number of icy pebbles drift across the CO snowline. Here we report the first direct evidence of an elevated C/H ratio in disk gas. We use two thermo-chemical codes to model the $^{13}$C$^{18}$O, C$^{17}$O, and C$^{18}$O (2-1) line spectra of the HD 163296 disk. We show that the gas inside the CO snowline ($\sim$70 au) has a C/H ratio of 1-2 times higher than the stellar value. This ratio exceeds the expected value substantially, as only 25-60% of the carbon should be in gas at these radii. Although we cannot rule out the case of a normal C/H ratio inside 70 au, the most probable solution is an elevated C/H ratio of 2-8 times higher than the expectation. Our model also shows that the gas outside 70 au has a C/H ratio of 0.1$\times$ the stellar value. This picture of enriched C/H gas at the inner region and depleted gas at the outer region is consistent with numerical simulations of icy pebble growth and drift in protoplanetary disks. Our results demonstrate that the large-scale drift of icy pebble can occur in disks and may significantly change the disk gas composition for planet formation.

Recent theoretical, numerical, and observational work have suggested that when a growing planet opens a gap in its disk the flow of gas into the gap is dominated by gas falling vertically from a height of at least one gas scale height. Our primary objective is to include, for the first time, the chemical impact that accreting gas above the midplane will have on the resulting C/O. We compute the accretion of gas onto planetary cores beginning at different disk radii and track the chemical composition of the gas and small icy grains to predict the resulting carbon-to-oxygen ratio (C/O) in their atmospheres. In our model, all of the planets which began their evolution inward of 60 AU open a gap in the gas disk, and hence are chemically affected by the vertically accreting gas. Two important conclusions follow from this vertical flow: (1) more oxygen rich icy dust grains become available for accretion onto the planetary atmosphere. (2) The chemical composition of the gas dominates the final C/O of planets in the inner ($<$ 20 AU) part of the disk. This implies that with the launch of the James Webb Space Telescope we can trace the disk material that sets the chemical composition of exoplanetary atmospheres.

Context. The gas kinetic temperature (TK) determines the physical and chemical evolution of the Interestellar Medium (ISM). However, obtaining reliable TK estimates usually requires expensive observations including the combination of multi-line analysis and dedicated radiative transfer calculations. Aims. This work explores the use of HCN and HNC observations, and particularly its I(HCN)/I(HNC) intensity ratio of their J=1-0 lines, as direct probe of the gas kinetic temperature in the molecular ISM. Methods. We obtained a new set of large-scale observations of both HCN and HNC (1-0) lines along the Integral Shape Filament (ISF) in Orion. In combination with ancillary gas and dust temperature measurements, we find a systematic temperature dependence of the observed I(HCN)/I(HNC) intensity ratio across our maps. Additional comparisons with chemical models demonstrate that these observed I(HCN)/I(HNC) variations are driven by the effective destruction and isomerization mechanisms of HNC under low energy barriers. Results. The observed variations of I(HCN)/I(HNC) with TK can be described with a two-part linear function. This empirical calibration is then used to create a temperature map of the entire ISF. Comparisons with similar dust temperature measurements in this cloud, as well as in other regions and galactic surveys, validate this simple technique to obtain direct estimates of the gas kinetic temperature in a wide range of physical conditions and scales with an optimal working range between 15 K and 40 K. Conclusions. Both observations and models demonstrate the strong sensitivity of the I(HCN)/I(HNC) ratio to the gas kinetic temperature. Since these lines are easily obtained in observations of local and extragalactic sources, our results highlight the potential use of this observable as new chemical thermometer for the ISM.

The region around the H$_2$O ice line, due to its higher surface density, seems to be the ideal location to form planets. The core of Jupiter, as well as the cores of close in gas giants are thus thought to form in this region of the disk. Actually constraining the formation location of individual planets has proven to be difficult, however. We aim to use the Nitrogen abundance in Jupiter, which is around 4 times solar, in combination with \textit{Juno} constraints on the total mass of heavy elements in Jupiter, to narrow down its formation scenario. Different pathways of enrichment of Jupiter's atmosphere, such as the accretion of enriched gas, pebbles or planetesimals are considered and their implications for the oxygen abundance of Jupiter is discussed. The super solar Nitrogen abundance in Jupiter necessitates the accretion of extra N$_2$ from the proto-solar nebula. The only location of the disk that this can happen is outside, or just inside the N$_2$ ice line. These constraints favor a pebble accretion origin of Jupiter, both from the composition as well as from a planet formation perspective. We predict that Jupiter's oxygen abundance is between 3.6 and 4.5 times solar.

Gas giants accrete their envelopes from the gas and dust of proto-planetary disks, so it is important to determine the composition of the inner few AU, where most giant planets are expected to form. We aim to constrain the elemental carbon and oxygen abundance in the inner disk ($R<$2.3 AU) of TW Hya and compare with the outer disk ($R>2.3$ AU) where carbon and oxygen appear underabundant by a factor of $\sim$50. Archival infrared observations of TW Hya are compared with a detailed thermo-chemical model, DALI. The inner disk gas mass and elemental C and O abundances are varied to fit the infrared CO, H$_2$ and H$_2$O line fluxes. Best fitting models have an inner disk that has a gas mass of $ 2 \times 10^{-4} M_\odot$ with C/H $\approx 3 \times 10^{-6}$ and O/H $\approx 6 \times 10^{-6}$. The elemental oxygen and carbon abundances of the inner disk are $\sim 50$ times underabundant compared to the ISM and are consistent with those found in the outer disk. The uniformly low volatile abundances imply that the inner disk is not enriched by ices on drifting bodies that evaporate. This indicates that drifting grains are stopped in a dust trap outside the water ice line. Such a dust trap would also form a cavity as seen in high resolution sub-millimeter continuum observations. If CO is the major carbon carrier in the ices, dust needs to be trapped efficiently outside the CO ice line of $\sim$20 AU. This would imply that the shallow sub-millimeter rings in the TW Hya disk outside of 20 AU correspond to very efficient dust traps. The more likely scenario is that more than 98\% of the CO has been converted into less volatile species, e.g. CO$_2$ and CH$_3$OH. A giant planet forming in the inner disk would be accreting gas with low carbon and oxygen abundances as well as very little icy dust, potentially leading to a planet atmosphere with strongly substellar C/H and O/H ratios.

[abridged]CO rovibrational lines are efficient probes of warm molecular gas and can give unique insights into the inner 10 AU of proto-planetary disks. Recent studies have found a relation between the ratio of lines originating from the second and first vibrationally excited state, denoted as $v2/v1$, and the emitting radius of CO. In disks around Herbig Ae stars the vibrational excitation is low when CO lines come from close to the star, and high when lines only probe gas at large radii (more than 5 AU). We aim to find explanations for the observed trends between CO vibrational ratio, emitting radii, and NIR excess, and identify their implications in terms of the physical and chemical structure of inner disks around Herbig stars. Slab models and full disk thermo chemical models are calculated. Simulated observations from the models are directly compared to the data. Broad CO lines with low vibrational ratios are best explained by a warm (400-1300 K) inner disk surface with gas-to-dust ratios below 1000; no CO is detected within/at the inner dust rim, due to dissociation at high temperatures. In contrast, explaining the narrow lines with high vibrational ratios requires an inner cavity of a least 5 AU in both dust and gas, followed by a cool (100-300 K) molecular gas reservoir with gas-to-dust ratios greater than 10000 at the cavity wall. In all cases the CO gas must be close to thermalization with the dust. The high gas-to-dust ratios needed to explain high $v2/v1$ in narrow CO lines for a subset of group I disks can naturally be interpreted as due to the dust traps that have been proposed to explain millimeter dust cavities. The broad lines seen in most group II objects indicate a very flat disk in addition to inner disk substructures within 10 AU that can be related to the substructures recently observed with ALMA.

CO has long been thought to be the best tracer to measure gas masses as it is readily detected at (sub)mm wavelengths in many disks. Inferred gas masses from CO in recent ALMA observations of large samples of disks seem inconsistent with their inferred dust masses. The derived gas-to-dust mass ratios from CO are 1-2 orders of magnitude lower than $\sim$100 even if photodissociation and freeze-out are included. Herschel measurements of HD line emission imply gas masses in line with gas-to-dust mass ratios of 100. This suggests that at least one additional mechanism is removing gaseous CO. Here we test the suggestion that the bulk of the CO is chemically processed and that the carbon is sequestered into less volatile species such as CO2, CH3OH and CH4 in the dense, shielded midplane regions of the disk. Using our gas-grain chemical code we perform a parameter exploration and follow the CO abundance evolution over a range of conditions representative of shielded disk midplanes. We find that no chemical processing of CO takes place on 1-3 Myr timescales for low cosmic-ray ionisation rates, $< 5\times 10^{-18}$ s$^{-1}$. Assuming an ionisation rate of $10^{-17}$ s$^{-1}$, more than 90% of the CO is destroyed, but only in the parts of the disk below 30 K. This order of magnitude destruction of CO is robust against the choice of grain-surface reaction rate parameters except for temperatures below 20 K, where there is a strong dependence on the H tunnelling barrier. The low temperatures needed for CO chemical processing indicate that the exact disk temperature structure is important, with warm disks around luminous Herbig stars expected to have little to no CO conversion. For cold disks around T Tauri stars, a large fraction of the emitting CO layer is affected unless the disks are young (< 1 Myr). This can lead to inferred gas masses that are up to two orders of magnitude too low.

The efficiency of radial transport of icy solid material from outer disk to the inner disk is currently unconstrained. Efficient radial transport of icy dust grains could significantly alter the composition of the gas in the inner disk. Our aim is to model the gaseous CO$_2$ abundance in the inner disk and use this to probe the efficiency of icy dust transport in a viscous disk. Features in the simulated CO$_2$ spectra are investigated for their dust flux tracing potential. We have developed a 1D viscous disk model that includes gas and grain motions as well as dust growth, sublimation and freeze-out and a parametrisation of the CO$_2$ chemistry. The thermo-chemical code DALI was used to model the mid-infrared spectrum of CO$_2$, as can be observed with JWST-MIRI. CO$_2$ ice sublimating at the iceline increases the gaseous CO$_2$ abundance to levels equal to the CO$_2$ ice abundance of $\sim 10^{-5}$, which is three orders of magnitude more than the gaseous CO$_2$ abundances of $\sim 10^{-8}$ observed by Spitzer. Grain growth and radial drift further increase the gaseous CO$_2$ abundance. A CO$_2$ destruction rate of at least $10^{-11}$ s$^{-1}$ is needed to reconcile model prediction with observations. This rate is at least two orders of magnitude higher than the fastest known chemical destruction rate. A range of potential physical mechanisms to explain the low observed CO$_2$ abundances are discussed. Transport processes in disks can have profound effects on the abundances of species in the inner disk. The discrepancy between our model and observations either suggests frequent shocks in the inner 10 AU that destroy CO$_2$, or that the abundant midplane CO$_2$ is hidden from our view by an optically thick column of low abundance CO$_2$ in to the disk surface XDR/PDR. Other molecules, such as CH$_4$ or NH$_3$, can give further handles on the rate of mass transport.

[Abridged] The infrared ro-vibrational emission lines from organic molecules in the inner regions of protoplanetary disks are unique probes of the physical and chemical structure of planet forming regions and the processes that shape them. The non-LTE excitation effects of carbon dioxide (CO2) are studied in a full disk model to evaluate: (i) what the emitting regions of the different CO2 ro-vibrational bands are; (ii) how the CO2 abundance can be best traced using CO2 ro-vibrational lines using future JWST data and; (iii) what the excitation and abundances tell us about the inner disk physics and chemistry. CO2 is a major ice component and its abundance can potentially test models with migrating icy pebbles across the iceline. A full non-LTE CO2 excitation model has been built. The characteristics of the model are tested using non-LTE slab models. Subsequently the CO2 line formation has been modelled using a two-dimensional disk model representative of T-Tauri disks. The CO2 gas that emits in the 15 $\mu$m and 4.5 $\mu$m regions of the spectrum is not in LTE and arises in the upper layers of disks, pumped by infrared radiation. The v$_2$ 15 $\mu$m feature is dominated by optically thick emission for most of the models that fit the observations and increases linearly with source luminosity. Its narrowness compared with that of other molecules stems from a combination of the low rotational excitation temperature (~250 K) and the inherently narrower feature for CO2. The inferred CO2 abundances derived for observed disks are more than two orders of magnitude lower than those in interstellar ices (~10$^5$), similar to earlier LTE disk estimates. Line-to-continuum ratios are low, of order a few %, thus high signal-to-noise (S/N > 300) observations are needed for individual line detections. Prospects of accurate abundance retreival with JWST-MIRI and JWST-NIRSpec are discussed.

A new collection of photodissociation and photoionisation cross sections for 102 atoms and molecules of astrochemical interest has been assembled, along with a brief review of the basic processes involved. These have been used to calculate dissociation and ionisation rates, with uncertainties, in a standard ultraviolet interstellar radiation field (ISRF) and wavelength-dependent radiation fields. The new ISRF rates generally agree within 30% with our previous compilations, with a few notable exceptions. The reduction of rates in shielded regions was calculated as a function of dust, molecular and atomic hydrogen, atomic C, and self-shielding column densities. The relative importance of shielding types depends on the species in question and the dust optical properties. The new data are publicly available from the Leiden photodissociation and ionisation database. Sensitivity of rates to variation of temperature and isotope, and cross section uncertainties, are tested. Tests were conducted with an interstellar-cloud chemical model, and find general agreement (within a factor of two) with the previous iteration of the Leiden database for the ISRF, and order-of-magnitude variations assuming various kinds of stellar radiation. The newly parameterised dust-shielding factors makes a factor-of-two difference to many atomic and molecular abundances relative to parameters currently in the UDfA and KIDA astrochemical reaction databases. The newly-calculated cosmic-ray induced photodissociation and ionisation rates differ from current standard values up to a factor of 5. Under high temperature and cosmic-ray-flux conditions the new rates alter the equilibrium abundances of abundant dark cloud abundances by up to a factor of two. The partial cross sections for H2O and NH3 photodissociation forming OH, O, NH2 and NH are also evaluated and lead to radiation-field-dependent branching ratios.