Assessment of Double-Hybrid Density Functional Theory for Magnetic Exchange Coupling in Manganese Complexes

Molecular systems containing magnetically interacting (exchange-coupled) manganese ions are important in catalysis, biomimetic chemistry, and molecular magnetism. The reliable prediction of exchange coupling constants with quantum chemical methods is key for tracing the relationships between structure and magnetic properties in these systems. Density functional theory (DFT) in the broken-symmetry approach has been employed extensively for this purpose and hybrid functionals with moderate levels of Hartree–Fock exchange admixture have often been shown to perform adequately. Double-hybrid density functionals that introduce a second-order perturbational contribution to the Kohn–Sham energy are generally regarded as a superior approach for most molecular properties, but their performance remains unexplored for exchange-coupled manganese systems. An assessment of various double-hybrid functionals for the prediction of exchange coupling constants is presented here using a set of experimentally characterized dinuclear manganese complexes that cover a wide range of exchange coupling situations. Double-hybrid functionals perform more uniformly compared to conventional DFT methods, but they fail to deliver improved accuracy or reliability in the prediction of exchange coupling constants. Reparametrized double-hybrid density functionals (DHDFs) perform no better, and most often worse, than the original B2-PLYP double-hybrid method. All DHDFs are surpassed by the hybrid-meta-generalized gradient approximation (GGA) TPSSh functional. Possible directions for future methodological developments are discussed.


