Srinivasan S. Iyengar

Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals

We propose and implement an alternative approach to the original Car–Parrinello method where the density matrix elements (instead of the molecular orbitals) are propagated together with the nuclear degrees of freedom. Our new approach has the advantage of leading to an O(N) computational scheme in the large system limit. Our implementation is based on atom-centered Gaussian orbitals, which are especially suited to deal effectively with general molecular systems. The methodology is illustrated by applications to the three-body dissociation of triazine and to the dynamics of a cluster of a chloride ion with 25 water molecules.

Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals. III. Comparison with Born-Oppenheimer dynamics

In a recently developed approach to ab initio molecular dynamics (ADMP), we used an extended Lagrangian to propagate the density matrix in a basis of atom centered Gaussian functions. Results of trajectory calculations obtained by this method are compared with the Born–Oppenheimer approach (BO), in which the density is converged at each step rather than propagated. For NaCl, the vibrational frequency with ADMP is found to be independent of the fictitious electronic mass and to be equal to the BO trajectory result. For the photodissociation of formaldehyde, H2CO→H2+CO, and the three body dissociation of glyoxal, C2H2O2→H2+2CO, very good agreement is found between the Born–Oppenheimer trajectories and the extended Lagrangian approach in terms of the rotational and vibrational energy distributions of the products. A 1.2 ps simulation of the dynamics of chloride ion in a cluster of 25 water molecules was used as a third test case. The Fourier transform of the velocity–velocity autocorrelation function showed ...

Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals. II. Generalizations based on mass-weighting, idempotency, energy conservation and choice of initial conditions

A generalization is presented here for a newly developed approach to ab initio molecular dynamics, where the density matrix is propagated with Gaussian orbitals. Including a tensorial fictitious mass facilitates the use of larger time steps for the dynamics process. A rigorous analysis of energy conservation is presented and used to control the deviation of the fictitious dynamics trajectory from the corresponding Born–Oppenheimer dynamics trajectory. These generalizations are tested for the case of the Cl−(H2O)25 cluster. It is found that, even with hydrogen atoms present in the system, no thermostats are necessary to control the exchange of energy between the nuclear and the fictitious electronic degrees of freedom.

Further analysis of solutions to the time‐independent wave packet equations of quantum dynamics. II. Scattering as a continuous function of energy using finite, discrete approximate Hamiltonians

We consider further how scattering information (the S‐matrix) can be obtained, as a continuous function of energy, by studying wave packet dynamics on a finite grid of restricted size. Solutions are expanded using recursively generated basis functions for calculating Green’s functions and the spectral density operator. These basis functions allow one to construct a general solution to both the standard homogeneous Schrodinger’s equation and the time‐independent wave packet, inhomogeneous Schrodinger equation, in the non‐interacting region (away from the boundaries and the interaction region) from which the scattering solution obeying the desired boundary conditions can be constructed. In addition, we derive new expressions for a ‘‘remainder or error term,’’ which can hopefully be used to optimize the choice of grid points at which the scattering information is evaluated. Problems with reflections at finite boundaries are dealt with using a Hamiltonian which is damped in the boundary region as was done by ...

The properties of ion-water clusters. II. Solvation structures of Na+, Cl−, and H+ clusters as a function of temperature

Ion-water-cluster properties are investigated both through the multistate empirical valence bond potential and a polarizable model. Equilibrium properties of the ion-water clusters H+(H2O)100, Na+(H2O)100, Na+(H2O)20, and Cl-(H2O)17 in the temperature region 100-450 K are explored using a hybrid parallel basin-hopping and tempering algorithm. The effect of the solid-liquid phase transition in both caloric curves and structural distribution functions is investigated. It is found that sodium and chloride ions largely reside on the surface of water clusters below the cluster melting temperature but are solvated into the interior of the cluster above the melting temperature, while the solvated proton was found to have significant propensity to reside on or near the surface in both the liquid- and solid-state clusters.

Hydrogen Tunneling in an Enzyme Active Site: A Quantum Wavepacket Dynamical Perspective

