Search NASA⌕ Search

SEARCH · Search NASA

Results for “electron and proton transfer processes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

Diffusion of Brønsted acidic dopants in conjugated polymers

Many semiconductor devices (e.g., light emitting diodes and photovoltaics) utilize heterojunctions of doped and undoped layers or depend on gradients of electronic doping to control charge transport. Understanding of the formation and stability of gradients in doping requires an understanding of diffusion of dopants and the complex changes in polymer properties that arise during doping. Conjugated polymers can be electrically doped by strong acids, but the details of the reaction mechanism and subsequent stability are not understood. Here, we show a clear kinetic isotope effect in the doping of thin films of poly(3-hexylthiophene) (P3HT) by bis(trifluoromethane)sulfonimide (HTFSI) from solution indicating that this doping process is limited by proton transfer to the polymer. Complementary X-ray photoelectron spectroscopy and dynamic secondary ion mass spectrometry (DSIMS) depth profiling of dopant concentrations show definitive evidence of dopant enrichment at the P3HT surface. These surface-limited concentration profiles suggest that diffusivity of dopants vary inversely with dopant concentration due to doping-induced changes to the structure of the conjugated polymer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding the Competition between Alcohol Formation and Dimerization during Electrochemical Reduction of Aromatic Carbonyl Compounds

The electrochemical reductive dimerization of small aromatic carbonyl compounds derived from lignocellulosic biomass is a crucial C−C coupling reaction for upgrading small molecules to longchain hydrocarbons, particularly in the synthesis of drop-in sustainable aviation fuels. Although other electrochemical reduction reactions of these reactants (i.e., hydrogenation and hydrogenolysis) have undergone extensive mechanistic investigation, the understanding of dimerization remains relatively underdeveloped. Most importantly, there is a lack of understanding of the selectivity-determining step between dimerization and monomer reduction and critical factors that can affect this step. In this study, we provide a comprehensive mechanistic model to explain the competition between dimerization and monomer reduction of benzaldehyde under various conditions. Our model proposes that the selectivity between dimerization and monomer reduction depends on the competition between desorption of a ketyl radical from the electrode, necessary for dimerization, and further reduction of the ketyl radical to an alcohol on the electrode by proton-coupled electron transfer (PCET). Computationally comparing the adsorption/desorption energy and PCET activation barrier energy is challenging because conventional DFT calculations substantially underestimate the PCET kinetic barriers. In this study, we employed constrained DFTbased configuration interaction (CDFT-CI) to obtain a reliable comparison of these energies. Our mechanistic model was tested and supported by experimental results obtained with four electrodes (Cu, Pb, Bi, graphite), three pH conditions (acidic, neutral, basic), and three potentials. Our study offers a coherent mechanistic foundation that can explain how each of these conditions impacts the desorption and PCET processes and the selectivities for dimerization and alcohol production.

09 BIOMASS FUELS↗

Mechanism of Catechol Oxidation by the Coupled Bicupric Active Site of Tyrosinase: Completion of the Oxygenase/Oxidase Reaction Cycle

Tyrosinase contains a coupled binuclear copper (CBC) active site, which in its bicuprous form (deoxy-Ty) binds O 2 to form a side-on peroxide [Cu(II) 2 O 2 ] 2+ intermediate (oxy-Ty) that performs the regioselective monooxygenation of monophenols to catechols and their subsequent 2e − oxidation to quinones. Previously, we used spectroscopic, kinetic, and computational methods to elucidate the mechanism of the initial steps in the monooxygenation reaction. Herein, we investigated the final step in catecholate formation and provide experimental and computational results elucidating the catechol oxidation reaction, formally a proton-coupled 2e − process. We employed single-turnover stoppedflow absorption to observe the elusive catecholate-Cu(II) 2 OH intermediate via the anaerobic reaction of resting Cu(II) 2 OH Ty (met-Ty) with the catecholate substrate. This intermediate was cryo-trapped and characterized by electron paramagnetic resonance and X-ray absorption spectroscopies. These experimental results were correlated to quantum mechanics/molecular mechanics (QM/MM) and QM calculations to describe the complete catalytic cycle of tyrosinase, revealing that the 2e − oxidation of catechol to quinone involves two steps: (i) coupled 1H + /1e − transfer from the bound monoanionic catecholate to the bridged hydroxide of the met-Ty active site, followed by (ii) the second 1e − transfer. The latter is calculated to be the rate-limiting step in catechol oxidation, confirmed by experimental solvent kinetic isotope effect studies. Our computational data suggest that quinone release from the protein provides the necessary driving force for this reaction. These results establish a detailed description of this oxidase cycle of the tyrosinase mechanism and more broadly provide molecular-level insights into the diverse reactivity of CBC sites in biology.

