Strong-correlation density functionals made simple.

Recent work on incorporating strong-correlation (sc) corrections into the scLH22t local hybrid functional [A. Wodynski and M. Kaupp, J. Chem. Theory Comput. 18, 6111-6123 (2022)] used a hybrid procedure, applying a strong-correlation factor derived from the reverse Becke-Roussel machinery of the KP16/B13 and B13 functionals to the nonlocal correlation term of a local hybrid functional. Here, we show that adiabatic-connection factors for strong-correlation-corrected local hybrids (scLHs) can be constructed in a simplified way based on a comparison of semi-local and exact exchange-energy densities only, without recourse to exchange-hole normalization. The simplified procedure is based on a comparative analysis of Becke's B05 real-space treatment of nondynamical correlation and that in LHs, and it allows us to use, in principle, any semi-local exchange-energy density in the variable used to construct local adiabatic connections. The derivation of competitive scLHs is demonstrated based on either a modified Becke-Roussel or a simpler Perdew-Burke-Ernzerhof (PBE) energy density, leading to the scLH23t-mBR and scLH23t-tPBE functionals, which both exhibit low fractional spin errors while retaining good performance for weakly correlated situations. We also report preliminary attempts toward more detailed modeling of the local adiabatic connection, allowing a reduction of unphysical local maxima in spin-restricted bond-dissociation energy curves (scLH23t-mBR-P form). The simplified derivations of sc-factors reported here provide a basis for future constructions and straightforward implementation of exchange-correlation functionals that escape the zero-sum game between low self-interaction and static-correlation errors.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science.

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Local hybrid functionals augmented by a strong-correlation model.

The strong-correlation factor of the recent KP16/B13 exchange-correlation functional has been adapted and applied to the framework of local hybrid (LH) functionals. The expression identifiable as nondynamical (NDC) and dynamical (DC) correlations in LHs is modified by inserting a position-dependent KP16/B13-style strong-correlation factor qAC(r) based on a local version of the adiabatic connection. Different ways of deriving this factor are evaluated for a simple one-parameter LH based on the local density approximation. While the direct derivation from the LH NDC term fails due to known deficiencies, hybrid approaches, where the factor is determined from the B13 NDC term as in KP16/B13 itself, provide remarkable improvements. In particular, a modified B13 NDC expression using Patra's exchange-hole curvature showed promising results. When applied to the simple LH as a first attempt, it reduces atomic fractional-spin errors and deficiencies of spin-restricted bond dissociation curves to a similar extent as the KP16/B13 functional itself while maintaining the good accuracy of the underlying LH for atomization energies and reaction barriers in weakly correlated situations. The performance of different NDC expressions in deriving strong-correlation corrections is analyzed, and areas for further improvements of strong-correlation corrected LHs and related approaches are identified. All the approaches evaluated in this work have been implemented self-consistently into a developers' version of the Turbomole program.

Hyperfine-Coupling Tensors from Projected Hartree-Fock Theory.

We assess the calculation of hyperfine coupling (HFC) tensors by different variants of Projected Hartree-Fock (PHF) theory. For a set of small main-group S = 1/2 radicals (BO, CO+, CN, AlO, vinyl, methyl, ethynyl), spin-symmetry as well as complex-conjugation and point-group symmetry are first broken in a reference determinant, and then variationally restored, in the frame of the modern formulation of PHF theory. Historically, PHF theory was basically restricted to the restoration of spin symmetry from an unrestricted HF determinant (conserving Sz symmetry). This afforded unsatisfactory HFCs. We obtain far better results for isotropic (and anisotropic) HFCs when the variational energy is further lowered by working with generalized determinants that completely break spin symmetry, and when additional symmetries are used. Specifically, complex-conjugation projection recovers a substantial fraction of the dynamical correlation energy in small molecules, and the detailed equations for combined complex-conjugation, spin- and point-group projection in the density-matrix/diagonalization formulation of PHF theory are here reported for the first time. The compact representation of the PHF wave function allows for a straightforward evaluation of the spin-density matrix and of HFC tensors with little effort. The promising performance of PHF theory may motivate the application of post-PHF methods to the calculation of HFC tensors.

