Pith. sign in

REVIEW 3 major objections 5 minor 41 references

Polarization and dynamic phases of aligning active matter in periodic obstacle arrays

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

Pith's one-line read Periodic obstacle arrays leave the polarization transition threshold unchanged but lock the polarized flow to the lattice symmetry directions, and anisotropic arrays split it into coupled and then independent one-dimensional lanes.

desk verdict Likely real new effect—substrate geometry steering collective polarization in aligning active matter—but the lane-state phase diagram needs finite-size checks before the coupled-lane phase is fully convincing. read the letter →

arxiv 2411.16882 v1 pith:TVY4YYVE submitted 2024-11-25 cond-mat.soft

classification cond-mat.soft
keywords activematterpolarizationtransitionobstaclearraysdirectionallockinglaneformationaligningtorquesBrownianparticlesflocking
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

This paper uses numerical simulations of self-propelled particles with steric repulsion and torque-based alignment to ask whether a periodic obstacle substrate changes how a flock forms. It finds that the threshold for the polarization transition is essentially unchanged by obstacles, but the direction of the polarized flow becomes locked to the symmetry axes of a dense square array. When the array is made anisotropic, global two-dimensional polarization gives way first to a coupled lane state with dense and sparse lanes, then to a decoupled state in which each lane polarizes independently left or right. The result matters because it shows that passive geometry can control where active flows go, which could be used to guide collective motion in biological or synthetic active matter.

What carries the argument

The central mechanism is the competition between two time scales: the polarization time $\tau_P \sim \mu L^2/(\beta N \sigma^2 v_0)$, the time collisions need to align orientations, and the persistence time $\tau_D = 1/D$, the time noise decorrelates orientation. The transition occurs when $\tau_P < \tau_D$, giving $\beta_c \sim \mu D/(v_0 f)$ with $f = N\sigma^2/L^2$. The obstacle array does not enter this balance, which explains why the threshold is unchanged; instead it breaks translational symmetry, locking the polarization direction and, for anisotropic arrays, confining motion to channels whose connectivity determines whether lanes couple or act independently.

What would settle it

Repeat the simulations in larger boxes (e.g., $L = 100\sigma$ or $200\sigma$) and with more obstacle rows (16 or 32) at the same filling fraction, and measure the velocity-interface probability and the locations of the quasi-isotropic, coupled, and uncoupled boundaries. If the nearly 0.2 and 0.5 interface probabilities or the phase boundaries shift systematically with system size, the identification of these as generic substrate-controlled phases would be undermined.

Watch

Extended reading notes

Core claim

In the model studied here, active Brownian particles interacting through short-range repulsion and an aligning torque $\beta(\hat{n}_i \times F_i)\cdot\hat{z}$, the transition to a polarized state at critical $\beta/D$ or $\beta N$ is not shifted by adding obstacles. Instead, for a dense enough square array the polarization vector and individual velocities become locked to the substrate's $0^\circ$, $90^\circ$, $180^\circ$, and $270^\circ$ symmetry directions. For rectangular arrays with fixed 8 rows in $y$, increasing the number of columns produces three regimes: a quasi-isotropic globally polarized state, a coupled lane state in which a few dense lanes flow in one direction while sparse lanes coexist and neighboring lanes usually share polarity (interface probability about 0.2), and a decoupled lane state in which particles cannot hop between lanes, each lane flows at uniform density, and the interface probability rises to about 0.5, indicating random per-lane left/right choice.

Load-bearing premise

The phase diagram and the lane-state interface probabilities are established in a single simulation box of side $50\sigma$ with exactly eight obstacle rows in the $y$ direction, and the paper reports no systematic variation of box size or row number, so the claim that these constitute generic substrate-controlled phases rests on this small fixed geometry representing an infinite array.

Editorial extensions

