Insights into the Spin-State Transitions in [Fe(tpy)2]2+: Importance of the Terpyridine Rocking Motion

Activating a high-spin iron(II) complex to thermal spin-crossover with an inert non-isomorphous molecular dopant

[Fe(bpp)2][ClO4]2 (bpp = 2,6-bis{pyrazol-1-yl}pyridine; monoclinic, C2/c) is high-spin between 5-300 K, and crystallises with a highly distorted molecular geometry that lies along the octahedral-trigonal prismatic distortion pathway. In contrast, [Ni(bpp)2][ClO4]2 (monoclinic, P21) adopts a more regular, near-octahedral coordination geometry. Gas phase DFT minimisations (ω-B97X-D/6-311G**) of [M(bpp)2]2+ complexes show the energy penalty associated with that coordination geometry distortion runs as M2+ = Fe2+ (HS) ≈ Mn2+ (HS) < Zn2+ ≈ Co2+ (HS) ≲ Cu2+ ≪ Ni2+ ≪ Ru2+ (LS; HS = high-spin, LS = low-spin). Slowly crystallised solid solutions [FexNi1-x(bpp)2][ClO4]2 with x = 0.53 (1a) and 0.74 (2a) adopt the P21 lattice, while x = 0.87 (3a) and 0.94 (4a) are mixed-phase materials with the high-spin C2/c phase as the major component. These materials exhibit thermal spin-transitions at T½ = 250 ± 1 K which occurs gradually in 1a, and abruptly and with narrow thermal hysteresis in 2a-4a. The transition proceeds to 100% completeness in 1a and 2a; that is, the 26% Ni doping in 2a is enough to convert high-spin [Fe(bpp)2][ClO4]2 into a cooperative, fully SCO-active material. These results were confirmed crystallographically for 1a and 2a, which revealed similarities and differences between these materials and the previously published [FexNi1-x(bpp)2][BF4]2 series. Rapidly precipitated powders with the same compositions (1b-4b) mostly resemble 1a-4a, except that 2b is a mixed-phase material; 2b-4b also contain a fraction of amorphous solid in addition to the two crystal phases. The largest iron fraction that can be accommodated by the P21 phase in this system is 0.7 ± 0.1.

A Survey of the Angular Distortion Landscape in the Coordination Geometries of High-Spin Iron(II) 2,6-Bis(pyrazolyl)pyridine Complexes

Reaction of 2,4,6-trifluoropyridine with sodium 3,4-dimethoxybenzenethiolate and 2 equiv of sodium pyrazolate in tetrahydrofuran at room temperature affords 4-(3,4-dimethoxyphenylsulfanyl)-2,6-di(pyrazol-1-yl)pyridine (L), in 30% yield. The iron(II) complexes [FeL2][BF4]2 (1a) and [FeL2][ClO4]2 (1b) are high-spin with a highly distorted six-coordinate geometry. This structural deviation from ideal D2d symmetry is common in high-spin [Fe(bpp)2]2+ (bpp = di{pyrazol-1-yl}pyridine) derivatives, which are important in spin-crossover materials research. The magnitude of the distortion in 1a and 1b is the largest yet discovered for a mononuclear complex. Gas-phase DFT calculations at the ω-B97X-D/6-311G** level of theory identified four minimum or local minimum structural pathways across the distortion landscape, all of which are observed experimentally in different complexes. Small distortions from D2d symmetry are energetically favorable in complexes with electron-donating ligand substituents, including sulfanyl groups, which also have smaller energy penalties associated with the lowest energy distortion pathway. Natural population analysis showed that these differences reflect greater changes to the Fe-N{pyridyl} σ-bonding as the distortion proceeds, in the presence of more electron-rich pyridyl donors. The results imply that [Fe(bpp)2]2+ derivatives with electron-donating pyridyl substituents are more likely to undergo cooperative spin transitions in the solid state. The high-spin salt [Fe(bpp)2][CF3SO3]2, which also has a strong angular distortion, is also briefly described and included in the analysis.

On the use of vibronic coherence to identify reaction coordinates for ultrafast excited-state dynamics of transition metal-based chromophores

