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.     

On the curvature dependence of the interfacial tension in a symmetric three-component interface

We consider a symmetric interface between two polymers A(N) and B(N) in a common monomeric solvent S using the mean-field Scheutjens-Fleer self-consistent field theory and focus on the curvature dependence of the interfacial tension. In multi-component systems there is not one unique scenario to curve such an interface. We elaborate on this by keeping either the chemical potential of the solvent or the bulk concentration of the solvent fixed, that is we focus on the semi-grand canonical ensemble case. Following Helfrich, we expand the surface tension as a Taylor series in the curvature parameters and find that there is a non-zero linear dependence of the interfacial tension on the mean curvature in both cases. This implies a finite Tolman length. In a thermodynamic analysis we prove that the non-zero Tolman length is related to the adsorption of solvent at the interface. Similar, but not the same, correlations between the solvent adsorption and the Tolman length are found in the two scenarios. This result indicates that one should be careful with symmetry arguments in a Helfrich analysis, in particular for systems that have a finite interfacial tension: one not only should consider the structural symmetry of the interface, but also consider the constraints that are enforced upon imposing the curvature. The volume fraction of solvent, the chain length N as well as the interaction parameter chi(AB) in the system can be used to take the system in the direction of the critical point. The usual critical behavior is found. Both the width of the interface and the Tolman length diverge, whereas the density difference between the two phases, adsorbed amount of solvent at the interface, interfacial tension, spontaneous curvature, mean bending modulus as well as the Gaussian bending modulus vanish upon approach of the critical point.     

Numerical analysis of the radial-axial structure transition with an applied field in a nematic droplet

In this paper the director configurations and the free energies of a nematic droplet with a surface normal anchoring condition are calculated numerically. For this surface anchoring, a transition occurs between the radial and axial structures with respect to an applied field. In the calculation of the director configurations, the position of a disclination has been fixed. Comparing the free energies for different disclinations, the stable position which gives the minimum free energy is found. In calculating the free energy of a droplet, it is assumed that the free energy density of the nematic phase does not exceed the isotropic free energy density, so that the large distortion in the vicinity of the disclination causes a nematic-isotropic transition and the free energy density of the disclination core becomes equal to the isotropic free energy density. The director configuration in a droplet is calculated as a function of an applied field for different isotropic free energy densities, elastic constant ratios and droplet shapes. The relation between the radial-axial structure transition and these factors are clarified.     

Twist viscosities and flow alignment of biaxial nematic liquid crystal phases of a soft ellipsoid-string fluid studied by molecular dynamics simulation

We have calculated the twist viscosity and the alignment angle between the director and the stream lines in shear flow of a liquid crystal model system, which forms biaxial nematic liquid crystals, as functions of the density, from the Green-Kubo relations by equilibrium molecular dynamics simulation and by a nonequilibrium molecular dynamics algorithm, where a torque conjugate to the director angular velocity is applied to rotate the director. The model system consists of a soft ellipsoid-string fluid where the ellipsoids interact according a repulsive version of the Gay-Berne potential. Four different length-to-width-to-breadth ratios have been studied. On compression, this system forms discotic or calamitic uniaxial nematic phases depending on the dimensions of the molecules, and on further compression a biaxial nematic phase is formed. In the uniaxial nematic phase there is one twist viscosity and one alignment angle. In the biaxial nematic phase there are three twist viscosities and three alignment angles corresponding to the rotation around the various directors and the different alignments of the directors relative to the stream lines, respectively. It is found that the smallest twist viscosity arises by rotation around the director formed by the long axes, the second smallest one arises by rotation around the director formed by the normals of the broadsides, and the largest one by rotation around the remaining director. The first twist viscosity is rather independent of the density whereas the last two ones increase strongly with density. One finds that there is one stable director alignment relative to the streamlines, namely where the director formed by the long axes is almost parallel to the stream lines and where the director formed by the normals of the broadsides is almost parallel to the shear plane. The relative magnitudes of the components of the twist viscosities span a fairly wide interval so this model should be useful for parameterisation experimental data.     

Nematic-isotropic interfaces under shear: a molecular-dynamics simulation

