A Non-Born–Oppenheimer Self-consistent Field Method

This communication outlines the development of a novel method for the ab initio computation of molecular systems wherein the Born–Oppenheimer approximation is not invoked. In the current method, a common Hamiltonian is employed to operate on the electrons and nuclei simultaneously. In addition, an inseparable wavefunction is generated to describe the non-Born–Oppenheimer behaviour. The physical implications of the new method are discussed.     

Quantum wave packet ab initio molecular dynamics: an approach to study quantum dynamics in large systems

A methodology to efficiently conduct simultaneous dynamics of electrons and nuclei is presented. The approach involves quantum wave packet dynamics using an accurate banded, sparse and Toeplitz representation for the discrete free propagator, in conjunction with ab initio molecular dynamics treatment of the electronic and classical nuclear degree of freedom. The latter may be achieved either by using atom-centered density-matrix propagation or by using Born-Oppenheimer dynamics. The two components of the methodology, namely, quantum dynamics and ab initio molecular dynamics, are harnessed together using a time-dependent self-consistent field-like coupling procedure. The quantum wave packet dynamics is made computationally robust by using adaptive grids to achieve optimized sampling. One notable feature of the approach is that important quantum dynamical effects including zero-point effects, tunneling, as well as over-barrier reflections are treated accurately. The electronic degrees of freedom are simultaneously handled at accurate levels of density functional theory, including hybrid or gradient corrected approximations. Benchmark calculations are provided for proton transfer systems and the dynamics results are compared with exact calculations to determine the accuracy of the approach.     

Ab initio quantum mechanical charge field (QMCF) molecular dynamics: a QM/MM – MD procedure for accurate simulations of ions and complexes

A new formalism for quantum mechanical / molecular mechanical (QM/MM) dynamics of chemical species in solution has been developed, which does not require the construction of any other potential functions except those for solvent–solvent interactions, maintains all the advantages of large simulation boxes and ensures the accuracy of ab initio quantum mechanics for all forces acting in the chemically most relevant region. Interactions between solute and more distant solvent molecules are incorporated by a dynamically adjusted force field corresponding to the actual molecular configuration of the simulated system and charges derived from the electron distribution in the solvate. The new formalism has been tested with some examples of hydrated ions, for which accurate conventional ab initio QM/MM simulations have been previously performed, and the comparison shows equivalence and in some aspects superiority of the new method. As this simulation procedure does not require any tedious construction of two-and three-body interaction potentials inherent to conventional QM/MM approaches, it opens the straightforward access to ab initio molecular dynamics simulations of any kind of solutes, such as metal complexes and other composite species in solution.     

The chemist’s concept of molecular structure

The concept of molecular structure is fundamental to the practice and understanding of chemistry, but the meaning of this term has evolved and is still evolving. The Born–Oppenheimer separation of electronic and nuclear motions lies at the heart of most modern quantum chemical models of molecular structure. While this separation introduces a great computational and practical simplification, it is neither essential to the conceptual formulation of molecular structure nor universally valid. Going beyond the Born–Oppenheimer approximation introduces new paradigms, bringing fresh insight into the chemistry of fluxional molecules, proteins, superconductors and macroscopic dielectrics, thus opening up new avenues for exploration. But it requires that our ideas of molecular structure need to evolve beyond simple ball-and-stick-type models.

Atoms in molecules: beyond Born–Oppenheimer paradigm

