Abstract
Mesoscale turbulence is an innate phenomenon, distinct from inertial turbulence, that spontaneously occurs at low Reynolds number in fluidized biological systems. This spatiotemporal disordered flow radically changes nutrient and molecular transport in living fluids and can strongly affect the collective behaviour in prominent biological processes, including biofilm formation, morphogenesis and cancer invasion. Despite its crucial role in such physiological processes, understanding mesoscale turbulence and any relation to classical inertial turbulence remains obscure. Here we show how the motion of active matter along a microchannel transitions to mesoscale turbulence through the evolution of locally disordered patches (active puffs) from an ordered vortexlattice flow state. We demonstrate that the stationary critical exponents of this transition to mesoscale turbulence in a channel coincide with the directed percolation universality class. This finding bridges our understanding of the onset of lowReynoldsnumber mesoscale turbulence and traditional scaleinvariant turbulence in confinement.
Introduction
LowReynoldsnumber turbulence is established through continuous energy injection from the constituent elements of an active fluid in many biological systems, including bacterial suspensions^{1,2,3,4,5,6}, cellular monolayers^{7,8,9} or subcellular filament/motor protein mixtures^{10,11}. Although the inertia is negligible (Reynolds number ≫1) in such systems, active turbulence is characterized by a highly disordered distribution of vortices^{12,13}. However, mesoscale turbulence in living fluids possesses a characteristic vortex length scale, which distinguishes it from scaleinvariant inertial turbulence^{14}, and it is considered a new class of turbulent flow^{2,12,15}.
Despite extensive implications for diverse fluid dynamical systems and more than a century of research, the transition from pressuredriven laminar flow to inertial turbulence in even the simplest geometries remains one of the major unresolved problems in fluid mechanics. Elaborate experiments have recently shed new light on the nature of this transition by measuring the decay and splitting of local turbulent domains/clusters (puffs) in pipe flows and have determined the critical Reynolds number—ratio of inertial to viscous forces—at which the transition occurs^{16,17}. Shortrange interactions between the locally turbulent puffs, which feed on surrounding laminar flow as an absorbing state, drive a continuous transition to a fully turbulent flow. Recent experimental evidence from channel and circular Couette flows^{18,19}, together with direct numerical simulation studies and predator–prey models^{20}, have provided evidence that the transition at the critical Reynolds number is characterized by the directed percolation universality class.
Strikingly, here we show that for a profoundly distinct class of turbulence at low Reynolds number, the transition in a channel can also be characterized by the emergence of puffs created by microscopic activity of biological fluids. Even for this low Reynolds number class of turbulentlike flows that apparently lacks a perfectly unique absorbing state, we find that the critical exponents correspond to the directed percolation universality class subject to a small nonordering, conjugated and intrinsic field.
Results
Vortex lattice as an intermediate state
To study the transition to lowReynoldsnumber mesoscale turbulence, we computationally solve the continuum equations of active nematics in microchannels, which have successfully reproduced the patterns of bacterial ordering in bulk^{5} and in confinement^{21}, the flow structure and correlation lengths of microtuble bundles^{10,11,22} and the flow patterns of dividing cells^{8,23} (see Methods for the details of the model). Through this continuum description, the transition to turbulence occurs by increasing the amount of local energy injection (activity) in the living fluids. In a confined environment, the activity leads to spontaneous symmetry breaking and the generation of unidirectional flow^{24}, which is followed by an oscillatory regime characterized by distorted streamlines^{25,26}, upon increasing the activity. Further increase in the activity leads to the emergence of a stable lattice of vortices throughout the channel^{27} (Fig. 1a), and this transitions to mesoscale turbulence at higher activities (Fig. 1b). The emergence of the intermediate vortex lattice in active matter has been observed experimentally in motility assays of microtubles^{28}, in bacterial suspension in a channel confinement^{6}, and also numerically by shortrange attraction of selfpropelled particles^{29} and hydrodynamic screening of activityinduced flows due to frictional damping^{30}. To focus on the effect of the confining channel on the transition to mesoscale turbulence, we consider the ideal ‘wet’ limit and neglect additional frictional damping^{11,31,32}. This is because in the experimental systems studied so far, there is no obvious qualitative effect on the active turbulence, friction appears to be a small effect, and the mesoscale turbulent state we consider here is unaffected^{2,6,10,11}. We find that the intermediate vortex lattice is stable to impulsive perturbations, which do not produce growing modes (see Methods). In stark contrast to inertial turbulence, the Reynolds number is irrelevant here and the transition between flow regimes is governed by the dimensionless activity number (Fig. 2a). This parameter characterizes the ratio of the channel height h, which here is equivalent to the hydrodynamic screening length, to the characteristic activityinduced length scale , which represents the relative importance of the intrinsic activity ζ and the orientational elasticity K of the nematic fluid^{33,34}.
The marked difference between the various flow states is clearly seen in the structure of the vorticity. Therefore, to characterize the transition between the regimes, we measure the distribution of the local enstrophy , averaged across the channel. This quantity represents the strength of vortices in the flow, and has also been used for determining the nature of inertial turbulence. The vortexlattice state possesses a welldefined peak in the enstrophy (Fig. 2b). As the active flow transitions at higher activities, the enstrophy distribution broadens, demonstrating that vorticity cascades down into mesoscale turbulence. The gradual disappearance of the peak in the enstrophy distribution (Fig. 2b) suggests a continuous transition from the vortex lattice to mesoscale turbulence. But how does the active turbulence develop from the vortex lattice?
Critical behaviour
Figure 1c shows a snapshot of the vorticity field in a long channel close to the transition. The vortex lattice predominantly occupies the entire channel. Locally, however, we can identify regions of the channel where vortex pairs split into smaller nonordered vortices (Fig. 1d). This coexistence of the global vortex lattice and clusters of local active turbulence controls the transition to turbulence in the channel. We term these localized domains of nonordered vorticity active puffs, in analogy to the inertial puffs observed in the experiments on scaleinvariant turbulence in long tubes^{16}. Unlike the inertial puffs that are externally initiated by perturbations to the flow field (such as induced pressure jumps), active puffs are intrinsic to our simulations.
The behaviour of these active puffs is clearly characterized in the spacetime kymograph of enstrophy (Fig. 3). An active puff can split, giving birth to new puffs, or decay into the ordered vortexlattice state. Below some critical activity number, puffs tend to decay back to the inactive vortexlattice state (Fig. 3a). On the other hand, active puffs span the entire system when splitting occurs at a high enough rate to produce a statistical steady state at the critical activity number (Fig. 3b). At higher activities yet, the competition between splitting and decaying of puffs results in a welldefined turbulence fraction within the channel (Fig. 3c). The active flow approaches the fully turbulent state with active puffs ultimately occupying the entire channel when the decay time far exceeds the splitting time^{16}.
We thus measure the turbulence fraction, the area fraction occupied by active puffs in the channel, as a function of the activity number (Fig. 4a). Well below the critical point, active puffs have a short lifetime and rarely split (Fig. 3a), leading to a negligible turbulence fraction in the steady state (Fig. 4a). However, as the critical value of the activity is approached, puff decay becomes less likely and splitting time decreases substantially (Fig. 3b). Above the critical point, the puff population does not die out, producing a steadystate, nonzero turbulence fraction (Fig. 3c), and we find the turbulence fraction continuously increases with a powerlaw dependence ∼(A−A_{cr})^{β} (Fig. 4a). We measure the stationary exponent to be β=0.275±0.043, which closely matches the universal critical exponent of the (1+1) directed percolation process (β=0.276)^{35} and is in agreement with the value that has recently been measured for inertial turbulence in Couette flow (β=0.28±0.03)^{19}.
This strong agreement is striking, as it draws a parallel between the lowReynoldsnumber mesoscale turbulence in living fluids, which possesses a characteristic vorticity length scale, and highReynoldsnumber inertial turbulence, which is scaleinvariant. Furthermore, the exponent is particularly surprising as active puffs are generated at a small activityindependent rate. In this way, the inactive vortex lattice is not a perfectly absorbing state, and yet there is remarkable agreement between the measured stationary exponent β and the critical exponent of the directed percolation universality class. The agreement with the universal critical exponent from directed percolation can be understood by recognizing that the creation of puffs corresponds to a weak, nonordering field conjugated to the turbulence fraction, which is known to have no detectable effect on the stationary exponents sufficiently close to the critical point^{36,37}.
To further scrutinize the critical behaviour at the transition point, we measure the spatial and temporal distributions of vortexlattice gaps (see Methods). These distributions of the inactive state characterize correlations of the active puffs^{38} and obey power laws with exponents , μ_{} for space and time, respectively (Fig. 4b,c). The temporal exponent is measured to be μ_{}=1.84±0.04 and the spatial exponent is =1.8±0.1. These values also correspond to the critical exponents for (1+1) directed percolation (μ_{}=1.84,=1.748)^{35}. The values of the exponents obtained from our measurements for mesoscale turbulence in a channel and for (1+1) directed percolation with spontaneous site activation are summarized in Table 1 and are compared with the experimentally measured exponents for the inertial turbulence in simple shear experiments in onedimensional geometries^{19}. It would be of interest to see if the critical exponents of the directed percolation universality class will continue to be found in geometries with higher effective dimensionality, as in experiments on inertial turbulence in quasitwodimensional passive liquid crystals^{39,40}, in channel flows^{18} and in threedimensional quantum turbulence^{41}.
Discussion
Our findings present a first concrete connection between turbulence in living fluids and classical scaleinvariant turbulence, beyond a superficial visual similarity, by showing that the transitions to these two profoundly distinct types of spatiotemporal disorder in channel flows are characterized by the same scaling behaviour, corresponding to a critical directed percolation process. While the transition to scaleinvariant turbulence corresponds to a critical absorbing phase transition, our results suggest that mesoscale turbulence is driven away from criticality by rare active puff creation that corresponds to directed percolation in the presence of a conjugated field. This opens new possibilities for further investigation of the nature of mesoscale turbulence and using tools from nonequilibrium statistical mechanics to explain transitionary behaviours in biological systems. Future research should investigate the transitions between ordered flow states, the nature of the weak conjugated field and the possibility of nonuniversal dynamical scaling behaviour or superexponential dependence of puff lifetime on activity number.
Methods
Active nematohydrodynamics simulations
The spatiotemporal evolution of a living fluid is described by active nematohydrodynamics equations based on the theory of liquid crystals. This formulation has been extensively applied to biological systems including bacterial suspensions^{21}, microtuble/motor protein mixtures^{10,22,42} and cellular monolayers^{8,43}. The total density ρ and the velocity field u of the active matter obey the incompressible Navier–Stokes equations
where Π is the stress tensor. While several studies of mesoscale turbulence have characterized the dynamics of the flow using only the velocity field as the relevant order parameter^{2,15}, an additional order parameter field is required to account for the orientational order of active fluids. This is particularly important since several experiments have established the existence and pivotal role of the orientational order in the dynamics of bacterial suspensions^{5,21}, microtuble bundles^{10,11}, assemblies of fibroblast cells^{44} and more recently in stem cell cultures^{9}. To account for the macroscopic orientational order of microscopic active and anisotropic particles, the nematic tensor is considered, where q denotes the coarsegrained magnitude of the orientational order, n is the director and I the identity tensor. The nematic tensor evolves as
where Γ is a rotational diffusivity and the corotation term
accounts for the response of the orientation field to the extensional and rotational components of the velocity gradients, as characterized by the strain rate E=(∇ u^{T}+∇ u)/2 and vorticity Ω=(∇ u^{T}−∇ u)/2 tensors, and weighted by the tumbling parameter λ. The relaxation of the orientational order is determined by the molecular field,
where denotes the free energy. We use the Landau–de Gennes bulk free energy^{45},
and , which describes the cost of spatial inhomogeneities in the order parameter, assuming a single elastic constant K.
In addition to the viscous stress Π^{visc}=2η E, equation (2) must account for contributions to the stress Π from the nematic elasticity and the activity. The nematic contribution to the stress is
which includes the pressure P (ref. 46). The active contribution to the stress takes the form Π^{act}=−ζ Q (ref. 47), such that any gradient in Q generates a flow field, with strength determined by the activity coefficient, ζ.
The equations of active nematohydrodynamics (equations (1, 2, 3)) are solved using a hybrid lattice Boltzmann and finite difference method^{48,49,50}. We model a twodimensional velocity field and nematic director^{22}. While the director field can theoretically develop outofplane components, for active experiments and the parameter values used in this study, it does not move out of plane. Discrete space and time steps are chosen as unity and all quantities can be converted to physical units in a materialdependent manner^{22,51,52}. Simulations are performed with the parameters A=0, B=0.3, C=−0.3, Γ=0.34, K=0.04, λ=0.3, ρ=1 and μ=2/3, in lattice Boltzmann units. Parameter fitting of the continuum equations to physical active systems remains a topic of research; therefore, we consider a generic parameter set that has been shown to reproduce the active turbulent state observed in bacteria and microtuble/motor protein suspensions^{10,22}. With the largest observed average speed and vortex size, these parameters indicate a small Reynolds number Re≲O(10^{−1}). Additional details can be found in refs 48, 49, 50, 53.
We use a channel with a height h=25 and length L=3,000. Noslip boundary conditions are applied to channel walls and periodic boundary conditions are used at the channel extremities. The results reported here are for strong homogeneous boundary conditions for the director field on the channel walls. In addition, we have performed simulations with homeotropic and weak anchoring boundary conditions and find that the transitions described in the main text are independent of the anchoring boundary conditions on the walls. To obtain sufficient statistics for the turbulence fraction (Fig. 4a), several instances (5–10) of each set of parameters with initially randomized fields are simulated. A warmup of 3 × 10^{5} time steps is allowed to reach the stationary state. Above the critical point, this produces sample sizes of ≳100 puffs.
Directed percolation model with spontaneous activation
To examine the behaviour of the (1+1) directed percolation universality class, we utilized a Domany–Kinzel cellular automaton^{54} and chose probabilities to correspond to sitedirected percolation. This stochastic model is a discrete system on a diagonal square lattice of linear spatial size L with periodic boundary conditions. The state s(i,t) of site i at time t can be inactive (or empty) with s(i,t)=0, or the site can be activated (or occupied) s(i,t)=1 (ref. 35). In confined active nematic flows, the inactive state corresponds to the vortexlattice state, and the activated phase corresponds to the active puffs of mesoscale turbulence. At time t, each site is occupied with some probability P_{2} if both backward sites (at time t−1) are occupied, and with probability P_{1} if only one backward site is occupied. Sitedirected percolation is recovered with the choice P=P_{1}=P_{2} (ref. 55). When the probability P reaches the critical probability P_{c}, the system transitions from the absorbing phase of entirely inactive states to one in which the stationary density of active sites is nonzero and this transition is known to belong to the directed percolation universality class^{35}.
In the confined active nematic, we find a small but nonzero rate of puff creation P_{0}. Puff creation destroys the absorbing state, driving the system away from criticality. However, directed percolation in the presence of spontaneous site activation has been well studied, and it is known that the transition point is only weakly dependent on P_{0} (ref. 36), that spontaneous site activation is equivalent to a weak external field conjugated to the order parameter^{35}, and that the stationary universal scalings of the directed percolation universality class continue to be found for sufficiently small spontaneous site activation (P_{0}≲10^{−3})^{56}.
Our directed percolation simulations employ periodic boundary conditions and a lattice size of 10^{4} sites in the spatial dimension to coincide with the lattice Boltzmann system. Data are obtained from 10^{3} runs of 5 × 10^{3} time steps each. We consider P_{0}={0, 10^{−9}, 10^{−8}, 10^{−7}, 10^{−6}} and find the critical probability P_{c}=0.64470±0.00002 as expected (0.6447001(1) (refs 35, 57)). We have found that the rate of puff creation does not depend on the activity coefficient ξ. Thus, at least near the critical point, there does not appear to be a suitable control parameter for P_{0}. From the lattice Boltzmann simulations, we find that the rate of puff formation is (0.54±0.10) × 10^{−7}∼10^{−7}, which is the value used in the main text. Although P_{c} moves the system slightly away from the critical point, the directed percolation scaling exponents are not observed to change and are known to be independent of P_{0} in this weak field limit^{32}. Measuring and N_{}, as for the lattice Boltzmann simulations, supplies the critical exponents reported in the main text.
Calculating the turbulence fraction
The enstrophy field ɛ(x, y, t)=Ω·Ω is calculated from the vorticity field Ω(x, y, t). The field is averaged across the channel to create a 1+1 dimensional kymographical spacetime signal. As described in the main text, the channelaveraged enstrophy in the vortexlattice phase shows regular periodic oscillations, while local active turbulence domains (the active puffs) exhibit fluctuating, noisy enstrophy signals (Fig. 5a). To produce the visual traces of the active puffs (Fig. 3), the kymographs are Fourier transformed in both time and space. The primary peaks are masked in reciprocal spacetime using Gaussian fits to produce the kymographs without the structured oscillations of the periodic background of the vortex lattice (Fig. 5b).
To quantitatively determine the turbulence fraction, the total simulation time is divided into n windows of size τ, creating a L × n array of spacetime intervals, with each window i ranging from t_{i × τ} to t_{(i+1) × τ}. Each of the temporal intervals is analysed for periodicity using a discretized temporal autocorrelation function , where k runs from 0 to τ, , and is the enstrophy signal averaged over the interval. While aperiodic signals decay to zero, oscillating functions possess periodic peaks. An interval is defined to be periodic if the peak amplitude exceeds a threshold value set by the 95 per cent confidence interval for a normal distribution with standard deviation τ^{−1/2}. The time interval τ is chosen between three to five full oscillations of the vortex lattice, which is long enough to detect periodicity but minimizes the number of intervals containing both periodic and aperiodic regions. To reduce the possibility of falsepositives, we recognize that the vortex unit cells have a finite size and duration. We filter the L × n kymograph of periodic/aperiodic windows, demanding that each periodic interval must possess at least one neighbouring periodic interval in time or space; otherwise, it is rejected.
The spatial interval distribution of the inactive state (the vortex lattice) is measured from the processed kymographs by recording the length intervals between active puff regions for fixed temporal coordinates. Similarly, the time interval distribution N_{} between puffs is found by recording the temporal duration of the regions for fixed spatial coordinates. At short spatial intervals L, exhibits oscillations, which represent the characteristic size of the repeating vortexlattice state.
Assessing finitesize effects
To test whether the critical exponents are affected by finitesize effects and physical parameters, we measure the critical exponents for various smaller lengths of the channel than those used in the main text, as well as increased nematic elasticity (Fig. 6). Within the uncertainty of the measurements, the stationary exponent β and the temporal correlations exponent μ_{} are independent of these parameters.
Stability of the vortex lattice
To probe the stability of the vortex lattice, the flow is subjected to instantaneous pulses at different times and locations along the channel. Even near the transition point, the perturbation to the mean enstrophy signal rapidly decays and the vortex lattice is not destroyed. The pulses only instantaneously shift the vortexlattice configuration in space, and are not observed to grow into chaotic regions (active turbulence). Moreover, Fourier analysis of the perturbed enstrophy fields only exhibits the welldefined peaks associated with the vortex lattice and does not reveal growing modes that might indicate linear instability (below, at or above the critical point).
Data availability
Source data and simulation materials are available from the authors upon request.
Additional information
How to cite this article: Doostmohammadi, A. et al. Onset of mesoscale turbulence in active nematics. Nat. Commun. 8, 15326 doi: 10.1038/ncomms15326 (2017).
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
 1
