Freezing on a sphere

Diffusion, density, and defects on spheres

We simulate and model diffusion of spherical colloids of radius, a, on spherical surfaces of radius, R, as a function of relative size and surface concentration. Using Brownian dynamics simulations, we quantify diffusion and microstructure at different concentrations ranging from single particles to dense crystalline states. Self-diffusion and structural metrics (pair distribution, local density, and topological charge) are indistinguishable between spheres and planes for all concentrations up to dense liquid states. For concentrations approaching and greater than the freezing transition, smaller spheres with higher curvature show increased diffusivities and nonuniform density/topological defect distributions, which differ qualitatively from planar surfaces. The total topological charge varies quadratically with sphere radius for dense liquid states and linearly with sphere radius for dense crystals with icosahedrally organized grain scars. Between the dense liquid and dense crystal states on spherical surfaces is a regime of fluctuating and interacting defect clusters. We show local density governs self-diffusion in dense liquids on flat and spherical surfaces via the pair distribution. In contrast, dynamic topological defects couple to finite diffusivities through freezing and in low density crystal states on spherical surfaces, where neither exist on flat surfaces.

Self-limiting stacks of curvature-frustrated colloidal plates: Roles of intraparticle versus interparticle deformations

In geometrically frustrated assemblies local intersubunit misfits propagate to intra-assembly strain gradients, giving rise to anomalous self-limiting assembly thermodynamics. Here we use theory and coarse-grained simulation to study a recently developed class of "curvamer" particles, flexible shell-like particles that exhibit self-limiting assembly due to the build up of curvature deformation in cohesive stacks. To address a generic, yet poorly understood aspect of frustrated assembly, we introduce a model of curvamer assembly that incorporates both intraparticle shape deformation as well as compliance of interparticle cohesive gaps, an effect we can attribute to a finite range of attraction between particles. We show that the ratio of intraparticle (bending elasticity) to interparticle stiffness not only controls the regimes of self-limitation but also the nature of frustration propagation through curvamer stacks. We find a transition from uniformly bound, curvature-focusing stacks at small size to gap opened, uniformly curved stacks at large size is controlled by a dimensionless measure of inter- versus intracurvamer stiffness. The finite range of interparticle attraction determines the range of cohesion in stacks that are self-limiting, a prediction which is in strong agreement with numerical studies of our coarse-grained colloidal model. These predictions provide critical guidance for experimental realizations of frustrated particle systems designed to exhibit self-limitation at especially large multiparticle scales.

Nuclear size-regulated emergence of topological packing order on growing human lung alveolospheres

Within multicellular living systems, cells coordinate their positions with spatiotemporal accuracy to form various structures, setting the clock to control developmental processes and trigger maturation. These arrangements can be regulated by tissue topology, biochemical cues, as well as mechanical perturbations. However, the fundamental rules of how local cell packing order is regulated in forming three-dimensional (3D) multicellular architectures remain unclear. Furthermore, how cellular coordination evolves during developmental processes, and whether this cell patterning behavior is indicative of more complex biological functions, is largely unknown. Here, using human lung alveolospheres as a model system, by combining experiments and numerical simulations, we find that, surprisingly, cell packing behavior on alveolospheres resembles hard-disk packing but with increased randomness; the stiffer cell nuclei act as the hard disks surrounded by deformable cell bodies. Interestingly, we observe the emergence of topological packing order during alveolosphere growth, as a result of increasing nucleus-to-cell size ratio. Specifically, we find more hexagon-concentrated cellular packing with increasing bond orientational order, indicating a topological gas-to-liquid transition. Additionally, by osmotically changing the compactness of cells on alveolospheres, we observe that the variations in packing order align with the change of nucleus-to-cell size ratio. Together, our findings reveal the underlying rules of cell coordination and topological phases during human lung alveolosphere growth. These static packing characteristics are consistent with cell dynamics, together suggesting that better cellular packing stabilizes local cell neighborhoods and may regulate more complex biological functions such as organ development and cellular maturation.

Programmable Potentials Choreograph Defects in a Colloidal Crystal Shell

Crystallization on spherical surfaces is obliged by topology to induce lattice defects. But controlling the organization of such defects remains a great challenge due to the long-range constraints of the curved geometry. Here, we report on DNA-coated colloids whose programmable interaction potentials can be used to regulate the arrangement of defects and even achieve perfect icosahedral order on a sphere. Combined simulations and theoretical analysis show how the potential can be tuned by changing the temperature, thereby controlling the number of defects. An explicit expression for the effective potential is derived, allowing us to distinguish the effects of entropic repulsion and enthalpic attraction. Altogether, the present findings provide insights into the physics of crystallization on curved spaces and may be used for designing desired crystal geometries.

Thermal stability and secondary aggregation of self-limiting, geometrically frustrated assemblies: Chain assembly of incommensurate polybricks

In geometrically frustrated assemblies, equilibrium self-limitation manifests in the form of a minimum in the free energy per subunit at a finite, multisubunit size which results from the competition between the elastic costs of frustration within an assembly and the surface energy at its boundaries. Physical realizations-from ill-fitting particle assemblies to self-twisting protein superstructures-are capable of multiple mechanisms of escaping the cumulative costs of frustration, resulting in unlimited equilibrium assembly, including elastic modes of "shape flattening" and the formation of weak, defective bonds that screen intra-assembly stresses. Here we study a model of one-dimensional chain assembly of incommensurate "polybricks" and determine its equilibrium assembly as a function of temperature, concentration, degree of shape frustration, elasticity, and interparticle binding, notably focusing on how weakly cohesive, defective bonds give rise to strongly temperature-dependent assembly. Complex assembly behavior derives from the competition between multiple distinct local minima in the free-energy landscape, including self-limiting chains, weakly bound aggregates of self-limiting chains, and strongly bound, elastically defrustrated assemblies. We show that this scenario, in general, gives rise to anomalous multiple aggregation behavior, in which disperse subunits (stable at low concentration and high temperature) first exhibit a primary aggregation transition to self-limiting chains (at intermediate concentration and temperature) which are ultimately unstable to condensation into unlimited assembly of finite-chains through weak binding beyond a secondary aggregation transition (at low temperature and high concentration). We show that window of stable self-limitation is determined both by the elastic costs of frustration in the assembly as well as energetic and entropic features of intersubunit binding.

