Integral motion

Integrated motion of molecular machines (rotaxane- and catenane-type) in supramolecular polymeric scaffolds 13PC2395. 

Figure Al.2.7. Trajectory of two coupled stretches, obtained by integrating Hamilton s equations for motion on a PES for the two modes. The system has stable anhamionic synmretric and antisyimnetric stretch modes, like those illustrated in figrne Al.2.6. In this trajectory, semiclassically there is one quantum of energy in each mode, so the trajectory corresponds to a combination state with quantum numbers nj = [1, 1]. The woven pattern shows that the trajectory is regular rather than chaotic, corresponding to motion in phase space on an invariant torus.

If the adiabatic work is independent of the path, it is the integral of an exact differential and suffices to define a change in a function of the state of the system, the energy U. (Some themiodynamicists call this the internal energy , so as to exclude any kinetic energy of the motion of the system as a whole.)... 

In the example of the previous section, the release of the stop always leads to the motion of the piston in one direction, to a final state in which the pressures are equal, never in the other direction. This obvious experimental observation turns out to be related to a mathematical problem, the integrability of differentials in themiodynamics. The differential Dq, even is inexact, but in mathematics many such expressions can be converted into exact differentials with the aid of an integrating factor. 

Wlien H has reached its minimum value this is the well known Maxwell-Boltzmaim distribution for a gas in themial equilibrium with a unifomi motion u. So, argues Boltzmaim, solutions of his equation for an isolated system approach an equilibrium state, just as real gases seem to do. Up to a negative factor (-/fg, in fact), differences in H are the same as differences in the themiodynamic entropy between initial and final equilibrium states. Boltzmaim thought that his //-tiieorem gave a foundation of the increase in entropy as a result of the collision integral, whose derivation was based on the Stosszahlansatz. 

To generalize what we have just done to reactive and inelastic scattering, one needs to calculate numerically integrated trajectories for motions in many degrees of freedom. This is most convenient to develop in space-fixed Cartesian coordinates. In this case, the classical equations of motion (Hamilton s equations) are given... 

Electronic spectra are almost always treated within the framework of the Bom-Oppenlieimer approxunation [8] which states that the total wavefiinction of a molecule can be expressed as a product of electronic, vibrational, and rotational wavefiinctions (plus, of course, the translation of the centre of mass which can always be treated separately from the internal coordinates). The physical reason for the separation is that the nuclei are much heavier than the electrons and move much more slowly, so the electron cloud nonnally follows the instantaneous position of the nuclei quite well. The integral of equation (BE 1.1) is over all internal coordinates, both electronic and nuclear. Integration over the rotational wavefiinctions gives rotational selection rules which detemiine the fine structure and band shapes of electronic transitions in gaseous molecules. Rotational selection rules will be discussed below. For molecules in condensed phases the rotational motion is suppressed and replaced by oscillatory and diflfiisional motions. 

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],... 

Ryckaert J-P, Ciccotti G and Berendsen H J C 1977 Numerical integration of the Cartesian equations of motion of a system with constraints molecular dynamics of n-alkanes J. Comput. Phys. 23 327-41... 

Molecular dynamics tracks tire temporal evolution of a microscopic model system tlirough numerical integration of tire equations of motion for tire degrees of freedom considered. The main asset of molecular dynamics is tliat it provides directly a wealtli of detailed infonnation on dynamical processes. 

The center of the wavepacket thus evolves along the trajectory defined by classical mechanics. This is in fact a general result for wavepackets in a hannonic potential, and follows from the Ehrenfest theorem [147] [see Eqs. (154,155) in Appendix C]. The equations of motion are straightforward to integrate, with the exception of the width matrix, Eq. (44). This equation is numerically unstable, and has been found to cause problems in practical applications using Morse potentials [148]. As a result, Heller inboduced the P-Z method as an alternative propagation method [24]. In this, the matrix A, is rewritten as a product of matrices... 

We may conclude that the matter of optimal algorithms for integrating Newton s equations of motion is now nearly settled however, their optimal and prudent use [28] has not been fully exploited yet by most programs and may still give us an improvement by a factor 3 to 5. 

In order to solve the classical equations of motion numerically, and, thus, to t)btain the motion of all atoms the forces acting on every atom have to be computed at each integration step. The forces are derived from an energy function which defines the molecular model [1, 2, 3]. Besides other important contributions (which we shall not discuss here) this function contains the Coulomb sum... 

This discussion suggests that even the reference trajectories used by symplectic integrators such as Verlet may not be sufficiently accurate in this more rigorous sense. They are quite reasonable, however, if one requires, for example, that trajectories capture the spectral densities associated with the fastest motions in accord to the governing model [13, 15]. Furthermore, other approaches, including nonsymplectic integrators and trajectories based on stochastic differential equations, can also be suitable in this case when carefully formulated. 

The LIN method ( Langevin/Implicit/Normal-Modes ) combines frequent solutions of the linearized equations of motions with anharmonic corrections implemented by implicit integration at a large timestep. Namely, we express the collective position vector of the system as X t) = Xh t) + Z t). (In LN, Z t) is zero). The first part of LIN solves the linearized Langevin equation for the harmonic reference component of the motion, Xh t)- The second part computes the residual component, Z(t), with a large timestep. 

The solution Xh(t) of the linearized equations of motion can be solved by standard NM techniques or, alternatively, by explicit integration. We have experimented with both and found the second approach to be far more efficient and to work equally well. Its handling of the random force discretization is also more straightforward (see below). For completeness, we describe both approaches here. 

We further discuss how quantities typically measured in the experiment (such as a rate constant) can be computed with the new formalism. The computations are based on stochastic path integral formulation [6]. Two different sources for stochasticity are considered. The first (A) is randomness that is part of the mathematical modeling and is built into the differential equations of motion (e.g. the Langevin equation, or Brownian dynamics). The second (B) is the uncertainty in the approximate numerical solution of the exact equations of motion. 

Another difference is related to the mathematical formulation. Equation (1) is deterministic and does not include explicit stochasticity. In contrast, the equations of motion for a Brownian particle include noise. Nevertheless, similar algorithms are adopted to solve the two differential equations as outlined below. The most common approach is to numerically integrate the above differential equations using small time steps and preset initial values. 

The constrained equations of motion in cartesian eoordinates can be solved by the SHAKE or (the essentially equivalent) RATTLE method (see [8]) which requires the solution of a non-linear system of equations in the Lagrange multiplier funetion A. The equivalent formulation in local coordinates ean still be integrated by using the explicit Verlet method. 

One property of the exact trajectory for a conservative system is that the total energy is a constant of the motion. [12] Finite difference integrators provide approximate solutions to the equations of motion and for trajectories generated numerically the total energy is not strictly conserved. The exact trajectory will move on a constant energy surface in the 61V dimensional phase space of the system defined by. 

---