Dombrowski, C., Cisneros, L., Chatkaew, S., Goldstein, R. E. & Kessler, J. O. Selfconcentration and largescale coherence in bacterial dynamics. Phys. Rev. Lett. 93, 098103 (2004).
 2
Wensink, H. H. et al. Mesoscale turbulence in living fluids. Proc. Natl Acad. Sci. USA 109, 14308–14313 (2012).
 3
Sokolov, A. & Aranson, I. S. Physical properties of collective motion in suspensions of bacteria. Phys. Rev. Lett. 109, 248109 (2012).
 4
Zhou, S., Sokolov, A., Lavrentovich, O. D. & Aranson, I. S. Living liquid crystals. Proc. Natl Acad. Sci. USA 111, 1265–1270 (2014).
 5
Nishiguchi, D., Nagai, K. N., Chaté, H. & Sano, M. Longrange nematic order and anomalous fluctuations in suspensions of swimming filamentous bacteria. Phys. Rev. E 95, 020601 (2017).
 6
Wioland, H., Lushi, E. & Goldstein, R. E. Directed collective motion of bacteria under channel confinement. New J. Phys. 18, 075002 (2016).
 7
Vedula, S. K. et al. Emerging modes of collective cell migration induced by geometrical constraints. Proc. Natl Acad. Sci. USA 109, 12974–12979 (2012).
 8
