Next Article in Journal
Preface to the Special Issue on “Applications of Fuzzy Optimization and Fuzzy Decision Making”
Previous Article in Journal
“Realistic Choice of Annual Matrices Contracts the Range of λS Estimates” under Reproductive Uncertainty Too
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shearing Effects on the Phase Coarsening of Binary Mixtures Using the Active Model B

by
Antonio Lamura
1,* and
Adriano Tiribocchi
2
1
Istituto Applicazioni Calcolo, CNR, Via Amendola 122/D, 70126 Bari, Italy
2
Istituto per le Applicazioni del Calcolo, CNR, Via dei Taurini 19, 00185 Rome, Italy
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(23), 3008; https://doi.org/10.3390/math9233008
Submission received: 8 November 2021 / Revised: 20 November 2021 / Accepted: 22 November 2021 / Published: 24 November 2021
(This article belongs to the Section Mathematical Physics)

Abstract

:
The phase separation of a two-dimensional active binary mixture is studied under the action of an applied shear through numerical simulations. It is highlighted how the strength of the external flow modifies the initial shape of growing domains. The activity is responsible for the formation of isolated droplets which affect both the coarsening dynamics and the morphology of the system. The characteristic dimensions of domains along the flow and the shear direction are modulated in time by oscillations whose amplitudes are reduced when the activity increases. This induces a broadening of the distribution functions of domain lengths with respect to the passive case due to the presence of dispersed droplets of different sizes.

1. Introduction

Over the last decades much research has been dedicated to studying active materials, i.e., non-equilibrium systems in which internal units continuously consume energy, usually stored in the environment, to self-propel [1,2,3,4,5]. Examples of active matter range from bacterial and algal suspensions [6,7,8,9] to colloidal particles acquiring motion through chemical reactions [10,11,12], up to the cytoskeleton of living cells [13,14,15], to name but a few. Importantly, to convert their energy into motion, particles violate time-reversal symmetry (TRS) at the micro-scale.
An active process of particular relevance to us is the motility-induced phase separation (MIPS), where a suspension of repulsive motile particles can phase separate into bulk dense (liquid) and dilute (vapor) regions [16,17,18]. This essentially occurs because active particles aggregate where they move slowly, thus leading to a local density increase causing a slowing down and further accumulation. The resulting clusters show highly dynamical particle exchange leading to macroscopic phase separation [19,20,21]. The kinetics of MIPS has been found to share many features with that of passive systems with attractive interactions, thus suggesting that TRS could be restored at the macroscale level [16,17,22,23,24,25,26,27]. However, large scale simulations of active Brownian particles [20,21,28,29] as well as experiments [10,11,12] reveal that MIPS may also exhibit further non-equilibrium features, such as the formation of mesoscopic vapor bubbles within the particles aggregates and microphase separation [20,28,29,30]. This scenario would indeed support the view that time-reversal symmetry is manifestly broken macroscopically.
Alongside particle-based simulations, the phenomenology of MIPS can be also described using continuum field theories, either by an explicit coarse-graining of the microscopic dynamics or via symmetry arguments and conservation laws [31]. The latter strategy has led to the construction of the Active Model B (AMB) [24,28,32], a field theory that, in the absence of solvent, essentially extends the passive model B [33] by incorporating an active gradient term that cannot be obtained from a free-energy functional, thus breaking time-reversal symmetry. The passive Model B (MB) falls within a class of phase-field models describing the kinetics of phase separation through a conserved scalar order parameter ϕ capturing, for example, the density of colloids or the concentration in a binary mixture. If hydrodynamics can be neglected, the evolution of ϕ is governed by a Cahn-Hilliard equation [34,35], where the thermodynamic force driving the relaxation of ϕ stems from the functional derivative of a ϕ 4 free energy with square gradient terms. Such simplified description captures, for example, the result R t α , for the dependence of the domain size R on time t with a characteristic growth exponent α = 1 / 3 [36].
MB provides a robust framework to describe the passive phase separation subject to an external shearing [37,38,39,40,41,42,43]. In that case the growth is anisotropic, being characterized by domains stretched and tilted along the flow direction as observed in simulations [42,44] and experiments [45]. Indeed, numerical solutions of the dynamical equation with shear [42] confirm the presence of anisotropy of growing domains, which are found to exhibit two typical lengths. Theoretical calculations based on a renormalization group analysis shows that the growth exponent α y in the shear direction, perpendicular to the flow, is not changed by the applied velocity profile, while the exponent α x along the flow direction is increased by 1 [41]. Moreover, relevant physical quantities are found to oscillate on a logarithmic time scale, due to a periodic stretching and breaking-up of domains. The interest towards the role of shear in phase ordering is still vivid, as witnessed by recent studies [46,47].
In the context of active matter, while efforts have been addressed to investigate the active phase separation in the absence of an external driving [20,21,28,48,49], much less is known about the dynamic response observed when a shear flow is applied. Though it was shown [50] that in AMB the transition from homogeneous to bulk phase separation belongs to the same universality class of equilibrium MB, numerical studies of coarsening in AMB provide different scenarios. In Refs. [32,51], for example, it was argued that in AMB without shear the size of growing domains is still compatible with a growth exponent α = 1 / 3 (i.e., activity has negligible effects on the coarsening dynamics) with a change in the static phase diagram owing to a pressure jump, even in the case of a flat interface. On the contrary, very recent simulations [52,53] put forward the existence of a late-time growth exponent α = 1 / 4 , a result akin to that obtained using particle-based simulations [20,54].
In this work we aim at characterizing, in the AMB theory, the shear-induced morphology of growing domains in two spatial dimensions focusing on the role played by the external flow and the activity. Our results show that, in agreement with the AMB without shear, the activity promotes the formation of isolated drops. This significantly affects the dynamic behavior of the mixture subject to a shear flow. Indeed, unlike the passive counterpart, the amplitude of the time oscillations of domain sizes diminishes for increasing values of activity, although their growth exponents are only weakly modified. In addition, the probability distribution functions of the lengths of patterns computed along the two spatial directions are found to broaden as the activity augments.
The paper is structured as follows. In Section 2 the model is properly defined. The numerical results are presented and discussed in Section 3. The phenomenology of phase separation is discussed in weak and strong regimes, remarking the differences with respect to the passive model. The role played by isolated droplets appearing during coarsening is highlighted in connection with the time evolution of the size of domains along the two directions. Finally, some conclusions are drawn.

