Mechanistic elucidation of O2 production from tBuOOH in water using the Mn(II) catalyst [Mn2(mcbpen)2(H2O)2]2+: a DFT study

This study employs density functional theory at the SMD/B3LYP-D3/6-311+G(2d,p),def2-TZVPP//SMD/B3LYP-D3/6-31G(d),SDD level of theory to explore the mechanistic details of O2 generation from tBuOOH, using H218O as the solvent, in the presence of the Mn(II) catalyst [Mn2(mcbpen)2(H2O)2]2+. Since this chemistry was reported to occur through the reaction of Mn(III)(μ-O)Mn(IV)-O˙ with water, we first revaluated this proposal and found that it occurs with an activation barrier greater than 36 kcal mol-1, ruling out the functioning of such a dimer as the active catalyst. Experimental evidence has shown that the oxidation of [Mn2(mcbpen)2(H2O)2]2+ by tBuOOH in H218O produces the Mn(IV) species [Mn(18O)(mcbpen)]+. Our investigations revealed a plausible mechanism for this observation in which [Mn (18O)(mcbpen)]+ acts as the active catalyst, generating the tert-butyl peroxyl radical (tBuOO˙) through its reaction with tBuOOH. In this proposed mechanism, the O-O bond is formed through the interaction of tBuOO˙ with another [Mn(18O)(mcbpen)]+, finally leading to the formation of the 16O18O product. Our findings underscore the pivotal role of [Mn(18O)(mcbpen)]+ in both generating the active species tBuOO˙ and consuming it to produce 16O18O. With activation barriers as low as about 9 kcal mol-1, these elementary steps highlight the feasibility of our proposed mechanism. Moreover, this mechanism elucidates why, experimentally, one of the oxygen atoms in the released O2 comes from water, while the other originates from tBuOOH. This research broadens our understanding of high oxidation state manganese chemistry, setting the stage for the development of more efficient Mn-based catalysts, aimed at improving processes in both renewable energy and synthetic chemistry.

Elucidating the catalytic mechanisms of O2 generation by [Mn2(μ-O)2(terpy)2(OH2)2]3+ using DFT calculations: a focus on ClO- as oxidant

The experimentally reported Mn(IV)Mn(III) complex [Mn2(μ-O)2(terpy)2(OH2)2]3+ has been observed catalyzing O2 generation with oxidants like ClO- and HSO5-. Previous mechanistic studies primarily focused on O2 generation with HSO5-, concluding that Mn(IV)Mn(III) acts as a catalyst, generating a Mn(IV)Mn(IV)-oxyl species as a key intermediate responsible for O-O bond formation. This computational study employs DFT calculations to investigate whether the catalytic generation of O2 using ClO- follows the same mechanism previously identified with HSO5- as the oxidant, or if it proceeds through an alternate pathway. To this end, we explored multiple pathways using ClO- as the oxidant. Interestingly, our findings confirm that in the case of ClO- as the oxidant, similar to what was observed with HSO5-, the Mn(IV)Mn(IV)-oxyl species indeed plays a crucial role in driving the catalytic evolution of O2 with the potential formation of the binuclear complexes Mn(IV)Mn(IV)-oxy and Mn(IV)Mn(IV)-OH during the reaction. These complexes are reactive in producing O2, with activation free energies of 15.9 and 14.3 kcal mol-1, respectively. However, our calculations revealed that the Mn(IV)Mn(IV)-oxyl complex is significantly more reactive in producing O2 than Mn(IV)Mn(IV)-oxy and Mn(IV)Mn(IV)-OH, with a lower free energy barrier of 8.1 kcal mol-1. Consequently, even though Mn(IV)Mn(IV)-oxyl is predicted to be present in much lower concentrations than Mn(IV)Mn(IV)-oxy and Mn(IV)Mn(IV)-OH, it emerges as the species acting as the active catalyst for catalytic O2 generation. This study enhances our knowledge of high oxidation state (+3 and +4) manganese chemistry, highlighting its key role in catalysis and paving the way for more efficient Mn-based catalysts with broad applications.

A Review of Recent Developments in Molecular Dynamics Simulations of the Photoelectrochemical Water Splitting Process

In this review, we provide a short overview of the Molecular Dynamics (MD) method and how it can be used to model the water splitting process in photoelectrochemical hydrogen production. We cover classical non-reactive and reactive MD techniques as well as multiscale extensions combining classical MD with quantum chemical and continuum methods. Selected examples of MD investigations of various aqueous semiconductor interfaces with a special focus on TiO2 are discussed. Finally, we identify gaps in the current state-of-the-art where further developments will be needed for better utilization of MD techniques in the field of water splitting.

Quantum Chemical Modeling of Homogeneous Water Oxidation Catalysis

The design of efficient and robust water oxidation catalysts has proven challenging in the development of artificial photosynthetic systems for solar energy harnessing and storage. Tremendous progress has been made in the development of homogeneous transition-metal complexes capable of mediating water oxidation. To improve the efficiency of the catalyst and to design new catalysts, a detailed mechanistic understanding is necessary. Quantum chemical modeling calculations have been successfully used to complement the experimental techniques to suggest a catalytic mechanism and identify all stationary points, including transition states for both O-O bond formation and O₂ release. In this review, recent progress in the applications of quantum chemical methods for the modeling of homogeneous water oxidation catalysis, covering various transition metals, including manganese, iron, cobalt, nickel, copper, ruthenium, and iridium, is discussed.

First-Principles Molecular Dynamics Studies of Organometallic Complexes and Homogeneous Catalytic Processes

Computational chemistry is a valuable aid to complement experimental studies of organometallic systems and their reactivity. It allows probing mechanistic hypotheses and investigating molecular structures, shedding light on the behavior and properties of molecular assemblies at the atomic scale. When approaching a chemical problem, the computational chemist has to decide on the theoretical approach needed to describe electron/nuclear interactions and the composition of the model used to approximate the actual system. Both factors determine the reliability of the modeling study. The community dedicated much effort to developing and improving the performance and accuracy of theoretical approaches for electronic structure calculations, on which the description of (inter)atomic interactions rely. Here, the importance of the model system used in computational studies is highlighted through examples from our recent research focused on organometallic systems and homogeneous catalytic processes. We show how the inclusion of explicit solvent allows the characterization of molecular events that would otherwise not be accessible in reduced model systems (clusters). These include the stabilization of nascent charged fragments via microscopic solvation (notably, hydrogen bonding), transfer of charge (protons) between distant fragments mediated by solvent molecules, and solvent coordination to unsaturated metal centers. Furthermore, when weak interactions are involved, we show how conformational and solvation properties of organometallic complexes are also affected by the explicit inclusion of solvent molecules. Such extended model systems may be treated under periodic boundary conditions, thus removing the cluster/continuum (or vacuum) boundary, and require a statistical mechanics simulation technique to sample the accessible configurational space. First-principles molecular dynamics, in which atomic forces are computed from electronic structure calculations (namely, density functional theory), is certainly the technique of choice to investigate chemical events in solution. This methodology is well established and thanks to advances in both algorithms and computational resources simulation times required for the modeling of chemical events are nowadays accessible, though the computational requirements use to be high. Specific applications reviewed here include mechanistic studies of the Shilov and Wacker processes, speciation in Pd chemistry, hydrogen bonding to metal centers, and the dynamics of agostic interactions.