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

PBE-GGA predicts the B8↔B2 phase boundary of FeO at Earth’s core conditions

FeO is a crucial component of the Earth’s core, and its thermodynamic properties are essential to developing more accurate core models. It is also a notorious correlated insulator in the NaCl-type (B1) phase at ambient conditions. It undergoes two polymorphic transitions at 300 K before it becomes metallic in the NiAs-type (B8) structure at ~100 GPa. Although its phase diagram is not fully mapped, it is well established that the B8 phase transforms to the CsCl-type (B2) phase at core pressures and temperatures. Here, we report a successful ab initio calculation of the B8↔B2 phase boundary in FeO at Earth’s core pressures. We show that fully anharmonic free energies computed with the Perdew–Burke–Ernzerhof-generalized gradient approximation coupled with thermal electronic excitations reproduce the experimental phase boundary within uncertainties at P > 255 GPa, including the largely negative Clapeyron slope of –52 MPa/K. This study validates the applicability of a standard density functional theory functional to FeO under Earth’s core conditions and demonstrates the theoretical framework that enables complex predictive studies of this region.

58 GEOSCIENCES↗
DOE OSTI2025

Micropolar Elastoplasticity Using a Fast Fourier Transform‐Based Solver

ABSTRACT This work presents a micromechanical spectral formulation for obtaining the full‐field and homogenized response of elastoplastic micropolar composites. A closed‐form radial‐return mapping is derived from thermodynamics‐based micropolar elastoplastic constitutive equations to determine the increment of plastic strain necessary to return the generalized stress state to the yield surface, and the algorithm implementation is verified using the method of numerically manufactured solutions. Then, size‐dependent material response and micro‐plasticity are shown as features that may be efficiently simulated in this micropolar elastoplastic framework. The computational efficiency of the formulation enables the generation of large datasets in reasonable computing times.

42 ENGINEERING↗
DOE OSTI2022

Solvation Thermodynamics of Solutes in Water and Ionic Liquids Using the Multiscale Solvation-Layer Interface Condition Continuum Model

We report molecular assembly processes are generally driven by thermodynamic properties in solutions. Atomistic modeling can be very helpful in designing and understanding complex systems, except that bulk solvent is very inefficient to treat explicitly as discrete molecules. The SLIC/CDC multiscale model combines continuum solvent electrostatics based on the solvent layer interface condition (SLIC) with statistical thermodynamic models for hydrogen bonding and nonpolar modes: cavity formation, dispersion interactions, combinatorial mixing (CDC). The SLIC/CDC model predicts Gibbs energies of solvation for a database of 500 solutes in water with average accuracy better than 1 kcal/mol both for experimental measurements and for explicit-solvent molecular dynamics simulations. The separate SLIC/CDC energy mode values agree quantitatively with those computed from explicit-solvent molecular dynamics. The SLIC/SASA multiscale model combines the SLIC continuum electrostatic model with the solvent accessible surface area (SASA) nonpolar energy mode. The SLIC/SASA model predicts Gibbs energies of solvation with better than 1.4 kcal/mol average accuracy in aqueous systems and better than 1.6 kcal/mol average accuracy in ionic liquids. Both models predict solvation entropies, and are the first implicit-solvation models capable of predicting solvation heat capacities.

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

Techno-Economic Analysis of CSP Incorporating sCO2 Brayton Power Cycles: Trade-Off Between Cost and Performance

Concentrating solar power (CSP) plants, thanks to the implementation of cost-competitive thermal energy storage, represent a dispatchable zero-emission alternative to traditional fossil fuel power plants. Next generation solar towers are expected to adopt high temperature receivers (>700 degrees C) coupled to sCO2-based power blocks, which optimal design is generally pushed towards the maximum cycle efficiency, often neglecting the economic impact with the justification that the main share of the capital cost is represented by the heliostat field. As result, the scientific literature lacks in comprehensive studies on techno-economic evaluation of CSP+sCO2 power plants addressing the important correlation that exists between system cost and performances. This work provides a preliminary techno-economic analysis of a solar power tower comparing four different cycle configurations for the sCO2 power block. Results have been reported on a Pareto front, highlighting the tradeoff between the plant investment cost and the solar-to-electricity plant efficiency. The trends of the optimization variables and cycle results have been reported to give useful insights about proper assumptions for the sCO2 power block design. The recompressed cycle with intercooling resulted as the most promising configuration and it has been further analyzed through a comparison of different solutions on the Pareto front. The cost breakdown of the sCO2 power block has been reported to highlight which components have the greatest impact on the overall plant cost and how they vary along the optimal solutions front. Eventually, the optimization has been repeated introducing a correlation to compute the turbomachinery isentropic efficiencies, to investigate their effect on the techno- economic analysis.

