Complex Crystals from Size-Disperse Spheres

Nanocrystal Assemblies: Current Advances and Open Problems

We explore the potential of nanocrystals (a term used equivalently to nanoparticles) as building blocks for nanomaterials, and the current advances and open challenges for fundamental science developments and applications. Nanocrystal assemblies are inherently multiscale, and the generation of revolutionary material properties requires a precise understanding of the relationship between structure and function, the former being determined by classical effects and the latter often by quantum effects. With an emphasis on theory and computation, we discuss challenges that hamper current assembly strategies and to what extent nanocrystal assemblies represent thermodynamic equilibrium or kinetically trapped metastable states. We also examine dynamic effects and optimization of assembly protocols. Finally, we discuss promising material functions and examples of their realization with nanocrystal assemblies.

Extinction and coexistence in a binary mixture of proliferating motile disks

A binary mixture of two-different-size proliferating motile disks is studied. As growth is space limited, we focus on the conditions such that there is a coexistence of both large and small disks, or dominance of the larger disks. The study involves systematically varying some system parameters, such as diffusivities, growth rates, and self-propulsion velocities. In particular, we demonstrate that diffusing faster confers a competitive advantage, so that larger disks can in the long time coexist or even dominate the smaller ones. In the case of self-propelled disks, a coexistence regime is induced by the activity where the two types of disks show the same spatial distribution: both particles are phase separated or both are homogeneously distributed in the whole system.

Theory and simulation of ligand functionalized nanoparticles - a pedagogical overview

Synthesizing reconfigurable nanoscale synthons with predictive control over shape, size, and interparticle interactions is a holy grail of bottom-up self-assembly. Grand challenges in their rational design, however, lie in both the large space of experimental synthetic parameters and proper understanding of the molecular mechanisms governing their formation. As such, computational and theoretical tools for predicting and modeling building block interactions have grown to become integral in modern day self-assembly research. In this review, we provide an in-depth discussion of the current state-of-the-art strategies available for modeling ligand functionalized nanoparticles. We focus on the critical role of how ligand interactions and surface distributions impact the emergent, pre-programmed behaviors between neighboring particles. To help build insights into the underlying physics, we first define an "ideal" limit - the short ligand, "hard" sphere approximation - and discuss all experimental handles through the lens of perturbations about this reference point. Finally, we identify theories that are capable of bridging interparticle interactions to nanoscale self-assembly and conclude by discussing exciting new directions for this field.

The entropy-controlled strategy in self-assembling systems

Self-assembly of various building blocks has been considered as a powerful approach to generate novel materials with tailorable structures and optimal properties. Understanding physicochemical interactions and mechanisms related to structural formation and transitions is of essential importance for this approach. Although it is well-known that diverse forces and energies can significantly contribute to the structures and properties of self-assembling systems, the potential entropic contribution remains less well understood. The past few years have witnessed rapid progress in addressing the entropic effects on the structures, responses, and functions in the self-assembling systems, and many breakthroughs have been achieved. This review provides a framework regarding the entropy-controlled strategy of self-assembly, through which the structures and properties can be tailored by effectively tuning the entropic contribution and its interplay with the enthalpic counterpart. First, we focus on the fundamentals of entropy in thermodynamics and the entropy types that can be explored for self-assembly. Second, we discuss the rules of entropy in regulating the structural organization in self-assembly and delineate the entropic force and superentropic effect. Third, we introduce the basic principles, significance and approaches of the entropy-controlled strategy in self-assembly. Finally, we present the applications where this strategy has been employed in fields like colloids, macromolecular systems and nonequilibrium assembly. This review concludes with a discussion on future directions and future research opportunities for developing and applying the entropy-controlled strategy in complex self-assembling systems.

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.

Microstructural diversity, nucleation paths, and phase behavior in binary mixtures of charged colloidal spheres

