A Computational Study of Carbene Ligand Stabilization of Biomimetic Models of the Rotated H red State of [FeFe]‐Hydrogenase

Probing the performance of DFT in the structural characterization of [FeFe] hydrogenase models

In this work, a series of DFT and DFT-D methods is combined with double-ζ basis sets to benchmark their performance in predicting the structures of five newly synthesized hexacarbonyl diiron complexes with a bridging ligand featuring a μ-S2C3 motif in a ring-containing unit functionalized with aromatic groups. Such complexes have been considered as [FeFe] hydrogenase catalytic site models with potential for eco-friendly energetic applications. According to this assessment, r2SCAN is identified as the density functional recommended for the reliable description of the molecular and crystal structures of the herein studied models. However, the butterfly (μ-S)2Fe2 core of the models demonstrates a minor deformation of its optimized geometry obtained from both molecular and periodic calculations. The FeFe bond length is slightly underestimated while the FeS bonds tend to be too long. Adding the D3(BJ) correction to r2SCAN does not lead to any improvement in the calculated structures.

Can carbene decorated [FeFe]-hydrogenase model complexes catalytically produce dihydrogen? An insight from theory

Cyclic alkyl amino carbene (CAAC) anchored [FeFe]-hydrogenase model complex featuring rotated conformation at one of the iron centers are found to be promising candidate for effective production of dihydrogen. A stepwise comparison of the complete mechanism using the CAAC stabilized model complex [1]⁰ has been performed with that of an experimentally isolated one ([2]⁰). Interestingly, the reduction events involved in the catalytic cycles are found to be more favorable than those previously reported for a similar experimentally known system. Furthermore, the computed ΔpK_a values indicate that the distal iron center with a vacant coordination site is more basic compared to the amino nitrogen atom of the azadithiolate bridge. We also made an attempt to determine the oxidation states of the iron centers for the intermediates involved in the catalytic cycles on the basis of the computed Mössbauer isomer shift and Mulliken spin density values.