Molecular dynamics tunneling

Mavri, J., Berendsen, H.J.C., Van Gunsteren, W.F. Influence of solvent on intramolecular proton transfer in hydrogen malonate. Molecular dynamics study of tunneling by density matrix evolution and nonequilibrium solvation. J. Phys. Chem. 97 (1993) 13469-13476. 

The development of an ab initio quantum molecular dynamics method is guided by the need to overcome two main obstacles. First, one needs to develop an efficient, yet accurate, method for solving the electronic Schrodinger equation for both ground and excited electronic states. Second, the quantum mechanical character of the nuclear dynamics must be addressed. (This is necessary for the description of photochemical and tunneling processes.) This section provides a detailed discussion of the approaches we have taken to solve these two problems. 

Semiclassical techniques like the instanton approach [211] can be applied to tunneling splittings. Finally, one can exploit the close correspondence between the classical and the quantum treatment of a harmonic oscillator and treat the nuclear dynamics classically. From the classical trajectories, correlation functions can be extracted and transformed into spectra. The particular charm of this method rests in the option to carry out the dynamics on the fly, using Born Oppenheimer or fictitious Car Parrinello dynamics [212]. Furthermore, multiple minima on the hypersurface can be treated together as they are accessed by thermal excitation. This makes these methods particularly useful for liquid state or other thermally excited system simulations. Nevertheless, molecular dynamics and Monte Carlo simulations can also provide insights into cold gas-phase cluster formation [213], if a reliable force field is available [189]. 

The concept of tunneling has recently been invoked to explain the mechanism of photodissociation of matrix-isolated molecules. Previously, photodissociation was customarily accounted for by the fact that transla-tionally hot photofragments escape from the cage and stabilize in separate matrix sites, thereby avoiding recombination. Using the time-dependent self-consistent field approximation for molecular dynamics simulations,... 

Thermochemistry DSC (see Box 9.1) TGA (see Box 9.1) Vapour pressure measurements Solid-state NMR Ultramicroscopy Scanning electron microscopy (SEM) Tunnelling electron microscopy (TEM) Atomic force microscopy (AFM) Molecular dynamics... 

Nonpolar and dipolar altitudinal rotors (compounds 2 and 3 in Fig. 17.3) have been synthesized. 19F NMR spectroscopy showed that the barrier to rotation in 3 was extremely low in solution. Both systems have then been immobilized on Au(l 11) surfaces and studied with a variety of techniques.57 The results obtained indicated that for a fraction of molecules the static electric field from the scanning tunneling microscopy (STM) tip could induce an orientation change in the dipolar rotor but not in the nonpolar analog (for a recent example of an azimuthal molecular rotor controlled by the STM tip, see Reference 58). Compound 3 can exist as three pairs of helical enantiomers because of the propeller-like conformation of the tetra-arylcyclobutadienes. For at least one out of the three diastereomers, an asymmetric potential energy surface can be predicted by molecular dynamics simulations on application of an alternating electric field.55... 

When the H- H dipole-dipole interaction can be measured for a specific pair of H nuclei, studies of the temperature dependence of both the H NMR line-shape and the H NMR relaxation provide a powerful way of probing the molecular dynamics, even in very low temperature regimes at which the dynamics often exhibit quantum tunnelling behaviour. In such cases, H NMR can be superior to quasielastic neutron scattering experiments in terms of both practicality and resolution. The experimental analysis can be made even more informative by carrying out H NMR measurements on single crystal samples. In principle, studies of both the H NMR lineshape and relaxation properties can be used to derive correlation times (rc) for the motion in practice, however, spin-lattice relaxation time (T measurements are more often used to measure rc as they are sensitive to the effects of motion over considerably wider temperature ranges. 

At even lower temperatures, some unusual properties of matter are displayed. Consequently, new experimental and theoretical methods are being created to explore and describe chemistry in these regimes. In order to account for zero-point energy effects and tunneling in simulations, Voth and coworkers developed a quantum molecular dynamics method that they applied to dynamics in solid hydrogen. In liquid helium, superfluidity is displayed in He below its lambda point phase transition at 2.17 K. In the superfluid state, helium s thermal conductivity dramatically increases to 1000 times that of copper, and its bulk viscosity drops effectively to zero. Apkarian and coworkers have recently demonstrated the disappearance of viscosity in superfluid helium on a molecular scale by monitoring the damped oscillations of a 10 A bubble as a function of temperature. These unique properties make superfluid helium an interesting host for chemical dynamics. 

---