Domain size polydispersity effects on the structural and dynamical properties in lipid monolayers with phase coexistence. | |

Elena Rufeil Fiori and Adolfo J. Banchio | |

Lipid monolayers with phase coexistence are a frequently used model for lipid membranes. In these systems, domains of the liquid–condensed phase always present size polydispersity. However, very few theoretical works consider size distribution effects on the monolayer properties. Because of the difference in surface densities, domains have excess dipolar density with respect to the surrounding liquid expanded phase, originating a dipolar inter–domain interaction. This interaction depends on the domain area, and hence the presence of a domain size distribution is associated with interaction polydispersity. Inter–domain interactions are fundamental to understanding the structure and dynamics of the monolayer. For this reason, it is expected that polydispersity significantly alters monolayer properties. By means of Brownian dynamics simulations, we study the radial distribution function (RDF), the average mean square displacement and the average time–dependent self–diffusion coefficient, , of lipid monolayers with normal distributed size domains. For this purpose, we vary the relevant system parameters, polydispersity and interaction strength, within a range of experimental interest. We also analyze the consequence of using a monodisperse model for determining the interaction strength from an experimental RDF. It was found that polydispersity strongly affects the value of the interaction strength obtained, which is greatly underestimated if polydispersity is not considered. However, within a certain range of parameters, the RDF obtained from a polydisperse model can be well approximated by that of a monodisperse model, suitably fitting the interaction strength, even for polydispersities. For small interaction strengths or small polydispersities, the polydisperse systems obtained from fitting the experimental RDF have an average mean square displacement and in good agreement with that of the monodisperse system. | |

^{0}

^{0}footnotetext: Universidad Nacional de Córdoba, Facultad de Matemática, Astronomía, Física y Computación, Córdoba, Argentina.

^{0}

^{0}footnotetext: Consejo Nacional de Investigaciones Científicas y Técnicas, CONICET, IFEG, Córdoba, Argentina.

^{0}

^{0}footnotetext: e-mail: rufeil@famaf.unc.edu.ar

## 1 Introduction

Most biologically relevant monolayers
present phase coexistence characterized by domains formed by lipids in an ordered phase state, dispersed in a continuous, disordered phase ^{1, 2, 3, 4}.
The domains interact with each other ^{5, 6, 7}, and these interactions affect their own movement ^{8, 9} as well as that of other species present in the monolayer ^{10, 11}.
Inter–domain interaction may be related to electrostatic forces (dipolar or Coulombic repulsions), forces related to the spontaneous curvature of the coexisting phases and hydrodynamic forces that appear when domains are in motion.

Dipolar inter–domain interaction, which is always present, arises from the excess dipolar density of the ordered phase with respect to the continuous phase. Hence, the dipolar strength is proportional to the domain area.

Lipid monolayers typically exhibit domain size polydispersity
^{12, 13, 14, 15}.
In particular, Langmuir monolayers at the air–water interface show a wide domain size distribution ^{10, 9, 16, 14}.
Due to the origin of the dipolar interaction, size polydispersity leads to interaction polydispersity, which usually turns out to be the most important, as a consequence of the quadratic dependence of dipolar strength on domain size.

Structural and dynamical properties of monolayers are mainly determined by inter–domain interactions and hence, for systems where the dipolar interaction is dominant, are expected to be strongly affected by the presence of domains of different sizes.

To the best of our knowledge, the effects of size polydispersity on the structural and dynamical properties of lipid monolayers have not been studied. However, for the determination of physical parameters of the constituting lipids, size distribution has been taken into account.
Mulder ^{12} studied the use of size distributions of circular domains in Langmuir monolayers for determining physical parameters of surfactants. He approximates the exact size distribution by a Gaussian distribution and uses a simplified theoretical analysis, in which the inter–domain interactions are treated approximately by considering equally sized domains arranged in a regular hexagonal array.
Lee et al. ^{14} obtain the excess dipolar density by fitting the size distribution with an equilibrium thermodynamic expression. Their scheme assumes no interactions between domains, and hence it is valid for sufficiently diluted systems.

