Computational modelling of nematic phase ordering by film and droplet growth over heterogeneous substrates

This paper presents a computational study of defect nucleation associated with the kinetics of the isotropic-to-nematic phase ordering transition over heterogeneous substrates, as it occurs in new liquid crystal biosensor devices, based on the Landau-de Gennes model for rod-like thermotropic nematic liquid crystals. Two regimes are identified due to interfacial tension inequalities: (i) nematic surface film nucleation and growth normal to the heterogeneous substrate, and (ii) nematic surface droplet nucleation and growth. The former, known as wetting regime, leads to interfacial defect shedding at the moving nematic-isotropic interface. The latter droplet regime, involves a moving contact line, and exhibits two texturing mechanisms that also lead to interfacial defect shedding: (a) small and large contact angles of drops spreading over a heterogeneous substrate, and (b) small drops with large curvature growing over homogeneous patches of the substrate. The numerical results are consistent with qualitative defect nucleation models based on the kinematics of the isotropic-nematic interface and the substrate-nematic-isotropic contact line. The results extend current understanding of phase ordering over heterogeneous substrates by elucidating generic defect nucleation processes at moving interfaces and moving contact lines.     

Interfacial nematodynamics of heterogeneous curved isotropic-nematic moving fronts

The early stages of liquid crystal phase ordering upon thermal quenches of isotropic phases into unstable and metastable temperature ranges is studied using two-dimensional (2D) computational solutions of the governing Landau-de Gennes (L-dG) equations for low molar mass nematic liquid crystals and analysis based on the corresponding interfacial nematodynamic model. The early phase ordering stage, for both unstable and metastable quenches of the isotropic phase, is shown to lead to highly textured nematic spherulites through a mechanism of interfacial defect nucleation. The underlying mechanisms of interface-driven texturing are elucidated using complementary 2D computational parametric studies of the bulk L-dG equation and analysis of the IN model. It is shown that for highly curved nanodomains and realistic elastic anisotropy, sharp interfacial transitions between uniaxial and biaxial states arise and are resolved by interfacial defect nucleation, which upon subsequent migration into the spherulite's interior leads to strong texturing. This paper shows that texture formation in the early stages of phase ordering is interface driven, and due to low interface tension, elastic anisotropy, and large curvature. Interfacial defect shedding in highly curved, low tension, anisotropic interfaces is a significant defect nucleation mechanism that needs to be taken into account when considering texturing processes.     

Preliminary communication Nematic-substrate repulsion in the nematic-isotropic phase coexistence region

Observations of two types of nematic droplet in the nematic-isotropic phase coexistence region are reported. One type contains topological defects and is free to move within a thin, homeotropically treated cell; the other is defect free and appears to be pinned at the substrates. The freely moving droplet represents an apparently new liquid crystal-substrate repulsion, which depends on the director alignments at the substrate and at the surface of the nematic droplet.     

Texture formation under phase ordering and phase separation in polymer-liquid crystal mixtures

Computational modeling of texture formation in coupled phase separation-phase ordering processes in polymer/liquid crystal mixtures is performed using a unified model based on the nematic tensor order parameter and gradient orientation elasticity. The computational methods are able to resolve defect nucleation, defect-defect interactions, and defect-particle interactions, as well as global and local morphological features in the concentration and order parameter spatiotemporal behavior. Biphasic structures corresponding to polymer dispersed liquid crystals (PDLCs), crystalline filled nematic (CFNs), and random filled nematics (RFNs) are captured and analyzed using liquid crystal defect physics and structure factors. Under spinodal decomposition due to concentration fluctuations, the PDLC structure emerges, and the nucleation and repulsive interaction of defects within nematic droplets leads to bipolar nematic droplets. Under spinodal decomposition due to ordering fluctuations, the CFNs structure emerges, and the stable polymer droplet crystal is pinned by a lattice of topological defects. For intermediate cases, where the mixture is unstable to both concentration and nematic order fluctuations, the RFN structure emerges, and polymer droplets and fibrils are pinned by a defect network, whose density increases with the curvature of the polymer-liquid crystal interface. The simulations provide an information of the role of topological defects on phase separation-phase ordering processes in polymer-liquid crystal mixtures.     

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.     

