Universal sequence of ordered structures obtained from mesoscopic description of self-assembly

Spontaneous pattern formation in monolayers of binary mixtures with competing interactions

A model for a monolayer of two types of particles spontaneously forming ordered patterns is studied using a mesoscopic theory and MC simulations. We assume hard-cores of the same size a for both components. For r > a, like particles attract and repel each other at short and large distances, respectively, with the same potential u(r) for both species, and the cross-interaction is -u(r). The model is inspired by oppositely charged particles or macromolecules with preferential solubility in different components of a solvent that is close to a miscibility critical point, in particular by inclusions in biological membranes. We obtain the phase diagram in the chemical potentials and temperature variables as well as in the concentration, density and temperature variables, using the mean-field one-shell approximation. We find that the presence of the second component significantly extends the temperature range of stability of the ordered phases. We obtain three stable phases with periodic concentration: the lamellar L phase with alternating stripes of the two components for similar chemical potentials, and a hexagonal arrangement of the clusters of the minority component in the liquid of the majority component. The latter two phases, however, are stable only at relatively high temperatures. At lower temperatures, the L phase coexists with a disordered one-component fluid or with very dilute gas with mixed components. At still lower temperatures, the one-component phase coexisting with the L phase can be disordered or ordered, depending on the chemical potentials. The theoretical results are confirmed by MC simulations for selected thermodynamic states.

Phase diagram of SALR fluids on spherical surfaces

We investigate the phase diagram of a fluid of hard-core disks confined to the surface of a sphere and whose interaction potential contains a short-range attraction followed by a long-range repulsive tail (SALR). Based on previous work in the bulk we derive a stability criterion for the homogeneous phase of the fluid, and locate a region of instability linked to the presence of a negative minimum in the spherical harmonics expansion of the interaction potential. The inhomogeneous phases contained within this region are characterized using a mean-field density functional theory. We find several inhomogeneous patterns that can be separated into three broad classes: cluster crystals, stripes, and bubble crystals, each containing topological defects. Interestingly, while the periodicity of inhomogeneous phases at large densities is mainly determined by the position of the negative minimum of the potential expansion, the finite size of the system induces a richer behavior at smaller densities.

Structural and Dynamical Behaviour of Colloids with Competing Interactions Confined in Slit Pores

Systems with short-range attractive and long-range repulsive interactions can form periodic modulated phases at low temperatures, such as cluster-crystal, hexagonal, lamellar and bicontinuous gyroid phases. These periodic microphases should be stable regardless of the physical origin of the interactions. However, they have not yet been experimentally observed in colloidal systems, where, in principle, the interactions can be tuned by modifying the colloidal solution. Our goal is to investigate whether the formation of some of these periodic microphases can be promoted by confinement in narrow slit pores. By performing simulations of a simple model with competing interactions, we find that both the cluster-crystal and lamellar phases can be stable up to higher temperatures than in the bulk system, whereas the hexagonal phase is destabilised at temperatures somewhat lower than in bulk. Besides, we observed that the internal ordering of the lamellar phase can be modified by changing the pore width. Interestingly, for sufficiently wide pores to host three lamellae, there is a range of temperatures for which the two lamellae close to the walls are internally ordered, whereas the one at the centre of the pore remains internally disordered. We also find that particle diffusion under confinement exhibits a complex dependence with the pore width and with the density, obtaining larger and smaller values of the diffusion coefficient than in the corresponding bulk system.

Adsorption in Mixtures with Competing Interactions

A binary mixture of oppositely charged particles with additional short-range attraction between like particles and short-range repulsion between different ones in the neighborhood of a substrate preferentially adsorbing the first component is studied by molecular dynamics simulations. The studied thermodynamic states correspond to an approach to the gas–crystal coexistence. Dependence of the near-surface structure, adsorption and selective adsorption on the strength of the wall–particle interactions and the gas density is determined. We find that alternating layers or bilayers of particles of the two components are formed, but the number of the adsorbed layers, their orientation and the ordered patterns formed inside these layers could be quite different for different substrates and gas density. Different structures are associated with different numbers of adsorbed layers, and for strong attraction the thickness of the adsorbed film can be as large as seven particle diameters. In all cases, similar amount of particles of the two components is adsorbed, because of the long-range attraction between different particles.

