## Abstract

Emulsions comprising isotropic fluid drops within a nematic host are of interest for applications ranging from biodetection to smart windows, which rely on changes of molecular alignment structures around the drops in response to chemical, thermal, electric, and other stimuli. We show that absorption or desorption of trace amounts of common surfactants can drive continuous transformations of elastic multipoles induced by the droplets within the uniformly aligned nematic host. Out-of-equilibrium dynamics of director structures emerge from a controlled self-assembly or desorption of different surfactants at the drop-nematic interfaces, with ensuing forward and reverse transformations between elastic dipoles, quadrupoles, octupoles, and hexadecapoles. We characterize intertransformations of droplet-induced surface and bulk defects, probe elastic pair interactions, and discuss emergent prospects for fundamental science and applications of the reconfigurable nematic emulsions.

## INTRODUCTION

Emulsions are common soft matter systems associated with various foods and drinks, health and personal care products, bodily fluids essential for functioning of humans and animals, and many other (*1*). Nematic emulsions comprise nematic liquid crystals (LCs) as at least one of the fluid emulsion components (*2*, *3*). They have recently attracted a great deal of applied and fundamental interest because of bringing about facile response of the LC to highly enrich the physical behavior of these systems (*4*–*18*). Nematic drops within isotropic fluids and nematic host materials with embedded isotropic droplets exhibit fascinating structures of defects and distortions in the LC molecular alignment field (*2*–*18*), which comply with the constraints of topological theorems (*13*). Such structures of molecular alignment within or around the droplets can serve as reporters of biological microbes or chemical sensors (*4*, *5*) and can drive various forms of reconfigurable elasticity-mediated self-assembly (*2*, *7*, *11*–*13*). In emulsions comprising foreign fluid drops within an aligned nematic host, a variety of defect and structures have been reported, giving origin to particle-induced deformations of the far-field alignment that range from elastic dipoles (*2*) to quadrupoles (*7*, *16*), to hexadecapoles, and to even higher-order multipoles (*17*–*19*). Some of these structures, which typically can be characterized by leading-order elastic multipoles, were found intertransforming in response to applied electric fields or changes of temperature and other experimental conditions (*7*, *18*). Elasticity-mediated self-assembly of the droplets within the nematic host medium is highly dependent on the elastic multipolar nature of the director (**n** ≡ −**n**) distortions around the droplets, as well as on the structure and geometry of the corresponding particle-induced topological defects, and could be used in designing mesostructured composite materials (*2*–*20*). However, limitations of controlling the occurrence of different elastic multipoles and corresponding defect structures hinder the utility of LC emulsions in fundamental science and technological applications.

Nematic emulsions are of great fundamental interest. How a sphere with well-defined boundary conditions interacts with a vector or director field is understood within the predictions of topological theorems, such as the hairy ball theorem (*21*). However, topological theorems define the rules for the field-surface interactions only up to the topological characteristics so that a number of topologically analogous but energetically similar structures can arise while having very different structures of particle-induced defects (*2*, *3*, *22*, *23*). In the context of LC colloids and emulsions, this manifests itself via various interesting topology-compliant director configurations appearing under different conditions, with a classic example being formation of either a hedgehog point defect or a Saturn ring disclination that both allow embedding a drop with perpendicular boundary conditions into a uniformly aligned nematic background (*2*, *8*, *13*, *16*). Transformations between some of these elastic multipole structures have been demonstrated, including when prompted by external fields (*7*, *18*, *24*, *25*), although these studies so far were limited to only few types of elastic multipoles. Mechanisms of transformations between various topologically admissible configurations of defects and elastic multipoles as transiently changing or stable structures are known for only limited types of surface boundary conditions (*7*), even though fascinating nematic emulsions and colloids with patchy surface anchoring have been demonstrated (*26*–*30*).

Here, we report spontaneous transformations between elastic dipoles, quadrupoles, octupoles, and hexadecapoles occurring in different forward and reverse sequences in response to absorption or desorption of surfactant molecules at the LC-isotropic fluid interfaces of droplets. While the elastic dipoles and quadrupoles are commonly studied (*2*, *8*, *13*) and high-order elastic multipoles have been reported to occur at very specific predesigned boundary conditions on surfaces of solid particles (*31*, *32*), our study reveals how these different multipoles can arise within an emergent hierarchical self-assembly process occurring at length scales from ~1 nm to tens of micrometers. These highly controlled transformations of elastic multipoles can occur in a forward and reverse order, in agreement with our computer simulations based on Landau–de Gennes free energy minimization, when the system is driven to different equilibria corresponding to different types of behavior exhibited by surfactant molecules at the LC-isotropic fluid interfaces. We demonstrate how a systematic variation of patchy surface boundary conditions through controlling the surfactant molecule self-assembly triggers fusion, splitting, and various other transformations of singular nematic defects at surfaces and in the bulk of the nematic host medium. We show that the ensuing formation and transformation of various elastic multipoles can progress under out-of-equilibrium conditions, being driven by absorption and self-assembly of surfactant molecules, or the reverse process. The newly observed elastic multipoles of different orders and signs, each with a characteristic configuration of droplet-induced topological defects, may enable new forms of colloidal self-assembly (*13*), nanoparticle templating by defects (*33*), and biodetection (*4*), whereas their fascinating intertransformation sequences can be used to design programmable reconfigurability of the mesostructured nematic emulsions.