Theoretical simulation of thermally induced phase separation in a main‐chain liquid‐crystalline polymer solution

Phase diagrams of main‐chain liquid‐crystalline polymer (MCLCP) solutions have been calculated self‐consistently on the basis of a simple addition of the Flory–Huggins free energy for isotropic mixing, the Maier–Saupe free energy for nematic ordering, and the Flory free energy for chain rigidity of the MCLCP backbone. The calculated phase diagram is an upper critical solution type overlapping with the nematic–isotropic transition. The phase diagram consists of liquid–liquid, liquid–nematic, and pure nematic regions. Subsequently, the dynamics of thermally induced phase separation and morphology development have been investigated by the incorporation of the combined free energy density into the coupled time‐dependent Ginzburg–Landau (model C) equations, which involve conserved compositional and nonconserved orientational order parameters. The numerical calculations reveal a variety of the morphological patterns arising from the competition between liquid–liquid phase separation and nematic ordering of the liquid‐crystalline polymer. Of particular interest is the observation of an inflection in the growth dynamic curve, which may be attributed to the nematic ordering of the MCLCP component, which leads to the breakdown of the interconnected domains. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 41: 913–926, 2003     

Regimes of water droplet evaporation on copper substrates

Distilled water droplet evaporation has been studied on copper substrate surfaces with different degrees of roughness. Data on variations in the contact diameter have been employed to distinguish between the regimes of distilled water droplet spreading over the copper surfaces that proceed after the viscous regime. For each isolated regime, the duration has been determined as a fraction of the total evaporation time and the main physical processes have been described. Variations in contact angles have been analyzed as depending on copper surface temperature. It has been established that, as the substrate temperature is elevated, wetting becomes better, while the adhesion work remains almost unchanged, thereby indicating the absence of chemical and structural transformations at the liquid–substrate interface.     

Study of the surface and interfacial tensions in systems containing a low molar mass liquid crystal

The surface tension of a low molar mass liquid crystal (LMMLC), 4-cyano-4'-n-heptyloxybiphenyl (70CB), was measured as a function of temperature using the pendant drop method, forming drops of different volumes ranging from 5 to 11 mm³. Contact angles formed by drops of 70CB in the nematic and isotropic phases on plates of polystyrene (PS) and of a liquid crystal polymer (LCP), VECTRA A910, were also measured. Only large drops could be used for surface tension analysis. It was shown that in the nematic phase the surface tension of 70CB decreases with increasing temperature, and that in the isotropic phase the surface tension increases with increasing temperature. Using the values of contact angle and of surface tension of 7OCB it was possible to evaluate the interfacial energy between 7OCB and PS and between 7OCB and VECTRA. The interfacial energy between 7OCB and PS, and between 7OCB and VECTRA, decreased with increasing temperature for ranges of temperatures corresponding to both phases of 70CB.     

Study of the surface and interfacial tensions in systems containing a low molar mass liquid crystal

The surface tension of a low molar mass liquid crystal (LMMLC), 4-cyano-4'-n-heptyloxybiphenyl (70CB), was measured as a function of temperature using the pendant drop method, forming drops of different volumes ranging from 5 to 11 mm³. Contact angles formed by drops of 70CB in the nematic and isotropic phases on plates of polystyrene (PS) and of a liquid crystal polymer (LCP), VECTRA A910, were also measured. Only large drops could be used for surface tension analysis. It was shown that in the nematic phase the surface tension of 70CB decreases with increasing temperature, and that in the isotropic phase the surface tension increases with increasing temperature. Using the values of contact angle and of surface tension of 7OCB it was possible to evaluate the interfacial energy between 7OCB and PS and between 7OCB and VECTRA. The interfacial energy between 7OCB and PS, and between 7OCB and VECTRA, decreased with increasing temperature for ranges of temperatures corresponding to both phases of 70CB.     

Confinement of a liquid crystal to small droplets and its effect on nuclear magnetic relaxation

Abstract

