Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grand Potential”

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 37 records · Page 2

Mesoscale Modeling for Restructuring and Fragmentation in High Burnup UO 2

This report summarizes the mesoscale modeling conducted in fiscal year 2025 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, focusing on the microstructural evolution and restructuring in high burnup UO 2 nuclear fuel and its impact on fuel fragmentation. We developed a pioneering phase-field model to simulate restructuring behavior across different regions of high burnup fuel, including the dark zone and rim region. A grand-potential-based phase-field model is employed to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of temperature and burnup rate were studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to initiate around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. The restructuring occurs faster with higher burnup rate. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. The predicted microstructures are consistent with experimental observations of the restructured regions. Based on the mesoscale simulations, a mechanistic model for restructuring and grain size evolution was implemented in BISON. Thus, this work provides a first-of-its-kind restructuring model for different regions of high-burnup fuel to BISON. It was shown that the model predicts appropriate grain size evolution along fuel radius as those observed in experiments. This model enables BISON to predict the effect of restructuring on the fission gas release. Finally, the phase-field fracture simulation with dark zone specific microstructures were presented to provide the fragmentation criteria for the dark zone. This work introduces a first-of-its-kind restructuring model for high burnup fuel, significantly enhancing BISON's predictive capabilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Phase-field modeling of radiation-induced segregation for multicomponent alloys

Structural alloys under irradiation are known to undergo radiation-induced solute redistribution (RIS) at grain boundaries, leading to detrimental effects such as intergranular corrosion and stress-assisted cracking. To better understand the phenomenon, improved models of RIS applicable to concentrated, multicomponent alloys, and mesoscale microstructures are needed. In this talk, we present a novel grand-potential-based phase-field model to account for the complete set of multicomponent kinetic and thermodynamic couplings between atoms and point defects in the Onsager transport equations. We demonstrate multiscale modeling capability by deriving the Onsager coefficient matrix from atomistic-based Kinetic Monte Carlo simulations. Model predictions and validations of RIS and the effect of defect production, grain boundary sink strength and density are demonstrated for a model FCC FeCrNi system. We also demonstrate the novel capability to describe RIS in the presence of equilibrium segregation described using a density-based CALPHAD thermodynamic framework. This modeling approach overcomes certain limitations in conventional RIS models and provides a step closer towards multiscale modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulating water adsorption in metal–organic frameworks with open metal sites using the 12-6-4 Lennard–Jones potential

Metal-organic frameworks (MOFs) with coordinatively unsaturated open metal sites (OMSs) are promising sorbent materials in atmospheric water harvesting (AWH) systems at low relative humidity (RH) due to their strong interactions with water molecules. However, accurate computational modelling of water adsorption in those materials are challenging, as standard force fields (FFs) based on the 12-6 Lennard-Jones (L-J) potential cannot properly describe the water-OMS interactions. In biomolecular simulations, the 12-6-4 L-J potential has been successfully used to model metal ion-water interactions in which the extra 1/r 4 -term is for charge-induced dipole electrostatics. In this work, we adopted this strategy and used the 12-6-4 L-J potential to model water-OMS interactions in MOFs for low RH AWH applications. Notably, without any modifications, the parameters used for aqueous metal ions were able to greatly improve the accuracy in predicting water adsorption isotherms in MOF-74, which highlights the simplicity and transferability of this method.

74 ATOMIC AND MOLECULAR PHYSICS↗

Dynamical tuning of the chemical potential to achieve a target particle number in grand canonical Monte Carlo simulations

We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system. Corrections to the chemical potential are made according to time-averaged estimates of the mean and variance of the particle number, with the latter being proportional to thermodynamic compressibility. Here, we perform a variety of tests, and in all cases find rapid convergence of the chemical potential—inexactness of the tuning algorithm contributes only a minor part of the total measurement error for realistic simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

Mu2e Experiment

The goal of the Mu2e experiment is right in the name. Mu2e (pronounced “mew-to-ee”) is shorthand for “muon to electron.” Scientists are looking for the one-to-one conversion of a particle called the muon into its more familiar relative, the electron. New theories predict this never observed transformation, and Mu2e will be 10,000 times more sensitive than previous attempts to discover it. A successful sighting could hint at undiscovered particles and potentially illuminate a grand unification theory of nature’s four forces. The experiment complements research at the Large Hadron Collider in Europe, potentially clarifying the origins of new particles.

Fermilab, Fermilab↗

Grand unification at the cosmological collider with chemical potential

We introduce a tree-level chemical potential mechanism for spin-1 particles within cosmological collider physics, allowing them to be detected in primordial non-Gaussianities for masses above the inflationary Hubble scale. We apply this mechanism to orbifold grand unification and the massive unification partners of the standard model gauge bosons. Our mechanism requires at least a pair of massive vector fields which are singlets of the standard model, a condition which is satisfied in the classic “trinification” scenario. Assuming that the gauge hierarchy problem is solved by supersymmetry, gauge coupling running points to unification partners at ~ 10$^{15}$ GeV. We show that, within high-scale inflation, chemical potential enhancement can lead to observably strong signals for trinification partners in future cosmological surveys.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Understanding protein-complex assembly through grand canonical maximum entropy modeling