## RESULTS

### Droplet-induced elastic multipoles

Far away from droplets and both bulk or surface defects induced around the isotropic droplet surfaces (at distances much larger than the droplet dimensions), the director distortions with respect to the uniform far-field background can be represented as a summation over the elastic multipoles (*31*, *32*)*l* determines the order of a multipole as the 2* ^{l}*-th pole, −

*l ≤ m ≤ l*, and

*r*

_{0}is a radius of the droplet. One can immediately notice the analogy between the multipole expansion of the droplet-induced director distortions in the far-field and the corresponding expansion of electrostatic charge distributions (

*31*,

*32*,

*34*). This analogy stems from the fact that the nematic LC elastic free energy density within the one-elastic-constant approximation reads

*K*is the average Frank elastic constant and

*n*

_{μ}(μ =

*x*,

*y*) are director components perpendicular to the far-field director

**n**

_{0}(along the

*z*direction). We note that the use of this one-constant approximation in this case is crucial to the electrostatics analogy and is justified by the fact that the different elastic constants often have comparable values (

*34*). The structures of elastic colloidal multipoles minimize this free energy density that gives Euler-Lagrange equations in the form of the Laplace equations, bringing about the electrostatic analogy (

*34*)

Multipole moments *n _{x}*,

*n*<< 1, and

_{y}*n*≈ 1 far from the particle surface. While elastic colloidal dipoles (

_{z}*l*= 1) and quadrupoles (

*l*= 2) in emulsions and solid particle dispersions have been widely studied for over two decades (

*2*,

*13*), higher-order multipoles such as hexadecapoles (

*l*= 4) were observed only recently and typically require rather special shape or surface functionalization to exist (

*17*,

*18*,

*31*,

*32*). While symmetry-allowed higher-order multipoles are always present in the above expansion, the question often is if they need to be considered, which depends on relative strengths of different multipole moments (

*13*). Our findings below show how self-assembly and desorption of surfactant molecules at LC-glycerol interfaces drive transformations between highly diverse director structures that, up to leading order, can be described as multipoles with

*l*= 1, 2, 3, and 4.

### Topological transformations from boojum quadrupole to hedgehog dipole

Micrometer-sized spherical droplets of isotropic fluid in a nematic LC were obtained by mixing a small amount of glycerol (~1 to 5 volume %; fig. S1A) with a nematic LC and vigorous agitation via a combination of flicking and sonication (see Materials and Methods). These droplets maintain their spherical shape when mixed with the LC because of the high interfacial tension. The emulsions are infiltrated into glass cells of about 30- to 60-μm gap thickness, much larger than the diameter of the droplets. To satisfy the boundary conditions at the droplet surface, the LC director field **n**(**r**) deviates from the far-field director **n**_{0} defined by boundary conditions at the confining glass plates and forms characteristic patterns with accompanying topological defects. These different director perturbations are characterized by the ensuing elastic colloidal multipoles (*13*, *31*, *32*, *34*). When pure glycerol is used to form droplets dispersed in 4-cyano-4′-pentylbiphenyl (5CB) or ZLI-2806 nematic LCs, planar azimuthally degenerate boundary conditions are naturally promoted at the LC-glycerol interface (*3*). This results in a quadrupolar configuration of **n**(**r**) around the droplet (Fig. 1A, a) with two surface point defects (called “boojums”) at its poles on diametrically opposite sides of the droplet along **n**_{0} [Fig. 1, B (a) and C (a)]. Such droplets form an elastic quadrupole with *l* = 2 and *m* = +1 [Fig. 1, A (a) to D (a)] (*7*, *31*), which we will refer to as a boojum quadrupole. Boojum quadrupoles obtained with a pure glycerol are stable over long periods of time without changing the **n**(**r**)-configuration. However, adding a trace amount of surfactant to glycerol before mixing with the LC alters boundary conditions and changes this droplet-induced elastic multipole behavior.

