Search NASA⌕ Search

SEARCH · Search NASA

Results for “Free energy”

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 OSTI2023

Minimum Free-Energy Shapes of Ag Nanocrystals: Vacuum vs Solution

Here, we use two variants of replica-exchange molecular dynamics (MD) simulations, parallel tempering MD and partial replica exchange MD, to probe the minimum free-energy shapes of Ag nanocrystals containing 100–200 atoms in a vacuum, ethylene glycol (EG) solvent, and EG solvent with a PVP polymer containing 100 repeat units. Our simulations reveal a shape intermediate between a Dh and an Ih, a Dh-Ih, that has distinct structural signatures and magic sizes. We find several prominent features associated with entropy: pure FCC nanocrystals are less common than FCC crystals containing stacking faults, and crystals with the minimum potential energy are not always preferred over the range of relevant temperatures. The shapes of the nanocrystals in solution are influenced by the chemical identities of the solution-phase molecules. Comparing Ag nanocrystal shapes in EG to those in an EG+PVP solution, we find more icosahedra in EG and more decahedra in EG+PVP across all of the nanocrystal sizes probed in this study. At certain critical sizes, nanocrystal shapes can change dramatically with the addition and removal of a single atom or with a change in temperature at a fixed size. The information in our study could be useful in efforts to devise processing routes to achieve selective nanocrystal shapes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Neural Network Approach to Predict Gibbs Free Energy of Ternary Solid Solutions

Here, we present a data-centric deep learning (DL) approach using neural networks (NNs) to predict the thermodynamics of ternary solid solutions. We explore how NNs can be trained with a dataset of Gibbs free energies computed from a CALPHAD database to predict ternary systems as a function of composition and temperature. We have chosen the energetics of the FCC solid solution phase in 226 binaries consisting of 23 elements at 11 different temperatures to demonstrate the feasibility. The number of binary data points included in the present study is 102,000. We select six ternaries to augment the binary dataset to investigate their influence on the NN prediction accuracy. We examine the sensitivity of data sampling on the prediction accuracy of NNs over selected ternary systems. It is anticipated that the current DL workflow can be further elevated by integrating advanced descriptors beyond the elemental composition and more curated training datasets to improve prediction accuracy and applicability.

42 ENGINEERING↗
DOE OSTI2020

Subleading corrections to the free energy in a theory with N 5/3 scaling

We numerically investigate the sphere partition function of a Chern-Simons-matter theory with SU(N) gauge group at level k coupled to three adjoint chiral multiplets that is dual to massive IIA theory. Beyond the leading order N 5/3 behavior of the free energy, we find numerical evidence for a term of the form (2/9) log N. We conjecture that this term may be universal in theories with N 5/3 scaling in the large-N limit with the Chern-Simons level k held fixed.

1/N Expansion↗
DOE OSTI2021

Potent Noncovalent Inhibitors of the Main Protease of SARS-CoV-2 from Molecular Sculpting of the Drug Perampanel Guided by Free Energy Perturbation Calculations

Starting from our previous finding of 14 known drugs as inhibitors of the main protease (M pro ) of SARS-CoV-2, the virus responsible for COVID-19, we have redesigned the weak hit perampanel to yield multiple noncovalent, nonpeptidic inhibitors with ca. 20 nM IC 50 values in a kinetic assay. Free-energy perturbation (FEP) calculations for M pro -ligand complexes provided valuable guidance on beneficial modifications that rapidly delivered the potent analogues. The design efforts were confirmed and augmented by determination of high-resolution X-ray crystal structures for five analogues bound to M pro . Results of cell-based antiviral assays further demonstrated the potential of the compounds for treatment of COVID-19. In addition to the possible therapeutic significance, the work clearly demonstrates the power of computational chemistry for drug discovery, especially FEP-guided lead optimization.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An ab initio free energy study of the reaction mechanism and rate-limiting steps of CO 2 capture by aqueous glycine

