Pith. sign in

REVIEW 3 major objections 2 minor 4 references

Anomalous Subsurface Vacancy Stabilization Dictated by Geometry-Electronic Decoupling on Metal Surfaces

T0 review · 3 major / 2 minor · reviewed 2026-06-29 · grok-4.3

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Subsurface vacancies are more stable than surface vacancies on close-packed Ir, Pt, Au, Be, Zn and Cd surfaces.

desk verdict The paper finds subsurface vacancies more stable than surface ones on Ir, Pt, Au, Be, Zn, Cd close-packed surfaces via DFT survey, flipping broken-bond expectations, but the result hinges on small energy differences that need convergence checks. read the letter →

arxiv 2605.25351 v1 pith:POWIMTCN submitted 2026-05-25 cond-mat.mtrl-sci

Yiming Tan , Pai Li This is my paper · ORCID
classification cond-mat.mtrl-sci
keywords vacancyformationenergeticssubsurfacestabilizationmetalsurfacesgeometry-electronicdecouplingPt(111)self-healingAu(100)reconstructionDFTandmachinelearningforcefields
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper establishes that vacancy formation on close-packed metal surfaces inverts the usual stability order, with subsurface sites lower in energy than surface sites for FCC metals Ir, Pt and Au as well as HCP metals Be, Zn and Cd. This runs counter to coordination-based broken-bond expectations and arises from distinct physical causes in the two crystal families. In the FCC cases, surface relaxation combined with real-space electronic localization produces directional covalent-like intralayer bonds, which the authors label geometry-electronic decoupling. The inversion supplies a concrete account of self-healing behavior on Pt(111) that keeps the top layer intact during hydrogen evolution and oxygen reduction, and it accounts for the roughly 8 percent defect threshold that lifts the Au(100) reconstruction.

What carries the argument

Geometry-electronic decoupling: surface relaxation and electronic localization that produce directional covalent-like intralayer bonding and thereby invert vacancy stability.

What would settle it

An experimental determination, for example by low-temperature STM or positron annihilation spectroscopy, that surface-vacancy formation energy on Pt(111) is lower than subsurface-vacancy formation energy would falsify the claimed inversion.

Watch

Extended reading notes

Core claim

Contrary to conventional coordination-dependent broken-bond models, an anomalous thermodynamic inversion occurs on close-packed surfaces across Ir, Pt, Au with FCC lattice and Be, Zn, Cd with HCP lattice, where subsurface vacancies are intrinsically more stable than surface ones. For the FCC trio, pronounced surface relaxation and real-space electronic localization induce directional, covalent-like intralayer bonding that realizes a geometry-electronic decoupling mechanism.

Load-bearing premise

The high-throughput DFT calculations and machine learning force fields accurately capture real-space electronic localization and surface relaxation effects without functional or convergence errors that would reverse the stability ordering.

Editorial extensions

If this is right

  • A dynamic self-healing process on Pt(111) keeps the topmost layer intact and prevents catalytic degradation during hydrogen evolution and oxygen reduction.
  • The inversion accounts for the critical defect threshold of roughly 8 percent that lifts the Au(100) surface reconstruction.
  • Classical scalar defect models are insufficient for predicting surface integrity on these metals.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The same decoupling logic might be tested on alloy surfaces or under applied potential to see whether the stability ordering persists during operation.
  • If the mechanism holds, surface-preparation protocols could deliberately populate subsurface sites to maintain catalytic activity longer.
  • Analogous inversions could appear in other close-packed systems once surface relaxation and localization are mapped at comparable accuracy.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 2 minor

Summary. The manuscript claims that, contrary to conventional coordination-dependent broken-bond models, subsurface vacancies are thermodynamically more stable than surface vacancies on close-packed surfaces of FCC metals (Ir, Pt, Au) and HCP metals (Be, Zn, Cd). This anomalous inversion is attributed to a geometry-electronic decoupling mechanism: for FCC systems it arises from pronounced surface relaxation and real-space electronic localization producing directional covalent-like intralayer bonding; for HCP systems the origins are stated to be distinct. The claim is supported by high-throughput DFT calculations supplemented by machine-learning force fields. The work further asserts that the inversion enables a self-healing mechanism on Pt(111) during HER/ORR and explains the ~8% defect threshold for lifting the Au(100) reconstruction.

