An accurate non-Born–Oppenheimer calculation of the first purely vibrational transition in LiH molecule

Gaussians for Electronic and Rovibrational Quantum Dynamics

The assumptions underpinning the adiabatic Born-Oppenheimer (BO) approximation are broken for molecules interacting with attosecond laser pulses, which generate complicated coupled electronic-nuclear wave packets that generally will have components of electronic and dissociation continua as well as bound-state contributions. The conceptually most straightforward way to overcome this challenge is to treat the electronic and nuclear degrees of freedom on equal quantum-mechanical footing by not invoking the BO approximation at all. Explicitly correlated Gaussian (ECG) basis functions have proved successful for non-BO calculations of stationary molecular states and energies, reproducing rovibrational absorption spectra with very high accuracy. In this Article, we present a proof-of-principle study of the ability of fully flexible ECGs (FFECGs) to capture the intricate electronic and rovibrational dynamics generated by short, high-intensity laser pulses. By fitting linear combinations of FFECGs to accurate wave function histories obtained on a large real-space grid for a regularized 2D model of the hydrogen atom and for the 2D Morse potential, we demonstrate that FFECGs provide a very compact description of laser-driven electronic and rovibrational dynamics.

Nuclear-electronic all-particle density matrix renormalization group

We introduce the Nuclear-Electronic All-Particle Density Matrix Renormalization Group (NEAP-DMRG) method for solving the time-independent Schrödinger equation simultaneously for electrons and other quantum species. In contrast to the already existing multicomponent approaches, in this work, we construct from the outset a multi-reference trial wave function with stochastically optimized non-orthogonal Gaussian orbitals. By iterative refining of the Gaussians' positions and widths, we obtain a compact multi-reference expansion for the multicomponent wave function. We extend the DMRG algorithm to multicomponent wave functions to take into account inter- and intra-species correlation effects. The efficient parameterization of the total wave function as a matrix product state allows NEAP-DMRG to accurately approximate the full configuration interaction energies of molecular systems with more than three nuclei and 12 particles in total, which is currently a major challenge for other multicomponent approaches. We present the NEAP-DMRG results for two few-body systems, i.e., H₂ and H₃ ⁺, and one larger system, namely, BH₃.

On the physical interpretation of the nuclear molecular orbital energy

Recently, several groups have extended and implemented molecular orbital (MO) schemes to simultaneously obtain wave functions for electrons and selected nuclei. Many of these schemes employ an extended Hartree-Fock approach as a first step to find approximate electron-nuclear wave functions and energies. Numerous studies conducted with these extended MO methodologies have explored various effects of quantum nuclei on physical and chemical properties. However, to the best of our knowledge no physical interpretation has been assigned to the nuclear molecular orbital energy (NMOE) resulting after solving extended Hartree-Fock equations. This study confirms that the NMOE is directly related to the molecular electrostatic potential at the position of the nucleus.

Electric transition dipole moment in pre-Born-Oppenheimer molecular structure theory

This paper presents the calculation of the electric transition dipole moment in a pre-Born-Oppenheimer framework. Electrons and nuclei are treated equally in terms of the parametrization of the non-relativistic total wave function, which is written as a linear combination of basis functions constructed from explicitly correlated Gaussian functions and the global vector representation. The integrals of the electric transition dipole moment are derived corresponding to these basis functions in both the length and the velocity representation. The calculations are performed in laboratory-fixed Cartesian coordinates without relying on coordinates which separate the center of mass from the translationally invariant degrees of freedom. The effect of the overall motion is eliminated through translationally invariant integral expressions. The electric transition dipole moment is calculated between two rovibronic levels of the H2 molecule assignable to the lowest rovibrational states of the X (1)Σ(g)(+) and B (1)Σ(u)(+) electronic states in the clamped-nuclei framework. This is the first evaluation of this quantity in a full quantum mechanical treatment without relying on the Born-Oppenheimer approximation.

Accurate ab initio potential energy curve of F2. III. The vibration rotation spectrum

An analytical expression is found for the accurate ab initio potential energy curve of the fluorine molecule that has been determined in the preceding two papers. With it, the vibrational and rotational energy levels of F(2) are calculated using the discrete variable representation. The comparison of this theoretical spectrum with the experimental spectrum, which had been measured earlier using high-resolution electronic spectroscopy, yields a mean absolute deviation of about 5 cm(-1) over the 22 levels. The dissociation energy with respect to the lowest vibrational energy is calculated within 30 cm(-1) of the experimental value of 12 953+/-8 cm(-1). The reported agreement of the theoretical spectrum and dissociation energy with experiment is contingent upon the inclusion of the effects of core-generated electron correlation, spin-orbit coupling, and scalar relativity. The Dunham analysis [Phys. Rev. 41, 721 (1932)] of the spectrum is found to be very accurate. New values are given for the spectroscopic constants.

