## Topology of three-dimensional active nematic turbulence confined to droplets

Simon Čopar, 1, ∗ Jure Aplinc, 1, ∗ Žiga Kos, 1, ∗ Slobodan Žumer, 1, 2 and Miha Ravnik 1, 2

1 Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, 1000 Ljubljana, Slovenia

2 J. Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia

(Dated: August 5, 2019)

## Abstract

Active nematics contain topological defects which under sufficient activity move, create and annihilate in a chaotic quasi-steady state, called active turbulence. However, understanding active defects under confinement is an open challenge, especially in three-dimensions. Here, we demonstrate the topology of three-dimensional active nematic turbulence under the spherical confinement, using numerical modelling. In such spherical droplets, we show the three-dimensional structure of the topological defects, which due to closed confinement emerge in the form of closed loops or surface-to-surface spanning line segments. In the turbulent regime, the defects are shown to be strongly spatially and time varying, with ongoing transformations between positive winding, negative winding and twisted profiles, and with defect loops of zero and non-zero topological charge. The timeline of the active turbulence is characterised by four types of bulk topology-linked events - breakup, annihilation, coalescence and cross-over of the defects - which we discuss could be used for the analysis of the active turbulence in different three-dimensional geometries. The turbulent regime is separated by a first order structural transition from a low activity regime of a steady-state vortex structure and an offset single point defect. We also demonstrate coupling of surface and bulk topological defect dynamics by changing from strong perpendicular to inplane surface alignment. More generally, this work is aimed to provide insight into three-dimensional active turbulence, distinctly from the perspective of the topology of the emergent three-dimensional topological defects.

∗ These authors contributed equally to this work.

## I. INTRODUCTION

The ability of the material to employ externally or internally stored energy to spontaneously organize, flow, move, or change shape - i.e. be active - is found in a variety of materials [32]. A major class of active materials are active nematics [9], which exhibit strong collective behaviour and self organisation, emergent as local orientational order of material building blocks in material systems such as bacteria [4] and microtubule-kinesin mixtures [37, 45]. Active nematic formalism could be also applied to describe other biological systems, such as cytoskeleton [5, 43] and biofilm dynamics [24, 46]. A recurring behaviour in active nematics is the formation of topological defects, i.e. regions of broken orientational order, which are subject to topological invariant conservation [1]. Due to strong localized deformation of the orientational order, some active nematic defects can perform as strong effective sources of material flow, as dependent on the symmetry of the defects [7, 15]. A recurring property of topological defects in different active nematic systems is that at larger activities, they undergo into a regime of irregular and chaotic motion at low Reynolds numbers - called active turbulence. In the turbulent regime, generally, the topological defects are constantly created and annihilated [4, 45] which can be characterised by using statistical tools from classical turbulence [14, 53]. The threshold between regular and irregular dynamics is strongly dependent on the confinement [10, 21], higher order force multipoles [31] and friction, which can even lead to stabilization of ordered defect phases in active nematics [7, 44].

Geometrical confinements is today seen as one of the prime mechanisms for controlling the dynamics of active nematics [29, 50, 58, 59], notably also in the context of possible energy harvesting [49]. Active nematics interact with their confinement through hydrodynamic boundary conditions, such as slip, no-slip and friction, as well as through boundary conditions on the surface imposed orientational ordering, also called surface anchoring. Active nematic droplets based on the mixture of microtubules and molecular motors were shown to periodically morph their shape [25]. Active nematic layers organised at the surface of spherical droplets are also explored from the perspective of topological defect trajectories and their mutual interactions [21, 25, 26, 60]. When interfacing an ordered passive nematic fluid, active nematics can also exhibit directional streaming along the passive nematic director [17], and actuate the defects in the passive nematic host [18]. Fluid flow inside a droplet

can produce net propulsion of the droplet, which itself acts as an individual active swimmer particle [13, 30, 51].

Realising and controlling three-dimensional active nematic materials is an open emergent challenge in experiments [2, 8, 39, 40, 47, 61] and numerical simulations [19, 48, 53, 57], where one of the standing major complexities is how to characterise the complex 3D spatially and time-varying flow and orientational fields, that moreover have embedded topological defect regions [53]. Analogous systems exist that are dominated by topological defects in three-dimensions, distinctly passive (non-active) liquid crystal colloids [33, 52] and chiral ferromagnets [12, 41], where structures are controlled and characterised by using topological approaches. Indeed, due to the apolar nature of nematic director, the nematic (passive or active) in three dimensions can form not only point defects (as in 2D), but also defect lines and loops, which can be characterised with different topological invariants, including winding number, topological charge and self linking number [1, 62]. In passive nematics, coarsening structures during relaxation of a quenched state from isotropic to nematic also show notable analogy in view of defect loop topology to 3D active turbulent defect profiles [38]. To generalize, the idea of this paper is to apply concepts of topology to characterize structural events in 3D active nematic turbulence.

In this paper, we explore the topological defects regimes of three-dimensional active nematic under the spherical confinement and no-slip surface, specifically focusing on the topology-affecting events in the active turbulence. We show that for moderate activities, the active nematic assumes an oriented state with an off-center shifted single point defect in the form of a small loop, and vortex flow with direction-set angular momentum. Upon increasing activity, we observe the onset of topological turbulence, characterised by a firstorder structural transition and hysteresis between the offset point defect and the turbulent regime, also showing the corresponding phase diagram as dependent on activity and droplet size. In the turbulent regime, we show that the active turbulent dynamics can be interpreted as a time-series of topological defect affecting events: defect splitting, annihilation, merging and crossovers. Depending on surface anchoring conditions, the defects are in the form of closed loops (perpendicular anchoring) or surface-to-surface lines that effectively wet the confining surface (inplane anchoring), resulting in interesting surface-to-bulk conditioned turbulent dynamics. More generally, the results demonstrate dynamics of strongly confined active nematic and point towards using concepts of topology for controlling the

three-dimensional active nematic turbulence.

## II. RESULTS

The defect phenomena of the active nematic is explored in an elementary confinement of a droplet with fixed spherical shape by using mescoscopic numerical modelling of active nematodynamics [9, 15, 18, 20], which was shown to give good agreement with experiments on dense active nematic systems [9, 14, 56]. The approach relies on the dynamic coupling between the material flow and the mesocospic order parameter tensor Q ij that covers the orientational ordering of active nematic. The activity is described by dipolar-like forcing via the active stress contribution. The confining spherical surface is set to impose strong perpendicular (homeotropic; in Figs. 1-3) or inplane (degenerate planar; in Fig. 4) alignment of the active nematic at the surface, which experimentally, would correspond to different surface functionalities [34]. Using strong anchoring and fixed spherical shape allows us to discern the effects of shape from the effects of topology that are otherwise inherently intertwined. Generally, such fixed shape regime corresponds to having materials with large surface tension or a background medium that is rigid enough to support the shape of the active nematic droplet (e.g. gel-like [11]). We use no-slip boundary condition at the surfaces to simulate the host medium that resists flow. All distances are measured in units of the nematic correlation length ξ N , time is measured in units of the intrinsic nematic time scale τ N , and activity ζ in the units of L/ξ 2 N where L is the single elastic constant of the material. Equally, activity can be also described with the active length ξ ζ = √ L/ζ . Note that the nematic correlation length measures the effective size (thickness) of the defects and is given as a relative strength of the nematic elasticity vs. variations in nematic order, which are beside the activity and confinement - the key energetic mechanisms that affect the formation, structure and dynamics of topological defects. For more on the approach, please see Methods. Finally, note that this work focuses on the defect phenomena in three-dimensional active nematic, i.e. active nematic in the whole bulk of the spherical droplet, which is different to many current works which consider a thin layer of active nematic material, for example at the surface of a droplet [18, 60]

