Generalization of Wertheim's theory for the assembly of various types of rings

Breakdown of the law of rectilinear diameter and related surprises in the liquid-vapor coexistence in systems of patchy particles

The phase diagram of molecular or colloidal systems depends strongly on the range and angular dependence of the interactions between the constituent particles. For instance, it is well known that the critical density of particles with "patchy" interactions shifts to lower values as the number of patches is decreased [see Bianchi et al. Phys. Rev. Lett. 97, 168301 (2006)]. Here, we present simulations that show that the phase behavior of patchy particles is even more interesting than had been appreciated. In particular, we find that, upon cooling below the critical point, the width of the liquid-vapor coexistence region of a system of particles with tetrahedrally arranged patches first increases, then decreases, and finally increases again. In other words, this system exhibits a doubly re-entrant liquid-vapor transition. As a consequence, the system exhibits a very large deviation from the law of rectilinear diameter, which assumes that the critical density can be obtained by linear extrapolation of the averages of the densities of the coexisting liquid and vapor phases. We argue that the unusual behavior of this system has the same origin as the density maximum in liquid water and is not captured by the Wertheim theory. The Wertheim theory also cannot account for our observation that the phase diagram of particles with three patches depends strongly on the geometrical distribution of the patches and on the degree to which their position on the particle surface is rigidly constrained. However, the phase diagram is less sensitive to small angular spreads in the patch locations. We argue that the phase behavior reported in this paper should be observable in experiments on patchy colloids and may be relevant for the liquid-liquid equilibrium in solutions of properly functionalized dendrimers.

How to simulate patchy particles⋆

Patchy particles is the name given to a large class of systems of mesoscopic particles characterized by a repulsive core and a discrete number of short-range and highly directional interaction sites. Numerical simulations have contributed significantly to our understanding of the behaviour of patchy particles, but, although simple in principle, advanced simulation techniques are often required to sample the low temperatures and long time-scales associated with their self-assembly behaviour. In this work we review the most popular simulation techniques that have been used to study patchy particles, with a special focus on Monte Carlo methods. We cover many of the tools required to simulate patchy systems, from interaction potentials to biased moves, cluster moves, and free-energy methods. The review is complemented by an educationally oriented Monte Carlo computer code that implements all the techniques described in the text to simulate a well-known tetrahedral patchy particle model.

Valence, loop formation and universality in self-assembling patchy particles

Patchy particles have emerged as an attractive model to mimic phase separation and self-assembly of globular proteins solutions, colloidal patchy particles, and molecular fluids where directional interactions are operative. In our previous work, we extensively explored the coupling of directional and isotropic interactions on both the phase separation and self-assembly in a system of patchy particles with five spots. Here, we extend this work to consider different patch valences and isotropic interaction strengths with an emphasis on self-assembly. Although the location of self-assembly transition lines in the temperature-density plane depend on a number of parameters, we find universal behavior of cluster size that is dependent only on the probability of a spot being bound, the patch valence, and the density. Using these principles, we quantify both the mass distribution and the shape for all clusters, as well as clusters containing loops. Following the logical implications of these results, combined with a simplified version of a mean-field theory that incorporates Flory-Stockmayer theory, we find a universal curve for the temperature dependence of cluster mass and a universal curve for the fraction of clusters that contain loops. As the curves are dependent on the particle valence, such results provide a method for parameterizing patchy particle models using experimental data.

A theory for the effect of patch/non-patch attractions on the self-assembly of patchy colloids

In this paper, we develop a thermodynamic perturbation theory to describe the self-assembly of patchy colloids which exhibit both patch-patch attractions as well as patch/non-patch attractions. That is, patches attract other patches as well as the no patch region (we call this region Ψ). In general, the patch-patch and patch-Ψ attractions operate on different energy scales allowing for a competition between different modes of attraction. This competition may result in anomalous thermodynamic properties. As an application, we tune the patch parameters to reproduce the liquid density (suitably scaled) maximum of water. It is then shown that the liquid branch of the colloids phase diagram has liquid densities consistent with both saturated and super-cooled liquid water. Finally, it is shown that the colloids reproduce water's anomalous minimum in isothermal compressibility and negative volume expansivity.

Limiting the valence: advancements and new perspectives on patchy colloids, soft functionalized nanoparticles and biomolecules

Artistic representation of limited valance units consisting of a soft core (in blue) and a small number of flexible bonding patches (in orange).

Temperature (de)activated patchy colloidal particles

We present a new model of patchy particles in which the interaction sites can be activated or deactivated by varying the temperature of the system. We study the thermodynamics of the system by means of Wertheim's first order perturbation theory, and use Flory-Stockmayer theory of polymerization to analyse the percolation threshold. We find a very rich phase behaviour including lower critical points and reentrant percolation.