Anisotropic Nematic Liquid Crystals

• Anisotropic nematic liquid crystals are orientationally ordered fluids characterized by directional elastic and viscous properties that influence phase behavior.
• Continuum theories like Landau–de Gennes and Oseen–Frank model their anisotropic moduli, guiding experimental quantification of defect dynamics and phase transitions.
• Tunable anisotropy in these materials enables precise control over tactoid morphology, defect core splitting, and practical applications in optics and microfluidics.

Anisotropic nematic liquid crystals are orientationally ordered soft-matter phases defined by a directionally dependent (anisotropic) elasticity, viscosity, or related physical property. The anisotropy manifests through distinct moduli, transport coefficients, defect energetics, and hydrodynamics that couple to the molecular director field, and has far-reaching implications for equilibrium phase transitions, interface morphologies, topological singularities, and dynamic response functions. This entry reviews the continuum statistical, molecular, and hydrodynamic theories of anisotropic nematics, emphasizing how elastic and viscous anisotropy are represented, their impact on defect and interface structure, and quantitative connections to experiment.

1. Theoretical Formulation of Anisotropy

Anisotropic nematic liquid crystals are most generally described by continuum order parameters capturing the broken rotational symmetry and spatially varying molecular arrangement. The Landau–de Gennes (LdG) Q-tensor theory posits a symmetric, traceless $3\times3$ tensor $Q(\mathbf{x})$, representing the second moment of the molecular axis distribution. In uniaxial nematics, $Q = S(n\otimes n - I/3)$, with $n\in S^2$ the director and $S$ the scalar order parameter. The isotropic state corresponds to $Q=0$.

The free energy functional comprises both bulk (thermodynamic) and elastic (gradient) contributions,

$$E[Q] = \int_\Omega \left\{ f_\text{elastic}(Q,\nabla Q) + W_\text{bulk}(Q) \right\} dV,$$

where the bulk energy is a fourth-order polynomial in $Q$ invariants, e.g.,

$$W(Q) = 3a\,\mathrm{Tr}(Q^2) - 2b\,\mathrm{Tr}(Q^3) + \tfrac14[\mathrm{Tr}(Q^2)]^2,$$

and the elastic part is expressed as a sum of quadratic invariants in gradients of $Q$: $$f_\text{elastic}(Q,\nabla Q) = \tfrac12 L_1\,\partial_k Q_{ij}\,\partial_k Q_{ij} + \tfrac12 L_2\,(\partial_j Q_{ij})(\partial_k Q_{ik}) + \tfrac12 L_3\,(\partial_j Q_{ik})(\partial_k Q_{ij}).$$ Anisotropy enters via the three elastic constants $L_1$, $L_2$, and $L_3$. For $L_2\gg L_1,L_3$, splay deformations are strongly penalized relative to bend and twist, modeling highly anisotropic materials such as chromonic nematics (Golovaty et al., 2019).

The Oseen–Frank director theory is recovered in the uniaxial case and frames director distortions in terms of splay ($K_1$), twist ($K_2$), and bend ($K_3$) elastic constants, with anisotropy when $K_1\neq K_2\neq K_3$ (Borthagaray et al., 2020). The relationship between $L$ parameters and $K_i$ is material- and symmetry-dependent, and mapping from LdG to Oseen–Frank under uniaxial reduction yields explicit correspondences.

2. Phase Transitions and Tactoid Morphology

Anisotropic elasticity plays a decisive role in the first-order isotropic-to-nematic phase transition as well as the geometry of phase coexistence domains (tactoids). The LdG bulk potential admits coexisting isotropic ($Q=0$) and nematic ($|Q|>0$) minima for $0 < a < b^2/4$, with equal well depth at $a_0=2b^2/9$ (Golovaty et al., 2019). For $a<a_0$ (finite undercooling), nucleation and growth of nematic tactoids in an isotropic matrix or vice versa are observed.

In the presence of strong elastic anisotropy ($L_2\gg L_1,L_3$), tactoid interfaces exhibit noncircular, elongated shapes with aspect ratios fixed by the balance between splay energy and interfacial tension. The director couples to the interface via the $L_2$ term, resulting in equilibrium shapes where the director lies tangent to the boundary. Numerical solutions to the gradient-flow LdG equations (in the thin-film limit) naturally recover these morphologies, with interface width scaling as $O(\varepsilon)$, $\varepsilon = L_1/L_2 \ll 1$.

