Arrested states in colloidal fluids with competing interactions: A static replica study.

We present the first systematic application of the integral equation implementation of the replica method to the study of arrested states in fluids with microscopic competing interactions (short-range attractive and long-range repulsive, SALR), as exemplified by the prototype Lennard-Jones-Yukawa model. Using a wide set of potential parameters, we provide as many as 11 different phase diagrams on the density (ρ)-temperature (T) plane, embodying both the cluster-phase boundary, TC(ρ), and the locus below which arrest takes place, TD(ρ). We describe how the interplay between TC and TD-with the former falling on top of the other, or the other way around, depending on thermodynamic conditions and potential parameters-gives rise to a rich variety of non-ergodic states interspersed with ergodic ones, of which both the building blocks are clusters or single particles. In a few cases, we find that the TD locus does not extend all over the density range subtended by the TC envelope; under these conditions, the λ-line is within reach of the cluster fluid, with the ensuing possibility to develop ordered microphases. Whenever a comparison is possible, our predictions favorably agree with previous numerical results. Thereby, we demonstrate the reliability and effectiveness of our scheme to provide a unified theoretical framework for the study of arrested states in SALR fluids, irrespective of their nature.

Supersolid Phases of Bosonic Particles in a Bubble Trap.

Confinement can have a considerable effect on the behavior of particle systems and is therefore an effective way to discover new phenomena. A notable example is a system of identical bosons at low temperature under an external field mimicking an isotropic bubble trap, which constrains the particles to a portion of space close to a spherical surface. Using path integral Monte Carlo simulations, we examine the spatial structure and superfluid fraction in two emblematic cases. First, we look at soft-core bosons, finding the existence of supersolid cluster arrangements with polyhedral symmetry; we show how different numbers of clusters are stabilized depending on the trap radius and the particle mass, and we characterize the temperature behavior of the cluster phases. A detailed comparison with the behavior of classical soft-core particles is provided too. Then, we examine the case, of more immediate experimental interest, of a dipolar condensate on the sphere, demonstrating how a quasi-one-dimensional supersolid of clusters is formed on a great circle for realistic values of density and interaction parameters. Crucially, this supersolid phase is only slightly disturbed by gravity. We argue that the predicted phases can be revealed in magnetic traps with spherical-shell geometry, possibly even in a lab on Earth. Our results pave the way for future simulation studies of correlated quantum systems in curved geometries.

A density functional theory and simulation study of stripe phases in symmetric colloidal mixtures.

In a binary mixture, stripes refer to a one-dimensional periodicity of the composition, namely, a regular alternation of layers filled with particles of mostly one species. We have recently introduced [Munaò et al., Phys. Chem. Chem. Phys. 25, 16227 (2023)] a model that possibly provides the simplest binary mixture endowed with stripe order. The model consists of two species of identical hard spheres with equal concentration, which mutually interact through a square-well potential. In that paper, we have numerically shown that stripes are present in both liquid and solid phases when the attraction range is rather long. Here, we study the phase behavior of the model in terms of a density functional theory capable to account for the existence of stripes in the dense mixture. Our theory is accurate in reproducing the phases of the model, at least insofar as the composition inhomogeneities occur on length scales quite larger than the particle size. Then, using Monte Carlo simulations, we prove the existence of solid stripes even when the square well is much thinner than the particle diameter, making our model more similar to a real colloidal mixture. Finally, when the width of the attractive well is equal to the particle diameter, we observe a different and more complex form of compositional order in the solid, where each species of particle forms a regular porous matrix holding in its holes the other species, witnessing a surprising variety of emergent behaviors for a very basic model of interaction.

Microphase versus macrophase separation in the square-well-linear fluid: A theoretical and computational study.

Due to the presence of competing interactions, the square-well-linear fluid can exhibit either liquid-vapor equilibrium (macrophase separation) or clustering (microphase separation). Here we address the issue of determining the boundary between these two regimes, i.e., the Lifshitz point, expressed in terms of a relationship between the parameters of the model. To this aim, we carry out Monte Carlo simulations to compute the structure factor of the fluid, whose behavior at low wave vectors accurately captures the tendency of the fluid to form aggregates or, alternatively, to phase separate. Specifically, for a number of different combinations of attraction and repulsion ranges, we make the system go across the Lifshitz point by increasing the strength of the repulsion. We use simulation results to benchmark the performance of two theories of fluids, namely, the hypernetted chain (HNC) equation and the analytically solvable random phase approximation (RPA); in particular, the RPA theory is applied with two different prescriptions as for the direct correlation function inside the core. Overall, the HNC theory proves to be an appropriate tool to characterize the fluid structure and the low-wave-vector behavior of the structure factor is consistent with the threshold between microphase and macrophase separation established through simulation. The structural predictions of the RPA theory turn out to be less accurate, but this theory offers the advantage of providing an analytical expression of the Lifshitz point. Compared to simulation, both RPA schemes predict a Lifshitz point that falls within the macrophase-separation region of parameters: in the best case, barriers roughly twice higher than predicted are required to attain clustering conditions.

