Interstellar Chemistry and Molecule Formation

The interstellar medium is a chemistry set run at the wrong temperature. Number densities span $\sim1\,\mathrm{cm^{-3}}$ in the diffuse gas to $\gtrsim10^6\,\mathrm{cm^{-3}}$ in prestellar cores, and temperatures sit between $10$ and $100\,\mathrm{K}$ — cold enough that any reaction with an activation barrier is exponentially dead, and dilute enough that the three-body encounters which stabilize products in a flask or a stellar atmosphere essentially never happen. Yet space is full of molecules: over three hundred distinct species, from $\mathrm{H_2}$ and CO to methanol, glycolaldehyde, and long carbon chains. The whole subject is the story of how a rich chemistry survives under conditions that forbid almost every reaction we learn on Earth. Two loopholes carry the entire enterprise. The first is the ion-molecule reaction: a charge polarizes a neutral from far away, creating a long-range $-1/r^4$ attraction that pulls the pair together with no barrier to climb — so these reactions run just as fast at $10\,\mathrm{K}$ as at room temperature. The second is the dust grain, a solid surface that acts as the third body a two-body gas collision lacks, soaking up the binding energy released when two atoms join. Grains are why molecular hydrogen — the coolant, shield, and feedstock for everything else — exists at all, since two H atoms meeting in the gas cannot radiate away their bond energy and simply fly apart. Layer cosmic rays on top and the network ignites. Cosmic-ray protons plunge into clouds that starlight cannot reach, ionize $\mathrm{H_2}$ , and seed $\mathrm{H_3^+}$ — the proton donor that launches oxygen, carbon, and nitrogen chemistry. The abundances we observe are never an equilibrium; they are a time-dependent record of a cloud's history, which is what makes interstellar chemistry a clock and a thermometer for star formation.

Key takeaways

At $10$ – $100\,\mathrm{K}$ only barrierless reactions proceed, so interstellar chemistry runs on two channels: gas-phase ion-molecule reactions and catalysis on dust-grain surfaces.
Ion-molecule reactions are fast and temperature-independent because the charge-induced-dipole (Langevin) attraction gives a capture rate $k_L=2\pi q\sqrt{\alpha/\mu}\sim10^{-9}\,\mathrm{cm^3\,s^{-1}}$ with no explicit $T$ .
Molecular hydrogen cannot form in the gas — radiative association of two H atoms is dipole-forbidden ( $\sim10^{-29}\,\mathrm{cm^3\,s^{-1}}$ ) — so essentially all $\mathrm{H_2}$ forms on grains at $R\approx3\times10^{-17}\,\mathrm{cm^3\,s^{-1}}$ .
Cosmic rays ionize $\mathrm{H_2}$ to $\mathrm{H_2^+}$ , which makes $\mathrm{H_3^+}$ , the universal proton donor that drives chains to $\mathrm{H_2O}$ , $\mathrm{OH}$ , $\mathrm{NH_3}$ , hydrocarbons, and $\mathrm{HCO^+}$ .
Abundances are set by competition between formation, destruction by photons and cosmic rays, freeze-out onto grains, and desorption — a non-equilibrium network whose composition encodes cloud history.
Key idea: at ISM temperatures only barrierless reactions proceed — Langevin-fast ion-molecule reactions build the gas-phase network from cosmic-ray-made $\mathrm{H_3^+}$ , while dust grains catalyse $\mathrm{H_2}$ formation and the hydrogenation reactions gas-phase chemistry cannot achieve.

The core insight

At ISM temperatures only barrierless reactions proceed: ion-molecule reactions (Langevin-fast, temperature-independent) build the gas-phase network from cosmic-ray-produced $\mathrm{H_3^+}$ , while dust grains catalyse $\mathrm{H_2}$ formation and hydrogenation reactions that gas-phase chemistry cannot achieve.

Where this lives today

Astrochemistry is now a precision observational science, and its center of gravity is ALMA — the Atacama Large Millimeter/submillimeter Array — which resolves molecular line emission in protoplanetary disks and hot cores and has turned "over 300 detected species" into a working periodic table of the ISM, including prebiotic molecules such as glycolaldehyde (a simple sugar) and the first aromatic interstellar organics. The James Webb Space Telescope now measures ice-mantle abundances — $\mathrm{H_2O}$ , CO, $\mathrm{CO_2}$ , $\mathrm{CH_3OH}$ — directly in absorption toward embedded protostars, tying the grain-surface chemistry of this article to the material delivered into planet-forming disks. Molecular ions serve as cosmic-ray dosimeters: infrared absorption by $\mathrm{H_3^+}$ and submillimeter lines of $\mathrm{OH^+}$ and $\mathrm{H_2O^+}$ pin down $\zeta$ along individual sightlines and reveal its order-of-magnitude variation across the Galaxy. And deuterium fractionation — the enhancement of species like $\mathrm{N_2D^+}$ in the coldest cores — is now the standard chemical clock for prestellar cores near collapse. The open frontiers: the origin of complex organics (warm gas-phase versus grain-surface radical routes), the true value and variability of $\zeta$ , and how much of a comet's volatile inventory is inherited chemically unchanged from the parent cloud.