If this is right

  • Below the MIPS threshold, the polarization transition (near $\beta/D = 25$ or $\beta N = 70$) is essentially unaffected by adding obstacles, so the substrate does not alter the onset of collective order, only its direction.
  • For square arrays with $8 \times 8$ or more obstacles, both individual particle velocities and the system's mean polarization lock to $0^\circ$, $90^\circ$, $180^\circ$, and $270^\circ$; the locking sharpens as obstacle density rises.
  • For anisotropic arrays, increasing the $x$-direction obstacle count from about 15 to 25 turns the globally polarized quasi-isotropic state into a coupled lane state with a small number of dense lanes and interface probability near 0.2, and above about 25 obstacles per row into uncoupled lanes with random per-lane left/right polarity and interface probability near 0.5.
  • In the decoupled lane state, the transient time to reach steady flow decreases with increasing particle number $N$ and angular mobility $\beta$, because more collisions help lanes escape transient jamming.
  • A triangular substrate should analogously lock polarization to multiples of $60^\circ$, as the paper states.

Reading between the lines

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

  • Because the obstacle array only selects the direction of an already-formed polarization rather than changing its onset, the natural design principle that follows is to use the array geometry as a passive switch for collective flow: one could route a flock along desired axes without strengthening the alignment interaction.
  • The coupled-to-uncoupled transition in the simulations coincides with the obstacle gap becoming smaller than the particle diameter, so the lane-state boundaries should be controlled by this geometric ratio; testing this in other steric models or experiments would turn the phase diagram into a predictive engineering rule.
  • In the decoupled lane state each lane's left/right choice appears independent, so the 2D system effectively becomes an array of coupled quasi-1D aligning systems; varying the anisotropy tunes the interlane coupling continuously, which should make it possible to study the 1D-to-2D crossover in flocking order in one experimental setup.
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

3 major / 5 minor

Summary. The manuscript numerically studies a two-dimensional system of active Brownian particles with short-range steric repulsion and aligning torques, in the presence of periodic square and rectangular obstacle arrays. In the obstacle-free case the authors report a polarization transition near β/D ≈ 25, with a data collapse suggesting βc ∼ μD/(v0f). For square obstacle arrays they find that the polarization transition persists at essentially the same parameters, and that dense arrays lock the polarization direction to the substrate symmetry axes (0°, 90°, 180°, 270°). For anisotropic arrays, obtained by increasing the number of obstacles along x while keeping Np,y = 8 rows, they identify a crossover from a quasi-isotropic globally polarized state to coupled one-dimensional lanes and then to uncoupled lanes, characterized by the mean square velocities, the probability of velocity interfaces between neighboring lanes, transient times, and heat maps in the (N, Np,x) plane.

Significance. If the central claims hold, the paper establishes a new and potentially useful mechanism for controlling the direction and spatial organization of collective active flows: passive obstacle geometry can lock polarization to substrate symmetry directions and, for anisotropic arrays, induce lane states with controlled inter-lane coupling. The manuscript's strengths are the use of multiple independent order parameters (individual and collective polarization histograms, mean square velocity components, interface probabilities, transient times), the large number of independent realizations (up to 360 in Fig. 4), and the inclusion of error bars in the main order-parameter plots. The scaling collapse in Fig. 2 over D and N is a genuine attempt to make the transition estimate quantitative. However, the lane-state phase diagram is established for a single box size with exactly eight obstacle rows, which is a load-bearing gap that must be addressed before the phase diagram can be regarded as generic rather than a finite-size artifact.

major comments (3)
  1. [Section III.C, Figs. 5, 6, and 8] The identification of the coupled-lane and decoupled-lane states as distinct phases rests entirely on simulations with exactly Np,y = 8 obstacle rows, so the interface probability in Fig. 6(b) is averaged over only seven nearest-neighbor lane pairs. The reported coupled-lane value of about 0.2 and decoupled value of about 0.5 therefore correspond to roughly one or two and three or four interfaces, respectively, a small integer range in which finite-size effects could easily mimic a systematic crossover. Since the heat maps in Fig. 8 vary N and Np,x while holding Np,y = 8, they cannot establish that the phase boundaries or the interface probability are independent of the number of lanes. I request simulations at larger system sizes and at least Np,y = 12 and 16 (or, alternatively, a scaling argument showing that eight rows are sufficient) before the lane-state phase diagram can be accepted as a generic substrate-controlled phase diagram.
  2. [Section III.C, Fig. 7(a)] The transient-time measurements in Fig. 7(a) include only lanes whose particle density is at least 80% of N/8, a post hoc cutoff that is not varied or justified. Because the number of dense lanes changes with N and Np,x, this filter can systematically alter the reported τ and the conclusion that higher N reduces jamming; the text also defines the steady state only verbally as 'having a lane collectively move to the left or right.' Please report the number of excluded lanes for each parameter set and show that the trends in τ survive reasonable variations of the density cutoff and of the steady-state criterion.
  3. [Section III.A, Eq. (2) and Fig. 3] The claim that the obstacle array leaves the polarization transition threshold unchanged is supported only by Fig. 3, which shows curves for a single case (N = 240, D = 0.01) without a quantitative estimate of βc or error bars on the transition location. The scaling estimate βc ∼ μD/(v0f) is an order-of-magnitude argument, and the collapse in Fig. 2 is not a finite-size scaling analysis. I ask for a quantitative measure of the transition (e.g., susceptibility peak or Binder cumulant crossing) as a function of obstacle density, or a clear statement that the observed shift is within the resolution of the present measurements.