The question of whether one can use information from quantum coherence as a means of identifying vibrational degrees of freedom that are active along an excited-state reaction coordinate is discussed. Specifically, we are exploring the notion of whether quantum oscillations observed in single-wavelength kinetics data exhibiting coherence dephasing times that are intermediate between that expected for either pure electronic or pure vibrational dephasing are vibronic in nature and therefore may be coupled to electronic state-to-state evolution. In the case of a previously published Fe(II) polypyridyl complex, coherences observed subsequent to ¹A₁ → ¹MLCT excitation were linked to large-amplitude motion of a portion of the ligand framework; dephasing times on the order of 200-300 fs suggested that these degrees of freedom could be associated with ultrafast (∼100 fs) conversion from the initially formed MLCT excited state to lower-energy, metal-centered ligand-field excited state(s) of the compound. Incorporation of an electronically benign but sterically restrictive Cu(I) ion into the superstructure designed to interfere with this motion yielded a compound exhibiting a ∼25-fold increase in the compound's MLCT lifetime, a result that was interpreted as confirmation of the initial hypothesis. However, new data acquired on a different chemical system - Cr(acac')₃ (where acac' represents various derivatives of acetylacetonate) - yielded results that call into question this same hypothesis. Coherences observed subsequent to ⁴A₂ → ⁴T₂ ligand-field excitation on a series of molecules implicated similar vibrational degrees of freedom across the series, but exhibited dephasing times ranging from 340 fs to 2.5 ps without any clear correlation to the dynamics of excited-state evolution in the system. Taken together, the results obtained on both of these chemical platforms suggest that while identification of coherences can indeed point to degrees of freedom that should be considered as candidate modes for defining reaction trajectories, our understanding of the factors that determine the interplay across coherences, dephasing times, and electronic and geometric structure is insufficient at the present time to view this parameter as a robust metric for differentiating active versus spectator modes for ultrafast dynamics.

Chiral control of spin-crossover dynamics in Fe(II) complexes

Iron-based spin-crossover complexes hold tremendous promise as multifunctional switches in molecular devices. However, real-world technological applications require the excited high-spin state to be kinetically stable-a feature that has been achieved only at cryogenic temperatures. Here we demonstrate high-spin-state trapping by controlling the chiral configuration of the prototypical iron(II)tris(4,4'-dimethyl-2,2'-bipyridine) in solution, associated for stereocontrol with the enantiopure Δ- or Λ-enantiomer of tris(3,4,5,6-tetrachlorobenzene-1,2-diolato-κ²O¹,O²)phosphorus(V) (P(O₂C₆Cl₄)₃^- or TRISPHAT) anions. We characterize the high-spin-state relaxation using broadband ultrafast circular dichroism spectroscopy in the deep ultraviolet in combination with transient absorption and anisotropy measurements. We find that the high-spin-state decay is accompanied by ultrafast changes of its optical activity, reflecting the coupling to a symmetry-breaking torsional twisting mode, contrary to the commonly assumed picture. The diastereoselective ion pairing suppresses the vibrational population of the identified reaction coordinate, thereby achieving a fourfold increase of the high-spin-state lifetime. More generally, our results motivate the synthetic control of the torsional modes of iron(II) complexes as a complementary route to manipulate their spin-crossover dynamics.

Accurate calculation of spin-state energy gaps in Fe(III) spin-crossover systems using density functional methods

Fe(III) complexes are receiving ever-increasing attention as spin crossover (SCO) systems because they are usually air stable, as opposed to Fe(II) complexes, which are prone to oxidation. Here, we present the first systematic study exclusively devoted to assess the accuracy of several exchange-correlation functionals when it comes to predicting the energy gap between the high-spin (S = 5/2) and the low-spin (S = 1/2) states of Fe(III) complexes. Using a dataset of 24 different Fe(III) hexacoordinated complexes, it is demonstrated that the B3LYP* functional is an excellent choice not only for predicting spin-state energy gaps for Fe(III) complexes undergoing spin-transitions but also for discriminating Fe(III) complexes that are either low- or high-spin in the whole range of temperatures. Our benchmark study has led to the identification of a very versatile Fe(III) compound whose SCO properties can be engineered upon changing a single axial ligand. Overall, this work demonstrates that B3LYP* is a reliable functional for screening new spin-crossover systems with tailored properties.

Resonance Raman Excitation Profiles of Fe(II)-Terpyridine Complexes: Electronic Effects of Ligand Modifications

Metal 2,2':6',2″-terpyridine (tpy) complexes are readily used as building blocks in metallo-supramolecular polymers that stand out for their photophysical properties in solar energy assemblies. Furthermore, Resonance Raman (RR) excitation profiles are sensitive indicators of the electronic properties of chromophores. Previously, using RR spectroscopy, we studied the [Fe(tpy)₂]²⁺ complex and metallo-supramolecular polymers formed by tpy derivatives and Fe(II) ions. Here, we compare RR spectra of iron (Fe(II)) complexes with 4'-substituted tpy ligands─[Fe(4'-R-tpy)₂]²⁺, with R = H (1a), Cl (2a), 4-chlorophenyl (3a), and 2-thienyl (4a) to describe changes in their electronic structure after functionalization. By combining theoretical calculations, RR, and UV/vis spectra, we elucidated differences in the RR excitation profiles of 1a, 2a, and 4a complexes. In all Raman modes, complexes 1a and 2a showed maximal enhancement only at 532 nm excitation, whereas complex 4a exhibited maximal enhancement selectively at either 532 or 633 nm excitations. Based on our calculations, the mixed metal/ligand character of the highest occupied molecular orbital (HOMO) of 4a complex manifests itself through selective enhancement of vibration modes, mainly localized on the 2-thienyl unit at 633 nm excitation, which may explain the unique behavior of this complex. Therefore, complex 4a is a prospective candidate for further detailed photophysical explorations toward developing sensitizers for solar cells.

