Pith. sign in

REVIEW 4 major objections 4 minor 39 references

Rolling, sliding and trapping of driven particles in square obstacle lattices

T0 review · 4 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read For a magnetically driven microparticle in a square lattice of cylindrical pillars, the ratio of pillar edge-to-edge spacing to particle radius determines whether it rolls forward, slides backward, or becomes trapped in vertical…

desk verdict A genuinely new experimental observation—rolling-to-sliding reversal and trapping in obstacle lattices—is supported by a plausible but not fully quantitative simulation mechanism; the geometric range inconsistency and missing error bars need fixing. read the letter →

arxiv 2505.21702 v1 pith:AFOPKF7C submitted 2025-05-27 cond-mat.soft

Galor Geva , Arin Escobar , Paula Magrinya , Pablo Llombart , Alfredo Alexander-Katz , Laura R. Arriaga , Juan L. Aragones This is my paper
classification cond-mat.soft
keywords rotatingmagneticmicroparticlessquareobstaclelatticerolling-to-slidingtransitionpressureforcesroto-translationalcouplingStokesiandynamicshydrodynamictrappingmicrofluidiccargotransport
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 establishes that the direction of travel of a magnetically driven microparticle rolling over a flat floor can be reversed by surrounding it with a square lattice of cylindrical pillars, with the pillar spacing relative to the particle radius acting as the control knob. At large spacing the particle rolls forward in the direction of its rotation because substrate shear forces dominate. At small spacing the rotational flow of the particle collides with the pillars, creating pressure gradients that overpower the shear and push the particle backward in a sliding motion. In between, the particle is hydrodynamically trapped, oscillating up and down with a frequency locked to the rotation rate. Because forward, backward, and trapped states are all selected simply by changing the geometry of the environment, the result offers a direct design rule for steering microscale cargo.

What carries the argument

The central object is the normalized roto-translational coupling parameter $\xi/\xi_0$, which compares a particle's in-lattice translation speed to its free-rolling speed on a bare substrate. The mechanism that carries the argument is the competition between shear lubrication forces from the substrate, which drive rolling, and pressure forces generated when the particle's rotational flow field is reflected by the obstacle cylinders, breaking the fore-aft symmetry of the flow and driving backward sliding. The simulations isolate the pressure contribution by holding the particle at fixed height using a sphere-based Rotne\textendash Prager\textendash Yamakawa model of a 12-bead dodecahedral roller of effective radius $R = 0.782a$, and show that the same non-monotonic velocity profile across the unit cell appears even when the particle cannot move vertically, with the strongest pressure force at the unit-cell center.

What would settle it

Reverse the magnetic field's rotation direction in a lattice with small $l_p/R$: if pressure forces drive sliding, the particle should keep sliding opposite to the rolling direction, flipping its travel direction, and the critical spacing at which rolling switches to sliding should remain unchanged. A failure of either prediction would show that something other than fore-aft pressure forces sets the direction of motion.

Watch

Extended reading notes

Core claim

The central claim is a geometry-controlled reversal of the translational direction of an externally driven rotating particle. For a freely rotating particle near a flat substrate, the measured dimensionless roto-translational coupling, $\xi_0 = v_p/(2\pi f R)$, is positive, meaning the particle rolls in the direction of its rotation. Inside a square lattice of cylindrical obstacles, the normalized coupling $\xi/\xi_0$ falls linearly as the edge-to-edge spacing $l_p/R$ shrinks; when $l_p/R$ becomes small enough that $\xi/\xi_0$ crosses zero, the particle switches from rolling to sliding opposite the rolling direction. The paper locates an intermediate trapping regime near $l_p/R \approx 4.0$\textendash$4.4$ where $\xi/\xi_0 \approx 0$ and particles execute a small vertical oscillation whose frequency scales linearly with the rotation frequency. Using simulations of a bead-resolved dodecahedral roller at fixed height, the paper shows that the position-dependent sliding velocity across the unit cell is reproduced without any vertical motion, identifying the pressure force from the reflected rotational flow as the dominant cause of sliding; heavier particles with smaller gaps, which have stronger shear, roll at spacings where lighter particles slide.

Load-bearing premise

The conclusion rests on the assumption that the computer model reproduces the true near-wall fluid forces well enough that the backward push it predicts is a real pressure effect rather than an artifact of how the particle and pillars are built from small spheres; if that push is a modeling artifact, the mechanism for backward sliding is unproven.

Editorial extensions

