Combining Renewable Electricity and Renewable Carbon: Understanding Reaction Mechanisms of Biomass-Derived Furanic Compounds for Design of Catalytic Nanomaterials.

ConspectusDespite the growing deployment of renewable energy conversion technologies, a number of large industrial sectors remain challenging to decarbonize. Aviation, heavy transport, and the production of steel, cement, and chemicals are heavily dependent on carbon-containing fuels and feedstocks. A hopeful avenue toward carbon neutrality is the implementation of renewable carbon for the synthesis of critical fuels, chemicals, and materials. Biomass provides an opportune source of renewable carbon, naturally capturing atmospheric CO2 and forming multicarbon linkages and useful chemical functional groups. The constituent molecules nonetheless require various chemical transformations, often best facilitated by catalytic nanomaterials, in order to access usable final products.Catalyzed transformations of renewable biomass compounds may intersect with renewable energy production by offering a means to utilize excess intermittent electricity and store it within chemical bonds. Electrochemical catalytic processes can often offer advantages in energy efficiency, product selectivity, and modular scalability compared to thermal-driven reactions. Electrocatalytic reactions with renewable carbon feedstocks can further enable related processes such as water splitting, where value-adding organic oxidation reactions may replace the evolution of oxygen. Organic electroreduction reactions may also allow desirable hydrogenations of bonds without intermediate formation of H2 and need for additional reactors.This Account highlights recent work aimed at gaining a fundamental understanding of transformations involving biomass-derived molecules in electrocatalytic nanomaterials. Particular emphasis is placed on the oxidation of biomass derived furanic compounds such as furfural and 5-hydroxymethylfurfural (HMF), which can yield value-added chemicals, including furoic acid (FA), maleic acid (MA), and 2,5-furandicarboxylic acid (FDCA) for renewable materials and other commodities. We highlight advanced implementations of online electrochemical mass spectrometry (OLEMS) and vibrational spectroscopies such as attenuated total reflectance surface enhanced infrared reflection absorption spectroscopy (ATR-SEIRAS), combined with microkinetic models (MKMs) and quantum chemical calculations, to shed light on the elementary mechanistic pathways involved in electrochemical biomass conversion and how these paths are influenced by catalytic nanomaterials. Perspectives are given on the potential opportunities for materials development toward more efficient and selective carbon-mitigating reaction pathways.

Tight-binding model predicts exciton energetics and structure for photovoltaic molecules.

Conjugated molecules and polymers are being designed as acceptor and donor materials for organic photovoltaic (OPV) cells. OPV performance depends on generation of free charge carriers through dissociation of excitons, which are electron-hole pairs created when a photon is absorbed. Here, we develop a tight-binding model to describe excitons on homo-oligomers, alternating co-oligomers, and a non-fullerene acceptor - IDTBR. We parameterize our model using density functional theory (DFT) energies of neutral, anion, cation, and excited states of constituent moieties. A symmetric molecule like IDTBR has two ends where an exciton can sit; but the product wavefunction approximation for the exciton breaks symmetry. So, we introduce a tight-binding model with full correlation between electron and hole, which allows the exciton to coherently explore both ends of the molecule. Our approach predicts optical singlet excitation energies for oligomers of varying length as well as IDTBR in good agreement with time-dependent DFT and spectroscopic results.

Negative Dielectric Constant of Water at a Metal Interface.

Water polarizability at a metal interface plays an essential role in electrochemistry. We devise a classical molecular dynamics approach with an efficient description of metal polarization and a novel ac field method to measure the local dielectric response of interfacial water. Water adlayers next to the metal surface exhibit higher-than-bulk in-plane and negative out-of-plane dielectric constants, the latter corresponding physically to overscreening of the applied field. If we account for the gap region at the interface, the average out-of-plane dielectric constant is quite low (ε_{⊥}≈2), in agreement with reported measurements on confined thin films.

Stabilization of Dinuclear Rhodium and Iridium Clusters on Layered Titanate and Niobate Supports.