The physical picture

Picture chemistry as a board game played in a deep freeze, where the only legal moves release energy and need no push to start. Two neutral atoms drifting together feel almost nothing until they touch, and joining them across even a modest barrier is forbidden — the Boltzmann factor $e^{-E_a/kT}$ with $kT\approx10^{-3}\,\mathrm{eV}$ at $10\,\mathrm{K}$ crushes any barrier of a fraction of an electron-volt to nothing. Ions rewrite the rules: a charge polarizes a passing neutral, and the attraction curves distant trajectories inward, guaranteeing a reaction. Dust grains rewrite them again — a grain is a solid third body that absorbs the recoil energy a lone pair of colliding atoms has nowhere to dump. The first figure shows these two channels side by side, the only two ways to make molecules in the cold.

Figure

Two barrierless channels of cold chemistry. On a grain, two adsorbed H atoms meet and the new H₂ carries off part of the 4.48 eV bond energy while the lattice absorbs the rest, at rate R ≈ 3×10⁻¹⁷ cm³ s⁻¹. In the gas, a cosmic ray ionizes H₂ and launches the H₃⁺ ion-molecule network at Langevin rate k_L ≈ 10⁻⁹ cm³ s⁻¹, independent of temperature.

Move from a single reaction to a whole cloud and a second organizing idea appears: depth. Ultraviolet starlight photodissociates $\mathrm{H_2}$ and CO but is absorbed in the outer skin of a cloud, so molecules survive only where the dust column — measured by the visual extinction $A_V$ — is thick enough to shield them. Because different molecules need different shielding, a cloud is layered like an onion: atomic hydrogen gives way to $\mathrm{H_2}$ at $A_V\approx1$ , while ionized carbon $\mathrm{C^+}$ becomes neutral C and then locks into CO only deeper, near $A_V\approx2$ – $4$ . This stratification is the most important structural fact about molecular gas, and it is why a layer of "CO-dark" molecular gas exists — rich in $\mathrm{H_2}$ but with no CO to trace it.

Figure

Chemical stratification of a cloud with depth. As the shielding column A_V increases into the cloud, hydrogen turns molecular at A_V ≈ 1 (blue, fraction 2n(H₂)/n_H) and carbon locks into CO only deeper, near A_V ≈ 2 to 4 (green). Between the two thresholds lies CO-dark molecular gas: H₂-rich but CO-poor.

The last picture to hold is that abundances are a story, not a state. Reaction timescales in the cold ISM run to $10^5$ – $10^7$ years, comparable to how long a cloud lives, so the chemistry never relaxes to equilibrium — what we read in a spectrum is a snapshot of a running clock, and chemistry becomes a diagnostic of physics.

The core equations

Interstellar chemistry is a set of coupled rate equations, but a handful of relations carry the physics: the master balance for $\mathrm{H_2}$ , the two microphysical rates that appear everywhere (grain catalysis and Langevin capture), the cosmic-ray trigger, and the freeze-out clock. The abundance of any species is a race between formation and destruction. For $\mathrm{H_2}$ the formation term is grain catalysis and the dominant destruction is photodissociation, and their balance controls the atomic-to-molecular transition.

\frac{\mathrm{d}n(\mathrm{H_2})}{\mathrm{d}t}=R\,n_{\mathrm{H}}\,n(\mathrm{H})-\zeta_{\mathrm{diss}}\,n(\mathrm{H_2}) \qquad (1)

where: $R$ — grain-surface $\mathrm{H_2}$ formation rate coefficient, $\approx3\times10^{-17}\,\mathrm{cm^3\,s^{-1}}$ ; $n_{\mathrm{H}}$ — total hydrogen-nucleon density, $n_{\mathrm{H}}=n(\mathrm{H})+2n(\mathrm{H_2})$ ; $n(\mathrm{H})$ — atomic-hydrogen density; $\zeta_{\mathrm{diss}}$ — photodissociation rate of $\mathrm{H_2}$ per molecule. Reading: molecules are made on grains at a rate set by how often an H atom meets a grain, and destroyed by UV photons; where shielding kills $\zeta_{\mathrm{diss}}$ , the ratio $n(\mathrm{H_2})/n(\mathrm{H})$ climbs steeply — the $\mathrm{H\,I}$ -to- $\mathrm{H_2}$ transition. The formation coefficient $R$ is not fundamental; it is a product of microphysical factors — how fast H atoms move, how much grain surface they can hit, and how efficiently a stuck atom actually leaves as $\mathrm{H_2}$ .

