Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Structurally Constrained Evolutionary Algorithm for the Discovery and Design of Metastable Phases

Metastable materials are abundant in nature and technology, showcasing remarkable properties that inspire innovative materials design. However, traditional crystal structure prediction methods, which rely solely on energetic factors to determine a structure’s fitness, are not suitable for predicting the vast number of potentially synthesizable phases that represent a local minimum corresponding to a state in thermodynamic equilibrium. Here, we present a new approach for the prediction of metastable phases with specific structural features, and interface this method with the XTALOPT evolutionary algorithm. Our method relies on structural features that include the local crystalline order (e.g., the coordination number or chemical environment), and symmetry (e.g., Bravais lattice and space group) to filter the breeding pool of an evolutionary crystal structure search. The effectiveness of this approach is benchmarked on three known metastable systems: XeN 8 , with a two-dimensional polymeric nitrogen sublattice, brookite TiO 2 , and a high pressure BaH 4 phase that was recently characterized. Additionally, a newly predicted metastable melaminate salt, P1¯WC 3 N 6 , was found to possess an energy that is lower than two phases proposed in a recent computational study. Here, the method presented here could help in identifying the structures of compounds that have already been synthesized, and developing new synthesis targets with desired properties.

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

Nonadiabatic Molecular Dynamics Study of the Relaxation Pathways of Photoexcited Cyclooctatetraene

In the current study, we present nonadiabatic (NAMD) and adiabatic molecular dynamics simulations of the transition-state dynamics of photoexcited cyclooctatetraene (COT). The equilibrium-state structure and absorption spectra are analyzed using the semiempirical Austin Model 1 potential. The NAMD simulations are obtained by a surface-hopping algorithm. We analyzed in detail an active excited to ground state relaxation pathway accompanied by an S 2 /S 3 (D 2d ) → S 1 (D 8h ) → S o (D 4h ) → S o (D 2d ) double-bond shifting mechanism. The simulated excitation lifetime is in good agreement with experiment. The first excited singlet state S1 plays a crucial role in the photochemistry. Here, the obtained critical molecular conformations, energy barrier, and transition-state lifetime results will provide a basis for further investigations of the bond-order inversion and photoswitching process of COT.

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

Spin-Polarization Simulations for the Future Circular Collider $E^+E^-$ using BMAD

The high precision measurement of the centre-of-mass energy in the Future Circular Collider e+e- (FCC-ee) at Z and W energies can be realized through resonant spin depolarization utilizing transversely polarized beams. This requires a guaranteed sufficiently-high spin polarization in the presence of lattice imperfections. Investigations of the impact of misalignments on the equilibrium polarization are conducted using analytical and Monte-Carlo spin simulations with Bmad. Potential optimization schemes to ensure high polarization using orbit bumps have been explored.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Depletion-driven thermochemistry of molten salt reactors: review, method, and analysis

Molten salt reactors (MSRs) are innovative advanced nuclear reactors that utilize nuclear fuel by dissolving it in a high-temperature liquid salt. This unique feature differentiates MSRs from other types of reactors and allows for enhanced safety and economic performance. The liquid fuel also entails several multiphysics effects that can complicate reactor design and operation. One primary effect termed here as depletion-driven thermochemistry is a driving force in altering the multiphysics behavior of the reactor. Essentially, depletion-driven thermochemistry is the effect that fuel depletion has on changing the chemical redox potential of the fuel salt over time. As the fuel is consumed, the redox potential shifts toward a more oxidizing state. Without active control, the changing chemistry due to depletion increases corrosion thereby limiting reactor component lifetimes. Additionally, the changing redox potential of the fuel salt alters the vapor pressures of chemical species dissolved in the fuel salt. Changing vapor pressures of species in the fuel salt is an important parameter to understand when off-gassing volatile species during normal reactor operation, and for source term characterization during accident scenario transients. The present work represents a fundamental step toward modeling and coupling the driving physics (i.e., neutronics and chemistry) involved in altering the redox potential in an MSR. Here, the neutronic code Griffin models the depletion of the fuel-salt system, while the chemical equilibrium code Thermochimica calculates the thermochemical state of the isotopic inventory, using the Molten Salt Thermodynamic Database - Thermochemical (MSTDB-TC). These two codes are tightly coupled to predict the impact of fuel depletion in altering the chemistry in MSR systems. Redox potential control methods are discussed and can be modeled using this multiphysics approach. The vapor pressures of chemical species that could be extracted to an off-gas system, as determined by the reactor’s thermochemical state, are examined. The neutronics-chemistry coupling developed in this work is expected to have potential application for analyzing corrosion, source term evolution, and material safeguards in MSR systems. Lastly, suggestions for areas of further improvements of the models to expand these capabilities by incorporating other coupled physics effects is provided.

