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

High-throughput calculations of charged point defect properties with semi-local density functional theory—performance benchmarks for materials screening applications

Abstract Calculations of point defect energetics with Density Functional Theory (DFT) can provide valuable insight into several optoelectronic, thermodynamic, and kinetic properties. These calculations commonly use methods ranging from semi-local functionals with a-posteriori corrections to more computationally intensive hybrid functional approaches. For applications of DFT-based high-throughput computation for data-driven materials discovery, point defect properties are of interest, yet are currently excluded from available materials databases. This work presents a benchmark analysis of automated, semi-local point defect calculations with a-posteriori corrections, compared to 245 “gold standard” hybrid calculations previously published. We consider three different a-posteriori correction sets implemented in an automated workflow, and evaluate the qualitative and quantitative differences among four different categories of defect information: thermodynamic transition levels, formation energies, Fermi levels, and dopability limits. We highlight qualitative information that can be extracted from high-throughput calculations based on semi-local DFT methods, while also demonstrating the limits of quantitative accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Statistical perspective on embrittling potency for intergranular fracture

Embrittling potency is a thermodynamic metric that assesses the influence of solute segregation to a grain boundary (GB) on intergranular fracture. Historically, authors of studies have reported embrittling potency as a single scalar value, assuming a single segregation site of importance at a GB and a particular cleavage plane. However, the topography of intergranular fracture surfaces is not generally known a priori. Accordingly, we, in this study, present a statistical ensemble approach to compute embrittling potency, where many free surface (FS) permutations are systematically considered to model fracture of a GB. The result is a statistical description of the thermodynamics of GB embrittlement. As a specific example, embrittling potency distributions are presented for Cr segregation to sites at two Ni $\langle 111 \rangle$ symmetric tilt GBs using atomistic simulations. We show that the average embrittling potency for a particular GB site, considering an ensemble of FS permutations, is not equal to the embrittling potency computed using the lowest energy pair of FSs. A mean GB embrittlement is proposed, considering both the likelihood of formation of a particular FS and the probability of solute occupancy at each GB site, to compare the relative embrittling behavior of two distinct GBs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Destabilizing high-capacity high entropy hydrides via earth abundant substitutions: From predictions to experimental validation

The vast chemical space of high entropy alloys (HEAs) makes trial-and-error experimental approaches for materials discovery intractable and often necessitates data-driven and/or first principles computational insights to successfully target materials with desired properties. In the context of materials discovery for hydrogen storage applications, a theoretical prediction-experimental validation approach can vastly accelerate the search for substitution strategies to destabilize high-capacity hydrides based on benchmark HEAs, e.g. TiVNbCr alloys. Here, in this study, machine learning predictions, corroborated by density functional theory calculations, predict substantial hydride destabilization with increasing substitution of earth-abundant Fe content in the (TiVNb) 75 Cr 25-x Fe x system. The as-prepared alloys crystallize in a single-phase bcc lattice for limited Fe content x < 7, while larger Fe content favors the formation of a secondary C14 Laves phase intermetallic. Short range order for alloys with x < 7 can be well described by a random distribution of atoms within the bcc lattice without lattice distortion. Hydrogen absorption experiments performed on selected alloys validate the predicted thermodynamic destabilization of the corresponding fcc hydrides and demonstrate promising lifecycle performance through reversible absorption/desorption. This demonstrates the potential of computationally expedited hydride discovery and points to further opportunities for optimizing bcc alloy ↔ fcc hydrides for practical hydrogen storage applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A guided ion beam investigation of UO 2 + thermodynamics and f orbital participation: Reactions of U + + CO 2 , UO + + O 2 , and UO + + CO