This contribution presents the first atoms in molecules study that goes beyond the Born–Oppenheimer paradigm employing the newly developed two-component quantum theory of atoms in molecules (TC-QTAIM). The LiH, LiD, and LiT systems containing quantum instead of clamped hydrogen nuclei are used as typical examples. The computational analysis that is done on non-adiabatic wavefunctions derived from the fully variational multicomponent molecular orbital approach (FV-MC-MO) results in hydrogen atomic basins without any clamped nucleus. The topological analysis of the Γ-field, the field that replaces the usual one-electron density used in the orthodox topological analysis, reveals delicate differences among the considered systems. The calculation of basin properties also demonstrates that the TC-QTAIM differentiates among atomic basins containing isotopes. Since the nuclear dynamics is contained intrinsically in non-adiabatic wavefunctions, the nuclear contribution to both topological analysis and basin properties naturally emerges from the TC-QTAIM analysis resolving the long-standing obstacle of consistent incorporation of nuclear dynamics within the context of the orthodox QTAIM. Also, a similar analysis is done on non-adiabatic wavefunctions describing excited instead of ground nuclear vibrations of the considered systems demonstrating the fact that TC-QTAIM is capable of being employed for both ground and excited nuclear vibrational states.     

Incompleteness of atomic and molecular eigenfunctions

It is well known that the set of bound state hydrogenic eigenfunctions is not complete. Here we show that the set of bound state eigenfunctions for all atomic and molecular systems is incomplete in the Born–Oppenheimer approximation. A trivial argument that there are no positive eigenvalues is also presented which results from the virial theorem. This revised version was published online in July 2006 with corrections to the Cover Date.     

Theoretical study of vitamin properties from combined QM-MM methods: Comparison of chemical shifts and energy

The combination of quantum mechanics (QM) and molecular mechanics (MM) methods has become an alternative tool for many applications for which pure QM and MM are not suitable. The QM-MM method has been used for different types of problems, for example, structural biology, surface phenomena, and the liquid phase. In this paper, we have implemented these methods for vitamins, an important kind of biological molecule, and then compared results. The calculations were done by the full ab initio method (HF/3–21 g and HF/6–31 g) and QM-MM (ONIOM) method with HF(3–21 g)/AM1/UFF; then, we found that the geometry obtained by the QM-MM method is very accurate and this rapid method can be used in place of time consuming ab initio methods for large molecules. A comparison of energy values in the QM-MM and QM methods is given. We compare chemical shifts and conclude that the QM-MM method is a perturbed full QM method. The text was submitted by the authors in English.     

Molecular structure and Bader's theory

 It is shown that a supposed catastrophe of Bader's theory of atoms in molecules, suggested by Cassam-chena? and Jayatilaka [Theor Chem Acc (2001) 105: 213] is merely a consequence of the approximate character of the adiabatic Born–Oppenheimer theory of molecular structure, and that nonadiabatic approaches could be in accordance with Bader's ideas. Received: 4 April 2001 / Accepted: 5 September 2001 / Published online: 3 June 2002     

Non-Markovian dynamics and entanglement in quantum Brownian motion

Dynamical aspects of quantum Brownian motion in a low temperature environment are investigated. We give a systematic calculation of quantum entanglement among two Brownian oscillators without invoking Born–Markov approximation widely used for the study of open systems. Our approach is suitable to probe short time dynamics at cold temperatures where many experiments on quantum information processing are performed.     

Theoretical study on isotope and temperature effect in hydronium ion using ab initio path integral simulation

We have applied the ab initio path integral molecular dynamics simulation to study hydronium ion and its isotopes, which are the simplest systems for hydrated proton and deuteron. In this simulation, all the rotational and vibrational degrees of freedom are treated fully quantum mechanically, while the potential energies of the respective atomic configurations are calculated "on the fly" using ab initio quantum chemical approach. With the careful treatment of the ab initio electronic structure calculation by relevant choices in electron correlation level and basis set, this scheme is theoretically quite rigorous except for Born-Oppenheimer approximation. This accurate calculation allows a close insight into the structural shifts for the isotopes of hydronium ion by taking account of both quantum mechanical and thermal effects. In fact, the calculation is shown to be successful to quantitatively extract the geometrical isotope effect with respect to the Walden inversion. It is also shown that this leads to the isotope effect on the electronic structure as well as the thermochemical properties.     