Like aggregation from unlike attraction: stripes in symmetric mixtures of cross-attracting hard spheres.

Self-assembly of colloidal particles into striped phases is at once a process of relevant technological interest-just think about the possibility to realise photonic crystals with a dielectric structure modulated along a specific direction-and a challenging task, since striped patterns emerge in a variety of conditions, suggesting that the connection between the onset of stripes and the shape of the intermolecular potential is yet to be fully unravelled. Hereby, we devise an elementary mechanism for the formation of stripes in a basic model consisting of a symmetric binary mixture of hard spheres that interact via a square-well cross attraction. Such a model would mimic a colloid in which the interspecies affinity is of longer range and significantly stronger than the intraspecies interaction. For attraction ranges shorter enough than the particle size the mixture behaves like a compositionally-disordered simple fluid. Instead, for wider square-wells, we document by numerical simulations the existence of striped patterns in the solid phase, where layers of particles of one species are interspersed with layers of the other species; increasing the attraction range stabilises the stripes further, in that they also appear in the bulk liquid and become thicker in the crystal. Our results lead to the counterintuitive conclusion that a flat and sufficiently long-ranged unlike attraction promotes the aggregation of like particles into stripes. This finding opens a novel way for the synthesis of colloidal particles with interactions tailored at the development of stripe-modulated structures.

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.

Condensation and Crystal Nucleation in a Lattice Gas with a Realistic Phase Diagram.

We reconsider model II of Orban et al. (J. Chem. Phys. 1968, 49, 1778−1783), a two-dimensional lattice-gas system featuring a crystalline phase and two distinct fluid phases (liquid and vapor). In this system, a particle prevents other particles from occupying sites up to third neighbors on the square lattice, while attracting (with decreasing strength) particles sitting at fourth- or fifth-neighbor sites. To make the model more realistic, we assume a finite repulsion at third-neighbor distance, with the result that a second crystalline phase appears at higher pressures. However, the similarity with real-world substances is only partial: Upon closer inspection, the alleged liquid−vapor transition turns out to be a continuous (albeit sharp) crossover, even near the putative triple point. Closer to the standard picture is instead the freezing transition, as we show by computing the free-energy barrier relative to crystal nucleation from the "liquid".

Clustering in Mixtures of SALR Particles and Hard Spheres with Cross Attraction.

Self-assembling complex fluids are often modeled as particles with effective competing isotropic interactions, combining a short-range attraction (SA) followed by a longer-range repulsion (LR). For moderately low temperatures and densities, SALR particles form clusters in equilibrium, at least provided that the potential parameters are appropriate. Here we inquire into the possibility that cluster formation in SALR fluids might be pushed by a foreign species even under thermodynamic conditions that would not allow for clusterization of the pure system. To this aim, we study by Monte Carlo simulations a mixture of hard-sphere two-Yukawa particles and hard spheres, with a cross interaction modeled by a square-well attraction, and we investigate the conditions of clustering in terms of strength of attraction and relative concentration of the two species. We find that clusters can occur in the mixture for the same temperature and density where the pure SALR fluid is almost structureless. In particular, we single out a cross attraction such that clusters are formed with a SALR concentration as low as 5%. We also find a situation where nearly pure droplets of hard spheres are held together by a shell of SALR particles. Conversely, we show that clustering can be undermined in the mixture under conditions for which this process takes place in the parent SALR fluid. Using a simple criterion, based on the second virial coefficients of the attractive part of interaction potentials (the so-called "reference attractive fluids"), we are able to predict accurately whether clustering is favored (or hindered) in the mixture, as compared to the pure SALR fluid.

Early stages of aggregation in fluid mixtures of dimers and spheres: a theoretical and simulation study.

We use Monte Carlo simulation and the Reference Interaction Site Model (RISM) theory of molecular fluids to investigate a simple model of colloidal mixture consisting of dimers, made up of two tangent hard monomers of different size, and hard spheres. In addition to steric repulsion, the two species interact via a square-well attraction only between small monomers and spheres. Recently, we have characterized the low-temperature regime of this mixture by Monte Carlo, reporting on the spontaneous formation of a wide spectrum of supramolecular aggregates [Prestipino et al., J. Phys. Chem. B, 2019, 123, 9272]. Here we focus on a regime of temperatures where, on cooling, the appearance of local inhomogeneties first, and the early stages of aggregation thereafter, are observed. In particular, we find signatures of aggregation in the onset of a low-wavevector peak in the structure factors of the mixture, as computed by both theory and simulation. Then, we link the structural information to the microscopic arrangement through a detailed cluster analysis of Monte Carlo configurations. In this regard, we devise a novel method to compute the maximum distance for which two spheres can be regarded as bonded together, a crucial issue in the proper identification of fluid aggregates. The RISM theory provides relatively accurate structural and thermodynamic predictions in comparison with Monte Carlo, but with slightly degrading performances as the fluid progresses inside the locally inhomogeneous phase. Our study certifies the efficacy of the RISM approach as a useful complement to numerical simulation for a reasoned analysis of aggregation properties in colloidal mixtures.