Aqueous amino acids are promising absorbents for direct air capture (DAC) of CO 2 . Herein, we investigate the possibility of kinetic control of CO 2 absorption with aqueous anionic glycine (GLY - ) by employing extensive ab initio molecular dynamics simulations, free energy analysis, and reaction rate theory. We find that first GLY - binds to CO 2 by overcoming a barrier (7.4 kcal/mol) to form a zwitterion intermediate, which then releases a proton by overcoming a similar barrier. Despite the similarity in the barrier, zwitterion formation appears to be the rate-limiting step because it is two orders of magnitude slower (microseconds) than the proton release step. This is predominantly due to stronger nonequilibrium solvent effects for the former that cause many barrier-recrossing events and effectively slow down the reaction rate. Such a detailed fundamental understanding of the amino acid-based CO 2 -absorption mechanism and rates is key to improving the kinetic efficiency of DAC technology.

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

Surface Gibbs free energy analyses of Sr segregation in lanthanum strontium iron oxide

Sr surface segregation in mixed-cation oxides often causes degradation of chemical or electrochemical performance. In this work, the thermodynamic driving force for Sr segregation is investigated to explain the Sr segregation behavior in La0.5Sr0.5FeO3 (LSF). Specifically, Density Functional Theory is used to calculate the surface Gibbs free energy of various surface terminations of cubic (100) LSF as a function of chemical potential and temperature. Here, SrO-terminated surfaces are predicted to be dominant, promoting Sr segregation at high temperatures and low oxygen partial pressure.

08 HYDROGEN↗
DOE OSTI2021

Proton-Coupled Defects Impact O–H Bond Dissociation Free Energies on Metal Oxide Surfaces

Proton-coupled electron transfer (PCET) reactions on metal oxides require coupling between proton transfer at the solid-liquid interface and electron transfer involving defects at or near the band edge. Herein, hybrid functional periodic density functional theory is used to elucidate the impact of proton-coupled defects on the bond dissociation free energies (BDFEs) of O-H bonds on anatase TiO 2 surfaces. These O-H BDFEs are directly related to interfacial PCET thermochemistry. Comparison between geometrically similar O-H bonds associated with different defect types, namely conduction d-band electrons or valence p-band holes, reveals that the BDFEs differ by ~81 kcal/mol (3.50 eV), comparable to the wide TiO 2 band gap. These differences are shown to be determined primarily by differences in electron transfer driving forces, which are analyzed by using band energies and inner-sphere reorganization energies within a Marcus theory framework. Furthermore, these fundamental insights about the impact of proton-coupled defects on PCET thermochemistry at semiconductor surfaces have broad implications for electrocatalysis.

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

Correlating Solvation Free Energy to Electrolyte Properties for Lithium Metal Batteries

The electrolyte plays a critical role in lithium metal batteries. In particular, ion solvation profoundly impacts key electrolyte properties and battery performance. Here, in this study, we systematically investigate solvation-property relationships in a series of electrolytes with different solvent-diluent ratios. Through potentiometric techniques that measure the relative solvation free energies of electrolytes, we find that weaker solvation correlates with larger ion clusters, lower ionic conductivity and diffusion coefficient, and superior electrochemical stability. Weaker solvation leads to the formation of a small number of Li clusters with large hydrodynamic radii, which lowers the Li + diffusivity and ionic conductivity of the electrolyte. Concurrently, weaker solvation leads to improved electrochemical stability at both the cathode and anode interfaces. Understanding these solvation-property relationships and trade-offs is important to designing electrolytes for optimized lithium metal battery performance.

25 ENERGY STORAGE↗
DOE OSTI2023

Systematic modification of functionality in disordered elastic networks through free energy surface tailoring

