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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Large-scale application of free energy perturbation calculations for antibody design

Abstract Alchemical free energy perturbation (FEP) is a rigorous and powerful technique to calculate the free energy difference between distinct chemical systems. Here we report our implementation of automated large-scale FEP calculations, using the Amber software package, to facilitate antibody design and evaluation. In combination with Hamiltonian replica exchange, our FEP simulations aim to predict the effect of mutations on both the binding affinity and the structural stability. Importantly, we incorporate multiple strategies to faithfully estimate the statistical uncertainties in the FEP results. As a case study, we apply our protocols to systematically evaluate variants of the m396 antibody for their conformational stability and their binding affinity to the spike proteins of SARS-CoV-1 and SARS-CoV-2. By properly adjusting relevant parameters, the particle collapse problems in the FEP simulations are avoided. Furthermore, large statistical errors in a small fraction of the FEP calculations are effectively reduced by extending the sampling, such that acceptable statistical uncertainties are achieved for the vast majority of the cases with a modest total computational cost. Finally, our predicted conformational stability for the m396 variants is qualitatively consistent with the experimentally measured melting temperatures. Our work thus demonstrates the applicability of FEP in computational antibody design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory

Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)-an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Finally, analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration-physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.

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

Protein–Ligand Binding Free-Energy Calculations with ARROW–A Purely First-Principles Parameterized Polarizable Force Field

Protein–ligand binding free-energy calculations using molecular dynamics (MD) simulations have emerged as a powerful tool for in silico drug design. Here, we present results obtained with the ARROW force field (FF)–a multipolar polarizable and physics-based model with all parameters fitted entirely to high-level ab initio quantum mechanical (QM) calculations. ARROW has already proven its ability to determine solvation free energy of arbitrary neutral compounds with unprecedented accuracy. The ARROW FF parameterization is now extended to include coverage of all amino acids including charged groups, allowing molecular simulations of a series of protein–ligand systems and prediction of their relative binding free energies. We ensure adequate sampling by applying a novel technique that is based on coupling the Hamiltonian Replica exchange (HREX) with a conformation reservoir generated via potential softening and nonequilibrium MD. ARROW provides predictions with near chemical accuracy (mean absolute error of ~0.5 kcal/mol) for two of the three protein systems studied here (MCL1 and Thrombin). The third protein system (CDK2) reveals the difficulty in accurately describing dimer interaction energies involving polar and charged species. Overall, for all of the three protein systems studied here, ARROW FF predicts relative binding free energies of ligands with a similar accuracy level as leading nonpolarizable force fields.

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

Anisotropic Cahn-Hilliard free energy and interfacial energies for binary alloys with pairwise interactions

The original Cahn-Hilliard derivation of the contribution of compositional inhomogeneity to the free energy of a binary alloy with pairwise interactions is extended to include higher-order inhomogeneity terms. For alloys on a cubic lattice, the coefficient of the first inhomogeneity is a second-rank tensor and reduces to a scalar, but it is shown that the second order and the third order inhomogeneity terms are weighted by fourth-rank and sixth-rank tensors, thus resulting in anisotropic contributions. Furthermore, each interaction shell generates a unique set of inhomogeneity coefficients that is determined by the set of vectors connecting an atom to its neighbors on that shell. These coefficients are calculated for fcc and bcc alloys with interactions up fourth nearest neighbors. Phase field simulations based on these extended Cahn-Hilliard free energies are performed to measure interface free energies along specific crystallographic directions as a function of temperature, and to obtain the equilibrium shape of precipitates. Interface free energies, and the resulting anisotropies, are compared to those obtained by discrete models and Monte Carlo simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Studying the Collective Functional Response of a Receptor in Alchemical Ligand Binding Free Energy Simulations with Accelerated Solvation Layer Dynamics

