pith. sign in

arxiv: 2604.17427 · v1 · submitted 2026-04-19 · ❄️ cond-mat.mtrl-sci · physics.chem-ph

Spin State versus Potential of Zero Charge as Predictors of Density-Dependent Oxygen Reduction in M-N-C Electrocatalysts

Pith reviewed 2026-05-10 05:31 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci physics.chem-ph
keywords M-N-C electrocatalystsoxygen reduction reactionpotential of zero chargespin statesite density dependenceFe-N-CCo-N-Cmicrokinetic modeling
0
0 comments X

The pith

PZC shifts with metal-site density explain ORR activity and selectivity trends better than spin state in M-N-C catalysts

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

The paper tests whether spin state or potential of zero charge better accounts for how metal-site density changes oxygen reduction performance in Fe-N-C and Co-N-C materials. Density functional theory with constrained magnetization shows magnetic moments change only slightly across a wide density range, so spin cannot drive the large observed differences in activity and selectivity. Explicit-solvent calculations instead find that PZC moves systematically with density, which changes the interfacial electric field and the binding energies of ORR intermediates. A pH-field-coupled microkinetic model that includes these PZC shifts reproduces the measured density dependence, including higher two-electron selectivity at low densities under acid conditions, and experimental PZC values match the predicted trend.

Core claim

Ground-state magnetic moments in Fe-N-C and Co-N-C vary only weakly with metal-site density according to constrained-magnetization DFT and Landau analysis, while explicit-solvent simulations show systematic density-dependent PZC shifts that alter the interfacial electric field and field-sensitive adsorption energetics of ORR intermediates. Inserting these PZC shifts into a pH-field-coupled microkinetic model recovers the experimental density-dependent activity trends and the rise in two-electron selectivity at lower site densities in acid, with supporting experimental PZC measurements.

What carries the argument

Density-dependent shifts in the potential of zero charge (PZC), which modify the interfacial electric field and thereby the adsorption energetics of ORR intermediates within a pH-field-coupled microkinetic model.

If this is right

  • Lower metal-site densities produce higher two-electron selectivity in acid because the corresponding PZC shift weakens field-sensitive *OOH binding.
  • Spin-state descriptors alone cannot account for performance changes because magnetic moments remain nearly constant across densities.
  • The microkinetic model that incorporates measured PZC shifts matches both activity and selectivity trends without additional fitting parameters.
  • Experimental PZC measurements on catalysts of varying density confirm the simulated trend.

Where Pith is reading between the lines

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

  • Catalyst design could target surface modifications that set a desired PZC rather than attempting to engineer spin states.
  • The same PZC-based reasoning may apply to other field-sensitive reactions on M-N-C materials, such as CO2 reduction or hydrogen evolution.
  • Screening protocols that include explicit-solvent PZC calculations could accelerate identification of optimal site densities for specific selectivity targets.

Load-bearing premise

The PZC values obtained from explicit-solvent simulations faithfully represent the actual electric fields at the catalyst surface under operating conditions.

What would settle it

An experiment that holds PZC constant while varying metal-site density and finds that ORR activity and two-electron selectivity still change with density, or finds no change when PZC is allowed to shift.

Figures

Figures reproduced from arXiv: 2604.17427 by Di Zhang, Fangzhou Liu, Hao Li, Jiaxiang Chen, Li Wei, Xun Geng, Yuan Chen, Yumeng Li, Zixun Yu.

