Fine-tuning of the spin-crossover properties of Fe(III) complexes via ligand design.

Exploring the chemical space of a given ligand aiming to modulate its ligand field strength is a versatile strategy for the fine-tuning of physical properties such as the transition temperature (T1/2) of spin-crossover (SCO) complexes. The computational study presented herein aims at systematically exploring the extent to which the ligand substituent effects can modulate T1/2 in two families of Fe(III) SCO systems with a N4O2 coordination environment and at identifying the best descriptors for fast and accurate prediction of changes in T1/2 upon ligand functionalization. B3LYP* calculations show that the attachment of substituents to ß-ketoiminato fragments (L1) leads to drastic changes in T1/2, while functionalization of phenolato moieties (L2) allows for a finer degree of control over T1/2. Natural Bond Orbital (NBO) charges of the donor atoms, Hammett parameters for both para and meta-functionalization of L2, and Swain-Lupton parameters for L1 and para-functionalization of L2 have been found to be the suitable descriptors for predicting the changes in T1/2. Further analysis of the ligand-field splitting in such systems rationalizes the observed trends and shows that ligand substituents modify both the σ and π bonds between the Fe(III) center and the ligands. Thus, we provide simple yet reliable guidelines for the rational design of new SCO systems with specific values of T1/2 based on their ligand design.

Piano-Stool Ruthenium(II) Complexes with Delayed Cytotoxic Activity: Origin of the Lag Time.

We have recently reported a series of piano-stool ruthenium(II) complexes of the general formula [RuCl2(η6-arene)(P(1-pyrenyl)R2R3)] showing excellent cytotoxic activities (particularly when R2 = R3 = methyl). In the present study, new members of this family of compounds have been prepared with the objective to investigate the effect of the steric hindrance of a bulky phosphane ligand, namely diisopropyl(1-pyrenyl)phosphane (L), on exchange reactions involving the coordinated halides (X = Cl, I). Two η6-arene rings were used, i.e. η6-methyl benzoate (mba) and η6-p-cymene (p-cym), and four complexes were synthesized, namely [RuCl2(mba)(L)] (1Cl2iPr), [RuI2(mba)(L)] (1I2iPr), [RuCl2(p-cym)(L)] (2Cl2iPr), and [RuI2(p-cym)(L)] (2I2iPr). Unexpectedly, all of the complexes exhibited poor cytotoxic activities after 24 h of incubation with cells, in contrast to the related compounds previously reported. However, it was observed that aged DMSO solutions of 2I2iPr (from 2 to 7 days) exhibited better activities in comparison to freshly prepared solutions and that the activity improved over "aging" time. Thorough studies were therefore performed to uncover the origin of this lag time in the cytotoxicity efficiency. The data achieved clearly demonstrated that compounds 2I2iPr and 2Cl2iPr were undergoing a series of transformation reactions in DMSO (with higher rates for the iodido complex 2I2iPr), ultimately generating cyclometalated species through a mechanism involving DMSO as a coordinated proton abstractor. The cyclometalated complexes detected in solution were subsequently prepared; hence, pure [RuCl(p-cym)(κ2C-diisopropyl(1-pyrenyl)phosphane)] (3CliPr), [RuI(p-cym)(κ2C-diisopropyl(1-pyrenyl)phosphane)] (3IiPr), and [Ru(p-cym)(κS-dmso)(κ2C-diisopropyl(1-pyrenyl)phosphane)]PF6 (3dmsoiPr) were synthesized and fully characterized. Remarkably, 3CliPr, 3IiPr, and 3dmsoiPr are all very efficient cytotoxic agents, exhibiting slightly better activities in comparison to the chlorido noncyclometalated complexes [RuCl2(η6-arene)(P(1-pyrenyl)R2R3)] described in an earlier report. For comparison purposes, the iodido compounds [RuI2(mba)(dimethyl(1-pyrenyl)phosphane)] (1I2Me) and [RuI2(p-cym)(dimethyl(1-pyrenyl)phosphane)] (2I2Me), bearing the less hindered dimethyl(1-pyrenyl)phosphane ligand, have also been prepared. The cytotoxic and chemical behaviors of 1I2Me and 1I2Me were comparable to those of their chlorido counterparts reported previously.