A guided ion beam tandem mass spectrometer was employed to study the reactions of U + + CO 2 , UO + + O 2 , and the reverse of the former, UO + + CO. Reaction cross sections as a function of kinetic energy over about a three order of magnitude range were studied for all systems. The reaction of U + + CO 2 proceeds to form UO + + CO with an efficiency of 118% ± 24% as well as generating UO 2 + + C and UCO + + O. The reaction of UO + + O 2 forms UO 2 + in an exothermic, barrierless process and also results in the collision-induced dissociation of UO + to yield U + . In the UO + + CO reaction, the formation of UO 2 + in an endothermic process is the dominant reaction, but minor products of UCO + + O and U + + (O + CO) are also observed. Analysis of the kinetic energy dependences observed provides the bond energies, D 0 (U + –O) = 7.98 ± 0.22 and 8.05 ± 0.14 eV, D 0 (U + –CO) = 0.73 ± 0.13 eV, and D 0 (OU + –O) = 7.56 ± 0.12 eV. The values obtained for D 0 (U + –O) and D 0 (OU + –O) agree well with the previously reported literature values. To our knowledge, this is the first experimental measurement of D 0 (U + –CO). Furthermore, an analysis of the oxide bond energies shows that participation of 5f orbitals leads to a substantial increase in the thermodynamic stability of UO 2 + relative to ThO 2 + and especially transition metal dioxide cations.

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

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Thermal Extraction and Delivery from Advanced Reactor Systems

This document presents preliminary thermodynamic analysis of thermal extraction and delivery from advanced nuclear reactors. Light-water, high temperature gas, sodium fast, and molten salt reactors are templated and analyzed for thermal extraction from nominal steam production conditions. Uncontrolled, controlled, and partial admission steam extraction methods are described and analyzed to show the differences in system impact from the various methods. Thermal delivery modeling results show comparisons that have become possible through the computational capabilities of the thermal delivery model. Due to the significant degrees of freedom of the problem, specific cases are used to compare the thermodynamic costs of heat transport using various heat transfer fluids. All this work is done using Modelica models, allowing for the future application of models in dynamic analyses beyond the steady-state calculations presented in this report. Appendices describe additional capabilities and methods that can help users analyze complex problems using HYBRID models developed in Modelica.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A Simple Criterion for Feasibility of Heat Integration between Distillation Streams Based on Relative Volatilities

In a multicomponent distillation configuration, there are numerous sources and sinks of heat, and a potential way to reduce the heat duty requirement is to perform heat integration. Unfortunately, an algorithmic search of the optimal heat integration opportunities is intractable when the required temperatures of intermediate streams are computed via complex models. Instead, in this work, we introduce pressure-scaled pseudo relative volatility, a new metric to compare stream temperatures. Here, we justify the use of pseudo relative volatility by proving that this variable is a monotonically increasing function of the liquid fraction in a saturated mixture stream. Using this metric, we derive a shortcut criterion to check the feasibility of various heat integration opportunities, such as thermal coupling via heat transfer (TCH). The advantage of this approach is that it circumvents the need for explicit temperatures and instead relies on composition, component relative volatilities, and pressure—quantities that are readily available in shortcut models for optimization of distillation configurations. Leveraging this fact, we propose a new optimization framework that identifies feasible TCHs that we consider within the formulation while minimizing the total heat duty of a distillation configuration. We demonstrate, on a few examples, that our formulation can identify heat duty efficient configurations, some of which are multieffect configurations. Using this methodology, we discover configurations that are not only simpler than the fully thermally coupled (FTC) configuration but also have a much lower heat duty.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Molecular Characterization of Organosulfate-Dominated Aerosols over Agricultural Fields from the Southern Great Plains by High-Resolution Mass Spectrometry

