Implicit Mass-Matrix Penalization of Hamiltonian dynamics with application to exact sampling of stiff systems

Details Implicit Mass-Matrix Penalization of Hamiltonian dynamics with application to exact sampling of stiff systems Implicit Mass-Matrix Penalization of Hamiltonian dynamics with application to exact sampling of stiff systems

2009-05-28

arxiv.org/abs/0905.4737v1

Mathematics

Numerical Analysis

39 pages

Scientific paper

An implicit mass-matrix penalization (IMMP) of Hamiltonian dynamics is proposed, and associated dynamical integrators, as well as sampling Monte-Carlo schemes, are analyzed for systems with multiple time scales. The penalization is based on an extended Hamiltonian with artificial constraints associated with some selected DOFs. The penalty parameters enable arbitrary tuning of timescales for the selected DOFs. The IMMP dynamics is shown to be an interpolation between the exact Hamiltonian dynamics and the dynamics with rigid constraints. This property translates in the associated numerical integrator into a tunable trade-off between stability and dynamical modification. Moreover, a penalty that vanishes with the time-step yields order two convergent schemes for the exact dynamics. Moreover, by construction, the resulting dynamics preserves the canonical equilibrium distribution in position variables, up to a computable geometric correcting potential, leading to Metropolis-like unbiased sampling algorithms. The algorithms can be implemented with a simple modification of standard geometric integrators with algebraic constraints imposed on the selected DOFs, and has no additional complexity in terms of enforcing the constraints and force evaluations. The properties of the IMMP method are demonstrated numerically on the $N$-alkane model, showing that the time-step stability region of integrators and the sampling efficiency can be increased with a gain that grows with the size of the system. This feature is mathematically analyzed for a harmonic atomic chain model. When a large stiffness parameter is introduced, the IMMP method is shown to be asymptotically stable and to converge towards the heuristically expected Markovian effective dynamics on the slow manifold.

Affiliated with

Also associated with

No associations

Rating

Implicit Mass-Matrix Penalization of Hamiltonian dynamics with application to exact sampling of stiff systems does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Implicit Mass-Matrix Penalization of Hamiltonian dynamics with application to exact sampling of stiff systems, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Implicit Mass-Matrix Penalization of Hamiltonian dynamics with application to exact sampling of stiff systems will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-526599