pith. sign in

arxiv: 2604.26987 · v1 · submitted 2026-04-28 · ❄️ cond-mat.mtrl-sci

The effect of Van der Waals interaction on the microstructure of EPD deposits: a simulation study

R\'emi Martin , Sandrine Duluard , C\'eline Merlet This is my paper

Pith reviewed 2026-05-07 15:40 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords electrophoretic depositionVan der Waals interactionparticle simulationmicrostructurecoatingsaggregationelectric field
0
0 comments X

The pith

Simulations show Van der Waals self-cohesion stops shaping electrophoretic deposit microstructures above a critical electric field.

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

The paper uses particle-based simulations to test how attractive forces between particles alter the structure of coatings formed by electrophoretic deposition at different electric field strengths. These forces produce distinct arrangements near the substrate and through the deposit bulk, with measurable effects on mechanical properties. At high fields, beyond a critical strength, the attractions cease to matter and structures are set solely by particles not overlapping. This identifies a practical switch between aggregation-controlled and packing-controlled deposition regimes.

Core claim

Particle-based simulations of electrophoretic deposition that include Van der Waals self-cohesion between particles produce different microstructures near the substrate and in the bulk compared with cohesion-free cases, together with altered mechanical signatures. At high applied electric fields the influence of self-cohesion vanishes beyond a critical field value, after which the resulting microstructures are governed by volume exclusion effects alone.

What carries the argument

Particle-based simulations that add Van der Waals self-cohesion to model aggregation under an applied electric field.

If this is right

  • At low and moderate fields, self-cohesion produces clustered or less dense deposits that differ measurably from cohesion-free cases.
  • Mechanical signatures of the deposits change systematically with the shift from aggregation to volume-exclusion control.
  • Above the critical field, microstructure becomes insensitive to the strength of particle attractions.
  • Process design can target the high-field regime to obtain predictable packing-dominated coatings.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Coating engineers could operate above the critical field to reduce sensitivity of final structure to small changes in particle surface chemistry.
  • The location of the critical field likely depends on particle size, shape, and material, offering a calibration target for different suspensions.
  • Extending the simulations to include polydisperse particles would show how size variation interacts with the identified transition.

Load-bearing premise

The particle simulation model with added Van der Waals self-cohesion accurately reproduces the microscopic mechanisms and mechanical outcomes of real electrophoretic deposits across the studied field range.

What would settle it

Direct experimental comparison of real deposit microstructures and mechanical properties with and without modeled cohesion, measured while ramping the electric field until the two cases converge at the predicted critical value.

Figures

Figures reproduced from arXiv: 2604.26987 by C\'eline Merlet, R\'emi Martin, Sandrine Duluard.

