SEARCH · Engineering Papers

Results for “Molecular simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Learning local equivariant representations for large-scale atomistic dynamics

Abstract A simultaneously accurate and computationally efficient parametrization of the potential energy surface of molecules and materials is a long-standing goal in the natural sciences. While atom-centered message passing neural networks (MPNNs) have shown remarkable accuracy, their information propagation has limited the accessible length-scales. Local methods, conversely, scale to large simulations but have suffered from inferior accuracy. This work introduces Allegro, a strictly local equivariant deep neural network interatomic potential architecture that simultaneously exhibits excellent accuracy and scalability. Allegro represents a many-body potential using iterated tensor products of learned equivariant representations without atom-centered message passing. Allegro obtains improvements over state-of-the-art methods on QM9 and revMD17. A single tensor product layer outperforms existing deep MPNNs and transformers on QM9. Furthermore, Allegro displays remarkable generalization to out-of-distribution data. Molecular simulations using Allegro recover structural and kinetic properties of an amorphous electrolyte in excellent agreement with ab-initio simulations. Finally, we demonstrate parallelization with a simulation of 100 million atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Exploring the impacts of conformer selection methods on ion mobility collision cross section predictions

The prediction of structure dependent molecular properties, such as collision cross sections as measured using ion mobility spectrometry, are crucially dependent on the selection of the correct population of molecular conformers. Here, we report an in-depth evaluation of multiple conformation selection techniques, including simple averaging, Boltzmann weighting, lowest energy selection, low energy threshold reductions, and similarity reduction. Generating 50,000 conformers each for 18 molecules, we used the In Silico Chemical Library Engine (ISiCLE) to calculate the collision cross sections for the entire dataset. First, we employed Monte Carlo simulations to understand the variability between conformer structures as generated using simulated annealing. Then we employed Monte Carlo simulations to the aforementioned conformer selection techniques applied on the simulated molecular property—the ion mobility collision cross section. Based on our analyses, we found Boltzmann weighting to be a good tradeoff between precision and theoretical accuracy. Combining multiple techniques revealed that energy thresholds and root-mean-squared deviation-based similarity reductions can save considerable computational expense while maintaining property prediction accuracy. Molecular dynamic conformer generation tools like AMBER can continue to generate new lowest energy conformers even after tens of thousands of generations, decreasing precision between runs. This reduced precision can be ameliorated and theoretical accuracy increased by running density functional theory geometry optimization on carefully selected conformers.

Nielson, Felicity F.↗

DOE OSTI2025

Carbon mineralization pathways in interfacial adsorbed water nanofilms

Carbon mineralization in humidified carbon dioxide offers a promising route to mitigate anthropogenic emissions in a world stressed by water security. Despite its technological importance, our understanding of carbonation in water-poor environments lags, as traditional dissolution-precipitation pathways struggle to explain the adsorbed water nanofilm-mediated reactivity. Here, we utilize in operando X-ray diffraction (XRD) and advanced molecular simulations to investigate nanoconfined reactions driving forsterite carbonation, the magnesium-rich olivine. By examining magnesium ion dissolution and transport in atomistic simulations of the forsterite-water-carbon dioxide interface and comparing these with the in operando XRD activation energies, we identify both processes as rate-limiting at saturation. Our simulations reveal a mechanistic view of interfacial carbonation, where dissolution and precipitation are mediated by anomalous quasi two-dimensional diffusion. The transport process involves intermittent diffusive hopping in the desorbed state, separated by crawling events that are spatially short but temporally long. This understanding transcends carbon mineralization, with implications for understanding the transport of contaminants in geosystems, the design of multifunctional materials, water desalination, and molecular recognition systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Controlling Structure and Properties of Vapor-Deposited Glasses of Organic Semiconductors: Recent Advances and Challenges

The past decade has seen great progress in manipulating the structure of vapor-deposited glasses of organic semiconductors. Upon varying the substrate temperature during deposition, glasses with a wide range of density and molecular orientation can be prepared from a given molecule. Here, we review recent studies that show the structure of vapor-deposited glasses can be tuned to significantly improve the external quantum efficiency and lifetime of organic light-emitting diodes (OLEDs). We highlight the ability of molecular simulations to reproduce experimentally observed structures, setting the stage for in silico design of vapor-deposited glasses in the coming decade. Finally, we identify research opportunities for improving the properties of organic semiconductors by controlling the structure of vapor-deposited glasses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Elucidating molecular level interfacial interactions between a de novo protein and nucleated calcite with solid-state NMR

