Computational self-assembly of a one-component icosahedral quasicrystal

Nearly Ideal Random Tiling with Dodecagonal Symmetry from Pentablock Quarterpolymers of the AB1CB2D Type

Two-dimensional tiling manners as cross-sectional views of cylindrical domain assembly formed by pentablock quarterpolymers of the AB1CB2D type in bulk were investigated. Several binary and ternary blends from three mother polymers having different ϕB1/ϕB2 ratios (ϕB1 and ϕB2 are the volume fractions of the B1 and B2 blocks, respectively) represent nonperiodic but ordered triangle/square tilings, where the N3/N4 ratios (N3 and N4 are the numbers of triangles and squares in the observed area, respectively) are all close enough to the theoretical value of 4/√3 ≑ 2.31 for the dodecagonal quasicrystalline (DDQC) state, irrespective of the total number of polygons. The TEM images, having almost the same N3/N4 ratios, were proved to show 4- and 6-fold symmetries in terms of the angular appearance of equilateral polygon sides via image analyses. Among them, a ternary blend showed a nearly ideal random tiling pattern that is almost equivalent to the theoretically predicted tiling by SCFT. Moreover, the magnitude of phason strain estimated for a TEM image from the ternary blend was proved to be quite small when the observing area is narrow, while it deviates from the ideal quasicrystalline tiling with an increasing number of vertices in the observing area.

Local Structural Features Elucidate Crystallization of Complex Structures

Complex crystal structures are composed of multiple local environments, and how this type of order emerges spontaneously during crystal growth has yet to be fully understood. We study crystal growth across various structures and along different crystallization pathways, using self-assembly simulations of identical particles that interact via multiwell isotropic pair potentials. We apply an unsupervised machine learning method to features from bond-orientational order metrics to identify different local motifs present during a given structure's crystallization process. In this manner, we distinguish different crystallographic sites in highly complex structures. Tailoring this order parameter to structures of varying complexity and coordination number, we study the emergence of local order along a multistep crystal growth pathway─from a low-density fluid to a high-density, supercooled amorphous liquid droplet and to a bulk crystal. We find a consistent under-coordination of the liquid relative to the average coordination number in the bulk crystal. We use our order parameter to analyze the geometrically frustrated growth of a Frank-Kasper phase and discover how structural defects compete with the formation of crystallographic sites that are more high-coordinated than the liquid environments. The method presented here for classifying order on a particle-by-particle level has broad applicability to future studies of structural self-assembly and crystal growth, and they can aid in the design of building blocks and for targeting pathways of formation of soft-matter structures.

Quasicrystal of Binary Hard Spheres on a Plane Stabilized by Configurational Entropy

Because of their aperiodic nature, quasicrystals are one of the least understood phases in statistical physics. One significant complication they present in comparison to their periodic counterparts is the fact that any quasicrystal can be realized as an exponentially large number of different tilings, resulting in a significant contribution to the quasicrystal entropy. Here, we use free-energy calculations to demonstrate that it is this configurational entropy which stabilizes a dodecagonal quasicrystal in a binary mixture of hard spheres on a plane. Our calculations also allow us to quantitatively confirm that in this system all tiling realizations are essentially equally likely, with free-energy differences less than 0.0001k_{B}T per particle-an observation that could be related to the observation of only random tilings in soft-matter quasicrystals. Owing to the simplicity of the model and its available counterparts in colloidal experiments, we believe that this system is an excellent candidate to achieve the long-awaited quasicrystal self-assembly on the micron scale.

Learning to learn by using nonequilibrium training protocols for adaptable materials

Evolution in time-varying environments naturally leads to adaptable biological systems that can easily switch functionalities. Advances in the synthesis of environmentally responsive materials therefore open up the possibility of creating a wide range of synthetic materials which can also be trained for adaptability. We consider high-dimensional inverse problems for materials where any particular functionality can be realized by numerous equivalent choices of design parameters. By periodically switching targets in a given design algorithm, we can teach a material to perform incompatible functionalities with minimal changes in design parameters. We exhibit this learning strategy for adaptability in two simulated settings: elastic networks that are designed to switch deformation modes with minimal bond changes and heteropolymers whose folding pathway selections are controlled by a minimal set of monomer affinities. The resulting designs can reveal physical principles, such as nucleation-controlled folding, that enable such adaptability.

Targeted Discovery of Low-Coordinated Crystal Structures via Tunable Particle Interactions

