Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic dynamics”

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

Fouling behavior of zwitterionic membranes compared to polyamide membranes

Membrane fouling remains a critical bottleneck for reverse osmosis (RO) desalination, driving energy consumption and reducing membrane lifetime. Here, we employ all-atom molecular dynamics simulations to investigate the antifouling behavior of random zwitterionic amphiphilic copolymer (r-ZAC) membranes composed of sulfobetaine methacrylate (SBMA) and allyl methacrylate (AMA), benchmarked against conventional polyamide (PA) RO membranes. Structural and dynamical analyses—including radial distribution functions, coordination numbers, tetrahedral order parameters, vector orientation, and residence-time correlation functions—reveal that r-ZAC surfaces sustain tightly bound, long-lived hydration layers with preserved tetrahedrality and anisotropic water orientation, in sharp contrast to the weak and disordered hydration of PA. Steered molecular dynamics simulations demonstrate that r-ZAC membranes impose substantial free-energy barriers to foulant approach (alginate ≈ 90 kcal/mol, sucrose ≈ 35 kcal/mol, humic acid ≈ 15 kcal/mol), whereas PA membranes exhibit negligible barriers (< 1 kcal/mol) and thermodynamically favorable adsorption. Detailed foulant–surface interaction analyses show that zwitterionic hydration and electrostatic heterogeneity in r-ZAC suppress adhesion, except in the case of amphiphilic humic acid, which exploits multiple binding modes. Together, these results establish molecular-level design principles for antifouling membranes: the combination of zwitterionic hydration, structured interfacial water, and controlled amphiphilic balance in r-ZAC membranes provides superior resistance to organic fouling relative to PA.

Cross-linked polyamide↗
DOE OSTI2025

Effects of high-pressure hydrogen exposure on filler-elastomer adhesion

Elastomers are known to gain enhanced mechanical properties through compounding with nanosized filler particles such as silica or carbon black. Filler dispersion and filler-polymer interfacial strength are key contributing factors to this improvement. The interfacial strength is critical to part lifetime in pressurized gas sealing applications such as O-rings, where weak binding between the filler particle and polymer matrix can lead to internal void structures. With the aim to build a fundamental understanding of precursors to pressurized hydrogen-induced failure in elastomers, we use all-atom molecular dynamics simulations to study the impact of hydrogen oversaturation on filler-polymer interaction strength. We systematically study the interface between a commonly used elastomer, ethylene-propylene-diene monomer (EPDM) and silica by varying gas concentration, crosslink density, and surface chemistry. Our simulations predict that decompression leads to a localization of excess gas near the interface. In conclusion, we demonstrate that this localized gas can weaken interfacial adhesion and quantify the interaction using thermodynamic approaches.

EPDM↗
DOE OSTI2024

A coarse-grained model of clay colloidal aggregation and consolidation with explicit representation of the electrical double layer

The aggregation of clay minerals in liquid water exemplifies colloidal self-assembly in nature. These negatively charged aluminosilicate platelets interact through multiple mechanisms with different sensitivities to particle shape, surface charge, aqueous chemistry, and interparticle distance and exhibit complex aggregation structures. Experiments have difficulty resolving the associated colloidal assemblages at the scale of individual particles. Conversely, all-atom molecular dynamics (MD) simulations provide detailed insight on clay colloidal interaction mechanisms, but they are limited to systems containing a few particles. We develop a new coarse-grained (CG) model capable of representing assemblages of hundreds of clay particles with accuracy approaching that of MD simulations, at a fraction of the computational cost. Our CG model is parameterized based on MD simulations of a pair of smectite clay particles in liquid water. A distinctive feature of our model is that it explicitly represents the electrical double layer (EDL), i.e., the cloud of charge-compensating cations that surrounds the clay particles. Our model captures the simultaneous importance of long-range colloidal interactions (i.e., interactions consistent with simplified analytical models, already included in extant clay CG models) and short-range interactions such as ion correlation and surface and ion hydration effects. The resulting simulations correctly predict, at low solid-water ratios, the existence of ordered arrangements of parallel particles separated by water films with a thickness up to ~10 nm and, at high solid-water ratios, the coexistence of crystalline and osmotic swelling states, in agreement with experimental observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Shape memory effect in metallic glasses

