Pith. sign in

REVIEW 3 major objections 5 minor 79 references

Simulations of Three-dimensional Nematic Guidance of Microswimmers

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

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

Pith's one-line read A hybrid simulation coupling active Brownian particles to a lattice-Boltzmann nematic model shows that the winding profile of a 3D loop defect dictates whether confined microswimmers expand, shrink, or buckle the loop.

desk verdict A useful hybrid ABP-LBM method for living nematics with credible 2D validation, but the 3D winding-profile claims rest on single trajectories from hand-built director fields and need a stronger control. read the letter →

arxiv 2501.07816 v2 pith:3FV4F6C2 submitted 2025-01-14 cond-mat.soft physics.bio-phphysics.flu-dyn

Zeyang Mou , Yuan Li , Zhihong You , Rui Zhang This is my paper · ORCID
classification cond-mat.softphysics.bio-phphysics.flu-dyn MSC 76A15 PACS 83.80.Xz47.63.mf
keywords livingnematicactiveBrownianparticleslatticeBoltzmannmethodliquidcrystalhydrodynamics3Dtopologicaldefectsloopdefectdynamicsmicroswimmerguidancebacterialreversaltime
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

The paper develops a simulation method for living nematics, where bacteria or synthetic swimmers are dispersed in a liquid crystal, and claims that the method is faithful enough to explain existing 2D experiments and to predict new 3D behavior. The method couples active Brownian particle dynamics for the swimmers to a lattice-Boltzmann hydrodynamic model of the nematic director and flow. It reproduces stable swirling, polar jets, and their bend-induced undulation in patterned cells, with quantitative wavelengths and thresholds matching experiments. When applied to 3D loop defects, the simulations show that a loop's fate depends on its winding profile: +1/2 wedge windings and radial twist windings control expansion, shrinkage, and buckling, and the swimmer distribution feeds back on these dynamics. The paper also predicts that with stochastic swimming-direction reversals, bacteria do not necessarily accumulate in splay regions, challenging a common expectation.

What carries the argument

The carrying object is the hybrid method itself: point-like active Brownian particles (overdamped translation, orientation coupling to the director and to the flow via a Jeffery/Bretherton torque, and orientational noise) whose active stress feeds back through an anisotropic-Gaussian kernel into a lattice-Boltzmann solver of the Beris-Edwards nematodynamic equations, including the $Q$-tensor, Landau-de Gennes free energy, and Navier-Stokes flow. For the 3D loop predictions, the load-bearing construction is the local director profile of a disclination winding, Eq. (13), with tangent $\mathbf{t}$, rotation axis $\boldsymbol{\Omega}$, phase offset $\alpha$, and twist angle $\beta$; $\beta=0$ or $\pi$ gives wedge windings, $\beta=\pi/2$ pure-twist windings, $\alpha=0$ radial twist, and $\alpha=\pi/2$ tangential twist. The argument proceeds by imposing these profiles as initial conditions and letting the coupled particle-nematic dynamics evolve the loop, then attributing expansion, shrinkage, and buckling to the active flow generated by the +1/2 wedge or radial twist windings.

What would settle it

Prepare a living nematic with a single, well-characterized 3D loop defect of known winding profile and track its lateral size over time across a range of bacterial densities; for example, a homeotropic-cell wedge-twist loop that shrinks at a density where Fig. 4(b) predicts expansion would falsify the claim. A second check is to measure bacterial density near the hybrid wall as a function of reversal time, since the claim predicts that for large $t_{\mathrm{rev}}$ bacteria accumulate at bend-dominated regions, so observing splay-only accumulation regardless of reversal time would disprove it.

Watch

Extended reading notes

Core claim

The central claim is that the morphodynamics and destiny of a 3D disclination loop in a living nematic are governed by the loop's local winding profile as well as by swimmer activity, effective size, and initial distribution. Using the director-field construction $\mathbf{n} = \cos(\phi/2)\mathbf{m} + \sin(\phi/2)(\cos\beta\,\mathbf{m}\times\mathbf{t} + \sin\beta\,\mathbf{t})$ (Eq. 13) to prescribe wedge and twist windings, the authors find that in a wedge-twist loop the +1/2 wedge drives expansion or shrinkage depending on its orientation, while in a pure-twist loop radial twist windings act like the wedge and other windings generate out-of-plane flows that buckle the loop. They also find that swimmers migrate to splay regions, which can accelerate loop motion and, when the +1/2 wedge forms a dense cluster, can eject child loops; with large swimmer spacing the reduced active stress fails to stabilize the loop. When swimmers reverse direction stochastically with characteristic time $t_{\mathrm{rev}}$, the accumulation pattern changes, and for long $t_{\mathrm{rev}}$ swimmers may accumulate where bend deformation dominates, so splay accumulation is not universal.

Load-bearing premise

The 3D loop-defect predictions rest on hand-constructed initial director fields (Eq. 13); if real living-nematic loops are not well represented by those prescribed wedge or twist profiles, the finding that winding profile controls loop dynamics could be an artifact of the initialization rather than a property of the system.

Editorial extensions