Figure 1
Figure 1. Figure 1: FIG. 1. Snapshots extracted from a particle-based simulation of EPD. Starting from an homogeneous sus view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Structure of the deposit at the end of the simulation showing a layered region, close to the planar view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Pair potentials and wall-colloid interaction potential, as a function of the surface-surface separa view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Pair potentials (thick lines) and first rdf peak, averaged over the bulk (thin lines). The rdf represented view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6 view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Average volume fraction in the core of the deposit (12 view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Average pore diameter in the core region of the deposits (12 view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Core-averaged tensile strengths as a function of electric field for view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Surface density of particles in the first two layers as a function of electric field. view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Layer bond surface density view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Voronoi tilings of first layers projected in the view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Average Voronoi tile surface, normalized by particle surface area view at source ↗
read the original abstract

Electrophoretic deposition is a method of choice for generating coatings thanks to its ease of implementation and its ability to produce coatings of relatively large thicknesses in a single step process. While this process also benefits from a large number of tunable parameters to adapt the coating to each application (applied electric field, particle concentration, viscosity of the suspension, etc...), such a freedom can lead to the selection of parameters being an overwhelming task. A better fundamental understanding of the microscopic phenomena and mechanisms at play during deposition can provide clues for the more efficient design of optimized coatings. Particle-based models, allowing for the simulation of deposit microstructures for various process parameters, are particularly interesting to get insights in such systems. Nevertheless, such studies are rare and usually do not involve the possibility of self-cohesion between particles, while it seems crucial for the final structure of the deposit. Here, we use particle-based simulations to study the influence of aggregation on the deposit formed for different applied electric fields. We show that the self-cohesion indeed leads to different microstructures, both in the close vicinity of the substrate and in the bulk of the deposit, and relate this to the mechanical signature of the deposits. Our results reveal that at high electric field, the influence of self-cohesion on resulting microstructures essentially vanishes beyond a critical field strength. This marks the transition between a deposition regime affected by aggregation to a regime largely dominated by volume exclusion effects.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 3 minor

Summary. The manuscript reports particle-based simulations of electrophoretic deposition (EPD) that incorporate Van der Waals self-cohesion between particles. By comparing runs with and without this term across a range of electric field strengths, the authors find that self-cohesion alters local packing near the substrate and bulk microstructure at moderate fields, with corresponding changes in mechanical response; above a critical field the differences vanish and volume-exclusion effects dominate.

Significance. If the numerical results hold, the work supplies a clear mechanistic picture of a field-driven crossover from aggregation-controlled to volume-exclusion-controlled deposition regimes. This supplies a practical handle for tuning EPD coating density and mechanical properties via field strength alone. The explicit inclusion of self-cohesion in an otherwise standard electrophoretic model is a constructive modeling choice that previous simulation studies often omit.

major comments (1)
  1. [Results] Results section: the critical field strength is introduced as the point at which self-cohesion effects 'essentially vanish,' yet no quantitative threshold (e.g., relative difference in coordination number, packing fraction, or radial distribution function falling below X %) is stated. Because the central claim is the existence and location of this transition, an explicit, reproducible definition is required.
minor comments (3)
  1. [Methods] Methods: the precise functional form, cutoff, and Hamaker constant of the added Van der Waals term are not given; these parameters directly set the magnitude of self-cohesion and therefore the location of the reported transition.
  2. [Figures and text] Figure captions and text: the mechanical signature is repeatedly invoked but never defined (e.g., whether it is obtained from a simulated uniaxial compression test, from contact-number statistics, or from another observable).
  3. [Methods and Results] The manuscript would benefit from a short statement of the number of independent runs, system size, and statistical uncertainty on the microstructure metrics shown in the figures.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive assessment of the work and for the constructive comment on the definition of the critical field. We address this point below and have prepared a revised manuscript.

read point-by-point responses

  1. Referee: [Results] Results section: the critical field strength is introduced as the point at which self-cohesion effects 'essentially vanish,' yet no quantitative threshold (e.g., relative difference in coordination number, packing fraction, or radial distribution function falling below X %) is stated. Because the central claim is the existence and location of this transition, an explicit, reproducible definition is required.

    Authors: We agree that an explicit, quantitative definition of the critical field is required to make the central claim reproducible. In the revised manuscript we now define the critical field as the lowest applied field at which the relative difference in bulk packing fraction between the cohesive and non-cohesive runs falls below 2 % and the relative difference in average coordination number falls below 5 %; both quantities are required to satisfy the criterion simultaneously. These thresholds were chosen because they correspond to the field value at which the two data sets converge within the statistical uncertainty obtained from our ensemble of independent runs. We have added a dedicated paragraph in the Results section stating this definition, updated the figure captions to mark the critical field with this criterion, and briefly justified the choice of thresholds in the text. revision: yes

Circularity Check

0 steps flagged

No significant circularity: results are direct simulation outputs

full rationale

The paper reports outcomes from particle-based simulations that compare microstructures with and without an added Van der Waals self-cohesion term under varying electric fields. The claimed transition (self-cohesion effects vanishing above a critical field, leaving volume-exclusion dominant) is an observed difference between simulation runs rather than a quantity defined in terms of itself or fitted to a subset and then re-predicted. No mathematical derivation chain, self-referential definitions, or load-bearing self-citations appear in the argument; the model inputs (forces, fields, particle interactions) are independent of the reported microstructural metrics.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review yields no explicit free parameters, axioms, or invented entities; the model implicitly assumes standard colloidal interaction potentials and a particular representation of the electric field driving force.

pith-pipeline@v0.9.0 · 5562 in / 1052 out tokens · 33996 ms · 2026-05-07T15:40:34.158340+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

29 extracted references · 29 canonical work pages

  1. [1]

    Pair potential The pair potential models inter-colloid interactions using the classical DLVO framework as the sum of a repulsive long-range screened electrostatic potentialU e, a mid-range attractive vdW potentialU vdW , and a close-range steric repulsive potentialU r: UDL V O=U r +UvdW +Ue.(1) The screened electrostatic pair potential is modeled as an ex...

    work page

  2. [2]

    External potential A similar Hamaker potential is used for the interaction between particles and the substrate by replacing the radius of one of the spheres with infinity. The potential then becomes a function of thezcoordinate of the particle: uwall r = Acc σ 6 LJ· 7560 7a−z (z−a) 7 + 7a+z (z+a) 7 ,(7) U wall vdW =− Acc 6 · 2za z2 −a 2 +ln z−a z+a .(8) 7...

    work page

  3. [3]

    Dissipative interactions The interaction with the implicit solvent is modeled through the Fast Lubrication Dynamics (FLD) algorithm[20], an approximation of Stokesian dynamics[21]. This method has been devel- oped for modelling pairwise short-range lubrication and long-range hydrodynamic interaction, as well as the associated noise, efficiently and consis...

    work page

  4. [4]

    A. R. Boccaccini and I. Zhitomirsky, Current Opinion in Solid State and Materials Science6, 251 (2002)

    work page 2002

  5. [5]

    O. O. V . der Biest and L. J. Vandeperre, Annual Review of Materials Science29, 327 (1999)

    work page 1999

  6. [6]

    Besra and M

    L. Besra and M. Liu, Progress in Materials Science52, 1 (2007)

    work page 2007

  7. [7]

    Amrollahi, J

    P. Amrollahi, J. S. Krasinski, R. Vaidyanathan, L. Tayebi, and D. Vashaee, inHandbook of Nano- electrochemistry: Electrochemical Synthesis Methods, Properties, and Characterization Techniques, edited by M. Aliofkhazraei and A. S. H. Makhlouf (Springer International Publishing, Cham, 2016) pp. 561–591

    work page 2016

  8. [8]

    Corni, M

    I. Corni, M. P. Ryan, and A. R. Boccaccini, Journal of the European Ceramic Society28, 1353 (2008)

    work page 2008

  9. [9]

    Prioux, S

    M. Prioux, S. Duluard, F. Ansart, G. Pujol, P. Gomez, and L. Pin, Journal of the American Ceramic Society103, 6724 (2020)

    work page 2020

  10. [10]

    Verma, H

    K. Verma, H. Cao, P. Mandapalli, and R. Wille, Journal of Computational Science41, 101075 (2020). 32

    work page 2020

  11. [11]

    M. A. Salazar de Troya, J. R. Morales, B. Giera, A. J. Pascall, M. A. Worsley, R. Landingham, W. L. Du Frane, and J. D. Kuntz, Materials & Design209, 110000 (2021)

    work page 2021

  12. [12]

    Giera, L

    B. Giera, L. A. Zepeda-Ruiz, A. Luis, A. J. Pascall, and T. H. Weisgraber, Langmuir33, 652 (2017)

    work page 2017

  13. [13]

    Derjaguin and L

    B. Derjaguin and L. Landau, Acta Physicochim. URSS14, 633 (1941)

    work page 1941

  14. [14]

    Verwey and G

    E. Verwey and G. Overbeek,Theory of the stability of lyophobic colloids(Elsevier, 1948)

    work page 1948

  15. [15]

    J. N. Israelachvili,Intermolecular and Surface F orces, 3rd ed. (Academic Press, 2011)

    work page 2011

  16. [16]

    J. J. Karnes, A. P. Pascall, C. Rehbock, V . Ramesh, M. A. Worsley, S. Barcikowski, E. Lee, and B. Giera, Computer Physics Communications297, 109062 (2024)

    work page 2024

  17. [17]

    Sarkar and P

    P. Sarkar and P. S. Nicholson, Journal of the American Ceramic Society79, 1987 (1996)

    work page 1987

  18. [18]

    J. D. Weeks, D. Chandler, and H. C. Andersen, The Journal of Chemical Physics54, 5237 (1971)

    work page 1971

  19. [19]

    Filion, R

    L. Filion, R. Ni, D. Frenkel, and M. Dijkstra, The Journal of Chemical Physics134, 134901 (2011)

    work page 2011

  20. [20]

    Giera, L

    B. Giera, L. A. Zepeda-Ruiz, A. J. Pascall, J. D. Kuntz, C. M. Spadaccini, and T. H. Weisgraber, Journal of The Electrochemical Society162, D3030 (2015)

    work page 2015

  21. [21]

    Elimelech, J

    M. Elimelech, J. Gregory, X. Jia, and R. williams,Particle Deposition and Aggregation: Measure- ment, Modelling, and Simulation(Butterworth-Heinemann, 1998)

    work page 1998

  22. [22]

    Everaers and M

    R. Everaers and M. R. Ejtehadi, Phys. Rev. E67, 041710 (2003)

    work page 2003

  23. [23]

    Kumar,Microscale dynamics in suspensions of non-spherical particles, Ph.D

    A. Kumar,Microscale dynamics in suspensions of non-spherical particles, Ph.D. thesis, University of Illinois at Urbana-Champaign (2010)

    work page 2010

  24. [24]

    J. F. Brady and G. Bossis, Annual Review of Fluid Mechanics20, 111 (1988)

    work page 1988

  25. [25]

    Kumar and J

    A. Kumar and J. L. Higdon, Physical Review E82, 051401 (2010)

    work page 2010

  26. [26]

    B. J. Alder and T. E. Wainwright, The Journal of Chemical Physics27, 1208 (1957)

    work page 1957

  27. [27]

    S.-E. Phan, W. B. Russel, Z. Cheng, J. Zhu, P. M. Chaikin, J. H. Dunsmuir, and R. H. Ottewill, Physical Review E54, 6633 (1996)

    work page 1996

  28. [28]

    Parisi and F

    G. Parisi and F. Zamponi, The Journal of Chemical Physics123, 144501 (2005)

    work page 2005

  29. [29]

    G. L. Hunter and E. R. Weeks, Reports on Progress in Physics75, 066501 (2012). 33

    work page 2012