Computational electrochemistry: prediction of liquid-phase reduction potentials

High-valent metal-oxo species transformation and regulation by co-existing chloride: Reaction pathways and impacts on the generation of chlorinated by-products

High-valent metal-oxo species (HMOS) have been extensively recognized in advanced oxidation processes (AOPs) owing to their high selectivity and high chemical utilization efficiency. However, the interactions between HMOS and halide ions in sewage wastewater are complicated, leading to ongoing debates on the intrinsic reactive species and impacts on remediation. Herein, we prepared three typical HMOS, including Fe(IV), Mn(V)-nitrilotriacetic acid complex (Mn(V)NTA) and Co(IV) through peroxymonosulfate (PMS) activation and comparatively studied their interactions with Cl- to reveal different reactive chlorine species (RCS) and the effects of HMOS types on RCS generation pathways. Our results show that the presence of Cl- alters the cleavage behavior of the peroxide OO bond in PMS and prohibits the generation of Fe(IV), spontaneously promoting SO4•- production and its subsequent transformation to secondary radicals like Cl• and Cl2•-. The generation and oxidation capacity of Mn(V)NTA was scarcely influenced by Cl-, while Cl- would substantially consume Co(IV) and promote HOCl generation through an oxygen-transfer reaction, evidenced by density functional theory (DFT) and deuterium oxide solvent exchange experiment. The two-electron-transfer standard redox potentials of Fe(IV), Mn(V)NTA and Co(IV) were calculated as 2.43, 2.55 and 2.85 V, respectively. Due to the different reactive species and pathways in the presence of Cl-, the amounts of chlorinated by-products followed the order of Co(II)/PMS > Fe(II)/PMS > Mn(II)NTA/PMS. Thus, this work renovates the knowledge of halide chemistry in HMOS-based systems and sheds light on the impact on the treatment of salinity-containing wastewater.

Limiting factors in the accuracy of DFT calculation for redox potentials

In the present study, we have investigated factors affecting the accuracy of computational chemistry calculation of redox potentials, namely the gas-phase ionization energy (IE) and electron affinity (EA), and the continuum solvation effect. In general, double-hybrid density functional theory methods yield IEs and EAs that are on average within ~0.1 eV of our high-level W3X-L benchmark, with the best performing method being DSD-BLYP/ma-def2-QZVPP. For lower-cost methods, the average errors are ~0.2-0.3 eV, with ωB97X-3c being the most accurate (~0.15 eV). For the solvation component, essentially all methods have an average error of ~0.3 eV, which shows the limitation of the continuum solvation model. Curiously, the directly calculated redox potentials show errors of ~0.3 eV for all methods. These errors are notably smaller than what can be expected from error propagation with the two components (IE and EA, and solvation effect). Such a discrepancy can be attributed to the cancellation of errors, with the lowest-cost GFN2-xTB method benefiting the most, and the most accurate ωB97X-3c method benefiting the least. For organometallic species, the redox potentials show large deviations exceeding ~0.5 eV even for DSD-BLYP. The large errors are attributed to those for the gas-phase IEs and EAs, which represents a major barrier to the accurate calculation of redox potentials for such systems.

Development of a multi-step screening procedure for redox active molecules in organic radical polymer anodes and as redox flow anolytes

Benzo[d]-X-zolyl-pyridinyl (XO, S, NH) radicals represent a promising class of redox-active molecules for organic batteries. We present a multistep screening procedure to identify the most promising radical candidates. Experimental investigations and highly correlated wave function-based calculations are performed to determine benchmark redox potentials. Based on these, the accuracies of different methods (semi-empirical, density functional theory, wave function-based), solvent models, dispersion corrections, and basis sets are evaluated. The developed screening procedure consists of three steps: First, a conformer search is performed with CREST. The molecules are selected based on the redox potentials calculated using GFN2-xTB. Second, HOMO energies calculated with reparametrized B3LYP-D3(BJ) and the def2-SVP basis set are used as selection criteria. The final molecules are selected based on the redox potentials calculated from Gibbs energies using BP86-D3(BJ)/def2-TZVP. With this multistep screening approach, promising molecules can be suggested for synthesis, and structure-property relationships can be derived.

Redox potentials in ionic liquids: Anomalous behavior?

Redox potentials depend on the nature of the solvent/electrolyte through the solvation energies of the ionic solute species. For concentrated electrolytes, ion solvation may deviate significantly from the Born model predictions due to ion pairing and correlation effects. Recently, Ghorai and Matyushov [J. Phys. Chem. B 124, 3754-3769 (2020)] predicted, on the basis of linear response theory, an anomalous trend in the solvation energies of room temperature ionic liquids, with deviations of hundreds of kJ/mol from the Born model for certain size solutes/ions. In this work, we computationally evaluate ionic solvation energies in the prototypical ionic liquid, 1-butyl-3-methylimidazolium tetrafluoroborate (BMIM/BF4), to further explore this behavior and benchmark several of the approximations utilized in the solvation energy predictions. For comparison, we additionally compute solvation energies within acetonitrile and molten NaCl salt to illustrate the limiting behavior of purely dipolar and ionic solvents. We find that the overscreening effect, which results from the inherent charge oscillations of the ionic liquid, is substantially reduced in magnitude due to screening from the dipoles of the molecular ions. Therefore, for the molten NaCl salt, for which the ions do not have permanent dipoles, modulation of ionic solvation energies from the overscreening effect is most significant. The conclusion is that ionic liquids do indeed exhibit unique solvation behavior due to peak(s) in the electrical susceptibility caused by the ion shell structure; redox potential shifts for BMIM/BF4 are of more modest order ∼0.1 V, but may be larger for other ionic liquids that approach molten salt behavior.

Machine Learning-Based Prediction of Reduction Potentials for PtIV Complexes

Some of the well-known drawbacks of clinically approved PtII complexes can be overcome using six-coordinate PtIV complexes as inert prodrugs, which release the corresponding four-coordinate active PtII species upon reduction by cellular reducing agents. Therefore, the key factor of PtIV prodrug mechanism of action is their tendency to be reduced which, when the involved mechanism is of outer-sphere type, is measured by the value of the reduction potential. Machine learning (ML) models can be used to effectively capture intricate relationships within PtIV complex data, leading to highly accurate predictions of reduction potentials and other properties, and offering significant insights into their electrochemical behavior and potential applications. In this study, a machine learning-based approach for predicting the reduction potentials of PtIV complexes based on relevant molecular descriptors is presented. Leveraging a data set of experimentally determined reduction potentials and a diverse range of molecular descriptors, the proposed model demonstrates remarkable predictive accuracy (MSE = 0.016 V2, RMSE = 0.13 V, R2 = 0.92). Ab initio calculations and a set of different machine learning algorithms and feature engineering techniques have been employed to systematically explore the relationship between molecular structure and similarity and reduction potential. Specifically, it has been investigated whether the reduction potential of these compounds can be described by combining ML models across different combinations of constitutional, topological, and electronic molecular descriptors. Our results not only provide insights into the crucial factors influencing reduction potentials but also offer a rapid and effective tool for the rational design of PtIV complexes with tailored electrochemical properties for pharmaceutical applications. This approach has the potential to significantly expedite the development and screening of novel PtIV prodrug candidates. The analysis of principal components and key features extracted from the model highlights the significance of structural descriptors of the 2D Atom Pairs type and the lowest unoccupied molecular orbital energy. Specifically, with just 20 appropriately selected descriptors, a notable separation of complexes based on their reduction potential value is achieved.

