Lamellar order, microphase structures, and glassy phase in a field theoretic model for charged colloids

Non-equilibrium view of the amorphous solidification of liquids with competing interactions

The interplay between short-range attractions and long-range repulsions (SALR) characterizes the so-called liquids with competing interactions, which are known to exhibit a variety of equilibrium and non-equilibrium phases. The theoretical description of the phenomenology associated with glassy or gel states in these systems has to take into account both the presence of thermodynamic instabilities (such as those defining the spinodal line and the so called λ line) and the limited capability to describe genuine non-equilibrium processes from first principles. Here, we report the first application of the non-equilibrium self-consistent generalized Langevin equation theory to the description of the dynamical arrest processes that occur in SALR systems after being instantaneously quenched into a state point in the regions of thermodynamic instability. The physical scenario predicted by this theory reveals an amazing interplay between the thermodynamically driven instabilities, favoring equilibrium macro- and micro-phase separation, and the kinetic arrest mechanisms, favoring non-equilibrium amorphous solidification of the liquid into an unexpected variety of glass and gel states.

Competition between clustering and phase separation in binary mixtures containing SALR particles

We investigate by Monte Carlo simulations a mixture of particles with competing interactions (hard-sphere two-Yukawa, HSTY) and hard spheres (HS), with same diameters σ and a square-well (SW) cross attraction. In a recent study [G. Munaò et al., J. Phys. Chem. B, 2022, 126, 2027-2039], we have analysed situations-in terms of relative concentration and attraction strength-where HS promote the formation of clusters involving particles of both species under thermodynamic conditions that would not allow for clustering of the pure HSTY fluid. Here, we focus on the role played by the range of cross attraction in determining the equilibrium structure of the mixture, starting from a homogeneous low-density state. When the width of the well exceeds approximately σ, clustering takes place in the system, with aggregates characterised by various sizes and shapes. Only for low HSTY concentrations (less than 10%) a single big cluster appears, anticipating the behaviour observed for a wider well, around 1.2σ. In the latter case, a spherical cluster encompassing almost all particles is the stable structure at equilibrium. We interpret this outcome as a macrophase, liquid-vapour separation where the spherical cluster is just the form taken at low density by the liquid phase inside the vapour phase: indeed, when the density takes larger values, periodic boundary conditions select liquid-vapour interfaces with other non-spherical shapes, similarly as found for a finite sample of simple fluid going through the liquid-vapour coexistence region. For still higher densities we document the existence of a solid phase characterized by the alternation of bilayers filled with particles of one species and bilayers of the other species, giving the solid a peculiar wafer structure.

Hierarchical microphase separation in non-conserved active mixtures

Non-equilibrium phase separating systems with reactions, such as biomolecular condensates and bacteria colonies, can break time-reversal symmetry (TRS) in two distinct ways. Firstly, the conservative and non-conservative sectors of the dynamics can be governed by incompatible free energies; when both sectors are present, this is the leading-order TRS violation, captured in its simplest form by 'Model AB'. Second, the diffusive dynamics can break TRS in its own right. This happens only at higher order in the gradient expansion (but is the leading behaviour without reactions present) and is captured by 'Active Model B+' (AMB+). Each of the two mechanisms can lead to microphase separation, by quite different routes. Here we introduce Model AB+, for which both mechanisms are simultaneously present, and show that for slow reaction rates the system can undergo a new type of hierarchical microphase separation, which we call 'bubbly microphase separation'. In this state, small droplets of one fluid are continuously created and absorbed into large droplets, whose length-scales are controlled by the competing reactive and diffusive dynamics.

Solution of disordered microphases in the Bethe approximation

