Discontinuous structural transition in a thin hybrid liquid crystal film

Ordering of Oblate Hard Particles between Hybrid Penetrable Walls

We report a Monte Carlo (MC) simulation study of a model discotic liquid crystal (DLC) confined between hybrid walls with controllable penetrability. The model consists of oblate hard Gaussian overlap (HGO) particles. Particle-substrate interactions are modeled as follows: each substrate sees a particle as a disc of zero thickness and diameter D less than or equal to that of the actual particle, σ0, embedded inside the particle and located halfway along, and perpendicular to, its minor axis. This allows us to control the anchoring properties of the substrates, from planar (edge-on) for D ≈ 0 to homeotropic (face-on) for D ≈ σ0, which can be done independently at either substrate. Depending on the values of Ds ≡ D/σ0 at the top (Dst) and bottom (Dsb) substrates, we find domains in (Dsb, Dst) space in which particle alignment is uniform planar (UP), is uniform homeotropic (UH), or varies linearly from planar at one substrate to homeotropic at the other (Lin). These domains are separated by regions of bistability (P-Lin and H-Lin), which appear to be wider than for prolate HGOs, and there may be also a small tristable (P-H-Lin) region. Results are compared with the predictions of density functional theory, implemented at the level of Onsager's second-virial approximation with Parsons-Lee rescaling. As in the case of symmetric confinement studied previously, the agreement between theory and simulation is substantially less good than for prolate HGOs: in particular, for the investigated substrate separation L = 6σ0, the Lin configuration is never predicted. These discrepancies are likely a consequence of the fact that Onsager's theory is less accurate for discs than for rods.

Ordering transitions of weakly anisotropic hard rods in narrow slitlike pores

The effect of strong confinement on the positional and orientational ordering is examined in a system of hard rectangular rods with length L and diameter D (L>D) using the Parsons-Lee modification of the second virial density-functional theory. The rods are nonmesogenic (L/D<3) and confined between two parallel hard walls, where the width of the pore (H) is chosen in such a way that both planar (particle's long axis parallel to the walls) and homeotropic (particle's long axis perpendicular to the walls) orderings are possible and a maximum of two layers is allowed to form in the pore. In the extreme confinement limit of H≤2D, where only one-layer structures appear, we observe a structural transition from a planar to a homeotropic fluid layer with increasing density, which becomes sharper as L→H. In wider pores (2D<H<3D) planar order with two layers, homeotropic order, and even combined bilayer structures (one layer is homeotropic, while the other is planar) can be stabilized at high densities. Moreover, first-order phase transitions can be seen between different structures. One of them emerges between a monolayer and a bilayer with planar orders at relatively low packing fractions.

Computational aspects of the smectization process in liquid crystals: An example study of a perfectly aligned two-dimensional hard-boomerang system

A replica method for calculation of smectic liquid crystal properties within the Onsager theory has been presented and applied to an exemplary case of two-dimensional perfectly aligned needlelike boomerangs. The method allows one to consider the complete influence of the interaction terms in contrast to the Fourier expansion method which uses mostly first or second order terms of expansion. The program based on the replica algorithm is able to calculate a single representative layer as an equivalent set of layers, depending on the size of the considered width of the sample integration interval. It predicts successfully smectic density distributions, energies, and layer thicknesses for different types of layer arrangement-of the antiferroelectric or of the smectic A order type. Specific features of the algorithm performance and influence of the numerical accuracy on the physical properties are presented. Future applications of the replica method to freely rotating molecules are discussed.

Competition of lattice and basis for alignment of nematic liquid crystals

Due to elastic anisotropy, two-dimensional patterning of substrates can promote weak azimuthal alignment of adjacent nematic liquid crystals. Here we consider how such alignment can be achieved using a periodic square lattice of circular or elliptical motifs. In particular, we examine ways in which the lattice and motif can combine to favor differing orientations. Using Monte Carlo simulation and continuum elasticity we find, for circular motifs, that the coverage fraction controls both the polar anchoring angle and a transition in the azimuthal orientation. If the circles are generalized to ellipses, arbitrary control of the effective easy axis and effective anchoring potential becomes achievable by appropriate tuning of the ellipse motif relative to the periodic lattice patterning. This has possible applications in both monostable and bistable liquid crystal device contexts.

Neural-network approach to modeling liquid crystals in complex confinement

Finding the structure of a confined liquid crystal is a difficult task since both the density and order parameter profiles are nonuniform. Starting from a microscopic model and density-functional theory, one has to either (i) solve a nonlinear, integral Euler-Lagrange equation, or (ii) perform a direct multidimensional free energy minimization. The traditional implementations of both approaches are computationally expensive and plagued with convergence problems. Here, as an alternative, we introduce an unsupervised variant of the multilayer perceptron (MLP) artificial neural network for minimizing the free energy of a fluid of hard nonspherical particles confined between planar substrates of variable penetrability. We then test our algorithm by comparing its results for the structure (density-orientation profiles) and equilibrium free energy with those obtained by standard iterative solution of the Euler-Lagrange equations and with Monte Carlo simulation results. Very good agreement is found and the MLP method proves competitively fast, flexible, and refinable. Furthermore, it can be readily generalized to the richer experimental patterned-substrate geometries that are now experimentally realizable but very problematic to conventional theoretical treatments.

Independent control of polar and azimuthal anchoring

Monte Carlo simulation, experiment, and continuum theory are used to examine the anchoring exhibited by a nematic liquid crystal at a patterned substrate comprising a periodic array of rectangles that, respectively, promote vertical and planar alignment. It is shown that the easy axis and effective anchoring energy promoted by such surfaces can be readily controlled by adjusting the design of the pattern. The calculations reveal rich behavior: for strong anchoring, as exhibited by the simulated system, for rectangle ratios ≥2 the nematic aligns in the direction of the long edge of the rectangles, the azimuthal anchoring coefficient changing with pattern shape. In weak anchoring scenarios, however, including our experimental systems, preferential anchoring is degenerate between the two rectangle diagonals. Bistability between diagonally aligned and edge-aligned arrangement is predicted for intermediate combinations of anchoring coefficient and system length scale.

Competing alignments of nematic liquid crystals on square-patterned substrates

A theoretical analysis is presented of a nematic liquid crystal confined between substrates patterned with squares that promote vertical and planar alignment. Two approaches are used to elucidate the behavior across a wide range of length scales: Monte Carlo simulation of hard particles and Frank-Oseen continuum theory. Both approaches predict bistable degenerate azimuthal alignment in the bulk along the edges of the squares; the continuum calculation additionally reveals the possibility of an anchoring transition to diagonal alignment if the polar anchoring energy associated with the pattern is sufficiently weak. Unlike the striped systems previously analyzed, the Monte Carlo simulations suggest that there is no "bridging" transition for sufficiently thin cells. The extent to which these geometrically patterned systems resemble topographically patterned substrates, such as square wells, is also discussed.

Density functional theory of liquid crystals and surface anchoring: hard Gaussian overlap-sphere and hard Gaussian overlap-surface potentials

This article applies the density functional theory to confined liquid crystals, comprised of ellipsoidal shaped particles interacting through the hard Gaussian overlap (HGO) potential. The extended restricted orientation model proposed by Moradi and co-workers [J. Phys.: Condens. Matter 17, 5625 (2005)] is used to study the surface anchoring. The excess free energy is calculated as a functional expansion of density around a reference homogeneous fluid. The pair direct correlation function (DCF) of a homogeneous HGO fluid is approximated, based on the optimized sum of Percus-Yevick and Roth DCF for hard spheres; the anisotropy introduced by means of the closest approach parameter, the expression proposed by Marko [Physica B 392, 242 (2007)] for DCF of HGO, and hard ellipsoids were used. In this study we extend an our previous work [Phys. Rev. E 72, 061706 (2005)] on the anchoring behavior of hard particle liquid crystal model, by studying the effect of changing the particle-substrate contact function instead of hard needle-wall potentials. We use the two particle-surface potentials: the HGO-sphere and the HGO-surface potentials. The average number density and order parameter profiles of a confined HGO fluid are obtained using the two particle-wall potentials. For bulk isotropic liquid, the results are in agreement with the Monte Carlo simulation of Barmes and Cleaver [Phys. Rev. E 71, 021705 (2005)]. Also, for the bulk nematic phase, the theory gives the correct density profile and order parameter between the walls.

Nematic liquid-crystal alignment on stripe-patterned substrates

Here, we use molecular simulation to consider the behavior of a thin nematic film confined between two identical nanopatterned substrates. Using patterns involving alternating stripes of homeotropic-favoring and homogeneous-favoring substrates, we investigate the influence of the relative stripe width and the film thickness. From this, we show that the polar anchoring angle can be varied continuously from planar to homeotropic by appropriate tuning of these parameters. For very thin films with equal stripe widths, we observe orientational bridging, the surface patterning being written in domains which traverse the nematic film. This dual-bridging-domain arrangement breaks down with increase in film thickness, however, being replaced by a single tilted monodomain. Strong azimuthal anchoring in the plane of the stripe boundaries is observed for all systems.

Competition between capillarity, layering and biaxiality in a confined liquid crystal

The effect of confinement on the phase behaviour and structure of fluids made of biaxial hard particles (cuboids) is examined theoretically by means of Onsager second-order virial theory in the limit where the long particle axes are frozen in a mutually parallel configuration. Confinement is induced by two parallel planar hard walls (slit-pore geometry), with particle long axes perpendicular to the walls (perfect homeotropic anchoring). In bulk, a continuous nematic-to-smectic transition takes place, while shape anisotropy in the (rectangular) particle cross-section induces biaxial ordering. As a consequence, four bulk phases, uniaxial and biaxial nematic and smectic phases, can be stabilised as the cross-sectional aspect ratio is varied. On confining the fluid, the nematic-to-smectic transition is suppressed, and either uniaxial or biaxial phases, separated by a continuous transition, can be present. Smectic ordering develops continuously from the walls for increasing particle concentration (in agreement with the supression of nematic-smectic second-order transition at confinement), but first-order layering transitions, involving structures with n and n + 1 layers, arise in the confined fluid at high concentration. Competition between layering and uniaxial-biaxial ordering leads to three different types of layering transitions, at which the two coexisting structures can be both uniaxial, one uniaxial and another biaxial, or both biaxial. Also, the interplay between molecular biaxiality and wall interactions is very subtle: while the hard wall disfavours the formation of the biaxial phase, biaxiality is against the layering transitions, as we have shown by comparing the confined phase behaviour of cylinders and cuboids. The predictive power of Onsager theory is checked and confirmed by performing some calculations based on fundamental-measure theory.

Simulation and theory of hybrid aligned liquid crystal films

We present a study of the effects of nanoconfinement on a system of hard Gaussian overlap particles interacting with planar substrates through the hard-needle-wall potential, extending earlier work by two of us [D. J. Cleaver and P. I. C. Teixeira, Chem. Phys. Lett. 338, 1 (2001)]. Here, we consider the case of hybrid films, where one of the substrates induces strongly homeotropic anchoring, while the other favors either weakly homeotropic or planar anchoring. These systems are investigated using both Monte Carlo simulation and density-functional theory, the latter implemented at the level of Onsager's second-virial approximation with Parsons-Lee rescaling. The orientational structure is found to change either continuously or discontinuously depending on substrate separation, in agreement with earlier predictions by others. The theory is seen to perform well in spite of its simplicity, predicting the positional and orientational structure seen in simulations even for small particle elongations.

Structure and stability of isotropic states of hard platelet fluids

We study the thermodynamics and the pair structure of hard, infinitely thin, circular platelets in the isotropic phase. Monte Carlo simulation results indicate a rich spatial structure of the spherical expansion components of the direct correlation function, including nonmonotonical variation of some of the components with density. Integral equation theory is shown to reproduce the main features observed in simulations. The hypernetted chain closure, as well as its extended versions that include the bridge function up to second and third order in density, perform better than both the Percus-Yevick closure and Verlet bridge function approximation. Using a recent fundamental measure density functional theory, an analytic expression for the direct correlation function is obtained as the sum of the Mayer bond and a term proportional to the density and the intersection length of two platelets. This is shown to give a reasonable estimate of the structure found in simulations, but to fail to capture the nonmonotonic variation with density. We also carry out a density functional stability analysis of the isotropic phase with respect to nematic ordering and show that the limiting density is consistent with that where the Kerr coefficient vanishes. As a reference system, we compare to simulation results for hard oblate spheroids with small, but nonzero elongations, demonstrating that the case of vanishingly thin platelets is approached smoothly.

Control of the nematic-isotropic phase transition by an electric field

