## Kinetics of motile solitons in fluid nematics

Satoshi Aya  ,# , Fumito Araoka 

Physicochemical Soft Matter Research Team, RIKEN Center for Emergent Matter Science (CEMS), Japan

Solitary waves, dubbed 'solitons', are special types of waves that propagate for an infinite  distance  under  ideal  conditions.  These  waves  are  ubiquitously  found  in nature such as typhoon or neuron signals. Yet, their artificial generation and the control  of  their  propagation  remain  outstanding  challenges  in  materials  science owing  to an  insufficient understanding  of  the experimental  conditions and theoretical aspects. Herein, a generic strategy for forming particle-like solitons and controlling  their  kinetics  in  nematic  fluid  media  is  reported.  The  key  to  the realisation  of  the  generation  of  solitons  and  the  control  of  their  kinetics  is  the coupling between the fluid elasticity and the background flow flux, as evidenced by experimental observations and theoretical approaches. The findings of this study enable the exploration of solitons in a wide range of materials and have technological ramifications for the lossless transport of energy or structures.

Solitons are packets of spatially localised moving waves that convey energy with very little or no mass transport. Since the first observations of the solitonic behaviour of water waves by Russell in 1834 and the introduction of the term of 'soliton' by Zabusky and Kruskal 1 , solitons have been identified in various natural phenomena. However, scientists

 To whom correspondence should be addressed, E-mail: satoshi.aya@riken.jp, fumito.araoka@riken.jp

# Current contact address: South China University of Technology, satoshiaya@scut.edu.cn

have little access to artificial solitons, despite many  dedicated theoretical and computational studies. From a mathematical perspective, solitons are recognised by their particle-like nature such as nondissipative and elastic collisional properties as a result of balancing  their  dispersive  effects  with  the  nonlinear  counterparts  described  by  a Korteweg-de-Vries or nonlinear Schrö dinger equation . This nature of solitons enables 2 the  distortion-free  long-distance  transport  of  waves  or  structures,  which  has  attracted considerable  interest  from  both  fundamental  and  technological viewpoints  in  many branches of science, including atmospheric circulation 3,4 , fluid  and mechanical waves 5 68 organisms 9,10 condensates 11-13 and

- , living , Bose-Einstein , light confinement, propagation 14-16 .

Very  recently,  there  have  been  several  groundbreaking  studies  observing solitonic analogues for physical systems such as spintronic 17-19 and fluid systems 20-24 . In most of these cases, the solitons possess topological singularities in spin-like vector fields that  are  distinct  from  the  far-field  background.  In  other  words,  the  structures  of  the solitons  are  topologically  protected.  Although  it  is  expected  that  such  phenomena  are generally triggered by the spontaneous breaking of either the rotational or translational symmetry or both, the primary unsolved questions lie in the understanding of the physics underpinning both the generation and transport of solitons in real materials.

In  this  article,  we  describe  many  unknown  dynamic  solitonic  phenomena closely involved in many systems with a particle-like nature that can be directly generated and  even  controlled  in  nematic  fluids.  A  swarm  of  a  novel  type  of  motile  solitons  is stabilised and 'fueled' by the coupling of the dielectric and hydrodynamic convectional responses  to  low-voltage  electric  fields.  This  discovery  establishes  the  fundamental physics and strategy of manipulating the kinetics of solitons in an experimental setting and permits the observation of their time evolution, facilitating further experimental and theoretical studies of dynamic soliton excitations in soft-matter systems.

## Soliton creation

As a soliton medium, we prepare thin films of frustrative mixtures comprising two nematics, E7 and 4'-butyl-4-heptyl-bicyclohexyl-4-carbonitrile (CCN47) (Fig. 1a). We uniformly align them between two glass plates with the average molecular orientation, denoted  by  the  director 𝒏 ,  parallel  to  the  plates.  To  create  solitons,  a  sinusoidal  or rectangular  AC  electric  field  is  applied  normal  to  the  plates  (Fig.  1a).  The  electric response significantly changes because of the frustration of the dielectric response and the  hydrodynamic  instability  incurred  by  changing  the  combination  of  the  dielectric anisotropy, Δε  = ε∥ -ε⊥ , and conductivity anisotropy, Δσ = σ∥ -σ⊥ (the subscripts ∥ and ⊥ indicate the directions parallel and perpendicular to 𝒏 , respectively). Generally, the  combination  of  dielectricity  and  conductivity  can  be  classified  using  de  Gennes' notation  as (ΔεΔσ) 25 .  As  E7  and  CCN47  exhibit  positive (++) and  negative (--) anisotropies,  respectively,  mixing  E7  with  CCN47  enables  smooth  sign  tuning  from (--) to (++) via (-+) .  As  the  mixture  of  anisotropies  corresponds  to  specific parameter strengths, the effects of this tuning can be investigated.

In Figs. 1b-d, we show time-dependent microscopic texture images acquired during soliton creation in a mixture containing 16 wt% E7 for a 10-V electric field at 12 Hz at the specific times designated in Fig. 1e. Two processes-nucleation and packingare clearly visible.

Nucleation . Independent cruciform solitons are randomly excited in space and time (Fig. 1b) as a result of the localisation of ions when an electric field is applied, while they disappear when the field is absent (Supplementary section 1 and 2, Supplementary Figs. S1 and S2). These solitons fluctuate around and remain relatively near their nucleation positions  with  a  space-filling  ratio  of 𝑋 = 𝜋𝑅 2 /𝐴 &lt; 0.6, where 𝑅 and 𝐴 are  the soliton radius and the total area of the region, respectively. The time dependence of X is fitted by the Kolmogorov-Johnson-Mehl-Avrami (KJMA) expression, which might give

insights to the procedures of the migration and the localisation of the ions during the soliton generation (Supplementary Fig. S3), ln[ln 1/{1 - 𝑋}] = 𝑚 ln 𝑡 + ln 𝐾 , where 𝑋 , 𝑡 , 𝐾 , and 𝑚 are  the  space  occupancy  of  the  director-deformed  region,  the  time,  the temperature-dependent Avrami coefficient, and the exponent, respectively 26 . The value obtained  for 𝑚 is  2.9,  suggesting  the  involvement  of  two-dimensional  homogeneous nucleation and growth processes triggered by the localisation of ions 27 .

Packing. Centered rectangular packing of the solitons occurs when the filling ratio 𝑋 &gt; 0.6 (Fig. 1d, Supplementary Fig. S4), a value consistent with the close-packed ratio of perfectly hexagonally packed solid spheres in two dimensions (0.604). For such packing states,  the  solitons  are  fusion-resistant.  Similar  packing  structures  occur  in  some  selfdriven systems 24,28,29

Structure of solitons. We investigate the director structure of the solitons.  Since the solitons are dynamic with time, magnified views of solitons obtained with an exposure time of 200 ms  under crossed polarisers and with a standard tint waveplate (retardation = 530 nm) taken by a conventional camera (Figs. 1f and g) just reflect a timeaveraged structure, an array of concentric director looking like patterns with one s =+1 topological defect is clear in the midplane of the nematic film. In order to elucidate the dynamic structure, we perform high-frame-rate polarising microscopy at 100 fps. Figure 1h shows a time-dependent oscillation of the structure of the solitons. Surprisingly, unlike the static structure observed in Figs. 1f and g, the director field of the solitons transforms between two defect-less structures, i.e., biconvex and biconcave structures, synchronised with  the  applied  electric  field.  It  is  noted  that,  by  analysing  the  director  fields,  the biconvex structure has  a larger deviation angle  of  the director in  the  midplane  of the nematic film than the biconcave structure (Supplementary Fig. S5). This explains why the time-averaged structure looks like a concentric domain with a s =+1 defect.

