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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Relationship between wave energy and free energy from pickup ions in the Comet Halley environment

The free energy available from the implanted heavy ion population at Comet Halley is calculated by assuming that the initial unstable velocity space ring distribution of the ions evolves toward a bispherical shell. Ultimately this free energy adds to the turbulence in the solar wind. Upstream and downstream free energies are obtained separately for the conditions observed along the Giotto spacecraft trajectory. The results indicate that the waves are mostly upstream propagating in the solar wind frame. The total free energy density always exceeds the measured wave energy density because, as expected in the nonlinear process of ion scattering, the available energy is not all immediately released. An estimate of the amount which has been released can be obtained from the measured oxygen ion distributions and again it exceeds that observed. The theoretical analysis is extended to calculate the k spectrum of the cometary-ion-generated turbulence.

Huddleston, D. E.↗
DOE OSTI2025

Free-energy orbital-free density functional theory: recent developments, perspective, and outlook

By summarizing the constraint-based development of orbital-free free-energy density functional approximations, we provide a perspective on progress over the last 15 years, the limitations of existing functionals, and the challenges awaiting resolution. We outline the chronology of the development of noninteracting and exchange-correlation free-energy orbital-free functionals and summarize the theoretical basis of existing local density approximation, second-order approximation, generalized gradient approximation (GGA), and meta-GGAs. We discuss limitations and challenges such as problems with thermodynamic derivatives, free-energy nonadditivity and the closely related issue of all-electron versus valence-only local pseudo-potential performance.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

CHARMM-GUI Free Energy Calculator for Absolute and Relative Ligand Solvation and Binding Free Energy Simulations