We study low-salt, binary aqueous suspensions of charged colloidal spheres of size ratio Γ = 0.57, number densities below the eutectic number density n_E, and number fractions of p = 1.00-0.40. The typical phase obtained by solidification from a homogeneous shear-melt is a substitutional alloy with a body centered cubic structure. In strictly gas-tight vials, the polycrystalline solid is stable against melting and further phase transformation for extended times. For comparison, we also prepare the same samples by slow, mechanically undisturbed deionization in commercial slit cells. These cells feature a complex but well reproducible sequence of global and local gradients in salt concentration, number density, and composition as induced by successive deionization, phoretic transport, and differential settling of the components, respectively. Moreover, they provide an extended bottom surface suitable for heterogeneous nucleation of the β-phase. We give a detailed qualitative characterization of the crystallization processes using imaging and optical microscopy. By contrast to the bulk samples, the initial alloy formation is not volume-filling, and we now observe also α- and β-phases with low solubility of the odd component. In addition to the initial homogeneous nucleation route, the interplay of gradients opens various further crystallization and transformation pathways leading to a great diversity of microstructures. Upon a subsequent increase in salt concentration, the crystals melt again. Wall-based, pebble-shaped β-phase crystals and facetted α-crystals melt last. Our observations suggest that the substitutional alloys formed in bulk experiments by homogeneous nucleation and subsequent growth are mechanically stable in the absence of solid-fluid interfaces but thermodynamically metastable.

Densest packings from size segregation of particles in geometric confinement

A correlation between density maximization and size segregation for packings of polydisperse particles in geometric confinement has been discovered, through the derivation of a general solution for the densest-packed zigzag arrangements of polydisperse particles. This solution is a size-graded structure in which the larger a particle the closer it is located to either end of the system, such that the smaller particles in the interior are encapsulated by the larger ones away from it. Any pair of different-sized adjacent particles is a grain-boundary-like configuration that reduces the overall packing efficiency of the system, and this solution corresponds to a minimization of excess-volume contributions from grain-boundary-like configurations of different-sized particles as a result of the clustering of equal- or like-sized particles. Our findings provide new insights into how structural order of polydisperse particles emerges in confined settings.

Tuning Stoichiometry to Promote Formation of Binary Colloidal Superlattices

The self-assembly of binary nanoparticle superlattices from colloidal mixtures is a promising method for the fabrication of complex colloidal cocrystal structures. However, binary mixtures often form amorphous or metastable phases instead of the thermodynamically stable phase. Here we show that in binary mixtures of differently sized spherical particles, an excess of the smaller component can promote-and, in some cases, may be necessary for-the self-assembly of a binary cocrystal. Using computer simulations, we identify two mechanisms responsible for this phenomenon. First, excess small particles act like plasticizers and enable systems to reach a greater supersaturation before kinetic arrest occurs. Second, they can disfavor competing structures that may interfere with the growth of the target structure. We find the phase behavior of simulated mixtures of nearly hard spheres closely matches published experimental results. We demonstrate the generality of our findings for mixtures of particles of arbitrary shape by presenting a binary mixture of hard shapes that only self-assembles with an excess of the smaller component.

Efficient event-driven simulations of hard spheres

 Hard spheres are arguably one of the most fundamental model systems in soft matter physics, and hence a common topic of simulation studies. Event-driven simulation methods provide an efficient method for studying the phase behavior and dynamics of hard spheres under a wide range of different conditions. Here, we examine the impact of several optimization strategies for speeding up event-driven molecular dynamics of hard spheres and present a light-weight simulation code that outperforms existing simulation codes over a large range of system sizes and packing fractions. The presented differences in simulation speed, typically a factor of five to ten, save significantly on both CPU time and energy consumption and may be a crucial factor for studying slow processes such as crystal nucleation and glassy dynamics.

Effect of polydispersity on the dynamics of active Brownian particles

