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.     

Contracted polarization functions for the atoms helium through neon

Contracted Gaussian-type function sets are proposed for polarization functions of the atoms helium through neon, with the exception of lithium. A segmented contraction scheme is used for its compactness and computational efficiency. The contraction coefficients and orbital exponents are fully optimized to minimize the difference from accurate atomic natural orbitals. The present polarization functions yield more than 99% of atomic correlation energies predicted by accurate natural orbitals of the same size. Received: 16 April 1997 / Accepted: 26 June 1997     

Contracted polarization functions for the atoms Ca, Ga–Kr, Sr, and In–Xe

 Contracted Gaussian-type function sets are proposed for polarization functions of the atoms Ga–Kr and In–Xe. We also report polarization functions for Ca and Sr. A segmented contraction scheme is used for its compactness and computational efficiency. The contraction coefficients and orbital exponents are fully optimized to minimize the deviation from accurate atomic natural orbitals. The present polarization functions yield more than 99% of atomic correlation energies predicted by accurate natural orbitals of the same size. Received: 23 February 2001 / Accepted: 19 April 2001 / Published online: 13 June 2001     

Systematically convergent correlation consistent basis sets for molecular core-valence correlation effects: the third-row atoms gallium through krypton

The family of correlation consistent polarized valence basis sets has been extended in order to account for core-core and core-valence correlation effects within the third-row, main group atoms gallium through krypton. Construction of the basis sets is similar to that of the atoms boron through argon, where either the difference between core-correlated and valence-only correlation energies were calculated via configuration interaction (CISD) computations on the ground electronic states of the atoms (named cc-pCVnZ) or the sets were optimized with respect to the core-valence correlation energy and a small weight of core-core correlation energy (cc-pwCVnZ). Due to the correlation of 3d orbitals, added shells of higher angular momentum exponents compared to the valence sets are necessary to describe the core region. The pattern of added core-correlating functions is (1s1p1d1f) for double-zeta, (2s2p2d2f1g) for triple-zeta, (3s3p3d3f2g1h) for quadruple-zeta, and (4s4p4d4f3g2h1i) for quintuple-zeta. Atomic and molecular results show good convergence to the CBS limit, with the cc-pwCVnZ sets showing improved convergence compared to the cc-pCVnZ ones for molecular core-valence correlation effects. After testing the basis sets on the homonuclear diatomics Ga2-Kr2 with coupled cluster wave functions, it is concluded that a treatment of core-valence correlation effects is essential for high-accuracy ab initio investigations of third-row-containing molecules. Though the basis sets are optimal for 3s3p3d correlation, preliminary atomic and molecular results show the basis sets to be efficient with respect to 3d-only correlation, and these potentially could be used with 3d-only correlation for more qualitative studies on larger species.     

Contracted basis sets for density functional calculations: segmented versus general contraction

The differences between segmented and general contracted basis sets of double and triple zeta quality are analyzed for first and second row elements. Based on coverage of the exponent space and the performance for molecular properties, it is shown that a segmented contraction requires duplication of one primitive function compared to a general contraction for double zeta type basis sets. For triple zeta basis sets, segmentation necessitates either addition of one primitive function and expanding to a quadruple valence space, or addition of two primitive functions. For molecular properties depending on the valence orbitals, such as atomization energies, equilibrium distances, and vibrational frequencies, some of the inner functions describing the core orbitals can be removed without significantly affecting the accuracy. Several of the popular basis sets in common use correspond to such core-pruned basis sets.     

Linear Approximation of Basis Set Orbital Exponents for Isoelectronic Series of Atoms

The orbital exponents of Slater type atomic orbitals (AOs) in isoelectronic series of atoms may be approximated by the linear dependence on the nuclear charge using a technique developed for optimization of AO basis sets in Hartree–Fock–Roothaan calculations. This approach yields the analytical Hartree–Fock wave functions for any ion in the isoelectronic atomic series without optimization of orbital exponents. The approximated linear equations for atomic orbital basis sets of B, C, O, and F in the ground state are presented as an example.     

Molecular states of atoms. II

Orbital energy parameters, previously obtained from atomic valence state energies, are used in calculating approximate wave functions for their orbitals. The radial factors of these wave functions are expressed as linear combinations of three Gaussian type orbitals with selected exponents, the coefficients being determined by normalisation and reproduction of the kinetic energy and interelectron repulsion parameters. Wave functions of universal form are obtained for the non-transition elements up to xenon. Each calculated s orbital wave function (except 1s) has a radial node, as is appropriate if there is a p orbital in the same shell with none.     

