Phase behavior of polydisperse spheres: Simulation strategies and an application to the freezing transition

Emerging exotic compositional order on approaching low-temperature equilibrium glasses

The ultimate fate of a glass former upon cooling has been a fundamental problem in condensed matter physics and materials science since Kauzmann. Recently, this problem has been challenged by a model with an extraordinary glass-forming ability effectively free from crystallisation and phase separation, two well-known fates of most glass formers, combined with a particle-size swap method. Thus, this system is expected to approach the ideal glass state if it exists. However, we discover exotic compositional order as the coexistence of space-spanning network-like structures formed by small-large particle connections and patches formed by medium-size particles at low temperatures. Therefore, the glass transition is accompanied unexpectedly by exotic compositional ordering inaccessible through ordinary structural or thermodynamic characterisations. Such exotic compositional ordering is found to have an unusual impact on structural relaxation dynamics. Our study thus raises fundamental questions concerning the role of unconventional structural ordering in understanding glass transition.

Packing Polydisperse Colloids into Crystals: When Charge-Dispersity Matters

Monte Carlo simulations, fully constrained by experimental parameters, are found to agree well with a measured phase diagram of aqueous dispersions of nanoparticles with a moderate size polydispersity over a broad range of salt concentrations, c_{s}, and volume fractions, ϕ. Upon increasing ϕ, the colloids freeze first into coexisting compact solids then into a body centered cubic phase (bcc) before they melt into a glass forming liquid. The surprising stability of the bcc solid at high ϕ and c_{s} is explained by the interaction (charge) polydispersity and vibrational entropy.

Different routes into the glass state for soft thermo-sensitive colloids

We report an experimental and theoretical investigation of glass formation in soft thermo-sensitive colloids following two different routes: a gradual increase of the particle number density at constant temperature and an increase of the radius in a fixed volume at constant particle number density. Confocal microscopy experiments and the non-equilibrium self-consistent generalized Langevin equation (NE-SCGLE) theory consistently show that the two routes lead to a dynamically comparable state at sufficiently long aging times. However, experiments reveal the presence of moderate but persistent structural differences. Successive cycles of radius decrease and increase lead instead to a reproducible glass state, indicating a suitable route to obtain rejuvenation without using shear fields.

Composition inversion in mixtures of binary colloids and polymer

Understanding the phase behaviour of mixtures continues to pose challenges, even for systems that might be considered "simple." Here, we consider a very simple mixture of two colloidal and one non-adsorbing polymer species, which can be simplified even further to a size-asymmetrical binary mixture, in which the effective colloid-colloid interactions depend on the polymer concentration. We show that this basic system exhibits surprisingly rich phase behaviour. In particular, we enquire whether such a system features only a liquid-vapor phase separation (as in one-component colloid-polymer mixtures) or whether, additionally, liquid-liquid demixing of two colloidal phases can occur. Particle-resolved experiments show demixing-like behaviour, but when combined with bespoke Monte Carlo simulations, this proves illusory, and we reveal that only a single liquid-vapor transition occurs. Progressive migration of the small particles to the liquid phase as the polymer concentration increases gives rise to composition inversion-a maximum in the large particle concentration in the liquid phase. Close to criticality, the density fluctuations are found to be dominated by the larger colloids.

Local order and crystallization of dense polydisperse hard spheres

Computer simulations give precious insight into the microscopic behavior of supercooled liquids and glasses, but their typical time scales are orders of magnitude shorter than the experimentally relevant ones. We recently closed this gap for a class of models of size polydisperse fluids, which we successfully equilibrate beyond laboratory time scales by means of the swap Monte Carlo algorithm. In this contribution, we study the interplay between compositional and geometric local orders in a model of polydisperse hard spheres equilibrated with this algorithm. Local compositional order has a weak state dependence, while local geometric order associated to icosahedral arrangements grows more markedly but only at very high density. We quantify the correlation lengths and the degree of sphericity associated to icosahedral structures and compare these results to those for the Wahnström Lennard-Jones mixture. Finally, we analyze the structure of very dense samples that partially crystallized following a pattern incompatible with conventional fractionation scenarios. The crystal structure has the symmetry of aluminum diboride and involves a subset of small and large particles with size ratio approximately equal to 0.5.

Disappearance of the Hexatic Phase in a Binary Mixture of Hard Disks

