pith. sign in

arxiv: 2512.20289 · v2 · submitted 2025-12-23 · ❄️ cond-mat.soft · cond-mat.stat-mech· physics.bio-ph· physics.comp-ph

Spatiotemporal Chaos and Defect Proliferation in Polar-Apolar Active Mixture

Partha Sarathi Mondal , Tamas Vicsek , Shradha Mishra This is my paper

Pith reviewed 2026-05-16 20:21 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mechphysics.bio-phphysics.comp-ph
keywords spatiotemporal chaosactive mattertopological defectspolar-apolar mixturehydrodynamic equationsband structuresLyapunov exponentdefect proliferation
0
0 comments X

The pith

In a polar-apolar active particle mixture, an intermediate regime produces a chaotic phase of evolving high-density bands with continual creation and annihilation of half-integer topological defects.

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

The paper solves the hydrodynamic equations for a mixture of self-propelled polar and apolar particles. It finds a non-monotonic response of the apolar component to changes in polar density and activity. At intermediate values a disordered state appears with high-density bands that evolve chaotically and with ongoing creation and annihilation of half-integer defects. The authors quantify the resulting spatiotemporal chaos through the spectrum of density fluctuations and the maximal Lyapunov exponent. This regime is positioned as richer than the disorder seen when the apolar species is passive.

Core claim

Numerical solution of the coupled hydrodynamic equations for the polar-apolar mixture reveals that, in an intermediate regime of polar activity and density, the system enters a dynamically disordered phase marked by high-density, chaotically evolving band-like structures together with continuous creation and annihilation of half-integer topological defects; this phase is spatiotemporally chaotic, as shown by the spectral properties of density fluctuations and a positive maximal Lyapunov exponent.

What carries the argument

The coupled hydrodynamic equations for the densities and order parameters of the polar and apolar components, solved numerically across parameter ranges to identify the intermediate chaotic regime.

If this is right

  • The apolar species exhibits a non-monotonic response to polar density and activity.
  • Spatiotemporal chaos appears through an inverse energy cascade typical of active systems rather than a direct cascade.
  • Continual defect creation and annihilation occurs together with chaotic band evolution.
  • The observed states broaden the set of possible complex transitions beyond those in living liquid crystal systems.

Where Pith is reading between the lines

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

  • The chaotic regime could be used to control defect statistics by tuning only polar activity in an experimental mixture.
  • Similar band-and-defect chaos may appear in other active mixtures such as colloidal rods combined with swimming bacteria.
  • The Lyapunov exponent provides a practical diagnostic that could be tracked in future numerical studies of active mixtures to locate the onset of chaos.

Load-bearing premise

The chosen hydrodynamic equations fully capture the microscopic interactions and couplings between the polar and apolar components without requiring extra higher-order terms.

What would settle it

Direct measurement of the maximal Lyapunov exponent remaining zero (instead of positive) in either simulation or experiment at the reported intermediate polar densities and activities would falsify the spatiotemporal-chaos claim.

Figures

Figures reproduced from arXiv: 2512.20289 by Partha Sarathi Mondal, Shradha Mishra, Tamas Vicsek.