As the soliton filling ratio increases, mutual repulsion, separating the solitons by distances of about twice the soliton radius, is observed. As a result of this repulsion, the solitons become resistant to the external stress without coalescence, confirming the pseudoparticle nature of the solitons.

Conditions of soliton creation: topological state diagram. The topological pattern that includes the solitons varies depending on the voltage and frequency (Figs. 2a-e). On the basis of the director deformations and flow profiles (Figs. 2f-j), we classify the frustrated states into five categories: (I) ground planar alignment (G-state), (II) ∥ -roll pattern ( ∥ -state), (III) rectangular pattern (R-state), (IV) soliton state (Soliton-state), and (V) ⊥ -roll pattern ( ⊥ -state).  Fig.  2k  shows  a  state  diagram  for  a  mixture  containing  16  wt%  E7, which exhibits (-+) anisotropy.

In  the  G-state,  neither  elastic  deformation  nor  hydrodynamic  effects  occur; therefore,  the  initial  uniform  uniaxial  alignment  is  maintained.  In  the ∥ -state,  which appears  in  the  lower-frequency  region  below  9  Hz,  only  an  in-plane  splay  elastic deformation  of  the  director  occurs  because  of  the  coupling  with  the  flexion-induced polarisation,  the  so-called  flexoelectric  effect 30 .  Because  the  strong  surface  anchoring effect  forces  fix  the  surface  directors  in  the  initial  state,  twist  deformation  along  the electric  field  direction  ( z direction)  is  induced,  yielding  a  one-dimensional  periodic pattern. Particle tracking confirms the absence of flow (see Methods). An increase in the voltage in the ∥ -state causes a field-induced hydrodynamic convective flow parallel to the initial director ( 𝒏 0). The coupling between the flow and the flexoelectric effect induces a  two-dimensional  version  of  the ∥ -state,  i.e.,  the  R-state.  The ⊥ -state  occurs  in  a dielectric  regime  in  which  ions  are  effectively  immobile,  but  the  nematic  director oscillates with the applied electric field. The convective flow field parallel to 𝒏 0 arising from an electrohydrodynamic effect (Carr-Helfrich instability) yields a stripe pattern in which the director has a wavy modulation in space, as observed in nematics with (-+)

anisotropy 26 .  The  appearances  of  the  respective  roll  patterns  in  the ∥ -  and ⊥ -states significantly differ in terms of the roll direction: rolls in the ∥ -state run parallel to the director; those in the ⊥ -state are perpendicular. The Soliton-state occurs in the frequency range  between  the ∥ -/R-  and ⊥ -states  and  spans  a  conductive  regime  in  which  ionic convection dominates the dielectric oscillation of the director 31 . Accordingly, the soliton system  has  features  of  both  the ∥ and ⊥ -states  with  a  compresence  of  elastic  and convective effects, as discussed below.

We next turn to the major focus of this study: how the controlled anisotropies of dielectricity and conductivity affect the soliton stabilisation mechanism. To systematically  examine  the  effects  of  the  material  parameters,  we  either  altered  the concentration of E7 or added excesses of the ionic species TBABE to the CCN47. Fig. 2l shows a generalised scheme of the procedural conditions. Solitons are observed only in (-+) mixtures at low electric fields ( ≦ 10 V), a typical condition in which the CarrHelfrich electrohydrodynamic mechanism, where a convective flow occurs parallel to the background  director,  becomes  active 25 . In  contrast,  no  observable  flow  occurs  for (- -) and (+ +) mixtures, with the flexoelectric effect appearing for the former and Freederikz-effect-driven  uniform  reorientation  appearing  for  the  latter.  Although  the present system does not realise (+ -) mixtures, the Freederikz effect dominates (+ -) mixtures  in  conventional  systems.  The  strength  of  the  conductivity  is  also  vital  in determining the soliton stability. Fig. 2m shows that solitons appear in the intermediate strength  range 8 × 10 -9 &lt; σ &lt; 4 × 10 -8  Ω -1 ∙ m -1 ,  where  the  localisation  of  ions  is realised. This once again confirms that the Carr-Helfrich electrohydrodynamic mechanism is necessary for the presence of solitons (Supplementary Information 8).

## Soliton kinetics

The soliton system described herein exhibits two types of kinetic regimesswimming and proliferation-which appear in ascending order as the frequency of the

Soliton-state  increases  at  a  fixed  voltage.  Worth  noting  is  that  voltage  or  frequency contributes differently to the kinetics of the solitons as discussed below.

Swimming  regime:  collision  and  elastic  reflection. An  increase  in  the  voltage  or frequency in the oscillation regime (the transition from the orange to purple-blue regions in Supplementary Fig. S6) produces steady directional motion of the solitons in which the initial cruciform configuration changes to a bug-eye-like pattern because of asymmetric elastic  deformation  (Figs.  3a  and  b).  The  direction  of  this  motion  can  be  tuned  by changing either the amplitude or frequency of the electric field. Figs. 3c-h show the timedependent trajectories of solitons under the control of the tuned amplitude of the electric field at a characteristic frequency of 20 Hz, where the solitons are most stable in their shape. At low voltages (8.2 &lt; V &lt; 8.5 V, Fig. 3c), the solitons move directionally parallel to 𝒏 0 in the xy plane. Notably, they randomly move to the left or right because of the previously  discussed  dielectric  oscillation.  As  plotted  in  Fig.  3i,  at V &gt;  8.5  V,  the directional  angle  of  motion  with  respect  to 𝒏 0  continuously  varies  as  the  voltage  is increased (Figs. 3d-h), corresponding to the increasingly asymmetric deformation of the bug-eye-like pattern (Fig. 3b). High-speed camera observations reveal that this movement is  driven  by  the  oscillation  of  the  director  field  between  distorted  versions  of  two orientational states, as shown in Fig. 1h. The directional motion ends at about 10 V, after which the soliton motion becomes somewhat chaotic, leading to a fractalised proliferation process that will be discussed later. The average velocity of the motion of the solitons is measured through a trajectory tracking as shown in Figs. 3c-d and plotted as a function of voltage and frequency in Figs. 3j, k. Although the velocity is nearly constant with respect  to  the  voltage  because  the  background  flow  field  is  mostly  unaffected  in  the limited voltage range, a significant frequency dependence is observed. The solitons are resonantly active at about 20 Hz, as shown in Fig. 3k. This frequency corresponds to a critical frequency that separates the conductive and dielectric regimes 32 , as described by the  competition  between  the  relaxation  process  of  mobile  charge  in  space  and  the

dielectric relaxation of the director, 𝜎 // ɛ0ɛ//. Given 𝜎 // = 10 -9 1/Ω∙m and ɛ// = 4, as measured experimentally, we obtain the critical frequency of 28 Hz by calculation, consistent with the experimentally observed value.

