---
title: Landau–de Gennes Flow Dynamics
url: https://www.emergentmind.com/topics/landau-de-gennes-flow
type: topic
---

# Landau–de Gennes Flow Dynamics

Landau–de Gennes flow is the evolution of a nematic \(Q\)-tensor field driven by a Landau–de Gennes free energy. The order parameter is a symmetric, traceless matrix field,
\[
Q \in S_0^{3\times 3}=\{Q\in \mathbb{R}^{3\times 3}:Q^T=Q,\ \operatorname{tr}Q=0\},
\]
or, in equivalent notation, \(\mathcal S=\{Q\in \mathbb{R}^{3\times 3}:Q^\top=Q,\ \operatorname{tr}Q=0\}\). Uniaxial states have the form \(Q=s(n\otimes n-I/3)\), the isotropic state is \(Q=0\), and biaxial states have three distinct eigenvalues. In the PDE literature, the term usually refers to an \(L^2\)-type gradient flow or a closely related relaxational dynamics for \(Q\); in broader continuum settings it also appears as the no-flow limit of Beris–Edwards models on fixed or evolving surfaces [1505.06143], [2311.06240], [2409.09837].

## 1. Energetic formulation and admissible states

The classical Landau–de Gennes framework combines an elastic energy with a bulk potential. In the one-constant setting,
\[
F[Q]=\int_\Omega \left(\frac{L}{2}|\nabla Q|^2+f_B(Q)\right)\,dx,
\]
with
\[
f_B(Q)=\frac{A}{2}\operatorname{tr}(Q^2)-\frac{B}{3}\operatorname{tr}(Q^3)+\frac{C}{4}\big(\operatorname{tr}(Q^2)\big)^2.
\]
At the isotropic–nematic transition temperature, \(A=B^2/(27C)\), and the bulk potential is minimized both by \(Q=0\) and by the continuum of uniaxial states \(Q=s_+(n\otimes n-I/3)\) with \(s_+=B/(3C)\) [1505.06143].

More general Landau–de Gennes flows arise by replacing the isotropic elastic density with anisotropic quadratic or quartic elastic terms. One anisotropic quadratic elastic energy is
\[
G[Q]=\int_\Omega \big[L_1Q_{ij,k}Q_{ij,k}+L_2Q_{ik,j}Q_{ij,k}+L_3Q_{ij,k}Q_{ik,j}\big]\,dx,
\]
with the structural condition \(L_1>3|L_2+L_3|\), which yields strong Legendre ellipticity and a coercive lower bound. A quartic elastic variant introduces tensorial factors \(S_1(Q)=\frac{s_0}{3}I+Q\), \(S_2(Q)=\frac{2s_0}{3}I-Q\) and the additional term \(\frac{L_5}{2}|Q|^2|\nabla Q|^2\); when \(L_5=0\), it reduces to the Golovaty–Novack–Sternberg quartic elastic model [2011.09541], [2409.09837].

A distinct class of models replaces the quartic bulk polynomial by the Ball–Majumdar singular potential. In that case the physically admissible set is characterized by eigenvalue constraints
\[
\lambda_i(Q)\in\left(-\frac13,\frac23\right),
\]
and the bulk potential is finite only on this set, convex and lower semicontinuous on \(S_0^{3\times 3}\), smooth in the interior, and divergent at the physical boundary. This enforces physical admissibility directly at the continuum level rather than through a posteriori estimates [1207.1643], [1310.8474].

This variety of energies implies that “Landau–de Gennes flow” is not a single canonical PDE. It is a variational class of \(Q\)-tensor evolutions whose analytical and physical properties depend strongly on the elastic sector, the bulk potential, the boundary conditions, and whether hydrodynamic or thermal couplings are retained.

## 2. Relaxational dynamics and dissipation structure

The basic gradient-flow structure is
\[
\partial_t Q=-\frac{\delta F}{\delta Q},
\]
up to mobility, viscosity, or projection factors. For the classical one-constant energy, the \(L^2\) gradient flow is
\[
\gamma \partial_t Q=-\frac{\delta F}{\delta Q},
\]
which becomes
\[
\gamma Q_t=L\Delta Q-AQ+B\left(Q^2-\frac{|Q|^2}{3}I\right)-C|Q|^2Q.
\]
This is a coupled system of five nonlinear parabolic PDEs for the independent components of \(Q\in S_0\), and it satisfies the dissipation identity
\[
\frac{d}{dt}F[Q(t)]=-\frac1\gamma\int_\Omega |Q_t|^2\,dx\le 0
\]
[1505.06143].