Exploring the Cooperation of the Redox Non-Innocent Ligand and Di-Cobalt Center for the Water Oxidation Reaction Catalyzed by a Binuclear Complex

Water oxidation is a crucial reaction in the artificial photosynthesis system. In the present work, density functional calculations were employed to decipher the mechanism of water oxidation catalyzed by a binuclear cobalt complex, which was disclosed to be a homogeneous water oxidation catalyst in pH=7 phosphate buffer. The calculations showed that the catalytic cycle starts from the CoIII,III-OH2 species. Then, a proton-coupled electron transfer followed by a one-electron transfer process leads to the generation of the formal CoIV,IV-OH intermediate. The subsequent PCET produces the active species, namely the formal CoIV,V=O intermediate (4). The oxidation processes mainly occur on the ligand moiety, including the coordinated water moiety, implying a redox non-innocent behavior. Two cobalt centers keep their oxidation states and provide one catalytic center for water activation during the oxidation process. 4 triggers the O-O bond formation via the water nucleophilic attack pathway, in which the phosphate buffer ion functions as the proton acceptor. The O-O bond formation is the rate-limiting step with a calculated total barrier of 17.7 kcal/mol. The last electron oxidation process coupled with an intramolecular electron transfer results in the generation of O2.

Mesoporous Single Atom-Cluster Fe-N/C Oxygen Evolution Electrocatalysts Synthesized with Bottlebrush Block Copolymer-Templated Rapid Thermal Annealing

Current electrocatalysts for oxygen evolution reaction (OER) are either expensive (such as IrO2, RuO2) or/and exhibit high overpotential as well as sluggish kinetics. This article reports mesoporous earth-abundant iron (Fe)-nitrogen (N) doped carbon electrocatalysts with iron clusters and closely surrounding Fe-N4 active sites. Unique to this work is that the mechanically stable mesoporous carbon-matrix structure (79 nm in pore size) with well-dispersed nitrogen-coordinated Fe single atom-cluster is synthesized via rapid thermal annealing (RTA) within only minutes using a self-assembled bottlebrush block copolymer (BBCP) melamine-formaldehyde resin composite template. The resulting porous structure and domain size can be tuned with the degree of polymerization of the BBCP backbone, which increases the electrochemically active surface area and improves electron transfer and mass transport for an effective OER process. The optimized electrocatalyst shows a required potential of 1.48 V (versus RHE) to obtain the current density of 10 mA/cm2 in 1 M KOH aqueous electrolyte and a small Tafel slope of 55 mV/decade at a given overpotential of 250 mV, which is significantly lower than recently reported earth-abundant electrocatalysts. Importantly, the Fe single-atom nitrogen coordination environment facilitates the surface reconstruction into a highly active oxyhydroxide under OER conditions, as revealed by X-ray photoelectron spectroscopy and in situ Raman spectroscopy, while the atomic clusters boost the single atoms reactive sites to prevent demetalation during the OER process. Density functional theory (DFT) calculations support that the iron nitrogen environment and reconstructed oxyhydroxides are electrocatalytically active sites as the kinetics barrier is largely reduced. This work has opened a new avenue for simple, rapid synthesis of inexpensive, earth-abundant, tailorable, mechanically stable, mesoporous carbon-coordinated single-atom electrocatalysts that can be used for renewable energy production.

Pivotal Role of Geometry Regulation on O-O Bond Formation Mechanism of Bimetallic Water Oxidation Catalysts

In this study, we highlight the impact of catalyst geometry on the formation of O-O bonds in Cu2 and Fe2 catalysts. A series of Cu2 complexes with diverse linkers are designed as electrocatalysts for water oxidation. Interestingly, the catalytic performance of these Cu2 complexes is enhanced as their molecular skeletons become more rigid, which contrasts with the behavior observed in our previous investigation with Fe2 analogs. Moreover, mechanistic studies reveal that the reactivity of the bridging O atom results in distinct pathways for O-O bond formation in Cu2 and Fe2 catalysts. In Cu2 systems, the coupling takes place between a terminal CuIII -OH and a bridging μ-O⋅ radical. Whereas in Fe2 systems, it involves the coupling of two terminal Fe-oxo entities. Furthermore, an in-depth structure-activity analysis uncovers the spatial geometric prerequisites for the coupling of the terminal OH with the bridging μ-O⋅ radical, ultimately leading to the O-O bond formation. Overall, this study emphasizes the critical role of precisely adjusting the spatial geometry of catalysts to align with the O-O bonding pathway.

Divergent Bimetallic Mechanisms in Copper(II)-Mediated C-C, N-N, and O-O Oxidative Coupling Reactions

Copper-catalyzed aerobic oxidative coupling of diaryl imines provides a route for conversion of ammonia to hydrazine. The present study uses experimental and density functional theory computational methods to investigate the mechanism of N-N bond formation, and the data support a mechanism involving bimolecular coupling of Cu-coordinated iminyl radicals. Computational analysis is extended to CuII-mediated C-C, N-N, and O-O coupling reactions involved in the formation of cyanogen (NC-CN) from HCN, 1,3-butadiyne from ethyne (i.e., Glaser coupling), hydrazine from ammonia, and hydrogen peroxide from water. The results reveal two different mechanistic pathways. Heteroatom ligands with an uncoordinated lone pair (iminyl, NH2, OH) undergo charge transfer to CuII, generating ligand-centered radicals that undergo facile bimolecular radical-radical coupling. Ligands lacking a lone pair (CN and CCH) form bridged binuclear diamond-core structures that undergo C-C coupling. This mechanistic bifurcation is rationalized by analysis of spin densities in key intermediates and transition states, as well as multiconfigurational calculations. Radical-radical coupling is especially favorable for N-N coupling owing to energetically favorable charge transfer in the intermediate and thermodynamically favorable product formation.

Study of the microscopic mechanism of stepwise charge injection in co-sensitive DSSCs in the framework of a D-π-A dye and chlorophyll