Consistent gaussian basis sets of double‐ and triple‐zeta valence with polarization quality of the fifth period for solid‐state calculations

Consistent basis sets of double‐ and triple‐zeta valence with polarization quality for the fifth period have been derived for periodic quantum‐chemical solid‐state calculations with the crystalline‐orbital program CRYSTAL. They are an extension of the pob‐TZVP basis sets, and are based on the full‐relativistic effective core potentials (ECPs) of the Stuttgart/Cologne group and on the def2‐SVP and def2‐TZVP valence basis of the Ahlrichs group. We optimized orbital exponents and contraction coefficients to supply robust and stable self‐consistent field (SCF) convergence for a wide range of different compounds. The computed crystal structures are compared to those obtained with standard basis sets available from the CRYSTAL basis set database. For the applied hybrid density functional PW1PW, the average deviations of calculated lattice constants from experimental references are smaller with pob‐DZVP and pob‐TZVP than with standard basis sets. © 2018 Wiley Periodicals, Inc.     

6‐31G* basis set for third‐row atoms

Medium basis sets based upon contractions of Gaussian primitives are developed for the third‐row elements Ga through Kr. The basis functions generalize the 6‐31G and 6‐31G* sets commonly used for atoms up to Ar. A reexamination of the 6‐31G* basis set for K and Ca developed earlier leads to the inclusion of 3d orbitals into the valence space for these atoms. Now the 6‐31G basis for the whole third‐row K through Kr has six primitive Gaussians for 1s, 2s, 2p, 3s, and 3p orbitals, and a split‐valence pair of three and one primitives for valence orbitals, which are 4s, 4p, and 3d. The nature of the polarization functions for third‐row atoms is reexamined as well. The polarization functions for K, Ca, and Ga through Kr are single set of Cartesian d‐type primitives. The polarization functions for transition metals are defined to be a single 7f set of uncontracted primitives. Comparison with experimental data shows good agreement with bond lengths and angles for representative vapor‐phase metal complexes. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 976–984, 2001     

Hartree–Fock High-Precision Analytical Functions of Open-Shell Atoms

The algorithm of high-precision optimization of basis functions suggested previously for calculating the analytical Hartree–Fock orbitals of closed-shell atoms is generalized to open-shell systems described by the Roothaan method (1960). Expressions for the first (free gradient) and second (Hesse matrix) derivatives of the system's energy with respect to the nonlinear parameters (orbital exponents) of the basis functions are derived in terms of density matrices for the filled and open shells. An algorithm is proposed for high-precision optimization of the nonlinear parameters using these equations based on Murtagh–Sargent and Newton minimization procedures. To illustrate the application of this algorithm, we give optimization of the basis sets of Slater type functions for atoms from the second row, as well as for Al, Si, P, K, Sc, and Fe atoms. The analytical Hartree–Fock orbitals giving nearly Hartree–Fock energies are calculated with a high degree of accuracy.     

An approximate frozen core approach to valence only molecular calculations

A simple valence electron-only theory based on an approximate frozen core approach and an exact core-valence strong orthogonality condition is developed for atomic and molecular systems. A unique reduced basis is introduced in which both core and valence orbitals are expanded. The core representation is roughly approximated, and the valence orbital overlap with the corresponding all-electron reference functions is nearly exact. The size of the reduced basis in terms of primitive functions is practically the same as that adopted by effective core potential methods in which the valence orbitals have the correct nodal properties. Results obtained with the present approach are presented for LiO, BeO and CaO molecules, and compared with the corresponding all-electron frozen core calculations. In addition, a detailed investigation on Li _n Be clusters (n=1,..., 6) is carried out.Dedicated to Professor J. Koutecký on the occasion of his 65th birthday     

The Bond Order of C2 from a Strictly N‐Representable Natural Orbital Energy Functional Perspective

