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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Hotspot’s Better Half: Non-Equilibrium Intra-Molecular Strain in Shock Physics

Shockwave interactions with a material’s microstructure localizes energy into hotspots, which act as nucleation sites for complex processes such as phase transformations and chemical reactions. To date, hotspots have been described via their temperature fields. Nonreactive, all-atom molecular dynamics simulations of shock-induced pore collapse in a molecular crystal show that more energy is localized as potential energy (PE) than can be inferred from the temperature field and that PE localization persists beyond thermal diffusion. Furthermore, the origin of the PE hotspot is traced to large intramolecular strains, storing energy in modes readily available for chemical decomposition.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Extreme Metastability of Diamond and its Transformation to the BC8 Post-Diamond Phase of Carbon

Diamond possesses exceptional physical properties due to its remarkably strong carbon–carbon bonding, leading to significant resilience to structural transformations at very high pressures and temperatures. Despite several experimental attempts, synthesis and recovery of the theoretically predicted post-diamond BC8 phase remains elusive. Through quantum-accurate multimillion atom molecular dynamics (MD) simulations, we have uncovered the extreme metastability of diamond at very high pressures, significantly exceeding its range of thermodynamic stability. We predict the post-diamond BC8 phase to be experimentally accessible only within a narrow high pressure–temperature region of the carbon phase diagram. The diamond to BC8 transformation proceeds through premelting followed by BC8 nucleation and growth in the metastable carbon liquid. Here we propose a double-shock compression pathway for BC8 synthesis, which is currently being explored in experiments at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Molecular Picture of the Transient Nature of Lipid Rafts

In biological membranes, lipid rafts are now thought to be transient and nanoscopic. However, the mechanism responsible for these nanoscopic assemblies remains poorly understood, even in the case of model membranes. Therefore, it has proven extremely challenging to probe the physicochemical properties of lipid rafts at the molecular level. Here, we use all-atom molecular dynamics (MD) simulations and inelastic X-ray scattering (IXS), an intrinsically nanoscale technique, to directly probe the energy transfer and collective short-wavelength dynamics (phonons) of biologically relevant model membranes. We show that the nanoscale propagation of stress in lipid rafts takes place in the form of collective motions made up of longitudinal (compression waves) and transverse (shear waves) molecular vibrations. Importantly, we provide a molecular picture for the so-called van der Waals mediated "force from lipid" [Anishkin, A. et al. Proc. Natl. Acad. Sci. U.S.A. 2014, 111, 7898], a key parameter for the ionic channel mechano-transduction and the mechanism for the lipid transfer of molecular level stress [Aponte-Santamaría, C. et al. J. Am. Chem. Soc. 2017, 139, 13588]. Specifically, we describe how lipid rafts are formed and maintained through the propagation of molecular stress, lipid raft rattling dynamics, and a relaxation process. Eventually, the rafts dissipate through the self-diffusion of lipids making up the rafts. We also show that the molecular stress and viscoelastic properties of transient lipid rafts can be modulated through the use of hydrophobic biomolecules such as melatonin and tryptophan. Overall, the herein proposed mechanism describing the molecular interactions for the formation and dissolution of lipid rafts may offer insights as to how lipid rafts enable biological function.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Adsorption, Orientation, and Speciation of Amino Acids at Air–Aqueous Interfaces for the Direct Air Capture of CO 2

Amino acids make up a promising family of molecules capable of direct air capture (DAC) of CO 2 from the atmosphere. Under alkaline conditions, CO 2 reacts with the anionic form of an amino acid to produce carbamates and deactivated zwitterionic amino acids. The presence of the various species of amino acids and reactive intermediates can have a significant effect on DAC chemistry, the role of which is poorly understood. In this study, all-atom molecular dynamics (MD) based computational simulations and vibrational sum frequency generation (vSFG) spectroscopy studies were conducted to understand the role of competitive interactions at the air–aqueous interface in the context of DAC. We find that the presence of potassium bicarbonate ions, in combination with the anionic and zwitterionic forms of amino acids, induces concentration and charge gradients at the interface, generating a layered molecular arrangement that changes under pre- and post-DAC conditions. In parallel, an enhancement in the surface activity of both anionic and zwitterionic forms of amino acids is observed, which is attributed to enhanced interfacial stability and favorable intermolecular interactions between the adsorbed amino acids in their anionic and zwitterionic forms. The collective influence of these competitive interactions, along with the resulting interfacial heterogeneity, may in turn affect subsequent capture reactions and associated rates. Finally, these effects underscore the need to consider dynamic changes in interfacial chemical makeup to enhance DAC efficiency and to develop successful negative emission and storage technologies.

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

Relaxor Ferroelectric-Like Spatiotemporal Memory in Field-Driven Lipid Bilayers

