Abstract
Topology has an increasingly important role in the physics of condensed matter, quantum systems, material science, photonics and biology, with spectacular realizations of topological concepts in liquid crystals. Here we report on longlived hidden topological states in thermally quenched, chiral nematic droplets, formed from stringlike, triangular and polyhedral constellations of monovalent and polyvalent singular point defects. These topological defects are regularly packed into a spherical liquid volume and stabilized by the elastic energy barrier due to the helical structure and confinement of the liquid crystal in the microsphere. We observe, for the first time, topological threedimensional point defects of the quantized hedgehog charge q=−2, −3. These highercharge defects act as ideal polyvalent artificial atoms, binding the defects into polyhedral constellations representing topological molecules.
Introduction
Topology currently emerges as a major research theme in a number of subfields of physics, and many new topological phenomena have been observed in very different contexts. For example, nontrivial topology of the electronic structure of crystalline materials leads to exotic material properties, such as surface conductivity of topological insulators^{1,2} and peculiar electronic structure of Weyl and Dirac semimetals^{3,4,5,6}. Topologically important phenomena can nowadays be traced in condensed matter^{1}, quantum systems^{7}, material science^{8}, photonics^{9}, active matter^{10,11,12} and even cell division processes in biology^{13}. Understanding the topological properties of matter and applying topological concepts across various subfields of physics not only helps us understand the physics of these systems, but can inspire design of new materials with unusual material properties.
Liquid crystals (LCs) are well known for their richness of topological defects and phenomena, which are observed and analysed on the micrometre scale using optical methods. During the past decade, several spectacular realizations of topological concepts have been demonstrated in LCs, such as topological charge creation and manipulation in nematic liquid crystals (NLCs)^{14,15,16}, and fascinating defect motion in active NLCs^{11,12}. By using laser tweezers and fluorescent confocal polarized microscopy (FCPM) it is possible to analyse the topological charges, knot and link the tensorial nematic ordering field^{17,18,19} and experimentally prove fundamental theorems, such as the PoincaréHopf theorem^{20,21}.
Closely related to topology are hidden states of matter, which can be stabilized by topology and refer to states that are not accessible under equilibrium conditions, but can be created if the system is rapidly quenched^{22}. Optically induced hidden states were reported in metallic glasses^{23} and layered, quasitwodimensional electronic crystals^{24} by applying strong laser pulses. In LCs, transient hidden states of topological defects are obtained by a rapid temperature or pressure quench across the isotropictonematic phase transition. This quench generates a random constellation of topological defects via the KibbleZurek mechanism^{25,26,27}, consisting of mutually compensating topological charges, which are unstable and annihilate into vacuum. Topological defects of unit charge can be stabilized in LCs by colloidal inclusions, appearing as micrometrescale, looplike structural imperfections, accompanying spheres^{17}, fibres^{14}, handlebodies^{20} and knotted particles^{19}. In chiral nematic liquid crystals, topological monopoles of unit charge are stabilized by the spontaneous winding of the nematic orientational field, which forms skyrmionlike twisted threedimensional (3D) structures called torons^{28,29,30}. Topological defects of higher than unit charge have not been observed before in nematic LCs, as an equivalent number of smaller charges usually achieves a lower free energy due to a lower elastic distortion^{31,32}.
NLCs are characterized by spontaneous orientational ordering of long axes of rodlike molecules, which is combined with complete positional disorder of the centres of gravity of molecules. This orientation order is described with a tensorial orientational field, where the largest eigenvalue of this Qtensor is pointing into a direction called the director, also referred to as a headless vector n. If the constituent molecules are chiral, the nematic phase spontaneously twists along a direction perpendicular to the director and forms twisted or chiral nematic liquid crystalline phase. When such a spontaneously twisted LC is confined into a sphere, which enforces perpendicular orientation of the LC molecules on the interface, the orientational field is frustrated and can form a variety of different topological states. It has recently been demonstrated experimentally that these states are rich in topological point defects^{30}, where the core of the defect is molten and the orientational order is strongly depressed. Theory also predicted the existence of linked and knotted loop defects in chiral nematic droplets^{33}, but such topological states have not yet been confirmed.
Here we show that longlived topological states in thermally quenched chiral nematic droplets are formed from stringlike, triangular and polyhedral constellations of monovalent and polyvalent singular point defects. Topological defects are regularly packed into a spherical liquid volume and stabilized by the elastic energy barrier due to the helical structure and confinement of the LC in the microsphere. We observe for the first time the quantization of topological charge of 3D point defects in terms of a unit topological charge. In addition to the q=−1 hedgehogs, which are usually stable in LC droplets and colloids, we observe higher charges with a multiple hedgehog charge q=−2, and −3. These monster charges can be regarded as polyvalent artificial atoms, and are able to bind the surrounding unit charge defects into polyhedral constellations representing topological molecules.
Results
Quenched chiral nematic droplets
To obtain droplets with various topological constellations of point defects we mix chiral LCs with chiral pitch in the range 6 to 12 μm into a viscous liquid, which promotes perpendicular anchoring to form droplets with diameters 10 to 20 μm (see Methods for details). After the droplets are formed, the sample is heated to the isotropic phase and cooled at a rate of several degrees per second back to the chiral nematic phase. This produces topologically complex and higher free energy director structures, which are stable for several days^{33,34}. The complexity of the topological states in chiral nematic droplets depends on the chirality parameter N=2d/p_{0}, where d is the diameter of the droplet and p_{0} is the intrinsic pitch of the chiral nematic LC. By far the most common are layered cholesteric structures with a single point defect, which is due to the topology of spherical confinement, but for N>2.5 structures with more point defects can appear. The droplets are examined by widefield optical and FCPM, and director fields are reconstructed from FCPM data using a recently developed method^{30} (details are described in Methods).
In samples with a chirality N∼2.5–4 we obtain two +1 radial hedgehogs residing at the surface of the droplet and facing each other on a symmetry axis, as shown in Fig. 1a–h. Their topological charge amounts to +2 and another hyperbolic hedgehog with negative charge −1 appears in between, adding up the total topological charge inside the sphere to +1. This is because the PoincaréHopf theorem requires the net hedgehog charge q of the director field n(r) with perpendicular surface anchoring along the closed bounding surface S with genus g to be equal to (up to sign) q=±(1−g) (refs 31, 32). For a sphere with g=0, having no handles, there should be a total hedgehog charge of +1, where we choose the sign according to the convention that the director field, represented as a vector field, points outwards from the surface of the droplet.
The space between the +1 and −1 hedgehogs is filled with a twisted LC, having the characteristic shape of an elongated cholesteric bubble, anchored with one end to the +1 defect, as shown in Fig. 1f. The cholesteric bubbles can be easily recognized in a streamline representation of the director field in Fig. 1g which shows only the projection of the director field on a plane. The white areas where the director is mostly perpendicular to the crosssection indicate the locations of the cholesteric bubbles. The bubbles are nonsingular, axially symmetric, smoothly winding structures, which induce elastic repulsion between the attracting pairs of topological charges and stabilize the defect constellations. The crosssection of a cholesteric bubble is a twodimensional skyrmion^{29} of the Bloch type^{35}, as shown in Fig. 1h.
At a higher chirality parameter N>3, we observe in different droplets many diverse constellations for the same value of N, which can be classified into two species: (i) strings of alternating point defects with unit hedgehog charges, (ii) constellations with a mixture of unit and multiple (that is, higher) hedgehog charges. In all cases, the +1 hedgehogs are found next to the surface, each anchored to its corresponding cholesteric bubble, while the negatively charged hedgehogs arrange in the bulk between the bubbles.
Stringlike constellations
Strings of alternating unitcharge hedgehogs are presented in Fig. 1i–o for 5 and in Fig. 2a,b and Fig. 2c,d for 7 and 9 unitcharge hedgehogs, respectively. Figure 1i shows a nonpolarized micrograph of a droplet with a chirality parameter N=3.1 and Fig. 1j–m show the corresponding experimental FCPM data. The reconstructed director field in Fig. 1n reveals a Vshaped constellation of five unit, chargealternating, point defects and the streamline representation in Fig. 1o indicates the positions of three cholesteric bubbles between the point defects as is schematically shown in the inset. The three hedgehogs close to the surface are of the +1 radial type, whereas the two hedgehogs in the midst of both arms of the Vshaped constellation are the −1 hyperbolic type, so the total topological charge is +1. Similarly as in the droplet with three point defects (Fig. 1a–h), the cholesteric bubbles, anchored to the +1 hedgehogs, face the hyperbolic −1 hedgehogs between them.
Even longer strings with 7 and 9 unit charge hedgehogs are found at higher chirality (Fig. 2, Supplementary Movies 1 and 2), always containing an odd number of hedgehogs, with each added pair of +1 and −1 hedgehogs forming another bubble and extending the string of hedgehogs with a new arm as shown in the schematic representations in insets to Fig. 2b,d. Whereas the 5 unitcharge hedgehogs in Fig. 1i–o can be neatly packed into the equatorial plane of the droplet, the constellations with 7 and 9 unit hedgehogs in Fig. 2b,d pack in stringlike structures in 3D.
Higher topological charges
Some hidden states are not formed only of unit topological charge hedgehogs, but also contain double (q=−2) and triple (q=−3) hedgehog charges as the binding elements. A double hedgehog charge appears in a planar configuration in Fig. 3a–c and Supplementary Movie 3 with a threefold optical symmetry containing three +1 hedgehogs close to the surface of the droplet and the fourth hedgehog (cyan dot) located exactly in the centre of the droplet. The three nearsurface +1 hedgehogs together with the central hedgehog must add up to the +1 total topological charge of the droplet, which means that the central hedgehog has a −2 topological charge. Its threefold symmetry of binding to the surface unit hedgehogs (Fig. 3b) resembles a trivalent atom with three symmetrically positioned orbitals (for example, a carbon atom with sp^{2} hybridized orbitals). This structure does not match any of the previously experimentally observed point defects in 3D nematics.
An even larger topological charge, carrying a hedgehog charge q=−3, is shown in Fig. 3d,e and Supplementary Movie 4. Here four unitcharge q=+1 hedgehogs form a tetrahedral constellation around the tetravalent topological charge in the centre of the droplet (green dot). Its topological charge q=−3 follows from the conservation of the overall topological charge. The q=−3 hedgehog forms a tightly squeezed director structure at the centre of the droplet, surrounded by four neighbouring cholesteric bubbles, which are compressing it towards the centre. The structure of this q=−3 hedgehog defect is fundamentally threedimensional, reminiscent of a sp^{3} hybridized atomic orbital.
Like in molecular chemistry, a single charged unit can be replaced by a compound structure with the same charge, while not changing the valence of the configuration. For example, in Fig. 4a,b and Supplementary Movie 5, one of the +1 hedgehogs from Fig. 3c is replaced by a stringlike constellation of three hedgehogs, similar to the threehedgehog string in Fig. 1a–g and likewise carrying a +1 topological charge. The tetravalent configuration from Fig. 3d,e also has a less symmetric variant, shown in Fig. 4c,d and Supplementary Movie 6 with six +1 hedgehogs in the vertices of an octahedron. Three of them are a part of the tetravalent configuration with a q=−3 point defect at the centre. The other three, together with two q=−1 hedgehogs, form a Vshaped string with a total charge q=+1, which acts as the fourth +1 hedgehog of the tetrahedral configuration. The stability of these structures demonstrates that the highcharge defects do not occur strictly in a confinement field with just the right symmetry, but can also withstand exchanges of the surrounding hedgehogs with larger constellations.
Discussion
The experiments shown in Figs 3 and 4 and Supplementary Movies 3–6 constitute the firstever observation of higher point topological charges in 3D. In NLCs, the topological charge is measured by counting the number of patches of the director field that are piercing the sphere which is enveloping the defect^{36}. Figure 5a,b show these patches, formed by the reconstructed director field around the −2 defect, which is decorated continuously with arrows; a patch is a region where the director field points outwards (red regions) on a continuous background of inwardpointing director (blue regions; note that this choice of arrows is opposite to the one we chose for the calculation of topological charge and its sign to make patches visually more discernible). Each patch can be understood as a bundle of hyperbolic director streamlines in Fig. 5c that converge towards a neighbouring +1 hedgehog. A defect with M patches has a hedgehog charge of q=1−M (ref. 36); therefore, the defect in Fig. 5a,b,d,e has a q=−2 charge and Fig. 5f shows a q=−3 charge. The patch structure reveals that q=−2 and q=−3 defects are generalizations of the hyperbolic hedgehog q=−1, and form at the confluence of streamlines from the +1 charges below the droplet surface, similar to how saddle points occur in an electric field formed by a set of equally charged point charges.
Highly symmetric constellations are obtained for the higher charges, which reside in the centre of the droplet: the q=−2 charge stabilizes the triangular molecule in Fig. 5g, the q=−3 charge binds the tetrahedral molecule in Fig. 5h, and the q=−5 the octahedral one in a model in Fig. 5i. The valence properties of the higher topological charges in chiral nematic spheres are similar to the tetrahedral structure of the monovalent topological defects in nematic shells^{37} which were proposed to be used for realization of tetravalent colloids^{38}.
It is quite interesting that we find the +1 hedgehogs close to the surface in all our experiments on topological hidden states in cholesteric droplets. This can be understood by considering the growth of the cholesteric regions (the cholesteric bubbles) after the quench. It is energetically favourable to expel the +1 hedgehogs towards the surface, where the helical order cannot exist because of the n·(∇ × n)=0 surface constraint. The homeotropic surface itself naturally prefers a general radiallike structure of the nearsurface defects, which is why the defects at the surface all have the same form and the same unit charge +1. Defects with positive topological charge greater than +1 cannot form close to the interface—they are not equivalent to hyperbolic saddle points, which occur naturally between other defects.
In quasitwodimensional systems, the higher topological charges are, in most cases, unstable, as an equivalent number of smaller defects or disclinations achieves a lower free energy due to a lower elastic distortion, but were observed in specific experiments^{39,40,41}. It was shown theoretically that in 3D achiral nematics the hedgehogs with higher charges than ±1 are unstable^{42}. Here they are stabilized by chirality (Fig. 6), which acts as a stabilizing spring inside the spherical confinement, preventing dissociation of the highercharge defects. As we can see from the stability chart Fig. 6, several different structures can appear at a given relative chirality N. This implies that the state which will emerge in a droplet after the quench is not determined by the helical pitch and droplet size; instead, the metastable states appear with probabilities which depend on their relative energies^{17,33}. The most complex hidden states appear only in the middle of the stability range of structures with the same number of positive charges (Fig. 6). This shows there is an optimal ratio between chirality and confinement, which promotes the formation of higher charge topological defects.
The higher charge defects could in general be true point defects, or have the internal structure of a defect loop, analogous to recently demonstrated defect loops in the cores of +1 defects^{43}, but with a more complex structure of its crosssection^{44}. Alternatively, higher charge defects could be composed of several tightly packed topological point defects, for example, the q=−2 could be composed of three −1 defects and one +1 defect. However, the FCPM image in Fig. 3b shows that the q=−2 defect is confined to a volume comparable to a −1 defect. Any internal structure would have to be confined to a volume smaller than the resolution of our microscope, which is ∼300 nm. This is an order of magnitude smaller than the helical period of the LC mixture, which sets the minimum separation between the defects for which the helical twist still acts as a repulsive force^{30}. We can, therefore, safely conclude that the central topological defect is a single point defect, carrying the topological charge q=−2.
Our experiments in chiral nematic droplets clearly reveal, for the first time, the full anatomy of the hedgehog defects, as we are able to reconstruct the topology of the 3D director and perform a patch analysis of the streamlines. Higher topological charges, which were previously considered unstable are in fact strongly stabilized by the chirality and spherical confinement of the NLC droplets on a time scale of several days. Surprisingly, the topological defects can be arranged in a perfect polyhedral constellation inside a liquid sphere, which has implications far beyond the field of LCs. Our strategy, which combines spherical confinement and a chiral field, offers guidance for the formation of hidden skyrmionic constellations in chiral magnets, cold atoms, and polyvalent softmatter colloidal science. Regular packing of the topological defects into a sphere can be regarded as the synthesis of topological molecules, formed of polyvalent defects, where the skyrmionlike structures have the role of topological bonds.
From a practical perspective, chiral nematic droplets with higher topological charges in the centre could serve as polyvalent colloidal particles, providing directional bonds to the surrounding colloidal particles. For example, a q=−3 topological charge in the centre of the chiral nematic droplet provides a tetravalent coordination of unit charge q=+1 hedgehogs located close to the surface of the droplet. If the surroundings of these q=+1 hedgehogs are decorated with polymer linkers, similar to the idea of Nelson for the nematic shells^{38}, a cholesteric droplet would be similar to a 4valent atom like carbon, silicon and germanium. Here the polymer linkers could provide artificial bonds to neighbouring droplets for selfassembly into colloidal crystals with a diamond lattice.
Methods
Preparation of droplets
The LC droplets were prepared with a lowbirefringence LC mixture of 1:1 weight ratio of 4′butyl4heptylbicyclohexyl4carbonitrile (CCN47) and 4,4′dipentylbicyclohexyl4carbonitrile (CCN55, both purchased from Nematel; refractive indices of the mixture n_{o}=1.47 and n_{e}=1.50) doped with 1–2% of chiral dopant S811 (Merck) to get a mixture with a pitch in the range ∼6–12 μm. A small amount of dye N,N′bis(2,5ditertbutylphenyl)3,4,9,10perylenedicarboximide (BTBP, Sigma Aldrich) was added to the LC mixture to enable FCPM by first dissolving the dye in acetone, adding the dye/acetone solution to the LC mixture and evaporating the solvent at room temperature to achieve a homogeneous mixture. Droplets with diameters in the range ∼10–20 μm were produced by mixing the LC mixture in a glycerol medium (n_{glycerol}=1.47) doped with 4% wt. Lαphosphatidylcholine (lecithin, Avanti Polar Lipids) to achieve perpendicular orientation of LC molecules on the droplet surface. The medium with the droplets was sandwiched between a microscopic cover glass of 150 μm thickness and a thicker 1 mm glass, which were separated by 30 μm mylar spacers. The cell was sealed along the perimeter with a fast curing twocomponent epoxy glue to prevent droplet movements because of glycerol flow. Constellations of point defects formed after the sample was quenched from the isotropic phase at a rate of several K s^{−1} to room temperature.
FCPM microscopy and director reconstruction
The structures in the droplets were imaged by fluorescence confocal polarizing microscopy (FCPM)^{45} on a Leica TCS SP5 X confocal microscope with a Leica WLL laser light source. The microscope was modified by inserting a quarterwave plate for the selected excitation wavelength (488 nm) in the slot below the objective to transform the polarization of the microscope from linear to circular. Linear polarizers were inserted above the waveplate into the same slot as the waveplate to select different linear polarizations of excitation. A full xyzscan was preformed at four linear polarizations in the xy plane, separated by 45° to obtain intensities I_{0}, I_{π/4}, I_{π/2} and I_{3π/4}. The scan at the first polarization was repeated in the end, to be used for linear correction of intensities because of bleaching between the scans. The xyzstacks were then deconvolved with SVI Huygens Professional software. Orientation of the director in the xy plane was calculated from the intensities as^{21}:
Total intensity I_{exp} was calculated as: I_{exp}=(I_{0}+I_{π/4}+I_{π/2}+I_{3π/4}) and zdependent background I_{offset} and normalization I_{norm} corrections were applied to calculate the outofplane angle θ from equation^{30}: . The director field was calculated from θ and φ as: n=(cos θ cos φ, cos θ sin φ, ±sin θ). A simulated annealing algorithm^{30} was used to determine the sign of the z component of the director field by randomly selecting a point in the droplet and calculating the elastic energy of director deformation in that point from:
where L is the elastic constant in the one constant approximation and q=2π/p_{0} is the inverse cholesteric pitch, which makes the elastic energy chiralitydependent. The tensor Q is calculated from the director as
where the scalar order parameter S, which is not available from the FCPM experiment, is taken to be constant. This elastic energy is compared with the elastic energy in that point with flipped zcomponent of director field (n_{x}, n_{y}, −n_{z}). If the energy of the flipped state is lower than the starting energy, this state is kept, but if it is higher, the flip is accepted with a probability given by a Boltzmann weight exp(−ΔE/t), where ΔE is the energy difference between the two states and t is a free parameter. The algorithm starts at a large t, iterating the flipping procedure over the whole volume of the droplet until a stable energy of the whole structure is reached, and then t is lowered and the procedure is repeated until thermalization of the structure is achieved. This procedure scans all possible configurations of signs of n_{z} and finds the ones with the lowest elastic energy within the restraint of experimental data for θ and φ, which means it is the best approximation of the observed structure that can be found from experimental data^{30}.
To generate the streamline representations of the director, the director field in a crosssection of the sample was projected on the crosssection plane and smoothed to indicate the average local direction. Areas where the projection of the normalized director is smaller than 1/3 (where director is mostly perpendicular to the crosssection) are excluded from the final streamline presentation.
Data availability
All the data relevant to the findings of this study are available from the authors upon request.
Additional information
How to cite this article: Posnjak, G. et al. Hidden topological constellations and polyvalent charges in chiral nematic droplets. Nat. Commun. 8, 14594 doi: 10.1038/ncomms14594 (2017).
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
 1
