Theoretical Investigation of the Electronic Structure of Fe(II) Complexes at Spin-State Transitions

Femtosecond Core-Level Spectroscopy Reveals Involvement of Triplet States in the Gas-Phase Photodissociation of Fe(CO)5

Excitation of iron pentacarbonyl [Fe(CO)5], a prototypical photocatalyst, at 266 nm causes the sequential loss of two CO ligands in the gas phase, creating catalytically active, unsaturated iron carbonyls. Despite numerous studies, major aspects of its ultrafast photochemistry remain unresolved because the early excited-state dynamics have so far eluded spectroscopic observation. This has led to the long-held assumption that ultrafast dissociation of gas-phase Fe(CO)5 proceeds exclusively on the singlet manifold. Herein, we present a combined experimental-theoretical study employing ultrafast extreme ultraviolet transient absorption spectroscopy near the Fe M2,3-edge, which features spectral evolution on 100 fs and 3 ps time scales, alongside high-level electronic structure theory, which enables characterization of the molecular geometries and electronic states involved in the ultrafast photodissociation of Fe(CO)5. We assign the 100 fs evolution to spectroscopic signatures associated with intertwined structural and electronic dynamics on the singlet metal-centered states during the first CO loss and the 3 ps evolution to the competing dissociation of Fe(CO)4 along the lowest singlet and triplet surfaces to form Fe(CO)3. Calculations of transient spectra in both singlet and triplet states as well as spin-orbit coupling constants along key structural pathways provide evidence for intersystem crossing to the triplet ground state of Fe(CO)4. Thus, our work presents the first spectroscopic detection of transient excited states during ultrafast photodissociation of gas-phase Fe(CO)5 and challenges the long-standing assumption that triplet states do not play a role in the ultrafast dynamics.

Infrared Ion Spectroscopic Characterization of the Gaseous [Co(15-crown-5)(H2O)]2+ Complex

We report fingerprint infrared multiple-photon dissociation spectra of the gaseous monohydrated coordination complex of cobalt(II) and the macrocycle 1,4,7,10,13-pentaoxacyclopentadecane (or 15-crown-5), [Co(15-crown-5)(H2O)]2+. The metal-ligand complexes are generated using electrospray ionization, and their IR action spectra are recorded in a quadrupole ion trap mass spectrometer using the free-electron laser FELIX. The electronic structure and chelation motif are derived from spectral comparison with computed vibrational spectra obtained at the density functional theory level. We focus here on the gas-phase structure, addressing the question of doublet versus quartet spin multiplicity and the chelation geometry. We conclude that the gas-phase complex adopts a quartet spin state, excluding contributions of doublet species, and that the chelation geometry is pseudo-octahedral with the six oxygen centers of 15-crown-5 and H2O coordinated to the metal ion. We also address the possible presence of higher-energy conformers based on the IR spectral evidence and calculated thermodynamics.

How complex-surface interactions modulate the spin transition of Fe(II) SCO complexes supported on metallic surfaces?

The deposition of a prototypical spin-crossover [Fe(phen)2(NCS)2] complex on Au(111), Cu(111) and Ag(111) surfaces has been investigated by means of periodic DFT+U calculations, with the aim of understanding how different metallic surfaces affect the spin state switching. Our results show that adsorption is metal- and spin-dependent, with different preferred adsorption sites for the different surfaces and spin states. For the three considered surfaces adsorption energies are larger in the LS state than in the HS one, which increases the transition enthalpy by 58.7 kJ mol-1 for Cu(111), 14.6 kJ mol-1 for Au(111) and 9.6 kJ mol-1 for Ag(111) with respect to the free molecule. There is a clear correlation between this effect and the extent of the complex-surface interaction, which can be established from adsorption energies, surface-complex distances and charge density difference plots as: Cu(111) > Au(111) > Ag(111). Therefore, a stronger interaction with the surface produces a larger energy difference between two spin states, making the spin transition less probable to occur. Finally, our calculations show that it would be possible to probe the spin-state of the deposited molecules from the STM images, in line with the recent experimental results.

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.

Branching mechanism of photoswitching in an Fe(II) polypyridyl complex explained by full singlet-triplet-quintet dynamics

