---
title: Landau–de Gennes Model
url: https://www.emergentmind.com/topics/landau-de-gennes-model
type: topic
---

# Landau–de Gennes Model

Searching arXiv for recent and foundational Landau–de Gennes papers relevant to the requested encyclopedia article.
The Landau–de Gennes model is a continuum theory of liquid crystals in which the local state is represented by a symmetric, traceless tensor \(Q(x)\), and the free energy combines elastic penalties for spatial distortion with a bulk potential that selects orientationally ordered phases below a temperature threshold. In the cited literature, this framework appears in its classical nematic form, in asymptotic reductions to Oseen–Frank and thin-film models, and in modified versions that incorporate chirality, smectic layering, free boundaries, anisotropic elasticity, and homogenized microstructure effects [2601.12653] [1501.07339] [1501.05236].

## 1. Order parameter, bulk potential, and vacuum structure

The basic Landau–de Gennes state space is
\[
S_0:=\{Q\in \mathbb{R}^{3\times 3}: Q^T=Q,\ \operatorname{tr}Q=0\}.
\]
The isotropic state is \(Q=0\). A uniaxial state has the form
\[
Q=s\left(n\otimes n-\frac{I_3}{3}\right),
\]
where \(n\in S^2\) is the director and \(s\) is the scalar degree of order; biaxial states have three distinct eigenvalues [2601.12653] [2410.09930].

In the standard quartic theory, the bulk potential is
\[
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,
\]
with \(A=\alpha_1(T-T_1^*)\) and \(B,C>0\). In the three-dimensional setting summarized for the modified smectic model, the isotropic state minimizes \(f_b\) if \(A>\frac{B^2}{27C}\), whereas for \(A<\frac{B^2}{27C}\) the minimizers are uniaxial with
\[
s_+=\frac{B+\sqrt{B^2-24AC}}{4C}.
\]
This is the canonical thermotropic mechanism by which the bulk term selects the nematic well [2408.03343].

Several cited works use equivalent normalizations. In the vanishing-elasticity literature, the bulk density is written as
\[
f(Q)=k-\frac{a}{2}\operatorname{tr}(Q^2)-\frac{b}{3}\operatorname{tr}(Q^3)+\frac{c}{4}\big(\operatorname{tr}(Q^2)\big)^2,
\]
with \(k\) chosen so that \(\inf f=0\), and the vacuum manifold becomes
\[
\mathcal N=\left\{s_*\left(n\otimes n-\frac{I}{3}\right):n\in S^2\right\},
\qquad
s_*=\frac{b+\sqrt{b^2+24ac}}{4c}.
\]
Because \(n\equiv -n\), this vacuum manifold is diffeomorphic to \(\mathbb{RP}^2\), which is the topological source of half-integer disclinations in the classical nematic theory [1501.05236].