The periodic microphases that self-assemble in systems with competing short-range attractive and long-range repulsive (SALR) interactions are structurally both rich and elegant. Significant theoretical and computational efforts have thus been dedicated to untangling their properties. By contrast, disordered microphases, which are structurally just as rich but nowhere near as elegant, have not been as carefully considered. Part of the difficulty is that simple mean-field descriptions make a homogeneity assumption that washes away all of their structural features. Here, we study disordered microphases by exactly solving a SALR model on the Bethe lattice. By sidestepping the homogenization assumption, this treatment recapitulates many of the key structural regimes of disordered microphases, including particle and void cluster fluids as well as gelation. This analysis also provides physical insight into the relationship between various structural and thermal observables, between criticality and physical percolation, and between glassiness and microphase ordering.

Hunting mermaids in real space: known knowns, known unknowns and unknown unknowns

We review efforts to realise so-called mermaid (or short-ranged attraction/long ranged repulsion) interactions in 3d real space. The repulsive and attractive contributions to these interactions in charged colloids and colloid-polymer mixtures, may be accurately realised, by comparing particle-resolved studies with colloids to computer simulation. However, when we review work where these interactions have been combined, despite early indications of behaviour consistent with predictions, closer analysis reveals that in the non-aqueous systems used for particle-resolved studies, the idea of summing the attractive and repulsive components leads to wild deviations with experiment. We suggest that the origin lies in the weak ion dissociation in these systems with low dielectric constant solvents. Ultimately this leads even to non-centro-symmetric interactions and a new level of complexity in these systems.

Clustering and assembly dynamics of a one-dimensional microphase former

Both ordered and disordered microphases ubiquitously form in suspensions of particles that interact through competing short-range attraction and long-range repulsion (SALR). While ordered microphases are more appealing materials targets, understanding the rich structural and dynamical properties of their disordered counterparts is essential to controlling their mesoscale assembly. Here, we study the disordered regime of a one-dimensional (1D) SALR model, whose simplicity enables detailed analysis by transfer matrices and Monte Carlo simulations. We first characterize the signature of the clustering process on macroscopic observables, and then assess the equilibration dynamics of various simulation algorithms. We notably find that cluster moves markedly accelerate the mixing time, but that event chains are of limited help in the clustering regime. These insights will inspire further study of three-dimensional microphase formers.

Equilibrium Phase Behavior of a Continuous-Space Microphase Former

Periodic microphases universally emerge in systems for which short-range interparticle attraction is frustrated by long-range repulsion. The morphological richness of these phases makes them desirable material targets, but our relatively coarse understanding of even simple models hinders controlling their assembly. We report here the solution of the equilibrium phase behavior of a microscopic microphase former through specialized Monte Carlo simulations. The results for cluster crystal, cylindrical, double gyroid, and lamellar ordering qualitatively agree with a Landau-type free energy description and reveal the nontrivial interplay between cluster, gel, and microphase formation.

Cooperativity flows and shear-bandings: a statistical field theory approach

Cooperativity effects have been proposed to explain the non-local rheology in the dynamics of soft jammed systems. Based on the analysis of the free-energy model proposed by L. Bocquet, A. Colin and A. Ajdari, Phys. Rev. Lett., 2009, 103, 036001, we show that cooperativity effects resulting from the non-local nature of the fluidity (inverse viscosity) are intimately related to the emergence of shear-banding configurations. This connection materializes through the onset of inhomogeneous compact solutions (compactons), wherein the fluidity is confined to finite-support subregions of the flow and strictly zero elsewhere. The compacton coexistence with regions of zero fluidity ("non-flowing vacuum") is shown to be stabilized by the presence of mechanical noise, which ultimately shapes up the equilibrium distribution of the fluidity field, the latter acting as an order parameter for the flow-noflow transitions occurring in the material.

Polydispersity and gelation in concentrated colloids with competing interactions

