How much energy does it take to open a bubble in stretched liquid?
Averaging the work from many fast, non-equilibrium simulations recovers the full energy barrier for opening a nanoscale cavity in a stretched liquid, even under mild conditions that ordinary simulations cannot reach.
Source
Free Energy Evaluation of Cavity Formation in Metastable Liquid Based on Stochastic Thermodynamics
Study at a glance
- Design
- Computational / modelling — Molecular dynamics of a Lennard-Jones liquid in which an expanding repulsive field carves out a cavity; the work done in 100 repeated runs is averaged with the Jarzynski equality to obtain free energy versus cavity radius at four negative pressures.
- N
- No participants; each pressure condition used 100 independent simulations from different equilibrium starting configurations.
- Population
- Simulated Lennard-Jones (argon-like) liquid under negative pressure at reduced temperature 0.8
- Outcome
- Free-energy change versus cavity radius; critical cavity size and activation energy; fitted surface tension and bulk free-energy difference
Structured fields used in claim comparison tables when every cited study has a complete layer.
What they did
The authors simulated a simple Lennard-Jones liquid held at negative pressure, where it is metastable and could boil or cavitate. A repulsive spherical field was grown steadily at the centre of the box to push atoms out and make a cavity, and the work done on the atoms was recorded. Repeating this 100 times per condition, they used the Jarzynski equality, which relates an exponential average of non-equilibrium work to the equilibrium free-energy difference, to obtain free energy as a function of cavity radius at four pressures.
What they found
Each pressure gave a smooth free-energy curve with a maximum, the critical nucleus, as predicted by classical nucleation theory; in one case the critical radius was about 6.5 atomic diameters with an activation energy of 190 plus or minus 10 energy units. The barrier and critical size both shrank as the liquid was stretched more, and at the mildest pressure the barrier of about 400 units implies a probability too tiny for direct simulation. Fitting the curves to the classical cubic-plus-square form gave a surface tension that rose with the degree of non-equilibrium.
The limits
What it doesn't show
The cavity is a vacuum forced open by an external field, not a real vapour bubble that nucleates spontaneously, so the results describe cavity formation and only approximately bubble nucleation. The authors warn that 100 samples may be too few, because very few runs had work below the estimated free energy and more sampling could lower the barrier. The method assumes a spherical cavity, limiting it to low temperatures, and the expansion speed, although tested, was still fast. Only one simple model liquid was studied.
Key terms
- Metastable liquid
- A liquid held beyond its phase boundary, here under tension, that has not yet transformed because forming the new phase requires crossing an energy barrier.
- Classical nucleation theory (CNT)
- A model in which the free energy of a bubble is a negative volume term plus a positive surface-tension term, giving a barrier at a critical radius.
- Jarzynski equality
- An exact result stating that the exponential average of work done in many non-equilibrium processes equals the exponential of the equilibrium free-energy difference.
- Critical nucleus
- The bubble size at the top of the free-energy barrier; smaller bubbles tend to shrink and larger ones grow spontaneously.
- Lennard-Jones potential
- A simple pair interaction with short-range repulsion and weaker long-range attraction, commonly used to model argon-like atoms.
- Binodal and spinodal
- The binodal marks where two phases coexist in equilibrium; the spinodal marks where the metastable phase becomes unstable and transforms without a barrier.
Flashcards
0 of 10 answers reviewed
Research intelligence for this paper
See its role on concept claims, tensions it is part of, placement history, and related discoveries.
Quiz yourself
What does the Jarzynski equality relate?
Common questions
Why not simply wait for bubbles to form in an ordinary simulation?
Near the binodal the barrier is so high that spontaneous nucleation would essentially never occur in a simulation's time span; previous studies therefore worked near the spinodal.
Why is an exponential average of work needed instead of the plain average?
Fast processes dissipate energy, so the mean work overestimates the free-energy change; the exponential weighting emphasises rare low-work runs and removes that bias.
How does this relate to real vapour bubbles?
A correction for vapour pressure inside the bubble is small at this low temperature, so the cavity results approximate bubble energetics, but the simulation itself creates vacuum.
More on Phase transitions