Microphase separation versus the vapor-liquid transition in systems of spherical particles

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.

Phase behavior of mixtures of hard colloids and soft coarse-grained macromolecules

Effective "soft" interactions between macromolecules such as polymers, amphiphilic dendrimers, and suitably designed DNA based dendritic molecules have been shown to be purely repulsive and non-diverging. We report the structure and phase behavior of a mixture of hard colloids and soft coarse-grained macromolecules. Through the use of Reference Interaction Site Model theory and molecular dynamics simulations we find that hard colloids and soft macromolecules act as depletants toward each other, generating a medium-induced effective attraction. This effective attraction leads to the formation of non-dispersed phases at high densities. At low and high fractions of hard colloids the system macrophase separates into two disparate regions of hard colloids and soft macromolecules. However, this system microphase separates into a hard-rich and soft-rich self-assembled domains at intermediate compositions. The formation of microphase separated structure in this system of isotropic, disconnected, and purely repulsive colloids is surprising and quite novel. This behavior is likely due to a softening of the interface between hard-rich and soft-rich self-assembled domains. Molecular dynamics simulations have revealed that the microphase separated state has an overall disordered bicontinuous morphology. The hard-rich domain forms an ordered FCC structure and the soft-rich domain forms a disordered cluster-fluid, making the structure simultaneously ordered and disordered.

Dioctyl sodium sulfosuccinate surfactant self-assembly dependency of solvent hydrophilicity: a modelling study

The self-assembly of dioctyl sodium sulfosuccinate (AOT) model surfactant in solvent environments of differing polarity is examined by means of dissipative particle dynamics (DPD) bead model parametrized against Hildebrand solubility parameters from atomistic molecular dynamics (MD) simulations. The model predicts that in hydrophobic solvents (e.g. dodecane) the surfactant forms small (Nagg ∼ 8) reverse micellar aggregates, while in a solvent corresponding to water lamellar assembly takes place, in good agreement with literature structural parameters. Interestingly, solvents of intermediate polarity lead to formation of large, internally structured aggregates. In these, the surfactant headgroups cluster within the aggregate, surrounded by a continuous phase formed by the hydrocarbon tails. We show that the partitioning of the headgroups between the aggregate surface layer and the inner clustered phase depends primarily on solvent polarity, and can be controlled by the solvent, but also system composition. Finally, we compare the DPD assembly response to simplified effective interaction potentials derived at dilute concentration limit for the interactions. The comparison reveals that the simplified effective potential descriptions provide good level of insight on the assembly morphologies, despite drastic, isotropic interactions simplification involved.

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.

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.

On the dynamically arrested states of equilibrium and non-equilibrium gels: a comprehensive Brownian dynamics study

In this work a systematic study over a wide number of final thermodynamic states for two gel-forming liquids was performed. Such two kind of gel formers are distinguished by their specific interparticle interaction potential. We explored several thermodynamic states determining the thermodynamic, structural and dynamic properties of both liquids after a sudden temperature change. The thermodynamic analysis allows to identify that the liquid with short range attraction and long range repulsion lacks of a stable gas-liquid phase separation liquid, in contrast with the liquid with short range attractions. Thus, although for some thermodynamic states the structural behavior, measured by the static structure factor, is similar to and characteristic of the gel phase, for the short range attractive fluid the gel phase is a consequence of a spinodal decomposition process. In contrast, gelation in the short range attraction and long range repulsion liquid is not due to a phase separation. We also analyze the similarities and differences of the dynamic behavior of both systems through the analysis of the mean square displacement, the self part of the intermediate scattering function, the diffusion coefficient and theαrelaxation time. Finally, using one of the main results of the non-equilibrium self-consistent generalized Langevin equation theory (NE-SCGLE), we determine the dynamic arrest phase diagram in the volume fraction and temperature (φvsT) plane.

Inverse design of equilibrium cluster fluids applied to a physically informed model