For quartic elastic energies, the symmetric-traceless constraint can be enforced explicitly by projection. One formulation is
\[
\partial_t Q=-M\mathcal P\!\left(\frac{\delta \mathcal F(Q)}{\delta Q}\right),
\qquad
\mathcal P(A)=(A+A^\top)-\frac{\operatorname{tr}A}{3}I,
\]
with a variational derivative containing nonlinear \(\operatorname{div}Q\), \(\operatorname{curl}Q\), and \(|Q|^2|\nabla Q|^2\) contributions. The associated continuous energy law is
\[
\frac{d}{dt}\mathcal F(Q)=-M\int_\Omega |\partial_tQ|^2\,dx,
\]
so the dynamics are purely relaxational and converge to local minimizers of the chosen energy [2409.09837].

With anisotropic quadratic elasticity and singular bulk potential, the flow remains gradient-driven but the parabolic operator becomes non-diagonal because of the \(L_2\) and \(L_3\) terms. The full PDE contains
\[
2L_1\Delta Q_{ij}+2(L_2+L_3)\partial_j\partial_kQ_{ik}-(L_2+L_3)\delta_{ij}\partial_k\partial_\ell Q_{\ell k},
\]
together with the projected singular bulk force and the linear destabilizing term \(2aQ\). The convexity of the singular potential is central in the dissipation and comparison theory [2011.09541].

On surfaces, suppressing hydrodynamics and fixing the geometry yields the pure Landau–de Gennes flow limit
\[
\partial_t Q=-\Gamma\,\frac{\delta F_{LdG}}{\delta Q},
\]
or, in conforming surface variables,
\[
\partial_t q=-\Gamma(h_{EL}+h_{TH}),\qquad \partial_t \beta=-\Gamma(\omega_{EL}+\omega_{TH}),
\]
where curvature enters explicitly through the surface mean curvature \(H\), Gaussian curvature \(K\), and shape operator \(I\) [2311.06240].

## 3. Phase fronts, mean-curvature motion, and geometric limits

At the isotropic–nematic transition, the coexistence of the isotropic well \(Q=0\) and the uniaxial nematic well creates a distinguished front-propagation regime. In a three-dimensional droplet with radial hedgehog boundary data, radially symmetric uniaxial dynamics reduce to a scalar equation for the radial order parameter \(h(r,t)\). For “efficient interface” initial data and sufficiently small elastic constant, the isotropic–nematic interface is well defined and, for small times, its radius \(\rho(t)\) satisfies
\[
\frac{d\rho}{dt}=-\frac{2}{\rho},
\]
which is motion by mean curvature for spheres. Numerical experiments also show rapid front formation for broader uniaxial, biaxial, and non-radially symmetric initial data [1505.06143].

A separate sharp-interface scaling uses
\[
E_\varepsilon(Q)=\int_\Omega \frac{\varepsilon}{2}|\nabla Q|^2+\frac{F(Q)}{\varepsilon}\,dx,
\qquad
\partial_tQ^\varepsilon=\Delta Q^\varepsilon-\frac{1}{\varepsilon^2}\partial_QF(Q^\varepsilon).
\]
As \(\varepsilon\to 0\), the energy concentrates near a codimension-one interface \(I_t\). On the isotropic side, \(Q\to 0\); on the nematic side, \(Q\to s_c(u\otimes u-I/3)\), where the director satisfies
\[
\partial_tu=\Delta u+|\nabla u|^2u,
\]
with homogeneous Neumann boundary condition on the evolving interface, while the interface itself satisfies
\[
V_n=H.
\]
The rigorous convergence is proved under well-prepared initial data and as long as the limiting mean curvature flow remains smooth [2010.07786].

On a disc, the long-time behavior depends strongly on planarity and boundary forcing. Planar initial conditions remain planar for all time and preserve an isotropic core, whereas non-planar perturbations can trigger escape into the third dimension and convergence to a smooth globally ordered uniaxial state. Under non-minimal biaxial Dirichlet data, boundary layers form and the interior may become largely ordered or almost entirely isotropic depending on the initial interface location [1505.06143].

