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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Non-equilibrium states and interactions in the topological insulator and topological crystalline insulator phases of NaCd4As3

Topological materials are of great interest because they can support metallic edge or surface states that are robust against perturbations, with the potential for technological applications. Here, we experimentally explore the light-induced non-equilibrium properties of two distinct topological phases in NaCd4As3: a topological crystalline insulator (TCI) phase and a topological insulator (TI) phase. This material has surface states that are protected by mirror symmetry in the TCI phase at room temperature, while it undergoes a structural phase transition to a TI phase below 200 K. After exciting the TI phase by an ultrafast laser pulse, we observe a leading band edge shift of >150 meV that slowly builds up and reaches a maximum after ∼0.6 ps and that persists for ∼8 ps. The slow rise time of the excited electron population and electron temperature suggests that the electronic and structural orders are strongly coupled in this TI phase. It also suggests that the directly excited electronic states and the probed electronic states are weakly coupled. Both couplings are likely due to a partial relaxation of the lattice distortion, which is known to be associated with the TI phase. In contrast, no distinct excited state is observed in the TCI phase immediately or after photoexcitation, which we attribute to the low density of states and phase space available near the Fermi level. Our results show how ultrafast laser excitation can reveal the distinct excited states and interactions in phase-rich topological materials.

Angle-resolved photoemission spectroscopy↗
DOE OSTI2024

Voltage-Dependent First-Principles Barriers to Li Transport within Li-Ion Battery Solid Electrolyte Interphases

Charging a Li-ion battery requires Li-ion transport between the cathode and the anode. This Li-ion transport is dependent on (among other factors) the electrostatic environment that the ion encounters within the solid electrolyte interphase (SEI), which separates the anode from the surrounding electrolyte. A previous first-principles work has illuminated the reaction barriers through likely atomistic SEI environments but has had difficulty accurately reflecting the larger electrostatic potential landscape that an ion encounters moving through the SEI. In this work, we apply the recently developed quantum continuum approximation (QCA) technique to provide an equilibrium electronic potentiostat for first-principles interface calculations. Using QCA, we calculate the potential barrier for Li-ion transport through LiF, Li 2 O, and Li 2 CO 3 SEIs along with LiF–LiF and LiF–Li 2 O grain boundaries, all paired with Li metal anodes. Here, we demonstrate that the SEI potential barrier is dependent on the electrochemical potentials of the anode in each system. Finally, we use these techniques to estimate the change in the diffusion barrier for a Li ion moving in a LiF SEI as a function of the anode potential. We find that properly accounting for interface and electronic voltage effects significantly lowers reaction barriers compared with previous literature results.

25 ENERGY STORAGE↗
DOE OSTI2023

Neural Network Water Model Based on the MB-Pol Many-Body Potential

The MB-pol many-body potential accurately predicts many properties of water, including cluster, liquid phase, and vapor–liquid equilibrium properties, but its high computational cost can make applying it in large-scale simulations quite challenging. Here, in order to address this limitation, we developed a “deep potential” neural network (DPMD) model based on the MB-pol potential for water. We find that a DPMD model trained on mostly liquid configurations yields a good description of the bulk liquid phase but severely underpredicts vapor–liquid coexistence densities. By contrast, adding cluster configurations to the neural network training set leads to a good agreement for the vapor coexistence densities. Liquid phase densities under supercooled conditions are also represented well, even though they were not included in the training set. These results confirm that neural network models can combine accuracy and transferability if sufficient attention is given to the construction of a representative training set for the target system.

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

Reaching the equilibrium state of the frustrated triangular Ising magnet Ca 3 Co 2 O 6

Ca 3 Co 2 O 6 is a frustrated magnet consisting of a triangular arrangement of chains of Ising spins. It shows regular magnetization steps versus magnetic field every 1.2 T that are metastable with very slow dynamics. This has puzzled the community for many years and given rise to numerous potential theories. Here we approach the problem by seeking the elusive magnetic equilibrium state at T=2 K. To this end, we explore two approaches: (1) bypassing the slow dynamics produced by changing fields by instead field-cooling directly to the target temperature and (2) quantum annealing in transverse magnetic fields. While we observe no measurable effect of the quantum annealing in fields up to 7 T, which is likely due to the large Ising anisotropy of Co spins in this material, we find that for the field cooling in longitudinal fields we achieve the predicted equilibrium 1/3 magnetization. We perform Monte Carlo simulations of the ground state phase diagram and we also simulate the quantum annealing process and find good agreement between experiment and theory. Thus we present an investigation of the elusive ground state properties of the canonical frustrated triangular system Ca 3 Co 2 O 6 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Optimization Methodology of Pebble Bed HTGR Start-Up and Running-in Strategy