Biomineralization is the process by which organisms use biomolecules to produce hierarchically structured organic–inorganic composites. Using biology as inspiration, a protein construct (FD31) was previously designed to accelerate formation of nano-calcite with an unconventional {110} face. Here, to understand the molecular interactions essential for protein aided calcite nucleation, solid-state nuclear magnetic resonance (ssNMR) spectroscopy was used in this work to characterize the FD31–calcite interface at the atomic level. Glutamic acid side chains designed to interact directly with calcium ions on the surface were found to have dynamics on the sub-millisecond timescale, indicating possible interactions between the protein and surface waters that were not included in the original model. Dipolar ssNMR recoupling techniques also showed that the protein backbone is ∼2 Å closer to the surface than in the original docking model. Refined molecular simulations were done in the presence of explicit waters, which resulted in the protein backbone closer to the surface than in the original docking structure, providing better agreement with experiment and highlighting the important role played by water in FD31–calcite interactions. While this work provides the first experimental confirmation that FD31 interactions with calcite are localized to the surface of the protein designed to serve as a template, these studies do indicate a more dynamic binding and closer binding mode between FD31 and the nucleated surface than originally proposed. In all, this enhanced molecular insight into the FD31–calcite interface has advanced our fundamental understanding of the atomic interactions at the organic–inorganic interface and will aid in the design of biological templates for the nucleation of inorganic crystals.