## A. Active regimes in spherical droplet

Depending on the activity of the material, two active regimes with distinct behaviour of topological defects are found for the three-dimensional bulk active nematic in the confinement of a spherical droplet with perpendicular alignment: for low activity, structure with an offset stationary point defect is observed (Fig. 1a) whereas for higher activities, a regime of three-dimensional active turbulence with spatially- and time-varying defect loops is observed (Fig. 1c). These two regimes are separated by a (hysteretic) structural transition (Fig. 1b). In the limit of no activity (i.e. also no material flow), the considered spherical confinement would exhibit a radial nematic director structure with a single +1 radial hedgehog point defect at the center, as imposed by the confining surface conditions; indeed, such profile is observed for activities below the activity ζ onset = 1 6 . × 10 -3 L/ξ 2 n (see also Fig. 2). Note that the hedgehog point defect appears in the form of a small ring of local winding number +1 2 / , with diameter of the order of nematic correlation length, and is topologically equivalent to a point [55]; the exact small-scale structure will in practice depend on the microscopic building blocks of the material. Therefore, in view of our work, loops with radius of the order of the nematic correlation length will be called point defects. Above the activity ζ onset a self-sustained structure of a twisted director profile and a vortex flow appears, characterised by the orientation and magnitude of the angular momentum vector (Fig. 1a). Similar as in (passive) cholesteric nematic droplets [42], twisted director field pushes the point defect away from the droplet center. Since our system has no preferred chirality, the shift of the point defect can be both along or against the angular momentum vector. Additionally, there is a velocity field component of roughly ∼ 2 5 . -times smaller magnitude then the circular vortex flow which goes along the center of the droplet towards the point defect and then near the surface back to the opposite end of the spherical droplet confinement. For increasing activity, the magnitude of the vortex flow increases and the defect shifts more and more toward the surface of the droplet, and as a critical activity is reached, the small defect ring opens into a larger loop (see Fig. 1b). At this value of the activity, the activity-induced flow can overcome the elastic barrier of the nematic field to transform the point defect into a topologically equivalent large defect loop. This process involves contributions of defect core line tension, nematic elasticity, advection, backflow, shear active flow as well as coupling to confining geometry and actual local and global defect topology. Above the critical activity,

FIG. 1. Active nematic regimes in spherical droplet confinement with perpendicular

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00904v1-figure-c001d886-6870-44e0-b1c2-9a9c6b95c677.jpg)

surface alignment. Snapshots of dynamics for different activities, showing time evolution, velocity streamlines and director field. (a) Active nematic regime at low activities with single point defect (in red, in the form of a small ring) displaced from the droplet center and flow vortex (gray streamlines) with the angular momentum vector Γ . Note that the flow velocity component in the direction of -Γ is ∼ 2 5 . -times smaller than the flow inplane component in the xy plane. Right panel shows director profiles in the plane of the vortex and perpendicular to the vortex. Red color corresponds to the director along x axis, green color along y axis and blue color along Γ . Bottom panel shows how with increasing activity the point defect (in red) moves further away from the droplet center and flow magnitude (blue arrows) increases. (b) Transition of point defect regime into the turbulent regime upon an increase of activity from 0 035 . L/ξ 2 N to 0 036 . L/ξ 2 N ; note stretching and deformation of the defect loop, and the disintegration of the regular velocity profile. (c) Turbulent active nematic regime. Defect loops are shown as isosurfaces of reduced degree of order (in red) and flow field with gray streamlines. Right panel shows the corresponding director in the given plane; the local variability of the defect cross section is illustrated with splay-bend parameter (yellow and blue isosurfaces; see Methods). (d) Selected time-line of active turbulence in the active nematic droplet. Red isosurface of the scalar order parameter is drawn at S = 0 45 . ; blue arrows and bright green isosurfaces represent velocity field above selected cut-off magnitudes (see Methods).

FIG. 2. Analysis of active nematic regimes in confinement of spherical droplet. (a)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00904v1-figure-aae63b99-35a0-47f4-a2f8-8e0df493bf79.jpg)

Multiple material characteristics of the three-dimensional active nematic in spherical droplet confinement of radius R = 136 ξ N with perpendicular surface alignment as dependent on the activity. The calculations are performed for increasing or decreasing activity, clearly showing a pronounced hysteresis. At each activity change, previous structure is taken as the initial condition and simulation is run for 23000 τ N for each data point. Top panel shows the total length of defect loop cores at a given activity, with the bars indicating the standard deviation of the natural oscillation of the defect length over time. Second panel in (a) shows average radial distance of the defect core regions from the droplet center. Third panel gives the total angular momentum set by the flow and its standard deviation in time. The bottom panel gives the probability for finding given number of defect loops within the droplet. Note three active regimes - point defect, hysteretic regime and developed active turbulence - each with different response of the active material characteristics. Linear fit of the angular moment dependence on the activity, in the point defect regime, is used to set the onset activity below which the vortex flow profile has negligible magnitude and the point defect remains at the droplet center. Graphs of angular momentum and probability are shown only for hysteresis behaviour at decreasing activity. Hysteretic region is delimited by activity numbers Ac ≈ 22 1 . and Ac ≈ 25 5 . , for the considered droplet size. (b) Phase diagram of active regimes in spherical droplets with perpendicular surface alignment, as given by droplet radius R and activity

the active nematic transitions into a chaotic irregular behaviour - i.e. the active turbulence regime. There can be strong material flow due to large elastic deformation around the core of the loop, and overall, strong dynamic deformations of the loop are observed, including topology-modifying events, such as breakup of a loop into multiple loops, as shown in Fig. 1c,d.

Performing a more quantitative study of the active nematic regimes in the spherical droplet confinement, several distinct material characteristics are analysed (as shown in Fig. 2), that can assess the three-dimensional nematic orientational ordering and material flow, as well as the topology of the defects. For each simulation, we changed the activity between two neighbouring data points in Fig. 2a, let the system evolve for 4500 τ N to reach a dynamic steady state, and then over a time interval of 18500 τ N collected data on total defect length, (average) defect position, total angular momentum, and probability for observing N-number of defect loops (Fig. 2a). Note that, in the analysis we employ direct 3D tracking of defect line cores based on variations in the nematic degree of order, as further explained in Methods. The defect length inside the droplets quantifies that at low activities there is a regime of a single point defect, but upon increasing activity, a structural transition occurs to the active turbulent regime, in which the single defect loop or multiple loops dynamically transform and vary in-time, with total defect line length changing in time at the order of multiple 10% from the average value. The position of the centre of the point defect changes with increasing activity and shifts from the droplet centre towards the surface of the droplet. As the active nematic transitions to the turbulent regime, the active defect loops explore and transverse the available volume irrespective of actual activity. The average radial position of defect segments from the droplet center settles in the turbulent regime to ≈ 0 3 . R , as at high activities the defect segments become more uniformly distributed. The average defect position also shows hysteretic behaviour. The average defect core position from the droplet center well measures the point defect position in the low activity state, but in the turbulent regime, it is more meaningful to observe the probability density of defect line segments per unit volume, with respect to relative distance from the center of the droplet, as shown in Supplementary Fig. 1. The angular momentum grows within the point defect regime continuously with activity, according to a linear dependence ( ζ -ζ onset ) . ζ onset is the activity for the onset of vortex flow, which we determine numerically by a linear fit of the angular momentum dependence on activity. Note that such emergence of active flows in confined

and defect systems is known to depend on the symmetry and profile of the equilibirum (noactivity) structure [16, 27, 54]. When the active nematic passes to the active turbulence regime, the average magnitude of the angular momentum gradually saturates with activity, which again signifies the break up of the regular defect dynamics. At this transition between point defect and active turbulence regimes, also the process of defect loop breakup (i.e. creation of new loops) and annihilation of emerging loops begins, resulting in a growing fraction of time when there is more than one loop present in the spherical droplet.

