Physical Review E 90, 013019 (2014)

Casimir Effect in Active Matter Systems

D. Ray1,2, C. Reichhardt1, and C. J. Olson Reichhardt1

1Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
2Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556, USA

(Received 11 March 2014; published 23 July 2014)

We numerically examine run-and-tumble active matter particles in Casimir geometries composed of two finite parallel walls. We find that there is an attractive force between the two walls of a magnitude that increases with increasing run length. The attraction exhibits an unusual exponential dependence on the wall separation, and it arises due to a depletion of swimmers in the region between the walls by a combination of the motion of the particles along the walls and a geometric shadowing effect. This attraction is robust as long as the wall length is comparable to or smaller than the swimmer run length, and is only slightly reduced by the inclusion of steric interactions between swimmers. We also examine other geometries and find regimes in which there is a crossover from attraction to repulsion between the walls as a function of wall separation and run length.
I. INTRODUCTION
II. SIMULATION
III. RESULTS
A. Attractive force on walls
B. Reverse Casimir effect
C. Infinite walls
IV. SUMMARY
References

I.  INTRODUCTION

The Casimir geometry consists of two finite parallel plates placed at a fixed distance from each other that experience an attraction due to the confinement of fluctuations in the media between the plates [1]. In the original Casimir effect calculation, electromagnetic vacuum fluctuations produce an attractive force between two metal places in a vacuum. Much later, Casimir forces were experimentally measured [2], and they have recently been studied in a variety of systems in order to understand how to control their magnitude [3] or polarity [4]. Confined classical fluctuations can produce the so-called critical Casimir effect [5,6,7], first proposed by Fisher and de Gennes near critical demixing in bulk mixtures. The critical Casimir effect has been directly measured [8] and studied in various colloidal systems where it can produce colloidal aggregation [8,9,10]. Casimir type effects have also been studied in granular media, where attractive forces arise between objects or plates placed in vibrated or flowing sand [11]. Critical Casimir effects have also been proposed to occur near percolation thresholds [12] and in biological systems such as near fluctuating cellular membranes [13]. The ability to enhance or control such forces can lead to a wide variety of applications in self-assembly, particle transport, and the creation of novel devices.
Strong fluctuations appear in active matter or self-driven particle systems [14] such as swimming bacteria undergoing run-and-tumble dynamics [15,16]. Recently a number of non-biological active matter systems have been realized experimentally, including artificial swimmers [17], self-driven colloids [18,19,20], or light-activated colloidal particles performing a directed random walk [19,20,21]. Particles undergoing run-and-tumble or active Brownian motion can phase separate at large run length or high density, forming dense regions separated by a dilute active gas [22,23,24,21]. Monodisperse active particles, such as self-propelled disks, can form intermittent dense patches with crystalline order, termed "living crystals" [20]. It has been shown that run-and-tumble dynamics and active Brownian motion can be mapped on to each other when the motility parameters depend on the particle density, so that results obtained with one class of system should be generalizable to the other [25].
Here we demonstrate that a Casimir-like attractive force arises between two plates placed in a bath of active particles and can be controlled by changing the plate geometry. There is already some evidence that active matter can induce forces on objects. One example is the ratchet effect observed for swimming bacteria and run-and-tumble particles in an array of asymmetric funnels in dilute or semi-dense regimes [16,26,27,28,29]. Here, when the active particles run along the funnel walls they can escape through the easy direction of the funnel or become trapped in a funnel tip. Their continued swimming produces forces on the funnel walls they contact. Such induced forces were more clearly demonstrated in systems where an untethered asymmetric sawtooth gear rotates in a preferred direction when placed in a bacterial bath [30,31], but only if the bacteria are actively swimming. In simulations, asymmetric objects placed in an active matter bath are pushed to produce an active matter-driven micro-shuttle [32,33].
In the Casimir geometry, there are no sharp corners in which the particles can accumulate [36]; however, we show that confinement effects alone are sufficient to create fluctuation-induced forces. The magnitude of the attractive force increases with increasing run length and can vary over several orders of magnitude. For fixed run length, the force F between the plates as a function of plate spacing d obeys F(d) ∝ exp(−d/λ), rather than a power law which is usually observed in Casimir geometries [1,2,4,6]. The attractive force arises due to the depletion of the particle density between the plates by a combination of the motion of the particles along the walls and geometric shadowing, both of which increase with increasing run length. It is not an excluded volume effect or depletion interaction of the types studied in Refs. [34,10], since in our model the particles and walls are of vanishing width. We also examine other geometries and show that it is possible to control the magnitude of the forces and induce a crossover from an attractive to a repulsive force. For the case of infinite walls, we analytically derive the force as a function of the distance between the walls.
Fig1.png
Figure 1: A schematic of the system containing run-and-tumble particles (balls) with some particle trajectories indicated by lines and arrows. The run length is Rl. The two parallel walls (bars) of length l are separated by a distance d. When a particle moves along a wall it imparts a force against the wall. We measure the net force per particle ∆F/N=(|Fout|−|Fin|)/N.

