Multiconfiguration Pair-Density Functional Theory Predicts Spin-State

Apr 4, 2017 - PSL Research University, Institut de Recherche de Chimie Paris ... Department of Chemistry, Chemical Theory Center, and Minnesota Superc...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of Newcastle, Australia

Letter

Multiconfiguration Pair-Density Functional Theory Predicts Spin-State Ordering in Iron Complexes with the Same Accuracy as Complete Active Space Second-Order Perturbation Theory at a Significantly Reduced Computational Cost Liam Wilbraham, Donald G. Truhlar, Laura Gagliardi, and Ilaria Ciofini J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.7b00570 • Publication Date (Web): 04 Apr 2017 Downloaded from http://pubs.acs.org on April 8, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 13

The Journal of Physical Chemistry Letters for JPC Lett., April 4, 2017

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Multiconfiguration Pair-Density Functional Theory Predicts Spin-State Ordering in Iron Complexes with the Same Accuracy as Complete Active Space Second-Order Perturbation Theory at a Significantly Reduced Computational Cost Liam Wilbraham,a Donald G. Truhlar,b* Laura Gagliardi,b*and Ilaria Ciofinia* a

PSL Research University, Institut de Recherche de Chimie Paris IRCP, CNRS – Chimie ParisTech, 11 rue Pierre et Marie Curie, F-75005 Paris, France b Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States

ABSTRACT: The spin-state orderings in nine Fe(II) and Fe(III) complexes with ligands of diverse ligand-field strength were investigated with multiconfiguration pair-density functional theory, MC-PDFT. The performance of this method was compared to that of complete active space second-order perturbation theory, CASPT2, and Kohn-Sham density functional theory, KSDFT. We also investigated the dependence of CASPT2 and MC-PDFT results on the size of the active-space. MC-PDFT reproduces the CASPT2 spin-state ordering, the dependence on the ligand field strength, and the dependence on active space at a computational cost significantly reduced as compared to CASPT2.

TOC image

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

Page 2 of 13

2 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Accurate description of spin-state energetics in molecular complexes containing transition metals is essential for the understanding and prediction of catalytic reactivity,1,2 including biological reactions, and of magnetic properties, such as spin-crossover.3,4 The quantitative calculation of spin-state ordering in transition metal complexes has been highlighted repeatedly as a significant challenge for modern computational methods.5,6 Furthermore, in many cases the key issue is not the quantitative magnitude of the spin state splitting but simply the correct ordering of the spin states, which is important for understanding the reactivity and catalytic power of transition metal complexes, because these properties are strongly dependent on the spin state of the system. Within this context, Kohn-Sham density functional theory (KS-DFT) suffers from sensitivity to the chosen approximation of the exchange-correlation density functional. Local approximations, such as the generalised gradient approximation (GGA), are biased towards lowspin ground states due to their over-delocalization of bonding molecular orbitals.7,8 This can be countered by the admixture of Hartree-Fock (HF) exchange hybrid density functionals, which are known to favor high-spin ground states;9,10 the optimal fraction of HF exchange, however, is system-dependent and may depend on the ligand field strength and,11,14 due to the inherently multiconfiguration nature of many of the spin states, a balanced and controlled description of static and dynamic correlation is required. Consequently some commonly-used functionals, such as B3LYP12 and PBE0,13 sometimes fail to correctly predict the ground spin state. Another method that has been applied to this problem16–18 is the complete active space selfconsistent field15 (CASSCF) method. However, CASSCF calculations recover only a small portion of dynamic correlation energy. Dynamic correlation energy may be added by multireference perturbation theory, such as CASPT2.19 The CASPT2 method, however, is generally limited to small- to medium-sized systems because its computer time and computer memory requirements increase rapidly with system size. The more-recently-developed multiconfiguration pair-density functional theory20,21 (MCPDFT) seeks to combine the explicitly multiconfigurational CASSCF with the low cost of KSDFT; and it provides a promising alternative to CASPT2 in terms of the required computational resources.22 The formulation of MC-PDFT has two key differences from KS-DFT: (i) the reference wave function is multiconfigurational, and (ii) the density functional employed (called the on-top density functional) depends on the total density and the on-top pair density (in contrast to the exchange-correlation density functionals of KS-DFT, which depend on the up-spin and the down-spin densities). Furthermore, MC-PDFT has recently been extended to work with generalised-active-space self-consistent field (GASSCF) wave functions,24,25,26 which have lower cost than CASSCF.

