## Coarse-Grained Molecular Dynamics Simulations of Depletion-Induced Interactions for Soft Matter Systems

Tyler N. Shendruk, 1, ∗ Martin Bertrand, 2 James L. Harden, 2 Gary W. Slater, 2 and Hendrick W. de Haan 3

1 The Rudolf Peierls Centre for Theoretical Physics,

Department of Physics, Theoretical Physics, University of Oxford,

1 Keble Road, Oxford, OX1 3NP, United Kingdom

2 University of Ottawa, Department of Physics, 150 Louis-Pasteur, Ottawa, ON, K1N 6N5, Canada

3 University of Ontario Institute of Technology, Faculty of Science, 2000 Simcoe St. (Dated: August 5, 2014)

North, Oshawa, ON, Canada

Given the ubiquity of depletion effects in biological and other soft matter systems, it is desirable to have coarse-grained Molecular Dynamics simulation approaches appropriate for the study of complex systems. This paper examines the use of two common truncated Lennard-Jones (WCA) potentials to describe a pair of colloidal particles in a thermal bath of depletants. The shiftedWCA model is the steeper of the two repulsive potentials considered, while the combinatorial-WCA model is the softer. It is found that the depletion-induced well depth for the combinatorial-WCA model is significantly deeper than the shifted-WCA model because the resulting overlap of the colloids yields extra accessible volume for depletants. For both shifted- and combinatorial-WCA simulations, the second virial coefficients and pair potentials between colloids are demonstrated to be well approximated by the Morphometric Thermodynamics (MT) model. This agreement suggests that the presence of depletants can be accurately modelled in MD simulations by implicitly including them through simple, analytical MT forms for depletion-induced interactions. Although both WCA potentials are found to be effective generic coarse-grained simulation approaches for studying depletion effects in complicated soft matter systems, combinatorial-WCA is the more efficient approach as depletion effects are enhanced at lower depletant densities. The findings indicate that for soft matter systems that are better modelled by potentials with some compressibility, predictions from hard-sphere systems could greatly underestimate the magnitude of depletion effects at a given depletant density.

## 1. INTRODUCTION

form lock-and-key complexes [22-24] and to control particle motion and modify curvature in vesicles [25, 26].

Depletion forces of entropic origin are universal in soft matter systems, where they dictate the formation of a variety of soft solid phases [1-3]. In particular, for biological systems, these forces have been shown to play a prominent role in a myriad of phenomena [4, 5]. For example, we know depletion induced attractions promote fiber bundling [6], aggregation of red blood cells [7], and DNA collapse [8]. Reaction kinetics can also be influenced such as in DNA loop formation [5], protein folding [9, 10], and modified ligand-protein binding within nuclei [11]. Furthermore, these short-ranged depletion forces have been used to control biological systems, such as form active bundles [12]) and improved polymerase chain reactions [13].

Additionally, depletion forces occur in countless nonbiological soft matter contexts. They are often used to engineer novel microscopic systems by inducing ordered structures [14], as is done to align tobacco mosaic viruses [15-17]. Depletion forces can be significant enough to cause separation in colloid-polymer mixtures [3] and determine the phase diagrams of hardsphere mixtures [18, 19]. The presence of depletants can also alter the conformations of polymer chains [20, 21],

∗ Electronic address: tyler.shendruk@physics.ox.ac.uk

Yet, despite their experimental significance, simulations of complex situations that include depletioninduced interactions due to a thermal bath of smaller depletant particles are relatively rare. Computational research primarily focuses on the nature of the depletion forces themselves. A rich literature on binary hardsphere fluids has benefited from Density Functional Theories (DFT) (see Refs. [19, 27]). The majority of the work using DFT has been focused on understanding the statistical state of hard-sphere mixtures [18, 28-36]. Likewise, Monte Carlo (MC) simulations have played an essential role in verifying various depletion theories [37-40]. The computational efficiency of MC simulations is continually improving [41, 42] and they are being applied to more exotic mixtures than binary hard-spheres [17, 43, 44] and material properties [45].

Since such numerical calculations are computationally expensive, they have primarily been limited to simple situations. Furthermore, it is not always obvious how they can be used in conjunction with the generic particlebased computational tools that soft matter physicists have found to be useful for studying complex microscopic matter, such as Molecular Dynamics (MD) [46]. In order to study the physical phenomena that the depletants induce in complex soft matter systems, coarse-grained computational models that can seamlessly couple to those used by the soft matter community are needed. Previ-

ous work has used hard-core Yukawa potentials to study how additional interactions between depletants can impact depletion-induced effects [47], while molecular simulations have been applied to better understand gelation in depletant-colloid systems [39] with large size ratios [48] and the effects of using polymers as depletants [49-51].

Kim and Szleifer [52] studied the attraction between two approaching polymer chains in a bath of depletants via MD simulations using a truncated Lennard-Jones potential, which is also called the Week-Chandler Andersen (WCA) potential. The interaction was characterized by calculation of the potential of mean force (PMF), but as the system was quite extensive in terms of the number of particles, only a small number of volume fractions of depletants could be studied. In this work, we focus on an elementary system: two colloids in a bath of depletants (Fig. 1). This allows us to explore depletant interactions across a wide range of volume fractions for two common WCA implementations. In addition, the results are compared to theoretical predictions for the depletion potentials, including the Morphometric Thermodynamics model [36, 53, 54].

Our work demonstrates that the details of the excluded volume interactions used in the MD simulations greatly impacts depletion-induced forces. Indeed, an appropriate choice of parameters can significantly increase the depth of the attractive well, which in turn decreases the second virial coefficient and effective volume of the colloids, in a predictable manner. As a result, simulations using 'soft' yet steep potentials can reveal important entropic effects at lower volume fractions of depletants than hard sphere simulations and thus reduce computational costs for future studies into depletion effects in complex soft matter systems. Further, for soft matter systems in which the most appropriate coarse-grained model allows for some softness in the interaction between two constituents, the magnitude of depletion effects would be significantly underestimated by hard sphere models and better approximated by potentials that are not infinitely steep at the nominal edge of the particle.

## 1.1. Depletion-Induced Attraction

Because they can produce 'order from disorder' [55], at first glance entropic forces can feel counter-intuitive. In fact, because they are purely entropic in nature, depletion interactions are universal in biology and common in soft matter systems. Before discussing the simulation techniques or results found in this work, let us first review depletion-induced pair potentials.

As shown in Fig. 1, consider two large, hard, colloidal spheres of radius R c and volume V c suspended in a thermal bath of small, hard depletant spheres of radius R d and volume V d with thermal energy k B T = β -1 . The large colloids each have a thin sheath of radius R = R c + R d from which the centres of the depletants are excluded. When the two large spheres are brought

FIG. 1: Two relatively large spherical colloidal particles (of radius R c ) in a solution of depletants (radius R d &lt; R c ); each possesses an excluded volume of radius R = R c + R d . When the colloids are near (2 R c &lt; r &lt; 2 R ) an overlap volume V o exists causing the number of states available to the depletants to be greater than when the colloids are separate ( r ≥ 2 R ). This leads to an effective force that favours bringing the colloids together.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0309v1-figure-5ad0fd1e-8831-4eff-9213-0b49bc28e795.jpg)

together such that their centre-to-centre separation r is less than 2 R , their excluded volumes overlap and the free space available to the depletants increases. This results in a free energy change as a function of available volume for the depletants giving rise to an osmotic pressure pushing the colloids together.

## 1.2. Asakura-Oosawa Pair Potential

Asakura and Oosawa were the first to estimate this depletion-induced pair potential [56]. They discussed the attraction in terms of an entropic, depletant-induced (and therefore statistical in nature), short-ranged pair potential u ρ ( dep , r ). The Asakura-Oosawa (AO) model has also been referred to as the penetrating hard spheres (PHS) model [1] because the depletants are assumed to be hard with respect to the large colloids but an ideal gas amongst themselves, which is only appropriate for small volume fractions of depletants. In the thermodynamic, dilute limit of depletants ( N d → ∞ while φ d glyph[lessmuch] 1), one can approximate the osmotic pressure to be given by a van't Hoff relation Π ≈ -( φ /V d d ) β -1 . Explicitly, the

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0309v1-figure-0e940b2a-fc89-4ec2-a405-75f0a9012b39.jpg)

-

d

