Hybrid Density Functional Approximations Research Articles

Efficient all-electron hybrid density functionals for atomistic simulations beyond 10 000 atoms.

Hybrid density functional approximations (DFAs) offer compelling accuracy for abinitio electronic-structure simulations of molecules, nanosystems, and bulk materials, addressing some deficiencies of computationally cheaper, frequently used semilocal DFAs. However, the computational bottleneck of hybrid DFAs is the evaluation of the non-local exact exchange contribution, which is the limiting factor for the application of the method for large-scale simulations. In this work, we present a drastically optimized resolution-of-identity-based real-space implementation of the exact exchange evaluation for both non-periodic and periodic boundary conditions in the all-electron code FHI-aims, targeting high-performance central processing unit (CPU) compute clusters. The introduction of several new refined message passing interface (MPI) parallelization layers and shared memory arrays according to the MPI-3 standard were the key components of the optimization. We demonstrate significant improvements of memory and performance efficiency, scalability, and workload distribution, extending the reach of hybrid DFAs to simulation sizes beyond ten thousand atoms. In addition, we also compare the runtime performance of the PBE, HSE06, and PBE0 functionals. As a necessary byproduct of this work, other code parts in FHI-aims have been optimized as well, e.g., the computation of the Hartree potential and the evaluation of the force and stress components. We benchmark the performance and scaling of the hybrid DFA-based simulations for a broad range of chemical systems, including hybrid organic-inorganic perovskites, organic crystals, and ice crystals with up to 30 576 atoms (101 920 electrons described by 244 608 basis functions).

The role of the coupling matrix elements in time-dependent density functional theory on the simulation of core-level spectra of transition metal complexes

Time-dependent density functional theory (TD-DFT) stands out as an efficient tool for computing core-level spectra in large molecules, particularly transition metal complexes. However, despite their relatively moderate computational demands, TD-DFT methods can still pose challenges for typical computations involving transition metal complexes with over a thousand basis functions. In this study, we investigate the role of the Coulomb, Hartree–Fock exchange, and exchange-correlation kernel contributions to the TD-DFT coupling matrix elements when simulating core-level spectra in transition metal complexes. Our observations reveal that the exchange-correlation kernel contribution, responsible for more than 50% of the computational time in a hybrid TD-DFT calculation, surprisingly has no discernible impact on the qualitative aspects of the calculated spectra. While the Coulomb term plays a crucial role in describing L2,3 -edge spectra, its significance becomes negligible when considering K, L 1, and M4,5 edges. In contrast, the scaled Hartree–Fock exchange is demonstrated to be the most influential term, underscoring the necessity for hybrid density functional approximations in accurately simulating core-level spectra. These trends hold irrespective of the chosen basis set and exchange-correlation functional, providing valuable insights for the development of approximate methods for incorporating two-electron interactions within the realm of core-level spectroscopies.

First-Principles Anharmonic Infrared and Raman Vibrational Spectra of Materials: Fermi Resonance in Dry Ice.

We introduce a computational tool for the quantum-mechanical simulation of anharmonic infrared and Raman vibrational spectra of materials. The approach, implemented in the CRYSTAL software, stems from Taylor's expansion of the potential energy surface (PES) on the basis of normal modes up to cubic and quartic terms. The PES can be sampled with four different numerical schemes at the level of density functional theory (DFT), with local, generalized-gradient, and hybrid density functional approximations. Anharmonic states are obtained by solving Shrödinger's nuclear equation with either the vibrational self-consistent field (VSCF) or vibrational configuration interaction (VCI) methods. Nuclear quantum effects (NQEs) are thus fully accounted for. Infrared intensities are computed numerically through a Berry phase approach or analytically through a coupled-perturbed (CP) approach. Raman intensities are computed analytically via the CP approach. A variety of anharmonic features of vibrational spectra of materials can be simulated, including band shifts, combination bands, overtones, resonances (first-order Fermi, second-order Darling-Dennison), and hot bands. We showcase the effectiveness of the approach on the description of a first-order Fermi resonance (FR) in CO2 dry ice: a challenging test-case given that the FR occurs in the Raman spectrum, requires NQEs, and involves two- and three-mode couplings. Fundamental mechanistic differences with respect to the well-known FR in molecular CO2 are addressed. This application represents the first quantum-mechanical, periodic description of FR in dry ice.