Ab initio centroid path integral molecular dynamics: application to vibrational dynamics of diatomic molecular systems

An ab initio centroid molecular dynamics (CMD) method is developed by combining the CMD method with the ab initio molecular orbital method. The ab initio CMD method is applied to vibrational dynamics of diatomic molecules, H2 and HF. For the H2 molecule, the temperature dependence of the peak frequency of the vibrational spectral density is investigated. The results are compared with those obtained by the ab initio classical molecular dynamics method and exact quantum mechanical treatment. It is shown that the vibrational frequency obtained from the ab initio CMD approaches the exact first excitation frequency as the temperature lowers. For the HF molecule, the position autocorrelation function is also analyzed in detail. The present CMD method is shown to well reproduce the exact quantum result for the information on the vibrational properties of the system.     

Quantum wavepacket ab initio molecular dynamics: an approach for computing dynamically averaged vibrational spectra including critical nuclear quantum effects

We have introduced a computational methodology to study vibrational spectroscopy in clusters inclusive of critical nuclear quantum effects. This approach is based on the recently developed quantum wavepacket ab initio molecular dynamics method that combines quantum wavepacket dynamics with ab initio molecular dynamics. The computational efficiency of the dynamical procedure is drastically improved (by several orders of magnitude) through the utilization of wavelet-based techniques combined with the previously introduced time-dependent deterministic sampling procedure measure to achieve stable, picosecond length, quantum-classical dynamics of electrons and nuclei in clusters. The dynamical information is employed to construct a novel cumulative flux/velocity correlation function, where the wavepacket flux from the quantized particle is combined with classical nuclear velocities to obtain the vibrational density of states. The approach is demonstrated by computing the vibrational density of states of [Cl-H-Cl]-, inclusive of critical quantum nuclear effects, and our results are in good agreement with experiment. A general hierarchical procedure is also provided, based on electronic structure harmonic frequencies, classical ab initio molecular dynamics, computation of nuclear quantum-mechanical eigenstates, and employing quantum wavepacket ab initio dynamics to understand vibrational spectroscopy in hydrogen-bonded clusters that display large degrees of anharmonicities.     

A heuristic potential function for intramolecular and intrasupermolecular chemical reactions

A heuristic fitting procedure to obtain an analytical potential function for describing a reactive potential energy surface in the neighborhood of the intrinsic reaction coordinate (IRC ) has been developed. For discussion, the pairwise potential function form, ∑a_nr^?n, is assumed in order to fit ab initio quantum mechanical calculations of intramolecular (or intrasupermolecular) interaction energies and its use is found advantageous because all the calculation can be carried out by the linear least squares method. Normal modes perpendicular to IRC are utilized to prepare an initial data base for the potential fitting in the neighborhood of IRC . Some trial molecular dynamics (MD ) simulations are performed in order to check the fitted potential function and, unless they lead to reasonable energies within the tolerance assumed, their results are utilized to construct an improved data base (the dynamic sampling). The present systematic optimization procedure has been applied to the proton transfer reaction of the formamidine–water (FW ) system. The normal mode analysis in both the transition state (TS ) and the stable state (SS ) regions suggests that the present fitted potential function can reproduce satisfactorily the Born–Oppenheimer (BO ) adiabatic surface obtained by ab initio molecular orbital (MO ) calculations. We conclude that our procedure works well for the chemical reaction molecular dynamics (CRMD ) simulation.     

Theoretical advances in the dissolution studies of mineral–water interfaces

This article presents theoretical advances in computational modeling of dissolution at mineral–water interfaces with specific emphasis on silicates. Two different Monte Carlo methods have been developed that target equilibrium properties and kinetics in silicate–water dissolution. The equilibrium properties are explored using the combined reactive Monte Carlo and configurational bias Monte Carlo (RxMC-CBMC) method. The new RxMC-CBMC method is designed to affordably simulate the three-dimensional structure of the mineral with explicit water molecules. The kinetics of the overall dissolution process is studied using a stochastic kinetic Monte Carlo method that utilizes rate constants obtained from accurate ab initio calculations. Both these methods provide important complementary perspective of the complex dynamics involving chemical and physical interactions at the mineral–water interface. The results are compared to experimental and previous computational data available in the literature.     