The dynamics of collision between solitons and their reflection by obstacles can be understood in terms of their pseudoparticle nature. Figs. 4a-c show that their head-on collisional behaviour varies depending on the degree of offset, ∆, between soliton cores. When R &lt; ∆ &lt; 2 R , interference occurs between the peripheries of soliton pairs with the same orientation. Consequently, the solitons repel each other at a small post-collision angle;  at  0 ≦ ∆  &lt; R ,  elastic  deformation  becomes  more  significant  because  of  the increased mismatch of the director field between solitons, leading to vertical repulsion between  the  solitons.  Following  this  collision,  the  solitons  separate  such  that  elastic deformation  is  avoided  and  then  continue  to  move  along  their  original  directions  of motion. Interestingly, the solitons are elastically reflected at an air-nematic interface in the manner of real particles colliding with a wall (Fig. 4d).

Proliferation  regime. Figs.  5a  and  b  show  the  fractalised  proliferation  trajectory  of destabilised solitons for a 10-V field at 40 Hz. Soliton fractalisation produces clones for movement  perpendicular  to 𝒏 0.  Upon  fractalisation,  the  solitons  bifurcate  and  move obliquely for a short distance. The fractalisation process corresponds to the continuous growth in soliton size. Considering that an optimum size exists for the solitons at a fixed film thickness  (Supplementary  Fig.  S7),  the  fractalisation  process results in the continuous accumulation of excess elastic energy. The solitons fractalise when the energy penalty  from  elastic  deformation  cannot  be  compensated  by  the  effects  of  surface anchoring and dielectric interaction. This is illustrated in Fig. 5c, which shows the soliton diameter, 2 R , as a function of time; fractalisation occurs when the soliton size exceeds a specific value. Notably, the fractalisation simultaneously occurs for all solitons (Figs. 5a and b). This implies that a synchronised growth mechanism is involved for the solitons,

resulting in a solitonic trajectory with a topological Hausdorff fractal dimension of D = 1.

## Discussion

Soliton kinetics. We first discuss the coupling of the far-field background flow field with elastic  deformation  in  the  Soliton-state  that  drives  the  kinetics  of  the  solitons.  In  the Soliton-state, convective flows with waves obliquely traveling to both the left and right relative to 𝒏 0  simultaneously appear, as confirmed by the observation of roll patterns coexisting with solitons in the ⊥ -state. The localised dynamics are assumed to consist of two oblique waves traveling at an angle 𝜃 with respect to 𝒏 0, as introduced elsewhere 3335 .  As  a  first-order  approximation,  the  amplitude  expansion  of  the  flow  field  is 𝑤 = [𝐴(𝑟, 𝑡)𝑒 𝑖𝑘𝐴∙𝑟 +𝐵(𝑟,𝑡)𝑒 -𝑖𝑘𝐵∙𝑟 ]𝑒 (𝑖𝑘 𝐶 ∙𝑟-𝜔𝑡) sin 𝑧 , where 𝐴 and 𝐵 are  the  respective spatiotemporal variable amplitudes of the two assumed traveling waves; 𝑘𝐴 , 𝑘𝐵 , and 𝑘𝐶 are  the  wave  numbers; 𝑟 is  the  spatial  coordinate; 𝑡 is  time;  and 𝜔 is  the  angular frequency.  For  small  amplitudes 𝐴(= 𝐴0𝑒 -𝑖(𝑄𝑥+𝑃𝑦-Ω𝑡) ) and 𝐵(= 𝐵0𝑒 -𝑖(𝑄𝑥-𝑃𝑦-Ω𝑡) ) , the description of the flow waves in the Ginzburg-Landau framework is valid, yielding 𝜕 𝜕𝑡 𝐴 = 𝜕 𝜕𝐴 [-∫𝑑𝑉(𝒖𝐴 ∙ ∇𝐴) + (𝜇 - 𝑐|𝐴| 2 -𝑔|𝐵| 2 )𝐴 2 +(∇𝐴) 2 ] and 𝜕 𝜕𝑡 𝐵 = 𝜕 𝜕𝐵 [-∫𝑑𝑉( 𝒖𝐵 ∙ ∇𝐵) + (𝜇 - 𝑐|𝐵| 2 -𝑔|𝐴| 2 )𝐵 2 +(∇𝐵) 2 ] , where 𝒖𝐴,𝐵 are  the  group velocities  of  the  waves, 𝜇 is  the  bifurcation  parameter,  and 𝑐 and 𝑔 are  coupling coefficients. Forming a symmetric superposition mode of the two waves in the xy plane with a finite Ω proportional to |𝐴0 | 2 produces a two-dimensional localised traveling rectangle solution corresponding to a solitonic state. If only the flow field is under action, two types of solitonic motion are possible: a directional motion parallel to 𝒏 0 and a twodimensional random walk within the xy plane. However, in our results, variable angles of solitonic motion with respect to 𝒏 0 are achieved. This motivates us to consider the elastic deformation  of  the  bulk  directors  by  introducing  a  spatially  modulated  imbalance  of amplitudes between waves 𝐴 and 𝐵 .  Although convective flow is considered in Ref.

[33], our solitons have splay deformation in the xy plane and twist-bend deformation in the xz plane  (Fig.  2i)  in  addition  to  a  convective  flow.  Owing  to  the  flow-alignment property of the director, its orientation tends to be parallel to the flow field. Since the orientation of the director oscillates with time, the local flow also changes its direction accordingly through a self-regulating procedure (Supplementary Fig. S9). Such processes can  trigger  the  local  motion  of  topological  defects 36 ,  thereby  providing  an  additional traveling mode out of the xz plane. Indeed, the electrohydrodynamic convection flow can be  modulated  by  twist  deformation  to  produce  a  spiral  flow 37 .  Because  the  elastic deformation  in  the xy plane  weakens  in  the  transition  from  the ∥ -state  through  the Soliton-state to the ⊥ -roll state as the frequency is increased or the voltage is decreased, the  inclination  of  the  oblique  solitonic  movement  decreases  with  respect  to 𝒏 0;  with neither  flow  nor  elastic  deformation  in  the xy plane,  no  two-dimensional  motion  is observed in the studied parameters (Figs. 2a-e). Thus, the local coupling between the elastic  deformation  and  the  flow  field  in  an  anisotropic  medium  permits  effective dictation of the solitonic motion directionality, which cannot be achieved through a flow field alone.

Pseudoparticle soliton properties. Herein, we discuss the pseudoparticle properties of solitons  in  anisotropic  flows.  For  characterisation,  we  analyse  a  radial  distribution function of packed solitons, 𝑔(𝑟 cc ) (Supplementary Fig. S9) and the thermal-energyactivated solitonic Brownian diffusion. Fig. 6a shows the pairwise potential 𝑉(𝑟 cc ) as a function of the center-to-center distance between solitons, 𝑟 𝑐𝑐 . At the dilute limit, 𝑉(𝑟 cc ) is directly calculated as 𝑉(𝑟 cc ) 𝑘B𝑇 = -ln[𝑔(𝑟 cc )] . The first potential closely corresponds to the mean nearest-neighbor distance between solitons, 𝑟 cc,o = 15.2 µm. The depth of the potential  well  slightly  exceeds  the  thermal  energy, 𝑘𝐵𝑇 ,  which  is  consistent  with  the observed  minimal  soliton  fluctuation  beyond 𝑟 cc,o .  The  solitonic  Brownian  diffusion process is calculated as the two-dimensional mean square displacement of the solitons over  time 𝑡 (Fig.  6b).  The  effective  diffusion  constant  of  the  solitons  is  1.4  µ m /s, 2

