Accuracy and efficiency of atomic basis set methods versus plane wave calculations with ultrasoft pseudopotentials for DNA base molecules

Recent results from Preuss et al. (J Comput Chem 2004, 25, 112) on DNA base molecules, obtained by plane wave density functional calculations using ultrasoft pseudopotentials, are compared with calculations using Gaussian basis sets. Bond lengths and angles agree closely, but dihedral angles and vibrational frequencies show significant differences. The Gaussian basis calculations are at least an order of magnitude more efficient than the plane wave/ultrasoft pseudopotential calculations at a similar level of accuracy; the advantage is even larger if the Fourier Transform Coulomb method is used. To obtain definite benchmark values, the geometries of the four DNA bases were optimized at the MP2 level with large basis sets, up to cc-pVQZ and aug-cc-pVTZ.     

Atomic energy levels from configuration interaction calculations with relativistic corrections

It is shown that configuration interaction calculations, with inclusion of the relativistic corrections, constitute an appropriate approach for the prediction of atomic energy levels and that results of experimental accuracy are possible given the availability of large-scale, fast computers. The results obtained for He through F emphasize both the practical difficulties to be encountered and the possibility of predictions with less than 1% error.     

Relativistic effects determined using the Douglas-Kroll contracted basis sets and correlation consistent basis sets with small-core relativistic pseudopotentials

The coupled cluster approximation with single, double, and quasiperturbative triple excitations [CCSD(T)] was used in combination with the Douglas-Kroll contracted correlation consistent basis sets [cc-pVnZ-DK, where n = D(2), T(3), Q(4), and 5] and small-core relativistic pseudopotentials (PP) with correlation consistent polarized valence basis sets (cc-pVnZ-PP and aug-cc-pVnZ-PP) to investigate the impact of scalar relativistic corrections on energetic and structural properties of small molecules containing third-row (Ga-Kr) atoms. These molecules were taken from the Gaussian-2 extended test set for third-row atoms. Atomization energies, ionization energies, electron affinities, and proton affinities for molecules in the test set were determined and compared with nonrelativistic results which were obtained in a recent study in which the standard and augmented correlation consistent basis sets were used in combination with CCSD(T). Several schemes were used to extrapolate the energies to the complete basis set limit.     

Reliability of atomic natural orbital basis sets in calculations involving pseudopotentials

The performance of Atomic Natural Orbital (ANO) basis sets for calculations involving nonempirical core pseudopotentials has been studied by comparing the results for atomic and molecular nitrogen obtained using contracted ANO basis sets with those obtained using both the primitive set and a segmented one. The primitive set has been optimized at the SCF level for atomic N treated as a five-electron pseudo-atom, and consists of 7s and 7p primitive GTOs supplemented by 2d and 1f GTOs optimized at the CI level. From this primitive set three contracted [3s 3p 2d 1f] sets have been obtained. The first one has been derived from the ANOs of the neutral atom, the second has been obtained from an averaged density matrix and the third one is a segmented set. For the atom, the segmented set gives a zero contraction error at the SCF level as it must be in valence-only calculations. The ANO basis sets show some small contraction error at the SCF level but perform better in CI calculations. However, for the diatomic N₂ molecule the ANO basis sets exhibit a rather large contraction error in the calculated SCF energy. A detailed analysis of the origin of this error is reported, which shows that the conventional strategy used to derive ANO basis sets does not work very well when pseudopotentials are involved.     

On ab-initio calculations of singly excited states of atoms with augmented fuess-type basis functions

It is suggested that by simply augmenting a given optimized basis set for the ground state wavefunction of an atom by Fuess-type functions, quality wavefunctions of the singly excited states of the atom could be obtained. The non-linear parameters of the Fuess-type functions are determined from the known quantum defect data. As an illustrative example we present the results for the singly excited ¹S and ¹P states of He. The calculated energies are within 0.2% of the experimental values. The dipole oscillator strengths are also calculated from both the length and the velocity formulas. The results compare well with the experimental values; the velocity form giving generally slightly closer agreement with observations.     

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.     

Highly accurate relativistic universal Gaussian basis set for Dirac–Fock–Breit calculations

