Engineering Papers⌕ Search

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.

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

At least 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Sequence-defined Pareto frontier of a copolymer structure

The correlations between the sequence of monomers in a macromolecule and its three-dimensional (3D) structure is a grand challenge in polymer science. The properties and functions of macromolecules depend on their 3D shape that has appeared to be dictated by their monomer sequence. However, the progress towards understanding the sequence–structure-property correlations and their utilization in materials engineering are slow because it is almost impossible to characterize an astronomically large number of possible sequences of a copolymer using traditional experimental and simulation methods. To address this problem, here, we combine evolutionary computing and coarse-grained molecular dynamics (CGMD) simulation and study the sequence-structure correlations of a model AB-type copolymer in a solution and assess the impact of sequence on the packing density in its bulk phase. The CGMD-based evolutionary algorithm (EA) screens the sequence space of a single chain copolymer efficiently and identifies a wide range of single-molecule structures including extremal radii of gyration. The data are utilized to estimate the Pareto front of the structure-space of a binary copolymer as a function of its composition. The monomer packings in single-molecule solution phase and multimolecular bulk phase are found to be identical. Finally, this work highlights the opportunities of sequence-specific control of macromolecular structure for designing target materials.

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

Isotropic-to-Nematic Transition in Salt-Free Polyelectrolyte Coacervates from Coarse-Grained Simulations

Recent interest in complex coacervation between oppositely charged polyelectrolytes (PEs) has been fueled by its relevance to biology in the context of membraneless organelle formation within living cells. For PEs with limited flexibility (such as double-stranded DNA), theoretical treatments and recent experiments have reported the emergence of liquid crystalline order (LCO) within the resulting coacervate phases. In this work, we study the underlying physics of this phenomenon using coarse-grained molecular dynamics simulations of symmetric semiflexible-semiflexible and asymmetric semiflexible-flexible coacervates. By comparing coacervates with the corresponding semidilute solutions of neutral polymers, we demonstrate that the presence of Coulomb interactions in coacervates facilitates orientational ordering, in agreement with theoretical predictions. Quantitative comparisons between our simulations and theory indicate that, for asymmetric nematic coacervates, the strong orientational ordering of stiff polyanions induces a weak ordering of the flexible polycations-an effect that was not anticipated by available theoretical studies. Simulations reveal that, for nematic coacervates, the preferred orientation of the PE chains at the liquid-liquid coacervate-supernatant interface is parallel, and the alignment of semiflexible PEs is homogeneous. The results presented here provide new molecular-level insights into the intracoacervate LCO and will help motivate further experimental and theoretical activities in this area.

36 MATERIALS SCIENCE↗
DOE OSTI2021

EZFF: Python library for multi-objective parameterization and uncertainty quantification of interatomic forcefields for molecular dynamics

Parameterization of interatomic forcefields is a necessary first step in performing molecular dynamics simulations. This is a non-trivial global optimization problem involving quantification of multiple empirical variables against one or more properties. We present EZFF, a lightweight Python library for parameterization of several types of interatomic forcefields implemented in several molecular dynamics engines against multiple objectives using genetic-algorithm-based global optimization methods. The EZFF scheme provides unique functionality such as the parameterization of hybrid forcefields composed of multiple forcefield interactions as well as built-in quantification of uncertainty in forcefield parameters and can be easily extended to other forcefield functional forms as well as MD engines.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Probing the state of hydrogen in δ − AlOOH at mantle conditions with machine learning potential

Hydrous and nominally anhydrous minerals are a fundamental class of solids of enormous significance to geophysics. They are the water carriers in the deep geological water cycle and impact structural, elastic, plastic, and thermodynamic properties and phase relations in Earth's forming aggregates (rocks). They play a critical role in the geochemical and geophysical processes that shape the planet. Their complexity has prevented predictive calculations of their properties, but progress in materials simulations ushered by machine-learning potentials is transforming this state of affairs. Here, we adopt a hybrid approach that combines deep learning potentials (DPs) with the strongly constrained and appropriately normed meta-generalized gradient approximation functional to simulate a prototypical hydrous system. We illustrate the success of this approach to simulate δ − AlOOH ( δ ) , a phase capable of transporting water down to near the core-mantle boundary of the Earth ( ∼ 2900 km depth and ∼ 135 GPa ) in subducting slabs. A high-throughput sampling of phase space using molecular dynamics simulations with DPs sheds light on the hydrogen-bond behavior and proton diffusion at geophysical conditions. These simulations provide a pathway for a deeper understanding of these crucial components that shape Earth's internal state. Published by the American Physical Society 2024

