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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2026

Diffusion Quantum Monte Carlo Calculation of the Austenite and Martensite Phases of NiTi

NiTi is a promising material for smart and active technologies due to its exhibition of the shape memory effect, superelasticity, and biocompatibility. The shape memory effect is tied to the reversible transition between the austenite and martensite phases. A major research direction is to alloy NiTi with Zr, Hf, Pd, Pt, etc., in order to tune the martensitic transition temperature (MTT). Modeling the MTT from first principles is challenging because the lattice dynamics is complicated by anharmonicity and various low-energy structures. Using density functional theory, the energy difference between the austenite and martensite phases of NiTi varies by up to 100 meV/atom depending on the choice of density functional, which is of the same order of the energy difference itself. Consequently, free energy calculations with different functionals can result in estimates of the MTT that vary by several hundred K. Using diffusion quantum Monte Carlo, we calculated the energy difference between the B2 and B19' structures of NiTi to be 70.9 +- 2.5 meV/atom.

Kevin K Ly↗
DOE OSTI2023

Impact of lanthanide ion complexation and temperature on the chemical reactivity of N , N , N ', N '-tetraoctyl diglycolamide (TODGA) with the dodecane radical cation

The impact of trivalent lanthanide ion complexation and temperature on the chemical reactivity of N,N,N',N'-tetraoctyl diglycolamide (TODGA) with the n-dodecane radical cation (RH˙+) has been measured by electron pulse radiolysis and evaluated by quantum mechanical calculations. Additionally, Arrhenius parameters were determined for the reaction of the non-complexed TODGA ligand with the RH˙ + from 10–40 °C, giving the activation energy (E a = 17.43 ± 1.64 kJ mol –1 ) and pre-exponential factor (A = (2.36 ± 0.05) × 10 13 M –1 s –1 ). The complexation of Nd(III), Gd(III), and Yb(III) ions by TODGA yielded [LnIII(TODGA)3(NO3)3] complexes that exhibited significantly increased reactivity (up to 9.3× faster) with the RH˙ + , relative to the non-complexed ligand: k([Ln III (TODGA) 3 (NO 3 ) 3 ] + RH˙ + ) = (8.99 ± 0.93) × 10 10 , (2.88 ± 0.40) × 10 10 , and (1.53 ± 0.34) × 10 10 M –1 s –1 , for Nd(III), Gd(III), and Yb(III) ions, respectively. The rate coefficient enhancement measured for these complexes exhibited a dependence on atomic number, decreasing as the lanthanide series was traversed. Preliminary reaction free energy calculations—based on a model [Ln III (TOGDA)] 3+ complex system—indicate that both electron/hole and proton transfer reactions are energetically unfavorable for complexed TODGA. Furthermore, complementary average local ionization energy calculations showed that the most reactive region of model N,N,N',N'-tetraethyl diglycolamide (TEDGA) complexes, [Ln III (TEGDA) 3 (NO 3 ) 3 ], toward electrophilic attack is for the coordinated nitrate (NO 3 – ) counter anions. Furthermore, it is possible that radical reactions with the complexed NO 3 – counter anions dominate the differences in rates seen for the [Ln III (TODGA) 3 (NO 3 ) 3 ] complexes, and are likely responsible for the reported radioprotection in the presence of TODGA complexes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Dual Enhancement of Thermostability and Activity of Xylanase through Computer-Aided Rational Design

In the realm of enzyme engineering, the dual enhancement of thermostability and activity remains a challenge. Herein, we employed a computer-aided approach integrating folding free energy calculations and evolutionary analysis to engineer Paecilomyces thermophila xylanase into a hyperthermophilic enzyme for application in the paper and pulp industry. Through the computational rational design, XynM9 with superior thermostability and enhanced activity was designed. Its optimal reaction temperature increases by 10 °C to 85 °C, its T m increases by 10 °C to 93 °C, and its half-life increases 11-fold to 5.8 h. Additionally, its catalytic efficiency improves by 57% to 3926 s –1 mM –1 . Molecular dynamics simulations revealed that XynM9 is stabilized by more hydrogen bonds and salt bridges than wild-type xylanase. The mutant’s narrower catalytic cleft enhances the substrate-binding affinity, thus improving the catalytic efficiency. In harsh conditions at 80 °C and pH 10, using XynM9 significantly reduced both hemicellulose and lignin, which makes it a good candidate for use in the paper and pulp process. Finally, our study presents an accurate and efficient strategy for the dual enhancement of enzyme properties, guiding further improvement of computational tools for protein stabilization.

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