Introduction
Oligonuclear complexes with magnetically interacting manganese ions have a long and rich history in chemistry, not only as synthetic analogues of essential biological clusters, such as the oxygen-evolving complex of photosystem II, but also in their own right, because of their magnetic, spectroscopic, and catalytic properties [1][2][3][4][5][6][7][8][9][10][11][12][13][14].From the perspective of quantum chemistry, a fundamental challenge in the description of the electronic structure of these systems is to model reliably the magnetic energy levels, that is, the energy levels associated with the magnetic "interaction" between the open-shell Mn ions [15].This forms the basis for obtaining the full range of spin-dependent observables for a given system.By definition, this is a genuine multireference problem and in principle should be addressed with appropriate multireference quantum chemical methods [16][17][18][19][20][21][22].However, the challenge of directly computing the full spectrum of energy states often can or has to be reduced for practical purposes to the more modest target of extracting pairwise exchange coupling constants that parameterize the magnetic energy levels in the framework of an effective spin Hamiltonian.This is often achieved at the level of Kohn-Sham density functional theory (DFT) using the broken-symmetry approach Inorganics 2019, 7, 57; doi:10.3390/inorganics7050057www.mdpi.com/journal/inorganics(BS-DFT) [23][24][25][26][27][28][29][30][31].Various spectroscopic properties can be subsequently obtained using spin-projection techniques applied to single-determinant broken-symmetry solutions.
The BS-DFT approach has been used extensively in the study of exchange-coupled transition metal complexes, with documented achievements as well as failures.The chemical nature of a given system, the DFT functional, and the methodological details of the approach interact in physically non-transparent ways that lead to inconsistent and unpredictable behavior.This has so far precluded the establishment of a universally applicable BS-DFT approach, encouraging instead the empirical choice of an "optimal" functional for a specific type of chemical system.For manganese complexes a series of benchmark studies that compared generalized gradient approximation (GGA), meta-GGA, hybrid, and hybrid-meta-GGA functionals have suggested that hybrid functionals perform better than non-hybrid functionals for the prediction of exchange-coupling constants [32][33][34][35][36][37][38][39][40][41][42][43].Pure GGA functionals tend to yield too large antiferromagnetic exchange couplings, a result attributed to excessive delocalization of spin density.Moderate admixture of Hartree-Fock (HF) exchange appears to provide an adequate counterbalance to this behavior for manganese systems, so that the calculated exchange coupling constants agree better with experiment when the BS-DFT results are used with spin projection.Still, there is no magic number for the exact exchange admixture and the best compromise for a given type of system depends on other features of the functional.Although the popular B3LYP hybrid functional (20% HF exchange) [44,45] has been used with reasonably good results, the hybrid-meta-GGA TPSSh functional (10% HF exchange) [46] is a better and more widely validated choice, not only for exchange coupling constants but for a range of molecular properties [34,[47][48][49][50].This choice of functional is not necessarily transferable, as demonstrated for example by the fact that TPSSh is no longer at the methodological sweet spot for the calculation of exchange coupling constants even for complexes of metal ions isoelectronic with manganese [51].The adequate performance of a given functional for a restricted chemical space and the emergence of inconsistent behavior or large deviations outside this restricted space are typical of standard BS-DFT approaches.
Double-hybrid density functionals (DHDFs) [52] are regarded as "higher-rung" DFT methods compared to hybrid and meta-GGA functionals.They mix standard DFT exchange (E DFT X ) and correlation (E DFT C ) with HF exchange (E HF X ) and an additional second-order perturbation theory contribution (E PT2 C ): The E PT2 C contribution is obtained through a Møller-Plesset perturbational term (MP2) based on Kohn-Sham orbitals that are self-consistently optimized with respect to the first three terms: The archetypal example is the B2-PLYP functional [53] that uses the B88 and LYP functionals for exchange and correlation, mixing HF exchange (α X = 0.53) and perturbational correction (α C = 0.27).The percentage of HF exchange in DHDFs is usually considerably higher than that in most standard hybrid functionals.DHDFs have been shown to perform robustly and with generally superior results compared to other DFT methods in numerous benchmarks involving various types of energetics [52,[54][55][56][57][58].
Studies involving exchange-coupled systems are comparatively scarce.Schwabe and Grimme [59] evaluated double-hybrid DFT on a set of organic compounds and a few simple copper complexes [60], supporting the good performance of the method.Rajaraman and coworkers [61] reported that the double-hybrid B2-PLYP functional performed better than hybrid functionals such as B3LYP, PBE0, or TPSSh for the description of magnetic coupling in Gd(III)-radical complexes.Vogiatzis et al. [62] also reported that B2-PLYP performs well within the BS-DFT approach for dinuclear paddlewheel complexes, however they also noted that it is inferior to B3LYP in the case of a trichromium system [62].
It is, therefore, not a given that the use of a double-hybrid functional automatically can ensure better results for exchange coupling problems across transition metal chemistry.The increase in the number of DHDFs over recent years further complicates the choice of method because new functionals are typically not tested for transition metal complexes with magnetically coupled ions.Double-hybrid DFT has similarly not been evaluated so far for manganese systems of relevance to biomimetic chemistry and molecular magnetism.The oxidation states Mn(III) and Mn(IV) are most relevant in this area, that is why previous assessments of theoretical methods for magnetic and spectroscopic properties of manganese complexes have made extensive use of oligonuclear complexes with Mn ions in these oxidation states [32][33][34][35][36][37][38][39][40][41][42][43].In the present study a set of dinuclear Mn complexes are used to evaluate the performance of double-hybrid DFT for the prediction of exchange coupling constants.The manganese complexes comprising the database are crystallographically characterized compounds with experimentally determined exchange coupling constants that span the range from moderate ferromagnetic to strong antiferromagnetic coupling.A series of traditional as well as more modern DHDF formulations are included in the evaluation.The results are compared with those obtained by standard GGA, meta-GGA, hybrid-GGA and hybrid-meta-GGA functionals.

Test Set of Dinuclear Manganese Complexes
Five dinuclear complexes were selected for the present study as representative of various exchange coupling situations in manganese containing systems (see Table 1 and Figure 1).It is noted that the exchange coupling constants J in the present study adhere to the following form of the Heisenberg-Dirac-van Vleck Hamiltonian for the isotropic bilinear coupling of two spins S 1 and S 2 : Inorganics 2019, 7, x FOR PEER REVIEW 4 of 14 Recent years have seen a surge in interest in DHDFs and, as a result, a rapid increase in the number of variants and modified forms [52,58,70].In the present study we aim at a modest selection of functionals that is nevertheless representative of all widely available types of double-hybrid methods.A few non-double-hybrid methods were included to facilitate comparisons with the DHDFs.These are the GGA functional BLYP, the meta-GGA functionals TPSS [71] and SCAN [72], the hybrid functionals B3LYP (20% HF exchange) [44,45] and PBE0 (25% HF exchange) [73], and the hybrid-meta-GGA functional TPSSh (10% HF exchange) [46].The selected double-hybrid functionals include B2-PLYP [53], which is the first modern functional where a PT2 contribution replaces part of Antiferromagnetic coupling is the most common situation in complexes comprising Mn(III) and Mn(IV) ions, so only one of the systems (complex 1) represents a case of ferromagnetic coupling, with a ground state spin of S = 3 [63].The other complexes have low-spin ground states.Complex 2 is a weakly antiferromagnetically coupled Mn(III,III) dimer (S = 0) [64].Complex 3 is a mixed-valence system with asymmetric ligation [65] that facilitates spin localization.This complex has a moderate antiferromagnetic coupling resulting in a spin doublet (S = 1/2) ground state and has been the subject of a recent study [20] that evaluated the use of the density matrix renormalization group [66] in the multireference treatment of exchange coupling [18,22].Complex 4 is a classic example of a strongly coupled bis-µ-oxo Mn(IV,IV) system [67].Finally, complex 5 reaches the far limit of strong antiferromagnetic coupling [68].The tris-µ-oxo ligation in complex 5 brings the manganese ions in such close proximity that, in addition to ligand-mediated superexchange, direct metal-metal interaction contributes significantly in stabilizing the low-spin state [39].This situation is common in face-sharing d 3 -d 3 systems [51,69].a Definition of ligand abbreviations: Me 3 tacn = 1,4,7-trimethyl-1,4,7-triazacyclononane; bpy = bipyridine; bpea = N,N-bis(2-pyridylmethyl)ethylamine.

Selection of Functionals
Recent years have seen a surge in interest in DHDFs and, as a result, a rapid increase in the number of variants and modified forms [52,58,70].In the present study we aim at a modest selection of functionals that is nevertheless representative of all widely available types of double-hybrid methods.A few non-double-hybrid methods were included to facilitate comparisons with the DHDFs.These are the GGA functional BLYP, the meta-GGA functionals TPSS [71] and SCAN [72], the hybrid functionals B3LYP (20% HF exchange) [44,45] and PBE0 (25% HF exchange) [73], and the hybrid-meta-GGA functional TPSSh (10% HF exchange) [46].The selected double-hybrid functionals include B2-PLYP [53], which is the first modern functional where a PT2 contribution replaces part of the DFT correlation; it combines B88 for DFT exchange with LYP for DFT correlation, with a HF exchange coefficient α X = 0.53 and a PT2 coefficient α C = 0.27.mPW2-PLYP is similar but uses mPW exchange with α X = 0.55 and α C = 0.25.A series of reparametrized versions of the B2-PLYP functional are B2GP-PLYP (α X = 0.65 and α C = 0.36), B2K-PLYP (α X = 0.72 and α C = 0.42), and B2T-PLYP (α X = 0.60 and α C = 0.31), which were respectively suggested as more suitable for general-purpose calculations, kinetics, and thermochemistry [74].
In the spin-component-scaled MP2 approach (SCS-MP2) [75] a separate scaling factor is used for the correlation energy contribution of opposite-spin (E OS−PT2 C ) and same-spin (E SS−PT2 C ) electron pairs.This is usually an improvement over MP2 because it can take account of the fact that same-spin correlation is already present at the HF level.SCS-MP2 was also used in the context of double-hybrid density functional theory.Here we test functionals from the dispersion-corrected spin-component-scaled double-hybrid family (DSD) [76][77][78], which follow the formula: DSD-PBEP86 employs PBE exchange and P86 correlation with α X = 0.70, c C = 0.43, c O = 0.53, and c S = 0.25, while DSD-PBEB95 combines PBE exchange with B95 correlation and has coefficients α X = 0.66, c C = 0.55, c O = 0.46, and c S = 0.09.It is noted that the scaling factors for the SCS-MP2 part were reparametrized and do not follow the original SCS-MP2 definition [75].Finally, the PWPB95 functional [55] uses reparametrized PW91 exchange and B95 correlation, but instead of SCS-MP2 it employs the scaled-opposite-spin approach (SOS-MP2) [79] that considers exclusively the correlation contribution of opposite-spin electron pairs.The factor α X in PWPB95 is fixed at 0.50 and the correlation terms are related as c C = 1 − c O , with c C = 0.731 and c O = 0.269.

Conventional Density Functionals
It is instructive to examine first the performance of conventional DFT methods, where the term encompasses GGA, meta-GGA, and hybrid-(meta)-GGA functionals.Table 2 lists exchange coupling constants computed with six representative functionals from these families.The trends are largely anticipated from past studies: functionals that do not incorporate any exact exchange strongly favor the broken-symmetry over the high-spin state, yielding exchange coupling constants that are too negative compared to experiment.In contrast, hybrid functionals stabilize the high-spin states to an extent that depends principally on the admixture of HF exchange.
The BLYP functional performs very poorly, leading to unrealistically strong antiferromagnetic coupling constants for all complexes.The meta-GGA functional TPSS moves to the right direction, albeit without being able to sufficiently correct the GGA bias.The recently proposed SCAN functional is the best performing non-hybrid approach.It still fails to predict ferromagnetic coupling for complex 1, but improves over TPSS for all other complexes.The SCAN results suggest possible non-linearity because the results are rather poor in the weak coupling cases but improve significantly toward the strong exchange-coupling regime (complexes 4 and 5).TPSSh and B3LYP predict correctly the sign of the exchange coupling constants for all complexes.TPSSh, with 10% HF exchange, shows the best overall performance among all methods: it achieves deviations from experiment of merely 3-6 cm −1 for complexes 1, 3, and 4, with the largest deviation being 25 cm −1 for the strongest antiferromagnetic coupling in the test set (complex 5).The good performance of TPSSh documented here is in agreement with previous studies on synthetic manganese complexes and bioinorganic model systems [34][35][36]38,[40][41][42][43][80][81][82][83][84][85].The increase in the percentage of HF exchange to 20% in the B3LYP functional leads to a slight overestimation of the stability of high-spin states and larger deviations from experiment.Only complex 2 appears to be better described by B3LYP compared with TPSSh.Further increase of HF exchange to 25% in PBE0 results in exaggerated ferromagnetic coupling for complex 1, too small antiferromagnetic coupling for 3-5, and qualitatively incorrect reversal of the ground spin state for 2 from low to high spin.The effect of additional diffuse functions [86] in the basis set was tested and found to be negligible (variation of less than 0.2 cm −1 in the computed exchange coupling constants) because they do not have a differential effect on the energies of the high-spin and broken-symmetry solutions.The conductor-like polarizable continuum model (CPCM) [87] was additionally tested with an infinite dielectric in order to investigate possible effects on the computed exchange coupling constants.Compared to the gas-phase results, the CPCM calculations show variations in the J values of less than 2 cm −1 for the antiferromagnetically coupled dimers and up to 9 cm −1 for complex 1.Given that these values were obtained under the extreme assumption of a perfect conductor, it is concluded that the continuum model has only a limited effect on the computed values as it does not strongly favor any particular solution.In conclusion, the various technical aspects of the calculations appear to be converged.In terms of the performance of individual functionals, even though exceptions exist at the quantitative level for specific complexes, TPSSh offers the most balanced performance.

Double-Hybrid Density Functionals
The exchange coupling constants obtained with double-hybrid functionals are listed in Table 3.Compared to the divergent behavior of the functionals discussed above, the variation among double-hybrids appears remarkably small, with the exception of a few erratic results that are discussed below.With the marginal exception of PWPB95, all DHDFs fail to predict antiferromagnetic coupling for complex 2, and all DHDFs overestimate the stability of the high-spin state.
A comparison between Tables 2 and 3 shows that the DHDFs do not replicate the behavior of any conventional functional.For example, B2-PLYP produces the same exchange coupling constant as TPSSh for complex 1, but a qualitatively different result for complex 2, a similar error albeit with opposite sign than TPSSh for complex 3, and significantly greater errors than TPSSh for the more strongly coupled complexes 4 and 5. mPW2-PLYP tracks closely the B2-PLYP results, but with uniformly increased errors for all complexes.The reparametrized versions of B2-PLYP, i.e., the general-purpose B2GP-PLYP and the other two functionals (B2K-PLYP and B2T-PLYP) that were optimized for specific applications similarly show no improvement.Complex 4 yields an outlier for B2K-PLYP, which overestimates the strength of the antiferromagnetic coupling.The two DSD functionals that make use of spin-component-scaled MP2 have similar problems with complex 4, yielding unrealistically large antiferromagnetic exchange coupling constants, particularly in the case of DSD-PBEP86.This is presumably related to the coefficients employed for the SCS-MP2 correction, to the overall MP2 contribution, or both.Moreover, both DSD functionals fail to predict the absolute and relative sign of the coupling for complexes 1 and 2.
PWPB95, which employs only an opposite-spin MP2 component, behaves much more reasonably.It is the only double-hybrid functional that might be considered an improvement over B2-PLYP.This is less in the sense of overall superior numerical agreement with experiment, and more because it is the only double-hybrid functional that marginally predicts weak antiferromagnetic coupling for complex 2.
An important conclusion from the above results is that there are only three well-behaved functionals among the selection used here: B2-PLYP, mPW2-PLYP, and PWPB95.Although the present test set is too small to support entirely general conclusions, the present results suggest that the B2x-PLYP and DSD functionals are likely inappropriate for the problem of exchange coupling.Presumably the parametrization of these functionals for specific uses or with training sets that target narrowly defined properties reduces their generality, resulting in erratic behavior and rendering them unsuitable for exchange-coupled systems.But an equally important conclusion is that despite the fact that the three well-behaved DHDFs do not display the wild fluctuations and inconsistencies observed with conventional functionals (as suggested by their MAD values in Table 3), there is no advantage in choosing any of them over TPSSh for the study of exchange coupling in manganese complexes.

Energetic Contributions to Exchange Coupling from Double-Hybrid Density Functionals
To gain insight into the behavior of the double-hybrid density functionals it is illuminating to decompose the final exchange coupling constant into the energy contribution from the Kohn-Sham orbitals and the perturbational theory component.This is presented in Table 4.
In general, energies based on the Kohn-Sham orbitals without the perturbational contribution are strongly in favor of the high-spin state (J DFT in Table 4).This is due to the high HF admixture in all DHDFs.Looking at these results alone, the usual and expected correlation between the percentage of HF exchange and the stabilization of the high-spin state that was discussed for the conventional hybrid functionals (Table 2) becomes immediately apparent.B2-PLYP, mPW2-PLYP, and PWPB95 appear practically identical here because the results are primarily defined by the fact that all three functionals have almost the same α X factor (0.50-0.55).Reparametrized functionals with higher α X factors yield even greater stabilization of the high-spin state.In terms of the methodological utility of these DHDFs in the prediction of exchange coupling constants, the question is, to what extent can the perturbational contribution correct the flawed J DFT picture.Table 4. Density functional theory (DFT)-only exchange coupling constants J DFT (in cm −1 ) obtained from the Kohn-Sham orbitals of the double hybrid functionals by excluding the perturbational energy component, and the corresponding perturbational contribution (∆J PT2 ) that leads to the final results of Table 3.The ∆J PT2 values in Table 4 are illuminating in this respect.Almost without exception the perturbational component correctly stabilizes the low-spin state, thereby enhancing antiferromagnetic coupling.For the functionals that were termed "well-behaved" above, the perturbational corrections are all very similar.What can be concluded on the basis of these results is that the perturbational correction for these functionals is well-controlled but never sufficient to fully recover the experimentally determined strength of antiferromagnetic coupling.Interestingly, the SOS-MP2 used in PWPB95 gives practically the same corrections as the MP2 component of B2-PLYP and mPW2-PLYP.
The B2x-PLYP functionals are less consistent in the magnitude of the MP2 corrections as a result of their higher and variable α C factors.Thus, in the case of B2K-PLYP (α C = 0.42) the perturbational corrections are more pronounced and it is precisely the MP2 component that causes the problem of a too negative exchange coupling constant for complex 4. The decomposition of the energetic contributions in Table 4 reveals clearly the origin of the weakness of the DSD functionals for the present application.Their weakness stems in large part from the SCS-MP2 approach used in these functionals, or more precisely from the associated c C , c O , and c S parameters that lead to very large and hence unreliable perturbational terms.These two families of DHDFs vividly demonstrate the pitfalls of pursuing property-specific functional parametrizations.
It is beyond the scope of this contribution to devise an improved double-hybrid functional for the present systems.Besides, the test set of complexes used here is too small and chemically restricted to form a meaningful basis for methodological developments of general utility.Nevertheless, the results presented above provide useful insights in this direction.It is suggested that the PWPB95 functional can be a starting point for future refinements: the rather moderate 50% HF exchange and the SOS-MP2 approach with high percentage of DFT correlation (c C = 0.731) appear to provide a reasonable balance between the different energetic components.In the context of BS-DFT applications, the implication is that it could be perfectly sufficient for the perturbational term to correct any deficiencies of the density functional contribution only for opposite-spin correlation.
The high admixture of HF exchange in some of the DHDFs tested here is known to be beneficial for specific applications, for example in predicting reaction barriers.However, it is detrimental for the present problem of exchange coupling, presumably because it places unrealistic demands on the perturbational contribution by creating poor Kohn-Sham orbitals.Several methodological developments not represented by the functionals tested in the present work are currently being pursued in the field of double-hybrid density functionals [52].These include, for example, investigations into the effect of orbital optimization [88,89] and of post-MP2 contributions [90].
Regardless of where such theoretical advances will lead, an important point that needs to be stressed is that the commonly used training or reference data sets, with their focus on thermochemistry and kinetics, are not necessarily testing for the ability of functionals to handle aspects of the electronic structure that relate to spin-state energetics and exchange coupling.To facilitate progress in this area, appropriate benchmark sets must be developed to explicitly target challenging situations characteristic of open-shell transition metal chemistry.Small sets of copper dimers that are sometimes used in the literature [59, 60,91] are insufficient for the purpose of evaluating general-purpose functionals for exchange coupling interactions.Such sets should encompass spin-state energetics of mononuclear transition metal complexes, spin-crossover systems, and diverse cases of magnetic exchange interactions.The present collection of complexes can be part of a large, unbiased, and inclusive future test set that should contain examples spanning many different transition metal ions with different unpaired electron counts.Ideally, of course, a "next generation" of DHDFs would rely less on reference-set parametrization and more on physically motivated improvements that would address fundamental challenges such as that of strong correlation in polynuclear exchange-coupled systems.

Computational Methods
The crystallographic coordinates were taken from the Cambridge Structural Database [92] (see Table 1 for Refcodes of the compounds).Counterions and solvent molecules were removed, hydrogens were added where necessary and their positions were optimized with ORCA [93] using the TPSS functional [71] and the D3 dispersion corrections proposed by Grimme [94].The high-spin state was assumed for these optimizations.Following past practice, the crystallographic positions of heavy atoms were maintained in order to avoid introducing additional errors or bias in our comparison of methods.The zeroth-order regular approximation (ZORA) [95,96] was employed for treating scalar relativistic effects, in combination with ZORA-recontracted versions [97] of the def2-TZVP basis sets [98].Fully decontracted def/J basis sets [99] were used for fitting the Coulomb integrals in the resolution of the identity (RI) approximation.Tight energy convergence criteria, enhanced integration grids (Grid6) and increased radial integration accuracy (IntAcc 6) were used throughout.
The calculation of exchange-coupling constants followed the standard broken-symmetry DFT procedure that utilizes the energy of the high-spin (ferromagnetic) solution and the energy of the broken-symmetry Kohn-Sham determinant that results from "flipping" the spin of the electrons at the metal site with the lowest local spin population in the high-spin solution.All broken-symmetry solutions were verified for correctness by examining the manganese spin populations to confirm the presence of the expected number of unpaired electrons on each metal site (all Mn ions in all complexes have a local high-spin configuration).Explicit inspection is crucial because the assigned spin multiplicity does not necessarily constrain the total number of unpaired electrons in a broken-symmetry calculation.In some cases, the orbital optimization procedure may converge to solutions where the electrons are locally spin-coupled at a metal ion, even though the correct broken-symmetry solution with local high-spin configurations is lower in energy.Such electron-coupled configuration represents an excited state that is not relevant for the exchange-coupling problem.In these cases, it is necessary to restart from the orbitals of the initial high-spin solution and to carefully monitor the self-consistent field procedure in the broken-symmetry calculation, adjusting appropriately the convergence algorithms to avoid changes in orbital occupations and accumulation of numerical errors.Exchange coupling constants were computed with the Yamaguchi formula [25] that smoothly covers weak to strong exchange coupling by scaling the energy difference between the high-spin and broken-symmetry solutions according to their spin expectation values: Methodological details in terms of scalar relativity and all-electron basis sets were the same as discussed above.For hybrid functionals the chain-of-spheres approximation to exact exchange (COSX) [100] was employed with increased grid settings (GridX6 in ORCA convention) to speed up the calculation of exchange integrals without loss of accuracy.The resolution of the identity approach was used in the MP2 (or SCS-MP2) part of the calculation for all double-hybrid functionals, in combination with appropriate correlation fitting def2/C basis sets [101].

Conclusions
A set of dinuclear manganese complexes with various metal oxidation states and a broad range of exchange coupling constants was used to test the performance of representative conventional and double-hybrid density functionals.Among standard approaches, TPSSh was confirmed to be the best choice in line with past studies.Among the double-hybrid functionals, there appears to be little reason to prefer any of the more recent variants over the original B2PLYP.This sobering observation leads to an important conclusion concerning the development of new double-hybrid functionals: the minor or major adjustments and developments reported since 2006 offer no advantage for the problem of exchange coupling.Several double-hybrid density functionals, particularly some that incorporate high admixture of Hartree-Fock exchange or use the spin-component-scaled MP2 approach perform so poorly, that their use is discouraged.Analysis of the different energetic contributions suggests that erratic behavior by double-hybrid functionals is typically associated with large perturbational contributions.The results are not promising for applications to systems of higher nuclearity, where any weaknesses are expected to become even more pronounced.In terms of overall performance, no double-hybrid functional can compete with TPSSh, which remains the method of choice for exchange-coupled manganese systems.In addition to methodological advances that will address physical deficiencies of present density functionals, it is necessary to develop varied reference sets of exchange-coupled transition metal complexes, representative of electronically, chemically, and magnetically diverse systems, on which to assess any new methods.

Figure 1 .
Figure 1.Structures of the manganese complexes included in this study.Hydrogen atoms bound to carbons are omitted for clarity (Mn: purple; C: grey; N: blue; O: red; B: yellow; Cl: green).

Figure 1 .
Figure 1.Structures of the manganese complexes included in this study.Hydrogen atoms bound to carbons are omitted for clarity (Mn: purple; C: grey; N: blue; O: red; B: yellow; Cl: green).

Table 1 .
Dinuclear manganese complexes considered in this study, with their crystallographic identifiers, Mn oxidation states, Mn•••Mn distance R (in Å), and exchange coupling constant J (in cm −1 ).

Table 2 .
Exchange coupling constants J (in cm −1 ) computed with selected conventional density functionals for the five manganese complexes studied in this work, compared with experimentally fitted values.Mean absolute deviations (MAD) in cm −1 .

Table 3 .
Exchange coupling constants J (in cm −1 ) computed with selected double-hybrid density functionals for the five manganese complexes studied in this work, compared with experimentally fitted values.Mean absolute deviations (MAD) in cm −1 .