Forces between elongated particles in a nematic colloid

Nontrivial topological defects of micro-rods immersed in nematics and their phototuning

Combinations of different geometries and surface anchoring conditions give rise to the diversity of topological structures in nematic colloid systems. Tuning these parameters in a single system offers possibilities for observing the evolution of the topological transformation and for manipulating colloids through topological forces. Here we investigate the nontrivial topological properties of micro-rods dispersed in nematic liquid crystals through experimental observation and computer simulation. The topological variation is driven by photodynamically changing the surface anchoring using azobenzene-based surface-commander molecules, the majority of which are localized on both the substrates and the surface of micro-rods. By comparing experimental and simulation results, we show previously unidentified topological properties of the two-body LC-rod-colloid system. Moreover, unlike the traditional photoresponsive liquid crystal systems, the localization of azobenzene molecules on the surfaces makes it possible to change only the direction of the surface orientation, not disordering of the bulk structures. The results assist in the development of photo-driven micro-robotics in fluids.

Field-triggered vertical positional transition of a microparticle suspended in a nematic liquid crystal cell

In this paper, based on the numerical calculation of total energy utilizing the Green's function method, we investigate how a field-triggered vertical positional transition of a microparticle suspended in a nematic liquid crystal cell is influenced by the direction of the applied field, surface anchoring feature, and nematic's dielectric properties. The new equilibrium position of the translational movement is decided via a competition between the buoyant force and the effective force built on the microparticle by the elastic energy gradient along the vertical direction. The threshold value of external field depends on thickness L and Frank elastic constant K and slightly on the microparticle size and density, in a Fréedericksz-like manner, but by a factor. For a nematic liquid crystal cell with planar surface alignment, a bistable equilibrium structure for the transition is found when the direction of the applied electric field is (a) perpendicular to the two plates of the cell with positive molecular dielectric anisotropy or (b) parallel to the two plates and the anchoring direction of the cell with negative molecular dielectric anisotropy. When the electric field applied is parallel to both plates and perpendicular to the anchoring direction, the microparticle suspended in the nematic liquid crystal tends to be trapped in the midplane, regardless of the sign of the molecular dielectric anisotropy. Such a phenomenon also occurs for negative molecular dielectric anisotropy if the external field is applied perpendicular to the two plates. Explicit formulas proposed for the critical electric field agree extremely well with the numerical calculation.

Defect-mediated colloidal interactions in a nematic-phase discotic solvent

Interactions between colloidal inclusions dispersed in a nematic discotic liquid-crystalline solvent were investigated for different solute-solvent coupling conditions. The solvent was treated at the level of Gay-Berne discogens. Colloidal inclusions were coupled to the solvent with a generalized sphere-ellipsoid interaction potential. Energy strengths were varied to promote either homeotropic or planar mesogenic anchoring. Colloid-colloid interactions were modeled using a soft, excluded-volume contribution. Single-colloid and colloid-pair samples were evolved with Molecular Dynamics simulations. Equilibrium trajectories were used to characterize structural and dynamical properties of topological defects arising in the mesomorphic phase due to colloidal inclusions. Boojums were observed with planar anchoring, whereas Saturn rings were obtained with homeotropic anchoring. The manner in which these topological defects drive colloidal interactions was assessed through a free energy analysis, taking into account the relative orientation between a colloidal dyad and the nematic-field director. The dynamical behavior of defects was qualitatively surveyed from equilibrium trajectories borne from computer simulations.

Segregation of liquid crystal mixtures in topological defects

The structure and physical properties of liquid crystal (LC) mixtures are a function of composition, and small changes can have pronounced effects on observables, such as phase-transition temperatures. Traditionally, LC mixtures have been assumed to be compositionally homogenous. The results of chemically detailed simulations presented here show that this is not the case; pronounced deviations of the local order from that observed in the bulk at defects and interfaces lead to significant compositional segregation effects. More specifically, two disclination lines are stabilized in this work by introducing into a nematic liquid crystal mixture a cylindrical body that exhibits perpendicular anchoring. It is found that the local composition deviates considerably from that of the bulk at the interface with the cylinder and in the defects, thereby suggesting new assembly and synthetic strategies that may capitalize on the unusual molecular environment provided by liquid crystal mixtures.

Liquid crystal mixtures are used in commercial applications and their composition affects their properties. Here Rahimi et al. use atomistic simulations to show that defects influence the molecular arrangement of the mixture components leading to a deviation of the local order from that of the bulk.

Theory of elastic interaction between arbitrary colloidal particles in confined nematic liquid crystals