The newly synthesized dye molecules TY6 and CXC22 were selected to explain the influence of anthracene and acetylene groups on the power conversion efficiency (PCE) of the molecules at the microscopic level. Theoretical simulation was carried out to understand the properties of the two molecules, including frontier molecular orbitals, absorption spectra, light absorption efficiency, intramolecular charge transfer (ICT), dye regeneration, I-V prediction, etc. The results suggest that for CXC22, adding an anthracene and acetylene group in the conjugate bridge greatly enhances the molecule's absorption wavelength and molar extinction coefficient; CXC22 also has significant advantages in the intramolecular charge transfer and comparatively better dye regeneration and electron injection. These parameters cause CXC22 to have a higher PCE. Subsequently, CXC22 and the chlorophyll molecule (CHL7) were selected for co-sensitization to regulate performance. The stable structure in the co-sensitization configuration was screened, and the absorption spectrum characteristics and charge transfer mechanisms were revealed for the co-sensitization system. The designed evaluation model predicted that the PCE of CO1 (the cosensitive system of CXC22 and TY6 in H-H configuration is referred to as CO1) could reach 16.78%. This work provides an idea for developing an efficient dye-sensitized solar cell system.

Mechanism of photocatalytic CO2 reduction to HCO2H by a robust multifunctional iridium complex

The tetradentate PNNP-type IrIII complex Mes-IrPCY2 ([Cl-IrIII-H]+) is reported to be an efficient catalyst for the reduction of CO2 to formate with excellent selectivity under visible light irradiation. Density functional calculations have been carried out to elucidate the mechanism and the origin of selectivity in the present work. Calculations suggest that the double-reduced complex 1-H (1[IrI-H]0) demonstrates higher activity than the single-reduced complex 2-H (2[IrIII(L˙-)-H]+), possibly owing to the higher hydride donor ability of the former compared to the latter; thus 1-H functions as the active species in the overall CO2 reduction reaction. In the HCOO- formation pathway, the hydride of 1-H performs a nucleophilic attack on CO2via an outer-sphere fashion to generate species 1-OCHO (1[IrI-OCHO]0), which then releases HCOO- to produce an IrI intermediate. A subsequent protonation and chloride coordination of the Ir center leads to the regeneration of catalyst 1[Cl-IrIII-H]+. For the CO production, a nucleophilic attack on CO2 takes place by the Ir atom of 1-Hvia an inner-sphere manner to afford complex O2C-3-H (1[O2C-IrIII-H]0), followed by a two-proton-one-electron reduction to furnish the OC-2-H complex (2[OC-IrIII(L˙-)-H]+) after liberating a H2O. Ultimately, CO is released to form 2-H. The stronger nucleophilicity as well as smaller steric hindrance of the hydride than the Ir atom of the active species 1-H (1[IrI-H]0) is found to account for the favoring of formate formation over CO formation. Meanwhile, the CO2 reduction reaction is calculated to be preferred over the hydrogen evolution reaction, and this is consistent with the experimental product distributions.

Prediction of Aqueous Reduction Potentials of X•, ChH•, and XO• Radicals with X = Halogen and Ch = Chalcogen

The aqueous electron affinity and aqueous reduction potentials for F•, Cl•, Br•, I•, OH•, SH•, SeH•, TeH•, ClO•, BrO•, and IO• were calculated using electronic structure methods for explicit cluster models coupled with a self-consistent reaction field (SMD) to treat the aqueous solvent. Calculations were conducted using MP2 and correlated molecular orbital theory up to the CCSD(T)-F12b level for water tetramer clusters and MP2 for octamer cluster. Inclusion of explicit waters was found to be important for accurately predicting the redox potentials in a number of cases. The calculated reduction potentials for X• and ChH• were predicted to within ∼0.1 V of the reported literature values. Fluorine is anomalous due to abstraction of a hydrogen from one of the surrounding water molecules to form a hydroxyl radical and hydrogen fluoride, so its redox potential was calculated using only an implicit model. Larger deviations from experiment were predicted for ClO• and BrO•. These deviations are due to the free energy of solvation of the anion being too negative, as found in the pKa calculations, and that for the neutral being too positive with the current approach.

One-Electron Oxidation Potentials and Hole Delocalization in Heterogeneous Single-Stranded DNA

The study of DNA processes is essential to understand not only its intrinsic biological functions but also its role in many innovative applications. The use of DNA as a nanowire or electrochemical biosensor leads to the need for a deep investigation of the charge transfer process along the strand as well as of the redox properties. In this contribution, the one-electron oxidation potential and the charge delocalization of the hole formed after oxidation are computationally investigated for different heterogeneous single-stranded DNA strands. We have established a two-step protocol: (i) molecular dynamics simulations in the frame of quantum mechanics/molecular mechanics (QM/MM) were performed to sample the conformational space; (ii) energetic properties were then obtained within a QM1/QM2/continuum approach in combination with the Marcus theory over an ensemble of selected geometries. The results reveal that the one-electron oxidation potential in the heterogeneous strands can be seen as a linear combination of that property within the homogeneous strands. In addition, the hole delocalization between different nucleobases is, in general, small, supporting the conclusion of a hopping mechanism for charge transport along the strands. However, charge delocalization becomes more important, and so does the tunneling mechanism contribution, when the reducing power of the nucleobases forming the strand is similar. Moreover, charge delocalization is slightly enhanced when there is a correlation between pairs of some of the interbase coordinates of the strand: twist/shift, twist/slide, shift/slide, and rise/tilt. However, the internal structure of the strand is not the predominant factor for hole delocalization but the specific sequence of nucleotides that compose the strand.

Toward a molecular understanding of the conductivity of lithium-ion conducting polyanion polymer electrolytes by molecular dynamics simulation

With the improved lithium-ion transference number near unity, the low conductivity of single lithium-ion conducting solid polymer electrolytes (SLIC-SPEs) still hinders their application in high-rate batteries. Though some empirical conclusions on the conducting mechanism of SLIC-SPEs have been obtained, a more comprehensive study on the quantitative relationship between the molecular structure factors and ionic conduction performance is expected. In this study, a model structure that contains adjustable main chain and anion groups in the polyethylene oxide (PEO) matrix was used to clarify the influence of molecular structural factors on ionic conductivity and electrochemical stability of SLIC-SPEs. The anionic group was further disassembled into the intermediate group and end group while the main chain structure was distinguished into different degrees of polymerization and various lengths of the spacers between anions. Therefore, a well-defined molecular structure was employed to describe its relationship with ionic conductivity. In addition, the dissociation degree of salts and mobility of ions changing with the molecular structure were also discussed to explore the fundamental causes of conductivity. It can be concluded that the anion group affects the conductivity mainly via the dissociation degree, while the main chain structure impacts the conductivity by both dissociation degree and mobility.

Tailoring Ag Electron Donating Ability for Organohalide Reduction: A Bilayer Electrode Design

