The Langevin/implicit-Euler/normal-mode scheme for molecular dynamics at large time steps

Guihua Zhang, Tamar Schlick

Research output: Contribution to journalArticle

Abstract

As molecular dynamics simulations continue to provide important insights into biomolecular structure and function, a growing demand for increasing the time span of the simulations is emerging. Our focus here is developing a new algorithm, LIN (Langevin/implicit-Euler/normal mode), that combines normal-mode and implicit-integration techniques, for large time step biomolecular applications. In the normal-mode phase of LIN, we solve an approximate linearized Langevin formulation to resolve the rapidly varying components of the motion. In the implicit phase, we resolve the remaining components of the motion by numerical integration with the implicit-Euler scheme. Developments of the normal-mode phase of LIN are discussed in this paper. Specifically, we solve two crucial issues of the method. The first involves how to choose and how often to update the Hessian approximation for the linearized Langevin equation. This approximation must be computationally feasible and physically reasonable to capture the motion in the higher end of the vibrational spectrum. Three such general Hessian approximations are discussed. The related issue - the frequency of the Hessian update - is analyzed by projecting the motion onto the different vibrational modes. This analysis demonstrates that a one-picosecond interval is reasonable for updating the Hessian in the model system examined here. In this connection, we illustrate that the high-frequency motions are highly localized while the low-frequency motions are delocalized. We also show rigorously that the mode amplitudes are inversely proportional to the frequency (consistent with the equipartition theorem), with 90% of the displacement fluctuations coming from a very small group of low-frequency modes. Anharmonic effects essentially influence the low-frequency modes. The second issue involves how to solve the linearized Langevin equation at large timesteps correctly, where the usual discretized formulation of the random force is invalid. This is accomplished by using analytic expressions for the distributions associated with positions and velocities of the individual oscillators as a function of frequency, obtained as the solution of the corresponding Fokker-Planck equation. We apply LIN with these developments to the nucleic acid component deoxycytidine with timesteps ranging from 100 to 1000 fs. We demonstrate that LIN is stable in these simulations, with energies fluctuating about the same values - and possessing overall similar dynamical features - in comparison to 1 fs explicit simulations, though the fluctuations are significantly larger at larger timesteps. Moreover, continuous dynamics is maintained, and pathway information can be obtained. Computational performance is competitive only at very large time steps: a gain factor of 3-4 is obtained for runs with 1000 fs time steps. Larger gains may be achieved for biomolecules, where sparsity and parallelization can be exploited significantly.

Original languageEnglish (US)
Pages (from-to)4995-5012
Number of pages18
JournalThe Journal of chemical physics
Volume101
Issue number6
StatePublished - 1994

Fingerprint

Molecular dynamics
molecular dynamics
Fokker Planck equation
Deoxycytidine
Biomolecules
Vibrational spectra
Nucleic Acids
Computer simulation
low frequencies
simulation
approximation
equipartition theorem
formulations
nucleic acids
Fokker-Planck equation
numerical integration
vibrational spectra
vibration mode
emerging
oscillators

ASJC Scopus subject areas

  • Atomic and Molecular Physics, and Optics

Cite this

The Langevin/implicit-Euler/normal-mode scheme for molecular dynamics at large time steps. / Zhang, Guihua; Schlick, Tamar.

In: The Journal of chemical physics, Vol. 101, No. 6, 1994, p. 4995-5012.

Research output: Contribution to journalArticle

@article{61faa902bee1440189d548f4b5a076ae,
title = "The Langevin/implicit-Euler/normal-mode scheme for molecular dynamics at large time steps",
abstract = "As molecular dynamics simulations continue to provide important insights into biomolecular structure and function, a growing demand for increasing the time span of the simulations is emerging. Our focus here is developing a new algorithm, LIN (Langevin/implicit-Euler/normal mode), that combines normal-mode and implicit-integration techniques, for large time step biomolecular applications. In the normal-mode phase of LIN, we solve an approximate linearized Langevin formulation to resolve the rapidly varying components of the motion. In the implicit phase, we resolve the remaining components of the motion by numerical integration with the implicit-Euler scheme. Developments of the normal-mode phase of LIN are discussed in this paper. Specifically, we solve two crucial issues of the method. The first involves how to choose and how often to update the Hessian approximation for the linearized Langevin equation. This approximation must be computationally feasible and physically reasonable to capture the motion in the higher end of the vibrational spectrum. Three such general Hessian approximations are discussed. The related issue - the frequency of the Hessian update - is analyzed by projecting the motion onto the different vibrational modes. This analysis demonstrates that a one-picosecond interval is reasonable for updating the Hessian in the model system examined here. In this connection, we illustrate that the high-frequency motions are highly localized while the low-frequency motions are delocalized. We also show rigorously that the mode amplitudes are inversely proportional to the frequency (consistent with the equipartition theorem), with 90{\%} of the displacement fluctuations coming from a very small group of low-frequency modes. Anharmonic effects essentially influence the low-frequency modes. The second issue involves how to solve the linearized Langevin equation at large timesteps correctly, where the usual discretized formulation of the random force is invalid. This is accomplished by using analytic expressions for the distributions associated with positions and velocities of the individual oscillators as a function of frequency, obtained as the solution of the corresponding Fokker-Planck equation. We apply LIN with these developments to the nucleic acid component deoxycytidine with timesteps ranging from 100 to 1000 fs. We demonstrate that LIN is stable in these simulations, with energies fluctuating about the same values - and possessing overall similar dynamical features - in comparison to 1 fs explicit simulations, though the fluctuations are significantly larger at larger timesteps. Moreover, continuous dynamics is maintained, and pathway information can be obtained. Computational performance is competitive only at very large time steps: a gain factor of 3-4 is obtained for runs with 1000 fs time steps. Larger gains may be achieved for biomolecules, where sparsity and parallelization can be exploited significantly.",
author = "Guihua Zhang and Tamar Schlick",
year = "1994",
language = "English (US)",
volume = "101",
pages = "4995--5012",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "6",

}

