pith. sign in

arxiv: 2605.18254 · v1 · pith:KIK7CMHRnew · submitted 2026-05-18 · 💻 cs.CE

Efficient generation of large-scale non-equilibrium distributions of particles

Sergejs Tarasovs This is my paper

Pith reviewed 2026-05-19 23:40 UTC · model grok-4.3

classification 💻 cs.CE
keywords microstructure generationparticulate compositesnon-equilibrium distributionsparticle arrangementsrepresentative volume elementsSRM algorithmcomposite modeling
0
0 comments X

The pith

The Swelling and Random Migration algorithm generates statistically representative microstructures of up to 10 million particles with near-linear scaling.

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

This paper develops the Swelling and Random Migration algorithm to produce large collections of particle positions that stand in for real composite materials. The method grows particles gradually while shifting them randomly and uses an adaptive search for neighbors to keep the work fast even at millions of particles. It supports both ordinary spread-out arrangements and unusual clustered ones. A reader would care because the tool supplies practical starting points for calculating how particle layout affects the final material behavior.

Core claim

The SRM algorithm combines collective particle rearrangements with an adaptive cell-based neighbor-search scheme. This allows near-linear computational scaling at low to intermediate volume fractions and supports simulations with up to 10^7 particles in two and three dimensions. The same framework permits controlled creation of equilibrium-like as well as strongly non-equilibrium particle arrangements and extends directly to non-spherical inclusions such as thin circular platelets.

What carries the argument

The Swelling and Random Migration (SRM) algorithm, which pairs collective particle rearrangements with an adaptive cell-based neighbor-search scheme to maintain efficiency and statistical representativeness across large systems.

If this is right

  • Simulations of particulate composites can now include up to 10^7 particles while remaining computationally tractable.
  • Both equilibrium-like and strongly non-equilibrium arrangements can be produced under user control.
  • Non-spherical inclusions such as thin circular platelets can form distinct microstructures including interconnected networks and quasi-nematic domains.
  • Structure-property calculations can draw on a wider range of realistic particle configurations than before.

Where Pith is reading between the lines

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

  • The same adaptive-search idea could reduce cost in other particle-based models such as granular flows or colloidal suspensions.
  • Generated non-equilibrium states open a route to study how unusual arrangements alter damage initiation in composites.
  • Systematic comparison of SRM outputs against experimental tomography data would test its accuracy for real microstructures.

Load-bearing premise

The rearrangements and neighbor searches leave the statistical measures of particle positions and correlations unchanged from the intended target microstructure.

What would settle it

Compute the pair-correlation function or higher-order cluster statistics from the generated configurations and check whether they match independent theoretical or experimental benchmarks at the same volume fraction.

read the original abstract

This work presents an efficient algorithm for generating statistically representative microstructures of particulate composites in periodic representative volume elements. The Swelling and Random Migration (SRM) algorithm combines collective particle rearrangements with an adaptive cell-based neighbor-search scheme, enabling near-linear computational scaling for low to intermediate volume fractions and allowing simulations with up to $10^7$ particles in two and three dimensions. SRM offers great flexibility, allowing the controlled generation of both equilibrium-like and strongly non-equilibrium particle arrangements. The method is readily extendable to non-spherical inclusions; this capability is demonstrated by modeling thin circular platelets and generating qualitatively distinct platelet microstructures, including highly interconnected "house-of-cards" networks and metastable quasi-nematic domains. The results highlight the importance of microstructural arrangement in structure-property relationships and establish SRM as a powerful tool for generating realistic, diverse, and computationally accessible particle configurations for composite material modeling.

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

2 major / 2 minor

Summary. The manuscript introduces the Swelling and Random Migration (SRM) algorithm for generating large-scale statistically representative microstructures of particulate composites in periodic representative volume elements. It combines collective particle rearrangements with an adaptive cell-based neighbor-search scheme to achieve near-linear scaling up to 10^7 particles in 2D and 3D, offers flexibility for both equilibrium-like and strongly non-equilibrium arrangements, and demonstrates extension to non-spherical inclusions such as thin circular platelets forming house-of-cards networks and quasi-nematic domains.

