Glassy dynamics, spinodal fluctuations, and the kinetic limit of nucleation in suspensions of colloidal hard rods

Diffusionless rotator-crystal transitions in colloidal truncated cubes

Upon osmotic compression, rotationally symmetric faceted colloidal particles can form translationally ordered, orientationally disordered rotator mesophases. This study explores the mechanism of rotator-to-crystal phase transitions where orientational order is gained in a translationally ordered phase, using rotator-phase forming truncated cubes as a testbed. Monte Carlo simulations were conducted for two selected truncations (s), one for s = 0.527 where the rotator and crystal lattices are dissimilar and one for s = 0.572 where the two phases have identical lattices. These differences set the stage for a qualitative difference in their rotator-crystal transitions, highlighting the effect of lattice distortion on phase transition kinetics. Our simulations reveal that significant lattice deviatoric effects could hinder the rotator-to-crystal transition and favor arrangements of lower packing fraction instead. Indeed, upon compression, it is found that for s = 0.527, the rotator phase does not spontaneously transition into the stable, densely packed crystal due to the high lattice strains involved but instead transitions into a metastable solid phase to be colloquially referred to as "orientational salt" for short, which has a similar lattice as the rotator phase and exhibits two distinct particle orientations having substitutional order, alternating regularly throughout the system. This study paves the way for further analysis of diffusionless transformations in nanoparticle systems and how lattice-distortion could influence crystallization kinetics.

Comparing brute force to transition path sampling for gas hydrate nucleation with a flat interface: comments on time reversal symmetry

Fluid to solid nucleation is often investigated with the rare event method transition path sampling (TPS). I claim that the inherent irreversibility of solid nucleation, even at stationary conditions, calls into question TPS's applicability for determining solid nucleation mechanisms, especially for pre-critical behavior. Even when applied to a phenomenon which displays time reversal asymmetry like solid nucleation, TPS is a good means of exploring phase space and giving trends in post-critical structure, and its ability to facilitate nucleation rate and free energy calculations remains outstanding. Forward-only splitting and ratcheting methods such as forward flux sampling are more attractive for understanding nucleation mechanisms as they do not require time reversal symmetry, but at low driving forces may suffer from the same limitations as brute force: they may never make it to the first ratchet. Here I briefly summarize the TPS method and gas hydrate nucleation simulation literature, focusing on topics within both to facilitate a comparison of brute force hydrate nucleation to transition path sampling of hydrate nucleation. Perhaps anecdotally, the brute force technique results in more crystalline trajectories despite having higher driving forces than TPS. I maintain this difference is because of the inherent irreversibility of hydrate nucleation, meaning its pre-critical behavior cannot accurately be determined by the melting trajectories that comprise approximately half of the configurations in TPS's path ensemble.

How does a hyperuniform fluid freeze?

All phase transitions can be categorized into two different types: continuous and discontinuous phase transitions. Discontinuous phase transitions are normally accompanied with significant structural changes, and nearly all of them have the kinetic pathway of nucleation and growth, if the system does not suffer from glassy dynamics. Here, in a system of barrier-controlled reactive particles, we find that the discontinuous freezing transition of a nonequilibrium hyperuniform fluid into an absorbing state does not have the kinetic pathway of nucleation and growth, and the transition is triggered by long-wavelength fluctuations. The transition rate decreases with increasing the system size, which suggests that the metastable hyperuniform fluid could be kinetically stable in an infinitely large system. This challenges the common understanding of metastability in discontinuous phase transitions. Moreover, we find that the "metastable yet kinetically stable" hyperuniform fluid features a scaling in the structure factor [Formula: see text] in 2D, which is the third dynamic hyperuniform state in addition to the critical hyperuniform state with [Formula: see text] and the nonequilibrium hyperuniform fluid with [Formula: see text].

Effect of particle anisotropy on the thermodynamics and kinetics of ordering transitions in hard faceted particles

