D/H & 12C/13C Ratios in Astrochemistry

• D/H and 12C/13C ratios are fundamental isotopic diagnostics used to gauge deuterium and carbon fractionation in astrochemistry, planetary science, and cosmochemistry.
• They are measured via radio/mm spectroscopy, mass spectrometry, and spatial mapping, enabling detailed insights into molecular formation and volatile evolution.
• Observed variations across environments—from molecular clouds to protostellar disks and comets—constrain chemical pathways and inform models of chemical inheritance and evolution.

The deuterium-to-hydrogen (D/H) and $^{12}$C/$^{13}$C ratios are critical isotopic diagnostics in astrochemistry, planetary science, and cosmochemistry. These ratios probe isotope fractionation processes, molecular formation environments, and the mixing histories of volatiles and organics in diverse astronomical contexts, including molecular clouds, protoplanetary disks, comets, and protostellar environments.

1. Fundamental Definitions and Measurement Approaches

The D/H ratio is defined as

$$\mathrm{D/H} = \frac{N(\mathrm{XD})}{N(\mathrm{XH})}$$

where $N$ represents the column density of a molecule X containing deuterium (D) relative to its hydrogen (H) analogue. For carbon, the $^{12}$C/$^{13}$C ratio is

$${}^{12}\mathrm{C}/{}^{13}\mathrm{C} = \frac{N(\text{species with }^{12}\mathrm{C})}{N(\text{species with }^{13}\mathrm{C})}$$

These ratios are typically derived via rotational or vibrational spectroscopy of isotopologue transitions, or gas-phase/in situ mass spectrometry in Solar System contexts (Biver et al., 2016, Araki et al., 2016, Müller et al., 2022, Jørgensen et al., 2018). Astrophysical observations typically rely on assumptions of local thermodynamic equilibrium (LTE), optically thin emission, and the use of partition functions and excitation temperatures explicitly referenced in the original analyses.

The measurement methodologies include:

• Radio/mm spectroscopy of molecular lines, applying LTE/non-LTE excitation analysis with optically thin approximation, e.g., for HC$_3$N and HCN isotopologues (Araki et al., 2016, Biver et al., 2016, Rampinelli et al., 4 Apr 2025).
• Mass spectrometry, e.g., ROSINA/DFMS on-board Rosetta for direct volatile composition, with rigorous gain and pixel correction protocols for precise D/H and $^{13}$C/$^{12}$C extraction in multiple volatile species (Müller et al., 2022).
• Radial and spatially resolved mapping using interferometric data, e.g., ALMA, to reconstruct isotopic gradients in protoplanetary disks (Rampinelli et al., 4 Apr 2025).

2. Observational Results Across Environments

D/H and $^{13}$0C/$^{13}$1C have been measured in a wide range of solar/galactic environments, revealing substantial diversity and diagnostic value.

Environment	D/H	$^{13}$2C/$^{13}$3C	Reference
ISM, Solar neighborhood (HCO$^{13}$4)	—	$^{13}$5	(Luo et al., 2024)
Dense cores (HC$^{13}$6N, L1527)	$^{13}$7	$^{13}$8	(Araki et al., 2016)
Protoplanetary disk (PDS 70, DCN/HCN)	$^{13}$9	$$\mathrm{D/H} = \frac{N(\mathrm{XD})}{N(\mathrm{XH})}$$0	(Rampinelli et al., 4 Apr 2025)
Protostar IRAS 16293B (CH$$\mathrm{D/H} = \frac{N(\mathrm{XD})}{N(\mathrm{XH})}$$1OCH$$\mathrm{D/H} = \frac{N(\mathrm{XD})}{N(\mathrm{XH})}$$2)	$$\mathrm{D/H} = \frac{N(\mathrm{XD})}{N(\mathrm{XH})}$$3 (4%)	$$\mathrm{D/H} = \frac{N(\mathrm{XD})}{N(\mathrm{XH})}$$4	(Jørgensen et al., 2018)
Comet 67P (H$$\mathrm{D/H} = \frac{N(\mathrm{XD})}{N(\mathrm{XH})}$$5O)	$$\mathrm{D/H} = \frac{N(\mathrm{XD})}{N(\mathrm{XH})}$$6	$$\mathrm{D/H} = \frac{N(\mathrm{XD})}{N(\mathrm{XH})}$$7 (in alkanes)	(Müller et al., 2022)
Comet Lovejoy (H$$\mathrm{D/H} = \frac{N(\mathrm{XD})}{N(\mathrm{XH})}$$8O, HCN)	$$\mathrm{D/H} = \frac{N(\mathrm{XD})}{N(\mathrm{XH})}$$9	$N$0 (HCN)	(Biver et al., 2016)
Comet Lemmon (H$N$1O)	$N$2	$N$3 (HCN)	(Biver et al., 2016)

These ratios often deviate substantially from elemental abundances (e.g., ISM D/H $N$4, $N$5C/$N$6C $N$7 69–70), revealing enrichment or dilution as a result of fractionation, local chemical pathways, and evolutionary conditions.

