Polarization consistent basis sets. 4: the elements He, Li, Be, B, Ne, Na, Mg, Al, and Ar

Polarization consistent basis sets, optimized for density functional calculations, are proposed for the elements He, Li, Be, B, Ne, Na, Mg, Al, and Ar. The basis sets for He, B, Ne, Al, and Ar are assigned based on the previously proposed basis sets for H, C-F, and Si-Ar. The basis sets for Li, Be, Na, and Mg are defined based on energetic analysis along the lines used in previous work and the performance for molecular systems. The performance for atomization energies is comparable to those for systems composed of the elements H, C-F, and Si-Ar.     

Medium-size Gaussian basis sets for hydrogen through argon

Summary Medium-sized Gaussian basis sets are reoptimized for the ground states of the atoms from hydrogen through argon. The composition of these basis sets is (4s), (5s), and (6s) for H and He, (9s5p) and (12s7p) for the atoms Li to Ne, and (12s8p) and (12s9p) for the atoms Na to Ar. Basis sets for the² P states of Li and Na, and the³ P states of Be and Mg are also constructed since they are useful in molecular calculations. In all cases, our energies are lower than those obtained previously with Gaussian basis sets of the same size.     

A systematic preparation of new contracted Gaussian-type orbitals. IX [54/5], [64/5], [64/6], [74/6], [74/7] and MAXI-1–MAXI-5 from Li to Ne

The new contracted Gaussian-type orbitals (CGTO s) for molecular calculations have been developed from Li to Ne. The CGTO s are minimal type, i.e. composed of two s-type CGTO s, s₁, s₂, and one p-type CGTO , p₁. They are new family of CGTO s given by Tatewaki and Huzinaga, and others. In the previous works three primitive GTO s are used for s₂, which is the main part of the 2s orbital, whereas four primitive GTO s are employed in the present work. The sets generated are [54/5], [64/5], [64/6], [74/6], and [74/7]. In almost all the cases the errors in the 2s and 2p orbital energies are smaller than those of DZ . The resulting 2s orbitals are close to the orbitals of the uncontracted GTO sets, (13/n) and (14/n) of Duijneveldt. It is found that the 2s and 2p orbitals given by [64/6], [74/6], and [74/7] are satisfactorily near to those of Hartree–Fock. The basis sets [54/5], [64/6], and [74/7] are applied to the N₂ molecule in the split valence forms of [5211/311], [6211/3111], and [7211/4111]. Adding the d-type polarization functions from one through three, the quality of the basis sets has been examined. All of the three sets show good behavior and the sets augmented with three d-type polarization functions give almost entirely the same results as the very extended basis set.     

Atomic SCF loge localized wave functions

We report in this work SCF atomic calculations for Li, Be, B, C, and Ne using a basis set of completely loge-localized functions. For these second row atoms the total volume R³ was partitioned into a spherical loge of radius R and its volume complement. The loge-localized basis functions were constructed as a product of Slater-type orbitals and a cut-off factor. The energy values obtained differ significantly from the Hartree–Fock ones indicating that the delocalization effects—not included in these calculations—are important.     

Calculation of Negative Ions of B,C, N,O and F Using Noninteger n Slater Type Orbitals

Combined Hartree‐Fock‐Roothaan calculations have been performed using noninteger n Slater type orbitals for the ground states of the lowest electron configurations 1s²2s²2pⁿ (2 ≤ n ≤ 6) for negative ions of B, C, N, O and F. These results are compared with the corresponding results obtained from the use of integer n Slater type orbitals. All of the nonlinear parameters are fully optimized. The results of calculation of coupling‐projection coefficients, orbital and total energies and virial ratios are presented. It is shown that the noninteger n Slater type orbitals, in general, improve the orbital energies.     

Model study of the impact of orbital choice on the accuracy of coupled-cluster energies. II. Valence-universal coupled-cluster method