Hasan, M. & Kane, C. Colloquium: topological insulators. Rev. Mod. Phys. 82, 3045–3067 (2010).
 2
Qi, X.L. & Zhang, S.C. Topological insulators and superconductors. Rev. Mod. Phys. 83, 1057–1110 (2011).
 3
Burkov, A. Topological semimetals. Nat. Mater. 15, 1145–1148 (2016).
 4
Xu, S.Y. et al. Discovery of a Weyl fermion semimetal and topological Fermi arcs. Science 349, 613–617 (2015).
 5
Liu, Z. K. et al. Discovery of a threedimensional topological Dirac semimetal, Na3Bi. Science 343, 864–867 (2014).
 6
Neupane, M. et al. Observation of a threedimensional topological Dirac semimetal phase in highmobility Cd3As2 . Nat. Commun. 5, 3786 (2014).
 7
Hall, D. S. et al. Tying quantum knots. Nat. Phys. 12, 478–483 (2016).
 8
McConney, M. E. et al. Topography from topology: Photoinduced surface features generated in liquid crystal polymer networks. Adv. Mater. 25, 5880–5885 (2013).
 9
Dennis, M. R., King, R. P., Jack, B., O’Holleran, K. & Padgett, M. J. Isolated optical vortex knots. Nat. Phys. 6, 118–121 (2010).
 10
