pith. sign in

arxiv: 1906.10605 · v1 · pith:FHWYPD5Ynew · submitted 2019-06-25 · ❄️ cond-mat.mtrl-sci · cond-mat.dis-nn

van der Waals forces stabilize low-energy polymorphism in B2O3: Implications for the crystallization anomaly

Pith reviewed 2026-05-25 16:24 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci cond-mat.dis-nn
keywords boron oxidepolymorphismvan der Waals forcesglass formationcrystallization anomalydensity functional theoryrandom phase approximation
0
0 comments X

The pith

Van der Waals forces make the known B2O3 crystal lower in energy than predicted competing polymorphs.

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

The paper computes cohesive energies of B2O3 structures with density functional theory plus high-accuracy many-body methods to compare the known crystal against never-synthesized polymorphs. It finds that dispersion forces are required to place the experimental B2O3-I structure at the bottom of the energy ladder, with other candidates only a few kcal/mol higher. Those metastable crystals turn out to match the glass in both energy and density while remaining mechanically soft and anisotropic. The result reframes the material's strong tendency to form glass rather than crystals as a consequence of near-degenerate low-energy polymorphs.

Core claim

The van der Waals forces play a key role in making the experimentally known polymorph (B2O3-I) the lowest in energy, with many competing metastable structures lying only a few kcal/mol above. All metastable crystals are comparable in energy and density to the glass, while having anisotropic and mechanically soft structures. The best metastable polymorph shares a structural motif found in both the glass and a recently synthesized borosulfate compound.

What carries the argument

Cohesive-energy comparisons from random-phase approximation and quantum Monte Carlo calculations that include dispersion forces, applied to a set of recently predicted B2O3 polymorph structures.

If this is right

  • Metastable polymorphs sit close enough in energy and density to the glass to compete with it during cooling.
  • The lowest-lying metastable structure shares a local motif already seen in both the glass and a borosulfate compound.
  • Anisotropic and soft character of the metastable crystals supplies a structural reason for kinetic hindrance to crystallization.
  • Inclusion of van der Waals forces is required to recover the experimental stability ordering among B2O3 polymorphs.

Where Pith is reading between the lines

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

  • Similar dispersion-driven near-degeneracies may operate in other network-forming oxides that also vitrify easily.
  • Targeted synthesis routes could be designed around the shared structural motif of the lowest metastable polymorph.
  • Mechanical softness of the competing crystals suggests they could be stabilized by epitaxial strain or pressure.

Load-bearing premise

The recently predicted but unsynthesized polymorph structures are the relevant competitors whose energies must be compared to the glass and to B2O3-I.

What would settle it

Observation or calculation of any additional B2O3 polymorph whose energy lies below B2O3-I once van der Waals contributions are included, or a reversal of the stability ordering when dispersion forces are omitted.

Figures

Figures reproduced from arXiv: 1906.10605 by Fran\c{c}ois-Xavier Coudert, Francesco Mauri, Guillaume Ferlat, Henri Hay, Maria Hellgren, Michele Casula.