ACS Paragon Plus Environment

Page 3 of 13

The Journal of Physical Chemistry Letters

3 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

In this light, we decided to test the performance of MC-PDFT on the prediction of spinstate ordering. To this end we tested MC-PDFT on nine Fe(II) and Fe(III) complexes with five ligands of diverse ligand-field strength, as shown in Figure 1. We shall compare the performance of this method with both CASPT2 and KS-DFT results for selected density functionals. Additionally, we investigated the dependence of CASPT2 and MC-PDFT on the size of the active-space and on the formulation of the on-top density functional. for the latter purpose we carried out calculations using four on-top density functionals, two of which are based on the PBE27 exchange-correlation functional with the other two based on revPBE28, a modified form of PBE where one parameter is recalibrated based on atomic exchange energies. These exchangecorrelation functionals are converted to on-top functionals by translation protocols explained previously, in particular by the original translation20 that does not include the gradient of the ontop pair density and by the full translation29 that does include it. The original translation is indicated by the prefix “t” and the full translation by the prefix “ft”. The resulting on-top density functionals are tPBE, ftPBE, trevPBE, and ftrevPBE. These functionals are examples of firstgeneration on-top density functionals that, in principle, could be further improved. In the main text we shall discuss results obtained with tPBE and trevPBE, with the fully translated results available in supporting information (SI). For clarity, we shall provide a brief review of MC-PDFT here; a more detailed description can be found elsewhere.20,21 MC-PDFT involves two steps: (1) a multiconfiguration selfconsistent-field calculation and (2) a post-SCF calculation with the on-top density functional. The MC-PDFT energy is , Π = ∑ ℎ + ∑     +  , Π + 

(1)

where  is the MCSCF density; Π is the MCSCF on-top pair density; p, q, r, and s are orbital indices; ℎ are the one-electron kinetic energy, nuclear attraction, and effective core potential integrals, is an element of the MCSCF one-electron density matrix,   are the twoelectron integrals,  is the on-top energy, and  is the internuclear repulsion. The on-top energy, akin to the exchange-correlation energy in KS-DFT, constitutes exchange and correlation additions to the MCSCF classical Coulomb energy (which is given by the four-index sum in eq 1). We note that the correlation energy explicitly accounts for the effect of dynamic correlation to minimize the electron-electron interaction energy, and it also implicitly contains a correction to the MCSCF kinetic energy. Details of the calculations are given in the SI. We calculate the spin-state ordering of the aforementioned complexes using a CASSCF wave function as the MCSCF wave function used as a basis for post-SCF energy calculations using both CASPT2 and MCPDFT. We use basis sets analogous to those previously suggested by Pierloot and co-workers.17 An exploratory

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

Page 4 of 13

4 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

investigation of the effect of active space size was performed, and we shall refer to the active spaces considered as (m/m-1, n), where m [m-1] is the number of active electrons in the active space for Fe(II) [Fe(III)], and n is the number of active orbitals. Five choices of active space are considered (the orbitals are pictured in Figure S2 in the SI): •

(6/5, 5) – Fe(II) valence d orbitals



(10/9, 7) – Fe(II) valence d orbitals plus two Fe-ligand -bonding orbitals



(6/5, 10) – Fe(II) valence d orbitals plus a set correlating d orbitals



(10/9, 12) – (6/5, 10) plus two -bonding, eg-like Fe-ligand orbitals



(12/11, 14) – (10/9, 12) plus Fe 3s and a correlating orbital for 3s.