The Role of Cu 3+ in the Oxygen Evolution Activity of Copper Oxides

Cu-based oxides and hydroxides represent an important class of materials from a catalytic and corrosion perspective. In this study, we investigate the formation of bulk and surface Cu 3+ species that are stable under water oxidation catalysis in alkaline media. So far, no direct evidence existed for the presence of hydroxides (CuOOH) or oxides, which were primarily proposed by theory. Here, this work directly places CuOOH in the oxygen evolution reaction (OER) Pourbaix stability region with a calculated free energy of −208.68 kJ/mol, necessitating a revision of known Cu–H 2 O phase diagrams. We also predict that the active sites of CuOOH for the OER are consistent with a bridge O* site between the two Cu 3+ atoms with onset at ≥1.6 V vs the reversible hydrogen electrode (RHE), aligning with experimentally observed Cu 2+/3+ oxidation waves in cyclic voltammetry of Fe-free and Fe-spiked copper in alkaline media. Trace amounts of Fe (2 μg/mL (ppm) to 5 μg/mL) in the solution measurably enhance the catalytic activity of the OER, likely due to the adsorption of Fe species that serve as the active sites . Importantly, modulation excitation X-ray absorption spectroscopy (ME-XAS) of a Cu thin-film electrode shows a distinct Cu 3+ fingerprint under OER conditions at 1.8 V vs RHE. Additionally, in situ Raman spectroscopy of polycrystalline Cu in 0.1 mol/L (M) KOH revealed features consistent with those calculated for CuOOH in addition to CuO. Overall, this work provides direct evidence of bulk electrochemical Cu 3+ species under OER conditions and expands our longstanding understanding of the oxidation mechanism and catalytic activity of copper.

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

Superoxide Dismutase 1 Folding Stability as a Target for Molecular Tweezers in SOD1-Related Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation are hallmarks of many severe neurodegenerative diseases including Alzheimer's, Parkinson's and Huntington's disease. As a supramolecular ligand that binds to lysine and arginine residues, the molecular tweezer CLR01 was found to modify the aggregation pathway of disease-relevant proteins in vitro and in vivo with beneficial effects on toxicity. However, the molecular mechanisms of how tweezers exert these effects remain mainly unknown, hampering further drug development. Here, we investigate the modulation mechanism of unfolding and aggregation pathways of SOD1, which are involved in amyotrophic lateral sclerosis (ALS), by CLR01. Using a truncated version of the wildtype SOD1 protein, SOD1 bar , we show that CLR01 acts on the first step of the aggregation pathway, the unfolding of the SOD1 monomer. CLR01 increases, by ~10 °C, the melting temperatures of the A4V and G41D SOD1 mutants, which are commonly observed mutations in familial ALS. Molecular dynamics simulations and binding free energy calculations as well as native mass spectrometry and mutational studies allowed us to identify K61 and K92 as binding sites for the tweezers to mediate the stability increase. The data suggest that the modulation of SOD1 conformational stability is a promising target for future developments of supramolecular ligands against neurodegenerative diseases.

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

PTM‐Psi : A python package to facilitate the computational investigation of p ost‐ t ranslational m odification on p rotein s tructures and their i mpacts on dynamics and functions

