Pith. sign in

REVIEW 2 major objections 3 minor 36 references

Stationary and transient correlations in driven electrolytes

T0 review · 2 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Under a steady electric field, all ion-ion correlations in a symmetric electrolyte acquire the same algebraic, conical long-range shape, and after a field quench they relax by anisotropic diffusion.

desk verdict A careful SDFT/RPA derivation of a genuinely new conical long-range correlation shape in driven electrolytes, with the Gaussian closure as the main caveat. read the letter →

arxiv 2411.17264 v2 pith:N4K2PRD2 submitted 2024-11-26 cond-mat.stat-mech

Haggai Bonneau , Vincent Démery , Elie Raphaël This is my paper · ORCID
classification cond-mat.stat-mech
keywords chargedfluidsdrivendiffusivesystemsfluctuatinghydrodynamicstransportpropertiescorrelationsstochasticdensityfunctionaltheorynon-equilibriumsteadystateelectrolytes
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper asks what pair correlations look like in a binary symmetric electrolyte (two ion species of opposite charge and equal mobility) driven by a uniform electric field, and how they evolve after the field is suddenly switched on or off. Working in stochastic density functional theory truncated at Gaussian order, the authors find that in the non-equilibrium steady state (NESS) all ion-ion correlation functions share the same algebraic long-range tail, whose zero contour is a cone around the field direction rather than the parabola found for short-range interacting driven fluids. After switching the field on, the correlations relax by anisotropic diffusion, with the cone appearing inside a growing length scale $\sqrt{T\kappa t}$ while equilibrium exponential correlations persist beyond it. These diffusive correlations reproduce the previously reported algebraic $t^{-3/2}$ relaxation of the total charge current on switching on, and the exponential relaxation on switching off.

What carries the argument

The carrying object is the singular part of the Fourier-space correlation, $\tilde{c}_{\rm sing}=s^2/[2(s^2+f^2 s_{\parallel}^2)]$, the term that contains the discontinuity at the wavevector origin after higher powers of $s_{\parallel}$ and $s_{\perp}$ are discarded. It is identical for all correlation elements, and its inverse Fourier transform is the Laplacian of the Green function of an anisotropic Poisson equation, which produces the conical shape. The same diffusive kernel explains the transient behavior: after a quench the correlations relax like solutions of the anisotropic diffusion equation with diffusion constant enhanced by $1+f^2$ along the field direction, and the mesoscopic number-charge equations in the long-time regime reduce to this diffusion equation.

What would settle it

Simulate the full stochastic density equations without discarding third-order fluctuations, or measure density-density correlations in a dilute driven electrolyte at separations of tens of Debye lengths: the zero contour of the correlation should be a cone with $\tan\theta=\sqrt{2/(1+f^2)}$ in three dimensions, whose angle tends to $\sin^{-1}\sqrt{2/3}\approx 54.7^\circ$ as the field goes to zero rather than flattening to $90^\circ$.

Watch

Extended reading notes

Core claim

The central discovery is that the non-equilibrium steady state of a binary symmetric electrolyte (two species, opposite charges, equal mobility) has a universal conical correlation structure. In Fourier space, the discontinuity at the origin is captured by the same singular part $\tilde{c}_{\rm sing}=s^2/[2(s^2+f^2 s_{\parallel}^2)]$ for every correlation element; in real space this is the Laplacian of the Green function of an anisotropic Poisson equation, so the correlation decays algebraically and vanishes on a cone of angle $\Theta_d=\sin^{-1}\sqrt{(d-1)/(d+f^2)}$. The cone angle tends to $\sin^{-1}\sqrt{1-1/d}$ as the field tends to zero, so the long-range anisotropy persists even for infinitesimal driving. For a quench in the field, the correlation obeys the diffusive scaling $c(x,\tau)=[1/(8\pi^{d/2}\tau^{d/2})][f^2/(1+f^2)^{3/2}]\Phi(x/\sqrt{\tau})$, and the same dynamics combined with the long-range Coulomb and hydrodynamic kernels yields the algebraic $t^{-3/2}$ relaxation of the charge current.