If this is right

  • The ratio $l_p/R$ alone selects rolling, sliding, or trapping, so microfluidic channel design can steer rotating particles without changing the magnetic field or the particle itself.
  • Particle density shifts the transition because denser particles sediment closer to the substrate, increasing shear; a lattice that slides a light particle will roll a heavy one at the same spacing.
  • Sliding is fastest at the center of the unit cell and weakest near obstacles, making the lattice a position-dependent velocity modulator that could sort particles by their response.
  • In the trapping regime the vertical oscillation frequency is locked linearly to the rotation frequency, providing a direct readout of the local shear-to-pressure balance.
  • Because shear forces depend logarithmically on gap height, even sub-micrometer vertical displacements can flip the rolling/sliding balance, giving a sensitive control handle.

Reading between the lines

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

  • If pressure-driven sliding is generic, then reversing the magnetic field's rotation should reverse the sliding direction while leaving the critical spacing unchanged; this is a clean experiment that separates pressure-driven sliding from shear-driven rolling.
  • The linear $\xi/\xi_0$ versus $l_p/R$ relation predicts a measurable critical spacing for each particle\textendash lattice pair; mapping this critical spacing across particle radii, densities, obstacle heights, and field frequencies would test whether the transition truly collapses onto a single geometric ratio.
  • Because the fixed-height simulations reproduce the sliding profile without vertical motion, an experiment that mechanically constrains the particle height should preserve backward sliding at small spacings while eliminating the oscillatory trapping regime, cleanly separating the two phenomena.
  • The frequency-locked vertical oscillation in the trapping regime could be used as a local rheological probe, with the amplitude and phase of the vertical motion reporting the local shear-to-pressure ratio without resolving sub-micron fluid flow.
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

4 major / 4 minor

Summary. The paper studies a single rotating magnetic microparticle near a substrate patterned with a square lattice of cylindrical pillars. The authors report that the direction of net translation depends on the ratio l_p/R of the edge-to-edge pillar spacing to particle radius: at large l_p/R the particle rolls in the direction of rotation, at small l_p/R it slides backward, and in an intermediate range it becomes hydrodynamically trapped with an oscillatory vertical motion. Rotne-Prager-Yamakawa bead-model simulations at fixed particle height are used to argue that pressure forces from the obstacle lattice, rather than reduced shear, dominate the backward sliding and its position dependence across the unit cell. A control experiment with denser particles is presented as evidence that increased shear restores rolling.

Significance. If the mechanism is correct, the paper offers a clean geometric design rule for reversing the direction of a microroller's transport in a porous medium, which would be of clear interest to microfluidics and active-matter transport. The experiments are compelling at the qualitative level: the trajectories and movies show a genuine direction reversal, the denser-particle control supports the shear-pressure competition, and the linear scaling of the trapping oscillation frequency with the actuation frequency is a good check that the vertical motion is actuation-driven. The fixed-height simulation design is also a sensible way to separate height-induced shear reduction from pressure effects. However, the quantitative support for the phase boundary and the hydrodynamic fidelity of the bead-resolved RPY model are not yet sufficient to establish the central pressure-dominance claim.

major comments (4)
  1. [Experimental setup and Fig. 2] The stated parameter ranges are inconsistent with the reported crossover. With l_p (edge-to-edge) between 6 and 11 µm and R between 3 and 8 µm, the maximum l_p/R is 11/3 ≈ 3.7, yet the text places the rolling-sliding transition at l_p/R ≈ 4.0–4.4. The definition of l_p, the particle radii used in Fig. 2, or the transition values must be reconciled; as written, the central phase boundary lies outside the accessible experimental range.
  2. [Fig. 2 and Fig. S2] The normalized coupling data are presented without error bars, track counts, or per-particle statistics, even though the free-roller ξ0 values are quoted with uncertainties in the preceding paragraph. The linear decrease of ξ/ξ0 with l_p/R and the zero crossing are load-bearing quantitative claims; the authors should report the number of trajectories, standard errors, and the uncertainty of the linear fit, and show that the collapse across particle radii and rotation frequencies is statistically justified.
  3. [Numerical simulations, §4 and Fig. 4] The RPY bead model is the only support for the pressure-dominance mechanism, but it omits near-field lubrication corrections at the small gaps where the shear-pressure balance is set. The paper itself notes (ref. 31) that shear forces grow logarithmically as the separation decreases, and the experimental gaps are in the lubricated regime; a 12-bead dodecahedron with RPY mobility cannot faithfully resolve that regime. In addition, the simulated normalizing baseline ξ0 = 0.07 is imported from ref. 24 and is not checked against the experimentally measured ξ0 = 0.039–0.048 for this system. A smooth-particle boundary-integral calculation with lubrication corrections, or at least a demonstration that the sign and magnitude of ξ/ξ0 are converged with bead resolution and gap size, is needed before the pressure-force attribution is established.
  4. [Trapping regime and Fig. 3] The interpretation assumes the particle rotates synchronously with the applied field in all regimes, but no measurement or simulation verifies that the particle's angular velocity equals 2πf inside the lattice, particularly in the trapping and sliding regimes. If the particle loses phase lock, the apparent reversal could arise from asynchronous rotation rather than from pressure forces; the authors should provide direct evidence of synchronous rotation, or discuss the consequences of partial phase slip.