The impact of the choice of molecular orbital sets on the results of the valence-universal coupled cluster method involving up to three-body amplitudes (VU-CCSDT) was studied for the H4 model. This model offers a straightforward way of representing all possible symmetry-adapted orbitals. Moreover, the degree of quasi-degeneracy of its lowest ¹A₁ states can be varied over a wide range by changing its geometry. Calculations were performed both for 13 sets of standard quantum chemical orbitals and for a vast variety of nonstandard orbital sets defined by nodes of a two-dimensional orbital grid. The performance of various standard orbital sets in VU-CCSDT calculations is compared. It is also documented that for every quasi-degeneracy region there exist nonstandard orbital sets which allow one to obtain more accurate VU-CCSDT energies than the standard orbital sets. In an attempt to provide a general interpretation for some of the alternative orbital sets, we defined a set of orbitals which maximize the proximity of the model and target spaces—maximum proximity orbitals (MPO). It is demonstrated that outside the strong quasi-degeneracy region the energies obtained for the VU-CCSDT approach based on the MPOs are more accurate than for the standard orbital sets. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 67: 221–237, 1998     

Valence p functions for alkali and alkaline-earth atoms

Contracted Gaussian-type function (CGTF) basis sets are reported for valence p orbitals of the six alkali and alkaline-earth atoms Li, Be, Na, Mg, K, and Ca for molecular applications. These sets are constructed by Roothaan–Hartree–Fock calculations for the ns → np excited states of atoms, in which both linear and nonlinear parameters of CGTFs are variationally optimized. The present CGTF sets reproduce well the numerical Hartree–Fock ns → np excitation energies: the largest error is 0.0009 hartrees for Li. New CGTFs are tested with diatomic Li₂, Na₂, K₂, and MH molecules, where M = Li, Be, Na, Mg, K, and Ca, by self-consistent-field (SCF) and multiconfiguration SCF calculations. The resultant spectroscopic constants compare well with those of more elaborate calculations and are sufficiently close to experimental values, supporting the efficiency of the present set for the valence p orbitals. Received: 9 July 1998 / Accepted: 17 September 1998 / Published online: 1 February 1999     

What do kinetic-energy anisotropies tell us about chemical bonding? I. Diatomic hydrides

The kinetic-energy anisotropies of fifteen diatomic hydrides AH with A = H, Li, Be, B, C, N, O, F, Na, Mg, Al, Si, P, S, Cl are calculated from self-consistent-field wave functions constructed from extended basis sets of Slater-type orbitals. It is found that there is no consistent ordering of the bond-parallel and bond-perpendicular components of the kinetic energy with respect to separated atom values. An analysis of the orbital contributions reveals that nonbonding π orbitals make large contributions to the total kinetic-energy anisotropy. This makes it difficult, if not impossible, to deduce anything about the nature of the chemical bond from the total anisotropy. However, certain dimensionless orbital kinetic-energy anisotropies are useful for interpretative studies because, in free atoms, these quantities have fixed values that depend only on the symmetry of the orbital.     

An ab initio comparison of dichalcogen hydrides

A theoretical study of HSH, HSeH, HTeH, HSSH, HSeSH, HSeSeH, HTeSH, HTeSeH, and HTeTeH was carried out by ab initio molecular orbital methods employing minimal Gaussian basis sets MINI-1 and MINI-1* of Huzinaga and his group. Both basis sets yield accurate estimates on the equilibrium geometries of monochalcogen hydrides. In the case of dichalcogen hydrides, however, the inclusion of d-polarization functions for sulfur, selenium, and tellurium greatly improve the accuracy of the geometry prediction. The unpolarized MINI-1 basis sets yield essentially correct orbital energies and therefore suffice for the comparative study on the electronic structures in similar molecules. The results with both basis sets imply close similarities in the electronic structures of SS, SeS, and SeSe bonds with more marked differences in bonds containing tellurium as a consequence of notably smaller orbital energy of the 5s- orbital of tellurium as compared to the corresponding orbitals in sulfur and selenium. The barriers to internal rotation about the chalcogen—chalcogen bond in all dichalcogen hydrides are similar. The cis- and trans-barrier heights are ca. 23 and 14 kJ mol⁻¹, respectively. The relative stabilities of different hydrides are discussed.     

