Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userUnderscreening and hidden ion structures in large scale simulations of concentrated electrolytes
56
0
0.0
(
0
)
Added by
Emily Krucker-Velasquez
Publication date
2021
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The electrostatic screening length predicted by Debye-Huckel theory decreases with increasing ionic strength, but recent experiments have found that the screening length can instead increase in concentrated electrolytes. This phenomenon, referred to as underscreening, is believed to result from ion-ion correlations and short-range forces such as excluded volume interactions among ions. We use Brownian Dynamics to simulate a version of the Restrictive Primitive Model for electrolytes over a wide range of ion concentrations, ionic strengths, and ion excluded volume radii for binary electrolytes. We measure the decay of the charge-charge correlation among ions in the bulk, and compare it against scaling trends found experimentally and determined in certain weak coupling theories of ion-ion correlation. Moreover, we find that additional large scale ion structures emerge at high concentrations. In this regime, the frequency of oscillations computed from the charge-charge correlation function is not dominated by electrostatic interactions but rather by excluded volume interactions and with oscillation periods on the order of the ion diameter. We also find that the nearest neighbor correlation of ions sharing the same charge transitions from negative at small concentrations to positive at high concentrations, representing the formation of small, like-charge ion clusters. We conclude that the increase in local charge density due to the formation of these clusters and the topological constraints of macroscopic charged surfaces can help explain the degree of underscreening observed experimentally.
rate research
Read More
Ion-pairing is commonly considered as a culprit for the reduced ionic conductivity in polymer electrolyte systems. However, this simple thermodynamic picture should not be taken literally, as ion-pairing is a dynamical phenomenon. Here we construct model PEO-LiTFSI systems with different degree of ion-pairing by tuning solvent polarity, and examine the relation between the cation-anion distinct conductivity $sigma^rm{d}_{+-}$ and the lifetime of ion-pairs $tau_{+-}$ using molecular dynamics simulations. It is found that there exist two distinct regimes where $sigma^rm{d}_{+-}$ scales with $1/tau_{+-}$ and $tau_{+-}$ respectively, and the latter is a signature of longer-lived ion-pairs which contribute negatively to the total ionic conductivity. This suggests that ion-pairs are kinetically different depending on the solvent polarity, which renders the ion-pair lifetime highly important when discussing its effect on ion transport in polymer electrolyte systems.
We investigate electric current structures in Short Large-Amplitude Magnetic Structures (SLAMS) in the terrestrial ion foreshock region observed by the Magnetospheric Multiscale mission. The structures with intense currents (|J|~1 {mu}A/m^2) have scale lengths comparable to the local ion inertial length (di). One current structure type is a current sheet due to the magnetic field rotation of the SLAMS, and a subset of these current sheets can exhibit reconnection features including the electron outflow jet and X-line-type magnetic topology. The di-scale current sheet near the edge of a SLAMS propagates much more slowly than the overall SLAMS, suggesting that it may result from compression. The current structures also exist as magnetosonic whistler waves with fci < f < flh, where fci and flh are the ion cyclotron frequency and the lower-hybrid frequency, respectively. The field rotations in the current sheets and whistler waves generate comparable |J| and energy conversion rates. Electron heating is clearly observed in one whistler packet embedded in a larger-scale current sheet of the SLAMS, where the parallel electric field and the curvature drift opposite to the electric field energize electrons. The results give insight about the thin current structure generation and energy conversion at thin current structures in the shock transition region.
A novel integral equations approach is applied for studying ion pairing in the restricted primitive model (RPM) electrolyte, i. e., the three point extension (TPE) to the Ornstein-Zernike integral equations. In the TPE approach, the three-particle correlation functions $g^{[3]}({bf r}_{1},{bf r}_{2},{bf r}_{3})$ are obtained. The TPE results are compared to molecular dynamics (MD) simulations and other theories. Good agreement between TPE and MD is observed for a wide range of parameters, particularly where standard integral equations theories fail, i. e., low salt concentration and high ionic valence. Our results support the formation of ion pairs and aligned ion complexes.
Inspired by recent experimental observations of anomalously large decay lengths in concentrated electrolytes, we revisit the Restricted Primitive Model (RPM) for an aqueous electrolyte. We investigate the asymptotic decay lengths of the one-body ionic density profiles for the RPM in contact with a planar electrode using classical Density Functional Theory (DFT), and compare these with the decay lengths of the corresponding two-body correlation functions in bulk systems, obtained in previous Integral Equation Theory (IET) studies. Extensive Molecular Dynamics (MD) simulations are employed to complement the DFT and IET predictions. Our DFT calculations incorporate electrostatic interactions between the ions using three different (existing) approaches: one based on the simplest mean field treatment of Coulomb interactions (MFC), whilst the other two employ the Mean Spherical approximation (MSA). The MSAc invokes only the MSA bulk direct correlation function whereas the MSAu also incorporates the MSA bulk internal energy. Although MSAu yields profiles that agree best with MD simulations in the near field, in the far field we observe that the decay lengths are consistent between IET, MSAc, and MD simulations, whereas those from MFC and MSAu deviate significantly. Using DFT we calculated the solvation force, which relates directly to surface force experiments. We find that its decay length is neither qualitatively nor quantitatively close to the large decay lengths measured in experiments and conclude that the latter cannot be accounted for by the primitive model. The anomalously large decay lengths found in surface force measurements require an explanation that lies beyond primitive models.
The non-Newtonian behavior of a monodisperse concentrated dispersion of spherical particles was investigated using a direct numerical simulation method, that takes into account hydrodynamic interactions and thermal fluctuations accurately. Simulations were performed under steady shear flow with periodic boundary conditions in the three directions. The apparent shear viscosity of the dispersions was calculated at volume fractions ranging from 0.31 to 0.56. Shear-thinning behavior was clearly observed at high volume fractions. The low- and high-limiting viscosities were then estimated from the apparent viscosity by fitting these data into a semi-empirical formula. Furthermore, the short-time motions were examined for Brownian particles fluctuating in concentrated dispersions, for which the fluid inertia plays an important role. The mean square displacement was monitored in the vorticity direction at several different Peclet numbers and volume fractions so that the particle diffusion coefficient is determined from the long-time behavior of the mean square displacement. Finally, the relationship between the non-Newtonian viscosity of the dispersions and the structural relaxation of the dispersed Brownian particles is examined.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