FIG. 2: Depletion-induced attractive pair potential between two spherical, hard colloids for a variety of colloid-depletant ratios. The inset shows the entropic pair potentials for various hard-sphere colloid and depletant sizes, while the main figure shifts these curves. The main abscissa is a shifted and scaled by the centre-to-centre separation emphasizing differnces in well depths. The hard-sphere Aakura-Oosawa (A0) model (Eq. (1) and Eq. (2)) is shown as dashed lines, while the Morphometric Thermodynamics (MT) predictions are shown as solid lines. Scaling the pair potentials removes φ d dependence from the AO model but not the MT model. The MT curves shown are for the example volume fraction of φ d = 0 2 .

AO entropic depletion pair potential is

$$\text{AO entropic depletion pair potential is} & & \text{ due to rep} \\ \beta u _ { \text{AO} } = \begin{cases} - \left ( \frac { \mathcal { R } } { R _ { d } } \right ) ^ { 3 } \left [ 1 - \frac { 3 } { 4 } \frac { r } { \mathcal { R } } + \frac { 1 } { 1 6 } \left ( \frac { r } { \mathcal { R } } \right ) ^ { 3 } \right ] \phi _ { d } & 2 R _ { c } \leq r \leq 2 \mathcal { R } \, & \text{Compu!} \\ r > 2 \mathcal { R } & \text{$\text{$\text{$\text{$\text{$\text{$by}$ and}$} } \\ ( 1 ) & \text{between $\cdot$} } \\ \text{Note that the AO model (and MT model, to be dis-} \\ \text{$therely s1$} \\ \text{$therely s2$} \end{cases}$$

Note that the AO model (and MT model, to be discussed below) is generally defined for hard spheres and must be supplemented with a steric-repulsion condition for r &lt; 2 R c . For hard spheres of radius R c ,

$$U _ { \text{HS} } = \begin{cases} \infty & r < 2 R _ { \text{c} } \\ 0 & r \geq 2 R _ { \text{c} }. \end{cases} \quad \begin{pmatrix} \P a \P a \\ ( 2 ) & \P ( \text{M} ) \\ \end{pmatrix}$$

Hence, the total interaction is u AO + U HS in this definition. As will be shown, we are free to relax this constraint when the particles are not strictly hard spheres.

The AO pair potential (Fig. 2; dashed lines) is characterized by four points: (i) the potential is strictly attractive; (ii) the range is set by the depletant size R d ; (iii) the minimum occurs when the large colloids are in contact

$$\beta u _ { A O } ^ { \min } = - \left [ 1 + \frac { 3 } { 2 } \frac { R _ { c } } { R _ { d } } \right ] \phi _ { d }, \quad \ ( 3 ) \quad \ u _ { \ } i s \\ \i _ { \dots \dots \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$$

which deepens as the size ratio R /R c d is increased (though is not zero in the limit R /R c d → 1 [57]); (iv)

via osmotic pressure, the volume fraction of depletants and temperature both scale the well depth.

The AO model is beautifully simple; however, it has two major limitations:

- · Depletants must be small R d glyph[lessmuch] R c .
- · The volume fraction must be small φ d glyph[lessmuch] 1.

The essential conclusion is that in order to consider the many interesting situations that do not meet these conditions, we must move beyond the AO limit.

## 1.3. Morphometric Pair Potential

To move beyond the φ d glyph[lessmuch] 1 limit, the free energy cost of inserting the large colloids must be greater than it is in the PHS assumption. Simply replacing the PHS osmotic pressure with a form that includes higher virial coefficients for the hard sphere gas of depletants allows the AO model to predict colloid pressures at higher φ d but is not sufficient to predict the form of the depletion pair potential [58, 59]. The missing ingredient is the realization that colloids that are close to one another mutually perturb the statistical arrangement of their surrounding shells of depletants and cause an inhomogeneous probability distribution of depletants. This results in the pair potential oscillating between attractive and repulsive as a function of separation [59]. The well effectively becomes deeper as it is now the difference between the maximum due to repulsion and the minimum due to attraction, but the minimum itself remains mostly unaffected.

Computationally intensive numerical techniques such as DFT can determine the local depletant particle density and so find the resulting entropic pair potential between colloids. Such techniques have proven extremely successful in exploring non-additive effects [30] but are not straight-forward analytically simple models. One conceptually simple model that can reproduce the first repulsive/anti-correlation component of the pair potential is called Morphometric Thermodynamics (MT) [36, 53, 54].

The MT model discusses free energy in terms of the geometric shapes involved and their thermodynamic conjugates. In the AO model, only the overlap volume V o and osmotic pressure Π associated with it are considered. Just as the osmotic pressure is a thermodynamic force that arises because of restrictions on the placement of depletants within the overlap volume, an entropic surface tension arises due to the cost of disallowing depletants to be at certain points on the shell of depletants surrounding each colloid due to the presence of the other colloid (Fig. 1). The entropic interaction pair potential u is then modelled as arising from changes to the accessible volume V o , to the surface area A restriction with resulting entropic surface tension γ and to the Gaussian curvatures of the overlap surface C 1 and C 2 with corre-

sponding entropic bending rigidities κ 1 and κ 2 :

$$u _ { \text{MT} } = \Pi V _ { o } + \gamma A + \kappa _ { 1 } C _ { 1 } + \kappa _ { 2 } C _ { 2 } \quad \ \ ( 4 ) \quad \ \widehat { \xi } ^ { \dagger }$$

where the geometric terms are functions of the colloidcolloid separation while the prefactors depend on the volume fraction of depletants (listed in Appendix B). As shown in Fig. 2 (solid lines), the MT framework is able to produce the attractive well and also the primary repulsive barrier, while the AO model only predicts the attractive well. Correlations and anti-correlations beyond the primary contributions cannot be reproduced within the MT model.

## 2. SIMULATION METHODS

Simplified analytical models such as the AO or MT potentials are appealing for straightforward systems but coarse-grained simulations offer a path to studying more complicated soft matter systems. In order to computationally study depletion effects via Molecular Dynamics simulations, a generic steep repulsive potential to produce hard steric effects between colloids and depletants is needed. To coarse-grain the system, we neglect all hydrodynamic interactions that may propagate through the fluid medium [60, 61] and any possible non-equilibrium effects [62]. We consider the solvent to be simply a dissipative thermal bath with Langevin Dynamics. Both the colloidal particles and the depletants are modelled as quasi-'hard spheres' in the framework of coarse-grained MD.

Because depletion forces arise due to excluded volumes between particles that are inaccessible to the smaller depletant particles, any repulsive potential can technically be used to represent the steric repulsion (though as we shall see, the details can have a significant impact on the magnitude of the effective interactions). In typical coarse-grained MD simulations, this is done using the Lennard-Jones (LJ) pairwise potential, which includes two terms, a steric repulsion that goes as r -12 and a dispersion-type interaction that goes as r -6 [46]. In order to neglect attractive dispersion-type attraction, the LJ potential is commonly truncated at its minimum and shifted to positive values, as will be explained in § 2.1 and § 2.2. This truncated LJ potential is referred to as the Weeks-Chandler-Andersen (WCA) potential. This is a simple and commonly employed potential; however, using a potential that is not infinitely steep alters the attractive entropic interactions and they must be carefully characterized. It will be seen in this work that an appropriate choice of parameters can significantly increase the depth of the depletion induced potential well, which in turn decreases the second virial coefficient and the effective volume of the colloids. We consider two versions of the WCA potential as coarse-grained generic simulation models for depletion effects. For both approaches, the characteristic energy for all interactions is glyph[epsilon1] and the

6 0

FIG. 3: Comparison of the shifted-WCA (blue) and combinatorial-WCA (red) potentials for R c = 2 5 . σ , ε = 1 k B T and ∆ ij = 4 σ for the shifted-WCA case. The combinatorialWCA extends further and diverges at a small separation and is thus the 'softer' of the two. However, both potentials are hard in that they rise as 1 /r 12 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0309v1-figure-c6cd97cd-ca65-4c61-bcd7-bd8f9ce880a3.jpg)

characteristic length scale is σ - which is equal to the diameter of the depletants.

## 2.1. Shifted WCA

The shifted-WCA model (sWCA) associates the same repulsive potential to each spherical particle but shifts the barrier depending on their radii. For a particle of nominal radius R i and one of R j , it is given by

$$\text{typical} \\ \text{ng the} \\ \text{icludes} \\ \text{and a} \\ \text{In or} - \\ \text{in, the} \\ \text{in and} \\ \mathfrak { I } \, \S \, 2. 1 \\ \mathfrak { I } \, \text{ to as} \quad. \quad. \quad.$$

where r ij is the centre-to-centre separation between particles i and j , ε is the depth of the potential well, σ is the nominal lengthscale of the potential, ∆ ij = R i + R j -σ is a shift of the potential and r cut ij = 2 1 / 6 σ +∆ ij is the cutoff distance. The repulsive potential thus rises in the same rapid manner, over the same distance of 2 2 / 6 σ , for all pairs regardless of the size of the particles. Here we have chosen σ and ε in Eq. (5) to explicitly be the characteristic values but one could conceivably chose some other distance and/or energy. As shown in Fig. 3, the shiftedWCA potential U sWCA rises from zero at r cut ij , crosses U sWCA = 1 ε at r ij = R i + R j , and diverges as r → ∆ . ij

## 2.2. Combinatorial WCA

On the other hand, the combinatorial-WCA model (cWCA) does not shift the WCA repulsion but rather sets an effective size between two MD particles. It is defined by where σ ij is no longer constant but instead is given by R i + R j and r cut ij = 2 1 / 6 σ . The resulting repulsion between particles U cWCA ( r ij ) is still hard - the potential still steeply rises as ∼ r -12 ij . However, the combinatorialWCA repulsion between two colloids is slightly softer than its shifted-WCA counterpart. As shown in Fig. 3, the repulsive potential extends further and does not rise as sharply. The same is true for colloid-depletant interactions. Since σ is equal to the size of the depletants, 2 R d , the depletant-depletant interactions are equivalent in the sWCA and cWCA potentials.

$$& \text{defined by} \\ & U = \begin{cases} 4 \varepsilon \left [ \left ( \frac { \sigma _ { i j } } { r _ { i j } } \right ) ^ { 1 2 } - \left ( \frac { \sigma _ { i j } } { r _ { i j } } \right ) ^ { 6 } \end{cases} \right ] + \varepsilon \quad r _ { i j } < r _ { i j } ^ { c u t } \quad & \text{$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$\quad$
& \text{where $\sigma_{ij}$ is no longer constant but instead is given by} \quad & \text{and} \quad r _ { i j } ^ { c u t } = 2 ^ { 1 / 6 } \sigma. \text{ The resulting repulsion be-} \quad & \text{and} \quad r _ { i j } ^ { c u t } = 2 ^ { 1 / 6 } \sigma. \text{ The resulting repulsion be-} \quad & \text{and} \quad r _ { i j } ^ { c u t } = 2 ^ { 1 / 6 } \sigma. \text{ The resulting repulsion be-} \quad & \text{and} \quad r _ { i j } ^ { c u t } \end{cases}$$

## 2.3. Effective Colloid Size

In using continuous potentials instead of perfectly hard spheres, the question arises as to what the effective size of the colloids are. Not only are the potentials not infinite at the nominal radius, but the colloid-colloid interaction extends beyond the 2 R c (Fig. 3). The second virial coefficient for a pairwise-interacting solution of colloids can be used to define an effective volume (and hence radius) for the colloids defined by either sWCA and cWCA:

Here U r ( ) is the potential that defines the colloids, which for hard spheres is given by Eq. (2). Hence for hardspheres, the second virial coefficient is 4 times the volume of the hard sphere B HS c 2 , = 4 ( 4 3 πR 3 C ) = 4 V HS .

$$r & \text{ the colloids defined by either swUCA and } c W U A \colon \\ B _ { c, 2 } = - 2 \pi \int _ { 0 } ^ { \infty } \left ( e ^ { - \beta U ( r ) } - 1 \right ) r ^ { 2 } d r. \quad ( 7 ) \quad T _ { c } \\ \text{ere} \, U \left ( r \right ) & \text{ is the potential that defines the colloids, which} \quad \text{coll} _ { c }$$

For the sWCA potential, we use Eq. (5) for colloids of nominal size 2 R c in Eq. (7) to solve for B c 2 , . Defining the effective volume of the sWCA colloids as V sWCA = B c 2 , / 4, we can calculate V sWCA /V HS , which depends on the nominal radius R c . For the case of R c = 2 5 . σ as studied in this work, we find that

$$\frac { V _ { s W C A } } { V _ { H S } } = 1. 0 1 0. \quad \ \ \ ( 8 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Hence, using the sWCA potential to define a colloid with a nominal diameter of 5 σ results in a colloid that has a 1% larger volume than the hard sphere equivalent. This corresponds to an effective radius of R sWCA = 2 508. .

Following a similar procedure for the cWCA potential but using Eq. (6), we obtain

$$\frac { V _ { s W C A } } { V _ { H S } } = 1. 0 5 3. \quad \ \ ( 9 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

The cWCA potential thus results in a colloid that has a 5.3% larger volume than the hard sphere result, corresponding to an effective radius of R cWCA = 2 543. . As opposed to the sWCA result, the length of the cWCA potential scales with the nominal size and thus this result is independent of R c .

The R sWCA and R cWCA values given above are used in theoretical calculations in the rest of this work. Although the corrections to the nominal radii appear to be small, they do have a significant effect on the results and yield much better agreement between the simulations and the theory for the pair potentials, well depths, and B c 2 , values. It is important to keep in mind that both of the WCA potentials are 'hard' in that they are both short-ranged, rapidly rising, repulsive potentials.

## 2.4. Net potential between two colloids

The net pair potential between two colloids is simply the sum of the excluded volume potentials U (Eq. (5) or (6)) and the entropic component u :

$$W = U \left ( r _ { i j }, R _ { c }, \sigma \right ) + u \left ( \phi _ { d }, r _ { i j }, R _ { c }, R _ { d } \right ). \quad ( 1 0 )$$

In this notation, W indicates the net potential, U denotes the repulsive, excluded volume potentials (WCA), and u specifies the entropic, depletion potentials. When the WCA potentials are subtracted off the net interaction potential then only the depletion-induced component remains.

## 3. SIMULATIONS OF DEPLETANT-INDUCED PAIR-POTENTIAL

To measure the entropic pair potential between two colloids, the mean force on each colloidal particle held at a fixed separation is measured. Integrating over the entire range of simulated separations produces the potential of mean force W sim between two colloidal particles.

## 3.1. Shifted-WCA

The net pair potential W sim sWCA ( r ), as measured by force integration, is shown in Fig. 4. The strong, shortranged, shifted-WCA repulsion U sWCA is clearly visible at small separations and the potential diverges at r = ∆ ij = 4 σ for these simulations. For clarity, The main figure shows low to moderately volume fractions ( φ d ≤ 0 25) and the inset shows the highest volume frac-. tions considered here ( φ d = 0 3 and 0 35). . . Subtracting U sWCA off leaves only the depletion-induced component u sWCA (Fig. 5). Persistently negative potentials for large centre-to-centre separations suggests that the force integrated results from simulations underpredict the net pair potentials at the highest volume fractions of depletants considered ( φ d &gt; . 0 25).

FIG. 4: Simulated net pair potentials W sim sWCA from integrated force measurements of shifted-WCA beads of size R c = 2 5 . σ in baths of various volume fractions of shifted-WCA depletants of size R d = 0 5 . σ . While the main figure shows low to moderate values, the inset shows the high volume fraction cases. Solid lines are the sum of the entropic pair potential predicted by the Morphometric Thermodynamics (MT) model and the shifted-WCA repulsion, u MT + U sWCA .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0309v1-figure-2cc6e7ab-4c9f-4943-b9de-31837062f736.jpg)

At low volume fractions, the simulated pair potentials are well approximated by both the AO and the MT models but as φ d glyph[greaterorsimilar] 0 1 the AO model begins to deviate from . the measured curves (Fig. 5; inset). On the other hand, the MT model continues to capture the behaviour of the primary attractive well and repulsive barrier (the region within approximately two depletants from the colloid's surface: 2 R c ≤ r glyph[lessorsimilar] 2 [ R c +2 R d ]). Of course, the MT model does not capture the secondary features in Fig. 5 because current functionals are not applicable beyond 2 R = 2( R c + R d ) [53].

Although, the AO model does not well represent the pair potential as a function of separation, Fig. 5 shows that it does estimate the contact energy well, i.e. u min sWCA ≈ u min AO . This can be seen quantitatively in Fig. 6. In fact, Fig. 6 shows that the linear AO minimum (Eq. (3)) agrees with the simulated min[ W sim sWCA ].

To better understand these predictions, we theoretically combine the MT model for the depletion potential and the sWCA excluded volume potential: u MT + U sWCA . This approach successfully replicates the net pair potential as compared to simulation for φ d glyph[lessorsimilar] 0 25 . (Fig. 4; main). This is because the depletion potential is well approximated by the MT model for all but the highest depletant volume fractions (Fig. 5). For the highest volume fractions considered (Fig. 4; inset), the MT model appears to overpredict the primary barrier height and underpredict the primary well depth for sWCA. The minima of the u MT + U sWCA data are shown in Fig. 6 as solid lines. While good agreement is found up to

FIG. 5: Simulated depletant-induced entropic pair potentials u sWCA from integrated force measurements of shifted-WCA beads of size R /R c d = 5. Solid lines show the induced pair potentials predicted by the Morphometric Thermodynamics (MT) model, while dashed lines in the inset show the AakuraOosawa (AO) model. These theoretical curves are zero for r &gt; 2 R = 6 σ and infinite for r &lt; 2 R c = 5 σ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0309v1-figure-261c4172-b828-43fa-bd27-4ca6eca1ffbf.jpg)

φ = 0 25, . min[ u MT + U sWCA ] actually begins to curve upwards thereafter. This trend is an artifact of the MT model where at high volume fractions, the first repulsive peak is exaggerated - which drags the curve upwards at r ≈ 2 R c and thus increases the minimum [54]. Due to this artifact, even the linear AO model does a better job than the MT model at predicting the minimum in this regime.

To correct this, the difference between the minimum of the contact attraction and the maximum of the anticorrelation barrier can be calculated. This calculation can be done for both the simulations (Fig. 6; red squares) and the u MT + U sWCA potential (this is not possible for the AO potential as no repulsive barrier is predicted). Excellent agreement is now found between the simulation results and the theoretical u MT + U sWCA potential across all studied volume fractions. This suggests that MT's overprediction of the primary barrier and underprediction of the primary well depth cancel.

The well depth competes with thermal energy in determining if colloids prefer to agglomorate in clusters or explore the entire available volume. Qualitatively speaking since the well depth W min sWCA in Fig. 6 is only a few k B T even as φ d → 0 4, . we expect that extremely high number densities of depletants would be required to overcome thermal jostling. This will be further quantified when we consider the second virial coefficients of colloids in thermal baths of depletants ( § 4.1).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0309v1-figure-64dc147c-c2d5-4606-ae1a-1e41cc66f646.jpg)

d

FIG. 6: Measures of the well depth of the pair potential as a function of depletant volume fraction for shifted-WCA particles of size R c = 2 5 . σ and depletants of size R d = 0 5 . σ . The minimum from simulations ( W sWCA ; black circles) is compared to the Aakura-Oosawa (A0) value ( u min AO from Eq. (3); dashed line) and the theoretical minima of the Morphometric Thermodynamics (MT) model plus the shifted-WCA potential (min[ u MT + U sWCA ]; solid black line). The effective well depths (primary minimum minus primary maximum) from simulations (∆ W sWCA ; red squares) and the theoretical net value (∆[ u MT + U sWCA ]; solid red line) are also shown.

## 3.2. Combinatorial-WCA

The depletant-induced pair potential is measured in the same manner as in § 3.1 for cWCA interactions; however, separations smaller than 2 R c are allowed by the cWCA repulsion and must now be considered as seen in Fig. 7. The competition between U cWCA and the entropic pair potential u cWCA results in a minimum in the net interaction W cWCA for all volume fractions. A comparison between Fig. 4 and Fig. 7 shows that the well is much deeper in the combinatorial-WCA model than for the shifted-WCA. This is because of the additional overlap that can occur. Since the U cWCA repulsion rises less dramatically and does not have an infinitely-hard repulsive core, the combinatorial-WCA wells are broader, extending to smaller separations. For high volume fractions, the combinatorial-WCA repulsion U cWCA does not dominate over the depletion-induced pair potential u cWCA until much smaller separations than in the shifted-WCA model. This can be seen in Fig. 8, which displays the location of the minimum of the net potentials W cWCA and W sWCA as measured in simulations as a function of volume fraction. Excellent agreement is again found between theory (MT) and simulation.

At low volume fraction, colloids are actually more likely to be found further apart for cWCA than sWCA since the cut-off for the former extends further. However,

FIG. 7: Simulated net pair potentials W cWCA from integrated force measurements of combinatorial-WCA beads of size R c = 2 5 . σ in baths of various volume fractions of combinatorial-WCA depletants of size R d = 0 5 . σ . The inset shows the high volume fraction cases and the main figure shows the low to moderate values. Solid lines are the sum of the MT model for the depletion potential and the combinatorial-WCA excluded volume potential, u MT + U cWCA .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0309v1-figure-e53e7698-e61c-43c8-9f77-44cb1ac7eb84.jpg)

due to the relative softness of u cWCA , the minimum in the potential decreases at a much faster rate and, at high volume fractions, colloids defined by cWCA are most likely found to be significantly closer together than for sWCA. It is interesting to note that the dependence of r min on φ d is approximately linear for both WCA potentials, with a cross-over around φ d ≈ 0 14 and . r min ≈ 5 1. .

Once again the WCA potential can be subtracted from the total potential W in order to produce the depletantinduced pair potential u cWCA (Fig. 9). The depletioninduced pair potential continues into separations r &lt; 2 R c . As discussed, the definitions of the AO and MT implicitly have infinitely hard cores at 2 R c . However, this is an imposed restriction and we can remove it without modifying the formalism.

As seen in Fig. 9 (inset), the AO potential crosses the simulated u potential quite near r ≈ 2 R c but when significant overlap occurs the AO potential u AO cannot reproduce the measured depletion-induced pair potential. On the other hand, allowing the MT pair potential u MT to extend to smaller separations does a remarkable job of reproducing the depletant-induced component of the pair potential except at the highest volume fractions of depletants considered (Fig. 7; inset) at which point the primary repulsive barrier is overpredicted. For this reason, summing u MT and U cWCA at low to moderate volume fractions once again reproduces the net interaction potentials.

Because the potentials are much softer, greater overlap

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0309v1-figure-b80841c9-4e0e-4a8a-96d6-de6b9cd67add.jpg)

d

FIG. 8: Position of the minimum in the net potential measured from simulations for both sWCA (black diamonds) and cWCA (red triangles). The prediction from the theoretical potentials u MT + U sWCA (solid black line) and u MT + U cWCA (solid red line) are also included.

between the colloids is allowed. This in turn yields a significant boost in the volume accessible to the depletants and thus a deep attractive well. The resulting minima (Fig. 10) are much deeper at high volume fractions than they were for the sWCA potentials (Fig. 6), though not at smaller volume fractions of φ d glyph[lessorsimilar] 0 2. . In fact, because the cWCA potential extends further (Fig. 3), the AO model is not even accurate at low volume fractions. At both high and low depletant volume fractions, the AO model for hard-spheres is not sufficient for these cWCA simulations. However, since summing the MT model and the cWCA potential reproduces the pair potentials so well for φ d glyph[lessorsimilar] 0 25 (Fig. 7), it also predicts the behaviour of . the minima and the well depths as functions of φ d quite well (Fig. 10). At the highest volume fractions considered ( φ d glyph[greaterorsimilar] 0 30), . the MT model once again begins to fail. While in the sWCA case, the primary repulsion was overpredicted but the primary well depth was underpredicted (Fig. 5 and Fig. 6), in the cWCA case both the repulsive barrier height and the well depth are overpredicted (Fig. 9 and Fig. 10).

## 4. VIRIAL COEFFICIENTS: IMPLICIT DEPLETANTS

Considering a reduced system of only two colloids in a bath of depletants has allowed us to characterize the depletant pair potential in detail. However, systems of interest generally consist of many large particles in a bath of depletants - a much more computationally expensive system to study. To circumvent this, the success of the MT potential in replicating the depletant pair potential

FIG. 9: Simulated depletant-induced entropic pair potentials u cWCA as a function of centre-to-centre separation for R /R c d = 5. Solid lines show the depletion-induced pair potentials predicted by the Morphometric Thermodynamics (MT) model, while dashed lines in the inset show the AakuraOosawa (AO) model. The AO model fails to accurately predict the entropic attraction for separations r &lt; 2 R c but the MT model remains highly accurate.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0309v1-figure-88f6586b-563f-4ee9-8837-37d58cb46e94.jpg)