Polarized basis sets of Slater-type orbitals: H to Ne atoms

We present three Slater-type atomic orbital (STO) valence basis (VB) sets for the first and second row atoms, referred to as the VB1, VB2, and VB3 bases. The smallest VB1 basis has the following structure: [3, 1] for the H and He atoms, [5, 1] for Li and Be, and [5, 3, 1] for the B to Ne series. For the VB2 and VB3 bases, both the number of shells and the number of functions per shell are successively increased by one with respect to VB1. With the exception of the H and Li atoms, the exponents for the VB1 bases were obtained by minimizing the sum of the Hartree-Fock (HF) and frozen-core singles and doubles configuration interaction (CISD FC) energies of the respective atoms in their ground state. For H and Li, we minimized the sum of the HF and CISD FC energies of the corresponding diatoms (i.e., of H(2) or Li(2)) plus the ground-state energy of the atom. In the case of the VB2 basis sets, the sum that was minimized also included the energies of the positive and negative ions, and for the VB3 bases, the energies of a few lowest lying excited states of the atom. To account for the core correlations, the VBx (x = 1, 2, and 3) basis sets for the Li to Ne series were enlarged by one function per shell. The exponents of these extended (core-valence, CV) basis sets, referred to, respectively, as the CVBx (x = 1, 2, and 3) bases, were optimized by relying on the same criteria as in the case of the VBx (x = 1, 2, and 3) bases, except that the full CISD rather than CISD FC energies were employed. We show that these polarized STO basis sets provide good HF and CI energies for the ground and excited states of the atoms considered, as well as for the corresponding ions.     

Energies of the first row atoms from quantum Monte Carlo

All-electron variational and diffusion quantum Monte Carlo calculations of the ground state energies of the first row atoms (from Li to Ne) are reported. The authors use trial wave functions of four types: single-determinant Slater-Jastrow wave functions, multideterminant Slater-Jastrow wave functions, single-determinant Slater-Jastrow wave functions with backflow transformations, and multideterminant Slater-Jastrow wave functions with backflow transformations. At the diffusion quantum Monte Carlo level and using their multideterminant Slater-Jastrow wave functions with backflow transformations, they recover 99% or more of the correlation energies for Li, Be, B, C, N, and Ne, 97% for O, and 98% for F.     

Compact wave functions for the beryllium isoelectronic series, Li- to Ne6+: a standard Hylleraas approach

Variational calculations have been carried out for the ground states of several members of the beryllium isoelectronic series using a standard Hylleraas approach involving Slater-type basis functions. The species examined are Li(-), Be, B(+), C(2+), N(3+), O(4+), F(5+), and Ne(6+). For each species, the nonrelativistic energy, the electronic density at the nucleus, the expectation value , the moments for n = -1, 1, 2, and 3, and for n = -1, 1, and 2, are reported. With relatively compact basis sets, the ground state energies are obtained with uncertainties ranging from 50 parts per million to just under 4 parts per million.     

Optimized effective potential method: is it possible to obtain an accurate representation of the response function for finite orbital basis sets?

The optimized effective potential (OEP) equations are solved in a matrix representation using the orbital products of occupied and virtual orbitals for the representation of both the local potential and the response function. This results in a direct relationship between the matrix elements of local and nonlocal operators for the exchange-correlation potential. The effect of the truncation of the number of such products in the case of finite orbital basis sets on the OEP orbital and total energies and on the spectrum of eigenvalues of the response function is examined. Test calculations for Ar and Ne show that rather large AO basis sets are needed to obtain an accurate representation of the response function.     

Orbital relaxation in the Rydberg series of the lithium atom. An excited state SCF calculation