Particles interacting via isotropic, multiwell pair potentials have been shown to self-assemble into a range of crystal structures, yet how the characteristics of the underlying interaction potential give rise to the resultant structure remains largely unknown. We have thus developed a functional form for the interaction potential in which all features can be tuned independently. We perform continuous parameter space searches by systematically changing pairs of parameters, controlling the various features of the interaction potential. By enforcing a repulsive first well (controlling particle interactions of the first neighbor shell), we stimulate the formation of low-coordinated assemblies. We report the self-assembly of 20 previously unknown crystal structure types, 14 of which have low coordination numbers. Despite limiting the search to a small region of the vast parameter space of possible particle interactions, a wealth of complexity and symmetry is apparent within these crystal structures, which include clathrates with empty cages and low-symmetry structures. Our findings suggest that an unknown number of previously undiscovered crystal structure configurations are possible through self-assembly, which can serve as interesting design targets for soft condensed matter synthesis.

Self-assembly of dodecagonal and octagonal quasicrystals in hard spheres on a plane

Hard spheres are one of the most fundamental model systems in soft matter physics, and have been instrumental in shedding light on nearly every aspect of classical condensed matter. Here, we add one more important phase to the list that hard spheres form: quasicrystals. Specifically, we use simulations to show that an extremely simple, purely entropic model system, consisting of two sizes of hard spheres resting on a flat plane, can spontaneously self-assemble into two distinct random-tiling quasicrystal phases. The first quasicrystal is a dodecagonal square-triangle tiling, commonly observed in a large variety of colloidal systems. The second quasicrystal has, to our knowledge, never been observed in either experiments or simulations. It exhibits octagonal symmetry, and consists of three types of tiles: triangles, small squares, and large squares, whose relative concentration can be continuously varied by tuning the number of smaller spheres present in the system. The observed tile composition of the self-assembled quasicrystals agrees very well with the theoretical prediction we obtain by considering the four-dimensional (lifted) representation of the quasicrystal. Both quasicrystal phases form reliably and rapidly over a significant part of parameter space. Our results demonstrate that entropy combined with a set of geometrically compatible, densely packed tiles can be sufficient ingredients for the self-assembly of colloidal quasicrystals.

A columnar liquid quasicrystal with a honeycomb structure that consists of triangular, square and trapezoidal cells

Quasicrystals are intriguing structures that have long-range positional correlations but no periodicity in real space, and typically with rotational symmetries that are 'forbidden' in conventional periodic crystals. Here, we present a two-dimensional columnar liquid quasicrystal with dodecagonal symmetry. Unlike previous dodecagonal quasicrystals based on random tiling, a honeycomb structure based on a strictly quasiperiodic tessellation of tiles is observed. The structure consists of dodecagonal clusters made up of triangular, square and trapezoidal cells that are optimal for local packing. To maximize the presence of such dodecagonal clusters, the system abandons periodicity but adopts a quasiperiodic structure that follows strict packing rules. The stability of random-tiling dodecagonal quasicrystals is often attributed to the entropy of disordering when strict tiling rules are broken, at the sacrifice of the long-range positional order. However, our results demonstrate that quasicrystal stability may rest on energy minimization alone, or with only minimal entropic intervention.

Formation and fluctuation of two-dimensional dodecagonal quasicrystals

The self-assembly of two-dimensional dodecagonal quasicrystals (DDQCs) from patchy particles is investigated by Brownian dynamics simulations. The patchy particle has a five-fold rotational symmetry pattern described by the spherical harmonics Y₅₅. From the formation of the DDQC obtained by an annealing process, we find the following mechanism. The early stage of the dynamics is dominated by hexagonal structures. Then, nucleation of dodecagonal motifs appears by particle rearrangement, and finally the motifs span the whole system. The transition from the hexagonal structure into the dodecagonal motif is coincident with the collective motion of the particles. The DDQC consists of clusters of dodecagonal motifs, which can be classified into several packing structures. By the analyses of the DDQC under fixed temperature, we find that the fluctuations are characterised by changes in the network of the dodecagonal motifs. Finally we compare the DDQCs assembled from the patchy particle system and isotropic particle system. The two systems share a similar mechanism of the formation and fluctuation of DDQCs.

Freezing in two-length-scale systems: complexity, universality and prediction

Two-length-scale pair potentials arise ubiquitously in condensed matter theory as effective interparticle interactions in molecular, metallic and soft matter systems. The existence of two different bond lengths generated by the shape of potential causes complicated behavior in even one-component systems: polymorphism in solid and liquid states, water-like anomalies, the formation of quasicrystals and high stability against crystallization. Here we address general properties of freezing in one-component two-length-scale systems and argue that solidification of a liquid during cooling is essentially determined by the radial distribution function (RDF) of the liquid. We show that different two-length-scale systems having similar RDFs freeze into the same solid phases. In some cases, the similarity between RDFs can be expressed by the proximity of two dimensionless effective parameters: the ratio between effective bond lengths,λ, and the fraction of short-bonded particlesφ. We validate this idea by studying the formation of different solid phases in different two-length-scale systems. The method proposed allows predicting effectively the formation of solid phases in both numerical simulations and self-assembling experiments in soft matter systems with tunable interactions.

Using Machine Learning to Predict and Understand Complex Self-Assembly Behaviors of a Multicomponent Nanocomposite

