Search NASA⌕ Search

SEARCH · Search NASA

Results for “ion-ion correlations”

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.

Direct Experimental Observations of Ion Distributions during Overcharging at the Muscovite–Water Interface by Adsorption of Rb + and Halides (Cl – , Br – , I – ) at High Salinity

Classical electric double layer (EDL) models have been widely used to describe ion distributions at charged solid-water interfaces in dilute electrolytes. However, the chemistry of EDLs remains poorly constrained at high ionic strength where ion-ion correlations control non-classical behavior such as overcharging, i. e., the accumulation of counter-ions in amounts exceeding the substrate's surface charge. Here, we provide direct experimental observations of correlated cation and anion distributions adsorbed at the muscovite (001)-aqueous electrolyte interface as a function of dissolved RbBr concentration ([RbBr]=0.01–5.8 M) using resonant anomalous X-ray reflectivity. Our results show alternating cation-anion layers in the EDL when [RbBr]≳100 mM, whose spatial extension (i. e., ~20 Å from the surface) far exceeds the dimension of the classical Stern layer. Comparison to RbCl and RbI electrolytes indicates that these behaviors are sensitive to the choice of co-ion. This new in-depth molecular-scale understanding of the EDL structure during transition from classical to non-classical regimes supports the development of realistic EDL models for technologies operating at high salinity such as water purification applications or modern electrochemical storage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanisms of proton transport in aqueous acid solutions

The fundamental understanding of proton transport, specifically the Grotthuss-like mechanism, is critical for many technological applications. The present study presents a comprehensive analysis of proton transport in aqueous solutions of sulfuric and phosphoric acids using quasielastic neutron and light scattering, broadband dielectric spectroscopy, rheology, pulsed-field gradient-NMR measurements, and ab initio molecular dynamic simulations. Results show that in all systems, proton transport occurs through short “jumps” of ~0.5 Å that are faster than structural relaxation. However, proton hopping appears to be coupled to structural relaxation in aqueous solutions of sulfuric acid, while these processes are decoupled in phosphoric acid. Neutron scattering indicates that all protons have the same fast mobility in phosphoric acid systems, while at least one proton per sulfuric molecule remains slower than other protons. Analysis reveals correlated proton jumps, but these correlations suppress conductivity, suggesting that expected Grotthuss-like enhancement of conductivity is unlikely in bulk liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of Solvent Rearrangement on Mg2+ Solvation Structures in Dimethoxyethane Solutions Using Multimodal NMR Analysis

One of the main impediments faced for predicting important emergent properties of a multivalent electrolyte (such as conductivity and electrochemical stability) is the lack of quantitative analysis of ion-ion and ion-solvent interactions, which manifest in solvation structures and dynamics. In particular, the role of ion-solvent interactions is still unclear in cases where the strong electric field from multivalent cations can influence structural evolution of lower permittivity organic solvent molecules on solvation structure. These influences are commonly refered to as electrostriction pheonmena. Using quantitative 1H, 19F, 17O, and 25Mg NMR together with 19F nuclear spin relaxation and diffusion measurments, we find an unusual correlation between ion concentration and solvation structure of Mg(TFSI)2 salt in dimethoxyethane (DME) solution. The domiant solvation structure evolves from contact ion pairs (i.e. [Mg(TFSI)(DME)1-2]+) to fully solvated clusters (i.e. [Mg(DME)3]2+) as salt concentration increases or as temperature decreases. This transition is coupled to a phase separation, which we study here between 0.06 M and 0.36 M. Subsequent analysis is based on an explanation of the electrostriction phenomena and the competion between solvent molecules and TFSI anions for cation coordination.

Chen, Ying↗

Ion Conductivity and Correlations in Model Salt-Doped Polymers: Effects of Interaction Strength and Concentration

Correlated anion and cation motion can greatly reduce the overall ion conductivity in electrolytes versus the ideal conductivity calculated based on the diffusion constants alone. Using coarse-grained molecular dynamics simulations, we calculate conductivity and the degree of uncorrelated ion motion in salt-doped homopolymers and block copolymers as a function of concentration and interaction strengths. Calculating conductivity from ion mobility under an applied electric field increases accuracy versus the typical use of fluctuation dissipation relationships in equilibrium simulations. In typical electrolytes, correlation in cation-anion motion is often expected to be reduced at low ion concentrations. Yet, for these polymer electrolytes with strong ion-polymer and ion-ion interactions, we find cation-anion motion is more correlated at lower concentrations when other variables are held constant. We show this phenomenon is related to the slower ion cluster relaxation rate at low concentrations rather than the static spatial state of ion aggregation or the fraction of free ions.

