Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electrostatics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Substituent and Heteroatom Effects on π–π Interactions: Evidence That Parallel-Displaced π-Stacking is Not Driven by Quadrupolar Electrostatics

Stacking interactions are a recurring motif in supramolecular chemistry and biochemistry, where a persistent theme is a preference for parallel-displaced aromatic rings rather than face-to-face π-stacking. This is usually explained in terms of quadrupole–quadrupole interactions between the arene moieties but that interpretation is inconsistent with accurate calculations, which reveal that the quadrupolar picture is qualitatively wrong. At typical π-stacking distances, quadrupolar electrostatics may differ in sign from an exact calculation based on charge densities of the interacting arenes. We apply symmetry-adapted perturbation theory to dimers composed of substituted benzene and various aromatic heterocycles, which display a wide range of electrostatic interactions, and we investigate the interplay of Pauli repulsion, dispersion, and electrostatics as it pertains to parallel-displaced π-stacking. Profiles of energy components along cofacial slip-stacking coordinates support a prominent role for the “van der Waals model” (dispersion in competition with Pauli repulsion), even for polar monomers where electrostatic interactions are significant. While electrostatic interactions are necessary to explain the optimal face-to-face π-stacking distance and to account for the relative orientation of one polar arene with respect to another, we find no evidence to support continued invocation of quadrupolar electrostatics as a basis for π-stacking. Our results suggest that a driving force for offset-stacking exists even in the absence of electrostatic interactions. Consequently, tuning electrostatics via functionalization does not guarantee that slip-stacking can be avoided. This has implications for rational design of soft materials and other supramolecular architectures.

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

Colloidal Dispersions of Sterically and Electrostatically Stabilized PbS Quantum Dots: Structure Factors, Second Virial Coefficients, and Film-Forming Properties

Electrostatically stabilized nanocrystals (NCs) and, in particular, quantum dots (QDs) hold promise for forming strongly coupled superlattices due to their compact and electronically conductive surface ligands. However, studies of the colloidal dispersion and interparticle interactions of electrostatically stabilized sub-10 nm NCs have been limited, hindering the optimization of their colloidal stability and self-assembly. In this study, we employed small-angle X-ray scattering (SAXS) experiments to investigate the interparticle interactions and arrangement of PbS QDs with thiostannate ligands (PbS-Sn2S64-) in polar solvents. The study reveals significant deviations from the ideal solution behavior in electrostatically stabilized QD dispersions. Our results demonstrate that PbS-Sn2S64- QDs exhibit long-range interactions within the solvent, in contrast to the short-range steric repulsion characteristic of PbS QDs with oleate ligands (PbS-OA). Introducing highly charged multivalent electrolytes screens electrostatic interactions between charged QDs, reducing the length scale of the repulsive interactions. Furthermore, we calculated the second virial (B2) coefficients from SAXS data, providing insights into how surface chemistry, solvent, and size influence pair potentials. Finally, we explore the influence of long-range interparticle interactions of PbS-Sn2S64- QDs on the morphology of films produced by drying or spin-coating colloidal solutions. The long-range repulsive term of PbS-Sn2S64- QDs promotes the formation of amorphous films, and screening the electrostatic repulsion by the addition of an electrolyte enables the formation of crystalline domains. These findings highlight the critical role of NC-NC interactions in tailoring the properties of functional materials made of colloidal NCs.

SAXS↗
DOE OSTI2024

Charging, aggregation, and electrostatic dispersion of radioactive and nonradioactive particles in the atmosphere

Electrostatic dispersion can significantly impact the microphysical behavior of charged particles and ions until reaching zero space charge. However, although radioactive particles can be strongly charged in air, the influence of electrostatic dispersion has been neglected in understanding their behavior. This study is aimed at investigating the time evolution of the charge and size distributions of radioactive and nonradioactive particles in air and developing simple approaches for applications. With processes involving charging, aggregation, and electrostatic dispersion, a comprehensive population balance model (PBM) has been developed to examine particle charge/size distribution dynamics. It is shown that compared to nonradioactive particles, the charge and size distributions of radioactive particles may evolve differently with time because radioactivity and electrostatic dispersion can significantly affect the charging and aggregation kinetics of the particles. It is found that, after the Fukushima accident, background aerosols in the pathway of radioactive plumes might be highly charged due to ionizing radiation, suggesting that radiation fields may strongly influence in situ measurements of charged atmospheric particles. The comprehensive PBM is simplified, and then the verification and application of the simplified PBMs are discussed. This study provides useful insight into how radioactivity can affect the dynamic behavior of particles in atmospheric systems including radiation sources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Diabatization with Electrostatic Embedding for Studying Photophysics in Organic Molecular Crystals