Activity is the main control parameter for the transition between the point defect regime and the active turbulence; however, this transition could be also driven for example by changing the active nematic elastic constant or the size of the confining cavity. Figure 2b presents a phase diagram of active regimes with respect to activity and droplet size, showing that larger droplets transition into active turbulent regime at lower activities, whereas smaller droplets require larger activites to undergo structural transition from regular to irregular dynamics. Indeed, such behaviour can be well understood by considering the dimensionless activity number Ac = √ ζR /L 2 which gives the ratio between the confinement size R and active length ξ ζ . In Fig. 2b, the transition from point defect to active turbulence is observed roughly at Ac ∼ 20 .

The structural transition between the point defect structure and the active turbulence exhibits a clear hysteresis, which is a signature of discontinuous (first-order) transition between the two structures (Fig. 2a). The ordered phase with the offset point defect can be 'overactivated' (i.e. with activity above the transition value), which we simulate by gradually increasing activity in small steps from the low activity state until the turbulence appeared, assuring for each step in the activity change that the system reaches the dynamic steady state. Alternatively, the turbulent state can be 'under-activated' (i.e. with activity below the transition value), before the turbulent defect loop collapses into a point defect, which we simulate by gradually decreasing activity in small steps. The turbulent state can only appear when the flow is strong enough to overcome the elastic tension force of the defect line. Once the defect line is extended, it can remain extended even at lower activities. In the turbulent regime the system behaves chaotically, where changes in the shape of defects can lead to a very different future evolution of the system. The system has no static equilibrium, but reaches a dynamical steady state in which macroscopic variables, such as the defect length, angular momentum and average number of loops settle to fluctuating around

an average value.

## B. Topology of defect loops in active turbulence

Topology is a natural tool for studying topological defects that goes beyond exact detailed material characteristics, but rather focuses on the overall structure, and we use topological analysis for the characterisation of active turbulent regime under the confinement of active nematic into spherical droplets. We focus on topological properties of the actual defect loops - their number, how they are connected to each other, and their topological charge - as well as on the local variation of the nematic director surrounding the defect loops, i.e. considering defect loops not only as lines but rather as ribbon-like objects. Figure 3 shows the dynamics of active nematic turbulence in the spherical droplet, analysed from the perspective of the topology of three-dimensional active defect loops. In Fig. 3a, a time series of defect topology-affecting events is presented - importantly, with time-axis in scale for a selected time interval of the active tubulence in an active nematic droplet ( R = 136 ξ N and ζ = 0 04 . L/ξ 2 N ). Four different topological events are observed to occur stochastically in time, with the average event rate set by (and increasing with) the activity. New loops of zero topological charge are observed to detach (pinch-off) from a preexisting loop (the existence of at least one loop is topologically enforced by the droplet perpendicular surface conditions), extra loops annihilate, or merge back into the loop with nonzero charge, which by itself cannot disappear. Crossovers of the defect loop with itself can rewire the loop, reversing the path orientation of one part relative to the other. More generally, the visualised map reflects the conservation laws that would formally be described with a Morse theory, where topological events are critical points (saddles, endpoints) of the manifold comprised of the defect loci over time (the world sheet) [36]. Examples of loop splitting, annihilation of a single loop, merger of two loops, and a crossover (i.e. self-intersection) are shown in snapshots in Fig. 3c-f and in more detail in Supplementary Fig. 2.

The director cross section of defect loops in the regime of active turbulence varies along the defect loop in a complex way, and also changes in time, as shown in Fig. 3b and Supplementary Fig. 3. The defect loops in the three-dimensional active turbulent regime exhibit sections of local director symmetries of +1 2 / and -1 / 2 winding number and twisted profile, which has profound effect also on the local dynamics and local active flow. It is well known

FIG. 3. Topological events in turbulent regime. (a) Graphically depicted topological flow diagram for the active nematic defects loops in a selected time interval. Note that the diagram time axis is in scale (black line) and is given in units of the intrinsic nematic time scale τ N . At all times, one of the loops carries topological charge +1 and other have zero charge. (b) A snapshot of active turbulence with three loops. The color-coding of the disks for the selected defect loop shows local winding number of the defect line, where blue refers to +1 2 / , red to -1 2 / winding number and purple to local twist defect line profile. The variation of the defect line cross section is also illustrated more qualitatively with the splay-bend parameter (yellow and blue isosurfaces). (c-f) Time-sequences of selected loop topology-changing events: splitting, annihilation, merging and crossover. Panel (c) shows splitting of a neutral loop (upper) from a loop with a +1 topological charge (lower). The neutral loop soon annihilates in panel (d). Droplet radius is R = 136 ξ N and activity ζ = 0 04 . L/ξ 2 N .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00904v1-figure-96847ae6-318c-4951-8776-b7d011bfb23f.jpg)

from 2D active nematic systems that +1 2 / defect profiles are generally strong generators of flow and strongly propel this type of (in 2D) point defects, whereas -1 2 / defect profiles do not generate flow and are less motile [15]. Indeed, it is this variability of the director

field that we now observe in the local cross-sections of the three-dimensional defect loops (see Supplementary Fig. 3), which then affects the dynamics of not only itself but all defect loops, through both nematic elasticity and hydrodynamics. The role of local cross-sections is known to affect the dynamics of reconnection events between defect line segments [22, 23]. In our system, immediately before and after topological reconnection events the director orientations of both defect line segments generally lie in the same plane, and rewire in a way that makes the director between the defect lines uniform (without introducing topological solitons). For example, two defect line segments with locally twisted profile may reconnect into another pair of twisted profiles, a pair of +1 2 / and -1 / 2 local profiles, or other intermediate profiles with a net neutral winding number (see Supplementary Fig. 2). Finally, the fact that the defect lines in spherical droplet confinement with perpendicular conditions are closed into loops imposes that this local director profile ( +1 2 / , -1 / 2 or twisted) must come back to the same orientation when encircling the loop, which encodes the total topological charge carried by the loop [62].

At the level of the topology of the loop as a whole with surrounding director, we observe that one of the loops always carries a +1 topological charge - set by the +1 topological charge imposed by the perpendicular (homeotropic) surfaces - whereas the rest are topologically neutral (i.e. charge zero) and can annihilate without contacting another loop (see Fig. 3d), or alternatively, also merge back with another loop. From purely topological perspective, also pairs of oppositely charged loops could be formed (not observed in our simulations), especially in larger systems and with larger activities, while still preserving the net charge conservation set by confinement. However, the formation of oppositely charged pairs is - at least in the considered material regime - energetically even more unfavourable as the formation of zero-charged defect loops. Finally, a possible further approach to classify the active defect loops topologically could involve quaternionic approach [62] or construction of Pontryagin-Thom surfaces [6].

## C. Surface-contacting topological defects in active turbulence of spherical droplets with inplane alignment

Active turbulence in active droplets of fixed spherical geometry with strong surface degenerate inplane alignment (i.e. planar degenerate anchoring) is explored, showing a system

FIG. 4. Surface and bulk active nematic topological defect transitions in spherical

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00904v1-figure-48eb1f9d-d31f-4236-a90e-e5f6f425747c.jpg)

