Regularized relativistic corrections for polyelectronic and polyatomic systems with explicitly correlated Gaussians.

Drachmann's regularization approach is implemented for floating explicitly correlated Gaussians (fECGs) and molecular systems. Earlier applications of drachmannized relativistic corrections for molecular systems were hindered due to the unknown analytic matrix elements of 1/rix1/rjy-type operators with fECGs. In the present work, one of the 1/r factors is approximated by a linear combination of Gaussians, which results in calculable integrals. The numerical approach is found to be precise and robust over a range of molecular systems and nuclear configurations, and thus, it opens the route toward an automated evaluation of high-precision relativistic corrections over potential energy surfaces of polyatomic systems. Furthermore, the newly developed integration approach makes it possible to construct the matrix representation of the square of the electronic Hamiltonian relevant for energy lower-bound as well as time-dependent computations of molecular systems with a flexible and high-precision fECG basis representation.

One-Particle Operator Representation over Two-Particle Basis Sets for Relativistic QED Computations.

This work is concerned with two-spin-1/2-fermion relativistic quantum mechanics, and it is about the construction of one-particle projectors using an inherently two-particle, "explicitly correlated" basis representation necessary for good numerical convergence of the interaction energy. It is demonstrated that a faithful representation of the one-particle operators, which appear in intermediate but essential computational steps, can be constructed over a many-particle basis set by accounting for the full Hilbert space beyond the physically relevant antisymmetric subspace. Applications of this development can be foreseen for the computation of quantum-electrodynamics corrections for a correlated relativistic reference state and high-precision relativistic computations of medium-to-high-Z helium-like systems, for which other two-particle projection techniques are unreliable.

Bound-State Relativistic Quantum Electrodynamics: A Perspective for Precision Physics with Atoms and Molecules.

Precision physics aims to use atoms and molecules to test and develop the fundamental theory of matter, possibly beyond the Standard Model. Most of the atomic and molecular phenomena are described by the quantum electrodynamics (QED) sector of the Standard Model. Do we have the computational tools, algorithms, and practical equations for the most possible complete computation of atoms and molecules within the QED sector? What is the fundamental equation to start with? Is it still Schrödinger's wave equation for molecular matter, or is there anything beyond that? This paper provides a concise overview of the relativistic QED framework and recent numerical developments targeting precision physics and spectroscopy applications with common features of the robust and successful relativistic quantum chemistry methodology.

QED corrections to the correlated relativistic energy: One-photon processes.

This work is a collection of initial calculations and formal considerations within the Salpeter-Sucher exact equal-time relativistic quantum electrodynamics framework. The calculations are carried out as preparation for the computation of pair, retardation, and radiative corrections to the relativistic energy of correlated two-spin-1/2-fermion systems. In this work, particular attention is paid to the retardation and the "one-loop" self-energy corrections, which are known to be among the largest corrections to the correlated relativistic energy. The theoretical development is supplemented with identifying formal connections to the non-relativistic quantum electrodynamics framework, which is based on a correlated but non-relativistic reference, as well as to the "1/Z approach," which is built on a relativistic but independent-particle zeroth order. The two complementary directions currently provide the theoretical framework for light atomic-molecular precision spectroscopy and heavy-atom phenomena. The present theoretical efforts pave the way for relativistic QED corrections to (explicitly) correlated relativistic computations.

Methane dimer rovibrational states and Raman transition moments.

Benchmark-quality rovibrational data are reported for the methane dimer from variational nuclear motion computations using an ab initio intermolecular potential energy surface reported by [M. P. Metz et al., Phys. Chem. Chem. Phys., 2019, 21, 13504-13525]. A simple polarizability model is used to compute Raman transition moments that may be relevant for future direct observation of the intermolecular dynamics. Non-negligible ΔK ≠ 0 transition moments arise in this symmetric top system due to strong rovibrational couplings.