2. The Model

We consider a system with a scalar order parameter ϕ ( r , t ) at position r and time t. In the framework of the present model, the order parameter can be considered as the deviation of the density with respect to a reference value. It satisfies a conserved dynamics given by the following equation
ϕ t + · ( ϕ v ) = · J .
Thermal fluctuations are not considered here as it is usually done when studying phase separation [36]. At the l.h.s. a convective term couples ϕ to an imposed linear shear flow. This has the form v = γ ˙ y e x where γ ˙ is the shear rate, y is the coordinate along the (shear) y direction and e x is the unit vector along the (flow) x direction. The addition of this term to the original equation of the AMB is the main novelty of this study which allows us to consider the effect of an external velocity profile on the active coarsening dynamics. Hydrodynamic effects are not taken into account since the evolution Equation (1) is not coupled to the Navier-Stokes equation. This is because we are interested in considering only diffusive phase separation. Indeed, to account for the fluid motion in an active medium, the dynamics of ϕ needs to be coupled to that of a momentum-conserving solvent, whose evolution is governed by the Navier-Stokes equation including an active contribution in the stress tensor [55].
The current J is proportional to the negative gradient of a chemical potential
J = Γ μ
where the mobility Γ is set to unity in the following. The chemical potential μ = μ P + μ A is the sum of passive and active contributions given, respectively, by
(3) μ P = a ϕ + b ϕ 3 κ 2 ϕ , (4) μ A = λ ( ϕ ) 2 .
The term (3), where a , b , κ are positive and of the order unity, corresponds to the passive contribution of the model B, and can be obtained from the derivative of the square gradient φ 4 free-energy functional. In the passive case, the system would demix in two coexisting states at binodal densities ϕ P 1 , P 2 = ± ϕ P with ϕ P = a / b , corresponding to the minima of the free-energy density. The energy cost for the formation of interfaces is proportional to κ . The active contribution (4), on the contrary, cannot be obtained from a proper free energy and is the simplest choice at second order in gradients. It was introduced in the so-called active model B [32] in order to explicitly break TRS. For completeness we note that this term is similar to the one appearing in other nonlinear partial differential equations. Among others we cite, e.g., the Kardar-Parisi-Zhang equation for nonlinear interfacial diffusion [56], the Hunter-Saxton equation for nematic liquid crystals [57], and the Kuramoto–Sivashinsky equation for instabilities in laminar flame front [58,59]. Such a term is controlled by the active parameter λ which can be adjusted to go from the passive case ( λ = 0 ) to the active one ( λ 0 ) with λ O ( 1 ) [32]. We remark that the presence of the λ term is such that the system is invariant under the transformation ( ϕ , λ ) ( ϕ , λ ) . For this reason we restrict our attention to the case with λ 0 . For the active model B the binodals ϕ 1 , 2 ( ϕ 1 < 0 < ϕ 2 ) depend on λ and can be calculated by using the method put forward in Ref. [32] and generalized in Ref. [51]. The spinodals do not depend on λ [51] and are located at ϕ S 1 , S 2 = ± ϕ S with ϕ S = a / ( 3 b ) as in the model B. The state with uniform density is locally unstable between the spinodals while it is metastable between the spinodals and the binodals. When λ increases, the binodal gets closer to the spinodal on the negative ϕ side and stays far on the other side.
Equation (1) is solved in two dimensions by using a finite-difference scheme. The field ϕ is discretized on the nodes ( x i , y j ) ( i , j = 1 , 2 , , L ) of a square lattice with L × L nodes and mesh size Δ x . The time is discretized in time steps Δ t with time values given by t n = n Δ t , n = 1 , 2 , 3 , . Any discretized function h at time t n on a node ( x i , y j ) ( i , j = 1 , 2 , , L ) of the lattice is denoted by h ( x i , y j , t n ) = h i j n . At each time step we update ϕ n ϕ n + 1 using an explicit first-order Euler algorithm for the time derivative [60]
ϕ n + 1 = ϕ n Δ t ( γ ˙ y x ϕ n + α J α n ) .
Central-difference schemes [61] are coded for the spatial derivatives. The x derivative is given by
x h | i j n = h ( i + 1 ) j n h ( i 1 ) j n 2 Δ x
and analogously for the y derivative. The Laplacian operator appearing in the chemical potential μ P is discretized as
2 h | i j n = h ( i + 1 ) j n + h ( i 1 ) j n + h i ( j + 1 ) n + h i ( j 1 ) n 4 h i j n Δ x 2
Periodic boundary conditions (BC) are adopted along the flow direction and Lees–Edwards BC [62] are used in the shear direction. The latter take into account the space shift γ ˙ L Δ x Δ t occurring in a time step, due to shear, between the lower and upper row of the lattice. This is essentially done by identifying a point ( x i , y 1 ) on the lower row of the lattice with the one placed on the upper row at ( x i + γ ˙ L Δ x Δ t , y L ) ( i = 1 , 2 , , L ).

3. Results

In the following, all results are obtained with Δ x = 1 and Δ t = 0.001 , values ensuring numerical stability. We checked that simulations are converged for these values of space and time discretization units by looking at the behavior of the binodal values and of averaged domain sizes for different initial conditions (see the following). We verified that results are stable upon decreasing the discretization units. To this purpose we considered some runs with Δ x = 0.5 and Δ t = 0.0005 , and no significant differences were found in the binodals as well as in the typical extensions of forming patterns during coarsening. The model parameters in Equation (1) are a = b = 1 / 4 and κ = 1 while λ is varied within the range [ 0 , 3 ] . The shear rate γ ˙ is changed to access different shear regimes in the phase-separation process. To this purpose we introduce a dimensionless shear rate γ ˙ ^ = γ ˙ t D where t D is the interface diffusion time. Weak and strong shear regimes are characterized by values of γ ˙ ^ lower and higher than unit, respectively [43].