TY - JOUR

T1 - The Langevin/implicit-Euler/normal-mode scheme for molecular dynamics at large time steps

AU - Zhang, Guihua

AU - Schlick, Tamar

PY - 1994

Y1 - 1994

N2 - As molecular dynamics simulations continue to provide important insights into biomolecular structure and function, a growing demand for increasing the time span of the simulations is emerging. Our focus here is developing a new algorithm, LIN (Langevin/implicit-Euler/normal mode), that combines normal-mode and implicit-integration techniques, for large time step biomolecular applications. In the normal-mode phase of LIN, we solve an approximate linearized Langevin formulation to resolve the rapidly varying components of the motion. In the implicit phase, we resolve the remaining components of the motion by numerical integration with the implicit-Euler scheme. Developments of the normal-mode phase of LIN are discussed in this paper. Specifically, we solve two crucial issues of the method. The first involves how to choose and how often to update the Hessian approximation for the linearized Langevin equation. This approximation must be computationally feasible and physically reasonable to capture the motion in the higher end of the vibrational spectrum. Three such general Hessian approximations are discussed. The related issue - the frequency of the Hessian update - is analyzed by projecting the motion onto the different vibrational modes. This analysis demonstrates that a one-picosecond interval is reasonable for updating the Hessian in the model system examined here. In this connection, we illustrate that the high-frequency motions are highly localized while the low-frequency motions are delocalized. We also show rigorously that the mode amplitudes are inversely proportional to the frequency (consistent with the equipartition theorem), with 90% of the displacement fluctuations coming from a very small group of low-frequency modes. Anharmonic effects essentially influence the low-frequency modes. The second issue involves how to solve the linearized Langevin equation at large timesteps correctly, where the usual discretized formulation of the random force is invalid. This is accomplished by using analytic expressions for the distributions associated with positions and velocities of the individual oscillators as a function of frequency, obtained as the solution of the corresponding Fokker-Planck equation. We apply LIN with these developments to the nucleic acid component deoxycytidine with timesteps ranging from 100 to 1000 fs. We demonstrate that LIN is stable in these simulations, with energies fluctuating about the same values - and possessing overall similar dynamical features - in comparison to 1 fs explicit simulations, though the fluctuations are significantly larger at larger timesteps. Moreover, continuous dynamics is maintained, and pathway information can be obtained. Computational performance is competitive only at very large time steps: a gain factor of 3-4 is obtained for runs with 1000 fs time steps. Larger gains may be achieved for biomolecules, where sparsity and parallelization can be exploited significantly.

AB - As molecular dynamics simulations continue to provide important insights into biomolecular structure and function, a growing demand for increasing the time span of the simulations is emerging. Our focus here is developing a new algorithm, LIN (Langevin/implicit-Euler/normal mode), that combines normal-mode and implicit-integration techniques, for large time step biomolecular applications. In the normal-mode phase of LIN, we solve an approximate linearized Langevin formulation to resolve the rapidly varying components of the motion. In the implicit phase, we resolve the remaining components of the motion by numerical integration with the implicit-Euler scheme. Developments of the normal-mode phase of LIN are discussed in this paper. Specifically, we solve two crucial issues of the method. The first involves how to choose and how often to update the Hessian approximation for the linearized Langevin equation. This approximation must be computationally feasible and physically reasonable to capture the motion in the higher end of the vibrational spectrum. Three such general Hessian approximations are discussed. The related issue - the frequency of the Hessian update - is analyzed by projecting the motion onto the different vibrational modes. This analysis demonstrates that a one-picosecond interval is reasonable for updating the Hessian in the model system examined here. In this connection, we illustrate that the high-frequency motions are highly localized while the low-frequency motions are delocalized. We also show rigorously that the mode amplitudes are inversely proportional to the frequency (consistent with the equipartition theorem), with 90% of the displacement fluctuations coming from a very small group of low-frequency modes. Anharmonic effects essentially influence the low-frequency modes. The second issue involves how to solve the linearized Langevin equation at large timesteps correctly, where the usual discretized formulation of the random force is invalid. This is accomplished by using analytic expressions for the distributions associated with positions and velocities of the individual oscillators as a function of frequency, obtained as the solution of the corresponding Fokker-Planck equation. We apply LIN with these developments to the nucleic acid component deoxycytidine with timesteps ranging from 100 to 1000 fs. We demonstrate that LIN is stable in these simulations, with energies fluctuating about the same values - and possessing overall similar dynamical features - in comparison to 1 fs explicit simulations, though the fluctuations are significantly larger at larger timesteps. Moreover, continuous dynamics is maintained, and pathway information can be obtained. Computational performance is competitive only at very large time steps: a gain factor of 3-4 is obtained for runs with 1000 fs time steps. Larger gains may be achieved for biomolecules, where sparsity and parallelization can be exploited significantly.

UR - http://www.scopus.com/inward/record.url?scp=0001299323&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0001299323&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0001299323

VL - 101

SP - 4995

EP - 5012

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 6

ER -