Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular dynamic (MD) 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 433 records · Page 24

Modeling Nanoconfinement Effects Using Active Learning

Predicting the spatial configuration of gas in nanopores of is relevant in applications such as fluid flow forecasting and hydrocarbon reserves estimation. For example, shale reservoirs have suffered from computationally intractable multiscale problems, since fluid properties such as viscosity, density, and adsorption must be calculated by using expensive molecular dynamics (MD) simulations within each nanopore, whereas flow through these connected nanopores must be simulated at the micrometer scale. We utilize machine learning techniques to quickly and accurately model nanoscale confinement effects as an important step toward bridging the nano and micro scales. Our workflow is based on building and training physics-based deep-neural-networks models by learning from a database of MD calculations. The model accounts for the adsorption phenomenon by predicting the statistical distribution of gas inside nanopores. Because large databases of MD calculations are expensive to create, we investigate active learning (AL) as a data set construction strategy. In this workflow, new data are selected based on the model uncertainty via the query-by-committee approach. We show that our workflow obtains accurate models that generalize to real scanning electron microscopy geometries with 1/10th of the number of MD calculations required vs random data set generation. Our method enables the possibility of modeling nanoconfinement effects at the mesoscale, where complex connected sets of nanopores affect flow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Role of Phase Separation on Rayleigh-Plateau Type Instabilities in Alloys

Classical molecular dynamics (MD) simulations are used to investigate the role of phase separation (PS) on the Rayleigh-Plateau (RP) instability. Ni–Ag bulk structures are created at temperatures (2000 K and 1400 K) that generate different PS length scales, λ PS , relative to the RP instability length scale, λRP. Rectanguloids are then cut from the bulk structures and patterned with a perturbation of certain amplitude and wavelength, λ RP . It is found that when λ PS << λ RP (2000 K), the patterned rectanguloids break up into nanoparticles in a manner consistent with classical RP theory, whereas when λ PS << λ RP (1400 K), soluto-capillarity affects the RP instability significantly. Specifically, since Ag has a lower surface energy than Ni, Ag migrates to cover neighboring Ni regions, therefore modifying the RP instability. Thus, we demonstrate that the phase separation length scale of an immiscible alloy can be exploited to direct the assembly of functional bimetallic alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electric Field Effects on Water and Ion Structure and Diffusion at the Orthoclase (001)–Water Interface

Understanding the electrochemical properties of mineral–water interfaces tends to rely upon electrical double layer (EDL) models, but these models are based on the assumption that electrostatic equilibrium is constantly maintained. In reality, interfacial reactions, ion diffusion, and their electrochemical signatures are based in nonequilibrium conditions of locally or globally imbalanced electrical fields where current EDL models have limited purview. In this work, we performed molecular dynamics (MD) simulations of the orthoclase (001) surface in contact with a 1 M NaCl aqueous solution under various electric fields, to explore the interplay between EDL structure and dynamics when perturbed by electric fields of different direction and strength, by confinement, and by different distributions of structural surface charge. The simulations showed that confinement between two opposing (001) surfaces led to the development of an induced field when the applied field was perpendicular to the surfaces and, as a result, to ionic diffusion coefficients that were independent of electric field strength. In contrast, when the applied field was parallel to the surfaces, confinement resulted in ionic diffusion coefficients that were more strongly dependent on the magnitude of the electric field than in bulk water. Differences in the density and distribution of aluminol groups on the two surfaces had a significant impact on how the interfacial structure and dynamics varied in the presence of an electric field. Notably, these differences resulted in an electro-osmotic flow with opposite directions at the two surfaces under parallel applied electric field. Overall, the MD simulations highlighted the importance of considering atomic-level structure and heterogeneities when developing models of the electrochemical properties of mineral–water interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Water Structure and Dynamics Near the Surfaces of Silicalite-1

Zeolites are crystalline microporous aluminosilicates that are commonly used as industrial sorbents and membranes. Solvent structuring and dynamics near zeolite crystal surfaces are thought to influence interfacial transport and molecular recognition processes, as well as fundamental aspects of their crystallization from solution. Here, in this work, we use molecular dynamics (MD) simulations to investigate the behavior of interfacial water near the exposed (010), (100), and (101) crystal faces of silicalite-1, one of the most widely studied zeolites. The MD simulations reveal that water’s translational and orientational order is strongly influenced by the distinct corrugations, pore apertures, and functional group distributions presented on each surface. Specifically, we observe two distinct hydration layers near each surface. Water molecules in the contact layer are hydrogen bonded to two or three exposed surface silanol groups. The relative populations of water molecules that are doubly and triply hydrogen bonded to the surface strongly depend on the distribution of the exposed silanols. The interactions with exposed surface silanols also influence water’s local orientational order and distribution across each surface. We also show that the structuring of the solvent near the different faces of silicalite-1 strongly impacts the interfacial dynamics. The translational and orientational relaxation dynamics of water are slowest in the contact layers and correlate with structural ordering near each face.

