Applying Nuclear Forward Scattering as In Situ and Operando Tool for the Characterization of FeN4 Moieties in the Hydrogen Evolution Reaction

Nuclear forward scattering (NFS) is a synchrotron-based technique relying on the recoil-free nuclear resonance effect similar to Mössbauer spectroscopy. In this work, we introduce NFS for in situ and operando measurements during electrocatalytic reactions. The technique enables faster data acquisition and better discrimination of certain iron sites in comparison to Mössbauer spectroscopy. It is directly accessible at various synchrotrons to a broad community of researchers and is applicable to multiple metal isotopes. We demonstrate the power of this technique with the hydrogen evolution mechanism of an immobilized iron porphyrin supported on carbon. Such catalysts are often considered as model systems for iron-nitrogen-carbon (FeNC) catalysts. Using in situ and operando NFS in combination with theoretical predictions of spectroscopic data enables the identification of the intermediate that is formed prior to the rate-determining step. The conclusions on the reaction mechanism can be used for future optimization of immobilized molecular catalysts and metal-nitrogen-carbon (MNC) catalysts.

Ultrafast photogeneration of a metal-organic nitrene from 1,1'-diazidoferrocene

Ferrocene and its derivatives have fascinated chemists for more than 70 years, not least due to the analogies with the properties of benzene. Despite these similarities, the obvious difference between benzene and ferrocene is the presence of an iron ion and hence the availability of d-orbitals for properties and reactivity. Phenylnitrene with its rich photochemistry can be considered an analogue of nitrenoferrocene. As with most organic and inorganic nitrenes, nitrenoferrocene can be obtained by irradiating the azide precursor. We study the photophysical and photochemical processes of dinitrogen release from 1,1'-diazidoferrocene to form 1-azido-1'-nitrenoferrocene with UV-pump-mid-IR-probe transient absorption spectroscopy and time-dependent density functional theory calculations including spin-orbit coupling. An intermediate with a bent azide moiety is identified that is pre-organised for dinitrogen release via a low-lying transition state. The photochemical decay paths on the singlet and triplet surfaces including the importance of spin-orbit coupling are discussed. We compare our findings with the processes discussed for photochemical dinitrogen activation and highlight implications for the photochemistry of azides more generally.

FeN4 Environments upon Reduction: A Computational Analysis of Spin States, Spectroscopic Properties, and Active Species

FeN4 motifs, found, for instance, in bioinorganic chemistry as heme-type cofactors, play a crucial role in man-made FeNC catalysts for the oxygen reduction reaction. Such single-atom catalysts are a potential alternative to platinum-based catalysts in fuel cells. Since FeNC catalysts are prepared via pyrolysis, the resulting materials are amorphous and contain side phases and impurities. Therefore, the geometric and electronic nature of the catalytically active FeN4 site remains to be clarified. To further understand the behavior of FeN4 centers in electrochemistry and their expected spectroscopic behavior upon reduction, we investigate two FeN4 environments (pyrrolic and pyridinic). These are represented by the model complexes [Fe(TPP)Cl] and [Fe(phen2N2)Cl], where TPP = tetraphenylporphyrin and phen = 1,10-phenanthroline. We predict their Mössbauer, UV-vis, and NRV spectral data using density functional theory as windows into their electronic structure differences. By varying the axial ligand, we further show how well small chemical changes in both complexes can be discerned. We find that the differences in ligand field strength in pyrrolic and pyridinic coordination result in different spin ground states, which in turn leads to distinct Mössbauer spectroscopic properties. As a result, pyrrolic nitrogen donors with a weaker ligand field are predicted to show more pronounced spectroscopic differences under in situ and operando conditions, while pyridinic nitrogen donors are expected to show less pronounced spectroscopic changes upon reduction and/or ligand loss. We therefore suggest that a weaker ligand field leads to better detectability of catalytic intermediates in in situ and operando experiments.

Water co-catalysis in aerobic olefin epoxidation mediated by ruthenium oxo complexes

We report the development of a versatile Ru-porphyrin catalyst system which performs the aerobic epoxidation of aromatic and aliphatic (internal) alkenes under mild conditions, with product yields of up to 95% and turnover numbers (TON) up to 300. Water is shown to play a crucial role in the reaction, significantly increasing catalyst efficiency and substrate scope. Detailed mechanistic investigations employing both computational studies and a range of experimental techniques revealed that water activates the RuVI di-oxo complex for alkene epoxidation via hydrogen bonding, stabilises the RuIV mono-oxo intermediate, and is involved in the regeneration of the RuVI di-oxo complex leading to oxygen atom exchange. Distinct kinetics are obtained in the presence of water, and side reactions involved in catalyst deactivation have been identified.

AOMadillo: A program for fitting angular overlap model parameters

