Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermodynamics of computing”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Bayesian Inference applied to the Solubility of Phase-3 and Phase-5

This Section covers an introduction to the objectives and techniques used in this analysis. The objectives of the report are given in Subsection 1.1. An introduction to aqueous thermodynamics and how variance might propagate through the relevant thermodynamic equations is given in Subsection 1.2. An introduction to Bayesian inference and its application to thermodynamic modeling is given in Subsection 1.3.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Interactive Exploration of High-Dimensional Phase Diagrams

High-dimensional thermodynamic phase stability databases are becoming increasingly common due to the convergence of three recent trends: (i) the widespread interest in so-called “high-entropy” alloys, (ii) the availability of high-throughput computational assessments of phase stability in broad composition spaces and (iii) the ongoing development of ever-increasingly broad, multicomponent, multiphase CALPHAD databases. Although automated computational tools can readily process such high-dimensional data, scientists are often unable to visualize the relevant phase relations, an ability that is crucial to gaining an intuitive understanding of the stability constraints governing materials design. The present work addresses this need by providing algorithms that enable the interactive exploration of phase equilibria in high-dimensional spaces. These algorithms concentrate the complex nonlinear nonsmooth optimization needed into a preprocessing step that generates a large number of high-dimensional yet elementary graphical primitives. Furthermore, these primitives can then be cross-sectioned to yield 3-dimensional views in a computationally efficient manner that enables an interactive exploration of high-dimensional spaces. All of these operations are highly parallelizable, thus facilitating scaling of this method to large data sets.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Recombination of 𝐵 𝑐 mesons in ultrarelativistic heavy-ion collisions

High-energy heavy-ion collisions have been suggested as a favorable environment for the production of 𝐵 𝑐 mesons, due to a much larger abundance of charm and bottom quarks compared to elementary reactions. Motivated by recent CMS data for 𝐵$^{+}_{𝑐}$ production in Pb-Pb(5.02 TeV) collisions at the Large Hadron Collider (LHC), we deploy a previously developed transport approach for charmonia and bottomonia to evaluate the kinetics of 𝐵 𝑐 mesons throughout the fireball formed in these reactions. The main inputs to our approach are two transport parameters: the 𝐵 𝑐 's reaction rate and equilibrium limit. Both quantities are determined by previous calculations via a combination of charm and bottom sectors. In-medium binding energies of 𝐵 𝑐 mesons are calculated from a thermodynamic 𝑇 matrix with a lattice-QCD constrained potential, and figure in their inelastic reaction rates. Temperature-dependent equilibrium limits include charm- and bottom-quark fugacities based on their initial production. We compute the centrality dependence of inclusive 𝐵 𝑐 production and transverse-momentum (𝑝 𝑇 ) spectra using two different recombination models: instantaneous coalescence and resonance recombination. The main uncertainty in the resulting nuclear modification factors, 𝑅 AA , is currently associated with the 𝐵 𝑐 cross section in elementary 𝑝⁢𝑝 collisions, caused by the uncertainty in the branching ratio for the 𝐵$^{−}_{𝑐}$ → 𝐽/𝜓⁢𝜇 − ⁢$\overline{𝜈}$ decay. Our results indicate a large enhancement of the 𝑅 AA at low 𝑝 𝑇 , with significant regeneration contributions up to 𝑝 𝑇 ≃ 20GeV. Comparisons to CMS data are carried out but firm conclusions will require a more accurate value of the branching ratio, or alternative channels to measure the 𝐵 𝑐 production in 𝑝⁢ collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Framework for phase transitions between the Maxwell and Gibbs constructions at finite temperature

