Exchange Coupling Interactions from the Density Matrix Renormalization Group and N-Electron Valence Perturbation Theory: Application to a Biomimetic Mixed-Valence Manganese Complex

Combined Multireference-Multiscale Approach to the Description of Photosynthetic Reaction Centers

A first-principles description of the primary photochemical processes that drive photosynthesis and sustain life on our planet remains one of the grand challenges of modern science. Recent research established that explicit incorporation of protein electrostatics in excited-state calculations of photosynthetic pigments, achieved for example with quantum-mechanics/molecular-mechanics (QM/MM) approaches, is essential for a meaningful description of the properties and function of pigment-protein complexes. Although time-dependent density functional theory has been used productively so far in QM/MM approaches for the study of such systems, this methodology has limitations. Here we pursue for the first time a QM/MM description of the reaction center in the principal enzyme of oxygenic photosynthesis, Photosystem II, using multireference wave function theory for the high-level QM region. We identify best practices and establish guidelines regarding the rational choice of active space and appropriate state-averaging for the efficient and reliable use of complete active space self-consistent field (CASSCF) and the N-electron valence state perturbation theory (NEVPT2) in the prediction of low-lying excited states of chlorophyll and pheophytin pigments. Given that the Gouterman orbitals are inadequate as a minimal active space, we define specific minimal and extended active spaces for the NEVPT2 description of electronic states that fall within the Q and B bands. Subsequently, we apply our multireference-QM/MM protocol to the description of all pigments in the reaction center of Photosystem II. The calculations reproduce the electrochromic shifts induced by the protein matrix and the ordering of site energies consistent with the identity of the primary donor (ChlD1) and the experimentally known asymmetric and directional electron transfer. The optimized protocol sets the stage for future multireference treatments of multiple pigments, and hence for multireference studies of charge separation, while it is transferable to the study of any photoactive embedded tetrapyrrole system.

A configuration-based heatbath-CI for spin-adapted multireference electronic structure calculations with large active spaces

This work reports on a spin-pure configuration-based implementation of the heatbath configuration interaction (HCI) algorithm for selective configuration interaction. Besides the obvious advantage of being spin-pure, the presented method combines the compactness of the configurational ansatz with the known efficiency of the HCI algorithm and a variety of algorithmic and conceptual ideas to achieve a high level of performance. In particular, through pruning of the selected configurational space after HCI selection by means of a more strict criterion, a more compact wavefunction representation is obtained. Moreover, the underlying logic of the method allows us to minimize the number of redundant matrix-matrix multiplications while making use of just-in-time compilation to achieve fast diagonalization of the Hamiltonian. The critical search for 2-electron connections within the configurational space is facilitated by a tree-based representation thereof as suggested previously by Gopal et al. Usage of a prefix-based parallelization and batching during the calculation of the PT2-correction leads to a good load balancing and significantly reduced memory requirements for these critical steps of the calculation. In this way, the need for a semistochastic approach to the PT2 correction is avoided even for large configurational spaces. Finally, several test-cases are discussed to demonstrate the strengths and weaknesses of the presented method.

Variational Active Space Selection with Multiconfiguration Pair-Density Functional Theory

The selection of an adequate set of active orbitals for modeling strongly correlated electronic states is difficult to automate because it is highly dependent on the states and molecule of interest. Although many approaches have shown some success, no single approach has worked well in all cases. In light of this, we present the "discrete variational selection" (DVS) approach to active space selection, in which one generates multiple trial wave functions from a diverse set of systematically constructed active spaces and then selects between these wave functions variationally. We apply this DVS approach to 207 vertical excitations of small-to-medium-sized organic and inorganic molecules (with 3 to 18 atoms) in the QUESTDB database by (i) constructing various sets of active space orbitals through diagonalization of parametrized operators and (ii) choosing the result with the lowest average energy among the states of interest. This approach proves ineffective when variationally selecting between wave functions using the density matrix renormalization group (DMRG) or complete active space self-consistent field (CASSCF) energy but is able to provide good results when variationally selecting between wave functions using the energy of the translated PBE (tPBE) functional from multiconfiguration pair-density functional theory (MC-PDFT). Applying this DVS-tPBE approach to selection among state-averaged DMRG wave functions, we obtain a mean unsigned error of only 0.17 eV using hybrid MC-PDFT. This result matches that of our previous benchmark without the need to filter out poor active spaces and with no further orbital optimization following active space selection of the SA-DMRG wave functions. Furthermore, we find that DVS-tPBE is able to robustly and effectively select between the new SA-DMRG wave functions and our previous SA-CASSCF results.

Linearized Pair-Density Functional Theory for Vertical Excitation Energies

Multiconfiguration pair-density functional theory (MC-PDFT) is a computationally efficient method that computes the energies of electronic states in a state specific or state average framework via an on-top functional. However, MC-PDFT does not include state interaction among these states since the final energies do not come from the diagonalization of an effective model-space Hamiltonian. Recently, multistate extensions such as linearized PDFT (L-PDFT) have been developed to accurately model the potentials near conical intersections and avoided crossings. However, there has not been any systematic study evaluating their performance for predicting vertical excitations at the equilibrium geometry of a molecule, when the excited states are generally well separated. In this paper, we report the performance of L-PDFT on the extensive QUESTDB data set of vertical excitations using a database of automatically selected active spaces. We show that L-PDFT performs well on all these excitations and successfully reproduces the performance of MC-PDFT. These results further demonstrate the potential of L-PDFT, as its scaling is constant with the number of states included in the state-average manifold, whereas MC-PDFT scales linearly in this regard.

Local Oxidation States in {FeNO}6-8 Porphyrins: Insights from DMRG/CASSCF-CASPT2 Calculations

A first DMRG/CASSCF-CASPT2 study of a series of paradigmatic {FeNO}6, {FeNO}7, and {FeNO}8 heme-nitrosyl complexes has led to substantial new insight as well as uncovered key shortcomings of the DFT approach. By virtue of its balanced treatment of static and dynamic correlation, the calculations have provided some of the most authoritative information available to date on the energetics of low- versus high-spin states of different classes of heme-nitrosyl complexes. Thus, the calculations indicate low doublet-quartet gaps of 1-4 kcal/mol for {FeNO}7 complexes and high singlet-triplet gaps of ≳20 kcal/mol for both {FeNO}6 and {FeNO}8 complexes. In contrast, DFT calculations yield widely divergent spin state gaps as a function of the exchange-correlation functional. DMRG-CASSCF calculations also help calibrate DFT spin densities for {FeNO}7 complexes, pointing to those obtained from classic pure functionals as the most accurate. The general picture appears to be that nearly all the spin density of Fe[P](NO) is localized on the Fe, while the axial ligand imidazole (ImH) in Fe[P](NO)(ImH) pushes a part of the spin density onto the NO moiety. An analysis of the DMRG-CASSCF wave function in terms of localized orbitals and of the resulting configuration state functions in terms of resonance forms with varying NO(π*) occupancies has allowed us to address the longstanding question of local oxidation states in heme-nitrosyl complexes. The analysis indicates NO(neutral) resonance forms [i.e., Fe(II)-NO0 and Fe(III)-NO0] as the major contributors to both {FeNO}6 and {FeNO}7 complexes. This finding is at variance with the common formulation of {FeNO}6 hemes as Fe(II)-NO+ species but is consonant with an Fe L-edge XAS analysis by Solomon and co-workers. For the {FeNO}8 complex {Fe[P](NO)}-, our analysis suggests a resonance hybrid description: Fe(I)-NO0 ↔ Fe(II)-NO-, in agreement with earlier DFT studies. Vibrational analyses of the compounds studied indicate an imperfect but fair correlation between the NO stretching frequency and NO(π*) occupancy, highlighting the usefulness of vibrational data as a preliminary indicator of the NO oxidation state.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Systematic determination of coupling constants in spin clusters from broken-symmetry mean-field solutions

Quantum-chemical calculations aimed at deriving magnetic coupling constants in exchange-coupled spin clusters commonly utilize a broken-symmetry (BS) approach. This involves calculating several distinct collinear spin configurations, predominantly by density-functional theory. The energies of these configurations are interpreted in terms of the Heisenberg model, H̃=∑i Collapse

Key Words Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Shadan Ghassemi Tabrizi Department of Chemistry, University of Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany

Collapse

High-Spin Blatter's Triradicals

Robust organic triradicals with high-spin quartet ground states provide promising applications in molecular magnets, spintronics, etc. In this context, a triradical based on Blatter's radical has been synthesized recently, having two low-lying non-degenerate doublet states with a quartet ground state. The traditional broken-symmetry (BS)-DFT computed doublet-quartet energy gaps are reported to be somewhat overestimated in comparison to the experimentally observed values. In this work, we have employed different ab initio methods on this prototypical system to obtain more accurate doublet-quartet energy gaps for this triradical. The spin-constraint broken-symmetry (CBS)-DFT method has been used to reduce the overestimation of energy gaps from BS-DFT. To address the issues of spin-contamination and the multireference nature of low-spin states affecting the DFT methods, we have computed the energy gaps using appropriately state-averaged CASSCF and NEVPT2 computations. Using a series of active spaces, our calculations are shown to provide quite accurate values in concordance with the experimentally observed results. Furthermore, we have proposed and modeled another two triradicals based on Blatter's radical, which are of interest for experimental synthesis and characterization. Our computations show that all these triradicals also have a quartet ground state with a similar energy difference between the excited doublet states.

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r2SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew-Burke-Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r2SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

Magnetic Interactions in a [Co(II)3Er(III)(OR)4] Model Cubane through Forefront Multiconfigurational Methods

Strong electron correlation effects are one of the major challenges in modern quantum chemistry. Polynuclear transition metal clusters are peculiar examples of systems featuring such forms of electron correlation. Multireference strategies, often based on but not limited to the concept of complete active space, are adopted to accurately account for strong electron correlation and to resolve their complex electronic structures. However, transition metal clusters already containing four magnetic centers with multiple unpaired electrons make conventional active space based strategies prohibitively expensive, due to their unfavorable scaling with the size of the active space. In this work, forefront techniques, such as density matrix renormalization group (DMRG), full configuration interaction quantum Monte Carlo (FCIQMC), and multiconfiguration pair-density functional theory (MCPDFT), are employed to overcome the computational limitation of conventional multireference approaches and to accurately investigate the magnetic interactions taking place in a [Co(II)₃Er(III)(OR)₄] (chemical formula [Co(II)₃Er(III)(hmp)₄(μ₂-OAc)₂(OH)₃(H₂O)], hmp = 2-(hydroxymethyl)-pyridine) model cubane water oxidation catalyst. Complete active spaces with up to 56 electrons in 56 orbitals have been constructed for the seven energetically lowest different spin states. Relative energies, local spin, and spin-spin correlation values are reported and provide crucial insights on the spin interactions for this model system, pivotal in the rationalization of the catalytic activity of this system in the water-splitting reaction. A ferromagnetic ground state is found with a very small, ∼50 cm^-1, highest-to-lowest spin gap. Moreover, for the energetically lowest states, S = 3-6, the three Co(II) sites exhibit parallel aligned spins, and for the lower states, S = 0-2, two Co(II) sites retain strong parallel spin alignment.

Reimagining the Wave Function in Density Functional Theory: Exploring Strongly Correlated States in Pancake-Bonded Radical Dimers

Pancake bonding between π-conjugated radicals challenges conventional electronic structure approximations, due to the presence of both dispersion (van der Waals) interactions and "strong" electron correlation. Here we use a reimagined wave function-in-density functional theory (DFT) approach to model pancake bonds. Our generalized self-interaction correction extends DFT's reference system of noninteracting electrons, by introducing electron-electron interactions within an active space. We show that a small variation on our previous derivation recovers a DFT-corrected complete active space method proposed by Pijeau and Hohenstein. Comparison of the two approaches shows that the latter provides reasonable dissociation curves for single bonds and pancake bonds, including excited states inaccessible to conventional linear response time-dependent DFT. The results motivate broader adoption of wavefunction-in-DFT approaches for modeling pancake bonds.

Magneto-structural maps and bridged-ligand effect for dichloro-bridged dinuclear copper(ii) complexes: a theoretical perspective

Theoretical understanding of magneto-structural correlations in dichloro-bridged dicopper(ii) complexes can guide the design of magnetic materials having broad-scale applications. However, previous reports suggest these correlations are complicated and unclear. To clarify possible correlations, magnetic coupling constants (J _calc) of variants of a representative {Cu-(μ-Cl)₂-Cu} complex A were calculated through BS-DFT. The variation of the Cu-(μ-Cl)-Cu angle (α), Cu⋯Cu distance (R ₀), and Cu-Cl-Cu-Cl dihedral angle (τ) followed by structural optimization and calculation of the magnetic coupling constant (J _calc) revealed several trends. J _calc increased linearly with R ₀ and τ, and initially increased and then decreased with α. Further, bridging ligand effects on J _calc for dicopper(ii) complexes were evaluated through BS-DFT; the results revealed that J _calc increased with increasing ligand field strength (I^- < Br^- < Cl^- < N₃ ^- < F^-). Furthermore, a linear relationship was found between the spin density of the bridging ligand and J _calc.

