ABSTRACT
Topological solitons are quasi-particles of field deformations that are topologically nontrivial. Conceptualised by celebrated mathematicians and scientists including Carl Friedrich Gauss, Lord Kelvin, and Tony Skyrme, the ideas of knots, topology, and solitons have found a wealth of experimental embodiments in chiral liquid crystals in recent years. This is largely due to the exceptional capacity of the material for stabilising, characterising, and controlling topological solitons. In this review, I report the recent progress in realising and characterising topological solitons in chiral liquid crystals. Following the homotopy classification, I begin with the 1D twist walls and subsequently introduce 2D skyrmions and skyrmion bags, followed by 3D Hopf solitons including both hopfions and heliknotons. Physical properties of topological solitons are discussed, which include stability, soliton–soliton interaction, field-controlled transformation, field-activated dynamics, and crystalline self-assembly. Finally, the review ends by illustrating how liquid-crystal solitons can provide insights into topological solitons in other physical systems, and perspectives of future studies on the subject. The purpose of this review is to provide a fundamental overview of topological solitons in chiral liquid crystals.
Introduction
Topology is the study of properties that are preserved under continuous deformations [Citation1]. Topological properties include genus, knot invariant, linking number, etc. (). Fields, such as continuous vector fields, can also be topological. This can be illustrated by the fact that certain field configurations cannot be deformed continuously into one another without introducing singularities that tear the continuum, making these fields topologically distinct (). Topological solitons are topologically nontrivial particle-like field configurations embedded in a far-field background [Citation5,Citation6].
The concepts of solitons and knots have their roots in the early nineteenth century when Carl Friedrich Gauss envisaged that knots in fields could behave like particles and Lord Kelvin proposed a theory of chemical elements as different knots in Aether [Citation7]. Tony Skyrme, in the second half of the twentieth century, developed a model that describes nucleons as topological solitons in a nonlinear field theory [Citation8]. The model was later demonstrated to be a low-energy effective model of quantum chromodynamics [Citation9]. Today, topological solitons have been discovered or realised in many branches of physics, including nuclear physics, condensed matters, optics, and cosmology [Citation10–19]. Not only are topological solitons of great fundamental interest, but they also promise technological applications because of their unique particle-like properties and the capability of encoding information in topology [Citation20].
Liquid crystals (LCs) are ordered fluids with long-range orientational order and partial or no positional order [Citation21,Citation22]. The long-range orientational order allows the director field to be defined in nematic LCs, which is a non-polar line field of the average molecular orientation (). Beyond the multi-billion-dollar industry of electro-optic devices including displays [Citation23], LCs have served as a platform for studying topological structures and defects [Citation18,Citation24,Citation25]. LCs are excellent for studying topological solitons due to their unparalleled accessibility in generating and manipulating topological solitons, as well as structurally analysing topological solitons they host [Citation26]. Moreover, solitons in chiral LCs have been shown to gain additional stability. In chiral LCs, has an intrinsic tendency to twist due to molecular chirality; the distance over which twists by at equilibrium is the helical pitch (). The extra chirality term in the energy functional allows the system to evade the prediction of the Derrick-Hobart theorem, which shows solitons cannot be stabilised in dimensions of three and beyond within a linear field theory [Citation27–29].
This review is organised as follows. Homotopy groups and LC free-energy models are first introduced for classifying the topology and modelling the field configurations of topological solitons. Following the classification of homotopy groups, I provide a brief and non-comprehensive review of the progress in topological solitons in chiral LCs. I begin with the simplest one-dimensional (1D) twist walls, and progress to 2D skyrmions. The highlights include the stability of skyrmions, different geometric configurations, field-induced dynamics, and novel 2D field configurations with arbitrary skyrmion number – skyrmion bags. Next, I introduce 3D Hopf solitons, also known as hopfions, with the focus on their field-induced topological transformation. This is followed by heliknotons, a new class of 3D topological solitons discovered in the helical background of chiral LCs. I detail the fascinating properties of heliknotons, including 3D localisation, stable knotted vortex lines, self-assembled crystals, electrostriction, etc. I also delve into the geometric connection between hopfions and heliknotons and how they can reversibly transform into one another. Additionally, the 3D hopping dynamics of Hopf solitons during repeated inter-transformation between the two topologically identical, but geometrically distinct configurations is discussed. Finally, this review concludes by discussing how insights gained from studying LC topological solitons can be applied to studying solitons in other systems.
Homotopy classification of topological solitons
How do we describe or classify the topology of a field? A field is a map from the physical configuration space to the target space , , where and While the configuration space can be intuitively understood as the space the field lives in, the target space (or the order-parameter space in physical systems) depends on the nature of the field. For example, when the field is a unit vector confined to 2D orientations, the target space is the space composed of all possible orientations of the 2D unit vector, which is a unit circle (). In topology, if a field can be continuously deformed into , and are homotopic and belong to the same homotopy class. Moreover, the set of homotopy classes can form a homotopy group. Homotopy groups are often defined in terms of spheres; the nth homotopy group of , or , summarises topologically distinct maps from the n-dimensional sphere to the m-dimensional sphere and can be used to classify topological solitons () [Citation1,Citation5]. As an example, a 2D unit vector field in the 1D Euclidean space is a map . In this case, can be compactified into by identifying the positive and negative infinity in as one point without introducing discontinuity if , thus allowing the implementation of homotopy groups for topological classification (). The integer number of rotations that goes through from to is invariant under smooth geometric perturbations and is referred to as the topological degree or charge of the field. Thus, identifies with , the group of integers under addition. In nematic LCs, the nonpolar nature of LC director leads to the order-parameter space being , a sphere with diametrically opposite points identified (). A major consequence is the existence of LC defect lines in 3D, since while is a trivial group. However, as most homotopy groups of are homeomoric to homotopy groups of , and a continuous director field in a simply connected manifold can always be ‘vectorised’ into a vector field [Citation30], a director field can be replaced with a vector field for simplicity in many cases.
Free energy of chiral LCs and chiral LC ferromagnets
An important tool for investigating the structure and stability of topological solitons is numerical modelling of field configurations by free-energy minimisation. In chiral LCs, the free energy associated with deformations in the director field (treated as a unit vector field here) can be described by the Frank–Oseen free energy functional [Citation3,Citation21]
where is the helical pitch and , , and are the Frank elastic constants describing the energetic costs of splay, twist, bend, and saddle-splay deformations, respectively. When the LC is confined and surface anchoring is present, Equation (1) is supplemented by a surface anchoring potential at the confining substrates,
where is the surface anchoring strength and is the easy-axis orientation. To model the effect of electric field on , the dielectric coupling term can be included as
where is vacuum permittivity and is the dielectric anisotropy of the LC. In a chiral LC ferromagnet where magnetically monodomain nanoplates are uniformly dispersed in the hosting chiral LC, responds to weak magnetic fields due to the coupling between the magnetic moment of the nanoplates and the magnetic field, and the coupling between and magnetic nanoplates [Citation31]. In the case of strong homeotropic (perpendicular) surface anchoring of LC molecules to the magnetic nanoplates, the local magnetisation is assumed to be collinear with and vectorises the director field (). The magnetisation unit vector field can thus replace in Equations (1)–(3), and an additional linear magnetic coupling term between and an applied magnetic field reads
where is vacuum permeability and is the magnetisation of the material (approximated to be uniform) as the product of the nanoplate density and the average magnetic moment of a nanoplate . The material parameters for a few commercially available LC mixtures are summarised in .
It is important to note that the vectoral approach in Frank–Oseen free energy does not account for the non-polar nature of LC director field or the reduction in scalar order parameter near defects. In these cases, the tensorial approach and the Landau-de Gennes expansion is advised for a more complete description of free-energy cost in the perturbation of order in LCs [Citation22,Citation32]. Nevertheless, Frank–Oseen free energy is adequate for modelling orientable field configurations, especially when the detailed structures around defects are not the primary focus.
1D twist walls
1D solitonic walls represent the simplest of topological solitons. In condensed matter systems such as LCs and magnets, solitonic walls are rather common [Citation21,Citation33–35]. In chiral materials, these solitonic walls are also referred to as twist walls and can be stabilised by the frustration between the material’s chirality, surface energy, and coupling to external fields. In material systems where the order-parameter is a unit vector, an integer number of full-revolution (2π) twists of the order-parameter field is embedded in the uniform far-field, characterised by the homotopy group (). This is different for systems with a nonpolar order-parameter field, such as LCs. In a solitonic LC twist wall, the nonpolar director is confined in a 2D plane () and an integer number of half-revolution (π) twists can be embedded in the non-polar far-field (). The integer number of π-twists identifies with the topological charge. The order-parameter space is , which is homeomorphic to a circle, or the 1-sphere (). Therefore, nonpolar LC twist walls are also classified by the first homotopy group of , .
In LC samples where boundary conditions are often imposed at confining substrates, the configuration of twist walls depends on the effective anchoring strength. With perpendicular boundary conditions imposed by strong homeotropic anchoring, LC twist walls exist in the form of cholesteric fingers of the 3rd type, or CF-3 [Citation36]. In a CF-3, The director orientations are confined approximately in the 2D plane perpendicular to the direction that the finger extends (, top), and a 1D topological soliton in terminates at two half-integer defects near the substrates (). The defects extend along the translationally invariant direction and form disclination lines. We note that, without the approximation that director orientations are confined in 2D, both twist walls and disclination lines are described by . Notably, the director field around the two twist disclinations is non-orientable and the tangent of the director forms Möbius strips [Citation37]. Therefore, numerical modelling of configurations including twist disclinations requires the tensorial approach in Landau-de Gennes free energy.
Under boundary conditions with weak anchoring energy, the defect lines in a CF-3 escape beyond the confining surfaces, rendering a twist wall configuration translationally invariant across the thickness (). The structural stability between CF-3, twist wall, and the uniform state depends on the sample thickness-to-pitch ratio () and the effective anchoring strength , where is the surface anchoring energy density and is the average elastic constant in Landau-de Gennes free energy (). In brief, CF-3s are stable at large thicknesses and anchoring strengths, while in the weak anchoring regime, twist walls become the structurally stable configuration. The numerically derived stability between CF-3 and twist wall is consistent with experimental observations (), where CF-3s are observed in common LC nematics () with µm and Jm−2, corresponding to strong effective anchoring, while 1D twist walls are ubiquitously observed when µm [Citation3,Citation12,Citation33,Citation36,Citation37]. The cross-sectional 3PEF-PM images of a twist wall and a CF-3 reveal their detailed structural difference across the thickness (). The absence of high-energy defect lines in twist walls also allow them to be more abundant within the sample, lowering the overall elastic energy through twisting, while a CF-3 has a tendency to shrink due to the defect-induced tension. Overall, these results show the stability between 1D twist wall, CF-3, and the uniform state can be controlled by a set of material and sample parameters.
2D skyrmions
2D skyrmions (sometimes referred to as baby skyrmions due to their reduced dimensionality compared to the ones originally proposed by Tony Skyrme [Citation8]) are the topological solitons labelled by elements of the second homotopy group, for systems described by unit vector fields and for systems described by nonpolar director fields. Since and are homeomorphic, unit vector fields may be used in place of the nonpolar LC director field for convenience. Conceptually, a 2D skyrmion in chiral materials can be considered as configured by looping a twist wall into a circle [Citation37]. The topological charge of a 2D skyrmion, also called the skyrmion number, counts the integer number of times the entire field wraps around its target space , and is defined by . shows an elementary skyrmion that covers exactly once. Observed and studied extensively in LCs with strong perpendicular boundary conditions, a toron is a 2D skyrmion dressed by two point defects at its ends, with the existence of the defects mandated by the transformation of topology from a skyrmion into the uniform trivial state near the confining substrates (). Torons can therefore be understood as fragments of skyrmions terminating on singular point defects characterised by the same homotopy group, Notably, the two oppositely charged point defects are also a monopole and anti-monopole pair in the emergent magnetic field [Citation38]. Torons are common in both LCs and magnets, where the termination of skyrmions at point defects is related to confinement and other stabilising mechanisms [Citation3,Citation12,Citation39–43]
While torons are commonly observed as particle-like solitons in confined LCs with a sample thickness comparable to the helical pitch () and with strong perpendicular conditions [Citation39,Citation44], skyrmions with translationally invariant topology emerge when the strength of the surface anchoring or the size of the overall structure is reduced () [Citation12,Citation33]. The stability of skyrmions and torons, and their potential inter-transformations have been systematically studied as depending on both the thickness-to-pitch ratio () and the dimensionless effective anchoring strength () [Citation3]. Torons are stable compared to skyrmions or topologically trivial states when the surface anchoring is sufficiently strong and the thickness is sufficiently thick. At a smaller thickness, torons’ stability is lost to the uniform state. Skyrmions, on the other hand, are stable at weak anchoring, and the transition between torons and skyrmions happens at . Depending on the anchoring strength, the transformation between torons and skyrmions can occur at cell thicknesses of 10s of µm or sub-µm. Out of the non-exhaustive competing states, skyrmion or toron lattices can be the lowest energy state within certain parameter regions (), pointing at potential topological phases of matter analogous to the skyrmion lattice phase is solid-state chiral magnets [Citation10]. In practice, torons and skyrmions can inter-transform and co-exist at µm and µm, under strong homeotropic anchoring in a chiral LC based on the nematic 5CB. The structural difference between torons and skyrmions can be revealed by the distinct 2D patterns in the polarizing optical microscopy images. Nonlinear fluorescence images further reveal the detailed difference in their 3D structures (). Notably, the two point defects near the substrates are present in a toron but not in a skyrmion.
Helical-background skyrmions
Skyrmions can also be embedded in a helical background, contrasting the more commonly known axisymmetric, whirling configuration in the uniform background. The configuration of a helical-background skyrmion can be considered as splitting and splicing the uniform-background skyrmion while extending the helical region into a uniformly helical background with constant helical axis , which defines the axis that twists around in the far field (). Helical-background skyrmions can be stabilised in LC cells with planar boundary condition by an electric field along . The stabilised skyrmion extends in the in-plane direction and terminates at point defects, or monopoles in the emergent magnetic field [Citation38]. The energetics of in-plane helical-background skyrmions can be controlled by electric fields through the quadratic dielectric coupling. Additionally, when the LC is doped with ferromagnetic nanoplates with strong perpendicular boundary condition and forms a monodomain LC ferromagnet, the magnetisation field and the collinear director field additionally couple linearly to the applied magnetic field. This allows the length of the skyrmion and the separation velocity of the defect (monopole-antimonopole) pair to be controlled by the external fields (). The controlled dynamics of skyrmions and monopoles can lead to applications utilising the regions around monopoles with high-degree distortions in as energetic traps for cargo transport of micro- and nanoparticles.
Skyrmion bags
Beyond elementary skyrmions with unity skyrmion number, high-charge topological field configurations with arbitrary skyrmion number have recently been realised in LCs and modelled in solid-state magnets [Citation4,Citation45]. Dubbed ‘skyrmion bags’, these stable composite structures consist of antiskyrmions contained within a larger skyrmion, with the latter acting like a confining bag (). Skyrmion bags demonstrate a similar pair-wise interaction to their elementary counterparts and can be used to encode information in their topological degree of freedom. Conventional proposals of racetrack memory involve encoding bits of information by separations between elementary skyrmions, which allows for limited density and is susceptible to data loss due to their repulsive interaction [Citation20]. On the other hand, skyrmion bags utilise topology-encoded information to overcome this problem. For example, a train of skyrmion bags, S(N1), S(N2), … , encode information by the integer Nj, the number of (anti)skyrmions in each bag, enabling a new paradigm-changing approach in data storage and processing devices [Citation46].
3D hopfions
In 1931, German mathematician Heinz Hopf discovered the Hopf fibration (also called Hopf map or Hopf bundle) and demonstrated how space can be composed of fibres of circles, with each circle mapping to a single point in [Citation47]. This shows that a continuous map from a higher-dimensional sphere to a lower-dimensional sphere can be nontrivial, an unexpected result to the mathematical community at the time. Indeed, the third homotopy group of is identified with the group of integers, i.e. (). 3D Hopf solitons, or hopfions, are physical realisations of Hopf fibration, with each closed-loop in the 3D configuration space corresponding to a single distinct orientation in the order-parameter space for systems described by unit vectors [Citation5,Citation6,Citation48–53] (). Here the 3D configuration space can be given that it can be compactified to without introducing discontinuities. An example is the one-point compactification by identify the infinity of as one point when the far-field is uniform. Realisation of stable 3D hopfions had been challenging, as the Derrick-Hobart theorem shows the absence of localised stable solutions in linear filed theories in dimensions three or higher [Citation27,Citation28]. After long-preceded theoretical interest, stable hopfions were recently demonstrated in LCs [Citation54], LC ferromagnets [Citation13], and solid-state magnets [Citation55]. In LCs where the director field is nonpolar and the order-parameter space is , is homeomorphic to , and unit vector fields may be used in place of the non-polar LC director field.
Static and stable hopfions were first experimentally realised in chiral LCs and chiral LC ferromagnets where chirality and confinement help evade the prediction of Derrick-Hobart theorem and contribute to stability [Citation13,Citation54]. Topologically equivalent to the Hopf fibration, in a hopfion, the preimage (inverse image) of each point in the order-parameter space forms a loop in the 3D configuration space, and the preimages of all points in fill the 3D space (). Here preimage indicates the inverse image of the map , and can be obtained by finding the regions in that correspond to a given unit vector orientation in (; note that is used interchangeably with in ). The associated topological charge, Hopf index, can be geometrically obtained by the linking number of a pair of preimages of any two points in the order-parameter space [Citation56]. For elementary Hopf solitons with unity charge, the preimages are linked exactly once. A numerical integration method can also be applied to obtain the Hopf index [Citation57,Citation58]. Hopfions of Hopf indices and spontaneously emerge in samples of chiral LCs and chiral LC ferromagnets where the confining substrates impose perpendicular boundary conditions and the thickness over pitch ratio , while they can also be controllably generated by laser tweezers, as detailed in Refs [Citation13]. and [Citation54] (). The configuration and topology of LC hopfions were unambiguously identified by 3PEF-PM imaging, where experimentally reconstructed and preimages are consistent with numerical simulations based on minimisation of free energy (). LC hopfions show hard-particle-like elastic interactions and form 2D hexagonal lattice in the presence of a lateral confinement ().
Topological transformation of hopfions
Topological solitons hosted in chiral LC ferromagnets with equally strong facile responses to electric and magnetic fields provide an ideal platform for studying the interplay between external fields and topological solitons [Citation31]. In LC display industry, switching of pixels involves topologically trivial states [Citation23]. Understanding switching of topological configurations could additionally lead to multistable energy-efficient electro-optic and data devices. To prepare for topological transformation of solitons controlled by external fields, two elementary doughnut-shaped hopfions in a chiral LC with negative dielectric anisotropy () were coaxially arranged to create a composite hopfion with [Citation58] (). The preimage of each point in is two separate unlinked loops (Alexander-Briggs notation of knots [Citation5]), and preimages of two distinct points form a pair of Hopf links, giving rise to the total linking number of −2 (). Upon application of an intermediate electric field perpendicular to the confining substrates, uniform far-field regions that separate the knotted regions of two elementary hopfions disappear, and smoothly morphs into a single solitonic configuration. The partial merging of two elementary hopfion is demonstrated by the change in the topology of preimages, where a critical polar angle in separates two subspaces: (a) single-loop () preimages of points with and (b) preimages of points with in the form of two separate unlinked loops (). This gives rise to a rich linking topology between preimage pairs depending on their respective polar angle in (). Importantly, while the topology of individual preimage depends on and changes with , the overall field topology of is preserved, as evidenced by the overall linking number of the reconstructed preimages from 3PEF-PM images and numerical integration of Hopf index (). Above a critical , the interplay between the elastic and dielectric coupling energy causes to transform discontinuously from to through a series of short-lived singular transient states (). These processes demonstrate both complex morphing of preimage induced by electric field, while preserving the field topology, and the field-induced topological transformation of . Notably, similar topological transformations are also demonstrated for a hopfion consists of two oppositely charged elementary hopfions [Citation58]. Magnetic field also causes topological transformation of solitons hosted in chiral LC ferromagnets. Unlike the quadratic dielectric coupling to , the response of magnetisation field to the external magnetic field is polar, leading to rich transformation of both preimage and soliton topology. A structural stability diagram shows a wealth of relaxed knotted field configurations when a complex soliton is subject to a variable magnetic field . The equilibrium solitonic structures include those with and without singular defects, i.e. torons and hopfions, respectively, and with simple and complex topology of preimage ().
Heliknotons
Heliknotons are a new class of 3D knot solitons realised in the helical (cholesteric) background of chiral LCs [Citation14]. The LC helical field is ubiquitous in chiral materials (often the ground state) and can be considered as comprising a triad of orthonormal fields, namely the material director field and two immaterial line fields, the helical axis field and , which is the field perpendicular to both and (). Heliknotons were found to display field topology resembling Hopf fibration in and knotted half-integer singular vortex lines in and (). In conventional LCs such as 5CB, heliknotons spontaneously form after quenching from isotropic to chiral nematic phase when a moderate electric field was applied along the far-field helical axis , though they were later demonstrated to be stable without in materials with proper elastic anisotropy [Citation59]. Heliknotons can also be controllably generated using laser tweezers. Heliknotons are stable in the bulk of chiral LCs and display properties of 3D-localised particles including Brownian motion and orientation-correlated motion along (). This contrasts with 3D LC hopfions in a uniform far-field background, in which case confinement is required for stability. The pair-wise interaction between the heliknotons is anisotropic and highly tunable from 10 s to 1000 s of , depending on the LC material, and sample thickness (). The inter-heliknoton interaction arises from sharing long-range perturbations in the background helical field and minimising the overall free energy, much like the interaction between nematic colloids [Citation60,Citation61]. The interactions give rise to heliknoton crystals with different symmetries through self or guided assembly, including 2D rhombic latices, 2D kagome lattices, and 3D triclinic lattices (). Since LCs respond strongly to electric fields, LC heliknoton crystals exhibit giant anisotropic electrostriction, where two lattice parameters extend by 44% and 4% with less than 1.5 V change in voltage ().
The structure of heliknotons were revealed by consistent results of 3PEF-PM imaging and numerical simulations (); in heliknotons from interlinked closed-loop preimages with linking number + 1 (), and and exhibit singular half-integer vortex lines that close into a trefoil torus knot, labelled as T(2,3) or 31. This demonstrates a heliknoton as dual embodiments of a Hopf soliton and a vortex knot. The knot of singular vortex line in and also co-localise with a knotted tube of meron (fractional skyrmion) in [Citation14]. Unlike transient vortex lines that shrink and vanish due to energetically costly cores and distorted order around them, vortex-meron knots in heliknotons are energetically favourable and stable because of being nonsingular in the material field . The twisted structures with the same handedness of the hosting chiral material also aids in stability. Under the same condition that stabilises the elementary heliknoton, solitonic structures consistent with higher-degree heliknotons were also observed (), where they display Hopf indices of and , and form T(2,5) and T(2,7) torus knots (51 and 71). Whether heliknotons with even higher charges can be systematically generated and properties pertaining to their stability and many-body interactions require further study.
Heliknoton-hopfion transformation and 3D hopping
Heliknotons were later found to be stabilised in bulk chiral LCs without any electric field applied in nematic mixtures of rod-like and bent-core molecules that exhibit reduced bend elastic constant [Citation59,Citation62,Citation63]. The tunable elastic anisotropy allows the background far-field to be continuously varied between uniform, conical, and helical states with a continuous cone angle by an electric field when [Citation64]. In the uniform background, the hosted Hopf solitons are conventionally referred to as hopfions [Citation13,Citation54,Citation58], and in the helical background, Hopf solitons are identified as heliknotons [Citation14] (). Numerical simulations reveal the structural stability diagrams of Hopf solitons depending on elastic anisotropy, electric field, and confinement conditions, in topologically trivial backgrounds with varying (). In the LC bulk without any confinement (), at zero electric field, Hopf solitons are metastable particle-like, spatially localized structures in the helical background with higher energy density than the background when . At larger , heliknotons are metastable in the presence of an electric field, consistent with the original report that heliknotons require a stabilising in chiral LCs with that were used in the study [Citation14]. Hopf solitons can also be metastable in the bulk conical background. Notably, between metastable regions and across the helical-uniform and helical-conical boundaries for the background field, stable Hopf solitons with lower energy than the background fill the computational volume and form stable crystalline assemblies. Importantly, regardless of the distinct embedding backgrounds with different , the Hopf indices of the stabilised solitons remain unchanged. In the presence of confinement and perpendicular boundary conditions (), additional metastability regions of Hopf solitons emerge in the uniform background with (hopfions), differing from the case of no hopfions in the fully unwound LC bulk (). This further demonstrates how confinement and boundary conditions helps stabilise hopfions [Citation54,Citation65,Citation66].
The structural stability diagram with confinement suggests an in situ pathway for inter-transformation between a heliknoton and a hopfion by varying (), which was verified experimentally. The continuous inter-transformation provides further evidence that heliknotons and hopfions are indeed geometric embodiments of the same field topology (). The heliknoton-hopfion inter-transformation draws analogy with the transformation between geometric shapes of identical topological invariants, e.g. the surface of a coffee mug and a doughnut both have genus one (). Remarkably, the reversible full-cycle transformation is accompanied by a displacement of approximately half a pitch along the long axis of the heliknotons (), enabled by the nonreciprocal evolution of . Such transformation and the concomitant displacement can be repeated with a periodic voltage switching, leading to an activated propelling motion, analogous to the squirming motion of 2D LC skyrmions () [Citation67].
A close inspection of the translational motion of Hopf solitons under periodic inter-transformation revealed deviation from a linear trajectory (). Such deviation arises from the contrasting energetic landscapes between hopfions in a uniform background and heliknotons in a helical background. The former is monotonically increasing when the soliton displaces away from the energy minimum at the sample midplane, while the latter contains multiple energy minima along the far-field helical axis and the midplane is an unstable local energy maximum. This is evident from the different long-axis orientations of heliknotons at static equilibrium with no applied voltage (), since the orientation of the soliton correlates with its position along in a helical background [Citation14]. Therefore, in each inter-transformation cycle, the combined propelling motion, complex energetic landscape along (for heliknoton), and the correlated orientation and position of heliknoton in a helical background, results in a combined squirming motion in lateral directions and vertical motion along (or ) – an effective hopping motion in the 3D space. Additionally, a Hopf soliton also demonstrates a propelling motion when a hopfion responds to a low-magnitude voltage modulation that perturbs the background between the uniform and the helical state, while not completely inter-transforming between a hopfion and a heliknoton (). This demonstrates the universality and richness of nonreciprocal evolution in activating the propelling motion of solitons in LCs.
Topological solitons beyond LCs
LCs are an ideal material system for studying topological solitons due to their accessibility in the generation, manipulation, and optical characterisation of solitons. To date, a plethora of topological solitons in different dimensions have been realised and their equilibrium and out-of-equilibrium properties analysed, as partially covered in this review. Remarkably, because of the universality of the underlying topology, the insights into topological solitons may be shared between differential materials and physical systems such as liquid crystals, magnets, Bose-Einstein condensates, optics, etc [Citation10–19]. An exceptional example is solid-state chiral magnets and chiral LCs [Citation68], where the micromagnetic Hamiltonian of magnets takes a similar functional form to the Frank–Oseen free energy of LCs under the one-constant approximation by having . Notably, the molecular chirality in chiral LCs plays a role analogous to that of the Dzyaloshinskii–Moriya interaction in solid-state chiral magnets [Citation29,Citation68]. The close connection between the mean-field energy functional of the two systems has allowed LCs and magnets to mutually inform and enrich the understanding of topological solitons in each system. Examples of topological solitons in both LCs and magnets include 2D skyrmions [Citation3,Citation12,Citation69], 2D skyrmion bags [Citation4,Citation45], 3D Hopf solitons (hopfions and heliknotons) [Citation41,Citation66,Citation70,Citation71], as well as solitonic structures containing defects such as torons and chiral bobbers [Citation3,Citation39,Citation42,Citation54]. However, it is important to note that fundamental differences between the two systems exist; the nonlocal dipole–dipole interaction, magneto-crystalline anisotropy, and the intrinsically polar order parameter are unique to solid-state magnets. Interesting properties pertaining to topological solitons in magnets include the skyrmion Hall effect and topological Hall effect [Citation72], Hopf fibration in the emergent magnetic field [Citation70,Citation73], and topological phases of matter [Citation68], though detailing topological solitons in solid-state magnets is beyond the scope of this review.
Conclusion
Topological solitons are particle-like field deformations that are topologically nontrivial. In this review, I give a brief introduction to topological solitons in the nematic alignment (director) field of chiral LCs based on their classification using homotopy groups. The stability of topological solitons arises from their distinct topology that ensues energetic barriers between topologically distinct states – the so-called topological protection. Material properties such as chirality and elastic anisotropy, as well as confinement, surface anchoring, coupling to external fields, etc., also affect the energetic landscape of the system and the stability of solitons. The fascinating properties of LC topological solitons can largely be attributed to the fact that the contrast between the topological quasi-particles and the background is in the field alignment rather than the chemical composition. This enables topological solitons the capacity to be reconfigured and morphed, which determines soliton–soliton interaction, self-assembly, electrostriction, stimuli-induced dynamics, etc. This contrasts with ‘real’ particles such as colloids in a liquid, which are chemically distinct from their backgrounds.
The unparalleled accessibility of LCs in generating, manipulating, and characterising a rich collection of nontrivial field configurations has allowed LCs to be an ideal testbed for topological defects, topological solitons, and the underlying theories. The facile response of LCs to externally applied electric fields and LC ferromagnets to magnetic fields further allow additional knobs of control. Together with their topologically protected particle-like properties, LC topological solitons can find useful technological applications in electro-optics, data storage, data processing, microfluidics, etc.
Several interesting subjects pertaining to LC topological solitons were not covered in this review. For example, it was recently shown that nanoparticles are enriched in the distorted regions of heliknotons due to free-energy minimisation, making topological solitons a potential building block of templates for nanoparticle assembly and plasmonic metamaterials [Citation73]. Topological solitons in birefringent and reconfigurable media such as LCs can also be used to control light propagation or interact with optical solitons [Citation74,Citation75]. Moreover, a hierarchical self-assembly between 1D twist walls and defect lines result in a newly discovered family of solitonic structures, Möbiusons, which display Möbius-strip-like director field topology and self-propelled translational and rotational motions when activated with electric pulses [Citation37].
Moving forward, a diverse array of new frontiers of topological solitons in LCs and beyond await further investigations. Can we go beyond the 2D target space or associated with uniaxial polar and nonpolar order parameters? A 3D target space such as or its quotient groups will enable the realisation of topological solitons in the 3D configuration space, which share the same topology as nucleons in Skyrme’s nonlinear field theory [Citation8]. It would be interesting to see if such topological solitons can be stably realised in biaxial LC systems [Citation22,Citation76,Citation77]. Also, since topological solitons and colloids can both form crystals, it is tantalising to consider the possibility of hybrid colloidal-solitonic crystals in LC colloidal systems [Citation14,Citation61,Citation78]. Moreover, recent developments in ferroelectric LCs raise questions about how long-range dipolar interactions alter the stability of topological field configurations, and how hydrodynamic effects driven by ultralow electric field interact with solitons [Citation79]. Interestingly, in contrast to inanimate systems, biological systems made up of anisotropically shaped cells, such as tissues, embryos, and bacterial communities exhibit LC ordering and out-of-equilibrium activity from energy consumption [Citation80–83]. It remains largely unknown if topological solitons emerge in these active nematic systems, and whether nontrivial alignment field configurations play a role in their physiological functions. Finally, from a technological point of view, further studies and engineering of the stability, self-assembly, and controlled dynamics of topological solitons in different material systems will facilitate the advent of topological-soliton-enabled technologies. Overall, topological solitons in LCs and beyond represent a burgeoning field with interest in fundamental science and potential for practical applications.
Acknowledgement
I am grateful to Ivan Smalyukh for his invaluable advice, guidance, and support, during my years as a PhD student. I also thank the past and present Smalyukh Lab members for fruitful discussions and collaborations. I extend my gratitude to Ivan Smalyukh, Noel Clark, Mark Dennis, and Jing Yan for nominating me for the 2022 Glenn H. Brown Prize, and to the International Liquid Crystal Society for bestowing me the award. My PhD dissertation research was supported by US Department of Energy, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering, Award No. ER46921 and No. DE-SC0019293, and National Science Foundation grant DMR-1810513. I am currently a Damon Runyon Fellow supported by the Damon Runyon Cancer Research Foundation (Grant No. DRG-2446-21).
Disclosure statement
No potential conflict of interest was reported by the author(s).
References
- Hatcher A. Algebraic topology. Cambridge, England: Cambridge University Press; 2002.
- Crane K, Segerman H. Topology_joke.Jpg (3648×2886). Available from: http://www.segerman.org/images/topology_joke.jpg.
- Tai J-SB, Smalyukh II. Surface anchoring as a control parameter for stabilizing torons, skyrmions, twisted walls, fingers, and their hybrids in chiral nematics. Phys Rev E. 2020;101(4):042702. doi: 10.1103/PhysRevE.101.042702
- Foster D, Kind C, Ackerman PJ, et al. Two-dimensional skyrmion bags in liquid crystals and ferromagnets. Nat Phys. 2019;15(7):655. doi: 10.1038/s41567-019-0476-x
- Manton N, Sutcliffe P. Topological solitons. Cambridge, England: Cambridge University Press; 2004.
- Shnir YM. Topological and non-topological solitons in scalar field theories. Cambridge, England: Cambridge University Press; 2018.
- Turner JC, Griend PVD, Warner C. History and science of knots. Singapore: World Scientific; 1996.
- Skyrme THR, Particle states of a quantized meson field. Proc R Soc Lond Ser A Math Phys Sci. 1961; 262:237.
- Adkins GS, Nappi CR, Witten E. Static properties of nucleons in the skyrme model. Nucl Phys B. 1983;228(3):552. doi: 10.1016/0550-3213(83)90559-X
- Yu XZ, Onose Y, Kanazawa N, et al. Real-space observation of a two-dimensional skyrmion crystal. Nature. 2010;465(7300):901–904. doi: 10.1038/nature09124
- Mühlbauer S, Binz B, Jonietz F, et al. Skyrmion lattice in a chiral magnet. Science. 2009;323(5916):915. doi: 10.1126/science.1166767
- Ackerman PJ, Trivedi RP, Senyuk B, et al. Two-dimensional skyrmions and other solitonic structures in confinement-frustrated chiral nematics. Phys Rev E. 2014;90(1):012505. doi: 10.1103/PhysRevE.90.012505
- Ackerman PJ, Smalyukh II. Static three-dimensional topological solitons in fluid chiral ferromagnets and colloids. Nat Mater. 2017;16(4):426–432. doi: 10.1038/nmat4826
- Tai J-SB, Smalyukh II. Three-dimensional crystals of adaptive knots. Science. 2019;365(6460):1449. doi: 10.1126/science.aay1638
- Tsesses S, Ostrovsky E, Cohen K, et al. Optical skyrmion lattice in evanescent electromagnetic fields. Science. 2018;361(6406):993. doi: 10.1126/science.aau0227
- Zhang X, Zhou Y, Ezawa M. Antiferromagnetic skyrmion: stability, creation and manipulation. Sci Rep. 2016;6(1). doi: 10.1038/srep24795
- Lee W, Gheorghe AH, Tiurev K, et al. Synthetic electromagnetic knot in a three-dimensional skyrmion. Sci Adv. 2018;4(3):eaao3820. doi: 10.1126/sciadv.aao3820
- Chuang I, Durrer R, Turok N, et al. Cosmology in the laboratory: defect dynamics in liquid crystals. Science. 1991;251(4999):1336. doi: 10.1126/science.251.4999.1336
- Wu J-S, Smalyukh II. Hopfions, Heliknotons, Skyrmions, Torons and both abelian and nonabelian vortices in chiral liquid crystals. Liq Cryst Rev. 2022;1:1–35.
- Fert A, Cros V, Sampaio J. Skyrmions on the track. Nat Nanotechnol. 2013;8(3):152–156. doi: 10.1038/nnano.2013.29
- Chaikin PM, Lubensky TC. Principles of condensed matter physics. Cambridge, England: Cambridge University Press; 1995.
- de Gennes PG, Prost J. The physics of liquid crystals. Oxford: Clarendon Press; 1993.
- Cheng K-T. Liquid crystal display—present status and emerging technology. In: Lee C-C, editor. The current trends of optics and photonics. Netherlands, Dordrecht: Springer; 2015. p. 289–307.
- Lavrentovich OD, Kleman M, Bahr C. Cholesteric liquid crystals: defects and topology. In: Kitzerow H-S, editor. Chirality in liquid crystals. New York, NY: Springer; 2001. p. 115–158.
- Smalyukh II. Review: knots and other new topological effects in liquid crystals and colloids. Rep Prog Phys. 2020;83(10):106601. doi: 10.1088/1361-6633/abaa39
- Lee T, Trivedi RP, Smalyukh II. Multimodal nonlinear optical polarizing microscopy of long-range molecular order in liquid crystals. Opt Lett. 2010;35(20):3447. doi: 10.1364/OL.35.003447
- Derrick GH. Comments on nonlinear wave equations as models for elementary particles. J Math Phys. 1964;5(9):1252. doi: 10.1063/1.1704233
- Hobart RH. On the instability of a class of unitary field models. Proc Phys Soc. 1963;82(2):201. doi: 10.1088/0370-1328/82/2/306
- Bogdanov A, Hubert A. The stability of vortex-like structures in uniaxial ferromagnets. J Magn Magn Mater. 1999;195(1):182. doi: 10.1016/S0304-8853(98)01038-5
- Bouligand Y, Derrida B, Poenaru V, et al. Distortions with double topological character: the case of cholesterics. J Phys France. 1978;39(8):863. doi: 10.1051/jphys:01978003908086300
- Mertelj A, Lisjak D, Drofenik M, et al. Ferromagnetism in suspensions of magnetic platelets in liquid crystal. Nature. 2013;504(7479):237–241. doi: 10.1038/nature12863
- Ravnik M, Žumer S. Landau–de Gennes modelling of nematic liquid crystal colloids. Liq Cryst. 2009;36(10–11):1201. doi: 10.1080/02678290903056095
- Tai J-SB, Smalyukh II. Super-resolution stimulated emission depletion microscopy of director structures in liquid crystals. Opt Lett. 2018;43(20):5158. doi: 10.1364/OL.43.005158
- Kléman M. Points, lines, and walls: in liquid crystals, magnetic systems, and various ordered media. Hoboken, New Jersey: J. Wiley; 1983.
- Parkin SSP, Hayashi M, Thomas L. Magnetic domain-wall racetrack memory. Science. 2008;320(5873):190. doi: 10.1126/science.1145799
- Oswald P, Baudry J, Pirkl S. Static and dynamic properties of cholesteric fingers in electric field. Phys Rep. 2000;337(1–2):67. doi: 10.1016/S0370-1573(00)00056-9
- Zhao H, Tai J-SB, Wu J-S, et al. Liquid crystal defect structures with möbius strip topology. Nat Phys. 2023;19(3). 451–459. doi: 10.1038/s41567-022-01851-1
- Tai J-SB, Hess AJ, Wu J-S, et al. Field-controlled dynamics of skyrmions and monopoles. Sci Adv. 2024;10:eadj9373. doi: 10.1126/sciadv.adj9373
- Smalyukh II, Lansac Y, Clark NA, et al. Three-dimensional structure and multistable optical switching of triple-twisted particle-like excitations in anisotropic fluids. Nat Mater. 2010;9(2):139–145. doi: 10.1038/nmat2592
- Leonov AO, Inoue K. Homogeneous and heterogeneous nucleation of skyrmions in thin layers of cubic helimagnets. Phys Rev B. 2018;98(5):054404. doi: 10.1103/PhysRevB.98.054404
- Liu Y, Lake RK, Zang J. Binding a hopfion in a chiral magnet nanodisk. Phys Rev B. 2018;98(17):174437. doi: 10.1103/PhysRevB.98.174437
- Zheng F, Rybakov FN, Borisov AB, et al. Experimental observation of chiral magnetic bobbers in B20-type FeGe. Nat Nanotechnol. 2018;13(6):451–455. doi: 10.1038/s41565-018-0093-3
- Milde P, Köhler D, Seidel J, et al. Unwinding of a skyrmion lattice by magnetic monopoles. Science. 2013;340(6136):1076. doi: 10.1126/science.1234657
- Ackerman PJ, van de Lagemaat J, Smalyukh II. Self-assembly and electrostriction of arrays and chains of hopfion particles in chiral liquid crystals. Nat Commun. 2015;6(1). doi: 10.1038/ncomms7012
- Rybakov FN, Kiselev NS. Chiral magnetic skyrmions with arbitrary topological charge. Phys Rev B. 2019;99(6):064437. doi: 10.1103/PhysRevB.99.064437
- Kind C, Foster D. Magnetic skyrmion binning. Phys Rev B. 2021;103(10):L100413. doi: 10.1103/PhysRevB.103.L100413
- Hopf H. Über die Abbildungen der dreidimensionalen Sphäre auf die Kugelfläche. Math Ann. 1931;104:637. doi:10.1007/BF01457962
- Faddeev L, Niemi AJ. Stable knot-like structures in classical field theory. Nature. 1997;387(6628):58–61. doi: 10.1038/387058a0
- Battye RA, Sutcliffe PM. KnoTs as stable soliton solutions in a three-dimensional classical field theory. Phys Rev Lett. 1998;81(22):4798. doi: 10.1103/PhysRevLett.81.4798
- Gladikowski J, Hellmund M. Static solitons with Nonzero Hopf Number. Phys Rev D. 1997;56(8):5194. doi: 10.1103/PhysRevD.56.5194
- Hall DS, Ray MW, Tiurev K, et al. Tying quantum knots. Nat Phys. 2016;12(5):478–483. doi: 10.1038/nphys3624
- Acus A, Norvaišas E, Shnir Y. Hopfions interaction from the viewpoint of the product Ansatz. Phys Lett B. 2014;733:15. doi:10.1016/j.physletb.2014.04.009
- Kobayashi M, Nitta M. Torus Knots as Hopfions. Phys Lett B. 2014;728:314. doi:10.1016/j.physletb.2013.12.002
- Ackerman PJ, Smalyukh II. Diversity of knot solitons in liquid crystals manifested by linking of preimages in torons and hopfions. Phys Rev X. 2017;7(1):011006. doi: 10.1103/PhysRevX.7.011006
- Kent N, Reynolds N, Raftrey D, et al. Creation and observation of hopfions in magnetic multilayer systems. Nat Commun. 2021;12(1). doi: 10.1038/s41467-021-21846-5
- Bott R, Tu LW. Differential forms in algebraic topology. Berlin: Springer Science & Business Media; 2013.
- Whitehead JHC, An expression of Hopf’s invariant as an Integral, Proc Natl Acad Sci. 1947; 33:117.
- Tai J-SB, Ackerman PJ, Smalyukh II, Topological transformations of Hopf Solitons in chiral ferromagnets and liquid crystals. Proc Natl Acad Sci. 2018. 115:921. doi: 10.1073/pnas.1716887115
- Tai J-SB, Wu J-S, Smalyukh II. Geometric transformation and three-dimensional hopping of Hopf Solitons. Nat Commun. 2022;13:2986. doi: 10.1038/s41467-022-30494-2
- Poulin P, Stark H, Lubensky TC, et al. Novel colloidal interactions in anisotropic fluids. Science. 1997;275(5307):1770. doi: 10.1126/science.275.5307.1770
- Mundoor H, Senyuk B, Smalyukh II. Triclinic nematic colloidal crystals from competing elastic and electrostatic interactions. Science. 2016;352(6281):69. doi: 10.1126/science.aaf0801
- Babakhanova G, Parsouzi Z, Paladugu S, et al. Elastic and viscous properties of the nematic dimer CB7CB. Phys Rev E. 2017;96(6):062704. doi: 10.1103/PhysRevE.96.062704
- Tuchband MR, Shuai M, Graber KA, et al. The twist-bend nematic phase of bent mesogenic dimer CB7CB and its mixtures. arXiv:1511.07523.2015.
- Meyer RB. Effects of electric and magnetic fields on the structure of cholesteric liquid crystals. Appl Phys Lett. 2003;12(9):281. doi: 10.1063/1.1651992
- Ackerman PJ, Smalyukh II. Static three-dimensional Topological solitons in fluid chiral ferromagnets and colloids. Nat Mater. 2017;16(4):426. doi: 10.1038/nmat4826
- Tai J-SB, Smalyukh II. Static hopf solitons and knotted emergent fields in solid-state noncentrosymmetric magnetic nanostructures. Phys Rev Lett. 2018;121(18):187201. doi: 10.1103/PhysRevLett.121.187201
- Ackerman PJ, Boyle T, Smalyukh II. Squirming motion of baby skyrmions in nematic fluids. Nat Commun. 2017;8(1). doi: 10.1038/s41467-017-00659-5
- Han JH. Skyrmions in condensed matter. Springer; 2017.
- Leonov AO. Surface anchoring as a control parameter for shaping skyrmion or toron properties in thin layers of chiral nematic liquid crystals and noncentrosymmetric magnets. Phys Rev E. 2021;104(4):044701. doi: 10.1103/PhysRevE.104.044701
- Voinescu R, Tai J-SB, Smalyukh II. Hopf solitons in helical and conical backgrounds of chiral magnetic solids. Phys Rev Lett. 2020;125:057201. doi:10.1103/PhysRevLett.125.057201
- Sutcliffe P. Hopfions in chiral magnets. J Phys A Math Theor. 2018;51(37):375401. doi: 10.1088/1751-8121/aad521
- Nagaosa N, Tokura Y. Topological properties and dynamics of magnetic skyrmions. Nat Nanotech. 2013;8(12):899–911. doi: 10.1038/nnano.2013.243
- Mai Z, Yuan Y, Tai JSB, et al. Nematic Order, Plasmonic Switching and Self‐Patterning of Colloidal Gold Bipyramids. Adv Sci. 2021;8(22):2102854. doi: 10.1002/advs.202102854
- Hess AJ, Poy G, Tai J-SB, et al. Control of light by topological solitons in soft chiral birefringent media. Phys Rev X. 2020;10(3):031042. doi: 10.1103/PhysRevX.10.031042
- Poy G, Hess AJ, Seracuse AJ, et al. Interaction and Co-assembly of optical and topological solitons. Nat Photon. 2022;1(6):454–461. doi: 10.1038/s41566-022-01002-1
- Liu Q, Ackerman PJ, Lubensky TC, et al., Biaxial ferromagnetic liquid crystal colloids. Proc Natl Acad Sci. 2016;113:10479. doi: 10.1073/pnas.1601235113
- Mundoor H, Park S, Senyuk B, et al. Hybrid molecular-colloidal liquid crystals. Science. 2018;360(6390):768. doi: 10.1126/science.aap9359
- Mundoor H, Wu J-S, Wensink HH, et al. Thermally reconfigurable monoclinic nematic colloidal fluids. Nature. 2021;590(7845):268–274. doi: 10.1038/s41586-021-03249-0
- Chen X, Korblova, E, Dong, D, et al., First-principles experimental demonstration of ferroelectricity in a thermotropic nematic liquid crystal: polar domains and striking electro-optics. Proc Natl Acad Sci. 2020; 117:14021. doi: 10.1073/pnas.2002290117
- Saw TB, Doostmohammadi A, Nier V, et al. Topological defects in epithelia govern cell death and extrusion. Nature. 2017;544(7649):212–216. doi: 10.1038/nature21718
- Maroudas-Sacks Y, Garion L, Shani-Zerbib L, et al. Topological defects in the nematic order of actin fibres as organization centres of hydra morphogenesis. Nat Phys. 2021;17(2):251–259. doi: 10.1038/s41567-020-01083-1
- Yaman YI, Demir E, Vetter R, et al. Emergence of active nematics in chaining bacterial biofilms. Nat Commun. 2019;10(1). doi: 10.1038/s41467-019-10311-z
- Nijjer J, Li C, Kothari M, et al. Biofilms as self-shaping growing nematics. Nat Phys. 2023;1(12):1936–1944. doi: 10.1038/s41567-023-02221-1