A Convenient DFT-Based Strategy for Predicting Transition Temperatures of Valence Tautomeric Molecular Switches

Predicting valence tautomerism in diverse cobalt-dioxolene complexes: elucidation of the role of ligands and solvent

The ability of molecular switches to reversibly interconvert between different forms promises potential applications at the scale of single molecules up to bulk materials. One type of molecular switch comprises cobalt-dioxolene compounds that exhibit thermally-induced valence tautomerism (VT) interconversions between low spin Co(iii)-catecholate (LS-CoIII-cat) and high spin Co(ii)-semiquinonate (HS-CoII-sq) forms. Two families of these compounds have been investigated for decades but have generally been considered separately: neutral [Co(diox)(sq)(N2L)] and cationic [Co(diox)(N4L)]+ complexes (diox = generic dioxolene, N2L/N4L = bidentate/tetradentate N-donor ancillary ligand). Computational identification of promising new candidate compounds prior to experimental exploration is beneficial for environmental and cost considerations but requires a thorough understanding of the underlying thermochemical parameters that influence the switching. Herein, we report a robust approach for the analysis of both cobalt-dioxolene families, which involved a quantitative density functional theory-based study benchmarked with reliable quasi-experimental references. The best-performing M06L-D4/def2-TZVPP level of theory has subsequently been verified by the synthesis and experimental investigation of three new complexes, two of which exhibit thermally-induced VT, while the third remains in the LS-CoIII-cat form across all temperatures, in agreement with prediction. Valence tautomerism in solution is markedly solvent-dependent, but the origin of this has not been definitively established. We have extended our computational approach to elucidate the correlation of VT transition temperature with solvent stabilisation energy and change in dipole moment. This new understanding may inform the development of VT compounds for applications in soft materials including films, gels, and polymers.

Low-spin to low-spin valence tautomeric transition in cobalt bis-dioxolenes

Cobalt dioxolenes are a well-known class of switchable coordination compounds showing intramolecular electron transfer, which is always accompanied by a spin state change at the cobalt center. Here, we present the first example of thermally switchable cobalt bis-dioxolenes where intramolecular electron transfer seems to take place, but the spin state change is suppressed. This leads to the detection of thermal transition between a common ls-CoIII(SQ˙-)(Cat2-) and an extremely rare ls-CoII(SQ˙-)2 electronic state (hs - high-spin, ls - low-spin, SQ˙- - benzosemiquinonate(1-) radical and Cat2- - catecholate(2-)). Parallel to the present work, a similar work but on cobalt mono-dioxolenes has just appeared (Chem. Eur. J., 2023, 29, e202300091), suggesting thermal transition between ls-CoIII(Cat2-) and ls-CoII(SQ˙-) electronic states.

Metal-ligand interactions in a redox active ligand system. Electrochemistry and spectroscopy of [M(dipyvd)2]n+ (M=Zn, Ni, n=0, 1, 2)