Interpolation Methods for Molecular Potential Energy Surface Construction

The concept of a potential energy surface (PES) is one of the most important concepts in modern chemistry. A PES represents the relationship between the chemical system's energy and its geometry (i.e., atom positions) and can provide useful information about the system's chemical properties and reactivity. Construction of accurate PESs with high-level theoretical methodologies, such as density functional theory, is still challenging due to a steep increase in the computational cost with the increase of the system size. Thus, over the past few decades, many different mathematical approaches have been applied to the problem of the cost-efficient PES construction. This article serves as a short overview of interpolative methods for the PES construction, including global polynomial interpolation, trigonometric interpolation, modified Shepard interpolation, interpolative moving least-squares, and the automated PES construction derived from these.

Halogenation affects driving forces, reorganization energies and "rocking" motions in strained [Fe(tpy)2]2+ complexes

Controlling the energetics of spin crossover (SCO) in Fe(II)-polypyridine complexes is critical for designing new multifunctional materials or tuning the excited-state lifetimes of iron-based photosensitizers. It is well established that the Fe-N "breathing" mode is important for intersystem crossing from the singlet to the quintet state, but this does not preclude other, less obvious, structural distortions from affecting SCO. Previous work has shown that halogenation at the 6 and 6'' positions of tpy (tpy = 2,2';6',2''-terpyridine) in [Fe(tpy)₂]²⁺ dramatically increased the lifetime of the excited MLCT state and also had a large impact on the ground state spin-state energetics. To gain insight into the origins of these effects, we used density functional theory calculations to explore how halogenation impacts spin-state energetics and molecular structure in this system. Based on previous work we focused on the ligand "rocking" motion associated with SCO in [Fe(tpy)₂]²⁺ by constructing one-dimensional potential energy surfaces (PESs) along the tpy rocking angle for various spin states. It was found that halogenation has a clear and predictable impact on ligand rocking and spin-state energetics. The rocking is correlated to numerous other geometrical distortions, all of which likely affect the reorganization energies for spin-state changes. We have quantified trends in reorganization energy and also driving force for various spin-state changes and used them to interpret the experimentally measured excited-state lifetimes.

Reduced-dimensional surface hopping with offline-online computations

Molecular dynamics simulations often classically evolve the nuclear geometry on adiabatic potential energy surfaces (PESs), punctuated by random hops between energy levels in regions of strong coupling, in an algorithm known as surface hopping. However, the computational expense of integrating the geometry on a full-dimensional PES and computing the required couplings can quickly become prohibitive as the number of atoms increases. In this work, we describe a method for surface hopping that uses only important reaction coordinates, performs all expensive evaluations of the true PESs and couplings only once before simulating dynamics (offline), and then queries the stored values during the surface hopping simulation (online). Our Python codes are freely available on GitHub. Using photodissociation of azomethane as a test case, this method is able to reproduce experimental results that have thus far eluded ab initio surface hopping studies.

