Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 16814.
Published online 2017 December 1. doi:  10.1038/s41598-017-16967-1
PMCID: PMC5711923

Uncovering different states of topological defects in schlieren textures of a nematic liquid crystal


Topological defects are ubiquitously found in physical systems and therefore have been an important research subject of not only condensed matter physics but also cosmology. However, their fine structures remain elusive because of the microscopic scales involved. In the case of a liquid crystal, optical microscopy, although routinely used for the identification of liquid crystal phases and associated defects, does not have resolution high enough to distinguish fine structures of topological defects. Here we show that polarised and fluorescence microscopy, with the aid of numerical calculations on the orientational order and resulting image distortions, can uncover the structural states of topological defects with strength m =  ±1 in a thin cell of a nematic liquid crystal. Particularly, defects with m = +1 exhibit four different states arising from chiral symmetry breaking and up-down symmetry breaking. Our results demonstrate that optical microscopy is still a powerful tool to identify fine states of liquid crystalline defects.


Liquid crystals (LCs)15 have attracted interests of physicists as ideal model systems which facilitate direct observation of the structures of topological defects by optical means; for example, LCs have been studied as a tabletop model system of cosmic strings6,7. Not only from such a fundamental point of view, LCs defects have been drawing attention also as tuneable or reconfigurable templates for nano/micro-patterning of small non-LC objects, e.g., micro particles815, polymers16,17, self-assembled molecules18,19 and low molecular weight solutes20,21. Among different phases of LCs, the simplest nontrivial phase is the nematic phase characterised by broken rotational symmetry, while retaining the translational symmetry. Its local orientational order is specified by a unit vector without head-tail distinction, known as the director commonly denoted by n (ref.1).

The head-tail symmetry of n reflects the apolar nature of the nematic liquid crystal (NLC), and thus the order parameter space of a two-dimensional (2D) NLC is S 1/Z 2, a circle (S 1) whose two ends of a diameter are identical. Topological structures of the order parameter space and associated defects are described by the homotopy groups of the order parameter space. The first homotopy group of S 1/Z 2 is Z (refs3,5), that is, defects in a 2D NLC with different “winding number” m are topologically distinct. Moreover, m can be a half-integer as well as an integer because of the head-tail symmetry (Fig. 1b). The easiest way to determine m is an observation using a polarised optical microscope (POM). Under a POM, a NLC in a flat cell imposing planar surface alignment with no preferential direction (degenerate planar alignment) exhibits schlieren textures22,23 with topological defects from which dark brushes of even number (k) emanate (Fig. 1a). In a dark brush, n is parallel or perpendicular to the polariser (or analyser) axis, and the winding number m is related to k simply by k = 4|m|. The sign of m can be distinguished by the rotation of the crossed polarisers; dark brushes from a positive (negative) defect rotate in the direction the same as (opposite to) that of the polarisers.

Figure 1
Types of topological defects in schlieren textures on a thin planar cell of NLCs. (a) A schlieren texture found under POM with crossed polarisers, showing four types of defects characterized by m =  ±1/2 and ±1. ...

The order parameter space of a 3D NLC is S 2/Z 2, and its first homotopy group is Z 2, different from that of a 2D NLC3,5. A direct and notable consequence of this difference is that a line defect is topologically stable only when its winding number m is half-integer; one with integer m can be smoothed out to a configuration with no singularity15,2430. This is achieved by out-of-plane distributions of n, often referred to as “escape to the third dimension” (Fig. 1c–e). This escape occurs naturally to relax the energetic cost of defect core regions, and a well-known example can be found in a cylindrical capillary imposing surface normal alignment1,2430. Analogous structures known as umbilics can be induced by applying an electric field to a NLC with negative dielectric anisotropy1,31.

Although some of the 3D distributions of n in the regions surrounding specific defects in thick cells have been clarified with confocal fluorescence polarised microscopy3234, fine defect structures in schlieren textures of a thin flat cell, especially those involving 3D profiles, remain unexplored experimentally because of the limitation in the resolution of optical microscopy. Moreover, POM does not provide direct information on the component of n normal to the cell, and therefore POM alone cannot capture the escaped configuration of the defects. Here, to uncover the structural states of defects, particularly those with m =  ±1, we present POM and fluorescence optical microscopy (FOM) observations of the defects. Complementary FOM using dyes with different fluorescence anisotropies can confirm the presence of escaped structures at the defects with m =  ±1. Combining experimental POM and FOM images with numerically calculated distributions of n and simple ray-tracing calculations to simulate the lensing effect due to the escaped structures, we uncover four different states of the defects with m =  +1, attributable to the up/down symmetry breaking of the escape, and also the chiral symmetry breaking of the distribution of n (Fig. 1e). Anisotropic configurations of the defects with m = −1 (refs35,36) are also identified (Fig. 1d). We show that the escaped structures cause different lensing effects depending on the direction of the escape and thus, different appearances of defects in the FOM images.


Anisotropic fluorescence emissions of dyes in aligned NLCs

We first investigate the basic properties of the anisotropic fluorescence emissions from the dyes, pyrromethene 597 (PMN) and c-Naphox (CNX)37, in a uniaxially aligned NLC. Figure 2 shows the angle dependent fluorescence emissions with or without an electric field normal to the planar NLC cell. Without the electric field, both dyes show anisotropic emissions. The maximum emissions from PMN and CNX are found at the directions parallel and perpendicular to n = (n x, n y, n z), respectively. The results suggest that the transition dipole moment of the dye, t, projected onto the cell plane lies along n for PMN and perpendicularly to n for CNX on average. With the electric field that aligns n normal to the cell plane, both dyes show isotropic emissions around the axis of n. The results support that t indeed lies along n for PMN (positive fluorescence dichroic property) and perpendicularly to n for CNX (negative fluorescence dichroic property) on average (Fig. 2b,c).

Figure 2
Positive and negative fluorescent dichroic properties of dyes in NLC. (a) Analyser-angle-dependence of relative fluorescence intensities of dye-doped NLCs in a planar cell with the thickness of 2 μm collected using the fluorescence microscopy ...

