Spectral estimation of NMR relaxation

In this paper, spectral estimation of NMR relaxation is constructed as an extension of Fourier Transform (FT) theory as it is practiced in NMR or MRI, where multidimensional FT theory is used. nD NMR strives to separate overlapping resonances, so the treatment given here deals primarily with monoexponential decay. In the domain of real error, it is shown how optimal estimation based on prior knowledge can be derived. Assuming small Gaussian error, the estimation variance and bias are derived. Minimum bias and minimum variance are shown to be contradictory experimental design objectives. The analytical continuation of spectral estimation is constructed in an optimal manner. An important property of spectral estimation is that it is phase invariant. Hence, hypercomplex data storage is unnecessary. It is shown that, under reasonable assumptions, spectral estimation is unbiased in the context of complex error and its variance is reduced because the modulus of the whole signal is used. Because of phase invariance, the labor of phasing and any error due to imperfect phase can be avoided. A comparison of spectral estimation with nonlinear least squares (NLS) estimation is made analytically and with numerical examples. Compared to conventional sampling for NLS estimation, spectral estimation would typically provide estimation values of comparable precision in one-quarter to one-tenth of the spectrometer time when S/N is high. When S/N is low, the time saved can be used for signal averaging at the sampled points to give better precision. NLS typically provides one estimate at a time, whereas spectral estimation is inherently parallel. The frequency dimensions of conventional nD FT NMR may be denoted D₁, D₂, etc. As an extension of nD FT NMR, one can view spectral estimation of NMR relaxation as an extension into the zeroth dimension. In nD NMR, the information content of a spectrum can be extracted as a set of n-tuples (ω₁, … ω_n), corresponding to the peak maxima. Spectral estimation of NMR relaxation allows this information content to be extended to a set of (n + 1)-tuples (λ, ω₁, … ω_n), where λ is the relaxation rate.