42 ENGINEERING↗

Structural and Spectroscopic Properties of Butanediol-Modified Boehmite Materials

Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, in this study, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O–H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure and Dynamics of Aqueous Electrolytes at Quartz (001) and (101) Surfaces

Here, understanding and describing reactivity at mineral-water interfaces such as ion adsorption, the kinetics of dissolution, or surface charge development depends on our ability to improve the accuracy of electrical double layer (EDL) models. While molecular dynamics (MD) simulations are routinely used to investigate the structure and energetics of adsorbed ions comprising the EDL, less attention is paid to their self-diffusion dynamics, which can uniquely inform on coupling to interfacial reactions. Here we use MD to investigate both the organization and diffusion dynamics of water and electrolyte ions (NaCl, KCl, CaCl 2 ) at hydroxylated quartz (001) and (101) surfaces, a comparison which allowed us to assess surface structural effects of corrugation and silanol density. We found that inner- versus outer-sphere complex formation depends on cation size and charge but not necessarily hydration energies. Participation of surface silanols in the hydration spheres of Na + and K + generally indicated their preference for inner-sphere complexation, but this depends strongly on the orientation of the surface considered through its influence over the organization and dynamics of adsorbed water layers. In particular, surface orientation substantially affects the diffusive behavior of the near-surface water. Na + was found to decrease the mobility of water in the first layer, consistent with an increasing frequency of hydrolysis implied by faster quartz dissolution rates observed in experiments via the well known salt effect. Our results are also in good agreement with the observed dissolution rate of quartz vs. surface adsorption strength measure by Dove and Nix. This study sets the stage for a forthcoming paper examining how the dynamics at quartz/electrolyte interfaces are influenced by externally applied electric fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scale Invariance of Hot Spot Formation in TATB High Explosives

Shock-induced detonation of insensitive high explosives based on 1,3,5-triamino-2,4,6-trinitrobenzene starts with formation of hot spots at microstructural defects but has eluded atomistic modeling treatment at micron length scales. To this end, we performed multimicron scale all-atom molecular dynamics (MD) simulations of hot spots that form during the collapse of cylindrical pores with diameters between 10 and 300 nm. Our MD simulations show that hot spots formed at pores larger than 20 nm exhibit temperature fields with scale-invariant features for sizes up to at least 300 nm. Through a continuum-based grain-scale modeling framework, we span and extend beyond the size scales currently accessible to MD and find that hot spot scale invariance is a general feature that arises when the mechanical strength is insensitive to strain rate. Finally, our results demonstrate the applicability of all-atom MD to simulate the complicated dynamical evolution of micron-sized systems and bolster confidence in insights from MD simulations of materials that exhibit strength with negligible rate dependence over the relevant intervals.

36 MATERIALS SCIENCE↗

Effect of Hydrophilicity and Interfacial Water Structure on Particle Attachment

Oriented attachment (OA) of nanoparticles is an important and common crystallization pathway. Efforts to isolate and quantify the interparticle forces that underlie OA have only recently been initiated and the role played by the interfacial solvent structure in the confined region between nanocrystal surfaces in enabling and controlling OA remains an open question. In this work, classical molecular dynamics (MD) simulations were performed to compute the free energy landscape of nanoparticles approaching mineral substrates in both vacuum and pure water. Periclase (MgO), muscovite (KAl 3 Si 3 O 10 (OH) 2 ), and pyrophyllite (Al 2 Si 4 O 10 (OH) 2 ) nanoparticles exposing the (100), (001), and (001) surfaces, respectively, were brought to contact with substrates exposing the same surfaces. These mineral surfaces were selected to investigate a range of hydrophilicities (periclase > muscovite > pyrophyllite). The MD simulations revealed that the interfacial water structure created free energy minima for nanoparticle approach that corresponded to integer numbers of intervening water layers. The depth of the free energy minima and the free energy barriers between minima correlated positively with the hydrophilicity of the mineral surface. While the range of solvent-mediated interactions was predicted to be short (a few intervening water layers at most for MgO), the free energy barriers were calculated to be large (fraction of an eV per nm 2 ) and the effect of misalignment diminished rapidly with the number of intervening water layers. The results help explain the formation of mesocrystals and the observed ability of nanoparticles to rotate before attachment while in a long-lived solvent-separated state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Squeezing Out Interfacial Solvation: The Role of Hydrogen-Bonding in the Structural and Orientational Freedom of Molecular Self-Assembly