Doostmohammadi, A. et al. Celebrating Soft Matter’s 10th anniversary: Cell Division: a source of active stress in cellular monolayers. Soft Matter 11, 7328–7336 (2015).
 9
Kawaguchi, K., Kageyama, R. & Sano, M. Topological defect launches 3D mound in the active nematic sheet of neural progenitors. Preprint at arXiv:1605.06470 (2016).
 10
Sanchez, T., Chen, D. T. N., DeCamp, S. J., Heymann, M. & Dogic, Z. Spontaneous motion in hierarchically assembled active matter. Nature 491, 431–434 (2012).
 11
Guillamat, P., IgnésMullol, J. & Sagués, F. Control of active liquid crystals with a magnetic field. Proc. Natl Acad. Sci. USA 113, 5498–5502 (2016).
 12
Bratanov, V., Jenko, F. & Frey, E. New class of turbulence in active fluids. Proc. Natl Acad. Sci. USA 112, 15048–15053 (2015).
 13
Giomi, L. Geometry and topology of turbulence in active nematics. Phys. Rev. X 5, 031003 (2015).
 14
Heidenreich, S., Dunkel, J., Klapp, S. H. L. & Bär, M. Hydrodynamic lengthscale selection in microswimmer suspensions. Phys. Rev. E 94, 020601 (2016).
 15
