Engineering PapersSearch

SEARCH · Engineering Papers

Results for “electrostatic interactions”

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 73 records · Page 4

A Universal Augmentation Framework for Long-Range Electrostatics in Machine Learning Interatomic Potentials

Most current machine learning interatomic potentials (MLIPs) rely on short-range approximations, without explicit treatment of long-range electrostatics. To address this, we recently developed the Latent Ewald Summation (LES) method, which infers electrostatic interactions, polarization, and Born effective charges (BECs), just by learning from energy and force training data. Here, in this study, we present LES as a standalone library, compatible with any short-range MLIP, and demonstrate its integration with methods such as MACE, NequIP, Allegro, CACE, CHGNet, and UMA. We benchmark LES-enhanced models on distinct systems, including bulk water, polar dipeptides, and gold dimer adsorption on defective substrates, and show that LES not only captures correct electrostatics but also improves accuracy. Additionally, we scale LES to large and chemically diverse data by training MACELES-OFF on the SPICE set containing molecules and clusters, making a universal MLIP with electrostatics for organic systems, including biomolecules. MACELES-OFF is more accurate than its short-range counterpart (MACE-OFF) trained on the same data set, predicts dipoles and BECs reliably, and has better descriptions of bulk liquids. By enabling efficient long-range electrostatics without directly training on electrical properties, LES paves the way for electrostatic foundation MLIPs.

