pith. sign in

arxiv: 2604.03193 · v2 · submitted 2026-04-03 · ⚛️ physics.flu-dyn · cond-mat.soft

Bifurcations in Stokes Flow Sedimentation

Elias Huseby , Pierre Mathier , Meera Das , Arjun Menezes , Theo Witkamp , Ziqi Wang , Bernhard Mehlig , Greg A. Voth This is my paper

Pith reviewed 2026-05-13 18:10 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn cond-mat.soft
keywords stokes sedimentationbifurcationshelical particlescenter of mass offsetlow reynolds numberhydrodynamic interactionslimit cyclesparity time-reversal symmetry
0
0 comments X

The pith

Sedimentation of helical particles switches from complex limit cycles to simple alignment when center-of-mass offset exceeds a threshold below one percent of particle length.

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

The paper shows that helical ribbons sedimenting at low Reynolds number exhibit a sharp change in behavior controlled by the position of their center of mass relative to their geometric center. Experiments with controlled offsets identify a transition from rich, cycling dynamics to a single stable settling state as the offset grows. Simulations map this change onto an alignment bifurcation surface in three-dimensional offset space, inside which limit cycles appear through Hopf and homoclinic bifurcations. The transition is so sensitive that it occurs for offsets smaller than one percent of the ribbon length. A reader would care because it demonstrates how minute geometric asymmetries dictate whether particles tumble irregularly or fall steadily.

Core claim

An alignment bifurcation surface exists in the three-dimensional space of center-of-mass offsets for sedimenting helical particles. Inside the surface, particles follow limit cycles that arise through Hopf and homoclinic bifurcations; outside it they converge to a single attracting state. Cocentered particles possess three parity time-reversal symmetries under reflections normal to the eigenvectors of the translation-rotation coupling tensor, and any offset that preserves at least one of these symmetries permits closed orbits.

What carries the argument

The alignment bifurcation surface in the three-dimensional space of center-of-mass offsets, which separates regimes of complex and simple sedimentation dynamics.

If this is right

  • Offsets large enough to break all PT symmetries drive particles to a single attracting orientation rather than sustained cycling.
  • Limit cycles inside the bifurcation surface appear specifically via Hopf and homoclinic bifurcations.
  • Preservation of even one PT symmetry by the offset permits closed orbits within the surface.
  • The reference cocentered case sits at the center of the surface, with all symmetries intact until an offset breaks them.

Where Pith is reading between the lines

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

  • The same offset-controlled bifurcation structure may appear for other chiral or asymmetric shapes sedimenting in Stokes flow.
  • Microfluidic devices could exploit precise center-of-mass placement to steer particles into desired orientations during sedimentation.
  • Natural particles such as microorganisms or sediment grains might display analogous sensitivity to tiny mass asymmetries.

Load-bearing premise

The immersed boundary method mobility tensors accurately represent the hydrodynamic interactions of the experimental helical ribbons without significant numerical artifacts, and center-of-mass offsets can be set experimentally with precision sufficient to locate the bifurcation threshold.

What would settle it

Sedimentation experiments that vary the center-of-mass offset of helical ribbons in steps around 0.5 percent of particle length and record whether trajectories switch from closed limit cycles to steady alignment at the predicted surface location.

Figures

Figures reproduced from arXiv: 2604.03193 by Arjun Menezes, Bernhard Mehlig, Elias Huseby, Greg A. Voth, Meera Das, Pierre Mathier, Theo Witkamp, Ziqi Wang.