Vibrational infrared and Raman spectra of HCOOH from variational computations.

All vibrational energies of the formic acid molecule in different forms (trans-, cis-, delocalized-) were converged up to 4500 cm-1 beyond the zero-point vibrational energy with the GENIUSH-Smolyak variational approach and using an ab initio potential energy surface [D. P. Tew and W. Mizukami, J. Phys. Chem. A, 120, 9815-9828 (2016)]. The full-dimensional dipole and polarizability surfaces were fitted to points computed at the CCSD/aug-cc-pVTZ level of theory. Then, body-fixed vibrational dipole and polarizability transition moments were evaluated and used to simulate jet-cooled infrared and Raman spectra of HCOOH. The benchmark-quality vibrational energy, transition moment, and wave function data will be used in further work for comparison with vibrational experiments, and in further rovibrational computations.

The Bethe-Salpeter QED Wave Equation for Bound-State Computations of Atoms and Molecules.

Interactions in atomic and molecular systems are dominated by electromagnetic forces and the theoretical framework must be in the quantum regime. The physical theory for the combination of quantum mechanics and electromagnetism, quantum electrodynamics has been "established" by the mid-twentieth century, primarily as a scattering theory. To describe atoms and molecules, it is important to consider bound states. In the nonrelativistic quantum mechanics framework, bound states can be efficiently computed using robust and general methodologies with systematic approximations developed for solving wave equations. With the sight of the development of a computational quantum electrodynamics framework for atomic and molecular matter, the field theoretic Bethe-Salpeter wave equation expressed in space-time coordinates, its exact equal-time variant, and emergence of a relativistic wave equation, is reviewed. A computational framework, with initial applications and future challenges in relation with precision spectroscopy, is also highlighted.

Relativistic two-electron atomic and molecular energies using LS coupling and double groups: Role of the triplet contributions to singlet states.

The triplet contribution is computed to the 1 and 2 S0e1 states of the He atom, to the 1S0e1 state of the Li+ and Be2+ ions, and to the X1Σg + ground state of the H2 molecule by extensive use of double-group symmetry (equivalent to LS coupling for the atomic systems) during the course of the variational solution of the no-pair Dirac-Coulomb-Breit (DCB) wave equation. The no-pair DCB energies are converged within sub-parts-per-billion relative precision, using an explicitly correlated Gaussian basis optimized to the non-relativistic energies. The α fine-structure constant dependence of the triplet sector contribution to the variational energy is α4Eh at leading order, in agreement with the formal perturbation theory result available from the literature.

Evaluation of the Bethe Logarithm: From Atom to Chemical Reaction.

A general computational scheme for the (nonrelativistic) Bethe logarithm is developed, opening the route to "routine" evaluation of the leading-order quantum electrodynamics correction (QED) relevant for spectroscopic applications for small polyatomic and polyelectronic molecular systems. The implementation relies on the Schwartz method and minimization of a Hylleraas functional. In relation with electronically excited states, a projection technique is considered, which ensures positive definiteness of the functional over the entire parameter (photon momentum) range. Using this implementation, the Bethe logarithm is converged to a relative precision better than 1:103 for selected electronic states of the two-electron H2 and H3+, and the three-electron He2+ and H+H2 molecular systems. The present work focuses on nuclear configurations near the local minimum of the potential energy surface, but the computations can be repeated also for other structures.

Exact quantum dynamics developments for floppy molecular systems and complexes.

Molecular rotation, vibration, internal rotation, isomerization, tunneling, intermolecular dynamics of weakly and strongly interacting systems, intra-to-inter-molecular energy transfer, hindered rotation and hindered translation over surfaces are important types of molecular motions. Their fundamentally correct and detailed description can be obtained by solving the nuclear Schrödinger equation on a potential energy surface. Many of the chemically interesting processes involve quantum nuclear motions which are 'delocalized' over multiple potential energy wells. These 'large-amplitude' motions in addition to the high dimensionality of the vibrational problem represent challenges to the current (ro)vibrational methodology. A review of the quantum nuclear motion methodology is provided, current bottlenecks of solving the nuclear Schrödinger equation are identified, and solution strategies are reviewed. Technical details, computational results, and analysis of these results in terms of limiting models and spectroscopically relevant concepts are highlighted for selected numerical examples.

