Search NASASearch

SEARCH · Search NASA

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 19 records

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

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

Inferring effective electrostatic interaction of charge-stabilized colloids from scattering using deep learning

In this article, an innovative strategy is presented that incorporates deep auto-encoder networks into a least-squares fitting framework to address the potential inversion problem in small-angle scattering. To evaluate the performance of the proposed approach, a detailed case study focusing on charged colloidal suspensions was carried out. The results clearly indicate that a deep learning solution offers a reliable and quantitative method for studying molecular interactions. The approach surpasses existing deterministic approaches with respect to both numerical accuracy and computational efficiency. Overall, this work demonstrates the potential of deep learning techniques in tackling complex problems in soft-matter structures and beyond.

36 MATERIALS SCIENCE

pyEF: A Python Framework for QM and QM/MM Atom-Wise Electric Field Analysis

We introduce pyEF, a software package for computing molecular electric fields, electrostatic interaction energies, and electrostatic potentials from quantum mechanical (QM) atom-centered multipole expansions with atom-wise decomposable contributions. We demonstrate the computational efficiency and accuracy of this QM-derived electric field evaluation tool through several tests. To assess the influence of the underlying QM method and charge partitioning scheme on these electrostatic quantities, we analyze over 250 configurations of an acetone solute molecule in five solvents of variable polarity. We find that electric field calculations are highly sensitive to the choice of charge partitioning method. Even among real-space charge schemes, acetone Stark tuning rates differ by up to a factor of 2. Benchmarking computed solvent dipole moments against experimental bulk values, we conclude that the CM5, ADCH, and Hirshfeld-I charge schemes most reliably capture solvent electrostatics and therefore provide a more faithful foundation for computing electric fields. When constructed from these real-space charges, electric fields are nearly insensitive to basis set size and monotonically increase in magnitude with higher Fock exchange. We also demonstrate efficient convergence of QM electrostatics when more distant molecules are represented solely by MM point charges, reducing computational overhead. Leveraging these findings, we demonstrate the use of pyEF to deduce environmental effects on a transition metal complex from a Ga 4 L 6 12– nanocage and quantify the dominant role of organic linkers in orchestrating electrostatic preorganization.

electric fields

Single-Molecule Detection of the Encounter and Productive Electron Transfer Complexes of a Photosynthetic Reaction Center

Small, diffusible redox proteins play an essential role in electron transfer (ET) in respiration and photosynthesis, sustaining life on Earth by shuttling electrons between membrane-bound complexes via finely tuned and reversible interactions. Ensemble kinetic studies show transient ET complexes form in two distinct stages: an “encounter” complex largely mediated by electrostatic interactions, which subsequently, through subtle reorganization of the binding interface, forms a “productive” ET complex stabilized by additional hydrophobic interactions around the redox-active cofactors. Here, using single-molecule force spectroscopy (SMFS) we dissected the transient ET complexes formed between the photosynthetic reaction center-light harvesting complex 1 (RC-LH1) of Rhodobacter sphaeroides and its native electron donor cytochrome c 2 (cyt c 2 ). Importantly, SMFS resolves the distribution of interaction forces into low (~150 pN) and high (~330 pN) components, with the former more susceptible to salt concentration and to alteration of key charged residues on the RC. Thus, the low force component is suggested to reflect the contribution of electrostatic interactions in forming the initial encounter complex, whereas the high force component reflects the additional stabilization provided by hydrophobic interactions to the productive ET complex. Employing molecular dynamics simulations, we resolve five intermediate states that comprise the encounter, productive ET and leaving complexes, predicting a weak interaction between cyt c 2 and the LH1 ring near the RC-L subunit that could lie along the exit path for oxidized cyt c 2 . The multimodal nature of the interactions of ET complexes captured here may have wider implications for ET in all domains of life.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

High-throughput measurements of CO 2 permeance and solubility in ionic liquid reveal a synergistic role of ionic interactions and void fractions