If this is right

  • In thin patterned cells, the hybrid method reproduces the experimentally observed transition from dilute director-following to dense stable circulation, and then to undulated swirling, with quantitatively captured dominant undulation wavelength.
  • On periodic C-patterns, the method captures polar jet formation in splay regions and the density-dependent bend undulation, including the $\lambda \propto 1/\sqrt{\rho_j^s - \rho_j^{sc}}$ wavelength scaling and wave rupture at high density.
  • For 3D wedge-twist loops, the +1/2 wedge direction determines whether the loop expands (outward wedge, homeotropic cell) or shrinks (inward wedge, planar cell); swimmer clustering near the wedge can create child loops, while large swimmer spacing weakens active stress and accelerates shrinkage.
  • For pure-twist loops, radial twist windings control expansion and shrinkage, non-radial windings generate out-of-plane flows that buckle the loop, and an initially nonuniform swimmer distribution can transform a pure-twist loop into a wedge-twist loop through torsional buckling.
  • Varying the reversal time of swimming direction changes where bacteria accumulate, so the common picture that bacteria concentrate in splay regions is not universal; for large $t_{\mathrm{rev}}$ they can pile up at bend-dominated walls.

Reading between the lines

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

  • Because the method retains particle-level resolution, it could serve as a testbed for designing surface-anchor patterns that steer bacterial jets or stabilize defect loops in 3D, for example by combining patterned anchoring with time-varying reversal times.
  • The finding that nonuniform swimmer distributions can convert a pure-twist loop into a wedge-twist loop suggests that local activity gradients, not just global density, can rewire defect topology, a design handle the paper does not fully develop.
  • The predicted reversal-time-dependent accumulation could be tested with bacteria whose reversal frequency is controlled chemically or genetically, giving a direct experimental route to validate the model beyond reproduction of existing patterns.
  • The point-force treatment of swimmers likely underestimates excluded-volume and anchoring effects at high local densities; near dense clusters at +1/2 wedges, a squirmer-based model might change quantitative expansion rates while preserving the qualitative winding dependence.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

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 paper develops a hybrid simulation method in which active Brownian particles (microswimmers) are coupled to a lattice-Boltzmann nematic hydrodynamics solver. The microswimmers are modeled as extensile point-force dipoles that align with the local director, while their active stress enters the nematic Navier-Stokes equation through a Gaussian kernel. The method is first validated against published quasi-2D experiments on spiral-patterned and C-patterned cells, reproducing dilute following, stable swirling, undulated jets, and a wavelength scaling law. The authors then apply the method to 3D disclination loops initialized with prescribed wedge-twist and pure-twist director profiles described by Eq. (13), and report that loop expansion, shrinkage, buckling, and tumbling depend on the loop's winding profile (β and α) as well as swimmer activity, size, and initial distribution. A final section adds stochastic swimming-direction reversals and predicts density pumping in a hybrid cell, including accumulation near homeotropic walls for long reversal times.

Significance. If the 3D predictions hold, the paper would provide a validated simulation platform for living nematics with particle-level resolution, extending prior continuum and particle-based approaches. The explicit benchmarks against independent experiments (spiral and C-pattern) are a strength, as is the demonstration of stabilized undulations and the wavelength scaling. The 3D loop-defect study goes beyond uniform-activity active nematics by letting the swimmer distribution be modulated by the nematic structure, which is a promising direction. However, the quantitative 2D agreement is only semi-quantitative (factor of about two in critical density; free critical concentration in the fit), and the 3D conclusions rest on hand-imposed initial director fields and single trajectories, so the predictive claims require additional controls.

major comments (3)
  1. [III.B, Eq. (13), Fig. 4(b)] The 3D loop-defect dynamics are initialized with hand-imposed director profiles that are not equilibrium states. The manuscript itself states that this configuration is inherently unstable and that at ρ0=0 the loop rapidly shrinks due to elastic forces (Fig. 4(b)). Because swimmers align to the local director (Eq. 5) and the active stress is proportional to their local distribution, an artificially planted splay/bend pattern can seed swimmer accumulation before the loop has relaxed. The central claim that winding profile dictates expansion, shrinkage, and buckling therefore requires a passive-relaxation control: either relax the bare loop for a time comparable to the swimmer redistribution time before introducing swimmers, or measure the β and α profiles at the time active stress becomes significant and show that they match the prescribed profiles. Without such a control, the observed dynamics may be a relaxation artifact of the initialization rather than an intrinsic property of the living nematic.
  2. [III.B, Figs. 4-7] Each parameter set is represented by a single deterministic trajectory, with no independent realizations or ensemble statistics. This is particularly concerning near threshold densities, e.g., the transition between shrinkage and expansion around ρ0=2.3×10^-3 in Fig. 4(b), where stochastic fluctuations (the Dr term in Eq. 5 and the random reversal process in III.C) could move a single realization across the boundary. The 'destiny' classifications (expansion, shrinkage, buckling, tumbling) and the density-pumping prediction should be supported by at least a few independent realizations, or the authors should explicitly state that the results are single-trajectory observations whose run-to-run variability has not been assessed.
  3. [III.A.2, Fig. 2(h)] The quantitative validation is weaker than the text suggests. The simulated critical concentration for the C-pattern undulation is reported as 4.45×10^10 m^-2, compared with 2.36×10^10 m^-2 from experiment, a factor of about 1.9 discrepancy. In addition, the wavelength scaling λ ∝ 1/sqrt(ρ_s^j - ρ_sc^j) uses ρ_sc^j = 0.34 as a free fitting parameter (Fig. 2(h)). Since the abstract and conclusion describe the method as a 'faithful tool' validated against experiments, the authors should either qualify the 2D validation as semi-quantitative or discuss the origins of the factor-of-two threshold mismatch and the role of the free critical concentration in the scaling fit.