Understanding the Exchange Interaction between Paramagnetic Metal Ions and Radical Ligands: DFT and Ab Initio Study on Semiquinonato Cu(II) Complexes

The exchange coupling, represented by the J parameter, is of tremendous importance in understanding the reactivity and magnetic behavior of open-shell molecular systems. In the past, it was the subject of theoretical investigations, but these studies are mostly limited to the interaction between metallic centers. The exchange coupling between paramagnetic metal ions and radical ligands has hitherto received scant attention in theoretical studies, and thus the understanding of the factors governing this interaction is lacking. In this paper, we use DFT, CASSCF, CASSCF/NEVPT2, and DDCI3 methods to provide insight into exchange interaction in semiquinonato copper(II) complexes. Our primary objective is to identify structural features that affect this magnetic interaction. We demonstrate that the magnetic character of Cu(II)-semiquinone complexes are mainly determined by the relative position of the semiquinone ligand to the Cu(II) ion. The results can support the experimental interpretation of magnetic data for similar systems and can be used for the in-silico design of magnetic complexes with radical ligands.

Consistent Evaluation of Magnetic Exchange Couplings in Multicenter Compounds in KS-DFT: The Recomposition Method

The use of broken-symmetry calculations in Kohn-Sham density functional theory has offered an affordable route to study magnetic exchange couplings in transition-metal-based compounds. However, computing this property in compounds exhibiting several couplings is still challenging and especially due to the difficulties to overcome the well-known problem of spin contamination. Here, we present a new and general method to compute magnetic exchange couplings in systems featuring several spin sites. To provide a consistent spin decontamination of J values, our strategy exploits the decomposition method of the magnetic exchange coupling proposed by Coulaud et al. and generalizes our previous work on diradical compounds where the overall magnetic exchange coupling is defined as the sum of its three main and properly extracted physical contributions (direct exchange, kinetic exchange, and spin polarization). In this aim, the generalized extraction of all contributions is presented to systems with multiple spin sites bearing one unpaired electron. This is done by proposing a new paradigm to treat the kinetic exchange contribution, which proceeds through monorelaxations of the magnetic orbitals. This method, so-called the recomposition method, is applied to a compound featuring three Cu(II) ions with a linear arrangement and to a recently synthesized complex containing a Cu₄O₄ cubane unit presenting an unusual magnetic behavior.

A Theoretical Account of the Coupling between Metal- and Ligand-centred Spins