Abstract Post‐translational modification (PTM) of a protein occurs after it has been synthesized from its genetic template, and involves chemical modifications of the protein's specific amino acid residues. Despite of the central role played by PTM in regulating molecular interactions, particularly those driven by reversible redox reactions, it remains challenging to interpret PTMs in terms of protein dynamics and function because there are numerous combinatorially enormous means for modifying amino acids in response to changes in the protein environment. In this study, we provide a workflow that allows users to interpret how perturbations caused by PTMs affect a protein's properties, dynamics, and interactions with its binding partners based on inferred or experimentally determined protein structure. This Python‐based workflow, called PTM‐Psi , integrates several established open‐source software packages, thereby enabling the user to infer protein structure from sequence, develop force fields for non‐standard amino acids using quantum mechanics, calculate free energy perturbations through molecular dynamics simulations, and score the bound complexes via docking algorithms. Using the S ‐nitrosylation of several cysteines on the GAP2 protein as an example, we demonstrated the utility of PTM‐Psi for interpreting sequence–structure–function relationships derived from thiol redox proteomics data. We demonstrate that the S ‐nitrosylated cysteine that is exposed to the solvent indirectly affects the catalytic reaction of another buried cysteine over a distance in GAP2 protein through the movement of the two ligands. Our workflow tracks the PTMs on residues that are responsive to changes in the redox environment and lays the foundation for the automation of molecular and systems biology modeling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Lpp positions peptidoglycan at the AcrA-TolC interface in the AcrAB-TolC multidrug efflux pump

The multidrug efflux pumps of Gram-negative bacteria are a class of complexes that span the periplasm, coupling both the inner and outer membranes to expel toxic molecules. The best-characterized example of these tripartite pumps is the AcrAB-TolC complex of Escherichia coli. However, how the complex interacts with the peptidoglycan (PG) cell wall, which is anchored to the outer membrane (OM) by Braun’s lipoprotein (Lpp), is still largely unknown. In this work, we present molecular dynamics simulations of a complete, atomistic model of the AcrAB-TolC complex with the inner membrane, OM, and PG layers all present. We find that the PG localizes to the junction of AcrA and TolC, in agreement with recent cryo-tomography data. Free-energy calculations reveal that the positioning of PG is determined by the length and conformation of multiple Lpp copies anchoring it to the OM. The distance between the PG and OM measured in cryo-electron microscopy images of wild-type E. coli also agrees with the simulation-derived spacing. Finally, sequence analysis of AcrA suggests a conserved role for interactions with PG in the assembly and stabilization of efflux pumps, one that may extend to other trans-envelope complexes as well.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗
DOE OSTI2021

High-salinity brine desalination with amine-based temperature swing solvent extraction: A molecular dynamics study

Solvent extraction desalination (SED) is a promising technique for efficiently treating high salinity brines. However, only a few studies in the literature have investigated the influence of the molecular characteristics of the organic solvents on the emergent performance of the separation process. Thus, in this work we investigate the influence of the molecular structure of four different secondary amines on the desalination process of a concentrated aqueous solution of NaCl. Molecular dynamics simulations and free energy calculations are used to obtain the thermodynamic and structural properties of the pure solvent bulk phase, as well as a homogeneous solvent-water bulk mixture. Furthermore, solvent-water and solvent-brine interfacial systems are used to explicitly analyze the water partitioning behavior in the relevant temperature range for SED. The water solubilization in these secondary amines is found to be an exothermic process. Additionally, although the hydrogen-bonding interaction between the water and solvent molecules are beneficial, this interaction is not necessarily the dominant factor for ranking the water solubility performance of the different solvents. Finally, the excess volume of mixing of the bulk solvent-water mixture, as well as the electrostatic surface features of the solvents are also found to provide important clues for screening solvent behavior.

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

Graphene-supported single atom catalysts for high performance lithium-oxygen batteries

The optimal choice of d-block metals in single atom catalysts (SACs) is crucial for designing efficient electrocatalysts for activating the Oxygen reduction reaction (ORR)/ Oxygen evolution reaction (OER) in lithiumoxygen batteries (LOBs). Herein, we used the Quantum Mechanics methods to understand the origin of reactivity for a series of 16 d-block metals supported on nitrogen-doped graphene as SACs for ORR and OER in LOBs. Based on the Gibbs free energy calculations, we found that among the 16 SACs investigated, Zn-SAC exhibits the highest electrochemical activity with the lowest overpotential of 0.17 V. Here we then used machine learning (ML) to develop an intrinsic descriptor, phi, that correlates the catalytic activity with electronic and chemical properties of the catalytic centers at the M-N 4 active site on graphene surface. We established a linear relationship between phi and the catalytic activity that provides guidance for designing efficient SACs for electrocatalysis in LOBs. To validate these predictions, we report electrochemical measurements showing that Zn-SAC exhibits an ultra-stable cyclability with reduced overpotentials over Mo-SAC and nitrogen-doped graphene (NG), confirming our theoretical prediction. This fundamental work provides a deep understanding on the rational design of efficient SACs for OER/ ORR in LOBs.