3. Physical and Chemical Fractionation Mechanisms

Deuterium Fractionation Processes

Enhanced D/H in molecules reflects low-temperature ion-molecule exchange reactions and grain-surface chemistry. Specifically:

• In dense cold clouds (T $N$8 30 K), D-enrichment occurs via H$N$9 + HD $^{12}$0 H$^{12}$1D$^{12}$2 + H$^{12}$3 and subsequent transfer of deuterium to complex molecules (e.g., DCN, DCCCN), often resulting in D/H ratios orders of magnitude above elemental values (Araki et al., 2016, Rampinelli et al., 4 Apr 2025).
• In cometary and protostellar ices, D/H in organics is further enhanced by grain-surface reactions during mantle warm-up, especially for complex molecules (e.g., $^{12}$4\% for CH$^{12}$5OH, HCOOH, CH$^{12}$6OCHO) (Jørgensen et al., 2018).

Carbon Isotope Fractionation

$^{12}$7C/$^{12}$8C is modulated by:

• Isotope-selective photodissociation of CO: In photon-dominated regions, $^{12}$9CO self-shields more effectively, causing $^{13}$0CO/$^{13}$1CO $^{13}$2 elemental value at cloud surfaces (Colzi et al., 2020).
• Low-temperature isotopic exchange: Exothermic reactions such as $^{13}$3C$^{13}$4 + CO $^{13}$5 $^{13}$6CO + C$^{13}$7 ($^{13}$8 K) drive $^{13}$9C into tightly-bound species, depleting $${}^{12}\mathrm{C}/{}^{13}\mathrm{C} = \frac{N(\text{species with }^{12}\mathrm{C})}{N(\text{species with }^{13}\mathrm{C})}$$0C from radicals and chains at 10–20 K. Fractionation is reversed via atomic C + C$${}^{12}\mathrm{C}/{}^{13}\mathrm{C} = \frac{N(\text{species with }^{12}\mathrm{C})}{N(\text{species with }^{13}\mathrm{C})}$$1 exchange ($${}^{12}\mathrm{C}/{}^{13}\mathrm{C} = \frac{N(\text{species with }^{12}\mathrm{C})}{N(\text{species with }^{13}\mathrm{C})}$$2 K) late in cloud evolution, transiently enhancing $${}^{12}\mathrm{C}/{}^{13}\mathrm{C} = \frac{N(\text{species with }^{12}\mathrm{C})}{N(\text{species with }^{13}\mathrm{C})}$$3C in C-chains and nitriles (Colzi et al., 2020).

Gas-grain models show that $${}^{12}\mathrm{C}/{}^{13}\mathrm{C} = \frac{N(\text{species with }^{12}\mathrm{C})}{N(\text{species with }^{13}\mathrm{C})}$$4C/$${}^{12}\mathrm{C}/{}^{13}\mathrm{C} = \frac{N(\text{species with }^{12}\mathrm{C})}{N(\text{species with }^{13}\mathrm{C})}$$5C in HCN, HNC, or CN may range from $${}^{12}\mathrm{C}/{}^{13}\mathrm{C} = \frac{N(\text{species with }^{12}\mathrm{C})}{N(\text{species with }^{13}\mathrm{C})}$$650 (strong $${}^{12}\mathrm{C}/{}^{13}\mathrm{C} = \frac{N(\text{species with }^{12}\mathrm{C})}{N(\text{species with }^{13}\mathrm{C})}$$7C enrichment) to $${}^{12}\mathrm{C}/{}^{13}\mathrm{C} = \frac{N(\text{species with }^{12}\mathrm{C})}{N(\text{species with }^{13}\mathrm{C})}$$8200 (dilution), depending on temperature, density, CO freeze-out stage, and cosmic-ray ionization rate. This variation strongly impacts isotopic interpretations in star-forming regions.

4. Spatial Variation and Environmental Diagnostics

Galactic Gradient and ISM Baseline

The $${}^{12}\mathrm{C}/{}^{13}\mathrm{C} = \frac{N(\text{species with }^{12}\mathrm{C})}{N(\text{species with }^{13}\mathrm{C})}$$9C/$_3$0C ratio increases with Galactocentric radius, following

$_3$1

with $_3$2 in the Solar neighborhood and $_3$340 in the Galactic Center (diffuse gas) (Luo et al., 2024). Dense-gas tracers in the inner Galaxy yield ratios as low as 11–24 due to opacity and environmental fractionation.

Protoplanetary Disks

Radially resolved ALMA data demonstrate:

• D/H (DCN/HCN) is approximately $_3$4 radially in the PDS 70 disk, i.e., two orders of magnitude above ISM elemental, and flat from 40–100 au.
• $_3$5C/$_3$6C (HCN/H$_3$7CN) profile is constant at $_3$869, matching ISM levels, with no significant radial trend (Rampinelli et al., 4 Apr 2025).
• In contrast, HCN/HC$_3$9N exhibits strong radial gradients, diagnostic of N-fractionation via isotope-selective photodissociation of N$^{13}$0, while carbon and hydrogen fractionation remain comparatively invariant in the outer disk.

