Structural ordering of trapped colloids with competing interactions

Structural evolution of particle configurations: Zero-temperature phases under increasing confinement

In this study, we investigate the phase behavior and structural organization of colloidal particles in a two-dimensional (2D) system under isotropic harmonic confinement using overdamped Langevin dynamics simulations. We employ a modified mermaid potential, which introduces an additional short-distance term resulting in a null-force region, distinct from the conventional mermaid potential. This modification facilitates a richer exploration of self-assembled structures, revealing a variety of phases influenced by the interplay between confinement strength V0 and the interaction potential. Our analysis spans a wide range of parameters, resulting in a detailed phase diagram that captures transitions from dispersed clusters to well-ordered patterns, including square, triangular, rhomboidal, and mixed configurations, as the confinement strength increases. The findings underscore the intricate balance of forces governing the self-assembly of colloidal systems and offer valuable insights for future experimental realizations.

Computational self-assembly of a six-fold chiral quasicrystal

Quasicrystals are unique materials characterized by long-range order without periodicity. They are observed in systems such as metallic alloys, soft matter, and particle simulations. Unlike periodic crystals, which are invariant under real-space symmetry operations, quasicrystals possess symmetry that requires description by a space group in reciprocal space. In this study, we report the self-assembly of a six-fold chiral quasicrystal using molecular dynamics simulations of a two-dimensional particle system. The particles interact via the Lennard-Jones-Gauss pair potential and are subjected to a periodic substrate potential. We confirm the presence of chiral symmetry through diffraction patterns and order parameters, revealing unique local motifs in both real and reciprocal space. The quasicrystal's properties, including the tiling structure and symmetry and the extent of diffuse scattering, are strongly influenced by substrate potential depth and temperature. Our results provide insights into the mechanisms of chiral quasicrystal formation and the role and potential of external fields in tailoring quasicrystal structures.

Quasi-two-dimensional dispersions of Brownian particles with competitive interactions: phase behavior and structural properties

Competing short-range attractive (SA) and long range repulsive (LR) particle interactions can be used to describe three-dimensional charge-stabilized colloid or protein dispersions at low added salt concentrations, as well as membrane proteins with interaction contributions mediated by lipid molecules. Using Langevin dynamics (LD) simulations, we determine the generalized phase diagram, cluster shapes and size distributions of a generic quasi-two-dimensional (Q2D) dispersion of spherical SALR particles confined to in-plane motion inside a bulk fluid. The SA and LR interaction parts are modelled by a generalized Lennard-Jones potential and a screened Coulomb potential, respectively. The microstructures of the detected equilibrium and non-equilibrium Q2D phases are distinctly different from those observed in three-dimensional (3D) SALR systems, by exhibiting different levels of hexagonal ordering. We discuss a thermodynamic perturbation theory prediction for the metastable binodal line of a reference system of particles with SA interactions only, which in the explored Q2D-SALR phase diagram region separates cluster from non-clustered phases. The transition from the high-temperature (small SA) dispersed fluid (DF) phase to the lower-temperature equilibrium cluster (EC) fluid phase is characterised by a low-wavenumber peak height of the static structure factor (corresponding to a thermal correlation length of about twice the particle diameter) featuring a distinctly smaller value (≈1.4) than in 3D SALR systems. With decreasing temperature (increasing SA), the cluster morphology changes from disk-like shapes in the equilibrium cluster phase, to double-stranded anisotropic hexagonal cluster segments formed in a cluster-percolated (CP) gel-like phase. This transition can be quantified by a hexagonal order parameter distribution function. The mean cluster size and coordination number of particles in the CP phase are insensitive to changes in the attraction strength.

Peak effect and dynamics of stripe- and pattern-forming systems on a periodic one-dimensional substrate