Formation and internal ordering of periodic microphases in colloidal models with competing interactions

Theory and simulations predict that colloidal particles with short-range attractive and long-range repulsive interactions form periodic microphases if there is a proper balance between the attractive and repulsive contributions. However, the experimental identification of such structures has remained elusive to date. Using molecular dynamics simulations, we investigate the phase behaviour of a model system that stabilizes a cluster-crystal, a cylindrical and a lamellar phase at low temperatures. Besides the transition from the fluid to the periodic microphases, we also observe the internal freezing of the clusters at a lower temperature. Finally, our study indicates that, for the chosen model parameters, the three periodic microphases are kinetically accessible from the fluid phase.

Self-assembly in mixtures with competing interactions

A binary mixture of particles interacting with spherically-symmetrical potentials leading to microsegregation is studied by theory and molecular dynamics (MD) simulations. We consider spherical particles with equal diameters and volume fractions. Motivated by the mixture of oppositely charged particles with different adsorption preferences immersed in a near-critical binary solvent, we assume short-range attraction long-range repulsion for the interaction between like particles, and short-range repulsion long-range attraction for the interaction between different ones. In order to predict structural and thermodynamic properties of such complex mixtures, we develop a theory combining the density functional and field-theoretical methods. We show that concentration fluctuations in mesoscopic regions lead to a qualitative change of the phase diagram compared to mean-field predictions. Both theory and MD simulations show coexistence of a low-density disordered phase with a high-density phase with alternating layers rich in the first and second components. In these layers, crystalline structure is present in the solid, and absent in the liquid crystals. The density and the degree of order of the ordered phase decrease with increasing temperature, up to a temperature where the theory predicts a narrow two-phase region with increasing density of both phases for increasing temperature. MD simulations show that monocrystals of the solid and liquid crystals have a prolate shape with the axis parallel to the direction of concentration oscillations, and the deviation from the spherical shape increases with increasing periodic order.

Statics and dynamics of a finite two-dimensional colloidal system with competing attractive critical Casimir and repulsive magnetic dipole interactions

We continue our theoretical study of a recently proposed two-dimensional colloidal system with attractive critical Casimir and repulsive magnetic dipole forces that can be tuned easily and independently from each other via the temperature and the strength of an external magnetic field, respectively [K. Marolt, M. Zimmermann, and R. Roth, Phys. Rev. E 100, 052602 (2019)2470-004510.1103/PhysRevE.100.052602]. Using this freedom, it is possible to construct a competing interaction potential that causes microphase separation featuring spatially inhomogeneous cluster, stripe, and bubble phases in the bulk, i.e., in an infinite system without an external potential. In the present work, we demonstrate by means of density functional theory that microphase separation can also occur in finite geometries. In a square cell with a side length of 20 or 30 colloid diameters, we observe the emergence of highly structured cluster and ring phases at intermediate bulk densities in addition to almost uniform fluid phases for lower and higher bulk densities. We then employ dynamic density functional theory to determine how the system reacts when the temperature and the magnetic field are altered over time, and we show how to induce a transition from the liquid to the cluster/ring phase and also from the cluster directly to the ring phase. We find that often a slowly varying and nontrivial path in parameter space is required to reach a stable state, whereas abrupt changes are prone to lead to metastable configurations.

Competing active and passive interactions drive amoebalike crystallites and ordered bands in active colloids

Swimmers and self-propelled particles are physical models for the collective behavior and motility of a wide variety of living systems, such as bacteria colonies, bird flocks, and fish schools. Such artificial active materials are amenable to physical models which reveal the microscopic mechanisms underlying the collective behavior. Here we study colloids in a dc electric field. Our quasi-two-dimensional system of electrically driven particles exhibits a rich and exotic phase behavior exhibiting passive crystallites, motile crystallites, an active gas, and banding. Amongst these are two mesophases, reminiscent of systems with competing interactions. At low field strengths activity suppresses demixing, leading to motile crystallites. Meanwhile, at high field strengths, activity drives partial demixing to traveling bands. We parametrize a particulate simulation model which reproduces the experimentally observed phases.

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.