Geometric frustration of hard-disk packings on cones

Conical surfaces pose an interesting challenge to crystal growth: A crystal growing on a cone can wrap around and meet itself at different radii. We use a disk-packing algorithm to investigate how this closure constraint can geometrically frustrate the growth of single crystals on cones with small opening angles. By varying the crystal seed orientation and cone angle, we find that-except at special commensurate cone angles-crystals typically form a seam that runs along the axial direction of the cone, while near the tip, a disordered particle packing forms. We show that the onset of disorder results from a finite-size effect that depends strongly on the circumference and not on the seed orientation or cone angle. This finite-size effect occurs also on cylinders, and we present evidence that on both cylinders and cones, the defect density increases exponentially as circumference decreases. We introduce a simple model for particle attachment at the seam that explains the dependence on the circumference. Our findings suggest that the growth of single crystals can become frustrated even very far from the tip when the cone has a small opening angle. These results may provide insights into the observed geometry of conical crystals in biological and materials applications.

Early-stage bifurcation of crystallization in a sphere

Bifurcations in kinetic pathways decide the evolution of a system. An example is crystallization, in which the thermodynamically stable polymorph may not form due to kinetic hindrance. Here, we use confined self-assembly to investigate the interplay of thermodynamics and kinetics in the crystallization pathways of finite clusters. We report the observation of decahedral clusters from colloidal particles in emulsion droplets and show that these decahedral clusters can be thermodynamically stable, just like icosahedral clusters. Our hard sphere simulations reveal how the development of the early nucleus shape passes through a bifurcation that decides the cluster symmetry. A geometric argument explains why decahedral clusters are kinetically hindered and why icosahedral clusters can be dominant even if they are not in the thermodynamic ground state.

Visualizing defect dynamics by assembling the colloidal graphene lattice

Graphene has been under intense scientific interest because of its remarkable optical, mechanical and electronic properties. Its honeycomb structure makes it an archetypical two-dimensional material exhibiting a photonic and phononic band gap with topologically protected states. Here, we assemble colloidal graphene, the analogue of atomic graphene using pseudo-trivalent patchy particles, allowing particle-scale insight into crystal growth and defect dynamics. We directly observe the formation and healing of common defects, like grain boundaries and vacancies using confocal microscopy. We identify a pentagonal defect motif that is kinetically favoured in the early stages of growth, and acts as seed for more extended defects in the later stages. We determine the conformational energy of the crystal from the bond saturation and bond angle distortions, and follow its evolution through the energy landscape upon defect rearrangement and healing. These direct observations reveal that the origins of the most common defects lie in the early stages of graphene assembly, where pentagons are kinetically favoured over the equilibrium hexagons of the honeycomb lattice, subsequently stabilized during further growth. Our results open the door to the assembly of complex 2D colloidal materials and investigation of their dynamical, mechanical and optical properties.

Packing and emergence of the ordering of rods in a spherical monolayer

Spatially ordered systems confined to surfaces such as spheres exhibit interesting topological structures because of curvature induced frustration in orientational and translational order. The study of these structures is important for investigating the interplay between the geometry, topology, and elasticity, and for their potential applications in materials science, such as engineering directionally binding particles. In this work, we numerically simulate a spherical monolayer of soft repulsive spherocylinders (SRSs) and study the packing of rods and their ordering transition as a function of the packing fraction. In the model that we study, the centers of mass of the spherocylinders (situated at their geometric centers) are constrained to move on a spherical surface. The spherocylinders are free to rotate about any axis that passes through their respective centers of mass. We show that, up to moderate packing fractions, a two dimensional liquid crystalline phase is formed whose orientational ordering increases continuously with increasing density. This monolayer of orientationally ordered SRS particles at medium densities resembles a hedgehog-long axes of the SRS particles are aligned along the local normal to the sphere. At higher packing fractions, the system undergoes a transition to the solid phase, which is riddled with topological point defects (disclinations) and grain boundaries that divide the whole surface into several domains.

Curvature-controlled geometrical lensing behavior in self-propelled colloidal particle systems

In many biological systems, the curvature of the surfaces cells live on influences their collective properties. Curvature should likewise influence the behavior of active colloidal particles. We show using molecular simulation of self-propelled active particles on surfaces of Gaussian curvature (both positive and negative) how curvature sign and magnitude can alter the system's collective behavior. Curvature acts as a geometrical lens and shifts the critical density of motility-induced phase separation (MIPS) to lower values for positive curvature and higher values for negative curvature, which we explain theoretically by the nature of parallel lines in spherical and hyperbolic space. Curvature also fluidizes dense MIPS clusters due to the emergence of defect patterns disrupting the crystalline order inside the clusters. Using our findings, we engineer three confining surfaces that strategically combine regions of different curvature to produce a host of novel dynamical behaviors, including cyclic MIPS on spherocylinders, directionally biased cyclic MIPS on spherocones, and position dependent cluster fluctuations on metaballs.

Self-assembly of emulsion droplets through programmable folding