Efficient Approximation of Potential Energy Surfaces with Mixed-Basis Interpolation

The potential energy surface (PES) describes the energy of a chemical system as a function of its geometry and is a fundamental concept in modern chemistry. A PES provides much useful information about the system, including the structures and energies of various stationary points, such as stable conformers (local minima) and transition states (first-order saddle points) connected by a minimum-energy path. Our group has previously produced surrogate reduced-dimensional PESs using sparse interpolation along chemically significant reaction coordinates, such as bond lengths, bond angles, and torsion angles. These surrogates used a single interpolation basis, either polynomials or trigonometric functions, in every dimension. However, relevant molecular dynamics (MD) simulations often involve some combination of both periodic and nonperiodic coordinates. Using a trigonometric basis on nonperiodic coordinates, such as bond lengths, leads to inaccuracies near the domain boundary. Conversely, polynomial interpolation on the periodic coordinates does not enforce the periodicity of the surrogate PES gradient, leading to nonconservation of total energy even in a microcanonical ensemble. In this work, we present an interpolation method that uses trigonometric interpolation on the periodic reaction coordinates and polynomial interpolation on the nonperiodic coordinates. We apply this method to MD simulations of possible isomerization pathways of azomethane between cis and trans conformers. This method is the only known interpolative method that appropriately conserves total energy in systems with both periodic and nonperiodic reaction coordinates. In addition, compared to all-polynomial interpolation, the mixed basis requires fewer electronic structure calculations to obtain a given level of accuracy, is an order of magnitude faster, and is freely available on GitHub.

Evidence for a lowest energy 3MLCT excited state in [Fe(tpy)(CN)3]. Chem Commun (Camb) 2021;57:4658-4661. [PMID: 33977987 DOI: 10.1039/d1cc01090e]

Transient absorption data of [FeII(tpy)(CN)3]- reveals spectroscopic signatures indicative of 3MLCT with a ∼10 ps kinetic component. These data are supported by DFT and TD-DFT calculations, which show that excited state ordering is responsive to the number of cyanide ligands on the complex.

Ultrafast photoinduced dynamics in Prussian blue analogues

A review on ultrafast photoinduced processes in molecule-based magnets with an emphasis on Prussian blue analogues.

SCC-DFTB Parameters for Fe-C Interactions

We present an optimized density-functional tight-binding (DFTB) parameterization for iron-based complexes based on the popular trans3d set of parameters. The transferability of the original and optimized parameterizations is assessed using a set of 50 iron complexes, which include carbonyl, cyanide, polypyridine, and cyclometalated ligands. DFTB-optimized structures predicted using the trans3d parameters show a good agreement with both experimental crystal geometries and density functional theory (DFT)-optimized structures for Fe-N bond lengths. Conversely, Fe-C bond lengths are systematically overestimated. We improve the accuracy of Fe-C interactions by truncating the Fe-O repulsive potential and reparameterizing the Fe-C repulsive potential using a training set of six isolated iron complexes. The new trans3d*-LANLFeC parameter set can produce accurate Fe-C bond lengths in both geometry optimizations and molecular dynamics (MD) simulations, without significantly affecting the accuracy of Fe-N bond lengths. Moreover, the potential energy curves of Fe-C interactions are considerably improved. This improved parameterization may open the door to accurate MD simulations at the DFTB level of theory for large systems containing iron complexes, such as sensitizer-semiconductor assemblies in dye-sensitized solar cells, that are not easily accessible with DFT approaches because of the large number of atoms.

Quantum-chemistry-aided ligand engineering for potential molecular switches: changing barriers to tune excited state lifetimes

Substitution of terpyridine at the 4' position with electron withdrawing and donating groups is used to tune the quintet lifetime of its iron(ii) complex. DFT calculations suggest that the energy barrier between the quintet and singlet states can be altered significantly upon substitution, inducing a large variation of the lifetime of the photoexcited quintet state. This prediction was experimentally verified by transient optical absorption spectroscopy and good agreement with the trend expected from the calculations was found. This demonstrates that the potential energy landscape can indeed be rationally tailored by relevant modifications based on DFT predictions. This result should pave the way to advancing efficient theory-based ligand engineering of functional molecules to a wide range of applications.