We present a large-scale molecular-dynamics study of nematic-paranematic interfaces under shear. We use a model of soft repulsive ellipsoidal particles with well-known equilibrium properties, and consider interfaces which are oriented normal to the direction of the shear gradient (common stress case). The director at the interface is oriented parallel to the interface (planar). A fixed average shear rate is imposed with moving periodic boundary conditions, and the heat is dissipated with a profile-unbiased thermostat. First, we study the properties of the interface at one particular shear rate in detail. The local interfacial profiles and the capillary wave fluctuations of the interfaces are calculated and compared with those of the corresponding equilibrium interface. Under shear, the interfacial width broadens and the capillary wave amplitudes at large wavelengths increase. The strain is distributed inhomogeneously in the system (shear banding), the local shear rate in the nematic region being distinctly higher than in the paranematic region. Surprisingly, we also observe (symmetry-breaking) flow in the vorticity direction, with opposite direction in the nematic and the paranematic state. Finally, we investigate the stability of the interface for other shear rates and construct a nonequilibrium phase diagram.     

Modelling of director equilibrium states in a nematic cell with relief surface

We study two-dimensional equilibrium configurations of nematic liquid crystal (NLC) director in a cell bounded by two parallel surfaces. One surface is planar and the other one is spatially modulated. The relief of the modulated surface is described by a smooth periodic sine-like function. The orientation of NLC director easy axis is assumed to be homeotropic at one cell surface and planar at the other one. The NLC director anchoring with cell surfaces is assumed to be strong. We consider the case where disclination lines occur in the bulk of NLC above the extrema of the modulated surface. These disclination lines run along the crests and troughs of the surface relief. If the orientation of director at both bounding surfaces is of the same type, then NLC director field is continuous. For both configurations mentioned above (with defects and without defects), we obtain analytical expressions for director distribution in the bulk of NLC in the approximation of planar director deformations. Equilibrium distances from disclination lines to the spatially modulated surface are calculated when the defects occur. The dependences of these equilibrium distances on the period and depth of surface relief and the cell thickness are investigated in detail.     

Magnetic field effects on tactoids of plate-like colloids

We investigate the effect of a magnetic field on the shape and director field of nematic droplets in dispersions of sterically stabilized and charge-stabilized colloidal gibbsite platelets with a negative diamagnetic anisotropy. Depending on the magnetic field strength and tactoid size, we observe with polarized light microscopy several interesting structures, with different shapes and director fields both with and without defects. In particular, our findings provide the first experimental evidence for the existence of the split-core defect structure predicted ten years ago by Mkaddem and Gartland [Phys. Rev. E 62, 6694 (2000)]. The split-core structure is a metastable director-field configuration that can be stabilized by a sufficiently strong externally applied magnetic field but only if the diamagnetic anisotropy of the particles is negative. To account for our observations, we present a calculation of the stability regions of different shapes and director-field structures as a function of tactoid size, anchoring conditions, surface tension, elastic constants, and magnetic field strength. By fitting the experimental data to the theoretically predicted structures, we are able to extract values for the splay elastic constant, interfacial tension, and anchoring strength. Remarkably, we find significant differences between the two systems studied: for sterically stabilized gibbsite in bromotoluene the anchoring strength is one order of magnitude larger than that of aqueous gibbsite, with the latter exhibiting weak and the former strong anchoring of the director field to the interface. The splay elastic constants that we obtain are in agreement with earlier experiments, simulations, and theory, while the interfacial tension and anchoring strength are considerably larger than what was found in earlier experiments.     

Streamline Image Analysis: a new tool for investigating defects in nematic liquid crystals

A bidimensional investigation of the director field near defects normal to the symmetry plane in a nematic liquid crystal is here presented, based upon a new 'streamline' image processing of the data obtained by polarized light microscopy. With a visualization of vector field streamlines related to the actual director configuration, the disclination topological strength is revealed. The streamline visualization is also suitable for detecting smooth local variations of the director of field as functions of both space and time. Moreover, the orientation of the director field projection in the cell plane is determined with high sensitivity ( 2%) and the possible presence of the out-of-plane escape of the director field is evaluated.     

Generalised technique for calculation of plane director profiles in bounded nematic liquid crystals

We propose an improved conformal mapping technique for analytical calculation of two-dimensional profiles of the nematic liquid crystal (NLC) director in a cell with a simply connected cross-sectional region. We consider the case of the strong anchoring and the piecewise constant director pretilt on the piecewise smooth curve which bounds the region. Obtained expressions for the director profile explicitly depend on the conformal mapping which maps the region onto the upper half plane of a complex plane. An advantage of our method in comparison with the standard conformal mapping technique is that it does not require the knowledge of the inverse mapping and the calculation of the integral in the Poisson formula. Proposed technique allows to take into account topological defects in the bulk of the NLC on the symmetry axis of the region. As an example of how the method can be used, we find an analytical expression for the director profile in a horizontal cylindrical groove partly filled with the NLC. We consider the case where a disclination line parallel to the axis of the groove occurs in the bulk of the NLC. The equilibrium position of the disclination line is found.     

