pith. sign in

arxiv: 2604.21773 · v1 · submitted 2026-04-23 · ⚛️ physics.chem-ph

Molecular dynamics simulations of Nafion thin films at a platinum catalyst surface: Correlating structure with charging behaviour

Dustin Vivod , Binny A. Davis , Tobias Binninger , Michael Eikerling This is my paper

Pith reviewed 2026-05-08 13:30 UTC · model grok-4.3

classification ⚛️ physics.chem-ph
keywords Nafionplatinum catalystmolecular dynamicsthin filmsdifferential capacitancehydronium ionsPEM fuel cellsVoronoi tessellation
0
0 comments X

The pith

Molecular dynamics simulations show water films under Nafion on platinum remain stable below 1.3 nm and tie charging behaviour to surface crowding by ions or polymer.

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

The paper builds atom-scale models of a platinum surface coated with Nafion and controlled water layers to examine the local environment that governs catalysis in fuel cells. Molecular dynamics runs with a Voronoi-based film assembly identify that water configurations thinner than 1.3 nm are energetically stable. When the platinum surface is charged, the simulations track shifts in electrostatic potential and differential capacitance, with the observed trends linked to accumulation of hydronium ions or Nafion segments directly at the metal. This structural insight matters because the immediate surroundings of the catalyst control reaction rates and ion transport in operating devices. The authors frame the entire protocol as reusable for testing alternative ionomer materials.

Core claim

Through construction of a dense Nafion film on platinum via Voronoi tessellation and molecular dynamics simulations, water film configurations with thickness less than 1.3 nm are found to be stable. Simulations of charged platinum surfaces show that electrostatic conditions and differential capacitance of the interface can be interpreted through crowding of hydronium ions or the Nafion film at the platinum surface. The workflow thereby connects molecular arrangement directly to charging behaviour at the catalyst-ionomer boundary.

What carries the argument

Voronoi tessellation construction that assembles a dense Nafion film fully covering the platinum substrate, enabling controlled addition of water layers for stability and charging analysis.

If this is right

  • Water films thinner than 1.3 nm remain energetically stable, indicating that limited hydration does not destabilize the Nafion-platinum interface.
  • Trends in differential capacitance and electrostatics under charging arise from whether hydronium ions or Nafion segments crowd the platinum surface.
  • The atomistic model directly correlates molecular structure with macroscopic electrical response at the catalyst surface.
  • The simulation protocol applies without major modification to other ionomer materials for PEMFC studies.

Where Pith is reading between the lines

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

  • The 1.3 nm stability threshold may help set ionomer thickness targets that balance hydration and avoid flooding or dry-out at the catalyst in devices.
  • Running the same models across a range of applied potentials could produce capacitance-voltage curves for direct comparison with electrochemical experiments.
  • Ionomer modifications that alter chain flexibility or ion distribution could be screened to reduce unwanted crowding and tune interface capacitance.
  • Repeating the setup on stepped or defective platinum surfaces would test whether the crowding mechanism holds for realistic catalyst particles.

Load-bearing premise

The Voronoi tessellation produces a realistic dense Nafion film that fully covers the platinum without artifacts from the assembly protocol or force-field inaccuracies at the interface.

What would settle it

Energy minimization calculations that find lower energies for water films thicker than 1.3 nm, or experimental capacitance data that fail to show the predicted dependence on hydronium versus Nafion crowding, would contradict the stability and interpretation results.

Figures

Figures reproduced from arXiv: 2604.21773 by Binny A. Davis, Dustin Vivod, Michael Eikerling, Tobias Binninger.

