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.     

Augmented Gaussian basis set of quintuple zeta valence quality for H and from Li to Ar: Applications in DFT calculations of molecular electric properties

Augmented Gaussian basis set of quintuple zeta valence quality plus polarization functions (A5ZP) for H and Li–Ar is presented. It was determined from the 5ZP basis set by addition of diffuse (s and p symmetries) and polarization (p, d, f, g, and h symmetries) functions that were optimized for the anion at the Hartree–Fock and Mller–Plesset second-order levels, respectively. It was shown that in general this basis set in combination with the density functional theory can be used with success to predict electric properties for a sample of molecules. Comparison with theoretical and experimental values available in the literature is done.     

Basis set and electron correlation effects on the polarizability and second hyperpolarizability of model open-shell pi-conjugated systems

The basis set and electron correlation effects on the static polarizability (alpha) and second hyperpolarizability (gamma) are investigated ab initio for two model open-shell pi-conjugated systems, the C(5)H(7) radical and the C(6)H(8) radical cation in their doublet state. Basis set investigations evidence that the linear and nonlinear responses of the radical cation necessitate the use of a less extended basis set than its neutral analog. Indeed, double-zeta-type basis sets supplemented by a set of d polarization functions but no diffuse functions already provide accurate (hyper)polarizabilities for C(6)H(8) whereas diffuse functions are compulsory for C(5)H(7), in particular, p diffuse functions. In addition to the 6-31G(*)+pd basis set, basis sets resulting from removing not necessary diffuse functions from the augmented correlation consistent polarized valence double zeta basis set have been shown to provide (hyper)polarizability values of similar quality as more extended basis sets such as augmented correlation consistent polarized valence triple zeta and doubly augmented correlation consistent polarized valence double zeta. Using the selected atomic basis sets, the (hyper)polarizabilities of these two model compounds are calculated at different levels of approximation in order to assess the impact of including electron correlation. As a function of the method of calculation antiparallel and parallel variations have been demonstrated for alpha and gamma of the two model compounds, respectively. For the polarizability, the unrestricted Hartree-Fock and unrestricted second-order M?ller-Plesset methods bracket the reference value obtained at the unrestricted coupled cluster singles and doubles with a perturbative inclusion of the triples level whereas the projected unrestricted second-order M?ller-Plesset results are in much closer agreement with the unrestricted coupled cluster singles and doubles with a perturbative inclusion of the triples values than the projected unrestricted Hartree-Fock results. Moreover, the differences between the restricted open-shell Hartree-Fock and restricted open-shell second-order M?ller-Plesset methods are small. In what concerns the second hyperpolarizability, the unrestricted Hartree-Fock and unrestricted second-order M?ller-Plesset values remain of similar quality while using spin-projected schemes fails for the charged system but performs nicely for the neutral one. The restricted open-shell schemes, and especially the restricted open-shell second-order M?ller-Plesset method, provide for both compounds gamma values close to the results obtained at the unrestricted coupled cluster level including singles and doubles with a perturbative inclusion of the triples. Thus, to obtain well-converged alpha and gamma values at low-order electron correlation levels, the removal of spin contamination is a necessary but not a sufficient condition. Density-functional theory calculations of alpha and gamma have also been carried out using several exchange-correlation functionals. Those employing hybrid exchange-correlation functionals have been shown to reproduce fairly well the reference coupled cluster polarizability and second hyperpolarizability values. In addition, inclusion of Hartree-Fock exchange is of major importance for determining accurate polarizability whereas for the second hyperpolarizability the gradient corrections are large.     

Valence basis sets for lanthanide 4f-in-core pseudopotentials adapted for crystal orbital ab initio calculations

Crystal orbital adapted Gaussian (4s4p3d), (5s5p4d) and (6s6p5d) valence primitive basis sets have been derived for calculating periodic bulk materials containing trivalent lanthanide ions modeled with relativistic energy-consistent 4f-in-core lanthanide pseudopotentials of the Stuttgart-Koeln variety. The calibration calculations of crystalline A-type Pm₂O₃ using different segmented contraction schemes (4s4p3d)/[2s2p2d], (4s4p3d)/[3s3p2d], (5s5p4d)/[2s2p2d], (5s5p4d)/[3s3p3d], (5s5p4d)/[4s4p3d], (6s6p5d)/[2s2p2d], (6s6p5d)/[3s3p3d] and (6s6p5d)/[4s4p4d] are discussed at both Hartree–Fock (HF) and density functional theory (DFT) levels for the investigation of basis set size effects. Applications to the geometry optimization of A-type Ln₂O₃ (Ln = La-Pm) show a satisfactory agreement with experimental data using the lanthanide valence basis sets (6s6p5d)/[4s4p4d] and the standard set 6-311G* for oxygen. The corresponding augmented sets (8s7p6d)/[6s5p5d] with additional diffuse functions for describing neutral lanthanide atoms were applied to calculate atomic energies of free lanthanide atoms for the evaluation of cohesive energies for A-Ln₂O₃ within both conventional Kohn-Sham DFT and the a posteriori-HF correlation DFT schemes.     

