Ironing out the photochemical and spin-crossover behavior of Fe(II) coordination compounds with computational chemistry

Impact of flexible hexyl chain ordering in a mononuclear spin crossover iron(III) complex

A novel FeIII complex [Fe(Hex-tnal)2]BPh4 (1) with a tridentate N2O ligand having an n-hexyl chain, Hex-Htnal (=1-((((1-hexyl-1H-1,2,3-triazol-4-yl)methyl)imino)methyl)naphthalen-2-ol), is reported. Temperature-dependent magnetic susceptibility measurements revealed that 1 exhibits a two-step spin crossover (SCO) transition in the 400-10 K temperature range, including an unusual gradual χMT change above RT (300-345 K) and a hysteretic χMT jump in a narrow temperature range of 345-357 K. These behaviors were also characterized by differential scanning calorimetry. Variable-temperature single-crystal X-ray diffraction studies revealed that the order-disorder transition and conformational change of the hexyl chains and the symmetry change associated with the re-entrant structural phase transition, namely triclinic P1̄ (100-275 K) ↔ monoclinic C2/c (296-340 K) ↔ triclinic P1̄ (360 K), are coupled to variations in intermolecular interactions and the N4O2 coordination environment, resulting in the occurrence of the unusual two-step SCO transition of 1. This study demonstrates that the flexible motion of alkyl substituents in the supramolecular lattice influences the occurrence of anomalous magnetic switching properties.

Controlling the Photophysical Properties of a Series of Isostructural d6 Complexes Based on Cr0, MnI, and FeII

Development of first-row transition metal complexes with similar luminescence and photoredox properties as widely used RuII polypyridines is attractive because metals from the first transition series are comparatively abundant and inexpensive. The weaker ligand field experienced by the valence d-electrons of first-row transition metals challenges the installation of the same types of metal-to-ligand charge transfer (MLCT) excited states as in precious metal complexes, due to rapid population of energetically lower-lying metal-centered (MC) states. In a family of isostructural tris(diisocyanide) complexes of the 3d6 metals Cr0, MnI, and FeII, the increasing effective nuclear charge and ligand field strength allow us to control the energetic order between the 3MLCT and 3MC states, whereas pyrene decoration of the isocyanide ligand framework provides control over intraligand (ILPyr) states. The chromium(0) complex shows red 3MLCT phosphorescence because all other excited states are higher in energy. In the manganese(I) complex, a microsecond-lived dark 3ILPyr state, reminiscent of the types of electronic states encountered in many polyaromatic hydrocarbon compounds, is the lowest and becomes photoactive. In the iron(II) complex, the lowest MLCT state has shifted to so much higher energy that 1ILPyr fluorescence occurs, in parallel to other excited-state deactivation pathways. Our combined synthetic-spectroscopic-theoretical study provides unprecedented insights into how effective nuclear charge, ligand field strength, and ligand π-conjugation affect the energetic order between MLCT and ligand-based excited states, and under what circumstances these individual states become luminescent and exploitable in photochemistry. Such insights are the key to further developments of luminescent and photoredox-active first-row transition metal complexes.

Iron(II) complexes of 2,6-bis(imidazo[1,2-a]pyridin-2-yl)pyridine and related ligands with annelated distal heterocyclic donors

Following a published synthesis of 2,6-bis(imidazo[1,2-a]pyridin-2-yl)pyridine (L1), treatment of α,α'-dibromo-2,6-diacetylpyridine with 2 equiv. 2-aminopyrimidine or 2-aminoquinoline in refluxing acetonitrile respectively gives 2,6-bis(imidazo[1,2-a]pyrimidin-2-yl)pyridine (L2) and 2,6-bis(imidazo[1,2-a]quinolin-2-yl)pyridine (L3). Solvated crystals of [Fe(L1)2][BF4]2 (1[BF4]2) and [Fe(L2)2][BF4]2 (2[BF4]2) are mostly high-spin, although one solvate of 1[BF4]2 undergoes thermal spin-crossover on cooling. The iron coordination geometry is consistently distorted in crystals of 2[BF4]2 which may reflect the influence of intramolecular, inter-ligand N⋯π interactions on the molecular conformation. Only 1 : 1 Fe : L3 complexes were observed in solution, or isolated in the solid state; a crystal structure of [FeBr(py)2L3]Br·0.5H2O (py = pyridine) is presented. A solvate crystal structure of high-spin [Fe(L4)2][BF4]2 (L4 = 2,6-di{quinolin-2-yl}pyridine; 4[BF4]2) is also described, which exhibits a highly distorted six-coordinate geometry with a helical ligand conformation. The iron(II) complexes are high-spin in solution at room temperature, but 1[BF4]2 and 2[BF4]2 undergo thermal spin-crossover equilibria on cooling. All the compounds exhibit a ligand-based emission in solution at room temperature. Gas phase DFT calculations mostly reproduce the spin state properties of the complexes, but show small anomalies attributed to intramolecular, inter-ligand dispersion interactions in the sterically crowded molecules.

Spin-state energetics and magnetic anisotropy in penta-coordinated Fe(III) complexes with different axial and equatorial ligand environments

The penta-coordinated trigonal-bi-pyramidal (TBP) Fe(III) complex (PMe₂Ph)₂FeCl₃ shows a reduced magnetic anisotropy in its intermediate-spin (IS) state as compared to its methyl-analog (PMe₃)₂Fe(III)Cl₃. In this work, the ligand environment in (PMe₂Ph)₂FeCl₃ is systematically altered by replacing the axial -P with -N and -As, the equatorial -Cl with other halides, and the axial methyl group with an acetyl group. This has resulted in a series of Fe(III) TBP complexes modelled in their IS and high-spin (HS) states. Lighter ligands -N and -F stabilize the complex in the HS state, while the magnetically anisotropic IS state is stabilized by -P and -As at the axial site, and -Cl, -Br, and -I at the equatorial site. Larger magnetic anisotropies appear for complexes with nearly degenerate ground electronic states that are well separated from the higher excited states. This requirement, largely controlled by the d-orbital splitting pattern due to the changing ligand field, is achieved with a certain combination of axial and equatorial ligands, such as -P and -Br, -As and -Br, and -As and -I. In most cases, the acetyl group at the axial site enhances the magnetic anisotropy compared to its methyl counterpart. In contrast, the presence of -I at the equatorial site compromises the uniaxial type of anisotropy of the Fe(III) complex leading to an enhanced rate of quantum tunneling of magnetization.