High-dimensional neural network potentials for accurate vibrational frequencies: the formic acid dimer benchmark.

In recent years, machine learning potentials (MLP) for atomistic simulations have attracted a lot of attention in chemistry and materials science. Many new approaches have been developed with the primary aim to transfer the accuracy of electronic structure calculations to large condensed systems containing thousands of atoms. In spite of these advances, the reliability of modern MLPs in reproducing the subtle details of the multi-dimensional potential-energy surface is still difficult to assess for such systems. On the other hand, moderately sized systems enabling the application of tools for thorough and systematic quality-control are nowadays rarely investigated. In this work we use benchmark-quality harmonic and anharmonic vibrational frequencies as a sensitive probe for the validation of high-dimensional neural network potentials. For the case of the formic acid dimer, a frequently studied model system for which stringent spectroscopic data became recently available, we show that high-quality frequencies can be obtained from state-of-the-art calculations in excellent agreement with coupled cluster theory and experimental data.

Variational vs perturbative relativistic energies for small and light atomic and molecular systems.

Variational and perturbative relativistic energies are computed and compared for two-electron atoms and molecules with low nuclear charge numbers. In general, good agreement of the two approaches is observed. Remaining deviations can be attributed to higher-order relativistic, also called non-radiative quantum electrodynamics (QED), corrections of the perturbative approach that are automatically included in the variational solution of the no-pair Dirac-Coulomb-Breit (DCB) equation to all orders of the α fine-structure constant. The analysis of the polynomial α dependence of the DCB energy makes it possible to determine the leading-order relativistic correction to the non-relativistic energy to high precision without regularization. Contributions from the Breit-Pauli Hamiltonian, for which expectation values converge slowly due the singular terms, are implicitly included in the variational procedure. The α dependence of the no-pair DCB energy shows that the higher-order (α4Eh) non-radiative QED correction is 5% of the leading-order (α3Eh) non-radiative QED correction for Z = 2 (He), but it is 40% already for Z = 4 (Be2+), which indicates that resummation provided by the variational procedure is important already for intermediate nuclear charge numbers.

On the Breit interaction in an explicitly correlated variational Dirac-Coulomb framework.

The Breit interaction is implemented in the no-pair variational Dirac-Coulomb (DC) framework using an explicitly correlated Gaussian basis reported in the previous paper [P. Jeszenszki, D. Ferenc, and E. Mátyus, J. Chem. Phys. 156, 084111 (2022)]. Both a perturbative and a fully variational inclusion of the Breit term are considered. The no-pair DC plus perturbative Breit and the no-pair DC-Breit energies are compared with perturbation theory results including the Breit-Pauli Hamiltonian and leading-order non-radiative quantum electrodynamics corrections for low Z values. Possible reasons for the observed deviations are discussed.

Variational Dirac-Coulomb explicitly correlated computations for atoms and molecules.

The Dirac-Coulomb equation with positive-energy projection is solved using explicitly correlated Gaussian functions. The algorithm and computational procedure aims for a parts-per-billion convergence of the energy to provide a starting point for further comparison and further developments in relation with high-resolution atomic and molecular spectroscopy. Besides a detailed discussion of the implementation of the fundamental spinor structure, permutation, and point-group symmetries, various options for the positive-energy projection procedure are presented. The no-pair Dirac-Coulomb energy converged to a parts-per-billion precision is compared with perturbative results for atomic and molecular systems with small nuclear charge numbers. Paper II [D. Ferenc, P. Jeszenszki, and E. Mátyus, J. Chem. Phys. 156, 084110 (2022).] describes the implementation of the Breit interaction in this framework.

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.

