The open-source Python package, MRSimulator, is presented as a simple-to-use, fast, versatile, and extendable package capable of simulating one- and higher-dimensional Nuclear Magnetic Resonance (NMR) spectra under static, magic-angle, and variable-angle conditions. High benchmarks in spectral simulations are achieved by assuming that there are no degeneracies in the energy eigenstates, i.e., all dipolar couplings are in the weak limit and that there are no rotational resonances during evolution periods. Under these assumptions, the symmetry pathway formalism is exploited to reduce an NMR method applied to a spin system into a sum of individual transition pathways, whose signals are more efficiently calculated individually than as part of a full-density matrix simulation. To increase numerical efficiencies further, our approach restricts coherence transfer among transitions to pure rotations about an axis in the x–y plane of the rotating frame or through an artificial total mixing operation between selected transitions of adjacent free evolution periods. The assumptions used in this approach are valid for most commonly used solid-state NMR methods. Details of the implementation, along with example code usage, are given, including a least-squares spectral analysis.
more »
« less
Simulating multipulse NMR spectra of polycrystalline solids in the frequency domain
An approach is presented for simulating multipulse nuclear magnetic resonance (NMR) spectra of polycrystalline solids directly in the frequency domain. The approach integrates the symmetry pathway concept for multipulse NMR with efficient algorithms for calculating spinning sideband amplitudes and performing interpolated finite-element numerical integration over all crystallite orientations in a polycrystalline sample. The numerical efficiency is achieved through a set of assumptions used to approximate the evolution of a sparse density matrix through a pulse sequence as a set of individual transition pathway signals. The utility of this approach for simulating the spectra of complex materials, such as glasses and other structurally disordered materials, is demonstrated.
more »
« less
- Award ID(s):
- 2107636
- PAR ID:
- 10519288
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 160
- Issue:
- 23
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Nuclear magnetic resonance (NMR) spectroscopy is a powerful tool for obtaining precise information about the local bonding of materials, but difficult to interpret without a well-vetted dataset of reference spectra. The ability to predict NMR parameters and connect them to three-dimensional local environments is critical for understanding more complex, long-range interactions. New computational methods have revealed structural information available from29Si solid-state NMR by generating computed reference spectra for solids. Such predictions are useful for the identification of new silicon-containing compounds, and serve as a starting point for determination of the local environments present in amorphous structures. In this study, we have used 42 silicon sites as a benchmarking set to compare experimentally reported29Si solid-state NMR spectra with those computed by CASTEP-NMR and Vienna Ab Initio Simulation Program (VASP). Data-driven approaches enable us to identify the source of discrepancies across a range of experimental and computational results. The information from NMR (in the form of an NMR tensor) has been validated, and in some cases corrected, in an effort to catalog these for the local spectroscopy database infrastructure (LSDI), where over 10,00029Si NMR tensors for crystalline materials have been computed. Knowledge of specific tensor values can serve as the basis for executing NMR experiments with precision, optimizing conditions to capture the elements accurately. The ability to predict and compare experimental observables from a wide range of structures can aid researchers in their chemical assignments and structure determination, since the computed values enables the extension beyond tables of typical chemical shift (or shielding) ranges.more » « less
-
Rapid and automated lipid profiling by nuclear magnetic resonance spectroscopy using neural networksAbstract Nuclear magnetic resonance (NMR) spectroscopy is a powerful tool for quantitative metabolomics; however, quantification of metabolites from NMR data is often a slow and tedious process requiring user input and expertise. In this study, we propose a neural network approach for rapid, automated lipid identification and quantification from NMR data. Multilayered perceptron (MLP) networks were developed with NMR spectra as the input and lipid concentrations as output. Three large synthetic datasets were generated, each with 55,000 spectra from an original 30 scans of reference standards, by using linear combinations of standards and simulating experimental‐like modifications (line broadening, noise, peak shifts, baseline shifts) and common interference signals (water, tetramethylsilane, extraction solvent), and were used to train MLPs for robust prediction of lipid concentrations. The performances of MLPS were first validated on various synthetic datasets to assess the effect of incorporating different modifications on their accuracy. The MLPs were then evaluated on experimentally acquired data from complex lipid mixtures. The MLP‐derived lipid concentrations showed high correlations and slopes close to unity for most of the quantified lipid metabolites in experimental mixtures compared with ground‐truth concentrations. The most accurate, robust MLP was used to profile lipids in lipophilic hepatic extracts from a rat metabolomics study. The MLP lipid results analyzed by two‐way ANOVA for dietary and sex differences were similar to those obtained with a conventional NMR quantification method. In conclusion, this study demonstrates the potential and feasibility of a neural network approach for improving speed and automation in NMR lipid profiling and this approach can be easily tailored to other quantitative, targeted spectroscopic analyses in academia or industry.more » « less
-
Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.more » « less
-
We present a methodology for simulating multidimensional electronic spectra of molecular aggregates with coupling of electronic excitation to a structured environment using the stochastic non-Markovian quantum state diffusion (NMQSD) method in combination with perturbation theory for the response functions. A crucial aspect of our approach is that we propagate the NMQSD equation in a doubled system Hilbert space but with the same noise. We demonstrate that our approach shows fast convergence with respect to the number of stochastic trajectories, providing a promising technique for numerical calculation of two-dimensional electronic spectra of large molecular aggregates.more » « less
An official website of the United States government