minor comments (5)
  1. [Introduction and III.B headings] There are typos in the text, for example 'protypical' in the first sentence of the Introduction and 'T o pological Defects' in the Section III.B header; please proofread the manuscript.
  2. [III.B, Fig. 6(e) caption] The caption of Fig. 6(e) refers to the 'left side of the simulation box (y < Ny/2)' containing 'triple as many microswimmers' as the right side, while the text says the '+y half' contains 3Np/4 and the '−y half' contains Np/4; please reconcile the notation so the nonuniform initial distribution is unambiguous.
  3. [III.A.2, Ref. [62]] Reference [62] is cited as 'in preparation' to support the traveling-wave interpretation of the stabilized undulation; such an unpublished reference cannot be verified and should be updated or removed before publication.
  4. [II.B, Eq. (12)] The parameter l is used in figure captions and in Eq. (12) but is never explicitly defined; please state that g⊥ = l and give its value in the parameter list.
  5. [II.B, passive stress tensor] There is a garbled sentence in Section III.B: 'with trev ∼ 2 − 12 sin a nematic DSCG solution' appears to be missing a word; in addition, the expression for Πp contains tensor-index notation that is difficult to parse and should be checked for missing parentheses and consistent index ordering.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the central 3D predictions are emergent simulation outcomes, and the 2D validation is benchmarked against independent experiments.

full rationale

The paper's chain is a hybrid ABP plus lattice-Boltzmann nematic simulation, not an analytic derivation, so circularity must be shown by an equation-level reduction or a fitted quantity presented as a prediction. No such reduction appears. The 2D validation (Figs. 1 and 2) compares against independent published experiments (Peng et al.; Koizumi et al.; Turiv et al.) and there is no statement that model parameters were fitted to those validation targets; reproducing stable swirling, undulated jets, thresholds, and scaling is therefore external support. The 3D loop-defect results initialize the director field via Eq. (13), but the winding-profile dependence of expansion, shrinkage, and buckling is an emergent output of integrating the Beris-Edwards equation and the active-stress coupling, not an identity with the prescribed profile. The paper even provides a passive control (rho0=0 in Figs. 4(b), 5(b), 6(b)) showing the handed-in loop shrinks under elasticity, allowing the active cases to be read against it. The reversal-time 'prediction' is a parameter sweep with prescribed reversal rates; the accumulation location near homeotropic walls is a simulation outcome, not a fitted target. The only self-citations are Ref. [53] for the LBM method and Ref. [62] as an in-preparation companion; neither is load-bearing because the method's fidelity is established by the external 2D benchmarks, and [62] only supplies additional details of a wave advance already compared with experiment [29]. The manuscript itself flags the main limitation in Sec. III.D: 'Further experimental validation is needed to investigate bacterial reversal near different types of anchoring surfaces and to test how this influences disclination lines or loops in the bulk.' That is an honest external-validity caveat, not circularity. A related robustness concern—whether the hand-constructed Eq. (13) profiles faithfully represent real living-nematic loops—is a physical/initial-condition issue, not a circular derivation.

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

The model introduces no new physical entities, but the central predictions depend on several hand-set dimensionless parameters and on modeling assumptions about swimmer stress, alignment, and defect initialization. The largest burden is the active-stress ansatz and the prescribed defect-loop profiles.

free parameters (8)
  • gamma_B (nematic alignment torque coefficient) = 0.004
    Controls how strongly each swimmer's orientation is pulled toward the local director; set by hand, not measured, and directly sets guidance behavior in Eqs. (5)-(6).
  • zeta_B (active stress coefficient, pusher strength) = 0.03 (2D), 0.1-0.3 (3D)
    Sets the strength of the active stress feedback onto the nematic; varies by setup and shows no experimental calibration.
  • v0 (self-propulsion speed) = 0 to 0.005
    Varied across simulations; the paper maps 0.005 to bacterial speeds of 0-16 um/s, but the mapping is not a fit.
  • sigma (swimmer excluded-volume size) = 1 to 3 lattice units
    Controls packing and activity density; l/sigma = 2 is fixed, and sigma is varied to change swimmer clustering.
  • rho_c^j (critical jet concentration for wavelength scaling) = 0.34
    Used in the fit lambda proportional to 1/sqrt(rho_s^j - rho_c^j) in Fig. 2(h); chosen to collapse simulation data, not derived independently.
  • D_r (rotational diffusion) = 1e-5
    Noise amplitude for swimmer orientation; chosen by hand.
  • LC model parameters (Gamma, A0, U, L, Ws, rho_f, eta, xi) = 0.1, 0.01, 0.35, 0.01, 1, 1, 0.33, 0.3
    Standard Beris-Edwards and lattice Boltzmann parameters chosen in dimensionless units; not fitted to experiments.
  • epsilon_LJ and gamma_LJ (Lennard-Jones repulsion) = 4.167e-4 and 0.1
    Excluded-volume interaction strength and mobility coefficient; chosen by hand.