Surfactant SDS (fig. S1B) uniformly dissolves within glycerol at studied small concentrations. At SDS concentrations >0.1 weight % (wt %), homeotropic boundary conditions are induced at droplet surfaces, giving the origins to a dipolar **n**(**r**)-configuration with a bulk point defect (Fig. 1A, h). For SDS at <0.1 wt %, the droplets initially appear as boojum quadrupoles with tangential boundary conditions [Fig. 1, A (a) to C (a)] right after the preparation, but **n**(**r**)-configurations and singular defects change with time. First, two boojums at the poles open into the quarter-integer surface-bound defect loops encircling the surface of the droplet at its poles. Later, these surface defect loops slowly drift away from the poles [Fig. 1A (a to d) and movie S1] and, upon reaching the equator, combine to form a single half-integer bulk defect loop called “Saturn ring,” as well as the corresponding quadrupolar **n**(**r**)-configuration with *l* = 2 and *m* = −1 (Fig. 1A, e) (*16*, *31*). Then, the Saturn ring starts moving toward one of the poles of the droplet (Fig. 1A, f and g), eventually transforming into a bulk point defect of elementary hedgehog charge; the corresponding dipolar **n**(**r**)-configuration (Fig. 1A, h) has *l* = 1 and *m* = +1. Polarizing and bright-field microscopy textures (Fig. 1, A and B, respectively) reveal details of these transformations, with the corresponding schematics and numerical visualizations of elastic multipoles shown in Fig. 1 (C and D). Note that, during the first phase of the transformations (Fig. 1A, b to d), when two quarter-integer surface defect loops are moving toward the equator, the polarizing micrographs feature eight bright lobes around the droplet (fig. S2A), as in the elastic hexadecapoles with conically degenerate surface anchoring (fig. S2E) studied recently (*17*–*19*). These experimentally observed topological transformations of elastic multipoles are caused by tangential-to-homeotropic changes of patchy surface anchoring boundary conditions at the LC-glycerol interface (Fig. 2, A to C). SDS molecules are initially evenly dispersed within the volume of the glycerol droplet (Fig. 2A), but with time, they diffuse to and self-organize at the LC-glycerol interface and cause the surface anchoring transition. Nucleation of surfactant monolayer islands with homeotropic anchoring typically starts at the poles of the droplet (Fig. 2B) because the ensuing local change of boundary conditions minimizes the corresponding elastic energy. Self-assembled surfactant islands with homeotropic boundary conditions at the poles open the boojums into quarter-integer surface defect loops and, with time, transform elastic quadrupoles into hexadecapoles with mixed anchoring [Figs. 1B (b) and 2B], different from hexadecapoles studied previously (*17*, *18*). Increasing the surfactant density leads to the growth of these patches with homeotropic boundary conditions, expanding dynamics of the surface defect loops that then move toward the equator (Fig. 1B, b to d). As the two patches with homeotropic boundary conditions meet at the equator and two quarter-integer surface disclination loops merge to form a half-integer bulk disclination loop, the elastic hexadecapole with mixed anchoring transforms into an elastic quadrupole with a Saturn ring [Figs. 1A (e) and 2C]. The transformation from a Saturn ring quadrupole to an elastic dipole is nearly twice faster than the transformation from a boojum quadrupole to a Saturn ring quadrupole (Fig. 2D). Figure 3 shows the corresponding results of numerical computer simulations, where boundary conditions are continuously changed from tangential to homeotropic starting from the poles, mimicking the experimental observations (Materials and Methods). Computer-simulated polarizing optical microscopy textures of transformations of the director structure around the droplet and corresponding elastic multipoles (Fig. 3 and movies S2) are in excellent agreement with the experimental observations.

By controlling the amount of the added surfactant, one can stop topological transformations at different desired stages with a specific director configuration around the droplet and a resulting elastic multipole. Although topological transformations from a boojum quadrupole to an elastic dipole are rather interesting, even more exciting are the intermediate states that correspond to the observed elastic hexadecapoles [Figs. 3A (b and c) and 4] and octupoles [Figs. 3D (b to d) and 5] formed by droplets with mixed patchy tangential-homeotropic surface anchoring. The growth of self-assembled surfactant islands is an out-of-equilibrium process that enriches accessible structures during transformation by prompting the appearance of states that cannot be found during conventional switching transitions, say when induced by electric fields (*7*, *25*). As we discuss in detail below, the observed transformations of elastic multipoles conserve topological charges of the induced defects like boojums, hedgehogs, and Saturn rings while they all compensate the effective topological charge due to uniform or mixed boundary conditions on the spherical droplet’s surface (*7*, *13*, *21*, *34*).

### Elastic hexadecapoles with mixed patchy anchoring

The configuration of **n**(**r**) around the droplet intermediate between a boojum quadrupole and a Saturn ring quadrupole is hexadecapolar [Figs. 3A (b and c) and 4 and fig. S2]. A polarizing optical micrograph (Fig. 4A) shows eight bright lobes around the perimeter of the droplet separated by narrow dark regions within which **n**(**r**) is parallel to polarizer or analyzer. Director field deviates away from **n**_{0} in the regions of bright lobes. In addition, insertion of the phase retardation plate into the polarizing microscope’s optical path reveals that the tilt of **n**(**r**) with respect to **n**_{0} is different within adjacent bright lobes (Fig. 4B) and switches between clockwise (blue color in the texture) and counterclockwise (yellow color in the texture) directions eight times along the perimeter. These textural features and **n**(**r**) are similar to the ones exhibited by an elastic hexadecapole (fig. S2, E to H) previously found forming around colloidal particles with conically degenerate surface anchoring (*17*–*19*). However, the **n**(**r**)-structure and defects in the hexadecapoles with tilted anchoring are different from that in hexadecapoles with mixed anchoring observed here. First, bright-field textures of hexadecapoles with mixed anchoring (Fig. 4C) reveal the presence of two quarter-integer surface defect loops, which separate areas with tangential and normal boundary conditions [Figs. 1A (b and c) and 3A (b and c)], whereas hexadecapoles with conically degenerate anchoring have two boojums at the poles and one defect loop at the equator [compare Fig. 4 (A to D) and fig. S2 (E to G)]. Second, the hexadecapolar structure observed in the current study has four bright lobes corresponding to the areas with tangential boundary conditions, and four bright lobes close to the poles correspond to the areas with homeotropic boundary conditions of the droplet. Experimental textures (Fig. 4, A to C) reveal the detailed **n**(**r**) around the droplet (Fig. 4D), which is in good agreement with numerical modeling (Fig. 3A, b or c) (*17*–*19*, *31*). A color map of the projection *n _{x}* of