Self-consistent-field (SCF ) calculations for a series of Rydberg states (1s²ns)²S of the Li atom are performed using the generalized Brillouin theorem (GBT) method. The calculated energy is a proper upper bound to the excited state energy. The SCF term values of the Rydberg states are almost the same as those of the frozen-core approximation ones. The orbital behavior shows that the core is slightly expanded by the penetration of the Rydberg orbitals, and the higher Rydberg orbitals can be very well represented by the modified hydrogen-like orbitals.     

Ground states of molecules. 53.MNDO calculations for molecules containing chlorine

MNDO calculations of heats of formation, dipole moments, ionization potentials, and structures are reported for a wide range of compounds containing chlorine in its characteristic valence state (Cl^I) and one or more of the elements H, B, Be, C, N, O, and F. The calculated errors in the heats of formation and the dipole moments are not significantly greater than those previously reported for compounds containing no chlorine. First vertical ionization potentials were on average 0.95 eV too high. The ordering of higher cationic states was found to be correct, even for species such as Cl₂O, Cl₂, and HOCl, where ab initio–Koopmans' theorem calculations predict the incorrect ordering. The calculated energies and geometries of compounds such as CIF₃ are qualitatively incorrect, probably because of the lack of 3d atomic orbitals in the orbital basis set.     

Compact contracted Gaussian-type basis sets for halogen atoms. Basis-set superposition effects on molecular properties

Compact, contracted Gaussian basis sets for halogen atoms are generated and tested in ab initio molecular calculations. These basis sets have similar structure to that of Huzinaga and co-workers' (HTS ) sets; however, they give both better atomic total energies and better properties of atomic valence orbitals. These sets, after splitting of valence orbitals and augmenting with polarization functions, provide molecular results that agree well with those given by extended calculations. Basis set superposition error (BSSE ) is calculated using the counterpoise method. BSSE has only slight influence on calculated equilibrium geometry, shape of potential curve, and electric properties (dipole and quadrupole moments) of molecules. However, atomization energies may be significantly changed by the BSSE .     

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.     

Intramolecular hydrogen bond in the hydroxycyclohexadienyl peroxy radicals

The hydroxycyclohexadienyl peroxy radicals (HO? C₆H₆? O₂) produced from the reaction of OH‐benzene adduct with O₂ were studied with density functional theory (DFT) calculations to determine their characteristics. The optimized geometries, vibrational frequencies, and total energies of 2‐hydroxycyclohexadienyl peroxy radical II_s and 4‐hydroxycyclohexadienyl peroxy radical III_s were calculated at the following theoretical levels, B3LYP/6‐31G(d), B3LYP/6‐311G(d,p), and B3LYP/6‐311+G(d,p). Both were shown to contain a red‐shifted intramolecular hydrogen bond (O? H … O? H bond). According to atoms‐in‐molecules (AIM) analysis, the intramolecular hydrogen bond in the 2‐hydroxycyclohexadienyl peroxy radical II_s is stronger than that one in 4‐hydroxycyclohexadienyl peroxy radical III_s, and the former is the most stable conformation among its isomers. Generally speaking, hydrogen bonding in these radicals plays an important role to make them more stable. Based on natural bond orbital (NBO) analysis, the stabilization energy between orbitals is the main factor to produce red‐shifted intramolecular hydrogen bond within these peroxy radicals. The hyperconjugative interactions can promote the transfer of some electron density to the O? H antibonding orbital, while the increased electron density in the O? H antibonding orbital leads to the elongation of the O? H bond and the red shift of the O? H stretching frequency. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Orbitals expanded in slater functions with single-exponent by shell and by subshell

Slater-type orbitals (STO s) with a single-exponent by shell or by subshell have been constructed to reduce the number of integrals evaluated in the electronic calculations. The expansion of orbitals in these new basis sets has been carried out in detail for the ground state of the Ne atom. We have carried out a study of STO basis sets with a different size for this atom that could help to propose empirical rules for the selection of these basis sets for other atoms. The usefulness of STO s with single-exponent by shell and subshell and the splitting of s and p functions are discussed. © 1994 John Wiley & Sons, Inc.