droplets with planar degenerate surface alignment. (a) Snapshot of active nematic droplet showing two bulk defect lines, which terminate at the surface (ends of one are encircled with magenta), and a boojum defect (encircled in cyan). The inplane director field at the spherical droplet surface is represented by grayscale streaks. Bottom panel shows Peirce quincuncial projection map of the surface director and defects from the entire sphere onto a square; the sphere can be reproduced by folding four corners of the square back together into one point where the yellow square represents the equator of the sphere. Note the pointed out defects encircled in magenta and cyan, as shown in the sphere and the projection square. (b) Annihilation sequence of a selected ± 1 2 / defect pair at the surface of the active nematic droplet shown using Peirce projection map. During the sequence a surface profile with five +1 2 / defects and one -1 2 / defect is transformed into four +1 2 / defect profile, conserving the +2 surface topological charge set by the Euler characteristic of the sphere. In the bulk of the droplet, the sequence corresponds to a line segment beeing annihilated on the surface, as shown in (d). The annihilating surface defect pair is encircled in orange. (c) At the droplet surface, +1 surface boojum breaks into two +1 2 / defect profiles (encircled in green). In the bulk, this surface process corresponds to a defect line segment interacting with the surface boojum, resulting in two bulk defect lines, as shown in (e). (d) Full view of the changes in the bulk defect structure (in red, drawn as isosurfaces of S = 0 51 . ) that corresponds to the defect annihilation shown in (b). (e) Full view of the changes in the bulk defect structure (in red, drawn as isosurfaces of S = 0 523 . ) that corresponds to the breaking of a boojum shown in (c)

of interesting surface-to-bulk coupled defect dynamics. Under degenerate planar anchoring, the defect lines within the bulk are not any more necessarily closed into loops, but can also terminate at the surface, as shown in Fig. 4. Indeed, in addition to closed defect loops observed for droplets with perpendicular surface alignment, defect lines with surface-to-surface spanning ends and surface boojum defects are observed, resulting in the turbulent defect dynamics which in the bulk is generally similar to the dynamics in homeotropic droplets but it is additionally coupled to the topological events at the droplet surface. Figure 4 shows two such surface events using the Peirce quincuncial projection of the surface director and defects from the sphere onto a square: in Figs. 4b,d, a pair of ± 1 2 / defects within the surface director field is annihilated, and in Figs. 4c,e, a surface boojum is split into a pair of +1 2 / defects. During the pair surface annihilation event, a structure with four +1 2 / defects and a ± 1 / 2 defect pair is transformed into a structure with only four +1 2 / defects. In the bulk, such event: (i) either corresponds to the annihilation of a defect surface-to-surface line segment on the surface of the droplet (which was the actual case in the presented simulation, as shown in Fig. 4d), or (ii) a merger of two ends of defect line segments into one segment and the following detachment of the defect segment from the surface. In Fig. 4c, a +1 surface boojum interacts (and merges) with a section of a bulk defect, which results in splitting of the boojum into two +1 2 / defects, transforming the surface director profile with two +1 2 / defects and the +1 boojum into a surface field with four +1 2 / defects. In the bulk of the droplet, the defect line that merged with the boojum splits into two parts, each of them terminating at one of the newly created +1 2 / defects at the surface, as shown in Fig. 4e. Additionally, not shown here, we observe creation of surface defect pairs when an existing bulk defect loop segment is pushed from the bulk to the surface, where it is split into two ends (seen at the surface as a pair of ± 1 / 2 defects), or also alternatively, a reconfiguration of the surface charges by merging of combinations of -1 2 / defects and +1 boojums.

Topologically, the observed active turbulent defect dynamics in droplets with inplane alignment is determined by the director on the surface where the total winding number (2D topological charge) must be equal to the Euler characteristic of the confining surface χ . Specifically, in our droplets, the net 2D surface topological surface charge must equal χ = 2 , which in the limit of zero activity, is realized by a pair of +1 boojums at the opposite poles of the droplet. However, for non-zero activity, as shown if Fig. 4, the surface boojums are only rarely observed and the surface charge is rather distributed between four +1 2 / defects

and additional pairs of ± 1 / 2 defects, which are actually endings of the bulk +1 2 / or -1 2 / defect lines touching the droplet surface. More generally, the topological dynamics of the defect lines (and loops) in the 3D bulk of active nematic induces a lower dimensional active defect dynamics at the confining surfaces.

The turbulent defect dynamics at the surface of the droplet shares similarities with the 2D active nematic turbulence (e.g. in films or shells); however, it is different from the point that it is actually coupled to the active flow and defect reconfiguration of the entire bulk. The no-slip boundary condition for the material flow, used in this work, effectively prevents surface defects to be driven by the activity directly on the droplet surface, but it is the active dynamics of bulk defects that cause apparent active motion, creation and annihilation of surface defects in the turbulent regime as provided via effective elastic line tension and splitting of bulk defect lines at the surface. In different realizations, though, we very much envisage that the velocity boundary condition could be different, from no slip, to partial slip or (full) slip, in which case one would expect similar defect dynamics, but now also propelled through the surface. Alternatively, the defect dynamics might be also affected by different surface and bulk rotational viscosities, effectively, causing the defects in the bulk and at the surface to move at different time scales. Finally, such combined actively-driven surface and bulk dynamics might lead to the emergence of effectively coupled or decoupled dynamics of defects at the surface and in the bulk.

## III. DISCUSSION

Confinement is today seen as one the major routes for controlling active matter - including active nematics - with multiple geometries explored [25, 29, 58, 59]. From the perspective of this work, we should emphasize that we purposely chose the regime where the droplet radius R is generally comparable but larger than the active length ξ ζ = √ L/ζ , i.e. R/ξ ζ &lt; ∼ 35 , and the surface interactions (surface anchoring) are strong, such that the confinement has a profound role. The main effect of confinement with the strong homeotropic boundary condition is that it enforces existence of at least one bulk topological defect (loop) by imposing non net-zero topological charge in the bulk of the droplet, and similarly strong planar anchoring imposes non-zero surface topological charge, thus enforcing existence of surface defects. This confinement-imposed existence of defects in the limit of low activity

has a notable effect on the onset of active turbulence as it circumvents the need for spontaneous creation of defect pairs. Instead of a spontaneous defect creation, existing defect or defects get transformed - e.g. existing defect splits into a defect pair. As a result of this lower energy barrier mechanism, in this work we do not observe events where new defects would emerge from an instability of the homogeneous nematic field and would involve spontaneous appearance of a single loop with zero topological charge, but novel defect loops are instead created by a breakup from an existing defect loop. In similar manner, one could envisage that some surface anchoring profiles or structures with sharp edges, that would cause strong local distortions of nematic could perform as nucleating regions for the emergence of the active nematic defects. Finally, the topology and geometry of the confining surfaces thus can play an important part in generating distinct active regimes.

The major body of current work in active nematics is for two-dimensional or quasi-twodimensional regimes and geometries [4, 45]. From the perspective of topological defects, the key difference between 2D and 3D active nematics is that the defects change from defect points in plane to defect lines and loops in three-dimensional space, which introduces a notable increase in the complexity of tracking and monitoring apporaches for the defects. Topologically, the key difference is that in 2D nematic fields only +1 2 / and -1 / 2 defects are possible, whereas in 3D nematics the defects can be in the from of points, loops or walls, and even lines if these can terminate appropriatelly like on surfaces. Also, solitonic solutions are starting to be seen in objects like active skyrmions [35]. These three-dimensional nematic defect structures can be characterised with different topological invariants (in 2D nematics there is only 2D topological charge) that account for the topology of the defects as a whole, such as 3D topological charge, linking and self-linking number. Defects can be characterized also locally through their cross-section, for example with the winding number. At the level of orientational order fields, active nematics are topologically equivalent to the passive (i.e. not active) nematics, in which 3D defect structures as complex as knots and links were demonstrated [33, 52]. It is an interesting question if similarly complex three-dimensional topological defect structures could emerge also in active nematics and moreover, if they could have some role in real living matter.