Size distribution effects, on the other hand, have long been studied in colloidal systems, mainly
for static properties ^{17, 18, 19},
phase transitions ^{20, 21, 22, 23, 24},
crystallization ^{25, 26, 27, 28, 29, 30, 31, 32},
glass transition ^{33, 34, 35},
self–assembly ^{36},
drying of colloidal dispersions ^{37},
determination of effective interactions ^{38, 39, 40, 41}
and dynamical properties ^{19}.
In particular, for quasi–two–dimensional systems, most of the works consider binary mixture ^{17, 18, 32}, while a few explore polydisperse colloidal suspensions ^{40, 41, 38, 25}.
These study the influence of polydispersity in the determination of an effective interaction potential, when the polydisperse system is regarded as a monodisperse system.

In a previous work ^{16}, we proposed a novel way of estimating dipolar repulsion, using a passive method that involves the analysis of images of the monolayer with phase coexistence.
The method is based on comparing the pair correlation function obtained from experiments with that obtained from simulations of systems of monodisperse domains interacting by a dipolar density pair potential.
We also studied the point dipole approximation for the dipolar density pair potential, and determined an effective point dipole interaction strength that reproduces the structural properties of a system with dipolar density pair potential.

In this work, we use Brownian dynamics simulations to investigate the effects of polydispersity on the structure and dynamics of lipid monolayer models, with parameters chosen within a range of experimental interest.
For this purpose, we study the pair correlation function, the average mean square displacement and the average time–dependent self–diffusion coefficient, as a function of polydispersity and dipolar interaction strength.
In addition, we study how polydispersity would affect determination of dipolar repulsion strength from experimental data, when the method proposed in Ref. ^{16} is used based on a monodisperse model.

## 2 Model and Theory

### 2.1 Polydisperse domain interactions

We consider a monolayer in its two–phase liquid–condensed (LC) and liquid–expanded (LE) coexistence region, where the LC phase forms domains in the LE phase, which occupies the larger area of the monolayer.
Because of the difference in surface densities, the LC domains possess an excess dipole density, , with respect to the surrounding LE phase ^{42}.
This originates dipolar repulsive interactions between the domains.
We model the mixed monolayer as a 2D–dispersion of domains which interact through a dipolar pair potential. Considering only dipole components perpendicularly oriented to the interface, and using the approximation of a point dipole in the center of each domain with an excluded area,
the resulting pair potential between domain and domain can be described by:

(1) |

where is a hard core repulsive potential, and

(2) |

where is the dipole moment of the domain representing the dipole density over the domain area , ,
center–to–center domain distance,
is the vacuum permittivity and is an effective permittivity ^{43} that considers the relative permittivities of the membrane, , water, , and air, :

(3) |

We consider polydisperse circular domains of radius , accordingly Eq. 2 becomes:

(4) |

This equation makes explicit how domain size distribution leads to interaction polydispersity.

For convenience, we define a dimensionless interaction strength:

(5) |

where is the Boltzmann’s constant, the absolute temperature and the mean value of the radii distribution. Then, the dimensionless dipolar pair potential takes the form:

(6) |

### 2.2 Domain size distribution

Most experiments on lipid monolayers present domain size polydispersity ^{13, 14, 8, 16, 4}. The functional form of the domain size distribution is clearly crucial to the structural and dynamical properties of the monolayer. Here, we approximate the size polydispersity by a truncated normal distribution function,

(7) |

where is the mean or expectation of the distribution, is the standard deviation, and is a normalization constant which arises from simetrically truncating the normal distribution function to exclude negative domain radii. The system polydispersity, , is characterized by the ratio of the width of the distribution to its mean; .

In our simulations we use a discrete counterpart of (see Fig. 1), where we chose a number of domain species , and a bin width , such that the radii distribution of the system is described by the set , , where is the number of domains of type and is the total number of domains. For the system with the largest polydispersity studied, , of the domain radii fall outside the histograms.

### 2.3 Radial distribution function

A key quantity for characterizing the structure of the monolayer is the radial distribution function (RDF), . Considering a distribution of domains in the monolayer plane, is related to the probability of finding a domain at a distance from another domain chosen as a reference point:

(8) |

where are the partial radial distribution functions, defined as:

(9) |

with the total monolayer area and the angular brackets indicating an equilibrium ensemble average. Note that is the probability density of finding a particle at a distance from an particle, where is the number density of domains with radius .

