Magnetic Coupling in Transition-Metal Binuclear Complexes by Spin-Flip Time-Dependent Density Functional Theory

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.

Natural orbitals and two-particle correlators as tools for the analysis of effective exchange couplings in solids

Using generalizations of spin-averaged natural orbitals and two-particle charge correlators for solids, we investigate the electronic structure of antiferromagnetic transition-metal oxides with a fully self-consistent, imaginary-time GW method. Our findings disagree with the Goodenough-Kanamori (GK) rules that are commonly used for the qualitative interpretation of such solids. First, we found a strong dependence of the natural orbital occupancies on momenta, contradicting GK assumptions. Second, along the momentum path, the character of natural orbitals changes. In particular, the contributions of oxygen 2s orbitals are important, which has not been considered in the GK rules. To analyze the influence of the electronic correlation on the values of effective exchange coupling constants, we use both natural orbitals and two-particle correlators and show that electronic screening modulates the degree of superexchange by stabilizing the charge-transfer contributions, which greatly affects these coupling constants. Finally, we give a set of predictions and recommendations regarding the use of density functional, Green's function, and wave-function methods for evaluating effective magnetic couplings in molecules and solids.

Revisiting magnetic exchange couplings in heterodinuclear complexes through the decomposition method in KS-DFT

Providing tools to understand the physical mechanisms governing magnetic properties in transition metal-based compounds is still of great interest. Here, the magnetic exchange coupling in a series of heterodinuclear complexes is investigated by means of the decomposition method. This work presents the first application of the decomposition method to systems where magnetic centres may bear more than one unpaired electron. By decomposing the coupling into three physical contributions (direct exchange, kinetic exchange, and spin polarisation), we provide numerical arguments to confirm or infirm the rationalisation allowed by the conceptual analysis of the magnetic d orbitals. We also take advantage of the recently proposed generalisation of the method [David et al., J. Chem. Theory Comput., 2023, 19, 157] to get more insights into the underlying mechanisms by disentangling the coupling between centres into its electron-electron interactions.

Magnetic exchange interactions in binuclear and tetranuclear iron(III) complexes described by spin-flip DFT and Heisenberg effective Hamiltonians

Low-energy spectra of single-molecule magnets (SMMs) are often described by Heisenberg Hamiltonians. Within this formalism, exchange interactions between magnetic centers determine the ground-state multiplicity and energy separation between the ground and excited states. In this contribution, we extract exchange coupling constants (J) for a set of iron (III) binuclear and tetranuclear complexes from all-electron calculations using non-collinear spin-flip time-dependent density functional theory (NC-SF-TDDFT). For 12 binuclear complexes with J-values ranging from -6 to -132 cm^-1 , our benchmark calculations using the short-range hybrid ωPBEh functional and 6-31G(d,p) basis set agree well with the experimentally derived values (mean absolute error of 4.7 cm^-1 ). For the tetranuclear SMMs, the computed J constants are within 6 cm^-1 from the experimentally derived values. We explore the range of applicability of the Heisenberg model by analyzing bonding patterns in these Fe(III) complexes using natural orbitals (NO), their occupations, and the number of effectively unpaired electrons. The results illustrate the efficiency of the spin-flip protocol for computing the exchange couplings and the utility of the NO analysis in assessing the validity of effective spin Hamiltonians.

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.

Validation of the Green's Function Approximation for the Calculation of Magnetic Exchange Couplings

In this work, we assess the potential of the Green's function approximation to predict isotropic magnetic exchange couplings and to reproduce the standard broken-symmetry energy difference approach for transition metal complexes. To this end, we have selected a variety of heterodinuclear, homodinuclear, and polynuclear systems containing 3d transition metal centers and computed the couplings using both the Green's function and energy difference methods. The Green's function approach is shown to have mixed results for the cases tested. For dinuclear complexes with large strength couplings (≳50 cm^-1), the Green's function method is unable to reliably reproduce the energy difference values. However, for weaker dinuclear couplings, the Green's function approach acceptably reproduces broken-symmetry energy difference couplings. In polynuclear cases, the Green's function approximation worked remarkably well, especially for Fe^III complexes. On the other hand, for a Ni^II polynuclear complex, qualitatively wrong couplings are predicted. Overall, the evaluation of exchange couplings from local rigid magnetization rotations offers a powerful alternative to time-consuming energy differences methods for large polynuclear transition metal complexes, but to achieve a quantitative agreement, some improvements to the method are needed.

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.

Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory (MRSF-TDDFT) as a Simple yet Accurate Method for Diradicals and Diradicaloids