36 MATERIALS SCIENCE↗

Impact of Ion–Ion Correlated Motion on Salt Transport in Solvated Ion Exchange Membranes

The influence of dynamical ion-ion correlations and ion pairing on salt transport in ion exchange membranes remain poorly understood. In this study, we use the framework of Onsager transport coefficients within atomistic molecular dynamics sim- ulations to study the impact of ion-ion correlated motion on salt transport in hydrated polystyrene sulfonate membranes and compare with the results from aqueous salt so- lutions. At sufficiently high salt concentrations, cation-anion dynamical correlations exert a significant influence on both salt diffusivities and conductivities. Anion-anion distinct correlations, arising from the imbalance between the concentration of free (mo- bile) cations and anions, and the retarding effect of the fixed charge groups on cations, proves to be an additional important feature for polymer membranes. Furthermore, our results demonstrate that dynamical correlations should become an important consideration in experimental measurements of salt diffusivities and conductivities for non-dilute salt solutions in polymer membranes.

36 MATERIALS SCIENCE↗

Ion correlations drive charge overscreening and heterogeneous nucleation at solid–aqueous electrolyte interfaces

Significance Ion distributions at charged solid–water interfaces, referred to as the electrical double layer (EDL), are poorly understood at high ion concentrations, in part due to the lack of molecular-scale descriptions of the interactions between adsorbed hydrated ions. Here, direct visualization of the salinity-dependent evolution of EDL structure reveals molecular origins of nonclassical transformation of the EDL, in which charge overscreening and heterogeneous nucleation are driven by ion–ion correlations at the interfaces. This manifestation of the atomistic basis of nonclassical behaviors provides a much-needed understanding of the impact of ion cooperativity at charged interfaces for the development of predictive models for element transport in natural environments and advanced technologies for material growth and synthesis in saline environments.

36 MATERIALS SCIENCE↗

Modeling electrokinetic flows with the discrete ion stochastic continuum overdamped solvent algorithm

In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the discrete ion stochastic continuum overdamped solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Herein, several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced-charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a “dry diffusion” approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy.

97 MATHEMATICS AND COMPUTING↗

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hard-photon-triggered jets in 𝑝−𝑝 and 𝐴−𝐴 collisions

An investigation of high-transverse-momentum (high-𝑝 𝑇 ) photon-triggered jets in proton-proton (𝑝−𝑝) and ion-ion (𝐴−𝐴) collisions at $\sqrt{s_{NN}}$=0.2 and 5.02TeV is carried out, using the multistage description of in-medium jet evolution. Monte Carlo simulations of hard scattering and energy loss in heavy-ion collisions are performed using parameters tuned in a previous study of the nuclear modification factor (𝑅 𝐴⁢𝐴 ) for inclusive jets and high-𝑝𝑇 hadrons. We obtain a good reproduction of the experimental data for photon-triggered jet 𝑅 𝐴⁢𝐴 , as measured by the ATLAS detector, the distribution of the ratio of jet to photon 𝑝 𝑇 (𝑋 𝐽⁢𝛾 ), measured by both CMS and ATLAS, and the photon-jet azimuthal correlation as measured by CMS. We obtain a moderate description of the photon-triggered jet 𝐼 𝐴⁢𝐴 , as measured by STAR. A noticeable improvement in the comparison is observed when one goes beyond prompt photons and includes bremsstrahlung and decay photons, revealing their significance in certain kinematic regions, particularly at 𝑋 𝐽⁢𝛾 >1. Moreover, azimuthal angle correlations demonstrate a notable impact of bremsstrahlung photons on the distribution, emphasizing their role in accurately describing experimental results. This work highlights the success of the multistage model of jet modification to straightforwardly predict (this set of) photon-triggered jet observables. This comparison, along with the role played by bremsstrahlung photons, has important consequences on the inclusion of such observables in a future Bayesian analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The ion-ion correlations in organic ionic plastic crystal