Monte Carlo simulations were used to study the influence of particle aspect ratio on the kinetics and phase behavior of hard gyrobifastigia (GBF). First, the formation of a highly anisotropic nucleus shape in the isotropic-to-crystal transition in regular GBF is explained by the differences in interfacial free energies of various crystal planes and the nucleus geometry predicted by the Wulff construction. GBF-related shapes with various aspect ratios were then studied, mapping their equations of state, determining phase coexistence conditions via interfacial pinning, and computing nucleation free-energy barriers via umbrella sampling using suitable order parameters. Our simulations reveal a reduction of the kinetic barrier for isotropic-crystal transition upon an increase in aspect ratio, and that for highly oblate and prolate aspect ratios, an intermediate nematic phase is stabilized. Our results and observations also support two conjectures for the formation of the crystalline state from the isotropic phase: that low phase free energies at the ordering phase transition correlate with low transition barriers and that the emergence of a mesophase provides a steppingstone that expedites crystallization.

Topological fine structure of smectic grain boundaries and tetratic disclination lines within three-dimensional smectic liquid crystals

Observing and characterizing the complex ordering phenomena of liquid crystals subjected to external constraints constitutes an ongoing challenge for chemists and physicists alike. To elucidate the delicate balance appearing when the intrinsic positional order of smectic liquid crystals comes into play, we perform Monte-Carlo simulations of rod-like particles in a range of cavities with a cylindrical symmetry. Based on recent insights into the topology of smectic orientational grain boundaries in two dimensions, we analyze the emerging three-dimensional defect structures from the perspective of tetratic symmetry. Using an appropriate three-dimensional tetratic order parameter constructed from the Steinhardt order parameters, we show that those grain boundaries can be interpreted as a pair of tetratic disclination lines that are located on the edges of the nematic domain boundary. Thereby, we shed light on the fine structure of grain boundaries in three-dimensional confined smectics.

Topology of Orientational Defects in Confined Smectic Liquid Crystals

We propose a general formalism to characterize orientational frustration of smectic liquid crystals in confinement by interpreting the emerging networks of grain boundaries as objects with a topological charge. In a formal idealization, this charge is distributed in pointlike units of quarter-integer magnitude, which we identify with tetratic disclinations located at the end points and nodes. This coexisting nematic and tetratic order is analyzed with the help of extensive Monte Carlo simulations for a broad range of two-dimensional confining geometries as well as colloidal experiments, showing how the observed defect networks can be universally reconstructed from simple building blocks. We further find that the curvature of the confining wall determines the anchoring behavior of grain boundaries, such that the number of nodes in the emerging networks and the location of their end points can be tuned by changing the number and smoothness of corners, respectively.

Multistep nucleation of anisotropic molecules

Phase transition of anisotropic materials is ubiquitously observed in physics, biology, materials science, and engineering. Nevertheless, how anisotropy of constituent molecules affects the phase transition dynamics is still poorly understood. Here we investigate numerically the phase transition of a simple model system composed of anisotropic molecules, and report on our discovery of multistep nucleation of nuclei with layered positional ordering (smectic ordering), from a fluid-like nematic phase with orientational order only (no positional order). A trinity of molecular dynamics simulation, machine learning, and molecular cluster analysis yielding free energy landscapes unambiguously demonstrates the dynamics of multistep nucleation process involving characteristic metastable clusters that precede supercritical smectic nuclei and cannot be accounted for by the classical nucleation theory. Our work suggests that molecules of simple shape can exhibit rich and complex nucleation processes, and our numerical approach will provide deeper understanding of phase transitions and resulting structures in anisotropic materials such as biological systems and functional materials.

Crystal nucleation in colloidal rod suspensions: The effect of depletant size

In order to better control the assembly of nanorods, knowledge of the pathways by which they form ordered structures is desirable. In this paper, we characterize crystal nucleation in suspensions of spherocylindrical rods with aspect ratio L/D = 2.3 in the presence of both small and large polymer depletants. Using a combination of Langevin dynamics and Monte Carlo simulations, together with biased sampling techniques, we show that the preferred pathway always involves the formation of monolayer assemblies irrespective of the volume fraction of the initial isotropic phase and the diameter of the depletants. This includes the previously neglected case of nucleation from the colloidal liquid phase and shows that the presence of depletion attraction can alter nucleation pathways even when the initial phase is dense.