We numerically study the dynamics and the phases of self-propelled disk-shaped particles of different sizes with soft repulsive potential in two dimensions. Size diversity is introduced by the polydispersity index (PDI) ε, which is the width of the uniform distribution of the particle's radius. The self-propulsion speed of the particles controls the activity v. We observe enhanced dynamics for large size diversity among the particles. We calculate the effective diffusion coefficient D_{eff} in the steady state. The system exhibits four distinct phases, jammed phase with small D_{eff} for small activity and liquid phase with enhanced D_{eff} for large activity. The number fluctuation is larger and smaller than the equilibrium limit in the liquid and jammed phases, respectively. Further, the jammed phase is of two types: solid jammed and liquid jammed for small and large PDI. Whereas the liquid phase is called motility induced phase separation (MIPS) liquid for small PDI and for large PDI, we find enhanced diffusivity and call it the pure liquid phase. The system is studied for three packing densities ϕ, and the response of the system for polydispersity is the same for all ϕ's. Our study can help understand the behavior of cells of various sizes in a tissue, artificial self-driven granular particles, or living organisms of different sizes in a dense environment.

Softness-driven complexity in supercrystals of gold nanoparticles

Many soft matter systems are composed of roughly spherical objects that can self-assemble in ordered structures. Unlike hard spheres, at high volume fraction these soft spheres adapt their shape to the local geometrical constraints and the question of space filling needs to be entirely revisited. Hydrophobically coated gold nanocrystals self-assemble in supercrystals and are good candidates to explore this question. When the soft coating is thin compared to the rigid core, a FCC structure is obtained, with a behaviour similar to that of hard spheres. In the opposite case, for a thick soft coating, a BCC structure is found instead. This paper focus on the intermediate region between these two classical structures. By varying the gold core radius R and the ligand fully extended length L, we establish a structure diagram based on a large experimental data set. The hexagonal Frank-Kasper C14 structure is observed for various values of R and L and can coexist with a FCC phase. Depending on the structure, values of the minimum thickness e of the ligand shell compared to L are different. These experimental results confirm that the C14 Frank-Kasper phase is a solution to the problem of filling the space with soft particles even with a rigid core and should help to establish pertinent models in order to predict the structures of the superlattices built by gold nanoparticles.

Monodisperse patchy particle glass former

Glass formers are characterized by their ability to avoid crystallization. As monodisperse systems tend to rapidly crystallize, the most common glass formers in simulations are systems composed of mixtures of particles with different sizes. Here, we make use of the ability of patchy particles to change their local structure to propose them as monodisperse glass formers. We explore monodisperse systems with two patch geometries: a 12-patch geometry that enhances the formation of icosahedral clusters and an 8-patch geometry that does not appear to strongly favor any particular local structure. We show that both geometries avoid crystallization and present glassy features at low temperatures. However, the 8-patch geometry better preserves the structure of a simple liquid at a wide range of temperatures and packing fractions, making it a good candidate for a monodisperse glass former.

Dynamical phase transitions and their relation to structural and thermodynamic aspects of glass physics

We review recent developments in structural-dynamical phase transitions in trajectory space based on dynamic facilitation theory. An open question is how the dynamic facilitation perspective on the glass transition may be reconciled with thermodynamic theories that posit collective reorganization accompanied by a growing static length scale and, eventually, a vanishing configurational entropy. In contrast, dynamic facilitation theory invokes a dynamical phase transition between an active phase (close to the normal liquid) and an inactive phase, which is glassy and whose order parameter is either a time-averaged dynamic or structural quantity. In particular, the dynamical phase transition in systems with non-trivial thermodynamics manifests signatures of a lower critical point that lies between the mode-coupling crossover and the putative Kauzmann temperature, at which a thermodynamic phase transition to an ideal glass state would occur. We review these findings and discuss such criticality in the context of the low-temperature decrease in configurational entropy predicted by thermodynamic theories of the glass transition.

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.

Tetrahedrality Dictates Dynamics in Hard Sphere Mixtures