Spin-Crossover Properties of an Iron(II) Coordination Nanohoop.

Addition of the bipyridyl-embedded cycloparaphenylene nanohoop bipy[9]CPP to [Fe{H2 B(pyz)2 }] (pyz=pyrazolyl) produces the distorted octahedral complex [Fe(bipy[9]CPP){H2 B(pyz)2 }2 ] (1). The molecular structure of 1 shows that the nanohoop ligand contains a non-planar bipy unit. Magnetic susceptibility measurements indicate spin-crossover (SCO) behaviour with a T1/2 of 130 K, lower than that of 160 K observed with the related compound [Fe(bipy){H2 B(pyz)2 }2 ] (2), which contains a conventional bipy ligand. A computational study of 1 and 2 reveals that the curvature of the nanohoop leads to the different SCO properties, suggesting that the SCO behaviour of iron(II) can be tuned by varying the size and diameter of the nanohoop.

Electronic Control of Spin-Crossover Properties in Four-Coordinate Bis(formazanate) Iron(II) Complexes.

The transition between spin states in d-block metal complexes has important ramifications for their structure and reactivity, with applications ranging from information storage materials to understanding catalytic activity of metalloenzymes. Tuning the ligand field (ΔO) by steric and/or electronic effects has provided spin-crossover compounds for several transition metals in the periodic table, but this has mostly been limited to coordinatively saturated metal centers in octahedral ligand environments. Spin-crossover complexes with low coordination numbers are much rarer. Here we report a series of four-coordinate, (pseudo)tetrahedral Fe(II) complexes with formazanate ligands and demonstrate how electronic substituent effects can be used to modulate the thermally induced transition between S = 0 and S = 2 spin states in solution. All six compounds undergo spin-crossover in solution with T1/2 above room temperature (300-368 K). While structural analysis by X-ray crystallography shows that the majority of these compounds are low-spin in the solid state (and remain unchanged upon heating), we find that packing effects can override this preference and give rise to either rigorously high-spin (6) or gradual spin-crossover behavior (5) also in the solid state. Density functional theory calculations are used to delineate the empirical trends in solution spin-crossover thermodynamics. In all cases, the stabilization of the low-spin state is due to the π-acceptor properties of the formazanate ligand, resulting in an "inverted" ligand field, with an approximate "two-over-three" splitting of the d-orbitals and a high degree of metal-ligand covalency due to metal → ligand π-backdonation. The computational data indicate that the electronic nature of the para-substituent has a different influence depending on whether it is present at the C-Ar or N-Ar rings, which is ascribed to the opposing effect on metal-ligand σ- and π-bonding.

[Co/Fe(α-Alkyl-tpdt)2]x-: Alkyl-Substituted Cobalt and Iron Bis-dithiolenethiophenic Complexes.

Tetraphenylphosphonium salts of Co and Fe complexes with alkyl-substituted, tert-butyl (tb), and isopropyl (dp) 2,3-thiophenedithiolate (α-tpdt) ligands, namely, TPP[Co(α-tb-tpdt)2] (3), TPP2[Fe(α-tb-tpdt)2]2 (4a-b), TPP[Co(α-dp-tpdt)2] (5), and TPP[Fe(α-dp-tpdt)2] (6) were prepared and characterized by cyclic voltammetry, single crystal X-ray diffraction, magnetic susceptibility measurements, and 57Fe Mössbauer spectroscopy. Compound 3 and 5 are isostructural with their Au and Ni analogues with a square-planar coordination geometry. Compound 4 presents two polymorphs (4a-b) both showing a Fe(III) bisdithiolene dimerization. The magnetic susceptibility of 3 and 5 exhibits behavior dominated by antiferromagnetic interactions, with room-temperature magnetic moments of 3.40 and 3.36 µB, respectively, indicating that these square-planar Co(III) complexes assume an intermediate spin electronic configuration (S = 1) as supported by multiconfigurational and DFT calculations.