Blends of nanoparticles, polymers, and small molecules can self-assemble into optical, magnetic, and electronic devices with structure-dependent properties. However, the relationship between a multicomponent nanocomposite's formulation and its assembled structure is complex and cannot be predicted by theory. The blends can be strongly influenced by processing conditions, which can introduce non-equilibrium states. Currently, nanocomposite devices are designed through cycles of experimental trial and error. Machine learning (ML) methods are a compelling alternative because they can use existing datasets to map high-dimensional spaces. These methods do not rely on known relationships between parameters, so they are suited to complex systems without a solid theoretical foundation. Here, a dataset of 595 microscopy images of nanocomposite thin films is used to train a series of ML models. Correlations between the input and output parameters are examined, providing new insights into the system. Finally, the most successful ML model is used to predict the structures of new nanocomposite compositions. The results confirm that ML techniques can be used to improve the efficiency of nanocomposite device design. More broadly, the current study suggests some of the advantages and challenges associated with applying ML to complex systems.

The Largest Quasicrystalline Tiling with Dodecagonal Symmetry from a Single Pentablock Quarterpolymer of the AB1CB2D Type

A quasicrystalline tiling pattern with tile size of ca. 60 nm has been discovered in the bulk state of a four-component pentablock polymer molecule of the AS₁IS₂P type, where A, S, I, and P denote poly(4-vinylbenzyldimethylamine), polystyrene, polyisoprene, and poly(2-vinylpyridine), respectively. The polymer samples used were prepared by anionic polymerizations and have narrow molecular weight distribution. The sample films were obtained by an extremely slow solvent-cast process from dilute solutions of tetrahydrofuran for 14 days. It has been found by TEM observation that the quarterpolymer, AS₁IS₂P-4 (M_n = 149 kg/mol, ϕ_A/ϕ_S1/ϕ_I/ϕ_S2/ϕ_P = 0.12/0.27/0.20/0.29/0.13), reveals two final stable structures, i.e., a 3.3.4.3.4 periodic tiling pattern as a minor component and a quasicrystalline (QC) tiling with dodecagonal symmetry as a major component, where the former includes a triangle/square number ratio of 2 and the latter has one of approximately 2.28, which is close enough to the ideal ratio, 4/ 3 ≑ 2.31, for the triangle/square random tiling of the dodecagonal QC tiling (DDQC). Two structures were also clearly proved by SAXS diffraction patterns. Here, it should be noted this QC structure having a tile side length of ca. 60 nm was created with a single block polymer molecule.

A bottom-up design strategy for controllable self-assembly based on the isotropic double-well potential

Controllable self-assembly of particles or atoms is still challenging in the synthesis of materials with desirable properties that are highly relevant to the microscopic structures determined by the interparticle interactions. To gain insight into how the interactions affect the self-assembly, we designed various kinds of isotropic double-well potentials and simulated the motion of the particles. By controlling the depth and location of the potential wells and the height of the barriers, we studied their effects on the aggregation structures and the related microscopic kinetic processes. We identified five aggregation patterns at different temperatures and eight kinds of crystals, including Frank-Kasper phases, and observed the expansion or contraction of crystals. We found that the system usually stays in a sparse configuration at very low or very high temperatures. The particles typically assemble into a loosely packed cluster at medium temperatures and then deplete into a tightly packed state with a specific pattern. These phenomena can be explained from the perspective of energy. In contrast, very few structures could be obtained for the system guided by a single-well potential under the same simulation conditions. Thus, the interparticle interactions driven by the double-well potential greatly enrich the possible packing morphology of the system. The information obtained from this work helps us to understand how to achieve a specific self-assembled architecture through a reasonable selection of materials.

Molecular dynamics investigation of a one-component model for the stacking motif in complex alloy structures

Developing realistic three-dimensional growth models for quasicrystals is a fundamental requirement. The present work employs classical molecular dynamics simulations to investigate the adsorption of Al on a close-packed Al layer containing atomic vacancies. Simulation results show that the adsorbed Al atoms are located preferentially above and below the atomic vacancies in the close-packed layer, and the results obtained from a one-component system of atoms interacting via an interatomic pair potential for Al–Al appropriately reproduce the stacking motif seen in complex alloys such as the μ-Al4Mn phase. The simulations also reveal the formation of a deformed icosahedron. These results provide new insights into the growth mechanism and origin of complex alloys and quasicrystals.

DNA self-organization controls valence in programmable colloid design