In recent years interest in advanced reactor technologies has increased significantly. However, the methods used for analysis of traditional nuclear reactors are insufficient to consider all the different and varied advanced reactor designs without further development. One promising advanced reactor design is the pebble bed reactor (PBR). PBRs possess unique operational and fuel cycle features that require the development of specific analysis methodologies to adequately design and analyze the systems. There is a need in PBR research for a capability to analyze and optimize the process of transitioning from the start-up reactor core to the equilibrium core (known as the “running-in” of the reactor). The start-up of a PBR and the transition to the equilibrium core is a complex, multi-physics challenge that has not yet been well researched and has many opportunities for design, analysis, and optimization of the process. In this research, a methodology is defined to consider the potential strategies in PBR start-up and run-in to the equilibrium core. Multiple candidate software are considered, and their pros and cons are discussed for the PBR-specific application in the methodology. A preliminary software selection for the physics engine is made, and initial verification of key modules is performed. Software to support optimization of the reactor run-in strategies through reduced order modeling (ROM) and machine learning are considered. Challenges for a full implementation of the methodology are discussed. Additional code selections and verification of models relevant to the application are needed before full demonstration of the methodology can be achieved and an optimal strategy determined.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Mean-field theory for the structure of strongly interacting active liquids

Active systems, which are driven out of equilibrium by local non-conservative forces, exhibit unique behaviors and structures with potential utility for the design of novel materials. An important and difficult challenge along the path toward this goal is to precisely predict how the structure of active systems is modified as their driving forces push them out of equilibrium. Here, we use tools from liquid-state theories to approach this challenge for a classic minimal active matter model. First, we construct a nonequilibrium mean-field framework that can predict the structure of systems of weakly interacting particles. Second, motivated by equilibrium solvation theories, we modify this theory to extend it with surprisingly high accuracy to systems of strongly interacting particles, distinguishing it from most existing similarly tractable approaches. Furthermore, our results provide insight into spatial organization in strongly interacting out-of-equilibrium systems.

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

Simulation of divertor heat flux width on EAST by BOUT++ transport code

The BOUT++ edge plasma transport code is applied to study the effects of neutral recycling, drifts and heating scheme on edge plasma profiles and the divertor heat flux width of two discharges of the experimental advanced superconducting tokamak (EAST) steady-state H-mode plasmas heated by low hybrid wave (LHW) and neutral beam (NB), respectively. Neutral recycling seems to have played an important role in the plasma density profile. The edge plasma density drops dramatically for the case w/o neutral recycling, while it can be sustained for the case w/ neutral recycling. Drifts are found to have significant influences on edge plasma profiles. Both the amplitude and width of the divertor heat flux are found to have increased a lot due to drifts. The simulated heat flux width w/ drifts for the two discharges shows reasonable agreement with the experiments. However, the width from the simulation and experiment for the LHW heated discharge is much larger than that of the NB heated discharge. Comparison with Goldston's drift-based model and more detailed analysis on the heat flux contributions from drifts versus turbulence show that drifts are the dominant factor in edge plasma transport for both LHW and NB heated discharges, while turbulence may have played a more important role in determining the heat flux width in LHW heated discharges than that in NB heated discharges. This may account for the larger heat flux width in the LHW heated discharge. The magnetic topology and the equilibrium change by the LHW power may also be a potential reason for the larger heat flux width in the LHW heated discharges, which still needs more evidence and studies to support it. Finally, additional SOLPS simulation w/o drifts for the two discharges show reasonable agreement with the counterpart from the BOUT++ simulation, suggesting the two are both suitable codes for EAST edge plasma simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Radiative pulsed L-mode operation in ARC-class reactors

