pith. sign in

arxiv: 2603.17486 · v2 · submitted 2026-03-18 · ❄️ cond-mat.soft · physics.chem-ph

Free-Energy Analysis of Bubble Nucleation on Electrocatalytic Surfaces

Qingguang Xie , Paolo Malgaretti , Othmane Aouane , Simon Thiele , Jens Harting This is my paper

Pith reviewed 2026-05-15 09:01 UTC · model grok-4.3

classification ❄️ cond-mat.soft physics.chem-ph
keywords bubble nucleationfree-energy modelelectrocatalysissupersaturationelectrolysisnucleation barrier
0
0 comments X

The pith

A free-energy model predicts activation energies and critical radii for bubbles on electrocatalyst surfaces.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a free-energy model that calculates the energy barrier and size of the critical nucleus for gas bubbles forming on catalyst surfaces during electrolysis. The model incorporates gas supersaturation, temperature, pressure, and surface wettability to forecast nucleation behavior. It derives power-law scalings where the activation energy falls as the inverse square of supersaturation and the critical radius falls as the inverse of supersaturation. Predictions of critical radii for hydrogen, oxygen, and nitrogen bubbles match experimental values. A coupled diffusion-reaction model further connects supersaturation levels to applied current density.

Core claim

The free-energy model shows that the maximum activation barrier ΔG_max scales as ζ^{-2} with supersaturation ζ while the critical nucleus radius R_c scales as ζ^{-1}. These relations, combined with surface wettability, yield quantitative agreement between predicted critical radii for H2, O2, and N2 bubbles and measured values across experimental conditions.

What carries the argument

Classical free-energy expression for bubble nucleation that balances supersaturation-driven volume energy against surface tension cost, adjusted for contact angle on the catalyst surface.

If this is right

  • Activation energies can be computed directly for any supersaturation, temperature, and wettability to guide electrolyzer operation.
  • Critical bubble sizes decrease inversely with supersaturation, setting thresholds for when bubbles appear at a given current.
  • Coupling gas diffusion to reaction kinetics gives the highest achievable supersaturation at a prescribed current density.
  • Catalyst layer design can target wettability and geometry to raise or lower nucleation barriers for better gas removal.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Operating currents could be chosen to stay below nucleation thresholds predicted by the scaling laws, reducing overpotentials from bubble blocking.
  • If real surfaces show strong heterogeneity, the model would underpredict the spread in nucleation sites and require statistical averaging over local contact angles.
  • The same free-energy framework could extend to other gas-evolving reactions such as chlorine or CO2 reduction by swapping gas properties into the same equations.

Load-bearing premise

Classical nucleation theory remains valid for nanoscale bubbles on real catalyst surfaces with given wettability and no significant non-classical effects or heterogeneity.

What would settle it

A set of measured critical radii for hydrogen bubbles at several supersaturations that deviate from the predicted inverse scaling with supersaturation.

Figures

Figures reproduced from arXiv: 2603.17486 by Jens Harting, Othmane Aouane, Paolo Malgaretti, Qingguang Xie, Simon Thiele.