Atomically dispersed organometallic clusters can provide well-defined nuclearity of active sites for both fundamental studies as well as new regimes of activity and selectivity in chemical transformations. More recently, dinuclear clusters adsorbed onto solid surfaces have shown novel catalytic properties resulting from the synergistic effect of two metal centers to anchor different reactant species. Difficulty in synthesizing, stabilizing, and characterizing isolated atoms and clusters without agglomeration challenges allocating catalytic performance to atomic structure. Here, we explore the stability of dinuclear rhodium and iridium clusters adsorbed onto layered titanate and niobate supports using molecular precursors. Both systems maintain their nuclearity when characterized using aberration-corrected high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM). Statistical analysis of HAADF-STEM images revealed that rhodium and iridium dimers had mean cluster-to-cluster distances very similar to what is expected from a random distribution of atoms over a large area, indicating that they are dispersed without aggregation. The stability of dinuclear rhodium clusters supported on titanate nanosheets was also investigated by X-ray absorption fine structure (EXAFS), DRIFTS, and first-principles calculations. Both X-ray absorption spectroscopy and HAADF-STEM simulations, guided by density functional theory (DFT)-optimized structure models, suggested that rhodium dimers adsorb onto the nanosheets in an end-on binding mode that is stable up to 100 °C under reducing conditions. This study highlights that crystalline nanosheets derived from layered metal oxides can be used as model supports to selectively stabilize dinuclear clusters, which could have implications for heterogeneous catalysis.

Hydrogen Bond Thermodynamics in Aqueous Acid Solutions: A Combined DFT and Classical Force-Field Approach.

The thermodynamics of hydrogen bonds in aqueous and acidic solutions significantly impacts the kinetics and thermodynamics of acid reaction chemistry. We utilize in this work a multiscale approach, combining density functional theory (DFT) with classical molecular dynamics (MD) to model hydrogen bond thermodynamics in an acidic solution. Using thermodynamic cycles, we split the solution phase free energy into its gas phase counterpart plus solvation free energies. We validate this DFT/MD approach by calculating the aqueous phase hydrogen bond free energy between two water molecules (H2O-···-H2O), the free energy to transform an H3O+ cation into an H5O2+ cation, and the hydrogen bond free energy of protonated water clusters (H3O+-···-H2O and H5O2+-···-H2O). The computed equilibrium hydrogen bond free energy of H2O-···-H2O is remarkably accurate, especially considering the large individual contributions to the thermodynamic cycle. Turning to cations, we find the ion to be more stable than H3O+ by roughly 1-2 kBT. This small free energy difference allows for thermal fluctuation between the two idealized motifs, consistent with spectroscopic and simulation studies. Lastly, hydrogen bonding free energies between either H+ cation and H2O in solution were found to be stronger than between two H2O, though much less so than in vacuum because of dielectric screening in solution. Altogether, our results suggest the DFT/MD approach is promising for application in modeling hydrogen bonding and proton transfer thermodynamics in condensed phases.

Accelerated Anti-Markovnikov Alkene Hydrosilylation with Humic-Acid-Supported Electron-Deficient Platinum Single Atoms.

The hydrosilylation reaction is one of the largest-scale applications of homogeneous catalysis, and Pt homogeneous catalysts have been widely used in this reaction for the commercial manufacture of silicon products. However, homogeneous Pt catalysts result in considerable problems, such as undesired side reactions, unacceptable catalyst residues and disposable platinum consumption. Here, we synthesized electron-deficient Pt single atoms supported on humic matter (Pt1 @AHA_U_400), and the catalyst was used in hydrosilylation reactions, which showed super activity (turnover frequency as high as 3.0×107  h-1 ) and selectivity (>99 %). Density functional theory calculations reveal that the high performance of the catalyst results from the atomic dispersion of Pt and the electron deficiency of the Pt1 atoms, which is different from conventional Pt nanoscale catalysts. Excellent performance is maintained during recycle experiments, indicating the high stability of the catalyst.

