Trajectory computation analysis

The term standard studies refers to those based upon tracking the time evolution of swarms of trajectories, so that correlation functions and spectra can be computed. Instead, it is recommended that the analysis be based upon stability parameters of trajectories. Stability analysis of the CH periodic orbit, as a function of stretch-wag coupling, was the subject of an earlier study by Garcia-Ayllon et al. (145). 

In a classic 1976 paper, Baldwin used the Biirgi-Dunitz trajectory to define the well-known rules for the design of cyclizations (the Baldwin rules). However, stereoelectronic factors for a bond formation to alkynes have been controversial. Originally, the rules for alkyne cyclizations were based on the assumption that nucleophiles add to alkynes along an acnte trajectory, instead of the obtuse Biirgi-Dunitz angle of attack. Subsequent experimental and computational analysis suggested that this trajectory is stereoelectronically... 

The measured ammonia and ammonium concentrations at the Finnish EMEP measurement stations of Virolahti and Ahtari show a maximum on March 21, 1989 this time corresponds to computed arrival times. However, such maximum values are not uncommon in springtime. A detailed analysis of backward trajectories, computed for these measurement stations, and the measured concentrations do not show conclusively that the measured maximum values would have been caused by the accident (Kukkonen et al., 1993). Most of the NH due to the accident may have escaped the available measurement stations. 

As a consequence of this observation, the essential dynamics of the molecular process could as well be modelled by probabilities describing mean durations of stay within different conformations of the system. This idea is not new, cf. [10]. Even the phrase essential dynamics has already been coined in [2] it has been chosen for the reformulation of molecular motion in terms of its almost invariant degrees of freedom. But unlike the former approaches, which aim in the same direction, we herein advocate a different line of method we suggest to directly attack the computation of the conformations and their stability time spans, which means some global approach clearly differing from any kind of statistical analysis based on long term trajectories. 

Another view of this theme was our analysis of spectral densities. A comparison of LN spectral densities, as computed for BPTI and lysozyme from cosine Fourier transforms of the velocity autocorrelation functions, revealed excellent agreement between LN and the explicit Langevin trajectories (see Fig, 5 in [88]). Here we only compare the spectral densities for different 7 Fig. 8 shows that the Langevin patterns become closer to the Verlet densities (7 = 0) as 7 in the Langevin integrator (be it BBK or LN) is decreased. 

In numerous cases an atomically detailed picture is required to understand function of biological molecules. The wealth of atomic information that is provided by the Molecular Dynamics (MD) method is the prime reason for its popularity and numerous successes. The MD method offers (a) qualitative understanding of atomic processes by detailed analysis of individual trajectories, and (b) comparison of computations to experimental data by averaging over a representative set of sampled trajectories. 

It is also possible to use normal mode analysis [7] to estimate the difference between the exact and the optimal trajectories. Yet another formula is based on the difference between the optimal and the exact actions 2a w [5[Yeract(t)] (f)]]- The action is computed (of course), employ-... 

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

Production. When the simulation is equilibrated, the dynamic simulation is considered reliable. Prom this point on, the trajectory generated is stored for further analysis. Typical production runs take from several hundred picoseconds up to tens of nanoseconds (depending on the size of the system and the available computer power). 

An important issue, the significance of which is sometime underestimated, is the analysis of the resulting molecular dynamics trajectories. Clearly, the value of any computer simulation lies in the quality of the information extracted from it. In fact, it is good practice to plan the analysis procedure before starting the simulation, as the goals of the analysis will often detennine the character of the simulation to be performed. 

In the MD/QM technique each tool is used separately, in an attempt to exploit their particular strengths. Classical molecular dynamics as a very fast sampling technique is first used for efficient sampling of the conformational space for the molecule of interest. A cluster analysis of the MD trajectory is then used to identify the main con-formers (clusters). Finally QM calculations, which provide a more accurate (albeit more computationally expensive) representation of the system, can be applied to just a small number of snapshots carefully extracted from each representative cluster from the MD-generated trajectory. 

Reaction rates are macroscopic averages of the number of microscopical molecules that pass from the reactant to the product valley in the potential hypersurface. An estimation of this rate can be obtained from the energy of the highest point in the reaction path, the transition state. This approach will however fail when the reaction proceeds without an enthalpic barrier or when there are many low frequency modes. The study of these cases will require the analysis of the trajectory of the molecule on the potential hypersurface. This idea constitutes the basis of molecular dynamics (MD) [96]. Molecular dynamics were traditionally too computationally demanding for transition metal complexes, but things seem now to be changing with the use of the Car-Parrinello (CP) method [97]. This approach has in fact been already succesfully applied to the study of the catalyzed polymerization of olefins [98]. 

For the description of a solution of alanine in water two models were compared and combined with one another (79), namely the continuum model approach and the cluster ansatz approach (148,149). In the cluster approach snapshots along a trajectory are harvested and subsequent quantum chemical analysis is carried out. In order to learn more about the structure and the effects of the solvent shell, the molecular dipole moments were computed. To harvest a trajectory and for comparison AIMD (here CPMD) simulations were carried out (79). The calculations contained one alanine molecule dissolved in 60 water molecules. The average dipole moments for alanine and water were derived by means of maximally localized Wannier functions (MLWF) (67-72). For the water molecules different solvent shells were selected according to the three radial pair distributions between water and the functional groups. An overview about the findings is given in Tables II and III. 

A theoretical analysis of the possible conformations of polylp-phenylene terephthalate) (PPTA) and polylp-phenylene isophthalate) (PPIA) is performed on the basis of molecular mechanics and molecular dynamics trajectories. The dependence of the persistence length on the fluctuations of the torsional angle around the ester bond is discussed for PPTA in the frame of the RIS model. Realistic parameters like bond length and bond angles are provided by computer simulations using MD. 

---