Accuracy of basis-set extrapolation schemes for DFT-RPA correlation energies in molecular calculations

We construct a reference benchmark set for atomic and molecular random phase approximation (RPA) correlation energies in a density functional theory framework at the complete basis-set limit. This set is used to evaluate the accuracy of some popular extrapolation schemes for RPA all-electron molecular calculations. The results indicate that for absolute energies, accurate results, clearly outperforming raw data, are achievable with two-point extrapolation schemes based on quintuple- and sextuple-zeta basis sets. Moreover, we show that results in good agreement with the benchmark can also be obtained by using a semiempirical extrapolation procedure based on quadruple- and quintuple-zeta basis sets. Finally, we analyze the performance of different extrapolation schemes for atomization energies.     

Effect of double bonds on the conducting properties of ciguatoxin 3C and tetrahydropyrane-based polymers: a theoretical study

The electronic structure of the ciguatoxin 3C is analyzed through the Kohn-Sham model by using two different kinds of basis sets: localized basis set (Gaussian functions) and nonlocalized basis set (plane wave functions). With the localized basis functions, two approximations are used for the exchange-correlation functional: the local density approximation and the generalized gradient approximation. With the nonlocalized basis set, just the local density approximation is used. The energy gap, obtained from the frontier molecular orbitals, for this molecule predicts that this system is a semiconductor, even when the number of double bonds is increased inside the structure. However, as large molecules built with the basic unit--the tetrahydropyrane--of the ciguatoxin 3C are found in nature, it suggests studying the gap in polymeric systems built with the basic unit of this molecule. It is demonstrated that the presence of double bonds reduces considerably the gap, indicating the possibility of forming conducting materials by introducing double bonds in this kind of molecular systems. Thus, molecules strongly linked with biological systems can be used as precursor to build electric conducting systems.     

Basis set convergence of molecular correlation energy differences within the random phase approximation

The basis set convergence of energy differences obtained from the random phase approximation (RPA) to the correlation energy is investigated for a wide range of molecular interactions. For dispersion bound systems the basis set incompleteness error is most pronounced, as shown for the S22 benchmark [P. Jurecka et al., Phys. Chem. Chem. Phys. 8, 1985 (2006)]. The use of very large basis sets (> quintuple-zeta) or extrapolation to the complete basis set (CBS) limit is necessary to obtain a reliable estimate of the binding energy for these systems. Counterpoise corrected results converge to the same CBS limit, but counterpoise correction without extrapolation is insufficient. Core-valence correlations do not play a significant role. For medium- and short-range correlation, quadruple-zeta results are essentially converged, as demonstrated for relative alkane conformer energies, reaction energies dominated by intramolecular dispersion, isomerization energies, and reaction energies of small organic molecules. Except for weakly bound systems, diffuse augmentation almost universally slows down basis set convergence. For most RPA applications, quadruple-zeta valence basis sets offer a good balance between accuracy and efficiency.     

Effects of different basis sets and donor‐acceptor groups on linear and second‐order nonlinear optical properties and molecular frontier orbital energies

The calculation of molecular hyperpolarizability, molecular frontier orbital energies of some donor‐acceptor oxadiazoles ( 5a – f , 8a – f , and 9a – f ) have been investigated using ab initio methods and different basis sets. Ab initio optimizations were performed at the Hartree–Fock (HF) and density functional (Beckee‐3–Lee–Yang–Parr; B3LYP) levels of theory with 6‐31G basis set. The polarizability (<α>), anisotropy of polarizability (Δα), and ground‐state dipole moment (μ), first hyperpolarizability (β), and molecular frontier orbital (HOMO, highest occupied molecular orbital and LUMO, lowest unoccupied molecular orbital) energies of 5a – f , 8a – f , and 9a – f have been calculated at the HF and B3LYP methods with 6‐31G, 6‐31G(d), 6‐31+G(d), 6‐31++G(d,p), 6‐311G, 6‐311G(d), 6‐311+G(d), and 6‐311++G(d,p) basis sets. Also, the molecular hardness (η) and electronegativity (χ) parameters have been obtained using molecular frontier orbital energies. The <α>, Δα, μ, β, HOMO, LUMO energies, η and χ parameters have been investigated as dependence on the choice of method and basis set. The variation graphics of <α>, Δα, μ, β, η, and χ parameters using HF and B3LYP methods with different basis sets are presented. We have examined the frontier molecular orbital pictures of 5a – f , 8a – f , and 9a – f using B3LYP/6‐31++G(d,p) level. The 5a – f , 8a – f , and 9a – f display significant linear, second‐order molecular nonlinearity, and molecular parameters and provide the basis for future design of efficient nonlinear optical materials having the 1,3,4‐oxadiazole core. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