A combined machine learning–physics–based approach is explored for molecular and materials engineering. Specifically, collective variables, akin to those used in enhanced sampled simulations, are constructed using a machine learning model trained on data gathered from a single system. Through the constructed collective variables, it becomes possible to identify critical molecular interactions in the considered system, the modulation of which enables a systematic tailoring of the system’s free energy landscape. To explore the efficacy of the proposed approach, we use it to engineer allosteric regulation and uniaxial strain fluctuations in a complex disordered elastic network. Its successful application in these two cases provides insights regarding how functionality is governed in systems characterized by extensive connectivity and points to its potential for design of complex molecular systems.

42 ENGINEERING↗
DOE OSTI2026

Ensemble Effects on Hydroxide Bond Dissociation Free Energies in Polyoxovanadate Clusters

Understanding structure-property relationships is foundational to numerous modern chemistries, such as proton-coupled electron transfer (PCET). However, an experimentally measured property is the result of the behavior from an ensemble of molecules. Neglecting ensemble effects, especially under complex chemical environments, may obfuscate these relationships and lead to discrepancies between theory and experiment. In this work, we demonstrate the impact of configurational entropy and local chemical environments on hydroxide bond dissociation free energies [BDFE- (O−H)] for a set of polyoxovanadate nanoclusters, at ambient conditions. The O−H bond strengths are investigated via density functional theory (DFT) coupled with statistical thermodynamic analysis and bilinear modeling, and compared with previous experimental results on the same systems, namely electrochemical solutions of: [V 6 O 13−x (OH) x (TRIOL R ) 2 ] −2 (x = 2, 4, 6; R = NO 2 , Me) and [V 6 O 11−x (OMe) 2 (OH) x (TRIOL NO 2 ) 2 ] −2 (x = 2, 4). Interestingly, we find that ensemble effects, even at room temperature, can account for a significant portion of the BDFE(O−H) trend with the degree of reduction via H atom binding, which cannot be fully captured by single-structure, static DFT calculations. Moreover, we find that the ensemble effects may be replicated statistically, requiring only enumeration of energetically accessible H-binding sites. With the ensemble effects resolved, we present a simple bilinear model to reconcile remaining biases between experiment and ensemble-informed theory, which corelate with clusterspecific electronic environment differences. The bilinear model achieves outstanding accuracy vs experiments with a root-mean squared error of 0.4 kcal/mol. Finally, based on the physicochemical characteristics of hydrogen interaction with polyoxometalates, we present a simple methodology that captures the BDFE(O−H) trend while dramatically reducing required DFT calculations by 98% and achieving accuracy within 1 kcal/mol. Overall, this work elucidates the roles and structural origins of configurational entropy and chemical effects on polyoxometalate hydroxide bond energies, with potential applicability to various atomically precise metal oxide systems. Importantly, it introduces models for rapid and highly accurate property calculations in connection with experiments.

Cluster chemistry↗
DOE OSTI2025

Morphology Effects on Free Energies of Proton-Coupled Electron Transfer in Polyoxotungstates

Polyoxotungstates have previously been established to facilitate the hydrogenation of small molecule substrates via hydrogen atom transfer from reactive hydroxyl groups formed at the assembly surface. Understanding structure−function relationships that dictate the thermochemistry and kinetics of protoncoupled electron transfer is key to controlling this chemistry. In this work, we combine comprehensive electrochemical experiments and density functional theory calculations to address how different polyoxotungstate morphologies, specifically W 6 O 19 −2 , W 10 O 32 −4 , SiW 12 O 40 −4 , and P 2 W 18 O 62 −6 , affect the bond dissociation free energies of surface hydroxides (BDFE(O−H)) formed upon reduction of the assembly in acidic media. Our results reveal increasing hydroxide bond strengths with increasing cluster size, and that anisotropic cluster geometries result in substantial thermodynamic differentiation of H-binding sites. We demonstrate an excellent agreement between theory and experiments on the reported BDFE(O−H) values and, importantly, we elucidate how cluster size and shape affect electronic properties (local charges and frontier molecular orbitals), giving rise to sites with increased preference for hydrogen binding, demonstrated in higher BDFE(O−H). Overall, this work aids the understanding and design of polyoxometalates exhibiting surface sites with tailored interaction strengths.