Dunkel, J. et al. Fluid dynamics of bacterial turbulence. Phys. Rev. Lett. 110, 228102 (2013).
 16
Avila, K. et al. The onset of turbulence in pipe flow. Science 333, 192–196 (2011).
 17
Barkley, D. et al. The rise of fully turbulent flow. Nature 526, 550–553 (2015).
 18
Sano, M. & Tamai, K. A universal transition to turbulence in channel flow. Nat. Phys. 12, 249–253 (2016).
 19
Lemoult, G. et al. Directed percolation phase transition to sustained turbulence in Couette flow. Nat. Phys. 12, 254–258 (2016).
 20
Shih, H. Y., Hsieh, T. L. & Goldenfeld, N. Ecological collapse and the emergence of travelling waves at the onset of shear turbulence. Nat. Phys. 12, 245–248 (2016).
 21
Volfson, D., Cookson, S., Hasty, J. & Tsimring, L. S. Biomechanical ordering of dense cell populations. Proc. Natl Acad. Sci. USA 105, 15346–15351 (2008).
 22
Thampi, S. P., Golestanian, R. & Yeomans, J. M. Velocity correlations in an active nematic. Phys. Rev. Lett. 111, 118101 (2013).
 23
Saw, T. B. et al. Topological defects in epithelia govern cell death and extrusion. Nature 544, 212–216 (2017).
 24