minor comments (5)
  1. [Section III.B, Fig. 4] The sentence 'the individual particle velocities in Fig. 4(h) exhibit small additional peaks' should refer to Fig. 4(g); Fig. 4(h) is the histogram of the mean polarization direction, not the individual particle velocities.
  2. [Section III.C and Fig. 8 caption] The schematic phase diagram is first introduced as Fig. 5(d), but the text and the Fig. 8 caption refer to it as Fig. 6(d). Please correct the cross-reference.
  3. [Section III.A] The filling fraction f = Nσ²/L² is used in the scaling estimate without being explicitly defined at first use; please define it when it is introduced.
  4. [Throughout] The terms 'coupled/decoupled' and 'coupled/uncoupled' lane states are used interchangeably (e.g., 'decoupled lane state' in the text versus 'uncoupled lanes' in Fig. 5 and the abstract); please choose one consistent nomenclature.
  5. [References] Reference [30] displays 'Soft Matter20, 2310' without a space between the journal name and the volume number; please fix the formatting.

Circularity Check

0 steps flagged · score 0.0 of 10

No material circularity: the phase diagram and scaling collapse are direct simulation measurements, and the self-citations are contextual rather than load-bearing.

full rationale

The paper's central claims are (i) that the polarization transition occurs at essentially the same alignment parameter with and without obstacles, (ii) that dense square obstacle arrays lock the polarization to the substrate symmetry directions, and (iii) that anisotropic arrays produce quasi-isotropic, coupled-lane, and decoupled-lane states. Each of these is established by direct simulation measurements reported in Figs. 3-8, not by construction from an assumed result. The scaling estimate beta_c ~ mu D/(v0 f) in Section III.A is an order-of-magnitude argument derived from the model equations, and the subsequent data collapse in Fig. 2 is a test of that scaling rather than a fit of the threshold. The lane-state order parameters (mean-square velocities, interface probability, number of velocity interfaces) are measured independently and then used to construct the schematic phase diagram, so the phase boundaries are not equal to the inputs by definition. Self-citations appear in two places: Ref. [30] is cited for the equations of motion, and Refs. [7,15,16] are cited for directional locking of run-and-tumble particles. Both are contextual or comparative; neither supplies a uniqueness theorem, a fitted parameter, or an ansatz to which the present conclusions reduce. The one notable caveat, the fixed box size L=50 sigma with Np,y=8 obstacle rows and no systematic variation in the number of lanes, is a finite-size or robustness concern rather than a circularity in the derivation chain. No circular step meets the evidentiary bar of reducing a prediction to an input by construction.

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

The central claims rest on the chosen active particle model, the WCA collision rule, and the specific periodic geometry. There are no fitted free parameters in the ledger: model parameters (beta, D, N, v0, mu, L) are inputs, and the measured thresholds are outputs rather than fit constants.