Mesoscopic theory for systems with competing interactions near a confining wall

Mesoscopic theory for self-assembling systems near a planar confining surface is developed. Euler-Lagrange equations and the boundary conditions (BCs) for the local volume fraction and the correlation function are derived from the density functional theory expression for the grand thermodynamic potential. Various levels of approximation can be considered for the obtained equations. The lowest-order nontrivial approximation [generic model (GM)] resembles the Landau-Brazovskii-type theory for a semi-infinite system. Unlike in the original phenomenological theory, however, all coefficients in our equations and BCs are expressed in terms of the interaction potential and the thermodynamic state. Analytical solutions of the linearized equations in the GM are presented and discussed on a general level and for a particular example of the double-Yukawa potential. We show exponentially damped oscillations of the volume fraction and the correlation function in the direction perpendicular to the confining surface. The correlations show oscillatory decay in directions parallel to this surface too, with the decay length increasing significantly when the system boundary is approached. The framework of our theory allows for a systematic improvement of the accuracy of the results.

Effect of aggregation on adsorption phenomena

Adsorption at an attractive surface in a system with particles self-assembling into small clusters is studied by molecular dynamics simulation. We assume Lennard-Jones plus repulsive Yukawa tail interactions and focus on small densities. The relative increase in the temperature at the critical cluster concentration near the attractive surface (CCCS) shows a power-law dependence on the strength of the wall-particle attraction. At temperatures below the CCCS, the adsorbed layer consists of undeformed clusters if the wall-particle attraction is not too strong. Above the CCCS or for strong attraction leading to flattening of the adsorbed aggregates, we obtain a monolayer that for strong or very strong attraction consists of flattened clusters or stripes, respectively. The accumulated repulsion from the particles adsorbed at the wall leads to a repulsive barrier that slows down the adsorption process, and the accession time grows rapidly with the strength of the wall-particle attraction. Beyond the adsorbed layer of particles, a depletion region of a thickness comparable with the range of the repulsive tail of interactions occurs, and the density in this region decreases with increasing strength of the wall-particle attraction. At larger separations, the exponentially damped oscillations of density agree with theoretical predictions for self-assembling systems. Structural and thermal properties of the bulk are also determined. In particular, a new structural crossover associated with the maximum of the specific heat and a double-peaked histogram of the cluster size distribution are observed.

Following the nucleation pathway from disordered liquid to gyroid mesophase

Mesophases have order intermediate between liquids and crystals and arise in systems with frustration, such as surfactants, block copolymers, and Janus nanoparticles. The gyroid mesophase contains two interpenetrated, nonintersecting chiral networks that give it properties useful for photonics. It is challenging to nucleate a gyroid from the liquid. Elucidating the reaction coordinate for gyroid nucleation could assist in designing additives that facilitate the formation of the mesophase. However, the complexity of the gyroid structure and the extreme weakness of the first-order liquid to gyroid transition make this a challenging quest. Here, we investigate the pathway and transition states for the nucleation of a gyroid from the liquid in molecular simulations with a mesogenic binary mixture. We find that the gyroid nuclei at the transition states have a large degree of positional disorder and are not compact, consistent with the low surface free energy of the liquid-gyroid interface. A combination of bond-order parameters for the minor component is best to describe the passage from liquid to gyroid, among those we consider. The committor analyses, however, show that this best coordinate is not perfect and suggests that accounting for the relative ordering of the two interpenetrated networks in infant nuclei, as well as for signatures of ordering in the major component of the mesophase, would improve the accuracy of the reaction coordinate for gyroid formation and its use to evaluate nucleation barriers. To our knowledge, this study is the first to investigate the reaction coordinate and critical nuclei for the formation of any mesophase from an amorphous phase.

Modeling micro-heterogeneity in mixtures: The role of many body correlations

