Scalable Prediction of Complex Surface Reconstructions under Operating Conditions via Harmony-Search-Based Global Optimization
Pith reviewed 2026-06-27 21:24 UTC · model grok-4.3
The pith
A harmony-search algorithm paired with machine learning potentials locates complex catalyst surface structures under reaction conditions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
HASGO integrates universal MLIPs with a harmony search algorithm, incorporating multi-head replay fine-tuning to overcome PES softening and stochastic structural perturbation for robust convergence. This enables identification of intricate surface oxide overlayers that align with atomic-resolution microscopy, resolving the square-pyramidal subsurface O5 motif on Ag(100) during ethylene epoxidation.
What carries the argument
Harmony-search-based Atomic Structural Global Optimization (HASGO), which maintains a harmony memory of candidate structures, applies stochastic perturbations, and evaluates energies via fine-tuned MLIPs to locate minimum-energy operando surface configurations.
If this is right
- Global optimization becomes feasible for expansive search spaces in low-symmetry operando systems
- Multi-head replay fine-tuning corrects PES inaccuracies without full retraining
- Stochastic perturbation provides fault-tolerant convergence to experimental structures
- The approach supplies atomic models that accelerate rational design of industrial catalysts
Where Pith is reading between the lines
- The same protocol could be tested on other metal surfaces and reactions to check whether similar oxide motifs appear
- Coupling HASGO outputs with in-situ spectroscopy might create closed-loop refinement of predicted structures
- Running the search at elevated temperature rather than zero Kelvin could reveal whether the O5 motif persists dynamically
Load-bearing premise
The fine-tuned machine learning interatomic potentials remain accurate on low-symmetry large-cell surface systems so that the located energy minima correspond to real physical structures.
What would settle it
Atomic-resolution microscopy images of the Ag(100) surface under ethylene epoxidation conditions that show a subsurface oxygen arrangement different from the predicted square-pyramidal O5 motif.
read the original abstract
The dynamic structural evolution of catalyst surfaces under operating conditions dictates catalytic performance, yet capturing these reconstructions atomically remains challenging. Global optimization based on machine learning interatomic potentials (MLIPs) is promising, but scaling to large-scale, low-symmetry operando systems is hindered by expansive search spaces and potential energy surface (PES) inaccuracies. Herein, we present Harmony-search-based Atomic Structural Global Optimization (HASGO), a framework integrating universal MLIPs with a harmony search algorithm. HASGO overcomes the problem of PES softening by incorporating a multi-head replay fine-tuning protocol. Moreover, the stochastic structural perturbation step in its algorithm offers a fault-tolerant strategy to enhance the robustness of global convergence. These enable HASGO to identify intricate surface oxide overlayers that align with atomic-resolution microscopy, thereby resolving the square-pyramidal subsurface O5 motif on Ag(100) during ethylene epoxidation. This scalable framework provides a robust approach for uncovering operando structures, accelerating the rational design of industrial catalysts.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces Harmony-search-based Atomic Structural Global Optimization (HASGO), a framework that integrates universal MLIPs with a harmony search algorithm and a multi-head replay fine-tuning protocol to scale global optimization to large, low-symmetry operando catalyst surfaces. It claims that this approach overcomes PES softening, enables robust convergence via stochastic perturbation, and identifies intricate surface oxide overlayers on Ag(100) during ethylene epoxidation that match atomic-resolution microscopy, thereby resolving the square-pyramidal subsurface O5 motif.
Significance. If the MLIP accuracy on the reported low-symmetry structures holds, the work provides a scalable computational route to predict operando reconstructions that are otherwise inaccessible, directly linking to experimental microscopy and supporting rational catalyst design.
major comments (1)
- [Results / Methods (protocol description)] The central claim that HASGO resolves the square-pyramidal subsurface O5 motif requires that the fine-tuned universal MLIP minima are physically accurate rather than artifacts. No section reports single-point DFT energies or forces on the final HASGO structures or on held-out large-cell, low-symmetry oxide configurations to quantify residual PES errors at the relevant scale and symmetry.
minor comments (2)
- [Methods] Clarify the exact definition and implementation of the multi-head replay fine-tuning protocol, including the replay buffer size, head architecture, and convergence criteria for the fine-tuning loss.
- [Results] The abstract states alignment with microscopy but the main text should include quantitative metrics (e.g., RMSD to experimental positions or simulated vs. observed image features) rather than qualitative statements.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The concern regarding explicit DFT validation of the final HASGO structures is well taken, and we address it directly below.
read point-by-point responses
-
Referee: [Results / Methods (protocol description)] The central claim that HASGO resolves the square-pyramidal subsurface O5 motif requires that the fine-tuned universal MLIP minima are physically accurate rather than artifacts. No section reports single-point DFT energies or forces on the final HASGO structures or on held-out large-cell, low-symmetry oxide configurations to quantify residual PES errors at the relevant scale and symmetry.
Authors: We agree that direct quantification of residual PES errors on the final large, low-symmetry structures is necessary to substantiate the claim. The multi-head replay fine-tuning was validated on smaller reference systems during development, and the identified O5 motif was cross-checked against atomic-resolution microscopy, but these do not replace single-point DFT benchmarks at the target scale. In the revised manuscript we will add single-point DFT energies and forces (using the same functional and settings as the original training data) for the converged HASGO O5 structure on Ag(100) as well as for a set of held-out large-cell oxide configurations. These results will be reported in a new subsection of the Results and discussed in the context of residual error. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper's core claim is that HASGO, combining universal MLIPs with harmony search and a multi-head replay fine-tuning protocol, enables discovery of operando surface reconstructions (including the O5 motif) that match microscopy data. No equations, fitting procedures, or self-citation chains are shown that reduce the reported structures or their identification to inputs by construction; the fine-tuning is presented as a general robustness measure rather than a fit to the target motifs, and experimental alignment functions as external corroboration. The derivation chain therefore remains self-contained and predictive rather than tautological.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Weaver, J. F. Surface chemistry of late transition metal oxides.Chemical Reviews113, 4164–4215 (2013)
2013
-
[2]
Polo-Garzon, F., Bao, Z., Zhang, X., Huang, W. & Wu, Z. Surface reconstructions of metal oxides and the consequences on catalytic chemistry.ACS Catalysis9, 5692–5707 (2019)
2019
-
[3]
Bossche, M. V. d. & Gronbeck, H. Methane oxidation over PdO (101) revealed by first-principles kinetic modeling.Journal of the American Chemical Society137, 12035–12044 (2015)
2015
-
[4]
& Iglesia, E
Chin, Y.-H., Garc´ ıa-Dieguez, M. & Iglesia, E. Dynamics and thermodynamics of Pd–PdO phase transi- tions: effects of Pd cluster size and kinetic implications for catalytic methane combustion.The Journal of Physical Chemistry C120, 1446–1460 (2016). 11
2016
-
[5]
& Jeong, H
Ryu, S.-H., Kim, S., Lee, H., Choi, J.-H. & Jeong, H. Robust palladium oxide nano-cluster cat- alysts using atomic ions and strong interactions for high-performance methane oxidation.Nature Communications15, 8348 (2024)
2024
-
[6]
Yang, W.et al.Active Palladium Structures on Ceria Obtained by Tuning Pd–Pd Distance for Efficient Methane Combustion.ACS Catalysis14, 16459–16468 (2024)
2024
-
[7]
Ammon, C.et al.Dissociation and oxidation of methanol on Cu(110).Surface Science507, 845–850 (2002)
2002
-
[8]
J.et al.Investigation of the evolution of Pd–Pt supported on ceria for dry and wet methane oxidation.Nature Communications13, 5080 (2022)
Divins, N. J.et al.Investigation of the evolution of Pd–Pt supported on ceria for dry and wet methane oxidation.Nature Communications13, 5080 (2022)
2022
-
[9]
Merrill, M. D. & Dougherty, R. C. Metal oxide catalysts for the evolution of O2 from H2O.The Journal of Physical Chemistry C112, 3655–3666 (2008)
2008
-
[10]
Hu, X.et al.Engineering nonprecious metal oxides electrocatalysts for two-electron water oxidation to H2O2.Advanced Energy Materials12, 2201466 (2022)
2022
-
[11]
Rønne, N.et al.Atomistic structure search using local surrogate model.The Journal of Chemical Physics157(2022)
2022
-
[12]
Dominguez, J. E. Q., Christiansen, M.-P. V., Neyman, K. M., Hammer, B. & Bruix, A. Efficient grand canonical global optimization with on-the-fly-trained machine-learning interatomic potentials.arXiv preprint arXiv:2509.19968(2025)
arXiv 2025
-
[13]
Shi, X.et al.Accessing complex reconstructed material structures with hybrid global optimization accelerated via on-the-fly machine learning.Chemical Science14, 8777–8784 (2023)
2023
-
[14]
& Liu, Z.-P
Chen, D., Shang, C. & Liu, Z.-P. Automated search for optimal surface phases (ASOPs) in grand canonical ensemble powered by machine learning.The Journal of Chemical Physics156(2022)
2022
-
[15]
Chen, D.et al.Square-pyramidal subsurface oxygen [Ag4OAg] drives selective ethene epoxidation on silver.Nature Catalysis7, 536–545 (2024)
2024
-
[16]
W., Oganov, A
Glass, C. W., Oganov, A. R. & Hansen, N. USPEX–Evolutionary crystal structure prediction.Computer Physics Communications175, 713–720 (2006)
2006
-
[17]
Wang, Y., Lv, J., Zhu, L. & Ma, Y. CALYPSO: A method for crystal structure prediction.Computer Physics Communications183, 2063–2070 (2012)
2063
-
[18]
& Dolg, M
Zhang, J. & Dolg, M. ABCluster: the artificial bee colony algorithm for cluster global optimization. Physical Chemistry Chemical Physics17, 24173–24181 (2015)
2015
-
[19]
Gruznev, D.et al.Atomic structure and electronic properties of the two-dimensional (Au, Al)/Si (111) 2×2 compound.Physical Review B92, 245407 (2015)
2015
-
[20]
V.et al.Two-dimensional In–Sb compound on silicon as a quantum spin Hall insulator
Gruznev, D. V.et al.Two-dimensional In–Sb compound on silicon as a quantum spin Hall insulator. Nano Letters18, 4338–4345 (2018)
2018
-
[21]
G., Kvashnin, D
Kvashnin, A. G., Kvashnin, D. G. & Oganov, A. R. Novel unexpected reconstructions of (100) and (111) surfaces of NaCl: Theoretical prediction.Scientific Reports9, 14267 (2019)
2019
-
[22]
& Gastegger, M
Sch¨ utt, K., Unke, O. & Gastegger, M. Equivariant message passing for the prediction of tensorial prop- erties and molecular spectra.Proceedings of the 38th International Conference on Machine Learning 9377–9388 (2021)
2021
-
[23]
P., Simm, G., Ortner, C
Batatia, I., Kovacs, D. P., Simm, G., Ortner, C. & Cs´ anyi, G. MACE: Higher order equivariant message passing neural networks for fast and accurate force fields.Advances in neural information processing systems35, 11423–11436 (2022). 12
2022
-
[24]
Neumann, M.et al.Orb: A fast, scalable neural network potential.arXiv preprint arXiv:2410.22570 (2024)
arXiv 2024
-
[25]
Liao, Y.-L. & Smidt, T. Equiformer: Equivariant graph attention transformer for 3d atomistic graphs. arXiv preprint arXiv:2206.11990(2022)
arXiv 2022
-
[26]
Xu, J., Xie, W., Han, Y. & Hu, P. Atomistic insights into the oxidation of flat and stepped platinum surfaces using large-scale machine learning potential-based grand-canonical monte carlo.ACS Catalysis 12, 14812–14824 (2022)
2022
-
[27]
Xu, J., Cao, X.-M. & Hu, P. Accelerating metadynamics-based free-energy calculations with adaptive machine learning potentials.Journal of Chemical Theory and Computation17, 4465–4476 (2021)
2021
-
[28]
& Ouyang, R
Zhang, Y., Li, Z., Han, Z.-K. & Ouyang, R. Global optimization of cation ordering in perovskites by recommendation-based basin-hopping.Journal of Chemical Theory and Computation20, 6971–6979 (2024)
2024
-
[29]
Journal of Chemical Theory and Computation20, 7370–7382 (2024)
Lee, C.et al.Galaxydock-dl: Protein–ligand docking by global optimization and neural network energy. Journal of Chemical Theory and Computation20, 7370–7382 (2024)
2024
-
[30]
Yang, C., Wu, C., Xie, W., Xie, D. & Hu, P. General reactive element-based machine learning potentials for heterogeneous catalysis.Nature Catalysis1–14 (2025)
2025
-
[31]
Liu, Y.et al.An automated framework for exploring and learning potential-energy surfaces.Nature Communications16, 7666 (2025)
2025
-
[32]
Deng, B.et al.Chgnet as a pretrained universal neural network potential for charge-informed atomistic modelling.Nature Machine Intelligence5, 1031–1041 (2023)
2023
-
[33]
Choudhary, K.et al.Unified graph neural network force-field for the periodic table: solid state applications.Digital Discovery2, 346–355 (2023)
2023
-
[34]
& Ong, S
Chen, C. & Ong, S. P. A universal graph deep learning interatomic potential for the periodic table. Nature Computational Science2, 718–728 (2022)
2022
-
[35]
Merchant, A.et al.Scaling deep learning for materials discovery.Nature624, 80–85 (2023)
2023
-
[36]
npj Computational Materials10, 94 (2024)
Zhang, D.et al.Pretraining of attention-based deep learning potential model for molecular simulation. npj Computational Materials10, 94 (2024)
2024
-
[37]
arXiv preprint arXiv:2405.04967(2024)
Yang, H.et al.Mattersim: A deep learning atomistic model across elements, temperatures and pressures. arXiv preprint arXiv:2405.04967(2024)
Pith/arXiv arXiv 2024
-
[38]
Gao, R.et al.A foundation machine learning potential with polarizable long-range interactions for materials modelling.Nature Communications16, 10484 (2025)
2025
-
[39]
Mazitov, A.et al.Pet-mad as a lightweight universal interatomic potential for advanced materials modeling.Nature Communications16, 10653 (2025)
2025
-
[40]
Deng, B.et al.Overcoming systematic softening in universal machine learning interatomic potentials by fine-tuning.arXiv preprint arXiv:2405.07105(2024)
arXiv 2024
-
[41]
Kulichenko, M.et al.Data generation for machine learning interatomic potentials and beyond.Chemical Reviews124, 13681–13714 (2024)
2024
-
[42]
Huo, Songze and Cao, Xiao-Ming. Accelerating High-Throughput Catalyst Screening by Direct Generation of Equilibrium Adsorption Structures.arXiv preprint arXiv:2512.15228(2025)
arXiv 2025
-
[43]
Kumar, A., Raghunathan, A., Jones, R., Ma, T. & Liang, P. Fine-tuning can distort pretrained features and underperform out-of-distribution.arXiv preprint arXiv:2202.10054(2022). 13
arXiv 2022
-
[44]
Ramasesh, V. V., Dyer, E. & Raghu, M. Anatomy of catastrophic forgetting: Hidden representations and task semantics.arXiv preprint arXiv:2007.07400(2020)
arXiv 2007
-
[45]
G., Allen, C
Radova, M., Stark, W. G., Allen, C. S., Maurer, R. J. & Bart´ ok, A. P. Fine-tuning foundation models of materials interatomic potentials with frozen transfer learning.npj Computational Materials11, 237 (2025)
2025
-
[46]
W., Kim, J
Geem, Z. W., Kim, J. H. & Loganathan, G. V. A new heuristic optimization algorithm: harmony search.Simulation76, 60–68 (2001)
2001
-
[47]
& Zeevi, Y
Eldar, Y., Lindenbaum, M., Porat, M. & Zeevi, Y. Y. The farthest point strategy for progressive image sampling.IEEE transactions on image processing6, 1305–1315 (1997)
1997
-
[48]
Batatia, Ilyes and Lin, Chen and Hart, Josephet al.Cross learning between electronic structure theories for unifying molecular, surface, and inorganic crystal foundation force fields.arXiv preprint arXiv:2510.25380(2025)
arXiv 2025
-
[49]
F., Zhang, F., Pan, L., Li, T
Weaver, J. F., Zhang, F., Pan, L., Li, T. & Asthagiri, A. Vacancy-mediated processes in the oxidation of CO on PdO (101).Accounts of Chemical Research48, 1515–1523 (2015)
2015
-
[50]
Zhang, F.et al.CO oxidation on PdO (101) during temperature-programmed reaction spectroscopy: role of oxygen vacancies.The Journal of Physical Chemistry C118, 28647–28661 (2014)
2014
-
[51]
A genetic algorithm tutorial.Statistics and computing4, 65–85 (1994)
Whitley, D. A genetic algorithm tutorial.Statistics and computing4, 65–85 (1994)
1994
-
[52]
& Eberhart, R
Kennedy, J. & Eberhart, R. Particle swarm optimization.Proceedings of ICNN’95 - International Conference on Neural Networks4, 1942–1948 (1995)
1942
-
[53]
& Dorigo, M
Socha, K. & Dorigo, M. Ant colony optimization for continuous domains.European journal of operational research185, 1155–1173 (2008)
2008
-
[54]
Mirjalili, S., Mirjalili, S. M. & Lewis, A. Grey wolf optimizer.Advances in Engineering Software69, 46–61 (2014)
2014
-
[55]
& Zakian, P
Kaveh, A. & Zakian, P. Improved GWO algorithm for optimal design of truss structures.Engineering with Computers34, 685–707 (2018)
2018
-
[56]
& Khiyabani, F
Nasiri, J. & Khiyabani, F. M. A whale optimization algorithm (WOA) approach for clustering.Cogent Mathematics & Statistics5, 1483565 (2018)
2018
-
[57]
Biogeography-based optimization.IEEE transactions on evolutionary computation12, 702–713 (2008)
Simon, D. Biogeography-based optimization.IEEE transactions on evolutionary computation12, 702–713 (2008)
2008
-
[58]
Mirjalili, S., Mirjalili, S. M. & Hatamlou, A. Multi-verse optimizer: a nature-inspired algorithm for global optimization.Neural Computing and Applications27, 495–513 (2016)
2016
-
[59]
& Goerigk, L
Grimme, S., Ehrlich, S. & Goerigk, L. Effect of the damping function in dispersion corrected density functional theory.Journal of Computational Chemistry32, 1456–1465 (2011)
2011
-
[60]
& Krieg, H
Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu.The Journal of Chemical Physics132(2010)
2010
-
[61]
Shipilin, M.et al.The influence of incommensurability on the long-range periodicity of the Pd (100)-(5×
-
[62]
R27◦-PdO (101).Surface Science660, 1–8 (2017)
2017
-
[63]
and Jia Chen, Z.et al.Unidirectional suppression of hydrogen oxidation on oxidized platinum clusters.Nature Communications4, 2500 (2013)
Hang Li, Y., Xing, J. and Jia Chen, Z.et al.Unidirectional suppression of hydrogen oxidation on oxidized platinum clusters.Nature Communications4, 2500 (2013)
2013
-
[64]
and Hu, P
Michaelides, A. and Hu, P. Catalytic water formation on platinum: A first-principles study.Journal of the American Chemical Society123, 4235–4242 (2001). 14
2001
-
[65]
and Gronbeck, H
Jørgensen, M. and Gronbeck, H. First-principles microkinetic modeling of methane oxidation over Pd (100) and Pd (111).ACS Catalysis6, 6730–6738 (2016)
2016
-
[66]
and Oh, J.et al.Sulfur atoms adsorbed on Cu (100) at low coverage: Characterization and stability against complexation.The Journal of Physical Chemistry B122, 963–971 (2018)
Walen, H., Liu, D. and Oh, J.et al.Sulfur atoms adsorbed on Cu (100) at low coverage: Characterization and stability against complexation.The Journal of Physical Chemistry B122, 963–971 (2018)
2018
-
[67]
& Altman, E
Zheng, G. & Altman, E. The oxidation mechanism of Pd (100).Surface Science504, 253–270 (2002)
2002
-
[68]
Kosteln´ ık, P.et al.The Pd (100)-(5×5) R27 ◦—O surface oxide: A leed, dft and stm study.Surface Science601, 1574–1581 (2007)
2007
-
[69]
T., Saidi, W
Li, M., Curnan, M. T., Saidi, W. A. & Yang, J. C. Uneven oxidation and surface reconstructions on stepped Cu (100) and Cu (110).Nano Letters22, 1075–1082 (2022)
2022
-
[70]
C., McFarlane, R
Zeng, H. C., McFarlane, R. A. & Mitchell, K. A. R. A LEED crystallographic investigation of some missing row models for the Cu(100)-(2 √ 2× √
-
[71]
R45 ◦-O surface structure.Surface Science Letters 208, L7–L14 (1989)
1989
-
[72]
W.et al.Theory of the oxygen-induced restructuring of Cu (110) and Cu (100) surfaces
Jacobsen, K. W.et al.Theory of the oxygen-induced restructuring of Cu (110) and Cu (100) surfaces. Physical Review Letters65, 1788 (1990)
1990
-
[73]
& McGaughey, A
Lee, M. & McGaughey, A. J. Energetics and kinetics of the c (2×2) to (2 √ 2× √
-
[74]
R45 ◦ transition during the early stages of Cu (100) oxidation.Physical Review Bondensed Matter and Materials Physics 83, 165447 (2011)
2011
-
[75]
& Henkelman, G
Lian, X., Xiao, P., Liu, R. & Henkelman, G. Calculations of oxygen adsorption-induced surface reconstruction and oxide formation on Cu (100).Chemistry of Materials29, 1472–1484 (2017)
2017
-
[76]
Monig, H.et al.Understanding scanning tunneling microscopy contrast mechanisms on metal oxides: a case study.ACS Nano7, 10233–10244 (2013)
2013
-
[77]
Calculations of oxygen adsorption-induced surface reconstruction and oxide formation on Cu(100).Chemistry of Materials2017,29(4), 1472–1484
Lian, X.; Xiao, P.; Liu, R.; Henkelman, G. Calculations of oxygen adsorption-induced surface reconstruction and oxide formation on Cu(100).Chemistry of Materials2017,29(4), 1472–1484
-
[78]
& Dosch, H
Stierle, A., Costina, I., Kumaragurubaran, S. & Dosch, H. In situ x-ray diffraction study ofag(100) at ambient oxygen pressures.The Journal of Physical Chemistry C111, 10998–11002 (2007)
2007
-
[79]
Rocca, M.et al.Phase transition of dissociatively adsorbed oxygen onag(001).Physical Review B61, 213 (2000)
2000
-
[80]
& Lin, X
Hua, M., Li, S., Tian, X., Li, G. & Lin, X. Independent switchable atomic silver quantum transistor via potential-driven surface reconstruction.The Journal of Chemical Physics163, 114101 (2025)
2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.