Figure 1
Figure 1. Figure 1: Spatial trajectories of five helical ribbons with identical shape and initial [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Orientation dynamics of a cocentered particle. (a) gravity vector trajectories in [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Helical ribbons with holes where metal spheres are epoxied to move the center [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Increasing 𝑦 (minor axis) offset of the center of mass. (a) through (e) show experimental phase diagrams as the center of mass moves along the 𝑦 axis. (f) shows a bifurcation diagram of the tilt (𝜃) position of the fixed points from numerical simulations. Saddles are green and centers are yellow. Vertical lines in (f) show the center of mass offset of each experimental data set. rotation by 𝜋. We compare e… view at source ↗
Figure 5
Figure 5. Figure 5: Increasing 𝑧 offset of the center of mass. (a) through (e) show experimental phase diagrams as the center of mass moves along the 𝑧 axis. (f) shows a bifurcation diagram of the spin (𝜓) position of the fixed points from numerical simulations. Saddles are green and centers are yellow. Vertical lines in (f) show the center of mass offset of each experimental data set. movies linked from Appendix 9.5 which de… view at source ↗
Figure 6
Figure 6. Figure 6: Increasing offset of the center of mass along the intermediate ( [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Computationally generated phase spaces for a helical ribbon a with small center [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Surfaces defining bifurcations as a function of the center of mass position. (a) [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: A schematic of dynamical regimes in a slice of the helical ribbon’s bifurcation [PITH_FULL_IMAGE:figures/full_fig_p018_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Phase spaces as the center of mass is moved relative to the center of mobility [PITH_FULL_IMAGE:figures/full_fig_p018_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Phase spaces as the center of mass is moved relative to the center of mobility [PITH_FULL_IMAGE:figures/full_fig_p019_11.png] view at source ↗
read the original abstract

Particles whose shapes couple translation to rotation display a rich array of behaviors as they sediment at low Reynolds number. We introduce a unifying perspective in which the possible dynamical regimes and bifurcations between them can be understood. We use experimental measurements of helical ribbons, with controlled center of mass offsets, to identify the key bifurcation from complex dynamics to a single attracting state as the magnitude of the offset increases. The sedimentation dynamics are very sensitive to small center of mass offsets, with the bifurcation occurring for offsets less than one percent of the particle length. Using mobility tensors obtained from immersed boundary method simulations, we simulate helical particle sedimentation and identify an alignment bifurcation surface, defined in the three dimensional space of center of mass offsets, that separates simple from complex sedimentation dynamics. Inside this surface we find limit cycles which emerge through Hopf and homoclinic bifurcations. Cocentered particles with coincident centers of force and mobility provide a reference case at the center of the bifurcation surface. We show how the geometric and dynamical symmetries of sedimenting cocentered particles are broken as the center of force offset moves away from the cocentered case. Three parity time-reversal (PT) symmetries exist for all cocentered particles under reflections normal to the eigenvectors of its translation-rotation coupling tensor. When a center of force offset preserves at least one of these PT symmetries, then there are closed orbits for particles inside the alignment bifurcation surface.

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 / 3 minor

Summary. The manuscript examines low-Reynolds-number sedimentation of particles whose shape couples translation and rotation, focusing on helical ribbons. Through controlled experiments with center-of-mass offsets and immersed-boundary-method-derived mobility tensors, it identifies an alignment bifurcation surface in three-dimensional offset space that separates simple attracting states from complex dynamics. Inside the surface, limit cycles arise via Hopf and homoclinic bifurcations; the transition occurs for offsets below 1% of particle length. The analysis centers on the breaking of three PT symmetries present in the cocentered reference case, with closed orbits preserved when at least one symmetry remains.

Significance. If the reported sensitivity and bifurcation structure hold, the work supplies a symmetry-based organizing principle for the phase space of sedimenting chiral particles, directly linking geometric offsets to qualitative changes in long-time behavior. The combination of precision experiments and mobility-tensor simulations provides a reproducible route to predict regimes, with potential relevance to microfluidic sorting, colloidal assembly, and biological sedimentation.

major comments (2)
  1. [Abstract / experimental results] Abstract and experimental results section: the central sensitivity claim (bifurcation for offsets <1% of particle length) is load-bearing yet lacks reported uncertainty quantification, calibration precision for offset placement, or explicit data-exclusion criteria; without these, the precise location of the bifurcation surface cannot be assessed.
  2. [Numerical methods] Numerical methods and validation: the assumption that IBM mobility tensors fully capture the hydrodynamic interactions for the experimental ribbons (including any higher-order or wall effects) is not cross-validated against the measured sedimentation trajectories at the reported offset thresholds; this directly affects the reliability of the simulated bifurcation surface.
minor comments (3)
  1. [Symmetry analysis] Notation for the translation-rotation coupling tensor and its eigenvectors should be introduced once with a clear definition before the PT-symmetry discussion.
  2. [Figures] Figure captions for the bifurcation surface and phase portraits should explicitly state the range of offsets shown and the integration time used to classify simple versus complex dynamics.
  3. [Introduction] A brief statement on the Reynolds-number range achieved in both experiment and simulation would strengthen the low-Re assumption.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the positive assessment of its significance. We address each major comment below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [Abstract / experimental results] Abstract and experimental results section: the central sensitivity claim (bifurcation for offsets <1% of particle length) is load-bearing yet lacks reported uncertainty quantification, calibration precision for offset placement, or explicit data-exclusion criteria; without these, the precise location of the bifurcation surface cannot be assessed.

    Authors: We agree that explicit uncertainty quantification and calibration details are needed to support the central claim. In the revised manuscript we will add a dedicated subsection in the experimental results reporting the precision of center-of-mass offset placement (via microfabrication tolerances and optical calibration), estimated uncertainties on the offset values, and the explicit criteria used to exclude trajectories (e.g., excessive noise or wall proximity). These additions will allow readers to assess the robustness of the reported bifurcation threshold. revision: yes

  2. Referee: [Numerical methods] Numerical methods and validation: the assumption that IBM mobility tensors fully capture the hydrodynamic interactions for the experimental ribbons (including any higher-order or wall effects) is not cross-validated against the measured sedimentation trajectories at the reported offset thresholds; this directly affects the reliability of the simulated bifurcation surface.

    Authors: We concur that direct cross-validation strengthens the link between simulation and experiment. The revised numerical methods section will include quantitative comparisons of simulated versus measured sedimentation trajectories for helical ribbons at offsets both inside and near the bifurcation surface. These comparisons will confirm that the IBM-derived mobility tensors reproduce the observed dynamics, including any wall effects present in the experimental setup, thereby validating the location of the simulated bifurcation surface. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper derives its central claims—the existence of an alignment bifurcation surface in center-of-mass offset space and the associated Hopf/homoclinic bifurcations—from independent experimental measurements of helical ribbons and from mobility tensors computed via immersed-boundary-method simulations. These inputs are not obtained by fitting parameters to the same sedimentation trajectories that are later “predicted”; the symmetry-breaking analysis follows directly from the translation-rotation coupling tensor of the cocentered reference case without self-referential definitions or load-bearing self-citations. No step reduces by construction to a quantity defined from the target result itself.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The analysis rests on the Stokes-flow approximation and the accuracy of mobility tensors computed by immersed boundary methods; no free parameters are introduced beyond the experimentally controlled offsets, and no new entities are postulated.

axioms (2)
  • domain assumption Low Reynolds number Stokes flow approximation holds for the sedimentation dynamics
    Invoked throughout the experimental and simulation analysis of particle trajectories.
  • domain assumption Mobility tensors from immersed boundary simulations accurately represent the hydrodynamic coupling for the helical ribbons
    Used to generate the simulated trajectories and identify the bifurcation surface.

pith-pipeline@v0.9.0 · 5569 in / 1342 out tokens · 39561 ms · 2026-05-13T18:10:19.186863+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

3 extracted references · 3 canonical work pages

  1. [1]

    & Eaton, John K.2010 Turbulent dispersed multiphase flow.Annual Review of Fluid Mechanics42(Volume 42, 2010), 111–133

    Balachandar, S. & Eaton, John K.2010 Turbulent dispersed multiphase flow.Annual Review of Fluid Mechanics42(Volume 42, 2010), 111–133. Barron, Laurence D2009Molecular light scattering and optical activity. Cambridge University Press. 25 Bender, Carl M. & Boettcher, Stefan1998 Real spectra in non-hermitian hamiltonians having PT symmetry.Phys. Rev. Lett.80...

    work page arXiv 2010

  2. [2]

    Hobden, M

    Springer Science & Business Media. Hobden, M. V.1967 Optical activity in a non-enantiomorphous crystal silver gallium sulphide.Nature 216(5116), 678–678. Huseby, Elias, Gissinger, Josephine, Candelier, Fabien, Pujara, Nimish, Verhille, Gautier, Mehlig, Bernhard & Voth, Greg2025 Helical ribbons: Simple chiral sedimentation.Physical Review Fluids10(2), 0241...

    work page 1967

  3. [3]

    Oehmke, Theresa B., Bordoloi, Ankur D., Variano, Evan & Verhille, Gautier2021 Spinning and tumbling of long fibers in isotropic turbulence.Physical Review Fluids6(4), 044610. Pasteur, Louis1848 Recherches sur les relations qui peuvent exister entre la forme cristalline, la composition chimique et le sens de la polarisation rotatoire.Comptes Rendus de l’Ac...

    work page arXiv 2026