Computer simulation of topological defects around a colloidal particle or droplet dispersed in a nematic host

Lock-key microfluidics: simulating nematic colloid advection along wavy-walled channels

Liquid crystalline media mediate interactions between suspended particles and confining geometries, which not only has potential to guide patterning and bottom-up colloidal assembly, but can also control colloidal migration in microfluidic devices. However, simulating such dynamics is challenging because nemato-elasticity, diffusivity and hydrodynamic interactions must all be accounted for within complex boundaries. We model the advection of colloids dispersed in flowing and fluctuating nematic fluids confined within 2D wavy channels. A lock-key mechanism between homeotropic colloids and troughs is found to be stronger for planar anchoring on the wavy walls compared to homeotropic anchoring on the wavy walls due to the relative location of the colloid-associated defects. Sufficiently large amplitudes result in stick-slip trajectories and even permanent locking of colloids in place. These results demonstrate that wavy walls not only have potential to direct colloids to specific docking sites but also to control site-specific resting duration and intermittent elution.

Reduced-variance orientational distribution functions from torque sampling

We introduce a method to sample the orientational distribution function in computer simulations. The method is based on the exact torque balance equation for classical many-body systems of interacting anisotropic particles in equilibrium. Instead of the traditional counting of events, we reconstruct the orientational distribution function via an orientational integral of the torque acting on the particles. We test the torque sampling method in two- and three-dimensions, using both Langevin dynamics and overdamped Brownian dynamics, and with two interparticle interaction potentials. In all cases the torque sampling method produces profiles of the orientational distribution function with better accuracy than those obtained with the traditional counting method. The accuracy of the torque sampling method is independent of the bin size, and hence it is possible to resolve the orientational distribution function with arbitrarily small angular resolutions.

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.

Topology of Orientational Defects in Confined Smectic Liquid Crystals

We propose a general formalism to characterize orientational frustration of smectic liquid crystals in confinement by interpreting the emerging networks of grain boundaries as objects with a topological charge. In a formal idealization, this charge is distributed in pointlike units of quarter-integer magnitude, which we identify with tetratic disclinations located at the end points and nodes. This coexisting nematic and tetratic order is analyzed with the help of extensive Monte Carlo simulations for a broad range of two-dimensional confining geometries as well as colloidal experiments, showing how the observed defect networks can be universally reconstructed from simple building blocks. We further find that the curvature of the confining wall determines the anchoring behavior of grain boundaries, such that the number of nodes in the emerging networks and the location of their end points can be tuned by changing the number and smoothness of corners, respectively.

Phase-field model for elastocapillary flows of liquid crystals

We propose a phase-field model to study interfacial flows of nematic liquid crystals that couple the capillary forces on the interface with the elastic stresses in the nematic phase. The theoretical model has two key ingredients: A tensor order parameter that provides a consistent description of the molecular and distortional elasticity, and a phase-field formalism that accurately represents the interfacial tension and the nematic anchoring stress by approximating a sharp-interface limit. Using this model, we carry out finite-element simulations of drop retraction in a surrounding fluid, with either component being nematic. The results are summarized by eight representative steady-state solutions in planar and axisymmetric geometries, each featuring a distinct configuration for the drop and the defects. The dynamics is dominated by the competition between the interfacial tension and the distortional elasticity in the nematic phase, mediated by the anchoring condition on the drop surface. As consequences of this competition, the steady-state drop deformation and the clearance between the defects and the drop surface both depend linearly on the elastocapillary number.

Topological defects in an unconfined nematic fluid induced by single and double spherical colloidal particles

We present numerical solutions to the Landau-de Gennes free-energy model under the one-constant approximation for systems of single and double spherical colloidal particles immersed in an otherwise uniformly aligned nematic liquid crystal. A perfect homeotropic surface anchoring of liquid-crystal molecules on the spherical surface is considered. A large parameter space is carefully examined, including those in the free-energy model and those describing the dimer configurations and the background liquid-crystal orientation. The stability of the resulting liquid-crystal defects appearing in the neighborhood of the colloidal dimer pair is analyzed in light of the numerical results for their free energies. A number of scenarios are considered: a free dimer pair in a nematic fluid where the free-energy ground states are described in terms of a phase diagram, and a constrained dimer pair where the interparticle distance and the relative orientation of the distance vector to the nematic director can be manipulated. We pay particular attention to the nonsymmetric solutions, which yield several metastable defect states that can be observed in real systems. The high-precision numerical calculations are based on a spectral method, which is an enabling factor that allows us to compare the subtle difference in the free energies of different defect structures.

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.