Surface formulations introduce an additional geometric layer. In the conforming decomposition \(Q=C_{(q,\beta)}\), the elastic molecular field contains terms such as
\[
\Delta_s q-(H^2-2K)q+3\beta H\Big(I-\frac{H}{2}g^{-1}\Big),
\qquad
\Delta_s\beta+2\,I:q-3\beta(H^2-2K),
\]
so curvature can directly bias alignment and defect structure. This identifies curvature-driven reorientation as an intrinsic part of surface Landau–de Gennes flow rather than an external correction [2311.06240].

## 4. Existence, regularity, and physical admissibility

For the classical quartic-bulk gradient flow on bounded smooth domains with smooth Dirichlet data, standard parabolic theory yields a unique smooth solution for all \(t>0\). In that setting the flow admits an \(L^\infty\) bound
\[
|Q(\cdot,t)|\le \sqrt{\frac23}\,\frac{B}{3C},
\]
which is used repeatedly in the analysis of fronts and long-time dynamics [1505.06143].

With anisotropic elasticity and singular potential, the analytical framework is more delicate. For periodic domains \(\Omega=T^n\), \(n=2,3\), and initial data \(Q_0\in L^2(\Omega;Q_{\mathrm{phy}})\) with \(\psi_{\mathrm{sing}}(Q_0)\in L^1\), the energy is proper, bounded from below, lower semicontinuous, and \(\lambda\)-convex with \(\lambda=-2a\). This yields a unique global EVI gradient-flow solution. For every \(t_0>0\), the solution is strong, belongs to \(L^\infty(t_0,\infty;H^2)\), and satisfies an exact energy identity. Under the stronger coercivity condition
\[
L_1-3|L_2+L_3|-aC_{TR}>0,
\]
the solution eventually detaches from the physical boundary: after some \(T_0\), all eigenvalues satisfy
\[
-\frac13+\kappa\le \lambda_i(Q(t,x))\le \frac23-\kappa.
\]
Before that time, the contact set \(E_t=\{x:Q(t,x)\in \partial Q_{\mathrm{phy}}\}\) satisfies \(\dim_H(E_t)\le 2\) in three dimensions and \(\dim_H(E_t)=0\) in two dimensions [2011.09541].

Non-isothermal extensions couple \(Q\) to incompressible flow and temperature. In one periodic three-dimensional model with singular bulk potential, the unknowns are \((u,Q,\theta)\), the molecular field is
\[
H=[\Delta Q]-[\partial f(Q)]+U(\theta)[\partial G(Q)],
\]
and the total energy and entropy satisfy thermodynamically consistent balances. Global-in-time weak solutions exist for arbitrary physically relevant initial data, and the temperature remains strictly positive:
\[
\theta(t,\cdot)\ge c\,e^{-At}\theta_0
\]
[1207.1643].

A related nonisothermal Beris–Edwards-type model on the torus employs the free energy
\[
\mathcal F(Q,\nabla Q,\vartheta)=\frac12|\nabla_xQ|^2+\vartheta f(Q)-\frac{\lambda}{2}|Q|^2-\vartheta\log\vartheta-a\vartheta^m
\]
and proves global weak solutions under integrated entropy inequalities. A key estimate is a Hessian coercivity inequality for the Ball–Majumdar potential, used to control \(\mathrm{L}[\partial f/\partial Q]\) and \(\Delta Q\) in \(L^2\) [1310.8474].

## 5. Finite element discretization and computational analysis

A recent finite element treatment of quartic elastic Landau–de Gennes flow uses quasi-uniform triangulations and the discrete spaces
\[
\mathcal T_h^0,\qquad \mathcal T_h^g=\mathcal T_h^0+g_h,
\]
consisting of continuous, piecewise linear, symmetric traceless \(3\times 3\) fields. A basis is built from five symmetric-traceless generators \(R_\alpha\) and nodal hat functions \(\psi_i\), so that
\[
\mathcal B_h=\{R_\alpha\psi_i:\alpha\in[5],\, i\in[N]\}
\]
is a basis of \(\mathcal T_h^0\) [2409.09837].

