Take the 2-minute tour ×
Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. It's 100% free, no registration required.

Okumura et al., J. Chem. Phys. 2007 states that the Gear predictor-corrector integration scheme, used in particular in some molecular dynamics packages for the dynamics of rigid bodies using quaternions to represent molecular orientations, is not a symplectic scheme. My question is: how can one prove that? Does it follow from the fact that the Gear integrator is not time-reversible (and if so, how can one show that)? If not, how do you prove that an integration scheme is not symplectic?

share|improve this question
add comment

1 Answer

up vote 15 down vote accepted

Take a look at the notes on lectures 1 and 2 of Geometric Numerical Integration found here. Quoting from Lecture 2

A numerical one-step method $y_{n+1} = \Phi_h(y_n)$ is called symplectic if, when applied to a Hamiltonian system, the discrete flow $y \mapsto \Phi_h(y)$ is a symplectic transformation for all sufficiently small step sizes.

From your link you have $$x(t+h) = x(t) + h \dot{x}(t) + h^2 \left\{\frac{3}{24}f(t+h) +\frac{10}{24}f(t) -\frac{1}{24}f(t-h) \right\}$$ and $$\dot{x}(t+h) = \frac{x(t+h) - x(t)}{h} + h \dot{x}(t) + h \left\{\frac{7}{24}f(t+h) +\frac{6}{24}f(t) -\frac{1}{24}f(t-h) \right\}$$

Now take $\omega(\xi,\eta) = \xi^T J \eta$ where $J = \left(\begin{array}{cc} 0 & \mathbb{I} \\ \mathbb{I} & 0 \end{array}\right)$. Then the integrator is symplectic if and only if $\omega(x(t),\dot{x}(t))=\omega(x(t+h),\dot{x}(t+h))$ for sufficiently small $h$.

All that you need to do is to fill in the values of $x(t+h)$ and $\dot{x}(t+h)$ from the integrator, and show that this condition does not hold.

share|improve this answer
    
Dear Joe Fitzsimons, you should insert a minus sign into the matrix $J$ representing the symplectic form $\omega$. –  Giuseppe Jan 21 '12 at 15:09
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.