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
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.
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
- 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
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.
Referee Report
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)
- 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.
- 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'.
- 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)
- Notation for the PZC and the electric-field coupling in the microkinetic model should be defined consistently between the methods and results sections.
- Figure captions should state the exact metal-site densities used in both the DFT and microkinetic calculations.
Simulated Author's Rebuttal
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
-
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
-
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
-
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
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
axioms (3)
- domain assumption Constrained-magnetization DFT and Landau analysis correctly identify ground-state magnetic moments across densities
- domain assumption Explicit-solvent simulations accurately capture density-dependent PZC shifts and interfacial fields
- domain assumption The pH-field-coupled microkinetic model validly represents ORR kinetics and selectivity
Reference graph
Works this paper leans on
-
[1]
(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]
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]
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, ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.