Performance of a black-box-type rovibrational method in comparison with a tailor-made approach: Case study for the methane-water dimer.

The present work intends to join and respond to the excellent and thoroughly documented rovibrational study of X. G. Wang and T. Carrington, Jr. [J. Chem. Phys. 154, 124112 (2021)] that used an approach tailored for floppy dimers with an analytic dimer Hamiltonian and a non-product basis set including Wigner D functions. It is shown in the present work that the GENIUSH black-box-type rovibrational method can approach the performance of the tailor-made computation for the example of the floppy methane-water dimer. Rovibrational transition energies and intensities are obtained in the black-box-type computation with a twice as large basis set and in excellent numerical agreement in comparison with the more efficient tailor-made approach.

All-order explicitly correlated relativistic computations for atoms and molecules.

A variational solution procedure is reported for the many-particle no-pair Dirac-Coulomb and Dirac-Coulomb-Breit Hamiltonians aiming at a parts-per-billion (ppb) convergence of the atomic and molecular energies, described within the fixed nuclei approximation. The procedure is tested for nuclear charge numbers from Z = 1 (hydrogen) to 28 (iron). Already for the lowest Z values, a significant difference is observed from leading-order Foldy-Woythusen perturbation theory, but the observed deviations are smaller than the estimated self-energy and vacuum polarization corrections.

Fingerprint region of the formic acid dimer: variational vibrational computations in curvilinear coordinates.

Curvilinear kinetic energy models are developed for variational nuclear motion computations including the inter- and the low-frequency intra-molecular degrees of freedom of the formic acid dimer. The coupling of the inter- and intra-molecular modes is studied by solving the vibrational Schrödinger equation for a series of vibrational models, from two up to ten active vibrational degrees of freedom by selecting various combinations of active modes and constrained coordinate values. Vibrational states, nodal assignment, and infrared vibrational intensity information is computed using the full-dimensional potential energy surface (PES) and electric dipole moment surface developed by Qu and Bowman [Phys. Chem. Chem. Phys., 2016, 18, 24835; J. Chem. Phys., 2018, 148, 241713]. Good results are obtained for several fundamental and combination bands in comparison with jet-cooled vibrational spectroscopy experiments, but the description of the ν8 and ν9 fundamental vibrations, which are close in energy and have the same symmetry, appears to be problematic. For further progress in comparison with experiment, the potential energy surface, and in particular, its multi-dimensional couplings representation, requires further improvement.

Orientational decoherence within molecules and emergence of the molecular shape.

The question of classicality is addressed in relation with the shape of the nuclear skeleton of molecular systems. As the most natural environment, the electrons of the molecule are considered as continuously monitoring agents for the nuclei. For this picture, an elementary formalism of decoherence theory is developed and numerical results are presented for few-particle systems. The numerical examples suggest that the electron-nucleus Coulomb interaction is sufficient for inducing a blurred shape with strong quantum coherences in compounds of the lightest elements, H2, D2, T2, and HeH+.

Nonadiabatic, Relativistic, and Leading-Order QED Corrections for Rovibrational Intervals of

The rovibrational intervals of the ^{4}He_{2}^{+} molecular ion in its X ^{2}Σ_{u}^{+} ground electronic state are computed by including the nonadiabatic, relativistic, and leading-order quantum-electrodynamics corrections. Good agreement of theory and experiment is observed for the rotational excitation series of the vibrational ground state and the fundamental vibration. The lowest-energy rotational interval is computed to be 70.937 69(10) cm^{-1} in agreement with the most recently reported experimental value, 70.937 589(23)(60)_{sys} cm^{-1} [L. Semeria et al., Phys. Rev. Lett. 124, 213001 (2020)PRLTAO0031-900710.1103/PhysRevLett.124.213001].