Inverse design strategies have proven highly useful for the discovery of interaction potentials that prompt self-assembly of a variety of interesting structures. However, often the optimized particle interactions do not have a direct relationship to experimental systems. In this work, we show that Relative Entropy minimization is able to discover physically meaningful parameter sets for a model interaction built from depletion attraction and electrostatic repulsion that yield self-assembly of size-specific clusters. We then explore the sensitivity of the optimized interaction potentials with respect to deviations in the underlying physical quantities, showing that clustering behavior is largely preserved even as the optimized parameters are perturbed.

Adsorption time scales of cluster-forming systems

A microscopic model of adsorption in cluster forming systems with competing interaction is considered. The adsorption process is described by the master equation and modelled by a kinetic Monte Carlo method. The evolution of the particle concentration and interaction energy during the adsorption of particles on a plane triangular lattice is investigated. The simulation results show a diverse behavior of the system time evolution depending on the temperature and chemical potential and finally on the formation of clusters in the system. The characteristic relaxation times of adsorption vary in several orders of magnitude depending on the thermodynamic parameters of the final equilibrium state of the adsorbate. A very fast adsorption of particles is observed for highly ordered adsorbate equilibrium states.

Equilibrium mechanisms of self-limiting assembly

Self-assembly is a ubiquitous process in synthetic and biological systems, broadly defined as the spontaneous organization of multiple subunits (e.g. macromolecules, particles) into ordered multi-unit structures. The vast majority of equilibrium assembly processes give rise to two states: one consisting of dispersed disassociated subunits, and the other, a bulk-condensed state of unlimited size. This review focuses on the more specialized class of self-limiting assembly, which describes equilibrium assembly processes resulting in finite-size structures. These systems pose a generic and basic question, how do thermodynamic processes involving non-covalent interactions between identical subunits "measure" and select the size of assembled structures? In this review, we begin with an introduction to the basic statistical mechanical framework for assembly thermodynamics, and use this to highlight the key physical ingredients that ensure equilibrium assembly will terminate at finite dimensions. Then, we introduce examples of self-limiting assembly systems, and classify them within this framework based on two broad categories: self-closing assemblies and open-boundary assemblies. These include well-known cases in biology and synthetic soft matter - micellization of amphiphiles and shell/tubule formation of tapered subunits - as well as less widely known classes of assemblies, such as short-range attractive/long-range repulsive systems and geometrically-frustrated assemblies. For each of these self-limiting mechanisms, we describe the physical mechanisms that select equilibrium assembly size, as well as potential limitations of finite-size selection. Finally, we discuss alternative mechanisms for finite-size assemblies, and draw contrasts with the size-control that these can achieve relative to self-limitation in equilibrium, single-species assemblies.

Local order and cluster formation in model fluids with competing interactions: a simulation and theoretical study

In SALR fluids, theory and simulation predict tiny morphological changes in the density profile occurring at the onset of clustering.

Microphase separation in a two-dimensional colloidal system with competing attractive critical Casimir and repulsive magnetic dipole interactions

We propose and study theoretically a colloidal system in two dimensions with attractive critical Casimir and repulsive magnetic dipole forces, wherein the strength of attraction and repulsion can be easily and independently tuned by adjusting the temperature and an external magnetic field, respectively. We expect this setup to be experimentally accessible and are confident that it can serve to deepen our understanding of the mechanisms behind microphase separation due to competing interactions. We develop a density functional theory for our model and present first results of our calculations in the form of a phase diagram for fixed temperature, but varying magnetic fields and bulk densities. For certain values of these parameters, we are able to confirm the existence of thermodynamically stable inhomogeneous density profiles in the bulk, such as parallel lamellar stripes, as well as clusters and voids on a hexagonal lattice.

Self-assembly of spiral patterns in confined systems with competing interactions

