Strong anchoring boundary conditions in nematic liquid crystals: Higher-order corrections to the Oseen-Frank limit and a revised small-domain theory

Strong anchoring boundary conditions are conventionally modelled by imposing Dirichlet conditions on the order parameter in Landau-de Gennes theory, neglecting the finite surface energy of realistic anchoring. This work revisits the strong anchoring limit for nematic liquid crystals in confined two-dimensional domains. By explicitly retaining a Rapini-Papoular surface energy and adopting a scaling where the extrapolation length $l_{ex}$ is comparable to the coherence length $\xi$, we analyse both the small-domain ($\ep = h/\xi\to 0$; $h$ is the domain size) and Oseen-Frank $(\ep \to \infty$) asymptotic regimes. In the small-domain limit, the leading-order equilibrium solution is given by the average of the boundary data, which can vanish in symmetrically frustrated geometries, leading to isotropic melting. In the large-domain limit, matched asymptotic expansions reveal that surface anchoring introduces an $O(1/\ep)$ correction to the director field, in contrast to the $O(1/\ep^2)$ correction predicted by Dirichlet conditions. The analysis captures the detailed structure of interior and boundary defects, showing that mixed (Robin-type) boundary conditions yield smoother defect cores and more physical predictions than rigid Dirichlet conditions. Numerical solutions for square and circular wells with tangential anchoring illustrate the differences between the two boundary condition treatments, particularly in defect morphology. These results demonstrate that a consistent treatment of anchoring energetics, together with stability considerations, is essential for accurate modelling of nematic equilibria in micro- and nano-scale confined geometries.