Abstract:To set up the general framework for relativistic explicitly correlated wave function methods, the electron-electron coalescence conditions are derived for the wave functions of the Dirac-Coulomb (DC), Dirac-Coulomb-Gaunt (DCG), Dirac-Coulomb-Breit (DCB), modified Dirac-Coulomb (MDC), and zeroth-order regularly approximated (ZORA) Hamiltonians. The manipulations make full use of the internal symmetries of the reduced two-electron Hamiltonians such that the asymptotic behaviors of the wave functions emerge naturally. The results show that, at the coalescence point of two electrons, the wave functions of the DCG Hamiltonian are regular, while those of the DC and DCB Hamiltonians have weak singularities of the type r(12)(ν) with ν being negative and of O(α(2)). The behaviors of the MDC wave functions are related to the original ones in a simple manner, while the spin-free counterparts are somewhat different due to the complicated electron-electron interaction. The behaviors of the ZORA wave functions depend on the chosen potential in the kinetic energy operator. In the case of the nuclear attraction, the behaviors of the ZORA wave functions are very similar to those of the nonrelativistic ones, just with an additional correction of O(α(2)) to the nonrelativistic cusp condition. However, if the Coulomb interaction is also included, the ZORA wave functions become close to the large-large components of the DC wave functions. Note that such asymptotic expansions of the relativistic wave functions are only valid within an extremely small convergence radius R(c) of O(α(2)). Beyond this radius, the behaviors of the relativistic wave functions are still dominated by the nonrelativistic limit, as can be seen in terms of direct perturbation theory (DPT) of relativity. However, as the two limits α → 0 and r(12) → 0 do not commute, DPT is doomed to fail due to incorrect descriptions of the small-small component Ψ(SS) of the DC wave function for r(12) < R(c). Another deduction from the possible divergence of Ψ(SS) at r(12) = R(c) is that the DC Hamiltonian has no bound electronic states, although the last word cannot be said. These findings enrich our understandings of relativistic wave functions. On the practical side, it is shown that, under the no-pair approximation, relativistic explicitly correlated wave function methods can be made completely parallel to the nonrelativistic counterparts, as demonstrated explicitly for MP2-F12. Yet, this can only be achieved by using an extended no-pair projector.

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

Gaussians for Electronic and Rovibrational Quantum Dynamics

Density Matrix Renormalization Group Approach Based on the Coupled-Cluster Downfolded Hamiltonians

Orbital-optimized Density Functional Calculations of Molecular Rydberg Excited States with Real Space Grid Representation and Self-Interaction Correction

Ehrenfest dynamics with localized atomic-orbital basis sets within the projector augmented-wave method

A regularized second-order correlation method from Green's function theory

Relativistic coupled cluster with completely renormalized and perturbative triples corrections

Formulation and Implementation of Frequency-Dependent Linear Response Properties with Relativistic Coupled Cluster Theory for GPU-accelerated Computer Architectures

Relativistic Explicit Correlation: Coalescence Conditions and Practical Suggestions.

Electric-field gradients at the nuclei from all-electron, four-component relativistic density-functional theory using Gaussian-type orbitals

Renormalized Internally-Contracted Multireference Coupled Cluster with Perturbative Triples

A functional renormalization group approach to electronic structure calculations for systems without translational symmetry

Relativistic Linear Response in Quantum-Electrodynamical Density Functional Theory

One-particle operator representation over two-particle basis sets for relativistic QED computations

Accurate variational electronic structure calculations with the density matrix renormalization group

Multidimensional quantum dynamics with explicitly correlated Gaussian wave packets using Rothe's method

Efficient treatment of local meta-generalized gradient density functionals via auxiliary density expansion: the density fitting (DF) J+X approximation

Explicitly correlated electronic structure calculations with transcorrelated matrix product operators

Challenges with relativistic GW calculations in solids and molecules

Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation.

Cholesky Decomposition-Based Implementation of Relativistic Two-Component Coupled-Cluster Methods for Medium-Sized Molecules