Figure 1
Figure 1. Figure 1: Global characteristics of active nematics containing an active polar component. Panel (a) shows the phase diagram of the system based on the steady-state characteristics of the active nematic. In the homogeneous regime, both ρn and Q are homogeneous, whereas in the IN-regime ρn and Q exhibit significant spatial fluctuations. The insets show the probability density function of Q = |Q|, P(Q). In the homogene… view at source ↗
Figure 2
Figure 2. Figure 2: Snapshots of the magnitude of Q-field (i.e. Q) for the apolar species for ρp0 = 0.09 in Phase-II for different values of polar activity : (a) vp = 0.001, (b) vp = 0.005, (c) vp = 0.05, (d) vp = 0.10, (e) vp = 0.20, and (f) vp = 0.50. The heatmap depicts the magnitude of the nematic order parameter field of active nematics Q = |Q|. Parameters: System size, L = 512. The rest of the parameters are the same as… view at source ↗
Figure 3
Figure 3. Figure 3: Variation of structural properties of the bands with changing control parameters in the IN regime. Panel (a) shows the variation of averaged correlation length, < lρn >, with polar density, ρp0, in Phase-II for two different system sizes L = 512, & 1024 for vp = 0.25. The inset shows the time series of lρn (t) for vp = 0.25 and ρp = 0.07; Panel (b) presents lρn vs. ρp0 plot for different values of vp. The … view at source ↗
Figure 4
Figure 4. Figure 4: Plot of the autocorrelation function of the fluctuations in lρn (t), Cacf (t),vs.time, t, for a different set of control parameters. Panel (a) shows Cacf (t)vs.t for different values of ρp0 for vp = 0.20. The inset shows the plot of the correlation time, τc, vs. ρp0. Panel (b) shows Cacf (t)vs.t for different values of vp for ρp0 = 0.07. The inset shows the plot of the correlation time, τc, vs. vp. The err… view at source ↗
Figure 5
Figure 5. Figure 5: Mechanism of the formation of chaotic bands following an instantaneous change in the vp. The system is prepared to have an initial condition that corresponds to the steady-state for vp = 0.001, as shown in panel (a), which shows the configuration of the nematic order parameter field |Q(r)|. The system is then instantaneously changed to vp = 0.07. The subsequent snapshots in the top row (b-e) show the tempo… view at source ↗
Figure 6
Figure 6. Figure 6: Variation of the mean density fluctuation of the polar species, ∆ρp (in panel (a)) and mean stress, σ (in panel (b)) with time. Parameters : L = 400. The rest of the parameters are the same as in figure 1. 3.2.2 Mechanism of formation of bands : Before closing this section, we look into the mechanism of formation of modulating bands at high activity of polar species. At relatively lower vp, the steady stat… view at source ↗
Figure 7
Figure 7. Figure 7: Visual description of the mechanism of formation of defects. The panels (a)-(e) present a zoomed-in view of the system at subsequently increasing times. The colour represents the magnitude of the local nematic order parameter (Q, given by Eq.5) according to the colour bar, and the lines represent the orientations of the local nematic director. The panels (f)-(j) depict the plot of the stress (∇Qxy) at time… view at source ↗
Figure 8
Figure 8. Figure 8: Variation of mean number of 1/2-integer topological defects in the steady state of the system, ⟨nd⟩, on tuning the control parameters. Panel (a) shows the plot ⟨nd⟩ vs. ρp0 for different vp; Panel (b) shows the plot of ⟨nd⟩vs.vp for two different ρp0. The error bar represents the standard deviation calculated over time in steady state as well as independent realisations; Panel (c) presents the phase diagra… view at source ↗
Figure 9
Figure 9. Figure 9: Characteristics of the frequency spectrum of the fluctuations in lρn (t) in the steady state of the system. Panel (a) illustrates the piecewise stationary characteristics of the frequency spectrum. The frequency spectrum exhibits exponential characteristics at low frequencies and power-law characteristics at high frequencies; Panel (b) showcases the effect of a finite number of points in the time series on… view at source ↗
Figure 10
Figure 10. Figure 10: Variation of maximal Lyapunov exponent (MLE) on control parameters (ρp0, vp) for both NT and T S methods. Panel (a) shows the variation of MLE with ρp0 for a fixed vp. Panel (b) shows the variation of MLE with vp at a fixed ρp0. The error bars in (a) and (b) represent the standard deviation calculated over 10 independent realisations. Thus, the comprehensive analysis in this section establishes that the d… view at source ↗
read the original abstract

Chaotic transitions in inertial fluids typically proceed through a direct energy cascade from large to small scales. In contrast, active systems, composed of self propelled units, inject energy at microscopic scales and therefore exhibit an inverse cascade, giving rise to distinctly unconventional flow patterns. Here, we investigate an active mixture consisting of both apolar and polar self driven components, a setting expected to display richer behaviours than those found in living liquid crystal (LLC) systems, where the apolar constituent is passive. Using numerical solutions of the corresponding hydrodynamic equations, we uncover a variety of complex dynamical states. Our results reveal a non-monotonic response of the apolar species to changes in the density and activity of the polar component. In an intermediate regime, reminiscent of LLC-induced disorder, the system develops a dynamically disordered phase characterised by high-density, chaotically evolving band-like structures and by the continual creation and annihilation of half integer topological defects. We show that this regime exhibits spatiotemporal chaos, which we quantify through two complementary measures: the spectral properties of density fluctuations and the maximal Lyapunov exponent. Together, these findings broaden the understanding of complex transitions in active matter and suggest potential experimental realisations in bacterial suspensions or synthetic microswimmer assemblies.

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 numerically solves hydrodynamic equations for a polar-apolar active mixture and reports a non-monotonic response of the apolar component to polar density and activity. In an intermediate regime the system enters a disordered phase with high-density chaotically evolving bands, continual creation/annihilation of half-integer defects, and spatiotemporal chaos diagnosed by density-fluctuation spectra and a positive maximal Lyapunov exponent.