Control of Molecular Bonding Strength on Metal Catalysts with Organic Monolayers for CO2 Reduction.

The development of separate levers for controlling the bonding strength of different reactive species on catalyst surfaces is challenging but essential for the design of highly active and selective catalysts. For example, during CO2 reduction, production of CO often requires balancing a trade-off between the adsorption strength of the reactant and product states: weak binding of CO is desirable from a selectivity perspective, but weak binding of CO2 leads to low activity. Here, we demonstrate a new method of controlling both CO2 adsorption and CO desorption over supported metal catalysts by employing a single self-assembly step where organic monolayer films were deposited on the catalyst support. Binding of phosphonic acid monolayers on supported Pt and Pd catalysts weakened CO binding via a through-support effect. The weakened CO adsorption was generally accompanied by decreased adsorption and reactivity of CO2. However, by the incorporation of basic amine functions at controlled positions in the modifying film, strong CO2 adsorption and hydrogenation reactivity could be restored. Thus, both through-surface and through-space interactions could be manipulated by design of the organic modifiers. After surface modification, the catalysts exhibited significantly improved selectivity (up to ∼99% at conversions near 50%) and activity toward CO production. Moreover, the rate of deactivation was notably reduced due to prevention of CO poisoning.

Tight binding models accurately predict band structures for copolymer semiconductors.

Conjugated polymers possess a wide range of desirable properties including accessible band gaps, plasticity, tunability, mechanical flexibility and synthetic versatility, making them attractive for use as active materials in organic photovoltaics (OPVs). In particular, push-pull copolymers, consisting of alternating electron-rich and electron-deficient moieties, offer broad optical absorption, tunable band gaps, and increased charge transfer between monomer units. However, the large number of possible monomer combinations to explore means screening OPV copolymers by first-principles quantum calculations is computationally intensive. If copolymer band structures could be rapidly computed from homopolymer data, potential materials could be screened more efficiently. In this work, we construct tight binding models of copolymer band structures with parameters determined by density functional theory (DFT) calculations on homopolymers. We use these models to predict copolymer valence and conduction bands, which compare well to direct DFT calculations of copolymer band structures.

Polaron hopping barriers and rates in semiconducting polymers.

Conjugated polymers are potential next-generation materials for organic electronic devices. The ability of these materials to transport charges is a key factor limiting their performance. Charge carriers in conjugated polymers are localized by disorder and polaronic effects. Charge transport in these materials is often described by thermally activated hopping, with a rate given by Marcus theory. The polaron hopping activation energy determines the temperature dependence of the Marcus rate. This energy barrier is dictated by the transition state, in which the charge carrier is equally divided between the initial and final locations. The prefactor for the polaron hopping rate is set by the charge tunneling rate between the initial and final locations. We use a tight-binding polaron model, in which charge carriers are stabilized by both nuclear reorganization and polarization of the surrounding dielectric, to compute the activation energy, charge tunneling rate and overall rate constant for intrachain and interchain charge hopping processes in poly(3-hexylthiophene) (P3HT) crystalline lamellae and amorphous melts. Charge transport in these two environments is limited by interchain hopping processes. Both hopping barriers and rates predicted by the model are in good agreement with experiments on a variety of crystalline and amorphous P3HT materials. Qualitatively, the barriers largely depend on how well the transition state is stabilized by polarization effects, and on the hopping integral between the initial and final locations, both of which penalize hopping over longer distances.

Thermal Fluctuations Lead to Cumulative Disorder and Enhance Charge Transport in Conjugated Polymers.