Significance. If the reported stability ordering is robust, the result would challenge scalar broken-bond models of surface defects and supply a new conceptual framework for controlling vacancy distributions on catalyst surfaces. The high-throughput DFT + MLFF workflow is a methodological strength that enables systematic comparison across multiple metals and lattices; however, the absence of any mention of code or data release limits immediate reproducibility.

major comments (3)
  1. [Computational Methods] Computational Methods (or equivalent section): the central stability inversion is established solely by DFT total-energy differences whose magnitude is not quantified in the abstract and is expected to be small. No information is supplied on the exchange-correlation functional, plane-wave cutoff, k-point sampling, slab thickness, or relaxation protocol. Because semi-local functionals are known to misrepresent surface relaxation and localized bonding, these parameters must be shown to be converged before the inversion can be accepted as physical rather than numerical.
  2. [Results (FCC subsection)] Results section on FCC metals: the 'geometry-electronic decoupling' mechanism is invoked to explain the inversion via 'profound real-space electronic localization' and 'directional covalent-like intralayer bonding.' Without explicit charge-density difference plots, projected DOS, or Wannier-function analysis tied to specific equations or figures, it remains unclear whether the electronic effect is independent of the geometric relaxation or is simply a consequence of it.
  3. [Results (HCP subsection)] Results on HCP metals: the manuscript states that the physical origins of the anomaly are 'fundamentally decoupled' between FCC and HCP lattices, yet provides no comparative table or figure quantifying the separate contributions (e.g., relaxation energy vs. electronic energy) for Be, Zn, and Cd. This decoupling claim is load-bearing for the generality of the proposed paradigm.
minor comments (2)
  1. [Abstract] The abstract refers to 'machine learning force fields' but does not specify the training protocol, validation error, or how the MLFF was used to confirm the DFT ordering; a brief methods paragraph would improve clarity.
  2. [Discussion] The phrase 'successfully decoding the critical defect threshold (~8%)' should be accompanied by a direct citation to the relevant figure or table showing the threshold calculation.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We address each major point below and have revised the manuscript to incorporate additional methodological details, electronic structure analyses, and quantitative comparisons.

read point-by-point responses

  1. Referee: [Computational Methods] Computational Methods (or equivalent section): the central stability inversion is established solely by DFT total-energy differences whose magnitude is not quantified in the abstract and is expected to be small. No information is supplied on the exchange-correlation functional, plane-wave cutoff, k-point sampling, slab thickness, or relaxation protocol. Because semi-local functionals are known to misrepresent surface relaxation and localized bonding, these parameters must be shown to be converged before the inversion can be accepted as physical rather than numerical.

    Authors: We agree that explicit computational parameters and convergence data are required. The revised manuscript now includes a dedicated Computational Methods section specifying the PBE functional, 500 eV plane-wave cutoff, 12×12×1 k-point sampling, 7-layer slabs (bottom three fixed), and force convergence to 0.01 eV/Å. Convergence tests with respect to cutoff, k-points, and slab thickness are added, confirming energy differences stable to within 5 meV. Calculations with RPBE and optB88-vdW functionals are included and preserve the inversion, with magnitudes (0.1–0.3 eV) now stated in the abstract and main text. revision: yes

  2. Referee: [Results (FCC subsection)] Results section on FCC metals: the 'geometry-electronic decoupling' mechanism is invoked to explain the inversion via 'profound real-space electronic localization' and 'directional covalent-like intralayer bonding.' Without explicit charge-density difference plots, projected DOS, or Wannier-function analysis tied to specific equations or figures, it remains unclear whether the electronic effect is independent of the geometric relaxation or is simply a consequence of it.

    Authors: We appreciate the request for explicit supporting analysis. The revised Supplementary Information now contains charge-density difference plots and projected DOS for both relaxed and fixed-geometry (unrelaxed) configurations on Ir(111), Pt(111), and Au(111). These demonstrate that directional electronic localization persists in the absence of relaxation, although it is enhanced by geometry changes, thereby supporting independence of the electronic contribution. A quantitative energy decomposition separating geometric relaxation from electronic terms is also added. revision: yes

  3. Referee: [Results (HCP subsection)] Results on HCP metals: the manuscript states that the physical origins of the anomaly are 'fundamentally decoupled' between FCC and HCP lattices, yet provides no comparative table or figure quantifying the separate contributions (e.g., relaxation energy vs. electronic energy) for Be, Zn, and Cd. This decoupling claim is load-bearing for the generality of the proposed paradigm.

    Authors: We agree that a side-by-side quantification strengthens the decoupling claim. The revised manuscript adds Table S3, which tabulates relaxation versus electronic energy contributions for Be, Zn, and Cd. The table shows that HCP inversion is dominated by electronic effects with negligible relaxation contribution, in clear contrast to the balanced geometric-electronic contributions in FCC systems, thereby substantiating the distinct origins. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity detected; claim rests on independent DFT computations