The angular overlap model (AOM) is an established parameterization scheme within ligand field theory (LFT). In principle, its application is fairly straightforward, but can be tedious and involve a trial-and-error approach to identify and judge the best set of parameters. With the availability of quantum chemical methods to predict d-d transitions in transition metal complexes, a rich source of computational spectroscopic data with unambiguous assignments to electronic states is available. Herein, we present AOMadillo, a software package that is designed to interface the output of ab initio LFT calculations from the ORCA suite of programs and performs a least-squares fit for a chosen AOM parameterization. Many steps of the AOM parameterization are automated, so that scans of geometric parameters and evaluations of sets of similar complexes are convenient. The fitting routine is highly configurable, allowing the efficient evaluation of different parameter sets.

Vibrational Coherences of the Photoinduced Mixed-Valent Creutz-Taube Ion Revealed by Excited State Dynamics

A recent study of photoinduced mixed-valency in the one-electron reduced form (μ-pz)[RuII(NH3)5]24+ of the Creutz-Taube ion used transient absorption spectroscopy with vis-NIR broadband detection to uncover a mixed-valent excited state with a typical intervalence charge transfer band and a nanosecond lifetime [Pieslinger et al. Angew. Chem., Int. Ed. 2022, 61, e202211747]. Herein, we use excited state dynamics simulations with implicit solvation to elucidate the electronic and vibrational evolution in the first 10 ps after the optical excitation. A manifold of excited states with weak interaction between the metal centers is populated already at time zero due to the breakdown of the Condon approximation and dominates the population of electronic states at short time scales (<0.5 ps). A long-lived vibrational wave packet mostly confined to oscillations of the metal center-bridge distances is observed. The oscillations are traced to the electronic structure properties of states with weak metal-metal coupling. The long-lived mixed-valent excited state of the Creutz-Taube ion analogue is formed vibrationally cold and has a more compact geometry. While experimentally, intersystem crossing and vibrational relaxation were deduced to be completed within 1 ps, our analysis indicates that both processes might persist at longer times.

Nature of S-States in the Oxygen-Evolving Complex Resolved by High-Energy Resolution Fluorescence Detected X-ray Absorption Spectroscopy

Photosystem II, the water splitting enzyme of photosynthesis, utilizes the energy of sunlight to drive the four-electron oxidation of water to dioxygen at the oxygen-evolving complex (OEC). The OEC harbors a Mn4CaO5 cluster that cycles through five oxidation states Si (i = 0-4). The S3 state is the last metastable state before the O2 evolution. Its electronic structure and nature of the S2 → S3 transition are key topics of persisting controversy. Most spectroscopic studies suggest that the S3 state consists of four Mn(IV) ions, compared to the Mn(III)Mn(IV)3 of the S2 state. However, recent crystallographic data have received conflicting interpretations, suggesting either metal- or ligand-based oxidation, the latter leading to an oxyl radical or a peroxo moiety in the S3 state. Herein, we utilize high-energy resolution fluorescence detected (HERFD) X-ray absorption spectroscopy to obtain a highly resolved description of the Mn K pre-edge region for all S-states, paying special attention to use chemically unperturbed S3 state samples. In combination with quantum chemical calculations, we achieve assignment of specific spectroscopic features to geometric and electronic structures for all S-states. These data are used to confidently discriminate between the various suggestions concerning the electronic structure and the nature of oxidation events in all observable catalytic intermediates of the OEC. Our results do not support the presence of either peroxo or oxyl in the active configuration of the S3 state. This establishes Mn-centered storage of oxidative equivalents in all observable catalytic transitions and constrains the onset of the O-O bond formation until after the final light-driven oxidation event.

The Complex Reactivity of [(salen)Fe]2(μ-O) with HBpin and Its Implications in Catalysis

We report a detailed study into the method of precatalyst activation during alkyne cyclotrimerization. During these studies we have prepared a homologous series of Fe(III)-μ-oxo(salen) complexes and use a range of techniques including UV-vis, reaction monitoring studies, single crystal X-ray diffraction, NMR spectroscopy, and LIFDI mass spectrometry to provide experimental evidence for the nature of the on-cycle iron catalyst. These data infer the likelihood of ligand reduction, generating an iron(salan)-boryl complex as a key on-cycle intermediate. We use DFT studies to interrogate spin states, connecting this to experimentally identified diamagnetic and paramagnetic species. The extreme conformational flexibility of the salan system appears connected to challenges associated with crystallization of likely on-cycle species.

The Marcus dimension: identifying the nuclear coordinate for electron transfer from ab initio calculations

