arxiv: 2602.10521 · v2 · submitted 2026-02-11 · ❄️ cond-mat.soft

Defect binding-unbinding transition in active nematic membranes

Yuki Hirota , Nariya Uchida This is my paper

Pith reviewed 2026-05-16 03:50 UTC · model grok-4.3

classification ❄️ cond-mat.soft

keywords active nematicsdeformable membranestopological defectsactive turbulencecurvature couplingdefect bindingnonequilibrium dynamics

0 comments

The pith

Active nematic membranes undergo a continuous transition from trapped defects to active turbulence at a critical activity scaling as coupling strength squared over bending stiffness.

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

The paper establishes that the balance between active stress and anisotropic curvature coupling in nematic liquid crystals on deformable membranes produces a binding-unbinding transition for topological defects. In the low-activity regime, local membrane deformations trap defects and suppress their motion. Above a threshold, activity dominates and drives turbulence while preserving correlations between the director field and membrane shape through wave-like curvature profiles induced by orientational walls. A reader would care because the model supplies a physical mechanism for defect-driven deformations in nonequilibrium biological membranes using only a minimal set of equations.

Core claim

Simulations of the coupled evolution of the nematic order parameter and membrane height field reveal a continuous transition from a curvature-dominated regime, in which topological defects remain trapped by local deformations, to an activity-dominated regime that exhibits active turbulence. Scaling analysis shows the critical activity threshold ζ_c scales as α²/κ, where α is the coupling constant and κ is the bending stiffness; this relation is confirmed numerically. Even in the turbulent regime, walls in the director field generate characteristic wave-like curvature profiles that maintain significant orientational-geometry correlations.

What carries the argument

The minimal coupled dynamical equations for the nematic order parameter and membrane height, incorporating active stress and anisotropic curvature coupling.

Load-bearing premise

The minimal model of coupled nematic order parameter and membrane height evolution captures the dominant physics of active stress and anisotropic curvature coupling without requiring additional biological or hydrodynamic details.

What would settle it

Numerical or experimental data showing either a discontinuous jump in defect mobility with activity or a critical threshold that deviates from the α²/κ scaling would falsify the continuous transition and the reported scaling relation.

Figures

Figures reproduced from arXiv: 2602.10521 by Nariya Uchida, Yuki Hirota.