A Local Hybrid Functional with Wide Applicability and Good Balance between (De)Localization and Left-Right Correlation.

A new local hybrid functional, LH20t, with a position-dependent exact-exchange admixture governed by a simple local mixing function (g(r) = b·τW(r)/τ(r)), combined with gradient-corrected (PBE) exchange and meta-GGA (B95) correlation, as well as a second-order GGA-based pig2 calibration function to address the ambiguity of exchange-energy densities, has been constructed. The adjustable parameters of LH20t have been optimized in a multistep procedure based on thermochemical kinetics data and measures of spurious nondynamical correlation. LH20t has subsequently been evaluated for the full GMTKN55 main-group energetics test suite, with and without an added DFT-D4 dispersion correction. Performance of the new functional in the GMTKN55 tests is excellent, better than any global hybrid so far, approaching the best results for any rung-4 functional, without any noticeable artifacts due to the gauge ambiguity. The robust performance across the board is combined with enhanced exact-exchange admixtures of >70% near the nuclei and asymptotically (but low admixture in bonds). This helps to provide excellent performance for a wide variety of excitation classes (core, valence singlet and triplet, Rydberg, short-range intervalence charge-transfer) in TDDFT evaluations. Notably, LH20t is the first functional that provides simultaneously the correct description for the most extreme localized and delocalized cases of the MVO-10 test set of gas-phase mixed-valence systems. This outstanding performance for mixed-valence systems, which signals a very fine balance between reduced delocalization errors and a reasonable description of left-right correlation, is corroborated by tests on ground- and excited-state properties for organic and organometallic mixed-valence systems in solution.

Hubbard Trimer with Spin-Orbit Coupling: Hartree-Fock Solutions, (Non)Collinearity, and Anisotropic Spin Hamiltonian.

We present unrestricted and generalized Hartree-Fock solutions (UHF and GHF, respectively) for the single-band Hubbard model of an equilateral triangle. Spin-orbit coupling (SOC) is treated self-consistently, and HF stability and properties of different spin structures are studied in detail. The GHF solution switches from noncollinear to collinear when crossing a high-symmetry point in parameter space (spanned by the amplitudes of spin-conserving and spin-dependent hopping, i.e., kinetic energy and SOC, respectively). The collinear GHF solution represents a simple example to disprove the notion that a collinear vector spin density in a Slater determinant necessarily entails a defined spin projection. Spin Hamiltonian parameters for the anisotropic interaction between three spin-1/2 centers are extracted from HF energies and subsequently compared to exact results from effective Hamiltonian theory. This provides an unambiguous benchmark for interpreting broken-symmetry mean-field solutions in terms of spin configurations and puts this semiclassical approach (frequently applied in broken-symmetry density functional theory) on a firmer basis.

Exact Mapping from Many-Spin Hamiltonians to Giant-Spin Hamiltonians.

