Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chemical 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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The effect of blue and infrared laser melting frequency on oxide morphology in 304L

Continuous Wave Laser Beam (LB) melting offers control over the localised heating and cooling of melt pools for welding, brazing, additive manufacturing and solidification. Research into laserliquid metal interactions has primarily focused on the heat and mass transport under large thermal gradients imposed by the localised melting conditions. Recent commercial availability of blue (450 nm) wavelength lasers has enabled the printing of material with poor IR (1064) nm attenuation like Cu. However, little research has been conducted into varying input laser frequency and scan rate at fixed absorption to understand the effects of laser light on surface oxide formation and the underlying liquid metal thermodynamic response. Here, this article conducts laser melting of 304L under Blue (450 nm) and Infrared (IR, 1064 nm) laser frequencies and examines their impact on laser oxide thickness, chemistry and coloration. We find that laser frequencies induce changes in the oxide layer thickness and chemistry that cannot be explained using conventional thermal absorption shifts and fluid dynamics. We suggest that light coupling may have thermodynamic implications for the chemical potential of the liquid metal, which may drive the observed phase behaviour. Identifying mechanisms for the tuning of the chemical potential via laser frequency control could enable the ability to control solidification microstructure and atomic ordering in laser welding and additive manufacturing processes.

304L↗
DOE OSTI2025

Multifunctional Catalysts for the Tandem Reactions of Oxygenates

Industrially-relevant catalytic reactions rarely consist of a simple sequence of elementary steps. Moreover, kinetic coupling of multiple reactions on a catalyst surface is highly desired for process intensification and improved energy efficiency for large scale chemical transformations. The shifting landscape of hydrocarbon chemical feedstocks in the US also motivates research on the selective conversion of more complex molecules. One desirable type of catalytic reaction is the reduction of carboxylic acids that are produced from biomass feedstocks to their corresponding alcohols. The proposed research explores the fundamental importance of hydrogen spillover on a multifunctional catalyst for carboxylic acid reduction with H 2 composed of metal particles coupled to metal oxide particles. Recent work has demonstrated the excellent performance of supported tungsten oxide clusters for carboxylic acid reduction, but only after they are promoted with a late transition metal such as palladium. Elucidating the active state of the catalyst and the associated reaction mechanism for acid reduction on that active state are the overall goals of the proposed project and successful completion will enable future design of efficient multifunctional catalysts. The critically important role of the metal promoter is hypothesized to be its ability to dissociate H2 and spillover atomic H to the support. Although hydrogen spillover is a well-recognized phenomenon in catalysis, its role in both catalyst activation and catalytic turnover are still unresolved. The study combined materials synthesis, characterization, reactivity testing, and molecular simulations, to explore the effect of hydrogen chemical potential on the formation of the active catalytic sites and on the steady state catalytic reduction of carboxylic acid. Varying the hydrogen chemical potential through modification of the gas conditions, support composition, and metal loading to modulated the structure and catalytic performance of the reducible metal oxide. Dual function catalysts containing supported Pd and WO x species co-located on a non-reducible carrier (silica) and a reducible carrier (titania) were synthesized and characterized by electron microscopy, temperature-programmed reduction, and chemisorption. Spectroscopic methods such as X-ray absorption and UV-vis were also used to evaluate the catalysts, which were used in the reduction of propionic acid to aldehyde and alcohol. Quantum chemical calculations, including ab initio phase diagrams provided molecular insights into the H spillover phenomenon.

09 BIOMASS FUELS↗
DOE OSTI2026

Li + /H + Exchange in Solid-State Oxide Li-Ion Conductors

Understanding the moisture stability of oxide Li-ion conductors is important for their practical applications in solid-state batteries. Unlike sulfide or halide conductors, oxide conductors generally better resist degradation when in contact with water but can still undergo topotactic Li + /H + exchange (LHX). Here, we combine density functional theory (DFT) calculations with a machine-learning interatomic potential model to investigate the thermodynamic driving force of the LHX reaction for two representative oxide Li-ion conductor families: garnets and NASICONs. Li-stuffed garnets exhibit a strong driving force for proton exchange due to their high Li chemical potential. In contrast, NASICONs demonstrate a higher resistance against proton exchange due to the lower Li chemical potential and the lower O–H bond covalency for polyanion-bonded oxygens. Our findings reveal a critical trade-off: Li stuffing enhances conductivity but increases moisture susceptibility. This study underscores the importance of designing Li-ion conductors that possess both high conductivity and high stability in practical environments.

