Molecular dynamics study of the isotropic-nematic quench

Long term phase separation dynamics in liquid crystal-enriched microdroplets obtained from binary fluid mixtures

The dynamics of long term phase separation in binary liquid mixtures remains a subject of fundamental interest. Here, we study a binary liquid mixture, where the minority phase is confined to a liquid crystal (LC)-rich droplet, by investigating the evolution of size, defect and mesogen alignment over time. We track the binary liquid mixture evolving towards equilibrium by visualising the configuration of the liquid crystal droplet through polarisation microscopy. We compare our experimental findings with computational simulations and elucidate differences between bulk phases and confined droplets based on the respective thermodynamics of phase separation. Our work provides insights on how phase transitions on the microscale can deviate from bulk phase diagrams with relevance to other material systems, such as the liquid-liquid phase separation of polymer and protein solutions.

History-dependent phase transition character

We consider history-dependent behavior in domain-type configurations in orientational order that are formed in configurations reached via continuous symmetry-breaking phase transitions. In equilibrium, these systems exhibit in absence of impurities a spatially homogeneous order. We focus on cases where domains are formed via (i) Kibble-Zurek mechanism in fast enough quenches or by (ii) Kibble mechanism in strongly supercooled phases. In both cases, domains could be arrested due to pinned topological defects that are formed at domain walls. In systems exhibiting polar or quadrupolar order, point and line defects (disclinations) dominate, respectively. In particular, the disclinations could form complex entangled structures and are more efficient in stabilizing domains. Domain patterns formed by fast quenches could be arrested by impurities imposing a strong enough random-field type disorder, as suggested by the Imry-Ma theorem. On the other hand, domains formed in supercooled systems could be also formed if large enough energy barriers arresting domains are established due to large enough systems' stiffness. The resulting effective interactions in established domain-type patterns could be described by random matrices. The resulting eigenvectors reveal expected structural excitations formed in such structures. The most important role is commonly played by the random matrix largest eigenvector. Qualitatively different behavior is expected if this eigenvector exhibits a localized or extended character. In the former case, one expects a gradual, non-critical-type transition into a glass-type structure. However, in the latter case, a critical-like phase behavior could be observed.

Structural properties and ring defect formation in discotic liquid crystal nanodroplets

In this work, we performedNpTMonte Carlo simulations of a Gay-Berne discotic liquid crystal confined in a spherical droplet under face-on anchoring and fixed pressure. We find that, in contrast to the unbounded system, a plot of the order parameter as function of temperature does not show a clear evidence of a first-order isotropic-nematic transition. We also find that the impossibility of simultaneously satisfy the uniform director field requirement of a nematic phase with the radial boundary conditions, results in the appearance of a ring disclination line as a stress release mechanism in the interior of the droplet. Under further cooling, a columnar phase appears at the center of the droplet.

Stochastic kinetics reveal imperative role of anisotropic interfacial tension to determine morphology and evolution of nucleated droplets in nematogenic films

For isotropic fluids, classical nucleation theory predicts the nucleation rate, barrier height and critical droplet size by ac- counting for the competition between bulk energy and interfacial tension. The nucleation process in liquid crystals is less understood. We numerically investigate nucleation in monolayered nematogenic films using a mesoscopic framework, in par- ticular, we study the morphology and kinetic pathway in spontaneous formation and growth of droplets of the stable phase in the metastable background. The parameter κ that quantifies the anisotropic elastic energy plays a central role in determining the geometric structure of the droplets. Noncircular nematic droplets with homogeneous director orientation are nucleated in a background of supercooled isotropic phase for small κ. For large κ, noncircular droplets with integer topological charge, accompanied by a biaxial ring at the outer surface, are nucleated. The isotropic droplet shape in a superheated nematic background is found to depend on κ in a similar way. Identical growth laws are found in the two cases, although an unusual two-stage mechanism is observed in the nucleation of isotropic droplets. Temporal distributions of successive events indi- cate the relevance of long-ranged elasticity-mediated interactions within the isotropic domains. Implications for a theoretical description of nucleation in anisotropic fluids are discussed.