The characteristics of the hadron-to-quark first-order phase transition differ depending on whether charge neutrality is locally or globally fulfilled. In 𝛽-equilibrated matter, these two possibilities correspond to the Maxwell and Gibbs constructions. Recently, we presented a new framework in which a continuously varying parameter allows one to describe a first-order phase transition in intermediate scenarios to the two extremes of fully local and fully global charge neutrality. In this work, we extend the previous framework to finite temperatures and out-of-𝛽 equilibrium conditions, making it available for simulations of core-collapse supernovae and binary neutron star mergers. We investigate its impact on key thermodynamic quantities across a range of baryon densities, temperatures, and electron fractions. We find that when matter is not in 𝛽 equilibrium, the pressure in the mixed phase is not constant even for the case of fully local charge neutrality. Moreover, we compute the thermal index using three different approaches, demonstrating that the finite-temperature extension of an equation of state using a constant thermal index can be ill defined when applied to the mixed phase.

QCD phase transitions↗
DOE OSTI2021

Watch and learn—a generalized approach for transferrable learning in deep neural networks via physical principles

Transfer learning refers to the use of knowledge gained while solving a machine learning task and applying it to the solution of a closely related problem. Such an approach has enabled scientific breakthroughs in computer vision and natural language processing where the weights learned in state-of-the-art models can be used to initialize models for other tasks which dramatically improve their performance and save computational time. Here we demonstrate an unsupervised learning approach augmented with basic physical principles that achieves fully transferrable learning for problems in statistical physics across different physical regimes. By coupling a sequence model based on a recurrent neural network to an extensive deep neural network, we are able to learn the equilibrium probability distributions and inter-particle interaction models of classical statistical mechanical systems. Our approach, distribution-consistent learning, DCL, is a general strategy that works for a variety of canonical statistical mechanical models (Ising and Potts) as well as disordered interaction potentials. Using data collected from a single set of observation conditions, DCL successfully extrapolates across all temperatures, thermodynamic phases, and can be applied to different length-scales. This constitutes a fully transferrable physics-based learning in a generalizable approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Atomic structure considerations for the low-temperature opacity of Xenon

Here, we have begun a preliminary investigation into the opacity of Xe at low temperatures (<50 eV). The emissivity and opacity of Xe is a crucial factor in determining the utility of Xe in EUV lithography, with numerous industrial applications. To this end, we have been exploring the accuracy of some approximations used in opacity models for the relevant ion stages of Xe (6 times ionized through 20 times ionized). Due to the complexity of Xe atomic structure, one needs to use full configuration-interaction to properly describe the strong mixing in the various n=4 sub-shells that give rise to the Δn=0 and Δn=1 transitions that dominate the opacity spectrum at low temperatures. Since calculations that include full configuration-interaction for large numbers of configurations quickly become computationally prohibitive, we have explored hybrid calculations, in which full configuration-interaction is retained for the most important transitions, while intermediate coupling is employed for all other transitions. After verifying atomic structure properties, local-thermodynamic-equilibrium (LTE) opacities are generated using the ATOMIC code at selected temperatures and densities and compared to various experimental results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

SAP-X2C: Optimally-Simple Two-Component Relativistic Hamiltonian with Size-Intensive Picture Change

We present a simple relativistic exact 2-component (X2C) Hamiltonian that models two-electron picture-change effects using Lehtola’s superposition of atomic potentials (SAP) [S. Lehtola, J. Chem. Theory Comput. 15, 1593−1604 (2019)]. The SAP-X2C approach retains the low cost and technical simplicity of the popular 1-electron X2C (1eX2C) predecessor but is significantly more accurate and has a well-defined thermodynamic limit, making it applicable to extended systems (such as large molecules and periodic crystals). The assessment of the SAP-X2C-based Hartree−Fock total and spinor energies, spin−orbit splittings, equilibrium bond distances, and harmonic vibrational frequencies suggests that SAP-X2C is similar to the more complex atomic meanfield (AMF) X2C counterparts in its ability to approximate the 4-component Dirac−Hartree−Fock reference.

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

“Magic” Orientation Angles to Suppress Spin-Driven Hall Currents in Anisotropic 2D Materials with an Ideal Skyrmion Gas