Alchemical free energy simulations have long been utilized to predict free energy changes for binding affinity and solubility of small molecules. However, while the theoretical foundation of these methods is well established, seamlessly handling many of the practical aspects regarding the preparation of the different thermodynamic end states of complex molecular systems and the numerous processing scripts often remains a burden for successful applications. Here, we present CHARMM-GUI Free Energy Calculator (http://www.charmm-gui.org/input/fec) that provides various alchemical free energy perturbation molecular dynamics (FEP/MD) systems with input and post-processing scripts for NAMD and GENESIS. Four submodules are available: Absolute Ligand Binder (for absolute ligand binding FEP/MD), Relative Ligand Binder (for relative ligand binding FEP/MD), Absolute Ligand Solvator (for absolute ligand solvation FEP/MD), and Relative Ligand Solvator (for relative ligand solvation FEP/MD). Each module is designed to build multiple systems of a set of selected ligands at once for high-throughput FEP/MD simulations. The capability of Free Energy Calculator is illustrated by absolute and relative solvation FEP/MD of a set of ligands and absolute and relative binding FEP/MD of a set of ligands for T4-lysozyme in solution and the adenosine A 2A receptor in a membrane. The calculated free energy values are overall consistent with the experimental and published free energy results (within ~1 kcal/mol). We hope that Free Energy Calculator is useful to carry out high-throughput FEP/MD simulations in the field of biomolecular sciences and drug discovery.

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

Calculated solvent reorganization entropies, free energies, and fluctuations of the energy gaps for intramolecular electron transfer and excitation of the solvatochromic dye B30

Intramolecular electron transfer between two biphenyl groups linked by an androstane spacer and excitation of the pyridinium-N-phenolate betaine dye B30 to the first excited singlet state are studied by quantum/classical molecular-dynamics simulations at temperatures between 150 and 300 K in solvents with a range of polarities. Temperature dependences of the solvent reorganization energies, free energies, entropies, and the inhomogeneous broadening of B30’s absorption band are examined. The variances of fluctuations of the energy gap between the reactant and product states do not have the direct proportionality to temperature that often is assumed to hold. An explanation for the observations is suggested.

Chemistry↗
DOE OSTI2020

Combined Force-Frequency Sampling for Simulation of Systems Having Rugged Free Energy Landscapes.

An adaptive, machine learning-based sampling method is presented for simulation of systems having rugged, multidimensional free energy landscapes. The method's main strength resides in its ability to learn both from the frequency of visits to distinct states and the generalized force estimates that arise in a system as it evolves in phase space. This is accomplished by introducing a self-integrating artificial neural network, which generates an estimate of the free energy directly from its derivatives. The usefulness of the proposed combined approach is examined in the context of two concrete examples, namely, an alanine dipeptide molecule in water and a polymer diffusing through a narrow pore. This new method is found to be robust, faster, and more accurate than approaches that rely only on frequency-based or generalized force-based estimations. After combining the proposed approach with overfill protection and support for sparse data storage and training, the method is shown to be more effective than comparable, previously available techniques and capable of scaling efficiently to larger numbers of collective variables.

Sevgen, Emre↗
DOE OSTI2021

Molecular Dynamics Free Energy Simulations Reveal the Mechanism for the Antiviral Resistance of the M66I HIV-1 Capsid Mutation

While drug resistance mutations can often be attributed to the loss of direct or solvent-mediated protein−ligand interactions in the drug-mutant complex, in this study we show that a resistance mutation for the picomolar HIV-1 capsid (CA)-targeting antiviral (GS-6207) is mainly due to the free energy cost of the drug-induced protein side chain reorganization in the mutant protein. Among several mutations, M66I causes the most suppression of the GS-6207 antiviral activity (up to ~84,000-fold), and only 83- and 68-fold reductions for PF74 and ZW-1261, respectively. To understand the molecular basis of this drug resistance, we conducted molecular dynamics free energy simulations to study the structures, energetics, and conformational free energy landscapes involved in the inhibitors binding at the interface of two CA monomers. To minimize the protein−ligand steric clash, the I66 side chain in the M66I−GS-6207 complex switches to a higher free energy conformation from the one adopted in the apo M66I. In contrast, the binding of GS-6207 to the wild-type CA does not lead to any significant M66 conformational change. Based on an analysis that decomposes the absolute binding free energy into contributions from two receptor conformational states, it appears that it is the free energy cost of side chain reorganization rather than the reduced protein−ligand interaction that is largely responsible for the drug resistance against GS-6207.

Sun, Qinfang (ORCID:0000000173621543)↗
DOE OSTI2018

Free-Energy Landscape of the Dissolution of Gibbsite at High pH

The individual elementary reactions involved in the dissolution of a solid into solution remain mostly speculative due to a lack of direct experimental probes. In this regard, in this paper we have applied atomistic simulations to map the free-energy landscape of the dissolution of gibbsite from a step edge as a model of metal hydroxide dissolution. The overall reaction combines kink formation and kink propagation. Two individual reactions were found to be rate-limiting for kink formation, that is, the displacement of Al from a step site to a ledge adatom site and its detachment from ledge/terrace adatom sites into the solution. Finally, as a result, a pool of mobile and labile adsorbed species, or adatoms, exists before the release of Al into solution. Because of the quasi-hexagonal symmetry of gibbsite, kink site propagation can occur in multiple directions. Overall, our results will enable the development of microscopic mechanistic models of metal oxide dissolution.

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

Molecular Level Understanding of the Free Energy Landscape in Early Stages of Metal-Organic Framework Nucleation

The assembly mechanism of hierarchical materials controlled by the choice of solvent and presence of spectator ions. In this paper, we use enhanced sampling molecular dynamics methods to investigate these effects on the configurational landscape of metal-linker interactions in the early stages of synthesis, using MIL-101(Cr) as a prototypical example. Microsecond-long well-tempered metadynamics (WTmetaD) uncover a complex free energy structural landscape, with distinct crystal (C) and non-crystal (NC) like configurations and their equilibrium population. In presence of ions (Na+, F-), we observe a complex effect on the crystallinity of secondary building units (SBUs), by encouraging/suppressing salt bridges between C configurations and consequently controlling the percentage of defects. Solvent effects are assessed by introducing N, N-dimethylformamide (DMF) instead of water, where SBU adducts are appreciably more stable and compact. These results shed light on how solvent nature and ionic strength impact the free energy of nucleation phenomena that ultimately control materials synthesis and defect formation. The work described in this publication was supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences and performed at Pacific Northwest National Laboratory (PNNL). The authors acknowledge the use of the UCL Legion High Performance Computing Facility (Legion @ UCL) and associated support services in the completion of this work. This research used resources of the National Energy Research Scientific Computing Centre, a DOE Office of Science User Facility. Battelle operates PNNL for DOE under Contract DE-AC05-76RL01830.

Kollias, Loukas↗
DOE OSTI2021

Toward Automated Sampling of Polymorph Nucleation and Free Energies with the SGOOP and Metadynamics

Understanding the driving forces behind the nucleation of different polymorphs is of great importance for material sciences and the pharmaceutical industry. This includes understanding the reaction coordinate that governs the nucleation process and correctly calculating the relative free energies of different polymorphs. Here, we demonstrate, for the prototypical case of urea nucleation from the melt, how one can learn such a one-dimensional reaction coordinate as a function of prespecified order parameters and use it to perform efficient biased all-atom molecular dynamics simulations. The reaction coordinate is learnt as a function of the generic thermodynamic and structural order parameters using the “spectral gap optimization of order parameters (SGOOP)” approach [Tiwary, P. and Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. (2016)] and is biased using well-tempered metadynamics simulations. The reaction coordinate gives insights into the role played by different structural and thermodynamics order parameters, and the biased simulations obtain accurate relative free energies for different polymorphs. This includes an accurate prediction of the approximate pressure at which urea undergoes a phase transition and one of the metastable polymorphs becomes the most stable conformation. Here, we believe the ideas demonstrated in this work will facilitate efficient sampling of nucleation in complex, generic systems.

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

Adsorbate Free Energies from DFT-Derived Translational Energy Landscapes

Adsorption free energies are fundamental to surface chemistry and catalysis. Standard models combine some assumed analytical form of the translational potential energy surface, often parameterized against density functional theory (DFT) calculations, with an analytical expression for the resultant translational densities of states (DOS), free energy, and entropy. Here we compare the performance of such models against numerical evaluations of the DOS and thermodynamic functions derived from solutions to the translational Schrodinger equation. We compare results for a translational potential energy surface (PES) derived from nudged eleastic band calculations with those obtained from adsorbate rastering across a series of monatomic (O, S, C, N, H) and polyatomic (NH x ) adsorbates on (100) Pt and Au facets. We find that analytical models as commonly parameterized have mixed performance for describing the translational PES and that the consequences for computed free energies are modest but potentially significant in microkinetic models. Here, numerical solutions are possible for modest to no additional computational cost over analytical models and thus should be considered when reliable free energy estimates are needed or translational potential energy surfaces are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Temperature-dependent solvent reorganization entropies, free energies, and transition dipole strengths for the photoexcitation of Reichardt’s dye B30

Absorption spectra of the solvatochromic dye 2,6-diphenyl-4-2,4,6-triphenyl-1-pyridinophenolate (B30) were measured in seven solvents of varying polarity over temperature ranging from each solvent’s freezing point to 300 K. The excitation energies and their variances allowed calculations of the solvent reorganization energies, reorganization free energies and reorganization entropies as functions of temperature. The entropies of solvent packing around the chromophore are found to make major contributions to the reorganization free energies. The variances of the excitation energies depend only weakly on temperature, in disagreement with an expression that is often used for solvent reorganization free energies. Polar solvents reduce the transition dipole strength of B30’s long-wavelength absorption band, probably because interactions with the solvent enhance the charge-transfer character of the transition. The dipole strength drops further at low temperatures.

Chemistry↗
NASA NTRS2012

Non-Equilibrium Properties from Equilibrium Free Energy Calculations

Calculating free energy in computer simulations is of central importance in statistical mechanics of condensed media and its applications to chemistry and biology not only because it is the most comprehensive and informative quantity that characterizes the eqUilibrium state, but also because it often provides an efficient route to access dynamic and kinetic properties of a system. Most of applications of equilibrium free energy calculations to non-equilibrium processes rely on a description in which a molecule or an ion diffuses in the potential of mean force. In general case this description is a simplification, but it might be satisfactorily accurate in many instances of practical interest. This hypothesis has been tested in the example of the electrodiffusion equation . Conductance of model ion channels has been calculated directly through counting the number of ion crossing events observed during long molecular dynamics simulations and has been compared with the conductance obtained from solving the generalized Nernst-Plank equation. It has been shown that under relatively modest conditions the agreement between these two approaches is excellent, thus demonstrating the assumptions underlying the diffusion equation are fulfilled. Under these conditions the electrodiffusion equation provides an efficient approach to calculating the full voltage-current dependence routinely measured in electrophysiological experiments.

Pohorille, Andrew↗