REVIEW 3 major objections 6 minor 33 references
Accurate Prediction of the $\alpha \to \beta$ Phase Transformation Temperature in Tin via Full Anharmonic Treatment
T0 review · 3 major / 6 minor · reviewed 2026-08-01 · deepseek-v4-flash
Pith's one-line read Anharmonicity, not the functional, sets tin's α→β transition near 288 K.
desk verdict Solid paper showing beta-Sn is strongly anharmonic, but the headline number mixes classical anharmonicity with quantum effects—qualitatively right, quantitatively overstated. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is a machine-learned Atomic Cluster Expansion (ACE) interatomic potential fitted to PBE density-functional theory data, which reproduces the potential energy surface of both α- and β-Sn and enables large-scale molecular dynamics. Its role is to make thermodynamic integration feasible while keeping the electronic reference fixed. The decisive comparison is between two free energies computed on this same potential: the quasi-harmonic free energy, built from volume-dependent phonon frequencies via Eq. (1), and the thermodynamic-integration free energy, which couples the ACE potential to an Einstein crystal and integrates over the coupling parameter λ (Eqs. 2–3). Their differe
What would settle it
Compute the same QHA-vs-TI comparison using a potential trained on a functional with a smaller static energy difference, such as PBE+U with U≈1–1.5 eV; if the TI transition temperature moves away from 286 K by more than about 10 K, the agreement with experiment partly reflects a cancellation between the static-energy error and the anharmonic correction. More directly, a quantitative calorimetric measurement of β-Sn's heat capacity above 200 K, compared with the TI prediction including the electronic and dilation terms, would falsify the magnitude of the anharmonic free energy if the mismatch e
Extended reading notes
Core claim
The central discovery is that β-Sn is strongly anharmonic in a temperature-dependent way that stabilizes it relative to α-Sn. The explicit anharmonic free energy, defined as the difference between the thermodynamic-integration and quasi-harmonic free energies (ΔF_anh = F_TI − F_QHA), reaches roughly −11 meV/atom at 300 K for β-Sn and is negligible for α-Sn. This anharmonic stabilization lowers the predicted α→β transformation temperature from 377 K (QHA) to 288 K (TI), almost exactly the experimental 286 K. The claim is corroborated by four independent observations: excess heat capacity above the Dulong–Petit limit, a temperature-driven collapse of a 3–3.5 THz feature in the vibrational dens
Load-bearing premise
The interpretation of ΔF_anh = F_TI − F_QHA as purely 'explicit anharmonic free energy' assumes the two free energies are directly comparable, yet TI is run with classical nuclei while QHA includes zero-point quantum motion, so nuclear-quantum effects are not separately accounted for and may contribute to the 85 K shift.
Editorial extensions
If this is right
- Predictions of phase boundaries in soft polymorphic metals should include explicit anharmonicity; quasi-harmonic treatments can be off by tens to hundreds of kelvin.
- The observed stability of white tin at ambient conditions is partly a temperature-driven dynamical effect, not just a static-energy preference.
- With anharmonicity accounted for, the residual uncertainty in the transition temperature is dominated by the 0 K static energy difference, so improving the exchange-correlation description (e.g., beyond PBE) becomes the next bottleneck.
- The combination of machine-learned potentials and thermodynamic integration offers a practical route to fully anharmonic free energies at near-ab initio accuracy for materials with shallow potential energy landscapes.
- The TI–QHA difference provides a quantitative decomposition of heat capacity and free energy into harmonic, dilation, and explicit phonon–phonon contributions, enabling direct comparison with calorimetric data.
Reading between the lines
- The near-exact agreement with 286 K may partly result from a cancellation between PBE's overestimated static energy difference and the anharmonic correction, as the paper itself concedes; a test with a functional giving a smaller static ΔE (e.g., PBE+U) would reveal whether the anharmonic mechanism is robust or the match is fortuitous.
- The finding suggests that other soft metallic phases with shallow potential wells may harbor large hidden anharmonic contributions, implying that quasi-harmonic screening of phase diagrams could be systematically biased for such systems.
- One could decompose the explicit anharmonic free energy by phonon mode (e.g., via temperature-dependent phonon self-energies) to identify which vibrations drive the stabilization of β-Sn, offering a microscopic design lever for alloying or control of tin pest.
- The classical-nuclei nature of the thermodynamic integration leaves a small unquantified nuclear-quantum correction; quantifying it (e.g., via path-integral methods) would sharpen the attribution of the 85 K shift purely to anharmonicity.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper constructs an Atomic Cluster Expansion (ACE) machine-learned interatomic potential for tin trained on PBE-DFT data and uses it to compute the free energies of α-Sn and β-Sn by two routes: the quasi-harmonic approximation (QHA, Eq. 1) and classical thermodynamic integration (TI, Eqs. 2–3). On the same ACE potential, the QHA transformation temperature is 377 K, whereas TI gives 288 K, in close agreement with the experimental 286 K. The authors attribute the 85 K shift to explicit phonon–phonon anharmonicity, reporting ΔF_anh ≈ −11 meV/atom for β-Sn at 300 K and a much smaller effect for α-Sn. The anharmonicity of β-Sn is further supported by heat capacity, velocity-autocorrelation spectra, and force-deviation metrics. The central claim is that explicit anharmonicity is essential for the phase stability of tin, while the absolute transition temperature remains limited by the underlying 0 K PBE energetics.
Significance. If the central attribution is quantitatively reliable, the work is significant: it demonstrates a transferable MLIP-based strategy for isolating explicit anharmonic contributions to a phase boundary and resolves a long-standing discrepancy for tin without fitting to the experimental transition temperature. The design is sound in that both free-energy routes use the same potential, so the comparison isolates methodological differences rather than potential errors. The microscopic corroboration from heat capacity, VACF spectra, and force residuals is a real strength and makes the qualitative conclusion—that β-Sn is strongly anharmonic—credible. However, the quantitative claim that the −11 meV/atom term is purely anharmonic requires an additional correction or justification because the two free energies being subtracted are not computed at the same statistical-mechanical level.
major comments (3)
- [Thermodynamic integration, Eqs. (2)–(3) versus QHA Eq. (1)] The subtraction ΔF_anh = F_TI − F_QHA does not isolate explicit anharmonicity alone. F_QHA (Eq. 1) is a quantum harmonic free energy with zero-point motion and Bose–Einstein occupation, while F_TI from classical MD is a classical free energy. Thus ΔF_anh = (F_TI_classical − F_classical_harmonic) − (F_QHA_quantum − F_classical_harmonic). The second bracket is a nuclear-quantum correction. For Sn's 3–5 THz modes at 286 K, the high-T per-mode correction is roughly (ħω)^2/(24 k_B T) ≈ 0.3–0.7 meV, and the summed phase-dependent correction could be ~1 meV/atom—exactly the quantity entering the claimed 85 K shift. Please quantify this correction (e.g., by adding the harmonic quantum correction to the TI result or by path-integral MD) or explicitly report ΔF_anh as the TI−QHA difference with the nuclear-quantum contribution separated. As written, the title and conclusions state that −11 meV/ato
- [Methods, Thermodynamic integration; Fig. 4] It is not specified whether the TI free energies are Helmholtz free energies at fixed 0 K volumes or Gibbs free energies minimized over volume. The QHA curves are explicitly min_V F(V,T) (Methods, QHA), while the TI description lists only temperatures and λ values. If TI was performed at a single volume per phase, the TI−QHA difference includes the thermal-expansion contribution to the free energy in addition to explicit anharmonicity. Please specify the volume protocol (e.g., NPT simulations, or multiple NVT volumes followed by minimization) and report the volumes used at each temperature. Without this, the comparison of the two free-energy curves is not well-defined.
- [Results and Fig. 4] No statistical uncertainties are reported for F_TI, ΔF_anh, or the resulting T_TI. The TI result is a numerical quadrature over λ combined with finite MD sampling, and the transition temperature is a crossing of two free-energy curves. The statement that TI gives 288 K, 'within 2 K of experiment,' has no error bar. Please report standard errors (e.g., from block averaging, multiple independent TI runs, or bootstrap over the integration) and propagate them to the transition temperature. This is necessary to assess whether the agreement with experiment is meaningful or fortuitous.
minor comments (6)
- [Abstract] The phrase 'existing theoretical approaches over- or underestimating' is grammatically awkward; suggest 'overestimating or underestimating' or a similar rewrite.
- [Fig. 6] The VACF spectra are shown at 10 K and 300 K. Since the transition is at 286 K, adding a spectrum just above Tαβ would more directly illustrate the claimed collapse of the 3–3.5 THz feature at the transition.
- [Methods, Eqs. (9)–(10)] The harmonic potential energy E_harm = (1/2)u^T Φ u uses displacements u from MD snapshots. Please clarify whether global translations/rotations are projected out before evaluating the harmonic energy, as these can artificially inflate the harmonic term.
- [Table 3] The 'Force cosine' column is not defined in the main text or caption. Since it is used as a metric for anharmonicity, a definition (e.g., average cosine of the angle between F_MD and F_harm) should be provided.
- [Title and Conclusions] The phrase 'Full Anharmonic Treatment' may be misleading because the TI/MD simulations treat nuclei classically. Consider adding the qualifier 'classical' or noting in the Conclusions that nuclear quantum effects are not included in the TI result.
- [Discussion] The comparison with Legrain and Manzhos (Ref. 6) is useful, but the quoted 1 meV/atom estimate is from a different electronic-structure and fitting context. A sentence noting the methodological differences beyond run length (e.g., reference potential and fixed-volume mapping) would improve the comparison.
Circularity Check
No significant circularity: the central result is a model-internal QHA-vs-TI prediction on a DFT-trained potential, not a fitted or self-citational reduction.
full rationale
The derivation chain is: (1) DFT-PBE training data are used to fit an ACE potential; (2) the same ACE potential is used for both quasi-harmonic phonon free energies (Eq. 1) and thermodynamic-integration free energies (Eqs. 2-3); (3) the transformation temperature is identified as the crossing of the two free-energy curves. No experimental transition temperature enters as a fitting target; the paper explicitly contrasts its approach with the Ravelo-Baskes MEAM, where 'the experimental value was used as an explicit fitting target.' The QHA and TI free energies are evaluated on the same potential energy surface by design, so their difference is a well-defined model-internal quantity; calling this difference the 'explicit anharmonic free energy' is an operational definition and a physical interpretation, not a circular derivation. The main substantive concern—that F_TI is computed classically while F_QHA is a quantum phonon free energy including zero-point and Bose-Einstein occupations—is a statistical-mechanics/accuracy limitation in the decomposition, not a circularity: no equation reduces to its own input, and no parameter is fitted to the resulting 288 K prediction. The self-citations to ACE/Pacemaker implementation references are methodological, not load-bearing uniqueness claims, and the paper's own caveat that 'we cannot exclude that the quantitative agreement benefits in part from a cancellation between an overestimated static energy and the anharmonic contribution' is a limitation on the absolute temperature, not evidence that the prediction is equivalent to the inputs used to construct the model. The agreement with the experimental 286 K is therefore a genuine prediction rather than a fitted or definitionally forced outcome.
Assumptions & free parameters
free parameters (1)
- ACE potential parameters =
Fitted to 4,427 DFT configurations; RMSE 5.7 meV/atom and 33.0 meV/Å on validation set
assumptions (4)
- domain assumption PBE-DFT with PAW pseudopotentials gives a sufficiently accurate 0 K energy difference and potential energy surface for tin
- domain assumption Classical thermodynamic-integration free energies are directly comparable to quantum quasi-harmonic free energies at 286 K
- domain assumption The ACE potential transfers accurately to finite-temperature anharmonic configurations of β-Sn, not just training configurations
- domain assumption Metastable α-Sn can be thermodynamically integrated above 286 K without transforming to β-Sn during the λ path
Cite this review
Pith. "Pith review of Accurate Prediction of the $\alpha \to \beta$ Phase Transformation Temperature in Tin via Full Anharmonic Treatment." pith.science (2026). https://pith.science/paper/7SJ64WJV
@misc{pith2026260725978,
author = {Pith},
title = {Pith review of: Accurate Prediction of the $\alpha \to \beta$ Phase Transformation Temperature in Tin via Full Anharmonic Treatment},
year = {2026},
howpublished = {\url{https://pith.science/paper/7SJ64WJV}},
note = {Machine review of arXiv:2607.25978}
}
abstract
Predicting the $\alpha \to \beta$ (grey-to-white) transition temperature in tin presents a longstanding challenge for atomistic simulations, with existing theoretical approaches over- or underestimating the experimental boundary (286 K) by up to several hundred Kelvin. In this work, we construct an Atomic Cluster Expansion (ACE) potential trained on density functional theory data to evaluate the finite-temperature free energies of both phases. Evaluated on the same potential energy surface, the quasi-harmonic approximation predicts a transformation temperature of 377 K, whereas full thermodynamic integration, which accounts for explicit vibrational anharmonicity, yields 288 K. This shift directly quantifies the explicit anharmonic free energy, which is substantial for metallic $\beta$-Sn but negligible for semiconducting $\alpha$-Sn. The anisotropic anharmonicity in $\beta$-Sn is corroborated by its excess heat capacity, temperature-driven renormalization of its vibrational spectrum, and deviations of its atomic forces and displacements from the harmonic reference. Our results demonstrate that capturing full lattice anharmonicity is essential for predicting the phase stability of tin, while the absolute transition temperature remains limited by the accuracy of the underlying 0 K energetics.
Figures
Figures from the paper (4 more)
Reference graph
Works this paper leans on
-
[1]
Farrow, R.et al.The growth of metastable, heteroepitaxial films of α-sn by metal beam epitaxy.Journal of Crystal Growth 54, 507–518 (1981)
1981
-
[2]
& Smith, R
Vnuk, F., De Monte, A. & Smith, R. The effect of pressure on the semiconductor-to-metal transition temperature in tin and in dilute sn-ge alloys.Journal of applied physics55, 4171–4176 (1984). 12/14
1984
-
[3]
Ihm, J. & Cohen, M. L. Equilibrium properties and the phase transition of grey and white tin.Phys. Rev. B23, 1576–1579 (1981). URLhttps://link.aps.org/doi/10.1103/PhysRevB.23.1576
-
[4]
& de Gironcoli, S
Pavone, P., Baroni, S. & de Gironcoli, S. α↔β phase transition in tin: A theoretical study based on density-functional perturbation theory.Phys. Rev. B57, 10421–10423 (1998). URL https://link.aps.org/doi/10.1103/ PhysRevB.57.10421
1998
-
[5]
& Park, C.-H
Na, S.-H. & Park, C.-H. First-principles study of structural phase transition of sn.Journal of the Korean Physical Society 56, 494–497 (2010)
2010
-
[6]
& Manzhos, S
Legrain, F. & Manzhos, S. Understanding the difference in cohesive energies between alpha and beta tin in dft calculations. AIP Advances6(2016)
2016
-
[7]
Friák, M.et al.Multi-methodological study of temperature trends in mössbauer effect in sn.Computational Materials Science215, 111780 (2022)
2022
-
[8]
Physical Review Materials7, 053603 (2023)
Chen, T.et al.Modeling the high-pressure solid and liquid phases of tin from deep potentials with ab initio accuracy. Physical Review Materials7, 053603 (2023). 9.Nitol, M. S.et al.Hybrid interatomic potential for Sn.Physical Review Materials7, 043601 (2023)
2023
Show all 33 references
-
[10]
J.et al.Tin-pest problem as a test of density functionals using high-throughput calculations.Physical Review Materials5, 083608 (2021)
Mehl, M. J.et al.Tin-pest problem as a test of density functionals using high-throughput calculations.Physical Review Materials5, 083608 (2021)
2021
-
[11]
Houben, K.et al.In situ study of the α-sn to β-sn phase transition in low-dimensional systems: Phonon behavior and thermodynamic properties.Phys. Rev. B100, 075408 (2019). URL https://link.aps.org/doi/10.1103/ PhysRevB.100.075408
2019
-
[12]
P., Burke, K
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple.Physical Review Letters77, 3865–3868 (1996)
1996
-
[13]
& Lee, M
Ko, W.-S., Kim, D.-H., Kwon, Y .-J. & Lee, M. Atomistic simulations of pure tin based on a new modified embedded-atom method interatomic potential.Metals8, 900 (2018)
2018
-
[14]
& Pasquarello, A
Haxhijaj, A., Riemelmoser, S. & Pasquarello, A. Lattice dynamics and structural phase stability of group-iv elemental solids with the r2SCAN functional.Phys. Rev. B113, 104105 (2026). URL https://link.aps.org/doi/10. 1103/znkf-3g37
2026
-
[15]
Khvan, A.et al.Thermodynamic properties of tin: Part I experimental investigation, ab-initio modelling of α-, β-phase and a thermodynamic description for pure metal in solid and liquid state from 0 K.Calphad65, 50–72 (2019)
2019
-
[16]
& Järvinen, M
Merisalo, M. & Järvinen, M. Anharmonicity of lattice vibrations in white tin by x-ray measurement of almost-forbidden reflections.Philosophical Magazine B37, 233–240 (1978)
1978
-
[17]
& Baskes, M
Ravelo, R. & Baskes, M. Equilibrium and thermodynamic properties of grey, white, and liquid tin.Physical Review Letters 79, 2482–2485 (1997)
1997
-
[18]
R.et al.Structural and dynamic properties of liquid tin from a new modified embedded-atom method force field
Vella, J. R.et al.Structural and dynamic properties of liquid tin from a new modified embedded-atom method force field. Physical Review B95, 064202 (2017)
2017
-
[19]
A., Baskes, M
Etesami, S. A., Baskes, M. I., Laradji, M. & Asadi, E. Thermodynamics of solid sn and pbsn liquid mixtures using molecular dynamics simulations.Acta Materialia161, 320–330 (2018)
2018
-
[20]
Chesser, I.et al.The structure and migration of twin boundaries in tetragonal β-Sn: An application of machine learning based interatomic potentials.Acta Materialia299, 121390 (2025)
2025
-
[21]
Atomic cluster expansion for accurate and transferable interatomic potentials.Physical Review B99, 014104 (2019)
Drautz, R. Atomic cluster expansion for accurate and transferable interatomic potentials.Physical Review B99, 014104 (2019)
2019
-
[22]
H., Forslund, A., Srinivasan, P
Jung, J. H., Forslund, A., Srinivasan, P. & Grabowski, B. Dynamically stabilized phases with full ab initio accuracy: Thermodynamics of Ti, Zr, Hf with a focus on the hcp-bcc transition.Physical Review B108(2023)
2023
-
[23]
L., Rowe, J
Price, D. L., Rowe, J. M. & Nicklow, R. M. Lattice dynamics of grey tin and indium antimonide.Phys. Rev. B3, 1268–1279 (1971). URLhttps://link.aps.org/doi/10.1103/PhysRevB.3.1268
1971 doi
-
[24]
M., Brockhouse, B
Rowe, J. M., Brockhouse, B. N. & Svensson, E. C. Lattice dynamics of white tin.Phys. Rev. Lett.14, 554–556 (1965). URLhttps://link.aps.org/doi/10.1103/PhysRevLett.14.554. 25.Zdetsis, A. D. On the lattice dynamics of grey tin.Journal of Physics and Chemistry of Solids38, 1113–1...
1965 doi
-
[26]
Rayne, J. A. & Chandrasekhar, B. S. Elastic constants of β-tin from 4.2 K to 300 K.Physical Review122, 1962–1962 (1961)
1962
-
[27]
& Hafner, J
Kresse, G. & Hafner, J. Ab initiomolecular dynamics for open-shell transition metals.Physical Review B48, 13115–13118 (1993)
1993
-
[28]
& Furthmüller, J
Kresse, G. & Furthmüller, J. Efficient iterative schemes forab initiototal-energy calculations using a plane-wave basis set. Physical Review B54, 11169–11186 (1996)
1996
-
[29]
& Joubert, D
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method.Physical Review B59, 1758–1775 (1999)
1999
-
[30]
Lysogorskiy, Y .et al.Performant implementation of the atomic cluster expansion (PACE) and application to copper and silicon.npj Computational Materials7(2021)
2021
-
[31]
& Drautz, R
Lysogorskiy, Y ., Bochkarev, A., Mrovec, M. & Drautz, R. Active learning strategies for atomic cluster expansion models. Physical Review Materials7, 043801 (2023)
2023
-
[32]
P.et al.LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales.Computer Physics Communications271, 108171 (2022)
Thompson, A. P.et al.LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales.Computer Physics Communications271, 108171 (2022)
2022
-
[33]
R.et al.Kokkos 3: Programming model extensions for the exascale era.IEEE Transactions on Parallel and Distributed Systems33, 805–817 (2022)
Trott, C. R.et al.Kokkos 3: Programming model extensions for the exascale era.IEEE Transactions on Parallel and Distributed Systems33, 805–817 (2022)
2022
-
[34]
& Tanaka, I
Togo, A., Chaput, L., Tadano, T. & Tanaka, I. Implementation strategies in phonopy and phono3py.J. Phys. Condens. Matter35, 353001 (2023). 35.Togo, A. First-principles phonon calculations with phonopy and phono3py.J. Phys. Soc. Jpn.92, 012001 (2023)
2023
-
[36]
Materials and technologies for sustainable development
Menon, S., Lysogorskiy, Y ., Rogal, J. & Drautz, R. Automated free-energy calculation from atomistic simulations.Phys. Rev. Mater .5, 103801 (2021). URLhttps://link.aps.org/doi/10.1103/PhysRevMaterials.5.103801. Acknowledgements MF acknowledges financial support from the Czech...
2021 doi
Reviewed August 1, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.