Shape memory effect (SME), mainly present in crystalline Ti-Ni alloys, is basically missing in metallic glasses (MGs) that lack the long-range periodic order of crystals. In this work, we report experimental results of SME in annealed MGs, in which the low-energy configuration state recovery is observed by both differential scanning calorimetry and X-ray photon correlation spectroscopy. We elucidate the origin of SME in MGs under the potential energy landscape framework, i.e., after annealing, the energy of MGs enters into a deep basin and atoms are located in the low-energy configuration state. Albeit deviating from their relative stable configurations by temperature changes, atoms in the annealed MGs tend to return to the low-energy atomic configurations along certain trajectories as pre-annealing temperature is approached, These results could extend the application of MGs as functional materials by directionally manipulating their energy states via annealing and rejuvenation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Room temperature sodium diffusion in sodium iron phosphate investigated with neutron scattering

Sodium iron phosphate, Na x FePO 4 , is currently one of the most promising cathode materials for sodium ion batteries. Here we employ a systematic approach using three neutron scattering techniques to investigate temperature dependent properties with a focus on the charge transfer processes. Two distinct Na + diffusive motions are observed with quasielastic neutron scattering, including room temperature motion. These processes correspond to rapid double jumps through extrinsic vacancy sites of ∼6 Å along with the more fundamental nearest neighbor jumps of 3.1 Å. Experimental observations of room temperature ionic diffusion are highly desirable, as the diffusion pathways are naturally well suited for efficient transport. Inelastic neutron scattering was used to probe the vibrational density of states that assist in the ionic and electronic (polaronic) transport processes, which are further supported with density functional theory. Neutron diffraction examines the disordering of the sodium sublattice upon transition to the solid-solution phase, which is believed to enhance the overall kinetics during sodium insertion/extraction. Finally, 57 Fe Mössbauer spectroscopy measurements probes polaron hopping and the local structural evolution that modifies the electron density surrounding the redox active iron sites. This study sheds light on the temperature induced atomic scale dynamics that occur in olivine Na x FePO 4 .