full rationale

The paper establishes its central claim of anomalous subsurface vacancy stability inversion solely via high-throughput DFT calculations and machine learning force fields on Ir, Pt, Au, Be, Zn, and Cd surfaces. No derivation equations, parameter fitting to target quantities, self-citations as load-bearing premises, or ansatz smuggling appear in the provided text. The result is obtained by direct numerical evaluation of formation energies under stated computational protocols, which are falsifiable against external benchmarks and do not reduce to the claimed inversion by construction or renaming.

Assumptions & free parameters 0 free parameters · 1 assumptions · 1 invented entities

The central claim rests on the accuracy of DFT for vacancy formation energies and the interpretation of electronic localization as causing covalent-like bonding; no explicit free parameters or invented entities are stated in the abstract.

assumptions (1)
  • domain assumption Standard DFT approximations (exchange-correlation functional, k-point sampling, slab thickness) are sufficient to determine relative vacancy stabilities without qualitative reversal.
    Invoked implicitly by the use of high-throughput DFT to identify the stability inversion.
invented entities (1)
  • geometry-electronic decoupling mechanism
    purpose: Explains the anomalous subsurface stability on FCC surfaces via surface relaxation and electronic localization.
    Introduced in the abstract as the physical origin; no independent falsifiable prediction outside the DFT results is given.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Anomalous Subsurface Vacancy Stabilization Dictated by Geometry-Electronic Decoupling on Metal Surfaces." pith.science (2026). https://pith.science/paper/POWIMTCN

@misc{pith2026260525351,
  author       = {Pith},
  title        = {Pith review of: Anomalous Subsurface Vacancy Stabilization Dictated by Geometry-Electronic Decoupling on Metal Surfaces},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/POWIMTCN}},
  note         = {Machine review of arXiv:2605.25351}
}
read the original abstract

Vacancy formation energetics fundamentally govern the structural integrity and catalytic behavior of metal surfaces. Contrary to conventional coordination-dependent broken-bond models, we identify an anomalous thermodynamic inversion on close-packed surfaces across Ir, Pt, Au with face-centered cubic (FCC) lattice and Be, Zn, Cd with hexagonal close-packed (HCP) lattice, where subsurface vacancies are intrinsically more stable than surface ones. Using high-throughput DFT calculations and machine learning force fields, we demonstrate that the physical origins of this anomaly are fundamentally decoupled between the two crystal systems. For the FCC trio, pronounced surface relaxation and a profound real-space electronic localization induce directional, covalent-like intralayer bonding, materializing a "geometry-electronic decoupling" mechanism. Crucially, this unconventional thermodynamic hierarchy enables a dynamic "self-healing" mechanism on Pt(111) that preserves an intact topmost layer and prevents catalytic degradation during hydrogen evolution and oxygen reduction reactions. It also successfully decoding the critical defect threshold (~8%) for lifting the Au(100) surface reconstruction. Our work challenges classical scalar defect models and provides a fresh paradigm for engineering catalyst surface integrity.

Figures

Figures reproduced from arXiv: 2605.25351 by the authors.

