🔬 condensed matter

Elucidating the mechanism of helium evaporation from liquid water

Using molecular dynamics simulations and a Fokker-Planck-based continuum theory, this study elucidates that helium evaporation from liquid water produces a super-Maxwellian kinetic energy distribution driven by anomalously low friction at the deformable liquid-vapor interface.

Original authors: Kritanjan Polley, Kevin R. Wilson, David T. Limmer

Published 2026-09-18

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Kritanjan Polley, Kevin R. Wilson, David T. Limmer

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

The boundary where air meets water is far more than a simple line; it is a bustling frontier where molecules constantly cross back and forth, driving the chemistry of our atmosphere and the climate of our planet. When a gas molecule from the air dissolves into the ocean, or when a dissolved substance escapes back into the sky, it must navigate this transition zone. For most substances, this journey is a predictable shuffle, governed by the random jostling of heat that keeps everything in a state of thermal balance. In this balanced state, the speeds of escaping molecules follow a standard pattern known as the Maxwell-Boltzmann distribution, where the average speed is directly tied to the temperature of the liquid. However, nature sometimes breaks these rules, and one particular gas, helium, has long been known to behave strangely when it tries to leave the water.

Researchers have long observed that when helium evaporates from liquid water, it does not follow the standard rules of thermal balance. Instead of the expected distribution of speeds, the escaping helium atoms carry significantly more kinetic energy than the surrounding water temperature would suggest. This phenomenon, known as a super-Maxwellian distribution, has been seen in pure water, salty solutions, and even in the presence of surfactants, yet the exact mechanism driving this extra speed has remained elusive. Previous theories suggested that a specific force, derived from the average energy landscape of the helium atom, was pushing the gas out of the water like a spring releasing a compressed object. While this force certainly exists, it tells only part of the story, as it assumes the atom moves infinitely slowly and ignores the friction it encounters while moving.

To uncover the full story of this energetic escape, a team of scientists at Lawrence Berkeley National Laboratory and the University of California, Berkeley, turned to detailed computer simulations. They constructed a virtual slab of liquid water containing a single helium atom and watched its journey from the deep liquid to the vapor above. By tracking the atom's position and speed over time, they confirmed the experimental findings: the helium atom indeed leaves the water with a mean kinetic energy of 1.50 kcal/mol at 300 Kelvin, which is 0.3 kcal/mol higher than what standard thermal physics predicts. The researchers then sought to explain why this extra energy persists, moving beyond simple force calculations to include the resistance the atom feels as it moves.

The key to the mystery lies in how the water surface behaves around the escaping helium. The team developed a mathematical model that treats the evaporation process as a particle moving through a fluid environment where both the forces and the resistance change rapidly over tiny distances. They discovered that the resistance, or friction, the helium atom feels drops dramatically just before the water density itself begins to fall. This happens because the liquid-vapor interface is not a rigid wall but a flexible surface that can deform. As the helium atom, which hates being surrounded by water, approaches the surface, the water molecules rearrange themselves to create a small pocket or valley around the atom. This deformation allows the helium to experience a gas-like environment even while it is still technically below the average surface level of the water.

Because the helium atom finds itself in this gas-like pocket, the friction holding it back vanishes almost instantly, long before it has fully escaped the liquid. In a normal scenario, friction would act like a brake, slowing a fast-moving particle down until it matches the average speed of its surroundings. Here, however, the friction disappears so quickly that the helium atom is never slowed down. It is pushed out by the repulsive forces of the water but faces almost no resistance to dampen its speed. The result is that the atom retains the extra energy it gained during its escape, shooting out of the water with a velocity that defies the thermal expectations of the liquid it just left.

This new understanding allows the researchers to predict how this process changes with temperature. By adjusting their model to account for how water's viscosity and the energy required to dissolve helium change as the water gets colder, they calculated that the excess kinetic energy of the evaporating helium should remain relatively constant at lower temperatures, around 250 Kelvin, before rising again as the water warms. While this specific temperature dependence has not yet been measured in a laboratory, the model provides a clear, testable prediction for future experiments. The work demonstrates that the evaporation of trace gases is not just a simple diffusion process but a complex interplay where the flexibility of the water surface itself plays a critical role in determining how fast and how energetically a molecule can escape into the air.

