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.

Correct scaling of the correlation length from a theory for concentrated electrolytes.

Self-consistent theory for concentrated electrolytes is developed. Oscillatory decay of the charge-charge correlation function with the decay length that shows perfect agreement with the experimentally discovered and so far unexplained scaling is obtained. For the density-density correlations, monotonic asymptotic decay with the decay length comparable with the decay length of the charge correlations is found. We show that the correlation lengths in concentrated electrolytes depend crucially on the local variance of the charge density.

Triangular lattice models for pattern formation by core-shell particles with different shell thicknesses.

Triangular lattice models for pattern formation by hard-core soft-shell particles at interfaces are introduced and studied in order to determine the effect of the shell thickness and structure. In model I, we consider particles with hard-cores covered by shells of cross-linked polymeric chains. In model II, such inner shell is covered by a much softer outer shell. In both models, the hard cores can occupy sites of the triangular lattice, and nearest-neighbor repulsion following from overlapping shells is assumed. The capillary force is represented by the second or the fifth neighbor attraction in model I or II, respectively. Ground states with fixed chemical potential µ or with fixed fraction of occupied sites c are thoroughly studied. For T > 0, the µ(c) isotherms, compressibility and specific heat are calculated by Monte Carlo simulations. In model II, 6 ordered periodic patterns occur in addition to 4 phases found in model I. These additional phases, however, are stable only at the phase coexistence lines at the (µ, T) diagram, which otherwise looks like the diagram of model I. In the canonical ensemble, these 6 phases and interfaces between them appear in model II for large intervals of c and the number of possible patterns is much larger than in model I. We calculated line tensions for different interfaces, and found that the favorable orientation of the interface corresponds to its smoothest shape in both models.

Adsorption anomalies in a two-dimensional model of cluster-forming systems.

Adsorption on a boundary line confining a monolayer of particles self-assembling into clusters is studied by Monte Carlo simulations. We focus on a system of particles interacting via competing interaction potential in which effectively short-range attraction is followed by long-range repulsion. For the chemical potential values below the order-disorder phase transition the adsorption isotherms were shown to undergo nonstandard behavior, i.e., the adsorption exhibits a maximum on structural transition between structureless and disordered cluster fluid. In particular, we have found that the adsorption decreases for increasing chemical potential when (i) clusters dominate over monomers in the bulk, (ii) the density profile in the direction perpendicular to the confining line exhibits an oscillatory decay, and (iii) the correlation function in the layer near the adsorption wall exhibits an oscillatory decay in the direction parallel to this wall. Our report indicates striking differences between simple and complex fluid adsorption processes.

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.

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.

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.

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.

Orientational ordering of lamellar structures on closed surfaces.

Self-assembly of particles with short-range attraction and long-range repulsion interactions on a flat and on a spherical surface is compared. Molecular dynamics simulations are performed for the two systems having the same area and the density optimal for formation of stripes of particles. Structural characteristics, e.g., a cluster size distribution, a number of defects, and an orientational order parameter (OP), as well as the specific heat, are obtained for a range of temperatures. In both cases, the cluster size distribution becomes bimodal and elongated clusters appear at the temperature corresponding to the maximum of the specific heat. When the temperature decreases, orientational ordering of the stripes takes place and the number of particles per cluster or stripe increases in both cases. However, only on the flat surface, the specific heat has another maximum at the temperature corresponding to a rapid change of the OP. On the sphere, the crossover between the isotropic and anisotropic structures occur in a much broader temperature interval; the orientational order is weaker and occurs at significantly lower temperature. At low temperature, the stripes on the sphere form spirals and the defects resemble defects in the nematic phase of rods adsorbed at a sphere.

Exactly solvable model for self-assembly of hard core-soft shell particles at interfaces.

A generic model for self-assembly of a monolayer of hybrid core-shell particles at an interface is developed. We assume that for distances larger than the size of the incompressible core a soft repulsion appears, and the repulsion is followed by an attraction at larger separations. The model is solved exactly in a one-dimensional lattice version. One, two or three periodic structures and variety of shapes of the pressure-density isotherms may occur in different versions of the model. For strong interactions the isotherm consists of nearly vertical segments at densities optimal for the periodic structures that are connected by segments with a small slope. The range of order depends very strongly on the strength of attraction and on the density. Our results agree with experimental observations.

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).

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.

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.

Effects of rigid or adaptive confinement on colloidal self-assembly. Fixed vs. fluctuating number of confined particles.

