Molecular properties from perturbation theory: A Unified treatment of energy derivatives

The definition of a molecular property as a derivative of the electronic energy with respect to one or more applied perturbations is reviewed. The explicit enumeration of terms entering the derivative formulas is performed by considering in turn the various parameter spaces on which the energy and wave function depend. After deriving general expressions for first, second, and third derivatives for different types of perturbation, the parameter spaces involved in MCSCF and CI cases are identified and used to obtain expressions for the first and second derivatives. An example of an MCSCF third derivative is also given. In addition, the various equation systems defining the perturbed wave functions in each order are derived. Some attention is given to the efficient computer implementation of derivative calculations, and the present work is compared with that of other authors.     

Exact solution of multidimensional hyper‐radial Schrödinger equation for many‐electron quantum systems

In quantum theory, solving Schrödinger equation analytically for larger atomic and molecular systems with cluster of electrons and nuclei persists to be a tortuous challenge. Here, we consider, Schrödinger equation in arbitrary N‐dimensional space corresponding to inverse‐power law potential function originating from a multitude of interactions participating in a many‐electron quantum system for exact solution within the framework of Frobenius method via the formulation of an ansatz to the hyper‐radial wave function. Analytical expressions for energy spectra, and hyper‐radial wave functions in terms of known coefficients of inverse‐power potential function, and wave function parameters have been obtained. A generalized two‐term recurrence relation for power series expansion coefficients has been established. © 2016 Wiley Periodicals, Inc.     

Addressing an instability in unrestricted density functional theory direct dynamics simulations

In Density Functional Theory (DFT) direct dynamics simulations with Unrestricted Hartree Fock (UHF) theory, triplet instability often emerges when numerically integrating a classical trajectory. A broken symmetry initial guess for the wave function is often used to obtain the unrestricted DFT potential energy surface (PES), but this is found to be often insufficient for direct dynamics simulations. An algorithm is described for obtaining smooth transitions between the open-shell and the closed-shell regions of the unrestricted PES, and thus stable trajectories, for direct dynamics simulations of dioxetane and its •O CH₂-CH₂ O• singlet diradical. © 2018 Wiley Periodicals, Inc.     

On the correspondence between classical and quantum mechanics in macromolecular systems: Too much classical chaos

In a recent study we found the classical dynamics of a polyethylene (PE) chain to exhibit low dimensional chaos at temperatures as low as a few Kelvin. These results strongly suggest that classical molecular dynamic simulations in polymer systems can grossly overestimate vibrational motion, which consequently results in disordered structures. In contrast, quantum mechanical calculations using Internal Coordinate Quantum Monte Carlo (an improved method with an initial conjecture for the correct wave function) indicate that the quantum ground state for a three-dimensional model PE chain is far more rigid than determined from molecular dynamics (MD) simulations, even at energies as low as a small fraction of the ground state energy. This result casts uncertainty on the reliability of MD estimates of dynamical or structural quantities relevant to the study of some macromolecular systems.     

Full optimization of Jastrow-Slater wave functions with application to the first-row atoms and homonuclear diatomic molecules

We pursue the development and application of the recently introduced linear optimization method for determining the optimal linear and nonlinear parameters of Jastrow-Slater wave functions in a variational Monte Carlo framework. In this approach, the optimal parameters are found iteratively by diagonalizing the Hamiltonian matrix in the space spanned by the wave function and its first-order derivatives, making use of a strong zero-variance principle. We extend the method to optimize the exponents of the basis functions, simultaneously with all the other parameters, namely, the Jastrow, configuration state function, and orbital parameters. We show that the linear optimization method can be thought of as a so-called augmented Hessian approach, which helps explain the robustness of the method and permits us to extend it to minimize a linear combination of the energy and the energy variance. We apply the linear optimization method to obtain the complete ground-state potential energy curve of the C(2) molecule up to the dissociation limit and discuss size consistency and broken spin-symmetry issues in quantum Monte Carlo calculations. We perform calculations for the first-row atoms and homonuclear diatomic molecules with fully optimized Jastrow-Slater wave functions, and we demonstrate that molecular well depths can be obtained with near chemical accuracy quite systematically at the diffusion Monte Carlo level for these systems.     

Multireference self-consistent-field energies without the many-electron wave function through a variational low-rank two-electron reduced-density-matrix method

