Pith. sign in

REVIEW 2 major objections 4 minor 40 references

Thin active nematohydrodynamic layers: asymptotic theories and instabilities

T0 review · 2 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash

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

Pith's one-line read Allowing a film to bend and change thickness flips the sign rule for activity-driven ordering: on a cylinder, both extensile and contractile stress can induce nematic order from an isotropic state, unlike fixed-surface active nematics.

desk verdict A well-built asymptotic framework for active nematic films, but the cylindrical isotropic-phase claim rests on an unverified dispersion relation that fails internal consistency checks. read the letter →

arxiv 2506.16523 v1 pith:OSJPTOPA submitted 2025-06-19 cond-mat.soft physics.bio-ph

Mehrana R. Nejad , L. Mahadevan This is my paper · ORCID
classification cond-mat.softphysics.bio-ph
keywords activenematicsthinfilmsasymptoticexpansionlubricationtheorylinearstabilitynematohydrodynamicsfilmthicknesscurvature
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

This paper derives thin-film equations for active nematic layers whose thickness and center-surface shape can change, starting from the full three-dimensional nematohydrodynamic equations with an asymptotic long-wavelength expansion. Using these equations, the authors show that deformability breaks the usual rules of two-dimensional active nematics: in a flat film the nematic phase becomes unstable in all perturbation directions for both extensile and contractile activity, and in a cylindrical film contractile activity drives a thickness instability that deforms the cylinder. The headline result is that in the isotropic phase of a cylindrical film, both extensile and contractile activity can induce nematic order, whereas active nematics on fixed surfaces only order under extensile activity. If correct, this means that thickness changes and curvature are not passive spectators but active players in morphogenetic processes such as gastrulation.

What carries the argument

The central object is the effective film force balance obtained by integrating the three-dimensional Stokes equations across the thin dimension. In the nearly flat case this produces an in-plane tension tensor (Eqs. (40)–(42)) whose entries combine the active stress, proportional to the activity coefficient $m$, the film thickness $h$, and the in-plane nematic components $Q_{ij}$, with viscous shear terms; a thickness evolution equation $\partial_t h = -\nabla\cdot(h\mathbf{u})$ (Eq. (33)); and a torque-balance shape equation $T_{11}\partial_1^2 H + 2T_{12}\partial_1\partial_2 H + T_{22}\partial_2^2 H = 0$ (Eq. (43)). For a curved film the same expansion yields a perpendicular force balance (Eq. (59)) that ties curvature, stretch rates, thickness, and active stress together, along with compatibility relations for the surface metric. These effective equations carry the argument: every new instability follows from the sign of the activity coefficient in the resulting growth rates, and that sign changes for circumferential perturbations because of the curvature terms in the effective equations.

What would settle it

A direct numerical solution of the full three-dimensional equations (1)–(3) for a contractile active nematic cylinder in the isotropic phase ($B=2/3$) at parameters where Eq. (76) predicts growth for circumferential mode $P=2$ and small axial wavenumber would settle the claim: if the order parameter $S$ does not grow and no circumferential thickness or shape bands appear, the sign-reversal result is refuted.

Watch

Extended reading notes

Core claim

On its own terms, the paper establishes that a deformable thin active nematic film supports a richer set of activity-driven instabilities than a fixed two-dimensional layer. In the flat-film nematic phase, the growth rates (Eqs. (44)–(45)) show contractile activity producing a uniform, direction-independent instability and extensile activity destabilizing perturbations along the order direction, so the classical bend-only-for-extensile and splay-only-for-contractile dichotomy disappears when thickness can vary. In the isotropic phase of a cylindrical film, the growth rate for the order parameter $S$ (Eq. (76)) has an activity coefficient whose sign depends on the circumferential mode number $P$: long-axis perturbations ($P=0$) grow only in extensile systems, while circumferential modes $P\ge2$ grow in contractile systems. Because the thickness and center-surface perturbations are slaved to $\delta S$ through Eqs. (77)–(78), any ordering instability immediately produces thickness and shape changes. The authors also show that on a non-deformable cylinder the sign of the activity coefficient never changes (Eq. (80)), isolating deformability as the cause of the reversal.