Protosolar and Planetary Building Block Contexts

Cometary measurements show:

• D/H in water spans a wide range: as low as $^{13}$1 (VSMOW-like; C/2014 Q2 Lovejoy) to as high as $^{13}$2 (C/2012 F6 Lemmon). In comet 67P, D/H in water is $^{13}$3, with higher D/H in organics (alkanes: $^{13}$4–$^{13}$5) (Müller et al., 2022, Biver et al., 2016).
• $^{13}$6C/$^{13}$7C in HCN or alkanes is solar/terrestrial within uncertainties ($^{13}$889–124), but $^{13}$9C-enrichment is noted in specific protostellar organics (Jørgensen et al., 2018).
• In protostars, the $^{12}$0C/$^{12}$1C ratio in O-bearing organics (dimethyl ether, glycolaldehyde, methyl formate) can be as low as 25–40, half the local ISM value, possibly due to selective retention of $^{12}$2CO in ices or UV-driven fractionation (Jørgensen et al., 2018).

5. Implications for Volatile Evolution, Chemical Pathways, and Cosmochemistry

• Fractionation as a probe of origin: Variation of D/H and $^{12}$3C/$^{12}$4C enables reconstruction of volatile delivery to planet-forming disks and planetary atmospheres, tracing inheritance vs. disk chemistry (Rampinelli et al., 4 Apr 2025).
• Constraints on molecular formation: Nearly equal H$^{12}$5CCCN and HC$^{12}$6CCN abundances in HC$^{12}$7N confirm a two-equivalent-carbon atom progenitor (C$^{12}$8H$^{12}$9), while $^{13}$00C-enrichment at the central C position arises from CN supplied by a low $^{13}$01C/$^{13}$02C reservoir. Post-formation isotopic scrambling is ruled out (Araki et al., 2016).
• Limitations and caveats: Optically thick tracers, local excitation, and unknown spatial source structure can bias derived ratios, particularly in emission-line studies of dense regions (Luo et al., 2024, Colzi et al., 2020).
• Disk-averaged vs. spatially resolved fractionation: Direct measurement of radial profiles (disk, core, envelope) resolves degenerate interpretations associated with averaged values, revealing active spatially varying fractionation (Rampinelli et al., 4 Apr 2025).
• Solar System connection: D/H and $^{13}$03C/$^{13}$04C in cometary ices overlap or bracket terrestrial and protosolar values, with organics generally more enriched in deuterium than water, reflecting inheritance from, and chemical processing in, low-temperature presolar and protosolar environments (Müller et al., 2022, Biver et al., 2016).

6. Theoretical and Modeling Perspectives

• Detailed gas-grain chemical networks, incorporating both gas-phase fractionation channels (ion-molecule, neutral-neutral exchange, PDR photodissociation) and ice-surface processes (accretion, photodesorption, binding energy effects), reproduce the observed diversity in observed ratios. Critical dependence on temperature, density, and evolutionary time is demonstrated (Colzi et al., 2020).
• The newly proposed C$^{13}$05 exchange reaction ($^{13}$06C + C$^{13}$07 $^{13}$08 $^{13}$09CC$^{13}$10 + C) is efficient post-CO-freezeout and drives late-stage $^{13}$11C enhancement in C-chains (Colzi et al., 2020).
• The "double-isotope method"—using observed ratios of H$^{13}$12CN/H$^{13}$13CN and an assumed $^{13}$14C/$^{13}$15C to infer $^{13}$16N/$^{13}$17N—can yield errors up to a factor of 3.5 if actual $^{13}$18C/$^{13}$19C is not explicitly measured and is variable (Colzi et al., 2020).

7. Comparative Synthesis and Outlook

D/H and $^{13}$20C/$^{13}$21C ratios are foundational isotopic metrics for deciphering the chemical evolution of astrophysical matter. Their measured values and radial profiles in environments from cold molecular clouds to protostellar systems, protoplanetary disks, and comets encode the interplay of fractionation processes, inheritance, spatial mixing, and disk chemistry. The emerging observational paradigm, supported by refined gas-grain models and spatially resolved measurements, emphasizes environment-specific, temporally evolved, and process-dependent isotopic signatures that inform solar and exoplanetary volatile origins. Ongoing improvements in sensitivity, coverage, and spatial resolution—particularly with ALMA and in situ solar system probes—will further quantify isotopic variability, elucidating the chemical and dynamical history of solid and volatile material from interstellar to planetary scales (Luo et al., 2024, Rampinelli et al., 4 Apr 2025, Jørgensen et al., 2018, Araki et al., 2016, Müller et al., 2022, Biver et al., 2016, Colzi et al., 2020).