Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular dynamics 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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Solvent Isotherms and Structural Transitions in Nanoparticle Superlattice Assembly

Here, we introduce a Molecular Theory for Compressible Fluids (MOLT-CF) that enables us to compute free energies and other thermodynamic functions for nanoparticle superlattices with any solvent content, including the dry limit. Quantitative agreement is observed between MOLT-CF and united-atom molecular dynamics simulations performed to assess the reliability and precision of the theory. Among other predictions, MOLT-CF shows that the amount of solvent within the superlattice decreases approximately linearly with its vapor pressure and that in the late stages of drying, solvent-filled voids form at lattice interstitials. Applied to single-component superlattices, MOLT-CF predicts fcc-to-bcc Bain transitions for decreasing vapor pressure and for increasing ligand length, both in agreement with experimental results. We explore the stability of other single-component phases and show that the C14 Frank-Kasper phase, which has been reported in experiments, is not a global free-energy minimum. Implications for precise assembly and prediction of multicomponent nanoparticle systems are discussed.

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

Pivotal Roles of Triple Screening-Topological, Electrostatic, and Hydrodynamic-On Dynamics in Semidilute Polyelectrolyte Solutions

For semidilute polyelectrolyte solutions, it is generally assumed that topological, electrostatic, and hydrodynamic interactions are screened (called triple screening). Despite a large body of research focused on polyelectrolyte solutions, the concept of triple screening has never been rigorously verified. In this work, we test the concept by probing concentration fluctuations in aqueous solutions containing a well-studied polyelectrolyte, sodium poly(styrenesulfonate) (NaPSS) with neutron scattering, theory, and molecular dynamics simulations. Neutron spin–echo (NSE) and small-angle neutron scattering (SANS) data from semidilute solutions of NaPSS are presented at different polymer and salt (NaCl) concentrations. A combined theory for structure (J. Chem. Phys. 105, 5183 (1996)) and dynamics (J. Chem. Phys. 107, 2619 (1997)), which captures effects of hydrodynamic, topological, and electrostatic screening, is used to interpret the experimental results. Here, the theory quantitatively predicts the decay rate obtained from the NSE measurements while capturing the shape and concentration dependencies of the polyelectrolyte peak observed in the SANS spectra. Detailed comparisons of the theory and the experiments reveal that the wavevector-dependent decay rate of concentration fluctuations in semidilute solutions of polyelectrolytes is dictated by the screening of hydrodynamic, topological, and electrostatic interactions. This conclusion is corroborated by coarse-grained molecular dynamics simulations, executed without any hydrodynamic interactions, which fail to capture the correct wavevector dependence of the decay rate. These results highlight that the theories based on the concept of triple screening provide a quantitative framework for predicting a relation between the structure and dynamics of polyelectrolyte solutions.

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

Solution Structure and Scaling Laws of Cylindrical and Tapered Bottlebrush Polymers

Bottlebrush polymers are a unique class of macromolecular architectures with a plethora of potential industrial and pharmaceutical applications that critically depend on the bottlebrush shape and dimensions. Here, a systematic series of 12 cylindrical and 12 cone-shaped (tapered) bottlebrush polymers with poly(tert-butyl acrylate) (PtBA) or polystyrene (PS) side chains were synthesized using the sequential addition of macromonomers ring-opening metathesis polymerization (SAM-ROMP) grafting-through method. Small-angle neutron scattering (SANS) studies on dilute solutions of the two types of bottlebrush polymers provided noninvasive characterization of their structural dimensions and chain conformations. Simulated SANS traces, generated using coarse-grained molecular dynamics simulations, reproduced the distinctive features observed in the experimental SANS signals and provided necessary validation for data modeling. The combined analysis of experimental and simulated SANS signals yielded key structural and conformational parameters, including the bottlebrush radius, length, and Kuhn length as well as the excluded volume parameter and the correlation length of the polymer side chains. Importantly, the obtained structural parameters followed well-defined scaling laws as a function of the backbone and side chain degrees of polymerization, as predicted by mean field theories. In conclusion, these findings provide clear experimental and computational evidence of the interdependence of structural and conformational properties in an important class of polymer architectures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unveiling the Mechanism of Mn Dissolution Through a Dynamic Cathode‐Electrolyte Interphase on LiMn2O4