Thermal spin crossover in Fe(ii) and Fe(iii). Accurate spin state energetics at the solid state.

The thermal spin crossover (SCO) phenomenon refers to an entropy-driven spin transition in some materials based on d6-d9 transition metal complexes. While its molecular origin is well known, intricate SCO behaviours are increasingly common, in which the spin transition occurs concomitantly to e.g. phase transformations, solvent absorption/desorption, or order-disorder processes. The computational modelling of such cases is challenging, as it requires accurate spin state energies in the solid state. Density Functional Theory (DFT) is the best framework, but most DFT functionals are unable to balance the spin state energies. While a few hybrid functionals perform better, they are still too expensive for solid-state minima searches in moderate-size systems. The best alternative is to dress cheap local (LDA) or semi-local (GGA) DFT functionals with a Hubbard-type correction (DFT+U). However, the parametrization of U is not straightforward due to the lack of reference values, and because ab initio parametrization methods perform poorly. Moreover, SCO complexes undergo notable structural changes upon transition, so intra- and inter-molecular interactions might play an important role in stabilizing either spin state. As a consequence, the U parameter depends strongly on the dispersion correction scheme that is used. In this paper, we parametrize U for nine reported SCO compounds (five based on FeII, 1-5 and four based on FeIII, 6-9) when using the D3 and D3-BJ dispersion corrections. We analyze the impact of the dispersion correction treatments on the SCO energetics, structure, and the unit cell dimensions. The average U values are different for each type of metal ion (FeIIvs. FeIII), and dispersion correction scheme (D3 vs. D3-BJ) but they all show excellent transferability, with mean absolute errors (MAE) below chemical accuracy (i.e. MAE <4 kJ mol-1). This enables a better description of SCO processes and, more generally, of spin state energetics, in materials containing FeII and FeIII ions.

Assessment of the SCAN Functional for Spin-State Energies in Spin-Crossover Systems.

The strongly constrained and appropriately normed (SCAN) functional has been tested toward the calculation of spin-state energy differences in a data set of 20 spin-crossover (SCO) systems, ranging from d4 to d7. Results show that the SCAN functional is able to correctly predict the low-spin state as the ground state for all systems, and the energy window provided by the calculations falls in the approximate range of energies that will allow for SCO to occur. Moreover, the SCAN functional can be used in periodic boundary condition calculations, accounting for the effect of collective crystal vibrations and counterions in the thermochemistry of the spin transition. Our results validate this functional as a potential method for in silico screening of new SCO systems at both, molecular and crystal-packed levels.

Electronic and Steric Control of the Spin-Crossover Behavior in [(CpR)2Mn] Manganocenes.

A computational study of the spin-crossover behavior in the family [(CpR)2Mn] (R = Me, iPr, tBu) is presented. Using the OPBE functional, the different electronic and steric effects over the metal's ligand field are studied, and trends in the spin-crossover-temperature (T1/2) behavior are presented in terms of the cyclopentadienyl (Cp) ligand functionalization. Our calculations outlined a delicate balance between both electronic and steric effects. While an increase in the number of electron-donating groups increases the spin-crossover temperature (T1/2) to the point that the transition is suppressed and only the low-spin state is observed, steric effects play an opposite role, increasing the distance between the Cp rings, which in turns shifts T1/2 to lower values, eventually stabilizing the high-spin state. Both effects can be rationalized by exploring the electronic structure of such systems in terms of the relevant d-based molecular orbitals.

Benchmarking Density Functional Methods for Calculation of State Energies of First Row Spin-Crossover Molecules.

