Search NASA⌕ Search

SEARCH · Search NASA

Results for “electrostatic interaction”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2019

Connection between Mott physics and crystal structure in a series of transition metal binary compounds

Abstract The choice that a solid system “makes” when adopting a crystal structure (stable or metastable) is ultimately governed by the interactions between electrons forming chemical bonds. Here we analyze six prototypical binary transition metal compounds and shed light on the connection between Mott physics and the behavior of the energy as a function of the spatial arrangement of the atoms in these systems. Remarkably, we find that the main qualitative features of this complex behavior in the Mott phase of these systems can be traced back to the fact that the strong d -electron correlations influence substantially the charge transfer mechanism, which, in turn, controls the electrostatic interactions. This result advances our understanding of the influence of strong correlations on the crystal structure, opens a new avenue for extending structure prediction methodologies to strongly correlated materials, and paves the way for predicting and studying metastability and polymorphism in these systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Decay of Electrostatic Force of Dust Particles on Photovoltaic Modules

Photovoltaic (PV) energy yield loss due to solar module soiling has become increasingly important as solar module deployment is now at the hundreds of gigawatts scale and continues to grow rapidly. We have reported on direct measurements-using atomic force microscopy (AFM)-of strong electric-field-induced attraction and adhesion force (Fes) of dust particles onto solar panel that are 1 to 2 orders of magnitude stronger than the van der Waals and water capillary forces, corroborated by observing the increase in system voltage-induced soiling rate. Here, we report another characteristic of Fes on soiling-long lasting or slow decay after turning off the high voltage applied to solar panels. The Fes decay time varies in a wide time range of 1 to 10 hours, depending on two factors: 1) either/both the cell or/and particle were charged with high voltage before the voltages were turned off; and 2) how the cell was connected to the ground after the voltage was turned off-either connected through the power supply electronics, directly connected to the ground, or electrically floated. The Fes decay is understood in terms of 1) net electrical charge dissipations in both particle and cell, 2) thermal disordering of dipole polarization in the module glass dielectrics, and 3) charge redistribution by the electrostatic interaction of particle and module glass. This long-lasting Fes for hours can significantly affect the solar panel soiling after sunset, especially combined with water condensation.

14 SOLAR ENERGY↗
DOE OSTI2025

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

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

Adsorption↗
DOE OSTI2020

Effect of Dispersion Medium Composition and Ionomer Concentration on the Microstructure and Rheology of Fe–N–C Platinum Group Metal-free Catalyst Inks for Polymer Electrolyte Membrane Fuel Cells

In this paper, we present an investigation of the microstructure and rheological behavior of catalyst inks consisting of Fe–N–C platinum group metal-free catalysts and a perfluorosulfonic acid ionomer in a dispersion medium (DM) of water and 1-propanol ( n PA). The effects of the ionomer-to-catalyst (I/C) ratio and weight percentage of water (H 2 O %) in the DM on the ink microstructure were studied. Steady-shear and dynamic-oscillatory-shear rheology, in combination with synchrotron X-ray scattering, was utilized to understand interparticle interactions and the level of agglomeration of the inks. In the absence of the ionomer, the inks were significantly agglomerated, approaching a gel-like microstructure for catalyst concentrations as low as 2 wt %. The effect of H 2 O % in the DM on particle agglomeration was found to vary with particle concentration. In concentrated inks (≥2 wt % catalyst), increasing H 2 O % was found to increase agglomeration because of the hydrophobic nature of the catalysts. In dilute inks (<1 wt % catalyst), the trend was reversed with increasing H 2 O %, suggesting that electrostatic interactions are dominating the behavior. In inks with 5 wt % catalyst, the addition of an ionomer was found to significantly stabilize the catalyst against agglomeration. Maximum stability was observed at 0.35 I/C for all DM H 2 O % studied. At high ionomer concentrations (I/C > 0.35), interesting differences were observed between n PA-rich inks (H 2 O % ≤ 50%) and H 2 O-rich (82% H 2 O) inks. The n PA-rich inks remained predominantly stable—ink viscosity only weakly increased with I/C and the Newtonian behavior was maintained for I/C up to 0.9. In contrast, the H 2 O-rich inks exhibited a significant increase in viscoelasticity with increasing I/C, suggesting flocculation of the catalyst by the ionomer. These differences suggest that the nature of the interactions between the ionomer and catalyst is highly dependent on the H 2 O % in the DM.

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

Partitioning of Seven Different Classes of Antibiotics into LPS Monolayers Supports Three Different Permeation Mechanisms through the Outer Bacterial Membrane