We report bioinspired membrane molecules with improved physical properties and enhanced stability can serve as functional models for conventional lipid or amphiphilic species. Importantly, these molecules can also provide new insights into emergent phenomena that manifest during self-assembly at interfaces. Here, we elucidate the structural response and mechanistic steps underlying the self-assembly of the amphiphilic, charged oligodimethylsiloxane imidazolium cation (ODMS-MIM+) at the air–aqueous interface using Langmuir trough methods with coincident surface-specific vibrational sum-frequency generation (SFG) spectroscopy. We find evidence for a new compression-induced desolvation step that precedes commonly known disordered-to-ordered phase transitions to form nanoscopic assemblies. The experimental data was supported by atomistic molecular dynamics (MD) simulations to provide a detailed mechanistic picture underlying the assembly and the role of water in these phase transitions. The sensitivity of the hydrophobic ODMS tail conformations to compression-owing to distinct water–ODMS interactions and tail–tail solvation properties-offers new strategies for the design of interfaces that can be further used to develop soft-matter electronics and low-dimensional materials using physical and chemical controls.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Binding Interactions between Receptor-Binding Domain of Spike Protein and Human Angiotensin Converting Enzyme-2 in Omicron Variant

The emergence of new SARS-CoV-2 Omicron variant of concern (OV) has exacerbated the COVID-19 pandemic because of a large number of mutations in the spike protein, particularly in the receptor-binding domain (RBD), resulting in highly contagious and/or vaccine-resistant strains. Herein, we present a systematic analysis based on detailed molecular dynamics (MD) simulations in order to understand how the OV RBD mutations affect the ACE2 binding. We show that the OV RBD binds to ACE2 more efficiently and tightly predominantly because of strong electrostatic interactions, thereby promoting increased infectivity and transmissibility compared to other strains. Some of the OV RBD mutations are predicted to affect the antibody neutralization either through their role in the S-protein conformational changes, such as S371L, S373P, and S375F, or through changing its surface charge distribution, such as G339D, N440K, T478K, and E484A. Other mutations, such as K417N, G446S, and Y505H, decrease the ACE2 binding, whereas S447N, Q493R, G496S, Q498R, and N501Y tend to increase it.

60 APPLIED LIFE SCIENCES↗

The Role of Ligand–Ligand Interactions in Multimodal Ligand Conformational Equilibria and Surface Pattern Formation

Multimodal chromatography uses multiple modes of interaction such as charge, hydrophobic, or hydrogen bonding to separate proteins. Recently, we used molecular dynamics (MD) simulations to show that ligands immobilized on surfaces can interact and associate with neighboring ligands to form hydrophobic and charge patches, which may have important implications for the nature of protein–surface interactions. Here, we study interfacial systems of increasing complexity—from a single immobilized multimodal ligand to high density surfaces—to better understand how ligand behavior is affected by the presence of a surface and the presence of other ligands in the vicinity, and how this behavior scales to larger systems. Furthermore, we find that tethering a ligand to a surface restricts its conformations to a subset of those observed in free solution, yet the ligand maintains flexibility in the plane of the surface and can form contacts with neighboring ligands. We find that although the formation of a contact between two neighboring ligands is slightly unfavorable, three neighboring ligands exhibit a preference for the formation of a fully connected cluster. To explore how these trends in ligand association extend to a larger surface with high density of ligands, we performed coarse-grained Monte Carlo (MC) simulations of a 132-ligand surface using ligand interactions parametrized based on free energies obtained from the three-ligand MD simulations. Despite their simplicity, the coarse-grained simulations qualitatively capture the cluster size distribution of ligands observed in detailed MD simulations. Quantitative differences between the two suggest opportunities for improvements in the coarse-grained energy function for efficient predictions of cluster and pattern formations. Our approach presents a promising route to the engineering of multimodal patterns for future chromatographic resin design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluating the Structural Response of Amphiphilic Monolayers to Environmental Stimuli