Load-bearing premise

The predictions assume surface tension is negligible (capillary number $\mathrm{Ca}\gg1$) so the film surfaces are stress-free, and assume the nematic tensor is uniform across the film thickness; if either assumption fails, the instabilities could be suppressed or altered.

Editorial extensions

If this is right

  • In a flat film with variable thickness, the nematic phase is unstable for perturbations in every direction under both extensile and contractile activity; contractile activity gives a leading-order growth rate independent of perturbation angle (Eq. (44)).
  • On a cylinder, contractile activity creates a thickness instability with growth rate $\omega_h = m/(2\mu)$ that deforms the cylinder even when its radius is below the orientational-instability threshold (Eq. (71)).
  • In the isotropic phase of a cylinder, contractile activity can produce nematic order through circumferential modes $P\ge2$, and once order appears the film develops thickness bands and shape changes (Eqs. (76)–(78)).
  • The reversal is caused by deformability: for a non-deformable cylinder the activity coefficient in the isotropic growth rate keeps a single sign, recovering the extensile-only rule (Eq. (80)).
  • For extensile systems, the cylinder's nematic phase still shows the familiar bend instability, but only above a critical radius $R_c = (4\mu K_Q/(-3m\gamma))^{1/2}$; below it the ordered cylinder is stable (Eq. (70)).

Reading between the lines

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

  • By extension, the same sign-reversal mechanism should operate on other curved films, such as spheres or tubes with varying mean curvature, because the paper's general curved-film equations (Table II) contain curvature-weighted active terms; a spherical analogue would predict contractile activity ordering the isotropic phase in high-order angular modes.
  • The contractile thickness instability is a concrete target for experiments with contractile actomyosin gels or microtubule-kinesin suspensions coated on cylindrical substrates: one should observe periodic axial thickenings without prior nematic order, provided surface tension is weak ($\mathrm{Ca}\gg1$).
  • If the flat-film result persists beyond linear order, the standard classification of active nematic defects by bend versus splay activity would need revision in free-standing films; the paper does not analyze defects, but its equations are set up for nonlinear simulations that could check this.
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

2 major / 4 minor

Summary. The manuscript develops asymptotic thin-film models for active nematic layers by reducing the full three-dimensional nematohydrodynamic equations in two regimes: nearly flat films with small curvature, and generally curved films with O(1) curvature. For flat films it derives effective equations for thickness, in-plane velocity, and nematic order and performs linear stability analysis of the nematic and isotropic phases. For cylindrical films it performs a similar analysis and identifies couplings between thickness, shape, and order. The paper's headline claim is that in the isotropic phase on a deformable cylinder, both extensile and contractile activity can induce nematic order, in contrast to fixed-surface active nematics.

Significance. The derivation from the stated 3D equations is mostly systematic, contains no fitted parameters, and the flat-film results recover known 2D active nematic limits in the appropriate reductions, which is a clear strength. The paper also produces explicit, falsifiable dispersion relations and phase diagrams. However, the central claim about contractile ordering in the isotropic cylindrical phase is not supported by the analysis as written: it rests on a quoted dispersion relation, Eq. (76), that has an unphysical long-wavelength divergence and no derivation. The value of the paper is therefore contingent on the authors being able to substantiate or correct Eq. (76).