3.1. Planar Interface

In the case of the model B one has t D = ξ 3 ϕ P 2 Γ σ [63], where ξ = 2 2 κ a is the width of the planar interface between two coexisting phases described by the function ϕ ( x ) = ϕ P tanh ( 2 x / ξ ) , and σ = 2 3 ϕ P 2 2 κ a is the interface tension [36]. For our choice of the parameters, it is t D = 384 in model units. However, for the active model B with λ > 0 there are no explicit expressions for ξ and σ as well as for the binodals, though it can be shown that there is a solution with a planar interface between the coexisting phases at densities ϕ 1 , 2 [32]. For this reason we first compute numerically the binodals ϕ 1 , 2 ( λ ) by considering the relaxation of a flat interface between two states initially set at the values ± ϕ P . Once the values ϕ 1 , 2 ( λ ) are found, a sharp profile between the initial states set at ϕ 1 , 2 ( λ ) is let to evolve to the steady interface. The stationary profiles are shown in Figure 1 for the values λ = 0 , 1 , 2 , 3 . Numerical data are successfully fitted by a kink profile (see the Appendix A for further details). The values of the interface width ξ obtained from fits for different values of the activity are presented in Figure 2 and show a linear dependence of ξ on λ for λ 1 . The binodal densities ϕ 1 and ϕ 2 increase with λ as ϕ 1 approaches the spinodal ϕ S 1 = ϕ S , being ϕ S = 3 / 3 when a = b (like in our case).
The diffusion time t D is computed as the time for the initial sharp profile to relax to the steady interface (see the Appendix A). Numerical results give t D 4 × 10 2 at λ = 0 in agreement with the theoretical estimate of the model B. When λ > 0 , it is found that t D increases with λ reaching the highest value t D 1.1 × 10 3 when λ = 3 . The values of t D are plotted in the right panel of Figure 2 and show that t D increases linearly with λ up to λ 1.5 , before slowing down its growth. In the following we consider a range of values of the shear rate between γ ˙ w = 0.9 × 10 3 and γ ˙ s = 3.6 × 10 3 . The values γ ˙ w and γ ˙ s are such to access weak ( γ ˙ w t D ( λ ) 1 ) and strong ( γ ˙ s t D ( λ ) > 1 ) shear regimes, respectively, for all considered values of λ . This guarantees that the shear regime is not affected when changing activity while keeping fixed the strength of applied flow.

3.2. Phase Separation under Weak and Strong Shear

Now we move on to the study of the phase separation under an external shear flow by varying the activity parameter λ for γ ˙ = γ ˙ w , γ ˙ s . We consider a system initially prepared in a symmetric disordered state with φ ( r , 0 ) = ω where ω is a random number in the range [ 0.01 , 0.01 ] . This state corresponds to a critical composition of the system. The size of the lattice is L = 1024 and measures are averaged over five independent runs. The different values of the dimensionless shear rate γ ˙ ^ influence the initial morphology of the forming domains. This can be seen in Figure 3 and Figure 4 where snapshots of systems at consecutive times are shown for the cases at weak ( γ ˙ t D = 1.0 ) and strong shear ( γ ˙ t D = 4.0 ), respectively, with λ = 3 .
At weak shear, domains form and grow while the density ϕ attains the steady binodal values before shear effects come into play at values of strain γ ˙ t > 1 . This occurs since the relaxation time scale t D ( λ ) is smaller than the shear characteristic time 1 / γ ˙ . Once isolated drops at density ϕ 1 are produced ( γ ˙ t = 1 ), the applied flow advects them along the x direction and favours their merging, thus causing the formation of elongated domains along the y direction ( γ ˙ t = 1.75 ). Afterwards the shear stretches these structures which are tilted as well, being characterized by different thicknesses. Isolated droplets at density ϕ 1 can be observed inside domains at density ϕ 2 ( γ ˙ t = 5.3 ). Later on, the shear further deforms and tilts the domains, which may eventually break up. After the overstretching, domains retract forming ϕ 2 phases with a larger thickness where isolated droplets of the other phase get trapped ( γ ˙ t = 19.4 ). Such stretching and bursting persist though this periodic behavior cannot be observed over very long periods of time due to the finite size of the system.
In the strong shear regime, the initial morphology is different with respect to the previous case, since the interface diffusion time is larger than the inverse of the shear rate. Indeed, while domains grow (as quantified later) and start to be deformed by the flow, the two phases are still far from the steady binodal densities (see the panel at γ ˙ t = 2 in Figure 4), an effect promoting an initial thinning of domains along the shear direction. Afterwards, the elongations and ruptures of domains previously described take place once again, a phenomenon also observed in simulations of passive model B with strong shear [41,42]. However, unlike such cases, the main feature here is that, under shear, isolated droplets of the ϕ 1 phase survive and are dispersed in the ϕ 2 matrix (see in particular the last snapshots of Figure 3 and Figure 4). This is not the case if λ = 0 , as demonstrated in Figure 5 where instantaneous configurations of the passive case taken at equal times are shown.
A scenario akin to that just described in the AMB occurs when considering sheared binary mixtures with surface diffusion [64]. However, in the present work the physical explanation is different. The fraction β of the ϕ 1 phase is larger than the one of the ϕ 2 phase when λ > 0 despite the fact that we consider a critical system with a symmetric initial composition such that < ϕ > = 0 (the symbol < > denotes an average over the lattice). Since the dynamics described by Equation (1) is conserved, the average value of the order parameter ϕ is preserved and it has to be β ϕ 1 + ( 1 β ) ϕ 2 = 0 . In a passive mixture the values of the binodals ϕ 1 , ϕ 2 are symmetric with respect to the value ϕ = 0 and it results β = 0.5 . In the AMB the binodals are not symmetric anymore with respect to the value ϕ = 0 and one has β > 0.5 , since | ϕ 1 | < ϕ 2 when λ > 0 . As a result, the ϕ 1 phase is more abundant than the ϕ 2 one. In our study it is ϕ 1 0.71 and ϕ 2 1.15 when λ = 3 so that β 0.62 . Therefore, in AMB, it is the activity that produces an effective off-symmetric mixture though the initial state is not, in agreement with previous studies of AMB [32]. Such off-symmetric mixture consists of droplets of the majority ( ϕ 1 ) phase dispersed in the minority ( ϕ 2 ) one, a situation persisting under shear. We finally note that a different morphology would be observed in the coarsening of off-symmetric passive mixtures under shear with a similar ratio between the two phases [65]. In this case, at low shear rate small droplets evaporate due to Ostwald ripening, while for strong flows it is the minority phase to be dispersed in a large number of droplets.