Voituriez, R., Joanny, J. F. & Prost, J. Spontaneous flow transition in active polar gels. Europhys. Lett. 70, 404–410 (2005).
 25
Giomi, L., Mahadevan, L., Chakraborty, B. & Hagan, M. F. Banding, excitability and chaos in active nematic suspensions. Nonlinearity 25, 2245 (2012).
 26
Ramaswamy, R. & Jülicher, F. Activity induces traveling waves, vortices and spatiotemporal chaos in a model actomyosin layer. Sci. Rep. 6, 20838 (2016).
 27
Shendruk, T. N., Dootmohammadi, A., Thijssen, K. & Yeomans, J. M. Dancing disclinations in confined active nematics. Soft Matter doi: 10.1039/c6sm02310j (2017).
 28
Sumino, Y. et al. Largescale vortex lattice emerging from collectively moving microtubules. Nature 483, 448–452 (2012).
 29
Großmann, R., Romanczuk, P., Bär, M. & SchimanskyGeier, L. Vortex arrays and mesoscale turbulence of selfpropelled particles. Phys. Rev. Lett. 113, 258104 (2014).
 30
Doostmohammadi, A., Adamer, M. F., Thampi, S. P. & Yeomans, J. M. Stabilization of active matter by flowvortex lattices and defect ordering. Nat. Commun. 7, 10557 (2016).
 31