Inside a cell, heterotypic proteins assemble in inhomogeneous, crowded systems where the abundance of these proteins vary with cell types. While some protein complexes form putative structures that can be visualized with imaging, there are far more protein complexes that are yet to be solved because of their dynamic associations with one another. Nevertheless, it is possible to infer these protein complexes through a physical model. However, it is often not clear to physicists what kind of data from biology is necessary for such a modeling endeavor. Here, we aim to model these clusters of coarse-grained protein assemblies from multiple subunits through the constraints of interactions among the subunits and the chemical potential of each subunit. We obtained the constraints on the interactions among subunits from the known protein structures. We inferred the chemical potential that dictates the particle number distribution of each protein subunit from the knowledge of protein abundance from experimental data. Guided by the maximum entropy principle, we formulated an inverse statistical mechanical method to infer the distribution of particle numbers from the data of protein abundance as chemical potentials for a grand canonical multicomponent mixture. Using grand canonical Monte Carlo simulations, we captured a distribution of high-order clusters in a protein complex of succinate dehydrogenase with four known subunits. The complexity of hierarchical clusters varies with the relative protein abundance of each subunit in distinctive cell types such as lung, heart, and brain. When the crowding content increases, we observed that crowding stabilizes emergent clusters that do not exist in dilute conditions. We, therefore, proposed a testable hypothesis that the hierarchical complexity of protein clusters on a molecular scale is a plausible biomarker of predicting the phenotypes of a cell.

59 BASIC BIOLOGICAL SCIENCES↗

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.

boron-doped copper↗

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗

Electrocatalytic Hydrogen Evolution at Full Atomic Utilization over ITO-Supported Sub-nano-Pt n Clusters: High, Size-Dependent Activity Controlled by Fluxional Pt Hydride Species

A combination of density functional theory (DFT) and experiments with atomically size-selected Pt n clusters deposited on indium-tin oxide (ITO) electrodes was used to examine the effects of applied potential and Pt n size on the electrocatalytic activity of Pt n (n = 1, 4, 7, 8) for the hydrogen evolution reaction (HER). Activity is found to be negligible for isolated Pt atoms on ITO, increasing rapidly with Pt n size, such that Pt 7 /ITO and Pt 8 /ITO have roughly double the activity per Pt atom compared to atoms in the surface layer of polycrystalline Pt. Both DFT and experiment find that hydrogen under-potential deposition (H upd ) results in Pt n /ITO (n = 4, 7, 8) adsorbing ~2 H atoms/Pt atom at the HER threshold potential, equal to ca. double the Hupd observed for Pt bulk or nanoparticles. Here, the cluster catalysts under electrocatalytic conditions are hence best described as a Pt hydride compound, significantly departing from a metallic Pt cluster. The exception is Pt 1 /ITO, where H adsorption at the HER threshold potential is energetically unfavorable. Theory combines global optimization with grand canonical approaches for the influence of potential, uncovering that several metastable structures contribute to HER, changing with the applied potential. It is hence critical to include reactions of the ensemble of energetically accessible Pt n H x /ITO structures to correctly predict the activity vs. Pt n size and applied potential. For the small clusters, spillover of H ads from the clusters to the ITO support is significant, resulting in a competing channel for loss of H ads , particularly at slow potential scan rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The current status of high temperature electrochemistry-based CO 2 transport membranes and reactors for direct CO 2 capture and conversion

The concept of direct CO 2 capture and conversion has attracted significant interest from industries and academia in recent decades due to its potential to address the current grand challenge of global warming/climate change, rapid depletion of fossil fuels and realization of a future carbon neutral ecosystem. The incumbent benchmark technology for CO 2 capture is the post-combustion flue-gas “amine washing”, which is energy intensive and costly for large-scale commercial implementation. The CO 2 conversion technologies, on the other hand, are still at their infancy with many technical challenges to overcome, but primarily being explored in laboratory-scale, low-temperature, solution-based and high-temperature, solid-oxide-based electrochemical cells with renewable electricity perceived as the energy input. In this article, we provide a comprehensive overview on an emergent class of high-temperature electrochemical CO 2 transport membranes that can capture and convert CO 2 into valuable chemicals in single catalytic reactor fashion. The review starts with the chemistry and transport theory of three basic types of membranes purposely designed for different CO 2 feedstocks and downstream conversions. A range of key functional materials used in these membranes and their microstructural/electrochemical properties important to the CO 2 transport are then thoroughly discussed in conjunction with the effects of surface modifications and operating conditions. Several types of combined CO 2 capture and conversion catalytic reactors based on these membranes are also assessed with a focus on their working principles, system configurations and performance demonstrations. Finally, challenges and prospective of these electrochemical CO 2 transport membranes and their associated conversion reactors are candidly discussed for future development.

10 SYNTHETIC FUELS↗