Abstract
Coexistence of order and fluidity in soft matter often mimics that in biology, allowing for complex dynamics and applications-like displays. In active soft matter, emergent order can arise because of such dynamics. Powered by local energy conversion, this behavior resembles motions in living systems, like schooling of fish. Similar dynamics at cellular levels drive biological processes and generate macroscopic work. Inanimate particles capable of such emergent behavior could power nanomachines, but most active systems have biological origins. Here we show that thousands-to-millions of topological solitons, dubbed âskyrmionsâ, while each converting macroscopically-supplied electric energy, exhibit collective motions along spontaneously-chosen directions uncorrelated with the direction of electric field. Within these âschoolsâ of skyrmions, we uncover polar ordering, reconfigurable multi-skyrmion clustering and large-scale cohesion mediated by out-of-equilibrium elastic interactions. Remarkably, this behavior arises under conditions similar to those in liquid crystal displays and may enable dynamic materials with strong emergent electro-optic responses.
Similar content being viewed by others
Introduction
Soft matter and living systems are commonly described as close cousins1, both with properties stemming from interactions between the constituent building blocks that are comparable in strength to thermal fluctuations. Active soft matter systems2,3 are additionally out-of-equilibrium in nature, like everything alive. They exhibit emergent collective dynamics that closely mimic such behavior in living systems2,3,4. For example, mechanically agitated fluidized monolayers of rods form a dynamic granular liquid crystal (LC)5. Coherent motion emerges in many systems where particles communicate through collisions or short-range interactions like screened electrostatic repulsions6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22. However, these interactions typically cannot be tuned in strength or switched from attractive to repulsive. Moreover, with a few exceptions5,6,7,10,11,19,20,21,22, including the ones in which electric energy is used to power motions6,7,11, most active matter systems have biological origins and either chemical or mechanical energy conversion within the constituent building blocks. This poses the grand challenge to develop versatile reconfigurable active matter formed by inanimate, man-made particles both as models of biological systems and for technological uses21.
We describe an emergent collective dynamic behavior of skyrmions23,24,25, particle-like two-dimensional (2D) topological analogs of Skyrme solitons used to model atomic nuclei with different baryon numbers26,27. In LCs1,28, these skyrmions are elements of the second homotopy group29 and contain smooth but topologically nontrivial, spatially localized structures in the alignment field of constituent rod-like molecules, the director field n(r). They are characterized by integer-valued topological invariants, the skyrmion numbers29. Depending on the applied voltage, the internal n(r) structures within our skyrmions adopt different orientations relative to the 2D sample plane and the far-field alignment. Thousands to millions of skyrmions start from random orientations and motions while each individually converting energy due to oscillating voltage, but then synchronize motions and develop polar ordering within seconds. The ensuing schools of topological solitons differ from all previously studied systems2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22 because skyrmions have no physical boundaries, membranes, chemical composition or density gradients, or singularities in the order parameter at the level of the host fluid30,31,32,33,34,35,36, even though they exhibit giant-number fluctuations in terms of the skyrmionic, topologically protected, localized structures of n(r). Although the LC medium is nonpolar, spontaneous symmetry breaking and many-body dynamic interactions lead to polar ordering of skyrmions, which is characterized by near-unity values of polar and velocity order parameters. Electrically tunable interactions stemming from the orientational elasticity of LCs30,31 provide a versatile means of controlling this behavior while probing order and giant-number fluctuations within the schools. The dynamic multi-skyrmion assemblies echo formation of high-baryon-number skyrmions in nuclear physics due to the addition of charge-one topological invariants, which in equilibrium condensed matter could be only achieved when forming skyrmion bags29, very differently from the behavior of singular active matter topological defects that conserve the net winding number to always add to zero2,33,34,35,36. Our findings highlight the interplay between nonsingular topology of field configurations and out-of-equilibrium behavior and promise a host of technological uses.
Results
Skyrmion schooling
Schooling of fish (Fig. 1a), like many other forms of collective motions2, is accompanied by inhomogeneities and dynamic local clustering. Similar behavior is observed in our rather unusual schools formed by thousands to millions of localized particle-like skyrmionic orientational structures of n(r) within LCs (Fig. 1bâf). While moving and bypassing obstacles, these skyrmions exhibit dynamically self-reconfigurable assembly (Fig. 1f and Supplementary Movie 1). In our experiments, under conditions and sample preparation similar to that in LC displays, such skyrmions are controllably mass-produced at different initial densities (see Methods)32 and also generated one-by-one using laser tweezers23,29. Skyrmion stability is enhanced by soft perpendicular boundary conditions on the inner surfaces of confining glass plates (Fig. 1d)23 and the LCâs chirality, which prompts n(r) twisting23,24. At no fields, the structure of each skyrmion is axisymmetric (Fig. 1g)23, with Ï-twist of n(r) from the center to periphery in all radial directions and containing all possible n(r) orientations within it. Due to the used LCâs negative dielectric anisotropy, electric field E applied across the cell tends to align n(r)â¥E (Fig. 1d), so that n(r) around the skyrmions progressively tilts away from the cell normal with increasing voltage U, whereas n(r) within the skyrmions morphs from an originally axisymmetric structure (Fig. 1g) to a highly asymmetric one (Fig. 1j, k) that matches this tilted director surrounding. Since the tilting of n(r) in E with respect to the sample plane breaks the nonpolar symmetry of the resulting effectively 2D structure, we vectorize n(r) and visualize it with arrows colored by orientations and corresponding points in the two-sphere \({\Bbb S}^2\)-order parameter space (Fig. 1g)29. The asymmetric skyrmion is described by a preimage vector connecting preimages of the south and north poles of \({\Bbb S}^2\) (Fig. 1k), the regions where n(r) points into and out-of-the sample plane, respectively. The skyrmion number, a topological invariant describing how many times n(r) within the single skyrmion wraps around \({\Bbb S}^2\), remains equal to unity, indicating topological stability with respect to smooth n(r) deformations prompted by changing U. This is confirmed by experimental polarizing optical micrographs of individual skyrmions at different U (Fig. 1i, m) that closely match their computer-simulated counterparts (Fig. 1h, l). Skyrmions exhibit only Brownian motion at no fields and at high frequencies (like 1âkHz) of applied field, at which n(r) cannot follow the temporal changes of E and responds to its time average.24 At frequencies for which the voltage oscillation period TUâ=â1/f is comparable or larger than the LC response time, spatial translations of individual skyrmions24 in an oscillating electric field arise because the temporal evolution of asymmetric n(r) is not invariant upon reversal of time with turning U on and off within each TU (Figs. 1gâm and 2a, b). Because of the initial axisymmetric structure of skyrmions and homeotropic n(r) background around them, the spontaneous symmetry breaking leads to random motion directions of individual skyrmions within the sample plane. This feature of our system indicates that the coherent unidirectional motion of many skyrmions within the thousands-to-million schools is an emergent phenomenon with a physical mechanism that relies on inter-skyrmion interactions, which we explore in detail below.
Collective motions of two-to-hundred skyrmions
To gain insights into the emergent schooling with electrically reconfigurable clusters of topological solitons (Fig. 1gâm), we first probe how the dynamics of skyrmions change with increasing number density (Fig. 2). Using laser tweezers, we set up a âraceâ by arranging skyrmion pairs along a straight line together with single skyrmions (Fig. 2a) and then start oscillating U at fâ=â2âHz by effectively turning it on and off every 0.5âs while using a 1âkHz carrier frequency electric source. The single solitons move faster than pairs (Fig. 2a) and trio-chains (Supplementary Movie 2), whereas large assemblies of 40 and more skyrmions barely move (Supplementary Fig. 1). Since spatial translations of skyrmions arise from temporal evolution of asymmetric n(r) not invariant upon reversal of time with turning U on and off, skyrmions within large assemblies tend to share these asymmetric distortions to reduce the elastic free energy30,31, which impedes their translations while in tightly packed reduced-asymmetry assemblies. However, with the oscillating fieldâs modulation changing to fâ=â1âHz, the solitons tend to spread apart, so that the motion of well-spread pairs is then nearly as fast as that of single skyrmions (Fig. 2b). The dynamic evolution of n(r) during motions of skyrmion pairs is enriched by elastic interactions that tend to minimize the ensuing elastic free energy costs due to arranging skyrmions at different relative positions and can be controlled from attractive to repulsive (Fig. 2c). Tuning U and f alters these inter-skyrmion interactions and reconfigures larger kinetic assemblies of tens (Supplementary Movie 3) and hundreds of skyrmions (Supplementary Movie 4), which in turn alters their dynamics. For example, with tuning f within 2â8âHz, a cluster of 29 skyrmions shown in Supplementary Movie 3 re-forms into a long chain, which breaks into smaller clusters and then rearranges again and again, multiple times. Hundreds of skyrmions at initial packing fractions <0.01 (Fig. 2dâf) form chains meandering like snakes (Fig. 2dâf and Supplementary Fig. 2 and Movie 4) and bounce from each other within a dynamic cluster-like region. Inaccessible to skyrmions at equilibrium29, such self-reconfigurable behavior emerges in both right- and left-handed chiral LCs obtained with different chiral additives (see Methods). Moreover, over long periods of time, there is no repetition of assemblies and trajectories of the randomly directed motion, which we verify by analyzing the complex dynamics (Fig. 2f and Supplementary Movies 3 and 4). Although the polarizing video microscopy reveals how the asymmetric periodically changing n(r) evolution powers motions of multi-skyrmion assemblies (Fig. 2),24 this mechanism alone cannot explain self-reconfigurable randomly directed dynamics of two-to-hundreds skyrmion assemblies and schooling of thousands to millions of skyrmions with tunable clustering within the schools. Systematic analysis of pair interactions (Fig. 2c) reveals that this complex behavior can be understood by taking into account the out-of-equilibrium elastic interactions between skyrmions that arise from minimizing the elastic free energy due to partial sharing of n(r) distortions associated with multiple moving skyrmions, as we detail below.
Out-of-equilibrium elastic interactions
Elastic interactions between skyrmions emerge to reduce the free energy costs of n(r) distortions around these topological solitons, like in nematic colloids29,30,31,32, albeit typically without the dynamic n(r) fully reaching equilibrium because of the voltage modulation and soliton motions. The elastic interactions between skyrmions confined to a 2D plane have dipolar nature29, although the complex temporal evolution of n(r) in periodically modulated U makes these elastic dipoles effectively change their tilt relative to the 2D sample plane within TU and self-propel while they interact. Such dynamic dipolar skyrmions mutually repel at small U, but exhibit anisotropic interactions (Fig. 2c), including attractions, when oscillating E prompts their n(r) symmetry breaking (Fig. 1j, k) and motions. Oscillating E rotates preimage dipoles from pointing orthogonally to the sample plane at Uâ=â0 (Fig. 1g) to being tilted or in-plane (Fig. 1j, k) when U increases, with the effective tilt periodically changing with a voltage modulation period (Fig. 1d) typically comparable to the LCâs response time. When released at different relative initial positions using laser tweezers32, skyrmions with parallel preimage vectors perpendicular to substrates always repel, whereas skyrmions with in-plane preimage dipoles attract when placed head to tail and repel when side by side (insets of Fig. 2c). Depending on U, f, and relative skyrmion positions, the strength of reconfigurable elastic pair interactions (Fig. 2c) varies within (1â10,000)kBT, where kB is the Boltzmann constant and T is the absolute temperature. Since the response of n(r) to oscillating U is fast on the timescales of skyrmion motions at ~1âμm per second, tuning n(r) by U and f modifies elastic forces between parallel dipoles by changing the effective tilt (averaged over TU) of the dipole moments relative to the sample plane.
Emergence of polar order and coherent motions
In the presence of thousands to millions of skyrmions (Fig. 1bâf), applied E initially induces random tilting of the director around individual skyrmions, so that their southânorth preimage unit vectors piâ=âPi/|Pi | point in random in-plane directions (Fig. 3). Individual skyrmions exhibit translational motions with velocity vectors vi roughly antiparallel to their pi. With time, coherent directional motions emerge (Supplementary Movies 1 and 5â7), with schooling of skyrmions either individually dispersed (Fig. 4 and Supplementary Movie 5) or in various assemblies (Fig. 5 and Supplementary Movies 1, 6, and 7). Velocity and polar order parameters \(S = \left| {\mathop {\sum}\nolimits_i^N {{\mathbf{v}}_i} } \right|/(N{{\mathbf{v}}_s})\) and \(Q = \left| {\mathop {\sum}\nolimits_i^N {{\mathbf{p}}_i} } \right|/N\) characterize degrees of ordering of vi and pi within the moving schools19, where N is the number of skyrmionic particles and vs is the absolute value of velocity of a coherently moving school. Both S and Q increase from 0 to ~0.9 within seconds (Fig. 4d), indicating the emergence of coherent unidirectional motion of polar skyrmionic particles, like that of fish in schools2,3. At relatively low initial packing fractions (~0.1 by area), we observe no clustering of skyrmions as they move coherently within the schools, repelling each other at short distances and weakly attracting at larger distances (Fig. 4c). This emergent behavior is different from pair interactions and dynamics at similar voltages (Fig. 2c), where moving skyrmions tend to attract to form chains at shorter inter-skyrmion distances. The presence of such short-range repulsive and long-range attractive interactions is consistent with the formation of coherently moving schools and results from many-body interactions (Fig. 4c), where elastic interactions between skyrmions with periodically evolving n(r) are further enriched by backflows and electro-kinetic effects.28 As the elastic interactions vary from attractive to repulsive within TU, the effective time averaging of these interactions localizes skyrmionic particles at distances roughly corresponding to the distance at which pair interactions are comparable to kBT (Figs. 2c and 4a, c). The many-body interactions between skyrmions then lead to effective cohesion within the school and their coherent collective motion. Using video microscopy, we also analyze the mean <N> and root mean square ÎNâ=â<(Nâ<N>)2 >1/2 of particles within different sample areas (see Methods and Fig. 4e). Unlike in the case of random Brownian motion of colloidal particles or the same skyrmions with αâ=â0.5, when ÎNâââ<N>â1/2, skyrmions in schools exhibit giant-number fluctuations with ÎNâââ<N>α, where αâ=â0.763 (Fig. 4e), as well as fluctuations in the local number density probed by counting the numbers of skyrmions within a selected sample area versus time (Fig. 4f).
Tunable clustering, edges, and cohesion in skyrmion schools
We alter the skyrmion schooling behavior by inducing formations of clusters (Fig. 5). At moderately large 0.1â0.4 skyrmion packing fractions within the schools (Fig. 5), we observe motions of dynamically self-assembled clusters at Uâ=â2.5â3.75âV and linear chains at Uâ=â3.75â4.5âV (Fig. 5aâd and Supplementary Movies 1, 6, and 7). Like individual skyrmions in schools of lower density (Fig. 4), moving clusters and chains remain separated at distances ~30âμm corresponding to effective pairwise interactions ~kBT (Fig. 4c). During this schooling, small clusters and linear polar chains exhibit giant-number fluctuations with varying values of αâ=â0.61â0.85 (Fig. 5e). Within the clusters and chains, skyrmions are kept at separation distances comparable to their lateral size (Fig. 5aâm) and roughly consistent with the separation distances corresponding to minima of potentials of pair interactions at similar conditions (Fig. 2c), which can be tuned by U and f through tuning the temporal evolution of n(r), as we show for the case of chains in Fig. 5n. This dynamic assembly of multi-skyrmions echoes nuclear physics models, where subatomic particles with high-baryon numbers can be modeled as clusters of elementary skyrmions27. Each skyrmion cluster can be characterized by a net skyrmion number corresponding to a sum of topological invariants of elementary skyrmions within it (e.g., clusters in Fig. 5f, j have net skyrmion numbers of 7 and 5, respectively). Tuning packing fractions, U and f, allow emergence of a broad range of this collective behavior (Fig. 5). The direction of collective motion within inch-square cells (Figs. 1bâd, 4, and 5) is selected spontaneously and emerges only at sufficiently large number densities of skyrmions, although gradients of cell gap thickness and external fields could potentially be used to control it.
Edges of schools are well defined regardless of the internal clustering within the schools (Fig. 6aâd). Individual skyrmions, which happen to be slowed down by imperfections (Fig. 6e), move faster than the clusters and thus âcatch upâ to the schoolâs edge. Although the speed decreases as the number of skyrmions within the clusters increases, this reduction is <50% even for very large clusters containing over 100 skyrmions (Fig. 6f), also showing how the dynamic behavior under schooling conditions differs from that of individual skyrmions and their small-to-medium clusters (Supplementary Fig. 1).
Diagram of dynamic and static states
We summarize the schooling behavior of skyrmions using a structural diagram (Fig. 7), where we present results for square-wave electric fields oscillating at frequency f (we note that the diagram of states changes when other electric signal waveforms and various carrier frequencies are used, although exploration of all these parameter spaces is outside the scope of our present work). Being unstable at high U and low f (Fig. 7), skyrmions exhibit static self-assemblies at low U and high f and dynamic structures at intermediate U and f. The intermediate-strength E is needed to asymmetrically morph the axisymmetric skyrmions observed at low U, without destroying the skyrmions by the strong electric alignment taking place at high U. Using intermediate frequencies avoids various electro-kinetic instabilities at low f (at which skyrmions become unstable due to spatial re-distributions of ions that further alter the director field) while still allowing for out-of-equilibrium temporal evolution of n(r) not invariant under turning the instantaneous voltage on and off. This is because the electric field oscillation period TUâ=â1/f within this frequency range is comparable to the LCâs rising and falling response times (within 20â100âms for our samples). Formation of clusters and chains within schools at different voltages is consistent with the nature of out-of-equilibrium elastic interactions between skyrmions at oscillating U revealed by numerical modeling (Fig. 5fâm) and qualitatively discussed above by using the dipolar elastic interactions. Fine details of clusters, like inter-particle distances (Fig. 5n) vary along the f-axis and would be difficult to capture within a single diagram of states, but the simplified three-dimensional diagram in Fig. 7 overviews the tendencies and helps to emphasize the physical underpinnings of the observed rich out-of-equilibrium behavior, which we summarize below.
Discussion
Playing a key role in skyrmion schooling, many-body elastic interactions (Figs. 4 and 5) minimize the energetic costs of periodically varying n(r) within schools by tending to position individual skyrmions such that they share the dynamic distortions and reduce the overall free energy. Motion of clusters and chains is impeded as compared to the fastest individual skyrmions (Fig. 6e, f) because of the very same sharing of asymmetric dynamic distortions between individual skyrmions, the non-reciprocal evolution of which is the source of motion. During each TUâ=â10â20âms within schooling (Figs. 4 and 5), n(r) never fully relaxes because the LCâs 20â100âms response time is longer than TU, so that the skyrmions are always asymmetric with periodically changing preimage dipole tilts (Fig. 1j, k). When the instantaneous U within TU drops to zero, the elastic torque tends to relax the n(r) of all skyrmions to an axially symmetric state shown in Fig. 1g, but well before this happens, competing electric and elastic torques re-morph the skyrmions back to the highly asymmetric structures (Fig. 1j, k). While a viscous torque resists changes of n(r), torque balances are different in the presence of E and without it, making the responses to turning instantaneous U on and off highly asymmetric and non-reciprocal24. Consequently, asymmetric skyrmions translate within a tilted director background in response to oscillating voltage, roughly antiparallel to their pi, much like polar granular particles translate in response to mechanical vibrations19. Asymmetric skyrmions synchronize pi and vi even before colliding by sensing each other through the long-range many-body elastic interactions.29 Remarkably, this elasticity-enhanced synchronization can take place at packing fractions ~0.01 and inter-skyrmion distances ~10 times larger than the solitonâs lateral size, consistent with the long-range nature of elastic forces. Given that collective motions can arise at carrier frequencies ~1âkHz, at which ions are too slow to follow oscillating fields and their dynamics can be neglected, electro-kinetics is not a pre-requisite for the studied effects. Since numerical modeling reproduces motions of individual skyrmions24 and their chains and clusters when using only the rotational viscosity/torque (Fig. 5fâm and see Methods), it appears that flows are not essential for the collective dynamics of skyrmions, like in the âdryâ types of active matter2 (e.g., herds of cows and biological cells crawling on substrates), although such flows are locally present23. Backflow and electro-kinetic effects enrich the collective dynamics, although their detailed study and uses are beyond the scope of our current work. Skyrmion clusters also dynamically interact to exchange and re-arrange elementary skyrmions both spontaneously and during interactions with obstacles, such as other skyrmions pinned to substrates using laser tweezers (Fig. 8 and Supplementary Movie 1). Importantly, this local bypassing of obstacles does not alter the direction of schooling, but could potentially be a useful tool in controlling collective dynamics in schools of skyrmions, as well as could potentially be extended to other active matter systems37.
To conclude, we have demonstrated active matter formed by solitonic particle-like field configurations, with salient features of energy conversion at the individual particle level and synchronization of initially random motion directions that leads to skyrmion schools. This schooling displays voltage-controlled self-reconfigurations of coherent motions with and without clustering. While much of the recent excitement in active matter has been generated by topological defects, which exhibit fascinating dynamics33,34 and play key roles in living tissues35,36, our findings demonstrate that not just singular defects, but also topological solitons can behave like active particles. Skyrmion schooling can allow for modeling diverse forms of non-equilibrium behavior, benefiting from non-biological origins and on-demand creation/elimination of skyrmions using laser tweezers29 and providing insights into the role of topology and orientational elasticity in active matter. For example, it will be of interest to explore how giant-number fluctuations arise despite of and in the presence of the self-aligning nematic fluid hosts of skyrmions with orientational elasticity and long-range inter-skyrmion interactions. Our system comprises commercially available ingredients, with design and preparation techniques benefiting from display industry developments. It also connects the topology29 and active matter2 paradigms, potentially resulting in fertile new research directions at their interfaces. From a materials-applications perspective, one can envisage photonic and electro-optic materials, including displays and privacy windows, with built-in emergent responses capable of controlling light, effectively expanding potential exciting applications of more common active matterials2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,33,34,35,36. Being compatible with the touch-screen displays and related technologies, skyrmion schooling can be coupled to external stimuli responses and interactions with humans, potentially yielding active matter art and computer games invoking emergent behavior. Since these skyrmions can carry nanoparticle cargo25,32, their schooling can yield active self-reconfigurable metamaterials and nanophotonic devices. Our approach can be also extended to nematic colloids30,31, where one can potentially achieve electrically powered schooling of colloidal particles within the LC. These findings call for the development of novel active matter modeling approaches capable of handling collective behaviors of thousands to millions of schooling skyrmions, each with periodically morphing complex structures of the molecular alignment field and with temporal director evolution coupled to flows and electro-kinetic effects.
Methods
Sample preparation
Chiral nematic LC mixtures with negative dielectric anisotropies were prepared by mixing a chiral additive (CB15, ZLI-811, or QL-76) with a room-temperature nematic host (MLC-6609 or ZLI-2806). CB15 and ZLI-2806 were purchased from EM Chemicals. ZLI-811 and MLC-6609 were purchased from Merck. The QL-76 chiral additive38,39 was obtained from the Air Force Research Laboratory (Dayton, OH). Pitch, p, of the mixtures was controlled by varying the concentration, c, of the chiral additive with known helical twisting power (Table 1), hHTP, according to the relation pâ=â1/(hHTP·c). Studied samples had pâ=â3â10âµm and d/pâââ1.23 The LCs were mixed with ~0.1âwt% of cationic surfactant hexadecyltrimethylammonium bromide (CTAB, purchased from Sigma-Aldrich) in order to allow for inducing electrohydrodynamic instability with low-frequency applied field32. The CTAB doping allowed for facile generation of large numbers (thousands to millions) of skyrmions at different initial packing fractions upon relaxation of the cells from electrohydrodynamic instability (Fig. 9), although the presence of CTAB is not required for skyrmion schooling as similar dynamics could also be obtained in samples without CTAB and with skyrmions generated by laser tweezers. Furthermore, in chiral nematic cells with weak perpendicular boundary conditions and d/pâââ1, skyrmions could be formed spontaneously upon quenching samples from isotropic to the LC phase and exhibited similar dynamic behavior. LC cells were constructed using glass substrates with transparent indium tin oxide conductive layers and spin-coated with polyimide coatings (SE-1211 purchased from Nissan Chemical) to impose the finite-strength perpendicular surface boundary conditions for the LC director. Spin coating was done at 2700âr.p.m. for 30âs. The substrates were then baked for 5âmin at 90â°C and for 1âh at 190â°C to induce cross-linking of the alignment layer. The substrates were glued together, with the treated surfaces facing inward, and the cell gap was set with glass fiber segments dispersed in the ultraviolet-curable glue. The glue was cured for 60âs with an OmniCure UV lamp, Series 2000. Commercially available homeotropic cells (purchased from Instec) were also used. Electrical connections for voltage application across the depth of the LC cell were achieved by soldering leads to the ITOâelectrode surfaces. Finally, the LC was heated to the isotropic phase, infiltrated into the constructed cells via capillary action and sealed with 5-min fast-setting epoxy.
Generation, manipulation, pinning, and control of skyrmions
Skyrmions can be formed spontaneously upon thermally quenching the sample from the isotropic to the chiral nematic phase, by relaxing the LC from electrohydrodynamic instability (Fig. 9), or by a direct one-by-one optical generation using holographic laser tweezers23,32. Spontaneous formation of large densities of skyrmions was achieved by inducing an electrohydrodynamic instability in cells doped with CTAB by applying low-frequency voltage of Uâ=â5â25âV at frequencies within 2â10âHz. The control of this voltage and carrier frequency allows for the selection of the initial skyrmion packing fraction within 0.01â0.40 (Fig. 9). Laser-induced generation of individual skyrmions in all other cells without CTAB was done using optical tweezers comprised of a 1064ânm Ytterbium-doped fiber laser (YLR-10-1064, IPG Photonics) and a phase-only spatial light modulator (P512-1064, Boulder Nonlinear Systems). Using this setup, we can controllably produce arbitrary, dynamically evolving 3D patterns of laser light intensity within the sample and generate twisted structures by means of optically induced local reorientation of the director field known as optical Fredericks transition. Upon focusing this laser beam of power >50âmW in the midplane of the cell, the local LC director realigns away from the far-field background by coupling to the optical frequency electric field of the laser beam. Skyrmions were individually generated by laser tweezers at ~50âmW power and selectively pinned to the substrate surface to act as obstacles in desired locations using powers of 70â150âmW. The LCâs tendency to twist makes skyrmions energetically favorable at given confinement conditions, so that they form spontaneously after the uniform background of homeotropic cells is distorted by electrohydrodynamic instabilities or laser-induced realignment (after voltage or laser light are turned off, as shown in Fig. 9), whereas control of distortions in these two cases allows for defining initial densities and locations of skyrmions. In order to morph skyrmions and power their motions via macroscopically supplied energy, electric field was applied across the cell using a homemade MATLAB-based voltage-application program coupled with a data-acquisition board (NIDAQ-6363, National Instruments)24.
The means of controlling dynamics of skyrmions include voltage driving schemes, selection of LCs with different material parameters and chiral additives (Table 1), design of LC cell geometry, and strength of surface boundary conditions, and so on. Collective motion effects could be obtained both when simply using low-frequency oscillating field (yielding field oscillations at timescales comparable to the LCâs response time) and when modulating high-frequency carrier signals (e.g., at 1âkHz) at modulation frequencies that again yield modulation periods comparable to the LCâs response time. The large range of possibilities exists in terms of controlling the skyrmion schooling by varying electric signal waveforms, modulation, and carrier frequencies, and so on, but detailed exploration of all these possibilities is outside the scope of our present study.
Numerical modeling
We computer-simulated structures of skyrmions and their self-assemblies at experimental conditions (Figs. 1 and 5) by using LC free energy with elastic and electric coupling terms24,40,41:
where Frank elastic constants K11, K22, K33, and K24 represent the elastic costs for splay, twist, bend, and saddle-splay deformations of n(r), respectively. The chiral wavenumber of the ground-state chiral nematic mixture is defined as q0â=â2Ï/p and Îε is the dielectric anisotropy. We take K24â=âK22, as in previous studies32, whereas all other material parameters used correspond to the experimental values (Table 1). As the applied voltage is modulated, the competing electric and elastic torques are balanced by a viscous torque associated with rotational viscosity, γ, that opposes the fast rotation of the director1. The resulting director dynamics is governed by a torque balance equation1, γâni/âtâ=ââδW/δni, from which both the equilibrium n(r) and the effective temporal evolution of the director field towards equilibrium are obtained for the director, ni(t), where ni is the component of n along the ith axis (iâ=âx, y, z). As in experiments, we start our computer simulations from skyrmions embedded in a homeotropic LC background, for which the equilibrium director structure is obtained by minimizing free energy in Eq. (1) at no external fields (Fig. 1g), as detailed in our previous studies.24,29 Then, we minimize free energy to obtain skyrmionâs field configurationâs at various applied voltages (Fig. 1j, k) when starting from the skyrmion structure at no fields (Fig. 1g) within the computational volume as the initial condition. In a similar way, to obtain multi-skyrmion clusters or chains, we start from minimizing free energy at Uâ=â0 for seven (Fig. 5f) or five (Fig. 5j) axisymmetric skyrmions embedded in the homeotropic background of the computational volume and then use these structures as initial conditions to obtain clusters in the corresponding applied fields. Similar to experiments, skyrmions self-organize into clusters and chains (Fig. 5fâm) at corresponding voltages and translate laterally by about p/2 each time as we effectively turn voltage on and off and minimize free energy at the corresponding conditions. Periodic turning voltage on and off that corresponds to TU results in a periodic non-reciprocal director-field evolution that yields an asymmetric shifting of the skyrmions between the voltage-on and voltage-off states, resulting in a displacement within the computational volume24, similar to that seen in experiments. Such periodic displacements add to yield lateral translations of both individual skyrmions and their clusters and chains (Fig. 5fâm). Similar to experiments, the velocity vectors that we obtain from analyzing displacements of skyrmionic n(r) structures are antiparallel to the preimage vectors pi. Dynamic evolution of n(r) has been used to derive the intermediate states between the voltage-on and voltage-off states by taking snapshots of the director field during the process of evolution towards equilibrium within each voltage modulation period TU24. Once the director structures are obtained, we utilize a Jones matrix method23,24 to generate polarizing optical micrographs for experimental parameters such as sample thickness, optical refractive index anisotropy, and p (Table 1).
Optical microscopy, video characterization, and data analysis
Images and videos were obtained using charge-coupled device cameras Grasshopper (purchased from Point Grey Research, Inc.) or SPOT 14.2 Color Mosaic (purchased from Diagnostic Instruments, Inc.), which were mounted on an upright BX-51 Olympus microscope. Dry Ã2, Ã4, Ã10, and Ã20 objectives (with numerical apertures ranging from 0.3 to 0.9) were used, with different relative orientations of polarizers adjusted to increase contrast between the skyrmionic structures and the background. This contrast was integral to successfully analyze the polarizing micrographs for skyrmion motion using the open-source ImageJ/FIJI software (obtained from the National Institute of Health). Built-in particle-tracking tools were applied, through which skyrmion positional information and skyrmion number density were extracted for each frame. Then, the data analysis and plotting in MATLAB software (obtained from MathWorks) were performed to characterize trajectory pathways, velocity and polar order parameters, giant-number fluctuation scaling, and density fluctuations. The temporal evolution of both polar and velocity order was characterized by analyzing the positional data for skyrmion motion between frames of the videos. The velocity vector for an individual skyrmion vi was defined by drawing a vector between the skyrmionâs positions in consecutive frames of the video, pointing along the direction of motion4. The polar preimage vector for each skyrmion pi was defined by drawing a normalized vector between the south-pole and north-pole preimages. The giant-number fluctuations and scaling trends were analyzed using the skyrmion number density data obtained by means of the ImageJ/FIJI particle-counting features15. This was done by doing skyrmion number density analysis with time for 15 areas of different sizes, ranging from 25âµmâÃâ25âµm to 1250âµmâÃâ1250âµm, for each experimental video. Figure 4f represents the number density fluctuation analysis for one representative 400âµmâÃâ400âµm region as an example. Various regions within the samples were probed for each experimental video and a composite of 3â5 videos were analyzed for each case (individual, clusters, and chains of skyrmions), resulting in ~120 data points each to represent the different sample areas. The time period over which the fluctuations were characterized for each video was within 150â180âs. The density data points by area were compiled and plotted as logâlog plots of the mean particles <N> and root mean square ÎNâ=â<â(Nâââ<N>â)2â>â1/2 in Figs. 4e and 5e.
Data availability
All datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.
Code availability
MATLAB codes generated and analyzed during the current study are available from the corresponding author on request.
Change history
29 January 2020
An amendment to this paper has been published and can be accessed via a link at the top of the paper.
References
Chaikin, P. M. & Lubensky, T. C. Principles of Condensed Matter Physics (Cambridge University Press, Cambridge, 1995).
Ramaswamy, S. The mechanics and statistics of active matter. Annu. Rev. Condens. Matter Phys. 1, 323â345 (2010).
Marchetti, M. C. et al. Hydrodynamics of soft active matter. Rev. Mod. Phys. 85, 1143â1189 (2012).
Vicsek, T., Czirók, A., Ben-Jacob, E., Cohen, I. & Shochet, O. Novel type of phase transition in a system of self-driven particles. Phys. Rev. Lett. 75, 1226 (1995).
Narayan, V., Ramaswamy, S. & Menon, N. Long-lived giant number fluctuations in a swarming granular nematic. Science 317, 105â108 (2007).
Bricard, A., Caussin, J.-B., Desreumaux, N., Dauchot, O. & Bartolo, D. Emergence of macroscopic directed motion in populations of motile colloids. Nature 503, 95â98 (2013).
Shelley, M. Applied physics: on a roll. Nature 503, 43 (2013).
Needleman, D. & Dogic, Z. Active matter at the interface between materials science and cell biology. Nat. Rev. Mater. 2, 17048 (2017).
Sanchez, T., Chen, D. T. N., DeCamp, S. J., Heymann, M. & Dogic, Z. Spontaneous motion in hierarchically assembled active matter. Nature 491, 431 (2012).
Peruani, F., Ginelli, F., Bär, M. & Chaté, H. Polar vs. apolar alignment in systems of polar self-propelled particles. J. Phys. Conf. Ser. 297, 012014 (2011).
Yan, J. et al. Reconfiguring active particles by electrostatic imbalance. Nat. Mater. 15, 1095 (2016).
Ristroph, L. & Zhang, J. Anomalous hydrodynamic drafting of interacting flapping flags. Phys. Rev. Lett. 101, 194502 (2008).
Ellis, P. W. et al. Curvature-induced defect unbinding and dynamics in active nematic toroids. Nat. Phys. 14, 85 (2018).
Zhang, H. P., Beâer, A., Florin, E.-L. & Swinney, H. L. Collective motion and density fluctuations in bacterial colonies. Proc. Natl Acad. Sci. USA 107, 13626â13630 (2010).
Toner, J. & Tu, Y. Long-range order in a two-dimensional dynamical XY model: how birds fly together. Phys. Rev. Lett. 75, 4326 (1995).
Kudrolli, A., Lumay, G., Volfson, D. & Tsimring, L. S. Swarming and swirling in self-propelled polar granular rods. Phys. Rev. Lett. 100, 058001 (2008).
Peruani, F. et al. Collective motion and nonequilibrium cluster formation in colonies of gliding bacteria. Phys. Rev. Lett. 108, 098102 (2012).
Schaller, V., Weber, C., Semmrich, C., Frey, E. & Bausch, A. Polar patterns of driven filaments. Nature 467, 73â77 (2010).
Kumar, N., Soni, H., Ramaswamy, S. & Sood, A. K. Flocking at a distance in active granular matter. Nat. Commun. 5, 4688 (2014).
Palacci, J., Sacanna, S., Steinberg, A. P., Pine, D. J. & Chaikin, P. M. Living crystals of light-activated colloidal surfers. Science 339, 936â940 (2013).
Kernbach, S. Handbook of Collective RoboticsâFundamentals and Challenges (Pan Stanford Publishing, Boca Raton, 2013).
Werfel, J., Petersen, K. & Nagpal, R. Designing collective behavior in a termite-inspired robot construction team. Science 343, 754â758 (2014).
Ackerman, P. J., Trivedi, R. P., Senyuk, B., van de Lagemaat, J. & Smalyukh, I. I. Two-dimensional skyrmions and other solitonic structures in confinement-frustrated chiral nematics. Phys. Rev. E 90, 012505 (2014).
Ackerman, P. J., Boyle, T. & Smalyukh, I. I. Squirming motion of baby skyrmions in nematic fluids. Nat. Commun. 8, 673 (2017).
Sohn, H. R. O. et al. Dynamics of topological solitons, knotted streamlines, and transport of cargo in liquid crystals. Phys. Rev. E 97, 052701 (2018).
Skyrme, T. H. R. A non-linear field theory. Proc. R. Soc. Lond. Ser. A 260, 127â138 (1961).
Witten, E. Global aspects of current algebra. Nucl. Phys. B 223, 422â432 (1983).
Gennes, de, P. G. & Prost, J. The Physics of Liquid Crystals 2nd edn (Clarendon Press,: Oxford, 1993).
Foster, D. et al. Two-dimensional skyrmion bags in liquid crystals and ferromagnets. Nat. Phys. 15, 655â659 (2019).
Poulin, P., Stark, H., Lubensky, T. C. & Weitz, D. A. Novel colloidal interactions in anisotropic fluids. Science 275, 1770â1773 (1997).
Smalyukh, I. I. Liquid crystal colloids. Annu. Rev. Condens. Matter Phys. 9, 207â226 (2018).
Ackerman, P. J., van de Lagemaat, J. & Smalyukh, I. I. Self-assembly and electrostriction of arrays and chains of hopfion particles in chiral liquid crystals. Nat. Commun. 6, 6012 (2015).
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).
Shankar, S., Ramaswamy, S., Marchetti, M. C. & Bowick, M. J. Defect unbinding in active nematics. Phys. Rev. Lett. 121, 108002 (2018).
Saw, T. B. et al. Topological defects in epithelia govern cell death and extrusion. Nature 544, 212 (2017).
Kawaguchi, K., Kageyama, R. & Sano, M. Topological defects control collective dynamics in neural progenitor cell cultures. Nature 545, 327 (2017).
Morin, A., Desreumaux, N., Caussin, J. B. & Bartolo, D. Distortion and destruction of colloidal flocks in disordered environments. Nat. Phys. 13, 63 (2017).
Li, Q. et al. Reversible photoswitchable axially chiral dopants with high helical twisting power. J. Am. Chem. Soc. 129, 12908â12909 (2007).
White, T. J. et al. Phototunable azobenzene cholesteric liquid crystals with 2000 nm range. Adv. Funct. Mater. 19, 3484â3488 (2009).
Smalyukh, I. I., Lansac, Y., Clark, N. A. & Trivedi, R. P. Three-dimensional structure and multistable optical switching of triple-twisted particle-like excitations in anisotropic fluids. Nat. Mater. 9, 139 (2010).
Ackerman, P. J. & Smalyukh, I. I. Diversity of knot solitons in liquid crystals manifested by linking of preimages in torons and hopfions. Phys. Rev. X 7, 011006 (2017).
Acknowledgements
This research was supported by the National Science Foundation through grants DMR-1810513 (research), DGE-1144083 (Graduate Research Fellowship to H.R.O.S.) and ACI-1532235 and ACI-1532236 (RMACC Summit supercomputer used for the numerical modeling). We thank P. Ackerman, M. Bowick, M. Cates, A. Hess, T. Lubensky, D. Marenduzzo, S. Ramaswamy, M. Ravnik, M. Tasinkevych, J. Toner, J. Yeomans, and S. Zumer for discussions and P. Ackerman for technical assistance. We thank Corbin Sohn for the assistance with taking the photograph of a school of fish shown in Fig. 1a, which is part of the Sohn family collection.
Author information
Authors and Affiliations
Contributions
H.R.O.S. performed experiments (with assistance from C.D.L.) and numerical modeling. H.R.O.S., C.D.L., and I.I.S. analyzed data. I.I.S. wrote the manuscript (with the input from all authors), conceived the project, designed experiments, and provided funding.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work.
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Sohn, H.R.O., Liu, C.D. & Smalyukh, I.I. Schools of skyrmions with electrically tunable elastic interactions. Nat Commun 10, 4744 (2019). https://doi.org/10.1038/s41467-019-12723-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-019-12723-3