### 2.4 Diffusion

In order to evaluate the effects of polydispersity on the system dynamics, we studied the mean square displacement (MSD) and self–diffusion of domains. The mean square displacement of a domain of type with its center at position at time is given by:

(10) |

where the angular brackets indicate an equilibrium ensemble average. However, many experimental results do not differentiate between domain radii. Hence, we evaluate the average MSD as a representative quantity:

(11) |

where is the molar fraction of domains of type .

The time–dependent self–diffusion coefficient, , is defined as the time derivative of . The short–time limit of corresponds to the free domain diffusion coefficient, which in this work is approximated by the diffusion coefficient of a disk in a two–dimensional simple fluid:

(12) |

where is the viscosity of the fluid. It is important to note that, in our simulations, only the ratios between the different are relevant.

The long–time self–diffusion coefficient is defined as:

(13) |

Analogously, the average time–dependent self–diffusion coefficient, , is obtained from the time derivative of , and the corresponding limits are:

(14) |

(15) |

### 2.5 Simulations

We model the mixed monolayer as a two–dimensional Brownian suspension of interacting hard disks with polydisperse radii immersed in an effective fluid, each disk representing an idealized lipid domain.
The inter–domain interactions are described by the point dipole pair interaction Eq. (6) plus a hard disk repulsive part.
Hydrodynamics interactions are disregarded.
To study the static and dynamical properties of the mixed monolayer model, we performed Brownian dynamics (BD) simulations.
In this scheme, the finite difference equation describing the in-plane displacement of Brownian disks immersed in a fluid during the time step is given by ^{44}

(16) |

where is the position of domain of type at time , is the force on disk of type due to the disk of type and a random displacement vector of domain of type originating from solvent particle collisions. is sampled from a Gaussian distribution with zero mean and covariance matrix:

(17) |

where is the identity matrix, and the Kronecker delta.

The simulated systems consisted of disks with radii distribution under periodic boundary conditions, using the minimum image convention.

The size of the simulation box, , was determined from the condensed area fraction, , defined as:

(18) |

The system is completely characterized by the following parameters: the interaction strength (or equivalently, the dipolar density ), the size distribution and the total area fraction .

Throughout this work, the area fraction is fixed at the value , which was chosen as a typical value of experimental monolayer micrographs.

In Figure 1, we show the distributions used in our studies and, in Table S1 of the electronic supplementary information (ESI), the size distribution of each system studied can be found. We have carried out simulation studies with four distinct polydispersities: and . These systems are compared with simulations of perfectly monodisperse domains.

The bin width was chosen as . This choice leads to a small number of domain types for each polydispersity, and at the same time reproduces the distribution shape qualitatively well. For a typical experimental average radius of px, the bin width px is in the range of typical optical microscopy experimental error. We verified that, for the systems studied, the Gaussian distribution is well–described by the selected discretization . The RDF and the average MSD do not change substantially if half of the bin width value is used.

Finally, the interaction strength was varied such that the system remains in its fluid phase and has experimental interest.

The time scale used in the simulations is . In our simulations, we used and for monodisperse, and , respectively. These values were selected so that the dynamics of the smallest domain in each system is well–resolved. For all systems we used domains. The number of domains of each type for the systems considered are specified in Table S1 of the ESI. For the systems studied, we verified that there is no system size dependency.

## 3 Results and Discussions

### 3.1 Effects of polydispersity on and

With the aim of analyzing the effect of polydispersity, we consider systems with four different interaction strengths, which correspond to the liquid regime for the selected area fraction (), and for each interaction strength we vary the polydispersity from (monodisperse) to . Here, we present the results for and , and, in the ESI, the intermediate values of and are shown (Fig. S1).

Figure 2 shows the pair correlation function for (a) and (b), for and . In both cases, we observe a decrease and a broadening of the first peak as the polydispersity grows. Furthermore, there is a shift of the peak position to greater distances and the start to show correlations for shorter distances. In general, polydispersity softens the peaks and minima of the . It is remarkable that, even for of polydispersity, there is still a well–developed first minimum and second maximum, for highly interacting systems. This can be attributed to the fact that the interaction strength grows with the fourth power of the domain size.

