Basis Set Of Slater-type Orbitals Research Articles

Optimization of the r2SCAN-3c Composite Electronic-Structure Method for Use with Slater-Type Orbital Basis Sets.

The “Swiss army knife” composite density functional electronic-structure method r2SCAN-3c (J. Chem. Phys.2021, 154, 064103) is extended and optimized for the use with Slater-type orbital basis sets. The meta generalized-gradient approximation (meta-GGA) functional r2SCAN by Furness et al. is combined with a tailor-made polarized triple-ζ Slater-type atomic orbital (STO) basis set (mTZ2P), the semiclassical London dispersion correction (D4), and a geometrical counterpoise (gCP) correction. Relativistic effects are treated explicitly with the scalar-relativistic zeroth-order regular approximation (SR-ZORA). The performance of the new implementation is assessed on eight geometry and 74 energy benchmark sets, including the extensive GMTKN55 database as well as recent sets such as ROST61 and IONPI19. In geometry optimizations, the STO-based r2SCAN-3c is either on par with or more accurate than the hybrid density functional approximation M06-2X-D3(0)/TZP. In energy calculations, the overall accuracy is similar to the original implementation of r2SCAN-3c with Gaussian-type atomic orbitals (GTO), but basic properties, intermolecular noncovalent interactions, and barrier heights are better described with the STO approach, resulting in a lower weighted mean absolute deviation (WTMAD-2(STO) = 7.15 vs 7.50 kcal mol–1 with the original method) for the GMTKN55 database. The STO-optimized r2SCAN-3c outperforms many conventional hybrid/QZ approaches in most common applications at a fraction of their cost. The reliable, robust, and accurate r2SCAN-3c implementation with STOs is a promising alternative to the original implementation with GTOs and can be generally used for a broad field of quantum chemical problems.

Potential-energy curve for the a3Σu+ state of a lithium dimer with Slater-type orbitals

We report state-of-the-art ab initio calculations of the potential energy curve for the $a^3\Sigma_u^+$ state of the lithium dimer conducted to achieve spectroscopic accuracy ($<$1cm$^{-1}$) without any prior adjustment to fit the corresponding experimental data. The nonrelativistic clamped-nuclei component of the interaction energy is calculated with a composite method involving six-electron coupled cluster and full configuration interaction theories combined with basis sets of Slater-type orbitals ranging in quality from double- to sextuple-zeta. To go beyond the nonrelativistic Born-Oppenheimer picture we include both the leading-order relativistic and adiabatic corrections, and find both of these effects to be non-negligible within the present accuracy standards. The potential energy curve developed by us allowed to calculate molecular parameters ($D_e$, $D_0$, $\omega_e$ etc.) for this system, as well as the corresponding vibrational energy levels, with an error of only a few tenths of a wavenumber ($0.2-0.4\,$cm$^{-1}$). We also report an ab initio value for the scattering length of two $^2S$ lithium atoms which determines the stability of the related Bose-Einstein condensate.

An efficient hybrid scheme for time dependent density functional theory.

A hybrid approach able to perform Time Dependent Density Functional Theory (TDDFT) simulations with the same accuracy as that of hybrid exchange-correlation (xc-) functionals but at a fraction of the computational cost is developed, implemented, and validated. The scheme, which we name Hybrid Diagonal Approximation (HDA), consists in employing in the response function a hybrid xc-functional (containing a fraction of the non-local Hartree-Fock exchange) only for the diagonal elements of the omega matrix, while the adiabatic local density approximation is employed for the off-diagonal terms. HDA is especially (but not exclusively) advantageous when using Slater type orbital basis sets and allows one to employ them in a uniquely efficient way, as we demonstrate here by implementing HDA in a local version of the Amsterdam Density Functional code. The new protocol is tested on NH3, C6H6, and the [Au25(SCH3)18]- cluster as prototypical cases ranging from small molecules to ligand-protected metal clusters, finding excellent agreement with respect to both full kernel TDDFT simulations and experimental data. Additionally, a specific comparison test between full kernel and HDA is considered at the Casida level on seven other molecular species, which further confirm the accuracy of the approach for all investigated systems. For the [Au25(SCH3)18]- cluster, a speedup by a factor of seven is obtained with respect to the full kernel. The HDA, therefore, promises to provide a quantitative description of the optical properties of medium-sized systems (nanoclusters) at an affordable cost, thanks to its computational efficiency, especially in combination with a complex polarization algorithm version of TDDFT.

A Quadratic Pair Atomic Resolution of the Identity Based SOS-AO-MP2 Algorithm Using Slater Type Orbitals

We report a production level implementation of pair atomic resolution of the identity (PARI) based second-order Møller–Plesset perturbation theory (MP2) in the Slater type orbital (STO) based Amsterdam Density Functional (ADF) code. As demonstrated by systematic benchmarks, dimerization and isomerization energies obtained with our code using STO basis sets of triple-ζ-quality show mean absolute deviations from Gaussian type orbital, canonical, basis set limit extrapolated, global density fitting (DF)-MP2 results of less than 1 kcal/mol. Furthermore, we introduce a quadratic scaling atomic orbital based spin-opposite-scaled (SOS)-MP2 approach with a very small prefactor. Due to a worst-case scaling of , our implementation is very fast already for small systems and shows an exceptionally early crossover to canonical SOS-PARI-MP2. We report computational wall time results for linear as well as for realistic three-dimensional molecules and show that triple-ζ quality calculations on molecules of several hundreds of atoms are only a matter of a few hours on a single compute node, the bottleneck of the computations being the SCF rather than the post-SCF energy correction.

Ab initio Potential Energy Curve for the Ground State of Beryllium Dimer.

This work concerns ab initio calculations of the complete potential energy curve and spectroscopic constants for the ground state X1Σ g+ of the beryllium dimer, Be2. High accuracy and reliability of the results is one of the primary goals of the paper. To this end, we apply large basis sets of Slater-type orbitals combined with high-level electronic structure methods including triple and quadruple excitations. The effects of the relativity are also fully accounted for in the theoretical description. For the first time the leading-order quantum electrodynamics effects are fully incorporated for a many-electron molecule. Influence of the finite nuclear mass corrections (post-Born-Oppenheimer effects) turns out to be completely negligible for this system. The predicted well-depth ( De = 934.6 ± 2.5 cm-1) and the dissociation energy ( D0 = 807.7 cm-1) are in a very good agreement with the most recent experimental data. We confirm the existence of the weakly bound twelfth vibrational level [Patkowski et al. Science 2009, 326, 1382] that it lies just below the onset of the continuum.

Jj-Coupling-based atomic self-consistent-field calculations with relativistic effective core potentials and two-component spinors

A self-consistent-field (SCF) program for the calculation of atomic energies and wave functions defined in jj-coupling using two-component atomic spinors and relativistic effective core potentials (RECPs) is described. The code is based on the linear combination of atomic orbitals SCF algorithm for atomic states defined in LS-coupling developed by Roothaan and Bagus. Hamiltonian matrix elements with respect to one- and two-electron operators, including RECPs, are calculated for two-component atomic spinor basis functions of either Gaussian-type orbitals (GTOs) or Slater-type orbitals (STOs). Electronic states are defined as eigenfunctions of the total angular momentum squared operator and tables of the required vector coupling coefficients that define such pure states are provided. In addition, one or more GTO expansions of large- and small-core RECPs and their corresponding GTO basis sets are supplied for all elements Z=3 through Z=118. Optimized two-component basis sets of STOs or GTOs can be calculated for use in molecular structure codes based on RECPs and relativistic electronic structure theory. Atomic asymptotic state energies at the SCF level of theory for analysis of molecular dissociation limits can be studied. Program summaryProgram title: jjatomProgram Files doi:http://dx.doi.org/10.17632/zs4twp8r67.1Licensing provisions: GPLv3Programming language: Fortran 90Nature of problem: Relativistic atomic energies and wave functions; optimization of two component spinor basis sets.Solution method: Atomic spinors defined in jj-coupling and expansions in two-component spinor basis functions of Gaussian- or Slater-type functions. Self-consistent field optimization of the total energy. Used for optimization of two-component atomic spinor basis sets and calculations of valence spectra of heavy elements.

Combining Slater-type orbitals and effective core potentials

We present a general methodology to evaluate matrix elements of the effective core potentials (ECPs) within one-electron basis set of Slater-type orbitals (STOs). The scheme is based on translation of individual STO distributions in the framework of Barnett-Coulson method. We discuss different types of integrals which naturally appear and reduce them to few basic quantities which can be calculated recursively or purely numerically. Additionally, we consider evaluation of the STOs matrix elements involving the core polarisation potentials (CPP) and effective spin-orbit potentials. Construction of the STOs basis sets designed specifically for use with ECPs is discussed and differences in comparison with all-electron basis sets are briefly summarised. We verify the validity of the present approach by calculating excitation energies, static dipole polarisabilities and valence orbital energies for the alkaline earth metals (Ca, Sr, Ba). Finally, we evaluate interaction energies, permanent dipole moments and ionisation energies for barium and strontium hydrides, and compare them with the best available experimental and theoretical data.