II.  SIMULATION

In Fig. 1 we show a schematic of our system which resides in a two-dimensional (2D) simulation box of size Lx ×Ly with periodic boundary conditions in the x- and y-directions. Unless otherwise noted, we take Lx=120 and Ly=60. Within the box are N active particles and two parallel walls of length l separated by a distance d. The run-and-tumble particles move in a fixed randomly chosen direction during a running time τ before undergoing a tumbling event and running in a new randomly chosen direction. Since we consider the dilute limit where particle-particle interactions are rare, for efficiency we use event-driven (ED) dynamics simulations and neglect particle-particle interactions. The particles run with speed v=1 and travel a run length of Rlv between tumbles. A particle contacting a wall at an angle θ moves with velocity sinθ along the wall and exerts a force cosθ on the wall. A particle that reaches the end of the wall before tumbling resumes its original swimming direction. The assumption that run-and-tumble particles move along or accumulate at walls was confirmed in ratchet geometry [16] and asymmetric gear experiments [30,31] and observed in simulations [26,27,28]. We measure the time-averaged forces Fout (Fin) imparted on the outer (inner) surfaces of the walls. The net force per particle acting to bring the walls together or apart is ∆F/N=(|Fout|−|Fin|)/N. We confirm our results using molecular dynamics (MD) simulations with and without particle-particle interactions [26], where we implement the steric interactions via a short-range harmonic repulsion between the particles which have a radius of rd = 0.25. We find that all of our results are robust for the particle densities we consider.
Fig2.png
Figure 2: (a) ∆F/N, the attractive force between the walls, vs d, the inter-wall distance, from ED simulations with l=20, for various Rl, showing an exponential behavior. (b) ∆F/N vs Rl for the same system at various d. For large Rl, ∆F/N saturates. Inset: ∆F/N vs d for l=20, Rl=40, and N=400 in ED simulations (•) and MD simulations without (+) and with (⧫) steric particle-particle interactions.

III.  RESULTS

A.  Attractive Force on Walls