Walker, Samuel A.↗
DOE OSTI2025

Machine learned potential for high-throughput phonon calculations of metal—organic frameworks

Metal–organic frameworks (MOFs) are highly porous and versatile materials studied extensively for applications such as carbon capture and water harvesting. However, computing phonon-mediated properties in MOFs, like thermal expansion and mechanical stability, remains challenging due to the large number of atoms per unit cell, making traditional Density Functional Theory (DFT) methods impractical for high-throughput screening. Recent advances in machine learning potentials have led to foundation atomistic models, such as MACE-MP-0, that accurately predict equilibrium structures but struggle with phonon properties of MOFs. In this work, we developed a workflow for computing phonons in MOFs within the quasi-harmonic approximation with a fine-tuned MACE model, MACE-MP-MOF0. The model was trained on a curated dataset of 127 representative and diverse MOFs. The fine-tuned MACE-MP-MOF0 improves the accuracy of phonon density of states and corrects the imaginary phonon modes of MACE-MP-0, enabling high-throughput phonon calculations with state-of-the-art precision. The model successfully predicts thermal expansion and bulk moduli in agreement with DFT and experimental data for several well-known MOFs. These results highlight the potential of MACE-MP-MOF0 in guiding MOF design for applications in energy storage and thermoelectrics.

Elena, Alin Marin↗
DOE OSTI2020

Correlations Between Fission Fragment and Neutron Anisotropies in Neutron-Induced Fission

Several sources of angular anisotropy for fission fragments and prompt neutrons have been studied in neutron-induced fission reactions. These include kinematic recoils of the target from the incident neutron beam and the fragments from the emission of the prompt neutrons, preferential directions of the emission of the fission fragments with respect to the beam axis due to the population of particular transition states at the fission barrier, and forward-peaked angular distributions of pre-equilibrium neutrons which are emitted before the formation of a compound nucleus. In addition, there are several potential sources of angular anisotropies that are more difficult to disentangle: the angular distributions of prompt neutrons from fully accelerated fragments or from scission neutrons, and the emission of neutrons from fission fragments that are not fully accelerated. Here, we study the effects of the first group of anisotropy sources, particularly exploring the correlations between the fission-fragment anisotropy and the resulting neutron anisotropy. While kinematic effects were already accounted for in our Hauser-Feshbach Monte Carlo code cgmf, anisotropic angular distributions for the fission fragments and pre-equilibrium neutrons resulting from neutron-induced fission on 233,234,235,238 U, 239,241 Pu, and 237 Np have been introduced for the first time. The effects of these sources of anisotropy are examined over a range of incident neutron energies, from thermal to 20 MeV, and compared with experimental data from the Chi-Nu liquid scintillator array. The anisotropy of the fission fragments is reflected in the anisotropy of the prompt neutrons, especially as the outgoing energy of the prompt neutrons increases, allowing for an extraction of the fission-fragment anisotropy to be made from a measurement of the neutrons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Ultrasonic-Assisted Extrusion Processing for Enhancing Physical Properties of High-Density Polyethylene by Flow-Induced Crystallization

The evolution of crystallinity resulting from stress imposed on a melt, known as flow-induced crystallinity, can strongly influence the mechanical and physical properties of semicrystalline polymers. This study investigates shear-induced crystallization by applying an ultrasonic field to the melt flow as it passes through dies with various geometries. A custom-built sonication die is employed for controlling the dynamic temperature and shear environment, resulting in molecular alignment and potential for flow-induced crystallization. Application of both conventional and ultrasonic shear rates at the equilibrium melt temperature of high-density polyethylene (HDPE) was investigated to accelerate crystallinity and manipulate the crystal morphology across the film in pursuit of improved mechanical and gas barrier properties without the need for additives or other polymer layers. The relationships among ultrasonic-assisted extrusion processing, polymer structure, and performance were analyzed using wide- and small-angle X-ray scattering (WAXS and SAXS), tensile testing, and oxygen transmission rate (OTR) analysis. Multiple linear regression models were implemented to predict the correlation among HDPE structure, process, and properties. Structural analysis revealed that both conventional and ultrasonic shear rates had the most significant influence on lamellar spacing and redistribution of rigid and soft amorphous fractions within the crystalline domains, ultimately dictating the mechanical and physical properties of the films. The goal is to explore the potential of the ultrasonic-assisted high crystallinity monolayer that can replace some of the functionality of complex, heterogeneous multilayer packaging with a single-material film having enhanced oxygen barrier properties.