In the realm of particle self-assembly, it is possible to reliably construct nearly arbitrary structures if all the pieces are distinct^1-3, but systems with fewer flavours of building blocks have so far been limited to the assembly of exotic crystals^4-6. Here we introduce a minimal model system of colloidal droplet chains⁷, with programmable DNA interactions that guide their downhill folding into specific geometries. Droplets are observed in real space and time, unravelling the rules of folding. Combining experiments, simulations and theory, we show that controlling the order in which interactions are switched on directs folding into unique structures, which we call colloidal foldamers⁸. The simplest alternating sequences (ABAB...) of up to 13 droplets yield 11 foldamers in two dimensions and one in three dimensions. Optimizing the droplet sequence and adding an extra flavour uniquely encodes more than half of the 619 possible two-dimensional geometries. Foldamers consisting of at least 13 droplets exhibit open structures with holes, offering porous design. Numerical simulations show that foldamers can further interact to make complex supracolloidal architectures, such as dimers, ribbons and mosaics. Our results are independent of the dynamics and therefore apply to polymeric materials with hierarchical interactions on all length scales, from organic molecules all the way to Rubik's Snakes. This toolbox enables the encoding of large-scale design into sequences of short polymers, placing folding at the forefront of materials self-assembly.

Competition between clustering and phase separation in binary mixtures containing SALR particles

We investigate by Monte Carlo simulations a mixture of particles with competing interactions (hard-sphere two-Yukawa, HSTY) and hard spheres (HS), with same diameters σ and a square-well (SW) cross attraction. In a recent study [G. Munaò et al., J. Phys. Chem. B, 2022, 126, 2027-2039], we have analysed situations-in terms of relative concentration and attraction strength-where HS promote the formation of clusters involving particles of both species under thermodynamic conditions that would not allow for clustering of the pure HSTY fluid. Here, we focus on the role played by the range of cross attraction in determining the equilibrium structure of the mixture, starting from a homogeneous low-density state. When the width of the well exceeds approximately σ, clustering takes place in the system, with aggregates characterised by various sizes and shapes. Only for low HSTY concentrations (less than 10%) a single big cluster appears, anticipating the behaviour observed for a wider well, around 1.2σ. In the latter case, a spherical cluster encompassing almost all particles is the stable structure at equilibrium. We interpret this outcome as a macrophase, liquid-vapour separation where the spherical cluster is just the form taken at low density by the liquid phase inside the vapour phase: indeed, when the density takes larger values, periodic boundary conditions select liquid-vapour interfaces with other non-spherical shapes, similarly as found for a finite sample of simple fluid going through the liquid-vapour coexistence region. For still higher densities we document the existence of a solid phase characterized by the alternation of bilayers filled with particles of one species and bilayers of the other species, giving the solid a peculiar wafer structure.

Entropic control of nanoparticle self-assembly through confinement

Entropy can be the sole driving force for the construction and regulation of ordered structures of soft matter systems. Specifically, under confinement, the entropic penalty could induce enhanced entropic effects which potentially generate visually ordered structures. Therefore, spatial confinement or a crowding environment offers an important approach to control entropy effects in these systems. Here, we review how spatial confinement-mediated entropic effects accurately and even dynamically control the self-assembly of nanoscale objects into ordered structures, focusing on our efforts towards computer simulations and theoretical analysis. First, we introduce the basic principle of entropic ordering through confinement. We then introduce the applications of this concept to various systems containing nanoparticles, including polymer nanocomposites, biological macromolecular systems and macromolecular colloids. Finally, the future directions and challenges for tailoring nanoparticle organization through spatial confinement-mediated entropic effects are detailed. We expect that this review could stimulate further efforts in the fundamental research on the relationship between confinement and entropy and in the applications of this concept for designer nanomaterials.

Observation of two-step melting on a sphere

Melting in two-dimensional flat space is typically two-step and via the hexatic phase. How melting proceeds on a curved surface, however, is not known. Topology mandates that crystalline particle assemblies on these surfaces harbor a finite density of defects, which itself can be ordered, like the icosahedral ordering of 5-coordinated disclination defects on a sphere. Thus, melting even on a sphere, the simplest closed surface, involves the loss of both crystalline and defect order. Probing the interplay of these two forms of order, however, requires a system in which melting can be performed in situ, and this has not been achieved hitherto. Here, by tuning interparticle interactions in situ, we report an observation of an intermediate hexatic phase during the melting of colloidal crystals on a sphere. Remarkably, we observed a precipitous drop in icosahedral defect order in the hexatic phase where the shear modulus is expected to vanish. Furthermore, unlike in flat space, where disorder can fundamentally alter the nature of the melting process, on the sphere, we observed the signature characteristics of ideal melting. Our findings have profound implications for understanding, for instance, the self-assembly and maturation dynamics of viral capsids and also phase transitions on curved surfaces.

Interfacial self-assembly of SiO2-PNIPAM core-shell particles with varied crosslinking density

Spherical particles confined to liquid interfaces generally self-assemble into hexagonal patterns. It was theoretically predicted by Jagla two decades ago that such particles interacting via a soft repulsive potential are able to form complex, anisotropic assembly phases. Depending on the shape and range of the potential, the predicted minimum energy configurations include chains, rhomboid and square phases. We recently demonstrated that deformable core-shell particles consisting of a hard silica core and a soft poly(N-isopropylacrylamide) shell adsorbed at an air/water interface can form chain phases if the crosslinker is primarily incorporated around the silica core. Here, we systematically investigate the interfacial self-assembly behavior of such SiO₂-PNIPAM core-shell particles as a function of crosslinker content and core size. We observe chain networks predominantly at low crosslinking densities and smaller core sizes, whereas higher crosslinking densities lead to the formation of rhomboid packing. We correlate these results with the interfacial morphologies of the different particle systems, where the ability to expand at the interface and form a thin corona at the periphery depends on the degree of crosslinking close to the core. We perform minimum energy calculations based on Jagla-type pair potentials with different shapes of the soft repulsive shoulder. We compare the theoretical phase diagram with experimental findings to infer to which extent the interfacial interactions of the experimental system may be captured by Jagla pair-wise interaction potentials.