Low Interfacial Free Energy Describes the Bulk Ordering Transition in Colloidal Cubes

Many hard faceted nanoparticles are known to undergo disorder-to-order phase transitions following a classical nucleation and growth mechanism. In a previous study [J. Phys. Chem. B 2018, 122, 9264-9273], it was shown that hard cubes undergo a nonclassical phase transition with a bulk character instead of originating from consolidated nuclei. Significantly, an unusually high fraction of ordered particles was observed in the metastable basin of the disordered phase, even for very low degrees of supersaturation. This work aims to substantiate the conjecture that these unique properties originate from a comparatively low interfacial free energy between the disordered and ordered phases for hard cubes relative to other hard particle systems. Using the cleaving wall method to directly measure the interfacial free energy for cubes, it is found that its values are indeed small; e.g., at phase coexistence conditions, it is only one-fifth that for hard spheres. A theoretical nucleation model is used to explore the broader implications of low interfacial tension values and how this could result in a bulk ordering mechanism.

Translational and rotational critical-like behaviors in the glass transition of colloidal ellipsoid monolayers

Critical-like behaviors have been found in translational degrees of freedom near the glass transition of spherical particle systems mainly with local polycrystalline structures, but it is not clear if criticality exists in more general glassy systems composed of nonspherical particles without crystalline structures. Here, through experiments and simulations, we show critical-like behaviors in both translational and rotational degrees of freedom in monolayers of monodisperse colloidal ellipsoids in the absence of crystalline orders. We find rich features of the Ising-like criticality in structure and slow dynamics at the ideal glass transition point ϕ₀, showing the thermodynamic nature of glass transition at ϕ₀ A dynamic criticality is found at the mode-coupling critical point ϕ_c for the fast-moving clusters whose critical exponents increase linearly with fragility, reflecting a dynamic glass transition. These results cast light on the glass transition and explain the mystery that the dynamic correlation lengths diverge at two different temperatures.

Hyperspectral Imaging of Photonic Cellulose Nanocrystal Films: Structure of Local Defects and Implications for Self-Assembly Pathways

Cellulose nanocrystals (CNCs) can spontaneously assemble into chiral nematic films capable of reflecting circularly polarized light in the visible range. As many other photonic materials obtained by bottom-up approaches, CNC films often display defects that greatly impact their visual appearance. Here, we study the optical response of defects in photonic CNC films, coupling optical microscopy with hyperspectral imaging, and we compare it to optical simulations of discontinuous cholesteric structures of increasing complexity. Cross-sectional SEM observations of the film structure guided the choice of simulation parameters and showed excellent agreement with experimental optical patterns. More importantly, it strongly suggests that the last fraction of CNCs to self-assemble, upon solvent evaporation, does not undergo the typical nucleation and growth pathway, but a spinodal decomposition, an alternative self-assembly pathway so far overlooked in cast films and that can have far-reaching consequences on choices of CNC sources and assembly conditions.

Supercrystallographic Reconstruction of 3D Nanorod Assembly with Collectively Anisotropic Upconversion Fluorescence

Constructing three-dimensional (3D) metamaterials from functional nanoparticles endows them with emerging collective properties tailored by the packing geometries. Herein, we report 3D supercrystals self-assembled from upconversion nanorods (NaYF₄:Yb,Er NRs), which exhibit both translational ordering of NRs and orientational ordering between constituent NRs in the superlattice (SL). The construction of 3D reciprocal space mappings (RSMs) based on synchrotron-based X-ray scattering measurements was developed to uncover the complex structure of such an assembly. That is, the two main orthogonal sets of hexagonal close-packing (hcp)-like SLs share the [110]_SL axis, and NRs within the SL possess orientational relationships of [120]_NR//[100]_SL, [210]_NR//[010]_SL, and [001]_NR//[001]_SL. Notably, these supercrystals containing well-aligned NRs exhibit collectively anisotropic upconversion fluorescence in two perpendicular directions. This study not only demonstrates novel crystalline superstructures and functionality of NR-based 3D assemblies but also offers a unique tool for deciphering a wide range of complex nanoparticle supercrystals.