R=\tfrac12\,\langle v_{\mathrm{H}}\rangle\left(\frac{n_{\mathrm{gr}}\sigma_{\mathrm{gr}}}{n_{\mathrm{H}}}\right)S(T)\,\epsilon \qquad (2)

where: $\langle v_{\mathrm{H}}\rangle=\sqrt{8kT/\pi m_{\mathrm{H}}}$ — mean thermal speed of H atoms; $n_{\mathrm{gr}}\sigma_{\mathrm{gr}}/n_{\mathrm{H}}$ — dust geometric cross section per H nucleon, $\approx10^{-21}\,\mathrm{cm^2}$ ; $S(T)$ — sticking probability of an impinging H atom; $\epsilon$ — probability an adsorbed atom forms $\mathrm{H_2}$ before it desorbs. Reading: rate $=\tfrac12\times$ (H thermal flux) $\times$ (grain area per H) $\times S\times\epsilon$ , the $\tfrac12$ because two atoms make one molecule. With $\langle v_{\mathrm{H}}\rangle\sim1.5\times10^5\,\mathrm{cm\,s^{-1}}$ and $S\epsilon\sim0.3$ this recovers $R\approx3\times10^{-17}\,\mathrm{cm^3\,s^{-1}}$ . The engine of the gas phase is the ion-molecule reaction, whose rate follows from the classical capture of a neutral by the long-range field of an ion. The remarkable feature is what is absent from the answer.

k_L=2\pi q\sqrt{\frac{\alpha}{\mu}} \qquad (3)

where: $q$ — ion charge; $\alpha$ — polarizability volume of the neutral; $\mu$ — reduced mass of the colliding pair. In Gaussian units with $q=4.8\times10^{-10}\,\mathrm{esu}$ and $\alpha\sim10^{-24}\,\mathrm{cm^3}$ , $k_L\sim10^{-9}\,\mathrm{cm^3\,s^{-1}}$ . Reading: there is no temperature and no barrier here. Because the capture cross section scales as $1/v$ , the product $\langle\sigma v\rangle$ is velocity-independent — ion-molecule chemistry runs at the same brisk rate at $10\,\mathrm{K}$ as in a laboratory. Every gas-phase chain has the same ignition: a cosmic ray ionizes $\mathrm{H_2}$ , and the ion instantly reacts with another $\mathrm{H_2}$ to make the proton donor $\mathrm{H_3^+}$ .

\mathrm{H_2}+\mathrm{c.r.}\to\mathrm{H_2^+}+e^-,\qquad \mathrm{H_2^+}+\mathrm{H_2}\to\mathrm{H_3^+}+\mathrm{H} \qquad (4)

where: $\mathrm{c.r.}$ — cosmic-ray proton or fast secondary electron; $\mathrm{H_3^+}$ — triatomic hydrogen ion, the universal proton donor of interstellar chemistry. Reading: the second reaction is Langevin-fast, so $\mathrm{H_2^+}$ never accumulates; the rate-limiting step is the first, set by the cosmic-ray ionization rate $\zeta\sim10^{-16}\,\mathrm{s^{-1}}$ (a few $\times10^{-17}$ in dense cores, $\sim10^{-15}$ in diffuse gas, measured from $\mathrm{H_3^+}$ and $\mathrm{OH^+}$ , $\mathrm{H_2O^+}$ ). Control $\zeta$ and you control the whole network. Finally, in dense gas molecules do not just react — they stick. The freeze-out clock says how fast a gas-phase species accretes onto cold grains, and it is short enough in cores to matter.

t_{\mathrm{freeze}}\approx\frac{1}{n_{\mathrm{gr}}\sigma_{\mathrm{gr}}\,\langle v\rangle\,S}\approx\frac{10^{9}}{n_{\mathrm{H}}\,[\mathrm{cm^{-3}}]}\ \mathrm{yr} \qquad (5)

where: $t_{\mathrm{freeze}}$ — timescale for a molecule to accrete onto grain surfaces; $\langle v\rangle$ — thermal speed of the accreting species; $n_{\mathrm{H}}$ — hydrogen density. Reading: at $n_{\mathrm{H}}\sim10^5\,\mathrm{cm^{-3}}$ this is only $\sim10^4$ years — far shorter than a core's life — so CO and other molecules deplete out of the gas onto icy mantles, reshaping the chemistry we observe.

Derivation AWhy molecular hydrogen cannot form in the gas phase

The goal is to show quantitatively that the two gas-phase routes to $\mathrm{H_2}$ — radiative association of two hydrogen atoms and three-body recombination — are both hopelessly slow under interstellar conditions, so that grain-surface catalysis is not merely helpful but strictly necessary. We estimate each gas-phase route in turn, then show how a grain rescues the reaction.