3.3. Domain Size

The coarsening dynamics can be investigated by looking at evolution of the typical measures of domains. The sizes R x and R y along the flow and the shear directions, respectively, are shown in Figure 6 and are computed as the inverse of the first moments of the structure factor
R x , y ( t ) = π d k C ( k , t ) d k | k x , y | C ( k , t ) .
C ( k , t ) = ϕ ( k , t ) ϕ ( k , t ) is the structure factor averaged over different realizations of the system and ϕ ( k , t ) is the Fourier transform of the order parameter.
At weak shear, it results that R x R y until γ ˙ t 1 , then R y grows initially faster than R x due to the shear-induced dragging of droplets along the flow direction. The typical size along the shear direction attains a maximum at γ ˙ t 1.75 and then decreases, while the growth rate along the flow direction increases due to shear stretching. The quantity R y oscillates on a logarithmic time scale between minimum and maximum values corresponding to elongated and broken domains, respectively. The amplitude of such oscillations diminishes when increasing the activity parameter λ and is related to the existence of small droplets that cannot grow any further. The radii R x for different values of λ show a similar trend compatible with an asymptotic power-law growth R x A ( λ ) t α x . Fitting the power law for R x to numerical data shows that the amplitudes A decrease with the activity as illustrated in Table 1, while the exponent α x 1 does not seem to be affected by λ .
Due to the limited size of the simulated system and the superimposed oscillations, it is difficult to estimate the growth exponent α y along the shear direction although, at γ ˙ t 5 , one observes a regime with α x α y 1 . This result is due to the advection of domains by the applied flow, and has been observed in the passive model B under shear as well [41,42].
At increasing values of shear rate, the growth along the flow direction proceeds faster than that in the shear direction, with a local maximum of R x at γ ˙ t 3 corresponding to a minimum for R y , since domains are strongly deformed by the shear. Note that the height of the maximum of R x is reduced by increasing the activity, a feature associated to the presence of droplets which can be deformed only slightly by the flow. When γ ˙ t 7 , R x shows a power-law growth with an exponent α x 1.1 not depending on λ , while the size R y exhibits a periodic behavior on a logarithmic time scale with amplitude that shrinks with λ . By using power counting, it might be expected that the net effect of the convective term is to increase the growth exponent of the unsheared system by 1 in the flow direction. However, as previously discussed, in AMB with no external flow the value of the growth exponent has not been definitely determined yet. Our results seem to point towards a picture in which the activity has mild effect on the coarsening dynamics, as in the unsheared case [32], since we find that α x 1 as in the passive case and α x α y 1 . Additional simulations run on systems of size L = 512 do not provide further insights for evaluating α x and α y as well as for attempting a finite-size scaling analysis. Indeed, a more accurate estimate of the growth exponents would likely require much larger systems, thus dramatically increasing the computational resources necessary to simulate their late-time dynamics.
To elucidate the role of flow strength, simulations with shear rate γ ˙ = 2.8 × 10 3 are also considered. This value is such that γ ˙ t D ( λ ) 1.1 for 0 λ 3 so that the system is in the strong regime. The time behavior of R x and R y is similar to the one shown in Figure 6, with no significant effect on the growth exponents. We only observe a reduction in the amplitudes of oscillations at the smallest values of activity along the x-direction, since domains are less affected by flow. Being less deformed, they can grow along the shear direction, a result witnessed by wider amplitudes of R y with respect to the case with γ ˙ = 3.6 × 10 3 . These effects reduce when λ increases, since, as previously discussed, the dynamics is deeply modified by the presence of droplets for the highest value of activity.
The size distribution of domains can be analyzed by calculating the normalized probability distributions P ( L x , y ) , having domains of lengths L x and L y along the flow and the shear directions, respectively. By moving along all rows (the flow direction) of the computational domain, L x is computed as the size of unidimensional domains with equal composition (same sign of the density ϕ ). From the registered values of L x , the function P ( L x ) is derived. The same procedure is adopted along the columns (the shear direction) of the lattice to determine L y and then P ( L y ) . In Figure 7 and Figure 8 we plot P ( L x ) and P ( L y ) for λ = 0.3 in the weak and strong shear regime. At weak shear (see Figure 7), P ( L x ) and P ( L y ) show two peaks in the active case and a single one in the passive situation. The peaks at the smaller values of L x and L y (with L x L y 8 ) correspond to the presence of isolated droplets, while the other peaks indicate elongated domains. In particular, P ( L x ) broadens when going from the maximum of R y at γ ˙ t 1.75 to the minimum at γ ˙ t 5.3 and, then, shrinks when R y grows again. In the passive counterpart, the distribution functions are less broad. When looking at the shear direction, in the active case the peak of P ( L y ) at the larger value of L y oscillates in height, with a local maximum attained when R y is at the minimum ( γ ˙ t 5.3 ). Later P ( L y ) broadens. In the passive case the distribution P ( L y ) has a single peak (narrower than in the active case) that oscillates in height, following the cyclic dynamics of elongations and ruptures.
Moving to the strong shear regime (Figure 8), we observe that P ( L x ) is characterized by the presence of two peaks both in the passive and in the active case, roughly at the same positions. When λ = 3 , the peak at L x 10 is always the prevailing one, owing to the presence of droplets spanning the system. The value of L x , corresponding to the second peak, increases in time attaining its maximum at γ ˙ t 17.4 , when domains are highly stretched ( R y is at its minimum). The subsequent breaking-up of domains determines the shift in position of the second peak to smaller values. Moreover, it can be seen that the distributions are broader in the active case. This feature also holds for P ( L y ) , whose main peak oscillates in height during the time evolution. We note here that in the active system the first peak of P ( L y ) at the smaller value of L y , corresponding to isolated droplets, is less pronounced with respect to that of the weak shear regime. This can be attributed to the fact that, as previously discussed, initially forming droplets are fused by shear before growing in size.

