J. Phys. Chem. A, 107 (24), 4869 -4876, 2003. 10.1021/jp0261758 S1089-5639(02)06175-3
Web Release Date: May 23, 2003

Copyright © 2003 American Chemical Society

Application of the Hilbert-Huang Transform to the Analysis of Molecular Dynamics Simulations

Stephen C. Phillips, Robert J. Gledhill, and Jonathan W. Essex*

Department of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, U.K.

Colin M. Edge

GlaxoSmithKline, New Frontiers Science Park (North), Coldharbour Road, Harlow CM19 5AD, U.K.

Received: May 24, 2002

In Final Form: April 15, 2003

Abstract:

The Hilbert-Huang transform (HHT) is a new method for the analysis of nonstationary signals that allows a signal's frequency and amplitude to be evaluated with excellent time resolution. In this paper, the HHT method is described, and its performance is compared with the Fourier methods of spectral analysis. The HHT is then applied to the analysis of molecular dynamics simulation trajectories, including enhanced sampling trajectories produced by reversible digitally filtered molecular dynamics. Amplitude-time, amplitude-frequency, and amplitude-frequency-time spectra are all produced with the method and compared to equivalent results obtained using wavelet analysis. The wavelet and HHT analysis yield qualitatively similar results, but the HHT provides a better match to physical intuition than the wavelet transform. Moreover the HHT method is able to show the flow of energy into low-frequency vibrations during conformational change events and is able to identify frequencies appropriate for amplification by digital filters including the observation of a 10 cm-1 shift in target frequency.


Download the full text: PDF | HTML