The Impact of Colloidal Surface-Anchoring on the Smectic A Phase

Liquid-crystalline phases are known for their unique properties, i.e., the combination of fluidity and long-range orientational and/or positional order. The presence of a colloidal particle gives rise to perturbations of this order locally. These perturbations are the origin of intercolloidal forces driving the colloidal self-assembly in a directed manner. Hence, the understanding of these perturbations is the first step in understanding and controlling the self-assembly process. Here, we perform Monte Carlo simulations to investigate the perturbations of orientational and positional order in a smectic A phase caused by a spherical colloid. We model the host phase via an interaction potential that reproduces characteristic features of phase behavior, structure, dynamics, and elasticity [S. Püschel-Schlotthauer et al. J. Chem. Phys. 2016, 145, 164903]. For strong homeotropic anchoring conditions, we find a Saturn ring defect and an onion structure in the smectic A phase; the latter has never been reported for colloids so far. For strong planar anchoring conditions, we find Boojum defects that become elongated at low temperature, similar to what is observed in experiments. However, for weak planar anchoring conditions, a double surface ring defect is exhibited in the smectic A phase.

Shear dynamics of an inverted nematic emulsion

Here we study theoretically the dynamics of a 2D and a 3D isotropic droplet in a nematic liquid crystal under a shear flow. We find a large repertoire of possible nonequilibrium steady states as a function of the shear rate and of the anchoring of the nematic director field at the droplet surface. We first discuss homeotropic anchoring. For weak anchoring, we recover the typical behaviour of a sheared isotropic droplet in a binary fluid, which rotates, stretches and can be broken by the applied flow. For intermediate anchoring, new possibilities arise due to elastic effects in the nematic fluid. We find that in this regime the 2D droplet can tilt and move in the flow, or tumble incessantly at the centre of the channel. For sufficiently strong anchoring, finally, one or both of the topological defects which form close to the surface of the isotropic droplet in equilibrium detach from it and get dragged deep into the nematic state by the flow. In 3D, instead, the Saturn ring associated with the normal anchoring disclination line can be deformed and shifted downstream by the flow, but remains always localized in the proximity of the droplet, at least for the parameter range we explored. Tangential anchoring in 2D leads to a different dynamic response, as the boojum defects characteristic of this situation can unbind from the droplet under a weaker shear with respect to the normal anchoring case. Our results should stimulate further experiments with inverted liquid crystal emulsions under shear, as most of the predictions can be testable in principle by monitoring the evolution of liquid crystalline orientation patterns or by tracking the position and shape of the droplet over time.

Analytical description of the Saturn-ring defect in nematic colloids

We derive an analytical formula for the Saturn-ring configuration around a small colloidal particle suspended in nematic liquid crystal. In particular we obtain an explicit expression for the ring radius and its dependence on the anchoring energy. We work within Landau-de Gennes theory: Nematic alignment is described by a tensorial order parameter. For nematic colloids this model had previously been used exclusively to perform numerical computations. Our method demonstrates that the tensorial theory can also be used to obtain analytical results, suggesting a different approach to the understanding of nematic colloidal interactions.

Orientational ordering of Janus colloids in cholesteric liquid crystals

In this paper we show that Janus colloids, which are spherical particles with hybrid anchoring conditions, have preferable orientations in cholesteric media depending on the cholesteric wave vector. Simulations reveal that the tilt angle of a particle varies greatly with variation of the particle diameter to the cholesteric pitch ratio, which makes it possible to stabilize the appropriate particle orientation and to control it by variation of the cholesteric pitch.

Disclination lines at homogeneous and heterogeneous colloids immersed in a chiral liquid crystal

In the present work we perform Monte Carlo simulations in the isothermal-isobaric ensemble to study defect topologies formed in a cholesteric liquid crystal due to the presence of a spherical colloidal particle. Topological defects arise because of the competition between anchoring at the colloidal surface and the local director. We consider homogeneous colloids with either local homeotropic or planar anchoring to validate our model by comparison with earlier lattice Boltzmann studies. Furthermore, we perform simulations of a colloid in a twisted nematic cell and discuss the difference between induced and intrinsic chirality on the formation of topological defects. We present a simple geometrical argument capable of describing the complex three-dimensional topology of disclination lines evolving near the surface of the colloid. The presence of a Janus colloid in a cholesteric host fluid reveals a rich variety of defect structures. Using the Frank free energy we analyze these defects quantitatively indicating a preferred orientation of the Janus colloid relative to the cholesteric helix.

Defect topologies in a nematic liquid crystal near a patchy colloid

Using isothermal-isobaric Monte Carlo simulations we investigate defect topologies due to a spherical colloidal particle immersed in a nematic liquid crystal. Defects arise because of the competition between the preferential orientation at the colloid's surface and the far-field director ̂n(0). Considering a chemically homogeneous colloid as a special case we observe the well-known surface and saturn ring defect topologies for weak and strong perpendicular anchoring, respectively; for homogeneous, strong parallel anchoring we find a boojum defect topology that has been seen experimentally [see P. Poulin and D. A. Weitz, Phys. Rev. E 57, 626 (1998)] but not in computer simulations. We also consider a heterogeneous, patchy colloid where the liquid-crystal molecules anchor either preferentially planar or perpendicular at the surface of the colloid. For a patchy colloid we observe a boojum ring defect topology in agreement with recent experimental studies [see M. Conradi, M. Ravnik, M. Bele, M. Zorko, S. Žumer, and I. Muševič, Soft Matter 5, 3905 (2009)]. We also observe two other novel defect topologies that have not been reported thus far neither experimentally nor theoretically.

Defect structures in nematic liquid crystals around charged particles

We numerically study the orientation deformations in nematic liquid crystals around charged particles. We set up a Ginzburg-Landau theory with inhomogeneous electric field. If the dielectric anisotropy epsilon 1 is positive, Saturn-ring defects are formed around the particles. For epsilon 1< 0 , novel "ansa" defects appear, which are disclination lines with their ends on the particle surface. We find unique defect structures around two charged particles. To lower the free energy, oppositely charged particle pairs tend to be aligned in the parallel direction for epsilon 1> 0 and in the perpendicular plane for epsilon 1< 0 with respect to the background director. For identically charged pairs the preferred directions for epsilon 1> 0 and epsilon 1< 0 are exchanged. We also examine competition between the charge-induced anchoring and the short-range anchoring. If the short-range anchoring is sufficiently strong, it can be effective in the vicinity of the surface, while the director orientation is governed by the long-range electrostatic interaction far from the surface.

Different mechanisms of nucleation and self-organization of droplets in ferroelectric smectic membranes

New mechanisms of droplet nucleation and self-organization in ferroelectric membranes are described. The droplets may be accompanied by different number of topological defects (zero, one, two) whose location may be on the droplet boundary or in the membrane. Nucleation and self-organization of droplets with total topological charge S = 0 , S = + 1 and S = - 1 were investigated. We found that an S = - 1 topological defect may be the center of both droplet nucleation and chain formation. This mechanism of chaining drastically differs from the droplet self-organization described earlier which is realized by attraction of droplet-defect pairs. Our observations demonstrate new possibilities for manipulating the inclusions and their self-organization in smectic membranes.

Liquid-crystal-mediated force between a cylindrical nanoparticle and substrate

Using classical density functional theory, the structure of a molecular fluid around a cylindrical nanoparticle near a solid substrate is studied. The solvent-mediated force between the nanoparticle and the substrate is calculated in both the nematic and isotropic phases of the solvent. In the nematic phase, the force is short ranged and arises due to interaction between high-density regions near the substrate and nanoparticle. In the isotropic phase, the formation of a nematic bridge between the substrate and nanoparticle gives rise to an attractive force between them. The potential between the nanoparticle and substrate as a function of separation calculated numerically is compared to that calculated from the Derjaguin approximation. In the isotropic phase these are found to be in reasonable agreement at low separations, while the agreement is poorer in the nematic phase.

Director-configurational transitions around microbubbles of hydrostatically regulated size in liquid crystals

A high-pressure technique is introduced which allows a continuous variation of the inclusion size in liquid crystal colloids. We use a nematic liquid crystal host into which micrometer-sized gas bubbles are injected. By applying hydrostatic pressures, the diameter of these gas bubbles can be continuously decreased via compression and absorption of gas into the host liquid crystal, so that the director configurations around a single bubble can be investigated as a function of the bubble size. The theoretically predicted transition from a hyperbolic hedgehog to a Saturn-ring configuration, on reduction of the particle size below a certain threshold, is confirmed to occur at the radius of a few micrometers.

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.

Lattice Boltzmann algorithm to simulate isotropic-nematic emulsions

We present lattice Boltzmann simulations of the dynamical equations of motion of a drop of isotropic fluid in a nematic liquid crystal solvent, both in the absence and in the presence of an electric field. The coupled equations we solve are the Beris-Edward equations for the dynamics of the tensor order parameter describing the nematic solvent, the Cahn-Hilliard equation for the concentration evolution, and the Navier-Stokes equations for the determination of the instantaneous velocity field. We implement the lattice Boltzmann algorithm to ensure that spurious velocities are close to zero in equilibrium. We first study the effects of the liquid crystal elastic constant, K, anchoring strength, W, and surface tension, sigma, on the shape of the droplet and on the director field texture in equilibrium. We then consider how the drop behaves as the director field is switched by an applied electric field. We also show that the algorithm allows us to follow the motion of a drop of isotropic fluid placed in a liquid crystal cell with a tilted director field at the boundaries.

Molecular dynamics simulations studies of nanoparticles in an isotropic liquid crystal matrix: Single particle behavior and pairwise interactions

We report results of molecular dynamics simulation studies of the behavior of spherical nanoparticles (NPs) in a dense isotropic nematogen matrix comprised of soft spherocylinders (SSCs). The SSCs exhibit a tendency for frustrated planar anchoring at the NP surface that results in a long-range (compared to the size of the NPs and SSCs) reduction in local orientational ordering and increased fluctuations in local orientational ordering compared to the pure isotropic phase of the SSCs. The potential of mean force between two nanoparticles exhibits a novel long-range repulsive tail separated from short-range molecular packing peaks by a shallow local minimum in free energy. The long-range repulsion is caused by NP-induced ordering fluctuations while the shallow minimum results from increased local ordering within the confinement region in between two NPs. The influence of the NPs on local orientational order in the nematogen matrix and the nematogen-induced interaction between NPs are found to depend strongly on the size of the NPs.

Charged inclusion in nematic liquid crystals

We present a general theory of liquid crystals under inhomogeneous electric field in a Ginzburg-Landau scheme. The molecular orientation can be deformed by electric field when the dielectric tensor is orientation dependent. We then investigate the influence of a charged particle on the orientation order in a nematic state. The director is aligned either along or perpendicular to the local electric field around the charge, depending on the sign of the dielectric anisotropy. The deformation becomes stronger with increasing the ratio Ze/R, where Ze is the charge and R is the radius of the particle. Numerical analysis shows the presence of defects around the particle for large Ze/R. They are nanometer-scale defects for microscopic ions. If the dielectric anisotropy is positive, a Saturn ring defect appears. If it is negative, a pair of point defects appear apart from the particle surface, each being connected to the surface by a disclination line segment.

Structure of a liquid crystalline fluid around a macroparticle: Density functional theory

The structure of a molecular liquid, in both the nematic liquid crystalline and isotropic phases, around a cylindrical macroparticle, is studied using density functional theory. In the nematic phase the structure of the fluid is highly anisotropic with respect to the director, in agreement with results from simulation and phenomenological theories. On going into the isotropic phase the structure becomes rotationally invariant around the macroparticle with an oriented layer at the surface.

Measurement of elastic forces between iron colloidal particles in a nematic liquid crystal

The forces that arise between two iron particles in a nematic liquid crystal with a strong homeotropic anchoring were studied. For the first time, the short range repulsive force resulting from the presence of a hedgehog defect between two particles was precisely determined thanks to application of a small magnetic field and observation of the equilibrium position resulting from the balance between the elastic and magnetic forces. Above a given threshold force, the particles stuck together whereas the hedgehog defect was expelled and transformed into a Saturn ring located between the particles. The attractive part of the interparticle force was determined with the same method on the entire range of separation distances; we found that the equilibrium distance between two particles was r = 1.19 +/- 0.05 <d> (<d> was the average diameter of the pair of particles).