Lipid membranes are often regarded as passive barriers, yet their nonlinear dielectric response remains poorly understood. Using all-atom molecular dynamics, we show that fully hydrated dipalmitoylphosphatidylcholine bilayers exhibit relaxor ferroelectric-like behavior under time-dependent electric fields. Unlike crystalline relaxors, which are bipolar and display little remanent polarization, lipid bilayers exhibit a unipolar polarization response: even an alternating current field produces persistent, asymmetric polarization. Furthermore, the underlying free-energy landscape contains two distinct minima, a nonpolarized state and a unipolarly polarized state, between which stochastic thermally activated transitions occur. Directionally resolved Van Hove analysis reveals pronounced anisotropy arising from out-of-plane electric dipole alignment, interleaflet coupling, and lateral polarization domains. Each field cycle nucleates polarization at distinct sites and monitors their relaxation, marking a crossover from thermal fluctuations to field-sustained polarization. Remarkably, these polarized domains persist after field removal, generating long-lived, spatially coherent dipolar patterns that encode nanoscale polarization memory. Potassium chloride amplifies these effects via dielectric screening and a modified hydration structure, enhancing electric dipole flexibility and cooperativity. Together, these results establish protein-free bilayers as nonlinear, history-dependent dielectrics capable of sustaining field-tunable electromechanical coupling, providing an emergent physical foundation for nanoscale information storage and memory phenomena reminiscent of short- and long-term plasticity in soft neuromorphic systems.

Insulators↗
DOE OSTI2021

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↗
DOE OSTI2022

Charge-Density-Specific Response of Grafted Polyelectrolytes to Electric Fields: Bending or Tilting?

The response of polyelectrolytes (PEs) to applied electric fields drive applications in energetics, diagnostics, materials development, and many more. Here we employ all-atom molecular dynamics (MD) simulations to probe the response of grafted PE brushes to axial electric fields. For PEs with large charge densities, the electric field triggers a left-right asymmetric distribution of counterions around the PE backbone: consequently, depending on the location (left or right), there is an unequal screening and an unequal force on the PE functional groups causing a bending-driven brush height reduction. However, for the weakly charged PEs, the electric field causes a uniform distribution of the counterions across the brush enforcing a uniform (without left-right asymmetry) partially unscreened PE brush: therefore, all the brush segments experience similar force causing a brush tilting-driven brush height reduction. As a result, such bending versus tilting responses is commensurate with the electric-field-driven increase (decrease) in the flexibility of strongly (weakly) charged PEs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2022

Specific Ion Effects in Lanthanide–Amphiphile Structures at the Air–Water Interface and Their Implications for Selective Separation

The use of surfactants to attract dissolved ions to water surfaces and interfaces is an essential step in both solvent-based and solvent-free separation processes. Here, we have studied the interactions of lanthanide ions in the aqueous subphase with monolayers of dihexadecyl phosphate at air–water interfaces. With heavier lanthanides (atomic number Z ≥ 65) in the subphase, the floating layer can be compressed to an area/molecule of about half the molecular cross section, indicating bilayer formation. X-ray fluorescence and reflectivity data support this conclusion. In the presence of lighter lanthanides (Z < 65), only monolayers are observed. Subphase-concentration-dependent studies using Er 3+ (heavier) and Nd 3+ (lighter) lanthanides show a stepwise progression, with ions attaching to the monolayer only when the solution concentration is >3 × 10 –7 M. Above ~10 –5 M, bilayers form but only in the presence of the heavier lanthanide. Grazing incidence X-ray diffraction shows evidence of lateral ion–ion correlations in the bilayer structure but not in monolayers. Explicit solvent all-atom molecular dynamics simulations confirm the elevated ion–ion correlation in the bilayer system. This bilayer structure isolates heavier lanthanides but not lighter lanthanides from an aqueous solution and is therefore a potential mechanism for the selective separation of heavier lanthanides.

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

Hydration and Restructuring of Polar Polymer Interfaces: Implications in Antifouling and Responsive Materials

Manipulating polymer interfaces is crucial for understanding how structure influences function in applications spanning biofouling prevention to energy storage. Moreover, observing how polymers adapt their microscopic structure to changes in their local environment can reveal essential properties that govern their performance in such applications, providing key insights into how to design more effective interfaces. Here, in this study, a series of “grafting-from” polymer brushes with side chains varying in charge, are probed by sum frequency generation (SFG) and modeled using all-atom molecular dynamics (MD) simulations to elucidate how chemical makeup and charge mediate interfacial restructuring in dry versus hydrated states. Results show that charge, in progressing from nonpolar to cationic to zwitterionic, results in dramatic changes in interfacial structure and overall hydration. While net neutral systems, regardless of bulk phase polarity, show minimal interfacial water structuring, the cationic species exhibits strong bulk water signals from the surface potential. Meanwhile, the polymer brushes themselves restructure in water differently independent of charge, impacting the functional groups that are presented to the aqueous phase. Nonpolar and cationic species for instance undergo a change in alkyl group orientations to accommodate hydrating water molecules, whereas the zwitterionic polymer becomes completely disordered in water. Overall, the structure-based behavior trends presented herein have implications in antifouling applications and responsive material interfaces.

