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

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

Mechanochemical synthesis of iron aluminyl complexes

A series of iron aluminyl complexes have been synthesised in good crystalline yields from reactions between bulky diamido aluminium iodide complexes and K[Fe(CO)2Cp] in the solid state. The series of metal-metal bonded complexes have been characterised by X-ray crystallography and were investigated using density functional theory to probe the effects of ligand substitution on the Al-Fe bond.

Exploring the interaction between a fluorescent Ag(I)-biscarbene complex and non-canonical DNA structures: a multi-technique investigation

Silver compounds are mainly studied as antimicrobial agents, but they also have anticancer properties, with the latter, in some cases, being better than their gold counterparts. Herein, we analyse the first example of a new Ag(I)-biscarbene that can bind non-canonical structures of DNA, more precisely G-quadruplexes (G4), with different binding signatures depending on the type of G4. Moreover, we show that this Ag-based carbene binds the i-motif DNA structure. Alternatively, its Au(I) counterpart, which was investigated for comparison, stabilises mitochondrial G4. Theoretical in silico studies elucidated the details of different binding modes depending on the geometry of G4. The two complexes showed increased cytotoxic activity compared to cisplatin, overcoming its resistance in ovarian cancer. The binding of these new drug candidates with other relevant biosubstrates was studied to afford a more complete picture of their possible targets. In particular, the Ag(I) complex preferentially binds DNA structures over RNA structures, with higher binding constants for the non-canonical nucleic acids with respect to natural calf thymus DNA. Regarding possible protein targets, its interaction with the albumin model protein BSA was also tested.

Evaluation of DFT methods for predicting geometries and NMR spectra of Bi(III) dithiocarbamate complexes with antitumor properties

CONTEXT

Bismuth complexes with dithiocarbamate ligands have attracted attention because of their biological applications, such as antimicrobial, antileishmanial, and anticancer properties. These complexes have high cytotoxic activity against cancer cells, being more active than the standard drugs cisplatin, doxorubicin, and tamoxifen. In the present study, we investigated the ability of some DFT methods to reproduce the geometries and NMR spectra of the Bi(III) dithiocarbamate complexes, selected based on their proven antitumor activity. Our investigation revealed that the M06-L/def2-TZVP/ECP/CPCM method presented good accuracy in predicting geometries, while the TPSSh/def2-SVP/ECP/CPCM method proved effective in analyzing the 13C NMR spectra of these molecules. In general, all examined methods exhibited comparable performance in predicting 1H NMR signals.

METHODS

Calculations were performed with the Gaussian 09 program using the def2-SVP and def2-TZVP basis sets, employing relativistic effective core potential (ECP) for Bi and using the CPCM solvent model. The exchange-correlation functionals BP86, PBE, OLYP, M06-L, B3LYP, B3LYP-D3, M06-2X, TPSSh, CAM-B3LYP, and ωB97XD were used in the study. Geometry optimizations were started from crystallographic structures available at the Cambridge Structural Database. The theoretical results were compared with experimental data using the mean root-mean-square deviation (RMSD), mean absolute deviations (MAD), and linear correlation coefficient (R2).

Deciphering the Origin of Regioselectivity in Ru(II)-Catalyzed C-H Annulation of N-Chlorobenzamides with 1,3-Diynes

Understanding the reaction mechanism and origin of regioselectivity in transition metal-catalyzed C-H activation/annulation reactions with 1,3-diynes has remained an intriguing challenge. In this article, to establish the mechanism and decipher the origin of regioselectivity, we report a detailed computational density functional theory-based mechanistic investigation on the recently developed Ru(II)-catalyzed [4 + 2] annulation of N-chlorobenzamides with 1,3-diynes for the synthesis of 3-alkynylated isoquinolone derivatives. Our calculations reveal a redox-neutral pathway for the annulation reaction. The stepwise analysis of the reaction channels indicates the migratory insertion step and the concerted reductive elimination/oxidative addition of the Ru(p-cymene) moiety to form the N-C bond leading to the 3-alkynylated product to be the selectivity- and rate-determining steps, respectively. Finally, the distortion/interaction analysis using the activation-strain model suggests the steric effect as the determining factor for the observed regioselectivity for the formation of the 3-alkynylated product. Overall, the computationally obtained key insights into the catalytic mechanism and the origin of regioselectivity in the C-H activation/annulation reaction can be used as a guide to rationally design and develop novel transformation strategies for heterocycle synthesis.

Electron beam-induced demetallation of Fe, Co, Ni, Cu, Zn, Pd, and Pt metalloporphyrins: insights in e-beam chemistry and metal cluster formations

Electron beams are versatile tools for nanoscale fabrication processes, however, the underlying e-beam chemistry remains in its infancy. Through operando transmission electron microscopy investigations, we elucidate a redox-driven cargo release of individual metal atoms triggered by electron beams. The chosen organic delivery molecule, tetraphenylporphyrin (TPP), proves highly versatile, forming complexes with nearly all metals from the periodic table and being easily processed in solution. A comprehensive cinematographic analysis of the dynamics of single metal atoms confirms the nearly instantaneous ejection of complexed metal atoms under an 80 kV electron beam, underscoring the system's broad versatility. Providing mechanistic insights, we employ density functional theory to support the proposed reductive demetallation pathway facilitated by secondary electrons, contributing novel perspectives to electron beam-mediated chemical reaction mechanisms. Lastly, our findings demonstrate that all seven metals investigated form nanoclusters once ejected from TPP, highlighting the method's potential for studying and developing sustainable single-atom and nanocluster catalysts.

Mechanism of Action of Antitumor Au(I) N-Heterocyclic Carbene Complexes: A Computational Insight on the Targeting of TrxR Selenocysteine

The targeting of human thioredoxin reductase is widely recognized to be crucially involved in the anticancer properties of several metallodrugs, including Au(I) complexes. In this study, the mechanism of reaction between a set of five N-heterocyclic carbene Au(I) complexes and models of the active Sec residue in human thioredoxin reductase was investigated by means of density functional theory approaches. The study was specifically addressed to the kinetics and thermodynamics of the tiled process by aiming at elucidating and explaining the differential inhibitory potency in this set of analogous Au(I) bis-carbene complexes. While the calculated free energy profile showed a substantially similar reactivity, we found that the binding of these Au(I) bis-carbene at the active CysSec dyad in the TrxR enzyme could be subjected to steric and orientational restraints, underlining both the approach of the bis-carbene scaffold and the attack of the selenol group at the metal center. A new and detailed mechanistic insight to the anticancer activity of these Au(I) organometallic complexes was thus provided by consolidating the TrxR targeting paradigm.

