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     

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.     

Relativistic correlating basis sets for the sixth-period d-block atoms from Lu to Hg

Contracted Gaussian-type function sets to describe valence correlation are developed for the sixth-period d-block atoms Lu through Hg. A segmented contraction scheme is employed for their compactness and efficiency. Contraction coefficients and exponents are determined by minimizing the deviation from accurate natural orbitals generated from configuration interaction calculations, in which relativistic effects are incorporated through the third-order Douglas-Kroll approximation. The present basis sets yield more than 99% of atomic correlation energies predicted by accurate natural orbital sets of the same size. Relativistic model core potential calculations with the present correlating sets give the spectroscopic constants of the AuH molecule in excellent agreement with experimental results.     

Correspondence between GTO and STO molecular basis sets

We present a method for the characterization of the distance between two spaces: one generated by a Gaussian basis set, and another by a Slater basis set. The method is an extension of one previously developed for atoms that has been modified to cover molecular problems. The current version enables us to obtain Slater basis sets capable of reproducing the results (multielectronic wave functions and orbitals) obtained with Gaussian basis sets. The interest of this result arises from the fact that we will be able to profit from the effort invested in the optimization of high‐quality Gaussian basis sets. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1655–1665, 2001     

Relativistic Gaussian basis sets for molecular calculations: Fully optimized single‐family exponent basis sets for H?Hg

Relativistic single‐family exponent Gaussian basis sets for molecular calculations are presented for the 80 atoms ₁H through ₈₀Hg. The exponent parameters shared by Gaussian basis functions of all symmetry species are fully optimized. Two nucleus models of uniformly charged sphere and Gaussian charge distribution are considered and two kinds of basis sets are generated accordingly. The total energy errors are less than 2 mhartree in any atoms. Some of the present basis sets include small variational collapse (or prolapse), but test calculations show that they could be reliably applied to molecular calculations. © 2005 Wiley Periodicals, Inc. J Comput Chem 27: 48–52, 2006     

Spatially restricted Double Z-Simplified Box Orbital basis sets: Optimization and comparison with some standard basis sets

As part of previous studies, we introduced a new type of basis function named Simplified Box Orbitals (SBOs) that belong to a class of spatially restricted functions which allow the zero differential overlap (ZDO) approximation to be applied with complete accuracy. The original SBOs and their Gaussian expansions SBO-3G form a minimal basis set, which was compared to the standard Slater-type orbital basis set (STO-3G). In the present paper, we have developed the SBO basis functions at double-zeta (DZ) level, and we have assessed the option of expanding the SBO-DZ as a combination of Gaussian functions. Finally, we have determined the quality of the new basis set by comparing the molecular properties calculated with SBO-nG with those achieved with some standard basis sets.     

Relativistic segmented contraction basis sets with core-valence correlation effects for atoms 57La through 71Lu: Sapporo-DK-nZP sets (n?=?D, T, Q)

For the 15 lanthanide atoms ₅₇La through ₇₁Lu, we report Sapporo-DK-nZP sets (n?=?D, T, Q), which are natural extensions of the Sapporo-(DK)-nZP sets for lighter atoms and efficiently incorporate the correlation among electrons in the N through P shells as well as the relativistic effect. The present sets well describe the correlation among the 4s and 4p electrons, which are important in the excitation of 4f electrons. Atomic test calculations of ₅₇La, ₅₈Ce, ₅₉Pr, and ₆₀Nd at configuration interaction with the Davidson correction level of theory confirm high performance of the present basis sets. Molecular test calculations are carried out for ₅₇LaF and ₇₀YbF diatomics at the coupled-cluster level of theory. The calculated spectroscopic constants approach smoothly to the experimental values as the quality of the basis set increases.     

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 relativistic adapted gaussian basis sets for francium through ununoctium without variational prolapse and to be used with both uniform sphere and gaussian nucleus models

Accurate relativistic adapted Gaussian basis sets (RAGBSs) for ₈₇Fr up to ₁₁₈Uuo atoms without variational prolapse were developed here with the use of a polynomial version of the Generator Coordinate Dirac‐Fock method. Two finite nuclear models have been used, the Gaussian and uniform sphere models. The largest RAGBS error, with respect to numerical Dirac‐Fock results, is 15.4 miliHartree for Ununoctium with a basis set size of 33s30p19d14f functions. © 2013 Wiley Periodicals, Inc.     

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     

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.     

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.     

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