Keber, F. C. et al. Topology and dynamics of active nematic vesicles. Science 345, 1135–1139 (2014).
 11
DeCamp, S. J., Redner, G. S., Baskaran, A., Hagan, M. F. & Dogic, Z. Orientational order of motile defects in active nematics. Nat. Mater. 14, 1110–1115 (2015).
 12
Doostmohammadi, A., Adamer, M. F., Thampi, S. P. & Yeomans, J. M. Stabilization of active matter by flowvortex lattices and defect ordering. Nat. Commun. 7, 10557 (2016).
 13
Yeomans, J. M., Doostmohammadi, A. & Thampi, S. Defectmediated morphologies in growing cell colonies. Phys. Rev. Lett. 117, 048102 (2016).
 14
Nikkhou, M. et al. Lightcontrolled topological charge in a nematic liquid crystal. Nat. Phys. 11, 183–187 (2015).
 15
Nikkhou, M., Škarabot, M. & Muševič, I. Topological binding and elastic interactions of microspheres and fibres in a nematic liquid crystal. Eur. Phys. J. E 38, 23 (2015).
 16
Nikkhou, M., Škarabot, M. & Muševič, I. Annihilation dynamics of topological monopoles on a fiber in nematic liquid crystals. Phys. Rev. E 93, 062703 (2016).
 17
