pith. sign in

arxiv: 2606.24097 · v1 · pith:X33OJ7JRnew · submitted 2026-06-23 · ❄️ cond-mat.mtrl-sci

Reversible non-equilibrium phase transformation in amorphous germanium

Pith reviewed 2026-06-25 23:45 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords germaniumamorphousliquid-glass transitionfirst-order phase transitionfree energymolecular dynamicsmetastable phaseshysteresis
0
0 comments X

The pith

First-principles simulations show amorphous germanium undergoes a reversible first-order transition to a metastable liquid below the melting point.

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

The paper simulates germanium with molecular dynamics and observes a reversible switch between its amorphous solid and liquid forms at temperatures below the normal melting point, accompanied by a wide hysteresis loop. Direct computation of free energies for both states, made possible by absolute entropy calculations, establishes that the switch occurs as a discontinuous first-order phase change between two metastable phases. This finding supplies thermodynamic grounding for models in which amorphous germanium thin films recrystallize explosively by passing through a transient liquid before reaching the stable crystal. The result matters for understanding how non-equilibrium phases can interconvert without reaching the equilibrium crystal structure.

Core claim

Direct calculation of the liquid and amorphous state free energies, enabled by absolute entropy calculations, verify that the transition is first order in character between two metastable phases.

What carries the argument

Absolute entropy calculations within first-principles molecular dynamics simulations that permit direct free-energy comparison between the liquid and amorphous states.

If this is right

  • The calculated free-energy ordering supports models of explosive recrystallization in amorphous Ge films that proceed through a metastable liquid intermediate.
  • Both the amorphous solid and the liquid remain metastable relative to the crystal below the equilibrium melting temperature.
  • The wide hysteresis loop indicates kinetic barriers that allow the liquid-glass transition to be observed reversibly on simulation timescales.
  • The same absolute-entropy method supplies a route to map non-equilibrium phase boundaries for other covalently bonded amorphous materials.

Where Pith is reading between the lines

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

  • Similar absolute-entropy techniques could test whether other group-IV amorphous semiconductors also host first-order liquid-glass lines.
  • The existence of a first-order line between two metastable phases implies that glass-transition diagrams may contain sharp boundaries rather than purely continuous crossovers.
  • Calorimetric or in-situ diffraction experiments that track latent heat or order-parameter jumps during rapid heating of amorphous Ge would provide direct experimental tests.

Load-bearing premise

The first-principles molecular dynamics simulations and absolute entropy calculations provide accurate thermodynamic descriptions of the liquid and amorphous states.

What would settle it

A continuous rather than discontinuous change in the free-energy difference or structural order parameter across the observed transition temperature would falsify the first-order claim.

Figures

Figures reproduced from arXiv: 2606.24097 by Marek Mihalkovi\v{c}, Michael Widom, Yang Huang.

