Periodic ordering of clusters in a one-dimensional lattice model

From disorder to order: A dynamic approach to mesophase formation in soft sphere model

This study explores the dynamics of self-assembly and mesophase formation through molecular dynamics simulations of hexagonal and lamellar systems using a simplified coarse-grained model. We focus on characterizing the order-disorder transitions driven by temperature variations and emphasize the often overlooked disordered regime, which serves as a precursor to periodic mesoscale ordering. Our findings not only underscore the morphological richness of the disordered regime, comparable to that of its periodic counterparts, but also reveal the presence of clustering regimes within isotropic phases, thus corroborating prior experimental and theoretical observations. By employing the dynamic correlation coefficient, this work introduces a novel approach to understanding the fundamental mechanisms of mesophase formation, providing new insights into the complex dynamics of self-assembly.

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.

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.

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.

Clustering and assembly dynamics of a one-dimensional microphase former

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

Thermodynamic signatures and cluster properties of self-assembly in systems with competing interactions

Colloidal particles, amphiphiles and functionalized nanoparticles are examples of systems that frequently exhibit short-range attraction coupled with long-range repulsion. We vary the ratio of attraction and repulsion in a simple isotropic model with competing interactions, using molecular simulations, and observe significant differences in the properties of the self-assembled clusters that form. We report conditions that lead to the self-assembly of clusters of a preferred size, accompanied by a change in the slope of the pressure with respect to density, similar to micelles formed by amphiphilic molecules. We also report conditions where repulsion dominates, clusters of a preferred size form and the pressure vs. density slope is unaffected by self-assembly. We investigate cluster structure by calculating the size distributions, free colloid density, cluster shape and density profiles. The system dynamics are characterized by cluster life-times. We do not find qualitative differences in structure or dynamics of the clusters, regardless the pressure behavior. Therefore, thermodynamic and structural quantities are required to classify the different clustering characteristics that are observable in systems with competing interactions. Our results have implications in terms of development of design principles for stable cluster self-assembly.

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.

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.

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

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.

Probing the existence of phase transitions in one-dimensional fluids of penetrable particles

Phase transitions in one-dimensional classical fluids are usually ruled out by using van Hove's theorem. A way to circumvent the conclusions of the theorem is to consider an interparticle potential that is everywhere bounded. Such is the case of, e.g., the generalized exponential model of index 4 (GEM-4 potential), which in three dimensions gives a reasonable description of the effective repulsion between flexible dendrimers in a solution. An extensive Monte Carlo simulation of the one-dimensional GEM-4 model [S. Prestipino, Phys. Rev. E 90, 042306 (2014)] has recently provided evidence of an infinite sequence of low-temperature cluster phases, however, also suggesting that upon pushing the simulation forward what seemed a true transition may eventually prove to be only a sharp crossover. We hereby investigate this problem theoretically by use of three different and increasingly sophisticated approaches (i.e., a mean-field theory, the transfer matrix of a lattice model of clusters, and the exact treatment of a system of point clusters in the continuum) to conclude that the alleged transitions of the one-dimensional GEM-4 system are likely just crossovers.

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.

Nematic phase in the J(1)-J(2) square-lattice Ising model in an external field

The J(1)-J(2) Ising model in the square lattice in the presence of an external field is studied by two approaches: the cluster variation method (CVM) and Monte Carlo simulations. The use of the CVM in the square approximation leads to the presence of a new equilibrium phase, not previously reported for this model: an Ising-nematic phase, which shows orientational order but not positional order, between the known stripes and disordered phases. Suitable order parameters are defined, and the phase diagram of the model is obtained. Monte Carlo simulations are in qualitative agreement with the CVM results, giving support to the presence of the new Ising-nematic phase. Phase diagrams in the temperature-external field plane are obtained for selected values of the parameter κ=J(2)/|J(1)| which measures the relative strength of the competing interactions. From the CVM in the square approximation we obtain a line of second order transitions between the disordered and nematic phases, while the nematic-stripes phase transitions are found to be of first order. The Monte Carlo results suggest a line of second order nematic-disordered phase transitions in agreement with the CVM results. Regarding the stripes-nematic transitions, the present Monte Carlo results are not precise enough to reach definite conclusions about the nature of the transitions.

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.