Electrochemical reduction of organohalides provides a green approach in the reduction of environmental pollutants, the synthesis of new organic molecules, and many other applications. The presence of a catalytic electrode can make the process more energetically efficient. Ag is known to be a very good electrode for the reduction of a wide range of organohalides. Herein, we examine the elementary adsorption and reaction steps that occur on Ag and the changes that result from changes in the Ag-coated metal, strain in Ag, solvent, and substrate geometry. The results are used to develop an electrode design strategy that can possibly be used to further increase the catalytic activity of pure Ag electrodes. We have shown how epitaxially depositing one to three layers of Ag on catalytically inert or less active support metal can increase the surface electron donating ability, thus increasing the adsorption of organic halide and the catalytic activity. Many factors, such as molecular geometry, lattice mismatches, work function, and solvents, contribute to the adsorption of organic halide molecules over the bilayer electrode surface. To isolate and rank these factors, we examined three model organic halides, namely, halothane, bromobenzene (BrBz), and benzyl bromide (BzBr) adsorption on Ag/metal (metal = Au, Bi, Pt, and Ti) bilayer electrodes in both vacuum and acetonitrile (ACN) solvent. The different metal supports offer a range of lattice mismatches and work function differences with Ag. Our calculations show that the surface of Ag becomes more electron donating and accessible to adsorption when it forms a bilayer with Ti as it has a lower work function and almost zero lattice mismatch with Ag. We believe this study will help to increase the electron donating ability of the Ag surface by choosing the right metal support, which in turn can improve the catalytic activity of the working electrode.

Toward Benchmark-Quality Ab Initio Predictions for 3d Transition Metal Electrocatalysts: A Comparison of CCSD(T) and ph-AFQMC

Generating accurate ab initio ionization energies for transition metal complexes is an important step toward the accurate computational description of their electrocatalytic reactions. Benchmark-quality data is required for testing existing theoretical methods and developing new ones but is complicated to obtain for many transition metal compounds due to the potential presence of both strong dynamical and static electron correlation. In this regime, it is questionable whether the so-called gold standard, coupled cluster with singles, doubles, and perturbative triples (CCSD(T)), provides the desired level of accuracy─roughly 1-3 kcal/mol. In this work, we compiled a test set of 28 3d metal-containing molecules relevant to homogeneous electrocatalysis (termed 3dTMV) and computed their vertical ionization energies (ionization potentials) with CCSD(T) and phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) in the def2-SVP basis set. A substantial effort has been made to converge away the phaseless bias in the ph-AFQMC reference values. We assess a wide variety of multireference diagnostics and find that spin-symmetry breaking of the CCSD wave function and the PBE0 density functional correlate well with our analysis of multiconfigurational wave functions. We propose quantitative criteria based on symmetry breaking to delineate correlation regimes inside of which appropriately performed CCSD(T) can produce mean absolute deviations from the ph-AFQMC reference values of roughly 2 kcal/mol or less and outside of which CCSD(T) is expected to fail. We also present a preliminary assessment of density functional theory (DFT) functionals on the 3dTMV set.

N-H Bond Activation Catalyzed by an Anderson-Type Polyoxometalate-Based Compound: Key Role of Transition-Metal Heteroatom

Polyoxometalates (POMs) have a broad array of applied platforms with well-characterized catalysis to achieve N-H bond activation. Herein, the mechanism of the Anderson-type POM-based catalyst [FeIIIMoVI6O18{(OCH2)3CNH2}2]3- ([TrisFeIIIMoVI6O18]3-, Tris = {(OCH2)3CNH2}2) for the N-H bond activation of hydrazine (PhHNNHPh) was investigated by density functional theory calculations. The results reveal that [TrisFeIIIMoVI6O18]3- as the active species is responsible for the continuous abstraction of two electrons and two protons of PhHNNHPh via a proton-coupled electron transfer pathway, resulting in the activation of two N-H bonds in PhHNNHPh and thus the product PhNNPh. H2O2 acts as an oxidant to regulate catalyst regeneration. Based on the proposed catalytic mechanism, the key role of the heteroatom FeIII in [TrisFeIIIMoVI6O18]3- was disclosed. The d-orbital of FeIII in [TrisFeIIIMoVI6O18]3- acts as an electron receptor to promote the electron transfer (ET) in the rate-determining step (RDS) of the catalytic cycle. The substitution of the heteroatom FeIII of [TrisFeIIIMoVI6O18]3- with CoIII, RuIII, or MnIII is expected to improve the catalytic activity for several reasons: (i) the unoccupied molecular orbitals of POM-based compounds containing CoIII or RuIII are low, which is beneficial for the ET of RDS; (ii) For N-H bond activation catalyzed by the MnIII-containing POM-based compound, the transition state of RDS is stable because the d-orbital of its active site is half-filled, which results in a low free-energy barrier.

Approaching Coupled Cluster Accuracy with Density Functional Theory Using the Generalized Connectivity-Based Hierarchy

This Perspective reviews connectivity-based hierarchy (CBH), a systematic hierarchy of error-cancellation schemes developed in our group with the goal of achieving chemical accuracy using inexpensive computational techniques ("coupled cluster accuracy with DFT"). The hierarchy is a generalization of Pople's isodesmic bond separation scheme that is based only on the structure and connectivity and is applicable to any organic and biomolecule consisting of covalent bonds. It is formulated as a series of rungs involving increasing levels of error cancellation on progressively larger fragments of the parent molecule. The method and our implementation are discussed briefly. Examples are given for the applications of CBH involving (1) energies of complex organic rearrangement reactions, (2) bond energies of biofuel molecules, (3) redox potentials in solution, (4) pK_a predictions in the aqueous medium, and (5) theoretical thermochemistry combining CBH with machine learning. They clearly show that near-chemical accuracy (1-2 kcal/mol) is achieved for a variety of applications with DFT methods irrespective of the underlying density functional used. They demonstrate conclusively that seemingly disparate results, often seen with different density functionals in many chemical applications, are due to an accumulation of systematic errors in the smaller local molecular fragments that can be easily corrected with higher-level calculations on those small units. This enables the method to achieve the accuracy of the high level of theory (e.g., coupled cluster) while the cost remains that of DFT. The advantages and limitations of the method are discussed along with areas of ongoing developments.

Ferrocene/ferrocenium, cobaltocene/cobaltocenium and nickelocene/nickelocenium: from gas phase ionization energy to one-electron reduction potential in solvated medium

We propose a theoretical procedure for accurate determination of reduction potentials for three metallocene couples, Cp₂M⁺/Cp₂M, where M = Fe, Co and Ni. This procedure first computes the gas phase ionization energy (IE) using the explicitly correlated CCSD(T)-F12 method and includes the zero-point energy correction, core-valence electronic correlation, and relativistic and spin-orbit coupling effects. By means of Born-Haber thermochemical cycle, the one-electron reduction potential is obtained as the sum of the gas phase IE and the corresponding Gibbs free energies of solvation (ΔG_solv) for both the neutral and cationic species. Among the three solvent models (PCM, SMD and uESE) investigated here, it turns out that only the SMD model (computed at the DFT level) gives the best estimation of the value for "ΔG_solv(cation) - ΔG_solv(neutral)" and thus, combining with the accurate IE values, the theoretical protocol is capable of yielding reliable values (in V) for , and . These predictions compare favorably with the available experimental data (in V): , , and . We show that our theoretical procedure is reliable for accurate reduction potential predictions of Cp₂Fe⁺/Cp₂Fe, Cp₂Co⁺/Cp₂Co and Cp₂Ni⁺/Cp₂Ni redox couples in aqueous and non-aqueous media; the maximum absolute deviation is as small as ≈120 mV, which outperforms those of the existing theoretical methods.