Load-bearing premise

The derivation keeps only Gaussian (second-order) fluctuations in the stochastic density equations; the cone's algebraic tail and angle would change if nonlinear density fluctuations or finite-size effects become relevant at the large distances where the cone is predicted.

Editorial extensions

If this is right

  • In any $d$-dimensional binary symmetric electrolyte, all correlation elements (same-charge and opposite-charge) share the same algebraic long-range tail; only the region within a Debye length distinguishes species.
  • The cone angle tends to $\sin^{-1}\sqrt{1-1/d}$ as the field goes to zero, so the long-range anisotropic correlation survives arbitrarily weak driving.
  • After a field quench, correlations obey diffusive scaling with length scale $\sqrt{T\kappa t}$: inside that scale the NESS cone is established, outside it correlations decay exponentially in space.
  • The diffusive correlation dynamics yields algebraic $t^{-3/2}$ relaxation of the total charge current after switch-on and exponential relaxation after switch-off, matching the previously reported current behavior.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the cone is universal, measuring the zero-contour angle at fixed large distance could serve as a non-invasive probe of the effective field strength $f=q\lambda_D E/T$, and deviations might expose charge renormalization.
  • The same singular Fourier structure should appear in other long-range interacting driven systems described by the same linearized field theory, such as charged colloids with screened electrostatics.
  • Because the derivation is for a bulk system, a natural test is confinement: charged walls or nanochannels should shift the cone angle by modifying the anisotropic Poisson Green function, distinguishing this mechanism from short-range parabolic correlations.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 3 minor

Summary. This paper uses stochastic density functional theory (SDFT/Dean–Kawasaki) with a linearized, Debye–Hückel/RPA closure to compute density–density correlations in a binary symmetric electrolyte after a sudden change of the applied electric field. In the non-equilibrium stationary state the authors isolate a singular Fourier part c_sing = s^2/[2(s^2+f^2 s_parallel^2)] (Eq. 15), whose inverse transform yields algebraically decaying conical correlations with zero-angle Theta_d = sin^{-1}(sqrt((d-1)/(d+f^2))) (Eq. 19). For a field quench on, they show that the transient correlations obey diffusive scaling c(x,tau) ~ tau^{-d/2} Phi(x/sqrt(tau)) (Eq. 26), with a front at length sqrt(T kappa t); for a quench off, the correlations relax to the equilibrium Yukawa form. The transient correlation dynamics is then connected to the algebraic tau^{-3/2} relaxation of the charge current previously reported in Ref. [23].

Significance. If the results hold, they provide a sharp and testable distinction between driven electrolytes and driven systems with short-range interactions: the long-range NESS correlations are conical rather than parabolic, with a universal angle that survives the zero-field limit. The paper is analytically self-contained: the singular part is extracted explicitly, the cone angle follows from a closed-form scaling function, and the tau^{-3/2} current relaxation is re-derived from the correlation dynamics rather than imported from earlier work. The main value is therefore as a concrete prediction of SDFT in the RPA regime, and as a conceptual link between transient correlations, diffusive spreading, and algebraic current relaxation. The strength of this significance is conditional on the status of the linearization, which is the central caveat discussed below.

major comments (2)
  1. [Section 3, Eq. (9)] The conical result and the diffusive-front scaling rest entirely on the linearized Dean–Kawasaki equation (9), where the nonlinear term nabla·(n nabla V*n) is discarded. No small parameter controls this truncation at the wave numbers that determine the long-range tail. In particular, the one-loop self-energy correction is of order kappa ∫ d^d q V(q) C(q); with V(q) ~ q^{-2} and the RPA C(q) of order unity for q→0 (Eq. 14), this integral is cutoff-sensitive in d=3 and not parametrically small. If nonzero, such a correction renormalizes the coefficient f^2 in the singular part (15), which would change the cone angle (19). The authors should either show that this one-loop correction leaves the form of (15) invariant, or explicitly restrict the universal-cone claim to the RPA model and state this qualification in the abstract and conclusion.
  2. [Section 4.2, Eqs. (29), (31), (33)] As written, these switch-off expressions do not satisfy the correct initial condition. At tau=0, Eq. (29) tends to -1/(2A) - f^2 u^2/(2ABC), hence to a negative value in the s→0 limit, whereas the initial state is the NESS correlation (15), which tends to +1/2. The sign of the first term should be +1/(2A); the same sign error is repeated in (31) and (33). This contradicts the statement in the text that the equilibrium value 1/2 is recovered as tau→∞, which is indeed what the c_{+-} expression (34) gives. Since these equations are used to justify the exponential relaxation after switch-off, they must be corrected.