The Marcus model forms the foundation for all modern discussion of electron transfer (ET). In this model, ET results in a change in diabatic potential energy surfaces, separated along an ET nuclear coordinate. This coordinate accounts for all nuclear motion that promotes electron transfer. It is usually assumed to be dominated by a collective asymmetric vibrational motion of the redox sites involved in the ET. However, this coordinate is rarely quantitatively specified. Instead, it remains a nebulous concept, rather than a tool for gaining true insight into the ET pathway. Herein, we describe an ab initio approach for quantifying the ET coordinate and demonstrate it for a series of dinitroradical anions. Using sampling methods at finite temperature combined with density functional theory calculations, we find that the electron transfer can be followed using the energy separation between potential energy surfaces and the extent of electron localization. The precise nuclear motion that leads to electron transfer is then obtained as a linear combination of normal modes. Once the coordinate is identified, we find that evolution along it results in a change in diabatic state and optical excitation energy, as predicted by the Marcus model. Thus, we conclude that a single dimension of the electron transfer described in Marcus-Hush theory can be described as a well-defined nuclear motion. Importantly, our approach allows the separation of the intrinsic electron transfer coordinate from other structural relaxations and environmental influences. Furthermore, the barrier separating the adiabatic minima was found to be sufficiently thin to enable heavy-atom tunneling in the ET process.

The π-interactions of ammonia ligands evaluated by ab initio ligand field theory

Ammonia and amine ligands are commonly assumed to be σ-only ligands in coordination chemistry, i.e. they are not expected to interact significantly with a metal via a π path. Ligand field analyses employing the Angular Overlap Model resulted in good fits to experimental data without a π parameter for ammonia ligands, thereby supporting this assumption. In this work, we challenge this assumption and suggest that it is an oversimplification. We use complete active space calculations for electronic structure analyses of copper ammine complexes that are in good agreement with the transitions observed in experimental UV-vis spectra. These findings lead to a reinterpretation of the experimental spectra that necessitates a significant π interaction of the ammonia ligands. The strength of the ammonia π interaction is evaluated by parameterizing the ligand field splittings of a series of metal hexammine complexes ([M(NH₃)₆]ⁿ⁺ with M = Cr, Mn, Fe, Co, Ni, Ru, Os and n = 2, 3) and selected tetrammine complexes ([M(NH₃)₄]ⁿ⁺ with M = Cr, Mn, Fe, Co, Ni and n = 2 or 3) with the Angular Overlap Model. The resulting π parameters show that ammonia is a π donor of similar strength as chloride.

Fine-Tuning Redox Properties of Heteroleptic Molybdenum Complexes through Ligand-Ligand-Cooperativity

Heteroleptic molybdenum complexes bearing 1,5-diaza-3,7-diphosphacyclooctane (P2N2) and non-innocent dithiolene ligands were synthesized and electrochemically characterized. The reduction potentials of the complexes were found to be fine-tuned by a synergistic effect identified by DFT calculations as ligand-ligand cooperativity via non-covalent interactions. This finding is supported by electrochemical studies combined with UV-vis spectroscopy and temperature-dependent NMR spectroscopy. The observed behavior is reminiscent of enzymatic redox modulation using second ligand sphere effects.

How Metal Nuclearity Impacts Electrocatalytic H2 Production in Thiocarbohydrazone-Based Complexes

Thiocarbohydrazone-based catalysts feature ligands that are potentially electrochemically active. From the synthesis point of view, these ligands can be easily tailored, opening multiple strategies for optimization, such as using different substituent groups or metal substitution. In this work, we show the possibility of a new strategy, involving the nuclearity of the system, meaning the number of metal centers. We report the synthesis and characterization of a trinuclear nickel-thiocarbohydrazone complex displaying an improved turnover rate compared with its mononuclear counterpart. We use DFT calculations to show that the mechanism involved is metal-centered, unlike the metal-assisted ligand-centered mechanism found in the mononuclear complex. Finally, we show that two possible mechanisms can be assigned to this catalyst, both involving an initial double reduction of the system.

Substituent Effects in Iron Porphyrin Catalysts for the Hydrogen Evolution Reaction

For a future hydrogen economy, non-precious metal catalysts for the water splitting reactions are needed that can be implemented on a global scale. Metal-nitrogen-carbon (MNC) catalysts with active sites constituting a metal center with fourfold coordination of nitrogen (MN₄ ) show promising performance, but an optimization rooted in structure-property relationships has been hampered by their low structural definition. Porphyrin model complexes are studied to transfer insights from well-defined molecules to MNC systems. This work combines experiment and theory to evaluate the influence of porphyrin substituents on the electronic and electrocatalytic properties of MN₄ centers with respect to the hydrogen evolution reaction (HER) in aqueous electrolyte. We found that the choice of substituent affects their utilization on the carbon support and their electrocatalytic performance. We propose an HER mechanism for supported iron porphyrin complexes involving a [Fe^II (P⋅)]^- radical anion intermediate, in which a porphinic nitrogen atom acts as an internal base. While this work focuses on the HER, the limited influence of a simultaneous interaction with the support and an aqueous electrolyte will likely be transferrable to other catalytic applications.

An Iron‐Catalyzed Route to Dewar 1,3,5‐Triphosphabenzene and Subsequent Reactivity

