An official website of the United States government
We present an approach to approximating static properties of glasses without experimental inputs rooted in the first-principles random structure sampling. In our approach, the glassy system is represented by a collection (composite) of periodic, small-cell (few 10 s of atoms) local minima on the potential energy surface. These are obtained by generating a set of periodic structures with random lattice parameters and random atomic positions, which are then relaxed to their closest local minima on the potential energy surface using the first-principles methods. Using vitreous SiO2 as an example, we illustrate and discuss how well various atomic and electronic structure properties calculated as averages over the set of such local minima reproduce experimental data. The practical benefit of our approach, which can be rigorously thought of as representing an infinitely quickly quenched liquid, is in that it transfers the computational burden to linear scaling and easy to converge averages of properties computed on small-cell structures, rather than simulation cells with 100 s if not 1000 s of atoms while retaining a good overall predictive accuracy. Because of this, it enables the future use of high-cost/high-accuracy electronic structure methods, thereby bringing the modeling of glasses and amorphous phases closer to the state of modeling of crystalline solids.
more »
« less
The glassy solid as a statistical ensemble of crystalline microstates
Abstract We present an alternative and, for the purpose of non-crystalline materials design, a more suitable description of covalent and ionic glassy solids as statistical ensembles of crystalline local minima on the potential energy surface. Motivated by the concept of partially broken ergodicity, we analytically formulate the set of approximations under which the structural features of ergodic systems such as the radial distribution function (RDF) and powder X-ray diffraction (XRD) intensity can be rigorously expressed as statistical ensemble averages over different local minima. Validation is carried out by evaluating these ensemble averages for elemental Si and SiO2over the local minima obtained through the first-principles random structure sampling that we performed using relatively small simulation cells, thereby restricting the sampling to a set of predominantly crystalline structures. The comparison with XRD and RDF from experiments (amorphous silicon) and molecular dynamics simulations (glassy SiO2) shows excellent agreement, thus supporting the ensemble picture of glasses and opening the door to fully predictive description without the need for experimental inputs.
more »
« less
- Award ID(s):
- 1945010
- PAR ID:
- 10150391
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- npj Computational Materials
- Volume:
- 6
- Issue:
- 1
- ISSN:
- 2057-3960
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Efforts to reduce the carbon footprint associated with cement and concrete production have resulted in a number of promising lower‐emission alternatives. Still, research has emphasized a small subset of potentially useful precursor materials. With the goal of expanding the precursor pool, this work presents results of parallel literature mining and rate modeling activities. As a result of literature mining, materials with appropriate SiO2, Al2O3, and CaO concentrations were assembled into a comprehensive, representative ternary diagram. 23 000+ materials were extracted from 7000 journal articles, and 7500 materials from 6000 articles with 80 ≤ SiO2 + Al2O3 + CaO ≤105 wt% automatically classified. Both supervised and semi‐supervised models were used for dissolution rate prediction of glassy materials with all models pulling from a single data set (n = 802 reported dissolution rates from 105 different glasses). Supervised modeling utilized linear and decision tree regressions to determine features most predictive of dissolution rate, resulting in log‐linear relationships between rate and pH, inverse temperature (1/K), and non‐bridging oxygen per tetrahedron (NBO/T). Semi‐supervised modeling was observed to be more robust to broader feature inclusion, providing similar predictive ability with a relatively larger set of descriptive features. Most importantly, results indicated that models trained on data from disparate scientific communities were adequately predictive (RMSE ≈ 1), particularly under pH ≥7 conditions relevant to the cement and alkali activation communities.more » « less
-
Abstract Solid‐state NMR measurements coupled with density functional theory (DFT) calculations demonstrate how hydrogen positions can be refined in a crystalline system. The precision afforded by rotational‐echo double‐resonance (REDOR) NMR to interrogate13C–1H distances is exploited along with DFT determinations of the13C tensor of carbonates (CO32−). Nearby1H nuclei perturb the axial symmetry of the carbonate sites in the hydrated carbonate mineral, hydromagnesite [4 MgCO3⋅Mg(OH)2⋅4 H2O]. A match between the calculated structure and solid‐state NMR was found by testing multiple semi‐local and dispersion‐corrected DFT functionals and applying them to optimize atom positions, starting from X‐ray diffraction (XRD)‐determined atomic coordinates. This was validated by comparing calculated to experimental13C{1H} REDOR and13C chemical shift anisotropy (CSA) tensor values. The results show that the combination of solid‐state NMR, XRD, and DFT can improve structure refinement for hydrated materials.more » « less
-
Abstract We present a new method to construct an internally consistent thermodynamic model using a compilation of high‐pressure melting experiments. The steepest descent method and Monte Carlo sampling are combined to constrain all model parameters simultaneously instead of determining each parameter sequentially from relevant experiments. Our approach is applied to the published melting experiments on mantle materials to obtain the thermodynamic parameters of the MgO‐FeO‐SiO2ternary system. Inversion with the subsets of experimental data is conducted as well to investigate the source of discrepancy among existing studies, and the key parameters are found to be the thermal expansivity of SiO2and the excess volume of mixing between MgO and SiO2. Mixing between FeO and SiO2is only constrained with large uncertainty, which could also imply that oxides with low concentrations have minimal effects on melting. Constraining the thermodynamics of MgO and SiO2will be important for a better understanding of mantle melting at high pressures.more » « less
-
Abstract In this study, high yields of CO are reported from CO2using the silica (SiO2) supported perovskite oxide, La0.75Sr0.25FeO3(LSF), composites in the reverse water gas shift chemical looping (RWGS‐CL) process. XRD patterns of materials formed upon adding SBA‐15 to the perovskite sol‐gel precursor solution indicated successful formation of an orthorhombic perovskite oxide structure in the composites. The total surface area increased by ∼300 % with the addition of 50 % LSF to SBA‐15 by mass and surface accessibility of perovskite oxide crystallites was verified by CO2chemisorption and XPS measurements. Composite materials achieved up to a factor of 10 increases in CO yields (∼3.5 vs 0.35 mmol CO/gLSF) compared to pure LSF through six consecutive RWGS‐CL cycles at 700 °C. Following these RWGS‐CL cycles, XRD Scherrer analyses showed that the perovskite oxide in the composite material decreased in crystallite size. This approach to synthesis of supported perovskite oxides is expected to be valuable for large‐scale CO2conversion by RWGS‐CL.more » « less
