Quantum molecular dynamics calculations

Figure 14. (a) Potential-energy surfaces, with a trajectory showing the coherent vibrational motion as the diatom separates from the I atom. Two snapshots of the wavepacket motion (quantum molecular dynamics calculations) are shown for the same reaction at / = 0 and t = 600 fs. (b) Femtosecond dynamics of barrier reactions, IHgl system. Experimental observations of the vibrational (femtosecond) and rotational (picosecond) motions for the barrier (saddle-point transition state) descent, [IHgl] - Hgl(vib, rot) + I, are shown. The vibrational coherence in the reaction trajectories (oscillations) is observed in both polarizations of FTS. The rotational orientation can be seen in the decay of FTS spectra (parallel) and buildup of FTS (perpendicular) as the Hgl rotates during bond breakage (bottom). 

The setup of these calculations is very similar for both quantum and molecular mechanics. In practice, molecular dynamics calculation s using the nl) initio and semi-empirical quantum mechanical SCFmethods are limited to relatively small systems. Each time step requires a complete calculation of the wave function and the forces. 

If the inverse in Eq. (2.8) does not exist then the metric is singular, in which case the parameterization of the manifold of states is redundant. That is, the parameters are not independent, or splitting of the manifold occurs, as in potential curve crossing in quantum molecular dynamics. In both cases, the causes of the singularity must be studied and revisions made to the coordinate charts on the manifold (i.e. the way the operators are parameterized) in order to proceed with calculations. 

Fig. 13 Results from the quantum calculations on the duplex sequence 5 -GAGG-3. In a, the sodium ions and their solvating water molecules are located at positions near the phosphate anions of the DNA backbone. In b, one sodium ion is moved from near a phosphate anion to N-7 of a guanine, which molecular dynamics calculations show to be a preferred site. The balloons represent the hole density on the GAGG sequences with the two different sodium ion orientations. The radical cation clearly changes its average location with movement of the sodium ion...

This chapter focuses on the calculation of atomic charges in molecules. It discusses why atomic charges can neither be measured nor calculated unambiguously and provides two different recipes for obtaining atomic charges from quantum chemical calculations. The chapter concludes with a discussion about generating atomic charges for use in molecular mechanics/molecular dynamics calculations. 

To better understand the possible weak links in the decanio-bate ion that might give rise to oxygen exchange, a molecular dynamics calculation was carried out with the PQS quantum chemistry code (www.pqs-chem.com) of the decaniobate ion embedded in the COSMO continuum solvent model (55) (dielectric constant = 80, probe radius = 1.4 A). The time step was set at... 

Induced dipole autocorrelation functions of three-body systems have not yet been computed from first principles. Such work involves the solution of Schrodinger s equation of three interacting atoms. However, classical and semi-classical methods, especially molecular dynamics calculations, exist which offer some insight into three-body dynamics and interactions. Very useful expressions exist for the three-body spectral moments, with the lowest-order Wigner-Kirkwood quantum corrections which were discussed above. 

The chapter also examined three molecular methods (1) ab initio quantum mechanical calculations, which are typically used to get better interatomic potentials, (2) MC calculations, and (3) molecular dynamic calculations. The latter two molecular methods are most useful to probe the behavior of a small number of molecules, in which experimental capability is constrained by either space or time. 

Cimiraglia, R. (1992) Adiabatic and diabatic sets in molecular calculations, in Broeckhove. J. and Latbouwers, L. (eds), Time-dependent quantum molecular dynamics, Plenum Press, New York,pp. 11-26. 

Besides charge transfer interactions, dipolar coupling between ttk transitions of bases may lead to delocalization of the excited states. In order to obtain some guidelines for our experimental studies, we have undertaken the calculation of excited Frank-Condon states within the framework of the exci-ton theory [26]. These studies were enriched by combining data from quantum chemistry and molecular dynamics calculations in collaboration with Krystyna Zakrzewska and Richard Lavery [26,27,27-29]. The general formalism is described in the Chapter by E. Bittner and A. Czader in the present volume. 

Molecular dynamics calculations apply the laws of motion to molecules. Thus one can simulate the motion of an enzyme as it changes shape on binding to a substrate, or the motion of a swarm of water molecules around a molecule of solute quantum mechanical molecular dynamics also allows actual chemical reactions to be simulated. 

The solute charge distribution obtained from the quantum calculation is then used as input in the molecular dynamics calculation. The solute-solvent Lennard-Jones parameters and the complete solvent-solvent force field are obtained from the literature. 

Two methods are in common use for simulating molecular liquids the Monte Carlo method (MC) and molecular dynamics calculations (MD). Both depend on the availability of reasonably accurate potential energy surfaces and both are based on statistical classical mechanics, taking no account of quantum effects. In the past 10-15 years quantum Monte Carlo methods (QMC) have been developed that allow intramolecular degrees of freedom to be studied, but because of the computational complexity of this approach results have only been reported for water clusters. 

F. Barocchi, M. Zoppi, and M. Neumann. First-order quantum corrections to depolarized interaction induced light scattering spectral moments Molecular dynamics calculation. Phys. Rev. A, 27 1587-1593 (1983). 

---