Significance. If the numerical diagnostics prove robust under refinement, the work extends chaos studies in active matter beyond passive-apolar LLCs by demonstrating defect proliferation and inverse-cascade chaos in a fully active mixture, with potential experimental relevance to bacterial or microswimmer suspensions.

major comments (2)
  1. [Numerical results (Lyapunov exponent subsection)] Numerical results section on Lyapunov exponent: the reported positive maximal Lyapunov exponent is presented without grid-refinement studies, domain-size scaling, or checks against artificial viscosity and perturbation norm. In continuum active-matter PDEs the sign of this exponent is known to be sensitive to spatial discretization; absence of these controls leaves open the possibility that the observed chaos is a numerical artifact of under-resolved inverse cascades.
  2. [Results on density spectra] Section on spectral properties of density fluctuations: the power spectra are shown without ensemble averaging, error bars, or explicit convergence tests with respect to system size and integration time. This weakens the claim that the spectra independently confirm spatiotemporal chaos, especially given the non-monotonic parameter dependence highlighted in the abstract.
minor comments (2)
  1. [Abstract] Abstract: the statement that the disordered phase is 'reminiscent of LLC-induced disorder' would benefit from a one-sentence quantitative comparison to the relevant LLC literature (e.g., defect density or spectral exponent).
  2. [Figure captions] Figure captions: several panels lack explicit parameter values or grid resolution, making it difficult for readers to reproduce the reported states.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed report. We address each major comment below and will revise the manuscript to incorporate additional numerical controls that strengthen the evidence for spatiotemporal chaos.

read point-by-point responses

  1. Referee: [Numerical results (Lyapunov exponent subsection)] Numerical results section on Lyapunov exponent: the reported positive maximal Lyapunov exponent is presented without grid-refinement studies, domain-size scaling, or checks against artificial viscosity and perturbation norm. In continuum active-matter PDEs the sign of this exponent is known to be sensitive to spatial discretization; absence of these controls leaves open the possibility that the observed chaos is a numerical artifact of under-resolved inverse cascades.

    Authors: We agree that systematic convergence tests are necessary to exclude numerical artifacts in continuum models of active matter. In the revised manuscript we will add grid-refinement studies (halving the grid spacing), domain-size scaling (doubling the system size while keeping resolution fixed), and explicit checks varying artificial viscosity and perturbation norm. Additional simulations performed since submission confirm that the maximal Lyapunov exponent remains positive and quantitatively stable under these refinements, indicating that the reported chaos is not an artifact of under-resolution. revision: yes

  2. Referee: [Results on density spectra] Section on spectral properties of density fluctuations: the power spectra are shown without ensemble averaging, error bars, or explicit convergence tests with respect to system size and integration time. This weakens the claim that the spectra independently confirm spatiotemporal chaos, especially given the non-monotonic parameter dependence highlighted in the abstract.

    Authors: We acknowledge that ensemble averaging and convergence diagnostics would make the spectral evidence more robust. In the revised version we will include ensemble averages over at least five independent realizations, with error bars, together with explicit tests demonstrating convergence of the spectral shape with increasing system size and integration time. These additions will confirm that the inverse-cascade features and the non-monotonic dependence on polar density and activity are reproducible and not sensitive to finite-size or finite-time effects. revision: yes

Circularity Check

0 steps flagged

Numerical study applies independent diagnostics to simulated fields with no self-referential derivation

full rationale

The paper solves a set of hydrodynamic PDEs numerically for varying parameters and directly computes two independent diagnostics (Fourier spectra of density fluctuations and the maximal Lyapunov exponent) on the resulting spatiotemporal fields. No step fits a parameter to a data subset and then renames the output as a prediction, no self-citation supplies a load-bearing uniqueness theorem, and no ansatz or definition is smuggled in via prior work by the same authors. The central claim of an intermediate chaotic regime follows from the observed behavior of the simulated solutions rather than from any tautological reduction to the input equations or fitted quantities.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard hydrodynamic description of active matter with added polar-apolar coupling terms. Several activity and density parameters are chosen to place the system in the reported regime; the hydrodynamic approximation itself is taken as given.

free parameters (2)
  • polar activity coefficient
    Tuned to explore the intermediate regime where non-monotonic response and chaos appear.
  • apolar density and activity parameters
    Varied to demonstrate non-monotonic dependence on polar component.
axioms (1)
  • domain assumption Hydrodynamic equations with polar-apolar couplings accurately describe the mixture at the chosen scales.
    Invoked by choosing to solve those equations numerically without additional microscopic validation.