Amphiphilic monolayers composed of end groups with distinct polar and nonpolar functional groups offer rapid and reversible interfacial adaptation in response to environmental stimuli such as a change in interfacial medium polarity. We have synthesized and characterized a suite of monolayers with functional groups of competing polarity designed to reconfigure their interfacial chemical composition in response to solvent polarity. In these films, the end group is designed to be able to reorient and expose the functional groups that minimize the interfacial free energy between the film and the environment. Using a combination of spectroscopic, computational, and wettability studies, we have investigated the responsive interfacial behavior of different end groups upon exposure to environments with varying polarities. Contact angle measurements across a series of polar and dispersive probe liquids reveal trends that reflect the underlying molecular flexibility and composition. Vibrational sum frequency generation (SFG) spectroscopy and atomistic molecular dynamics (MD) simulations confirm solvent-driven reorientation of the end groups, with restructuring observed at the interface. To quantify these effects, we have developed a surface energy calculation model that incorporates solvent-induced surface rearrangements into the estimations. Our findings reveal a strong dependence of surface energy and switching behavior on the length and flexibility of the functionalities in the end group, which affects the exposure of certain interfacial compositions under different solvents. These results offer new insights into the design of adaptive monolayers and provide a framework for evaluating solvent-responsive surfaces.

functional groups↗

Surface Energies and Structure of Salt–Brine Interfaces

Permeability of salt formations is controlled by the equilibrium between the salt–brine and salt–salt interfaces described by the dihedral angle, which can change with the composition of the intergranular brine. In this work, classical molecular dynamics (MD) simulations were used to investigate the structure and properties of the salt–brine interface to provide insight into the stability of salt systems. Mixed NaCl–KCl brines were investigated to explore differences in ion size on the surface energy and interface structure. Nonlinearity was noted in the salt–brine surface energy with increasing KCl concentration, and the addition of 10% KCl increased surface energies by 2–3 times (5.0 M systems). Size differences in Na+ and K+ ions altered the packing of dissolved ions and water molecules at the interface, impacting the surface energy. Additionally, ions at the interface had lower numbers of coordinating water molecules than those in the bulk and increased hydration for ions in systems with 100% NaCl or 100% KCl brines. Ultimately, small changes in brine composition away from pure NaCl altered the structure of the salt–brine interface, impacting the dihedral angle and the predicted equilibrium permeability of salt formations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Percolated Ionic Aggregate Morphologies and Decoupled Ion Transport in Precise Sulfonated Polymers Synthesized by Ring-Opening Metathesis Polymerization

Here, we describe a set of precise single-ion conducting polymers that form self-assembled percolated ionic aggregates in glassy polymer matrices and have decoupled transport of metal cations. These precise single-ion conductors (SICs), synthesized by a scalable ring-opening metathesis polymerization, consist of a polyethylene backbone with a sulfonated phenyl group pendant on every fifth carbon and are fully neutralized by a counterion X + (Li + , Na + , or Cs + ). Experimental X-ray scattering measurements and fully atomistic molecular dynamics (MD) simulations are in good agreement. The MD simulations show that the ionic groups nanophase separate from the polymer backbone to form percolating ionic aggregates. Using graph theory, we find that within the Li + - and Na + -neutralized polymers the percolated aggregates exhibit planar and ribbon-like configurations at intermediate length scales, while the percolated aggregates within the Cs + -neutralized polymers are more isotropic. Electrical impedance spectroscopy measurements show that the ionic conductivities exhibit Arrhenius behavior, with conductivities of 10 –7 to 10 –6 S/cm at 180 °C. In the MD simulations, the cations move between sulfonate groups in the percolated aggregates, larger ions travel further, and overall cations travel further than the polymer backbones, indicating a decoupled ion-transport mechanism. Thus, the percolated ionic aggregates in these polymers can serve as pathways to facilitate decoupled ion motion through a glassy polymer matrix.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of Molecular Architecture on Ion Transport in Ethylene oxide-Based Polymer Electrolytes

This work aims to develop a detailed mechanistic understanding of the role of a graft polymer architecture on lithium ion (Li + ) transport in poly(ethylene oxide)-based polymer electrolytes. Specifically, we compare Li + transport in poly-(ethylene oxide) (PEO) versus poly(oligo oxyethylene methacrylate) (POEM) polymers doped with lithium bis( trifluoromethanesulfonyl) (LiTFSI) salts, using both experimental electrochemical characterization and molecular dynamics (MD) simulations. Furthermore, our results indicate that POEM exhibits a range of relaxation processes that cannot be interpreted solely in terms of glass-transition temperature (T g ) effects. Due to its side-chain architecture, the segmental relaxation of POEM is nonuniform across ether oxygens (EOs) and shows a more pronounced sensitivity to temperature above T g compared to PEO. Moreover, POEM also exhibits a nonuniform Li + coordination behavior, in which Li + is primarily solvated by two different chains in POEM, compared to a single chain in PEO. Li + transport in POEM occurs via two events with distinct characteristic times: a fast intrachain hopping along side chains and a slow interchain hopping between side chains. Taken together, the relaxation processes and ion transport mechanisms identified in POEM provide useful insights into design of more effective solid polymer electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantification of Mono- and Multivalent Counterion-Mediated Bridging in Polyelectrolyte Brushes