A two-component interaction model is introduced herein, which allows us to describe macroscopic miscibility with various modes of tunable micro-segregation, ranging from phase separation to micro-segregation, and is in excellent agreement with structural quantities obtained from simulations and the liquid state hypernetted-chain like integral equation theory. The model is based on the conjecture that the many-body correlation bridge function term in the closure relation can be divided into one part representing the segregation effects, which are modeled herein, and the usual part representing random many body fluctuations. Furthermore, the model allows us to fully neglect these second contributions, thus increasing the agreement between the simulations and the theory. The analysis of the retained part of the many body correlations gives important clues about how to model the many body bridge functions for more realistic systems exhibiting micro-segregation, such as aqueous mixtures.

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.

Nonequilibrium phase transitions of sheared colloidal microphases: Results from dynamical density functional theory

By means of classical density functional theory and its dynamical extension, we consider a colloidal fluid with spherically symmetric competing interactions, which are well known to exhibit a rich bulk phase behavior. This includes complex three-dimensional periodically ordered cluster phases such as lamellae, two-dimensional hexagonally packed cylinders, gyroid structures, or spherical micelles. While the bulk phase behavior has been studied extensively in earlier work, in this paper we focus on such structures confined between planar repulsive walls under shear flow. For sufficiently high shear rates, we observe that microphase separation can become fully suppressed. For lower shear rates, however, we find that, e.g., the gyroid structure undergoes a kinetic phase transition to a hexagonally packed cylindrical phase, which is found experimentally and theoretically in amphiphilic block copolymer systems. As such, besides the known similarities between the latter and colloidal systems regarding the equilibrium phase behavior, our work reveals further intriguing nonequilibrium relations between copolymer melts and colloidal fluids with competing interactions.

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.

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.

Microphase Separation in Oil-Water Mixtures Containing Hydrophilic and Hydrophobic Ions

We develop a lattice-based Monte Carlo simulation method for charged mixtures capable of treating dielectric heterogeneities. Using this method, we study oil-water mixtures containing an antagonistic salt, with hydrophilic cations and hydrophobic anions. Our simulations reveal several phases with a spatially modulated solvent composition, in which the ions partition between water-rich and water-poor regions according to their affinity. In addition to the recently observed lamellar phase, we find tubular and droplet phases, reminiscent of those found in block copolymers and surfactant systems. Interestingly, these structures stem from ion-mediated interactions, which allows for tuning of the phase behavior via the concentrations, the ionic properties, and the temperature.

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.

Self-Assembly of Mesophases from Nanoparticles

A growing number of crystalline and quasi-crystalline structures have been formed by coating nanoparticles with ligands, polymers, and DNA. The design of nanoparticles that assemble into mesophases, such as those formed by block copolymers, would combine the order, mobility, and stimuli responsive properties of mesophases with the electronic, magnetic, and optical properties of nanoparticles. Here we use molecular simulations to demonstrate that binary mixtures of unbound particles with simple short-ranged pair interactions produce the same mesophases as block copolymers and surfactants, including lamellar, hexagonal, gyroid, body-centered cubic, face-centered cubic, perforated lamellar, and semicrystalline phases. The key to forming the mesophases is the frustrated attraction between particles of different types, achieved through control over interparticle size and over strength and softness of the interaction. Experimental design of nanoparticles with effective interactions described by the potentials of this work would provide a distinct, robust route to produce ordered tunable liquid crystalline mesophases from nanoparticles.

Purely entropic self-assembly of the bicontinuous Ia3d gyroid phase in equilibrium hard-pear systems

We investigate a model of hard pear-shaped particles which forms the bicontinuous Ia[Formula: see text]d structure by entropic self-assembly, extending the previous observations of Barmes et al. (2003 Phys. Rev. E68, 021708. (doi:10.1103/PhysRevE.68.021708)) and Ellison et al. (2006 Phys. Rev. Lett.97, 237801. (doi:10.1103/PhysRevLett.97.237801)). We specifically provide the complete phase diagram of this system, with global density and particle shape as the two variable parameters, incorporating the gyroid phase as well as disordered isotropic, smectic and nematic phases. The phase diagram is obtained by two methods, one being a compression-decompression study and the other being a continuous change of the particle shape parameter at constant density. Additionally, we probe the mechanism by which interdigitating sheets of pears in these systems create surfaces with negative Gauss curvature, which is needed to form the gyroid minimal surface. This is achieved by the use of Voronoi tessellation, whereby both the shape and volume of Voronoi cells can be assessed in regard to the local Gauss curvature of the gyroid minimal surface. Through this, we show that the mechanisms prevalent in this entropy-driven system differ from those found in systems which form gyroid structures in nature (lipid bilayers) and from synthesized materials (di-block copolymers) and where the formation of the gyroid is enthalpically driven. We further argue that the gyroid phase formed in these systems is a realization of a modulated splay-bend phase in which the conventional nematic has been predicted to be destabilized at the mesoscale due to molecular-scale coupling of polar and orientational degrees of freedom.

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.