Calibrating Reaction Enthalpies: Use of Density Functional Theory and the Correlation Consistent Composite Approach in the Design of Photochromic Materials.

Acquisition of highly accurate energetic data for chromium-containing molecules and various chromium carbonyl complexes is a major step toward calibrating bond energies and thermal isomerization energies from mechanisms for Cr-centered photochromic materials being developed in our laboratories. The performance of six density functionals in conjunction with seven basis sets, utilizing Gaussian-type orbitals, has been evaluated for the calculation of gas-phase enthalpies of formation and enthalpies of reaction at 298.15 K on various chromium-containing systems. Nineteen molecules were examined: Cr(CO)6, Cr(CO)5, Cr(CO)5(C2H4), Cr(CO)5(C2ClH3), Cr(CO)5(cis-(C2Cl2H2)), Cr(CO)5(gem-(C2Cl2H2)), Cr(CO)5(trans-(C2Cl2H2)), Cr(CO)5(C2Cl3H), Cr(CO)5(C2Cl4), CrO2, CrF2, CrCl2, CrCl4, CrBr2, CrBr4, CrOCl2, CrO2Cl2, CrOF2, and CrO2F2. The performance of 69 density functionals in conjunction with a single basis set utilizing Slater-type orbitals (STO) and a zeroth-order relativistic approximation was also evaluated for the same test set. Values derived from density functional theory were compared to experimental values where available, or values derived from the correlation consistent composite approach (ccCA). When all reactions were considered, the functionals that exhibited the smallest mean absolute deviations (MADs, in kcal mol-1) from ccCA-derived values were B97-1 (6.9), VS98 (9.0), and KCIS (9.4) in conjunction with quadruple-ζ STO basis sets and B97-1 (9.3) in conjunction with cc-pVTZ basis sets. When considering only the set of gas-phase reaction enthalpies (ΔrH°gas), the functional that exhibited the smallest MADs from ccCA-derived values were B97-1 in conjunction with cc-pVTZ basis sets (9.1) and PBEPBE in conjunction with polarized valence triple-ζ basis set/effective core potential combination for Cr and augmented and multiple polarized triple-ζ Pople style basis sets (9.5). Also of interest, certainly because of known cancellation of errors, PBEPBE with the least-computationally expensive basis set combination considered in the present study (valence double-ζ basis set/effective core potential combination for Cr and singly-polarized double-ζ Pople style basis sets) also provided reasonable accuracy (11.1). An increase in basis set size was found to have an improvement in accuracy for the best performing functional (B97-1).

Simulating Third-Order Nonlinear Optical Properties Using Damped Cubic Response Theory within Time-Dependent Density Functional Theory.

A general implementation for damped cubic response properties is presented using time-dependent density functional theory (TDDFT) and Slater-type orbital basis sets. To directly calculate two-photon absorption (TPA) cross sections, we also present an implementation of a reduced damped cubic response approach. Validation of the implementations includes a detailed comparison between response theory and the sum-over-states approach for calculating the nonlinear optical properties of LiH, as well as a comparison between the simulated and experimental TPA and third-harmonic generation (THG) spectra for the dimethylamino-nitrostilbene (DANS) molecule. The study of LiH demonstrates the incorrect pole structure obtained in response theory due to the adiabatic approximation typically employed for the exchange-correlation kernel. For DANS, we find reasonable agreement between simulated and experimental TPA and THG spectra. Overall, this work shows that care must be taken when calculating higher-order response functions in the vicinity of one-photon poles due to the approximate kernels typically used in the simulations.

Reexamination of the calculation of two-center, two-electron integrals over Slater-type orbitals. III. Case study of the beryllium dimer

In this paper we present results of ab initio calculations for the beryllium dimer with a basis set of Slater-type orbitals (STOs). Nonrelativistic interaction energy of the system is determined using the frozen-core full configuration interaction calculations combined with high-level coupled-cluster correction for inner-shell effects. We have developed STO basis sets, ranging in quality from double to sextuple $\ensuremath{\zeta}$, which are used in these computations. Principles of their construction are discussed and several atomic benchmarks are presented. Relativistic effects of order ${\ensuremath{\alpha}}^{2}$ are calculated perturbatively by using the Breit-Pauli Hamiltonian and are found to be significant. We also estimate the leading-order QED effects. Influence of the adiabatic correction is found to be negligible. Finally, the electronic binding energy of the beryllium dimer is determined to be $929.0\phantom{\rule{0.16em}{0ex}}\ifmmode\pm\else\textpm\fi{}\phantom{\rule{0.16em}{0ex}}1.9\phantom{\rule{0.16em}{0ex}}{\text{cm}}^{\ensuremath{-}1}$, in a very good agreement with the recent experimental value.