Colloidal particles in polymer solutions and functionalized nanoparticles often exhibit short-range attraction coupled with long-range repulsion (SALR) leading to the spontaneous formation of symmetric patterns. Chiral nanostructures formed by thin films of SALR particles have not been reported yet. In this study, we observe striking topological transitions from a symmetric pattern of concentric rings to a chiral structure of a spiral shape, when the system is in hexagonal confinement. We find that the spiral formation can be induced either by breaking the system symmetry with a wedge, or by melting of the rings. In the former case, the chirality of the spiral is determined by the orientation of the wedge and thus can be controlled. In the latter, the spiral arises due to thermally induced defects and is absent in the average particle distribution, which forms highly regular hexagonal patterns in the central part of the system. These hexagonal patterns can be explained by interference of planar density waves. Thermodynamic considerations indicate that equilibrium spirals can appear spontaneously in any stripe-forming system confined in a hexagon with a small wedge, provided that certain conditions are satisfied by a set of phenomenological parameters.

Some general features of mesophase formation in hard-core plus tail potentials

The formation of mesophases in fluids with hard-core plus tail interactions is investigated and compared with the occurrence of cluster crystals in ultra-soft repulsive potentials by using a simple variational expression for the Helmholtz free energy. The purpose of this study is mostly qualitative, i.e., to explain the origin of the different behavior of these systems, and the reason why, in the hard-core case, interactions which are apparently quite different display a common pattern for the phase diagram, featuring spheres, cylinders, lamellae, inverted cylinders, and inverted spheres as the density is increased. In the limit of zero temperature, our approach also yields some simple predictions for the densities at which the transitions between different mesophases are expected to take place, as well as for the size of their clusters at the transitions. We find that these results compare favorably with those obtained in a former study of a model fluid with competing attractive and repulsive interactions by density-functional theory with numerical minimization.

Complex structures generated by competing interactions in harmonically confined colloidal suspensions

We investigate the structural properties of colloidal particle systems interacting via an isotropic pair potential and confined by a three-dimensional harmonic potential. The interaction potential has a repulsive-attractive-repulsive profile that varies with the interparticle distance (also known as a 'mermaid' potential). We performed Langevin dynamics simulations to find the equilibrium configurations of the system. We show that particles can self-assemble in complex structural patterns, such as compact disks, fringed disks, rods, spherical clusters with superficial entrances among others. Also, for particular values of the parameters of the interaction potential, we could identify that some configurations were formed by quasi two-dimensional (2D) structures which are stable for 2D systems.

Combined density functional and Brazovskii theories for systems with spontaneous inhomogeneities

The low-T part of the phase diagram in self-assembling systems is correctly predicted by known versions of density functional theory (DFT). The high-T part obtained in DFT, however, does not agree with simulations even on the qualitative level. In this work, a new version of DFT for systems with spontaneous inhomogeneities on a mesoscopic length scale is developed. The contribution to the grand thermodynamic potential associated with mesoscopic fluctuations is explicitly taken into account. The expression for this contribution is obtained by methods known from the Brazovskii field theory. Apart from developing the approximate expression for the grand thermodynamic potential that contains the fluctuation contribution and is ready for numerical minimization, we develop a simplified version of the theory valid for weakly ordered phases, i.e. for the high-T part of the phase diagram. The simplified theory is verified by comparison with the results of simulations for a particular version of the short-range attraction long-range repulsion (SALR) interaction potential. Except for the fact that in our theory the ordered phases are stable at lower T than in simulations, a good agreement for the high-T part of the phase diagram is obtained for the range of density that was considered in simulations. In addition, the equation of state and compressibility isotherms are presented. Finally, the physical interpretation of the fluctuation-contribution to the grand potential is discussed in detail.

Structural crossover in a model fluid exhibiting two length scales: Repercussions for quasicrystal formation

We investigate the liquid state structure of the two-dimensional model introduced by Barkan et al. [Phys. Rev. Lett. 113, 098304 (2014)10.1103/PhysRevLett.113.098304], which exhibits quasicrystalline and other unusual solid phases, focusing on the radial distribution function g(r) and its asymptotic decay r→∞. For this particular model system, we find that as the density is increased there is a structural crossover from damped oscillatory asymptotic decay with one wavelength to damped oscillatory asymptotic decay with another distinct wavelength. The ratio of these wavelengths is ≈1.932. Following the locus in the phase diagram of this structural crossover leads directly to the region where quasicrystals are found. We argue that identifying and following such a crossover line in the phase diagram towards higher densities where the solid phase(s) occur is a good strategy for finding quasicrystals in a wide variety of systems. We also show how the pole analysis of the asymptotic decay of equilibrium fluid correlations is intimately connected with the nonequilibrium growth or decay of small-amplitude density fluctuations in a bulk fluid.