major comments (2)
  1. [§V.A.2, Eq. (76)] The abstract's claim that contractile activity can induce nematic order in the isotropic phase of a cylindrical film rests on Eq. (76), but the active term in that equation is singular as Q→0 for all circumferential modes P≥2. For fixed P≥2 and contractile activity m>0, the factor (Q²−2P²+2)(Q²+P²)/Q² tends to P²(2−2P²)/Q², so the term −λm(Q²−2P²+2)(Q²+P²)/(9Q²µ) diverges to +∞. This predicts an arbitrarily large growth rate at arbitrarily long axial wavelength for any nonzero contractile activity, with no threshold; such a divergence is unphysical in a Stokes-flow thin-film model. Eq. (76) is quoted without derivation, so the divergence cannot be traced to a justified term in the asymptotic expansion. The authors must either provide the full derivation and show that the singularity is a true physical effect, or correct the formula and re-examine the contractile-ordering claim.
  2. [§V.A.2, Eqs. (76)–(80) and Fig. 4] Eq. (76) is also not consistent with the non-deformable limit that the paper uses as its benchmark. For P=0 and Q→0, Eq. (76) gives an active coefficient −2λm/(9µ), whereas Eq. (80) for a non-deformable cylinder and Eq. (47) for a flat layer at φ=0 both give −λm/(6µ). In addition, the Fig. 4 caption and the surrounding text attribute the contractile-ordering instability to Eq. (80), but Eq. (80) has a sign-definite active coefficient and can only be destabilized by extensile activity; the contractile result must come from Eq. (76). These internal inconsistencies indicate that either Eq. (76) is an erroneous transcription of the dispersion relation or the text describes a different calculation. As written, the central result of the paper is not supported.
minor comments (4)
  1. [§IV, Eq. (34)] The sentence immediately after Eq. (34) says that flows with positive divergence increase the thickness, but the right-hand side of Eq. (34) has a minus sign before h(∂1u+∂2v); positive divergence decreases the thickness. Please correct the wording.
  2. [Fig. 4 caption] The caption states that contractile systems exhibit growth of nematic order in the circumferential direction 'see Eq. (80)', but Eq. (80) is the non-deformable result that has a sign-definite active coefficient and predicts only extensile-driven growth. The citation likely should be to Eq. (76), and the sentence should be reworded after Eq. (76) is corrected.
  3. [End of Sec. III] The list of 13 unknown fields ends with 'a2,a2'; the second entry should be 'a1'.
  4. [SM, Sec. VIII B, Eqs. (S26)-(S27)] The expressions for σ0_12 and σ0_22 in the SM do not match the effective stresses in Eqs. (40)-(42): the shear component should read µ(∂1v+∂2u)+mQ12 and the second normal component should read 2µ(∂1u+2∂2v)+m(Q22−Q33). Please correct these typos or check the torque-balance derivation.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the reduced models and instabilities are derived from stated 3D active-nematic equations with no fitted inputs; self-citations serve only as external benchmarks.

full rationale

The derivation chain starts from the stated 3D incompressible active-nematohydrodynamic equations (Eqs. 1-3) and the nematic free energy (Eq. 4), together with stress-free and kinematic boundary conditions (Eqs. 21-23). The flat-film effective equations (Eqs. 38-43) and the curved-film effective equations (Eqs. 55, 59-61) are obtained by asymptotic expansion in the slenderness parameter, not by imposing the target instability results. The linear stability growth rates, including the cylindrical isotropic-phase result Eq. (76), are computed by Fourier-analyzing these effective equations; no parameter is fitted to the predicted growth rates or to the order-formation criteria. Citations [26,27,35,37] are used only to state the conventional 2D active-nematic benchmarks against which the new film results are compared, and Eqs. (46) and (80) reproduce those benchmarks; they are not used as inputs to the derivation. The assumptions Ca >> 1 and through-thickness homogeneity of Q are stated modeling assumptions rather than hidden redescriptions of the results. The discrepancy between Eq. (76) and the non-deformable limit Eq. (80) noted in the skeptic summary is a potential algebraic or consistency issue, not circularity: if Eq. (76) is erroneous, the central claim would be unsupported, but it would not make the claim equivalent to its own inputs. No fitted-input-as-prediction, self-citation load-bearing, uniqueness-imported-from-authors, ansatz-smuggled-via-citation, or known-result-renaming pattern is present.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new physical entities, fitted parameters, or ad hoc constants. The free parameters are all standard material parameters (A, B, K_Q, mu, m, lambda, gamma); B is set to -2/3 or +2/3 to fix the equilibrium order parameter to S=1 or S=0, which is a normalization, not a fit. The axioms are explicit modeling choices about surface tension, through-thickness order, in-plane director, and the free energy form.