The Hall effect depending on conduction electron spin projection becomes very different for anisotropic 2D crystals. In this case the spin-dependent electron current strongly determined by the orientation angle, θ, of the sample with respect to an applied electric eld. The spin-up and -down components of the direct and Hall charge currents oscillate with the angle 2&theta;. The direct and Hall components of the current have the structure where there are the angle-independent part, oscillation amplitude, and phase shift. All three quantitates strongly depend on the electron spin projection, electron mass ratio, and skyrmion size. We nd that there are "magic" orientation angles where the spin-up, spin-down, and total Hall currents vanish. There is a great interest in computations based on 2D materials with skyrmions. Such properties can be useful for computer logic operations based on skyrmions.

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

An Atomistic Study of Reactivity in Solid-State Electrolyte Interphase Formation for Li/Li7P3S11

Lithium metal batteries offer superior volumetric and gravimetric specific capacities compared to those based on traditional graphite anodes. Although advancements in solid-state electrolytes address safety concerns, challenges remain, particularly regarding interphase formation in lithium metal anodes. This work presents a computational framework based on high-throughput first-principles density functional theory and machine-learning interatomic potentials (MLIPs) including automated iterative, active learning to enable robust computational exploration of interphase formation between lithium metal anodes and an inorganic solid-state electrolyte. As a demonstration, we apply the framework to a Li/Li7P3S11 interface and find that it accurately identifies the experimentally observed, thermodynamically stable interphase products as well as their overall spatial arrangement within a heterogeneous, amorphous layered structure, with Li2S domains of nanocrystallinity. Our simulations show two stages, a fast and slow diffusion reaction regime, that corroborate the relative phase formation rate of Li x P, Li2S, and Li3P. Using the Onsager transport theory, we capture time-dependent ionic diffusion within the reacting interface, including cross-correlation effects. We found that cross-correlation effects between Li-P and P-S ionic motion significantly influence P-ion diffusion, making it highly sensitive to the local environment and potentially leading to "kinetic trapping" of Li-P phases. The passivation of the interface is shown as the ionic fluxes all approach zero, effectively halting interphase growth.

Diffusion↗
DOE OSTI2025

Theoretical and experimental quantification of Suzuki segregation enthalpy and strengthening mechanisms in a binary alloy

Solute segregation to planar defects in metallic alloys has been shown to drastically alter mechanical properties. While various works using first-principles and thermodynamic calculations have studied the fundamental driving forces for solute segregation via the Suzuki criterion, planar defect energy, or a comparison of energies of the HCP-like phase and FCC matrix, a quantitative experimental and computational comparison of equilibrium composition and segregation enthalpies has not yet been reported. In this work, we predict the equilibrium composition and segregation enthalpy to intrinsic stacking faults in a Ni-60Co (at.%) alloy and compare the results to two independent experimental methods. We observed that Co segregates to the innermost two planes of the intrinsic stacking fault, and we found that the experimental segregation enrichment, measured from transmission electron microscopy energy dispersive X-ray spectroscopy, of the faults is 6.8 at.% Co, which is 2.2 at.% less than the predicted value at the same temperature. We also find that the segregation enthalpy measured from the composition profile is −21.1 ± 6.4 meV/atom and separately from differential scanning calorimetry segregation enthalpy is −33.2 meV/atom, whereas the predicted enthalpy is −31 ± 1 meV/atom. Based on these results, we determine that segregation occurs very rapidly, within 8 min at temperatures as low as 36% of the homologous solidus temperature. Furthermore, this analysis provides an overview of the possible dislocation mechanisms responsible for strengthening effects due to solute segregation, and concludes that changes in room temperature hardness from local phase transformation is likely tied to post-segregation room temperature equilibrium partial separation distance.

Ab initio calculation↗
DOE OSTI2024

Wave steepening and shock formation in ultracold neutral plasmas