Lowest excitation energy of 9Be

Variational calculations employing explicitly correlated Gaussian functions and explicitly including the nuclear motion [i.e., without assuming the Born-Oppenheimer (BO) approximation] have been performed to determine the lowest singlet transition energy in the 9Be atom. The non-BO wave functions were used to calculate the alpha2 relativistic corrections (alpha=1/137.035,999,679). With those corrections and with the alpha3 and alpha4 QED corrections determined previously by others, we obtained 54,677.35 cm(-1) for the 3(1)S-->2(1)S transition energy. This result falls within the error bracket for the experimental transition of 54,677.26(10) cm(-1). This is the first time an electronic transition of Be has been calculated from first principles with the experimental accuracy.

Calculations of the ground states of BeH and BeH+ without the Born-Oppenheimer approximation

Non-Born-Oppenheimer variational calculations employing explicitly correlated Gaussian basis functions have been performed for the ground states of the beryllium monohydride molecule (BeH) and its ion (BeH+), as well as for the beryllium atom (Be) and its ion (Be+). An approach based on the analytical energy gradient calculated with respect to the Gaussian exponential parameters was employed. The calculated energies were used to determine the ionization potential of BeH and the dissociation energies of BeH and BeH+. Also, the generated wave functions were used to compute various expectation values, such as the average interparticle distances and the nucleus-nucleus correlation functions.

Darwin and mass-velocity relativistic corrections in non-Born-Oppenheimer variational calculations

The Pauli approach to account for the mass-velocity and Darwin relativistic corrections has been applied to the formalism for quantum mechanical molecular calculations that does not assume the Born-Oppenheimer (BO) approximation regarding separability of the electronic and nuclear motions in molecular systems. The corrections are determined using the first order perturbation theory and are derived for the non-BO wave function of a diatomic system expressed in terms of explicitly correlated Gaussian functions with premultipliers in the form of even powers of the internuclear distance. As a numerical example we used calculations of the transition energies for pure vibrational states of the HD(+) ion.

Non-Born-Oppenheimer variational calculation of the ground-state vibrational spectrum of LiH+

Very accurate, rigorous, variational, non-Born-Oppenheimer (non-BO) calculations have been performed for the fully symmetric, bound states of the LiH(+) ion. These states correspond to the ground and excited vibrational states of LiH(+) in the ground (2)Sigma(+) electronic state. The non-BO wave functions of the states have been expanded in terms of spherical N-particle explicitly correlated Gaussian functions multiplied by even powers of the internuclear distance and 5600 Gaussians were used for each state. The calculations that, to our knowledge, are the most accurate ever performed for a diatomic system with three electrons have yielded six bound states. Average interparticle distances and nucleus-nucleus correlation function plots are presented.

Darwin and mass-velocity relativistic corrections in the non-Born-Oppenheimer calculations of pure vibrational states of H2

The Darwin and mass-velocity relativistic corrections have been calculated for all pure vibrational states of the H2 using the perturbation theory and very accurate variational wave functions obtained without assuming the Born-Oppenheimer (BO) approximation. Expansions in terms of explicitly correlated Gaussians with premultipliers in the form of even powers of the internuclear distance were used for the wave functions. With the inclusion of the two relativistic corrections to the non-BO energies the transition energies for the highest states agree more with the experimental results.

Matrix elements of N-particle explicitly correlated Gaussian basis functions with complex exponential parameters

In this work we present analytical expressions for Hamiltonian matrix elements with spherically symmetric, explicitly correlated Gaussian basis functions with complex exponential parameters for an arbitrary number of particles. The expressions are derived using the formalism of matrix differential calculus. In addition, we present expressions for the energy gradient that includes derivatives of the Hamiltonian integrals with respect to the exponential parameters. The gradient is used in the variational optimization of the parameters. All the expressions are presented in the matrix form suitable for both numerical implementation and theoretical analysis. The energy and gradient formulas have been programmed and used to calculate ground and excited states of the He atom using an approach that does not involve the Born-Oppenheimer approximation.