The application of an alkyne cyclotrimerization regime with an [Fe(salen)]₂‐μ‐oxo (1) catalyst to triphenylmethylphosphaalkyne (2) yields gram‐scale quantities of 2,4,6‐tris(triphenylmethyl)‐Dewar‐1,3,5‐triphosphabenzene (3). Bulky lithium salt LiHMDS facilitates a rearrangement of 3 to the 1,3,5‐triphosphabenzene valence isomer (3′), which subsequently undergoes an intriguing phosphorus migration reaction to form the ring‐contracted species (3′′). Density functional theory calculations provide a plausible mechanism for this rearrangement. Given the stability of 3, a diverse array of unprecedented transformations was investigated. We report novel crystallographically characterized products of successful nucleophilic/electrophilic addition and protonation/oxidation reactions.

Identification of the Catalytically Dominant Iron Environment in Iron- and Nitrogen-Doped Carbon Catalysts for the Oxygen Reduction Reaction

For large-scale utilization of fuel cells in a future hydrogen-based energy economy, affordable and environmentally benign catalysts are needed. Pyrolytically obtained metal- and nitrogen-doped carbon (MNC) catalysts are key contenders for this task. Their systematic improvement requires detailed knowledge of the active site composition and degradation mechanisms. In FeNC catalysts, the active site is an iron ion coordinated by nitrogen atoms embedded in an extended graphene sheet. Herein, we build an active site model from in situ and operando ⁵⁷Fe Mössbauer spectroscopy and quantum chemistry. A Mössbauer signal newly emerging under operando conditions, D4, is correlated with the loss of other Mössbauer signatures (D2, D3a, D3b), implying a direct structural correspondence. Pyrrolic N-coordination, i.e., FeN₄C₁₂, is found as a spectroscopically and thermodynamically consistent model for the entire catalytic cycle, in contrast to pyridinic nitrogen coordination. These findings thus overcome the previously conflicting structural assignments for the active site and, moreover, identify and structurally assign a previously unknown intermediate in the oxygen reduction reaction at FeNC catalysts.

Photoelectrochemical Conversion of Dinitrogen to Benzonitrile: Selectivity Control by Electrophile‐ versus Proton‐Coupled Electron Transfer

Nitride complexes are key species in homogeneous nitrogen fixation to NH₃ via stepwise proton‐coupled electron transfer (PCET). In contrast, direct generation of nitrogenous organic products from N₂‐derived nitrides requires new strategies to enable efficient reductive nitride transfer in the presence of organic electrophiles. We here present a 2‐step protocol for the conversion of dinitrogen to benzonitrile. Photoelectrochemical, reductive N₂ splitting produces a rhenium(V) nitride with unfavorable PCET thermochemistry towards ammonia generation. However, N‐benzoylation stabilizes subsequent reduction as a basis for selective nitrogen transfer in the presence of the organic electrophile and Brønsted acid at mild reduction potentials. This work offers a new strategy for photoelectrosynthetic nitrogen fixation beyond ammonia—to yield nitrogenous organic products.

Revisiting the Fundamental Nature of Metal‐Ligand Bonding: An Impartial and Automated Fitting Procedure for Angular Overlap Model Parameters

The properties and reactivities of transition metal complexes are often discussed in terms of Ligand Field Theory (LFT), and with ab initio LFT a direct connection to quantum chemical wavefunctions was recently established. The Angular Overlap Model (AOM) is a widely used, ligand‐specific parameterization scheme of the ligand field splitting that has, however, been restricted by the availability and resolution of experimental data. Using ab initio LFT, we present here a generalised, symmetry‐independent and automated fitting procedure for AOM parameters that is even applicable to formally underdetermined or experimentally inaccessible systems. This method allows quantitative evaluations of assumptions commonly made in AOM applications, for example, transferability or the relative magnitudes of AOM parameters, and the response of the ligand field to structural or electronic changes. A two‐dimensional spectrochemical series of tetrahedral halido metalates ([M^IIX₄]²⁻, M=Mn−Cu) served as a case study. A previously unknown linear relationship between the halide ligands’ chemical hardness and their AOM parameters was found. The impartial and automated procedure for identifying AOM parameters introduced here can be used to systematically improve our understanding of ligand–metal interactions in coordination complexes.

Extreme -Tensor Anisotropy and Its Insensitivity to Structural Distortions in a Family of Linear Two-Coordinate Ni(I) Bis-N-heterocyclic Carbene Complexes

We report a new series of homoleptic Ni(I) bis-N-heterocyclic carbene complexes with a range of torsion angles between the two ligands from 68° to 90°. Electron paramagnetic resonance measurements revealed a strongly anisotropic g-tensor in all complexes with a small variation in g_∥ ∼ 5.7-5.9 and g_⊥ ∼ 0.6. The energy of the first excited state identified by variable-field far-infrared magnetic spectroscopy and SOC-CASSCF/NEVPT2 calculations is in the range 270-650 cm^-1. Magnetic relaxation measured by alternating current susceptibility up to 10 K is dominated by Raman and direct processes. Ab initio ligand-field analysis reveals that a torsion angle of <90° causes the splitting between doubly occupied d_xz and d_yz orbitals, which has little effect on the magnetic properties, while the temperature dependence of the magnetic relaxation appears to have no correlation with the torsion angle.