The bond order of the ground electronic state of the carbon dimer has been analyzed in the light of natural orbital functional theory calculations carried out with an approximate, albeit strictly N‐representable, energy functional. Three distinct solutions have been found from the Euler equations of the minimization of the energy functional with respect to the natural orbitals and their occupation numbers, which expand upon increasing values of the internuclear coordinate. In the close vicinity of the minimum energy region, two of the solutions compete around a discontinuity point. The former, corresponding to the absolute minimum energy, features two valence natural orbitals of each of the following symmetries, σ, σ*, π and π*, and has three bonding interactions and one antibonding interaction, which is very suggestive of a bond order large than two but smaller than three. The latter, features one σ–σ* linked pair of natural orbitals and three degenerate pseudo‐bonding like orbitals, paired each with one triply degenerate pseudo‐antibonding orbital, which points to a bond order larger than three. When correlation effects, other than Hartree–Fock for example, between the paired natural orbitals are accounted for, this second solution vanishes yielding a smooth continuous dissociation curve. Comparison of the vibrational energies and electron ionization energies, calculated on this curve, with their corresponding experimental marks, lend further support to a bond order for C ₂ intermediate between acetylene and ethylene.     

Self‐consistent shielding in atoms and molecules

The orbital exponents of trial wave functions for simple systems can be found from the potential energy terms alone. Shielding of the nuclear charge by one electron on another is determined by the relative values of the nuclear–electron attraction and the electron–electron repulsion. For two electrons in the same orbital, the shielding is divided equally. For different orbitals, only the inner electron shields the outer. The systems tested are first‐row atoms, using Slater orbitals. It appears that if this approach can be generalized, it may not be necessary to calculate kinetic energies in chemical systems, since they will be determined by the orbital exponents. This would be useful if trial wave functions were not available, but trial electron density functions were. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Correlating basis sets for the H atom and the alkali-metal atoms from Li to Rb

 Contracted Gaussian-type function sets are developed for correlating p, d, and f functions for a valence electron of the hydrogen atom and alkali-metal atoms from Li to Rb. A segmented contraction scheme is used for its compactness and efficiency. Contraction coefficients and exponents are determined by minimizing the deviation from the K orbitals of the atoms. The present basis sets yield an accuracy comparable to the correlation-consistent basis set for the hydrogen atom and also give a similar high accuracy for the alkali-metal atoms. In the calculations of spectroscopic constants of alkali hydrides, the decontraction of the p function plays an important role, especially for LiH. The contributions of d and f functions are nontrivial for KH and RbH. Received: 6 September 2002 / Accepted: 13 November 2002 / Published online: 19 March 2003 Acknowledgements. This work was supported in part by a Grant-in-Aid for Scientific Research from the Ministry of Education of Japan. Correspondence to: T. Noro e-mail: tashi@sci.hokudai.ac.jp     

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.     

Compact and efficient basis sets of s- and p-block elements for model core potential method

We propose compact and efficient valence-function sets for s- and p-block elements from Li to Rn to appropriately describe valence correlation in model core potential (MCP) calculations. The basis sets are generated by a combination of split MCP valence orbitals and correlating contracted Gaussian-type functions in a segmented form. We provide three types of basis sets. They are referred to as MCP-dzp, MCP-tzp, and MCP-qzp, since they have the quality comparable with all-electron correlation consistent basis sets, cc-pVDZ, cc-pVTZ, and cc-pVQZ, respectively, for lighter atoms. MCP calculations with the present basis sets give atomic correlation energies in good agreement with all-electron calculations. The present MCP basis sets systematically improve physical properties in atomic and molecular systems in a series of MCP-dzp, MCP-tzp, and MCP-qzp. Ionization potentials and electron affinities of halogen atoms as well as molecular spectroscopic constants calculated by the best MCP set are in good agreement with experimental values.     

Basis set selections for relativistic self-consistent field calculations

Practical methods of generating reliable and economic basis sets for relativistic self-consistent fields (RSCF) calculations are developed. Large component basis sets are generated from constrained optimizations of exponents in the nonrelativistic atomic calculations for light atoms. For heavy atoms, large component basis sets for inner core orbitals are generated by fitting numerical atomic spinors of Dirac-Hartree-Fock calculations with appropriate number of Slater-type functions. Small component basis sets are obtained by using the kinetic balance condition and other computational criteria. With judicious selections of the basis sets, virtual orbitals in RSCF calculations become very similar to those in nonrelativistic calculations, implying that relativistic virtual orbitals can be used in electron correlation calculations in the same manner as the conventional nonrelativistic virtual orbitals. It is also evident that the Koopmans' theorem is also valid in RSCF results.