Generalized elimination of the global translation from explicitly correlated Gaussian functions

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.

Deformed explicitly correlated Gaussians

Deformed explicitly correlated Gaussian (DECG) basis functions are introduced, and their matrix elements are calculated. All matrix elements can be calculated analytically in a closed form, except the Coulomb one, which has to be approximated by a Gaussian expansion. The DECG basis functions can be used to solve problems with nonspherical potentials. One example of such potential is the dipole self-interaction term in the Pauli-Fierz Hamiltonian. Examples are presented showing the accuracy and necessity of deformed Gaussian basis functions to accurately solve light-matter coupled systems in cavity QED.

Lower Bounds for Nonrelativistic Atomic Energies

A recently developed lower bound theory for Coulombic problems (E. Pollak, R. Martinazzo, J. Chem. Theory Comput. 2021, 17, 1535) is further developed and applied to the highly accurate calculation of the ground-state energy of two- (He, Li⁺, and H^-) and three- (Li) electron atoms. The method has been implemented with explicitly correlated many-particle basis sets of Gaussian type, on the basis of the highly accurate (Ritz) upper bounds they can provide with relatively small numbers of functions. The use of explicitly correlated Gaussians is developed further for computing the variances, and the necessary modifications are here discussed. The computed lower bounds are of submilli-Hartree (parts per million relative) precision and for Li represent the best lower bounds ever obtained. Although not yet as accurate as the corresponding (Ritz) upper bounds, the computed bounds are orders of magnitude tighter than those obtained with other lower bound methods, thereby demonstrating that the proposed method is viable for lower bound calculations in quantum chemistry applications. Among several aspects, the optimization of the wave function is shown to play a key role for both the optimal solution of the lower bound problem and the internal check of the theory.

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₃.

H3 + as a five-body problem described with explicitly correlated Gaussian basis sets

Various explicitly correlated Gaussian (ECG) basis sets are considered for the solution of the molecular Schrödinger equation with particular attention to the simplest polyatomic system, H₃ ⁺. Shortcomings and advantages are discussed for plain ECGs, ECGs with the global vector representation, floating ECGs and their numerical projection, and ECGs with complex parameters. The discussion is accompanied with particle density plots to visualize the observations. In order to be able to use large complex ECG basis sets in molecular calculations, a numerically stable algorithm is developed, the efficiency of which is demonstrated for the lowest rotationally and vibrationally excited states of H₂ and H₃ ⁺.

Effective non-adiabatic Hamiltonians for the quantum nuclear motion over coupled electronic states

The quantum mechanical motion of the atomic nuclei is considered over a single- or a multidimensional subspace of electronic states which is separated by a gap from the rest of the electronic spectrum over the relevant range of nuclear configurations. The electron-nucleus Hamiltonian is block-diagonalized up to O(εⁿ⁺¹) through a unitary transformation of the electronic subspace, and the corresponding nth-order effective Hamiltonian is derived for the quantum nuclear motion. Explicit but general formulas are given for the second- and the third-order corrections. As a special case, the second-order Hamiltonian corresponding to an isolated electronic state is recovered which contains the coordinate-dependent mass-correction terms in the nuclear kinetic energy operator. For a multidimensional, explicitly coupled electronic band, the second-order Hamiltonian contains the usual Born-Oppenheimer terms and nonadiabatic corrections, but generalized mass-correction terms appear as well. These, earlier neglected terms, perturbatively account for the outlying (discrete and continuous) electronic states not included in the explicitly coupled electronic subspace.

Non-adiabatic mass-correction functions and rovibrational states of 4 He 2 + ( X 2 Σ u + )

The mass-correction functions in the second-order non-adiabatic Hamiltonian are computed for the ⁴ He 2 + molecular ion using the variational method, floating explicitly correlated Gaussian functions, and a general coordinate-transformation formalism. When non-adiabatic rovibrational energy levels are computed using these (coordinate-dependent) mass-correction functions and a highly accurate potential energy and diagonal Born-Oppenheimer correction curve, significantly improved theoretical results are obtained for the nine rotational and two rovibrational intervals known from high-resolution spectroscopy experiments.

Explicitly correlated Gaussian functions with shifted-center and projection techniques in pre-Born-Oppenheimer calculations

Numerical projection methods are elaborated for the calculation of eigenstates of the non-relativistic many-particle Coulomb Hamiltonian with selected rotational and parity quantum numbers employing shifted explicitly correlated Gaussian functions, which are, in general, not eigenfunctions of the total angular momentum and parity operators. The increased computational cost of numerically projecting the basis functions onto the irreducible representations of the three dimensional rotation-inversion group is the price to pay for the increased flexibility of the basis functions. This increased flexibility allowed us to achieve a substantial improvement for the variational upper bound to the Pauli-allowed ground-state energy of the H 3 + = { p + , p + , p + , e - , e - } molecular ion treated as an explicit five-particle system. We compare our pre-Born-Oppenheimer result obtained for this molecular ion with rotational-vibrational calculations carried out on a potential energy surface.

Fundamental vibration frequency and rotational structure of the first excited vibrational level of the molecular helium ion ( He 2 + )

The term values of the rotational levels of the first excited vibrational state of the electronic ground state of He 2 + with a rotational quantum number N ⁺ ≤ 13 have been determined with an accuracy of 1.2 × 10^-3 cm^-1 (∼35 MHz) by multichannel-quantum-defect-theory-assisted Rydberg spectroscopy of metastable He₂. Comparison of the experimental term values with the most accurate ab initio results for He 2 + available in the literature [W.-C. Tung, M. Pavanello, and L. Adamowicz, J. Chem. Phys. 136, 104309 (2012)] reveals inconsistencies between the theoretical and experimental results that increase with increasing rotational quantum numbers. The fundamental vibrational wavenumber of He 2 + was determined to be 1628.3832(12) cm^-1 by fitting effective molecular constants to the obtained term values.