pith. sign in

arxiv: 2511.21359 · v2 · pith:RBUTCXBZnew · submitted 2025-11-26 · ❄️ cond-mat.soft · physics.bio-ph

Spatiotemporal Control of Charge +1 Topological Defects in Polar Active Matter

Birte C. Geerds , Abhinav Singh , Mathieu Dedenon , Daniel J. G. Pearce , Frank J\"ulicher , Ivo F. Sbalzarini , Karsten Kruse This is my paper

Pith reviewed 2026-05-25 07:37 UTC · model grok-4.3

classification ❄️ cond-mat.soft physics.bio-ph
keywords topological defectsactive polar fluiddefect controlproportional integral controllerannular activitypolar boundary conditionsactive liquid crystals
0
0 comments X

The pith

A proportional integral controller steers +1 topological defects along complex paths in an active polar fluid by varying active annulus size and boundary orientation.

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

The paper shows how to direct the motion of a +1 topological defect inside a disk-shaped active polar fluid. Localizing activity to an annulus makes the defect circle the center at a speed and radius set by the annulus dimensions and the polar boundary angle. An ansatz for the polar field yields explicit formulas for those speeds and radii, which are then used to build a proportional-integral feedback loop that continuously changes the annulus size and boundary angle. The loop can therefore drive the defect along any prescribed trajectory inside the disk.

Core claim

If activity is localized in an annulus within the disk, the defect moves on a circular trajectory around the center of the disk. Using an ansatz for the polar field, the dependence of the angular speed and the circle radius on the boundary orientation of the polar field and the active annulus is determined. Using a proportional integral controller, the defect is guided along complex trajectories by changing the active annulus size and the boundary orientation.

What carries the argument

Proportional-integral feedback that adjusts the radius of the active annulus and the polar boundary orientation to enforce a target defect path.

If this is right

  • The defect can be held at a chosen location for a chosen duration.
  • Arbitrary smooth trajectories inside the disk become reachable by suitable time-dependent annulus and boundary inputs.
  • Defect motion can be controlled without direct forcing on the defect core itself.
  • The same inputs can be used to switch between different circular orbits or to reverse direction.

Where Pith is reading between the lines

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

  • The approach may generalize to other defect charges or to non-circular confinement if a comparable ansatz exists.
  • Real-time optical or chemical patterning could implement the required annulus and boundary changes in an experiment.
  • Multi-defect interactions might be managed by extending the controller to a vector of defect positions.
  • The method supplies a concrete testbed for whether biological systems already use analogous activity patterning to position defects during development.

Load-bearing premise

The ansatz for the polar field remains a valid approximation when activity is localized to an annulus and the boundary orientation is varied dynamically.

What would settle it

Apply the proportional-integral controller with a chosen target trajectory and record the defect position over time; if the measured path deviates substantially from the target, the control method fails.

Figures

Figures reproduced from arXiv: 2511.21359 by Abhinav Singh, Birte C. Geerds, Daniel J. G. Pearce, Frank J\"ulicher, Ivo F. Sbalzarini, Karsten Kruse, Mathieu Dedenon.