In conclusion, using numerical modelling, we demonstrate the topology of three-dimensional active turbulence in droplets with fixed spherical shape of active nematic, considering both droplets with perpendicular and inplane imposed surface orientation. Structural transi-

tion with developed hysteresis between bulk regular offset point defect regime and irregular regime of active turbulence is shown for homeotropic droplets. Defect length, defect position, angular momentum and number of defect loops are used as bulk material characteristics that can account for the coupled material flow, orientation and topological variability of the three-dimension active nematic in turbulent and non-turbulent regime. The bulk active defect lines and loops are shown to exhibit local profiles with spatially and in-time varying cross-section, with winding number of +1 2 / , -1 2 / and twisted profile. Individually, defect loops are generally topological-charge neutral, except for compensating the charge imposed by the confining surface. Topological flow diagram is shown as a strong tool for topological characterisation of active turbulence, identifying distinct set of topology-affecting or changing events - defect splitting, annihilation, merging and crossovers - with their conservation laws drawing on analogy with Morse theory. For droplets with degenerate inplane surface alignment, we show that active defects can touch (wet) the surface, creating active defect line segments that span between different surface regions, topologically coupling the surface and bulk defect dynamics. More generally, the aim of this work is to approach - in full three-dimensions - the complex temporal and spatial variability of active matter using a topological tool-box and moreover, to explain the observed continual dynamics in view of distinct transitions between selected well-characterized topological structures.

## ACKNOWLEDGMENTS

Authors acknowledge discussion on 3D active defect lines with D. Beller. The work was supported by the Slovenian Research Agency (ARRS) under contracts P1-0099, J1-9149 and L1-8135, and COST action EUTOPIA (CA17139). M.R. also acknowledges support under EPSRC grant EP/R014604/1 at Isaac Newton Institute, University of Cambridge (programme The mathematical design of new materials).

## AUTHOR CONTRIBUTIONS

S. Č., J. A. and Ž. K. contributed equally to the work. J. A. and Ž. K. performed numerical simulations and analysed the results. S. Č. performed geometric and topological analysis. Ž. K. developed numerical code. M. R. and S. Ž. supervised the work. All authors

contributed to writing of the manuscript. M.R. initiated and led the research.

## APPENDIX: METHODS

## 1. Modelling of active nematics

We apply mesoscopic continuum description of (dense) active nematic, that is based on coupled dynamic equations for orientational order and the material flow field [9, 18, 20]. Roughly, this approach is based on the material flow generation via the active forces caused by the distortions in the orientational ordering of the active nematic, and the back-coupled response of the active nematic orientation to the generated flow. The nematic ordering is described by a traceless tensor order parameter Q ij whose largest eigenvalue is the degree of order S (also called scalar order parameter) and the corresponding eigenvector gives the main ordering axis, called the director n . The dynamics of the Q-tensor is given by the adapted Beris-Edwards equation [3]

$$\left ( \partial _ { t } + u _ { k } \partial _ { k } \right ) Q _ { i j } - S _ { i j } = \Gamma ^ { \text{bulk} } H _ { i j },$$

where ∂ t is a derivative over time t , ∂ k is a derivative over k -th spatial coordinate, u is fluid velocity, and Γ bulk the rotational viscosity coefficient. Molecular field H ij drives the system towards the equilibrium of Q ij and can be interpreted from the free energy F (but also alternatively, see Ref. [9]):

$$H _ { i j } = - \frac { \delta F } { \delta Q _ { i j } } + \frac { \delta _ { i j } } { 3 } \text{Tr} \frac { \delta F } { \delta Q _ { k l } },$$

where F is written in the Landau-de Gennes form as:

$$F & = \int _ { V } \left \{ \frac { A } { 2 } Q _ { i j } Q _ { j i } + \frac { B } { 3 } Q _ { i j } Q _ { j k } Q _ { k i } + \frac { C } { 4 } ( Q _ { i j } Q _ { j i } ) ^ { 2 } \\ & \quad + \frac { L } { 2 } \frac { \partial Q _ { i j } } { \partial x _ { k } } \frac { \partial Q _ { i j } } { \partial x _ { k } } \right \} \text{d} V.$$

Derivatives of Q ij in Eq. 3 describe the effective elastic behavior of the director field, where L is the elastic constant. Summation over the repeated indices is implied. A B , , and C are parameters that can be used to tune the nematic phase behavior.

The advection term S ij couples the velocity and nematic ordering:

$$S _ { i j } & = ( \chi D _ { i k } - \Omega _ { i k } ) \left ( Q _ { k j } + \frac { \delta _ { k j } } { 3 } \right ) \\ & \quad + \left ( Q _ { i k } + \frac { \delta _ { i k } } { 3 } \right ) ( \chi D _ { k j } + \Omega _ { k j } ) \\ & \quad - 2 \chi \left ( Q _ { i j } + \frac { \delta _ { i j } } { 3 } \right ) Q _ { k l } W _ { l k },$$

where D ij is the symmetric, and Ω ij the antisymmetric part of the velocity gradient tensor W ij = ∂ u i j . The alignment parameter χ depends on the molecular shape and defines the flow-aligning or flow-tumbling regime.

Surface alignment is modelled as follows. For droplets with perpendicular surface alignment (homeotropic anchoring), strong anchoring with fixed radial director at the droplet surface is assumed. Droplets with inplane degenerate planar alignment (degenerate planar anchoring) are modelled by the surface free energy density:

$$f _ { \deg } = W _ { \deg } \left ( \tilde { Q } _ { i j } - \tilde { Q } _ { i j } ^ { \perp } \right ) ^ { 2 },$$

where ˜ Q ij = Q ij + Sδ ij / 2 , ˜ Q ⊥ ij = ( δ ik -ν ν i k ) ˜ Q kl ( δ kj -ν ν k j ) , and ν is the surface normal. Surface Q-tensor field follows the dynamics

$$\dot { Q } _ { i j } = \Gamma ^ { \text{surf} } \left ( H _ { i j } ^ { \text{surf} } - \frac { \delta _ { i j } } { 3 } \text{Try} H _ { k l } ^ { \text{surf} } \right ),$$

where

$$H _ { i j } ^ { \text{surf} } = - \frac { \partial f _ { \text{vol} } } { \partial ( \partial _ { k } Q _ { i j } ) } \nu _ { k } - \frac { \partial f _ { \text{deg} } } { \partial Q _ { i j } }$$

and f vol is the free energy density expressed in Eq. 3.

The fluid velocity obeys the incompressibility condition and the Navier-Stokes equation

$$\nabla \cdot u = 0,$$

$$\rho \left ( \partial _ { t } + u _ { k } \partial _ { k } \right ) u _ { i } = \partial _ { j } \Pi _ { i j },$$

where ρ is the fluid density and Π ij the stress tensor, which consists of a passive and an

active part Π ij = Π passive ij +Π active ij , where

$$; \Pi _ { i j } = \Pi _ { i j } ^ { \dots } \cdots + \Pi _ { i j } ^ { \dots } \cdots, \text{ where} \\ \Pi _ { i j } ^ { \text{passive} } & = - P \delta _ { i j } + 2 \chi \left ( Q _ { i j } + \frac { \delta _ { i j } } { 3 } \right ) Q _ { k l } H _ { k l } \\ & \quad - \chi H _ { i k } \left ( Q _ { k j } + \frac { \delta _ { k j } } { 3 } \right ) & & ( 1 0 ) \\ & \quad - \chi \left ( Q _ { i k } + \frac { \delta _ { i k } } { 3 } \right ) H _ { k j } - \partial _ { i } Q _ { k l } \frac { \delta F } { \delta \partial _ { j } Q _ { k l } } \\ & \quad + Q _ { i k } H _ { k j } - H _ { i k } Q _ { k j } + 2 \eta D _ { i j }, \\ \Pi _ { i j } ^ { \text{active} } & = - \zeta Q _ { i j }. & & ( 1 1 ) \\ \end{cases}$$