Robotic swimming in curved space via geometric phase

Locomotion by shape changes or gas expulsion is assumed to require environmental interaction, due to conservation of momentum. However, as first noted in [J. Wisdom, Science 299, 1865-1869 (2003)] and later in [E. Guéron, Sci. Am. 301, 38-45 (2009)] and [J. Avron, O. Kenneth, New J. Phys, 8, 68 (2006)], the noncommutativity of translations permits translation without momentum exchange in either gravitationally curved spacetime or the curved surfaces encountered by locomotors in real-world environments. To realize this idea which remained unvalidated in experiments for almost 20 y, we show that a precision robophysical apparatus consisting of motors driven on curved tracks (and thereby confined to a spherical surface without a solid substrate) can self-propel without environmental momentum exchange. It produces shape changes comparable to the environment's inverse curvatures and generates movement of [Formula: see text] cm per gait. While this simple geometric effect predominates over short time, eventually the dissipative (frictional) and conservative forces, ubiquitous in real systems, couple to it to generate an emergent dynamics in which the swimming motion produces a force that is counter-balanced against residual gravitational forces. In this way, the robot both swims forward without momentum and becomes fixed in place with a finite momentum that can be released by ceasing the swimming motion. We envision that our work will be of use in a broad variety of contexts, such as active matter in curved space and robots navigating real-world environments with curved surfaces.

Coloration in Supraparticles Assembled from Polyhedral Metal-Organic Framework Particles

Supraparticles are spherical colloidal crystals prepared by confined self‐assembly processes. A particularly appealing property of these microscale structures is the structural color arising from interference of light with their building blocks. Here, we assemble supraparticles with high structural order that exhibit coloration from uniform, polyhedral metal–organic framework (MOF) particles. We analyse the structural coloration as a function of the size of these anisotropic building blocks and their internal structure. We attribute the angle‐dependent coloration of the MOF supraparticles to the presence of ordered, onion‐like layers at the outermost regions. Surprisingly, even though different shapes of the MOF particles have different propensities to form these onion layers, all supraparticle dispersions show well‐visible macroscopic coloration, indicating that local ordering is sufficient to generate interference effects.

Data-driven criterion for the solid-liquid transition of two-dimensional self-propelled colloidal particles far from equilibrium

We establish an explicit data-driven criterion for identifying the solid-liquid transition of two-dimensional self-propelled colloidal particles in the far from equilibrium parameter regime, where the transition points predicted by different conventional empirical criteria for melting and freezing diverge. This is achieved by applying a hybrid machine learning approach that combines unsupervised learning with supervised learning to analyze a huge amount of the system's configurations in the nonequilibrium parameter regime on an equal footing. Furthermore, we establish a generic data-driven evaluation function, according to which the performance of different empirical criteria can be systematically evaluated and improved. In particular, by applying this evaluation function, we identify a new nonequilibrium threshold value for the long-time diffusion coefficient, based on which the predictions of the corresponding empirical criterion are greatly improved in the far from equilibrium parameter regime. These data-driven approaches provide a generic tool for investigating phase transitions in complex systems where conventional empirical ones face difficulties.

Entropy-Driven Unconventional Crystallization of Spherical Colloidal Nanocrystals Confined in Wide Cylinders

The densest packings of identical spherical colloidal nanocrystals in a thin cylinder generally give rise to confinement-induced chiral ordering. Here, we demonstrate that entropy can invalidate Pauling's packing rules for the nanocrystals confined in wide cylinders and novel ordered phases, where chiral ordering is broken, emerge. The nucleation and growth of spherical colloidal nanocrystals in the wide cylinders exhibit unique mechanisms which are distinctly different from that of thin ones. Furthermore, theoretical models which capture the essential physics of the ordering transitions are developed to reproduce the achiral ordering and reveal that the ordered phases are thermodynamically stable and stabilized through confinement-mediated entropic effect. These findings demonstrate that entropy arising from thermal motion can invalidate Pauling's packing rules of spherical colloidal nanocrystals confined in cylinders, which provides new insights into confinement physics of colloidal particles and might inspire nonintuitive design rules for the fabrication of novel ordered phases through confinement.

Curved colloidal crystals of discoids at near-critical liquid-liquid interface

The assembly of disc-shaped particles at curved liquid-liquid interfaces was studied by using confocal microscopy. The interface is formed by a phase-separating critical liquid mixture of 2,6-lutidine and heavy water, where the colloids spontaneously assembled forming a dome. The novelty of this system is three-fold. First, the domes can be constructed and annihilated remotely and reversibly, which allows dynamic control of the colloidal assembly. Second, the effect of curvature can be investigated by analyzing domes of different radii ranging from 5 μm to 125 μm. Third, the slow dynamics due to hydrodynamic interaction among the particles can be utilized to investigate the time-evolution of defect morphology. Unlike the widely studied repulsive colloids, the interparticle potential near the critical point has an attractive component. I contrasted the packing and defects morphology of a solid-like and liquid-like dome differing in particle number density. In the solid-like dome, a chain of 5- and 7-fold coordinated particles was observed. The analysis of trajectories showed that particles were bound in a potential well of a depth of about ten times the thermal energy, which matched well with the calculation of the pair-potential by considering the attractive critical Casimir force among the particles. In the liquid-like dome, 6-fold particles separated by clusters of 5- and 7-coordinated particles were observed, which is suggestive of liquid-solid coexistence. The uniqueness of this system will open up a new research avenue to investigate the effect of varying curvature on the crystallization, defects, and phase diagram of colloidal assemblies.

