pith. sign in

arxiv: 2605.26588 · v1 · pith:NUPT7L33new · submitted 2026-05-26 · ⚛️ physics.comp-ph

Local Surrogates for Harmonic Vibrational Entropy in Multilattices

Pith reviewed 2026-06-29 15:04 UTC · model grok-4.3

classification ⚛️ physics.comp-ph
keywords vibrational entropylocal surrogatesmultilatticesharmonic approximationdefect thermodynamicscutoff error estimatessublattice resolution
0
0 comments X

The pith

Finite-range models let harmonic entropy be calculated from local sublattice surrogates with bounded error.

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

Harmonic vibrational entropy calculations normally require diagonalizing a large Hessian matrix, which costs time that grows with the cube of the system size. The paper proves that when atom interactions have finite range, the entropy contributions are local to each sublattice site, allowing replacement by a regression model trained on small clusters. This local approach scales linearly and maintains accuracy for repeated evaluations in defect studies of multilattices like alloys and semiconductors.

Core claim

For finite-range or screened atomistic models, sublattice-resolved locality and cutoff-error estimates justify replacing the global entropy calculation by a local, symmetry-respecting regression problem with controlled truncation error.

What carries the argument

Sublattice-resolved locality proofs and local regression surrogates that respect species and internal shift symmetries.

If this is right

  • Harmonic entropy can be evaluated at linear cost in system size for fixed cutoff radius.
  • Surrogates transfer across different supercell sizes in multispecies crystals.
  • Explicit controls on truncation and regression error are available for stability checks.
  • Method applies to semiconductors, ordered alloys, and crystals with multi-atom bases.

Where Pith is reading between the lines

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

  • Such surrogates could enable high-throughput screening of defect formation energies including vibrational contributions.
  • Extension to anharmonic effects might require similar locality arguments for higher-order force constants.

Load-bearing premise

The atomistic models have finite interaction range or are screened.

What would settle it

Observation of significant entropy contributions from atoms beyond the cutoff distance in a model with strictly finite-range interactions.

Figures

Figures reproduced from arXiv: 2605.26588 by Jiale Linghu, Tina Torabi, Yangshuai Wang.

Figure 1
Figure 1. Figure 1: Representative two-sublattice crystal structures ( [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Envelope of |𝜕𝑆𝑁 𝓁0 ,𝛼 ∕𝜕𝑢𝛽 (𝑛)| as a function of |𝑛 − 𝓁0 | in 𝑑 = 1, 2, 3. The 1D and 2D panels show the median over independent random base configurations, with the 10%–90% interquantile band shaded; the 3D panel shows a single finite-difference run due to cost. Dashed reference slopes correspond to (1 + 𝑟) −2𝑑 in each panel. 100 101 Cutoff radius rcut (cells) 10−5 10−4 10−3 10−2 Truncation error |S N 0,… view at source ↗
Figure 3
Figure 3. Figure 3: Truncation error |𝑆 𝑁 0,0 (𝑢) − 𝑆̃ 0,0 (𝑢; 𝑟cut)| as a function of cutoff radius 𝑟cut in 𝑑 = 1, 2, 3. Dashed reference slopes 𝑟 −𝑑 cut. The empirical decay sits below the upper bound; in 1D and 2D the extra decay reflects the localized displacement tails, while the 3D modulation is closer to the worst-case envelope. 5.2 Sublattice- and Species-Resolved Entropy Regression The next question is whether the si… view at source ↗
Figure 4
Figure 4. Figure 4: Multispecies ACE site model on Si SW. Held-out parity at body order [PITH_FULL_IMAGE:figures/full_fig_p022_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Species-resolved local-polynomial test on CdTe SW. (a) Held-out site parity for a body-order- [PITH_FULL_IMAGE:figures/full_fig_p023_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: (a) Wall-time scaling of direct diagonalization versus surrogate evaluation on the 1D diatomic chain [PITH_FULL_IMAGE:figures/full_fig_p024_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Si SW: cross-supercell transfer and 3D timing. (a) Held-out site-entropy parity for the ACE site [PITH_FULL_IMAGE:figures/full_fig_p025_7.png] view at source ↗
read the original abstract

Harmonic vibrational entropy is a key finite-temperature contribution to defect thermodynamics, but direct evaluation by dense Hessian diagonalization scales cubically with atom count and is too costly for supercell convergence, migration-path sampling, and high-throughput defect studies. We develop local surrogate models for harmonic entropy in multilattices, including semiconductors, ordered alloys, and multispecies crystals with multi-atom bases and internal degrees of freedom. Unlike Bravais lattices, multilattices contain internal-shift degrees of freedom and optical phonon modes coupled to acoustic strain; entropy models must therefore resolve sublattice and species labels. For finite-range or screened atomistic models, we prove sublattice-resolved locality and cutoff-error estimates that justify replacing the global entropy calculation by a local, symmetry-respecting regression problem with controlled truncation error. This turns vibrational entropy from a global spectral calculation into a reusable local site model with linear evaluation cost at fixed cutoff. Numerical tests confirm the predicted locality behavior and show that sublattice/species-resolved surrogates achieve accurate regression, transfer across supercell sizes, and linear-scaling evaluation on Stillinger--Weber Si and CdTe benchmarks. The resulting method enables repeated harmonic-entropy evaluations in multispecies defect calculations while retaining explicit stability, truncation, and surrogate-error controls.

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

2 major / 2 minor

Summary. The manuscript develops local surrogate models for harmonic vibrational entropy in multilattices (semiconductors, ordered alloys, multispecies crystals). For finite-range or screened atomistic models it proves sublattice-resolved locality together with cutoff-error estimates; these justify replacing global Hessian diagonalization by a local, symmetry-respecting regression problem whose evaluation cost is linear at fixed cutoff. Numerical benchmarks on Stillinger-Weber Si and CdTe are reported to confirm the predicted locality scaling, regression accuracy, and transfer across supercell sizes.

Significance. If the conditional proofs and error bounds hold, the work supplies an explicit, reusable site model with controlled truncation and surrogate error for a quantity that is otherwise cubic-cost. This directly enables repeated harmonic-entropy evaluations inside defect thermodynamics, migration-path sampling, and high-throughput studies of multilattices. Credit is due for (i) stating the finite-range/screened assumption up front, (ii) deriving sublattice/species-resolved locality rather than assuming it, and (iii) validating on two independent short-range benchmarks rather than a single toy system.

major comments (2)
  1. [Theoretical development / locality proof] The central claim rests on the proof of sublattice-resolved locality and the cutoff-error estimates. The abstract states that these hold precisely when the underlying model is finite-range or screened; the manuscript must therefore contain an explicit derivation (presumably in the theoretical section) showing how the finite-range assumption produces the stated truncation bound. Without that derivation the error-control guarantee cannot be verified.
  2. [Numerical results / benchmarks] Table or figure reporting regression accuracy (R², MAE, or similar) on the Stillinger-Weber Si and CdTe benchmarks should be accompanied by the exact regression procedure, feature set, and how the predicted scaling with cutoff is quantified. The current high-level description leaves open whether the numerical confirmation actually matches the analytic cutoff-error estimate.
minor comments (2)
  1. Notation for sublattice and species labels in the surrogate model equations should be introduced once and used consistently; a small table summarizing the symbols would improve readability.
  2. The abstract mentions 'explicit stability, truncation, and surrogate-error controls'; a short dedicated paragraph or subsection listing the three controls and where each appears in the manuscript would help readers locate them.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of the work's significance and for the constructive major comments. We address each point below and will revise the manuscript to strengthen the presentation of the theoretical derivation and the numerical procedures.

read point-by-point responses
  1. Referee: [Theoretical development / locality proof] The central claim rests on the proof of sublattice-resolved locality and the cutoff-error estimates. The abstract states that these hold precisely when the underlying model is finite-range or screened; the manuscript must therefore contain an explicit derivation (presumably in the theoretical section) showing how the finite-range assumption produces the stated truncation bound. Without that derivation the error-control guarantee cannot be verified.

    Authors: We agree that an explicit, self-contained derivation is required to allow verification of the error-control guarantee. The current manuscript states the finite-range/screened assumption and asserts the resulting locality and truncation bounds, but the step-by-step derivation from the assumption to the sublattice-resolved bounds is not presented with sufficient detail. In the revised version we will add a dedicated subsection (or appendix) that derives the locality result and the cutoff-error estimate directly from the finite-range hypothesis, including the relevant matrix-element decay arguments and the symmetry-respecting truncation. revision: yes

  2. Referee: [Numerical results / benchmarks] Table or figure reporting regression accuracy (R², MAE, or similar) on the Stillinger-Weber Si and CdTe benchmarks should be accompanied by the exact regression procedure, feature set, and how the predicted scaling with cutoff is quantified. The current high-level description leaves open whether the numerical confirmation actually matches the analytic cutoff-error estimate.

    Authors: We will expand the numerical section to include (i) the precise regression algorithm and hyperparameters, (ii) the explicit feature construction (local atomic neighborhoods with species and sublattice labels up to the chosen cutoff), and (iii) a direct quantitative comparison between the observed regression errors and the analytic cutoff-error bounds derived in the theory section. A new table or supplementary figure will report R², MAE, and scaling of the error with cutoff radius for both benchmarks, confirming consistency with the theoretical predictions. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is conditional on stated model assumptions

full rationale

The central claim is a conditional proof of sublattice-resolved locality and cutoff-error estimates that hold precisely when the atomistic model is finite-range or screened; this assumption is stated explicitly and required for the result. No load-bearing self-citation, self-definitional step, or fitted input renamed as prediction appears. Numerical benchmarks on independent short-range models (Stillinger-Weber Si, CdTe) provide external validation of the predicted scaling. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the domain assumption of finite-range or screened interactions; no free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption Finite-range or screened atomistic models exhibit sublattice-resolved locality for vibrational entropy contributions.
    Invoked to justify replacing global calculation with local surrogates.

pith-pipeline@v0.9.1-grok · 5754 in / 1035 out tokens · 12949 ms · 2026-06-29T15:04:28.456000+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

57 extracted references · 44 canonical work pages · 1 internal anchor

  1. [1]

    Atomic cluster expansion: Completeness, efficiency and stability

    Bachmayr, M., Csanyi, G., Dusson, G., Drautz, R., Etter, S., van der Oord, C., Ortner, C., 2022. Atomic cluster expansion: Completeness, efficiency and stability. Journal of Computational Physics 454, 110946. doi:10.1016/j.jcp.2022.110946

  2. [2]

    Polynomial approximation of symmetric functions

    Bachmayr, M., Dusson, G., Ortner, C., Drautz, R., 2024. Polynomial approximation of symmetric functions. Mathematics of Computation 93, 811–839. doi:10.1090/mcom/3868

  3. [3]

    Boundsforthetraceoftheinverseandthedeterminantofsymmetricpositive definite matrices

    Bai,Z.,Golub,G.H.,1996. Boundsforthetraceoftheinverseandthedeterminantofsymmetricpositive definite matrices. Annals of Numerical Mathematics 4, 29–38

  4. [4]

    On representing chemical environments

    Bartók, A.P., Kondor, R., Csányi, G., 2013. On representing chemical environments. Physical Review B 87, 184115. doi:10.1103/PhysRevB.87.184115

  5. [5]

    Gaussian approximation potentials: The accuracy of quantum mechanics, without the electrons

    Bartók, A.P., Payne, M.C., Kondor, R., Csányi, G., 2010. Gaussian approximation potentials: The accuracy of quantum mechanics, without the electrons. Physical Review Letters 104, 136403. doi:10. 1103/PhysRevLett.104.136403

  6. [6]

    Higher-order far-field boundary conditions for crys- talline defects

    Braun, J., Ortner, C., Wang, Y., Zhang, L., 2025. Higher-order far-field boundary conditions for crys- talline defects. SIAM Journal on Numerical Analysis 63, 520–541. doi:10.1137/24M165836X

  7. [7]

    Computer modelling of defects in solids

    Catlow, C.R.A., 1986. Computer modelling of defects in solids. Annual Review of Materials Science 16, 517–548. doi:10.1146/annurev.ms.16.080186.002505

  8. [8]

    QM/MMmethodsforcrystallinedefects.part3: Machine-learned MM models

    Chen,H.,Ortner,C.,Wang,Y.,2022. QM/MMmethodsforcrystallinedefects.part3: Machine-learned MM models. Multiscale Modeling & Simulation 20, 1490–1518. doi:10.1137/21M1441122

  9. [9]

    On the mathematical foundations of learning

    Cucker, F., Smale, S., 2002. On the mathematical foundations of learning. Bulletin of the American Mathematical Society 39, 1–49. doi:10.1090/S0273-0979-01-00923-5

  10. [10]

    Embedded-atom method: Derivation and application to impurities, surfaces, andotherdefectsinmetals

    Daw, M.S., Baskes, M.I., 1984. Embedded-atom method: Derivation and application to impurities, surfaces, andotherdefectsinmetals. PhysicalReviewB29, 6443–6453. doi:10.1103/PhysRevB.29. 6443

  11. [11]

    Atomicclusterexpansionforaccurateandtransferableinteratomicpotentials

    Drautz,R.,2019. Atomicclusterexpansionforaccurateandtransferableinteratomicpotentials. Physical Review B 99, 014104. doi:10.1103/PhysRevB.99.014104

  12. [12]

    Uniform accuracy of the quasicontinuum method

    E, W., Lu, J., Yang, J.Z., 2006. Uniform accuracy of the quasicontinuum method. Physical Review B 74, 214115. doi:10.1103/PhysRevB.74.214115

  13. [13]

    Anharmonic free energies and phonon dispersions from the stochastic self-consistent harmonic approximation

    Errea, I., Calandra, M., Mauri, F., 2014. Anharmonic free energies and phonon dispersions from the stochastic self-consistent harmonic approximation. Physical Review B 89, 064302. doi:10.1103/ PhysRevB.89.064302

  14. [14]

    Method to extract anharmonic force constants from first principles calculations

    Esfarjani, K., Stokes, H.T., 2008. Method to extract anharmonic force constants from first principles calculations. Physical Review B 77, 144112. doi:10.1103/PhysRevB.77.144112

  15. [15]

    Evaluation of harmonic methods for calculating the free energy of defects in solids

    Foiles, S.M., 1994. Evaluation of harmonic methods for calculating the free energy of defects in solids. Physical Review B 49, 14930–14938. doi:10.1103/PhysRevB.49.14930

  16. [16]

    New Monte Carlo method to compute the free energy of arbitrary solids

    Frenkel, D., Ladd, A.J.C., 1984. New Monte Carlo method to compute the free energy of arbitrary solids. application to the fcc and hcp phases of hard spheres. Journal of Chemical Physics 81, 3188–

  17. [17]

    doi:10.1063/1.448024

  18. [18]

    Calculation of the vacancy formation energy in aluminium

    Gillan, M.J., 1989. Calculation of the vacancy formation energy in aluminium. Journal of Physics: Condensed Matter 1, 689–711. doi:10.1088/0953-8984/1/4/005. 31

  19. [19]

    Breakdown of the Arrhenius law in de- scribing vacancy formation energies: The importance of local anharmonicity revealed byab initiother- modynamics

    Glensk, A., Grabowski, B., Hickel, T., Neugebauer, J., 2014. Breakdown of the Arrhenius law in de- scribing vacancy formation energies: The importance of local anharmonicity revealed byab initiother- modynamics. Physical Review X 4, 011018. doi:10.1103/PhysRevX.4.011018

  20. [20]

    Reinforcing materials modelling by encoding the structures of defects in crystalline solids into distortion scores

    Goryaeva, A.M., Maréchal, C., Marinica, M.C., 2020. Reinforcing materials modelling by encoding the structures of defects in crystalline solids into distortion scores. Nature Communications 11, 4691. doi:10.1038/s41467-020-18282-2

  21. [21]

    Physical Review B 79, 134106

    Grabowski, B., Ismer, L., Hickel, T., Neugebauer, J., 2009.Ab initioup to the melting point: Anhar- monicity and vacancies in aluminum. Physical Review B 79, 134106. doi:10.1103/PhysRevB.79. 134106

  22. [22]

    Calculation of dislocation binding to helium-vacancy defects in tungsten using hybrid ab initio-machine learning methods

    Grigorev, P., Goryaeva, A.M., Marinica, M.C., Kermode, J.R., Swinburne, T.D., 2023. Calculation of dislocation binding to helium-vacancy defects in tungsten using hybrid ab initio-machine learning methods. Acta Materialia 247, 118734. doi:10.1016/j.actamat.2023.118734

  23. [23]

    Approximating spectral sums of large-scale matrices usingstochasticChebyshevapproximations

    Han, I., Malioutov, D., Avron, H., Shin, J., 2017. Approximating spectral sums of large-scale matrices usingstochasticChebyshevapproximations. SIAMJournalonScientificComputing39,A1558–A1585. doi:10.1137/16M1078148

  24. [24]

    Approximation properties of ACE basis func- tions

    Hashemi, M., Wittwer, S., Suldér, B., Drautz, R., 2023. Approximation properties of ACE basis func- tions. Multiscale Modeling & Simulation 21, 978–1010. doi:10.1137/22M1535040

  25. [25]

    Temperaturedependenteffectivepotential method for accurate free energy calculations of solids

    Hellman,O.,Steneteg,P.,Abrikosov,I.A.,Simak,S.I.,2013. Temperaturedependenteffectivepotential method for accurate free energy calculations of solids. Physical Review B 87, 104111. doi:10.1103/ PhysRevB.87.104111

  26. [26]

    Improved tangent estimate in the nudged elastic band method for findingminimumenergypathsandsaddlepoints

    Henkelman, G., Jónsson, H., 2000. Improved tangent estimate in the nudged elastic band method for findingminimumenergypathsandsaddlepoints. JournalofChemicalPhysics113,9978–9985. doi:10. 1063/1.1323224

  27. [27]

    doi:10.1137/1.9780898717778

    Higham,N.J.,2008.FunctionsofMatrices: TheoryandComputation.SocietyforIndustrialandApplied Mathematics, Philadelphia. doi:10.1137/1.9780898717778

  28. [28]

    In operando active learning of interatomic interaction during large- scale simulations

    Hodapp, M., Shapeev, A., 2020. In operando active learning of interatomic interaction during large- scale simulations. Machine Learning: Science and Technology 1, 045005. doi:10.1088/2632-2153/ aba373

  29. [29]

    Machine-learning potentials enable predictive and tractable high- throughput screening of random alloys

    Hodapp, M., Shapeev, A., 2021. Machine-learning potentials enable predictive and tractable high- throughput screening of random alloys. Physical Review Materials 5, 113802. doi:10.1103/ PhysRevMaterials.5.113802

  30. [30]

    Analysisofstablescrewdislocationconfigurationsinanantiplanelattice model

    Hudson,T.,Ortner,C.,2014a. Analysisofstablescrewdislocationconfigurationsinanantiplanelattice model. SIAM Journal on Mathematical Analysis 46, 291–320. doi:10.1137/130941318

  31. [31]

    Existence and stability of a screw dislocation under anti-plane deforma- tion

    Hudson, T., Ortner, C., 2014b. Existence and stability of a screw dislocation under anti-plane deforma- tion. ArchiveforRationalMechanicsandAnalysis213,887–929. doi:10.1007/s00205-014-0746-9

  32. [32]

    A stochastic estimator of the trace of the influence matrix for Lapla- cian smoothing splines

    Hutchinson, M.F., 1990. A stochastic estimator of the trace of the influence matrix for Lapla- cian smoothing splines. Communications in Statistics — Simulation and Computation 19, 433–450. doi:10.1080/03610919008812866

  33. [33]

    Statisticalmechanicsoffluidmixtures

    Kirkwood,J.G.,1935. Statisticalmechanicsoffluidmixtures. JournalofChemicalPhysics3,300–313. doi:10.1063/1.1749657

  34. [34]

    Crack propagation in b.c.c

    Kohlhoff, S., Gumbsch, P., Fischmeister, H.F., 1991. Crack propagation in b.c.c. crystals studied with a combined finite-element and atomistic model. Philosophical Magazine A 64, 851–878. doi:10.1080/ 01418619108213953. 32

  35. [35]

    Ma- chine learning surrogate models for strain-dependent vibrational properties and migration rates of point defects

    Lapointe, C., Swinburne, T.D., Proville, L., Becquart, C.S., Mousseau, N., Marinica, M.C., 2022. Ma- chine learning surrogate models for strain-dependent vibrational properties and migration rates of point defects. Physical Review Materials 6, 113803. doi:10.1103/PhysRevMaterials.6.113803

  36. [36]

    Machine learning surrogate models for prediction of point defect vibrational entropy

    Lapointe,C.,Swinburne,T.D., Thiry,L.,Mallat,S.,Proville,L.,Becquart, C.S.,Marinica,M.C.,2020. Machine learning surrogate models for prediction of point defect vibrational entropy. Physical Review Materials 4, 063802. doi:10.1103/PhysRevMaterials.4.063802

  37. [37]

    Theatomicsimulationenvironment—aPythonlibraryfor working with atoms

    Larsen,A.H.,Mortensen,J.J.r.,Blomqvist,J.,Castelli,I.E.,Christensen,R.,Dułak,M.,Friis,J.,Groves, M.N.,Hammer,B.r.,Hargus,C.,etal.,2017. Theatomicsimulationenvironment—aPythonlibraryfor working with atoms. Journal of Physics: Condensed Matter 29, 273002. doi:10.1088/1361-648X/ aa680e

  38. [38]

    Fine-tuning universal machine-learned interatomicpotentials: Atutorialonmethodsandapplications

    Liu, X., Zeng, K., Luo, Z., Wang, Y., Zhao, T., Xu, Z., 2026. Fine-tuning universal machine-learned interatomicpotentials: Atutorialonmethodsandapplications. JournalofAppliedPhysics139,041101. doi:10.1063/5.0299305

  39. [39]

    Atomistic-to-continuum coupling

    Luskin, M., Ortner, C., 2013. Atomistic-to-continuum coupling. Acta Numerica 22, 397–508. doi:10. 1017/S0962492913000068

  40. [40]

    Performant implementation of the atomic cluster expansion (PACE) and application to copper and silicon

    Lysogorskiy,Y.,vanderOord,C.,Bochkarev,A.,Menon,S.,Rinaldi,M.,Hammerschmidt,T.,Mrovec, M., Thompson, A., Csányi, G., Ortner, C., Drautz, R., 2021. Performant implementation of the atomic cluster expansion (PACE) and application to copper and silicon. npj Computational Materials 7, 97. doi:10.1038/s41524-021-00559-9

  41. [41]

    Theory of Lattice Dynamics in the Harmonic Approximation

    Maradudin, A.A., Montroll, E.W., Weiss, G.H., Ipatova, I.P., 1971. Theory of Lattice Dynamics in the Harmonic Approximation. 2nd ed., Academic Press, New York

  42. [42]

    The role of the activation barriers and the diffusion path on the rate constants for vacancy diffusion in bcc Fe

    Marinica, M.C., Becquart, C.S., Domain, C., 2013. The role of the activation barriers and the diffusion path on the rate constants for vacancy diffusion in bcc Fe. Journal of Nuclear Materials 433, 491–497. doi:10.1016/j.jnucmat.2012.04.024

  43. [43]

    Elastic far-field decay from dislocations in multilat- tices

    Olson, E., Ortner, C., Wang, Y., Zhang, L., 2023. Elastic far-field decay from dislocations in multilat- tices. Multiscale Modeling & Simulation 21, 1379–1409. doi:10.1137/22M1502021

  44. [44]

    Moment tensor potentials: A class of systematically improvable interatomic po- tentials

    Shapeev, A.V., 2016. Moment tensor potentials: A class of systematically improvable interatomic po- tentials. Multiscale Modeling & Simulation 14, 1153–1173. doi:10.1137/15M1054183

  45. [45]

    Computer simulation of local order in condensed phases of silicon

    Stillinger, F.H., Weber, T.A., 1985. Computer simulation of local order in condensed phases of silicon. Physical Review B 31, 5262–5271. doi:10.1103/PhysRevB.31.5262

  46. [46]

    Unsupervised calculation of free energy barriers in large crys- talline systems

    Swinburne, T.D., Marinica, M.C., 2018. Unsupervised calculation of free energy barriers in large crys- talline systems. Physical Review Letters 120, 135503. doi:10.1103/PhysRevLett.120.135503

  47. [47]

    Modeling solid-state chemistry: Interatomic potentials for multicomponent systems

    Tersoff, J., 1989. Modeling solid-state chemistry: Interatomic potentials for multicomponent systems. Physical Review B 39, 5566–5568. doi:10.1103/PhysRevB.39.5566

  48. [48]

    LAMMPS – a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales

    Thompson, A.P., Aktulga, H.M., Berger, R., Bolintineanu, D.S., Brown, W.M., Crozier, P.S., in ’t Veld, P.J., Kohlmeyer, A., Moore, S.G., Nguyen, T.D., et al., 2022. LAMMPS – a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Computer Physics Communications 271, 108171. doi:10.1016/j.cpc.2021.108171

  49. [49]

    Surrogate models for vibrational entropy based on a spatial decomposition

    Torabi, T., Wang, Y., Ortner, C., 2025. Surrogate models for vibrational entropy based on a spatial decomposition. Multiscale Modeling & Simulation 23, 514–544. doi:10.1137/24M165168X

  50. [50]

    Latticegreenfunctionforextendeddefectcalculations: Computationanderrores- timationwithlong-rangeforces

    Trinkle, D.R., 2008. Latticegreenfunctionforextendeddefectcalculations: Computationanderrores- timationwithlong-rangeforces. PhysicalReviewB78, 014110. doi:10.1103/PhysRevB.78.014110. 33

  51. [51]

    Fastestimationoftr(𝑓(𝑎))viastochasticLanczosquadrature

    Ubaru,S.,Chen,J.,Saad,Y.,2017. Fastestimationoftr(𝑓(𝑎))viastochasticLanczosquadrature. SIAM Journal on Matrix Analysis and Applications 38, 1075–1099. doi:10.1137/16M1104974

  52. [52]

    Frequency factors and isotope effects in solid state rate processes

    Vineyard, G.H., 1957. Frequency factors and isotope effects in solid state rate processes. Journal of Physics and Chemistry of Solids 3, 121–127. doi:10.1016/0022-3697(57)90059-8

  53. [53]

    Thermodynamics of Crystals

    Wallace, D.C., 1972. Thermodynamics of Crystals. Wiley, New York

  54. [54]

    A posteriori error estimate and adaptivity for QM/MMmodelsofcrystallinedefects

    Wang, Y., Kermode, J.R., Ortner, C., Zhang, L., 2024. A posteriori error estimate and adaptivity for QM/MMmodelsofcrystallinedefects. ComputerMethodsinAppliedMechanicsandEngineering428, 117097. doi:10.1016/j.cma.2024.117097

  55. [55]

    Scattered Data Approximation

    Wendland, H., 2004. Scattered Data Approximation. Cambridge University Press, Cambridge

  56. [56]

    ACEpotentials.jl: A Julia implementationoftheatomicclusterexpansion

    Witt, W.C., van der Oord, C., Gelzinyte, E., Jarvinen, T., Ross, A., Darby, J.P., Ho, C.H., Baldwin, W.J., Sachs, M., Kermode, J., Bernstein, N., Csanyi, G., Ortner, C., 2023. ACEpotentials.jl: A Julia implementationoftheatomicclusterexpansion. TheJournalofChemicalPhysics159,164101. doi:10. 1063/5.0158783

  57. [57]

    Directcalculationofthelatticegreenfunctionwitharbitraryinteractions for general crystals

    Yasi,J.A.,Trinkle,D.R.,2012. Directcalculationofthelatticegreenfunctionwitharbitraryinteractions for general crystals. Physical Review E 85, 066706. doi:10.1103/PhysRevE.85.066706. 34