Trajectory complex

Three-dimensional Q aU ions maintained in stable trajectory (complex sinusoidal path) within the trap... 

In Chemical Kinetics one considers reactions evolving close to equilibrium. Attention is focused on detailed reaction mechanisms by which given initial species are transformed into final products, and equilibrate with them. The dynamics is rather simple, not in term in reaction mechanism, but in term of trajectory complexity. Most of the time, one essentially observes relaxation towards equilibrium. It is precisely the outcome of Thermodynamics that equilibrium is not only stable, but also that, because of the principle of detailed balance, which is part of the equilibrium concept, it is stable in a way that precludes any really complex behaviour. Relaxation towards equilibrium has to be monotonous. ... 

The Mathieu equation for the quadnipole ion trap again has stable, bounded solutions conesponding to stable, bounded trajectories inside the trap. The stability diagram for the ion trap is quite complex, but a subsection of the diagram, correspondmg to stable trajectories near the physical centre of the trap, is shown in figure Bl.7.15. The interpretation of the diagram is similar to that for tire quadnipole mass filter. 

The simplest way to add a non-adiabatic correction to the classical BO dynamics method outlined above in Section n.B is to use what is known as surface hopping. First introduced on an intuitive basis by Bjerre and Nikitin [200] and Tully and Preston [201], a number of variations have been developed [202-205], and are reviewed in [28,206]. Reference [204] also includes technical details of practical algorithms. These methods all use standard classical trajectories that use the hopping procedure to sample the different states, and so add non-adiabatic effects. A different scheme was introduced by Miller and George [207] which, although based on the same ideas, uses complex coordinates and momenta. 

To exemplify both aspects of the formalism and for illustration purposes, we divide the present manuscript into two major parts. We start with calculations of trajectories using approximate solution of atomically detailed equations (approach B). We then proceed to derive the equations for the conditional probability from which a rate constant can be extracted. We end with a simple numerical example of trajectory optimization. More complex problems are (and will be) discussed elsewhere [7]. 

The model consists of a two dimensional harmonic oscillator with mass 1 and force constants of 1 and 25. In Fig. 1 we show trajectories of the two oscillators computed with two time steps. When the time step is sufficiently small compared to the period of the fast oscillator an essentially exact result is obtained. If the time step is large then only the slow vibration persists, and is quite accurate. The filtering effect is consistent (of course) with our analytical analysis. Similar effects were demonstrated for more complex systems [7]. 

Computational issues that are pertinent in MD simulations are time complexity of the force calculations and the accuracy of the particle trajectories including other necessary quantitative measures. These two issues overwhelm computational scientists in several ways. MD simulations are done for long time periods and since numerical integration techniques involve discretization errors and stability restrictions which when not put in check, may corrupt the numerical solutions in such a way that they do not have any meaning and therefore, no useful inferences can be drawn from them. Different strategies such as globally stable numerical integrators and multiple time steps implementations have been used in this respect (see [27, 31]). 

MacroModel (we tested Version 6.5) is a powerful molecular mechanics program. The program can be run from either its graphic interface or an ASCII command file. The command file structure allows very complex simulations to be performed. The XCluster utility permits the analysis and filtering of a large number of structures, such as Monte Carlo or dynamics trajectories. The documentation is very thorough. 

Most potential energy surfaces are extremely complex. Fiber and Karplus analyzed a 300 psec molecular dynamics trajectory of the protein myoglobin. They estimate that 2000 thermally accessible minima exist near the native protein structure. The total number of conformations is even larger. Dill derived a formula to calculate the upper bound of thermally accessible conformations in a protein. Using this formula, a protein of 150 residues (the approx-... 

Unlike simple deflections or accelerations of ions in magnetic and electric fields (Chapter 25), the trajectory of an ion in a quadrupolar field is complex, and the equations of motion are less easy to understand. Accordingly, a simplified version of the equations is given here, with a fuller discussion in the Appendix at the end of this chapter. 

Time reversibility. Newton s equation is reversible in time. Eor a numerical simulation to retain this property it should be able to retrace its path back to the initial configuration (when the sign of the time step At is changed to —At). However, because of chaos (which is part of most complex systems), even modest numerical errors make this backtracking possible only for short periods of time. Any two classical trajectories that are initially very close will eventually exponentially diverge from one another. In the same way, any small perturbation, even the tiny error associated with finite precision on the computer, will cause the computer trajectories to diverge from each other and from the exact classical trajectory (for examples, see pp. 76-77 in Ref. 6). Nonetheless, for short periods of time a stable integration should exliibit temporal reversibility. 

From the obtained trajectories the following structural and electronic properties were extracted (i) structure factors, (ii) stability of Sn complexes, in particular, the Zintl anions, and (iii) the electrical conductivities - which yield answers to the questions mentioned in the first paragraph. 

Lloyd and Pagels show that these three requirements lead uniquely to an average complexity of a state proportional to the Shannon entropy of the set of (experimentally determined) trajectories leading to the given state (= EiPi oSzPi)- The thermodynamic depth of a state S to which the system S has evolved via the possible trajectory is equal to the amount of information required to specify that trajectory, or Djj S) Hamiltonian systems, Lloyd and... 

Not every collision, not every punctilious trajectory by which billiard-ball complexes arrive at their calculable meeting places leads to reaction. 

---