The assessment and application of an approach to noncovalent interactions: the energy decomposition analysis (EDA) in combination with DFT of revised dispersion correction (DFT-D3) with Slater-type orbital (STO) basis set

An assessment study is presented about energy decomposition analysis (EDA) in combination with DFT including revised dispersion correction (DFT-D3) with Slater-type orbital (STO) basis set. There has been little knowledge about the performance of the EDA + DFT-D3 concerning STOs. In this assessment such an approach was applied to calculate noncovalent interaction energies and their corresponding components. Complexes in S22 set were used to evaluate the performance of EDA in conjunction with four representative types of GGA-functionals of DFT-D3 (BP86-D3, BLYP-D3, PBE-D3 and SSB-D3) with three STO basis sets ranging in complexity from DZP, TZ2P to QZ4P. The results showed that the approach of EDA + BLYP-D3/TZ2P has a better performance not only in terms of calculating noncovalent interaction energy quantitatively but also in analyzing corresponding energy components qualitatively. This approach (EDA + BLYP-D3/TZ2P) was thus applied further to two representative large-system complexes including porphine dimers and fullerene aggregates to gain a better insight into binding characteristics.

Open-system electronic dynamics and thermalized electronic structure

We propose and implement a novel computational method for simulating open-system electronic dynamics and obtaining thermalized electronic structures within an open quantum system framework. The system-bath interaction equation of motion is derived and modeled from the local harmonic oscillator description for electronic density change. The nonequilibrium electronic dynamics in a thermal bath is simulated using first-order kinetics. The resultant electronic densities are temperature-dependent and can take characteristics of the ground and excited states. We present results of calculations performed on H(2) and 1,3-butadiene performed at the Hartree-Fock level of theory using a minimal Slater-type orbital basis set.

Cusps and derivatives for wave‐functions expanded in Slater orbitals: A density study

AbstractThis article recalls basic properties of Slater type Orbital (STO) basis sets, in particular, with respect to satisfying Kato's electron‐nucleus cusp condition. It is shown that a suitable starting point is the set of hydrogen‐like orbitals or, more flexibly, their generalization as Coulomb Sturmians. A case study is presented on the Hartree–Fock density obtained for a water molecule, N2, and CO2 showing that the nuclear cusp condition may be easily accounted for near the heavier nuclei, but only with difficulties at the hydrogen center. In particular, bond stretching affects the expected cusp condition on hydrogen atoms. Some indications are given on the improvement of STO basis sets. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009

Importance of the Basis Set for the Spin-State Energetics of Iron Complexes

