Chirality, anisotropic viscosity and elastic anisotropy in three-dimensional active nematic turbulence

Various active materials exhibit strong spatio-temporal variability of their orientational order known as active turbulence, characterised by irregular and chaotic motion of topological defects, including colloidal suspensions, biofilaments, and bacterial colonies.In particular in three dimensions, it has not yet been explored how active turbulence responds to changes in material parameters and chirality.Here, we present a numerical study of three-dimensional (3D) active nematic turbulence, examining the influence of main material constants: (i) the flow-alignment viscosity, (ii) the magnitude and anisotropy of elastic deformation modes (elastic constants), and (iii) the chirality. Specifically, this main parameter space covers contractile or extensile, flow-aligning or flow tumbling, chiral or achiral elastically anisotropic active nematic fluids. The results are presented using time- and space-averaged fields of defect density and mean square velocity. The results also discuss defect density and mean square velocity as possible effective order parameters in chiral active nematics, distinguishing two chiral nematic states -- active nematic blue phase and chiral active turbulence. This research contributes to the understanding of active turbulence, providing a numerical main phase space parameter sweep to help guide future experimental design and use of active materials.


INTRODUCTION
Active fluids are diverse synthetic and biological materials [1][2][3][4][5][6][7][8][9][10].Active nematics are characterized by orientational order along the director vector n and propelled by anisotropic active stress [11,12].A notable and rather ubiquitous dynamic state exhibited by active nematics is active turbulence [13], characterized by continuous spatio-temporal defect proliferation and annihilation.This phenomenon is governed by the intricate interplay between defects within the director field and the structures of the velocity field [14].In three dimensions, the structure of active turbulence comprises a dynamic rewiring network of defect lines and loops [15][16][17][18], where the diverse range of possible defect shapes and their interplay with the velocity field pose significant challenges for its understanding, control and possible applicability.
Approaches for studying three-dimensional (3D) active turbulence include theoretical models of individual defect loops [19][20][21], spectral analysis [16,17], and defect tracking and extracting the statistics of linked loops [22], defect curvature and length [19,23].Another possible approach is to extract the mean-field observables such as average defect density and squared velocity.The scaling of such observables with the activity was determined for single elastic constant and constant viscosity parameters [23,24].However, a systematic study of 3D active * ziga.kos@fmf.uni-lj.sinematic turbulence for different elasticity, viscosity, and chirality material parameters also at different values of activities has not been performed yet.The role of chirality in active turbulence is also relevant as the active matter materials are often weakly chiral [25,26].As shown also for 2D active nematics [27][28][29], a systematic numerical study can lead to better experimental control of active mater with different material properties and serve as a benchmark for theoretical models of activity-dependent dynamic regimes, instabilities, and coupled structures in the flow and orientation field.
In this paper, we show three-dimensional active nematic turbulence at different values of the main material parameters -i.e.(i) the flow-alignment viscosity parameter, (ii) the elastic constants of splay, twist, and bend deformation modes, and (iii) the intrinsic chirality.Specifically, chirality is introduced via chiral elastic energy contribution.From numerical simulations of active dynamics, defect density and mean square velocity are extracted.The simulations are performed for contractile and extensile active materials and show distinct scalings with alignment parameter and average elastic constant, while elastic anisotropy has little effect on the mean-field averaged observables of active turbulence.For chiral active nematic turbulence, we show that increasing inverse chiral pitch q 0 increases the defect density and decreases the mean square velocity.At high values of q 0 , a transition between the active nematic blue phase and chiralityaffected active turbulence is observed.Distinctly, in the active nematic blue phase we observe effective jamming of the defect lattice and a strong drop in the magnitude of the flow field.