The variational two-electron reduced-density-matrix (2-RDM) method allows for the computation of accurate ground-state energies and 2-RDMs of atoms and molecules without the explicit construction of an N-electron wave function. While previous work on variational 2-RDM theory has focused on calculating full configuration-interaction energies, this work presents the first application toward approximating multiconfiguration self-consistent-field (MCSCF) energies via low-rank restrictions on the 1- and 2-RDMs. The 2-RDM method with two- or three-particle N-representability conditions reduces the exponential active-space scaling of MCSCF methods to a polynomial scaling. Because the first-order algorithm [Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)] represents each form of the 1- and 2-RDMs by a matrix factorization, the RDMs are readily defined to have a low rank rather than a full rank by setting the matrix factors to be rectangular rather than square. Results for the potential energy surfaces of hydrogen fluoride, water, and the nitrogen molecule show that the low-rank 2-RDM method yields accurate approximations to the MCSCF energies. We also compute the energies along the symmetric stretch of a 20-atom hydrogen chain where traditional MCSCF calculations, requiring more than 17x10(9) determinants in the active space, could not be performed.     

On the computation of excited states with MCSCF methods

We discuss the theoretical and practical problems arising when trying to compute excited states of nonrelativistic electrons in a molecular system, by multiconfiguration (MCSCF) methods. These nonlinear models approximate the linear Schrödinger theory and are a generalization of the well-known Hartree–Fock approach. Due to the MCSCF nonlinearity, a theoretical definition of what should be a MCSCF excited state is not clear at all, contrarily to the ground state case. We compare various definitions used in Quantum Chemistry. We in particular stress that some defects may lead to important computational problems, already observed in Quantum Chemistry (root flipping). We then present a definition of MCSCF excited states based on a solid mathematical ground and compare it with the most used methods. This new definition leads to a completely new algorithm for computing the first excited state, which was proposed and tested in a collaboration with Cancès and Galicher. Numerical results are provided for the simple case of two-electron systems, as an illustration of the possible issues which can arise as consequences of the nonlinearity of the MCSCF method.     

Geometrical energy derivative evaluation with MRCI wave functions

The theory of MCSCF and CI energy derivatives with respect to geometrical variations is briefly reviewed with special attention given to the MCSCF and MRCI energy gradients. A computational procedure is proposed for MRCI energy gradients that does not require the solution to any “coupled-perturbed MCSCF ” equations, it does not require any expensive direct-CI matrix-vector products involving derivative integrals, and it does not require any derivative integrals to be transformed from the AO basis to the MO basis. An additional feature is that it does not require any changes to existing MCSCF gradient evaluation programs in order to compute MRCI gradients. The only difference in the two cases is the exact nature of the data passed to the gradient evaluation program from the previous steps in the computational procedure. The additional effort required to compute the entire MRCI energy gradient vector is approximately that required for one additional iteration of the MRCI diagonalization procedure and for one additional MCSCF iteration. For large scale MRCI wave functions, the MRCI energy gradient evaluation should only require about 10% of the effort of computing the wave function itself. This computational procedure removes a major computational botleneck of potential energy surface evaluation.     

Hybrid quantum/classical path integral approach for simulation of hydrogen transfer reactions in enzymes

A hybrid quantum/classical path integral Monte Carlo (QC-PIMC) method for calculating the quantum free energy barrier for hydrogen transfer reactions in condensed phases is presented. In this approach, the classical potential of mean force along a collective reaction coordinate is calculated using umbrella sampling techniques in conjunction with molecular dynamics trajectories propagated according to a mapping potential. The quantum contribution is determined for each configuration along the classical trajectory with path integral Monte Carlo calculations in which the beads move according to an effective mapping potential. This type of path integral calculation does not utilize the centroid constraint and can lead to more efficient sampling of the relevant region of conformational space than free-particle path integral sampling. The QC-PIMC method is computationally practical for large systems because the path integral sampling for the quantum nuclei is performed separately from the classical molecular dynamics sampling of the entire system. The utility of the QC-PIMC method is illustrated by an application to hydride transfer in the enzyme dihydrofolate reductase. A comparison of this method to the quantized classical path and grid-based methods for this system is presented.     

Energy conserving approximations to the quantum potential: dynamics with linearized quantum force

Solution of the Schrodinger equation within the de Broglie-Bohm formulation is based on propagation of trajectories in the presence of a nonlocal quantum potential. We present a new strategy for defining approximate quantum potentials within a restricted trial function by performing the optimal fit to the log-derivatives of the wave function density. This procedure results in the energy-conserving dynamics for a closed system. For one particular form of the trial function leading to the linear quantum force, the optimization problem is solved analytically in terms of the first and second moments of the weighted trajectory distribution. This approach gives exact time-evolution of a correlated Gaussian wave function in a locally quadratic potential. The method is computationally cheap in many dimensions, conserves total energy and satisfies the criterion on the average quantum force. Expectation values are readily found by summing over trajectory weights. Efficient extraction of the phase-dependent quantities is discussed. We illustrate the efficiency and accuracy of the linear quantum force approximation by examining a one-dimensional scattering problem and by computing the wavepacket reaction probability for the hydrogen exchange reaction and the photodissociation spectrum of ICN in two dimensions.     

Semiclassical quantum unimolecular reaction rate theory revisited