The outer membrane (OM) of Gram-negative (G-) bacteria presents a barrier for x classes of antibacterial agents. Lipopolysaccharide (LPS), present in the outer leaflet of the OM, is stabilized by divalent cations and is considered to be the major impediment for antibacterial agent permeation. However, the actual affinities of major antibiotic classes toward LPS have not yet been determined. In the present work, we use Langmuir monolayers formed from E. coli Re and Rd types of LPS to record pressure-area isotherms in the presence of antimicrobial agents. Our observations suggest three general types of interactions. Firstly, some antimicrobials demonstrated no measurable interactions with LPS. This lack of interaction in the case of cefsulodin, a third-generation cephalosporin antibiotic, correlates with its low efficacy against G- bacteria. Ampicillin and ciprofloxacin also show no interactions with LPS, but in contrast to cefsulodin, both exhibit good efficacy against G- bacteria, indicating permeation through common porins. Secondly, we observe substantial intercalation of the more hydrophobic antibiotics novobiocin, rifampicin, azithromycin, and telithromycin into relaxed LPS monolayers. These largely re-partition back to the subphase with monolayer compression. We find that the hydrophobic area, charge, and dipole all show correlations with both the mole fraction of antibiotic retained in the monolayer at the monolayer-bilayer equivalence pressure and the efficacies of these antibiotics against G- bacteria. Thirdly, amine-rich gentamicin and the cationic antimicrobial peptides polymyxin B and colistin show no hydrophobic insertion but are instead strongly driven into the polar LPS layer by electrostatic interactions in a pressure-independent manner. Their intercalation stably increases the area per molecule (by up to 20%), which indicates massive formation of defects in the LPS layer. These defects support a self-promoted permeation mechanism of these antibiotics through the OM, which explains the high efficacy and specificity of these antimicrobials against G- bacteria.

lateral pressure↗
DOE OSTI2021

ChAdOx1 interacts with CAR and PF4 with implications for thrombosis with thrombocytopenia syndrome

Vaccines derived from chimpanzee adenovirus Y25 (ChAdOx1), human adenovirus type 26 (HAdV-D26), and human adenovirus type 5 (HAdV-C5) are critical in combatting the severe acute respiratory coronavirus 2 (SARS-CoV-2) pandemic. As part of the largest vaccination campaign in history, ultrarare side effects not seen in phase 3 trials, including thrombosis with thrombocytopenia syndrome (TTS), a rare condition resembling heparin-induced thrombocytopenia (HIT), have been observed. This study demonstrates that all three adenoviruses deployed as vaccination vectors versus SARS-CoV-2 bind to platelet factor 4 (PF4), a protein implicated in the pathogenesis of HIT. We have determined the structure of the ChAdOx1 viral vector and used it in state-of-the-art computational simulations to demonstrate an electrostatic interaction mechanism with PF4, which was confirmed experimentally by surface plasmon resonance. These data confirm that PF4 is capable of forming stable complexes with clinically relevant adenoviruses, an important step in unraveling the mechanisms underlying TTS.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

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

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

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

Molecular Dynamics and Free Energy Calculations of Dicyclohexano-18-crown-6 Diastereoisomers with Sm 2+ , Eu 2+ , Dy 2+ , Yb 2+ , Cf 2+ , and Three Halide Salts in Tetrahydrofuran and Acetonitrile Using the AMOEBA Force Field

With the continual development of lanthanides (Ln) in current technological devices, an efficient separation process is needed that can recover greater amounts of these rare elements. Dicyclohexano-18-crown-6(DCH18C6) is a crown ether that may be a promising candidate for Ln separation, but additional research is required. As such, molecular dynamics(MD) simulations have been performed on four divalent lanthanide halide salts(Sm 2+ , Eu 2+ , Dy 2+ , and Yb 2+ ) and one divalent actinide halide salt (Cf 2+ ) bound to three diastereoisomers of DCH18C6. Dy 2+ , Yb 2+ , Cf 2+ , DCH18C6, and tetrahydrofuran (THF) solvent were parameterized for the AMOEBA polarizable force field for the first time, whereas existing parameters for Sm 2+ and Eu 2+ were utilized from our previous efforts. A coordination number (CN) of six for Ln 2+ /An 2+ –O solvated in THF indicated that the cations interacted almost entirely with the oxygens of the polyether ring. A CN of one for Ln 2+ /An 2+ -N solvated in acetonitrile for systems containing iodide suggested that theN atom of acetonitrile was competitive with I – for cation interactions. Fluctuation between five and six CNs for Dy 2+ and Yb 2+ suggested that although the cations remained in the polyether ring, the size of the ring may not be an ideal fit as these cations possess comparatively smaller ionic radii. Gibbs binding free energies of Sm 2+ in all DCH18C6 diastereoisomers solvated in THF were calculated. The binding free energy of the cis-syn-cis diastereoisomer was the most favorable, followed by cis-anti-cis, and then trans-anti-trans. Lastly, two major types of conformation were observed for each diastereoisomer that were related to the electrostatic interactions and charge density of the cations.

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

