Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

Hybrid Quantum Mechanical, Molecular Mechanical, and Machine Learning Potential for Computing Aqueous-Phase Adsorption Free Energies on Metal Surfaces

Performing reliable computer simulations of elementary processes occurring at metal–water interfaces is pivotal for novel catalyst design in sustainable energy applications. Computational catalyst design hinges on the ability to reliably and efficiently compute the potential energy surface (PES) of the system. Here, due to the large system sizes needed for studying processes at liquid water–metal interfaces, these systems can currently not be described using density functional theory (DFT). In this work, we used a hybrid quantum mechanical, molecular mechanical, and machine learning potential for studying the adsorption behavior of phenol, atomic hydrogen, 2-butanol, and 2-butanone on the (0001) facet of Ru under reducing conditions when Ru is not oxidized. Specifically, we describe the adsorbate and the surrounding metal atoms at the DFT level of theory. Here, we also considered the electrostatic field effect of the water molecules on adsorbate–metal interactions. Next, for the water–water and water–adsorbate interactions, we used established classical force fields. Finally, for the water–Ru surface interaction, for which no reliable force fields have been published, we used Behler–Parrinello high-dimensional neural network potentials (HDNNPs). Employing this setup, we used our explicit solvation for metal surface (eSMS) approach to compute the aqueous-phase effect on the low-coverage adsorption of selected molecules and atoms on the (0001) facet of Ru. In agreement with previous experimental and computational studies of oxygenated molecules over transition metal facets, we found that liquid water destabilizes the tested adsorbates on Ru(0001). Interestingly, our findings indicate that adsorbates on Ru are less affected by the presence of an aqueous phase than on other transition metals (e.g., Pt), highlighting the necessity of experimental investigations of Ru-based catalytic systems in liquid water.

Adsorption↗
NASA NTRS1996

Approximate scaling properties of RNA free energy landscapes

RNA free energy landscapes are analysed by means of "time-series" that are obtained from random walks restricted to excursion sets. The power spectra, the scaling of the jump size distribution, and the scaling of the curve length measured with different yard stick lengths are used to describe the structure of these "time series". Although they are stationary by construction, we find that their local behavior is consistent with both AR(1) and self-affine processes. Random walks confined to excursion sets (i.e., with the restriction that the fitness value exceeds a certain threshold at each step) exhibit essentially the same statistics as free random walks. We find that an AR(1) time series is in general approximately self-affine on timescales up to approximately the correlation length. We present an empirical relation between the correlation parameter rho of the AR(1) model and the exponents characterizing self-affinity.

Signal Processing, Computer-Assisted↗
NASA NTRS1977

Variational energy principle for compressible, baroclinic flow. 2: Free-energy form of Hamilton's principle

The first and second variations are calculated for the irreducible form of Hamilton's Principle that involves the minimum number of dependent variables necessary to describe the kinetmatics and thermodynamics of inviscid, compressible, baroclinic flow in a specified gravitational field. The form of the second variation shows that, in the neighborhood of a stationary point that corresponds to physically stable flow, the action integral is a complex saddle surface in parameter space. There exists a form of Hamilton's Principle for which a direct solution of a flow problem is possible. This second form is related to the first by a Friedrichs transformation of the thermodynamic variables. This introduces an extra dependent variable, but the first and second variations are shown to have direct physical significance, namely they are equal to the free energy of fluctuations about the equilibrium flow that satisfies the equations of motion. If this equilibrium flow is physically stable, and if a very weak second order integral constraint on the correlation between the fluctuations of otherwise independent variables is satisfied, then the second variation of the action integral for this free energy form of Hamilton's Principle is positive-definite, so the action integral is a minimum, and can serve as the basis for a direct trail and error solution. The second order integral constraint states that the unavailable energy must be maximum at equilibrium, i.e. the fluctuations must be so correlated as to produce a second order decrease in the total unavailable energy.

Schmid, L. A.↗
NASA NTRS1985

Statistical mechanics of light elements at high pressure. VII - A perturbative free energy for arbitrary mixtures of H and He

A model free energy is presented which accurately represents results from 45 high-precision Monte Carlo calculations of the thermodynamics of hydrogen-helium mixtures at pressures of astrophysical and planetophysical interest. The free energy is calculated using free-electron perturbation theory (dielectric function theory), and is an extension of the expression given in an earlier paper in this series. However, it fits the Monte Carlo results more accurately, and is valid for the full range of compositions from pure hydrogen to pure helium. Using the new free energy, the phase diagram of mixtures of liquid metallic hydrogen and helium is calculated and compared with earlier results. Sample results for mixing volumes are also presented, and the new free energy expression is used to compute a theoretical Jovian adiabat and compare the adiabat with results from three-dimensional Thomas-Fermi-Dirac theory. The present theory gives slightly higher densities at pressures of about 10 megabars.

Hubbard, W. B.↗
NASA NTRS2013

Good Practices in Free-energy Calculations

As access to computational resources continues to increase, free-energy calculations have emerged as a powerful tool that can play a predictive role in drug design. Yet, in a number of instances, the reliability of these calculations can be improved significantly if a number of precepts, or good practices are followed. For the most part, the theory upon which these good practices rely has been known for many years, but often overlooked, or simply ignored. In other cases, the theoretical developments are too recent for their potential to be fully grasped and merged into popular platforms for the computation of free-energy differences. The current best practices for carrying out free-energy calculations will be reviewed demonstrating that, at little to no additional cost, free-energy estimates could be markedly improved and bounded by meaningful error estimates. In energy perturbation and nonequilibrium work methods, monitoring the probability distributions that underlie the transformation between the states of interest, performing the calculation bidirectionally, stratifying the reaction pathway and choosing the most appropriate paradigms and algorithms for transforming between states offer significant gains in both accuracy and precision. In thermodynamic integration and probability distribution (histogramming) methods, properly designed adaptive techniques yield nearly uniform sampling of the relevant degrees of freedom and, by doing so, could markedly improve efficiency and accuracy of free energy calculations without incurring any additional computational expense.

Pohorille, Andrew↗
DOE OSTI2025

Nonadiabatic Force Matching for Alchemical Free-Energy Estimation

We propose a method to compute free-energy differences from nonadiabatic alchemical transformations by using flow-based generative models. The method, nonadiabatic force matching, hinges on estimating the dissipation along an alchemical switching process in terms of a nonadiabatic force field that can be learned through stochastic flow matching. The learned field can be used in conjunction with short-time trajectory data to evaluate upper and lower bounds on the alchemical free energy that variationally converge to the exact value if the field is optimal. Applying the method to evaluate the alchemical free energy of atomistic models shows that it can substantially reduce the simulation cost of a free-energy estimate at a negligible loss of accuracy when compared with thermodynamic integration.

Computational chemistry↗
NASA NTRS1994

Direct measurement of free-energy barrier to nucleation of crystallites in amorphous silicon thin films

A method is introduced to measure the free-energy barrier W(sup *), the activation energy, and activation entropy to nucleation of crystallites in amorphous solids, independent of the energy barrier to growth. The method allows one to determine the temperature dependence of W(sup *), and the effect of the preparation conditions of the initial amorphous phase, the dopants, and the crystallization methds on W(sup *). The method is applied to determine the free-energy barrier to nucleation of crystallites in amorphous silicon (a-Si) thin films. For thermally induced nucleation in a-Si thin films with annealing temperatures in the range of from 824 to 983 K, the free-energy barrier W(sup *) to nucleation of silicon crystals is about 2.0 - 2.1 eV regardless of the preparation conditions of the films. The observation supports the idea that a-Si transforms into an intermediate amorphous state through the structural relaxation prior to the onset of nucleation of crystallites in a-Si. The observation also indicates that the activation entropy may be an insignificant part of the free-energy barrier for the nucleation of crystallites in a-Si. Compared with the free-energy barrier to nucleation of crystallites in undoped a-Si films, a significant reduction is observed in the free-energy barrier to nucleation in Cu-doped a-Si films. For a-Si under irradiation of Xe(2+) at 10(exp 5) eV, the free-energy barrier to ion-induced nucleation of crystallites is shown to be about half of the value associated with thermal-induced nucleation of crystallites in a-Si under the otherwise same conditions, which is much more significant than previously expected. The present method has a general kinetic basis; it thus should be equally applicable to nucleation of crystallites in any amorphous elemental semiconductors and semiconductor alloys, metallic and polymeric glasses, and to nucleation of crystallites in melts and solutions.

Shi, Frank G.↗
NASA NTRS2019

Determining the Transport of Magnetic Helicity and Free Energy in the Sun's Atmosphere

The most important factors determining solar coronal activity are believed to be the availability of magnetic free energy and the constraint of magnetic helicity conservation. Direct measurements of the helicity and magnetic free energy in the coronal volume are difficult, but their values may be estimated from measurements of the helicity and free energy transport rates through the photosphere. We examine these transport rates for a topologically open system such as the corona, in which the magnetic fields have a nonzero normal component at the boundaries, and derive a new formula for the helicity transport rate at the boundaries. In addition, we derive new expressions for helicity transport due to flux emergence/submergence versus photospheric horizontal motions. The key feature o four formulas is that they are manifestly gauge invariant. Our results are somewhat counterintuitive in that only the lamellar electric field produced by the surface potential transports helicity across boundaries, and the solenoidal electric field produced by a surface stream function does not contribute to the helicity transport. We discuss the physical interpretation of this result. Furthermore, we derive an expression for the free energy transport rate and show that a necessary condition for free energy transport across a boundary is the presence of a closed magnetic field at the surface, indicating that there are current systems within the volume. We discuss the implications of these results for using photospheric vector magnetic and velocity field measurements to derive the solar coronal helicity and magnetic free energy, which can then be used to constrain and drive models for coronal activity.

Schuck, Peter W.↗
NASA NTRS2001

Calculating Free Energies Using Average Force

A new, general formula that connects the derivatives of the free energy along the selected, generalized coordinates of the system with the instantaneous force acting on these coordinates is derived. The instantaneous force is defined as the force acting on the coordinate of interest so that when it is subtracted from the equations of motion the acceleration along this coordinate is zero. The formula applies to simulations in which the selected coordinates are either unconstrained or constrained to fixed values. It is shown that in the latter case the formula reduces to the expression previously derived by den Otter and Briels. If simulations are carried out without constraining the coordinates of interest, the formula leads to a new method for calculating the free energy changes along these coordinates. This method is tested in two examples - rotation around the C-C bond of 1,2-dichloroethane immersed in water and transfer of fluoromethane across the water-hexane interface. The calculated free energies are compared with those obtained by two commonly used methods. One of them relies on determining the probability density function of finding the system at different values of the selected coordinate and the other requires calculating the average force at discrete locations along this coordinate in a series of constrained simulations. The free energies calculated by these three methods are in excellent agreement. The relative advantages of each method are discussed.

Darve, Eric↗
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↗
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.↗