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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Evaluation of Methods for Molecular Dynamics Simulation of Ionic Liquid Electric Double Layers

We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.

Haskins, Justin B.↗

NASA NTRS2001

Carbon Nanotubes in Water: MD Simulations of Internal and External Flow, Self Organization

We have developed computational tools, based on particle codes, for molecular dynamics (MD) simulation of carbon nanotubes (CNT) in aqueous environments. The interaction of CNTs with water is envisioned as a prototype for the design of engineering nano-devices, such as artificial sterocillia and molecular biosensors. Large scale simulations involving thousands of water molecules are possible due to our efficient parallel MD code that takes long range electrostatic interactions into account. Since CNTs can be considered as rolled up sheets of graphite, we expect the CNT-water interaction to be similar to the interaction of graphite with water. However, there are fundamental differences between considering graphite and CNTs, since the curvature of CNTs affects their chemical activity and also since capillary effects play an important role for both dynamic and static behaviour of materials inside CNTs. In recent studies Gordillo and Marti described the hydrogen bond structure as well as time dependent properties of water confined in CNTs. We are presenting results from the development of force fields describing the interaction of CNTs and water based on ab-initio quantum mechanical calculations. Furthermore, our results include both water flows external to CNTs and the behaviour of water nanodroplets inside heated CNTs. In the first case (external flows) the hydrophobic behaviour of CNTs is quantified and we analyze structural properties of water in the vicinity of CNTs with diagnostics such as hydrogen bond distribution, water dipole orientation and radial distribution functions. The presence of water leads to attractive forces between CNTs as a result of their hydrophobicity. Through extensive simulations we quantify these attractive forces in terms of the number and separation of the CNT. Results of our simulations involving arrays of CNTs indicate that these exhibit a hydrophobic behaviour that leads to self-organising structures capable of trapping water clusters. In the second case (internal flows) we study the behaviour of water droplets confined inside CNTs. Constant temperature simulations allow us to capture structural properties such as the contact angles and density profiles of the equilibrated drops. By heating and subsequently cooling of the CNT, we are able to measure the evaporation and the condensation rate of the entrapped water.

Jaffe, Richard L.↗

DOE OSTI2025

Energy-efficient, Large-scale Molecular Dynamics Simulations via Hardware- and Algorithm-level Optimization

This work aims to develop a framework for energy-efficient computing that will enable molecular dynamics (MD) simulations of large-scale phenomena with atomic precision and simultaneously remove computational bottlenecks limiting the speed of MD simulations. We seek to implement such an approach through the development of surrogate models for the interatomic force calculation combined with the use of mixed numerical precision formats. For a model system of neutral atoms (only pairwise interactions), significant force calculation efficiency improvements were achieved, without detrimental effects on atomic structures or average energies, using single precision, by developing a surrogate model (deep neural network), and by quantizing this surrogate model. For a model system of charged atoms, the reciprocal-space calculation of electrostatic interactions was identified as the main bottleneck, and the development of a surrogate model should be pursued to achieve an estimated one-order-of-magnitude additional speedup.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2026

Valence-free open nanoparticle superlattices

A cornerstone of advanced materials design is establishing a framework for assembling nanoparticle superstructures with tailored symmetries. A longstanding challenge has been assembling diamond-like superstructures for photonic devices. Traditionally, such open superstructures require functionalized nanoparticles with directional or anisotropic interactions, reminiscent of valence bonding in a diamond. Here, we present a robust strategy for assembling valence-free nanoparticles into a broad array of cubic superstructures. By grafting nanoparticles with oppositely charged, end-functionalized water-soluble polymers of adjustable molecular weight, we gain control over electrostatic interactions and conformational constraints. This unified approach yields lattices analogous to rock salt, CsCl, zinc-blende, diamond, and the rare simple cubic phase, with tunable lattice constants. Theoretical models and simulations elucidate the underlying interactions, providing a framework for engineering valence-free nanoparticle superlattices.

36 MATERIALS SCIENCE↗

NASA NTRS2001

Electrostatic Charging of Polymers by Particle Impact at Low Pressures

Studies of the electrostatic interaction between micrometer-sized particles and polymer surfaces are of great interest to NASA's planetary exploration program. The unmanned landing missions to Mars planned for this decade as well as the possible manned missions that might take place during the second decade of this century require a better understanding of the electrostatic response of the materials used in landing crafts and equipment when exposed to wind-blown dust or to surface dust and sand particles. We report on preliminary experiments designed to measure the electrostatic charge developed on five polymer surfaces as they are impacted simultaneously by Mars simulant particles less than 5 micrometers in diameter moving at 20 m/s. Experiments were performed in a CO2 atmosphere at 10 mbars of pressure using a particle delivery method that propels the particles with contact. Experiments were also performed in dry air at atmospheric pressures using a pressurized particle delivery system. The five polymer surfaces, commonly used in space applications, were chosen so that they span the triboelectric series.