Long-lasting strong electrostatic attraction and adhesion forces of dust particles on photovoltaic modules

Photovoltaic (PV) energy-yield loss due to solar module soiling has become increasingly important as solar module deployment is now at the hundreds of gigawatts scale and continues to grow rapidly. The electrostatic attraction and adhesion force (F es ) of dust particles by the high voltages from solar panels have been reported to be 1 to 2 orders of magnitude stronger than the van der Waals and water capillary forces, which has been corroborated by observing the increase in system voltage-induced soiling rate by setting up an outdoor test. Here, we report another characteristic of F es on soiling - long-lasting or slow decay after turning off the high voltage applied to solar panels. The F es decay time varies across a wide time range of 1 to 10 h, depending on two factors: 1) whether the cell, the particle, or both were charged with high voltage before the voltages were turned off; and 2) how the cell was connected to the ground after the voltage was turned off - either connected through the power-supply electronics, directly connected to the ground, or electrically floated. The Fes decay is understood in terms of 1) net electrical charge dissipations in both particle and cell, 2) thermal disordering of dipole polarization in the module glass dielectrics, and 3) charge redistribution by the electrostatic interaction of particle and module glass. This long-lasting F es for hours can affect the solar panel soiling after sunset, and it can have an even greater effect when combined with water condensation at night.

14 SOLAR ENERGY↗
DOE OSTI2021

Fast and Accurate Machine Learning Strategy for Calculating Partial Atomic Charges in Metal–Organic Frameworks

Computational high-throughput screening using molecular simulations is a powerful tool for identifying top-performing metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges are often required to model the electrostatic interactions between the MOF and the adsorbate, especially when the adsorption involves molecules with dipole or quadrupole moments such as water and CO 2 . Although ab initio methods can be used to calculate accurate partial atomic charges, these methods are impractical for screening large material databases because of the high computational cost. We developed a random forest machine learning model to predict the partial atomic charges in MOFs using a small yet meaningful set of features that represent both the elemental properties and the local environment of each atom. The model was trained and tested on a collection of about 320 000 density-derived electrostatic and chemical (DDEC) atomic charges calculated on a subset of the Computation-Ready Experimental Metal–Organic Framework (CoRE MOF-2019) database and separately on charge model 5 (CM5) charges. The model predicts accurate atomic charges for MOFs at a fraction of the computational cost of periodic density functional theory (DFT) and is found to be transferable to other porous molecular crystals and zeolites. In conclusion, a strong correlation is observed between the partial atomic charge and the average electronegativity difference between the central atom and its bonded neighbors.

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

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

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

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2022

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

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

Lunar dust↗
DOE OSTI2022

Noncoordinating Secondary Sphere Ion Modulates Supramolecular Clustering of Lanthanides

The role of counterions in molecular recognition of lanthanides is underexplored, especially when they exhibit weak interactions with the metal cations. Here, we report a complementary and comprehensive investigation integrating theoretical calculations with X-ray absorption fine structure spectroscopy, dynamic light scattering, and small-angle X-ray scattering to reveal atomic-scale structural features beyond the immediate coordination sphere of a system used for rare-earth element separations. Here, our results indicate the formation of an unusual T-shaped outer-sphere lanthanide complex, containing two ligands and two nitrate ions in the first coordination sphere, whereas the third nitrate is weakly coordinated and resides in the second shell. This unique structural arrangement causes inhomogeneous charge distribution, leading to self-assembly of the complexes into larger nanoclusters through sterically directed electrostatic interactions in the nonpolar medium. Our findings point to the importance of "noncoordinating" anions in defining the degree of supramolecular aggregation and ion cluster assembly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

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

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

Banin, A.↗
DOE OSTI2023

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics

The chemical equilibrium between self-ionized and molecular water dictates the acid–base chemistry in aqueous solutions, yet understanding the microscopic mechanisms of water self-ionization remains experimentally and computationally challenging. Herein, Density Functional Theory (DFT)–based deep neural network (DNN) potentials are combined with enhanced sampling techniques and a global acid–base collective variable to perform extensive atomistic simulations of water self-ionization for model systems of increasing size. The explicit inclusion of long-range electrostatic interactions in the DNN potential is found to be crucial to accurately reproduce the DFT free energy profile of solvated water ion pairs in small (64 and 128 H 2 O) cells. The reversible work to separate the hydroxide and hydronium to a distance S is found to converge for simulation cells containing more than 500 H 2 O, and a distance of ∼ 8 Å is the threshold beyond which the work to further separate the two ions becomes approximately zero. The slow convergence of the potential of mean force with system size is related to a restructuring of water and an increase of the local order around the water ions. Calculation of the dissociation equilibrium constant illustrates the key role of long-range electrostatics and entropic effects in the water autoionization process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2013

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

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

Johnson, Scott↗