Anisotropic laser trapping in nematic colloidal dispersion

The interaction between a colloidal particle and a focused laser beam in a nematic liquid crystal reveals an unusual anisotropic Coulomb-like character. Experiments demonstrate two opposite directions in which the particle is attracted to and repelled from the nematic region deformed by the light-induced director reorientation. In this work we present analytical analysis of such behavior and derive the energy of interaction between colloidal particle and deformed director field. The analytical solution is in good agreement with recent results obtained by computer simulation.

Laser trapping of low refractive index colloids in a nematic liquid crystal

We describe and analyze laser trapping of small colloidal particles in a nematic liquid crystal, where the index of refraction of colloids is smaller compared to the indices of the liquid crystal. Two mechanisms are identified that are responsible for this anomalous trapping: (i) below the optical Fréedericksz transition, the trapping is due to the anisotropic dielectric interaction of the polarized light with the inhomogeneous director field around the colloid, (ii) above the optical Fréedericksz transition, the optical trapping is accompanied by the elasticity-mediated interaction between the optically distorted region of a liquid crystal and the colloid. In the majority of the experiments, the trapping above the Fréedericksz transition is highly anisotropic. Qualitative agreement is found with a numerical analysis, considering the nematic director elastic distortion, dielectric director-light field coupling and optical repulsion due to low refraction index colloid in high index surroundings.

Colloid-induced structure in liquid crystal media

The structural perturbations induced by colloidal particles immersed in a model nematic subjected to an external field are calculated employing integral equation methods. Maps of the density-orientational distribution about a colloidal particle are obtained, and these provide a microscopic picture of the colloid's nematic coat. We focus on colloidal particles that favor homeotropic anchoring, but planar anchoring cases are also considered. The range and structure of the nematic coat is shown to be significantly influenced by the nature of the anchoring, the size of the colloidal particle, the range and strength of the colloid-nematogen interaction, and the external field strength. All of these factors are discussed.

Colloidal structures from bulk demixing in liquid crystals

This experimental paper deals with phase separations of binary mixtures composed of a continuous liquid crystal phase and an isotropic dispersed phase. In contrast to isotropic binary mixtures, the investigated mixtures do not lead to a full phase separation but to a self-ordering of colloidal particles, as reported earlier (Loudet, J. C. et al. Nature 2000, 407, 611). We present here further aspects of such phase separations which include the kinetics of the phase separation, the origin of the formation of dislocation-like patterns, the influence of surfactants, chiral additives, and temperature on the formed colloidal structures. The present results show that (i) the dislocations in chain arrays can be seen as kinetically frozen defects, (ii) temperature can be used to control the size of the domains formed upon demixing, (iii) a slight change in surface chemistry, via the addition of surfactants, profoundly alters the formed colloidal structures, and (iv) chiral additives allow the formation of unique helical pearl chains which reflect the symmetry of the liquid crystal phase they are embedded in.

Interactions between spherical colloids mediated by a liquid crystal: A molecular simulation and mesoscale study

Monte Carlo simulations and dynamic field theory (DyFT) are used to study the interactions between dilute spherical particles, dispersed in nematic and isotropic phases of a liquid crystal. A recently developed simulation method (expanded ensemble density of states) was used to determine the potential of mean force (PMF) between the two spheres as a function of their separation and size. The PMF was also calculated by a dynamic field theory that describes the evolution of the local tensor order parameter. Both methods reveal an overall attraction between the colloids in the nematic phase; in the isotropic phase, the overall attraction between the colloids is much weaker, whereas the repulsion at short range is stronger. In addition, both methods predict a new topology of the disclination lines, which arises when the particles approach each other. The theory is found to describe the results of simulations remarkably well, down to length scales comparable to the size of the molecules. At separations corresponding to the width of individual molecular layers on the particles' surface, the two methods yield different defect structures. We attribute this difference to the neglect of density inhomogeneities in the DyFT. We also investigate the effects of the size of spherical colloids on their interactions.

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.

Nematic liquid crystal around a spherical particle: Investigation of the defect structure and its stability using adaptive mesh refinement