In colloids with competing short-range attractions and long-range repulsions, microcrystalline gels are experimentally formed under conditions where computer simulations point to a lamellar phase as the ground state. Here, upon applying a low-frequency alternating electric field, we bring the system from an initial gel state to a columnar-like state. While molecular dynamics simulations on monodisperse colloids reveal that a columnar structure spontaneously evolves towards a lamellar phase, the columnar-like state in experiments relaxes back to the initial disordered gel state once the electric field is switched off. Similarly, a columnar phase in molecular dynamics simulations decomposes into finite-size crystalline clusters as the relative polydispersity of the colloids is around 1.0%. We conclude that the experimentally observed melting of the columnar structure is driven by polydispersity. Moreover, further simulations reveal that the critical polydispersity required to destabilize a long-range ordered structure increases with the attraction range, pointing to the possibility of observing periodic structures in experiments if the attraction range is sufficiently long compared to the polydispersity of the colloids.

Phase separation and self-assembly of colloidal dimers with tunable attractive strength: from symmetrical square-wells to Janus dumbbells

We numerically investigate colloidal dimers with asymmetric interaction strengths to study how the interplay between molecular geometry, excluded volume effects and attractive forces determines the overall phase behavior of such systems. Specifically, our model is constituted by two rigidly-connected tangent hard spheres interacting with other particles in the first instance via identical square-well attractions. Then, one of the square-well interactions is progressively weakened, until only the corresponding bare hard-core repulsion survives, giving rise to a "Janus dumbbell" model. We investigate structure, thermodynamics and phase behavior of the model by means of successive umbrella sampling and Monte Carlo simulations. In most of the cases, the system behaves as a standard simple fluid, characterized by a gas-liquid phase separation, for sufficiently low temperatures. In these conditions we observe a remarkable linear scaling of the critical temperature as a function of the interaction strength. But, as the interaction potential approaches the Janus dumbbell limit, we observe the spontaneous formation of self-assembled lamellar structures, preempting the gas-liquid phase separation. Comparison with previous studies allows us to pinpoint the role of the interaction range in controlling the onset of ordered structures and the competition between the formation of these structures and gas-liquid condensation.

Controlling local packing and growth in calcium-silicate-hydrate gels

We investigate the development of gels under out-of-equilibrium conditions, such as calcium-silicate-hydrate (C-S-H) gels that form during cement hydration and are the major factor responsible for cement mechanical strength. We propose a new model and numerical approach to follow the gel formation upon precipitation and aggregation of nano-scale colloidal hydrates, whose effective interactions are consistent with forces measured in experiments at fixed lime concentrations. We use Grand Canonical Monte Carlo to mimic precipitation events during Molecular Dynamics simulations, with their rate corresponding to the hydrate production rate set by the chemical environment. Our results display hydrate precipitation curves that indeed reproduce the acceleration and deceleration regime typically observed in experiments and we are able to correctly capture the effect of lime concentration on the hydration kinetics and the gel morphology. Our analysis of the evolution of the gel morphology indicates that the acceleration is related to the formation of an optimal local crystalline packing that allows for large, elongated aggregates to grow and that is controlled by the underlying thermodynamics. The defects produced during precipitation favor branching and gelation that end up controlling the deceleration. The effects on the mechanical properties of C-S-H gels are also discussed.

Reversible gels of patchy particles: role of the valence

We simulate a binary mixture of colloidal patchy particles with two and three patches, respectively, for several relative concentrations and hence relative average valences. For these limited-valence systems, it is possible to reach low temperatures, where the lifetime of the patch-patch interactions becomes longer than the observation time without encountering phase separation in a colloid-poor (gas) and a colloid rich (liquid) phase. The resulting arrested state is a fully connected long-lived network where particles with three patches provide the branching points connecting chains of two-patch particles. We investigate the effect of the valence on the structural and dynamic properties of the resulting gel and attempt to provide a theoretical description of the formation and of the resulting gel structure based on a combination of the Wertheim theory for associated liquids and the Flory-Stockmayer approach for modeling chemical gelation.

Self-tuning phase separation in a model with competing interactions inspired by biological cell polarization

