🔬 mesoscale physics

Molecular hydrogen formation on dust: The impact of gas-dust drift on formation efficiency

This study utilizes kinetic Monte Carlo simulations to demonstrate that while gas-dust drift increases H-atom collision rates, the resulting higher impact energies can significantly suppress molecular hydrogen formation efficiency by reducing sticking probabilities, thereby challenging the assumption that relative motion always enhances H $_2$ production in dynamically active interstellar environments.

Original authors: Stefan Reissl, Simon C. O. Glover, Ralf S. Klessen, Liam S. Morrissey, Mordecai-Mark Mac Low

Published 2026-10-07

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Stefan Reissl, Simon C. O. Glover, Ralf S. Klessen, Liam S. Morrissey, Mordecai-Mark Mac Low

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the vast, cold spaces between the stars, a quiet but vital transformation is constantly taking place. Hydrogen, the most abundant element in the universe, exists mostly as single atoms drifting through the void. Yet, to build the stars and planets that light up the cosmos, these atoms must pair up to form molecules. This pairing is difficult to achieve in open space because the atoms are too far apart to bump into each other often enough, and when they do, they usually bounce apart rather than sticking together. Nature has found a clever workaround: tiny specks of cosmic dust. These grains, made of materials like silicate rock or carbon, act as microscopic workbenches. Single hydrogen atoms land on them, slide around the surface, and eventually meet a partner to form a molecule of hydrogen gas, which then flies off into space. This process is the starting point for almost all chemistry in the galaxy and helps cool down the clouds of gas where new stars are born.

For decades, astronomers have understood that the temperature of these dust grains and the properties of their surfaces are critical to how well this molecular factory works. However, a new study by researchers in Germany, Canada, and the United States suggests that the story is incomplete without considering how the dust moves relative to the gas. In many active regions of space, such as outflows from young stars or turbulent clouds, the dust grains do not sit still; they drift through the gas at high speeds. The researchers wanted to know how this relative motion changes the efficiency of hydrogen molecule creation. They built a sophisticated computer simulation that tracks individual hydrogen atoms as they collide with dust grains, accounting for the fact that the grains are moving through the gas, which changes how hard the atoms hit and whether they stick.

The team developed a detailed model of spherical dust grains, mimicking both silicate and carbon-based surfaces. They simulated the arrival of hydrogen atoms, which can either stick to the grain and wander across it, or react immediately with an atom already waiting there. A key part of their work was testing two different ideas about how atoms stick to the surface. In one scenario, they assumed the chance of an atom sticking is constant, regardless of how fast it hits. In the other, they used a more realistic rule where the sticking probability drops if the atom hits the surface too hard. This distinction is crucial because the drift of the grain through the gas increases the speed at which atoms strike the surface.

When the researchers ran their simulations with a constant chance of sticking, they found that faster drift actually helped. The moving grain swept up more hydrogen atoms per second, leading to a higher rate of molecule formation. However, the results changed dramatically when they included the rule that harder impacts make atoms less likely to stick. In this more realistic scenario, strong drift became a hindrance. As the grains moved faster through the gas, the hydrogen atoms hit the surface with such force that they bounced off instead of staying. The formation efficiency, which measures how many collisions successfully create a molecule, plummeted. For grains drifting at high speeds, the efficiency dropped to a mere one to three percent, compared to twenty to fifty percent when the drift was slow.

The study also revealed that the location where molecules form on the grain is not uniform. When a grain drifts, it creates a wind of hydrogen atoms hitting one side more than the other. At low speeds, the distribution of reactions is fairly even. But at high speeds, the pattern becomes highly uneven. If the sticking probability depends on impact speed, the most active region for forming molecules is not the point where the gas hits the grain most directly. Instead, it shifts to a ring around the side of the grain. At the very front, the atoms hit so hard they bounce off. At the very back, the wind is too weak to bring many atoms. The sweet spot is in between, where the collision rate is still high, but the impact speed is low enough for the atoms to stick.

Another major finding concerns the type of dust grain. Carbon-based grains proved to be much more resilient than silicate ones. They could continue to form hydrogen molecules efficiently at much higher temperatures, up to 150 Kelvin, whereas silicate grains stopped working effectively around 110 Kelvin. This is because carbon surfaces hold onto hydrogen atoms more tightly, preventing them from evaporating back into space as the grain heats up. Despite this advantage, the negative effect of high-speed drift applied to both types of grains. In almost all cases, the reaction that involves a gas atom hitting a stationary atom on the surface dominated the process, while the reaction where two wandering atoms meet became less important as the drift increased and the surface population thinned out.

The implications of these findings extend beyond just counting molecules. The formation of hydrogen on dust grains releases energy that can spin the grains, helping them align with magnetic fields in space. This alignment is what allows astronomers to measure magnetic fields by observing the polarization of starlight. If the drift of the grains changes how often molecules form, it could also change how the grains spin and align. This suggests that in dynamic environments like stellar outflows or shock waves, where dust and gas move at different speeds, the standard models used to predict molecular formation might be significantly off. The researchers conclude that to truly understand the chemistry of the universe, we must account for the relative motion between gas and dust, as it fundamentally alters the efficiency of the cosmic machinery that builds the building blocks of stars.

Technical Summary: The Impact of Gas–Dust Drift on Molecular Hydrogen Formation Efficiency