Abstract Understanding the formation and evolution of the cathode‐electrolyte interphase (CEI), which forms at the interface between the cathode and electrolyte, is crucial for revealing degradation mechanisms in cathode materials, especially for developing strategies to stabilize the interphase in the strongly oxidizing conditions that evolve at high operating voltages in next‐generation Li‐ion batteries. However, The present understanding of the CEI is challenged by its complex and dynamic nature. In this work, near‐edge X‐ray absorption fine structure spectroscopy, electrochemical characterization, and reactive molecular dynamics simulations are combined to reveal a mechanism for CEI formation and evolution above model LiMn 2 O 4 (LMO) thin‐film electrodes in contact with conventional carbonate‐based electrolytes. It is found that Mn dissolution from LMO can be understood in terms of repetitive Mn 3 O 4 formation and dissolution behavior during cycling, which is closely connected to electrolyte decomposition and a key aspect of the CEI formation and growth. The behavior of the CEI in this model system offers detailed insight into the dynamic chemistry of the interphase, underscoring the important role of electrolyte composition and cathode surface structure in interphase degradation.

Ou, Wenhan↗
DOE OSTI2024

Effect of K + Force Fields on Ionic Conductivity and Charge Dynamics of KOH in Ethylene Glycol

Predicting ionic conductivity is crucial for developing efficient electrolytes for energy storage and conversion and other electrochemical applications. An accurate estimate of ionic conductivity requires understanding complex ion–ion and ion–solvent interactions governing the charge transport at the molecular level. Molecular simulations can provide key insights into the spatial and temporal behavior of electrolyte constituents. However, such insights depend on the ability of force fields to describe the underlying phenomena. In this work, molecular dynamics simulations were leveraged to delineate the impact of force field parameters on ionic conductivity predictions of potassium hydroxide (KOH) in ethylene glycol (EG). Four different force fields were used to represent the K + ion. Diffusion-based Nernst–Einstein and correlation-based Einstein approaches were implemented to estimate the ionic conductivity, and the predicted values were compared with experimental measurements. The physical aspects, including ion-aggregation, charge distribution, cluster correlation, and cluster dynamics, were also examined. A force field was identified that provides reasonably accurate Einstein conductivity values and a physically coherent representation of the electrolyte at the molecular level.

25 ENERGY STORAGE↗
DOE OSTI2026

Free Energy and Flexibility Analysis of Autoinhibited Human BRAF

The RAF serine/threonine protein kinases function as direct effectors of RAS in the intracellular transmission of extracellular growth signals, and they are key targets for drug discovery, given the high incidence of oncogenic mutations in RAF and other components of this signaling pathway. In its inactive state, RAF is held in an autoinhibited conformation in the cytosol through a combination of intramolecular interactions and binding to a regulatory 14−3−3 protein dimer. Activation of RAF is initiated by its interaction with membrane-localized GTP-bound RAS, which induces conformational changes that release RAF from its autoinhibited state. However, the molecular mechanisms governing RAF activation remain incomplete, largely due to the challenges in experimentally capturing the intermediate conformational states in this process. To address this gap, we developed a comprehensive all-atom model of BRAF based on existing cryo-EM structures. Using this model, we performed extensive molecular dynamics simulations to evaluate the stability and free energy landscape of autoinhibited BRAF in solution. Our analysis reveals conformational flexibility within the autoinhibited complex, suggesting that this dynamic behavior may play a role in facilitating BRAF activation upon engagement with the membrane-bound RAS.

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

Roles of kaolinite-oil-gas molecular interactions in hydrogen storage within depleted reservoirs