Density functional theory optimized basis sets for gradient corrected functionals: 3d transition metal systems

Density functional theory optimized basis sets for gradient corrected functionals for 3d transition metal atoms are presented. Double zeta valence polarization and triple zeta valence polarization basis sets are optimized with the PW86 functional. The performance of the newly optimized basis sets is tested in atomic and molecular calculations. Excitation energies of 3d transition metal atoms, as well as electronic configurations, structural parameters, dissociation energies, and harmonic vibrational frequencies of a large number of molecules containing 3d transition metal elements, are presented. The obtained results are compared with available experimental data as well as with other theoretical data from the literature.     

Quasirelativistic energy-consistent 5f-in-core pseudopotentials for trivalent actinide elements

Quasirelativistic energy-consistent 5f-in-core pseudopotentials modelling trivalent actinides, corresponding to a near-integral 5f ⁿ occupation (n = 0–14 for Ac–Lr), have been generated. Energy-optimized (6s5p4d), (7s6p5d), and (8s7p6d) primitive valence basis sets contracted to polarized double to quadruple zeta quality as well as 2f1g correlation functions have been derived. Corresponding smaller basis sets (4s4p3d), (5s5p4d), and (6s6p5d) suitable for calculations on actinide(III) ions in crystalline solids form subsets of these basis sets designed for calculations on neutral molecules. Results of Hartree–Fock test calculations for actinide(III) monohydrates and actinide trifluorides show a satisfactory agreement with corresponding calculations using 5f-in-valence pseudopotentials. Even in the beginning of the actinide series, where the 5f shell is relatively diffuse, only quite acceptable small deviations occur as long as the 5f-shell does not participate significantly in covalent bonding. Electronic Supplementary Material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Molecular structure of diatomic lanthanide compounds

The molecular constants of selected diatomic lanthanide compounds (LaH, LaO, LaF, EuH, EuO, EuF, EuS, GdO, GdF, GdH, YbH, YbO, YbF, YbS, LuH, LuO and LuF) have been calculated by using relativistic small-core pseudopotentials and optimized (14s13p10d8f6g)/ [6s6p5d4f3g] valence basis sets. The results are in good agreement with available experimental data, with exception of YbO and LuF. The reasons for the discrepancies in case of YbO are due to a complicated mixing of configurations in the ground state, whereas in case of LuF the binding energy estimated by experimentalists appears to be too low.     

Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy

Gaussian basis sets of quadruple zeta valence quality for Rb-Rn are presented, as well as bases of split valence and triple zeta valence quality for H-Rn. The latter were obtained by (partly) modifying bases developed previously. A large set of more than 300 molecules representing (nearly) all elements-except lanthanides-in their common oxidation states was used to assess the quality of the bases all across the periodic table. Quantities investigated were atomization energies, dipole moments and structure parameters for Hartree-Fock, density functional theory and correlated methods, for which we had chosen M?ller-Plesset perturbation theory as an example. Finally recommendations are given which type of basis set is used best for a certain level of theory and a desired quality of results.     

Optimized Slater-type basis sets for the elements 1-118