P is fluid pressure, η the isotropic viscosity contribution, and ζ is the activity parameter characterizing the strength of force dipoles of contractile ( ζ &lt; 0 ) and extensile ( ζ &gt; 0 ) objects.

Coupled equations for the fluid velocity u i and the nematic order Q ij are solved numerically by the hybrid lattice-Boltzmann algorithm [10, 18]. The hybrid algorithm consists of an explicit finite difference method for the Q-tensor evolution (Eq. 1), and the D3Q19 lattice Boltzmann model for the Navier-Stokes equation and the compressibility condition (Eqs. 8, 9). The nematic stress tensor is implemented in the lattice Boltzmann algorithm as a force contribution with the half-force correction [28]. The spherical droplet cavity is allocated on a rectangular grid and a no-slip boundary condition is implemented by the bounce-back rule. A radial surface normal is allocated on the surface nodes, and it is used for the calculation of the surface Q-tensor field either for homeotropic or for planar degenerate anchoring. We performed multiple test to verify that results are not affected by spurious velocities and numerical method induced symmetry or symmetry breaking. For example, in Fig. 1a, for different random initial conditions, the defect is displaced from the center in different (arbitrary) directions that are typically not along the (rectangular) simulation grid axes. The results of the simulations are expressed in the units of elastic constant L , nematic correlation length ξ N = √ L/ A ( + BS eq + 9 2 CS 2 eq ) where S eq is the equilibrium nematic degree of order, and nematic intrinsic time scale τ N = ξ 2 N / Γ bulk L . The phase parameters are set to A = -0 190 . L/ξ 2 N , B = -2 34 . L/ξ 2 N and C = 1 91 . L/ξ 2 N , the nematic is in alignment regime with χ = 1 , isotropic viscosity contribution equals η = 1 38 . ξ 2 N /Lτ N , and the strength of the planar degenerate anchoring is W deg = 6 6 . · 10 -4 L/ξ N with surface rotational viscosity parameter Γ surf = 0 67 Γ . bulk /ξ N . Grid resolution is set to ∆ x = 1 5 . ξ N and time resolution to

∆ = 0 057 t . τ N . Velocity field of the active nematic can be analysed also by calculating the total angular momentum Γ = ∫ ρ r × u d V , where r is the distance from the droplet centre and integration is performed over the whole volume of the active nematic droplet.

The parameters of the considered 3D active nematic system can be summarised by introducing two selected dimensionless numbers (according to Eqs. 1 and 9): Ericksen number Er = uR/ Γ bulk L comparing the viscous terms to the elastic terms and activity number Ac = √ ζR /L 2 characterizing the relative strength of activity vs. confinement, where u is maximum velocity in droplets of radius R and activity ζ . In our simulations we take activities of up to ζ &lt; ∼ 0 06 . L/ξ 2 N and droplet radius of R = 136 ξ N which corresponds to Er &lt; ∼ 100 and Ac &lt; ∼ 35 .

## 2. Tracing and visualisation of defects

Modelling of active nematics gives a continuous three-dimensional Q-tensor field and a velocity field on a discrete grid at each time step. To perform the topological analysis, geometry of defect lines is extracted from the Q field. Defects are detected by first finding the point in the sample that has the lowest degree of order S . A small nearby sphere with 3 data points in radius is searched for the lowest S and the line is propagated until it closes the defect loop by meeting the initial point. Visited points are removed from the search space, and the algorithm is repeated until no regions with sufficiently low order are found. This procedure produces polygonal lines that approximate the defect loops, and an associated coordinate frame is constructed, consisting of the tangent ˆ t and two perpendicular vectors ˆ t 1 and ˆ t 2 . The choice of these two vectors is in principle arbitrary, as long as they vary continuously around the loop. For each defect loop segment, the surrounding director field is analysed on a small circle spanned by the perpendicular vectors ˆ t 1 and ˆ t 2 . Two mutually perpendicular normalized vectors ˆ n 1 and ˆ n 2 are found which are used to parametrise the circle on the unit sphere of all directions, visited by the director as traversing the circle. These directions ˆ n 1 and ˆ n 2 are flipped head-to-tail if needed so that they vary continuously along the loop. Keeping the signs and orientation of the local coordinate frame consistent with respect to the initial reference point overcomes the ambiguities inherent in the line field topology. Note that we trace two coordinate frames over the entire loop. The first frame consists of ˆ t , ˆ t 1 , and ˆ t 2 and the second frame consists of ˆ n 1 , ˆ n 2 , and ˆ n 1 × ˆ n 2 . The

orientations of these frames are represented with quaternions instead of rotational matrices, to retain the full topological information in a closed loop [62]. The pair ˆ n 1 and ˆ n 2 define the discs, drawn in Fig. 3b to represent the cross section, while the colour represents the quantity ˆ t · (ˆ n 1 × ˆ ) n 2 / 2 that distinguishes the winding number, from pure -1 / 2 profile in red, to pure +1 2 / in blue. Splay-bend parameter is used to visualise local director around the active defect loops [38].

Velocity field in Fig. 1 is shown by blue arrows above selected cut-off magnitudes of v = 0 0320 . ξ N /τ N in (a) at ζ = 0 003 . L/ξ 2 N , v = 0 127 . ξ N /τ N in (a) at ζ = 0 013 . L/ξ 2 N , v = 0 195 . ξ N /τ N in (a) at ζ = 0 022 . L/ξ 2 N , v = 0 304 . ξ N /τ N in (b), and v = 0 326 . ξ N /τ N in (d). Velocity magnitude is additionally represented by a bright green isosurfaces at similar cutoff values of v = 0 0305 . ξ N /τ N in (a) at ζ = 0 003 . L/ξ 2 N , v = 0 125 . ξ N /τ N in (a) at ζ = 0 013 . L/ξ 2 N , v = 0 0320 . ξ N /τ N in (a) at ζ = 0 192 . L/ξ 2 N , v = 0 297 . ξ N /τ N in (b), and v = 0 316 . ξ N /τ N in (d).

## APPENDIX: SUPPLEMENTARY FIGURES

1

SUPPLEMENTARY FIG. 1. Distribution of defect line segments within the droplet.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00904v1-figure-b1283017-d739-4940-bca7-5dee03c99843.jpg)

Probability density of defect line segments per unit volume, with respect to relative distance from the center of the droplet. The probability is normalized to unity when integrated over the whole volume of the droplet. All three graphs are in the active turbulent regime. At lower activities, the defects are more likely found in the middle of the droplet, whereas with increasing activity the probability becomes more uniform. Immediately adjacent to the surface, probability is somewhat larger due to the flow spreading and extending the defects against the surface.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00904v1-figure-afcda86b-e200-4b65-ad8e-2b8a068c07ec.jpg)

SUPPLEMENTARY FIG. 2. Defect cross sections during two selected rewiring events. Splits (a), merges (b) and cross-overs are all based on local rewiring of two defect line segments. The disks showing the plane of directions the director takes around the defect line are shown to be coplanar between the segments both before and after the rewiring. The cross sections are also topologically complementary - a pair of -1 2 / and +1 2 / profiles (red and blue), a pair of oppositely twisted profiles (purple), or a general intermediate pair. The insets show the director in the plane that separates the defect line segments. In both presented events, the director is planar to the separation plane before the rewiring and perpendicular to the separation plane after the rewiring. Note that each panel shows all defects present in a droplet at a given time, so the total topological charge in each panel must equal +1 (enclosing droplet is not shown). In panel (a), the smaller split-off loop has zero topological charge. In panel (b), the bottom loop before merging has zero topological charge.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00904v1-figure-831502a0-e233-47bf-8cc7-e0bcb2b8f81b.jpg)