A systematic study of the performance of several density functional methodologies to study spin-crossover (SCO) on first row transition metal complexes is reported. All functionals have been tested against several mononuclear systems containing first row transition metal complexes and exhibiting spin-crossover. Among the tested functionals, the hybrid meta-GGA functional TPSSh with a triple-ζ basis set including polarization functions on all atoms provides the best results across different metals and oxidation states, and its performance in both predicting the correct ground state and the right energy window for SCO to occur is quite satisfactory. The effect of some additional contributions, such as zero-point energies, relativistic effects, and intramolecular dispersion interactions, has been analyzed. The reported strategy thus expands the use of the TPSSh functional to other metals and oxidation states other than FeII, making it the method of choice to study SCO in first row transition metal complexes. Additionally, the presented results validate the potential use of the TPSSh functional for virtual screening of new molecules with SCO, or its use in the study of the electronic structure of such systems.

Mercurophilic interactions: a theoretical study on the importance of ligands.

In this work, a theoretical analysis of intermolecular HgHg contacts in the presence of different ligands is presented. A survey of structural databases to explore the geometrical preferences among experimental structures presenting short Hg(ii)Hg(ii) contacts reveals the main interaction topologies depending on the nature of the ligand. A benchmark study of several dispersion corrected-density functional methods is carried out to determine the optimal computational methodology for the theoretical study of such interaction. This methodology is later used in the study of several dimers of dicoordinated mercury compounds [HgX2]2 (X = Cl, Br, I, CN, CH3, CF3, SiH3, SiMe3) as well as to evaluate the stability of different interacting topologies in such dimers. It is found that, whenever possible, two mercury-containing molecules will preferentially join through intermolecular HgX donor-acceptor interactions and only ligands devoid of atoms with lone pairs available for such interactions can form dimers attached predominantly via mercurophilic interactions.

Molecular Mechanisms of Spin Crossover in the {Fe(pz)[Pt(CN)4]} Metal-Organic Framework upon Water Adsorption.

The rational design of multifunctional materials with properties that can be selectively controlled at the molecular level is key to the development and application of nanoscale devices. In this study, molecular dynamics simulations using ligand-field molecular mechanics are performed to elucidate, for the first time, the molecular mechanisms responsible for the variation of the spin-crossover properties of the {Fe(pz)[Pt(CN)4]} metal-organic framework upon water adsorption. The simulations demonstrate a direct relationship between the water loading adsorbed in the pores and the variation of the spin-crossover transition temperature, with the high-spin state of the material becoming gradually more stabilized as the number of adsorbed water molecules increases. The decrease of the spin-crossover temperature of {Fe(pz)[Pt(CN)4]} upon water adsorption predicted by the simulations is in agreement with the available experimental data and is traced back to the elongation of the bonds between the Fe(II) atoms and the organic linkers induced by interactions of the adsorbed water molecules with the framework.

Structure of the Reduced Copper Active Site in Preprocessed Galactose Oxidase: Ligand Tuning for One-Electron O2 Activation in Cofactor Biogenesis.

Galactose oxidase (GO) is a copper-dependent enzyme that accomplishes 2e- substrate oxidation by pairing a single copper with an unusual cysteinylated tyrosine (Cys-Tyr) redox cofactor. Previous studies have demonstrated that the post-translational biogenesis of Cys-Tyr is copper- and O2-dependent, resulting in a self-processing enzyme system. To investigate the mechanism of cofactor biogenesis in GO, the active-site structure of Cu(I)-loaded GO was determined using X-ray absorption near edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) spectroscopy, and density-functional theory (DFT) calculations were performed on this model. Our results show that the active-site tyrosine lowers the Cu potential to enable the thermodynamically unfavorable 1e- reduction of O2, and the resulting Cu(II)-O2•- is activated toward H atom abstraction from cysteine. The final step of biogenesis is a concerted reaction involving coordinated Tyr ring deprotonation where Cu(II) coordination enables formation of the Cys-Tyr cross-link. These spectroscopic and computational results highlight the role of the Cu(I) in enabling O2 activation by 1e- and the role of the resulting Cu(II) in enabling substrate activation for biogenesis.

Electronic Structure Modulation in an Exceptionally Stable Non-Heme Nitrosyl Iron(II) Spin-Crossover Complex.