The link between local structure and dynamical slowdown in glassy fluids has been the focus of intense debate for the better part of a century. Nonetheless, a simple method to predict the dynamical behavior of a fluid purely from its local structural features is still missing. Here, we demonstrate that the diffusivity of perhaps the most fundamental family of glass formers-hard sphere mixtures-can be accurately predicted based on just the packing fraction and a simple order parameter measuring the tetrahedrality of the local structure. Essentially, we show that the number of tetrahedral clusters in a hard sphere mixture is directly linked to its global diffusivity. Moreover, the same order parameter is capable of locally pinpointing particles in the system with high and low mobility. We attribute the power of the local tetrahedrality for predicting local and global dynamics to the high stability of tetrahedral clusters, the most fundamental building and densest-packing building blocks for a disordered fluid.

Tuning the Glass Transition: Enhanced Crystallization of the Laves Phases in Nearly Hard Spheres

Colloidal crystals with a diamond and pyrochlore structure display wide photonic band gaps at low refractive index contrasts. However, these low-coordinated and open structures are notoriously difficult to self-assemble from colloids interacting with simple pair interactions. To circumvent these problems, one can self-assemble both structures in a closely packed MgCu₂ Laves phase from a binary mixture of colloidal spheres and then selectively remove one of the sublattices. Although Laves phases have been proven to be stable in a binary hard-sphere system, they have never been observed to spontaneously crystallize in such a fluid mixture in simulations nor in experiments of micron-sized hard spheres due to slow dynamics. Here we demonstrate, using computer simulations, that softness in the interparticle potential suppresses the degree of 5-fold symmetry in the binary fluid phase and enhances crystallization of Laves phases in nearly hard spheres.

Grand canonical inverse design of multicomponent colloidal crystals

Inverse design methods are powerful computational approaches for creating colloidal systems which self-assemble into a target morphology by reverse engineering the Hamiltonian of the system. Despite this, these optimization procedures tend to yield Hamiltonians which are too complex to be experimentally realized. An alternative route to complex structures involves the use of several different components, however, conventional inverse design methods do not explicitly account for the possibility of phase separation into compositionally distinct structures. Here, we present an inverse design scheme for multicomponent colloidal systems by combining active learning with a method to directly compute their ground state phase diagrams. This explicitly accounts for phase separation and can locate stable regions of Hamiltonian parameter space which grid-based surveys are prone to miss. Using this we design low-density, binary structures with Lennard-Jones-like pairwise interactions that are simpler than in the single component case and potentially realizable in an experimental setting. This reinforces the concept that ground states of simple, multicomponent systems might be rich with previously unappreciated diversity, enabling the assembly of non-trivial structures with only few simple components instead of a single complex one.

Packing Polydisperse Colloids into Crystals: When Charge-Dispersity Matters

Monte Carlo simulations, fully constrained by experimental parameters, are found to agree well with a measured phase diagram of aqueous dispersions of nanoparticles with a moderate size polydispersity over a broad range of salt concentrations, c_{s}, and volume fractions, ϕ. Upon increasing ϕ, the colloids freeze first into coexisting compact solids then into a body centered cubic phase (bcc) before they melt into a glass forming liquid. The surprising stability of the bcc solid at high ϕ and c_{s} is explained by the interaction (charge) polydispersity and vibrational entropy.

Influence of Softness on the Stability of Binary Colloidal Crystals