The largest of these is the active space proposed by Hauser18 and co-workers. For comparison, and to examine the effects of Hartree-Fock exchange on spin-state ordering, we also performed calculations at KS-DFT level using PBE26 and PBE013 functionals, the latter containing 25% Hartree-Fock exchange. In all cases, D2h symmetry was imposed for both geometries and energy calculations (except for complexes 4 and 9, for which Ci symmetry was used). We used symmetry mainly because it retains consistency of the active space across different complexes with ligands of varying ligand-field strength. Generally, we discuss results in terms of the energy difference between high- and low-spin forms of each complex (Δ/ ). Positive (negative) values of Δ/ indicate that the low (high) spin state is energetically favored. We consider the difference in Δ/ for different ligands, and we note that strong-field ligands are expected to produce higher (more positive or less negative) values of Δ/ than weak-field ligands. We consider five different ligands with the ordering expected from the spectrochemical series being H2O < NH3 < NCH < CO < CNH, in terms of Δ/ . Based on diffusion Monte Carlo results from Droghetti31 and on ligand field theory, we expect Fe complexes with CNH and CO will usually have low-spin ground states, and complexes with NCH, NH3 and H2O will usually have high-spin ground states. Figure 2 shows ΔEH/L values calculated by KS-DFT (using PBE and PBE0), by MC-PDFT (using tPBE and trevPBE), and by CASPT2. For MC-PDFT and CASPT2, the (12/11, 14) active space was used. For both Fe(II) and Fe(III) complexes, we note that PBE overestimates ΔEH/L for all ligands with respect to CASPT2. To stress the absolute spin-state ordering obtained in each case, an analogous plot is given in the SI. As anticipated above, based on previous observations for GGA exchange-correlation functionals,11,32 this is interpreted as stemming from an overdelocalization of bonding orbitals, and we find that it severely overstabilizes low-spin states in

ACS Paragon Plus Environment

Page 5 of 13

The Journal of Physical Chemistry Letters

5 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

the case of strong-field ligands. Introducing HF exchange in the form of PBE0 generally stabilizes high-spin states, reducing the value of ΔEH/L for all complexes with respect to PBE. However, the effect of HF exchange is not uniform for ligands of different ligand field strengths, with the induced high-spin stabilization of strong-field ligands (CNH, CO) more severe than for weak field ligands (NCH, NH3, H2O), echoing the previous investigation of the effect of HF exchange on spin-state splitting of octahedral Fe complexes by Ioannidis and Kulik.11This effect is exemplified by the performance of PBE0, which underestimates (overestimates) ΔEH/L for strong-field (weak-field) ligands. The imbalance in the PBE0 treatment of ligands of different ligand field strengths can be rationalized by considering the different types of ligand–metal bonding present in each complex. Ligands that are good π-acceptors (CNH, CO) are most affected by introducing HF exchange. As the self-interaction error32 is reduced in orbitals of π-symmetry, these orbitals are destabilized and the overall ligand-field strength is reduced, resulting in lower values of Δ/ . Conversely, ligands that are good π-donors (NCH, H2O) are less affected due to the increased population of more localized, antibonding d-orbitals upon ligand-to-metal donation. Lastly, ligands which only participate in σ-donation (NH3) are the least affected by the introduction of HF exchange, due to the more local nature of the σ-symmetry bonding orbitals. For both Fe(II) and Fe(III) complexes, the difference in Δ/ computed between each of the two strong-field ligands remains largely constant, irrespective of the functional used. The same holds true for weak-field ligands. This indicates the need for an effective and balanced description of different types of metal-ligand bonding. In analogous work, a recent analysis of the effects of density localization in octahedral metal complexes computed with GGA, global hybrid, and range-separated hybrid density functionals compared with CASPT2 reference values indicated the need for a separate treatment of density delocalization and energetic errors to recover properties that rely on simultaneous treatment of density and orbital energies.33 In this context, we shall now discuss the performance of the MC-PDFT approach compared to CASPT2. Despite the systematic overestimation of Δ/ by tPBE with respect to CASPT2 (Figure 2), we see that the trend of Δ/ values with respect to ligand-field strength at MCPDFT level traces that of CASPT2. Consequently, and in contrast to the results obtained at KSDFT level, the difference in Δ/ as a function of ligand type is largely preserved when using MC-PDFT, demonstrating that using an on-top density functional does not compromise a balanced treatment of strong- and weak-field ligand bonding provided by CASPT2. This is especially striking when one considers that the CASSCF step that precedes the Post-SCF step in both CASPT2 and MC-PDFT does not yield good results, as shown in Fig. S3 of the SI.

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

Page 6 of 13