Luo, Chenxing (ORCID:0000000341166851)↗
DOE OSTI2026

From 2D to 4D: a containerized workflow and browser to explore dynamic chromatin architecture

Background Characterizing the physical organization of the genome is essential for understanding long-range gene regulation, chromatin compartmentalization, and epigenetic accessibility. Hi-C experiments generate two-dimensional (2D) genome-wide contact maps of chromatin interactions by capturing the spatial proximity between genomic loci, which reveal interaction frequencies but lack the spatial resolution needed to interpret the three-dimensional (3D) genome structure(s). Emerging evidence suggests that epigenetic regulation is closely linked to 3D genome architecture, and that structural changes over time (4D) drive key biological processes in development, disease, and environmental response. Thus, integrating 3D structure with functional data is critical for a more complete understanding of genome regulation. Previous work, most notably the 4DHiC chromosome modeling framework, has shown that physical multi-dimensional modeling approaches rooted in polymer physics and molecular dynamics can resolve these structures at biologically meaningful resolutions by integrating temporal Hi-C data with physical constraints to uncover dynamic chromosome reorganization. Thus, molecular dynamics simulations, constrained by Hi-C contact matrices, can resolve fine-scale structural changes and reveal functionally significant transitions in chromatin conformation. Results Herein, we present the 4D Genome Browser Workflow (4DGBWorkflow) and the 4D Genome Browser (4DGB). The algorithm is based on the 4DHiC method, and the containerized tool is an end-to-end workflow that can transform, filter, and view 4D epigenomics and chromatin datasets, allowing non-specialists to apply three-dimensional modeling principles to diverse datasets and experimental conditions. The software executes on a laptop running macOS, Linux or Windows. From input Hi-C files (.hic), the 4DGBWorkflow produces 3D reconstructions of chromosomes, integrates the reconstruction with track data (e.g., epigenetic marks, transcriptome profiles), and provides comparative visualization of the results in a single workflow. Conclusions The 4DGBWorkflow and 4D Genome Browser are open-source tools for comparative analysis and visualization of 4D chromosome datasets, including chromatin architecture and epigenomic signals. Automatic integration of Hi-C data with molecular dynamics democratizes the construction of time resolved 3D genome structures, simplifying complex simulations and data integration schemes.

3D Genome Browser↗
DOE OSTI2020

Liquid Structure and Transport Properties of the Deep Eutectic Solvent Ethaline

A range of techniques including physical property measurements, neutron scattering experiments, ab initio molecular dynamics, and classical molecular dynamics simulations are used to probe the structural, thermodynamic, and transport properties of a deep eutectic solvent comprised of a 1:2 molar ratio of choline chloride and ethylene glycol. This mixture, known as Ethaline, has many desirable properties for use in a range of applications, and therefore, understanding its liquid structure and transport properties is of interest. Simulation results are able to capture experimental densities, diffusivities, viscosities, and structure factors extremely well. Here, the solvation environment is dynamic and dominated by different hydrogen bonding interactions. Dynamic heterogeneities resulting from hydrogen bonding interactions are quantified. Rotational dynamics of molecular dipole moments of choline and ethylene glycol are computed and found to exhibit a fast and slow mode.

25 ENERGY STORAGE↗
DOE OSTI2022

Probing electrolyte–silica interactions through simulations of the infrared spectroscopy of nanoscale pores