It has long been known that irradiation with visible light converts Fe(II) polypyridines from their low-spin (singlet) to high-spin (quintet) state, yet mechanistic interpretation of the photorelaxation remains controversial. Herein, we simulate the full singlet-triplet-quintet dynamics of the [Fe(terpy)₂]²⁺ (terpy = 2,2':6',2"-terpyridine) complex in full dimension, in order to clarify the complex photodynamics. Importantly, we report a branching mechanism involving two sequential processes: a dominant ³MLCT→³MC(³T_2g)→³MC(³T_1g)→⁵MC, and a minor ³MLCT→³MC(³T_2g)→⁵MC component. (MLCT = metal-to-ligand charge transfer, MC = metal-centered). While the direct ³MLCT→⁵MC mechanism is considered as a relevant alternative, we show that it could only be operative, and thus lead to competing pathways, in the absence of ³MC states. The quintet state is populated on the sub-picosecond timescale involving non-exponential dynamics and coherent Fe-N breathing oscillations. The results are in agreement with the available time-resolved experimental data on Fe(II) polypyridines, and fully describe the photorelaxation dynamics.

Toward High-Temperature Light-Induced Spin-State Trapping in Spin-Crossover Materials: The Interplay of Collective and Molecular Effects

Spin-crossover (SCO) materials display many fascinating behaviors including collective phase transitions and spin-state switching controlled by external stimuli, e.g., light and electrical currents. As single-molecule switches, they have been fêted for numerous practical applications, but these remain largely unrealized-partly because of the difficulty of switching these materials at high temperatures. We introduce a semiempirical microscopic model of SCO materials combining crystal field theory with elastic intermolecular interactions. For realistic parameters, this model reproduces the key experimental results including thermally induced phase transitions, light-induced spin-state trapping (LIESST), and reverse-LIESST. Notably, we reproduce and explain the experimentally observed relationship between the critical temperature of the thermal transition, T_1/2, and the highest temperature for which the trapped state is stable, T_LIESST, and explain why increasing the stiffness of the coordination sphere increases T_LIESST. We propose strategies to design SCO materials with higher T_LIESST: optimizing the spin-orbit coupling via heavier atoms (particularly in the inner coordination sphere) and minimizing the enthalpy difference between the high-spin (HS) and low-spin (LS) states. However, the most dramatic increases arise from increasing the cooperativity of the spin-state transition by increasing the rigidity of the crystal. Increased crystal rigidity can also stabilize the HS state to low temperatures on thermal cycling yet leave the LS state stable at high temperatures following, for example, reverse-LIESST. We show that such highly cooperative systems offer a realistic route to robust room-temperature switching, demonstrate this in silico, and discuss material design rationale to realize this.

Trend in light-induced excited-state spin trapping in Fe(II)-based spin crossover systems

A computational study of the light-induced excited spin-state trapping (LIESST) in a number of Fe(II) spin crossover complexes, coordinated by monodentate, bidentate and multidentate ligands is carried out, with the goal to uncover the trend in the low temperature relaxation rate. A nine order of magnitude change in low temperature relaxation rate is observed among the complexes. The trend is rationalized in terms of the change in metal-ligand covalency, numerically estimated by the crystal orbital Hamiltonian population, thus influencing the back donation or delocalization of the electrons from the low-lying Fe(II)-centered molecular orbital to the empty low-lying ligand-centered π* antibonding molecular orbitals.

Toward Simulation of Fe(II) Low-Spin → High-Spin Photoswitching by Synergistic Spin-Vibronic Dynamics

A new theoretical approach is presented and applied for the simulation of Fe(II) low-spin (LS, singlet, t_2g⁶e_g⁰) → high-spin (HS, quintet, t_2g⁴e_g²) photoswitching dynamics of the octahedral model complex [Fe(NCH)₆]²⁺. The utilized synergistic methodology heavily exploits the strengths of complementary electronic structure and spin-vibronic dynamics methods. Specifically, we perform 3D quantum dynamics (QD) and full-dimensional trajectory surface hopping (TSH, in conjunction with a linear vibronic coupling model), with the modes for QD selected by TSH. We follow a hybrid approach which is based on the application of time-dependent density functional theory (TD-DFT) excited-state potential energy surfaces (PESs) and multiconfigurational second-order perturbation theory (CASPT2) spin–orbit couplings (SOCs). Our method delivers accurate singlet–triplet–quintet intersystem crossing (ISC) dynamics, as assessed by comparison to our recent high-level ab initio simulations and related time-resolved experimental data. Furthermore, we investigate the capability of our simulations to identify the location of ISCs. Finally, we assess the approximation of constant SOCs (calculated at the Franck–Condon geometry), whose validity has central importance for the combination of TD-DFT PESs and CASPT2 SOCs. This efficient methodology will have a key role in simulating LS → HS dynamics for more complicated cases, involving higher density of states and varying electronic character, as well as the analysis of ultrafast experiments.

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.

Dependence of the Fe(II)-Gallic Acid Coordination Compound Formation Constant on the pH

One important property of tannins involves their ability to form coordination compounds with metal ions, which is vital for the bioavailability of these ions, as well as for the antibacterial and antioxidative activities of tannins. In this study, the pH dependence of interactions between gallic acid, one of the basic building blocks of tannins, and Fe(II) ions, was investigated using UV/Vis spectroscopy, in conjunction with density functional theory (DFT) calculations. Moreover, two models were developed to explain the processes taking place in the solution. The first model treated the reaction as a simple bimolecular process while the second also considered the protolytic equilibrium, which was proven very successful in discerning the pH dependence of formation constants, and whose assumptions were well supported by DFT calculations. We showed that the two-time deprotonated gallic acid species forms the coordination compound with Fe(II) ions in a 1:1 molar ratio. To gain better insight into the process, the coordination compound formation was also studied using various DFT functionals, which further supported the model results. Furthermore, due to the relatively low sample amounts needed, the methodology developed here will be useful to study compounds that are more difficult to isolate.

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.

Theoretical inspection of the spin-crossover [Fe(tzpy)2(NCS)2] complex on Au(100) surface

We explore the deposition of the spin-crossover [Fe(tzpy)₂(NCS)₂] complex on the Au(100) surface by means of density functional theory (DFT) based calculations. Two different routes have been employed: low-cost finite cluster-based calculations, where both the Fe complex and the surface are maintained fixed while the molecule approaches the surface; and periodic DFT plane-wave calculations, where the surface is represented by a four-layer slab and both the molecule and surface are relaxed. Our results show that the bridge adsorption site is preferred over the on-top and fourfold hollow ones for both spin states, although they are energetically close. The LS molecule is stabilized by the surface, and the HS-LS energy difference is enhanced by about 15%-25% once deposited. The different Fe ligand field for LS and HS molecules manifests on the composition and energy of the low-lying bands. Our simulated STM images indicate that it is possible to distinguish the spin state of the deposited molecules by tuning the bias voltage of the STM tip. Finally, it should be noted that the use of a reduced size cluster to simulate the Au(100) surface proves to be a low-cost and reliable strategy, providing results in good agreement with those resulting from state-of-the-art periodic calculations for this system.

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.

Spin-state energetics of metallocenes: How do best wave function and density functional theory results compare with the experimental data?

We benchmark the accuracy of quantum-chemical methods, including wave function theory methods [coupled cluster theory at the CCSD(T) level, multiconfigurational perturbation-theory (CASPT2, NEVPT2) and internally contracted multireference configuration interaction (MRCI)] and 30 density functional theory (DFT) approximations, in reproducing the spin-state splittings of metallocenes. The reference values of the electronic energy differences are derived from the experimental spin-crossover enthalpy for manganocene and the spectral data of singlet-triplet transitions for ruthenocene, ferrocene, and cobaltocenium. For ferrocene and cobaltocenium we revise the previous experimental interpretations regarding the lowest triplet energy; our argument is based on the comparison with the lowest singlet excitation energy and herein reported, carefully determined absorption spectrum of ferrocene. When deriving vertical energies from the experimental band maxima, we go beyond the routine vertical energy approximation by introducing vibronic corrections based on simulated vibrational envelopes. The benchmarking result confirms the high accuracy of the CCSD(T) method (in particular, for UCCSD(T) based on Hartree-Fock orbitals we find for our dataset: maximum error 0.12 eV, weighted mean absolute error 0.07 eV, weighted mean signed error 0.01 eV). The high accuracy of the single-reference method is corroborated by the analysis of a multiconfigurational character of the complete active space wave function for the triplet state of ferrocene. On the DFT side, our results confirm the non-universality problem with approximate functionals. The present study is an important step toward establishing an extensive and representative benchmark set of experiment-derived spin-state energetics for transition metal complexes.

A model electronic Hamiltonian to describe low-lying d-d and metal-to-ligand charge-transfer excited states of [Fe(bpy)3 ]2

A simple practical method to compute both d-d and metal-to-ligand charge-transfer (MLCT) excited states of iron(II) polypyridyl complexes is proposed for use in simulation studies. Specifically, a model electronic Hamiltonian developed previously for d-d excited states of [Fe(bpy)₃ ]²⁺ is extended to deal with low-lying MLCT excited states simultaneously by including the MLCT electronic configurations into the basis functions of the model Hamiltonian. As a first attempt, parameters in the model Hamiltonian matrix elements are determined by using density functional theory (DFT) and time-dependent (TD-)DFT calculation results as benchmarks. To examine the performance of the model Hamiltonian, the potential energy curves along the interpolation between the lowest singlet and quintet state structures are compared to those from the (TD-)DFT calculations and to those from CASPT2 calculations in literature. The electronic absorption spectrum computed through molecular dynamics simulation is compared to the experimental spectrum. The spin-orbit couplings at the ground state structure are also compared to those from wavefunction-based ab initio electronic structure calculations. The results indicate that the constructed model Hamiltonian provides reasonable information on both the low-lying d-d and MLCT excited states of [Fe(bpy)₃ ]²⁺ .

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.

Benchmarking quantum chemistry methods for spin-state energetics of iron complexes against quantitative experimental data

The accuracy of relative spin-state energetics predicted by selected quantum chemistry methods: coupled cluster theory at the CCSD(T) level, multiconfigurational perturbation theory (CASPT2, NEVPT2), multireference configuration interaction at the MRCISD+Q level, and a number of DFT methods, is quantitatively evaluated by comparison with the experimental data of four octahedral iron complexes. The available experimental data, either spin-forbidden transition energies or spin crossover enthalpies, are corrected for relevant environmental effects in order to derive the quantitative benchmark set of iron spin-state energetics. Comparison of theory predictions with the resulting reference data: (1) validates the high accuracy of the CCSD(T) method, particularly when based on Kohn-Sham orbitals, giving the maximum error below 2 kcal mol-1 and the mean absolute error (MAE) below 1 kcal mol-1; (2) corroborates the tendency of CASPT2 to systematically overstabilize higher-spin states by up to 5.5 kcal mol-1; (3) confirms that the latter problem is partly remedied by the recently proposed CASPT2/CC approach [Phung et al., J. Chem. Theory Comput., 2018, 14, 2446-2455]; (4) demonstrates that NEVPT2 performs worse than CASPT2, by giving errors up to 7 kcal mol-1; (5) shows that the accuracy of MRCISD+Q spin-state energetics strongly depends on the size-consistency correction: the Davidson-Silver and Pople corrections perform best (MAE < 3 kcal mol-1), whereas the standard Davidson correction is not recommended (MAE of 7 kcal mol-1). Only a few DFT methods (including the best performing ones identified in this study: B2PLYP-D3 and OPBE) are able to provide a balanced description of the spin-state energetics for all four studied iron complexes simultaneously, corroborating the non-universality problem of approximate density functionals.

Visualizing the coordination-spheres of photoexcited transition metal complexes with ultrafast hard X-rays

The concept of coordination sphere (CS) is central to the rational development of hierarchical molecular assemblies in modern chemistry.

Excited-state solvation structure of transition metal complexes from molecular dynamics simulations and assessment of partial atomic charge methods

Excited-state solvation structure (radial distribution function) of transition metal complexes by classical and mixed quantum-classical (QM/MM) molecular dynamics simulations.

Observation of the First Spin Crossover in an Iron(II) Complex with an S6 Coordination Environment: Tris[bis(N,N-diethylamino)carbeniumdithiocarboxylato]iron(II) Hexafluorophosphate

For the first time, the spin-crossover (SCO) phenomenon has been observed in an Fe^II -S₆ system in a tris(chelate)-type iron(II) complex with a zwitterionic sulfur donor bidentate, bis(N,N-diethylamino)carbeniumdithiocarboxylate (EtL), [Fe^II (EtL)₃ ](PF₆ )₂ (1), as synthesized by the reaction of a precursor complex [Fe^II (CH₃ CN)₆ ](PF₆ )₂ with EtL. In the solid state, the high-spin (HS) d⁶ state at ambient temperature and the low-spin (LS) d⁶ state at temperatures lower than approximately 240 K were evidenced by magnetic measurements with SQUID and Mössbauer spectra in the temperature range 4-290 K. X-ray analyses of the crystals at various temperatures disclosed that the distorted trigonal prismatic coordination environments essentially do not change; however, contraction of Fe-S distances by approximately 10 % (0.22 Å), ordering of alkyl groups in EtL and PF₆ ^- counteranions, and formation of significant intermolecular S⋅⋅⋅S interactions between adjacent molecules (average distances of 3.59 Å) take place during the transition from the HS to the LS state. A large decrease in the volume of the formula unit (78.1 ų ) might be responsible for the large activation barrier, thereby resulting in a slow phase transition upon cooling.

Monitoring Ultrafast Spin Crossover Intermediates in an Iron(II) Complex by Broad Band Stimulated X-ray Raman Spectroscopy

The photoinduced spin crossover dynamics of transition metal complexes is of fundamental scientific importance and is used for sensor device applications and solar energy harvesting. Current X-ray and optical spectroscopy experiments for [Fe^II(bpy)₃], an archetypal earth-abundant metal complex, show conflicting spin dynamics. We have simulated the broad band transient X-ray absorption and hybrid (broad + narrow band) X-ray stimulated Raman signals at the N and Fe K-edges of the key excited state intermediates involved in the spin crossover process of this complex. We find that these signals are much more sensitive to electron and spin populations than transition absorption and may be useful in the design of photovoltaic and artificial photosynthetic systems.

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.

Computational Spectroscopic Characterization of a Bistable Binuclear Complex [(CO)2(benzoate)FeII/III(terephthalate)CoIII/II(benzoate)(CO)2]+

Electric dipole moments, polarizabilities, and IR, Raman, optical rotatory dispersion, and electronic and vibrational circular dichroism spectra of the four cis–trans isomers of the proposed [(CO)2(benzoate)FeII/III(terephthalate)CoIII/II(benzoate)(CO)2]+ binuclear complex, having bistablity due to intramolecular charge transfer (IMCT), is investigated using the time-dependent density functional theory ((TD)DFT) B3LYP/6–31G(d,p)[LanL2DZ] method. Results show that the two FeII–CoIII and FeIII–CoII IMCT states of this binuclear complex have distinctly different spectroscopic, optical, and electric response properties, and are sensitive to the cis–trans arrangement of the ligands around the two metallic centres. Furthermore, intrinsic reaction coordinates inter-connecting the two IMCT states are identified using the Duschinsky matrix method. Only one or two of the normal coordinates remain almost (above 80 %) intact during the IMCT reaction which denotes global changes in the bonding strengths and potential energy hypersurface of this bistable binuclear complex. Analysis of the calculated spin densities characterizes the IMCT transition state structures of the trans–trans, cis–cis, and trans–cis isomers as early, early, and late transition states, respectively.

Computational characterization of competing energy and electron transfer states in bimetallic donor-acceptor systems for photocatalytic conversion

The rapidly growing interest in photocatalytic systems for direct solar fuel production such as hydrogen generation from water splitting is grounded in the unique opportunity to achieve charge separation in molecular systems provided by electron transfer processes. In general, both photoinduced and catalytic processes involve complicated dynamics that depend on both structural and electronic effects. Here the excited state landscape of metal centered light harvester-catalyst pairs is explored using density functional theory calculations. In weakly bound systems, the interplay between structural and electronic factors involved can be constructed from the various mononuclear relaxed excited states. For this study, supramolecular states of electron transfer and excitation energy transfer character have been constructed from constituent full optimizations of multiple charge/spin states for a set of three Ru-based light harvesters and nine transition metal catalysts (based on Ru, Rh, Re, Pd, and Co) in terms of energy, structure, and electronic properties. The complete set of combined charge-spin states for each donor-acceptor system provides information about the competition of excited state energy transfer states with the catalytically active electron transfer states, enabling the identification of the most promising candidates for photocatalytic applications from this perspective.

Diffusion Monte Carlo Perspective on the Spin-State Energetics of [Fe(NCH)6](2.)

The energy difference between the high spin and the low spin state of the model compound [Fe(NCH)6](2+) is investigated by means of Diffusion Monte Carlo (DMC), where special attention is dedicated to analyzing the effect of the fix node approximation on the accuracy of the results. For this purpose, we compare several Slater-Jastrow and multireference Slater-Jastrow trial wave functions. We found that a Slater-Jastrow trial wave function constructed with the generalized Kohn-Sham orbitals from hybrid DFT represents the optimal choice. This is understood by observing that hybrid functionals account for the subtle balance between exchange and correlation effects and the respective orbitals accurately describe the ligand-metal hybridization as well as the charge reorganization accompanying the spin transition. Finally the DMC results are compared with those obtained by Hartree-Fock, DFT, CASSCF, and CASPT2. While there is no clear reference value for the high spin-low spin energy difference, DMC and high level CCSD(T) calculations agree within around 0.3 eV.

Manipulating charge transfer excited state relaxation and spin crossover in iron coordination complexes with ligand substitution

Optical and X-ray free-electron laser measurements reveal ligand substitution in an Fe(ii)-centered complex extends its MLCT lifetime.

Developing light-harvesting and photocatalytic molecules made with iron could provide a cost effective, scalable, and environmentally benign path for solar energy conversion. To date these developments have been limited by the sub-picosecond metal-to-ligand charge transfer (MLCT) electronic excited state lifetime of iron based complexes due to spin crossover – the extremely fast intersystem crossing and internal conversion to high spin metal-centered excited states. We revitalize a 30 year old synthetic strategy for extending the MLCT excited state lifetimes of iron complexes by making mixed ligand iron complexes with four cyanide (CN^–) ligands and one 2,2′-bipyridine (bpy) ligand. This enables MLCT excited state and metal-centered excited state energies to be manipulated with partial independence and provides a path to suppressing spin crossover. We have combined X-ray Free-Electron Laser (XFEL) Kβ hard X-ray fluorescence spectroscopy with femtosecond time-resolved UV-visible absorption spectroscopy to characterize the electronic excited state dynamics initiated by MLCT excitation of [Fe(CN)₄(bpy)]^2–. The two experimental techniques are highly complementary; the time-resolved UV-visible measurement probes allowed electronic transitions between valence states making it sensitive to ligand-centered electronic states such as MLCT states, whereas the Kβ fluorescence spectroscopy provides a sensitive measure of changes in the Fe spin state characteristic of metal-centered excited states. We conclude that the MLCT excited state of [Fe(CN)₄(bpy)]^2– decays with roughly a 20 ps lifetime without undergoing spin crossover, exceeding the MLCT excited state lifetime of [Fe(2,2′-bipyridine)₃]²⁺ by more than two orders of magnitude.

Fe N-Heterocyclic Carbene Complexes as Promising Photosensitizers

The photophysics and photochemistry of transition metal complexes (TMCs) has long been a hot field of interdisciplinary research. Rich metal-based redox processes, together with a high variety in electronic configurations and excited-state dynamics, have rendered TMCs excellent candidates for interconversion between light, chemical, and electrical energies in intramolecular, supramolecular, and interfacial arrangements. In specific applications such as photocatalytic organic synthesis, photoelectrochemical cells, and light-driven supramolecular motors, light absorption by a TMC-based photosensitizer and subsequent excited-state energy or electron transfer constitute essential steps. In this context, TMCs based on rare and expensive metals, such as ruthenium and iridium, are frequently employed as photosensitizers, which is obviously not ideal for large-scale implementation. In the search for abundant and environmentally benign solutions, six-coordinate Fe(II) complexes (Fe(II)L6) have been widely considered as highly desirable alternatives. However, not much success has been achieved due to the extremely short-lived triplet metal-to-ligand charge transfer ((3)MLCT) excited state that is deactivated by low-lying metal-centered (MC) states on a 100 fs time scale. A fundamental strategy to design useful Fe-based photosensitizers is thus to destabilize the MC states relative to the (3)MLCT state by increasing the ligand field strength, with special focus on making eg σ* orbitals on the Fe center energetically less accessible. Previous efforts to directly transplant successful strategies from Ru(II)L6 complexes unfortunately met with limited success in this regard, despite their close chemical kinship. In this Account, we summarize recent promising results from our and other groups in utilizing strongly σ-donating N-heterocyclic carbene (NHC) ligands to make strong-field Fe(II)L6 complexes with significantly extended (3)MLCT lifetimes. Already some of the first homoleptic bis(tridentate) complexes incorporating (CNHC^Npyridine^CNHC)-type ligands gratifyingly resulted in extension of the (3)MLCT lifetime by more than 2 orders of magnitude compared to the parental [Fe(tpy)2](2+) (tpy = 2,2':6',2″-terpyridine) complex. Quantum chemical (QC) studies also revealed that the (3)MC instead of the (5)MC state likely dictates the deactivation of the (3)MLCT state, a behavior distinct from traditional Fe(II)L6 complexes but rather resembling Ru analogues. A heteroleptic Fe(II) NHC complex featuring mesoionic bis(1,2,3-triazol-5-ylidene) (btz) ligands also delivered a 100-fold elongation of the (3)MLCT lifetime relative to its parental [Fe(bpy)3](2+) (bpy = 2,2'-bipyridine) complex. Again, a Ru-like deactivation mechanism of the (3)MLCT state was indicated by QC studies. With a COOH-functionalized homoleptic complex, a record (3)MLCT lifetime of 37 ps was recently observed on an Al2O3 nanofilm. As a proof of concept, it was further demonstrated that the significant improvement in the (3)MLCT lifetime indeed benefits efficient light harvesting with Fe(II) NHC complexes. For the first time, close-to-unity electron injection from the lowest-energy (3)MLCT state to a TiO2 nanofilm was achieved by a stable Fe(II) complex. This is in complete contrast to conventional Fe(II)L6-derived photosensitizers that could only make use of high-energy photons. These exciting results significantly broaden the understanding of the fundamental photophysics and photochemistry of d(6) Fe(II) complexes. They also open up new possibilities to develop solar energy-converting materials based on this abundant, inexpensive, and intrinsically nontoxic element.

Molecular and Interfacial Calculations of Iron(II) Light Harvesters

Iron-carbene complexes show considerable promise as earth-abundant light-harvesters, and adsorption onto nanostructured TiO2 is a crucial step for developing solar energy applications. Intrinsic electron injection capabilities of such promising Fe(II) N-heterocyclic complexes (Fe-NHC) to TiO2 are calculated here, and found to correlate well with recent experimental findings of highly efficient interfacial injection. First, we examine the special bonding characteristics of Fe-NHC light harvesters. The excited-state surfaces are examined using density functional theory (DFT) and time-dependent DFT (TD-DFT) to explore relaxed excited-state properties. Finally, by relaxing an Fe-NHC adsorbed on a TiO2 nanocluster, we show favorable injection properties in terms of interfacial energy level alignment and electronic coupling suitable for efficient electron injection of excited electrons from the Fe complex into the TiO2 conduction band on ∼100 fs time scales.

Characterizing the Solvated Structure of Photoexcited [Os(terpy)₂](2+) with X-ray Transient Absorption Spectroscopy and DFT Calculations

Characterizing the geometric and electronic structures of individual photoexcited dye molecules in solution is an important step towards understanding the interfacial properties of photo-active electrodes. The broad family of “red sensitizers” based on osmium(II) polypyridyl compounds often undergoes small photo-induced structural changes which are challenging to characterize. In this work, X-ray transient absorption spectroscopy with picosecond temporal resolution is employed to determine the geometric and electronic structures of the photoexcited triplet state of [Os(terpy)₂]²⁺ (terpy: 2,2′:6′,2″-terpyridine) solvated in methanol. From the EXAFS analysis, the structural changes can be characterized by a slight overall expansion of the first coordination shell [OsN₆]. DFT calculations supports the XTA results. They also provide additional information about the nature of the molecular orbitals that contribute to the optical spectrum (with TD-DFT) and the near-edge region of the X-ray spectra.

Ammonia-modified Co(II) sites in zeolites: spin and electron density redistribution through the Co(II)-NO bond

Electronic factors essential for the bonding of a non-innocent NO ligand to ammonia-modified Co(2+) sites in cobalt-exchanged zeolites are examined for small cluster models using DFT and advanced correlated wave function calculations. The analysis of charge transfer processes between the NO ligand and the cobalt center involves two protocols: valence-bond expansion of the multiconfiguration CASSCF wave function (in terms of fragment-localized active orbitals) and spin-resolved natural orbitals for chemical valence (SR-NOCV). Applicability of SR-NOCV analysis to transition metal complexes involving non-innocent fragments is critically assessed and the approach based on the CASSCF wave function turns out to be much more robust and systematic for all studied models. It is shown that the character and direction of electron density redistribution through the Co-N-O bond, quantified by relative share of the Co(II)-NO(0), Co(III)-NO(-), and Co(I)-NO(+) resonance structures in the total wave function, fully rationalize the activation of the N-O bond upon NH3 co-ligation (evidenced by calculated and measured red-shift of the NO stretching frequency and commonly ascribed to enhanced backdonation). The huge red-shift of νN-O is attributed to an effective electron transfer between the ammonia-modified Co(ii) centers and the NO antibonding π*-orbitals (related to the increased share of the Co(III)-NO(-) form). Unexpectedly, the effect is stronger for the singlet complex with three NH3 ligands than for that with five NH3 ligands bound to the cobalt center. Our results also indicate that high-efficiency electron transfers between the Co(ii) center and the NO ligand may be enabled for the selected spin state and disabled for the other spin state of the adduct. This illustrates how the cobalt center may serve to fine-tune the electronic communication between the NO ligand and its binding site.

Understanding and tuning the properties of redox-accumulating manganese helicates

The multiple redox transitions of pentanuclear Mn clusters and the tuning of their redox potentials by ligand design are investigated computationally.

TD-DFT study of the light-induced spin crossover of Fe(iii) complexes

Two light-induced spin-crossover Fe(iii) compounds have been studied with time-dependent density functional theory (TD-DFT) to investigate the deactivation mechanism and the role of the ligand-field states as intermediates in this process.

Insight into the light-induced spin crossover of [Fe(bpy)3]2+ in aqueous solution from molecular dynamics simulation of d–d excited states

Lifetimes of triplet d–d states were evaluated through molecular dynamics simulations to gain insight into relaxation dynamics of aqueous [Fe(bpy)₃]²⁺.

Sub-50-fs photoinduced spin crossover in [Fe(bpy)₃]²⁺

It is known that excitation by visible light of the singlet metal-to-ligand charge-transfer ((1)MLCT) states of Fe(II) complexes leads to population of the lowest-lying high-spin quintet state ((5)T) with unity quantum yield. Here we investigate this so-called spin crossover (SCO) transition in aqueous iron(II)tris(bipyridine). We use pump-probe transient absorption spectroscopy with a high time resolution of <60 fs in the ultraviolet probe range, in which the (5)T state absorbs, and of <40 fs in the visible probe range, in which both the hot MLCT state and the (5)T state absorb. Our results show that the (5)T state is impulsively populated in less than 50 fs, which is the time we measured for the depopulation of the MCLT manifold. We propose that non-totally-symmetric modes mediate the process, possibly high-frequency modes of the bipyridine (bpy) ligand. These results show that even though the SCO process in Fe(II) complexes represents a strongly spin-forbidden (ΔS = 2) two-electron transition, spin flipping occurs at near subvibrational times and is intertwined with the electron and structural dynamics of the system.

Fe(II)-Polypyridines as Chromophores in Dye-Sensitized Solar Cells: A Computational Perspective

Over the past two decades, dye-sensitized solar cells (DSSCs) have become a viable and relatively cheap alternative to conventional crystalline silicon-based systems. At the heart of a DSSC is a wide band gap semiconductor, typically a TiO2 nanoparticle network, sensitized with a visible light absorbing chromophore. Ru(II)-polypyridines are often utilized as chromophores thanks to their chemical stability, long-lived metal-to-ligand charge transfer (MLCT) excited states, tunable redox potentials, and near perfect quantum efficiency of interfacial electron transfer (IET) into TiO2. More recently, coordination compounds based on first row transition metals, such as Fe(II)-polypyridines, gained some attention as potential sensitizers in DSSCs due to their low cost and abundance. While such complexes can in principle sensitize TiO2, they do so very inefficiently since their photoactive MLCT states undergo intersystem crossing (ISC) into low-lying metal-centered states on a subpicosecond time scale. Competition between the ultrafast ISC events and IET upon initial excitation of Fe(II)-polypyridines is the main obstacle to their utilization in DSSCs. Suitability of Fe(II)-polypyridines to serve as sensitizers could therefore be improved by adjusting relative rates of the ISC and IET processes, with the goal of making the IET more competitive with ISC. Our research program in computational inorganic chemistry utilizes a variety of tools based on density functional theory (DFT), time-dependent density functional theory (TD-DFT) and quantum dynamics to investigate structure-property relationships in Fe(II)-polypyridines, specifically focusing on their function as chromophores. One of the difficult problems is the accurate determination of energy differences between electronic states with various spin multiplicities (i.e., (1)A, (1,3)MLCT, (3)T, (5)T) in the ISC cascade. We have shown that DFT is capable of predicting the trends in the energy ordering of these electronic states in a set of structurally related complexes with the help of appropriate benchmarks, based either on experimental data or higher-level ab initio calculations. Models based on TD-DFT and quantum dynamics approaches have proven very useful in understanding IET processes in Fe(II)-polypyridine-TiO2 assemblies. For example, they helped us to elucidate the origin of "band selective" sensitization in the [Fe(bpy-dca)2(CN)2]-TiO2 assembly (bpy-dca = 2,2'-bipyridine-4,4'-dicarboxylic acid), first observed by Ferrere and Gregg [ Ferrere , S. ; Gregg , B. A. J. Am. Chem. Soc. 1998 , 120 , 843 . ]. They also shed light on the relationship between the linker group that anchors Fe(II)-polypyridines onto the TiO2 surface and the speed of IET in Fe(II)-polypyridine-TiO2 assemblies. More interestingly, our results show that the IET efficiency is strongly correlated with the amount of electron density on the linker group and that one can obtain insights into the IET in dye-semiconductor assemblies based on ground state electronic structure calculations alone. This may be useful for quick screening of a large number of complexes for use as potential sensitizers in DSSCs, especially if followed up by TD-DFT and quantum dynamics simulations for selected target compounds to confirm efficient sensitization. While our focus over the past few years has been exclusively on Fe(II)-polypyridines, the computational strategies outlined in this Account are applicable to a wide variety of sensitizers.

Detailed Characterization of a Nanosecond-Lived Excited State: X-ray and Theoretical Investigation of the Quintet State in Photoexcited [Fe(terpy)2]2

Theoretical predictions show that depending on the populations of the Fe 3d _xy , 3d _xz , and 3d _yz orbitals two possible quintet states can exist for the high-spin state of the photoswitchable model system [Fe(terpy)₂]²⁺. The differences in the structure and molecular properties of these ⁵B₂ and ⁵E quintets are very small and pose a substantial challenge for experiments to resolve them. Yet for a better understanding of the physics of this system, which can lead to the design of novel molecules with enhanced photoswitching performance, it is vital to determine which high-spin state is reached in the transitions that follow the light excitation. The quintet state can be prepared with a short laser pulse and can be studied with cutting-edge time-resolved X-ray techniques. Here we report on the application of an extended set of X-ray spectroscopy and scattering techniques applied to investigate the quintet state of [Fe(terpy)₂]²⁺ 80 ps after light excitation. High-quality X-ray absorption, nonresonant emission, and resonant emission spectra as well as X-ray diffuse scattering data clearly reflect the formation of the high-spin state of the [Fe(terpy)₂]²⁺ molecule; moreover, extended X-ray absorption fine structure spectroscopy resolves the Fe-ligand bond-length variations with unprecedented bond-length accuracy in time-resolved experiments. With ab initio calculations we determine why, in contrast to most related systems, one configurational mode is insufficient for the description of the low-spin (LS)-high-spin (HS) transition. We identify the electronic structure origin of the differences between the two possible quintet modes, and finally, we unambiguously identify the formed quintet state as ⁵E, in agreement with our theoretical expectations.