This study addresses the magnetic interaction between paramagnetic metal ions and the radical ligands taking the [Cu^II (hfac)₂ (imVDZ)] and [M^II (hfac)₂ (pyDTDA)] (imVDZ=1,5-dimethyl-3-(1-methyl-2-imidazolyl)-6-oxoverdazyl; hfac=(1,1,1,5,5,5)hexafluroacetylacetonate; pyDTDA=4-(2'-pyridyl)-1,2,3,5-dithiadiazolyl), (M=Cu, Ni, Co, Fe, Mn) compounds as reference systems. The coupling between the metal and ligand spins is quantified in terms of the exchange coupling constant (J) in the platform of density functional theory (DFT) and the wave function-based complete active space self-consistent field (CASSCF) method. Application of DFT and broken symmetry (BS) formalism results ferromagnetic coupling for all the transition metal complexes except the Mn(II) complex. This DFT-BS prediction of magnetic nature matches with the experimental finding for all the complexes other than the Fe(II)-pyDTDA complex, for which an antiferromagnetic coupling between high spin iron and the thiazyl ligand has been reported. However, evaluation of spin state energetics through the multiconfigurational wave function-based method produces the S=3/2 ground spin state for the iron-thiazyl in parity with experiment. Electronic structure analyses find the overlap between the metal- and ligand-based singly occupied molecular orbitals (SOMOs) to be one of the major reasons attributing to different extent of exchange coupling in the systems under investigation.

Water oxidation in oxygenic photosynthesis studied by magnetic resonance techniques

The understanding of light-induced biological water oxidation in oxygenic photosynthesis is of great importance both for biology and (bio)technological applications. The chemically difficult multistep reaction takes place at a unique protein-bound tetra-manganese/calcium cluster in photosystem II whose structure has been elucidated by X-ray crystallography (Umena et al. Nature 2011, 473, 55). The cluster moves through several intermediate states in the catalytic cycle. A detailed understanding of these intermediates requires information about the spatial and electronic structure of the Mn₄ Ca complex; the latter is only available from spectroscopic techniques. Here, the important role of Electron Paramagnetic Resonance (EPR) and related double resonance techniques (ENDOR, EDNMR), complemented by quantum chemical calculations, is described. This has led to the elucidation of the cluster's redox and protonation states, the valence and spin states of the manganese ions and the interactions between them, and contributed substantially to the understanding of the role of the protein surrounding, as well as the binding and processing of the substrate water molecules, the O-O bond formation and dioxygen release. Based on these data, models for the water oxidation cycle are developed.

Localized Active Space-State Interaction: a Multireference Method for Chemical Insight

Multireference electronic structure methods, like the complete active space (CAS) self-consistent field model, have long been used to characterize chemically interesting processes. Important work has been done in recent years to develop modifications having a lower computational cost than CAS, but typically these methods offer no more chemical insight than that from the CAS solution being approximated. In this paper, we present the localized active space-state interaction (LASSI) method that can be used not only to lower the intrinsic cost of the multireference calculation but also to improve interpretability. The localized active space (LAS) approach utilizes the local nature of the electron-electron correlation to express a composite wave function as an antisymmetrized product of unentangled wave functions in local active subspaces. LASSI then uses these LAS states as a basis from which to express complete molecular wave functions. This not only makes the molecular wave function more compact but also permits flexibility in choosing those states to be included in the basis. Such selective inclusion of states translates to the selective inclusion of specific types of interactions, thereby allowing a quantitative analysis of these interactions. We demonstrate the use of LASSI to study charge migration and spin-flip excitations in multireference organic molecules. We also compute the J coupling parameter for a bimetallic compound using various LAS bases to construct the Hamiltonian to provide insights into the coupling mechanism.

Computational design of magnetic molecules and their environment using quantum chemistry, machine learning and multiscale simulations

Having served as a playground for fundamental studies on the physics of d and f electrons for almost a century, magnetic molecules are now becoming increasingly important for technological applications, such as magnetic resonance, data storage, spintronics and quantum information. All of these applications require the preservation and control of spins in time, an ability hampered by the interaction with the environment, namely with other spins, conduction electrons, molecular vibrations and electromagnetic fields. Thus, the design of a novel magnetic molecule with tailored properties is a formidable task, which does not only concern its electronic structures but also calls for a deep understanding of the interaction among all the degrees of freedom at play. This Review describes how state-of-the-art ab initio computational methods, combined with data-driven approaches to materials modelling, can be integrated into a fully multiscale strategy capable of defining design rules for magnetic molecules.

Exploring optimal multimode vibronic pathways in singlet fission of azaborine analogues of perylene

The development of new singlet fission chromophores is a vibrant area of research to explore the possibility of efficient photovoltaic devices. Using high-level ab-initio density matrix renormalization group calculations, we present a systematic analysis of BN-doped perylenes for their potential application as singlet fission candidates. Four singlet fission chromophores are identified considering the monomer-based properties and their excitonic characters are further analyzed in a dimer configuration optimized in a six-dimensional space for local maxima of fission rates. Furthermore, a multistate multimode vibronic Hamiltonian is employed to identify intra- and interstate vibrational pathways for excitation energy modulation. Several photophysical properties such as Davydov splitting, activation energy and vibronic admixture of multiexcitonic and charge-transfer states are calculated for physically accessible dimers. The optimal dimer packing results in appropriate vibrational relaxation of singlet fission states and promotes significant population transfer which would be more attenuated without such couplings. This work not only identifies potential singlet fission systems with favorable electronic properties but also highlights the sensitivity of dimer packings with respect to the substitution patterns in singlet fission chromophores.

Harnessing Colossal Magnetic Anisotropy in Sandwiched 3d2-Metallocenes

Single-molecule magnets are gaining attention in recent years with the growing focus on achieving higher barriers of magnetization reversal. Metallocenes, owing to their unique sandwiched structure, assure themselves as plausible molecular systems for the development of novel single-molecule magnets (SMMs). Here in this work, we have explicitly investigated metallocenes of first-row transition elements, along with their one-electron-oxidized (cationic) and -reduced (anionic) analogues, for their magnetic anisotropies by adopting multireference ab initio calculations. Herein, we report a high magnetic anisotropy for 3d² systems among all 3d-metallocenes.

On the fly swapping algorithm for ordering of degrees of freedom in density matrix renormalization group

Density matrix renormalization group (DMRG) and its time-dependent variants have found widespread applications in quantum chemistry, includingab initioelectronic structure of complex bio-molecules, spectroscopy for molecular aggregates, and charge transport in bulk organic semiconductors. The underlying wavefunction ansatz for DMRG, matrix product state (MPS), requires mapping degrees of freedom (DOF) into a one-dimensional topology. DOF ordering becomes a crucial factor for DMRG accuracy. In this work, we propose swapping neighboring DOFs during the DMRG sweeps for DOF ordering, which we term 'on the fly swapping' (OFS) algorithm. We show that OFS is universal for both static and time-dependent DMRG with minimum computational overhead. Examples are given for one dimensional antiferromagnetic Heisenberg model,ab initioelectronic structure of N₂molecule, and the S₁/S₂internal conversion dynamics of pyrazine molecule. It is found that OFS can indeed improve accuracy by finding better DOF ordering in all cases.

Post-Density Matrix Renormalization Group Methods for Describing Dynamic Electron Correlation with Large Active Spaces

The ab initio density matrix renormalization group (DMRG) method has been well-established and has become one of the most accurate numerical methods for the precise electronic structure solution of large active spaces. In the past few years, to capture the missing dynamic correlation, various post-DMRG approaches have been proposed through the combination of DMRG and multireference quantum chemical methods or density functional theory. With this in mind, this work provides a brief overview of ab initio DMRG principles and the new developments within post-DMRG methods. For clarity, post-DMRG methods are classified into two main categories depending on whether high-order n-electron reduced density matrices are used, and their merits and disadvantages are properly discussed. Finally, we conclude by discussing unsolved bottlenecks and giving development perspectives of post-DMRG approaches, which are expected to yield quantitative descriptions of complex electronic structures in large strongly correlated molecules and materials.

Antiferromagnetic Coupling Supported by Metallophilic Interactions: Theoretical View

The antiferromagnetic coupling supported by metallophilic interactions has been studied in the framework of the broken symmetry approach (BS) and multiconfigurational calculations (CASSCF). A series of heterobimetallic complexes of the form [PtCo(X)₄(Y)]₂ (X = tba thiobenzoate, SAc thioacetate, and Y = H₂O, NO₂py, py), previously reported, have been used as model systems. Magnetic coupling constants were found in good agreement with the experimental reports, and it could be concluded that axial ligands with a pure σ-donor character have a marked effect on the J value strengthening the antiferromagnetic coupling, as shown for [PtCo(SAc)₄(H₂O)]₂ and [PtNi(SAc)₄(H₂O)]₂. The latter complex, included for comparative purposes, also made it possible to evidence that the interaction between magnetic orbitals and low-level excitation in the Pt···Pt region is also relevant favoring the stronger antiferromagnetic coupling found in this case. A careful analysis of the energetic components involved in Pt···Pt interaction suggests that the stabilization arises from a combination of favorable orbital contributions, which allows a weak covalent Pt···Pt σ(dz²...dz²) bond. Theoretical tools evidence that the weak σ-bond found between monomeric units is responsible for a spin polarization mechanism resulting in the observed antiferromagnetic interaction. Multiconfigurational calculations finally allowed us to establish that the spin polarization mechanism involves not only the d_z² orbitals in the M-Pt···Pt-M bond direction but also the empty 6p_z orbitals of Pt atoms. The inclusion of these orbitals favors a correlation-induced delocalization of magnetic orbitals and therefore a better balance among direct and kinetic exchange. The results shown in this work are relevant in the molecular design of systems supported by metallophilic interactions not only between platinum atoms but also could be extended to other cases with similar interactions.

A DMRG/CASPT2 Investigation of Metallocorroles: Quantifying Ligand Noninnocence in Archetypal 3d and 4d Element Derivatives

Hybrid density functional theory (B3LYP) and density matrix renormalization group (DMRG) theory have been used to quantitatively compare the degree of ligand noninnocence (corrole radical character) in seven archetypal metallocorroles. The seven complexes, in decreasing order of corrole noninnocent character, are Mn[Cor]Cl > Fe[Cor]Cl > Fe[Cor](NO) > Mo[Cor]Cl2 > Ru[Cor](NO) ≈ Mn[Cor]Ph ≈ Fe[Cor]Ph ≈ 0, where [Cor] refers to the unsubstituted corrolato ligand. DMRG-based second-order perturbation theory calculations have also yielded detailed excited-state energetics data on the compounds, shedding light on periodic trends involving middle transition elements. Thus, whereas the ground state of Fe[Cor](NO) (S = 0) is best described as a locally S = 1/2 {FeNO}7 unit antiferromagnetically coupled to a corrole A' radical, the calculations confirm that Ru[Cor](NO) may be described as simply {RuNO}6-Cor3-, that is, having an innocent corrole macrocycle. Furthermore, whereas the ferromagnetically coupled S = 1{FeNO}7-Cor•2- state of Fe[Cor](NO) is only ∼17.5 kcal/mol higher than the S = 0 ground state, the analogous triplet state of Ru[Cor](NO) is higher by a far larger margin (37.4 kcal/mol) relative to the ground state. In the same vein, Mo[Cor]Cl2 exhibits an adiabatic doublet-quartet gap of 36.1 kcal/mol. The large energy gaps associated with metal-ligand spin coupling in Ru[Cor](NO) and Mo[Cor]Cl2 reflect the much greater covalent character of 4d-π interactions relative to analogous interactions involving 3d orbitals. As far as excited-state energetics is concerned, DMRG-CASPT2 calculations provide moderate validation for hybrid density functional theory (B3LYP) for qualitative purposes, but underscore the possibility of large errors (>10 kcal/mol) in interstate energy differences.

Magnetic exchange coupling in Cu dimers studied with modern multireference methods and broken-symmetry coupled cluster theory

Spin-state energetics of exchange-coupled copper complexes pose a persistent challenge for applied quantum chemistry. Here, we provide a comprehensive comparison of all available theoretical approaches to the problem of exchange coupling in two antiferromagnetically coupled bis-μ-hydroxo Cu(II) dimers. The evaluated methods include multireference methods based on the density matrix renormalization group (DMRG), multireference methods that incorporate dynamic electron correlation either perturbatively, such as the N-electron valence state perturbation theory, or variationally, such as the difference-dedicated configuration interaction. In addition, we contrast the multireference results with those obtained using broken-symmetry approaches that utilize either density functional theory or, as demonstrated here for the first time in such systems, a local implementation of coupled cluster theory. The results show that the spin-state energetics of these copper dimers are dominated by dynamic electron correlation and represent an impossible challenge for multireference methods that rely on brute-force expansion of the active space to recover correlation energy. Therefore, DMRG-based methods even at the limit of their applicability cannot describe quantitatively the antiferromagnetic exchange coupling in these dimers, in contrast to dinuclear complexes of earlier transition metal ions. The convergence of the broken-symmetry coupled cluster approach is studied and shown to be a limiting factor for the practical application of the method. The advantages and disadvantages of all approaches are discussed, and recommendations are made for future developments.

Modern multireference methods and their application in transition metal chemistry

Transition metal chemistry is a challenging playground for quantum chemical methods owing to the simultaneous presence of static and dynamic electron correlation effects in many systems. Wavefunction based multireference (MR) methods constitute a physically sound and systematically improvable Ansatz to deal with this complexity but suffer from some conceptual difficulties and high computational costs. The latter problem partially arises from the unfavorable scaling of the Full Configuration Interaction (Full-CI) problem which in the majority of MR methods is solved for a subset of the molecular orbital space, the so-called active space. In the last years multiple methods such as modern variants of selected CI, Full-CI Quantum Monte Carlo (FCIQMC) and the density matrix renormalization group (DMRG) have been developed that solve the Full-CI problem approximately for a fraction of the computational cost required by conventional techniques thereby significantly extending the range of applicability of modern MR methods. This perspective review outlines recent advancements in the field of MR electronic structure methods together with the resulting chances and challenges for theoretical studies in the field of transition metal chemistry. In light of its emerging importance a special focus is put on the selection of adequate active spaces and the concomitant development of numerous selection aides in recent years.

Improved Description and Efficient Implementation of Spin-Projected Perturbation Theory for Practical Applications

In this study, we continue to develop the recently proposed second-order perturbation theory for the spin-projected Hartree-Fock method [Tsuchimochi, T.; Ten-no, S. L. J. Chem. Theory Comput. 2019, 15, 6688] in various aspects. A new, stable imaginary level-shift scheme is derived to obtain a well-conditioned equation, enabling a significantly faster convergence. To achieve a further speed-up, we propose a preconditioning scheme considering the pair character on a spin-projected basis. We also eliminate the computational memory bottleneck in solving the linear equation for large systems using a distributed memory parallel implementation. Finally, for the description of open-shell molecules, several modified zeroth-order Hamiltonians are introduced and tested using the Mn₂O₂(NHCHCO₂)₄ complex. These developments enable practical calculations of a second-order perturbation theory with improved accuracy at a reduced computational cost.

Magnetic Anisotropy in a Cubane-like Ni4 Complex: An Ab Initio Perspective

Magnetic anisotropy, in the absence of an external magnetic field, relates to the degeneracy lift of energy levels. In the standard case of transition metal complexes, this property is usually modeled by an anisotropic spin Hamiltonian and one speaks of "zero-field splitting" (ZFS) of spin states. While the case of mononuclear complexes has been extensively described by means of ab initio quantum mechanical calculations, the literature on polynuclear complexes studied with these methodologies is rather scarce. In this work, advanced multiconfigurational wave function theory methods are applied to compute the ZFS of the ground S = 4 state of an actual tetranickel(II) complex, displaying a magnet behavior below 0.5 K. First, the isotropic couplings are computed in the absence of the spin-orbit coupling operator, in the full complex and also in clusters with only two active nickel(II) centers, confirming the occurrence of weak ferromagnetic couplings in this system. Second, the single-site magnetic anisotropies are computed on a cluster bearing only one active nickel(II) site, showing that the single-site anisotropy axes are not oriented in an optimal fashion for generating a large uniaxial molecular anisotropy. Furthermore, the possibility for involving only a few local orbital excited states in the calculation is assessed, actually opening the way for a consistent and manageable treatment of the ZFS of the ground S = 4 state. Third, multiconfigurational calculations are performed on the full complex, confirming the weak uniaxial anisotropy occurring for this state and also, interestingly, revealing a significant contribution of the lowest-lying orbitally excited S = 3 states. Overall, by comparison with the experiment, the reported results question the common habit of using only one structure, in particular derived from a crystallography experiment, to compute magnetic anisotropy parameters.

First-Principles Investigations of Magnetic Anisotropy and Spin-Crossover Behavior of Fe(III)-TBP Complexes

With the ongoing effort to obtain mononuclear 3d-transition-metal complexes that manifest slow relaxation of magnetization and, hence, can behave as single-molecule magnets (SMMs), we have modeled 14 Fe(III) complexes based on an experimentally synthesized (PMe₃)₂FeCl₃ complex [J. Am. Chem. Soc. 2017, 139 (46), 16474-16477], by varying the axial ligands with group XV elements (N, P, and As) and equatorial halide ligands from F, Cl, Br, and I. Out of these, nine complexes possess large zero field splitting (ZFS) parameter D in the range of -40 to -60 cm^-1. The first-principles investigation of the ground-spin state applying density functional theory (DFT) and wave function-based multiconfigurations methods, e.g., SA-CASSCF/NEVPT2, are found to be quite consistent except for few delicate cases with near-degenerate spin states. In such cases, the hybrid B3LYP functional is found to be biased toward high-spin (HS) state. Altering the percentage of exact exchange admixed in the B3LYP functional leads to intermediate-spin (IS) ground state consistent with the multireference calculations. The origin of large zero field splitting (ZFS) in the Fe(III)-based trigonal bipyramidal (TBP) complexes is investigated. Furthermore, a number of complexes are identified with very small ΔG_HS-IS^adia. values indicating the possible spin-crossover phenomenon between the bistable spin states.

Beyond the Classical Contributions to Exchange Coupling in Binuclear Transition Metal Complexes

Complexes with two or more magnetically coupled metal ions have attracted considerable attention as catalysts of many vital processes, single-molecule magnets, or spin-crossover compounds. Elucidation of their electronic structures is essential for understanding their catalytic and magnetic properties. Here, we provide an unprecedented insight into exchange-coupling mechanisms between the magnetic centers in six prototypical bis-μ-oxo bimetallic M₂O₂ complexes, including two biologically relevant models of non-heme iron enzymes. Employing multiconfigurational/multireference methods and related orbital entanglement analysis, we revealed the essential and counterintuitive role of predominantly unoccupied valence metal d orbitals in their strong antiferromagnetic coupling. We found that the participation of these orbitals is twofold. First, they enhance the superexchange between the singly occupied d orbitals. Second, they become substantially occupied and thus directly magnetically active, which we perceive as a new mechanism of the exchange interaction between the magnetic transition metal centers.

Interpretation of Exchange Interaction through Orbital Entanglement

Recently, the analysis of single-orbital entropy and mutual information has been introduced as a tool for the investigation of contributions to the exchange (J) coupling between open-shell metal ions [Stein et al. J. Phys. Chem. Lett. 2019, 10, 6762-6770]. Here, we show that this analysis may lead to an incorrect interpretation of the J-coupling mechanism. Instead, we propose an orbital-entanglement analysis that is based on the two-electron density and that provides a coherent picture of the contributing exchange pathways, which seems fully consistent with the available J values. For this purpose, we used a prototypical bis-μ-oxo binuclear manganese complex ([Mn₂O₂(NH₃)₈]⁴⁺) and demonstrated that its antiferromagnetism (J < 0), calculated by using the active space composed of all valence p_O and d_Mn orbitals, correlates well with the largest elements in the differential low-spin vs high-spin entanglement map. These elements correspond to interactions between the pairs of d_Mn orbitals mediated by the oxo-bridging out-of-plane p orbitals, representing the π superexchange pathway. We also show that the reduction of active space to manifold of the singly occupied magnetic orbitals does not lead to discrepancy between the calculated J values and entanglement maps. This contrasts to analysis of mutual information, which suggests the "direct" d_Mn-d_Mn interactions to play a dominant role for the J coupling, irrespective of the size of active space as well as of the antiferromagnetism expected. The failure is attributed to the large contribution of spin entanglement contained in the mutual information of the low-spin state, which may be regarded as the origin of the different complexity of its wave function and electron density.

Organic Mixed-Valence Compounds and the Overhauser Effect in Insulating Solids

Recent experiments have shown that the organic free radical 1,3-bisdiphenylene-2-phenylallyl (BDPA) can induce an Overhauser effect dynamic nuclear polarization in insulating solids, a feat previously considered not to be possible. Here, we establish that this peculiar ability of the BDPA radical stems from its mixed-valence nature and the ensuing intramolecular charge transfer. Using state-of-the-art DMRGSCF calculations, we confirm the class II mixed-valence nature of BDPA with the characteristic double-well potential energy surface, and we investigate the mechanism of the consequent electron hopping. A two-component vibronic Hamiltonian is then employed to compute the rate of electron hopping from a quantum dynamical time-propagation of the density matrix. The predicted hyperfine coupling oscillations indeed fall within the frequency range required for an Overhauser effect. The paradigm of mixed-valence compounds as a mining source opens many possibilities for the development and fine tuning of novel polarizing agents.

Successes, challenges, and opportunities for quantum chemistry in understanding metalloenzymes for solar fuels research

Overview of the rich and diverse contributions of quantum chemistry to understanding the structure and function of the biological archetypes for solar fuel research, photosystem II and hydrogenases.

Extending the ASS1ST Active Space Selection Scheme to Large Molecules and Excited States

Multireference electronic structure methods based on the CAS (complete active space) ansatz are well-established as a means to provide reliable predictions of physical properties of strongly correlated systems. A critical aspect of every CAS calculation is the selection of an adequate active space, in particular as the boundaries for tractable active spaces have been shifted significantly with the emergence of efficient approximations to the Full-CI problem like the density matrix renormalization group and full-CI quantum Monte Carlo. Recently, we proposed an active space selection based on first-order perturbation theory (ASS1ST) that yields satisfactory results for the electronic ground state of a variety of strongly correlated systems. In this work, we present a state-averaged extension of ASS1ST (SA-ASS1ST) that determines suitable active spaces when electronically excited states are targeted. Furthermore, the computational costs of the single state and state-averaged variants are significantly reduced by a simple approximation that avoids the most expensive step of the original method, the evaluation of active space four-electron reduced density matrices, altogether. After the applicability of the approximation is established, test calculations on a biomimetic Mn₄O₄ cluster demonstrate the enhanced range of ASS1ST in terms of system size and complexity. Furthermore, calculations on [VOCl₄]₂^-, MeMn(CO)₃-α-diimine, and anthracene show that SA-ASS1ST suggests well-suited active spaces to describe d → d and charge-transfer excitations in transition-metal complexes as well as π → π* excitations in aryl compounds. Finally, the application of ASS1ST on multiple points of the potential energy surface of Cr₂ illustrates the applicability of the method even when extremely complicated bonding patterns are met. More importantly, however, it highlights the necessity to use special strategies when different points of a potential energy surface are investigated, e.g., during chemical reactions.

An ab initio multireference study of reductive eliminations from organoferrates(iii) in the gas-phase: it is all about the spin state

The reductive elimination reaction from organoferrates(iii) of the composition [FeR₃R′]⁻ is studied by state-of-the-art multireference electronic structure calculations.

Electronic Communication as a Transferable Property of Molecular Bridges?

Electronic communication through molecular bridges is important for different types of experiments, such as single-molecule conductance, electron transfer, superexchange spin coupling, and intramolecular singlet fission. In many instances, the chemical structure of the bridge determines how the two parts it is connecting communicate, and does so in ways that are transferable between these different manifestations (for example, high conductance often correlates with strong antiferromagnetic spin coupling, and low conductance due to destructive quantum interference correlates with ferromagnetic coupling). Defining electronic communication as a transferable property of the bridge can help transfer knowledge between these different areas of research. Examples and limits of such transferability are discussed here, along with some possible directions for future research, such as employing spin-coupled and mixed-valence systems as structurally well-controlled proxies for understanding molecular conductance and for validating first-principles theoretical methodologies, building conceptual understanding for the growing experimental work on intramolecular singlet fission, and developing measures for the transferability of electronic communication as a bridge property.

Orbital Entanglement Analysis of Exchange-Coupled Systems

A new tool for the interpretation of multiconfigurational wave functions representing the spin states of exchange-coupled transition metal complexes is introduced. Based on orbital entanglement measures, herein derived from multiconfigurational density matrix renormalization group calculations, the complexity of the wave function is reduced, thus facilitating a connection with established concepts for the interpretation of magnetically coupled systems. We show that the entanglement of localized orbitals with a small basis set is a good representation of the magnetic coupling topology and that it is sensitive to chemical changes in homologous complexes. Furthermore, we introduce a measure for the magnetic relevance of orbitals in the active subspace and a concept for the quantitative comparison of different chemical species. The approach presented here will be easily applicable to higher nuclearity clusters, providing a direct insight into all states of the Heisenberg spin ladder for systems previously accessible only by single-configurational methods.

Carotenoid Nuclear Reorganization and Interplay of Bright and Dark Excited States

We report quantum chemical calculations using multireference perturbation theory (MRPT) with the density matrix renormalization group (DMRG) plus photothermal deflection spectroscopy measurements to investigate the manifold of carotenoid excited states and establish their energies relative to the bright state (S₂) as a function of nuclear reorganization. We conclude that the primary photophysics and function of carotenoids are determined by interplay of only the bright (S₂) and lowest-energy dark (S₁) states. The lowest-lying dark state, far from being energetically distinguishable from the lowest-lying bright state along the entire excited-state nuclear reorganization pathway, is instead computed to be either the second or first excited state depending on what equilibrium geometry is considered. This result suggests that, rather than there being a dark intermediate excited state bridging a non-negligible energy gap from the lowest-lying dark state to the lowest-lying bright state, there is in fact no appreciable energy gap to bridge following photoexcitation. Instead, excited-state nuclear reorganization constitutes the bridge from S₂ to S₁, in the sense that these two states attain energetic degeneracy along this pathway.

Implications of structural heterogeneity for the electronic structure of the final oxygen-evolving intermediate in photosystem II

Heterogeneity in intermediate catalytic states of the oxygen-evolving complex (OEC) of Photosystem II is known from a wide range of experimental and theoretical data, but its potential implications for the mechanism of water oxidation remain unexplored. We delineate the consequences of structural heterogeneity for the final step of the catalytic cycle by tracing the evolution of three spectroscopically relevant and structurally distinct components of the last metastable S₃ state to the transient O₂-evolving S₄ state of the OEC. Using quantum chemical calculations, we show that each S₃ isomer leads to a different electronic structure formulation for the active S₄ state. Crucially, in addition to previously hypothesized Mn(IV)-oxyl species, we establish for the first time, how a genuine Mn(V)-oxo can be obtained in the catalytically active S₄ state: this takes the form of a five-coordinate and locally high-spin (S_Mn = 1) Mn(V) site. This formulation for the S₄ state evolves naturally from a preceding S₃-state structural intermediate that contains a quasi-trigonal-bipyramidal Mn(IV) ion. The results strongly suggest that water binding in the S₃ state is not prerequisite for reaching the oxygen-evolving S₄ state of the complex, supporting the notion that both substrates are preloaded at the beginning of the catalytic cycle. This scenario allows true four-electron metal-centered hole accumulation to precede OO bond formation and hence the latter can proceed via a genuine even-electron mechanism. This can occur as intramolecular nucleophilic coupling of two oxo units synchronously with the binding of a water substrate for the next catalytic cycle.

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.

The S3 State of the Oxygen-Evolving Complex: Overview of Spectroscopy and XFEL Crystallography with a Critical Evaluation of Early-Onset Models for O–O Bond Formation

The catalytic cycle of the oxygen-evolving complex (OEC) of photosystem II (PSII) comprises five intermediate states Si (i = 0–4), from the most reduced S0 state to the most oxidized S4, which spontaneously evolves dioxygen. The precise geometric and electronic structure of the Si states, and hence the mechanism of O–O bond formation in the OEC, remain under investigation, particularly for the final steps of the catalytic cycle. Recent advances in protein crystallography based on X-ray free-electron lasers (XFELs) have produced new structural models for the S3 state, which indicate that two of the oxygen atoms of the inorganic Mn4CaO6 core of the OEC are in very close proximity. This has been interpreted as possible evidence for “early-onset” O–O bond formation in the S3 state, as opposed to the more widely accepted view that the O–O bond is formed in the final state of the cycle, S4. Peroxo or superoxo formation in S3 has received partial support from computational studies. Here, a brief overview is provided of spectroscopic information, recent crystallographic results, and computational models for the S3 state. Emphasis is placed on computational S3 models that involve O–O formation, which are discussed with respect to their agreement with structural information, experimental evidence from various spectroscopic studies, and substrate exchange kinetics. Despite seemingly better agreement with some of the available crystallographic interpretations for the S3 state, models that implicate early-onset O–O bond formation are hard to reconcile with the complete line of experimental evidence, especially with X-ray absorption, X-ray emission, and magnetic resonance spectroscopic observations. Specifically with respect to quantum chemical studies, the inconclusive energetics for the possible isoforms of S3 is an acute problem that is probably beyond the capabilities of standard density functional theory.

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.

Spin, Valence, and Structural Isomerism in the S3 State of the Oxygen-Evolving Complex of Photosystem II as a Manifestation of Multimetallic Cooperativity

Photosynthetic water oxidation is catalyzed by a Mn₄CaO₅-cluster in photosystem II through an S-state cycle. Understanding the roles of heterogeneity in each S-state, as identified recently by the EPR spectroscopy, is very important to gain a complete description of the catalytic mechanism. We performed herein hybrid DFT calculations within the broken-symmetry formalism and associated analyses of Heisenberg spin models to study the electronic and spin structures of various isomeric structural motifs (hydroxo-oxo, oxyl-oxo, peroxo, and superoxo species) in the S₃ state. Our extensive study reveals several factors that affect the spin ground state: (1) (formal) Mn oxidation state; (2) metal-ligand covalency; (3) coordination geometry; and (4) structural change of the Mn cluster induced by alternations in Mn···Mn distances. Some combination of these effects could selectively stabilize/destabilize some spin states. We found that the high spin state ( S_total = 6) of the oxyl-oxo species can be causative for catalytic function, which manifests through mixing of the metal-ligand character in magnetic orbitals at relatively short O5···O6 distances (<2.0 Å) and long Mn_A···O5 distances (>2.0 Å). These results will serve as a basis to conceptually identify and rationalize the physicochemical synergisms that can be evoked by the unique "distorted chair" topology of the cluster through cooperative Jahn-Teller effects on multimetallic centers.

Benchmarking quantum chemistry methods for spin-state energetics of iron complexes against quantitative experimental data

The accuracy of relative spin-state energetics predicted by selected quantum chemistry methods: coupled cluster theory at the CCSD(T) level, multiconfigurational perturbation theory (CASPT2, NEVPT2), multireference configuration interaction at the MRCISD+Q level, and a number of DFT methods, is quantitatively evaluated by comparison with the experimental data of four octahedral iron complexes. The available experimental data, either spin-forbidden transition energies or spin crossover enthalpies, are corrected for relevant environmental effects in order to derive the quantitative benchmark set of iron spin-state energetics. Comparison of theory predictions with the resulting reference data: (1) validates the high accuracy of the CCSD(T) method, particularly when based on Kohn-Sham orbitals, giving the maximum error below 2 kcal mol-1 and the mean absolute error (MAE) below 1 kcal mol-1; (2) corroborates the tendency of CASPT2 to systematically overstabilize higher-spin states by up to 5.5 kcal mol-1; (3) confirms that the latter problem is partly remedied by the recently proposed CASPT2/CC approach [Phung et al., J. Chem. Theory Comput., 2018, 14, 2446-2455]; (4) demonstrates that NEVPT2 performs worse than CASPT2, by giving errors up to 7 kcal mol-1; (5) shows that the accuracy of MRCISD+Q spin-state energetics strongly depends on the size-consistency correction: the Davidson-Silver and Pople corrections perform best (MAE < 3 kcal mol-1), whereas the standard Davidson correction is not recommended (MAE of 7 kcal mol-1). Only a few DFT methods (including the best performing ones identified in this study: B2PLYP-D3 and OPBE) are able to provide a balanced description of the spin-state energetics for all four studied iron complexes simultaneously, corroborating the non-universality problem of approximate density functionals.

Meeting the Challenge of Magnetic Coupling in a Triply-Bridged Chromium Dimer: Complementary Broken-Symmetry Density Functional Theory and Multireference Density Matrix Renormalization Group Perspectives

Face-sharing octahedral dinuclear Cr(III) compounds with d³–d³ electronic configurations represent nontrivial examples of electronic complexity, posing particular challenges for theoretical and computational studies. A tris-hydroxy-bridged Cr(III)–Cr(III) system has proven to be a richly rewarding target for studies of magnetism and electron paramagnetic resonance spectroscopy. It was also reported to be a peculiarly difficult system to treat with density functional theory (DFT). In this work the magnetic coupling problem for this dimer is approached with broken-symmetry (BS)-DFT and multireference calculations that utilize the density matrix renormalization group (DMRG) to handle full-valence active spaces. BS-DFT is shown to recover the correct ordering and energy spacing of Heisenberg spin states if used in conjunction with appropriate spin projection procedures, albeit with pronounced functional sensitivity. The contrasting conclusions of previous studies are traced to incorrect inclusion of electronically excited configurations. Analysis of the direct and differential overlap of corresponding orbital pairs from the BS-DFT solution indicates that metal–metal through-space interaction is the dominant contributor to antiferromagnetic coupling. At the DFT level a procedure that utilizes pseudopotential substitution is demonstrated that allows evaluation of the direct exchange vs superexchange contributions. A complementary description is obtained with DMRG-SCF calculations that enable state-averaged CASSCF calculations with both metal and bridge orbitals in the active space. A localized orbital subspace analysis supports the DFT conclusions that in contrast to doubly bridged isoelectronic analogues, antiferromagnetic coupling in the chromium dimer arises primarily from direct metal–metal interaction but is significantly enhanced by ligand-mediated superexchange.

UNO DMRG CAS CI calculations of binuclear manganese complex Mn(IV)2 O2 (NHCHCO2 )4 : Scope and applicability of Heisenberg model

Both direct exchange and super-exchange interactions cooperate to realize inter-spin magnetic interaction in binuclear manganese complex Mn(IV)₂ O₂ (NHCHCO₂ )₄ with a di-μ-oxo path. We revisited this spin system using DMRG CAS methods and CAS selection procedures. Our results indicate that our previous "dynamically extended spin polarization" (DE-SP) procedure for organic polyradicals and so forth does not work well. Thus, we have examined another selection procedure, the "dynamically extended super-exchange" (DE-SE) procedure. DMRG CASCI [18,18] by UB3LYP(HS)-UNO(DE-SE) can realize antiferromagnetic J values similar to experimental ones (-87 cm^-1 ). In addition, all J values between all spin states (HS[septet],IS[quintet],IS[triplet],LS[singlet])were also shown to be correct under sufficiently large M values. © 2018 Wiley Periodicals, Inc.