Thermodynamic and spectroscopic data of exchange-coupled molecular spin clusters (e.g. single-molecule magnets) are routinely interpreted in terms of two different models: the many-spin Hamiltonian (MSH) explicitly considers couplings between individual spin centers, while the giant-spin Hamiltonian (GSH) treats the system as a single collective spin. When isotropic exchange coupling is weak, the physical compatibility between both spin Hamiltonian models becomes a serious concern, due to mixing of spin multiplets by local zero-field splitting (ZFS) interactions ('S-mixing'). Until now, this effect, which makes the mapping MSH→GSH ('spin projection') non-trivial, had only been treated perturbationally (up to third order), with obvious limitations. Here, based on exact diagonalization of the MSH, canonical effective Hamiltonian theory is applied to construct a GSH that exactly matches the energies of the relevant (2S+1) states comprising an effective spin multiplet. For comparison, a recently developed strategy for the unique derivation of effective ('pseudospin') Hamiltonians, now routinely employed in ab initio calculations of mononuclear systems, is adapted to the problem of spin projection. Expansion of the zero-field Hamiltonian and the magnetic moment in terms of irreducible tensor operators (or Stevens operators) yields terms of all ranks k (up to k=2S) in the effective spin. Calculations employing published MSH parameters illustrate exact spin projection for the well-investigated [Ni(hmp)(dmb)Cl]4 ('Ni4 ') single-molecule magnet, which displays weak isotropic exchange (dmb=3,3-dimethyl-1-butanol, hmp- is the anion of 2-hydroxymethylpyridine). The performance of the resulting GSH in finite field is assessed in terms of EPR resonances and diabolical points. The large tunnel splitting in the M=± 4 ground doublet of the S=4 multiplet, responsible for fast tunneling in Ni4 , is attributed to a Stevens operator with eightfold rotational symmetry, marking the first quantification of a k=8 term in a spin cluster. The unique and exact mapping MSH→GSH should be of general importance for weakly-coupled systems; it represents a mandatory ultimate step for comparing theoretical predictions (e.g. from quantum-chemical calculations) to ZFS, hyperfine or g-tensors from spectral fittings.

Understanding Thermodynamic and Spectroscopic Properties of Tetragonal Mn12 Single-Molecule Magnets from Combined Density Functional Theory/Spin-Hamiltonian Calculations.

We apply broken-symmetry density functional theory to determine isotropic exchange-coupling constants and local zero-field splitting (ZFS) tensors for the tetragonal Mn12(t)BuAc single-molecule magnet. The obtained parametrization of the many-spin Hamiltonian (MSH), taking into account all 12 spin centers, is assessed by comparing theoretical predictions for thermodynamic and spectroscopic properties with available experimental data. The magnetic susceptibility (calculated by the finite-temperature Lanczos method) is well approximated, and the intermultiplet excitation spectrum from inelastic neutron scattering (INS) experiments is correctly reproduced. In these respects, the present parametrization of the 12-spin model represents a significant improvement over previous theoretical estimates of exchange-coupling constants in Mn12, and additionally offers a refined interpretation of INS spectra. Treating anisotropic interactions at the third order of perturbation theory, the MSH is mapped onto the giant-spin Hamiltonian describing the S = 10 ground multiplet. Although the agreement with high-field EPR experiments is not perfect, the results clearly point in the right direction and for the first time rationalize the angular dependence of the transverse-field spectra from a fully microscopic viewpoint. Importantly, transverse anisotropy of the effective S = 10 manifold is explicitly shown to arise largely from the ZFS-induced mixing of exchange multiplets. This effect is given a thorough analysis in the approximate D2d spin-permutational symmetry group of the exchange Hamiltonian.

Construction of Giant-Spin Hamiltonians from Many-Spin Hamiltonians by Third-Order Perturbation Theory and Application to an Fe3 Cr Single-Molecule Magnet.

A general giant-spin Hamiltonian (GSH) describing an effective spin multiplet of an exchange-coupled metal cluster with dominant Heisenberg interactions was derived from a many-spin Hamiltonian (MSH) by treating anisotropic interactions at the third order of perturbation theory. Going beyond the existing second-order perturbation treatment allows irreducible tensor operators of rank six (or corresponding Stevens operator equivalents) in the GSH to be obtained. Such terms were found to be of crucial importance for the fitting of high-field EPR spectra of a number of single-molecule magnets (SMMs). Also, recent magnetization measurements on trigonal and tetragonal SMMs have found the inclusion of such high-rank axial and transverse terms to be necessary to account for experimental data in terms of giant-spin models. While mixing of spin multiplets by local zero-field splitting interactions was identified as the major origin of these contributions to the GSH, a direct and efficient microscopic explanation had been lacking. The third-order approach developed in this work is used to illustrate the mapping of an MSH onto a GSH for an S=6 trigonal Fe3 Cr complex that was recently investigated by high-field EPR spectroscopy. Comparisons between MSH and GSH consider the simulation of EPR data with both Hamiltonians, as well as locations of diabolical points (conical intersections) in magnetic-field space. The results question the ability of present high-field EPR techniques to determine high-rank zero-field splitting terms uniquely, and lead to a revision of the experimental GSH parameters of the Fe3 Cr SMM. Indeed, a bidirectional mapping between MSH and GSH effectively constrains the number of free parameters in the GSH. This notion may in the future facilitate spectral fitting for highly symmetric SMMs.