General orbital invariant MP2-F12 theory

A general form of orbital invariant explicitly correlated second-order closed-shell Moller-Plesset perturbation theory (MP2-F12) is derived, and compact working equations are presented. Many-electron integrals are avoided by resolution of the identity (RI) approximations using the complementary auxiliary basis set approach. A hierarchy of well defined levels of approximation is introduced, differing from the exact theory by the neglect of terms involving matrix elements over the Fock operator. The most accurate method is denoted as MP2-F12/3B. This assumes only that Fock matrix elements between occupied orbitals and orbitals outside the auxiliary basis set are negligible. For the chosen ansatz for the first-order wave function this is exact if the auxiliary basis is complete. In the next lower approximation it is assumed that the occupied orbital space is closed under action of the Fock operator [generalized Brillouin condition (GBC)]; this is equivalent to approximation 2B of Klopper and Samson [J. Chem. Phys. 116, 6397 (2002)]. Further approximations can be introduced by assuming the extended Brillouin condition (EBC) or by neglecting certain terms involving the exchange operator. A new approximation MP2-F12/3C, which is closely related to the MP2-R12/C method recently proposed by Kedzuch et al. [Int. J. Quantum Chem. 105, 929 (2005)] is described. In the limit of a complete RI basis this method is equivalent to MP2-F12/3B. The effect of the various approximations (GBC, EBC, and exchange) is tested by studying the convergence of the correlation energies with respect to the atomic orbital and auxiliary basis sets for 21 molecules. The accuracy of relative energies is demonstrated for 16 chemical reactions. Approximation 3C is found to perform equally well as the computationally more demanding approximation 3B. The reaction energies obtained with smaller basis sets are found to be most accurate if the orbital-variant diagonal Ansatz combined with localized orbitals is used for the first-order wave function. This unexpected result is attributed to geminal basis set superposition errors present in the formally more rigorous orbital invariant methods.     

Compact basis sets for LCAO-LSD calculations. Part II: Tests for Cr2 and Ni4

The compact orbital and auxiliary basis sets for LCAO-LSD calculations introduced in Part I are tested in molecular calculations on Cr₂ and Ni₄. The present results for spectroscopic constants and valence orbital energies obtained using medium size orbital expansions with a double-zeta representation for valence orbitals are in very good agreement with those previously calculated with very extended sets. Since the computational time of the present calculations is reduced severalfold compared with the extended basis set calculations, the present basis sets allow increased efficiency of the LCAO-LSD calculations and allow the method to be extended to larger systems.     

Electronic-structure studies of solides. V. Rigorous Hartree-Fock treatment of metallic hydrogen using a plane-wave basis

We have performed Hartree-Fock calculations for simple cubic metallic hydrogen crystals using Bloch functions expanded in plane waves, All integrals were evaluated accurately including exchangeterms. Increasingly larger basis sets were used, and the total Hartree-Fock energy obtained with the maximum number of plane waves (27) was ?0.4770 hartrees/atom. This total energy is believed to be within a few millihartress of the Hartree-Fock limit results. The deficiency of a plane-wave expansion to represent the atomic cusps, however, makes it difficult to obtain the exact Hartree-Fock limit with a plane-wave expansion. When the correlation energy (calculated in the random-phase approximation with Hartree-Fock bands and functions as zeroth order states) is added, and upper limit of ?0.501 hartrees/atom is found for the total energy of this system. The Fermi surface was found to touch the Brillouin zone boundaries around the X points due to an appreciable depression of the band energies in that part of the Brillouin zone. The equilibrium lattice spacing (a = 2.705 bohrs) was slightly smaller than that obtained earlier with an atomic orbital basis.     