corresponding to a mean solitonic fluctuation velocity of about 1 µ m/s. The effective mass of the solitons, 𝑚𝑒𝑓𝑓 ,  is  expressed  as 𝑈 = 1/2𝑚𝑒𝑓𝑓𝑣 2 ,  where 𝑈 is  the  kinetic energy and 𝑣 the mean velocity of solitonic fluctuation. By balancing the thermal energy, 𝑘B𝑇 , and the kinetic energy of the solitons, 𝑚𝑒𝑓𝑓 is estimated to be 1 pg. In contract, the real mass of molecules constituting each soliton is estimated to be 𝜌𝜋𝑅 2 ℎ ≈ 1000 pg, where 𝜌 and ℎ are the material density and nematic film thickness, respectively. The effective mass of the soliton is significantly smaller than the real mass of the materials constituting the solitons. Notably, 𝑚𝑒𝑓𝑓 is independent of the film thickness irrespective of  the  change  in  the  real  mass  (Supplementary  Fig.  S12).  Because  the  elastic  energy density decreases as the film thickness increases, the total elastic energy does not vary. This means that 𝑚𝑒𝑓𝑓 reflects the elastic energy stored by the elastic deformation in the director field.

## Summary

We  developed  and  demonstrated  a  strategy  for  manipulating  the  complex spatiotemporal behaviour of solitons by fine-tuning several parameters upon which the coherent  interactions  between  the  traveling  flow  and  elastic  deformation  are  highly dependent. This system provides a model for observing and controlling soliton dynamics and suggests that solitary behaviours reflect a more universal nature applicable to other soft-matter systems. Further work is required to determine whether additional parameters such as chirality and local ordering drive new solitonic behavioural features.

## References

- 1. Zabusky, N. J. &amp; Kruskal, M. D. Interaction of solitons in a collisionless plasma and recurrence of initial states. Phys. Rev. Lett. 15 , 240-243 (1965).
- 2. Manton, N. &amp; Sutcliffe, P. Topological Solitons (Cambridge Univ. Press, 2004).
- 3. Doviak, R. J. &amp; Ge, R. S. An atmospheric solitary gust observed with a Doppler radar, a tall tower and a surface network. J. Atmos. Sci. 41 , 2559-2573 (1984).
- 4. Doviak, R. J. &amp; Christie, D. R. Thunderstorm-generated solitary waves - A wind shear hazard. J. Aircr. 26 , 423-431 (1989).
- 5. Lakshmanan, M. Solitons, tsunamis and oceanographical applications of. in Extreme

Environmental Events (ed. Meyers, R. A.) 873-888 (Springer, 2011).

- 6. Deng,  B.,  Raney,  J.  R.,  Tournat,  V .  &amp;  Bertoldi,  K.  Elastic  vector  solitons  in  soft architected materials. Phys. Rev. Lett. 118(20) , 204102 (2017).
- 7. Alden, J. S. et al. Strain solitons and topological defects in bilayer graphene. Proc. Natl. Acad. Sci. 110 , 11256-11260 (2013).
- 8. Chen,  B.  G.,  Upadhyaya,  N.  &amp; Vitelli,  V .  Nonlinear  conduction  via  solitons  in  a topological mechanical insulator. Proc. Natl. Acad. Sci. 111 , 13004-13009 (2014).
- 9. Kuwayama, H. &amp; Ishida, S. Biological soliton in multicellular movement. Sci. Rep. 3 , 2272 (2013).
- 10. Ward, A. et al. Solid friction between soft filaments. Nat. Mater. 14 , 583-588 (2015).
- 11. Lenz, G., Meystre, P. &amp; Wright, E. M. Nonlinear atom optics. Phys. Rev. Lett. 71 , 3271-3274 (1993).
- 12. Denschlag,  J.  et  al.  Generating  solitons  by  phase  engineering  of  a  Bose-Einstein condensate. Science 287 , 97-100 (2000).
- 13. Khaykovich, L. et al. Formation of a matter-wave bright soliton. Science 296 , 12901293 (2002).
- 14. Mitchell, M. &amp; Segev, M. Self-trapping of incoherent white light. Nature 387 , 880883 (1997).
- 15. Piccardi, A., Alberucci, A., Bortolozzo, U., Residori, S. &amp; Assanto, G. Soliton gating and switching in liquid crystal light valve. Appl. Phys. Lett. 96(7) , 071104 (2010).
- 16. Assanto,  G.,  Fratalocchi,  A.  &amp;  Peccianti,  M.  Spatial  solitons  in  nematic  liquid crystals: from bulk to discrete. Opt. Express 15 , 5248-5259 (2007).
- 17. Mochizuki, M., Yu, X. Z., Seki, S., Kanazawa, N., Koshibae, W., Zang, J., Mostovoy, M.,  Tokura,  Y.  &amp;  Nagaosa,  N.  Thermally  driven  ratchet  motion  of  a  skyrmion microcrystal and topological magnon Hall effect. Nat. Mater. 13 , 241-246 (2014).
- 18. Okamura, Y., Kagawa, F., Seki &amp; Tokura, Y. Transition to and from the skyrmion lattice  phase  by  electric  fields  in  a  magnetoelectric  compound. Nat.  Commun. 7 , 12669 doi: 10.1038/ncomms12669 (2016).
- 19. Takagi, R., White, J. S., Hayami, S., Arita, R., Honecker, D., Rønnow, H. M., Tokura, Y. &amp; Seki, S. Multipleq noncollinear magnetism in an itinerant hexagonal magnet. Sci. Adv. 4(11), eaau3402 (2018).
- 20. Ackerman, P. J. &amp; Smalyukh, I. I. Static three-dimensional topological solitons in fluid chiral ferromagnets and colloids. Nat. Mater. 16 , 426-432 (2017).
- 21. Ackerman, P. J., Boyle, T. &amp; Smalyukh, I. I. Squirming motion of baby skyrmions in nematic fluid. Nat. Commun. 8(1) , 673 (2017).
- 22. 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).
- 23. Li, B.-X. et al. Electrically driven three-dimensional solitary waves as director bullets in nematic liquid crystals. Nat. Commun. 9(1) , 2912 (2018).
- 24. Nych, A., Fukuda, J., Ognysta, U., Žumer, S. &amp; Muševič, I. Spontaneous formation and dynamics of half-skyrmions in a chiral liquid-crystal film. Nat. Phys. 13 , 12151220 (2017).
- 25. de Gennes, P. G. &amp; Prost, J. The Physics of Liquid Crystals . (Clarendon Press, Oxford, 1993).
- 26. Avrami, M. Kinetics of phase change. I. General theory. J. Chem. Phys. 7 , 1103-1112 (1939); Avrami, M. Kinetics of phase change. II. Transformation-time relations for random distribution of nuclei. J. Chem. Phys. 8 , 212-224 (1940); Avrami, M. Kinetics of phase change III. Granulation, phase change, and microstructure. J. Chem. Phys. 9 ,

