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.     

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.     

Explicitly correlated Gaussian functions with r factors for calculations of the ground state of the helium atom

Explicitly correlated Gaussian functions with r exp( ? β r ) factors have been used in variational calculations of the ground state of the helium atom. Additional correlation factors in the form of even powers of r _ij were introduced to the Gaussian functions with exponential correlation components by differentiating these functions with respect to the correlation exponent β. The algorithm of this method and its computational implementation is described. A number of calculations were performed for the ground state of helium atom to test the performance of the basis sets comparising different numbers of the Gaussians with the exp( ? β r ) and r exp( ? β r ) correlation factors. The numerical results indicate that including functions with r factors does not lead to improved results, contrary to what was anticipated initially. © 1994 by John Wiley & Sons, Inc.     

Extended Hartree–Fock calculations for the helium ground state

The extended Hartree–Fock (EHF) wave function of an n-electron system is defined (Löwdin, Phys. Rev. 97 , 1509 (1955)) as the best Slater determinant built on one-electron spin orbitals having a complete flexibility and projected onto an appropriate symmetry subspace. The configuration interaction equivalent to such a wavefunction for the ¹S state of a two-electron atom is discussed. It is shown that there is in this case an infinite number of solutions to the variational problem with energies lower than that of the usual Hartree–Fock function, and with spin orbitals satisfying all the extremum conditions. Two procedures for obtaining EHF spin orbitals are presented. An application to the ground state of Helium within a basic set made up of 4(s), 3(p₀), 2(d₀) and 1 (f₀) Slater orbitals has produced 90% of the correlation energy.     

Collision-induced dipole polarizability of helium dimer from explicitly correlated calculations

Large expansions in basis sets of explicitly correlated Gaussian functions and the variation-perturbation technique were used to calculate the static dipole polarizability of the helium dimer at 16 different internuclear separations from 1.0 to 9.0 bohrs. The convergence towards the complete basis set limit was analyzed in order to estimate uncertainties of all the calculated values. The results are significantly more accurate than literature data. Asymptotically correct analytic fits for the trace and anisotropy of collision-induced polarizability were obtained.     

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.     

Variational functions of the padé type—application to the ground state of the helium atom

It is shown that trial functions involving Padé approximants yield satisfactory results for the ground state of the helium atom. In particular, the five-parameter form reproduces the best variational function of the type Ψ = e^–zs?(u) obtained numerically, to a remarkable extent.     

Accurate calculations of intermolecular interaction energies using explicitly correlated wave functions

Explicitly correlated second-order M?ller-Plesset (MP2-F12) calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented and compared with standard MP2 results. The MP2 complete basis set limits are estimated using basis set extrapolation and augmented quadruple-zeta and quintuple-zeta basis sets. Already with augmented double-zeta basis sets the MP2-F12 interaction energies are found to be closer to the complete basis set limits than standard MP2 calculations with augmented quintuple-zeta basis sets. Various possible approximations in the MP2-F12 method are systematically tested. Best results are obtained with localized orbitals and the diagonal MP2-F12/C(D) ansatz. Hybrid approximations, in which some contributions of the auxiliary basis set are neglected and which considerably reduce the computational cost, have a negligible effect on the interaction energies. Also the orbital-invariant fixed-amplitude approximation of Ten-no leads to only slightly less accurate results. Preliminary results for the neon and benzene dimers, obtained with the recently proposed CCSD(T)-F12a approximation, indicate that the CCSD(T) basis set limits can also be very closely approached using augmented triple-zeta basis sets.     

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.     

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.     

Variational calculations in a halfspace for the interaction of a ground and excited state helium atom with its mirror image