Experimental measurement of the azimuthal anchoring energy function at a SiO-nematic interface

The functional form of the azimuthal anchoring energy, i.e. the anisotropic part of the interfacial free energy, at the interface between the nematic liquid crystal 4-n-pentyl-4'-cyanobiphenyl and an obliquely evaporated SiO substrate is measured for the first time by using a reflectometric method. The anchoring energy function is obtained by measuring the director rotation on the interface caused by an external magnetic field ranging from 0 to 2·3 T. The dependence of the anchoring energy on the director azimuthal angle is found to be well fitted by the function W_a(ϕ) = W_asin² ϕ in agreement with the predictions of the Berreman model for the anchoring at a grooved interface.     

On the compatibilization of polymer blends

By means of statistical thermodynamics we consider the effect of a diblock copolymer on the interface of a demixed homopolymer blend. In contrast to the usual case where the blocks are identical with the components of the homopolymers, we investigate the use of arbitrary blocks XY respective to the A-B blend. As examples, we calculate interfacial properties such as the interfacial tension and the width of the interface as well as the concentration profiles of the blocks. We expect a strong compatibilization effect if the blocks are sufficiently long and have a preferential repulsive interaction with one of the two homopolymers.     

Water molecule clusters measured at water/air interfaces using atomic force microscopy

During the tip approach to hydrophobic surfaces like the water/air interface, the measured interaction force reveals a strong attraction with a range of approximately 250 nm at some points along the interface. The range of this force is approximately 100 times larger than the measured for gold (approximately 3 nm) and 10 times larger than the one for hydrophobic silicon surfaces (approximately 25 nm). At other points the interface exerts a medium range repulsive force growing stepwise as the tip approaches the interface plane, consequently the hydrophobic force is a strong function of position. To explain these results we propose a model where the force on the tip is associated with the exchange of a small volume of the interface with a dielectric permittivity epsilon(int) by the tip with a dielectric permittivity epsilon(tip). By assuming a oscillatory spatial dependence for the dielectric permittivity it is possible to fit the measured force profiles. This dielectric spatial variation was associated with the orientation of the water molecules arrangement in the interfacial region. Small nanosized hydrogen-bond connected cages of water molecules present in bulk water at the interface are oriented by the interfacial electric field generated by the water molecules broken bonds, one broken hydrogen bond out of every four. This interfacial field orients small clusters formed by approximately 100 water molecules into larger clusters (approximately 100 nm). In the limit of small (less than 5 nm thick) water molecule cages we have modeled the static dielectric permittivity (epsilon) as the average response of those cages. In these regions the dielectric permittivity for water/air interfaces decreases monotonically from the bulk value epsilon approximately 80 to approximately 2 at the interface. For regions filled with medium size cages, the tip senses the structure of each cage and the static dielectric permittivity is matched to the geometrical features of these cages sized approximately 25 to 40 nm. Interfacial electric energy density values were calculated using the electric field intensity and the dielectric permittivity obtained by the fitting of the experimental points. The integration of the electric energy density along the interfacial region gives a value of 0.072 J m(-2) for interfacial energy density for points where the hydrophobic force has a range of approximately 250 nm. Regions formed by various clusters result in lower values of the interfacial energy density.     

Effect of quenched size polydispersity on the fluid-solid transition in charged colloidal suspensions

We study the effect of quenched size polydispersity on the phase behavior of charged colloidal suspensions using free-energy calculations in Monte Carlo simulations. The colloids are assumed to interact with a hard-core repulsive Yukawa (screened-Coulomb) interaction with constant surface potential, so that the particles are polydisperse both in size and charge. In addition, we take the size distribution to be fixed in both the fluid and crystal phase (no size fractionation is allowed). We study the fluid-solid transition for various screening lengths and surface potentials, finding that upon increasing the size polydispersity the freezing transition shifts toward higher packing fractions and the density discontinuity between the two coexisting phases diminishes. Our results provide support for a terminal polydispersity above which the freezing transition disappears.     

Effect of an electric field on defects in a nematic liquid crystal with variable surface anchoring