Combined Machine Learning, Computational, and Experimental Analysis of the Iridium(III) Complexes with Red to Near-Infrared Emission

Various coordination complexes have been the subject of experimental and theoretical studies in recent decades because of their fascinating photophysical properties. In this work, a combined experimental and computational approach was applied to investigate the optical properties of monocationic Ir(III) complexes. An interpretative machine learning-based quantitative structure-property relationship (ML/QSPR) model was successfully developed that could reliably predict the emission wavelength of the Ir(III) complexes and provide a foundation for the theoretical evaluation of the optical properties of Ir(III) complexes. A hypothesis was proposed to explain the differences in the emission wavelengths between structurally different individual Ir(III) complexes. The efficacy of the developed model was demonstrated by high R2 values of 0.84 and 0.87 for the training and test sets, respectively. It is worth noting that a relationship between the N-N distance in the diimine ligands of the Ir(III) complexes and emission wavelengths is detected. This effect is most probably associated with a degree of distortion in the octahedral geometry of the complexes, resulting in a perturbed ligand field. This combined experimental and computational approach shows great potential for the rational design of new Ir(III) complexes with the desired optical properties. Moreover, the developed methodology could be extended to other transition-metal complexes.

One Collision-Two Substituents: Gas-Phase Preparation of Xylenes under Single-Collision Conditions

The fundamental reaction pathways to the simplest dialkylsubstituted aromatics-xylenes (C6 H4 (CH3 )2 )-in high-temperature combustion flames and in low-temperature extraterrestrial environments are still unknown, but critical to understand the chemistry and molecular mass growth processes in these extreme environments. Exploiting crossed molecular beam experiments augmented by state-of-the-art electronic structure and statistical calculations, this study uncovers a previously elusive, facile gas-phase synthesis of xylenes through an isomer-selective reaction of 1-propynyl (methylethynyl, CH3 CC) with 2-methyl-1,3-butadiene (isoprene, C5 H8 ). The reaction dynamics are driven by a barrierless addition of the radical to the diene moiety of 2-methyl-1,3-butadiene followed by extensive isomerization (hydrogen shifts, cyclization) prior to unimolecular decomposition accompanied by aromatization via atomic hydrogen loss. This overall exoergic reaction affords a preparation of xylenes not only in high-temperature environments such as in combustion flames and around circumstellar envelopes of carbon-rich Asymptotic Giant Branch (AGB) stars, but also in low-temperature cold molecular clouds (10 K) and in hydrocarbon-rich atmospheres of planets and their moons such as Triton and Titan. Our study established a hitherto unknown gas-phase route to xylenes and potentially more complex, disubstituted benzenes via a single collision event highlighting the significance of an alkyl-substituted ethynyl-mediated preparation of aromatic molecules in our Universe.

Double Centrosymmetric Si···π Tetrel Bonds as New Synthons─Evidence from Crystal Structures and DFT Calculations

The crystal structure of bis((μ2-ethynylsilyloxo)-dichloro-aluminum), BEDCA, and a few related structures are characterized by the occurrence of tetrel bonds that link molecules. Particularly, centosymmetric dimers in such structures occur that are connected by two equivalent Si···π tetrel bonds. The dimer of BEDCA and dimers of other model species that similarly are linked by two equivalent Si···π tetrel bonds are analyzed theoretically. Some of the complexes calculated here are also characterized by the occurrence of triel bonds. Thus, ωB97XD/aug-cc-pVTZ calculations are performed and these DFT results are further supported by calculations with the use of other theoretical approaches: the quantum theory of atoms in molecules, QTAIM; the natural bond orbital, NBO; the energy decomposition analysis, EDA; and the noncovalent interactions method, NCI. The results show that the tetrel bonds analyzed here are rather weak, and they are not detected by the QTAIM approach; however, they are detected by other approaches, like NBO, for example. On the other hand, the triel bonds that occur in a few complexes discussed here are very strong and possess characteristics of covalent bonds.

Chiral-at-Metal Racemization Unraveled for MX2 (a-chel)2 by means of a Computational Analysis of MoO2 (acnac)2

Octahedral chiral-at-metal complexes MX2 (a-chel)2 (a-chel=asymmetric chelate) can rearrange their ligands by four mechanisms known as the Bailar (B), Ray-Dutt (RD), Conte-Hippler (CH), and Dhimba-Muller-Lammertsma (DML) twists. Racemization involves their interconnections, which were computed for MoO2 (acnac)2 (acnac=β-ketoiminate) using density functional theory at ωB97x-D with the 6-31G(d,p) and 6-311G(2d,p) basis sets and LANL2DZ for molybdenum. Racemizing the cis(NN) isomer, being the global energy minimum with trans oriented imine groups, is a three step process (DML-CH-DML) that requires 17.4 kcal/mol, while all three cis isomers (cis(NN), cis(NO), and cis(OO)) interconvert at ≤17.9 kcal/mol. The B and RD twists are energetically not competitive and neither are the trans isomers. The interconnection of all enantiomeric minima and transition structures is summarized in a graph that also visualizes valley ridge inflection points for two of the three CH twists. Geometrical features of the minima and twists are given. Lastly, the influence of N-substitution on the favored racemization pathway is evaluated. The present comprehensive study serves as a template for designing chiral-at-metal MX2 (a-chel)2 catalysts that may retain their chiral integrity.

An Anthracene-Thiolate-Ligated Ruthenium Complex: Computational Insights into Z-Stereoselective Cross Metathesis

Stereoselective control of the cross metathesis of olefins is a crucial aspect of synthetic procedures. In this study, we utilized density functional theory methods to calculate thermodynamic and kinetic descriptors to explore the stereoselectivity of cross metathesis between allylbenzene and 2-butene-1,4-diyl diacetate. A ruthenium-based complex, characterized primarily by an anthracene-9-thiolate ligand, was designed in silico to completely restrict the E conformation of olefins through a bottom-bound mechanism. Our investigation of the kinetics of all feasible propagation routes demonstrated that Z-stereoisomers of metathesis products can be synthesized with an energy cost of only 13 kcal/mol. As a result, we encourage further research into the synthetic strategies outlined in this work.

