Theoretical Investigation of the Potential-Dependent CO Adsorption on Copper Electrodes.

Despite the importance of CO adsorption in many electrocatalytic reaction mechanisms, there has been little investigation of the dependence of the free energy of CO adsorption on the applied potential. Herein, we report on the potential-dependent adsorption of CO on Cu electrodes using a grand-canonical density functional theory approach. We demonstrate that, within the working potential range of electrocatalytic CO2 reduction on Cu(111) and Cu(100), the CO adsorption strength can change by over 0.1 eV. Our analyses explain the potential dependence through an interfacial capacitance loss upon CO adsorption as well as orbital relaxation induced by the electrode potential. Via sensitivity analysis with respect to two electrolyte model parameters (solvent dielectric constant and Debye screening length), we find that the surface excess charge density is a useful descriptor of the CO adsorption free energy.

Substituent Impact on Quinoxaline Performance and Degradation in Redox Flow Batteries.

Aqueous redox flow batteries (RFBs) are attractive candidates for low-cost, grid-scale storage of energy from renewable sources. Quinoxaline derivatives represent a promising but underexplored class of charge-storing materials on account of poor chemical stability in prior studies (with capacity fade rates >20%/day). Here, we establish that 2,3-dimethylquinoxaline-6-carboxylic acid (DMeQUIC) is vulnerable to tautomerization in its reduced form under alkaline conditions. We obtain kinetic rate constants for tautomerization by applying Bayesian inference to ultraviolet-visible spectroscopic data from operating flow cells and show that these rate constants quantitatively account for capacity fade measured in cycled cells. We use density functional theory (DFT) modeling to identify structural and chemical predictors of tautomerization resistance and demonstrate that they qualitatively explain stability trends for several commercially available and synthesized derivatives. Among these, quinoxaline-2-carboxylic acid shows a dramatic increase in stability over DMeQUIC and does not exhibit capacity fade in mixed symmetric cell cycling. The molecular design principles identified in this work set the stage for further development of quinoxalines in practical, aqueous organic RFBs.

High-Performance Iridium-Molybdenum Oxide Electrocatalysts for Water Oxidation in Acid: Bayesian Optimization Discovery and Experimental Testing.

Ir oxides are costly and scarce catalysts for oxygen evolution reaction (OER) in acid. There has been extensive interest in developing alternatives that are either Ir-free or require smaller amounts of Ir to drive the reactions at acceptable rates. One design strategy is to identify Ir-based mixed oxides that achieve similar performance while requiring smaller amounts of Ir. The obstacle to this strategy has been a very large phase space of the Ir-based mixed metal oxides, in terms of the metals combined with Ir and the different crystallographic structures of the mixed oxides, which prevents a thorough exploration of possible materials. In this work, we developed a workflow that uses machine-learning-aided Bayesian optimization in combination with density functional theory to make the exploration of this phase space plausible. This screening identified Mo as a promising dopant for forming acid-tolerant Ir-based oxides for the OER. We synthesized and characterized the Ir-Mo mixed oxides in the form of thin-film electrocatalysts with a known surface area. We show that these mixed oxides exhibited overpotentials ∼30 mV lower than a pure Ir control while maintaining 24% lower Ir dissolution rates than the Ir control. These findings suggest that Mo is a promising dopant and highlight the promise of machine learning to guide the experimental exploration and optimization of catalytic materials.

Understanding capacity fade in organic redox-flow batteries by combining spectroscopy with statistical inference techniques.

Organic redox-active molecules are attractive as redox-flow battery (RFB) reactants because of their low anticipated costs and widely tunable properties. Unfortunately, many lab-scale flow cells experience rapid material degradation (from chemical and electrochemical decay mechanisms) and capacity fade during cycling (>0.1%/day) hindering their commercial deployment. In this work, we combine ultraviolet-visible spectrophotometry and statistical inference techniques to elucidate the Michael attack decay mechanism for 4,5-dihydroxy-1,3-benzenedisulfonic acid (BQDS), a once-promising positive electrolyte reactant for aqueous organic redox-flow batteries. We use Bayesian inference and multivariate curve resolution on the spectroscopic data to derive uncertainty-quantified reaction orders and rates for Michael attack, estimate the spectra of intermediate species and establish a quantitative connection between molecular decay and capacity fade. Our work illustrates the promise of using statistical inference to elucidate chemical and electrochemical mechanisms of capacity fade in organic redox-flow battery together with uncertainty quantification, in flow cell-based electrochemical systems.