indicates that it could be used to include the effects of the depletants implicitly . In this approach, explicit depletants would not be included in the system but their effects would be included by imposing the MT potential between all pairs of particles. We use the second virial coefficient as a quantitative measure to predict the efficacy of such an approach.

## 4.1. Second Virial Coefficient

Recall that the second virial coefficient was used to calculate the effective size of the colloids according to the potentials that define them. The effective size as modified by the presence of the depletants can be calculated by a similar method where the potential u r ( ) in Eq. (7) can now be: i) the simulated pair potentials W sim sWCA and W sim cWCA , ii) the sum of u MT and U (either U sWCA or U cWCA ), iii) the net AO potential u AO + U HS .

Although virial coefficient expansions of the equation of state breakdown in non-thermalized systems, predictions of B c 2 , &lt; 0 suggests that attraction dominates over thermal jostling and so demixing of colloidal and depletant particles will occur [54, 63]. In other words, the colloids are expected to agglomorate for negative enough B c 2 , values. The depletion-induced interactions demand a correction such that the virial coefficient is reduced from the hard-sphere value and so we shall discuss the reduced second virial coefficient B ∗ c 2 , = B c 2 , /V c .

FIG. 10: The same as Fig. 6 but for the combinatorial-WCA simulation model. Unlike the shifted-WCA case (Fig. 6), the minimum of the Aakura-Oosawa (AO; dashed line) model does not well-approximate the value from simulations at any volume fraction. The theoretical minimum from the MT model plus the combinatorial-WCA potential (solid black line) and the effective well depth are both found to be in good agreement with the corresponding simulation result.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0309v1-figure-423a65dd-c987-4e02-a1cb-8c53ddfa1291.jpg)