Novak, Eric [Swarthmore College, Swarthmore, PA (U↗
DOE OSTI2024

Partial Solvation of Lithium Ions Enhances Conductivity in a Nanophase-Separated Polymer Electrolyte

In this study, we demonstrate that a multiblock lithium-ion-conducting polymer can be swollen with ethylene carbonate solvent to increase the conductivity relative to the dry polymer material by nearly 4 orders of magnitude. This increase is due to the partial solvation of lithium ions by ethylene carbonate, which leads to Li + diffusion along the solvent–polymer interface. This differs from the vehicular transport mechanism for lithium ions in pure solvent. We use a combination of broadband dielectric spectroscopy, X-ray scattering, and all-atom molecular dynamics simulations to probe the effect of the solvent on the polymer morphology and to elucidate the mechanism of lithium ion transport.

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

Solid 3D Li–S Battery Design via Stacking 2D Conductive Microporous Coordination Polymers and Amorphous Li–S Layers

To make a lithium-sulfur (Li-S) battery practical, not only high gravimetric energy capacity is important, but also high volumetric energy capacity will be required. The currently explored Li-S cathode designs often deploy systems with liquid electrolyte infiltration, hence with relatively low volumetric capacity. In the current study, we theoretically test a compact solid three-dimensional (3D) design (more like a Li-ion battery cathode than a conventional Li-S cathode) consisted of a sandwich structure alternating between the two-dimensional (2D) Mn-hexaaminobenzene-based coordination polymer (2D Mn-HAB-CP) layer and the amorphous Li-S layer. We study the theoretical limits for both its gravimetric and volumetric energy capacity, as well as its structural stability and Li diffusion within the cathode system. To study the Li diffusion within an amorphous system, we also develop a pull-atom molecular dynamics (PA-MD) to calculate the barrier heights of such disordered systems. We reveal the mechanism that determines the Li diffusion in the amorphous layer of the system. Overall, we find such a 3D solid Li-S cathode can be practical, with sufficient large gravimetric and volumetric energy capacity, as well as the Li diffusion constant. It also solves many other common Li-S cathode problems, from Li polysulfide dissolution to electrical insulating, and structure instabilities.

25 ENERGY STORAGE↗
DOE OSTI2020

Persistent Protein Motions in a Rugged Energy Landscape Revealed by Normal Mode Ensemble Analysis

Proteins are allosteric machines that couple motions at distinct, often distant, sites to control biological function. Low-frequency structural vibrations are a mechanism of this long-distance connection and are often used computationally to predict correlations, but experimentally identifying the vibrations associated with specific motions has proved challenging. Spectroscopy is an ideal tool to explore these excitations, but measurements have been largely unable to identify important frequency bands. The result is at odds with some previous calculations and raises the question what methods could successfully characterize protein structural vibrations. Here we show the lack of spectral structure arises in part from the variations in protein structure as the protein samples the energy landscape. However, by averaging over the energy landscape as sampled using an aggregate 18.5 μs of all-atom molecular dynamics simulation of hen egg white lysozyme and normal-mode analyses, we find vibrations with large overlap with functional displacements are surprisingly concentrated in narrow frequency bands. These bands are not apparent in either the ensemble averaged vibrational density of states or isotropic absorption. However, in the case of the ensemble averaged anisotropic absorption, there is persistent spectral structure and overlap between this structure and the functional displacement frequency bands. We systematically lay out heuristics for calculating the spectra robustly, including the need for statistical sampling of the protein and inclusion of adequate water in the spectral calculation. The results show the congested spectrum of these complex molecules obscures important frequency bands associated with function and reveal a method to overcome this congestion by combining structurally sensitive spectroscopy with robust normal mode ensemble analysis.

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

Paired Simulations and Experimental Investigations into the Calcium-Dependent Conformation of Albumin

Serum albumin is the most abundant protein in blood plasma, and it is involved in multiple biological processes. Serum albumin has recently been adapted for improving biomaterial integration with bone tissue, and studies have shown the importance of this protein in bone repair and regeneration. However, the mechanism of action is not yet clear. In stark contrast, other studies have demonstrated that albumin blocks cell adhesion to surfaces, which is seen as a limitation to its bone healing role. These apparent contradictions suggest that the conformation of albumin facilitates its bioactivity, leading to enhanced bone repair. Serum albumin is known to play a major role in maintaining the calcium ion concentration in blood plasma. Due to the prevalence of calcium at bone repair and regeneration sites, it has been hypothesized that calcium binding to serum albumin triggers a conformational change, leading to bioactivity. In the current study, molecular modeling approaches including molecular docking, atomic molecular dynamics (MD) simulation, and coarse-grained MD simulation were used to test this hypothesis by investigating the conformational changes induced in bovine serum albumin by interaction with calcium ions. The computational results were qualitatively validated with experimental Fourier-transform infrared spectroscopy analysis. We find that free calcium ions in solution transiently bind with the three major loops in albumin, triggering a conformational change where N-terminal and C-terminal domains separate from each other in a partial unfolding process. The separation distance between these domains was found to correlate with the calcium ion concentration. The experimental data support the simulation results showing that albumin has enhanced conformational heterogeneity upon exposure to intermediate levels of calcium, without any significant secondary structure changes.

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

All-Atom Simulations Uncover Structural and Dynamical Properties of STING Proteins in the Membrane System

Recent studies have shown that the stimulator of interferon gene (STING) protein plays a central role in the immune system by facilitating the production of type I interferons in cells. The STING signaling pathway is also a prominent activator of cancer-killing T cells that initiate a powerful adaptive immune response. Since biomolecular signaling pathways are complicated and not easily identified through traditional experiments, molecular dynamics (MD) has often been used to study structural and dynamical responses of biological pathways. In this work, we carried out MD simulations for full-length chicken and human STING (chSTING and hSTING) proteins. Specifically, we investigated ligand-bound closed (holo) and ligand-unbound open (apo) forms of STING in the membrane system by comparing their conformational and dynamical differences. Our research provides clues for understanding the mechanism of the STING signaling pathway by uncovering detailed insights for the examined systems: the residues from each chain in the binding pocket are strongly correlated to one another in the open STING structure compared with those in the closed STING structure. Ligand-bound closed STING displays ~174° rotation of the ligand-binding domain (LBD) relative to the open STING structure. The dynamical analysis of residue Cys148 located in the linker region of hSTING does not support the earlier hypothesis that Cys148 can form disulfide bonds between adjacent STING dimers. We also demonstrate that using the full-length proteins is critical, since the MD simulations of the LBD portion alone cannot properly describe the global conformational properties of STING.

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

Machine Learning a Simple Interpretable Short-Range Potential for Silica

A wide array of models, spanning from computationally expensive ab initio methods to a spectrum of force-field approaches, have been developed and employed to probe silica polymorphs and understand growth processes and atomic-level dynamical transitions in silica. However, the quest for a model capable of making accurate predictions with high computational efficiency for various silica polymorphs is still ongoing. Recent developments in short-range machine-learned models, such as GAP and NNPScan, have shown promise in providing reasonable descriptions of silica, but their computational cost remains high compared to force fields such as BKS which are based on simple interpretable functional forms. Here, in this study, we build on the recent success of our reinforcement learning (RL) workflow to derive a new set of optimal parameters for a promising short-range BKS-based model proposed by Soules. We use RL to navigate the eight-dimensional parameter space of the Soules potential using an experimental training data set that includes both local and global structural features from approximately 21 experimentally realized silica polymorphs, including high density phases and porous zeolites. We compare the performance of our machine-learned ML-Soules model with other high quality models including our recent machine-learned parametrization of BKS (ML-BKS), a machine-learned potential (GAP), as well as predictions of ab initio calculations with the highly fidelity SCAN functional. The ML-Soules accurately captures the relative energetic ordering of various polymorphs as well as their structural features at a significantly reduced computational expense. The ML-Soules model also reasonably captures the structure, density, and elastic constants of quartz, as well as metastable silica polymorphs. We further discuss the limitations of the Soules functional form and propose potential enhancements, including the incorporation of additional three-body terms and/or the utilization of different short-ranged functional forms to achieve greater accuracy for both global and local features in the modeling of silica while retaining low computational cost.

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

Toward Automated Sampling of Polymorph Nucleation and Free Energies with the SGOOP and Metadynamics

Understanding the driving forces behind the nucleation of different polymorphs is of great importance for material sciences and the pharmaceutical industry. This includes understanding the reaction coordinate that governs the nucleation process and correctly calculating the relative free energies of different polymorphs. Here, we demonstrate, for the prototypical case of urea nucleation from the melt, how one can learn such a one-dimensional reaction coordinate as a function of prespecified order parameters and use it to perform efficient biased all-atom molecular dynamics simulations. The reaction coordinate is learnt as a function of the generic thermodynamic and structural order parameters using the “spectral gap optimization of order parameters (SGOOP)” approach [Tiwary, P. and Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. (2016)] and is biased using well-tempered metadynamics simulations. The reaction coordinate gives insights into the role played by different structural and thermodynamics order parameters, and the biased simulations obtain accurate relative free energies for different polymorphs. This includes an accurate prediction of the approximate pressure at which urea undergoes a phase transition and one of the metastable polymorphs becomes the most stable conformation. Here, we believe the ideas demonstrated in this work will facilitate efficient sampling of nucleation in complex, generic systems.

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

Accurate Calculation of Solvation Properties of Lithium Ions in Nonaqueous Solutions

Here, we perform all-atom molecular dynamics simulations of lithium triflate in 1,2-dimethoxyethane using six different literature force fields. This system is representative of many experimental studies of lithium salts in solvents and polymers. We show that multiple historically common force fields for lithium ions give qualitatively incorrect results when compared with those from experiments and quantum chemistry calculations. We illustrate the importance of correctly selecting force field parameters and give recommendations on the force field choice for lithium electrolyte applications.

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

Is the Local Ion Density Sufficient to Drive NaCl Nucleation from the Melt and Aqueous Solution?

Even though nucleation is ubiquitous in different science and engineering problems, investigating nucleation is extremely difficult due to the complicated ranges of time and length scales involved. In this work, we simulate NaCl nucleation in both molten and aqueous environments using enhanced sampling of all-atom molecular dynamics with deep-learning-based estimation of reaction coordinates. By incorporating various structural order parameters and learning the reaction coordinate as a function thereof, we achieve significantly improved sampling relative to traditional ad hoc descriptions of what drives nucleation, particularly in an aqueous medium. Our results reveal a one-step nucleation mechanism in both environments, with reaction coordinate analysis highlighting the importance of local ion density in distinguishing solid and liquid states. However, although fluctuations in the local ion density are necessary to drive nucleation, they are not sufficient. Our analysis shows that near the transition states, descriptors such as enthalpy and local structure become crucial. Furthermore, our protocol proposed here enables robust nucleation analysis and phase sampling and could offer insights into nucleation mechanisms for generic small molecules in different environments.

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

Side-Chain and Ring-Size Effects on Permeability in Artificial Water Channels

Artificial water channels (AWCs) have emerged as a promising framework for stable water permeation, with water transport rates comparable to aquaporins (3.4–40.3 × 10 8 H 2 O/channel/s). In this study, we probe the influence of ring-size and side-chain length on the water permeability observed within a class of AWCs termed ligand-appended pillar[n]arenes (LAPs) that have an adjustable ring-size (m) and side-chain length (n). Through all-atom molecular dynamics simulations, we calculate the permeability of these channels using the collective diffusion model and find their permeabilities. We characterize the mechanistic influence of pillar[n]arene ring-size and side-chain length on the channel water permeability by analyzing the characteristics of the internal permeating water-wire and the surrounding channel structure. We observe that water permeability decreases as a function of increasing ring-size due to increases in hydrophilic contacts between the permeating water-wire and the oxygen groups on the channel wall. Further, we observe an increase in water permeability as a function of side-chain length due to increased partitioning of the channel terminal groups into the hydrophilic blocks of the surrounding bilayer. For the LAP6 channel, with increase in side-chain length, the distance between terminal groups increases and leads to an increase in pore size, thereby enhancing water permeability. In the case of LAP5, as side-chain length increases, the channel displays a compensatory effect between tilt and bend angle due to the flexible side-chains. Such flexibility leads to higher terminal group partitioning in the hydrophilic blocks of the bilayer and extends the permeating water-wire. Furthermore, this increase in water-wire length and hydrophilic block access overcomes the nonmonotonic pore size trend in pillar[5]arene channels.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Self-Assembled Oligomers Facilitate Amino Acid-Driven CO 2 Capture at the Air–Aqueous Interface

Direct air capture of CO 2 using amino acid absorbents, such as glycine or sarcosine, is constrained by the relatively slow mass transfer of CO 2 through the air–aqueous interface. Our recent study showed a marked improvement in CO 2 capture by introducing CO 2 -permeable oligo-dimethylsiloxane (ODMS-MIM + ) oligomers with cationic (imidazolium, MIM + ) headgroups. Here, in this work, we have employed all-atom molecular dynamics simulations in combination with subensemble analysis using network theory to provide a detailed molecular picture of the behavior of CO 2 and the glycinate anions (Gly – ) at the ODMS-MIM + decorated air–aqueous interfaces. We show that the cationic head groups of the surfactants enhance the concentration and lifetime of Gly – in the interfacial region, while ODMS tails promote the physisorption of CO 2 in the interfacial region. Together, these two factors increase the effective region of contact and the probability of interactions between CO 2 and Gly – compared to that of the pure air–aqueous interface. The fundamental insights gained in this work establish essential foundations for developing hybrid systems with oligomer-decorated interfaces to maximize the overall CO 2 capture rates.

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

Single-Molecule Electron Transport in Peptoids

Peptoids are structural analogs of peptides in which side chains are appended to the backbone nitrogen rather than the α-carbon. The sequence-defined modularity of peptoids enables precise control over structure−function relationships, enabling applications in energy storage and biomedical materials. Despite recent progress, the role of sequence and conformation on electron transport in peptoid molecules is not fully understood. Here, we synthesize a library of peptoid oligomers and characterize their molecular electronic properties using the scanning tunneling microscope-break junction (STM-BJ) technique. Our results show well-defined electron transport behavior for peptoid sequences containing aromatic side groups lacking hydrogen bonds (H-bonds) and without chemical substitutions at the N−C α position. This behavior fundamentally differs from electron transport in peptides, where H-bond interactions give rise to higher conductance states. All-atom molecular dynamics (MD) simulations are used to understand the conformational heterogeneity of peptoids, and molecular conformations obtained from MD simulations are used in quantum mechanical calculations based on the nonequilibrium Green’s function−density functional theory (NEGF-DFT) formalism. In all cases, computational results are in reasonable qualitative agreement with experiments. Our work demonstrates that the conductance behavior of peptoids depends on monomer identity, including side-chain aromaticity and substitution at the N−C α position. Overall, this work provides new insights into the structure−function relationships governing electron transport in peptoid-based materials and establishes design rules for peptoid-based molecular junctions.

Charge transport↗