anions↗
DOE OSTI2025

HEPOM: Using Graph Neural Networks for the Accelerated Predictions of Hydrolysis Free Energies in Different pH Conditions

Hydrolysis is a fundamental family of chemical reactions where water facilitates the cleavage of bonds. The process is ubiquitous in biological and chemical systems, owing to water’s remarkable versatility as a solvent. However, accurately predicting the feasibility of hydrolysis through computational techniques is a difficult task, as subtle changes in reactant structure like heteroatom substitutions or neighboring functional groups can influence the reaction outcome. Furthermore, hydrolysis is sensitive to the pH of the aqueous medium, and the same reaction can have different reaction properties at different pH conditions. In this work, we have combined reaction templates and high-throughput ab initio calculations to construct a diverse data set of hydrolysis free energies. The developed framework automatically identifies reaction centers, generates hydrolysis products, and utilizes a trained graph neural network (GNN) model to predict ΔG values for all potential hydrolysis reactions in a given molecule. The long-term goal of the work is to develop a data-driven, computational tool for high-throughput screening of pH-specific hydrolytic stability and the rapid prediction of reaction products, which can then be applied in a wide array of applications including chemical recycling of polymers and ion-conducting membranes for clean energy generation and storage.

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

Prediction of α $IIb$ $β$ 3 integrin structures along its minimum free energy activation pathway

The adhesion protein integrin is a transmembrane heterodimer that plays a pivotal role in cellular processes such as cell signaling and cell migration. To execute its function, integrin undergoes extensive conformational changes from a bent-closed to an extended-open state. Resolving the structures across these changes remains a challenge with both experimental and computational methods, but it is crucial for understanding the activation mechanism of integrin. We address this challenge for the platelet integrin α IIb β 3 by employing finite temperature string method with structures of the images along the initial guess path generated by a multiscale data-driven framework. The full-length all-atom structures along the resulting minimum free energy path between the inactive bent-closed and active extended-open states of α IIb β 3 integrin are consistent with a variety of experimentally resolved structures. Changes in these predicted structures along the path show that the extension and separation of the α and β subunits from the bent-closed to the extended-open state require correlated movements between the subdomain pairs in α IIb β 3 . Furthermore, these results provide new insights into integrin activation mechanism, and the predicted structures have potential applications in guiding the design of integrin-targeting therapeutics.