In Fig. 2(a) we plot the normalized force ∆F/N versus d for ED simulations of a system with l=20 at various run lengths. ∆F/N increases with increasing Rl, and we find ∆F(d) ∝ Aexp(−d/λ), with λ increasing for increasing Rl. At small run lengths such as Rl=1, ∆F/N becomes very small, and in the Brownian limit of infinitesimal Rl, ∆F/N=0. In Fig. 2(b) we plot ∆F versus Rl for the same system at various d, showing that the attractive force increases with increasing Rl and saturates at large Rl.
To verify the robustness of the results, we perform three different types of simulations highlighted in the inset of Fig. 2(b) where we plot ∆F/N vs d for systems with N=400, l=20, and Rl=40. We find nearly identical results for ED simulations and MD simulations with and without steric interactions, indicating that our results are robust for interacting particles in the dilute limit. We expect that at higher densities, phase separation or clustering will become important if steric particle-particle interactions are included [22,23,24].
Fig3.png
Figure 3: ρv(r), the rate at which locations are visited by particles, for the system in Fig. 2(a) with Lx=60, l=20, and d=10 at Rl=1 (a), 5 (b), 20 (c), 50 (d). ρv between the plates drops with increasing Rl, giving a larger attractive force.
In order to better understand the origin of the attractive force, in Fig. 3 we plot ρv(r), the rate at which individual spatial locations are visited by particles, from ED simulations with Lx=60, l=20, and d=10 for different values of Rl. We use the normalization that ρv=1.0 in a homogeneous system without walls. In Fig. 3(a), for short run lengths Rl=1, ρv ≈ 1 everywhere and ∆F/N ≈ 0. At larger Rl in Figs. 3(b,c), two trends appear. Firstly, particles near the walls strike the walls and run along them rather than being reflected, raising ρv along the wall surfaces and diminishing ρv nearby; since these particles increasingly reach the end of the wall before tumbling as Rl/l increases, the regions at the ends of the walls are traversed very frequently by particles, forming spots of high ρv. Secondly, a geometric shadowing effect depresses ρv in the inner region between the walls: this region can only be replenished by bulk particles that enter the narrow gap between walls. This will occur with diminishing frequency as the gap size d decreases. In contrast, the region just outside the walls is replenished from the bulk, whose population is continually refreshed by particles reaching the end of either side of a wall. Since the inner and outer wall surfaces are populated by particles arriving from the inner and outer regions respectively, an imbalance occurs with more particles accumulating on the outer wall surfaces than the inner surfaces. The result is a nonzero ∆F/N which increases as Rl is increased or d is decreased. For very long Rl=50 in Fig. 3(d), additional features appear in ρv at the wall ends due to the periodic boundary conditions.
These results show that Casimir attractive forces can arise in active matter systems and that they can be controlled by modifying the run length of the swimming particles. The Casimir effect that we observe is not a critical Casimir effect since we are not at a critical point. Instead, the fluctuations or straight runs of the run-and-tumble particles are cut off by the walls, producing a relative drop in the number of particles on the inner walls.
Fig4.png
Figure 4: (a) ∆F/N vs d for a system with Rl=40 and various l. Negative forces appear when l approaches Ly. (b) ∆F/N vs l for fixed d=15 and Rl=40. (c) ∆F/N vs d for a system with l=59 at various Rl. (d) ρv(r) for a system with l=58, d=15, and Rl=40, in the negative force regime.

B.  Reverse Casimir effect

We next consider the effect of increasing l, the wall length. Since we have periodic boundary conditions, when l approaches Ly the system is better described as having apertures of width w=Lyl=60−l, as shown in Fig. 4(d) for a system with w=2. In Fig. 4(a) we plot ∆F/N versus d for increasing l, showing that for l ≥ 53, ∆F/N becomes negative over a range of d, while at larger d it becomes positive again. The magnitude of the negative force depends on l, and is largest for the smallest aperture sizes, as shown in Fig. 4(b). Here ∆F/N is nearly zero at small l and increases with increasing l up to l ≈ 40 due to the shadowing effect. For l > 40, ∆F/N decreases rapidly and becomes negative due to a particle trapping effect that occurs as the system enters the aperture limit. The particles moving along an outside wall return to the bulk upon reaching the wall end when l is small, but for large l they instead are trapped by the interior region, raising the interior density and producing a net repulsive force. This trapping effect becomes more prominent with decreasing w. In Fig. 4(a) regarding the curves with l ≥ 53 where negative forces can occur, for d < 30 the particles in the interior region spend most of their time running along the interior walls, producing a repulsive force. For d > 30 the particles have a chance to turn away from the inside wall during their run and avoid striking it, so the force becomes attractive again.
In Fig. 4(c) we plot ∆F/N versus d for systems with different Rl. For small Rl=10, ∆F/N is always positive, but as Rl increases, the force drops below zero over a region that increases with increasing Rl, and returns to positive values at higher values of d. The resulting shape of ∆F/N(d) resembles an atomic interaction force curve, with a stable characteristic distance determined by the point at which ∆F/N crosses zero with positive slope, suggesting that freely moving walls could be stabilized at this spacing. These results show that by varying the geometry it is possible to achieve detailed control over the magnitude of the fluctuation-induced forces which could be useful for self-assembly.
Fig5.png
Figure 5: Outward force −∆F/N vs d/Rl for particles confined between infinite walls at fixed particle density. (a) 1D systems with varied Rl (in legend) and d. The curves collapse to the functional form ∆F(d/Rl) = 1/(1+d/Rl). (b) 2D systems with various d and Rl. The curves collapse to a similar functional form, with altered coefficients.