We present observations of wave steepening and signatures of shock formation during expansion of ultracold neutral plasmas formed with an initial density distribution that is centrally peaked and decays exponentially with distance. The plasma acceleration and velocity decrease at large distance from the plasma center, leading to central ions overtaking ions in the outer regions and the development of a steepening front that is narrow compared to the size of the plasma. The density and velocity change dramatically across the front, and significant heating of the ions is observed in the region of steepest gradients. For a reasonable estimate of electron temperature, the relative velocity of ions on either side of the front modestly exceeds the local sound speed (Mach number M≳1). This indicates that by sculpting steep density gradients, it is possible to create the conditions for shock formation, or very close to it, opening a new avenue of research for ultracold neutral plasmas.

Computational fluid dynamics↗
DOE OSTI2021

Computational discovery of stable and metastable ternary oxynitrides

Materials design from first principles enables exploration of uncharted chemical spaces. Extensive computational searches have been performed for mixed-cation ternary compounds, but mixed-anion systems are gaining increased interest as well. Central to computational discovery is the crystal structure prediction, where the trade-off between reliance on prototype structures and size limitations of unconstrained sampling has to be navigated. We approach this challenge by letting two complementary structure sampling approaches compete. We use the kinetically limited minimization approach for high-throughput unconstrained crystal structure prediction in smaller cells up to 21 atoms. On the other hand, ternary—and, more generally, multinary—systems often assume structures formed by atomic ordering on a lattice derived from a binary parent structure. Thus, we additionally sample atomic configurations on prototype lattices with cells up to 56 atoms. Using this approach, we searched 65 different charge-balanced oxide–nitride stoichiometries, including six known systems as the control sample. The convex hull analysis is performed both for the thermodynamic limit and for the case of synthesis with activated nitrogen sources. We identified 34 phases that are either on the convex hull or within a viable energy window for potentially metastable phases. We further performed structure sampling for “missing” binary nitrides whose energies are needed for the convex hull analysis. Among these, we discovered metastable Ce3N4 as a nitride analog of the tetravalent cerium oxide, which becomes stable under slightly activated nitrogen condition µN > +0.07 eV. Given the outsize role of CeO2 in research and application, Ce3N4 is a potentially important discovery.

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

Adaptive simulations enable computational design of electron beam processing of nanomaterials with supersonic micro-jet precursor

Focused Electron Beam Induced Processing (FEBIP) is a powerful tool for the “direct-write” of nanomaterials with the possibility of atomistic control on suspended 2D material substrates. FEBIP capabilities have been significantly expanded by using a localized jet-based delivery of precursors, and especially when a thermally energized supersonic micro-jet enhances the delivery of mass flux along with controlling the far-from-equilibrium thermodynamic state of adsorbed adatoms. The possibilities of growing nanomaterials with “dialed-in” composition with ultra-high growth rates and aspect ratios up to 100:1 have been demonstrated using the supersonic micro-jet-FEBIP. Bringing this scientific discovery to the level of maturity required for practical applications in additive nanomanufacturing requires simulation tools that are capable of capturing the complex flow physics of micro-jet-substrate interactions that bridge a wide range of flow regimes from the high-density gas micro-jet expanding into the vacuum environment of FEBIP. To address this significant computational challenge, a new approach has been developed and described in this work for a multiscale adaptive DSMC (Direct Simulation Monte Carlo) algorithm to predict the micro-jet gas dynamics in an FEBIP environment, spanning the full range of flow regimes from low Knudsen (Kn) number O(0.01) continuum flow to high Kn of O(10) for the molecular flow within a unified computational framework. Here, the fundamental principles of the adaptive DSMC algorithm are described, its viability as a DSMC technique is demonstrated, and computational improvements for the benchmark cases are discussed in the context of advancing 3D nanofabrication with FEBIP. Ultimately, combining the first principle simulations via adaptive DSMC with the complementary experimental data will enable the creation of the powerful CAD tools for in silico design and optimal operation of micro-jet-FEBIP.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Degrees of rate control and AutoDiff-driven direct sensitivity analysis in heterogeneous catalysis

