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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Ammonia Synthesis under Ambient Conditions: Insights into Water–Nitrogen–Magnetite Interfaces

New routes for transforming nitrogen into ammonia at ambient conditions would be a milestone toward an energy efficient and economically attractive production route in comparison to the traditional Haber-Bosch process. Recently, the synthesis of ammonia from water and nitrogen at room temperature and atmospheric pressure has been reported to be catalyzed by Fe 3 O 4 at the air-water interface. By integrating ambient pressure X-ray photoelectron spectroscopy and ab initio molecular dynamics and free energy calculations, we investigate the underlying thermodynamic mechanisms governing ammonia and hydrazine formation at the water-Fe 3 O 4 -nanoparticle interface. Here, we find that, unlike pure Fe 3 O 4 where N 2 can only interact with a limited number of Fe sites, hydroxylated species introduce large and diverse adsorption geometries where N 2 can bind through either Fe sites or Fe-OH groups, each of which are capable of independently facilitating proton-coupled electron transfer.

Chandy, Sruthy K. [University of California, Berke↗
DOE OSTI2021

A comparative study of Bi, Sb, and BiSb for electrochemical nitrogen reduction leading to a new catalyst design strategy

Recent studies identified Bi as one of the most promising non-noble metal elements that can promote the electrochemical N 2 reduction reaction (ENRR) to produce NH 3 . The electronic features that make Bi a promising ENRR catalyst may also be owned by Sb that belongs to the same group as Bi. Thus, the ENRR properties of Bi, Sb, and a BiSb alloy were investigated comparatively to identify common characteristics that facilitate the ENRR. These catalysts were prepared as uniform coating layers on high surface area carbon felt electrodes, which could serve as both regular electrodes and pseudo-gas diffusion electrodes. The experimental results demonstrated that while Bi and Sb show comparable ENRR performances, the formation of a BiSb alloy distinctively increases the faradic efficiency for NH 3 production. Additionally, the X-ray photoelectron spectroscopy results revealed that Bi in BiSb possesses a partial positive charge while Sb in BiSb possesses a partial negative charge, which can impact the way the catalyst surface interacts with the reactants and reaction intermediates of the ENRR and hydrogen evolution reaction (HER), the major competing reaction with the ENRR. Computational investigations including the Bader charge analysis and Gibbs free energy calculations for the elemental steps of the ENRR and HER provided an explanation of how the formation of a BiSb alloy can change the selectivity for the ENRR. The combined experimental and theoretical results and discussion contained in this study lead to a new strategy for designing efficient metal catalysts for the ENRR. Additionally, this study investigated how the use of gas phase and dissolved N 2 affected the ENRR performances of the Bi, Sb, and BiSb catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nanoconfinement matters in humidified CO 2 interaction with metal silicates

With enigmatic observations of enhanced reactivity of wet CO 2 -rich fluids with metal silicates, the mechanistic understanding of molecular processes governing carbonation proves critical in designing secure geological carbon sequestration and economical carbonated concrete technologies. Here, we use the first principle and classical molecular simulations to probe the impact of nanoconfinement on physicochemical processes at the rock–water–CO 2 interface. We choose nanoporous calcium–silicate–hydrate (C–S–H) and forsterite (Mg 2 SiO 4 ) as model metal silicate surfaces that are of significance in cement chemistry and geochemistry communities, respectively. We show that while a nanometer-thick interfacial water film persists at unsaturated conditions consistent with in situ infrared spectroscopy, the phase behavior of the water–CO 2 mixture changes from its bulk counterpart depending on the surface chemistry and nanoconfinement. We also observe enhanced solubility at the interface of water and CO 2 phases, which could amplify the CO 2 speciation rate. Additionally, through free energy calculations, we show that CO 2 could be found in a metastable state near the C–S–H surface, which can potentially react with surface water and hydroxyl groups to form carbonic acid and bicarbonate. These findings support the explicit consideration of nanoconfinement effects in reactive and non-reactive pore-scale processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Confinement-induced clustering of H 2 and CO 2 gas molecules in hydrated nanopores

Gas molecule clustering within nanopores holds significance in the fields of nanofluidics, biology, gas adsorption/desorption, and geological gas storage. However, the intricate roles of nanoconfinement and surface chemistry that govern the formation of gas clusters remain inadequately explored. In this study, through free energy calculation in molecular simulations, we systematically compared the tendencies of H 2 and CO 2 molecules to aggregate within hydrated hydrophobic pyrophyllite and hydrophilic gibbsite nanopores. The results indicate that nanoconfinement enhances gas dimer formation in the nanopores, irrespective of surface chemistry. However, surface hydrophilicity prohibits the formation of gas clusters larger than dimers, while large gas clusters form easily in hydrophobic nanopores. Despite H 2 and CO 2 both being non-polar, the larger quadrupole moment of CO 2 leads to a stronger preference for dimer/cluster formation compared to H 2 . Here, our results also indicate that gases prefer to enter the nanopores as individual molecules, but exit the nanopores as dimers/clusters. This investigation provides a mechanistic understanding of gas cluster formation within nanopores, which is relevant to various applications, including geological gas storage.

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

Monte Carlo Simulations of Crystal Defects in Open Ensembles

Zero- and two-dimensional crystal defects form in open statistical ensembles, such as the grand canonical, that are usually inaccessible with conventional simulation techniques. This longstanding challenge is overcome with a new Hamiltonian Monte Carlo method that samples energy-biased gradual transformations. In conclusion, the method enables free energy calculations for nonideal point defects and the direct prediction of finite-temperature interface structures.

Grain boundaries↗
DOE OSTI2024

De novo design of drug-binding proteins with predictable binding energy and specificity

The de novo design of small molecule–binding proteins has seen exciting recent progress; however, high-affinity binding and tunable specificity typically require laborious screening and optimization after computational design. We developed a computational procedure to design a protein that recognizes a common pharmacophore in a series of poly(ADP-ribose) polymerase–1 inhibitors. One of three designed proteins bound different inhibitors with affinities ranging from <5 nM to low micromolar. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free energy calculations performed directly on the designed models were in excellent agreement with the experimentally measured affinities. We conclude that de novo design of high-affinity small molecule–binding proteins with tuned interaction energies is feasible entirely from computation.

Science & Technology - Other Topics↗
DOE OSTI2024

Zero and Finite Temperature Quantum Simulations Powered by Quantum Magic

We introduce a quantum information theory-inspired method to improve the characterization of many-body Hamiltonians on near-term quantum devices. We design a new class of similarity transformations that, when applied as a preprocessing step, can substantially simplify a Hamiltonian for subsequent analysis on quantum hardware. By design, these transformations can be identified and applied efficiently using purely classical resources. In practice, these transformations allow us to shorten requisite physical circuit-depths, overcoming constraints imposed by imperfect near-term hardware. Importantly, the quality of our transformations is t u n a b l e : we define a 'ladder' of transformations that yields increasingly simple Hamiltonians at the cost of more classical computation. Using quantum chemistry as a benchmark application, we demonstrate that our protocol leads to significant performance improvements for zero and finite temperature free energy calculations on both digital and analog quantum hardware. Specifically, our energy estimates not only outperform traditional Hartree-Fock solutions, but this performance gap also consistently widens as we tune up the quality of our transformations. In short, our quantum information-based approach opens promising new pathways to realizing useful and feasible quantum chemistry algorithms on near-term hardware.

Physics↗
DOE OSTI2022

Thermodynamics Modeling for Actinide Monocarbides and Mononitrides from First Principles

The high-temperature thermodynamical properties for the actinide monocarbides and mononitrides ThC, ThN, UC, UN, PuC, and PuN are calculated from first-principles electronic-structure theory. The electronic structure is modeled with density-functional theory (DFT) and is fully relativistic, including the spin-orbit interaction. Furthermore, the DFT is extended to account for orbital–orbital interactions, by means of a parameter-free orbital-polarization (OP) technique, that has proven to be essential for the 5f electrons in plutonium. Strong anharmonicity and the temperature dependence of the lattice vibrations are captured with the self-consistent ab initio lattice dynamics (SCAILD) method. The calculated free energies and heat capacities are compared to published results from quasi-harmonic (QH) theory, and experiments, where available. For the uranium and plutonium compounds, we make use of CALPHAD assessments to help evaluate the theory. Generally, our anharmonic relativistic approach compares well with both CALPHAD and experiments. For the thorium compounds, our theory is in good accord with QH modeling of the free energy at lower temperatures but for the heat capacity the comparison is less favorable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Enhancing Paraoxon Binding to Organophosphorus Hydrolase Active Site

Organophosphorus hydrolase (OPH) is a metalloenzyme that can hydrolyze organophosphorus agents resulting in products that are generally of reduced toxicity. The best OPH substrate found to date is diethyl p-nitrophenyl phosphate (paraoxon). Most structural and kinetic studies assume that the binding orientation of paraoxon is identical to that of diethyl 4-methylbenzylphosphonate, which is the only substrate analog co-crystallized with OPH. In the current work, we used a combined docking and molecular dynamics (MD) approach to predict the likely binding mode of paraoxon. Then, we used the predicted binding mode to run MD simulations on the wild type (WT) OPH complexed with paraoxon, and OPH mutants complexed with paraoxon. Additionally, we identified three hot-spot residues (D253, H254, and I255) involved in the stability of the OPH active site. We then experimentally assayed single and double mutants involving these residues for paraoxon binding affinity. The binding free energy calculations and the experimental kinetics of the reactions between each OPH mutant and paraoxon show that mutated forms D253E, D253E-H254R, and D253E-I255G exhibit enhanced substrate binding affinity over WT OPH. Interestingly, our experimental results show that the substrate binding affinity of the double mutant D253E-H254R increased by 19-fold compared to WT OPH.

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

Structure and energetics of model amphiphilic molecules at the water liquid-vapor interface - A molecular dynamics study

A molecular dynamics study of adsorption of p-n-pentylphenol at infinite dilution at the water liquid-vapor interface is reported. The calculated free energy of adsorption is -8.8 +/- 0.7 kcal/mol, in good agreement with the experimental value of -7.3 kcal/mol. The transition between the interfacial region and the bulk solution is sharp and well-defined by energetic, conformational, and orientational criteria. At the water surface, the phenol head group is mostly immersed in aqueous solvent. The most frequent orientation of the hydrocarbon tail is parallel to the interface, due to dispersion interactions with the water surface. This arrangement of the phenol ring and the alkyl chain requires that the chain exhibits a kink. As the polar head group is being moved into the solvent, the chain length increases and the tail becomes increasingly aligned toward the surface normal, such that the nonpolar part of the molecule exposed to water is minimized. The same effect was achieved when phenol was replaced by a more polar head group, phenolate.

Pohorille, Andrew↗
NASA NTRS1987

Chemistry Research

Philip Morris research center scientists use a computer program called CECTRP, for Chemical Equilibrium Composition and Transport Properties, to gain insight into the behavior of atoms as they progress along the reaction pathway. Use of the program lets the scientist accurately predict the behavior of a given molecule or group of molecules. Computer generated data must be checked by laboratory experiment, but the use of CECTRP saves the researchers hundreds of hours of laboratory time since experiments must run only to validate the computer's prediction. Philip Morris estimates that had CECTRP not been available, at least two man years would have been required to develop a program to perform similar free energy calculations.

Source record↗
NASA NTRS2010

Interactions of Ku70/80 with Double-Strand DNA: Energetic, Dynamics, and Functional Implications

Space radiation is a proficient inducer of DNA damage leading to mutation, aberrant cell signaling, and cancer formation. Ku is among the first responding proteins in nucleus to recognize and bind the DNA double strand breaks (DSBs) whenever they are introduced. Once loaded Ku works as a scaffold to recruit other repair factors of non-homologous end joining and facilitates the following repair processes. The crystallographic study of the Ku70/80 heterodimer indicate the core structure of this protein shows virtually no conformational change after binding with DNA. To investigate the dynamical features as well as the energetic characteristics of Ku-DNA binding, we conduct multi-nanosecond molecular dynamics simulations of a modeled Ku70/80 structure and several complexes with two 24-bp DNA duplexes. Free energy calculations show significant energy differences between the complexes with Ku bound at DSBs and those with Ku associated at an internal site of a chromosome. The results also reveal detailed interactions between different nucleotides and the amino acids along the DNA-binding cradle of Ku, indicating subtle binding preference of Ku at specific DNA sequences. The covariance matrix analyses along the trajectories demonstrate the protein is stimulated to undergo correlated motions of different domains once bound to DNA ends. Additionally, principle component analyses identify these low frequency collective motions suitable for binding with and translocation along duplex DNA. It is proposed that the modification of dynamical properties of Ku upon binding with DSBs may provide a signal for the further recruitment of other repair factors such as DNA-PKcs, XLF, and XRCC4.

Hu, Shaowen↗
NASA NTRSDate not supplied

Special Grain Boundaries in NiTi Shape Memory Alloys as Sites for Preferential Martensite Nucleation

Shape memory alloys (SMAs) exhibit unique thermomechanical properties due to a reversible martensitic transformation which can be controlled via chemical composition and microstructural features. Among the latter, grain boundaries (GBs) are key in determining how the transformation nucleates and propagates. We performed molecular dynamics simulations to examine the roles of a few special GBs in the austenite phase of NiTi. The GBs can act as preferential host sites for martensite nuclei which can substantially lower the nucleation barrier for transformation and thereby result in reduced thermal hysteresis, an important SMA property for cyclic actuation applications. Free energy calculations show that a characteristic of these GBs is a negative entropy which drops sharply close to the transformation temperature, an anomalous behavior for GBs with fixed composition which is a direct result of the martensitic transformation. We discuss implications of these results with respect to SMA processing techniques for achieving improved properties.

Gabriel Plummer↗
NASA NTRS1984

Theoretical study of miscibility and glass-forming trends in mixtures of polystyrene spheres

A theoretical study of glass-forming trends and miscibility in mixtures of polystyrene spheres (polyballs) of different diameters, suspended in an aqueous solution, is presented. The polyballs are assumed to be charged and to interact via a Debye-Hueckel screened Coulomb potential. The Helmholtz free energy is calculated from a variational principle based on the Gibbs-Bogoliubov inequality, in which a mixture of hard spheres of different diameters is chosen as the reference system. It is found that when the charges of the two types of polyballs are sufficiently different, the variationally determined ratio of hard-sphere diameters differs substantially, leading to packing difficulties characteristic of glass formation. The experimentally observed range of glass formation corresponds to a ratio of hard-sphere diameters of 0.8 or less. Calculations of the free energy as a function of concentration indicate that the liquid polyball mixture is stable against the phase separation, even for widely different polyball charges.

Shih, W.-H.↗
DOE OSTI2021

Using a Coarse-Grained Modeling Framework to Identify Oligomeric Motifs with Tunable Secondary Structure

Coarse-grained modeling can be used to explore general theories that are independent of specific chemical detail. In this paper, we present cg_openmm, a Python-based simulation framework for modeling coarse-grained hetero-oligomers and screening them for structural and thermodynamic characteristics of cooperative secondary structures. cg_openmm facilitates the building of coarse-grained topology and random starting configurations, setup of GPU-accelerated replica exchange molecular dynamics simulations with the OpenMM software package, and features a suite of postprocessing thermodynamic and structural analysis tools. In particular, native contact analysis, heat capacity calculations, and free energy of folding calculations are used to identify and characterize cooperative folding transitions and stable secondary structures. In this work, we demonstrate the capabilities of cg_openmm on a simple 1–1 Lennard-Jones coarse-grained model, in which each residue contains 1 backbone and 1 side-chain bead. By scanning both nonbonded and bonded force-field parameter spaces at the coarse-grained level, we identify and characterize sets of parameters which result in the formation of stable helices through cooperative folding transitions. Furthermore, we show that the geometries and stabilities of these helices can be tuned by manipulating the force-field parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Chemomechanics in alloy phase stability

We describe a first-principles statistical mechanics method to calculate the free energies of crystalline alloys that depend on temperature, composition, and strain. The approach relies on an extension of the alloy cluster expansion to include an explicit dependence on homogeneous strain in addition to site occupation variables that track the degree of chemical ordering. The method is applied to the Si-Ge binary alloy and is used to calculate free energies that describe phase stability under arbitrary epitaxial constraints. We find that while the incoherent phase diagram (in which coexisting phases are not affected by coherency constraints) hosts a miscibility gap, coherent phase equilibrium predicts ordering and negative enthalpies of mixing. Instead of chemical instability, the chemomechanical free energy exhibits instabilities along directions that couple the composition of the alloy with a volumetric strain order parameter. Furthermore, this has fundamental implications for phase field models of spinodal decomposition as it indicates the importance of gradient energy coefficients that couple gradients in composition with gradients in strain.

Materials Science↗