Multivalent counterion-induced bridging interactions have been identified as the key mechanism of drastic collapse of the height of polyelectrolyte (PE) brushes. In this article, we employ all-atom molecular dynamics (MD) simulations to quantify the bridging interactions in PE brushes for counterions of different sizes and valences. We identify that unlike the current notion, bridging interactions are not the sole function of the counterion valence. Rather the bridging interactions depend on the fraction of counterions (of a given type) that get physically condensed on the PE backbone as well as the size of the counterion solvation shell. These mechanisms ensure that certain monovalent counterions demonstrate much stronger bridging interactions than that witnessed for certain divalent and trivalent counterions, while certain counterions of identical valences show drastically different bridging. Here, we argue that these counterion-specific bridging interactions eventually enable not only the significant reduction of the PE brush height in presence of certain multivalent screening counterions, but may also give rise to scenarios where the brush height reduction for certain monovalent counterions is larger than certain divalent and trivalent counterions. This latter observation contradicts the experimental findings where the multivalent counterions invariably led to a larger decrease in the height of the PE brushes: we argue that this discrepancy stems from the fact that in our simulations we only consider densely grafted and short (and hence less flexible) PE brushes that hinder the formation of different laterally inhomogeneous structures (like pinned micelles and cylindrical bundles) that would have led to a larger brush height reduction (in experiments, which invariably consider longer and less densely grafted brushes, the formation of such inhomogeneous structures are primarily responsible for larger brush height reduction in presence of multivalent counterions). Finally, we also probe the dynamic properties of the counterions (i.e., their time-dependent displacements) and their bridging interactions (i.e., lifetime of bridging interactions).

42 ENGINEERING↗

Atomic Edge-Guided Polyethylene Crystallization on Monolayer Two-Dimensional Materials

In this report we combine an advanced synthesis of two-dimensional (2D) materials (MoSe 2 ) having well-defined atomic edge configurations with ab initio and atomistic molecular dynamics (MD) simulations to study how atomic edges interact with polyethylene (HDPE) chains in a dilute solution assembly process. Our results reveal that Mo-terminated zigzag (Mo-ZZ) edges act as preferred nucleation sites and strongly interact with HDPE chains. The HDPE chains align in parallel with the Mo-ZZ edges and form arrays of lamellae that are perpendicular to the edges. Interestingly, atomic edge configurations are observed to dramatically change such interactions. The crystallization discrepancy at different edges was demonstrated on the same piece of MoSe 2 with different types of edges. The ab initio and MD simulations between n-alkane (n = 5 and 25), a segment of HDPE, and MoSe 2 suggest that the atomic structures of MoSe 2 can affect their interactions with n-alkane chains. Following the Mo-ZZ edge preferred nucleation principle, controlled long-range alignment of HDPE lamellae can be realized by creating multilayer MoSe 2 with parallel atomic steps. This research opens a pathway toward an atomic level understanding of polymer–2D nanomaterial interactions. It also bridges the gap between atomic-level and long-range mesoscopic structures and introduces a novel strategy for long-range structural control.

36 MATERIALS SCIENCE↗

Effects of Ionic Group Distribution on the Structure and Dynamics of Amorphous Polymer Melts

Ionizable groups tethered to a polymer backbone often associate to form ionic clusters, whose features are determined by the balance of the polymer backbone and electrostatics. These assemblies impact the dynamics of the macromolecules and their ability to transport ions. Here, using fully atomistic molecular dynamics (MD) simulations, we investigate the effects of the distribution of ionizable groups along the polymer backbone on cluster characteristics and the resulting impacts on the structure and dynamics of amorphous polymers. Particularly, we probe polystyrene sulfonates (PSS) with random, precise, and block configurations of the SO 3 – sulfonate groups along the backbone with Na + as the counterion. We find that the distribution of the ionic groups affects the shape and distribution of the clusters as well as the internal packing of the ionizable groups in the cluster and the number of unique chains that participate in each cluster, affecting the structure and the dynamics of the polymers. Furthermore, the signature of ionic clusters, observed in the static structure factor S(q) for all three distributions, is significantly more pronounced for the precise and blocky polymers compared to the random one. Remarkably, we find that the local mobility of the polymer segments is not only affected by the number and size of the clusters but also by the number of polymer chains associated with clusters.

36 MATERIALS SCIENCE↗