Dasetty, Siva [University of Chicago, IL (United S↗
DOE OSTI2022

First-principles free energies by hybrid thermodynamic integration for phase equilibria and fission product solubility in molten salts

The overarching goal of this research project is quantitative microscopic modeling of phase diagrams and solubility of actinides and lanthanide-fission products in molten salts for molten salt reactors. Application of molten salts as coolants and liquid fuels for next-generation reactors requires quantitative understanding of phase diagrams, equations of state and transport properties of binary and ternary molten salt mixtures, as well as solubility of fuel, fission products and reactor-wall corrosion products, but these remain beyond the reach of current first-principles computational techniques. We will address this fundamental capability gap using a novel hybrid framework for thermodynamic integration to directly predict absolute free energies of solid and molten salts at any specified conditions of composition, temperature and pressure. This framework will combine a hierarchy of first-principles, molecular dynamics and continuum techniques to deliver first-principles-level accuracy at a fraction of the computational cost of conventional ab initio molecular dynamics approaches. In particular, this technique is expected to provide first-principles accuracy, while off-loading most of the computational work on much less expensive classical molecular dynamics and continuum techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Free energy difference to create the M-OH* intermediate of the oxygen evolution reaction by time-resolved optical spectroscopy

Theoretical descriptors differentiate the catalytic activity of materials for the oxygen evolution reaction by the strength of oxygen binding in the reactive intermediate created upon electron transfer. Recently, time-resolved spectroscopy of a photo-electrochemically driven oxygen evolution reaction followed the vibrational and optical spectra of this intermediate, denoted M-OH*. However, these inherently kinetic experiments have not been connected to the relevant thermodynamic quantities. Here we discover that picosecond optical spectra of the Ti-OH* population on lightly doped SrTiO 3 are ordered by the surface hydroxylation. A Langmuir isotherm as a function of pH extracts an effective equilibrium constant relatable to the free energy difference of the first oxygen evolution reaction step. Thus, time-resolved spectroscopy of the catalytic surface reveals both kinetic and energetic information of elementary reaction steps, which provides a critical new connection between theory and experiment by which to tailor the pathway of water oxidation and other surface reactions.

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

Conformational free energy landscape of a glutamate transporter and microscopic details of its transport mechanism

Removing glutamate from the synaptic cleft is vital for proper function of the brain. Excitatory amino acid transporters mediate this process by uptaking the neurotransmitter from the synaptic cleft back to the cell after its release. The archaeal homolog, Glt Ph , an aspartate transporter fromPyrococcus horikoshii, presents the best structurally characterized model for this family of transporters. In order to transport, Glt Ph undergoes elevator-like conformational changes between inward-facing (IF) and outward-facing (OF) states. Here, we characterize, at an atomic level, the OF⇌IF transition of Glt Ph in differentapo/bound states using a combination of ensemble-based enhanced sampling techniques, employing more than two thousand of coupled simulation replicas of membrane-embedded Glt Ph . The resulting free-energy profiles portray the transition ofapo/bound states as a complex four-stage process, while sodium binding alone locks the structure in one of its states. Along the transition, the transport domain (TD) disengages from the scaffold domain (SD), allowing it to move as a piston sliding vertically with respect to the membrane during the elevator-like motion of TD. Lipid interactions with residues comprising the SD–TD interface directly influence the large-scale conformational changes and, consequently, the energetics of transport. Structural intermediates formed during the transition leak water molecules and may correlate to the uncoupled Cl − ion conductance observed experimentally in both prokaryotic and mammalian glutamate transporters. Mechanistic insights obtained from our study provide a structural framework for better development of therapeutic for neurological disorders.

Science & Technology - Other Topics↗
DOE OSTI2021

Strong coupling expansion of free energy and BPS Wilson loop in $$ \mathcal{N} $$ = 2 superconformal models with fundamental hypermultiplets

As a continuation of the study (in arXiv:2102.07696 and arXiv:2104.12625 ) of strong-coupling expansion of non-planar corrections in \( \mathcal{N} \) = 2 4d superconformal models we consider two special theories with gauge groups SU( N ) and Sp(2 N ). They contain N -independent numbers of hypermultiplets in rank 2 antisymmetric and fundamental representations and are planar-equivalent to the corresponding \( \mathcal{N} \) = 4 SYM theories. These \( \mathcal{N} \) = 2 theories can be realised on a system of N D3-branes with a finite number of D7-branes and O7-plane; the dual string theories should be particular orientifolds of AdS 5 × S 5 superstring. Starting with the localization matrix model representation for the \( \mathcal{N} \) = 2 partition function on S 4 we find exact differential relations between the 1 /N terms in the corresponding free energy F and the \( \frac{1}{2} \) -BPS Wilson loop expectation value \( \left\langle \mathcal{W}\right\rangle \) and also compute their large ’t Hooft coupling ( λ » 1) expansions. The structure of these expansions is different from the previously studied models without fundamental hypermultiplets. In the more tractable Sp(2 N ) case we find an exact resummed expression for the leading strong coupling terms at each order in the 1 /N expansion. We also determine the exponentially suppressed at large λ contributions to the non-planar corrections to F and \( \left\langle \mathcal{W}\right\rangle \) and comment on their resurgence properties. We discuss dual string theory interpretation of these strong coupling expansions.

1/N expansion↗