Abstract A new ARC-class, highly-radiative, pulsed, L-mode, burning plasma scenario is developed and evaluated as a candidate for future tokamak reactors. Pulsed inductive operation alleviates the stringent current drive requirements of steady-state reactors, and operation in L-mode affords ELM-free access to ∼ 90 % core radiation fractions, significantly reducing the divertor power handling requirements. In this configuration the fusion power density can be maximized despite L-mode confinement by utilizing high-field to increase plasma densities and current. This allows us to obtain high gain in robust scenarios in compact devices with P fus > 1000 MW despite low confinement. We demonstrate the feasibility of such scenarios here; first by showing that they avoid violating 0D tokamak limits, and then by performing self-consistent integrated simulations of flattop operation including neoclassical and turbulent transport, magnetic equilibrium, and radiofrequency current drive models. Finally we examine the potential effect of introducing negative triangularity with a 0D model. Our results show high-field radiative pulsed L-mode scenarios are a promising alternative to the typical steady state advanced tokamak scenarios which have dominated tokamak reactor development.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Application of Machine Learning and Data Augmentation Algorithms in the Discovery of Metal Hydrides for Hydrogen Storage

The development of efficient and sustainable hydrogen storage materials is a key challenge for realizing hydrogen as a clean and flexible energy carrier. Among various options, metal hydrides offer high volumetric storage density and operational safety, yet their application is limited by thermodynamic, kinetic, and compositional constraints. In this work, we investigate the potential of machine learning (ML) to predict key thermodynamic properties—equilibrium plateau pressure, enthalpy, and entropy of hydride formation—based solely on alloy composition using Magpie-generated descriptors. We significantly expand an existing experimental dataset from ~400 to 806 entries and assess the impact of dataset size and data augmentation, using the PADRE algorithm, on model performance. Models including Support Vector Machines and Gradient Boosted Random Forests were trained and optimized via grid search and cross-validation. Results show a marked improvement in predictive accuracy with increased dataset size, while data augmentation benefits are limited to smaller datasets and do not improve accuracy in underrepresented pressure regimes. Furthermore, clustering and cross-validation analyses highlight the limited generalizability of models across different material classes, though high accuracy is achieved when training and testing within a single hydride family (e.g., AB2). The study demonstrates the viability and limitations of ML for accelerating hydride discovery, emphasizing the importance of dataset diversity and representation for robust property prediction.

augmentation↗
DOE OSTI2024

Current potential patches

A novel form of the current potential, a mathematical tool for the design of stellarators and stellarator coils, is developed. Specifically, these are current potentials with a finite-element basis, called current potential patches. Current potential patches leverage the relationship between distributions of magnetic dipoles and current potentials to explore limits of the access properties of stellarator coil sets. An example calculation is shown using the Helically Symmetric Experiment (HSX) equilibrium, demonstrating the method's use in coil design and understanding the limits of the access properties of coil sets. Current potential patches have additional desirable properties such as of promoting sparse current sheet solutions and identifying crucial locations of shaping current placement. A result is found for the HSX equilibrium that shaping currents covering only 22% of the winding surface are sufficient to produce the equilibrium to a good accuracy, provided a toroidal field is generated by an exterior coil set.

Elder, Todd↗
DOE OSTI2021

Effect of interfacial structures on phonon transport across atomically precise Si/Al heterojunctions

Phonons are important carriers of energy and information in many cryogenic devices used for quantum information science and in fundamental physics experiments such as dark matter detectors. In these systems phonon behavior can be dominated by interfaces, hence there is increasing demand for a more detailed understanding of interfacial phonon transport in relevant material systems. Previous studies have focused on understanding thermal transport over the entire phonon spectrum at and above room temperature. At ultra-low temperatures, however, knowledge is missing regarding athermal phonon behavior due to the challenge in modeling the extreme conditions in micro-scale, heterogeneous cryogenic systems, as well as extracting single-phonon information from a large ensemble. In this work, we delineate interfacial phonon transport using a combination of classical molecular dynamics (MD) and phonon wave packet simulations, to illustrate the consistency and differences between the ensemble- and single-phonon dynamics. We consider three single-crystal Si surface reconstructions --- (1×1), (v3×v3) and (7×7) --- and model both experimentally observed Si(1×1)/Al interfaces and hypothesized Si(v3×v3)/Al and Si(7×7)/Al interfaces. The overall interfacial thermal conductance calculated from non-equilibrium MD shows that for the Si(1×1)/Al system, the presence of Al twin grain boundaries can hinder phonon transport and reduce thermal conductance by 2--12% relative to single-crystal Al; whereas the Si (v3×v3) and (7×7) reconstructions can enhance it by 6--19%. Normal mode decomposition reveals that both the increase and decrease in conductance are related to inelastic phonon scattering. Single-phonon wave packet simulations show results consistent with NEMD, while further suggesting that phonon polarization conversion is significant even when elastic transmission dominates, and that the interfacial structures have anisotropic impacts on atomic vibrations along different lattice directions. Our findings suggest alternative perspectives in achieving selective phonon transport through controlling interfacial structures of materials with atomically-precise fabrication techniques, and that the phonon wave packet formalism is a potentially powerful method for developing a detailed understanding of non-equilibrium phenomena in the low-temperature limit.