assumptions (5)
  • domain assumption Surface tension is negligible, Ca >> 1, so free surfaces are stress-free (Sec. III B).
    Used to set boundary condition Eq. (21); if violated, a pressure jump across the interface would enter the normal force balance and could stabilize thickness and shape instabilities.
  • domain assumption Active stress is large compared to passive elastic stress, so passive elastic stress is omitted from the momentum balance Eq. (1).
    Stated in Sec. II; the momentum equation includes only viscous, pressure, and active stress, not the Ericksen stress from nematic elasticity.
  • domain assumption The nematic tensor is homogeneous across the thickness, Q = Q(x1,x2,t) (Sec. III B).
    Justified by the claim that free-surface conditions force the order to be parallel to the surface; this removes all n-dependence of Q and restricts the model to thin layers with uniform order through the thickness.
  • domain assumption The nematic orientation lies in the plane of the film, Q13 = Q23 = 0 (Sec. II).
    Assumed for both flat and curved films; out-of-plane tilt of the director is not considered.
  • domain assumption The Landau-de Gennes free energy with a single elastic constant, Eq. (4), describes the passive nematic relaxation.
    Used to write the molecular field Eq. (5); the model does not treat anisotropic elasticity or surface anchoring energies.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Thin active nematohydrodynamic layers: asymptotic theories and instabilities." pith.science (2026). https://pith.science/paper/OSJPTOPA

@misc{pith2026250616523,
  author       = {Pith},
  title        = {Pith review of: Thin active nematohydrodynamic layers: asymptotic theories and instabilities},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OSJPTOPA}},
  note         = {Machine review of arXiv:2506.16523}
}
read the original abstract

Starting from a three-dimensional description of an active nematic layer, we employ an asymptotic theory to derive a series of low-dimensional continuum models that capture the coupled dynamics of flat and curved films, including variations in film thickness, shape deformations, internal velocity fields, and the dynamics of orientational order. Using this asymptotic theory, we investigate instabilities driven by activity in both the nematic and isotropic phases for cylindrical and flat films. In the flat case, we demonstrate that incorporating shape and thickness variations fundamentally alters the bend and splay nature of instabilities compared to conventional two dimensional nematic instabilities. In the isotropic phase, we find that both extensile and contractile activity can induce nematic order, in contrast with active nematics on fixed surfaces, where only extensile activity leads to ordering. For the case of curved geometries such as a cylindrical film, we reveal that thickness and shape instabilities are inherently coupled. In the isotropic phase, the emergence of nematic order triggers both thickness and shape instabilities. In the nematic phase, contractile activity induces thickness instabilities, which in turn drive geometric deformations. Our results highlight the crucial interplay between activity, thickness variations, and curvature, providing new insights into the behavior of active nematic films beyond the conventional two dimensional paradigm that has been studied to date.

Figures

Figures reproduced from arXiv: 2506.16523 by the authors.