We use a relatively simple continuum model to investigate the effects of dielectric inhomogeneity within confined liquid-crystal cells. Specifically, we consider, in planar, cylindrical, and spherical geometries, the stability of a nematic-isotropic interface subject to an applied voltage when the nematic liquid crystal has a positive dielectric anisotropy. Depending on the magnitude of this voltage, the temperature, and the geometry of the cell, the nematic region may shrink until the material is completely isotropic within the cell, grow until the nematic phase fills the cell, or, in certain geometries, coexist with the isotropic phase. For planar geometry, no coexistence is found, but we are able to give analytical expressions for the critical voltage for an electric-field-induced phase transition as well as the critical wetting layer thickness for arbitrary applied voltage. In cells with cylindrical and spherical geometries, however, locally stable nematic-isotropic coexistence is predicted, the thickness of the nematic region being controllable by alteration of the applied voltage.

Density functional theory of liquid crystals and surface anchoring

This paper applies the density functional theory to confined liquid crystals comprising ellipsoidal shaped particles interacting through the hard Gaussian overlap (HGO) potential. The restricted orientation model proposed by Rickayzen [Mol. Phys. 95, 393 (1998)] is extended to study the surface anchoring. The excess free energy is calculated as a functional expansion of density around a reference homogeneous fluid. The pair direct correlation function (DCF) of a homogeneous HGO fluid is approximated, based on the Percus-Yevick DCF for hard spheres; the anisotropy is introduced by means of the closest approach parameter. The average number density and orientational order parameter profiles of a HGO fluid confined in between planar walls are obtained using a hard needle-wall potential to represent the particle-wall interactions. For short and long needle lengths, the homeotropic and planar anchoring are observed, respectively. For the bulk isotropic phase the calculated density and order parameter profiles are in agreement with the Monte Carlo simulation of Barmes and Cleaver [Phys. Rev. E 69, 61705 (2004)]. However, for the bulk nematic phase the theory gives the correct density profile between the walls. The correct order parameters are obtained close to the walls whereas for the region in the middle of the walls, the agreement is less satisfactory.

Density profile and order parameter of a hard ellipsoidal fluid confined to a slit

The density profile and order parameter of a fluid of hard axially symmetric ellipsoids confined in between two parallel hard walls is obtained by using the density functional theory. The required input direct correlation function of the homogeneous fluid is calculated by the variational method introduced by Marko (1989 Phys. Rev. 39 2050) and the modified closest approach method proposed by Rickayzen (1998 Mol. Phys. 95 393). Here the restricted orientation model, ROM, is extended to study a fluid comprising molecules which can be aligned in more than six directions, making it more representative of a normal fluid. The density profiles, the average number density and order parameter are obtained for different values of density and elongations. The results are in agreement with the previous theory and available Monte Carlo simulation results.

Using particle shape to induce tilted and bistable liquid crystal anchoring

We use Monte Carlo simulations of hard Gaussian overlap (HGO) particles symmetrically confined in slab geometry to investigate the role of particle-substrate interactions on liquid crystalline anchoring. Despite the restriction here to purely steric interactions and smooth substrates, a range of behaviors are captured, including tilted anchoring and homeotropic-planar bistability. These macroscopic behaviors are all achieved through appropriate tuning of the microscopics of the HGO-substrate interaction, based upon nonadditive descriptions for the HGO-substrate shape parameter.

Computer simulation of a liquid-crystal anchoring transition

We present a study of the effects of confinement on a system of hard Gaussian overlap particles interacting with planar substrates through the hard-needle-wall potential. Using geometrical arguments to calculate the molecular volume absorbed at the substrates, we show that both planar and homeotropic arrangements can be obtained using this model. Monte Carlo simulations are then used to perform a systematic study of the model's behavior as a function of the system density and the hard-needle-wall interaction parameter. As well as showing the homeotropic to planar anchoring transition, the anchoring phase diagrams computed from these simulations indicate regions of bistability. This bistable behavior is examined further through the explicit simulation of field-induced two-way switching between the two arrangements.

Structures and transitions in thin hybrid nematic films: a Monte Carlo study

We confirm by Monte Carlo simulations of a Lebwohl-Lasher lattice spin model the existence of a biaxially ordered nonbent structure in a liquid-crystalline cell subject to opposing boundary conditions. We report on the observation of the bending transition from the biaxial to the bent-director structure when the temperature of the system is lowered. The structural transition is monitored both by the change of the order parameters and by heat capacity. We discuss the thickness dependence of the transition temperature by means of wetting-induced phenomena and elastic deformations. We propose the correspondence to the phenomenological description, which agrees well without any fitting parameters.