Kim, Dongjin [University of California, Berkeley,

Probing specific ion effects at air-aqueous dibutyl phosphate interfaces using vibrational sum frequency generation spectroscopy

Molecular properties at air–liquid and liquid–liquid interface hold the key to many processes involving molecular transport across phase boundaries from aerosol formation to carbon cycling and material separation using solvent extraction techniques. Using dibutyl phosphate (DBP) as a representative for partially aqueous soluble surfactants, the specific ion effect (SIE) of the Hofmeister series cations Cs + , Na + , Li + , and Mg 2+ on the partition and interaction between surfactant molecules and water molecules in the air–aqueous interface are investigated using vibrational sum frequency generation spectroscopy and surface tension measurements. In the presence of 1 mM and 1M bulk aqueous phase ionic strength salt concentrations, fundamental qualitative relationships are observed for the salting out of DBP relative to bulk aqueous phase nitrate salt concentrations and the specific cations species. At 1 mM ionic strength, the interfacial charge and hence the interfacial potential modulates the electrostatic interactions; in particular, the counter cations partially screen the negatively charged interface induced by the DBP in a direct Hofmeister order. At 1M ionic strength, the electric field at the interface or interfacial potential is effectively neutralized, and the counter cations promote the partitioning of DBP to the interface depending on their specific interaction with the DBP head group and metal ion hydration properties. The present results lay a foundation to study SIEs of heavier metals on hydrophobic-aqueous DBP interfaces.

Louie, Christina [Washington State University, Pul

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

Predicting Partial Atomic Charges in Metal–Organic Frameworks: An Extension to Ionic MOFs

Molecular simulation is an invaluable tool to predict and understand the usage of metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges, commonly obtained from density functional theory (DFT) calculations, are often required to model the electrostatic interactions between the MOF and adsorbates, especially when the adsorbates have dipole or quadrupole moments, such as water and CO 2 . Machine learning (ML) models have been previously employed to predict partial charges and avoid the computational cost associated with DFT calculations. However, previous ML models suffer from small training data sets, which limit their scope of application. In this work, we introduce two novel machine learning models, PACMOF2-neutral and PACMOF2-ionic, aimed at predicting the density-derived electrostatic and chemical (DDEC6) partial atomic charges for both neutral and ionic MOFs. These models not only yield DFT-level accuracy at a fraction of the computational cost but also demonstrate a remarkable improvement in prediction of adsorption, as validated with grand canonical Monte Carlo simulations. Furthermore, the robustness and fast computational time of the PACMOF2 models, along with their transferability to other porous materials such as covalent organic frameworks and zeolites, underscores their potential in high-throughput screening of MOFs for diverse applications.

36 MATERIALS SCIENCE

Partial Model of Insulator/Insulator Contact Charging

Two papers present a two-phase equilibrium model that partly explains insulator/ insulator contact charging. In this model, a vapor of ions within a gas is in equilibrium with a submonolayer of ions of the same species that have been adsorbed on the surface of an insulator. The surface is modeled as having localized states, each with a certain energy of adsorption for an ion. In an earlier version of the model described in the first paper, the ions do not interact with each other. Using the grand canonical ensemble, the chemical potentials of both vapor and absorbed phases are derived and equated to determine the vapor pressure. If a charge is assigned to the vapor particles (in particular, if single ionization is assumed), then the surface charge density associated with adsorbed ions can be calculated as a function of pressure. In a later version of the model presented in the second paper, the submodel of the vapor phase is extended to include electrostatic interactions between vapor ions and adsorbed ones as well as the screening effect, at a given distance from the surface, of ions closer to the surface. Theoretical values of this model closely match preliminary experimental data on the discharge of insulators as a function of pressure.

Hogue, Michael

Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water

Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.

Adsorption

Interconnected nanoconfining pore networks enhance catalyst CO 2 interaction in electrified reactive capture

Systems that sequentially capture and upgrade CO 2 from air to fuels/fuel-intermediates, such as syngas and ethylene, rely on an energy-intensive CO 2 release process. Electrified reactive capture systems transform CO 2 obtained directly from carbonate capture liquid into products. Previous reactive capture systems show a decline in Faradaic efficiencies (FE) at current densities above 200 mA/cm 2 . Here we show the chemical origins of this problem, finding that prior electrocatalyst designs failed to arrest, activate, and reduce in situ-generated CO 2 (i-CO 2 ) before it traversed the catalyst layer and entered the tailgas stream. We develop a templated synthesis to define pore structures and the sites of Ni single atoms, and find that carbon-nitrogen-based nanopores are effective in accumulating i-CO 2 via short-range, non-electrostatic interactions between CO 2 molecules and the nanochannel walls. These interactions confine and enrich i-CO 2 within the pores, enhancing its binding and activation. We report as a result carbonate electrolysis at 300 mA/cm 2 with FE to CO of 50% ± 3%, and with <1% CO 2 in the tailgas outlet stream. This corresponds to a projected energy efficiency (EE) to 2:1 syngas of 46% at 300 mA/cm 2 when H 2 is added using a water electrolyzer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Shear-Induced Conformations of Salt-Free Polyelectrolytes in Semidilute Solutions

Here, we present coarse-grained molecular dynamics simulations of salt-free polyelectrolyte chains in semidilute solutions under simple shear flow, with full hydrodynamic interactions and explicit dipolar solvent. At equilibrium, chain orientation statistics follow a pseudo-Voigt distribution, and the structural correlation length and chain end-to-end vector autocorrelation function exhibit scaling behavior consistent with theoretical predictions for polyelectrolytes. Under shear, chains transition from coiled to stretched states and the end-to-end vector autocorrelation function reveals oscillatory dynamics at high Weissenberg numbers. Analysis of the gyration tensor and shear strain distributions identified three distinct chain populations with directional alignment along and against the shear gradient. Compared with their neutral polymer counterparts, polyelectrolytes exhibit stronger shear thinning and enhanced chain alignment under the same Weissenberg number, which is attributed to electrostatic interactions and shear-induced counterion release. These findings provide molecular insight into the distinct flow response of charged polymers and have implications for tailoring the rheological properties of polyelectrolyte-based materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Investigating Particle-Particle Electrostatic Effects on Charged Lunar Dust Transport via Discrete Element Modeling

NASA surface exploration missions have always seen negative effects of dust including the Apollo missions. The astronaut-witnessed unusual behavior of the dust particles that surround the vehicle after engine cutoff has the potential to have more of an influence on surface systems dust loading than the high velocity lunar rocket plume ejecta in the landing process. The levitation and transport of the fine components of regolith on lunar surface has been linked to electrostatic effects and electric field, but so far there is no accurate model considering the inter-particle electrostatic interactions, especially when the particles are charged by rocket plume or other mechanical interactions due to exploration activities. This study is proposed to investigate the dynamics of charged lunar regolith with a discrete element method (DEM) approach focusing on the inter-particle interactions and contact charge transfer. The grain dynamics is coupled with mechanical and electrical particle interactions, and both short- and long-range interactions between spherical particles are incorporated. A tribo-charging model based on instantaneous collisions between particles is adopted and validated by comparing the simulation results to existing experimental data. Sensitivity analysis is conducted to quantify the effects of initial charge, tribo-charging, and E-field on transport of lunar dust based on JSC-1 simulants with a radius of 50 lm. DEM simulations are also conducted in a near realistic lunar environment with the estimations of initial conditions that shows the difference in position and velocity distributions between charged particles and uncharged particles. The results indicate that the charged dust particles have higher dispersion of position and velocity by several orders of magnitude due to electrostatic effects. This provides a potential explanation for the phenomena of the approximately 30 s dust lofting following Apollo Lunar Module landing.

Lunar dust

pH profile of the adsorption of nucleotides onto montmorillonite. II - Adsorption and desorption of 5-prime-AMP in iron-calcium montmorillonite systems

The interaction of 5-prime-AMP with montmorillonite saturated with various ratios of two metals found ubiquitously on the surface of earth, that is, iron and calcium, is investigated. Adsorption and desorption of the nucleotide were studied in the pH range of 2-12 at three levels of addition: 0.080, 0.268 and 0.803 mmole 5-prime-AMP per gram of clay. Two desorption stages were employed - H2O wash and NaOH extraction (pH = 12.0). 5-prime-AMP was preferentially adsorbed on the Fe-containing clays relative to the Ca clay. The nucleotide was fully recovered by the two desorption stages, mostly by the NaOH extraction. The evidence at hand indicates that 5-prime-AMP reaction with clay is affected by electrostatic interactions involving both attraction and repulsion forces. Some specific adsorption, possibly the result of covalent bonding and complex formation with the adsorbed ion, cannot be ruled out for iron but does not appear to operate for calcium. Changes in pH cause varying degrees of attaction and repulsion of 5-prime-AMP and may have been operating on the primitive earth, leading to sequences of adsorption and release of this biomolecule.

Banin, A.

PowderSim: Lagrangian Discrete and Mesh-Free Continuum Simulation Code for Cohesive Soils

PowderSim is a calculation tool that combines a discrete-element method (DEM) module, including calibrated interparticle-interaction relationships, with a mesh-free, continuum, SPH (smoothed-particle hydrodynamics) based module that utilizes enhanced, calibrated, constitutive models capable of mimicking both large deformations and the flow behavior of regolith simulants and lunar regolith under conditions anticipated during in situ resource utilization (ISRU) operations. The major innovation introduced in PowderSim is to use a mesh-free method (SPH-based) with a calibrated and slightly modified critical-state soil mechanics constitutive model to extend the ability of the simulation tool to also address full-scale engineering systems in the continuum sense. The PowderSim software maintains the ability to address particle-scale problems, like size segregation, in selected regions with a traditional DEM module, which has improved contact physics and electrostatic interaction models.

Johnson, Scott

Cation selectivities in zwitterion grafted nanopores: effect of zwitterion architecture

Selective separation of monovalent cations is a critical challenge in applications such as water purification and lithium recovery from salt brines. Cross-linked zwitterionic amphiphilic copolymer (ZAC-X) membranes have gained attention for their exceptional anion permselectivity, attributed to self-assembled zwitterion-lined nanodomains that interact preferentially with anions according to their hydrated radii r hyd . However, these membranes show minimal selectivity among monovalent cations, despite significant differences in their hydration structures, motivating studies on the underlying mechanisms of cation transport and selectivity in this family of materials. In this study, we conducted molecular dynamics simulations of aqueous salt solutions within zwitterion-functionalized nanopores to elucidate the influence of dipole orientation of the zwitterionic (ZI) ligands on cation diffusivities, partitioning, and permeabilities. To this end, we examined two contrasting ZI ligand architectures: Motif A (surface–cation–anion, S–ZI + –ZI − ) and Motif B (surface–anion–cation, S–ZI − –ZI + ). Our results show that in Motif A, the sulfonate groups of the ZI ligands are localized near the pore center radially, leading to strong electrostatic interactions with small bare cations (Mg 2+ and Li + ). This configuration results in high cation partitioning but low cation diffusion, maintaining solution-diffusion tradeoff typical of functionalized membranes. In contrast, Motif B show that sulfonate groups shift radially toward the mid-region of the pore. This shift, especially for small bare cations, introduces steric constraints that weaken their interactions with the sulfonate groups, thereby enhancing hydration and lowering partitioning, while still maintaining their low self-diffusivity. These findings establish zwitterion dipole orientation as a powerful design lever for tuning cation selectivity in membrane systems and offer molecular-level insights for engineering next-generation ion separation materials.

Morishita, Kazuya [Univ. of Texas, Austin, TX (Uni

Modifying Specific Ion Effects: Studies of Monovalent Ion Interactions with Amines

Specific ion effects in the interactions of monovalent anions with amine groups-one of the hydrophilic moieties found in proteins-were investigated using octadecylamine monolayers floating at air–aqueous solution interfaces. Here, we find that at solution pH 5.7, larger monovalent anions induce a nonzero pressure starting at higher areas/molecules, i.e., a wider “liquid expanded” region in the monolayer isotherms. Using X-ray fluorescence at near total reflection (XFNTR), an element- and surface-specific technique, ion adsorption to the amines at pH 5.7 is confirmed to be ion-specific and to follow the conventional Hofmeister series. However, at pH 4, this ion specificity is no longer observed. We propose that at the higher pH, the amine headgroups are only partially protonated, and large polarizable ions such as iodine are better able to boost amine protonation. At the lower pH, on the other hand, the monolayer is fully protonated, and electrostatic interactions dominate over ion specificity. These results demonstrate that ion specificity can be modified by changing the experimental conditions.

36 MATERIALS SCIENCE

Charge regulation effects on colloidal mixture nanoparticles

Changes in pH within a system containing dissociable sites affect the protonation and deprotonation of these groups, thereby influencing their physical properties. In response, the system modifies their surface charge, affecting electrostatic interactions, aggregation, stability, and structural behavior. Although the pH can be tuned in experiments, it is difficult to model this phenomenon using simulations or theoretical approaches. Here, we perform hybrid Monte Carlo-molecular dynamics simulations to model charge regulation effects in an equimolar colloidal charged system. We compare charge regulation effects with those of a system in which the charges of colloidal nanoparticles are not dissociable. The comparison between the two cases modifies the phase diagram, and it changes the volume fraction where a percolation network of nanoparticles is found. Charge regulation is found to destroy network formation, as the charge in the nanoparticles is modified because of the cooperativity dependency of the degree of charge dissociation sites among the nanoparticles favoring cluster formation. Furthermore, our work suggests that the ionic and/or electronic conductivity in functionalized nanoparticles can be modified by changing pH values. It also guides the experimental design of oppositely charged nanoparticles as inks for 3D printing processes.

Classical statistical mechanics

Stabilizing Nickel‐Rich Cathodes in Aqueous Process through Nanocellulose as Water Barrier

Nickel-rich LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NMC 811) cathode offers high voltage and high specific capacity, making it promising for high energy density batteries. However, large-scale manufacturing of aqueous-processed NMC 811 electrodes remains challenging due to proton exchange causing material decomposition and capacity loss. This work addresses this issue by constructing an in situ nanocellulose protective layer for NMC 811 particles via electrostatic interactions during the slurry preparation. For the first time, the interatomic spacing between inter-chains of nanocellulose is measured through wide-angle X-ray scattering and demonstrate the ability to effectively confine interlayer water using atomistic simulations. Moreover, this nanocellulose coverage simultaneously minimizes Li + surface segregation and mitigates water infiltration. Owing to less material decomposition during the aqueous processing, nanocellulose-protected NMC electrodes exhibit higher initial coulombic efficiency (83% vs 62% at 0.1C) and capacity (133 vs 59 mAh g −1 at 6C) than unprotected electrodes. Additionally, optimized aqueous-processed NMC electrodes offer comparable or even superior electrochemical properties compared to the electrodes fabricated using the conventional toxic organic solvent, N-methyl-2-pyrrolidone. Consequently, the developed approach enables affordable, sustainable aqueous processing for Nickel-rich NMC 811 cathodes with excellent electrochemical performances.

25 ENERGY STORAGE

Self‐Assembled Membranes for High Ion Selectivity and Proton Blocking in Electrochemical Applications

Anion-exchange membranes (AEMs) with high anion/cation selectivity and exceptional proton-blocking ability are critical for applications such as bipolar membrane electrodialysis and electrochemical acid recovery. However, existing AEMs are constrained by a trade-off between ionic conductivity and selectivity, largely due to the intrinsic coupling between charge density and water content, and they suffer from excessive proton leakage facilitated by the Grotthuss hopping mechanism. In this work, poly(vinylimidazolium) membranes functionalized with long alkyl side chains that self-assemble into well-defined microphase-separated morphologies stabilized by hydrophobic and electrostatic interactions are reported. These unique structures localize the charge density along the polymer backbone to promote fast and selective ion transport. As a result, these membranes exhibit ionic conductivities and counter-ion diffusivities surpassing those of conventional homogeneous membranes, along with unprecedented counter-ion/co-ion selectivity and proton-blocking ability. These results establish a new design paradigm for high-performance, phase separated charged polymer membranes that overcome the limitations of homogeneous membranes, with broad implications for advanced electrochemical technologies.

36 MATERIALS SCIENCE

The Anisotropic Adsorption of De Novo Allosteric Two‐Component Protein Fibers on Mica Surfaces

Protein adsorption at solid–liquid interfaces underlies many biomedical and materials applications, yet the mechanisms governing adsorption of proteins and their assemblies remain poorly understood. Here we investigate de novo–designed proteins that self-assemble into three fibrillar morphologies — small (S), large (L), and helical (H)—on muscovite mica. S-fibers are metastable, forming first but diminishing as L- and H-fibers develop and deposit. Adsorption of both fibers and monomers depends on fiber morphology and solvent environment, which is modulated by the substrate. The anisotropic surface features of the fiber types—long axis of S- and L-fibers and helical grooves of H-fibers—correlate with specific crystallographic directions on mica (001). S- and L-fibers align along the unique lattice axis, while the substrate-facing groove of H-fibers preferentially aligns along the remaining symmetry-related directions. Increasing potassium chloride (KCl) concentrations to molar levels alters adsorption, decreasing monomer coverage relative to fibers. These observations, interpreted through the effect of interfacial interface, indicate that ion-mediated modulation of hydration layers governs electrostatic interactions and alignment. This study reveals how coupling between protein topography, substrate crystallography, and hydration structure dictates selective adsorption and orientation of protein assemblies, offering insights for rational biomolecular material design.

AFM

Electrostatically Enhanced Buried Interface Binding of Self‐Assembled Monolayers for Efficient And Stable Inverted Perovskite Solar Cells

Inverted p‐i‐n structure perovskite solar cells (PSCs) have outperformed traditional n‐i‐p PSCs in recent years. A key advancement is the use of self‐assembled monolayers (SAMs) as hole transport layers. One class of widely used SAMs is carbazole‐based phosphonic acids. However, it is found that these SAMs lack strong binding with transparent conducting oxides (TCO) and perovskite. The weak binding strength results in suboptimal interfacial adhesion of the buried interface, which limits the device's stability. Here, interfacial binding is enhanced by increasing the dipole moment that creates a strong interfacial electric field that enhances electrostatic interactions at the TCO/perovskite interface, while incorporating tailored functional groups in SAMs to improve chemical anchoring to TCO and binding to perovskite. Specifically, the donor‐acceptor SAM molecule 4‐(7‐(4‐(bis(4‐methoxyphenyl)amino)‐2,5‐difluorophenyl)benzo[c][1,2,5]thiadiazol‐4‐yl)benzoic acid (PAFTB) is employed, which features an enhanced dipole moment along with electron‐donating and electron‐withdrawing functional groups to optimize interfacial interactions. Compared to extensively used [2‐(9H‐carbazol‐9‐yl)ethyl]phosphonic acid (2PACz), PAFTB enhances total interfacial adhesion by 2.8 times, thereby improving the thermal stability of the layer. Using this approach, PSCs are demonstrated with a certified quasi‐steady‐state power conversion efficiency of 24.9% and maintain 80% of the initial efficiency after 900 h of maximum power point tracking at 85 °C.

Interfacial adhesion