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 577 records · Page 32

The performance of ensemble-based free energy protocols in computing binding affinities to ROS1 kinase

Optimization of binding affinities for compounds to their target protein is a primary objective in drug discovery. Herein we report on a collaborative study that evaluates a set of compounds binding to ROS1 kinase. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to rank the binding free energies. The predicted binding free energies from ESMACS simulations show good correlations with experimental data for subsets of the compounds. Consistent binding free energy differences are generated for TIES and ESMACS. Although an unexplained overestimation exists, we obtain excellent statistical rankings across the set of compounds from the TIES protocol, with a Pearson correlation coefficient of 0.90 between calculated and experimental activities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Entropy Metal Sulfide Nanoparticles Promise High-Performance Oxygen Evolution Reaction

Transition metal sulfide with a multi-elemental nature represents a class of promising catalysts for oxygen evolution reaction (OER) owing to its good catalytic activity. However, its synthesis remains a grand challenge due to the thermodynamic immiscibility of the constituent multimetallic elements in a sulfide structure. Herein, we for the first time report a high-entropy metal sulfide (HEMS, i.e., (CrMnFeCoNi)Sx) solid solution nanoparticles. Computational and X-ray photoelectron spectroscopy analysis suggest that the (CrMnFeCoNi)S x exhibits a synergistic effect among metal atoms that leads to desired electronic states to enhance OER activity. The (CrMnFeCoNi)Sx nanoparticles showed one of the best activities (low overpotential 295 mV at 100 mA cm -2 in 1 M KOH solution) and good durability (only slight polarization after 10 hours by chronopotentiometry) compared with the unary, binary, ternary, and quaternary sulfide counterparts. This work opens up a new synthesis paradigm for high-entropy compound nanoparticles for highly efficient electrocatalysis applications.

25 ENERGY STORAGE↗

Solvent Stabilization of Protic Oxonium/Ammonium Intermediates in Cation Radical Cyclization Reactions Investigated via Computational Approaches

Anodic intramolecular cyclization reactions have substantial synthetic utility for formation of cyclic carbon–carbon or carbon–heteroatom bonds. For cases of intramolecular trapping of a cation radical by a protic nucleophile, the cyclization step coincides with a substantial increase in substrate acidity and thus may exhibit particularly pronounced solvent effects. In this computational work, we employ both quantum chemical (QM) and quantum mechanics/molecular mechanics (QM/MM) methods to compute solvent effects on free energy profiles for cyclization and deprotonation reaction steps for cation radical intermediates of substrates representative for anodic intramolecular cyclizations. We find substantial solvent contribution to the thermodynamic driving force for cation radical cyclization; for example, methanol and tetrahydrofuran solvents provide ∼30–35 kJ/mol driving force to form cyclic oxonium cation radicals and ∼15–25 kJ/mol driving force to form cyclic ammonium cation radicals, compared to baseline reactions in dichloromethane solvent. Given that these solvent shifts are on par with the innate cyclization reaction thermodynamics, the choice of solvent plays a crucial role in promoting/driving the cation radical cyclization step. Methanol is particularly effective at facilitating rapid deprotonation of the cyclic cation radical intermediate, which may lead to the full electrochemical process (e.g., second electron transfer) proceeding heterogeneously at the anode.

Khan, Shahriar [ORNL] (ORCID:0000000289138430)↗

CeTi 2 O 6 —A Promising Oxide for Solar Thermochemical Hydrogen Production