minor comments (4)
  1. [Throughout] There are several typographical and grammatical errors, including 'as illustrated schematically Fig. 1c' (missing 'in'), 'discontinous' in the Fig. 2 caption, and 'See fig. S8' for 'See Fig. S8'.
  2. [Fig. 4 references] The main text says the largest negative ξ/ξ0 at Δy/l_p = 0 appears in Fig. 4a, but the caption describes Fig. 4a as height dependence and Fig. 4b as position dependence; the in-text references to these panels appear to be swapped.
  3. [Fig. 2 caption and main text] The Fig. 2 caption says the black square corresponds to a denser particle with nominal radius R = 6 µm, while the main text gives R = 5.4 µm for the same particles; this nominal-versus-measured radius discrepancy should be clarified.
  4. [Intensity-to-height calibration] The conversion of the brightest-pixel intensity variations into vertical displacements of order 1 µm is not described in the main text; a brief description of the calibration or a reference to the relevant SM section would let the reader assess the magnitude of the out-of-plane motion.

Circularity Check

0 steps flagged · score 2.0 of 10

No load-bearing circularity; only a minor self-citation for the free-roller normalization baseline.

full rationale

The derivation chain is self-contained on the key point: the rolling-to-sliding reversal is observed experimentally and reproduced by fixed-height RPY simulations computed in this paper, with no parameter fitted to the reversal curve. The experimental ξ0 values (0.039–0.048) are measured baselines used only to normalize ξ, and the sign of ξ/ξ0 is set by the measured/simulated velocity direction, not by the baseline magnitude. The fixed-height simulations use standard RPY mobility (refs 35–38) and are not imported from the authors' prior work. The only self-citation is ref 24, which supplies the free-roller normalization ξ0=0.07 and the earlier pressure-force concept; this is not load-bearing because the reversal sign and the non-monotonic unit-cell velocity profile are independent of that normalization, and the heavier-particle experiment independently tests the shear-pressure balance. No uniqueness theorem or ansatz is smuggled in via citation. The RPY bead-model coarse-graining is a modeling-accuracy concern, not circularity.

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

The central mechanism relies on standard Stokes-flow hydrodynamics plus a set of modeling choices. The experiments contribute measured ξ0 baselines and trajectory data; the simulations contribute the pressure-force attribution. The main uncertainty is whether the coarse-grained RPY model and the assumed synchronous rotation capture the near-wall lubrication physics at the small sedimented gaps, and whether the chosen simulation geometry (R_o/R=2, h/R=6.4, z/R=1.75) is representative of the experimental parameter range. No new physical entities are introduced.

free parameters (4)
  • Simulation fixed particle height z/R = 1.75
    Chosen operating height at which the free roller has ξ0=0.07; the sign and magnitude of ξ/ξ0 in the lattice depend on this height.
  • Simulation obstacle-to-roller radius ratio R_o/R = 2
    Chosen geometry; experimental particles span R_o/R about 0.63 to 1.67, so the simulated transition is not a quantitative match to Fig. 2.
  • Simulation obstacle height h/R = 6.4
    Chosen geometry; experiments use obstacle height 11.5 um with particle radii 3-8 um, a ratio that is not fixed.
  • Slope and intercept of linear ξ/ξ0 versus lp/R fit = not reported
    Used to locate the rolling-sliding crossover; no fit parameters or confidence intervals are given in the main text.