biointerface↗
DOE OSTI2025

Antielectrophoretic Response-Driven Bending–Tilting Deformation of Cationic Polyelectrolyte Brushes Drives Nonlinear Electroosmotic Transport in Brush-Grafted Nanochannels

In this paper, we use all-atom molecular dynamics (MD) simulations to describe a non-linearly enhanced electroosmotic (EOS) flow, where, in a nanochannel grafted with cationic PMETAC ([Poly(2-(Methacryloyloxy)Ethyl) Trimethylammonium Chloride]) brushes, a two-fold increase in the electric field strength leads to a several-fold (more than two-fold) increase in the EOS flow strength and volume flow rate. The electric field enforces the PMETAC brushes to undergo a bending-tilting driven deformation with a significant portion of the brush layer becoming parallel to the grafting surface. In response, a substantial fraction of the counterions leave the brush layer (hence become more mobile), but instead of going into the bulk, accumulate at the brush-bulk interface, i.e., stay in proximity of the brush segments aligned parallel to the grafting surface. This creates an interesting situation, where the counterions are not completely within the brush layer, yet they fully screen the brush charges. Such “freer” conditions enable the counterions to achieve very high velocity, thereby ensuring that the water solvating the counterions themselves move very fast triggering the significantly augmented EOS transport. Probing deeper we can identify that the bending-tilting driven brush deformation, enforcing the brushes to align parallel to the substrate, results from the anti-electrophoretic behavior of the brushes, where despite being positively charged, the brushes move against the electric field direction. Such an anti-electrophoretic behavior of the PE brushes, which has not been reported before, can be associated with the very fast velocities of the negatively charged counterions and the electrostatic and hydrodynamic coupling of the counterions with the positive functional groups of the brushes. Here, we anticipate that the findings of this paper will shed light on strategies for nanochannel flow, the anti-electrophoretic response of charged polymer chains, and the significance of capturing the detailed chemical architecture of polyelectrolytes in nanoscale science and engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Melatonin-Induced Modulation of Cholesterol-Enriched Model Neuronal Membranes

Melatonin, a hormone primarily produced by the brain’s pineal gland, not only regulates circadian rhythms, but also influences the structural and biophysical properties of neuronal membranes. Its amphiphilic nature enables direct incorporation into lipid bilayers and preferential interactions with cholesterol-rich lipid rafts, critical hubs for cellular signaling and membrane organization. Despite increasing recognition of its membrane activity, the molecular basis of melatonin’s interactions with coexisting liquid-ordered (L o ) and liquid-disordered (L d ) phases remains unclear. Here, in this study, we combine small-angle neutron scattering (SANS) and all-atom molecular dynamics simulations to examine model neuronal membranes composed of DSPC, DOPC, POPC, and cholesterol. Our results show that melatonin preserves domain morphology while adopting distinct orientations within the bilayer and at the membrane interface, allowing both lateral and transmembrane bridging across lipid phases. These findings establish the molecular underpinnings of melatonin’s modulation of membrane heterogeneity and provide strong support for its receptor-independent actions.

atomistic simulations↗
DOE OSTI2021

Overscreening, Coion-Dominated Electroosmosis, and Electric Field Strength Mediated Flow Reversal in Polyelectrolyte Brush Functionalized Nanochannels

Controlling the direction and strength of nanofluidic electrohydrodyanmic transport in the presence of an externally applied electric field is extremely important in a number of nanotechnological applications. Here, we employ all-atom molecular dynamics simulations to discover the possibility of changing the direction of electroosmotic (EOS) liquid flows by merely changing the electric field strength in a nanochannel functionalized with polyelectrolyte (PE) brushes. In exploring this, we have uncovered three facets of nanoconfined PE brush behavior and resulting EOS transport. First, we identify the onset of an overscreening effect: such overscreening refers to the presence of more counterions (Na + ) within the brush layer than needed to neutralize the negative brush charges. Accordingly, as a consequence of the overscreening, in the bulk liquid outside the brush layer, there is a greater number of co-ions (Cl – ) than counterions in the presence of an added salt (NaCl). Second, this specific ion distribution ensures that the overall EOS flow is along the direction of motion of the co-ions. Such co-ion-dictated EOS transport directly contradicts the notion that EOS flow is always dictated by the motion of the counterions. Finally, for large-enough electric fields, the brush height reduces significantly, causing some of the excess overscreening-inducing counterions to squeeze out of the PE brush layer into the brush-free bulk. As a result, the overscreening effect disappears and the number of co-ions and counterions outside the PE brush layer become similar. Despite that there is an EOS transport, this EOS transport, unlike the standard EOS transport that occurs due to the imbalance of the co-ions and counterions, occurs since a larger residence time of the water molecules in the first solvation shell of the counterions (Na + ) ensures a water transport in the direction of motion of the counterions. In this work, the net effect is the reversal of the direction of the EOS transport by merely changing the strength of the electric field.

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