Statistical Mechanics and Thermodynamics of Liquids and Crystals.

Thermodynamic phases are the most prominent manifestation of emergent behavior [...].

Self-Assembled Structures of Colloidal Dimers and Disks on a Spherical Surface.

We study self-assembly on a spherical surface of a model for a binary mixture of amphiphilic dimers in the presence of guest particles via Monte Carlo (MC) computer simulation. All particles had a hard core, but one monomer of the dimer also interacted with the guest particle by means of a short-range attractive potential. We observed the formation of aggregates of various shapes as a function of the composition of the mixture and of the size of guest particles. Our MC simulations are a further step towards a microscopic understanding of experiments on colloidal aggregation over curved surfaces, such as oil droplets.

Entropy Multiparticle Correlation Expansion for a Crystal.

As first shown by H. S. Green in 1952, the entropy of a classical fluid of identical particles can be written as a sum of many-particle contributions, each of them being a distinctive functional of all spatial distribution functions up to a given order. By revisiting the combinatorial derivation of the entropy formula, we argue that a similar correlation expansion holds for the entropy of a crystalline system. We discuss how one- and two-body entropies scale with the size of the crystal, and provide fresh numerical data to check the expectation, grounded in theoretical arguments, that both entropies are extensive quantities.

Ultracold Bosons on a Regular Spherical Mesh.

Here, the zero-temperature phase behavior of bosonic particles living on the nodes of a regular spherical mesh ("Platonic mesh") and interacting through an extended Bose-Hubbard Hamiltonian has been studied. Only the hard-core version of the model for two instances of Platonic mesh is considered here. Using the mean-field decoupling approximation, it is shown that the system may exist in various ground states, which can be regarded as analogs of gas, solid, supersolid, and superfluid. For one mesh, by comparing the theoretical results with the outcome of numerical diagonalization, I manage to uncover the signatures of diagonal and off-diagonal spatial orders in a finite quantum system.

Complex Self-Assembly from Simple Interaction Rules in Model Colloidal Mixtures.

Building structures with hierarchical order through the self-assembly of smaller blocks is not only a prerogative of nature, but also a strategy to design artificial materials with tailored functions. We explore in simulation the spontaneous assembly of colloidal particles into extended structures, using spheres and size-asymmetric dimers as solute particles, while treating the solvent implicitly. Besides rigid cores for all particles, we assume an effective short-range attraction between spheres and small monomers to promote, through elementary rules, dimer-mediated aggregation of spheres. Starting from a completely disordered configuration, we follow the evolution of the system at low temperature and density, as a function of the relative concentration of the two species. When spheres and large monomers are of same size, we observe the onset of elongated aggregates of spheres, either disconnected or cross-linked, and a crystalline bilayer. As spheres grow bigger, the self-assembling scenario changes, getting richer overall, with the addition of flexible membrane sheets with crystalline order and monolayer vesicles. With this wide assortment of structures, our model can serve as a viable template to achieve a better control of self-assembly in dilute suspensions of microsized particles.

Molecular dynamics determination of liquid-vapor coexistence in molten alkali halides.

We show by extensive molecular dynamics simulations that accurate predictions of liquid-vapor coexistence in molten alkali halides can be achieved in terms of a rigid ion potential description in which temperature-dependent ionic diameters are employed. The new ionic sizes result from the fitting of the experimental isothermal compressibilities, a condition whose physical implications and consequences are illustrated. The same diameters also allow us to formulate confident predictions for the compressibilities of salts in cases where the experimental data are lacking. The extension of the present approach to molten alkali-halide mixtures and to other classes of molten salts is discussed.

The barrier to ice nucleation in monatomic water.

Crystallization from a supercooled liquid initially proceeds via the formation of a small solid embryo (nucleus), which requires surmounting an activation barrier. This phenomenon is most easily studied by numerical simulation, using specialized biased-sampling techniques to overcome the limitations imposed by the rarity of nucleation events. Here, I focus on the barrier to homogeneous ice nucleation in supercooled water, as represented by the monatomic-water model, which in the bulk exhibits a complex interplay between different ice structures. I consider various protocols to identify solidlike particles on a computer, which perform well enough for the Lennard-Jones model, and compare their respective impact on the shape and height of the nucleation barrier. It turns out that the effect is stronger on the nucleus size than on the barrier height. As a by-product of the analysis, I determine the structure of the nucleation cluster, finding that the relative amount of ice phases in the cluster heavily depends on the method used for classifying solidlike particles. Moreover, the phase which is most favored during the earlier stages of crystallization may happen, depending on the nucleation coordinate adopted, to be different from the stable polymorph. Therefore, the quality of a reaction coordinate cannot be assessed simply on the basis of the barrier height obtained. I explain how this outcome is possible and why it just points out the shortcoming of collective variables appropriate to simple fluids in providing a robust method of particle classification for monatomic water.