pith-pipeline@v0.9.0 · 5531 in / 1397 out tokens · 26892 ms · 2026-05-16T20:21:52.897782+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.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    The evolution equations for the density and the symmetry-broken variable for both species are discussed below. ... ∂tρp = ∇·(Dρp ∇ρp − vp ρp P) ... ∂tQij = −ΓQ ∂FQ/∂Q + ... + γ(Pi Pj − ½ δij Pk Pk)

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

60 extracted references · 60 canonical work pages

  1. [1]

    Vicsek T, Czir´ ok A, Ben-Jacob E, Cohen I and Shochet O 1995Physical Review Letters75 1226

    work page

  2. [2]

    Bechinger C, Di Leonardo R, L¨ owen H, Reichhardt C, Volpe G and Volpe G 2016Reviews of Modern Physics88045006

    work page

  3. [3]

    B¨ ar M, Großmann R, Heidenreich S and Peruani F 2020Annual Review of Condensed Matter Physics11441–466

    work page

  4. [4]

    thesis Indian Institute of Technology Hyderabad

    Ganguly C 2013Aspects of Self-PropulsionPh.D. thesis Indian Institute of Technology Hyderabad

    work page

  5. [5]

    Vicsek T and Zafeiris A 2012Physics reports51771–140

    work page

  6. [6]

    4th edition(Garland Science)

    Alberts B, Johnson A, Lewis J, Raff M, Roberts K and Walter P 2002 The self-assembly and dynamic structure of cytoskeletal filamentsMolecular Biology of the Cell. 4th edition(Garland Science)

    work page 2002

  7. [7]

    Kolomeisky A B and Fisher M E 2007Annu. Rev. Phys. Chem.58675–695

    work page

  8. [8]

    Niu B, Wang Het al.2012Discrete Dynamics in Nature and Society2012

    work page

  9. [9]

    Shapiro J A 1995BioEssays17597–607

    work page

  10. [10]

    Jolles J W, Boogert N J, Sridhar V H, Couzin I D and Manica A 2017Current Biology27 2862–2868

    work page

  11. [11]

    Moreno Garc´ ıa C A, Maxwell T M, Hickford J and Gregorini P 2020Frontiers in Veterinary Science774

    work page

  12. [12]

    Gueron S, Levin S A and Rubenstein D I 1996Journal of Theoretical Biology18285–98

    work page

  13. [13]

    Wang X and Lu J 2019IEEE Circuits and Systems Magazine196–22

    work page

  14. [14]

    Sokolov A and Aranson I S 2012Physical review letters109248109 20

    work page

  15. [15]

    Dunkel J, Heidenreich S, Drescher K, Wensink H H, B¨ ar M and Goldstein R E 2013Physical review letters110228102

    work page

  16. [16]

    Section A

    Gruler H, Schienbein M, Franke K and Boisfleury-chevance A d 1995Molecular Crystals and Liquid Crystals Science and Technology. Section A. Molecular Crystals and Liquid Crystals 260565–574

    work page

  17. [17]

    Zhou S, Sokolov A, Lavrentovich O D and Aranson I S 2014Biophysical Journal106420a

    work page

  18. [18]

    Turiv T, Koizumi R, Thijssen K, Genkin M M, Yu H, Peng C, Wei Q H, Yeomans J M, Aranson I S, Doostmohammadi Aet al.2020Nature Physics16481–487

    work page

  19. [19]

    Kumar A, Galstian T, Pattanayek S K and Rainville S 2013Molecular Crystals and Liquid Crystals57433–39

    work page

  20. [20]

    Doostmohammadi A, Shendruk T N, Thijssen K and Yeomans J M 2017Nature communications815326

    work page

  21. [21]

    de la Cotte A, Pearce D J, Nambisan J, Puggioni L, Levy A, Giomi L and Fernandez-Nieves A 2025Proceedings of the National Academy of Sciences122e2512147122

    work page

  22. [22]

    Shankar S, Souslov A, Bowick M J, Marchetti M C and Vitelli V 2022Nature Reviews Physics 4380–398

    work page

  23. [23]

    Alert R, Casademunt J and Joanny J F 2022Annual Review of Condensed Matter Physics13 143–170

    work page

  24. [24]

    Genkin M M, Sokolov A, Lavrentovich O D and Aranson I S 2017Physical Review X7011029

    work page

  25. [25]

    Vats A, Yadav P K, Banerjee V and Puri S 2023Physical Review E108024701

    work page

  26. [26]

    Kemkemer R, Kling D, Kaufmann D and Gruler H 2000The European Physical Journal E1 215–225

    work page

  27. [27]

    Balasubramaniam L, M` ege R M and Ladoux B 2022Current Opinion in Genetics & Development73101897

    work page

  28. [28]

    Saw T B, Xi W, Ladoux B and Lim C T 2018Advanced materials301802579

    work page

  29. [29]

    Sampat P B and Mishra S 2021Physical Review E104024130

    work page

  30. [30]

    Mondal P S, Mishra P K, Vicsek T and Mishra S 2025Physica A: Statistical Mechanics and its Applications659130338 ISSN 0378-4371 URL https://www.sciencedirect.com/science/article/pii/S0378437124008483

    work page

  31. [31]

    Sharma A and Soni H 2025Soft Matter

    work page

  32. [32]

    De Gennes P G and Prost J 1993The physics of liquid crystals83 (Oxford university press)

    work page

  33. [33]

    Brochard F and De Gennes P 1970Journal de Physique31691–708

    work page

  34. [34]

    Lopatina L M and Selinger J V 2009Physical review letters102197802

    work page

  35. [35]

    Bisht K, Wang Y, Banerjee V and Majumdar A 2020Physical Review E101022706

    work page

  36. [36]

    Aditi Simha R and Ramaswamy S 2002Physical review letters89058101

    work page

  37. [37]

    Ramaswamy S, Simha R A and Toner J 2003Europhysics Letters62196

    work page

  38. [38]

    Bertin E, Chat´ e H, Ginelli F, Mishra S, Peshkov A and Ramaswamy S 2013New journal of physics15085032

    work page

  39. [39]

    Goldenfeld N 2018Lectures on phase transitions and the renormalization group(CRC Press)

    work page

  40. [40]

    Toner J 2022Active Matter and Non-Equilibrium Statistical Physics: Lecture Notes of the Les Houches Summer School, September 2018 112, eds. G52–101

    work page 2018

  41. [41]

    Hohenberg P and Shraiman B I 1989Physica D: Nonlinear Phenomena37109–115 21

    work page

  42. [42]

    Wenzel D, Nestler M, Reuther S, Simon M and Voigt A 2021Computational Methods in Applied Mathematics21683–692

    work page

  43. [43]

    Delmarcelle T and Hesselink L 1994 The topology of symmetric, second-order tensor fields Proceedings Visualization’94(IEEE) pp 140–147

    work page 1994

  44. [44]

    Cross M and Hohenberg P 1994Science2631569–1570

    work page

  45. [45]

    Packard N H, Crutchfield J P, Farmer J D and Shaw R S 1980Physical review letters45712

    work page

  46. [46]

    Lai Y C and Ye N 2003International Journal of Bifurcation and Chaos131383–1422

    work page

  47. [47]

    Skokos C H, Gottwald G A and Laskar J 2016Chaos detection and predictabilityvol 1 (Springer)

    work page

  48. [48]

    Berera A and Ho R D 2018Physical review letters120024101

    work page

  49. [49]

    Boffetta G and Musacchio S 2017Physical review letters119054102

    work page

  50. [50]

    Valsakumar M, Satyanarayana S and Sridhar V 1997Pramana4869–85

    work page

  51. [51]

    Marwan N and Braun T 2023Chaos: An Interdisciplinary Journal of Nonlinear Science33

    work page

  52. [52]

    Adenij A E 2024Communication in Physical Sciences11

    work page

  53. [53]

    Mishra S, Baskaran A and Marchetti M C 2010Physical Review E—Statistical, Nonlinear, and Soft Matter Physics81061916

    work page

  54. [54]

    Marchetti M C, Joanny J F, Ramaswamy S, Liverpool T B, Prost J, Rao M and Simha R A 2013Reviews of modern physics851143–1189

    work page

  55. [55]

    Mishra P K, Mondal P S, Jena P and Mishra S 2025New Journal of Physics27074602 URL https://doi.org/10.1088/1367-2630/adef71

    work page doi:10.1088/1367-2630/adef71

  56. [56]

    Kantz H and Schreiber T 2003Nonlinear time series analysis(Cambridge university press)

    work page

  57. [57]

    Casdagli M 1992 Nonlinear forecasting, chaos and statisticsModeling Complex Phenomena: Proceedings of the Third Woodward Conference, San Jose State University, April 12–13, 1991 (Springer) pp 131–152

    work page 1992

  58. [58]

    Fraser A M and Swinney H L 1986Physical review A331134

    work page

  59. [59]

    Kennel M B, Brown R and Abarbanel H D 1992Physical review A453403

    work page

  60. [60]

    Kaplan D T and Glass L 1992Physical review letters68427 22

    work page