All conjugated polymers examined to date exhibit significant cumulative lattice disorder, although the origin of this disorder remains unclear. Using atomistic molecular dynamics (MD) simulations, the detailed structures for single crystals of a commonly studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) are obtained. It is shown that thermal fluctuations of thiophene rings lead to cumulative disorder of the lattice with an effective paracrystallinity of about 0.05 in the π-π stacking direction. The thermal-fluctuation-induced lattice disorder can in turn limit the apparent coherence length that can be observed in diffraction experiments. Calculating mobilities from simulated crystal structures demonstrates that thermal-fluctuation-induced lattice disorder even enhances charge transport in P3HT. The mean inter-chain charge transfer integral is enhanced with increasing cumulative lattice disorder, which in turn leads to pathways for fast charge transport through crystals.

Energetics of exciton binding and dissociation in polythiophenes: a tight binding approach.

Organic photovoltaics offer a potential low-cost alternative to inorganic solar cells. Crucial to the performance of these devices is the generation of free charges, which occurs through the dissociation of excitons. Here we study excitons in polythiophenes, their stability and energetics of dissociation and separation into charge carriers. Excitons are excited electron and hole pairs bound by Coulomb interactions. To separate into unbound charges, the exciton binding energy must be overcome. We use a tight binding Hamiltonian to describe the exciton binding energy and its dissociation potential, for an exciton confined to a single polymer chain. Our model accounts for polaronic effects, arising from reorganization of nuclei and from polarization of the surrounding dielectric, which stabilize the separated carriers and thereby affect the exciton dissociation potential. We examine the effects of an applied electric field on the dissociation potential, and relate the field strength necessary to unbind the hole-electron pair to the maximum attractive Coulomb force between them. We apply our model to study the exciton at a donor-acceptor interface on a block-copolymer. Interfacial polarization alters the exciton binding potential, rendering the hole-electron pair easier to unbind.

Elementary kinetics of nitrogen electroreduction on Fe surfaces.

Electrochemical ammonia synthesis could provide a sustainable and efficient alternative to the energy intensive Haber-Bosch process. Development of an active and selective N2 electroreduction catalyst requires mechanism determination to aid in connecting the catalyst composition and structure to performance. Density functional theory (DFT) calculations are used to examine the elementary step energetics of associative N2 reduction mechanisms on two low index Fe surfaces. Interfacial water molecules in the Heyrovsky-like mechanism help lower some of the elementary activation barriers. Electrode potential dependent barriers show that cathodic potentials below -1.5 V-RHE (reversible hydrogen electrode) are necessary to give a significant rate of N2 electroreduction. DFT barriers suggest a larger overpotential than expected based on elementary reaction free energies. Linear Brønsted-Evans-Polanyi relationships do not hold across N-H formation steps on these surfaces, further confirming that explicit barriers should be considered in DFT studies of the nitrogen reduction reaction.

BCC-Phased PdCu Alloy as a Highly Active Electrocatalyst for Hydrogen Oxidation in Alkaline Electrolytes.

Anion-exchange membrane fuel cells hold promise to greatly reduce cost by employing nonprecious metal cathode catalysts. More efficient anode catalysts are needed, however, to improve the sluggish hydrogen oxidation reaction in alkaline electrolytes. We report that BCC-phased PdCu alloy nanoparticles, synthesized via a wet-chemistry method with a critical thermal treatment, exhibit up to 20-fold HOR improvement in both mass and specific activities, compared with the FCC-phased PdCu counterparts. HOR activity of the BCC-phased PdCu is 4 times or 2 times that of Pd/C or Pt/C, respectively, in the same alkaline electrolyte. In situ HE-XRD measurements reveal that the transformation of PdCu crystalline structure favors, at low annealing temperature (<300 °C), the formation of FCC structure. At higher annealing temperatures (300-500 °C), a BCC structure dominates the PdCu NPs. Density functional theory (DFT) computations unravel a similar H binding strength and a much stronger OH binding of the PdCu BCC surface (cf. FCC surface), both of which are simultaneously close to those of Pt surfaces. The synergistic optimization of both H and OH binding strengths is responsible for the enhancement of HOR activity on BCC-phased PdCu, which could serve as an efficient anode catalyst for anion-exchange membrane fuel cells. This work might open a new route to develop efficient HOR catalysts from the perspective of crystalline structure transformation.