QM/MM modeling↗

Investigation of the background in coherent J/ψ production at the EIC

Understanding various fundamental properties of nucleons and nuclei is among the most important scientific goals at the upcoming Electron-Ion Collider (EIC). With the unprecedented opportunity provided by the next-generation machine, the EIC might provide definitive answers to many standing puzzles and open questions in modern nuclear physics. We investigate one of the golden measurements proposed at the EIC, which is to obtain the spatial gluon density distribution within a lead (Pb) nucleus. The proposed experimental process is the exclusive J/ψ vector-meson production off the Pb nucleus: e + Pb → e' + J/ψ + Pb'. The Fourier transformation of the momentum transfer |t| distribution of the coherent diffraction is the transverse gluon spatial distribution. In order to measure it, the experiment has to overcome an overwhelmingly large background arising from the incoherent diffractive production, where the nucleus Pb' mostly breaks up into fragments of particles in the far-forward direction close to the hadron-going beam rapidity. We systematically study the rejection of incoherent J/ψ production by vetoing products from these nuclear breakups—protons, neutrons, and photons—which is based on the BeAGLE event generator and the most up-to-date EIC Far-forward Interaction Region design. The achieved vetoing efficiency, the ratio between the numbers of vetoed events and total incoherent events, ranges from about 80% to 99% depending on |t|. Assuming a 5% smearing applied to the reconstructed |t| resolution in the Sar t re model, this vetoing efficiency can suppress the incoherent background to at least the first minimum of the coherent |t| distribution. Experimental and accelerator machine challenges as well as potential improvements are discussed herein.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Deciphering Reaction Mechanisms of Molecular Proton Reduction Catalysts with Cyclic Voltammetry: Kinetic vs Thermodynamic Control

The kinetics and thermodynamics of elementary reaction steps involved in the catalytic reduction of protons to hydrogen define the reaction landscape for catalysis. The mechanisms can differ in the order of the elementary proton transfer, electron transfer, and bond-forming steps and can be further differentiated by the sites at which protons and electrons localize. Access to fully elucidated mechanistic, kinetic, and thermochemical details of molecular catalysts is crucial to facilitate the development of new catalysts that operate with optimal efficiency, selectivity, and durability. The mechanism by which a catalyst operates, as well as the kinetics and thermodynamics associated with the individual steps, can often be accessed through electroanalytical studies. Here, this Account details the application of cyclic voltammetry to interrogate reaction mechanisms and quantify the kinetics and thermodynamics of elementary reaction steps for a series of molecular catalysts that mediate electrochemical proton reduction. I distinguish the limiting scenarios wherein a catalyst operates under kinetic control vs thermodynamic control, with a focus on detecting how cyclic voltammetry features shift with proton source strength and concentration, as well as scan rate. For systems that operate under kinetic control, catalytic currents are observed at, or slightly positive toward, the formal potential for the redox process that triggers catalysis. Under thermodynamic control, catalytic responses shift as a function of the proton source pKa and effective pH of the solution. After drawing this distinction, we introduce the appropriate voltammetry experiments and accompanying analytical expressions for extracting key metrics from the data. To illustrate analytical strategies to quantify elementary reaction steps of catalysts operating under kinetic control, I describe our studies of proton reduction catalysts Co(dmgBF 2 ) 2 (CH 3 CN) 2 (dmgBF 2 = difluoroboryl-dimethylglyoxime) and [Ni(P 2 Ph N 2 Ph ) 2 ] 2+ (P 2 Ph N 2 Ph = 1,5-phenyl-3,7-phenyl-1,5-diaza-3,7-diphosphacyclooctane). Here, peak shift analysis, foot-of-the-wave analysis, and plateau current analysis are applied to data sets wherein voltammetric response are recorded as a function of catalyst concentration, proton source concentration, proton source strength, and scan rate to quantify rate constants for elementary proton transfer and bond-forming steps in a catalytic cycle. Further, the case study of [Ni(P 2 Ph N 2 Ph ) 2 ] 2+ illustrates how complementary spectroscopic methods can bolster the mechanistic assignment. Collectively, these two studies showcase how detailed mechanistic studies inform on rate-limiting elementary steps in catalysis and other key processes underpinning catalysis. Second, I present analytical strategies to interrogate catalysts operating under thermodynamic control, centered on the case study of [Ni II (P 2 Ph N 2 Bn ) 2 ] 2+ (P 2 Ph N 2 Bn = 1,5-dibenzyl-3,7-diphenyl-1,5-diaza-3,7-diphosphacyclooctane). Here, the application of nonaqueous Pourbaix theory to extract thermodynamic information is introduced, and the construction of a coupled Pourbaix diagram is detailed. This study identifies ligand-based protonation as the key process that places catalysis under thermodynamic control and influences the reaction mechanism. Together, the work detailed in this Account showcases the utility of electroanalytical methods to disentangle complex reaction mechanisms and extract key thermochemical and kinetic parameters for elementary steps of catalysis. Through detailed presentation of the key analytical expressions that underpin these analyses, this Account seeks to facilitate the adoption of cyclic voltammetry by the community to fully extract kinetic, thermochemical, and mechanistic information on electrochemical small-molecule activation.