We examine the ordering, pinning, and dynamics of two-dimensional pattern-forming systems interacting with a periodic one-dimensional substrate. In the absence of the substrate, particles with competing long-range repulsion and short-range attraction form anisotropic crystal, stripe, and bubble states. When the system is tuned across the stripe transition in the presence of a substrate, we find that there is a peak effect in the critical depinning force when the stripes align and become commensurate with the substrate. Under an applied drive, the anisotropic crystal and stripe states can exhibit soliton depinning and plastic flow. When the stripes depin plastically, they dynamically reorder into a moving stripe state that is perpendicular to the substrate trough direction. We also find that when the substrate spacing is smaller than the widths of the bubbles or stripes, the system forms pinned stripe states that are perpendicular to the substrate trough direction. The system exhibits multiple reentrant pinning effects as a function of increasing attraction, with the anisotropic crystal and large bubble states experiencing weak pinning but the stripe and smaller bubble states showing stronger pinning. We map out the different dynamic phases as a function of filling, the strength of the attractive interaction term, the substrate strength, and the drive, and demonstrate that the different phases produce identifiable features in the transport curves and particle orderings.

Local orderings in the melting process of a square crystal and in the resulting liquid

Using Brownian dynamics simulations we investigate the melting processes of a square crystalline lattice of colloidal particles interacting via an isotropic potential, which comprises both a hard-core repulsion and an additional softened square-well potential. For temperatures slightly lower than the transition one, we found a proliferation of small liquid clusters surrounded by the square lattice. These clusters are not static, quite the opposite, they have an intense dynamics and are continuously formed and destroyed over time. However, no unbound topological defects are observed. At the transition temperature, one of these liquid clusters starts to grow, until the entire system becomes in the liquid phase, then, characterizing a first-order phase transition. The tetratic intermediate phase, as given by the KTHNY theory, was not observed. Moreover, the liquid phase exhibits a considerable number of crystalline clusters having square and triangular orderings, which remain present even when increasing temperature by an order of magnitude. As the temperature increases, structural changes within the liquid phase are analyzed by evaluating the number and sizes of the square and triangular clusters. A transition of the dominant clusters is observed.

Structural transitions in two-dimensional modulated systems under triangular confinement

We study numerically the structural transitions of two-dimensional systems of classic particles with competing interactions under a triangular confinement with two different types of soft-wall potentials. We observe a variety of novel confinement-induced equilibrium configurations as a function of particle density and confinement steepness for each considered confinement potential. The specific role played by the confining potentials on the ordering of the particle clusters is revealed. These findings allow us to control the self-organization of modulated systems through using external confinements.

Ordering and Dynamics of Interacting Colloidal Particles under Soft Confinement

Confinement can induce substantial changes in the physical properties of macromolecules in suspension. Soft confinement is a particular class of restriction where the boundaries that constraint the particles in a region of the space are not well-defined. This scenario leads to a broader structural and dynamical behavior than observed in systems enclosed between rigid walls. In this contribution, we study the ordering and diffusive properties of a two-dimensional colloidal model system subjected to a one-dimensional parabolic trap. Increasing the trap strength makes it possible to go through weak to strong confinement, allowing a dimensional transition from two- to one-dimension. The non-monotonic response of the static and dynamical properties to the gradual dimensionality change affects the system phase behavior. We find that the particle dynamics are connected to the structural transitions induced by the parabolic trap. In particular, at low and intermediate confinement regimes, complex structural and dynamical scenarios arise, where the softness of the external potential induces melting and freezing, resulting in faster and slower particle diffusion, respectively. Besides, at strong confinements, colloids move basically along one direction, and the whole system behaves structurally and dynamically similar to a one-dimensional colloidal system.

Nonequilibrium pattern formation in circularly confined two-dimensional systems with competing interactions

We numerically investigate the nonequilibrium behaviors of classic particles with competing interactions confined in a two-dimensional logarithmic trap. We reveal a quench-induced surprising dynamics exhibiting rich dynamic patterns depending upon confinement strength and trap size, which is attributed to the time-dependent competition between interparticle repulsions and attractions under a circular confinement. Moreover, in the collectively diffusive motions of the particles, we find that the emergence of dynamic structure transformation coincides with a diffusive mode transition from superdiffusion to subdiffusion. These findings are likely useful in understanding the pattern selection and evolution in various chemical and biological systems in addition to modulated systems, and add a new route to tailoring the morphology of pattern-forming systems.

Dynamics of Magnus-dominated particle clusters, collisions, pinning, and ratchets