Step 1. Two ground-state H atoms colliding must shed the $4.48\,\mathrm{eV}$ bond energy during the encounter to end up bound. But $\mathrm{H_2}$ is homonuclear: it has no permanent electric dipole moment and therefore no allowed electric-dipole transition to radiate that energy away. Stabilization can only proceed through a forbidden quadrupole transition with radiative rate $A\sim10^{-7}\,\mathrm{s^{-1}}$ . Why: the selection rule, not the energetics, is the killer — the energy is available, but the molecule has no fast way to emit it.

Step 2. The transient collision complex lasts only about a vibrational period, $\tau\sim10^{-13}\,\mathrm{s}$ . The probability that it radiates before flying apart is $A\tau\sim10^{-7}\times10^{-13}\sim10^{-20}$ . Folding this tiny probability into a gas-kinetic collision rate gives an effective radiative-association coefficient $k_{\mathrm{ra}}\sim10^{-29}\,\mathrm{cm^3\,s^{-1}}$ . Why: multiplying a normal collision rate by a $10^{-20}$ survival probability is what turns a fast process into a negligible one.

Step 3. Three-body recombination, $\mathrm{H+H+H}\to\mathrm{H_2+H}$ , uses the third atom to carry off the energy and needs no radiation, but its rate scales as $n(\mathrm{H})^3$ ; at ISM densities of $10^2$ – $10^6\,\mathrm{cm^{-3}}$ triple encounters are astronomically rare, becoming competitive only above $n\sim10^8\,\mathrm{cm^{-3}}$ (protostellar-disk interiors). Why: the cube of a small density is vanishingly small, so a route that dominates in a dense flask is irrelevant in a cloud.

Step 4. A dust grain supplies the missing third body. An H atom strikes a grain, sticks, hops across the surface, and finds a second adsorbed H; the new $\mathrm{H_2}$ shares its $4.48\,\mathrm{eV}$ with the whole lattice, which acts as a heat sink and often ejects the molecule in excited rovibrational states (seen in the near-infrared). Standard grain abundances give $R\approx3\times10^{-17}\,\mathrm{cm^3\,s^{-1}}$ — twelve orders of magnitude above the gas-phase coefficient. Why: the solid absorbs the recoil a two-body gas collision cannot, converting a forbidden gas reaction into a routine surface one.

Go deeper · Derivation B · 1 minThe Langevin capture rate and why ion-molecule chemistry ignores temperature

Derivation BThe Langevin capture rate and why ion-molecule chemistry ignores temperature

The goal is to derive from classical mechanics the capture cross section of an ion by a polarizable neutral, and to show that when it is thermally averaged the temperature drops out entirely, leaving the constant Langevin rate coefficient that makes cold interstellar gas chemically active. We build the potential, add the centrifugal barrier, find the critical impact parameter, and average over velocity to expose the cancellation that removes any $T$ dependence.

Step 1. The ion's field $E=q/r^2$ induces a dipole $p=\alpha E$ in the neutral, and the interaction energy is $U(r)=-\tfrac12\alpha E^2=-\alpha q^2/2r^4$ — a long-range attraction far softer than a bare Coulomb potential but reaching well beyond the geometric size of either partner. Why: the $-1/r^4$ tail is what lets an ion "grab" a neutral from a distance, giving a cross section much larger than $\pi d^2$ .

Step 2. For a collision with speed $v$ and impact parameter $b$ , angular momentum $L=\mu v b$ adds a centrifugal term, so the effective radial potential is $U_{\mathrm{eff}}(r)=\tfrac12\mu v^2 b^2/r^2-\alpha q^2/2r^4$ . This has a maximum — a centrifugal barrier — whose height grows with $b$ . Why: whether the pair spirals together or is flung apart is decided by this barrier, and the barrier is tuned by the impact parameter.

Step 3. Setting $\mathrm{d}U_{\mathrm{eff}}/\mathrm{d}r=0$ and evaluating $U_{\mathrm{eff}}$ at the barrier gives a barrier height $E^2 b^4/(4\alpha q^2)$ , with $E=\tfrac12\mu v^2$ . Capture occurs when $E$ exceeds it, i.e. for $b^4<b_{\mathrm{crit}}^4=4\alpha q^2/(\mu v^2)$ , so $\sigma_L=\pi b_{\mathrm{crit}}^2=2\pi q\sqrt{\alpha/\mu}\,/\,v$ . Why: this is the crucial result — $\sigma_L$ scales as $1/v$ , because slower collisions are captured from larger impact parameters.

