SEARCH · Search NASA
Results for “grand canonical potential kinetics”
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.
Laser-Irradiated Holey Graphene-Supported Single-Atom Catalyst towards Hydrogen Evolution and Oxygen Reduction
Single-atom catalysts (SAC) can boost the intrinsic catalytic activity of hydrogen evolution reaction (HER) and oxygen reduction reaction (ORR). However, the challenge remains due to the complex synthesis process and insufficient stability. A sustainable approach is applied to synthesizing SACs through laser irradiation and gaining mesoporous graphene oxide (MGO). The surface dangling bonds of nitrogen-doped MGO (NMGO) extract metal atoms species from Co or Fe metal foams and convert them to SAC via an appropriate synthesis approach. Notably, the Co-NMGO electrocatalyst requires low potentials of 146 mV to convey a current density of 10 mA cm -2 towards HER. Similarly, the Fe-NMGO electrocatalyst offers an onset of 0.79 V towards ORR in acidic solution. The individual metal atoms are confirmed via aberration-corrected scanning transmission electron microscopy, and X-ray absorption near-edge structure and extended X-ray absorption fine structure. Density functional theory calculations by applying the grand canonical potential kinetics model revealed that Co-NMGO shows the optimum free reaction energy of -0.17 eV at -0.1 V for HER, and Fe-NMGO has less limiting potential than that of Co-NMGO for ORR case. Furthermore, this work opens a new approach towards the synthesis of SAC and its mechanistic understandings.
Reaction mechanism and kinetics for N 2 reduction to ammonia on the Fe–Ru based dual-atom catalyst
Environmental and energy considerations demand that the Haber-Bosch process for reducing N 2 to NH 3 be replaced with electrochemical ammonia synthesis where the H atoms come from water instead of from H 2 . But a practical realization of electrochemical N 2 reduction reaction (NRR) requires the development of new generation electrocatalysts with low overpotential and high Faraday efficiency (FE). A major problem here is that the hydrogen evolution reaction (HER) competes with NRR. Herein, we consider new generation dual-site catalysts involving two different metals incorporated into a novel two-dimensional C 3 N–C 2 N heterostructure that provides a high concentration of well-defined but isolated active sites that bind two distinct metal atoms in a framework that facilitates electron transfer. We report here the mechanism and predicted kinetics as a function of applied potential for both NRR and HER for the (Fe–Ru)/C 3 N–C 2 N dual atom catalyst. These calculations employ the grand canonical potential kinetics (GCP-K) methodology to predict reaction free energies and reaction barriers as a function of applied potential. The rates are then used in a microkinetic model to predict the turn-over-frequencies (TOF) as a function of applied potential. At U = 0 V, the FE for NRR is 93%, but the current is only 2.0 mA cm –2 . The onset potential (at 10 mA cm –2 ) for ammonia on Fe–Ru/C 3 N–C 2 N is –0.22 V RHE . This leads to a calculated TOF of 434 h –1 per Fe–Ru site. In conclusion, we expect that the mechanisms for NRR and HER developed here will help lead to new generations of NRR with high TOF and FE.
Reaction mechanism and kinetics for CO 2 reduction on nickel single atom catalysts from quantum mechanics
Experiments have shown that graphene-supported Ni-single atom catalysts (Ni-SACs) provide a promising strategy for the electrochemical reduction of CO 2 to CO, but the nature of the Ni sites (Ni-N 2 C 2 , Ni-N 3 C 1 , Ni-N 4 ) in Ni-SACs has not been determined experimentally. Here, we apply the recently developed grand canonical potential kinetics (GCP-K) formulation of quantum mechanics to predict the kinetics as a function of applied potential (U) to determine faradic efficiency, turn over frequency, and Tafel slope for CO and H 2 production for all three sites. We predict an onset potential (at 10 mA cm –2 ) U onset = –0.84 V (vs. RHE) for Ni-N 2 C 2 site and U onset = –0.92 V for Ni-N 3 C 1 site in agreement with experiments, and U onset = –1.03 V for Ni-N 4 . We predict that the highest current is for Ni-N 4 , leading to 700 mA cm –2 at U = –1.12 V. To help determine the actual sites in the experiments, we predict the XPS binding energy shift and CO vibrational frequency for each site.
CO 2 reduction on pure Cu produces only H 2 after subsurface O is depleted: Theory and experiment
Significance Electrochemical reduction of CO 2 to fuels and valuable chemicals is a global imperative to address energy and climate challenges. Despite enormous efforts, there is not yet a commercial process. The best catalyst remains Cu metal, but a great deal of controversy revolves around the role of the presence of oxygen in the material. We combine theory and experiments to show that pure Cu without subsurface oxygen leads to no organic products. Thus, maintaining substantial oxygen at Cu surfaces is essential for long-term reduction of CO 2 to organics.
Metal dissolution from first principles: Potential-dependent kinetics and charge transfer
Metal dissolution and its inverse process are integral to both corrosion and electrodeposition; however, many mechanistic details regarding the dissolution process are challenging to decipher. These include how ion dissolution kinetics and charge transfer are influenced by the competition between metal and solvent interactions under an electrode potential. Here, in this work, we introduce a computational framework based on density functional theory with grand-canonical treatment of electrons to directly predict the potential energy landscape for metal dissolution at a constant potential. Using aluminum as an example, we demonstrate that dissolution kinetics is governed by competing kinetics between two physical processes associated with metal–metal bond breaking and ion-migration within the electrical double layer, respectively. We identify a kinetic transition between regimes dominated by each of these processes and show that this transition depends on the operating electrode potential, among other key factors. It is further found that kinetics and thermodynamics of these processes can be described with a simple, one-parameter Marcus-theory-type model. Beyond offering new understanding of charge transfer during dissolution, our simulation protocol provides a recipe for directly predicting other important quantities in electrochemical reactions from first principles that are difficult to measure, such as the symmetry factor.
Correlating surface adsorbate configuration and electrochemical performance of IrO 2 during seawater-relevant electrolysis
Seawater electrolysis alleviates freshwater demand to produce clean hydrogen while eliminating the need for water purification steps. The anodic process, seawater oxidation, typically requires high overpotentials and yields low selectivity to oxygen via the oxygen evolution reaction (OER), primarily due to the competing chlorine evolution reaction (CER) and hypochlorite evolution reaction (HCER) in pH-neutral conditions. Here, combining in situ surface-enhanced Raman characterization, grand canonical density functional theory-based calculations, and kinetic Monte Carlo simulations, we report the evolution of surface adsorbate configurations driven by applied potential and pH during seawater-relevant OER over IrO 2 , a highly OER-active and chloride-corrosion-resistant catalyst. As a result, the chemical properties of active sites, and thereby the kinetics of OER and CER/HCER, are effectively tuned. However, it is revealed that there is no optimal combination of potential and pH to achieve both high activity and high selectivity for seawater-relevant OER. To address this limitation, we establish a correlation between activity/selectivity and surface adsorbate configurations, enabling the optimization of highly active and OER-selective IrO 2 -based catalysts in seawater-relevant oxidation by modulating the local adsorbate environment of active sites.
CO (2) Electroreduction on Borated Copper Surfaces: Boron Active Sites, Not Copper
Boron-doped copper has recently emerged as an active and stable catalyst for the electrochemical reduction of CO 2 to value-added C 2 products. Here, in this work, we develop a realistic model of CO electroreduction on surface borides of copper under operational conditions, taking into account the effects of electrode potential, electrolyte environment, and pH. We study the possible reconstruction of the electrocatalyst surface using grand canonical DFT and global optimization to obtain a potential-dependent grand canonical ensemble description of metastable, hydrogen-covered catalyst surfaces. Two key surface configurations, low H coverage (LC, −0.6 V SHE ) and high H coverage (HC, −0.8 V SHE ), dominate this ensemble, with the former being kinetically persistent and C 2 selective under strongly reducing conditions. Nonmetallic boron sites on the surface copper boride are found to bind CO more strongly than copper sites, and mechanistic investigation of CO electroreduction pathways presents a surprisingly unconventional case of boron-centered reactivity in contrast to typical copper-centered reactivity. Neighboring boron sites present along boron chains on the surface copper boride are found to facilitate C−C coupling, thereby driving the high C 2 selectivity of this electrocatalyst.
Oxygen evolution reaction over catalytic single-site Co in a well-defined brookite TiO 2 nanorod surface
Efficient electrocatalysts for the oxygen evolution reaction (OER) are paramount to the development of electrochemical devices for clean energy and fuel conversion. However, the structural complexity of heterogeneous electrocatalysts makes it a great challenge to elucidate the surface catalytic sites and OER mechanisms. Here, we report that catalytic single-site Co in a well-defined brookite TiO 2 nanorod (210) surface (Co-TiO 2 ) presents turnover frequencies that are among the highest for Co-based heterogeneous catalysts reported to date, reaching 6.6 ± 1.2 and 181.4 ± 28 s -1 at 300 and 400 mV overpotentials, respectively. Based on grand canonical quantum mechanics calculations and the single-site Co atomic structure validated by in situ and ex situ spectroscopic probes, we have established a full description of the catalytic reaction kinetics for Co-TiO 2 as a function of applied potential, revealing an adsorbate evolution mechanism for the OER. The computationally predicted Tafel slope and turnover frequencies exhibit exceedingly good agreement with experiment.
Nitrate-To-Ammonia Electroconversion at Neutral pH on Polycrystalline Vanadium Sulfide Derived from Vanadium Disulfide
The electrochemical nitrate reduction reaction (NO3RR) offers a pathway to produce NH3 for fuel and fertilizer from waste NO3-. In this work, a polycrystalline vanadium sulfide (VSx), which is derived from solvothermally grown and annealed VS2, is shown to exhibit excellent NO3RR activity (2.3 +- 0.6 mg.cm-2 geo..h-1 @ -0.92 VRHE) and Faradaic efficiency to NH4+ (69 +- 6% at -0.69 VRHE) in buffered neutral pH electrolyte containing 0.1 M NO3-. A variety of characterization techniques are leveraged to support the VSx assignment, including X-ray photoelectron spectroscopy, near-edge X-ray absorption fine structure spectroscopy, selected area electron diffraction, and X-ray diffraction measurements. The VS2 annealing step reduces the oxide character and generates VSx, which, based on the improved NO3RR activity, results in the creation of active sites for NO3- binding. To help shed light on NO3RR on VSx, VS2 is used as a model system, and a grand-canonical density functional theory (GC-DFT) investigation of VS2 shows strong evidence that S vacancies are active sites for NO3RR, where NO3- outcompetes H+ for adsorption at the S-vacancy sites. Moreover, GC-DFT results highlight a thermodynamically favorable reaction to generate NH4+ in an aqueous electrolyte at relevant cathodic potentials. As an annealed material, VSx may contain undersaturated V sites, which show an electronic structure similar to the theoretically calculated S-vacancy site of VS2, and these sites may contribute to the observed increase in NO3RR activity and selectivity for NH4+ on VSx versus unannealed VS2. Finally, kinetic isotope effect measurements suggest that the kinetic rate-limiting step of the NO3RR on VSx is not proton-coupled, indicating it may be the first electron transfer to adsorbed NO3*.
Synergic Effects of Surface Chemistry and Applied Potentials on the Kinetics of the Electrocatalytic Oxygen Evolution Reaction in IrO 2
Improving the efficacy of the oxygen evolution reaction (OER) through water oxidation is critical for advancing photoelectrochemical water splitting. Among the catalysts, IrO 2 exhibits a high OER catalytic activity and stability under acidic conditions. The OER mechanism in this system has been a topic of intense research; however, many mechanistic understandings are lacking. In this work, we carried out first-principles calculations with an implicit solvation model at a constant potential to study the OER behavior on the IrO 2 (110) surface. We find that the surface hydrogen coverage has a significant effect on OER kinetics and transition states. We then develop a microkinetics model that accounts for the continuous evolution of both OER activation energy and hydrogen coverage as a function of an applied potential. We show that this inclusion leads to significant improvement in the simulated Tafel plot compared to available experiments. In conclusion, our results point to a complex interplay between surface chemistry and the applied potential on OER kinetics.
Modeling stochastic fluctuations in relativistic kinetic theory
Using the information current, we develop a Lorentz-covariant framework for modeling equilibrium fluctuations in relativistic kinetic theory in the grand-canonical ensemble. The resulting stochastic theory is proven to be causal and covariantly stable, and its predictions do not depend on the choice of spacetime foliation used to define the grand-canonical probabilities. As expected, in a box containing N > 5 particles, Boltzmann’s molecular chaos postulate is broken with (almost exact) probability N -1/2 , leading to a breakdown of the Boltzmann equation in small systems. Here, we also verify that, in ultrarelativistic gases, transient hydrodynamics already accounts for at least 80% of the equilibrium fluctuations of the stress-energy tensor at a given time. Finally, we compute the correlators at nonequal times for two selected collision kernels: that of a chemically active diluted solution, and that of ultrarelativistic scalar particles self-interacting via a quartic potential. For the former, we compute the density-density correlators analytically in real space, and dehydrodynamization of the stochastic theory is proven to occur whenever the mean free path diverges at high energy.
Extended kinetic lattice grand canonical Monte Carlo simulation method for transport of multicomponent ion mixtures through a model nanopore system
Abstract An extended version of the original kinetic lattice grand canonical Monte Carlo simulation method combined with mean field theory (KLGCMC/MF) ( J. Chem. Phys. 2007 , 127, 024706) is presented for the study of transport of multicomponent ion mixtures through a model nanopore. Comparison of the extended KLGCMC/MF (eKLGCMC/MF) simulation results with Poisson–Nernst–Planck (PNP) calculations is also made to confirm the validity of the extended simulation approach. Unlike the original version of KLGCMC/MF simulation method that treats only a binary ionic solution with one cation and one anion species, this extended version can deal with a system that includes ternary ion mixtures. A diffusion probability algorithm is also added to the extended version of the simulation method to describe the inhomogeneous diffusivity of ions that is often observed in the ion permeation through nanopores. Both Legendre and Chebyshev polynomials of the second kind were tested as a basis set for the basis set expansion (BSE) method with which to calculate the reaction field energy in the eKLGCMC/MF simulation. It turned out that the Legendre polynomials perform better than the Chebyshev polynomials, and as a result, the Legendre polynomials were implemented in the current version of eKLGCMC/MF simulation algorithm. The presented eKLGCMC/MF simulation method with new features finds its potential applications in nanopore systems where the correlation between ion species with the same sign of charges plays a key role such as oscillating ion currents or anomalous mole fraction effects.