Correlated trial wave-functions (depending explicitly upon the inter-electronic separation in the helium atom) subjected to boundary conditions excluding the helium electrons from the solid surface are used in half space, to evaluate in a variational approach the holding potential between an helium atom in its ground and metastable³S state with a perfectly imaging metallic conductor limited by a plane surface. Inter-electronic correlation effects are pointed out, and are seen to be non negligible. The holding potential for the helium³S displays a well of nearly an order of magnitude deeper than for the ground state. These results could be useful to determine the distance of closest approach of thermal atoms, and an approximate wave function for an helium atom interacting with a surface.     

High‐precision Hy–CI variational calculations for the ground state of neutral helium and helium‐like ions

Hylleraas–configuration interaction (Hy–CI) method variational calculations with up to 4648 expansion terms are reported for the ground ¹S state of neutral helium. Convergence arguments are presented to obtain estimates for the exact nonrelativistic energy of this state. The nonrelativistic energy is calculated to be ?2.9037 2437 7034 1195 9829 99 a.u. Comparisons with other calculations and an energy extrapolation give an estimated nonrelativistic energy of ?2.9037 2437 7034 1195 9830(2) a.u., which agrees well with the best previous variational energy, ?2.9037 2437 7034 1195 9829 55 a.u., of Korobov (Phys Rev A 2000, 61, 64503), obtained using the universal (exponential) variational expansion method with complex exponents (Frolov, A. M.; Smith, V. H. Jr. J Phys B Atom Mol Opt Phys 1995, 28, L449). In addition to He, results are also included for the ground ¹S states of H^?, Li⁺, Be⁺⁺, and B⁺³. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

An empirically adjusted Newton–Raphson algorithm for finding local minima on molecular potential energy surfaces

A simple extension of the Newton–Raphson method is proposed that approximately accounts for anharmonicity in bond-stretching coordinates. By modeling each bonded distance in a polyatomic molecule as a Morse oscillator with no anharmonic stretch-stretch or stretch-bend coupling, a multiplicative correction factor to the Newton–Raphson step is derived. Representative examples suggest that the rate of convergence of the proposed scheme is typically faster than that of the standard Newton–Raphson method.     

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.     

Variational Monte Carlo calculations for some cations and anions of the first‐row atoms using explicitly correlated wave functions

The variational Monte Carlo method is applied to calculate ground‐state energies of some cations and anions of the first‐row atoms. Accurate values providing between 80 and 90% of the correlation energy are obtained. Explicitly correlated wave functions including up to 42 variational parameters are used. The nondynamic correlation due to the 2s ? 2p near degeneracy effect is included by using a multideterminant wave function. The variational free parameters have been fixed by minimizing the energy that has shown to be a more convenient functional than the variance of the local energy, which is the most commonly employed method in variational Monte Carlo calculations. The energies obtained improve previous works using similar wave functions. © 2002 Wiley Periodicals, Inc.; DOI 10.1002/qua.10125     

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.     

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.     

The use of explicitly correlated,partially antisymmetric wave functions in atomic and molecular calculations

Trial wave functions, written as the sum of a configuration interaction expansion and an explicitly correlated term which is not antisymmetric, are proposed for use in calculating the electronic properties of atoms and molecules. A variational principle, modified to allow the use for such partially antisymmetric wave functions, is developed. It is shown that the consequences of partial antisymmetry on calculated expectation values can be estimated. The method avoids difficult three-electron integrals which arise in other theories.     

Explicitly correlated Gaussian calculations of the (2)D Rydberg states of the boron atom

Accurate non-relativistic variational calculations are performed for the seven lowest members of the (2)D Rydberg series (1s(2)2s2p(2), and 1s(2)2s(2)nd, n = 3, [ellipsis (horizontal)], 8) of the boron atom. The wave functions of the states are expanded in terms of all-electron explicitly correlated Gaussian basis functions and the effect of the finite nuclear mass is directly included in the calculations allowing for determining the isotopic shifts of the energy levels. The Gaussian basis is optimized independently for each state with the aid of the analytic energy gradient with respect to the Gaussian parameters. The calculations represent the highest accuracy level currently achievable for the considered states. The computed energies are compared with the available experimental data.