concentrated solar power↗
DOE OSTI2025

Collisional Excitation of HCN by CO to Refine the Modeling of Cometary Comae

Here, we present the first dataset of collisional (de)-excitation rate coefficients of HCN induced by CO, one of the main perturbing gases in cometary atmospheres. The dataset spans the temperature range of 5–50 K. It includes both state-to-state rate coefficients involving the lowest ten and nine rotational levels of HCN and CO, respectively, and the so-called “thermalized” rate coefficients over the rotational population of CO at each kinetic temperature. The derivation of these coefficients exploited the good performance of the statistical adiabatic channel model (SACM) on top of an accurate interaction potential computed at the CCSD(T)-F12b/CBS level of theory. The reliability of the SACM approach was validated by comparison with full quantum calculations restricted at the lowest total angular momentum of the system. These results provide essential input to accurately model the distribution among the rotational energy levels and the abundance of HCN in cometary atmospheres, accounting for deviations from local thermodynamic equilibrium that typically occurs in such environments.

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

Collisional excitation of C2H− by H2: New interaction potential and scattering calculations

Interstellar anions play an important role in astrochemistry as being tracers of the physical and chemical conditions in cold molecular clouds and circumstellar gas. The local thermodynamic equilibrium is generally not fulfilled in media where anions are detected and radiative and collisional data are required to model the observed lines. The C2H− anion has not yet been detected in the interstellar medium; however, collisional data could be used for non-LTE models that would help in identifying the most intense lines. For this purpose, we have computed the first 4D potential energy surface (PES) of the C2H−–H2 complex using an explicitly correlated coupled-cluster approach. The PES is characterized by a single deep minimum with a well-depth of 924.96 cm−1. From this interaction potential, we derived excitation cross sections and rate coefficients of C2H− induced by collisions with para- and ortho-H2. The results obtained for collisions with para-H2 are compared to previous calculations performed using a 2D-PES obtained from an average over H2 rotations.

Chemistry↗
DOE OSTI2025

Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.

Cage effect↗
DOE OSTI2021

Inelastic rate coefficients for collisions of C 4 H - with H 2

Carbon-chain anions were recently detected in the interstellar medium. These very reactive species are used as tracers of the physical and chemical conditions in a variety of astrophysical environments. However, the local thermodynamic equilibrium conditions are generally not fulfilled in these environments. Therefore, collisional as well as radiative rates are needed to accurately model the observed emission lines. We determine in this work the state-to-state rate coefficients of C 4 H- in collision with both ortho- and para-H 2 . A new ab initio 4D potential energy surface was computed using explicitly correlated coupled-cluster procedures. This surface was then employed to determine rotational excitation and de-excitation cross-sections and rate coefficients for the first 21 rotational levels (up to rotational level j 1 = 20) using the close-coupling method, while the coupled-state approximation was used to extend the calculations up to j 1 = 30. State-to-state rate coefficients were obtained for the temperature range 2–100K. The differences between the ortho- and para-H 2 rate coefficients are found to be small.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Mechanistic Insights into Molecular Copper Hydride Catalysis: the Kinetic Stability of CuH Monomers toward Aggregation is a Critical Parameter for Catalyst Performance

The activity of molecular copper hydride (CuH) complexes towards the selective insertion of unsaturated hydrocarbons under mild conditions has contributed significantly to versatile methodologies for upgrading these feedstocks. However, these catalysts are particularly susceptible to deleterious aggregation, leading to the depletion of active CuH species. Little is known about the mechanisms of CuH aggregation, how it influences overall catalyst performance, and how it can be controlled. We address these challenges with mechanistic studies on a model reaction of unactivated alkene hydroboration catalyzed by (IPr*CPh 3 )CuH (LCuH). Here, we report a comprehensive mechanistic investigation of this system, identifying an aggregation pathway that continuously depletes catalytically active LCuH to form inactive CuH clusters during turnover. Deactivation of LCuH is controlled primarily by the competition between the kinetics of the initial LCuH dimerization step and that of alkene insertion. We therefore propose that a more comprehensive understanding of CuH catalyst performance must account for the kinetics of the initial LCuH dimerization step, revising a previously explored thermodynamic understanding of CuH aggregation, where the concentration of active species is controlled by equilibria established between CuH dimers and monomers. With a series of (NHC)CuH congeners (NHC = N-heterocyclic carbene), we demonstrate that ostensibly minor structural modifications to the ligand peripheries can drastically affect the LCuH dimerization kinetics, while maintaining reactivity towards on–cycle alkene insertion. We employed a computational approach based on molecular dynamics simulations to provide an in-depth understanding of how specific structural ligand modifications can substantially increase the kinetic stability of monomeric CuH catalysts. Our combined experimental and computational studies suggest strategies for rational ligand design that can be broadly applied to molecular catalyst systems that are susceptible to deactivation via aggregation pathways.