Problem Statement
Molecular hydrogen ( $H_2$ ) is the most abundant molecule in the interstellar medium (ISM) and a critical coolant for star-forming regions. Its formation is predominantly catalyzed on dust grain surfaces via the Langmuir-Hinshelwood (LH) and Eley-Rideal (ER) mechanisms. While the dependence of formation efficiency on grain temperature, surface composition (silicate vs. carbonaceous), and morphology is well-studied, the impact of relative motion between gas and dust (gas–dust drift) has been largely neglected in microscopic simulations. In dynamically active environments—such as protostellar outflows, shocks, and turbulent clouds—drift velocities can become comparable to or exceed the thermal velocity of atomic hydrogen. This relative motion alters the collision rate and the impact energy of incident atoms, potentially modifying the sticking probability and the subsequent formation efficiency. This paper investigates how gas–dust drift modifies $H_2$ formation, specifically addressing the competition between drift-enhanced collision rates and the reduced sticking probability at higher impact energies.

Methodology
The authors developed an event-driven Kinetic Monte Carlo (KMC) model to simulate individual hydrogen atoms on spherical silicate and amorphous carbon grains. Key methodological components include:

Grain Model: The model utilizes a Fibonacci distribution scheme to place adsorption sites (physisorption and chemisorption) uniformly on the grain surface. It distinguishes between weakly bound physisorption sites (p-sites) and strongly bound chemisorption sites (c-sites), incorporating transitions between them.
Drift Implementation: Gas–dust drift is modeled by shifting the Maxwellian velocity distribution of incident H atoms by a drift velocity $U$ . A dimensionless drift parameter $s = U/v_{th}$ (where $v_{th}$ is the thermal velocity) is used. The collision rate is corrected using a factor $\kappa(s)$ derived from the shifted distribution.
Sticking Probability: The study compares two scenarios:
1. A constant sticking probability ( $S=0.5$ ).
2. A drift-dependent sticking probability, where the sticking function $S(T_g, T_d)$ is evaluated using an effective gas temperature $T_{eff} = T_g + m_H U_\perp^2 / k_B$ . Here, $U_\perp$ is the component of drift velocity normal to the surface, accounting for the increased kinetic energy of impact.
Reaction Channels: The simulation tracks adsorption, surface migration (thermal hopping and quantum tunneling), thermal desorption, and $H_2$ formation via both LH (two adsorbed atoms reacting) and ER (incident atom reacting with an adsorbed atom) mechanisms.
Parameters: Simulations were run for grain radii of 10 nm, covering a range of dust temperatures ( $T_d$ ), gas temperatures ( $T_g = 100$ K), and drift parameters ( $s = 0.1, 1.0, 10$ ).

Key Results

Anisotropy of Formation: Gas–dust drift induces a strong anisotropy in the spatial distribution of adsorbed H atoms and $H_2$ formation.
- At low drift ( $s=0.1$ ), the distribution is nearly isotropic.
- At high drift ( $s=10$ ) with constant sticking, formation concentrates on the upstream hemisphere due to the enhanced collision rate.
- At high drift ( $s=10$ ) with drift-dependent sticking, the formation pattern shifts. The highest formation rates occur in a ring-like region perpendicular to the drift vector. This is because the leading edge (maximum impact energy) suffers from a drastically reduced sticking probability, while the upstream hemisphere receives fewer collisions at angles where sticking remains viable.
Formation Efficiency ( $\epsilon_{H2}$ ):
- Constant Sticking: Increasing drift generally enhances $H_2$ formation efficiency (up to $\epsilon_{H2} \approx 0.3\text{--}0.4$ ) due to the increased collision rate.
- Drift-Dependent Sticking: Strong drift ( $s=10$ ) significantly suppresses $H_2$ formation for both silicate and carbon grains, reducing efficiency to $\epsilon_{H2} \approx 0.01\text{--}0.03$ . The increased impact energy reduces the probability of atoms sticking to the surface, outweighing the benefit of higher collision rates.
Material Dependence: Carbonaceous grains maintain efficient formation up to higher dust temperatures ( $T_d \approx 150$ K) compared to silicate grains ( $T_d \approx 110$ K) due to deeper chemisorption binding energies. However, strong drift suppresses efficiency for both materials.
Reaction Mechanisms: The ER mechanism dominates $H_2$ formation across the entire parameter space. As drift increases, the surface population of adsorbed atoms decreases (due to reduced sticking), which further suppresses the LH mechanism. Consequently, ER becomes even more dominant at high drift velocities.
Timescales: The characteristic timescale for the grain surface to reach equilibrium occupation ( $\tau_{occ}$ ) increases with drift when using the drift-dependent sticking function, as the reduced sticking slows the accumulation of surface species.

Significance and Conclusions
The paper concludes that gas–dust drift is a critical factor in $H_2$ formation that cannot be ignored in dynamically active astrophysical environments. The primary finding is that an enhanced collision rate does not necessarily imply a higher formation rate; the increased impact energy associated with drift can reduce sticking probabilities to the point where formation is strongly suppressed.

The authors assert that prescriptions for $H_2$ formation in models of turbulent clouds, shocks, or outflows must account for the relative gas–dust velocity. Specifically, the drift velocity affects both the collision rate and the adsorption probability, and these effects can oppose one another. The study suggests that variations in inferred $H_2$ formation rates in observational data (e.g., in Photodissociation Regions or shocks) may reflect gas–dust drift rather than solely variations in density, radiation fields, or grain temperature.

Finally, the paper notes that while the current model provides a first step, significant uncertainties remain regarding the precise sticking functions for amorphous materials at high impact energies. The authors suggest that future work should rely on first-principles molecular dynamics simulations to constrain these sticking probabilities more rigorously.

More like this