Improving XYG3-type doubly hybrid approximation using self-interaction corrected SCAN density and orbitals via the PZ-SIC framework: The xDH@SCAN(SIC) approach

XYG3-type doubly hybrid (xDH) approximations have gained widespread recognition for their accuracy in describing a diverse range of chemical and physical interactions. However, a recent study [Song et al., J. Phys. Chem. Lett. 12, 800–807 (2021)] has highlighted the limitation of xDH methods in calculating the dissociation of NaCl molecules. This issue has been related to the density and orbitals used for evaluating the energy in xDH methods, which are obtained from lower-rung hybrid density functional approximations (DFAs) and display substantial density errors in the dissociation limit. In this work, we systematically investigate the influence of density on several challenging datasets and find that xDH methods are less sensitive to density errors compared to semi-local and hybrid DFAs. Furthermore, we demonstrate that the self-interaction corrected SCAN density approach offers superior accuracy compared to the self-consistent SCAN density and Hartree–Fock density approaches, as evidenced by performing charge analysis on the dissociation of heterodimers, such as NaCl and LiF. Building on these insights, we propose a five-parameter xDH method using the SCAN density and orbitals corrected by the PZ-SIC scheme. This new xDH@SCAN(SIC) method provides a balanced and accurate description across a wide range of challenging systems.

The role of electron correlations in the electronic structure of putative Chern magnet TbMn6Sn6

A member of the RMn6Sn6 rare-earth family materials, TbMn6Sn6, recently showed experimental signatures of the realization of a quantum-limit Chern magnet. In this work, we use quantum Monte Carlo (QMC) and density functional theory with Hubbard U (DFT + U) calculations to examine the electronic structure of TbMn6Sn6. To do so, we optimize accurate, correlation-consistent pseudopotentials for Tb and Sn using coupled-cluster and configuration–interaction (CI) methods. We find that DFT + U and single-reference QMC calculations suffer from the same overestimation of the magnetic moments as meta-GGA and hybrid density functional approximations. Our findings point to the need for improved orbitals/wavefunctions for this class of materials, such as natural orbitals from CI, or for the inclusion of multi-reference effects that capture the static correlations for an accurate prediction of magnetic properties. DFT + U with Mn magnetic moments adjusted to the experiment predict the Dirac crossing in bulk to be close to the Fermi level, within ~120 meV, in agreement with the experiments. Our non-stoichiometric slab calculations show that the Dirac crossing approaches even closer to the Fermi level, suggesting the possible realization of Chern magnetism in this limit.

Enabling Molecular-Level Computational Description of Redox and Proton-Coupled Electron Transfer Reactions of Samarium Diiodide.

