## Flow-induced nonequilibrium self-assembly in suspensions of stiff, apolar, active filaments

Ankita Pandey, 1 P. B. Sunil Kumar, 1 and R. Adhikari 2

1 Indian Institute of Technology, Madras, Chennai 600036, India 2 The Institute of Mathematical Sciences, CIT Campus, Chennai 600113, India (Dated: August 20, 2015)

Active bodies in viscous fluids interact hydrodynamically through self-generated flows. Here we study spontaneous aggregation induced by hydrodynamic flow in a suspension of stiff, apolar, active filaments. Lateral hydrodynamic attractions in extensile filaments lead, independent of volume fraction, to anisotropic aggregates which translate and rotate ballistically. Lateral hydrodynamic repulsion in contractile filaments lead, with increasing volume fractions, to microstructured states of asters, clusters, and incipient gels where, in each case, filament motion is diffusive. Our results demonstrate that the interplay of active hydrodynamic flows and anisotropic excluded volume interactions provides a generic nonequilibrium mechanism for hierarchical self-assembly of active soft matter.

Filamentous structures along which chemical energy is converted to mechanical motion are found in many biological contexts such as actin filaments, molecular motors walking on microtubules, flagellar bundles and ciliary hairs. These cellular components maintain structure, influence signaling, and provide motile forces. Much recent effort has been directed at synthesising biomimetic analogues of such components. Motor-microtubule mixtures [1] and self-assembled motor-microtubules bundles capable of motility [2] are two examples where dynamic structures remarkably similar to those occurring in vivo have been synthesized in the laboratory. These emergent structures are maintained by nonequilibrium 'active' forces generated by the consumption of chemical energy. Hydrodynamic flow provides a way of transmitting nonequilibrium forces and its role in creating and maintaining dynamic structures remains relatively unexplored.

Our present theoretical understanding of hydrodynamics in nonequilibrium energy-converting systems is based largely on continuum [3] or kinetic [4] theories that prescribe, respectively, the spatiotemporal evolution of coarse-grained fields or distribution functions. These theories predict long wavelength collective phenomena like instabilities [5, 6], the onset of flows [7] through spontaneous symmetry breaking [8] and the continuous generation of topological defects [9]. The complementary domain of short wavelength phenomena where near-field hydrodynamics, the shape of the active body and excluded volume interactions are important, has been much less studied [10, 11]. This domain is relevant in situations where particle numbers are not thermodynamically large and a hydrodynamic limit is not necessarily attained, as for example in the cell. Theoretical descriptions that do not presume the existence of continuum and hydrodynamic limits, then, are required to provide insight into the physics of short wavelength phenomena in active matter.

In this Letter we study short wavelength aggregation phenomena in suspensions of stiff, apolar, active filaments using a description that resolves individual filaments and includes their hydrodynamic and excluded volume interactions. Our model for an active filament consists of rigidly connected active beads that produce force-free, torque-free flows in the fluid [10]. We approximate the active flow of each bead by its slowest decaying component, the irreducible stresslet [12]. The stresslet principal axis is always aligned along the filament tangent and the principal values are either positive (extensile) or negative (contractile). This ensures that an isolated straight filament produces active flows but does not actively translate or rotate and is, thus, individually apolar. However, as we show, translational and rotational motion is possible in the combined flow of two or more stiff, apolar filaments. This inherently collective and nonequilibrium source of motion, when constrained by the anisotropic excluded volume interactions between filaments, leads to non-trivial states of aggregation in suspension.

The apolar filament suspension is an example of a system where forces that break time reversal invariance, and which are usually associated with dissipation and the cessation of motion, drive the system and lead to nonquilibrium steady states (NESS) with non-zero fluxes. A wellknown paradigm of NESS is the driven dissipative system, where external forces at the system boundary produce nonequlibrium. In contrast, here internal forces distributed continuously within the system boundary produce nonequilibrium. Therefore, apolar active suspensions are novel model systems for studying NESS both theoretically and experimentally. Many recent developments in the statistical physics of NESS, such as fluctuation theorems, macroscopic fluctuation theory, and motility-induced phase separation can be tested in the system we study here.