Seven different types of Slater type basis sets for the elements H (Z = 1) up to E118 (Z = 118), ranging from a double zeta valence quality up to a quadruple zeta valence quality, are tested in their performance in neutral atomic and diatomic oxide calculations. The exponents of the Slater type functions are optimized for the use in (scalar relativistic) zeroth-order regular approximated (ZORA) equations. Atomic tests reveal that, on average, the absolute basis set error of 0.03 kcal/mol in the density functional calculation of the valence spinor energies of the neutral atoms with the largest all electron basis set of quadruple zeta quality is lower than the average absolute difference of 0.16 kcal/mol in these valence spinor energies if one compares the results of ZORA equation with those of the fully relativistic Dirac equation. This average absolute basis set error increases to about 1 kcal/mol for the all electron basis sets of triple zeta valence quality, and to approximately 4 kcal/mol for the all electron basis sets of double zeta quality. The molecular tests reveal that, on average, the calculated atomization energies of 118 neutral diatomic oxides MO, where the nuclear charge Z of M ranges from Z = 1-118, with the all electron basis sets of triple zeta quality with two polarization functions added are within 1-2 kcal/mol of the benchmark results with the much larger all electron basis sets, which are of quadruple zeta valence quality with four polarization functions added. The accuracy is reduced to about 4-5 kcal/mol if only one polarization function is used in the triple zeta basis sets, and further reduced to approximately 20 kcal/mol if the all electron basis sets of double zeta quality are used. The inclusion of g-type STOs to the large benchmark basis sets had an effect of less than 1 kcal/mol in the calculation of the atomization energies of the group 2 and group 14 diatomic oxides. The basis sets that are optimized for calculations using the frozen core approximation (frozen core basis sets) have a restricted basis set in the core region compared to the all electron basis sets. On average, the use of these frozen core basis sets give atomic basis set errors that are approximately twice as large as the corresponding all electron basis set errors and molecular atomization energies that are close to the corresponding all electron results. Only if spin-orbit coupling is included in the frozen core calculations larger errors are found, especially for the heavier elements, due to the additional approximation that is made that the basis functions are orthogonalized on scalar relativistic core orbitals.     

Quasirelativistic energy-consistent 5f-in-core pseudopotentials for divalent and tetravalent actinide elements

Quasirelativistic energy-consistent 5f-in-core pseudopotentials modeling divalent (5f ⁿ⁺¹ occupation with n = 5–13 for Pu–No) respectively tetravalent (5f ^n-1 occupation with n = 1–9 for Th–Cf) actinides together with corresponding core-polarization potentials have been generated. Energy-optimized (6s5p4d) and (7s6p5d) valence basis sets as well as 2f1g correlation functions have been derived and contracted to polarized double, triple, and quadruple zeta quality. Corresponding smaller (4s4p) and (5s5p) respectively (4s4p3d) and (5s5p4d) basis sets suitable for calculations on actinide(II) respectively actinide(IV) ions in crystalline solids form subsets of these basis sets designed for calculations on molecules. Results of Hartree–Fock test calculations for actinide di- and tetrafluorides show a satisfactory agreement with calculations using 5f-in-valence pseudopotentials. Electronic Supplementary Material The online version of this article doi: contains supplementary material, which is available to authorized users.     

Density functional study of platinum polyyne monomer,oligomer, and polymer: Ground state geometrical and electronic structures

Geometries of monomers and oligomers of a platinum polyyne and its free ligands were optimized using density functional theory with B3LYP hybrid functional. The LANL2DZ basis set was used for Pt and the 6‐31G* for other atoms in geometry optimizations. The electronic structures of these compounds were analyzed using Stuttgart/Dresden ECPs (SDD) basis set for metal atoms and 6‐311G* for others. The polymerization has very little effect on the bond lengths and by introducing the metal, the acetylide bond length increases slightly. The strong overlap between metal sp_x orbitals and σp_x orbitals of acetylides results in localized σ bonding. The hybridization between the ligand pπ orbitals and the platinum dπ orbital resulted in the π‐conjugation enhancement. This conjugation enhancement causes some effects such as the highest‐occupied molecular orbital–lowest‐unoccupied molecular orbital gap reduction and charge transfer characteristic of low‐energy vertical transitions. © 2013 Wiley Periodicals, Inc.     

DFT Calculations of the Ionization Potential and Electron Affinity of Alaninamide

Adiabatic and vertical ionization potentials (IPs) and valence electron affinities (EAs) of alaninamide in gas phase have been determined using density functional theory (BLYP,B3LYP,B3P86) methods with 6-311++G(d,p) basis set,respectively. IPs and EAs of alaninamide in solutions have been calculated at the B3LYP/6-311++G(d,p) level. Five possible conformers of alaninamide and their charged states have been optimized employing density functional theory B3LYP method with 6-311++(d,p) basis set,respectively.     

Augmented Gaussian basis sets of double and triple zeta valence qualities for the atoms K and Sc–Kr: Applications in HF, MP2, and DFT calculations of molecular electric properties

Augmented Gaussian basis sets of double and triple zeta valence qualities plus polarization functions for the atoms K and from Sc to Kr are presented. They were generated from the all-electron unaugmented sets by addition of diffuse functions (s, p, d, f, and g symmetries) that were optimized for the anion ground states. From these sets, Hartree–Fock, second-order Møller–Plesset perturbation theory, and density functional theory electric dipole moment and dipole polarizability calculations for a sample of molecules were carried out. Comparison with theoretical and experimental values available in the literature was done.