The highly stable nitrosyl iron(II) mononuclear complex [Fe(bztpen)(NO)](PF6 )2 (bztpen=N-benzyl-N,N',N'-tris(2-pyridylmethyl)ethylenediamine) displays an S=1/2↔S=3/2 spin crossover (SCO) behavior (T1/2 =370 K, ΔH=12.48 kJ mol(-1) , ΔS=33 J K(-1) mol(-1) ) stemming from strong magnetic coupling between the NO radical (S=1/2) and thermally interconverted (S=0↔S=2) ferrous spin states. The crystal structure of this robust complex has been investigated in the temperature range 120-420 K affording a detailed picture of how the electronic distribution of the t2g -eg orbitals modulates the structure of the {FeNO}(7) bond, providing valuable magneto-structural and spectroscopic correlations and DFT analysis.

Theoretical Modeling of the Ligand-Tuning Effect over the Transition Temperature in Four-Coordinated Fe(II) Molecules.

Spin-crossover molecules are systems of great interest due to their behavior as molecular level switches, which makes them promising candidates for nanoscale memory devices, among other applications. In this paper, we report a computational study for the calculation of the transition temperature (T(1/2)), a key physical quantity in the characterization of spin-crossover systems, for the family of tetracoordinated Fe(II) transition-metal complexes of generic formula [PhB(MesIm)3FeNPR1R2R3]. Our calculations correctly reproduce the experimentally reported decrease in the T(1/2) with an increasing size of the phosphine and allow for the prediction of the T(1/2) in new members of the family that are not reported so far. More importantly, further insight into the factors that control the fine-tuning of the T(1/2) can be obtained by direct analysis of the underlying electronic structure in terms of the relevant molecular orbitals.

Non-Switching 1,2-Dithienylethene-based Diplatinum(II) Complex Showing High Cytotoxicity.

A diplatinum(II) complex was prepared from a new 1,2-dithienylethene-based ligand containing N-methylimidazole groups as metal-binding units. Reaction of the ligand 1,2-bis[2-methyl-5-(1-methyl-1H-imidazol-2-yl)-3-thienyl]-cyclopentene (L2(H)) with cis-dichlorobis(dimethylsulfoxido)platinum(II) generated the bimetallic complex trans-[Pt2Cl4(DMSO)2(L2(H))] (DMSO = dimethyl sulfoxide), whose DNA-interacting properties were investigated using different techniques. Cytotoxicity assays with various cancer cell lines showed that this compound is active, with IC50 values in the micromolar range. Surprisingly, the diplatinum(II) complex does not exhibit the anticipated photoswitching properties; indeed, UV irradiation does not lead to the photocyclization of the ligand L2(H) or of the metal complex. Computational studies were performed and revealed significant differences in the electronic structure of L2(H) compared with L1(H) (i.e., 1,2-bis[2-methyl-5-(4-pyridyl)-3-thienyl]-cyclopentene, which exhibits photoswitching properties), in terms of the relevant molecular orbitals involved in the UV-vis absorption features, which ultimately is responsible for the inertia of L2(H) toward photocyclization.

Theoretical modeling of spin crossover in metal-organic frameworks: [Fe(pz)2Pt(CN)4] as a case study.

Metal-organic frameworks (MOFs) with spin-crossover behavior are promising materials for applications in memory storage and sensing devices. A key parameter that characterizes these materials is the transition temperature T1/2, defined as the temperature with equal populations of low-spin and high-spin species. In this study, we describe the development, implementation, and application of a novel hybrid Monte Carlo/molecular dynamics method that builds upon the Ligand Field Molecular Mechanics approach and enables the modeling of spin-crossover properties in bulk materials. The new methodology is applied to the study of a spin-crossover MOF with molecular formula [Fe(pz)2Pt(CN)4] (pz = pyrazine). The total magnetic moment of the material is determined as a function of the temperature from direct calculations of the relative equilibrium populations of both low-spin and high-spin states of each Fe(II) center of the framework. The T1/2 value, calculated from the temperature dependence of the magnetization curve, is in good agreement with the available experimental data. A comparison between the spin-crossover behavior of the isolated secondary building block of the framework and the bulk material is presented, which reveals the origin of the different spin-crossover properties of the isolated molecular system and corresponding MOF structure.