Just like atoms combine into molecules, colloids can self-organize into predetermined structures according to a set of design principles. Controlling valence-the number of interparticle bonds-is a prerequisite for the assembly of complex architectures. The assembly can be directed via solid "patchy" particles with prescribed geometries to make, for example, a colloidal diamond. We demonstrate here that the nanoscale ordering of individual molecular linkers can combine to program the structure of microscale assemblies. Specifically, we experimentally show that covering initially isotropic microdroplets with N mobile DNA linkers results in spontaneous and reversible self-organization of the DNA into Z(N) binding patches, selecting a predictable valence. We understand this valence thermodynamically, deriving a free energy functional for droplet-droplet adhesion that accurately predicts the equilibrium size of and molecular organization within patches, as well as the observed valence transitions with N Thus, microscopic self-organization can be programmed by choosing the molecular properties and concentration of binders. These results are widely applicable to the assembly of any particle with mobile linkers, such as functionalized liposomes or protein interactions in cell-cell adhesion.

Effective potentials induced by mixtures of patchy and hard co-solutes

The addition of co-solutes to colloidal suspensions is often employed to induce tunable depletion interactions. In this work, we investigate effective colloidal interactions arising from binary co-solute mixtures of hard spheres and patchy particles. By changing the relative concentration of the two species, we show that the resulting effective potential V_eff continuously changes from the one obtained for a single-component hard sphere co-solute to that mediated by the single-component patchy particle co-solute. Interestingly, we find that, independent of the relative concentration of the two components, the resulting V_eff is additive, i.e., it is well-described by the linear combination of the effective interactions mediated by respective pure co-solutes. However, a breakdown of the additivity occurs when the co-solute mixture is close to the onset of a demixing transition. These results represent a step forward in understanding and predicting colloidal behavior in complex and crowded environments and for exploiting this knowledge to design targeted colloidal superstructures.

How to design an icosahedral quasicrystal through directional bonding

Icosahedral quasicrystals (IQCs) are materials that exhibit long-range order but lack periodicity in any direction. Although IQCs were the first reported quasicrystals¹, they have been experimentally observed only in metallic alloys², not in other materials. By contrast, quasicrystals with other symmetries (particularly dodecagonal) have now been found in several soft-matter systems^3-5. Here we introduce a class of IQCs built from model patchy colloids that could be realized experimentally using DNA origami particles. Our rational design strategy leads to systems that robustly assemble in simulations into a target IQC through directional bonding. This is illustrated for both body-centred and primitive IQCs, with the simplest systems involving just two particle types. The key design feature is the geometry of the interparticle interactions favouring the propagation of an icosahedral network of bonds, despite this leading to many particles not being fully bonded. As well as furnishing model systems in which to explore the fundamental physics of IQCs, our approach provides a potential route towards functional quasicrystalline materials.

Entropic formation of a thermodynamically stable colloidal quasicrystal with negligible phason strain

Quasicrystals have been discovered in a variety of materials ranging from metals to polymers. Yet, why and how they form is incompletely understood. In situ transmission electron microscopy of alloy quasicrystal formation in metals suggests an error-and-repair mechanism, whereby quasiperiodic crystals grow imperfectly with phason strain present, and only perfect themselves later into a high-quality quasicrystal with negligible phason strain. The growth mechanism has not been investigated for other types of quasicrystals, such as dendrimeric, polymeric, or colloidal quasicrystals. Soft-matter quasicrystals typically result from entropic, rather than energetic, interactions, and are not usually grown (either in laboratories or in silico) into large-volume quasicrystals. Consequently, it is unknown whether soft-matter quasicrystals form with the high degree of structural quality found in metal alloy quasicrystals. Here, we investigate the entropically driven growth of colloidal dodecagonal quasicrystals (DQCs) via computer simulation of systems of hard tetrahedra, which are simple models for anisotropic colloidal particles that form a quasicrystal. Using a pattern recognition algorithm applied to particle trajectories during DQC growth, we analyze phason strain to follow the evolution of quasiperiodic order. As in alloys, we observe high structural quality; DQCs with low phason strain crystallize directly from the melt and only require minimal further reduction of phason strain. We also observe transformation from a denser approximant to the DQC via continuous phason strain relaxation. Our results demonstrate that soft-matter quasicrystals dominated by entropy can be thermodynamically stable and grown with high structural quality--just like their alloy quasicrystal counterparts.

Programming patchy particles to form three-dimensional dodecagonal quasicrystals

Model patchy particles have been shown to be able to form a wide variety of structures, including symmetric clusters, complex crystals, and even two-dimensional quasicrystals. Here, we investigate whether we can design patchy particles that form three-dimensional quasicrystals, in particular targeting a quasicrystal with dodecagonal symmetry that is made up of stacks of two-dimensional quasicrystalline layers. We obtain two designs that are able to form such a dodecagonal quasicrystal in annealing simulations. The first is a one-component system of seven-patch particles but with wide patches that allow them to adopt both seven- and eight-coordinated environments. The second is a ternary system that contains a mixture of seven- and eight-patch particles and is likely to be more realizable in experiments, for example, using DNA origami. One interesting feature of the first system is that the resulting quasicrystals very often contain a screw dislocation.