The total fluorescence intensity without the analyser given by the integral of the plot of Fig. 2a over the angle is much lower under an electric field than without the field in the case of PMN. This suggests that in the FOM images of schlieren textures the region with the larger |n z| will appear darker. In the case of CNX, the relationship between the FOM intensity and |n z| is opposite to that of PMN; the region with larger |n z| will appear brighter. Therefore, using both dyes with positive and negative fluorescence dichroic properties, the distribution of |n z| in the schlieren textures can be qualitatively estimated in a complementary way. (Since the incident angle of the light through the objective lens is not perfectly parallel to the optical axis z in the FOM setup, the degree of polarization is not strictly maintained. Therefore, the present fluorescence dichroic properties should be regarded as qualitative.)

Schlieren textures at different thicknesses

Figure 3a–d shows typical schlieren textures observed at different local cell thicknesses d in a wedge cell (see Methods). The local thickness can be estimated by the colours arising from the thickness dependent retardation of the NLC (5CB) with planar alignment (see Methods). Four regions having local thickness of approximately 1, 2, 3, and 4 μm were analysed in the following analysis.

Figure 3
Four regions with different cell thickness showing typical schlieren textures. Schlieren textures at the approximate thickness d of (a) 1, (b) 2, (c) 3, and (d) 4 μm. The colours arise from the interference due to the thickness dependent retardation ...

The number densities of defects with half-integer and integer m are plotted in Fig. 3e, showing decrease with increasing d. In general, a pair of defects with the same |m| and opposite signs tends to be annihilated after they attract each other3841, which gives rise to coarsening of local nematic domains. In the present system, the transient and local alignment direction at the polyimide surfaces is memorized and frozen in short time scales (several seconds) after the liquid crystal is quenched to the nematic phase from the isotropic phase (surface memory effect4246). Therefore, the observed frozen defects should be a consequence of the competition between the defect annihilation and the surface memorizing processes.

Since the surface memory grows only at the interface, the cell thickness should not affect the characteristic time scale of the growth. On the other hand, the coarsening process of local nematic domains that reduces the number of defects has more chances to proceed in a thicker cell. This may be a qualitative reason for the decrease in the total number of defects with increasing d. In contrast to the sharp decline of number density of defects with half-integer m with d, those with integer m decay more slowly. The defects with integer m, which emerge with the lower density, may have less chance to meet and be annihilated. Although the detail of the survival dynamics of defects is an interesting problem, it is beyond the scope of this paper. For the evaluation of the stable and frozen defects specifically with integer m, the present range of d is technically favourable because a sufficient number of defects for analysis can be acquired from single images captured at one time.

Observation of defects in schlieren textures

Figures 4 and and55 show all types of defects recognised in the schlieren textures of NLCs with PMN and CNX, respectively, at different thicknesses, d. Examples of the original images are shown in Supplementary Figs 1 and 2. The defects with m = ±1/2 appear as bright spots in the FOM images irrespective to the sign of the fluorescent dichroicity, suggesting the accumulation of the dye molecules at the defect cores20,21. Here we focus on the defects with m = ±1, and in FOM images they appear as dark and bright spots when the host NLC is doped with PMN and CNX, respectively. The spot size increases with increasing d. The results indicate that n z is nonzero at the defect centre, and thus, escaped structures exist and the extent of the escaping region increases with increasing d. Defects with m = −1, appear as spots elongated in one axis in FOM images. The estimated 2D distribution of n around the defect reconstructed by images of the POM and FOM with a polariser indicates that the axis of elongation is parallel to one of the axes with radial alignment of (n x, n y).

Figure 4
Discernible states of the defects found in schlieren textures of NLC doped with PMN. A set of microscope images (10 × 10 μm2 squares) of discernible defects with m =  ±1/2 and ±1. ...
Figure 5
Discernible states of the defects found in schlieren textures of NLC doped with CNX. A set of microscope images (10 × 10 μm2 squares) of discernible defects with m =  ±1/2 and ±1. ...

Defects with m =  +1 exhibit four discernible types of patterns in the optical images of the NLC doped with PMN (Fig. 4a–d). The first difference is that dark spots with two different sizes appear at each thickness. Smaller and larger dark spots are denoted by “S” and “L”, respectively, in the brackets in Fig. 4a–d. Examples of the fluorescence profiles at d ~2 μm shown in Fig. 4e indicate the difference. The size of both “S” and “L” increase with increasing d and are comparable to d (Fig. 4f). On the other hand, in the optical images of the NLC doped with CNX, the size difference is not clearly recognised, although the tendency of the increase of the bright spot size with d is confirmed (Fig. 5). The second difference is that the streamline of (n x, n y) from the centre in the plane winds in the clockwise (CW), or counter-clockwise (CCW) manner, which is also shown in the brackets in Fig. 4a–d. It is also noted that the pure radial alignment of (n x, n y) shown in Fig. 1e is practically never found. Such curved orientation profiles as curved dark brushes in POM have been reported in previous literature1,2830.

The presence of two different winding direction is indicative of chiral symmetry breaking. To confirm it, we observed the schlieren textures of a NLC doped with PMN and a small amount of CB15 that primarily induces right-handed twist distortions. Only two of the four discernible types appear for the defects with m =  +1 (Supplementary Fig. 3); smaller dark spots with CCW vortex (S,CCW), and larger dark spots with CW vortex (L,CW). This result indicates that the pair (S,CCW) and (L,CW) is associated with one handedness of chirality, and the other pair (S,CW) and (L,CCW) with the other handedness.

We also obtained confocal laser scanning (CLS) FOM images for defects with m =  ±1 (Fig. 6). Since the excitation laser is linearly polarised in our setup, obtained images strongly reflect the 2D distribution of n, and thus, resemble FOM ones with a polariser. Because of the circumferential contrast around the defect centre on these images, it was difficult to distinguish the difference in the size of the dark spots at the defect centre found for the defect with m =  +1 under FOM without the polariser. Nevertheless, the defect centres with PMN and CNX appear dark and bright, respectively, as in FOM images with a polariser, supporting the existence of the escaped structures.