The factors that govern CO 2 solubility in ionic liquids (ILs) are of great interest for the development of new materials for CO 2 capture and utilization. The cationic functional group (i.e., imidazolium, pyrrolidinium, pyridinium, etc.), alkyl chain length of cation, degree of fluorination of anion, anion size, and the void fraction in IL are known to influence CO 2 solubility. However, a comprehensive explanation of how these factors collectively affect CO 2 solubility has not been developed yet. This knowledge gap is largely attributed to the lack of CO 2 solubility data for IL structures other than imidazolium based ILs. We report here an automated high-throughput (HT) setup for the measurement of CO 2 solubility in room-temperature ILs (RTILs) combining six different anions and nine different cations for a total of 19 different specific ranges of RTILs. The HT setup first dispenses up to 200 µL of RTILs in a 96-well microtiter plate and then utilizes a robotic arm to measure cyclic voltammogram (CV) in each well using maneuverable Ag electrodes. The Cottrell analysis of the CO 2 reduction CV peak provides a direct measurement of CO 2 permeance in RTILs, which yields Henry’s constant from the estimated diffusion coefficient of CO 2 . Henry’s constants thus obtained are in very good agreement with those reported earlier. The measured CO 2 permeance and Henry’s constant of all RTILs seem to follow a first-order dependence on void fraction and a second-order dependence on electrostatic interaction between anion and cation of IL, with some synergistic dependence on the product of a void fraction and electrostatic interaction, making them two important descriptors for the design of novel ILs.

CO2 Solubility

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

Interfacial Adhesion Mechanism of Anionic Polyelectrolyte Brushes Induced by Oppositely Charged Macromolecular Counterions

Abstract Polyelectrolyte brushes are widely used as model systems for investigating electrostatic interactions at soft interfaces and to achieve exceptional lubrication properties. Previous studies have primarily focused on their behavior in the presence of multivalent counterions, which induce brush collapse, ionic crosslinking, and pronounced changes in interfacial structures. However, interactions between polyelectrolyte brushes and oppositely charged macromolecular counterions, such as polycations, remain poorly understood. Here, we prepare well-defined polystyrene sulfonate (PSS) brushes fabricated via surface-initiated grafting and employ surface forces apparatus measurements to investigate their interactions with oppositely charged polycations. In contrast to multivalent ions, polycations do not induce noticeable brush collapse, but instead generate significant adhesion between symmetric PSS brush layers. This adhesion increases with both contact time and applied load, eventually reaching a steady plateau, indicating that the interaction is governed not simply by electrostatic screening but by the gradual formation of polycation-mediated bridging under confinement. Furthermore, the introduction of monovalent Na+ ions disrupts the adhesive interaction even at very low concentrations, suggesting that the bridging function of the adsorbed polycation is highly sensitive to competitive ionic screening. In comparison, measurements with trivalent counterions reveal the expected brush collapse behavior but minimal dependence of adhesion on contact time, highlighting a clear mechanistic distinction from the polymeric counterion case. Collectively, these results demonstrate that molecular size, configurational restriction, and confinement-induced rearrangement, rather than charge valency alone, govern the interaction behavior of macromolecular counterions at brush interfaces. This work provides new insight into the molecular origins of adhesion and the regulation of interfacial interactions in charged polymer brush systems.