assumptions (6)
  • domain assumption Beris-Edwards Q-tensor equation (Eq. 7) with one-elastic-constant free energy (Eq. 9) accurately describes the nematic host dynamics.
    The entire LC response, including defect dynamics, is computed from this model; no comparison to three-constant nematic elasticity is made.
  • ad hoc to paper Microswimmers act as point force dipoles with extensile sign and active stress -zeta_B sum f(r - r_i) Q^B_i (Eq. 12 and the active stress tensor).
    This is the only coupling from swimmers to the nematic; higher multipoles, surface anchoring on the swimmer body, and swimmer-induced flow disturbances beyond the dipole are neglected.
  • domain assumption Swimmer orientation evolves by Eq. (5), combining nematic alignment torque with Jeffery's equation and white noise; no explicit boundary anchoring of the swimmer body is included.
    Used in all 2D and 3D results; the alignment torque coefficient gamma_B is set by hand.
  • ad hoc to paper The local director field of a defect loop is represented by Eq. (13) with winding angle beta and phase offset alpha, and loops are initialized by prescribing such profiles.
    The 3D loop conclusions in Figs. 4-6 are generated from hand-imposed initial windings rather than from loops that nucleate spontaneously.
  • domain assumption Swimmer direction reversals occur as independent Poisson events with mean interval t_rev.
    Used in Section III.C for density pumping; real reversal may depend on local anchoring, flagellar state, and environment, as the authors themselves note.
  • domain assumption Lattice Boltzmann flow solver without thermal noise in the LC is sufficient; fluctuations come only from swimmer orientation noise.
    The paper acknowledges LBM does not include thermal fluctuations in the nematic; these may matter for defect dynamics.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Simulations of Three-dimensional Nematic Guidance of Microswimmers." pith.science (2026). https://pith.science/paper/3FV4F6C2

@misc{pith2026250107816,
  author       = {Pith},
  title        = {Pith review of: Simulations of Three-dimensional Nematic Guidance of Microswimmers},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3FV4F6C2}},
  note         = {Machine review of arXiv:2501.07816}
}
read the original abstract

It has been shown that an anisotropic liquid crystalline (LC) environment can be used to guide the self-propulsion dynamics of dispersed microswimmers, such as bacteria. This type of composite system is named "living nematic" (LN). In the dilute limit, bacteria are found to mainly follow the local director field. Beyond the dilute limit, however, they exhibit novel dynamical behaviors, from swirling around a spiral +1 defect pattern to forming undulating waves, and to active turbulence. Our current knowledge of how these different behaviors emerge at different population densities remains limited. Here we develop a hybrid method to simulate the dynamics of microswimmers dispersed in a nematic LC. Specifically, we model the microswimmers using active Brownian dynamics method, which is coupled to a hydrodynamic model of nematic LCs to describe the evolution of the flow field and the LC structure. Our method is validated by comparing to existing quasi-two-dimensional (2D) experiments, including undulated swirling around a spiral pattern and stabilized undulated jets on a periodic C-pattern. We further extend our method to three-dimensional (3D) systems by examining loop-defect dynamics. We find that the morphodynamics and destiny of a loop defect not only depend on the activity (self-propulsion velocity), effective size, and the initial distribution of the swimmers, but also rely on its winding profile. Specifically, +1/2 wedge and radial twist winding can dictate the dynamics of loop defects. By varying the characteristic reversal time, we predict that microswimmers not necessarily accumulate in splay regions. Taken together, our hybrid method provides a faithful tool to explain and guide the experiments of LNs in both 2D and 3D, sheds light on the interplay between microswimmer distribution and defect dynamics, and unravels the design principles of using LCs to control active matter.

Figures

Figures reproduced from arXiv: 2501.07816 by the authors.