Ion mobilities of H₂O⁺ drifting in helium are calculated and compared with experiment. These calculations employ global potential energy surfaces of the H₂O⁺–He complex, which in the present case were calculated ab initio at the unrestricted MP2 level of theory using a basis set of aug‐cc‐pVTZ quality, and treating the ion as a rigid body. Details are presented of the general characteristics of both the ground and first‐excited electronic states of the complex. Although only the ground‐state surface was used for the mobility calculations, the ab initio determination of the ground state necessitated the inclusion of the first‐excited state owing to the presence of a crossing between the two. This crossing is also described. Mobilities calculated from the global surfaces are in good agreement with experiment. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Fully relativistic pseudopotentials for alkaline atoms: Dirac-Hartree-Fock and configuration interaction calculations of alkaline monohydrides

Fully relativistic four-component energy-adjusted pseudopotentials and corresponding valence basis sets have been derived for the alkaline atoms Li through Cs, treating them as one-valence electron systems. Core-valence correlation effects are accounted for by a core-polarization potential, deviations of the core-nucleus repulsion from a point charge model by a suitable correction. The results of Dirac-Hartree-Fock and configuration interaction calculations are presented for atomic properties not used in the pseudopotential adjustment, i.e. electron affinities and dipole polarizabilities, as well as for the spectroscopic constants of the ground states of the alkaline monohydrides. The analytic form of the cut-off function for the electric field in the core-polarization term and its effects on atomic and molecular properties is discussed.     

Simplified ab-initio calculations for molecular systems

A method is described for reducing a large part of the arithmetic of exact ab-initio SCF molecularorbital calculations based on Slater-type-orbitals without noticeable loss of numerical accuracy. The procedure involves the transformation to Löwdin orthogonalized orbitals and then invoking the NDDO approximation. The three- and four-centre two-electron integrals required are estimated by a truncated Ruedenberg expansion. All one-electron integrals are evaluated exactly. No empirical parameters are employed. Numerical tests on CO, OF₂, O₃ and ONF show that the NDDO approximation is very accurate for Löwdin functions and that the Ruedenberg expansion is arithmetically satisfactory for the SCF MO calculations.     

Simplified ab-initio calculations on hydrogen-containing molecules

The simplified ab-initio method described in an earlier paper is tested on some hydrogen-containing molecules. The performance is slightly below that found previously for molecules composed entirely of first-row atoms but should be suitable for applications where limited numerical accuracy is sufficient. The hope of improved performance through limited expansion of the basis, especially on hydrogen, is not realised and so alternative treatments of the two-electron many-centre integrals should be sought if greater numerical accuracy is required.     

Ab initio CI calculations on benzene with an extended basis set

An extended basis set of triple zeta plus polarization quality is employed to carry out configuration interaction (CI ) calculations of the three lowest singlet and triplet excited states of benzene. The CI calculation is carried out by taking into account single and double excitations of π and σ electrons. In the CI , composite natural orbitals (CNO s), which are constructed from the natural orbitals of the ground state of ethylene, are used as virtual orbitals. The aim of using CNOs is to reduce the number of virtual orbitals to be used in constructing configuration-state functions, thus cutting down CI dimensions without losing reasonable accuracy. The excitation energies resulting from the CI are in fairly good agreement with experiment. The root mean square of the deviation is 0.22 eV for the six calculated energies and the largest disagreement is 0.37 eV for the third singlet excited state. To obtain better excitation energies by an ab initio calculation, it seems likely that we need to take into account more electron correlation than in the present calculation. © 1994 John Wiley & Sons, Inc.     

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.     

Energy-consistent pseudopotentials for quantum Monte Carlo calculations

The authors present scalar-relativistic energy-consistent Hartree-Fock pseudopotentials for the main-group elements. The pseudopotentials do not exhibit a singularity at the nucleus and are therefore suitable for quantum Monte Carlo (QMC) calculations. They demonstrate their transferability through extensive benchmark calculations of atomic excitation spectra as well as molecular properties. In particular, they compute the vibrational frequencies and binding energies of 26 first- and second-row diatomic molecules using post-Hartree-Fock methods, finding excellent agreement with the corresponding all-electron values. They also show their pseudopotentials give superior accuracy than other existing pseudopotentials constructed specifically for QMC. Finally, valence basis sets of different sizes (VnZ with n=D,T,Q,5 for first and second rows, and n=D,T for third to fifth rows) optimized for our pseudopotentials are also presented.     

Improved relativistic energy-consistent pseudopotentials for 3d-transition metals