crystallinity↗
DOE OSTI2021

Evidence for equilibrium exciton condensation in monolayer WTe 2

We present evidence that the two-dimensional bulk of monolayer WTe2 contains electrons and holes bound by Coulomb attraction—excitons—that spontaneously form in thermal equilibrium. On cooling from room temperature to 100 K, the conductivity develops a V-shaped dependence on electrostatic doping, while the chemical potential develops a step at the neutral point. These features are much sharper than is possible in an independent-electron picture, but they can be accounted for if electrons and holes interact strongly and are paired in equilibrium. Our calculations from first principles show that the exciton binding energy is larger than 100 meV and the radius as small as 4 nm, explaining their formation at high temperature and doping levels. Below 100 K, more strongly insulating behaviour is seen, suggesting that a charge-ordered state forms. The observed absence of charge density waves in this state is surprising within an excitonic insulator picture, but we show that it can be explained by the symmetries of the exciton wavefunction. Therefore, in addition to being a topological insulator, monolayer WTe2 exhibits strong correlations over a wide temperature range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Actinium–DOTA coordination in water from hybrid ML/MM: Structure, free energies, and water-exchange pathways

Quantitative simulation of trivalent ƒ-block chelates in water remains challenging because bonded and non-bonded force-field models make different approximations for coordination structure, exchange dynamics, and ion–ligand interactions in highly charged systems. Here, we develop a hybrid machine-learning/molecular-mechanics (ML/MM) framework for Ac 3+ –DOTA in explicit solvent by training an E(3)-equivariant neural network potential (MACELES) on mechanically embedded QM/MM data for Ac aquo and Ac–DOTA species and coupling it to NAMD 2.14 with particle-mesh Ewald electrostatics. Nanosecond ML/MM trajectories remain numerically stable and preserve chelate integrity, yielding a compact DOTA inner shell with an inner-sphere water coordination number of CN Ac,O w ≈ 1.7 arising from a dynamic equilibrium between one- and two-water states (37.5% and 59.9% of frames; three waters 2.5%). A 5 ns potential of mean force shows two low-lying basins at CN Ac,O w ≈ 1 and CN Ac,O w ≈ 2. DFT end-state free energies are consistent with the ML/MM profile, and DFT minimum-energy paths provide a qualitative electronic-structure reference for the observed basin connectivity. State-resolved kinetics reveal picosecond water-exchange pathways that couple hydration changes to transient DOTA arm fluctuations, and training-set comparisons show that temperature-matched Ac–DOTA data optimize energy/force accuracy while more diverse solvated data improve charge prediction. Overall, the present hybrid ML/MM model provides a practical description of Ac 3+ –DOTA hydration thermodynamics and short-time exchange behavior in explicit water at MD-like cost.

Actinium↗
DOE OSTI2021

Inverse design of equilibrium cluster fluids applied to a physically informed model

Inverse design strategies have proven highly useful for the discovery of interaction potentials that prompt self-assembly of a variety of interesting structures. However, often the optimized particle interactions do not have a direct relationship to experimental systems. In this work, we show that Relative Entropy minimization is able to discover physically meaningful parameter sets for a model interaction built from depletion attraction and electrostatic repulsion that yield self-assembly of size-specific clusters. Here, we then explore the sensitivity of the optimized interaction potentials with respect to deviations in the underlying physical quantities, showing that clustering behavior is largely preserved even as the optimized parameters are perturbed.

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

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗
DOE OSTI2025

Operando Characterization of Fe in Doped Ni x (Fe 1– x )O y H z Catalysts for Electrochemical Oxygen Evolution