An electric-field induced dynamical state in dispersions of highly charged colloidal rods: comparison of experiment and theory

Concentrated dispersions of highly charged rod-like colloids (fd-virus particles) in isotropic-nematic coexistence exhibit a dynamical state when subjected to low-frequency electric fields [Soft Matter, 2010, 6, 273]. This dynamical state consists of nematic domains which persistently melt and form on time scales typically of the order of seconds. The origin of the dynamical state has been attributed to a field-induced, cyclic dissociation and association of condensed ions [Soft Matter, 2014, 10, 1987, Soft Matter, 2015, 11, 2893]. The ionic strength increases on dissociation of condensed ions, rendering the nematic domains unstable, while the subsequent decrease of the ionic strength due to association of condensed ions leads to a recurrent stabilization of the nematic state. The role of dissociation/association of condensed ions in the phase/state behaviour of charged colloids in electric fields has not been addressed before. The electric field strength that is necessary to dissociate sufficient condensed ions to render a nematic domain unstable, depends critically on the ambient ionic strength of the dispersion without the external field, as well as the rod-concentration. The aim of this paper is to compare experimental results for the location of transition lines and the dynamics of melting and forming of nematic domains at various ionic strengths and rod-concentrations with the ion-dissociation/association model. Phase/state diagrams in the field-amplitude versus frequency plane at two different ambient ionic strengths and various rod-concentrations are presented, and compared to the theory. The time scale on which melting and forming of the nematic domains occurs diverges on approach of the transition line where the dynamical state appears. The corresponding critical exponents have been measured by means of image time-correlation spectroscopy [Eur. Phys. J. E, 2009, 30, 333], and are compared to the theoretical values predicted by the ion-dissociation/association model.

An electric-field induced dynamical state in dispersions of charged colloidal rods

The response of concentrated dispersions of charged colloids to low-frequency electric fields is governed by field-induced inter-colloidal interactions resulting from the polarization of electric double layers and the layer of condensed ions, association and dissociation of condensed ions, as well as hydrodynamic interactions through field-induced electro-osmotic flow. The phases and states that can be formed by such field-induced interactions are an essentially unexplored field of research. Experiments on concentrated suspensions of rod-like colloids (fd-virus particles), within the isotropic-nematic phase coexistence region, showed that a number of phases/states are induced, depending on the field amplitude and frequency [Soft Matter, 2010, 6, 273]. In particular, a dynamical state is found where nematic domains form and melt on a time scale of the order of seconds. We discuss the microscopic origin of this dynamical state, which is attributed to the cyclic, electric-field induced dissociation and association of condensed ions. A semi-quantitative theory is presented for the dynamics of melting and formation of nematic domains, including a model for the field-induced dissociation/association of condensed ions. The resulting equation of motion for the orientational order parameter is solved numerically for parameters complying with the fd-virus system. A limit-cycle is found, with a cycling-time that diverges at the transition line in the field-amplitude versus frequency plane where the dynamical state first appears, in accord with experimental findings.

Symmetry breaking in nematic liquid crystals: analogy with cosmology and magnetism

Universal behavior related to continuous symmetry breaking in nematic liquid crystals is studied using Brownian molecular dynamics. A three-dimensional lattice system of rod-like objects interacting via the Lebwohl-Lasher interaction is considered. We test the applicability of predictions originally derived in cosmology and magnetism. In the first part we focus on coarsening dynamics following the temperature driven isotropic-nematic phase transition for different quench rates. The behavior in the early coarsening regime supports predictions made originally by Kibble in cosmology. For fast enough quenches, symmetry breaking and causality give rise to a dense tangle of defects. When the degree of orientational ordering is large enough, well defined protodomains characterized by a single average domain length are formed. With time subcritical domains gradually vanish and supercritical domains grow with time, exhibiting a universal scaling law. In the second part of the paper we study the impact of random-field-type disorder on a range of ordering in the (symmetry broken) nematic phase. We demonstrate that short-range order is observed even for a minute concentration of impurities, giving rise to disorder in line with the Imry-Ma theorem prediction only for the appropriate history of systems.

