Lattice Model Results for Pattern Formation in a Mixture with Competing Interactions.

A monolayer consisting of two types of particles, with energetically favored alternating stripes of the two components, is studied by Monte Carlo simulations and within a mesoscopic theory. We consider a triangular lattice model and assume short-range attraction and long-range repulsion between particles of the same kind, as well as short-range repulsion and long-range attraction for the cross-interaction. The structural evolution of the model upon increasing temperature is studied for equal chemical potentials of the two species. We determine the structure factor, the chemical potential-density isotherms, the specific heat, and the compressibility, and show how these thermodynamic functions are associated with the spontaneous formation of stripes with varying degrees of order.

Pattern Formation in Two-Component Monolayers of Particles with Competing Interactions.

Competing interactions between charged inclusions in membranes of living organisms or charged nanoparticles in near-critical mixtures can lead to self-assembly into various patterns. Motivated by these systems, we developed a simple triangular lattice model for binary mixtures of oppositely charged particles with additional short-range attraction or repulsion between like or different particles, respectively. We determined the ground state for the system in contact with a reservoir of the particles for the whole chemical potentials plane, and the structure of self-assembled conglomerates for fixed numbers of particles. Stability of the low-temperature ordered patterns was verified by Monte Carlo simulations. In addition, we performed molecular dynamics simulations for a continuous model with interactions having similar features, but a larger range and lower strength than in the lattice model. Interactions with and without symmetry between different components were assumed. We investigated both the conglomerate formed in the center of a thin slit with repulsive walls, and the structure of a monolayer adsorbed at an attractive substrate. Both models give the same patterns for large chemical potentials or densities. For low densities, more patterns occur in the lattice model. Different phases coexist with dilute gas on the lattice and in the continuum, leading to different patterns in self-assembled conglomerates ('rafts').

Short-time dynamics in active systems: the Vicsek model.

We study the short-time dynamics (STD) of the Vicsek model (VM) with vector noise. The study of STD has proved to be very useful in the determination of the critical point, critical exponents and spinodal points in equilibrium phase transitions. Here we aim is to test its applicability in active systems. We find that, despite the essential non-equilibrium characteristics of the VM (absence of detailed balance, activity), the STD presents qualitatively the same phenomenology as in equilibrium systems. From the STD one can distinguish whether the transition is continuous or discontinuous (which we have checked also computing the Binder cumulant). When the transition is continuous, one can determine the critical point and the critical exponents.

Self-propelled Vicsek particles at low speed and low density.

We study through numerical simulation the Vicsek model for very low speeds and densities. We consider scalar noise in two and three dimensions and vector noise in three dimensions. We focus on the behavior of the critical noise with density and speed, trying to clarify seemingly contradictory earlier results. We find that, for scalar noise, the critical noise is a power law in both density and speed, but although we confirm the density exponent in two dimensions, we find a speed exponent different from earlier reports (we consider lower speeds than previous studies). On the other hand, for the vector noise case we find that the dependence of the critical noise cannot be separated as a product of power laws in speed and density. Finally, we study the dependence of the relaxation time with speed. At the critical point we find a power law, with the same exponent in two and three dimensions.

Unexpected crossover dynamics of single polymer in a corrugated tube.

We present molecular dynamics study of a generic (coarse-grained) model for single-polymer diffusion confined in a corrugated cylinder. For a narrow tube, i.e., diameter of the cylinder δ < 2.3, the axial diffusion coefficient D(∣∣) scales as D(∣∣) ∝ N(-3∕2), with chain length N, up to N ≈ 100 and then crosses over to Rouse scaling for the larger N values. The N(-3∕2) scaling is due to the large fluctuation of the polymer chain along its fully stretched equilibrium conformation. The stronger scaling, namely N(-3∕2), is not observed for an atomistically smooth tube and/or for a cylinder with larger diameter.

From capillary condensation to interface localization transitions in colloid-polymer mixtures confined in thin-film geometry.

Monte Carlo simulations of the Asakura-Oosawa model for colloid-polymer mixtures confined between two parallel repulsive structureless walls are presented and analyzed in the light of current theories on capillary condensation and interface localization transitions. Choosing a polymer-to-colloid size ratio of q=0.8 and studying ultrathin films in the range of D=3 to D=10 colloid diameters thickness, grand canonical Monte Carlo methods are used; phase transitions are analyzed via finite size scaling, as in previous work on bulk systems and under confinement between identical types of walls. Unlike the latter work, inequivalent walls are used here: While the left wall has a hard-core repulsion for both polymers and colloids, at the right-hand wall an additional square-well repulsion of variable strength acting only on the colloids is present. We study how the phase separation into colloid-rich and colloid-poor phases occurring already in the bulk is modified by such a confinement. When the asymmetry of the wall-colloid interaction increases, the character of the transition smoothly changes from capillary condensation type to interface localization type. For very thin films (i.e., for D=3 ) and a suitable choice of the wall-colloid interactions, evidence is found that the critical behavior falls in the universality class of the two-dimensional Ising model. Otherwise, we observe crossover scaling between different universality classes (namely, the crossover from the three-dimensional to the two-dimensional Ising model universality class). The colloid and polymer density profiles across the film in the various phases are discussed, as well as the correlation of interfacial fluctuations in the direction parallel to the confining walls. The broadening of the interface between the coexisting colloid-rich and polymer-rich phases (located parallel to the confining walls) is understood in terms of capillary wave fluctuations. The experimental observability of all these phenomena is briefly discussed.

Confinement effects on phase behavior of soft matter systems.

When systems that can undergo phase separation between two coexisting phases in the bulk are confined in thin film geometry between parallel walls, the phase behavior can be profoundly modified. These phenomena shall be described and exemplified by computer simulations of the Asakura-Oosawa model for colloid-polymer mixtures, but applications to other soft matter systems (e.g. confined polymer blends) will also be mentioned. Typically a wall will prefer one of the phases, and hence the composition of the system in the direction perpendicular to the walls will not be homogeneous. If both walls are of the same kind, this effect leads to a distortion of the phase diagram of the system in thin film geometry, in comparison with the bulk, analogous to the phenomenon of "capillary condensation" of simple fluids in thin capillaries. In the case of "competing walls", where both walls prefer different phases of the two phases coexisting in the bulk, a state with an interface parallel to the walls gets stabilized. The transition from the disordered phase to this "soft mode phase" is rounded by the finite thickness of the film and is not a sharp phase transition. However, a sharp transition can occur where this interface gets localized at (one of) the walls. The relation of this interface localization transition to wetting phenomena is discussed. Finally, an outlook to related phenomena is given, such as the effects of confinement in cylindrical pores on the phase behavior, and more complicated ordering phenomena (lamellar mesophases of block copolymers or nematic phases of liquid crystals under confinement).