Organic ionic plastic crystals (OIPCs) are emerging as promising electrolyte materials for solid-state batteries. However, despite the fast ionic diffusion, OIPCs exhibit relatively low DC conductivity in solid phases caused by strong ion-ion correlations that suppress charge transport. To understand the origin of this suppression, we performed a study of ion dynamics in the OIPC 1-Ethyl-1-methylpyrrolidinium bis (trifluoromethyl sulfonyl) imide [P 12 ][TFSI] utilizing dielectric spectroscopy, light scattering, and Nuclear Magnetic Resonance diffusometry. Comparison of the results obtained in this study with the published earlier results on an OIPC with a completely different structure (Diethyl(methyl)(isobutyl)phosphonium Hexafluorophosphate [P 1,2,2,4 ][PF 6 ]) revealed strong similarities in ion dynamics in both systems. Unlike DC conductivity, which may drop more than ten times between melted and solid phases, diffusion of anions and cations remains high and does not show strong changes at phase transition. The conductivity spectra in the broad frequency range demonstrate unusual shapes in solid phases with an additional step separating fast local ion motions from suppressed long-range charge diffusion controlling DC conductivity. We suggested that in solid phases, anions and cations can jump only between the specific ion sites defined by the crystalline structure. These constraints lead to strong cation-cation and anion-anion correlations strongly suppressing long-range charge transport.

36 MATERIALS SCIENCE↗

Verification Testing of OLI Systems Mixed Solvent Electrolyte Model for the Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 ) System to High Ionic Strength at 25°C

This technical report summarizes model verification results and summary statistics for 41 evaporite mineral solubility cases evaluated by Savannah River National Laboratory using OLI Systems’ aqueous electrolyte thermodynamic modeling software. The 41 verification cases containing a total of 60 solubility curves comprise mineral solubility data from low to high ionic strength at 25°C for the eight-component system Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 O as reported by Harvie et al. (1984). Thermodynamic calculations were executed using OLI Systems’ Stream Analyzer computation module within the OLI Studio software platform (Ver. 11.0, Rev. 11.0.1.9). The Mixed Solvent Electrolyte (MSE) thermodynamic framework was chosen for this investigation because of its superiority in modeling high ionic-strength inorganic salt solutions and actinide redox chemistry and solubility, both of which are relevant to the geological repository conditions at the Waste Isolation Pilot Plant in Carlsbad, New Mexico. Mineral solubility data in various inorganic salt solutions were digitized and extracted from figures generated by Harvie et al. (1984). For each of the 60 solubility curves, a case-specific chemistry model and input file were generated in OLI Studio using OLI Stream Analyzer and the MSE (H 3 O + ion) public databank provided by OLI Systems. Model simulation results were exported to Microsoft Excel to calculate summary statistics and to generate graphs comparing the OLI model predictions to the solubility data. Summary statistics include residuals (model – data) and concordance (accuracy × precision, where precision is indicated by the Pearson correlation coefficient and accuracy accounts for bias and scale differential). Private databanks were not developed, and activity coefficient model regressions were not performed to improve OLI model fits to the data. Of the 41 model verification plots, 83% have a mean of the percent residuals less than or equal to 25%. Similarly, 75% display a concordance greater than or equal to 0.75. Only seven of the 41 verification plots fail to show good agreement between the model and data. Of these seven, three are relevant to the WIPP repository because they involve the Mg-OH-Cl-SO 4 -CO 3 aqueous system. The remaining four address salt solubilities at the pH extremes (strong acid and strong base). It should be noted that in two of the three Mg-OH-Cl-SO 4 -CO 3 system cases, the regressed Harvie et al. (1984) solubility curve also deviated from the data. Lack of agreement between the OLI model-predicted solubility curves and the data is attributable to one or more of the following: specific solid species are not included in the OLI MSE databank; there is significant variation among the different solubility datasets chosen by Harvie et al. (1984); the OLI MSE model’s thermodynamic parameters were determined using different solubility datasets; and the activity coefficient parameters for certain relevant ion-ion and ion-molecule pairs have not been optimized via data regression. Two recommendations for future work are to (1) evaluate solubility data for the Mg-OH-Cl-SO 4 -CO 3 system at high ionic strength and, if necessary, develop a private OLI MSE database that includes missing species and, where necessary, regressed standard state properties and interaction parameters; (2) perform similar verification testing of the OLI model for actinide solubility data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effect of Ionic Core on the Properties of Non-Isothermal Plasmas

