Hydrodynamic liquid crystal models for lipid bilayers

Axel Voigt; Ingo Nitschke; Jan Magnus Sischka

arxiv: 2512.14374 · v2 · submitted 2025-12-16 · ❄️ cond-mat.soft · math-ph· math.MP· physics.flu-dyn

Hydrodynamic liquid crystal models for lipid bilayers

Ingo Nitschke , Jan Magnus Sischka , Axel Voigt This is my paper

Pith reviewed 2026-05-16 22:04 UTC · model grok-4.3

classification ❄️ cond-mat.soft math-phmath.MPphysics.flu-dyn

keywords lipid bilayershydrodynamic modelsLandau-Helfrich modelBeris-Edwards modelsurface liquid crystalsorder parametermembrane viscosityasymmetric bilayers

0 comments

The pith

Lipid bilayers are described by hydrodynamic surface liquid crystal models with a scalar order parameter for molecular alignment along the surface normal.

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

The paper shows how to derive refined hydrodynamic models for lipid bilayers by starting from surface liquid crystal hydrodynamics and introducing a scalar order parameter that tracks lipid alignment perpendicular to the membrane. This adds molecular detail to earlier approaches that treated the bilayer as a homogeneous fluid while minimizing Helfrich energy. The new models are a surface Landau-Helfrich version for asymmetric bilayers and a surface Beris-Edwards version for symmetric bilayers. In the limit of perfect order both recover the known surface Navier-Stokes-Helfrich equations, giving an alternative derivation of those standard models.

Core claim

Starting from hydrodynamic surface liquid crystal models, we obtain a hydrodynamic surface Landau-Helfrich model for asymmetric lipid bilayers and a surface Beris-Edwards model for symmetric lipid bilayers. The scalar order parameter represents the molecular alignment of the lipids along the surface normal, with other degrees of freedom averaged out. When the alignment is complete, both models reduce exactly to the known surface Navier-Stokes-Helfrich models.

What carries the argument

The scalar order parameter that encodes lipid molecular alignment along the surface normal inside hydrodynamic surface liquid crystal models.

If this is right

Membrane viscosity and lipid alignment dynamics are treated explicitly rather than implicitly.
Asymmetric bilayers obey a Landau-Helfrich-type hydrodynamic system while symmetric bilayers obey a Beris-Edwards-type system.
The models are consistent with existing surface Navier-Stokes-Helfrich descriptions when the order parameter reaches its maximum value.
Continuous descriptions of lipid bilayers now include a controllable degree of molecular ordering.

Where Pith is reading between the lines

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

The same starting liquid-crystal framework could be used to add vector or tensor order parameters for tilted or gel phases.
Numerical implementations would allow direct simulation of how local alignment modulates membrane bending and flow coupling.
The derivation supplies a template for obtaining hydrodynamic models of other surface-confined liquid-crystal systems such as nematic shells.

Load-bearing premise

The lipid bilayer can be represented as a hydrodynamic surface liquid crystal whose molecular alignment is captured by a single scalar order parameter along the surface normal, with all other molecular degrees of freedom averaged out.

What would settle it

A direct measurement of lipid tilt or flow response in an asymmetric bilayer under shear that cannot be reproduced by the evolution equation for the single scalar order parameter.

read the original abstract

Coarse-grained continuous descriptions for lipid bilayers are typically based on minimizing the Helfrich energy. Such models consider the fluid properties of these structures only implicitly and have been shown to nicely reproduce equilibrium properties. Model extensions that also address the dynamics of these structures are surface (Navier--)Stokes--Helfrich models. They explicitly account for membrane viscosity. However, these models also usually treat the lipid bilayer as a homogeneous continuum, neglecting the molecular degrees of freedom of the lipids. Here, we derive refined models which consider in addition a scalar order parameter representing the molecular alignment of the lipids along the surface normal. Starting from hydrodynamic surface liquid crystal models, we obtain a hydrodynamic surface Landau--Helfrich model for asymmetric lipid bilayers and a surface Beris--Edwards model for symmetric lipid bilayers. The fully ordered case for both models leads to the known surface (Navier--)Stokes--Helfrich models. Besides more detailed continuous models for lipid bilayers, we therefore also provide an alternative derivation of surface (Navier--)Stokes--Helfrich models.

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 derives hydrodynamic bilayer models with a scalar order parameter for normal alignment and recovers the standard Stokes-Helfrich limits in the ordered case, but the averaging of other molecular degrees of freedom is the assumption that needs verification.

read the letter