Park, Jinwoo [Argonne National Laboratory , , , ,;

Interactions between Per- and Polyfluoroalkyl Substances (PFAS) at the Water–Air Interface

Per- and polyfluoroalkyl substances (PFAS)--so-called “forever chemicals”--contaminate the drinking water of about 100 million people in the U.S. alone and are inefficiently removed by standard treatment techniques. A key property of these compounds that underlies their fate and transport and the efficacy of several promising remediation approaches is that they accumulate at the water−air interface. This phenomenon remains incompletely understood, particularly under conditions relevant to natural and treatment systems where water−air interfaces often carry significant loads of other organic contaminants or natural organic matter. To understand the impact of organic loading on PFAS adsorption, we carried out molecular dynamics simulations of PFAS at varying interfacial densities. We find that adsorbed PFAS form strong mutual interactions (attraction between perfluoroalkyl chains and electrostatic interactions among charged head groups) that give rise to ordered interfacial coatings. Furthermore, these interactions often involve near-cancellation of hydrophobic attraction and Coulomb repulsion. Our findings explain an apparent paradox whereby PFAS adsorption isotherms often suggest minimal mutual interactions while simultaneously displaying a high sensitivity to the composition and density of interfacial coatings. Consideration of the compounds present with PFAS at the interface has the potential to allow for more accurate predictions of fate and transport and the design of more efficient remediation approaches.

54 ENVIRONMENTAL SCIENCES

Unravelling the Glycan Code: Molecular Dynamics and Quantum Chemistry Reveal How O-Glycan Functional Groups Govern OgpA Selectivity in Mucin Degradation by Akkermansia muciniphila

Mucins, heavily O-glycosylated glycoproteins, are a key component of mucus, and certain gut microbiota, including Akkermansia muciniphila , can utilise mucin glycans as a carbon source. Akkermansia muciniphila produces the O-glycopeptidase enzyme OgpA, which cleaves peptide bonds at the N-terminus of serine (Ser) or threonine (Thr) residues carrying O-glycan substitutions, with selectivity influenced by the O-glycan functional groups. Using molecular dynamics (MD) simulations and quantum chemistry calculations, we explored how different O-glycan groups affect OgpA's selectivity. Our results show that peptides bind to the enzyme via hydrogen bonds, π–π interactions, van der Waals forces and electrostatic interactions, with key residues, including Tyr90, Val138, Gly176, Tyr210 and Glu91, playing important roles. The primary determinant of selectivity is the interaction between the peptide's functional group and the enzyme's binding cavity, while peptide–enzyme interface interactions are secondary. Quantum chemistry calculations reveal that OgpA prefers peptides with a lower electrophilic character. This study provides new insights into mucin degradation by gut microbiota enzymes, advancing our understanding of this critical biological process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Computational Investigation of the Chemical Bond between An(III) Ions and Soft-Donor Ligands

The chemical bonding of actinide ions with arene and borohydride ligands is explored via quantum chemical methods to understand how the transuranium elements interact with softdonor ligands. Specifically, the [An(C 6 Me 6 )(BH 4 ) 3 ] complexes (An = U, Np, and Pu) and their reduced congeners are studied. Density functional theory (DFT) shows that the metal–ligand interactions in the neutral complexes are governed by electrostatic interactions. Both DFT and complete active space (CASSCF) results show that as one moves from U to Pu, the 5f-orbitals are stabilized leading to a poorer energy match with the ligand orbitals. This contributes to progressively weaker metal-arene and metal-borohydride interactions across the series due to a decrease in energy-driven covalency. A reduction in orbital contributions to bonding is obtained for the transuranium-arene interactions as well. Upon reduction, the arene is reduced, forming a δ-bond. This causes the An–arene distances to contract by 0.1–0.2 Å compared to the neutral complexes. The ground state is assigned as the intermediate-spin state where the arene radical is antiferromagnetically coupled to the metal-centered f-electrons in Np and Pu. On the other hand, the ferromagnetically and antiferromagnetically coupled states are close in energy in the uranium complex, but do not mix when spin– orbit coupling is included using a state-interaction approach (SO-CASPT2). The population of the CASSCF δ*-antibonding natural orbital increases from U to Pu consistent with the increased An–arene distances, weaker interactions, and decreasing covalency across the series. Although the An–B distance increases by ca. 0.06 Å upon reduction, both the neutral and reduced species involve an An(III)–borohydride bond and as such are qualitatively similar. The Np complexes can be assigned to have slightly weaker bonding than the uranium analogs but are overall “uranium-like”. The Pu complexes are predicted to have less covalent contributions to bonding in both the Pu–arene and Pu–borohydride interactions; however, the Pu–arene interaction is predicted to be particularly weak.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Influence of shape on crystal structure and optical properties of heterocyclic conjugated molecules

Organic optoelectronics are increasingly important due to their tunablilty, flexibility, and solution processability. Tuning optical properties of these materials as solids relies on the balance of weak non-covalent interactions that dictate crystal structure, but are difficult to predict. Our research aims to improve our understanding of how electrostatic interactions can direct and facilitate intramolecular interactions that dictate emergent properties of crystalline materials. This paper focuses on exploring how multi-fused thiophene ring systems that are popular in modern organic optoelectronic materials impact intramolecular interactions, while also investigating the role of molecular shape. In these examples, the shape of heterocyclic systems correlate with the crystal structures: while the bent heterocyclics show no discrete and discernable intramolecular interactions, those with bent shapes interact cofacially with one of the electron poor ArF pendants by twisting the arylene ethynylene backbone. Two of the control molecules, which bear non-fluorinated benzyl ester substituents, show intramolecular edge-face interactions, and several of these molecules show clear polymorphic behavior. Furthermore, these findings further our understanding of how discrete interactions can be altered not only by electrostatics, but also by shape, allowing for increasingly nuanced control over the crystal structures and optical properties of optoelectronic materials.

Guzmán, Elisa

Solvent‐Phobic and Ionophilic Carboxylated Polythiophene Layer for Fluoride‐Rich Cathode Electrolyte Interphase

Abstract One focal area of contemporary organic mixed ionic‐electronic conductor (OMIEC) research relates to utilization of dual‐conductive properties to enhance the ion/electron transfer kinetics for energy storage applications. Insight regarding OMIEC response toward the electrolyte anion and solvent used in lithium‐ion batteries (LIBs), however, is limited. Here, for the first time, the solvent‐phobic and ionophilic (SP‐IP) properties of the OMIEC, poly[3‐(potassium‐4‐butanoate)thiophene‐2,5‐diyl] (P3KBT), are revealed through comprehensive evaluation and characterization. The solvent‐phobic characteristics arise from the cooperation of dispersive interaction, polar interaction, and hydrogen‐bonding between P3KBT and electrolyte solvent. The ionophilic nature is driven by electrostatic interactions between P3KBT side chain carboxylate groups and LiPF 6 , and the reversible electrochemical doping/de‐doping of the polythiophene backbone with PF 6 ⁻ . The SP‐IP properties induce formation of a LiF‐ rich , Li 2 CO 3 ‐ limited cathode electrolyte interphase (CEI) layer when a P3KBT coating layer is applied to the active material surface, significantly improving half‐cell life to over 1500 cycles at 2C.

Ren, Haoze [Department of Chemical and Bimolecular

Pillared Laminar Vermiculite Membranes with Tunable Monovalent and Multivalent Ion Selectivity

Effective membrane separation of Li + from Na + and Mg 2+ is crucial for lithium extraction from water yet challenging for conventional polymeric membranes. Two dimensional (2D) membranes with ordered laminar structures and tunable physicochemical properties offer distinctive ion-sieving capabilities promising for lithium extraction. Recently, phyllosilicates are introduced as abundant and cost-effective source materials for such membranes. However, their water instability and low inherent ion transport selectivity hinder practical applications. Herein, a new class of laminar membranes with excellent stability and tunable ion sieving is reported by incorporating inorganic alumina pillars into vermiculite interlayers. Crosslinking vermiculite flakes with alumina pillars significantly strengthens interlamellar interactions, resulting in robust water stability. Doping of Na + before the pillaring process reverses the membrane's surface charge, substantially boosting Li + separation from multivalent cations via electrostatic interactions. Lithium extraction is often complicated by the presence of co-existing monovalent cations (e.g., Na + ) at higher concentrations. Here, by introducing excess Na + into the membrane after the pillaring process, the separation of Li + from monovalent cations is enhanced through steric effects. This work realizes both monovalent/multivalent and monovalent/monovalent selective ion sieving with the same membrane platform. A separation mechanism is proposed based on Donnan exclusion and size exclusion, providing new insights for membrane design for resource recovery applications.

2D materials

Multi-scale computational screening and mechanistic insights of cyclic amines as solvents for improved lignocellulosic biomass processing

A computational screening workflow for the efficient deconstruction of cellulose, lignin and hemicellulose fractions of lignocellulosic biomass using cyclic amines as solvents. Lignocellulosic biomass is a promising feedstock for production of affordable fuels and chemicals from renewable resources. Effective solubilization and subsequent deconstruction of its cellulose, hemicellulose, and lignin fractions is essential for the viability of future biorefineries. This study used quantum chemistry-based equilibrium thermodynamics methods to evaluate the potential of 650 cyclic amines to solubilize cellulose, hemicellulose, and lignin. The activity coefficients of solvent - biopolymer interactions were predicted using the COSMO-RS (COnductor-like Screening MOdel for Real Solvents) method and used to identify cyclic amines that can efficiently dissolve and extract selective fractions of biopolymers during biomass pretreatment. Among the 650 cyclic amines, 1-piperazineethanmaine was predicted to be an effective solvent for extracting all three polymers and was experimentally shown to achieve the highest lignin removal (97.1%). Non-covalent interaction, reduced density gradient and quantum chemical calculations were performed to elucidate the dissolution mechanism of lignin, cellulose and hemicellulose and gain further molecular level insights into the interactions between the cyclic amines and biomass polymers that promote efficient solubilization and extraction. These analyses indicated that 1-piperazineethanmaine and 1-methylimidazole make noncovalent van der Waals, electrostatic interactions and hydrogen bonding with lignin, leading to enhanced lignin removal, while the strong intramolecular hydrogen bonding interactions in cellulose and hemicellulose result in weaker solvent-biopolymer interactions. Overall, the computational approach provided an efficient method for identifying cyclic amines tailored for optimal biomass pretreatment and resulted in the identification of a potential new class of solvents for effective biomass pretreatment.

Kumar, Nikhil

Binding Free Energy Analysis of Colicin D, E3 and E8 to Their Respective Cognate Immunity Proteins Using Computational Simulations

Colicins are antimicrobial proteins produced by bacteria for the purpose of destroying neighboring bacteria. Colicin activity is neutralized by a specific cognate immunity protein in order to protect the host. This study investigates the structural and binding mechanisms underlying the interaction of colicin-D, -E3 and -E8 to their respective immunity proteins (ImD, Im3 and Im8) using structure prediction, molecular dynamics (MD) simulations and MM-PBSA approach of free energy calculations. High-confidence colicin-immunity (Col-Im) complex structures predicted using AlphaFold2 were subjected to MD simulations of 150 ns with GROMACS and were analyzed for the binding free energy calculation using gmx_MMPBSA. Results showed that the complex of Col_E3-Im3 exhibited the most favorable binding free energy, driven by strong van der Waals and electrostatic interactions. Col_D-ImD and Col_E8-Im8 also showed the favorable binding. Electrostatics and hydrogen bonding emerged as a key factor driving binding and stability, while polar solvation acted as a destabilizing factor across all systems. These outcomes provide an understanding of the molecular mechanisms of Col-Im systems, with potential applications for developing natural antimicrobials for food safety.

Biochemistry & Molecular Biology

Interacting dust grains in complex plasmas: Ion wake formation and the electric potential

Dust grains have been used as minimally invasive probes to determine plasma parameters including the plasma density, temperature, and electric field in a plasma discharge. However, the dust grains in a plasma generate local potential disturbances due to the collection of charge and the subsequent electrostatic interactions between the dust and charged plasma particles. Dust grains in close proximity to one another exhibit interesting non-reciprocal interactions and self-organize into structures such as one-dimensional filamentary chains, two-dimensional “zigzags,” and three-dimensional helices, among others. The formation of these structures suggests that although the dust grains may be less invasive than traditional plasma probes, the disturbance to the local plasma environment introduced by dust grains is non-trivial. Commonly used analytic forms of the electric potential describing complex plasmas have failed to resolve the near-dust region, and as a result are insufficient to provide insight about the formation of complex dust structures. Here, we use an N-body simulation to compute the electric potential from ion densities near various dust grain configurations. We provide an alternative description to the standard analytic model for the electric potential of dust and ion wakes based on a Gaussian shaped cloud of ions. The electric potential obtained from simulations is used to identify minimum energy configurations for two and three dust grains. It is further demonstrated that the minimum potential region identified for N dust grains and their associated ion wakes does not predict the minimum-energy configuration of N + 1 dust grains.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY