Hydrodynamic Surface Landau–Helfrich Model

• The hydrodynamic surface Landau–Helfrich model is a continuum framework that integrates membrane hydrodynamics, curvature elasticity, and lipid order to capture biomembrane dynamics.
• It employs techniques such as enforcing local area preservation and Allen–Cahn evolution to analyze vesicle morphogenesis and dynamic ordering transitions.
• The model underpins investigations of asymmetric bilayers using scalar or tensorial order parameters, bridging classical Helfrich mechanics with liquid-crystalline descriptions.

The hydrodynamic surface Landau–Helfrich model is a continuum theoretical framework describing the coupled dynamics of lipid bilayers, incorporating membrane hydrodynamics, curvature elasticity, and internal liquid crystalline order. It generalizes the traditional Helfrich energy minimization by dynamically coupling a surface velocity field, membrane bending, and a scalar or tensorial order parameter representing local lipid molecular orientation. Primary applications include the study of biomembrane mechanics, vesicle morphogenesis, and self-organization phenomena in soft and biological matter (Nitschke et al., 16 Dec 2025).

1. Physical Setup and Modeling Assumptions

The model is formulated for a smooth, closed, evolving surface $\Gamma(t)\subset\mathbb{R}^3$ representing the mid-surface of a lipid bilayer. At each $x\in\Gamma$, the unit surface normal is $n(x)$, and the orthogonal projection operator onto the tangent plane is $P=I-n\otimes n$. The material velocity is $V:\Gamma\to\mathbb{R}^3$, with normal ($V\cdot n = v_n$) and tangential ($v = P V$) components. Inextensibility (local area conservation) is strictly enforced by $\operatorname{div}_\Gamma v = v_n H$, where $H$ is the surface mean curvature.

Beyond the classical homogeneous continuum assumption, a scalar order parameter $\beta:\Gamma\to[0,2/3]$ quantifies the local degree of lipid alignment along $n$, with $\beta=2/3$ indicating full order (perpendicular alignment) and $\beta=0$ isotropy. This scalar (or, in more general models, a Q-tensor) encodes internal molecular order absent in standard Helfrich models, enabling refined treatment of asymmetric and liquid-crystalline bilayer structures (Nitschke et al., 16 Dec 2025, Nitschke et al., 2023).

2. Landau–Helfrich Free Energy Functional

The total free energy functionally couples membrane shape and order via: $$F[\beta] = \int_\Gamma \left( \frac{3L}{4} |\nabla_\Gamma \beta|^2 + \frac{9}{4}\beta^2\left[\frac{\kappa}{2}\left(H-(\hat{H}_0+\frac{3}{2}(H_0-\hat{H}_0)\beta)\right)^2 - \bar{\kappa} K\right] + f_{\text{dw}}(\beta) \right) dA$$ where $L$ is the order-parameter elastic constant, $H$ and $K$ are the mean and Gaussian curvatures, $\kappa$ and $\bar{\kappa}$ the bending and Gaussian-modulus coefficients, and $H_0$, $\hat{H}_0$ the spontaneous curvatures in fully ordered and isotropic states. $f_{\text{dw}}(\beta)$ is a double-well potential with minima at $\beta=0$ and $\beta=2/3$, defined as

$$f_{\text{dw}}(\beta) = \frac{9}{4}\beta^2\left(\frac{\hat{a}}{9}(2/3-\beta)^2 - \frac{3\varpi}{2}\beta(4 - \frac{9}{2}\beta)\right)$$

with $\varpi,\hat{a}>0$ for boundedness (Nitschke et al., 16 Dec 2025).

This energy penalizes spatial inhomogeneity in $\beta$ (gradient energy), energetically favors ordered/disordered domains (double well), and couples order to bending rigidity and spontaneous curvature, with a term $H_0(\beta)$ interpolating between $H_0$ and $\hat{H}_0$ for asymmetric bilayers.

3. Dynamic Equations and Constitutive Laws

3.1 Momentum Balance

The surface Navier–Stokes equation with additional order and curvature-coupling stresses reads: $$\rho[D_t V] = -\nabla_\Gamma p + F_{\text{viscous}} + F_{\text{IM}} + F_{\text{LH}}$$ with $\rho$ the surface mass density and $D_t$ the material derivative on $\Gamma$. The tangential and normal force components are:

• Tangential: anisotropic surface viscous stress, order-elastic stress, curvature–order and material-immobility forces.
• Normal: pressure, normal viscous stress, normal material-immobility, classical Helfrich force, and higher-order curvature–order couplings.

3.2 Inextensibility Constraint

Local area preservation is enforced: $\operatorname{div}_\Gamma v = v_n H$

3.3 Order Parameter Evolution