## 4.1.1. Simulations

Consider first the harder, shifted-WCA simulations (Fig. 11; black diamonds). The simulated shifted-WCA second virial coefficient is linear with volume fraction of depletants until approximately φ d ≈ 0 2 after . which point the curvature increases and B ∗ c 2 , begins to drop more rapidly. The second virial coefficient becomes negative somewhere in the range of φ d between 0 25 and 0 3. . .

The results for the combinatorial-WCA simulations (Fig. 11; red triangles) are similar to the sWCA results at low volume fractions. However, the combinatorialWCA B ∗ c 2 , falls quite rapidly. It becomes negative around φ d = 0 2 and drops to dramatically more negative values . at high volume fraction.

## 4.1.2. AO Approximation: φ d glyph[lessmuch] 1

Assuming that the colloids are hard and φ d glyph[lessmuch] 1 simplifies matters immensely. This is the AO model and, as shown in Appendix A, the depletant-induced attractions appear as a correction on the hard-sphere result, specifically taking the form

$$\frac { B _ { c, 2 } ^ { * } } { 4 } \approx 1 - g \phi _ { d }, \text{ \quad \ \ } ( 1 1 ) \text{ \ \ } \text{$\stackrel{----} { \text{bo} } \\ \text{At} \\ \text{um} \text{ }$$

where g ( R /R d c ) is a third order polynomial given by equation Eq. (A2). The AO approximation (Eq. (11))

FIG. 11: Reduced second virial coefficient B ∗ c 2 , = B c 2 , /V c for colloidal particles with R /R c d = 5 as a function of φ d . Results for the linear Aakura-Oosawa (AO) model is shown as a dashed line. The MT based predictions are also included: u MT + U sWCA (solid black line) and u MT + U cWCA (solid red line)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0309v1-figure-c823c003-4e91-41ea-b4c7-ae94b9f35e1d.jpg)