Effective Pair Interactions and Structure in Liquid Noble Metals within Wills-Harrison and Bretonnet-Silbert Models

Recently, for calculating the effective pair interactions in liquid transition metals, we have developed an approach which includes the Wills-Harrison and Bretonnet-Silbert models as limit cases. Here, we apply this approach to noble liquid metals. The dependencies of pair potentials and corresponding MD-simulated pair correlation functions in pure liquid Cu, Ag and Au on the portion of the non-diagonal (with respect to the magnet quantum number) d-d-electron couplings in the metal under consideration are studied. The model provides a good agreement with experimental and ab initio data for pair correlation functions, structure factors and velocity autocorrelation functions.

Self-assembly of shape-tunable oblate colloidal particles into orientationally ordered crystals, glassy crystals and plastic crystals

The shapes of colloidal particles are crucial to self-assembled superstructures. Understanding the relationship between the shapes of building blocks and the resulting crystal structures is an important fundamental question. Here, we demonstrate that, by using particles whose shape interpolates between a flat disc and a sphere, not only are self-assembled superstructures but also their orientations sensitively dependent on the particle shape. By changing the shape gradually from a flat disc to a spherical shape, a crystal sequence from orientationally ordered crystals to orientationally disordered crystals with frozen and more free rotations are found. The latter two phases are identified as a glassy crystal and a plastic crystal, respectively. By combining theoretical model calculations, the formed crystal structures and the occurring transitions are found to be dictated by the interplay between particle shape and particle-particle interaction as well as particle-wall interaction. In particular, for quasi-spherical shapes, when the strong attraction dominates, a glassy crystal forms, or otherwise a plastic crystal forms. These results demonstrate that the interplay between the particle shape and the interaction can be used to tune crystallization and further fabricate colloid-based new structured and dynamic materials.

Two-dimensional oxide quasicrystal approximants with tunable electronic and magnetic properties

Recently, the discovery of the quasiperiodic order in ultra-thin perovskite films reinvigorated the field of 2-dimensional oxides on metals, and raised the question of the reasons behind the emergence of the quasiperiodic order in these systems. The effect of size-mismatch between the two separate systems has been widely reported as a key factor governing the formation of new oxide structures on metals. Herein, we show that electronic effects can play an important role as well. To this end, the structural, thermodynamic, electronic and magnetic properties of freestanding two-dimensional oxide quasicrystalline approximants and their characteristics when deposited over metallic substrates are systematically investigated to unveil the structure-property relationships within the series. Our thermodynamic approach suggests that the formation of these aperiodic systems is likely for a wide range of compositions. In addition, the magnetic properties and work functions of the thin films can be controlled by tuning their chemical composition. This work provides well-founded general insights into the driving forces behind the emergence of the quasiperiodic order in ternary oxides grown on elemental metals and offers guidelines for the discovery of new oxide quasicrystalline ultra-thin films with interesting physical properties.

Moving beyond the constraints of chemistry via crystal structure discovery with isotropic multiwell pair potentials

The rigid constraints of chemistry-dictated by quantum mechanics and the discrete nature of the atom-limit the set of observable atomic crystal structures. What structures are possible in the absence of these constraints? Here, we systematically crystallize one-component systems of particles interacting with isotropic multiwell pair potentials. We investigate two tunable families of pairwise interaction potentials. Our simulations self-assemble a multitude of crystal structures ranging from basic lattices to complex networks. Sixteen of the structures have natural analogs spanning all coordination numbers found in inorganic chemistry. Fifteen more are hitherto unknown and occupy the space between covalent and metallic coordination environments. The discovered crystal structures constitute targets for self-assembly and expand our understanding of what a crystal structure can look like.

Long-range order in quadrupolar systems on spherical surfaces

The interplay between curvature, confinement and ordering on curved manifolds, with anisotropic interactions between building blocks, takes a central role in many fields of physics. In this paper, we investigate the effects of lattice symmetry and local positional order on orientational ordering in systems of long-range interacting point quadrupoles on a sphere in the zero temperature limit. Locally triangular spherical lattices show long-range ordered quadrupolar configurations only for specific symmetric lattices as strong geometric frustration prevents general global ordering. Conversely, the ground states on Caspar-Klug lattices are more diverse, with many different symmetries depending on the position of quadrupoles within the fundamental domain. We also show that by constraining the quadrupole tilts with respect to the surface normal, which models interactions with the substrate, and by considering general quadrupole tensors, we can manipulate the ground state configuration symmetry.

Interaction Patterns for Staggered Assembly of Fibrils from Semiflexible Chains