Accurate calculations of intermolecular interaction energies using explicitly correlated wave functions

Explicitly correlated second-order M?ller-Plesset (MP2-F12) calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented and compared with standard MP2 results. The MP2 complete basis set limits are estimated using basis set extrapolation and augmented quadruple-zeta and quintuple-zeta basis sets. Already with augmented double-zeta basis sets the MP2-F12 interaction energies are found to be closer to the complete basis set limits than standard MP2 calculations with augmented quintuple-zeta basis sets. Various possible approximations in the MP2-F12 method are systematically tested. Best results are obtained with localized orbitals and the diagonal MP2-F12/C(D) ansatz. Hybrid approximations, in which some contributions of the auxiliary basis set are neglected and which considerably reduce the computational cost, have a negligible effect on the interaction energies. Also the orbital-invariant fixed-amplitude approximation of Ten-no leads to only slightly less accurate results. Preliminary results for the neon and benzene dimers, obtained with the recently proposed CCSD(T)-F12a approximation, indicate that the CCSD(T) basis set limits can also be very closely approached using augmented triple-zeta basis sets.     

Approaching the theoretical limit in periodic local MP2 calculations with atomic-orbital basis sets: the case of LiH

The atomic orbital basis set limit is approached in periodic correlated calculations for solid LiH. The valence correlation energy is evaluated at the level of the local periodic second order M?ller-Plesset perturbation theory (MP2), using basis sets of progressively increasing size, and also employing "bond"-centered basis functions in addition to the standard atom-centered ones. Extended basis sets, which contain linear dependencies, are processed only at the MP2 stage via a dual basis set scheme. The local approximation (domain) error has been consistently eliminated by expanding the orbital excitation domains. As a final result, it is demonstrated that the complete basis set limit can be reached for both HF and local MP2 periodic calculations, and a general scheme is outlined for the definition of high-quality atomic-orbital basis sets for solids.     

The ellipsoidal Gaussian basis in molecular orbital theory

The ellipsoidal Gaussian basis function used in a minimal valence atomic orbital representation is compared with the double-zeta spherical Gaussian basis orbital representation for some seventeen molecules made up of first row atoms and hydrogen. Except for acetylene the double-zeta basis gives consistently better total electronic energies and generally better property values than the optimized ellipsoidal single zeta basis. Difference molecular density contour maps comparing the two basis sets, as well as other one-electron property values, indicate that the ellipsoidal basis exaggerates the transfer of charge from the atomic regions to the interatomic and lone pair regions of molecules. Apparently, the forced complete elliptization of the valence atomic orbital in the single-zeta representation does not allow the basis set sufficient flexibility to simultaneously represent both the basically spherical atomic part of these orbitals and the non-spherical molecular bond formation. Other properties and aspects of the ellipsoidal Gaussian basis are also discussed.     

Study on the maximum accuracy of the pseudopotential density functional method with localized atomic orbitals versus plane-wave basis sets

A detailed study on the accuracy attainable with numerical atomic orbitals in the context of pseudopotential first-principles density functional theory is presented. Dimers of first- and second-row elements are analyzed: bond lengths, atomization energies, and Kohn-Sham eigenvalue spectra obtained with localized orbitals and with plane-wave basis sets are compared. For each dimer, the cutoff radius, the shape, and the number of the atomic basis orbitals are varied in order to maximize the accuracy of the calculations. Optimized atomic orbitals are obtained following two routes: (i) maximization of the projection of plane wave results into atomic orbital basis sets and (ii) minimization of the total energy with respect to a set of primitive atomic orbitals as implemented in the OPENMX software package. It is found that by optimizing the numerical basis, chemical accuracy can be obtained even with a small set of orbitals.     

Hartree-Fock complete basis set limit properties for transition metal diatomics