minor comments (3)
  1. [Introduction] There are several typographical errors, including "fluctuacting hydrodynamics", "hovewer", "swich off", "Nondimensionlizing", and "algebraiclly". These should be corrected in a final pass.
  2. [Section 2, after Eq. (8)] The sentence ending "...and is referred to as the hydrodynamic correction. It." contains a stray fragment "It." that should be removed or completed.
  3. [Section 5.2, Eqs. (49)-(50)] The integrals leading from (47)-(48) to the closed-form prefactors in (49) and (50) are not shown. Since these prefactors are central to the comparison with Ref. [23], including the intermediate integration steps (or a short appendix) would improve verifiability.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the conical stationary correlations and diffusive transient are derived from the linearized SDFT equations, and the current-relaxation comparison is a consistency check rather than an imported prediction.

full rationale

The paper's derivation is self-contained from the linearized SDFT dynamics. Equation (9) is the input model; equation (12) is the exact correlation dynamics for a binary symmetric electrolyte; the stationary solution (13) is obtained by solving the algebraic steady-state equation, and the singular part (15) is the leading small-wavevector behavior of that solution, not an assumed ansatz. The real-space cone (17)-(19) follows by Fourier inversion and elementary algebra. The transient result (26) is likewise obtained by taking the long-range limit s→0 of the exact quench solution (20)-(23) and inverting the resulting diffusion kernel; no diffusive front is assumed. Section 5 re-derives the current relaxation from the same correlation dynamics: equations (47)-(50) are explicit integrals over the computed correlations, and the agreement with Ref. [23] is a backward consistency check, not a load-bearing import. References [9] and [23] are by overlapping authors, but the equations used from them are re-derived here or are framework-level model assumptions. The RPA/Gaussian closure and the quasi-stationary elimination (43) are physical approximations; their range of validity is a correctness question, not a circularity, because none of the target results—cone shape, zero-angle, diffusive scaling, or tau^{-3/2} relaxation—is fed into the derivation as an input. No fitted parameter is renamed as a prediction. Hence no circular step is identified.

Assumptions & free parameters 0 free parameters · 6 assumptions · 0 invented entities

The central derivations use no free parameters fitted to data. External field, temperature, mobility, charge, density, and dielectric constant are physical inputs; f is their dimensionless combination. The nontrivial premises are the linearized SDFT/RPA closure and the quasi-static treatment of charge fluctuations at long time. No new entities are introduced.