Figure 1
Figure 1. Figure 1: Chemical structure of Nafion with backbone and side-chain parts highlighted view at source ↗
Figure 2
Figure 2. Figure 2: Exemplary snapshots of the two initial Nafion layers used within this work. view at source ↗
Figure 3
Figure 3. Figure 3: Exemplary snapshot of the dense Nafion layer. Hydronium ions were added for view at source ↗
Figure 4
Figure 4. Figure 4: Snapshots during production runs of systems with different water content. (a): view at source ↗
Figure 5
Figure 5. Figure 5: Density distributions of water oxygen from different levels of theory.[43, 44] view at source ↗
Figure 6
Figure 6. Figure 6: Distributions of water OH bond angles to surface normal from this study (blue) view at source ↗
Figure 7
Figure 7. Figure 7: Densities along Z direction from the Pt surface for selected atom species (col view at source ↗
Figure 8
Figure 8. Figure 8: Side-chain angles histograms of the ‘dense’ layer systems with different water view at source ↗
Figure 9
Figure 9. Figure 9: Comparison of charge densities (left) and electrostatic potentials (right) be view at source ↗
Figure 10
Figure 10. Figure 10: Centers of molecular excess charge (cf. Eq. 7) per water content with respect view at source ↗
Figure 11
Figure 11. Figure 11: (a) Electric surface potentials ϕs as a function of surface charge for each water content. The potential reference is chosen such that the uncharged systems have ϕs=0. (b) Differential capacitances Cdiff as a function of electric surface potential for each water content. Calculated using Eq. 4 using neighbouring charge steps and ϕs values as described in Sec. 2.2. the presence of the second adsorption lay… view at source ↗
Figure 12
Figure 12. Figure 12: Comparison of Z densities between uncharged and charged systems of hydro view at source ↗
Figure 13
Figure 13. Figure 13: Temporal evolution of the number of atoms that are within 4 view at source ↗
Figure 14
Figure 14. Figure 14: Visual example of side-chain angles w.r.t. the platinum surface. view at source ↗
Figure 15
Figure 15. Figure 15: Exemplary Voronoi diagrams for the square and hex Nafion layers (col view at source ↗
Figure 16
Figure 16. Figure 16: Exemplary Voronoi diagrams for the square and hex Nafion layers (column view at source ↗
Figure 17
Figure 17. Figure 17: Exemplary Voronoi diagrams for the square and hex Nafion layers (column view at source ↗
Figure 18
Figure 18. Figure 18: Exemplary Voronoi diagrams for the square and hex Nafion layers (column view at source ↗
Figure 19
Figure 19. Figure 19: Exemplary Voronoi diagrams for the dense Nafion layers at both investigated view at source ↗
Figure 20
Figure 20. Figure 20: Left: Temporal evolution of the potentials at the platinum surface during the view at source ↗
Figure 21
Figure 21. Figure 21: Comparison of charge densities (left) and potentials (right) between un view at source ↗
Figure 22
Figure 22. Figure 22: Distribution of water dipole angles w.r.t. the surface normal for the charged view at source ↗
read the original abstract

Electrocatalysis is greatly influenced by the local reaction environment, which is governed by the structure of the catalyst, the distribution of the electrolyte, and the local electric field. In catalytic systems comprised of complex molecular species like ionomers, the distribution of electrolyte can vary substantially, resulting in divers local reaction environments. In order to gain atom-scale insight into this micro-environment we construct a model system consisting of a platinum surface, varying levels of water, and a Nafion thin film and conduct molecular dynamics simulations. We employ a construction based on Voronoi tesselation to assemble a dense film of ionomer that fully covers the platinum substrate. An energy analysis reveals that water film configurations with thickness of less then 1.3 nm are stable. Simulations with charged platinum surfaces are analysed in view of electrostatic conditions and differential capacitance of the interface configuration. Trends observed in these properties can be interpreted in view of the crowding of hydronium ions or the Nafion film at the platinum surface. The presented workflow can be easily applied to investigate novel ionomers for use in PEMFCs.

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

3 major / 2 minor

Summary. The manuscript presents molecular dynamics simulations of Nafion thin films on a platinum catalyst surface, using a Voronoi tessellation construction to assemble a dense ionomer film that fully covers the substrate. An energy analysis indicates that water film configurations thinner than 1.3 nm are stable. Simulations with charged platinum surfaces are used to examine electrostatic conditions and differential capacitance, with trends interpreted in terms of hydronium ion or Nafion crowding at the interface. The authors propose that the workflow can be applied to investigate novel ionomers for PEMFCs.

Significance. If the results hold after addressing validation gaps, this work would provide useful atomistic insights into the local structure and charging behavior at complex catalyst-ionomer interfaces in fuel cells. The MD approach to multi-component thin films and the suggested workflow for screening ionomers represent a constructive contribution to understanding electrocatalytic microenvironments.

major comments (3)
  1. [Methods (Voronoi tessellation)] Methods section describing the Voronoi tessellation construction: No quantitative validation of the assembled Nafion film is provided (e.g., mass density profiles perpendicular to the Pt surface, chain end-to-end distances, or sulfonate distributions compared to bulk Nafion or experimental thin-film data). This is load-bearing for the central claims, as the reported water-film stability threshold and capacitance trends via crowding both assume the film is a realistic, dense, artifact-free layer.
  2. [Energy analysis of water films] Results section on energy analysis: The specific claim that water films with thickness less than 1.3 nm are stable is stated without error bars, raw energy data, force-field details, or equilibration protocols, preventing assessment of the threshold's robustness.
  3. [Charged platinum surfaces and capacitance] Results section on charged surfaces: The analysis of differential capacitance and its interpretation in terms of hydronium/Nafion crowding lacks details on the capacitance calculation method from the trajectories and reports no statistical uncertainties, weakening the link between structure and electrostatic properties.
minor comments (2)
  1. [Abstract] Abstract: 'less then 1.3 nm' should read 'less than 1.3 nm'; 'divers local reaction environments' should read 'diverse local reaction environments'.
  2. [Methods] The manuscript would benefit from explicit statements of the force field, thermostat/barostat settings, and simulation lengths to improve reproducibility.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their careful and constructive review of our manuscript. We agree that additional quantitative validation, methodological details, and uncertainty reporting will improve the clarity and robustness of our findings. We address each major comment below and will incorporate the suggested revisions.

read point-by-point responses

  1. Referee: Methods section describing the Voronoi tessellation construction: No quantitative validation of the assembled Nafion film is provided (e.g., mass density profiles perpendicular to the Pt surface, chain end-to-end distances, or sulfonate distributions compared to bulk Nafion or experimental thin-film data). This is load-bearing for the central claims, as the reported water-film stability threshold and capacitance trends via crowding both assume the film is a realistic, dense, artifact-free layer.

    Authors: We agree that quantitative validation of the Nafion film is necessary to substantiate the model. In the revised manuscript, we will add mass density profiles perpendicular to the Pt surface, comparisons of chain end-to-end distances and sulfonate distributions against bulk Nafion reference simulations, and discussion of consistency with experimental thin-film data. These additions will confirm that the Voronoi-based assembly yields a dense, realistic layer without significant artifacts. revision: yes

  2. Referee: Results section on energy analysis: The specific claim that water films with thickness less than 1.3 nm are stable is stated without error bars, raw energy data, force-field details, or equilibration protocols, preventing assessment of the threshold's robustness.

    Authors: We acknowledge that the energy analysis requires more supporting information. The revised version will report error bars on the energy values (from block averaging or independent runs), provide full details on the force field (including Pt-interface parameters), and describe the equilibration protocols. Raw energy data will be included in the supplementary information to allow assessment of the 1.3 nm threshold. revision: yes

  3. Referee: Results section on charged surfaces: The analysis of differential capacitance and its interpretation in terms of hydronium/Nafion crowding lacks details on the capacitance calculation method from the trajectories and reports no statistical uncertainties, weakening the link between structure and electrostatic properties.

    Authors: We will expand the description of the differential capacitance calculation, specifying the method used to extract it from the trajectories (including relevant equations and assumptions about constant potential or charge). Statistical uncertainties will be added using appropriate averaging techniques. These changes will better support the structural interpretations regarding hydronium and Nafion crowding. revision: yes

Circularity Check

0 steps flagged

No circularity: results are direct MD outputs

full rationale

The paper constructs a Nafion film via Voronoi tessellation, runs molecular dynamics trajectories, and reports energy minima for water films <1.3 nm plus capacitance trends from electrostatic analysis. These quantities are computed outputs of the simulations rather than quantities redefined or fitted to match prior inputs within the paper. No self-citations, ansatzes, or uniqueness theorems are invoked as load-bearing steps, and no equations reduce the reported stability or crowding interpretations to the assembly protocol by construction. The derivation chain is therefore self-contained.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claims depend on classical molecular dynamics assumptions and a custom film-assembly protocol whose validity is not independently verified in the abstract.

free parameters (2)
  • water film thickness threshold
    Identified via energy analysis as less than 1.3 nm for stability; treated as an observed outcome rather than an input fit.
  • Nafion coverage density
    Set by the Voronoi tessellation procedure to achieve full substrate coverage.
axioms (2)
  • domain assumption Classical force fields sufficiently capture Pt-Nafion-water interactions and charging behavior
    Invoked throughout the MD simulations without mention of quantum corrections or experimental validation.
  • ad hoc to paper Voronoi tessellation yields a physically realistic dense ionomer film
    The construction method is introduced to assemble the model system.

pith-pipeline@v0.9.0 · 5503 in / 1372 out tokens · 69466 ms · 2026-05-08T13:30:45.891810+00:00 · methodology

discussion (0)

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