Search NASA⌕ Search

SEARCH · Search NASA

Results for “Free energy calculations”

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 conditions for accurate modeling of urea production via co-electrolysis

Co-electrolysis of carbon oxides and nitrogen oxides promise to simultaneously help restore the balance of the C and N cycles while producing valuable chemicals such as urea. However, co-electrolysis processes are still largely inefficient and numerous knowledge voids persist. Here, we provide a solid thermodynamic basis for modelling urea production via co-electrolysis. First, we determine the energetics of aqueous urea produced under electrochemical conditions based on experimental data, which enables an accurate assessment of equilibrium potentials and overpotentials. Next, we use density functional theory (DFT) calculations to model various co-electrolysis reactions producing urea. The calculated reaction free energies deviate significantly from experimental values for well-known GGA, meta-GGA and hybrid functionals. These deviations stem from errors in the DFT-calculated energies of molecular reactants and products. In particular, the error for urea is approximately -0.25 ± 0.10 eV. Finally, we show that all these errors introduce large inconsistencies in the calculated free-energy diagrams of urea production via co-electrolysis, such that gas-phase corrections are strongly advised.

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

Mechanistic Insights on Permeation of Water over Iron Cations in Nanoporous Silicon Oxide Films for Selective H 2 and O 2 Evolution

Electrocatalysts encapsulated by an ultrathin and semipermeable oxide layer offer a promising avenue for efficient, selective, and cost-effective production of hydrogen through photoelectrochemical water splitting. This architecture is especially attractive for Z-scheme water splitting, for which a nanoporous oxide film can be leveraged to mitigate undesired, yet kinetically facile, reactions involving redox shuttles, such as aqueous iron cations, by limiting transport of these species to catalytically active sites. In this work, molecular dynamics simulations were combined with electrochemical measurements to provide a mechanistic understanding of permeation of water and Fe(III)/Fe(II) redox shuttles through nanoporous SiO 2 films. It is shown that even for SiO 2 pores with a width as small as 0.8 nm, water does not experience any energy barrier for permeating into the pores due to a favorable interaction with hydrophilic silanol groups on the oxide surface. In contrast, permeation of Fe(III) and Fe(II) into microporous SiO 2 pores is limited due to high energy barriers, which stem from a combination of distortion and dehydration of the second and third ion solvation shells. Furthermore, our simulations and experimental results show that SiO 2 coatings can effectively mitigate undesired Fe(III)/Fe(II) redox reactions at underlying electrodes by attenuating permeation of iron cations, while allowing water to permeate and thus participate in water splitting reactions. In a broader context, our study demonstrates that selectivity of solvated cations can be manipulated by controlling the pore size and surface chemistry of oxide films.

14 SOLAR ENERGY↗
NASA NTRS1996

Redox Disproportionation of Glucose as a Major Biosynthetic Energy Source

Previous studies have concluded that very little if any energy is required for the microbial biosynthesis of amino acids and lipids from glucose -- processes that yield almost as much ATP (adenosine triphosphate) as they consume. However, these studies did not establish the strength nor the nature of the energy source driving these biological transformations. To identify and estimate the strength of the energy source behind these processes, we calculated the free energy change due to the redox disproportionation of substrate carbon of (a) 26 redox-balanced fermentation reactions, and (b) the biosynthesis of amino acids, lipids, and nucleotides of E. coli from glucose. A plot of the negative free energy of these reactions per mmole of carbon as a function of the number of disproportionative electron transfers per mmol of carbon showed that the energy yields of these fermentations and biosyntheses were directly proportional to the degree of redox disproportionation of carbon. Since this linear relationship showed that redox disproportionation was the dominant energy source of these reactions, we were able to establish that amino acid and lipid biosynthesis obtained most of their energy from redox disproportionation (greater than 94%). In contrast nucleotide biosynthesis was not driven by redox disproportionation of carbon, and consequently depended completely on ATP for energy. This crucial and previously unrecognized role of sugars as an energy source of biosynthesis suggests that sugars were involved at the earliest stage in the origin of anabolic metabolism.

