Nematic order parameter ($Q_{ij}$)

A nematic order parameter measures orientational alignment without distinguishing the two ends of a molecular axis. For unit molecular axes $\mathbf{u}$ in three dimensions, a common tensor convention isIt is symmetric and traceless, and is unchanged by $\mathbf{u}\mapsto -\mathbf{u}$. An isotropic distribution gives $Q=0$ even though individual molecules retain definite axes.

For a uniaxial nematic liquid crystal, this convention gives $Q=S(\mathbf{nn}-I/3)$, where $\mathbf{n}$ is the unoriented nematic director and $S=⟨(3{\cos }^2\theta -1)/2⟩$. Perfect alignment has $S=1$ and eigenvalues $2/3,-1/3,-1/3$. Other normalizations multiply $Q$ by a constant, so coefficients in a Landau-Ginzburg theory must use the same convention. The tensor representation handles head–tail symmetry and permits alignment to vanish at defect cores.

For a narrow nematic liquid crystal orientation distribution about its unoriented nematic director, expand $S=⟨(3{\cos }^2\theta -1)/2⟩=1-3⟨{\theta }^2⟩/2+O(⟨{\theta }^4⟩)$. This proves the displayed root-mean-square width. It does not determine the full distribution or a full width at half maximum. If the two small transverse tilt components are independent equal-variance Gaussians with variance ${\sigma }^2$, then $⟨{\theta }^2⟩\simeq 2{\sigma }^2$, so $\sigma \simeq \sqrt{(1-S)/3}$.

For scalar nematic order parameter $S=⟨(3{\cos }^2\theta -1)/2⟩$, the displayed angle has ${\cos }^2{\theta }_{\mathrm{eff}}=⟨{\cos }^2\theta ⟩$. It summarizes a second moment, not the literal mean angle. For example an isotropic orientation distribution has $S=0$, giving ${\theta }_{\mathrm{eff}}=\arccos (1/\sqrt{3})$; the average acute angle to its unoriented axis is instead one radian, obtained by integrating $\theta \sin \theta$ over $[0,\pi /2]$.

A uniaxial nematic order parameter has one distinguished nematic director and equal eigenvalues in its perpendicular plane. In the displayed normalization the eigenvalues are $\lambda ,-\lambda /2,-\lambda /2$, giving $\mathrm{Tr}{Q}^2=3{\lambda }^2/2$ and $\mathrm{Tr}{Q}^3=3{\lambda }^3/4$. Positive and negative $\lambda$ describe prolate and oblate orientational anisotropy, respectively; reversing the nematic director leaves the tensor unchanged.

For $S=⟨3\mathbf{s}{\mathbf{s}}^{\mathsf{T}}-I⟩=Q(3\mathbf{n}{\mathbf{n}}^{\mathsf{T}}-I)$ and unit molecular orientations $\mathbf{s}$, contraction with $\mathbf{n}$ gives $2Q=3⟨(\mathbf{s}\cdot \mathbf{n}{)}^2⟩-1$. Since the squared projection is between zero and one, the displayed bounds follow. The eigenvalues in this normalization are $2Q,-Q,-Q$. A differently normalized nematic order parameter has correspondingly rescaled bounds; its numerical coefficients cannot be transferred unchanged.

Under rotation by a rotation matrix, a symmetric second-rank tensor transforms by conjugation, so its traces of powers are invariant. For a three-dimensional traceless second-rank tensor, the Cayley-Hamilton theorem implies $\mathrm{Tr}{Q}^4=(\mathrm{Tr}{Q}^2{)}^2/2$. In two dimensions its eigenvalues are $s,-s$, so $\mathrm{Tr}{Q}^3=0$. These identities organize the Landau-de Gennes free energy.

The Landau-de Gennes free energy is a rotational symmetry preserving Landau free energy of the nematic order parameter. One three-dimensional normalization through fourth order isThe cubic invariant generically produces a first-order phase transition. Its absence in two dimensions is an algebraic identity, not a guarantee that a mean-field approximation captures all fluctuations.

The leading scalar coupling of an electric field to the nematic order parameter is quadratic in the field and linear in the tensor. For uniaxial nematic order, it is $-3\chi \lambda E^2({\cos }^2\vartheta -1/3)/2$. When $\chi \lambda >0$, minimization aligns the nematic director with the field axis and adds $-\chi E^2\lambda$ to the scalar Landau free energy.

For a scalar Landau-de Gennes free energy with $\bar{b}>0$, $\bar{c}<0$, and coupling $-\chi E^2\lambda$ with $\chi >0$, the critical endpoint of a quartic Landau free energy is reached at the displayed field strength and ${\bar{a}}_c=3{\bar{c}}^2/(8\bar{b})$, ${\lambda }_c=-\bar{c}/(4\bar{b})$. It terminates the field-biased first-order phase transition when both field and quadratic coefficient can be tuned. This is a quartic small-field mean-field approximation; higher terms may alter its location.