We develop the method proposed by Chernyshuk and Lev [Phys. Rev. E 81, 041701 (2010)] for theoretical investigation of elastic interactions between colloidal particles of arbitrary shape and chirality (polar as well as azimuthal anchoring) in the confined nematic liquid crystal (NLC). General expressions for six different types of multipole elastic interactions are obtained in the confined NLC: monopole-monopole (Coulomb type), monopole-dipole, monopole-quadrupole, dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole interactions. The obtained formulas remain valid in the presence of the external electric or magnetic fields. The exact equations are found for all multipole coefficients for the weak anchoring case. For the strong anchoring coupling, the connection between the symmetry of the shape or director and multipole coefficients is obtained, which enables us to predict which multipole coefficients vanish and which remain nonzero. The particles with azimuthal helicoid anchoring are considered as an example. Dipole-dipole interactions between helicoid cylinders and cones are found in the confined NLC. In addition, the banana-shaped particles in homeotropic and planar nematic cells are considered. It is found that the dipole-dipole interaction between banana-shaped particles differs greatly from the dipole-dipole interaction between the axially symmetrical particles in the nematic cell. There is a crossover from attraction to repulsion between banana particles along some directions in nematic cells. It is shown that monopoles do not "feel" the type of nematic cell: monopole-monopole interaction turns out to be the same in homeotropic and planar nematic cells and converges to the Coulomb law as thickness increases, L→∞.

Elastic and hydrodynamic torques on a colloidal disk within a nematic liquid crystal

The orientationally dependent elastic energy and hydrodynamic behavior of colloidal disks with homeotropic surface anchoring suspended in the nematic liquid crystal 4-cyano-4'-pentylbiphenyl (5CB) have been investigated. In the absence of external torques, the disks align with the normal of the disk face â parallel to the nematic director n[over ^]. When a magnetic field is applied, the disks rotate â by an angle θ so that the magnetic torque and the elastic torque caused by distortion of the nematic director field are balanced. Over a broad range of angles, the elastic torque increases linearly with θ in quantitative agreement with a theoretical prediction based on an electrostatic analogy. When the disks are rotated to angles θ>π/2, the resulting large elastic distortion makes the disk orientation unstable, and the director undergoes a topological transition in which θ→π-θ. In the transition, a defect loop is shed from the disk surface, and the disks spin so that â sweeps through π radians as the loop collapses back onto the disk. Additional measurements of the angular relaxation of disks to θ=0 following removal of the external torque show a quasi-exponential time dependence from which an effective drag viscosity for the nematic can be extracted. The scaling of the angular time dependence with disk radius and observations of disks rotating about â indicate that the disk motion affects the director field at surprisingly modest Ericksen numbers.

Theory of elastic interaction between colloidal particles in a nematic cell in the presence of an external electric or magnetic field

The Green's function method developed previously [S. B. Chernyshuk and B. I. Lev, Phys. Rev. E 81, 041701 (2010)] is used to describe elastic interactions between axially symmetric colloidal particles in a nematic cell in the presence of an external electric or magnetic field. Formulas for dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole interactions in the homeotropic and planar nematic cells with parallel and perpendicular field orientations are obtained. A set of predictions has been made: (1) The deconfinement effect for dipole particles in the homeotropic nematic cell when an electric field is approaching its Freedericksz threshold value E⇒E(t). This means cancellation of the confinement effect found in [M. Vilfan et al., Phys. Rev. Lett. 101, 237801 (2008)] near the Freedericksz transition. In the planar nematic cell this deconfinement effect exists for both dipole and quadrupole particles and depends on the field orientation as well as on the sign of dielectric anisotropy Δε. (2) The effect of tunable stabilization of the particles is predicted. The equilibrium distance between two particles, which are attracted along the electric field parallel to the planes of a homeotropic nematic cell with Δε<0, depends on the strength of the field. (3) Attraction and repulsion zones for all elastic interactions are changed dramatically under the action of the external field.

Theory of elastic interaction of colloidal particles in nematic liquid crystals near one wall and in the nematic cell

We apply the method developed [Chernyshuk and Lev, Phys. Rev. E 81, 041701 (2010)] for theoretical investigation of colloidal elastic interactions between axially symmetric particles in the confined nematic liquid crystal near one wall and in the nematic cell with thickness L. Both cases of homeotropic and planar director orientations are considered. Particularly, dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole interactions of the one particle with the wall and within the nematic cell are found as well as corresponding two particle elastic interactions. A set of results has been predicted: The effective power of repulsion between two dipole particles at height h near the homeotropic wall is reduced gradually from inverse 3 to 5 with an increase of dimensionless distance r / h; near the planar wall, the effect of dipole-dipole isotropic attraction is predicted for large distances r > r(dd) = 4.76 h; maps of attraction and repulsion zones are crucially changed for all interactions near the planar wall and in the planar cell; and one dipole particle in the homeotropic nematic cell was found to be shifted by the distance δ(eq) from the center of the cell. The proposed theory fits very well with experimental data for the confinement effect of elastic interaction between spheres in the homeotropic cell [Vilfan et al., Phys. Rev. Lett. 101, 237801 (2008)] in the range 1-1000 kT. The influence of the K(24) and K(13) terms as well as connection with other theoretical approaches are discussed.

