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 di ﬀ erence 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 ﬂ uid 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.


Introduction
It is well-known that liquid-crystal molecules (mesogens) in the nematic phase exhibit no long-range positional order of their centers-of-mass but tend to align their longer axes parallel to a preferred direction specied by the global director n0 . 1 If a colloid is immersed in such a homogeneously ordered nematic phase the order is perturbed locally.This is because mesogens in the vicinity of the colloid's surface may orient themselves in such a way that their orientation is in contradiction with n0 . 2 As a consequence of the competition between n0 and the local orientation of mesogens near the surface of a colloid, defects arise in the liquid-crystal host phase.][4][5][6] Depending on the specic preparation of the colloid the preferred local orientation of mesogens may be either parallel or perpendicular with respect to the colloid's surface normal.These specic ways of orienting mesogens can be realized, for example, by varying the chemical composition of the colloid's surface. 7f more than a single colloidal particle is suspended in the liquid-crystal host phase, their specic defect topologies interact if the separation of the colloids is small enough.This gives rise to an effective interaction between the colloids which may be repulsive or attractive, such that in the latter case specically structured assemblies can be formed. 8For example, if a chemically homogeneous spherical colloid with locally homeotropic surface anchoring (i.e., anchoring of mesogens parallel with respect to the colloid's local surface normal) is placed in a nematic host phase the well known Saturn ring defect topology forms around the colloid's equator. 6If two such colloids come sufficiently close, the defect lines interact and can entangle the colloidal dimer in several ways. 9,10ne may measure these effective intercolloidal forces by trapping them in a certain conguration with optical tweezers which can then be used to move the colloids through the host uid in a controlled manner. 11More recently, optical tweezers have also been used for reconnecting defects of already linked colloids to achieve a more stable conguration. 12Furthermore, optical tweezers are important tools for analyzing complex defect structures. 13The possibility of formation and controlled manipulation of these structures eventually allows one to fabricate novel ordered arrangements of an assembly of colloidal particles which already opened a wide range of applications in optics and photonics. 14riginally, most research in the eld was focused on spherical colloids suspended in nematic liquid crystals but in today's research there is no limit regarding the colloid's shape.These shapes range, for example, from rod-like and polygonal shapes to star-or ring-like structures. 86][17][18] This is primarily because recent advances in chemical synthesis allow one to prepare inhomogeneous colloids in a controlled fashion even down to particle sizes in the nanometer range. 19Generally speaking, the surface of a Janus colloid is composed of antithetic materials.For example, one may envision a spherical Janus colloid where mesogens are anchored in a locally planar fashion at one hemisphere whereas the anchoring may be locally homeotropic at the other.Clearly, more complex defect topologies are expected if such a Janus colloid is immersed in a nematic host phase. 15,16he complexity of the colloidal suspension may be further enhanced by replacing the nematic host phase by another ordered but structurally more sophisticated one.An interesting candidate in this respect is the cholesteric phase forming if the mesogens are chiral.The cholesteric phase may be thought of as a nematic phase in which the director becomes local and rotates continuously around a distinguished axis.Hence, the director forms a helix characterized by the pitch p which is the distance to be travelled along the preferred axis if the director eld rotates by an angle of 2p.
Consequently, to characterize order in the cholesteric phase requires n0 to be replaced by n0 (z) if we take the z-axis as that axis around which the helix evolves.It is then clear that in the cholesteric phase n0 (z)$ê z ¼ 0 where êz is a unit vector pointing in the z-direction.Hence, in planes orthogonal to êz the cholesteric phase may be perceived as a stacked sequence of twodimensional nematics where n0 (z) is rotated in the x-y plane between neighboring nematic planes.
If a colloid is now immersed in such a cholesteric host phase it has been demonstrated 20,21 that the ratio of pitch and colloidal size inuences the formation of defect structures.Hence, compared with the nematic phase, where defect topologies are usually independent of the colloid's size, the ratio of p and the colloidal size enters as a new parameter in the formation of defect topologies in the cholesteric phase.Hence, more complex topologies are anticipated in a liquid-crystal host phase composed of chiral mesogens.This is even more so if two or more colloids are suspended in a cholesteric host phase.Because of the inherent twist of the cholesteric helix the defect lines have even more possibilities to entangle the colloids. 12,131][22] As far as computer simulations of a colloid suspended in a cholesteric phase are concerned, no study exists to date to the best of our knowledge.The reason is perhaps twofold.First, because p is generally quite large, fairly large systems need to be employed in molecular simulations which can quickly render such simulations computationally prohibitive.Second, most standard models for liquid crystals composed of chiral mesogens are rather costly from a computational perspective because of pronounced shape anisotropy of the mesogens.For example, Memmer et al. 23 suggested an interaction potential for a pair of chiral Gay-Berne molecules.They could only study fairly small systems comprising a few hundred to a few thousand molecules which casts some doubt on the signicance of the results obtained for the cholesteric 24 and even more so for the much richer structures characteristic of blue phases. 23n contrast, we have recently demonstrated 25 that large systems of up to 4 Â 10 4 chiral mesogens can be studied at moderate computational expense if one employs mesogens with a much smaller aspect ratio of about 1.30. 257][28][29][30][31] In the achiral Hess-Su model a pair of mesogens interact via a Lennard-Jones potential where the strength of the attraction depends on the orientation of the mesogens.Because the chiral version of this model reproduces cholesteric and even blue phases unambiguously 25 we feel that it is timely to apply it to investigate defect topologies arising if a colloid is immersed in a cholesteric phase.Our work is based upon Monte Carlo (MC) simulations in the isothermal-isobaric ensemble.
We organized our manuscript as follows.In Section. 2 we introduce our model system.Section 3 is given to key theoretical quantities and to a presentation of our results which we discuss in the concluding Section.4.