Step 4. The rate coefficient is the velocity average $k_L=\langle\sigma_L v\rangle$ . But $\sigma_L v=2\pi q\sqrt{\alpha/\mu}$ is already velocity-independent, so the average is trivial: $k_L=2\pi q\sqrt{\alpha/\mu}$ , with no temperature anywhere. Why: the $1/v$ in $\sigma_L$ exactly cancels the $v$ in the flux, which is the deep reason ion-molecule reactions run undiminished at $10\,\mathrm{K}$ and thus dominate cold gas-phase chemistry.

Go deeper · Derivation C · 1 minThe ion-molecule chain from H₃⁺ to interstellar water

Derivation CThe ion-molecule chain from H₃⁺ to interstellar water

The goal is to trace the pathway by which cosmic-ray-produced $\mathrm{H_3^+}$ delivers oxygen into water and hydroxyl, showing how one initiating ionization propagates through a sequence of fast, exothermic ion-molecule steps and one dissociative recombination, and why the identical template with C or N substituted builds the rest of the network. We follow the proton from $\mathrm{H_3^+}$ to oxygen and up the hydrogenation ladder to the electron that closes the cycle.

Step 1. $\mathrm{H_3^+}$ hands a proton to atomic oxygen, $\mathrm{H_3^+}+\mathrm{O}\to\mathrm{OH^+}+\mathrm{H_2}$ . (In practice O is frequently first ionized by charge exchange with $\mathrm{H^+}$ , after which $\mathrm{O^+}+\mathrm{H_2}\to\mathrm{OH^+}+\mathrm{H}$ ; either entry point yields $\mathrm{OH^+}$ .) Why: this is the gateway — it converts an unreactive neutral atom into a molecular ion that can climb the hydrogenation ladder using the most abundant reactant available, $\mathrm{H_2}$ .

Step 2. Successive hydrogen abstractions from $\mathrm{H_2}$ build up the hydronium ion: $\mathrm{OH^+}+\mathrm{H_2}\to\mathrm{H_2O^+}+\mathrm{H}$ , then $\mathrm{H_2O^+}+\mathrm{H_2}\to\mathrm{H_3O^+}+\mathrm{H}$ . Each step is a fast, exothermic Langevin reaction with $\mathrm{H_2}$ , so the chain runs to completion wherever $\mathrm{H_2}$ is abundant. Why: because $\mathrm{H_2}$ is everywhere in molecular gas, these steps are effectively unlimited, and the bottleneck stays at Step 1, tying the water yield to the oxygen budget and to $\zeta$ .

Step 3. The chain terminates when $\mathrm{H_3O^+}$ recombines with a free electron: $\mathrm{H_3O^+}+e^-\to\mathrm{OH}+2\mathrm{H}$ or $\to\mathrm{H_2O}+\mathrm{H}$ , with roughly a quarter of recombinations yielding water and most of the rest hydroxyl. Why: dissociative recombination is how a charged network hands its products back as neutrals — it is the exit door that populates the observable molecules OH and $\mathrm{H_2O}$ .

Step 4. The same template, with C or N in place of O, builds hydrocarbons and ammonia, while proton transfer to CO makes the tracer ion $\mathrm{HCO^+}$ — so a single cosmic-ray ionization cascades through the whole oxygen, carbon, and nitrogen chemistry. Why: this universality is what lets one parameter, the cosmic-ray ionization rate, govern the chemical activity of a shielded cloud.

Limiting cases

Which chemistry operates depends sharply on density, temperature, and shielding. The same gas passes through several of these regimes as it collapses from diffuse cloud to hot core.

Case	Result	Meaning
Diffuse, unshielded ( $n\sim10\,\mathrm{cm^{-3}}$ )	$A_V\lesssim1$ , gas mostly $\mathrm{H\,I}$	UV photodissociates $\mathrm{H_2}$ and CO; only trace molecules survive in self-shielded pockets; chemistry is sparse.
Dense, shielded core ( $n\gtrsim10^4\,\mathrm{cm^{-3}}$ )	$A_V\gtrsim5$ , $T\approx10\,\mathrm{K}$	UV excluded; cosmic-ray-driven ion-molecule chemistry rules; rich inventory; heavy CO freeze-out onto grains.
Protostellar hot core	$T\gtrsim100\,\mathrm{K}$	Ice mantles sublimate, releasing $\mathrm{H_2O}$ , $\mathrm{CH_3OH}$ , and complex organics; warm reactions with small barriers switch on.
Shock / high-temperature	$T\gtrsim300$ – $1000\,\mathrm{K}$	Endothermic neutral-neutral reactions (e.g. $\mathrm{O}+\mathrm{H_2}\to\mathrm{OH}+\mathrm{H}$ ) unlock, rapidly driving oxygen into water behind shocks.

Common trapDo not assume interstellar chemistry is in equilibrium