The spatial dependence of the orientation of the molecular director and of the nematic order parameter is obtained by minimization of the Landau–de Gennes free energy of the nematic liquid crystal confined in a spherical droplet. Special attention is given to the vicinity of the nematic–isotropic transition. The influence of the resulting nematic structure, large liquid crystal–polymer interface and restricted molecular diffusion on the nuclear magnetic relaxation is analysed. The translationally-induced molecular reorientation and the liquid crystal–polymer cross relaxation are discussed in particular. The possibility of an indirect study of the molecular anchoring on the polymer surface is demonstrated.     

Nematic Homopolymers: From Segmented to Wormlike Chains

We present a density functional approach to orientational ordering in homopolymeric systems. The polymers are modeled as chains of identical rodlike segments connected via a simple generic bending potential. The segments are impenetrable to each other, and it is their mutual excluded volume that drives the transition from the orientationally disordered isotropic phase to the orientationally ordered nematic fluid. These excluded volume effects are accounted for within the so‐called Onsager approximation at the chain–chain level and in an independent pairwise overlap approximation at the segment–segment level. The Khokhlov and Semenov formalism for nematic wormlike polymers is shown to be an exact limiting case of our treatment. The ordering transition is studied analytically by using a linear stability analysis of the isotropic phase yielding the properties of the system at the isotropic‐nematic (I–N) bifurcation point. Using a numerical scheme, the equilibrium distribution functions in the nematic phase are calculated, and the location of the thermodynamic I–N transition is determined. For stiff bending potentials, chains with a relatively small number of segments are found to behave like wormlike chains, and we determine the regime of model parameters for which this identification holds.     

Heterogeneous nucleation of crystallization of high polymers from the melt. III. Nucleation kinetics and interfacial energies

In the melt crystallization of isotactic polypropylene, poly(ethylene oxide) and poly(butene-1) in contact with substrates, the existence of a fixed number of nucleating sites on the substrate surfaces has been established. When these sites become active successively (the transient in the number of nuclei is long) during crystallization, pseudohomogeneous nucleation on the substrate occurs. Nucleation rates for poly(butene-1) and poly(ethylene oxide) on substrates and in bulk have been measured. These data can be used for comparing the nucleating ability of substrates. Estimates of the variation of bulk nucleation rates from one volume element to another as well as for repeated crystallization within a given volume element have been included. Finally, the temperature coefficients of heterogeneous nucleation rates have been combined with the temperature coefficient of spherulitic growth rate of poly(butene-1), to yield values of the interfacial energy parameters appearing in the theory of heterogeneous nucleation. The quantitative characterization of the nucleating ability of substrates by this method is an improvement over the mere use of nucleation densities or nucleation rates.     

Topologically non-equivalent textures generated by the nematic electrohydrodynamics

Some classes of nematic liquid crystals can be driven through turbulent regimes when forced by an external electric field. In contrast to isotropic fluids, a turbulent nematic exhibits a transition to a stochastic regime that is characterised by a network of topological defects. We study the deformations arising after the electric field has been switched-off. In contrast to the turbulent regime, the relaxation of this topological-defect regime involves the annihilation of an interlacement of defect lines. We show that these defect lines separate regions of the nematic having topologically non-equivalent textures.     

Nematic–isotropic phase behaviors of polydisperse polymer/liquid‐crystal systems: Chain‐length‐dependent interaction parameter

Phase behaviors of polydisperse polystyrene (PS)/nematic liquid‐crystal systems [P‐ethoxy ‐ benzylidene ‐ p ‐ n‐butylaniline (EBBA)] are investigated with a thermo‐optical analysis technique. We also develop a thermodynamic framework to describe the phase behaviors of polydisperse PS/EBBA systems. The proposed model is based on a modified double‐lattice model to describe isotropic mixing and Maier–Saupe theory for anisotropic ordering. To correlate the polymer chain length and energy parameters in a nematic–isotropic biphasic region and to apply the primary interaction parameter in an isotropic–isotropic phase‐transition behaviors of polydisperse PS/EBBA systems. The proposed model shows remarkable agreement with experimental data for the model systems in comparison with an existing model. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 1031–1039, 2006     