Ryan, David E. [Pacific Northwest National Laborat↗
DOE OSTI2024

The Bell-Evans-Polanyi relation for hydrogen evolution reaction from first-principles

The versatile Bell-Evans-Polanyi (BEP) relation stipulates the kinetics of a reaction in terms of thermodynamics. Herein, we establish the BEP relation for the hydrogen evolution reaction (HER) from fundamental electrochemical principles leveraging the Butler-Volmer relation for a one-step, one-electron process and the transition state theory. Based on first-principles investigations of HER mechanisms on fourteen metal electrodes, we firmly justify the BEP relation solely using an easy-to compute hydrogen adsorption free energy and universal electrochemical constants.

42 ENGINEERING↗
DOE OSTI2023

Dimensionally reduced machine learning model for predicting single component octanol–water partition coefficients

Abstract MF-LOGP, a new method for determining a single component octanol–water partition coefficients ( $$LogP$$ LogP ) is presented which uses molecular formula as the only input. Octanol–water partition coefficients are useful in many applications, ranging from environmental fate and drug delivery. Currently, partition coefficients are either experimentally measured or predicted as a function of structural fragments, topological descriptors, or thermodynamic properties known or calculated from precise molecular structures. The MF-LOGP method presented here differs from classical methods as it does not require any structural information and uses molecular formula as the sole model input. MF-LOGP is therefore useful for situations in which the structure is unknown or where the use of a low dimensional, easily automatable, and computationally inexpensive calculations is required. MF-LOGP is a random forest algorithm that is trained and tested on 15,377 data points, using 10 features derived from the molecular formula to make $$LogP$$ LogP predictions. Using an independent validation set of 2713 data points, MF-LOGP was found to have an average $$RMSE$$ RMSE = 0.77 ± 0.007, $$MAE$$ MAE = 0.52 ± 0.003, and $${R}^{2}$$ R 2 = 0.83 ± 0.003. This performance fell within the spectrum of performances reported in the published literature for conventional higher dimensional models ( $$RMSE$$ RMSE = 0.42–1.54, $$MAE$$ MAE = 0.09–1.07, and $${R}^{2}$$ R 2 = 0.32–0.95). Compared with existing models, MF-LOGP requires a maximum of ten features and no structural information, thereby providing a practical and yet predictive tool. The development of MF-LOGP provides the groundwork for development of more physical prediction models leveraging big data analytical methods or complex multicomponent mixtures. Graphical Abstract

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

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗
DOE OSTI2020

MARTINI-Compatible Coarse-Grained Model for the Mesoscale Simulation of Peptoids

Peptoids (poly-N-substituted glycines) are a class of synthetic polymers that are regioisomers of peptides (poly-C-substituted glycines), in which the point of side-chain connectivity is shifted from the backbone C to the N atom. Peptoids have found diverse applications as peptidomimetic drugs, protein mimetic polymers, surfactants, and catalysts. Computational modeling is valuable in the understanding and design of peptoid-based nanomaterials. In this work, we report the bottom-up parameterization of coarse-grained peptoid force fields based on the MARTINI peptide force field against all-atom peptoid simulation data. Our parameterization pipeline iteratively refits coarse-grained bonded interactions using iterative Boltzmann inversion and nonbonded interactions by matching the potential of mean force for chain extension. We assure good sampling of the amide bond cis/trans isomerizations in the all-atom simulation data using parallel bias metadynamics. We develop coarse-grained models for two representative peptoids—polysarcosine (poly(N-methyl glycine)) and poly(N-((4-bromophenyl)ethyl)glycine)—and show their structural and thermodynamic properties to be in excellent accord with all-atom calculations but up to 25-fold more efficient and compatible with MARTINI force fields. Here, this work establishes a new rigorously parameterized coarse-grained peptoid force field for the understanding and design of peptoid nanomaterials at length and time scales inaccessible to all-atom calculations.

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

Quasi‐Binary Transition Metal Dichalcogenide Alloys: Thermodynamic Stability Prediction, Scalable Synthesis, and Application

Abstract Transition metal dichalcogenide (TMDCs) alloys could have a wide range of physical and chemical properties, ranging from charge density waves to superconductivity and electrochemical activities. While many exciting behaviors of unary TMDCs have been demonstrated, the vast compositional space of TMDC alloys has remained largely unexplored due to the lack of understanding regarding their stability when accommodating different cations or chalcogens in a single‐phase. Here, a theory‐guided synthesis approach is reported to achieve unexplored quasi‐binary TMDC alloys through computationally predicted stability maps. Equilibrium temperature–composition phase diagrams using first‐principles calculations are generated to identify the stability of 25 quasi‐binary TMDC alloys, including some involving non‐isovalent cations and are verified experimentally through the synthesis of a subset of 12 predicted alloys using a scalable chemical vapor transport method. It is demonstrated that the synthesized alloys can be exfoliated into 2D structures, and some of them exhibit: i) outstanding thermal stability tested up to 1230 K, ii) exceptionally high electrochemical activity for the CO 2 reduction reaction in a kinetically limited regime with near zero overpotential for CO formation, iii) excellent energy efficiency in a high rate Li–air battery, and iv) high break‐down current density for interconnect applications. This framework can be extended to accelerate the discovery of other TMDC alloys for various applications.