The liquid-crystal host phase
In this work we consider a liquid-crystal host phase in which a colloidal particle is immersed.The liquid crystal is composed of N mesogens interacting with each other in a pairwise additive fashion.The interactions are described by the uid-uid (ff) potential function where r ij ¼ r i À r j is the distance vector connecting the centersof-mass of mesogens i and j located at r i and r j , respectively, r ij ¼ |r ij |, and rij ¼ r ij /r ij .Unit vectors ûi and ûj specify the orientations of mesogens i and j in space.Hence, u ff is a Lennard-Jones potential where the attractive term has been modied to account for different orientations of the mesogens.Accordingly, s is the "diameter" of a spherical reference molecule and 3 ff is the depth of the attractive well.As in our previous work 25 we take the anisotropy function to be given by where P 2 ðxÞ ¼ 1 2 ð3x 2 À 1Þ is the second Legendre polynomial.
Terms proportional to 3 1 and 3 2 arise from an expansion of the orientation dependence of the interaction between two achiral mesogens in terms of rotational invariants. 32,33The pseudo-scalar proportional to 3 3 accounts for the chirality of the interaction between a pair of mesogens where the degree of chirality is controlled by the magnitude of 3 3 and the handedness by its sign.The achiral version of the interaction potential was originally suggested by Hess and Su 26 whereas the introduction of chirality by a pseudo-scalar follows in spirit an earlier suggestion by Memmer et al. for a Gay-Berne-like model potential. 23ence, the anisotropy function preserves the head-tail symmetry of the interaction between pairs of mesogens characteristic of many liquid-crystalline materials, that is the interaction potential in eqn ( 1) is invariant with respect to a change in sign of ûi and/or ûj . 1 The model has already been demonstrated to be capable of reproducing properties of liquid crystals in a variety of contexts ranging from the formation of nematic phases [27][28][29][30] to the development of defect topologies arising near the surface of spherical colloidal particles immersed in a nematic liquidcrystal host phase composed of achiral mesogens. 16In its chiral version, it has recently been shown to correctly reproduce cholesteric and blue phases. 25