Recent Advancements in Mechanistic Studies of Palladium- and Nickel-Catalyzed Ethylene Copolymerization with Polar Monomers

The introduction of polar functional groups into polyolefin chain structures creates opportunities to enhance specific properties, such as adhesion, dyeability, printability, compatibility, thermal stability, and electrical conductivity, which widen the range of potential applications for these modified materials. Transition metal catalysts, especially late transition metals, have proven to be highly effective in copolymerization processes due to their reduced Lewis acidity and electrophilicity. However, when compared to the significant progress and summary of synthetic methods, there is a distinct lack of a comprehensive summary of mechanistic studies pertaining to the catalytic systems involved in ethylene copolymerization catalyzed by palladium and nickel catalysts. In this review, we have provided a comprehensive summary of the latest developments in mechanistic studies of ethylene copolymerization with polar monomers catalyzed by late-transition-metal complexes. Experimental and computational methods were employed to conduct a detailed investigation of these organic and organometallic systems. It is mainly focused on ligand substitution, changes in binding modes, ethylene/polar monomer insertion, chelate opening, and β-H elimination. Factors that control the catalytic activity, molecular weight, comonomer incorporation ratios, and branch content are analyzed, these include steric repulsions between ligands and monomers, electronic effects arising from both ligands and monomers, and so on.

DFT Study of Imine-Exchange Reactions in Iron(II)-Coordinated Pincers

The imine bond is among the most applied motifs in dynamic covalent chemistry. Although its uses are varied and often involve coordination to a transition metal for stability, mechanistic studies on imine exchange reactions so far have not included metal coordination. Herein, we investigated the condensation and transimination reactions of an Fe2+ -coordinated diimine pyridine pincer, employing wB97XD/6-311G(2d,2p) DFT calculations in acetonitrile. We first experimentally confirmed that Fe2+ is strongly coordinated by these pincers, and is thus a justified model ion. When considering a four-membered ring-shaped transition state for proton transfers, the required activation energies for condensation and transimination reaction exceeded the values expected for reactions known to be spontaneous at room temperature. The nature of the incoming and exiting amines and the substituents on the para-position of the pincer had no effect on this. Replacing Fe2+ with Zn2+ or removing it altogether did not reduce it either. However, the addition of two ethylamine molecules lowered the energy barriers to be compatible with experiment (19.4 and 23.2 kcal/mol for condensation and transimination, respectively). Lastly, the energy barrier of condensation of a non-coordinated pincer was significantly higher than found for Fe2+ -coordinating pincers, underlining the catalyzing effect of metal coordination on imine exchange reactions.

Molecular Structure of Gaseous Oxopivalate Co(II): Electronic States of Various Multiplicities

Synchronous electron diffraction/mass spectrometry was used to study the composition and structure of molecular forms existing in a saturated vapor of cobalt(II) oxopivalate at T = 410 K. It was found that monomeric complexes Co4O(piv)6 dominate in the vapor. The complex geometry possesses the C3 symmetry with bond lengths Co-Oc = 1.975(5) Å and Co-O = 1.963(5) Å, as well as bond angles Oc-Co-O = 111.8(3)°, Co-Oc-Co = 110.4(6)°, O-Co-O = 107.1(3)° in the central OcCo4 fragment and four OcCoO3 fragments. The presence of an open 3d shell for each Co atom leads to the possibility of the existence of electronic states of the Co4O(piv)6 complex with Multiplicities 1, 3, 5, 7, 9, 11, and 13. For them, the CASSCF and XMCQDPT2 calculations predict similar energies, identical shapes of active orbitals, and geometric parameters, the difference between which is comparable with the error of determination by the electron diffraction experiment. QTAIM and NBO analysis show that the Co-Oc and Co-O bonds can be attributed to ionic (or coordination) bonds with a significant contribution of the covalent component. The high volatility and simple vapor composition make it possible to recommend cobalt (II) oxopivalate as precursors in the preparation of oxide films or coatings in the CVD technologies. The features of the electronic and geometric structure of the Co4O(piv)6 complex allows for the conclude that only a very small change in energy is required for the transition from antiferromagnetically to ferromagnetically coupled Co atoms.

Assessment of Density Functional Theory Methods for the Structural Prediction of Transition and Post-Transition Metal-Nucleic Acid Complexes

Understanding the structure of metal-nucleic acid systems is important for many applications such as the design of new pharmaceuticals, metal detection platforms, and nanomaterials. Herein, we explore the ability of 20 density functional theory (DFT) functionals to reproduce the crystal structure geometry of transition and post-transition metal-nucleic acid complexes identified in the Protein Data Bank and Cambridge Structural Database. The environmental extremes of the gas phase and implicit water were considered, and analysis focused on the global and inner coordination geometry, including the coordination distances. Although gas-phase calculations were unable to describe the structure of 12 out of the 53 complexes in our test set regardless of the DFT functional considered, accounting for the broader environment through implicit solvation or constraining the model truncation points to crystallographic coordinates generally afforded agreement with the experimental structure, suggesting that functional performance for these systems is likely due to the models rather than the methods. For the remaining 41 complexes, our results show that the reliability of functionals depends on the metal identity, with the magnitude of error varying across the periodic table. Furthermore, minimal changes in the geometries of these metal-nucleic acid complexes occur upon use of the Stuttgart-Dresden effective core potential and/or inclusion of an implicit water environment. The overall top three performing functionals are ωB97X-V, ωB97X-D3(BJ), and MN15, which reliably describe the structure of a broad range of metal-nucleic acid systems. Other suitable functionals include MN15-L, which is a cheaper alternative to MN15, and PBEh-3c, which is commonly used in QM/MM calculations of biomolecules. In fact, these five methods were the only functionals tested to reproduce the coordination sphere of Cu²⁺-containing complexes. For metal-nucleic acid systems that do not contain Cu²⁺, ωB97X and ωB97X-D are also suitable choices. These top-performing methods can be utilized in future investigations of diverse metal-nucleic acid complexes of relevance to biology and material science.

Bracelet-like Complexes of Lithium Fluoride with Aromatic Tetraamides, and Their Potential for LiF-Mediated Self-Assembly: A DFT Study