Mechanistic Insights into Photochemical CO2 Reduction to CH4 by a Molecular Iron-Porphyrin Catalyst

Iron tetraphenylporphyrin complex modified with four trimethylammonium groups (Fe-p-TMA) is found to be capable of catalyzing the eight-electron eight-proton reduction of CO₂ to CH₄ photochemically in acetonitrile. In the present work, density functional theory (DFT) calculations have been performed to investigate the reaction mechanism and to rationalize the product selectivity. Our results revealed that the initial catalyst Fe-p-TMA ([Cl-Fe(III)-LR₄]⁴⁺, where L = tetraphenylporphyrin ligand with a total charge of -2, and R₄ = four trimethylammonium groups with a total charge of +4) undergoes three reduction steps, accompanied by the dissociation of the chloride ion to form [Fe(II)-L^••2-R₄]²⁺. [Fe(II)-L^••2-R₄]²⁺, bearing a Fe(II) center ferromagnetically coupled with a tetraphenylporphyrin diradical, performs a nucleophilic attack on CO₂ to produce the ¹η-CO₂ adduct [CO₂^•--Fe(II)-L^•-R₄]²⁺. Two intermolecular proton transfer steps then take place at the CO₂ moiety of [CO₂^•--Fe(II)-L^•-R₄]²⁺, resulting in the cleavage of the C-O bond and the formation of the critical intermediate [Fe(II)-CO]⁴⁺ after releasing a water molecule. Subsequently, [Fe(II)-CO]⁴⁺ accepts three electrons and one proton to generate [CHO-Fe(II)-L^•-R₄]²⁺, which finally undergoes a successive four-electron-five-proton reduction to produce methane without forming formaldehyde, methanol, or formate. Notably, the redox non-innocent tetraphenylporphyrin ligand was found to play an important role in CO₂ reduction since it could accept and transfer electron(s) during catalysis, thus keeping the ferrous ion at a relatively high oxidation state. Hydrogen evolution reaction via the formation of Fe-hydride ([Fe(II)-H]³⁺) turns out to endure a higher total barrier than the CO₂ reduction reaction, therefore providing a reasonable explanation for the origin of the product selectivity.

Intramolecular and intermolecular hole delocalization rules the reducer character of isolated nucleobases and homogeneous single-stranded DNA

The use of DNA strands as nanowires or electrochemical biosensors requires a deep understanding of charge transfer processes along the strand, as well as of the redox properties. These properties are computationally assessed in detail throughout this study. By applying molecular dynamics and hybrid QM/continuum and QM/QM/continuum schemes, the vertical ionization energies, adiabatic ionization energies, vertical attachment energies, one-electron oxidation potentials, and delocalization of the hole generated upon oxidation have been determined for nucleobases in their free form and as part of a pure single-stranded DNA. We show that the reducer ability of the isolated nucleobases is explained by the intramolecular delocalization of the positively charged hole, while the enhancement of the reducer character when going from aqueous solution to the strand correlates very well with the intermolecular hole delocalization. Our simulations suggest that the redox properties of DNA strands can be tuned by playing with the balance between intramolecular and intermolecular charge delocalization.

Electrocatalytic and green system coupling strategy for simultaneous recovery and purification of uranium from uranium-containing wastewater

The recycling of uranium in wastewater is not only beneficial to the protection of ecological safety but also has great significance for the sustainable development of nuclear energy. However, there is no satisfactory method to recover and reuse uranium efficiently up to now. Here, we have developed an efficient and economical strategy that can achieve uranium recovery and direct reuse in wastewater. The feasibility analysis verified that the strategy still had good separation and recovery ability in acidic, alkaline, and high-salinity environments. The purity of uranium recovered from the separated liquid phase after electrochemical purification was up to about 99.95%. Ultrasonication could greatly increase the efficiency of this strategy, and 99.00% of high-purity uranium could be recovered within 2 h. We further improved the overall recovery rate by recovering the residual solid-phase uranium, and the overall recovery of uranium was increased to 99.40%. Moreover, the concentration of impurity ions in the recovered solution met the World Health Organization guidelines. In summary, the development of this strategy is of great importance for the sustainable use of uranium resources and environmental protection.

A Systematic Study on the Redox Potentials of Phenazine-Derivatives in Aqueous Media: A Combined Computational and Experimental Work

Phenazines are an emerging class of organic compounds that have been recently utilized in aqueous redox flow batteries, a promising technology for large-scale energy storage. A virtual screening based on density functional theory calculations is used to investigate the redox potentials of around 100 phenazine derivatives in aqueous media containing various electron-donating or electron-withdrawing groups at different positions. The calculations identify the crucial positions that should be functionalized with multiple hydroxy groups to design new anolytes. The combined experimental-computational methodology reported herein guides the development of a new molecule with a record low reversible redox potential as a potential anolyte for aqueous redox flow batteries.

QM Calculations Revealed that Outer-Sphere Electron Transfer Boosted O-O Bond Cleavage in the Multiheme-Dependent Cytochrome bd Oxygen Reductase

The cytochrome bd oxygen reductase catalyzes the four-electron reduction of dioxygen to two water molecules. The structure of this enzyme reveals three heme molecules in the active site, which differs from that of heme-copper cytochrome c oxidase. The quantum chemical cluster approach was used to uncover the reaction mechanism of this intriguing metalloenzyme. The calculations suggested that a proton-coupled electron transfer reduction occurs first to generate a ferrous heme b₅₉₅. This is followed by the dioxygen binding at the heme d center coupled with an outer-sphere electron transfer from the ferrous heme b₅₉₅ to the dioxygen moiety, affording a ferric ion superoxide intermediate. A second proton-coupled electron transfer produces a heme d ferric hydroperoxide, which undergoes efficient O-O bond cleavage facilitated by an outer-sphere electron transfer from the ferrous heme b₅₉₅ to the O-O σ* orbital and an inner-sphere proton transfer from the heme d hydroxyl group to the leaving hydroxide. The synergistic benefits of the two types of hemes rationalize the highly efficient oxygen reduction repertoire for the multi-heme-dependent cytochrome bd oxygen reductase family.

Theoretical study on the mechanism of water oxidation catalyzed by a mononuclear copper complex: important roles of a redox non-innocent ligand and HPO4 2- anion

