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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dynamics of disordered mechanical systems with large connectivity, free probability theory, and quasi-Hermitian random matrices

Disordered mechanical systems with high connectivity represent a limit opposite to the more familiar case of disordered crystals. Individual ions in a crystal are subjected essentially to nearest-neighbor interactions. In contrast, the systems studied in this paper have all their degrees of freedom coupled to each other. Thus, the problem of linearized small oscillations of such systems involves two full positive-definite and non-commuting matrices, as opposed to the sparse matrices associated with disordered crystals. Consequently, the familiar methods for determining the averaged vibrational spectra of disordered crystals, introduced many years ago by Dyson and Schmidt, are inapplicable for highly connected disordered systems. In this paper we apply random matrix theory (RMT) to calculate the averaged vibrational spectra of such systems, in the limit of infinitely large system size. At the heart of our analysis lies a calculation of the average spectrum of the product of two positive definite random matrices by means of free probability theory techniques. We also show that this problem is intimately related with quasi-hermitian random matrix theory (QHRMT), which means that the ‘hamiltonian’ matrix is hermitian with respect to a non-trivial metric. This extends ordinary hermitian matrices, for which the metric is simply the unit matrix. The analytical results we obtain for the spectrum agree well with our numerical results. The latter also exhibit oscillations at the high-frequency band edge, which fit well the Airy kernel pattern. We also compute inverse participation ratios of the corresponding amplitude eigenvectors and demonstrate that they are all extended, in contrast with conventional disordered crystals. Finally, we compute the thermodynamic properties of the system from its spectrum of vibrations. In addition to matrix model analysis, we also study the vibrational spectra of various multi-segmented disordered pendula, as concrete realizations of highly connected mechanical systems. A universal feature of the density of vibration modes, common to both pendula and the matrix model, is that it tends to a non-zero constant at vanishing frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Corrosion of 316 L stainless steel under the natural circulation of molten NaCl-MgCl 2 salt

The corrosion behavior of 316 L stainless steel (SS) was studied via the natural circulation of molten eutectic NaCl-MgCl 2 salt through a microloop. The post-corrosion tested 316 L SS microloop sections were characterized with microscopy techniques to determine the microstructural and microchemical changes that occurred at the alloy/salt interface. It was found that 316 L SS showed heterogeneous dissolution at the hot-leg, whereas deposition of corrosion products occurred at the cold-leg. For the first time, experimentally obtained molten salt flow-induced corrosion of 316 L SS results were combined with computational thermodynamic-kinetic models to validate the dissolution and deposition in terms of elemental compositional changes at the alloy/salt interface. The thermodynamic-kinetic modeling predicted that the heterogeneous dissolution of Cr and Fe from the hot-leg section of 316 L SS persisted throughout the salt circulation. The model also estimated that, despite Cr deposition starting earlier than Fe, the total redeposition of Fe is expected to be significantly greater than that of Cr over the circulation of salt. Furthermore, the modeling accurately predicted the subsurface enrichment of Mo which is attributed to the reduced Cr activity and the relatively higher diffusion rate of Mo within the alloy matrix. Here, the agreement between modeling and experimental results confirms that Fe chlorides dissolve at the hot-leg and subsequently deposit at the cold-leg due to activity changes driven by the thermal gradient. By contrast, Cr was not detected in the cold-leg deposits, which is attributed to its weaker temperature dependence on activity, limiting its redeposition under these conditions.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

In-situ formation of stable interface towards Li-in anode for halide solid-state electrolyte

Halide-based solid-state electrolytes (SSEs) are promising candidates for next-generation all-solid-state lithium batteries (ASSLBs) due to their high ionic conductivity and chemical stability. However, their poor interfacial compatibility with lithium metal anode and Li-In alloy significantly hinder practical application due to the requirement for a protective interlayer. In this study, a novel approach to overcome this limitation is presented by introducing iron (Fe) doping into Li 3 InCl 6 (LIC), which enables direct and stable contact with lithium-indium (Li-In) metal without a protective interlayer. Thermodynamic and computational analyses identified Fe 3+ as a suitable dopant based on its similar reduction potential to In 3+ and structural compatibility within the halide lattice. The synthesized 10 at. % Fe-doped LIC exhibits high phase purity, retained ionic conductivity, and notably improved interfacial stability. Full-cell tests using Fe-LIC achieve over 300 cycles with 80 % capacity retention. At the same time, symmetric Li-In/ Fe-LIC/ Li-In cells sustain over 500 h of operation, representing the first reported long-term cycling of LIC-based ASSLB without a protective interlayer. In conclusion, this work establishes Fe doping as an effective strategy to stabilize halide SSEs of In system against Li-In alloy, thereby simplifying cell architecture and advancing the development of safer, high-performance halide-based solid-state electrolytes.

Halide-based solid-state electrolytes↗
DOE OSTI2021