Figure 1
Figure 1. Figure 1: FIG. 1. Effect of heterogeneous activity patterning. [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Controlling defect motion with boundary angle. [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Dynamic feedback control of the defect location. [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
read the original abstract

Topological defects are a conspicuous feature of active liquid crystals that have been associated with important morphogenetic transitions in organismal development. Robust development thus requires a tight control of the motion and placement of topological defects. In this manuscript, we study a mechanism to control +1 topological defects in an active polar fluid confined to a disk. If activity is localized in an annulus within the disk, the defect moves on a circular trajectory around the center of the disk. Using an ansatz for the polar field, we determine the dependence of the angular speed and the circle radius on the boundary orientation of the polar field and the active annulus. Using a proportional integral controller, we guide the defect along complex trajectories by changing the active annulus size and the boundary orientation.

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

1 major / 0 minor

Summary. The manuscript studies control of +1 topological defects in a polar active fluid confined to a disk. When activity is localized to an annulus, the defect follows a circular trajectory. An ansatz for the polar field is introduced to derive analytic expressions for the defect's angular speed and orbital radius as functions of the active annulus size and boundary orientation. These relations are then used to implement a proportional-integral controller that dynamically adjusts the annulus radius and boundary angle to steer the defect along prescribed complex trajectories.

Significance. If the ansatz remains accurate under time-dependent parameter variation, the work demonstrates a concrete, analytically grounded route to spatiotemporal defect control in active polar fluids. This is relevant to biological morphogenesis where defect positioning is thought to matter. The combination of a reduced analytic model with a standard PI controller is a clear methodological strength; however, the absence of any reported comparison to full time-dependent simulations leaves the practical utility of the control law unconfirmed.

major comments (1)
  1. [Abstract (ansatz and controller sections)] The derivation of angular speed and radius (described in the abstract as obtained from an ansatz for the polar field) is used directly as the plant model for the PI controller. No numerical integration of the underlying time-dependent active polar fluid equations is reported to test whether the quasi-static ansatz remains valid when the annulus radius and boundary orientation are varied in real time by the controller. This validation step is load-bearing for the central claim that the controller successfully guides the defect along complex trajectories.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We are grateful to the referee for their thorough review and valuable feedback on our manuscript. We address the major comment point by point below.

read point-by-point responses

  1. Referee: [Abstract (ansatz and controller sections)] The derivation of angular speed and radius (described in the abstract as obtained from an ansatz for the polar field) is used directly as the plant model for the PI controller. No numerical integration of the underlying time-dependent active polar fluid equations is reported to test whether the quasi-static ansatz remains valid when the annulus radius and boundary orientation are varied in real time by the controller. This validation step is load-bearing for the central claim that the controller successfully guides the defect along complex trajectories.

    Authors: We agree that direct numerical validation of the quasi-static ansatz under time-varying control parameters would strengthen the manuscript. However, the present work is centered on deriving the analytic relations from the ansatz and demonstrating the feasibility of the PI control strategy within this reduced model. The ansatz is constructed to satisfy the boundary conditions and the steady-state structure of the polar field for fixed parameters. For the controller, we assume that the parameter variations are sufficiently slow for the quasi-static approximation to hold, which is a standard approach in reduced-order modeling for control design. Full time-dependent simulations of the active polar fluid equations with dynamic control are computationally intensive and beyond the scope of this theoretical study, but we have confirmed the ansatz accuracy for static cases. We can add a discussion of this assumption in the revised manuscript. revision: partial

Circularity Check

0 steps flagged

No significant circularity; ansatz-based derivation and standard controller are independent of inputs

full rationale

The paper states it uses an ansatz for the polar field to derive angular speed and radius as functions of annulus size and boundary orientation, then applies a standard PI controller to vary those parameters for trajectory control. No quoted equations reduce a claimed prediction to a fitted parameter by construction, no self-citations are invoked as load-bearing uniqueness theorems, and the controller is an external tool rather than a self-derived relation. The derivation chain therefore remains self-contained; any concern about ansatz validity under dynamics is a modeling assumption, not circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of a polar-field ansatz under localized activity and on the assumption that activity can be confined exactly to an annulus without additional hydrodynamic effects.

axioms (1)
  • domain assumption The polar field can be approximated by a specific ansatz when activity is localized to an annulus.
    Explicitly invoked in the abstract as the basis for determining angular speed and radius.

pith-pipeline@v0.9.0 · 5687 in / 1302 out tokens · 28791 ms · 2026-05-25T07:37:02.398788+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

36 extracted references · 36 canonical work pages · 1 internal anchor

  1. [1]

    The initial conditions forpare obtained by fixing the position of the defect at somer̸= 0 and minimizing the elastic energyF, Fig

    The boundary conditions impose the presence of a topological defect of charge +1 in the polar field. The initial conditions forpare obtained by fixing the position of the defect at somer̸= 0 and minimizing the elastic energyF, Fig. S1a and Sect. S IV A [28]. We then add a small random angular perturbation topwith a magnitude that vanishes at the boundary ...

    work page

  2. [2]

    The minimum results from two competing aspects: As the defect approaches the disk boundary, the gradient terms in the elastic energy increase

    Above this valuer min increases monotonically with ri, Fig 3b. The minimum results from two competing aspects: As the defect approaches the disk boundary, the gradient terms in the elastic energy increase. In contrast, as the defect enters the passive region, the term in the elastic energy proportional toχincreases, becauseP 0 is larger in the passive tha...

    work page 2023

  3. [3]

    Prost, F

    J. Prost, F. J¨ ulicher, and J. F. Joanny, Active gel physics, Nature Physics11, 111 (2015)

    work page 2015

  4. [4]

    Balasubramaniam, R.-M

    L. Balasubramaniam, R.-M. M` ege, and B. Ladoux, Ac- tive nematics across scales from cytoskeleton organiza- tion to tissue morphogenesis, Current Opinion in Genet- ics & Development73, 101897 (2022)

    work page 2022

  5. [5]

    Spontaneous motion in hierarchically assembled active matter

    T. Sanchez, D. T. N. Chen, S. J. DeCamp, M. Hey- mann, and Z. Dogic, Spontaneous motion in hierarchi- cally assembled active matter, Nature491, 431 (2012), 1301.1122

    work page internal anchor Pith review Pith/arXiv arXiv 2012

  6. [6]

    Gruler, U

    H. Gruler, U. Dewald, and M. Eberhardt, Nematic liquid crystals formed by living amoeboid cells, The European Physical Journal B - Condensed Matter and Complex Systems11, 187 (1999)

    work page 1999

  7. [7]

    Duclos, S

    G. Duclos, S. Garcia, H. G. Yevick, and P. Silberzan, Perfect nematic order in confined monolayers of spindle- shaped cells, Soft Matter10, 2346 (2014)

    work page 2014

  8. [8]

    T. B. Saw, A. Doostmohammadi, V. Nier, L. Kocgozlu, S. Thampi, Y. Toyama, P. Marcq, C. T. Lim, J. M. Yeo- mans, and B. Ladoux, Topological defects in epithelia govern cell death and extrusion, Nature544, 212 (2017)

    work page 2017

  9. [9]

    Kawaguchi, R

    K. Kawaguchi, R. Kageyama, and M. Sano, Topological defects control collective dynamics in neural progenitor cell cultures, Nature545, 327 (2017)

    work page 2017

  10. [10]

    Guillamat, C

    P. Guillamat, C. Blanch-Mercader, G. Pernollet, K. Kruse, and A. Roux, Integer topological defects or- ganize stresses driving tissue morphogenesis, Nature Ma- terials21, 588 (2022)

    work page 2022

  11. [11]

    Endresen, A

    K. Endresen, A. Murali, B. Geerds, D. J. G. Pearce, and F. Serra, Actuation of cell sheets in 3d (2025), arXiv:2411.17834 [cond-mat.soft]

    work page arXiv 2025

  12. [12]

    Maroudas-Sacks, L

    Y. Maroudas-Sacks, L. Garion, L. Shani-Zerbib, A. Livshits, E. Braun, and K. Keren, Topological de- fects in the nematic order of actin fibres as organization centres of Hydra morphogenesis, Nature Physics17, 251 (2021)

    work page 2021

  13. [13]

    Ravichandran, M

    Y. Ravichandran, M. Vogg, K. Kruse, D. J. G. Pearce, and A. Roux, Topology changes of Hydra define actin ori- entation defects as organizers of morphogenesis, Science Advances11, eadr9855 (2025)

    work page 2025

  14. [14]

    V. A. Webster-Wood, M. Guix, N. W. Xu, B. Behkam, H. Sato, D. Sarkar, S. Sanchez, M. Shimizu, and K. K. Parker, Biohybrid robots: recent progress, chal- lenges, and perspectives, Bioinspiration & Biomimetics 18, 015001 (2022)

    work page 2022

  15. [15]

    D. J. G. Pearce, S. Gat, G. Livne, A. Bernheim- Groswasser, and K. Kruse, Programming active meta- materials using topological defects, arXiv (2020), 2010.13141

    work page arXiv 2020

  16. [16]

    M. M. Norton, P. Grover, M. F. Hagan, and S. Fraden, Optimal Control of Active Nematics, Physical Review Letters125, 178005 (2020), 2007.14837

    work page arXiv 2020

  17. [17]

    Zhang, S

    R. Zhang, S. A. Redford, P. V. Ruijgrok, N. Kumar, A. Mozaffari, S. Zemsky, A. R. Dinner, V. Vitelli, Z. Bryant, M. L. Gardel, and J. J. d. Pablo, Spatiotem- poral control of liquid crystal structure and dynamics through activity patterning, Nature Materials20, 875 (2021), 1912.01630

    work page arXiv 2021

  18. [18]

    Shankar, L

    S. Shankar, L. V. D. Scharrer, M. J. Bowick, and M. C. Marchetti, Design rules for controlling active topological defects, Proceedings of the National Academy of Sciences 121, 10.1073/pnas.2400933121 (2024), 2212.00666

    work page doi:10.1073/pnas.2400933121 2024

  19. [19]

    Floyd, A

    C. Floyd, A. R. Dinner, and S. Vaikuntanathan, Tailor- ing interactions between active nematic defects with re- inforcement learning, Soft Matter21, 4488 (2025)

    work page 2025

  20. [20]

    Nishiyama, J

    K. Nishiyama, J. Berezney, M. M. Norton, A. Aggar- wal, S. Ghosh, M. F. Hagan, Z. Dogic, and S. Fraden, Closed-loop control of active nematic flows, arXiv 10.48550/arxiv.2408.14414 (2024), 2408.14414

    work page doi:10.48550/arxiv.2408.14414 2024

  21. [21]

    Ghosh, C

    S. Ghosh, C. Joshi, A. Baskaran, and M. F. Hagan, Spa- tiotemporal control of structure and dynamics in a polar 6 active fluid, Soft Matter20, 7059 (2024)

    work page 2024

  22. [22]

    Kruse, J.-F

    K. Kruse, J.-F. Joanny, F. J¨ ulicher, J. Prost, and K. Seki- moto, Asters, vortices, and rotating spirals in active gels of polar filaments., Phys. Rev. Lett.92, 078101 (2004)

    work page 2004

  23. [23]

    Kruse, J

    K. Kruse, J. F. Joanny, F. J¨ ulicher, J. Prost, and K. Seki- moto, Generic theory of active polar gels: a paradigm for cytoskeletal dynamics., The European physical journal. E, Soft matter16, 5 (2005)

    work page 2005

  24. [24]

    Schrader, S

    B. Schrader, S. Reboux, and I. F. Sbalzarini, Discretiza- tion correction of general integral PSE operators for par- ticle methods, Journal of Computational Physics229, 4159 (2010)

    work page 2010

  25. [25]

    Singh, A

    A. Singh, A. Foggia, P. Incardona, and I. F. Sbalzarini, A Meshfree Collocation Scheme for Surface Differential Op- erators on Point Clouds, Journal of Scientific Computing 96, 89 (2023)

    work page 2023

  26. [26]

    Incardona, A

    P. Incardona, A. Leo, Y. Zaluzhnyi, R. Ramaswamy, and I. F. Sbalzarini, OpenFPM: A scalable open framework for particle and particle-mesh codes on parallel comput- ers, Comput. Phys. Commun.241, 155 (2019)

    work page 2019

  27. [27]

    Singh, P

    A. Singh, P. Incardona, and I. F. Sbalzarini, A C++ expression system for partial differential equations en- ables generic simulations of biological hydrodynamics, Eur. Phys. J. E44, 117 (2021)

    work page 2021

  28. [28]

    Singh, P

    A. Singh, P. H. Suhrcke, P. Incardona, and I. F. Sbalzarini, A numerical solver for active hydrodynamics in three dimensions and its application to active turbu- lence, Phys. Fluids35, 105155 (2023)

    work page 2023

  29. [29]

    Singh, L

    A. Singh, L. Kraatz, S. Yaskovets, P. Incardona, and I. F. Sbalzarini, Integrating odeint time stepping into openfpm for distributed and gpu accelerated numerical solvers (2025), arXiv:2309.05331 [cs.MS]

    work page arXiv 2025

  30. [30]

    [18, 20, 22– 27, 29, 30]

    See supplemental material for additional figures, tables, calculations and computational details, [URL will be pro- vided by publisher], which includes Ref. [18, 20, 22– 27, 29, 30]

    work page

  31. [31]

    Giomi, M

    L. Giomi, M. J. Bowick, P. Mishra, R. Sknepnek, and M. Cristina Marchetti, Defect dynamics in active nemat- ics, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences372, 20130365 (2014)

    work page 2014

  32. [32]

    Di Leonardo, S

    R. Di Leonardo, S. Keen, F. Ianni, J. Leach, M. Pad- gett, and G. Ruocco, Hydrodynamic interactions in two dimensions, Physical Review E—Statistical, Nonlinear, and Soft Matter Physics78, 031406 (2008)

    work page 2008

  33. [33]

    H. P. Geering,Regelungstechnik: Mathematische Grund- lagen, Entwurfsmethoden, Beispiele, 6th ed. (Springer, 2013)

    work page 2013

  34. [34]

    Z. Zhao, H. Li, Y. Yao, Y. Zhao, F. Serra, K. Kawaguchi, H. Zhang, and M. Sano, Integer topological defects offer a methodology to quantify and classify active cell mono- layers, Nature Communications16, 2452 (2025)

    work page 2025

  35. [35]

    Coyle, B

    S. Coyle, B. Doss, Y. Huo, H. R. Singh, D. Quinn, K. Jimmy Hsia, and P. R. LeDuc, Cell alignment modu- lated by surface nano-topography – Roles of cell-matrix and cell-cell interactions, Acta Biomaterialia142, 149 (2022)

    work page 2022

  36. [36]

    K. D. Endresen, M. Kim, M. Pittman, Y. Chen, and F. Serra, Topological defects of integer charge in cell monolayers, Soft Matter17, 5878 (2021)

    work page 2021