Highly emissive organic molecular crystals find applications in several areas, such as organic electronics, solar cells, and sensors. Understanding the excited-state mechanisms underlying these applications is essential for optimizing and controlling them effectively. Exciton models coupled with nonadiabatic dynamics, particularly quantum dynamics, provide crucial insights into photochemical and photophysical processes in molecular crystals. Nevertheless, there remains a lack of general tools and automated workflows to facilitate such simulations. In this paper, we present a computational strategy to investigate the photoactivated dynamics of organic molecular crystals, bridging methodologies traditionally used for molecular systems and materials science, with a particular focus on the interplay between local excitations and charge transfer (CT) processes. We have implemented an interface between the fromage and Overdia programs, enabling the construction of vibronic Hamiltonians for molecular crystals within an excited-state ONIOM(QM:QM′) framework, incorporating long-range electrostatics through a RESP-based Ewald summation. Fragment-based diabatization provides a route to quantum dynamics simulations in weak-to-intermediate coupling regimes. The method was applied to the photophysics of dibenzo[g,p]chrysene (DBC) crystals using time-dependent DFT. The fromage/ Overdia interface was employed to compute the couplings of local excitations and CT states for 18 unique DBC dimers in the crystal and to quantify the influence of electrostatic embedding, which was found to be modest (10−20%). Simulations on π-stacked dimers reproduced the small red shift observed experimentally from solution to crystal, attributed to electronic interactions among fixed monomers rather than crystal electrostatics. Quantum dynamics simulations revealed ultrafast population transfer from bright local excitations to CT states. This approach establishes a robust framework linking molecular and solid-state excited-state dynamics, with potential applications for studying excitations, defects, and impurities in molecular crystals.

Crystals↗
DOE OSTI2026

Electrostatically driven pattern formation in mixed charged–neutral multicomponent elastic membranes

Multicomponent crystalline and amorphous elastic shells exhibit heterogeneous surface patterns that provide distinctive functionalities in cellular environments. Such patterning typically arises from the competition between short-range attractive and long-range repulsive interactions in membranes. Here, we demonstrate that the intrinsic competition between electrostatic repulsion and elastic deformation is sufficient to drive spontaneous surface patterning in elastic shells, requiring no additional attractive interactions. Using numerical simulations, we demonstrate pattern formation in mechanically homogeneous membranes with heterogeneous surface charge composition across different topologies, including spheres, discs, and flat periodic membranes. We also examine patterns in crystalline and amorphous shells of coassembled charged and neutral components with different bending rigidities. At low charge fraction, discrete charged surface domains form. At intermediate charge fraction, the competition between electrostatics and elasticity leads to elongated domains (rods) of the charged component, which results in lamellar patterns at nearly equal fraction of the charged and neutral components. At high charge fraction, nanodomains of the neutral component form. Amorphous shells exhibit similar progressions but with disordered structures rather than ordered lamellar patterns. These pattern morphologies are observed in both the closed shells and flat membranes. As salt concentration increases, all patterns coarsen due to the screening of electrostatic interactions.

charge heterogeneity↗
DOE OSTI2025

Two-Dimensional Binary Superlattice of BNNT-Surfactant Vesicle Complex Induced by Electrostatic Interaction

For a wide range of practical applications of boron nitride nanotubes (BNNTs), it is essential to achieve their highly ordered self-assembled structures. This study reports on a two-dimensional (2D) binary superlattice of individually exfoliated BNNTs with a negative surface charge (p-BNNT25) and cationic surfactant vesicles (CTAT/SDBS vesicles, prepared by mixing cetyltrimethylammonium tosylate (CTAT) and sodium dodecylbenzenesulfonate (SDBS)) complexes through electrostatic interactions. Depending on the surface charge density of the CTAT/SDBS vesicles and the mass ratio between the CTAT/SDBS vesicle and p-BNNT25, the CTAT/SDBS-BNNT complexes formed highly ordered superstructures. These structures include an intercalated lamellar phase with a centered rectangular structure (ICLP), in which a 2D array of p-BNNT25 is inserted into the multilamellar structure, and an AB 3 structure, in which the BNNTs are surrounded by surfactant micelles in a triangular arrangement. To the best of our knowledge, this is the first demonstration of the fabrication of highly ordered superstructures of individually exfoliated and negatively charged BNNTs with positively charged surfactant vesicles through electrostatic interactions. This approach for the 2D binary superlattices of CTAT/SDBS-BNNT complexes induced by electrostatic interactions is expected to be beneficial for a wide range of one-dimensional (1D) nanoparticle applications.

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