177-184 (1941).

- 27. Torrens-Serra, J., Venkataraman, S., Stoica, M., Kuehn, U., Roth, S. &amp; Eckert, J. Nonisothermal kinetic analysis of the crystallization of metallic glasses using the master curve method. Materials 4 , 2231-2243 (2011).
- 28. Sumino, Y., Nagai, K. H., Shitaka, Y., Tanaka, D., Yoshikawa, K., Chate, H. &amp; Oiwa, K. Large-scale vortex lattice emerging from collectively moving microtubules. Nature 483 , 448-452 (2012).
- 29. Doostmohammadi, A., Adamer, M. F., Thampi, S. P. &amp; Yeomans, J. M. Stabilization of active matter by flow-vortex lattices and defect ordering. Nat. Commun. 7 , 10557 (2016).
- 30. Blinov, L. M. Domain instabilities in liquid crystals. J. Phys. Colloq. 40(C3) , C3-247 (1979).
- 31. Kramer, L. &amp; Pesch, W. Convection instabilities in nematic liquid crystals. Ann. Rev. Fluid. Mech. 27 , 515-539 (1995).
- 32. Éber, N., Salamon, P. &amp; Buka, Á. Electrically induced patterns in nematics and how to avoid them. Liq. Cry. Rev. 4 , 101-134 (2016).
- 33. Riecke, H. &amp; Granzow, G. D. Localization of waves without bistability: Worms in nematic electroconvection. Phys. Rev. Lett. 81 , 333-336 (1998).
- 34. Dennin, M., Ahlers, G. &amp; Cannell, D. S. Spatiotemporal chaos in electroconvection. Science 272 , 388-390 (1996).
- 35. Treiber, M. &amp; Kramer, L. Coupled complex Ginzburg-Landau equations for the weak electrolyte model of electroconvection. Phys. Rev. E 58 , 1973-1982 (1998).
- 36. Riecke, H. Localized structures in pattern-forming systems. in: Pattern Formation in Continuous  and  Coupled  Systems.  The  IMA  Volumes  in  Mathematics  and  its Applications (Eds. Golubitsky, M., Luss, D., Strogatz, S. H.) vol. 115 (Springer, New York, NY, 1999)
- 37. Hertrich,  A.,  Krekhov,  A.  &amp;  Scaldin,  O.  The  electrohydrodynamic  instability  in twisted nematic liquid crystals. J. Phys. II 4 , 239-252 (1994).

## Methods

Sample preparation. The frustrative nematic mixtures were formed by mixing E7 (a mixture of cyanobiphenyl and cyanoterphenyl components, Wako, Aldrich, and made in house) with 4ʹ-butyl-4-heptyl-bicyclohexyl-4-carbonitrile (CCN47, Nematel GmbH) by weight. These molecules form a uniaxial nematic phase in which the average molecular orientation is characterised by a vectorial director, n . Because of the macroscopic and nonpolar nature of the nematic phase, the head and tail of n are equivalent. The liquid crystals are uniformly aligned between two glass plates with a controlled thickness in the range of 2-50 µ m. In the manuscript, if not specified, we show solitonic behaviours

in a 5.2-µ m-thick film. Although E7 has positive anisotropies (𝚫𝛆𝚫𝛔) = (+ +) , those of CCN47 are negative (𝚫𝛆𝚫𝛔) = (- -) . Tetrabutylammonium benzoate (TBABE) was added to the E7-CCN47 mixture to tune its net strength of conductivity. Films comprising these mixtures were sandwiched between pairs of substrates coated with rubbed polyimide on indium-tin-oxide electrode layers, with the film thickness adjusted by adding silica particles to obtain a 2-12-µ m-thick film (micromer, Micromod) or by adding polyester film spacers to obtain a 20-50-µ m-thick film (Mylar, Dupont). In many other liquid crystalline systems, including Schiff bases and cyanobiphenyl mesogens, solitons were also found.

Conductivity measurement. The electrical and dielectric responses in the liquid crystal mixtures were characterised via dielectric spectroscopy (Solartron Analytical, impedance/gain-phase analyser 1260A with dielectric interface system 1296A). An AC electric field at 0.05 Vrms sweeping at a frequency of 10 -2 -10 4 Hz was applied to the films, and the resulting impedance response was recorded and analysed (Supplementary Fig. S13).

Director mapping. The orientation of the directors was determined optically using polarising microscopy. The highest birefringence of the mixtures was sufficiently low (0.03 &lt; 𝚫𝒏  = 𝒏𝐞 - 𝒏𝐨 &lt; 0.08, where 𝒏𝐞 𝐚𝐧𝐝 𝒏𝐨 are the extraordinary and ordinary refractive indices, respectively) to carry out this task. On the basis of the measured transmittance of light through the films, we calculated the orientation of the directors. Our method was the same as that presented in Ref. [21] . We also conducted fluorescence confocal polarising microscopy (TCS SP8 STED, Leica) to probe the three-dimensional distribution of the soliton directors. Although volume scanning was performed during soliton movement, owing to their fast fluctuation that was synchronised with the electric field, the solitons were not clear enough to allow their analysis. Instead, we measured the emission intensities at single planes on the surfaces

of the slab and at the middle of the cell. The intensity on the surfaces was nearly constant during the application of the electric field, confirming a strong surface anchoring effect. The signal was maximum at the middle the slab, corresponding to elastic deformation of nearly zero on the surface that reached a maximum at the middle of the slab.

Characterisation of flow field. The particle tracking method was used to visualise the effective flow in the mixtures in the presence of an electric field. In this process, a small number of 2-µ m particles were dispersed into the nematic mixtures, and their motions were tracked by both bright-field and polarising microscopy.

Time-dependent trajectory analysis. The recorded dynamics of the solitons were analysed using custom-made ImageJ macros. For the precise image analysis needed to extract sharp shapes and calculate the centers of gravity of the solitons, a dark but slightly structured time-dependent background was averaged over time.

## Acknowledgments

The authors thank H. Orihara (Hokkaido University), A. Buka (Hungarian Academy of Science), P. Salamon (Hungarian Academy of Science), H. Yoshida (Osaka University), I. I. Smalyukh (University of Colorado Boulder), and J. Yoshioka (Ritsumeikan University, ex. RIKEN CEMS) for their valuable input. We would like to thank Editage (www.editage.jp) for English language editing.

## Author contributions

S.A. and F.A. designed and directed the research, discussed the mechanisms, and wrote the manuscript. S.A. performed all the experiments and analyses.

## Additional information

Supplementary information is available in the online version of this paper. Reprints and permissions information is available online at www.nature.com/reprints. Correspondence and requests for materials should be addressed to S.A. or F.A.

## Competing financial interests

The authors declare no competing financial interests.

Figure 1. Creation of ordered solitons. (a) Slab geometry and chemical structures of the main materials-cyanobiphenyl (CB) and 4'-butyl-4-heptyl-bicyclohexyl-4-carbonitrile (CCN47). (b-d) Nucleation and growth processes of electrically pumped solitons visualised using polarising microscopy with crossed polarisers for specific times indicated by the time evolution of the soliton space-filling ratio in (e). Scale bar, 100 µ m. (e) Time evolution of the space-filling ratio of solitons. (f) Polarising microscopy images of packed solitons. The soliton radius, R , is defined as the radius of the primary area undergoing elastic deformation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-6be95c39-3773-4f56-9190-3219da892ab8.jpg)