Elongation and percolation of defect motifs in anisotropic packing problems

We examine the regime between crystalline and amorphous packings of anisotropic objects on surfaces of different genus by continuously varying their size distribution or shape from monodispersed spheres to bidispersed mixtures or monodispersed ellipsoidal particles; we also consider an anisotropic variant of the Thomson problem with a mixture of charges. With increasing anisotropy, we first observe the disruption of translational order with an intermediate orientationally ordered hexatic phase as proposed by Nelson, Rubinstein and Spaepen, and then a transition to amorphous state. By analyzing the structure of the disclination motifs induced, we show that the hexatic-amorphous transition is caused by the growth and connection of disclination grain boundaries, suggesting this transition lies in the percolation universality class in the scenarios considered.

Spontaneous organization of supracolloids into three-dimensional structured materials

Periodic nano- or microscale structures are used to control light, energy and mass transportation. Colloidal organization is the most versatile method used to control nano- and microscale order, and employs either the enthalpy-driven self-assembly of particles at a low concentration or the entropy-driven packing of particles at a high concentration. Nonetheless, it cannot yet provide the spontaneous three-dimensional organization of multicomponent particles at a high concentration. Here we combined these two concepts into a single strategy to achieve hierarchical multicomponent materials. We tuned the electrostatic attraction between polymer and silica nanoparticles to create dynamic supracolloids whose components, on drying, reorganize by entropy into three-dimensional structured materials. Cryogenic electron tomography reveals the kinetic pathways, whereas Monte Carlo simulations combined with a kinetic model provide design rules to form the supracolloids and control the kinetic pathways. This approach may be useful to fabricate hierarchical hybrid materials for distinct technological applications.

Ultracold Bosons on a Regular Spherical Mesh

Here, the zero-temperature phase behavior of bosonic particles living on the nodes of a regular spherical mesh ("Platonic mesh") and interacting through an extended Bose-Hubbard Hamiltonian has been studied. Only the hard-core version of the model for two instances of Platonic mesh is considered here. Using the mean-field decoupling approximation, it is shown that the system may exist in various ground states, which can be regarded as analogs of gas, solid, supersolid, and superfluid. For one mesh, by comparing the theoretical results with the outcome of numerical diagonalization, I manage to uncover the signatures of diagonal and off-diagonal spatial orders in a finite quantum system.

Cooperatively rearranging regions change shape near the mode-coupling crossover for colloidal liquids on a sphere

The structure and dynamics of liquids on curved surfaces are often studied through the lens of frustration-based approaches to the glass transition. Competing glass transition theories, however, remain largely untested on such surfaces and moreover, studies hitherto have been entirely theoretical/numerical. Here we carry out single particle-resolved imaging of dynamics of bi-disperse colloidal liquids confined to the surface of a sphere. We find that mode-coupling theory well captures the slowing down of dynamics in the moderate to deeply supercooled regime. Strikingly, the morphology of cooperatively rearranging regions changed from string-like to compact near the mode-coupling crossover—a prediction unique to the random first-order theory of glasses. Further, we find that in the limit of strong curvature, Mermin–Wagner long-wavelength fluctuations are irrelevant and liquids on a sphere behave like three-dimensional liquids. A comparative evaluation of competing mechanisms is thus an essential step towards uncovering the true nature of the glass transition.

The static and dynamic behavior of condensed phases residing on curved surfaces can be fundamentally different from their counterparts in Euclidean space. Singh et al. test several competing glass theories on colloidal liquids confined to the surface of a sphere and show they behave like 3D bulk liquids.

Colloidal assembly in droplets: structures and optical properties

Colloidal assembly in emulsion drops provides fundamental tools for studying optimum particle arrangement under spherical confinement and practical means for producing photonic microparticles. Recent progress has revealed that energetically favored cluster configurations are different from conventional supraballs, which could enhance optical performance. This paper reviews state-of-the-art emulsion-templated colloidal clusters, and particularly focuses on recently reported novel structures such as icosahedral, decahedral, and single-crystalline face-centered cubic (fcc) clusters. We classify the clusters according to the number of component particles as small (N < O(10²)), medium (O(10²) ≤N≤O(10⁴)), and large (N≥O(10⁵)). For each size of clusters, we discuss the detailed structures, mechanisms of cluster formation, and optical properties and potential applications. Finally, we outline current challenges and questions that require further investigation.

Formation and stability of conformal spirals in confined 2D crystals

We investigate the ground-state and dynamical properties of nonuniform two-dimensional (2D) clusters of long-range interacting particles. We demonstrate that, when the confining external potential is designed to produce an approximate 1/r²density profile, the particles crystallize into highly ordered structures featuring spiral crystalline lines. Despite the strong inhomogeneity of the observed configurations, most of them are characterized by small density of topological defects, typical of conformal crystals, and the net topological charge induced by the simply-connected geometry of the system is concentrated near the cluster center. These crystals are shown to be robust with respect to thermal fluctuations up to a certain threshold temperature, above which the net charge is progressively redistributed from the center to the rest of the system and the topological order is lost. The crystals are also resilient to the shear stress produced by a small nonuniform azimuthal force field, rotating as a rigid body (RB). For larger forces, topological defects proliferate and the RB rotation gives place to plastic flow.

Topological defects of dipole patchy particles on a spherical surface

We investigate the assembly of dipole-like patchy particles confined to a spherical surface by Brownian dynamics simulations. The surface property of the spherical particle is described by the spherical harmonic Y10, and the orientation of the particle is defined as the uniaxial axis. On a flat space, we observe a defect-free square lattice with nematic order. On a spherical surface, defects appear due to the topological constraint. As for the director field, four defects of winding number +1/2 are observed, satisfying the Euler characteristic. We have found many configurations of the four defects lying near a great circle. Regarding the positional order for the square lattice, eight grain boundary scars proliferate linearly with the sphere size. The positions and orientations of the eight grain boundary scars are strongly related to the four +1/2 defect cores.