The Allen–Cahn-type evolution for $\beta$ is: $$(M+\frac{\xi^2}{2})D_t \beta = L\Delta_\Gamma \beta - \frac{3}{2}\kappa\beta(...)+3\bar{\kappa} K\beta - ...$$ where $M$ is the order mobility, $\xi$ an anisotropy parameter, and $\Delta_\Gamma$ the Laplace–Beltrami operator. The detailed nonlinear term couples $\beta$ both to geometry and the double-well structure, enabling capture of dynamic ordering transitions driven by curvature and flow (Nitschke et al., 16 Dec 2025).

3.4 Constitutive Relations

Constitutive Effect	Mathematical Form	Physical Role
Anisotropic viscous stress	$$\sigma_{\text{viscous}} = (1+\frac{\xi}{2}\beta)^2[2D_\Gamma(v)-(\text{div}_\Gamma v)I]$$	Direction-dependent membrane viscosity
Order-elastic stress	$$\propto L\nabla_\Gamma\beta\otimes\nabla_\Gamma\beta$$	Opposes gradients in molecular order
Curvature–order coupling	$$\propto \kappa\beta^2(H-\dots)P + \bar{\kappa}\beta^2 P$$	Bending modulus and spontaneous curvature coupled to $\beta$
Double-well “pressure”	$\propto f_{\text{dw}}(\beta)I$	Isotropic domain-forming pressure

(Nitschke et al., 16 Dec 2025)

4. Special Cases and Relationship to Other Models

In the limit $\beta\equiv 2/3$ (“fully ordered”), the evolution equation collapses to a Lagrange multiplier enforcing $\beta=2/3$, all immobility terms vanish with $M\to 0$, and the model reduces to the classical surface (Navier–)Stokes–Helfrich system: $$\rho D_t V = -\nabla_\Gamma p + \mu \Delta_\Gamma v - \kappa \nabla_\Gamma \cdot[(H-H_0)P] + ...$$ with inextensibility $\operatorname{div}_\Gamma v = v_n H$.

Tensorial generalizations employ a surface Landau–de Gennes (LdG) functional for a symmetric traceless Q-tensor order parameter, resulting in a coupled Navier–Stokes–LdG–Helfrich system. Beris–Edwards-type models on surfaces extend this framework to model nematic liquid crystals with full Q-tensor dynamics and surface anchoring constraints, covering both symmetric (Q-tensor) and scalar (uniaxial) cases (Nitschke et al., 2023).

If order is replaced by a phase field $c:\Gamma\to[-1,1]$, the model connects to Navier–Stokes–Cahn–Hilliard-type surface systems capturing two-phase coexistence and membrane domain formation (Bachini et al., 2023).

5. Boundary and Initial Conditions

For closed vesicles, the default is no-flux for $\beta$ and stress-free or prescribed membrane traction, i.e., $\partial_\Gamma = \emptyset$. For $\partial_\Gamma\neq\emptyset$, boundary conditions specify either normal velocities or tangential stress, with Neumann-type conditions for $\beta$ ($\partial_n\beta = 0$).

Anchoring and other constraints may be imposed via Lagrange multipliers, including tangential anchoring for the Q-tensor or conservation of particular eigenvalues (Nitschke et al., 2023, Nitschke et al., 16 Dec 2025).

6. Dimensionless Parameters and Scaling Laws

Key nondimensional groups governing the model dynamics include:

Symbol	Name	Expression	Interpretation
$\mathrm{Re}$	Reynolds number	$\rho U R / \mu_\mathrm{eff}$, $\mu_\mathrm{eff}=(1+\xi/3)^2\mu$	Inertia/viscosity
$\mathrm{Ca}_b$	Bending capillary number	$\mu_\mathrm{eff} U R^2/\kappa$	Bending/viscous dissipation
$\mathrm{Pe}$	Order Péclet number	$MU R / L$	Advective vs diffusive transport of $\beta$
$C_o$	Curvature–order coupling	$H_0 R$	Spontaneous curvature effect
$A$	Anisotropy number	$\xi\beta$	Hydrodynamic anisotropy

These parameters regulate the interplay of hydrodynamics, elasticity, and phase/ordering dynamics, controlling the dominant relaxation and driven mechanisms (Nitschke et al., 16 Dec 2025).

7. Generalizations and Connections

The hydrodynamic surface Landau–Helfrich model forms a basis for refined continuous descriptions of biomembranes, linking to established models such as:

• Classical surface (Navier–)Stokes–Helfrich for fully ordered (homogeneous) regimes.
• Beris–Edwards Q-tensor models for nematic/symmetric bilayers on evolving surfaces, incorporating constraints such as tangential anchoring, constant normal eigenvalue, or uniaxial/biaxial conditions (Nitschke et al., 2023).
• Navier–Stokes–Cahn–Hilliard-type models for phase-separating, two-phase (raft-like) membranes (Bachini et al., 2023).

A plausible implication is that these frameworks together enable rigorous and thermodynamically consistent modeling of a spectrum of membrane phenomena, ranging from reversible shape fluctuations and order-driven curvature sorting to dynamic phase separation and flow-coupled reorganization in biological interfaces.