We report experimental studies on defects in a nematic liquid crystal with negative dielectric anisotropy mounted in a cell with perfluoropolymer-coated surfaces. The sample exhibits a discontinuous anchoring transition from planar to homeotropic on cooling at zero or a small electric field, and above a cross-over voltage a continuous ‘inverse Freedericksz transition’, at which the director starts tilting in opposite directions at the two surfaces. Defects of strength ±1/2 are either annihilated or expelled when the director tilts. On the other hand, disclination lines of ±1 which end in partial point defects (boojums) at the surfaces in the planar alignment regime acquire point defects of strength ±1 at the midplane of the cell when the director tilts. At a low enough temperature, the homeotropic anchoring becomes strong, and an electric field above the Freedericksz threshold generates the usual umbilic defects, which follow the dynamic scaling laws found in earlier studies.     

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.     

Direct decoration of disclinations by solidification-induced band texture for a nematic side chain liquid crystalline polymer

For a nematic polymethacrylate side chain liquid crystalline polymer, g 154 N 298 I (°C), the solidification-induced band texture has been observed aligned along the disclination under a polarizing optical microscope, when the specimen was quenched from 280°C to room temperature. The decoration technique of solidification-induced band texture, which is usually reported for main chain liquid crystalline polymers, was then introduced to reveal the director field pattern along a disclination for this side chain liquid crystalline polymer. It was found by infra-red dichroism measurements that the director orientation is parallel with the direction of the band. On this basis, disclinations with strength s=±1/2 and s=±1 were mapped according to the corresponding pattern of solidification-induced band texture. In addition, two types of inversion wall, loop-like and splay-type walls, were also found to be decorated by the solidification-induced band texture.     

Defect controlled dynamics of nematic liquids

The long time dynamical response of a nematic liquid exhibiting banded textures (inversion walls) during the twist Freedericksz transition is presented. A dynamical model of approach to equilibrium through defect interaction and the resulting dissolution of the banded textures is presented. A linear stability analysis shows that splay-bend inversion wall defects are unstable to two dimensional infinitesimal perturbations. A model of inversion wall segment collapse with production of a disclination line pair is given. The energy-momentum tensor gives the force of interaction between inversion walls and disclination lines. A perturbation analysis gives the evolution of the director field in closed form. Entropy production gives the velocity of each line. The growth law governing the wall dissolution is given.     

Direct decoration of disclinations by solidification-induced band texture and focal-conic texture for a low molar mass liquid crystal

Two decoration methods, solidification-induced band texture decoration and focal-conic texture decoration, were established to map the director field of disclinations in a low molar mass methacrylate liquid crystal. It was found that when the specimen film is quenched from the nematic melt to room temperature, solidification-induced band texture is observed arrayed along with the schlieren texture, and the orientation of the director field can be displayed. Moreover, when the specimen is cooled from the nematic melt to 63°C and annealed, the focal-conic texture of the smectic A phase is found to grow around the corresponding disclination core with good orientation to reveal the director field. By the two decoration techniques, the director fields of disclinations with strength s = 1/2 and s = ±1 were revealed. Two types of inversion wall, loop-like wall and splay-type wall, were found by both solidification-induced band texture decoration and focal-conic texture decoration.     

Influence of perpendicular external magnetic field on microstructures of monolayer composed of ferromagnetic particles: analysis by means of quasi-two-dimensional Monte Carlo simulation

We investigated the influences of the magnetic field strength and particle areal density on the microstructure of a quasi-two-dimensional monolayer composed of ferromagnetic particles by means of a Monte Carlo simulation. The magnetic field was applied along a direction perpendicular to the plane of the monolayer. Microstructures of the monolayer obtained in the simulations were analyzed in terms of radial distribution and orientational distribution functions. Formation of the microstructures is discussed from the perspective of particle-particle interaction energy and the perpendicular magnetic susceptibility of the monolayer was calculated from simulated magnetization curves. The obtained results are summarized as follows. For small areal density of particles, formation of chain-like structures is prevented by the repulsive magnetic interaction between particles due to orientations of the magnetic moments in the particles along the magnetic field direction. For intermediate areal density of particles, the chain-like structures remain even when a relatively strong magnetic field is applied, because contributions of the attractive magnetic interactions increase. For large areal density of particles, mixtures of chain-like and locally ordered structures appear due to the anisotropic attractive magnetic interactions in the absence of the magnetic field. However, when a sufficiently strong magnetic field is applied, the magnetic interactions between particles change to isotropic repulsive interactions, which results in the short-range repulsive steric interactions between particles becoming dominant with the appearance of hexagonal close packed structures.