Assessing the Calculation of Exchange Coupling Constants and Spin Crossover Gaps Using the Approximate Projection Model To Improve Density Functional Calculations

Electronic Structure and Magneto-Structural Correlations Study of Cu2UL Trinuclear Schiff Base Complexes: A 3d-5f-3d Case

The magnetic properties of trinuclear Schiff base complexes M₂AnLⁱ (M^II = Zn, Cu; An^IV = Th, U; Lⁱ = Schiff base; i = 1-4, 6, 7, 9), exhibiting the [M(μ-O)₂]₂U core structure with adjacent M1···U and M2···U and next-adjacent M1···M2 interactions, featuring 3d-5f-3d subsystems, have been investigated theoretically using relativistic ZORA/B3LYP computations combined with the broken symmetry (BS) approach. Bond order and natural population analyses reveal that the covalent contribution to the bonding within the Cu-O-U coordination is important thus favoring superexchange coupling between the transition metal and the uranium magnetic centers. The calculated coupling constants J_CuU between the Cu and U atoms, agree with the observed shift from the antiferromagnetic (AF) character of the L^1,2,3,4 complexes to the ferromagnetic (ferro) of the L^6,7,9 ones. The structural parameters, i.e., the Cu···U distances and the Cu-O-U angles, as well as the electronic factors driving the magnetic couplings are discussed. The analyses are supported by the study of the mixed ZnCuULⁱ and Cu₂ThLⁱ systems, where in the first complex the Cu^II (3d⁹) ion is replaced by the diamagnetic Zn^II (3d¹⁰) one, whereas in the second complex the U^IV (5f²) paramagnetic center is replaced by the diamagnetic Th^IV (5f⁰) one.

Good Vibrations: Calculating Excited-State Frequencies Using Ground-State Self-Consistent Field Models

The use of Δ-self-consistent field (SCF) approaches for studying excited electronic states has received a renewed interest in recent years. In this work, the use of this scheme for calculating excited-state vibrational frequencies is examined. Results from Δ-SCF calculations for a set of representative molecules are compared with those obtained using configuration interaction with single substitutions (CIS) and time-dependent density functional theory (TD-DFT) methods. The use of an approximate spin purification model is also considered for cases where the excited-state SCF solution is spin-contaminated. The results of this work demonstrate that an SCF-based description of an excited-state potential energy surface can be an accurate and cost-effective alternative to CIS and TD-DFT methods.

Cooperative redox and spin activity from three redox congeners of sulfur-bridged iron nitrosyl and nickel dithiolene complexes

The synthesis of sulfur-bridged Fe-Ni heterobimetallics was inspired by Nature's strategies to "trick" abundant first row transition metals into enabling 2-electron processes: redox-active ligands (including pendant iron-sulfur clusters) and proximal metals. Our design to have redox-active ligands on each metal, NO on iron and dithiolene on nickel, resulted in the observation of unexpectedly intricate physical properties. The metallodithiolate, (NO)Fe(N2S2), reacts with a labile ligand derivative of [NiII(S2C2Ph2)]0, NiDT, yielding the expected S-bridged neutral adduct, FeNi, containing a doublet {Fe(NO)}7. Good reversibility of two redox events of FeNi led to isolation of reduced and oxidized congeners. Characterization by various spectroscopies and single-crystal X-ray diffraction concluded that reduction of the FeNi parent yielded [FeNi]-, a rare example of a high-spin {Fe(NO)}8, described as linear FeII(NO-). Mössbauer data is diagnostic for the redox change at the {Fe(NO)}7/8 site. Oxidation of FeNi generated the 2[FeNi]+⇌[Fe2Ni2]2+ equilibrium in solution; crystallization yields only the [Fe2Ni2]2+ dimer, isolated as PF6- and BArF- salts. The monomer is a spin-coupled diradical between {Fe(NO)}7 and NiDT+, while dimerization couples the two NiDT+ via a Ni2S2 rhomb. Magnetic susceptibility studies on the dimer found a singlet ground state with a thermally accessible triplet excited state responsible for the magnetism at 300 K (χMT = 0.67 emu·K·mol-1, µeff = 2.31 µB), and detectable by parallel-mode EPR spectroscopy at 20 to 50 K. A theoretical model built on an H4 chain explains this unexpected low energy triplet state arising from a combination of anti- and ferromagnetic coupling of a four-radical molecular conglomerate.

Comparison of Linear Response Theory, Projected Initial Maximum Overlap Method, and Molecular Dynamics-Based Vibronic Spectra: The Case of Methylene Blue

The simulation of optical spectra is essential to molecular characterization and, in many cases, critical for interpreting experimental spectra. The most common method for simulating vibronic absorption spectra relies on the geometry optimization and computation of normal modes for ground and excited electronic states. In this report, we show that the utilization of such a procedure within an adiabatic linear response (LR) theory framework may lead to state mixings and a breakdown of the Born-Oppenheimer approximation, resulting in a poor description of absorption spectra. In contrast, computing excited states via a self-consistent field method in conjunction with a maximum overlap model produces states that are not subject to such mixings. We show that this latter method produces vibronic spectra much more aligned with vertical gradient and molecular dynamics (MD) trajectory-based approaches. For the methylene blue chromophore, we compare vibronic absorption spectra computed with the following: an adiabatic Hessian approach with LR theory-optimized structures and normal modes, a vertical gradient procedure, the Hessian and normal modes of maximum overlap method-optimized structures, and excitation energy time-correlation functions generated from an MD trajectory. Because of mixing between the bright S₁ and dark S₂ surfaces near the S₁ minimum, computing the adiabatic Hessian with LR theory and time-dependent density functional theory with the B3LYP density functional predicts a large vibronic shoulder for the absorption spectrum that is not present for any of the other methods. Spectral densities are analyzed and we compare the behavior of the key normal mode that in LR theory strongly couples to the optical excitation while showing S₁/S₂ state mixings. Overall, our study provides a note of caution in computing vibronic spectra using the excited-state adiabatic Hessian of LR theory-optimized structures and also showcases three alternatives that are less sensitive to adiabatic state mixing effects.

A computational mechanistic study of CH hydroxylation with mononuclear copper–oxygen complexes

A computational study of methane hydroxylation by oxygen-bound monocopper complexes.

Revealing the nature of electron correlation in transition metal complexes with symmetry breaking and chemical intuition

In this work, we provide a nuanced view of electron correlation in the context of transition metal complexes, reconciling computational characterization via spin and spatial symmetry breaking in single-reference methods with qualitative concepts from ligand-field and molecular orbital theories. These insights provide the tools to reliably diagnose the multi-reference character, and our analysis reveals that while strong (i.e., static) correlation can be found in linear molecules (e.g., diatomics) and weakly bound and antiferromagnetically coupled (monometal-noninnocent ligand or multi-metal) complexes, it is rarely found in the ground-states of mono-transition-metal complexes. This leads to a picture of static correlation that is no more complex for transition metals than it is, e.g., for organic biradicaloids. In contrast, the ability of organometallic species to form more complex interactions, involving both ligand-to-metal σ-donation and metal-to-ligand π-backdonation, places a larger burden on a theory's treatment of dynamic correlation. We hypothesize that chemical bonds in which inter-electron pair correlation is non-negligible cannot be adequately described by theories using MP2 correlation energies and indeed find large errors vs experiment for carbonyl-dissociation energies from double-hybrid density functionals. A theory's description of dynamic correlation (and to a less important extent, delocalization error), which affects relative spin-state energetics and thus spin symmetry breaking, is found to govern the efficacy of its use to diagnose static correlation.