Sliding states of a soft-colloid cluster crystal: Cluster versus single-particle hopping

We study a two-dimensional model for interacting colloidal particles which displays spontaneous clustering. Within this model we investigate the competition between the pinning to a periodic corrugation potential and a sideways constant pulling force which would promote a sliding state. For a few sample particle densities and amplitudes of the periodic corrugation potential we investigate the depinning from the statically pinned to the dynamically sliding regime. This sliding state exhibits the competition between a dynamics where entire clusters are pulled from a minimum to the next and a dynamics where single colloids or smaller groups leave a cluster and move across the corrugation energy barrier to join the next cluster downstream in the force direction. Both kinds of sliding states can occur either coherently across the entire sample or asynchronously: the two regimes result in different average mobilities. Finite temperature tends to destroy separate sliding regimes, generating a smoother dependence of the mobility on the driving force.

Cluster Crystals Stabilized by Hydrophobic and Electrostatic Interactions

Cluster crystals are crystalline materials in which each site is occupied by multiple identical particles, atoms, colloids, or polymers. There are two classes of systems that make cluster crystals. One is composed of particles that interact via potentials that are bound at the origin and thus are able to penetrate each other. The other consists of non-interpenetrating particles whose interaction potential diverges at the origin. The goal of this work is to find which systems of the second class can make cluster crystals that are stable at room temperature. First, the general properties of the required potentials are established using an analytical model and Monte Carlo simulations. Next, we ask how such potentials can be constructed by combining hydrophobic attraction and electrostatic repulsion. A colloid model with a hard-sphere core and a repulsive wall is introduced to mimic the hydrophobic interaction. Charge is added to create long-range repulsion. A search in the parameter space of the colloid size, counterion type, and charge configuration uncovers several models for which effective colloid-colloid interaction, determined in explicit solvent as a potential of mean force, has the necessary shape. For the effective potential, cluster crystals are confirmed as low free-energy configurations in replica-exchange molecular dynamics simulations, which also generate the respective transition temperature. The model that exhibits a transition above room temperature is further studied in explicit solvent. Simulations on a 10 ns time scale show that crystalline conformations are stable below the target temperature but disintegrate rapidly above it, supporting the idea that hydrophobic and electrostatic interactions are sufficient to induce an assembly of cluster crystals. Finally, we discuss which physical systems are good candidates for experimental observations of cluster crystals.

Pattern formation and self-assembly driven by competing interactions

Colloidal fluids interacting via effective potentials which are attractive at the short range and repulsive at the long range have long been raising considerable attention because such an instance provides a simple mechanism leading to pattern formation even for isotropic interactions. If the competition between attraction and repulsion is strong enough, the gas-liquid phase transition is suppressed, and replaced by the formation of mesophases, i.e., inhomogeneous phases displaying periodic density modulations whose length, although being larger than the particle size, cannot nevertheless be considered macroscopic. We describe a fully numerical implementation of density-functional theory in three dimensions, tailored to periodic phases. The results for the equilibrium phase diagram of the model are compared with those already obtained in previous investigations for the present system as well as for other systems which form mesophases. The phase diagram which we find shows a strong similarity with that of block copolymer melts, in which self-assembly also results from frustration of a macroscopic phase separation. As the inhomogeneous region is swept by increasing the density from the low-density side, one encounters clusters, bars, lamellae, inverted bars, and inverted clusters. Moreover, a bicontinuous gyroid phase consisting of two intertwined percolating networks is predicted in a narrow domain between the bar and lamellar phases.

Phase behavior and bulk structural properties of a microphase former with anisotropic competing interactions: A density functional theory study