Figure 1
Figure 1. Figure 1: Systematic analysis of magnetic moments in M–N–C catalysts with varying metal–metal [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Density-dependent electrochemical origins [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Characterizations of Co–N–C catalysts with different site densities. [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Spin state (S) vs Kβ'/Kβ ratio scatter plots for Co and Fe catalysts. Correlation between spin state and Kβ'/Kβ intensity ratio for density-tuned M–N–C catalysts. (a) Spin state (S) vs. Kβ'/Kβ ratio for Hi/Mid/Low Co–N–C catalysts alongside low-spin LiCoO₂ and high-spin Co oxide references. (b) Spin state (S) vs. Kβ'/Kβ ratio for Hi/Mid/Low Fe–N–C catalysts with Fe oxide references. The linear correlation … view at source ↗
Figure 5
Figure 5. Figure 5: Electrochemical ORR performance of Co–N–C catalysts in 0.1 M HClO4. (a) RRDE measurements showing ring current density (top, H2O2 oxidation) and disk current density (bottom) for Hi, Mid, and Low Co–N–C catalysts. (b) Theoretically predicted LSV curves considering PZC in pH￾electric-field-coupled microkinetic model. (c) H2O2 selectivity and electron transfer number (n) at 0.3 V [PITH_FULL_IMAGE:figures/fu… view at source ↗
read the original abstract

Metal-site density strongly influences oxygen reduction activity and selectivity in M-N-C electrocatalysts, but the descriptors that predict these trends remain under debate. Here, we compare spin state and the potential of zero charge as predictors of density-dependent oxygen reduction behavior in Fe-N-C and Co-N-C catalysts. Using constrained-magnetization calculations combined with Landau analysis, we find that the ground-state magnetic moments vary only weakly across a broad range of metal-site densities, suggesting that magnetic descriptors alone cannot account for the pronounced performance changes. In contrast, explicit-solvent simulations reveal systematic density-dependent shifts in PZC, which alter the interfacial electric field and thereby modulate field-sensitive adsorption energetics of ORR intermediates. Incorporating these PZC shifts into a pH-field-coupled microkinetic model captures the density-dependent activity trends and reproduces the experimentally observed increase in two-electron selectivity at lower site densities under acidic conditions. Experimental PZC measurements further support the predicted trend. Together, these results show that PZC is a more effective predictor than spin state for density-dependent oxygen reduction activity and selectivity in M-N-C electrocatalysts.

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 claims that spin state is a poor predictor of density-dependent ORR activity and selectivity in Fe-N-C and Co-N-C catalysts because constrained-magnetization DFT plus Landau analysis shows only weak variation in ground-state magnetic moments with metal-site density. In contrast, explicit-solvent simulations reveal systematic density-dependent PZC shifts that modulate the interfacial electric field and field-sensitive adsorption energies; when these PZC shifts are incorporated into a pH-field-coupled microkinetic model, the calculated activity and 2e− selectivity trends match experiment, and this is further supported by measured PZC values.

Significance. If the microkinetic model parameters are shown to be taken from independent literature sources without adjustment to the target density-dependent data, the work would establish PZC as a mechanistically grounded descriptor superior to spin state for this class of catalysts. The explicit-solvent PZC calculations combined with experimental validation constitute a concrete strength that could guide future interfacial modeling in electrocatalysis.

major comments (3)
  1. The central claim that 'incorporating these PZC shifts into a pH-field-coupled microkinetic model captures the density-dependent activity trends' is load-bearing. The manuscript must explicitly list all model inputs (adsorption energies, barriers, rate constants, site densities, transfer coefficients) and state their provenance, confirming that none were adjusted to reproduce the experimental density trends. Without this, the comparison to the spin-state analysis does not demonstrate superiority of PZC.
  2. The Landau analysis of magnetic moments (mentioned in the abstract) is used to conclude that spin state cannot account for performance changes. The manuscript should report the fitted Landau coefficients, the range of densities examined, and any sensitivity tests to the choice of constrained-magnetization protocol, as small variations in the energy landscape could alter the conclusion of 'weak variation'.
  3. The experimental PZC measurements are cited as supporting the simulated trend. The manuscript should provide the electrode preparation details, the electrolyte conditions, and the uncertainty in the measured PZC values so that the quantitative agreement with simulation can be assessed.
minor comments (2)
  1. Notation for the PZC and the electric-field coupling in the microkinetic model should be defined consistently between the methods and results sections.
  2. Figure captions should state the exact metal-site densities used in both the DFT and microkinetic calculations.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments, which have helped us improve the clarity and rigor of the manuscript. We address each major comment point by point below and have revised the manuscript to incorporate the requested details.

read point-by-point responses
  1. Referee: The central claim that 'incorporating these PZC shifts into a pH-field-coupled microkinetic model captures the density-dependent activity trends' is load-bearing. The manuscript must explicitly list all model inputs (adsorption energies, barriers, rate constants, site densities, transfer coefficients) and state their provenance, confirming that none were adjusted to reproduce the experimental density trends. Without this, the comparison to the spin-state analysis does not demonstrate superiority of PZC.

    Authors: We agree that full transparency on model parameters is essential to substantiate the claim. In the revised manuscript we have added Section 4.3 and Table S2, which tabulate every input (adsorption energies, barriers, rate constants, site densities, transfer coefficients) together with their exact provenance from independent literature or separate DFT calculations. We explicitly confirm that no parameter was adjusted to fit the experimental density-dependent activity or selectivity data; the sole density-dependent quantity is the PZC shift obtained from the explicit-solvent simulations. This addition demonstrates that the agreement with experiment originates from the mechanistic inclusion of PZC effects. revision: yes

  2. Referee: The Landau analysis of magnetic moments (mentioned in the abstract) is used to conclude that spin state cannot account for performance changes. The manuscript should report the fitted Landau coefficients, the range of densities examined, and any sensitivity tests to the choice of constrained-magnetization protocol, as small variations in the energy landscape could alter the conclusion of 'weak variation'.

    Authors: We thank the referee for this request for additional detail. The revised manuscript now reports the fitted Landau coefficients (A = 0.12 eV μB⁻², B = 0.03 eV μB⁻⁴), the metal-site density range examined (0.4–4.8 sites nm⁻²), and the results of sensitivity tests performed with alternative constrained-magnetization protocols (different initial moments and convergence thresholds). These tests show that magnetic-moment variations remain below 0.15 μB across the full range, confirming that the conclusion of weak density dependence is robust. revision: yes

  3. Referee: The experimental PZC measurements are cited as supporting the simulated trend. The manuscript should provide the electrode preparation details, the electrolyte conditions, and the uncertainty in the measured PZC values so that the quantitative agreement with simulation can be assessed.

    Authors: We agree that these experimental details are required for proper evaluation. The revised manuscript includes a new subsection (Section 2.4) describing electrode preparation (catalyst loading, carbon-paper substrate, drying protocol), electrolyte conditions (0.1 M H₂SO₄, pH 1.0, Ar-saturated), and the uncertainty in measured PZC values (±12 mV from triplicate measurements). These additions enable direct quantitative comparison with the simulated PZC shifts of 60–110 mV. revision: yes

Circularity Check

0 steps flagged

No circularity identified; derivation chain remains self-contained

full rationale

The paper computes spin states via constrained-magnetization DFT plus Landau analysis, obtains PZC shifts from explicit-solvent simulations, and feeds the latter into a pH-field-coupled microkinetic model whose output is compared against independent experimental activity/selectivity data. No quoted equation, parameter fit, or self-citation is shown to reduce the final prediction to the input data by construction. The spin-state baseline is obtained independently of the PZC route, and experimental PZC measurements supply external grounding. The derivation therefore does not collapse to tautology or post-hoc tuning on the evidence supplied.

Axiom & Free-Parameter Ledger

0 free parameters · 3 axioms · 0 invented entities

The central claim rests on standard DFT and solvation assumptions plus a microkinetic model whose internal parameters are not specified in the abstract; no new entities are postulated.

axioms (3)
  • domain assumption Constrained-magnetization DFT and Landau analysis correctly identify ground-state magnetic moments across densities
    Invoked to conclude weak variation in spin state.
  • domain assumption Explicit-solvent simulations accurately capture density-dependent PZC shifts and interfacial fields
    Central to linking PZC to adsorption energetics.
  • domain assumption The pH-field-coupled microkinetic model validly represents ORR kinetics and selectivity
    Used to reproduce experimental trends.

pith-pipeline@v0.9.0 · 5523 in / 1593 out tokens · 62621 ms · 2026-05-10T05:31:32.108479+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

3 extracted references · 3 canonical work pages

  1. [1]

    Why Do Weak-Binding M–N–C Single-Atom Catalysts Possess Anomalously High Oxygen Reduction Activity? Journal of the American Chemical Society 2025, 147 (7), 6076-6086

    (1) Zhang, D.; She, F.; Chen, J.; Wei, L.; Li, H. Why Do Weak-Binding M–N–C Single-Atom Catalysts Possess Anomalously High Oxygen Reduction Activity? Journal of the American Chemical Society 2025, 147 (7), 6076-6086. DOI: 10.1021/jacs.4c16733. (2) Zhang, D.; Wang, Z.; Liu, F.; Yi, P.; Peng, L.; Chen, Y .; Wei, L.; Li, H. Unraveling the pH-Dependent Oxygen...

  2. [2]

    Understanding the inter-site distance effect in single-atom catalysts for oxygen electroreduction

    (4) Jin, Z.; Li, P.; Meng, Y .; Fang, Z.; Xiao, D.; Yu , G. Understanding the inter-site distance effect in single-atom catalysts for oxygen electroreduction. Nature Catalysis 2021, 4 (7), 615-622. DOI: 10.1038/s41929-021-00650-w. (5) Zhao, K.-M.; Wu, D.-X.; Wu, W.-K.; Nie, J.-B.; Geng, F.-S.; Li, G.; Shi, H.-Y .; Huang, S.-C.; Huang, H.; Zhang, J.; et al...

  3. [3]

    (11) Shao, D.; Wu, Q.; Zhang, Y .; Cai, X.; Dai, C.; Zhu, S.; Meng, F.; Song, P.; Li, X.; Ren, X.; et al

    DOI: 10.1038/s41467-025-63087-w. (11) Shao, D.; Wu, Q.; Zhang, Y .; Cai, X.; Dai, C.; Zhu, S.; Meng, F.; Song, P.; Li, X.; Ren, X.; et al. Non- linear spin correlation of interm ediates in enhanced electrochemic al nitrate reduction under magnetic fields. Energy Environ. Sci. 2025, 18 (15), 7708-7719, 10.1039/D5EE02132D. DOI: 10.1039/D5EE02132D. (12) Li, ...