Tkalec, U., Ravnik, M., Čopar, S., Žumer, S. & Muševič, I. Reconfigurable knots and links in chiral nematic colloids. Science 333, 62–65 (2011).
 18
Jampani, V. S. R. et al. Colloidal entanglement in highly twisted chiral nematic colloids: Twisted loops, Hopf links, and trefoil knots. Phys. Rev. E 84, 031703 (2011).
 19
Martinez, A. et al. Mutually tangled colloidal knots and induced defect loops in nematic fields. Nat. Mater. 13, 258–263 (2014).
 20
Senyuk, B. et al. Topological colloids. Nature 493, 200–205 (2012).
 21
Chen, B. G.g., Ackerman, P. J., Alexander, G. P., Kamien, R. D. & Smalyukh, I. I. Generating the Hopf fibration experimentally in nematic liquid crystals. Phys. Rev. Lett. 110, 237801 (2013).
 22
Koshihara, S.Y., Tokura, Y., Mitani, T., Saito, G. & Koda, T. Photoinduced valence instability in the organic molecular compound tetrathiafulvalenepchloranil (TTFCA). Phys. Rev. B 42, 6853 (1990).
 23
Takubo, N. et al. Persistent and reversible alloptical phase control in a manganite thin film. Phys. Rev. Lett. 95, 017404 (2005).
 24