To analyze how polydispersity affects the dynamical properties of the monolayer, we calculated the average MSD and the average time–dependent diffusion coefficient. Results for the selected systems are presented in Fig. 3. For the weak interaction (Fig. 3a), we observe that the average short–time diffusion coefficients of polydisperse systems are larger than that of monodisperse system. This is a direct consequence of the inverse radius dependence of the short–time coefficient diffusion, , and of the symmetry of the radii distribution. Here, for the more polydisperse case, , the difference is about . For the average long–time diffusion coefficients we find a similar ordering, but with slightly larger differences. In particular, for the difference reaches . However, the relative decrease of the average long–time diffusion constant with respect to the short–time limit, , differs by less than . For all polydispersities, diminish roughly at the long–time limit. In the inset of Fig. 3a, the respective average MSDs are shown. Note that, within the logarithmic scale, the polydisperse systems appear very similar to the monodisperse system, except for a slight translation to larger values.

Figure 3b shows the average diffusion quantities for . For these systems, the same qualitative behavior is observed as in previous ones. In particular, the short–time limits are identical since they have identical radii distribution. In addition, the average long–time diffusion coefficients are smaller, because of the stronger interactions.

Comparing the most polydisperse system, , with the monodisperse one, it is found that the average long–time diffusion coefficient is of the order of larger, while the ratio is roughly larger. In general, an enhancement of average diffusion is observed as polydispersity increases. This effect is stronger in the long–time regime and is also more pronounced for systems with stronger interactions.

To describe the structure of the polydisperse system in more detail, we show, for the system with and , the partial radial distribution functions in Fig. 4a and in Fig. 4b. These show that the partial RDFs start to differentiate from zero at distances larger than the respective contact values, i.e., domains do not come into contact. Note that the distance between the first and the second neighbor shell and the first minimum depth are very similar for different domain sizes. This is clearly seen in the inset of Fig. 4a, where we have plotted horizontally shifted by their respective peak positions, . Considering the spatial correlation between the smaller domains and the other domain types, as shown in Fig. 4b, it is observed that the probability of finding the smallest domains (type ) around domains of other types increases with domain sizes. However, the first minimum depth of decreases with domain type . The other , not shown here, also indicates that small domains are more probably found close to larger ones.

Focusing now on the time dependent diffusivities, in Fig. 5 we show the time dependent diffusion coefficient for each domain type, , and the corresponding average, . Here, to compare the relative slowdown of the dynamics between the different types, the diffusivities are normalized by their corresponding short time limit, . It is observed that, except for the smallest domain type, the self–diffusion of the different domain sizes slows down similarly, reaching a long–time limit of roughly . There is a remarkably large difference in the behavior of the smallest domains, which show a much smaller slowdown than the others, i.e., they are able to leave the neighbor cages more easily. This could be attributed to the fact that the dipolar strength of the domains grows quadratically with the radius.

### 3.2 Monodisperse vs. polydisperse models

In lipid monolayers with phase coexistence, the dipolar repulsion () is usually not a known parameter.
One method to determine it is to fit the experimentally measured RDF with one obtained from simulations, using as the only adjustable parameter^{16}.
This assumes monodisperse distribution of domain radii, and hence monodisperse inter–domain interactions.
Therefore, it is important to assess the effects of polydispersity on the determination of using this method.
For this purpose, we used the from the monodisperse systems studied in this work as the reference RDF, as if they were previously fitted to experimental data sets.
Then, we fitted the reference with the polydisperse systems shown in Fig. 1 (as was already stated for all systems) and we analyzed how the interaction strength varies with .
Note that, as in our previous work ^{16}, only the first peak height of the reference is considered in the fitting procedure.

Figure 6 shows the results for systems with (a) and (b) (monodisperse models with and , respectively). In the ESI, systems with and ( and , respectively) are shown (Fig. S3).

The RDF from the monodisperse system with (Fig. 6a) is qualitatively well captured by the polydisperse systems, except for short distances, where for more polydisperse systems the start to grow at shorter distances, as expected. In the inset, we have plotted the horizontally shifted so that the first peak position coincides. This clearly shows that the distance between the different neighbor shells and the depth of minima are very similar for all polydispersities. This agreement is striking, since we are comparing systems with polydispersities as large as .