Figure 1
Figure 1. Figure 1: Left: 3x3x3 supercell of B2O3-I and T10 seen in the (a,c) plane. Right: layers of T0 and T0-0.5b, the latter based on a mixed decoration of BO3 triangle and boroxol units. All the polymorphs considered in this work are 3D networks with the exception of T0, T0-0.5b, and T0-b which are layered structures. II. METHODS Highly accurate QMC simulations have been carried out on three polymorphs (B2O3-I and two pr… view at source ↗
Figure 2
Figure 2. Figure 2: Density and lattice parameters of B2O3-I from different ab initio schemes. Errors on lattice parameters are expressed relatively to the experimental values [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Density and lattice parameters of T0 and T10 from [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 6
Figure 6. Figure 6: Corrections to the PBE density (left panel) and rel [PITH_FULL_IMAGE:figures/full_fig_p004_6.png] view at source ↗
Figure 5
Figure 5. Figure 5: Energies, relative to B2O3-I (in kcal mol−1 ) obtained from the different ab initio schemes, including RPA, for three polymorphs (T0-0.5b, T0-b, and T3). The energies for T0 and T10 are reported in [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: Left: Energy, relative to B2O3-I, and density for all polymorphs using DFT-D2 (red) and PBE (gray symbols). Symbols refer to the relative proportions of structural units (triangles:boroxol) in the polymorphs: 4 only BO3 triangles (1:0), 9 only boroxol rings (0:1), mixed decoration (3:1), ♦ mixed decoration (1:1). Experimental data (liquid and glass) are from Ref. [51]. Right: energy and bulk modulus using … view at source ↗
read the original abstract

The cohesive energies and structural properties of recently predicted, and never synthesized, B$_2$O$_3$ polymorphs are investigated from first principles using density functional theory and high-accuracy many-body methods, namely, the random phase approximation and quantum Monte Carlo. We demonstrate that the van der Waals forces play a key role in making the experimentally known polymorph (B$_2$O$_3$-I) the lowest in energy, with many competing metastable structures lying only a few kcal/mol above. Remarkably, all metastable crystals are comparable in energy and density to the glass, while having anisotropic and mechanically soft structures. Furthermore, the best metastable polymorph according to our stability criteria has a structural motif found in both the glass and a recently synthesized borosulfate compound. Our findings provide new perspectives for understanding the B$_2$O$_3$ anomalous behavior, namely, its propensity to vitrify in a glassy structure drastically different from the known crystal.

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 uses density functional theory supplemented by random phase approximation and quantum Monte Carlo calculations to compute cohesive energies of B2O3 polymorphs. It claims that inclusion of van der Waals forces reverses the energy ordering to make the experimentally known B2O3-I the ground state, with several unsynthesized metastable polymorphs lying only a few kcal/mol higher—energies and densities comparable to the glass—while exhibiting soft, anisotropic structures. The work links this near-degeneracy to the material’s propensity for vitrification rather than crystallization.

Significance. If the reported energy ordering holds, the results supply a first-principles rationale for the long-standing crystallization anomaly in B2O3 by showing how dispersion forces stabilize the known crystal while leaving multiple metastable competitors nearly degenerate with the amorphous phase. The explicit use of RPA and QMC for dispersion-dominated energy differences on the scale of a few kcal/mol is a methodological strength that goes beyond standard DFT. The identification of a structural motif shared with both the glass and a recently synthesized borosulfate compound offers a concrete, testable connection between crystal and amorphous states.

major comments (2)
  1. [§3 (selection of polymorphs)] §3 (selection of polymorphs): The central claim that B2O3-I becomes the lowest-energy structure once vdW is included rests on the assumption that the finite set of recently predicted polymorphs exhausts the relevant low-energy competitors. No global structure search, exhaustive enumeration, or explicit argument for completeness of this set is provided; any missed lower-energy structure would alter the reported stability ordering and the interpretation of the crystallization anomaly.
  2. [§4–5 (energy differences and convergence)] §4–5 (energy differences and convergence): The key result is that vdW corrections reverse the ordering by only a few kcal/mol. However, the manuscript does not report basis-set convergence tests, k-point sampling details, or statistical error bars from the QMC runs that would establish whether the reported differences are robust to the level of numerical precision required for such small energy scales.
minor comments (2)
  1. [Table 1] Table 1: the caption should explicitly state whether the listed densities are calculated with or without vdW corrections, as this affects direct comparison to the glass density.
  2. [Abstract and Table 2] The abstract states that metastable crystals are 'comparable in energy and density to the glass'; a quantitative statement of the density differences (e.g., in a dedicated column of Table 2) would strengthen this claim.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their constructive comments, which help clarify the scope and robustness of our findings. We respond to each major comment below.

read point-by-point responses
  1. Referee: [§3 (selection of polymorphs)] The central claim that B2O3-I becomes the lowest-energy structure once vdW is included rests on the assumption that the finite set of recently predicted polymorphs exhausts the relevant low-energy competitors. No global structure search, exhaustive enumeration, or explicit argument for completeness of this set is provided; any missed lower-energy structure would alter the reported stability ordering and the interpretation of the crystallization anomaly.

    Authors: We agree that our study relies on the set of polymorphs predicted in recent literature rather than performing a new global structure search. These candidates represent the low-energy structures identified by prior computational work, and our calculations demonstrate that vdW forces stabilize B2O3-I while leaving multiple metastable structures near-degenerate with the glass. We will add an explicit discussion of this limitation and its implications for the crystallization anomaly interpretation. revision: partial

  2. Referee: [§4–5 (energy differences and convergence)] The key result is that vdW corrections reverse the ordering by only a few kcal/mol. However, the manuscript does not report basis-set convergence tests, k-point sampling details, or statistical error bars from the QMC runs that would establish whether the reported differences are robust to the level of numerical precision required for such small energy scales.

    Authors: The referee is correct that the original manuscript omitted detailed convergence information. Although the calculations used converged settings (as indicated by consistency across DFT, RPA, and QMC), we will add basis-set convergence tests, k-point sampling details, and QMC statistical error bars (with the latter already computed but not reported) to the revised manuscript, preferably in the supplementary information, to confirm the robustness of the few kcal/mol differences. revision: yes

standing simulated objections not resolved
  • A new global structure search to exhaustively confirm the completeness of the polymorph set is beyond the scope of the current revision and would require substantial additional computational effort.

Circularity Check

0 steps flagged

No circularity: direct first-principles energy ordering

full rationale

The derivation computes cohesive energies of B2O3 polymorphs via DFT, RPA and QMC on a fixed set of input structures. The claim that vdW forces make B2O3-I lowest in energy follows from explicit numerical differences between those computed energies; no parameter is fitted to the target ordering, no result is defined in terms of itself, and no uniqueness theorem or ansatz is imported via self-citation to force the conclusion. The completeness of the structure list is an external modeling choice, not a definitional reduction. The paper is therefore self-contained against its computational benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the accuracy of RPA and QMC for dispersion-dominated energy differences in boron-oxygen networks and on the completeness of the set of predicted polymorphs; no explicit free parameters are named in the abstract, but standard DFT dispersion corrections implicitly contain fitted coefficients.

axioms (1)
  • domain assumption RPA and QMC provide sufficiently accurate relative energies for B2O3 polymorphs differing by a few kcal/mol
    Invoked when the abstract states that these methods demonstrate the vdW role in stability ordering.

pith-pipeline@v0.9.0 · 5720 in / 1417 out tokens · 20440 ms · 2026-05-25T16:24:53.489618+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

71 extracted references · 71 canonical work pages · 1 internal anchor

  1. [1]

    Ferlat, Rings in Network Glasses: The B2O3 Case (Springer, Switzerland, 2015), chap

    G. Ferlat, Rings in Network Glasses: The B2O3 Case (Springer, Switzerland, 2015), chap. 14, p. 367

  2. [2]

    A. C. Wright, Phys. Chem. Glasses: Eur. J. Glass Sci. Technol. B59, 65 (2018)

  3. [3]

    Nicholas, S

    J. Nicholas, S. Sinogeikin, J. Kieffer, and J. Bass, Phys. Rev. Lett. 92, 215701 (2004)

  4. [4]

    S. K. Lee, K. Mibe, Y. Fei, G. D. Cody, and B. O. Mysen, Phys. Rev. Lett.94, 165507 (2005)

  5. [5]

    Trachenko, V

    K. Trachenko, V. V. Brazhkin, G. Ferlat, M. T. Dove, and E. Artacho, Phys. Rev. B78, 172102 (2008)

  6. [6]

    Zeidler, K

    A. Zeidler, K. Wezka, D. A. J. Whittaker, S. P. Salmon, A. Baroni, S. Klotz, H. E. Fischer, M. C. Wilding, C. L. Bull, M. G. Tucker, et al., Phys. Rev. B90, 024206 (2014)

  7. [7]

    S. K. Lee, Y.-H. Kim, P. Chow, Y. Xiao, J. Cheng, and G. Shen, Proc. Natl. Acad. Sci. USA115, 5855 (2018)

  8. [8]

    V. V. Brazhkin, I. Farnan, K. I. Funakoshi, M. Kanzaki, Y. Katayama, A. G. Lyapin, and H. Saitoh, Phys. Rev. Lett. 105, 115701 (2010)

  9. [9]

    O. L. G. Alderman, G. Ferlat, A. Baroni, M. Salanne, M. Micoulaut, C. J. Benmore, A. Lin, A. Tamalonis, and J. K. R. Weber, J. Phys. Condens. Matter27, 455104 (2015)

  10. [10]

    G. E. Gurr, P. W. Montgomery, C. D. Knutson, and B. T. Gorres, Acta Cristallogr.B26, 906 (1970)

  11. [11]

    S. S. Cole and N. W. Taylor, J. Am. Ceram. Soc.18, 55 (1935)

  12. [12]

    P. M. Piccione, C. Laberty, S. Yang, M. A. Camblor, A. Navrotsky, and M. E. Davis, J. Phys. Chem. B104, 10001 (2000)

  13. [13]

    R. E. Youngman, S. T. Haubrich, J. W. Zwanziger, M. T. Janicke, and B. F. Chmelka, Science269, 1416 (1995)

  14. [14]

    Umari and A

    P. Umari and A. Pasquarello, Phys. Rev. Lett.95, 137401 (2005)

  15. [15]

    E. D. Zanotto and D. R. Cassar, Sci. Rep.7, 43022 (2017)

  16. [16]

    D. R. Ulhmann, J. F. Hays, and D. Turnbull, Phys. Chem. Glasses 8, 1 (1967)

  17. [17]

    Takada, C

    A. Takada, C. R. A. Catlow, and G. D. Price, Phys. Chem. Glasses 44, 147 (2003)

  18. [18]

    Huang, M

    L. Huang, M. Durandurdu, and J. Kieffer, J. Phys. Chem. C 111, 13712 (2007)

  19. [19]

    Ferlat, A

    G. Ferlat, A. P. Seitsonen, M. Lazzeri, and F. Mauri, Nature Mat. 11, 925 (2012)

  20. [20]

    Claeyssens, J

    F. Claeyssens, J. N. Hart, N. C. Norman, and N. L. Allan, Adv. Funct. Mater.23, 5887 (2013)

  21. [21]

    Daub and H

    M. Daub and H. Hillebrecht, Eur. J. Inorg. Chem.2015, 4176 (2015)

  22. [22]

    Sorella, M

    S. Sorella, M. Casula, and D. Rocca, J. Chem. Phys.127, 014105 (2007)

  23. [23]

    C. J. Umrigar, J. Toulouse, C. Filippi, S. Sorella, and R. G. Hennig, Phys. Rev. Lett.98, 110201 (2007)

  24. [24]

    Barborini, S

    M. Barborini, S. Sorella, and L. Guidoni, J. Chem. Theory Comput. 8, 1260 (2012)

  25. [25]

    Casula, C

    M. Casula, C. Filippi, and S. Sorella, Phys. Rev. Lett.95, 100201 (2005)

  26. [26]

    Casula, S

    M. Casula, S. Moroni, S. Sorella, and C. Filippi, J. Chem. Phys. 132, 154113 (2010)

  27. [27]

    Dagrada, S

    M. Dagrada, S. Karakuzu, V. L. Vildosola, M. Casula, and S. Sorella, Phys. Rev. B94, 245108 (2016)

  28. [28]

    Sorella, Turborvb, quantum monte carlo software for electronic structure calculations , http://people.sissa.it/ sorella/web/

    S. Sorella, Turborvb, quantum monte carlo software for electronic structure calculations , http://people.sissa.it/ sorella/web/

  29. [29]

    D. C. Langreth and J. P. Perdew, Solid State Comm.17, 1425 (1975)

  30. [30]

    Gunnarsson and B

    O. Gunnarsson and B. I. Lundqvist, Phys. Rev. B13, 4274 (1976)

  31. [31]

    Marini, P

    A. Marini, P. García-González, and A. Rubio, Phys. Rev. Lett. 96, 136404 (2006)

  32. [32]

    Lebègue, J

    S. Lebègue, J. Harl, T. Gould, J. G. Ángyán, G. Kresse, and J. F. Dobson, Phys. Rev. Lett.105, 196401 (2010)

  33. [33]

    Kresse and J

    G. Kresse and J. Furthmüller, Phys. Rev. B54, 11169 (1996)

  34. [34]

    Kresse and D

    G. Kresse and D. Joubert, Phys. Rev. B59, 1758 (1999)

  35. [35]

    Grimme, J

    S. Grimme, J. Comput. Chem.27, 1787 (2006)

  36. [36]

    Tkatchenko and M

    A. Tkatchenko and M. Scheffler, Phys. Rev. Lett.102, 073005 (2009)

  37. [37]

    K. Lee, E. D. Murray, L. Kong, B. I. Lundqvist, and D. C. Langreth, Phys. Rev. B82, 081101(R) (2010)

  38. [38]

    Thonhauser, S

    T. Thonhauser, S. Zuluaga, C. A. Arter, K. Berland, E. Schröder, and P. Hyldgaard, Phys. Rev. Lett.115, 136402 (2015)

  39. [39]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996)

  40. [40]

    S. J. Clark, M. D. Segall, C. J. Pickard, P. J. Hasnip, M. J. Probert, K. Refson, and M. C. Payne, Z. Kristallogr. 220, 567 (2005)

  41. [41]

    Giannozzi, S

    P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococcioni, I. Dabo, et al., J. Phys. Condens. Matter 21, 395502 (2009), URL http://www.quantum-espresso.org

  42. [42]

    Coudert, Phys

    F.-X. Coudert, Phys. Chem. Chem. Phys. 15, 16012 (2013). 7

  43. [43]

    Hay, Ph.D

    H. Hay, Ph.D. thesis, Université Pierre & Marie Curie, Paris (2016), URL https://tel.archives-ouvertes. fr/tel-0147013a1

  44. [44]

    P. Haas, F. Tran, and P. Blaha, Phys. Rev. B79, 085104 (2009)

  45. [45]

    H. Hay, G. Ferlat, M. Casula, A. P. Seitsonen, and F. Mauri, Phys. Rev. B92, 144111 (2015)

  46. [46]

    W. Sun, S. T. Dacek, S. P. Ong, G. Hautier, A. Jain, W. D. Richards, A. C. Gamst, K. A. Persson, and G. Ceder, Sci. Adv. 2, e1600225 (2016)

  47. [47]

    7) is read- ily available in the CRYSTAL14 code [58] for D2 only (among vdW schemes)

    The calculation of the mechanical moduli (Fig. 7) is read- ily available in the CRYSTAL14 code [58] for D2 only (among vdW schemes)

  48. [48]

    10.8 kcal/(mol B2O3)

    For binary oxides, the range of observed polymorphs (as defined by the 90th percentile of a statistical analysis [46]) is 94 meV/atom, i.e. 10.8 kcal/(mol B2O3). Taking silica as a close parent system, the highest energy above quartz (among SiO2 polymorphs for which calorimetric data are available) is 6.88 kcal/(mol 2SiO2) (ISV zeolite) [12]

  49. [49]

    Wang, G.-C

    M.-S. Wang, G.-C. Guo, W.-T. Chen, G. Xu, W.-W. Zhou, K.-J. Wu, and J.-S. Huang, Angew. Chem. Int. Ed.46, 3909 (2007)

  50. [50]

    M.-C. Liu, P. Zhou, H. G. Yao, S.-H. Ji, R.-C. Zhang, M. Ji, and Y.-L. An, Eur. J. Inorg. Chem.31, 4622 (2009)

  51. [51]

    N. E. Shmidt, Russ. J. Inorg. Chem.11, 241 (1966)

  52. [52]

    However, the criterion on anisotropy likely does not apply to layered structures (T0, T0-0.5b and T0-b)

    Minimal value of the Young modulusEmin≥ 20 GPa and elastic anisotropyη≤ 5, whereη is max (Emax Emin ,Gmax Gmin ) . However, the criterion on anisotropy likely does not apply to layered structures (T0, T0-0.5b and T0-b)

  53. [53]

    Perim, D

    E. Perim, D. Lee, Y. Liu, C. Toher, P. Gong, Y. Li, W. N. Simmons, O. Levy, J. J. Vlassak, J. Schroers, et al., Nat. Commun. 7, 12315 (2016)

  54. [54]

    C. H. L. Goodman, Nature257, 370 (1975)

  55. [55]

    G. G. Naumis, Phys. Rev. E85, 061505 (2012)

  56. [56]

    Ronceray and P

    P. Ronceray and P. Harrowell, Phys. Rev. E96, 042602 (2017)

  57. [57]

    Russo, F

    J. Russo, F. Romano, and H. Tanaka, Phys. Rev. X8, 021040 (2018)

  58. [58]

    Dovesi, R

    R. Dovesi, R. Orlando, A. Erba, C. M. Zicovich-Wilson, B. Civalleri, S. Casassa, L. Maschio, M. Ferrabone, M. De La Pierre, P. D’Arco, et al., Int. J. Quantum Chem. 114, 1287 (2014)

  59. [59]

    Winkler, C

    B. Winkler, C. J. Pickard, V. Milman, and G. Thimm, Chem. Phys. Lett.337, 36 (2001)

  60. [60]

    Vanderbilt, Phys

    D. Vanderbilt, Phys. Rev. B41, 7892 (1990)

  61. [61]

    Richet, Y

    P. Richet, Y. Bottinga, L. Denielou, J. Petitet, and C. Tequi, Geochim. Cosmochim. Acta46, 2639 (1982)

  62. [62]

    Mouhat and F.-X

    F. Mouhat and F.-X. Coudert, Phys. Rev. B90, 224104 (2014)

  63. [63]

    Gaillac, P

    R. Gaillac, P. Pullumbi, and F.-X. Coudert, J. Phys. Condens. Matter 28, 275201 (2016)

  64. [64]

    Burkatzki, C

    M. Burkatzki, C. Filippi, and M. Dolg, J. Chem. Phys. 126, 234105 (2007)

  65. [65]

    Casula, C

    M. Casula, C. Attaccalite, and S. Sorella, J. Chem. Phys. 121, 7110 (2004)

  66. [66]

    Marchi, S

    M. Marchi, S. Azadi, M. Casula, and S. Sorella, J. Chem. Phys. 131, 154116 (2009)

  67. [67]

    Devaux, M

    N. Devaux, M. Casula, F. Decremps, and S. Sorella, Phys. Rev. B 91, 081101 (2015)

  68. [68]

    Calandra Buonaura and S

    M. Calandra Buonaura and S. Sorella, Phys. Rev. B57, 11446 (1998)

  69. [69]

    Assaraf and M

    R. Assaraf and M. Caffarel, J. Chem. Phys.119, 10536 (2003)

  70. [70]

    Filippi and C

    C. Filippi and C. J. Umrigar, Phys. Rev. B61, R16291 (2000)

  71. [71]

    van der Waals forces stabilize low-energy polymorphism in B2O3: Implications for the crystallization anomaly

    S. Sorella and L. Capriotti, J. Chem. Phys.133, 234111 (2010). 1 Supporting Material van der Waals forces stabilize low-energy polymorphism in B 2O3: implications for the crystallization anomaly Guillaume Ferlat1, Maria Hellgren1, François-Xavier Coudert2, Henri Hay1, Francesco Mauri1,3, Michele Casula1 1Sorbonne Université, MNHN, UMR CNRS 7590, IRD, IMPM...