catalysts↗

Hydrolysis of Small Oxo/Hydroxo Molecules Containing High Oxidation State Actinides (Th, Pa, U, Np, Pu): A Computational Study

The energetics of hydrolysis reactions for high oxidation states of oxo/hydroxo monomeric actinide species (Th IV O 2 , Pa IV O 2 , U IV O 2 , Pa V O 2 (OH), U V O 2 (OH), U VI O 3 , Np VI O 3 , Np VII O 3 (OH), and Pu VII O 3 (OH)) were calculated at the CCSD(T) level. The first step is the formation of a Lewis acid/base adduct with H 2 O (hydration), followed by a proton transfer to form a dihydroxide molecule (hydrolysis); this process is repeated until all oxo groups are hydrolyzed. The physisorption (hydration) for each H 2 O addition was predicted to be exothermic, ca. -20 kcal/mol. The hydrolysis products are preferred energetically over the hydration products for the +IV and +V oxidation states. The compounds with An VI are a turning point in terms of favoring hydration over hydrolysis. For An VII O 3 (OH), hydration products are preferred, and only two waters can bind; the complete hydrolysis process is now endothermic, and the oxidation state for the An in An(OH) 7 is +VI with two OH groups each having one-half an electron. The natural bond order charges and the reaction energies provide insights into the nature of the hydrolysis/hydration processes. The actinide charges and bond ionicity generally decrease across the period. Here, the ionic character decreases as the oxidation state and coordination number increase so that covalency increases moving to the right in the actinide period.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optical Absorption Spectroscopy Probes Water Wire and Its Ordering in a Hydrogen-Bond Network

Water wires, quasi-one-dimensional chains composed of hydrogen-bonded (H-bonded) water molecules, play a fundamental role in numerous chemical, physical, and physiological processes. Yet direct experimental detection of water wires has been elusive so far. Based on advanced many-body theory that includes electron-hole interactions, we report that optical absorption spectroscopy can serve as a sensitive probe of water wires and their ordering. In both liquid and solid water, the main peak of the spectrum is discovered to be a charge-transfer exciton. In water, the charge-transfer exciton is strongly coupled to the H-bonding environment where the exciton is excited between H-bonded water molecules with a large spectral intensity. In regular ice, the spectral weight of the charge-transfer exciton is enhanced by a collective excitation occurring on proton-ordered water wires, whose spectral intensity scales with the ordering length of water wire. The spectral intensity and excitonic interaction strength reaches its maximum in ice XI, where the long-range ordering length yields the most pronounced spectral signal. Our findings suggest that water wires, which widely exist in important physiological and biological systems and other phases of ice, can be directly probed by this approach. Published by the American Physical Society 2025