We study the hydrogen tunneling problem in a model system that represents the active site of the biological enzyme, soybean lipoxygenase-1. Toward this, we utilize quantum wavepacket dynamics performed on potential surfaces obtained by using hybrid density functional theory under the influence of a dynamical active site. The kinetic isotope effect is computed by using the transmission amplitude of the wavepacket, and the experimental value is reproduced. By computing the hydrogen nuclear orbitals (eigenstates) along the reaction coordinate, we note that tunneling for both hydrogen and deuterium occurs through the existence of distorted, spherical s-type proton wave functions and p-type polarized proton wave functions for transfer along the donor-acceptor axis. In addition, there is also a significant population transfer through distorted p-type proton wave functions directed perpendicular to the donor-acceptor axis (via intervening pi-type proton eigenstate interactions) which underlines the three-dimensional nature of the tunneling process. The quantum dynamical evolution indicates a significant contribution from tunneling processes both along the donor-acceptor axis and along directions perpendicular to the donor-acceptor axis. Furthermore, the tunneling process is facilitated by the occurrence of curve crossings and avoided crossings along the proton eigenstate adiabats.

Ab initio molecular dynamics: Propagating the density matrix with gaussian orbitals. IV. Formal analysis of the deviations from born-oppenheimer dynamics

In the context of the recently developed Atom-centered Density Matrix Propagation (ADMP) approach to ab initio molecular dynamics, a formal analysis of the deviations from the Born—Oppenheimer surface is conducted. These deviations depend on the fictitious mass and on the magnitude of the commutator of the Fock and density matrices. These quantities are found to be closely interrelated and the choice of the fictitious mass provides a lower bound on the deviations from the Born—Oppenheimer surface. The relations are illustrated with an example calculation for the Cl−(H2O)25 cluster. We also show that there exists a direct one-to-one correspondence between approximate Born—Oppenheimer dynamics, where SCF convergence is restricted by a chosen threshold value for the commutator of the Fock and density matrices, and extended Lagrangian dynamics performed using a finite value for the fictitious mass. The analysis is extended to the nuclear forces used in the ADMP approximation. The forces are shown to be more general than those standardly used in Born—Oppenheimer dynamics, with the addition terms in the nuclear forces depending on the commutator mentioned above.

Insights from first principles molecular dynamics studies toward infrared multiple-photon and single-photon action spectroscopy: Case study of the proton-bound dimethyl ether dimer

We have used finite temperature ab initio molecular dynamics simulations in conjunction with computation of critical quantum nuclear effects to probe the differences between single-photon argon tagged action spectral results and infrared multiple-photon dissociation experiments for a proton bound molecular ion system. We find that the principal difference between the results in these experimental techniques is essentially that of cluster temperature. The multiple-photon dissociation experiments conducted using room temperature ions reflect a larger degree of conformational freedom compared to the colder single-photon argon tagged action spectral results. Our ab initio molecular dynamics simulation techniques accurately capture the effects of conformational sampling, adequately reproduce both spectra, and can be utilized to assign the dynamically averaged finite temperature spectra.

Effect of time-dependent basis functions and their superposition error on atom-centered density matrix propagation (ADMP): connections to wavelet theory of multiresolution analysis.

We present a rigorous analysis of the primitive Gaussian basis sets used in the electronic structure theory. This leads to fundamental connections between Gaussian basis functions and the wavelet theory of multiresolution analysis. We also obtain a general description of basis set superposition error which holds for all localized, orthogonal or nonorthogonal, basis functions. The standard counterpoise correction of quantum chemistry is seen to arise as a special case of this treatment. Computational study of the weakly bound water dimer illustrates that basis set superposition error is much less for basis functions beyond the 6-31+G(*) level of Gaussians when structure, energetics, frequencies, and radial distribution functions are to be calculated. This result will be invaluable in the use of atom-centered Gaussian functions for ab initio molecular dynamics studies using Born-Oppenheimer and atom-centered density matrix propagation.

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.