The space-filling ratio X is defined as the ratio of the total area of N solitons to the sample area, A . (g) Counterpart to (f) taken with a tint plate. (h) Monochromatic textures with a 𝜆 /4 plate and schematics of the alternating director field between two orientational states, taken with a high-frame-rate camera.

Figure 2. Conditions for generating solitons. (a-e) Polarising microscope textures of the topological states in the presence of an electric field. (f-j) Schematics of the director field by nail notation and motion trajectories of particle motion. The director field in the xy plane is on the middle plane of the cell. (k) State diagram as a function of the voltage and frequency. (l) Primary phenomena occurring because of the changes in the signs of the anisotropies of dielectricity (Δε) and conductivity (Δσ). The solid line indicates a route for tuning the

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-7179dc97-1c32-4db9-b171-d8792f099c1b.jpg)

combination of anisotropies by mixing E7 at different weight percentages, which are shown in square brackets, with CCN47. (m) Change in the specimen conductivity obtained by changing the concentration of E7 in CCN47.

Figure 3. Time-dependent dynamics versus electric-field parameters. (a, b) Director distortion of the kinetic solitons at 8.3 V and 8.9 V. Schematics of the alternating distorted convex and distorted concave director fields are presented with magnified patterns of them. (c-h) Typical time-dependent trajectories of solitons during different motions at various voltages at 20 Hz, stacked over 8 s. Scale bar = 50 µ m. (i) Angles of soliton motion with respect to 𝒏 0 at 20 Hz. (j) Velocity of the solitons at 20 Hz plotted as a function of the

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-6b9814c8-66ff-4321-ab79-5fad09a46496.jpg)

voltage. (k) Angles of soliton motion with respect to 𝒏 0 at 10 V Hz as a function of the frequency. (l) Velocity of solitons at 10 V as a function of the frequency.

Figure 4. Elastic collisional trajectory tracking. (a-c) Snapshots of collisions between solitons (with a tint-plate) and (d) reflection of solitons by an air boundary (with crossed polarisers) with schematics of solitonic topologies. In the second column from the left, the schematics of distorted biconvex director fields before collisions are drawn. Instantaneously, the director field alternates between distorted versions of two orientational state, as shown in Fig. 3a. Dashed lines in the schematics are used to clarify the positional offset, ∆. The overall trajectories are shown with colour-code traces and arrows in the right-most column. Scale bar, 10 µ m.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-2a06cb98-6b0f-4d98-be37-733cc04612df.jpg)

Figure 5. Fractal cloning of solitons. (a) Time-dependent trajectories of solitons during fractal cloning, stacked over (a) 10 s and (b) 5 s. Scale bars, 100 µ m. (c) Fractalisation occurring when the soliton size is nearly doubled.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-e4f96adc-b5af-43a0-8bac-c0aee5004946.jpg)

Figure 6. Pseudoparticle properties of solitons. (a) Pairwise soliton potential. (b) MSD of

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-564e0fd6-5fcf-4ea1-9a6f-b081552a57a4.jpg)

solitons as a function of time.

## Supplementary Information

Satoshi Aya  ,# and Fumito Araoka 

## Table of Contents

## Methods and Discussion

- 1. Ion localisation at soliton generation
- 2. Analysis of the soliton growth process
- 3. Packing of solitons
- 4. Structure of solitons
- 5. Conditions of soliton creation: comparison between the current solitonic systems and conventional convective systems
- 6. Kinetics of solitons
- 7. Go-and-stop motion of swimming solitons
- 8. Size of solitons
- 9. Feedback processes of the director and flow field through a flowalignment process
- 10. Discussion of soliton stabilisation
- 11. Radial distribution of solitons in the centered rectangular packing state
- 12. Relationship between the solitonic trajectories and the pair-wise potential curve
- 13. Characterisation of the electrical and dielectric properties

Supplementary Figures (Supplementary Figs S1-S12)

Supplementary References (Supplementary Figs S1-S9)

* To whom correspondence should be addressed.

E-mail: satoshi.aya@riken.jp (S.A.); fumito.araoka@riken.jp (F.A.)

#

Current contact address: South China University of Technology, E-mail: satoshiaya@scut.edu.cn

## 1. Ion localisation in solitons

Ions play important roles in stabilising solitons. Here, we show an essential experimental result regarding ion localisation in solitons, i.e., the observation of the redistribution of ions. In the experiment, we first switch on an electric field to create solitons. At an arbitrary time, we remove the electric field, and the solitons instantly disappear. Then, we wait for X seconds for the next application of the electric field to create solitons. If a localisation process of ions is involved in soliton generation, the diffusion or redistribution of ions should occur after the electric field is removed. We observe this by tracking the location at which the solitons are regenerated while tuning the waiting time X . As a result, we find that solitons are regenerated at the same location as before the electric field is removed if X is less than about 20 s (Fig. S1); otherwise, they are randomly regenerated in space. This means that some ions are indeed accumulated at the center of the solitons upon their generation. Once the diffusion of ions is complete after the electric field is removed, new solitons are randomly created in space at the second time of generation because there are new seeds for the solitons due to the localisation of ions. Taking a typical diffusion coefficient of the ions in liquid crystals of 2× 10 -9 m/s , it is calculated that the delocalisation of the ions takes about 20 2 s, which is consistent with the observed time scale (Fig. S2).

## 2. Analysis of the soliton growth process

A Kolmogorov-Johnson-Mehl-Avrami (KJMA) analysis of the structureforming events provides information on how homogenously these events unfold and how their nucleation and growth processes proceed S1-S4 . To estimate the occurrence frequencies of the electrically pumped solitons, we plotted the time evolution of the

soliton volume ratio during their creation using the Avrami expression, 𝑙𝑛[𝑙𝑛 1/{1 - 𝛷}] = 𝑚 𝑙𝑛 𝑡 + 𝑙𝑛 𝐾 , where 𝑋 , 𝑡 , 𝐾 , and 𝑚 are the soliton volume fraction, time, temperature-dependent Avrami coefficient, and Avrami exponent, respectively (Fig. S3). We found that the overall creation process presents a good fit with the Avrami scale with an Avrami exponent of 2.9. Considering the random incidence of each soliton within space, as revealed by polarising microscopy, the creation of solitons can be seen to be mediated by the homogeneous nucleation of soliton seeds, which is initiated by a localisation of ions, as discussed later, and the subsequent two-dimensional spreading of their network.

Now, we turn to consider why the KJMA model fits our data well. When creating solitons, switching on an electric field initiates the localisation of ions, and the localised ions becomes seeds that have potential to trigger the creation of the solitons. Let us assume the number of seeds at time 0 is N0 and these can be activated to create solitons through a firstorder-like reaction process. Then, the number inactivated seeds at time   can be written as t N0exp(-a ), and the activation rate of the seeds to be nucleus of the solitons per unit time t aN0exp(-a ). As a result, the emerging number of solitons is proportional to N0-aN0exp(-a ) t t ∝ 1exp(-k t n ), which is in a same expression to that of the KJMA model and validates the usage of the model to evaluate the growth process of the solitons.