Upon tactoid shrinkage, integer-strength ($\pm 1$) disclination defects do not persist; instead, the Q-tensor formalism predicts splitting into pairs of half-charge ($\pm 1/2$) disclinations, in quantitative agreement with experiment—a phenomenon inaccessible to pure director models. The core radii of emergent $\pm 1/2$ disclinations scale as $C \varepsilon^{1/2}$ near onset. Coalescence of nematic islands leads to net half-charge defect trapping, a signature of the Kibble mechanism (Golovaty et al., 2019).

3. Defects, Biaxial Cores, and Elastic Anisotropy

Topological defects in nematic order—disclinations—display rich structure in the presence of elastic anisotropy. Both LdG and self-consistent molecular-field theories accommodate anisotropic elastic energies of the form (Schimming et al., 2020): $$f_e(Q,\nabla Q) = L_1\,\partial_k Q_{ij}\partial_k Q_{ij} + L_2\,(\partial_j Q_{ij})(\partial_k Q_{ik}) + L_3\,Q_{kl}\,\partial_k Q_{ij}\partial_l Q_{ij}.$$ Mapping to Oseen–Frank’s elastic constants yields $K_{33}-K_{11} = (16/3) S^3 L_3$ and $K_{11}+K_{33}\propto S^2(2L_1 + L_2)$. The $L_3$ term removes degeneracy between bend and splay.

Near defect cores, the Q-tensor is necessarily spatially inhomogeneous, and the solution exhibits a crossover: isotropic ($S=P=0$) at the core, a region of finite biaxiality $P(r)$ at intermediate distances, and asymptotic uniaxial order far from the core. The anisotropy modifies the director distortion near the core as predicted by Dzyaloshinsky’s solution: the director angle for a $+1/2$ defect reads $$\phi(\theta)\approx \tfrac12\theta + (\varepsilon/6)\sin\theta + \mathcal{O}(\varepsilon^2)$$, $\varepsilon=(K_{33}-K_{11})/(K_{33}+K_{11})$ (Schimming et al., 2020).

The size of cores is controlled primarily by the surface tension term $|\nabla S|^2$ (through $2L_1+L_2$), and increasing $L_2$ increases the core size without introducing further anisotropy. In active nematics, the anisotropic polarization $p_i = \partial_j Q_{ij}$ generated by spatially inhomogeneous $Q$ is proportional to the active force density, and elastic anisotropy thus imparts nontrivial vector-flow patterns around defects, directly influencing active-stress distributions (Schimming et al., 2020).

4. Two-Dimensional Elastic Anisotropy and Solution Landscapes

In geometrically confined nematic systems, elastic anisotropy substantially enriches the space of critical points and stable textures. In the two-dimensional LdG model, anisotropy is typically parametrized by $L_2$: $$\mathcal{F}_\lambda[Q] = \int_{\Omega} \left\{ \tfrac12 |\nabla Q|^2 + \tfrac{L_2}{2} (\mathrm{div} Q)^2 \right\} + \frac{\lambda^2}{L} \left(\tfrac{A}{2}\mathrm{tr} Q^2 - \tfrac{B}{3}\mathrm{tr} Q^3 + \tfrac{C}{4}(\mathrm{tr} Q^2)^2\right) dA$$ (Han et al., 2021).

Depending on $L_2$ and the domain size $\lambda$, at least five classes of symmetric critical points exist: WORS (well-order-reconstruction solution), Ring$^+$ (central $+1$ defect), Ring$^-$, Constant (out-of-plane order), and pWORS (periodic alternation of $\pm 1/2$ defects). Increasing $L_2$ stabilizes Ring$^+$ and the Constant solution, and in the large-$\lambda$, large-$L_2$ regime the system transitions to a uniform out-of-plane texture. Bifurcation diagrams reveal a hierarchy of bifurcations and multistability as $L_2$ is tuned, providing a mechanism for device-level control over defect and alignment states (Han et al., 2021).

