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     

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     

Adapted Gaussian basis sets for atoms Cs to Lr based on the generator coordinate Hartree–Fock method

We have applied a discretized version of the generator coordinate Hartree–Fock method to generate adapted Gaussian basis sets for atoms Cs (Z=55) to Lr (Z=103). Our Hartree–Fock total energy results, for all atoms studied, are better than the corresponding Hartree–Fock energy results attained with previous Gaussian basis sets. For the atoms Cs to Lr we have obtained an energy value within the accuracy of 10⁻⁴ to 10⁻³ hartree when compared with the corresponding numerical Hartree–Fock total energy results. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 858–865, 1998     

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.     

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     

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     

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     

Accurate atomic Gaussian basis functions for second-row atoms: Small split-valence 3-21SP and 4-22SP basis sets

Small split-valence Gaussian 3-21SP and 4-22SP basis sets, previously reported for the first-row atoms [Chem. Phys. Lett., 229 , 151 (1996)], have been extended for the second-row elements of the Periodic Table. The total energies of the ground states of the second-row atoms calculated with the new basis sets are significantly lower than those obtained with the well-known 3-21G (J. Am. Chem. Soc., 104 , 2797 (1982)] and 4-31G [J. Chem. Phys., 56 , 5255 (1972)] basis sets. This is because, as first noted in our previous work for first-row atoms, that the 3-21G and 4-31G basis sets only correspond to a local minimum of the Hartree–Fock energy functional, which is relatively far from its global minimum. The proposed basis sets have been tested by performing geometry optimizations and calculations of normal frequencies in the harmonic approximation of some diatomic and polyatomic molecules at the Hartree–Fock level. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1200–1210     

Towards a complete basis set limit of Hartree–Fock method: correlation-consistent versus polarized-consistent basis sets

In this paper the convergence pattern of correlation-consistent (cc-pVxZ) and polarized-consistent (PC-n) hierarchies relative to the complete basis set limit have been considered in a small set of diatomic molecules. Using the sequence of these basis sets it was demonstrated that potential energy surfaces derived from basis-set-dependent solution of the Hartree–Fock equations achieves the exact numerical derived potential energy surfaces (PESs) in an ordered manner. So it was possible to compute the spectroscopic parameters in the complete basis set limit with considerable accuracy using the most extended members of both hierarchies. On the other hand, for the first time the detailed convergence patterns of total energies in three separate inter-nuclear distances have been considered in these molecules and it was demonstrated that the total energies arrive at microhartree accuracy at a considerable rate. Possible performance of extrapolation schemes is discussed and it was demonstrated that reliable extrapolation procedures indeed exist. A successful test of the proposed extrapolation method, using the three most extended members of polarized-consistent basis sets, has been accomplished on selected polyatomic molecules.     

Quantum supercharger library: Hyper‐parallelism of the Hartree–Fock method

We present here a set of algorithms that completely rewrites the Hartree–Fock (HF) computations common to many legacy electronic structure packages (such as GAMESS‐US, GAMESS‐UK, and NWChem) into a massively parallel compute scheme that takes advantage of hardware accelerators such as Graphical Processing Units (GPUs). The HF compute algorithm is core to a library of routines that we name the Quantum Supercharger Library (QSL). We briefly evaluate the QSL's performance and report that it accelerates a HF 6‐31G Self‐Consistent Field (SCF) computation by up to 20 times for medium sized molecules (such as a buckyball) when compared with mature Central Processing Unit algorithms available in the legacy codes in regular use by researchers. It achieves this acceleration by massive parallelization of the one‐ and two‐electron integrals and optimization of the SCF and Direct Inversion in the Iterative Subspace routines through the use of GPU linear algebra libraries. © 2015 Wiley Periodicals, Inc.     

Use of noninteger n‐Slater type orbitals in combined Hartree–Fock–Roothaan theory for calculation of isoelectronic series of atoms Be to Ne

The ground state calculations in the combined Hartree–Fock–Roothaan approach are performed for the neutral and the first 20 cationic members of the isoelectronic series of atoms from Be to Ne using noninteger n‐Slater type orbitals. For the total energies obtained, only a small deviation has been found. At the same time, the size of the present noninteger n‐Slater type orbitals is smaller than that of the usual extended integer n‐Slater functions in literature. All of the nonlinear parameters are fully optimized. The relationship between optimized parameters and atomic number Z is also investigated. For each atom, the total energies are given in tables. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Local Hartree–Fock orbitals using a three‐level optimization strategy for the energy

Using the three‐level energy optimization procedure combined with a refined version of the least‐change strategy for the orbitals—where an explicit localization is performed at the valence basis level—it is shown how to more efficiently determine a set of local Hartree–Fock orbitals. Further, a core–valence separation of the least‐change occupied orbital space is introduced. Numerical results comparing valence basis localized orbitals and canonical molecular orbitals as starting guesses for the full basis localization are presented. The results show that the localization of the occupied orbitals may be performed at a small computational cost if valence basis localized orbitals are used as a starting guess. For the unoccupied space, about half the number of iterations are required if valence localized orbitals are used as a starting guess compared to a canonical set of unoccupied Hartree–Fock orbitals. Different local minima may be obtained when different starting guesses are used. However, the different minima all correspond to orbitals with approximately the same locality. © 2013 Wiley Periodicals, Inc.     

Complete one‐ and two‐center partitioning scheme for the total energy in the Hartree–Fock theory

We proposed a complete calculation scheme for attributing the total energy by the Hartree–Fock theory to atoms (E_A) and the region between two atoms (E_AB). It was pointed out that the conventional method using the Fock matrix includes a large amount of mutual contamination in both E_A and E_AB. The new scheme was derived from the basic expression of the total energy. Calculated results by the new scheme satisfy the theoretical requirements. The scaling effect on partitioned energies was also examined. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 71: 35–46, 1999     

Solution of periodic Poisson's equation and the Hartree–Fock approach for solids with extended electron states: Application to linear augmented plane wave method

We formulate a Hartree–Fock‐LAPW method for electronic band structure calculations. The method is based on the Hartree–Fock–Roothaan approach for solids with extended electron states and closed core shells where the basis functions of itinerant electrons are linear augmented plane waves. All interactions within the restricted Hartree–Fock approach are analyzed and in principle can be taken into account. In particular, we obtained the matrix elements for the exchange interactions of extended states and the crystal electric field effects. To calculate the matrix elements of exchange for extended states, we first introduce an auxiliary potential and then integrate it with an effective charge density corresponding to the electron exchange transition under consideration. The problem of finding the auxiliary potential is solved by using the strategy of the full potential LAPW approach, which is based on the general solution of periodic Poisson's equation. Here, we use an original technique for the general solution of periodic Poisson's equation and multipole expansions of electron densities. We apply the technique to obtain periodic potentials of the face‐centered cubic lattice and discuss its accuracy and convergence in comparison with other methods. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

An instability condition for the Hartree–Fock solution of the infinite one‐dimensional system with two‐crossing bands. I. Singlet‐instability check of metallic carbon nanotube

An instability condition is derived for the Hartree–Fock solution so that it can be applied to the system in which the highest occupied and the lowest unoccupied bands cross at the in‐between point in the Brillouin zone. The instability check developed here is further applied to a metallic single‐walled carbon nanotube having the two‐crossing bands toward prediction of its instability. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 574–582, 2000