The water oxidation reaction is the bottleneck problem of the artificial photosynthetic system. In this work, the mechanism of water oxidation catalyzed by a mononuclear copper complex in alkaline conditions was studied by density functional calculations. Firstly, a water molecule coordinating with the copper center of the complex (Cuii, 1) generates Cuii-H2O (2). 2 undergoes two proton-coupled electron transfer processes to produce intermediate (4). The oxidation process occurs mainly on the ligand moiety, and 4 (˙L-Cuii-O˙) can be described as a Cuii center interacting with a ligand radical antiferromagnetically and an oxyl radical ferromagnetically. 4 is the active species that can trigger O-O bond formation via the water nucleophilic attack mechanism. This process occurs in a step-wise manner. The attacking water transfers one of the protons to the HPO4 2- coupled with an electron transfer to the ligand radical, which generates a transient OH˙ interacting with the oxyl radical and H2PO4 -. Then the O-O bond is formed through the direct coupling of the oxo radical and the OH radical. The triplet di-oxygen could be released after two oxidation processes. According to the Gibbs free energy diagram, the O-O bond formation was suggested to be the rate-limiting step with a calculated total barrier of 19.5 kcal mol-1. More importantly, the copper complex catalyzing water oxidation with the help of a redox non-innocent ligand and HPO4 2- was emphasized.

An Efficient Multilayer Approach to Model DNA-Based Nanobiosensors

In this work, we present a full computational protocol to successfully obtain the one-electron reduction potential of nanobiosensors based on a self-assembled monolayer of DNA nucleobases linked to a gold substrate. The model is able to account for conformational sampling and environmental effects at a quantum mechanical (QM) level efficiently, by combining molecular mechanics (MM) molecular dynamics and multilayer QM/MM/continuum calculations within the framework of Marcus theory. The theoretical model shows that a guanine-based biosensor is more prone to be oxidized than the isolated nucleobase in water due to the electrostatic interactions between the assembled guanine molecules. In addition, the redox properties of the biosensor can be tuned by modifying the nature of the linker that anchor the nucleobases to the metal support.

Redox Potentials with COSMO-RS: Systematic Benchmarking with Different Databases

Recent techniques of computational electrochemistry can yield redox potentials with accuracy as good as 0.1 V. Yet, many such methods are not universal, easy to use, or computationally efficient. Herein, we provide a systematic benchmarking of a relatively cheap and straightforward computational approach for fairly accurate computations of redox potentials. It is based on a combination of the conductor-like screening model for real solvents (COSMO-RS) and the density functional theory (DFT). The benchmarking is done with databases covering diverse redox systems, including transition-metal aquacomplexes and various organic and inorganic compounds. In addition, we also present our own test set aiming at maximum chemical diversity and maximum range of redox potential values. We assess the performance of the fairly efficient computational protocol combining the COSMO-RS with the BP86 DFT functional. This is done by calibrating it against different high-level state-of-the-art techniques, in particular, polarizable continuum model (PCM) connected to composite G3(MP2,CC)(+) method, domain-based pair natural orbital implementation of coupled cluster theory, or complete basis set CBS-QB3 method. We also elaborate on the absolute reduction potential value of standard hydrogen electrode to be used with COSMO-RS, and we propose the value of approx. 4.4 V. The COSMO-RS/BP86-D3(BJ) combination outperforms other methods on a wide range of redox systems. However, we show that its accuracy is based on a balanced error cancelation and, therefore, it cannot be further systematically improved. As a result, the proposed procedure represents a pragmatic choice for large-scale screening, yet it could be combined with more advanced methods for detailed studies.

MoBioTools: A toolkit to setup quantum mechanics/molecular mechanics calculations

We present a toolkit that allows for the preparation of QM/MM input files from a conformational ensemble of molecular geometries. The package is currently compatible with trajectory and topology files in Amber, CHARMM, GROMACS and NAMD formats, and has the possibility to generate QM/MM input files for Gaussian (09 and 16), Orca (≥4.0), NWChem and (Open)Molcas. The toolkit can be used in command line, so that no programming experience is required, although it presents some features that can also be employed as a python application programming interface. We apply the toolkit in four situations in which different electronic-structure properties of organic molecules in the presence of a solvent or a complex biological environment are computed: the reduction potential of the nucleobases in acetonitrile, an energy decomposition analysis of tyrosine interacting with water, the absorption spectrum of an azobenzene derivative integrated into a voltage-gated ion channel, and the absorption and emission spectra of the luciferine/luciferase complex. These examples show that the toolkit can be employed in a manifold of situations for both the electronic ground state and electronically excited states. It also allows for the automatic correction of the active space in the case of CASSCF calculations on an ensemble of geometries, as it is shown for the azobenzene derivative photoswitch case.

Quantitative structure-activity relationship for the photooxidation of aromatic micro-pollutants induced by graphene oxide in water