Weber, Arthur L.↗
DOE OSTI2024

Thermodynamics of Tritium Trapping by Point Defects in Intermetallic Al 12 (TM) 2.35 Aluminide Coating Phase

Density functional theory simulations have been carried out to investigate the potential for tritium trapping by metal vacancies in intermetallic Al 12 (TM) 2.35 phase (TM = Fe, Cr, and Ni) as function of temperature and tritium partial pressure. It was found that tritium could be favorably trapped by Fe and Ni vacancies and not favorably trapped by Al and Cr vacancies. However, due to the presence of partially occupied Al sites in bulk Al 12 (TM) 2.35 , leading to the approximate number of ~255 Al atoms in the unit cell, 86 sites were found energetically favorable to the creation of an Al vacancy. While adding a tritium atom in an Al vacancy is not energetically favorable, the tritiated defect still has a negative Gibbs free energy because the energy gain for creating an Al vacancy overcome the energy cost of adding the tritium species. Based on the calculated Gibbs free energy, the first tritiation of a metal vacancy, at conditions relevant to in-reactor operations, should be more favorable for Al, followed Fe, Ni, and Cr vacancies. By comparing the behavior of tritium in Al 12 (TM) 2.35 with previously studied Fe-Al coating phases (i.e., FeNiAl 5 , Fe 4 Al 13 , and Fe 2 Al 5.6 ), we found that there is a correlation between interstitial tritium solubility and the potential for vacancy trapping. The current trend suggests that if the insertion of an interstitial tritium cost more than 0.3 eV, then trapping by metal vacancies should be preferred. By combining the simulations results obtained to date, we noticed different trapping mechanisms of tritium in the Al coating. Tritium is mostly trapped by Fe and Ni vacancies in the outer Fe-Al coating phase Al 12 (TM) 2.35 while tritium should be preferentially trapped by Al and Fe vacancies for the inner Fe-Al coating phases (FeNiAl 5 , Fe 4 Al 13 , Fe 2 Al 5.6 ). Altogether, these studies show that tritium interacts differently with the various Fe-Al aluminide phases, they also suggest that tritium trapping and retention could be more efficient if metal defects are present and if the solubility of interstitial tritium in the different phases is low.

36 MATERIALS SCIENCE↗
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

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 OSTI2023

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1997

Energy from Redox Disproportionation of Sugar Carbon Drives Biotic and Abiotic Synthesis

To identify the energy source that drives the biosynthesis of amino acids, lipids, and nucleotides from glucose, we calculated the free energy change due to redox disproportionation of the substrate carbon of: (1) 26-carbon fermentation reactions and (2) the biosynthesis of amino acids and lipids of E. coli from glucose. The free energy (cal/mmol of carbon) of these reactions was plotted as a function of the degree of redox disproportionation of carbon (disproportionative electron transfers (mmol)/mmol of carbon). The zero intercept and proportionality between energy yield and degree of redox disproportionation exhibited by this plot demonstrate that redox disproportionation is the principal energy source of these redox reactions (slope of linear fit = -10.4 cal/mmol of disproportionative electron transfers). The energy and disproportionation values of E. coli amino acid and lipid biosynthesis from glucose lie near this linear curve fit with redox disproportionation accounting for 84% and 96% (and ATP only 6% and 1 %) of the total energy of amino acid and lipid biosynthesis, respectively. These observations establish that redox disproportionation of carbon, and not ATP, is the primary energy source driving amino acid and lipid biosynthesis from glucose. In contrast, we found that nucteotide biosynthesis involves very little redox disproportionation, and consequently depends almost entirely on ATP for energy. The function of sugar redox disproportionation as the major source of free energy for the biosynthesis of amino acids and lipids suggests that sugar disproportionation played a central role in the origin of metabolism, and probably the origin of life.