## 3. Packing of solitons

We analysed images of packed solitons to gain a better understanding of their packing states. Figure S4a shows a swarm of solitons packed in a centered rectangular manner. A fast Fourier transform of Fig. S3a is shown in Fig. S3b, and the corresponding direct observation via the conoscopy of light diffraction from the packing

state of the solitons is shown in Fig. S3c. It is obvious that the fast Fourier and conoscopy results are consistent.

## 4. Structure of solitons

As stated in the main manuscript, the solitonic structure is dynamic. When we observe the structure of the solitons under a POM by a conventional low-speed camera, which has a capturing rate lower than the oscillating frequency of the dynamic solitons, the solitons appear as a time-averaged structure. In this case, they look like bearing topological defects near the their center. In order to clarify the time sequence of the structure, we perform high-frame-rate polarising microscopy with a 𝜆 /4 plate at 100 fps. In Fig. S5, the time evolutions of the transmittance in two domains of the observed textures are plotted as a function of time along with a schematic director field and the textures at specified times. T is the periodicity of the applied electric field. Clearly, it is shown that the transmittance of the biconvex structure is stronger than that of the biconcave one. This suggests that the maximum deviation angle of the director in the biconvex structure is larger than that in the biconcave one. Another key feature is that the inclination angle of the brushes in the textures between the biconvex and biconcave structures is very different, indicated by the white dotted lines. Lastly, both the switching times from the biconvex to the biconcave and from the biconcave to the biconvex are of the order of several millisecond, and the switching time from the biconvex to the biconcave is shorter. These suggest the structural oscillation is elasticity-driven. All of these evidences confirm that the biconvex and biconcave structures are not a simple symmetric inversion of the director field, but the director field is more distorted in the biconvex structure.

## 5. Conditions of soliton creation: comparison between the current solitonic systems and conventional convective systems.

The state diagram (Fig. 2a-e, Figs. S6,S7) infers an important connection between the current solitonic systems and conventional convective systems. We compare the known state diagram in a review by Eber et al. S5 (Fig. S6, Left) with our state diagram (Fig. S6, Right). It turns out that the G-, ||-, R- and ⊥ -states appear in similar ranges with the ones in conventional convective systems. The significant difference is that the Soliton state replaces the travelling waves and dielectric rolls regimes, and suppresses the chaotic regime.

## 6. Kinetics of solitons

In order to show the details of the kinetics of solitons, we show a state diagram in Fig. S7, focusing on the soliton state. The dynamic solitons just fluctuate around their equilibrium position at low frequencies (orange regime). At low frequencies and high voltages (Purple regime), the solitons swim parallel to the direction of the background director. At higher frequencies and high voltages (Blue regime), the solitons swim perpendicular to the direction of the background director. Since this regime is right next to the electrohydrodynamic regime ((V) ⊥ -state) where the direction of background flow is parallel to the direction of the background director, the swimming direction is dominantly determined by the flow as we stated in the Discussion part. In the narrow shaded regime between the blue and orange regimes, the solitons change their swimming direction between parallel and perpendicular to the direction of the background director, as shown in Figs. 3i-l. In this regime, the

background flow becomes slightly oblique, suggesting an additional travelling mode out of the xz plane becomes stronger.

## 7. Go-and-stop motion of swimming solitons

We observe the trajectory of the solitons upon their swimming, taken by a high-speed camera. This allows us to analyse the accurate dynamics of the solitons. Fig. S8a demonstrates an increment of x -coordinate of a soliton as a function of time when the soliton is moving along x direction. It is clear that x -coordinate of the solitons averagely-increases, and oscillates. Fig. S8b shows a power spectrum of Fig. S8a, revealing the rate of the go-and-stop motion is about 16.5 Hz, consistent with the frequency of the electric field, 16 Hz. This fact suggests that the go-and-stop motion is caused by the switching between biconvex and biconcave structures upon solitonic swimming.

## 8. Size of solitons

The size of solitons is determined by the balance between the elastic energy and the surface anchoring. As discussed in the manuscript, the solitons carry out-ofplane deformation along the   axis. Therefore, it is expected that the size of the solitons z is dependent on the sample thickness; that is, a thicker sample results in larger solitons. In Fig. S9, the linear dependence of the size of the solitons on the sample thickness is shown.

- 9.  Feedback  processes  of  the  director  and  flow  field  through  a  flow-alignment process

The orientation of the director tends to be parallel to the flow field as a result of the flow-alignment property. Since the flow field deviates from the initial state after the reorientation of the director, the orientations of both the flow field and director oscillate as a result of a self-regulating procedure (Supplementary Fig. S10).

## 10. Discussion of soliton stabilisation

Herein, we discuss why dynamic solitons have been so difficult to observe, despite the lengthy history of research on the electrohydrodynamics of liquid crystals, and, correspondingly, why our system is a special case. We again emphasise that the key property of the soliton system examined in this study is the frustration of dielectricity and conductivity. As shown in Fig. 2l, the Soliton-state appears only for (- +) mixtures at moderately low conductivities in the range of 8 × 10 -9 &lt; σ &lt; 4 × 10 -8 Ω -1 ∙ m -1 , which is two or three orders of magnitude lower than that of typical materials used in electrohydrodynamic studies but similar to previously reported values S6,S7 . The importance of this middle conductivity region can be understood by considering the coupling between the space charge Q and the electric field E . In the lower-conductivity region ( ∝ Q ), the number of free ions is spatially limited; as a result, the viscous torque proportional to the time-averaged quantity &lt; QE&gt; T S8 must destabilise the ground planar state, but is too small to overcome the total stabilising force produced by surface anchoring and dielectric alignment in the framework of the Carr-Helfrich electrohydrodynamic effect. In contrast, in the high-conductivity regime in which Q is high and uniformly distributed in space, a viscous torque dominates and produces a turbulent flow. In the intermediate region, there is a sufficient number of ions for the Carr-Helfrich electrohydrodynamic effect to apply, but these can still be spatially

localised. This consideration is well supported by the numerical calculation of electroconvective dynamics by Treiber and Kramer S9 . Overall, the value of &lt; QE&gt; T effectively reflects the appearance of spatially localised waves corresponding to the Soliton-state.

## 11. Radial distribution of solitons in the centered rectangular packing state

The radial distribution function of the solitons in the centered rectangular lattice, 𝑔(𝑟 𝑐𝑐 ) , was calculated using imaging analysis. The radial distribution function is calculated as the ratio of the average soliton density, ρ(𝑟 𝑐𝑐 ) , at a distance of (𝑟 𝑐𝑐 ) to the average density of solitons over the entire structure, ρ , as 𝑔(𝑟 𝑐𝑐 ) = 2ρ(𝑟 𝑐𝑐 )/ρ (Fig. S11).

## 12. Relationship between the solitonic trajectories and the pair-wise potential curve

In Fig. 4, it is shown that the solitonic dynamic behaviour upon a collision is changed depending on the positional offset, ∆, of the solitons. Especially, when the positional offset, ∆, of two solitons in a collision event is about zero, i.e., a head-on collision, the solitons repels and make a vertical repulsion by a displacement of about 15 µ m. This displacement is consistent with the value of the mean nearest-neighbor distance between solitons that is calculated by the pairwise potential curve (Fig. 6a). From this perspective, the collision is surely triggered by an elastic repulsion between the solitons.