**n**(

**r**) onto the

*x*axis orthogonal to

**n**

_{0}highlights positive, near-zero, and negative

*n*[Figs. 3B (b and c) and 4E] around the droplet and is also consistent with a hexadecapolar

_{x}**n**(

**r**)-configuration (

*17*–

*19*,

*31*,

*32*). Mapping of

*n*at some distance from the spherical surface reveals that the tilting of

_{x}**n**(

**r**) with respect to

**n**

_{0}in the hexadecapole with mixed anchoring has an opposite sequence as compared to that in the hexadecapole with conic anchoring [compare Fig. 4E, fig. S2 (B and D), and fig. S2 (F and H)]. Thus, both hexadecapoles have the same

*l*= 4 leading multipole order but differ in that we have

*m*= +1 for the hexadecapoles with conic anchoring and

*m*= −1 for the hexadecapoles with mixed anchoring (

*31*).

Sharing of the far-field deformations of the director between different elastic multipoles gives rise to elasticity-mediated interactions. At center-to-center distances much larger than the droplet size, the colloidal pair interaction potential, derived from an electrostatic analogy of the far-field distortions (*17*), reads*R* is a distance between the centers of two interacting droplets, *a _{l}* =

*b*

_{l}r_{0}

^{l + 1}is an elastic multipole moment of the

*l*th order,

*P*

_{l + l′}(cosθ) are the Legendre polynomials, and θ is an angle between

**n**

_{0}and the center-to-center droplet separation vector

**R**. Although Eq. 5 is accurate only at large distances, it gives qualitative insights into the pair interaction potential behavior even at distances corresponding to only few particle sizes, although singular defects and near-field effects start playing important roles when

*R*is becoming comparable to the droplet diameter. The θ-dependent elastic pair interactions of hexadecapoles have eight zones of attraction separated by eight zones of repulsion (

*17*–

*19*), which correlate with eight bright lobes in polarizing optical micrographs. Two droplets attract if zones of the same color (Fig. 4B) in the polarizing textures with a phase retardation plate radially face each other. Pair interactions of hexadecapoles with mixed anchoring depend on relative areas of patches with dissimilar anchoring, with quadrupolar interactions dominating when their size is very different [Fig. 3, A (b) and B (b)]. The quadrupolar moment is at minimum and hexadecapole moment is maximized when the sizes of the bright (Fig. 4A) or blue and yellow (Fig. 4B) lobes along the droplet perimeter in polarizing textures are comparable [see also Fig. 3 (A and B, c)]. For example, the ratio between the hexadecapolar (

*b*

_{4}) and quadrupolar (

*b*

_{2}) moments obtained from fitting (Fig. 4F) is

*b*

_{4}/

*b*

_{2}≈ 1.93 at θ ≈ 38°. While the far-field director structure and interactions of such droplets are consistent with their leading-order hexadecapole nature, note that the near-field effects and proximity of singular defects alter this behavior at small distances.

### Elastic octupoles with mixed anchoring

Another interesting **n**(**r**)-configuration around the droplet and corresponding elastic multipole can be obtained when one of the boojums opens into the growing surface defect loop but the other remains at the pole (Fig. 5 and movies S3). In our experiments, because of the associated nucleation and growth process, such asymmetric evolution of structures often happens spontaneously, but this process can be controlled if a small foreign inclusion is trapped by one of the boojums but not the other, thus resulting in different energetic barriers for the surfactant island self-assembly at the two poles. In this case, polarizing micrographs (Fig. 5A) show only six bright lobes around the droplet perimeter. Tilt of **n**(**r**) with respect to **n**_{0} is different within adjacent bright lobes (Fig. 5D) and switches between clockwise (blue color in the texture) and counterclockwise (yellow color in the texture) directions with respect to **n**_{0}, like in the hexadecapole with mixed anchoring, but only six times. The bright-field micrograph (Fig. 5B) reveals the presence of one quarter-integer surface disclination in the lower droplet’s hemisphere and the boojum defect at the top pole. The experimentally reconstructed **n**(**r**) around such a droplet shown in Fig. 5C is in good agreement with numerical calculations (Fig. 3D, b to d). A color map of the projection *n _{x}* of