Mixtures of two types of nanoparticles can self-assemble into a wide variety of binary colloidal crystals (also called binary nanoparticle superlattices), which are interesting for their structural diversity and potential applications. Although so-called packing models-which usually treat the particles as "hard" with only excluded volume interactions-seem to explain many reported dense crystalline phases, these models often fail to predict the right structure. Here, we examine the role of soft repulsive interparticle interactions on binary colloidal crystals comprising two sizes of spherical particles; such "softness" can arise due to ligand shells or screened electrostatics. We determine the ground state phase diagram of binary systems of particles interacting with an additive inverse power law potential using a basin hopping algorithm to calculate the enthalpy of an extremely large pool of candidate structures. We find that a surprisingly small amount of softness can destabilize dense packings in favor of less densely packed structures, which provides further evidence that considerations beyond packing are necessary for describing many of the observed phases of binary colloidal crystals. Importantly, we find that several of the phases stabilized by softness, which are characterized by relatively few interparticle contacts and a tendency for local icosahedral order, are more likely to be observed experimentally than those predicted by packing models. We also report a previously unknown dense AB₄ phase and conduct free energy calculations to examine how the stability of several crystals will vary with temperature. Our results further our understanding of why particular binary colloidal crystals form and will be useful as a reference for experimentalists working with softly repulsive colloids.

Precise modulation of molecular weight distribution for structural engineering

As one of the most critical molecular parameters, molecular weight distribution has a profound impact on the structure and properties of polymers. Quantitative and comprehensive understanding, however, has yet to be established, mainly due to the challenge in the precise control and regulation of molecular weight distribution. In this work, we demonstrated a robust and effective approach to artificially engineer the molecular weight distribution through precise recombination of discrete macromolecules. The width, symmetry, and other characteristics of the distribution can be independently manipulated to achieve absolute control, serving as a model platform for highlighting the importance of chain length heterogeneity in structural engineering. Different from their discrete counterparts, each individual component in dispersed samples experiences a varied degree of supercooling at a specific crystallization temperature. Non-uniform crystal nucleation and growth kinetics lead to distinct molecular arrangements. This work could bridge the gap between discrete and dispersed macromolecules, providing fundamental perspectives on the critical role of molecular weight distribution.

Engineering self-assembly of giant molecules in the condensed state based on molecular nanoparticles

In biological systems, it is well-known that the activities and functions of biomacromolecules are dictated not only by their primary chemistries, but also by their secondary, tertiary, and quaternary hierarchical structures. Achieving control of similar levels in synthetic macromolecules is yet to be demonstrated. Most of the critical molecular parameters associated with molecular and hierarchical structures, such as size, composition, topology, sequence, and stereochemistry, are heterogenous, which impedes the exploration and understanding of structure formation and manipulation. Alternatively, in the past few years we have developed a unique giant molecule system based on molecular nanoparticles, in which the above-mentioned molecular parameters, as well as interactions, are precisely defined and controlled. These molecules could self-assemble into a myriad of unconventional and unique structures in the bulk, thin films, and solution. Giant molecules thus offer a robust platform to manipulate the hierarchical structures via precise and modular assemblies of building blocks in an amplified size level compared with small molecules. It has been found that they are not only scientifically intriguing, but also technologically relevant.

Free Energy Landscape of Colloidal Clusters in Spherical Confinement

The structure of finite self-assembling systems depends sensitively on the number of constituent building blocks. Recently, it was demonstrated that hard sphere-like colloidal particles show a magic number effect when confined in emulsion droplets. Geometric construction rules permit a few dozen magic numbers that correspond to a discrete series of completely filled concentric icosahedral shells. Here, we investigate the free energy landscape of these colloidal clusters as a function of the number of their constituent building blocks for system sizes up to several thousand particles. We find that minima in the free energy landscape, arising from the presence of filled, concentric shells, are significantly broadened, compared to their atomic analogues. Colloidal clusters in spherical confinement can flexibly accommodate excess particles by ordering icosahedrally in the cluster center while changing the structure near the cluster surface. In between these magic number regions, the building blocks cannot arrange into filled shells. Instead, we observe that defects accumulate in a single wedge and therefore only affect a few tetrahedral grains of the cluster. We predict the existence of this wedge by simulation and confirm its presence in experiment using electron tomography. The introduction of the wedge minimizes the free energy penalty by confining defects to small regions within the cluster. In addition, the remaining ordered tetrahedral grains can relax internal strain by breaking icosahedral symmetry. Our findings demonstrate how multiple defect mechanisms collude to form the complex free energy landscape of colloidal clusters.