## 13. Characterisation of the electrical and dielectric properties

The electric and dielectric responses of the liquid-crystal mixtures were characterised using dielectric spectroscopy (Solartron Analytical, impedance/gain-phase analyser 1260A with dielectric interface system 1296A). We focused on the lowfrequency range (0.01-10  Hz), in which the main response originates from the charge 4 polarisation of the mobile ions S10 . The mixtures were confined between parallel glass plates with thin layers of indium tin oxide and were subjected to a 0.05-Vrms AC electric field sweeping a frequency range of 0.01-10 4 Hz. Fig. S13 shows the dielectric spectra produced at various weight ratios of E7 to the host CCN47. Clearly, an increase in the ratio of E7 increases both the ionic relaxation frequency and dielectric strength as a result of the increased number of movable ions. Using the ionic characterisation method introduced by Sawada et al. S10 , we calculated the conductivities of the mixtures, as shown in Fig. 2m in the main manuscript.

## Supplementary Figures

Figure S1. Trajectories of solitons. At the time (indicated by orange) where the first application of the electric field is removed, the solitons instantaneously disappear. Within 5 s (&lt; X ), the electric field is applied a second time, and the solitons appear at

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-1f664a70-f6bd-4dfa-b571-1cc704544c4f.jpg)

same locations at which they disappeared when first application of the electric field was stopped.

Figure S2. Simulation results for the diffusion of ions in space. At an initial state at   = 0 t s, the ions are assumed to be localised at 12.5 µ m upon the generation of solitons. The curves of the spatial distribution of the ions at 1, 2, 4, and 20 s are given in blue, black, red and green colours.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-323a0dbb-0aef-491b-8d84-6312a2573ef0.jpg)

Figure S3. Avrami analysis of soliton creation with the time-dependent soliton volume ratio plotted on the Avrami scale.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-8570e158-7dbb-4064-8fdd-5e142d53bf28.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-f8f18a27-0954-4f48-aa0e-f8da16bc99c3.jpg)

b

Figure S4. (a) Centered rectangularly packed solitons. (b) Fast Fourier transform of (a).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-6ea3937c-a8ee-47da-a99a-9fb04c6e7fed.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-981a5b6c-adf6-4192-9d96-dce132a5d566.jpg)

(c) Conoscopy image of (a).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-5f49e5dd-171a-4c4f-9065-4ee889ad32ea.jpg)

Time

Figure S5. Time evolution of the structural variation in the solitons observed with highframe-rate polarising microscopy with a 𝜆 /4 plate.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-40d19844-ee35-48bd-8649-7f57ed6ded06.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-b4d0095f-485e-476b-8ff2-d68c08cecaa3.jpg)

2

Figure S6. Comparison between the state diagrams of the current soliton system (Left) and conventional convective systems (Right, Reproduced from É ber et al. s5 ).

Figure S7. State diagram showing different regimes in the soliton state.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-6df85ccb-3349-45d7-ab3e-948c51e65fb6.jpg)

Figure S8. Analysis of the go-and-stop motion of the solitons at an electric field of 16 Hz.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-3daa4156-0d2b-4581-8f40-7ca4a98ff098.jpg)

Figure S9. Thickness dependence of the diameter of the solitons.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-01c15dc4-f7e2-4747-a2c9-0af56096f99f.jpg)

Figure S10. (a) Roll pattern in relation to the wave vector and the direction of director.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-1118870a-d53d-4858-9e5e-3eca0db16784.jpg)

- (b) Schematic of the feedback mechanism of oscillation of the director and flow field.

Figure S11. Radial distribution function of centered rectangularly packed solitons.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-1e57af2d-280c-49ca-8480-73e86735c03a.jpg)

Figure S12. Thickness dependence of the effective mass of the solitons.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-a41b1195-da1b-4a7a-b9d8-a59a53092062.jpg)

Figure S13. Dielectric spectra of CCN47-E7 mixtures. Top: Real part of the complex dielectric constant as a function of the frequency. Bottom: Imaginary part of the complex dielectric constant as a function of the frequency.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00316v1-figure-8590f9b5-7229-4a18-ba33-33e5365dc4e9.jpg)

## Supplementary References

[S1] Avrami, M. Kinetics of phase change. I. General theory. J. Chem. Phys. 7 , 11031112 (1939).

[S2] Avrami, M. Kinetics of phase change. II. Transformation-time relations for random distribution of nuclei. J. Chem. Phys. 8 , 212-224 (1940).

[S3] Avrami, M. Kinetics of phase change III. Granulation, phase change, and microstructure. J. Chem. Phys. 9 , 177-184 (1941).

[S4] Torrens-Serra, J., Venkataraman, S., Stoica, M., Kuehn, U., Roth, S. &amp; Eckert, J.

Non-isothermal kinetic analysis of the crystallization of metallic glasses using the master curve method. Materials 4 , 2231-2243 (2014).

[S5] É ber, N., Salamon, P. &amp; Buka, Á. Electrically induced patterns in nematics and how to avoid them. Liq. Cry. Rev. 4, 101-134 (2016).

[S6] Brand, H. R., Fradin, C., Finn, P. L., Pesch, W. &amp; Cladis, P. E. Electroconvection in nematic liquid crystals: Comparison between experimental results and the hydrodynamic model. Phys. Lett. A. 235 , 508-514 (1997).

[S7] Kramer, L. &amp; Pesch, W. Convection instabilities in nematic liquid crystals. Ann. Rev. Fluid. Mech. 27 , 515-539 (1995).

[S8] Treiber, M. &amp; Kramer, L. Coupled complex Ginzburg-Landau equations for the weak electrolyte model of electroconvection. Phys. Rev. E 58 , 1973-1982 (1998). [S9] Sawada, A., Nakazono, Y., Tarumi, K. &amp; Naemura, S. Complex dielectric constant of liquid crystal materials containing ionic impurities in low frequency region. Mol. Cryst. Liq. Cryst. 318 , 225-242 (1998); Sawada, A., Nakazono, Y., Tarumi, K. &amp; Naemura, S. Space charge polarization by ions in a liquid crystal material. Mol. Cryst. Liq. Cryst. 331 , 457-464 (1999); Sawada, A., Tarumi, K. &amp; Naemura, S. Novel

Characterization method of ions in liquid crystal materials by complex dielectric constant measurements. Jpn. J. Appl. Phys. 38 , 1423-1427 (1999).
Satoshi Aya Fumito Araoka
Kinetics of motile solitons in fluid nematics
Solitary waves, dubbed "solitons", are special types of waves that propagate
for an infinite distance under ideal conditions. These waves are ubiquitously
found in nature such as typhoon or neuron signals. Yet, their artificial
generation and the control of their propagation remain outstanding challenges
in materials science owing to an insufficient understanding of the experimental
conditions and theoretical aspects. Herein, a generic strategy for forming
particle-like solitons and controlling their kinetics in nematic fluid media is
reported. The key to the realisation of the generation of solitons and the
control of their kinetics is the coupling between the fluid elasticity and the
background flow flux, as evidenced by experimental observations and theoretical
approaches. The findings of this study enable the exploration of solitons in a
wide range of materials and have technological ramifications for the lossless
transport of energy or structures.