Recent studies of melting in hard disks have confirmed the existence of a hexatic phase occurring in a narrow window of density which is separated from the isotropic liquid phase by a first-order transition, and from the solid phase by a continuous transition. However, little is known concerning the melting scenario in mixtures of hard disks. Here we employ tailored Monte Carlo simulations to elucidate the phase behavior of a system of large (l) and small (s) disks with diameter ratio σ_{l}/σ_{s}=1.4. We find that as small disks are introduced to a system of large ones, the stability window of the hexatic phase shrinks progressively until the line of continuous transitions terminates at an end point beyond which melting becomes a first-order liquid-solid transition. This occurs at surprisingly low concentrations of the small disks, c≲1%, emphasizing the fragility of the hexatic phase. We speculate that the change to the melting scenario is a consequence of strong fractionation effects, the nature of which we elucidate.

Effect of Size Polydispersity on the Nature of Lennard-Jones Liquids

Polydisperse fluids are encountered everywhere in biological and industrial processes. These fluids naturally show a rich phenomenology exhibiting fractionation and shifts in critical point and freezing temperatures. We study here the effect of size polydispersity on the basic nature of Lennard-Jones (LJ) liquids, which represent most molecular liquids without hydrogen bonds, via two- and three-dimensional molecular dynamics computer simulations. A single-component liquid constituting spherical particles and interacting via the LJ potential is known to exhibit strong correlations between virial and potential energy equilibrium fluctuations at constant volume. This correlation significantly simplifies the physical description of the liquid, and these liquids are now known as Roskilde-simple (RS) liquids. We show that this simple nature of the single-component LJ liquid is preserved even for very high polydispersities (above 40% polydispersity for the studied uniform distribution). We also investigate isomorphs of moderately polydisperse LJ liquids. Isomorphs are curves in the phase diagram of RS liquids along which structure, dynamics, and some thermodynamic quantities are invariant in dimensionless units. We find that isomorphs are a good approximation even for polydisperse LJ liquids. The theory of isomorphs thus extends readily to size polydisperse fluids and can be used to improve even further the understanding of these intriguing systems.

Solid-liquid transition in polydisperse Lennard-Jones systems

We study melting of a face-centered crystalline solid consisting of polydisperse Lennard-Jones spheres with Gaussian polydispersity in size. The phase diagram reproduces the existence of a nearly temperature invariant terminal polydispersity (δ(t) =/~ 0.11), with no signature of reentrant melting. The absence of reentrant melting can be attributed to the influence of the attractive part of the potential upon melting. We find that at terminal polydispersity the fractional density change approaches zero, which seems to arise from vanishingly small compressibility of the disordered phase. At constant temperature and volume fraction the system undergoes a sharp transition from crystalline solid to the disordered amorphous or fluid state with increasing polydispersity. This has been quantified by second- and third-order rotational invariant bond orientational order, as well as by the average inherent structure energy. The translational order parameter also indicates a similar sharp structural change at δ =/~ 0.09 in case of T(*) = 1.0, φ = 0.58. The free energy calculation further supports the sharp nature of the transition. The third-order rotationally invariant bond order shows that with increasing polydispersity, the local cluster favors a more icosahedral arrangement and the system loses its local crystalline symmetry. Interestingly, the value of structure factor S(k) of the amorphous phase at δ =/~ 0.10 (just beyond the solid-liquid transition density at T(*) = 1) becomes 2.75, which is below the value of 2.85 required for freezing given by the empirical Hansen-Verlet rule of crystallization, well known in the theory of freezing.

Effect of nanoparticle polydispersity on the self-assembly of polymer tethered nanospheres

Recent simulations predict that aggregating nanospheres functionalized with polymer "tethers" can self-assemble to form a cylinder, perforated lamellae, lamellae, and even the double gyroid phase, which are phases also seen in block copolymer and surfactant systems. Nanoparticle size polydispersity is likely to be a characteristic of these systems. If too high, polydispersity may destabilize a phase. Using multiple thermodynamic paths to explore the phase diagram as a function of temperature and polydispersity, we explore the effect of nanosphere size polydispersity on the phase diagram. We show that in the portions of the phase diagram characterized by an icosahedral local nanoparticle packing motif, a low amount of polydispersity lowers the energy and a large amount of polydispersity raises the energy of the system by disrupting the icosahedral packing. In general, regions of the phase diagram characterized by liquid-like icosahedral packing have high terminal polydispersities from 15% to more than 30%. In the regions of the phase diagram characterized by crystalline local packing, polydispersity raises the energy of the system and induces a phase transition from crystalline to liquid-like ordering within the nanosphere rich regions of the microphase. We find the bilayer crystalline lamellae phase has a terminal polydispersity of 6%, but may still be partially crystalline up to 12%.