Topological defects in a two-dimensional liquid crystal confined in a circular nanocavity

Nematic-isotropic phase transitions in thin slabs of liquid crystals with topological defect arrays

We study the nematic-to-isotropic phase transitions in thin slabs of nematic liquid crystals with photopatterned director fields of topological defect arrays at constant heating rates and show that the transition kinetics is significantly impacted by both the heating rate and the topological strengths of these defects. Specifically, with ±1/2 defect arrays, the isotropic domains emerge from the defect cores when the heating rate is high, while from random places when the heating rate is low. With ±1 defect arrays, the isotropic domains always emerge from the defect cores regardless of the heating rate. Furthermore, the isotropic domains show significant movements at slow heating rates, and the total area of the isotropic domains grows with the temperature T following a simple power law (T - T')γ, where the exponent γ is approximately 1 in most cases and is 2/3 for the ±1 defect arrays at low heating rates when the isotropic domains are pinned on the defect cores. We attribute this phenomenon to an interplay between the surface tension and bulk free energy.

Theoretical Understanding of Polar Topological Phase Transitions in Functional Oxide Heterostructures: A Review

The exotic topological phase is attracting considerable attention in condensed matter physics and materials science over the past few decades due to intriguing physical insights. As a combination of "topology" and "ferroelectricity," the ferroelectric (polar) topological structures are a fertile playground for emergent phenomena and functionalities with various potential applications. Herein, the review starts with the universal concept of the polar topological phase and goes on to briefly discuss the important role of computational tools such as phase-field simulations in designing polar topological phases in oxide heterostructures. In particular, the history of the development of phase-field simulations for ferroelectric oxide heterostructures is highlighted. Then, the current research progress of polar topological phases and their emergent phenomena in ferroelectric functional oxide heterostructures is reviewed from a theoretical perspective, including the topological polar structures, the establishment of phase diagrams, their switching kinetics and interconnections, phonon dynamics, and various macroscopic properties. Finally, this review offers a perspective on the future directions for the discovery of novel topological phases in other ferroelectric systems and device design for next-generation electronic device applications.

Topological fine structure of smectic grain boundaries and tetratic disclination lines within three-dimensional smectic liquid crystals

Observing and characterizing the complex ordering phenomena of liquid crystals subjected to external constraints constitutes an ongoing challenge for chemists and physicists alike. To elucidate the delicate balance appearing when the intrinsic positional order of smectic liquid crystals comes into play, we perform Monte-Carlo simulations of rod-like particles in a range of cavities with a cylindrical symmetry. Based on recent insights into the topology of smectic orientational grain boundaries in two dimensions, we analyze the emerging three-dimensional defect structures from the perspective of tetratic symmetry. Using an appropriate three-dimensional tetratic order parameter constructed from the Steinhardt order parameters, we show that those grain boundaries can be interpreted as a pair of tetratic disclination lines that are located on the edges of the nematic domain boundary. Thereby, we shed light on the fine structure of grain boundaries in three-dimensional confined smectics.

Defect patterns of two-dimensional nematic liquid crystals in confinement

A two-dimensional or quasi-two-dimensional nematic liquid crystal refers to a surface-confined system. When such a system is further confined by external line boundaries or excluded from internal line boundaries, the nematic directors form a deformed texture that may display defect points or defect lines, for which winding numbers can be clearly defined. Here, a particular attention is paid to the case when the liquid crystal molecules prefer to form a boundary nematic texture in parallel to the wall surface (i.e., following the homogeneous boundary condition). A general theory, based on geometric argument, is presented for the relationship between the sum of all winding numbers in the system (the total winding number) and the type of confinement angles and curved segments. The conclusion is validated by comparing the theoretical defect rule with existing nematic textures observed experimentally and theoretically in recent years.

Confinement effects on the random sequential adsorption packings of elongated particles in a slit

The behavior of a system of two-dimensional elongated particles (discorectangles) packed in a slit between the two parallel walls was analyzed using a simulation approach. The packings were produced using the random sequential adsorption model with continuous positional and orientational degrees of freedom. The aspect ratio (length-to-width ratio, ɛ=l/d) of the particles was varied within the range ɛ∈[1;32] while the distance between the walls was varied within the range h/d∈[1;80]. The properties of deposits in jammed state [the coverage, the order parameter, and the long-range (percolation) connectivity between particles] were studied numerically. The values of ɛ and h significantly affected the structure of the packings and the percolation connectivity. Particularly, the observed nontrivial dependencies of the jamming coverage φ(ɛ) or φ(h) were explained by the interplay of the different geometrical factors related to confinement, particle orientation degrees of freedom and excluded volume effects.

Particle-resolved topological defects of smectic colloidal liquid crystals in extreme confinement

Confined samples of liquid crystals are characterized by a variety of topological defects and can be exposed to external constraints such as extreme confinements with nontrivial topology. Here we explore the intrinsic structure of smectic colloidal layers dictated by the interplay between entropy and an imposed external topology. Considering an annular confinement as a basic example, a plethora of competing states is found with nontrivial defect structures ranging from laminar states to multiple smectic domains and arrays of edge dislocations, which we refer to as Shubnikov states in formal analogy to the characteristic of type-II superconductors. Our particle-resolved results, gained by a combination of real-space microscopy of thermal colloidal rods and fundamental-measure-based density functional theory of hard anisotropic bodies, agree on a quantitative level.

Rodlike molecules in extreme confinement

A unique feature of colloid particles and biopolymers is the molecule's intrinsic rigidity characterized by a molecular-level length scale. Under extreme confinement conditions at cellular scales or in nanodevices, these molecules can display orientational ordering accompanied by severe density depletion. Conventional liquid-crystal theories, such as the Oseen-Frank or Landau-de Gennes theories, cannot capture the essential molecular-level properties: the boundary effects, which extend to a distance of the rigidity length scale, and the drastic variations of the inhomogeneous molecular density. Here we show, based on a simple interpretation of the Onsager model, that rodlike molecules in extreme annular confinement produce unusual liquid-crystal defect structures that are independent phases from the patterns usually seen in a weaker confinement environment.

Highly confined mixtures of parallel hard squares: A density-functional-theory study

Using the fundamental-measure density-functional theory, we study theoretically the phase behavior of extremely confined mixtures of parallel hard squares in slit geometry. The pore width is chosen such that configurations consisting of two consecutive big squares, or three small squares, in the transverse direction, perpendicular to the walls, are forbidden. We analyze two different mixtures with edge lengths of species selected so as to allow or forbid one big plus one small square to fit into the channel. For the first mixture we obtain first-order transitions between symmetric and asymmetric packings of particles: Small and big squares are preferentially adsorbed at different walls. Asymmetric configurations are shown to lead to more efficient packing at finite pressures. We argue that the stability region of the asymmetric phase in the pressure-composition plane is bounded so that the symmetric phase is stable at low and very high pressure. For the second mixture, we observe strong demixing between phases which are rich in different species. Demixing occurs in the lateral direction, i.e., the dividing interface is perpendicular to the walls, and phases exhibit symmetric density profiles. The possible experimental realization of this behavior (which in practical terms is precluded by jamming) in strictly two-dimensional systems is discussed. Finally, the phase behavior of a mixture with periodic boundary conditions is analyzed and the differences and similarities between the latter and the confined system are discussed. We claim that, although exact calculations exclude the existence of true phase transitions in (1+ε)-dimensional systems, density-functional theory is still successful in describing packing properties of large clusters of particles.

Topological defects in two-dimensional liquid crystals confined by a box

When a spatially uniform system that displays a liquid-crystal ordering on a two-dimensional surface is confined inside a rectangular box, the liquid crystal direction field develops inhomogeneous textures accompanied by topological defects because of the geometric frustrations. We show that the rich variety of nematic textures and defect patterns found in recent experimental and theoretical studies can be classified by the solutions of the rather fundamental, extended Onsager model. This is critically examined based on the determined free energies of different defect states, as functions of a few relevant, dimensionless geometric parameters.

Finite particle size drives defect-mediated domain structures in strongly confined colloidal liquid crystals

When liquid crystals are confined to finite volumes, the competition between the surface anchoring imposed by the boundaries and the intrinsic orientational symmetry-breaking of these materials gives rise to a host of intriguing phenomena involving topological defect structures. For synthetic molecular mesogens, like the ones used in liquid-crystal displays, these defect structures are independent of the size of the molecules and well described by continuum theories. In contrast, colloidal systems such as carbon nanotubes and biopolymers have micron-sized lengths, so continuum descriptions are expected to break down under strong confinement conditions. Here, we show, by a combination of computer simulations and experiments with virus particles in tailor-made disk- and annulus-shaped microchambers, that strong confinement of colloidal liquid crystals leads to novel defect-stabilized symmetrical domain structures. These finite-size effects point to a potential for designing optically active microstructures, exploiting the as yet unexplored regime of highly confined liquid crystals.

Liquid crystals confined to micrometre-sized geometries can be well described by a continuum theory, where the size effect of constituent mesogens is negligible. Here, the authors show how the continuum theory breaks down in colloidal liquid crystal, leading to the formation of defect-mediated domains.

Nematic ordering of semiflexible polymers confined on a toroidal surface

We study the isotropic-like and nematic states of wormlike liquid-crystal polymers embedded on the surface of a torus. The role played by surface curvature, which couples to the molecular rigidity, is reported as the main reason that causes the weak nematic ordering in an otherwise ordinary isotropic phase. The same coupling has a profound effect on the nematic states as well, which are stabilized by the Onsager excluded-volume interaction; the latter has been frequently used to study lyotropic liquid crystal polymers and is used here as an example of the physical mechanisms that drive the system to make orientational ordering. We identify important parameters in the system which are used as axes of the four-dimensional phase diagram. The numerical study demonstrates a strong correlation between the liquid-crystal defect-free and defect structures and the geometry of the liquid-crystal embedded surface.

Ordering in granular-rod monolayers driven far from thermodynamic equilibrium

The orientational order in vertically agitated granular-rod monolayers is investigated experimentally and compared quantitatively with equilibrium Monte Carlo simulations and density functional theory. At sufficiently high number density, short rods form a tetratic state and long rods form a uniaxial nematic state. The length-to-width ratio at which the order changes from tetratic to uniaxial is around 7.3 in both experiments and simulations. This agreement illustrates the universal aspects of the ordering of rod-shaped particles across equilibrium and nonequilibrium systems. Moreover, the assembly of granular rods into ordered states is found to be independent of the agitation frequency and strength, suggesting that the detailed nature of energy injection into such a nonequilibrium system does not play a crucial role.

Defect structures mediate the isotropic-nematic transition in strongly confined liquid crystals

Using Monte Carlo simulations, we study rod-like lyotropic liquid crystals confined to a square slab-like geometry with lateral dimensions comparable to the length of the particles. We observe that this system develops linear defect structures upon entering the planar nematic phase. These defect structures flank a lens-shaped nematic region oriented along a diagonal of the square box. We interpret these structures as a compromise between the 2-fold order of the bulk nematic phase and the 4-fold order imposed by the lateral boundaries. A simple Onsager-type theory that effectively implements these competing tendencies is used to model the phase behavior in the center of the box and shows that the free-energy cost of forming the defect structures strongly offsets the transition-inducing effects of both the transverse and lateral confinement.

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.

Texture and shape of two-dimensional domains of nematic liquid crystals

We present a generalized approach to compute the shape and internal structure of two-dimensional nematic domains. By using conformal mappings, we are able to compute the director field for a given domain shape that we choose from a rich class, which includes drops with large and small aspect ratios and sharp domain tips as well as smooth ones. Results are assembled in a phase diagram that for given domain size, surface tension, anchoring strength, and elastic constant shows the transitions from a homogeneous to a bipolar director field, from circular to elongated droplets, and from sharp to smooth domain tips. We find a previously unaccounted for regime, where the drop is nearly circular, the director field bipolar, and the tip rounded. We also find that bicircular director fields, with foci that lie outside the domain, provide a remarkably accurate description of the optimal director field for a large range of values of the various shape parameters.

Spontaneous circulation of confined active suspensions

Many active fluid systems encountered in biology are set in total geometric confinement. Cytoplasmic streaming in plant cells is a prominent and ubiquitous example, in which cargo-carrying molecular motors move along polymer filaments and generate coherent cell-scale flow. When filaments are not fixed to the cell periphery, a situation found both in vivo and in vitro, we observe that the basic dynamics of streaming are closely related to those of a nonmotile stresslet suspension. Under this model, it is demonstrated that confinement makes possible a stable circulating state; a linear stability analysis reveals an activity threshold for spontaneous autocirculation. Numerical analysis of the longtime behavior reveals a phenomenon akin to defect separation in nematic liquid crystals and a high-activity bifurcation to an oscillatory regime.

Nematic order in small systems: measuring the elastic and wall-anchoring constants in vibrofluidized granular rods

We investigate nematic order in vibrated granular rods confined to a small quasi-2D container less than 10 rod lengths in diameter. As rod density ρ increases, patterning shifts from bipolar to uniform alignment. We find that a continuum liquid crystal free energy functional captures key patterning features down to almost the particle size. By fitting theory to experiments, we estimate the relative values of bend and splay elastic constants and wall anchoring. We find that splay is softer than bend for all ρ and rod lengths tested, while the ratio of the average elastic constant to wall anchoring increases with ρ.

Nematic liquid crystal dynamics under applied electric fields

In this paper we investigate the coarsening dynamics of liquid crystal textures in a two-dimensional nematic under applied electric fields, using numerical simulations performed using a publicly available liquid crystal algorithm developed by the authors. We consider both positive and negative dielectric anisotropies and two different possibilities for the orientation of the electric field (parallel and perpendicular to the two-dimensional lattice). We determine the effect of an applied electric field pulse on the evolution of the characteristic length scale and other properties of the liquid crystal texture network. In particular, we show that different types of defects are produced after the electric field is switched on, depending on the orientation of the electric field and the sign of the dielectric anisotropy.