Homogeneous and Janus colloid
In this work the N mesogens interact with a colloidal particle via a Yukawa-like potential given by assuming that the colloid is placed at the origin of the (Cartesian) coordinate system at the center of the simulation cell and where the subscript "fc" refers to the interaction between a uid molecule (i.e., mesogen) and the colloid.In eqn (3), r i ¼ |r i |, ri ¼ r i /r i , r 0 is the hard-core radius of the colloid, 3 fc determines the strength of the interaction, h À1 is the screening length of the attraction, and g is the anchoring function that we will discuss further below.Notice that because the colloid is centered at the origin of the coordinate system, ri points from the center-ofmass of the colloid to the center-of-mass of mesogen i. Parameters are introduced to guarantee that the minimum of the potential remains at r min ¼ r 0 + s and that the depth of the attractive well remains unchanged if one varies h at xed g.
In this work we also focus on Janus colloids with surfaces composed of antithetic materials.In the present case the conictive properties of the surface of a Janus colloid are different ways of anchoring a mesogen locally.Anchoring is described by the anchoring function 0 # g # 1 in eqn (3) which serves to discriminate energetically the desired (undesired) local orientation of a mesogen relative to the surface normal of the colloid.
More specically, we take the anchoring function to be given as 16 g(r i , ûi where the dimensionless parameter 0 # c # 1 has to be introduced to weaken the homeotropic local alignment relative to the parallel one because in the Hess-Su model, mesogens always exhibit a rather pronounced tendency to favor the former alignment over the latter.The tendency to align homeotropically at at surfaces with no specic surface anchoring is known for real liquid crystals 34 as well as for model systems i.e. for Gay-Berne 35 or Hess-Su 29 mesogens.This effect can be explained by the strong side-side attraction of two mesogens and has also been observed for strongly curved surfaces. 16Beacuse of the shape anisotropy homeotropic alignment at a surface enables a higher areal density of mesogens.Because of the stronger side-side attraction this is also reduces the congurational potential energy compared with a planar arrangement of mesogens at the surface plane.Even though the latter would be entropically favored, the energatically favored homeotropic alignment apparently wins.In eqn (7), (8a) where g k anchors mesogen i locally such that its orientation ûi is parallel to the surface of the colloid (i.e., perpendicular to ri ) whereas g t serves to align mesogen i in a locally homeotropic fashion (i.e., parallel to ri ).Associated with g k and g t are weighting functions where êx is a unit vector pointing along the (positive) x-axis of the space-xed Cartesian coordinate system.The function w k is largest at the North Pole of the colloid dened by ri $ê x ¼ +1 whereas w t assumes its maximum at the colloid's South Pole where ri $ê x ¼ À1.To be physically meaningful both weighting functions in eqn (9) have to be positive semidenite which is . Specic values of g k and g t control the sharpness of the transition between the antithetic surface parts of the colloid.
Notice also that our model is capable of describing a chemically homogeneous colloid where mesogens are aligned either parallel or homeotropically at each point on the colloid's surface.To that end one needs to set one of the two weighting functions in eqn (9) equal to one and the other one equal to zero depending on the desired anchoring scenario (see also Table 1).

Fixing the far-eld directors
To study experimentally defect topologies evolving around colloidal particles in an ordered liquid-crystal host phase, it is customary to place the colloid and its host between solid substrates. 36Through a specic anchoring scenario at the substrates, this setup allows one to x the (nonlocal) far-eld director n(k) 0 in regions of the liquid crystal that are unperturbed by the presence of the colloid.Superscript k indicates that the far-eld director at the lower substrate (k ¼ 1) may differ from its counterpart at the upper substrate (k ¼ 2).
In the vicinity of the colloid, the mismatch between the local alignment at the colloid's surface and the alignment of mesogens at the solid substrates will generally give rise to a local director eld n0 (r) where n0 (r) does not necessarily coincide with n(k) 0 .Therefore, we follow the experimental setup in spirit and place our liquid-crystalline host phase plus the dissolved colloid between plane parallel, structureless solid surfaces separated along the z-axis by a distance s z .This distance is always chosen large enough so that a direct interaction with the colloid is precluded and such that in the cholesteric host phase a helix of a certain pitch can be realized based upon earlier experience (see below for specic values of s z employed here). 25enceforth, we assume the substrates to be located at z ¼ Às z /2 (k ¼ 1) and s z /2 (k ¼ 2).The interaction between a mesogen and the substrate is described by the potential function where Dz (k) i h z i AE s z /2 and the sign is chosen depending on whether a mesogen interacts with the lower (k ¼ 1, +) or upper substrate (k ¼ 2, À), respectively.Analogous to eqn (1) and ( 3) the parameter 3 fs determines the strength of the interaction between a mesogen and the solid substrates.Throughout this work we maintain 3 fs /3 ff ¼ 5.0.
To x the far-eld directors in space we introduce the anchoring function where we assume each far-eld director to lie in the x-y plane and q (k) 0 is an angle associated with projections of n(k) 0 onto the x-and y-axes, respectively (see Section. 3.2 for specic choices of q (k) 0 ).
Fixing the in-plane unidirectional alignment at the surface is a widely used operational technique in the experimental setup. 6,13,22For example, a thin layer of polyimide on at glass substrates may be "rubbed unidirectionally to obtain an excellent planar alignment". 13In this spirit we are following the experimental approach by anchoring the mesogens in the plane of the surfaces.

Properties and numerical details
In this work we employ MC simulations carried out predominantly in a specialized isothermal-isobaric ensemble.As explained elsewhere 16 in this ensemble the thermodynamic state of the liquid crystal is specied by N, volume s z A, the ratio s x /s y (s a is the side length of the simulation cell in a-direction), and the transverse component tensor P where P aa (a ¼ x or y) is a diagonal element of P; offdiagonal elements of P vanish as we do not expose the liquid crystal to any shear strain.For xed values of these variables we generate a distribution of congurations in the 5N-dimensional space spanned by the 3N Cartesian coordinates of the N centers of mass and the 2N orientations of the N linear mesogens.This distribution is proportional to exp Numerically, this distribution can be realized by employing a standard Metropolis algorithm adapted to the specialized isothermal-isobaric ensemble. 37This algorithm proceeds in a sequence of two steps.Initially, one of the N mesogens is picked sequentially and it is decided with equal probability whether to displace its center-of-mass or to rotate it around one of the three axes of the space-xed Cartesian coordinate system.The specic axis is also picked randomly with equal probability for each of the axes.Once a displacement or a rotation has been attempted for all mesogens, one attempt is made to change the area A ¼ s x s y .The sequence of N displacement or rotation attempts plus the one attempt to change A constitutes an MC cycle.Our results are typically based upon MC simulations employing 5.0 Â 10 4 cycles for equilibration followed by 1.5-3.0Â 10 5 cycles during which ensemble averages are computed.
Throughout this work we use the customary dimensionless (i.e., "reduced") units, that is lengths are given in units of s, energy in units of 3 ff , temperature in units of 3 ff /k B , and pressure in units of 3 ff /s 3 .Based upon previous work 16,25 we x T ¼ 0.95 and P k ¼ 1.80 because under these thermodynamic conditions a Boojum and a Saturn ring defect topology is stable if a homogeneous colloid with respective planar and homeotropic surface anchoring is immersed in a nematic phase of the achiral version of our liquid crystal. 16Using instead a Janus colloid, a Boojum ring defect topology would form under the same thermodynamic conditions. 16ur systems typically accommodate 2 Â 10 4 to 7 Â 10 4 mesogens interacting with the various components of the model via eqn (1), (3), and (10) where we set 3 1 ¼ À3 2 /2 ¼ 0.04 throughout this work.Under these conditions the liquid crystal is in the nematic phase if 3 3 ¼ 0.00 and in the cholesteric phase if 3 3 ¼ À0.10. 25These are the only values of the chirality parameter considered in this work unless otherwise stated.
In selecting this value of the coupling constant 3 3 some caution is advisable.First, one wishes to make |3 3 | as large as possible in order to minimize the pitch of the cholesteric phase.This is because we focus on a high ratio of the colloid's size to p.One could, of course, maximize the ratio d 0 /p (d 0 ¼ 2r 0 ) by increasing the size of the colloid.However, this would not be sensible as it would require host phases comprising many more mesogens.Consequently, this would render the MC simulations computationally more demanding.Second, it is known 38 and also observed in our simulations (not shown here) that in the presence of a colloid a transition to a blue phase is favored even before this transition occurs in the absence of the colloid, i.e. in a bulk liquid crystal. 25Hence, to avoid the transition, an upper threshold exists for |3 3 | in our model which is slightly larger than the value chosen.
To ensure that effects caused by the presence of the colloidal particle vanish at the boundaries of the simulation cell we calculated the local density prole where hN(z)i is the average number of molecules located in a slab of volume Adz centred on z.Plots in Fig. 1 reveal that r(z) ¼ 0 in the immediate vicinity of the substrate because of the diverging repulsive mesogen-substrate interactions as z / Às z /2 (Dz (1) / 0) [see eqn (10)].As z increases, r(z) turns out to be damped oscillatory function reecting the formation of layers of mesogens.These layers become increasingly less pronounced as z increases because u (1) fs decays to zero as z becomes larger.It is particularly gratifying that despite this nontrivial structure r(z) with and without a colloid are nearly identically the same for |z| T 7 indicating that in this region the presence of the colloid does not affect the substrate-induced structure.As one move towards the center of the simulation cell (z ¼ 0) where the colloid is located a region exists where r(z) without a colloid is slightly larger than that in the presence of a colloid which is to be expected.
The key property analyzed in this work is the local director eld n0 (r) which we obtain as a numerical solution of the eigenvalue equation 16,39 where l(r) is the local eigenvalue associated with the local eigenvector n(r) and is the local alignment tensor.In eqn (14), angular brackets denote an average in the isothermal-isobaric ensemble, r is the local density, the operator 5 represents the tensor product, 1 is the unit tensor, and d is the Dirac d-function.Hence, Q is a real, symmetric, traceless, second-rank tensor which we represent by a 3 Â 3 matrix.We obtain three eigenvalues and -vectors numerically using Jacobi's method. 40As usual we take the largest eigenvalue l + (r) as the local nematic order parameter and its associated eigenvector as the local nematic director eld.In a similar fashion we obtain n0 from the nonlocal analogues of eqn ( 13) and ( 14).Topological defects, which are the main subject of this work are singularities of the local director eld n0 (r).However, singularities in the mathematical sense cannot arise in practice but are instead signalled by nite regions of lower order.In the subsequent sections we are particularly interested in identifying such regions characterized by small values of l + (r).To that end we discretize our simulation cell into small cubic boxes of side length ds ¼ 0.3 and for all visualizations, mark a box with a blue   dot if l + (r) # 0.25.Whereas the choice of this threshold is admittedly somewhat arbitrary, it turned out in practice to optimize the visibility of the rather complex three-dimensional variation of disclination lines.
In the cholesteric phase the characteristic helix evolves along the z-axis as shown in Fig. 3. Fixing s z ¼ 27, we obtain a stable helical structure comprising half a pitch in the z-direction of length p x 50 as can be inferred from Fig. 2. Because of the excluded volume near the substrate s z À 2dz x p/2 (see Fig. 1).The helix is stabilized through the anchoring of mesogens at the solid substrate where we choose q (1) 0 ¼ q (2) 0 ¼ 0 in eqn (11).The other case considered in this work is a system in which s z ¼ 40 such that about three quarters of an unstrained full pitch are accommodated.In this case, q (2) 0 ¼ 0 as before but q (1) 0 ¼ 90 .Other parameters of our model system are summarized in Table 1.
At this stage it seams noteworthy that the cholesteric helix can withstand quite a substantial strain upon compressing or stretching it p/4.For example, in a system in which an unstrained helix of total length p would form this helix can be compressed all the way down to a total length of p/4 by reducing s z before the helix is destroyed.Destruction of the helix is signalled by a decay of l + (z) to a value close to zero at z ¼ 0 in the middle between the substrates.Were the helix still intact, l + (z) would be constant across the pore at a relatively high value of about 0.4.
The problem with straining the cholesteric helix on account of an improperly chosen value of s z can in principle be avoided by replacing the monostable directional anchoring conditions (one easy axis) by degenerate planar conditions (innitely many easy axes) so that n(1) 0 and n(2) 0 are free to adjust themselves in the x-y plane so as to remove the spurious strain.In this case an unstrained helix of a total length corresponding to some fraction of p can be accommodated irrespective of s z .
However, in view of the Janus colloid, where the angle between its orientation and the direction of the local director eld enters as a new parameter, the monostable directional anchoring of mesogens at both substrates is to be preferred over the degenerate planar one.This is because in the later case any orientation of n(1) 0 and n(2) 0 on the unit circle in the x-y plane is energetically permissible as long as the angle between n(1) 0 and n (2)  0 remains xed such that at any given s z an unstrained helix of constant length can be accommodated.Unfortunately, this way one loses control over the orientation of the local director eld at the position of the Janus colloid.
To ensure that the cholesteric order at the boundaries is not perturbed by the presence of the colloid we compute the local director eld n0 (z) as a function of z.Plots of the components n x (z) and n y (z) of n0 (z) and of the local nematic order parameter l + (z) reveal that an unperturbed helix exists in the host phase (see Fig. 3).This helix is apparently affected only negligibly by the presence of the colloid at the simulation cell's boundaries as one concludes from Fig. 3.
With n0 (r) from MC we also compute the Frank free energy as a quantitative measure of deformation of the local director eld.In the one-constant approximation, the Frank free energy can be cast as 41 In eqn (15), the constant q 0 h 2p/p is the wave number of the cholesteric helix and K x 1.63 ff /s (ref.42) is a material-specic constant, for our model turns out to be of the same order of magnitude as that characteristic of the Gay-Berne model. 43oreover, assuming T ¼ 285 K and a particle size of s ¼ 3.4 Â 10 À10 m we obtain 3 ff z 3.9 Â 10 À21 J.With these numbers the elastic constant for our model turns out to be K z 19pN in SI units.Thus, this order of magnitude is realistic for a typical liquid crystal such as p-azoxyanisole. 44In ref. 42 it was also shown that the one-constant approximation is valid here on account of the small aspect ratio of the mesogens.As one can see from eqn (15) q 0 modies the twist contribution to F .In the nematic phase, where q 0 ¼ 0 because p / N, eqn (15) reduces to the standard form 1 such that F is composed of the usual splay, twist, and bend contributions.
In general, the contribution of each term in eqn ( 15) is weighted separately by material-specic constants K 1 , K 2 , and K 3 which can be attributed to bend, splay, and twist deformations of the director eld, respectively.However, the one-constant approximation K 1 x K 2 x K 3 ¼ K is very well satised for the Hess-Su model. 42This is because of the relatively small aspect ratio of 1.30 of mesogens in the Hess-Su model.In practice, we compute F numerically by differentiating n(r) according to eqn (15).
At this stage it seems noteworthy that the approach taken here is different from the one usually applied in that we obtain n0 (r) as an ensemble average corresponding either to a state at thermodynamic equilibrium or to a rather stable, "long lived" metastable state.That is we obtain n0 (r) as result of intermolecular interactions.The conventional procedure is to take a mean-eld expression for the free energy and minimize this functional with respect to n0 (r) in principle.Both approaches should provide consistent results as long as the mean-eld approximation remains valid.

Homogeneous colloid in a cholesteric host phase
To validate the model system as being sufficiently realistic we begin our presentation of results by considering a chemically homogeneous colloid suspended in a cholesteric host phase.In this case the resulting defect topology is well known both theoretically using lattice Boltzmann calculations 20,21 and experimentally. 13he cholesteric phase is characterized by a helix evolving along the z-axis.As shown in Fig. 3 the unperturbed helix (i.e., sufficiently far away from the surface of the colloid) is characterized by a local director eld that can be tted almost perfectly by the expression To accommodate a helix of the desired length between the substrates we introduce the phase shi 4 0 which can be obtained by simply tting the MC data in Fig. 3 with eqn (16).Alternatively, 4 0 can be dertermined as follows.
Because of the form of the uid-substrate potential function given in eqn (10) the minimum of that potential is located at z min ¼ AEs z /2 H dz depending on whether the upper or lower substrate is considered.To accommodate a helix of pitch p between substrates separated by a distance s z we need to have where q (k) 0 (k ¼ 1 or k ¼ 2) is the anchoring angle either at the lower (k ¼ 1, +) or at the upper substrate (k ¼ 2, À).For our present choice 3 3 ¼ À0.10, we determine p x 50. 25 With these values and a given anchoring angle either at the lower or at the upper substrate, the phase shi 4 0 can easily be determined.
The mesogens are anchored directionally at the lower substrate such that their longer axes point along the y-axis, and hence q in eqn (11).Under these conditions 4 0 z 5 4 p is obtained from eqn (17) for s z ¼ 40.Because we have chosen both p and s z such that three quarters of a full pitch can be accommodated without exposing the helix to any spurious strain, we also anchor our mesogens at the upper substrate directionally such that their longer axes now point along the x-axis which is realized by setting q (2) 0 ¼ 0 in eqn ( 11) leading, of course, to the same 4 0 from eqn (17).
We begin our discussion of results by considering a chemically homogeneous colloid with r 0 ¼ 15 (d 0 /p ¼ 0.6) which anchors the mesogens in a locally homeotropic fashion.If this colloid were immersed in a nematic rather than a cholesteric host phase the well-known Saturn ring defect topology would evolve. 16It consists of a closed disclination line encircling the colloid in the vicinity of its equator if the line connecting the colloid's North and South Pole is parallel to the far-eld director.
Because of a possible mismatch between n0 (z) in a cholesteric phase and the local orientation of the mesogens at the surface of the colloid, defects also arise in this case.However, in the present case their topology is much more complex than for a nematic host uid.For the visualization of this topology we adopt the procedure already described in Section.3.1.The plot in Fig. 4 shows resulting disclination lines entangling a colloid suspended in a cholesteric host phase.The change of n0 (z) along the z-axis is indicated by the small double arrows rotating around a line parallel to the z-axis in the x-y plane.The resulting disclination line is closed and wrapped around the colloid in a rather complex fashion.In particular, the plot in Fig. 4 agrees very nicely with earlier lattice Boltzmann studies 21 and experimental results 13 and thus validates the model used here.
To rationalize the complex spatial variation of the disclination line one may invoke a simple geometrical argument.Consider a vector r pointing from the center of the colloid in an outward direction such that locally r is collinear to the orientation of a mesogen that is homeotropically aligned.Hence, for an ideal defect to arise, r needs to be orthogonal to the local director eld satisfying the expression where y and z are restricted to the interval [Àr eff , r eff ].The effective radius r eff arises because depending on the anchoring strength, the presence of repulsive mesogen-colloid interaction, and details of the thermodynamic conditions, disclination lines will in general be removed to some extent from the colloids's surface, thus r eff > r 0 .In practice, we shall take r eff as an adjustable parameter that we choose to match the MC data to the maximum extent possible.Eqn (18) can easily be solved and one obtains the disclination-line path in parametric form.Notice, that only the "++" and "ÀÀ" combinations of the vector components of r t are solutions of eqn (18).By tting our MC data we obtain r eff x 20.With this value the red line in Fig. 4 represents the disclination-line path calculated from eqn (19) indicating that the above reasoning is well justied.It is also interesting to note that larger colloids are encircled by disclination lines with more than single loop.This is effected in eqn ( 19) by increasing r eff .Disclination lines with more than a single loop have also been observed in lattice Boltzmann simulations presented in ref. 20.However, if the spherical, chemically homogeneous colloid anchors the mesogens in a locally planar fashion, the defect topology changes dramatically.Instead of one closed disclination loop as in Fig. 4 we observe two individual disclination lines of nite length wrapped around the colloid as shown in Fig. 5. Again, our results are in accord with earlier lattice Boltzmann results. 20he geometric concept behind eqn (18) can also be applied to the case of planar anchoring but requires additional modi-cations.In this case, the ideal defect arises if the preferred, planar anchoring is perpendicular to the local director, thus n0 (z) is parallel to r and satises the expression The solution of this expression leads to two defect points due to the fact that n0 (z) is rotating in the x-y plane and can only be parallel to r if its z-components vanish.For that reason one expects the formation of these defect points at z ¼ 0.
However, our simulation data reveal two line defects instead.In this case, one may wonder how the plot in Fig. 5 depends on our criterion to identify disclination lines.To that end we varied our threshold.It turned out that the length of the disclination lines depicted in Fig. 5 is independent of the threshold used for l + (r); however, the thickness of the blue lines vanishes completely (and thus the disclination lines) if one uses l + (r) # 0.05 to identify regions of lower order in the cholesteric host phase.Because the threshold value for l + (r) does not affect the length of the disclination lines shown in Fig. 5, we propose the following line of arguments to unravel the apparent discrepancy between MC data and the point defects that would form ideally.
In general, the formation of topological defects (which in a mathematical sense are singularities in the director-eld) costs free energy.However, it is conceivable that in comparison with the formation of point defects, a strong local deformation of the (continuous) director eld is even more costly in terms of free energy.If that is so, the formation of defects is advantageous and entire lines of defects may form.Eventually, however, the deformation of the director eld wins energetically such that the disclination lines must end and therefore be of nite length.
Based upon this qualitative energetic argument one would also argue that a disclination line forms as long as the mismatch between n0 (z) and the local alignment of mesogens at the colloid's surface is maximized.Mathematically speaking, the disclination lines satisfy for each xed z.Differentiating the right side of eqn ( 18) and solving eqn (21) one easily veries that the disclination line path is given in this case by Notice though that r k (z) in eqn ( 22) again turns out to encircle the colloid in a closed loop as compared to eqn (19), only sine and cosine functions are interchanged.However, to rationalize that disclination lines are of nite extent in reality, we argue that the largest mismatch between n0 (z) and the locally planar alignment of mesogens at the surface of the colloid arise at z ¼ 0. The mismatch then decreases with increasing |z|, vanishing eventually at |z| ¼ r eff .However, there has to exist a certain threshold value n0 (z)$r k ¼ dr for the formation of topological defects.Beyond this value an elastic deformation of the director eld is energetically favored compared to the defect formation.Clearly, the value of the threshold is a priori unknown but may be adjusted like r eff to obtain an optimal representation of the MC data.This representation is possible for r eff x 22 and dr x 0.8r eff as the plot in Fig. 5 illustrates where again the "++" and "ÀÀ" combinations of the components of r k in the interval [Àr eff , r eff ] have been used.
The similarity between eqn ( 19) and ( 22) can be rationalized as follows.Consider a cross section of the colloid in the x-y plane that is located at some z ¼ z 0 .With respect to the local director n0 (z 0 ) but irrespective of both the actual anchoring of mesogens at the surface of the colloid and the direction of n0 (z 0 ), four distinguished points exist (see Fig. 6).At two of them n0 (z 0 ) is locally parallel to the cross section whereas at the other two n0 (z 0 ) is locally orthogonal to it.Hence, at the latter two one anticipates the formation of defects.The position of defects alternates with position of the other two points at which the local orientation of the cross section agrees perfectly with n0 (z 0 ).Because of eqn ( 16) all four points rotate in the x-y plane in a counter clockwise fashion with increasing z.

The twisted nematic cell
In the preceding section an intrinsic cholesteric host phase could form because of the chirality of the interactions between a pair of mesogens.However, a quasi-cholesteric phase may be induced if one places achiral mesogens (3 3 ¼ 0) between planar solid surfaces with hybrid anchoring for which the energetically most favorable alignment of mesogens at one substrate differs from that at the other one.
Consider rst a homogeneous anchoring scenario realized by setting q (1) ¼ q (2) in eqn (11), such that n(1) 0 ¼ n(2) 0 .If thermodynamic conditions have been chosen properly, a nematic Fig. 5 As Fig. 4, but for a planar anchoring.The fit of eqn (22) to the MC data is obtained for r eff ¼ 22 and dr ¼ 0.2r eff .
phase will form between the substrates characterized by n(1) If one now maintains q (2) 0 but changes q (1) 0 the mesogens have to change their orientation from the one imposed by the lower substrate to the favorable one at the upper substrate.This is possible by gradually switching the local director eld from n0 (Às z /2) ¼ n (1)  0 to n0 (s z /2) ¼ n(2) 0 along the z-axis normal to the substrates' planes.As a result a helical structure evolves along the z-axis similar to the one formed when the mesogens themselves are chiral.
However, in the present case the quasi-cholesteric phase is induced solely by the substrates: The important difference is that in the case of intrinsic chirality, the pitch of the cholesteric helix can be modied by changing 3 3 in the mesogen-mesogen interaction potential.A similar tuning of interactions is, of course, impossible in real materials.Here the pitch is a quantity characteristic of any given material.In the quasi-cholesteric phase, however, p is a function of the angle increment Dq 0 h q (1) 0 À q (2) 0 and the handedness of the helix is determined by the sign of Dq 0 (see Fig. 7).Notice, in particular, the different pitch in the quasi-cholesteric helix (see Fig. 7) compared with typical values of p in an intrinsically cholesteric liquid crystal (see Fig. 2).
The formation of a quasi-cholesteric, substrate-induced phase can also be realized experimentally in a so called twisted nematic cell (TNC).However, experimentally it is generally impossible to gradually vary the anchoring at one substrate relative to the other.In an experimental setup one would simply employ homogeneously anchoring substrates where one of them is rotated in the x-y plane relative to the other.
If a chemically homogeneous colloid is immersed in a quasicholesteric phase, one anticipates the formation of disclination lines.Starting with a colloid anchoring mesogens in a locally homeotropic fashion, we see in Fig. 8(a) that the well-known Saturn ring defect topology arises for the case q (1) 0 ¼ q (2) 0 ¼ 0 where the host phase is nematic with n0 pointing in the xdirection.Because of the direction of n0 the maximum mismatch between n0 and the local alignment of mesogens at the colloidal surface arises around the equator of the colloid as one can see from Fig. 8(a).
Increasing now q (1) 0 in steps of 10 sequentially from one MC run to the next one realizes from Fig. 7 that p decreases from innity (nematic phase) to a large but nite value in the quasicholesteric phase.Plots in Fig. 8(b)-(f) show that in the quasicholesteric phase the initial Saturn ring [see Fig. 8(a)] becomes increasingly deformed as p decays (i.e., as q (1) 0 increases).However, it is particularly noteworthy from the plots in Fig. 8 that both undeformed and deformed disclination lines can be represented by the simple analytic expression given in eqn (19).The projection of disclination lines onto the x-y plane in Fig. 8(f) in conjunction with its full three-dimensional representation in Fig. 8(e) reveals that the disclination line is a closed loop wrapped around the colloid for q (1) 0 ¼ 110 .At this point another feature of a quasi-cholesteric helix in a TNC setup should perhaps be emphasized.For q (1) 0 > 90 the helix has two options of responding to a further increase of q (1) 0 .On the one hand, p may decrease further such that the already metastable quasi-cholesteric helix becomes more twisted or, on the other hand, the handedness of the helix may ip at the expense of an increase of p.The latter is possible because of the head-tail symmetry of the achiral mesogens in our model.In practice, we observe in our MC simulations that the change in helix handedness and a simultaneous sudden increase of p occur under TNC conditions if q (1) 0 > 110 .At this stage another important difference between a quasicholesteric phase in a TNC setup and a true cholesteric phase in a liquid crystal with intrinsic mesogenic chirality arises.Starting with a quasi-cholesteric and a cholesteric phase of the same p we are able to reduce p continuously in the latter case by increasing 3 3 .At some critical value of 3 3 the cholesteric phase undergoes a phase transition to a blue phase as we demonstrated elsewhere. 25This blue phase is distinguished from the cholesteric phase structurally because it exhibits doubleinstead of single-twist helices.A similar phase transition is impossible in the TNC with a liquid crystal composed of achiral mesogens.
Turning now to a discussion of a TNC where a colloid with planar local surface anchoring is immersed in a quasi-cholesteric host phase, we investigate in Fig. 9 the effect of gradually decreasing p as q (1) 0 increases.Starting again with the nematic phase of innite p (e.g., q (1) 0 ¼ q (2) 0 ¼ 0 ) we see from Fig. 9(a) that the well-known Boojum defect topology forms where two defects are located at the North and South Pole of the colloid.Fig. 7 The red line represents the pitch of a TNC with s z ¼ 27 as a function of the angle q (1) 0 at constant q (2) 0 ¼ 0 (see eqn 11) whereas the circles mark the relevant pitches used in Fig. 8 and 9 in a TNC.
Again, as q (1) 0 increases from its initial value of 0 the localized defects of the Boojum topology become short lines which are stretched the larger q (1) 0 (the smaller p) becomes [see Fig. 9(b)-(e)].For q (1) 0 ¼ 120 the three-dimensional representation in Fig. 9(e) and its projection onto the x-y plane in Fig. 9(f) show that the disclination lines are partially wrapped around the colloid but remain disconnected.This is apropos of the situation depicted in Fig. 5 for a colloid suspended in a cholesteric phase composed of intrinsically chiral mesogens.However, compared with the plot in Fig. 5 the disclination lines depicted in Fig. 9(f) are shorter, which is primarily ascribed to the smaller size of the colloid in the latter case.However, the reader should realize that it is not the absolute value of the colloid's size that matters but only the ratio between this size and the pitch.Typical values for this ratio in the present case [see Fig. 8(f)] are 0.16 as opposed to, for instance, 0.6 for data presented in Fig. 4 or 5.
For anchoring angles q (1) 0 > 120 the same sudden change in handedness of the already metastable quasi-cholesteric helix Fig. 8 As Fig. 4, but for a spherical colloid of radius r 0 ¼ 6 suspended in a TNC.The angle q in the top-right corner of the plots (a)-(f) indicates the rotation of the directional anchoring vector in eqn (11) and is used for calculation of the pitch.The ratio d 0 /p in parts (a)-(e) is 0.00, 0.04, 0.08, 0.12, and 0.15 respectively.and increase of p occurs that was already described above for a colloid with locally homeotropic surface anchoring of the mesogens.Nevertheless, the plots presented in Fig. 9(a)-(f) also show that the simple analytic expression for disclination lines given in eqn ( 22) is capable of representing the MC data very well.Hence we conclude that as far as disclination lines are concerned there is essentially no difference between a cholesteric phase formed by intrinsically chiral mesogens and a quasi-cholesteric phase composed of achiral mesogens in a TNC.

Janus colloid
If, instead of a colloid with homogeneous surface anchoring, a Janus colloid is immersed in a cholesteric phase, the already quite complex spatial variation of disclination lines becomes even more complicated.This is because in the case of a Janus colloid the orientation of its antithetic surfaces with respect to the local orientation of the director eld at the position of the colloid's center of mass enters the analysis as a new parameter w which we introduce via n0 (0)$û J ¼ cos w where ûJ is a unit vector describing the orientation of the Janus colloid.Hence, we focus on the rotation of the Janus colloid around the z-axis, such that ûJ $ê z ¼ 0.
For all the results presented below we x the position of the Janus colloid such that its center of mass coincides with the origin of the Cartesian coordinate system located at the center of the simulation cell.To analyze the effect of different orientations of the Janus colloid with respect to the cholesteric helix we proceed as follows.
We x ûJ such that is parallel to the x-axis.By changing the angles q (k) 0 (k ¼ 1, 2) in eqn (11) simultaneously one can then expose the Janus colloid to different orientations of the director eld.Therefore, disclination lines form over different parts of the Janus colloid as w is varied.Moreover, because we x the substrate distance s z ¼ 27 and maintain the coupling constant 3 3 ¼ À0.10 the pitch is always xed such that half a pitch ts between the substrates without spurious strain.
Results plotted in Fig. 10 illustrate the rather complex spatial variation of disclination lines as w varies between 0 and 90 .In all cases we have only a single line of disclinations with disconnected ends such that the line neither forms a closed loop as in the case of a chemically homogeneous colloid with homeotropic surface anchoring nor form two disconnected lines as is the case for locally planar surface anchoring.However, as before in the case of a chemically homogeneous colloid with planar anchoring the disclination line we are dealing with here is of nite length.The fact that we are dealing only with a single disclination line emphasizes the interplay between the locally planar and homeotropic surface anchoring over two hemispheres of the Janus colloid.Generally speaking, the plots in Fig. 10 This variation of disclination lines with w also has another interesting consequence which becomes apparent as one changes the size of the Janus colloid.As it is mainly n0 (r) that changes as w is varied, it seems plausible to discriminate energetically between states for different w via the Frank free energy which is a quantitative measure of the free energy density associated with different defect topologies.Based upon a quantitative analysis of the Frank free energy [see eqn (15)] plots in Fig. 11 illustrate the effect.Starting with a small colloid of hard-core radius r 0 ¼ 3 (d 0 /p ¼ 0.12) we see in Fig. 11(a) that the Frank free energy is apparently insensitive to a variation of w except the weak minimum at w ¼ 0. Comparing Fig. 10(a) and 12(a) reveals the size independence of the defect structure at w ¼ 90 as far as the overall topology of the disclination line is concerned.This changes dramatically at w ¼ 0 .The disclination line in Fig. 12(b) is much closer to the surface and is less curved than the one shown in Fig. 10(b).These differences between disclination lines presented in Fig. 10(b) and 12(b) might be the reason for the different behavior of F (w).For a vanishing ratio d 0 /p (p / N, nematic host phase) one anticipates the formation of a Boojum ring defect topology as conrmed by plots in Fig. 10(c) and in agreement with experiments, Landau-de Gennes theory, and molecular simulations. 7,15,16Thus, the defect topology in a cholesteric phase, where d 0 /p ¼ 0.12 [see Fig. 12(b)] is intermediate between the one shown in Fig. 10(b) (d 0 /p ¼ 0.24) and the Boojum ring topology shown Fig. 12(c).
Returning now to a more detailed analysis of the Frank free energy in Fig. 11, we notice that for d 0 /p ¼ 0.16 the plot of F in Fig. 11(a) reveals that F becomes a monotonically decaying function of w with the largest values arising at small angles.However, F quickly levels off as w increases beyond 30 .Going nally to a Janus colloid of r 0 ¼ 6 (d 0 /p ¼ 0.24) we notice that the plot of F in Fig. 11(a) exhibits a maximum at w ¼ 0 and then decays monotonically until a constant value is assumed.The different values in the plateaus of all three curves in Fig. 11(a) are related to the increase in the length of the disclination lines as the colloid becomes bigger.Dividing the curves plotted in Fig. 11(a) by the respective plateau value F 0 reveals that the width of the free energy barrier as well as the relative height of this barrier increases the larger the colloid becomes.This can be seen in Fig. 11(b) where the plot of F /F 0 approaches the common value of 1 at increasingly larger w as the radius of the colloid increases.The same results are expected for a xed size of the Janus colloid but increasingly shorter pitch.The crucial quantity here is, again, the between the colloid's size and the pitch and not individual values of these quantities.
Because of the periodicity of the cholesteric helix the results shown in Fig. 11 will also be repeated periodically.For example, for the curves shown in Fig. 11 there will be another free-energy barrier at angles approaching 180 .These observations offer the possibility for the following Gedankenexperiment.Suppose, one were able to x the orientation of a Janus colloid which might be conceivable by using a homogeneous magnetic eld. 45For Janus colloids much bigger then the ones that could be studied here on account of computational limitations most angles would correspond to energetically unfavorable regions because the width of the free-energy barrier increases with d 0 /p.Hence, the Janus colloid would be more or less restricted to positions somewhere in the x-y plane so that it might be feasible to form a quasi two-dimensional layer of Janus colloids all with the same orientation through orientational self-assembly.

Conclusions
In this work we employ MC simulations carried out in the isothermal-isobaric ensemble to calculate the local director eld of a cholesteric liquid-crystal host phase in the presence of an immersed spherical particle.We employ a simple, modied Lennard-Jones model for the liquid crystal following earlier suggestions by Hess and Su. 26 To form a cholesteric phase the mesogens need to be chiral.We introduce chirality through a pseudo-scalar coupling in addition to the other orientation dependent terms in the Hess-Su potential.This follows in spirit a suggestion by Memmer et al. for the interaction between a mesogenic pair in the classical Gay-Berne model. 23y considering homogeneous colloids with either locally planar or homeotropic surface anchoring of the mesogens we determine disclination lines that had already been obtained earlier in either experiments or in lattice-Boltzmann calculations. 13,20,21Besides the validation of our model through this agreement with earlier works we also offer a simple geometrical argument that explains the complex spatial variation of disclination lines quite nicely.
Based upon the present study one might consider a system with two colloids immersed in a cholesteric host to study the  effective pair-wise interactions.This has already been done in experiments 22 and lattice-Boltzmann calculations 46,47 and could, in principle, also be realized in our model system.In another work we disperse two homogeneous colloids with planar anchoring in a nematic phase and obtain the Frank free energy as a function of the distance between two colloids and the angle formed by the distance vector connecting their centers and the nematic director. 48It turns out that the Frank free energy is wellsuited to study the effective interactions between a colloidal pair that are mediated by a nematic host phase.The Frank free energy allowed us 48 to explain quantitatively the "magic angle" of about 30 formed between the center-of-mass distance vector of the colloidal pair and n0 which is obtained experientially 6,11 for a stable conguration of the colloids.From this agreement we conclude that F (w) is equally well-suited to determine the stability of a Janus colloid immersed in a cholesteric host phase.
However, to apply the approach of ref.48 to study effective interactions between a pair of Janus colloids in such a cholesteric phase comprehensively would be much more demanding.Here the Frank free energy is a function of all three coordinates (e.g.spherical coordinates) and the ratio of d 0 /p.In the case of Janus colloids one has to take into account the relative orientation of the Janus colloid with respect to the cholesteric helix in addition.Given large and sufficient computational resources one can, of course, use this model system to study a colloidal interactions mediated by the cholesteric phase.However, at least at present the dimension of the parameter space corresponding to such a setup seems prohibitive.
In addition to the Hess-Su potential, where the intrinsic chirality is introduced via the interaction potential, we also study a quasi-cholesteric phase where the mesogen-mesogen interaction is achiral and the formation of a cholesteric helix is introduced through the anchoring of mesogens at two planar solid substrates.This setup follows in spirit that characteristic of a TNC.By comparing disclination lines arising in the intrinsically cholesteric and the quasi-cholesteric liquid crystal, we demonstrate that the topology of the disclination lines apparently only depends on the presence of a helically varying director eld irrespective of the precise origin of the helix.Again, the analytic expression derived here from a geometrical argument is capable of describing the complex topology of the resulting disclination lines.
Although the topology of disclination lines evolving around a spherical colloid with homogeneous surface anchoring is quite complex, even more complicated structures occur in the vicinity of a Janus colloid consisting of two hemispheres at which the mesogens are anchored in a locally planar and homeotropic fashion, respectively.Because of the presence of the two hemispheres of antithetic surface anchoring, the orientation of the Janus colloid relative to the local director eld of the cholesteric phase matters.Depending on this orientation the topology of disclination lines varies.
Associated with these different topologies is a higher (or lower) free energy which we determine by calculating the Frank free energy in the one-constant approximation.The oneconstant approximation is nearly exact because of the small aspect ratio of mesogens in the Hess-Su model. 26 varying the size of the colloid or the pitch of the cholesteric helix, we show that barriers in free energy arise for certain orientations of the Janus colloid relative to the local director eld.These barriers increase in magnitude with growing size of the Janus colloids.Simultaneously, the range of angles over which the free-energy barriers are observed becomes wider as the colloid becomes larger (or the pitch smaller).As the size of the colloid as well as the pitch are controllable parameters in corresponding experiments, 7,13,49 it seems conceivable that similar free-energy barriers might be used to form self-assembled quasi two-dimensional layers of Janus colloids experimentally, provided the orientation of the Janus colloid can be maintained.The latter may be possible by using homogeneous magnetic elds such that orientationally ordered quasi twodimensional arrays of Janus colloids could be formed that might have interesting properties.However, at this stage more in-depth research is clearly required to investigate the potential usefulness of the effect reported here for future applications.
Finally, it seems worthwhile emphasizing that the disclination lines obtained in this work have been computed as ensemble averages from relatively long Markov chains of congurations.In none of these individual simulations we observed "jumps" between structurally or thermodynamically distinct states similar to what is occasionally encountered in MC simulations of coexisting phases near the coexistence line.Hence, we conclude that the data presented here correspond to thermodynamically stable or metastable states of a lifetime exceeding that of the simulations considerably.

Fig. 1
Fig. 1 Plot of the local density r(z) as a function of position z in the lower part (z # 0) of our system; ( ) without and () with a colloidal particle.

Fig. 3
Fig. 3 Plots of components n x ( ) and n y ( ) of the local director field and of the local nematic order parameter l + ( ) as functions of z.Data are obtained for 3 3 ¼ À0.10 and s z ¼ 27, where lines refer to the presence and symbols to the absence of the colloid.

Fig. 2
Fig. 2 Pitch p as a function of the magnitude of the chiral coupling constant, |3 3 |.The red line is intended to guide the eye.The transition from the cholesteric to the blue phase occurs at |3 3 | ¼ 0.2. 25

Fig. 4
Fig. 4 Disclination line represented by a set of blue dots (see text) in the vicinity of a spherical colloidal particle of radius r 0 ¼ 15 (d 0 /p ¼ 0.6) suspended in a cholesteric host phase.Mesogens are anchored in a locally homeotropic fashion at each point of the colloid's surface.The red line is the fit of the expression in eqn (19) to the MC data taking r eff x 20 as a fit parameter.The little double arrows indicate the local director and the coloured plane visualizes their continuous variation along the z axis.The colors green and red indicate one or the other side of this plane.

Fig. 6
Fig. 6 (a) Sketch of a cross-section of the simulation cell in the x-y plane at z ¼ z 0 .Colloid and n0 (z) are represented by the gray shaded area and the dashes respectively.Red and green areas refer to respective regions in which locally planar anchored mesogens are orthogonal or collinear to n0 (z).(b) Corresponding data for n0 (z) (dashes) and the local order parameter at z ¼ 0.2 (see attached color bar).
(a) and (b) show that the way in which the disclination line is wrapped around the Janus colloid changes signicantly with w.

Fig. 10
Fig. 10 As Fig. 4, but for a Janus colloid (r 0 ¼ 6, d 0 /p ¼ 0.24) with planar surface anchoring at the upper (black) hemisphere and a homeotropic anchoring at the lower (grey) hemisphere.The arrow through the colloid indicates its orientation ûJ .The configurations are related to (a) w ¼ 90 and (b) w ¼ 0 .

Fig. 12
Fig.12 As Fig.10, but for r 0 ¼ 3, (d 0 /p ¼ 0.12).(a) w ¼ 90 , (b) w ¼ 0 , and (c) Boojum ring defect topology forming in a nematic host fluid (d 0 /p ¼ 0.).Notice that only a fraction of the entire system is shown.Therefore, the local director fields at the upper and lower boundary of the z-axis are not orthogonal as they would be at z ¼ AEs z /2 on account of the anchoring at the substrates.

Table 1
16del parameters used to realize specific defect topologies near a spherical colloid immersed in a nematic host phase composed of achiral mesogens (see Section. 2.2)16