Samarium diiodide (SmI2, Kagan's reagent) is a one-electron reductant with applications ranging from organic synthesis to nitrogen fixation. Highly inaccurate relative energies of redox and proton-coupled electron transfer (PCET) reactions of Kagan's reagent are predicted by pure and hybrid density functional approximations (DFAs) when only scalar relativistic effects are accounted for. Calculations including spin-orbit coupling (SOC) show that the SOC-induced differential stabilization of the Sm(III) versus the Sm(II) ground state is little affected by ligands and solvent, and a standard SOC correction derived from atomic energy levels is thus included in the reported relative energies. With this correction, selected meta-GGA and hybrid meta-GGA functionals predict Sm(III)/Sm(II) reduction free energies to within 5 kcal/mol of the experiment. Considerable discrepancies remain, however, in particular for the PCET-relevant O-H bond dissociation free energies, for which no regular DFA is within 10 kcal/mol of the experiment or CCSD(T). The main cause behind these discrepancies is the delocalization error, which leads to excess ligand-to-metal electron donation and destabilizes Sm(III) versus Sm(II). Fortunately, static correlation is unimportant for the present systems, and the error may be reduced by including information from virtual orbitals via perturbation theory. Contemporary, parametrized double-hybrid methods offer promise as companions to experimental campaigns in the further development of the chemistry of Kagan's reagent.

CRYSTAL23: A Program for Computational Solid State Physics and Chemistry.

The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.

Quasi-Relativistic Calculation of EPR g Tensors with Derivatives of the Decoupling Transformation, Gauge-Including Atomic Orbitals, and Magnetic Balance

We present an exact two-component (X2C) ansatz for the EPR g tensor using gauge-including atomic orbitals (GIAOs) and a magnetically balanced basis set expansion. In contrast to previous X2C and four-component relativistic ansätze for the g tensor, this implementation results in a gauge-origin-invariant formalism. Furthermore, the derivatives of the relativistic decoupling matrix are incorporated to form the complete analytical derivative of the X2C Hamiltonian. To reduce the associated computational costs, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The quasi-relativistic X2C and DLU-X2C Hamiltonians accurately reproduce the results of the parent four-component relativistic theory when accounting for two-electron picture-change effects with the modified screened nuclear spin-orbit approximation in the respective one-electron integrals and integral derivatives. According to our benchmark studies, the uncontracted Dyall and segmented-contracted Karlsruhe x2c-type basis sets perform well when compared to large even-tempered basis sets. Moreover, (range-separated) hybrid density functional approximations such as LC-ωPBE and ωB97X-D are needed to match the experimental findings. The impact of the GIAOs depends on the distribution of the spin density, and their use may change the Δg shifts by 10-50% as shown for [(C5Me5)2Y(μ-S)2Mo(μ-S)2Y(C5Me5)2]-. Routine calculations of large molecules are possible with widely available and comparably low-cost hardware as demonstrated for [Pt(C6Cl5)4]- with 3003 basis functions and three spin-(1/2) La(II) and Lu(II) compounds, for which we observe good agreement with the experimental findings.

Metal-organic frameworks properties from hybrid density functional approximations.

The chemical versatility and modular nature of Metal-Organic Frameworks (MOFs) make them unique hybrid inorganic-organic materials for several important applications. From a computational point of view, ab initio modeling of MOFs is a challenging and demanding task, in particular, when the system reaches the size of gigantic MOFs as MIL-100 and MIL-101 (where MIL stands for Materials Institute Lavoisier) with several thousand atoms in the unit cell. Here, we show how such complex systems can be successfully tackled by a recently proposed class of composite electronic structure methods revised for solid-state calculations. These methods rely on HF/density functional theory hybrid functionals (i.e., PBEsol0 and HSEsol) combined with a double-zeta quality basis set. They are augmented with semi-classical corrections to take into account dispersive interactions (D3 scheme) and the basis set superposition error (gCP). The resulting methodologies, dubbed "sol-3c," are cost-effective yet reach the hybrid functional accuracy. Here, sol-3c methods are effectively applied to predict the structural, vibrational, electronic, and adsorption properties of some of the most common MOFs. Calculations are feasible even on very large MOFs containing more than 2500 atoms in the unit cell as MIL-100 and MIL-101 with reasonable computing resources. We propose to use our composite methods for the routine in silico screening of MOFs targeting properties beyond plain structural features.

Assessing Density Functional Theory for Chemically Relevant Open-Shell Transition Metal Reactions.

Due to the principle lack of systematic improvement possibilities of density functional theory, careful assessment of the performance of density functional approximations (DFAs) on well-designed benchmark sets, for example, for reaction energies and barrier heights, is crucial. While main-group chemistry is well covered by several available sets, benchmark data for transition metal chemistry is sparse. This is especially the case for larger, chemically relevant molecules. Addressing this issue, we recently introduced the MOR41 benchmark which covers chemically relevant reactions of closed-shell complexes. In this work, we extend these efforts to single-reference open-shell systems and introduce the "reactions of open-shell single-reference transition metal complexes" (ROST61) benchmark set. ROST61 includes accurate coupled-cluster reference values for 61 reaction energies with a mean reaction energy of -42.8 kcal mol-1. Complexes with 13-93 atoms covering 20 d-block elements are included, but due to the restriction to single-reference open-shell systems, important elements such as iron or platinum could not be taken into account, or only to a small extent. We assess the performance of 31 DFAs in combination with three London dispersion (LD) correction schemes. Further, DFT-based composite methods, MP2, and a few semiempirical quantum chemical methods are evaluated. Consistent with the results for the MOR41 closed-shell benchmark, we find that the ordering of DFAs according to Jacob's ladder is preserved and that adding an LD correction is crucial, clearly improving almost all tested methods. The recently introduced r2SCAN-3c composite method stands out with a remarkable mean absolute deviation (MAD) of only 2.9 kcal mol-1, which is surpassed only by hybrid DFAs with low amounts of Fock exchange (e.g., 2.3 kcal mol-1 for TPSS0-D4/def2-QZVPP) and double-hybrid (DH) DFAs but at a significantly higher computational cost. The lowest MAD of only 1.6 kcal mol-1 is obtained with the DH DFA PWPB95-D4 in the def2-QZVPP basis set approaching the estimated accuracy of the reference method. Overall, the ROST61 set adds important reference data to a sparsely sampled but practically relevant area of chemistry. At this point, it provides valuable orientation for the application and development of new DFAs and electronic structure methods in general.

Assessing the Accuracy of Local Hybrid Density Functional Approximations for Molecular Response Properties.

A comprehensive overview of the performance of local hybrid functionals for molecular properties like excited states, ionization potentials within the GW framework, polarizabilities, magnetizabilities, NMR chemical shifts, and NMR spin-spin coupling constants is presented. We apply the generalization of the kinetic energy, τ, with the paramagnetic current density to all magnetic properties and the excitation energies from time-dependent density functional theory. This restores gauge invariance for these properties. Different ansätze for local mixing functions such as the iso-orbital indicator, the correlation length, the Görling-Levy second-order limit, and the spin polarization are compared. For the latter, we propose a modified version of the corresponding hyper-generalized gradient approximation functional of Perdew, Staroverov, Tao, and Scuseria (PSTS) [Phys. Rev. A 2008, 78, 052513] to allow for a numerically stable evaluation of the exchange-correlation kernel and hyperkernel. The PSTS functional leads to a very consistent improvement compared to the related TPSSh functional. It is further shown that the "best" choice of the local mixing function depends on the studied property and molecular class. While functionals based on the iso-orbital indicator lead to rather accurate excitation energies and ionization energies, the results are less impressive for NMR properties, for which a considerable dependence on the considered molecular test set and nuclei is observed. Johnson's local hybrid functional based on the correlation length yields remarkable results for NMR shifts of compounds featuring heavy elements and also for the excitation energies of organic compounds.

Doubly Hybrid Functionals Close to Chemical Accuracy for Both Finite and Extended Systems: Implementation and Test of XYG3 and XYGJ-OS.

While being widely used to understand the chemical reactions in heterogeneous catalysis or other multidisciplinary systems, a great challenge that semilocal and hybrid density functional approximations (DFAs) are facing is to deliver a uniformly accurate description for both finite and extended systems. Herein, we perform reliable and well-converged periodic calculations of two doubly hybrid approximations (DHAs), XYG3 and XYGJ-OS, and demonstrate that the good accuracy of DHAs achieved for molecules is transferable to the semiconductors and insulators. Such an accuracy is not only for energetic properties but also for the first- and second-order response properties, which is general for different kinds of chemical environments, including simple cubic bulks, perovskite-type transition metal oxides like TiO2, and heterogeneous systems like CO adsorption on the NaCl(100) surface. The present finding has strengthened the predictive power of DFT, which not only will inspire the future development of the top-rung DFAs but also will boost their applications in multidisciplinary studies with high accuracy and efficiency.

An improved seminumerical Coulomb and exchange algorithm for properties and excited states in modern density functional theory.

A seminumerical algorithm capable of performing large-scale (time-dependent) density functional theory (TD-DFT) calculations to extract excitation energies and other ground-state and excited-state properties is outlined. The algorithm uses seminumerical integral techniques for evaluating Coulomb and exchange parts for a set of density matrices as occurring in standard TD-DFT or similar methods for the evaluation of vibrational frequencies. A suitable optimized de-aliasing procedure is introduced. The latter does not depend on further auxiliary quantities and retains the symmetry of a given density matrix. The algorithm is self-contained and applicable to any orbital basis set available without the need for further auxiliary basis sets or optimized de-aliasing grids. Relativistic two-component excited-state TD-DFT calculations are reported for the first time using the developed seminumerical algorithm for standard and local hybrid density functional approximations. Errors are compared with the widely used "resolution of the identity" (RI) approximations for Coulomb (RI-J) and exchange integrals (RI-K). The fully seminumerical algorithm does not exhibit an enlarged error for standard DFT functionals compared to the RI approximation. For the more involved local hybrid functionals and within strong external fields, accuracy is even considerably improved.

Extended Koopmans' theorem in the adiabatic connection formalism: Applied to doubly hybrid density functionals.

A rigorous framework that combines the extended Koopmans' theorem (EKT) with the adiabatic connection (AC) formalism of density functional theory is developed here, namely, EKT-AC, to calculate the vertical ionization potentials (IPs) of molecular systems. When applied to the doubly hybrid density functional approximations (DH-DFAs), the EKT-DH approach is established for the B2PLYP-type DHs with one-parameter and two-parameters, as well as the XYG3-type DHs. Based on EKT-DH, an approximation of the KT type is introduced, leading to the KT-DH approach. The IP-condition that the calculated vertical IPs with EKT-DH or KT-DH are to reproduce the experimental IPs closely is applied to investigate the commonly used DH-DFAs for such a purpose and is utilized as a principle for DH-DFA developments. Considering the systematic improvements, as well as its numeric stability, we recommend the KT-B2GPPLYP approach as a pragmatic way for vertical IP calculations.

Accuracy of Localized Resolution of the Identity in Periodic Hybrid Functional Calculations with Numerical Atomic Orbitals.

We present an implementation of hybrid density functional approximations for periodic systems within a pseudopotential-based, numerical atomic orbital (NAO) framework. The two-electron Coulomb repulsion integrals (ERIs) are evaluated using the localized resolution-of-the-identity (LRI) approximation. The accuracy of the LRI approximation is benchmarked unambiguously against independent reference results obtained via a computational scheme whereby the ERIs are accurately evaluated by expanding the products of NAOs in terms of plane waves. An alternative strategy for constructing auxiliary basis sets is proposed, and its accuracy is assessed and compared to the previously used procedure. Finally, the reliability of our algorithm and implementation is benchmarked against other established implementations within different numerical frameworks in terms of the calculated band gap values of a set of semiconductors and insulators.

Accuracy of density functional theory methods for the calculation of magnetic exchange couplings in binuclear iron(III) complexes

We assess the performance of popular hybrid density functional approximations for the prediction of magnetic exchange couplings (J) in FeIII complexes. To this end, eleven dinuclear oxo-bridged iron(III) complexes with accurately determined J values were chosen to cover cases of weak and strong antiferromagnetic couplings ranging from −6 cm-1 to −132 cm-1. Seven representative exchange-correlation functionals ranging from the hybrids PBEh and B3LYP to meta hybrids (M06, M062X), range separated functionals (HSE, LC-ωPBE), and LSDA with 30% of HF exchange were used for this assessment. These functionals have been suggested as good performers for this particular property in transition metal complexes and hence this assessment aims to narrow the interval of confidence for the particular case of Fe complexes. In addition, we examined the performance of an alternative method based on differential local spin rotations to calculate magnetic exchange coupling parameters in these complexes. Our results show that, although some of these functionals perform on par, the M06 and HSE functionals yield magnetic exchange couplings in closer agreement with experimental data, with mean absolute percentage errors of 8.6 and 9.2%, respectively.

Renormalized Singles Green's Function for Quasi-Particle Calculations beyond the G0 W0 Approximation.

Quasi-particle energies and band gaps in particular are critical for investigating novel materials. Commonly used density functional approximations (DFAs) systematically underestimate band gaps, and GW approximation is the established method of choice for good accuracy and reliability. However, G0 W0 has some undesired dependence on the DFA, while self-consistent GW (sc GW) is expensive and not consistent in accuracy improvement. Here a simple and efficient GRS W0 approach has been developed: a subspace diagonalization of the Hartree-Fock (HF) Hamiltonian with the DFA density matrix provides the new reference Green's function GRS that incorporates the effect of all single excitation contributions to the self-energy, thereby essentially eliminating the starting-point dependence. Calculations for molecules and large band gap solids demonstrate the significant improvement over G0 W0 and greatly reduced dependence on the initial DFA. GRS W0 approach also improve results for other bulks over G0 W0, but to a lesser extent, which could be due to the limitations in current implementation for bulks. The results demonstrate that to achieve good accuracy, it is not necessary to use hybrid DFA, which is expensive for bulks. This work should be greatly significant in making GW a more robust approach.

Computation of indirect nuclear spin-spin couplings with reduced complexity in pure and hybrid density functional approximations.

We present a (sub)linear-scaling algorithm to determine indirect nuclear spin-spin coupling constants at the Hartree-Fock and Kohn-Sham density functional levels of theory. Employing efficient integral algorithms and sparse algebra routines, an overall (sub)linear scaling behavior can be obtained for systems with a non-vanishing HOMO-LUMO gap. Calculations on systems with over 1000 atoms and 20 000 basis functions illustrate the performance and accuracy of our reference implementation. Specifically, we demonstrate that linear algebra dominates the runtime of conventional algorithms for 10 000 basis functions and above. Attainable speedups of our method exceed 6 × in total runtime and 10 × in the linear algebra steps for the tested systems. Furthermore, a convergence study of spin-spin couplings of an aminopyrazole peptide upon inclusion of the water environment is presented: using the new method it is shown that large solvent spheres are necessary to converge spin-spin coupling values.

On the accuracy of commonly used density functional approximations in determining the elastic constants of insulators and semiconductors.

We have performed density functional calculations using a range of local and semi-local as well as hybrid density functional approximations of the structure and elastic constants of 18 semiconductors and insulators. We find that most of the approximations have a very small error in the lattice constants, of the order of 1%, while the errors in the elastic constants and bulk modulus are much larger, at about 10% or better. When comparing experimental and theoretical lattice constants and bulk modulus we have included zero-point phonon effects. These effects make the experimental reference lattice constants 0.019 Å smaller on average while making the bulk modulus 4.3 GPa stiffer on average. According to our study, the overall best performing density functional approximations for determining the structure and elastic properties are the PBEsol functional, the two hybrid density functionals PBE0 and HSE (Heyd, Scuseria, and Ernzerhof), as well as the AM05 functional.

Improved Electronic Excitation Energies from Shape-Corrected Semilocal Kohn-Sham Potentials

We propose a general method for obtaining accurate valence and Rydberg excitation energies from standard density-functional approximations in adiabatic linear-response time-dependent density-functional theory. The method consists in modeling the sum of Hartree (Coulomb) and exchange-correlation potentials, v(HXC)(r), by the Hartree-exchange-correlation potential of the corresponding partially ionized system in which a fraction of electron charge (δ = 0.15 to 0.30, depending on the functional) is removed from the highest occupied Kohn-Sham orbital level. The model potential is less repulsive and closer to exact in valence and near asymptotic regions, so it yields more accurate Kohn-Sham orbitals and orbital eigenvalues. By applying this scheme to conventional local, semilocal, and hybrid density-functional approximations, we improve their accuracy for Rydberg excitations by almost an order of magnitude without sacrificing the already good performance for valence transitions.