Packing, entropic patchiness, and self-assembly of non-convex colloidal particles: A simulation perspective

Directed assembly of small binary clusters of magnetizable ellipsoids

We report the effect of shape anisotropy and material properties on the directed assembly of binary suspensions composed of magnetizable ellipsoids. In a Monte Carlo simulation, we implement the ellipsoid-dipole model to calculate the pairwise dipolar interaction energy as a function of position and orientation. The analysis explores dilute suspensions of paramagnetic and diamagnetic ellipsoids with different aspect ratios in a superparamagnetic medium. We analyze the local order of binary structures as a function of particle aspect ratio, medium permeability, and dipolar interaction strength. Our results show that local order and symmetry are tunable under the influence of a uniform magnetic field when one component of the structure is dilute with respect to the other. The simulation results match previously reported experiments on the directed assembly of binary suspension of spheres. Additionally, we report the conditions on particle aspect ratios and medium properties for various structures with rotational symmetries, as well as open and enclosed structures under the influence of a uniform magnetic field.

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.

Mesoscale simulations of diffusion and sedimentation in shape-anisotropic nanoparticle suspensions

We determine the long-time self-diffusion coefficient and sedimentation coefficient for suspensions of nanoparticles with anisotropic shapes (octahedra, cubes, tetrahedra, and spherocylinders) as a function of nanoparticle concentration using mesoscale simulations. We use a discrete particle model for the nanoparticles, and we account for solvent-mediated hydrodynamic interactions between nanoparticles using the multiparticle collision dynamics method. Our simulations are compared to theoretical predictions and experimental data from existing literature, demonstrating good agreement in the majority of cases. Further, we find that the self-diffusion coefficient of the regular polyhedral shapes can be estimated from that of a sphere whose diameter is the average of their inscribed and circumscribed sphere diameters.

Porous Self-Assemblies Mediated by Dumbbell Particles as Cross-Linking Agents

Self-assembly of colloidal particles is emerging as a promising approach for producing novel materials. These colloidal particles can be synthesized with protrusions (lobes) on their surfaces that allow the formation of porous structures with a wide range of applications. Using Langevin dynamics simulations, we studied self-assembly in the binary mixtures of lobed colloidal particles with variations in their lobe sizes to investigate the feasibility of using dumbbell particles (with two lobes) as cross-linkers to increase the porosity in self-assembled morphologies. Each binary system was formed by mixing the dumbbell particles with one of the following types of particles: trigonal planar (three lobes), tetrahedral (four lobes), trigonal bipyramidal (five lobes), and octahedral (six lobes). We observed that the lobe size on each particle can be tuned to favor the formation of random aggregates and spherical aggregates when the lobes are larger and well-ordered crystalline structures when the lobes are smaller. We also observed that these polydisperse systems form self-assembled structures characterized by porosities higher than those of the structures formed by the monodisperse systems. These results indicate that the lobe size is an important design feature that can be optimized to achieve desired structures with distinct morphologies and porosities, and the dumbbell particles are effective cross-linking agents to enhance the porosity in self-assembled structures.

Predicting Outcomes of Nanoparticle Attachment by Connecting Atomistic, Interfacial, Particle, and Aggregate Scales

Predicting nanoparticle aggregation and attachment phenomena requires a rigorous understanding of the interplay among crystal structure, particle morphology, surface chemistry, solution conditions, and interparticle forces, yet no comprehensive picture exists. We used an integrated suite of experimental, theoretical, and simulation methods to resolve the effect of solution pH on the aggregation of boehmite nanoplatelets, a case study with important implications for the environmental management of legacy nuclear waste. Real-time observations showed that the particles attach preferentially along the (010) planes at pH 8.5 and the (101) planes at pH 11. To rationalize these results, we established the connection between key physicochemical phenomena across the relevant length scales. Starting from molecular-scale simulations of surface hydroxyl reactivity, we developed an interfacial-scale model of the corresponding electrostatic potentials, with subsequent particle-scale calculations of the resulting driving forces allowing successful prediction of the attachment modes. Finally, we scaled these phenomena to understand the collective structure at the aggregate-scale. Our results indicate that facet-specific differences in surface chemistry produce heterogeneous surface charge distributions that are coupled to particle anisotropy and shape-dependent hydrodynamic forces, to play a key role in controlling aggregation behavior.

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.

Phase behavior of hard circular arcs

By using Monte Carlo numerical simulation, this work investigates the phase behavior of systems of hard infinitesimally thin circular arcs, from an aperture angle θ→0 to an aperture angle θ→2π, in the two-dimensional Euclidean space. Except in the isotropic phase at lower density and in the (quasi)nematic phase, in the other phases that form, including the isotropic phase at higher density, hard infinitesimally thin circular arcs autoassemble to form clusters. These clusters are either filamentous, for smaller values of θ, or roundish, for larger values of θ. Provided the density is sufficiently high, the filaments lengthen, merge, and straighten to finally produce a filamentary phase while the roundels compact and dispose themselves with their centers of mass at the sites of a triangular lattice to finally produce a cluster hexagonal phase.

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.

Chiral self-sorting of active semiflexible filaments with intrinsic curvature

Many-body interactions in systems of active matter can cause particles to move collectively and self-organize into dynamic structures with long-range order. In cells, the self-assembly of cytoskeletal filaments is critical for cellular motility, structure, intracellular transport, and division. Semiflexible cytoskeletal filaments driven by polymerization or motor-protein interactions on a two-dimensional substrate, such as the cell cortex, can induce filament bending and curvature leading to interesting collective behavior. For example, the bacterial cell-division filament FtsZ is known to have intrinsic curvature that causes it to self-organize into rings and vortices, and recent experiments reconstituting the collective motion of microtubules driven by motor proteins on a surface have observed chiral symmetry breaking of the collective behavior due to motor-induced curvature of the filaments. Previous work on the self-organization of driven filament systems have not studied the effects of curvature and filament structure on collective behavior. In this work, we present Brownian dynamics simulation results of driven semiflexible filaments with intrinsic curvature and investigate how the interplay between filament rigidity and radius of curvature can tune the self-organization behavior in homochiral systems and heterochiral mixtures. We find a curvature-induced reorganization from polar flocks to self-sorted chiral clusters, which is modified by filament flexibility. This transition changes filament transport from ballistic to diffusive at long timescales.

Dense packings of hard circular arcs

This work investigates dense packings of congruent hard infinitesimally thin circular arcs in the two-dimensional Euclidean space. It focuses on those denotable as major whose subtended angle θ∈(π,2π]. Differently than those denotable as minor whose subtended angle θ∈[0,π], it is impossible for two hard infinitesimally thin circular arcs with θ∈(π,2π] to arbitrarily closely approach once they are arranged in a configuration, e.g., on top of one another, replicable ad infinitum without introducing any overlap. This makes these hard concave particles, in spite of being infinitesimally thin, most densely pack with a finite number density. This raises the question as to what are these densest packings and what is the number density that they achieve. Supported by Monte Carlo numerical simulations, this work shows that one can analytically construct compact closed circular groups of hard major circular arcs in which a specific, θ-dependent, number of them (counter) clockwise intertwine. These compact closed circular groups then arrange on a triangular lattice. These analytically constructed densest-known packings are compared to corresponding results of Monte Carlo numerical simulations to assess whether they can spontaneously turn up.

Role of Entropy in Colloidal Self-Assembly

Entropy plays a key role in the self-assembly of colloidal particles. Specifically, in the case of hard particles, which do not interact or overlap with each other during the process of self-assembly, the free energy is minimized due to an increase in the entropy of the system. Understanding the contribution of entropy and engineering it is increasingly becoming central to modern colloidal self-assembly research, because the entropy serves as a guide to design a wide variety of self-assembled structures for many technological and biomedical applications. In this work, we highlight the importance of entropy in different theoretical and experimental self-assembly studies. We discuss the role of shape entropy and depletion interactions in colloidal self-assembly. We also highlight the effect of entropy in the formation of open and closed crystalline structures, as well as describe recent advances in engineering entropy to achieve targeted self-assembled structures.

When does entropy promote local organization?

Crowded soft-matter and biological systems organize locally into preferred motifs. Locally-organized motifs in soft systems can, paradoxically, arise from a drive to maximize overall system entropy. Entropy-driven local order has been directly confirmed in model, synthetic colloidal systems, however similar patterns of organization occur in crowded biological systems ranging from the contents of a cell to collections of cells. In biological settings, and in soft matter more broadly, it is unclear whether entropy generically promotes or inhibits local organization. Resolving this is difficult because entropic effects are intrinsically collective, complicating efforts to isolate them. Here, we employ minimal models that artificially restrict system entropy to show that entropy drives systems toward local organization, even when the model system entropy is below reasonable physical bounds. By establishing this bound, our results suggest that entropy generically promotes local organization in crowded soft and biological systems of rigid objects.

Multivalent Patchy Colloids for Quantitative 3D Self-Assembly Studies

We report methods to synthesize sub-micron- and micron-sized patchy silica particles with fluorescently labeled hemispherical titania protrusions, as well as routes to efficiently characterize these particles and self-assemble these particles into non-close-packed structures. The synthesis methods expand upon earlier work in the literature, in which silica particles packed in a colloidal crystal were surface-patterned with a silane coupling agent. Here, hemispherical amorphous titania protrusions were successfully labeled with fluorescent dyes, allowing for imaging by confocal microscopy and super-resolution techniques. Confocal microscopy was exploited to experimentally determine the numbers of protrusions per particle over large numbers of particles for good statistical significance, and these distributions were compared to simulations predicting the number of patches as a function of core particle polydispersity and maximum separation between the particle surfaces. We self-assembled these patchy particles into open percolating gel networks by exploiting solvophobic attractions between the protrusions.

Colloidal molecules and patchy particles: complementary concepts, synthesis and self-assembly

About the latest developments regarding self-assembly of textured colloids and its prospects.

The entropic bond in colloidal crystals

A vast array of natural phenomena can be understood through the long-established schema of chemical bonding. Conventional chemical bonds arise through local gradients resulting from the rearrangement of electrons; however, it is possible that the hallmark features of chemical bonding could arise through local gradients resulting from nonelectronic forms of mediation. If other forms of mediation give rise to "bonds" that act like conventional ones, recognizing them as bonds could open new forms of supramolecular descriptions of phenomena at the nano- and microscales. Here, we show via a minimal model that crowded hard-particle systems governed solely by entropy exhibit the hallmark features of bonding despite the absence of chemical interactions. We quantitatively characterize these features and compare them to those exhibited by chemical bonds to argue for the existence of entropic bonds. As an example of the utility of the entropic bond classification, we demonstrate the nearly equivalent tradeoff between chemical bonds and entropic bonds in the colloidal crystallization of hard hexagonal nanoplates.

Preparation of Dimpled Polystyrene-Silica Colloidal Nanocomposite Particles

Preparation of polymer-silica colloidal nanocomposite particles with concave shape is challenging and seldom reported. This paper presents a novel and facile method to prepare dimpled polymer-silica nanocomposite particles with a thin silica shell through the judicious combination of alcoholic dispersion polymerization and the decane evaporation method. Submicrometer-sized polystyrene-silica (PS-SiO₂) nanocomposite particles were first prepared by dispersion polymerization of styrene in methanol in the presence of a methanolic silica sol, and then dimpled PS-SiO₂ particles were prepared by heating the near-spherical PS-SiO₂ particles dispersed in methanol/water media in the presence of decane and subsequent cooling. The effects of different heat treatment parameters, such as methanol/water ratio, stirring temperature, and stirring rate on the formation of the nanocomposite particles were investigated. Optimization of the heating conditions allowed ∼100% of dimpled particles to be achieved with one dimple on each particle. Moreover, calcination of the dimpled PS-SiO₂ nanocomposite particles led to the formation of hollow dimpled particles with a thin silica shell. This method is expected to enrich the shapes of polymer-silica nanocomposite particles.

Phase Behavior of Bowl-Shaped Colloids: Order and Dynamics in Plastic Crystals and Glasses

Charged fluorescent bowl-shaped colloids consisting of a polystyrene core surrounded by a poly(N-isopropylmethacrylamide) shell are obtained by nanoengineering spherical composite microgels. The phase diagram of these soft bowl-shaped colloids interacting through long-range Yukawa-type interactions is investigated using confocal laser scanning microscopy. The bowl-shaped structure leads to marked differences in phase-behavior compared to their spherical counterpart. With increasing number density, a transition from a fluid to a plastic crystal phase, with freely rotating particles, followed by a glass-like state is observed. It is found that the anisotropic bowl shape frustrates crystallization and slows down crystallization kinetics and causes the glass-like transition to shift to a significantly lower volume fraction than for the spheres. Quantitative analysis of the positional and orientational order demonstrates that the plastic crystal phase exhibits quasi-long range translational order and orientational disorder, while in the disordered glass-like phase the long-range translational order vanishes and short-range rotational order appears, dictated by the specific bowl shape. It is further shown that the different structural transitions are characterized by decoupling of the translational and orientational dynamics.

Hard convex lens-shaped particles: metastable, glassy and jammed states

We generate and study dense positionally and/or orientationally disordered, including jammed, monodisperse packings of hard convex lens-shaped particles (lenses). Relatively dense isotropic fluid configurations of lenses of various aspect ratios are slowly compressed via a Monte Carlo method based procedure. Under this compression protocol, while 'flat' lenses form a nematic fluid phase (where particles are positionally disordered but orientationally ordered) and 'globular' lenses form a plastic solid phase (where particles are positionally ordered but orientationally disordered), 'intermediate', neither 'flat' nor 'globular', lenses do not form either mesophase. In general, a crystal solid phase (where particles are both positionally and orientationally ordered) does not spontaneously form during lengthy numerical simulation runs. In correspondence to those volume fractions at which a transition to the crystal solid phase would occur in equilibrium, a 'downturn' is observed in the inverse compressibility factor versus volume fraction curve beyond which this curve behaves essentially linearly. This allows us to estimate the volume fraction at jamming of the dense non-crystalline packings so generated. These packings are nematic for 'flat' lenses and plastic for 'globular' lenses, while they are robustly isotropic for 'intermediate' lenses, as confirmed by the calculation of the τ order metric, among other quantities. The structure factors S(k) of the corresponding jammed states tend to zero as the wavenumber k goes to zero, indicating they are effectively hyperuniform (i.e., their infinite-wavelength density fluctuations are anomalously suppressed). Among all possible lens shapes, 'intermediate' lenses with aspect ratio around 2/3 are special because they are those that reach the highest volume fractions at jamming while being positionally and orientationally disordered and these volume fractions are as high as those reached by nematic jammed states of 'flat' lenses and plastic jammed states of 'globular' lenses. All of their attributes, taken together, make such 'intermediate' lens packings particularly good glass-forming materials.

A Shape-Induced Orientation Phase within 3D Nanocrystal Solids

When nanocrystals self assemble into ordered superstructures they form functional solids that may inherit the electronical properties of the single nanocrystals. To what extent these properties are enhanced depends on the positional and orientational order of the nanocrystals within the superstructure. Here, the formation of micrometer-sized free-standing supercrystals of faceted 20 nm Bi nanocrystals is investigated. The self-assembly process, induced by nonsolvent into solvent diffusion, is probed in situ by synchrotron X-ray scattering. The diffusion-gradient is identified as the critical parameter for controlling the supercrystal-structure as well as the alignment of the supercrystals with respect to the substrate. Monte Carlo simulations confirm the positional order of the nanocrystals within these superstructures and reveal a unique orientation phase: the nanocrystal shape, determined by the atomic Bi crystal structure, induces a total of 6 global orientations based on facet-to-facet alignment. This parallel alignment of facets is a prerequisite for optimized electronic and optical properties within designed nanocrystal solids.

Many Body Effects and Icosahedral Order in Superlattice Self-Assembly

We elucidate how nanocrystals "bond" to form ordered structures. For that purpose we consider nanocrystal configurations consisting of regular polygons and polyhedra, which are the motifs that constitute single component and binary nanocrystal superlattices, and simulate them using united atom models. We compute the free energy and quantify many body effects, i.e., those that cannot be accounted for by pair potential (two-body) interactions, further showing that they arise from coalescing vortices of capping ligands. We find that such vortex textures exist for configurations with local coordination number ≤6. For higher coordination numbers, vortices are expelled and nanocrystals arrange in configurations with tetrahedral or icosahedral order. We provide explicit formulas for the optimal separations between nanocrystals, which correspond to the minima of the free energies. Our results quantitatively explain the structure of superlattice nanocrystals as reported in experiments and reveal how packing arguments, extended to include soft components, predict ordered nanocrystal aggregation.

Relevance of packing to colloidal self-assembly

Since the 1920s, packing arguments have been used to rationalize crystal structures in systems ranging from atomic mixtures to colloidal crystals. Packing arguments have recently been applied to complex nanoparticle structures, where they often, but not always, work. We examine when, if ever, packing is a causal mechanism in hard particle approximations of colloidal crystals. We investigate three crystal structures composed of their ideal packing shapes. We show that, contrary to expectations, the ordering mechanism cannot be packing, even when the thermodynamically self-assembled structure is the same as that of the densest packing. We also show that the best particle shapes for hard particle colloidal crystals at any finite pressure are imperfect versions of the ideal packing shape.

Packing of sedimenting equiaxed dendrites

The packing of free-floating crystal grains during solidification has a strong impact on the phase-change process as well as on the structure and the defects in the solidified material. The packing fraction is affected by the particular dendritic morphology of the grains and by their low inertia resulting from the small density difference between solid and liquid. Understanding the grain packing phenomenon during metal alloy solidification is not experimentally possible since packing is coupled to many other phenomena. We therefore investigate the packing of equiaxed dendrites on a model system, consisting of fixed-shape nonconvex model particles sedimenting in conditions hydrodynamically similar to those encountered in solidifying metals. We perform numerical simulations by a discrete-element model and experiments with transparent liquids in a sedimentation column. The combination of experiments and simulations enables us to determine the packing fraction as a function of (i) the grain morphology, expressed by a shape parameter, and (ii) the hydrodynamic conditions, expressed by the particle Stokes number.

Macroscopic chiral symmetry breaking in monolayers of achiral nonconvex platelets

The fabrication of chiral structures using achiral building blocks is a fundamental problem that remains a challenge in materials science. In this work we present a molecular dynamics simulation study of nonconvex polygonal platelets, interacting via soft-repulsive interactions, that are confined in two-dimensional space. These particle models are designed to promote, even at moderate densities, a natural offset displacement between the edges of neighbouring particles. In particular we demonstrate that nonconvex platelets exhibit macroscopic chiral symmetry breaking when the symmetry of the particles equals (or is multiple of) the number of nearest neighbours in the condensed crystalline phase, corresponding to the situation of platelets with 4-, 6-, and 12-fold symmetries.