New approaches for the calibration of exchange-energy densities in local hybrid functionals.

The ambiguity of exchange-energy densities is a fundamental challenge for the development of local hybrid functionals, or of other functionals based on a local mixing of exchange-energy densities. In this work, a systematic construction of semi-local calibration functions (CFs) for adjusting the exchange-energy densities in local hybrid functionals is provided, which directly links a given CF to an underlying semi-local exchange functional, as well as to the second-order gradient expansion of the exchange hole. Using successive steps of integration by parts allows the derivation of correction terms of increasing order, resulting in more and more complicated but also more flexible CFs. We derive explicit first- and second-order CFs (pig1 and pig2) based on B88 generalized-gradient approximation (GGA) exchange, and a first-order CF (tpig1) based on τ-dependent B98 meta-GGA exchange. We combine these CFs with different long-range damping functions and evaluate them for calibration of LDA, B88 GGA, and TPSS meta-GGA exchange-energy densities. Based on a minimization of unphysical nondynamical correlation contributions in three noble-gas dimer potential-energy curves, free parameters in the CFs are optimized, and performance of various approaches in the calibration of different exchange-energy densities is compared. Most notably, the second-order pig2 CF provides the largest flexibility with respect to the diffuseness of the damping function. This suggests that higher-order CFs based on the present integration-by-parts scheme may be particularly suitable for the flexible construction of local hybrid functionals.

Validation of local hybrid functionals for TDDFT calculations of electronic excitation energies.

The first systematic evaluation of local hybrid functionals for the calculation of electronic excitation energies within linear-response time-dependent density functional theory (TDDFT) is reported. Using our recent efficient semi-numerical TDDFT implementation [T. M. Maier et al., J. Chem. Theory Comput. 11, 4226 (2015)], four simple, thermochemically optimized one-parameter local hybrid functionals based on local spin-density exchange are evaluated against a database of singlet and triplet valence excitations of organic molecules, and against a mixed database including also Rydberg, intramolecular charge-transfer (CT) and core excitations. The four local hybrids exhibit comparable performance to standard global or range-separated hybrid functionals for common singlet valence excitations, but several local hybrids outperform all other functionals tested for the triplet excitations of the first test set, as well as for relative energies of excited states. Evaluation for the combined second test set shows that local hybrids can also provide excellent Rydberg and core excitations, in the latter case rivaling specialized functionals optimized specifically for such excitations. This good performance of local hybrids for different excitation types could be traced to relatively large exact-exchange (EXX) admixtures in a spatial region intermediate between valence and asymptotics, as well as close to the nucleus, and lower EXX admixtures in the valence region. In contrast, the tested local hybrids cannot compete with the best range-separated hybrids for intra- and intermolecular CT excitation energies. Possible directions for improvement in the latter category are discussed. As the used efficient TDDFT implementation requires essentially the same computational effort for global and local hybrids, applications of local hybrid functionals to excited-state problems appear promising in a wide range of fields. Influences of current-density dependence of local kinetic-energy dependent local hybrids, differences between spin-resolved and "common" local mixing functions in local hybrids, and the effects of the Tamm-Dancoff approximation on the excitation energies are also discussed.

A Relativistic Quantum-Chemical Analysis of the trans Influence on (1)H NMR Hydride Shifts in Square-Planar Platinum(II) Complexes.