Li, Zhuohan [Lawrence Berkeley National Laboratory↗
DOE OSTI2026

Viscous Gubser flow with conserved charges to benchmark fluid simulations

We present semi-analytical solutions for the evolution of both the temperature and chemical potentials for viscous Gubser flow with conserved charges. Such a solution can be especially useful in testing numerical codes intended to simulate relativistic fluids with large chemical potentials. The freeze-out hypersurface profiles for constant energy density are calculated, along with the corresponding normal vectors, and presented as a new unit test for numerical codes. We also compare the influence of the equation of state on the semi-analytical solutions. We benchmark the newly developed smoothed particle hydrodynamics code ccake that includes both shear viscosity and three conserved charges. Here, the numerical solutions are in excellent agreement with the semi-analytical solution and are also able to accurately reproduce the hypersurface at freeze-out.

Hydrodynamic models↗
DOE OSTI2026

Confinement and String Breaking in the Compact Abelian Higgs Model

While real-time simulation of Quantum Chromodynamics remains technologically out of reach, simplified models for studying elements of QCD phenomenology abound. This work presents a simple model, a spin-1 truncation of the Compact Abelian Higgs Model simulated on qutrit sites, in which confinement and string breaking is accessible to current simulation methods. In the low-energy regime of 1+1D scalar electrodynamics, the heavy modes are integrated out, producing a spin chain effective Hamiltonian in which Gauss' law is implicitly satisfied. We study the spectrum of string-like excitations using DMRG methods on the order of 100 sites. We demonstrate that an added, local chemical potential, playing a role analogous to external charges, permits parameter-dependent measurements of physical features of interest like the string tension and effective meson mass. Varying the chemical potential also permits a characterization of string stability not assessed in prior studies of confining lattice models.

Senseman, Blake [Iowa U.]↗
DOE OSTI2026

Machine-Learning-Driven Discovery of Water Splitting BaFe 2 O 4 and Human-in-the-Loop Improvement via Al-Substitution for Increased Thermal Stability

Thermochemical hydrogen (TCH) production offers a promising method for converting thermal energy into hydrogen fuel through heat-driven redox cycles of metal oxides. Here, in this work a defect graph neural network (dGNN) was used to predict oxygen vacancy formation energies ΔH V O combined with Materials Project predictions of oxygen chemical potential stability to screen candidate oxides via high-throughput database analysis. BaFe 2 O 4 was identified as a promising material for experimental validation based on its predicted ΔH V O , oxygen chemical potential stability range, and potential for tunable substitutions to improve thermal properties. Experimental validation using thermogravimetric analysis (TGA), stagnation flow reactor (SFR), X-ray diffraction (XRD), and electron microscopy confirmed positive water-splitting behavior but also revealed limitations in thermal stability under aggressive reduction conditions. To address this, a human-in-the-loop modification strategy was employed introducing Al substitution in BaFe 2–x Al x O 4 ; this modification improves thermal stability, alters the crystal structure and enhances overall performance. These results demonstrate a combined computational and experimental workflow in which machine learning accelerates identification of promising candidates, while targeted experimental design enables optimization of functional performance. This approach advances the development of robust, cost-effective TCH materials and highlights the importance of integrating data-driven discovery with human-guided materials design in paving the way for scalable hydrogen production technologies.

organic↗
DOE OSTI2026

Locating the QCD critical point through contours of constant entropy density

We propose a new method to investigate the existence and location of the conjectured high-temperature critical point of strongly interacting matter via contours of constant entropy density. By approximating these lines as a power series in the baryon chemical potential 𝜇 𝐵 , one can extrapolate them from first-principle results at zero net-baryon density, and use them to locate the quantum chromodynamics (QCD) critical point, including the associated first-order and spinodal lines. As a proof of principle, we employ currently available continuum-extrapolated lattice data from the Wuppertal-Budapest collaboration to find a critical point at a temperature and a baryon chemical potential of 𝑇 𝑐 = 114.3 ± 6.9 MeV and 𝜇 𝐵,𝑐 = 602.1 ± 62.1 MeV, respectively, at expansion order 𝒪⁡(𝜇$^{2}_{𝐵}$). We advocate for a more precise determination of the required expansion coefficients via lattice QCD simulations as a means of pinpointing the location of the critical endpoint in the phase diagram of strongly interacting matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Lattice QCD constraints on the critical point from an improved precision equation of state

In this paper we employ lattice simulations to search for the critical point of QCD. We search for the onset of a first-order QCD transition on the phase diagram by following contours of constant entropy density from imaginary to real chemical potentials under conditions of strangeness neutrality. We scan the phase diagram and investigate whether these contours meet to determine the probability that the critical point is located in a certain region on the 𝑇−𝜇 𝐵 plane. To achieve this we introduce a new, continuum extrapolated equation of state at zero density with improved precision using lattices with 𝑁 𝜏 =8, 10, 12, 16 time slices, and supplement it with new data at imaginary chemical potential. The current precision allows us to exclude, at the 2⁢𝜎 level, the existence of a critical point at 𝜇 𝐵 <450 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Modeling the influence of the solid electrolyte interphase on the sand’s time and dendrite formation on lithium metal electrodes

Lithium metal is a sought after battery material for its high energy density due to the low electrochemical potential and density. However, lithium metal is also highly reactive, which results in a strong propensity for dendrite formation. The Sand’s time has previously been used to predict the time of dendrite initiation on metals that do not form a solid-electrolyte interphase (SEI), but it has been shown that the Sand’s time is not accurate for lithium electrodes when using transport parameters associated with the electrolyte. Thus, we built a numerical model to simulate lithium ion transport through a growing SEI to predict the Sand’s time. The numerical model is shown to be more accurate than previous analytical solutions, especially for low current densities. We then analyze the sensitivity of the Sand’s time to different SEI properties and the chemical potential gradients present in the SEI, driving lithium transport. The results showed that high lithium concentration has a greater impact at high current density, while fast diffusivity is more important at low current density. Lastly, we modeled the influence of surface roughness on the plating evolution and chemical potential gradients when an SEI is present in comparison to the electrolyte. As a result, we demonstrate that the SEI plays a critical role in lithium electrode stability, and that improved characterization techniques are needed to better understand transport through the SEI and increase lithium metal utilization in energy storage devices.

Chemistry↗
DOE OSTI2025

Revisiting the First, Second and Combined Laws of Thermodynamics

The present study revisits the first, second, and combined laws of thermodynamics by introducing partial internal energy in the first law, entropy production and partial entropy in entropy change, and partial volume in the context of work exchange. These developments yield a more rigorous and internally consistent thermodynamic framework, within which chemical potential is defined as the partial derivative of internal energy with respect to composition at constant entropy and volume, and is related to partial internal energy, partial entropy, and partial volume.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling Europium (II/III) ion solvation in the LiCl-KCl eutectic mixture with polarizable force fields

Here, the solvation processes of Europium (II/III) ions within the molten salt eutectic mixture, 3LiCl-2KCl, are investigated over temperatures ranging from 673 K to 1173 K. New polarizable ion force fields are proposed to model Europium ions in a molten salt eutectic mixture with a goal of accurately capturing the underlying physics in the solutions. In contrast to rigid-ion models, the polarizable Drude model with an adjusted WBK (Wang-Buckingham) force field significantly improves predictions for diffusion coefficients, the average diffusion activation energy, and changes in excess ion chemical potential and partial molar entropy, producing good agreement with experimental measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Thermal quasiparticle theory

The widely used thermal Hartree–Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is postulated in consultation with the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasiparticle) energies in such a way that all fundamental thermodynamic relations are obeyed. The associated density matrix is of a one-electron type, whose diagonal elements take the form of the Fermi–Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi–Dirac or thermal HF theory. The theory thus stipulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green’s function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperatures, the theory approaches finite-temperature MBPT of the same order, but it may outperform the latter at intermediate temperatures by including additional electron-correlation effects through orbital energies. Here, a physical meaning of these thermal orbital energies is proposed (encompassing that of thermal HF orbital energies, which has been elusive) as a finite-temperature version of Janak’s theorem.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Computing with a Chemical Reservoir

