pith. sign in

arxiv: 2604.21033 · v1 · submitted 2026-04-22 · ❄️ cond-mat.soft · physics.flu-dyn

Orientation Dynamics of Gyrotactic Microswimmers in Turbulent Flows

Suraj Kumar Nayak , Vishwanath Shukla , Akshay Bhatnagar This is my paper

Pith reviewed 2026-05-09 22:32 UTC · model grok-4.3

classification ❄️ cond-mat.soft physics.flu-dyn
keywords verticalalignmentswimmersorientationtimedecaydimensionlessdisplacement
0
0 comments X

The pith

DNS of gyrotactic microswimmers in turbulence shows strong gyrotaxis promotes vertical alignment with shape-dependent distributions at high activity, exponential orientation autocorrelation decaying as 1/(2ψ), and transition from ballistic to diffusive motion.

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

Gyrotactic microswimmers are tiny organisms that swim while also trying to point in a preferred direction due to gravity or similar torques. In turbulent water full of random swirls, their orientations result from competition between swimming, flow shear, and reorientation. Researchers ran detailed computer simulations of fluid turbulence and tracked thousands of such swimmers, changing three properties: how elongated they are, how fast they swim relative to the flow, and how quickly they reorient. When reorientation is fast, swimmers align vertically. When slow, orientations become nearly random. Shape effects appear mainly in the tails of the orientation distributions at higher swimming speeds, with rods showing stronger vertical alignment at weak gyrotaxis. Orientation memory fades exponentially over time, with the decay rate scaling directly with the reorientation parameter. Swimmers travel straight at short times but spread randomly at long times. A reduced two-dimensional model captured the main patterns seen in the full three-dimensional runs. These behaviors help explain how microorganisms navigate chaotic flows.

Core claim

Strong gyrotaxis (smaller ψ) promotes vertical alignment of the swimmers, while weak gyrotaxis leads to nearly isotropic orientations. At small ψ the rod-shaped swimmers respond to shear by aligning with the stretching direction of the strain-rate tensor, while at large ψ the alignment with the vorticity vector is preferred. The orientation autocorrelation decays exponentially with a decay rate that scales as 1/(2ψ).

Load-bearing premise

The assumption that the background flow is homogeneous and isotropic turbulence and that swimmer dynamics are fully captured by the three non-dimensional parameters without additional effects such as biological variability or non-spherical flow interactions.

Figures

Figures reproduced from arXiv: 2604.21033 by Akshay Bhatnagar, Suraj Kumar Nayak, Vishwanath Shukla.

