Completeness criteria for explicitly correlated Gaussian geminal bases of axial symmetry

The completeness criteria for the basis set of explicitly correlated Gaussian-type geminals adapted to C_∞v symmetry are given. Specifically, we show that any pair function of Σ⁺ symmetry can be expanded in terms of products involving two spherical Gaussian orbitals located on the internuclear axis and a Gaussian correlating factor with a positive exponent. Pair functions corresponding to other irreducible representations of C_∞v can be expressed as linear combinations of products of a σ⁺ function and an angular factor depending on the azimuthal angles. The minimal set of the angular factors needed for completeness is given. These factors are relevant also for other explicitly correlated bases. © 1997 John Wiley & Sons, Inc.     

Analysis of the errors in explicitly correlated electronic structure theory

The explicitly correlated second order M?ller-Plesset (MP2-R12) methods perform well in reproducing the last detail of the correlation cusp, allowing higher accuracy than can be accessed through conventional means. Nevertheless in basis sets that are practical for calculations on larger systems (i.e., around triple- or perhaps quadruple-zeta) MP2-R12 fails to bridge the divide between conventional MP2 and the MP2 basis set limit. In this contribution we analyse the sources of error in MP2-R12 calculations in such basis sets. We conclude that the main source of error is the choice of the correlation factor r12. Sources of error that must be avoided for accurate quantum chemistry include the neglect of exchange commutators and the extended Brillouin condition. The generalized Brillouin condition is found not to lead to significant errors.     

SF-[2]R12: a spin-adapted explicitly correlated method applicable to arbitrary electronic states

The [2](R12) method [M. Torheyden and E. F. Valeev, J. Chem. Phys. 131, 171103 (2009)] is an explicitly correlated perturbative correction that can greatly reduce the basis set error of an arbitrary electronic structure method for which the two-electron density matrix is available. Here we present a spin-adapted variant (denoted as SF-[2](R12)) that is formulated completely in terms of spin-free quantities. A spin-free cumulant decomposition and multi-reference generalized Brillouin condition are used to avoid three-particle reduced density matrix completely. The computational complexity of SF-[2](R12) is proportional to the sixth power of the system size and is comparable to the cost of the single-reference MP2-R12 method. The SF-[2](R12) method is shown to decrease greatly the basis set error of multi-configurational wave functions.     

New correlation factors for explicitly correlated electronic wave functions

We have investigated the correlation factors exp(-zetar12), r12 exp(-zetar12), erfc(zetar12), and r12 erfc(zetar12) in place of the linear-r12 term for use in explicitly correlated electronic-structure methods. The accuracy obtained with all of these correlation factors is significantly greater than that obtained with the plain correlation factor r12. Polarization functions that are more diffuse than those of standard basis sets give even better results. The correlation factor exp(-zetar12) is very close to the optimum correlation factor for helium and outperforms the others.     

Analytic first derivatives for explicitly correlated,multicenter, Gaussian geminals

Variational calculations utilizing the analytic gradient of explicitly correlated Gaussian molecular integrals are presented for the ground state of the hydrogen molecule. Preliminary results serve to motivate the need for general formulas for analytic first derivatives of molecular integrals involving multicenter, explicitly correlated Gaussian geminals with respect to Gaussian exponents and coordinates of the orbital centers. Explicit formulas for analytic first derivatives of Gaussian functions containing correlation factors of the form exp(-βr_ij²) are derived and discussed. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 991–999, 1997     

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.     

Frequency-dependent nonlinear optical properties with explicitly correlated coupled-cluster response theory using the CCSD(R12) model

Response theory up to infinite order is combined with the explicitly correlated coupled-cluster singles and doubles model including linear-r(12) corrections, CCSD(R12). The additional terms introduced by the linear-r(12) contributions, not present in the conventional CCSD calculation, are derived and discussed with respect to the extra costs required for their evaluation. An implementation is presented up to the cubic response function for one-electron perturbations, i.e., up to frequency-dependent second hyperpolarizabilities. As first applications the authors computed the electronic polarizabilities and second hyperpolarizabilities of BH, N(2), and formaldehyde and show that the improvement in the one-electron basis set convergence known from the R12 method for ground state energies is retained for higher-order optical properties. Frequency-dependent results are presented for the second hyperpolarizability of N(2).     

Second order explicitly correlated R12 theory revisited: a second quantization framework for treatment of the operators' partitionings