Problem Statement
The evaporation of trace gases from liquid water is a critical process in atmospheric chemistry, yet the kinetic behavior of inert species like helium at the air-water interface defies standard thermal expectations. While the evaporation of water and most solutes follows a Maxwell-Boltzmann distribution consistent with thermal equilibrium, experimental and simulation data have consistently shown that evaporating helium exhibits a "super-Maxwellian" distribution of kinetic energies. Previous theoretical attempts to explain this anomaly relied primarily on the potential of mean force (PMF), which describes the average conservative force on a molecule. However, the authors argue that this description is incomplete because it assumes infinitesimally slow motion and neglects dissipative, frictional forces that are essential for accurately describing particle motion at finite velocities. The central problem addressed is the lack of a mechanistic framework that incorporates both the spatially varying conservative forces and the position-dependent friction to explain the excess kinetic energy observed in helium evaporation.

Methodology
The study employs a multi-scale approach combining explicit molecular dynamics (MD) simulations with a reduced continuum theory.

Molecular Dynamics Simulations: Simulations were conducted using LAMMPS with a slab geometry containing 768 water molecules and one helium atom at 300 K. The SWM4-NDP polarizable water model and a Drude oscillator model for helium were used to capture polarization effects. Helium was initialized in the bulk liquid, and trajectories were propagated until the atom evaporated (defined as reaching 20 Å from the local instantaneous interface).
Parameterization of Reduced Model: To bridge the gap between atomistic detail and macroscopic transport, the authors coarse-grained the dynamics into a Fokker-Planck equation. This required two key parameters:
- Potential of Mean Force (PMF), $F(z)$ : Calculated via umbrella sampling to determine the reversible work required to move helium across the interface.
- Position-Dependent Friction, $\gamma(z)$ : Determined by comparing short-time MD trajectories against numerical solutions of the Fokker-Planck equation. An optimization procedure (simulated annealing) was used to fit the friction profile to the MD data, allowing for variations beyond simple sigmoidal forms.
Continuum Modeling: The resulting Fokker-Planck equation, incorporating the computed $F(z)$ and $\gamma(z)$ , was solved to predict the velocity and kinetic energy distributions of evaporating helium.

Key Results

Super-Maxwellian Distribution: The MD simulations reproduced the experimental finding that evaporating helium possesses a mean kinetic energy of 1.50 kcal/mol at 300 K, an excess of 0.3 kcal/mol over the expected $2k_BT$ for a flux-weighted Maxwell-Boltzmann distribution.
Anomalously Low Interfacial Friction: The computed friction profile $\gamma(z)$ reveals a precipitous drop in friction that occurs before the water density profile begins to decrease significantly. This drop is attributed to the ability of the liquid-vapor interface to deform around the helium atom. Because helium has a large, positive solvation free energy, it is thermodynamically favorable for the interface to deform, creating a "valley" where the helium experiences a gas-phase-like environment despite being below the Gibbs dividing surface.
Mechanistic Explanation: The continuum model successfully reproduces the non-Maxwellian distribution. The mechanism is identified as the motion of an underdamped particle in a potential subject to a rapidly varying viscous environment. The PMF acts to repel helium from the liquid, while the friction near the interface is so low (approaching gas-phase values) that the particle cannot relax its velocity back to a thermal equilibrium distribution before escaping.
Temperature Dependence: Using the continuum model, the authors predicted the temperature dependence of the excess kinetic energy from 230 K to 300 K. The model predicts that the excess energy reaches a plateau at lower temperatures where friction increases sharply, and increases with temperature in the range where viscosity is roughly constant, correlating with changes in solvation energy.

Significance and Claims
The paper claims to provide a mechanistic interpretation of trace gas evaporation that resolves the anomaly of super-Maxwellian distributions by explicitly including the spatial variation of friction. The authors assert that the standard PMF-only description is insufficient because it ignores the critical role of dissipative forces near the interface. By demonstrating that the interface deforms around the gas particle, leading to anomalously low friction, the study offers a unified explanation for the excess kinetic energy.

The significance of this work lies in its validation of a reduced description—using a Fokker-Planck equation parameterized by molecularly resolved forces—as a tool that can quantitatively reproduce detailed simulation results. This approach allows for the exploration of conditions, such as varying temperatures, where explicit molecular simulation might be computationally difficult. The authors modestly note that while their classical force field with polarization is an approximation (ignoring potential nuclear quantum effects), the framework provides a robust basis for understanding interfacial mass transfer. The predicted temperature dependence of the excess kinetic energy is presented as a testable hypothesis for future experimental verification.

More like this