Symmetry breaking and structure of a mixture of nematic liquid crystals and anisotropic nanoparticles

Orientational ordering of a homogeneous mixture of uniaxial liquid crystalline (LC) molecules and magnetic nanoparticles (NPs) is studied using the Lebwohl-Lasher lattice model. We consider cases where NPs tend to be oriented perpendicularly to LC molecules due to elastic forces. We study domain-type configurations of ensembles, which are quenched from the isotropic phase. We show that for large enough concentrations of NPs the long range uniaxial nematic ordering is replaced by short range order exhibiting strong biaxiality. This suggests that the impact of NPs on orientational ordering of LCs for appropriate concentrations of NPs is reminiscent to the influence of quenched random fields which locally enforce a biaxial ordering.

Orientational quenched disorder of a nematic liquid crystal

By means of direct imaging, we map the surface heterogeneities of the nematic director orientation on a SiOx anchoring layer. The spatial correlations of surface director orientations are well fitted with a compressed exponential with exponent of 1.5 and typical correlation length of few microns. To discuss these results a formal analogy is established between the equation governing the nematic surface torques and the Langevin equation. Based on this analogy we prove that the disorder is spatially correlated orientational quenched disorder. The measured correlation length is discussed in terms of substrate morphology and molecular adsorption.

Thermotropic uniaxial nematic liquid crystal with ± 1/2 strength disclinations: defect-antidefect production and correlation

We report the observation of mostly wedge disclinations of strength ± 1/2 in a uniaxial nematic liquid crystal. The dense network of the disclination is stabilized in cells of thickness (∼6 µm) treated with perfluoropolymer. The distribution of net defect strength was found to be Gaussian. The scaling exponent obtained from the width of the distribution is ν = 0.29 ± 0.08 and is in good agreement with the prediction based on the cosmological Kibble mechanism. The exponent and the constant of the scaling relation are compared with the values reported for lyotropic liquid crystal with half-strength disclination by Mukai and co-workers.

The influence of disorder on thermotropic nematic liquid crystals phase behavior

We review the theoretical research on the influence of disorder on structure and phase behavior of condensed matter system exhibiting continuous symmetry breaking focusing on liquid crystal phase transitions. We discuss the main properties of liquid crystals as adequate systems in which several open questions with respect to the impact of disorder on universal phase and structural behavior could be explored. Main advantages of liquid crystalline materials and different experimental realizations of random field-type disorder imposed on liquid crystal phases are described.

Isotropic-to-nematic phase transition in a liquid-crystal droplet

In this paper, we focus on the isotropic-to-nematic phase transition in a liquid-crystal droplet. We present the results of an experiment to measure the growth of the nematic phase within an isotropic phase liquid-crystal droplet. Experimentally, we observe two primary phase transition regimes. At short time scales, our experimental results (R(t) approximately t0.51) show good agreement with a Stefan-type model of the evolution of the nematic phase within the isotropic phase of a liquid crystal. As time progresses, the growth of the nematic phase is restricted by increased confinement of the droplet boundary. During this stage of growth, the nematic phase grows at a slower rate of R(t) approximately t0.31. The slower growth at later stages might be due to a variety of factors such as confinement-induced latent heat reduction; a change of defect strength during its evolution; or interactions between the defect and the interface between the liquid crystal and oil or between adjacent defects. The presence of two growth regimes is also consistent with the molecular simulations of Bradac et al. (Bradac, Z.; Kralj, S.; Zumer, S. Phys. Rev. E 2002, 65, 021705) who identify an early stage domain regime and a late stage confinement regime. For the domain and confinement regimes, Bradac et al. (Bradac, Z.; Kralj, S.; Zumer, S. Phys. Rev. E 2002, 65, 021705) obtained growth exponents of 0.49 +/- 0.05 and 0.25 +/- 0.05. These are remarkably close to the values 0.51 and 0.31 observed in our experiments.

Defect-enhanced nematic surface order reconstruction