Second order R12 theory is presented and derived alternatively using the second quantized hole-particle formalism. We have shown that in order to ensure the strong orthogonality between the R12 and the conventional part of the wave function, the explicit use of projection operators can be easily avoided by an appropriate partitioning of the involved operators to parts which are fully describable within the computational orbital basis and complementary parts that involve imaginary orbitals from the complete orbital basis. Various Hamiltonian splittings are discussed and computationally investigated for a set of nine molecules and their atomization energies. If no generalized Brillouin condition is assumed, with all relevant partitionings the one-particle contribution arising in the explicitly correlated part of the first order wave function has to be considered and has a significant role when smaller atomic orbital basis sets are used. The most appropriate Hamiltonian splitting results if one follows the conventional perturbation theory for a general non-Hartree-Fock reference. Then, no couplings between the R12 part and the conventional part arise within the first order wave function. The computationally most favorable splitting when the whole complementary part of the Hamiltonian is treated as a perturbation fails badly. These conclusions also apply to MP2-F12 approaches with different correlation factors.     

Benchmark calculations for He2 and LiH molecules using explicitly correlated Gaussian functions

Explicitly correlated Gaussian (ECG) functions with carefully optimized non-linear parameters are used to calculate the electronic energies of He₂⁺ and LiH at their equilibrium internuclear distances. The obtained variational upper bounds (−4.99464392 and −8.070538 hartree, respectively) are the lowest reported to date. By extrapolating results obtained with various expansion lengths, the estimations of the Born–Oppenheimer limits are made.     

Molecular structure calculations: a unified quantum mechanical description of electrons and nuclei using explicitly correlated Gaussian functions and the global vector representation

We elaborate on the theory for the variational solution of the Schro?dinger equation of small atomic and molecular systems without relying on the Born-Oppenheimer paradigm. The all-particle Schro?dinger equation is solved in a numerical procedure using the variational principle, Cartesian coordinates, parameterized explicitly correlated Gaussian functions with polynomial prefactors, and the global vector representation. As a result, non-relativistic energy levels and wave functions of few-particle systems can be obtained for various angular momentum, parity, and spin quantum numbers. A stochastic variational optimization of the basis function parameters facilitates the calculation of accurate energies and wave functions for the ground and some excited rotational-(vibrational-)electronic states of H(2) (+) and H(2), three bound states of the positronium molecule, Ps(2), and the ground and two excited states of the (7)Li atom.     

Confinement effects on the electronic structure of M‐shell atoms: A study with explicitly correlated wave functions

The ground and some excited states of Na and Mg atoms confined at the center of a spherical box with impenetrable walls are studied. Variational wave functions including dynamic correlations and configuration mixing have been obtained. Level crossings induced by confinement have been analyzed in terms of the energy of the occupied orbitals of the M shell and the weight of the different configurations. Confinement effects on the correlation energy have been studied. The parameterized optimized effective potential and the variational Monte Carlo methods have been employed. A cut off‐factor has been included to account for the hard wall confinement.     

Application of explicitly correlated Gaussian functions for calculations of the ground state of the beryllium atom

The electronic energy of atoms and molecules may be evaluated accurately by the use of wave functions where the interelectronic distances are explicitly present. In particular, explicitly correlated Gaussian-type functions make these types of calculations feasible and computationally tractable even for more extended systems. The resulting multielectron integrals may be reduced to standard one- and two-electron integrals that are readily evaluated. Initial calculations have been made for the Be atom where all four electrons were correlated at the same time. The preliminary results show that accurate results may be obtained. © 1993 John Wiley & Sons, Inc.     

Energy and energy gradient matrix elements with N-particle explicitly correlated complex Gaussian basis functions with L=1

In this work we consider explicitly correlated complex Gaussian basis functions for expanding the wave function of an N-particle system with the L=1 total orbital angular momentum. We derive analytical expressions for various matrix elements with these basis functions including the overlap, kinetic energy, and potential energy (Coulomb interaction) matrix elements, as well as matrix elements of other quantities. The derivatives of the overlap, kinetic, and potential energy integrals with respect to the Gaussian exponential parameters are also derived and used to calculate the energy gradient. All the derivations are performed using the formalism of the matrix differential calculus that facilitates a way of expressing the integrals in an elegant matrix form, which is convenient for the theoretical analysis and the computer implementation. The new method is tested in calculations of two systems: the lowest P state of the beryllium atom and the bound P state of the positronium molecule (with the negative parity). Both calculations yielded new, lowest-to-date, variational upper bounds, while the number of basis functions used was significantly smaller than in previous studies. It was possible to accomplish this due to the use of the analytic energy gradient in the minimization of the variational energy.     

Simple coupled-cluster singles and doubles method with perturbative inclusion of triples and explicitly correlated geminals: The CCSD(T)R12 model