Numerical Hartree-Fock (HF) energies accurate to at least 1 microhartree are reported for 27 diatomic transition-metal-containing species. The convergence of HF energies toward this numerical limit upon increasing the basis set size has been investigated, where standard nonrelativistic all-electron correlation consistent basis sets and augmented basis sets, developed by Balabanov and Peterson [J. Chem. Phys. 123, 064107 (2005)], were employed. Several schemes which enable the complete basis set (CBS) limit to be determined have been investigated, and the resulting energies have been compared to the numerical Hartree-Fock energies. When comparing basis set extrapolation schemes, those in the form of exponential functions perform well for our test set, with mean absolute deviations from numerical HF energies of 234 and 153 microE(h), when the CBS limit has been determined using a two-point fit as proposed by Halkier et al. [Chem. Phys. Lett. 302, 437 (1999)] on calculations of triple- and quadruple-zeta basis set qualities and calculations of quadruple- and quintuple-zeta basis set qualities, respectively. Overall, extrapolation schemes in the form of a power series are not recommended for the extrapolation of transition metal HF energies. The impact of basis set superposition error has also been examined.     

Calibration study of the CCSD(T)-F12a/b methods for C2 and small hydrocarbons

Explicitly correlated CCSD(T)-F12a/b methods combined with basis sets specifically designed for this technique have been tested for their ability to reproduce standard CCSD(T) benchmark data covering 16 small molecules composed of hydrogen and carbon. The standard method calibration set was obtained with very large one-particle basis sets, including some aug-cc-pV7Z and aug-cc-pV8Z results. Whenever possible, the molecular properties (atomization energies, structures, and harmonic frequencies) were extrapolated to the complete basis set limit in order to facilitate a direct comparison of the standard and explicitly correlated approaches without ambiguities arising from the use of different basis sets. With basis sets of triple-ζ quality or better, the F12a variant was found to overshoot the presumed basis set limit, while the F12b method converged rapidly and uniformly. Extrapolation of F12b energies to the basis set limit was found to be very effective at reproducing the best standard method atomization energies. Even extrapolations based on the small cc-pVDZ-F12/cc-pVTZ-F12 combination proved capable of a mean absolute deviation of 0.20 kcal/mol. The accuracy and simultaneous cost savings of the F12b approach are such that it should enable high quality property calculations to be performed on chemical systems that are too large for standard CCSD(T).     

A natural linear scaling coupled-cluster method

It is shown that using an appropriate localized molecular orbital (LMO) basis, one is able to calculate coupled-cluster singles and doubles (CCSD) wave functions and energies for very large systems by performing full CCSD calculations on small subunits only. This leads to a natural linear scaling coupled-cluster method (NLSCC), in which total correlation energies of extended systems are evaluated as the sum of correlation energy contributions from individual small subunits within that system. This is achieved by defining local occupied orbital correlation energies. These are quantities, which in the LMO basis become transferable between similar molecular fragments. Conventional small scale existing molecular CCSD codes are all that is needed, the local correlation effect being simply transmitted via the appropriate LMO basis. Linear scaling of electronic correlation energy calculations is thus naturally achieved using the NLSCC approach, which in principle can treat nonperiodic extended systems of infinite basis set size. Results are shown for alkanes and several polyglycine molecules and the latter compared to recent results obtained via an explicit large scale LCCSD calculation. (c) 2004 American Institute of Physics.     

Systematic regularities of molecular SCF energies

New formulae for the approximate computation of total molecular energies are developed based on ab initio calculations of n-alkanes. Their application to various kinds of molecules reveals that good expectation values for total molecular energies can be obtained by considering only the one-electron terms h _i and the nuclear repulsion energy. It is further shown that very good agreement with SCF total energies is obtained by a relationship which connects the total energy with the sum of inner-shell (core) orbital energies. The results turn out to be better than those obtained using Ruedenberg's approximation, which takes both inner-shell and valence-shell orbital energies into account.     

On the limits of highest-occupied molecular orbital driven reactions: the frontier effective-for-reaction molecular orbital concept