Time discretization is nonlinearly implicit and midpoint-like. Given \(Q_h^n\in \mathcal T_h^g\), one solves for \(Q_h^{n+1}\in \mathcal T_h^g\) through
\[
\frac{Q_h^{n+1}-Q_h^n}{\Delta t}=-M\,H^{n+},
\]
with \(H^{n+}=H_1^{n+}+\cdots+H_6^{n+}\) representing the quartic elastic and bulk contributions evaluated at temporal midpoints. The scheme satisfies the discrete energy dissipation law
\[
\frac{\mathcal F(Q_h^{n+1})-\mathcal F(Q_h^n)}{\Delta t}
=
-\frac1M\left|\frac{Q_h^{n+1}-Q_h^n}{\Delta t}\right|^2,
\]
which gives unconditional energy stability. Solvability at each time step is obtained by a fixed-point iteration starting from \(Q_h^{n+1,0}=Q_h^n\); for sufficiently small \(\Delta t\), the fixed-point map is contractive and the solution is unique [2409.09837].

The same analysis proves coercivity and \(\Gamma\)-convergence of the discrete energies when \(L_5>0\). In weak \(H^1(\Omega;\mathcal S)\), cluster points of discrete global minimizers converge to global minimizers of the continuum energy, and isolated local minimizers are approximated by discrete minimizers. The \(\Gamma\)-convergence argument uses truncation, mollification, and piecewise linear interpolation. For \(L_5=0\), the numerical scheme remains energy stable and exhibits similar dynamics, but the \(\Gamma\)-convergence proof is not established [2409.09837].

The numerical experiments reproduce isotropic-to-nematic phase transitions and tactoid dynamics. On \(\Omega=[0,2]^2\), the reported convergence is approximately second order in \(h\) for \(Q\) and the energy, and second order in \(\Delta t\). A practical CFL-like restriction is observed: the maximum stable \(\Delta t\) scales like \(h^2\). On the unit disk, degree \(+1\), degree \(-1\), and degree \(0\) tactoids show, respectively, splitting into \(\pm 1/2\) or \(-1/2\) defect pairs, or collapse to a uniform nematic state; in all cases the discrete energy decays monotonically [2409.09837].

## 6. Relation to Oseen–Frank theory, hydrodynamics, and other uses of “flow”

A central structural property is the reduction to Oseen–Frank elasticity on the uniaxial manifold. For
\[
Q=s_0\left(n\otimes n-\frac{I}{3}\right),
\]
the quartic elastic model of Golovaty et al. reduces to the Oseen–Frank energy for suitable parameter identifications. The additional quartic term in the modified model is motivated by the identity
\[
\frac{K_1}{2}|\nabla n|^2
=
\frac{3K_1}{8s_0^4}|Q|^2|\nabla Q|^2,
\]
which justifies including \(|Q|^2|\nabla Q|^2\) and links its coefficient directly to \(K_1\) and \(s_0\) [2409.09837].

Landau–de Gennes flow is not synonymous with the Beris–Edwards system. Pure gradient-flow models suppress advection, corotation, velocity-pressure coupling, and Ericksen stresses. By contrast, bulk non-isothermal Beris–Edwards models use
\[
\partial_tQ+\nabla\cdot(Qu)-S(\nabla u,Q)=\Gamma(\theta)H
\]
or closely related forms, and surface Beris–Edwards models use objective derivatives \(D_t^mQ\) or \(D_t^JQ\), inextensible surface flow, and coupled momentum balances. The pure Landau–de Gennes flow is recovered in the no-flow regime, for example by setting \(v=0\), \(v_n=0\), or imposing the NF constraint on surfaces [1207.1643], [2311.06240].

A separate use of the word “flow” appears in functional renormalization-group studies of the Landau–de Gennes model. There the flow variable is the coarse-graining scale \(k\), with \(t=-\ln(k/k_0)\), and the effective average action satisfies the Wetterich equation rather than a time-dependent relaxational PDE. In that framework, the numerical solution of the RG flow indicates a first-order nematic–isotropic transition and yields the estimate
\[
T_c-T^*=5.85\,\mathrm{K}
\]
for the transition-temperature difference in the studied parametrization [1803.03683].

Taken together, these developments establish Landau–de Gennes flow as a broad variational framework for nematic evolution. It accommodates isotropic and nematic phases, biaxiality, defect cores, front propagation, geometric coupling to curvature, singular physical constraints on eigenvalues, and structure-preserving numerical discretization within a common \(Q\)-tensor formalism.

Source: https://www.emergentmind.com/topics/landau-de-gennes-flow