The effects of confinement on colloidal self-assembly in the case of fixed number of confined particles are studied in the one dimensional lattice model solved exactly in the grand canonical ensemble (GCE) in PÈ©kalski et al. [J. Chem. Phys. 142, 014903 (2015)]. The model considers a pair interaction defined by a short-range attraction plus a longer-range repulsion. We consider thermodynamic states corresponding to self-assembly into clusters. Both fixed and adaptive boundaries are studied. For fixed boundaries, there are particular states in which, for equal average densities, the number of clusters in the GCE is larger than in the canonical ensemble. The dependence of pressure on density has a different form when the system size changes with fixed number of particles and when the number of particles changes with fixed size of the system. In the former case, the pressure has a nonmonotonic dependence on the system size. The anomalous increase of pressure for expanding system is accompanied by formation of a larger number of smaller clusters. In the case of elastic confining surfaces, we observe a bistability, i.e., two significantly different system sizes occur with almost the same probability. The mechanism of the bistability in the closed system is different to that of the case of permeable walls, where the two equilibrium system sizes correspond to a different number of particles.

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.

Periodic ordering of clusters and stripes in a two-dimensional lattice model. II. Results of Monte Carlo simulation.

The triangular lattice model with nearest-neighbor attraction and third-neighbor repulsion, introduced by PÈ©kalski, Ciach, and Almarza [J. Chem. Phys. 140, 114701 (2014)] is studied by Monte Carlo simulation. Introduction of appropriate order parameters allowed us to construct a phase diagram, where different phases with patterns made of clusters, bubbles or stripes are thermodynamically stable. We observe, in particular, two distinct lamellar phases-the less ordered one with global orientational order and the more ordered one with both orientational and translational order. Our results concern spontaneous pattern formation on solid surfaces, fluid interfaces or membranes that is driven by competing interactions between adsorbing particles or molecules.

Periodic ordering of clusters and stripes in a two-dimensional lattice model. I. Ground state, mean-field phase diagram and structure of the disordered phases.

The short-range attraction and long-range repulsion between nanoparticles or macromolecules can lead to spontaneous pattern formation on solid surfaces, fluid interfaces, or membranes. In order to study the self-assembly in such systems we consider a triangular lattice model with nearest-neighbor attraction and third-neighbor repulsion. At the ground state of the model (T = 0) the lattice is empty for small values of the chemical potential µ, and fully occupied for large µ. For intermediate values of µ periodically distributed clusters, bubbles, or stripes appear if the repulsion is sufficiently strong. At the phase coexistences between the vacuum and the ordered cluster phases and between the cluster and the lamellar (stripe) phases the entropy per site does not vanish. As a consequence of this ground state degeneracy, disordered fluid phases consisting of clusters or stripes are stable, and the surface tension vanishes. For T > 0 we construct the phase diagram in the mean-field approximation and calculate the correlation function in the self-consistent Brazovskii-type field theory.

Periodic ordering of clusters in a one-dimensional lattice model.

A generic lattice model for systems containing particles interacting with short-range attraction long-range repulsion (SALR) potential that can be solved exactly in one dimension is introduced. We assume attraction J1 between the first neighbors and repulsion J2 between the third neighbors. The ground state of the model shows existence of two homogeneous phases (gas and liquid) for J2/J1 <1/3. In addition to the homogeneous phases, the third phase with periodically distributed clusters appears for J2/J1 > 1/3. Phase diagrams obtained in the self-consistent mean-field approximation for a range of values of J2/J1 show very rich behavior, including reentrant melting, and coexistence of two periodic phases (one with strong and the other one with weak order) terminated at a critical point. We present exact solutions for the equation of state as well as for the correlation function for characteristic values of J2/J1. Based on the exact results, for J2/J1 > 1/3 we predict pseudo-phase transitions to the ordered cluster phase indicated by a rapid change of density for a very narrow range of pressure, and by a very large correlation length for thermodynamic states where the periodic phase is stable in mean field. For 1/9 < J2/J1 < 1/3 the correlation function decays monotonically below certain temperature, whereas above this temperature exponentially damped oscillatory behavior is obtained. Thus, even though macroscopic phase separation is energetically favored and appears for weak repulsion at T = 0, local spatial inhomogeneities appear for finite T. Monte Carlo simulations in canonical ensemble show that specific heat has a maximum for low density ρ that we associate with formation of living clusters, and if the repulsion is strong, another maximum for ρ = 1/2.

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.

Separation of components in lipid membranes induced by shape transformation.

Vesicles composed of a two component membrane with each component characterized by different spontaneous curvature are investigated by minimization of the free energy consisting of Helfrich elastic energy and entropy of mixing. The results show that mixing and demixing of membrane components can be induced by elongating a vesicle or changing its volume, if one of the components forms a complex with macromolecules on the outer monolayer. The influence of elastic coefficients on the separation of components is also examined.