Quantum dynamics, complex coordinate

So far, no exact 6D quantum dynamics calculation has been reported for diatomic dissociation on surface. But with modern computer power, such numerical endeavor will undoubtly be realized soon. We also note here recent mixed quantum/classical studies of Jackson who treated three COM coordinates classically and three internal molecular coordinates quantum mechanically for H2/Cu(100) (120). Such treatment seems quite promising for more complex systems. 

A review of the Journal of Physical Chemistry A, volume 110, issues 6 and 7, reveals that computational chemistry plays a major or supporting role in the majority of papers. Computational tools include use of large Gaussian basis sets and density functional theory, molecular mechanics, and molecular dynamics. There were quantum chemistry studies of complex reaction schemes to create detailed reaction potential energy surfaces/maps, molecular mechanics and molecular dynamics studies of larger chemical systems, and conformational analysis studies. Spectroscopic methods included photoelectron spectroscopy, microwave spectroscopy circular dichroism, IR, UV-vis, EPR, ENDOR, and ENDOR induced EPR. The kinetics papers focused on elucidation of complex mechanisms and potential energy reaction coordinate surfaces. 

In contrast to the subsystem representation, the adiabatic basis depends on the environmental coordinates. As such, one obtains a physically intuitive description in terms of classical trajectories along Born-Oppenheimer surfaces. A variety of systems have been studied using QCL dynamics in this basis. These include the reaction rate and the kinetic isotope effect of proton transfer in a polar condensed phase solvent and a cluster [29-33], vibrational energy relaxation of a hydrogen bonded complex in a polar liquid [34], photodissociation of F2 [35], dynamical analysis of vibrational frequency shifts in a Xe fluid [36], and the spin-boson model [37,38], which is of particular importance as exact quantum results are available for comparison. 

However, details of this process including the mode of urea binding, the protonation state of individual surround protein residues, and the exact identity of the nucleophile are still under debate. Cyanate also was proposed as a possible intermediate in the urease mechanism (33). Recent quantum chemical calculations and molecular dynamics simulations indicated that hydrolytic and ehmination mechanisms might indeed compete, and that both are viable reaction channels for urease (34—37). Finally, an important issue is Why does urease require nickel as the metal of choice, whereas most other metallohydrolases use zinc While it was speculated that, inter alia, the relatively rigid and stable coordination environment around the Ni(II) ions as opposed to the higher kinetic lability and lower thermodynamic stability of Zn(II) complexes might play a role (31), this fundamental question has not yet been answered. 

---