The Interplay of Pauli Repulsion, Electrostatics, and Field Inhomogeneity for Blueshifting and Redshifting Vibrational Probe Molecules

Many molecules' vibrational frequencies are sensitive to intermolecular electric fields, enabling them to probe the field in complex molecular environments. However, it is often unclear whether the probe is responding to the local electric field or other types of intermolecular interactions, inhibiting interpretation of the frequency and effectiveness as probes. This is especially true for molecules whose vibrational frequencies blueshift instead of the more typical redshift in hydrogen bonding configurations. In this work, we computationally investigate the causes of redshifting versus blueshifting over a range of vibrational reporters. First, we apply adiabatic energy decomposition analysis to a paradigmatic set of probes, finding that redshifting only occurs when electrostatic interactions are strong enough to overcome the dominant and large blueshifting contribution of Pauli repulsion. Furthermore, we demonstrate that field inhomogeneity can further shift the frequency of many probes substantially to either reinforce or counteract the shift expected from a homogeneous field. We find that redshifting is reinforced by electric field inhomogeneity, otherwise field inhomogeneity further weakens the electrostatic contribution relative to Pauli repulsion, leading to blueshifting. Further calculations indicate that the probe's response to field inhomogeneity can be understood by considering the mass of the atoms involved in the stretching mode and sign of the electric field. In explaining the interplay of different intermolecular interactions and field inhomogeneity for many probes, our results should enable the use and interpretation of spectroscopic probes and their connection to electric fields in more complex systems.

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

Direct Imaging of Asymmetric Interfaces and Electrostatic Potentials inside a Hafnia–Zirconia Ferroelectric Nanocapacitor

In hafnia-based thin-film ferroelectric devices, chemical phenomena during growth and processing, such as oxygen vacancy formation and interfacial reactions, appear to strongly affect device performance. However, the correlation between the structure, chemistry, and electrical potentials at the nanoscale in these devices is not fully known, making it difficult to understand their influence on device properties. Here, we directly image the composition and electrostatic potential with nanometer resolution in the cross section of a nanocrystalline W/Hf 0.5 Zr 0.5 O 2−δ (HZO)/W ferroelectric capacitor using multimodal electron microscopy. This reveals a 1.4 nm wide tungsten suboxide interfacial layer formed at the bottom interface during fabrication, which introduces a potential dip and leads to asymmetric switching fields. Additionally, we compare the measured potentials to DFT calculations and find it is nearly 3 V lower than expected in the HZO, which appears to be caused by oxygen vacancies and a resulting negative built-in potential. In conclusion, these chemical and electrostatic details are important to characterize and tune to achieve high-performance ferroelectric devices.

Defects in solids↗
DOE OSTI2025

A self-consistent electrostatic electrified shallow water model for charged liquid surface dynamics

The dynamics of an electrified liquid surface are investigated using a shallow water model that is self-consistently coupled with an electrostatic solver. To account for the spatiotemporal variation of a curved liquid surface, the electric field on curved surfaces is calculated by solving a two-dimensional electrostatic equation using the weighted least squares (WLSQ) interpolation method. First, the WLSQ implementation is verified with analytical theory obtained from a Laplace solution assuming a half-cylinder liquid surface profile. Then, the coupled electrified shallow water model is used to study the instability of liquid surface perturbation as a function of the potential drop between an electrode and conducting liquid, surface tension, and gravity. We present a linear dispersion theory of liquid surface instability for long-wavelength perturbations, including the effects of liquid viscosity. Furthermore, the effects of multiple sinusoidal surface perturbations on the electrified liquid surface instability are investigated.

Capillary waves↗
DOE OSTI2025

The influence of protein electrostatics on potential inversion in flavoproteins