Within the Landau-de Gennes phenomenological theory, we study the influence of an applied electric field with average strength E{a} on the position of a nematic line defect with topological charge M=+/-1/2 in a hybrid cell. We explore the biaxial structure of the defect core and we describe its expulsion from the cell upon increasing E{a}. We show that prior to the expulsion the defect core displays dramatic changes for strong enough surface anchorings. At a critical value of E{a}, the core broadens and merges into a surface layer with a large biaxiality. This transition corresponds to the reconstruction of the nematic order already observed in the bulk in response to an applied electric field. A similar order reconstruction could take place even in the absence of defects, but at a higher threshold.

Annihilation of nematic point defects: pre-collision and post-collision evolution

The annihilation of the nematic hedgehog and anti-hedgehog within an infinite cylinder of radius R is studied. The semi-microscopic lattice-type model and Brownian molecular dynamics are used. We distinguish among the i) early pre-collision, ii) late pre-collision, iii) early post-collision, and iv) late post-collision stages. In the pre-collision stage our results agree qualitatively with the existing experimental observations and also continuum-type simulations. The core of each defect exhibits a ring-like structure, where the ring axis is set perpendicular to the cylinder symmetry axis. For xi(0)d/(2R) > 1 the interaction between defects is negligible, where xi(0)d describes the initial separation of defects. Consequently, the defects annihilate within the simulation time window for xi(0)d/(2R) < 1. For close enough defects their separation scales as xi(d) [see text] (t(c)- t)(0.4+/-0.1), where t(c) stands for the collision time. In elastically anisotropic medium the hedgehog is faster than the anti-hedgehog. In the early pre-collision stage the defects can be treated as point-like particles, possessing inherent core structure, that interact via the nematic director field. In the late pre-collision stage the cores reflect the interaction between defects. After the collision a charge-less ring structure is first formed. In the early post-collision stage the ring adopts an essentially untwisted circular structure of the radius xi(r). In the late post-collision stage we observe two qualitatively different scenarios. For mu = xi(r)/R < mu(c) approximately 0.25 the ring collapses leading to the escaped radial equilibrium structure. For mu > mu(c) the chargeless ring triggers the nucleation growth into the planar polar structure with line defects.

Random anisotropy nematic model: nematic-non-nematic mixture

The influence of a random-anisotropy- (RA-) type disorder on the phase separation of the nematogen-non-nematogen mixture is studied. A combination of the phenomenological Landau-de Gennes and Flory-Huggins theories is used. We assume that the non-nematogen component (i.e., impurity) enforces the RA disorder to the enclosing thermotropic liquid-crystal (LC) phase. The Imry-Ma argument is used according to which the lower-temperature phase exhibits a domain-type pattern. The disorder strength is measured in terms of the dimensionless parameter Lambda. We consider the case in which the LC molecules and impurities mix in the isotropic phase for Lambda=0. The impurities enforce a finite degree of orientational ordering even in the high-temperature paranematic phase. In the low-temperature phase they give rise to a domain-type structure, resulting in the distorted nematic (speronematic) phase. We show that the onset of orientational ordering increases the phase separation tendency. The RA field, however, opposes this tendency. With increasing value of Lambda the difference between the paranematic and speronematic ordering decreases. Consequently the structure of the phase-separated pattern can be much more complex in comparison to the Lambda=0 case.

Annihilation of nematic point defects: postcollision scenarios

We perform a study of the annihilation of a nematic radial and hyperbolic point defects with the main focus on the confinement induced collision and postcollision scenarios. Brownian molecular dynamics on a semimicroscopic lattice is used. Initially a pair of defects, separated for 1.4-1.7 radii, is induced at the axis of the cylindrical capillary. In such a configuration defects start to approach slowly. In the early stage, their cores are negligibly influenced by the mutual interaction. When the distance becomes comparable to the nematic correlation length, the cores significantly deform. In the collision regime, defects gradually merge. We observe two qualitatively different scenarios in the postcollision regime, depending on the degree of (meta) stability of the initially imposed escaped structure with point defects.