It almost never is. Formation and destruction timescales in cold gas run to $10^5$ – $10^7$ years — comparable to cloud lifetimes — so abundances depend on history, not a steady state: the age since the gas became molecular, its irradiation and shock record, its freeze-out cycles. Two clouds at identical density and temperature can show very different spectra because they sit at different points along their chemical evolution. A related trap is to estimate the $\mathrm{H_2}$ formation rate from gas-phase radiative association: that coefficient, $\sim10^{-29}\,\mathrm{cm^3\,s^{-1}}$ , is twelve orders of magnitude too small, so neglecting grain catalysis means no cloud would ever turn molecular. Always put $\mathrm{H_2}$ on the grains.

Worked examples

1 · The Langevin rate for H₃⁺ + CO → HCO⁺ + H₂ (fully worked)

Estimate the ion-molecule rate coefficient for the reaction that makes the classic tracer ion $\mathrm{HCO^+}$ , using the polarizability of CO, $\alpha=1.95\times10^{-24}\,\mathrm{cm^3}$ , and compare with the measured value $\approx1.7\times10^{-9}\,\mathrm{cm^3\,s^{-1}}$ .

Step 1. Find the reduced mass. The partners are $\mathrm{H_3^+}$ (mass $3\,\mathrm{u}$ ) and CO (mass $28\,\mathrm{u}$ ): $\mu=\dfrac{3\times28}{3+28}\,\mathrm{u}=2.71\,\mathrm{u}=2.71\times1.66\times10^{-24}\,\mathrm{g}=\mathbf{4.50\times10^{-24}\,g}$ .

Step 2. Assemble the Langevin formula in Gaussian units with $q=4.80\times10^{-10}\,\mathrm{esu}$ . First the square root: $\sqrt{\alpha/\mu}=\sqrt{1.95\times10^{-24}/4.50\times10^{-24}}=\sqrt{0.433}=0.658\ \mathrm{cm^{3/2}\,g^{-1/2}}$ .

Step 3. Then $k_L=2\pi q\sqrt{\alpha/\mu}=6.283\times(4.80\times10^{-10})\times0.658=\mathbf{2.0\times10^{-9}\,cm^3\,s^{-1}}$ . (The units work: $\mathrm{esu}=\mathrm{g^{1/2}\,cm^{3/2}\,s^{-1}}$ , so $\mathrm{esu}\times\mathrm{cm^{3/2}\,g^{-1/2}}=\mathrm{cm^3\,s^{-1}}$ .)

Step 4. Compare. The Langevin estimate $2.0\times10^{-9}$ sits just above the laboratory value $1.7\times10^{-9}\,\mathrm{cm^3\,s^{-1}}$ , as it should — Langevin theory is an upper bound (it counts every capture as a reaction), and real reactions occasionally bounce back apart. Agreement to $\sim20\%$ with no adjustable parameters, and crucially no temperature: this rate is the same at $10\,\mathrm{K}$ and $300\,\mathrm{K}$ .

2 · How long does it take gas to become molecular?

A given H atom converts to $\mathrm{H_2}$ on the timescale $t_{\mathrm{H_2}}=1/(R\,n_{\mathrm{H}})$ , since $R\,n_{\mathrm{H}}$ is the rate at which it finds partners on grains. Evaluate it in a diffuse cloud ( $n_{\mathrm{H}}=100\,\mathrm{cm^{-3}}$ ) and a dense clump ( $n_{\mathrm{H}}=10^4\,\mathrm{cm^{-3}}$ ), with $R=3\times10^{-17}\,\mathrm{cm^3\,s^{-1}}$ .

Step 1. Diffuse: $t_{\mathrm{H_2}}=1/(3\times10^{-17}\times100)=1/(3\times10^{-15})=3.3\times10^{14}\,\mathrm{s}$ . Converting with $1\,\mathrm{yr}=3.15\times10^{7}\,\mathrm{s}$ gives $\mathbf{1.1\times10^{7}\,yr}\approx10\,\mathrm{Myr}$ .

Step 2. Dense: the rate is $100\times$ faster, so $t_{\mathrm{H_2}}=\mathbf{1.1\times10^{5}\,yr}$ .

Step 3. Interpret. Ten million years is comparable to the lifetime of a diffuse cloud, which is exactly why diffuse gas is only partly molecular and sits perpetually near its $\mathrm{H\,I}$ -to- $\mathrm{H_2}$ transition, while dense gas turns fully molecular in a small fraction of its life. (By the same token the freeze-out clock $t_{\mathrm{freeze}}\approx10^9/n_{\mathrm{H}}\,\mathrm{yr}$ gives only $\sim10^4\,\mathrm{yr}$ at $n_{\mathrm{H}}=10^5\,\mathrm{cm^{-3}}$ , so CO depletes onto grains long before such a core collapses.) The molecular fraction of a cloud is a clock reading — it measures how long the gas has been dense.

