An official website of the United States government
The glass transition temperatures of common binary oxides, including those with low glass-forming ability, are estimated using pair distribution functions (PDFs) from ab initio molecular dynamics simulations. The computed glass transition temperatures for good glass-formers such as silica (SiO2), germania (GeO2), and boron oxide (B2O3) are in agreement with measured values. These calculations are then used to compute the glass transition temperatures of alumina (Al2O3), tantala (Ta2O5), and telluria (TeO2), which are known to exhibit low glass-forming ability. For Al2O3 and Ta2O5, we also compute the simulated caloric curve from molecular dynamics simulations using two-body empirical force fields. Finally, we discuss the possibility of extracting the glass transition temperature by measuring the thermal broadening of the PDFs from scattering measurements.
more »
« less
How close are the classical two-body potentials to ab initio calculations? Insights from linear machine learning based force matching
In this work, we propose a linear machine learning force matching approach that can directly extract pair atomic interactions from ab initio calculations in amorphous structures. The local feature representation is specifically chosen to make the linear weights a force field as a force/potential function of the atom pair distance. Consequently, this set of functions is the closest representation of the ab initio forces, given the two-body approximation and finite scanning in the configurational space. We validate this approach in amorphous silica. Potentials in the new force field (consisting of tabulated Si–Si, Si–O, and O–O potentials) are significantly different than existing potentials that are commonly used for silica, even though all of them produce the tetrahedral network structure and roughly similar glass properties. This suggests that the commonly used classical force fields do not offer fundamentally accurate representations of the atomic interaction in silica. The new force field furthermore produces a lower glass transition temperature (Tg ∼ 1800 K) and a positive liquid thermal expansion coefficient, suggesting the extraordinarily high Tg and negative liquid thermal expansion of simulated silica could be artifacts of previously developed classical potentials. Overall, the proposed approach provides a fundamental yet intuitive way to evaluate two-body potentials against ab initio calculations, thereby offering an efficient way to guide the development of classical force fields.
more »
« less
- Award ID(s):
- 2309000
- PAR ID:
- 10503145
- Publisher / Repository:
- AIP Publishing
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 160
- Issue:
- 5
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Vibrational spectroscopy has been widely used to investigate various structural aspects of the glass network, and there are a plethora of papers reporting subtle but consistent changes in infrared and Raman spectral features of glass upon alterations of glass compositions, thermal histories, mechanical stresses, or surface treatments. However, interpretations of such spectral features are still obscured due to the lack of well‐established physical principles accurately describing vibrational modes of the non‐crystalline glass network. Due to the non‐equilibrium nature of the glass network, three‐dimensionally connected without any long‐range orders, vibrational spectral features of glass cannot be interpreted using the analogy with those of isolated molecular moieties or crystalline counterparts. This feature article explains why such comparisons are outdated and describes the recent advances made from theoretical calculations of vibrational spectral features of amorphous networks or comparisons of computational results with experimental data. For the interpretation of vibrational spectral features of silica and silicate glasses, the following empirical relationships are suggested: (i) the intensity‐weighted peak position of the Si‐O‐Si stretch mode negatively correlates with the weighted average of the Si‐O bond length distribution, and (ii) the broad band of the Si‐O‐Si bending mode negatively correlates with the Si‐O‐Si bond angle distribution. Selected examples of vibrational spectroscopic imaging of surface defects are discussed to deliberate the implication of these findings in the structure‐property relationship of silica and silicate glass materials. Unanswered questions and continuing research challenges are identified.more » « less
-
Ab initio design of polymer nanocomposite materials for high breakdown strength requires prediction of localized trap states at the polymer–filler interface. Systematic first-principles calculations of realistic interfaces can be challenging, particularly for amorphous polymers and fillers that necessitate the calculation of ensembles of large unit cells with hundreds of atoms. We present a computational approach for automatically generating reasonable structures for amorphous polymer–filler interfaces, combining classical molecular dynamics and Monte Carlo simulations. We identify trap states by analyzing the localization of electronic eigenstates calculated using density functional theory on ensembles of interface structures, clearly distinguishing shallow trap states from delocalized band-edge states. Applying this approach to silica–polyethylene interfaces as an initial example, we find under-coordination and distorted coordination structures at amorphous silica surfaces contribute a combination of deep and shallow traps at these interfaces, whereas polyethylene does not generate localized interfacial states.more » « less
-
Transferability evaluation of the deep potential model for simulating water-graphene confined systemMachine learning potentials (MLPs) are poised to combine the accuracy of ab initio predictions with the computational efficiency of classical molecular dynamics (MD) simulation. While great progress has been made over the last two decades in developing MLPs, there is still much to be done to evaluate their model transferability and facilitate their development. In this work, we construct two deep potential (DP) models for liquid water near graphene surfaces, Model S and Model F, with the latter having more training data. A concurrent learning algorithm (DP-GEN) is adopted to explore the configurational space beyond the scope of conventional ab initio MD simulation. By examining the performance of Model S, we find that an accurate prediction of atomic force does not imply an accurate prediction of system energy. The deviation from the relative atomic force alone is insufficient to assess the accuracy of the DP models. Based on the performance of Model F, we propose that the relative magnitude of the model deviation and the corresponding root-mean-square error of the original test dataset, including energy and atomic force, can serve as an indicator for evaluating the accuracy of the model prediction for a given structure, which is particularly applicable for large systems where density functional theory calculations are infeasible. In addition to the prediction accuracy of the model described above, we also briefly discuss simulation stability and its relationship to the former. Both are important aspects in assessing the transferability of the MLP model.more » « less
-
Ab initio molecular dynamics liquid-quench simulations and hybrid density functional calculations are performed to model the effects of room-temperature atomic fluctuations and photo-illumination on the structural and electronic properties of amorphous sub-stoichiometric In2O2.96. A large configurational ensemble is employed to reliably predict the distribution of localized defects as well as their response to the thermal and light activation. The results reveal that the illumination effects on the carrier concentration are greater in amorphous configurations with shorter In–O bond length and reduced polyhedral sharing as compared to the structures with a more uniform morphology. The obtained correlation between the photo-induced carrier density and the reduction in the number of fully coordinated In-atoms implies that metal oxides with a significant fraction of crystalline/amorphous interfaces would show a more pronounced response to illumination. Photo-excitation also produces In–O2–In defects that have not been previously found in sub-stoichiometric amorphous oxides; these defects are responsible for carrier instabilities due to overdoping.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0175756