**n**(

**r**) onto the

*x*axis orthogonal to

**n**

_{0}(Fig. 5E) closely resembles the electrostatic octupole (

*13*,

*31*), although here, its elastic octupolar counterpart is formed by induced director distortions around a colloidal inclusion undergoing anisotropic diffusion with different coefficients (

*D*

_{||}/

*D*

_{⊥}= 1.35) measured parallel and perpendicular to

**n**

_{0}(fig. S3, G and H) (

*35*,

*36*).

Two elastic octupoles attract (Fig. 5F) if zones of the same color in the polarizing micrographs with a retardation plate (Fig. 5D) align radially more or less against each other, so that the angular dependence of their elastic pair interactions has six zones of attraction separated by six zones of repulsion (Fig. 5H and inset of Fig. 5G). This behavior correlates with the angular diagram of elastic interactions of two antiparallel octupoles (Fig. 5H and inset of Fig. 5G) calculated using elastic moments (*b*_{1},*b*_{3},*b*_{5}) = (0.182, 0.499, −0.005), obtained from fitting the potential (Fig. 5G) extracted from the experimental attraction between two droplets (Fig. 5F). We note that the dependence described by Eq. 5 is, strictly speaking, expected to hold only at large distances, but it provides qualitative insights even at *R* corresponding to several droplet sizes; here, like for hexadecapoles discussed above, singular defects and near-field effects modify the short-range elastic interactions when droplets closely approach each other. The angular dependence of interactions (Fig. 5H) will change for parallel octupoles: Zones of attraction will become zones of repulsion and vice versa. The elastic interaction potential of octupoles amounts to hundreds of *k*_{B}*T* (Fig. 5G) and is an order of magnitude weaker than that for elastic dipoles and quadrupoles formed by colloidal inclusions of comparable dimensions (*13*). The strength of the octupolar moment depends on the ratio between the area of patches with homeotropic boundary conditions and the area of patches with tangential boundary conditions between the equator and a surface defect (Figs. 5, A to E, and 3, D to F). The quadrupolar moment is dominant when the area with homeotropic boundary conditions is small, and a dipolar moment increases when the surface defect line is close to the equator and homeotropic boundary conditions cover nearly half of the droplet. The octupolar moment is at maximum (*b*_{3}/*b*_{1} ≈ 2.7) when the quarter-integer surface defect line is somewhere in the midway between the pole and the equator (Fig. 5, D and F).

The elastic octupole can transform into the elastic dipole if there are enough SDS molecules to cover the whole LC-glycerol interface (Fig. 5I) and change boundary conditions to homeotropic all over the droplet’s surface. In this case, the quarter-integer surface disclination loop formed by opening the boojum at the bottom pole travels all the way from the bottom pole to the top pole, where it merges with the boojum surface point defect, which, in turn, produces a bulk hedgehog point defect.

### Topological transformations from hedgehog dipole to boojum quadrupole

Topological transformations in the opposite direction, from an elastic dipole to a boojum quadrupole, also occur within the inverse LC emulsions. To demonstrate them, as well as to show that the studied transformations are rather universal for LCs and not specific to a certain nematic compound or a mixture, we used a nematic 5CB as a host medium for the glycerol droplets with trace amounts of water, different from a multicomponent nematic mixture studied above. An alignment agent *N*,*N*-dimethyl-*N*-octadecyl-3-aminopropyl-trimethoxysilyl chloride (DMOAP) (fig. S1C) was then used as a surfactant mixed with glycerol at <0.1 wt %, sufficient to initially cover the entire droplet’s surface by the surfactant monolayer. Figure 6A shows the sequence of bright-field textures during this transition. Polarizing micrographs reveal that a hyperbolic hedgehog point defect at one of the poles and the ensuing dipolar **n**(**r**)-configuration right after the sample preparation (Fig. 6A, a), when strong homeotropic boundary conditions (Fig. 6B, a) are defined by DMOAP molecules at the interface. With time, however, this dipolar configuration becomes unstable and the hedgehog opens into the half-integer disclination loop, the Saturn ring, which then expands and moves toward the equator [Fig. 6A (b) and movie S4]. This transformation of the hedgehog into the Saturn ring configuration can be explained by weakening of the homeotropic surface anchoring, which might be caused by a slow dissociation of DMOAP molecules from the interface (*37*). While full understanding of the nature of this anchoring weakening at the interface will require separate detailed studies, a possible reason is that the trace amounts of water within glycerol cause DMOAP molecules to dimerize/polymerize and be less effective in inducing the homeotropic interfacial boundary conditions for the director and even eventually to dissociate from the interface. The change of the surface anchoring strength alters the interplay between surface and bulk elastic free energy contribution, causing the out-of-equilibrium structural transformations that allow the system to reduce the overall free energy via transformation to the quadrupolar structure (*24*, *25*, *38*). After the Saturn ring of a half-integer disclination reaches the equator, it splits into two quarter-integer surface disclinations moving toward the droplet poles [Fig. 6A (c and d) and movie S4]. The Saturn ring splits into two disclinations when boundary conditions at the equator locally change from homeotropic to tangential, being a reverse process of what we described above during the surfactant self-assembly causing the growth of homeotropic anchoring patches. The two quarter-integer surface disclinations separate regions of dissimilar, tangential and homeotropic, boundary conditions and move toward poles as an area with tangential boundary conditions continuously grows out from the near-equator region. The area with tangential boundary conditions nucleates at the equator because a change of boundary conditions from homeotropic to tangential decreases director deformations [**n**(**r**) becomes parallel to **n**_{0}] at the equator and corresponding elastic energy. The loops of quarter-integer surface disclinations eventually contract into two boojums [Fig. 6, A (e) and C, and movie S4] after they reach the poles of the droplet. These transformations of elastic multipoles from an elastic dipole to a boojum quadrupole undergo with conservation of topological charges and become possible because of the change of boundary conditions at the LC-glycerol interface from homeotropic to tangential. The patchy surface anchoring transition likely happens because of hydrolysis of silane groups and polymerization of individual DMOAP molecules initially located at the LC-glycerol interface but eventually desorbing from it. This patchy out-of-equilibrium anchoring transition drives transformation of multipolar structures that is remarkably similar (but in opposite direction) to that of a boojum quadrupole to a hedgehog dipole transformation discussed above. Numerical modeling once again reproduces the entire transformation sequence while additionally revealing the final stage of transformation when boojums inherit the half-ring handle-like geometry (Fig. 6C and movie S5), similar to that previously studied for solid colloidal particles with different genus (*22*, *39*).

### Fusion, splitting, and other transformations of singular defects

Structures of dipolar and quadrupolar director distortions, with accompanying hedgehog, Saturn ring, and boojum defects, are widely studied and well understood [Fig. 1C (a, e, and h)] (*2*, *3*, *13*, *22*, *23*). In all these cases, the droplet-induced defects appear to match well-defined tangential or homeotropic boundary conditions at the LC-droplet interface to the uniform far-field alignment. A spherical glycerol droplet has a genus equal to zero and Euler characteristic χ_{Euler} = 2. Consistent with topological theorems (*2*, *3*, *13*, *22*, *23*), the hedgehog charge of the three-dimensional topological defects π_{2}(^{2}/ℤ_{2}) = ℤ induced in the nematic LC bulk [the hyperbolic hedgehog point defect and the equivalent to it Saturn ring of π_{1}(S^{2}/ℤ_{2}) = ℤ_{2} disclination] that compensates the hedgehog charge of surface boundary conditions at the drop surface is equal ±χ_{Euler}/2 = ±1; the sign of the induced defect can be determined upon decorating the nonpolar director field with a unit vector field along one of the two antiparallel directions (*23*, *40*). The net winding number of π_{1}(S^{1}) LC defects in the two-dimensional interfacial director field at the LC-glycerol interface is equal to χ_{Euler} = 2, with the two boojums each having a winding of unity [Figs. 1C (a) and 7]. When the core of a boojum splits into a handle-like disclination loop (Fig. 7A), the interfacial director field then contains four defects with half-integer winding numbers. The intermediate structures of topological defects are subject to the same topological theorem constraints, albeit now for the mixed tangential-homeotropic patchy boundary conditions (*2*, *3*, *13*, *22*, *23*). An important feature of the intermediate structures observed during transformations is the presence of fractional surface defect lines. Unlike in the bulk, the winding number of disclinations is not constrained to be an integer multiplied by ±^{1}/_{2} so that fractional disclinations are also allowed (*41*). The ones observed in our study are quarter-integer disclinations with the winding number of −^{1}/_{4} (Figs. 1C, 4D, 5C, and 6B), where this winding number is found by locally circumnavigating the surface defect line and measuring the director rotation within the LC in the proximity of this defect, divided by 2π. While past studies demonstrated junctions of such quarter-integer surface defect lines and the bulk half-integer defect lines (*41*), our current study shows how two quarter-integer surface disclinations can combine to form a single loop of the bulk half-integer disclination π_{1}(S^{2}/ℤ_{2}) = ℤ_{2} (Fig. 1C), as well as that a reverse process can take place (Fig. 6B) while preserving conservation of topological characteristics. The forward-directed transformations of elastic multipoles involve fusion of quarter-integer surface disclinations into a single half-integer bulk defect loop (Fig. 1C, c to e), and its splitting into two quarter-integer surface disclination loops is observed during the reverse process (Fig. 6B, b to d).

Another interesting aspect related to topological defects is associated with the nature of nucleation of surfactant monolayer islands (Fig. 7, A to N). Consistently with experiments, the numerical modeling shows that the free energy minimization of the LC host medium favors nucleation of surfactant islands at the droplet poles, at the location of point-like or split-core handle-like surface boojums (Fig. 7O), although even the islands that nucleate away from poles are strongly attracted to the poles. The diverse set of possible structures in our system includes boojums with point-like cores and handle-like split cores and quarter-integer disclination loops around surfactant islands replacing the boojums or coexisting and interacting with them (Fig. 7, A to N). A local minimum and a nonmonotonic behavior of free energy versus displacement is observed for the islands of surfactant with homeotropic patch of anchoring nucleated away from the poles; this behavior arises as a result of interaction and fusion of distinct singular cores of a boojum defect and a disclination loop surrounding the island. While the structures of cores of defects with locally reduced scalar order parameter look more complex when surfactant islands nucleate away from the top/bottom poles (Fig. 7, D, E, I, M, N, and inset of O), their net charge in the interfacial director field still adds to χ_{Euler} = 2, with a unity winding number per top/bottom hemisphere. Well-defined defect structures are important for controlling spatial localization of various nanoparticles and even self-assembled molecular structures (*13*), which makes our finding particularly interesting because the observed transformations of defect structures could be potentially used for reconfiguration of defect-entrapped nanostructures.

## DISCUSSION

Our results indicate that self-assembly of surfactant molecules can drive transformation of defects and director structures surrounding isotropic liquid droplets within LC emulsions formed from both single-compound LCs such as 5CB and various nematic mixtures. Our study can be extended to explore the differences between such transformations in emulsions with single- and multicomponent LCs, where various forms of interfacial phase separation could play a role in the latter case. Our observations for large droplets of size *r*_{0}≫ *K*/*W* are consistent with the regime of strong polar anchoring imposed by patches of the droplet surface with or without surfactant monolayers. The polar anchoring surface boundary conditions for micrometer-sized drops and *K* ~ 10^{−11} N and *W* ~ 10^{−5} J/m^{2} can be considered as strong because the minimization of the total free energy mainly occurs through relaxation of the director structure within the LC bulk while keeping the director at LC-droplet interfaces fixed along the spatially varying easy orientation determined by the patchy distribution of surfactant molecules at the LC-droplet interface. One can envisage even richer behavior occurring as the droplet size would be reduced first to the regime of soft boundary conditions, where surface and bulk terms could be comparable and compete, and then to the regime of small drops when the surface energy would be relaxed while keeping bulk director distortions at minimum. The latter regime could be analogous to what was already observed for LC nanodroplets (*42*) so that the LC surrounding could drive patterning of surfactants within the droplet surfaces without elastic director distortions around the droplet. Isotropic fluid droplets, surfactant monolayers, and the surrounding LC could all be polymerized by fixing respective configurations through cross-linking and could be decorated with various nanoparticles, potentially offering new routes for hierarchical assembly of structured materials. When using microfluidics, droplet dimensions and amounts of accessible surfactants could be controlled rather precisely, potentially allowing for fabrication of emulsions with monodisperse droplet sizes, surfactant coverage, and types of elastic multipoles induced by the droplets within the emulsions (*13*). This could allow for probing elasticity-mediated colloidal self-assembly due to high-order elastic multipoles and formation of various mesostructured composite materials enabled by self-assembly of reconfigurable elastic colloidal multipoles. Nematic emulsions with out-of-equilibrium structural transformations could be also engineered for various biodetection applications (*4*, *14*, *15*), where different elastic multipoles and defect structures could potentially serve as reporters of biological molecules and microbes at the LC-fluid interfaces (note that water could be used instead of glycerol).

## MATERIALS AND METHODS

### Samples preparation and imaging techniques

Two different LCs were used as hosts for the glycerol (Sigma-Aldrich; fig. S1A) droplets: 5CB (from Frinton Laboratories Inc.) and ZLI-2806 (EMD Performance Materials Corp.). A small amount (<0.1 wt %) of surfactant had been mixed with glycerol to alter a surface anchoring for LC molecules at the interface between LC and glycerol droplets. We mixed glycerol with an SDS surfactant (fig. S1B) for droplets in ZLI-2806 and an alignment agent DMOAP (fig. S1C) for droplets in 5CB. To obtain micrometer-sized spherical droplets in a nematic host, the glycerol mixture was added to the LC and then vigorously stirred using a combination of flicking and sonication. The ensuing diluted LC emulsions with spherical droplets (*r*_{0} ≈ 1 to 10 μm) suspended evenly in a LC host were infiltrated into ~30- to 60-μm-thick cells made of two glass plates separated by glass spacers and sealed with a ultraviolet-curable glue. A polyimide PI2555 (HD Microsystems) was spin-coated on the glass plates, baked at 270°C, and unidirectionally rubbed with a velvet cloth to promote a tangential alignment of the LC in the cell. Experimental samples were stable over at least several days, and LC structures around droplets were studied by means of bright-field and polarized light optical microscopy with an inverted Olympus IX81 microscope and a ×100 (numerical aperture, 1.42) oil objective. Optical trapping and manipulations of the droplets were performed by holographic optical tweezers operating at 1064 nm (*43*). Microscopy image acquisition and analysis were performed using a charge-coupled device camera (Flea, Point Grey) and ImageJ software, respectively.

### Details of computer simulations

Numerical continuum simulations were performed to study dynamics of nematic tensorial order parameter **Q**. For uniaxial systems, the nematic order parameter is written in the form of **Q** = *S*(**nn** − **I**/3), where unit vector **n** is the nematic director field and *S*, the largest eigen value of **Q**, quantifies the degree of uniaxial alignment. Within nematodynamics framework, the evolution of nematic tensorial order parameter follows the Beris-Edwards formalism (*44*)**S** is the generalized advective term that accounts for the response of the nematic order parameter to the symmetric and antisymmetric parts of velocity gradient. Note that the experimental observations reveal the quasi-static nature of nematic director evolution that this slow dynamic justifies, simplifying the hydrodynamic equations and neglecting the reorientation of the director with backflow effects. Hence, Eq. 6 reduces to the Ginzburg-Landau equation

The molecular field

The free energy density, *f*_{LdG}, is the sum of bulk (describing isotropic/nematic phase transition) and elastic energies (accounting for penalizing nematic distortion from the uniform configuration) given by

The relaxation rate is controlled by the collective rotational diffusion constant Γ. The phenomenological coefficient *A*_{0} sets the energy scale, *U* controls the magnitude of the order parameter, and *L* is the elastic constant in the one-elastic constant approximation. At the droplet surface with unit normal **ν**, the anchoring condition is imposed by adding a surface term to the free energy*V*. Note that hybrid types of anchoring of the director field is enforced at the droplet surface. In a portion of the droplet with the homeotropic anchoring, the Rapini-Papoular–like surface free energy density of the form**Q**^{0} = *S _{eq}*(

**νν**−

**I**/3). The fourth-order Fournier-Galatola free energy density is adopted to apply the degenerate planar anchoring boundary condition on the rest of droplet surface

*W*controls the anchoring strength,

**p**=

**δ**−

**νν**. Assuming that all of the transitions are axisymmetric (i.e., rotationally symmetric with respect to the director far field), patchy boundary conditions with different regions having either tangential or homeotropic anchoring were imposed on the droplet surface. While mimicking experimental observations, these boundary conditions were gradually changed axisymmetrically by advancing the location of the quarter-integer disclination along the polar angle–changing direction.

The computer simulations were performed on a cubic mesh of 240 × 240 × 240 uniformly distributed points. The total free energy of the system was minimized, corresponding to the nulls of first variation of F, and the resultant Euler-Lagrange partial differential equations were solved using an explicit Euler relaxation finite difference scheme (*45*–*51*). The droplet radius was chosen to be *r*_{0} = 50.25 in simulation unit. Formation of the hyperbolic hedgehog point defect was ensured by initializing the director field with the ansatz of the following form (*46*)*P* = 1.5 determines the separation distance of the defect and the droplet center, and **d** and **d**_{0} are the position vector and the position of the center of mass of the droplet, respectively.

The optical polarization micrograph textures were simulated using the method of Jones matrices and experimental material parameters (*52*, *53*). In all of the simulated textures, the incident polarized wave vector propagated along the *y* axis (perpendicular to the image), while polarizer and analyzer were crossed with the polarizer aligned along the far-field director. The following numerical parameters have been used throughout the simulations: *A*_{0} = 6 × 10^{5} J m^{−3} and *U* = 3.0, yielding equilibrium value of *S _{eq}* = 0.5,

*L*= 6 pN, and the anchoring strength in the range of

*W*= 6 × 10

^{−5}to 6 × 10

^{−3}J m

^{−2}.

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/7/25/eabg0377/DC1

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is **not** for commercial advantage and provided the original work is properly cited.

## REFERENCES AND NOTES

**Acknowledgments:**We thank Q. Liu, B. Fleury, and J. B. ten Hove for the discussions.

**Funding:**This work was partially supported by Department of Energy, Office of Basic Energy Sciences, under contract DE-SC0019293 with the University of Colorado, Boulder (B.S., K.C., and I.I.S.). K.C. acknowledges financial support from NSF and the Soft Materials Research Center (SMRC) Research Experience for Undergraduates (REU) program at the University of Colorado, Boulder. The work of J.J.d.P., A.M., and R.Z. is supported by the Department of Energy, Basic Energy Sciences, Division of Materials Science and Engineering, grant DE-SC0019762. The development of software for simulation of morphology and dynamics in liquid crystals is supported by DOE under the Midwest Integrated Center for Computational Materials, MICCOM, at Argonne National Laboratory.

**Author contributions:**B.S. and K.C. conducted the experiments. A.M. and J.J.d.P. developed the models and performed the simulations. R.Z. contributed to the theoretical analysis. I.I.S. conceived and directed the project. B.S. and I.I.S. analyzed the results and wrote the manuscript, with input from all authors.

**Competing interests:**The authors declare that they have no competing interests.

**Data and materials availability:**All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

- Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).