Despite the wide application and benefits of the degree of rate control (DRC) analysis, several details remain argued, particularly about the conservation of DRCs at transient (TR) and steady-state (SS) conditions, especially for complex reaction networks. This work argues that previous proofs about the conservation properties of DRCs have been incomplete, and we provide new mathematical proofs at TR and SS conditions. In addition, we use both analytical (automatic differentiation) and numerical (finite difference) approaches to compute DRCs for the case study of ethane hydrogenolysis (EH) over Pt(111). This work confirms that at both TR and SS conditions, the sum of all DRCs, i.e., sum of the degrees of kinetic (DKRC) and thermodynamic rate control (DTRC), is conserved at zero. At SS conditions, the sum of DKRC is conserved at 1 while the sum of DTRC is conserved at −1. In corroboration of previous works, we show that the DTRC for any adsorbate at SS is equal to the product of the species coverage and a constant. In contrast, at TR conditions, the individual sums of both DTRC and DKRC are not conserved and can be any real number, with potential implications for the novel field of dynamic catalysis. Finally, we show that the conventional finite difference (FD) approach, only useful at SS, is prone to inaccuracy and very sensitive to the value of the differential change applied. The optimal differential value also varies significantly with system and rate definition. Consequently, we describe and illustrate in this work the application of the automatic differentiation (AD) approach for the more accurate determination of DRCs at both TR and SS conditions.

Automatic differentiation↗
DOE OSTI2023

From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstraction by $CH_3O\dot{O}$ radicals

Hydrogen atom abstraction by methyl peroxy ($CH_3O\dot{O}$) radicals can play an important role in gasoline/ethanol interacting chemistry for fuels that produce high concentrations of methyl radicals. Detailed kinetic reactions for hydrogen atom abstraction by $CH_3O\dot{O}$ radicals from the components of FGF-LLNL (a gasoline surrogate) including cyclopentane, toluene, 1-hexene, n -heptane, and isooctane have been systematically studied in this work. Here, the M06-2X/6-311 ++ G(d,p) level of theory was used to obtain the optimized structure and vibrational frequency for all stationary points and the low-frequency torsional modes. The 1-D hindered rotor treatment for low-frequency torsional modes was treated at M06-2X/6-31G level of theory. The UCCSD(T)-F12a/cc-pVDZ-F12 and QCISD(T)/CBS level of theory were used to calculate single point energies for all species. High pressure limiting rate constants for all hydrogen atom abstraction channels were performed using conventional transition state theory with unsymmetric tunneling corrections. Individual rate constants are reported in the temperature range from 298.15 to 2000 K. Our computed results show that the abstraction of allylic hydrogen atoms from 1-hexene is the fastest at low temperatures. When the temperature increases, the hydrogen atom abstraction reaction channel at the primary alkyl site gradually becomes dominant. Thermodynamics properties for all stable species and high-pressure limiting rate constants for each reaction pathway obtained in this work were incorporated into the latest gasoline surrogate/ethanol model to investigate the influence of the rate constants calculated here on model predicted ignition delay times.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Thermodynamically consistent incorporation of the Langmuir adsorption model into compressible fluctuating hydrodynamics

For a gas–solid interfacial system where chemical species undergo reversible adsorption, we develop a mesoscopic stochastic modeling method that simulates both gas-phase hydrodynamics and surface coverage dynamics by coupling the Langmuir adsorption model with compressible fluctuating hydrodynamics. To this end, we derive a thermodynamically consistent mass–energy update scheme that accounts for how the mass and energy variables in the gas and surface subsystems should be updated according to the changes in the number of molecules of each species in each subsystem due to adsorption and desorption events. By performing a stochastic analysis for the ideal Langmuir model and the full hydrodynamic system, we analytically confirm that our mass–energy update scheme captures thermodynamic equilibrium predicted by equilibrium statistical mechanics. We find that an internal energy correction term is needed, which is attributed to the difference in the mean kinetic energy of gas molecules colliding with the surface from that computed from the Maxwell–Boltzmann distribution. By performing an equilibrium simulation study for an ideal gas mixture of CO and Ar, with CO undergoing reversible adsorption, we validate our overall simulation method and implementation.

Adsorption↗
DOE OSTI2022

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