6 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2 also includes the trevPBE results. We see that the over-stabilization of low-spin states with tPBE is reduced, yet – as before – the trend in Δ/ with respect to ligand field strength is retained. This is consistent with the formulation of the underlying revPBE functional, in which one parameter of the PBE functional has been recalibrated to improve agreement of atomic exchange energies with known exact results. The stability of the good results obtained with MC-PDFT to functional variation extends to other on-top density functionals (ftPBE, ftrevPBE) as well, as seen in the SI. Crucially, a change in the functional formulation shows no difference in the treatment of strong- and weak-field ligands. We also emphasise that these results represent the application of first-generation on-top density functionals and are, in this respect, especially encouraging. The MC-PDFT results are also stable with respect to the size of the chosen active space, as is shown in Figures S3 and S4 of the SI, underlining the even treatment of different metal-ligand interactions by the CASSCF wave function. This is very encouraging. In order to substantiate the claim that MC-PDFT is significantly less computationally expensive than CASPT2, we measured the wall-time taken for each approach. Figure 3 shows the results of this analysis for the Fe(II) complexes, expressed as a ratio between calculation times for CASPT2 and MC-PDFT with respect to active space size. Analogous plots for Fe(III) complexes as well as the full set of absolute wall clock times for each complex are given in the SI. For very small active spaces, we see that there is a negligible increase in speed only in the case of ligands with a reduced number of heavy atoms – i.e. a smaller basis set (H2O, NH3). For all other complexes (CNH, CO, NCH) we observe a significant decrease in calculation time even at the smallest active spaces. The time-saving benefit of MC-PDFT becomes more evident with increasing active space size and for complexes with a larger basis sets. The most striking speed improvement is for [Fe(CO)6]2+, using the (12/11, 14) active space, with MC-PDFT energies for the quintet state calculated approximately 18 times faster than at CASPT2 level. Especially when combining with GASSCF spaces, as mentioned above, we anticipate being able to treat much larger problems with MC-PDFT than with CASPT2. By testing the performance of MC-PDFT when calculating the spin-state ordering of nine Fe(II) and Fe(III) complexes, we have found that this method able to qualitatively reproduce the spin-state ordering calculated at CASPT2 level; the ordering is reproduced in all but one case, and the relative treatment of ligands with diverse ligand field strengths is retained. We conclude that MC-PDFT is able to take advantage of the apparently even-handed even treatment of different ligand-metal bonding types by the underlying CASSCF wave function, and at the same time it benefits from the reduced computational cost associated with calculating exchange and correlation energy using a density functional. MC-PDFT also reacts similarly to CASPT2 as a

ACS Paragon Plus Environment

Page 7 of 13

The Journal of Physical Chemistry Letters

7 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

function of active space size, and shows modest, yet consistent, sensitivity to the on-top density functional.  ASSOCIATED CONTENT Supporting Information. Supporting Information. Computed structural details; active space dependence; on-top density functional dependence; computational details. This material is available free of charge via the Internet at http://pubs.acs.org.  AUTHOR INFORMATION Corresponding Authors *E-mail: [email protected] (IC) *E-mail: [email protected] (LG) *E-mail: [email protected] (DGT) ORCiD

Liam Wilbraham: 0000-0003-4522-6380 Donald G. Truhlar: 0000-0002-7742-7294 Laura Gagliardi: 0000-0001-5227-1396

Ilaria Ciofini: 0000-0002-5391-4522 Notes The authors declare no competing financial interest.  ACKNOWLEDGMENT The authors are grateful to Pragya Verma for helpful interactions. This work was supported in part (L.G. and D.G.T.) by the U. S. National Science Foundation by Grant CHE14-64536. I.C. received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement No 648558, STRIGES CoG grant).

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

Page 8 of 13

8 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. Complexes considered in this study.

Figure 2. ΔEH/L calculated for Fe(II) (left) and Fe(III) (right) complexes in ascending order of computed ligand field strength. Results obtained using PBE, PBE0, MC-PDFT (with two different on-top functionals tPBE and trevPBE), and CASPT2 are shown, with the CASPT2 and MC-PDFT results based on the (12/11, 14) active space. Positive (negative) values of ΔEH/L correspond to a low-spin (high-spin) ground state.

ACS Paragon Plus Environment

Page 9 of 13

The Journal of Physical Chemistry Letters

9 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. Speed increase of MC-PDFT (tPBE) calculations with respect to CASPT2, expressed as the ratio between wall clock times of CASPT2 (tCASPT2) and MC-PDFT (tMC-PDFT) calculations. Squares (triangles) indicate singlet (quintet) spin states. Only results for Fe(II) complexes are shown, with analogous plots for Fe(III) complexes given in the supporting information.

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

Page 10 of 13

10 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC Graphic

ACS Paragon Plus Environment

Page 11 of 13

The Journal of Physical Chemistry Letters

11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

 REFERENCES (1)