The photocatalytic degradation behavior of aromatic micro-pollutants (AMPs) exhibits complexity and uncertainty, which mainly depends on the properties of different substituents on benzene. And with similar catalytic reaction substrates, the reaction rate constant could reveal the influence of different characteristics of molecular structure in a specific system. Therefore, the photooxidation pseudo first-order kinetic rate constants (k_obs) of 30 AMPs were experimentally determined in Photo-GO system. A quantitative structure-activity relationship (QSAR) model for predicting the photooxidation reaction of AMPs has been developed by stepwise multiple linear regression (MLR) method, based on the lg k_obs and representative molecule descriptors (20 in total) including physicochemical, quantum chemical and electrostatic descriptors. Afterwards, R_adj², Q_LOO², and Q_ext² were calculated as 0.870, 0.841, and 0.732 respectively, which exhibited the excellent goodness-of-fit, robustness, and predictability of the QSAR model, indicating its great prediction ability for photooxidation behavior of AMPs. Meanwhile, during the photooxidation process of AMPs with GO, the model revealed that the one-electron oxidation potential (E_ox), molecular dipole moment (μ), and number of hydrogen bond donors (#HD) were the most important molecular structural parameters, which showed that the single electron transfer pathway and adsorption were as the significant steps. Additionally, the Hammett correlation showed that photooxidation of AMPs in Photo-GO system is of typical electrophilic reactions, which demonstrated that the electron-donating substituents could promote the photooxidation of AMPs. The QSAR model was constructed and evaluated to perform the prediction of AMPs reaction kinetics, which provided a guidance for the study of the mechanism and selective oxidation of AOPs photooxidation system based on GO.

Mode of action of p-quinone derivatives with trypanocidal activity studied by experimental and in silico models

Quinones are attractive pharmacological scaffolds for developing new agents for the treatment of different transmissible and non-transmissible human diseases due to their capacity to alter the cell redox homeostasis. The bioactivity and potential mode of action of 19 p-quinone derivatives fused to different aromatic rings (carbo or heterocycles) and harboring distinct substituents were investigated in infective Trypanosoma brucei brucei. All the compounds, except for a furanequinone (EC₅₀₌38 μM), proved to be similarly or even more potent (EC₅₀ = 0.5-5.5 μM) than the clinical drug nifurtimox (EC₅₀ = 5.3 μM). Three furanequinones and one thiazolequinone displayed a higher selectivity than nifurtimox. Two of these selective hits resulted potent inhibitors of T. cruzi proliferation (EC₅₀₌0.8-1.1 μM) but proved inactive against Leishmania infantum amastigotes. Most of the p-quinones induced a rapid and marked intracellular oxidation in T. b. brucei. DFT calculations on the oxidized quinone (Q), semiquinone (Q^•-) and hydroquinone (QH₂) suggest that all quinones have negative ΔG for the formation of Q^•-. Qualitative and quantitative structure-activity relationship analyses in two or three dimensions of different electronic and biophysical descriptors of quinones and their corresponding bioactivities (killing potency and oxidative capacity) were performed. Charge distribution over the quinone ring carbons of Q and Q^.- and the frontier orbitals energies of SUMO (Q^.-) and LUMO (Q) correlate with their oxidative and trypanocidal activity. QSAR analysis also highlighted that both bromine substitution in the p-quinone ring and a bulky phenyl group attached to the furane and thiazole rings (which generates a negative charge due to the π electron system polarized by the nearby heteroatoms) are favorable for activity. By combining experimental and in silico procedures, this study disclosed important information about p-quinones that may help to rationally tune their electronic properties and biological activities.

Protocols for Understanding the Redox Behavior of Copper-Containing Systems

Suitability of single-reference density functional theory (DFT) methods for the calculation of redox potentials of copper-containing macrocycle complexes was confirmed by the use of T ₁ diagnostics along with a verification of negligible spin contamination or wave function instability. When examining the effect of improvement in the cc-pVnZ basis set series on calculated redox potentials, the results readily converged at the cc-pVTZ level. The all-electron Def2-TZVPP basis set is shown to be a suitable choice of a basis set for the calculation of redox potentials when utilizing a cc-pVTZ geometry. The best-performing model chemistries are determined to be the M06/polarizable continuum model (PCM); therefore, a scheme for redox potential calculations of copper macrocycles using either M06/cc-pVTZ with PCM solvation is proposed to reliably reproduce experimental trends.

Explicit and Hybrid Solvent Models for Estimates of Parameters Relevant to the Reduction Potential of Ethylene Carbonate

Using ethylene carbonate as a sample solvent, we investigated two molecular parameters used to estimate the reduction potential of the solvent: electron affinity, and the energy of the lowest unoccupied molecular orbital (LUMO). The results showed that the values of these parameters are inconsistent for a single ethylene carbonate molecule in vacuum calculations and in the continuous effective solvent. We performed a series of calculations employing explicit or hybrid (explicit/continuous) solvent models for aggregates of solvent molecules or solvated salt ions. In the hybrid solvent model, values of the two estimates extrapolated to an infinite system size converged to one common value, whereas the difference of 1 eV was calculated in the purely explicit solvent. The values of the gap between the highest occupied molecular orbital (HOMO) and the LUMO obtained in the hybrid model were significantly larger than those resulting from the explicit solvent calculations. We related these differences to the differences in frontier orbitals and changes of electron density obtained in the two solvent models. In the hybrid solvent model, the location of the additional electron in the reduced system usually corresponds to the LUMO orbital of the oxidized system. The presence of salt ions in the solvent affects the extrapolated values of the electron affinity and LUMO energy.

The role of tridentate ligands on the redox stability of anticancer gold(III) complexes

Gold(III) complexes are promising compounds for cancer chemotherapy, whose action depends on their redox stability. In this context, the choice of ligands is crucial to adjust their reactivity and biological response. The present study addressed the effect of the gold coordination sphere on the reduction potential (E^o) for ten gold(III) complexes containing five or six-membered rings tridentate ligands - [Au^III(trident)Cl]³⁺ⁿ (trident = N^N^N, C^N^N, C^C^N, C^N^C, and N^C^N). The calculated E^o covered a broad range of 2500 mV with the most stable complexes containing two AuC bonds (E^o = -1.85 V for [Au^III(C^C^N)Cl] - f). For complexes with one AuC bond, the N^C^N ligands stabilize the gold(III) complex more efficiently than N^N^C; however, the inclusion of the non-innocent ligand bipy (2,2'-bipyridine) in N^N portion provides an extra stabilization effect. Among the derivatives with one AuC bond, [Au^III(N^N^C)Cl]⁺ (N^N = bipy) (a) showed E^o = -1.20 V. For the complexes with N^N^N ligands, E^o was positive and almost constant (+0.60 V). Furthermore, the kinetics for ligand exchange reactions (Cl^-/H₂O, H₂O/Cys and Cl^-/Cys) were monitored for the most stable compounds and the energy profiles compared to the reduction pathways.

Hydrogen Atom Abstraction from an OsII(NH3)2 Complex Generates an OsIV(NH2)2 Complex: Experimental and Computational Analysis of the N-H Bond Dissociation Free Energies and Reactivity

Double hydrogen atom abstraction from (TMP)Os^II(NH₃)₂ (TMP = tetramesitylporphyrin) with phenoxyl or nitroxyl radicals leads to (TMP)Os^IV(NH₂)₂. This unusual bis(amide) complex is diamagnetic and displays an N-H resonance at 12.0 ppm in its ¹H NMR spectrum. ¹H-¹⁵N correlation experiments identified a ¹⁵N NMR spectroscopic resonance signal at -267 ppm. Experimental reactivity studies and density functional theory calculations support relatively weak N-H bonds of 73.3 kcal/mol for (TMP)Os^II(NH₃)₂ and 74.2 kcal/mol for (TMP)Os^III(NH₃)(NH₂). Cyclic voltammetry experiments provide an estimate of the pK_a of [(TMP)Os^III(NH₃)₂]⁺. In the presence of Barton's base, a current enhancement is observed at the Os(III/II) couple, consistent with an ECE event. Spectroscopic experiments confirmed (TMP)Os^IV(NH₂)₂ as the product of bulk electrolysis. Double hydrogen atom abstraction is influenced by π donation from the amides of (TMP)Os^IV(NH₂)₂ into the d orbitals of the Os center, favoring the formation of (TMP)Os^IV(NH₂)₂ over N-N coupling. This π donation leads to a Jahn-Teller distortion that splits the energy levels of the d_xz and d_yz orbitals of Os, results in a low-spin electron configuration, and leads to minimal aminyl character on the N atoms, rendering (TMP)Os^IV(NH₂)₂ unreactive toward amide-amide coupling.

Trade-Off between Redox Potential and the Strength of Electrochemical CO2 Capture in Quinones

Electrochemical carbon dioxide capture recently emerged as a promising alternative approach to conventional energy-intensive carbon-capture methods. A common electrochemical capture approach is to employ redox-active molecules such as quinones. Upon electrochemical reduction, quinones become activated for the capture of CO₂ through a chemical reaction. A key disadvantage of this method is the possibility of side-reactions with oxygen, which is present in almost all gas mixtures of interest for carbon capture. This issue can potentially be mitigated by fine-tuning redox potentials through the introduction of electron-withdrawing groups on the quinone ring. In this article, we investigate the thermodynamics of the electron transfer and chemical steps of CO₂ capture in different quinone derivatives with a range of substituents. By combining density functional theory calculations and cyclic voltammetry experiments, we support a previously described trade-off between the redox potential and the strength of CO₂ capture. We show that redox potentials can readily be tuned to more positive values to impart stability to oxygen, but significant decreases in CO₂ binding free energies are observed as a consequence. Our calculations support this effect for a large series of anthraquinones and benzoquinones. Different trade-off relationships were observed for the two classes of molecules. These trade-offs must be taken into consideration in the design of improved redox-active molecules for electrochemical CO₂ capture.

Water Oxidation Catalyzed by a Bioinspired Tetranuclear Manganese Complex: Mechanistic Study and Prediction

Density functional theory calculations were utilized to elucidate the water oxidation mechanism catalyzed by polyanionic tetramanganese complex a [Mn^III ₃ Mn^IV O₃ (CH₃ COO)₃ (A-α-SiW₉ O₃₄ )]^6- . Theoretical results indicated that catalytic active species 1 (Mn₄ ^{III,III,IV,IV} ) was formed after O₂ formation in the first turnover. From 1, three sequential proton-coupled electron transfer (PCET) oxidations led to the Mn^IV -oxyl radical 4 (Mn₄ ^IV,IV,IV,IV -O⋅). Importantly, 4 had an unusual butterfly-shaped Mn₂ O₂ core for the two substrate-coordinated Mn sites, which facilitated O-O bond formation via direct coupling of the oxyl radical and the adjacent Mn^IV -coordinated hydroxide to produce the hydroperoxide intermediate Int1 (Mn₄ ^III,IV,IV,IV -OOH). This step had an overall energy barrier of 24.9 kcal mol^-1 . Subsequent PCET oxidation of Int1 to Int2 (Mn₄ ^III,IV,IV,IV -O₂ ⋅) enabled the O₂ release in a facile process. Furthermore, apart from the Si-centered complex, computational study suggested that tetramanganese polyoxometalates with Ge, P, and S could also catalyze the water oxidation process, where those bearing P and S likely present higher activities.

Feasible Cluster Model Method for Simulating the Redox Potentials of Laccase CueO and Its Variant

Laccases are regarded as versatile green biocatalysts, and recent scientific research has focused on improving their redox potential for broader industrial and environmental applications. The density functional theory (DFT) quantum mechanics approach, sufficiently rigorous and efficient for the calculation of electronic structures, is conducted to better comprehend the connection between the redox potential and the atomic structural feature of laccases. According to the crystal structure of wild type laccase CueO and its variant, a truncated miniature cluster model method was established in this research. On the basic of thermodynamic cycle, the overall Gibbs free energy variations before and after the one-electron reduction were calculated. It turned out that the trends of redox potentials to increase after variant predicted by the theoretical calculations correlated well with those obtained by experiments, thereby validating the feasibility of this cluster model method for simulating the redox potentials of laccases.

The Role of Redox-Inactive Metals in Modulating the Redox Potential of the Mn4CaO4 Model Complex

Photosynthetic oxygen-evolving center (OEC), the "engine of life", is a unique Mn₄CaO₅ cluster catalyzing the water oxidation. The role of redox-inactive component Ca²⁺, which can only be functionally replaced by Sr²⁺ in a biological environment, has been under debate for a long time. Recently, its modulating effect on the redox potential of native OEC and artificial structural OEC model complex has received great attention, and linear relationship between the potential and the Lewis acidity of the redox-inactive metal has been proposed for the MMn₃O₄ model complex. In this work, the modulating effect has been studied in detail using the Mn₄CaO₄ model complex, which is the closest structural model to OEC to date and has a similar redox potential at the S₁-S₂ transition. We found the redox-inactive metal only has a weak modulating effect on the potential, which is comparable in strength to that of the ligand environments. Meanwhile, the net charge of the complex, which could be changed along with the redox-inactive metal, has a high impact on the potential and can be unified by protonation, deprotonation, or ligand modification. Although the modulating effect of the redox-inactive metal is not very strong, the linear relationship between the potential and the Lewis acidity is still valid for Mn₄MO₄ complexes. Our results of strong modulating effects for net charge and weak modulating effects for redox-inactive metal fit with the previous experimental observations on Mn₄MO₄ (M = Ca²⁺, Y³⁺, and Gd³⁺) model complexes, and suggest that Ca²⁺ can be structurally and electrochemically replaced with other metal cations, together with proper ligand modifications.

Computation of Oxidation Potentials of Solvated Nucleobases by Static and Dynamic Multilayer Approaches

The determination of the redox properties of nucleobases is of paramount importance to get insight into the charge-transfer processes in which they are involved, such as those occurring in DNA-inspired biosensors. Although many theoretical and experimental studies have been conducted, the value of the one-electron oxidation potentials of nucleobases is not well-defined. Moreover, the most appropriate theoretical protocol to model the redox properties has not been established yet. In this work, we have implemented and evaluated different static and dynamic approaches to compute the one-electron oxidation potentials of solvated nucleobases. In the static framework, two thermodynamic cycles have been tested to assess their accuracy against the direct determination of oxidation potentials from the adiabatic ionization energies. Then, the introduction of vibrational sampling, the effect of implicit and explicit solvation models, and the application of the Marcus theory have been analyzed through dynamic methods. The results revealed that the static direct determination provides more accurate results than thermodynamic cycles. Moreover, the effect of sampling has not shown to be relevant, and the results are improved within the dynamic framework when the Marcus theory is applied, especially in explicit solvent, with respect to the direct approach. Finally, the presence of different tautomers in water does not affect significantly the one-electron oxidation potentials.

Ligand-Centered Hydrogen Evolution with Ni(II) and Pd(II)DMTH

In this study, we report a pair of electrocatalysts for the hydrogen evolution reaction (HER) based on the noninnocent ligand diacetyl-2-(4-methyl-3-thiosemicarbazone)-3-(2-pyridinehydrazone) (H₂DMTH, H₂L¹). The neutral complexes NiL¹ and PdL¹ were synthesized and characterized by spectroscopic and electrochemical methods. The complexes contain a non-coordinating, basic hydrazino nitrogen that is protonated during the HER. The pK_a of this nitrogen was determined by spectrophotometric titration in acetonitrile to be 12.71 for NiL¹ and 13.03 for PdL¹. Cyclic voltammograms of both NiL¹ and PdL¹ in acetonitrile exhibit diffusion-controlled, reversible ligand-centered events at -1.83 and -1.79 V (vs ferrocenium/ferrocene) for NiL¹ and PdL¹, respectively. A quasi-reversible, ligand-centered event is observed at -2.43 and -2.34 V for NiL¹ and PdL¹, respectively. The HER activity in acetonitrile was evaluated using a series of neutral and cationic acids for each catalyst. Kinetic isotope effect (KIE) studies suggest that the precatalytic event observed is associated with a proton-coupled electron transfer step. The highest turnover frequency values observed were 6150 s^-1 at an overpotential of 0.74 V for NiL¹ and 8280 s^-1 at an overpotential of 0.44 V for PdL¹. Density functional theory (DFT) computations suggest both complexes follow a ligand-centered HER mechanism where the metals remain in the +2 oxidation state.