shows that, in this limit, the reduced second virial coefficient decreases as the volume fraction of depletants is increased, reflecting the deepening of the implicit depletion attraction. As is expected for the AO model, the correction from the hard-sphere value is linear in φ d . The linear AO prediction compares well to both the shiftedand the combinatorial-WCA simulations for R /R c d = 5 when the volume fraction of depletants is low (Fig. 11). The AO approximation remains relatively accurate for larger values of φ d for sWCA. For φ d glyph[lessorsimilar] 0 15, . the approximation predicts the second virial coefficient of the shifted-WCA simulations quite well. However, as it must, the AO approximation for B ∗ c 2 , fails at high volume fractions and does not cross to zero until φ d = 0 525 and . 0 527 for sWCA and cWCA, respectively. .

## 4.1.3. MT plus WCA Approximation

While discussing the simulation models, it was shown that summing the MT depletant potential and the WCA excluded volume interactions successfully predicts the total pair potentials. These theoretical pair potentials can be integrated over to predict the second virial coefficient as a function of volume fraction of depletants.

The result for u MT + U sWCA is shown as a solid black line in Fig. 11. The agreement between the AO curve and both simulation results is reasonable until φ d ≈ 0 20. . At this point, B ∗ c 2 , for u MT + U sWCA begins to curve upwards. This again is an artifact of the exaggeration of the repulsive barrier in the MT model at high volume fractions. On the other hand, the B ∗ c 2 , for the cWCA

simulations behaves physically, beginning to drop more rapidly at this point.

At low volume fractions, the result for u MT + U cWCA (solid red line) agrees with the data and in fact all other curves. However, as φ d increases, the B ∗ c 2 , for u MT + U cWCA begins to drop quite dramatically. The results are in remarkable agreement with the values derived from the simulations. Note that since the cWCA potential allows for the position of the minimum to be significantly less than 2 R c and for correspondingly deeper wells, the exaggeration of the repulsive barrier at higher volume fractions has a much smaller impact on the calculation of B ∗ c 2 , .

## 4.1.4. Escape Time

The second virial coefficient represents two-body interactions and is often conveniently interpreted as an effective size. To give an alternative physical interpretation that is appropriate when the second virial coefficient is a large, negative number, a set of simulations was performed in which two colloids are initially in contact. Setting R /R c d = 5, the volume fraction of depletants was varied between 0 and 0.35. Since the depletion force increases with the volume fraction, the particles are held together more tightly with increasing φ d . To quantify this, the average time it takes for the particles to separate was measured as an escape time, τ esc . The colloids are considered to have escaped the depletion potential well when their centre-to-centre distance is larger than 2( R c + R d ).

As shown in Fig. 12, the τ esc curve mirrors the B ∗ c 2 , results. When the volume fraction of depletants is large the primary attractive well is deep and the primary repulsive barrier is high (Fig. 10) such that colloids that are in contact must diffusively overcome the effective barrier of height ∆ W in order to escape contact. In this way, reaching r = 2 R is a Kramers' process and τ esc ∝ e β ∆ W [64, 65]. As the volume fraction of depletants is increased, the effective barrier height becomes greater (Fig. 10), the escape time rises rapidly as seen in Fig. 12 and the colloids effectively become stuck together. Concurrently, the second virial coefficient B ∗ c 2 , is plummeting to large negative values reflecting the strong attractive interactions as the volume fraction increases. The pair potentials (Fig. 7) demonstrate that the primary attractive well and repulsive barrier dominate over the secondary correlations and are indeed the significant contribution to the second virial coefficient. The inset to Fig. 12 shows that B ∗ c 2 , ∼ e β ∆ W - as for τ esc - when the effective barrier is greater than the thermal energy β ∆ W &gt; 1. The relationship τ esc ∼ B ∗ c 2 , suggests a simple physical picture: when primary barrier is large the magnitude of B ∗ c 2 , reflects the relative amount of time that colloids are held together by depletant forces τ esc .

FIG. 12: Overlay of the escape time τ esc onto the reduced second virial coefficient B ∗ c 2 , = B c 2 , /V c (triangles denoted cWCA simulations and the solid red line presents the corresponding MT-based predictions) as a function of φ d . The right-hand vertical axis corresponds to the escape time in simulation units. The inset shows the reduced second virial coefficient as an exponential function of the effective well depth, ∆ W .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0309v1-figure-412b084c-2527-4e52-b472-c7675aad6746.jpg)

## 5. CONCLUDING REMARKS

In order to explore coarse-grained methodologies for modelling depletion interactions in soft matter systems, we charaterized the use of two types of truncated Lennard-Jones (WCA) potentials in Molecular Dynamics simulations of large colloids immersed in a bath of smaller depletants.

