Improved generator coordinate Hartree–Fock method applied to generate Gaussian basis sets for the isoelectronic series of the atoms He to Ne

We have generated Gaussian basis sets (GBSs) for the neutral and the first 20 cations members of the isoelectronic series of each ground state atom from He to Ne with the improved generator coordinate Hartree–Fock (IGCHF) method. For all atomic species studied here, our total energy errors are smaller than those calculated by the original GCHF method using GBSs of the same sizes. The largest difference between our total energy results and those computed with a numerical Hartree–Fock approach is equal to 215 μhartree for Co¹⁷⁺. We also compare the ionization potentials obtained with the IGCHF with the corresponding experimental values. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem 88: 252–262, 2002     

An improved generator coordinate Hartree–Fock method applied to the choice of contracted Gaussian basis sets for first‐row diatomic molecules

Accurate Gaussian basis sets (18s for Li and Be and 20s11p for the atoms from B to Ne) for the first‐row atoms, generated with an improved generator coordinate Hartree–Fock method, were contracted and enriched with polarization functions. These basis sets were tested for B₂, C₂, BeO, CN⁻, LiF, N₂, CO, BF, NO⁺, O₂, and F₂. At the Hartree–Fock (HP), second‐order Møller–Plesset (MP2), fourth‐order Møller–Plesset (MP4), and density functional theory (DFT) levels, the dipole moments, bond lengths, and harmonic vibrational frequencies were studied, and at the MP2, MP4, and DFT levels, the dissociation energies were evaluated and compared with the corresponding experimental values and with values obtained using other contracted Gaussian basis sets and numerical HF calculations. For all diatomic molecules studied, the differences between our total energies, obtained with the largest contracted basis set [6s5p3d1f], and those calculated with the numerical HF methods were always less than 3.2 mhartree. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 15–23, 2000     

Accurate Gaussian basis sets for the H2O molecule generated with the molecular improved generator coordinate Hartree–Fock method

The molecular improved generator coordinate Hartree–Fock (MIGCHF) method is used to generate accurate basis sets of primitive Gaussian-type functions for the H₂O molecule. Sequences of increasing size atom centered basis sets are employed to explore the accuracy that can be achieved with this method. Using the O(24s14p8d5f2g1h);H(22s9p5d2f1g) basis set, the HF and second-order electron correlation energies of the H₂O ground state at the experimental geometry are computed as −76.0674680 and −0.3491935 hartree, respectively. The HF energy is in error by 20 μhartree and the second-order correlation energy corresponds to 96.5% of an estimate of the limiting value. The relevance of the present calculations is to show the accuracy that can be achieved in studies of small polyatomic molecules with the MIGCHF method.     

Generator coordinate Hartree–Fock method applied to the choice of a contracted Gaussian basis for the second-row atoms

The generator coordinate Hartree–Fock (GCHF) method is employed as a criterion for the selection of a 18s12p Gaussian basis for the atoms Na–Ar. The role of the weight functions in the assessment of the numerical integration range of the GCHF equations is shown. The extended basis is then contracted to (10s6p) by a standard procedure and in combination with the previously contracted (7s5p) Gaussian basis for the atoms Li–Ne is enriched with polarization functions. This basis is tested for AlF, SiO, PN, BCl, and P₂. The properties of interest were HF total energies, MP2 dipolar moments, bond distances, and dissociation energies. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 927–934, 1997     

The generator coordinate method in the unrestricted Hartree‐Fock formalism

The generator coordinate method was implemented in the unrestricted Hartree‐Fock formalism. Weight functions were built from Gaussian generator functions for 1s, 2s, and 2p orbitals of carbon and oxygen atoms. These weight functions show a similar behavior to those found in the generator coordinate restricted Hartree‐Fock method, i.e., they are smooth, continuous, and tend to zero in the limits of integration. Moreover, the weight functions obtained are different for spin‐up and spin‐down electrons what is a result from spin polarization. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Relativistic Gaussian basis sets for the atoms hydrogen to neon

Gaussian basis sets for use in relativistic molecular calculations are developed for atoms and ions with one to ten electrons. A relativistic radial wavefunction coupled to an angular function of l-symmetry is expanded into a linear combination of spherical Gaussians of the form r ^l exp (–r ²). One set of basis functions is used for all large and small components of the same angular symmetry. The expansion coefficients and the orbital exponents have been determined by minimizing the integral over the weighted square of the deviation between the Dirac or Dirac-Fock radial wavefunctions and their analytical approximations. The basis sets calculated with a weighting function inversely proportional to the radial distance are found to have numerical constants very similar to those of their energy-optimized non-relativistic counterparts. Atomic sets are formed by combining l-subsets. The results of relativistic and non-relativistic calculations based on these sets are analyzed with respect to different criteria, e.g. their ability to reproduce the relativistic total energy contribution and the spin-orbit splitting. Contraction schemes are proposed.Dedicated to Prof. Dr. A. Neckel on occasion of his 60th birthday     

Accurate relativistic adapted Gaussian basis sets for hydrogen through xenon without variational prolapse and to be used with both uniform sphere and Gaussian nucleus models

Accurate relativistic adapted Gaussian basis sets (RAGBSs) from H (Z = 1) through Xe (Z = 54) without variational prolapse have been developed by employing a polynomial version of the Generator Coordinate Dirac‐Fock (p‐GCDF) method. Two nuclear models have been used in this work: (1) the finite nucleus of uniform proton‐charge distribution, and (2) the finite nucleus with a Gaussian proton–charge distribution. The largest errors observed are only 1.5 mHartree (silver and cadmium) and the RAGBS sizes are much smaller than previous accurate relativistic Gaussian basis sets that were shown to be free of variational prolapse. © 2005 Wiley Periodicals, Inc. J Comput Chem 27: 61–71, 2006     