Lu, Zexi↗
DOE OSTI2023

Development of a Systematic and Extensible Force Field for Peptoids (STEPs)

Peptoids (N-substituted glycines) are a class of biomimetic polymers that have attracted significant attention due to their accessible synthesis and enzymatic and thermal stability relative to their naturally occurring counterparts (polypeptides). While these polymers provide the promise of more robust functional materials via hierarchical approaches, they present a new challenge for computational structure prediction for material design. The reliability of calculations hinges on the accuracy of interactions represented in the force field used to model peptoids. For proteins, structure prediction based on sequence and de novo design has made dramatic progress in recent years; however, these models are not readily transferable for peptoids. Current efforts to develop and implement peptoid-specific force fields are spread out, leading to replicated efforts and a fragmented collection of parameterized sidechains. Here, we developed a peptoid-specific force field containing 70 different side chains, using GAFF2 as starting point. The new model is validated based on the generation of Ramachandran-like plots from DFT optimization compared against force field reproduced potential energy and free energy surfaces as well as the reproduction of equilibrium cis/trans values for some residues experimentally known to form helical structures. In conclusion, equilibrium cis/trans distributions (Kct) are estimated for all parameterized residues to identify which residues have an intrinsic propensity for cis or trans states in the monomeric state.

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

Ionic-Based Electrochemical Gas Sensors for Low-Cost, High-Sensitivity SO2 Detection

Sulfur dioxide (SO2) is a toxic gas associated with adverse health and environmental effects that necessitate reliable monitoring techniques. Here, we report the development of an all-solid-state electrochemical sensor utilizing a lithium borate (Li3BO3) solid electrolyte capable of subppm of SO2 detection. While subppm of SO2 sensing has been previously demonstrated in other solid-state electrolyte systems─such as stabilized zirconia, natrium super ionic conductors (NASICON) under mixed-potential conditions─here we establish Li3BO3 as an alternative solid electrolyte enabling equilibrium potentiometric sensing in an all-solid architecture. This sensor demonstrates a detection limit of at least 0.25 ppm, surpassing the human-olfactory threshold and meeting the rigorous requirements for industrial and personal monitoring applications. The sensing mechanism relies on the formation of Li2SO4 on the electrode surface, as evidenced by multimodal characterization techniques, including Raman spectroscopy, scanning electron microscopy (SEM), and scanning transmission electron microscopy (STEM). The strong linear correlation between the open-circuit potential (OCV) and the logarithm of SO2 concentration between 0.25 and 2 ppm indicates that the response is Nernstian in nature.

Lagunas, Francisco (ORCID:000000026377683X)↗
DOE OSTI2022

Ultrafast perturbation of magnetic domains by optical pumping in a ferromagnetic multilayer

Ultrafast optical pumping of systems with spatially nonuniform magnetic textures is known to cause far-from-equilibrium spin transport effects, such as the broadening of domain-walls. Here, we study the dynamics of labyrinth domain networks in ferromagnetic CoFe/Ni multilayers subject to a femtosecond optical pump and find an ultrafast domain dilation by 6% within 1.6 ps. This surprising result is based on the unambiguous determination of a harmonically-related shift of ultrafast magnetic X-ray diffraction for the first- and third-order rings. Domain dilation is plausible from conservation of momentum arguments, whereby inelastic scattering from a hot, quasi-ballistic, radial current transfers momentum to the magnetic domains. Our results suggest a potentially rich variety of unexpected physical phenomena associated with far-from-equilibrium inelastic electron-magnon scattering processes in the presence of spin textures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effects of emitting surfaces and trapped ions on the sheath physics and current flow in multidimensional plasma systems