Calle, Carlos I.↗

DOE OSTI2026

Restoring Symmetry and Enhancing Exchange via Chiral Molecular-Magnetic Hexagons

Modified behaviors of molecules, designed for device or energy applications, can occur due to lattice-molecule incompatibilities in point-group symmetries (PGS), long-range interactions between neighboring molecules, or latticemolecule charge transfer. Such sensitivities are prerequisite for the use of quantum molecules as devices, but removing certain broken symmetries is desirable from the standpoint of simplifying their behaviors. Here, we demonstrate symmetry restoration to molecular-magnetic lattices and show that it leads to strengthened exchange-coupling due to dipolar electrostatic interactions between neighboring molecules and due to manifestations of an earlier identification of a spin-sensitive slightly mobile electron by Hooshmand et al. Starting with a broken symmetry molecular magnet of current interest, we demonstrate a more robust structure composed of six molecular magnets on a h-BN-like surface. The broken point-group symmetry is healed by constructing interpenetrating equilateral triangular structures that are appropriately ratcheted to preserve three-fold rotational symmetry and inversion symmetry.

Pederson, Mark R.↗

NASA NTRS1972

Elastic properties of alpha quartz and the alkali halides based on an interatomic force model.

A two-body central-force atomic model can be used to describe accurately the elastic properties of alpha quartz if the nontetrahedral O:O forces are included. The strength of the Si:O interaction has little effect on the bulk modulus. The technique is sufficiently general to allow calculations of the elastic properties of a specified structure under arbitrary pressure from a complete description of the interatomic forces. The elastic constants for the NaCl structure and the CsCl structure are examined. Our model includes two-body, central, anion-anion, anion-cation, and electrostatic interactions.

Weidner, D. J.↗

NASA NTRS1999

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗

NASA NTRS2013

Physical and Chemical Toeholds for Exoplanet Bioastronomy

If a search for exoplanet life were mounted today, the likely focus would be to detect oxygen (or ozone) in the atmosphere of a water-bearing rocky planet orbiting roughly 1AU from a G-type star. This appropriately conservative and practical default is necessary in large part because biological input on the question of where and how to look for life has progressed little beyond a purely empirical reliance on the example of terrestrial biology. However, fundamental physical and chemical considerations may impose significant yet universal constraints on biological potential. The liquid water + oxygen paradigm will be considered as an example, with a focus on the question, is liquid water a prerequisite for life? . Life requires a solvent to mediate interactions among biological molecules. A key class of these interactions is molecular recognition with high specificity, which is essential for high fidelity catalysis and (especially) information processing. For example, to correctly reproduce a string consisting of 600,000 units of information (e.g., 600 kilobases, equivalent to the genome of the smallest free living terrestrial organisms) with a 90% success rate requires specificity greater than 10(exp 7):1 for the target molecule vs. incorrect alternatives. Such specificity requires (i) that the correct molecular association is energetically stabilized by at least 40 kJ/mol relative to alternatives, and (ii) that the system is able to sample among possible states (alternative molecular associations) rapidly enough to allow the system to fall under thermodynamic control and express the energetic stabilization. We argue that electrostatic interactions are required to confer the necessary energetic stabilization vs. a large library of molecular alternatives, and that a solvent with polarity and dielectric properties comparable to water is required for the system to sample among possible states and express thermodynamic control. Electrostatic associations can be made in non-polar solvents, but the resulting complexes are too stable to be "unmade" with sufficient frequency to confer thermodynamic control on the system. Such considerations suggest that water, or a solvent with properties very like water, is necessary to support high-fidelity information processing a feature that must be common to all biology and can therefore be considered a critical prerequisite for life.

Hoehler, Tori↗

NASA NTRSDate not supplied

Design of Electrostatic Dust Lofting Suborbital Flight Experiment Examining Photoionization under Lunar Gravity

Dust on the lunar surface electrostatically charges due to the plasma environment surrounding the Moon, causing grains to become lofted and adhere to nearby surfaces including landers and astronauts. Studying the behaviors of these charged particles in the lunar environment is essential to plan around the deleterious effects of dust to future Moon missions. Models attempt to predict the amount of dust loading that can be expected in many of these scenarios, but they require experimental validation to be predictive. This physics cannot be fully studied on Earth due to the six times larger gravitational force obscuring the electrostatic interactions, so it is necessary to run experiments in a more relevant environment, including vacuum and near-lunar gravitational effects. An experiment has been designed to fly on the Lunar Gravity Acceleration (LGA) mission aboard the Blue Origin New Shepard suborbital rocket. This experiment will perform photoionization charging of lunar regolith simulant grains under the illumination of an ultraviolet (UV) source. As a result, the charged grains will then electrostatically repel one another and loft in the reduced gravity environment; their trajectories will be imaged via a high-speed camera. Preliminary laboratory results influencing the design of this experiment will be presented, including characterization of several UV sources, measurements of photoionization currents under various vacuum conditions, and examination of lunar simulant dust lofting under terrestrial gravity. Results from this flight will be compared with ground-based testing and the laboratory results outlined above to examine the dependence on gravity and will be fed into the dust charging and lofting models currently under development.

electrostatics↗

DOE OSTI2026

Cosolvent-tuned interactions in ionic liquids: A vibrational and quantum-chemical study of ethylene glycol ratio effects

Ionic liquids (ILs) are attractive media for CO 2 capture but remain limited by viscosity and cost. Blending ILs with ethylene glycol (EG) is a practical route to mitigate these constraints, yet the molecular origins of cosolvent effects and their dependence on composition are not well resolved. We combine Fourier-transform infrared (FT-IR) spectroscopy with quantum-chemical (DFT) analysis to elucidate how the IL:EG molar ratio modulates intermolecular interactions and electronic structure. Computed vibrational frequencies enable mode assignment and deconvolution of overlapping bands, revealing systematic, ratio-dependent shifts and broadenings in (i) EG O–H stretching, (ii) cation and EG C–H stretchings (imidazolium C2–H, C4–H, C5–H, methyl and ethyl groups, -CH2 of EG), (iii) anion signature modes (e.g., CN motifs), and (iv) EG C–O and C–C stretchings, consistent with the redistribution of hydrogen-bonding networks. Molecular electrostatic potential (MESP) maps quantify attenuation of extreme potential regions with increasing EG, indicating progressive screening of cation–anion electrostatic interactions. Quantum Theory of Atoms in Molecules (QTAIM) identifies emergent bond critical points between EG and the IL ions, while Reduced Density Gradient–Noncovalent Interaction (RDG–NCI) analysis differentiates strong directional hydrogen bonds from dispersive contacts across compositions. Together, these results show that EG fraction controls a switch from predominantly ion–ion to mixed ion–EG coordination, altering local polarity and polarizability that underlie the observed FT-IR trends. The framework provides composition–structure–spectrum relationships that can guide rational selection of IL:EG ratios to balance favorable molecular interactions with practical performance targets in scalable CO 2 capture systems.

DAC↗

NASA NTRS2011

The Importance of Water for High Fidelity Information Processing and for Life

Is water an absolute prerequisite for life? Life depends on a variety of non-covalent interactions among molecules, the nature of which is determined as much by the solvent in which they occur as by the molecules themselves. Catalysis and information processing, two essential functions of life, require non-covalent molecular recognition with very high specificity. For example, to correctly reproduce a string consisting of 600,000 units of information (e.g ., 600 kilobases, equivalent to the genome of the smallest free living terrestrial organisms) with a 90% success rate requires specificity > 107 : 1 for the target molecule vs. incorrect alternatives. Such specificity requires (i) that the correct molecular association is energetically stabilized by at least 40 kJ/mol relative to alternatives, and (ii) that the system is able to sample among possible states (alternative molecular associations) rapidly enough to allow the system to fall under thermodynamic control and express the energetic stabilization. We argue that electrostatic interactions are required to confer the necessary energetic stabilization vs. a large library of molecular alternatives, and that a solvent with polarity and dielectric properties comparable to water is required for the system to sample among possible states and express thermodynamic control. Electrostatic associations can be made in non-polar solvents, but the resulting complexes are too stable to be "unmade" with sufficient frequency to confer thermodynamic control on the system. An electrostatic molecular complex representing 3 units of information (e.g., 3 base pairs) with specificity > 107 per unit has a stability in non-polar solvent comparable to that of a carbon-carbon bond at room temperature. These considerations suggest that water, or a solvent with properties very like water, is necessary to support high-fidelity information processing, and can therefore be considered a critical prerequisite for life.

Hoehler, Tori M.↗

NASA NTRS2001

Membrane rejection of nitrogen compounds

Rejection characteristics of nitrogen compounds were examined for reverse osmosis, nanofiltration, and low-pressure reverse osmosis membranes. The rejection of nitrogen compounds is explained by integrating experimental results with calculations using the extended Nernst-Planck model coupled with a steric hindrance model. The molecular weight and chemical structure of nitrogen compounds appear to be less important in determining rejection than electrostatic properties. The rejection is greatest when the Donnan potential exceeds 0.05 V or when the ratio of the solute radius to the pore radius is greater than 0.8. The transport of solute in the pore is dominated by diffusion, although convective transport is significant for organic nitrogen compounds. Electromigration contributes negligibly to the overall solute transport in the membrane. Urea, a small organic compound, has lower rejection than ionic compounds such as ammonium, nitrate, and nitrite, indicating the critical role of electrostatic interaction in rejection. This suggests that better treatment efficiency for organic nitrogen compounds can be obtained after ammonification of urea.

NASA Discipline Life Support Systems↗

DOE OSTI2025

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, ↗

DOE OSTI2025

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↗

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↗

DOE OSTI2024

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↗

NASA NTRS2005

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↗