Accurate relativistic adapted Gaussian basis sets for Cesium through Radon without variational prolapse and to be used with both uniform sphere and Gaussian nucleus models

Accurate relativistic adapted Gaussian basis sets (RAGBSs) from Cs (Z = 55) through Rn (Z = 86) without variational prolapse were developed by using the polynomial version of the Generator Coordinate Dirac-Fock method. The RAGBSs presented here can be used with any of two popular finite nucleus models, the uniform sphere and the Gaussian models. The largest RAGBS error is 4.5 mHartree for Radon with a size of 30s27p17d11f.     

A polynomial version of the generator coordinate Dirac-Fock method

A polynomial version of the Generator Coordinate Dirac-Fock (p-GCDF) method is introduced and applied to develop Adapted Gaussian Basis Sets (AGBS) for helium- and beryllium-like atomic species (He, Ne +8, Ar +16, Sn +48, Be, Ne +6, Ar +14, and Sn +46) and for Kr and Xe atoms. The Dirac-Fock-Coulomb and Dirac-Fock-Breit energies obtained with these basis sets are in excellent agreement with numerical finite-difference calculations. Moreover, the sizes of the AGBS generated here with the p-GCDF method are significantly smaller than the size of previous relativistic Gaussian basis sets.     

High-quality Gaussian basis sets for fourth-row atoms

Summary Energy-optimized Gaussian basis sets of triple-zeta quality for the atoms Rb-Xe have been derived. Two series of basis sets are developed; (24s 16p 10d) and (26s 16p 10d) sets which we expand to 13d and 19p functions as the 4d and 5p shells become occupied. For the atoms lighter than Cd, the (24s 16p 10d) sets with triple-zeta valence distributions are higher in energy than the corresponding double-zeta distribution. To ensure a triple-zeta distribution and a global energy minimum the (26s 16p 10d) sets were derived. Total atomic energies from the largest basis sets are between 198 and 284E _H above the numerical Hartree-Fock energies.     

A study on universal Gaussian basis sets for first-row atoms

Analysis of various optimum and non-optimum Gaussian basis sets for firstrow elements have indicated that with a minimum increase of the basis set size and without loss of accuracy of the calculated total energy, a single universal Gaussian basis set may replace individually optimized Gaussian basis sets for a series of atoms. Such a universal Gaussian basis set may substantially reduce the computational work required for the calculation of molecular integrals in ab initio MO calculations.     

Gaussian basis sets for first‐, third‐, and fourth‐row positive and negative ions

An improved generator coordinate Hartree–Fock (HF) method is used to generate accurate triple‐optimized Gaussian basis sets for the cations from He⁺ (Z=2) through Ne⁺ (Z=10) and from K⁺ (Z=19) through Xe⁺ (Z=54), and for the anions from H⁻ (Z=1) through F⁻ (Z=9) and from K⁻ (Z=19) through I⁻ (Z=53). For all ions here studied, our ground‐state HF total energies are better than those calculated with the generator coordinate HF method, using optimized Gaussian basis sets of the same size. For all ions studied, the largest difference between our total energy values and the corresponding results obtained with a numerical HF method is equal to 3.434 mhartrees for Te⁺. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 126–130, 2001     

Gaussian basis sets for calculation of spin densities in first-row atoms

Summary The suitability of Gaussian basis sets for ab initio calculation of Fermi contact spin densities is established by application to the prototype first-row atoms B-F having open shell p electrons. Small multiconfiguration self-consistent-field wave functions are used to describe relevant spin and orbital polarization effects. Basis sets are evaluated by comparing the results to highly precise numerical grid calculations previously carried out with the same wave function models. It is found that modest contracted Gaussian basis sets developed primarily for Hartree-Fock calculations can give semiquantitative results if augmented by diffuse functions and if further uncontracted in the outer core-inner valence region.     

Hartree–Fock wave functions with a modified GTO basis for atoms

We have solved the atomic Hartree–Fock equations by using the algebraic approach, expanding the single-particle radial wave function in terms of a modified Gaussian type orbitals (GTOs) basis. Several atomic properties such as Kato's cusp condition for the electron density or the correct asymptotic behavior of the electron momentum density distribution are accurately verified. Additionally the energy of the atomic ground state can be obtained by using a smaller number of basis functions than in standard GTO expansions. This study has been performed for several atoms of the first three rows. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 65 : 59–64, 1997     

Compatibility of correlation‐consistent basis sets with a hybrid Hartree–Fock/density functional method

The present study examines the feasibility of combining the correlation‐consistent basis sets developed by Dunning and coworkers with the hybrid Hartree–Fock/density functional method B3LYP. Furthermore, extrapolation to the complete basis set (CBS) limit minimizes errors due to the presence of an incomplete basis set and can act as a rigorous test of the limitations of the B3LYP method. Equilibrium geometries, energies, and harmonic vibrational frequencies were determined for a series of well‐studied, yet computationally challenging, small inorganics and their respective ions. The results were then extrapolated to the CBS limit, where applicable, and compared to experiment. It was found that a union between the hybrid Hartree–Fock/density functional B3LYP method and Dunning's augmented correlation‐consistent basis sets gave results that were comparable to molecular orbital methods that explicitly account for electron correlation. Furthermore, the minimum basis set necessary to attain reasonable results for the systems studied was aug‐cc‐pVTZ. Upgrading to the aug‐cc‐pVQZ level and subsequent extrapolation to the CBS limit further improved the overall agreement with the experiment. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 207–216, 1999