25 ENERGY STORAGE↗
DOE OSTI2025

Design, Synthesis, and Biological Activity of Novel Ornithine Decarboxylase (ODC) Inhibitors

We here describe the design, synthesis, and biological activity of novel ornithine decarboxylase (ODC) inhibitors that show significantly higher potency in vitro than α-difluoromethylornithine (DFMO), a U.S. Food and Drug Administration (FDA) approved drug. We report two X-ray structures of ODC complexed with new ODC inhibitors, computational docking, molecular dynamics, and binding free energy calculations to validate the experimental models. The X-ray structures reveal that covalent adducts with pyridoxal phosphate (PLP) are formed in the active site of the human ODC enzyme, as verified by their preparation and enzymatic testing. Finally, we verified that the cellular activity of endogenous ODC was inhibited, and polyamine levels were reduced. Given that ODC is a clinically validated target, combined with the fact that DFMO is currently the only ODC inhibitor in clinical use for several indications, the further development of more potent ODC inhibitors with superior activity and physical properties is warranted.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Structure, Stability, and Electronic Properties of Boron Suboxide: A Density Functional Theory Study

Boron suboxide (B 6 O) is a boron-rich compound derived from the α-rhombohedral boron lattice with extreme hardness and unusual semiconducting properties. In this work, density functional theory (DFT) was used to show that unit cell volume, mechanical strength, band gaps, and thermodynamic stabilities of B 6 O were influenced by the interstitial elements and point defects at the icosahedral sites. While the hexagonal unit cell volume (HUCV) varies with interstitial occupancy, it is the icosahedral defect that weakens the intrinsic bulk modulus of B 6 O. Using the hybrid HSE functional, we confirmed that the perfect B 6 O bulk is a p-type semiconductor with a direct band gap of 2.8 eV. Furthermore, by screening α-boron compounds systematically, we found that a simple octet rule may offer a consistent explanation for the variations in the computed electronic structures. The formation free energies calculated over a wide range of temperatures (0–2500 K) and pressures (0–80 GPa) predict that formations of interstitial defects become favorable only at higher temperatures (ca. 1800 K) in bulk B 6 O lattices. The nudged elastic band (NEB) method was employed to identify the minimum energy pathways for the diffusions of dislocated B and O atoms. The diffusion of icosahedral B atoms has an energy barrier of 0.16 eV. More complex B diffusion paths involving the reorganization of icosahedral boron atoms incur higher barriers (>1 eV). In contrast, the diffusion of interstitial O atoms is facile with a barrier of 0.4 eV. Lastly, successive O insertions into the α-B lattice were performed using DFT to generate a basic understanding of the oxidation process. Furthermore, these calculations provide fundamental atomistic insights into the growth of B 6 O crystals and control of their point defects.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Direct Simulations of H–He Mixtures at Planetary Interior Conditions: Demixing, Insulator–Metal Transition and Miscibility Boundaries

Accurate knowledge of the electrical and thermal conductivities and structural properties of hydrogen–helium mixtures under thermodynamic conditions within and beyond the immiscibility range is very important to predict the thermal evolution and internal structure of gas giant planets like Jupiter and Saturn. Here, we propose a novel method to determine the immiscibility boundary accurately without the need for free energy calculations, while providing consistent insights into structural and transport properties of mixtures. We show with direct large-scale ab initio simulations that the insulator–metal transition (IMT) of the hydrogen subsystem is strongly affected by an admixture with a small fraction of helium and occurs at temperatures significantly higher than those of pure hydrogen. At pressures below 150 GPa, the IMT boundary is not related anymore to the H 2 subsystem dissociation, the system remains insulating even after the full dissociation of H 2 molecules and its transition to an H–He mixture. The offset of the IMT in the H–He mixture relative to the dissociation region in the hydrogen subsystem and the significant reduction of static electrical and thermal conductivity by a factor between two and a few thousand relative to pure hydrogen found in mixtures have consequences for Jupiter and Saturn’s thermal evolution, internal structure, and dynamo action, affecting a large fraction of the interior of both planets.

Helium↗
DOE OSTI2023

Relationships between Water’s Structure and Solute Affinity at Polypeptoid Brush Surfaces

Excellent antifouling surfaces are generally thought to create a tightly bound layer of water that resists solute adsorption, and highly hydrophilic surfaces such as those with zwitterionic functionalities are of significant current interest as antifoulant strategies. Furthermore, despite significant proofs-of-concept, we still lack a fundamental understanding of how the nanoscopic structure of this hydration layer translates to reduced fouling, how surface chemistry can be tuned to achieve antifouling through hydration water, and why, in particular, zwitterionic surfaces seem so promising. Here, we use molecular dynamics simulations and free energy calculations to investigate the molecular relationships among surface chemistry, hydration water structure, and surface–solute affinity across a variety of surface-decorated chemistries. Specifically, we consider polypeptoid-decorated surfaces that display well-known experimental antifouling capabilities and that can be synthesized sequence specifically, with precise backbone positioning of, e.g., charged groups. Through simulations, we calculate the affinities of a range of small solutes to polypeptoid brush surfaces of varied side-chain chemistries. We then demonstrate that measures of the structure of surface hydration water in response to a particular surface chemistry signal solute–surface affinity; specifically, we find that zwitterionic chemistries produce solute–surface repulsion through highly coordinated hydration water while suppressing tetrahedral structuring around the solute, in contrast to uncharged surfaces that show solute–surface affinity. Based on the relationship of this structural perturbation to the affinity of small-molecule solutes, we propose a molecular mechanism by which zwitterionic surface chemistries enhance solute repulsion, with broader implications for the design of antifouling surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular Mechanisms of Sorbed Ion Effects during Boehmite Particle Aggregation

Classical theories of particle aggregation, such as Derjaguin–Landau–Verwey–Overbeek (DLVO), do not explain recent observations of ion-specific effects or the complex concentration dependence for aggregation. Thus, here, we probe the molecular mechanisms by which selected alkali nitrate ions (Na + , K + , and NO 3 –) influence aggregation of the mineral boehmite (γ-AlOOH) nanoparticles. Nanoparticle aggregation was analyzed using classical molecular dynamics (CMD) simulations coupled with the metadynamics rare event approach for stoichiometric surface terminations of two boehmite crystal faces. Calculated free energy landscapes reveal how electrolyte ions alter aggregation on different crystal faces relative to pure water. Consistent with experimental observations, we find that adding an electrolyte significantly reduces the energy barrier for particle aggregation (~3–4×). However, in this work, we show this is due to the ions disrupting interstitial water networks, and that aggregation between stoichiometric (010) basal–basal surfaces is more favorable than between (001) edge–edge surfaces (~5–6×) due to the higher interfacial water densities on edge surfaces. The interfacial distances in the interlayer between aggregated particles with electrolytes (~5–10 Å) are larger than those in pure water (a few Ångströms). Together, aggregation/disaggregation in salt solutions is predicted to be more reversible due to these lower energy barriers, but there is uncertainty on the magnitudes of the energies that lead to aggregation at the molecular scale. By analyzing the peak water densities of the first monolayer of interstitial water as a proxy for solvent ordering, we find that the extent of solvent ordering likely determines the structures of aggregated states as well as the energy barriers to move between them. Further, the results suggest a path for developing a molecular-level basis to predict the synergies between ions and crystal faces that facilitate aggregation under given solution conditions. Such fundamental understanding could be applied extensively to the aggregation and precipitation utilization in the biological, pharmaceutical, materials design, environmental remediation, and geological regimes.

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

Effect of Liquid Electrolyte Soaking on the Interfacial Resistance of Li 7 La 3 Zr 2 O 1 2 for All-Solid-State Lithium Batteries