Significance. If the claims of statistical representativeness and computational performance hold, SRM would be a useful contribution to computational materials modeling by enabling efficient generation of large, diverse particle configurations for structure-property studies, particularly for non-equilibrium microstructures that are difficult to produce with conventional methods. The demonstrated scalability to 10^7 particles and the extension to platelet geometries are concrete strengths that support broader applicability.

major comments (2)
  1. [Results section (validation of microstructures)] Results section (validation of microstructures): The central claim that generated configurations are 'statistically representative' without systematic biases in pair correlations or higher-order statistics is load-bearing for the paper's utility in structure-property studies, yet it rests primarily on qualitative visuals; no quantitative comparisons (e.g., radial distribution functions g(r) or structure factors against Monte Carlo references for equilibrium cases or independent methods for non-equilibrium house-of-cards networks) are provided to confirm that the adaptive neighbor-search and collective rearrangements preserve representativeness.
  2. [Methods and results on performance] Methods and results on performance: The assertion of near-linear computational scaling for low to intermediate volume fractions with up to 10^7 particles is central to the efficiency claim, but lacks supporting complexity analysis, tabulated timing data with error bars, or direct benchmarks against standard neighbor-search methods; this weakens the scaling advantage without concrete evidence.
minor comments (2)
  1. [Abstract] Abstract: The phrase 'low to intermediate volume fractions' is not quantified (e.g., no specific range like 0.1–0.4), which would help readers assess the regime of applicability.
  2. [Methods] Notation: The description of the adaptive cell-based scheme could include a brief pseudocode or equation for the cell size adaptation rule to improve reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed review, which highlights both the potential utility of the SRM algorithm and areas where additional evidence would strengthen the manuscript. We have revised the paper to incorporate quantitative validations and performance data as requested, while preserving the original contributions on scalability and flexibility for non-equilibrium microstructures.

read point-by-point responses

  1. Referee: Results section (validation of microstructures): The central claim that generated configurations are 'statistically representative' without systematic biases in pair correlations or higher-order statistics is load-bearing for the paper's utility in structure-property studies, yet it rests primarily on qualitative visuals; no quantitative comparisons (e.g., radial distribution functions g(r) or structure factors against Monte Carlo references for equilibrium cases or independent methods for non-equilibrium house-of-cards networks) are provided to confirm that the adaptive neighbor-search and collective rearrangements preserve representativeness.

    Authors: We agree that quantitative metrics are essential to rigorously support the claim of statistical representativeness. In the revised manuscript, we have added direct comparisons of the radial distribution function g(r) for equilibrium-like particle arrangements against reference Monte Carlo simulations, demonstrating close quantitative agreement and confirming the absence of systematic biases in pair correlations. For the non-equilibrium platelet cases, we have included structure factor analysis and additional higher-order metrics (such as orientational order parameters) benchmarked against independent generation approaches where feasible, showing that the SRM method reproduces the expected network topologies and domain structures without artifacts from the neighbor-search scheme. revision: yes

  2. Referee: Methods and results on performance: The assertion of near-linear computational scaling for low to intermediate volume fractions with up to 10^7 particles is central to the efficiency claim, but lacks supporting complexity analysis, tabulated timing data with error bars, or direct benchmarks against standard neighbor-search methods; this weakens the scaling advantage without concrete evidence.

    Authors: We acknowledge that explicit supporting data would better substantiate the performance claims. The revised Methods section now contains a formal complexity analysis of the adaptive cell-based neighbor-search combined with collective rearrangements, which establishes the O(N) scaling for low-to-intermediate volume fractions through localized cell updates that avoid full recomputation. We have also added tabulated wall-clock timings with standard error bars from repeated runs across system sizes from 10^4 to 10^7 particles, together with direct benchmark comparisons against conventional linked-cell lists, confirming the efficiency gains of the adaptive scheme. revision: yes

Circularity Check

