skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Theory and simulation of electrokinetic fluctuations in electrolyte solutions at the mesoscale
Electrolyte solutions play an important role in energy storage devices, whose performance relies heavily on the electrokinetic processes at sub-micron scales. Although fluctuations and stochastic features become more critical at small scales, the long-range Coulomb interactions pose a particular challenge for both theoretical analysis and simulation of fluid systems with fluctuating hydrodynamic and electrostatic interactions. Here, we present a theoretical framework based on the Landau–Lifshitz theory to derive closed-form expressions for fluctuation correlations in electrolyte solutions, indicating significantly different decorrelation processes of ionic concentration fluctuations from hydrodynamic fluctuations, which provides insights for understanding transport phenomena of coupled fluctuating hydrodynamics and electrokinetics. Furthermore, we simulate fluctuating electrokinetic systems using both molecular dynamics (MD) with explicit ions and mesoscopic charged dissipative particle dynamics (cDPD) with semi-implicit ions, from which we identify that the spatial probability density functions of local charge density follow a gamma distribution at sub-nanometre scale (i.e. $$0.3\,{\rm nm}$$ ) and converge to a Gaussian distribution above nanometre scales (i.e. $$1.55\,{\rm nm}$$ ), indicating the existence of a lower limit of length scale for mesoscale models using Gaussian fluctuations. The temporal correlation functions of both hydrodynamic and electrokinetic fluctuations are computed from all-atom MD and mesoscale cDPD simulations, showing good agreement with the theoretical predictions based on the linearized fluctuating hydrodynamics theory.  more » « less
Award ID(s):
2103967
PAR ID:
10328419
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
942
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract A fundamental principle of chaotic quantum dynamics is that local subsystems eventually approach a thermal equilibrium state. The corresponding timescales increase with subsystem size as equilibration is limited by the hydrodynamic build-up of fluctuations on extended length scales. We perform large-scale quantum simulations that monitor particle-number fluctuations in tunable ladders of hard-core bosons and explore how the build-up of fluctuations changes as the system crosses over from integrable to fully chaotic dynamics. Our results indicate that the growth of large-scale fluctuations in chaotic, far-from-equilibrium systems is quantitatively determined by equilibrium transport coefficients, in agreement with the predictions of fluctuating hydrodynamics. This emergent hydrodynamic behaviour of subsystem fluctuations provides a test of fluctuation–dissipation relations far from equilibrium and allows the accurate determination of equilibrium transport coefficients using far-from-equilibrium quantum dynamics. 
    more » « less
  2. We report a Monte Carlo simulation study of length-scale dependent density fluctuations in cavities in the coarse-grained mW representation of water at ambient conditions. Specifically, we use a combination of test particle insertion and umbrella sampling techniques to examine the full range of water occupation states in spherical cavities up to 6.3 Å in radius in water. As has previously been observed, water density fluctuations are found to be effectively Gaussian in nature for atomic-scale cavities, but as the cavities get larger they exhibit a non-Gaussian “fat-tail” distribution for lower occupancy states. We introduce a new statistical thermodynamic approach to analyze non-Gaussian fluctuations based on the radial distribution of waters about cavities with varying numbers of waters within its boundaries. It is shown that the onset on these non-Gaussian fluctuations is a result of the formation of a bubble within the cavity as it is emptied that is accompanied by the adsorption of waters onto its interior surface. We revisit a theoretical framework we previously introduced to describe Gaussian fluctuations within cavities to now incorporate bubble formation by including surface tension contributions. This modified theory accurately describes density fluctuations within both atomic and meso-scale cavities. Moreover, the theory predicts the transition from Gaussian to non-Gaussian fluctuations at a specific cavity occupancy in excellent agreement with simulation observations. 
    more » « less
  3. Abstract The arrangement of solvent molecules and ions at solid–liquid interfaces determines electrochemical properties that are important in separations platforms, sensing technologies, and energy‐storage systems. Here we show that single glass and polymer pores in contact with propylene carbonate (PC) solutions of LiClO4exhibit an effective surface potential that is modulated by the enantiomeric excess of the solvent. In particular, electrochemical and electrokinetic measurements of ionic transport through glass pipettes and polymer pores reveal that the effective surface potential is significantly lower in solutions prepared using enantiomerically pure PC than in solutions prepared using racemic PC. Both pore systems became positively charged in all racemic solutions examined in the range of LiClO4concentrations between 1 mM and 100 mM, whereas solutions in (R)‐(+)‐PC induced a positive surface potential only at concentrations above ~5 mM. The effective surface potential is quantified through asymmetry in current–voltage curves and zeta‐potential measurements. Vibrational sum‐frequency‐generation experiments on LiClO4solutions in racemic and enantiomerically pure PC indicate that the surface lipid‐bilayer‐like region in the former is more strongly organized than in the latter, dictating the favorable positions for lithium and perchlorate ions in each case. The more ordered molecular packing in the racemic liquid leads to accumulation of lithium ions on the outside of the bilayer, creating a higher effective positive charge. Our results highlight the extreme sensitivity of the interfacial potential on molecular organization of the solvent, and the relatively unexplored role that chirality can play in electrokinetic phenomena. 
    more » « less
  4. We develop a theoretical framework to quantify how active forces renormalize the effective bending rigidity, Gaussian modulus, and surface tension of thermally fluctuating membranes. Building on classical statistical mechanics, we extend the analysis to include nonequilibrium active forces, both direct forces and those coupled to membrane curvature, within a nonlinear continuum formulation. Our model also incorporates hydrodynamic interactions mediated by the surrounding viscous fluid, which significantly alter the fluctuation spectrum. We find that direct active forces enhance long-wavelength undulations, leading to a substantial reduction in both the effective bending rigidity and surface tension, with the extent of softening strongly modulated by fluid viscosity. In contrast, curvature-coupled active forces primarily influence intermediate and short-wavelength fluctuations and show minimal sensitivity to viscosity. Together, these findings provide key insights into the nonequilibrium mechanics of active membranes and yield testable predictions for interpreting fluctuation spectra in both biological contexts and engineered membrane systems. 
    more » « less
  5. Abstract We study solar wind turbulence anisotropy in the inertial and energy-containing ranges in the inbound and outbound directions during encounters 1–9 by the Parker Solar Probe (PSP) for distances between ∼21 and 65R. Using the Adhikari et al. approach, we derive theoretical equations to calculate the ratio between the 2D and slab fluctuating magnetic energy, fluctuating kinetic energy, and the outward/inward Elsässer energy in the inertial range. For this, in the energy-containing range, we assume a wavenumberk−1power law. In the inertial range, for the magnetic field fluctuations and the outward/inward Elsässer energy, we consider that (i) both 2D and slab fluctuations follow a power law ofk−5/3, and (ii) the 2D and slab fluctuations follow the power laws withk−5/3andk−3/2, respectively. For the velocity fluctuations, we assume that both the 2D and slab components follow ak−3/2power law. We compare the theoretical results of the variance anisotropy in the inertial range with the derived observational values measured by PSP, and find that the energy density of 2D fluctuations is larger than that of the slab fluctuations. The theoretical variance anisotropy in the inertial range relating to thek−5/3andk−3/2power laws between 2D and slab turbulence exhibits a smaller value in comparison to assuming the same power lawk−5/3between 2D and slab turbulence. Finally, the observed turbulence energy measured by PSP in the energy-containing range is found to be similar to the theoretical result of a nearly incompressible/slab turbulence description. 
    more » « less