We present a theoretical study of a system with competing short-range ferromagnetic attraction and a long-range antiferromagnetic repulsion, in the presence of a uniform external magnetic field. The interplay between these interactions, at sufficiently low temperature, leads to the self-tuning of the magnetization to a value which triggers phase coexistence, even in the presence of the external field. The investigation of this phenomenon is performed using a Ginzburg-Landau functional in the limit of an infinite number of order parameter components (large N model). The scalar version of the model is expected to describe the phase separation taking place on a cell surface when this is immersed in a uniform concentration of chemical stimulant. A phase diagram is obtained as a function of the external field and the intensity of the long-range repulsion. The time evolution of the order parameter and of the structure factor in a relaxation process is studied in different regions of the phase diagram.

Two-dimensional fluid with competing interactions exhibiting microphase separation: theory for bulk and interfacial properties

Colloidal particles that are confined to an interface such as the air-water interface are an example of a two-dimensional fluid. Such dispersions have been observed to spontaneously form cluster and stripe morphologies in certain systems with isotropic pair potentials between the particles, due to the fact that the pair interaction between the colloids has competing attraction and repulsion over different length scales. Here we present a simple density functional theory for a model of such a two-dimensional fluid. The theory predicts a bulk phase diagram exhibiting cluster, stripe, and bubble modulated phases, in addition to homogeneous fluid phases. Comparing with simulation results for this model from the literature, we find that the theory is qualitatively reliable. The model allows for a detailed investigation of the structure of the fluid and we are able to obtain simple approximate expressions for the static structure factor and for the length scale characterizing the modulations in the microphase separated phases. We also investigate the behavior of the system under confinement between two parallel hard walls. We find that the confined fluid phase behavior can be rather complex.

Rheology of colloidal microphases in a model with competing interactions

We study the rheological properties of colloidal microphases in two dimensions simulating a model of colloidal particles with competing interactions. Due to the competition between short-range attraction and long-range repulsion, as a function of the density the model exhibits a variety of microphases such as clusters, stripes, or crystals with bubbles. We prepare the system in a confined microphase employing Monte Carlo simulations and then shear the resulting configurations by applying a drag force profile. We integrate numerically the equation of motion for the particles and analyze the dynamics as a function of the density and the applied strain rate. We measure the stress-strain curves and characterize the yielding of the colloidal microphases. The results depend on the type of microphase. (i) Clusters are easily sheared along layers and the relative motion is assisted by rotations. (ii) Stripes shear easily when they are parallel to the flow and tend to jam when they are perpendicular to it. Under a sufficiently strong shear rate perpendicular stripes orient in the flow direction. (iii) Crystals with bubbles yield by fracturing along the bubbles and eventually forming stripes. We discuss the role of dislocations, emitted by the bubbles, in the yielding process. Finally, we analyze the effect of thermal fluctuations on the rheological properties.

Phase behavior of a fluid with competing attractive and repulsive interactions

Fluids in which the interparticle potential has a hard core, is attractive at moderate separations, and repulsive at large separations are known to exhibit unusual phase behavior, including stable inhomogeneous phases. Here we report a joint simulation and theoretical study of such a fluid, focusing on the relationship between the liquid-vapor transition line and any new phases. The phase diagram is studied as a function of the amplitude of the attraction for a certain fixed amplitude of the long ranged repulsion. We find that the effect of the repulsion is to substitute the liquid-vapor critical point and a portion of the associated liquid-vapor transition line, by two first-order transitions. One of these transitions separates the vapor from a fluid of spherical liquidlike clusters; the other separates the liquid from a fluid of spherical voids. At low temperature, the two transition lines intersect one another and a vapor-liquid transition line at a triple point. While most integral equation theories are unable to describe the new phase transitions, the Percus-Yevick approximation does succeed in capturing the vapor-cluster transition, as well as aspects of the structure of the cluster fluid, in reasonable agreement with the simulation results.