The structural and dynamical properties of nanoconfined solutions can differ dramatically from those of the corresponding bulk systems. Understanding the changes induced by confinement is central to controlling the behavior of synthetic nanostructured materials and predicting the characteristics of biological and geochemical systems. A key outstanding issue is how the molecular-level behavior of nanoconfined electrolyte solutions is reflected in different experimental, particularly spectroscopic, measurements. This is addressed here through molecular dynamics simulations of the OH stretching infrared (IR) spectroscopy of NaCl, NaBr, and NaI solutions in isotopically dilute HOD/D 2 O confined in hydroxylated amorphous silica slit pores of width 1–6 nm and pH ~2. In addition, the water reorientation dynamics and spectral diffusion, accessible by pump–probe anisotropy and two-dimensional IR measurements, are investigated. In this work, the aim is to elucidate the effect of salt identity, confinement, and salt concentration on the vibrational spectra. It is found that the IR spectra of the electrolyte solutions are only modestly blue-shifted upon confinement in amorphous silica slit pores, with both the size of the shift and linewidth increasing with the halide size, but these effects are suppressed as the salt concentration is increased. This indicates the limitations of linear IR spectroscopy as a probe of confined water. However, the OH reorientational and spectral diffusion dynamics are significantly slowed by confinement even at the lowest concentrations. The retardation of the dynamics eases with increasing salt concentration and pore width, but it exhibits a more complex behavior as a function of halide.

Mesoporous material↗
DOE OSTI2020

Resolving Heterogeneous Dynamics of Excess Protons in Aqueous Solution with Rate Theory

Rate theories have found great utility across the chemical sciences by providing a physically transparent way to analyze dynamical processes. Here we demonstrate the benefits of using transition state theory and Marcus theory to study the rate of proton transfer in HCl solutions. By using long ab initio molecular dynamics simulations we show that good agreement is obtained between these two different formulations of rate theory and how they can be used to study the pathways and life-time of proton transfer in aqueous solution. Since both rate theory formulations utilize identical sets of molecular data, this provides a self-consistent theoretical picture of the rates of aqueous phase proton transfer. Specifically, we isolate and quantify the rates of proton transfer, ion-pair dissociation, and solvent exchange in concentrated HCl solutions. Our analysis predicts a concentration dependence to both proton transfer and ion-pairing. Moreover, our estimate of the life-time for the Zundel species is 0.8 ps and 1.3 ps for 2M and 8M HCl, respectively. In this paper, we demonstrate that concentration effects, can indeed be quantified through the combination of state-of-the-art simulation and theory and provides a picture of the important correlations between the cation (hydronium) and the counter-ion in acid solutions.

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

Structural basis of agonist specificity of α 1A -adrenergic receptor

α 1 -adrenergic receptors (α 1 -ARs) play critical roles in the cardiovascular and nervous systems where they regulate blood pressure, cognition, and metabolism. However, the lack of specific agonists for all α 1 subtypes has limited our understanding of the physiological roles of different α 1 -AR subtypes, and led to the stagnancy in agonist-based drug development for these receptors. Here we report cryo-EM structures of α 1A -AR in complex with heterotrimeric G-proteins and either the endogenous common agonist epinephrine or the α 1A -AR-specific synthetic agonist A61603. These structures provide molecular insights into the mechanisms underlying the discrimination between α 1A -AR and α 1B -AR by A61603. Guided by the structures and corresponding molecular dynamics simulations, we engineer α 1A -AR mutants that are not responsive to A61603, and α 1B -AR mutants that can be potently activated by A61603. Together, these findings advance our understanding of the agonist specificity for α 1 -ARs at the molecular level, opening the possibility of rational design of subtype-specific agonists.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Thermodynamic driving forces in contact electrification between polymeric materials