C.  Infinite Walls

We have also considered the case of two infinite walls that confine the particles to an interior chamber, producing only an outward force. Such a system could be created by trapping bacteria or active particles in a confined geometry. We consider both one-dimensional (1D or linear) and 2D samples with different wall spacings d and constant particle density. For Brownian particles, ∆F/N is independent of d. In Fig. 5(a) we plot the outward force −∆F/N versus d/Rl for a 1D system. When l is infinite, the only length scales are d and Rl. The fraction of time a particle spends traveling along the walls is unchanged if both these lengths are scaled by the same factor, so the force must be a function of d/Rl. In Fig. 5, the curves for various Rl collapse to the form F(x)=1/(1+x) with xd/Rl, which can be derived from elementary considerations.
For simplicity, let xd/Rl be an integer, and measure length in units of Rl and time in units of τ. The particle moves in 1D along y with walls located at y=0 and y=x. During each run time, the particle moves a distance ±1 along y. To calculate the relative amount of time a particle spends at a wall, we note that a particle emerges into the bulk by tumbling away from a wall; taking this to be the y=0 wall, the particle takes one step to y=1. We then need to calculate the expected time tb required for the particle to either return to y=0 or reach y=x. This is the "gambler's ruin" problem, with the well-known solution [37] of (x−1). Counting the first step into the bulk, the expected number of steps in the bulk between wall motions is just x, so tb=x. To determine tw, the amount of time the particle spends at the wall after reaching it, we note that each tumble provides a 50:50 chance of escaping the wall, giving a probability p=[1/2] that tw=0, p=[1/4] that tw=1, and so on. The expected value of tw is easily found to be tw=1. Therefore, the ratio of tw to the total time is tw/(tw+tb)=1/(1+x). This equals the time-averaged force on the walls since in the 1D scenario, the force on a wall is 1 when a particle is present at the wall and 0 otherwise. For noninteger x, it is not too difficult to derive a more general formula [(2−[(x)/(ceil(x))])/(1+ceil(x))] where ceil(x) is the smallest integer ≥ x.
In Fig. 5(b) we plot −∆F/N for a 2D system. Here we could not derive a fully analytic solution; however, we obtain a fit of F(x) = a/(1+bx) with a=2/π and b ≈ 1.25. The value for a is exact and arises because in 2D all angles of incidence are possible rather than just normal incidence. Since the angle distribution is uniform, the force on the wall decreases by a factor of [1/([(π)/2]−0)]∫0π/2 cosθ dθ = 2/π. We cannot calculate b exactly, but it can be understood heuristically: compared to the 1D derivation, the particle spends more time in the bulk because only the x component of its motion moves it towards a wall. This increases the number of steps required to reach a wall, giving b > 1.
When steric particle-particle interactions are included, the forces on the walls will be reduced; however, in the dilute limit these results should hold and should be readily testable in experiments. It would also be interesting to study forces between two spheres or a sphere and a wall. In this case it would be important to understand how particles move along curved surfaces [36].

IV.  SUMMARY