Non-Phenomenological Description of the Time-Resolved Emission in Solution with Quantum-Classical Vibronic Approaches-Application to Coumarin C153 in Methanol

We report a joint experimental and theoretical work on the steady-state spectroscopy and time-resolved emission of the coumarin C153 dye in methanol. The lowest energy excited state of this molecule is characterized by an intramolecular charge transfer thus leading to remarkable shifts of the time-resolved emission spectra, dictated by the methanol reorganization dynamics. We selected this system as a prototypical test case for the first application of a novel computational protocol aimed at the prediction of transient emission spectral shapes, including both vibronic and solvent effects, without applying any phenomenological broadening. It combines a recently developed quantum-classical approach, the adiabatic molecular dynamics generalized vertical Hessian method (Ad-MD|gVH), with nonequilibrium molecular dynamics simulations. For the steady-state spectra we show that the Ad-MD|gVH approach is able to reproduce quite accurately the spectral shapes and the Stokes shift, while a ∼0.15 eV error is found on the prediction of the solvent shift going from gas phase to methanol. The spectral shape of the time-resolved emission signals is, overall, well reproduced, although the simulated spectra are slightly too broad and asymmetric at low energies with respect to experiments. As far as the spectral shift is concerned, the calculated spectra from 4 ps to 100 ps are in excellent agreement with experiments, correctly predicting the end of the solvent reorganization after about 20 ps. On the other hand, before 4 ps solvent dynamics is predicted to be too fast in the simulations and, in the sub-ps timescale, the uncertainty due to the experimental time resolution (300 fs) makes the comparison less straightforward. Finally, analysis of the reorganization of the first solvation shell surrounding the excited solute, based on atomic radial distribution functions and orientational correlations, indicates a fast solvent response (≈100 fs) characterized by the strengthening of the carbonyl-methanol hydrogen bond interactions, followed by the solvent reorientation, occurring on the ps timescale, to maximize local dipolar interactions.

Photoactive Metal-to-Ligand Charge Transfer Excited States in 3d6 Complexes with Cr0, MnI, FeII, and CoIII

Many coordination complexes and organometallic compounds with the 4d⁶ and 5d⁶ valence electron configurations have outstanding photophysical and photochemical properties, which stem from metal-to-ligand charge transfer (MLCT) excited states. This substance class makes extensive use of the most precious and least abundant metal elements, and consequently there has been a long-standing interest in first-row transition metal compounds with photoactive MLCT states. Semiprecious copper(I) with its completely filled 3d subshell is a relatively straightforward and well explored case, but in 3d⁶ complexes the partially filled d-orbitals lead to energetically low-lying metal-centered (MC) states that can cause undesirably fast MLCT excited state deactivation. Herein, we discuss recent advances made with isoelectronic Cr⁰, Mn^I, Fe^II, and Co^III compounds, for which long-lived MLCT states have become accessible over the past five years. Furthermore, we discuss possible future developments in the search for new first-row transition metal complexes with partially filled 3d subshells and photoactive MLCT states for next-generation applications in photophysics and photochemistry.

Can range-separated functionals be optimally tuned to predict spectra and excited state dynamics in photoactive iron complexes?

Density functional theory is an efficient computational tool to investigate photophysical and photochemical processes in transition metal complexes, giving invaluable assistance in interpreting spectroscopic and catalytic experiments. Optimally tuned range-separated functionals are particularly promising, as they were created to address some of the fundamental deficiencies present in approximate exchange-correlation functionals. In this paper, we scrutinize the selection of optimally tuned parameters and its influence on the excited state dynamics, using the example of the iron complex [Fe(cpmp)₂]²⁺ with push-pull ligands. Various tuning strategies are contemplated based on pure self-consistent DFT protocols, as well as on the comparison with experimental spectra and multireference CASPT2 results. The two most promising sets of optimal parameters are then employed to carry out nonadiabatic surface-hopping dynamics simulations. Intriguingly, we find that the two sets lead to very different relaxation pathways and timescales. While the set of optimal parameters from one of the self-consistent DFT protocols predicts the formation of long-lived metal-to-ligand charge transfer triplet states, the set in better agreement with CASPT2 calculations leads to deactivation in the manifold of metal-centered states, in better agreement with the experimental reference data. These results showcase the complexity of iron-complex excited state landscapes and the difficulty of obtaining an unambiguous parametrization of long-range corrected functionals without experimental input.

Pseudo-Octahedral Iron(II) Complexes with Near-Degenerate Charge Transfer and Ligand Field States at the Franck-Condon Geometry