The shifted-WCA model is the 'harder' of the two models while the combinatorial-WCA potential both extends further and rise less steeply - and is thus 'softer' (while still being a hard potential). We find that this seemingly small difference between the excluded volume potentials has a dramatic impact on the depletion effects that arise in the system; the well depths are deeper, the minimum separations are closer, and the second virial coefficients are much lower for cWCA than sWCA. Hence, modelling the system with a cWCA approach yields significantly enhanced depletion effects for a fixed volume fraction of depletants.

Comparison was also made between two simple theories for the depletion potentials. The simplest approach, the Asakura-Oosawa model for penetrating hardsphere depletants, was found to give good results in the limit of low volume fractions. However, as φ d increased, this agreement soon broke down for both WCA potentials. The geometric-based Morphometric Thermodynamics model not only reproduces the minimum in the depletion potential, but also the first repulsive barrier. Adding u MT to either U sWCA or U cWCA yielded excellent agreement for the net pair potentials, the depletion potentials, and the effective well depths at low and mod-

erate volume fractions ( φ d glyph[lessorsimilar] 0 25) but begins to fail for . the highest volume fractions considered here. However, the second virial coefficient calculated from u MT + U cWCA closely matched the cWCA simulation results across all volume fractions, although the predictions for sWCA exhibited known artifacts at these size ratios. This result suggests that using the MT depletion potential along side the cWCA excluded volume interaction may be an effective model for implicitly incorporating depletion effects into future coarse-grained MD simulations.

The relative softness of cWCA significantly enhances depletion effects and thus presents a computationally efficient simulation approach since accessible volume fractions generate large primary attractive wells and repulsive barriers. Moreover, not only does the use of WCA potentials allow for MD-based exploration of depletant effects but also is likely to yield more physical results for a host of soft matter systems in which interactions are softer than the hard-sphere model. In many soft matter systems the effective excluded volume interaction has a softness associated with it; examples include polymer coated colloids, solutes that are compressible (e.g. micro or nanogels, dendrimers), charged objects in an electrolyte solution where the interaction is governed by overlapping Debye layers or easily deformed materials that can fit together more effectively (e.g., protein folding, chromosome collapse [66]). The results from hard

- [1] H. Lekkerkerker and R. Tuinier. Colloids and the Depletion Interaction . Springer, 2011.
- [2] J. Bibette, F. Leal Calderon, and P. Poulin. Emulsions: basic principles. Reports on Progress in Physics , 62(6):969, 1999.
- [3] R. Tuinier, J. Rieger, and C.G. de Kruif. Depletioninduced phase separation in colloidpolymer mixtures. Advances in Colloid and Interface Science , 103(1):1 - 31, 2003.
- [4] D. Marenduzzo, K. Finan, and P. R. Cook. The depletion attraction: an underappreciated force driving cellular organization. The Journal of Cell Biology , 175(5):681-686, 2006.
- [5] D. Marenduzzo, C. Micheletti, and P. R. Cook. Entropydriven genome organization. Biophysical Journal , 90(10):3712 - 3721, 2006.
- [6] M. Hosek and J. X. Tang. Polymer-induced bundling of f actin and the depletion force. Physical Review E , 69:051907, May 2004.
- [7] P. Steffen, C. Verdier, and C. Wagner. Quantification of depletion-induced adhesion of red blood cells. Physical Review Letters , 110:018102, Jan 2013.
- [8] J. Pelletier, K. Halvorsen, B.-Y. Ha, R. Paparcone, S. J. Sandler, C. L. Woldringh, W. P. Wong, and S. Jun. Physical manipulation of the escherichia coli chromosome reveals its soft nature. Proceedings of the National Academy of Sciences , 2012.
- [9] Y. Zhou and C. K. Hall. Solute excluded-volume effects on the stability of globular proteins: A statistical thermodynamic theory. Biopolymers , 38(2):273-284, 1996.

sphere calculations thus represent the lower bound for estimating depletion effects at a given volume fraction and may significantly underestimate their impact for such soft matter systems. While the correct modelling will be system specific, this possibility of underestimating depletant effects is important when considering predictions from hard sphere models. Our results have shown that even a small amount of softness leading to a marginal amount of overlap can result in a dramatic increase in the depletion potential, resulting in 'soft' colloids bound together by depletion-induced attraction that would be insufficient for 'harder' colloids.

## Acknowledgements

We gratefully acknowledge support through NSERC Discovery Grants to G.W.S. and J.L.H. and EMBO funding to T.N.S (ALTF181-2013). High performance computational resources were graciously provided by Prof. Laura Ramuno at the University of Ottawa and Sharcnet.

## References

- [10] Y. Snir and R. D. Kamien. Entropically driven helix formation. Science , 307(5712):1067, 2005.
- [11] K. Richter, M. Nessling, and P. Lichter. Macromolecular crowding and its potential impact on nuclear function. Biochimica et Biophysica Acta (BBA) - Molecular Cell Research , 1783(11):2100 - 2107, 2008.
- [12] T. Sanchez, D. T. N. Chen, S. J. DeCamp, M. Heymann, and Z. Dogic. Spontaneous motion in hierarchically assembled active matter. Nature , 491(7424):431-434, 2012.
- [13] R. R. Lareu, K. S. Harve, and M. Raghunath. Emulating a crowded intracellular environment in vitro dramatically improves rt-pcr performance. Biochemical and Biophysical Research Communications , 363(1):171 - 177, 2007.
- [14] A. G. Yodh, K. Lin, J. C. Crocker, A. D. Dinsmore, R. Verma, and P. D. Kaplan. Entropically driven selfassembly and interaction in suspension. 359(1782):921937, 2001.
- [15] M. Adams, Z. Dogic, S. L. Keller, and S. Fraden. Entropically driven microphase transitions in mixtures of colloidal rods and spheres. Nature , 393(6683):349-352, 1998.
- [16] S. C. Glotzer and M. J. Solomon. Anisotropy of building blocks and their assembly into complex structures. Nature materials , 6(8):557-562, 2007.
- [17] W. Li, T. Yang, and H.-r. Ma. Depletion potentials in colloidal mixtures of hard spheres and rods. The Journal of Chemical Physics , 128(4), 2008.
- [18] R. Roth, R. Evans, and S. Dietrich. Depletion potential in hard-sphere mixtures:theory and applications. Physical Review E , 62:5360-5377, Oct 2000.