Due to their multiconfigurational nature featuring strong electron correlation, accurate description of diradicals and diradicaloids is a challenge for quantum chemical methods. The recently developed mixed-reference spin-flip (MRSF)-TDDFT method is capable of describing the multiconfigurational electronic states of these systems while avoiding the spin-contamination pitfalls of SF-TDDFT. Here, we apply MRSF-TDDFT to study the adiabatic singlet-triplet (ST) gaps in a series of well-known diradicals and diradicaloids. On average, MRSF displays a very high prediction accuracy of the adiabatic ST gaps with the mean absolute error (MAE) amounting to 0.14 eV. In addition, MRSF is capable of accurately describing the effect of the Jahn-Teller distortion occurring in the trimethylenemethane diradical, the violation of the Hund rule in a series of the didehydrotoluene diradicals, and the potential energy surfaces of the didehydrobenzene (benzyne) diradicals. A convenient criterion for distinguishing diradicals and diradicaloids is suggested on the basis of the easily obtainable quantities. In all of these cases, which are difficult for the conventional methods of density functional theory (DFT), MRSF shows results consistent with the experiment and the high-level ab initio computations. Hence, the present study documents the reliability and accuracy of MRSF and lays out the guidelines for its application to strongly correlated molecular systems.

Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory

Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe²⁺, Fe³⁺, and Co³⁺, including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. We find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method-in terms of computational cost and in terms of accuracy on monatomic systems-and therefore we chose this method for molecular benchmarks. Then Kohn-Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF₂, FeCl₂, and CoF₂. We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.

Spin-flip methods in quantum chemistry

This Perspective discusses salient features of the spin-flip approach to strong correlation and describes different methods that sprung from this idea. The spin-flip treatment exploits the different physics of low-spin and high-spin states and is based on the observation that correlation is small for same-spin electrons. By using a well-behaved high-spin state as a reference, one can access problematic low-spin states by deploying the same formal tools as in the excited-state treatments (i.e., linear response, propagator, or equation-of-motion theories). The Perspective reviews applications of this strategy within wave function and density functional theory frameworks as well as the extensions for molecular properties and spectroscopy. The utility of spin-flip methods is illustrated by examples. Limitations and proposed future directions are also discussed.

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.

Quantitative Assessment of rPM6 for Fluorine- and Chlorine-Containing Metal Complexes: Comparison with Experimental, First-Principles, and Other Semiempirical Results

We report a reparameterization of PM6 parameters for fluorine and chlorine using our training set containing transition metal complexes. Spin unrestricted calculations with the resulting rPM6 (UrPM6) were examined quantitatively using two test sets: (i) the description of magnetic interactions in 25 dinuclear metal complexes and (ii) the prediction of barrier heights and reaction energies for epoxidation and fluorination reactions catalyzed by high-valent manganese-oxo species. The conventional UPM6 and UPM7 methods were also evaluated for comparison on the basis of either experimental or computational (the UB3LYP/SVP level) outcomes. The merits of UrPM6 are highlighted by both the test sets. As regards magnetic exchange coupling constants, the UrPM6 method had the smallest mean absolute errors from the experimental data (19 cm-1), followed by UPM7 (119 cm-1) and UPM6 (373 cm-1). For the epoxidation and fluorination reactions, all of the transition state searches were successful using UrPM6, while the success rates obtained by UPM6 and UPM7 were only 50%. The UrPM6-optimized stationary points also agreed well with the reference UB3LYP-optimized geometries. The accuracy for estimating reaction energies was also greatly remedied.

Fermi-Löwdin orbital self-interaction correction to magnetic exchange couplings

We analyze the effect of removing self-interaction error on magnetic exchange couplings using the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method in the framework of density functional theory (DFT). We compare magnetic exchange couplings obtained from self-interaction-free FLOSIC calculations with the local spin density approximation (LSDA) with several widely used DFT realizations and wave function based methods. To this end, we employ the linear H-He-H model system, six organic radical molecules, and [Cu₂Cl₆]^2- as representatives of different types of magnetic interactions. We show that the simple self-interaction-free version of LSDA improves calculated couplings with respect to LSDA in all cases, even though the nature of the exchange interaction varies across the test set, and in most cases, it yields results comparable to modern hybrids and range-separated approximate functionals.

Singlet-triplet energy gaps and the degree of diradical character in binuclear copper molecular magnets characterized by spin-flip density functional theory

Molecular magnets, defined here as organic polyradicals, can be used as building blocks in the fabrication of novel and structurally diverse magnetic light-weight materials. We present a theoretical investigation of the lowest spin states of several binuclear copper diradicals. In contrast to previous studies, we consider not only the energetics of the low-lying states (which are related to the exchange-coupling parameter within the Heisenberg-Dirac-van-Vleck model), but also the character of the diradical states themselves. We use natural orbitals, their occupations, and the number of effectively unpaired electrons to quantify bonding patterns in these systems. We compare the performance of spin-flip time-dependent density functional theory (SF-TDDFT) using various functionals and effective core potentials against the wave function based approach, equation-of-motion spin-flip coupled-cluster method with single and double substitutions (EOM-SF-CCSD). We find that SF-TDDFT paired with the PBE50 and B5050LYP functionals performs comparably to EOM-SF-CCSD, with respect to both singlet-triplet gaps and states' characters. Visualization of frontier natural orbitals shows that the unpaired electrons are localized on copper centers, in some cases exhibiting slight through-bond interaction via copper d-orbitals and p-orbitals of neighboring ligand atoms. The analysis reveals considerable interactions between the formally unpaired electrons in the antiferromagnetic diradicaloids, meaning that they are poorly described by the Heisenberg-Dirac-van-Vleck model. Thus, for these systems the experimentally derived exchange-coupling parameters are not directly comparable with the singlet-triplet gaps. This explains systematic discrepancies between the computed singlet-triplet energy gaps and the exchange-coupling parameters extracted from experiment.

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