Long-range elastic-mediated interaction between nanoparticles adsorbed on free-standing smectic films

We determine the elastic-mediated interaction between colloidal nanoparticles adsorbed on the surface of free-standing smectic films. In contrast with the short-range character of the elastic-mediated force between particles adsorbed on smectic films supported by a solid substrate, the effective force acquires a long-range character in free-standing films, decaying with the particles distance R as slow as 1/R . We also discuss the dependence of the effective interaction potential on the surface tension gamma and film thickness. We show that it decays as 1/gamma in the regime of large surface tensions and becomes independent of the film thickness at a characteristic surface tension.

Forces between cylindrical nanoparticles in a liquid crystal

Using classical density functional theory, the forces between two cylindrical nanoparticles in a liquid crystal solvent are calculated. Both the nematic and isotropic phases of the solvent are considered. In the nematic phase, the interaction is highly anisotropic. At short range, changes in the defect structure around the cylinders leads to a complex interaction between them. In the isotropic phase, an attractive interaction arises due to overlap between halos of ordered fluid adsorbed on the surfaces of the cylinders.

Effective triplet interactions in nematic colloids

Three-body effective interactions emerging between parallel cylindrical rods immersed in a nematic liquid crystals are calculated within the Landau-de Gennes free-energy description. Collinear, equilateral and midplane configurations of the three colloidal particles are considered. In the last two cases the effective triplet interaction is of the same magnitude and range as the pair one.

Depletion effects in smectic phases of hard-rod-hard-sphere mixtures

It is known that when hard spheres are added to a pure system of hard rods the stability of the smectic phase may be greatly enhanced, and that this effect can be rationalised in terms of depletion forces. In the present paper we first study the effect of orientational order on depletion forces in this particular binary system, comparing our results with those obtained adopting the usual approximation of considering the rods parallel and their orientations frozen. We consider mixtures with rods of different aspect ratios and spheres of different diameters, and we treat them within Onsager theory. Our results indicate that depletion effects, and consequently smectic stability, decrease significantly as a result of orientational disorder in the smectic phase when compared with corresponding data based on the frozen-orientation approximation. These results are discussed in terms of the tau parameter, which has been proposed as a convenient measure of depletion strength. We present closed expressions for tau, and show that it is intimately connected with the depletion potential. We then analyse the effect of particle geometry by comparing results pertaining to systems of parallel rods of different shapes (spherocylinders, cylinders and parallelepipeds). We finally provide results based on the Zwanzig approximation of a fundamental-measure density-functional theory applied to mixtures of parallelepipeds and cubes of different sizes. In this case, we show that the tau parameter exhibits a linear asymptotic behaviour in the limit of large values of the hard-rod aspect ratio, in conformity with Onsager theory, as well as in the limit of large values of the ratio of rod breadth to cube side length, d, in contrast to Onsager approximation, which predicts tau approximately d (3). Based on both this result and the Percus-Yevick approximation for the direct correlation function for a hard-sphere binary mixture in the same limit of infinite asymmetry, we speculate that, for spherocylinders and spheres, the tau parameter should be of order unity as d tends to infinity.

Anisotropic nanoparticles immersed in a nematic liquid crystal: defect structures and potentials of mean force

We report results for the potential of mean force (PMF) and the defect structures that arise when spherocylindrical nanoparticles are immersed in a nematic liquid crystal. Using a dynamic field theory for the tensor order parameter Q of the liquid crystal, we analyzed configurations, including one, two, and three elongated particles, with strong homeotropic anchoring at their surfaces. For systems with one nanoparticle, the most stable configuration is achieved when the spherocylinder is placed with its long axis perpendicular to the far-field director, for which the defect structure consists of an elongated Saturn ring. For systems with two or three nanoparticles with their long axes placed perpendicular to the far-field director, at small separations the defect structures consist of incomplete Saturn rings fused with new disclination rings orthogonal to the original ones, in analogy to results previously observed for spherical nanoparticles. The shape of these orthogonal rings depends on the nanoparticles' configuration, i.e., triangular, linear, or parallel with respect to their long axis. A comparison of the PMFs indicates that the latter configuration is the most stable. The stability of the different arrays depends on whether orthogonal disclination rings form or not, their size, and the curvature effects in the interparticle regions. Our results suggest that the one-elastic-constant approximation is valid for the considered systems; similar results were obtained when a three-constant expression is used to represent the elastic free energy. The attractive interactions between the elongated particles were compared to those observed for spheres of similar diameters. Similar interparticle energies were observed for linear arrays; in contrast, parallel and triangular arrays of spherocylinders yielded interactions that were up to 3.4 times stronger than those observed for spherical particles.