Tang, Fujie (ORCID:0000000197615359)↗

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗

Nature of the Reactive Biferric Peroxy Intermediate P′ in the Arylamine Oxygenases and Related Binuclear Fe Enzymes

Binuclear nonheme iron enzymes activate O 2 to perform a wide range of chemical transformations. The process of O 2 activation typically involves a biferric peroxy-level intermediate P. It has been previously found that this intermediate undergoes further activation, either protonation or rearrangement to form P′ or further oxidation to form high-valent intermediates Q or X. Here, this study defines the structure of the P′ intermediate in the N-oxygenases CmlI (and AurF based on previous data) using nuclear resonance vibrational spectroscopy (NRVS) in conjugation with density functional theory (DFT) calculations. These results, combined with variable temperature variable field (VTVH) magnetic circular dichroism (MCD) spectroscopy on the 1-electron cryoreduced P′, define the structure of the P′ intermediate as a μ-1,2-hydroxoperoxo biferric site with a second hydroxide bridge. Reaction coordinate calculations demonstrate that single electron transfer (SET) is facilitated by protonation of the peroxo, activating its reductive cleavage, and that the additional hydroxide bridge does not impact this reaction. VTVH MCD studies further reveal that the hydroxide bridge is absent in the biferrous site, suggesting that during the O 2 reaction with the biferrous site, a water molecule forms the hydroxide bridge in providing the proton that activates the peroxide in P′ for reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessment of Hydroxyl Radical Reactivity in Sulfur-Containing Amino Acid Models Under Acidic pH

Methionine residues in proteins and peptides are frequently oxidized by losing one electron. The presence of nearby amide groups is crucial for this process, enabling methionine to participate in long-range electron transfer. Hydroxyl radical (HO•) plays an important role being generated in aerobic organisms by cellular metabolisms as well as by exogenous sources such as ionizing radiations. The reaction of HO• with methionine mainly affords the one-electron oxidation of the thioether moiety through two consecutive steps (HO• addition to the sulfur followed by HO− elimination). We recently investigated the reaction of HO• with model peptides mimicking methionine and its cysteine-methylated counterpart, i.e., CH3C(O)NHCHXC(O)NHCH3, where X = CH2CH2SCH3 or CH2SCH3 at pH 7. The reaction mechanism varied depending on the distance between the sulfur atom and the peptide backbone, but, for a better understanding of various suggested equilibria, the analysis of the flux of protons is required. We extended the previous study to the present work at pH 4 using pulse radiolysis techniques with conductivity and optical detection of transient species, as well as analysis of final products by LC-MS and high-resolution MS/MS following γ-radiolysis. Comparing all the data provided a better understanding of how the presence of nearby amide groups influences the one-electron oxidation mechanism.

Biochemistry & Molecular Biology↗

Ammonia Synthesis under Ambient Conditions: Insights into Water–Nitrogen–Magnetite Interfaces

New routes for transforming nitrogen into ammonia at ambient conditions would be a milestone toward an energy efficient and economically attractive production route in comparison to the traditional Haber-Bosch process. Recently, the synthesis of ammonia from water and nitrogen at room temperature and atmospheric pressure has been reported to be catalyzed by Fe 3 O 4 at the air-water interface. By integrating ambient pressure X-ray photoelectron spectroscopy and ab initio molecular dynamics and free energy calculations, we investigate the underlying thermodynamic mechanisms governing ammonia and hydrazine formation at the water-Fe 3 O 4 -nanoparticle interface. Here, we find that, unlike pure Fe 3 O 4 where N 2 can only interact with a limited number of Fe sites, hydroxylated species introduce large and diverse adsorption geometries where N 2 can bind through either Fe sites or Fe-OH groups, each of which are capable of independently facilitating proton-coupled electron transfer.