Hemmat, Zahra↗
DOE OSTI2025

Unveiling the thermite-driven lithium fire ignition in solid-state batteries

Here, this study challenges the assumption of the non-flammability of lithium metal all-solid-state batteries (LiSSBs) and other lithium metal batteries without flammable electrolytes. Through thermodynamic calculations and ex situ experiments, we reveal for the first time the risk of thermite reactions between lithium metal and LiFePO 4 in both charged and discharged states. Reactivity is worsened by excess lithium metal in the cell, reaching final maximum adiabatic temperatures of 2,500°C in the charged state, which is hot enough to boil lithium. The thermite reaction triggers spontaneously at 500°C, with poor surface contact, while increasing surface contact through mixing initiates the reaction at room temperature in an inert environment. Despite its fast kinetics, this reaction is transport limited due to lithium passivation, leading to long burn times and reignition risks. Given the risk of lithium metal contacting the cathode during failure, understanding these reactions is crucial for ensuring the safe deployment of LiSSBs.

LiFePO4↗
DOE OSTI2025

Site-specific surface reactivity on SnO 2 : Evaluating selective atomic layer deposition processes

Area selective atomic layer deposition (AS-ALD) is a bottom-up synthesis approach with potential for deposition with molecular level precision. Here, the site-specific hydration of metal oxide substrates, combined with surface H 2 O-selective ALD processes, provides a potentially powerful path to targeted synthesis. Density functional theory (DFT) calculations are used to predict the thermodynamics of ALD precursor reactivity and hydration for (001), (101), (110), and (100) rutile SnO 2 facets as a function of temperature. Trimethylaluminum (TMA) and dimethyl aluminum isopropoxide (DMAI) dimers are predicted to react with both dehydrated and hydrated SnO 2 (001), (101), and (110) facets at ALD-relevant temperatures, while the SnO 2 (100) facet is predicted to be uniquely unreactive with TMA and DMAI monomers as well as dehydrate near 177 °C making this facet more amenable to targeted ALD. In situ ellipsometric studies of Al 2 O 3 ALD on polycrystalline SnO 2 at 150 °C are consistent with the computational predictions of rapid and unselective nucleation, in stark contrast to inhibited and selective ALD on isostructural rutile TiO 2 .

Atomic Layer Deposition↗
DOE OSTI2021

Analytic and numerical vibronic spectra from quasi-classical trajectory ensembles

The truncated Wigner approximation to quantum dynamics in phase space is explored in the context of computing vibronic line shapes for monomer linear optical spectra. We consider multiple model potential forms including a shifted harmonic oscillator with both equal and unequal frequencies on the ground and excited state potentials as well as a shifted Morse potential model. For the equal-frequency shifted harmonic oscillator model, we derive an analytic expression for the exact vibronic line shape that emphasizes the importance of using a quantum mechanical distribution of phase space initial conditions. For the unequal-frequency shifted harmonic oscillator model, we are no longer able to obtain an exact expression for the vibronic line shape in terms of independent deterministic classical trajectories. We show how one can rigorously account for corrections to the truncated Wigner approximation through nonlinear responses of the line shape function to momentum fluctuations along a classical trajectory and demonstrate the qualitative improvement in the resulting spectrum when the leading-order quantum correction is included. Finally, we numerically simulate absorption spectra of a highly anharmonic shifted Morse potential model. We find that, while finite quantization and the dissociation limit are captured with reasonable accuracy, there is a qualitative breakdown of the quasi-classical trajectory ensemble’s ability to describe the vibronic line shape when the relative shift in Morse potentials becomes large. The work presented here provides clarity on the origin of unphysical negative features known to contaminate absorption spectra computed with quasi-classical trajectory ensembles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