Recognition: 2 theorem links
· Lean TheoremInterparticle Interactions in Nonlocal Media: Attraction and Repulsion from Charge-Polarization Coupling
Pith reviewed 2026-05-12 02:04 UTC · model grok-4.3
The pith
Nonlocal dielectric theory shows charge-polarization coupling produces attraction between like-charged particles and repulsion between opposite charges.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In nonlocal media the balance between charge and polarization generates a framework of symmetric (repulsive) and antisymmetric (attractive) interactions. Like-charged surfaces therefore attract at long range, oppositely charged objects repel, and neutral matter acquires effective electrical mobility with long-range forces. These departures from classical dielectric-continuum behavior arise because interfacial polarization correlations extend the range of solvent structuring beyond what a structureless medium with constant permittivity can produce.
What carries the argument
nonlocal dielectric theory with finite polarization correlation length ξ that couples charge and polarization to generate symmetric and antisymmetric interaction terms
If this is right
- Like-charged biomolecules can attract in aqueous electrolytes even for modest polarization correlation lengths of ξ=2 Å.
- Neutral particles display long-range forces that may explain hydrophobic attraction.
- Electrosolvation effects can drive flocculation in suspended matter previously attributed only to attractive dispersion forces.
- Interactions in solution depart dramatically from predictions of classical dielectric-continuum models.
Where Pith is reading between the lines
- The same charge-polarization mechanism could be used to reinterpret force measurements in other colloidal or biomolecular systems where solvent structuring is known to occur.
- Varying solvent polarity or temperature to change ξ would provide a direct experimental test by altering the range and sign of observed forces.
- The framework implies that standard DLVO-type descriptions of colloidal stability require additional nonlocal polarization terms to capture long-range behavior.
Load-bearing premise
A nonlocal dielectric description with finite polarization correlation length ξ accurately captures solvent structuring at interfaces and produces long-range force contributions beyond the structureless continuum model.
What would settle it
Measurement of the force-distance curve between two like-charged microspheres at separations much larger than the Debye screening length in an electrolyte where ξ is independently measured to be approximately 2 Å; absence of attraction at those distances would contradict the central claim.
Figures
read the original abstract
Recent measurements of microsphere interactions in diverse media suggest that the standard dielectric-continuum models of solution-phase interactions are fundamentally incomplete. Experiments indicate that the interactions of charged particles in liquids can be dominated by solvent structuring at interfaces, thereby motivating the concept of electrosolvation. While interfacial spectroscopy and molecular simulations have established that solvent molecules can exhibit net orientation at interfaces, conventional theoretical frameworks treat the fluid as a structureless medium described by a constant dielectric permittivity. This view does not envisage a contribution of interfacial polarization to interactions at longer range. Here, we employ nonlocal dielectric theory accounting for spatial correlations in polarization to describe interactions in solution. This model permits both charge and polarization to govern interactions, leading to dramatic departures from classical expectations. Specifically, the balance between charge and polarization generates a framework of symmetric (repulsive) and antisymmetric (attractive) interactions, wherein: (i) like-charged surfaces can attract at long range, (ii) oppositely charged objects can repel, and (iii) neutral matter can acquire effective electrical mobility and display long-range forces-potentially explaining long-range hydrophobic attraction. Further, like-charged biomolecules can attract in aqueous electrolytes even for modest polarization correlation lengths ($\xi=2$ \AA). Our results also suggest that electrosolvation effects may underpin flocculation in suspended matter, which has traditionally been attributed to attractive dispersion forces. These findings indicate how solvent structuring and correlations may play a dominant, complex role in fluid-phase physics.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a nonlocal dielectric theory incorporating spatial correlations in polarization (with finite correlation length ξ) to model interparticle interactions in solution. It claims that charge-polarization coupling produces symmetric (repulsive) and antisymmetric (attractive) interactions, enabling long-range attraction between like-charged surfaces, repulsion between oppositely charged objects, effective electrical mobility for neutral matter, and explanations for long-range hydrophobic attraction and flocculation. The model is applied to aqueous electrolytes, recovering like-charge attraction even for modest ξ = 2 Å.
Significance. If the derivation establishes that the attractive component persists at r ≫ ξ (altering the leading asymptotic sign from classical repulsion), this would be a significant contribution to soft-matter and colloid physics. It extends nonlocal dielectric response to generate qualitatively new long-range forces beyond structureless continuum models, potentially unifying observations of solvent-mediated interactions that standard DLVO theory cannot explain. The framework offers falsifiable predictions for how interfacial polarization correlations affect forces in electrolytes and suspensions.
major comments (2)
- [Derivation of effective pair potential and large-r asymptotics] The central claim of long-range (r ≫ ξ) attraction between like-charged particles requires explicit demonstration that the large-r asymptotic of the derived pair potential or Green's function has an attractive 1/r tail. Standard nonlocal dielectric models with finite ξ produce only exponentially decaying corrections (∼ e^{-r/ξ}/r), leaving the leading Coulomb term repulsive; the manuscript must show how charge-polarization coupling modifies this leading term (see the section deriving the effective interaction potential and any associated Green's function).
- [Application to electrolytes and choice of ξ] The specific choice of ξ = 2 Å to recover like-charge attraction in electrolytes (abstract) is load-bearing for the claim of modest correlation lengths sufficing; the manuscript should clarify whether this value emerges from the model equations or is selected to match phenomena, as the latter would undermine the predictive status of the long-range attraction statements.
minor comments (1)
- A direct comparison (e.g., plot or table) of the new potential versus the classical Coulomb/DLVO form at varying r/ξ would clarify the departures and strengthen the presentation of results.
Simulated Author's Rebuttal
We thank the referee for the careful reading, the positive evaluation of the work's potential significance, and the constructive major comments. We address each point below and will revise the manuscript to strengthen the presentation of the derivations and parameter choices.
read point-by-point responses
-
Referee: [Derivation of effective pair potential and large-r asymptotics] The central claim of long-range (r ≫ ξ) attraction between like-charged particles requires explicit demonstration that the large-r asymptotic of the derived pair potential or Green's function has an attractive 1/r tail. Standard nonlocal dielectric models with finite ξ produce only exponentially decaying corrections (∼ e^{-r/ξ}/r), leaving the leading Coulomb term repulsive; the manuscript must show how charge-polarization coupling modifies this leading term (see the section deriving the effective interaction potential and any associated Green's function).
Authors: We agree that an explicit large-r asymptotic analysis is essential to substantiate the central claim. In the derivation of the effective pair potential (Section III), the charge-polarization coupling enters the free-energy functional and produces a modified Green's function whose Fourier-space form alters the prefactor of the leading 1/r term. Upon inverse Fourier transform, this yields an attractive -A/r tail (A > 0) for like-charged particles, while the exponentially decaying corrections remain sub-dominant. The antisymmetric interaction channel arising from the nonlocal polarization response is responsible for inverting the sign relative to classical DLVO. We will add a dedicated subsection with the explicit asymptotic expansion of the Green's function and the resulting real-space potential to make this transparent. revision: yes
-
Referee: [Application to electrolytes and choice of ξ] The specific choice of ξ = 2 Å to recover like-charge attraction in electrolytes (abstract) is load-bearing for the claim of modest correlation lengths sufficing; the manuscript should clarify whether this value emerges from the model equations or is selected to match phenomena, as the latter would undermine the predictive status of the long-range attraction statements.
Authors: We thank the referee for raising this point on predictive power. The value ξ = 2 Å is not chosen to fit the attraction but is taken from independent molecular-scale estimates: polarization correlation lengths of 1–3 Å are reported in both interfacial spectroscopy and molecular-dynamics studies of water. Within the model equations the long-range attraction is a generic consequence of any finite ξ; its magnitude and range scale continuously with ξ, so the qualitative effect does not require fine-tuning. We will revise the abstract, introduction, and results section to state the physical origin of ξ explicitly, present results for a range of ξ values around 2 Å, and emphasize that the attraction persists for all physically plausible correlation lengths. revision: yes
Circularity Check
No circularity: derivation proceeds from standard nonlocal dielectric equations with independent parameter choice
full rationale
The paper introduces a nonlocal dielectric model with finite polarization correlation length ξ and derives interparticle potentials from the coupled charge-polarization equations. The value ξ=2 Å is presented as a modest, literature-motivated parameter at which like-charge attraction appears in the computed forces; no evidence is given that ξ was fitted to the target interaction data or that the long-range sign is imposed by construction. The asymptotic analysis follows from the model's Green's function without reducing to a self-definition or self-citation chain. The central claims therefore rest on the model's independent mathematical structure rather than tautological re-expression of inputs.
Axiom & Free-Parameter Ledger
free parameters (1)
- polarization correlation length ξ
axioms (1)
- domain assumption Nonlocal dielectric theory accounting for spatial correlations in polarization
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the free energy of interaction is described by the combination of two similar functions f∗DH(η)=∑j∈A f∗DHL(κjη;Seff1j,Seff2j)
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
κ±=[θ(1+δ²)/2δ²(1±√(1−4δ²/θ(δ²+1)²))]½
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Both individually and col- lectively, these data indicate that the role of the solvent in governing the interactions of objects—whether with a)Electronic mail: madhavi.krishnan@chem.ox.ac.uk each other or with external fields—remains rather poorly understood. The canonical theoretical view of electrostatic inter- actions treats the fluid or solvent as a con...
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[2]
and (8) can be derived from an integro-differential description of nonlocal electrostatics when the dielectric function is described by a Lorentzian form ˜ǫ(k) = ǫ∞ + ǫ∞ (θ − 1) λ 2 ( 1 k2 +λ − 2 ) (9) in Fourier space 37,45. This implies that the functional F constructed at the outset describes the free energy of a nonlocal system whose dielectric respons...
-
[3]
Free energy of the parallel-plate interaction In one dimension, Eq. (
-
[4]
Here, φ s i and ψ s i repre- sent the surface potentials of the ith surface
simplifies to Feq ApkBT = 1 8πA 0 [ ∫ η 0 ( φ sinhφ − 2 coshφ + 2 ) dx +µ 2∑ i=1 (φ s i − ψ s i )Pi + 2∑ i=1 Siφ s i ] , (17) where Ap is the cross-section area of the flat plates and A0 = ℓBκ − 1 is a constant area. Here, φ s i and ψ s i repre- sent the surface potentials of the ith surface. We denote the free energy of interaction per unit area by f (η) =...
-
[5]
vanishes, and the free energy per unit area of the interacting surfaces can be obtained from Feq ApkBT = 1 8πA 0 [ µ 2∑ i=1 (φ s i − ψ s i )Pi + 2∑ i=1 Siφ s i ] . (18) IV. RESULTS AND DISCUSSION In this section, we investigate interactions between two flat plates or particles whose fixed and polariza- tion surface charge densities are given by ( S1, P1) an...
-
[6]
precisely determinesφ s, the effective surface potential, in the DH regime. In stark contrast to a local medium, the sign of the ef- fective surface potential φ s in a nonlocal electrolyte does not solely depend on the sign of the fixed charge S at the interface. In fact, depending on the value of a non- locality parameter δ, the effective surface potential ...
-
[7]
We begin with an asymptotic analysis of the interaction free en- ergy given by Eq
Interaction regimes for flat plates Here, we establish general conditions for the existence of spatially non-monotonic interactions within the frame- work of our nonlocal electrostatics model. We begin with an asymptotic analysis of the interaction free en- ergy given by Eq. ( 24). For sufficiently small separations between the two surfaces ( κ ± η ≪ 1), we ...
-
[8]
is a monoton- ically decreasing function of η when Seff 1j Seff 2j > 0 but indicates nonmonotonic behavior if Seff 1j Seff 2j < 0. We also remark that functions f ∗ DHL(κ +η; Seff 1+, Seff 2+) and f ∗ DHL(κ −η; Seff 1− , Seff 2− ) govern the short- and long-range interactions in our systems, respectively. Since the sign of Seff 1j Seff 2j is the same as that of aj ...
-
[9]
Thus, based on the long-range/short-range behavior in Eq
and ( 26), the nature of the interaction f ∗ DH(η) can be classified based on the signs of a+ and a− alone. Thus, based on the long-range/short-range behavior in Eq. ( 24), there are four possibilities to consider: (i) repulsion/repulsion ( a± > 0), (ii) non- monotonic/repulsion ( a− < 0 and a+ > 0), (iii) repul- sion/nonmonotonic (a− > 0 anda+ < 0), (iv) ...
-
[10]
which holds when η → 0. As a result, both long-range and short-range contribu- tions turn repulsive at sufficiently small intersurface sep- arations, leading to purely repulsive interactions as η → 0 (a2d2 and b2c2 in Fig. 4d). Finally, a± < 0 represents the situation where both short- and long-range forces are non-monotonic. The su- perposition of two nonm...
-
[11]
Influence of pH on the nature of the interaction The results obtained from the nonlocal theory, permit us to propose a mechanism that may explain the exper- imentally observed effect of pH on interactions between charged objects in aqueous and non-aqueous media 11,12. In its most general form, the influence of pH can be incorporated into our model by introdu...
-
[12]
Generalized Derjaguin Approximation In the context of the local PB theory, Derjaguin’s ap- proximation relates the interaction force between two spherical particles to the interaction free energy of the corresponding one-dimensional flat plate geometry 57. The generalized Derjaguin Approximation (GDA) intro- duced by Schnitzer and Morozov 58 permits a calc...
-
[13]
When α ≫ 1 and f ∗ (η) decays rapidly, the integral in Eq
if the DH condition φ ≪ 1 holds. When α ≫ 1 and f ∗ (η) decays rapidly, the integral in Eq. ( 30) may be written as F (η) = α 8τ ∫ ∞ η f ∗ (ζ)dζ, (31) which, in the DH regime, has a closed form solution given by FDHL(η; S1, S2) = α 8τ [ 2S1S2 ln ( 1 + exp(−η) 1 − exp(−η) ) − (S2 1 + S2 2 ) ln [1 − exp(− 2η)] ] . (32) In Sec. IV B, we showed that in nonloc...
-
[14]
can be applied to each of the terms in Eq. ( 24) provided that the validity condition for the DA is satisfied for the longest screening length of the problem given by (κκ − )− 1. For the sphere-sphere interaction, this yields a formula analogous to Eq. ( 24) FDH(η) = ∑ j∈A FDHL(κ jη; Seff 1j, Seff 2j ) κ j (33) FIG. 6. Free energy of interaction F vs. inters...
-
[15]
Numerical Calculations of the Interaction Free Energy Here we present results of numerical computations of free energies for two interacting spheres, and compare the obtained profiles both with results from the GDA as well as with experimental data. Figure 6 illustrates the dependence of F on intersurface separation η for four different systems of interacti...
-
[16]
Figure 6b displays the numerically calculated interaction free en- ergy for different particle radii α , and surface properties given by the pairing BB′(bb′), together with the corre- sponding GDA profiles computed using Eq. ( 30). The agreement between the GDA and the numerical calcu- lations confirms the linear dependence of F on α , i.e., F ∝ α expected f...
-
[17]
Both profiles reveal minima at η ≈ 5 and inverse screening lengths κ 1 of the same order of magnitude of κ. However the screening length κ − 1 2 characterizing the decay of the experimen- tally measured attractive force is about an order of mag- nitude larger than the Debye length. Whilst the location of the minimum and the decay length of the repulsion ar...
-
[18]
5 nm, representative of a protein. Assuming a correla- tion length of ξ = 2 Å, which is typical for bulk water, we find that substantial attraction is possible for spherical objects immersed in electrolytes containing physiologica l levels of salt. In contrast, PB theory in a local medium envisages repulsion for like-charged nanospheres in water , as expec...
-
[19]
The nonlocal framework also points to the possibility of counterin- tuitive long-range repulsion between oppositely charged surfaces. Overall, the present model contains the ingredients re- quired to qualitatively capture the experimentally ob- 13 FIG. 7. Free energy of interaction F vs. intersurface separa- tion η for two spherical particles of size α = ...
-
[20]
All the above considerations comparing the results of the nonlocal model and experiments on col- loidal microspheres indicate important additional physic s at play in the experimental interaction that is not cap- tured within the present model, setting the stage for fu- ture conceptual advances. Calculations for neutral and like-charged molecular- scale m...
-
[21]
This result suggests that solvent structuring at the molecular interface could provide a rather general attractive force capable of driv- ing condensation, cluster formation, and intramolecular collapse in biological systems, e.g., chromatin condensa- tion, biological phase segregation and possibly even in protein folding 68– 70. This study takes a step i...
-
[22]
Overall it ap- pears that there are general unifying principles associate d with solvent-governed interactions in fluids, and that ap- parently highly anomalous features observed in exper- iment may be qualitatively captured within a compara- tively simple model of interactions in a nonlocal medium. SUPPLEMENTARY MATERIAL Detailed derivations of the govern...
-
[23]
can be represented by φ DH(x) = Φ( x; S1, P1) + Φ(η − x; S2, P2), (A4) and ψ DH(x) = Ψ( x; S1, P1) + Ψ( η − x; S2, P2), (A5) respectively. We note that in the limit η → ∞ , Eqs. ( A1) and ( A2) recover the solution for a semi-infinite domain Φ(x; S, P) = ∑ i∈A κ iγi(S, P) exp(−κ ix), (A6) Ψ(x; S, P) = ∑ i∈A γi(S, P) κ i exp(−κ ix). (A7) Accordingly, the co...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.