Reduction induced S-nucleophilicity in mono-dithiolene molybdenum complexes - in situ generation of sulfonium ligands

The reduction of a molybdenum monodithiolene complex, [Mo(CO)₂(dt)(dppe)], in the presence of dichloromethane leads to the transfer of CH₂ to sulfur and respective sulfonium species. Detailed analytical and mechanistical spectroscopic and electrochemical studies reveal the reasons for the unexpected formation and composition of the very unusual resultant complexes to be electronic-energetic in nature.

Picomolar FKBP inhibitors enabled by a single water-displacing methyl group in bicyclic [4.3.1] aza-amides

Methyl groups can have profound effects in drug discovery but the underlying mechanisms are diverse and incompletely understood. Here we report the stereospecific effect of a single, solvent-exposed methyl group in bicyclic [4.3.1] aza-amides, robustly leading to a 2 to 10-fold increase in binding affinity for FK506-binding proteins (FKBPs). This resulted in the most potent and efficient FKBP ligands known to date. By a combination of co-crystal structures, isothermal titration calorimetry (ITC), density-functional theory (DFT), and 3D reference interaction site model (3D-RISM) calculations we elucidated the origin of the observed affinity boost, which was purely entropically driven and relied on the displacement of a water molecule at the protein-ligand-bulk solvent interface. The best compounds potently occupied FKBPs in cells and enhanced bone morphogenic protein (BMP) signaling. Our results show how subtle manipulation of the solvent network can be used to design atom-efficient ligands for difficult, solvent-exposed binding pockets.

Structurally Characterized μ-1,2-Peroxo/Superoxo Dicopper(II) Pair

Superoxo complexes of copper are primary adducts in several O₂-activating Cu-containing metalloenzymes as well as in other Cu-mediated oxidation and oxygenation reactions. Because of their intrinsically high reactivity, however, isolation of Cu_x(O₂^•-) species is challenging. Recent work (J. Am. Chem. Soc. 2017, 139, 9831; 2019, 141, 12682) established fundamental thermochemical data for the H atom abstraction reactivity of dicopper(II) superoxo complexes, but structural characterization of these important intermediates was so far lacking. Here we report the first crystallographic structure determination of a superoxo dicopper(II) species (3) together with the structure of its 1e^- reduced peroxo congener (2; a rare cis-μ-1,2-peroxo dicopper(II) complex). Interconversion of 2 and 3 occurs at low potential (-0.58 V vs Fc/Fc⁺) and is reversible both chemically and electrochemically. Comparison of metric parameters (d(O-O) = 1.441(2) Å for 2 vs 1.329(7) Å for 3) and of spectroscopic signatures (ν̃(¹⁶O-¹⁶O) = 793 cm^-1 for 2 vs 1073 cm^-1 for 3) reflects that the redox process occurs at the bridging O₂-derived unit. The Cu^II-O₂^•--Cu^II complex has an S = 1/2 spin ground state according to magnetic and EPR data, in agreement with density functional theory calculations. Computations further show that the potential associated with changes of the Cu-O-O-Cu dihedral angle is shallow for both 2 and 3. These findings provide a structural basis for the low reorganization energy of the kinetically facile 1e^- interconversion of μ-1,2-superoxo/peroxo dicopper(II) couples, and they open the door for comprehensive studies of these key intermediates in Cu_x/O₂ chemistry.

Cyanate Formation via Photolytic Splitting of Dinitrogen

Light-driven N₂ cleavage into molecular nitrides is an attractive strategy for synthetic nitrogen fixation. However, suitable platforms are rare. Furthermore, the development of catalytic protocols via this elementary step suffers from poor understanding of N-N photosplitting within dinitrogen complexes, as well as of the thermochemical and kinetic framework for coupled follow-up chemistry. We here present a tungsten pincer platform, which undergoes fully reversible, thermal N₂ splitting and reverse nitride coupling, allowing for experimental derivation of thermodynamic and kinetic parameters of the N-N cleavage step. Selective N-N splitting was also obtained photolytically. DFT computations allocate the productive excitations within the {WNNW} core. Transient absorption spectroscopy shows ultrafast repopulation of the electronic ground state. Comparison with ground-state kinetics and resonance Raman data support a pathway for N-N photosplitting via a nonstatistically vibrationally excited ground state that benefits from vibronically coupled structural distortion of the core. Nitride carbonylation and release are demonstrated within a full synthetic cycle for trimethylsilylcyanate formation directly from N₂ and CO.

Spectroscopic discernibility of dopants and axial ligands in pyridinic FeN4 environments relevant to single-atom catalysts