Figure 6
Confocal laser scanning (CLS) FOM images. Images (10 × 10 μm2 squares) of systems with (a) PMN and (b) CNX at d ~2 μm. From left to right in each panel, transmittance image with crossed-nicols-like configuration, ...


Structure of defects

Let us discuss the origin of the four different states found for defects with m =  +1. We first recall that the escape direction can be either upward or downward as shown in Fig. 1c, and that the configurations with different escape directions are of the same free energy because the two confining substrates are equivalent.

The other important contribution is that of chiral symmetry breaking as mentioned above. The chiral symmetry breaking is attributed to twist distortions due to the twist elastic constant (K 22) smaller than those for splay and bend distortions (K 11 and K 33, respectively), and give rise to configurations with different twist sense. The up/down symmetry breaking and the chiral symmetry breaking result in four distinct configurations, shown as [U,L,CW], [U,R,CCW], [D,R,CW] and [D,L,CCW] in Fig. 1e. Here we denote the escape direction by U and D (Up/Down), the sense of the twist of n by L and R (Left-handed/Right-handed), and the sense of the spiral profile of n projected onto the cell plane by CW and CCW (ClockWise/Counter-ClockWise). It is important to note that [U,L,CW] and [D,L,CCW] are the same profile viewed from the different directions ([U,R,CCW] and [D,R,CW] are also the same). A chiral dopant inducing right-handed twist was used in the experiment shown in Supplementary Fig. 3, and therefore the profiles (L,CCW) and (S,CW) found there should be assigned to [U,R,CCW] or [D,R,CW], respectively. From the same argument for the structures with left-handed twist, [D,L,CCW] or [U,L,CW] should correspond to (S,CCW) and (L,CW). Thus the director profiles of the experimentally observed four distinct states are successfully assigned, and the size of the dark spot of the defect (S/L) is associated to the direction of the escape (U/D).

To show that this chiral symmetry breaking indeed occurs for defects with m =  +1, we calculate the optimised profiles of orientational order within the framework of the Landau-de Gennes formulation. The elastic energy in our Landau-de Gennes formulation corresponds to choosing K 22 < K 11 = K 33 in the Frank elastic energy (see Methods for details). Although we do not focus on the effect of unequal K 11 and K 33, we note that radial far-field orientation profiles with splay distortions are realised when K 11 < K 33 as is the case for 5CB, and circular profiles with bend distortions when K 11 > K 33.

We relax the orientation profile from the initial profile, an ideal escaped structure without chiral symmetry breaking shown in Fig. 1e framed by a black square. Chiral symmetry breaking occurs yielding the profile shown in Fig. 7a despite the absence of randomness in the initial profile and the relaxation process, which indicates the robustness of the chiral symmetry breaking (rounding numerical errors are sufficient for the symmetry breaking) and the instability of the profiles without chiral symmetry breaking ([U,N,R] and [D,N,R] in Fig. 1e). The profile in Fig. 7a exhibits right-handed twist distortions, and that with left-handed twist distortions can emerge with equal probability because of the inversion symmetry of the free energy and the boundary conditions.

Figure 7
3D distribution of n around the defects with integer m. Calculated structures of defects with (a) m =  +1 and (b) m = −1. Here the director n is the unit eigenvector of the tensor order parameter Q ij with ...

Now let us focus on the structure of a defect with m = −1. The projected 2D profile of a defect with m = −1 (Fig. 1b) appears to have C4v symmetry. However, the actual 3D profile does not when the escape is appropriately taken into account (Fig. 1d)27,35,36; the escape is upward along one lateral axis (“Up-axis” in Fig. 1d) and downward along an axis perpendicular to the Up-axis (“Down-axis”). The symmetry of the m = −1 defect is therefore C2v, and the elongated spots observed in experiments are consistent with the C2v symmetry rather than the C4v symmetry. The direction of the elongation is parallel to one of the escape axes (Up-axis or Down-axis). The escaped structure shown in Fig. 1d intrinsically involves twist distortions; e.g., those with different handedness can be identified along the diagonal directions of the green square in Fig. 1d. Therefore, the defect with m = 1 cannot break the chiral symmetry in the manner that the defect with m =  +1 does because lower K 22 would simply boost both twist distortions with different handedness locally to minimise the total energy without breaking the original C2v symmetry. In Fig. 7b, we show the numerically calculated structure of a defect with m = −1, and the initial condition is an ideal one shown schematically in Fig. 1d. We also perform different calculations with chiral perturbations introduced in the ideal initial condition. These perturbations turned out to decay resulting in an achiral profile, which indicates the robustness of the configuration shown in Fig. 1d.

From the numerical calculations, we also find that the width of the escape remains finite, although in the case of a cylindrical capillary1,2426, it is equal to the capillary radius (that is, the larger the radius, the larger the escape width). In our case of a thin cell, cell thickness d instead of the capillary radius enters as the characteristic length, and therefore the escape width should be proportional to d. This speculation is consistent with the observation in Fig. 4f that the dark spot size corresponding to the escape width increases almost linearly with d.

Optics of defects

Our remaining task is to clarify why the direction of the escape (up/down) can manifest itself as the size of the dark spot. For this purpose, we investigate how light rays are bent through the LC cell near the defects with a non-uniform distribution of refractive indices47. The bend of light through a LC has long been known4850 and applied to the LC lens technology51. If the lensing effect were absent, the up/down-states would be indistinguishable because the intensity of the FOM would depend only on (n x,n y) and therefore insensitive to the escape direction. Our observations motivate us to evaluate the lensing effect here.

To simplify the problem for our qualitative arguments, we consider a defect with m =  +1 without twist deformation shown in Fig. 1e (see Methods for the calculation detail). We restrict our discussion to rays confined in the 2D plane (r,z) that includes the z axis (r is the distance from the z axis). Since the director n is also confined in the (r,z) plane in our setup, the extraordinary light (where the electric field E is in the (r,z) plane) and the ordinary light (where E is perpendicular to the (r,z) plane) are decoupled, and the latter propagates without bending. Therefore, we discuss only the extraordinary component. Our calculations presented below are based on geometric optics and the principle of Fermat5,52.

