Title: Machine learning of 27Al NMR electric field gradient tensors for crystalline structures from DFT
NMR crystallography has emerged as a promising technique for the determination and refinement of atomic coordinates in crystal structures. The crystal structure of compounds containing quadrupolar nuclei, such as 27Al, can be improved by directly comparing solid-state NMR measurements to DFT computations of the electric field gradient (EFG) tensor. The non-negligible computational cost of these first-principles calculations limits the applicability of this method to all but the most well-defined structures. We developed a fast, low-cost machine learning model to predict EFG parameters based on local structural motifs and elemental parameters. We computed 8081 EFG tensors from 1681 27Al crystalline solids using DFT and benchmarked them against 105 experimentally measured 27Al sites. Surprisingly, simple local geometric features dominate the predictive performance of the resulting random-forest model, yielding an R2 value of 0.98 and an RMSE of 0.61 MHz for CQ, the quadrupolar coupling constant. This model accuracy should enable pre-refining future structural assignments before finally validating with first-principles calculations. Such a catalogue of 27Al NMR tensors can serve as a tool for researchers assigning complex NMR spectra influenced by the nuclear electric quadrupole interaction.  more » « less
Award ID(s):
2004915 1640899
PAR ID:
10629814
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Scientific Reports (Nature)
Date Published:
Journal Name:
Scientific Reports
Volume:
15
ISSN:
2045-2322
Page Range / eLocation ID:
26456
Subject(s) / Keyword(s):
27Al NMR NMR crystallography electric field gradient tensor density functional theory machine learning
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT This work implements a methodology for studying quadrupolar nuclear spin relaxation in ionic liquids. The dynamic behavior of the ions in the liquid is described by ab initio molecular dynamics (aiMD) with forces obtained from density functional theory (DFT) calculations with periodic boundary conditions and a non‐hybrid functional. The electric field gradient (EFG) driving the quadrupolar relaxation was calculated with free boundary conditions, using clusters that contained the ion of interest surrounded by two coordination shells treated quantum mechanically and augmented with a solvation model. Tests showed that EFG calculations using only the first coordination shell, containing five nearest neighbors, also provide a suitable model, because the relaxation rates differ by no more than 4% from the results from the two‐shell solvation. The results of this study show that the relaxation of the deuterated ethylammonium nitrate () occurs within the extreme narrowing regime for a spectrometer magnetic field T and is therefore characterized by the ensemble variance of the EFG and the correlation time associated with the EFG autocorrelation function. The quadrupolar relaxation of demanded molecular dynamics production times longer than 330 ps and averaging over multiple ions, as well as independent trajectories to get suitably converged relaxation rates. The calculated relaxation rate is Hz, about 60% above the rate reported experimentally. However, the approach utilized in the present study has an accuracy similar to, or better than, what has been previously reported for systems involving non‐ionic solvents that required simulations of 100 ps duration or less. 
    more » « less
  2. null (Ed.)
    Rechargeable aluminum–graphite batteries are an emerging energy storage technology with great promise: they exhibit high rate performance, cyclability, and a discharge potential of approximately 2 V, while both electrodes are globally abundant, low cost, and inherently safe. The batteries use chloroaluminate-containing electrolytes and store charge in the graphite electrodes when molecular AlCl 4 − anions electrochemically intercalate within them. However, much remains to be understood regarding the ion intercalation mechanism, in part due to challenges associated with characterizing the chloroaluminate anions themselves. Here, we use solid-state 27 Al nuclear magnetic resonance (NMR) spectroscopy to probe the molecular-level electronic and magnetic environments of intercalated chloroaluminate anions at different states-of-charge. The results reveal broad 27 Al NMR signals associated with intercalated AlCl 4 − anions, reflecting high extents of local disorder. The intercalated anions experience a diversity of local environments, many of which are far from the ideal crystalline-like structures often depicted in graphite staging models. Density functional theory (DFT) calculations of the total 27 Al isotropic shifts enable the contributions of chemical shift, ring-current effects, and electric quadrupolar interactions to be disentangled quantitatively. In combination, the solid-state NMR and DFT results reveal the molecular geometries and environments of intercalated AlCl 4 − anions and capture the significant disorder present in intercalated graphite battery electrodes. 
    more » « less
  3. Abstract The deviation of the electron density around the nuclei from spherical symmetry determines the electric field gradient (EFG), which can be measured by various types of spectroscopy. Nuclear Quadrupole Resonance (NQR) is particularly sensitive to the EFG. The EFGs, and by implication NQR frequencies, vary dramatically across materials. Consequently, searching for NQR spectral lines in previously uninvestigated materials represents a major challenge. Calculated EFGs can significantly aid at the search’s inception. To facilitate this task, we have applied high-throughput density functional theory calculations to predict EFGs for 15187 materials in the JARVIS-DFT database. This database, which will include EFG as a standard entry, is continuously increasing. Given the large scope of the database, it is impractical to verify each calculation. However, we assess accuracy by singling out cases for which reliable experimental information is readily available and compare them to the calculations. We further present a statistical analysis of the results. The database and tools associated with our work are made publicly available by JARVIS-DFT ( https://www.ctcms.nist.gov/~knc6/JVASP.html ) and NIST-JARVIS API ( http://jarvis.nist.gov/ ). 
    more » « less
  4. We present a first-principles study based on density functional theory (DFT) on the electronic and structural properties of Ta2⁢NiSe5, a layered transition metal chalcogenide that has been considered as a possible candidate for an excitonic insulator. Our systematic DFT results however provide a nonexcitonic mechanism for the experimentally observed electronic and structural phase transitions in Ta2⁢NiSe5, in particular explaining why sulfur substitution of selenium reduces the distortion angle in the low-temperature phase and potassium dosing closes the gap in the electronic structure. Moreover, the calculations show that these two effects couple to each other. Further, our first-principles calculations predict several changes in both the crystal structure and electronic structure under the effects of uniform charge dosing and uniaxial strain, which could be tested experimentally. 
    more » « less
  5. This paper reports the derivation and implementation of the electric dipole-magnetic dipole and electric dipole-electric quadrupole polarizability tensors at the density functional theory level with periodic boundary conditions (DFT-PBC). These tensors are combined to evaluate the Buckingham/Dunn tensor that describes the optical rotation (OR) in oriented chiral systems. We describe several aspects of the derivation of the equations and present test calculations that verify the correctness of the tensor formulation and their implementation. The results show that the full OR tensor is completely origin invariant as for molecules and that PBC calculations match molecular cluster calculations on 1D chains. A preliminary investigation on the choice of density functional, basis set, and gauge indicates a similar dependence as for molecules: the functional is the primary factor that determines the OR magnitude, followed by the basis set and to a much smaller extent the choice of gauge. However, diffuse functions may be problematic for PBC calculations even if they are necessary for the molecular case. A comparison with experimental data of OR for the tartaric acid crystal shows reasonable agreement given the level of theory employed. The development presented in this paper offers the opportunity to simulate the OR of chiral crystalline materials with general-purpose DFT-PBC methods, which, in turn, may help to understand the role of intermolecular interactions on this sensitive electronic property. 
    more » « less