Physics of the portrait

Interstellar chemistry had no theory for a long time. Molecules were spotted in space decades before anyone could explain how the cold, near-vacuum ISM assembles them; the explanation, when it came, was ion-molecule chemistry — and its architect turned a spectroscopist's puzzle into a whole discipline.

Origins · 1 minWilliam Klemperer and how this was discovered

William Klemperer

1927–2017 · founder of chemical astrophysics and the ion-molecule theory of interstellar chemistry

William Klemperer was a Harvard physical chemist and master molecular spectroscopist who became, almost against the grain of his field, the founder of interstellar chemistry. In 1970 radio astronomers Buhl and Snyder detected a mysterious $89\,\mathrm{GHz}$ emission line they dubbed "X-ogen." Klemperer proposed it came from a molecular ion — $\mathrm{HCO^+}$ , protonated carbon monoxide — a radical suggestion when few believed ions could be abundant in space and no laboratory spectrum of $\mathrm{HCO^+}$ yet existed to check against. He was right, and the identification was confirmed in the laboratory a few years later. With his student Eric Herbst, Klemperer laid out in 1973 the first comprehensive gas-phase ion-molecule network for dense clouds, showing how cosmic-ray ionization and the proton donor $\mathrm{H_3^+}$ build the observed molecular inventory through barrierless reactions that run even at $10\,\mathrm{K}$ . That framework — now called chemical astrophysics — is the backbone of every astrochemical model today, from prestellar cores to protoplanetary disks. His insight was to trust the physics of the charge-induced-dipole attraction over the prejudice that cold, thin gas must be chemically inert.

History

1940

— McKellar identifies interstellar CH and CN in optical absorption, the first known molecules in space

1963

— Gould and Salpeter argue that H₂ forms catalytically on the surfaces of dust grains

1970

— Klemperer identifies the unidentified X-ogen radio line as the molecular ion HCO⁺

1973

— Herbst and Klemperer construct the first comprehensive ion-molecule reaction network

1996

— Geballe and Oka detect H₃⁺ by infrared absorption in dense interstellar clouds

Going further

Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 12 (The Interstellar Medium and Star Formation) · Tielens, The Physics and Chemistry of the Interstellar Medium (2005) · Draine, Physics of the Interstellar and Intergalactic Medium (2011), Ch. 31–33 · van Dishoeck & Blake, Ann. Rev. Astron. Astrophys. 36, 317 (1998) · Herbst & van Dishoeck, Ann. Rev. Astron. Astrophys. 47, 427 (2009).

Provenance. The values $R\approx3\times10^{-17}\,\mathrm{cm^3\,s^{-1}}$ , $E_b(\mathrm{H_2})=4.48\,\mathrm{eV}$ , and $D_0(\mathrm{CO})=11.09\,\mathrm{eV}$ are standard textbook figures (Carroll & Ostlie §12.1; Draine Ch. 31–33; Tielens 2005). A verifier should confirm: the Langevin coefficient $k_L=2\pi q\sqrt{\alpha/\mu}$ from the capture condition $b_{\mathrm{crit}}^4=4\alpha q^2/\mu v^2$ and its $1/v$ cancellation; the worked $k_L\approx2\times10^{-9}$ for $\mathrm{H_3^+}+\mathrm{CO}$ against the measured $1.7\times10^{-9}$ ; the grain rate from Eq. (2) with grain area per H $\approx10^{-21}\,\mathrm{cm^2}$ , $\langle v_{\mathrm{H}}\rangle\sim1.5\times10^5\,\mathrm{cm\,s^{-1}}$ , $S\epsilon\sim0.3$ ; and $t_{\mathrm{H_2}}=1/(Rn_{\mathrm{H}})$ . Reaction pathways follow van Dishoeck & Blake (1998) and Herbst & van Dishoeck (2009); the radiative-association estimate $\sim10^{-29}\,\mathrm{cm^3\,s^{-1}}$ follows the forbidden-transition argument in Draine Ch. 31.

Prerequisites. Phases of the Interstellar Medium · Interstellar Dust Grains · Interstellar Gas Heating and Cooling · Molecular Clouds and Giant Molecular Clouds. Leads to. CO as a Molecular Gas Tracer and the X-Factor · The 21 cm Line of Neutral Hydrogen · Thermal Dust Emission and Infrared Radiation.

Ask the tutor: Want to walk through how cosmic-ray ionization, $\mathrm{H_3^+}$ , and grain catalysis combine to set the abundance of any given molecule in a cloud — or how deuterium fractionation turns that chemistry into a clock for star formation?

If H₂ is so hard to make, why is it the most abundant molecule in the universe?