Figure 1
Figure 1. Figure 1: Orientation distribution in pxpypz space. Rod-like microswimmers for (a) ψ = 10 and (b) ψ = 0.5. Spheroids for (c) ψ = 10 and (d) ψ = 0.5. Spheres for (e) ψ = 10 and (f) ψ = 0.5. Swimming number is kept fixed at ϕ = 10. Projections of the orientation vector along xˆ, yˆ and zˆ are given by px ≡ p · xˆ, py ≡ p · yˆ and pz ≡ p · zˆ, respectively. of p and d is the width. In our study, we neglect the inertia … view at source ↗
Figure 2
Figure 2. Figure 2: Probability distribution functions of cos θ3(≡ p · zˆ) for swimming numbers: (a) ϕ = 0, (b) ϕ = 1, and (c) ϕ = 10. Behavior is shown for the stability numbers ψ = 0.5 (orange curves), 1 (green curves) and 10 (blue curves). Each (ϕ, ψ) combination is explored for three shapes: rods (γ = 1, square markers), spheroids (γ = 0.5, diamond markers) and spheres (γ = 0, circle markers). The inset shows the magnifie… view at source ↗
Figure 3
Figure 3. Figure 3: Probability distribution functions (PDFs) of the projections of [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Plots of orientation autocorrelation function [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: (a) Trajectories of three randomly chosen swimmers, where [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Probability distribution functions of: (a) the vertical displacement [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Reduced model for spherical microswimmers. (a) Probability distribution of orientation about the [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Reduced model for spherical microswimmers. Probability distribution functions of: (a) the vertical dis [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
read the original abstract

We study the dynamics of gyrotactic microswimmers suspended in homogeneous and isotropic turbulence by using direct numerical simulations (DNS). The swimmers are characterized by three non-dimensional parameters: their aspect ratio ($\gamma$), a dimensionless swimming speed ($\phi$), and a dimensionless reorientation time ($\psi$). Strong gyrotaxis (smaller $\psi$) promotes vertical alignment of the swimmers, while weak gyrotaxis leads to nearly isotropic orientations. At low swimming numbers, the orientation distribution is largely shape-independent with spheres and spheroids showing marginally greater vertical alignment than rods, whereas at higher activity the peaks of the distributions exhibit largely shape-independent behavior and the tails show a clear dependence on particle shape. However, at large $\psi$ rods exhibit a stronger alignment along the vertical. We observe that at small $\psi$ the rod-shaped swimmers respond to shear by aligning with the stretching direction of the strain-rate tensor, while at large $\psi$ the alignment with the vorticity vector is preferred. The orientation autocorrelation is found to decay exponentially, with a decay rate that scales as $1/(2\psi)$. Analysis of the mean-squared displacement (MSD) reveals a transition from a ballistic motion at short times to a diffusive regime at longer times. To assess the efficiency of vertical migration, we compute the probability distributions of vertical displacement over a fixed time interval and the time taken to migrate a specific vertical distance. Furthermore, we use a simplified two-dimensional model for spherical swimmers that qualitatively reproduces the key trends observed in the full three-dimensional (3D) simulations.

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 reports direct numerical simulations of gyrotactic microswimmers in homogeneous isotropic turbulence, parameterized by aspect ratio γ, dimensionless swimming speed φ, and reorientation time ψ. It finds that smaller ψ enhances vertical alignment, with shape-dependent orientation distributions at varying φ; rods align preferentially with the strain-rate stretching eigenvector at small ψ and with vorticity at large ψ. Orientation autocorrelation decays exponentially with rate scaling as 1/(2ψ). Mean-squared displacement shows a ballistic-to-diffusive transition, and vertical migration efficiency is assessed via displacement statistics. A simplified 2D model for spheres qualitatively reproduces key 3D trends.

Significance. If validated, the work offers concrete statistics on how gyrotaxis strength and shape modulate swimmer orientation and migration in turbulence, relevant to plankton ecology and bioconvection. The DNS-based conditional alignments and the observed 1/(2ψ) scaling provide testable predictions, while the reduced 2D model adds practical value. The three-parameter formulation is standard, but the significance hinges on whether the small-scale flow diagnostics are numerically reliable.

major comments (2)
  1. [DNS methodology and rod alignment results] DNS methodology and results on rod alignments: The central claims that rod-shaped swimmers align with the stretching eigenvector of the strain-rate tensor at small ψ but prefer the vorticity vector at large ψ are diagnosed from local flow gradients. No grid resolution (relative to Kolmogorov scale), Reynolds number, or convergence tests for these eigenvector-based conditional statistics are reported. Because strain and vorticity are concentrated at dissipative scales, insufficient resolution risks distorting the reported alignment preferences.
  2. [Orientation autocorrelation analysis] Orientation autocorrelation and scaling: The exponential decay with rate 1/(2ψ) is presented as a robust outcome across the parameter space. However, without reported statistical uncertainties, ensemble sizes, or sensitivity to the underlying orientation time series length, it is unclear whether this scaling holds uniformly or is influenced by finite-time sampling in the DNS.
minor comments (2)
  1. [Abstract and parameter definitions] The abstract refers to 'low swimming numbers' while the parameter is defined as dimensionless swimming speed φ; consistent terminology and explicit ranges for γ, φ, and ψ should be stated early.
  2. [Results figures] Figures showing orientation distributions and MSD curves lack error bars or indication of statistical convergence; adding these would strengthen the presentation of the shape-dependent tails and ballistic-diffusive transition.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper's central results on swimmer orientations, alignment with strain-rate eigenvectors versus vorticity, exponential autocorrelation decay scaling as 1/(2ψ), and MSD transitions are obtained from direct numerical simulations of the three-parameter model (γ, φ, ψ) in homogeneous isotropic turbulence. These quantities are computed outputs from particle trajectories and local flow gradients rather than algebraic reductions, parameter fits, or self-citations that collapse back to the inputs by construction. The simplified 2D model is introduced only to qualitatively reproduce trends and does not serve as the source of the 3D claims. No load-bearing self-citations, uniqueness theorems, or ansatzes smuggled via prior work appear in the derivation chain.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The work rests on standard fluid-dynamics assumptions and a conventional gyrotaxis model; three non-dimensional parameters are varied but not fitted to new data.

free parameters (3)
  • γ (aspect ratio)
    Non-dimensional shape parameter varied across spheres, spheroids, and rods.
  • φ (dimensionless swimming speed)
    Non-dimensional activity parameter controlling swimming relative to flow.
  • ψ (dimensionless reorientation time)
    Non-dimensional gyrotaxis strength parameter controlling reorientation rate.
axioms (2)
  • domain assumption Homogeneous and isotropic turbulence
    Background flow field assumed in the DNS setup.
  • domain assumption Standard gyrotactic torque model
    Swimmer reorientation follows conventional torque balance with gravity and flow.

pith-pipeline@v0.9.0 · 5586 in / 1406 out tokens · 27908 ms · 2026-05-09T22:32:08.513347+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [1]

    Lampert, Functional ecology , 21 (1989)

    W. Lampert, Functional ecology , 21 (1989)

    work page 1989

  2. [2]

    W. M. Durham, E. Climent, M. Barry, F. De Lillo, G. Boffetta, M. Cencini, and R. Stocker, Nature communications4, 2148 (2013)

    work page 2013

  3. [3]

    De Lillo, M

    F. De Lillo, M. Cencini, W. M. Durham, M. Barry, R. Stocker, E. Climent, and G. Boffetta, Physical review letters112, 044502 (2014)

    work page 2014

  4. [4]

    Rusconi and R

    R. Rusconi and R. Stocker, Current opinion in microbiology25, 1 (2015)

    work page 2015

  5. [5]

    Gustavsson, F

    K. Gustavsson, F. Berglund, P. Jonsson, and B. Mehlig, Physical review letters116, 108104 (2016)

    work page 2016

  6. [6]

    Sengupta, F

    A. Sengupta, F. Carrara, and R. Stocker, Nature543, 555 (2017)

    work page 2017

  7. [7]

    Gupta, A

    A. Gupta, A. Roy, A. Saha, and S. S. Ray, Physical Review Fluids5, 052601 (2020)

    work page 2020

  8. [8]

    J. K. Alageshan, A. K. Verma, J. Bec, and R. Pandit, Physical Review E101, 043110 (2020)

    work page 2020

  9. [9]

    Si and L

    X. Si and L. Fang, Physical Review Fluids7, 094501 (2022)

    work page 2022

  10. [10]

    J. Qiu, C. Marchioli, and L. Zhao, Acta Mechanica Sinica38, 722323 (2022)

    work page 2022

  11. [11]

    W. M. Durham, E. Climent, and R. Stocker, Physical Review Letters106, 238102 (2011)

    work page 2011

  12. [12]

    Borgnino, G

    M. Borgnino, G. Boffetta, F. De Lillo, and M. Cencini, Journal of Fluid Mechanics856, R1 (2018)

    work page 2018

  13. [13]

    Cencini, G

    M. Cencini, G. Boffetta, M. Borgnino, and F. De Lillo, The European Physical Journal E42, 1 (2019)

    work page 2019

  14. [14]

    J. D. Wheeler, E. Secchi, R. Rusconi, and R. Stocker, Annual review of cell and developmental biology35, 213 (2019)

    work page 2019

  15. [15]

    Z. Liu, L. Jiang, and C. Sun, Physics of Fluids34(2022)

    work page 2022

  16. [16]

    J. O. Kessler, Nature313, 218 (1985)

    work page 1985

  17. [17]

    Kamykowski, E

    D. Kamykowski, E. Milligan, and R. Reed, Journal of Plankton Research20, 1781 (1998)

    work page 1998

  18. [18]

    W. M. Durham, J. O. Kessler, and R. Stocker, science323, 1067 (2009)

    work page 2009

  19. [19]

    Santamaria, F

    F. Santamaria, F. De Lillo, M. Cencini, and G. Boffetta, Physics of Fluids26(2014)

    work page 2014

  20. [20]

    J. O. Kessler, inNonequilibrium cooperative phenomena in physics and related fields(Springer, 1984) pp. 241–248. 11

    work page 1984

  21. [21]

    G. B. Jeffery, Proceedings of the Royal Society of London. Series A, Containing papers of a mathematical and physical character102, 161 (1922)

    work page 1922

  22. [22]

    Roberts, Journal of Experimental Biology53, 687 (1970)

    A. Roberts, Journal of Experimental Biology53, 687 (1970)

    work page 1970

  23. [23]

    J. G. Mitchell, A. Okubo, and J. A. Fuhrman, Limnology and oceanography35, 123 (1990)

    work page 1990

  24. [24]

    J. S. Guasto, R. Rusconi, and R. Stocker, Annual Review of Fluid Mechanics44, 373 (2012)

    work page 2012

  25. [25]

    C. Zhan, G. Sardina, E. Lushi, and L. Brandt, Journal of fluid mechanics739, 22 (2014)

    work page 2014

  26. [26]

    J. R. Taylor and R. Stocker, Science338, 675 (2012)

    work page 2012

  27. [27]

    O. A. Croze, G. Sardina, M. Ahmed, M. A. Bees, and L. Brandt, Journal of The Royal Society Interface10, 20121041 (2013)

    work page 2013

  28. [28]

    T. J. Pedley and J. Kessler, Proceedings of the Royal Society of London. Series B. Biological Sciences231, 47 (1987)

    work page 1987

  29. [29]

    Pumir and M

    A. Pumir and M. Wilkinson, New journal of physics13, 093030 (2011)

    work page 2011

  30. [30]

    Pujara, M

    N. Pujara, M. Koehl, and E. Variano, Journal of Fluid Mechanics838, 356 (2018)

    work page 2018

  31. [31]

    Borgnino, K

    M. Borgnino, K. Gustavsson, F. De Lillo, G. Boffetta, M. Cencini, and B. Mehlig, Physical review letters123, 138003 (2019)

    work page 2019

  32. [32]

    Pujara, J.-A

    N. Pujara, J.-A. Arguedas-Leiva, C. C. Lalescu, B. Bramas, and M. Wilczek, Journal of Fluid Mechanics922, R6 (2021)

    work page 2021

  33. [33]

    Lewis, Proceedings of the Royal Society of London

    D. Lewis, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences459, 1293 (2003)

    work page 2003

  34. [34]

    P. E. Hamlington, J. Schumacher, and W. J. Dahm, Physics of Fluids20(2008)

    work page 2008

  35. [35]

    Lovecchio, E

    S. Lovecchio, E. Climent, R. Stocker, and W. M. Durham, Science advances5, eaaw7879 (2019)

    work page 2019

  36. [36]

    Fouxon and A

    I. Fouxon and A. Leshansky, Physical Review E92, 013017 (2015)

    work page 2015

  37. [37]

    Marcos, J. R. Seymour, M. Luhar, W. M. Durham, J. G. Mitchell, A. Macke, and R. Stocker, Proceedings of the National Academy of Sciences108, 3860 (2011)

    work page 2011