Ab initio excited-state dynamics of the photoactive yellow protein chromophore

The photoisomerization mechanism of the neutral form of the photoactive yellow protein (PYP) chromophore is investigated using ab initio quantum chemistry and first-principles nonadiabatic molecular dynamics (ab initio multiple spawning or AIMS). We identify the nature of the two lowest-lying excited states, characterize the short-time behavior of molecules excited directly to S2, and explain the origin of the experimentally observed wavelength-dependent photoisomerization quantum yield.     

On the unimolecular elimination of gas-phase iodine monobromide following excitation in the visible wavelength region

Dissociation dynamics of iodine monobromide has been theoretically investigated at numerous excitation wavelengths in the visible wavelength region (450–540 nm) by an accurate time-dependent dynamical simulation technique. The interatomic potential energy curves are derived from ab initio methods which are then employed in a Landau–Zener treatment to calculate the spin-orbit branching ratios in the two halogen products. The extinction coefficients, in arbitrary units, for the discrete electron transition are obtained through a numerical algorithm. Based on these elastic scattering calculations, the results of Br*: Br branching ratios are determined at several simulation wavelengths. Also the electronic properties are evaluated to elucidate the mechanism of chemical dynamics of the system. It is found that the branching dynamics can be clarified via the Landau–Zener model, which treats an avoided crossing between two states of a given nature and evaluates a nonadiabatic effect relevant to the dissociation of IBr molecule.     

Ab initio potential-energy surfaces for complex, multichannel systems using modified novelty sampling and feedforward neural networks

A neural network/trajectory approach is presented for the development of accurate potential-energy hypersurfaces that can be utilized to conduct ab initio molecular dynamics (AIMD) and Monte Carlo studies of gas-phase chemical reactions, nanometric cutting, and nanotribology, and of a variety of mechanical properties of importance in potential microelectromechanical systems applications. The method is sufficiently robust that it can be applied to a wide range of polyatomic systems. The overall method integrates ab initio electronic structure calculations with importance sampling techniques that permit the critical regions of configuration space to be determined. The computed ab initio energies and gradients are then accurately interpolated using neural networks (NN) rather than arbitrary parametrized analytical functional forms, moving interpolation or least-squares methods. The sampling method involves a tight integration of molecular dynamics calculations with neural networks that employ early stopping and regularization procedures to improve network performance and test for convergence. The procedure can be initiated using an empirical potential surface or direct dynamics. The accuracy and interpolation power of the method has been tested for two cases, the global potential surface for vinyl bromide undergoing unimolecular decomposition via four different reaction channels and nanometric cutting of silicon. The results show that the sampling methods permit the important regions of configuration space to be easily and rapidly identified, that convergence of the NN fit to the ab initio electronic structure database can be easily monitored, and that the interpolation accuracy of the NN fits is excellent, even for systems involving five atoms or more. The method permits a substantial computational speed and accuracy advantage over existing methods, is robust, and relatively easy to implement.     

Charge-Tagged DNA Radicals in the Gas Phase Characterized by UV/Vis Photodissociation Action Spectroscopy

Adenosine radicals tagged with a fixed-charge group were generated in the gas phase and structurally characterized by tandem mass spectrometry, deuterium labeling, and UV/Vis action spectroscopy. Experimental results in combination with Born–Oppenheimer molecular dynamics, ab initio, and excited-state calculations led to unambiguous assignment of adenosine radicals as N-7 hydrogen atom adducts. The charge-tagged radicals were found to be electronically equivalent to natural DNA nucleoside radicals.