Figure 1
Figure 1. Figure 1: FIG. 1: Sketch of a system with a mixture of gas molecules and water without a bubble (a) and with a bubble (b). [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: (a) Energy difference ∆ [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: (a) Energy difference ∆ [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: a) Schematic illustration of the anode side of PEMWE. The gas (red circles) diffuses away through the [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
read the original abstract

Bubble nucleation at catalyst surfaces plays a critical role in the operation of electrolyzers. However, achieving controlled bubble nucleation remains challenging due to limited understanding of the underlying mechanisms. Here, we present a free-energy model that quantitatively predicts both the activation energy and critical nucleus size of bubbles at given supersaturation, temperature, pressure, and surface wettability. We find that the activation energy $\Delta G_{max}$ decreases with increasing supersaturation $\zeta$, following a power-law scaling of $\Delta G_{max} \sim \zeta^{-2}$, while the critical nucleus radius $R_c$ scales as $R_c\sim \zeta^{-1}$. Our theoretical predictions for the critical nucleus radius of hydrogen, oxygen and nitrogen bubbles are in quantitative agreement with experimental measurements. Finally, we present a simple model that couples gas diffusion and electrochemical reaction kinetics to determine the maximum gas supersaturation at a given current density. Our results advance the fundamental understanding of bubble nucleation at catalyst surfaces and provide practical guidelines for catalyst layer design to improve the performance of electrolyzers.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper develops a free-energy model based on classical nucleation theory for bubble nucleation on electrocatalytic surfaces. It derives power-law scalings ΔG_max ~ ζ^{-2} for the activation barrier and Rc ~ ζ^{-1} for the critical nucleus radius as functions of supersaturation ζ, incorporating temperature, pressure, and surface wettability via an effective contact angle. The central claim is that predicted Rc values for hydrogen, oxygen, and nitrogen bubbles show quantitative agreement with experimental measurements. A secondary model couples gas diffusion and electrochemical reaction kinetics to predict maximum supersaturation at given current density.

Significance. If the quantitative agreement holds after detailed validation, the work supplies a predictive, parameter-light framework linking supersaturation to nucleation barriers and sizes on catalyst surfaces. The simple power-law scalings and the diffusion-kinetics coupling offer practical design rules for electrolyzer catalyst layers to mitigate bubble-induced losses, advancing both fundamental understanding in soft-matter electrochemistry and device optimization.

major comments (2)
  1. [Abstract] Abstract: the claim of quantitative agreement between predicted Rc and experimental critical radii for H2, O2, and N2 bubbles is load-bearing yet unsupported by any reported contact-angle values, supersaturation ranges, error bars, or data-exclusion criteria, preventing independent verification of the match.
  2. [Free-energy model] Free-energy model (implied continuum functional): the derivation of Rc ~ ζ^{-1} assumes size-independent surface tension, negligible line tension at the three-phase line, and spatially uniform wettability; for the nanoscale regime (Rc typically 10-100 nm) any of these can shift predicted Rc by amounts comparable to experimental scatter, directly undermining the quantitative-agreement claim.
minor comments (1)
  1. [Abstract] The power-law scalings are stated as results but would benefit from an explicit derivation step (e.g., from the standard spherical-cap free-energy expression) to clarify the origin of the exponents.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major point below and have revised the manuscript accordingly where possible.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim of quantitative agreement between predicted Rc and experimental critical radii for H2, O2, and N2 bubbles is load-bearing yet unsupported by any reported contact-angle values, supersaturation ranges, error bars, or data-exclusion criteria, preventing independent verification of the match.

    Authors: We agree that the abstract would benefit from greater specificity to support the central claim. The full manuscript reports the contact angles (drawn from experimental literature for each gas-surface pair), the supersaturation ranges corresponding to the compared experiments, and the quantitative comparisons in the results section, but these were not explicitly summarized in the abstract. In the revision we will update the abstract to state the contact angles used, the supersaturation range, and that agreement holds within the experimental uncertainties reported in the literature. We will also add a compact table in the main text compiling all parameters, predicted versus measured Rc values, and data sources to enable independent verification. revision: yes

  2. Referee: [Free-energy model] Free-energy model (implied continuum functional): the derivation of Rc ~ ζ^{-1} assumes size-independent surface tension, negligible line tension at the three-phase line, and spatially uniform wettability; for the nanoscale regime (Rc typically 10-100 nm) any of these can shift predicted Rc by amounts comparable to experimental scatter, directly undermining the quantitative-agreement claim.

    Authors: The referee correctly notes the standard assumptions of classical nucleation theory on which the scaling Rc ~ ζ^{-1} is derived. These assumptions are widely used for bubble nucleation modeling, and our quantitative comparisons are performed within that framework using effective wettability parameters. We acknowledge that line tension and surface-tension curvature effects can become non-negligible at 10–100 nm scales. In the revised manuscript we have added a short subsection discussing the magnitude of these corrections (estimated <10% for the contact angles and radii in our study, consistent with prior literature) and have included a brief sensitivity analysis showing that the reported agreement remains within experimental scatter even after plausible corrections. The scaling itself is robust under the stated assumptions, which we now state more explicitly. revision: partial

Circularity Check

0 steps flagged

No significant circularity; scalings are standard CNT derivations

full rationale

The paper constructs its free-energy model from the classical nucleation theory continuum functional, yielding the standard power-law scalings ΔG_max ~ ζ^{-2} and Rc ~ ζ^{-1} directly from supersaturation, pressure, temperature and a single effective contact angle. These relations are not redefined in terms of the target experimental Rc values, nor are parameters fitted to the H2/O2/N2 bubble data and then relabeled as predictions. No self-citation chain is invoked to justify uniqueness or to smuggle an ansatz; the quantitative agreement is presented as external validation. The derivation chain therefore remains independent of its own outputs and is self-contained against established CNT benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The model rests on classical nucleation theory applied to surface bubbles and standard physical inputs for supersaturation, temperature, pressure, and contact angle; no new entities are introduced.

axioms (1)
  • domain assumption Classical nucleation theory applies to bubble formation on surfaces
    The free-energy expressions for activation barrier and critical radius follow the standard CNT form used for droplets and bubbles.

pith-pipeline@v0.9.0 · 5494 in / 1162 out tokens · 40056 ms · 2026-05-15T09:01:21.198546+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

42 extracted references · 42 canonical work pages

  1. [1]

    S. Yuan, C. Zhao, X. Cai, L. An, S. Shen, X. Yan, and J. Zhang. Bubble evolution and transport in PEM water electrolysis: Mechanism, impact, and management.Prog. Energy Combust. Sci., 96:101075, 2023

    work page 2023

  2. [2]

    L. Deng, L. Jin, L. Yang, C. Feng, A. Tao, X. Jia, Z. Geng, C. Zhang, X. Cui, and J. Shi. Bubble evolution dynamics in alkaline water electrolysis.eScience, 5:100353, 2025

    work page 2025

  3. [3]

    Vincent and D

    I. Vincent and D. Bessarabov. Low cost hydrogen production by anion exchange membrane electrolysis: A review.Renew. Sustain. Energy Rev., 81:1690–1704, 2018

    work page 2018

  4. [4]

    M. E. Ivanova, R. Peters, M. M¨ uller, S. Haas, M. F. Seidler, G. Mutschke, K. Eckert, P. R¨ ose, S. Calnan, R. Bagacki, R. Schlatmann, C. Grosselindemann, L.-A. Sch¨ afer, N. H. Menzler, A. Weber, R. Van De Krol, F. Liang, F. F. Abdi, S. Brendelberger, N. Neumann, J. Grobbel, M. Roeb, C. Sattler, I. Duran, B. Dietrich, M. E. C. Hofberger, L. Stoppel, N. ...

    work page 2023

  5. [5]

    S. Yuan, C. Zhao, C. Fu, J. Li, Y. Su, R. Xue, S. Shen, J. Yin, X. Yan, and J. Zhang. Discovery of bubble accumulation behavior in catalyst layer of proton exchange membrane water electrolyzer.Int. J. Heat Mass Transf., 227:125552, 2024

    work page 2024

  6. [6]

    J. Mo, Z. Kang, S. T. Retterer, D. A. Cullen, T. J. Toops, J. B. Green, M. M. Mench, and F. Zhang. Discovery of true electrochemical reactions for ultrahigh catalyst mass activity in water splitting.Sci. Adv., 2:e1600690, 2016

    work page 2016

  7. [7]

    J. Mo, Z. Kang, G. Yang, Y. Li, S. T. Retterer, D. A. Cullen, T. J. Toops, G. Bender, B. S. Pivovar, J. B. Green Jr, and F.-Y. Zhang. In situ investigation on ultrafast oxygen evolution reactions of water splitting in proton exchange membrane electrolyzer cells.J. Mater. Chem. A, 5:18469–18475, 2017

    work page 2017

  8. [8]

    Maheshwari, M

    S. Maheshwari, M. van der Hoef, X. Zhang, and D. Lohse. Stability of surface nanobubbles: A molecular dynamics study. Langmuir, 32:11116–11122, 2016

    work page 2016

  9. [9]

    Y. A. P. Sirkin, E. D. Gadea, D. A. Scherlis, and V. Molinero. Mechanisms of nucleation and stationary states of electrochemically generated nanobubbles.J. Am. Chem. Soc., 141:10801–10811, 2019

    work page 2019

  10. [10]

    Zhang, X

    Y. Zhang, X. Zhu, J. A. Wood, and D. Lohse. Threshold current density for diffusion-controlled stability of electrolytic surface nanobubbles.Proc. Natl. Acad. Sci. U.S.A., 121:e2321958121, 2024

    work page 2024

  11. [11]

    S. Zhan, X. Wang, R. Cheng, T. Zhou, W. Zhang, and J. Wang. Investigation of electrolytic hydrogen nanobubbles behavior on heterogeneous wettability surface by using molecular dynamics simulation.Int. J. Hydrog. Energy, 112:160–171, 2025

    work page 2025

  12. [12]

    Volmer and A

    M. Volmer and A. Weber. Keimbildung in ¨ ubers¨ attigten gebilden.Zeitschrift f¨ ur Physikalische Chemie, 119U:277–301, 1926

    work page 1926

  13. [13]

    L. Farkas. Keimbildungsgeschwindigkeit in ¨ ubers¨ attigten d¨ ampfen.Zeitschrift f¨ ur Physikalische Chemie, 125U:236–242, 1927

    work page 1927

  14. [14]

    Kelton and A.L

    K.F. Kelton and A.L. Greer. Beyond the classical theory. In K.F. Kelton and A.L. Greer, editors,Nucleation in Condensed Matter, volume 15 ofPergamon Materials Series, pages 85–123. Pergamon, 2010

    work page 2010

  15. [15]

    S. D. Lubetkin. The fundamentals of bubble evolution.Chem. Soc. Rev., 24:243–250, 1995

    work page 1995

  16. [16]

    K. J. Vachaparambil and K. E. Einarsrud. Explanation of bubble nucleation mechanisms: A gradient theory approach.J. Electrochem. Soc., 165:E504, 2018

    work page 2018

  17. [17]

    Meloni, A

    S. Meloni, A. Giacomello, and C. M. Casciola. Focus article: Theoretical aspects of vapor/gas nucleation at structured surfaces.J. Chem. Phys., 145:211802, 2016

    work page 2016

  18. [18]

    D. W. Oxtoby and R. Evans. Nonclassical nucleation theory for the gas–liquid transition.J. Chem. Phys., 89:7521–7530, 1988

    work page 1988

  19. [19]

    V. K. Shen and P. G. Debenedetti. Density-functional study of homogeneous bubble nucleation in the stretched lennard- jones fluid.J. Chem. Phys., 114:4149–4159, 2001. 9

    work page 2001

  20. [20]

    Zargarzadeh and J

    L. Zargarzadeh and J. A. W. Elliott. Thermodynamics of surface nanobubbles.Langmuir, 32:11309–11320, 2016

    work page 2016

  21. [21]

    L. Lan, Y. Pan, L. Zhao, and B. Wen. Quantitative thermodynamic analyses of nucleation, evolution, and stabilization of surface nanobubbles.Langmuir, 42:5967–5976, 2026

    work page 2026

  22. [22]

    E. D. Gadea, V. Molinero, and D. A. Scherlis. Nanobubble stability and formation on solid–liquid interfaces in open environments.Nano Lett., 23:7206–7212, 2023

    work page 2023

  23. [23]

    M. A. Edwards, H. S. White, and H. Ren. Voltammetric determination of the stochastic formation rate and geometry of individualH 2,N 2, andO 2 bubble nuclei.ACS Nano, 13:6330–6340, 2019

    work page 2019

  24. [24]

    ´A. M. Soto, S. R. German, H. Ren, D. van der Meer, D. Lohse, M. A. Edwards, and H. S. White. The nucleation rate of singleO 2 nanobubbles at Pt nanoelectrodes.Langmuir, 34:7309–7318, 2018

    work page 2018

  25. [25]

    Solymosi, M

    T. Solymosi, M. Geißelbrecht, S. Mayer, M. Auer, P. Leicht, M. Terlinden, P. Malgaretti, A. B¨ osmann, P. Preuster, J. Harting, M. Thommes, N. Vogel, and P. Wasserscheid. Nucleation as a rate-determining step in catalytic gas generation reactions from liquid phase systems.Sci. Adv., 8:eade3262, 2022

    work page 2022

  26. [26]

    B. T. Sangtam and H. Park. Review on bubble dynamics in proton exchange membrane water electrolysis: Towards optimal green hydrogen yield.Micromachines, 14:2234, 2023

    work page 2023

  27. [27]

    Selamet, U

    O.F. Selamet, U. Pasaogullari, D. Spernjak, D.S. Hussey, D.L. Jacobson, and M.D. Mat. Two-phase flow in a proton exchange membrane electrolyzer visualized in situ by simultaneous neutron radiography and optical imaging.Int. J. Hydrogen Energy, 38:5823–5835, 2013

    work page 2013

  28. [28]

    Jamnongwong, K

    M. Jamnongwong, K. Loubiere, N. Dietrich, and G. H´ ebrard. Experimental study of oxygen diffusion coefficients in clean water containing salt, glucose or surfactant: Consequences on the liquid-side mass transfer coefficients.Chem. Eng. J., 165:758–768, 2010

    work page 2010

  29. [29]

    J. Mo, Z. Kang, G. Yang, S. T. Retterer, D. A. Cullen, T. J. Toops, J. B. Green, and F.-Y. Zhang. Thin liquid/gas diffusion layers for high-efficiency hydrogen production from water splitting.Applied Energy, 177:817–822, 2016

    work page 2016

  30. [30]

    Shiva Kumar and V

    S. Shiva Kumar and V. Himabindu. Hydrogen production by PEM water electrolysis – A review.Mater. Sci. Energy Technol., 2:442–454, 2019

    work page 2019

  31. [31]

    Sakuma, Y

    G. Sakuma, Y. Fukunaka, and H. Matsushima. Nucleation and growth of electrolytic gas bubbles under microgravity.Int. J. Hydrogen Energy, 39:7638–7645, 2014

    work page 2014

  32. [32]

    Matsushima, D

    H. Matsushima, D. Kiuchi, and Y. Fukunaka. Measurement of dissolved hydrogen supersaturation during water electrolysis in a magnetic field.Electrochimica Acta, 54:5858–5862, 2009

    work page 2009

  33. [33]

    X. Zhao, H. Ren, and L. Luo. Gas bubbles in electrochemical gas evolution reactions.Langmuir, 35:5392–5408, 2019

    work page 2019

  34. [34]

    H. Lv, S. Wang, Y. Sun, J. Chen, W. Zhou, and C. Zhang. Anode catalyst layer with hierarchical pore size distribution for highly efficient proton exchange membrane water electrolysis.J. Power Sources, 564:232878, 2023

    work page 2023

  35. [35]

    Bernt and H

    M. Bernt and H. A. Gasteiger. Influence of ionomer content in IrO 2 /TiO2 electrodes on PEM water electrolyzer perfor- mance.J. Electrochem. Soc., 163:F3179–F3189, 2016

    work page 2016

  36. [36]

    T.-C. Ma, A. Hutzler, R. Hanke-Rauschenbach, and S. Thiele. Anode catalyst layer optimization in polymer electrolyte membrane water electrolysis: Modeling catalyst layer properties and interface effects.Electrochem. Commun., 178:107965, 2025

    work page 2025

  37. [37]

    Bernt, C

    M. Bernt, C. Schramm, J. Schr¨ oter, C. Gebauer, J. Byrknes, C. Eickes, and H. A. Gasteiger. Effect of the IrOx conductivity on the anode electrode/porous transport layer interfacial resistance in PEM water electrolyzers.J. Electrochem. Soc., 168:084513, 2021

    work page 2021

  38. [38]

    A. O. Maksimov, A. M. Kaverin, and V. G. Baidakov. Heterogeneous vapor bubble nucleation on a rough surface.Langmuir, 29:3924–3934, 2013

    work page 2013

  39. [39]

    Jones, G.M

    S.F. Jones, G.M. Evans, and K.P. Galvin. Bubble nucleation from gas cavities — a review.Adv. Colloid Interface Sci., 80:27–50, 1999

    work page 1999

  40. [40]

    A. A. Atchley and A. Prosperetti. The crevice model of bubble nucleation.J. Acoust. Soc. Am., 86:1065–1084, 1989

    work page 1989

  41. [41]

    L. S. Bartell. Tolman’sδ, surface curvature, compressibility effects, and the free energy of drops.J. Phys. Chem. B, 105:11615–11618, 2001

    work page 2001

  42. [42]

    S. R. German, M. A. Edwards, Q. Chen, and H. S. White. Laplace pressure of individual H 2 nanobubbles from pressure– addition electrochemistry.Nano Lett., 16:6691–6694, 2016

    work page 2016