Chandy, Sruthy K. [University of California, Berke↗

Photon-odderon interference in exclusive χ c charmonium production at the Electron-Ion Collider

Exclusive C = + 1 scalar, axial-vector, and tensor quarkonium production in high-energy electron-proton scattering requires a C -odd t -channel exchange of a photon or a three-gluon ladder. We derive the expressions for the corresponding amplitudes. The relative phase of the photon vs three-gluon exchange amplitudes is determined by the sign of the light-front matrix element of the eikonal color current operator d a b c J + a J + b J + c at moderate x , and is not affected by small- x QCD evolution. Model calculations predict constructive interference, which is particularly strong for momentum transfer | t | ∼ 1 GeV 2 where the cross section for χ c J production exceeds that for pure photon exchange by up to a factor of 4. Exclusive χ c J electroproduction at the high luminosity Electron-Ion Collider occurs with well measurable rates and measurements of these processes should find evidence for the perturbative odderon exchange. We also compute the total electroproduction cross section as a function of energy and provide first estimates of the number of χ c J events per month at the Electron-Ion Collider design luminosity. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Productivity and bioproduct formation in phototropic knock/out mutants in micro algae

Phototropin is a blue light receptor, which mediates a variety of blue-light elicited physiological processes in plants and algae. In higher plants these processes include phototropism, chloroplast movement and stomatal opening. In the green alga Chlamydomonas reinhardtii, phototropin plays a vital role in progression of the sexual life cycle and in the control of the eye spot size and light sensitivity. Phototropin is also involved in blue-light mediated changes in the synthesis of chlorophylls, carotenoids, chlorophyll binding proteins. We compared the transcriptome of phototropin knock out (PHOT KO) mutant and wild-type parent to analyze differences in gene expression in high light grown cultures (500 μmol photons m⁻²s⁻¹). Our results indicate the up-regulation of genes involved in photosynthetic electron transport chain, carbon fixation pathway, starch, lipid, and cell cycle control genes. With respect to photosynthetic electron transport genes, genes encoding proteins of the cytochrome b6f and ATP synthase complex were up regulated potentially facilitating proton-coupled electron transfer. In addition genes involved in limiting steps in the Calvin cycle Ribulose-1,5-biphosphate carboxylase/oxygenase (RuBisCO), Sidoheptulose 1,7 biophosphatase (SBPase), Glyceraldehyde-3-phosphate dehydrogenase (3PGDH) and that mediate cell-cycle control (CDK) were also up regulated along with starch synthase and fatty acid biosynthesis genes involved in starch and lipid synthesis. In addition, transmission electron micrographs show increased accumulation of starch granules in PHOT mutant compared to wild type, which is consistent with the higher expression of starch synthase genes. Collectively, the altered patterns of gene expression in the PHOT mutants were associated with a two-fold increase in growth and biomass accumulation compared to wild type when grown in environmental photobioreactors (Phenometrics) that simulate a pond environment. In conclusion, our studies suggest that phototropin may be a master gene regulator that suppresses rapid cell growth and promotes gametogenesis and sexual recombination in wild type strains.

59 BASIC BIOLOGICAL SCIENCES↗

Light‐driven catalysis with engineered enzymes and biomimetic systems

Abstract Efforts to drive catalytic reactions with light, inspired by natural processes like photosynthesis, have a long history and have seen significant recent growth. Successfully engineering systems using biomolecular and bioinspired catalysts to carry out light‐driven chemical reactions capitalizes on advantages offered from the fields of biocatalysis and photocatalysis. In particular, driving reactions under mild conditions and in water, in which enzymes are operative, using sunlight as a renewable energy source yield environmentally friendly systems. Furthermore, using enzymes and bioinspired systems can take advantage of the high efficiency and specificity of biocatalysts. There are many challenges to overcome to fully capitalize on the potential of light‐driven biocatalysis. In this mini‐review, we discuss examples of enzymes and engineered biomolecular catalysts that are activated via electron transfer from a photosensitizer in a photocatalytic system. We place an emphasis on selected forefront chemical reactions of high interest, including CH oxidation, proton reduction, water oxidation, CO 2 reduction, and N 2 reduction.

Edwards, Emily H.↗

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗

Gluon Sivers function from forward exclusive χ c 1 photoproduction on unpolarized protons

Exclusive production of a χ c 1 axial vector quarkonia in photon-proton scattering at high energies requires a C -odd t -channel exchange. In the limit of vanishing momentum transfer this occurs either via the exchange of a photon, the Primakoff process, where the spin of the proton does not change. For axial-vector meson production, as a consequence of the Landau-Yang theorem, the Primakoff cross section is finite as t → 0 . Alternatively, a C -odd spin dependent Odderon can be exchanged, which involves a spin flip of the proton. The resulting cross section is related to the square of the collinear trigluon correlator or the k ⊥ -moment of the gluon Sivers function. Using two models for the gluon Sivers function from the literature we compute the ratio of Sivers to Primakoff cross sections and the angular coefficient λ θ governing the angular distribution of the χ c 1 → J / ψ + γ decay as functions of x . We point out that these observables constrain the magnitude of the gluon Sivers function at small x which could be accessed in electron-proton scattering and ultraperipheral proton-proton and nucleus-proton collisions. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electrocatalytic Reduction of CO 2 to CO over Ag(110) and Cu(211) Modeled by Grand-Canonical Density Functional Theory

We report the results of modeling CO 2 reduction (CO 2 R) to CO over Ag(110) and Cu(211) surfaces at different applied potentials using grand-canonical density functional theory (GC-DFT), a method specifically designed to accurately model electrochemical systems. In addition to demonstrating GC-DFT’s ability to accurately model electrochemical processes, we also compare it with the computational hydrogen electrode (CHE) approach. GC-DFT predicts that the geometries of these reacting systems strongly depend on the applied potential and the Helmholtz free energies vary nonlinearly with the applied potential, which contradicts a central assumption of the CHE approach. The CHE approach neglects the change in the number of electrons on the electrode surface at different applied potentials, which reduces its accuracy as the potential changes from the potential of zero charge. Our results further demonstrate that the grand free energies of the reaction intermediates not only depend on the value of the applied potential but also on the metal surface type, adsorption site, and adsorbate. GC-DFT’s ability to predict the effect of the applied potential on adsorbate geometry enables it to evaluate different possible reaction mechanisms at different applied potentials. For instance, GC-DFT predicts that the first step of CO2R likely switches from proton-coupled electron transfer to sequential electron transfer and then proton transfer at more reducing potentials, a result that cannot be determined by the CHE because it assumes that all electron transfers are coupled to proton transfers and neglects the effect of the applied potential on the adsorbate geometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

[Fe 4 S 4 ] cubane in sulfite reductases: new insights into bonding properties and reactivity

The dissimilatory sulfite reductase enzyme has very characteristic active site where the substrate binds to an iron site, ligated by a siroheme macrocycle and a thiol directly connected to a [Fe 4 S 4 ] cluster. This arrangement gives the enzyme remarkable efficiency in reducing sulfite and nitrite all the way to hydrogen sulfide and ammonia. For the first time we present a theoretical study where substrate binding modalities and activation are elucidated using active site models containing proton supply side chains and the [Fe 4 S 4 ] cluster. Density functional theory (DFT) was deployed in conjunction with the energy decomposition scheme (as implemented in AMS), the quantum theory of atoms in molecules (QTAIM), and conceptual DFT (cDFT) descriptors. We quantified the role of the electrostatic interactions inside the active site created by the side chains as well as the influence of the [Fe4S4] cluster on the substrate binding. Furthermore, using conceptual DFT results we shed light of the activation process, thus, laying foundation for further mechanistic studies. We found that the bonding of the ligands to the iron complex is dominated by electrostatic interactions, but the presence of the [Fe 4 S 4 ] cubane leads to substantial changes in electronic interaction. The spin state of the cubane, however, affects the binding energy only marginally. The conceptual DFT results show that the presence of the [Fe 4 S 4 ] cubane affects the reactivity of the active site as it is involved in electron transfer. This is corroborated by an increase in the electrophilicity index, thus making the active site more prone to react with the ligands. The interaction energies between the ligand and the siroheme group are also increased upon the presence of the cubane group, thus, suggesting that the siroheme group is not an innocent spectator but plays an active role in the reactivity of the dSIR active site.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