Motivated by the recent work in skyrmions and active chiral matter systems, we examine pairs and small clusters of repulsively interacting point particles in the limit where the dynamics is dominated by the Magnus force. We find that particles with the same Magnus force can form stable pairs, triples, and higher ordered clusters or exhibit chaotic motion. For mixtures of particles with opposite Magnus force, particle pairs can combine to form translating dipoles. Under an applied drive, particles with the same Magnus force translate; however, particles with different or opposite Magnus force exhibit a drive-dependent decoupling transition. When the particles interact with a repulsive obstacle, they can form localized orbits with depinning or unwinding transitions under an applied drive. We examine the interaction of these particles with clusters or lines of obstacles and find that the particles can become trapped in orbits that encircle multiple obstacles. Under an ac drive, we observe a series of ratchet effects, including ratchet reversals, for particles interacting with a line of obstacles due to the formation of commensurate orbits. Finally, in assemblies of particles with mixed Magnus forces of the same sign, we find that the particles with the largest Magnus force become localized in the center of the cluster, while for mixtures with opposite Magnus forces, the motion is dominated by transient local pairs or clusters, where the translating pairs can be regarded as a form of active matter.

Colloidal systems with a short-range attraction and long-range repulsion: Phase diagrams, structures, and dynamics

Colloidal systems with both a short-range attraction and long-range repulsion (SALR) have rich phases compared with the traditional hard sphere systems or sticky hard sphere systems. The competition between the short-range attraction and long-range repulsion results in the frustrated phase separation, which leads to the formation of intermediate range order (IRO) structures and introduces new phases to both equilibrium and nonequilibrium phase diagrams, such as clustered fluid, cluster percolated fluid, Wigner glass, and cluster glass. One hallmark feature of many SALR systems is the appearance of the IRO peak in the interparticle structure factor, which is associated with different types of IRO structures. The relationship between the IRO peak and the clustered fluid state has been careful investigated. Not surprisingly, the morphology of clusters in solutions can be affected and controlled by the SALR potential. And the effect of the SALR potential on the dynamic properties is also reviewed here. Even though much progress has been made in understanding SALR systems, many future works are still needed to have quantitative comparisons between experiments and simulations/theories and understand the differences from different experimental systems. Owing to the large parameter space available for SALR systems, many exciting features of SALR systems are not fully explored yet. Because proteins in low-salinity solutions have SALR interactions, the understanding of SALR systems can greatly help understand protein behavior in concentrated solutions or crowded conditions.

2D melting of confined colloids with a mixture of square and triangular order

We implement Brownian dynamics to investigate the melting processes of colloidal particles confined isotropically and interacting via a potential which can be tailored in a repulsive-attractive-respulsive fashion as the interparticle distance increases. The stable configurations of such a system is composed of a large diversity of structures, which includes quasicrystalline, triangular, square, and mixed orderings, as well as the presence of fringes and holes, which are located, respectively, at the border and interior of the clusters. Our simulations demonstrate that during the melting process particles are able to swing between different micro phases. This intermediary stage, present in a finite range of temperature, precedes the melting in all cases investigated and is different from the hexatic phase of the KTNHY framework. We also test the fringes stability and find it to be higher than the one found in compact clusters. Finally, we show that, at the high temperature regime, the system loses its angular ordering while still preserves its radial interparticle confinement, which, ultimately, causes the proliferation of small subclusters.

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.

Static and dynamic properties of two-dimensional Coulomb clusters

We study the temperature dependence of static and dynamic responses of Coulomb interacting particles in two-dimensional confinements across the crossover from solid- to liquid-like behaviors. While static correlations that investigate the translational and bond orientational order in the confinements show the footprints of hexatic-like phase at low temperatures, dynamics of the particles slow down considerably in this phase, reminiscent of a supercooled liquid. Using density correlations, we probe long-lived heterogeneities arising from the interplay of the irregularity in the confinement and long-range Coulomb interactions. The relaxation at multiple time scales show stretched-exponential decay of spatial correlations in irregular traps. Temperature dependence of characteristic time scales, depicting the structural relaxation of the system, show striking similarities with those observed for the glassy systems, indicating that some of the key signatures of supercooled liquids emerge in confinements with lower spatial symmetries.

Structural transitions and hysteresis in clump- and stripe-forming systems under dynamic compression