Iron-doped nickel oxyhydroxides, Ni x (Fe 1– x )O y H z , are among the most promising oxygen evolution reaction (OER) electrocatalysts in alkaline environments. Although iron (Fe) significantly enhances the catalytic activity, there is still no clear consensus on whether Fe directly participates in the reaction or merely acts as a promoter. To elucidate the Fe’s role, we performed operando X-ray spectroscopy studies supported by DFT on Ni x (Fe 1– x )O y H z electrocatalysts. We probed the reversible changes in the structure and electronic character of Ni x (Fe 1– x )O y H z as the electrode potential is cycled between the resting (here at 1.10 V RHE ) and operational states (1.66 V RHE ). DFT calculations and XAS simulations on a library of Fe structures in various NiO y H z environments are in favor of a distorted local octahedral Fe(III)O 3 (OH) 3 configuration at the resting state with the NiO y H z scaffold going from α-Ni(OH) 2 to γ-NiOOH as the potential is increased. Under catalytic conditions, EXAFS and HERFD spectra reveal changes in p-d mixing (covalency) relative to the resting state between O/OH ligands and Fe leading to a shift from octahedral to square pyramidal coordination at the Fe site. XES measurements and theoretical simulations further support that the Fe equilibrium structure remains in a formal Fe(III) state under both resting and operational conditions. These spectral changes are attributed to potential dependent structural rearrangements around Fe. The results suggest that ligand dissociation leads to the C 4v symmetry as the most stable intermediate of the Fe during OER. This implies that Fe has a weakly coordinated or easily dissociable ligand that could serve to coordinate the O–O bond formation and, tentatively, play an active role in the Ni x (Fe 1– x )O y H z electrocatalyst.

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

A short review of defect superlattice formation in metals and alloys under irradiation

Irradiation damage drives complex and coupled phenomena in materials at far-from-equilibrium conditions. The self-organization of nanoscale defects in materials under irradiation shows great potential to tailor the physical properties of materials by controlling nanopatterned microstructures. Irradiation-induced gas bubble and void superlattices are two important ordered nanostructures of great scientific interest. Although both types of superlattices have been investigated extensively, a consensus has yet to be reached on their formation mechanisms. In this review article, the current research status of gas bubble and void superlattices in metals and alloys and their characterization, structural stability, and mechanistic modeling are summarized. The fundamental research goals to advance the mechanistic understanding of gas bubble and void superlattices are outlined.

36 MATERIALS SCIENCE↗
DOE OSTI2025

On the Sampling-Based Computation of Nash Equilibria Under Uncertainty via the Nikaido–Isoda Function

We consider the computation of an equilibrium of a stochastic Nash equilibrium problem, where the player objectives are assumed to be L 0 -Lipschitz continuous and convex, given rival decisions with convex and closed player-specific feasibility sets. To address this problem, we consider minimizing a suitably defined value function defined using the Nikaido–Isoda function. Such an avenue does not necessitate either monotonicity properties of the concatenated gradient map or potentiality requirements on the game but does require a suitable regularity requirement under which a stationary point is a Nash equilibrium. We design and analyze a sampling-enabled projected-gradient-response method, reliant on inexact resolution of a player-level best-response subproblem. Here, by deriving suitable Lipschitzian guarantees on the value function, we derive both asymptotic guarantees for the sequence of generated iterates as well as rate and complexity guarantees for computing a stationary point by appropriate choices of the sampling rate and inexactness sequence.

Nikaido-Isoda function↗
DOE OSTI2021

Calculation and interpretation of classical turning surfaces in solids

Abstract Classical turning surfaces of Kohn–Sham potentials separate classically allowed regions (CARs) from classically forbidden regions (CFRs). They are useful for understanding many chemical properties of molecules but need not exist in solids, where the density never decays to zero. At equilibrium geometries, we find that CFRs are absent in perfect metals, rare in covalent semiconductors at equilibrium, but common in ionic and molecular crystals. In all materials, CFRs appear or grow as the internuclear distances are uniformly expanded. They can also appear at a monovacancy in a metal. Calculations with several approximate density functionals and codes confirm these behaviors. A classical picture of conduction suggests that CARs should be connected in metals, and disconnected in wide-gap insulators, and is confirmed in the limits of extreme compression and expansion. Surprisingly, many semiconductors have no CFR at equilibrium, a key finding for density functional construction. Nonetheless, a strong correlation with insulating behavior can still be inferred. Moreover, equilibrium bond lengths for all cases can be estimated from the bond type and the sum of the classical turning radii of the free atoms or ions.