The molecular composition of organic aerosols, especially for day/nighttime variations of organosulfates above agricultural fields, is not well understood despite profound impacts on regional climate, crop production, air quality, and human health. Here, nanospray desorption electrospray ionization with high-resolution mass spectrometry (nano-DESI-HRMS) is used to interrogate the molecular composition of organic aerosols collected at the Southern Great Plains, located in an agricultural region of Oklahoma. Identified molecular formulae featured carbon, hydrogen, oxygen (CHO), nitrogen (CHNO), and/or sulfur (CHOS, CHNOS), with higher organosulfate proportions during daytime (41%) compared to nighttime (30%). In this work, nighttime aerosols featured increases in CHO, CHNO, and extremely low volatility organic carbon (ELVOC) species. However, due to high relative humidity, the nighttime aerosols phase state was found to be more liquid-like than daytime aerosols using parametrized glass transition temperatures. Aerosol molecular composition from an anthropogenically influenced plume (southerly winds) showed significant increases in CHOS and ELVOC species. By comparison with chamber studies, CHOS species are suspected to be of mixed biogenic and anthropogenic origin, whereas CHNOS species (not identified in the southerly winds) are suggested to predominately be of biogenic origin. Overall, this study provides key insight into organosulfates above agricultural fields, demonstrating dependence upon day/night cycles and episodic anthropogenic emissions.

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

Mechanistic Insights into the Conversion of Biorenewable Levoglucosanol to Dideoxysugars

Here, a molecular understanding of the conversion of biorenewable threo- and erythro-levoglucosanol (LGOL) to 3,4-dideoxysugars in aqueous medium is provided based on first-principles simulations. The synthetic importance of this transformation is that these intermediates can be quantitatively hydrogenated to (S,S)/(S,R) hexane-1,2,5,6-tetrol (tetrol), whose stereochemistry depends on which dideoxy sugar intermediates are formed during LGOL conversion. The thermodynamic and kinetic feasibility of the acetal (R 2 C(OR) 2 ) hydrolysis in LGOL is investigated via computing the free energy profile. In aqueous medium, the rate-determining step of LGOL hydrolysis is the protonation of the anhydro-bridge oxygen atom of LGOL concurrent with ring opening, yielding the cyclic forms of 3,4-dideoxymannose (DDM) and 3,4-dideoxyglucose (DDG) from threo- and erythro-LGOL, respectively. The measured activation energies of LGOL hydrolysis are 20.5 and 23.6 kcal/mol for DDM and DDG formation, respectively. These values are in agreement with the computed protonation free energies of 17.1 and 18.2 kcal/mol, respectively. Based on the simulations, a Bronsted base-catalyzed isomerization from DDG or DDM to 3,4-dideoxy fructose (DDF) is preferred with lower apparent activation free energy barriers compared to the acid-catalyzed isomerization. In summary, this study provides mechanistic information about the conversion of the biomass-derived anhydro-sugar LGOL to 3,4-dideoxy sugars, which are precursors to renewable high-value chemicals.

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

LoTToR: An Algorithm for Missing-Wedge Correction of the Low-Tilt Tomographic 3D Reconstruction of a Single-Molecule Structure

A single-molecule three-dimensional (3D) structure is essential for understanding the thermal vibrations and dynamics as well as the conformational changes during the chemical reaction of macromolecules. Individual-particle electron tomography (IPET) is an approach for obtaining a snap-shot 3D structure of an individual macromolecule particle by aligning the tilt series of electron tomographic (ET) images of a targeted particle through a focused iterative 3D reconstruction method. The method can reduce the influence on the 3D reconstruction from large-scale image distortion and deformation. Due to the mechanical tilt limitation, 3D reconstruction often contains missing-wedge artifacts, presented as elongation and an anisotropic resolution. Here, we report a post-processing method to correct the missing-wedge artifact. This low-tilt tomographic reconstruction (LoTToR) method contains a model-free iteration process under a set of constraints in real and reciprocal spaces. A proof of concept is conducted by using the LoTToR on a phantom, i.e., a simulated 3D reconstruction from a low-tilt series of images, including that within a tilt range of ±15°. The method is validated by using both negative-staining (NS) and cryo-electron tomography (cryo-ET) experimental data. A significantly reduced missing-wedge artifact verifies the capability of LoTToR, suggesting a new tool to support the future study of macromolecular dynamics, fluctuation and chemical activity from the viewpoint of single-molecule 3D structure determination.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A framework to evaluate machine learning crystal stability predictions