This paper starts from existing hydrodynamic surface liquid crystal models and reduces them by keeping only a scalar order parameter that tracks lipid alignment along the surface normal. For asymmetric bilayers the result is a hydrodynamic surface Landau-Helfrich model; for symmetric bilayers it is a surface Beris-Edwards model. Both recover the known surface Navier-Stokes-Helfrich equations when the order parameter is fully saturated. That reduction serves as a consistency check and also supplies an alternative route to the standard models. The explicit addition of the scalar order parameter is the concrete new step beyond homogeneous continuum treatments. It lets the equations carry information about local ordering while still treating the bilayer as a surface. That could matter for simulations of vesicles or membrane processes where alignment influences bending or flow resistance. The approach is straightforward and stays within established liquid-crystal hydrodynamics, which is a strength. The main soft spot is the modeling decision to average out all other molecular degrees of freedom and retain only the normal scalar. If in-plane director components or tilt fluctuations couple to the hydrodynamic flow at the scales of interest, the resulting stress tensor and dissipation may not be complete. The abstract presents the reduction as direct, but without the intermediate equations it is impossible to check whether the simplifications preserve the correct physics or introduce hidden assumptions. The circularity burden stays low because the new models are built from prior surface liquid-crystal equations rather than fitted to the target limits. This work is for people already working on continuum membrane models in soft-matter physics or biophysics. A reader who wants refined descriptions that include ordering would get value from the new equations, once the derivations are confirmed. It is not a wide-ranging advance but a useful incremental refinement. I would bring it to a reading group to discuss the averaging step. It deserves a serious referee because the grounding in liquid-crystal hydrodynamics is solid and the consistency check with known models is clean, even though the details will need close scrutiny.

Referee Report

2 major / 1 minor

Summary. The manuscript derives refined hydrodynamic models for lipid bilayers starting from established surface liquid crystal hydrodynamics. For asymmetric bilayers it obtains a hydrodynamic surface Landau-Helfrich model that incorporates a scalar order parameter representing molecular alignment along the surface normal; for symmetric bilayers it obtains a surface Beris-Edwards model with the same scalar order parameter. Both models are shown to reduce to the known surface (Navier-)Stokes-Helfrich models in the fully ordered limit, thereby supplying an alternative derivation of those standard models while including additional molecular degrees of freedom.

Significance. If the reductions are free of gaps, the work supplies a systematic route from liquid-crystal hydrodynamics to membrane models that retain a scalar order parameter for lipid alignment. This extends homogeneous continuum treatments by making alignment effects on stress and dissipation explicit, while the ordered-limit consistency check provides a useful validation. The approach could improve dynamical predictions for processes sensitive to local lipid ordering in asymmetric or symmetric bilayers.

major comments (2)

[Derivation sections (reduction from surface LC hydrodynamics)] The central modeling step that replaces the full director field by a single scalar order parameter (averaging out in-plane components and tilt fluctuations) is load-bearing for the stress tensor and dissipation. The manuscript must show explicitly how this averaging is performed and confirm that the resulting hydrodynamic equations remain consistent with the parent liquid-crystal model; without those steps the claim that the reduced models correctly capture alignment-hydrodynamic coupling cannot be verified.
[Fully ordered limit subsection] The statement that the fully ordered limit recovers the known surface Stokes-Helfrich models is presented as a consistency check, yet the intermediate algebra (how the order-parameter contributions to the stress and the Beris-Edwards or Landau-Helfrich terms vanish or integrate) is not shown. This reduction is central to the alternative-derivation claim and must be written out in detail.

minor comments (1)

[Abstract] The abstract introduces the two new models without a one-sentence distinction of their symmetry assumptions or the precise role of the scalar order parameter; a brief clarifying clause would help readers.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation of our manuscript and for the constructive comments on the derivation sections. We have revised the paper to supply the explicit averaging procedure and the full algebra for the ordered limit, thereby strengthening the consistency claims.

read point-by-point responses

Referee: [Derivation sections (reduction from surface LC hydrodynamics)] The central modeling step that replaces the full director field by a single scalar order parameter (averaging out in-plane components and tilt fluctuations) is load-bearing for the stress tensor and dissipation. The manuscript must show explicitly how this averaging is performed and confirm that the resulting hydrodynamic equations remain consistent with the parent liquid-crystal model; without those steps the claim that the reduced models correctly capture alignment-hydrodynamic coupling cannot be verified.

Authors: We agree that the averaging step from the full director field to the scalar order parameter requires explicit derivation. In the revised manuscript we have inserted a dedicated subsection (now Section 3.2) that performs the averaging: we decompose the director into its normal component (the scalar order parameter S) and in-plane/tilt fluctuations, impose statistical isotropy of the fluctuations on the tangent plane, and substitute the resulting moments into the Ericksen stress and dissipation function. The resulting hydrodynamic equations are shown to be consistent with the parent surface LC model by direct substitution and by verifying that the variational structure is preserved. These steps are now written out with intermediate expressions for the averaged stress tensor. revision: yes
Referee: [Fully ordered limit subsection] The statement that the fully ordered limit recovers the known surface Stokes-Helfrich models is presented as a consistency check, yet the intermediate algebra (how the order-parameter contributions to the stress and the Beris-Edwards or Landau-Helfrich terms vanish or integrate) is not shown. This reduction is central to the alternative-derivation claim and must be written out in detail.