Self-consistent theory for systems with mesoscopic fluctuations

We have developed a theory for inhomogeneous systems that allows for the incorporation of the effects of mesoscopic fluctuations. A hierarchy of equations relating the correlation and direct correlation functions for the local excess [Formula: see text] of the volume fraction of particles ζ has been obtained, and an approximation leading to a closed set of equations for the two-point functions has been introduced for the disordered inhomogeneous phase. We have numerically solved the self-consistent equations for one-dimensional (1D) and three-dimensional (3D) models with short-range attraction and long-range repulsion. Predictions for all of the qualitative properties of the 1D model agree with the exact results, but only semi-quantitative agreement is obtained in the simplest version of the theory. The effects of fluctuations in the two 3D models considered are significantly different, despite the very similar properties of these models in the mean-field approximation. In both cases we obtain the sequence of large-small-large compressibility for increasing ζ. The very small compressibility is accompanied by the oscillatory decay of correlations with correlation lengths that are orders of magnitude larger than the size of particles. In one of the two models considered, the small compressibility becomes very small and the large compressibility becomes very large with decreasing temperature, and eventually van der Waals loops appear. Further studies are necessary in order to determine the nature of the strongly inhomogeneous phase present for intermediate volume fractions in 3D.

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.

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.

Bistability in a self-assembling system confined by elastic walls: exact results in a one-dimensional lattice model

The impact of confinement on self-assembly of particles interacting with short-range attraction and long-range repulsion potential is studied for thermodynamic states corresponding to local ordering of clusters or layers in the bulk. Exact and asymptotic expressions for the local density and for the effective potential between the confining surfaces are obtained for a one-dimensional lattice model introduced by J. Pȩkalski et al. [J. Chem. Phys. 138, 144903 (2013)]. The simple asymptotic formulas are shown to be in good quantitative agreement with exact results for slits containing at least 5 layers. We observe that the incommensurability of the system size and the average distance between the clusters or layers in the bulk leads to structural deformations that are different for different values of the chemical potential μ. The change of the type of defects is reflected in the dependence of density on μ that has a shape characteristic for phase transitions. Our results may help to avoid misinterpretation of the change of the type of defects as a phase transition in simulations of inhomogeneous systems. Finally, we show that a system confined by soft elastic walls may exhibit bistability such that two system sizes that differ approximately by the average distance between the clusters or layers are almost equally probable. This may happen when the equilibrium separation between the soft boundaries of an empty slit corresponds to the largest stress in the confined self-assembling system.

How ions in solution can change the sign of the critical Casimir potential

We show that hydrophilic ions present in a confined, near-critical aqueous mixture can lead to an attraction between like charge surfaces with opposing preferential adsorption of the two species of the mixture, even though the corresponding Casimir potential in uncharged systems is repulsive. This prediction agrees with a recent experiment [Nellen et al., Soft Matter, 2011, 7, 5360]. We also show that oppositely charged hydrophobic surfaces can repel each other, although the Casimir potential between uncharged surfaces with like preferential adsorption (selectivity) is attractive. This behavior is expected when the electrostatic screening length is larger than the correlation length, and one of the confining surfaces is strongly selective and weakly charged, whereas the other confining surface is weakly selective and strongly charged. The Casimir potential can change sign because the hydrophilic ions near the weakly hydrophobic surface can overcompensate the effect of hydrophobicity, and this surface can act as a hydrophilic one. We also predict a more attractive interaction between charged, hydrophilic surfaces and a more repulsive interaction between charged, hydrophobic surfaces than given by the sum of the Casimir and Debye-Hückel potentials. Our theory is derived systematically from a microscopic approach, and combines the Landau-type and Debye-Hückel theories with an additional contribution of an entropic origin.