Simultaneous Energy Generation and Flow Enhancement (Electroslippage Effect) in Polyelectrolyte Brush Functionalized Nanochannels

Energy generation through nanofluidics is a topic of great nanotechnological relevance. In this work, we conduct all-atom molecular dynamics (MD) simulations of the transport of water and ions in a pressure-driven flow in nanochannels grafted with charged polyelectrolyte (PE) brushes and discover the possibility of simultaneous electrokinetic energy generation and flow enhancement (henceforth denoted as the electroslippage effect). Such PE-brush-functionalized nanochannels have been recently shown to demonstrate an overscreening (OS) effect (characterized by the presence of a greater number of screening counterions within the PE brush layer than needed to screen the PE brush charges), a consequent presence of excess coions within the PE brush-free bulk, and a coion-driven electroosmotic (EOS) transport in the presence of small to moderate applied axial electric fields. Here, however, we find that the streaming current, which represents the current generated by the flow-driven downstream advection of the charge imbalance present within the electric double layer (EDL) that screen the PE brush charges, is governed by the migration of the counterions. This stems from the fact that the highest contribution to the overall streaming current arises from the region near the PE brush-water interface (where there is an excess of counterions), while the brush-free bulk yields a hitherto unreported, but small, coion-dictated streaming current. This downstream advection of the charge imbalance (and the resultant counterion-driven streaming current) eventually leads to the development of an electric field (streaming electric field) in the direction that is opposite to the direction of the counterion-driven streaming current. The streaming current and the streaming electric field interact to generate the electrokinetic energy. Equally importantly, this streaming electric field induces an EOS transport, which becomes coion-driven, due to the presence of excess coions in the brush-free bulk. For the case of nanochannels grafted with negatively charged PE brushes, the streaming electric field will be in a direction that is opposite to the pressure-driven transport, and hence the coion (or anion) driven EOS flow, will be in the same direction as the pressure-driven transport. On the other hand, for the case of nanochannels grafted with positively charged PE brushes, the streaming electric field will be in the same direction as the pressure-driven flow, and hence the coion (or cation) driven EOS flow, will again be in the same direction as the pressure-driven flow. Therefore, whenever there occurs the presence of the OS and the resulting coion-driven EOS transport in PE brush grafted nanochannels, regardless of the sign of the charges of the PE brushes, this EOS transport will always aid the pressure-driven transport and will cause the most fascinating increase in the net volume flow rate across the nanochannel cross section, which is the electroslippage effect.

overscreening↗
DOE OSTI2020

Evolution of Metastable Structures at Bimetallic Surfaces from Microscopy and Machine-Learning Molecular Dynamics

The restructuring of interfaces plays a crucial role in materials science and heterogeneous catalysis. Bimetallic systems, in particular, often adopt very different compositions and morphologies at surfaces compared to the bulk. For the first time, we reveal a detailed atomistic picture of long-time scale restructuring of Pd deposited on Ag using microscopy, spectroscopy, and novel simulation methods. By developing and performing accelerated machine-learning molecular dynamics followed by an automated analysis method, we discover and characterize previously unidentified surface restructuring mechanisms in an unbiased fashion, including Pd-Ag place exchange and Ag pop-out as well as step ascent and descent. Remarkably, layer-by-layer dissolution of Pd into Ag is always preceded by an encapsulation of Pd islands by Ag, resulting in a significant migration of Ag out of the surface and a formation of extensive vacancy pits within a period of microseconds. These metastable structures are of vital catalytic importance, as Ag-encapsulated Pd remains much more accessible to reactants than bulk-dissolved Pd. We report our approach is broadly applicable to complex multimetallic systems and enables the previously intractable mechanistic investigation of restructuring dynamics at atomic resolution.

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

Attosecond science based on high harmonic generation from gases and solids

Recent progress in high power ultrafast short-wave and mid-wave infrared lasers has enabled gas-phase high harmonic generation (HHG) in the water window and beyond, as well as the demonstration of HHG in condensed matter. In this Perspective, we discuss the recent advancements and future trends in generating and characterizing soft X-ray pulses from gas-phase HHG and extreme ultraviolet (XUV) pulses from solid-state HHG. Then, we discuss their current and potential usage in time-resolved study of electron and nuclear dynamics in atomic, molecular and condensed matters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