Periodic trajectories, calculation

Quasi-Newton methods may be used instead of our full Newton iteration. We have used the fast (quadratic) convergence rate of our Newton algorithm as a numerical check to discriminate between periodic and very slowly changing quasi-periodic trajectories the accurate computed elements of the Jacobian in a Newton iteration can be used in stability computations for the located periodic trajectories. There are deficiencies in the use of a full Newton algorithm, such as its sometimes small radius of convergence (Schwartz, 1983). Several other possibilities for continuation methods also exist (Doedel, 1986 Seydel and Hlavacek, 1986). The pseudo-arc length continuation was sufficient for our calculations. 

One example of non-IRC trajectory was reported for the photoisomerization of cA-stilbene.36,37 In this study trajectory calculations were started at stilbene in its first excited state. The initial stilbene structure was obtained at CIS/6-31G, and 2744 argon atoms were used as a model solvent with periodic boundary conditions. In order to save computational time, finite element interpolation method was used, in which all degrees of freedom were frozen except the central ethylenic torsional angle and the two adjacent phenyl torsional angles. The solvent was equilibrated around a fully rigid m-stilbene for 20 ps, and initial configurations were taken every 1 ps intervals from subsequent equilibration. The results of 800 trajectories revealed that, because of the excessive internal potential energy, the trajectories did not cross the barrier at the saddle point. Thus, the prerequisites for common concepts of reaction dynamics such TST or RRKM theory were not satisfied. 

H-keto G in liquid water A ground state simulation of 9H-keto G embedded in 60 H20 molecules in a periodic setup at 300 K has been performed from which six configurations have been randomly selected as input for 6 nonadiabatic surface hopping trajectory calculations starting in the S1 excited state. Comparison with the simulations in the gas phase (see Section 10.3.3.2.1) permits analysis of the effects of the water solvent on the mechanism of radiationless decay. 

The 3-phospholene potential energy surface has an energy barrier between isomers with height Ej, = 5083 cm. Results from direct trajectory calculations by De Leon and Marston [23,63] are available for one energy, namely, 5133 cm . The PSS for this energy is shown in Fig. 35. It is seen that although overall the system displays characteristics of chaotic motion, a considerable portion of the PSS supports quasi-periodic motion. 

A Poincare surface of section may be used to identify the chaotic and quasi-periodic regions of phase space for a two-dimensional Hamiltonian. An ensemble of trajectories, chosen to randomly sample the phase space, are calculated and for each trajectory a point is plotted in the (9i,Pi)-plane every time Q2 = 0 for p2 > 0. A quasi-periodic trajectory lies on an invariant curve, while the points are scattered for a chaotic trajectory with no pattern. Figure 44 shows an example for a two-dimensional model for HOCl the HO bond distance is frozen in these calculations [351]. It clearly illustrates how the phase space becomes gradually more chaotic as the energy increases. 

Exponential divergence in systems that are chaotic prevents accurate long-time trajectory calculations of their dynamics. That is, numerical errors18 propagate exponentially during the dynamics so that accuracy beyond 100 characteristic periods of motion is extraordinarily difficult to achieve thus, accurate long-time dynamics is essentially uncomputable for chaotic classical systems. This serves as additional... 

Maergoiz et al. [28-30] performed classical trajectory calculations for ion-dipole, ion-quadrupole and dipole-dipole collisions, deriving results for capture rate coefficients and expressing them in terms of two reduced parameters, the reduced temperature, d, closely related to the Tr parameter in the work of Chesnavich and co-workers, and the Massey parameter, which is equal to the ratio of the coUisional timescale to the rotational period of the neutral. 3> 1 corresponds to the adiabatic limit (see below). They gave parameterized expressions for k /ki, similar to those given by Su and Chesnavich [26], but extending the range of validity. 

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

One drawback to a molecular dynamics simulation is that the trajectory length calculated in a reasonable time is several orders of magnitude shorter than any chemical process and most physical processes, which occur in nanoseconds or longer. This allows yon to study properties that change w ithin shorter time periods (such as energy finctnations and atomic positions), but not long-term processes like protein folding. 

It was necessary periodically to generate an adiabatic trajectory in order to obtain the odd work and the time correlation functions. In calculating E (t) on a trajectory, it is essential to integrate E)(t) over the trajectory rather than use the expression for E (T(f)) given earlier. This is because is insensitive to the periodic boundary conditions, whereas j depends on whether the coordinates of the atom are confined to the central cell, or whether the itinerant coordinate is used, and problems arise in both cases when the atom leaves the central cell on a trajectory. 

---