Energy-consistent relativistic pseudopotentials for 3d-transition metals Sc to Ni based on modified valence energies are proposed. The pseudopotentials are adjusted at the finite difference level within the intermediate coupling scheme with respect to multi-configuration Dirac–Hartree–Fock data based on the Dirac–Coulomb Hamiltonian with an estimate of the Breit contributions in quasidegenerate perturbation theory. Typically a few hundred to thousand J levels arising from about 35 to 40 configurations ranging from the anion down to the highly charged cation are considered as references. It is shown that introducing a small common energetic shift of all valence energies reduces the errors in the parameter adjustment considerably. Results of highly correlated atomic and molecular test calculations using large basis sets and basis set extrapolation techniques are presented. To be submitted to Theoretical Chemistry Accounts (special volume on the occasion of Prof. Dr. H. Stoll's 60th birthday)     

The ORP basis set designed for optical rotation calculations

Details of generation of the optical rotation prediction (ORP) basis set developed for accurate optical rotation (OR) calculations are presented. Specific rotation calculations carried out at the density functional theory (DFT) level for model chiral methane molecule, fluorooxirane, methyloxirane, and dimethylmethylenecyclopropane reveal that the ORP set outperforms larger basis sets, among them the aug‐cc‐pVTZ basis set of Dunning (J. Chem. Phys. 1989, 90, 1007) and the aug‐pc‐2 basis set of Jensen (J. Chem. Phys. 2002, 117, 9234; J. Chem. Theory Comput. 2008, 4, 719). It is shown to be an attractive choice also in the case of larger systems, namely norbornanone, β‐pinene, trans‐pinane, and nopinone. The ORP basis set is further used in OR calculations for 24 other systems, and the results are compared to the aug‐cc‐pVDZ values. Whenever large discrepancies of results are observed, the ORP values are in an excellent agreement with the aug‐cc‐pVTZ results. The ORP basis set enables accurate specific rotation calculations at a reduced cost and thus can be recommended for routine DFT OR calculations, also for large and conformationally flexible molecules. © 2013 Wiley Periodicals, Inc.     

Ab initio calculations of relativistic and electron correlation effects in polyatomics using the universal Gaussian basis set: XeF2

Ab initio accurate all-electron relativistic molecular orbital Dirac–Fock self-consistent field calculations are reported for the linear symmetric XeF₂ molecule at various internuclear distances with our recently developed relativistic universal Gaussian basis set. The nonrelativistic limit Hartree–Fock calculations were also performed for XeF₂ at various internuclear distances. The relativistic correction to the electronic energy of XeF₂ was calculated as ～ ?215 hartrees (?5850 eV) by using the Dirac–Fock method. The dominant magnetic part of the Breit interaction correction to the nonrelativistic interelectron Coulomb repulsion was included in our calculations by both the Dirac–Fock–Breit self-consistent field and perturbation methods. The calculated Breit correction is ～6.5 hartrees (177 eV) for XeF₂. The relativistic Dirac–Fock as well as the nonrelativistic HF wave functions predict XeF₂ to be unbound, due to neglect of electron correlation effects. These effects were incorporated for XeF₂ by using various ab initio post Hartree–Fock methods. The calculated dissociation energy obtained using the MP 2(full) method with our extensive basis set of 313 primitive Gaussians that included d and f polarization functions on Xe and F is 2.77 eV, whereas the experimental dissociation energy is 2.78 eV. The calculated correlation energy is ～ ?2 hartrees (?54 eV) at the predicted internuclear distance of 1.986 Å, which is in excellent agreement with the experimental Xe—F distance of 1.979 Å in XeF₂. In summary, electron correlation effects must be included in accurate ab initio calculations since it has been shown here that their inclusion is crucial for obtaining theoretical dissociation energy (D_e) close to experimental value for XeF₂. Furthermore, relativistic effects have been shown to make an extremely significant contribution to the total energy and orbital binding energies of XeF₂. © 1995 John Wiley & Sons, Inc.     

On the use of contracted basis functions in relativistic Hartree-Fock calculations

A procedure is suggested to build up contracted basis sets for relativistic atomic and molecular Hartree-Fock calculations when corresponding non-relativistic results are available or easy to obtain. Significant reductions of the size of relativistic Fock-Roothaan matrices are expected.     

A finite B-spline basis set for accurate diatomic molecule calculations

A finite basis set particularly adapted for solving the Hartree-Fock equation for diatomic molecules in prolate spheroidal coordinates has been constructed. These basis functions have been devised as products of B-splines times associated Legendre polynomials. Due to the large number of B-splines, the resulting set of eigenfunctions is amply distributed over excited states. This gives the possibility of using these basis sets to calculate sums over excited states, appearing in various orders of perturbation theory. As an illustration, the second-order corrections to the ground-state energy of some atoms and diatomic molecules with closed electron shells have been calculated.