Biology uses relatively few electron-transfer cofactors, tuning their potentials, electronic couplings, and reorganization energies to carry out the required chemistry. It is remarkable that the potential ordering of two-electron transfer active flavins can be normal (first oxidation at low potential and second oxidation at high potential) or inverted, and the gap between the potentials can be as large as one volt. Analysis based on structural bioinformatics and electrostatics indicates that the ordering of the flavin redox potential is influenced by protein electrostatics. In all 36 flavoproteins examined, the introduction of a negative charge near the flavin in silico increases the extent of potential inversion (by lowering the electrochemical potential of the second electron-transfer step); the introduction of a positive charge near the flavin favors normally ordered potentials. We also find that the addition of positive charges increases the electrochemical potential for the naturally occurring one-electron transition in flavodoxins (between deprotonated hydroquinone and neutral semiquinone) and also increases the second one-electron transition in bifurcating flavins (between anionic semiquinone and fully oxidized flavin). Finally, we find that proximity of a proton acceptor, notably conserved arginine, supports proton-coupled electron transfer because it may act as a proton acceptor, promoting potential inversion. This key arginine residue may enable two-electron transfer chemistry by promoting the proton-coupled electron transfer process over the pure electron transfer process, suggesting how a protein's flavin environment may influence one- or two-electron chemistry in flavoproteins.