Universal behavior of soft-core fluids near the threshold of thermodynamic stability.

We study, by using liquid-state theories and Monte Carlo simulation, the behavior of systems of classical particles interacting through a finite pair repulsion supplemented with a longer range attraction. Any such potential can be driven Ruelle-unstable by increasing the attraction at the expense of repulsion, until the thermodynamic limit is lost. By examining several potential forms, we find that all systems exhibit a qualitatively similar behavior in the fluid phase as the threshold of thermodynamic stability is approached (and possibly surpassed). The general feature underlying the approach to Ruelle instability is a pronounced widening of the liquid-vapor binodal (and spinodal) line at low temperatures, to such an extent that at the stability threshold a vanishing-density vapor would coexist with a diverging-density liquid. We attempt to rationalize the universal pathway to Ruelle instability in soft-core fluids by appealing to a heuristic argument.

Two-dimensional mixture of amphiphilic dimers and spheres: Self-assembly behaviour.

The emergence of supramolecular aggregates from simple microscopic interaction rules is a fascinating feature of complex fluids which, besides its fundamental interest, has potential applications in many areas, from biological self-assembly to smart material design. We here investigate by Monte Carlo simulation the equilibrium structure of a two-dimensional mixture of asymmetric dimers and spheres (disks). Dimers and disks are hard particles, with an additional short-range attraction between a disk and the smaller monomer of a dimer. The model parameters and thermodynamic conditions probed are typical of colloidal fluid mixtures. In spite of the minimalistic character of the interaction, we observe-upon varying the relative concentration and size of the two colloidal species-a rich inventory of mesoscale structures at low temperature, such as clusters, lamellæ (i.e., polymer-like chains), and gel-like networks. For colloidal species of similar size and near equimolar concentrations, a dilute fluid of clusters gives way to floating lamellæ upon cooling; at higher densities, the lamellæ percolate through the simulation box, giving rise to an extended network. A crystal-vapour phase-separation may occur for a mixture of dimers and much larger disks. Finally, when the fluid is brought in contact with a planar wall, further structures are obtained at the interface, from layers to branched patterns, depending on the nature of wall-particle interactions.

Self-assembly in a model colloidal mixture of dimers and spherical particles.

We investigate the structure of a dilute mixture of amphiphilic dimers and spherical particles, a model relevant to the problem of encapsulating globular "guest" molecules in a dispersion. Dimers and spheres are taken to be hard particles, with an additional attraction between spheres and the smaller monomers in a dimer. Using the Monte Carlo simulation, we document the low-temperature formation of aggregates of guests (clusters) held together by dimers, whose typical size and shape depend on the guest concentration χ. For low χ (less than 10%), most guests are isolated and coated with a layer of dimers. As χ progressively increases, clusters grow in size becoming more and more elongated and polydisperse; after reaching a shallow maximum for χ≈50%, the size of clusters again reduces upon increasing χ further. In one case only (χ=50% and moderately low temperature) the mixture relaxed to a fluid of lamellae, suggesting that in this case clusters are metastable with respect to crystal-vapor separation. On heating, clusters shrink until eventually the system becomes homogeneous on all scales. On the other hand, as the mixture is made denser and denser at low temperature, clusters get increasingly larger until a percolating network is formed.

Encapsulation of spherical nanoparticles by colloidal dimers.

We study by Monte Carlo simulation the coating process of colloidal dimers onto spherical nanoparticles. To this end we investigate a simplified mixture of hard spheres (the guest particles) and hard dimers formed by two tangent spheres of different sizes (the encapsulating agents) in an implicit-solvent representation; in our scheme, the range of effective interactions between the smaller particle in a dimer and a guest sphere depends on their relative size. By tuning the size and concentration of guests, under overall dilute conditions a rich phase behavior emerges: for small sizes and/or low concentrations, the preferred arrangement is compact aggregates (capsules) of variable sizes, where one or few guest particles are coated with dimers; for larger sizes and moderate guest concentrations, other scenarios are realized, including equilibrium separation between a guest-rich and a guest-poor phase. Our results serve as a framework for a more systematic investigation of self-assembled structures of functionalized dimers capable of encapsulating target particles, like for instance bioactive substances in a colloidal dispersion.