Recent one-dimensional simulations of planar sheaths with strong electron emission have shown that trapping of charge-exchange ions causes transitions from space-charge limited (SCL) to inverse sheaths. However, multidimensional emitting sheath phenomena with collisions remained unexplored, due in part to high computational cost. As such, we developed a novel continuum kinetic code to study the sheath physics, current flow and potential distributions in two-dimensional unmagnetized configurations with emitting surfaces. For small negatively biased thermionic cathodes in a plasma, the cathode sheath can exist in an equilibrium SCL state. The SCL sheath carries an immense density of trapped ions, neutralized by thermoelectrons, within the potential well of the virtual cathode. For further increases of emitted flux, the trapped ion cloud expands in space. The trapped ion space charge causes an increase of thermionic current far beyond the saturation limit predicted by conventional collisionless SCL sheath models without ion trapping. For sufficiently strong emission, the trapped ion cloud consumes the entire 2D plasma domain, forming a mode with globally confined ions and an inverse sheath at the cathode. In situations where the emitted flux is fixed and the bias is swept (e.g. emissive probe), the trapped ions cause a large thermionic current to escape for all biases below the plasma potential. Strong suppression of the thermionic emission, required for the probe to float, only occurs when the probe is above the plasma potential.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

From PAH to fullerenes: Closing a nano bowl with a nano lid

Quantum chemical DLPNO-CCSD(T)/cc-pVDZ//B3LYP/6–31G(d) calculations have been performed to unravel the formation mechanism of buckminsterfullerene (C 60 ) via the reaction of the C 40 nano bowl (C 40 H 10 ) and corannulene (C 20 H 10 ). The generated potential energy surfaces and molecular properties were further utilized to evaluate equilibrium constants and rate constants of elementary chemical reactions involved in the C 60 synthesis. The overall C 40 H 10 + C 20 H 10 → C 60 + 10H 2 reaction is shown to be highly exoergic and hence thermodynamically favorable at temperatures above 700 K and a kinetically feasible pathway under high-temperature conditions was revealed. This route involves hydrogen atom abstraction steps to activate closed-shell reactants/intermediates alternating with cyclodehydrogenation reactions eventually zipping the edges between the reacting C 40 and C 20 units by consecutive closures of six- and five-membered rings between them. The reaction is driven/catalyzed by hydrogen atoms which carry out the hydrogen abstractions from C 60 H x and are later recovered on the cyclodehydrogenation stages. The initial association reaction, C 40 H 9 + C 20 H 10 → C 60 H 18 + H, linking the two units by a covalent C-C bond appears to represent the kinetic bottleneck of the whole multistep process. This reaction is fast at low temperatures but slows down and has its equilibrium shifted toward the reactants at high temperatures. Here, the feasibility of the proposed mechanism corroborates the hypothesis that C 60 can be produced through a series of elementary reactions typical for the PAH growth.

Equilibrium constant↗
DOE OSTI2024

Imaging moiré excited states with photocurrent tunnelling microscopy

Moiré superlattices provide a highly tuneable and versatile platform to explore novel quantum phases and exotic excited states ranging from correlated insulators to moiré excitons. Scanning tunnelling microscopy has played a key role in probing microscopic behaviours of the moiré correlated ground states at the atomic scale. However, imaging of quantum excited states in moiré heterostructures remains an outstanding challenge. In this report we develop a photocurrent tunnelling microscopy technique that combines laser excitation and scanning tunnelling spectroscopy to directly visualize the electron and hole distribution within the photoexcited moiré exciton in twisted bilayer WS 2 . The tunnelling photocurrent alternates between positive and negative polarities at different locations within a single moiré unit cell. This alternating photocurrent originates from the in-plane charge transfer moiré exciton in twisted bilayer WS 2 , predicted by our GW-Bethe–Salpeter equation calculations, that emerges from the competition between the electron–hole Coulomb interaction and the moiré potential landscape. Our technique enables the exploration of photoexcited non-equilibrium moiré phenomena at the atomic scale.

36 MATERIALS SCIENCE↗