Trajectory calculations, multiple

The criterion for association in the classical trajectory calculations is the presence of multiple inner turning points between Cl and the CH3Clb center of mass. Two types of Cl -CH3Clb complexes, short-range and long-range, were observed in the trajectories. The two different complexes were identified by the Cl -C distance... 

Quantum chemical methods, exemplified by CASSCF and other MCSCF methods, have now evolved to an extent where it is possible to routinely treat accurately the excited electronic states of molecules containing a number of atoms. Mixed nuclear dynamics, such as swarm of trajectory based surface hopping or Ehrenfest dynamics, or the Gaussian wavepacket based multiple spawning method, use an approximate representation of the nuclear wavepacket based on classical trajectories. They are thus able to use the infoiination from quantum chemistry calculations required for the propagation of the nuclei in the form of forces. These methods seem able to reproduce, at least qualitatively, the dynamics of non-adiabatic systems. Test calculations have now been run using duect dynamics, and these show that even a small number of trajectories is able to produce useful mechanistic infomiation about the photochemistry of a system. In some cases it is even possible to extract some quantitative information. 

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

A calculation of f in a two-dimensional lattice that takes into account multiple return vacancy trajectories appears in Exercise 8.8. 

The protein-solvent interface was studied in an explicit solvent environment of 3182 water molecules by MD simulations performed on metmyoglobin [31].Both the structure and dynamics of the hydrated surface of myoglobin are similar to that obtained by experimental methods calculating three-dimensional density distributions, temperature factors and occupancy weights of the solvent molecules. On the basis of trajectories they identified multiple solvation layers around the protein surface including more than 500 hydration sites. Properties of theoretically calculated hydration clusters were compared to that obtained from neutron and X-ray data. This study indicates that the simulation unified the hydration picture provided by X-ray and neutron diffraction experiments. 

---