Figure 1
Figure 1. FIG. 1. Examples of thickness changes during gastrulation in (a) chick embryo, and (b) chameleon (adapted from Ref. [ [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p012_4.png]

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

40 extracted references · 35 canonical work pages

  1. [1]

    Nematic Phase Using Eqs. (9)–(10), (34), and (38)–(39) we can study the linear stability of a phase characterized by the steady- state solutions us = vs = θs = 0, hs = h0 and Ss = 1, corresponding to a nematic phase aligned along the e1 axis with no flow. For a non-zero flow aligning param- eter, the growthrate of the thickness and the nematic order are c...

    work page

  2. [2]

    Isotropic Phase Using Eqs. (9)–(10), (34), and (38)–(39) we can study the linear stability of a phase characterized by the steady- state solutions us = vs = θs = 0, hs = h0 and Ss = 0, corresponding to an isotropic phase with the director aligned along the x1 axis with no flow. Substituting the fourier decompositionδ ˜f(q,ω ) for the fields (S,θ,h,u,v ), ...

    work page

  3. [3]

    In this subsection, we study the stability of this phase

    Instability of the ordered phase As we explained for the linear stability analysis on the film, dynamics of the nematic tensor has a nematic steady state solution (S = 1) forB =−2/3. In this subsection, we study the stability of this phase. By substituting the Fourier decomposition of the fields ( δf (qx, P,ω )) into equations (9)–(10), (15)–(20), (55), (...

    work page

  4. [4]

    Instability of the isotropic phase To study the stability of an isotropic phase, we use B = 2/3, which leads to a steady state solution ( S = 0). As before, substituting the Fourier decomposition of the fields, the growth rates of the perturbations in the nematic order parameter S, thickness h, and center- surface V3, are respectively given by ωS =−γ−1(3K...

    work page

  5. [5]

    Hatwalne, S

    Y. Hatwalne, S. Ramaswamy, M. Rao, and R. A. Simha, Phys. Rev. Lett. 92, 118101 (2004)

    2004

  6. [6]

    M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha, Rev. Mod. Phys. 85, 1143 (2013)

    2013

  7. [7]

    Chandrakar, M

    P. Chandrakar, M. Varghese, S. Aghvami, A. Baskaran, Z. Dogic, and G. Duclos, Phys. Rev. Lett. 125, 257801 (2020)

    2020

  8. [8]

    Bricard, J.-B

    A. Bricard, J.-B. Caussin, D. Das, C. Savoie, V. Chikkadi, K. Shitara, O. Chepizhko, F. Peruani, D. Saintillan, and D. Bartolo, Nat. Commun. 6, 7470 (2015)

    work page 2015

Show all 40 references
  1. [9]

    Opathalage, M

    A. Opathalage, M. M. Norton, M. P. Juniper, B. Langes- lay, S. A. Aghvami, S. Fraden, and Z. Dogic, PNAS 116, 4788 (2019)

    2019

  2. [10]

    Thijssen, M

    K. Thijssen, M. R. Nejad, and J. M. Yeomans, Phys. Rev. Lett. 125, 218004 (2020)

    2020

  3. [11]

    S. J. DeCamp, G. S. Redner, A. Baskaran, M. F. Hagan, and Z. Dogic, Nat. Mater. 14, 1110 (2015)

    2015

  4. [12]

    H. Xu, M. R. Nejad, J. M. Yeomans, and Y. Wu, PNAS 120, e2219708120 (2023)

    2023

  5. [13]

    H. Li, X. Shi, M. Huang, X. Chen, M. Xiao, C. Liu, H. Chat´ e, and H. Zhang, PNAS 116, 777 (2019)

    2019

  6. [14]

    Y. I. Yaman, E. Demir, R. Vetter, and A. Kocabas, Nat. Commun. 10, 2285 (2019)

    2019

  7. [15]

    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, Nature 544, 212 (2017)

    2017

  8. [16]

    T. B. Saw, W. Xi, B. Ladoux, and C. T. Lim, Adv. Mater. 30, 1802579 (2018)

    2018

  9. [17]

    ˇSkult´ ety, D

    V. ˇSkult´ ety, D. B´ ardfalvy, J. Stenhammar, C. Nardini, and A. Morozov, J. Fluid Mech. 980, A28 (2024)

    2024

  10. [18]

    Maitra, Nat

    A. Maitra, Nat. Phys. 19, 733 (2023)

    2023

  11. [19]

    Mietke, F

    A. Mietke, F. J¨ ulicher, and I. F. Sbalzarini, PNAS 116, 29 (2019)

    2019

  12. [20]

    Teranishi, M

    A. Teranishi, M. Mori, R. Ichiki, S. Toda, G. Shioi, and S. Okuda, bioRxiv , 2023 (2023)

    2023

  13. [21]

    Pearce, New J

    D. Pearce, New J. Phys. 22, 063051 (2020)

    2020

  14. [22]

    Alaimo, C

    F. Alaimo, C. K¨ ohler, and A. Voigt, Sci. Rep. 7, 5211 (2017)

    2017

  15. [23]

    Nitschke, M

    I. Nitschke, M. Nestler, S. Praetorius, H. L¨ owen, and A. Voigt, Proc. R. Soc. A: Math. Phys. Eng. Sci. 474, 20170686 (2018)

    2018

  16. [24]

    Ravichandran, M

    Y. Ravichandran, M. Vogg, K. Kruse, D. J. Pearce, and A. Roux, Sci. Adv. 11, eadr9855 (2025)

    2025

  17. [25]

    Maroudas-Sacks, S

    Y. Maroudas-Sacks, S. Suganthan, L. Garion, Y. Ascoli- Abbina, A. Westfried, N. Dori, I. Pasvinter, M. Popovi´ c, and K. Keren, Devel. 152 (2025)

    2025

  18. [26]

    S. C. Al-Izzi and R. G. Morris, J. Fluid Mech. 957, A4 (2023)

    2023

  19. [27]

    Godr` eche and P

    C. Godr` eche and P. Manneville, Hydrodynamics and nonlinear instabilities (2005)

    2005

  20. [28]

    A. Oron, S. H. Davis, and S. G. Bankoff, Reviews of modern physics 69, 931 (1997)

    1997

  21. [29]

    Howell, European Journal of Applied Mathematics 7, 321 (1996)

    P. Howell, European Journal of Applied Mathematics 7, 321 (1996)

    1996

  22. [30]

    Santhosh, M

    S. Santhosh, M. R. Nejad, A. Doostmohammadi, J. M. Yeomans, and S. P. Thampi, J. Stat. Phys. 180, 699 (2020)

    2020

  23. [31]

    Joshi, S

    C. Joshi, S. Ray, L. M. Lemma, M. Varghese, G. Sharp, Z. Dogic, A. Baskaran, and M. F. Hagan, Phys. Rev. Lett. 129, 258001 (2022)

    2022

  24. [32]

    Chuai, G

    M. Chuai, G. Serrano N´ ajera, M. Serra, L. Mahadevan, and C. J. Weijer, Sci. Adv. 9, eabn5429 (2023)

    2023

  25. [33]

    Van De Fliert, P

    B. Van De Fliert, P. Howell, and J. Ockenden, J. Fluid Mech. 292, 359 (1995)

    1995

  26. [34]

    J. A. Teichman, Wrinkling and sagging of viscous sheets, Ph.D. thesis, Massachusetts Institute of Technology (2002). 17

    2002

  27. [35]

    A. C. Slim, J. Teichman, and L. Mahadevan, Journal of fluid mechanics 694, 5 (2012)

    2012

  28. [36]

    Srinivasan, Z

    S. Srinivasan, Z. Wei, and L. Mahadevan, Physical Re- view Fluids 2, 074103 (2017)

    2017

  29. [37]

    Narayanaswamy, J

    O. Narayanaswamy, J. Am. Ceram. Soc. 60, 1 (1977)

    1977

  30. [38]

    Buckmaster, A

    J. Buckmaster, A. Nachman, and L. Ting, Journal of Fluid Mechanics 69, 1 (1975)

    1975

  31. [39]

    Aditi Simha and S

    R. Aditi Simha and S. Ramaswamy, Phys. Rev. Lett. 89, 058101 (2002)

    2002

  32. [40]

    M. R. Nejad and J. M. Yeomans, PRX Life 1, 023008 (2023)

    2023

Pith tools

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