Contemporary computation is expensive, with large language models and artificial intelligence becoming more common in daily life. However, high-performance computing is reaching the limits in speed and energy expenditure, and domain science requires ever-increasing computational capacity, with simulations and data analysis pipelines ever-growing in complexity. As we progress towards post-exascale computation, with the associated high energy costs, new methods of energy-conscious computation are required. Novel analog and hybrid digital-analog systems can overcome these challenges, and chemical reactions offer a promising avenue. Computers based on chemistry can provide compact desktop devices with immense computational power. These devices are readily scalable by considering greater reaction systems or vessels, meeting the high-performance requirements for scientific workflows. In this article, we present ChemComp, a compilation pipeline for the conversion of ordinary differential equations into implementable chemical reactions. We then demonstrate the solving capabilities of ChemComp by emulating a potential chemical reservoir device. We leverage the multi-layer intermediate representation (MLIR) compiler framework to implement an expressive chemical reaction abstraction and propose a path for chemical reaction networks (CRNs) to represent mathematical problems effectively. Combined, we demonstrate a potential workflow that can harness chemistry’s computing power to create energy-efficient, high-performance computation systems for contemporary computing needs.

artificial intelligence↗
DOE OSTI2025

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↗
DOE OSTI2025

Chemical Process Safety at TRISO-Based, Metal-Based, and Salt-Based Fuel Fabrication Facilities: Technical Assessment and Guidance Assessment