Stojchevska, L. et al. Ultrafast switching to a stable hidden quantum state in an electronic crystal. Science 344, 177–180 (2014).
 25
Kibble, T. W. Topology of cosmic domains and strings. J. Phys. A 9, 1387–1398 (1976).
 26
Chuang, I., Durrer, R., Turok, N. & Yurke, B. Cosmology in the laboratory: defect dynamics in liquid crystals. Science 251, 1336–1342 (1991).
 27
Zurek, W. H. Cosmological experiments in condensed matter systems. Phys. Rep. 276, 177–221 (1996).
 28
Smalyukh, I. I., Lansac, Y., Clark, N. A. & Trivedi, R. P. Threedimensional structure and multistable optical switching of tripletwisted particlelike excitations in anisotropic fluids. Nat. Mater. 9, 139–145 (2010).
 29
Ackerman, P. J., Trivedi, R. P., Senyuk, B., van de Lagemaat, J. & Smalyukh, I. I. Twodimensional skyrmions and other solitonic structures in confinementfrustrated chiral nematics. Phys. Rev. E 90, 012505 (2014).
 30
Posnjak, G., Čopar, S. & Muševič, I. Points, skyrmions and torons in chiral nematic droplets. Sci. Rep. 6, 26361 (2016).
 31
Mermin, N. D. The topological theory of defects in ordered media. Rev. Mod. Phys. 51, 591–648 (1979).
 32
