Commun. Comput. Phys., 15 (2014), pp. 1343-1351. A Frequency Determination Method for Digitized NMR Signals H. Yan 1*, K. Li 1, R. Khatiwada 1, E. Smith 1, W. M. Snow 1, C. B. Fu 2, P.-H. Chu 3, H. Gao 3, W. Zheng 31 Indiana University, Bloomington, Indiana 47408, USA; Center for Exploration of Energy and Matter, Indiana University, Bloomington, IN 47408, USA. 2 Department of Physics, Shanghai Jiaotong University, Shanghai,200240, China; Indiana University, Bloomington, Indiana 47408, USA; Center for Exploration of Energy and Matter, Indiana University, Bloomington, IN 47408, USA. 3 Triangle Universities Nuclear Laboratory and Department of Physics, Duke University, Durham, North Carolina 27708, USA. Received 11 June 2013; Accepted (in revised version) 27 September 2013 Available online 21 February 2014 doi:10.4208/cicp.110613.270913a Abstract We present a high precision frequency determination method for digitized NMR FID signals. The method employs high precision numerical integration rather than simple summation as in many other techniques. With no independent knowledge of the other parameters of a NMR FID signal (phase $\phi$, amplitude $A$, and transverse relaxation time $T_{2}$) this method can determine the signal frequency $f_{0}$ with a precision of $1/(8\pi^{2}f_{0}^{2}T_{2}^{2})$ if the observation time $T\gg T_{2}$. The method is especially convenient when the detailed shape of the observed FT NMR spectrum is not well defined. When $T_{2}$ is $+\infty$ and the signal becomes pure sinusoidal, the precision of the method is $3/(2\pi^{2}f_{0}^{2}T^{2})$ which is one order more precise than the $\pm1$ count error induced precision of a typical frequency counter. Analysis of this method shows that the integration reduces the noise by bandwidth narrowing as in a lock-in amplifier, and no extra signal filters are needed. For a pure sinusoidal signal we find from numerical simulations that the noise-induced error in this method reaches the Cramer-Rao Lower Band (CRLB) on frequency determination. For the damped sinusoidal case of most interest, the noise-induced error is found to be within a factor of 2 of CRLB when the measurement time $T$ is 2 or 3 times larger than $T_{2}$.We discuss possible improvements for the precision of this method. Notice: Undefined variable: ams in /var/www/html/readabs.php on line 164 PACS: 95.75.Wx, 02.70.-c, 43.50.+y Key words: NMR, FID, numerical integration, frequency determination, noise reduction. *Corresponding author. Email: haiyan@umail.iu.edu (H. Yan)