Abstract Contact electrification, or contact charging, refers to the process of static charge accumulation after rubbing, or even simple touching, of two materials. Despite its relevance in static electricity, various natural phenomena, and numerous technologies, contact charging remains poorly understood. For insulating materials, even the species of charge carrier may be unknown, and the direction of charge-transfer lacks firm molecular-level explanation. Here, we use all-atom molecular dynamics simulations to investigate whether thermodynamics can explain contact charging between insulating polymers. Based on prior work suggesting that water-ions, such as hydronium and hydroxide ions, are potential charge carriers, we predict preferred directions of charge-transfer between polymer surfaces according to the free energy of water-ions within water droplets on such surfaces. Broad agreement between our predictions and experimental triboelectric series indicate that thermodynamically driven ion-transfer likely influences contact charging of polymers. Furthermore, simulation analyses reveal how specific interactions of water and water-ions proximate to the polymer-water interface explain observed trends. This study establishes relevance of thermodynamic driving forces in contact charging of insulators with new evidence informed by molecular-level interactions. These insights have direct implications for future mechanistic studies and applications of contact charging involving polymeric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atoms to fibers: Identifying novel processing methods in the synthesis of pitch-based carbon fibers

Understanding and optimizing the key mechanisms used in the synthesis of pitch-based carbon fibers (CFs) are challenging, because unlike polyacrylonitrile-based CFs, the feedstock for pitch-based CFs is chemically hetero- geneous, resulting in complex fabrication leading to inconsistency in the final properties. In this work, we use molecular dynamics simulations to explore the processing and chemical phase space through a framework of CF models to identify their effects on elastic performance. The results are in excellent agreement with experiments. We find that density, followed by alignment, and functionality of the molecular constituents dictate the CF mechanical properties more strongly than their size and shape. Last, we propose a previously unexplored fabrication route for high-modulus CFs. Unlike graphitization, this results in increased sp 3 fraction, achieved via generating high-density CFs. In addition, the high sp 3 fraction leads to the fabrication of CFs with isometric compressive and tensile moduli, enabling their potential applications for compressive loading.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Advancing Multiscale Simulation of Plasma-Surface Interfaces

We report the development of an atomistic-informed, surface-state-dependent predictive model for particle exchange in a carbon-tungsten plasma-surface interface. The predictive model uses machine learning (ML) techniques to learn the energy and angular distributions for particle exchange and rate functions for surface state evolution from molecular dynamics simulations of cumulative bombardment of tungsten by energetic carbon ions. Each predictive component is sensitive to the energy and trajectory of incident plasma species and the surface state. The surface state is represented by a set of surface state descriptors, which were derived from the atomistic surface state for each independent carbon bombardment event. These descriptors are representative of the composition and degree of amorphization of the outermost angstrom of surface material and were chosen to optimize predictive performance for particle exchange at the interface. The distributions for particle exchange (reflection/sputtering) are demonstrated to vary with each surface state descriptor, motivating the development of surface-state-dependent particle exchange models for plasma simulations. The performance of various ML methods was compared, including polynomial quantile regression, artificial neural networks, k-nearest neighbors, and random forest algorithms, with polynomial regression performing the best for interpolation and extrapolation of learned relationships. In addition to the particle exchange model, a neutral network was developed and used to identify data sufficiency throughout surface descriptor space, which will enable real-time feedback during future data production to ensure data is produced where it is most needed, and we provide commentary on improvements to the data production workflow for future endeavors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Experimental and Computational Insights into the Structural Dynamics of the Fc Fragment of IgG1 Subtype from Biosimilar VEGF‐Trap

The constant fragment (Fc) of the immunoglobulin G1 (IgG1) subtype is increasingly recognized as a crucial scaffold in the development of advanced therapeutics due to its enhanced specificity, efficacy, and extended half‐life. A prime example is VEGF‐Trap (Aflibercept), a recombinant fusion protein that merges the Fc region of the IgG1 subtype with the binding domains of vascular endothelial growth factor receptors (VEGFR)‐1 and VEGFR‐2. The Fc region's role in N‐glycosylation is particularly important, as it significantly influences protein stability. Herein, the first near‐physiological temperature structures of the N‐glycan‐bound Fc fragment of IgG1 subtype from a biosimilar VEGF‐Trap are presented, determined using the SPring‐8 Angstrom Compact free electron LAser (SACLA) and the Turkish Light Source (Turkish DeLight). Comparative analysis with cryogenic structures, including existing data, reveals alternate conformations within the glycan‐binding pocket. Furthermore, molecular dynamics simulations indicate the presence of a high degree of structural plasticity, explaining how the protein adapts its structure through conformational changes. The observed structural fluctuations/conformational changes demonstrate the effect of N‐glycans on protein stability. These findings offer new insights into the molecular basis of Fc‐mediated functions and provide valuable information for the design of next‐generation therapeutics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Mechanistic Insights into Molecular Copper Hydride Catalysis: the Kinetic Stability of CuH Monomers toward Aggregation is a Critical Parameter for Catalyst Performance

