Elimination of the translational kinetic energy contamination in pre-Born–Oppenheimer calculations

In this paper, we present a simple strategy for the elimination of the translational kinetic energy contamination of the total energy in pre-Born–Oppenheimer calculations carried out in laboratory-fixed Cartesian coordinates (LFCCs). The simple expressions for the coordinates and the operators are thus preserved throughout the calculations, while the mathematical form and the parametrisation of the basis functions are chosen so that the translational and rotational invariances are respected. The basis functions are constructed using explicitly correlated Gaussian functions (ECGs) and the global vector representation. First, we observe that it is not possible to parametrise the ECGs so that the system is at rest in LFCCs and at the same time the basis functions are square integrable with a non-vanishing (positive-definite) norm. Then, we work out a practical strategy to circumvent this problem by making use of the properties of the linear transformation between the LFCCs and translationally invariant and centre-of-mass Cartesian coordinates as well as the transformation properties of the corresponding basis function parameter matrices. By exploiting these formal mathematical relationships, we can identify and separate the translational contamination terms in the matrix representation of the kinetic energy operator in the LFCC formalism. We present numerical examples for the translational contamination and its elimination for the two lowest rotational energy levels of the singlet hydrogen molecule, corresponding to para- and ortho-H2, respectively, treated as four-particle quantum systems.

[1]  Ludwik Adamowicz,et al.  Born-Oppenheimer and non-Born-Oppenheimer, atomic and molecular calculations with explicitly correlated Gaussians. , 2013, Chemical reviews.

[2]  M. Reiher,et al.  Molecular structure calculations: a unified quantum mechanical description of electrons and nuclei using explicitly correlated Gaussian functions and the global vector representation. , 2012, The Journal of chemical physics.

[3]  Shant Shahbazian,et al.  The two-component quantum theory of atoms in molecules (TC-QTAIM): foundations , 2012, Theoretical Chemistry Accounts.

[4]  I. Mayer,et al.  Internal coordinates of quantum-mechanical systems , 2012 .

[5]  J. Hutter,et al.  Extracting elements of molecular structure from the all-particle wave function. , 2011, The Journal of chemical physics.

[6]  J. Hutter,et al.  On the emergence of molecular structure , 2011 .

[7]  Conrad Sanderson,et al.  Armadillo: An Open Source C++ Linear Algebra Library for Fast Prototyping and Computationally Intensive Experiments , 2010 .

[8]  Erwan Faou,et al.  Computing Semiclassical Quantum Dynamics with Hagedorn Wavepackets , 2009, SIAM J. Sci. Comput..

[9]  A. Császár,et al.  Toward black-box-type full- and reduced-dimensional variational (ro)vibrational computations. , 2009, The Journal of chemical physics.

[10]  J. Komasa,et al.  Nonadiabatic corrections to rovibrational levels of H2. , 2008, The Journal of chemical physics.

[11]  A. Chakraborty,et al.  Inclusion of explicit electron-proton correlation in the nuclear-electronic orbital approach using Gaussian-type geminal functions. , 2008, The Journal of chemical physics.

[12]  M. Quack,et al.  Rovibrational analysis of the ν 4,2ν 6 Fermi resonance band of CH35ClF2 by means of a polyad Hamiltonian involving the vibrational levels ν 4, 2ν 6,ν 6+ν 9 and 2ν 9, and comparison with ab initio calculations , 2006 .

[13]  B. Sutcliffe Comment on "elimination of translational and rotational motions in nuclear orbital plus molecular orbital theory" [J. Chem. Phys. 122, 164101 (2005)]. , 2005, The Journal of chemical physics.

[14]  S. Hyodo,et al.  Elimination of translational and rotational motions in nuclear orbital plus molecular orbital theory. , 2005, The Journal of chemical physics.

[15]  Edward F. Valeev,et al.  The electron and nuclear orbitals model: current challenges and future prospects , 2004 .

[16]  L. Adamowicz,et al.  Non-Born–Oppenheimer calculations of atoms and molecules , 2003 .

[17]  B. Kuhn,et al.  A new six-dimensional analytical potential up to chemically significant energies for the electronic ground state of hydrogen peroxide , 1999 .

[18]  K. Varga,et al.  Stochastic Variational Approach to Quantum-Mechanical Few-Body Problems , 1998 .

[19]  K. Varga,et al.  NEW DESCRIPTION OF ORBITAL MOTION WITH ARBITRARY ANGULAR MOMENTA , 1997, nucl-th/9701030.

[20]  D. Kinghorn Integrals and derivatives for correlated Gaussian functions using matrix differential calculus , 1996 .

[21]  Ed Anderson,et al.  LAPACK Users' Guide , 1995 .

[22]  J. Rychlewski,et al.  Many‐electron explicitly correlated Gaussian functions. I. General theory and test results , 1993 .

[23]  G. Hagedorn Semiclassical quantum mechanics , 1980 .

[24]  K. Szalewicz,et al.  High-accuracy Compton profile of molecular hydrogen from explicitly correlated Gaussian wave function , 1979 .

[25]  K. Singer,et al.  The use of Gaussian (exponential quadratic) wave functions in molecular problems - I. General formulae for the evaluation of integrals , 1960, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[26]  S. F. Boys,et al.  The integral formulae for the variational solution of the molecular many-electron wave equation in terms of Gaussian functions with direct electronic correlation , 1960, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[27]  D. Jackson High Resolution Spectroscopy , 1947, Nature.

[28]  Stephen Wilson,et al.  Handbook of molecular physics and quantum chemistry , 2003 .

[29]  Jacek Rychlewski,et al.  Explicitly correlated wave functions in chemistry and physics : theory and applications , 2003 .

[30]  B. T. Sutcliffe,et al.  Chapter 31 Coordinate Systems and Transformations , 2002 .

[31]  H. Nakai Simultaneous determination of nuclear and electronic wave functions without Born–Oppenheimer approximation: Ab initio NO+MO/HF theory , 2002 .

[32]  D. Luckhaus,et al.  Coupling Across Bonds: Ab Initio Calculations for the Anharmonic Vibrational Resonance Dynamics of the Coupled OH and CH Chromophores in Trans Formic Acid HCOOH , 2000 .