Exploring the light-induced dynamics in solvated metallogrid complexes with femtosecond pulses across the electromagnetic spectrum

Oligonuclear complexes of d⁴-d⁷ transition metal ion centers that undergo spin-switching have long been developed for their practical role in molecular electronics. Recently, they also have appeared as promising photochemical reactants demonstrating improved stability. However, the lack of knowledge about their photophysical properties in the solution phase compared to mononuclear complexes is currently hampering their inclusion into advanced light-driven reactions. In the present study, the ultrafast photoinduced dynamics in a solvated [2 × 2] iron(II) metallogrid complex are characterized by combining measurements with transient optical-infrared absorption and x-ray emission spectroscopy on the femtosecond time scale. The analysis is supported by density functional theory calculations. The photocycle can be described in terms of intra-site transitions, where the Fe^II centers in the low-spin state are independently photoexcited. The Franck-Condon state decays via the formation of a vibrationally hot high-spin (HS) state that displays coherent behavior within a few picoseconds and thermalizes within tens of picoseconds to yield a metastable HS state living for several hundreds of nanoseconds. Systematic comparison with the closely related mononuclear complex [Fe(terpy)₂]²⁺ reveals that nuclearity has a profound impact on the photoinduced dynamics. More generally, this work provides guidelines for expanding the integration of oligonuclear complexes into new photoconversion schemes that may be triggered by ultrafast spin-switching.

Photophysics and Photochemistry of Iron Carbene Complexes for Solar Energy Conversion and Photocatalysis

Earth-abundant first row transition metal complexes are important for the development of large-scale photocatalytic and solar energy conversion applications. Coordination compounds based on iron are especially interesting, as iron is the most common transition metal element in the Earth’s crust. Unfortunately, iron-polypyridyl and related traditional iron-based complexes generally suffer from poor excited state properties, including short excited-state lifetimes, that make them unsuitable for most light-driven applications. Iron carbene complexes have emerged in the last decade as a new class of coordination compounds with significantly improved photophysical and photochemical properties, that make them attractive candidates for a range of light-driven applications. Specific aspects of the photophysics and photochemistry of these iron carbenes discussed here include long-lived excited state lifetimes of charge transfer excited states, capabilities to act as photosensitizers in solar energy conversion applications like dye-sensitized solar cells, as well as recent demonstrations of promising progress towards driving photoredox and photocatalytic processes. Complementary advances towards photofunctional systems with both Fe(II) complexes featuring metal-to-ligand charge transfer excited states, and Fe(III) complexes displaying ligand-to-metal charge transfer excited states are discussed. Finally, we outline emerging opportunities to utilize the improved photochemical properties of iron carbenes and related complexes for photovoltaic, photoelectrochemical and photocatalytic applications.

Approximating Periodic Potential Energy Surfaces with Sparse Trigonometric Interpolation

The potential energy surface (PES) describes the energy of a chemical system as a function of its geometry and is a fundamental concept in computational chemistry. A PES provides much useful information about the system, including the structures and energies of various stationary points, such as local minima, maxima, and transition states. Construction of full-dimensional PESs for molecules with more than 10 atoms is computationally expensive and often not feasible. Previous work in our group used sparse interpolation with polynomial basis functions to construct a surrogate reduced-dimensional PESs along chemically significant reaction coordinates, such as bond lengths, bond angles, and torsion angles. However, polynomial interpolation does not preserve the periodicity of the PES gradient with respect to angular components of geometry, such as torsion angles, which can lead to nonphysical phenomena. In this work, we construct a surrogate PES using trigonometric basis functions, for a system where the selected reaction coordinates all correspond to the torsion angles, resulting in a periodically repeating PES. We find that a trigonometric interpolation basis not only guarantees periodicity of the gradient but also results in slightly lower approximation error than polynomial interpolation.

Molecular Dynamics Simulations on Relaxed Reduced-Dimensional Potential Energy Surfaces