Distinct Mechanisms and Hydricities of Cp*Ir-Based CO 2 Hydrogenation Catalysts in Basic Water

Transition metal-catalyzed reversible hydrogenation of CO 2 to formate in aqueous solutions under ambient conditions is an attractive environmentally friendly strategy for the storage and transportation of H 2 in a liquid chemical form. Here, mechanistic details for the CO 2 hydrogenation by a series of pentamethylcyclopentadienyl iridium(III) (Cp*Ir) complexes with picolinamidate ligands are investigated through density functional theory and kinetic isotope effects (KIEs) calculations. Additionally, adopting a speciation approach at a pH of 8.2, CO 2 hydrogenation mechanisms of 37 distinct protonation states and conformers of Cp*Ir type complexes are investigated. This imposes a challenging test for the computational modeling in terms of providing a consistent correlation with experimental kinetics and KIE studies of multiple catalysts instead of a single complex. Overall, H 2 heterolysis to generate an iridium hydride intermediate was demonstrated to be the rate determining step, which proceeds via distinct pathways depending on the stereoelectronic properties of the ligand. A peculiar iridium dihydride route was also uncovered, which could be optimized to accelerate the H 2 heterolysis. We have further computed the thermodynamic and kinetic hydricities of 74 complexes with [Cp*Ir(L)(H)] q ( q = -1, 0, 1 depending on the ligand charge by deprotonation or protonation events) general formula. Hydricity values do not show any noticeable correlation with the thermodynamics (Δ G ) of [Cp*Ir(L)(HCO 2 – )] formation but exhibit linear correlation with the kinetics (Δ G ‡ ) of electrophilic CO 2 attack to iridium hydride species, especially when the charges and local structural effects of the metal hydrides are considered. Insights gained in this study on (i) the factors determining the preferred pathway for the rate-limiting H 2 heterolysis step, (ii) the correlation between thermodynamic hydricity and the kinetics of CO 2 insertion to iridium hydride species, and (iii) mechanistic analysis via combination of free energy profiles and KIE studies for a series of iridium catalysts will provide guidelines for the design of next-generation catalysts for the reversible H 2 storage through CO 2 utilization.

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

A Machine Learning Framework for Modeling Ensemble Properties of Atomically Disordered Materials

Atomic disorder can strongly influence material properties such as charge transport, optical response, and catalytic activity. However, efficiently modeling these disorder effects remains challenging for first-principles methods due to the cost of sampling large configurational spaces and computing complex physical quantities. Recent advances of machine learning techniques, particularly graph neural networks (GNNs), has enabled the efficient and accurate predictions of complex material properties, offering promising tools for studying disordered systems. In this work, we present a general machine-learning-assisted computational framework that integrates equivariant GNNs with Monte Carlo simulations to compute the thermodynamic and ensemble-averaged functional properties of disordered materials. Using the surface-termination-disordered MXene monolayer Ti 3 C 2 T 2–x as a representative system, we find that electrical conductivity exhibits an emergent peak near the order–disorder phase transition temperature due to the interplay between electron scattering and doping. In contrast, optical conductivity remains largely insensitive to local atomic disorder and reflects the global surface chemical composition. These results highlight the role of atomic disorder in affecting material properties and demonstrate the potential of our approach for statistically modeling disorder effects in a wide range of materials such as high-entropy alloys and spin liquids.

MXene↗
DOE OSTI2026

Concerted Proton–Electron Transfer Minimizes Substituent Effects on Adsorbed Phthalocyanine Electrocatalysis

Molecularly modified electrodes (MMEs) are potent electrocatalysts, but few principles exist for their rational design. Electrocatalysis by soluble molecules depends strongly on substituents that tune the catalyst redox potential (E 1/2 ), but it is unclear if this parameter similarly impacts MME catalysis. Herein, we employ the hydrogen evolution reaction (HER) as a test case for comparing carbon-adsorbed cobalt phthalocyanine (CoPc/C) and cobalt hexadecafluoro-phthalocyanine (CoFPc/C). By correlating HER activity and voltammetric data to total Co surface concentration across a wide range of catalyst loadings, we find that only 5–25% of adsorbed Co sites contribute to the Co(II/I) redox wave and that this subpopulation poorly correlates with catalytic activity. Instead, in the low-loading limit, catalytic activity correlates linearly with the majority Co(II/I)-silent Co population, revealing per-site turnover frequency (TOF) values for HER. Despite a 230 mV difference in Co(II/I) redox potentials, CoPc/C and CoFPc/C display TOF values differing by less than a factor of 3 when compared over a wide potential range. Mechanistic studies point to an inner-sphere concerted proton–electron transfer step as rate-determining, suggesting that the Co–H bond dissociation free energy (BDFE) rather than the Co(II/I) E 1/2 is thermodynamically relevant. Computational studies indicate that the fluoro-substituents lead to compensatory changes in Co(II/I) E 1/2 and Co(I) basicity, leaving the Co–H BDFE largely unchanged between CoPc and CoFPc and thereby manifesting in similar catalytic rates. Furthermore, these results highlight the limited effect of E 1/2 -tuning on MME catalytic activity and motivate the development of methods to directly alter active site–substrate BDFE.