The design of colloidal interactions to achieve target self-assembled structures has especially been done for compact objects such as spheres with isotropic interaction potentials, patchy spheres and other compact objects with patchy interactions. Inspired by the self-assembly of collagen-I fibrils and intermediate filaments, we here consider the design of interaction patterns on semiflexible chains that could drive their staggered assembly into regular (para)crystalline fibrils. We consider semiflexible chains composed of a finite number of types of interaction beads (uncharged hydrophilic, hydrophobic, positively charged and negatively charged) and optimize the sequence of these interaction beads with respect to the interaction energy of the semiflexible chains in a number of target-staggered crystalline packings. We find that structures with the lowest interaction energies, that form simple lattices, also have low values of L/D (where L is chain length and D is stagger). In the low interaction energy sequences, similar types of interaction beads cluster together to form stretches. Langevin Dynamics simulations confirm that semiflexible chains with optimal sequences self-assemble into the designed staggered (para)crystalline fibrils. We conclude that very simple interaction patterns should suffice to drive the assembly of long semiflexible chains into staggered (para)crystalline fibrils.

Dynamic Observation of Dendritic Quasicrystal Growth upon Laser-Induced Solid-State Transformation

We report the laser-induced solid-state transformation between a periodic "approximant" and quasicrystal in the Al-Cr system during rapid quenching. Dynamic transmission electron microscopy allows us to capture in situ the dendritic growth of the metastable quasicrystals. The formation of dendrites during solid-state transformation is a rare phenomenon, which we attribute to the structural similarity between the two intermetallics. Through ab initio molecular dynamics simulations, we identify the dominant structural motif to be a 13-atom icosahedral cluster transcending the phases of matter.

Designing disorder into crystalline materials

Crystals are a state of matter characterized by periodic order. Yet, crystalline materials can harbour disorder in many guises, such as non-repeating variations in composition, atom displacements, bonding arrangements, molecular orientations, conformations, charge states, orbital occupancies or magnetic structure. Disorder can sometimes be random but, more usually, it is correlated. Frontier research into disordered crystals now seeks to control and exploit the unusual patterns that persist within these correlated disordered states in order to access functional responses inaccessible to conventional crystals. In this Review, we survey the core design principles that guide targeted control over correlated disorder. We show how these principles - often informed by long-studied statistical mechanical models - can be applied across an unexpectedly broad range of materials, including organics, supramolecular assemblies, oxide ceramics and metal-organic frameworks. We conclude with a forward-looking discussion of the exciting link between disorder and function in responsive media, thermoelectrics and topological phases.

Fabrication of magnetic colloidal bicones and their switchable self-assembly

The creation and fabrication of novel shapes are of paramount importance for colloids to self-assemble into desired structured materials and in order to use them as model systems for fundamental studies. Here, we demonstrate a simple yet effective method to fabricate bicone-shaped colloids by using SiO₂/PS core/shell particles through a thermal stretching route. The shape parameters are tunable. The aspect ratio is found to be predictable, and the experimental results are consistent with the model calculations. The yields are high enough to be used for self-assembly studies. As an illustration, we investigate the phase behavior of particles when interacting with depletion forces and show that the particles can form hexagonal or non-hexagonal crystal lattices in quasi-two-dimensional space. The method also allows functional nanoparticles to be integrated into the cores, resulting in responsive colloidal bicones. Moreover, the magnetically responsive self-assembly of particles is demonstrated.

Evolution of local motifs and topological proximity in self-assembled quasi-crystalline phases

Using methods from the field of topological data analysis, we investigate the self-assembly and emergence of three-dimensional quasi-crystalline structures in a single-component colloidal system. Combining molecular dynamics and persistent homology, we analyse the time evolution of persistence diagrams and particular local structural motifs. Our analysis reveals the formation and dissipation of specific particle constellations in these trajectories, and shows that the persistence diagrams are sensitive to nucleation and convergence to a final structure. Identification of local motifs allows quantification of the similarities between the final structures in a topological sense. This analysis reveals a continuous variation with density between crystalline clathrate, quasi-crystalline, and disordered phases quantified by 'topological proximity', a visualization of the Wasserstein distances between persistence diagrams. From a topological perspective, there is a subtle, but direct connection between quasi-crystalline, crystalline and disordered states. Our results demonstrate that topological data analysis provides detailed insights into molecular self-assembly.

Phase behavior of a two-dimensional core-softened system: new physical insights

In this work, we report results of extensive computer simulations regarding the phase behavior of a core-softened system. By using structural and thermodynamic descriptors, as well as self-diffusion coefficients, we provide a comprehensive view of the rich phase behavior displayed by the particular instance of the model studied in here. Our calculations agree with previously published results focused on a smaller region in the temperature-density parameter space (Dudalov et al 2014 Soft Matter 10 4966). In this work, we explore a broader region in this parameter space, and uncover interesting fluid phases with low-symmetry local order, that were not reported by previous works. Solid phases were also found, and have been previously characterized in detail by (Kryuchkov et al 2018 Soft Matter 14 2152). Our results support previously reported findings, and provide new physical insights regarding the emergence of order as disordered phases transform into solids by providing radial distribution function maps and specific heat data. Our results are summarized in terms of a phase diagram.