We have performed a systematic investigation of the influence of the basis set on relative spin-state energies for a number of iron compounds. In principle, with an infinitely large basis set, both Slater-type orbital (STO) and Gaussian-type orbital (GTO) series should converge to the same final answer, which is indeed what we observe for both vertical and relaxed spin-state splittings. However, we see throughout the paper that the STO basis sets give consistent and rapidly converging results, while the convergence with respect to the basis set size is much slower for the GTO basis sets. For example, the large GTO basis sets that give good results for the vertical spin-state splittings of compounds 1-3 (6-311+G**, Ahlrichs VTZ2D2P) fail for the relaxed spin-state splittings of compound 4 (where 1 is Fe-(PyPepS)2 (PyPepSH 2 = N-(2-mercaptophenyl)-2-pyridinecarboxamide), 2 is Fe(tsalen)Cl (tsalen = N, N'-ethylenebis-(thiosalicylideneiminato)), 3 is Fe(N(CH2-o-C6H4S) 3)(1-Me-imidazole), and 4 is FeFHOH). Very demanding GTO basis sets like Dunning's correlation-consistent (cc-pVTZ, cc-pVQZ) basis sets are needed to achieve good results for these relaxed spin states. The use of popular (Pople-type) GTO, effective core potentials basis set (ECPB), or mixed ECPB(Fe):GTO(rest) basis sets is shown to lead to substantial deviations (2-10 kcal/mol, 14-24 kcal/mol for 3-21G), in particular for the high spin states that are typically placed at too low energy. Moreover, the use of an effective core potential in the ECPB basis sets results in spin-state splittings that are systematically different from the STO-GTO results.

Augmenting basis set for time-dependent density functional theory calculation of excitation energies: Slater-type orbitals for hydrogen to krypton

Attention is given to the dual nature of the time-dependent density functional theory approach for predicting excitation spectra, namely its theoretical basis in dynamic polarizability and its practical implementation in computer programs as a single configuration interaction. A procedure for generating diffuse functions to be added to standard Slater-type orbital basis sets for H to Kr is proposed and tested on ten close-shell molecules. The database for comparing the performance of standard and augmented basis sets consists of the 15 lowest transitions to valence and Rydberg-like states for each of the molecules. The results of this computational study are very encouraging. The addition of new augmenting functions improves the performance of standard basis sets significantly. For example, the new augmented TZP basis set, about the same size as the standard TZ2P set (and considerably smaller than the large QZ4P set), led to an average absolute deviation of 150 predicted excitation energies of only 0.12 eV from those obtained with the very large basis set called QZ3P-3DIFFUSE, compared to 0.83, 0.79, and 0.19 eV, for the standard TZP, TZ2P, and QZ4P sets, respectively. Similar augmenting functions for Gaussian-type orbital basis sets for H to Kr are also suggested.

A sparse algorithm for the evaluation of the local energy in quantum Monte Carlo

A new algorithm is presented for the sparse representation and evaluation of Slater determinants in the quantum Monte Carlo (QMC) method. The approach, combined with the use of localized orbitals in a Slater-type orbital basis set, significantly extends the size molecule that can be treated with the QMC method. Application of the algorithm to systems containing up to 390 electrons confirms that the cost of evaluating the Slater determinant scales linearly with system size.

Density-functional generalized-gradient and hybrid calculations of electromagnetic properties using Slater basis sets.

In this paper we extend our density-functional theory calculations, with generalized gradient approximation and hybrid functionals, using Slater-type orbitals (STOs), to the determination of second-order molecular properties. The key to the entire methodology involves the fitting of all STO basis function products to an auxiliary STO basis, through the minimization of electron-repulsion integrals. The selected properties are (i) dipole polarizabilities, (ii) nuclear magnetic shielding constants, and (iii) nuclear spin-spin coupling constants. In all cases the one-electron integrals involving STOs were evaluated by quadrature. The implementation for (ii) involved some complexity because we used gauge-including atomic orbitals. The presence of two-electron integrals on the right-hand side of the coupled equations meant that the fitting procedure had to be implemented. For (iii) in the hybrid case, fitting procedures were again required for the exchange contributions. For each property we studied a number of small molecules. We first obtained an estimate of the basis set limit using Gaussian-type orbitals (GTOs). We then showed how it is possible to reproduce these values using a STO basis set. For (ii) a regular TZ2P quality STO basis was adequate; for (i) the addition of one set of diffuse functions (determined by Slater's rules) gave the required accuracy; for (iii) it was necessary to add a set of 1s functions, including one very tight function, to give the desired result. In summary, we show that it is possible to predict second-order molecular properties using STO basis sets with an accuracy comparable with large GTO basis sets. We did not encounter any major difficulties with either the selection of the bases or the implementation of the procedures. Although the energy code (especially in the hybrid case) may not be competitive with a regular GTO code, for properties we find that STOs are more attractive.

Even-tempered Slater-type orbitals revisited: from hydrogen to krypton.

Even-tempered Slater-type orbital basis sets were developed in 1973, based on total atomic energy optimization. Here, we revisit ET STOs and propose new sets based on past experience and recent computational studies. From preliminary atomic and molecular tests, these sets are shown to be very well balanced and to perform, at lower cost, almost as well as a very large (close to complete) basis set.

Electric field integrals for Slater‐type orbitals

AbstractGeneral algorithms for electric field integrals with Slater‐type orbitals (STO) are reported, and the efficiency of the computational procedure based in these algorithms is analyzed. An accuracy of nine decimal places, at least, is estimated for the molecular electric field computed on any point. Timing was examined for calculations with STO basis sets ranging from 30 to 300 basis functions. Computational cost was moderate, even for the largest basis set. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004

DFT calculation of core-electron binding energies

A total of 59 core-electron binding energies (CEBEs) were studied with the Amsterdam Density Functional Program (ADF) program and compared with the observed values. The results indicate that a polarized triple-zeta basis set of Slater-type orbitals is adequate for routine assessment of the performance of each method of computation. With such a basis set, seven density functionals were tested. In addition, the performance of 21 energy density functionals were computed from the density calculated with the statistical average of orbital potentials (SAOP). Among all the choices tested, the best density functional for core-electron binding energies of C to F turns out to be the combination of Perdew-Wang (1986) functional for exchange and the Perdew-Wang (1991) functional for correlation, confirming earlier studies based on contracted Gaussian-type orbitals. For this best functional, five Slater-type orbital basis sets were examined, ranging from polarized double-zeta quality to the largest set available in the ADF package. For the best functional with the best basis set, the average absolute deviation (AAD) of the calculated value from experiment is only 0.16 eV.