The accurate description of magnetic level energetics in oligonuclear exchange-coupled transition-metal complexes remains a formidable challenge for quantum chemistry. The density matrix renormalization group (DMRG) brings such systems for the first time easily within reach of multireference wave function methods by enabling the use of unprecedentedly large active spaces. But does this guarantee systematic improvement in predictive ability and, if so, under which conditions? We identify operational parameters in the use of DMRG using as a test system an experimentally characterized mixed-valence bis-μ-oxo/μ-acetato Mn(III,IV) dimer, a model for the oxygen-evolving complex of photosystem II. A complete active space of all metal 3d and bridge 2p orbitals proved to be the smallest meaningful starting point; this is readily accessible with DMRG and greatly improves on the unrealistic metal-only configuration interaction or complete active space self-consistent field (CASSCF) values. Orbital optimization is critical for stabilizing the antiferromagnetic state, while a state-averaged approach over all spin states involved is required to avoid artificial deviations from isotropic behavior that are associated with state-specific calculations. Selective inclusion of localized orbital subspaces enables probing the relative contributions of different ligands and distinct superexchange pathways. Overall, however, full-valence DMRG-CASSCF calculations fall short of providing a quantitative description of the exchange coupling owing to insufficient recovery of dynamic correlation. Quantitatively accurate results can be achieved through a DMRG implementation of second order N-electron valence perturbation theory (NEVPT2) in conjunction with a full-valence metal and ligand active space. Perspectives for future applications of DMRG-CASSCF/NEVPT2 to exchange coupling in oligonuclear clusters are discussed.

Helical Folding-Induced Stabilization of Ferromagnetic Polyradicals Based on Triarylmethyl Radical Derivatives

Magnetic ordering in purely organic π-conjugated materials is a challenging, rare, and desirable event. The interest lies on the unique magnetic properties derived from high-spin carbon-based polymers/macromolecules tailored through appropriate synthetic routes. Ground-breaking achievements have been reported regarding magnetic ordering in an organic polymer using spin clusters as building blocks. This strategy leads to two-dimensional extended polyradicals with a concomitant loss of appealing macroscopic properties such as expected magnetic anisotropy in elongated shaped macromolecules containing carbon-bearing radicals. Here we provide compelling evidence of a secondary structure-induced stabilization of ferromagnetic polyradicals with robust magnetic properties and strongly suggest revisiting a discarded attempt to obtain polymeric linear-like radicals. An alternative synthetic approach is also proposed, based on polyradicals obtained from discrete molecular precursors (oligomers) long enough to ensure a secondary structure, rather than from polymerization processes.

Valence excitation energies of alkenes, carbonyl compounds, and azabenzenes by time-dependent density functional theory: linear response of the ground state compared to collinear and noncollinear spin-flip TDDFT with the Tamm-Dancoff approximation

Time-dependent density functional theory (TDDFT) holds great promise for studying photochemistry because of its affordable cost for large systems and for repeated calculations as required for direct dynamics. The chief obstacle is uncertain accuracy. There have been many validation studies, but there are also many formulations, and there have been few studies where several formulations were applied systematically to the same problems. Another issue, when TDDFT is applied with only a single exchange-correlation functional, is that errors in the functional may mask successes or failures of the formulation. Here, to try to sort out some of the issues, we apply eight formulations of adiabatic TDDFT to the first valence excitations of ten molecules with 18 density functionals of diverse types. The formulations examined are linear response from the ground state (LR-TDDFT), linear response from the ground state with the Tamm-Dancoff approximation (TDDFT-TDA), the original collinear spin-flip approximation with the Tamm-Dancoff (TD) approximation (SF1-TDDFT-TDA), the original noncollinear spin-flip approximation with the TDA approximation (SF1-NC-TDDFT-TDA), combined self-consistent-field (SCF) and collinear spin-flip calculations in the original spin-projected form (SF2-TDDFT-TDA) or non-spin-projected (NSF2-TDDFT-TDA), and combined SCF and noncollinear spin-flip calculations (SF2-NC-TDDFT-TDA and NSF2-NC-TDDFT-TDA). Comparing LR-TDDFT to TDDFT-TDA, we observed that the excitation energy is raised by the TDA; this brings the excitation energies underestimated by full linear response closer to experiment, but sometimes it makes the results worse. For ethylene and butadiene, the excitation energies are underestimated by LR-TDDFT, and the error becomes smaller making the TDA. Neither SF1-TDDFT-TDA nor SF2-TDDFT-TDA provides a lower mean unsigned error than LR-TDDFT or TDDFT-TDA. The comparison between collinear and noncollinear kernels shows that the noncollinear kernel drastically reduces the spin contamination in the systems considered here, and it makes the results more accurate than collinear spin-flip TDDFT for functionals with a low percentage of Hartree-Fock exchange and sometimes for functionals with a higher percentage of Hartree-Fock exchange, but it yields less accurate results than ground-state TDDFT.