Single-atom catalysts (SACs) activate small molecules, e.g. the oxygen reduction reaction is catalysed by FeNC materials. Because the nature of active site(s) in this type of SAC is unclear, spectroscopic and computational insights are needed to clarify the atomistic composition and electronic structure. Using quantum chemistry, we show that key features of [Fe{phen2A2}L]n+ complexes (A = CH, N with n = 0, A = O with n = 0, 2; L = OH-, Cl-) can be differentiated spectroscopically.

Orbital Entanglement Analysis of Exchange-Coupled Systems

A new tool for the interpretation of multiconfigurational wave functions representing the spin states of exchange-coupled transition metal complexes is introduced. Based on orbital entanglement measures, herein derived from multiconfigurational density matrix renormalization group calculations, the complexity of the wave function is reduced, thus facilitating a connection with established concepts for the interpretation of magnetically coupled systems. We show that the entanglement of localized orbitals with a small basis set is a good representation of the magnetic coupling topology and that it is sensitive to chemical changes in homologous complexes. Furthermore, we introduce a measure for the magnetic relevance of orbitals in the active subspace and a concept for the quantitative comparison of different chemical species. The approach presented here will be easily applicable to higher nuclearity clusters, providing a direct insight into all states of the Heisenberg spin ladder for systems previously accessible only by single-configurational methods.

Implications of structural heterogeneity for the electronic structure of the final oxygen-evolving intermediate in photosystem II

Heterogeneity in intermediate catalytic states of the oxygen-evolving complex (OEC) of Photosystem II is known from a wide range of experimental and theoretical data, but its potential implications for the mechanism of water oxidation remain unexplored. We delineate the consequences of structural heterogeneity for the final step of the catalytic cycle by tracing the evolution of three spectroscopically relevant and structurally distinct components of the last metastable S₃ state to the transient O₂-evolving S₄ state of the OEC. Using quantum chemical calculations, we show that each S₃ isomer leads to a different electronic structure formulation for the active S₄ state. Crucially, in addition to previously hypothesized Mn(IV)-oxyl species, we establish for the first time, how a genuine Mn(V)-oxo can be obtained in the catalytically active S₄ state: this takes the form of a five-coordinate and locally high-spin (S_Mn = 1) Mn(V) site. This formulation for the S₄ state evolves naturally from a preceding S₃-state structural intermediate that contains a quasi-trigonal-bipyramidal Mn(IV) ion. The results strongly suggest that water binding in the S₃ state is not prerequisite for reaching the oxygen-evolving S₄ state of the complex, supporting the notion that both substrates are preloaded at the beginning of the catalytic cycle. This scenario allows true four-electron metal-centered hole accumulation to precede OO bond formation and hence the latter can proceed via a genuine even-electron mechanism. This can occur as intramolecular nucleophilic coupling of two oxo units synchronously with the binding of a water substrate for the next catalytic cycle.

Steric Switching From Photochemical to Thermal N2 Splitting: A Computational Analysis of the Isomerization Reaction {(Cp*)(Am)Mo}2(μ-η1:η1-N2) → {(Cp*)(Am)Mo}2(μ-N)2

A μ-η1:η1-N2-bridged Mo dimer, {(η5-C5Me5)[N(Et)C(Ph)N(Et)]Mo}2(μ-N2), cleaves dinitrogen thermally resulting in a crystallographically characterized bis-μ-N-bridged dimer, {(η5-C5Me5)[N(Et)C(Ph)N(Et)]Mo}2(μ-N)2. A structurally related Mo dimer with a bulkier amidinate ligand, ([N(iPr)C(Me)N(iPr)]), is only capable of photochemical dinitrogen activation. These opposing reactivities were rationalized as steric switching between the thermally and photochemically active species. A computational analysis of the geometric and electronic structures of intermediates along the isomerization pathway from Mo2(μ-η1:η1-N2) to Mo2(μ-η2:η1-N2) and Mo2(μ-η2:η2-N2), and finally Mo2(μ-N)2, is presented here. The extent to which dispersion affects the thermodynamics of the isomers is evaluated, and it is found that dispersion interactions play a significant role in stabilizing the product and making the reaction exergonic. The concept of steric switching is further explored with theoretical models with sterically even less demanding ligands, indicating that systematic ligand modifications could be used to rationally design the N2 activation energy landscape. An analysis of electronic excitations in the computed UV-vis spectra of the two complexes shows that a particular type of asymmetric excitations is only present in the photoactive complex.

Mono- and dinuclear Ni(i) products formed upon bromide abstraction from the Ni(i) ring-expanded NHC complex [Ni(6-Mes)(PPh3)Br]

New T- and Y-shaped Ni(i) complexes are reported and analysed by DFT and EPR.

Dinitrogen photoactivation: status quo and future perspectives

New perspectives for dinitrogen activation: an overview of photochemical pathways to cleave the strong N–N bond.