In this work, the impact of liquid electrolyte soaking on the interfacial resistance between the garnet-structured Li 7 La 3 Zr 2 O 12 (LLZO) solid electrolyte and metallic lithium has been studied. Lithium carbonate (Li 2 CO 3 ) formed by inadvertent exposure of LLZO to ambient conditions is generally known to increase interfacial impedance and decrease lithium wettability. Soaking LLZO powders and pellets in the electrolyte containing lithium tetrafluoroborate (LiBF 4 ) shows a significantly reduced interfacial resistance and improved contact between lithium and LLZO. Raman spectroscopy, X-ray diffraction, and soft X-ray absorption spectroscopy reveal how Li 2 CO 3 is continuously removed with increasing soaking time. On-line mass spectrometry and free energy calculations show how LiBF 4 reacts with surface carbonate to form carbon dioxide. Using a very simple and scalable process that does not involve heat-treatment and expensive coating techniques, we show that the Li–LLZO interfacial resistance can be reduced by an order of magnitude.

25 ENERGY STORAGE↗
DOE OSTI2020

Molecular Dynamics Simulation of Atomic Interactions in the Vancomycin Binding Site

Vancomycin is a glycopeptide antibiotic produced by Amycolaptopsis orientalis used to treat serious infections by Grampositive pathogens including methicillin-resistant Staphylococcus aureus. Vancomycin inhibits cell wall biosynthesis by targeting lipid II, which is the membrane-bound peptidoglycan precursor. The heptapeptide aglycon structure of vancomycin binds to the D-Ala-D-Ala of the pentapeptide stem structure in lipid II. The third residue of vancomycin aglycon is asparagine, which is not directly involved in the dipeptide binding. Nonetheless, asparagine plays a crucial role in substrate recognition, as the vancomycin analogue with asparagine substituted by aspartic acid (V D ) shows a reduction in antibacterial activities. To characterize the function of asparagine, binding of vancomycin and its aspartic-acid-substituted analogue V D to L-Lys-D-Ala-D-Ala and L-Lys-D-Ala-D-Lac was investigated using molecular dynamic simulations. Binding interactions were analyzed using root-mean-square deviation (RMSD), two-dimensional (2D) contour plots, hydrogen bond analysis, and free energy calculations of the complexes. The analysis shows that the aspartate substitution introduced a negative charge to the binding cleft of V D , which altered the aglycon conformation that minimized the repulsive lone pair interaction in the binding of a depsipeptide. Our findings provide new insight for the development of novel glycopeptide antibiotics against the emerging vancomycin-resistant pathogens by chemical modification at the third residue in vancomycin to improve its binding affinity to the D-Ala-D-Lac-terminated peptidoglycan in lipid II found in vancomycin-resistant enterococci and vancomycinresistant S. aureus.

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

Understanding Reaction Networks through Controlled Approach to Equilibrium Experiments Using Transient Methods

We report a combined experimental/theoretical approach to studying heterogeneous gas/solid catalytic processes using low-pressure pulse response experiments achieving a controlled approach to equilibrium that combined with quantum mechanics (QM)-based computational analysis provides information needed to reconstruct the role of the different surface reaction steps. We demonstrate this approach using model catalysts for ammonia synthesis/decomposition. Polycrystalline iron and cobalt are studied via low-pressure TAP (temporal analysis of products) pulse response, with the results interpreted through reaction free energies calculated using QM on Fe-BCC(110), Fe-BCC(111), and Co-FCC(111) facets. In TAP experiments, simultaneous pulsing of ammonia and deuterium creates a condition where the participation of reactants and products can be distinguished in both forward and reverse reaction steps. This establishes a balance between competitive reactions for D* surface species that is used to observe the influence of steps leading to nitrogen formation as the nitrogen product remains far from equilibrium. Here, the approach to equilibrium is further controlled by introducing delay timing between NH 3 and D 2 which allows time for surface reactions to evolve before being driven in the reverse direction from the gas phase. The resulting isotopic product distributions for NH 2 D, NHD 2 , and HD at different temperatures and delay times and NH 3 /D 2 pulsing order reveal the role of the N 2 formation barrier in controlling the surface concentration of NH x * species, as well as providing information on the surface lifetimes of key reaction intermediates. Conclusions derived for monometallic materials are used to interpret experimental results on a more complex and active CoFe bimetallic catalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