Spontaneous Crystallization in Systems of Binary Hard Sphere Colloids

Computer simulations of the fluid-to-solid phase transition in the hard sphere system were instrumental for our understanding of crystallization processes. But while colloid experiments and theory have been predicting the stability of several binary hard sphere crystals for many years, simulations were not successful to confirm this phenomenon. Here, we report the growth of binary hard sphere crystals isostructural to Laves phases, AlB_{2}, and NaZn_{13} in simulation directly from the fluid. We analyze particle kinetics during Laves phase growth using event-driven molecular dynamics simulations with and without swap moves that speed up diffusion. The crystallization process transitions from nucleation and growth to spinodal decomposition already deep within the fluid-solid coexistence regime. Finally, we present packing fraction-size ratio state diagrams in the vicinity of the stability regions of three binary crystals.

Statistical reprogramming of macroscopic self-assembly with dynamic boundaries

Self-assembly is a ubiquitous process that can generate complex and functional structures via local interactions among a large set of simpler components. The ability to program the self-assembly pathway of component sets elucidates fundamental physics and enables alternative competitive fabrication technologies. Reprogrammability offers further opportunities for tuning structural and material properties but requires reversible selection from multistable self-assembling patterns, which remains a challenge. Here, we show statistical reprogramming of two-dimensional (2D), noncompact self-assembled structures by the dynamic confinement of orbitally shaken and magnetically repulsive millimeter-scale particles. Under a constant shaking regime, we control the rate of radius change of an assembly arena via moving hard boundaries and select among a finite set of self-assembled patterns repeatably and reversibly. By temporarily trapping particles in topologically identified stable states, we also demonstrate 2D reprogrammable stiffness and three-dimensional (3D) magnetic clutching of the self-assembled structures. Our reprogrammable system has prospective implications for the design of granular materials in a multitude of physical scales where out-of-equilibrium self-assembly can be realized with different numbers or types of particles. Our dynamic boundary regulation may also enable robust bottom-up control strategies for novel robotic assembly applications by designing more complex spatiotemporal interactions using mobile robots.

The atomic structure of the Bergman-type icosahedral quasicrystal based on the Ammann-Kramer-Neri tiling

In this study, the atomic structure of the ternary icosahedral ZnMgTm quasicrystal (QC) is investigated by means of single-crystal X-ray diffraction. The structure is found to be a member of the Bergman QC family, frequently found in Zn-Mg-rare-earth systems. The ab initio structure solution was obtained by the use of the Superflip software. The infinite structure model was founded on the atomic decoration of two golden rhombohedra, with an edge length of 21.7 Å, constituting the Ammann-Kramer-Neri tiling. The refined structure converged well with the experimental diffraction diagram, with the crystallographic R factor equal to 9.8%. The Bergman clusters were found to be bonded by four possible linkages. Only two linkages, b and c, are detected in approximant crystals and are employed to model the icosahedral QCs in the cluster approach known for the CdYb Tsai-type QC. Additional short b and a linkages are found in this study. Short interatomic distances are not generated by those linkages due to the systematic absence of atoms and the formation of split atomic positions. The presence of four linkages allows the structure to be pictured as a complete covering by rhombic triacontahedral clusters and consequently there is no need to define the interstitial part of the structure (i.e. that outside the cluster). The 6D embedding of the solved structure is discussed for the final verification of the model.

Dodecagonal quasicrystal silicene: preparation, mechanical property, and friction behaviour

In this study, we obtained dodecagonal monolayer silicene with three-fold and four-fold coordination by melt quenching via molecular dynamics (MD) simulations. Stretching simulation of the pre-strained dodecagonal silicene showed lower critical stress than the honeycomb silicene and resulted in an increase in six-fold rings during the plastic deformation since the four-coordinated atom sites are less mechanically favoured than the three-coordinated sites. The friction behaviours with an AFM tip sliding on the dodecagonal and honeycomb surfaces under different loads and tip sizes were simulated and compared. For all the investigated cases, the dodecagonal surface always showed a lower mean friction force than the honeycomb surface. The lower friction of the quasicrystal was observed, and the mechanism was illuminated successfully for the first time by MD simulations. The reduced friction of dodecagonal silicene can be explained by the morphology of the one-dimensional potential energy surface (PES). The 1D PES of dodecagonal silicene has longer potential corrugation lengths than honeycomb silicene, which induce mild motion of the tip in the stick process and lower friction force. Considering the close density of the employed dodecagonal and honeycomb structure, the longer potential corrugation length is a consequence of the quasiperiodic morphology rather than the interspace between atoms. Besides, with a larger tip size, the 1D PES on the dodecagonal surface has a flatter area, which contributes further to the reduced friction force on the dodecagonal surface.

Hierarchical complex self-assembly in binary nanoparticle mixtures