Chemistry↗
DOE OSTI2026

Nonlinear thermodynamic computing out of equilibrium

We present the design for a thermodynamic computer that can perform arbitrary nonlinear calculations in or out of equilibrium. Simple thermodynamic circuits, fluctuating degrees of freedom in contact with a thermal bath and confined by a quartic potential, display an activity that is a nonlinear function of their input. Such circuits can therefore be regarded as thermodynamic neurons, and can serve as the building blocks of networked structures that act as thermodynamic neural networks, universal function approximators whose operation is powered by thermal fluctuations. We simulate a digital model of a thermodynamic neural network, and show that its parameters can be adjusted by genetic algorithm to perform nonlinear calculations at specified observation times, regardless of whether the system has attained thermal equilibrium. This work expands the field of thermodynamic computing beyond the regime of thermal equilibrium, enabling fully nonlinear computations, analogous to those performed by classical neural networks, at specified observation times.

Whitelam, Stephen [Lawrence Berkeley National Labo↗
DOE OSTI2025

Equilibrium expectations for non-Gaussian fluctuations near a QCD critical point

With the highly anticipated results from the Beam Energy Scan II program at RHIC being recently revealed, an understanding of particle-number fluctuations and their significance as a potential signature of a possible QCD critical point is crucial. Early works that embarked on this endeavor sought to estimate the fluctuations due to the presence of a critical point assuming they stay in equilibrium. From these results came the proposal to focus efforts on higher, non-Gaussian, moments of the event-by-event distributions, in particular of the number of protons. These non-Gaussian moments are especially sensitive to critical fluctuations, as their magnitudes are proportional to high powers of the critical correlation length. As the equation of state provides key input for hydrodynamical simulations of heavy-ion collisions, we estimate equilibrium fluctuations from the BEST equation of state (EoS) that includes critical features from the 3D Ising Model. In particular, the proton factorial cumulants and their dependence on non-universal mapping parameters is investigated within the BEST EoS. Furthermore, the correlation length, as a central quantity for the assessment of fluctuations in the vicinity of a critical point, is also calculated in a consistent manner with the scaling equation of state. An understanding of the equilibrium estimates of proton factorial cumulants will be useful for further comparison to estimates of out-of-equilibrium fluctuations in order to determine the magnitude of the observable fluctuations to be expected in heavyion collision experiments, in which the time spent near a critical point is short.

Karthein, Jamie M. [Massachusetts Institute of Tec↗
DOE OSTI2023

The nature and suppression strategies of interfacial reactions in all-solid-state batteries

Solid-state Li batteries are promising energy storage devices owing to their high safety and high theoretical energy density. However, the serious interfacial reaction between solid state electrolytes and cathodes deteriorates the battery performance, impeding the realization of long-term cyclability. The buried nature of the interphase layer presents a significant challenge in achieving a comprehensive understanding of the underlying interfacial reaction mechanisms. Herein, we systematically explore the interfacial reaction evolutions and interphase compositions and electronic properties between the popular oxide cathodes and sulfide solid electrolytes (SSEs). This includes analysis of the chemical and electrochemical reactions between cathodes/coatings and SSEs, as well as the electrochemical self-decomposition of SSEs by thermodynamic phase equilibrium analysis. Herein, we disclose that the driving force of the electrochemical reaction at the chemical potential of μ Li is much stronger than that of the chemical reaction, which dominates the interfacial reaction. Preventing the formation of an electronically conductive interphase is crucial in inhibiting the continuous interfacial degradation during long-term cycling, which can be achieved through the optimized combination of cathodes and SSEs, as well as the introduction of functional coatings between them. Based on these findings, the percentage of molar fraction ( f ) of electronically conductive species in the formed interphase is proposed as a key factor for indicating the interfacial stability for the first time. Furthermore, we propose a specific high-throughput screening scheme to filter the functional coating materials by comprehensively evaluating their functionality. The tiered screening identifies 48 coating materials with optimal properties. The work highlights the significant roles of rational coupling of the cathodes and SSEs, and optimizing interfacial coating materials for solid-state batteries. It opens new avenues for engineering an interphase with improved interfacial compatibility to realize long-term cyclability.

25 ENERGY STORAGE↗