Using classical density functional theory, we investigate systems exhibiting interactions where a short-range anisotropic attractive force competes with a long-range spherically symmetric repulsive force. The former is modelled within Wertheim's first-order perturbation theory for patchy particles, and the repulsive part is assumed to be a Yukawa potential which is taken into account via a mean-field approximation. From previous studies of systems with spherically symmetric competing interactions, it is well known that such systems can exhibit stable bulk cluster phases (microphase separation) provided that the repulsion is sufficiently weak compared to the attraction. For the present model system, we find rich phase diagrams including both reentrant clustering and liquid-gas binodals. In particular, the model predicts inhomogeneous bulk phases at extremely low packing fractions, which cannot be observed in systems with isotropic competing interactions.

Interactions and design rules for assembly of porous colloidal mesophases

Porous mesophases, where well-defined particle-depleted 'void' spaces are present within a particle-rich background fluid, can be self-assembled from colloidal particles interacting via isotropic pair interactions with competing attractions and repulsions. While such structures could be of wide interest for technological applications (e.g., filtration, catalysis, absorption, etc.), relatively few studies have investigated the interactions that lead to these morphologies and how they compare to those that produce other micro-phase-separated structures, such as clusters. In this work, we use inverse methods of statistical mechanics to design model isotropic pair potentials that form porous mesophases. We characterize the resulting porous structures, correlating features of the pair potential with the targeted pore size and the particle packing fraction. The former is primarily encoded by the amplitude and range of the repulsive barrier of the designed pair potential and the latter by the attractive well depth. We observe a trade-off with respect to the packing fraction of the targeted morphology: greater values support more spherical and monodisperse pores that themselves organize into periodic structures, while lower values yield more mobile pores that do not assemble into ordered structures but remain stable over a larger range of packing fraction. We conclude by commenting on the limitations of targeting a specific pore diameter within the present inverse design approach as well as by describing future directions to overcome these limitations.

Tiny changes in local order identify the cluster formation threshold in model fluids with competing interactions

In two-Yukawa fluids, the increment of attraction causes a tiny reversal of trend in the local density, within specific ranges of interparticle distances. This event precisely signals the onset of a clustered state.

Effects of confinement on pattern formation in two dimensional systems with competing interactions

Template-assisted pattern formation in monolayers of particles with competing short-range attraction and long-range repulsion interactions (SALR) is studied by Monte Carlo simulations in a simple generic model [N. G. Almarza et al., J. Chem. Phys., 2014, 140, 164708]. We focus on densities corresponding to formation of parallel stripes of particles and on monolayers laterally confined between straight parallel walls. We analyze both the morphology of the developed structures and the thermodynamic functions for broad ranges of temperature T and the separation L2 between the walls. At low temperature stripes parallel to the boundaries appear, with some corrugation when the distance between the walls does not match the bulk periodicity of the striped structure. The stripes integrity, however, is rarely broken for any L2. This structural order is lost at T = TK(L2) depending on L2 according to a Kelvin-like equation. Above the Kelvin temperature TK(L2) many topological defects such as breaking or branching of the stripes appear, but a certain anisotropy in the orientation of the stripes persists. Finally, at high temperature and away from the walls, the system behaves as an isotropic fluid of elongated clusters of various lengths and with various numbers of branches. For L2 optimal for the stripe pattern the heat capacity as a function of temperature takes the maximum at T = TK(L2).

Alternating strings and clusters in suspensions of charged colloids

We report the formation of alternating strings and clusters in a binary suspension of repulsive charged colloids with double layers larger than the particle size. Within a binary cell model we include many-body and charge-regulation effects under the assumption of a constant surface potential, and consider their repercussions on the two-particle interaction potential. We find that the formation of induced dipoles close to a charge-reversed state may explain the formation of these structures. Finally, we will touch upon the formation of dumbbells and small clusters in a one-component system, where the effective electrostatic interaction is always repulsive.

Gyroid phase of fluids with spherically symmetric competing interactions

We study the phase diagram of a fluid with spherically symmetric competing pair interactions that consist of a short-ranged attraction and a longer-ranged repulsion in addition to a hard core. To this end we perform free minimizations of three-dimensional triple periodic structures within the framework of classical density functional theory. We compare our results to those from Landau theory. Our main finding is that the double gyroid phase can exist as a thermodynamically stable phase.