The activity of molecular copper hydride (CuH) complexes towards the selective insertion of unsaturated hydrocarbons under mild conditions has contributed significantly to versatile methodologies for upgrading these feedstocks. However, these catalysts are particularly susceptible to deleterious aggregation, leading to the depletion of active CuH species. Little is known about the mechanisms of CuH aggregation, how it influences overall catalyst performance, and how it can be controlled. We address these challenges with mechanistic studies on a model reaction of unactivated alkene hydroboration catalyzed by (IPr*CPh 3 )CuH (LCuH). Here, we report a comprehensive mechanistic investigation of this system, identifying an aggregation pathway that continuously depletes catalytically active LCuH to form inactive CuH clusters during turnover. Deactivation of LCuH is controlled primarily by the competition between the kinetics of the initial LCuH dimerization step and that of alkene insertion. We therefore propose that a more comprehensive understanding of CuH catalyst performance must account for the kinetics of the initial LCuH dimerization step, revising a previously explored thermodynamic understanding of CuH aggregation, where the concentration of active species is controlled by equilibria established between CuH dimers and monomers. With a series of (NHC)CuH congeners (NHC = N-heterocyclic carbene), we demonstrate that ostensibly minor structural modifications to the ligand peripheries can drastically affect the LCuH dimerization kinetics, while maintaining reactivity towards on–cycle alkene insertion. We employed a computational approach based on molecular dynamics simulations to provide an in-depth understanding of how specific structural ligand modifications can substantially increase the kinetic stability of monomeric CuH catalysts. Our combined experimental and computational studies suggest strategies for rational ligand design that can be broadly applied to molecular catalyst systems that are susceptible to deactivation via aggregation pathways.

Ryan, David E. [Pacific Northwest National Laborat↗
NASA NTRS2017

Virtual Design of a 4-Bed Molecular Sieve for Exploration

Simulations of six new 4-Bed Molecular Sieve configurations have been performed using a COMSOL (COMSOL Multiphysics - commercial software) model. The preliminary results show that reductions in desiccant bed size and sorbent bed size when compared to the International Space Station configuration are feasible while still yielding a process that handles at least 4.0 kilograms a day CO2. The results also show that changes to the CO2 sorbent are likewise feasible. Decreasing the bed sizes was found to have very little negative effect on the adsorption process; breakthrough of CO2 in the sorbent bed was observed for two of the configurations, but a small degree of CO2 breakthrough is acceptable, and water breakthrough in the desiccant beds was not observed. Both configurations for which CO2 breakthrough was observed still yield relatively high CO2 efficiency, and future investigations will focus on bed size in order to find the optimum configuration.

Giesy, Timothy J.↗
NASA NTRS2017

Virtual Design of a 4-Bed Molecular Sieve for Exploration

Simulations of six new 4-Bed Molecular Sieve configurations have been performed using a COMSOL model. The preliminary results show that reductions in desiccant bed size and sorbent bed size when compared to the International Space Station configuration are feasible while still yielding a process that handles at least 4.0 kg/day CO2. The results also show that changes to the CO2 sorbent are likewise feasible. Decreasing the bed sizes was found to have very little negative effect on the adsorption process; breakthrough of CO2 in the sorbent bed was observed for two of the configurations, but water breakthrough in the desiccant beds was not observed. Nevertheless, both configurations for which CO2 breakthrough was observed still yield relatively high CO2 efficiency, and future investigations will focus on bed size in order to find the optimum configuration.

Giesy, Timothy J.↗
DOE OSTI2020

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. In this study, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