Hierarchical self-assembly of soft matter provides a powerful route to create complex materials with enhanced physical properties. The understanding of the fundamental processes leading to such organization can provide design rules to create new functional materials. In this work, we use a simple model of polymer-grafted nanoparticles to explore the self-assembly of binary mixtures. By using Monte Carlo simulations we study the interplay of composition, density and particle sizes on the self-organization of such nanoparticle systems. It is found that complex hierarchical organization can take place for conditions where one-component systems form simple lattices. In particular, a mixture where one component forms a structure with 18-fold symmetry in a sea of an apparent disordered phase of the second component is observed to emerge for certain parameter combinations.

Anomalous behavior and structure of a liquid of particles interacting through the harmonic-repulsive pair potential near the crystallization transition

A characteristic property of many soft matter systems is an ultrasoft effective interaction between their structural units. This softness often leads to complex behavior. In particular, ultrasoft systems under pressure demonstrate polymorphism of complex crystal and quasicrystal structures. Therefore, it is of interest to investigate how different can be the structure of the fluid state in such systems at different pressures. Here we address this issue for a model liquid composed of particles interacting through the harmonic-repulsive pair potential. This system can form different crystal structures as the liquid is cooled. We find that, at certain pressures, the liquid exhibits unusual properties, such as a negative thermal expansion coefficient. Besides, the volume and the potential energy of the system can increase during crystallization. At certain pressures, the system demonstrates high stability against crystallization and it is hardly possible to crystallize it on the timescales of the simulations. To address the liquid's structure at high pressures, we consider the scaled pair distribution function (PDF) and the bond-orientational order (BOO) parameters. The marked change happening with the PDF, as pressure increases, is the splitting of the first peak which is caused by the appearance of non-negligible interactions with the second neighbors and the following rearrangement of the structure. Our findings suggest that non-trivial effects, usually explained by different interactions at different spatial scales, can also be observed in one-component systems with simple one-length-scale ultrasoft repulsive interactions.

Complex Self-Assembly from Simple Interaction Rules in Model Colloidal Mixtures

Building structures with hierarchical order through the self-assembly of smaller blocks is not only a prerogative of nature, but also a strategy to design artificial materials with tailored functions. We explore in simulation the spontaneous assembly of colloidal particles into extended structures, using spheres and size-asymmetric dimers as solute particles, while treating the solvent implicitly. Besides rigid cores for all particles, we assume an effective short-range attraction between spheres and small monomers to promote, through elementary rules, dimer-mediated aggregation of spheres. Starting from a completely disordered configuration, we follow the evolution of the system at low temperature and density, as a function of the relative concentration of the two species. When spheres and large monomers are of same size, we observe the onset of elongated aggregates of spheres, either disconnected or cross-linked, and a crystalline bilayer. As spheres grow bigger, the self-assembling scenario changes, getting richer overall, with the addition of flexible membrane sheets with crystalline order and monolayer vesicles. With this wide assortment of structures, our model can serve as a viable template to achieve a better control of self-assembly in dilute suspensions of microsized particles.

Which Wave Numbers Determine the Thermodynamic Stability of Soft Matter Quasicrystals?

For soft matter to form quasicrystals an important ingredient is to have two characteristic length scales in the interparticle interactions. To be more precise, for stable quasicrystals, periodic modulations of the local density distribution with two particular wave numbers should be favored, and the ratio of these wave numbers should be close to certain special values. So, for simple models, the answer to the title question is that only these two ingredients are needed. However, for more realistic models, where in principle all wave numbers can be involved, other wave numbers are also important, specifically those of the second and higher reciprocal lattice vectors. We identify features in the particle pair interaction potentials that can suppress or encourage density modes with wave numbers associated with one of the regular crystalline orderings that compete with quasicrystals, enabling either the enhancement or suppression of quasicrystals in a generic class of systems.

Assembly of clathrates from tetrahedral patchy colloids with narrow patches

Here, we revisit the assembly of colloidal tetrahedral patchy particles. Previous studies have shown that the crystallization of diamond from the fluid phase depends more critically on patch width than on the interaction range: particles with patches narrower than 40° crystallize readily and those with wide patches form disordered glass states. We find that the crystalline structure formed from the fluid also depends on the patch width. Whereas particles with intermediate patches assemble into diamond (random stacking of cubic and hexagonal diamond layers), particles with narrow patches (with width ≈20° or less) crystallize frequently into clathrates. Free energy calculations show that clathrates are never (in the pressure-temperature plane) thermodynamically more stable than diamond. The assembly of clathrate structures is thus attributed to kinetic factors that originate from the thermodynamic stabilization of pentagonal rings with respect to hexagonal ones as patches become more directional. These pentagonal rings present in the fluid phase assemble into sII clathrate or into large clusters containing 100 particles and exhibiting icosahedral symmetry. These clusters then grow by interpenetration. Still, the organization of these clusters into extended ordered structures was never observed in the simulations.