To approach the complete basis set limit of the "gold-standard" coupled-cluster singles and doubles plus perturbative triples [CCSD(T)] method, we extend the recently proposed perturbative explicitly correlated coupled-cluster singles and doubles method, CCSD(2)(R12) [E. F. Valeev, Phys. Chem. Chem. Phys. 8, 106 (2008)], to account for the effect of connected three-electron correlations. The natural choice of the zeroth-order Hamiltonian produces a perturbation expansion with rigorously separable second-order energy corrections due to the explicitly correlated geminals and conventional triple and higher excitations. The resulting CCSD(T)(R12) energy is defined as a sum of the standard CCSD(T) energy and an amplitude-dependent geminal correction. The method is technically very simple: Its implementation requires no modification of the standard CCSD(T) program and the formal cost of the geminal correction is small. We investigate the performance of the open-shell version of the CCSD(T)(R12) method as a possible replacement of the standard complete-basis-set CCSD(T) energies in the high accuracy extrapolated ab initio thermochemistry model of Stanton et al. [J. Chem. Phys. 121, 11599 (2004)]. Correlation contributions to the heat of formation computed with the new method in an aug-cc-pCVXZ basis set have mean absolute basis set errors of 2.8 and 1.0 kJmol when X is T and Q, respectively. The corresponding errors of the standard CCSD(T) method are 9.1, 4.0, and 2.1 kJmol when X=T, Q, and 5. Simple two-point basis set extrapolations of standard CCSD(T) energies perform better than the explicitly correlated method for absolute correlation energies and atomization energies, but no such advantage found when computing heats of formation. A simple Schwenke-type two-point extrapolation of the CCSD(T)(R12)aug-cc-pCVXZ energies with X=T,Q yields the most accurate heats of formation found in this work, in error on average by 0.5 kJmol and at most by 1.7 kJmol.     

Matrix elements for Ĵ2 and Ĵz operators over explicitly correlated Cartesian Gaussian functions

General formalism for evaluation of multiparticle integrals involving J?² and J?_z operators over explicitly correlated Cartesian Gaussian functions is presented. The integrals are expressed in terms of the general overlap integrals. An explicitly correlated Cartesian Gaussian function is a product of spherical orbital Gaussian functions, powers of the Cartesian coordinates of the particle, and exponential Gaussian factors, which depend on interparticular distances. This development is relevant to both adiabatic and nonadiabatic calculations of energy and properties of multiparticle systems. © 1995 John Wiley & Sons, Inc.     

Newton–Raphson optimization of the explicitly correlated Gaussian functions for the ground state of the beryllium atom

Explicitly correlated Gaussian functions have been used in variational calculations on the ground state of the beryllium atom. In such calculations on systems with more electrons, it becomes imminent and essential to develop effective strategies for optimizing the parameters involved in the basis functions. The theory of analytical first and second derivatives of the variational functional with respect to the Gaussian exponents and its computational implementation in conjunction with the Newton–Raphson optimization technique is described. Some numerical results are presented to illustrate the performance of the method. © 1995 John Wiley & Sons, Inc.     

Multicenter and multiparticle integrals for explicitly correlated cartesian gaussian-type functions

An analytical derivation of multicenter and multiparticle integrals for explicitly correlated Cartesian Gaussian-type cluster functions is demonstrated. The evaluation method is based on the application of raising operators that transform spherical cluster Gaussian functions into Cartesian Gaussian functions.     

Coupled-cluster and explicitly correlated perturbation-theory calculations of the uracil anion

A valence-type anion of the canonical tautomer of uracil has been characterized using explicitly correlated second-order Moller-Plesset perturbation theory (RI-MP2-R12) in conjunction with conventional coupled-cluster theory with single, double, and perturbative triple excitations. At this level of electron-correlation treatment and after inclusion of a zero-point vibrational energy correction, determined in the harmonic approximation at the RI-MP2 level of theory, the valence anion is adiabatically stable with respect to the neutral molecule by 40 meV. The anion is characterized by a vertical detachment energy of 0.60 eV. To obtain accurate estimates of the vertical and adiabatic electron binding energies, a scheme was applied in which electronic energy contributions from various levels of theory were added, each of them extrapolated to the corresponding basis-set limit. The MP2 basis-set limits were also evaluated using an explicitly correlated approach, and the results of these calculations are in agreement with the extrapolated values. A remarkable feature of the valence anionic state is that the adiabatic electron binding energy is positive but smaller than the adiabatic electron binding energy of the dipole-bound state.     

Accurately solving the electronic Schrodinger equation of atoms and molecules using explicitly correlated (r12-) multireference configuration interaction. VII. The hydrogen fluoride molecule

We compute the potential-energy curve of the hydrogen fluoride molecule (HF) using a novel variant of the explicitly correlated multireference averaged coupled-pair functional method with a carefully selected basis set and reference space. After correcting for scalar relativistic effects and spin-orbit coupling, the potential is used to compute the dissociation energy, the equilibrium bond distance, the harmonic frequency, the anharmonicity, and the vibrational levels up to the dissociation limit. The errors in the equilibrium geometry constants compare favorably with the most elaborate (single reference) calculations of the literature. Starting at the region of RA/angstroms approximately 2,...,3, where the covalent HF bond begins to break and where single-reference methods become impractical, our potential begins to slightly underestimate the atomic interaction, which is reflected in an estimated error in the well depth of -0.2 kcal/mol.