Exchange Coupling Interactions from the Density Matrix Renormalization Group and N-Electron Valence Perturbation Theory: Application to a Biomimetic Mixed-Valence Manganese Complex

The accurate description of magnetic level energetics in oligonuclear exchange-coupled transition-metal complexes remains a formidable challenge for quantum chemistry. The density matrix renormalization group (DMRG) brings such systems for the first time easily within reach of multireference wave function methods by enabling the use of unprecedentedly large active spaces. But does this guarantee systematic improvement in predictive ability and, if so, under which conditions? We identify operational parameters in the use of DMRG using as a test system an experimentally characterized mixed-valence bis-μ-oxo/μ-acetato Mn(III,IV) dimer, a model for the oxygen-evolving complex of photosystem II. A complete active space of all metal 3d and bridge 2p orbitals proved to be the smallest meaningful starting point; this is readily accessible with DMRG and greatly improves on the unrealistic metal-only configuration interaction or complete active space self-consistent field (CASSCF) values. Orbital optimization is critical for stabilizing the antiferromagnetic state, while a state-averaged approach over all spin states involved is required to avoid artificial deviations from isotropic behavior that are associated with state-specific calculations. Selective inclusion of localized orbital subspaces enables probing the relative contributions of different ligands and distinct superexchange pathways. Overall, however, full-valence DMRG-CASSCF calculations fall short of providing a quantitative description of the exchange coupling owing to insufficient recovery of dynamic correlation. Quantitatively accurate results can be achieved through a DMRG implementation of second order N-electron valence perturbation theory (NEVPT2) in conjunction with a full-valence metal and ligand active space. Perspectives for future applications of DMRG-CASSCF/NEVPT2 to exchange coupling in oligonuclear clusters are discussed.

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.

Redox potential tuning by redox-inactive cations in nature's water oxidizing catalyst and synthetic analogues

Fundamental differences between synthetic manganese clusters and the biological water oxidizing catalyst are demonstrated in the modulation of their redox potential by redox-inactive cations.

A five-coordinate Mn(iv) intermediate in biological water oxidation: spectroscopic signature and a pivot mechanism for water binding

Among the four photo-driven transitions of the water-oxidizing tetramanganese-calcium cofactor of biological photosynthesis, the second-last step of the catalytic cycle, that is the S₂ to S₃ state transition, is the crucial step that poises the catalyst for the final O-O bond formation. This transition, whose intermediates are not yet fully understood, is a multi-step process that involves the redox-active tyrosine residue and includes oxidation and deprotonation of the catalytic cluster, as well as the binding of a water molecule. Spectroscopic data has the potential to shed light on the sequence of events that comprise this catalytic step, which still lacks a structural interpretation. In this work the S₂-S₃ state transition is studied and a key intermediate species is characterized: it contains a Mn₃O₄Ca cubane subunit linked to a five-coordinate Mn(iv) ion that adopts an approximately trigonal bipyramidal ligand field. It is shown using high-level density functional and multireference wave function calculations that this species accounts for the near-infrared absorption and electron paramagnetic resonance observations on metastable S₂-S₃ intermediates. The results confirm that deprotonation and Mn oxidation of the cofactor must precede the coordination of a water molecule, and lead to identification of a novel low-energy water binding mode that has important implications for the identity of the substrates in the mechanism of biological water oxidation.

How Accurately Can Extended X-ray Absorption Spectra Be Predicted from First Principles? Implications for Modeling the Oxygen-Evolving Complex in Photosystem II

First principle calculations of extended X-ray absorption fine structure (EXAFS) data have seen widespread use in bioinorganic chemistry, perhaps most notably for modeling the Mn4Ca site in the oxygen evolving complex (OEC) of photosystem II (PSII). The logic implied by the calculations rests on the assumption that it is possible to a priori predict an accurate EXAFS spectrum provided that the underlying geometric structure is correct. The present study investigates the extent to which this is possible using state of the art EXAFS theory. The FEFF program is used to evaluate the ability of a multiple scattering-based approach to directly calculate the EXAFS spectrum of crystallographically defined model complexes. The results of these parameter free predictions are compared with the more traditional approach of fitting FEFF calculated spectra to experimental data. A series of seven crystallographically characterized Mn monomers and dimers is used as a test set. The largest deviations between the FEFF calculated EXAFS spectra and the experimental EXAFS spectra arise from the amplitudes. The amplitude errors result from a combination of errors in calculated S0(2) and Debye-Waller values as well as uncertainties in background subtraction. Additional errors may be attributed to structural parameters, particularly in cases where reliable high-resolution crystal structures are not available. Based on these investigations, the strengths and weaknesses of using first-principle EXAFS calculations as a predictive tool are discussed. We demonstrate that a range of DFT optimized structures of the OEC may all be considered consistent with experimental EXAFS data and that caution must be exercised when using EXAFS data to obtain topological arrangements of complex clusters.

Ca K-edge XAS as a probe of calcium centers in complex systems