As part of efforts to prepare for potential and ongoing safety reviews for licensing of advanced non-light-water reactor fuel cycles, the U.S. Nuclear Regulatory Commission (NRC) tasked Pacific Northwest National Laboratory to prepare an assessment on the state of knowledge of potential chemical processes at fuel cycle facilities supporting the front end of these fuel cycles, and to assess the associated regulatory guidance. This report provides a technical assessment of chemical process safety considerations to support NRC licensing reviews of fabrication processes for tri-structural isotropic (TRISO) based, metallic-based, and salt-based fuels. The assessments involved collecting publicly available information on the fuel fabrication processes to (i) identify the operational process steps, characteristics and chemicals involved, (ii) identify the physical safety considerations and health safety considerations during licensing reviews of the various process steps, and (iii) collect information to support assessments of severity of accidents and potential mitigative measures to be implemented. The assessment provides a foundational basis on chemical process safety considerations for advanced fuel fabrication activities, although it is recognized that licensing reviews may necessitate design-specific considerations. The specific conditions under which chemical hazards emerge will require process-specific considerations, highlighting the importance of process-informed interpretation. The assessment also determined that exposure guidelines and limits to assess the consequences of acute exposures are limited for some chemicals, although alternative limits and supplementary information from databases or safety data sheets provide sufficient information to evaluate consequences of acute exposures. In addition, it was identified that metallic and salt fuel fabrication processes may involve beryllium, which is an exposure hazard. The regulatory framework for the licensing of advanced fuel cycle facilities, per 10 CFR Part 70 Domestic Licensing of Special Nuclear Material, is deemed robust and flexible to address the chemical safety considerations in this report. A review was conducted on various regulatory guidance and technical basis documents. This included reviewing NUREG-1520, Revision 2, Standard Review Plan for Fuel Cycle Facilities License Applications – Final Report and the process descriptions in Appendix A of NUREG/CR-6410, Nuclear Fuel Cycle Facility Accident Analysis Handbook, to address advanced fuel types. As new fuels will involve process-specific chemical uses, process-specific considerations are provided in this report. Additionally, it is noted that the U.S. Department of Energy protective action criteria database includes Temporary Emergency Exposure Limits (TEELs) for process-specific chemicals. This report provides technical information to support chemical safety assessments of new advanced fuel cycle facilities and identifies technical and safety information to support licensing reviews. No regulatory barriers were identified for the licensing of advanced fuel cycle facilities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Quantitative Analysis of the Semiconductor–Electrolyte Interface Using Cyclic Voltammetry Measurements