assumptions (6)
  • standard math Low-Reynolds-number, quasi-steady Stokes flow
    Viscous forces dominate; used in both experiment interpretation and RPY simulations.
  • standard math No-slip boundary conditions on substrate, particle, and obstacles
    Fundamental to shear-induced rolling and pressure generation.
  • domain assumption RPY mobility with a 12-bead dodecahedron adequately captures hydrodynamic coupling including lubrication
    The central mechanistic conclusion depends on this; RPY is a far-field approximation and lubrication is only approximate.
  • domain assumption Particle rotates synchronously with the applied magnetic field at all frequencies and confinements
    Not explicitly verified in the lattice and trapping regimes; phase slip would change the effective rotation rate and force balance.
  • domain assumption Brightest-pixel intensity is a monotonic proxy for particle vertical position
    Used to infer z-oscillations and height at the unit-cell center; no independent calibration is given.
  • domain assumption Heavier silica particles sediment closer to the substrate
    Used to explain why dense particles roll where light particles slide.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Rolling, sliding and trapping of driven particles in square obstacle lattices." pith.science (2026). https://pith.science/paper/AFOPKF7C

@misc{pith2026250521702,
  author       = {Pith},
  title        = {Pith review of: Rolling, sliding and trapping of driven particles in square obstacle lattices},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/AFOPKF7C}},
  note         = {Machine review of arXiv:2505.21702}
}
read the original abstract

Transport phenomena in complex and dynamic microscopic environments are fundamentally shaped by hydrodynamic interactions. In particular, microparticle transport in porous media is governed by the delicate interplay between particle-substrate friction and pressure forces. Here, we systematically investigate the motion of externally driven rotating magnetic microparticles near a substrate patterned with a square lattice of cylindrical obstacles, a model porous medium. Remarkably, we observe a reversal in the direction of particle translation as obstacle spacing decreases, highlighting a sensitive competition between shear-induced forward rolling and pressure-driven backward sliding due to flow-field symmetry breaking. These results demonstrate the crucial role of structured environments in determining microscale active particle transport, offering novel strategies for microfluidic design, targeted cargo delivery, and tunable active materials.

Figures

Figures reproduced from arXiv: 2505.21702 by the authors.