Because although each formation event is slow, hydrogen is overwhelmingly abundant and grains provide a huge, ever-present surface. The per-atom formation time $1/(R n_{\mathrm{H}})$ is about ten million years in diffuse gas but only $\sim10^5$ years in dense clumps, and once formed, $\mathrm{H_2}$ in a shielded cloud is protected by dust extinction and by self-shielding (the $\mathrm{H_2}$ already present absorbs the UV that would dissociate more). So any gas that stays dense and shielded for a few hundred thousand years turns essentially fully molecular — the difficulty is overcome by sheer abundance and time, which is what makes the molecular fraction such a useful clock.

What makes H₃⁺ so central — why not just start chemistry with H₂ itself?

$\mathrm{H_2}$ is stable and unreactive with most neutrals at low temperature, so on its own it does not start chemistry. $\mathrm{H_3^+}$ is different: a weakly bound proton carrier that readily donates $\mathrm{H^+}$ to almost anything it meets, and every such transfer is a fast, barrierless, exothermic ion-molecule reaction. Cosmic rays make $\mathrm{H_3^+}$ continuously (ionize $\mathrm{H_2}$ , then $\mathrm{H_2^+}+\mathrm{H_2}\to\mathrm{H_3^+}+\mathrm{H}$ ), so it is always available deep inside clouds where UV cannot reach. By protonating O, C, N, and CO it opens the hydrogenation ladders that build water, hydrocarbons, ammonia, and $\mathrm{HCO^+}$ — in effect converting the raw energy of cosmic rays into chemical initiation.

Does the Langevin rate really have no temperature dependence at all?

For the idealized capture it does not: the $1/v$ scaling of the cross section cancels the $v$ in the flux, so $k_L=2\pi q\sqrt{\alpha/\mu}$ contains no $T$ . This is an excellent approximation for ions reacting with nonpolar neutrals, which is why so many interstellar rate coefficients are quoted as constants near $\sim10^{-9}\,\mathrm{cm^3\,s^{-1}}$ . Two refinements matter. If the neutral is polar, the extra attraction makes the rate rise toward low temperature (as $T^{-1/2}$ ), so cold gas reacts even faster. And Langevin theory is an upper bound — it counts every capture as a reaction, whereas some complexes redissociate — so measured rates run a factor of a few below $k_L$ . The headline conclusion, fast and essentially $T$ -independent, survives both.

Why is CO used to trace molecular gas if H₂ is the dominant molecule?

Because $\mathrm{H_2}$ is nearly invisible. Being homonuclear it has no permanent dipole moment and hence no dipole rotational lines; its lowest transitions are weak quadrupole lines needing warm gas ( $\sim100\,\mathrm{K}$ ) to excite, so cold clouds emit essentially no $\mathrm{H_2}$ radiation. CO, by contrast, has a small dipole moment and a bright rotational line at $2.6\,\mathrm{mm}$ easily excited at $10\,\mathrm{K}$ and observable across the Galaxy and beyond. So we observe CO and convert to $\mathrm{H_2}$ mass using an empirical ratio (the X-factor). The catch is the CO-dark layer between the H/ $\mathrm{H_2}$ and $\mathrm{C^+}$ /CO transitions — molecular gas with little CO that surveys miss, recovered instead by ionized carbon and dust emission.

How do complex organic molecules form if the gas is only 10 K?

Mostly not in the gas — on the grains, then during warm-up. In a cold core, atoms and radicals (H, O, C, CO) freeze onto grains and are hydrogenated: sequential H additions turn CO into $\mathrm{H_2CO}$ and then methanol, $\mathrm{CH_3OH}$ . As a nearby protostar warms the ice toward $\sim30$ – $100\,\mathrm{K}$ , radicals gain mobility to diffuse and recombine, building larger molecules such as glycolaldehyde; further warming sublimates the whole mantle, producing the rich "hot core" spectra. Some warm gas-phase routes contribute too, and disentangling the two is an active question — but the cold, patient chemistry on ice mantles is where most of the complexity is first assembled.

What exactly does the visual extinction A_V measure, and why does chemistry care?

$A_V$ is the dimming of background starlight in magnitudes caused by dust along the line of sight, and because dust and gas are well mixed it is a direct proxy for the total column — roughly $N_{\mathrm{H}}\approx2\times10^{21}\,A_V\ \mathrm{cm^{-2}}$ . Chemistry cares because that same dust column absorbs the ultraviolet photons that photodissociate molecules. At $A_V\lesssim1$ the UV penetrates and keeps hydrogen atomic and carbon ionized; by $A_V\approx1$ enough UV is blocked (helped by $\mathrm{H_2}$ self-shielding) that hydrogen goes molecular; and only by $A_V\approx2$ – $4$ is CO shielded enough to survive. So $A_V$ is effectively the depth coordinate of the layered chemistry, turning a dust column into a map of which molecules exist where.