Stress driven fractionalization of vacancies in regular packings of elastic particles

Elucidating the interplay of defects and stress at the microscopic level is a fundamental physical problem that has a strong connection with materials science. Here, based on the two-dimensional crystal model, we show that the instability mode of vacancies with varying size and morphology conforms to a common scenario. A vacancy under compression is fissioned into a pair of dislocations that glide and vanish at the boundary. This neat process is triggered by the local shear stress around the vacancy. The remarkable fractionalization of vacancies creates rich modes of interaction between vacancies and other topological defects, and provides a new dimension for mechanical engineering of defects in extensive crystalline structures.

Stress-induced ordering of two-dimensional packings of elastic spheres

Packing of particles in confined environments is a common problem in multiple fields. Here, based on the two-dimensional Hertzian particle model, we study the packing of deformable spherical particles under compression and reveal the crucial role of stress as an ordering field in regulating particle arrangement. Specifically, under increasing compression, the squeezed particles spontaneously organize into regular packings in the sequence of triangular and square lattices, pentagonal tessellation, and the reentrant triangular lattice. The rich ordered patterns and complex structures revealed in this work suggest a fruitful organizational strategy based on the interplay of external stress and intrinsic elastic instability of particle arrays.

Dislocation screening in crystals with spherical topology

Whereas disclination defects are energetically prohibitive in two-dimensional flat crystals, their existence is necessary in crystals with spherical topology, such as viral capsids, colloidosomes, or fullerenes. Such a geometrical frustration gives rise to large elastic stresses, which render the crystal unstable when its size is significantly larger than the typical lattice spacing. Depending on the compliance of the crystal with respect to stretching and bending deformations, these stresses are alleviated either by a local increase of the intrinsic curvature in proximity of the disclinations or by the proliferation of excess dislocations, often organized in the form of one-dimensional chains known as "scars." The associated strain field of the scars is such as to counterbalance the one resulting from the isolated disclinations. Here we develop a continuum theory of dislocation screening in two-dimensional closed crystals with genus one. Upon modeling the flux of scars emanating from a given disclination as an independent scalar field, we demonstrate that the elastic energy of closed two-dimensional crystals with various degrees of asphericity can be expressed as a simple quadratic function of the screened topological charge of the disclinations, at both zero and finite temperature. This allows us to predict the optimal density of the excess dislocations as well as the minimal stretching energy attained by the crystal.

Properties of surface Landau-de Gennes Q-tensor models

Uniaxial nematic liquid crystals whose molecular orientation is subjected to tangential anchoring on a curved surface offer a non trivial interplay between the geometry and the topology of the surface and the orientational degree of freedom. We consider a general thin film limit of a Landau-de Gennes Q-tensor model which retains the characteristics of the 3D model. From this, previously proposed surface models follow as special cases. We compare fundamental properties, such as the alignment of the orientational degrees of freedom with principle curvature lines, order parameter symmetry and phase transition type for these models, and suggest experiments to identify suitable model assumptions.

Patchy particles by self-assembly of star copolymers on a spherical substrate: Thomson solutions in a geometric problem with a color constraint

Confinement or geometric frustration is known to alter the structure of soft matter, including copolymeric melts, and can consequently be used to tune structure and properties. Here we investigate the self-assembly of ABC and ABB 3-miktoarm star copolymers confined to a spherical shell using coarse-grained dissipative particle dynamics simulations. In bulk and flat geometries the ABC stars form hexagonal tilings, but this is topologically prohibited in a spherical geometry which normally is alleviated by forming pentagonal tiles. However, the molecular architecture of the ABC stars implies an additional 'color constraint' which only allows even tilings (where all polygons have an even number of edges) and we study the effect of these simultaneous constraints. We find that both ABC and ABB systems form spherical tiling patterns, the type of which depends on the radius of the spherical substrate. For small spherical substrates, all solutions correspond to patterns solving the Thomson problem of placing mobile repulsive electric charges on a sphere. In ABC systems we find three coexisting, possibly different tilings, one in each color, each of them solving the Thomson problem simultaneously. For all except the smallest substrates, we find competing solutions with seemingly degenerate free energies that occur with different probabilities. Statistically, an observer who is blind to the differences between B and C can tell from the structure of the A domains if the system is an ABC or an ABB star copolymer system.

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.

Precise Self-Positioning of Colloidal Particles on Liquid Emulsion Droplets

Decorating emulsion droplets by particles stabilizes foodstuff and pharmaceuticals. Interfacial particles also influence aerosol formation, thus impacting atmospheric CO₂ exchange. While studies of particles at disordered droplet interfaces abound in the literature, such studies for ubiquitous ordered interfaces are not available. Here, we report such an experimental study, showing that particles residing at crystalline interfaces of liquid droplets spontaneously self-position to specific surface locations, identified as structural topological defects in the crystalline surface monolayer. This monolayer forms at temperature T = T_s, leaving the droplet liquid and driving at T_d < T_s a spontaneous shape-change transition of the droplet from spherical to icosahedral. The particle's surface position remains unchanged in the transition, demonstrating these positions to coincide with the vertices of the sphere-inscribed icosahedron. Upon further cooling, droplet shape-changes to other polyhedra occur, with the particles remaining invariably at the polyhedra's vertices. At still lower temperatures, the particles are spontaneously expelled from the droplets. Our results probe the molecular-scale elasticity of quasi-two-dimensional curved crystals, impacting also other fields, such as protein positioning on cell membranes, controlling essential biological functions. Using ligand-decorated particles, and the precise temperature-tunable surface position control found here, may also allow using these droplets for directed supra-droplet self-assembly into larger structures, with a possible post-assembly structure fixation by UV polymerization of the droplet's liquid.

Interfacial aggregation of Janus rods in binary polymer blends and their effect on phase separation

Janus particles interfacially self-assemble into different structures when incorporated into multiphase systems. Dissipative particle dynamics simulations are employed herein to investigate the interplay between aggregation mechanisms and phase separation in polymer blends. Shorter rods with a standing configuration become increasingly "caged" or trapped in larger aggregates as weight fraction increases, which is reflected in the way that their diffusion is coupled to their aggregation rates. Janus rods of higher aspect ratios that are tilted at the interface aggregate side-by-side and are able to hinder phase separation kinetics. This is due to a combination of individual Janus rod conformations at the interface, their intrinsic aggregation mechanisms, aggregate fractal dimension, and aggregation rates, and can also be traced back to the scaling of the diffusion coefficient of aggregates with their size. Findings presented provide insight into the mechanisms governing two dimensionally growing colloidal aggregates at fluid interfaces, more specifically, those associated with Janus particles, and shed light on the potential of these systems in paving the way for designing new functional materials.

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.

Intermittency of Deformation and the Elastic Limit of an Icosahedral Virus under Compression

Viruses undergo mesoscopic morphological changes as they interact with host interfaces and in response to chemical cues. The dynamics of these changes, over the entire temporal range relevant to virus processes, are unclear. Here, we report on creep compliance experiments on a small icosahedral virus under uniaxial constant stress. We find that even at small stresses, well below the yielding point and generally thought to induce a Hookean response, strain continues to develop in time via sparse, randomly distributed, relatively rapid plastic events. The intermittent character of mechanical compliance only appears above a loading threshold, similar to situations encountered in granular flows and the plastic deformation of crystalline solids. The threshold load is much smaller for the empty capsids of the brome mosaic virus than for the wild-type virions. The difference highlights the involvement of RNA in stabilizing the assembly interface. Numerical simulations of spherical crystal deformation suggest intermittency is mediated by lattice defect dynamics and identify the type of compression-induced defect that nucleates the transition to plasticity.

Growth of curved crystals: competition between topological defect nucleation and boundary branching

Topological defect nucleation and boundary branching in crystal growth on a curved surface are two typical elastic instabilities driven by curvature induced stress, and have usually been discussed separately in the past. In this work they are simultaneously considered during crystal growth on a sphere. Phase diagrams with respect to sphere radius, size, edge energy and stiffness of the crystal for the equilibrium crystal morphologies are achieved by theoretical analysis and validated by Brownian dynamics simulations. The simulation results further demonstrate the detail of morphological evolution governed by these two different stress relaxation modes. Topological defect nucleation and boundary branching not only compete with each other but also coexist in a range of combinations of factors. Clarification of the interaction mechanism provides a better understanding of various curved crystal morphologies for their potential applications.

Building Reconfigurable Devices Using Complex Liquid-Fluid Interfaces

Liquid-fluid interfaces provide a platform both for structuring liquids into complex shapes and assembling dimensionally confined, functional nanomaterials. Historically, attention in this area has focused on simple emulsions and foams, in which surface-active materials such as surfactants or colloids stabilize structures against coalescence and alter the mechanical properties of the interface. In recent decades, however, a growing body of work has begun to demonstrate the full potential of the assembly of nanomaterials at liquid-fluid interfaces to generate functionally advanced, biomimetic systems. Here, a broad overview is given, from fundamentals to applications, of the use of liquid-fluid interfaces to generate complex, all-liquid devices with a myriad of potential applications.

Spherical structure factor and classification of hyperuniformity on the sphere

Understanding how particles are arranged on the surface of a sphere is not only central to numerous physical, biological, soft matter, and materials systems but also finds applications in computational problems, approximation theory, and analysis of geophysical and meteorological measurements. Objects that lie on a sphere experience constraints that are not present in Euclidean (flat) space and that influence both how the particles can be arranged as well as their statistical properties. These constraints, coupled with the curved geometry, require a careful extension of quantities used for the analysis of particle distributions in Euclidean space to distributions confined to the surface of a sphere. Here, we introduce a framework designed to analyze and classify structural order and disorder in particle distributions constrained to the sphere. The classification is based on the concept of hyperuniformity, which was first introduced 15 years ago and since then studied extensively in Euclidean space, yet has only very recently been considered also for spherical surfaces. We employ a generalization of the structure factor on the sphere, related to the power spectrum of the corresponding multipole expansion of particle density distribution. The spherical structure factor is then shown to couple with cap number variance, a measure of density variations at different scales, allowing us to analytically derive different forms of the variance pertaining to different types of distributions. Based on these forms, we construct a classification of hyperuniformity for scale-free particle distributions on the sphere and show how it can be extended to include other distribution types as well. We demonstrate that hyperuniformity on the sphere can be defined either through a vanishing spherical structure factor at low multipole numbers or through a scaling of the cap number variance-in both cases extending the Euclidean definition, while at the same time pointing out crucial differences. Our work thus provides a comprehensive tool for detecting global, long-range order on spheres and for the analysis of spherical computational meshes, biological and synthetic spherical assemblies, and ordering phase transitions in spherically distributed particles.

Moiré and honeycomb lattices through self-assembly of hard-core/soft-shell microgels: experiment and simulation

Moiré and honeycomb lattices result from the sequential double deposition of monolayers of core/shell microgels in dependence of the drying conditions.

Magic number colloidal clusters as minimum free energy structures

Clusters in systems as diverse as metal atoms, virus proteins, noble gases, and nucleons have properties that depend sensitively on the number of constituent particles. Certain numbers are termed ‘magic’ because they grant the system with closed shells and exceptional stability. To this point, magic number clusters have been exclusively found with attractive interactions as present between atoms. Here we show that magic number clusters exist in a confined soft matter system with negligible interactions. Colloidal particles in an emulsion droplet spontaneously organize into a series of clusters with precisely defined shell structures. Crucially, free energy calculations demonstrate that colloidal clusters with magic numbers possess higher thermodynamic stability than those off magic numbers. A complex kinetic pathway is responsible for the efficiency of this system in finding its minimum free energy configuration. Targeting similar magic number states is a strategy towards unique configurations in finite self-organizing systems across the scales.

Magic number cluster with closed shells and increased stability often result from potential energy minimization between attractive atoms or particles. Here, Wang et al. show that such magic number clusters can also result from entropy maximization in colloidal systems with negligible interactions.

Colloidal Metal-Organic Framework Hexapods Prepared from Postsynthesis Etching with Enhanced Catalytic Activity and Rollable Packing

Recent studies on the effect of particle shapes have led to extensive applications of anisotropic colloids as complex materials building blocks. Although much research has been devoted to colloids of convex polyhedral shapes, branched colloids remain largely underexplored because of limited synthesis strategies. Here we achieved the preparation of metal-organic framework (MOF) colloids in a hexapod shape, not directly from growth but from postsynthesis etching of truncated rhombic dodecahedron (TRD) parent particles. To understand the branch development, we used in situ optical microscopy to track the local surface curvature evolution of the colloids as well as facet-dependent etching rate. The hexapods show unique properties, such as improved catalytic activity in a model Knoevenagel reaction likely due to enhanced access to active sites, and the assembly into open structures which can be easily integrated with a self-rolled-up nanomembrane structure. Both the postsynthesis etching and the hexapod colloids demonstrated here show a new route of engineering micrometer-sized building blocks with exotic shapes and intrinsic functionalities originated from the molecular structure of materials.

Conforming nanoparticle sheets to surfaces with Gaussian curvature

Nanoparticle monolayer sheets are ultrathin inorganic-organic hybrid materials that combine highly controllable optical and electrical properties with mechanical flexibility and remarkable strength. Like other thin sheets, their low bending rigidity allows them to easily roll into or conform to cylindrical geometries. Nanoparticle monolayers not only can bend, but also cope with strain through local particle rearrangement and plastic deformation. This means that, unlike thin sheets such as paper or graphene, nanoparticle sheets can much more easily conform to surfaces with complex topography characterized by non-zero Gaussian curvature, like spherical caps or saddles. Here, we investigate the limits of nanoparticle monolayers' ability to conform to substrates with Gaussian curvature by stamping nanoparticle sheets onto lattices of larger polystyrene spheres. Tuning the local Gaussian curvature by increasing the size of the substrate spheres, we find that the stamped sheet morphology evolves through three characteristic stages: from full substrate coverage, where the sheet extends over the interstices in the lattice, to coverage in the form of caps that conform tightly to the top portion of each sphere and fracture at larger polar angles, to caps that exhibit radial folds. Through analysis of the nanoparticle positions, obtained from scanning electron micrographs, we extract the local strain tensor and track the onset of strain-induced dislocations in the particle arrangement. By considering the interplay of energies for elastic and plastic deformations and adhesion, we construct arguments that capture the observed changes in sheet morphology as Gaussian curvature is tuned over two orders of magnitude.

Length segregation in mixtures of spherocylinders induced by imposed topological defects

We explore length segregation in binary mixtures of spherocylinders of lengths L1 and L2 which are tangentially confined on a spherical surface of radius R. The orientation of the spherocylinders is constrained along an externally imposed direction field on the sphere which is either along the longitude or the latitude lines of the sphere. In both situations, integer orientational defects at the poles are imposed. Using computer simulations we show that these topological defects induce a complex segregation picture also depending on the length ratio factor γ = L2/L1 and the total packing fraction η of the spherocylinders. When the binary mixture is aligned along the longitude lines of the sphere, shorter rods tend to accumulate at the topological defects of the polar caps whereas longer rods occupy the central equatorial area of the spherical surface. In the reverse case of latitude ordering, a new state can emerge where longer rods are predominantly both in the cap and in the equatorial areas and shorter rods are localized in between. As a reference situation, we consider a defect-free situation in the flat plane and do not find any length segregation there at similar γ and η; hence, the segregation is purely induced by the imposed topological defects. We also develop an Onsager-like density functional theory which is capable of predicting length segregation in ordered mixtures. At low density, the results of this theory are in good agreement with the simulation data.

Formation of cluster crystals in an ultra-soft potential model on a spherical surface

We investigate the formation of cluster crystals with multiply occupied lattice sites on a spherical surface in systems of ultra-soft particles interacting via repulsive, bounded pair potentials. Not all interactions of this kind lead to clustering: we generalize the criterion devised in C. N. Likos et al., Phys. Rev. E, 2001, 63, 031206 to spherical systems in order to distinguish between cluster-forming systems and fluids which display reentrant melting. We use both DFT and Monte Carlo simulations to characterize the behavior of the system, and obtain semi-quantitative agreement between the two. We find that the number of clusters is determined by the ratio between the size σ of the ultra-soft particles and the radius R of the sphere in such a way that each stable configuration spans a certain interval of σ/R. Furthermore, we study the effect of topological frustration on the system due to the sphere curvature by comparing the properties of disclinations, i.e., clusters with fewer than six neighbors, and non-defective clusters. Disclinations are shown to be less stable, contain fewer particles, and be closer to their neighbors than other lattice points: these properties are explained on the basis of geometric and energetic considerations.