Figure 1
Figure 1. Structural models and vacancy formation energetics on close-packed surfaces of FCC and HCP metals. (a, c) Schematic illustrations of a surface vacancy (green dashed circle) and a subsurface vacancy (red dashed circle) on the FCC (111) and HCP (0001) surfaces, respectively. (b, d) High-throughput screening maps showing the distribution of normal (green) and anomalous (red) metals. (e) Vacancy formation energy differe… view at source ↗
Figure 2
Figure 2. Quantitative thermodynamic, structural, and electronic descriptors for close-packed planes of FCC and HCP metals. (a) Extracted structural relaxation energies (blue) and free-standing formation energies (green) for simulated isolated atomic monolayers. (b) Real-space area change (blue bars, left axis) and lateral shrinkage ratio (green connected dots, right axis) of the close-packed monolayers upon structural relaxa… view at source ↗
Figure 3
Figure 3. Real-space ELF analyses of perfect and defective FCC surfaces. Two-dimensional ELF contour maps are displayed for the top layer (upper row) and 2nd layer (middle row) of Pt, Ag, and Rh. Within each metal panel, from left to right, the subplots illustrate the perfect lattice, the system with a top-layer vacancy, and the system with a 2nd-layer (subsurface) vacancy, respectively. The white dashed lines indicate the li… view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: Microscopic site configurations, vacancy kinetic barriers, and electrochemical free energy landscapes for HER and ORR on Pt(111). (a) Top view of the perfect FCC (111) surface illustrating the high-symmetry adsorption sites. Lattice atoms are color-coded by depth: surf…
Figure 5
Figure 5. Figure 5: Schematic illustration of the reversible structural transition on the Au(100) surface and its underlying energetic mechanism. Upper row: Reversible structural evolution between the pristine square lattice and the densified hexagonal reconstruction under thermal anneali…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

4 extracted references · 4 canonical work pages

  1. [1]

    Journal of Catalysis 328, 36–42 (2015) https://doi.org/10.1016/j.jcat.2014.12.033

    (1) Medford, A. J.; Vojvodic, A.; Hummelshøj, J. S.; Voss, J.; Abild -Pedersen, F.; Studt, F.; Bligaard, T.; Nilsson, A.; Nørskov, J. K. From the Sabatier Principle to a Predictive Theory of Transition-Metal Heterogeneous Catalysis. Journal of Catalysis 2015, 328, 36 –42. https://doi.org/10.1016/j.jcat.2014.12.033. (2) Cui, X.; Li, W.; Ryabchuk, P .; Jung...

    work page doi:10.1016/j.jcat.2014.12.033 2015

  2. [2]

    E(3)-equivariant graph neural networks for data-efficient and accurate interatomic potentials

    https://doi.org/10.1038/s41467-022-29939-5. (27) Musaelian, A.; Batzner, S.; Johansson, A.; Sun, L.; Owen, C. J.; Kornbluth, M.; Kozinsky, B. Learning Local Equivariant Representations for Large -Scale Atomistic Dynamics. Nat Commun 2023, 14 (1),

    work page doi:10.1038/s41467-022-29939-5 2023

  3. [3]

    Learning local equivariant representations for large-scale atomistic dynamics

    https://doi.org/10.1038/s41467-023-36329-y. (28) Hammer, R.; Meinel, K.; Krahn, O.; Widdra, W. Surface Reconstruction of Pt(001) Quantitatively Revisited. Phys. Rev. B 2016, 94 (19), 195406. https://doi.org/10.1103/PhysRevB.94.195406. (29) Shiihara, Y .; Kohyama, M.; Ishibashi, S. Origin of Surface Stress on Late Transition Metal Surfaces: Ab Initio Local...

    work page doi:10.1038/s41467-023-36329-y 2016

  4. [4]

    2025 , month = apr, number =

    https://doi.org/10.48550/ARXIV.2504.16068. (47) Wang, H.; Zhang, L.; Han, J.; E, W. DeePMD -Kit: A Deep Learning Package for Many -Body Potential Energy Representation and Molecular Dynamics. Computer Physics Communications 2018, 228, 178–184. https://doi.org/10.1016/j.cpc.2018.03.016. (48) Li, P .; Ding, F. Origin of the Herringbone Reconstruction of Au(...

    work page doi:10.48550/arxiv.2504.16068 2018

Pith tools

Reviewed June 29, 2026 · model on record in the stance chip above.