The rapid adoption of machine learning in various scientific domains calls for the development of best practices and community agreed-upon benchmarking tasks and metrics. We present Matbench Discovery as an example evaluation framework for machine learning energy models, here applied as pre-filters to first-principles computed data in a high-throughput search for stable inorganic crystals. We address the disconnect between (1) thermodynamic stability and formation energy and (2) retrospective and prospective benchmarking for materials discovery. Alongside this paper, we publish a Python package to aid with future model submissions and a growing online leaderboard with adaptive user-defined weighting of various performance metrics allowing researchers to prioritize the metrics they value most. To answer the question of which machine learning methodology performs best at materials discovery, our initial release includes random forests, graph neural networks, one-shot predictors, iterative Bayesian optimizers and universal interatomic potentials. We highlight a misalignment between commonly used regression metrics and more task-relevant classification metrics for materials discovery. Accurate regressors are susceptible to unexpectedly high false-positive rates if those accurate predictions lie close to the decision boundary at 0 eV per atom above the convex hull. The benchmark results demonstrate that universal interatomic potentials have advanced sufficiently to effectively and cheaply pre-screen thermodynamic stable hypothetical materials in future expansions of high-throughput materials databases.

Riebesell, Janosh↗
DOE OSTI2026

Modeling the behavior of concentrated aqueous HNO 3 using machine learning interatomic potentials

We develop two multi-defect machine learning interatomic potentials (MLIPs) trained at the BLYP-D2 and PBE-D3 density functional theories using the DeepMD-kit, allowing for the investigation of structural and thermodynamic properties of nitric acid over a wide range of concentrations via molecular dynamics (MD) simulations. We directly compute the degree of dissociation, α, and pK a from MD simulations, revealing that HNO 3 behaves as a weaker acid at higher concentrations, noting that our standard-state pK a value is in excellent agreement with the experimental one. In general, good agreement is observed with experimental results such as α and density outside the training dataset, with only modest deviations at low-to-medium concentrations. We benchmark our custom multi-defect DeepMD MLIPs against foundational models MACE-MP0 and MACE-OFF23. The foundation models capture some aspects of HNO 3 /NO 3 − solvation in concentrated nitric acid but show noticeable density errors and miss subtle structural features relevant to spectroscopy, whereas the bespoke DeepMD MLIPs yield more compact solvation shells, reproduce density-concentration trends, and run ∼12–15× faster than MACE-MP0. Although classical FFs are still more efficient and match experimental densities better, they lack chemical reactivity and thus cannot predict α or pK a , underscoring the need for system-specific reactive MLIPs beyond universal MLIPs.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2025

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

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

pnnl/Pathway_Controlled_Optimization

A method for prediction of cell regulation to maximize cell growth rate while maintaining the solvent capacity of the cell. Prediction is formulated as an optimization problem using a thermodynamic framework that can leverage experimental data. We develop a formulation and variable initialization procedure that allows for computing solutions of the optimization with an interior point method

Central, PNNL Developer↗
DOE OSTI2021

Cooperative Research and Development Agreement among National Energy Technology Laboratory, Advanced Manufacturing LLC, and Connecticut Center for Advanced Technology [Abstract]

The National Energy Technology Laboratory (NETL), Advanced Manufacturing LLC (AMLLC) and Connecticut Center for Advanced Technology (CCAT) will collaborate to develop and manufacture cost-effective, oxide dispersion-strengthened, NiCrFeCo-rich high entropy alloys that are superior to Nickel-based superalloys (e.g. IN740) for repair or replacement service in extreme environments. The new alloys can result in higher efficiency and reduced cost of advanced power systems such as advanced ultra-supercritical (AUSC) boilers and gas turbines. High throughput, multiscale computational modeling will be performed to accelerate alloy discovery by interrogating the intrinsic properties of the alloys including thermodynamic, kinetic, and mechanical properties. Validation and evaluation will be performed, including additive manufacturing, microstructure characterization, and mechanical properties tests.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An approach to learn regulation to maximize growth and entropy production rates in metabolism