Empirical correlations between characteristic (1)H NMR shifts in Pt(II) hydrides with trans ligand influence series, Pt-H distances, and (195)Pt shifts are analyzed at various levels of including relativistic effects into density-functional calculations. A close examination of the trans ligand effects on hydride NMR shifts is shown to be dominated by spin-orbit shielding σ(SO). A rather complete understanding of the trends has been obtained by detailed molecular orbital (MO)-by-MO and localized MO analyses of the paramagnetic and spin-orbit (SO) contributions to the chemical shifts, noting that it is the perpendicular shift-tensor components that determine the trend of the (1)H hydride shifts. In contrast to previous assumptions, the change of the Pt-H distance in given complexes does not allow correlations between hydride shifts and metal-hydrogen bond length to be understood. Instead, variations in the polarization of metal 5d orbitals by the trans ligand affects the SO (and partly paramagnetic) shift contributions, as well as the Pt-H distances and the covalency of the metal-hydrogen bond (quantified, e.g., by natural atomic charges and delocalization indices from quantum theory atoms-in-molecules), resulting in a reasonable correlation of these structural/electronic quantities with hydride σ(SO) shieldings. Our analysis also shows that specific σ(p)- and σ(SO)-active MOs are not equally important across the entire series. This explains some outliers in the correlation for limited ranges of trans-influence ligands. Additionally, SO effects from heavy-halide ligands may further complicate trends, indicating some limitations of the simple one-parameter correlations. Strikingly, σ-donating/π-accepting ligands with a very strong trans influence are shown to invert the sign of the usually shielding σ(SO) contribution to the (1)H shifts, by a substantial reduction of the metal 5d orbital involvement in Pt-H bonding, and by involvement of metal 6p-type orbitals in the magnetic couplings, in violation of the Buckingham-Stephens "off-center ring-current" picture.

Towards improved local hybrid functionals by calibration of exchange-energy densities.

A new approach for the calibration of (semi-)local and exact exchange-energy densities in the context of local hybrid functionals is reported. The calibration functions are derived from only the electron density and its spatial derivatives, avoiding spatial derivatives of the exact-exchange energy density or other computationally unfavorable contributions. The calibration functions fulfill the seven more important out of nine known exact constraints. It is shown that calibration improves substantially the definition of a non-dynamical correlation energy term for generalized gradient approximation (GGA)-based local hybrids. Moreover, gauge artifacts in the potential-energy curves of noble-gas dimers may be corrected by calibration. The developed calibration functions are then evaluated for a large range of energy-related properties (atomization energies, reaction barriers, ionization potentials, electron affinities, and total atomic energies) of three sets of local hybrids, using a simple one-parameter local-mixing. The functionals are based on (a) local spin-density approximation (LSDA) or (b) Perdew-Burke-Ernzerhof (PBE) exchange and correlation, and on (c) Becke-88 (B88) exchange and Lee-Yang-Parr (LYP) correlation. While the uncalibrated GGA-based functionals usually provide very poor thermochemical data, calibration allows a dramatic improvement, accompanied by only a small deterioration of reaction barriers. In particular, an optimized BLYP-based local-hybrid functional has been found that is a substantial improvement over the underlying global hybrids, as well as over previously reported LSDA-based local hybrids. It is expected that the present calibration approach will pave the way towards new generations of more accurate hyper-GGA functionals based on a local mixing of exchange-energy densities.

Importance of the correlation contribution for local hybrid functionals: range separation and self-interaction corrections.

Local hybrid functionals with their position-dependent exact-exchange admixture are a conceptually simple and promising extension of the concept of a hybrid functional. Local hybrids based on a simple mixing of the local spin density approximation (LSDA) with exact exchange have been shown to be successful for thermochemistry, reaction barriers, and a range of other properties. So far, the combination of this generation of local hybrids with an LSDA correlation functional has been found to give the most favorable results for atomization energies, for a range of local mixing functions (LMFs) governing the exact-exchange admixture. Here, we show that the choice of correlation functional to be used with local hybrid exchange crucially influences the parameterization also of the exchange part as well as the overall performance. A novel ansatz for the correlation part of local hybrids is suggested based on (i) range-separation of LSDA correlation into short-range (SR) and long-range (LR) parts, and (ii) partial or full elimination of the one-electron self-correlation from the SR part. It is shown that such modified correlation functionals allow overall larger exact exchange admixture in thermochemically competitive local hybrids than before. This results in improvements for reaction barriers and for other properties crucially influenced by self-interaction errors, as demonstrated by a number of examples. Based on the range-separation approach, a fresh view on the breakdown of the correlation energy into dynamical and non-dynamical parts is suggested.