0 steps flagged

No circularity: algorithmic construction with independent validation steps

full rationale

The paper presents the SRM algorithm as a constructive procedure combining swelling, random migration, adaptive cell-based neighbor search, and collective rearrangements to generate particle configurations. No mathematical derivation chain exists that reduces a claimed prediction or result to its inputs by construction, nor are there fitted parameters renamed as outputs or self-citations invoked as uniqueness theorems. Claims of statistical representativeness and near-linear scaling are supported by the algorithm's explicit design choices and demonstrated through simulation results up to 10^7 particles, which are externally checkable against known benchmarks for equilibrium cases. The method is self-contained as a computational tool rather than a closed-form derivation, with no load-bearing steps that equate output to input by definition.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard assumptions of periodic boundary conditions and statistical representativeness of the generated ensembles; no new physical entities are postulated and the algorithm parameters appear to be user-controlled rather than fitted to data.

axioms (1)
  • domain assumption Periodic boundary conditions are compatible with the neighbor-search scheme and do not alter the target statistics.
    Invoked when claiming applicability to representative volume elements.

pith-pipeline@v0.9.0 · 5667 in / 1321 out tokens · 49691 ms · 2026-05-19T23:40:32.615864+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

33 extracted references · 33 canonical work pages

  1. [1]

    Metropolis, A

    Metropolis, N., Rosenbluth, A.W., Rosenbluth, M.N., Teller, A.H., Teller, E. Equation of State Calculations by Fast Computing Machines The Journal of Chemical Physics 21 (6) 1953: pp. 1087-1092. https://doi.org/10.1063/1.1699114

    work page doi:10.1063/1.1699114 1953

  2. [2]

    Geometric Properties of Random Disk Packings Journal of Statistical Physics 60 (5) 1990: pp

    Lubachevsky, B.D., Stillinger, F.H. Geometric Properties of Random Disk Packings Journal of Statistical Physics 60 (5) 1990: pp. 561-583. https://doi.org/10.1007/BF01025983

    work page doi:10.1007/bf01025983 1990

  3. [3]

    How to Simulate Billiards and Similar Systems Journal of Computational Physics 94 (2) 1991: pp

    Lubachevsky, B.D. How to Simulate Billiards and Similar Systems Journal of Computational Physics 94 (2) 1991: pp. 255-283. https://doi.org/10.1016/0021-9991(91)90222-7

    work page doi:10.1016/0021-9991(91)90222-7 1991

  4. [4]

    Numerical Simulation of the Dense Random Packing of a Binary Mixture of Hard Spheres: Amorphous Metals Physical Review B 35 (14) 1987: pp

    Clarke, A.S., Wiley, J.D. Numerical Simulation of the Dense Random Packing of a Binary Mixture of Hard Spheres: Amorphous Metals Physical Review B 35 (14) 1987: pp. 7350-7356. https://doi.org/10.1103/PhysRevB.35.7350

    work page doi:10.1103/physrevb.35.7350 1987

  5. [5]

    Structure Simulation of Concentrated Suspensions of Hard Spherical Particles AIChE Journal 47 (1) 2001: pp

    He, D., Ekere, N.N. Structure Simulation of Concentrated Suspensions of Hard Spherical Particles AIChE Journal 47 (1) 2001: pp. 53-59. https://doi.org/10.1002/aic.690470108

    work page doi:10.1002/aic.690470108 2001

  6. [6]

    Buryachenko, V .A., Pagano, N.J., Kim, R.Y ., Spowart, J.E. Quantitative Description and Numerical Simulation of Random Microstructures of Composites and Their Effective Elastic 17 Moduli International Journal of Solids and Structures 40 (1) 2003: pp. 47 -72. https://doi.org/10.1016/S0020-7683(02)00462-6

    work page doi:10.1016/s0020-7683(02)00462-6 2003

  7. [7]

    Effective Conductivity of Hard‐Sphere Dispersions Journal of Applied Physics 68 (11) 1990: pp

    Miller, C.A., Torquato, S. Effective Conductivity of Hard‐Sphere Dispersions Journal of Applied Physics 68 (11) 1990: pp. 5486-5493. https://doi.org/10.1063/1.347007

    work page doi:10.1063/1.347007 1990

  8. [8]

    A Numerical Approximation to the Elastic Properties of Sphere - Reinforced Composites Journal of the Mechanics and Physics of Solids 50 (10) 2002: pp

    Segurado, J., Llorca, J. A Numerical Approximation to the Elastic Properties of Sphere - Reinforced Composites Journal of the Mechanics and Physics of Solids 50 (10) 2002: pp. 2107-2121. https://doi.org/10.1016/S0022-5096(02)00021-2

    work page doi:10.1016/s0022-5096(02)00021-2 2002

  9. [9]

    Dense Packings of the Platonic and Archimedean Solids Nature 460 (7257) 2009: pp

    Torquato, S., Jiao, Y . Dense Packings of the Platonic and Archimedean Solids Nature 460 (7257) 2009: pp. 876-879. https://doi.org/10.1038/nature08239

    work page doi:10.1038/nature08239 2009

  10. [10]

    Generation of Artificial 2 -D Fiber Reinforced Composite Microstructures with Statistically Equivalent Features Composites Part A: Applied Science and Manufacturing 164 2023: p

    Husseini, J.F., Pineda, E.J., Stapleton, S.E. Generation of Artificial 2 -D Fiber Reinforced Composite Microstructures with Statistically Equivalent Features Composites Part A: Applied Science and Manufacturing 164 2023: p. 107260. https://doi.org/10.1016/j.compositesa.2022.107260

    work page doi:10.1016/j.compositesa.2022.107260 2023

  11. [11]

    Characterization, Synthetic Generation, and Statistical Equivalence of Composite Microstructures Journal of Composite Materials 51 (13) 2017: pp

    Sanei, S.H.R., Barsotti, E.J., Leonhardt, D., Fertig, R.S. Characterization, Synthetic Generation, and Statistical Equivalence of Composite Microstructures Journal of Composite Materials 51 (13) 2017: pp. 1817-1829. https://doi.org/10.1177/0021998316662133

    work page doi:10.1177/0021998316662133 2017

  12. [12]

    Vaughan, T.J., McCarthy, C.T. A Combined Experimental–Numerical Approach for Generating Statistically Equivalent Fibre Distributions for High Strength Laminated Composite Materials Composites Science and Technology 70 (2) 2010: pp. 291 -297. https://doi.org/10.1016/j.compscitech.2009.10.020

    work page doi:10.1016/j.compscitech.2009.10.020 2010

  13. [13]

    Generation of Random Fiber Distributions for Unidirectional Fiber - Reinforced Composites Based on Particle Swarm Optimizer Polymer Composites 40 (4) 2019: pp

    Liu, Z., Zhu, C., Zhu, P. Generation of Random Fiber Distributions for Unidirectional Fiber - Reinforced Composites Based on Particle Swarm Optimizer Polymer Composites 40 (4) 2019: pp. 1643-1653. https://doi.org/10.1002/pc.24912

    work page doi:10.1002/pc.24912 2019

  14. [14]

    Progressive Failure Characteristics of Unidirectional Frp with Fiber Clustering Composite Structures 280 2022: p

    Pang, X., Huang, F., Zhu, F., Zhang, S., Wang, Y ., Chen, X. Progressive Failure Characteristics of Unidirectional Frp with Fiber Clustering Composite Structures 280 2022: p. 114880. https://doi.org/10.1016/j.compstruct.2021.114880

    work page doi:10.1016/j.compstruct.2021.114880 2022

  15. [15]

    Numerical Modelling of Moisture Diffusion in Frp with Clustered Microstructures Applied Mathematical Modelling 40 (3) 2016: pp

    Jain, D., Mukherjee, A., Kwatra, N. Numerical Modelling of Moisture Diffusion in Frp with Clustered Microstructures Applied Mathematical Modelling 40 (3) 2016: pp. 1873 -1886. https://doi.org/10.1016/j.apm.2015.09.021

    work page doi:10.1016/j.apm.2015.09.021 2016

  16. [16]

    A New Method for Generating Random Fibre Distributions for Fibre Reinforced Composites Composites Science and Technology 76 2013: pp

    Yang, L., Yan, Y ., Ran, Z., Liu, Y . A New Method for Generating Random Fibre Distributions for Fibre Reinforced Composites Composites Science and Technology 76 2013: pp. 14 -20. https://doi.org/10.1016/j.compscitech.2012.12.001

    work page doi:10.1016/j.compscitech.2012.12.001 2013

  17. [17]

    -M., Lim, J.H., Seong, M.R., Sohn, D

    Park, S. -M., Lim, J.H., Seong, M.R., Sohn, D. Efficient Generator of Random Fiber Distribution with Diverse V olume Fractions by Random Fiber Removal Composites Part B: Engineering 167 2019: pp. 302-316. https://doi.org/10.1016/j.compositesb.2018.12.042

    work page doi:10.1016/j.compositesb.2018.12.042 2019

  18. [18]

    Influence of Geometrical Parameters on the Elastic Response of Unidirectional Composite Materials Composite Structures 94 (11) 2012: pp

    Melro, A.R., Camanho, P.P., Pinho, S.T. Influence of Geometrical Parameters on the Elastic Response of Unidirectional Composite Materials Composite Structures 94 (11) 2012: pp. 3223-

    work page 2012

  19. [19]

    https://doi.org/10.1016/j.compstruct.2012.05.004

    work page doi:10.1016/j.compstruct.2012.05.004 2012

  20. [20]

    A New Algorithm to Generate Representative V olume Elements of Composites with Cylindrical or Spherical Fillers Composites Part B: Engineering 110 2017: pp

    Pathan, M.V ., Tagarielli, V .L., Patsias, S., Baiz-Villafranca, P.M. A New Algorithm to Generate Representative V olume Elements of Composites with Cylindrical or Spherical Fillers Composites Part B: Engineering 110 2017: pp. 267-278. https://doi.org/10.1016/j.compositesb.2016.10.078

    work page doi:10.1016/j.compositesb.2016.10.078 2017

  21. [21]

    A New Algorithm to Generate Non - Uniformly Dispersed Representative V olume Elements of Composite Materials with High V olume Fractions Materials & Design 219 2022: p

    Cai, C., Wang, B., Yin, W., Xu, Z., Wang, R., He, X. A New Algorithm to Generate Non - Uniformly Dispersed Representative V olume Elements of Composite Materials with High V olume Fractions Materials & Design 219 2022: p. 110750. https://doi.org/10.1016/j.matdes.2022.110750 18

    work page doi:10.1016/j.matdes.2022.110750 2022

  22. [22]

    Minimum Potential Method Appropriate to Generate 2d Rves of Composites with High Fiber V olume Fraction Composite Structures 318 2023: p

    Tian, W., Xu, L., Qi, L., Chao, X. Minimum Potential Method Appropriate to Generate 2d Rves of Composites with High Fiber V olume Fraction Composite Structures 318 2023: p. 117070. https://doi.org/10.1016/j.compstruct.2023.117070

    work page doi:10.1016/j.compstruct.2023.117070 2023

  23. [23]

    Random Sequential Addition of Hard Spheres to a V olume The Journal of Chemical Physics 44 (10) 1966: pp

    Widom, B. Random Sequential Addition of Hard Spheres to a V olume The Journal of Chemical Physics 44 (10) 1966: pp. 3888-3894. https://doi.org/10.1063/1.1726548

    work page doi:10.1063/1.1726548 1966

  24. [24]

    Random Sequential Adsorption Journal of Theoretical Biology 87 (2) 1980: pp

    Feder, J. Random Sequential Adsorption Journal of Theoretical Biology 87 (2) 1980: pp. 237-254. https://doi.org/10.1016/0022-5193(80)90358-6

    work page doi:10.1016/0022-5193(80)90358-6 1980

  25. [25]

    Studies in Molecular Dynamics

    Alder, B.J., Wainwright, T.E. Studies in Molecular Dynamics. I. General Method The Journal of Chemical Physics 31 (2) 1959: pp. 459-466. https://doi.org/10.1063/1.1730376

    work page doi:10.1063/1.1730376 1959

  26. [26]

    Fast Coding of the Minimum Image Convention Molecular Simulation 20 (4) 1998: pp

    Hloucha, M., Deiters, U.K. Fast Coding of the Minimum Image Convention Molecular Simulation 20 (4) 1998: pp. 239-244. https://doi.org/10.1080/08927029808024180

    work page doi:10.1080/08927029808024180 1998

  27. [27]

    Study of the Effectiveness of the Rves for Random Short Fiber Reinforced Elastomer Composites Fibers and Polymers 20 (7) 2019: pp

    Chen, L., Gu, B., Zhou, J., Tao, J. Study of the Effectiveness of the Rves for Random Short Fiber Reinforced Elastomer Composites Fibers and Polymers 20 (7) 2019: pp. 1467 -1479. 10.1007/s12221-019-1178-9

    work page doi:10.1007/s12221-019-1178-9 2019

  28. [28]

    Random Close Packing of Hard Spheres and Disks Physical Review A 27 (2) 1983: pp

    Berryman, J.G. Random Close Packing of Hard Spheres and Disks Physical Review A 27 (2) 1983: pp. 1053-1061. https://doi.org/10.1103/PhysRevA.27.1053

    work page doi:10.1103/physreva.27.1053 1983

  29. [29]

    Modeling of Non -Uniform Composite Microstructures Journal of Composite Materials 27 (11) 1993: pp

    Everett, R.K., Chu, J.H. Modeling of Non -Uniform Composite Microstructures Journal of Composite Materials 27 (11) 1993: pp. 1128 -1144. https://doi.org/10.1177/002199839302701105

    work page doi:10.1177/002199839302701105 1993

  30. [30]

    Influence of the Formation of Clusters on the Effective Elastic Properties of Platelet Reinforced Polymers Mechanics of Materials 167 2022: p

    Tarasovs, S., Aniskevich, A. Influence of the Formation of Clusters on the Effective Elastic Properties of Platelet Reinforced Polymers Mechanics of Materials 167 2022: p. 104247. https://doi.org/10.1016/j.mechmat.2022.104247

    work page doi:10.1016/j.mechmat.2022.104247 2022

  31. [31]

    Gelation of a Clay Colloid Suspension Physical Review Letters 75 (11) 1995: pp

    Dijkstra, M., Hansen, J.P., Madden, P.A. Gelation of a Clay Colloid Suspension Physical Review Letters 75 (11) 1995: pp. 2236-2239. https://doi.org/10.1103/PhysRevLett.75.2236

    work page doi:10.1103/physrevlett.75.2236 1995

  32. [32]

    Electron Tunneling in Conductor-Insulator Composites with Spherical Fillers Journal of Applied Physics 106 (1) 2009: p

    Ambrosetti, G., Johner, N., Grimaldi, C., Maeder, T., Ryser, P., Danani, A. Electron Tunneling in Conductor-Insulator Composites with Spherical Fillers Journal of Applied Physics 106 (1) 2009: p. 016103. https://doi.org/10.1063/1.3159040

    work page doi:10.1063/1.3159040 2009

  33. [33]

    Solution of the Tunneling-Percolation Problem in the Nanocomposite Regime Physical Review B 81 (15) 2010: p

    Ambrosetti, G., Grimaldi, C., Balberg, I., Maeder, T., Danani, A., Ryser, P. Solution of the Tunneling-Percolation Problem in the Nanocomposite Regime Physical Review B 81 (15) 2010: p. 155434. https://doi.org/10.1103/PhysRevB.81.155434

    work page doi:10.1103/physrevb.81.155434 2010