Improved evaluation of spin-polarization energy contributions using broken-symmetry calculations

On Entangled Singlet Pure Diradicals

This work addresses a class of conjugated hydrocarbons that are expected to be singlet diradicals according to the topological Hückel Hamiltonian while possibly satisfying full on-bond electron pairing. These systems possess two degenerate singly occupied molecular orbitals (SOMOs), but aromaticity brought by properly positioned six-membered rings does prevent Jahn-Teller distortions. Density functional theory (DFT) calculations performed on two emblematic examples confirm the strong bond-length alternation in the closed-shell solutions and the clear spatial symmetry in the open-shell spin-unrestricted determinants, the latter solution always being found to have significantly lower energy. Since the SOMOs are here of different symmetry, the wave function is free from ionic valence-bond component, and spin decontamination of the unrestricted DFT solutions and wave function calculations at the CASSCF-plus-second-order-perturbation level confirm the expected pure diradical character of such molecules. In contrast to disjoint diradicals, the SOMOs of present systems have large amplitudes on neighbor atoms, and we propose to name them entangled pure diradicals, further providing some prescription rules for their design. Additional calculations point out the qualitative contrast between these molecules and the related diradicaloids.

Ultrafast Computational Screening of Molecules with Inverted Singlet-Triplet Energy Gaps Using the Pariser-Parr-Pople Semiempirical Quantum Chemistry Method

Molecules with an inverted energy gap between their first singlet and triplet excited states have promising applications in the next generation of organic light-emitting diode (OLED) materials. Unfortunately, such molecules are rare, and only a handful of examples are currently known. High-throughput virtual screening could assist in finding novel classes of these molecules, but current efforts are hampered by the high computational cost of the required quantum chemical methods. We present a method based on the semiempirical Pariser-Parr-Pople theory augmented by perturbation theory and show that it reproduces inverted gaps at a fraction of the cost of currently employed excited-state calculations. Our study paves the way for ultrahigh-throughput virtual screening and inverse design to accelerate the discovery and development of this new generation of OLED materials.

What governs magnetic exchange couplings in radical-bridged dinuclear complexes?

Coupling transition metal or lanthanide ions through a radical bridging ligand is a promising route to increase performances in the area of single molecular magnets. A better understanding of the underlying physical mechanisms governing the magnetic exchange couplings is thus of valuable importance to design future compounds. Here, couplings in three series of metal-radical-metal compounds based on transition metal ions are investigated by means of the decomposition/recomposition methods. This work presents the generalisation and first application of the method to systems with an arbitrary number of magnetic centres featuring several unpaired electrons. Thanks to the decomposition into the three main contributions (direct exchange, kinetic exchange, and spin polarisation) as well as a description in terms of electron-electron interactions, we study the influence of the nature of the metal centre and the radical ligand on the couplings. We combine the energetic contributions extracted with orbital and charge population analysis to rationalise the results.

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.

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.

Difficulty of the evaluation of the barrier height of an open-shell transition state between closed shell minima: The case of small C4n rings

C_4n cyclacenes exhibit strong bond-alternation in their equilibrium geometry. In the two equivalent geometries, the system keeps an essentially closed-shell character. The two energy minima are separated by a transition state suppressing the bond-alternation, where the wave function is strongly diradical. This paper discusses the physical factors involved in this energy difference and possible evaluations of the barrier height. The barrier given as the energy difference between the restricted density functional theory (DFT)/B3LYP for the equilibrium and the broken symmetry DFT/B3LYP of the transition state is either negative or small, in contradiction with the most reliable Wave Function Theory calculations. The minimal (two electrons in two molecular orbitals) Complete Active Space self-consistent field (CASSCF) overestimates the barrier, and the subsequent second-order perturbation cancels it. Due to the collective character of the spin-polarization effect, it is necessary to perform a full π CASSCF + second-order perturbation to reach a reasonable value of the barrier, but this type of treatment cannot be applied to large molecules. DFT procedures treating on an equal foot the closed-shell and open-shell geometries have been explored, such as Mixed-Reference Spin-Flip Time-dependent-DFT and a new spin-decontamination proposal, namely, DFT-dressed configuration interaction, but the results still depend on the density functional. M06-2X without or with spin-decontamination gives the best agreement with the accurate wave function results.

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.