Computation of hyperfine tensors for dinuclear Mn(III) Mn(IV) complexes by broken-symmetry approaches: anisotropy transfer induced by local zero-field splitting.

Based on broken-symmetry density functional calculations, the (55)Mn hyperfine tensors of a series of exchange-coupled, mixed-valence, dinuclear Mn(III) Mn(IV) complexes have been computed. We go beyond previous quantum chemical work by fully including the effects of local zero-field splitting (ZFS) interactions in the spin projection, following the first-order perturbation formalism of Sage et al. [J. Am. Chem. Soc. 1989, 111, 7239]. This allows the ZFS-induced transfer of hyperfine anisotropy from the Mn(III) site to the Mn(IV) site to be described with full consideration of the orientations of local hyperfine and ZFS tensors. After scaling to correct for systematic deficiencies in the quantum chemically computed local ZFS tensors, good agreement with experimental (55)Mn anisotropies at the Mn(IV) site is obtained. The hyperfine coupling anisotropies on the Mn(III) site depend sensitively on structural distortions for a d(4) ion. The latter are neither fully reproduced by using a DFT-optimized coordination environment nor by using experimental structures. For very small exchange-coupling constants, the perturbation treatment breaks down and a dramatic sensitivity to the scaling of the local ZFS tensors is observed. These results are discussed with respect to ongoing work to elucidate the structure of the oxygen-evolving complex of photosystem II by analysis of the EPR spectra.

Evaluation of a combination of local hybrid functionals with DFT-D3 corrections for the calculation of thermochemical and kinetic data.

Due to their position-dependent exact exchange admixture, local hybrid functionals offer a higher flexibility and thus the potential for more universal and accurate exchange correlation functionals compared to global hybrids with a constant admixture, as has been demonstrated in previous work. Yet, the local hybrid constructions used so far do not account for the inclusion of dispersion-type interactions. As a first exploratory step toward a more general approach that includes van der Waals-type interactions with local hybrids, the present work has added DFT-D3-type corrections to a number of simple local hybrid functionals. Optimization of only the s(8) and s(r,6) parameters for the S22 set provides good results for weak interaction energies but deteriorates the excellent performance of the local hybrids for G3 atomization energies and for classical reaction barriers. A combined optimization of the two DFT-D3 parameters with one of the two parameters of the spin-polarized local mixing function (LMF) of a local hybrid for a more general optimization set provides simultaneously accurate dispersion energies, improved atomization energies, and accurate reaction barriers, as well as excellent alkane protobranching ratios. For other LMFs, the improvements of such a combined optimization for the S22 energies have been less satisfactory. The most notable advantage of the dispersion-corrected local hybrids over, for example, a B3LYP-D3 approach, is in the much more accurate reaction barriers.

Density functional calculations of (55)Mn, (14)N and (13)C electron paramagnetic resonance parameters support an energetically feasible model system for the S(2) state of the oxygen-evolving complex of photosystem II.