Length-scales of dynamic heterogeneity in a driven binary colloid

Here we study the characteristic length scales in an aqueous suspension of a symmetric oppositely charged colloid subjected to a uniform electric field by Brownian dynamics simulations. We consider the in-plane structure in the presence of a sufficiently strong electric field where the like charges in the system form macroscopic lanes. We construct spatial correlation functions characterizing the structural order and that of particles of different mobilities in the plane transverse to the electric field at a given time. We call these functions equal time density correlation functions (ETDCFs). The ETDCFs between particles of different charges, irrespective of mobilities, are the structural ETDCFs, while those between particles of different mobilities are the dynamic ETDCFs. We extract the characteristic length of correlation by fitting the envelopes of the ETDCFs to exponential dependences. We find that the correlation length scales of the structural ETDCFs and the dynamic ETDCFs of the slow particles increase with time in a concurrent manner. This suggests that the clustering of particles tends to build up dynamically correlated slow particles in the plane transverse to the lanes. The ETDCFs can be measured for colloidal systems by directly following the particle motion by video-microscopy and may be useful to understand the patterns out of equilibrium.

Self-assembly of anisotropic nanoparticles into functional superstructures

Self-assembly of colloidal nanoparticles (NPs) into superstructures offers a flexible and promising pathway to manipulate the nanometer-sized particles and thus make full use of their unique properties. This bottom-up strategy builds a bridge between the NP regime and a new class of transformative materials across multiple length scales for technological applications. In this field, anisotropic NPs with size- and shape-dependent physical properties as self-assembly building blocks have long fascinated scientists. Self-assembly of anisotropic NPs not only opens up exciting opportunities to engineer a variety of intriguing and complex superlattice architectures, but also provides access to discover emergent collective properties that stem from their ordered arrangement. Thus, this has stimulated enormous research interests in both fundamental science and technological applications. This present review comprehensively summarizes the latest advances in this area, and highlights their rich packing behaviors from the viewpoint of NP shape. We provide the basics of the experimental techniques to produce NP superstructures and structural characterization tools, and detail the delicate assembled structures. Then the current understanding of the assembly dynamics is discussed with the assistance of in situ studies, followed by emergent collective properties from these NP assemblies. Finally, we end this article with the remaining challenges and outlook, hoping to encourage further research in this field.

Anisotropic viscoelastic phase separation in polydisperse hard rods leads to nonsticky gelation

Spinodal demixing into two phases having very different viscosities leads to viscoelastic networks-i.e., gels-usually as a result of attractive particle interactions. Here, however, we demonstrate demixing in a colloidal system of polydisperse, rod-like clay particles that is driven by particle repulsions instead. One of the phases is a nematic liquid crystal with a highly anisotropic viscosity, allowing flow along the director, but suppressing it in other directions. This phase coexists with a dilute isotropic phase. Real-space analysis and molecular-dynamics simulations both reveal a long-lived network structure that is locally anisotropic, yet macroscopically isotropic. We show that our system exhibits the characteristics of colloidal gelation, leading to nonsticky gels.

Dynamic arrest of adhesive hard rod dispersions

The phenomenon of dynamic arrest, more commonly referred to as gel and glass formation, originates as particle motion slows significantly. Current understanding of gels and glasses stems primarily from dispersions of spherical particles, but much less is known about how particle shape affects dynamic arrest transitions. To better understand the effects of particle shape anisotropy on gel and glass formation, we systematically measure the rheology, particle dynamics, and static microstructure of thermoreversible colloidal dispersions of adhesive hard rods (AHR). First, the dynamic arrest transitions are mapped as a function of temperature T, aspect ratio L/D≈ 3 to 7, and volume fraction φ≈ 0.1 to 0.5. The critical gel temperature T_gel and glass volume fraction φ_g are determined from the particle dynamics and rheology. Second, an effective orientation-averaged, short-range attraction between rods is quantified from small-angle scattering measurements and characterized by a reduced temperature τ. Similar τ is found at low rod concentrations, indicating that rod gelation occurs at similar effective attraction strength independent of L/D. Monte Carlo simulations reveal a similar convergence in τ when rods cluster and percolate with an average bond coordination number 〈n_c〉≈ 2.4, supporting the link between physical gelation and rigidity percolation. Lastly, AHR results are mapped onto a dimensionless state diagram to compare with previous predictions of attraction-driven gels, repulsion-driven glasses, and liquid crystal phases.