Ligand binding free energy simulations (LB-FES) that involve sampling of protein functional conformations have been longstanding challenges in research on molecular recognition. Particularly, modeling of the conformational transition pathway and design of the heuristic biasing mechanism are severe bottlenecks for the existing enhanced configurational sampling (ECS) methods. Inspired by the key role of hydration in regulating conformational dynamics of macromolecules, this report proposes a novel ECS approach that facilitates binding-associated structural dynamics by accelerated hydration transitions in combination with the lambda-exchange of free energy perturbation (FEP). Further, two challenging protein-ligand binding processes involving large configurational transitions of the receptor are studied, with hydration transitions at binding sites accelerated by Hamiltonian-simulated annealing of the hydration layer. Without the need for pathway analysis or ad hoc barrier flattening potential, LB-FES were performed with FEP/lambda-exchange molecular dynamics simulation at a minor overhead for annealing of the hydration layer. The LB-FES studies showed that the accelerated rehydration significantly enhances the collective conformational transitions of the receptor, and convergence of binding affinity calculations is obtained at a sweet-spot simulation time scale. Alchemical LB-FES with the proposed ECS strategy is free from the effort of trial and error for the setup and realizes efficient on-the-fly sampling for the collective functional response of the receptor and bound water and therefore presents a practical approach to high-throughput screening in drug discovery.

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

Maximizing Free Energy Gain

Maximizing the amount of work harvested from an environment is important for a wide variety of biological and technological processes, from energy-harvesting processes such as photosynthesis to energy storage systems such as fuels and batteries. Here, we consider the maximization of free energy—and by extension, the maximum extractable work—that can be gained by a classical or quantum system that undergoes driving by its environment. We consider how the free energy gain depends on the initial state of the system while also accounting for the cost of preparing the system. We provide simple necessary and sufficient conditions for increasing the gain of free energy by varying the initial state. We also derive simple formulae that relate the free energy gained using the optimal initial state rather than another suboptimal initial state. Finally, we demonstrate that the problem of finding the optimal initial state may have two distinct regimes, one easy and one difficult, depending on the temperatures used for preparation and work extraction. We illustrate our results on a simple model of an information engine.

Physics↗
DOE OSTI2023

Free energy and specific heat near a quantum critical point of a metal

We analyze the free energy and the specific heat for fermions interacting with a gapless boson at a quantum-critical point (QCP) in a metal. We use the Luttinger-Ward-Eliashberg formula for the free energy in the normal state, which includes contributions from bosons, fermions, and their interaction, all expressed via fully dressed fermionic and bosonic propagators. The sum of the last two contributions is the free energy F γ of an effective low-energy model of fermions with boson-mediated dynamical 4-fermion interaction V(Ω m ) ∝ 1/|Ω m | γ (the γ model). This purely electronic model has been used to analyze the interplay between non-Fermi liquid behavior and pairing near a QCP, which are both independent of the upper energy cutoff Λ. However, the specific heat C γ (T), obtained from F γ , does depend on Λ. We argue that this dependence is spurious and cancels out, once we include the contribution from bosons. We further argue that the full C(T) is the sum of the contribution from free fermions and the one from a critical boson, with the fully dressed propagator, other terms cancel out. We compare the full C(T) with the C γ (T), obtained using recently proposed regularization of F γ . We argue that for γ < 1, the full C(T) and the regularized C γ (T) differ by a γ-dependent prefactor, while for γ > 1, the full C(T) and C γ (T) differ by the positive contribution from free massless fermions (a positive constant for the electron-phonon case γ = 2). For these γ, C γ (T) is negative, but the full C(T) is positive. We argue that only the full C(T) matters as the positive and the negative contributions originate from the term in C(T) which contains the fully dressed bosonic propagator. Here, we then argue that the normal state remains stable until the pairing instability develops.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1989

Free energy surfaces in the superconducting mixed state