4. Conclusions

To summarize, we have investigated the phase separation of an active binary mixture subject to an applied shear flow. For this purpose we numerically solved the phenomenological equation of the active model B [32], supplemented by a convective term that couples the order parameter to the external velocity field. The initial morphology depends on the strength of shear. Later, growing domains are elongated, tilted, and burst by the flow. This is reflected in the typical sizes of domains R x and R y along the two spatial dimensions, which appear to be modulated by oscillations on a logarithmic time scale. However, the presence of activity is such that the fraction of one phase is larger than the other one, despite the initial symmetric composition. This induces the presence of droplets of the more abundant phase which span the system and are responsible for the observed reduction of the amplitudes of the oscillations of R x and R y when increasing the activity parameter. Though the limited size of the simulated system does not allow an estimate of the growth exponents α x and α y along the flow and the shear directions, respectively, we find that α x α y 1 as in the passive case [41,42]. The combined effect of activity and shear on the overall morphology has been studied by considering the probability distribution functions (PDFs) of the size of patterns along the two spatial directions. Our simulations suggest that such PDFs are characterized by a width that broadens with activity.
We hope that our results may stimulate further research on this system. It would be of interest, for example, understanding the role played by the shear in active binary mixtures where all terms breaking time-reversal symmetry to leading order in ∇ and ϕ are included [28], as well as how thermal noise is expected to impact on morphology and growth dynamics in the presence of an external driving.

Author Contributions

A.L. and A.T. designed the research, analysed the data and wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

A.T. acknowledges funding from the European Research Council under the European Union’s Horizon 2020 Framework Programme (No. FP/2014-2020) ERC Grant Agreement No.739964 (COPMAT).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon reasonable request from the corresponding author.

Acknowledgments

The authors wish to thank M. E. Cates for useful comments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Here we detail how the values of ξ and t D in Figure 2 are measured. The interface width ξ is obtained by fitting the steady interface profile along the x-direction with the following functional form
ϕ ( x , λ ) = ϕ 1 ( λ ) + ( ϕ 2 ( λ ) ϕ 1 ( λ ) ) 2 1 + tanh 2 x ξ ,
where ϕ 1 , 2 ( λ ) are the binodals and ξ is the only fit parameter.
The values of t D are evaluated by computing the ratio
Λ ( t ) = i | ϕ i ( t ) ϕ i ( t t 1 ) | i | ϕ i ( t t 1 ) | ,
over time, where ϕ i is the density along a unidimensional profile initialized between the binodals ϕ 1 and ϕ 2 and i is a discrete lattice index running over the whole planar interface ( i = 1 , , L ). This quantity provides an estimate of the relative difference between two profiles taken at a time interval equal to t 1 . The relaxation time t D is the time such that Λ ( t D ) < 10 2 . A reasonable estimate of t D is obtained by using t 1 = 50 . Although other values of t 1 can be adopted, our choice represents a reasonable compromise between the requirement t 1 t D to resolve t D and the need of having a sufficiently large time interval to avoid indistinguishable profiles at times t t 1 and t.