Our results are as follows. We find distinct fixed points for the non-linear dynamics of pairs of filaments translating and rotating in their mutual flow. Extensile fila-

FIG. 1. (Color online) Self-assembly of extensile filaments. (a) Flow field of two extensile filaments producing lateral attraction. (b) Lateral attraction drive rapid dimerization at early times. (c) Coagulation of k -mers at intermediate times. (d) Aggregation of k -mers into a single cluster at late times. The background color indicates the magnitude of the velocity normalized by its maximum. Snapshot b-d are for a φ = 0 05 system taken at . t = 0 1 . τ e , t = 0 2 . τ e and t = 1 8 . τ e respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0433v2-figure-414e1ae0-002d-4999-af48-186efd5ec10a.jpg)

ments, driven by dominant lateral attractions, join into stable rotating dimers. Contractile filaments, driven by dominant perpendicular attractions, join into translating and rotating T-shaped dimers. Extensile filament suspensions are always unstable and show a two-step aggregation kinetics. A rapid initial growth of small clusters is followed by a slower merging of these clusters into large, dynamic aggregates which translate and rotate ballistically. Contractile filament suspensions remain homogeneous and isotropic but develop short wavelength microstructure that depends on volume fraction. We observe such microstructures in a sequence of states of asters, clusters, and incipient gels with increasing volume fraction. The centres of mass of the filaments show a non-thermal diffusion which decreases with increasing volume fraction. We describe these results in detail below and then discuss their significance.

Active Zimm model: We briefly describe our model for an active filament [10]. We consider N b spherical active beads of radius a centred at r n . These are acted on by elastic forces F n = -∂U ( r 1 , . . . , r N b ) /∂ r n derived from the potential U ( r 1 , . . . , r N b ) which bonds neighbouring beads into a connected filament, penalises filament curvature, and enforces volume exclusion between beads. The equilibrium bond length is b . Each bead has a stresslet S n = S 0 ( t n t n -1 3 I ) aligned along the local tangent t n . The beads move in response to the elastic forces acting on them and the hydrodynamic flow which results, at low Reynolds numbers and small bead size, in the following equation of motion glyph[negationslash]

$$\dot { \mathbf r } _ { n } = \mu \mathbf F _ { n } + \sum _ { m \neq n } \left [ G ( \mathbf r _ { n } - \mathbf r _ { m } ) \cdot \mathbf F _ { m } + \nabla \mathbf G ( \mathbf r _ { n } - \mathbf r _ { m } ) \cdot \mathbf S _ { m } \right ]. \\ \mathbf r _ { \dots \cdots } \cdot \mathbf C / \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ \mathbf r _ { \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \ddots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots.$$

Here, G r ( ) is a Green function of the Stokes equation, µ is the bead mobility and S 0 is the principal value of the stresslet. The flow produced by a single bead is extensile or contractile when S 0 is, respectively, positive or negative. The stresslet self-term is absent as it cannot produce translational motion in an unbounded domain. These equations of motions may be thought of as a generalization, to active filaments, of the Zimm model of polymer dynamics where hydrodynamic interactions are treated at the level of the Kirkwood-Risemann superposition approximation. The accuracy of this approximation is ( a/b ) 3 for monopolar flows and, by an extension of the original study, ( a/b ) 4 for dipolar flows. With the choice of parameters made here, we make no more than 6% error in computing hydrodynamic interactions using the superposition approximation. We simulate a suspension of N such filaments using the lattice Boltzmann method to compute the active flows, as described in detail in [10].

Activity induced dimerization : As a prelude to collective dynamics we study the motion of a pair of ac-

FIG. 2. (Color online) Growth kinetics in extensile filaments. The mean of the k -mer distribution function 〈 k 〉 as a function of time t , for varying filament numbers N . Dimerization dominates the rapid growth for t glyph[lessmuch] τ e while coagulation of k -mers dominates the slow growth for t glyph[greatermuch] τ e . These regimes correspond, respectively, to panels (b) and (c) of Fig.(1).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0433v2-figure-c6ea55e7-edfc-4b9d-ab7f-08f9d1ec4d4c.jpg)

