Self-consistent field theory of semiflexible nematics: Density-nematic coupling, anisotropic elasticity, and defect core sizes

Technical Summary: Self-consistent field theory of semiflexible nematics

Problem Statement
While the linear response of small-molecule thermotropic liquid crystals is well-understood, the behavior of lyotropic liquid crystals—comprising polymers, micellar aggregates, and biofilaments—presents significant theoretical and experimental challenges. A central issue is the origin of large elastic anisotropy observed in lyotropic systems and the coupling between chain density and nematic order. Existing microscopic theories, largely based on the Onsager model, struggle to simultaneously address elasticity, topological defects, and density–order coupling in spatially varying semiflexible nematics. Furthermore, macroscopic continuum theories (Oseen–Frank and Landau–de Gennes) typically introduce phenomenological coefficients without explicit connection to underlying molecular statistics, leaving the microscopic origins of elastic anisotropy and defect core structures unexplored.

Methodology
The authors employ a Self-Consistent Field Theory (SCFT) framework modeling semiflexible nematics as wormlike chains. The system is defined by a Hamiltonian incorporating:

1. Bending energy: Governed by the persistence length ($l_p$) and contour length ($L$).
2. Maier–Saupe (MS) interactions: Anisotropic orientational interactions ($u_2$) driving nematic ordering.
3. Isotropic excluded volume: A solvent-mediated interaction ($u_0$) representing lyotropic conditions (implicitly treating solvent effects via Flory–Huggins parameters).

The study proceeds in three stages:

• Uniform State Analysis: The authors derive an effective free energy functional dependent on chain segment density ($\rho$) and the uniaxial order parameter ($S$). This allows for the determination of isotropic-nematic coexistence (binodal) and stability limits (spinodal) boundaries.
• Linear Response Analysis: Weak external fields are introduced to probe the system's response to perturbations in density, uniaxial order, and director orientation. This yields the inverse susceptibility matrix, from which Frank elastic constants ($K_1, K_2, K_3$) and correlation lengths are extracted.
• Defect Structure Simulation: The full SCFT equations are solved numerically for +1/2 disclinations to analyze director profiles and core sizes, comparing rigid-chain limits with semiflexible regimes.

Key Contributions and Results

• Density–Nematic Coupling and Stability:
  The derived effective free energy reveals a strong coupling between density ($\rho$) and the degree of order ($S$). The authors demonstrate that strong density–nematic coupling can destabilize the nematic state. The location of the spinodal instability depends critically on the ratio of excluded volume to nematic interaction strength ($u_0/u_2$). In contrast, director distortions couple to density and order variations only at higher orders, remaining effectively decoupled in the linear response regime.

• Elastic Anisotropy and Flexibility:
  The study quantifies how chain flexibility ($L/l_p$) alters Frank elastic constants.

  • Rigid Chains: $K_3$ (bend) increases rapidly with alignment strength, while $K_1$ (splay) and $K_2$ (twist) remain comparable.
  • Semiflexible Chains: Increasing flexibility strongly suppresses twist ($K_2$) and bend ($K_3$) elasticity, while splay elasticity ($K_1$) is affected comparatively weakly.
  • Crossover: This differential suppression leads to a crossover from bend-dominated to splay-dominated elasticity as flexibility increases, evidenced by a reversal in the ratio $K_1/K_3$ (from $<1$ to $>1$).

• Correlation Lengths and Anisotropy:
  Linear response analysis reveals that while the diagonal gradient coefficients of the inverse susceptibility suggest shorter correlation lengths parallel to the director ($K_{\parallel} < K_{\perp}$), the coupling between density and order renormalizes the effective coefficients. This results in an anisotropy where the correlation length parallel to the director is actually larger ($\xi_{\parallel} > \xi_{\perp}$). This finding aligns with experimental observations of isotropic-nematic interfaces and defect core structures.

• Defect Core Sizes:
  Numerical simulations of +1/2 disclinations show that the core size is anisotropic, with the radius being larger along the direction where the density gradient is parallel to the director. The study establishes that the characteristic defect core size is proportional to the equilibrium correlation length ($\xi_Q$), consistent with Landau–de Gennes scaling. The authors note that while lyotropic chromonics often exhibit anomalously large cores (potentially due to aggregate scission), their current model, which excludes polydispersity, captures the scaling relation between core size and correlation length driven by density fluctuations.

• Spinodal Behavior:
  As the system approaches the nematic spinodal (controlled by varying $u_0$), the correlation lengths diverge following the scaling $\xi \sim (u_0 - u_{0,i})^{-1/2}$. However, in defect configurations, this divergence is preempted by the local formation of an isotropic domain (tactoid) at the core, which suppresses nematic order before the spinodal is reached.

Significance
The paper provides a unified microscopic framework linking molecular flexibility, density–nematic order coupling, and elastic anisotropy in semiflexible nematics. By deriving Frank elastic constants and correlation lengths directly from chain statistics within the SCFT formalism, the work elucidates the microscopic origins of the large elastic anisotropy observed in lyotropic systems. The results validate the use of continuum theories (like Landau–de Gennes) for describing defect structures while highlighting the necessity of incorporating density fluctuations and chain flexibility to accurately predict elastic responses and stability limits. The authors conclude that their framework offers a basis for future extensions to self-assembled, polydisperse chains and reversible aggregation, though they do not propose specific new experimental protocols or applications beyond the scope of the current theoretical analysis.