References

  1. Ramaswamy, S. The Mechanics and Statistics of Active Matter. Annu. Rev. Condens. Mat. Phys. 2010, 1, 323–345. [Google Scholar] [CrossRef] [Green Version]
  2. Romanczuk, P.; Bär, M.; Ebeling, W.; Lindner, B.; Schimansky-Geier, L. Active Brownian particles. Eur. Phys. J. Spec. Top. 2012, 202, 1–62. [Google Scholar] [CrossRef] [Green Version]
  3. Marchetti, M.C.; Joanny, J.F.; Ramaswamy, S.; Liverpool, T.B.; Prost, J.; Rao, M.; Sima, R.A. Hydrodynamics of soft active matter. Rev. Mod. Phys. 2013, 85, 1143. [Google Scholar] [CrossRef] [Green Version]
  4. Bechinger, C.; Leonardo, R.D.; Löwen, H.; Reichhardt, C.; Volpe, G.; Volpe, G. Active particles in complex and crowded environments. Rev. Mod. Phys. 2016, 88, 045006. [Google Scholar] [CrossRef]
  5. Ramaswamy, S. Active fluids. Nat. Rev. Phys. 2019, 1, 640. [Google Scholar] [CrossRef]
  6. Leonardo, R.D.; Angelani, L.; Dell’Arciprete, D.; Ruocco, G.; Iebba, V.; Schippa, S.; Conte, M.P.; Mecarini, F.; Angelis, F.D.; Fabrizio, E.D. Bacterial ratchet motors. Proc. Nat. Acad. Sci. USA 2010, 107, 9541. [Google Scholar] [CrossRef] [Green Version]
  7. Sanchez, T.; Chen, D.T.N.; DeCamp, S.J.; Heymann, M.; Dogic, Z. Spontaneous motion in hierarchically assembled active matter. Nature 2012, 491, 431. [Google Scholar] [CrossRef] [Green Version]
  8. López, H.M.; Gachelin, J.; Douarche, C.; Auradou, H.; Clément, E. Turning Bacteria Suspensions into Superfluids. Phys. Rev. Lett. 2015, 115, 028301. [Google Scholar] [CrossRef] [Green Version]
  9. Dell’Arciprete, D.; Blow, M.L.; Brown, A.T.; Farrell, F.D.C.; Lintuvuori, J.S.; McVey, A.S.; Marenduzzo, D.; Poon, W.C.K. A growing bacterial colony in two dimensions as an active nematic. Nat. Commun. 2018, 9, 1. [Google Scholar] [CrossRef] [Green Version]
  10. Theurkauff, I.; Cottin-Bizonne, C.; Palacci, J.; Ybert, C.; Bocquet, L. Dynamic Clustering in Active Colloidal Suspensions with Chemical Signaling. Phys. Rev. Lett. 2012, 108, 268303. [Google Scholar] [CrossRef]
  11. Buttinoni, I.; Bialké, J.; Kümmel, F.; Löwen, H.; Bechinger, C.; Speck, T. Dynamical Clustering and Phase Separation in Suspensions of Self-Propelled Colloidal Particles. Phys. Rev. Lett. 2013, 110, 238301. [Google Scholar] [CrossRef] [Green Version]
  12. Palacci, J.; Sacanna, S.; Steinberg, A.P.; Pine, D.J.; Chaikin, P.M. Living crystals of light-activated colloidal surfers. Science 2013, 339, 936. [Google Scholar] [CrossRef] [PubMed]
  13. Needleman, D.; Dogic, Z. Active particles induce large shape deformations in giant lipid vesicles. Nat. Rev. Mater. 2017, 2, 17048. [Google Scholar] [CrossRef]
  14. Saw, T.B.; Doostmohammadi, A.; Nier, V.; Kocgozlu, L.; Thampi, S.; Toyama, Y.; Marcq, P.; Lim, C.T.; Yeomans, J.M.; Ladoux, B. Topological defects in epithelia govern cell death and extrusion. Nature 2017, 544, 212. [Google Scholar] [CrossRef] [PubMed]
  15. Shaebani, M.R.; Wysocki, A.; Winkler, R.G.; Gompper, G.; Rieger, H. Computational models for active matter. Nat. Rev. Phys. 2020, 2, 181. [Google Scholar] [CrossRef] [Green Version]
  16. Tailleur, J.; Cates, M.E. Statistical Mechanics of Interacting Run-and-Tumble Bacteria. Phys. Rev. Lett. 2008, 100, 218103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Cates, M.E.; Tailleur, J. Motility-Induced Phase Separation. Annu. Rev. Cond. Matt. Phys. 2015, 6, 219–244. [Google Scholar] [CrossRef] [Green Version]
  18. Gonnella, G.; Marenduzzo, D.; Suma, A.; Tiribocchi, A. Motility-induced phase separation and coarsening in active matter. Compt. Ren. Phys. 2015, 16, 316. [Google Scholar] [CrossRef] [Green Version]
  19. Fily, Y.; Marchetti, M.C. Athermal Phase Separation of Self-Propelled Particles with No Alignment. Phys. Rev. Lett. 2012, 108, 235702. [Google Scholar] [CrossRef] [Green Version]
  20. Stenhammar, J.; Tiribocchi, A.; Allen, R.J.; Marenduzzo, D.; Cates, M.E. Continuum Theory of Phase Separation Kinetics for Active Brownian Particles. Phys. Rev. Lett. 2013, 111, 145702. [Google Scholar] [CrossRef] [Green Version]
  21. Stenhammar, J.; Allen, R.J.; Marenduzzo, D.; Cates, M.E. Phase behaviour of active Brownian particles: The role of dimensionality. Soft Matter 2014, 10, 1489. [Google Scholar] [CrossRef] [Green Version]
  22. Speck, T.; Bialké, J.; Menzel, A.M.; Löwen, H. Effective Cahn-Hilliard Equation for the Phase Separation of Active Brownian Particles. Phys. Rev. Lett. 2014, 112, 218304. [Google Scholar] [CrossRef] [Green Version]
  23. Fodor, E.; Nardini, C.; Cates, M.E.; Tailleur, J.; Visco, P.; van Wijland, F. How Far from Equilibrium Is Active Matter? Phys. Rev. Lett. 2016, 117, 038103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Nardini, C.; Fodor, E.; Tjhung, E.; van Wijland, F.; Tailleur, J.; Cates, M.E. Entropy Production in Field Theories without Time-Reversal Symmetry: Quantifying the Non-Equilibrium Character of Active Matter. Phys. Rev. X 2017, 7, 021007. [Google Scholar] [CrossRef]
  25. Takatori, S.C.; Yan, W.; Brady, J.F. Swim Pressure: Stress Generation in Active Matter. Phys. Rev. Lett. 2014, 113, 028103. [Google Scholar] [CrossRef]
  26. Solon, A.P.; Chaté, H.; Tailleur, J. From Phase to Microphase Separation in Flocking Models. Phys. Rev. Lett. 2015, 114, 068101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Ginot, F.; Theuekauff, I.; Levis, D.; Ybert, C.; Bocquet, L.; Berthier, L.; Cottin-Bizonne, C. Nonequilibrium Equation of State in Suspensions of Active Colloids. Phys. Rev. X 2015, 5, 011004. [Google Scholar] [CrossRef]
  28. Tjhung, E.; Nardini, C.; Cates, M.E. Cluster Phases and Bubbly Phase Separation in Active Fluids: Reversal of the Ostwald Process. Phys. Rev. X 2018, 8, 031080. [Google Scholar] [CrossRef] [Green Version]
  29. Caporusso, C.B.; Digregorio, P.; Levis, D.; Cugliandolo, L.F.; Gonnella, G. Motility-Induced Microphase and Macrophase Separation in a Two-Dimensional Active Brownian Particle System. Phys. Rev. Lett. 2020, 125, 178004. [Google Scholar] [CrossRef]
  30. Shi, X.-Q.; Fausti, G.; Chaté, H.; Nardini, C.; Solon, A. Self-Organized Critical Coexistence Phase in Repulsive Active Particles. Phys. Rev. Lett. 2020, 125, 168001. [Google Scholar] [CrossRef]
  31. Cates, M.E. Active Field Theories. arXiv 2019, arXiv:1904.01330. [Google Scholar]
  32. Wittkowski, R.; Tiribocchi, A.; Stenhammar, J.; Allen, R.J.; Marenduzzo, D.; Cates, M.E. Scalar φ4 field theory for active-particle phase separation. Nat. Comm. 2014, 5, 4351. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Hohenberg, P.C.; Halperin, B.I. Theory of dynamic critical phenomena. Rev. Mod. Phys. 1977, 49, 435. [Google Scholar] [CrossRef]
  34. Cahn, J.W.; Hilliard, J.E. Free Energy of a Nonuniform System. I. Interfacial Free Energy. J. Chem. Phys. 1958, 28, 258. [Google Scholar] [CrossRef]
  35. Cahn, J.W.; Hilliard, J.E. Free Energy of a Nonuniform System. III. Nucleation in a Two-Component Incompressible Fluid. J. Chem. Phys. 1959, 31, 688. [Google Scholar] [CrossRef]
  36. Bray, J. Theory of phase-ordering kinetics. Adv. Phys. 1994, 43, 357. [Google Scholar] [CrossRef] [Green Version]
  37. Otha, T.; Nozaki, H.; Doi, M. Computer simulations of domain growth under steady shear flow. J. Chem. Phys. 1990, 93, 2664. [Google Scholar]
  38. Corberi, F.; Gonnella, G.; Lamura, A. Spinodal Decomposition of Binary Mixtures in Uniform Shear Flow. Phys. Rev. Lett. 1998, 81, 3852. [Google Scholar] [CrossRef] [Green Version]
  39. Corberi, F.; Gonnella, G.; Lamura, A. Structure and Rheology of Binary Mixtures in Shear Flow. Phys. Rev. E 2000, 61, 6621. [Google Scholar] [CrossRef] [Green Version]
  40. Berthier, L. Phase separation in a homogeneous shear flow: Morphology, growth laws, and dynamic scaling. Phys. Rev. E 2001, 63, 051503. [Google Scholar] [CrossRef] [Green Version]
  41. Corberi, F.; Gonnella, G.; Lamura, A. Two-Scale Competition in Phase Separation with Shear. Phys. Rev. Lett. 1999, 83, 4057. [Google Scholar] [CrossRef] [Green Version]
  42. Corberi, F.; Gonnella, G.; Lamura, A. Phase Separation of Binary Mixtures in Shear Flow: A Numerical Study. Phys. Rev. E 2000, 62, 8064. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Gonnella, G.; Lamura, A. Long-time behavior and different shear regimes in quenched binary mixtures. Phys. Rev. E 2007, 75, 011501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Rothman, D.H. Complex rheology in a model of a phase-separating fluid. Europhys. Lett. 1991, 14, 337. [Google Scholar] [CrossRef]
  45. Hashimoto, T.; Matsuzaka, K.; Moses, E.; Onuki, A. String Phase in Phase-Separating Fluids under Shear Flow. Phys. Rev. Lett. 1995, 74, 126. [Google Scholar] [CrossRef]
  46. Nakano, H.; Minami, Y.; Sasa, S. Long-Range Phase Order in Two Dimensions under Shear Flow. Phys. Rev. Lett. 2021, 126, 160604. [Google Scholar] [CrossRef]
  47. Saracco, G.P.; Gonnella, G. Critical behavior of the Ising model under strong shear: The conserved case. Phys. A 2021, 576, 126038. [Google Scholar] [CrossRef]
  48. Stenhammar, J.; Wittkowski, R.; Marenduzzo, D.; Cates, M.E. Activity-Induced Phase Separation and Self-Assembly in Mixtures of Active and Passive Particles. Phys. Rev. Lett. 2015, 114, 018301. [Google Scholar] [CrossRef]
  49. Cates, M.E.; Tjhung, E. Theories of binary fluid mixtures: From phase-separation kinetics to active emulsions. J. Fluid Mech. 2018, 836, P1. [Google Scholar] [CrossRef] [Green Version]
  50. Caballero, F.; Nardini, C.; Cates, M.E. From bulk to microphase separation in scalar active matter: A perturbative renormalization group analysis. J. Stat. Mech. 2018, 2018, 123208. [Google Scholar] [CrossRef] [Green Version]
  51. Solon, A.P.; Stenhammar, J.; Cates, M.E.; Kafri, Y.; Tailleur, J. Generalized thermodynamics of phase equilibria in scalar active matter. Phys. Rev. E 2018, 97, 020602R. [Google Scholar] [CrossRef] [Green Version]
  52. Pattanayak, S.; Mishra, S.; Puri, S. Ordering kinetics in the active model B. Phys. Rev. E 2021, 104, 014606. [Google Scholar] [CrossRef] [PubMed]
  53. Pattanayak, S.; Mishra, S.; Puri, S. Domain Growth in the Active Model B: Critical and Off-critical Composition. Soft Mater. 2021, 19, 286. [Google Scholar] [CrossRef]
  54. Redner, G.S.; Hagan, M.F.; Baskaran, A. Structure and Dynamics of a Phase-Separating Active Colloidal Fluid. Phys. Rev. Lett. 2013, 110, 055701. [Google Scholar] [CrossRef]
  55. Tiribocchi, A.; Wittkowski, R.; Marenduzzo, D.; Cates, M.E. Scalar Active Matter in a Momentum-Conserving Fluid. Phys. Rev. Lett. 2015, 115, 188302. [Google Scholar] [CrossRef] [Green Version]
  56. Kardar, M.; Parisi, G.; Zhang, Y.-C. Dynamic Scaling of Growing Interfaces. Phys. Rev. Lett. 1986, 56, 889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Hunter, J.K.; Saxton, R. Dynamics of Director Fields. Siam J. Appl. Math. 1991, 51, 1498. [Google Scholar] [CrossRef]
  58. Kuramoto, Y. Diffusion-Induced Chaos in Reaction Systems. Prog. Theor. Phys. Suppl. 1978, 64, 346. [Google Scholar] [CrossRef]
  59. Sivashinsky, G.I. On Flame Propagation Under Conditions of Stoichiometry. Siam J. Appl. Math. 1980, 39, 67. [Google Scholar] [CrossRef]
  60. Rogers, T.M.; Elder, K.R.; Desai, R.C. Numerical study of the late stages of spinodal decomposition. Phys. Rev. B 1988, 37, 9638. [Google Scholar] [CrossRef]
  61. Strikwerda, J.C. Finite Difference Schemes and Partial Differential Equations; Chapman and Hall: New York, NY, USA, 1989. [Google Scholar]
  62. Lees, A.W.; Edwards, S.F. The computer study of transport processes under extreme conditions. J. Phys. C 1972, 5, 1921. [Google Scholar] [CrossRef]
  63. Frischknecht, A. Effect of shear flow on the stability of domains in two-dimensional phase-separating binary fluids. Phys. Rev. E 1997, 56, 6970. [Google Scholar] [CrossRef] [Green Version]
  64. Gonnella, G.; Lamura, A. Sheared phase-separating binary mixtures with surface diffusion. J. Phys. A 2020, 53, 305002. [Google Scholar] [CrossRef]
  65. Lamura, A.; Gonnella, G. Computer simulations of domain growth in off-critical quenches of two-dimensional binary mixtures. Comput. Mater. 2002, 25, 531. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Steady density profiles for different values of the parameter λ .