We investigate the orientation profile and the structure of topological defects of a nematic liquid crystal around a spherical particle using an adaptive mesh refinement scheme developed by us previously. The previous work [J. Fukuda et al., Phys. Rev. E 65, 041709 (2002)] was devoted to the investigation of the fine structure of a hyperbolic hedgehog defect that the particle accompanies and in this paper we present the equilibrium profile of the Saturn ring configuration. The radius of the Saturn ring r(d) in units of the particle radius R(0) increases weakly with the increase of Epsilon, the ratio of the nematic coherence length to R(0). Next we discuss the energetic stability of a hedgehog and a Saturn ring. The use of adaptive mesh refinement scheme together with a tensor orientational order parameter Q (alpha, beta) allows us to calculate the elastic energy of a nematic liquid crystal without any assumption of the structure and the energy of the defect core as in the previous similar studies. The reduced free energy of a nematic liquid crystal, F= F/L1RO, with L(1) being the elastic constant, is almost independent of Epsilon in the hedgehog configuration, while it shows a logarithmic dependence in the Saturn ring configuration. This result clearly indicates that the energetic stability of a hedgehog to a Saturn ring for a large particle is definitely attributed to the large defect energy of the Saturn ring with a large radius.

Orientational defects near colloidal particles in a nematic liquid crystal

We study the interaction between a surface-anchoring colloidal particle and a liquid-crystalline host, and in particular the formation of orientational defects near the particle. A mean-field theory based on the nonlocal Marrucci-Greco nematic potential is used to represent molecular interactions in an inhomogeneous orientational field. An evolution equation for the molecular configuration tensor is solved numerically whose steady state minimizes the total free energy of the system. With strong homeotropic anchoring on the particle surface, three types of solutions may appear depending on initial conditions and particle size: Saturn rings, satellite point defects, and polar rings. The Saturn ring remains stable on micrometer-sized particles, contrary to previous calculations but consistent with experiments. A phase diagram is constructed for the three regimes. Based on the free energy, the most stable state is the Saturn ring for smaller particles and the satellite defect for larger ones.

Defect structure around two colloids in a liquid crystal

This Letter investigates the defect structures that arise between two colloidal spheres immersed in a nematic liquid crystal. Molecular simulations and a dynamic field theory are employed to arrive at molecular-level and mesoscopic descriptions of the systems of interest. At large separations, each sphere is surrounded by a Saturn ring defect. However, at short separations both theory and simulation predict that a third disclination ring appears in between the spheres, in a plane normal to the Saturn rings. This feature gives rise to an effective binding of the particles. The structures predicted by field theory and molecular simulations are consistent with each other.

Forces between elongated particles in a nematic colloid

Using molecular dynamics simulations we study the interactions between elongated colloidal particles (length to breath ratio >>1) in a nematic host. The simulation results are compared to the results of a Landau-de Gennes elastic free energy. We find that depletion forces dominate for the sizes of the colloidal particles studied. The tangential component of the force, however, allows us to resolve the elastic contribution to the total interaction. We find that this contribution differs from the quadrupolar interaction predicted at large separations. The difference is due to the presence of nonlinear effects, namely, the change in the positions and structure of the defects and their annihilation at small separations.

Dynamic interaction between suspended particles and defects in a nematic liquid crystal

Insertion of spherical particles into a uniform nematic liquid crystal gives rise to the formation of topological defects. In the present work, we investigate how a spherical particle accompanied by its topological defects interacts with neighboring disclination lines. We perform two- and three-dimensional dynamic simulations to analyze the effect of a particle on the annihilation process of two disclination lines. The dynamics of the liquid crystal is described by a time-dependent evolution equation on the symmetric traceless order parameter that includes some of the salient features of liquid crystalline materials: excluded volume effects, or equivalently, short-range order elasticity and long-range order elasticity. At the surface of the particle, the liquid crystal is assumed to exhibit strong homeotropic anchoring. The particle is located between two disclination lines of topological charges +1/2 and -1/2. Two-dimensional simulations indicate that the topological defects bound to the particle mediate an interaction between the two disclination lines which increases the attraction between them. This result is confirmed by three-dimensional simulations that provide a complete description of the director field and of the order parameter around the particle. These simulations indicate that a spherical particle between two disclination lines can be surrounded by a Saturn ring, and suggest that the dynamic behavior of disclination lines could be used to report the structure of a defect around the particle.