Quenched disorder in a liquid-crystal biosensor: Adsorbed nanoparticles at confining walls

We analyze the response of a nematic liquid-crystal film, confined between parallel walls, to the presence of nanoscopic particles adsorbed at the walls. This is done for a variety of patterns of adsorption (random and periodic) and operational conditions of the system that can be controlled in experimental liquid-crystal-based devices. We compute simulated optical textures and the total optical output of the sensor between crossed polars, as well as the correlation function for the liquid-crystal tensor order parameter; we use these observables to discuss the gradual destruction of the original uniform orientation. For large concentrations of particles adsorbed in random patterns, the liquid crystal at the center of the sensor adopts a multidomain state, characterized by a small correlation length of the tensor order parameter, and also by a loss of optical anisotropy under observation through crossed polars. In contrast, for particles adsorbed in periodic patterns, the nematic at the center of the cell can remain in a monodomain orientation state, provided the patterns in opposite walls are synchronized.

Key-lock mechanism in nematic colloidal dispersions

We consider the interaction between two-dimensional nematic colloids and planar or sculpted walls. The elastic interaction between colloidal disks and flat walls, with homeotropic boundary conditions, is always repulsive. These repulsions may be turned into strong attractions at structured or sculpted walls, with cavities that match closely the shape and size of the colloids. This key-lock mechanism is analyzed in detail for colloidal disks and spherocylindrical cavities of various length to depth ratios, by minimizing the Landau-de Gennes free energy functional of the nematic orientational order parameter. We find that the attractions occur only for walls with cavities within a small range of the colloidal size and a narrow range of orientations with respect to the cavity's symmetry axis.

Potential of mean force between a spherical particle suspended in a nematic liquid crystal and a substrate: sphere size effects

The expanded ensemble density of states method (ExEDOS) is used to investigate the effective interaction of a spherical colloidal particle suspended in a confined liquid crystal (LC) with a substrate. The potential of mean force (PMF) is determined as a function of the normal distance between the particle and the substrate's surface. The presence of the substrate induces a layered structure of the LC, which in turn greatly influences the PMF. We analyze the structure of the Saturn ring defect that accompanies the colloidal sphere, and find that the ring is displaced slightly towards the surface when the sphere is within the first LC surface layer. A transition occurs from an overall attraction of the colloid to the substrate to a global repulsion when the sphere's radius is roughly twice the length of the LC molecules.

Interaction between two spherical particles in a nematic liquid crystal

We numerically investigate the interaction between two spherical particles in a nematic liquid crystal mediated by elastic distortions in the orientational order. We pay attention to the cases where two particles with equal radii R0 impose rigid normal anchoring on their surfaces and carry a pointlike topological defect referred to as a hyperbolic hedgehog. To describe the geometry of our system, we use bispherical coordinates, which prove useful in the implementation of boundary conditions at the particle surfaces and at infinity. We adopt the Landau-de Gennes continuum theory in terms of a second-rank tensor order parameter Q(ij) for the description of the orientational order of a nematic liquid crystal. We also utilize an adaptive mesh refinement scheme that has proven to be an efficient way of dealing with topological defects whose core size is much smaller than the particle size. When the two "dipoles," composed of a particle and a hyperbolic hedgehog, are in parallel directions, the two-particle interaction potential is attractive for large interparticle distances D and proportional to D-3 as expected from the form of the dipole-dipole interaction, until the well-defined potential minimum at D approximately 2.46 R0 is reached. For the antiparallel configuration with no hedgehogs between the two particles, the interaction potential is repulsive and behaves as D-2 for D less than or approximately equal 10R0, which is stronger than the dipole-dipole repulsion (approximately D-3 ) expected theoretically as an asymptotic behavior for large D.

Interaction of colloids with a nematic-isotropic interface

The Landau-de Gennes free energy is used to calculate the interaction between long cylindrical colloids and the nematic-isotropic (NI) interface. This interaction has two contributions: one is specific of liquid crystals and results from the deformation of the director field close to the particles or to the interface, while the other is generic and results from wetting and surface tension effects. Deep in the nematic phase the director field of long cylindrical colloids, with strong homeotropic anchoring, exhibits two half-integer defect lines. As the colloid moves towards the interface, the director configuration changes through a series of discontinuous transitions, where one or two of the defects are annihilated. In addition, the NI interface bends towards the colloid in order to minimize the elastic free energy in the nematic. In the isotropic phase, the colloid is surrounded by a thin nematic layer that reduces the surface free energy under favorable wetting conditions. The interaction has a well-defined minimum near the interface. In this region the director and interfacial structures are complex and cannot be described analytically. Using the numerical results for the Landau-de Gennes free energy in the harmonic region, we obtained simple scaling laws for the (linear) force on the colloid.