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.     

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.     

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     

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.     

Relativistic Gaussian basis sets for radon through plutonium

Summary Relativistic Gaussian basis sets of neutral atoms Rn-Pu and ions Th⁺⁴, U⁺³ and Pu⁺³ in the configurations of average energies are presented. The exponent parameters of the basis sets are determined by least-squares fitting to the numerical Dirac-Fock wave functions. The total energies obtained are within 0.155 a.u. of the Dirac-Fock limits and the qualities of the basis sets are between double-zeta and triple-zeta in the valence parts. Using the exponent parameters the Breit interaction energies have been calculated by perturbation theory and the self-consistent field treatment.     

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     

Contracted Gaussian-type basis functions revisited II. Atoms Na through Ar

We report five minimal-type contracted Gaussian-type function (CGTF) basis sets of the second-row atoms, Na – Ar, for molecular applications. Three of the present CGTF sets are revised versions of those given by Huzinaga and co-workers and the other two are newly developed for more accurate calculations. Practical utility and improved reliability of the present basis sets, augmented by polarization functions, are confirmed by test calculations on the P atom and P₂ molecule both at the self-consistent field (SCF) and configuration interaction (CI) levels. Received: 10 February 1997 / Accepted: 23 April 1997     

The expansion of hydrogen states in Gaussian orbitals

The convergence properties of Gaussian orbitals are studied by considering a very simple system, the hydrogen atom. We have variationally optimized even-tempered basis sets containing up to 60 s functions for the ground state and the first excited S state of the hydrogen atom, to an accuracy of 10^–15E_h. In addition, we have freely optimized the exponents in basis sets containing up to 12 Gaussians. We have studied the convergence of the total energy, the kinetic energy, the extent of the atom as measured by r², and the Fermi-contact interaction at the nucleus in these basis sets as well as in basis sets augmented with additional diffuse or steep functions.     

Optimization of Gaussian basis sets for Dirac-Hartree-Fock calculations

We investigate the optimization of Gaussian basis sets for relativistic calculations within the framework of the restricted Dirac-Hartree-Fock (DHF) method for atoms. We compare results for Rn of nonrelativistic and relativistic basis set optimizations with a finite nuclear-size. Optimization of separate sets for each spin-orbit component shows that the basis set demands for the lower j component are greater than for the higher j component. In particular, the p _1/2 set requires almost as many functions as the s _1/2 set. This implies that for the development of basis sets for heavy atoms, the symmetry type for which a given number of functions is selected should be based on j, not on l, as has been the case in most molecular calculations performed to date.     

Choosing GTO basis sets for periodic HF calculations

We investigate numerical linear dependencies of Gaussian-type orbital basis sets employed in the framework of the Hartree-Fock self-consistent field method for periodic structures, which so far have hampered the use of extended basis sets for non-ionic crystals. These linear dependencies occur when diffuse basis functions are included in a basis set in an uncontrolled manner. We use the condition number of the overlap matrix to lead us in the construction of extended basis sets for periodic structures which avoid numerical linear dependencies. Extended basis sets of high quality are optimized for a number of periodic structures (fcc He, α-Be, α-BN, and B1 NaF) with respect to the energy of the constituent atoms or ions. The results obtained with our basis sets, which do not require reoptimization in the crystal environment, compare favorably with those obtained with other extended basis sets reported in the literature. Received: 20 July 1998/Accepted: 21 August 1998 / Published online: 19 October 1998     

Gaussian basis sets for molecular applications

The choice of basis set in quantum chemical calculations can have a huge impact on the quality of the results, especially for correlated ab initio methods. This article provides an overview of the development of Gaussian basis sets for molecular calculations, with a focus on four popular families of modern atom‐centered, energy‐optimized bases: atomic natural orbital, correlation consistent, polarization consistent, and def2. The terminology used for describing basis sets is briefly covered, along with an overview of the auxiliary basis sets used in a number of integral approximation techniques and an outlook on possible future directions of basis set design. © 2012 Wiley Periodicals, Inc.     

Gaussian basis sets for transition metals of the second series

A Gaussian basis set consisting of (15s, 9p, 8d) Gaussian functions has been optimized for the transition metal atoms of the second series (fourth-row atoms).     

Polarization functions for the modified m6‐31G basis sets for atoms Ga through Kr

The 2df polarization functions for the modified m6‐31G basis sets of the third‐row atoms Ga through Kr (Int J Quantum Chem, 2007, 107, 3028; Int J. Quantum Chem, 2009, 109, 1158) are proposed. The performances of the m6‐31G, m6‐31G(d,p), and m6‐31G(2df,p) basis sets were examined in molecular calculations carried out by the density functional theory (DFT) method with B3LYP hybrid functional, Møller‐Plesset perturbation theory of the second order (MP2), quadratic configuration interaction method with single and double substitutions and were compared with those for the known 6‐31G basis sets as well as with the other similar 641 and 6‐311G basis sets with and without polarization functions. Obtained results have shown that the performances of the m6‐31G, m6‐31G(d,p), and m6‐31G(2df,p) basis sets are better in comparison with the performances of the known 6‐31G, 6‐31G(d,p) and 6‐31G(2df,p) basis sets. These improvements are mainly reached due to better approximations of different electrons belonging to the different atomic shells in the modified basis sets. Applicability of the modified basis sets in thermochemical calculations is also discussed. © 2013 Wiley Periodicals, Inc.     

Reduced-size polarized basis sets for calculations of molecular electric properties. IV. First-row transition metals

Recent studies of the perturbation-dependent basis sets have indicated the possibility of a significant reduction of the size of the usual CGTO sets without considerable loss of accuracy in calculations of molecular electric properties. The resulting (ZPolX) basis sets have been developed for several atoms of the first and second row of the Periodic Table. The same method of the ZPolX basis set generation is extended for the first-row transition metals and the corresponding contracted ZPolX basis sets of the size [6s5p3d1f] are determined for both nonrelativistic and scalar relativistic calculations. The performance of the ZPolX basis sets is verified in calculations on the first-row transition metal oxides at the level of the ROHF, ROHF/CASPT2, and ROHF/CCSD(T) approximations. Also the study of the dipole polarizability of TiCl₄ confirms the excellent features of these very compact basis sets. The ZPolX basis sets for nonrelativistic and relativistic calculations of molecular electric properties are available on the web page http://www.chem.uni.torun.pl/zchk/basis-sets.html.     

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.