Elucidating cell regulation remains a challenging task due to the complexity of metabolism and the difficulty of experimental measurements. Here we present a method for prediction of cell regulation to maximize cell growth rate while maintaining the solvent capacity of the cell. Prediction is formulated as an optimization problem using a thermodynamic framework that can leverage experimental data. We develop a formulation and variable initialization procedure that allows for computing solutions of the optimization with an interior point method. The approach is applied to photoheterotrophic growth of Rhodospirilium rubrum using ethanol as a carbon source, which has applications to biosynthesis of ethylene production. Growth is captured as the rate of synthesis of amino acids into proteins, and synthesis of nucleotide triphoshaptes into RNA and DNA. The method predicts regulation that produces a high rate of protein and RNA synthesis while DNA synthesis is reduced close to zero in agreement with production of DNA being turned off for much of the cell cycle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Divergent carbon use efficiency-growth rate tradeoff in popular biological growth models

Carbon use efficiency (CUE) is an important trait emerging from processes regulating biological growth. CUE can be computed either based on the growth of structural biomass or total biomass divided by substrate uptake rate. Nonequilibrium thermodynamics and observations suggest that, for an exponentially growing population of cells, structural biomass CUE should first increase, then peak, and finally decrease with specific growth rate; meanwhile, total biomass CUE increases asymptotically with specific growth rate. We compared predictions from six popular models that are often used for plant and microbial growth in existing ecosystem models. We found that, for an exponentially growing population of biological cells, (1) the source-driven Pirt and Compromise models predict that structural biomass CUE increase asymptotically with growth rate; (2) the apparent sink-driven modified Droop model predicts that structural biomass CUE decreases with growth rate; and (3) the sink-driven variable internal storage model and two dynamic energy budget models predict that structural biomass CUE first increases, then peaks, and finally decreases with growth rate. Moreover, the modified Droop model predicts that total biomass CUE is constant with growth rate, while all other five models predict that total biomass CUE increases with growth rate asymptotically. For non-exponential biological growth, we show that there is no static relationship between total biomass CUE or structural biomass CUE with respect to either growth rate or temperature. Therefore, we contend that biological growth models should explicitly represent interactions between substrate acquisition, substate transformation, and maintenance respiration to better capture observed CUE dynamics, and the sink-driven model should be preferred for general ecosystem biogeochemistry modeling.

Tang, Jinyun [Lawrence Berkeley National Laborator↗
DOE OSTI2026

RidgeAlloy: A high pressure die casting alloy that captures the approaching wave of automotive body sheet scrap

Oak Ridge National Laboratory (ORNL) has recently launched an effort to economically expand the North American aluminum supply chain by designing and developing a new family of Al-Mg-Si-Fe-Mn structural die cast aluminum alloys that can be made from up to 100% mixed 5xxx and 6xxx series automotive post-consumer sheet scrap and do not require heat treatment. The goal is to develop an alloy capable of capturing the approaching scrap wave of aluminum sheet intensive vehicles (in the 2030’s) into integrated structural high pressure die castings, rather than down-cycling these high-grade sheet alloys into non-structural castings. ORNL discovered a significant gap in the commercial thermodynamic databases that were unable to predict the trends of a key primary intermetallic phase as a function of composition. Experiments were performed to correct this gap and to create a unique, corrected quaternary database for Al-Mg-Si-Fe. High throughput computational methods, including solidification simulations, were used to identify a composition volume capable of avoiding this embrittling primary intermetallic phase, even at higher Fe + Si contents. Laboratory scale castings validated the thermodynamic and solidification predictions and resulted in lab-cast alloys with little primary intermetallic, and which met or exceeded structural casting properties requirements for yield strength and ductility. An accelerated demonstration of an HPDC automotive part was conducted in collaboration with two U.S. small businesses. Ingots made from 100% recycled 5xxx and 6xxx series scrap (plus Fe) were used to die cast one variant of this new family of Al-Mg-Si-Fe-Mn alloys (with high Si + Fe content) into a mid-sized HPDC structural-type automotive component. The HPDC alloy demonstration part showed good properties (especially ductility) and good castability for a complex component.

Plotkowski, Alex [ORNL] (ORCID:0000000154718681)↗