On the Boroxol Ring Fraction in Melt-Quenched B${}_2$O${}_3$ Glass

Abstract

An atomistic structural model for melt-quenched B${}_2$O${}_3$ glass has eluded the simulation community so far. The difficulty lies in the abundance of the six-membered boroxol rings - an intermediate-range order motif suggested through Raman and NMR spectroscopy - which is challenging to obtain in atomistic molecular dynamics simulations. Here, we report the development of a DFT-accurate machine-learned potential for B${}_2$O${}_3$ and employ quench rates as low as 10${}^9$ K/s to obtain B${}_2$O${}_3$ glasses with more than 30% of boron atoms in boroxol rings. Also, we show that the pressure, and consequently the boroxol fraction, in the deep potential molecular dynamics (DPMD) simulations critically depends on the range of the geometry descriptor used in the embedding neural network, and at least a 9 $\unicode{x212B}$ range is required. The boroxol ring fraction increases with decreasing quench rate. Finally, amorphous B${}_2$O${}_3$ configurations display a minimum in energy at a boroxol fraction of 75%, intriguingly close to the experimental estimate in B${}_2$O${}_3$ glass.

References in corpus (6)

• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• DeePMD-kit v2: A software package for Deep Potential models
• Hidden polymorphs drive the vitrification in B2O3
• Understanding the Anomalous Diffusion of Water in Aqueous Electrolytes Using Machine Learned Potentials
• Monte Carlo simulations of glass-forming liquids beyond Metropolis
• Slowly Quenched, High Pressure Glassy B${}_2$O${}_3$ at DFT Accuracy