Geometries and binding energies of complexes between a LiF molecule and a model aromatic tetraamide are obtained using various DFT methods. The tetraamide consists of a benzene ring and four amides positioned so that the LiF molecule can bind via Li⋯O=C or N-H⋯F interactions. The complex with both interactions is the most stable one, followed by the complex with only N-H⋯F interactions. Doubling the size of the former resulted in a complex with a LiF dimer sandwiched between the model tetraamides. In turn, doubling the size of the latter resulted in a more stable tetramer with bracelet-like geometry having the two LiF molecules also sandwiched but far apart from each other. Additionally, all methods show that the energy barrier to transition to the more stable tetramer is small. The self-assembly of the bracelet-like complex mediated by the interactions of adjacent LiF molecules is demonstrated by all computational methods employed.

Mononuclear half-sandwich nd7 metallo drug complexes based on bidentate N∩N dendritic scaffolds: DFT (B3LYP; BP86 and B3PW91) examination

Through an use of three functionals (B3PW91, B3LYP and BP86) associated to a generic basis set LanL2DZ for transition metals (as well as halogen atoms) and 6-311+G (d,p) for others atoms, an examination of the bonding properties of a series of mononuclear half-sandwich nd⁷ transition metal (anticancer) complexes based on N∩N dendritic scaffolds (L) has been done. Collectively, complexes studied have adopted the piano-stool environment. An examination of the performance of each functional has shown that for the most reliable geometrical analysis of Metal-Nitrogen and Metal-Halogen bonds, the B3LYP and B3PW91 functionalities are suitable respectively. Regardless of the halogen ligand adopted, the B3LYP metal-nitrogen bond lengths are the most widely overestimated. A correlation has been built between the retained charge on each divalent transition metal cation and its metal ion affinity (MIA). Topological examinations reveal the higher instability of metal-N bonds compared to metal-X ones (X = Cl and Br). By the mean of the energy decomposition analysis, a predominant electrostatic character of metal … halogen and [LCP]^- … [MX]⁺ interaction has been demonstrated. The transition metal atom … (hydrophobic) surface (Cp*) interaction is most pronounced for the chloride rhodium complexes of rhodium (combined with (E)-N-(pyridin-2-ylmethylene) Propan-1-amine and 2,2'- dipyridylketone ligands and iridium combined with 2,2'- dipyridylketone ligand. The charge decomposition analysis displays the weakening of the [Formula: see text] bonds in the studied complexes.

Site-Averaged Ab Initio Kinetics: Importance Learning for Multistep Reactions on Amorphous Supports

Single-atom centers on amorphous supports include catalysts for polymerization, partial oxidation, metathesis, hydrogenolysis, and more. The disordered environment makes each site different, and the kinetics exponentially magnifies these differences to make ab initio site-averaged kinetics calculations extremely difficult. This work extends the importance learning algorithm for efficient and precise site-averaged kinetics estimates to ab initio calculations and multistep reaction mechanisms. Specifically, we calculate site-averaged proton transfer relaxation rates on an ensemble of cluster models representing Brønsted acid sites on silica-alumina. We include direct and water-assisted proton transfer pathways and simultaneously estimate the water adsorption and activation enthalpies for forward and backward proton transfers. We use density functional theory (DFT) to obtain a site-averaged rate, somewhat like a turnover frequency, for the proton transfer relaxation rate. Finally, we show that importance learning can provide orders-of-magnitude acceleration over standard sampling methods for site-averaged rate calculations in cases where the rate is dominated by a few highly active sites.

Improving a Methane C-H Activation Complex by Metal and Ligand Alterations from Computational Results

We present results for a series of complexes derived from a titanium complex capable of activating C-H bonds under mild conditions (PNP)Ti═CHtBu(CH2tBu), where PNP = N[2-PiPr2-4-methylphenyl]2-. In addition to the initial activation of methane, a tautomerization reaction to a terminal methylidene is also explored due to methylidene's potential use as a synthetic starting point. Analogous complexes with other low-cost 3d transition metals were studied, such as scandium, titanium, vanadium, and chromium as both isoelectronic and isocharged complexes. Our results predict that V(IV) and V(V) complexes are promising for methane C-H bond activation. The V(V) complex has a low rate-determining barrier for methane activation, specifically 16.6 kcal/mol, which is approximately 12 kcal/mol less than that for the Ti complex, as well as having a moderate tautomerization barrier of 29.8 kcal/mol, while the V(IV) complex has a methane activation barrier of 19.0 kcal/mol and a tautomerization barrier of 31.1 kcal/mol. Scandium and chromium complexes are much poorer for C-H bond activation; scandium has very high barriers, while chromium strongly overstabilizes the alkylidene intermediate, potentially stopping the further reaction. In addition to the original PNP ligand, some of the most promising ligands from a previous work were tested, although (as shown previously) modification of the ligand does not typically have large effects on the activity of the system. Our best ligand modification improves the performance of the V(V) complex via the substitution of the nitrogen in PNP by phosphorus, which reduces the tautomerization barrier by 5 to 24.4 kcal/mol.

The Tetrel Bond and Tetrel Halide Perovskite Semiconductors

The ion pairs [Cs+•TtX3-] (Tt = Pb, Sn, Ge; X = I, Br, Cl) are the building blocks of all-inorganic cesium tetrel halide perovskites in 3D, CsTtX3, that are widely regarded as blockbuster materials for optoelectronic applications such as in solar cells. The 3D structures consist of an anionic inorganic tetrel halide framework stabilized by the cesium cations (Cs+). We use computational methods to show that the geometrical connectivity between the inorganic monoanions, [TtX3-]∞, that leads to the formation of the TtX64- octahedra and the 3D inorganic perovskite architecture is the result of the joint effect of polarization and coulombic forces driven by alkali and tetrel bonds. Depending on the nature and temperature phase of these perovskite systems, the Tt···X tetrel bonds are either indistinguishable or somehow distinguishable from Tt-X coordinate bonds. The calculation of the potential on the electrostatic surface of the Tt atom in molecular [Cs+•TtX3-] provides physical insight into why the negative anions [TtX3-] attract each other when in close proximity, leading to the formation of the CsTtX3 tetrel halide perovskites in the solid state. The inter-molecular (and inter-ionic) geometries, binding energies, and charge density-based topological properties of sixteen [Cs+•TtX3-] ion pairs, as well as some selected oligomers [Cs+•PbI3-]n (n = 2, 3, 4), are discussed.

Hydrogen Dissociation Reaction on First-Row Transition Metal Doped Nanobelts

Zigzag molecular nanobelts have recently captured the interest of scientists because of their appealing aesthetic structures, intriguing chemical reactivities, and tantalizing features. In the current study, first-row transition metals supported on an H6-N3-belt[6]arene nanobelt are investigated for the electrocatalytic properties of these complexes for the hydrogen dissociation reaction (HDR). The interaction of the doped transition metal atom with the nanobelt is evaluated through interaction energy analysis, which reveals the significant thermodynamic stability of TM-doped nanobelt complexes. Electronic properties such as frontier molecular orbitals and natural bond orbitals analyses are also computed, to estimate the electronic perturbation upon doping. The highest reduction in the HOMO-LUMO energy gap compared to the bare nanobelt is seen in the case of the Zn@NB catalyst (4.76 eV). Furthermore, for the HDR reaction, the Sc@NB catalyst displays the best catalytic activity among the studied catalysts, with a hydrogen dissociation barrier of 0.13 eV, whereas the second-best catalytic activity is observed for the Zn@NB catalyst (0.36 eV). It is further found that multiple active sites, i.e., the presence of the metal atom and nitrogen atom moiety, help to facilitate the dissociation of the hydrogen molecule. These key findings of this study enhance the understanding of the relative stability, electronic features, and catalytic bindings of various TM@NB catalysts.

Mechanistic Investigation of the Synthesis of Dianionic In-Derived Coordination Polymers

The mechanism of formation of crystalline coordination polymers is as complex as the architectures themselves. In this Communication, we detail a three-tiered approach using density functional theory (DFT) analysis, synthesis, and in situ Raman spectroscopy to study the formation of coordination polymers. Specifically, the previously reported coordination polymers YCM-22 and YCM-51 containing the [In(CO₂R)₂X₃]^2- (X = halogen) molecular building unit (MBU) were investigated. DFT revealed two potential pathways of formation, involving the initial formation of either [InCl₄]^- or [In(CO₂R)Cl₃]^-. A molecular dimeric In species (8a) containing two [In(CO₂R)Cl₄]^2- centers bridged by 2,5-thiophenedicarboxylic acid was isolated. When a suspension of 8a was treated with a solution of 2,5-thiophenedicarboxylic acid, an isomer of the coordination polymer YCM-22 (denoted as YCM-22') was formed. In situ Raman analysis of the formation of YCM-22 confirms that [InCl₄]^- forms at the onset of the reaction and that the [In(CO₂R)₂X₃]^2- MBU forms at its expense. The totality of the data presented support a mechanism of formation of one-dimensional In-derived coordination polymers and present a roadmap for future investigations into the formation of other crystalline coordination polymers.

Halogen bonds with carbenes acting as Lewis base units: complexes of imidazol-2-ylidene: theoretical analysis and experimental evidence

ωB97XD/aug-cc-pVDZ and ωB97XD/aug-cc-pVTZ calculations were performed on complexes of imidazol-2-ylidene that are linked by halogen bonds. This singlet carbene acts as the Lewis base through a lone electron pair located at the carbon centre. The XCCH, XCN and X₂ units were chosen here as those that interact through the X Lewis acid halogen centre (X = Cl, Br and I); if X = F the complexes are linked by interactions which are not classified as halogen bonds. The properties of interactions that occur in complexes are analyzed using the results of DFT calculations which are supported by parameters derived from the Quantum Theory of Atoms in Molecules, QTAIM, and the Natural Bond Orbital, NBO, approaches. The energy decomposition analysis, EDA, applied here provided additional characteristics of interactions linking complexes analyzed. The majority of complexes are linked by the medium in strength and strong halogen bonds which often possess characteristics typical for covalent bonds. Searches through the Cambridge Structural Database were also performed and structures analogues to complexes analyzed theoretically were found, and these structures are also discussed in this study.

Benchmarking DFT functionals for photophysics of pyranoflavylium cations

An intense absorption, phosphorescence, a long triplet excited state lifetime and singlet oxygen generation capabilities are characteristics of pyranoflavylium cations, analogues to pyranoanthocyanidins originated in the maturation process of red wine. Such properties make these compounds potential photosensitizers to be applied in photodynamic therapy. In this context, the photophysical processes underlying that treatment critically depend on the electronic structure of the pyranoflavylium molecules. When employing density functional theory to describe the electronic structure of molecules, the choice of the most suitable functional is not trivial, and benchmark studies are needed to orient practitioners in the field. In this work, a benchmark of seven of the most commonly used density functionals in addressing the photophysical properties of a set of eight pyranoflavylium cations is reported. Ground and excited state geometries, molecular orbitals, and absorption, fluorescence and phosphorescence transition energies were calculated using density functional theory approaches, and evaluated and compared to experimental data and monoreferential wave function-based methodologies. Statistical analysis of the results indicates that global-hybrid functionals allow an excellent description of absorption and emission energies, with errors around 0.05 eV, while range-separated variants led to somewhat larger errors in the range 0.1-0.2 eV. In contrast, range-separated functionals display excellent phosphorescence energies with errors close to 0.05 eV, in this case global-hybrids showing increased discrepancies around 0.5-0.1 eV.

Interaction between [(η6-p-cym)M(H2O)3]2+ (MII = Ru, Os) or [(η5-Cp*)M(H2O)3]2+ (MIII = Rh, Ir) and Phosphonate Derivatives of Iminodiacetic Acid: A Solution Equilibrium and DFT Study

The pH-dependent binding strengths and modes of the organometallic [(η⁶-p-cym)M(H₂O)₃]²⁺ (M^II = Ru, Os; p-cym = 1-methyl-4-isopropylbenzene) or [(η⁵-Cp*)M(H₂O)₃]²⁺ (M^III = Rh, Ir; Cp* = pentamethylcyclopentadienyl anion) cations towards iminodiacetic acid (H₂Ida) and its biorelevant mono- and diphosphonate derivatives N-(phosphonomethyl)-glycine (H₃IdaP) and iminodi(methylphosphonic acid) (H₄Ida2P) was studied in an aqueous solution. The results showed that all three of the ligands form 1:1 complexes via the tridentate (O,N,O) donor set, for which the binding mode was further corroborated by the DFT method. Although with IdaP^3- and Ida2P^4- in mono- and bis-protonated species, where H⁺ might also be located at the non-coordinating N atom, the theoretical calculations revealed the protonation of the phosphonate group(s) and the tridentate coordination of the phosphonate ligands. The replacement of one carboxylate in Ida^2- by a phosphonate group (IdaP^3-) resulted in a significant increase in the stability of the metal complexes; however, this increase vanished with Ida2P^4-, which was most likely due to some steric hindrance upon the coordination of the second large phosphonate group to form (5 + 5) joined chelates. In the phosphonate-containing systems, the neutral 1:1 complexes are the major species at pH 7.4 in the millimolar concentration range that is supported by both NMR and ESI-TOF-MS.

Hybrid nucleobases as new and efficient unnatural genetic letters

To expand the existing genetic letters beyond the natural four nucleotides, such as G, C, A, and T, it is necessary to design robust nucleotides that can not only produce stable and unperturbed DNA but also function naturally in living cells. Although hydrophobic bases, such as d5SICS (2,6-dimethyl-2H-isoquiniline-1-thione) and dNaM (2-methoxy-3-methylnaphthalene) were shown to be replicated in bacterial cells, the d5SICS:dNaM base-pair was found to perturb the structure of the duplex DNA. Therefore, it is necessary to design nucleobases that can form base pairs like the natural G:C and A:T pairs. Here, a reliable dispersion-corrected density functional theory has been used to design several nucleobases that can produce three-hydrogen-bonded base pairs like the G:C pair. In doing so, the Watson-Crick faces of d5SICS and dNaM were modified by replacing the hydrophobic groups with hydrogen bond donors and acceptors. As dNaM contains an unnatural C-glycosidic bond (C-dNaM), it was also modified to contain the natural N-glycosidic bond (N-dNaM). This technique produced 91 new bases (N-d5SICS-X (X = 1-33), C-dNaM-X (X = 1-35), and N-dNaM-X (X = 1-23), where X is the different types of modifications applied to d5SICS and dNaM) and 259 base-pairs. Among these base pairs, 76 base pairs are found to be more stable than the G:C pair. Interestingly, the N-d5SICS-32:C-dNaM-32 and N-d5SICS-32:N-dNaM-20 pairs are found to be the most stable with binding energies of about -28.0 kcal/mol. The base-pair patterns of these pairs are also analogous to that of the G:C pair. Hence, it is proposed that N-d5SICS-32, C-dNaM-32, and N-dNaM-20 would act as efficient new genetic letters to produce stable and unperturbed artificial DNA.Communicated by Ramaswamy H. Sarma.

Electrophilicity of Hoveyda-Grubbs Olefin Metathesis Catalysts as the Driving Force that Controls Initiation Rates

The dissociative mechanism of initiation for a series of Hoveyda-Grubbs type metathesis catalysts modified at the para and meta positions in the isopropoxybenzylidene ligand is investigated by means of DFT calculations. The electron donating/withdrawing capacity of the ligand was screened through the incorporation of various substituents such as halogens, nitro, alkoxides, ketones, esters, amines, and amides. Variations in structural parameters, energy barriers for the Ru-O bond dissociation, and Ru-O bond strength were examined as a function of the Hammett constant. It was found that electronic properties of the catalysts such as chemical potential, hardness, and electrophilicity correlate linearly with the dissociative energy barriers. These findings enable a systematic rationalization and prediction of rate of precatalyst initiation through the calculation of only the HOMO-LUMO gap of catalysts, as the faster the initiation, the more electrophilic the catalyst.

Elusive intermediates in cisplatin reaction with target amino acids: Platinum(II)-cysteine complexes assayed by IR ion spectroscopy and DFT calculations

The reactivity of a widely used metal based antineoplastic drug, cisplatin, cis-PtCl₂(NH₃)₂, with L-cysteine (Cys) has been investigated using a combination of electrospray ionization mass spectrometry (ESI-MS), IRMPD gas phase ion spectroscopy and DFT calculations. The cysteine lateral chain represents one of the main platination sites in proteins, which is believed to be related to the resistance mechanisms to cisplatin. The vibrational features of the mass-selected substitution product cis-[PtCl(NH₃)₂(Cys)]⁺ and the intercepted cis-[PtCl(NH₃)₂(H₂O)(Cys)]⁺ intermediate complex were compared to calculated IR spectra, enabling the assessment of the sampled ions structures. In cis-[PtCl(NH₃)₂(Cys)]⁺, cysteine was found to bind platinum through the sulfur atom as a thiolate zwitterion, highlighting the enhanced acidity of the cysteine thiol group upon metal coordination. The cis-[PtCl(NH₃)₂(H₂O)(Cys)]⁺ structure complies with the non-covalent encounter complex, formed by cis-[PtCl(NH₃)₂(H₂O)]⁺ and neutral cysteine. This species is able to undergo the substitution process to produce cis-[PtCl(NH₃)₂(Cys)]⁺ when activated as a mass-isolated ion suggesting its participation in the reaction mechanism of cisplatin with cysteine in solution. Finally, the DFT-calculated energy profile for the substitution reaction was correlated with the peculiar gas-phase reactivity of this non-covalent complex, resulting to be 10-fold less reactive toward substitution than the corresponding methionine complex.

Racemization Pathway for MoO2(acac)2 Favored over Ray-Dutt, Bailar, and Conte-Hippler Twists

Chiral cis-MoO2(acac)2 racemizes via four pathways that agree with and extend upon Muetterties' topological analysis for dynamic MX2(chel)2 complexes. Textbook Ray-Dutt and Bailar twists are the least favored with barriers of 27.5 and 28.7 kcal/mol, respectively. Rotating both acac ligands of the Bailar structure by 90° gives the lower Conte-Hippler twist (20.0 kcal/mol), which represents a valley-ridge inflection that invokes the trans isomer. The most favorable is a new twist that was found by 90° rotation of only one acac ligand of the Bailar structure. The gas-phase barrier of 17.4 kcal/mol for this Dhimba-Muller-Lammertsma twist further decreases upon inclusion of the effects of solvents to 16.3 kcal/mol (benzene), 16.2 kcal/mol (toluene), and 15.4 kcal/mol (chloroform), which are in excellent agreement with the reported experimental values.

Enhancing Effects of the Cyano Group on the C-X∙∙∙N Hydrogen or Halogen Bond in Complexes of X-Cyanomethanes with Trimethyl Amine: CH3−n(CN)nX∙∙∙NMe3, (n = 0–3; X = H, Cl, Br, I)

In this paper, density functional theory and wave function theory calculations are carried out to investigate the strength and nature of the intermolecular C-X∙∙∙N bond interaction as a function of the number of cyano groups, CN, in the X-bond donor while maintaining the X-bond acceptor as fixed. Specifically, complexes of X-cyanomethanes with trimethyl amine CH_3−n(CN)_nX∙∙∙NMe₃ (n = 0–3; X = H, Cl, Br, I) are used as model systems. Geometrical parameters and vibrational C-X-stretching frequencies as well as interaction energies are used as relevant indicators to gauge hydrogen or halogen bond strength in the complexes. Additional characteristics of interactions that link these complexes, i.e., hydrogen or halogen bonds, are calculated with the use of the following theoretical tools: the atoms in molecules (AIM) approach, the natural bond orbital (NBO) method, and energy decomposition analysis (EDA). The results show that, for the specified X-center, the strength of C-X∙∙∙N interaction increases significantly and in a non-additive fashion with the number of CN groups. Moreover, the nature (noncovalent or partly covalent) of the interactions is revealed via the AIM approach.

Sandwich, Triple-Decker and Other Sandwich-like Complexes of Cyclopentadienyl Anions with Lithium or Sodium Cations

Density functional theory, DFT, calculations were carried out on complexes containing cyclopentadienyl anions and lithium or sodium cations; half-sandwich, sandwich and sandwich-like complexes (among them triple-decker ones) are analyzed. Searches performed through the Cambridge Structural Database revealed that crystal structures containing these motifs exist, mostly structures with sodium cations. The DFT calculations performed here include geometry optimization and frequency calculations of the complexes at the ωB97XD/aug-cc-pVTZ level, followed by the partitioning of the energy of interaction via the Energy Decomposition Analysis scheme, EDA, at the BP86-D3/TZ2P level. Additional calculations and analyses were performed using both the Quantum Theory of Atoms in Molecules, QTAIM, and the Natural Bond Orbital analyses, NBO. The results of this work show that the electrostatic interaction energy is the most important attractive contribution to the total interaction energy of each of the complex systems analyzed here, and that complexation itself leads to minor electron charge shifts.

Structure of Aqueous KNO3 Solutions by Wide-Angle X-ray Scattering and Density Functional Theory

The structure of aqueous KNO₃ solutions was studied by wide-angle X-ray scattering (WAXS) and density functional theory (DFT). The interference functions were subjected to empirical potential structure refinement (EPSR) modeling to extract the detailed hydration structure information. In aqueous KNO₃ solutions with a concentration from 0.50 to 2.72 mol·dm^-3, water molecules coordinate K⁺ with a mean coordination number (CN) from 6.6 ± 0.9 to 5.8 ± 1.2, respectively, and a K-O_W(H₂O) distance of 2.82 Å. To further describe the hydration properties of K⁺, a hydration factor (f_h) was defined based on the orientational angle between the water O-H vector and the ion-oxygen vector. The f_h value obtained for K⁺ is 0.792, 0.787, 0.766, and 0.765, and the corresponding average hydration numbers (HN) are 5.2, 5.1, 4.8, and 4.5. The reduced HN is compensated by the direct binding of oxygen atoms of NO₃^- with the average CN from 0.3 ± 0.7 to 2.6 ± 1.7, respectively, and the K-O_N(NO₃^-) distance of 2.82 Å. The average number of water molecules around NO₃^- slightly reduces from 10.5 ± 1.5 to 9.6 ± 1.7 with r_N-OW = 3.63 Å. K⁺-NO₃^- ion association is characterized by K-O_N and K-N pair correlation functions (PCFs). A K-N peak is observed at 3.81 Å as the main peak with a shoulder at 3.42 Å in g_K-N(r). This finding indicates that two occupancy sites are available for K⁺, i.e., the edge-shared bidentate (BCIP) and the corner-shared monodentate (MCIP) contact ion pairs. The structure and stability of the KNO₃(H₂O)₁₀ cluster were investigated by a DFT method and cross-checked with the results from WAXS.

Preparation of Deuterium-Labeled Armodafinil by Hydrogen–Deuterium Exchange and Its Application in Quantitative Analysis by LC-MS

Armodafinil, the R enantiomer of modafinil, was approved in 2007 by the US Food and Drug Administration as a wake-promoting agent for excessive sleepiness treatment. Due to its abuse by students and athletes, there is a need of its quantification. Quantitative analysis by liquid chromatography-mass spectrometry, however, though very common and sensitive, frequently cannot be performed without isotopically labeled standards which usually have to be specially synthesized. Here we reported our investigation on the preparation of deuterated standard of armodafinil based on the simple and inexpensive hydrogen–deuterium exchange reaction at the carbon centers. The obtained results clearly indicate the possibility of introduction of three deuterons into the armodafinil molecule. The introduced deuterons do not undergo back exchange under neutral and acidic conditions. Moreover, the deuterated and non-deuterated armodafinil isotopologues revealed co-elution during the chromatographic analysis. The ability to control the degree of deuteration using different reaction conditions was determined. The proposed method of deuterated armodafinil standard preparation is rapid, cost-efficient and may be successfully used in its quantitative analysis by LC-MS.

η6-Coordinated ruthenabenzenes from three-component assembly on a diruthenium μ-allenyl scaffold

The room temperature reactions with internal alkynes, RCCR, of the μ-allenyl acetonitrile complex [Ru₂Cp₂(CO)₂(NCMe){μ-η¹:η²-C¹HC²C³Me₂}]BF₄ (1-NCMe), freshly prepared from the tricarbonyl precursor [Ru₂Cp₂(CO)₃{μ-η¹:η²-C¹HC²C³Me₂}]BF₄, 1, proceeded with alkyne insertion into ruthenium-allenyl bond and allenyl-CO coupling, affording compounds [Ru₂Cp₂(CO)₂{μ-η²:η⁵-C(R)C(R)C¹HC²(C³MeCH₂)C(OH)}]BF₄ (R = Ph, 2; R = CO₂Me, 3; R = CO₂Et, 4) in 83-94% yields. Deprotonation of 2-4 by triethylamine gave [Ru₂Cp₂(CO)₂{μ-η²:η⁵-C(R)C(R)CHC(CMeCH₂)C(O)}] (R = Ph, 5; R = CO₂Me, 6; R = CO₂Et, 7) in 75-88% yields, and 2-4 could be recovered upon HBF₄·Et₂O addition to 5-7. All the products, 2-7, were fully characterized by elemental analysis, IR and multinuclear NMR spectroscopy. The structure of 2 was ascertained by single crystal X-ray diffraction and investigated by DFT calculations, revealing a six-membered ruthenacycle with Shannon aromaticity index in line with related compounds. The formation of ruthenium-coordinated ruthenabenzenes from a preexistent diruthenium scaffold is a versatile but underdeveloped approach exploiting cooperative effects typical of a dimetallic core.

Mechanisms, Challenges, and Opportunities of Dual Ni/Photoredox-Catalyzed C(sp2)-C(sp3) Cross-Couplings

The merging of photoredox and nickel catalysis has revolutionized the field of C-C cross-coupling. However, in comparison to the development of synthetic methods, detailed mechanistic investigations of these catalytic systems are lagging. To improve the mechanistic understanding, computational tools have emerged as powerful tools to elucidate the factors controlling reactivity and selectivity in these complex catalytic transformations. Based on the reported computational studies, it appears that the mechanistic picture of catalytic systems is not generally applicable, but is rather dependent on the specific choice of substrate, ligands, photocatalysts, etc. Given the complexity of these systems, the need for more accurate computational methods, readily available and user-friendly dynamics simulation tools, and data-driven approaches is clear in order to understand at the molecular level the mechanisms of these transformations. In particular, we anticipate that such improvement of theoretical methods will become crucial to advance the understanding of excited-state properties and dynamics of key species, as well as to enable faster and unbiased exploration of reaction pathways. Further, with greater collaboration between computational, experimental, and spectroscopic communities, the mechanistic investigation of photoredox/Ni dual-catalytic reactions is expected to thrive quickly, facilitating the design of novel catalytic systems and promoting our understanding of the reaction selectivity.

Investigation of the Nature of Intermolecular Interactions in Tetra(thiocyanato)corrolato-Ag(III) Complexes: Agostic or Hydrogen Bonded?

Tetra(thiocyanato)corrolato-Ag(III) complexes presented here constitute a new class of metallo-corrole complexes. The spectroscopic properties of these complexes are quite unusual and interesting. For example, the absorption spectra of these β-substituted corrolato-Ag(III) complexes are very different from those of the β-unsubstituted corrolato-Ag(III) derivatives. Single-crystal XRD analysis of a representative tetra(thiocyanato)corrolato-Ag(III) derivative reveals C-H···Ag interactions. The C-H···Ag interactions are rarely demonstrated in the crystal lattice of a discrete coordination/organometallic compound. Optimization of the hydrogen positions of the crystal structure discloses the geometrical parameters of the said interaction as a Ag···H distance of 2.597 Å and ∠C-H···Ag of 109.62°. The natural bond orbital analysis provides information about the donor-acceptor orbitals involved in the interactions and their interaction energies. It was observed that the σ_C-H orbital overlaps with the vacant d-orbital of Ag with an interaction energy of 17.93 kJ/mol. The filled d-orbital of Ag overlaps with the σ*_C-H orbital with an interaction energy of 4.79 kJ/mol. The highlights of this work are that the H···Ag distance is outside of the distance range for the typical agostic interaction but fitted with the weak H-bond distance. However, the ∠C-H···Ag angle is within the range of the agostic interaction. Both crystallographic data and electronic structure calculations reveal that these kinds of intermolecular interactions in square-planar d⁸ Ag(III) complexes are intermediate in nature. Thus, they cannot be categorically called either hydrogen bonding or agostic interaction.

Silaborative Assembly of Allenamides and Alkynes: Highly Regio- and Stereocontrolled Access to Bi- or Trimetallic Skipped Dienes

A highly stereo- and regiocontrolled multicomponent approach to skipped 1,4-dienes decorated with one boryl and two silyl functionalities is described. This Pd-catalyzed atom-economical union of allenamides, alkynes, and Me₂ PhSiBpin (or Et₃ SiBpin) proceeds without the use of phosphine ligands, instead relying on chelation through the internal amide group of the allenamide sulfonyl. A variety of alkynes, including those derived from complex bioactive molecules, can be efficiently coupled with allenamides and Me₂ PhSiBpin in good yields and with excellent selectivity. The synthetic potential was demonstrated through multiple valuable chemoselective transformations, establishing new disconnections for functionalized dienes. Density functional theory calculations revealed that the reaction first proceeded through borylation of the allenamide, followed by silylation of the alkyne and then reductive elimination, which convergently assemble the skipped 1,4-diene.

Formations of bimolecular barbituric acid complexes through hydrogen bonding interactions: DFT analyses of structural and electronic features

Formations of bimolecular barbituric acid (BA) complexes through hydrogen-bonding (HB) interactions were investigated in this work. BA has been known as a starting compound of pharmaceutical compounds developments, in which the molecular and atomic features of parent BA in homo-paring with another BA molecule were investigated here. The models were optimized to reach the stabilized structures and their properties were evaluated at the molecular and atomic scales. Density functional theory (DFT) calculations were performed to provide required information for achieving the goal of this work. Six dimer models were obtained finally according to examining all possible starting dimers configurations for involving in optimization calculations. N-H . . . O and C-H . . . O interactions were also involved in dimers formations besides participation of the X-center of parent BA in interaction. Molecular and atomic scales features were evaluated for characterizing the dimers formations. As a consequence, several configurations of BA dimers were obtained showing the importance of performing such structural analyses for developing further compounds from BA.

Elucidating Two Distinct Pathways for Electrocatalytic Hydrogen Production Using CoII Pincer Complexes

Hydrogen gas is a sustainable energy source with water as the sole combustion product. As a result, efforts to catalyze H₂ production are pertinent and widespread. The electrocatalytic H₂ generating capabilities of two Co^II complexes, [Co(κ³ -2,6-{Ph₂ PNR}₂ (NC₅ H₃ ))Br₂ ] with R=H (I) or R=Me (II), were presented for a variety of proton sources including trifluoroacetic acid (TFA), acetic acid (AA), and trifluoroethanol (TFE). Cyclic voltammetry and controlled potential coulometry demonstrated that electrocatalysis from I and II occurred at two different potentials and are associated with different reduction processes. Density functional theory analysis provided insight into the identities of the catalyst and supported two distinct reaction pathways for electrocatalytic proton reduction. Specifically, stronger acids (e. g., AA, TFA) proceeded at -1.31 to -1.45 V through a M^I /M^III pathway while sources with higher pK_a values (e. g., TFE, H₂ O) generated hydrogen at -2.4 V via M⁰ /M^II ligand-assisted metal-centered reduction.