Thermodynamics of Alkanethiol Self-Assembled Monolayer Assembly on Pd Surfaces.

We investigate the structure and binding energy of alkanethiolate self-assembled monolayers (SAMs) on Pd (111), Pd (100), and Pd (110) facets at different coverages. Dispersion-corrected density functional theory calculations are used to correlate the binding energy of alkanethiolates with alkyl chain length and coverage. The equilibrium coverage of thiolate layers strongly prefers 1/3 monolayer (ML) on the Pd (111) surface. The coverage of thiolates varies with chemical potential on Pd (100) and Pd (110), increasing from 1/3 to 1/2 ML on (100) and from 1/4 to 1/2 ML on (110) as the thiol chemical potential is increased. Higher coverages are driven by attractive dispersion interactions between the extended alkyl chains, such that transitions to higher coverages occur at lower thiol chemical potentials for longer chain thiolates. Stronger adsorption to the Pd (100) surface causes the equilibrium Wulff construction of Pd particles to take on a cubic shape upon saturation with thiols. The binding of H, O, and CO adsorbates is weakened as the thiolate coverage is increased, with saturation coverages causing unfavorable binding of O and CO on Pd (100) and weakened binding on other facets. Temperature-dependent CO diffuse reflectance infrared Fourier transform spectroscopy experiments are used to corroborate the weakened binding of CO in the presence of thiolate SAMs of varying surface density. Preliminary results of multiscale modeling efforts on the Pd-thiol system using a reactive force field, ReaxFF, are also discussed.

A first-principles study of stability of surface confined mixed metal oxides with corundum structure (Fe2O3, Cr2O3, V2O3).

Surface-confined mixed metal oxides can have different chemical properties compared to their host metal oxide support. For this reason, mixed transition metal oxides can offer tunable redox properties. Herein, we use density functional theory to predict the stability of the (0001) surface termination for mixed metal oxides consisting of Fe2O3, Cr2O3 and V2O3. We show that the pure oxide surface stability can predict the surface segregation preference of the surface-confined mixed metal oxides. We focus on substitution of Fe in the V2O3(0001) surface, for which we observe that Fe substitution increases the reducibility of the resulting mixed metal oxide surface. Our results suggest Fe is only stable on the surface under very high temperature and/or low-pressure conditions. Using thermodynamic relationships, we predict the transition points for these surface-confined mixed metal oxides at which exchange between surface/subsurface and subsurface/surface metal atoms occur due to changes in the oxygen chemical potential.

Applied Potentials in Variable-Charge Reactive Force Fields for Electrochemical Systems.

An atomic description of water dynamics and electrochemical properties at electrode-electrolyte interfaces is presented using molecular dynamics with the third generation of the charge-optimized many-body (COMB3) potential framework. Externally applied potentials in electrochemical applications were simulated by offsetting electronegativity on electrode atoms. This approach is incorporated into the variable charge scheme within COMB3 and is used to investigate electrochemical systems consisting of two Cu electrodes and a water electrolyte with varying concentrations of hydroxyls (OH-) and protons (H+). The interactions between the electronegativity offset method and the charge equilibration method in a variable charge scheme are analyzed. In addition, a charge equilibration method for electrochemical applications is proposed, where the externally applied potentials are treated by the electronegativity offset on the electrodes thus enforcing charge neutrality on the electrolyte. This method is able to qualitatively capture the relevant electrochemistry and predict consistently correct voltages with precalibration.

First-Principles Calculation of Pt Surface Energies in an Electrochemical Environment: Thermodynamic Driving Forces for Surface Faceting and Nanoparticle Reconstruction.