Molecular dynamics (MD) simulations with full-dimensional potential energy surfaces (PESs) obtained from high-level ab initio calculations are frequently used to model reaction dynamics of small molecules (i.e., molecules with up to 10 atoms). Construction of full-dimensional PESs for larger molecules is, however, not feasible since the number of ab initio calculations required grows rapidly with the increase of dimension. Only a small number of coordinates are often essential for describing the reactivity of even very large systems, and reduced-dimensional PESs with these coordinates can be built for reaction dynamics studies. While analytical methods based on transition-state theory framework are well established for analyzing the reduced-dimensional PESs, MD simulation algorithms capable of generating trajectories on such surfaces are more rare. In this work, we present a new MD implementation that utilizes the relaxed reduced-dimensional PES for standard microcanonical (NVE) and canonical (NVT) MD simulations. The method is applied to the pyramidal inversion of a NH₃ molecule. The results from the MD simulations on a reduced, three-dimensional PES are validated against the ab initio MD simulations, as well as MD simulations on full-dimensional PES and experimental data.

Designing air-stable cyclometalated Fe(ii) complexes: stabilization via electrostatic effects

Substitution of EWGs onto the cyclometelated iron complexes electrostatically stabilizes the Fe(ii) center while still preserving the increased ligand field strength.

Band-selective dynamics in charge-transfer excited iron carbene complexes

A combination of ultrafast spectroscopy and DFT/TD-DFT calculations of a recently synthesised iron carbene complex elucidates the ultrafast excited state evolution processes in these systems.

Theoretical study on the excited state decay properties of iron(ii) polypyridine complexes substituted by bromine and chlorine

Transition metal iron(ii) polypyridyl complexes with quintet ground states were deeply investigated by density functional theory (DFT) and time-dependent density functional theory (TDDFT). Compared with the parent complex [Fe(tpy)₂]²⁺ (tpy = 2,2′:6′,2′′-terpyridine), the ground states of the complexes substituted by halogen atoms changed from singlet states to quintet states with rare high spin excited state lifetimes. The substituted complex [Fe(dbtpy)₂]²⁺ (1) results in a high spin metal–ligand charge transfer lifetime of 17.4 ps, which is 1.4 ps longer than that of [Fe(dctpy)₂]²⁺ (2) with the substitution of chlorine atoms. The reason for this is explored by a combination of electronic structures, absorption spectra, extended transition state coupled with natural orbitals for chemical valence (ETS-NOCV) studies and potential energy curves (PECs). The distortion of 1 in the angles and dihedrals of the ligands is slightly larger than that in 2, although the average metal–ligand bond lengths of the latter are larger. The twisted octahedron decreases the interactions between the d orbitals of iron(ii) and the n/π orbitals of the ligands. Compared with 2, the enlarged energy gaps among the different PECs of 1 and the increased energy crossing points caused by the larger distortion result in the increase of its excited state lifetime. The different pairwise orbital interaction contributions between the metal center and the ligands in their singlet states are qualitatively estimated by ETS-NOCV. The results show that the substitution of bromine atoms will decrease the electrostatic attraction between the metal and ligands but not significantly impact the orbital interactions.

Transition metal iron(ii) halogen substituted polypyridyl complexes with quintet ground states were deeply investigated by density functional theory (DFT) and time-dependent density functional theory (TDDFT).

Insights into the excited state dynamics of Fe(ii) polypyridyl complexes from variable-temperature ultrafast spectroscopy