Figure 1. Steady density profiles for different values of the parameter λ .
Mathematics 09 03008 g001
Figure 2. Interface width ξ (left panel) and diffusion time t D (right panel) as functions of the activity λ .
Figure 2. Interface width ξ (left panel) and diffusion time t D (right panel) as functions of the activity λ .
Mathematics 09 03008 g002
Figure 3. Configurations of the system at consecutive times in the weak shear regime ( γ ˙ t D = 1.0 ) with λ = 3 . A central portion of size 512 × 512 of the whole lattice is shown.
Figure 3. Configurations of the system at consecutive times in the weak shear regime ( γ ˙ t D = 1.0 ) with λ = 3 . A central portion of size 512 × 512 of the whole lattice is shown.
Mathematics 09 03008 g003
Figure 4. Configurations of the system at consecutive times in the strong shear regime ( γ ˙ t D = 4.0 ) with λ = 3 . A central portion of size 512 × 512 of the whole lattice is shown.
Figure 4. Configurations of the system at consecutive times in the strong shear regime ( γ ˙ t D = 4.0 ) with λ = 3 . A central portion of size 512 × 512 of the whole lattice is shown.
Mathematics 09 03008 g004
Figure 5. Configurations of the system in the weak (left panel, γ ˙ t D = 0.4 ) and strong (right panel, γ ˙ t D = 1.4 ) shear regimes with λ = 0 . Snapshots are taken at the latest times of Figure 3 and Figure 4. A central portion of size 512 × 512 of the whole lattice is shown.
Figure 5. Configurations of the system in the weak (left panel, γ ˙ t D = 0.4 ) and strong (right panel, γ ˙ t D = 1.4 ) shear regimes with λ = 0 . Snapshots are taken at the latest times of Figure 3 and Figure 4. A central portion of size 512 × 512 of the whole lattice is shown.
Mathematics 09 03008 g005
Figure 6. Average measures of domains along the flow ( R x ) and the shear ( R y ) directions as functions of time at weak (left panel) and strong (right panel) shear for different values of λ .
Figure 6. Average measures of domains along the flow ( R x ) and the shear ( R y ) directions as functions of time at weak (left panel) and strong (right panel) shear for different values of λ .
Mathematics 09 03008 g006
Figure 7. Probability distribution functions P of domains of length L x (left panels) and L y (right panels) with λ = 3 (upper row) and λ = 0 (lower row) in the weak shear regime.
Figure 7. Probability distribution functions P of domains of length L x (left panels) and L y (right panels) with λ = 3 (upper row) and λ = 0 (lower row) in the weak shear regime.
Mathematics 09 03008 g007
Figure 8. Probability distribution functions P of domains of length L x (left panels) and L y (right panels) with λ = 3 (upper row) and λ = 0 (lower row) in the strong shear regime.
Figure 8. Probability distribution functions P of domains of length L x (left panels) and L y (right panels) with λ = 3 (upper row) and λ = 0 (lower row) in the strong shear regime.
Mathematics 09 03008 g008
Table 1. Fitted values of the amplitudes A for different values of the activity λ in the weak shear regime.
Table 1. Fitted values of the amplitudes A for different values of the activity λ in the weak shear regime.
λ A
04.40
12.75
20.99
30.82
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lamura, A.; Tiribocchi, A. Shearing Effects on the Phase Coarsening of Binary Mixtures Using the Active Model B. Mathematics 2021, 9, 3008. https://doi.org/10.3390/math9233008

AMA Style

Lamura A, Tiribocchi A. Shearing Effects on the Phase Coarsening of Binary Mixtures Using the Active Model B. Mathematics. 2021; 9(23):3008. https://doi.org/10.3390/math9233008

Chicago/Turabian Style

Lamura, Antonio, and Adriano Tiribocchi. 2021. "Shearing Effects on the Phase Coarsening of Binary Mixtures Using the Active Model B" Mathematics 9, no. 23: 3008. https://doi.org/10.3390/math9233008

APA Style

Lamura, A., & Tiribocchi, A. (2021). Shearing Effects on the Phase Coarsening of Binary Mixtures Using the Active Model B. Mathematics, 9(23), 3008. https://doi.org/10.3390/math9233008

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop