Structural phases of colloids interacting via a flat-well potential

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.

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.

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.

Crystal, Fivefold and Glass Formation in Clusters of Polymers Interacting with the Square Well Potential

We present results, from Monte Carlo (MC) simulations, on polymer systems of freely jointed chains with spherical monomers interacting through the square well potential. Starting from athermal packings of chains of tangent hard spheres, we activate the square well potential under constant volume and temperature corresponding effectively to instantaneous quenching. We investigate how the intensity and range of pair-wise interactions affected the final morphologies by fixing polymer characteristics such as average chain length and tolerance in bond gaps. Due to attraction chains are brought closer together and they form clusters with distinct morphologies. A wide variety of structures is obtained as the model parameters are systematically varied: weak interactions lead to purely amorphous clusters followed by well-ordered ones. The latter include the whole spectrum of crystal morphologies: from virtually perfect hexagonal close packed (HCP) and face centered cubic (FCC) crystals, to random hexagonal close packed layers of single stacking direction of alternating HCP and FCC layers, to structures of mixed HCP/FCC character with multiple stacking directions and defects in the form of twins. Once critical values of interaction are met, fivefold-rich glassy clusters are formed. We discuss the similarities and differences between energy-driven crystal nucleation in thermal polymer systems as opposed to entropy-driven phase transition in athermal polymer packings. We further calculate the local density of each site, its dependence on the distance from the center of the cluster and its correlation with the crystallographic characteristics of the local environment. The short- and long-range conformations of chains are analyzed as a function of the established cluster morphologies.

Two-dimensional crystals of star polymers: a tale of tails

The formation of non-hexagonal crystalline structures by the organisation of colloidal nanoparticles often involves the use of complex particles with anisotropic shape or interactions or the imposition of non-uniform external fields. Here we explore how unusual symmetries can be created using experimentally realistic particles that interact through isotropic and purely repulsive potentials. In particular, we use simulations to explore the phase behavior of two-dimensional systems of star polymers. We uncover how the tail of the pair potential has a large role in dictating the phase behavior. Star polymers interacting in the far field with a Gaussian potential only form hexagonal phases, while an exponential tail gives rise to stable primitive oblique and honeycomb lattices. We identify the ratio in strength between long and short range interactions and the nature of the transition between these regimes as crucial parameters to predict when non-hexagonal crystals of star polymers can be stable. This leads to experimental design rules for creating star polymers which should exhibit unusual lattice formation.

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.

On the formation of stripe, sigma, and honeycomb phases in a core-corona system

Using Monte Carlo simulations and free-energy calculations, we investigate the phase behaviour of a two-dimensional core-corona system. We model this system as particles consisting of an impenetrable hard core of diameter σ_HD surrounded by a purely repulsive soft corona of diameter δ = 1.95σ_HD. At low densities, we observe the spontaneous formation of a phase with a stripe texture as well as a honeycomb-like phase driven by both energy and entropy considerations. At high densities, we find that a two-dimensional analogue of the periodic sigma phase, considered as an approximant of dodecagonal quasicrystals, is energetically stabilised with respect to two distinct dodecagonal quasicrystals, namely, a square-triangle tiling and a square-triangle-shield tiling. We also find the formation of stable hexagonal phases at three distinct density ranges, which are energetically driven, i.e. by minimising the overlap of coronas. Furthermore, our calculations show that the low-density dodecagonal quasicrystal that was previously reported by Dotera et al., [Nature, 2014, 506, 208] is kinetically formed in the coexistence region between the honeycomb and the medium-density hexagonal phase.

The effect of temperature, interaction range, and pair potential on the formation of dodecagonal quasicrystals in core-corona systems

A two-dimensional dodecagonal quasicrystal was previously reported by Dotera et al (2014 Nature 506 208) in a system of particles interacting with a hard core of diameter σ and a repulsive square shoulder of diameter [Formula: see text]. In the current work, we examine the formation of this quasicrystal using bond orientational order parameters, correlation functions and tiling distributions. We find that this dodecagonal quasicrystal forms from a fluid phase. We further study the effect of the width of the repulsive shoulder by simulating the system over a range of values of δ. For the range of densities and temperatures considered, we observe the formation of the dodecagonal quasicrystal between [Formula: see text] and [Formula: see text]. We also study the effect of shape of the interaction potential by simulating the system using three other interaction potentials with two length scales, namely hard-core plus a linear ramp, modified exponential, or Buckingham (exp-6) potential. We observe the presence of the quasicrystal in all three systems. However, depending on the shape of the potential, the formation of the quasicrystal takes place at lower temperatures (or higher interaction strengths). Using free-energy calculations, we demonstrate that the quasicrystal is thermodynamically stable in the square-shoulder and linear-ramp system.

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.

Probing the limitations of isotropic pair potentials to produce ground-state structural extremes via inverse statistical mechanics

Inverse statistical-mechanical methods have recently been employed to design optimized short-range radial (isotropic) pair potentials that robustly produce novel targeted classical ground-state many-particle configurations. The target structures considered in those studies were low-coordinated crystals with a high degree of symmetry. In this paper, we further test the fundamental limitations of radial pair potentials by targeting crystal structures with appreciably less symmetry, including those in which the particles have different local structural environments. These challenging target configurations demanded that we modify previous inverse optimization techniques. In particular, we first find local minima of a candidate enthalpy surface and determine the enthalpy difference ΔH between such inherent structures and the target structure. Then we determine the lowest positive eigenvalue λ(0) of the Hessian matrix of the enthalpy surface at the target configuration. Finally, we maximize λ(0)ΔH so that the target structure is both locally stable and globally stable with respect to the inherent structures. Using this modified optimization technique, we have designed short-range radial pair potentials that stabilize the two-dimensional kagome crystal, the rectangular kagome crystal, and rectangular lattices, as well as the three-dimensional structure of the CaF(2) crystal inhabited by a single-particle species. We verify our results by cooling liquid configurations to absolute zero temperature via simulated annealing and ensuring that such states have stable phonon spectra. Except for the rectangular kagome structure, all of the target structures can be stabilized with monotonic repulsive potentials. Our work demonstrates that single-component systems with short-range radial pair potentials can counterintuitively self-assemble into crystal ground states with low symmetry and different local structural environments. Finally, we present general principles that offer guidance in determining whether certain target structures can be achieved as ground states by radial pair potentials.

Structural ordering of trapped colloids with competing interactions

The structure of colloids with competing interactions which are confined in a harmonic external trap potential is analyzed numerically by energy minimization in two spatial dimensions. A wealth of different cluster structures is found to be stable including clusters with a fringed outer rim (reminiscent to an ornamental border), clusters perforated with voids, as well as clusters with a crystalline core and a disordered rim. All cluster structures occur in a two-dimensional parameter space. The structural ordering can therefore be efficiently tuned by changing few parameters only providing access to a controlled fabrication of colloidal clusters.