Trust-region algorithm for the inversion of molecular diffusion NMR data.

Diffusion NMR experiments are very useful in studying structural and dynamical properties of molecules and in sorting out components from mixtures. A number of numerical approaches have been developed for the processing of diffusion NMR data. In this paper, numerical problems of the direct regularization methods such as CONTIN, MaxEnt, and the newly proposed ITAMeD approach are illustrated by analyzing simulated and experimental data. It allows us to further develop a new method to calculate the distribution of diffusion coefficients. Therefore, we present here an iterative regularization method based on the Trust-Region Algorithm for the Inversion (TRAIn) of molecular diffusion NMR data. It is demonstrated in this paper that our approach overcomes major numerical difficulties of the direct regularization methods mentioned above. Besides, this method reconstructs more reliable diffusion coefficient distributions, especially for real world samples of which the diffusion coefficients are nonsymmetrically distributed.

[1]  M. Nilsson,et al.  Speedy component resolution: an improved tool for processing diffusion-ordered spectroscopy data. , 2008, Analytical chemistry.

[2]  Ya-Xiang Yuan,et al.  Convergence and regularity of trust region methods for nonlinear ill-posed inverse problems , 2005 .

[3]  R Huo,et al.  Improved DOSY NMR data processing by data enhancement and combination of multivariate curve resolution with non-linear least square fitting. , 2004, Journal of magnetic resonance.

[4]  Charles S. Johnson,et al.  Diffusion ordered 2D NMR spectroscopy of phospholipid vesicles: determination of vesicle size distributions , 1993 .

[5]  Ya-Xiang Yuan,et al.  On the superlinear convergence of a trust region algorithm for nonsmooth optimization , 1985, Math. Program..

[6]  J. E. Tanner.  Use of the Stimulated Echo in NMR Diffusion Studies , 1970 .

[7]  Van Gorkom LCM,et al.  Analysis of DOSY and GPC-NMR Experiments on Polymers by Multivariate Curve Resolution , 1998, Journal of magnetic resonance.

[8]  Charles S. Johnson,et al.  Determination of Molecular Weight Distributions for Polymers by Diffusion-Ordered NMR , 1995 .

[9]  Lutgarde M. C. Buydens,et al.  Assessment of techniques for DOSY NMR data processing , 2003 .

[10]  T. Steihaug.  The Conjugate Gradient Method and Trust Regions in Large Scale Optimization , 1983 .

[11]  C. Vogel.  A constrained least squares regularization method for nonlinear ill-posed problems , 1990 .

[12]  W. Windig,et al.  Direct exponential curve resolution algorithm (DECRA): A novel application of the generalized rank annihilation method for a single spectral mixture data set with exponentially decaying contribution profiles , 1997 .

[13]  T. Malliavin,et al.  Maximum Entropy Processing of DOSY NMR Spectra , 1998 .

[14]  Charles S. Johnson.  Diffusion Ordered Nuclear Magnetic Resonance Spectroscopy: Principles and Applications , 1999 .

[15]  Charles S. Johnson,et al.  Resolution of discrete and continuous molecular size distributions by means of diffusion-ordered 2D NMR spectroscopy , 1993 .

[16]  D. Bernin,et al.  Iterative thresholding algorithm for multiexponential decay applied to PGSE NMR data. , 2013, Analytical chemistry.

[17]  Shanmin Zhang.  Pivotal steps towards quantification of molecular diffusion coefficients by NMR. , 2007, Chemphyschem : a European journal of chemical physics and physical chemistry.

[18]  S. Provencher.  A constrained regularization method for inverting data represented by linear algebraic or integral equations , 1982 .

[19]  I. J. Day.  On the inversion of diffusion NMR data: Tikhonov regularization and optimal choice of the regularization parameter. , 2011, Journal of magnetic resonance.

[20]  J. Jakeš.  Testing of the constrained regularization method of inverting Laplace transform on simulated very wide quasielastic light scattering autocorrelation functions , 1988 .

[21]  M. J. D. Powell,et al.  On the global convergence of trust region algorithms for unconstrained minimization , 1984, Math. Program..

[22]  J. E. Tanner,et al.  Spin diffusion measurements : spin echoes in the presence of a time-dependent field gradient , 1965 .

[23]  S. Provencher,et al.  Regularization Techniques for Inverse Problems in Molecular Biology , 1983 .

[24]  S. Provencher.  CONTIN: A general purpose constrained regularization program for inverting noisy linear algebraic and integral equations , 1984 .

[25]  P. Griffiths,et al.  Global Least-Squares Analysis of Large, Correlated Spectral Data Sets: Application to Component-Resolved FT-PGSE NMR Spectroscopy , 1996 .

[26]  Qun Chen,et al.  A Simple Access to the (Log/Normal) Molecular Weight Distribution Parameters of Polymers Using PGSTE NMR , 2011 .