Herein, Ca K-edge X-ray absorption spectroscopy (XAS) is developed as a means to characterize the local environment of calcium centers. The spectra for six, seven, and eight coordinate inorganic and molecular calcium complexes were analyzed and determined to be primarily influenced by the coordination environment and site symmetry at the calcium center. The experimental results are closely correlated to time-dependent density functional theory (TD-DFT) calculations of the XAS spectra. The applicability of this methodology to complex systems was investigated using structural mimics of the oxygen-evolving complex (OEC) of PSII. It was found that Ca K-edge XAS is a sensitive probe for structural changes occurring in the cubane heterometallic cluster due to Mn oxidation. Future applications to the OEC are discussed.

Metal oxidation states in biological water splitting

A central question in biological water splitting concerns the oxidation states of the manganese ions that comprise the oxygen-evolving complex of photosystem II.

A central question in biological water splitting concerns the oxidation states of the manganese ions that comprise the oxygen-evolving complex of photosystem II. Understanding the nature and order of oxidation events that occur during the catalytic cycle of five S_i states (i = 0–4) is of fundamental importance both for the natural system and for artificial water oxidation catalysts. Despite the widespread adoption of the so-called “high-valent scheme”—where, for example, the Mn oxidation states in the S₂ state are assigned as III, IV, IV, IV—the competing “low-valent scheme” that differs by a total of two metal unpaired electrons (i.e. III, III, III, IV in the S₂ state) is favored by several recent studies for the biological catalyst. The question of the correct oxidation state assignment is addressed here by a detailed computational comparison of the two schemes using a common structural platform and theoretical approach. Models based on crystallographic constraints were constructed for all conceivable oxidation state assignments in the four (semi)stable S states of the oxygen evolving complex, sampling various protonation levels and patterns to ensure comprehensive coverage. The models are evaluated with respect to their geometric, energetic, electronic, and spectroscopic properties against available experimental EXAFS, XFEL-XRD, EPR, ENDOR and Mn K pre-edge XANES data. New 2.5 K ⁵⁵Mn ENDOR data of the S₂ state are also reported. Our results conclusively show that the entire S state phenomenology can only be accommodated within the high-valent scheme by adopting a single motif and protonation pattern that progresses smoothly from S₀ (III, III, III, IV) to S₃ (IV, IV, IV, IV), satisfying all experimental constraints and reproducing all observables. By contrast, it was impossible to construct a consistent cycle based on the low-valent scheme for all S states. Instead, the low-valent models developed here may provide new insight into the over-reduced S states and the states involved in the assembly of the catalytically active water oxidizing cluster.

The protonation states of oxo-bridged Mn(IV) dimers resolved by experimental and computational Mn K pre-edge X-ray absorption spectroscopy

In nature, the protonation of oxo bridges is a commonly encountered mechanism for fine-tuning chemical properties and reaction pathways. Often, however, the protonation states are difficult to establish experimentally. This is of particular importance in the oxygen evolving complex of photosystem II, where identification of the bridging oxo protonation states is one of the essential requirements toward unraveling the mechanism. In order to establish a combined experimental and theoretical protocol for the determination of protonation states, we have systematically investigated a series of Mn model complexes by Mn K pre-edge X-ray absorption spectroscopy. An ideal test case for selective bis-μ-oxo-bridge protonation in a Mn dimer is represented by the system [Mn(IV)2(salpn)2(μ-OHn)2](n+). Although the three species [Mn(IV)2(salpn)2(μ-O)2], [Mn(IV)2(salpn)2(μ-O)(μ-OH)](+) and [Mn(IV)2(salpn)2(μ-OH)2](2+) differ only in the protonation of the oxo bridges, they exhibit distinct differences in the pre-edge region while maintaining the same edge energy. The experimental spectra are correlated in detail to theoretically calculated spectra. A time-dependent density functional theory approach for calculating the pre-edge spectra of molecules with multiple metal centers is presented, using both high spin (HS) and broken symmetry (BS) electronic structure solutions. The most intense pre-edge transitions correspond to an excitation of the Mn 1s core electrons into the unoccupied orbitals of local e(g) character (d(z)(2) and d(xy) based in the chosen coordinate system). The lowest energy experimental feature is dominated by excitations of 1s-α electrons, and the second observed feature is primarily attributed to 1s-β electron excitations. The observed energetic separation is due to spin polarization effects in spin-unrestricted density functional theory and models final state multiplet effects. The effects of spin polarization on the calculated Mn K pre-edge spectra, in both the HS and BS solutions, are discussed in terms of the strength of the antiferromagnetic coupling and associated changes in the covalency of Mn-O bonds. The information presented in this paper is complemented with the X-ray emission spectra of the same compounds published in an accompanying paper. Taken together, the two studies provide the foundation for a better understanding of the X-ray spectroscopic data of the oxygen evolving complex (OEC) in photosystem II.