assumptions (4)
  • domain assumption Overdamped Langevin dynamics with the aligning torque (n_i cross F_i) dot z captures the essential physics of aligning active matter.
    Eqs. (1)-(2) define the model; every simulation result depends on this choice.
  • domain assumption Purely repulsive WCA interactions between particles and obstacles are sufficient; hydrodynamic, adhesive, and nonlocal interactions are neglected.
    The paper states the interaction is modeled by a short-range WCA potential in Section II.
  • domain assumption A single periodic box of side 50 sigma with Np,y=8 fixed rows is representative of the infinite periodic array.
    The phase diagram and lane statistics are computed in this geometry; no finite-size or row-number scaling is reported.
  • domain assumption The polarization transition can be characterized by the average center-of-mass velocity and its scaling with beta, D, and N, without finite-size scaling analysis.
    Figs. 2 and 3 use average velocity, but the paper does not establish thermodynamic-limit behavior for the transition.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Polarization and dynamic phases of aligning active matter in periodic obstacle arrays." pith.science (2026). https://pith.science/paper/TVY4YYVE

@misc{pith2026241116882,
  author       = {Pith},
  title        = {Pith review of: Polarization and dynamic phases of aligning active matter in periodic obstacle arrays},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TVY4YYVE}},
  note         = {Machine review of arXiv:2411.16882}
}
read the original abstract

We numerically examine a system of monodisperse self-propelled particles interacting with each other via simple steric forces and aligning torques moving through a periodic array of obstacles. Without obstacles, this system shows a transition to a polarized or aligned state for critical alignment parameters. In the presence of obstacles, there is still a polarization transition, but for dense enough arrays, the polarization is locked to the symmetry directions of the substrate. When the obstacle array is made anisotropic, at low densities the particles can form a quasi-isotropic state where the system can be polarized in any of the dominant symmetry directions. For intermediate anisotropy, the particles self-organize into a coherent lane state with one-dimensional polarization. In this phase, a small number of highly packed lanes are adjacent to less dense lanes that have the same polarization, but lanes further away can have the opposite polarization, so that global polarization is lost. For the highest anisotropy, hopping between lanes is suppressed, and the system forms uniformly dense uncoupled but polarized lanes.

Figures

Figures reproduced from arXiv: 2411.16882 by the authors.

Figure 1
Figure 1. FIG. 1. Images of the particle locations and polarization [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. (a) The system polarization [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 4
Figure 4. FIG. 4. (a,c,e,g) Normalized histograms of the distribu [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figures from the paper (4 more)
Figure 5
Figure 5. Figure 5: FIG. 5. (a,b,c) Images of the particles and (d) a schematic phase diagram for a system with an anisotropic obstacle array [PITH_FULL_IMAGE:figures/full_fig_p005_5.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Heat maps as a function of [PITH_FULL_IMAGE:figures/full_fig_p006_8.png]
Figure 6
Figure 6. Figure 6: FIG. 6. (a) The mean square velocities [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]
Figure 7
Figure 7. Figure 7: (b) for a system with N = 528. At larger values of β, the lanes are better able to overcome transient jam￾ming due to the high angular mobility of the particles. This suggests that enhancing the collective alignment of the particles in the lane, either by increasing β …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

41 extracted references · 31 canonical work pages

  1. [1]

    M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha, Hydrody- namics of soft active matter, Rev. Mod. Phys. 85, 1143 (2013)

  2. [2]

    Bechinger, R

    C. Bechinger, R. Di Leonardo, H. L¨ owen, C. Reichhardt, G. Volpe, and G. Volpe, Active particles in complex and crowded environments, Rev. Mod. Phys. 88, 045006 (2016)

  3. [3]

    Fily and M

    Y. Fily and M. C. Marchetti, Athermal phase separation of self-propelled particles with no alignment, Physical Re- view Letters 108, 235702 (2012)

  4. [4]

    G. S. Redner, M. F. Hagan, and A. Baskaran, Structure and dynamics of a phase-separating active colloidal fluid, Biophysical Journal 104, 640a (2013)

  5. [5]

    Palacci, S

    J. Palacci, S. Sacanna, A. P. Steinberg, D. J. Pine, and P. M. Chaikin, Living crystals of light-activated colloidal surfers, Science 339, 936 (2013)

  6. [6]

    M. E. Cates and J. Tailleur, Motility-induced phase sep- aration, Annual Review of Condensed Matter Physics 6, 219–244 (2015)

  7. [8]

    H. E. Ribeiro, W. P. Ferreira, and F. Q. Potiguar, Trap- ping and sorting of active matter in a periodic back- ground potential, Phys. Rev. E 101, 032126 (2020)

  8. [9]

    Reichhardt and C

    C. Reichhardt and C. Reichhardt, Clogging, dynamics, and reentrant fluid for active matter on periodic sub- strates, Physical Review E 103, 062603 (2021)

Show all 41 references
  1. [10]

    C. M. Kjeldbjerg and J. F. Brady, Partitioning of active particles into porous media, Soft Matter 18, 2757 (2022)

  2. [11]

    K. J. Modica, A. K. Omar, and S. C. Takatori, Boundary design regulates the diffusion of active matter in hetero- geneous environments, Soft Matter 19, 1890 (2023). 8

  3. [12]

    Reichhardt and C

    C. Reichhardt and C. Reichhardt, Pattern formation and transport for externally driven active matter on periodic substrates (a), Europhysics Letters 142, 37001 (2023)

  4. [13]

    C. W. Chan, D. Wu, K. Qiao, K. L. Fong, Z. Yang, Y. Han, and R. Zhang, Chiral active particles are sensi- tive reporters to environmental geometry, Nature Com- munications 15, 1406 (2024)

  5. [14]

    Caprini, F

    L. Caprini, F. Cecconi, C. Maggi, and U. Marini Bet- tolo Marconi, Activity-controlled clogging and unclog- ging of microchannels, Physical Review Research 2, 043359 (2020)

  6. [15]

    Brun-Cosme-Bruny, A

    M. Brun-Cosme-Bruny, A. F¨ ortsch, W. Zimmermann, E. Bertin, P. Peyla, and S. Rafa ¨ ı, Deflection of phototac- tic microswimmers through obstacle arrays, Phys. Rev. Fluids 5, 093302 (2020)

  7. [16]

    Reichhardt and C

    C. Reichhardt and C. J. O. Reichhardt, Directional lock- ing effects for active matter particles coupled to a peri- odic substrate, Phys. Rev. E 102, 042616 (2020)

  8. [17]

    Nabil, A

    M. Nabil, A. Frankowski, A. Orosa, A. Fuller, and A. Nourhani, Modulating drift dynamics of circle swim- mers by periodic potentials, Physical Review E 105, 054610 (2022)

  9. [18]

    Pattanayak, R

    S. Pattanayak, R. Das, M. Kumar, and S. Mishra, En- hanced dynamics of active brownian particles in periodic obstacle arrays and corrugated channels, The European Physical Journal E 42, 1 (2019)

  10. [19]

    Fazelzadeh, Q

    M. Fazelzadeh, Q. Di, E. Irani, Z. Mokhtari, and S. Jabbari-Farouji, Active motion of tangentially driven polymers in periodic array of obstacles, The Journal of Chemical Physics 159 (2023)

  11. [20]

    Nayak, S

    S. Nayak, S. Das, P. Bag, T. Debnath, and P. K. Ghosh, Driven transport of active particles through arrays of symmetric obstacles, The Journal of Chemical Physics 159 (2023)

  12. [21]

    Vent´ ejou, I

    B. Vent´ ejou, I. Magniez-Papillon, E. Bertin, P. Peyla, and A. Dupont, Behavioral transition of a fish school in a crowded environment, Physical Review E 109, 064403 (2024)

  13. [22]

    M. J. Bowick, N. Fakhri, M. C. Marchetti, and S. Ra- maswamy, Symmetry, thermodynamics, and topology in active matter, Physical Review X 12, 010501 (2022)

  14. [23]

    Dauchot and V

    O. Dauchot and V. D´ emery, Dynamics of a self-propelled particle in a harmonic trap, Phys. Rev. Lett.122, 068002 (2019)

  15. [24]

    R. H. Damascena, L. R. E. Cabral, and C. C. d. S. Silva, Coexisting orbits and chaotic dynamics of a confined self- propelled particle, Phys. Rev. E 105, 064608 (2022)

  16. [25]

    R. H. Damascena and C. C. de Souza Silva, Noise-induced escape of a self-propelled particle from metastable orbits, Phys. Rev. E 108, 044605 (2023)

  17. [26]

    Giavazzi, M

    F. Giavazzi, M. Paoluzzi, M. Macchi, D. Bi, G. Scita, M. L. Manning, R. Cerbino, and M. C. Marchetti, Flock- ing transitions in confluent tissues, Soft matter 14, 3471 (2018)

  18. [27]

    Geyer, D

    D. Geyer, D. Martin, J. Tailleur, and D. Bartolo, Freezing a flock: Motility-induced phase separation in polar active liquids, Physical Review X 9, 031043 (2019)

  19. [28]

    Caprini, U

    L. Caprini, U. Marini Bettolo Marconi, and A. Puglisi, Spontaneous velocity alignment in motility-induced phase separation, Physical review letters 124, 078001 (2020)

  20. [29]

    Fazli and A

    Z. Fazli and A. Naji, Active particles with polar align- ment in ring-shaped confinement, Physical Review E 103, 022601 (2021)

  21. [30]

    Canavello, R

    D. Canavello, R. H. Damascena, L. R. E. Cabral, and C. C. de Souza Silva, Polar order, shear banding, and clustering in confined active matter, Soft Matter20, 2310 (2024)

  22. [31]

    Briand, M

    G. Briand, M. Schindler, and O. Dauchot, Spontaneously flowing crystal of self-propelled particles, Physical review letters 120, 208001 (2018)

  23. [32]

    Baconnier, D

    P. Baconnier, D. Shohat, C. H. L´ opez, C. Coulais, V. D´ emery, G. D¨ uring, and O. Dauchot, Selective and collective actuation in active solids, Nature Physics 18, 1234 (2022)

  24. [33]

    Baconnier, O

    P. Baconnier, O. Dauchot, V. D´ emery, G. D¨ uring, S. Henkes, C. Huepe, and A. Shee, Self-aligning polar active matter, arXiv preprint arXiv:2403.10151 (2024)

  25. [34]

    Briand and O

    G. Briand and O. Dauchot, Crystallization of self- propelled hard discs, Physical review letters 117, 098004 (2016)

  26. [35]

    J. R. Howse, R. A. L. Jones, A. J. Ryan, T. Gough, R. Vafabakhsh, and R. Golestanian, Self-motile colloidal particles: From directed propulsion to random walk, Physical Review Letters 99, 048102 (2007)

  27. [36]

    Volpe, I

    G. Volpe, I. Buttinoni, D. Vogt, H.-J. K¨ ummerer, and C. Bechinger, Microswimmers in patterned environ- ments, Soft Matter 7, 8810 (2011)

  28. [37]

    Drescher, J

    K. Drescher, J. Dunkel, L. H. Cisneros, S. Ganguly, and R. E. Goldstein, Fluid dynamics and noise in bacterial cell–cell and cell–surface scattering, Proceedings of the National Academy of Sciences 108, 10940 (2011), https://www.pnas.org/doi/pdf/10.1073/pnas.1019079108

  29. [38]

    Shimoyama, K

    N. Shimoyama, K. Sugawara, T. Mizuguchi, Y. Hayakawa, and M. Sano, Collective motion in a system of motile elements, Phys. Rev. Lett. 76, 3870 (1996)

  30. [39]

    Henkes, Y

    S. Henkes, Y. Fily, and M. C. Marchetti, Active jamming: Self-propelled soft particles at high density, Phys. Rev. E 84, 040301 (2011)

  31. [40]

    C. A. Weber, T. Hanke, J. Deseigne, S. L´ eonard, O. Dau- chot, E. Frey, and H. Chat´ e, Long-range ordering of vi- brated polar disks, Phys. Rev. Lett. 110, 208001 (2013)

  32. [41]

    K.-D. N. T. Lam, M. Schindler, and O. Dauchot, Self- propelled hard disks: implicit alignment and transition to collective motion, New Journal of Physics 17, 113056 (2015)

  33. [42]

    Martin-Roca, R

    J. Martin-Roca, R. Martinez, L. C. Alexander, A. L. Diez, D. G. A. L. Aarts, F. Alarcon, J. Ram ´ ırez, and C. Valeriani, Characterization of MIPS in a suspension of repulsive active Brownian particles through dynamical features, The Journal of Chemical Physics 154, 164901 (2021)

Pith tools

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