We describe a semiclassical quantum unimolecular reaction rate theory derived from the corresponding classical theory developed by Davis, Gray, Rice and Zhao (DGRZ). The analysis retains the intuitively useful mechanistic distinctions between intramolecular energy transfer and reaction, with the consequence that the semiclassical quantum theory version neglects some interference effects in the reaction dynamics. In the limiting case that intramolecular energy transfer is very fast compared to the rate of reaction we show that the DGRZ representation of the rate constant can be transformed, using the Weyl correspondence between quantum operators and classical variables, to the quantum flux–flux correlation function representation of the rate constant. In the more general case that the rate of intramolecular energy transfer influences the reaction dynamics, the semiclassical representation of the Wigner function for a classical system with both quasiperiodic and chaotic motion is used to obtain the reaction rate constant. Our analysis identifies the quantum analogue of the classical bottleneck to intramolecular energy transfer with the scars of unstable periodic orbits; it leads to a flux–flux correlation function representation of the rate constant for intramolecular energy transfer.     

Photoexcited tautomerization of vinyl alcohol to acetylaldehyde via a conical intersection from contracted Schrödinger theory

The photoexcited tautomerization of vinyl alcohol to acetylaldehyde via a conical intersection is analyzed through the direct calculation of two-electron reduced density matrices (2-RDMs) from solutions to the anti-Hermitian contracted Schr?dinger equation (ACSE). The study utilizes the recent generalization of the ACSE method for the treatment of excited states [G. Gidofalvi and D. A. Mazziotti, Phys. Rev. A, 2009, 80, 022507]. We computed absolute energies of the critical points as well as various intermediate points along the ground- and excited-state potential energy surface of vinyl alcohol and acetylaldehyde. The ACSE, MCSCF, second-order multireference many-body perturbation theory (MRMP2), and various coupled cluster methods all demonstrate the existence of a family of pathways from vinyl alcohol to acetylaldehyde via a conical intersection that are monotonically decreasing in energy. The conical intersection, proposed for the first time in this paper, is both structurally and energetically similar to the ground-state transition state. We observe a relationship between conical intersections and transition states both in this paper and in our previous work on bicyclobutane's ring conical intersection [J. W. Snyder, Jr. and D. A. Mazziotti, J. Chem. Phys., 2011, 135, 024107]. To treat multireference correlation, we seeded the ACSE with an initial 2-RDM guess from a multiconfiguration self-consistent field (MCSCF) calculation. The ACSE recovers more single-reference correlation energy than MRMP2 and more multireference correlation energy than comparable single-reference wave function methods. The 2-RDMs from the ACSE nearly satisfy all necessary N-representability conditions.     

Six-dimensional quantum dynamics of H2 dissociative adsorption on the Pt(211) stepped surface

Results of experimental studies, and theoretical calculations utilizing classical trajectories, have shown that dissociation of H2 on the Pt(211) stepped surface is enhanced at low energies by a molecular trapping mechanism. Because quantum effects can play a large role at the low energies and long lifetimes that characterize molecular trapping, we have undertaken quantum dynamics calculations for this system, the first to treat all molecular degrees of freedom of a gas molecule reacting on a stepped metallic surface. The calculations show that molecular trapping persists in the quantum system, but only at much lower energies than experimentally seen, pointing to possible deficiencies in the potential energy surface. Classical and quasiclassical trajectory calculations on the same potential provide a reasonable picture of reaction overall, but many of the finer details are inaccurate, and certain classical reaction mechanisms are entirely invalid. We conclude that some skepticism should be shown toward any classical study for which long-lived trapping states play a role.     

Performance of dynamically weighted multiconfiguration self-consistent field and spin-orbit coupling calculations of diatomic molecules of Group 14 elements

The efficacy of several multiconfiguration self-consistent field (MCSCF) methods in the subsequent spin-orbit coupling calculations was studied. Three MCSCF schemes to generate molecular orbitals were analyzed: state-specific, state-averaged, and dynamically weighted MCSCF. With Sn(2)(+) as the representative case, we show that the state-specific MCSCF orbitals lead to discontinuities in potential energy curves when avoided crossings of electronic states occur; this problem can be solved using the state-averaged or dynamically weighted MCSCF orbitals. The latter two schemes are found to give similar results when dynamic electron correlation is considered, which we calculated at the level of multiconfigurational quasidegenerate perturbation theory (MCQDPT). We employed the recently developed Douglas-Kroll spin-orbit adapted model core potential, ZFK3-DK3, and the dynamically weighted MCSCF scheme to calculate the spectroscopic constants of the mono-hydrides and compared them to the results obtained using the older set of potentials, MCP-TZP. We also showed that the MCQDPT tends to underestimate the dissociation energies of the hydrides and discussed to what extent coupled-cluster theory can be used to improve results.