Figure 1
Figure 1. FIG. 1. (a) Schematic of a magnetic particle rotating near [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Normalized roto-translational coupling parameter [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 4
Figure 4. FIG. 4. (a) Normalized roto-translational coupling parameter [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    C. P. Moerland, L. J. van IJzendoorn, and M. W. J. Prins, Lab Chip19, 919 (2019)

    work page 2019

  2. [2]

    W. Gao, D. Kagan, O. S. Pak, C. Clawson, S. Campuzano, E. Chuluun-Erdene, E. Shipton, E. E. Fullerton, L. E. Zhang, Liangfang, and J. Wang, Small 8, 460 (2012)

    work page 2012

  3. [3]

    J. Li, B. Esteban-Fern´ andez de´Avila, W. Gao, L. Zhang, and J. Wang, Science Robotics2, eaam6431 (2017)

    work page 2017

  4. [4]

    Palagi and P

    S. Palagi and P. Fischer, Nature Reviews Materials3, 113 (2018)

    work page 2018

  5. [5]

    Huang, S

    H. Huang, S. Yang, Y. Ying, X. Chen, J. Puigmart ´ ı-Luis, L. Zhang, and S. Pan´ e, Advanced Materials , 2305925 (2023)

    work page 2023

  6. [6]

    Ahmad, M

    B. Ahmad, M. Gauthier, G. J. Laurent, and A. Bolopion, IEEE Transactions on Robotics38, 646 (2021)

    work page 2021

  7. [7]

    Bozuyuk, P

    U. Bozuyuk, P. Wrede, E. Yildiz, and M. Sitti, Advanced Materials36, 2311462 (2024)

    work page 2024

  8. [8]

    E. M. Purcell, American journal of physics45, 3 (1977)

    work page 1977

Show all 39 references
  1. [9]

    T. A. Witten and H. Diamant, Reports on Progress in Physics83, 116601 (2020)

    2020

  2. [10]

    C. E. Sing, L. Schmid, M. F. Schneider, T. Franke, and A. Alexander-Katz, Proceedings of the National Academy of Sciences107, 535 (2010)

    2010

  3. [11]

    J. P. Steimel, J. L. Aragones, and A. Alexander-Katz, Physical Review Letters113, 178101 (2014)

    2014

  4. [12]

    L. W. Rogowski, J. Ali, X. Zhang, J. N. Wilking, H. C. Fu, and M. J. Kim, Nature Communications12, 1116 (2021)

    2021

  5. [13]

    Morimoto, T

    H. Morimoto, T. Ukai, Y. Nagaoka, N. Grobert, and T. Maekawa, Physical Review E—Statistical, Nonlinear, and Soft Matter Physics78, 021403 (2008)

    2008

  6. [14]

    Tierno, R

    P. Tierno, R. Golestanian, I. Pagonabarraga, and F. Sagu´ es, The Journal of Physical Chemistry B112, 16525 (2008)

    2008

  7. [15]

    Tierno, O

    P. Tierno, O. G¨ uell, F. Sagues, R. Golestanian, and I. Pagonabarraga, Physical Review E81, 011402 (2010)

    2010

  8. [16]

    Bozuyuk, E

    U. Bozuyuk, E. Yildiz, M. Han, S. O. Demir, and M. Sitti, Small19, 2303396 (2023)

    2023

  9. [17]

    Chamolly, E

    A. Chamolly, E. Lauga, and S. Tottori, Soft Matter16, 2611 (2020)

    2020

  10. [18]

    Z. Wu, Y. Zhang, N. Ai, H. Chen, W. Ge, and Q. Xu, Advanced Intelligent Systems4, 2100266 (2022)

    2022

  11. [19]

    C. J. Petell, K. Randene, M. Pappas, D. Sandoval, B. D. Strahl, J. S. Harrison, and J. P. Steimel, Elife10, e67525 (2021)

    2021

  12. [20]

    Bozuyuk, H

    U. Bozuyuk, H. Ozturk, and M. Sitti, Advanced Intelligent Systems , 2300099 (2023)

    2023

  13. [21]

    Alapan, U

    Y. Alapan, U. Bozuyuk, P. Erkoc, A. C. Karacakol, and M. Sitti, Science Robotics5, eaba5726 (2020)

    2020

  14. [22]

    X. Qi, S. Wang, S. Ma, K. Han, X. Bian, and X. Li, Physics of Fluids33(2021)

    2021

  15. [23]

    Y. Dou, P. M. Tzelios, D. Livitz, and K. J. Bishop, Soft Matter17, 1538 (2021)

    2021

  16. [24]

    Magrinya, P

    P. Magrinya, P. Palacios-Alonso, P. Llombart, R. Delgado-Buscalioni, A. Alexander-Katz, L. R. Arriaga, and J. L. Aragones, Proceedings of the National Academy of Sciences USA122, e2424236122 (2025)

    2025

  17. [25]

    E. B. Van Der Wee, B. C. Blackwell, F. Balboa Usabiaga, A. Sokolov, I. T. Katz, B. Delmotte, and M. M. Driscoll, Science Advances9, eade0320 (2023)

    2023

  18. [26]

    van Baalen, S

    C. van Baalen, S. Ketzetzi, A. Tintor, I. Gabay, and L. Isa, Soft Matter (2025)

    2025

  19. [27]

    Jiang, Z

    J. Jiang, Z. Yang, A. Ferreira, and L. Zhang, Advanced Intelligent Systems4, 2100279 (2022)

    2022

  20. [28]

    Xia and G

    Y. Xia and G. M. Whitesides, Angewandte Chemie International Edition37, 550 (1998)

    1998

  21. [29]

    E. E. Diel, J. W. Lichtman, and D. S. Richardson, Nature Protocols15, 2773 (2020)

    2020

  22. [30]

    H. Long, C. Lai, and C.-K. Chung, Surface and Coatings 6 Technology320, 315 (2017)

    2017

  23. [31]

    A. J. Goldmans, R. G. Cox, and H. Brenner, Chemical Engineering Science22, 637 (1967)

    1967

  24. [32]

    I. O. G¨ otze and G. Gompper, Europhysics Letters92, 64003 (2011)

    2011

  25. [33]

    J. C. Crocker and D. G. Grier, Journal of Colloid and Interface Science179, 298 (1996)

    1996

  26. [34]

    A. F. Demir¨ ors, A. Stauffer, C. Lauener, J. Cossu, S. N. Ramakrishna, J. De Graaf, C. C. Alcantara, S. Pan´ e, N. Spencer, and A. R. Studart, Soft Matter17, 1037 (2021)

    2021

  27. [35]

    Rotne and S

    J. Rotne and S. Prager, The Journal of Chemical Physics 50, 4831 (1969)

    1969

  28. [36]

    Yamakawa, The Journal of Chemical Physics53, 436 (1970)

    H. Yamakawa, The Journal of Chemical Physics53, 436 (1970)

    1970

  29. [37]

    A. M. Fiore and J. W. Swan, Journal of Fluid Mechanics 878, 544–597 (2019)

    2019

  30. [38]

    D. L. Ermak and J. A. McCammon, The Journal of Chemical Physics69, 1352 (1978)

    1978

  31. [39]

    F. B. Usabiaga, B. Kallemov, B. Delmotte, A. Bhalla, B. Griffith, and A. Donev, Communications in Applied Mathematics and Computational Science11, 217 (2016)

    2016

Pith tools

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