FIG. 3. (Color online) Self-assembly of contractile filaments. (a) Flow field of two contractile filaments producing perpendicular attraction. Microstructured steady states of asters at volume fraction φ = 0 11 in (b), clusters at . φ = 0 15 in (c) and incipient . gel at φ = 0 34 in (d). . A typical aster configuration is circled in (b). The background color indicates the magnitude of the velocity normalized by its maximum.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0433v2-figure-b353e14f-4010-47b1-bce6-f7b528f50ed2.jpg)

tive filaments in their mutual hydrodynamic flow. The flow points outwards along the long-axis for extensile filaments but inwards for contractile filaments. More generally, contractile flow is obtained from extensile flow by reversing all streamlines [10]. Extensile filament pairs attract laterally as is clear from the streamlines around a proximate pair shown in Fig.(1a). Extensile filament pairs repel perpendicularly as is clear from the reversed streamlines of Fig.(3a). The lateral attraction drives parallel filaments towards each other till steric forces prevent further motion. The shear flow along the filament long axis causes a small lateral motion, thus forming a stable rotating dimer (see movie in [13]). This is the stable state for a pair of extensile filaments and is independent of the special initial condition chosen in the simulation. Streamlines around a contractile filament flow in a direction opposite to that of extensile filaments. Thus, contractile filament pairs repel laterally and attract perpendicularly in their mutual flow as is clear from the streamlines in Fig.(1a) and Fig.(3a). The perpendicular attraction, combined with steric repulsion, drives contractile filaments to form T -shaped dimers which translate and rotate (see movie in [13]). Again, this stable state is independent of the initial condition used to generate it.

The dimerization described above is an essentially nonequilibrium phenomenon mediated by the active hydrodynamic flow. Although the microscopic processes that give rise to activity-induced flows and result in the dimerization phenomenon might vary greatly in origin and function, the resultant hydrodynamic flows and interactions can be completely parametrised using the irreducible expansion of the active flow [12], the leading term of which is retained in this analysis. This makes the nonequilibrium dimerisation both robust and generic.

Aggregation of extensile filaments : The nonequilibrium flow-induced forces described above produce selfassembled structures in suspensions of filaments which depend crucially on the sign of the stresslet. In exten- sile filaments, the lateral hydrodynamic attraction destabilises the homogeneous isotropic state of the suspension. Dimer pairs form rapidly and further aggregate to form k -mers. Asymptotically, the individual clusters collapse into a single dynamic cluster which translates and/or rotates according to its final shape. The qualitative features of the aggregation kinetics is shown in the panels (b-d) of Fig.1 and in the movie in [13].

To quantitatively understand the aggregation kinetics we compute the distribution P k, t ( ) of k -mers, which we define as k filaments that are within half a filament length of each other. The distribution is computed over an ensemble of 50 simulations with different initial conditions. The first moment 〈 k 〉 ( ) t = ∑ ∞ k =1 k P ( k, t ) provides a measure of the mean size of the aggregates and this quantity is plotted in Fig.(2). There is a rapid initial growth of the mean size for a duration τ e followed by a slower growth at later times. We estimate τ e from the intersection of a pair of straight lines that best fit the data. From this, we find that the duration τ e of rapid growth is weakly dependent or, possibly, independent of volume fraction while the duration of slow growth to saturation increases with volume fraction.

The nonequilibrium aggregation proceeds through a frequent coagulation process in which k -mers join j -mers to form ( k + )-mers and a rare fragmentation process j in which ( k + )-mers break into j j -mers and k -mers. A careful observation of the kinetics reveals a hierarchy of rates: dimerization is the most rapid process, followed by single filaments joining k -mers (with k ≥ 3), then dimers joining k -mers and, finally, k -mers joining j -mers (with j ≥ 3). The fragmentation rates show an opposite hierarchy: the breaking of a ( k + )-mer into j k -mers and j -mers with k ≈ j is the most rapid process, followed the breaking of dimers and, finally, monomers from ( k + )-mers. These rates are determined by the hydrodyj namic flow due to k -mers and j -mers and is sensitive to their shape. A direct calculation of the coagulation and

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0433v2-figure-d97d8250-f6f7-4c21-83c3-23e8b74a7f7e.jpg)