Alexander, G. P., Chen, B. G.g., Matsumoto, E. A. & Kamien, R. D. Colloquium: disclination loops, point defects, and all that in nematic liquid crystals. Rev. Mod. Phys. 84, 497–514 (2012).
 33
Seč, D., Čopar, S. & Žumer, S. Topological zoo of freestanding knots in confined chiral nematic fluids. Nat. Commun. 5, 3057 (2014).
 34
Orlova, T., Asshoff, S. J., Yamaguchi, T., Katsonis, N. & Brasselet, E. Creation and manipulation of topological states in chiral nematic microspheres. Nat. Commun. 6, 7603 (2015).
 35
Yu, X. Z. et al. Realspace observation of a twodimensional skyrmion crystal. Nature 465, 901–904 (2010).
 36
Čopar, S. & Žumer, S. Topological and geometric decomposition of nematic textures. Phys. Rev. E 85, 031701 (2012).
 37
LopezLeon, T., Koning, V., Devaiah, K., Vitelli, V. & FernándezNieves, A. Frustrated nematic order in spherical geometries. Nat. Phys. 7, 391–394 (2011).
 38
Nelson, D. R. Toward a tetravalent chemistry of colloids. Nano Lett. 2, 1125–1129 (2002).
 39
Madhusudana, N. & Pratibha, R. Studies on high strength defects in nematic liquid crystals? Mol. Cryst. Liq. Cryst. 103, 31–47 (1983).
 40
