Relativistic electronic-structure calculations employing a two-component no-pair formalism with external-field projection operators
TL;DR: It is found that correlation contribution and relativistic effects are nonadditive in the MRD-CI method.
Abstract: A no-pair formalism employing external-field projection operators correct to second order in the potential is used to calculate the 1s energies of one-electron atoms and ground-state properties of the bromine and silver atoms in the framework of the multireference double-excitation configuration-interaction (MRD-CI) method. It is found that the relativistic two-component method that has been used reproduces the one-particle energies of the Dirac equation to order (Z\ensuremath{\alpha}${)}^{3}$. The operator is bounded from below and can be used variationally in relativistic electron-structure calculations of many-electron atoms and molecules. The relativistic correction to the total energy recovers 97% of the relativistic correction of the Dirac-Hartree-Fock (DHF) result in the case of the bromine atom and more than 99% in the case of the silver atom. The relativistic correction of the ionization potential of silver has been calculated to be 0.47 eV at the CI level, in good agreement with DHF results, the correlation contribution in the relativistic case being 0.42 eV. The remaining discrepancy of the absolute value of 6.85 eV (DHF 6.34 eV) to experiment (7.57 eV) is attributed to basis-set deficiencies. The corresponding CI value of the electron affinity (relativistic CI value 1.05 eV, nonrelativistic 0.90 eV) is in much better agreement with experiment (1.30 eV). It is found that correlation contribution and relativistic effects are nonadditive.
Citations
More filters
TL;DR: An overview of NWChem is provided focusing primarily on the core theoretical modules provided by the code and their parallel performance, as well as Scalable parallel implementations and modular software design enable efficient utilization of current computational architectures.
4,666 citations
TL;DR: In this article, potential-dependent transformations are used to transform the four-component Dirac Hamiltonian to effective two-component regular Hamiltonians, which already contain the most important relativistic effects, including spin-orbit coupling.
Abstract: In this paper, potential‐dependent transformations are used to transform the four‐component Dirac Hamiltonian to effective two‐component regular Hamiltonians. To zeroth order, the expansions give second order differential equations (just like the Schrodinger equation), which already contain the most important relativistic effects, including spin–orbit coupling. One of the zero order Hamiltonians is identical to the one obtained earlier by Chang, Pelissier, and Durand [Phys. Scr. 34, 394 (1986)]. Self‐consistent all‐electron and frozen‐core calculations are performed as well as first order perturbation calculations for the case of the uranium atom using these Hamiltonians. They give very accurate results, especially for the one‐electron energies and densities of the valence orbitals.
3,585 citations
TL;DR: Seven different types of Slater type basis sets for the elements H (Z = 1) up to E118, ranging from a double zeta valence quality up to a quadruple zetavalence quality, are tested in their performance in neutral atomic and diatomic oxide calculations.
Abstract: 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.
2,112 citations
TL;DR: In this paper, the energy gradients in the zeroth order regular approximation (ZORA) to the Dirac equation were derived for the transition metal complexes W(CO), Os(CO)5, and Pt (CO)4.
Abstract: Analytical expressions are derived for the evaluation of energy gradients in the zeroth order regular approximation (ZORA) to the Dirac equation. The electrostatic shift approximation is used to avoid gauge dependence problems. Comparison is made to the quasirelativistic Pauli method, the limitations of which are highlighted. The structures and first metal-carbonyl bond dissociation energies for the transition metal complexes W(CO)6, Os(CO)5, and Pt(CO)4 are calculated, and basis set effects are investigated.
2,055 citations
TL;DR: HSE is a fast and accurate alternative to established density functionals, especially for solid state calculations, and correctly predicts semiconducting behavior in systems where pure functionals erroneously predict a metal, such as, for instance, Ge.
Abstract: This work assesses the Heyd-Scuseria-Ernzerhof (HSE) screened Coulomb hybrid density functional for the prediction of lattice constants and band gaps using a set of 40 simple and binary semiconductors. An extensive analysis of both basis set and relativistic effects is given. Results are compared with established pure density functionals. For lattice constants, HSE outperforms local spin-density approximation (LSDA) with a mean absolute error (MAE) of 0.037 A for HSE vs 0.047 A for LSDA. For this specific test set, all pure functionals tested produce MAEs for band gaps of 1.0–1.3 eV, consistent with the very well-known fact that pure functionals severely underestimate this property. On the other hand, HSE yields a MAE smaller than 0.3 eV. Importantly, HSE correctly predicts semiconducting behavior in systems where pure functionals erroneously predict a metal, such as, for instance, Ge. The short-range nature of the exchange integrals involved in HSE calculations makes their computation notably faster than regular hybrid functionals. The current results, paired with earlier work, suggest that HSE is a fast and accurate alternative to established density functionals, especially for solid state calculations.
1,564 citations