Weber, Arthur L.↗
NASA NTRS1997

Energy from Redox Disproportionation of Sugar Carbon Drives Biotic and Abiotic Synthesis

To identify the energy source that drives the biosynthesis of amino acids, lipids, and nucleotides from glucose, we calculated the free energy change due to redox disproportionation of the substrate carbon of: (1) 26-carbon fermentation reactions and (2) the biosynthesis of amino acids and lipids of E. coli from glucose. The free energy (cal/mmol of carbon) of these reactions was plotted as a function of the degree of redox disproportionation of carbon (disproportionative electron transfers (mmol)/mmol of carbon). The zero intercept and proportionality between energy yield and degree of redox disproportionation exhibited by this plot demonstrate that redox disproportionation is the principal energy source of these redox reactions (slope of linear fit = -10.4 cal/mmol of disproportionative electron transfers). The energy and disproportionation values of E. coli amino acid and lipid biosynthesis from glucose lie near this linear curve fit with redox disproportionation accounting for 84% and 96% (and ATP only 6% and 1%) of the total energy of amino acid and lipid biosynthesis, respectively. These observations establish that redox disproportionation of carbon, and not ATP, is the primary energy source driving amino acid and lipid biosynthesis from glucose. In contrast, we found that nucleotide biosynthesis involves very little redox disproportionation, and consequently depends almost entirely on ATP for energy. The function of sugar redox disproportionation as the major source of free energy for the biosynthesis of amino acids and lipids suggests that sugar disproportionation played a central role in the origin of metabolism, and probably the origin of life.

Weber, Arthur L.↗
DOE OSTI2021

Entropy Pair Functional Theory: Direct Entropy Evaluation Spanning Phase Transitions

We prove that, within the class of pair potential Hamiltonians, the excess entropy is a universal, temperature-independent functional of the density and pair correlation function. This result extends Henderson’s theorem, which states that the free energy is a temperature dependent functional of the density and pair correlation. The stationarity and concavity of the excess entropy functional are discussed and related to the Gibbs–Bugoliubov inequality and to the free energy. We apply the Kirkwood approximation, which is commonly used for fluids, to both fluids and solids. Approximate excess entropy functionals are developed and compared to results from thermodynamic integration. The pair functional approach gives the absolute entropy and free energy based on simulation output at a single temperature without thermodynamic integration. We argue that a functional of the type, which is strictly applicable to pair potentials, is also suitable for first principles calculation of free energies from Born–Oppenheimer molecular dynamics performed at a single temperature. In conclusion, this advancement has the potential to reduce the evaluation the free energy to a simple modification to any procedure that evaluates the energy and the pair correlation function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Deriving reliable nucleation rates from metadynamics simulations: Application to Yukawa fluids

In order to solidify the usefulness of metadynamics in studying nucleation of crystals from supercooled liquids, this work provides a specific procedure to calculate nucleation free energy barriers. After a pedagogical review of the important elements of classical nucleation theory and how metadynamics is used to find nucleation free energy barriers, we explain the benefits of local collective variables over more common global collective variables. We show how a metadynamics free energy barrier must be carefully post-processed so that classical nucleation theory can be applied to calculate nucleation rates. We apply our procedure to a Yukawa plasma and show that a particular physically motivated fit to metadynamics data reproduces low-temperature reference data, justifying the usefulness of metadynamics to predict nucleation rates.

42 ENGINEERING↗
DOE OSTI2024

Thermochemical Data for Furan-based Monomer Candidates for Frontal Ring-Opening Metathesis Polymerization (FROMP)

This dataset includes 471 furan-based monomer candidates for frontal ring-opening metathesis polymerization (FROMP) and relevant thermochemistry as calculated with density functional theory (DFT). The monomer candidates were combinatorically enumerated using Diels-Alder reactions of furan derivatives as dienes and four types of dienophiles (alkenes, alkynes, allenes, and benzynes). Common substituents were enumerated for the dienophile classes, and methyl substitution on the diene was explored. We used the SMILES arbitrary target specification (SMARTS) language to produce monomers and ring-opened structures from diene and dienophile precursor SMILES, and we studied the ring-opening reaction using a homodesmotic equation with ethene. RDKit conformers were initially generated from SMILES, then optimized with GFN2-xTB. The two conformers lowest in energy were then optimized with DFT using the wb97x-D3 functional, def2-TZVP basis set, and def2/J auxiliary basis set. Gibbs free energy corrections were obtained through frequency calculations. Structures with imaginary frequencies below -50 cm^{-1} were excluded from this work, and smaller imaginary modes were flipped to be positive for free energy calculations. Modes below 50 cm^{-1} were treated with the modified rigid rotor approximation, and all thermochemical values were calculated at T=200C. The CSV file contains the monomer SMILES, the free energy of reaction for Diels-Alder addition (G_DA_200), and the enthalpy of the ring-opening reaction (H_RO_200). All energies are given in kcal/mol. An interactive HTML is also included to visualize the monomers in this dataset.

Chua, Lauren↗
DOE OSTI2021

Neural Network Sampling of the Free Energy Landscape for Nitrogen Dissociation on Ruthenium

In heterogeneous catalysis, free energy profiles of reactions govern the mechanisms, rates, and equilibria. Energetics are conventionally computed using the harmonic approximation (HA), which requires determination of critical states a priori. Here, we use neural networks to efficiently sample and directly calculate the free energy surface (FES) of a prototypical heterogeneous catalysis reaction—the dissociation of molecular nitrogen on ruthenium—at density-functional-theory-level accuracy. We find that the vibrational entropy of surface atoms, often neglected in HA for transition metal catalysts, contributes significantly to the reaction barrier. The minimum free energy path for dissociation reveals an “on-top” adsorbed molecular state prior to the transition state. While a previously reported flat-lying molecular metastable state can be identified in the potential energy surface, it is absent in the FES at relevant reaction temperatures. Finally, these findings demonstrate the importance of identifying critical points self-consistently on the FES for reactions that involve considerable entropic effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2025

Site-specific surface reactivity on MgO for atomic layer deposition via selective hydration

Atomic layer deposition (ALD) is a powerful technique for thin film synthesis, offering atomic-scale precision and conformality. While ALD of MgO has been widely studied for applications in energy storage and microelectronics, its potential as surface on which deposition may be selective and defects repaired remains underexplored. Here, we present a combined theoretical and experimental investigation of MgO surface hydration and its implications for targeted ALD growth using water and dimethyl aluminum isopropoxide (DMAI) as reactants. We perform density functional theory (DFT) calculations to examine molecular and dissociative H 2 O adsorption on MgO (100) terraces and step-edge sites, including pristine surfaces and those with Mg/O vacancies. Reaction Gibbs free energies are calculated under various conditions to quantify surface reactivity. Our findings reveal facet- and defect-dependent hydration behaviors that align with experimental ALD growth trends on MgO (100). This study provides a molecular-level understanding of MgO surface chemistry critical for optimizing ALD processes for thin film growth and defect repair.

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

Quantum Calculations of VX Ammonolysis and Hydrolysis Pathways via Hydrated Lithium Nitride

Recently, lithium nitride (Li 3 N) has been proposed as a chemical warfare agent (CWA) neutralization reagent for its ability to produce nucleophilic ammonia molecules and hydroxide ions in aqueous solution. Quantum chemical calculations can provide insight into the Li 3 N neutralization process that has been studied experimentally. Here, we calculate reaction-free energies associated with the Li3N-based neutralization of the CWA VX using quantum chemical density functional theory and ab initio methods. We find that alkaline hydrolysis is more favorable to either ammonolysis or neutral hydrolysis for initial P-S and P-O bond cleavages. Reaction-free energies of subsequent reactions are calculated to determine the full reaction pathway. Notably, products predicted from favorable reactions have been identified in previous experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