Singh, Niven [Duke Univ., Durham, NC (United State↗
DOE OSTI2024

Drift kinetic electrostatic simulations of the edge localized mode heat pulse

In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Understanding the selectivity of nonsteroidal anti-inflammatory drugs for cyclooxygenases using quantum crystallography and electrostatic interaction energy

Quantum crystallography methods have been employed to investigate complex formation between nonsteroidal anti-inflammatory drugs (NSAIDs) and cyclooxygenase (COX) enzymes, with particular focus on the COX-1 and COX-2 isoforms. This study analyzed the electrostatic interaction energies of selected NSAIDs (flurbiprofen, ibuprofen, meloxicam and celecoxib) with the active sites of COX-1 and COX-2, revealing significant differences in binding profiles. Flurbiprofen exhibited the strongest interactions with both COX-1 and COX-2, indicating its potent binding affinity. Celecoxib and meloxicam showed a preference for COX-2, consistent with their known selectivity for this isoform, while ibuprofen showed comparable interaction energies with both isoforms, reflecting its nonselective inhibition pattern. Key amino-acid residues, including Arg120, Arg/His513 and Tyr355, were identified as critical determinants of NSAID selectivity and binding affinity. The findings highlight the complex interplay between interaction energy and selectivity, suggesting that while electrostatic interactions play a fundamental role, additional factors such as enzyme dynamics and the hydrophobic effect also contribute to the therapeutic efficacy and safety profiles of NSAIDs. These insights provide valuable guidance for the rational design of NSAIDs with enhanced therapeutic benefits and minimized adverse effects.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Local reduced-order modeling for electrostatic plasmas by physics-informed solution manifold decomposition

Despite advancements in high-performance computing and modern numerical algorithms, computational cost remains prohibitive for multi-query kinetic plasma simulations. Here, in this work, we develop data-driven reduced-order models (ROMs) for collisionless electrostatic plasma dynamics, based on the kinetic Vlasov-Poisson equation. Our ROM approach projects the equation onto a linear subspace defined by the proper orthogonal decomposition (POD) modes. We introduce an efficient tensorial method to update the nonlinear term using a precomputed third-order tensor. We capture multiscale behavior with a minimal number of POD modes by decomposing the solution manifold into multiple time windows and creating temporally local ROMs. We consider two strategies for decomposition: one based on the physical time and the other based on the electric field energy. Applied to the 1D1V Vlasov–Poisson simulations, that is, prescribed E-field, Landau damping, and two-stream instability, we demonstrate that our ROMs accurately capture the total energy of the system both for parametric and time extrapolation cases. The temporally local ROMs are more efficient and accurate than the single ROM. In addition, in the two-stream instability case, we show that the energy-windowing reduced-order model (EW-ROM) is more efficient and accurate than the time-windowing reduced-order model (TW-ROM). With the tensorial approach, EW-ROM solves the equation approximately 90 times faster than Eulerian simulations while maintaining a maximum relative error of 7.5% for the training data and 11% for the testing data.

Electrostatic plasmas↗
DOE OSTI2025

Scaling of energy delivered through an electrostatic discharge to a small series load

We study the energy delivered through a small-resistance series “victim” load during electrostatic discharge events in air. For gap lengths over 1 mm, the fraction of the stored energy delivered is mostly gap-length independent, with a slight decrease at larger gaps due to electrode geometry. The energy to the victim scales linearly with circuit capacitance and victim load resistance but is not strongly dependent on circuit inductance. This scaling leads to a simple approach to predicting the maximum energy that will be delivered to a series resistance for the case where the victim load resistance is lower than the spark resistance.

42 ENGINEERING↗
DOE OSTI2025

Ultra-sensitive radon assay using an electrostatic chamber in a recirculating system

Rare event searches such as neutrinoless double beta decay and Weakly Interacting Massive Particle detection require ultra-low background detectors. Radon contamination is a significant challenge for these experiments, which employ highly sensitive radon assay techniques to identify and select low-emission materials. This work presents the development of ultra-sensitive electrostatic chamber (ESC) instruments designed to measure radon emanation in a recirculating gas loop, for future lower background experiments. Unlike traditional methods that separate emanation and detection steps, this system allows continuous radon transport and detection. This is made possible with a custom-built recirculation pump. A Python-based analysis framework, PyDAn, was developed to process and fit time-dependent radon decay data. Radon emanation rates are given for various materials measured with this instrument. A radon source of known activity provides an absolute calibration, enabling statistically-limited minimal detectable activities of 20 µBq. These devices are powerful tools for screening materials in the development of low-background particle physics experiments.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Exploring Local Vibrational Structure in Protein-Bound Chlorophyll a: Isotope-Enrichment Experiments and Electrostatic Analysis

Local chlorophyll (Chl) vibrations play an essential role in biological photosynthesis by facilitating rapid energy transfer between pigments. In research studies, they also provide a useful spectroscopic probe of the local protein environment that surrounds each pigment. However, measuring the complete vibrational absorption spectrum of a protein-bound Chl molecule is much more difficult than, for example, Chl in neat solvent due to overlap with protein vibrations that typically drown out Chl vibrational features. Resonance Raman and fluorescence spectroscopies provide a way around this problem for Franck–Condon active vibrations, but these often rely on cryogenic measurement conditions and fail to capture vibrational signatures from, e.g., ester-group vibrations that lack coupling to a convenient electronic transition. In the present contribution, we use 13 C-enrichment of the protein backbone to shift protein background signals to lower frequency, providing a largely clean spectral window in which to study local Chl a C═O stretch modes. A room-temperature absorption spectrum for Chl a in the water-soluble chlorophyll protein (WSCP) of Lepidium virginicum is thus extracted as a difference between protein-plus-pigment and protein-only vibrational spectra. Excellent agreement in the molecular fingerprint region with the vibrational spectrum of Chl a in organic solvents confirms that the resulting spectrum represents the response of the protein-bound Chl a molecule. Furthermore, the ester group resonance is observed to shift in response to the S53P mutation that eliminates a 17 3 ester group hydrogen bond from the protein environment. Finally, we analyze these experimental results using MD-based electrostatic analysis, finding that electric-field mapping at the C atom of the ester group provides a satisfactory explanation for the observed frequency shifts between organic solvent and protein environment. MD analysis further suggests that a red-shifted ester peak observed experimentally for the S53P mutant results from solvation of the 17 3 Chl ester group due increase water penetration into the Chl-binding pocket relative to wild-type WSCP. Here, we anticipate that these results will prove useful both for benchmarking future simulation work and as a reference for interpreting Chl vibrational spectra as a probe of pigment–protein interactions.

14 SOLAR ENERGY↗
DOE OSTI2025

Accurate and Efficient Prediction of p K w in Aqueous Electrolytes Using Local Electrostatic Potentials

The K w =10 -pK w in aqueous electrolytes exhibit significant variations as a function of ion concentration. This behavior is not well-described by the Pitzer model, particularly at high concentrations. Here, we provide a molecular interpretation of the concentration dependent trends in p K w and develop a new method that accurately predicts the dissociation constant based upon the local electrostatic potential of the water oxygen that is determined by its nearest neighbors. The method is computationally efficient, relying only upon sampling of the local configurations (via molecular dynamics) and a small set of experimentally measured p K w . It is extensible to a range of salt compositions and the molecular understanding of the local potentials provides new routes toward the development of electrolytes with tailored physicochemical properties like p K w .

anions↗