assumptions (6)
  • domain assumption Electrolyte density dynamics are governed by the Dean-Kawasaki/SDFT equations with multiplicative Ito noise; correlations are computed from the linearized equations after discarding third-order fluctuations (RPA/Debye-Huckel closure).
    Invoked in Eqs. (1)-(12); all correlation results, including the cone shape and transient diffusion, are derived inside this Gaussian closure.
  • domain assumption Ions are overdamped Brownian particles in an incompressible, inertialess Stokes flow; hydrodynamic interactions enter through the Oseen tensor and fluctuating stress.
    Eqs. (2), (5)-(7); this excludes fluid inertia and memory, so the long-time tails discussed are not hydrodynamic long-time tails in the usual sense.
  • domain assumption The system is bulk, electroneutral, and subjected to a spatially uniform time-dependent field E(t)=E(t)e_x; electrode and boundary effects are ignored.
    Section 2; the quench protocol relies on a spatially uniform Heaviside field, and the bulk assumption underlies the algebraic long-range analysis.
  • domain assumption The explicit calculations are restricted to a binary symmetric electrolyte (two species, opposite unit charges, equal mobility, equal density).
    After Eq. (11); the closed-form Fourier solutions and the cone-angle formula are derived only under this symmetry.
  • domain assumption At large times and lengths, charge fluctuations are quasi-static, so the temporal derivative, diffusion, and noise in the charge equation are neglected, giving Delta = -(i/2) q E.k U.
    Section 5.1, Eq. (43); this approximation is load-bearing for the derivation of the current relaxation from the correlation dynamics.
  • standard math Distributional Fourier transforms and the split of a discontinuous Fourier function into a singular and a regular part are valid; the regular part decays faster in real space.
    Section 3.1; standard asymptotics, but the uniqueness of the chosen singular part is assumed.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Stationary and transient correlations in driven electrolytes." pith.science (2026). https://pith.science/paper/N4K2PRD2

@misc{pith2026241117264,
  author       = {Pith},
  title        = {Pith review of: Stationary and transient correlations in driven electrolytes},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/N4K2PRD2}},
  note         = {Machine review of arXiv:2411.17264}
}
read the original abstract

Particle-particle correlation functions in ionic systems control many of their macroscopic properties. In this work, we use stochastic density functional theory to compute these correlations, and then we analyze their long-range behavior. In particular, we study the system's response to a rapid change (quench) in the external electric field. We show that the correlation functions relax diffusively toward the non-equilibrium stationary state and that in a stationary state, they present a universal conical shape. This shape distinguishes this system from systems with short-range interactions, where the correlations have a parabolic shape. We relate this temporal evolution of the correlations to the algebraic relaxation of the total charge current reported previously.

Figures

Figures reproduced from arXiv: 2411.17264 by the authors.

Figure 1
Figure 1. The correlation function in the non-equilibrium stationary state for different values of the normalized external field f. Panels (a,b,c) are the anion￾cation correlation c−+, and panels (d,e,f) are the equal charge correlation cαα. One can observe that far from the origin, a cone develops and that at vanishing fields, the angle of the cone reaches a finite value, as predicted in equation (19). Surprisingly, the long… view at source ↗
Figure 2
Figure 2. The rescaled correlation function Cαβ|x∥| 3f−2 along slices of constant x∥ values. This allows us to compare it to the scaling function g3. The dotted lines correspond to f = 1 while the solid lines correspond to f = 5. We can see in solid black the evaluation of equation (18). unlike, for example, the Mach cone. From equation (17), we can see that the angle Θ of the cone where the value of the correlation function … view at source ↗
Figure 3
Figure 3. Sketch of the length scales at play and the functional dependence of the density-density correlation function on the distance from the origin. The behavior in the upper row (magenta) corresponds to the spatial relaxation of the correlation function from equilibrium to NESS (switching on). The behavior in the lower row (green) corresponds to the spatial relaxation of the correlation function from NESS to equilibrium … view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: A time series of the correlation element c−+ after a sudden switching on of the external field. In the upper row, c−+ along real axes. In the lower row, the rescaled correlation function c−+τ 3/2 along rescaled axes. The figure has been evaluated for f = 1. Near the ce…
Figure 5
Figure 5. Figure 5: A time series of the correlation element c−+ after a sudden switching off of the external field. In the upper row, c−+ along real axes. In the lower row, the rescaled correlation function c−+τ 3/2 along rescaled axes. The figure has been evaluated for f = 1. Away from …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

