A Landau-de Gennes Type Theory for Cholesteric-Helical Smectic-Smectic C* Liquid Crystal Phase Transitions

We present a rigorous mathematical analysis of a modified Landau-de Gennes (LdG) theory modeling temperature-driven phase transitions between cholesteric, helical smectic, and smectic C* phases. This model couples a tensor-valued order parameter (nematic orientational order) with a real-valued order parameter (smectic layer modulation). We establish the existence of energy minimizers of the modified LdG energy in three dimensions, subject to Dirichlet conditions, and rigorously analyze the energy minimizers in two asymptotic limits. First, in the Oseen–Frank limit, we show that the global minimizer strongly converges to a minimizer of the Landau-de Gennes bulk energy. Second, in the limit of dominant elastic constants, we prove that the global minimizers converge to a classical helical director profile. Finally, through stability analysis and bifurcation theory, we derive the complete sequence of symmetry-breaking transitions with decreasing temperature-from the cholesteric phase (with in-plane twist and no layering) to an intermediate helical smectic phase (with in-plane twist and layering), and ultimately to the smectic C* phase (with out-of-plane twist and layering). These theoretical results are supported by numerical simulations.

💡 Research Summary

This paper presents a rigorous mathematical treatment of a modified Landau‑de Gennes (LdG) theory designed to capture temperature‑driven phase transitions among three chiral liquid‑crystal phases: the cholesteric (chiral nematic), an intermediate helical smectic, and the smectic C* phase. The model couples a symmetric, traceless Q‑tensor (describing nematic orientational order) with a real‑valued scalar field δρ (describing the deviation of the mass density from its average, i.e., the smectic layering). The total free‑energy functional (Eq. 1) consists of five contributions: (i) an elastic term f_el that includes curl, divergence, and higher‑order gradient couplings with elastic constants η₁, η₂, η₂₄ and the chiral pitch σ; (ii) a bulk nematic term f_bn (the classic quartic polynomial in Q with temperature‑dependent coefficient A); (iii) a bulk smectic term f_bs (a quartic polynomial in δρ with temperature‑dependent coefficient d); (iv) a layer‑formation term f_layer that penalises deviations from a sinusoidal density modulation with wave number q; and (v) an angle‑coupling term f_angle that enforces a preferred angle θ₀ between the layer normal and the director.

The authors first establish the existence of global minimizers of the functional on a bounded, simply‑connected Lipschitz domain Ω⊂ℝ³ under strong (Dirichlet) anchoring conditions Q=Q_bc, δρ=δρ_bc on ∂Ω. By constructing appropriate Sobolev spaces (W¹,² for Q and W²,² for δρ) and proving three key lemmas—(1) coercivity of the elastic quadratic form under the parameter constraints η₁>0, 0<η₂₄<3η₁, 5η₁+10η₂−9η₂₄>0; (2) an L^p estimate for the second‑order operator governing δρ; and (3) weak lower semicontinuity of each energy contribution—the direct method of the calculus of variations yields a minimizer (Q*,δρ*).

Two asymptotic regimes are then analyzed. In the Oseen‑Frank limit (large domain, dominant bulk coefficients), the elastic contribution becomes negligible compared with the bulk term f_bn. The authors prove strong convergence of the minimizer’s Q‑field to a uniaxial configuration Q≈s₊(n⊗n−I/3), where s₊ is the temperature‑dependent scalar order parameter obtained from the minimisation of f_bn. This recovers the familiar “vector‑valued director” picture often employed in classical LdG analyses. In the opposite limit of dominant elastic constants (η_i→∞), the elastic energy forces the Q‑field to adopt the classical helical cholesteric profile n(z)=(cos (2πz/p), sin (2πz/p),0), confirming that the modified model contains the standard cholesteric theory as a special case.

The core of the paper is a detailed bifurcation and stability analysis of the temperature‑driven transitions. The smectic bulk coefficient d=α₂(T−T₂*) governs the onset of layering: for d>0 the smectic order parameter vanishes (cholesteric phase), while d<0 favours a non‑zero sinusoidal modulation (layered phases). Linear stability calculations identify a critical temperature T_c at which the cholesteric solution loses stability. Below T_c, the coupling terms f_layer and f_angle become active. When the preferred angle θ₀=0 (smectic A‑like), the system settles into an intermediate helical smectic state: the director retains the in‑plane twist of the cholesteric, but a sinusoidal density wave appears. When θ₀>0 (smectic C*), an additional symmetry breaking occurs: the director tilts out of the layer plane, producing the characteristic conical helix of the smectic C* phase. Using Lyapunov‑Schmidt reduction and the Crandall‑Rabinowitz theorem, the authors derive explicit critical values of the coupling parameters λ₁, λ₂ and the wave number q that dictate the sequence of bifurcations. The resulting phase diagram shows a cascade: cholesteric → helical smectic → smectic C* as temperature decreases.

Numerical simulations complement the analytical work. The authors discretise the functional on a three‑dimensional grid using finite differences and minimise it via a gradient‑descent algorithm while varying the temperature parameter. The computed bifurcation diagram matches the theoretical predictions, confirming the role of λ₁q² and λ₂cos²θ₀ in shifting the critical points. Visualisations of the minimisers illustrate the evolution of both Q‑tensor eigenvectors and the density modulation δρ across the three phases, highlighting defect‑free configurations in the cholesteric regime and the emergence of layered structures with tilted directors in the smectic C* regime.

In conclusion, the paper provides a comprehensive, mathematically rigorous foundation for a tensorial Landau‑de Gennes model that simultaneously captures orientational and positional ordering in chiral liquid crystals. It proves existence of minimisers, justifies common simplifying limits, and fully characterises the symmetry‑breaking cascade driven by temperature. The work bridges phenomenological modeling with rigorous analysis, opening pathways for future extensions that incorporate external fields, confinement effects, and dynamical (time‑dependent) Ginzburg‑Landau equations.