A large entropy of reduction is crucial in achieving materials capable of high-efficiency solar thermochemical hydrogen (STCH) production through two-step thermochemical water splitting cycles. We have recently demonstrated that the onsite electronic entropy of reduction attains an extreme value of 4.26 kB at 1500 K in Ce 4+ → Ce 3+ redox reactions, which explains the high performance and uniqueness of CeO 2 as an archetypal STCH material. However, ceria requires high temperatures (T > 1500 °C) to achieve a reasonable reduction extent because of its large reduction enthalpy, which is a major obstacle in practical applications. Therefore, new materials with a large entropy of reduction and lower reduction enthalpy are required. Here, we perform a systematic screening to search for Ce 4+ -based oxides which possess thermodynamics superior to CeO 2 for STCH production. We first search the Inorganic Crystal Structure Database (ICSD) and literature for Ce 4+ -based oxides and subsequently use density functional theory to compute their reduction enthalpies (i.e., oxygen vacancy formation energies). Here, we find that CeTi 2 O 6 with the brannerite structure is the most promising candidate for STCH because it possesses three essential characteristics of an STCH material: (i) a smaller reduction enthalpy compared to ceria yet large enough to split water, (ii) a high thermal stability, as reported experimentally, and (iii) a large entropy of reduction associated with Ce 4+ → Ce 3+ redox. Our proposed design strategy suggests that further exploration of Ce 4+ oxides for STCH production is warranted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Improved Precision Scaling for Simulating Coupled Quantum-Classical Dynamics

In this paper, we present a superpolynomial improvement in the precision scaling of quantum simulations for coupled quantum-classical systems. Such systems are found in, e.g., molecular-dynamics simulations within the Born-Oppenheimer approximation. By employing a framework based on the Koopman–von Neumann formulation of classical mechanics, we express the Liouville equation of motion as unitary dynamics and utilize phase kickback from a dynamical quantum simulation to calculate the quantum forces acting on classical particles. This approach allows us to simulate the dynamics of these classical particles without the overheads associated with measuring gradients and solving the equations of motion on a classical computer, resulting in a superpolynomial advantage at the price of increased space complexity. We demonstrate that these simulations can be performed in both microcanonical and canonical ensembles, enabling the estimation of thermodynamic properties from the prepared probability density. Published by the American Physical Society 2024

Simon, Sophia (ORCID:0000000332612338)↗

Wet-Bulb Temperature from Pressure, Relative Humidity, and Air Temperature

Abstract The algorithms currently in use for calculating the wet-bulb temperature from pressure, relative humidity, and air temperature have errors of a degree or more. Here, the equation for the thermodynamic wet-bulb temperature is derived using the Rankine–Kirchhoff approximations, which are highly accurate for meteorological and climatological applications. Likewise, equations are derived for the thermodynamic ice-bulb temperature and for the psychrometric wet-bulb and ice-bulb temperatures. The equations are fast to solve, requiring only a microsecond on a current laptop computer, and the results match empirical data to within hundredths of a degree. Furthermore, the equations reveal the existence of air temperatures and humidities that generate bistability: Evaporating liquid water is stable at a liquid temperature above freezing and, simultaneously, sublimating ice is stable at an ice temperature below freezing.

Romps, David M↗

Determining the Role of Synthesis and Reaction Conditions on the Exposed Surface and Shape of ..beta..-Mo2C Catalysts

The production of drop-in quality biofuels typically can be achieved through catalytic fast pyrolysis (CFP) of biomass, followed by hydrodeoxygenation (HDO) of the vapor-phase product to reduce the oxygen content and stabilize the bio-oil product. To favor selectivity toward desired deoxygenation products, bifunctional catalysts with both acidic and hydrogenation sites are often employed. Due to both its oxophilicity and its ability to promote hydrogenation reactions, molybdenum carbide catalysts have shown promise for promoting the HDO of CFP vapors. However, the surface structure and composition of carbide catalyst particles are not well understood, precluding the establishment of fundamental structure-function relationships needed for the rational development of molybdedum carbide for HDO, as well as other reactions. For orthorhombic ..beta..-Mo2C, the (100) facet is largely studied for computational surface-reaction analyses, possibly due to its similarities with the close-packed (111) facet of noble-metal catalysts. However, it has previously been postulated through ab initio thermodynamic studies that the (111), (110), and (011) facets are most stable, though there is no consensus. Furthermore, while a majority of previous surface stability analyses of ..beta..-Mo2C focus on bulk terminations, the highly carburizing synthesis conditions and the use of ..beta..-Mo2C in carbon-containing reaction applications could introduce non-stoichiometric coverages of surface carbon. As such, this study uses ab initio thermodynamics calculations to assess the relative surface stability of the low-index facets of ..beta..-Mo2C, including consideration of varying the coverage of surface carbon. From this information, Wulff constructions are utilized to predict the synthesis and reaction condition-dependent shape and surface carbon composition of ..beta..-Mo2C catalyst particles.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗

Alchemical Free Energy Estimators and Molecular Dynamics Engines: Accuracy, Precision, and Reproducibility

The binding free energy between a ligand and its target protein is an essential quantity to know at all stages of the drug discovery pipeline. Assessing this value computationally can offer insight into where efforts should be focused in the pursuit of effective therapeutics to treat a myriad of diseases. In this work, we examine the computation of alchemical relative binding free energies with an eye for assessing reproducibility across popular molecular dynamics packages and free energy estimators. The focus of this work is on 54 ligand transformations from a diverse set of protein targets: MCL1, PTP1B, TYK2, CDK2, and thrombin. These targets are studied with three popular molecular dynamics packages: OpenMM, NAMD2, and NAMD3 alpha. Trajectories collected with these packages are used to compare relative binding free energies calculated with thermodynamic integration and free energy perturbation methods. The resulting binding free energies show good agreement between molecular dynamics packages with an average mean unsigned error between them of 0.50 kcal/mol. The correlation between packages is very good, with the lowest Spearman’s, Pearson’s and Kendall’s tau correlation coefficients being 0.92, 0.91, and 0.76, respectively. Agreement between thermodynamic integration and free energy perturbation is shown to be very good when using ensemble averaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic and Kinetic Activity Descriptors for the Catalytic Hydrogenation of Ketones

Activity descriptors are a powerful tool for the design of catalysts than can efficiently utilize H 2 with minimal energy losses. In this study, we develop the use of hydricity and H - self-exchange rates as thermodynamic and kinetic descriptors for the hydrogenation of ketones by molecular catalysts. Two complexes with known hydricity, HRh(dmpe) 2 and HCo(dmpe) 2 , were investigated for the catalytic hydrogenation of ketones under mild conditions (1.5 atm, 25 °C). The rhodium catalyst proved to be an efficient catalyst for a wide range of ketones, whereas the cobalt catalyst could only hydrogenate electron-deficient ketones. Using a combination of experiment and electronic structure theory, thermodynamic hydricity values were established for 46 alkoxide/ketone pairs in both MeCN and THF solvent. Through comparison of the hydricities of the catalysts and substrates, it was determined that catalysis was only observed for catalyst/ketone pairs with an exergonic H - transfer step. Mechanistic studies revealed that H - transfer was rate-limiting step for catalysis, allowing for the experimental and computation construction of linear free-energy relationships (LFERs) for H - transfer. Further analysis revealed the LFERs could be reproduced using Marcus theory, in which the H - self-exchange rates for the HRh/Rh + and ketone/alkoxide pairs were used to predict the experimentally measured catalytic barriers within 2 kcal mol -1 . Finally, these studies significantly expand the scope of catalytic reactions that can be analyzed with a thermodynamic hydricity descriptor and firmly establish Marcus theory as a valid approach to develop kinetic descriptors for designing catalysts for H - transfer reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predictive coarse-grained (CG) modeling of morphologies in polymer nanocomposites with specific and directional intermolecular interactions (Final Report)