The bulk structure can also be generalized beyond the quartic case. A sextic potential analyzed in three dimensions takes the form
\[
f_b(Q)=a_1-a_2\operatorname{tr}(Q^2)+a_4\big(\operatorname{tr}(Q^2)\big)^2+a_6\big(\operatorname{tr}(Q^2)\big)^3+a_6' \big(\operatorname{tr}(Q^3)\big)^2,
\]
with \(a_2,a_4,a_6,a_6'>0\) and \(a_1\) chosen so that \(\min f_b=0\). Its vacuum manifold is no longer uniaxial:
\[
N=\{r_*(n\otimes n-m\otimes m):(n,m)\in M,\ n\cdot m=0\},
\]
where \(r_*>0\) solves \(4a_6r_*^4+2a_4r_*^2-a_2=0\). In that model, the well manifold is a compact smooth three-dimensional submanifold of \(S_0\) diffeomorphic to \(S^3/Q_8\) [2404.00677].

## 2. Elastic sector and variational formulation

The simplest Landau–de Gennes energy combines the bulk potential with a one-constant elastic term,
\[
E_\varepsilon(Q)=\int_\Omega \left\{\frac12 |\nabla Q|^2+\frac{1}{\varepsilon^2}f(Q)\right\}\,dx,
\]
or, equivalently, \(\mathcal F_L[Q]=\int_\Omega \{\frac{L}{2}|\nabla Q|^2+f_B(Q)\}\,dx\) with \(L=\varepsilon^2\) [1501.05236]. This form underlies much of the asymptotic defect theory.

A broader three-constant elastic density is
\[
f_e(\nabla Q)=\frac{L_1}{2}|\nabla Q|^2+\frac{L_2}{2}Q_{ij,j}Q_{ik,k}+\frac{L_3}{2}Q_{ik,j}Q_{ij,k},
\]
and occurs both in thin-film reduction and in homogenization settings [1501.07339] [2006.08776]. A divergence-penalized variant adds
\[
\frac{k}{2}|\operatorname{Div}Q|^2,
\]
leading, in the small-particle exterior-domain problem, to the linear anisotropic Euler–Lagrange system
\[
-\Delta Q_{ij}-\frac{k}{2}\left(\partial_j\operatorname{Div}Q_i+\partial_i\operatorname{Div}Q_j-\frac{2}{3}\operatorname{div}(\operatorname{Div}Q)\delta_{ij}\right)=0
\]
in \(\Omega=\mathbb R^3\setminus B_1(0)\), with homeotropic data on the colloid and uniaxial far-field alignment [2410.09930].

The elastic sector can also be reformulated to improve well-posedness for disparate Frank constants. A quartic elastic density proposed within generalized Landau–de Gennes theory is
\[
\begin{aligned}
\sigma(Q):={}&\frac{L_1}{2}\left|\left(\frac{s_0}{3}I+Q\right)\operatorname{div}Q\right|^2
+\frac{L_2}{2}\left|\left(\frac{s_0}{3}I+Q\right)\operatorname{curl}Q\right|^2\\
&+\frac{L_3}{2}\left|\left(\frac{2s_0}{3}I-Q\right)\operatorname{div}Q\right|^2
+\frac{L_4}{2}\left|\left(\frac{2s_0}{3}I-Q\right)\operatorname{curl}Q\right|^2.
\end{aligned}
\]
For \(L_i>0\), this density is nonnegative and coercive, and the associated minimization problem is well posed; the construction avoids the unboundedness from below associated with certain cubic elastic theories [1906.09232].

Variationally, existence results in the cited literature rely on the direct method under coercivity and lower semicontinuity assumptions. In the chiral–smectic model, for instance, the admissible space is \(W_Q\times W_{\delta\rho}\), with \(Q\in W^{1,2}_{S_0}(\Omega)\) and \(\delta\rho\in W^{2,2}(\Omega)\), together with Dirichlet anchoring on a bounded simply connected Lipschitz domain. Under explicit inequalities on \(\eta_1,\eta_2,\eta_{24}\), the energy is weakly lower semicontinuous and possesses a global minimizer [2601.12653].

## 3. Oseen–Frank correspondence and asymptotic limits

A persistent theme in the cited work is the reduction of Landau–de Gennes functionals to director theories when the bulk term strongly confines \(Q\) to the uniaxial well. In the chiral–smectic theory, the Oseen–Frank limit is obtained by rescaling \(A,B,C\to \lambda A_0,\lambda B_0,\lambda C_0\) with \(A_0<0\) and letting \(\lambda\to\infty\). Then \(Q_\lambda\) converges to the uniaxial manifold
\[
S_1=\left\{Q=s_+\left(n\otimes n-\frac{I_3}{3}\right)\right\},
\]
and the \(\Gamma\)-limit becomes an Oseen–Frank-type functional in \(n\), coupled to the smectic scalar field. In that reduction, the Frank constants are
\[
K_1=K_3=\frac{s_+^2(\eta_1+\eta_2)}{2},\qquad
K_2=s_+^2\eta_1,\qquad
K_4=\frac{s_+^2(\eta_{24}-\eta_1)}{2}
\]
[2601.12653].

The quartic elastic theory achieves an exact four-constant Oseen–Frank reduction on orientable uniaxial states. For \(Q=s_0(n\otimes n-I/3)\), the elastic density reduces to
\[
F_{OF}(n)=\int_\Omega \left[\frac{K_1}{2}(\operatorname{div}n)^2
+\frac{K_2}{2}((\operatorname{curl}n)\cdot n)^2
+\frac{K_3}{2}|(\operatorname{curl}n)\times n|^2
+\frac{K_2+K_4}{2}\big(\operatorname{tr}(\nabla n)^2-(\operatorname{div}n)^2\big)\right]dx,
\]
with the parameter correspondence
\[
s_0^4L_4=K_2+K_4,\qquad s_0^4(L_i+L_4)=K_i,\quad i=1,2,3.
\]
The corresponding dimensionless family \(\mathcal F_\varepsilon\) \(\Gamma\)-converges to \(\mathcal F_0\), and minimizers converge strongly in \(H^1\) [1906.09232].

Other asymptotic regimes emphasize special elastic hierarchies. In the cholesteric–smectic model, letting the twist modulus \(K_2\to\infty\) with \(K_2+K_4\) bounded forces
\[
\nabla\times n+\sigma n\to 0 \quad \text{in }L^2(\Omega),
\]
so the limiting director is helical,
\[
n_\sigma(z)=(\cos(\sigma z),\sin(\sigma z),0)^T.
\]
Depending on \(d,\lambda_1,\lambda_2,q,\theta_0\), the limiting state is cholesteric \((\delta\rho\equiv 0)\), helical smectic \((\delta\rho\neq 0,\ \theta\equiv 0)\), or SmC* \((\delta\rho\neq 0,\ \theta\neq 0)\) [2601.12653].

The same reduction logic appears in thin-film problems. When thickness and correlation length both vanish, the limiting energies contain leading-order perimeter terms and lower-order vortex terms, so the Landau–de Gennes theory simultaneously produces Allen–Cahn-type interface costs and Ginzburg–Landau-type singular structures [1801.04477].

## 4. Defects, singular sets, and convergence of minimizers

The defect theory of Landau–de Gennes minimizers is strongly shaped by the topology of the vacuum manifold and by the energy scaling regime. In two dimensions, for smooth boundary data valued in \(N\simeq \mathbb{RP}^2\), minimizers in the low-temperature regime are maximally biaxial near singularities: the biaxiality parameter
\[
\beta(Q)=1-\frac{6(\operatorname{tr}(Q^3))^2}{(\operatorname{tr}(Q^2))^3}
\]
reaches the value \(1\), while \(\inf_{\overline\Omega}|Q_\varepsilon|>0\), excluding isotropic melting in the core [1307.8065].

In three dimensions and under the logarithmic energy scaling
\[
E_\varepsilon(Q_\varepsilon)\le M(|\log \varepsilon|+1),
\]
the one-constant quartic model produces line defects. There exists a closed set \(S\subset \overline\Omega\) of finite length such that \(Q_{\varepsilon_n}\) converges strongly in \(H^1_{\mathrm{loc}}(\Omega\setminus S)\) to a locally minimizing harmonic map \(Q_0\in H^1_{\mathrm{loc}}((\Omega\setminus S),N)\). The support of the limiting energy measure is \(S\), the interior part defines a stationary \(1\)-varifold, the density is constant and equal to \(\kappa_*=\frac{\pi}{2}s_*^2\), and \(S\cap \Omega\) is locally a finite union of straight line segments [1501.05236].

Point defects have a different inner asymptotics. For global minimizers in three dimensions, if \(x_0\) is a singular point of the limiting harmonic map \(n_*\), then suitably chosen blow-ups \(Q_{\varepsilon_n}(x_n+\varepsilon_n y)\) converge in \(C^2_{\mathrm{loc}}\) to a local minimizer on \(\mathbb R^3\) whose tangent map at infinity is the standard hedgehog
\[
\Phi(x)=s_+\left(\frac{x}{|x|}\otimes \frac{x}{|x|}-\frac13 \mathrm{Id}\right).
\]
This implies that \(Q_{\varepsilon_n}\) is uniformly approximated by the Oseen–Frank minimizer outside an \(O(\varepsilon_n)\) neighborhood of the defect core [2207.04525].

The sextic theory changes the well topology and, with it, the defect taxonomy. In the bounded-energy regime, minimizers converge to locally minimizing harmonic maps into the biaxial vacuum manifold \(N\simeq S^3/Q_8\), and the singular set is locally finite. In the logarithmically divergent regime, the limiting line singular set \(S_{\mathrm{line}}\) is closed, countably \(\mathcal H^1\)-rectifiable, and in compact subsets is a finite union of closed straight line segments; the limiting density takes the quantized values \(K_*\) or \(2K_*\), where \(K_*=\pi/2\) [2404.00677].

Recent convergence refinements sharpen the classical vanishing-elasticity theory. For local minimizers with uniformly bounded energy and \(L^\infty\) norm, there exists a subsequence \(Q_{\varepsilon_i}\to Q_0\) strongly in \(H^1_{\mathrm{loc}}\), and the convergence is optimal in the sense that
\[
\lim_{i\to\infty}\|Q_{\varepsilon_i}-Q_0\|_{L^p(K)}=0
\qquad\text{for every }1<p<\infty
\]
and compact \(K\Subset \Omega\). The same analysis yields the sharp bulk-energy rate
\[
\int_K \frac{1}{\varepsilon_i^2}f(Q_{\varepsilon_i})\,dx \le C\varepsilon_i,
\]
with matching lower bounds in the hedgehog example [2507.14955].

## 5. Thin films, curved surfaces, and reduced geometries

When the domain is thin, Landau–de Gennes theories admit rigorous dimension reduction. For planar films occupying \(\Omega\times(0,h)\), with strong Dirichlet conditions on the lateral boundary and weak anchoring on the top and bottom, the rescaled energy \(F_\varepsilon\) \(\Gamma\)-converges to a two-dimensional functional \(F_0\) defined on
\[
H^1_g=\{Q\in H^1(\Omega;\mathcal D):Q|_{\partial\Omega}=g\},
\]
where \(\mathcal D\) is the set of \(Q\)-tensors minimizing the leading-order surface energy. In the regime emphasized in that work, \(\mathcal D\) encodes the requirement that the film normal \(\hat z\) is an eigenvector of \(Q\), and the limiting problem reduces, in suitable parameter regimes, to Allen–Cahn- or Ginzburg–Landau-type energies for an effective planar order parameter \(\mathbf p\) [1501.07339].

A different thin-film scaling, in which both thickness and nematic correlation length vanish, yields a singularly perturbed family whose \(\Gamma\)-limit consists of perimeter terms associated with the connected components of the zero set of an effective potential \(W\). In simply connected domains with planar anchoring and degree \(k\), minimizers satisfy the asymptotic law
\[
F_\varepsilon(Q_\varepsilon)=2\int_{\partial\Omega} d_{\sqrt W}(P_1,g(x))\,d\mathcal H^1(x)
+s_*^2\pi k\log\frac1\varepsilon+O(1),
\]
so the reduced model combines an interface energy with a vortex contribution [1801.04477].

On curved substrates, the thin-film limit retains a nontrivial remnant of the normal derivative. The limiting surface energy is
\[
F_0[Q]=\int_{\mathcal M}\left\{f_e^0(\nabla_{\mathcal M}Q,\nu)+\frac{1}{\delta^2}f_{LdG}(Q)+2f_s^{(1)}(Q,\nu)\right\}\,d\mathcal H^2,
\]
where
\[
f_e^0(\nabla_{\mathcal M}Q,\nu)=\min_{G\in\mathcal A} f_e(G\otimes \nu+\nabla_{\mathcal M}Q).
\]
This reduced energy depends explicitly on the surface normal and curvature. On a frustum, the equal-constant reduction leads to an effective director equation with curvature-induced preference for different winding numbers, and the minimizing class changes with the cone angle [1611.03011].

Surface Landau–de Gennes models derived from thin-film limits can preserve distinct three-dimensional features depending on how the normal eigenvalue is handled. In the general surface theory with tangential anchoring, the reduced tensor is decomposed as
\[
Q(q,\beta)=q-\frac{\beta}{2}g+\beta\,n\otimes n,
\]
and the surface energy contains the explicit curvature couplings
\[
L\|B\|^2\operatorname{tr}(q^2)-6LH\beta\langle B,q\rangle+\frac{9L}{2}\|B\|^2\beta^2.
\]
These terms align the in-plane director with principal curvature lines when \(H\beta\neq 0\), and different prescriptions for \(\beta\) produce either three-dimensional first-order or effectively two-dimensional second-order transitions [1912.08685].

## 6. Smectic, chiral, and free-boundary extensions

The classical Landau–de Gennes tensor model can be enlarged by coupling orientational order to positional order. A modified Landau–de Gennes theory for Smectic-A introduces a real scalar \(u\) representing density deviation from the average molecular density. In that model, the authors prove existence and regularity of global minimizers in three dimensions, show that the theory can capture the isotropic–nematic–smectic phase transition as a function of temperature under suitable assumptions, and analyze stable smectic phases on a square domain with tangent boundary conditions. Their asymptotics for \(\lambda\to 0\) and \(\lambda\to\infty\) relate layer number and thickness to the phenomenological parameters [2408.03343].

A more elaborate chiral theory couples the tensor \(Q\) to a real-valued smectic modulation \(\delta\rho(x)=\rho(x)-\rho_0\). Its total free energy is
\[
F(Q,\delta\rho)=\int_\Omega \{f_{el}(Q)+f_{bn}(Q)+f_{bs}(\delta\rho)+f_{layer}(Q,\delta\rho)+f_{angle}(Q,\delta\rho)\}\,dx,
\]
where the layering term is
\[
f_{layer}(Q,\delta\rho)=\lambda_1(\Delta \delta\rho+q^2\delta\rho)^2
\]
and the director–layer coupling is
\[
f_{angle}(Q,\delta\rho)=\lambda_2\left[\operatorname{tr}\!\left(D^2\delta\rho\left(Q+\frac{I_3}{3}\right)\right)+q^2\delta\rho\cos^2\theta_0\right]^2.
\]
The theory yields a complete sequence of temperature-driven symmetry-breaking transitions:
cholesteric \(\to\) helical smectic \(\to\) SmC*, with explicit stability criteria and a pitchfork bifurcation at
\[
d=-2\lambda_2 q^4\cos^4\theta_0
\]
[2601.12653].

Free boundaries can also be incorporated by coupling \(Q\) to a phase field \(\phi\). In the diffuse-interface theory for nematic droplets in isotropic liquid, the energy is
\[
E_\varepsilon=E^{\rm LdG}+\omega_p E_\varepsilon^{\rm mix}+\omega_a E_\varepsilon^{\rm anch}+\omega_v E^{\rm void},
\]
with
\[
E^{\rm mix}_\varepsilon[\phi]=\int_\Omega \left[\varepsilon |\nabla\phi|^2+\varepsilon^{-1}\phi^2(1-\phi)^2\right]dx,
\]
\[
E^{\rm anch}_\varepsilon[Q,\phi]=\int_\Omega \varepsilon \left|\left(Q+\frac{s_+}{3}I\right)\nabla\phi\right|^2 dx,
\]
and
\[
E^{\rm void}[Q,\phi]=\int_\Omega \frac12(1-\phi)^2|Q|^2\,dx.
\]
The model admits global minimizers, satisfies a uniform maximum principle under explicit assumptions, and has a sharp-interface \(\Gamma\)-limit in which the interfacial term blends perimeter and weak anchoring [2407.21437].

At still another level, a complex smectic order parameter \(\psi\) coupled to a director field \(n\) gives rise to surface smectic states. In the asymptotic regime \(q\tau=b\kappa^2\) with \(b>1\), the ground-state energy satisfies
\[
E_{\mathrm{g.st}}(\kappa,q,\tau,K_1,K_2,K_3)=\sqrt{b}\,\kappa\, e_0(b_{GL},\tau)+o(\kappa),
\]
and the smectic order parameter localizes near the boundary where the reduced half-space energy density \(E(b_{GL},\nu)\) is negative. This establishes the existence of a surface smectic state [1602.06241].

## 7. Selected variants, reduced models, and materials-scale applications

Beyond the standard bulk–elastic theory, the cited literature develops several specialized variants that modify defect morphology, solution landscapes, or effective material coefficients.

| Variant | Added structure | Reported consequence |
|---|---|---|
| Divergence-penalized exterior problem | \(\frac{k}{2}|\operatorname{Div}Q|^2\) | Saturn ring remains equatorial; biaxial region grows with \(k\) [2410.09930] |
| Cubic microlattice homogenization | Surface anchoring on a dilute scaffold | Effective bulk coefficients can be shifted, including the phase transition temperature [2006.08776] |
| 2D elastic anisotropy | \(L_2(\operatorname{div}Q)^2\) term | Stable \(WORS\), \(Ring^+\), and \(Constant\) branches, with \(Constant\) favored for large \(\lambda\) and large \(L_2\) [2105.10253] |
| 3D prism reduction | Reduced \(P\)-tensor on polygonal cross-sections | 3D critical points correspond to pathways on 2D solution landscapes [2211.07536] |
| Functional RG | Flow of \(U_k(\rho,\tau)\), cubic/quartic couplings, \(\eta\) | First-order NI transition; \(T_c-T^*=5.85\,K\) [1803.03683] |

The divergence-penalized model is formulated on \(\Omega=\mathbb R^3\setminus B_1(0)\) with homeotropic anchoring at the colloid and a uniform uniaxial far field. In that linearized small-particle regime, the energy is convex, minimizers are unique, and numerical results show that increasing \(k\) only marginally decreases the Saturn ring radius while enlarging the biaxial region around the ring [2410.09930].

Homogenization of a connected cubic microlattice scaffold embedded in a nematic host produces a surface contribution \(f_{\mathrm{hom}}(Q)\) in the effective bulk energy. In the cubic-symmetry case, suitable surface anchoring can realize
\[
f_{\mathrm{hom}}^{LDG}(Q)=(a'-a)\operatorname{tr}(Q^2)-(b'-b)\operatorname{tr}(Q^3)+(c'-c)(\operatorname{tr}(Q^2))^2,
\]
so the effective coefficients \(a',b',c'\) can be tuned through scaffold geometry and surface chemistry. In particular,
\[
T_*^{\mathrm{eff}}=T_*+\frac{a-a'}{a_*},
\]
which shifts the characteristic transition temperature [2006.08776].

In reduced two-dimensional models, elastic anisotropy reorganizes the entire equilibrium landscape. On a truncated square with tangent boundary conditions, the energy
\[
\mathcal F_\lambda[Q]=\int_\Omega \left\{\frac12 |\nabla Q|^2+\frac{L_2}{2}(\operatorname{div}Q)^2+\frac{\lambda^2}{L}f_B(Q)\right\}\,dA
\]
admits a symmetric critical point for every \(L_2>-1\), which is globally stable for small \(\lambda\). Numerical and asymptotic analysis reveal at least five classes of symmetric states: \(WORS\), \(Ring^\pm\), \(Constant\), and \(pWORS\). The \(Constant\) branch becomes energetically preferable for large domain size and large \(L_2\) [2105.10253].

A reduced three-dimensional prism model transfers this landscape logic into polygonal cross-sections. At the special temperature \(A=-B^2/(3C)\), the reduced functional
\[
F[P]=\int_{V_K}\left\{\frac12 |\nabla_{xy}P|^2+\frac{1}{2h^2}|P_{,z}|^2+\frac{\bar\lambda^2}{2C}\left(-\frac{B^2}{4C}\operatorname{tr}P^2+\frac{C}{4}(\operatorname{tr}P^2)^2\right)\right\}dV
\]
supports mixed three-dimensional critical points that correspond to paths between two-dimensional critical points on the cross-section. The numerical examples on cuboids and hexagonal prisms show explicitly how multistability in three dimensions can be tailored by two-dimensional solution landscapes [2211.07536].

At a more coarse-grained statistical level, a functional renormalization-group treatment of the Landau–de Gennes model uses the flow of the effective potential \(U_k(\rho,\tau)\), together with running cubic and quartic couplings and anomalous dimension, to study the nematic–isotropic transition. The computed effective potential displays a first-order transition, and the resulting estimate
\[
T_c-T^*=5.85\,K
\]
substantially improves on earlier values discussed in that work’s comparison [1803.03683]. A plausible implication is that the Landau–de Gennes model is not only a continuum variational theory but also a viable starting point for fluctuation analyses that probe the limits of mean-field predictions.

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