Bayesian Inference for NMR Spectroscopy with Applications to Chemical Quantification

Nuclear magnetic resonance (NMR) spectroscopy exploits the magnetic properties of atomic nuclei to discover the structure, reaction state and chemical environment of molecules. We propose a probabilistic generative model and inference procedures for NMR spectroscopy. Specically, we use a weighted sum of trigonometric functions undergoing exponential decay to model free induction decay (FID) signals. We discuss the challenges in estimating the components of this general model { amplitudes, phase shifts, frequencies, decay rates, and noise variances { and oer practical solutions. We compare with conventional Fourier transform spectroscopy for estimating the relative concentrations of chemicals in a mixture, using synthetic and experimentally acquired FID signals. We nd the proposed model is particularly robust to low signal to noise ratios (SNR), and overlapping peaks in the Fourier transform of the FID, enabling accurate predictions (e.g., 1% sensitivity at low SNR) which are not possible with conventional spectroscopy (5% sensitivity). 1. Introduction. Nuclear magnetic resonance (NMR) spectroscopy has greatly advanced our understanding of molecular properties, and is now widespread in analytical chemistry. The theory of nuclear magnetic resonance postulates that protons and neutrons behave like gyroscopes that spin about their axes, generating their own small magnetic elds. These concepts were rst described by Rabi et al. (1939), for which Isidor Rabi was awarded the 1944 Nobel prize in physics. Later, Bloch, Hanson and Packard (1946) and Purcell, Torrey and Pound (1946) showed how NMR could be used to understand the structure of molecules in liquids and solids, for which they shared the Nobel prize in physics in 1952. Richard Ernst then won the 1991 Nobel prize in chemistry for developing Fourier transform NMR spectroscopy 1 , which led to the prevalence of NMR as an analytic technique. NMR spectroscopy is well suited to studying both organic and inorganic molecules, including proteins, and other biochemical species (Barrett et al., 2013), and is routinely used to identify the structure of unknown chemical species or the composition of mixtures. NMR spectroscopy is quantitative, chemically specic and non-invasive and therefore

[1]  M. Maiwald,et al.  Process and reaction monitoring by low-field NMR spectroscopy. , 2012, Progress in nuclear magnetic resonance spectroscopy.

[2]  W. W. Hansen,et al.  The Nuclear Induction Experiment , 1946 .

[3]  José Domingo Salazar,et al.  Application of a Bayesian deconvolution approach for high-resolution (1)H NMR spectra to assessing the metabolic effects of acute phenobarbital exposure in liver tissue. , 2010, Analytical chemistry.

[4]  Andrew C. Miklos,et al.  Protein nuclear magnetic resonance under physiological conditions. , 2009, Biochemistry.

[5]  H. Prosper Bayesian Analysis , 2000, hep-ph/0006356.

[6]  Carl E. Rasmussen,et al.  Occam's Razor , 2000, NIPS.

[7]  M Andrec,et al.  A Metropolis Monte Carlo implementation of bayesian time-domain parameter estimation: application to coupling constant estimation from antiphase multiplets. , 1998, Journal of magnetic resonance.

[8]  Yu-Chung N. Cheng,et al.  Magnetic Resonance Imaging: Physical Principles and Sequence Design , 1999 .

[9]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[10]  G. L. Bretthorst Bayesian analysis. I. Parameter estimation using quadrature NMR models , 1990 .

[11]  E. Purcell,et al.  Resonance Absorption by Nuclear Magnetic Moments in a Solid , 1946 .

[12]  Richard R. Ernst,et al.  Nuclear magnetic resonance Fourier transform spectroscopy , 1992 .

[13]  William J. Astle,et al.  A Bayesian Model of NMR Spectra for the Deconvolution and Quantification of Metabolites in Complex Biological Mixtures , 2011, 1105.2204.

[14]  J. Keeler Understanding NMR Spectroscopy , 2005 .

[15]  Charles R Sanders,et al.  The quiet renaissance of protein nuclear magnetic resonance. , 2013, Biochemistry.

[16]  W. W. Hansen,et al.  Nuclear Induction , 2011 .

[17]  Lynn F. Gladden,et al.  Nuclear magnetic resonance in chemical engineering: Principles and applications , 1994 .

[18]  S. Godsill,et al.  Deterministic and statistical methods for reconstructing multidimensional NMR spectra , 2006, Magnetic resonance in chemistry : MRC.

[19]  I. I. Rabi,et al.  The Molecular Beam Resonance Method for Measuring Nuclear Magnetic Moments The Magnetic Moments of 3 Li 6 , 3 Li 7 and 9 F 19 , 1939 .

[20]  G Larry Bretthorst,et al.  High dynamic‐range magnetic resonance spectroscopy (MRS) time‐domain signal analysis , 2009, Magnetic resonance in medicine.

[21]  Bruce J Balcom,et al.  The internal magnetic field distribution, and single exponential magnetic resonance free induction decay, in rocks. , 2005, Journal of magnetic resonance.

[22]  J. Griffin,et al.  Time-domain Bayesian detection and estimation of noisy damped sinusoidal signals applied to NMR spectroscopy. , 2007, Journal of magnetic resonance.

[23]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[24]  Scott L. Whittenburg,et al.  Bayesian Estimation of NMR Spectral Parameters Under Low Signal-to-Noise Conditions , 1993 .

[25]  L. Dou,et al.  Bayesian inference and Gibbs sampling in spectral analysis and parameter estimation: II , 1995 .

[26]  M. F. Cardoso,et al.  The simplex-simulated annealing approach to continuous non-linear optimization , 1996 .

[27]  F. Malz,et al.  Validation of quantitative NMR. , 2005, Journal of pharmaceutical and biomedical analysis.

[28]  Elias Aboutanios,et al.  Instantaneous frequency based spectral analysis of nuclear magnetic resonance spectroscopy data , 2012, Comput. Electr. Eng..