Reaction of nickel and zinc triflates with the tridentate leucoverdazyl 1-isopropyl-3,5-di (2'-pyridyl)-6-oxo-2H-tetrazine (dipyvdH) and triethylamine resulted in the neutral coordination compounds M(dipyvd)2 (M = Ni,Zn). In acetonitrile, both compounds undergo two one electron oxidation processes, Zn (dipyvd)2 at -0.28 V and -0.12 V and Ni(dipyvd)2 at -0.32 V and -0.15 V vs ferrocene/ferricenium. Oxidations are ligand based resulting in an intermediate mixed valence species and a cationic bis(verdazyl) compound respectively. Oxidation of the ligand changes a localized, antiaromatic, non-planar 8π electron anion to a planar, delocalized 7π electron radical. The change in ligand structure results in an increase in the octahedral ligand field splitting from 10,500 cm-1 to ∼13,000 cm-1, suggesting an increase in the pi acceptor character of the ligand. In the mixed valence species, spectroscopic data suggests minimal interaction between ligands mediated by the metal center; i.e., these are class I-II systems in the Robin-Day classification.

Controlling Spin Crossover in a Family of Dinuclear Fe(III) Complexes via the Bis(catecholate) Bridging Ligand

Spin crossover (SCO) complexes can reversibly switch between low spin (LS) and high spin (HS) states, affording possible applications in sensing, displays, and molecular electronics. Dinuclear SCO complexes with access to [LS-LS], [LS-HS], and [HS-HS] states may offer increased levels of functionality. The nature of the SCO interconversion in dinuclear complexes is influenced by the local electronic environment. We report the synthesis and characterization of [{FeIII(tpa)}2spiro](PF6)2 (1), [{FeIII(tpa)}2Br4spiro](PF6)2 (2), and [{FeIII(tpa)}2thea](PF6)2 (3) (tpa = tris(2-pyridylmethyl)amine, spiroH4 = 3,3,3',3'-tetramethyl-1,1'-spirobi(indan)-5,5',6,6'-tetraol, Br4spiroH4 = 3,3,3',3'-tetramethyl-1,1'-spirobi(indan)-4,4',7,7'-tetrabromo-5,5',6,6'-tetraol, theaH4 = 2,3,6,7-tetrahydroxy-9,10-dimethyl-9,10-dihydro-9,10-ethanoanthracene), utilizing non-conjugated bis(catecholate) bridging ligands. In the solid state, magnetic and structural analysis shows that 1 remains in the [HS-HS] state, while 2 and 3 undergo a partial SCO interconversion upon cooling from room temperature involving the mixed [LS-HS] state. In solution, all complexes undergo SCO from [HS-HS] at room temperature, via [LS-HS] to mixtures including [LS-LS] at 77 K, with the extent of SCO increasing in the order 1 < 2 < 3. Gas phase density functional theory calculations suggest a [LS-LS] ground state for all complexes, with the [LS-HS] and [HS-HS] states successively destabilized. The relative energy separations indicate that ligand field strength increases following spiro4- < Br4spiro4- < thea4-, consistent with solid-state magnetic and EPR behavior. All three complexes show stabilization of the [LS-HS] state in relation to the midpoint energy between [LS-LS] and [HS-HS]. The relative stability of the [LS-HS] state increases with increasing ligand field strength of the bis(catecholate) bridging ligand in the order 1 < 2 < 3. The bromo substituents of Br4spiro4- increase the ligand field strength relative to spiro4-, while the stronger ligand field provided by thea4- arises from extension of the overlapping π-orbital system across the two catecholate units. This study highlights how SCO behavior in dinuclear complexes can be modulated by the bridging ligand, providing useful insights for the design of molecules that can be interconverted between more than two states.

Transition Temperature for Spin-Crossover Materials with the Mean Value Ensemble Hubbard-U Correction

Calculation of transition temperatures T1/2 for thermally driven spin-crossover in condensed phases is challenging, even with sophisticated state-of-the-art density functional approximations. The first issue is the accuracy of the adiabatic crossover energy difference ΔEHL between the low- and high-spin states of the bistable metal-organic complexes. The other is the proper inclusion of entropic contributions to the Gibbs free energy from the electronic and vibrational degrees of freedom. We discuss the effects of treatments of both contributions upon the calculation of thermochemical properties for a set of 20 spin-crossover materials using a Hubbard-U correction obtained from a reference ensemble spin-state. The U values obtained from a simplest bimolecular representation may overcorrect, somewhat, the ΔEHL values, hence giving somewhat excessive reduction of the T1/2 results with respect to their U = 0 values in the crystalline phase. We discuss the origins of the discrepancies by analyzing different sources of uncertainties. By use of a first-coordination-sphere approximation and the assumption that vibrational contributions from the outermost atoms in a metal-organic complex are similar in both low- and high-spin states, we achieve T1/2 results with the low-cost, widely used PBE generalized gradient density functional approximation comparable to those from the more costly, more sophisticated r2SCAN meta-generalized gradient approximation. The procedure is promising for use in high-throughput materials screening, because it combines rather low computational effort requirements with freedom from user manipulation of parameters.

Valence tautomerism in a cobalt-dioxolene complex containing an imidazolic ancillary ligand

This work reports the synthesis, structural, spectroscopic and magnetic investigation of two complexes, [Co(bmimapy)(3,5-DTBCat)]PF₆·H₂O (1) and [Co(bmimapy)(TCCat)]PF₆·H₂O (2), where bmimapy is an imidazolic tetradentate ancillary ligand and 3,5-DTBCat and TCCat are the 3,5-di-tert-butyl-catecholate and tetrachlorocatecholate anions, respectively. Their structures have been elucidated using single crystal X-ray diffraction, showing a pseudo-octahedral cobalt ion bound to a chelating dioxolene ligand and the ancillary bmimapy ligand in a folded conformation. Magnetometry displayed an entropy-driven, incomplete, Valence Tautomeric (VT) process for 1 in the 300-380 K temperature range, while 2 displayed a temperature independent, diamagnetic low-spin cobalt(iii)-catecholate charge distribution. This behaviour was interpreted on the basis of the cyclic voltammetric analysis, allowing the estimation of the free energy difference associated with the VT interconversion of +8 and +96 kJ mol^-1 for 1 and 2, respectively. A DFT analysis of this free energy difference highlighted the ability of the methyl-imidazole pendant arm of bmimapy favouring the onset of the VT phenomenon. This work introduces the imidazolic bmimapy ligand to the scientific community working in the field of valence tautomerism, increasing the library of ancillary ligands to prepare temperature switchable molecular magnetic materials.

One-Step versus Two-Step Valence Tautomeric Transitions in Tetraoxolene-Bridged Dinuclear Cobalt Compounds

The syntheses of valence tautomeric compounds with multistep transitions using new redox-active ligands are the long-term goal of the field of bistable materials. The redox-active tetraoxolene ligand, 2,7-di-tert-butylpyrene-4,5,9,10-tetraone (pyrene^Q-Q), is now developed to synthesize a pair of dinuclear compounds {[CoL₂]₂(pyrene^Sq-Sq)}[Co(CO)₄]₂·xCH₂Cl₂·2C₆H₅CH₃ (1, x = 2, L = 1,10-phenanthroline, phen; 2, x = 1.5, L = 2,2'-bipyridine, bpy). Variable-temperature magnetic susceptibilities and single-crystal X-ray diffraction measurements indicate a partial one-step valence tautomeric transition for 1 and a rare two-step valence tautomeric transition for 2, respectively. DFT calculation results are consistent with the experimental data, revealing the correlation between thermodynamic parameters and the one-step/two-step valence tautomeric behaviors.

Rich redox-activity and solvatochromism in a family of heteroleptic cobalt complexes

The combination of redox-active metals with redox-active ligands can lead to interesting charge transfer behaviours, including valence tautomerism and solvatochromism. With the aim of investigating a relatively underexplored redox-active metal/redox-active ligand combination, complexes [Co^II(acac)₂(X-BIAN)] (acac^- = acetylacetonate; X-BIAN = bis(4-X-phenyl)iminoacenaphthene; 1: X = -CF₃, 2: X = -Cl, 3: X = -H, 4: X = -Me) and [Co^III(acac)₂(Me-BIAN)]⁺ (5+) have been synthesised and characterised. At all temperatures investigated, and in both the solid and solution state, complexes 1-4 exist in a Co^II-BIAN⁰ charge distribution, while 5+ adopts a Co^III-BIAN⁰ charge distribution. In the case of 1-4, the potential Co^III-BIAN˙^- valence tautomer is inaccesible; the energy ordering between the ground Co^II-BIAN⁰ state and the excited Co^III-BIAN˙^- state must be reversed in order for an entropically driven interconversion to be possible. The energy gap between the states can be monitored via metal-to-ligand charge transfer bands in the visible region. We demonstrate tuning of this energy gap by varying the electronic properties of the BIAN ligand, as well as by controlling the molecular environment through solvent choice. Solvatochromic analysis, in combination with crystallographic evidence, allows elucidation of the specific solvent-solute interactions that govern the molecular behaviour of 1-4, affording insights that can inform potential future applications in sensing and switching.