Figure 1
Figure 1. FIG. 1. Microswimmer dynamics on a spiral [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Polar jets and their undulation over a periodic splay-bend “C”-pattern [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Schematic representation of the local director profile [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Dynamics of a wedge-twist loop defect in a homeotropic-anchoring cell [PITH_FULL_IMAGE:figures/full_fig_p008_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. Dynamics of a wedge-twist loop defect in a planar-anchoring cell [PITH_FULL_IMAGE:figures/full_fig_p009_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Dynamics of a pure-twist loop defect in a planar-anchoring cell [PITH_FULL_IMAGE:figures/full_fig_p010_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Interactions of microswimmers with a wedge-twist [PITH_FULL_IMAGE:figures/full_fig_p011_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

79 extracted references · 78 canonical work pages

  1. [62]

    B.-W. Park, J. Zhuang, O. Yasa, and M. Sitti,Multifunc- tional bacteria-driven microswimmers for targeted active drug delivery, ACS nano 11, 8910 (2017)

    work page 2017

  2. [1]

    Circulation and Undulation in Spiral-patterned Cell We first study the dynamics of microswimmers over a spiral-patterned surface in a thin nematic cell. Pre- vious literature [22, 26] show that in dilute dispersions, bacteria follow nonpolar spiral trajectories along the pre- imposedmolecularorientation; butastheirconcentration increases, they form unipol...

    work page

  3. [2]

    bacterial jet

    Polar Jets and Undulation in Periodic C-patterned Cell WenextexamineadifferentsystemwhereaperiodicC- patterned surface is used to guide a unidirectional bacte- rial flow (namely a “bacterial jet”) in an LN. As reported in Ref. [29], bacteria are expelled from the bend regions and condense into stable, unidirectional polar jets along the splay regions, unt...

    work page

  4. [3]

    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. modern physics 85, 1143 (2013)

    work page 2013

  5. [4]

    Saintillan, Rheology of active fluids, Annu

    D. Saintillan, Rheology of active fluids, Annu. review fluid mechanics 50, 563 (2018)

    work page 2018

  6. [5]

    Gompper, R

    G. Gompper, R. G. Winkler, T. Speck, A. Solon, C. Nar- dini, F. Peruani, H. Löwen, R. Golestanian, U. B. Kaupp, L.Alvarez, et al.,The 2020 motile active matter roadmap, J. Physics: Condens. Matter 32, 193001 (2020)

    work page 2020

  7. [6]

    Bechinger, R

    C. Bechinger, R. Di Leonardo, H. Löwen, C. Reichhardt, G. Volpe, and G. Volpe,Active particles in complex and crowded environments, Rev. modern physics 88, 045006 (2016)

    work page 2016

  8. [7]

    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, Nat. Commun. 15, 1406 (2024)

    work page 2024

Show all 79 references
  1. [8]

    S. S. Suarez and A. Pacey,Sperm transport in the female reproductive tract, Hum. reproduction update 12, 23 (2006)

    2006

  2. [9]

    Kantsler, J

    V. Kantsler, J. Dunkel, M. Blayney, and R. E. Goldstein, Correction: Rheotaxis facilitates upstream navigation of mammalian sperm cells, Elife 3, e03521 (2014)

    2014

  3. [10]

    I. I. Smalyukh, J. Butler, J. D. Shrout, M. R. Parsek, and G. C. L. Wong,Elasticity-mediated nematiclike bacterial organization in model extracellular DNA matrix, Phys. Rev. E 78, 030701 (2008)

    2008

  4. [11]

    Flemming and J

    H.-C. Flemming and J. Wingender,The biofilm matrix, Nat. reviews microbiology 8, 623 (2010)

    2010

  5. [12]

    D. J. Lemon, X. Yang, P. Srivastava, Y.-Y. Luk, and A. G. Garza, Polymertropism of rod-shaped bacteria: movement along aligned polysaccharide fibers, Sci. re- ports 7, 7643 (2017)

    2017

  6. [13]

    Repula, E

    A. Repula, E. Abraham, V. Cherpak, and I. I. Sma- lyukh, Biotropic liquid crystal phase transformations in cellulose-producing bacterial communities, Proc. Na- tional Acad. Sci. 119, e2200930119 (2022)

    2022

  7. [14]

    Gao and J

    W. Gao and J. Wang,Synthetic micro/nanomotors in drug delivery, Nanoscale 6, 10486 (2014)

    2014

  8. [15]

    K. J. Rao, F. Li, L. Meng, H. Zheng, F. Cai, and W. Wang,A force to be reckoned with: a review of syn- thetic microswimmers powered by ultrasound, Small 11, 2836 (2015)

    2015

  9. [16]

    Z. Wu, Y. Chen, D. Mukasa, O. S. Pak, and W. Gao, Medical micro/nanorobots in complex media, Chem. Soc. Rev. 49, 8088 (2020)

    2020

  10. [17]

    Masocha, B

    W. Masocha, B. Robertson, M. E. Rottenberg, J. Mh- langa, L. Sorokin, K. Kristensson,et al., Cerebral vessel laminins and IFN-γ define Trypanosoma brucei brucei penetration of the blood-brain barrier, The J. clinical investigation 114, 689 (2004)

    2004

  11. [18]

    Patteson, A

    A. Patteson, A. Gopinath, M. Goulian, and P. Arratia, Running and tumbling with E. coli in polymeric solutions, Sci. reports 5, 15761 (2015)

    2015

  12. [19]

    J. R. Gomez-Solano, A. Blokhuis, and C. Bechinger,Dy- namics of Self-Propelled Janus Particles in Viscoelastic Fluids, Phys. Rev. Lett. 116, 138301 (2016)

    2016

  13. [20]

    C.-k. Tung, C. Lin, B. Harvey, A. G. Fiore, F. Ardon, M. Wu, and S. S. Suarez,Fluid viscoelasticity promotes collective swimming of sperm, Sci.reports 7,3152(2017)

    2017

  14. [21]

    D. Yuan, Q. Zhao, S. Yan, S.-Y. Tang, G. Alici, J. Zhang, and W. Li,Recent progress of particle migration in vis- coelastic fluids, Lab on a Chip 18, 551 (2018)

    2018

  15. [22]

    Ishimoto and E

    K. Ishimoto and E. A. Gaffney,Hydrodynamic clustering of human sperm in viscoelastic fluids, Sci. Reports 8, 15600 (2018)

    2018

  16. [23]

    S. Liu, S. Shankar, M. C. Marchetti, and Y. Wu,Vis- coelastic control of spatiotemporal order in bacterial ac- tive matter, Nature 590, 80 (2021)

    2021

  17. [24]

    C. Peng, T. Turiv, Y. Guo, Q.-H. Wei, and O. D. Lavren- tovich, Command of active matter by topological defects and patterns, Science 354, 882 (2016)

    2016

  18. [25]

    Kumar, T

    A. Kumar, T. Galstian, S. K. Pattanayek, and S. Rainville, The motility of bacteria in an anisotropic liquid environment, Mol. Cryst. Liquid Cryst. 574, 33 (2013)

    2013

  19. [26]

    S. Zhou, A. Sokolov, O. D. Lavrentovich, and I. S. Aran- son, Living liquid crystals, Proc. National Acad. Sci. 111, 1265 (2014)

    2014

  20. [27]

    Sokolov, S

    A. Sokolov, S. Zhou, O. D. Lavrentovich, and I. S. Aran- son, Individual behavior and pairwise interactions be- tween microswimmers in anisotropic liquid, Phys. Rev. E 91, 013009 (2015)

    2015

  21. [28]

    Koizumi, T

    R. Koizumi, T. Turiv, M. M. Genkin, R. J. Lastowski, H. Yu, I. Chaganava, Q.-H. Wei, I. S. Aranson, and O. D. Lavrentovich,Control of microswimmers by spiral nematic vortices: Transition from individual to collective motion and contraction, expansion, and stable circulation of...

    2020

  22. [29]

    M. M. Genkin, A. Sokolov, O. D. Lavrentovich, and I. S. Aranson, Topological Defects in a Living Nematic En- snare Swimming Bacteria, Phys. Rev. X 7, 011029 (2017)

    2017

  23. [30]

    M. M. Genkin, A. Sokolov, and I. S. Aranson,Sponta- neous topological charging of tactoids in a living nematic, New J. Phys. 20, 43027 (2018)

    2018

  24. [31]

    Turiv, R

    T. Turiv, R. Koizumi, K. Thijssen, M. M. Genkin, H. Yu, C. Peng, Q.-H. Wei, J. M. Yeomans, I. S. Aranson, A. Doostmohammadi, and O. D. Lavrentovich, Polar jets of swimming bacteria condensed by a patterned liquid crystal, Nat. Phys. 16, 481 (2020)

    2020

  25. [32]

    P. C. Mushenheim, R. R. Trivedi, H. H. Tuson, D. B. Weibel, and N. L. Abbott, Dynamic self-assembly of motile bacteria in liquid crystals, Soft Matter 10, 88 (2013)

    2013

  26. [33]

    P.C.Mushenheim, R.R.Trivedi, S.S.Roy, M.S.Arnold, D. B. Weibel, and N. L. Abbott,Effects of confinement, surface-induced orientations and strain on dynamical be- haviors of bacteria in thin liquid crystalline films, Soft Matter 11, 6821 (2015)

    2015

  27. [34]

    R. R. Trivedi, R. Maeda, N. L. Abbott, S. E. Spagnolie, and D. B. Weibel,Bacterial transport of colloids in liquid crystalline environments, Soft Matter 11, 8404 (2015)

    2015

  28. [35]

    Lintuvuori, A

    J. Lintuvuori, A. Würger, and K. Stratford, Hydrody- namics Defines the Stable Swimming Direction of Spher- ical Squirmers in a Nematic Liquid Crystal, Phys. Rev. Lett. 119, 68001 (2017), publisher: American Physical Society

    2017

  29. [36]

    Gautam and J

    B. Gautam and J. S. Lintuvuori,Microswimmers Knead Nematics into Cholesterics, Phys. Rev. Lett. 132, 238301 (2024). 14

    2024

  30. [37]

    H. Chi, M. Potomkin, L. Zhang, L. Berlyand, and I. S. Aranson, Surface anchoring controls orientation of a mi- croswimmer in nematic liquid crystal, Commun. Phys. 3, 162 (2020)

    2020

  31. [38]

    S. Zhou, O. Tovkach, D. Golovaty, A. Sokolov, I. S. Aranson, and O. D. Lavrentovich, Dynamic states of swimming bacteria in a nematic liquid crystal cell with homeotropic alignment, New J. Phys. 19, 055006 (2017)

    2017

  32. [39]

    Daddi-Moussa-Ider and A

    A. Daddi-Moussa-Ider and A. M. Menzel,Dynamics of a simple model microswimmer in an anisotropic fluid: Im- plications for alignment behavior and active transport in a nematic liquid crystal, Phys. Review Fluids 3, 094102 (2018)

    2018

  33. [40]

    S. D. Ryan, L. Berlyand, B. M. Haines, and D. A. Karpeev, A Kinetic Model for Semidilute Bacterial Sus- pensions, Multiscale Model. Simul. 11, 1176 (2013)

    2013

  34. [41]

    Lee and M

    K.-W. Lee and M. G. Mazza,Stochastic rotation dynam- ics for nematic liquid crystals, The J. chemical physics 142, https://doi.org/10.1063/1.4919310 (2015)

    2015 doi

  35. [42]

    T. N. Shendruk and J. M. Yeomans,Multi-particle colli- sion dynamics algorithm for nematic fluids, Soft Matter 11, 5101 (2015)

    2015

  36. [43]

    Mandal and M

    S. Mandal and M. G. Mazza,Multiparticle collision dy- namics for tensorial nematodynamics, Phys. Review E 99, 063319 (2019)

    2019

  37. [44]

    Híjar,Dynamics of defects around anisotropic parti- cles in nematic liquid crystals under shear, Phys

    H. Híjar,Dynamics of defects around anisotropic parti- cles in nematic liquid crystals under shear, Phys. Review E 102, 062705 (2020)

    2020

  38. [45]

    Mandal and M

    S. Mandal and M. G. Mazza,Multiparticle collision dy- namics simulations of a squirmer in a nematic fluid, The Eur. Phys. J. E 44, 64 (2021)

    2021

  39. [46]

    Kozhukhov and T

    T. Kozhukhov and T. N. Shendruk,Mesoscopic simula- tions of active nematics, Sci. Adv. 8, eabo5788 (2022)

    2022

  40. [47]

    Macías-Durán, V

    J. Macías-Durán, V. Duarte-Alaniz, and H. Híjar,Active nematic liquid crystals simulated by particle-based meso- scopic methods, Soft Matter 19, 8052 (2023)

    2023

  41. [48]

    Kozhukhov, B

    T. Kozhukhov, B. Loewe, and T. N. Shendruk,Mitigat- ing density fluctuations in particle-based active nematic simulations, Commun. Phys. 7, 251 (2024)

    2024

  42. [49]

    Mandal, T

    S. Mandal, T. J. Mason, A. C. Croft, and M. G. Mazza, Cooperativity of Confined Nematic Microswim- mers: From One to Many, Phys. Review Lett. 134, 128302 (2025)

    2025

  43. [50]

    Denniston, E

    C. Denniston, E. Orlandini, and J. Yeomans, Lattice Boltzmann simulations of liquid crystal hydrodynamics, Phys. Review E 63, 056702 (2001)

    2001

  44. [51]

    Li, X.-q

    H. Li, X.-q. Shi, M. Huang, X. Chen, M. Xiao, C. Liu, H. Chaté, and H. Zhang,Data-driven quantitative model- ing of bacterial active nematics, Proc. National Acad. Sci. 116, 777 (2019)

    2019

  45. [52]

    H. Li, H. Chaté, M. Sano, X.-q. Shi, and H. P. Zhang,Ro- bust Edge Flows in Swarming Bacterial Colonies, Phys. Rev. X 14, 041006 (2024)

    2024

  46. [53]

    J. S. Dai, Euler–Rodrigues formula variations, quater- nion conjugation and intrinsic connections , Mech. Mach. Theory 92, 144 (2015)

    2015

  47. [54]

    Ishimoto, Jeffery’s orbits and microswimmers in flows: A theoretical review, J

    K. Ishimoto, Jeffery’s orbits and microswimmers in flows: A theoretical review, J. Phys. Soc. Japan 92, 062001 (2023)

    2023

  48. [55]

    Zhang, T

    R. Zhang, T. Roberts, I. S. Aranson, and J. J. De Pablo, Lattice Boltzmann simulation of asymmetric flow in ne- matic liquid crystals with finite anchoring , The J. Chem. Phys. 144 (2016)

    2016

  49. [56]

    A. N. Beris and B. J. Edwards,Thermodynamics of flow- ing systems: with internal microstructure, 36 (Oxford University Press, USA, 1994)

    1994

  50. [57]

    Darnton, L

    N. Darnton, L. Turner, K. Breuer, and H. C. Berg,Mov- ing fluid with bacterial carpets, Biophys. journal86, 1863 (2004)

    2004

  51. [58]

    Edwards and J

    S. Edwards and J. Yeomans,Spontaneous flow states in active nematics: a unified picture, Europhys. letters 85, 18008 (2009)

    2009

  52. [59]

    Drescher, R

    K. Drescher, R. E. Goldstein, N. Michel, M. Polin, and I. Tuval, Direct measurement of the flow field around swimming microorganisms, Phys. Review Lett. 105, 168101 (2010)

    2010

  53. [60]

    Hintsche, V

    M. Hintsche, V. Waljor, R. Großmann, M. J. Kühn, K. M. Thormann, F. Peruani, and C. Beta,A polar bun- dle of flagella can drive bacterial swimming by pushing, pulling, or coiling around the cell body, Sci. reports 7, 16771 (2017)

    2017

  54. [61]

    A. V. Singh, Z. Hosseinidoust, B.-W. Park, O. Yasa, and M. Sitti, Microemulsion-based soft bacteria-driven microswimmers for active cargo delivery, ACS nano 11, 9759 (2017)

    2017

  55. [63]

    Sourjik, and M

    Y.Alapan, O.Yasa, O.Schauer, J.Giltinan, A.F.Tabak, V. Sourjik, and M. Sitti, Soft erythrocyte-based bacte- rial microswimmers for cargo delivery, Sci. robotics 3, eaar4423 (2018)

    2018

  56. [64]

    J. Wu, M. Liu, Z. Mou, Y. Li, R. Wang, Z. Asilehan, Q. Shi, Z. You, R. Zhang, J. Jiang, and C. Peng,Pro- grammable double traveling waves in living liquid crystals (2025), in preparation

    2025

  57. [65]

    Ramaswamy,The mechanics and statistics of active matter, Annu

    S. Ramaswamy,The mechanics and statistics of active matter, Annu. Rev. Condens. Matter Phys. 1, 323 (2010)

    2010

  58. [66]

    Duclos, R

    G. Duclos, R. Adkins, D. Banerjee, M. S. Peterson, M. Varghese, I. Kolvin, A. Baskaran, R. A. Pelcovits, T. R. Powers, A. Baskaran,et al., Topological structure and dynamics of three-dimensional active nematics, Sci- ence 367, 1120 (2020)

    2020

  59. [67]

    Binysh, i

    J. Binysh, i. c. v. Kos, S. Čopar, M. Ravnik, and G.P.Alexander, Three-Dimensional Active Defect Loops, Phys. Rev. Lett. 124, 088001 (2020)

    2020

  60. [68]

    C. Long, X. Tang, R. L. Selinger, and J. V. Selinger,Ge- ometry and mechanics of disclination lines in 3D nematic liquid crystals, Soft Matter 17, 2265 (2021)

    2021

  61. [69]

    Kralj, M

    N. Kralj, M. Ravnik, and i. c. v. Kos,Defect Line Coars- ening and Refinement in Active Nematics, Phys. Rev. Lett. 130, 128101 (2023)

    2023

  62. [70]

    Goral, E

    M. Goral, E. Clement, T. Darnige, T. Lopez-Leon, and A. Lindner,Frustrated ‘run and tumble’of swimming Es- cherichia coli bacteria in nematic liquid crystals, Inter- face Focus 12, 20220039 (2022)

    2022

  63. [71]

    A. G. Prabhune, A. S. García-Gordillo, I. S. Aranson, T. R. Powers, and N. Figueroa-Morales,Bacteria Navi- gate Anisotropic Media using a Flagellar Tug-of-Oars, PRX Life 2, 033004 (2024)

    2024

  64. [72]

    H. C. Berg and D. A. Brown,Chemotaxis in Escherichia coli analysed by three-dimensional tracking, nature 239, 500 (1972)

    1972

  65. [73]

    Taute, S

    K. Taute, S. Gude, S. Tans, and T. Shimizu, High- throughput 3D tracking of bacteria on a standard phase contrast microscope, Nat. communications 6, 1 (2015). 15

    2015

  66. [74]

    Figueroa-Morales, R

    N. Figueroa-Morales, R. Soto, G. Junot, T. Darnige, C. Douarche, V. A. Martinez, A. Lindner, and E. Clé- ment, 3D Spatial Exploration by E. coli Echoes Motor Temporal Variability, Phys. Rev. X 10, 021004 (2020)

    2020

  67. [75]

    C. W. Chan, Z. Yang, Z. Gan, and R. Zhang, In- terplay of chemotactic force, Péclet number, and di- mensionality dictates the dynamics of auto-chemotactic chiral active droplets, The J. Chem. Phys. 161, https://doi.org/10.1063/5.0207355 (2024)

    2024 doi

  68. [76]

    S.Čopar, J.Aplinc, i.c.v.Kos, S.Žumer,andM.Ravnik, Topology of Three-Dimensional Active Nematic Turbu- lence Confined to Droplets, Phys. Rev. X 9, 031051 (2019)

    2019

  69. [77]

    L. N. Carenza, G. Gonnella, D. Marenduzzo, and G. Ne- gro, Rotation and propulsion in 3D active chiral droplets, Proc. National Acad. Sci. 116, 22065 (2019)

    2019

  70. [78]

    C. A. Whitfield, T. C. Adhyapak, A. Tiribocchi, G. P. Alexander, D. Marenduzzo, and S. Ramaswamy,Hydro- dynamic instabilities in active cholesteric liquid crystals, The Eur. Phys. J. E 40, 1 (2017)

    2017

  71. [79]

    Boule, S

    M.-A. Boule, S. Rainville, and T. Galstian, Dynamic guiding of bacteria in lyotropic chromonic liquid crystal using magnetic field, Mol. Cryst. Liquid Cryst. 712, 10 (2020)

    2020

Pith tools

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