First-principles NMR of oxide glasses boosted by machine learning

Thibault Charpentier
DOI: https://doi.org/10.1039/d4fd00129j
2024-06-27
Faraday Discussions
Abstract:Solid-State NMR has established itself as a cutting-egde spectroscopy for elucidating the structure of oxide glasses thanks to several decades of methodological and instrumental progresses. First-principles calculations of NMR properties combined with molecular dynamics (MD) simulations provides a powerful complementing approach for the interpretation of the NMR data although they still suffer from limitations in terms of size, time and high consumption of computational resources. We address this challenge by developing a machine-learning framework to boost predictive modelling of NMR spectra. We use kernel ridge regression techniques (least-square support vector regression and linear ridge regression) combined with the smooth overlap of atomic position (SOAP) atom-centered descriptors to efficiently predict NMR interactions: isotropic magnetic shielding and the electric field gradient (EFG) tensor. As illustrated in this work, this approach enables the simulation of MAS and MQMAS NMR spectra of very large models (more than 10000 atoms) and an efficient averaging of NMR properties over MD trajectories of nanoseconds for incorporating finite temperature effects, at computational cost of classical MD simulation. We illustrate these advances on sodium silicate glasses (SiO2-Na2O). NMR parameters (isotropic chemical shift and electric field gradient) could be predicted with an accuracy of 1 to 2% in terms of the total span of the NMR parameter values. To include vibrational effects, an approach is proposed by scaling the EFG tensor in NMR simulations with a factor obtained from the time auto-correlation functions computed on MD trajectory.
chemistry, physical
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is: **How to use machine - learning methods to accelerate and improve the prediction of nuclear magnetic resonance (NMR) properties of oxide glasses, especially for the research of large - scale models and dynamic effects**. Specifically, the goals of the paper include the following aspects: 1. **Improve the efficiency of NMR property prediction** - By developing methods based on the machine - learning (ML) framework, especially kernel ridge regression (KRR) and linear ridge regression (LRR), combined with atom - centered descriptors (such as SOAP), to efficiently predict NMR interactions (such as isotropic magnetic shielding and electric field gradient tensors). - This method can significantly reduce the computational cost, making it possible to simulate MAS and MQMAS NMR spectra of large - scale models containing more than 10,000 atoms. 2. **Overcome the limitations of traditional first - principles calculations** - Although first - principles calculations (such as DFT - GIPAW) are accurate, they are limited by size, time, and computational resources when dealing with large - scale models and long - time dynamic trajectories. - The method proposed in the paper can effectively average NMR properties on the time scale of classical molecular dynamics (MD) simulations, thereby considering the finite - temperature effect. 3. **Study the influence of vibration effects on NMR spectra** - A method is proposed to simulate vibration effects by scaling the electric field gradient (EFG) tensor, which is based on the time - auto - correlation function calculated from MD trajectories. 4. **Verify the accuracy and applicability of the method** - The effectiveness of this method has been verified in sodium silicate glass (SiO₂ - Na₂O), proving that it can predict NMR parameters (such as isotropic chemical shift and electric field gradient) with an accuracy of 1 - 2%. - Comparing the prediction results with experimental data further verifies the reliability of the method. ### Summary The core goal of the paper is to develop an efficient machine - learning framework for predicting the NMR properties of oxide glasses while overcoming the size and time limitations of traditional computational methods. This provides a new tool for studying the relationship between glass structure and properties and opens up the possibility of designing new glass materials.