Short and Soft: Multidomain Organization, Tunable Dynamics, and Jamming in Suspensions of Grafted Colloidal Cylinders with a Small Aspect Ratio

The yet virtually unexplored class of soft colloidal rods with a small aspect ratio is investigated and shown to exhibit a very rich phase and dynamic behavior, spanning from liquid to nearly melt state. Instead of the nematic order, these short and soft nanocylinders alter their organization with increasing concentration from isotropic liquid with random orientation to small domains with preferred local orientation and eventually a multidomain arrangement with a local orientational order. The latter gives rise to a kinetically suppressed state akin to structural glass with detectable terminal relaxation, which, on further increasing concentration, reveals features of hexagonally packed order as in ordered block copolymers. The respective dynamic response comprises four regimes, all above the overlapping concentration of 0.02 g/mL:(I) from 0.03 to 0.1 g/mol, the system undergoes a liquid-to-solidlike transition with a structural relaxation time that grows by 4 orders of magnitude. (II) From 0.1 to 0.2 g/mL, a dramatic slowing-down is observed and is accompanied by an evolution from isotropic to a multidomain structure. (III) Between 0.2 and 0.6 g/mol, the suspensions exhibit signatures of shell interpenetration and jamming, with the colloidal plateau modulus depending linearly on concentration. (IV) At 0.74 g/mL, in the densely jammed state, the viscoelastic signature of hexagonally packed cylinders from microphase-separated block copolymers is detected. These properties set short and soft nanocylinders apart from long colloidal rods (with a large aspect ratio) and provide insights for fundamentally understanding the physics in this intermediate soft colloidal regime and for tailoring the flow properties of nonspherical soft colloids.

A versatile simulation method for studying phase behavior and dynamics in colloidal rod and rod-polymer suspensions

Here, we present an implicit-solvent model for dynamic simulations of hard-rod and rod-polymer suspensions. Individual rods are represented by a rigid linear chain consisting of overlapping spheres which interact through a pseudohard-core potential based on the cut-and-shifted Mie (generalized Lennard-Jones) potential with exponents (50, 49). In the rod-polymer suspensions, the polymers are modeled as freely interpenetrable spheres with respect to each other, while there is the pseudohard-core repulsion between the polymer and rod spheres. Dynamic simulations with this model are carried out with a dissipative particle dynamics (DPD) thermostat-each sphere is put in a larger DPD sphere and thus interacts with others via additional pairwise frictional and random forces-which captures the effects of Brownian forces due to the solvent while conserving local momentum. The phase behavior of these models, obtained from continuous compression and expansion simulations, reproduces previous predictions based on theoretical calculations and Monte Carlo simulations. Our method is suited to study dynamic processes in these suspensions, including nucleation and self-assembly, and can be readily extended to colloidal particles of different shapes and chemistry.

Nucleation of pseudo hard-spheres and dumbbells at moderate metastability: appearance of A15 Frank–Kasper phase at intermediate elongations

Crystal nucleation of repulsive hard-dumbbells from the sphere to the two tangent spheres limit is investigated at moderately high metastability by brute-force molecular dynamics simulations.

Micro-structure and motion of two-dimensional dense short spherocylinder liquids