Density functional theory for systems with mesoscopic inhomogeneities

We study the effects of fluctuations on the mesoscopic length scale on systems with mesoscopic inhomogeneities. Equations for the correlation function and for the average volume fraction are derived in the self-consistent Gaussian approximation. The equations are further simplified by postulating the expression for the structure factor consistent with scattering experiments for self-assembling systems. Predictions of the approximate theory are verified by a comparison with the exact results obtained earlier for the one-dimensional lattice model with first-neighbor attraction and third-neighbor repulsion. We find qualitative agreement for the correlation function, the equation of state and the dependence of the chemical potential μ on the volume fraction ζ. Our results confirm also that strong inhomogeneities in the disordered phase are found only in the case of strong repulsion. The inhomogeneities are reflected in an oscillatory decay of the correlation function with a very large correlation length, three inflection points in the [Formula: see text] curve and a compressibility that for increasing ζ takes very large, very small and again very large values.

Assembly of nothing: equilibrium fluids with designed structured porosity

Controlled micro- to meso-scale porosity is a common materials design goal with possible applications ranging from molecular gas adsorption to particle size selective permeability or solubility. Here, we use inverse methods of statistical mechanics to design an isotropic pair interaction that, in the absence of an external field, assembles particles into an inhomogeneous fluid matrix surrounding pores of prescribed size ordered in a lattice morphology. The pore size can be tuned via modification of temperature or particle concentration. Moreover, modulating density reveals a rich series of microphase-separated morphologies including pore- or particle-based lattices, pore- or particle-based columns, and bicontinuous or lamellar structures. Sensitivity of pore assembly to the form of the designed interaction potential is explored.

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.

Effective interactions in active Brownian suspensions

Active colloids exhibit persistent motion, which can lead to motility-induced phase separation (MIPS). However, there currently exists no microscopic theory to account for this phenomenon. We report a first-principles theory, free of fit parameters, for active spherical colloids, which shows explicitly how an effective many-body interaction potential is generated by activity and how this can rationalize MIPS. For a passively repulsive system the theory predicts phase separation and pair correlations in quantitative agreement with simulation. For an attractive system the theory shows that phase separation becomes suppressed by moderate activity, consistent with recent experiments and simulations, and suggests a mechanism for reentrant cluster formation at high activity.

Structural orderings of anisotropically confined colloids interacting via a quasi-square-well potential

We implement Brownian dynamics to investigate the static properties of colloidal particles confined anisotropically and interacting via a potential which can be tailored in a repulsive-attractive-respulsive fashion as the interparticle distance increases. A diverse number of structural phases are self-assembled, which were classified according to two aspects, that is, their macroscopic and microscopic patterns. Concerning the microscopic phases we found the quasicrystalline, triangular, square, and mixed orderings, where this latter is a combination of square and triangular cells in a 3×2 proportion, i.e., the so-called (3(3),4(2)) Archimedian lattice. On the macroscopic level the system could self-organize in a compact or perforated single cluster surrounded or not by fringes. All the structural phases are summarized in detailed phases diagrams, which clearly show that the different phases are extended as the confinement potential becomes more anisotropic.

Structural ordering of trapped colloids with competing interactions

The structure of colloids with competing interactions which are confined in a harmonic external trap potential is analyzed numerically by energy minimization in two spatial dimensions. A wealth of different cluster structures is found to be stable including clusters with a fringed outer rim (reminiscent to an ornamental border), clusters perforated with voids, as well as clusters with a crystalline core and a disordered rim. All cluster structures occur in a two-dimensional parameter space. The structural ordering can therefore be efficiently tuned by changing few parameters only providing access to a controlled fabrication of colloidal clusters.

Effective isotropic potential for dipolar hard spheres

A new effective isotropic potential is proposed for the dipolar hard-sphere fluid, on the basis of recent results by others for its angle-averaged radial distribution function. The new effective potential is shown to exhibit oscillations even for moderately high densities and moderately strong dipole moments, which are absent from earlier effective isotropic potentials. The validity and significance of this result are briefly discussed.