{"id":"007a857e-b790-4426-9af9-8b6411ed6cb0","arxiv_id":"2607.25978","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Full anharmonic free-energy integration lowers the predicted α→β transition temperature of tin from 377 K (quasi-harmonic) to 288 K, matching experiment, because β-Sn carries a large explicit anharmonic free-energy contribution.","lead":"Using a machine-learning potential trained on density-functional-theory data, the authors compare quasi-harmonic and fully anharmonic free energies for the two solid phases of tin and obtain a grey-to-white transition temperature of 288 K, close to the measured 286 K. The central message is that anharmonic atomic vibrations in metallic white tin shift the phase boundary by about 85 K, so they cannot be ignored.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Classical TI vs quantum QHA: ΔF_anh mixes anharmonicity with nuclear quantum effects, so the 85 K shift is not purely anharmonic.","rationale":"The reader's weakest assumption correctly identifies the quantum-versus-classical mismatch between TI and QHA as the most load-bearing gap in the central claim. My analysis confirms that ΔF_anh, as defined, includes a nuclear quantum correction that is not accounted for anywhere in the paper. For Sn near 286 K, this correction is likely on the order of 0.3–1 meV/atom in the phase free-energy difference, which translates to roughly 10–30 K in transition temperature—a non-negligible fraction of the 85 K shift the authors attribute exclusively to anharmonicity. The paper's own corroborating evidence (excess heat capacity, VACF spectral collapse, anisotropic displacement distributions, and large residual force amplitude AF=0.58 in β-Sn) provides strong independent support for the qualitative conclusion that β-Sn is significantly anharmonic. Thus the concern does not invalidate the paper's central message; it only weakens the precise quantitative statement that the full 85 K shift is a direct measure of explicit anharmonicity. The authors also honestly acknowledge the PBE 0 K energy uncertainty and possible cancellation, which is a point in their favor. The proposed test—recomputing QHA with classical statistics, or running PIMD—would settle whether the quantum correction is small enough to ignore. Given the paper is already CONDITIONAL in the reader's verdict, no change is warranted; the conditionality should explicitly include this quantum-correction test.","tokens_in":13201,"tokens_out":9053,"duration_ms":91417,"concrete_test":"Recompute F_QHA for both phases using classical phonon statistics, i.e., replace the quantum occupation in Eq. 1 with kBT Σ ln(ħω/kBT) at the same volumes as in the QHA minimization, giving F_QHA_classical(T). Then compute the nuclear quantum correction ΔNQE = F_QHA_classical − F_QHA for each phase. If the β−α difference of ΔNQE at 286 K exceeds ~1 meV/atom, the reported ΔF_anh is substantially contaminated by quantum effects; if it is below ~0.3 meV/atom, the anharmonic attribution is safe. Alternatively, run path-integral MD thermodynamic integration at 286 K and compare with the classical TI result.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative claim is that the difference ΔF_anh = F_TI − F_QHA isolates the explicit anharmonic free energy, and that this anharmonic term alone lowers the transition temperature from 377 K to 288 K. However, the two free energies being subtracted are not computed at the same level of statistical mechanics. F_QHA (Eq. 1) is a quantum harmonic free energy including zero-point motion and Bose–Einstein occupation, while F_TI (Eqs. 2–3) is a classical thermodynamic integration from an MD simulation. Therefore ΔF_anh = (F_TI − F_classical-harmonic) + (F_classical-harmonic − F_QHA), where the first bracket is the true explicit anharmonic free energy (classical) and the second is a nuclear quantum correction. At 286 K, Sn's Debye temperature is only ~200 K, so the quantum correction is not negligible: for typical β-Sn modes (3–5 THz), the high-T expansion gives per-mode corrections of order (ħω)²/(24kBT) ≈ 0.3–0.7 meV/atom, and the phase difference could be ~1 meV/atom. This enters exactly the 85 K shift attributed entirely to anharmonicity, so the reported −11 meV/atom stabilization of β-Sn is likely an overestimate of the anharmonic contribution. The qualitative conclusion that β-Sn is strongly anharmonic is independently supported by heat capacity, VACF, and force metrics, so the core claim may survive, but the quantitative attribution and the 'within 2 K' agreement are not as clean as stated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":13537,"tokens_out":6297,"duration_ms":67558,"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":[{"comment":"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","section":"Thermodynamic integration, Eqs. (2)–(3) versus QHA Eq. (1)"},{"comment":"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.","section":"Methods, Thermodynamic integration; Fig. 4"},{"comment":"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.","section":"Results and Fig. 4"}],"minor_comments":[{"comment":"The phrase 'existing theoretical approaches over- or underestimating' is grammatically awkward; suggest 'overestimating or underestimating' or a similar rewrite.","section":"Abstract"},{"comment":"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.","section":"Fig. 6"},{"comment":"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.","section":"Methods, Eqs. (9)–(10)"},{"comment":"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.","section":"Table 3"},{"comment":"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.","section":"Title and Conclusions"},{"comment":"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.","section":"Discussion"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful thing here is a clean internal comparison: the authors build one ACE potential for Sn, then compute the α→β transition temperature two ways on the same potential energy surface. QHA gives 377 K, full thermodynamic integration gives 288 K, close to the experimental 286 K. The difference translates to an explicit anharmonic free energy of about −11 meV/atom for β-Sn at 300 K, with α-Sn staying nearly harmonic. That direct decomposition is new for the ambient-pressure Sn transition, and it makes a strong case that the QHA overshoot is not just a functional artifact.\n\nWhat the paper does well: the anharmonic signature is corroborated from several angles—β-Sn's heat capacity exceeding the Dulong–Petit limit, the temperature-driven collapse of its 3–3.5 THz VACF peak, the force residuals, and the anisotropic displacement distributions. No experimental transition temperature enters the fitting, so the 2 K agreement is a genuine prediction, not a fit. The potential itself looks well validated against DFT phonons and elastic constants.\n\nThe soft spot is the one the stress-test note flags: the TI free energy is classical, the QHA free energy is quantum, and the difference is treated as pure anharmonicity. Sn's Debye temperature is around 200 K, so nuclear quantum effects are not negligible at 300 K; the per-mode correction could be a few meV/atom, and the 85 K shift is therefore not purely anharmonic. The paper should say this explicitly, or better, estimate the quantum correction. Also missing: any statistical uncertainty on the free-energy curves or transition temperatures. The authors do acknowledge in the Discussion that the excellent agreement with experiment may partly come from cancellation with the PBE 0 K energy error, which is honest, but it means the 288 K number should be read as illustrative, not as a rigorous quantitative prediction.\n\nThese caveats do not break the core qualitative result. The evidence that β-Sn is far more anharmonic than α-Sn is independent of the classical/quantum bookkeeping, and the ordering of the transition temperature shift is robust. This paper deserves a serious referee and publication after the authors address the nuclear-quantum point and report error bars. I'd want it cited in future work on anharmonic phase stability, and I'd bring it to a reading group discussing MLIP free-energy methods.","headline":"Solid paper showing beta-Sn is strongly anharmonic, but the headline number mixes classical anharmonicity with quantum effects—qualitatively right, quantitatively overstated.","tokens_in":14037,"tokens_out":1853,"would_cite":true,"duration_ms":21509,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["64.70.K-","63.20.Ry","65.40.-b","71.15.Mb"],"model":"deepseek-v4-flash","headline":"Anharmonicity, not the functional, sets tin's α→β transition near 288 K.","keywords":["tin","α-β phase transition","anharmonicity","quasi-harmonic approximation","thermodynamic integration","machine-learned interatomic potential","atomic cluster expansion","free energy"],"falsifier":"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","tokens_in":13084,"feed_emoji":"🌡️","tokens_out":3789,"duration_ms":34171,"temperature":0.7,"pith_summary":"The paper argues that the long-standing overestimate of tin's grey-to-white (α→β) transition temperature in first-principles calculations comes mainly from neglected lattice anharmonicity in metallic β-Sn, not from errors in the electronic-structure functional. Using a machine-learned Atomic Cluster Expansion (ACE) potential trained on PBE density-functional data, the authors evaluate free energies of both phases two ways on the same potential energy surface: the quasi-harmonic approximation (QHA) and full thermodynamic integration (TI). QHA predicts a transition at 377 K; TI predicts 288 K, within 2 K of the experimental 286 K. The difference directly quantifies an explicit anharmonic free energy of about −11 meV/atom for β-Sn at room temperature, while α-Sn remains nearly harmonic. If correct, the result shows that capturing explicit phonon–phonon interactions is essential for predicting phase boundaries in soft, polymorphic metals, and it shifts the focus from exchange-correlation functionals to vibrational treatment.","feed_headline":"Anharmonicity, not the functional, sets tin's 286 K transition","feed_subtitle":"Switching from quasi-harmonic to full thermodynamic integration shifts the predicted tin transition by 89 K, matching experiment.","key_machinery":"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","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Full anharmonicity fixes tin transition prediction","Anharmonicity brings tin's predicted transition to 288 K","Tin's α→β transition demands full anharmonicity","Explicit anharmonicity narrows tin's transition gap to 2 K","Quasi-harmonic fails; full anharmonicity wins tin's 286 K"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Full anharmonicity fixes tin transition prediction","Anharmonicity brings tin's predicted transition to 288 K","Tin's α→β transition demands full anharmonicity","Explicit anharmonicity narrows tin's transition gap to 2 K","Quasi-harmonic fails; full anharmonicity wins tin's 286 K"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000774,"raw_usage":{"total_tokens":3274,"prompt_tokens":772,"completion_tokens":2502,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":516,"completion_tokens_details":{"reasoning_tokens":2428}},"tokens_in":516,"tokens_out":2502,"duration_ms":16374,"temperature":1.0,"reasoning_tokens":2428,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T00:55:49.669417+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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","supporting_citations":[],"review_version":1}