Oza, A. U. & Dunkel, J. Antipolar ordering of topological defects in active liquid crystals. New J. Phys. 18, 093006 (2016).
 32
Srivastava, P., Mishra, P. & Marchetti, M. C. Negative stiffness and modulated states in active nematics. Soft Matter 12, 8214–8225 (2016).
 33
Thampi, S. P., Golestanian, R. & Yeomans, J. M. Instabilities and topological defects in active nematics. Europhys. Lett. 105, 18001 (2014).
 34
Hemingway, E. J., Mishra, P., Marchetti, M. C. & Fielding, S. M. Correlation lengths in hydrodynamic models of active nematics. Soft Matter 12, 7943–7952 (2016).
 35
Hinrichsen, H. Nonequilibrium critical phenomena and phase transitions into absorbing states. Adv. Phys. 49, 815–958 (2000).
 36
Dickman, R., Rabêlo, W. R. M. & Ódor, G. Pair contact process with a particle source. Phys. Rev. E 65, 016118 (2001).
 37
Lübeck, S. & Willmann, R. D. Universal scaling behaviour of directed percolation and the pair contact process in an external field. J. Phys. A Math. Gen. 35, 10205 (2002).
 38
Kinzel, W. & Yeomans, J. M. Directed percolation: a finitesize renormalisation group approach. J. Phys. A Math. Gen. 14, L163 (1981).
 39