Explaining the structure sensitivity of Pt and Rh for aqueous-phase hydrogenation of phenol.

Phenol is an important model compound to understand the thermocatalytic (TCH) and electrocatalytic hydrogenation (ECH) of biomass to biofuels. Although Pt and Rh are among the most studied catalysts for aqueous-phase phenol hydrogenation, the reason why certain facets are active for ECH and TCH is not fully understood. Herein, we identify the active facet of Pt and Rh catalysts for aqueous-phase hydrogenation of phenol and explain the origin of the size-dependent activity trends of Pt and Rh nanoparticles. Phenol adsorption energies extracted on the active sites of Pt and Rh nanoparticles on carbon by fitting kinetic data show that the active sites adsorb phenol weakly. We predict that the turnover frequencies (TOFs) for the hydrogenation of phenol to cyclohexanone on Pt(111) and Rh(111) terraces are higher than those on (221) stepped facets based on density functional theory modeling and mean-field microkinetic simulations. The higher activities of the (111) terraces are due to lower activation energies and weaker phenol adsorption, preventing high coverages of phenol from inhibiting hydrogen adsorption. We measure that the TOF for ECH of phenol increases as the Rh nanoparticle diameter increases from 2 to 10 nm at 298 K and -0.1 V vs the reversible hydrogen electrode, qualitatively matching prior reports for Pt nanoparticles. The increase in experimental TOFs as Pt and Rh nanoparticle diameters increase is due to a larger fraction of terraces on larger particles. These findings clarify the structure sensitivity and active site of Pt and Rh for the hydrogenation of phenol and will inform the catalyst design for the hydrogenation of bio-oils.

Unveiling the Cerium(III)/(IV) Structures and Charge-Transfer Mechanism in Sulfuric Acid.

The Ce3+/Ce4+ redox couple has a charge transfer (CT) with extreme asymmetry and a large shift in redox potential depending on electrolyte composition. The redox potential shift and CT behavior are difficult to understand because neither the cerium structures nor the CT mechanism are well understood, limiting efforts to improve the Ce3+/Ce4+ redox kinetics in applications such as energy storage. Herein, we identify the Ce3+ and Ce4+ structures and CT mechanism in sulfuric acid via extended X-ray absorption fine structure spectroscopy (EXAFS), kinetic measurements, and density functional theory (DFT) calculations. We show EXAFS evidence that confirms that Ce3+ is coordinated by nine water molecules and suggests that Ce4+ is complexed by water and three bisulfates in sulfuric acid. Despite the change in complexation within the first coordination shell between Ce3+ and Ce4+, we show that the kinetics are independent of the electrode, suggesting outer-sphere electron-transfer behavior. We identify a two-step mechanism where Ce4+ exchanges the bisulfate anions with water in a chemical step followed by a rate-determining electron transfer step that follows Marcus theory (MT). This mechanism is consistent with all experimentally observed structural and kinetic data. The asymmetry of the Ce3+/Ce4+ CT and the observed shift in the redox potential with acid is explained by the addition of the chemical step in the CT mechanism. The fitted parameters from this rate law qualitatively agree with DFT-predicted free energies and the reorganization energy. The combination of a two-step mechanism with MT should be considered for other metal ion CT reactions whose kinetics have not been appropriately described.

Why halides enhance heterogeneous metal ion charge transfer reactions.

The reaction kinetics of many metal redox couples on electrode surfaces are enhanced in the presence of halides (i.e., Cl-, Br-, I-). Using first-principles metadynamics simulations, we show a correlation between calculated desorption barriers of V3+-anion complexes bound to graphite via an inner-sphere anion bridge and experimental V2+/V3+ kinetic measurements on edge plane pyrolytic graphite in H2SO4, HCl, and HI. We extend this analysis to V2+/V3+, Cr2+/Cr3+, and Cd0/Cd2+ reactions on a mercury electrode and demonstrate that reported kinetics in acidic electrolytes for these redox couples also correlate with the predicted desorption barriers of metal-anion complexes. Therefore, the desorption barrier of the metal-anion surface intermediate is a descriptor of kinetics for many metal redox couple/electrode combinations in the presence of halides. Knowledge of the metal-anion surface intermediates can guide the design of electrolytes and electrocatalysts with faster kinetics for redox reactions of relevance to energy and environmental applications.