- [19] R. Roth. Fundamental measure theory for hard-sphere mixtures: a review. Journal of Physics: Condensed Matter , 22(6):063102, 2010.
- [20] R. P. Sear. Coil-globule transition of a semiflexible polymer driven by the addition of spherical particles. Physical Review E , 58:724-728, Jul 1998.
- [21] B. Martnez-Haya and M. C. Gordillo. Effect of macromolecular crowding on the conformation of confined chain polymers. Macromolecular Theory and Simulations , 14(7):421-427, 2005.
- [22] S. Sacanna, W. T. M. Irvine, P. M. Chaikin, and D. J. Pine. Lock and key colloids. Nature , 464(7288):575-578, 2010.
- [23] T. Wang, R. Sha, R. Dreyfus, M. E. Leunissen, C. Maass, D. J. Pine, P. M. Chaikin, and N. C. Seeman. Selfreplication of information-bearing nanoscale patterns. Nature , 478(7368):225-228, 2011.
- [24] D. J. Ashton, Robert L. Jack, and N. B. Wilding. Selfassembly of colloidal polymers via depletion-mediated lock and key binding. Soft Matter , 9:9661-9666, 2013.
- [25] A. D. Dinsmore, A. G. Yodh, and D. J. Pine. Entropic control of particle motion using passive surface microstructures. Nature , 383(6597):239-242, 1996.
- [26] A. D. Dinsmore, D. T. Wong, P. Nelson, and A. G. Yodh. Hard spheres in vesicles: Curvature-induced forces and particle-induced curvature. Physical Review Letters , 80:409-412, Jan 1998.
- [27] M. Schmidt. Geometry-based density functional theory: an overview. Journal of Physics: Condensed Matter , 15(1):S101, 2003.
- [28] R. Roth, B. G¨ otzelmann, and S. Dietrich. Depletion forces near curved surfaces. Physical Review Letters , 83:448-451, Jul 1999.
- [29] M. Schmidt. Density functional for additive mixtures. Physical Review E , 62:3799-3802, 2000.
- [30] R. Roth and R. Evans. The depletion potential in non-additive hard-sphere mixtures. Europhysics Letters , 53(2):271, 2001.
- [31] M. Oettel. Depletion force between two large spheres suspended in a bath of small spheres: Onset of the derjaguin limit. Physical Review E , 69:041404, 2004.
- [32] S. A. Egorov. Effect of repulsive and attractive interactions on depletion forces in colloidal suspensions: A density functional theory treatment. Physical Review E , 70:031402, Sep 2004.
- [33] R. Roth and M. Kinoshita. Depletion potential between large spheres immersed in a multicomponent mixture of small spheres. The Journal of Chemical Physics , 125(8), 2006.
- [34] P.-M. K¨ onig, R. Roth, and S. Dietrich. Depletion forces between nonspherical objects. Physical Review E , 74:041404, 2006.
- [35] A. R. Herring and J. R. Henderson. Hard-sphere fluid adsorbed in an annular wedge: The depletion force of hardbody colloidal physics. Physical Review E , 75:011402, 2007.
- [36] M. Oettel, H. Hansen-Goos, P. Bryk, and R. Roth. Depletion interaction of two spheres full density functional theory vs. morphometric results. Europhysics Letters , 85(3):36003, 2009.
- [37] T. Biben, P. Bladon, and D. Frenkel. Depletion effects in binary hard-sphere fluids. Journal of Physics: Condensed Matter , 8(50):10799, 1996.
- [38] R. Dickman, P. Attard, and V. Simonian. Entropic forces
- in binary hard sphere mixtures: Theory and simulation. The Journal of Chemical Physics , 107(1), 1997.
- [39] L. Lue and L. V. Woodcock. Depletion effects and gelation in a binary hard-sphere fluid. Molecular Physics , 96(9):1435-1443, 1999.
- [40] A. Santos, S. B. Yuste, and M. L´pez de Haro. o Contact values of the particle-particle and wall-particle correlation functions in a hard-sphere polydisperse fluid. The Journal of Chemical Physics , 123(23), 2005.
- [41] J. Liu and E. Luijten. Rejection-free geometric cluster algorithm for complex fluids. Physical Review Letters , 92:035504, 2004.
- [42] D. J. Ashton, V. S´nchez-Gil, and N. B. Wilding. Monte a carlo methods for estimating depletion potentials in highly size-asymmetrical hard sphere mixtures. The Journal of Chemical Physics , 139(14), 2013.
- [43] G. Odriozola, F. Jim´ enez- ´ ngeles, A and M. LozadaCassou. Entropy driven key-lock assembly. The Journal of Chemical Physics , 129(11), 2008.
- [44] D. A. Triplett and K. A. Fichthorn. Entropic forces and directed alignment of hard squares in suspensions of rods and disks. The Journal of Chemical Physics , 133(14), 2010.
- [45] J. M. Rickman and H. D. Ou-Yang. Elastic response of binary hard-sphere fluids. Physical Review E , 84:012401, 2011.
- [46] G. W. Slater, C. Holm, M. V. Chubynsky, H. W. de Haan, A. Dub´ e, K. Grass, O. A. Hickey, C. Kingsburry, D. Sean, T. N. Shendruk, and L. Zhan. Modeling the separation of macromolecules: A review of current computer simulation methods. Electrophoresis , 30(5):792-818, 2009.
- [47] A. A. Louis, E. Allahyarov, H. L¨wen, and R. Roth. Efo fective forces in colloidal mixtures: From depletion attraction to accumulation repulsion. Physical Review E , 65:061407, 2002.
- [48] D. Henderson, A. Trokhymchuk, L. V. Woodcock, and K.-Y. Chan. Simulation and approximate formulae for the radial distribution functions of highly asymmetric hard sphere mixtures. Molecular Physics , 103(5):667674, 2005.
- [49] J. S. Smith, D. Bedrov, and G. D. Smith. A molecular dynamics simulation study of nanoparticle interactions in a model polymer-nanoparticle composite. Composites Science and Technology , 63(11):1599 - 1605, 2003. Modeling and Characterization of Nanostructured Materials.
- [50] S.-C. Kim, B.-S. Seong, and S.-H. Suh. Depletion interactions in two-dimensional colloid-polymer mixtures: Molecular dynamics simulations. Journal of Physics: Condensed Matter , 21(3):035101, 2009.
- [51] S. Kim. Molecular dynamics study of soft matter: Potential dependence of depletion forces. Journal of the Korean Physical Society , 61(7):1102-1107, 2012.
- [52] J. S. Kim and I. Szleifer. Depletion effect on polymers induced by small depleting spheres. The Journal of Physical Chemistry C , 114(48):20864-20869, 2010.
- [53] V. Bot ¸an, F. Pesth, T. Schilling, and M. Oettel. Hardsphere fluids in annular wedges: Density distributions and depletion potentials. Physical Review E , 79:061402, Jun 2009.
- [54] D. J. Ashton, N. B. Wilding, R. Roth, and R. Evans. Depletion potentials in highly size-asymmetric binary hardsphere mixtures: Comparison of simulation results with theory. Physical Review E , 84:061136, Dec 2011.

- [55] R. B. Phillips, J. Kondev, J. Theriot, N. Orme, and H. Garcia. Physical biology of the cell . Garland Science New York, 2009.
- [56] S. Asakura and F. Oosawa. On interaction between two bodies immersed in a solution of macromolecules. The Journal of Chemical Physics , 22(7):1255-1256, 1954.
- [57] A. Zaccone and E. M. Terentjev. Theory of molecular crowding in brownian hard-sphere liquids. Physical Review E , 85:061202, Jun 2012.
- [58] H. N. W. Lekkerkerker. Osmotic equilibrium treatment of the phase separation in colloidal dispersions containing non-adsorbing polymer molecules. Colloids and Surfaces , 51(0):419 - 426, 1990.
- [59] B. G¨ otzelmann, R. Evans, and S. Dietrich. Depletion forces in fluids. Phys. Rev. E , 57:6785-6800, Jun 1998.
- [60] J. Bafaluy, B. Senger, J.-C. Voegel, and P. Schaaf. Effect of hydrodynamic interactions on the distribution of adhering brownian particles. Physical Review Letters , 70:623-626, 1993.
- [61] J. T. Padding and A. A. Louis. Hydrodynamic interactions and brownian forces in colloidal suspensions: Coarse-graining over time and length scales. Physical Review E , 74:031402, Sep 2006.
- [62] J. Dzubiella, H. L¨wen, and C. N. Likos. Depletion forces o in nonequilibrium. Physical Review Letters , 91:248301, 2003.
- [63] M. G. Noro and D. Frenkel. Extended correspondingstates behavior for particles with variable range attractions. The Journal of Chemical Physics , 113(8), 2000.
- [64] H. A. Kramers. Brownian motion in a field of force and the diffusion model of chemical reactions. Physica , 7(4):284 - 304, 1940.
- [65] P. H¨ anggi, P. Talkner, and M. Borkovec. Reaction-rate theory: fifty years after Kramers. Reviews of Modern Physics , 62:251-341, Apr 1990.
- [66] T. N. Shendruk, M. Bertrand, H. W. de Haan, J. L. Harden, and G. W. Slater. Simulating the entropic collapse of coarse-grained chromosomes. http://arxiv.org/abs/1407.2850, 2014.
- [67] Y. Rosenfeld. Free-energy model for the inhomogeneous hard-sphere fluid mixture and density-functional theory of freezing. Physical Review Letters , 63:980-983, Aug 1989.
- [68] R. Roth, R. Evans, A. Lang, and G. Kahl. Fundamental measure theory for hard-sphere mixtures revisited: the white bear version. Journal of Physics: Condensed Matter , 14(46):12063, 2002.
- [69] Y.-X. Yu and J. Wu. Structures of hard-sphere fluids from a modified fundamental-measure theory. The Journal of Chemical Physics , 117(22):10156-10164, 2002.
- [70] H. Hansen-Goos and R. Roth. Density functional theory for hard-sphere mixtures: the white bear version mark ii. Journal of Physics: Condensed Matter , 18(37):8413, 2006.

## Appendix A: Second Virial Coefficient in the φ d glyph[lessmuch] 1 Limit:

## Depletant/Colloid van der Waals Gas

In the φ d glyph[lessmuch] 1 (or equivalently β u ∣ ∣ min AO ∣ ∣ glyph[lessmuch] 1) AO limit, the minimum value of the pair potential is well approximated by the AO model of hard-sphere colloids (Eq. (3))