36 extracted references · 26 canonical work pages

  1. [23]

    Péraud J P, Nonaka A J, Bell J B, Donev A and Garcia A L 2017Proceedings of the National Academy of Sciences114 10829–10833 ISSN 0027-8424 URLhttp://www.pnas.org/content/ 114/41/10829

    work page

  2. [1]

    The conductivity of electrolytes, for instance, has been explored with this approach in the early works of Debye and Hückel [1] and later Onsager [2]

    Introduction In many-particle systems with pair interactions, the transport properties, such as viscosity or conductivity, are related to the pair correlation in non-equilibrium steady states. The conductivity of electrolytes, for instance, has been explored with this approach in the early works of Debye and Hückel [1] and later Onsager [2]. There, the co...

    work page

  3. [2]

    Model We model an electrolyte as a system of charged Brownian particles of different species [9, 14, 16, 18]. The particles move in a homogeneous three-dimensional solution and are subjected to a uniform external electric field with a time-dependent amplitude E(t) = E(t)ˆex, where ˆex is the unit vector along thex axis. The particles interact via the elec...

    work page

  4. [3]

    The correlation elements in Fourier space for binary symmetric electrolytes can be calculated using equation (12)

    Correlations in electrolyte systems – Long range behavior at NESS In this section, we analyze the structure of the correlations between different ions at scales larger than the Debye length. The correlation elements in Fourier space for binary symmetric electrolytes can be calculated using equation (12). In the stationary state case, as shown in [9], the ...

    work page

  5. [4]

    The focus remains on lengths that are large compared to the Debye length and times that are large compared to the Debye time

    Time dependent correlations Now that we have characterized the correlations in the NESS, we investigate the temporal evolution of the correlations between the short-ranged correlations at equilibrium and the algebraic correlations in the NESS. The focus remains on lengths that are large compared to the Debye length and times that are large compared to the...

    work page

  6. [5]

    Mesoscopic density fields and conductivity The diffusive dynamics of the correlation suggest a relation to the algebraic relaxation of the charge current reported in [23]. To establish this relation, we start by normalizingtheexpressioninequation(8)togetthedimensionlesstotalchargecurrent J(τ ) σ0E0 = g(τ ) + 1 ¯ρλ3 D γel(τ, f) + g(τ ) rs λD γhyd(τ, f), (3...

    work page

  7. [6]

    At this stage, we can gain some insight into the relaxation rates of the conductivity corrections that we found in [23]

    It is easy to see that equation (24) gives the solution to this equation. At this stage, we can gain some insight into the relaxation rates of the conductivity corrections that we found in [23]. Stationary and transient correlations in driven electrolytes 13 5.2. Relaxation rates of the currents When we swich off the external field, the system transitions...

    work page

  8. [7]

    At NESS, the density- density correlation functions are anisotropic and decay algebraically with distance

    Conclusion In this work, we have characterized the behavior of the particle-particle correlation functions in the long-range regime, in the non-equilibrium stationary state (NESS), and in the transient regime as the system approaches NESS. At NESS, the density- density correlation functions are anisotropic and decay algebraically with distance. These prop...

    work page

Show all 36 references
  1. [8]

    Debye P and Hückel E 1923Physikalische Zeitschrift 24 305–325

  2. [9]

    Onsager L 1927Transactions of the Faraday Society23 341–349

  3. [10]

    Onsager L and Kim S K 1957The Journal of Physical Chemistry 61 198–215 URL https: //doi.org/10.1021/j150548a015

    doi

  4. [11]

    Onsager L and Fuoss R M 1932 The Journal of Physical Chemistry 36 2689–2778 URL https://doi.org/10.1021/j150341a001

    1932 doi

  5. [12]

    Brady J F and Morris J F 1997 Journal of Fluid Mechanics 348 103–139 URL https: //doi.org/10.1017/S0022112097006320

    1997 doi

  6. [13]

    Kawasaki K 1994 Physica A: Statistical Mechanics and its Applications 208 35–64 URL https://doi.org/10.1016/0378-4371(94)90533-9

    1994 doi

  7. [14]

    Dean D S 1996 Journal of Physics A: Mathematical and General29 L613–L617 URL http: //stacks.iop.org/0305-4470/29/i=24/a=001

    1996

  8. [15]

    Hansen J P and McDonald I R 2013Theory of simple liquids: with applications to soft mat- ter (Academic press) URL https://www.elsevier.com/books/theory-of-simple-liquids/ hansen/978-0-12-387032-2

  9. [16]

    Démery V and Dean D S 2016Journal of Statistical Mechanics: Theory and Experiment2016 023106 URL http://stacks.iop.org/1742-5468/2016/i=2/a=023106

    2016

  10. [17]

    Krüger M, Solon A, Démery V, Rohwer C M and Dean D S 2018The Journal of chemical physics 148 URL https://doi.org/10.1063/1.5019424

    doi

  11. [18]

    Robin P 2024 The Journal of Chemical Physics 160 064503 ISSN 0021-9606 URL https: //doi.org/10.1063/5.0188215

    2024 doi

  12. [19]

    Dean D S and Podgornik R 2014Phys. Rev. E89 032117 URL http://link.aps.org/doi/10. 1103/PhysRevE.89.032117

  13. [20]

    Dean D S, Lu B S, Maggs A C and Podgornik R 2016 Phys. Rev. Lett. 116 240602 URL http://link.aps.org/doi/10.1103/PhysRevLett.116.240602

    2016 doi

  14. [21]

    Mahdisoltani S and Golestanian R 2021Phys. Rev. Lett.126 158002 URL https://link.aps. org/doi/10.1103/PhysRevLett.126.158002

    doi

  15. [22]

    Du G, Dean D S, Miao B and Podgornik R 2024arXiv preprint arXiv:2404.06028

    arXiv

  16. [24]

    Donev A, Garcia A L, Péraud J P, Nonaka A J and Bell J B 2019 Current Opinion in Electrochemistry 13 1–10 ISSN 2451-9103 URL http://www.sciencedirect.com/science/ article/pii/S2451910318301960

    2019

  17. [25]

    Avni Y, Adar R M, Andelman D and Orland H 2022Phys. Rev. Lett. 128(9) 098002 URL https://link.aps.org/doi/10.1103/PhysRevLett.128.098002

    doi

  18. [26]

    Avni Y, Andelman D and Orland H 2022The Journal of Chemical Physics157 154502 URL https://doi.org/10.1063/5.0111645

    doi

  19. [27]

    Bernard O, Jardat M, Rotenberg B and Illien P 2023The Journal of Chemical Physics 159 164105 ISSN 0021-9606 URLhttps://doi.org/10.1063/5.0165533 Stationary and transient correlations in driven electrolytes 15

    doi

  20. [28]

    Poncet A, Bénichou O, Démery V and Oshanin G 2017 Phys. Rev. Lett. 118 118002 URL http://link.aps.org/doi/10.1103/PhysRevLett.118.118002

    2017 doi

  21. [29]

    Bain N and Bartolo D 2017Nature communications 8 15969 URL https://doi.org/10.1038/ ncomms15969

  22. [30]

    Bonneau H, Démery V and Raphaël E 2023 Journal of Statistical Mechanics: Theory and Experiment 2023 073205 URL https://dx.doi.org/10.1088/1742-5468/acdced

    2023 doi

  23. [31]

    Øksendal B 2000 Stochastic differential equations: an introduction with applications 5th ed (Springer) ISBN 9780387602431

    2000

  24. [32]

    De Zarate J M O and Sengers J V 2006Hydrodynamic fluctuations in fluids and fluid mixtures (Elsevier)

  25. [33]

    1 (Springer)

    Bockris J O, Reddy A K and Gamboa-Adelco M E 2006 Modern Electrochemistry Vol. 1 (Springer)

    2006

  26. [34]

    Onsager L 1939Journal of Physical Chemistry43 189–196

  27. [35]

    Mahdisoltani S and Golestanian R 2021 New Journal of Physics 23 073034 URL https: //dx.doi.org/10.1088/1367-2630/ac0f1a

    2021 doi

  28. [36]

    Minh THN , Rotenberg B and Marbach S 2023 Faraday Discussions 246 225–250 URL https://doi.org/10.1039/D3FD00031A

    2023 doi

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.