Data-based parameter estimation of generalized multidimensional Langevin processes.

The generalized Langevin equation is useful for modeling a wide range of physical processes. Unfortunately its parameters, especially the memory function, are difficult to determine for nontrivial processes. We establish relations between a time-discrete generalized Langevin model and discrete multivariate autoregressive (AR) or autoregressive moving average models (ARMA). This allows a wide range of discrete linear methods known from time series analysis to be applied. In particular, the determination of the memory function via the order of the respective AR or ARMA model is addressed. The method is illustrated on a one-dimensional test system and subsequently applied to the molecular dynamics time series of a biomolecule that exhibits an interesting relationship between the solvent method used, the respective molecular conformation, and the depth of the memory.

[1]  Gary P. Morriss,et al.  Statistical Mechanics of Nonequilibrium Liquids , 2008 .

[2]  Progress in Turbulence VI , 2005 .

[3]  H. Mori Transport, Collective Motion, and Brownian Motion , 1965 .

[4]  Causal relations in molecular dynamics from the multi-variate autoregressive model , 2006 .

[5]  R. Shanmugam Introduction to Time Series and Forecasting , 1997 .

[6]  A. Chorin,et al.  Stochastic Tools in Mathematics and Science , 2005 .

[7]  Dror Givon,et al.  Existence proof for orthogonal dynamics and the Mori-Zwanzig formalism , 2005 .

[8]  Oliver F. Lange,et al.  Collective Langevin dynamics of conformational motions in proteins. , 2006, The Journal of chemical physics.

[9]  T. DelSole,et al.  Optimally Persistent Patterns in Time-Varying Fields , 2001 .

[10]  L. Garrido,et al.  Systems far from equilibrium , 1980 .

[11]  S. Koda,et al.  Molecular dynamics simulation study on the transient response of solvation structure during the translational diffusion of solute. , 2005, The Journal of chemical physics.

[12]  M. Karplus,et al.  A Comprehensive Analytical Treatment of Continuum Electrostatics , 1996 .

[13]  Gwilym M. Jenkins,et al.  Time series analysis, forecasting and control , 1971 .

[14]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[15]  Masafumi Takahashi,et al.  Non-ideal Brownian motion, generalized Langevin Equation and its application to the security market , 1996 .

[16]  W. Jost,et al.  Physical Chemistry, An Advanced Treatise , 1974 .

[17]  Robert Zwanzig,et al.  Problems in nonlinear transport theory , 1980 .

[18]  Alexander Fischer,et al.  Automated Model Reduction for Complex Systems Exhibiting Metastability , 2006, Multiscale Model. Simul..

[19]  Illia Horenko,et al.  Likelihood-Based Estimation of Multidimensional Langevin Models and Its Application to Biomolecular Dynamics , 2008, Multiscale Model. Simul..

[20]  James Witkoskie,et al.  Basis set study of classical rotor lattice dynamics. , 2004, The Journal of chemical physics.

[21]  Richard A. Davis,et al.  Introduction to time series and forecasting , 1998 .

[22]  J. D. Doll,et al.  Generalized Langevin equation approach for atom/solid-surface scattering: General formulation for classical scattering off harmonic solids , 1976 .

[23]  R. Zwanzig Nonlinear generalized Langevin equations , 1973 .

[24]  Wasfy B. Mikhael,et al.  Multidimensional model based speech signal representations for automatic speaker identification , 2004, Circuits, Signals, and Systems.

[25]  James C. McWilliams,et al.  Langevin models of turbulence: Renormalization group, distant interaction algorithms or rapid distortion theory? , 2003 .

[26]  J. Deutch,et al.  Chapter 11 – Nonequilibrium Problems—Projection Operator Techniques , 1975 .

[27]  Andrew M. Stuart,et al.  Fitting SDE models to nonlinear Kac–Zwanzig heat bath models , 2004 .

[28]  Steven M. Kay,et al.  Vector space solution to the multidimensional Yule-Walker equations , 2003, 2003 IEEE International Conference on Acoustics, Speech, and Signal Processing, 2003. Proceedings. (ICASSP '03)..