This study investigates the effects of ionic cores on non-isothermal plasmas using a novel ion-ion interaction potential that incorporates screening effects from both ion cores and exchange-correlation interactions. Our findings indicate that with increasing distance, the effective potential approaches a Yukawa-like screening potential, while at shorter distances, strong electron binding weakens screening. The different values of the cutoff radius ' cut and core edge steepness α significantly influence the potential behavior and radial distribution functions (RDFs). Higher coupling parameters (Γ i ) strengthen the electron-ion interactions, leading to deeper potential wells and more pronounced non-ideality corrections. Increasing Θ decreases the absolute values of non-ideality corrections, indicating fewer interactions in the system. A larger cutoff radius r cut due to weaker screening effects. As Γ i at a fixed parameter also reduces corrections increases, non-ideality corrections grow, reflecting stronger coupling. The results show the importance of taking into account the ion core effects in dense non-isothermal plasma research.

dense plasmas↗

Perturbation theory of structure in classical liquid mixtures: Application to metallic systems near phase separation

The partial structure factors of classical simple liquid mixtures near phase separation are dicussed. The theory is developed for particles interacting through pair potentials, and is thus appropriate both to insulating fluids, and also to metallic systems if these may be described by an effective ion-ion pair interaction. The motivation arose from consideration of metallic liquid mixtures, in which resistive anomalies have been observed near phase separation. A mean field theory correction appropriate to 3 pair potential for the effects of correlated motions in the reference fluid is studied. The work is cast in terms of functions which are closely related to the direct correlation functions of Ornstein and Zernike. The results are qualitatively in accord with physical expectations. Quantitative agreement with experiment seems to turn on the selection of the hard core reference potential in terms of the metallic effective pair potential. It is suggested that the present effective pair potentials are perhaps not properly used to calculate the metallic structure factors at long wavelength.

Henderson, R. L.↗

Adhesion of a bimetallic interface

The Hohenberg-Kohn and Kohn-Sham formalisms are used to examine binding (binding energy as a function of separation) for combinations of the simple metals Al(111), Zn(0001), Mg(0001), and Na(110) in contact. Similar metal contacts between Al, Zn, Mg, and Na are examined self-consistently in an ab initio calculation using the Kohn-Sham formalism. Crystallinity is included using the Aschroft pseudopotential via first order perturbation theory for the electron-ion interaction; and the ion-ion interaction is included exactly via a lattice sum. Binding energy was determined both in the local-density approximation and including gradient corrections to the exchange and correlation energy. Binding was found in all cases. In dissimilar metal contacts, interfacial bonding was greater than that in the weaker material predicting the possibility of metallic transfer. The nonzero position of the energy minimum in like metal contacts is explained in terms of consistency between the Ashcroft pseudopotential and the bulk charge density. Good agreement with experimental surface energies is obtained in the self-consistent calculation when nonlocal terms are included.

Ferrante, J.↗

Theory and simulations of broadband electrostatic noise in the geomagnetic tail

The excitation mechanism for broadband electrostatic noise (BEN) and the effects of BEN particles in the geomagnetic tail are examined using the linear analysis theory and particle simulations. The linear theory for electrostatic instabilities is discussed. The plasma sheet particle population is simulated using counter-streaming cold ion beams, and warm ions and electrons. The ion-ion instability, ion-acoustic mode, and the electrostatic ion cyclotron harmonic waves are studied. The velocity distributions, electric field intensity, and electron plasma waves for the plasma sheet boundary layer are evaluated. The frequency wave spectrum and particle distributions are computed and analyzed. The conditions for the two simulations, which differ only in beam ion drift speed, are described; it is observed that in the first simulation the dominate modes propagate parallel to the magnetic field and in the second simulation the propagation modes are oblique. The simulation data reveal that when beam temperature is smaller than plasma sheet temperature ion-acoustic and ion-ion instabilities grow to large amplitudes heating both electrons and ions. The data are compared to ISEE-1 observations and good correlation is obtained.

Ashour-Abdalla, M.↗