Small changes in the chemical potential at a semiconductor interface can result in dramatic changes to the space-charge layer that underpins applications in the electronic and photovoltaic industries as well as in photoelectrochemical cells for fuel production. There has hence been great interest in techniques that directly probe the space-charge layer, yet many fail at the semiconductor–electrolyte interface due to the potential drop in the electric double-layer region of the electrolyte. This article demonstrates that photovoltages, obtained from straightforward cyclic voltammetry measurements, provide an experimental and quantitative approach for characterizing the semiconductor–electrolyte interface. Key parameters accessible through this approach include the flat-band potential ( E fb ), the fraction of the total potential that drops across the space-charge layer (γ sc ) and the electric double layer, as well as the surface recombination lifetime (τ s ). Here, we report photovoltage measurements for p -type Si(111) photoelectrodes in contact with electrolytes containing redox-active species with a range of known reduction potentials that exceed the 1.1 eV bandgap. In tetrabutylammonium [NBu 4 ] + electrolyte, the flat-band potential determined for hydrogen-terminated ( p -Si–H), methyl-terminated ( p -Si–CH 3 ), and chemically oxidized ( p -Si–cSiO x ) surfaces were −0.02, −0.31, and 0.30 V vs Fc +/0 , respectively, agreeing well with expected shifts arising from surface dipole modifications. The quantitative analysis also reveals that 67% of the applied bias drops across the space-charge layer for p -Si–H, 73% for p -Si–CH 3 , and only 44% for p -Si–cSiO x . The remaining potential drop is attributed to the interfacial surface layer, which consists of a molecular dipole or oxide overlayer, and the Helmholtz layer within the electrolyte. When the larger [NBu 4 ] + electrolyte was replaced with Li + , the flat-band position showed minimal changes, but the fraction of the potential drop across the space-charge layer increased significantly, consistent with the small cation altering the structure of the electric double layer.

electrolytes↗
DOE OSTI2025

Topological Mixed Valence Model for Twisted Bilayer Graphene

Song and Bernevig (SB) have recently proposed a topological heavy-fermion description of the physics of magic angle twisted bilayer graphene (MATBG), involving the hybridization of flat-band electrons with a relativistic conduction sea. Here, we explore the consequences of this model, seeking a synthesis of understanding drawn from heavy-fermion physics and MATBG experiments. Our work identifies a key discrepancy between measured and calculated on-site Coulomb interactions, implicating renormalization effects that are not contained in the current model. With these considerations in mind, we consider a SB model with a single, renormalized on-site interaction between the f electrons, containing a phenomenological heavy-fermion binding potential on the moiré A A sites. This feature allows the simplified model to capture the periodic reset of the chemical potential with filling and the observed stability of local moment behavior. We argue that a two-stage Kondo effect will develop in MATBG as a consequence of the relativistic conduction band: Kondo I occurs at high temperatures, establishing a coherent hybridization at the Γ points and a non-Fermi liquid of incoherent fermions at the moiré K points; at much lower temperatures, Kondo II leads to a Fermi liquid in the flat band. Utilizing an auxiliary-rotor approach, we formulate a mean-field treatment of MATBG that captures this physics, describing the evolution of the normal state across a full range of filling factors. By contrasting the relative timescales of phonons and valence fluctuations in bulk heavy-fermion materials with that of MATBG, we are led to propose a valley-polaron origin to the Coulomb renormalization and the heavy-fermion binding potential identified from experiment. We also discuss the possibility that the two-fluid, non-Fermi liquid physics of the relativistic Kondo lattice is responsible for the strange-metal physics observed in MATBG. Published by the American Physical Society 2025

Lau, Liam L. H. (ORCID:0000000166039088)↗
DOE OSTI2025

Machine‐Learning‐Driven Exploration of Surface Reconstructions of Reduced Rutile TiO 2

Abstract Titanium dioxide (TiO 2 ) is widely used as a catalyst support due to its stability, tunable electronic properties, and surface oxygen vacancies, which are crucial for catalytic processes such as the reverse water‐gas shift (RWGS) reaction. Reduced TiO 2 surfaces undergo complex surface reconstructions that endow unique properties but are computationally challenging to describe. In this study, we utilize machine‐learning interatomic potentials (MLIPs) integrated with an active‐learning workflow to efficiently explore reduced rutile TiO 2 surfaces. This approach enabled the prediction of a phase diagram as a function of oxygen chemical potential, revealing a variety of reconstructed phases, including a previously unreported subsurface shear plane structure. We further investigate the electronic properties of these surfaces and validate our results by comparing experimental and theoretical high‐resolution transmission electron microscopy (HRTEM). Our findings provide new insights into how extreme surface reductions influence the structural and electronic properties of TiO 2 , with potential implications for catalyst design.

Lee, Yonghyuk [Chemistry and Biochemistry Universi↗