RESULTS AND DISCUSSION
We performed numerical simulations of 3D active nematic turbulence using the Beris-Edwards model of nematodynamics and the active stress tensor (Methods).We numerically solve the Q-tensor and the velocity field dynamics on a periodic grid using finite difference and lattice-Boltzmann numerical approaches, respectively [17,24,[29][30][31][32].The main mechanisms of active turbulence that we focus on are visualised in Fig. 1.
Q-tensor field is shown in Fig. 1(a), where black rods show the director field as the main eigenvector of the Q-tensor and the gray isosurfaces show the degree of order representing the main eigenvalue.The degree of order is reduced near the core of disclination lines, which are string-like disordered structures spanning the simulation box. Figure 1(b) shows the velocity field and its magnitude for a selected snapshot of active turbulence.The nematic alignment has three main defor- and bend [n × (∇ × n)] 2 , each with its own elastic constant K 1 , K 2 , and K 3 , respectively.Twist mode can be spontaneously favoured in chiral nematic fluids by incorporating a finite intrinsic chiral pitch length p 0 (Fig. 1(d)).Active force is computed as −ζ∇ • Q and is generated by the splay and bend deformations of the active constituents alignment [26].The force direction in Fig. 1(c) is shown for extensile active materials (ζ > 0), opposite direction is expected for contractile materials (ζ < 0).Anisotropic viscosity of nematic fluids is in the Beris-Edwards model related to the flow-aligning parameter χ, which is typically dependent on the shape of nematic building blocks and describes if they are aligning in flow gradient at a fixed angle, or constantly tumbling (Fig. 1(e)).Note that in principle, a general (passive or active) nematic fluid has 6 viscosity coefficients, 5 of which are independent [33].While the concepts of elastic anistropy, chirality, and flow alignment are well understood in equilibrium and driven nematic fluids, their effects on the irregular state of active turbulence is less understood, particularly in three dimensions.Here, we show how the properties of active turbulence depend on the material coupling constants determining the strength or anisotropy of elasticity, viscosity, and chirality.

Extensile, contractile, aligning, and tumbling nematics
The defect density and the average flow magnitude in active turbulence are affected by the viscosity parameters of the nematic fluid.Here, we vary the flowaligning parameter χ, effectively modelling the flow-tumbling (χ < χ t ) and flow-aligning (χ > χ t ) nematic fluids [33], where the transition between flow-aligning and flow-tumbling regime in our simulations corresponds to χ t = 9Seq 4+3Seq = 0.86.At the highest value of the alignment parameter in the simulation (χ = 1.5), the crosssection of the director field and the velocity magnitude field show deformation at a larger scale compared to the lowest value of the aligning parameter at χ = 0.2 (Fig. 2(a)).During the simulation, both the defect density (Fig. 2(b)) and the mean square velocity (Fig. 2(c)) fluctuate in time due to finite size effect of the simulation box.We consider the effect of the alignment parameter both for the extensile (ζ > 0) and the contractile (ζ < 0) active nematics.Independently from the sign of the active stress, the mean defect density is decreasing with an increase of the alignment parameter χ (Fig. 2(d)).Differently, we observe that the mean velocity decreases with the alignment parameter for contractile active nematics and increases for extensile active nematics (Fig. 2(e)).Different values of the flow-aligning parameter result in different Ericksen-Leslie coefficients of the anisotropic viscosity tensor (see Methods).Similarly to 2D active nematics with substrate friction [34], some change of active behaviour is observed directly at the transition between flow-aligning and flow-tumbling regime at χ t = 0.86.Mean velocity dependence on the alignment parameter χ in Fig. 2(e) gradually changes its slope around χ t .

Role of elastic constants in active turbulence
Simulations at different values of the nematic elastic constants show two main results: (i) the defect density scales approximately linearly with the average elastic constant, and (ii) changing the elastic anisotropydifferent elastic penalties of the splay, twist, and bend director deformation modes-has a minor effect on the active nematic defect density and average velocity magnitude.
Figure 3 shows the defect density and mean square velocity when varying the magnitude of the average elastic constant.We observe that in 3D active turbulence the mean inverse defect density increases with increasing average Frank elastic constant K (Fig. 3(a)), in agreement with the scaling of ρ ∼ |ζ|γ1 ηK , predicted by considering the role of defect density on the line tension, drag force, and defect self-propulsion [24], where γ 1 is the rotational viscosity and η the effective isotropic viscosity.Likewise, the mean square velocity (Fig. 3(b)) also shows a slowly increasing trend with the average elastic constant, in line with the scaling ⟨v 2 ⟩ ∝ |ζ|K, obtained by considering the self-propulsion velocity [19,35] in the tubular neighborhood of a disclination with a radial size of ∼ 1/ √ ρ.The observed behaviour is in agreement with results from twodimensional active turbulence, where the inverse defect density and the mean square velocity are also reported to be proportional to the magnitude of the elastic constant [29,36].
Figure 4 shows the role of elastic anisotropy between the splay (K 1 ), twist (K 2 ) and bend (K 3 ) modes.We change the three elastic constants under the assumption of a fixed average elastic constant K = (K 1 + K 2 + K 3 )/3 (see also Methods).Different values of the elastic constants provide relatively comparable results as seen on colormaps for the mean defect density in Fig. 4(a) and mean velocity in Fig. 4(b).The results in Fig. 4(c) show that the elastic anisotropy in the considered range has little effect on the dynamics of the active turbulence since both mean defect density and mean velocity stay roughly equal under the condition of same K.With the condition of equal splay and twist modes (K 1 = K 2 , in Fig. 4(c), a slight increasing trend of the defect density is observed with increasing K 1 and K 2 and decreasing K 3 .Fig. 4(d) shows the role of anisotropy between splay and bend elastic constants for extensile and contractile systems.For extensile systems, a weakly increasing trend of mean defect density with increasing splay elastic constant can be observed.Contrary, for contractile active nematics, defect density gradually decreases with increasing K 1 at constant average Frank elastic constant.The main difference between extensile and contractile nematics is that the mean velocity is approximately 2-times larger in extensile systems for the same absolute value of activity.Similar behaviour was observed also in 2D active nematics [29].

Chiral active turbulence and transition to active blue phase
Material chirality in nematic systems can emerge as a result of chiral nematic building blocks, chiral dopants, or-in active systems-as a results of chiral dynamics of the active agents [25,26].In bulk 3D chiral active nematic, the intrinsic chirality increases the defect density and reduces the mean square velocity of active turbulence, notably already in the low chirality (i.e.large pitch) regime, as shown in Fig. 5.For example, for the chiral pitch equal to p 0 = 200 ∆x (i.e.89-times the active length l a = L/ζ) at activity ζ = 0.2 L/(∆x) 2 , the defect density is increased by 48% and the root mean square velocity decreased by 42% compared to the dynamic steady state value for achiral active turbulence (Fig. 2(d,e)).The effect of intrinsic chirality on defect density and mean square velocity is larger than the variation of ∼ 30 % that was observed for elastic anisotropy in Fig. 4. One possible explanation is that for elastic anisotropy in Fig. 4, the elastic deformation modes are energetically unfavorable; however, for intrinsic chirality in Fig. 5 nematic twist distortions are energetically preferred, which can have a greater effect on the emerged structure.
Upon further increasing chirality (i.e.reducing pitch) the passive chiral nematic is known to transition into 3D chiral orientational structures known as chiral blue phases I, II and III, and we observe a similar transition in active chiral nematics.Increasing chirality at fixed activity causes a steady increase of the defect density (Fig. 5(d,f,g)), up to a structure where defect lines of the active turbulence jam into a defect network, that one could identify as an effective active blue phase III (Fig. 5(a), first panel).Similarly as for the passive blue phase III [37], structural factor shows that the disclination network in active blue phase has no crystal-like order (Fig. 5(c)).Specifically, in active blue phase at pitch p 0 = 25 ∆x and ζ = 0.2 L/(∆x) 2 , the defect density is 5-times larger than in achiral nematic turbulence case (Fig. 2(d)).An even stronger indication of an effective structural transition is the root mean square velocity (Fig. 5(e,h,i)) that drops to 2% compared to the achiral case at ζ = 0.2 L/(∆x) 2 .For the numerical material parameters of our simulation, the transition to a passive blue phase is observed at a critical pitch length P c ≈ 40 ∆x.Activity reduces this critical pitch length to P c = 34 ∆x at ζ = 0.15 L/(∆x) 2 and even further to P c = 32 ∆x at ζ = 0.2 L/(∆x) 2 .How P c is obtained from the data is explained in Methods.
The mean square velocity drop-off and effective jamming of the disclination network can be explained by the difference of the structural features of the disclination network between active blue phase and active turbulence, as shown in Fig. 5(b).The local director crosssection profile of the disclination lines in the active blue phase are close to the −1/2 winding number, for which the active self-propulsion is well known to be zero [19].Selected double-twist cylinders are also shown in Fig. 5(b), which also do not generate a self-propulsion flow [38].Nodal points, where 4 disclination segments meet (Fig. 5(b)), are further characteristic features of the blue phase III [37].Dynamically, we observe that in the active blue phase at P 0 = 25 ∆x, the nodal points are stable for up to ∼ 150 (∆x) 2 /ΓL.Contrary, in the active chiral turbulence, nodal points appear during reconfiguration events between disclination lines and are much shorter lived -at P 0 = 150 ∆x each nodal point disappears on average after approximately ∼ 13 (∆x) 2 /ΓL.
We observe the distinct scaling of the disclination density and the mean square velocity near the transition between the active blue phase and the active turbulence.In Figure 5(d,e), we plot the inverse defect density and root mean square velocity dependence on the pitch length in proximity of the critical defect density ρ c , critical mean velocity v c and critical chiral pitch P c of the blue phase-active turbulence transition.The slope of the graph shows that inverse defect density scales roughly as (P 0 − P c ) 0.6 and mean velocity as (P 0 − P c ) 0.3 .Similar scaling is obtained for two different activities with a notable difference that at higher activity the blue phase-active turbulence transition occurs at smaller pitch lengths.

CONCLUSIONS
We explore dynamic reconfiguring network of disclination lines known as active turbulence using numerical simulations for selected main material parameters: chirality, flow alignement (anistropic viscosity), and elastic anistropy (different nematic elastic constants).The difference between extensile and contractile systems is observed in increasing or decreasing dependence of mean square velocity with the alignment parameter, showing the importance of the shape of active nematic building blocks [39].We confirm that defect density and mean square velocity are approximately inversely proportional to the magnitude of the average elastic constant, whereas the elastic anisotropy has only a small on the defect density and the mean square velocity of active turbulence.Though, we speculate that elastic anisotropy could affect the local structure of the defect lines (+1/2, −1/2 and twisted [19]) in the defect network.As the elastic instability is known to drive structural reconfigurations in passive liquid crystals [40], we expect that for active nematics elastic instability might have a greater effect in confined systems at the onset of active turbulence.While current experiments on 3D active turbulence explore bulk behaviour [15], the role of elastic anisotropy and alignment parameter is relevant also for possible future results in confined systems.
We performed simulations of an active chiral nematic and show the effective structural transition between chiral active turbulence and the active blue phase.The structures are distinct from each other and are separated by a continuous structural phase transition that we characterize by measuring the average inverse defect density and mean square velocity.Beyond the results reported here, the observed active blue phase dynamics would be (i) interesting to compare with blue phase dynamics due to thermal fluctuations [41] and (ii) explored in the context of driving with external field, such as electric or magnetic fields [42] or activity gradients [43].An interesting aspect is also the difference between the transition to active turbulence in blue phases and in modulated cholesteric phase, for which a linear hydrodynamic instability was predicted for extensile materials [26].While experimentally engineering 3D active blue phase materials can lead to very novel materials and phenomena, the effects of weak chirality that we show in the paper might be important also for present active nematic materials, since building blocks and processes in active matter are often chiral [25,26], which could be even further amplified by introduction of chiral dopants.Additionally, chiral symmetry breaking allows for additional active stresses [44][45][46][47].An open problem for future research is also the possibility of an active blue phase with symmetries of blue phase I or blue phase II [48].
More generally, this work is a contribution towards understanding the material-dependence of active nematic turbulence, to aid the experimental design and theoretical advances of active nematic phase [49][50][51].

Model equations of active nematodynamics
We simulate mesoscopic continuum description of active nematics using the adapted Beris-Edwards approach for active nematodynamics [30,32,52].The nematic order is described by a traceless tensor order parameter Q ij , with the director n as the main eigenvector, and evolves as where v is the fluid velocity and Γ is the rotational viscosity coefficient.The generalized advection term S ij couples the nematic order and fluid velocity where D ij and Ω ij are symmetric and antisymmetric part of the velocity gradient tensor W ij = ∂ i v j , and χ is the alignment parameter.The molecular field H ij drives sys- where F is the Landau-de Gennes free energy Here, A, B and C are material parameters and L 1 , L 2 , and L 3 are elastic constants in the tensorial formulation of the elastic free energy, and q 0 = 2π/P 0 is the inverse chiral pitch.L 1 , L 2 , and L 3 can be computed from the elastic constants of the splay K 1 , twist K 2 , and bend K 3 deformation modes that are formulated within the director-based Frank free energy.Often, single elastic constant approximation (L 1 = L, L 2 = L 3 = 0) is used, where splay, twist and bend elastic modes have equal contributions (K 1 = K 2 = K 3 = K).In Figs. 3 and 4, we explore the role of elastic anistropy, i.e. how individual elastic modes influence the active nematic dynamics and we use non-zero L 1 , L 2 and L 3 elastic constants, following the relations In the formulation of the elastic constants in Eqs. 2, K 24 is equal to K1  2 but is not relevant due to periodic boundary conditions.The flow field obeys the incompressibility condition and the Navier-Stokes equation, where ρ is the fluid density and Π ij is the stress tensor, which consists a passive and an active term , where Here, p is the fluid pressure, η is the isotropic viscosity, χ is the flow alignment parameter and ζ is the activity, which is positive in extensile materials and negative in contractile materials.The coupled equations for the nematic order Q ij and fluid velocity v i are numerically solved using the hybrid lattice-Boltzmann approach [17,24,[29][30][31][32], based on the finite difference method for solving the Q-tensor evolution (Eq.1), and the D3Q19 lattice Boltzmann method for the incompressibility and the Navier-Stokes equation (Eqs.3,4).

Material parameters
In the paper, we consider a single elastic constant approximation for the material constants (L 1 ̸ = 0, L 2 = 0, and L 3 = 0), except in Figs. 3 and 4. Changing the values of the elastic constants affects the nematic correlation length and in turn the resolution of the numerical mesh resolution.Accordingly, we compute the average Frank elastic constant from the mapping in Eqs. 2. To account for simulation results at different values of K in Fig. 3, we set L 1 = L, vary the elastic constant L 2 , and express the simulation results with a constant term Here, L represents the fixed value of the elastic constant as used in Figs. 2 and 5.In Figure 4, where we explore the role of the elastic anisotropy, we use the condition of a constant average Frank elastic constant K and choose the Frank elastic constants accordingly and express the simulation results in terms of K 0 .From a given set of the Frank elastic constants, the tensorial elastic constants L 1 , L 2 , and L 3 are determined from a mapping given by Eqs. 2 evaluated at S = S eq .
The overall results of the simulations are expressed in the units of the elastic constant L, rotational viscosity parameter Γ and mesh resolution ∆x.Mesh resolution is defined as ∆x = 1.5χ n , where χ n is nematic correlation length defined as χ n = L/(A + BS eq + 9 2 CS 2 eq ), where S eq = 0.533 is equilibrium degree of nematic order, A = −0.43L/(∆x) 2 , B = −5.3L/(∆x) 2 , C = 4.325L/(∆x) 2 and η = 1.38/Γ.The size of the simulation box is 201 × 201 × 201 mesh points and periodic boundary conditions are used in all three spatial directions.200 × 200 × 200 simulation box is used in Fig. 2. The time step in simulations is set to ∆t = 0.025(∆x) 2 /(LΓ).To recover the simulation results in the SI units, one can use typical parameter values for active systems L = 3 pN, ∆x = 2 µm, and Γ = 10 (Pa s) −1 , roughly estimated from active turbulence in bacterial and microtubule systems [2,54] and viscoelastic properties of lyotropic nematics [55].

Data analysis
We define the defect density as the length of defect lines over a unit volume and compute it from the defect volume fraction of the regions where scalar order parameter is S < 0.4 [24].In the next analysis step, we average either the defect density (Figs. 2 and 4) or the inverse defect density (Figs. 3 and 5) in time and calculate its standard deviation that is presented with error bars.For the velocity field, we compute the average of the velocity squared ⟨v 2 ⟩ V over the complete simulation volume V at given time.In Figs. 2, 4, and 5, we obtain the root mean square velocity as ⟨v⟩ = ⟨v 2 ⟩ V and the variability (presented by error bars) from its standard deviation.In Fig. 3, the mean inverse velocity squared ⟨1/v 2 ⟩ and its standard deviation are obtained from the time-dependence of the 1/⟨v 2 ⟩ V .To estimate the value of the critical pitch for the transition from the active chiral turbulence to the active blue phase, we performed nonlinear regression on simulation data for mean velocity and mean inverse defect density in Fig. 5.The value of P c with the precision of 1 ∆x was determined so that the mean velocity in the log-log plot (Fig. 5(i)) shows a power-law scaling for the longest range of the pitch values.Once P c was set, we employed a similar procedure to determine 1/ρ c based on the log-log plot of the mean inverse defect density (Fig. 5(g)).

FIG. 1 .
FIG. 1. Active nematic turbulence in a three dimensional system.(a,b) A snapshot of an active turbulence showing (a) a disclination network (gray isosurfaces) with a director field cross-section (black rods) and (b) a cross section of the velocity field.(c-e) A schematic representation of mechanisms of (chiral) active turbulence: (c) different elastic modes of director deformations-splay, twist and bend-with the direction of the active force generated by the splay and bend distortions, (d) chiral structure of the director field with pitch length P , and (e) flow alignment and flow tumbling regime of the director field in shear flow.

FIG. 2 .
FIG. 2. Dynamics of three-dimensional active nematic turbulence for extensile and contractile materials at different aligning parameters.(a) Cross-section of the velocity field magnitude (color map) and the director field (black rods) at two different activities (ζ = 0.2 L/(∆x) 2 and ζ = −0.2L/(∆x) 2 ) and two different alignment parameters (χ = 0.2 and χ = 1.5).(b,c) Defect density and volume-averaged root mean square velocity over time in a dynamic steady state for activity ζ = 0.1 L/(∆x) 2 and different alignment parameters, with the lowest plotted value of χ = 0.0 in dark blue colour and the highest value of χ = 2.0 in the dark red colour.(d,e) The dependence of the mean defect density and the mean velocity on the alignment parameter χ.For ζ = 0.1(∆x) 2 , the mean defect density and the mean velocity are directly computed from panels (b) and (c), respectively.The mean defect density decreases with the alignment parameter for extensile and contractile materials.The mean velocity shows an increasing trend with the alignment parameter for extensile materials (ζ > 0) and a decreasing trend for contractile materials (ζ < 0).Error bars in (d,e) represent the standard deviation originated by the time-averaged defect density and root mean square velocity, respectively.

FIG. 3 .
FIG. 3. Defect density and mean square velocity as dependent on the average elastic constant.(a) Steady state mean inverse defect density and (b) steady state mean inverse velocity squared, both for three different activities.The average elastic constant is expressed in units of K0 and is computed for different twist vs splay and bend elastic constants (see Methods).Error bars represent the standard deviation originated by the time-averaged quantities.

FIG. 4 .
FIG. 4. Effect of elastic anisotropy-different elastic constants-on three-dimensional active turbulence.(a) Mean defect density, (b) root mean square velocity and (c) selected cross-sections of the heat maps in (a,b), for the same average Frank elastic constant K = K0 and of fixed activity ζ = 0.2 L/(∆x) 2 .The data is represented using Frank elastic constants K1 and K2 (note that constant K implies K3 = 3K − K1 − K2).Small variability is observed, which is comparable to the fluctuations of the defect density and the mean velocity in time, as inherent to the chaotic active turbulent dynamics in finite volume.(d) Mean defect density and mean velocity of extensile (ζ = 0.2 L/(∆x) 2 ) and contractile (ζ = −0.2L/(∆x) 2 ) systems for constant twist mode (K2).Error bars in (c,d) represent the standard deviation originated by the time-averaged quantities.

FIG. 5 .
FIG. 5. Chiral active turbulence and active blue phase.(a) Defect network at three different values of chiral pitch P0.(b) Zoom of defect network at two different values of chiral pitch P0, with cross-section of director field in red.(c) Structure factor |S(k)| 2 at two different values of chiral pitch P0, on cuts along ky = 0.The structure factor was calculated from the Fourier transformation of the highest eigenvalue of the Q-tensor.(d) Mean inverse defect density dependence on the chiral pitch P0 for three different activities.Dashed line shows the value of the transition in a passive nematic at Pc ≈ 40 ∆x.(e) Root mean square velocity dependence on the chiral pitch P0. (f, h) A zoomed-in area of mean inverse defect density and root mean square velocity from P0 = 25 ∆x to P0 = 60 ∆x.(g, i) Log-log plot of mean inverse defect density and root mean square velocity, respectively.The data in the log-log plots is plotted relative to the critical pitch length Pc and the critical defect density.The values of Pc = 32 ∆x and 1/ρc = 160 ∆x were chosen for the activity ζ = 0.2 L/(∆x) 2 , and the values of Pc = 34 ∆x and 1/ρc = 170 ∆x for the activity ζ = 0.15 L/(∆x) 2 .The dashed line in (g) has a slope of x 0.6 and the dashed line in (i) a slope of x 0.3 .Error bars in (d-i) represent the standard deviation originated by the time-averaged quantities.