Platinum is a widely used catalyst in aqueous and electrochemical environments. The size and shape of Pt nanoparticles and the faceting (and roughness) of extended Pt surfaces change during use in these environments due to dissolution, growth, and reconstruction. Further, many Pt nanoparticle synthesis techniques are carried out in an aqueous environment. The surface structures formed are impacted by the relative surface energies of the low index facets in these environments. Density functional theory is used to calculate the surface energy of the low index facets of platinum as a function of electrochemical potential and coverage of adsorbed hydrogen, hydroxide, oxygen, and the formation of surface oxide in an aqueous environment. Whereas Pt(111) is the lowest energy bare surface in vacuum, the strong adsorption of hydrogen to Pt(100) at low potentials and of hydroxide to Pt(110) and oxygen to Pt(100) at high potentials drives these surfaces to be more stable in an electrochemical environment. We experimentally conditioned a polycrystalline platinum electrode by cycling the potential and find a growth in the total area as well as in the fraction of 110 and 100 sites, which are lower in energy at potentials where dissolved Pt is deposited or surface oxide is reduced. Further, we find that the lower surface energy of Pt(100) at low potentials may play a role in the growth of tetrahexahedral nanoparticles seen on square wave cycling of spherical Pt nanoparticles. Wulff constructions are presented as a function of Pt electrode potential.

Polaron formation mechanisms in conjugated polymers.

In semiconducting polymers, interactions with conformational degrees of freedom can localize charge carriers, and strongly affect charge transport. Polarons can form when charges induce deformations of the surrounding medium, including local vibrational modes or dielectric polarization. These deformations then interact attractively with the charge, tending to localize it. First we investigate vibrational polaron formation in poly(3-hexylthiophene) [P3HT], with a tight-binding model for charges hopping between adjacent rings, coupled to ring distortions. We use density functional theory (DFT) calculations to determine coupling constants, including the "spring constant" for ring distortions and the coupling to the charge carrier. On single chains, we find only broad, weakly bound polarons by this mechanism. In 2d crystalline layers of P3HT, even weak transverse hopping between chains destabilizes this polaron. Then, we consider polarons stabilized by dielectric polarization, described semiclassically with a polarizable continuum interacting with the carrier wavefunction. In contrast to vibrational polarons, we find dielectrically stabilized polarons in P3HT are narrower, more strongly bound, and stable in 2d crystalline layers.

Identification of Second Shell Coordination in Transition Metal Species Using Theoretical XANES: Example of Ti-O-(C, Si, Ge) Complexes.

X-ray absorption near-edge structure (XANES) is a common technique for elucidating oxidation state and first shell coordination geometry in transition metal complexes, among many other materials. However, the structural information obtained from XANES is often limited to the first coordination sphere. In this study, we show how XANES can be used to differentiate between C, Si, and Ge in the second coordination shell of Ti-O-(C, Si, Ge) molecular complexes based on differences in their Ti K-edge XANES spectra. Experimental spectra were compared with theoretical spectra calculated using density functional theory structural optimization and ab initio XANES calculations. The unique features for second shell C, Si, and Ge present in the Ti K pre-edge XANES are attributed to the interaction between the Ti center and the O-X (X = C, Si, or Ge) antibonding orbitals.

Co-adsorption of Cations as the Cause of the Apparent pH Dependence of Hydrogen Adsorption on a Stepped Platinum Single-Crystal Electrode.

The successful deployment of advanced energy-conversion systems depends critically on our understanding of the fundamental interactions of the key adsorbed intermediates (hydrogen *H and hydroxyl *OH) at electrified metal-aqueous electrolyte interfaces. The effect of alkali metal cations (Li+ , Na+ , K+ , Cs+ ) on the non-Nernstian pH shift of the step-related voltammetric peak of the Pt(553) electrode is investigated over a wide pH window (1 to 13) by means of experimental and computational methods. The co-adsorbed alkali cations along the step weaken the OH adsorption at the step sites, causing a positive shift of the potential of the step-related peak on Pt(553). Density functional calculations explain the observations on the identity and concentration of alkali cations on the non-Nernstian pH shift, and demonstrate that cation-hydroxyl co-adsorption causes the apparent pH dependence of "hydrogen" adsorption in the step sites of platinum electrodes.