Reversible integration

NVT, and in die course of the simulation the volume V of the simulation box is allowed to vary, according to the new equations of motion. A usefid variant allows the simulation box to change shape as well as size [89, 90], It is also possible to extend the Liouville operator-splitting approach to generate algoritlnns for MD in these ensembles examples of explicit, reversible, integrators are given by Martyna et al [91],... 

Martyna G J, Tuckerman M, Tobias D J and Klein M L 1996 Explicit reversible Integrators for extended systems dynamics Mol. Phys. 87 1117-57... 

Note that the same results have not been shown for symmetric (time-reversible) integration methods, although symmetric methods seem to perform quite well in practice. For a discussion of symmetric methods in the context of the QCMD model see [16, 17, 13]. 

Martyna GJ, Tuckerman ME, Tobias DJ, Klein ML (1996) Explicit reversible integrators for extended systems dynamics. Mol Phys 87(5) 1117-1157... 

After the momenta are selected from the distribution (8.39), the dynamics is propagated by a standard leapfrog algorithm (any symplectic and time-reversible integrator is suitable). The move is then accepted or rejected according to a criterion based on the detailed balance condition... 

An expanded formulation of the steady-state permeation model has been presented. Two numerical problems - stiffness and an ill-conditioned boundary value problem - are encountered in solving the system equations. These problems can be circumvented by matching forward and reverse integrations at a point near the inlet (n = 0) but outside the combustion zone. The model predicts a... 

The last question that must be addressed is this How does one obtain a stable integrator when P exceeds, say, 10 steps The answer to this question is due to a remarkable theorem in the theory of reversible integrators. We will simply quote the theorem here and interpret it. The interested reader can pursue the details in Refs. 48 and 49. If the Jacobian of the map (integrator)... 

Hamiltonian systems. Thus, one has to treat this non-volume-preserving piece of the integrator a bit more carefully. To ensure numerical stability, higher order reversible integration schemes in conjunction with multiple time step methods are preferred. The details of implementing this scheme are provided in Ref. 28. 

To construct a reversible integrator for Eqs. [218] with a NH thermostat, we write the total Liouvillian in the form... 

G. J. Martyna, M. E. Tuckerman, D. J. Tobias, andM. L. Klein, Mo/. Phys., 87,1117 (1996). Explicit Reversible Integrators for Extended System Molecular Dynamics. 

It is also possible to extend the Liouville operator-splitting approach to generate algorithms for MD in these ensembles examples of explicit, reversible, integrators are given by Martyna et al [91],... 

Of particular importance for molecular dynamics are the following properties a symplectic map will preserve volume, whereas a time-reversible map need not do so, and a symplectic integrator will approximately conserve energy due to the existence of the perturbed Hamiltonian, whereas a time-reversible integrator may... 

Cheng and Luyben (1985) compared the steady state designs for several multi-effect configurations for a benzene/toluene/xylene separation. The result showed 50% energy reduction and that the best configuration for the separation studied was a prefractionator/sidestream column arrangement with reverse integration. 

Nos developed a deterministic approach to constant temperature simulations based on an extended Lagrangian which does not disturb the dynamic properties of the system. Since then, modifications of this approach have been developed to generate constant pressure and temperature simulations. These extended system methods have been known to suffer from stability problems, as well as occasional failures in ergodicity. Klein and co-workers introduced the concept of Nos6-Hoover chains to overcome these problems, and used a Liouville operator formalism to obtain reversible integrators generating these chains. 

A major organisational change such as this needs to be discussed carefiiUy with staff the idea cannot be successfully presented at a one-off meeting. A phased programme of introduction should be considered and even reverse integration from the mainstream to the special school prior to the project as a preparation for staff and students of both schools. 

---