For a more structured system with , shown in Figure 6b, we observe similar behavior. However, for polydispersities larger than , the depth of the first minimum starts to differentiate, at the same time as the third peak position begins to dephase. This can be seen in the inset. Note that, for polydispersity, the fit is already not possible, i.e., there is no interaction strength for which the resulting reaches the maximum value of .

In general, to reach a certain value of a larger interaction strength is needed as the polydispersities increase. However, not all the values of can be obtained for systems with large polydispersity. Figure 7 shows as a function of the interaction strength for the different polydispersities studied here. It can be seen that the values tend to saturate for high interaction strength.

For experimental systems with and approximately Gaussian size distributions, Figure 7 can be used as a working curve to estimate (or ) directly from the experimental without implementing any simulation, generalizing the method presented in Ref. ^{16} to polydisperse systems.

At this point, we introduce a new parameter , which combines the number density, , and the interaction strength, , in one independent dimensionless parameter. Namely,

(19) |

where is the mean geometrical distance between domains.

Systems for which the hard disk interactions can be disregarded (i.e., strongly interacting or low density systems) are completely determined by , the scaled domain radii of the different species, , and their corresponding molar fractions, . Note that the parameter space is mapped to .

Given that, for the systems studied here, the hard disk interaction turned out to be irrelevant, Fig. 7 may be recast to show the maximum of the RDF as a function of the , as shown in Fig. 8. In this way, if the size polydispersity can be described by a Gaussian distribution with , and and are determined from the experiments, the value could be estimated directly from this figure, and subsequently the dipolar density can be obtained.

To further study how polydispersity affects the structure, Figure 9 shows the interaction strength needed to obtain a certain as a function of . The curves for and correspond to the RDFs shown in Fig. 6 a and b, respectively. Here, it can be clearly seen that the interaction strength for a given increases notably faster for more structured systems and, in particular, for and , and for and . The dotted line for the case indicates that no system with is able to reach this peak height.

Finally, we consider the dynamics of the systems studied shown in Figure 6. The corresponding average diffusion coefficient and MSD are shown in Figure 10.

It is remarkable that, for the less structured system (Fig. 10a), polydispersity does not much affect the intermediate and long–time average dynamical quantities. However, a completely different scenario occurs for the more structured system (Fig. 10b). In this case, up to , the average diffusion coefficient behaves similarly to that of the monodisperse system. On the other hand, for larger polydispersities, strongly deviates from the monodisperse system. The system reaches the sub–diffusive regime faster, as polydispersity increases. Besides, the long–time average diffusion coefficient is smaller for larger polydispersities. This is probably a consequence of the dependence of the interaction strength on the domain sizes, which is evident for strongly interacting systems.

## 4 Conclusions

In this work, we studied the influence of domain–size polydispersity on the structure and dynamics of model lipid monolayers with phase coexistence at the air–water interface. The size–polydispersity was modeled as a discretized Gaussian distribution.

Studying the monolayer structure, we found a decrease and a broadening of the first peak of the RDF as the polydispersity grows. Furthermore, a shift of the peak position to greater distances and the occurrence of correlations for shorter distances were observed. Notably, highly interacting systems () presented a well–developed first minimum and second maximum, even for of polydispersity. For all the systems studied, the partial RDF starts to differentiate from zero at distances larger than the respective contact values. This indicates that domains do not come into contact and hence that the hard–disk interaction is not relevant. Regarding the spatial correlation between the smaller domains and the other domain types, it was observed that the probability of finding the smallest domains around domains of other types increases with domain sizes.

Analyzing the domain dynamics, we found an enhancement of the average diffusion as polydispersity increases. This is more pronounced in the long–time regime. For systems with stronger interactions, the overall enhancement is more noticeable. It was also found that the self–diffusion of the smallest domains shows a much smaller slowdown than the other domain sizes, i.e., they are able to leave the neighbor cages more easily.

We also studied the effects of polydispersity in the determination of (or ) by the method proposed in Ref. ^{16}, where a monodisperse model is used to fit the experimental RDF.
For the systems considered, it was found that polydispersity strongly affects the value of obtained, which is greatly underestimated if polydispersity is not considered.
Only for the experimental system with small polydispersities and/or weak interactions is achieved a good approximation.