Hydrogen is a clean alternative to fossil fuels, emitting only water vapor during combustion. In a future hydrogen economy, large-scale storage will be an important component of the supply chain. Due to its low volumetric density, conventional surface storage methods are inadequate. Underground hydrogen storage (UHS) offers a viable solution, enabling the storage of millions of cubic meters. Among potential geological sites, including salt caverns, aquifers, and depleted gas reservoirs, depleted oil reservoirs show promise. Studying hydrogen interactions with residual oil and reservoir minerals is vital for understanding the properties of hydrogen and the reservoir post-injection. In this work, we employed molecular dynamics simulations to gain molecular-level insights into these interactions. We investigated hydrogen dissolution in oil, adsorption at kaolinite/oil interfaces, the role of CO 2 as a cushion gas, and the influence of kaolinite’s hydrophobicity on H 2 behavior. The main findings include: (1) hydrogen dissolves more in oil than in water, (2) the introduction of CO 2 suppresses hydrogen dissolution in oil and reduces the interfacial tension (IFT) between oil and gas, (3) CO 2 decreases H 2 partitioning near kaolinite surfaces due to its strong affinity for the hydrophilic gibbsite surface of kaolinite, and (4) CO 2 is more effective than H 2 in reducing IFT between kaolinite and the oil–gas mixture. These findings emphasize the effectiveness of CO 2 as a cushion gas and the important role of clay hydrophobicity in UHS, providing insights that are challenging to obtain experimentally.

08 HYDROGEN↗
DOE OSTI2026

Computational methods in solution-based plastics purification

Plastic waste can be recycled into resins with near-virgin properties by solution-based purification processes that selectively dissolve polymers, remove contaminants, or detach printing residues. Here, in this review, we examine computational methods for predicting the behavior governing solution-based plastic purification, motivated by the vast polymer–solvent–contaminant compositional space. We discuss thermodynamic and machine learning methods for predicting polymer–solvent and polymer–contaminant interaction and review physics-based molecular dynamics simulations that resolve molecular-scale phenomena within polymer matrices inaccessible to screening methods. We highlight how these methods have informed experimental design for dissolution-based recycling and solvent-based contaminant removal. Finally, we discuss the prospective role of agentic AI in integrating these computational tools with real-time sorting data to adapt purification conditions to the compositional variability of real post-consumer feedstocks. This review charts a path toward computationally guided solution-based purification workflows that can respond to the complexity inherent in plastic waste streams.