SUPPLEMENTARY FIG. 3. Time evolution of defect loop profile in active turbulence regime. At each time, the defect set is represented by vertical strips with length of each segment roughly corresponding to the defect loop length, and the color corresponding to the director profile variation along the loop. In time series, breakups of defects are seen as branching of a single 'stream' into two, coalescence is the opposite of that, annihilations are apex points, and crossovers are poorly represented in this flat projection. Note that defects are closed loops, but to represent them as vertical strips, they are cut at an arbitrary position and aligned with those at a previous time step for best matching.

- [1] G. P. Alexander, B. G. Chen, E. A. Matsumoto, and R. D. Kamien. Colloquium: Disclination loops, point defects, and all that in nematic liquid crystals. Rev. Mod. Phys. , 84:497, 2012.
- [2] D. Beller. Defect loops in 3D active nematics , January 2019. lecture in the workshop 'Optimal design of complex materials' at the Isaac Newton Institute, University of Cambridge, UK; https://www.newton.ac.uk/seminar/20190114161517001.
- [3] A. N. Beris and B. J. Edwards. Thermodynamics of Flowing Systems with Internal Microstructure . Oxford University Press, New York, 1994.
- [4] C. Blanch-Mercader, V. Yashunsky, S. Garcia, G. Duclos, L. Giomi, and P. Silberzan. Turbulent Dynamics of Epithelial Cell Cultures. Phys. Rev. Lett. , 120:208101, 2018.
- [5] J. Brugués and D. Needleman. Physical basis of spindle self-organization. Proc. Natl. Acad. Sci. , 111:18496, 2014.
- [6] B. G. Chen, P. J. Ackerman, G. P. Alexander, R. D. Kamien, and I. I. Smalyukh. Generating the hopf fibration experimentally in nematic liquid crystals. Phys. Rev. Lett. , 110:237801, 2013.
- [7] S. J. DeCamp, G. S. Redner, A. Baskaran, M. F. Hagan, and Z. Dogic. Orientational order of motile defects in active nematics. Nat. Mater. , 14:1110, 2015.
- [8] Z. Dogic. Structure and Dynamics of Three-Dimensional Active Nematics , July 2019. lecture at the GRC Liquid Crystals, NH, USA; https://www.grc.org/liquid-crystals-conference/2019/.
- [9] A. Doostmohammadi, J Ignés-Mullol, J. M. Yeomans, and F. Sagués. Active nematics. Nat. Commun. , 9:045006, 2018.
- [10] A. Doostmohammadi, T. N. Shendruk, K. Thijssen, and J. M. Yeomans. Onset of meso-scale turbulence in active nematics. Nat. Commun. , 8:15326, 2017.
- [11] Perry W. Ellis, Daniel J. G. Pearce, Ya-Wen Chang, Guillermo Goldsztein, Luca Giomi, and Alberto Fernandez-Nieves. Curvature-induced defect unbinding and dynamics in active nematic toroids. Nature Phys. , 14:85, 2018.
- [12] D. Foster, C. Kind, P. J. Ackerman, J.-S. B. Tai, M. R. Dennis, and I. I. Smalyukh. Twodimensional skyrmion bags in liquid crystals and ferromagnets. Nature Phys. , page 178, 2019.
- [13] Tong Gao and Zhaorui Li. Self-driven Droplet Powered By Active Nematics. Phys. Rev. Lett. , 119:564, 2017.

- [14] L. Giomi. Geometry and Topology of Turbulence in Active Nematics. Phys. Rev. X , 5:031003, 2015.
- [15] L. Giomi, M. J. Bowick, X. Ma, and M. C. Marchetti. Defect annihilation and proliferation in active nematics. Phys. Rev. Lett. , 110:228101, 2013.
- [16] Richard Green, John Toner, and Vincenzo Vitelli. Geometry of thresholdless active flow in nematic microfluidics. Phys. Rev. Fluids , 2:114, 2017.
- [17] P. Guillamat, J. Ignés-Mullol, and F. Sagués. Control of active liquid crystals with a magnetic field. Proc. Natl. Acad. Sci. , 113:5498, 2016.
- [18] P. Guillamat, Ž. Kos, J. Hardoüin, J. Ignés-Mullol, M. Ravnik, and F. Sagués. Active nematic emulsions. Science Advances , 4:eaao1470, 2018.
- [19] Raimo Hartmann, Praveen K. Singh, Philip Pearce, Rachel Mok, Boya Song, Francisco DíazPascual, Jörn Dunkel, and Knut Drescher. Emergence of three-dimensional order and structure in growing biofilms. Nature Phys. , 15:251, 2019.
- [20] Y. Hatwalne, S. Ramaswamy, M. Rao, and R. Simha. Rheology of active-particle suspensions. Phys. Rev. Lett. , 92:118101, 2004.
- [21] S. Henkes, M. C. Marchetti, and R. Sknepnek. Dynamical patterns in nematic active matter on a sphere. Phys. Rev. E , 97:3013, 2018.
- [22] M B Hindmarsh and T W B Kibble. Cosmic strings. Rep. Prog. Phys. , 58:477, 1995.
- [23] T Ishikawa and O. D Lavrentovich. Crossing of disclinations in nematic slabs. Europhys. Lett. , 41:171, 1998.
- [24] K. Kawaguchi, R. Kageyama, and M. Sano. Topological defects control collective dynamics in neural progenitor cell cultures. Nature , 545:327, 2017.
- [25] F. C. Keber, E. Loiseau, T. Sanchez, S. J. DeCamp, L. Giomi, M. J. Bowick, M. C. Marchetti, Z. Dogic, and A. R. Bausch. Topology and dynamics of active nematic vesicles. Science , 345:1135, 2014.
- [26] D. Khoromskaia and G. P. Alexander. Vortex formation and dynamics of defects in active nematic shells. New J. Phys. , 19:103043, 2017.
- [27] K. Kruse, J. F. Joanny, F. Jülicher, J. Prost, and K. Sekimoto. Asters, vortices, and rotating spirals in active gels of polar filaments. Phys. Rev. Lett. , 92:078101, Feb 2004.
- [28] T. Krüger, H. Kusumaatmaja, A. Kuzmin, O. Shardt, G. Silva, and E. M. Viggen. The Lattice Boltzmann Method . Springer Interantional Publishing, Switzerland, 2017.

- [29] Enkeleida Lushi, Hugo Wioland, and Raymond E. Goldstein. Fluid flows created by swimming bacteria drive self-organization in confined suspensions. Proc. Natl. Acad. Sci. , 111:9733, 2014.
- [30] C. C. Maass, C. Krüger, S. Herminghaus, and C. Bahr. Swimming Droplets. Annu. Rev. Condens. Matter Phys. , 7:171, 2016.
- [31] Ananyo Maitra, Pragya Srivastava, M. Cristina Marchetti, Juho S. Lintuvuori, Sriram Ramaswamy, and Martin Lenz. A nonequilibrium force can stabilize 2d active nematics. Proc. Natl. Acad. Sci. , 115:6934, 2018.
- [32] M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, Madan Rao, and R. Aditi Simha. Hydrodynamics of soft active matter. Rev. Mod. Phys. , 85:1143, 2013.
- [33] A. Martinez, M. Ravnik, B. Lucero, R. Visvanathan, S. Žumer, and I. I. Smalyukh. Mutually tangled colloidal knots and induced defect loops in nematic fields. Nat. Mater. , 13:258, 2014.
- [34] Berta Martínez-Prat, Jordi Ignés-Mullol, Jaume Casademunt, and Francesc Sagués. Selection mechanism at the onset of active turbulence. Nature Phys. , 15:362, 2019.
- [35] Luuk Metselaar, Amin Doostmohammadi, and Julia M. Yeomans. Topological states in chiral active matter: Dynamic blue phases and active half-skyrmions. J. Chem. Phys. , 150:064909, 2019.
- [36] J. Milnor. Morse Theory . Princeton University Press, Princeton, New Jersey, 1963.
- [37] Daniel Needleman and Zvonimir Dogic. Active matter at the interface between materials science and cell biology. Nat. Rev. Mater. , 2:17048, 2017.
- [38] Maryam Nikkhou, Miha Škarabot, Simon Čopar, Miha Ravnik, Slobodan Žumer, and Igor Muševič. Light-controlled topological charge in a nematic liquid crystal. Nature Phys. , 11:183, 2015.
- [39] D. Nishiguchi, I. S Aranson, A. Snezhko, and A. Sokolov. Engineering bacterial vortex lattice via direct laser lithography. Nat. Commun. , 9:15346, 2018.
- [40] A. Opathalage, M. M. Norton, M. P. N. Juniper, B. Langeslay, S. Ali Aghvami, S. Fraden, and Z. Dogic. Self-organized dynamics and the transition to turbulence of confined active nematics. Proc. Natl. Acad. Sci. , 116:4788, 2019.
- [41] C. Pfleiderer and A. Rosch. Condensed-matter physics: Single skyrmions spotted. Nature , 465:880, 2010.
- [42] G. Posnjak, S. Čopar, and I. Muševič. Hidden topological constellations and polyvalent charges in chiral nematic droplets. Nat. Commun. , 8:14594, 2017.

- [43] J. Prost, F. Jülicher, and J-F. Joanny. Active gel physics. Nature Phys. , 11:111, 2015.
- [44] E. Putzig, G. S. Redner, A. Baskaran, and A. Baskaran. Instabilities, defects, and defect ordering in an overdamped active nematic. Soft Matter , 12:3854, 2016.
- [45] T. Sanchez, D. T. N. Chen, S. J. DeCamp, M. Heymann, and Z. Dogic. Spontaneous motion in hierarchically assembled active matter. Nature , 491:431, 2012.
- [46] T. Beng Saw, A. Doostmohammadi, V. Nier, L. Kocgozlu, S. Thampi, Y. Toyama, P. Marcq, C. T. Lim, J. M. Yeomans, and B. Ladoux. Topological defects in epithelia govern cell death and extrusion. Nature , 544:212, 2017.
- [47] A. Scholich, S. Syga, H. Morales-Navarrete, F. Segovia-Miranda, H. Nonaka, K. Meyer, W. de Back, L. Brusch, Y. Kalaidzidis, M. Zerial, F. Jülicher, and B. M. Friedrich. Quantification of Nematic Cell Polarity in Three-dimensional Tissues. arXiv:1904.08886 , 2019.
- [48] Tyler N. Shendruk, Kristian Thijssen, Julia M. Yeomans, and Amin Doostmohammadi. Twistinduced crossover from two-dimensional to three-dimensional turbulence in active nematics. Phys. Rev. E , 98:010601, 2018.
- [49] S. P. Thampi, A. Doostmohammadi, T. N. Shendruk, R. Golestanian, and J. M. Yeomans. Active micromachines: Microfluidics powered by mesoscale turbulence. Science Advances , 2:e1501854, 2016.
- [50] Maxime Theillard, Roberto Alonso-Matilla, and David Saintillan. Geometric control of active collective motion. Soft Matter , 13:363, 2017.
- [51] E. Tjhung, D. Marenduzzo, and M. E. Cates. Spontaneous symmetry breaking in active droplets provides a generic route to motility. Proc. Natl. Acad. Sci. , 109:12381, 2012.
- [52] U. Tkalec, M. Ravnik, S. Čopar, S. Žumer, and I. Muševič. Reconfigurable Knots and Links in Chiral Nematic Colloids. Science , 333:62, 2011.
- [53] J. Urzay, A. Doostmohammadi, and J. M. Yeomans. Multi-scale statistics of turbulence motorized by active matter. J. Fluid Mech. , 822:762, 2017.
- [54] R Voituriez, J. F Joanny, and J Prost. Spontaneous flow transition in active polar gels. Europhys. Lett. , 70:404, 2005.
- [55] X. Wang, Y.-Ki Kim, E. Bukusoglu, B. Zhang, D. S. Miller, and N. L. Abbott. Experimental Insights into the Nanostructure of the Cores of Topological Defects in Liquid Crystals. Phys. Rev. Lett. , 116:147801, 2016.
- [56] H. H. Wensink, J. Dunkel, S. Heidenreich, K. Drescher, R. E. Goldstein, H. Löven, and J. M.

Yeomans. Meso-scale turbulence in living fluids. P. Natl. Acad. Sci. USA , 109, 2012.

- [57] Carl A. Whitfield, Tapan Chandra Adhyapak, Adriano Tiribocchi, Gareth P. Alexander, Davide Marenduzzo, and Sriram Ramaswamy. Hydrodynamic instabilities in active cholesteric liquid crystals. Eur. Phys. J. E , 40:1226, 2017.
- [58] Hugo Wioland, Francis G. Woodhouse, Jörn Dunkel, and Raymond E. Goldstein. Ferromagnetic and antiferromagnetic order in bacterial vortex lattices. Nature Phys. , 12:341, 2016.
- [59] Kun-Ta Wu, Jean Bernard Hishamunda, Daniel T. N. Chen, Stephen J. DeCamp, Ya-Wen Chang, Alberto Fernández-Nieves, Seth Fraden, and Zvonimir Dogic. Transition from turbulent to coherent flows in confined three-dimensional active fluids. Science , 355:eaal1979, 2017.
- [60] Rui Zhang, Ye Zhou, Mohammad Rahimi, and Juan J. de Pablo. Dynamic structure of active nematic shells. Nat. Commun. , 7:13483, 2016.
- [61] Shuang Zhou, Oleh Tovkach, Dmitry Golovaty, Andrey Sokolov, Igor S Aranson, and Oleg D Lavrentovich. Dynamic states of swimming bacteria in a nematic liquid crystal cell with homeotropic alignment. New J. Phys. , 19:055006, 2017.
- [62] Simon Čopar. Topology and geometry of nematic braids. Phys. Rep. , 538:1, 2014.
Simon Čopar Jure Aplinc Žiga Kos Slobodan Žumer Miha Ravnik
Topology of three-dimensional active nematic turbulence confined to droplets
Active nematics contain topological defects which under sufficient activity
move, create and annihilate in a chaotic quasi-steady state, called active
turbulence. However, understanding active defects under confinement is an open
challenge, especially in three-dimensions. Here, we demonstrate the topology of
three-dimensional active nematic turbulence under the spherical confinement,
using numerical modelling. In such spherical droplets, we show the
three-dimensional structure of the topological defects, which due to closed
confinement emerge in the form of closed loops or surface-to-surface spanning
line segments. In the turbulent regime, the defects are shown to be strongly
spatially and time varying, with ongoing transformations between positive
winding, negative winding and twisted profiles, and with defect loops of zero
and non-zero topological charge. The timeline of the active turbulence is
characterised by four types of bulk topology-linked events --- breakup,
annihilation, coalescence and cross-over of the defects --- which we discuss
could be used for the analysis of the active turbulence in different
three-dimensional geometries. The turbulent regime is separated by a first
order structural transition from a low activity regime of a steady-state vortex
structure and an offset single point defect. We also demonstrate coupling of
surface and bulk topological defect dynamics by changing from strong
perpendicular to inplane surface alignment. More generally, this work is aimed
to provide insight into three-dimensional active turbulence, distinctly from
the perspective of the topology of the emergent three-dimensional topological
defects.