(a) Stationary cluster size distribution

(b) Filament center of mass diffusion

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0433v2-figure-9506c030-16fa-40e6-b3f4-26edb6341f51.jpg)

FIG. 4. (Color online) (a) Steady state k -mer distribution function P glyph[star] ( k ) with varying volume fraction φ . The distribution shows two qualitative changes at φ = 0 15 and . φ = 0 34. . (b) Filament diffusion coefficient D as a function of volume fraction φ and mean square displacement (MSD) as a function of time t (inset). D is non-dimensionalised by L /τ 2 c , where L = ( N b -1) b is the filament length and τ c the typical time taken to reach stationarity from initially ordered states.

fragmentation rates is, therefore, a formidable problem. Nonetheless, constructing a tractable theoretical model for the rates of hydrodynamic coagulation and fragmentation is desirable as the rates can then be employed in a Smoluchowski-like mean-field description of aggregation kinetics.

Microstructure of contractile filaments: Nonequilib- rium self-assembly in contractile filaments proceeds in a manner very different from extensile filaments. Here, perpendicular hydrodynamic attraction induces dynamic microstructure in initially homogeneous and isotropic states of the suspension. Averaged over long times and large lengths, the suspension remains homogeneous and isotropic, but any instantaneous configuration reveals both positional inhomogeneities and orientational anisotropies. Three distinct microstructures emerge as a function of volume fraction, as shown in panels (bd) of Fig.(3). At the lowest volume fractions, filaments transiently cluster into star-shaped k -mers reminscent of asters (panel b). These k -mers frequently break and reform while their centers of mass diffuse as shown in the movie in [13]. Through this dynamics, the spatiotemporally averaged position and orientation become, respectively, uniform and isotropic. With an increase in volume fraction, the steric constraints limit the dynamic joining and breaking of k -mers and their diffusion (panel c). A percolating network of filaments is formed which, at short time scales, may present a solid-like response by being able to support forces and torques. With further increase in volume fraction, this network gets sterically constrained to a great extent and thus shows greatly reduced positional and orientational dynamics. States that are initially inhomogeneous and have a preferred orientation are unstable and evolve, rapidly, into homogeneous, isotropic microstructured states shown in the movie [13].

In order to distinguish between the distinct microstructural states, we compute the stationary k -mer distribution P glyph[star] ( k ) = lim t →∞ P k, t ( ) for different filament volume fractions and plot it in Fig. (4a). Asters correspond to φ = 0 05 0 11 0 15, clusters to . , . , . φ = 0 21 0 24 and the . , . incipient gel to φ = 0 34. . The distribution shows qualitatively different shapes for each microstructure and the rapid change in shape around φ = 0 3 provides a signa-. ture of gelation. The motion of individual filaments is diffusive in all three microstructured states as the linear variation of the mean square displacement with time shows in the inset of Fig.(4b). The diffusive constant D decreases with volume fraction as Fig.(4b) shows. The parallel hydrodynamic repulsion between filaments enhances their effective volume which increases the steric hindrance and so constrains filament motion, leading to smaller diffusion. It is important to emphasise that the diffusive motion occurs in the absence of thermal fluctuations and consequently D does not obey a Stokes-Einstein relation.

Discussion: Hierarchical self-assembly in thermodynamic equilibrium is driven by a competition between energy and entropy. In contrast, the self-assembly mechanism reported here is an essentially nonequilibrium phenomenon driven by dissipative, hydrodynamic flow produced by active processes. The relative translation of pairs of particles, which contributes to microstructure formation, is possible to due the absence of Onsager sym-

metry of the mobility matrix that prevents such motion in passive hydrodynamic flow. This short wavelength physics is resolved by our particle picture of active flow but is beyond the reach of current long wavelength continuum and kinetic theories.

Several recent experiments report dynamic selfassembly of anisotropic particles in hydrodynamic flow. The formation of 'living' colloidal crystals on planar substrates [14] bears a striking similarity with the aggregation of extensile filaments into liquid crystalline clusters. We suggest that extensile, dipolar hydrodynamic flows, modified by the presence of a no-slip boundary to induce attraction, play a crucial role in the self-assembly of such 'living' crystals. The dimerization of polar bi-metallic catalytic nanomotors [15] is, again, strikingly similar to the dimerization we report here. To directly verify some of our predictions we suggest experiments on apolar active particles made by adding a third metallic stripe of gold or platinum to bi-metallic catalytic nanomotors. Finally, we note that our work suggest a new nonequilibrium mechanism with which to control and manipulate hierarchical self-assembly of soft matter.