We numerically investigate the micro-structure and motion of 2D liquids composed of dense short spherocylinders, by reducing the shape aspect ratio from 3. It is found that reducing shape aspect ratio from 3 causes a smooth transition from heterogeneous structures composed of crystalline ordered domains with good tetratic alignment order to those with good hexagonal bond-orientational order at an aspect ratio equaling 1.35. In the intermediate regime, both structural orders are strongly deteriorated, and the translational hopping rate reaches a maximum due to the poor particle interlocking of the disordered structure. Shortening rod length allows easier rotation, induces monotonic increase of rotational hopping rates, and resumes the separation of rotational and translational hopping time scales at the small aspect ratio end, after the crossover of their rates in the intermediate regime. At the large shape aspect ratio end, the poor local tetratic order has the same positive effects on facilitating local rotational and translational hopping. In contrast, at the small shape aspect ratio end, the poor local bond orientational order has the opposite effects on facilitating local rotational and translational hopping.

Shearing Janus Nanoparticles Confined in Two-Dimensional Space: Reshaped Cluster Configurations and Defined Assembling Kinetics

The self-assembly of anisotropic nanoparticles (ANPs) possesses a wide array of potential applications in various fields, ranging from nanotechnology to material science. Despite intense research of the thermodynamic self-assembly of ANPs, elucidating their nonequilibrium behaviors under confinement still remains an urgent issue. Here, by performing simulation and theoretical justification, we present for the first time a study of the shear-induced behaviors of Janus spheres (the most elementary ANPs) confined in two-dimensional space. Our results demonstrate that the collective effects of shear and bonding structures can give rise to reshaped cluster configurations, featured by the chiral transition of clusters. Scaling analysis and numerical modeling are performed to quantitatively capture the assembling kinetics of dispersed Janus spheres, thereby suggesting an exotic way to bridge the gap between anisotropic and isotropic particles. The findings highlight confinement and shearing engineering as a versatile strategy to tailor the superstructures formed by ANPs toward unique properties.

Nonclassical Nucleation in a Solid-Solid Transition of Confined Hard Spheres

A solid-solid phase transition of colloidal hard spheres confined between two planar hard walls is studied using a combination of molecular dynamics and Monte Carlo simulation. The transition from a solid consisting of five crystalline layers with square symmetry (5□) to a solid consisting of four layers with triangular symmetry (4△) is shown to occur through a nonclassical nucleation mechanism that involves the initial formation of a precritical liquid cluster, within which the cluster of the stable 4△ phase grows. Free-energy calculations show that the transition occurs in one step, crossing a single free-energy barrier, and that the critical nucleus consists of a small 4△ solid cluster wetted by a metastable liquid. In addition, the liquid cluster and the solid cluster are shown to grow at the planar hard walls. We also find that the critical nucleus size increases with supersaturation, which is at odds with classical nucleation theory. The △-solid-like cluster is shown to contain both face-centered-cubic and hexagonal-close-packed ordered particles.

Cooling the two-dimensional short spherocylinder liquid to the tetratic phase: Heterogeneous dynamics with one-way coupling between rotational and translational hopping

We numerically demonstrate the transition from the isotropic liquid to the tetratic phase with quasilong-range tetratic alignment order (i.e., with nearly parallel or perpendicular aligned rods), for the cold two-dimensional (2D) short spherocylinder system before crystallization and investigate the thermal assisted heterogeneous rotational and translational micromotions. Comparing with the 2D liquid of isotropic particles, spherocylinders introduce extra rotational degrees of freedom and destroy packing isotropy and the equivalence between rotational and translational motions. It is found that cooling leads to the stronger dynamical heterogeneity with more cooperative hopping and the stronger retardations of rotational hopping than translational hopping. Under topological constraints from nearly parallel and perpendicular rods of the tetratic phase, longitudinal and transverse translational hopping can occur without rotational hopping, but not the reverse. The empty space trailing a neighboring translational hopping patch is needed for triggering the patch rotational hopping with its translational motion into the empty space. It is the origin for the observed increasing separation of hopping time scales and the one-way coupling between rotational and translational hopping. Strips of longitudinally or transversely aligned rods can be ruptured and reconnected with neighboring strips through buckling, kink formation, and patch rotation, under the unbalanced torques or forces from their neighboring rods and thermal kicks.

Localized orientational order chaperones the nucleation of rotator phases in hard polyhedral particles

The nucleation kinetics of the rotator phase in hard cuboctahedra, truncated octahedra, and rhombic dodecahedra is simulated via a combination of forward flux sampling and umbrella sampling. For comparable degrees of supersaturation, the polyhedra are found to have significantly lower free-energy barriers and faster nucleation rates than hard spheres. This difference primarily stems from localized orientational ordering, which steers polyhedral particles to pack more efficiently. Orientational order hence fosters here the growth of orientationally disordered nuclei.

Dual translocation pathways in smectic liquid crystals facilitated by molecular flexibility

We investigate translocation mechanisms in smectic A liquid crystals (LCs) by a realistic, coarse-grained model of a LC compound comprising a stiff azobenzene core with flexible tails. We observe that the molecules can permeate from one smectic layer to the next via two different mechanisms, with and without significant reorientation, the former being facilitated through transverse interlayer intermediates. This is possible due to the intrinsic flexibility of the molecules. The two processes lead to characteristic signatures in the Van Hove self-correlation function, which can also be observed experimentally.

Liquid crystal order in colloidal suspensions of spheroidal particles by direct current electric field assembly

DC electric fields are used to produce colloidal assemblies with orientational and layered positional order from a dilute suspension of spheroidal particles. These 3D assemblies, which can be visualized in situ by confocal microscopy, are achieved in short time spans (t < 1 h) by the application of a constant voltage across the capacitor-like device. This method yields denser and more ordered assemblies than had been previously reported with other assembly methods. Structures with a high degree of orientational order as well as layered positional order normal to the electrode surface are observed. These colloidal structures are explained as a consequence of electrophoretic deposition and field-assisted assembly. The interplay between the deposition rate and the rotational Brownian motion is found to be critical for the optimal ordering, which occurs when these rates, as quantified by the Peclet number, are of order one. The results suggest that the mechanism leading to ordering is equilibrium self-assembly but with kinetics dramatically accelerated by the application of the DC electric field. Finally, the crystalline symmetry of the densest structure formed is determined and compared with previously studied spheroidal assemblies.

Frustration of the isotropic-columnar phase transition of colloidal hard platelets by a transient cubatic phase

Using simulations and theory, we show that the cubatic phase is metastable for three model hard platelets. The locally favored structures of perpendicular particle stacks in the fluid prevent the formation of the columnar phase through geometric frustration resulting in vitrification. Also, we find a direct link between structure and dynamic heterogeneities in the cooperative rotation of particle stacks, which is crucial for the devitrification process. Finally, we show that the lifetime of the glassy cubatic phase can be tuned by surprisingly small differences in particle shape.

Glass transitions in quasi-two-dimensional suspensions of colloidal ellipsoids

We observed a two-step glass transition in monolayers of colloidal ellipsoids by video microscopy. The glass transition in the rotational degree of freedom was at a lower density than that in the translational degree of freedom. Between the two transitions, ellipsoids formed an orientational glass. Approaching the respective glass transitions, the rotational and translational fastest-moving particles in the supercooled liquid moved cooperatively and formed clusters with power-law size distributions. The mean cluster sizes diverge in power law as they approach the glass transitions. The clusters of translational and rotational fastest-moving ellipsoids formed mainly within pseudonematic domains and around the domain boundaries, respectively.

Crystallization and aging in hard-sphere glasses

We report new results from our programme of molecular dynamics simulation of hard-sphere systems, focusing on crystallization and glass formation at high concentrations. First we consider a much larger system than hitherto, N = 86 400 equal-sized particles. The results are similar to those obtained with a smaller system, studied previously, showing conventional nucleation and growth of crystals at concentrations near melting and crossing over to a spinodal-like regime at higher concentrations where the free energy barrier to nucleation appears to be negligible. Second, we investigate the dependence on the initial state of the system. We have devised a Monte Carlo 'constrained aging' method to move the particles in such a way that crystallization is discouraged. After a period of such aging, the standard molecular dynamics programme is run. For a system of N = 3200, we find that constrained aging encourages caging of the particles and slows crystallization somewhat. Nevertheless, both aged and unaged systems crystallize at volume fraction φ = 0.61 whereas neither system shows full crystallization in the duration of the simulation at φ = 0.62, a concentration still significantly below that of random close packing.