We carried out Hartree-Fock (HF) and density functional theory calculations for 61 compounds, the conjugated bases of carboxylic acids, phenols, and alcohols, and analyzed their acid-base behavior using molecular orbital (MO) energies and their dependence on solvent effects. Despite the well-known correlation between highest-occupied MO (HOMO) energies and pKa, we observed that HOMO energies are inadequate to describe the acid-base behavior of these compounds. Therefore, we established a criterion to identify the best frontier MO for describing pKa values and also to understand why the HOMO approach fails. The MO that fits our criterion provided very good correlations with pKa values, much better than those obtained by HOMO energies. Since they are the frontier molecular orbitals that drive the acid-base reactions in each compound, they were called frontier effective-for-reaction MOs, or FERMOs. By use of the FERMO concept, the reactions that are HOMO driven, and those that are not, can be better explained, independently from the calculation method used, as both HF and Kohn-Sham methodologies lead to the same FERMO.     

Density functional theory investigation of the polarizability and second hyperpolarizability of polydiacetylene and polybutatriene chains: treatment of exact exchange and role of correlation

The static polarizability and second hyperpolarizability of increasingly large polydiacetylene and polybutatriene (PBT) chains have been evaluated using the optimized effective potential for exact exchange (OEP-EXX) method developed by Yang and Wu [Phys. Rev. Lett. 89, 143002 (2002)], where the unknown part of the effective potential is expressed as a linear combination of Gaussian functions. Various conventional atomic orbital basis sets were employed for the exchange potential (X basis) as well as for the Kohn-Sham orbitals [molecular orbital (MO) basis]. Our results were compared to coupled-perturbed Hartree-Fock (CPHF) calculations and to ab initio correlated values obtained at various levels of approximation. It turns out that (a) small conventional basis sets are, in general, unsatisfactory for the X basis; (b) the performance of a given X basis depends on the MO basis and is generally improved when using a larger MO basis; (c) these effects are exaggerated for the second hyperpolarizability compared to the polarizability; (d) except for the second hyperpolarizability of PBT chains, using 6-311++G** for the X basis gives reasonable agreement with the CPHF results for all MO basis sets; (e) our results suggest that in the limit of a complete X basis the OEP-EXX values may approach the CPHF data; and (f) in general, the quality of a given conventional X basis degrades with the length of the oligomer, which correlates with the fact that the number of X basis functions becomes a smaller fraction of the number required to reproduce exactly the finite-basis-set Hartree-Fock energies. Linear and especially nonlinear electric field responses constitute a very stringent test for assessing the quality of functionals and potentials; appropriately tailored basis sets are needed to describe the latter. Finally, this study further highlights the importance of electron correlation effects on linear and nonlinear responses, for which correlated functionals with OEP are required.     

Basis-set extrapolation techniques for the accurate calculation of molecular equilibrium geometries using coupled-cluster theory

To reduce remaining basis-set errors in the determination of molecular equilibrium geometries, a basis-set extrapolation (BSE) scheme is suggested for the forces used in geometry optimizations. The proposed BSE scheme is based on separating the Hartree-Fock and electron-correlation contributions and uses expressions obtained by straightforward differentiation of well established extrapolation formulas for energies when using basis sets from Dunning's hierarchy of correlation-consistent basis sets. Comparison with reference data obtained at the R12 coupled-cluster level [CCSD(T)-R12] demonstrates that BSE significantly accelerates the convergence to the basis-set limit, thus leading to improvements comparable to or even better than those obtained by increasing the cardinal number in the used basis set by one. However, BSE alone is insufficient to improve agreement with experiment, even after additional consideration of inner-shell correlation and quadruple-excitation effects (mean error and standard deviation with extrapolation are -0.014 and 0.047 pm in comparison with mean error and standard deviation of -0.002 and 0.036 pm without extrapolation). Improvement is obtained only when other contributions of similar magnitude as the BSE contributions (e.g., pentuple-excitation effects and relativistic effects) are also considered. A rather large discrepancy (of the order of a few tenths of a picometer) is observed for the F(2) molecule indicating an enhanced basis-set requirement for the various contributions in this case.