We show that in a Casimir geometry of two parallel walls placed in a bath of run-and-tumble active particles, a robust attractive force arises between the walls due to a combination of particles moving along the walls and a geometric shadowing effect that depletes the particle density between the walls. The force increases when the depletion becomes more pronounced for increasing particle run length or reduced wall spacing. Our results are robust against the inclusion of steric particle-particle interactions in the dilute limit. For other geometries such as two walls containing small apertures, a particle trapping effect can produce repulsive forces between the walls. Our results show that active matter systems can exhibit a rich variety of fluctuation-induced forces between objects, which may be useful for applications such as self assembly or particle transport.
Note added- After submission we became aware of Ref. [38] on forces between two plates produced by active colloids with steric particle-particle interactions, where in the dilute limit a long-range attractive force arises of the same form that we observe. This indicates that our run-and-tumble results should also apply to active colloid systems. In the dense limit a repulsive force is caused by crystallization of the particles between the plates [38]. Additionally, Ref. [39] considers a Casimir type effect for swimming particles.
We thank D. Dalvit, L. Lopatina, and S. Redner for useful discussions. This work was carried out under the auspices of the NNSA of the U.S. DoE at LANL under Contract No. DE-AC52-06NA25396.

References

[1]
H.B.G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 (1948).
[2]
S.K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997).
[3]
F. Intravaia, S. Koev, I.W. Jung, A.A. Talin, P.S. Davids, R.S. Decca, V.A. Aksyuk, D.A.R. Dalvit, and D. Lopez, Nature Commun. 4, 3515 (2013).
[4]
J.N. Munday, F. Capasso, and V.A. Parsegian, Nature (London) 457, 170 (2009).
[5]
M.E. Fisher and P.G. de Gennes, C. R. Acad. Sci. Ser. B 287, 207 (1978).
[6]
C. Hertlein, L. Helden, A. Gambassi, S. Dietrich, and C. Bechinger, Nature (London) 451, 172 (2008).
[7]
A. Hanke, F. Schlesener, E. Eisenriegler, and S. Dietrich, Phys. Rev. Lett. 81, 1885 (1998); A. Gambassi, J. Phys.: Conf. Ser. 161, 012037 (2009).
[8]
D. Bonn, J. Otwinowski, S. Sacanna, H. Guo, G. Wegdam, and P. Schall, Phys. Rev. Lett. 103, 156101 (2009); A. Gambassi and S. Dietrich, Phys. Rev. Lett. 105, 059601 (2010).
[9]
S.J. Veen, O. Antoniuk, B. Weber, M.A.C. Potenza, S. Mazzoni, P. Schall, and G.H. Wegdam, Phys. Rev. Lett. 109, 248302 (2012); V.D. Nguyen, S. Faber, Z. Hu, G.H. Wegdam, and P. Schall, Nature Commun. 4, 1584 (2013).
[10]
S. Buzzaccaro, J. Colombo, A. Parola, and R. Piazza, Phys. Rev. Lett. 105, 198301 (2010).
[11]
C. Cattuto, R. Brito, U.M.B. Marconi, F. Nori, and R. Soto, Phys. Rev. Lett. 96, 178001 (2006); Y.Y. Villanueva, D.V. Denisov, S. de Man, and R.J. Wijngaarden, Phys. Rev. E 82, 041303 (2010).
[12]
N. Gnan, E. Zaccarelli, and F. Sciortino, Nature Commun. 5, 4267 (2014).
[13]
B.B. Machta, S.L. Veatch, and J.P. Sethna, Phys. Rev. Lett. 109, 138101 (2012).
[14]
S. Ramaswamy, Annu. Rev. Condens. Matter Phys. 1, 323 (2010); 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).
[15]
H.C. Berg, Random Walks in Biology (Princeton University Press, Princeton, 1983).
[16]
P. Galajda, J. Keymer, P. Chaikin, and R. Austin, J. Bacteriol. 189, 8704 (2007).
[17]
R. Dreyfus, J. Baudry, M.L. Roper, M. Fermigier, H.A. Stone, and J. Bibette, Nature (London) 437, 862 (2005).
[18]
W.F. Paxton, K.C. Kistler, C.C. Olmeda, A. Sen, S.K. St. Angelo, Y.Y. Cao, T.E. Mallouk, P.E. Lammert, and V.H. Crespi, J. Am. Chem. Soc. 126, 13424 (2004); J.R. Howse, R.A.L. Jones, A.J. Ryan, T. Gough, R. Vafabakhsh, and R. Golestanian, Phys. Rev. Lett. 99, 048102 (2007); R. Golestanian, Phys. Rev. Lett. 102, 188305 (2009).
[19]
G. Volpe, I. Buttinoni, D. Vogt, H.-J. Kummerer, and C. Bechinger, Soft Matter 7, 8810 (2011).
[20]
J. Palacci, S. Sacanna, A.P. Steinberg, D.J. Pine, and P.M. Chaikin, Science 339, 936 (2013).
[21]
I. Buttinoni, J. Bialke, F. Kümmel, H. Löwen, C. Bechinger, and T. Speck, Phys. Rev. Lett. 110, 238301 (2013).
[22]
J. Tailleur and M.E. Cates, Phys. Rev. Lett. 100, 218103 (2008); A.G. Thompson, J. Tailleur, M.E. Cates, and R.A. Blythe, J. Stat. Mech.: Theory Exp. 2011, P02029 (2011).
[23]
J. Bialke, T. Speck, and H. Löwen, Phys. Rev. Lett. 108, 168301 (2012).
[24]
Y. Fily and M.C. Marchetti, Phys. Rev. Lett. 108, 235702 (2012); G.S. Redner, M.F. Hagan, and A. Baskaran, Phys. Rev. Lett. 110, 055701 (2013); C. Reichhardt and C.J. Olson Reichhardt, Phys. Rev. E 90, 012701 (2014).
[25]
M.E. Cates and J. Tailleur, EPL 101, 20010 (2013).
[26]
M.B. Wan, C.J. Olson Reichhardt, Z. Nussinov, and C. Reichhardt, Phys. Rev. Lett. 101, 018102 (2008); C. Reichhardt and C.J. Olson Reichhardt, Phys. Rev. E 88, 062310 (2013).
[27]
J. Tailleur and M. E. Cates, Europhys. Lett. 86, 60002 (2009).
[28]
L. Angelani, A. Costanzo, and R. Di Leonardo, EPL 96, 68002 (2011); V. Kantsler, J. Dunkel, M. Polin, and R.E. Goldstein, Proc. Natl. Acad. Sci. (USA) 110, 1187 (2013); I. Berdakin, Y. Jeyaram, V.V. Moshchalkov, L. Venken, S. Dierckx, S.J. Vanderleyden, A.V. Silhanek, C.A. Condat, and V.I. Marconi, Phys. Rev. E 87, 052702 (2013).
[29]
M.E. Cates, Rep. Prog. Phys. 75, 042601 (2012).
[30]
L. Angelani, R. Di Leonardo, and G. Ruocco, Phys. Rev. Lett. 102, 048104 (2009); R. Di Leonardo, L. Angelani, D. Dell'Arciprete, G. Ruocco, V. Iebba, S. Schippa, M.P. Conte, F. Mecarini, F. De Angelis, and E. Di Fabrizio, Proc. Natl. Acad. Sci. U.S.A. 107, 9541 (2010).
[31]
A. Sokolov, M.M. Apodaca, B.A. Grzybowski, and I.S. Aranson, Proc. Natl. Acad. Sci. U.S.A. 107, 969 (2010).
[32]
N. Koumakis, A. Lepore, C. Maggi, and R. Di Leonardo, Nature Commun. 4, 2588 (2013).
[33]
L. Angelani and R. Di Leonardo, New J. Phys. 12, 113017 (2010).
[34]
L. Angelani, C. Maggi, M.L. Bernardini, A. Rizzo, and R. Di Leonardo, Phys. Rev. Lett. 107, 138302 (2011).
[35]
J.A. Drocco, C.J. Olson Reichhardt, and C. Reichhardt, Phys. Rev. E 85, 056102 (2012).
[36]
Y. Fily, A. Baskaran, and M.F. Hagan, arXiv:1402.5583.
[37]
S. Redner, A Guide to First-Passage Processes (Cambridge University Press, Cambridge, 2001).
[38]
R. Ni, M.A.C. Stuart, and P.G. Bolhuis, arXiv:1403.1533.
[39]
C. Parra-Rojas and R. Soto, Phys. Rev. E, in press.



File translated from TEX by TTHgold, version 4.00.
Back to Home