Takeuchi, K. A., Kuroda, M., Chaté, H. & Sano, M. Directed percolation criticality in turbulent liquid crystals. Phys. Rev. Lett. 99, 234503 (2007).
 40
Takeuchi, K. A., Kuroda, M., Chaté, H. & Sano, M. Experimental realization of directed percolation criticality in turbulent liquid crystals. Phys. Rev. E 80, 051116 (2009).
 41
Takahashi, M., Kobayashi, M. & Takeuchi, K. A. Universal critical behavior at a phase transition to quantum turbulence. Preprint at arXiv:1609.01561 (2016).
 42
Giomi, L., Bowick, M. J., Ma, X. & Marchetti, M. C. Defect annihilation and proliferation in active nematics. Phys. Rev. Lett. 110, 228101 (2013).
 43
Bittig, T., Wartlick, O., Kicheva, A., GonzálezGaitán, M. s. & Jülicher, F. Dynamics of anisotropic tissue growth. New J. Phys. 10, 063001 (2008).
 44
Duclos, G., Garcia, S., Yevick, H. G. & Silberzan, P. Perfect nematic order in confined monolayers of spindleshaped cells. Soft Matter 10, 2346–2353 (2014).
 45
de Gennes, P. G. & Prost, J. The Physics of Liquid Crystals Oxford Univ. Press (1995).
 46
Beris, A. N. & Edwards, B. J. Thermodynamics of Flowing Systems: With Internal Microstructure Oxford Univ. Press (1994).
 47
Simha, R. A. & Ramaswamy, S. Hydrodynamic fluctuations and instabilities in ordered suspensions of selfpropelled particles. Phys. Rev. Lett. 89, 058101 (2002).
 48
Marenduzzo, D., Orlandini, E., Cates, M. E. & Yeomans, J. M. Steadystate hydrodynamic instabilities of active liquid crystals: hybrid lattice Boltzmann simulations. Phys. Rev. E 76, 031921 (2007).
 49
Fielding, S. M., Marenduzzo, D. & Cates, M. E. Nonlinear dynamics and rheology of active fluids: simulations in two dimensions. Phys. Rev. E 83, 041910 (2011).
 50
Thampi, S. P., Golestanian, R. & Yeomans, J. M. Vorticity, defects and correlations in active turbulence. Philos. Trans. A Math. Phys. Eng. Sci. 372, 20130366 (2014).
 51
Cates, M. E., Fielding, S. M., Marenduzzo, D., Orlandini, E. & Yeomans, J. M. Shearing active gels close to the isotropicnematic transition. Phys. Rev. Lett. 101, 068102 (2008).
 52
Henrich, O., Stratford, K., Marenduzzo, D. & Cates, M. E. Ordering dynamics of blue phases entails kinetic stabilization of amorphous networks. Proc. Natl Acad. Sci. USA 107, 13212–13215 (2010).
 53
Denniston, C., Marenduzzo, D., Orlandini, E. & Yeomans, J. M. Lattice Boltzmann algorithm for three dimensional liquid crystal hydrodynamics. Philos. Trans. A Math. Phys. Eng. Sci. 362, 1745–1754 (2004).
 54
Domany, E. & Kinzel, W. Equivalence of cellular automata to Ising models and directed percolation. Phys. Rev. Lett. 53, 311–314 (1984).
 55
Lübeck, S. Crossover scaling in the DomanyKinzel cellular automaton. J. Stat. Mech. Theor. Exp. 2006, P09009 (2006).
 56
Lübeck, S. Scaling behavior of the order parameter and its conjugated field in an absorbing phase transition around the upper critical dimension. Phys. Rev. E 65, 046150 (2002).
 57
Jensen, I. Temporally disordered bond percolation on the directed square lattice. Phys. Rev. Lett. 77, 4988–4991 (1996).
Acknowledgements
This work was supported through funding from the ERC Advanced Grant 291234 MiCE and we acknowledge EMBO funding to T.N.S. (ALTF1812013). We thank also Paul van der Schoot.
Author information
Affiliations
Contributions
A.D., T.N.S. and J.M.Y. conceived the research and supervised the project. A.D., T.N.S. and K.T. carried out the simulations. All authors contributed to interpreting the results and writing the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Doostmohammadi, A., Shendruk, T., Thijssen, K. et al. Onset of mesoscale turbulence in active nematics. Nat Commun 8, 15326 (2017). https://doi.org/10.1038/ncomms15326
Received:
Accepted:
Published:
Further reading

Ultrafast viscosity measurement with ballistic optical tweezers
Nature Photonics (2021)

Autonomous materials systems from active liquid crystals
Nature Reviews Materials (2021)

Universal scaling of active nematic turbulence
Nature Physics (2020)

Collective Dynamics of Model PiliBased TwitcherMode Bacilliforms
Scientific Reports (2020)

Topological chaos in active nematics
Nature Physics (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.