The transition dipole moment t of PMN on average aligns along n (Fig. 2b), and thus, the extraordinary component dominates both excitation and emitted fluorescence rays (Fig. 8a). On the other hand, t is distributed perpendicularly to n for CNX (Fig. 2c). Therefore, excitation and emitted fluorescence rays contain both the ordinary and the extraordinary components depending on |n z|; the extraordinary component grows as |n z| varies from 0 (a peripheral part of the defect) to 1 (the defect centre) (Fig. 8b). This simply suggests that the systems with PMN and CNX show different lensing effects; a clearer effect is expected for the system with PMN than that with CNX, in agreement with the experimental result that the up/down-states are indistinguishable from the microscopy images for the latter. The following discussion on extraordinary rays therefore concerns microscope images on the system with PMN.

Figure 8
Ray-tracing around the ideal escaped structure. Schematics of the relationship between the nematic director, n, and the transition dipole moment of the dye, t, on a radial 2D plane with the escaped structure in the system with (a) PMN and (b) CNX. The ...

The calculated rays entering from different points near the defect centre at the bottom side of the cell with different incident angles are shown for both up/down-states in Fig. 8c,d. Although both states can bend the rays, the incident rays clearly tend to deviate more from the defect centre for the down-state. In the following, we discuss qualitatively how the bending of rays modulates 2D distributions of the excitation intensity and distort the resulting fluorescence images.

The excitation rays entering from the top side of the cell with the down- (or up-) state roughly correspond to those entering from the bottom side of the cell with the up- (or down-) state shown in Figs 8c (or 8d). Thus, for the down-state the excitation intensity is only slightly blurred and almost homogeneous. For the up-state, the excitation rays deviate toward the radial direction as rays travel deeper into the cell, resulting in the decreased intensity of excitation at the centre of the defect and the increased intensity roughly at the peripheral parts of the escaped region (r/r 0 ~1, where r 0 is the radius of the escaped region), in which the number density of the bent rays is highest at z/r 0 = 1~2 in Fig. 8d. Since the emitted fluorescence rays from the up-state propagate with the little bending following the traces shown in Fig. 8c, the resulting fluorescence image of the defect is expected to be a simple manifestation of the actual profile of |n z| in the radial direction with slightly increased intensity at the peripheral position.

On the other hand, the emitted fluorescence rays from the down-state, which has been excited almost homogeneously, bend toward the radial direction as shown in Fig. 8d and the deviation becomes larger for rays traveling from deeper parts. The ray from the point X at the vicinity of the defect centre in Fig. 8d bends in the radial direction, and thus, seems to originate at the point Y which is deviated from X towards the defect centre. As a result, the fluorescence intensity distribution is expected to be modulated and distorted so that the escaped region in the resulting fluorescence image appears as a smaller dark region. These qualitative arguments agree well with the experimental results, in which the dark region of the down-state appear smaller than that of the up-state and the fluorescence intensity in the peripheral region surrounding the dark spot is larger than that at a large distance from the defect (Fig. 4).

The same argument applies to the identification of the up- and down-axes of a defect with m = −1 (Fig. 1d) from the FOM images in Fig. 4; the direction of the elongation of the dark spot indicated by yellow arrows in Fig. 4 corresponds to the up-axis. Also in the system with CNX, although they are not so clear as the system with PMN, some of the bright spots in the FOM images of defects with m = 1 at larger d show elongated shapes (Fig. 5), in which the elongated direction should correspond to the up-axis.

Overall, the systematic modulations of the defect images due to the lensing effect are critical for the distinction of the up/down-states. For more precise discussions, rigorous optical treatment based on the solution of the full Maxwell equations will be required instead of the geometrical optics and the ray-tracing method53, which remains to be studied in future.

In summary, using a conventional optical microscope, we have uncovered the hidden states of the defects with integer strength, m, in the classical schlieren textures. For defects with m = +1, four states are identified as the result of up-down and chiral symmetry breaking. For those with m = −1, the directions of the two characteristic axes are distinguished. Two types of fluorescent dyes with different alignments of the transition dipole moment with respect to the local nematic director n enable us to confirm the escaped structures of defects with integer m. Through the Landau-de Gennes numerical modelling, we have calculated their detailed escaped structures and associated them with the observed different states. The simulation results also suggest that the chiral symmetry-broken structures showing 2D spiral-like streamline of n emerge only for defects with m =  +1 due to the lower elastic constant for the twist deformation than those of splay and bend ones, in agreement with the experimental results. Moreover, we have attributed the fluorescence images modulated depending on the escaped direction of n to the LC lensing effect due to the characteristic 3D distribution of the refractive indices around the defect. This suggests the importance of considering the lensing effect on the image analysis of more complex LC defects, e.g., defects with higher values of m (refs15,54), to reveal their hidden structural details. Clarification of the detailed structures of various defects in LCs could reveal their potential as unique templates for micro-patterning of asymmetric and/or chiral objects ranging from small molecules to colloidal particles. Since the escaped structures of the defects with m =  ±1 can locally and systematically modulate light depending on the size and the states, they also have potential as micro-light-modulators for gradient-index micro-lenses and a fine maskless direct patterning technique in photo-lithography.