Increasing the metal-to-ligand charge transfer (MLCT) excited state lifetime of polypyridine iron(II) complexes can be achieved by lowering the ligand's π* orbital energy and by increasing the ligand field splitting. In the homo- and heteroleptic complexes [Fe(cpmp)2 ]2+ (12+ ) and [Fe(cpmp)(ddpd)]2+ (22+ ) with the tridentate ligands 6,2''-carboxypyridyl-2,2'-methylamine-pyridyl-pyridine (cpmp) and N,N'-dimethyl-N,N'-di-pyridin-2-ylpyridine-2,6-diamine (ddpd) two or one dipyridyl ketone moieties provide low energy π* acceptor orbitals. A good metal-ligand orbital overlap to increase the ligand field splitting is achieved by optimizing the octahedricity through CO and NMe units between the coordinating pyridines which enable the formation of six-membered chelate rings. The push-pull ligand cpmp provides intra-ligand and ligand-to-ligand charge transfer (ILCT, LL'CT) excited states in addition to MLCT excited states. Ground and excited state properties of 12+ and 22+ were accessed by X-ray diffraction analyses, resonance Raman spectroscopy, (spectro)electrochemistry, EPR spectroscopy, X-ray emission spectroscopy, static and time-resolved IR and UV/Vis/NIR absorption spectroscopy as well as quantum chemical calculations.

Environmental Control of the Magnetic Behavior of Transition Metal Complexes: Density Functional Theory Study of Zeolite Y Embedded Complexes [M(bpy)3]2+@Y (M = Fe2+, Co2+)

Using the supramolecular approach developed for the study of the guest-host interactions in the zeolite Y encapsulated [Fe(bpy)₃]²⁺ compound: [Fe(bpy)₃]²⁺@Y (bpy = 2,2'-bipyridine) [Vargas et al., J. Chem. Theory Comput. 2009, 5, 97-115], we apply density functional theory (DFT) to the study of the influence of zeolite Y encapsulation on the structural and energetic properties of [Co(bpy)₃]²⁺ in the low-spin (LS) and high-spin (HS) states, while revisiting [Fe(bpy)₃]²⁺@Y. Although the accurate prediction of the HS-LS energy difference ΔE_HL^el remains challenging for current DFT methods, they give accurate estimates of its variation Δ(ΔE_HL^el) in a series of complexes of a given transition metal ion. Therefore, denoting [M(bpy)₃]²⁺@Y_SM as the supramolecular model of the inclusion compounds, the values of ΔE_HL^el for the bpy complexes in the gas phase and in the supercage of zeolite Y were determined by combining the DFT estimates of Δ(ΔE_HL^el) in the series {[M(NCH)₆]²⁺, [M(bpy)₃]²⁺, and [M(bpy)₃]²⁺@Y_SM}, with accurate CCSD(T) estimates of ΔE_HL^el in the benchmark complexes [M(NCH)₆]²⁺ (M = Fe, Co) [Lawson Daku et al., J. Chem. Theory Comput., 2012, 8, 4216-4231]. Generalized gradient approximations as well as global and range-separated hybrids were employed. In order to better account for the key role of dispersion, they were also augmented with the semiempirical D2, D3BJ, and D3BJM dispersion corrections when available. The use of the D3BJ and D3BJM corrections led to similar results, and this is only with the use of the D2 scheme that (i) the free and encapsulated [Fe(bpy)₃]²⁺ are correctly predicted as LS species and that (ii) the encapsulation of both complexes translates into a destabilization of their HS state with respect to their LS state. The increase of the HS-LS energy difference is smaller for [Co(bpy)₃]²⁺ than [Fe(bpy)₃]²⁺ because the HS-LS molecular volume difference ΔV_HL in [Co(bpy)₃]²⁺ is ∼50% smaller than in [Fe(bpy)₃]²⁺. Periodic DFT calculations performed on crystalline [M(bpy)₃]²⁺@Y show that the employed [M(bpy)₃]²⁺@Y_SM supramolecular model allows the influence of encapsulation on the geometry and the spin-state energetics of [M(bpy)₃]²⁺ (M = Fe, Co) to be quantitatively captured.

Adiabatic models for the quantum dynamics of surface scattering with lattice effects

In this contribution, we review models for the lattice effects in quantum dynamics calculations on surface scattering, which is important to modeling heterogeneous catalysis for achieving an interpretation of experimental measurements. Unlike dynamics models for reactions in the gas phase, those for heterogeneous reactions have to include the effects of the surface. For manageable computational costs in calculations, the effects of static surface (SS) are firstly modeled as this is simply and easily implemented. Then, the SS model has to be improved to include the effects of the flexible surface, that is the lattice effects. To do this, various surface models have been designed where the coordinates of the surface atoms are introduced in the Hamiltonian operator, especially those of the top surface atom. Based on this model Hamiltonian operator, extensive multi-dimension quantum dynamics calculations can be performed to recover the lattice effects. Here, we first review an overview of the techniques in constructing the Hamiltonian operator, which is a sum of the kinetic energy operator (KEO) and potential energy surface (PES). Since the PES containing the coordinates of the surface atoms in a cell is still expensive, the SS model is often accepted. We consider a mathematical model, called the coupled harmonic oscillator (CHO) model, to introduce the concepts of adiabatic and diabatic representations for separating the molecule and surface. Under the adiabatic model, we further introduce the expansion model where the potential function is Taylor expanded around the optimized geometry of the surface. By an expansion model truncated at the first and second order, various coupling surface models between the molecule and surface are derived. Moreover, by further and deeply understanding the adiabatic representation, an effective Hamiltonian operator is obtained by optimizing the total wave function in factorized form. By this factorized form of wave function and effective Hamiltonian operator, the geometry phase of the surface wave function is theoretically found. This theoretical prediction may be measured by carefully designing experiments. Finally, discussions on the adiabatic representation, the PES construction, and possibility of the classical-dynamics solutions are given. Based on these discussions, a simple outlook on the dynamics of photocatalytics is finally given.

Reconciling Local Coupled Cluster with Multireference Approaches for Transition Metal Spin-State Energetics

Spin-state energetics of transition metal complexes remain one of the most challenging targets for electronic structure methods. Among single-reference wave function approaches, local correlation approximations to coupled cluster theory, most notably the domain-based local pair natural orbital (DLPNO) approach, hold the promise of bringing the accuracy of coupled cluster theory with single, double, and perturbative triple excitations, CCSD(T), to molecular systems of realistic size with acceptable computational cost. However, recent studies on spin-state energetics of iron-containing systems raised doubts about the ability of the DLPNO approach to adequately and systematically approximate energetics obtained by the reference-quality complete active space second-order perturbation theory with coupled-cluster semicore correlation, CASPT2/CC. Here, we revisit this problem using a diverse set of iron complexes and examine several aspects of the application of the DLPNO approach. We show that DLPNO-CCSD(T) can accurately reproduce both CASPT2/CC and canonical CCSD(T) results if two basic principles are followed. These include the consistent use of the improved iterative (T₁) versus the semicanonical perturbative triple corrections and, most importantly, a simple two-point extrapolation to the PNO space limit. The latter practically eliminates errors arising from the default truncation of electron-pair correlation spaces and should be viewed as standard practice in applications of the method to transition metal spin-state energetics. Our results show that reference-quality results can be readily achieved with DLPNO-CCSD(T) if these principles are followed. This is important also in view of the applicability of the method to larger single-reference systems and multinuclear clusters, whose treatment of dynamic correlation would be challenging for multireference-based approaches.

Machine learning of material properties: Predictive and interpretable multilinear models

Machine learning models can provide fast and accurate predictions of material properties but often lack transparency. Interpretability techniques can be used with black box solutions, or alternatively, models can be created that are directly interpretable. We revisit material datasets used in several works and demonstrate that simple linear combinations of nonlinear basis functions can be created, which have comparable accuracy to the kernel and neural network approaches originally used. Linear solutions can accurately predict the bandgap and formation energy of transparent conducting oxides, the spin states for transition metal complexes, and the formation energy for elpasolite structures. We demonstrate how linear solutions can provide interpretable predictive models and highlight the new insights that can be found when a model can be directly understood from its coefficients and functional form. Furthermore, we discuss how to recognize when intrinsically interpretable solutions may be the best route to interpretability.

CASPT2 molecular geometries of Fe(II) spin-crossover complexes

Using fully internally contracted (FIC)-CASPT2 analytical gradients, geometry optimizations of spin-crossover complexes are reported. This approach is tested on a series of Fe(II) complexes with different sizes, ranging from 13 to 61 atoms. A combination of active space and basis set choices are employed to investigate their role in determining reliable molecular geometries. The reported strategy demonstrates that a wave function-based level of theory can be used to optimize the geometries of metal complexes in reasonable times and enables one to treat the molecular geometry and electronic structure of the complexes using the same level of theory. For a series of smaller Fe(II) SCO complexes, strong field ligands in the LS state result in geometries with the largest differences between DFT and CASPT2; however, good agreement overall is observed between DFT and CASPT2. For the larger complexes, moderate sized basis sets yield geometries that compare well with DFT and available experimental data. We recommend using the (10e,12o) active space since convergence to a minimum structure was more efficient than with truncated active spaces despite having similar Fe-ligand bond distances.

Heteroleptic iron(ii) complexes of chiral 2,6-bis(oxazolin-2-yl)-pyridine (PyBox) and 2,6-bis(thiazolin-2-yl)pyridine ligands – the interplay of two different ligands on the metal ion spin sate

The spin-crossover properties of [Fe(LR)L][ClO4]2 (LR = a chiral PyBox {L1R} or ThioPyBox {L2R} derivative) show subtle differences depending on the tridentate ‘L’ co-ligand.

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.

Luminescent First-Row Transition Metal Complexes

Precious and rare elements have traditionally dominated inorganic photophysics and photochemistry, but now we are witnessing a paradigm shift toward cheaper and more abundant metals. Even though emissive complexes based on selected first-row transition metals have long been known, recent conceptual breakthroughs revealed that a much broader range of elements in different oxidation states are useable for this purpose. Coordination compounds of V, Cr, Mn, Fe, Co, Ni, and Cu now show electronically excited states with unexpected reactivity and photoluminescence behavior. Aside from providing a compact survey of the recent conceptual key advances in this dynamic field, our Perspective identifies the main design strategies that enabled the discovery of fundamentally new types of 3d-metal-based luminophores and photosensitizers operating in solution at room temperature.

Fate of Sc-Ion Interaction With Water: A Computational Study to Address Splitting Water Versus Solvating Sc Ion

An exhaustive study of Sc-ion interaction with water molecules in all its possible oxidation and spin states has been carried out to delineate the relative propensity of Sc ions toward solvation and water splitting. Potential energy surface analysis of the Sc-ion reaction with water molecules, topological analysis of bonds, and the effect of sequential solvation up to 6 water molecules have been examined. Calculated values showed good agreement with the available experimental results. Close-shell systems such as singlet mono- and tricationic Sc ions prefer to split the water molecules. In contrast, the open-shell systems such as triplet mono- and doublet dicationic Sc ions prefer to get solvated than split the water molecule. Topological analysis of electron density predicted the Sc^+/2+–water bond as a noncovalent bond while Sc³⁺–OH₂, Sc²⁺–OH, and Sc⁺–H bonds as partially covalent in nature. Energy decomposition analysis revealed that Sc ion–water interactions are driven by electrostatic energy followed by polarization energy. The current study reveals that transition metal catalysis can be one of the most effective tools to employ in water splitting, by properly tuning the electrons, spin, and ligands around the catalytic center.

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.

Deciphering Cryptic Behavior in Bimetallic Transition-Metal Complexes with Machine Learning

We demonstrate an alternative, data-driven approach to uncovering structure-property relationships for the rational design of heterobimetallic transition-metal complexes that exhibit metal-metal bonding. We tailor graph-based representations of the metal-local environment for these complexes for use in multiple linear regression and kernel ridge regression (KRR) models. We curate a set of 28 experimentally characterized complexes to develop a multiple linear regression model for oxidation potentials. We achieve good accuracy (mean absolute error of 0.25 V) and preserve transferability to unseen experimental data with a new ligand structure. We also train a KRR model on a subset of 330 structurally characterized heterobimetallics to predict the degree of metal-metal bonding. This KRR model predicts relative metal-metal bond lengths in the test set to within 5%, and analysis of key features reveals the fundamental atomic contributions (e.g., the valence electron configuration) that most strongly influence the behavior of these complexes. Our work provides guidance for rational bimetallic design, suggesting that properties, including the formal shortness ratio, should be transferable from one period to another.

Structures and Spin States of Iron(II) Complexes of Isomeric 2,6-Di(1,2,3-triazolyl)pyridine Ligands

Iron(II) complex salts of 2,6-di(1,2,3-triazol-1-yl)pyridine (L¹) are unexpectedly unstable in undried solvent. This is explained by the isolation of [Fe(L¹)₄(H₂O)₂][ClO₄]₂ and [Fe(NCS)₂(L¹)₂(H₂O)₂]·L¹, containing L¹ bound as a monodentate ligand rather than in the expected tridentate fashion. These complexes associate into 4⁴ grid structures through O-H···N hydrogen bonding; a solvate of a related 4⁴ coordination framework, catena-[Cu(μ-L¹)₂(H₂O)₂][BF₄]₂, is also presented. The isomeric ligands 2,6-di(1,2,3-triazol-2-yl)pyridine (L²) and 2,6-di(1H-1,2,3-triazol-4-yl)pyridine (L³) bind to iron(II) in a more typical tridentate fashion. Solvates of [Fe(L³)₂][ClO₄]₂ are low-spin and diamagnetic in the solid state and in solution, while [Fe(L²)₂][ClO₄]₂ and [Co(L³)₂][BF₄]₂ are fully high-spin. Treatment of L³ with methyl iodide affords 2,6-di(2-methyl-1,2,3-triazol-4-yl)pyridine (L⁴) and 2-(1-methyl-1,2,3-triazol-4-yl)-6-(2-methyl-1,2,3-triazol-4-yl)pyridine (L⁵). While salts of [Fe(L⁵)₂]²⁺ are low-spin in the solid state, [Fe(L⁴)₂][ClO₄]₂·H₂O is high-spin, and [Fe(L⁴)₂][ClO₄]₂·3MeNO₂ exhibits a hysteretic spin transition to 50% completeness at T_1/2 = 128 K (ΔT_1/2 = 6 K). This transition proceeds via a symmetry-breaking phase transition to an unusual low-temperature phase containing three unique cation sites with high-spin, low-spin, and 1:1 mixed-spin populations. The unusual distribution of the spin states in the low-temperature phase reflects "spin-state frustration" of the mixed-spin cation site by an equal number of high-spin and low-spin nearest neighbors. Gas-phase density functional theory calculations reproduce the spin-state preferences of these and some related complexes. These highlight the interplay between the σ-basicity and π-acidity of the heterocyclic donors in this ligand type, which have opposing influences on the molecular ligand field. The Brønsted basicities of L¹-L³ are very sensitive to the linkage isomerism of their triazolyl donors, which explains why their iron complex spin states show more variation than the better-known iron(II)/2,6-dipyrazolylpyridine system.

Spin-States of Diastereomeric Iron(II) Complexes of 2,6-Bis(thiazolin-2-yl)pyridine (ThioPyBox) Ligands and a Comparison with the Corresponding PyBox Derivatives

This report investigates homoleptic iron(II) complexes of thiazolinyl analogues of chiral PyBox tridentate ligands: 2,6-bis(4-phenyl-4,5-dihydrothiazol-2-yl)pyridine (L¹Ph), 2,6-bis(4-isopropyl-4,5-dihydrothiazol-2-yl)pyridine (L¹iPr), and 2,6-bis(4-tert-butyl-4,5-dihydrothiazol-2-yl)pyridine (L¹t-Bu). Crystallographic data imply the larger and more flexible thiazolinyl rings reduce steric clashes between the R substituents in homochiral [Fe((R)-L¹R)₂]²⁺ or [Fe((S)-L¹R)₂]²⁺ (R = Ph, iPr, or t-Bu), compared to their PyBox (L²R) analogues. Conversely, the larger heterocyclic S atoms are in close contact with the R substituents in heterochiral [Fe((R)-L¹Ph)((S)-L¹Ph)]²⁺, giving it a more sterically hindered ligand environment than that in [Fe((R)-L²Ph)((S)-L²Ph)]²⁺ (L²Ph = 2,6-bis(4-phenyl-4,5-dihydrooxazol-2-yl)pyridine). Preformed [Fe((R)-L¹Ph)((S)-L¹Ph)]²⁺ and [Fe((R)-L¹iPr)((S)-L¹iPr)]²⁺ do not racemize by ligand redistribution in CD₃CN solution, but homochiral [Fe(L¹iPr)₂]²⁺ and [Fe(L¹t-Bu)₂]²⁺ both undergo partial ligand displacement in that solvent. Homochiral [Fe(L¹Ph)₂]²⁺ and [Fe(L¹iPr)₂]²⁺ exhibit spin-crossover equilibria in CD₃CN, centered at 344 ± 6 K and 277 ± 1 K respectively, while their heterochiral congeners are essentially low-spin within the liquid range of the solvent. These data imply that the diastereomers of [Fe(L¹Ph)₂]²⁺ and [Fe(L¹iPr)₂]²⁺ show a greater difference in their spin-state behaviors than was previous found for [Fe(L²Ph)₂]²⁺. Gas-phase DFT calculations (B86PW91/def2-SVP) of the [Fe(L¹R)₂]²⁺ and [Fe(L²R)₂]²⁺ complexes reproduce most of the observed trends, but they overstabilize the high-spin state of SCO-active [Fe(L¹iPr)₂]²⁺ by ca. 1.5 kcal mol^-1. This might reflect the influence of intramolecular dispersion interactions on the spin states of these compounds. Attempts to model this with the dispersion-corrected functionals B97-D2 or PBE-D3 were less successful than our original protocol, confirming that the spin states of sterically hindered molecules are a challenging computational problem.

Capturing photochemical and photophysical transformations in iron complexes with ultrafast X-ray spectroscopy and scattering

Light-driven chemical transformations provide a compelling approach to understanding chemical reactivity with the potential to use this understanding to advance solar energy and catalysis applications. Capturing the non-equilibrium trajectories of electronic excited states with precision, particularly for transition metal complexes, would provide a foundation for advancing both of these objectives. Of particular importance for 3d metal compounds is characterizing the population dynamics of charge-transfer (CT) and metal-centered (MC) electronic excited states and understanding how the inner coordination sphere structural dynamics mediate the interaction between these states. Recent advances in ultrafast X-ray laser science has enabled the electronic excited state dynamics in 3d metal complexes to be followed with unprecedented detail. This review will focus on simultaneous X-ray emission spectroscopy (XES) and X-ray solution scattering (XSS) studies of iron coordination and organometallic complexes. These simultaneous XES-XSS studies have provided detailed insight into the mechanism of light-induced spin crossover in iron coordination compounds, the interaction of CT and MC excited states in iron carbene photosensitizers, and the mechanism of Fe-S bond dissociation in cytochrome c.

Building an Iron Chromophore Incorporating Prussian Blue Analogue for Photoelectrochemical Water Oxidation

The replacement of traditional ruthenium-based photosensitizers with low-cost and abundant iron analogs is a key step for the advancement of scalable and sustainable dye-sensitized water splitting cells. In this proof-of-concept study, a pyridinium ligand coordinated pentacyanoferrate(II) chromophore is used to construct a cyanide-based CoFe extended bulk framework, in which the iron photosensitizer units are connected to cobalt water oxidation catalytic sites through cyanide linkers. The iron-sensitized photoanode exhibits exceptional stability for at least 5 h at pH 7 and features its photosensitizing ability with an incident photon-to-current conversion capacity up to 500 nm with nanosecond scale excited state lifetime. Ultrafast transient absorption and computational studies reveal that iron and cobalt sites mutually support each other for charge separation via short bridging cyanide groups and for injection to the semiconductor in our proof-of-concept photoelectrochemical device. The reorganization of the excited states due to the mixing of electronic states of metal-based orbitals subsequently tailor the electron transfer cascade during the photoelectrochemical process. This breakthrough in chromophore-catalyst assemblies will spark interest in dye-sensitization with robust bulk systems for photoconversion applications.

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.

A Ranked-Orbital Approach to Select Active Spaces for High-Throughput Multireference Computation

The past decade has seen a great increase in the application of high-throughput computation to a variety of important problems in chemistry. However, one area which has been resistant to the high-throughput approach is multireference wave function methods, in large part due to the technicalities of setting up these calculations and in particular the not always intuitive challenge of active space selection. As we look toward a future of applying high-throughput computation to all areas of chemistry, it is important to prepare these methods for large-scale automation. Here, we propose a ranked-orbital approach to select active spaces with the goal of standardizing multireference methods for high-throughput computation. This method allows for the meaningful comparison of different active space selection schemes and orbital localizations, and we demonstrate the utility of this approach across 1120 multireference calculations for the excitation energies of small molecules. Our results reveal that it is helpful to distinguish the method used to generate orbitals from the method of ranking orbitals in terms of importance for the active space. Additionally, we propose our own orbital ranking scheme that estimates the importance of an orbital for the active space through a pair-interaction framework from orbital energies and features of the Hartree-Fock exchange matrix. We call this new scheme the "approximate pair coefficient" (APC) method and we show that it performs quite well for the test systems presented.

Spin State Tunes Oxygen Atom Transfer towards FeIV O Formation in FeII Complexes

Oxoiron(IV) complexes bearing tetradentate ligands have been extensively studied as models for the active oxidants in non-heme iron-dependent enzymes. These species are commonly generated by oxidation of their ferrous precursors. The mechanisms of these reactions have seldom been investigated. In this work, the reaction kinetics of complexes [Fe^II (CH₃ CN)₂ L](SbF₆ )₂ ([1](SbF₆ )₂ and [2](SbF₆ )₂ ) and [Fe^II (CF₃ SO₃ )₂ L] ([1](OTf)₂ and [2](OTf)₂ (1, L=^Me,H Pytacn; 2, L=^nP,H Pytacn; ^R,R' Pytacn=1-[(6-R'-2-pyridyl)methyl]-4,7- di-R-1,4,7-triazacyclononane) with Bu₄ NIO₄ to form the corresponding [Fe^IV (O)(CH₃ CN)L]²⁺ (3, L=^Me,H Pytacn; 4, L=^nP,H Pytacn) species was studied in acetonitrile/acetone at low temperatures. The reactions occur in a single kinetic step with activation parameters independent of the nature of the anion and similar to those obtained for the substitution reaction with Cl^- as entering ligand, which indicates that formation of [Fe^IV (O)(CH₃ CN)L]²⁺ is kinetically controlled by substitution in the starting complex to form [Fe^II (IO₄ )(CH₃ CN)L]⁺ intermediates that are converted rapidly to oxo complexes 3 and 4. The kinetics of the reaction is strongly dependent on the spin state of the starting complex. A detailed analysis of the magnetic susceptibility and kinetic data for the triflate complexes reveals that the experimental values of the activation parameters for both complexes are the result of partial compensation of the contributions from the thermodynamic parameters for the spin-crossover equilibrium and the activation parameters for substitution. The observation of these opposite and compensating effects by modifying the steric hindrance at the ligand illustrates so far unconsidered factors governing the mechanism of oxygen atom transfer leading to high-valent iron oxo species.

Theoretical modeling of the singlet-triplet spin transition in different Ni(II)-diketo-pyrphyrin-based metal-ligand octahedral complexes

The structural stability, charge transfer effects and strength of the spin-orbit couplings in different Ni(ii)-ligand complexes have been studied at the DFT (B3LYP and CAM-B3LYP) and coupled cluster (DLPNO-CCSD(T)) levels of theory. Accordingly, two different, porphyrin- and diketo-pyrphyrin-based four-coordination macrocycles as planar ligands as well as pyridine (or pyrrole) and mesylate anion molecular groups as vertical ligands were considered in order to build metal-organic complexes with octahedral coordination configurations. For each molecular system, the identification of equilibrium geometries and the intersystem crossing (the minimum energy crossing) points between the potential energy surfaces of the singlet and triplet spin states is followed by computing the spin-orbit couplings between the two spin states. Structures, based on the diketo-pyrphyrin macrocycle as the planar ligand, show stronger six-coordination metal-organic complexes due to the extra electrostatic interaction between the positively charged central metal cation and the negatively charged vertical ligands. The results also show that the magnitude of the spin-orbit coupling is influenced by the atomic positions of deprotonations of the ligands, and implicitly the direction of the charge transfer between the ligand and the central metal ion.

A Photorobust Mo(0) Complex Mimicking [Os(2,2'-bipyridine)3]2+ and Its Application in Red-to-Blue Upconversion

Osmium(II) polypyridines are a well-known class of complexes with luminescent metal-to-ligand charge-transfer (MLCT) excited states that are currently experiencing a revival due to their application potential in organic photoredox catalysis, triplet-triplet annihilation upconversion, and phototherapy. At the same time, there is increased interest in the development of photoactive complexes made from Earth-abundant rather than precious metals. Against this background, we present a homoleptic Mo(0) complex with a new diisocyanide ligand exhibiting different bite angles and a greater extent of π-conjugation than previously reported related chelates. This new design leads to deep red emission, which is unprecedented for homoleptic arylisocyanide complexes of group 6 metals. With a ³MLCT lifetime of 56 ns, an emission band maximum at 720 nm, and a photoluminescence quantum yield of 1.5% in deaerated toluene at room temperature, the photophysical properties are reminiscent of the prototypical [Os(2,2'-bipyridine)₃]²⁺ complex. Under 635 nm irradiation with a cw-laser, the new Mo(0) complex sensitizes triplet-triplet annihilation upconversion of 9,10-diphenylanthracene (DPA), resulting in delayed blue fluorescence with an anti-Stokes shift of 0.93 eV. The photorobustness of the Mo(0) complex and the upconversion quantum yield are high enough to generate a flux of upconverted light that can serve as a sufficiently potent irradiation source for a blue-light-driven photoisomerization reaction. These findings are relevant in the greater contexts of designing new luminophores and photosensitizers for use in red-light-driven photocatalysis, photochemical upconversion, light-harvesting, and phototherapy.

Ultrafast processes: coordination chemistry and quantum theory

Coordination compounds, characterized by fascinating and tunable electronic properties, are capable of binding easily to proteins, polymers, wires and DNA. Upon irradiation, these molecular systems develop functions finding applications in solar cells, photocatalysis, luminescent and conformational probes, electron transfer triggers and diagnostic or therapeutic tools. The control of these functions is activated by the light wavelength, the metal/ligand cooperation and the environment within the first picoseconds (ps). After a brief summary of the theoretical background, this perspective reviews case studies, from 1st row to 3rd row transition metal complexes, that illustrate how spin-orbit, vibronic coupling and quantum effects drive the photophysics of this class of molecules at the early stage of the photoinduced elementary processes within the fs-ps time scale range.

On the Spin-State Dependence of Redox Potentials of Spin Crossover Complexes

Resistance switching properties of nanoscale junctions of spin crossover molecules have received recently much interest. In many cases, this property has been traced back to the variation of molecular orbital energies upon spin transition. However, one can also expect a substantial reorganization of the molecular structure due to charge localization, which calls for a better understanding of the relationship between the redox potential and the spin state of the molecule. To investigate this issue, we carried out a detailed density functional theory and variable temperature cyclic voltammetry investigation of the benchmark compound [Fe(HB(1,2,4-triazol-1-yl)₃)₂] in solution. We show that, for a correct thermodynamical picture, it is necessary to take into account the charge transfer-induced electronic and structural reorganization as well as spin equilibria in the oxidized and reduced species.

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.

Rapid Detection of Strong Correlation with Machine Learning for Transition-Metal Complex High-Throughput Screening

Despite its widespread use in chemical discovery, approximate density functional theory (DFT) is poorly suited to many targets, such as those containing open-shell, 3d transition metals that can be expected to have strong multireference (MR) character. For discovery workflows to be predictive, we need automated, low-cost methods that can distinguish the regions of chemical space where DFT should be applied from those where it should not. We curate more than 4800 open-shell transition-metal complexes up to hundreds of atoms in size from prior high-throughput DFT studies and evaluate affordable, finite-temperature DFT fractional occupation number (FON)-based MR diagnostics. We show that intuitive measures of strong correlation (i.e., the HOMO-LUMO gap) are not predictive of MR character as judged by FON-based diagnostics. Analysis of independently trained machine learning (ML) models to predict HOMO-LUMO gaps and FON-based diagnostics reveals differences in the metal and ligand sensitivity of the two quantities. We use our trained ML models to rapidly evaluate MR character over a space of ∼187000 theoretical complexes, identifying large-scale trends in spin-state-dependent MR character and finding small HOMO-LUMO gap complexes while ensuring low MR character.

Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory

Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe²⁺, Fe³⁺, and Co³⁺, including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. We find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method-in terms of computational cost and in terms of accuracy on monatomic systems-and therefore we chose this method for molecular benchmarks. Then Kohn-Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF₂, FeCl₂, and CoF₂. We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.

Detailed Pair Natural Orbital-Based Coupled Cluster Studies of Spin Crossover Energetics

In this work, a detailed study of spin-state splittings in three spin crossover model compounds with DLPNO-CCSD(T) is presented. The performance in comparison to canonical CCSD(T) is assessed in detail. It was found that spin-state splittings with chemical accuracy, compared to the canonical results, are achieved when the full iterative triples (T1) scheme and TightPNO settings are applied and relativistic effects are taken into account. Having established the level of accuracy that can be reached relative to the canonical results, we have undertaken a detailed basis set study in the second part of the study. The slow convergence of the results of correlated calculations with respect to basis set extension is particularly acute for spin-state splittings for reasons discussed in detail in this Article. In fact, for some of the studied systems, 5Z basis sets are necessary in order to come close to the basis set limit that is estimated here by basis set extrapolation. Finally, the results of the present work are compared to available literature. In general, acceptable agreement with previous CCSD(T) results is found, although notable deviations stemming from differences in methodology and basis sets are noted. It is noted that the published CASPT2 numbers are far away from the extrapolated CCSD(T) numbers. In addition, dynamic quantum Monte Carlo results differ by several tens of kcal/mol from the CCSD(T) numbers. A comparison to DFT results produced with a range of popular density functionals shows the expected scattering of results and showcases the difficulty of applying DFT to spin-state energies.

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.

Large-scale comparison of 3d and 4d transition metal complexes illuminates the reduced effect of exchange on second-row spin-state energetics

The origin of distinct 3d vs. 4d transition metal complex sensitivity to exchange is explored over a large data set.

Approximate Analytical Gradients and Nonadiabatic Couplings for the State-Average Density Matrix Renormalization Group Self-Consistent-Field Method

We present an approximate scheme for analytical gradients and nonadiabatic couplings for calculating state-average density matrix renormalization group self-consistent-field wave function. Our formalism follows closely the state-average complete active space self-consistent-field (SA-CASSCF) ansatz, which employs a Lagrangian, and the corresponding Lagrange multipliers are obtained from a solution of the coupled-perturbed CASSCF (CP-CASSCF) equations. We introduce a definition of the matrix product state (MPS) Lagrange multipliers based on a single-site tensor in a mixed-canonical form of the MPS, such that a sweep procedure is avoided in the solution of the CP-CASSCF equations. We apply our implementation to the optimization of a conical intersection in 1,2-dioxetanone, where we are able to fully reproduce the SA-CASSCF result up to arbitrary accuracy.

Iron(II) coordination complexes with panchromatic absorption and nanosecond charge-transfer excited state lifetimes

Replacing current benchmark rare-element photosensitizers with ones based on abundant and low-cost metals such as iron would help facilitate the large-scale implementation of solar energy conversion. To do so, the ability to extend the lifetimes of photogenerated excited states of iron complexes is critical. Here, we present a sensitizer design in which iron(II) centres are supported by frameworks containing benzannulated phenanthridine and quinoline heterocycles paired with amido donors. These complexes exhibit panchromatic absorption and nanosecond charge-transfer excited state lifetimes, enabled by the combination of vacant, energetically accessible heterocycle-based acceptor orbitals and occupied molecular orbitals destabilized by strong mixing between amido nitrogen atoms and iron. This finding shows how ligand design can extend metal-to-ligand charge-transfer-type excited state lifetimes of iron(II) complexes into the nanosecond regime and expand the range of potential applications for iron-based photosensitizers.

Long-Lived, Strongly Emissive, and Highly Reducing Excited States in Mo(0) Complexes with Chelating Isocyanides

Newly discovered tris(diisocyanide)molybdenum(0) complexes are Earth-abundant isoelectronic analogues of the well-known class of [Ru(α-diimine)₃]²⁺ compounds with long-lived ³MLCT (metal-to-ligand charge transfer) excited states that lead to rich photophysics and photochemistry. Depending on ligand design, luminescence quantum yields up to 0.20 and microsecond excited state lifetimes are achieved in solution at room temperature, both significantly better than those for [Ru(2,2'-bipyridine)₃]²⁺. The excited Mo(0) complexes can induce chemical reactions that are thermodynamically too demanding for common precious metal-based photosensitizers, including the widely employed fac-[Ir(2-phenylpyridine)₃] complex, as demonstrated on a series of light-driven aryl-aryl coupling reactions. The most robust Mo(0) complex exhibits stable photoluminescence and remains photoactive after continuous irradiation exceeding 2 months. Our comprehensive optical spectroscopic and photochemical study shows that Mo(0) complexes with diisocyanide chelate ligands constitute a new family of luminophores and photosensitizers, which is complementary to precious metal-based 4d⁶ and 5d⁶ complexes and represents an alternative to nonemissive Fe(II) compounds. This is relevant in the greater context of sustainable photophysics and photochemistry, as well as for possible applications in lighting, sensing, and catalysis.

DFT Study of the CNS Ligand Effect on the Geometry, Spin-State, and Absorption Spectrum in Ruthenium, Iron, and Cobalt Quaterpyridine Complexes

Geometry parameters, total energy of the system in different spin states, harmonic vibrational frequencies, and absorption spectra were computed for a range of mononuclear quaterpyridine Ru(II), Fe(III/II), and Co(III/II) complexes with two axial ambidentate CNS ligands by using density functional theory (DFT) and time-dependent DFT calculations. Both structural and electronic properties were found to be correlating with the type of the binding atom in the CNS ligand (isomerization differs by 4-13 kcal·mol-1). The N-bonding of CNS ligands is energetically favored. It was also found that the low spin (LS) state is the ground state for both Ru(II) and Co(III) complexes regardless of the CNS arrangement. The other complexes are the high-spin (HS) ground-state ones with the only exception of the S-bonded CNS isomer of the Fe(III) complex. The dependencies of energy differences between the HS and LS states versus C demonstrated stabilization of the HS state with an increasing amount of the exact exchange admixture (C) for iron and cobalt complexes. An opposite behavior was observed for ruthenium complexes. The best match in harmonic vibrational frequencies between the experimental and calculated values has been reached at C = 0.15 for all the complexes. The absorption profile of the Fe(II) complex with the alternatively bonded CNS ligands strongly depends on the angle between them. The light-harvesting efficiency of the Fe(II) complexes is very similar (∼0.4) and sufficiently close to that of the Ru(II) complexes. The iron-based coordination compounds are considered as a prospective dye for dye-sensitized solar cells. The results of calculations were completed with experimental reference data, thus providing a systematic compendium for practical use.