Close, Emily G. S. [Pacific Northwest National Lab↗

DOE OSTI2021

Adsorption space for microporous polymers with diverse adsorbate species

Abstract The enormous number of combinations of adsorbing molecules and porous materials that exist is known as adsorption space. The adsorption space for microporous polymers has not yet been systematically explored, especially when compared with efforts for crystalline adsorbents. We report molecular simulation data for the adsorptive and structural properties of polymers of intrinsic microporosity with a diverse set of adsorbate species with 345 distinct adsorption isotherms and over 240,000 fresh and swollen structures. These structures and isotherms were obtained using a sorption-relaxation technique that accounts for the critical role of flexibility of the polymeric adsorbents. This enables us to introduce a set of correlations that can estimate adsorbent swelling and fractional free volume dilation as a function of adsorbate uptake based on readily characterized properties. The separation selectivity of the 276 distinct binary molecular pairs in our data is reported and high-performing adsorbent systems are identified.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Visualizing group II intron dynamics between the first and second steps of splicing

Group II introns are ubiquitous self-splicing ribozymes and retrotransposable elements evolutionarily and chemically related to the eukaryotic spliceosome, with potential applications as gene-editing tools. Recent biochemical and structural data have captured the intron in multiple conformations at different stages of catalysis. Here, we employ enzymatic assays, X-ray crystallography, and molecular simulations to resolve the spatiotemporal location and function of conformational changes occurring between the first and the second step of splicing. We show that the first residue of the highly-conserved catalytic triad is protonated upon 5’-splice-site scission, promoting a reversible structural rearrangement of the active site (toggling). Protonation and active site dynamics induced by the first step of splicing facilitate the progression to the second step. Our insights into the mechanism of group II intron splicing parallels functional data on the spliceosome, thus reinforcing the notion that these evolutionarily-related molecular machines share the same enzymatic strategy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

DataSet for Elucidating molecular level interfacial interactions between a de novo protein and nucleated calcite with solid-state NMR

Biomineralization is the process by which organisms use biomolecules to produce hierarchically structured organic-inorganic composites. Using biology as inspiration, a protein construct (FD31) was previously designed to accelerate formation of nano-calcite with an unconventional {110} face. To understand the molecular interactions essential for protein aided calcite nucleation, solid-state nuclear magnetic resonance (ssNMR) spectroscopy was used in this work to characterize the FD31-calcite interface at the atomic level. Glutamic acid side chains designed to interact directly with calcium ions on the surface were found to have dynamics on the sub-millisecond timescale, indicating possible interactions between the protein and surface waters that were not included in the original model. Dipolar ssNMR recoupling techniques also showed that the protein backbone is ~2 Å closer to the surface than in the original docking model. Refined molecular simulations were done in the presence of explicit waters, which resulted in the protein backbone closer to the surface than in the original docking structure, providing better agreement with experiment and highlighting the important role played by water in FD31-calcite interactions. These studies provide the first experimental evidence to confirm that FD31 interactions with calcite are localized to the surface of the protein designed to serve as a template. However, these studies do indicate a more dynamic binding and closer binding mode between FD31 and the nucleated surface than originally proposed. In all, this enhanced molecular insight into the FD31-calcite interface has advanced our fundamental understanding of the atomic interactions at the organic-inorganic interface and will aid in the design of biological templates for the nucleation of inorganic crystals.

Saccuzzo Close, Emily Grace [Pacific Northwest Nat↗

DOE OSTI2026

Field-driven ion pairing dynamics in concentrated electrolytes

Here, we investigate ion pairing dynamics in electrolytes driven far from equilibrium using molecular simulations and nonequilibrium rate theory. Focusing on 0.5M LiPF 6 in water and acetonitrile under uniform electric fields, we compute transition path theory observables, including reactive fluxes and mean first-passage times of ion pairing. Moreover, we introduce a dynamical proxy of free-ion population, where its field-induced change is strongly correlated with the nonlinear enhancement of conductivity, yielding an increase of 40% at 50 mV/Å in acetonitrile, compared to that of less than 10% in aqueous electrolytes. Further kinetic analysis elucidates that Onsager’s classical theory substantially overestimates field-induced enhancement of ion pair dissociation in molecular electrolytes. This discrepancy arises from solvent-mediated dynamical pathways and field-induced dielectric decrement that suppress ion pair dissociation within explicit solvents, highlighting that a faithful description of molecular details is essential. Our results provide a molecular interpretation of nonlinear electrolyte transport beyond continuum theories and establish a general framework for quantifying nonequilibrium reaction kinetics in condensed phase systems.

Moon, Seokjin [University of California, Berkeley,↗

DOE OSTI2023

Structural and biochemical insight into a modular β-1,4-galactan synthase in plants

Rhamnogalacturonan I (RGI) is a structurally complex pectic polysaccharide with a backbone of alternating rhamnose and galacturonic acid residues substituted with arabinan and galactan side chains. Galactan synthase 1 (GalS1) transfers galactose and arabinose to either extend or cap the ..beta..-1,4-galactan side chains of RGI, respectively. Here we report the structure of GalS1 from Populus trichocarpa, showing a modular protein consisting of an N-terminal domain that represents the founding member of a new family of carbohydrate-binding module, CBM95, and a C-terminal glycosyltransferase family 92 (GT92) catalytic domain that adopts a GT-A fold. GalS1 exists as a dimer in vitro, with stem domains interacting across the chains in a 'handshake' orientation that is essential for maintaining stability and activity. In addition to understanding the enzymatic mechanism of GalS1, here we gained insight into the donor and acceptor substrate binding sites using deep evolutionary analysis, molecular simulations and biochemical studies. Combining all the results, a mechanism for GalS1 catalysis and a new model for pectic galactan side-chain addition are proposed.

09 BIOMASS FUELS↗

DOE OSTI2023

Molecular mechanisms of thickness-dependent water desalination in polyamide reverse-osmosis membranes

Here, using non-equilibrium molecular simulations, we systematically elaborate the relationships among synthesis (membrane's thickness and morphology), atomic-scale transport mechanism, and separation performance (permeability and selectivity) for 3D-printed polyamide (PA) membranes. Results indicated that water diffusion, swelling ratio, water flux, and water permeance proportionally decrease with increasing membrane thickness (4.0–32.5 nm). PA membranes with different thicknesses can achieve almost 100% salt rejection over the simulation time. Importantly, water permeability exponentially decreases with increasing thickness, and 15 nm is identified as the critical membrane thickness for efficient water transport. The discontinuous water-available space spreads all over PA membranes with thicknesses greater than 5 nm, allowing water molecules to jump by way of the temporary open-and-close pores. However, the connected water-useable space exists in PA membranes with thicknesses below 5 nm, offering the continuous channels to dominate water transport. More significantly, pore distribution is more homogeneous as the thickness increases. The applied high pressures can lead to membrane compaction during reverse osmosis and the thicker membranes show a lower compression ratio. In short, these investigations provide molecular insights for effectively designing and manufacturing PA membranes for water desalination and treatment at the molecular level.

42 ENGINEERING↗

DOE OSTI2025

Chromatin structures from integrated AI and polymer physics model

The physical organization of the genome in three-dimensional space regulates many biological processes, including gene expression and cell differentiation. Three-dimensional characterization of genome structure is critical to understanding these biological processes. Direct experimental measurements of genome structure are challenging; computational models of chromatin structure are therefore necessary. We develop an approach that combines a particle-based chromatin polymer model, molecular simulation, and machine learning to efficiently and accurately estimate chromatin structure fromindirectmeasures of genome structure. More specifically, we introduce a new approach where the interaction parameters of the polymer model are extracted from experimental Hi-C data using a graph neural network (GNN). We train the GNN on simulated data from the underlying polymer model, avoiding the need for large quantities of experimental data. The resulting approach accurately estimates chromatin structures across all chromosomes and across several experimental cell lines despite being trained almost exclusively on simulated data. The proposed approach can be viewed as a general framework for combining physical modeling with machine learning, and it could be extended to integrate additional biological data modalities. Ultimately, we achieve accurate and high-throughput estimations of chromatin structure from Hi-C data, which will be necessary as experimental methodologies, such as single-cell Hi-C, improve.

Biochemistry & Molecular Biology↗

DOE OSTI2021

Potential Foldamers Based on an ortho- Terphenyl Amino Acid

We describe the synthesis and characterization of a new class of oligomers built from a terphenyl-based amino acid. Here, these oligomeric amides are of interest because the adoption of specific conformations could potentially be driven by the coordinated formation of inter-residue hydrogen bonds and aromatic interactions. Although high-resolution structural data have proven inaccessible, circular dichroism and nuclear magnetic resonance studies suggest that the new oligomers fold concomitantly with discrete self-association in chloroform. Molecular simulations find no sign of folding as a monomer, and in fact suggest that hydrogen bond patterning is anti-cooperative, inhibiting longer helices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Multifunctional Catalysts for the Tandem Reactions of Oxygenates

Industrially-relevant catalytic reactions rarely consist of a simple sequence of elementary steps. Moreover, kinetic coupling of multiple reactions on a catalyst surface is highly desired for process intensification and improved energy efficiency for large scale chemical transformations. The shifting landscape of hydrocarbon chemical feedstocks in the US also motivates research on the selective conversion of more complex molecules. One desirable type of catalytic reaction is the reduction of carboxylic acids that are produced from biomass feedstocks to their corresponding alcohols. The proposed research explores the fundamental importance of hydrogen spillover on a multifunctional catalyst for carboxylic acid reduction with H 2 composed of metal particles coupled to metal oxide particles. Recent work has demonstrated the excellent performance of supported tungsten oxide clusters for carboxylic acid reduction, but only after they are promoted with a late transition metal such as palladium. Elucidating the active state of the catalyst and the associated reaction mechanism for acid reduction on that active state are the overall goals of the proposed project and successful completion will enable future design of efficient multifunctional catalysts. The critically important role of the metal promoter is hypothesized to be its ability to dissociate H2 and spillover atomic H to the support. Although hydrogen spillover is a well-recognized phenomenon in catalysis, its role in both catalyst activation and catalytic turnover are still unresolved. The study combined materials synthesis, characterization, reactivity testing, and molecular simulations, to explore the effect of hydrogen chemical potential on the formation of the active catalytic sites and on the steady state catalytic reduction of carboxylic acid. Varying the hydrogen chemical potential through modification of the gas conditions, support composition, and metal loading to modulated the structure and catalytic performance of the reducible metal oxide. Dual function catalysts containing supported Pd and WO x species co-located on a non-reducible carrier (silica) and a reducible carrier (titania) were synthesized and characterized by electron microscopy, temperature-programmed reduction, and chemisorption. Spectroscopic methods such as X-ray absorption and UV-vis were also used to evaluate the catalysts, which were used in the reduction of propionic acid to aldehyde and alcohol. Quantum chemical calculations, including ab initio phase diagrams provided molecular insights into the H spillover phenomenon.

09 BIOMASS FUELS↗

DOE OSTI2021

Understanding Solid Electrolyte Interphase Nucleation and Growth on Lithium Metal Surfaces

Experiments and theory are needed to decode the exact structure and distribution of components of a passivation layer formed at the anode surface of Li metal batteries, known as the Solid Electrolyte Interphase (SEI). Due to the inherent dynamic behavior as well as the lithium reactivity, the SEI structure and its growth mechanisms are still unclear. This study uses molecular simulation and computational chemistry tools to investigate the initial nucleation and growth dynamics of LiOH and Li2O that provide us with thermodynamics and structural information about the nucleating clusters of each species. Following the most favorable pathways for the addition of each of the components to a given nascent SEI cluster reveals their preferential nucleation mechanisms and illustrates different degrees of crystallinity and electron density distribution that are useful to understand ionic transport through SEI blocks.

25 ENERGY STORAGE↗

DOE OSTI2022

Production of β-ketoadipic acid from glucose in Pseudomonas putida KT2440 for use in performance-advantaged nylons

Biomass-derived chemicals can offer unique chemical functionality relative to petroleum-derived building blocks. To this end, here we report that β-ketoadipic acid (βKA), a C6 diacid with a β-ketone group, can be used as a performance-advantaged replacement for adipic acid in a nylon-6,6 analog. Building on our previous efforts to produce shikimate-derived products from carbohydrates, Pseudomonas putida KT2440 is engineered to produce βKA from glucose, achieving a 26 g/L titer. Following purification, βKA imparts an increase of 69 degrees C above the nylon-6,6 glass transition temperature and 20% reduced water permeability, equivalent to nylon-6,10. Molecular simulations predict that the enhanced thermal properties result from rigidity introduced by the β-ketone. Process analysis predicts that βKA can be produced for US$1.94/kg from sugars, requiring 63% less energy and emitting 43% less greenhouse gases than fossil-based adipic acid. Overall, this study illustrates the potential for βKA to serve as a useful building block for bio-based polymers.

09 BIOMASS FUELS↗

DOE OSTI2022

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