**Figure 2.** Figure 2: FIG. 2. Spatial correlation functions and time evolution of correlation lengths. (a) Orientational correlation function [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Phase diagram showing the time-averaged bending [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. Activity dependence of the bending energy and fluid [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗

**Figure 5.** Figure 5: FIG. 5. Snapshots of the director field and the mean curvature field for a fixed [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗

read the original abstract

We investigate the dynamics of active nematic liquid crystals on deformable membranes, focusing on the interplay between active stress and anisotropic curvature coupling. Using a minimal model, we simulate the coupled evolution of the nematic order parameter and membrane height. We demonstrate a continuous transition from a curvature-dominated regime, where topological defects are trapped by local deformation, to an activity-dominated regime exhibiting active turbulence. A scaling analysis reveals that the critical activity threshold $\zeta_c$ scales as $\alpha^2/\kappa$, where $\alpha$ and $\kappa$ are the coupling constant and bending stiffness, respectively; this relationship is confirmed by our numerical results. Furthermore, we find that significant correlations between the orientational pattern and membrane geometry persist even in the turbulent regime. Specifically, we identify that "walls" in the director field induce characteristic wave-like curvature profiles, providing a mechanism for dynamic coupling between order and shape. These results offer a physical framework for understanding defect-mediated deformation in nonequilibrium biological membranes.

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.

Desk Editor's Note private letter to a colleague

The paper identifies a continuous defect unbinding transition in active nematic membranes at activity scaling as α²/κ, with numerics confirming the scaling and persistent curvature waves from director walls.

read the letter

Colleague, the main thing here is a continuous transition in active nematics on deformable membranes: at low activity, curvature traps the defects, but above a threshold they unbind and the system goes turbulent. The scaling they extract is ζ_c ~ α²/κ from balancing active stress against the curvature coupling and bending rigidity, and the simulations reproduce it. They also note that director walls still drive wave-like curvature profiles even once turbulence sets in, so the order-shape link does not vanish completely. That combination of scaling plus the wave observation is the concrete addition beyond flat-substrate active nematic work. The model is deliberately minimal and drops hydrodynamics, which keeps the equations tractable and lets the balancing argument stay clean. Within those limits the numerics look direct: they track defect motion and flow to mark the transition without obvious fitting or circularity. The weakest part is that real membranes involve fluid flow and extra biology, so the exact location of ζ_c could shift, but the paper treats the result as a minimal-framework prediction rather than a universal one. No internal contradictions appear in the scaling or the reported checks. This is aimed at people working on active matter on curved or soft substrates, especially those thinking about biological membranes. The scaling and the mechanism are specific enough to be useful for guiding simulations or experiments in that niche. I would bring it to a reading group and would cite the scaling relation if I needed a concrete nonequilibrium membrane example. It deserves peer review; the analysis and numerical support are solid enough inside the stated model that referees can assess the details and the limitations without it being a desk reject.

Referee Report

2 major / 2 minor

Summary. The manuscript investigates the dynamics of active nematic liquid crystals on deformable membranes using a minimal coupled model for the nematic order parameter and membrane height. It reports a continuous transition from a curvature-dominated regime, in which topological defects are trapped by local deformations, to an activity-dominated regime exhibiting active turbulence. A scaling analysis predicts that the critical activity threshold ζ_c scales as α²/κ, with this relation confirmed by numerical simulations. Persistent correlations between the director field and membrane geometry are identified even in the turbulent regime, specifically wave-like curvature profiles induced by walls in the orientational pattern.

Significance. If the central scaling result holds, the work supplies a clean, minimal hydrodynamic-free framework for defect-mediated shape changes in active membranes, with direct relevance to biological systems. The dimensionally consistent balancing argument for ζ_c ~ α²/κ together with its direct numerical verification constitutes a falsifiable prediction that strengthens the claim; the additional observation of geometry-order correlations in the turbulent phase adds mechanistic insight beyond the transition itself.

major comments (2)

[§3] §3 (scaling analysis): the derivation of ζ_c ~ α²/κ balances active stress against the anisotropic curvature coupling and bending terms, but the precise prefactors and the assumed form of the curvature-coupling stress tensor are not written out explicitly; without these the scaling cannot be checked for hidden parameter dependence.
[§4.3] §4.3 (numerical results): the location of the transition is identified by the onset of defect unbinding and the appearance of turbulent flow, yet no quantitative criterion (e.g., defect-pair separation threshold or order-parameter correlation length) is stated; this makes the reported numerical confirmation of the scaling relation difficult to reproduce exactly.

minor comments (2)

[Figure 3] Figure 3: the color scale for membrane height h is not labeled with units or a reference value; adding this would clarify the amplitude of the reported wave-like profiles.
[§2] The model equations in §2 omit explicit mention of the Frank elastic constants; stating whether they are set to unity or retained as parameters would remove ambiguity in the scaling analysis.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their positive assessment and constructive comments on our manuscript. We address each major comment below and will revise the manuscript accordingly to improve clarity and reproducibility.

read point-by-point responses

Referee: [§3] §3 (scaling analysis): the derivation of ζ_c ~ α²/κ balances active stress against the anisotropic curvature coupling and bending terms, but the precise prefactors and the assumed form of the curvature-coupling stress tensor are not written out explicitly; without these the scaling cannot be checked for hidden parameter dependence.

Authors: We agree that the scaling derivation in §3 would benefit from explicit details. In the revised manuscript we will write out the full form of the curvature-coupling stress tensor and show the term-by-term balancing that yields ζ_c ~ α²/κ, including all prefactors, so that any hidden parameter dependence can be verified directly. revision: yes
Referee: [§4.3] §4.3 (numerical results): the location of the transition is identified by the onset of defect unbinding and the appearance of turbulent flow, yet no quantitative criterion (e.g., defect-pair separation threshold or order-parameter correlation length) is stated; this makes the reported numerical confirmation of the scaling relation difficult to reproduce exactly.

Authors: We acknowledge the need for a reproducible quantitative criterion. In the revised §4.3 we will define the transition point explicitly, for example as the activity at which the mean defect-pair separation exceeds a fixed fraction of the system size or at which the orientational correlation length falls below a stated threshold; this definition will then be used to locate the numerical points that confirm the ζ_c ~ α²/κ scaling. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper performs a scaling analysis on its minimal coupled equations for nematic order and membrane height to obtain the relation ζ_c ~ α²/κ by balancing active stress, anisotropic curvature coupling, and bending rigidity; this analytic result is then checked by direct numerical simulation of defect unbinding. No load-bearing step reduces to a self-definition, a fitted parameter renamed as a prediction, or a self-citation chain. The derivation remains independent of its own outputs and is externally falsifiable within the stated model assumptions.

Axiom & Free-Parameter Ledger

3 free parameters · 1 axioms · 0 invented entities

The central claim rests on a standard continuum description of active nematics plus membrane elasticity; no new entities are introduced and the free parameters are the usual activity, coupling, and stiffness coefficients.

free parameters (3)

activity strength zeta
Control parameter swept in simulations to locate the transition threshold.
coupling constant alpha
Strength of anisotropic curvature coupling between nematic order and membrane shape.
bending stiffness kappa
Membrane rigidity appearing in the scaling relation.

axioms (1)

domain assumption Continuum hydrodynamic description of nematic order parameter and membrane height remains valid across the transition
Invoked by the choice of minimal model for the coupled evolution equations.

pith-pipeline@v0.9.0 · 5461 in / 1411 out tokens · 80895 ms · 2026-05-16T03:50:50.987848+00:00 · methodology

discussion (0)

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.

A scaling analysis reveals that the critical activity threshold ζ_c scales as α²/κ ... F_curv = -α²/(2κ) ∫ |J(q)|²/q⁴
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

minimal model of coupled nematic order parameter and membrane height evolution

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

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

[1]

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)

work page 2013
[2]

Ramaswamy, Annu

S. Ramaswamy, Annu. Rev. Condens. Matter Phys.1, 323 (2010)

work page 2010
[3]

Bechinger, R

C. Bechinger, R. Di Leonardo, H. L¨ owen, C. Reichhardt, G. Volpe, and G. Volpe, Rev. Mod. Phys.88, 045006 (2016)

work page 2016
[4]

M. J. Bowick, N. Fakhri, M. C. Marchetti, and S. Ra- maswamy, Phys. Rev. X12, 010501 (2022)

work page 2022
[5]

Doostmohammadi, J

A. Doostmohammadi, J. Ign´ es-Mullol, J. M. Yeomans, and F. Sagu´ es, Nat. Commun.9, 3246 (2018)

work page 2018
[6]

Balasubramaniam, R.-M

L. Balasubramaniam, R.-M. M` ege, and B. Ladoux, Curr. Opin. Genet. Dev.73, 101897 (2022)

work page 2022
[7]

Prost, F

J. Prost, F. J¨ ulicher, and J.-F. Joanny, Nat. Phys.11, 111 (2015)

work page 2015
[8]

Sanchez, D

T. Sanchez, D. T. Chen, S. J. DeCamp, M. Heymann, and Z. Dogic, Nature491, 431 (2012)

work page 2012
[9]

H. H. Wensink, J. Dunkel, S. Heidenreich, K. Drescher, R. E. Goldstein, H. L¨ owen, and J. M. Yeomans, Proc. Natl. Acad. Sci. U.S.A.109, 14308 (2012)

work page 2012
[10]

L. M. Lemma, S. J. DeCamp, Z. You, L. Giomi, and Z. Dogic, Soft Matter15, 3264 (2019)

work page 2019
[11]

Giomi, Phys

L. Giomi, Phys. Rev. X5, 031003 (2015)

work page 2015
[12]

F. C. Keber, E. Loiseau, T. Sanchez, S. J. DeCamp, L. Giomi, M. J. Bowick, M. C. Marchetti, Z. Dogic, and A. R. Bausch, Science345, 1135 (2014)

work page 2014
[13]

P. W. Ellis, D. J. Pearce, Y.-W. Chang, G. Goldsztein, L. Giomi, and A. Fernandez-Nieves, Nat. Phys.14, 85 (2018)

work page 2018
[14]

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

work page 2017
[15]

Kawaguchi, R

K. Kawaguchi, R. Kageyama, and M. Sano, Nature545, 327 (2017)

work page 2017
[16]

Maroudas-Sacks, L

Y. Maroudas-Sacks, L. Garion, L. Shani-Zerbib, A. Livshits, E. Braun, and K. Keren, Nat. Phys.17, 251 (2021)

work page 2021
[17]

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, Development152, DEV204514 (2025)

work page 2025
[18]

X. Wang, M. Merkel, L. B. Sutter, G. Erdemci-Tandogan, M. L. Manning, and K. E. Kasza, Proc. Natl. Acad. Sci. U.S.A.117, 13541 (2020). 9

work page 2020
[19]

N. I. Petridou, S. Grigolon, G. Salbreux, E. Hannezo, and C.-P. Heisenberg, Nat. Cell Biol.21, 169 (2019)

work page 2019
[20]

T. Chen, Y. Zhao, X. Zhao, S. Li, J. Cao, J. Guo, W. Bu, H. Zhao, J. s. Du, Y. Cao,et al., Adv. Sci.9, 2104301 (2022)

work page 2022
[21]

Etournay, M

R. Etournay, M. Merkel, M. Popovi´ c, H. Brandl, N. A. Dye, B. Aigouy, G. Salbreux, S. Eaton, and F. J¨ ulicher, eLife5, e14334 (2016)

work page 2016
[22]

Erdemci-Tandogan and M

G. Erdemci-Tandogan and M. L. Manning, PLoS Com- put. Biol.17, e1009049 (2021)

work page 2021
[23]

Mietke, F

A. Mietke, F. J¨ ulicher, and I. F. Sbalzarini, Proc. Natl. Acad. Sci. U.S.A.116, 29 (2019)

work page 2019
[24]

B¨ acher and S

C. B¨ acher and S. Gekle, Phys. Rev. E99, 062418 (2019)

work page 2019
[25]

Vafa and L

F. Vafa and L. Mahadevan, Phys. Rev. Lett.129, 098102 (2022)

work page 2022
[26]

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

work page 2023
[27]

S. C. Al-Izzi and R. G. Morris, inSemin. Cell Dev. Biol., Vol. 120 (Elsevier, 2021) pp. 44–52, see this reference for a topical review of active, deformable surfaces in mor- phogenesis

work page 2021
[28]

Salbreux, F

G. Salbreux, F. J¨ ulicher, J. Prost, and A. Callan-Jones, Phys. Rev. Research4, 033158 (2022), this reference pro- vides a comprehensive theoretical treatment of active fluid surfaces

work page 2022
[29]

L. A. Hoffmann, L. N. Carenza, J. Eckert, and L. Giomi, Sci. Adv.8, eabk2712 (2022)

work page 2022
[30]

Z. Wang, M. C. Marchetti, and F. Brauns, Proc. Natl. Acad. Sci. U.S.A.120, e2220167120 (2023)

work page 2023
[31]

Nitschke and A

I. Nitschke and A. Voigt, Proc. R. Soc. A481(2025)

work page 2025
[32]

Mirza, A

W. Mirza, A. Torres-S´ anchez, G. Vilanova, and M. Ar- royo, New J. Phys.27, 043025 (2025)

work page 2025
[33]

C. Zhu, D. Saintillan, and A. Chern, Proc. R. Soc. A481 (2025)

work page 2025
[34]

Metselaar, J

L. Metselaar, J. M. Yeomans, and A. Doostmohammadi, Phys. Rev. Lett.123, 208001 (2019)

work page 2019
[35]

Khoromskaia and G

D. Khoromskaia and G. Salbreux, eLife12, e75878 (2023)

work page 2023
[36]

Uchida, Phys

N. Uchida, Phys. Rev. E66, 040902 (2002)

work page 2002
[37]

S. P. Thampi, R. Golestanian, and J. M. Yeomans, EPL 105, 18001 (2014)

work page 2014
[38]

Napoli and L

G. Napoli and L. Vergori, Phys. Rev. E85, 061701 (2012), see this reference for general treatments of spatial derivatives on curved nematic surfaces

work page 2012
[39]

Pollard, S

J. Pollard, S. C. Al-Izzi, and R. G. Morris, J. Fluid Mech. 1007, A32 (2025), see this reference for the choice of objective time derivatives and coordinate-invariance of the physics on fluid deformable surfaces

work page 2025
[40]

Arroyo and A

M. Arroyo and A. DeSimone, Phys. Rev. E79, 031915 (2009), see this reference for an exploration of couplings between flow and intrinsic curvature

work page 2009
[41]

M. L. Henle and A. J. Levine, Phys. Rev. E81, 011905 (2010), see this reference for the hydrodynamics of fluid membranes with curvature-flow coupling

work page 2010
[42]

R. G. Morris and M. S. Turner, Phys. Rev. Lett.115, 198101 (2015), see this reference for the influence of in- trinsic curvature on the mobility of membrane inclusions

work page 2015
[43]

Kinoshita and N

Y. Kinoshita and N. Uchida, Phys. Rev. E108, 014605 (2023)

work page 2023
[44]

V´ elez-Cer´ on, P

I. V´ elez-Cer´ on, P. Guillamat, F. Sagu´ es, and J. Ign´ es- Mullol, Proc. Natl. Acad. Sci. U.S.A.121, e2312494121 (2024)

work page 2024
[45]

Spontaneous rotation and propulsion of suspended capsules in active nematics

J. Santos, M. M. T. da Gama, and R. C. Coelho, arXiv:2510.17643 (2025)

work page internal anchor Pith review Pith/arXiv arXiv 2025
[46]

Shankar and M

S. Shankar and M. C. Marchetti, Phys. Rev. X9, 041047 (2019)

work page 2019
[47]

E. J. Hemingway, P. Mishra, M. C. Marchetti, and S. M. Fielding, Soft Matter12, 7943 (2016)

work page 2016
[48]

Hirota and N

Y. Hirota and N. Uchida, Data and code for: Defect binding-unbinding transition in active nematic mem- branes, https://doi.org/10.5281/zenodo.19394561

work page doi:10.5281/zenodo.19394561