Computing resources through HPCE, IIT Madras and Annapurna, IMSc are gratefully acknowledged. RA thanks the Indo-US Science and Technology Foundation for a fellowship. The authors thank M. E. Cates, P. Chaikin, S. Ghose, A. Laskar, T. E. Mallouk, H. Masoud, D. Pine, M. Shelley and R. Singh for helpful discussions, S. Ghose for help with figures and R. Manna for sharing simulation data on filament pairs.

- [1] F. J. N´d´ elec, e T. Surrey, A. C. Maggs, and S. Leibler, Nature 389 , 305 (1997).
- [2] T. Sanchez, D. Chen, S. DeCamp, M. Heymann, and Z. Dogic, Nature 491 , 431 (2012).
- [3] M. C. Marchetti, J. F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha, Rev. Mod. Phys. 85 , 1143 (2013).
- [4] D. Saintillan and M. J. Shelley, Comptes Rendus Physique 14 , 497 (2013).
- [5] R. A. Simha and S. Ramaswamy, Phys. Rev. Lett. 89 , 058101 (2002).
- [6] D. Saintillan and M. J. Shelley, Phys. Rev. Lett. 100 , 178103 (2008).
- [7] F. G. Woodhouse and R. E. Goldstein, Physical review letters 109 , 168105 (2012).
- [8] E. Tjhung, D. Marenduzzo, and M. E. Cates, Proceedings of the National Academy of Sciences 109 , 12381 (2012).
- [9] L. Giomi, M. J. Bowick, X. Ma, and M. C. Marchetti, Physical review letters 110 , 228101 (2013).
- [10] G. Jayaraman, S. Ramachandran, S. Ghose, A. Laskar, M. Saad Bhamla, P. B. Sunil Kumar, and R. Adhikari, Phys. Rev. Lett. 109 , 158302 (2012); A. Laskar, R. Singh, S. Ghose, G. Jayaraman, P. S. Kumar, and R. Adhikari, Scientific reports 3 (2013), 10.1038/srep01964.
- [11] R. Chelakkot, A. Gopinath, L. Mahadevan, and M. F. Hagan, Journal of The Royal Society Interface 11 , 20130884 (2014); H. Jiang and Z. Hou, Soft Matter 10 , 1012 (2014).
- [12] S. Ghose and R. Adhikari, Phys. Rev. Lett. 112 , 118102 (2014).
- [13] See Supplemental Material on YouTube for movies showing aggregation kinetics and microstructure formation and attached Supplementary Information for theoretical estimates of the flow-induced motion, and suggested designs of apolar nanorods.
- [14] J. Palacci, S. Sacanna, A. P. Steinberg, D. J. Pine, and P. M. Chaikin, Science 339 , 936 (2013).
- [15] W. Wang, W. Duan, A. Sen, and T. E. Mallouk, Proc. Natl. Acad. Sci. 110 , 17744 (2013).
Ankita Pandey P. B. Sunil Kumar R. Adhikari
Flow-induced nonequilibrium self-assembly in suspensions of stiff, apolar, active filaments
Active bodies in viscous fluids interact hydrodynamically through
self-generated flows. Here we study spontaneous aggregation induced by
hydrodynamic flow in a suspension of stiff, apolar, active filaments. Lateral
hydrodynamic attractions in extensile filaments lead, independent of volume
fraction, to anisotropic aggregates which translate and rotate ballistically.
Lateral hydrodynamic repulsion in contractile filaments lead, with increasing
volume fractions, to microstructured states of asters, clusters, and incipient
gels where, in each case, filament motion is diffusive. Our results demonstrate
that the interplay of active hydrodynamic flows and anisotropic excluded volume
interactions provides a generic nonequilibrium mechanism for hierarchical
self-assembly of active soft matter.