The overarching goal of the proposed work was to develop predictive models for investigating structure and dynamics in soft materials with chemistries that have specific and directional molecular interactions. The motivation behind studying materials with specific and directional interactions lies in the many desirable features these interactions provide when designing novel soft materials. Soft materials with specific and direction interactions (such as hydrogen bonds or H-bonds) can have a) thermally reversible phase behavior with different functions with varying temperature, b) precisely tuned nanostructure with desirable geometries that afford unique physical properties (e.g., color response, mechanical properties) and c) well-mixed/blended morphologies that are useful for variety of applications in energy field (e.g., materials for batteries require use of blended polymers where one domain gives superior mechanical properties and one domain promotes electrical conduction). Notably, biology makes extensive use of specific and directional interactions, in many cases based on H-bonds, to construct materials with precisely defined architectures and properties. Engineered soft materials with precisely tuned nanostructures and improved processiblity through thermoresponsive phase behavior are useful in numerous applications that are relevant to the Department of Energy (DOE) including high efficiency electronic devices, light-weight high-strength composite materials for batteries and fuel cells, and polymer membranes for separations, etc. While past computational studies have been tremendously useful in understanding molecular phenomena and guiding synthesis of new macromolecular soft materials for a wide variety of applications, the inability to capture small scale specific and directional interactions alongside macromolecular length and time scales represented a key limitation of most studies to date. Our work in this project addressed this grand challenge in computational materials chemistry, i.e., ability to model the anisotropic, directional, and specific interactions that govern the behavior of many macromolecular soft matter systems of interest, thus, has greatly expanded the predictive potential of simulations. Specifically the key outcomes were: successful development of new coarse-grained (CG) polymer models to study generic and specific polymer chemistries in which hydrogen-bonding interactions are dominant. These CG models were then used in molecular simulations to study structure and thermodynamics in polymer nanocomposites and blends; some studies were conducted in collaboration with experimentalists. We also published a perspective and a viewpoint which included some of the work we completed in this DOE project; we believe these perspective and viewpoint articles guide other researchers in the soft materials community on how to extend the computational approaches and models we have developed for the purposes of their studies.

36 MATERIALS SCIENCE↗

Modeling and Optimization of a Rotating Packed Bed Contactor with a Tetraamine-Appended Metal−Organic Framework for CO 2 Capture

A potential contactor technology for sorbent-based CO 2 capture is the rotating packed bed that contains separate sections for continuous adsorption and desorption. A heat exchanger can be embedded to remove heat in the adsorption section and add heat in the desorption section. In this work, we develop a two-dimensional (2D) model of a rotating packed bed for use in CO 2 capture applications. Mass and energy balances for the model are developed based on a Ljungström-type air preheater, which accounts for the counter-current axial flow of gas phases in separate sections of the bed and the rotation of a solid sorbent, which cycles between adsorption and desorption sections. The sorbent used for this analysis is the tetraamine-appended metal−organic framework Mg 2 (dobpdc)(3−4− 3), chosen for its stability and affinity for CO 2 at low partial pressures, such as those from a natural gas power plant source. An optimization problem is solved that considers the trade-off between maximizing the productivity of the bed and minimizing energy consumption. Maximum productivity and minimum energy are found to be 8.53 kg/h/m 3 and 3.84 MJ/kg, respectively, when these objectives are optimized independently. It is observed that the flue gas pressure and bed rotational speed are the desired operating variables to vary for model-based design of experiments to reduce uncertainty in parameter estimation, as these two variables yielded the most information content based on the Fisher information matrix.

20 FOSSIL-FUELED POWER PLANTS↗

Upstreamness and downstreamness in input–output analysis from local and aggregate information

Abstract Ranking sectors and countries within global value chains is of paramount importance to estimate risks and forecast growth in large economies. However, this task is often non-trivial due to the lack of complete and accurate information on the flows of money and goods between sectors and countries, which are encoded in input–output (I–O) tables. In this work, we show that an accurate estimation of the role played by sectors and countries in supply chain networks can be achieved without full knowledge of the I–O tables, but only relying on local and aggregate information, e.g., the total intermediate demand per sector. Our method, based on a rank-1 approximation to the I–O table, shows consistently good performance in reconstructing rankings (i.e., upstreamness and downstreamness measures for countries and sectors) when tested on empirical data from the world input–output database. Moreover, we connect the accuracy of our approximate framework with the spectral properties of the I–O tables, which ordinarily exhibit relatively large spectral gaps. Our approach provides a fast and analytically tractable framework to rank constituents of a complex economy without the need of matrix inversions and the knowledge of finer intersectorial details.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lithium oxidation and electrolyte decomposition at Li-metal/liquid electrolyte interfaces

In this work, we examine the evolution of events occurring when a Li metal surface is in contact with a 2 M solution of a Li salt in a solvent or mixture of solvents, via classical molecular dynamics simulations with a reactive force field allowing bond breaking and bond forming. The main events include Li oxidation and electrolyte reduction along with expansion of the Li surface layers forming a porous phase that is the basis for the formation of the solid-electrolyte interphase (SEI) components. Nucleation of the main SEI components (LiF, Li oxides, and some organics) is characterized. The analysis clearly reveals the details of these physical–chemical events as a function of time, during 20 nanoseconds. The effects of the chemistry of the electrolyte on Li oxidation and dissolution in the liquid electrolyte, and SEI nucleation and structure are identified by testing two salts: LiPF 6 and LiCF 3 SO 3 , and various solvents including ethers and carbonates and mixtures of them. The kinetics and thermodynamics of Li 6 F, the core nuclei in the LiF crystal, are studied by analysis of the MD trajectories, and via density functional theory calculations respectively. The SEI formed in this computational experiment is the “native” film that would form upon contact of the Li foil with the liquid electrolyte. As such, this work is the first in a series of computational experiments that will help elucidate the intricate interphase layer formed during battery cycling using metal anodes.

25 ENERGY STORAGE↗

The effects of multispecies Hohlraum walls on stimulated Brillouin scattering, Hohlraum dynamics, and beam propagation

We report experiments and simulations have been conducted to investigate the efficacy of Ta 2 O 5 -lined Hohlraum walls at reducing stimulated Brillouin backscattering (SBS) as well as any subsequent effects on the Hohlraum dynamics and capsule implosions in indirect drive experiments at the National Ignition Facility. Using a 1.1 MJ 400 TW, 351 nm, shaped laser pulse, we measure a 5× reduction in SBS power in the peak of the pulse from the wall on the outer 50° cone beams. The SBS spectrum indicates a reduction in the high-Z spectral signature when using multispecies wall materials. Detailed hydrodynamic simulations were performed using different heat conduction models with flux limiters. Additional simulations were performed on the plasma maps using the 3D parallel paraxial code pF3D to compare backscatter powers between the pure Au and Ta 2 O 5 -lined Hohlraums. Further analysis, using hydrodynamically equivalent plasmas, shows that the SBS reduction is clearly a result of the added ion Landau damping caused by the oxygen ions and not from differences in plasma conditions. The experimental and simulation results also show an increase in the wall plasma expansion when using the Ta 2 O 5 liner leading to a 70% more oblate implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Posterior Covariance Matrix Approximations

Here, the Davis equation of state (EOS) is commonly used to model thermodynamic relationships for high explosive (HE) reactants. Typically, the parameters in the EOS are calibrated, with uncertainty, using a Bayesian framework and Markov Chain Monte Carlo (MCMC) methods. However, MCMC methods are computationally expensive, especially for complex models with many parameters. This paper provides a comparison between MCMC and less computationally expensive Variational methods (Variational Bayesian and Hessian Variational Bayesian) for computing the posterior distribution and approximating the posterior covariance matrix based on heterogeneous experimental data. All three methods recover similar posterior distributions and posterior covariance matrices. This study demonstrates that for this EOS parameter calibration application, the assumptions made in the two Variational methods significantly reduce the computational cost but do not substantially change the results compared to MCMC.

97 MATHEMATICS AND COMPUTING↗

Screening and Development of Sacrificial Cathode Additives for Lithium‐Ion Batteries

Abstract This work presents a computational screening approach to identify Li‐rich transition‐metal oxide sacrificial cathode additives and provides experimental validation of antifluorite‐structured Li 6 MnO 4 as a potential candidate. Initial attempts to synthesize this compound result in low purity (≤40% by weight) owing to close thermodynamic competition with Li 2 O and MnO at low temperature. However, it is shown that a much higher purity of 85% by weight can be achieved by combining Li excess with rapid cooling from high temperature, which effectively stabilizes the Li 6 MnO 4 phase. The synthesized product delivers a high irreversible Li release capacity that exceeds 700 mAh g −1 by utilizing combined Mn oxidation (Mn 2+/3+ and Mn 3+/4+ ) and O oxidation. These results demonstrate that Li 6 MnO 4 may therefore be useful as a potential sacrificial cathode additive in Li‐ion batteries and motivate further investigation of other structurally‐related compounds. While attempts were made to synthesize two additional compounds among computationally screened candidates, it was not successful to experimentally realize the two candidates. The difficulty of experimental realization of the newly predicted materials remains a challenge and it is suggested that more efforts need to be devoted to developing computational techniques to precisely predict synthesizability and propose potential synthetic routes of the predicted materials.

Kim, Haegyeom↗

Neuromorphic heat transport effects in a molecular junction

Understanding energy transport at the nanoscale is an open and fundamental challenge in the molecular sciences with direct implications for the design of new electronics, computing devices, and materials. While nanoscale energy transport under steady-state conditions has been studied extensively, there is much less known about energy transport under time-dependent driving forces, particularly in the far-from-equilibrium regime. In this work, we use nonequilibrium molecular dynamics simulations and stochastic thermodynamics to investigate energy transport in a well-studied nanoscale system—a molecular junction—subjected to a time-periodic temperature gradient. The primary observation is that molecular junctions can exhibit heat transport hysteresis, a phenomenon in which the heat flux through a system depends not only on the instantaneous value of a time-dependent temperature bias but also on the temporal history of that bias. The presented findings illustrate that molecular junctions can exhibit the specific memory effect—heat transport hysteresis—that is essential for the design of thermal neuromorphic computers. This work elucidates a potential pathway toward the realization of such devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Disorderly Conduct of Benzamide IV: Crystallographic and Computational Analysis of High Entropy Polymorphs of Small Molecules

Benzamide, a simple derivative of benzoic acid and a common intermediate of pharmaceutical compounds, was reported to form two polymorphs in 1832, but the single crystal structure of the more stable form was not solved until 1959. Nearly 50 years later, the second form was characterized by powder diffraction, followed shortly thereafter by characterization of a third form, a polytype of the most thermodynamically stable Form I. These two new forms, Forms II and III, are metastable. Herein, we describe a fourth polymorph, Form IV, discovered by melt crystallization concurrently with its crystallization under confinement at small length scales (<10 nm), where it is stable indefinitely. Form III exists under confinement in larger pores, and melting point data for different pore sizes corroborate the existence of Form IV below 10 nm. Form IV is highly disordered, precluding indexing of powder diffraction data other than hk0 reflections. Nonetheless, a combination of powder X-ray diffraction and computational crystal structure prediction reveals that Form IV contains a 2D motif resembling that of Form II, but with longrange order in the third dimension masked by ubiquitous stacking faults. This approach relies on distilling a large number of candidate structures to a few possible disorder models based on benzamide tetrads that organize in 2D parquet-like tiles, with organization along the third dimension, that can be modeled with various stacking fault configurations having distinct intermolecular interactions and translations in the dimension orthogonal to the tiling planes. These observations reveal a bewildering crystallographic complexity for such a simple molecule. Nonetheless, the approach described herein demonstrates that challenging structures that may be abandoned prematurely because of poor crystallinity, twinning, or disorder can be solved.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