A NLC, 4-cyano-4′-pentylbiphenyl (5CB), whose twist elastic constant is significantly lower than those for splay and bend distortions55, was purchased from Wako Chemical. A fluorescent dye pyrromethene 597 (PMN, [[(4-tert-Butyl-3,5-dimethyl-1H-pyrrol-2-yl)(4-tert-butyl-3,5-dimethyl-2H-pyrrol-2-ylidene)methyl]methane](difluoroborane) was purchased from Sigma-Aldrich. Another fluorescent dye C-Naphox (CNX, racemic mixture) was synthesized according to the method previously reported by some of the authors37. Poly-amic acid (Poly(pyromellitic dianhydride-co-4,4′-oxydianiline), and N-methylpyroridone (NMP) for a precursor solution used for polyimide coating with planar alignment were obtained from Sigma-Aldrich. As a right-handed chiral reagent, 4-cyano-4′-((S)-2-methylbutyl)biphenyl (CB15, Merck) was used to investigate the effect of chirality.

Preparation of planar degenerate LC cells

We adopted a wedge cell with planar alignment without preferential direction (degenerate planar alignment) to observe very thin regions at the thickness of roughly 1 to 4 μm. A spin-coater (MSA-150, Mikasa) was used at 5,000 rpm for 60 s to coat cover glass slides with a poly-amic acid solution (1 wt%). The coated slides were then heated on a hot plate at 70 °C for 2 min to evaporate the solvent, and subsequently placed in an oven at 180 °C for 3 h to induce the transformation of poly(amic acid) to polyimide through the dehydration reaction. Then, two slides were assembled to a wedge cell using a Kapton® spacer placed at one side of the cell with a small amount of adhesive (liquid gasket1212, ThreeBond). The slides were held together using clips during curing of the adhesive. Although the wedge angle was ideally ~0.036 degree, the actual cell showed spatial fluctuation of the wedge angle (and the cell gap), probably due to the bending of the cover glass and inhomogeneity in the thickness of the polyimide layer. The NLC, 5CB, doped with 0.01 wt% PMN or 0.05 wt% CNX was injected into the wedge cell at an elevated temperature (40~45 °C) via capillary action. The sample was cooled to room temperature (22 ± 2 °C) and used for microscopic observations. Samples doped with CB15 at ~0.2 wt%, which gives the cholesteric pitch much larger than d, were prepared employing the same procedure.

For measurements of fluorescence dichroic properties, a planar cell with unidirectional planar alignment (KSRP-02/A111P1NSS05, EHC; cell gap of 2 ± 0.5 μm) filled with 5CB doped with dyes was used.

Optical microscopy

We observed the LC alignment using a transmitted polarizing optical microscope (POM) (BX-51P, Olympus) under crossed-nicols conditions.

To observe fluorescent images, we used a conventional fluorescence optical microscope (FOM), where the light source was mounted above the sample and the excitation light passed through the microscope objective lens on its way toward the sample. A Xe lamp (75 W) was used as the light source. To detect the emission from PMN and CNX, we used a fluorescence filter set (U-MWIB-3, Olympus) comprising an excitation filter that transmitted light with wavelengths between 460 nm and 495 nm, and an emission filter that transmitted light with wavelength larger than 510 nm. The images were collected at the pixel size of 0.154 μm using a Nikon DS-Qi1-Mc CCD (charge-coupled device) camera connected to a computer and controlled through imaging software (NIS-Elements, Nikon). An objective lens with a numerical aperture (NA) of 0.95 (UPLSAPO40 × 2, Olympus) was mainly used.

For the measurements of dichroic properties, the fluorescence intensity was measured using the FOM with a polarizer (analyser), through which the emitted fluorescence passed. The average intensities were measured at different angles, Φ, between the axis of the polarizer and the nematic director of the cell. When the maximum and minimum intensities were at Φ = 0 degree and ±90 degree, respectively, the transition moment of dye, t, was aligned to the nematic director, n, and the dye had positive dichroic fluorescence property. When the intensity variation behaves in an opposite manner, t was perpendicularly aligned to n, and the dye had negative dichroic fluorescence property. To further confirm the alignment of dyes in 5CB, the same measurements were conducted under the AC voltage to align n along the optical axis, that is, perpendicular to the cell.

We obtained microscope images in four regions exhibiting bluish white, yellow, blue, and light green colours in the POM under the crossed-nicols condition. Their thicknesses were estimated to be approximately 1, 2, 3, and 4 μm (error ± 0.5 μm) using the Michel-Levy interference chart [see, e.g., Olympus Microscopy Resource Centre homepage] with the birefringence of 5CB being ~0.18 at 23 °C56. For image acquisition, the focus was carefully adjusted to give tack-sharp peaks of defects with m =  ±1/2 (ref.20). The contrast of fluorescence images was linearly enhanced to maximize the visibility.

A confocal laser scanning (CLS) FOM (A1+ system, Nikon) was used to obtain the fluorescence images of defects with integer m. An optically pumped semiconductor laser (LU-N4 Laser Unit, Nikon, equipped with Sapphire 488, Coherent Inc.) was used to excite the fluorescent molecules at 488 nm and the emitted light between 525 and 595 nm was collected. An objective lens with an NA of 1.45 (PlanApoTIRF60 ×, oil, Nikon) was used. The CLS-FOM images were acquired at the pixel size of typically 60(x) × 60(y) nm2. The excitation laser was linearly polarized in y direction and all emitted fluorescence was collected.

Landau-de Gennes numerical modelling

A second rank tensor (Q ij) is used for the orientational order parameter, and the free energy of the system is formally expressed as


where d is the thickness of the cell and z is the coordinate along the cell normal. The bulk parts of the free energy density read15



where summations over repeated indices are implied, (∇×Q)ij ≡ ϵistsQtj, and ([nabla] [center dot] Q)j [equivalent] [nabla] i Q ij, where ϵist is the Levi-Civita antisymmetric symbol. We used typical material parameters for 5CB (ref.57); A = 0.044 × 106(T  T *) J m−3 K−1 with T * = 307 K and T = 296 K (room temperature), B = 0.816 × 106 J m−3, C = 0.45 × 106 J m−3, L 1 = 6 × 10−12 J m−1 and L 2 = 18 × 10−12 J m−1. Note that when a uniaxial form of the order parameter, Q ij = Q 0(n i n j  (1/3)δ ij), is assumed, f grad is reduced to the Frank elastic energy with K 11=K 33=(L 1 + L 2)Q 0 2 and K 22 = 2L 1 Q 0 2. Therefore, our choice of L 1 and L 2 yields K 22 smaller than K 11(=K 33) as is the case for 5CB.

The surface free energy density for two confining surfaces, f s0 and f sd were taken to be identical in our calculations, and we adopted the following form commonly used to model surfaces with planar degenerate anchoring58:


where Q s defines the strength of orientational order at the surface Q˜ijQij+(1/3)Qsδij, and Q˜ijPikQ˜klPlj with P ij = δ ij  ν i ν j being the projection operator where ν is a unit vector normal to the surface. We used W 1 = W 2 = 10−4 J m−2, a typical value for moderate anchoring strength.

The above-mentioned free-energy functional was minimised on a 201 × 201 × 51 cubic lattice with the lattice constant being Δx ≃ 20 nm (Therefore the cell thickness is approximately 1 μm). The minimisation was carried out by solving a simple rotational relaxation equation [partial differential]Q ij (r)/[partial differential]t = δF/δQ ij (r) + λδ ij, where the Lagrange multiplier λ ensures TrQ = 0.

The initial profile with an “escape” for m = +1 was set to be uniform along the z direction, and taken to be Qij(xyz) = Qinit(ni(xy)nj(xy)−(1/3)δij), where the director n(x,y) was set to the following:


where Q init = 1, r=(xxc)2+(yyc)2 with (x c, y c) being the coordinates of the centre of the cell, r 0 = 10Δx is the radius of the escape, e r = (x  x c,y  y c,0)/r is a radial unit vector, with and e z is a unit vector along the z direction. θ(r) = 2tan−1(r/r c) is the profile of escape obtained analytically for the “one-constant” case1,24,25, and this choice of θ(r) yields the down-state. We fixed n to e r at the outer boundaries in the (x,y) plane. Although this initial profile is achiral, it turned out that rounding errors were enough to induce the transition to a chiral state mentioned in the main text.

The initial profile for m = −1 is the same as that for m = +1, with n being flipped about the x axis locally. We performed another calculation with chiral perturbations introduced in the above-mentioned initial profile for m = −1, and these perturbations turned out to decay, resulting in an achiral profile mentioned in the main text.

Ray-tracing simulation

As noted in the main text, the main focus of our ray-tracing calculations is to provide a qualitative argument on the appearance of the escaped defects. Therefore, we consider a simple case of achiral configuration of the orientational order given in eq. (5) for rr 0, in which the director n lies in the (e r, e z) plane, and so does the light ray of interest. Now θ (r) = ± 2tan−1 (r/r 0), where + and  are for the down state and the up state, respectively.

According to Fermat’s principle, the path of extraordinary light from point A to B minimises the following integral5:


where dl is the line element of the path, n o and n e are the ordinary and the extraordinary refractive indices, respectively, and ϕ is the angle between the tangent of the path and n. In our calculation, we set n e = 1.7 and n o = 1.5, typical values for 5CB. When we denote the path by r(z), dl=dz1+(dr/dz)2. The unit vector tangent to the path is t=(dr/dz,1)/1+(dr/dz)2, and therefore cosφ=nt=(sinθ(r(z))dr/dz+cosθ(r(z)))/1+(dr/dz)2. The integral to be minimised is finally written as F=ABdzf{r(z),dr/dz} with


The Euler-Lagrangian equation [partial differential]f/[partial differential]r  (d/dz)([partial differential]f/[partial differential](dr/dz)) = 0 was solved numerically with the aid of Maple 2016.2 for given r(0). An additional initial condition must be given for dr/dz|z = 0, which is determined by appropriately taken into account the refraction at the interface (z = 0) between the liquid crystal and the outer medium (whose refractive index was set to 1.5).

Electronic supplementary material

SI-file(706K, pdf)


The authors would like to thank H. Yoshida and Y. Sato for useful discussions and T. Yamamoto for supplying the chiral reagent. T.O. and A.F. were supported by a grant for promotion of cooperation between Nagoya University and AIST. T.O. and A.F. were also supported by Grants for Fusion Emergent Research from Integrated Research Consortium on Chemical Sciences, Japan. J.F. was supported by JSPS KAKENHI Grant Number JP17H02947. K.K. was supported by JSPS KAKENHI Grant Numbers JP17H06417 and JP17H06413.

Author Contributions

Author Contributions

T.O. conceived the central idea, performed the experiments and analysed the results. J.F. carried out the numerical calculations. T.O. and J.F. contributed to the theoretical arguments. K.K. and T.O. arranged CLS-FOM and performed observations using CLS-FOM. C.W., A.F., and S.Y. designed and synthesised the dye, c-Naphox (CNX). All authors discussed the results and T.O. and J.F. wrote the manuscript.


Competing Interests

The authors declare that they have no competing interests.


Electronic supplementary material

Supplementary information accompanies this paper at 10.1038/s41598-017-16967-1.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributor Information

Takuya Ohzono, pj.og.tsia@ayukat-onozho.

Jun-ichi Fukuda,


1. de Gennes, P. G. & Prost, J. Physics of Liquid Crystals (Clarendon, Oxford, 1993).
2. Oswald, P. & Pieranski, P. Nematic and cholesteric liquid crystals (Taylor & Francis, New York, 2005).
3. Mermin ND. The topological theory of defects in ordered media. Rev. Mod. Phys. 1979;51:591–648. doi: 10.1103/RevModPhys.51.591. [Cross Ref]
4. Lavrentovich OD. Topological defects in dispersed words and worlds around liquid crystals, or liquid crystal drops. Liq. Cryst. 1998;16:117–124. doi: 10.1080/026782998207640. [Cross Ref]
5. Kleman, M. & Lavrentovich, O. D. Soft Matter Physics: An Introduction. (Springer–Verlag, New York, 2003).
6. Chuang I, Durrer R, Turok N, Yurke B. Cosmology in the laboratory: defect dynamics in liquid crystals. Science. 1991;251:1336–1342. doi: 10.1126/science.251.4999.1336. [PubMed] [Cross Ref]
7. Bowick MJ, Chandar L, Schiff EA, Srivastava AM. The Cosmological Kibble Mechanism in the Laboratory: String Formation in Liquid Crystals. Science. 1994;263:943–945. doi: 10.1126/science.263.5149.943. [PubMed] [Cross Ref]
8. Poulin P, Stark P, Lubensky HTC, Weitz DA. Novel Colloidal Interactions in Anisotropic Fluids. Science. 1997;275:1770–1773. doi: 10.1126/science.275.5307.1770. [PubMed] [Cross Ref]
9. Pires D, Fleury J-B, Galerne Y. Colloid particles in the interaction field of a disclination line in a nematic phase. Phys. Rev. Lett. 2007;98:247801. doi: 10.1103/PhysRevLett.98.247801. [PubMed] [Cross Ref]
10. Fleury J-B, Pires D, Galerne Y. Self-connected 3D architecture of microwires. Phys. Rev. Lett. 2009;103:267801. doi: 10.1103/PhysRevLett.103.267801. [PubMed] [Cross Ref]
11. Tkalec U, Ravnik M, Čopar S, Žumer S, Muševič I. Reconfigurable knots and links in chiral nematic colloids. Science. 2011;333:62–65. doi: 10.1126/science.1205705. [PubMed] [Cross Ref]
12. Ravnik M, Alexander GP, Yeomans JM, Žumer S. Three-dimensional colloidal crystals in liquid crystalline blue phases. Proc. Natl. Acad. Sci. USA. 2011;108:5188–5192. doi: 10.1073/pnas.1015831108. [PubMed] [Cross Ref]
13. Ohzono T, Fukuda J. Zigzag line defects and manipulation of colloids in a nematic liquid crystal in microwrinkle grooves. Nat. Commun. 2012;3:701. doi: 10.1038/ncomms1709. [PMC free article] [PubMed] [Cross Ref]
14. Sengupta A, Bahr C, Herminghaus S. Topological microfluidics for flexible micro-cargo concepts. Soft Matter. 2013;9:7251–7262. doi: 10.1039/c3sm50677k. [Cross Ref]
15. Yoshida H, Asakura K, Fukuda J, Ozaki M. Three-dimensional positioning and control of colloidal objects utilizing engineered liquid crystalline defect networks. Nat. Commun. 2015;6:7180. doi: 10.1038/ncomms8180. [PMC free article] [PubMed] [Cross Ref]
16. Kikuchi H, Yokota M, Hisakado Y, Yang H, Kajiyama T. Polymer-stabilized liquid crystal blue phases. Nat. Mater. 2002;1:64–68. doi: 10.1038/nmat712. [PubMed] [Cross Ref]
17. Higashiguchi K, Yasui K, Ozawa M, Odoi K, Kikuchi H. Spatial distribution control of polymer nanoparticles by liquid crystal disclinations. Polymer J. 2012;44:632–638. doi: 10.1038/pj.2012.44. [Cross Ref]
18. Wang X, Miller DS, Bukusoglu E, de Pablo JJ, Abbott NL. Topological defects in liquid crystals as templates for molecular self-assembly. Nat. Mater. 2016;15:106–112. doi: 10.1038/nmat4421. [PubMed] [Cross Ref]
19. Wang X, et al. Experimental insights into the nanostructure of the cores of topological defects in Liquid Crystals. Phys. Rev. Lett. 2016;116:147801. doi: 10.1103/PhysRevLett.116.147801. [PubMed] [Cross Ref]
20. Ohzono T, Katoh K, Fukuda J. Fluorescence microscopy reveals molecular localisation at line defects in nematic liquid crystals. Sci. Rep. 2016;6:36477. doi: 10.1038/srep36477. [PMC free article] [PubMed] [Cross Ref]
21. Fukuda J. Do small molecules assemble at topological defects of a liquid crystal? Proc. SPIE. 2017;10125(101250C):1–7.
22. Nehring J, Saupe A. On the schlieren texture in nematic and smectic liquid crystals. J. Chem. Soc. Faraday Trans. 1972;2(68):1–15. doi: 10.1039/f29726800001. [Cross Ref]
23. Chiccoli C, et al. Topological defects in schlieren textures of biaxial and uniaxial nematics. Phys. Rev. E. 2002;66:030701(R). doi: 10.1103/PhysRevE.66.030701. [PubMed] [Cross Ref]
24. Meyer RB. On the existence of even indexed disclinations in nematic liquid crystals. Phil. Mag. 1972;27:405–424. doi: 10.1080/14786437308227417. [Cross Ref]
25. Saupe A. Disclinations and properties of the directorfield in nematics and cholesteric liquid crystals. Mol. Cryst. Liq. Cryst. 1973;21:211–238. doi: 10.1080/15421407308083320. [Cross Ref]
26. Cladis PE, Kleman M. Non-singular disclinations of strength s =  +1 in nematics. J. Phys. (France) 1972;33:591–598. doi: 10.1051/jphys:01972003305-6059100. [Cross Ref]
27. Williams CE, Cladis PE, Kleman M. Screw disclinations in nematic samples with cylindrical symmetry. Mol. Cryst. Liq. Cryst. 1973;21:355–373. doi: 10.1080/15421407308083329. [Cross Ref]
28. Candau S, LeRoy P, Debeauvais F. Magnetic field effects in nematic and cholesteric droplets suspended in an isotropic liquid. Mol. Cryst. Liq. Cryst. 1973;23:283–297. doi: 10.1080/15421407308083378. [Cross Ref]
29. Press MJ, Arrott AS. Theory and experiments on configurations with cylindrical symmetry in liquid-crystal droplets. Phys. Rev. Lett. 1974;33:403–406. doi: 10.1103/PhysRevLett.33.403. [Cross Ref]
30. Press, M. J. & Arrott, A. S. Elastic energies and director fields in liquid crystal droplets, I cylindrical symmetry. J. Phys. (France) 36, C1–117–184 (1975).
31. Sasaki Y, et al. Large-scale self-organization of reconfigurable topological defect networks in nematic liquid crystals. Nat. Commun. 2016;7:13238. doi: 10.1038/ncomms13238. [PMC free article] [PubMed] [Cross Ref]
32. Smalyukh II, Shiyanovskii SV, Lavrentovich OD. Three-dimensional imaging of orientational order by fluorescence confocal polarizing microscopy. Chem. Phys. Lett. 2001;366:88–96. doi: 10.1016/S0009-2614(00)01471-8. [Cross Ref]
33. Trivedi RP, Klevets II, Senyuk B, Lee T, Smalyukh II. Reconfigurable interactions and three-dimensional patterning of colloidal particles and defects in lamellar soft media. Proc. Natl. Acad. Sci. USA. 2012;109:4744–4749. doi: 10.1073/pnas.1119118109. [PubMed] [Cross Ref]
34. Trivedi, R. P. & Smalyukh, I. I. 10. Optical microscopy of soft matter systems. Fluids, Colloids and Soft Materials: An Introduction to Soft Matter Physics, eds Fernandez-Nieves, A. & Puertas, A. M. (Wiley), pp 165–184 (2016).
35. Lubensky TC, Petty D, Currier N, Stark H. Topological defects and interactions in nematic emulsions. Phys. Rev. E. 1998;57:610–625. doi: 10.1103/PhysRevE.57.610. [Cross Ref]
36. Fukuda J, Yokoyama H. Separation-independent attractive force between like particles mediated by nematic-liquid-crystal distortions. Phys. Rev. Lett. 2005;94:148301. doi: 10.1103/PhysRevLett.94.148301. [PubMed] [Cross Ref]
37. Wang C, et al. A Phosphole oxide based fluorescent dye with exceptional resistance to photobleaching: a practical tool for continuous imaging in STED microscopy. Angew. Chem. Int. Ed. 2015;54:15213–15217. doi: 10.1002/anie.201507939. [PubMed] [Cross Ref]
38. Yurke B, Pargellis AN, Kovacs T, Huse DA. Coarsening dynamics of the XY model. Phys. Rev. E. 1993;47:1525–1531. doi: 10.1103/PhysRevE.47.1525. [PubMed] [Cross Ref]
39. Toth G, Denniston C, Yeomans JM. Hydrodynamics of topological defects in nematic liquid crystals. Phys. Rev. Lett. 2002;88:105504. doi: 10.1103/PhysRevLett.88.105504. [PubMed] [Cross Ref]
40. Dierking I, et al. Anisotropy in the annihilation dynamics of umbilic defects in nematic liquid crystals. Phys. Rev. E. 2012;85:021703. doi: 10.1103/PhysRevE.85.021703. [PubMed] [Cross Ref]
41. Digal S, Ray R, Saumia PS, Srivastave AM. Duality between the dynamics of line-like brushes of point defects in 2D and strings in 3D in liquid crystals. J. Phys. Condens. Matter. 2013;25:404204. doi: 10.1088/0953-8984/25/40/404204. [PubMed] [Cross Ref]
42. Cheng J, Boyd GD. The liquid-crystal alignment properties of photolithographic gratings. Appl. Phys. Lett. 1979;35:444. doi: 10.1063/1.91166. [Cross Ref]
43. Yokoyama H, Kobayashi S, Kamei H. Temperature dependence of the anchoring strength at a nematic liquid crystal-evaporated SiO interface. J. Appl. Phys. 1987;61:4501–4518. doi: 10.1063/1.338411. [Cross Ref]
44. Clark NL. Surface memory effects in liquid crystals: influence of surface composition. Phys. Rev. Lett. 1985;55:292–295. doi: 10.1103/PhysRevLett.55.292. [PubMed] [Cross Ref]
45. Berberi R, et al. Azimuthal anchoring of nematic on undulated substrate: elasticity versus memory. Eur. Phys. J. B. 1998;6:83–91. doi: 10.1007/s100510050529. [Cross Ref]
46. Ohzono T, Monobe H, Yamaguchi R, Shimizu Y, Yokoyama H. Dynamics of surface memory effect in liquid crystal alignment on reconfigurable microwrinkles. Appl. Phys. Lett. 2009;95:014101. doi: 10.1063/1.3167547. [Cross Ref]
47. Satiro C, Moraes F. Lensing effects in a nematic liquid crystal with topological defects. Eur. Phys. J. E. 2006;20:173–178. doi: 10.1140/epje/i2005-10127-2. [PubMed] [Cross Ref]
48. Nicholson TM. Computer simulation of liquid crystalline anisotropic structures. Mol. Cryst. Liq. Cryst. 1989;177:163–175.
49. Hwang DK, Rey AD. Computational modeling of light propagation in textured liquid crystals based on the finite-difference time-domain (FDTD) method. Liq. Cryst. 2005;32:487–497.
50. Mendoza CI, Reyes JA. Light propagation and transmission in hybrid-aligned nematic liquid crystal cells: geometrical optics calculations. Appl. Phys. Lett. 2006;89:091912. doi: 10.1063/1.2345042. [Cross Ref]
51. Sato S, Nose T, Yamaguchi R, Yanase S. Relationship between lens properties and director orientation in a liquid crystal lens. Liq. Cryst. 1989;5:1435–1442. doi: 10.1080/02678298908027781. [Cross Ref]
52. Born, M. & Wolf, M. Principles of Optics 7th ed. (Cambridge University Press, Cambridge, 1999).
53. Nych, A., Fukuda, J., Ognysta, U., Žumer, S. & Muševič, I. Spontaneous formation and dynamics of half-skyrmions in a chiral liquid-crystal film., Nat. Phys. 10.1038/nphys4245 (2017).
54. Kralj S, Murray BS, Rosenblatt C. Decomposition of strongly charged topological defects. Phys. Rev. E. 2017;95:042702. doi: 10.1103/PhysRevE.95.042702. [PubMed] [Cross Ref]
55. Madhusudana NV, Pratibha R. Elasticity and orientational order in some cyanobiphenyls: part IV. reanalysis of the data. Mol. Cryst. Liq. Cryst. 1982;89:249–257. doi: 10.1080/00268948208074481. [Cross Ref]
56. Horn RG. Refractive indices and order parameters of two liquid crystals. J. de Phys. 1978;39:105–109. doi: 10.1051/jphys:01978003901010500. [Cross Ref]
57. Coles H. Laser and Electric Field Induced Birefringence Studies on the Cyanobiphenyl Homologues. Mol. Cryst. Liq. Cryst. 1978;51:67–71. doi: 10.1080/00268947808070330. [Cross Ref]
58. Fournier J-B, Galatola P. Modeling planar degenerate wetting and anchoring in nematic liquid crystals. Europhys. Lett. 2005;72:403–409. doi: 10.1209/epl/i2005-10253-5. [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group