Figure 1
Figure 1. Figure 1: FIG. 1. Simulated structure of [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Hysteresis loops in (a) energy and (b) volume for [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Plots of [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: shows our results for several temperatures. Further details of these calculations are in the Appendix [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
read the original abstract

First principles molecular dynamics simulations of germanium reveal a reversible liquid-glass transition below the equilibrium melting point with a wide hysteresis loop. Direct calculation of the liquid and amorphous state free energies, enabled by absolute entropy calculations, verify that the transition is first order in character between two metastable phases. These results lend credence to models of explosive recrystallization from amorphous Ge thin films, through a metastable liquid state and finally reaching the low temperature crystalline structure.

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

1 major / 0 minor

Summary. The manuscript reports first-principles molecular dynamics simulations of germanium that identify a reversible liquid-glass transition below the equilibrium melting point, accompanied by a wide hysteresis loop. Absolute entropy calculations are used to compute the free energies of the liquid and amorphous states directly, confirming that the transition is first-order in character between two metastable phases. The results are presented as support for models of explosive recrystallization in amorphous Ge thin films proceeding through a metastable liquid intermediate.

Significance. If the underlying calculations hold, the work supplies direct thermodynamic evidence distinguishing a first-order transition from purely kinetic glass formation in a semiconductor system. The absolute-entropy route to free-energy comparison is a methodological strength that avoids reliance on fitting or reference states, lending weight to the interpretation of the observed hysteresis as a true phase coexistence between metastable liquid and amorphous phases.

major comments (1)
  1. [Abstract] Abstract and Methods (implied): the central claim that absolute-entropy-enabled free-energy calculations verify a first-order transition rests on the accuracy of the underlying first-principles MD trajectories. However, no information is supplied on the DFT functional, pseudopotentials, plane-wave cutoff, system sizes, equilibration times, or statistical uncertainties in the entropy and free-energy values. These omissions prevent assessment of whether the reported hysteresis and free-energy crossing are robust or sensitive to technical choices.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the detailed review and constructive feedback. The single major comment highlights a genuine omission in the manuscript that we will correct. Our point-by-point response follows.

read point-by-point responses
  1. Referee: [Abstract] Abstract and Methods (implied): the central claim that absolute-entropy-enabled free-energy calculations verify a first-order transition rests on the accuracy of the underlying first-principles MD trajectories. However, no information is supplied on the DFT functional, pseudopotentials, plane-wave cutoff, system sizes, equilibration times, or statistical uncertainties in the entropy and free-energy values. These omissions prevent assessment of whether the reported hysteresis and free-energy crossing are robust or sensitive to technical choices.

    Authors: We agree that the manuscript as submitted lacks the necessary technical specifications. In the revised version we will add a dedicated Computational Methods section that reports the DFT functional, pseudopotentials, plane-wave cutoff, system sizes, equilibration and production run lengths, and the statistical uncertainties (obtained via block averaging) on the entropy and free-energy values. These additions will enable readers to judge the robustness of the reported hysteresis and free-energy crossing. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The derivation rests on direct first-principles MD simulations yielding absolute entropies and free energies for liquid and amorphous Ge phases. These outputs are computed from the underlying potential and trajectories without reduction to fitted parameters, self-definitions, or self-citation chains. The first-order character follows from comparing the independently obtained G(T) curves; no equation or premise collapses to its own input by construction. This is the standard, self-contained computational route to the claimed result.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Review is abstract-only; ledger reflects assumptions stated or implied in the abstract with no access to full methods or parameter details.

axioms (2)
  • domain assumption First-principles molecular dynamics simulations accurately capture the atomic interactions and dynamics in both liquid and amorphous germanium.
    Invoked as the basis for observing the transition and hysteresis.
  • domain assumption Absolute entropy calculations can be reliably performed to enable direct free-energy comparison between the liquid and amorphous states.
    Required to verify the first-order character of the transition.

pith-pipeline@v0.9.1-grok · 5589 in / 1097 out tokens · 30301 ms · 2026-06-25T23:45:12.167431+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

36 extracted references

  1. [1]

    Pauling, The structure and entropy of ice and of other crystals with some randomness of atomic arrangement, J

    L. Pauling, The structure and entropy of ice and of other crystals with some randomness of atomic arrangement, J. Am. Chem. Soc.57, 2680 (1935)

  2. [2]

    P. H. Poole, F. Sciortino, U. Essmann, and H. E. Stanley, Phase behaviour of metastable water, Nature360, 324 (1992)

  3. [3]

    Galloet al., Water: A tale of two liquids, Chem

    P. Galloet al., Water: A tale of two liquids, Chem. Rev. 116, 7463 (2016)

  4. [4]

    van Thiel and F

    M. van Thiel and F. H. Ree, High-pressure liquid-liquid phase change in carbon, Phys. Rev. B48, 3591 (1993)

  5. [5]

    M. W. C. Dharma-Wardana and D. D. Klug, Evidence for multiple liquid–liquid phase transitions in carbon, and the friedel ordering of its liquid state, Phys. Plasmas29, 022108 (2022)

  6. [6]

    Sastry and C

    S. Sastry and C. A. Angell, Liquid–liquid phase transition in supercooled silicon, Nat. Mat.2, 739 (2003)

  7. [7]

    Ganesh and M

    P. Ganesh and M. Widom, Liquid-liquid transition in su- percooled silicon determined by first-principles simula- tion, Phys. Rev. Lett.102, 075701 (2009)

  8. [8]

    M. H. Bhat, E. S. V. Molinero, V. C. Solomon, S. Sas- try, J. L. Yarger, and C. A. Angel, Vitrification of a monatomic metallic liquid, Nature448, 787 (2007)

  9. [9]

    L. Xu, Z. Wang, J. Chen, S. Chen, W. Yang, Y. Ren, X. Zuo, J. Zeng, Q. Wu, and H. Sheng, Folded network and structural transition in molten tin, Nat. comm.13, 126 (2022)

  10. [10]

    Mishima, L

    O. Mishima, L. D. Calvert, and E. Whalley, An ap- parently first-order transition between two amorphous phases of ice induced by pressure, Nature314, 76 (1985)

  11. [11]

    P. F. McMillan, M. Wilson, D. Daisenberger, and D. Ma- chon, A density-driven phase transition between semicon- ducting and metallic polyamorphs of silicon, Nat. Mater. 4, 680 (2005)

  12. [12]

    Coppari, J

    F. Coppari, J. C. Chervin, A. Congeduti, M. Lazzeri, A. Polian, E. Principi, and A. Di Cicco, Pressure- induced phase transitions in amorphous and metastable crystalline germanium by raman scattering, x-ray spec- troscopy, and ab initio calculations, Phys. Rev. B80, 115213 (2009)

  13. [13]

    J. Tauc, R. Grigorovici, and A. Vancu, Optical properties and electronic structure of amorphous germanium, Phys. Stat. Solidi (1966)

  14. [14]

    explosively

    T. Takamori, R. Messier, and R. Roy, New noncrystalline germanium which crystallizes “explosively” at room tem- perature, Appl. Phys. Lett.20, 201 (1972)

  15. [15]

    R. L. Chapman, J. C. C. Fan, H. J. Zeiger, and R. P. Gale, Crystallization-front velocity during scanned laser crystallization of amorphous Ge films, Appl. Phys. Lett. 37, 292 (1980)

  16. [16]

    H. J. Leamy, W. L. Brown, G. K. Celler, G. Foti, G. H. Gilmer, and J. C. C. Fan, Explosive crystallization of amorphous germanium, Appl. Phys. Lett.38, 137 (1981)

  17. [17]

    Nikolova, M

    L. Nikolova, M. J. Stern, J. M. MacLeod, B. W. Reed, H. Ibrahim, G. H. Campbell, F. Rosei, T. LaGrange, and B. J. Siwick, In situ investigation of explosive crys- tallization in a-Ge: Experimental determination of the interface response function using dynamic transmission electron microscopy, J. Appl. Phys.116, 093512 (2014)

  18. [18]

    B. G. Bagley and H. S. Chen, A calculation of the thermo- dynamic first order amorphous semiconductor to metal- lic liquid transition temperature, AIP Conf. Proc.50, 97 (1979)

  19. [19]

    Kresse and J

    G. Kresse and J. Hafner, Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium, Phys. Rev. B49, 14251 (1994)

  20. [20]

    J. W. Gibbs,Elementary principles in statistical mechan- ics(Dover, New York, 1960)

  21. [21]

    C. E. Shannon, A Mathematical Theory of Communica- tion, Bell System Tech. J.27, 379 (1948)

  22. [22]

    R. E. Nettleton and M. S. Green, Expression in Terms of Molecular Distribution Functions for the Entropy Den- sity in an Infinite System, J. Chem. Phys.29, 1365 (1958)

  23. [23]

    H. J. Raveche, Entropy and Molecular Correlation Func- tions in Open Systems. I. Derivation, J. Chem. Phys.55, 2242 (1971)

  24. [24]

    Baranyai and D

    A. Baranyai and D. J. Evans, Direct entropy calculation from computer simulation of liquids, Phys. Rev. A40, 3817 (1989)

  25. [25]

    M. C. Gao and M. Widom, Information entropy of liquid metals, J. Phys. Chem. B122, 3550 (2018)

  26. [26]

    Widom and M

    M. Widom and M. Gao, First principles calculation of the entropy of liquid aluminum, Entropy21(2019)

  27. [27]

    Huang and M

    Y. Huang and M. Widom, Entropy approximations for simple fluids, Phys. Rev. E109, 034130 (2024)

  28. [28]

    L. F. Richardson, The approximate arithmetical solution by finite differences of physical problems involving dif- ferential equations, with an application to the stresses in a masonry dam, Phil. Trans. R. Soc. Lond. A210, 307 (1911)

  29. [29]

    Widom, High entropy alloys: fundamentals and ap- plications (Springer, Berlin, 2015) Chap

    M. Widom, High entropy alloys: fundamentals and ap- plications (Springer, Berlin, 2015) Chap. 8. Prediction of structure and phase transformations

  30. [30]

    R. L. C. Vink and G. T. Barkema, Configurational en- tropy of network-forming materials, Phys. Rev. Lett.89, 076405 (2002)

  31. [31]

    Kresse and D

    G. Kresse and D. Joubert, From ultrasoft pseudopoten- tials to the projector augmented-wave method, Phys. Rev. B59, 1758 (1999)

  32. [32]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)

  33. [33]

    Baldereschi, Mean-value point in the Brillouin zone, Phys

    A. Baldereschi, Mean-value point in the Brillouin zone, Phys. Rev. B7, 5212 (1973)

  34. [34]

    I. S. Galtsov, V. B. Fokin, A. V. Dorovatovsky, M. A. Paramonov, G. S. Demyanov, D. V. Minakov, M. A. Sheindlin, and P. R. Levashov, Thermophysical proper- ties of solid and liquid nickel near melting point, J. Appl. Phys.136, 145104 (2024)

  35. [35]

    Jinnouchi, F

    R. Jinnouchi, F. Karsai, and G. Kresse, On-the-fly ma- chine learning force field generation: Application to melt- ing points, Phys. Rev. B100, 014105 (2019)

  36. [36]

    B. D. McKay and A. Piperno, Practical graph isomor- phism, ii, J. Symb. Comp.60, 94 (2014). 6 Appendix A: End Matter first principles calculations—Our molecular dynamics and phonon calculations utilize the VASP code [31] in the PBE generalized gradient approximation [32] using the supplied projector augmented wave potentials for Ge with valence 4 at the d...