Thermodynamics of heterogeneous crystal nucleation in contact and immersion modes

One of the most intriguing problems of heterogeneous crystal nucleation in droplets is its strong enhancement in the contact mode (when the foreign particle is presumably in some kind of contact with the droplet surface) compared to the immersion mode (particle immersed in the droplet). Heterogeneous centers can have different nucleation thresholds when they act in contact or immersion modes. The underlying physical reasons for this enhancement have remained largely unclear. In this paper we present a model for the thermodynamic enhancement of heterogeneous crystal nucleation in the contact mode compared to the immersion one. To determine if and how the surface of a liquid droplet can thermodynamically stimulate its heterogeneous crystallization, we examine crystal nucleation in the immersion and contact modes by deriving and comparing with each other the reversible works of formation of crystal nuclei in these cases. The line tension of a three-phase contact gives rise to additional terms in the formation free energy of a crystal cluster and affects its Wulff (equilibrium) shape. As an illustration, the proposed model is applied to the heterogeneous nucleation of hexagonal ice crystals on generic macroscopic foreign particles in water droplets at T = 253 K. Our results show that the droplet surface does thermodynamically favor the contact mode over the immersion one. Surprisingly, the numerical evaluations suggest that the line tension contribution (from the contact of three water phases (vapor-liquid-crystal)) to this enhancement may be of the same order of magnitude as or even larger than the surface tension contribution.     

Fluctuating Interfaces in Liquid Crystals

Summary: We review and compare recent work on the properties of fluctuating interfaces between isotropic and nematic liquid-crystalline phases. Molecular dynamics and Monte Carlo simulations have been carried out for systems of ellipsoids and hard rods with aspect ratio 15:1, and the fluctuation spectrum of interface positions (the capillary wave spectrum) has been analyzed. In addition, the capillary wave spectrum has been calculated analytically within the Landau-de Gennes theory. The theory predicts that the interfacial fluctuations can be described in terms of a wave vector dependent interfacial tension, which is anisotropic at small wavelengths (stiff director regime) and becomes isotropic at large wavelengths (flexible director regime). After determining the elastic constants in the nematic phase, theory and simulation can be compared quantitatively. We obtain good agreement for the stiff director regime. The crossover to the flexible director regime is expected at wavelengths of the order of several thousand particle diameters, which was not accessible to our simulations.     

Evolution equation for a disclination line located between the uniaxial and isotropic phases of a nematic liquid crystal

We derive a supplemental evolution equation for a disclination line located on an interface between the uniaxial and isotropic phases of a nematic liquid crystal. This equation provides an additional kinetic relation accounting for the motion of an interfacial disclination line. In our treatment of the problem we neglect fluid motion. Our approach is based on the notion of configurational forces. To illustrate the role of our additional evolution equation, we consider two simple examples. We also identify an expression for the configurational force exerted by the uniaxial phase at the defect located on the phase interface.     

Interfacial forces and Marangoni flow on a nematic drop retracting in an isotropic fluid

Nematic-isotropic interfaces exhibit novel dynamics due to anchoring of the liquid crystal molecules on the interface. The objective of this study is to demonstrate the consequences of such dynamics in the flow field created by an elongated nematic drop retracting in an isotropic matrix. This is accomplished by two-dimensional flow simulations using a diffuse-interface model. By exploring the coupling among bulk liquid crystal orientation, surface anchoring and the flow field, we show that the anchoring energy plays a fundamental role in the interfacial dynamics of nematic liquids. In particular, it gives rise to a dynamic interfacial tension that depends on the bulk orientation. Tangential gradient of the interfacial tension drives a Marangoni flow near the nematic-isotropic interface. Besides, the anchoring energy produces an additional normal force on the interface that, together with the interfacial tension, determines the movement of the interface. Consequently, a nematic drop with planar anchoring retracts more slowly than a Newtonian drop, while one with homeotropic anchoring retracts faster than a Newtonian drop. The numerical results are consistent with prior theories for interfacial rheology and experimental observations.