Metal and ligand hyperfine couplings of a previously suggested, energetically feasible Mn(4)Ca model cluster (SG2009(-1)) for the S(2) state of the oxygen-evolving complex (OEC) of photosystem II (PSII) have been studied by broken-symmetry density functional methods and compared with other suggested structural and spectroscopic models. This was carried out explicitly for different spin-coupling patterns of the S=1/2 ground state of the Mn(III)(Mn(IV))(3) cluster. By applying spin-projection techniques and a scaling of the manganese hyperfine couplings, computation of the hyperfine and nuclear quadrupole coupling parameters allows a direct evaluation of the proposed models in comparison with data obtained from the simulation of EPR, ENDOR, and ESEEM spectra. The computation of (55)Mn hyperfine couplings (HFCs) for SG2009(-1) gives excellent agreement with experiment. However, at the current level of spin projection, the (55)Mn HFCs do not appear sufficiently accurate to distinguish between different structural models. Yet, of all the models studied, SG2009(-1) is the only one with the Mn(III) site at the Mn(C) center, which is coordinated by histidine (D1-His332). The computed histidine (14)N HFC anisotropy for SG2009(-1) gives much better agreement with ESEEM data than the other models, in which Mn(C) is an Mn(IV) site, thus supporting the validity of the model. The (13)C HFCs of various carboxylates have been compared with (13)C ENDOR data for PSII preparations with (13)C-labelled alanine.

Local hybrid functionals with an explicit dependence on spin polarization.

A new family of local hybrid functionals with a position-dependent mixture of local and exact exchange has been constructed. On the basis of conceptual similarities between local hybrids and explicit nondynamical correlation functionals, the relative spin polarization has been introduced as an additional variable into the local mixing function (LMF) that determines the position dependence of the exact exchange admixture. This leads to two new two-parameter LMFs, one based on the ratio of the local kinetic energy density to the von Weizsacker kinetic energy density and the other on the dimensionless density gradient. After optimization of the two free parameters for small test sets, significant improvements of atomization energies of the full G3 set are obtained (mean absolute errors of 2.96 and 3.18 kcal mol(-1), respectively), with minor deterioration for classical reaction barriers (HTBH38 and NHTBH38 test sets). The simultaneous description of two-center three-electron dimer cations remains a challenge. It is shown that the exact and local spin density exchange energy densities have closely related gauge definitions, so that gauge mismatch is only a minor problem for local hybrids based on local exchange.

On the self-consistent implementation of general occupied-orbital dependent exchange-correlation functionals with application to the B05 functional.

Occupied-orbital dependent (OOD) exchange-correlation functionals hold a particularly prominent place in current developments of density functional theory. Their self-consistent implementation is complicated by the fact that their orbital-dependent parts are not explicit but only implicit functionals of electron density, and the exchange-correlation potential may not be obtained straightforwardly by taking the functional derivative with respect to the density. A two-step procedure is required, in which initially the functional derivatives with respect to the orbitals (FDOs) are obtained, which may then be transformed into local and multiplicative potentials by techniques of the optimized-effective potential. In view of the rather large variety of OOD functionals under current study, we report here general, systematic, and transparent expressions of the FDOs of a generalized OOD functional and additionally a matrix-element version in a basis set of atomic orbitals. Explicit FDOs are for the first time derived and numerically tested for one of the currently most complex examples of an OOD functional, Becke's real-space model of nondynamical correlation (B05 functional) [J. Chem. Phys. 122, 064101 (2005)].

Coupled-Perturbed Scheme for the Calculation of Electronic g-Tensors with Local Hybrid Functionals.

A coupled-perturbed Kohn-Sham (CPKS) scheme for calculating second-order magnetic properties has been developed for the case of general occupied-orbital-dependent (OOD) exchange-correlation functionals involving the exact-exchange energy density. The origin of the coupling terms in the functional derivatives of OOD functionals with respect to the orbitals has been thoroughly analyzed, and general expressions for the resulting coupling terms have been obtained. The generalized CPKS scheme thus obtained has been implemented within the MAG-ReSpect code and tested in calculations of electronic g-tensors with local hybrid functionals. Compared to previously tested global hybrids, like B3LYP, thermochemically optimized local hybrids provide only little to moderate improvement for test sets of main-group radicals and paramagnetic transition-metal complexes. Closer analyses point to possible areas in which the fundamentally more flexible local hybrids may be improved for the property at hand.