The free energy surface for Tl2Ba2Ca2Cu3O1O has been measured as a function of temperature and magnetic field to determine the fundamental thermodynamic properties of the mixed state. The change in free energy, G(H)-G(O), is found to be linear in temperature over a wide range indicating that the specific heat is independent of field.

Finnemore, D. K.↗
DOE OSTI2025

Improving noisy free-energy measurements by adding more noise

Estimating free-energy differences using nonequilibrium work relations, such as the Jarzynski equality, is hindered by poor convergence when work fluctuations are large. For systems governed by overdamped Langevin dynamics, we propose the counterintuitive approach of adding noise in order to increase the precision of such calculations. Here, by introducing additional stochastic fluctuations to the system and rescaling its potential energy accordingly, we leave the thermodynamics of the system unchanged while increasing its relaxation rate. For a given time-dependent protocol this modification reduces the dissipated reduced work, leading to more accurate free-energy estimates. The method is designed to be used in experiment, and we illustrate its operation using computer simulations applied to two model systems. However, the regime of applicability of this strategy is likely limited, because it requires control of the system's potential energy in a way that is feasible in only a few experimental settings.

Whitelam, Stephen [Lawrence Berkeley National Labo↗
NASA NTRS2000

Apparent Minimum Free Energy Requirements for Methanogenic Archaea and Sulfate-Reducing Bacteria in an Anoxic Marine Sediment

Among the most fundamental constraints governing the distribution of microorganisms in the environment is the availability of chemical energy at biologically useful levels. To assess the minimum free energy yield that can support microbial metabolism in situ, we examined the thermodynamics of H2-consuming processes in anoxic sediments from Cape Lookout Bight, NC, USA. Depth distributions of H2 partial pressure, along with a suite of relevant concentration data, were determined in sediment cores collected in November (at 14.5 C) and August (at 27 C) and used to calculate free energy yields for methanogenesis and sulfate reduction. At both times of year, and for both processes, free energy yields gradually decreased (became less negative) with depth before reaching an apparent asymptote. Sulfate reducing bacteria exhibited an asymptote of -19.1 +/- 1.7 kj(mol SO4(2-)(sup -1) while methanogenic archaea were apparently supported by energy yields as small as -10.6 +/- 0.7 kj(mol CH4)(sup -1).

Hoehler, Tori M.↗
DOE OSTI2023

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

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

First principles free energy model with dynamic magnetism for δ -plutonium

We present an ab initio free energy model derived from a fully relativistic density functional theory (DFT) electronic structure with dynamic magnetism for δ -plutonium (face-centered cubic, fcc). The DFT model is extended with orbital-orbital interaction in a parameter free orbital polarization (OP) mechanism consistent with previous modeling of plutonium. Gibbs free energy is built from components associated with the temperature dependence of the electronic structure and the corresponding electronic entropy, lattice vibrations within an anharmonic lattice dynamics model, and dynamical fluctuations of the magnetization density, i.e. magnetic fluctuations. The fluctuation model consists of transverse and longitudinal modes driven by temperature induced excitations of the DFT + OP electronic structure. The ab initio model thus incorporates fluctuating states beyond the electronic ground state. Thanks to the dynamic magnetism, the theory predicts excellent thermodynamic properties and a Gibbs free energy in accord with CALPHAD and semi-empirical modeling developed from the thermodynamic observables. The magnetic fluctuations further explain anomalous behaviors of the thermal expansion in plutonium. Specifically, a thermal expansion for the δ -plutonium system turning from positive to negative at temperatures above room temperature, a tendency for gallium to reduce and remove the negative thermal expansion depending on composition, and a positive thermal expansion for the high temperature ϵ phase.

dynamic↗
NASA NTRS1983

Standard free energy of formation of iron iodide

An experiment is reported where silver iodide is used to determine the standard free energy of formation of iron iodide. By using silver iodide as a solid electrolyte, a galvanic cell, Ag/AgI/Fe-FeI2, is formulated. The standard free energy of formation of AgI is known, and hence it is possible to estimate the standard free energy of formation of FeI2 by measuring the open-circuit emf of the above cell as a function of temperature. The free standard energy of formation of FeI2 determined by this method is -38784 + 24.165T cal/mol. It is estimated that the maximum error associated with this method is plus or minus 2500 cal/mol.

Khandkar, A.↗
DOE OSTI2022

Ab initio free energies of liquid metal alloys: Application to the phase diagrams of Li-Na and K-Na

Comparison of free energies between different phases and different compositions underlies the prediction of alloy phase diagrams. To allow direct comparison, consistent reference points for the energies or enthalpies are required, and the entropy must be placed on an absolute scale, yielding absolute free energies. Here we derive absolute free energies of liquids from ab-initio molecular dynamics (AIMD) by combining the directly simulated enthalpies with an entropy derived from simulated densities and pair correlation functions. Additionally, as an example of the power of this method we calculate the phase diagrams of two binary alkali metal alloys, Li-Na and K-Na, revealing a critical point and liquid-liquid phase separation in the former case, and a deep eutectic in the latter. Good agreement with experimental data demonstrates the power of this simple method.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ligand structure and adsorption free energy of nanocrystals on solid substrates

In this paper, we present an investigation on the absorption of alkylthiolated nanocrystals on a solid substrate. We calculate adsorption free energies and report a number of effects induced by the substrate. Nearest neighbor distances and bonding free energies are significantly different than for a free floating case, there is a weakening of bonding free energies among nanocrystals, and the adsorption is manifestly anisotropic, i.e., stronger along certain directions of the nanocrystal core. We contend that this last result accounts for the Bain transition (fcc → bcc) observed in experimental results. We report the presence of vortices induced by the substrate, which explain the increased nearest neighbor distance among nanocrystals, which is in excellent quantitative agreement with experimental results and with the predictions of the Orbifold Topological Model. Implications for the assembly of nanostructures and future experiments are also discussed.

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

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↗
DOE OSTI2022

Extraction of the spin-glass free-energy landscape from $1/f$ noise measurements

The 1/f resistance noise has been measured in thin CuMn (13.5 at.%) spin-glass films. The temperature and frequency dependence have been analyzed in terms of the landscape of free-energy barriers. Here this analysis provides the full barrier distribution for various thin film thicknesses between 10 nm and 80 nm. The free-energy barrier height distribution's width and energy position have been determined. Contrary to previous models with fixed shape and energy, the free-energy landscape is described by a distribution of barriers that both shifts and changes shape as the temperature is reduced. The dependence of this distribution is in contrast with recent predictions. Using the fluctuation dissipation theorem, the 1/f measurements connect with χ", displaying agreement with other direct measurements of the latter.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Impact of magnetic and electric fields on the free energy to form a calcium carbonate ion-pair

Electromagnetic fields are used in water treatment and desalination to regulate scale formation and extend the lifetime of membranes. External electric and magnetic fields can promote or suppress mineral nucleation and growth. However, the molecular-scale mechanisms of such processes remain unknown. Computing the free energies needed to form ion pairs under external fields provides important insights into understanding the elemental steps during the initial formation of mineral scales. In this paper, we used molecular dynamics combined with metadynamics simulations to investigate the free energies of forming the [Ca–CO 3 ] 0 ion pair, a fundamental building block of carbonate scales, under a range of magnetic (up to 10 T) and electric (up to 10 V m −1 ) fields in water. The presence of constant magnetic or electric fields favored the ion pairing reaction and lowered the free energies by up to 3% to 6%. The internal energy and entropic components of the free energy showed significant changes and exhibited non-linear behavior with increasing field strength. The [Ca–CO 3 ] 0 ion pairing is an entropy-driven process in the absence of an external field, but the mechanism shifts to an internal energy-driven process under selected external fields, suggesting possible changes in the nucleation pathways.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