Many-body thermodynamics on quantum computers via partition function zeros
Quantum computers can study thermodynamics by finding zeros of functions in the complex plane.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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.
Physics-coupled data-driven design of high-temperature alloys
We present a materials design loop, which streamlines physics-coupled machine learning (ML) surrogate models to discover new alloy chemistries with improved properties. The efficacy is demonstrated by discovering a high-temperature alumina-forming austenitic (AFA) stainless steel with enhanced creep, followed by experimental validation. The ML models have been trained using a well-curated, highly consistent experimental dataset augmented with synthetic microstructural features from a computational thermodynamic approach. We have populated a large number of hypothetical AFA alloys to explore the high-dimensional composition space and have predicted their creep properties by providing the same synthetic input features obtained from the trained ML models. Uncertainties from the ML training were taken as thresholds for truncating predicted results to identify alloys with improved or deteriorated creep. Individual elemental compositions have been determined via probability density distribution analysis from the group of alloys at the top and bottom of the predicted creep values for further virtual and experimental validations. In conclusion, we anticipate that this workflow can be applied to screen desired conditions, such as chemistry and processing parameters, in high-dimensional space through physics-guided data analytics.
Alloy design↗
A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States
An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. Here we isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A novel design of transitional layer structure between reduced activation ferritic martensitic steels and tungsten for plasma facing materials
Plasma-facing components (PFCs) are among the most critical gaps for fusion energy to establish technical and economic feasibility. Tungsten as a first wall/blanket material in PFCs requires to be integrating with reduced activation ferritic martensitic (RAFM) steels as a structural component. Currently, major drawbacks are the requirement of brazing, the formation of a brittle interface, and a large difference between the coefficients of thermal expansion of tungsten and steel. Here, a novel transitional multilayer structure was designed and investigated to join tungsten and RAFM steels using three interlayers. The composition of each interlayer was selected based on computational thermodynamics and diffusion kinetics to ensure a body-centered cubic (bcc) single-phase structure and prevent the formation of a brittle intermetallic phase region in the temperature range of 600–1150 °C. Although the transitional layer structure was designed for additive manufacturing, spark plasma sintering (SPS) as proof of concept was used to bond the individual layers. Interfaces were investigated using scanning and transmission electron microscopy methods but no layered intermetallic phase was observed. Nanoindentation maps across the interface suggest major hardness differences at the interface between tungsten and the vanadium interlayer, as well as the interface between RAFM steel and the FeCrAl interlayer.
36 MATERIALS SCIENCE↗
Q-POP-Thermo: A general-purpose thermodynamics solver for ferroelectric materials
We report that Q-POP-Thermo is a program designed to compute thermodynamic monodomain equilibrium states and their properties for ferroelectric single crystals and thin films based on the Landau-Ginzburg-Devonshire (LGD) Theory. Utilizing symbolic manipulation with the SymPy Library, the governing equations along with appropriate boundary conditions are solved for speedy minimization of the free energy of a crystal. Utilizing the popular Differential Evolution algorithm, with appropriate hybridization, multiple phase diagrams, such as the pressure-temperature phase diagram for bulk single crystals and the common strain-temperature phase diagram for monodomain thin-film systems can be readily generated. Furthermore, a variety of material properties of stable ferroelectric phases, including dielectric, piezoelectric, and electrocaloric properties, can simultaneously be calculated. Validation studies are presented for both thin-film and single crystal systems to test the effectiveness and capability of the open-source program.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Generalized entanglement capacity of de Sitter space
Near horizons, quantum fields of low spin exhibit densities of states that behave asymptotically like 1 + 1 dimensional conformal field theories. In effective field theory, imposing some short-distance cutoff, one can compute thermodynamic quantities associated with the horizon, and the leading cutoff sensitivity of the heat capacity is found to equal to the leading cutoff sensitivity of the entropy. One can also compute contributions to the thermodynamic quantities from the gravitational path integral. For the cosmological horizon of the static patch of de Sitter space, a natural conjecture for the relevant heat capacity is shown to equal the Bekenstein-Hawking entropy. These observations allow us to extend the well-known notion of the generalized entropy to a generalized heat capacity for the static patch of de Sitter (dS). The finiteness of the entropy and the nonvanishing of the generalized heat capacity suggest it is useful to think about dS as a state in a finite dimensional quantum gravity model that is not maximally uncertain. Published by the American Physical Society 2024
Astronomy & Astrophysics↗
Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories
Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.
Iyer, Gopal R.↗
Optimal thermodynamic conditions to minimize kinetic by-products in aqueous materials synthesis
Abstract Phase diagrams offer substantial predictive power for materials synthesis by identifying the stability regions of target phases. However, thermodynamic phase diagrams do not offer explicit information regarding the kinetic competitiveness of undesired by-product phases. Here we propose a quantitative and computable thermodynamic metric to identify synthesis conditions under which the propensity to form kinetically competing by-products is minimized. We hypothesize that thermodynamic competition is minimized when the difference in free energy between a target phase and the minimal energy of all other competing phases is maximized. We validate this hypothesis for aqueous materials synthesis through two empirical approaches: first, by analysing 331 aqueous synthesis recipes text-mined from the literature; and second, by systematic experimental synthesis of LiIn(IO 3 ) 4 and LiFePO 4 across a wide range of aqueous electrochemical conditions. Our results show that even for synthesis conditions that are within the stability region of a thermodynamic Pourbaix diagram, phase-pure synthesis occurs only when thermodynamic competition with undesired phases is minimized.
Wang, Zheren (ORCID:0000000227429451)↗
Fundamental Issues Identified for Thermodynamic Description of Molten Salt Systems
Thermodynamic modeling of the molten salt systems is critical in the design of advanced nuclear power plants, materials recycling, and many other important clean energy applications. As one of the most capable computational thermodynamic approaches, the CALPHAD (Calculation of Phase Diagrams) method has been widely used for multicomponent thermodynamic database development and further support more physics-based modeling. Although considerable CALPHAD modeling efforts for molten salt systems have been made, we found that no evaluation has been performed regarding the unary thermodynamic models for the molten salt systems despite the fact that more than one pure substance database is available for the CALPHAD modeling. Therefore, the thermodynamic descriptions of the pure salts should be critically evaluated to ensure a high-fidelity molten salt thermodynamic database. In this work, we comprehensively analyzed some selected molten salts using different available thermodynamic descriptions from two databases as an example. One is the SSUB database released by the SGTE (Scientific Group Thermodata Europe), and the other is the FactPS database from the FactSage. Finally, the significant discrepancies observed in the thermodynamic description between the existing CALPHAD pure substance databases and experiments call for an urgent effort to improve.
42 ENGINEERING↗
Predicting Sintering Window of Binder Jet Additively Manufactured Parts Using a Coupled Data Analytics and CALPHAD Approach
Batch-to-batch variation in powder compositions for binder jet additive manufacturing (BJAM) can significantly deter defining an “ideal” sintering window for a given alloy. One way to overcome the problem is by running sintering experiments at various temperatures for each batch of the powder. However, such an approach increases the time required to achieve large-scale production of parts. The predictive capabilities of computational thermodynamic tools like CALPHAD can be leveraged to overcome the challenge, especially for binder jet additive manufacturing, since the process occurs under near-equilibrium conditions. However, calculating the sintering window using CALPHAD can be computationally expensive, considering many possible feedstock compositions within “specification”. Here, we generate high throughput CALPHAD data for nickel-based superalloys to develop machine learning models to predict the sintering window rapidly. The predictive capability of the models has been validated using published results on BJAM of Inconel 718 and 625. Further, validated models are lightweight and can be deployed in an industrial setting to get sintering window in an accelerated manner.
36 MATERIALS SCIENCE↗
Modernist materials synthesis: Finding thermodynamic shortcuts with hyperdimensional chemistry
Synthesis remains a challenge for advancing materials science. A key focus of this challenge is how to enable selective synthesis, particularly as it pertains to metastable materials. This perspective addresses the question: how can “spectator” elements, such as those found in double ion exchange (metathesis) reactions, enable selective materials synthesis? By observing reaction pathways as they happen (in situ) and calculating their energetics using modern computational thermodynamics, we observe transient, crystalline intermediates that suggest that many reactions attain a local thermodynamic equilibrium dictated by local chemical potentials far before achieving a global equilibrium set by the average composition. Furthermore, using this knowledge, one can thermodynamically “shortcut” unfavorable intermediates by including additional elements beyond those of the desired target, providing access to a greater number of intermediates with advantageous energetics and selective phase nucleation. Ultimately, data-driven modeling that unites first-principles approaches with experimental insights will refine the accuracy of emerging predictive retrosynthetic models for complex materials synthesis.
36 MATERIALS SCIENCE↗
Computationally Guided Design of Multiple Impurities Tolerant Electrode (Final Report)
The current project was based on a combined experimental and computational approach, which can help recommend better cathode materials under multiple impurities conditions. The PI will mainly take in charge of experimental and computational thermodynamics of the selected cathode materials for the SOFC applications under multiple impurities. While the co-PI will run the electrochemical tests of the cathodes recommended and eventually the long-term degradation tests. At the end of the project, a multiple tolerant cathode material based on the combined experimental and computational approach will be recommended and the reliability of the commonly used accelerated testing will be evaluated. It will address multiple impurities poisoning effect of SO2, CO2, Cr and H2O on the LSM, LSCF and LNO cathodes by identifying the formation of the detrimental secondary phases by XRD, SEM and TEM techniques. And further recommended cathode material will be subjected to electrochemical testing and the most promising ones will be applied to long-term tests. The hybrid approach the PI proposed will not only be applied to the design of multiple impurities tolerant cathodes but will also be considered in the future oxygen electrode applications in SOECs or reversible SOCs. This hybrid computational and experimental approach includes four sections: 1) Investigation of single impurity poisoning on LSM, LSCF and LNO cathodes in the presence of SO2, CO2, Cr. In this section, LSM, LSCF cathodes from FuelCellMaterials and LNO cathodes from Sol-Gel synthesis will be heat-treated in the above single impurity. And the formation of the secondary phases as well as the corresponding simulations will be cross compared, which shows good agreement between each other. Meanwhile, the accelerated testing approach will be evaluated in these systems compared with the previous published work to understand the reliability of the approach. 2) Investigation of multiple impurities poisoning on LSM, LSCF and LNO cathodes in the presence of Cr+H2O, SO2+Cr and SO2+Cr+H2O conditions. We have also applied these 3 candidate cathodes under these multiple impurities’ conditions and the long-term degradation mechanism of the multiple impurities will be understood with the help of the combined experimental and computational approach. Meanwhile, the synergistic effect of those impurities will be compared with the individual ones in the same cathode system to further reveal the actual operating conditions. 3) Electrochemical testing and polarization of the recommended cathode. We have demonstrated very low polarization resistance in LSCF (core)-LSM (shell) electrodes using MSD process and analyzed impedance spectra using DRT analysis and confirming that the low polarization resistance in LSM infiltrated MSD cells is due to reduction in polarization resistance associated with O2-adsorption process. 4). Long-term degradation testing. We have developed a versatile MSD-based process to deposit various continuous coatings onto porous scaffolds and established baseline for longer term Cr-impurity testing in future projects.
36 MATERIALS SCIENCE↗
A computational Tool for Thermodynamic Properties of Aqueous Species at HTHP
This computational tool provides the standard Gibbs energy of formation for aqueous species, the dielectric constant and ionization constant of water, over a wide range of temperatures, pressures, and densities both above and below the critical point of water.
Computational Tool↗
TChem v3.0: A Software Toolkit for the Analysis of Complex Kinetic Models
The TChem open-source software is a toolkit for computing thermodynamic properties, source term, and source term’s Jacobian matrix for chemical kinetic models that involve gas and surface reactions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Ab initio calculations of material properties for modeling debris (FY2021 Progress Report)
This report summarizes the work done during the first 6 months of the project “Ab initio calculations of material properties for modeling debris” that is part of the Rapid Response Research (R3) Venture. The main focus for this fiscal year (FY) has been the evaluation of classical interatomic potentials for the description of the condensation of uranium oxides in the regimes relevant to fireball conditions. We also delineate how these functions will be used within the molecular dynamics and Monte Carlo frameworks to compute thermodynamic properties of these materials. These descriptions are to be utilized within continuum modeling formulations to generate a more accurate characterization of the homogeneous nucleation process.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Solid-source metal-organic MBE for elemental Ir and Ru films
Thin films of elemental metals play a very important role in modern electronic nano-devices as conduction pathways, spacer layers, spin-current generators/detectors, and many other important functionalities. In this work, by exploiting the chemistry of solid metal-organic source precursors, we demonstrate the molecular beam epitaxy synthesis of elemental Ir and Ru metal thin films. The synthesis of these metals is enabled by thermodynamic and kinetic selection of the metal phase as the metal-organic precursor decomposes on the substrate surface. Film growth under different conditions was studied using a combination of in situ and ex situ structural and compositional characterization techniques. The critical role of substrate temperature, oxygen reactivity, and precursor flux in tuning film composition and quality is discussed in the context of precursor adsorption, decomposition, and crystal growth. Computed thermodynamics quantifies the driving force for metal or oxide formation as a function of synthesis conditions and changes in chemical potential. These results indicate that bulk thermodynamics are a plausible origin for the formation of Ir metal at low temperatures, while Ru metal formation is likely mediated by kinetics.
Materials Science↗
Thermodynamic Modeling of Solid Flux Interactions with Molten Aluminum
Abstract Oxide and dross formation during aluminum melt processing results in a considerable amount of loss of metallic aluminum. The total melt loss generated during melt processing can be greatly reduced by efficient use of flux, particularly for melting aluminum scrap or secondary alloys. Effective use of cover fluxes can significantly reduce dross generation by creating a barrier between the metallic aluminum and the atmosphere. The amount of metallic aluminum trapped within the dross layer can also be reduced by up to 50% using drossing fluxes. However, reducing the industrial average melt loss has remained difficult. To evaluate flux ingredients used in the casting industry, computational thermodynamic software ThermoCalc was used to calculate the driving force for reactions between twenty-one flux ingredients and eighteen common alloying and impurity elements in foundry alloys. The thermodynamic calculations, combined with other properties, were used to provide a desirable list of cover and drossing ingredients, which are being experimentally validated.
Moodispaw, Michael P. (ORCID:0000000267635882)↗
From electronic structure to model application of key reactions for gasoline/alcohol combustion: Hydrogen-atom abstractions by $CH_3\dot{O}$ radical
H-atom abstraction by methoxy radical ($CH_3\dot{O}$) plays an important role in capturing the kinetics of reactions between gasoline components and alcohols. This study focuses on determining the reaction rates and thermodynamic properties of methoxy radical reactions with five gasoline fuel components: n-heptane, iso-octane, 1-hexene, cyclopentane and toluene. Electronic structure calculations were per-formed for all the stationary points with M06-2X/6 -311 ++ g(d,p) method. G3 composite method with atomization method is used for determining Δ f H 0 of all the closed shell and radical species, using which the necessary thermodynamic data of all the species was determined. Coupled cluster theory QCISD(T)/cc-pVXZ (where X = D and T) and Møller-Plesset perturbation theory MP2/cc-pVXZ (where X = D, T and Q) were used to calculate single point energies. Subsequently, rate constants for all hydrogen atom ab-straction channels have been performed using conventional transition state theory with unsymmetric tunneling corrections. A systematic comparison of rates for abstraction from different sites within the same species and same site from different species is done in order to get insights into this reaction class. Here, the computed thermodynamic properties and rate constants were incorporated into a recent gasoline mechanism to investigate the impact of the calculations performed in this work. A shift in predicted NTC (negative temperature coefficient) behavior and significant reduction in model reactivity is observed upon incorporating the rates calculated herein.
33 ADVANCED PROPULSION SYSTEMS↗