Using numerical simulations, we study the dynamical evolution of particles interacting via competing long-range repulsion and short-range attraction in two dimensions. The particles are compressed using a time-dependent quasi-one dimensional trough potential that controls the local density, causing the system to undergo a series of structural phase transitions from a low density clump lattice to stripes, voids, and a high density uniform state. The compression proceeds via slow elastic motion that is interrupted with avalanche-like bursts of activity as the system collapses to progressively higher densities via plastic rearrangements. The plastic events vary in magnitude from small rearrangements of particles, including the formation of quadrupole-like defects, to large-scale vorticity and structural phase transitions. In the dense uniform phase, the system compresses through row reduction transitions mediated by a disorder-order process. We characterize the rearrangement events by measuring changes in the potential energy, the fraction of sixfold coordinated particles, the local density, and the velocity distribution. At high confinements, we find power law scaling of the velocity distribution during row reduction transitions. We observe hysteresis under a reversal of the compression when relatively few plastic rearrangements occur. The decompressing system exhibits distinct phase morphologies, and the phase transitions occur at lower compression forces as the system expands compared to when it is compressed.

Shell formation in short like-charged polyelectrolytes in a harmonic trap

Inspired by recent experiments and simulations on pattern formation in biomolecules by optical tweezers, a theoretical description based on the reference interaction site model (RISM) is developed to calculate the equilibrium density profiles of small polyelectrolytes in an external potential. The formalism is applied to the specific case of a finite number of Gaussian and rodlike polyelectrolytes trapped in a harmonic potential. The density profiles of the polyelectrolytes are studied over a range of lengths and numbers of polyelectrolytes in the trap, and the strength of the trap potential. For smaller polymers we recover the results for point charges. In the mean field limit the longer polymers, unlike point charges, form a shell at the boundary layer. When the interpolymer correlations are included, the density profiles of the polymers show sharp shells even at weaker trap strengths. The implications of these results are discussed.

Structural ordering of self-assembled clusters with competing interactions: transition from faceted to spherical clusters

The self-assembly of nanoparticles into clusters and the effect of the different parameters of the competing interaction potential on it are investigated. For a small number of particles, the structural organization of the clusters is almost unaffected by the attractive part of the potential, and for an intermediate number of particles the configuration strongly depends on the strength of it. The cluster size is controlled by the range of the interaction potential, and the structural arrangement is guided by the strength of the potential: i.e., the self-assembled cluster transforms from a faceted configuration at low strength to a spherical shell-like structure at high strength. Nonmonotonic behavior of the cluster size is found by increasing the interaction range. An approximate analytical expression is obtained that predicts the smallest cluster for a specific set of potential parameters. A Mendeleev-like table is constructed for different values of the strength and range of the attractive part of the potential in order to understand the structural ordering of the ground-state configuration of the self-assembled clusters.

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

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

Stripe systems with competing interactions on quasi-one dimensional periodic substrates

We numerically examine the two-dimensional ordering of a stripe forming system of particles with competing long-range repulsion and short-range attraction in the presence of a quasi-one-dimensional corrugated substrate. As a function of increasing substrate strength or period we show that a remarkable variety of distinct orderings can be realized, including modulated stripes, prolate clump phases, two dimensional ordered kink structures, crystalline void phases, and smectic phases. Additionally in some cases the stripes align perpendicular to the substrate troughs. Our results suggest that a new route to self assembly for systems with competing interactions can be achieved through the addition of a simple periodic modulated substrate.

Self-organization of highly symmetric nanoassemblies: a matter of competition

The properties and applications of metallic nanoparticles are inseparably connected not only to their detailed morphology and composition but also to their structural configuration and mutual interactions. As a result, the assemblies often have superior properties as compared to individual nanoparticles. Although it has been reported that nanoparticles can form highly symmetric clusters, if the configuration can be predicted as a function of the synthesis parameters, more targeted and accurate synthesis will be possible. We present here a theoretical model that accurately predicts the structure and configuration of self-assembled gold nanoclusters. The validity of the model is verified using quantitative experimental data extracted from electron tomography 3D reconstructions of different assemblies. The present theoretical model is generic and can in principle be used for different types of nanoparticles, providing a very wide window of potential applications.