{"id":"58348047-deac-4962-a9c5-5654af82410e","arxiv_id":"2411.17264","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In a driven symmetric electrolyte, ion density correlations form a universal conical long-range shape and relax diffusively after a field quench, explaining the algebraic relaxation of the current.","lead":"Using stochastic density functional theory, this paper shows that ions in an electrolyte driven by an electric field develop long-range correlations shaped like a cone, and that these correlations spread diffusively after the field is switched on. The result connects the microscopic correlation shape to the previously observed slow, algebraic relaxation of electric current after a field quench.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Gaussian closure is the load-bearing step; the dropped ∇·(n∇V*n) nonlinearity could renormalize the anisotropic diffusion in Eq. (15), so the cone angle is not yet established beyond RPA.","rationale":"The paper's internal logic is coherent: Eq. (12) is exact for the linearized model, Eqs. (15)–(19) follow from its small-s expansion, and the numerical inversions in Figs. 1, 2, 4, and 5 are consistent with the stated formulas. The reader's conditional verdict is appropriate. My stress-test focuses on the one unguarded physical input: the Gaussian/RPA closure of Eq. (9). This is not an internal inconsistency—the calculations are correct within the model—but a question of whether the model predicts the real electrolyte. Because the asymptotic tail is generated by the k→0 behavior, even a small nonlinear correction to the anisotropic diffusion constant changes the singular part and the cone angle, not just subleading terms. The absence of a small parameter and the absence of an external benchmark make this concern real. The proposed simulation check would settle it; until then, conditional acceptance is right. I recommend no change to the reader's verdict.","tokens_in":11718,"tokens_out":20706,"duration_ms":190606,"concrete_test":"Run Brownian-dynamics simulations of a 1:1 symmetric electrolyte at low packing fraction (e.g., λ_D/a ≈ 10, f = 1) in a periodic slab under a uniform field, and measure c_{++}(r) and c_{+-}(r) for separations 5–100 λ_D. Test the predicted scaling collapse of Fig. 2 and the zero-contour cone angle of Eq. (19). If the collapse fails or the angle shifts when the screening length is varied at fixed f, the Gaussian closure is not controlling the tail; if the collapse holds across two decades in λ_D, the concern is settled.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Every asymptotic statement in Sections 3 and 4 uses the linearized Eq. (9), obtained by dropping the ∇·(n∇V*n) nonlinearity from the current in Eq. (2). No small parameter is given that would control this truncation. Since the long-range behavior is controlled by k→0, the dropped term contributes at one loop a self-energy of order κ∫d^d q V(q)C(q). With V(q)∝q^{-2} and the Gaussian C(q) of order unity in the small-q region (Eq. 14), this integral is not parametrically small; in d=3 it is cutoff-sensitive. If nonzero, it renormalizes the anisotropic coefficient in the singular part c_sing = s^2/[2(s^2+f^2 s_parallel^2)] (Eq. 15), hence the cone angle Θ_d (Eq. 19). The paper contains no loop expansion or error estimate showing these corrections are subdominant at large distances. Thus the conical universality and the diffusive-front scaling are properties of the Gaussian model, not yet proven properties of the original stochastic dynamics.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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].","tokens_in":11874,"tokens_out":21914,"duration_ms":201818,"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":[{"comment":"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.","section":"Section 3, Eq. (9)"},{"comment":"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.","section":"Section 4.2, Eqs. (29), (31), (33)"}],"minor_comments":[{"comment":"There are several typographical errors, including \"fluctuacting hydrodynamics\", \"hovewer\", \"swich off\", \"Nondimensionlizing\", and \"algebraiclly\". These should be corrected in a final pass.","section":"Introduction"},{"comment":"The sentence ending \"...and is referred to as the hydrodynamic correction. It.\" contains a stray fragment \"It.\" that should be removed or completed.","section":"Section 2, after Eq. (8)"},{"comment":"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.","section":"Section 5.2, Eqs. (49)-(50)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the central derivation is coherent. The main risk is not internal consistency but the external validity of the RPA/linearization at large distances; the one-loop estimate in my report suggests this is a substantive concern that should be addressed by the authors, even if only by clearly delimiting the claims as RPA-level predictions. The sign errors in Section 4.2 are straightforward to fix, but they currently affect explicit equations used in the discussion. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Plainly: this paper contains a genuinely new result. It derives, within stochastic density functional theory and the RPA closure, that steady-state ion correlations in a driven symmetric electrolyte have a long-range conical shape, not the parabolic shape seen for short-range interactions. The singular Fourier part, the finite cone angle as the field goes to zero, and the diffusive transient scaling are new and internally consistent. I checked the directional limits, the zero of g_d, and the interpolation of the transient expression; they hold together. There are no fitted parameters, and the connection to the earlier algebraic current relaxation is a consistency check, not a circular import.\n\nThe main soft spot is the linearized SDFT/RPA closure itself. Equation (9) drops the nonlinear term ∇·(n∇V*n), and no small parameter controls that truncation. In d=3, a one-loop self-energy from that term is not parametrically small, and it could renormalize the anisotropic diffusion in Eq. (15), hence the cone angle. So the conical shape and transient front are rigorously properties of the Gaussian model, not yet proven for the full stochastic dynamics. The abstract and conclusion call the shape universal, which is overstated: the calculation is for a binary symmetric electrolyte inside RPA. A referee should ask for a one-loop estimate or simulation support. This is a genuine caveat, but not a fatal flaw.\n\nMinor: a few algebraic inversions are stated without derivation, which makes checking harder. That is a presentation issue.\n\nWho this is for: people working on electrolytes, driven soft matter, or fluctuating hydrodynamics. It deserves serious peer review. I would cite it, and I would bring it to a reading group.","headline":"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.","tokens_in":704,"tokens_out":2081,"would_cite":true,"duration_ms":36856,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["charged fluids","driven diffusive systems","fluctuating hydrodynamics","transport properties","correlations","stochastic density functional theory","non-equilibrium steady state","electrolytes"],"falsifier":"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$.","tokens_in":11462,"feed_emoji":"⚡","tokens_out":13572,"duration_ms":110392,"temperature":0.7,"pith_summary":"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.","feed_headline":"Driven electrolytes make ion correlations form a cone","feed_subtitle":"After a field quench the cone spreads by diffusion, explaining the algebraic current relaxation.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"supplies the exact stochastic density equation from which the linearized fluctuation dynamics is derived.","marker":"[7]"},{"why":"provides the SDFT formalism and the stationary-state correlation solution from which the singular part is extracted.","marker":"[9]"},{"why":"first points out long-ranged NESS correlations in driven electrolytes, motivating the shape analysis.","marker":"[14]"},{"why":"gives the treatment of hydrodynamic interactions in SDFT used to include solvent flow in the model.","marker":"[17]"},{"why":"supplies the parabolic correlation shape for driven short-range interacting particles that the conical shape is contrasted with.","marker":"[21]"},{"why":"reports the same parabolic long-range correlations in another short-range driven system, reinforcing the contrast class.","marker":"[22]"},{"why":"reports the algebraic relaxation of the total charge current that this paper reproduces and explains from diffusive correlations.","marker":"[23]"}],"fun_headline_variants":["Universal cone shape emerges in driven electrolytes","Quench in electric field spreads correlations diffusively","Algebraic current relaxation links to cone-shaped correlations","Driven electrolytes: correlations form a universal cone","Ion correlations relax diffusively after a field quench"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Universal cone shape emerges in driven electrolytes","Quench in electric field spreads correlations diffusively","Algebraic current relaxation links to cone-shaped correlations","Driven electrolytes: correlations form a universal cone","Ion correlations relax diffusively after a field quench"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000232,"raw_usage":{"total_tokens":1462,"prompt_tokens":887,"completion_tokens":575,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":503,"completion_tokens_details":{"reasoning_tokens":504}},"tokens_in":503,"tokens_out":575,"duration_ms":6152,"temperature":1.0,"reasoning_tokens":504,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:22:21.869015+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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$.","supporting_citations":[{"cited_title":"At NESS, the density- density correlation functions are anisotropic and decay algebraically with distance","cited_arxiv_id":null,"evidence_quote":"supplies the exact stochastic density equation from which the linearized fluctuation dynamics is derived."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the SDFT formalism and the stationary-state correlation solution from which the singular part is extracted."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"first points out long-ranged NESS correlations in driven electrolytes, motivating the shape analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the treatment of hydrodynamic interactions in SDFT used to include solvent flow in the model."},{"cited_title":"Correlation decoupling of Casimir interaction in an electrolyte driven by external electric fields","cited_arxiv_id":"2404.06028","evidence_quote":"reports the same parabolic long-range correlations in another short-range driven system, reinforcing the contrast class."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"reports the algebraic relaxation of the total charge current that this paper reproduces and explains from diffusive correlations."}],"review_version":1}