Altamimi, Ali [Univ. of Wisconsin, Madison, WI (Un↗
DOE OSTI2021

Molecular-level understanding of gibbsite particle aggregation in water

Using molecular dynamics simulations, we investigate the molecular scale origin of crystal face selectivity when one gibbsite particle attaches to another in water. A comparison of the free energy per unit surface area of particle–particle attachment indicates that particle attachment through edge surfaces, where the edge surfaces are either (1 0 0) or (1 1 0) crystal faces, is more energetically favorable compared to attachment between two basal surfaces (i.e., (0 0 1) crystal faces) or between the basal surface of one particle and the edge surface of another. This result suggests that gibbsite crystals with low basal/edge surface area ratio will preferentially attach through edge surfaces, potentially helping the crystals grow laterally. However, for larger gibbsite particles (high basal/edge surface area ratio) the total free energy, not normalized by surface area, of particle attachment through the basal surfaces is lower (more negative) than attachment through the edge surfaces, indicating that larger gibbsite particles will preferentially aggregate through basal surface attachments. The short-range electrostatic interactions including the interparticle hydrogen bonds from surface –OH groups drive particle attachment, and the dominant contribution to the free energy minimum is enthalpic rather than entropic. However, the enthalpy of basal-edge attachment is significantly offset by the entropy leading to a higher free energy (less negative) compared to that of basal-basal attachment. Study of the free energy for a few imperfect attachments of two particles indicates a higher free energy (i.e., less negative, less stable), compared to a perfect attachment

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

Microstructural impacts on ionic conductivity of oxide solid electrolytes from a combined atomistic-mesoscale approach

Abstract Although multiple oxide-based solid electrolyte materials with intrinsically high ionic conductivities have emerged, practical processing and synthesis routes introduce grain boundaries and other interfaces that can perturb primary conduction channels. To directly probe these effects, we demonstrate an efficient and general mesoscopic computational method capable of predicting effective ionic conductivity through a complex polycrystalline oxide-based solid electrolyte microstructure without relying on simplified equivalent circuit description. We parameterize the framework for Li 7- x La 3 Zr 2 O 12 (LLZO) garnet solid electrolyte by combining synthetic microstructures from phase-field simulations with diffusivities from molecular dynamics simulations of ordered and disordered systems. Systematically designed simulations reveal an interdependence between atomistic and mesoscopic microstructural impacts on the effective ionic conductivity of polycrystalline LLZO, quantified by newly defined metrics that characterize the complex ionic transport mechanism. Our results provide fundamental understanding of the physical origins of the reported variability in ionic conductivities based on an extensive analysis of literature data, while simultaneously outlining practical design guidance for achieving desired ionic transport properties based on conditions for which sensitivity to microstructural features is highest. Additional implications of our results are discussed, including a possible connection between ion conduction behavior and dendrite formation.

25 ENERGY STORAGE↗
DOE OSTI2025

In-Silico Design of Next Generation Cellulose-Derived Packaging Materials (CRADA Final Report)

Developing sustainable solutions for single-use packaging is an important objective to combat the environmental crisis of plastics pollution. Most embodiments of cellulose-based packaging materials, including CellophaneTM, are completely biodegradable in both terrestrial and marine environments. However, petroleum-derived alternatives offer some performance advantages for metrics such as moisture barriers and mechanical properties. This project leverages molecular dynamics simulation to investigate how molecular modifications to cellulose-based polymer assemblies impact their material properties. An important performance criterion for the modified materials was to retain biodegradability; thus, modifications by naturally occurring, biodegradable additives were the focus of this study. Specifically, we developed models with xylan and lignin of varying monomeric compositions into the cellulose matrix. The mechanical properties were investigated by performing stress-strain simulations, and the water barrier and hydrophobicity were investigated by simulating the water contact angle. Our findings indicate that the incorporation of xylan into the cellulose matrix tends to increase the mechanical properties with an optimal loading of ~27 wt%. We also predict that orienting the nanoscale directionality of the xylan chains such that they are perpendicular to the cellulose fibrils will dramatically increase mechanical strength. In contrast, the incorporation of lignin tends to weaken the composite at all loadings investigated. Simulations of water contact angle predicted that coating polymers on the surface of the cellulose assembly creates a more hydrophobic surface than incorporating them throughout the matrix. Of the coatings investigated, lignin resulted in the most hydrophobic surface, followed by pectin and keratin, which both imparted modest increases in hydrophobicity. Future experimental work done by Futamura will focus on designing material prototypes to capitalize on the predictions of performance enhancement obtained from molecular modeling. While substantial progress was made by the simulations performed in this project, there still exists a vast parameter space that we were unable to investigate, including branching, functional group decoration, and degree of polymerization of polymer additives. However, the methods developed in this initial investigation will facilitate more rapid evaluation of the impact of molecular characteristics on the performance of biopolymer composite materials and thereby accelerate future materials discovery efforts in this area.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Stockmayer fluid simulations for viscosity and glass transition temperature of ionic liquids

We develop a Stockmayer fluid model for molecular dynamics simulations of ionic liquids that captures molecular polarization, ionic conductivity, viscosity, and glass transition temperature, using ethylammonium nitrate (EAN) as an example. The ions in EAN are treated as spheres interacting via the Lennard-Jones potential with an embedded point charge and a permanent dipole moment. We show that our simulation results for EAN are consistent with experimental data and then explore the effects of the molecular parameters on the viscosity of ionic liquids. Our results indicate that viscosity monotonically increases with ionic charge and dipole moment but non-monotonically changes with ionic diameter (or molar volume). This non-monotonic trend arises from the competition among the electrostatic interactions, molecular packing, and size asymmetry between the cation and anion. In conclusion, our model also shows that long-lived ion pairs result in higher viscosities.

Coarse-grained simulations↗