Harvey, J. N.; Poli, R.; Smith, K. M. Understanding the Reactivity of Transition Metal Complexes Involving Multiple Spin States. Coord. Chem. Rev. 2003, 238, 347–361.

(2)

Bauer, I.; Knölker, H.-J. Iron Catalysis in Organic Synthesis. Chem. Rev. 2015, 115, 3170–3387.

(3)

Halcrow, M. Structure:function Relationships in Molecular Spin-Crossover Complexes. Chem. Soc. Rev. 2011, 40, 4119–4142.

(4)

Létard, J.-F.; Guionneau, P.; Goux-Capes, L. Towards Spin Crossover Applications. Top. Curr. Chem. 2004, 235, 221-249.

(5)

Kepp, K. P. Consistent Descriptions of Metal-Ligand Bonds and Spin-Crossover in Inorganic Chemistry. Coord. Chem. Rev. 2013, 257, 196-209.

(6)

Cramer, C. J.; Truhlar, D. G. Density Functional Theory for Transition Metals and Transition Metal Chemistry. Phys. Chem. Chem. Phys. 2009, 11, 10757–10816.

(7)

Ganzenmüller, G.; Berkaïne, N.; Fouqueau, A.; Casida, M. E.; Reiher, M. Comparison of Density Functionals for Differences between the High- (5T2g) and Low- (1A1g) Spin States of iron(II) Compounds. IV. Results for the Ferrous Complexes [Fe (L) (NHS4)]. J. Chem. Phys. 2005, 234321.

(8)

Swart, M.; Groenhof, A. R.; Ehlers, A. W.; Lammertsma, K. Validation of Exchange−Correlation Functionals for Spin States of Iron Complexes. J. Phys. Chem. A 2004, 108, 5479–5483.

(9)

Salomon, O.; Reiher, M.; Hess, B. A. Assertion and Validation of the Performance of the B3LYP* Functional for the First Transition Metal Row and the G2 Test Set. J. Chem. Phys. 2002, 117, 4729–4737.

(10)

Reiher, M.; Salomon, O.; Hess, B. A. Reparameterization of Hybrid Functionals Based on Energy Differences of States of Different Multiplicity. Theor. Chem. Acc. 2001, 107, 48– 55.

(11)

Ioannidis, E. I.; Kulik, H. J. Towards Quantifying the Role of Exact Exchange in Predictions of Transition Metal Complex Properties. J. Chem. Phys. 2015, 034104.

(12)

Stephens, P.; Devlin, F.; Chabalowski, C.; Frisch M. J. Ab-initio calculation of vibrational absorption and circular-dichroism spectra using density-functional force-fields. J. Phys. Chem. 1993, 98, 11623–11627.

(13)

Adamo, C.; Barone, V. Toward Reliable Density Functional Methods without Adjustable Parameters: The PBE0 Model. J. Chem. Phys. 1999, 110, 6158.

(14)

Ioannidis, E. I.; Kulik, H. J. Ligand-Field-Dependent Behavior of Meta-GGA Exchange in Transition-Metal Complex Spin-State Ordering. J. Phys. Chem. A 2017, 121, 874-884.

(15)

Roos, B. O.; Taylor, P. R.; Sigbahn, P. E. M. A Complete Active Space SCF Method (CASSCF) Using a Density Matrix Formulated Super-CI Approach. Chem. Phys. 1980, 48, 157–173.

(16)

Kepenekian, M.; Robert, V.; Le Guennic, B.; De Graaf, C. Energetics of [Fe(NCH)6]2+ Via CASPT2 Calculations: A Spin-Crossover Perspective. J. Comput. Chem. 2009, 30, 2327– 2333.

(17)

Pierloot, K.; Vancoillie, S. Relative Energy of the High-(5T2g) and Low-(1A1g) Spin States

ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

Page 12 of 13

12 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

of the Ferrous Complexes [Fe(L)(NHS4)]: CASPT2 versus Density Functional Theory. J. Chem. Phys. 2008, 128, 034104. (18)

Daku, L. M. L.; Aquilante, F.; Robinson, T. W.; Hauser, A. Accurate Spin-State Energetics of Transition Metal Complexes. 1. CCSD(T), CASPT2, and DFT Study of [M(NCH)6]2+ (M= Fe, Co). J. Chem. Theory. Comput. 2012, 8, 4216–4231.

(19)

