pith. machine review for the scientific record. sign in

arxiv: 2605.04625 · v1 · submitted 2026-05-06 · 🧮 math.AP

Global well-posedness and decay rates for the three dimensional incompressible active liquid crystals

Pith reviewed 2026-05-08 17:27 UTC · model grok-4.3

classification 🧮 math.AP
keywords active liquid crystalsglobal well-posednessQ-tensordecay estimatesNavier-StokesBeris-Edwardsmixing decayincompressible flow
0
0 comments X

The pith

The 3D active liquid crystals system admits global strong solutions for small initial data when activity exceeds a critical value.

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

The paper proves global existence and uniqueness of strong solutions to the coupled incompressible Navier-Stokes equations and parabolic Q-tensor system that models active liquid crystals in three dimensions. Solutions exist globally for initial data small in H^{s+1} times H^s when the constant activity parameter c is strictly larger than a critical threshold c_star. The proof also yields decay estimates on derivatives of the Q-tensor that combine the algebraic decay of the heat kernel with an extra exponential decay whose rate grows with the excess activity (c minus c_star). These estimates remain valid in the vanishing rotational viscosity limit and imply that sufficiently active nematics lose orientational order and become isotropic at an activity-dependent rate.

Core claim

For the three-dimensional incompressible active liquid crystal equations with constant activity, unique global strong solutions exist for sufficiently small initial data (Q_0, u_0) in H^{s+1} times H^s when the activity constant c exceeds a critical value c_star. When the data are further in L^1 and s is at least 4, the Q-tensor satisfies mixing decay bounds ||partial^k Q(t)||_L2 that combine the optimal algebraic decay of the heat kernel with an exponential factor exp(-c' (c-c_star) Gamma t) for k up to s-1. Sharp decay estimates also hold for the velocity field under an additional integrability assumption on the data.

What carries the argument

Refined commutator estimates combined with the Green's function method and time-weighted energy estimates applied to the forced Navier-Stokes system coupled to the Beris-Edwards Q-tensor evolution.

Load-bearing premise

Initial data must be sufficiently small in the indicated Sobolev norms and activity must strictly exceed the critical value c_star; the equations are assumed to take the standard Beris-Edwards form with constant activity coefficient.

What would settle it

A construction or numerical example of finite-time singularity or non-uniqueness for arbitrarily small initial data when activity equals or falls below c_star would falsify the global well-posedness statement.

read the original abstract

This paper investigates the global well-posedness and large-time behavior of 3D incompressible active liquid crystals under constant activity, modeled by a coupled system of forced incompressible Navier-Stokes equations for the velocity and a parabolic system for the $Q$-tensor order parameter. By employing refined commutator estimates, the existence and uniqueness of global strong solutions are proved for small initial data $(Q_0,u_0)\in H^{s+1}\times H^s$ $(s\geq 2)$ with activity $c>c_\star$, which improves a previous result in \cite{active-limit}. In addition, if the initial data further belong to $L^1$ and $s\geq 4$, we obtain a mixing decay estimate on $\|\partial^kQ(t)\|_{L^2}$ that combines both an extra exponential decay factor at a rate proportional to $(c-c_\star)\Gamma$ and the optimal algebraic decay rate that coincides with that of the heat kernel, where $k\leq s-1$. This result reveals that, in the high activity regime, active nematics become isotropic with an activity-dependent exponential convergence rate, and the estimate is stable in the infinite rotational viscosity limit, as $\Gamma\rightarrow 0$. Meanwhile, the sharp decay estimate on $\|\partial^ku(t)\|_{L^2}$ is also derived for $k\leq s-2$ with an additional initial assumption. The proof is established via a combination of the Green's function method and the time-weighted energy method. To the best of our knowledge, these results are the first reported for active/passive nematic liquid crystals within the Beris-Edwards framework, and the enhanced decay effect of the orientational field is essentially derived from the free energy. Furthermore, in the passive setting, our result implies the phase transition of thermotropic liquid crystals at high temperatures.

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

0 major / 4 minor

Summary. The manuscript establishes global well-posedness and uniqueness of strong solutions to the 3D incompressible active liquid crystal system (Beris-Edwards) for sufficiently small initial data (Q0, u0) in H^{s+1} × H^s (s ≥ 2) provided the constant activity c exceeds an implicit threshold c_star. For data additionally in L1 and s ≥ 4 it derives mixing decay estimates on ||∂^k Q(t)||_L2 that combine the optimal 3D heat-kernel algebraic rate with an extra exponential factor proportional to (c - c_star)Γ (k ≤ s-1), together with sharp decay for the velocity field; the estimates remain stable as Γ → 0. The proofs combine refined commutator estimates, Green's function representations, and time-weighted energy methods, and the results are claimed to be the first of their kind for active/passive nematics in this framework.

Significance. If the a-priori estimates close, the work supplies the first global strong-solution theory for active nematics within the Beris-Edwards model and introduces a novel activity-enhanced decay mechanism traceable to free-energy dissipation. The stability of the decay under Γ → 0 and the implication for high-temperature phase transitions in the passive case are of independent interest. The techniques employed (commutator estimates, Green's functions, time-weighted energies) are standard for coupled parabolic-hyperbolic systems and appear internally consistent.

minor comments (4)
  1. The critical threshold c_star is introduced in the abstract and main statements but its explicit construction or lower bound is not displayed; a brief remark on how it arises from the energy dissipation would improve readability.
  2. The rotational viscosity Γ appears in the exponential rate without being defined in the model section; it should be introduced together with the constitutive relations.
  3. The improvement over the cited result in [active-limit] is asserted but the precise additional assumptions removed (or the relaxed smallness condition) are not itemized; a short comparison paragraph would clarify the advance.
  4. In the decay statements the range k ≤ s-1 (for Q) and k ≤ s-2 (for u) is given; the reason for the difference in the velocity index should be explained in one sentence.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We sincerely thank the referee for the careful reading and positive assessment of our manuscript. The referee's summary accurately reflects the main results on global well-posedness for small data when the activity exceeds the threshold c_star and the activity-enhanced decay estimates that remain stable as Γ → 0. We appreciate the recognition of the novelty within the Beris-Edwards framework and the interest in the implications for the passive case. We will incorporate all minor revisions as recommended.

Circularity Check

0 steps flagged

No significant circularity detected in derivation chain

full rationale

The paper derives global well-posedness for small data and the stated decay rates directly from a-priori estimates on the Beris-Edwards system, using refined commutator estimates, Green's function representations, and time-weighted energies. These close for sufficiently small initial data in the indicated Sobolev spaces when activity exceeds the threshold c_star, with the exponential factor arising explicitly from the model's free-energy dissipation. No step reduces by construction to a fitted parameter, self-referential definition, or load-bearing self-citation whose validity is assumed rather than independently verified. The citation to prior work is used only to contextualize the improvement, not to justify the core estimates. The analysis is self-contained against the PDE system and standard parabolic theory.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The analysis relies on classical functional-analytic tools for 3D Navier-Stokes-type systems; no new free parameters, ad-hoc axioms, or invented entities are introduced beyond the standard PDE framework.

axioms (2)
  • standard math Sobolev embedding and commutator estimates in three dimensions for H^s spaces with s >= 2
    Invoked to close the nonlinear estimates for the coupled velocity-Q system.
  • standard math Existence of the Green's function for the linearized parabolic system
    Used to obtain the optimal algebraic decay rates matching the heat kernel.

pith-pipeline@v0.9.0 · 5640 in / 1549 out tokens · 38111 ms · 2026-05-08T17:27:01.214040+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

69 extracted references · 2 canonical work pages · 1 internal anchor

  1. [1]

    Abels, G

    H. Abels, G. Dolzmann and Y . N. Liu,Well-posedness of a fully coupled Navier-Stokes/Q-tensor system with inhomogeneous boundary data, SIAM Journal on Mathematical Analysis, 2014, 46(4): 3050–3077

  2. [2]

    Abels, G

    H. Abels, G. Dolzmann and Y . N. Liu,Strong solutions for the Beris-Edwards model for nematic liquid crystals with homogeneous Dirichlet boundary conditions, Advances in Differential Equations, 2016, 21(1/2): 109–152

  3. [3]

    De Anna,A global 2D well-posedness result on the order tensor liquid crystal theory, Journal of Differential Equations, 2017, 262(7): 3932–3979

    F. De Anna,A global 2D well-posedness result on the order tensor liquid crystal theory, Journal of Differential Equations, 2017, 262(7): 3932–3979

  4. [4]

    De Anna and A

    F. De Anna and A. Zarnescu,Uniqueness of weak solutions of the full coupled Navier-Stokes and Q-tensor system in 2D, Communications in Mathematical Sciences, 2016 14(8): 2127–2178

  5. [5]

    J. M. Ball and A. Majumdar,Nematic liquid crystals: from Maier-Saupe to a continuum theory, Molecular crystals and liquid crystals, 2010, 525(1): 1–11

  6. [6]

    Ballerini, N

    M. Ballerini, N. Cabibbo, R. Candelier and et al.,Interaction ruling animal collective behavior de- pends on topological rather than metric distance: Evidence from a field study, Proceedings of the national academy of sciences, 2008, 105(4): 1232–1237

  7. [7]

    Bahouri, J

    H. Bahouri, J. Chemin and R. Danchin,Fourier Analysis and Nonlinear Partial Differential Equa- tions, Grundlehren der Mathematischen Wissenschaften, 2011, 343

  8. [8]

    M. L. Blow, S. P. Thampi and J. M. Yeomans,Biphasic, lyotropic, active nematics, Physical review letters, 2014, 113(24): 248303

  9. [9]

    Chandrasekhar,Liquid Crystals, second edition, Cambridge University Press, 1992

    S. Chandrasekhar,Liquid Crystals, second edition, Cambridge University Press, 1992

  10. [10]

    I. .L. Chern and T. P. Liu,Convergence to diffusion waves of solutions for viscous conservation laws, Communications in mathematical physics, 1987, 110(3): 503–517

  11. [11]

    Cavaterra, E

    C. Cavaterra, E. Rocca, H. Wu and X. Xu,Global Strong Solutions of the Full Navier-Stokes and Q- Tensor System for Nematic Liquid Crystal Flows in Two Dimensions, SIAM Journal on Mathematical Analysis, 2016, 48(2): 1368–1399

  12. [12]

    Chakrabarti and P

    N. Chakrabarti and P. Das,Isotropic to nematic phase transition in f-actin, Journal of Surface Science and Technology, 2007, 23(3/4): 177

  13. [13]

    G.-Q. Chen, A. Majumdar, D. Wang D and R. Zhang,Global weak solutions for the compressible active liquid crystal system, SIAM Journal on Mathematical Analysis, 2018, 50(4): 3632–3675

  14. [14]

    Y . Chen, D. Wang and R. Zhang,On mathematical analysis of complex fluids in active hydrodynam- ics, Electronic Research Archive, 2021, 29(6)

  15. [15]

    G.-Q. Chen, A. Majumdar, D. Wang and R. Zhang,Global existence and regularity of solutions for active liquid crystals, Journal of Differential Equations, 2017, 263(1): 202–239. 41

  16. [16]

    Bae, Y .-P

    H. Bae, Y .-P. Choi and K. Kang,Well-Posedness and Asymptotic Stability of Solutions for the Incom- pressible Toner–Tu Model, SIAM Journal on Mathematical Analysis, 2025, 57(1): 637–660

  17. [17]

    Y .-P. Choi, K. Kang and W. Lee,Global existence and asymptotic stability for the Toner-Tu model of flocking, Communications on Pure and Applied Analysis, 2025, 24(7): 1296–1321

  18. [18]

    M. E. Cates and E. Tjhung,Theories of binary fluid mixtures: from phase-separation kinetics to active emulsions, Journal of Fluid Mechanics, 2018, 836: P1

  19. [19]

    M. Dai, E. Feireisl, E. Rocca and et al.,On asymptotic isotropy for a hydrodynamic model of liquid crystals, Asymptotic Analysis, 2016, 97(3–4): 189–210

  20. [20]

    Doostmohammadi, J

    A. Doostmohammadi, J. Ignés-Mullol, J. M. Yeomans and et al.,Active nematics, Nature communi- cations, 2018, 9(1): 3246

  21. [21]

    N. C. Darnton, L. Turner, S. Rojevsky and H. C. Berg,Dynamics of bacterial swarming, Biophysical journal, 2010, 98(10): 2082–2090

  22. [22]

    H. Du, X. Hu and C. Wang,Suitable weak solutions for the co-rotational Beris-Edwards system in dimension three, Archive for Rational Mechanics and Analysis, 2020, 238(2): 749–803

  23. [23]

    Feireisl, E

    E. Feireisl, E. Rocca, G. Schimperna and A. Zarnescu, Evolution of non-isothermal Landau–de Gennes nematic liquid crystals flows with singular potential, Communications in Mathematical Sci- ences, 2014, 12(2): 317–343

  24. [24]

    Feireisl, G

    E. Feireisl, G. Schimperna, E. Rocca and A. Zarnescu,Nonisothermal nematic liquid crystal flows with the Ball-Majumdar free energy, Annali di Matematica Pura ed Applicata, 2015, 194(5): 1269– 1299

  25. [25]

    Giomi, L

    L. Giomi, L. Mahadevan, B. Chakraborty and M. F. Hagan,Excitable patterns in active nematics, Physical review letters, 2011, 106(21): 218101

  26. [26]

    Giomi, L

    L. Giomi, L. Mahadevan, B. Chakraborty and M. F. Hagan,Banding, excitability and chaos in active nematic suspensions, Nonlinearity, 2012, 25(8): 2245

  27. [27]

    Giomi, M

    L. Giomi, M. J. Bowick, X. Ma and M. C. Marchetti,Defect annihilation and proliferation in active nematics, Physical review letters, 2013, 110(22): 228101

  28. [28]

    Giomi, M

    L. Giomi, M. J. Bowick, P. Mishra and et al.,Defect dynamics in active nematics, Philosophi- cal Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2014, 372(2029): 20130365

  29. [29]

    Guillén-González and M Á

    F. Guillén-González and M Á. Rodríguez-Bellido,Weak solutions for an initial-boundary Q-tensor problem related to liquid crystals, Nonlinear Analysis: Theory, Methods & Applications, 2015, 112: 84–104

  30. [30]

    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 Systems, 1999, 11: 187–192. 42

  31. [31]

    P. G. De Gennes and J. Prost,The physics of liquid crystals, second edition, Oxford university press, 1995

  32. [32]

    Hieber, A

    M. Hieber, A. Hussein and M. Wrona,Strong Well-Posedness of the Q-Tensor Model for Liquid Crystals: The Case of Arbitrary Ratio of Tumbling and Aligning Effectsξ, Archive for Rational Mechanics and Analysis, 2024, 248(3): 40

  33. [33]

    J. R. Huang and S. J. Ding,Global well-posedness for the dynamical Q-tensor model of liquid crys- tals, Science China Mathematics, 2015, 58: 1349–1366

  34. [34]

    Huang, Y

    J. Huang, Y . Wang, H. Wen and R. Zi, Optimal time-decay estimates for an Oldroyd-B model with zero viscosity, Journal of Differential Equations, 2022, 306: 456–491

  35. [35]

    Huang, Q

    J. Huang, Q. Liu and R. Zi,Global existence and decay rates of solutions to the Oldroyd-B model with stress tensor diffusion, Journal of Differential Equations, 2024, 389: 38–89

  36. [36]

    Jiang, L

    N. Jiang, L. Ke and X. Song,Well-posedness for the compressible active liquid crystal model, Journal of Mathematics (PRC), 2023, 43(2): 95–125

  37. [37]

    Jiang, S

    N. Jiang, S. Tang and B. Wang,Well-posedness and the zero activity limit for the active nematic liquid crystal model, Discrete and Continuous Dynamical Systems-B, 2024, 29(1): 68–109

  38. [38]

    Klainerman and A

    S. Klainerman and A. Majda,Singular limits of quasilinear hyperbolic systems with large parameters and the incompressible limit of compressible fluids, Communications on pure and applied Mathemat- ics, 1981, 34(4): 481–524

  39. [39]

    I. C. Khoo and S. T. Wu,Optics and nonlinear optics of liquid crystals, world scientific, 1993

  40. [40]

    Kagei and T

    Y . Kagei and T. Kobayashi,On large-time behavior of solutions to the compressible Navier-Stokes equations in the half space inR 3, Archive for Rational Mechanics and Analysis, 2002, 165(2): 89– 159

  41. [41]

    Kagei and T

    Y . Kagei and T. Kobayashi,Asymptotic behavior of solutions of the compressible Navier-Stokes equa- tions on the half space, Archive for Rational Mechanics and Analysis, 2005, 177(2): 231–330

  42. [42]

    Kierfeld, K

    J. Kierfeld, K. Frentzel, P. Kraikivski and et al.,Active dynamics of filaments in motility assays, The European Physical Journal Special Topics, 2008, 157: 123–133

  43. [43]

    Lian and R

    W. Lian and R. Zhang,Global weak solutions to the active hydrodynamics of liquid crystals, Journal of Differential Equations, 2020, 268(8): 4194–4221

  44. [44]

    T. B. Liverpool and M. C. Marchetti,Hydrodynamics and rheology of active polar filaments, In: Cell Motility. Biological and Medical Physics, Biomedical Engineering, Springer New York, NY , 2008: 177–206

  45. [45]

    Majda,Compressible fluid flow and systems of conservation laws in several space variables, Springer Science & Business Media, 1984

    A. Majda,Compressible fluid flow and systems of conservation laws in several space variables, Springer Science & Business Media, 1984

  46. [46]

    M. C. Marchetti, J. F. Joanny, S. Ramaswamy and et al.,Hydrodynamics of soft active matter, Re- views of modern physics, 2013, 85(3): 1143–1189. 43

  47. [47]

    Schonbek and Y

    M. Schonbek and Y . Shibata,Global well-posedness and decay for a Q tensor model of incompress- ible nematic liquid crystals inR N, Journal of Differential Equations, 2019, 266(6): 3034–3065

  48. [48]

    Matsumura,On the asymptotic behavior of solutions of semi-linear wave equations, Publications of the Research Institute for Mathematical Sciences, 1976, 12(1): 169–189

    A. Matsumura,On the asymptotic behavior of solutions of semi-linear wave equations, Publications of the Research Institute for Mathematical Sciences, 1976, 12(1): 169–189

  49. [49]

    Murata and Y

    M. Murata and Y . Shibata,Global well posedness for a Q-tensor model of nematic liquid crystals, Journal of Mathematical Fluid Mechanics, 2022, 24(2): 34

  50. [50]

    Nishida,Nonlinear hyperbolic equations and related topics in fluid dynamics, Publications Math- ematiques D’Orsay 78.02, 1978: 46–53

    T. Nishida,Nonlinear hyperbolic equations and related topics in fluid dynamics, Publications Math- ematiques D’Orsay 78.02, 1978: 46–53

  51. [51]

    T. J. Pedley and J. O. Kessler,Hydrodynamic phenomena in suspensions of swimming microorgan- isms, Annual Review of Fluid Mechanics, 1992, 24(1): 313–358

  52. [52]

    Paicu and A

    M. Paicu and A. Zarnescu,Global existence and regularity for the full coupled Navier–Stokes and Q-tensor system, SIAM journal on mathematical analysis, 2011, 43(5): 2009–2049

  53. [53]

    Paicu and A

    M. Paicu and A. Zarnescu,Energy dissipation and regularity for a coupled Navier–Stokes and Q- tensor system, Archive for Rational Mechanics and Analysis, 2012, 203: 45–67

  54. [54]

    Qiu and Y

    Z. Qiu and Y . Wang,Martingale solution for stochastic active liquid crystal system, Discrete and Continuous Dynamical Systems, 2021, 41(5): 2227–2268

  55. [55]

    Ravnik and J

    M. Ravnik and J. M. Yeomans,Confined active nematic flow in cylindrical capillaries, Physical review letters, 2013, 110(2): 026001

  56. [56]

    Sanchez, D

    T. Sanchez, D. T. N. Chen, S. J. DeCamp and et al.,Spontaneous motion in hierarchically assembled active matter, Nature, 2012, 491(7424): 431–434

  57. [57]

    C. K. Schmidt, M. Medina-Sánchez, R. J. Edmondson and et al.,Engineering microrobots for tar- geted cancer therapies from a medical perspective, Nature Communications, 2020, 11(1): 5618

  58. [58]

    Toner and Y

    J. Toner and Y . Tu,Long-range order in a two-dimensional dynamical XY model: how birds fly together, Physical review letters, 1995, 75(23): 4326

  59. [59]

    Wilkinson,Strictly physical global weak solutions of a Navier-Stokes Q-tensor system with sin- gular potential, Archive for Rational Mechanics and Analysis, 2015, 218(1): 487–526

    M. Wilkinson,Strictly physical global weak solutions of a Navier-Stokes Q-tensor system with sin- gular potential, Archive for Rational Mechanics and Analysis, 2015, 218(1): 487–526

  60. [60]

    D. Wang, X. Xu and C. Yu,Global weak solution for a coupled compressible Navier-Stokes and Q-tensor system, Communications in Mathematical Sciences, 2015, 13(1): 49–82

  61. [61]

    M. Y . Wong, C. Y . Tso, T. C. Ho and H. H. Lee,A review of state of the art thermal diodes and their potential applications, International Journal of Heat and Mass Transfer, 2021, 164: 120607

  62. [62]

    Xiao,Global strong solution to the three-dimensional liquid crystal flows of Q-tensor model, Journal of Differential Equations, 2017, 262(3): 1291–1316

    Y . Xiao,Global strong solution to the three-dimensional liquid crystal flows of Q-tensor model, Journal of Differential Equations, 2017, 262(3): 1291–1316

  63. [63]

    Yang and X

    F. Yang and X. Yang,Global well-posedness and decay estimates for the one-dimensional models of blood flow with a general parabolic velocity profile, Nonlinear Analysis: Real World Applications, 2024, 78: 104098. 44

  64. [64]

    Yang and C

    F. Yang and C. Li,Weak-strong uniqueness for three dimensional incompressible active liquid crys- tals, Acta Mathematica Scientia, 2024: 1–26

  65. [65]

    Yang and J

    F. Yang and J. Zhou,Weak-strong uniqueness of the full coupled Navier-Stokes and Q-tensor system in dimension three, arXiv preprint arXiv:2507.09281, 2025

  66. [66]

    Koch-Tataru theorem for 3D incompressible active nematic liquid crystals

    F. Yang,Koch-Tataru theorem for 3D incompressible active nematic liquid crystals, arXiv preprint arXiv:2604.04000, 2026

  67. [67]

    R. Zi, D. Fang and T. Zhang,Global solution to the incompressible Oldroyd-B model in the critical Lp framework: the case of the non-small coupling parameter, Archive for Rational Mechanics and Analysis, 2014, 213: 651–687

  68. [68]

    M. C. Zelati,Stable mixing estimates in the infinite Péclet number limit, Journal of Functional Anal- ysis, 2020, 279(4): 108562

  69. [69]

    Zhang, A

    R. Zhang, A. Mozaffari and J. J. de Pablo,Logic operations with active topological defects, Science advances, 2022, 8(8): eabg9060. 45