Authors: We accept that the ordered-limit reduction must be shown in full. The revised manuscript now contains an expanded subsection (Section 4.3) that carries out the algebra step by step: we set the scalar order parameter to its maximum value, demonstrate that all fluctuation-induced terms in the stress tensor integrate to zero under the surface divergence, and recover the standard surface Navier-Stokes-Helfrich equations together with the Helfrich bending energy. The same procedure is applied to both the Landau-Helfrich and Beris-Edwards cases, with explicit cancellation of the alignment-coupling contributions. revision: yes

Circularity Check

0 steps flagged

Derivation from established surface LC models via averaging approximation, with ordered limit as consistency check

full rationale

The paper starts from hydrodynamic surface liquid crystal models (treated as given inputs) and performs a modeling reduction by retaining only a scalar order parameter for normal alignment while averaging other molecular degrees of freedom. This produces the claimed Landau-Helfrich and Beris-Edwards surface models. The fully ordered limit is shown to recover the known Stokes-Helfrich models, presented explicitly as a consistency verification rather than a new prediction. No equation reduces an output quantity to a fitted parameter or self-citation by construction; the steps are standard coarse-graining approximations whose validity is independent of the target results. The derivation chain is therefore self-contained and non-circular.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of the starting hydrodynamic surface liquid crystal models and on the assumption that a single scalar order parameter suffices to capture lipid alignment; no free parameters or new invented entities are introduced in the abstract.

axioms (2)

domain assumption Lipid bilayers admit a hydrodynamic surface liquid crystal description
Invoked as the starting point for the derivation
domain assumption Molecular alignment of lipids can be represented by a single scalar order parameter along the surface normal
Core modeling choice that distinguishes the new models from homogeneous treatments

pith-pipeline@v0.9.0 · 5499 in / 1303 out tokens · 34340 ms · 2026-05-16T22:04:20.769794+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

[1]

Relaxation dynamics of fluid membranes

[AD09] Marino Arroyo and Antonio DeSimone. “Relaxation dynamics of fluid membranes”. In:Physical Review E79.3 (2009), p. 031915.doi:10.1103/physreve.79.031915. [AI93] A.L. Alexe-Ionescu. “Flexoelectric polarization and second order elasticity for nematic liquid crys- tals”. In:Physics Letters A180.6 (1993), pp. 456–460.doi:10.1016/0375-9601(93)90299-f. [A...

work page doi:10.1103/physreve.79.031915 2009
[2]

Tangential tensor fields on deformable surfaces — how to derive consistent𝐿 2-gradient flows

[NSV23] Ingo Nitschke, Souhayl Sadik, and Axel Voigt. “Tangential tensor fields on deformable surfaces — how to derive consistent𝐿 2-gradient flows”. In:IMA Journal of Applied Mathematics88.6 (2023), pp. 917–958.doi:10.1093/imamat/hxae006. [NV23] Ingo Nitschke and Axel Voigt. “Tensorial time derivatives on moving surfaces: General concepts and a specific ...

work page doi:10.1093/imamat/hxae006 2023

[1] [1]

Relaxation dynamics of fluid membranes

[AD09] Marino Arroyo and Antonio DeSimone. “Relaxation dynamics of fluid membranes”. In:Physical Review E79.3 (2009), p. 031915.doi:10.1103/physreve.79.031915. [AI93] A.L. Alexe-Ionescu. “Flexoelectric polarization and second order elasticity for nematic liquid crys- tals”. In:Physics Letters A180.6 (1993), pp. 456–460.doi:10.1016/0375-9601(93)90299-f. [A...

work page doi:10.1103/physreve.79.031915 2009

[2] [2]

Tangential tensor fields on deformable surfaces — how to derive consistent𝐿 2-gradient flows

[NSV23] Ingo Nitschke, Souhayl Sadik, and Axel Voigt. “Tangential tensor fields on deformable surfaces — how to derive consistent𝐿 2-gradient flows”. In:IMA Journal of Applied Mathematics88.6 (2023), pp. 917–958.doi:10.1093/imamat/hxae006. [NV23] Ingo Nitschke and Axel Voigt. “Tensorial time derivatives on moving surfaces: General concepts and a specific ...

work page doi:10.1093/imamat/hxae006 2023