Microphase separations of the fluids with spherically symmetric competing interactions

A density functional perturbation theory has been developed for studying the phase behaviors of a competing system in the spherical pores. The pore size as well as the intensity of competing interactions exerts a strong influence on the vapor-liquid, vapor-cluster, and cluster-liquid transitions of a competing system. The microdomain spacing (D) of the cluster is commensurate with the periodicity of modulation in the particle density distributions of a competing system in a spherical pore with the pore radius (R). For the cluster phase, we find that the multi-vaporlike void is formed depending on the periodicity of modulation by finite-size artifacts. For R < D, the competing system only shows the vapor-liquid transition at a high amplitude. For R > D, the vapor-cluster and cluster-liquid transitions are found at a high amplitude, whereas at a low amplitude, the cluster-liquid transition only occurs. The competing system exhibits two tricritical points, which are joined to one another by the line of second-order transitions at the low and high densities. A comparison with the result of a slit pore shows that (i) the tricritical points in a spherical pore, which has the highest symmetry, occur at a low amplitude compared with that of a slit pore because of the geometrical properties of the pores, and that (ii) the slit pore relatively shows the wide vapor-cluster and cluster-liquid coexistence regions compared with that of a spherical pore: the geometrical symmetry of a pore results in a weaker tendency for phase separation.

Spatial inhomogeneities in ionic liquids, charged proteins, and charge stabilized colloids from collective variables theory

Effects of size and charge asymmetry between oppositely charged ions or particles on spatial inhomogeneities are studied for a large range of charge and size ratios. We perform a stability analysis of the primitive model of ionic systems with respect to periodic ordering using the collective variables-based theory. We extend previous studies [Ciach et al., Phys. Rev. E 75, 051505 (2007)] in several ways. First, we employ a nonlocal approximation for the reference hard-sphere fluid which leads to the Percus-Yevick pair direct correlation functions for the uniform case. Second, we use the Weeks-Chandler-Anderson regularization scheme for the Coulomb potential inside the hard core. We determine the relevant order parameter connected with the periodic ordering and analyze the character of the dominant fluctuations along the λ lines. We show that the above-mentioned modifications produce large quantitative and partly qualitative changes in the phase diagrams obtained previously. We discuss possible scenarios of the periodic ordering for the whole range of size and charge ratios of the two ionic species, covering electrolytes, ionic liquids, charged globular proteins or nanoparticles in aqueous solutions, and charge-stabilized colloids.

Structure and phase behavior of two-Yukawa fluids with competing interactions in planar slit pores

A density functional perturbation theory, which is based both on the modified fundamental-measure theory and on the first-order mean-spherical approximation for long-range attractive and repulsive interactions, has been developed for studying the structure and phase behaviors of a competing system restricted to slit pores. The hysteresis loop for the adsorption and desorption curves indicates that the system exhibits vapor-cluster and cluster-liquid transitions which depend on the pair potential parameters and the slit width (H). The periodic spacing (D) of the cluster is commensurate with the periodicity of modulation in the particle density distribution and more closely related to the vapor-cluster and cluster-liquid phase transitions of the system. For the cluster phase, we find the transition from a single liquidlike slab to a multi-liquidlike slab with increasing the slit width. The multi-liquidlike slab is formed depending on the periodicity of modulation by finite-size artifacts. The cluster-related phase transitions, such as the vapor-cluster or cluster-liquid transitions occur for H>D, while for H<D the system only exhibits the vapor-liquid transition. At a low amplitude, the vapor-liquid transition disappears and the cluster-liquid transition only occurs for H<D. The coexistence curves for the confined phase diagram are contained within the corresponding bulk liquid-vapor coexistence curve. For a wide slit pore (H>D), the system exhibits two tricritical points, joined to one another by the line of second-order transition. The results support the conclusion that the confinement effect plays an important role in determining the equilibrium phase transition.