5. Dynamic Response and Hydrodynamic Anisotropy

The hydrodynamics of nematic liquid crystals incorporates viscous and elastic anisotropy at the level of both continuum and mesoscopic models. In continuum nematodynamics, the anisotropic viscous stress is given by the Leslie tensor $\nu_{ijkl}$, and friction may also be substrate-directional. The extension of the nematic hydrodynamic theory to include continuous update of the shear-stress-free configuration leads to dynamical models that yield the classical Ericksen–Leslie equations in the low-frequency limit (Biscari et al., 2015).

Importantly, in models incorporating anisotropy in relaxation times (shear modes parallel and perpendicular to the director), new algebraic relations among the Leslie viscosity coefficients emerge, such as

$$\frac{\alpha_2}{\alpha_3} = \frac{\alpha_4+\alpha_5}{\alpha_4+\alpha_6} = a^3,$$

where $a$ is the energetic ellipsoid aspect-ratio (Biscari et al., 2015). These predictions have been confirmed in molecular dynamics simulations and experimental data sets.

Anisotropic viscosity and friction profoundly affect defect dynamics in active nematics. In two-dimensional turbulence, introducing substrate-tuned viscosity anisotropy leads to global nematic order of $+1/2$ defects, while friction anisotropy does not. The degree of defect alignment is governed by both the anisotropy parameter and activity, with maximal ordering at intermediate values (Pearce, 2018).

6. Anisotropic Wave Propagation and Acousto-Optic Response

The propagation of acoustic waves in nematic liquid crystals is directionally dependent due to anisotropic elasticity and relaxation properties. A first-gradient continuum theory employing a hyperelastic anisotropic response with an evolving relaxed configuration successfully captures both the anisotropy of sound speed and its frequency dependence (Biscari et al., 2013). For small-amplitude longitudinal waves, the phase velocity

$$s(\omega, \theta) = v_0\left[1 + \eta\,f(\omega\tau)\left(2/3 - \alpha_1 + 3\alpha_1^2/2 + 3\alpha_1\cos^2\theta\right)\right],$$

with $\eta = \mu_0/p_I \ll 1$, $\tau$ the relaxation time, and $f(x) = x^2/(1 + x^2)$. The anisotropy of sound velocity and attenuation exhibits linear dependence on frequency in the 2–14 MHz band, matching experimental observations. Director oscillation amplitudes induced by acoustic waves are typically extremely small due to the weak molecular anisotropy, requiring high ultrasonic intensities for measurable acousto-optic effects (Turzi, 2014).

7. Applications and Device-Level Implications

The ability to program and tune elastic and viscous anisotropy in nematic devices enables a wide range of practical functions. In microfluidics, patterning the surface anchoring or applying an external field can create spatially varying anisotropic viscosity maps that robustly direct flow, enabling guided transport in Hele–Shaw geometries (Cousins et al., 2024). In the optical domain, devices exploiting the large birefringence (optical anisotropy) of nematics sandwiched with cholesteric layers operate as electrically tunable polarization rotators, with full modulation depth and linear-in-voltage control of output polarization (Gevorgyan et al., 2017).

From a computational perspective, mesoscopic stochastic rotation dynamics models faithfully capture non-Newtonian rheology (shear banding, anisotropic viscosities, flow alignment) and topological defect dynamics, as required to model realistic active and passive nematic fluids (Lee et al., 2015).

Anisotropic nematic elasticity and viscosity underlie the design and function of both passive display materials and active matter systems—the ability to control, predict, and exploit anisotropic energy landscapes is central to the architecture of advanced liquid crystal technologies.

---

Key References

• Landau–de Gennes tensorial approaches to anisotropic elasticity and phase transitions (Golovaty et al., 2019, Han et al., 2021)
• Molecular field theory and defect core structure (Schimming et al., 2020)
• Hydrodynamic relaxation and viscosity anisotropy (Biscari et al., 2015)
• Anisotropic Hele–Shaw flow applications (Cousins et al., 2024)
• Acoustic wave propagation and frequency-dependent anisotropy (Biscari et al., 2013)
• Active nematics and substrate-tuned anisotropy (Pearce, 2018)