and furthermore the well depth is less than the thermal energy ( β u ∣ ∣ min AO ∣ ∣ glyph[lessmuch] 1). This limit allows us to approximate the Mayer function f ( r ) = 1 -e -βu r ( ) ≈ βu r ( ), which can then be used to find an approximate, analytical form of the second virial coefficient:

$$\begin{array} { c } \text{can form on the second year} \text{elements} \\ \text{r} \\ B _ { c, 2 } = 2 \pi \left [ \frac { 2 } { \pi } \frac { 4 \pi R _ { c } ^ { 3 } } { 3 } - \int _ { 2 R _ { c } } ^ { 2 \mathcal { R } } \left ( e ^ { - \beta u } - 1 \right ) r ^ { 2 } d r + 0 \right ] \\ g \\ s, & \approx 2 \pi \left [ \frac { 2 } { \pi } \frac { 4 \pi R _ { c } ^ { 3 } } { 3 } + \int _ { 2 R _ { c } } ^ { 2 \mathcal { R } } \beta u _ { \mathrm A O } r ^ { 2 } d r \right ] \\ n & = 4 V _ { c } \left [ 1 - \phi _ { d } g \right ], & \text{(A1)} \\ \text{f-} \\ \text{f-} \end{array}$$

where we have defined

$$\underset { i c a l } { \text{in-} } \underset { g } { \text{s:} } \quad & g \left ( \frac { R _ { d } } { R _ { c } } \right ) \equiv \frac { 3 } { 2 } + \frac { 1 5 } { 8 } \left ( \frac { R _ { d } } { R _ { c } } \right ) + \frac { 3 } { 4 } \left ( \frac { R _ { d } } { R _ { c } } \right ) ^ { 2 } + \frac { 1 } { 8 } \left ( \frac { R _ { d } } { R _ { c } } \right ) ^ { 3 } \\ \underset { \text{c-es} } { \text{s:} }$$

for convenience. Notice that setting φ d = 0 automatically reproduces the hard-sphere value B c 2 , = 4 V c . The number density of depletants scales the correction, which otherwise is a function of R /R d c alone (encapsulated in g ). Less intuitively, the effective volume is predicted to reduce when the size of the depletants relative to the colloids is increased. This is a manifestation of the increased range of the attractive well but one must keep in mind the discussion in that the depletants must be significantly smaller than the colloidal particles. This condition remains.

The AO solution can be described as a van der Waals gas but, unlike a standard van der Waals gas in which only the first term has thermal energy dependence, in an AO solution the pressure scales with β -1 , reflecting the entropic basis of the attractive forces.

## Appendix B: Morphometric Pair Potential Details

For the MT potential defined in Eqn. 4, the geometric variables are simple functions of colloid-colloid separation r [36]:

$$V _ { 0 } & = \frac { 4 \pi \mathcal { R } ^ { 3 } } { 3 } \left [ 1 - \frac { 3 } { 2 } \left ( \frac { r } { 2 \mathcal { R } } \right ) + \frac { 1 } { 2 } \left ( \frac { r } { 2 \mathcal { R } } \right ) ^ { 3 } \right ], \quad \text{(B1a)} \\ A & = 4 \pi \mathcal { R } ^ { 2 } \left ( 1 - \frac { r } { \sim } \right ).$$

$$A & = 4 \pi \, \kappa \, \left ( 1 - \overline { 2 \mathcal { R } } \right ), & \text{(D1D)} \\, & C _ { 1 } & = \frac { A } { \mathcal { R } } + \pi \mathcal { R } \sqrt { 1 - \left ( \frac { r } { 2 \mathcal { R } } \right ) ^ { 2 } } \sin ^ { - 1 } \left ( \frac { r } { 2 \mathcal { R } } \right ), & \text{(B1c)} \\ | & C _ { 2 } & = 4 \pi. & \text{(B1d)}$$

$$\L ^ { \cdot \nu } & \quad 3 \quad \lfloor \, ^ { \cdot } \, \quad 2 \, \Lambda \, \Lambda \, \Lambda \, \quad 2 \, \Lambda \, \Lambda \, \quad \rfloor \,, \quad \langle \, \Lambda \, \Lambda \, \cdots \, \Lambda \, \rangle \\ A & = 4 \pi \mathcal { R } ^ { 2 } \left ( 1 - \frac { r } { 2 \mathcal { R } } \right ), \quad \ \ ( B 1 b ) \\ C. & = \frac { A } { 2 \, \perp \, \pi \mathcal { P }. } \sqrt { 1 \, - \, \left ( \frac { r } { 2 \, \perp \, \Lambda } \, \right ) ^ { 2 } } _ { \text{cin} ^ { - 1 } \left ( \frac { r } { 2 \, \perp \, \Lambda } \, \right ) } \quad \ \ ( \mathcal { R } 1 \in \Lambda )$$

$$C _ { 2 } & = 4 \pi.$$

However, the thermodynamic quantities are less straightforward. We utilize the Rosenfeld functionals [53, 67]:

$$\beta \Pi V _ { d } & = - \phi _ { d } \frac { 1 + \phi _ { d } + \phi _ { d } ^ { 2 } } { \left ( 1 - \phi _ { d } \right ) ^ { 3 } }, & \left ( B 2 a \right ) \left ( \begin{matrix} \beta \\ \phi _ { d }, \, w \\ \text{terms} \end{matrix} \right ) \\ \varepsilon \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{, } \text{ } \text{, } \text{ } \text{, } \text{, } \text{ } \text{ } \text{, } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text$$

$$\beta \kappa _ { 1 } V _ { \text{d} } = - \phi _ { \text{d} } R _ { \text{d} } ^ { 2 } & \frac { \lambda ^ { \prime \prime } _ { \text{d} } } { \left ( 1 - \phi _ { \text{d} } \right ) ^ { 3 } }, \\ & \quad \, \Big / \, - \varepsilon + 7 \theta \, - \, 1 1 \theta ^ { 2 } \quad \ln \, ( 1 - \theta \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot$$

$$\beta \gamma V _ { \text{d} } = \phi _ { \text{d} } R _ { \text{d} } \frac { \stackrel { \ell } { 3 \phi _ { \text{d} } } \stackrel { \ell \mu } { ( 1 + \phi _ { \text{d} } ) } } { 2 \left ( 1 - \phi _ { \text{d} } \right ) ^ { 3 } }, \quad \quad \quad ( \text{B2b} ) \\ \wedge \ \. \ \. \ \. \ \. \ \. \ \.$$

$$& \stackrel { r \cdots \iota } { r \cdots \iota } \quad \stackrel { r \cdots \iota } { \left ( 1 - \phi _ { d } \right ) ^ { 3 } } \,, \quad \stackrel { \vee \, \, \,$$

but other accurate forms are available in the binary hard sphere fluid literature such as the White Bear functionals [68-70]. The leading term in the osmotic pressure is φ d , while in both the surface tension and κ 2 , the leading terms are φ 2 d . In κ 1 the leading term is φ 3 d .
Tyler N. Shendruk Martin Bertrand James L. Harden Gary W. Slater Hendrick W. de Haan
Coarse-Grained Molecular Dynamics Simulations of Depletion-Induced Interactions for Soft Matter Systems
Given the ubiquity of depletion effects in biological and other soft matter
systems, it is desirable to have coarse-grained Molecular Dynamics simulation
approaches appropriate for the study of complex systems. This paper examines
the use of two common truncated Lennard-Jones (WCA) potentials to describe a
pair of colloidal particles in a thermal bath of depletants. The shifted-WCA
model is the steeper of the two repulsive potentials considered, while the
combinatorial-WCA model is the softer. It is found that the depletion-induced
well depth for the combinatorial-WCA model is significantly deeper than the
shifted-WCA model because the resulting overlap of the colloids yields extra
accessible volume for depletants. For both shifted- and combinatorial-WCA
simulations, the second virial coefficients and pair potentials between
colloids are demonstrated to be well approximated by the Morphometric
Thermodynamics (MT) model. This agreement suggests that the presence of
depletants can be accurately modelled in MD simulations by implicitly including
them through simple, analytical MT forms for depletion-induced interactions.
Although both WCA potentials are found to be effective generic coarse-grained
simulation approaches for studying depletion effects in complicated soft matter
systems, combinatorial-WCA is the more efficient approach as depletion effects
are enhanced at lower depletant densities. The findings indicate that for soft
matter systems that are better modelled by potentials with some
compressibility, predictions from hard-sphere systems could greatly
underestimate the magnitude of depletion effects at a given depletant density.