Andersson, K., Malmqvist, P. A., Roos, B. O., Sadlej, A. J., & Wolinski, K. Second-Order Perturbation Theory with a CASSCF Reference Function. J. Phys. Chem. 1990, 14, 5483– 5488.

(20)

Manni, G. L.; Carlson, R. K.; Luo, S.; Ma, D.; Olsen, J.; Truhlar, D. G.; Gagliardi, L. Multiconfigoration Pair-Density Functional Theory. J. Chem. Theory Comput. 2014, 10, 3669-3680

(21)

Gagliardi, L.; Truhlar, D. G.; Li Manni, G.; Carlson, R. K.; Hoyer, C. E.; Bao, J. L. Multiconfiguration Pair-Density Functional Theory: A New Way To Treat Strongly Correlated Systems. Acc. Chem. Res. 2016, 50, 66-73.

(22)

Sand, A. M.; Truhlar, D. G.; Gagliardi, L. Efficient Algorithm for Multiconfiguration PairDensity Functional Theory with Application to the Heterolytic Dissociation Energy of Ferrocene. J. Chem. Phys. 2017, 146, 034101.

(23)

Odoh, S. O.; Manni, G. L.; Carlson, R. K.; Truhlar, D. G.; Gagliardi, L. Separated-Pair Approximation and Separated-Pair Pair-Density Functional Theory. Chem. Sci. 2016, 7, 2399–2413.

(24)

Ghosh, S.; Cramer, C. J.; Truhlar, D. G.; Gagliardi, L. Generalized-Active-Space PairDensity Functional Theory : An Efficient Method to Study Large , Strongly Correlated , Conjugated Systems. Chem. Sci. 2017, DOI: 10.1039/c6sc05036k.

(25)

Bao, J. L.; Odoh, S. O.; Gagliardi, L.; Truhlar, D. G. Predicting Bond Dissociation Energies of Transition Metal Compounds by Multiconfiguration Pair-Density Functional Theory and Second-Order Perturbation Theory Based on Correlated Participating Orbitals and Separated Pairs. 2017, 16, 1–26.

(26)

Ma, D.; Li Manni, G.; Gagliardi, L. The generalized active space concept in multiconfigurational self-consistent field methods. J. Chem. Phys. 2011, 135, 044128.

(26)

Perdew, J. P.; Burke, K.; Ernzerhof, M.; Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868.

(27)

Zhang, Y.; Yang, W. Comment on “Generalized Gradient Approximation Made Simple.” Phys. Rev. Lett. 1998, 80, 890–890.

(28)

Carlson, R. K.; Truhlar, D. G.; Gagliardi, L. Multiconfiguration Pair-Density Functional Theory: A Fully Translated Gradient Approximation and Its Performance for Transition Metal Dimers and the Spectroscopy of Re2Cl82-. J. Chem. Theory Comput. 2015, 11, 4077– 4085.

(29)

Perdew, J. P.; Savin, A.; Burke, K. Escaping the symmetry dilemma through a pair-density interpretatuon of spin-density functional theory. Phys. Rev. A 1995, 51, 4531.

(30)

Droghetti, A.; Alfè, D.; Sanvito, S. Assessment of Density Functional Theory for iron(II) Molecules across the Spin-Crossover Transition. J. Chem. Phys. 2012, 137, 124303.

(31)

Fouqueau, A.; Casida, M. E.; Lawson Daku, L. M.; Hauser, A.; Neese, F. Comparison of Density Functionals for Energy and Structural Differences between the High-

ACS Paragon Plus Environment

Page 13 of 13

The Journal of Physical Chemistry Letters

13 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

[5T2g:(t2g)4(eg)2] and Low- [1A1g:(t2g)6(eg)0] Spin States of iron(II) Coordination Compounds II. More Functionals and the Hexaminoferrous Cation, [Fe(NH3)6]2+. J. Chem. Phys. 2005, 122, 44110. (32)

Pederson, M. R.; Heaton, R. A.; Lin, C. C. Density-functional theory with self interaction correction: Application to the lithium molecule. J. Chem. Phys. 1985, 82, 2688.

(33) Gani, T. C. H.; Kulik, H. J. Where Does the Density Localize? Convergent Behaviour for Global Hybrids, Range Separation, and DFT+U. J. Chem. Theory Comput. 2016, 12, 5931-5945.

ACS Paragon Plus Environment