Theory and simulation of the confined Lebwohl-Lasher model

Bistable director alignments of nematic liquid crystals confined in frustrated substrates

We studied in-plane bistable alignments of nematic liquid crystals confined by two frustrated surfaces by means of Monte Carlo simulations of the Lebwohl-Lasher spin model. The surfaces are prepared with orientational checkerboard patterns, on which the director field is locally anchored to be planar yet orthogonal between the neighboring blocks. We found the director field in the bulk tends to be aligned along the diagonal axes of the checkerboard pattern, as reported experimentally [J.-H. Kim et al., Appl. Phys. Lett. 78, 3055 (2001)APPLAB0003-695110.1063/1.1371246]. The energy barrier between the two stable orientations is increased, when the system is brought to the isotropic-nematic transition temperature. Based on an elastic theory, we found that the bistability is attributed to the spatial modulation of the director field near the frustrated surfaces. As the block size is increased and/or the elastic modulus is reduced, the degree of the director inhomogeneity is increased, enlarging the energy barrier. We also found that the switching rate between the stable states is decreased when the block size is comparable to the cell thickness.

Mean-field density functional theory of a nanoconfined classical, three-dimensional Heisenberg fluid. I. The role of molecular anchoring

In this work, we employ classical density functional theory (DFT) to investigate for the first time equilibrium properties of a Heisenberg fluid confined to nanoscopic slit pores of variable width. Within DFT pair correlations are treated at modified mean-field level. We consider three types of walls: hard ones, where the fluid-wall potential becomes infinite upon molecular contact but vanishes otherwise, and hard walls with superimposed short-range attraction with and without explicit orientation dependence. To model the distance dependence of the attractions, we employ a Yukawa potential. The orientation dependence is realized through anchoring of molecules at the substrates, i.e., an energetic discrimination of specific molecular orientations. If the walls are hard or attractive without specific anchoring, the results are "quasi-bulk"-like in that they can be linked to a confinement-induced reduction of the bulk mean field. In these cases, the precise nature of the walls is completely irrelevant at coexistence. Only for specific anchoring nontrivial features arise, because then the fluid-wall interaction potential affects the orientation distribution function in a nontrivial way and thus appears explicitly in the Euler-Lagrange equations to be solved for minima of the grand potential of coexisting phases.

Capillary and winding transitions in a confined cholesteric liquid crystal

We consider a Lebwohl-Lasher model of chiral particles confined in a planar cell (slit pore) under different boundary conditions, and solve it using mean-field theory. The phase behaviour of the system with respect to temperature and pore width is studied. Two phenomena are observed: (i) an isotropic-cholesteric transition, which exhibits an oscillatory structure with respect to pore width, and (ii) an infinite set of winding transitions caused by commensuration effects between cholesteric pitch and pore width. The latter transitions have been predicted and analysed by other authors for cholesterics confined in a fixed pore and subjected to an external field promoting the uniaxial nematic phase; here we induce winding transitions solely from geometry by changing the pore width at zero external field (a setup recently explored in atomic-force microscopy experiments). In contrast with previous studies, we obtain the phase diagram in the temperature vs. pore width plane, including the isotropic-cholesteric transition, the winding transitions and their complex relationship. In particular, the structure of winding transitions terminates at the capillary isotropic-cholesteric transition via triple points where the confined isotropic phase coexists with two cholesterics with different helix indices. For symmetric and asymmetric monostable plate anchorings the phase diagrams are qualitatively similar.

Domain walls in two-dimensional nematics confined in a small circular cavity

Using Monte Carlo simulation, we study a fluid of two-dimensional hard rods inside a small circular cavity bounded by a hard wall, from the dilute regime to the high-density, layering regime. Both planar and homeotropic anchoring of the nematic director can be induced at the walls through a free-energy penalty. The circular geometry creates frustration in the nematic phase and a polar-symmetry configuration with a distorted director field plus two +1/2 disclinations is created. At higher densities, a quasi-uniform structure is observed with a (minimal) director distortion which is relaxed via the formation of orientational domain walls. This novel structure is not predicted by elasticity theory and is similar to the step-like structures observed in three-dimensional hybrid slit pores. We speculate that the formation of domain walls is a general mechanism to relax elastic stresses under the conditions of strong surface anchoring and severe spatial confinement.

Liquid-crystal patterns of rectangular particles in a square nanocavity

Using density-functional theory in the restricted-orientation approximation, we analyze the liquid-crystal patterns and phase behavior of a fluid of hard rectangular particles confined in a two-dimensional square nanocavity of side length H composed of hard inner walls. Patterning in the cavity is governed by surface-induced order as well as capillary and frustration effects and depends on the relative values of the particle aspect ratio κ≡L/σ, with L the length and σ the width of the rectangles (L≥σ), and cavity size H. Ordering may be very different from bulk (H→∞) behavior when H is a few times the particle length L (nanocavity). Bulk and confinement properties are obtained for the cases κ=1, 3, and 6. In bulk the isotropic phase is always stable at low packing fractions η=Lσρ_{0} (with ρ_{0} the average density) and nematic, smectic, columnar, and crystal phases can be stabilized at higher η depending on κ: For increasing η the sequence of isotropic to columnar is obtained for κ=1 and 3, whereas for κ=6 we obtain isotropic to nematic to smectic (the crystal being unstable in all three cases for the density range explored). In the confined fluid surface-induced frustration leads to fourfold symmetry breaking in all phases (which become twofold symmetric). Since no director distortion can arise in our model by construction, frustration in the director orientation is relaxed by the creation of domain walls (where the director changes by 90^{∘}); this configuration is necessary to stabilize periodic phases. For κ=1 the crystal becomes stable with commensurate transitions taking place as H is varied. These transitions involve structures with different number of peaks in the local density. In the case κ=3 the commensurate transitions involve columnar phases with different number of columns. In the case κ=6 the high-density region of the phase diagram is dominated by commensurate transitions between smectic structures; at lower densities there is a symmetry-breaking isotropic to nematic transition exhibiting nonmonotonic behavior with cavity size. Apart from the present application in a confinement setup, our model could be used to explore the bulk region near close packing in order to elucidate the possible existence of disordered phases at close packing.