Tuning the spin-crossover properties of [Fe2] metal-organic cages.

A computational study on the interplay between ligand functionalization and guest effects on the transition temperature (T1/2) in the [Fe2(L1R)3]@X (L1 = 1,3-bis-(3-(pyridin-2-yl)-1H-pyrazol-5-yl)benzene, X = H-, F-, Cl-, Br-, I- and [BF4]-, R = H, F, or CH3) family of metal-organic cages (MOCs) is presented. Our results indicate that ligand functionalization with electron-donating or electron-withdrawing groups can significantly impact the T1/2 as expected, while the guest effect in lowering the T1/2 has a linear correlation with the increasing guest size. More importantly, small guests can move away from the center of the cavity, thus enhancing the two-step characteristic of the transition. All the data can be understood by analyzing the underlying electronic structure of the studied systems in terms of the relevant d-based molecular orbitals. These results can help in the rational design of new MOCs that can operate as sensors at specific temperatures, thus accelerating the discovery of new SCO devices with tailored properties.

Elucidating the exchange interactions in a {GdIIICuII4} propellor.

The multinucleating ligand 2,2'-(propane-1,3-diyldiimino)bis[2-(hydroxymethyl)-propane-1,3-diol] (bis-tris propane, H6L) is used in the design of a new family of 3d-4f complexes that display an unusual {LnCu4} four-blade propeller topology. We report the synthesis, structure and magnetic characterisation of [LnCu4(H4L)4](Cl)2(ClO4)·6CH3OH, where Ln = Gd (1), Tb (2), Dy (3), La (4). Previously we have used CH3COO- and NO3- as co-ligands with bis-tris propane, but here the use of Cl- and ClO4- leads to coordination of four {Cu(H4L)} units around the central Ln ion. A magneto-structural analysis reveals that the geometrical arrangement of the Cu(II) centres defined by the H4L2- ligands controls the magnetic communication between the different metal centres. DFT calculations performed on the isotropic (Gd) and diamagnetic (La) systems 1 and 4 help to unravel the intriguing exchange interactions.

Synthesis and characterization of highly diluted polyanionic iron(II) spin crossover systems.

Spin crossover (SCO) complexes, through their reversible spin transition under external stimuli, can work as switchable memory materials. Here, we present a protocol for the synthesis and characterization of a specific polyanionic iron SCO complex and its diluted systems. We describe steps for its synthesis and the determination of crystallographic structure of the SCO complex in diluted systems. We then detail a range of spectroscopic and magnetic techniques employed to monitor the spin state of the SCO complex in both diluted solid- and liquid-state systems. For complete details on the use and execution of this protocol, please refer to Galán-Mascaros et al.1.

Theoretical prediction of spin-crossover temperatures in ligand-driven light-induced spin change systems.

Spin-crossover compounds exhibit two alternative spin states with distinctive chemical and physical properties, a particular feature that makes them promising materials for nanotechnological applications as memory or display devices. A key parameter that characterizes these compounds is the spin-crossover temperature, T(1/2), defined as the temperature with equal populations of high and low-spin species. In this study, a theoretical/computational approach is described for the calculation of T(1/2) for the trans-[Fe(styrylpyridine)(4)(NCX)(2)] (X = S, Se, and BH(3), styrylpyridine in the trans configuration) ligand driven light-induced spin change (LD-LISC) complexes. In all cases, the present calculations provide an accurate description of both structural and electronic properties of the LD-LISC complexes and, importantly, predict spin-crossover temperatures in good agreement with the corresponding experimental data. Fundamental insights into the dependence of T(1/2) on the nature of the axial ligands are obtained from the direct analysis of the underlying electronic structure in terms of the relevant molecular orbitals.