Catalysts↗
DOE OSTI2020

A new magnesium sheet alloy with high tensile properties and room-temperature formability

Lightweight sheet alloys with superior mechanical performance such as high strength, ductility and formability at room temperature (RT) are desirable for high volume automotive applications. However, ductility or formability of metallic alloys at RT are generally inversely related to strength, thereby making it difficult to optimize all three simultaneously. Here we design a new magnesium sheet alloy-ZAXME11100 (Mg-1.0Zn-1.0Al-0.5Ca-0.4Mn-0.2Ce, wt. pct.) via CALPHAD (CALculation of PHAse Diagram) modeling and experimental validation. This new sheet alloy offers an excellent RT formability with a high Index Erichsen (I.E.) value of 7.8 mm in a solution-treated condition (T4), due to its weak and split basal texture and fine grain structure. The new ZAXME 11100 alloy also shows a rapid age-hardening response during post-forming artificial aging treatment at 210 °C for 1hour (T6), resulting in a significant increase of yield strength from 159 MPa (T4) to 270 MPa (T6). The excellent combination of T4 ductility (31%), T4 formability (7.8 mm) and T6 yield strength (270 MPa) in this new magnesium alloy is comparable to that of common 6xxx series aluminum sheet alloys. Thus, this new magnesium sheet alloy is highly attractive for sheet applications in automotive and other industries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nonequilibrium thermodynamics of erasure with superconducting flux logic

We implement a thermal-fluctuation-driven logical bit reset on a superconducting flux logic cell. We show that the logical state of the system can be continuously monitored with only a small perturbation to the thermally activated dynamics at 500 mK. We use the trajectory information to derive a single-shot estimate of the work performed on the system per logical cycle. We acquire a sample of 10 5 erasure trajectories per protocol and show that the work histograms agree with both microscopic theory and global fluctuation theorems. The results demonstrate how to design and diagnose complex, high-speed, and thermodynamically efficient computing using superconducting technology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Implementation of High Temperature and Pressure fluid property interpolation tables

This document describes the use of module property_interpolate by John Doherty (Doherty,2006). The document “Fast Lookup of CO 2 Properties” describes the numerical details associated with lookup table. Rajeh Pawar modified property_interpolate for easier implementation in the FEHM code. George Zyvoloski modified the package to produce tables for water and air at very high temperatures and pressures. He also modified the auxiliary codes(drivers) to use newer datasets from the National Institute of Standards and Technology. We note here that within the files associated with FEHM there three property_interpolate modules. They are interpolate_2a.f90 (water), interpolate_2b.f90 (air) , and interpolate_2c.f90(CO 2 ). The interpolate_2c module is different than the earlier module produced by Rajeh Pawar in that it includes high temperature and pressure data. While FEHM has a number of fluid physics modules, this document will describe only the software and algorithms associated with water, water vapor, and heat (WH).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
DOE OSTI2022

Computing Free Energies with Fluctuation Relations on Quantum Computers

As a central thermodynamic property, free energy enables the calculation of virtually any equilibrium property of a physical system, allowing for the construction of phase diagrams and predictions about transport, chemical reactions, and biological processes. Thus, methods for efficiently computing free energies, which in general is a difficult problem, are of great interest to broad areas of physics and the natural sciences. The majority of techniques for computing free energies target classical systems, leaving the computation of free energies in quantum systems less explored. Recently developed fluctuation relations enable the computation of free energy differences in quantum systems from an ensemble of dynamic simulations. While performing such simulations is exponentially hard on classical computers, quantum computers can efficiently simulate the dynamics of quantum systems. Here, we present an algorithm utilizing a fluctuation relation known as the Jarzynski equality to approximate free energy differences of quantum systems on a quantum computer. In this work we discuss under which conditions our approximation becomes exact, and under which conditions it serves as a strict upper bound. Furthermore, we successfully demonstrate a proof of concept of our algorithm using the transverse field Ising model on a real quantum processor. As quantum hardware continues to improve, we anticipate that our algorithm will enable computation of free energy differences for a wide range of quantum systems useful across the natural sciences.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

THERMODYNAMIC PROPERTIES OF SRF NIOBIUM

Bulk and thin films of niobium are the materials of choice in fabricating superconducting radio frequency (SRF) cavities for modern particle accelerators and quantum computing applications. The thermodynamic properties of Nb are of particular interest in heat management in cryogenic environments. Here, we report the results of measurements of the thermodynamic properties of niobium used in the fabrication of SRF cavities. The temperature and magnetic field dependence of thermal conductivity, Seebeck coefficient, and specific heat capacity was measured on bulk niobium samples.

Dhakal, Pashupati↗