In an effort to better define the nature of the nuclear coordinate associated with excited state dynamics in first-row transition metal-based chromophores, variable-temperature ultrafast time-resolved absorption spectroscopy has been used to determine activation parameters associated with ground state recovery dynamics in a series of low-spin Fe(ii) polypyridyl complexes. Our results establish that high-spin (⁵T₂) to low-spin (¹A₁) conversion in complexes of the form [Fe(4,4'-di-R-2,2'-bpy')₃]²⁺ (R = H, CH₃, or tert-butyl) is characterized by a small but nevertheless non-zero barrier in the range of 300-350 cm^-1 in fluid CH₃CN solution, a value that more than doubles to ∼750 cm^-1 for [Fe(terpy)₂]²⁺ (terpy = 2,2':6',2''-terpyridine). The data were analyzed in the context of semi-classical Marcus theory. Changes in the ratio of the electronic coupling to reorganization energy (specifically, H _ab ⁴/λ) reveal an approximately two-fold difference between the [Fe(bpy')₃]²⁺ complexes (∼1/30) and [Fe(terpy)₂]²⁺ (∼1/14), suggesting a change in the nature of the nuclear coordinate associated with ground state recovery between these two types of complexes. These experimentally-determined ratios, along with estimates for the ⁵T₂/¹A₁ energy gap, yield electronic coupling values between these two states for the [Fe(bpy')₃]²⁺ series and [Fe(terpy)₂]²⁺ of 4.3 ± 0.3 cm^-1 and 6 ± 1 cm^-1, respectively, values that are qualitatively consistent with the second-order nature of high-spin/low-spin coupling in a d⁶ ion. In addition to providing useful quantitative information on these prototypical Fe(ii) complexes, these results underscore the utility of variable-temperature spectroscopic measurements for characterizing ultrafast excited state dynamics in this class of compounds.

Impact of Ligand Substitutions on Multielectron Redox Properties of Fe Complexes Supported by Nitrogenous Chelates

Redox flow batteries (RFBs) have recently been recognized as a potentially viable technology for scalable energy storage. To take full advantage of RFBs, one possible approach for achieving high energy densities is to maximize a number of redox events by utilizing charge carriers capable of multiple one-electron transfers within the electrochemical window of solvent. However, past efforts to develop more efficient electrolytes for nonaqueous RFBs have mostly been empirical. In this manuscript, we shed light on design principles by theoretically investigating the effects of systematically substituting pyridyl moieties with imine ligands within a series of Fe complexes with some experimental validation. We found that such replacement is an effective strategy for reducing the molecular weight-to-charge ratios of these complexes. Simultaneously, calculations suggest that the reduction potentials and ligand-based redox activity of such substituted N-heterocyclic Fe compounds might be maintained within their +4 → -1 charge states. Additionally, by theoretically examining the role of coordination geometry, vis-à-vis reducing the number of redox noninnocent ligands within the first coordination sphere, we have demonstrated that Fe complexes with one such ligand were also capable of supporting multielectron reduction events and exhibited reduction potentials similar to their parent analogs supported by two or three of the same multidentate ligands. However, some differences in redox nature within the lower (+2 → -1) charge states were also noticed. Specifically, complexes containing two bidentate ligands, or one tridentate ligand, exhibited ligand-based reductions, whereas compounds with one bidentate ligand exhibited metal-centered reductions. The current results pave the way toward the design of the next-generation of Fe complexes with lower molecular weights and greater stored energy for redox flow batteries.

Iron-based trinuclear metal-organic nanostructures on a surface with local charge accumulation

Coordination chemistry relies on harnessing active metal sites within organic matrices. Polynuclear complexes-where organic ligands bind to several metal atoms-are relevant due to their electronic/magnetic properties and potential for functional reactivity pathways. However, their synthesis remains challenging; few geometries and configurations have been achieved. Here, we synthesise-via supramolecular chemistry on a noble metal surface-one-dimensional metal-organic nanostructures composed of terpyridine (tpy)-based molecules coordinated with well-defined polynuclear iron clusters. Combining low-temperature scanning probe microscopy and density functional theory, we demonstrate that the coordination motif consists of coplanar tpy's linked via a quasi-linear tri-iron node in a mixed (positive-)valence metal-metal bond configuration. This unusual linkage is stabilised by local accumulation of electrons between cations, ligand and surface. The latter, enabled by bottom-up on-surface synthesis, yields an electronic structure that hints at a chemically active polynuclear metal centre, paving the way for nanomaterials with novel catalytic/magnetic functionalities.

Synthesis and Computational Study of a Pyridylcarbene Fe(II) Complex: Unexpected Effects of fac/ mer Isomerism in Metal-to-Ligand Triplet Potential Energy Surfaces

The synthesis and the steady-state absorption spectrum of a new pyridine-imidazolylidene Fe(II) complex (Fe-NHC) are presented. A detailed mechanism of the triplet metal-to-ligand charge-transfer states decay is provided on the basis of minimum energy path (MEP) calculations used to connect the lowest-lying singlet, triplet, and quintet state minima. The competition between the different decay pathways involved in the photoresponse is assessed by analyzing the shapes of the obtained potential energy surfaces. A qualitative difference between facial ( fac) and meridional ( mer) isomers' potential energy surface (PES) topologies is evidenced for the first time in iron-based complexes. Indeed, the mer complex shows a steeper triplet path toward the corresponding ³MC minimum, which lies at a lower energy as compared to the fac isomer, thus pointing to a faster triplet decay of the former. Furthermore, while a major role of the metal-centered quintet state population from the triplet ³MC region is excluded, we identify the enlargement of iron-nitrogen bonds as the main normal modes driving the excited-state decay.

Tracking the picosecond deactivation dynamics of a photoexcited iron carbene complex by time-resolved X-ray scattering

Recent years have seen the development of new iron-centered N-heterocyclic carbene (NHC) complexes for solar energy applications. Compared to typical ligand systems, the NHC ligands provide Fe complexes with longer-lived metal-to-ligand charge transfer (MLCT) states. This increased lifetime is ascribed to strong ligand field splitting provided by the NHC ligands that raises the energy levels of the metal centered (MC) states and therefore reduces the deactivation efficiency of MLCT states. Among currently known NHC systems, [Fe(btbip)2]2+ (btbip = 2,6-bis(3-tert-butyl-imidazol-1-ylidene)pyridine) is a unique complex as it exhibits a short-lived MC state with a lifetime on the scale of a few hundreds of picoseconds. Hence, this complex allows for a detailed investigation, using 100 ps X-ray pulses from a synchrotron, of strong ligand field effects on the intermediate MC state in an NHC complex. Here, we use time-resolved wide angle X-ray scattering (TRWAXS) aided by density functional theory (DFT) to investigate the molecular structure, energetics and lifetime of the high-energy MC state in the Fe-NHC complex [Fe(btbip)2]2+ after excitation to the MLCT manifold. We identify it as a 260 ps metal-centered quintet (5MC) state, and we refine the molecular structure of the excited-state complex verifying the DFT results. Using information about the hydrodynamic state of the solvent, we also determine, for the first time, the energy of the 5MC state as 0.75 ± 0.15 eV. Our results demonstrate that due to the increased ligand field strength caused by NHC ligands, upon transition from the ground state to the 5MC state, the metal to ligand bonds extend by unusually large values: by 0.29 Å in the axial and 0.21 Å in the equatorial direction. These results imply that the transition in the photochemical properties from typical Fe complexes to novel NHC compounds is manifested not only in the destabilization of the MC states, but also in structural distortion of these states.

Bis-meridional Fe2+ spincrossover complexes of phenyl and pyridyl substituted 2-(pyridin-2-yl)-1,10-phenanthrolines

Fe²⁺ spincrossover complexes [Fe(L)₂]²⁺ (L = substituted (pyridin-2-yl)-1,10-phenanthroline) were prepared and SCO with changing coordination numbers was identified by ¹H NMR spectroscopy and in silico modeling.

Spin state variability in Fe2+ complexes of substituted (2-(pyridin-2-yl)-1,10-phenanthroline) ligands as versatile terpyridine analogues

Fe²⁺ spin crossover complexes [Fe(L)₂]²⁺ (L = substituted (pyridin-2-yl)-1,10-phenanthroline) were prepared and SCO properties were investigated in solution and in the solid state by an experiment and in silico.

Theoretical illumination of highly original photoreactive3MC states and the mechanism of the photochemistry of Ru(ii) tris(bidentate) complexes

Elucidation of the photoreactive mechanism of ruthenium(ii) complexes is reported along with identification of crucial and highly original metal-centred states.

Femtosecond X-Ray Scattering Study of Ultrafast Photoinduced Structural Dynamics in Solvated [Co(terpy)_{2}]^{2+}

We study the structural dynamics of photoexcited [Co(terpy)_{2}]^{2+} in an aqueous solution with ultrafast x-ray diffuse scattering experiments conducted at the Linac Coherent Light Source. Through direct comparisons with density functional theory calculations, our analysis shows that the photoexcitation event leads to elongation of the Co-N bonds, followed by coherent Co-N bond length oscillations arising from the impulsive excitation of a vibrational mode dominated by the symmetrical stretch of all six Co-N bonds. This mode has a period of 0.33 ps and decays on a subpicosecond time scale. We find that the equilibrium bond-elongated structure of the high spin state is established on a single-picosecond time scale and that this state has a lifetime of ∼7  ps.