Structures and Free Energies of Cerium Ions in Acidic Electrolytes.

The Ce3+/Ce4+ redox potential changes with the electrolyte, which could be due to unequal anion complexation free energies between Ce3+ and Ce4+ or a change in the solvent electrostatic screening. Ce complexation with anions and solvent screening also affect the solubility of Ce and charge transfer kinetics for electrochemical reactions involving waste remediation and energy storage. We report the structures and free energies of cerium complexes in seven acidic electrolytes based on Extended X-ray Absorption Fine Structure, UV-vis, and Density Functional Theory calculations. Ce3+ coordinates with nine water molecules as [Ce(H2O)9]3+ in all studied electrolytes. However, Ce4+ complexes with anions in all electrolytes except HClO4. Thus, our results suggest that Ce4+-anion complexation leads to the large shifts in standard redox potential. Long range screening effects are smaller than the anion complexation energies but could be responsible for changes in the Ce solubility with acid.

Inorganic Halide Double Perovskites with Optoelectronic Properties Modulated by Sublattice Mixing.

All-inorganic halide double perovskites have emerged as a promising class of materials that are potentially more stable and less toxic than lead-containing hybrid organic-inorganic perovskite optoelectronic materials. In this work, 311 cesium chloride double perovskites (Cs2BB'Cl6) were selected from a set of 903 compounds as likely being stable on the basis of a statistically learned tolerance factor (τ) for perovskite stability. First-principles calculations on these 311 double perovskites were then performed to assess their stability and identify candidates with band gaps appropriate for optoelectronic applications. We predict that 261 of the 311 Cs2BB'Cl6 compounds are likely synthesizable on the basis of a thermodynamic analysis of their decomposition to competing compounds (decomposition enthalpy <0.05 eV/atom). Of these 261 likely synthesizable compounds, 47 contain no toxic elements and have direct or nearly direct (within 100 meV) band gaps between 1 and 3 eV, as computed with hybrid density functional theory (HSE06). Within this set, we identify the triple-alkali perovskites Cs2[Alk]+[TM]3+Cl6, where Alk is a group 1 alkali cation and TM is a transition-metal cation, as a class of Cs2BB'Cl6 double perovskites with remarkable optical properties, including large and tunable exciton binding energies as computed by the GW-Bethe-Salpeter equation (GW-BSE) method. We attribute the unusual electronic structure of these compounds to the mixing of the Alk-Cl and TM-Cl sublattices, leading to materials with small band gaps, large exciton binding energies, and absorption spectra that are strongly influenced by the identity of the transition metal. The role of the double-perovskite structure in enabling these unique properties is probed through an analysis of the electronic structures and chemical bonding of these compounds in comparison with other transition-metal and alkali transition-metal halides.

Surpassing the single-atom catalytic activity limit through paired Pt-O-Pt ensemble built from isolated Pt1 atoms.

Despite the maximized metal dispersion offered by single-atom catalysts, further improvement of intrinsic activity can be hindered by the lack of neighboring metal atoms in these systems. Here we report the use of isolated Pt1 atoms on ceria as "seeds" to develop a Pt-O-Pt ensemble, which is well-represented by a Pt8O14 model cluster that retains 100% metal dispersion. The Pt atom in the ensemble is 100-1000 times more active than their single-atom Pt1/CeO2 parent in catalyzing the low-temperature CO oxidation under oxygen-rich conditions. Rather than the Pt-O-Ce interfacial catalysis, the stable catalytic unit is the Pt-O-Pt site itself without participation of oxygen from the 10-30 nm-size ceria support. Similar Pt-O-Pt sites can be built on various ceria and even alumina, distinguishable by facile activation of oxygen through the paired Pt-O-Pt atoms. Extending this design to other reaction systems is a likely outcome of the findings reported here.

New tolerance factor to predict the stability of perovskite oxides and halides.

Predicting the stability of the perovskite structure remains a long-standing challenge for the discovery of new functional materials for many applications including photovoltaics and electrocatalysts. We developed an accurate, physically interpretable, and one-dimensional tolerance factor, τ, that correctly predicts 92% of compounds as perovskite or nonperovskite for an experimental dataset of 576 ABX 3 materials (X = O2-, F-, Cl-, Br-, I-) using a novel data analytics approach based on SISSO (sure independence screening and sparsifying operator). τ is shown to generalize outside the training set for 1034 experimentally realized single and double perovskites (91% accuracy) and is applied to identify 23,314 new double perovskites (A 2 BB'X 6) ranked by their probability of being stable as perovskite. This work guides experimentalists and theorists toward which perovskites are most likely to be successfully synthesized and demonstrates an approach to descriptor identification that can be extended to arbitrary applications beyond perovskite stability predictions.

A Cu25 Nanocluster with Partial Cu(0) Character.

Atomically precise copper nanoclusters (NCs) are of immense interest for a variety of applications, but have remained elusive. Herein, we report the isolation of a copper NC, [Cu25H22(PPh3)12]Cl (1), from the reaction of Cu(OAc) and CuCl with Ph2SiH2, in the presence of PPh3. Complex 1 has been fully characterized, including analysis by X-ray crystallography, XANES, and XPS. In the solid state, complex 1 is constructed around a Cu13 centered-icosahedron and formally features partial Cu(0) character. XANES of 1 reveals a Cu K-edge at 8979.6 eV, intermediate between the edge energies of Cu(0) and Cu(I), confirming our oxidation state assignment. This assignment is further corroborated by determination of the Auger parameter for 1, which also falls between those recorded for Cu(0) and Cu(I).

Rate-Enhancing Roles of Water Molecules in Methyltrioxorhenium-Catalyzed Olefin Epoxidation by Hydrogen Peroxide.

Olefin epoxidation catalyzed by methyltrioxorhenium (MTO, CH3ReO3) is strongly accelerated in the presence of H2O. The participation of H2O in each of the elementary steps of the catalytic cycle, involving the formation of the peroxo complexes (CH3ReO2(η(2)-O2), A, and CH3ReO(η(2)-O2)2(H2O), B), as well as in their subsequent epoxidation of cyclohexene, was examined in aqueous acetonitrile. Experimental measurements demonstrate that the epoxidation steps exhibit only weak [H2O] dependence, attributed by DFT calculations to hydrogen bonding between uncoordinated H2O and a peroxo ligand. The primary cause of the observed H2O acceleration is the strong co-catalytic effect of water on the rates at which A and B are regenerated and consequently on the relative abundances of the three interconverting Re-containing species at steady state. Proton transfer from weakly coordinated H2O2 to the oxo ligands of MTO and A, resulting in peroxo complex formation, is directly mediated by solvent H2O molecules. Computed activation parameters and kinetic isotope effects, in combination with proton-inventory experiments, suggest a proton shuttle involving one or (most favorably) two H2O molecules in the key ligand-exchange steps to form A and B from MTO and A, respectively.

Synthesis and characterization of a Cu14 hydride cluster supported by neutral donor ligands.

The copper hydride clusters [Cu14 H12 (phen)6(PPh3)4][X]2 (X=Cl or OTf; OTf=trifluoromethanesulfonate, phen=1,10-phenanthroline) are obtained in good yields by the reaction of [(Ph3P)CuH]6 with phen, in the presence of a halide or pseudohalide source. The complex [Cu14H12 (phen)6(PPh3)4][Cl]2 reacts with CO2 in CH2Cl2 , in the presence of excess Ph3 P, to form the formate complex [(Ph3P)2Cu(κ(2)-O2CH)], along with [(phen)(Ph3 P)CuCl].

Water-catalyzed activation of H2O2 by methyltrioxorhenium: a combined computational-experimental study.

The formation of peroxorhenium complexes by activation of H2O2 is key in selective oxidation reactions catalyzed by CH3ReO3 (methyltrioxorhenium, MTO). Previous reports on the thermodynamics and kinetics of these reactions are inconsistent with each other and sometimes internally inconsistent. New experiments and calculations using density functional theory with the ωB97X-D and augmented def2-TZVP basis sets were conducted to better understand these reactions and to provide a strong experimental foundation for benchmarking computational studies involving MTO and its derivatives. Including solvation contributions to the free energies as well as tunneling corrections, we compute negative reaction enthalpies for each reaction and correctly predict the hydration state of all complexes in aqueous CH3CN. New rate constants for each of the forward and reverse reactions were both measured and computed as a function of temperature, providing a complete set of consistent activation parameters. New, independent measurements of equilibrium constants do not indicate strong cooperativity in peroxide ligand binding, as was previously reported. The free energy barriers for formation of both CH3ReO2(η(2)-O2) (A) and CH3ReO(η(2)-O2)2(H2O) (B) are predominantly entropic, and the former is much smaller than a previously reported value. Computed rate constants for a direct ligand-exchange mechanism, and for a mechanism in which a water molecule facilitates ligand-exchange via proton transfer in the transition state, differ by at least 7 orders of magnitude. The latter, water-assisted mechanism is predicted to be much faster and is consequently in much closer agreement with the experimentally measured kinetics. Experiments confirm the predicted catalytic role of water: the kinetics of both steps are strongly dependent on the water concentration, and water appears directly in the rate law.

Isolated catalyst sites on amorphous supports: a systematic algorithm for understanding heterogeneities in structure and reactivity.

Methods for modeling catalytic sites on amorphous supports lag far behind methods for modeling catalytic sites on metal surfaces, zeolites, and other crystalline materials. One typical strategy for amorphous supports uses cluster models with arbitrarily chosen constraints to model the rigid amorphous support, but these constraints arbitrarily influence catalyst site activity. An alternative strategy is to use no constraints, but this results in catalytic sites with unrealistic flexibility. We present a systematic ab initio method to model isolated active sites on insulating amorphous supports using small cluster models. A sequential quadratic programming framework helps us relate chemical properties, such as the activation energy, to active site structure. The algorithm is first illustrated on an empirical valence bond model energy landscape. We then use the algorithm to model an off-pathway kinetic trap in olefin metathesis by isolated Mo sites on amorphous SiO2. The cluster models were terminated with basis set deficient fluorine atoms to mimic the properties of an extended silica framework. We also discuss limitations of the current algorithm formulation and future directions for improvement.

Synthesis and micropatterning of photocatalytically reactive self-assembled monolayers covalently linked to Si(100) surfaces via a Si-C bond.

Selective generation of an amine-terminated self-assembled monolayer bound to silicon wafers via a silicon-carbon linkage was realized by photocatalytically reducing the corresponding azide-terminated, self-assembled monolayers (Az-SAMs). The Az-SAM was obtained by thermal deposition of 11-chloroundecene onto a hydrogen-terminated silicon wafer followed by nucleophilic substitution of the chloride with the azide ion in warm N,N'-dimethylformamide (DMF). The presence of the terminal azide group on the SAM was confirmed by reflection absorption infrared spectroscopy (RAIRS), by X-ray photoelectron spectroscopy (XPS), and by detecting the formation of a triazole upon reaction of the azide with an activated alkyne. The desired terminal amine groups were generated by photocatalytic reduction of the Az-SAM with cadmium selenide quantum dots (CdSe Qdots) using λ > 400 nm. Analysis of the reduced SAM by XPS gave results that were consistent with those obtained with an amine-terminated surface obtained by reducing the Az-SAM with triphenylphosphine. To demonstrate the feasibility of using the Az-SAM for surface patterning, a sample was coated with adsorbed CdSe Qdots and exposed to the output of a diode laser at λ = 407 nm through a micropatterned mask. Using a SEM, the pattern formed in this manner was revealed after removing the CdSe Qdots and subsequently adsorbing 10 nm gold nanoparticles (AuNPs) to the positively charged terminal-amine groups. The formation of the pattern by CdSe-photocatalyzed reduction of the azide demonstrates a novel route to create features by selective modification of organic monolayers on silicon wafers.