It is remarkable that, even for large polydispersities, the fitted RDFs result in good agreement with the reference one; they have the same second and third peak heights and neighbor shell distances. They mainly differ in the depth of the first minimum and for distances where the starts to grow.

With regard to the dynamics, on the other hand, only for weak interactions or small polydispersities does the average time–dependent diffusion coefficient agree with the reference system. For stronger interactions or larger polydispersities, a noticeable slow–down is observed in the average dynamics.

Finally, the method proposed in Ref. ^{16} may be straightforwardly generalized to include polydispersity by directly fitting the experimental with a model that accounts for the measured domain size distribution in the simulations.
Alternatively, using a Gaussian size–distribution model and for a selected range of , a set of figures, like Fig. 7, can be generated for different values of , for later use as working curves to estimate (or ) directly from the experimental data without implementing any simulation.
In particular, for a system in which the hard disk interaction can be disregarded, only one working curve is needed (Fig. 8).

### Acknowledgement

The authors acknowledge financial support from CONICET (Consejo Nacional de Investigaciones Científicas y Técnicas, Argentina), SECyT–UNC (Secretaría de Ciencia y Técnica de la Universidad Nacional de Córdoba, Argentina) and FonCyT (Fondo para la Investigación Científica y Tecnológica, Argentina).

## References

- Duncan et al. 2011 S. L. Duncan, I. S. Dalal and R. G. Larson, Biochimica et Biophysica Acta, 2011, 1808, 2450–2465.
- Dhar et al. 2012 P. Dhar, E. Eck, J. N. Israelachvili, D. W. Lee, Y. Min, A. Ramachandran, A. J. Waring and J. A. Zasadzinski, Biophysical Journal, 2012, 102, 56–65.
- Baoukina and Tieleman 2016 S. Baoukina and D. P. Tieleman, Biochimica et Biophysica Acta, 2016, 1858, 2431–2440.
- Min 2017 Y. Min, Current Opinion in Chemical Engineering, 2017, 15, 76–83.
- Andelman et al. 1985 D. Andelman, F. Brochard, P. G. de Gennes and J. F. Joanny, C. R. Acad. Sc. Paris, 1985, 301, 675–679.
- Andelman et al. 1986 D. Andelman, F. Brochard and J. F. Joanny, J. Chem. Phys., 1986, 86, 3673–3681.
- Ursell and W. S. Klug 2009 T. S. Ursell and R. P. W. S. Klug, Proc. Natl. Acad. Sci., 2009, 106, 13301–13306.
- Caruso et al. 2014 B. Caruso, M. Villareal, L. Reinaudi and N. Wilke, J. Phys. Chem. B, 2014, 118, 519–529.
- Wilke et al. 2010 N. Wilke, F. Vega Mercado and B. Maggio, Langmuir, 2010, 26, 11050–11059.
- Forstner et al. 2008 M. Forstner, D. Martin, F. Rückerl, J. A. Käs and C. Selle, Phys. Rev. E, 2008, 77, 051906–1–7.
- Rükerl et al. 2008 F. Rükerl, J. A. Käs and C. Selle, Langmuir, 2008, 24, 3365–3369.
- Mulder 2003 W. H. Mulder, J. Colloid Interface Sci., 2003, 264, 558–560.
- Hu et al. 2006 Y. Hu, K. Meleson and J. Israelachvili, Biophysical Journal, 2006, 91, 444–453.
- Lee et al. 2011 D. W. Lee, Y. Min, P. Dhar, A. Ramachandran, J. N. Israelachvili and J. A. Zasadzinski, Proc. Natl. Acad. Sci. U. S. A., 2011, 108, 9425–9430.
- Rosetti et al. 2017 C. M. Rosetti, A. Mangiarotti and N. Wilke, Biochimicia et Biophysica Acta, 2017, 1859, 789–802.
- Rufeil-Fiori et al. 2016 E. Rufeil-Fiori, N. Wilke and A. J. Banchio, Soft Matter, 2016, 12, 4769–4777.
- Hoffmann et al. 2006 N. Hoffmann, C. N. Likos and H. Löwen, J. Phys.: Condens. Matter, 2006, 18, 10193–10211.
- Assoud et al. 2009 L. Assoud, F. Ebert, P. Keim, R. Messina, G. Maret and H. Löwen, J. Phys.: Condens. Matter, 2009, 21, 464114.
- Yiannourakou et al. 2010 M. Yiannourakou, I. G. Economou and I. A. Bitsanis, J. Chem. Phys., 2010, 133, 224901.
- Yiannourakou et al. 2009 M. Yiannourakou, I. G. Economou and I. A. Bitsanis, J. Chem. Phys., 2009, 130, 194902.
- Fasolo and Sollich 2003 M. Fasolo and P. Sollich, Phys. Rev. Lett., 2003, 91, 068301.
- van der Linden et al. 2013 M. N. van der Linden, A. van Blaaderen and M. Dijkstra, J. Chem. Phys., 2013, 138, 114903.
- van den Pol et al. 2008 E. van den Pol, D. M. E. Thies-Weesie, A. V. Petukhov, G. J. Vroege and K. Kvashnina, J. Chem. Phys., 2008, 129, 164715.
- Sollich and Wilding 2011 P. Sollich and N. B. Wilding, Soft Matter, 2011, 7, 4472.
- Tretiakov and Wojciechowski 2012 K. V. Tretiakov and K. W. Wojciechowski, J. Chem. Phys., 2012, 136, 204506.
- Iacopini et al. 2009 S. Iacopini, T. Palberg and H. J. Schöpe, Phys. Rev. E, 2009, 79, 010601.
- Sear 1998 R. P. Sear, Europhys. Lett., 1998, 44, 531.
- Pusey 1987 P. N. Pusey, J. Phys. (Les Ulis, Fr.), 1987, 48, 709.
- McRae and Haymet 1988 R. McRae and A. D. J. Haymet, J. Chem. Phys., 1988, 88, 1114.
- Bolhuis and Kofke 1996 P. G. Bolhuis and D. A. Kofke, Phys. Rev. E, 1996, 54, 634.
- Auer and Frenkel 2001 S. Auer and D. Frenkel, Nature, 2001, 413, 711–713.
- Löwen et al. 2013 H. Löwen, T. Horn, T. Neuhaus and B. ten Hagen, Eur. Phys. J. Special Topics, 2013, 222, 2961–2972.
- Zaccarelli et al. 2015 E. Zaccarelli, S. M. Liddle and W. C. K. Poon, Soft Matter, 2015, 11, 324.
- Voigtmann 2003 T. Voigtmann, Phys. Rev. E, 2003, 68, 051401.
- Heckendorf et al. 2017 D. Heckendorf, K. J. Mutch, S. U. Egelhaaf and M. Laurati, Phys. Rev. Lett., 2017, 119, 048003.
- Cabane et al. 2016 B. Cabane, J. Li, F. Artzner, R. Botet, C. Labbez, G. Bareigts, M. Sztucki and L. Goehring, Phys. Rev. Lett., 2016, 116, 208001.
- Boulogne et al. 2016 F. Boulogne, Y. L. Kong, J. K. Nunes and H. A. Stone, Phys. Rev. Lett., 2016, 116, 238001.
- Frydel and Rice 2005 D. Frydel and S. A. Rice, Phys. Rev. E, 2005, 71, 041403.
- Pond et al. 2011 M. J. Pond, J. R. Errington and T. M. Truskett, J. Chem. Phys., 2011, 135, 124513–1–9.
- Pangburn and Bevan 2005 T. O. Pangburn and M. A. Bevan, J. Chem. Phys., 2005, 123, 174904.
- Pangburn and Bevan 2006 T. O. Pangburn and M. A. Bevan, J. Chem. Phys., 2006, 124, 054712.
- McConnell 1991 H. M. McConnell, Annu. Rev. Phys. Chem., 1991, 42, 171–195.
- Urbakh and Klafter 1993 M. Urbakh and J. Klafter, J. Phys. Chem., 1993, 97, 3344–3349.
- Ermak and McCammon 1978 D. Ermak and J. A. McCammon, J. Chem. Phys., 1978, 69, 1352–1360.