Lavrentovich, O. D. & Nastishin, Y. A. Defects in degenerate hybrid aligned nematic liquid crystals. Europhys. Lett. 12, 135 (1990).
 41
Lavrentovich, O. D. & Pergamenshchik, V. Patterns in thin liquid crystal films and the divergence (‘surfacelike’) elasticity. Int. J. Mod. Phys. B 9, 2389–2437 (1995).
 42
Brezis, H., Coron, J.M. & Lieb, E. H. Harmonic maps with defects. Commun. Math. Phys. 107, 649–705 (1986).
 43
Wang, X. et al. Experimental insights into the nanostructure of the cores of topological defects in liquid crystals. Phys. Rev. Lett. 116, 147801 (2016).
 44
Čopar, S. & Žumer, S. Quaternions and hybrid nematic disclinations. Proc. R. Soc. A 469, 20130204 (2013).
 45
Smalyukh, I. I. & Lavrentovich, O. D. Threedimensional director structures of defects in GrandjeanCano wedges of cholesteric liquid crystals studied by fluorescence confocal polarizing microscopy. Phys. Rev. E 66, 051703 (2002).
Acknowledgements
G.P. acknowledges the financial support of Slovenian Research Agency (ARRS) through contract PR05014. S.Č. acknowledges the financial support of ARRS through contracts Z16725 and P10099. I.M. acknowledges the financial support of ARRS through contracts J16723 and P10099.
Author information
Affiliations
Contributions
G.P. conducted the experiments, G.P. and S.Č. analysed the results, I.M. conceived and supervised the experiments and wrote the main manuscript. All authors contributed to the final version of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Movie 1
Droplet with a string constellation of 7 point defects. Experimental FCPM intensities in a xy plane for the four excitation/detection polarisations (left half of the screen, the polarisation direction is marked with an arrow in top right corner of each panel) and the corresponding director field, reconstructed from the experimental data (right half of the screen). Positions of point defects are marked with yellow (q = +1) and magenta (q = 1) spheres, connected with rods to help visualise the spatial organisation of the defects. The droplet is scanned from the bottom along the z direction, the z coordinate is in micrometres. The scale bar is 5 μm. (MOV 3909 kb)
Supplementary Movie 2
Droplet with a string constellation of 9 point defects. Experimental FCPM intensities in a xy plane for the four excitation/detection polarisations (left half of the screen, the polarisation direction is marked with an arrow in top right corner of each panel) and the corresponding director field, reconstructed from the experimental data (right half of the screen). Positions of point defects are marked with yellow (q = +1) and magenta (q = 1) spheres, connected with rods to help visualise the spatial organisation of the defects. The droplet is scanned 34 from the bottom along the z direction, the z coordinate is in micrometres. The scale bar is 5 μm. (MOV 5325 kb)
Supplementary Movie 3
Droplet with a q = 2 point defect. Experimental FCPM intensities in a xy plane for the four excitation/detection polarisations (left half of the screen, the polarisation direction is marked with an arrow in top right corner of each panel) and the corresponding director field, reconstructed from the experimental data (right half of the screen). Positions of point defects are marked with yellow (q = +1) and cyan (q = 2) spheres, connected with rods to help visualise the spatial organisation of the defects. The droplet is scanned from the bottom along the z direction, the z coordinate is in micrometres. The scale bar is 5 μm. (MOV 3756 kb)
Supplementary Movie 4
Droplet with a q = 3 point defect. Experimental FCPM intensities in a xy plane for the four excitation/detection polarisations (left half of the screen, the polarisation direction is marked with an arrow in top right corner of each panel) and the corresponding director field, reconstructed from the experimental data (right half of the screen). Positions of point defects are marked with yellow (q = +1) and green (q = 3) spheres, connected with rods to help visualise the spatial organisation of the defects. The droplet is scanned from the bottom along the z direction, the z coordinate is in micrometres. The scale bar is 5 μm. (MOV 4661 kb)
Supplementary Movie 5
Tetrahedral constellation in a droplet with a q = 2 point defect. Experimental FCPM intensities in a xy plane for the four excitation/detection polarisations (left half of the screen, the polarisation direction is marked with an arrow in top right corner of each panel) and the corresponding director field, reconstructed from the experimental data (right half of the screen). Positions of point defects are marked with yellow (q = +1), magenta (q = 1) and cyan (q = 2) spheres, connected with rods to help visualise the spatial organisation of the defects. The droplet is scanned from the bottom along the z direction, the z coordinate is in micrometres. The scale bar is 5 μm. (MOV 4431 kb)
Supplementary Movie 6
Octahedral constellation in a droplet with a q = 3 point defect. Experimental FCPM intensities in a xy plane for the four excitation/detection polarisations (left half of the screen, the polarisation direction is marked with an arrow in top right corner of each panel) and the corresponding director field, reconstructed from the experimental data (right half of the screen). Positions of point defects are marked with yellow (q = +1), magenta (q = 1) and green (q = 3) spheres, connected with rods to help visualise the spatial organisation of the defects. The droplet is scanned from the bottom along the z direction, the z coordinate is in micrometres. The scale bar is 5 μm. (MOV 6501 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Posnjak, G., Čopar, S. & Muševič, I. Hidden topological constellations and polyvalent charges in chiral nematic droplets. Nat Commun 8, 14594 (2017). https://doi.org/10.1038/ncomms14594
Received:
Accepted:
Published:
Further reading

Oblique light incidence method to study topological defects in nematic layers with conical boundary conditions
Scientific Reports (2021)

Selfassembly of fractal liquid crystal colloids
Nature Communications (2019)

Topologydependent selfstructure mediation and efficient energy conversion in heatfluxdriven rotors of cholesteric droplets
Nature Communications (2018)

Topologycommanded optical properties of bistable electricfieldinduced torons in cholesteric bubble domains
Scientific Reports (2017)

Bipolar configuration with twisted loop defect in chiral nematic droplets under homeotropic surface anchoring
Scientific Reports (2017)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.