BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Sharon C. Glotzer

Sharon C. Glotzer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.mtrl-scicond-mat.stat-mechphysics.comp-phcond-mat.mes-hallmath.MGphysics.chem-phBiological PhysicsComputational Engineering, Finance, and ScienceComputational GeometryComputer Visioncond-mat.dis-nncond-mat.othermath.OCMathematical Softwarenlin.AOphysics.optics

Catalog footprint

What is connected

44works
17topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2026arXiv

Entropic Colloidal Crystal Prediction: A Quantum Density Functional Theory Inspired Approach

In pursuit of a colloidal analogue to quantum density functional theory (DFT) predictions of atomic crystal structures, we report a new, classical DFT that predicts the relative thermodynamic stability of colloidal crystals of hard, convex particle shapes. In contrast to standard classical DFT approaches, our theory maps the hard particle system to an auxiliary system in which we treat the particles as fixed "nuclei" embedded in a fictitious, spatially varying density field that distributes throughout the auxiliary system. By minimizing the free energy of the auxiliary system, and through comparison with known equations of state and free energy calculations using thermodynamic integration, we show that the auxiliary system with the lowest free energy corresponds to the most probable crystal of hard shapes in the original system.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Using Particle Shape to Control Defects in Colloidal Crystals on Spherical Interfaces

Spherical particles confined to a sphere surface cannot pack densely into a hexagonal lattice without defects. In this study, we use hard particle Monte Carlo simulations to determine the effects of continuously deformable shape anisotropy and underlying crystal lattice preference on inevitable defect structures and their distribution within colloidal assemblies of hard rounded polyhedra confined to a closed sphere surface. We demonstrate that cube particles form a simple square assembly, overcoming lattice/topology incompatibility, and maximize entropy by distributing eight three-fold defects evenly on the sphere. By varying particle shape smoothly from cubes to spheres we reveal how the distribution of defects changes from square antiprismatic to icosahedral symmetry. Congruent studies of rounded tetrahedra reveal additional varieties of characteristic defect patterns within three, four, and six-fold symmetric lattices. This work has promising implications for programmable defect generation to facilitate different vesicle buckling modes using colloidal particle emulsions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Curvature-Controlled Geometrical Lensing Behavior in Self-Propelled Colloidal Particle Systems

In many biological systems, the curvature of the surfaces cells live on influence 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) that 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 dynamic phases, including cyclic MIPS on sphercylinders, wave-like MIPS on spherocones, and cluster fluctuations on metaballs.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Tuning Stoichiometry to Promote Formation of Binary Colloidal Superlattices

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

freud: A Software Suite for High Throughput Analysis of Particle Simulation Data

The freud Python package is a powerful library for analyzing simulation data. Written with modern simulation and data analysis workflows in mind, freud provides a Python interface to fast, parallelized C++ routines that run efficiently on laptops, workstations, and supercomputing clusters. The package provides the core tools for finding particle neighbors in periodic systems, and offers a uniform API to a wide variety of methods implemented using these tools. As such, freud users can access standard methods such as the radial distribution function as well as newer, more specialized methods such as the potential of mean force and torque and local crystal environment analysis with equal ease. While many comparable tools place a heavy emphasis on reading and operating on trajectory file formats, freud instead accepts numerical arrays of data directly as inputs. By remaining agnostic to its data source, freud is suitable for analyzing any coarse-grained particle simulation, regardless of the original data representation or simulation method. When used for on-the-fly analysis in conjunction with scriptable simulation software such as HOOMD-blue, freud enables smart simulations that adapt to the current state of the system, allowing users to study phenomena such as nucleation and growth.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Pinning Dislocations in Colloidal Crystals with Active Particles that Seek Stacking Faults

There is growing interest in functional, adaptive devices built from colloidal subunits of micron size or smaller. A colloidal material with dynamic mechanical properties could facilitate such microrobotic machines. Here we study via computer simulation how active interstitial particles in small quantities can be used to modify the bulk mechanical properties of a colloidal crystal. Passive interstitial particles are known to pin dislocations in metals, thereby increasing resistance to plastic deformation. We extend this tactic by employing anisotropic active interstitials that travel super-diffusively and bind strongly to stacking faults associated with partial dislocations. We find that: 1) interstitials that are effective at reducing plasticity compromise between strong binding to stacking faults and high mobility in the crystal bulk. 2) Reorientation of active interstitials in the crystal depends upon rotational transitions between high-symmetry crystal directions. 3) The addition of certain active interstitial shapes at concentrations as low as $60$ per million host particles ($0.006\%$) can create a shear threshold for dislocation migration.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Inverse Design of Simple Pair Potentials for the Self-Assembly of Complex Structures

The synthesis of complex materials through the self-assembly of particles at the nanoscale provides opportunities for the realization of novel material properties. However, the inverse design process to create experimentally feasible interparticle interaction strategies is uniquely challenging. Standard methods for the optimization of isotropic pair potentials tend toward overfitting, resulting in solutions with too many features and length scales that are challenging to map to mechanistic models. Here we introduce a method for the optimization of simple pair potentials that minimizes the relative entropy of the complex target structure while directly considering only those length scales most relevant for self-assembly. Our approach maximizes the relative information of a target pair distribution function with respect to an ansatz distribution function via an iterative update process. During this process, we filter high frequencies from the Fourier spectrum of the pair potential, resulting in interaction potentials that are smoother and simpler in real space, and therefore likely easier to make. We show that pair potentials obtained by this method assemble their target structure more robustly with respect to optimization method parameters than potentials optimized without filtering.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Crossover Behavior in the Packing and Assembly of Multivalent Lock-and-Key Colloids

Emergent behaviors occur in a vast array of systems across many scales, and are of fundamental physical importance because of the intrinsic difficulty in linking microscopic system properties to macroscopic behaviors. Here we study the emergent self-assembly behavior of model systems of recently synthesized families of concave dimpled hard spheres, or lock-and-key colloids. We find that as dimple size increases each family exhibits a crossover from a structure that does not reflect the particle symmetry to one that does and, surprisingly, the point at which this crossover occurs is approximately independent of the particle symmetry. Using a combination of numerical and analytic techniques we study systems at infinite and finite pressure, and find different common control parameters in each limit. Our results suggest there exists a set of experimentally realizable colloidal systems that exhibit complex emergent behaviors that can be traced to a common underlying microscopic control parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Digital Alchemy for Materials Design: Colloids and Beyond

Starting with the early alchemists, a holy grail of science has been to make desired materials by modifying the attributes of basic building blocks. Building blocks that show promise for assembling new complex materials can be synthesized at the nanoscale with attributes that would astonish the ancient alchemists in their versatility. However, this versatility means that making direct connection between building block attributes and bulk behavior is both necessary for rationally engineering materials, and difficult because building block attributes can be altered in many ways. Here we show how to exploit the malleability of the valence of colloidal nanoparticle "elements" to directly and quantitatively link building block attributes to bulk behavior through a statistical thermodynamic framework we term "digital alchemy". We use this framework to optimize building blocks for a given target structure, and to determine which building block attributes are most important to control for self assembly, through a set of novel thermodynamic response functions, moduli and susceptibilities. We thereby establish direct links between the attributes of colloidal building blocks and the bulk structures they form. Moreover, our results give concrete solutions to the more general conceptual challenge of optimizing emergent behaviors in nature, and can be applied to other types of matter. As examples, we apply digital alchemy to systems of truncated tetrahedra, rhombic dodecahedra, and isotropically interacting spheres that self assemble diamond, FCC, and icosahedral quasicrystal structures, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Self-assembly of a space-tessellating structure in the binary system of hard tetrahedra and octahedra

We report the formation of a binary crystal of hard polyhedra due solely to entropic forces. Although the alternating arrangement of octahedra and tetrahedra is a known space-tessellation, it had not previously been observed in self-assembly simulations. Both known one-component phases - the dodecagonal quasicrystal of tetrahedra and the densest-packing of octahedra in the Minkowski lattice - are found to coexist with the binary phase. No additional crystalline phases were observed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Role of Short-Range Order and Hyperuniformity in the Formation of Band Gaps in Disordered Photonic Materials

We study photonic band gap formation in two-dimensional high refractive index disordered materials where the dielectric structure is derived from packing disks in real and reciprocal space. Numerical calculations of the photonic density of states demonstrate the presence of a band gap for all polarizations in both cases. We find that the band gap width is controlled by the increase in positional correlation inducing short-range order and hyperuniformity concurrently. Our findings suggest that the optimization of short-range order, in particular the tailoring of Bragg scattering at the isotropic Brillouin zone, are of key importance for designing disordered PBG materials.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Using Depletion to Control Colloidal Crystal Assemblies of Hard Cuboctahedra

Depletion interactions arise from entropic forces, and their ability to induce aggregation and even ordering of colloidal particles through self-assembly is well established, especially for spherical colloids. We vary the size and concentration of penetrable hard sphere depletants in a system of cuboctahedra, and we show how depletion changes the preferential facet alignment of the colloids and thereby selects different crystal structures. Moreover, we explain the cuboctahedra phase behavior using perturbative free energy calculations. We find that cuboctahedra can form a stable simple cubic phase, and, remarkably, that the stability of this phase can only be rationalized by considering the effects of both the colloid and depletant entropy. We corroborate our results by analyzing how the depletant concentration and size affect the emergent directional entropic forces and hence the effective particle shape. We propose the use of depletants as a means of easily changing the effective shape of self-assembling anisotropic colloids.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A parallel algorithm for implicit depletant simulations

We present an algorithm to simulate the many-body depletion interaction between anisotropic colloids in an implicit way, integrating out the degrees of freedom of the depletants, which we treat as an ideal gas. Because the depletant particles are statistically independent and the depletion interaction is short-ranged, depletants are randomly inserted in parallel into the excluded volume surrounding a single translated and/or rotated colloid. A configurational bias scheme is used to enhance the acceptance rate. The method is validated and benchmarked both on multi-core CPUs and graphics processing units (GPUs) for the case of hard spheres, hemispheres and discoids. With depletants, we report novel cluster phases, in which hemispheres first assemble into spheres, which then form ordered hcp/fcc lattices. The method is significantly faster than any method without cluster moves and that tracks depletants explicitly, for systems of colloid packing fraction $ϕ_c<0.50$, and additionally enables simulation of the fluid-solid transition.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coarsening dynamics of binary liquids with active rotation

Active matter comprised of many self-driven units can exhibit emergent collective behaviors such as pattern formation and phase separation in both biologica and synthetic systems. While these behaviors are increasingly well understood for ensembles of linearly self-propelled particles, less is known about the collective behaviors of active rotating particles where energy input at the particle level gives rise to rotational particle motion. A recent simulation study revealed that active rotation can induce phase separation in mixtures of counter-rotating particles in 2D. In contrast to that of linearly self-propelled particles, the phase separation of counter-rotating fluids is accompanied by steady convective flows that originate at the fluid-fluid interface. Here, we investigate the influence of these flows on the coarsening dynamics of actively rotating binary liquids using a phenomenological, hydrodynamic model that combines a Cahn-Hilliard equation for the fluid composition with a Navier-Stokes equation for the fluid velocity. The effect of active rotation is introduced though an additional force within the Navier-Stokes equations that arises due to gradients in the concentrations of clockwise and counter-clockwise rotating particles. Depending on the strength of active rotation and that of frictional interactions with the stationary surroundings, we observe and explain new dynamical behaviors such as "active coarsening" via self-generated flows as well as the emergence of self-propelled vortex doublets. We confirm that many of the qualitative behaviors identified by the continuum model can also be found in discrete, particle-based simulations of actively rotating liquids. Our results highlight further opportunities for achieving complex dissipative structures in active materials subject to distributed actuation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Controlling Chirality of Entropic Crystals

Colloidal crystal structures with complexity and diversity rivaling atomic and molecular crystals have been predicted and obtained for hard particles by entropy maximization. However, so far homochiral colloidal crystals, which are candidates for photonic metamaterials, are absent. Using Monte Carlo simulations we show that chiral polyhedra exhibiting weak directional entropic forces self-assemble either an achiral crystal or a chiral crystal with limited control over the crystal handedness. Building blocks with stronger faceting exhibit higher selectivity and assemble a chiral crystal with handedness uniquely determined by the particle chirality. Tuning the strength of directional entropic forces by means of particle rounding or the use of depletants allows for reconfiguration between achiral and homochiral crystals. We rationalize our findings by quantifying the chirality strength of each particle, both from particle geometry and potential of mean force and torque diagrams.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Scalable Metropolis Monte Carlo for simulation of hard shapes

We design and implement HPMC, a scalable hard particle Monte Carlo simulation toolkit, and release it open source as part of HOOMD-blue. HPMC runs in parallel on many CPUs and many GPUs using domain decomposition. We employ BVH trees instead of cell lists on the CPU for fast performance, especially with large particle size disparity, and optimize inner loops with SIMD vector intrinsics on the CPU. Our GPU kernel proposes many trial moves in parallel on a checkerboard and uses a block-level queue to redistribute work among threads and avoid divergence. HPMC supports a wide variety of shape classes, including spheres / disks, unions of spheres, convex polygons, convex spheropolygons, concave polygons, ellipsoids / ellipses, convex polyhedra, convex spheropolyhedra, spheres cut by planes, and concave polyhedra. NVT and NPT ensembles can be run in 2D or 3D triclinic boxes. Additional integration schemes permit Frenkel-Ladd free energy computations and implicit depletant simulations. In a benchmark system of a fluid of 4096 pentagons, HPMC performs 10 million sweeps in 10 minutes on 96 CPU cores on XSEDE Comet. The same simulation would take 7.6 hours in serial. HPMC also scales to large system sizes, and the same benchmark with 16.8 million particles runs in 1.4 hours on 2048 GPUs on OLCF Titan.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Strong Orientational Coordinates and Orientational Order Parameters For Symmetric Objects

Recent advancements in the synthesis of anisotropic macromolecules and nanoparticles have spurred an immense interest in theoretical and computational studies of self-assembly. The cornerstone of such studies is the role of shape in self-assembly and in inducing complex order. The problem of identifying different types of order that can emerge in such systems can, however, be challenging. Here, we revisit the problem of quantifying orientational order in systems of building blocks with non-trivial rotational symmetries. We first propose a systematic way of constructing orientational coordinates for such symmetric building blocks. We call the arising tensorial coordinates strong orientational coordinates (SOCs) as they fully and exclusively specify the orientation of a symmetric object. We then use SOCs to describe and quantify local and global orientational order, and spatiotemporal orientational correlations in systems of symmetric building blocks. The SOCs and the orientational order parameters developed in this work are not only useful in performing and analyzing computer simulations of symmetric molecules or particles, but can also be utilized for the efficient storage of rotational information in long trajectories of evolving many-body systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Entropically Patchy Particles: Engineering Valence through Shape Entropy

Patchy particles are a popular paradigm for the design and synthesis of nanoparticles and colloids for self-assembly. In "traditional" patchy particles, anisotropic interactions arising from patterned coatings, functionalized molecules, DNA, and other enthalpic means create the possibility for directional binding of particles into higher-ordered structures. Although the anisotropic geometry of non-spherical particles contributes to the interaction patchiness through van der Waals, electrostatic, and other interactions, how particle shape contributes entropically to self-assembly is only now beginning to be understood. It has been recently demonstrated that, for hard shapes, entropic forces are directional. A newly proposed theoretical framework that defines and quantifies directional entropic forces demonstrates the anisotropic--that is, patchy--nature of these emergent, attractive forces. Here we introduce the notion of entropically patchy particles as the entropic counterpart to enthalpically patchy particles. Using three example "families" of shapes, we judiciously modify entropic patchiness by introducing geometric features to the particles so as to target specific crystal structures, which then assembled with Monte Carlo simulations. We quantify the emergent entropic valence via a potential of mean force and torque. We generalize these shape operations to shape anisotropy dimensions, in analogy with the anisotropy dimensions introduced for enthalpically patchy particles. Our findings demonstrate that entropic patchiness and emergent valence provide a way of engineering directional bonding into nanoparticle systems, whether in the presence or absence of additional, non-entropic forces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Excitations are localized and relaxation is hierarchical in glass-forming liquids

For several atomistic models of glass formers, at conditions below their glassy dynamics onset temperatures, ${T_\mathrm{o}}$, we use importance sampling of trajectory space to study the structure, statistics and dynamics of excitations responsible for structural relaxation. Excitations are detected in terms of persistent particle displacements of length $a$. At supercooled conditions, for $a$ of the order of or smaller than a particle diameter, we find that excitations are associated with correlated particle motions that are sparse and localized, occupying a volume with an average radius that is temperature independent and no larger than a few particle diameters. We show that the statistics and dynamics of these excitations are facilitated and hierarchical. Excitation energy scales grow logarithmically with $a$. Excitations at one point in space facilitate the birth and death of excitations at neighboring locations, and space-time excitation structures are microcosms of heterogeneous dynamics at larger scales. This nature of dynamics becomes increasingly dominant as temperature $T$ is lowered. We show that slowing of dynamics upon decreasing temperature below $T_\mathrm{o}$ is the result of a decreasing concentration of excitations and concomitant growing hierarchical length scales, and further that the structural relaxation time $τ$ follows the parabolic law, $\log(τ/ τ_\mathrm{o}) = J^2(1/T - 1/T_\mathrm{o})^2$, for $T<T_\mathrm{o}$, where $J$, $τ_\mathrm{o}$ and $T_\mathrm{o}$ can be predicted quantitatively from dynamics at short time scales. Particle motion is facilitated and directional, and we show this becomes more apparent with decreasing $T$. We show that stringlike motion is a natural consequence of facilitated, hierarchical dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Strong scaling of general-purpose molecular dynamics simulations on GPUs

We describe a highly optimized implementation of MPI domain decomposition in a GPU-enabled, general-purpose molecular dynamics code, HOOMD-blue (Anderson and Glotzer, arXiv:1308.5587). Our approach is inspired by a traditional CPU-based code, LAMMPS (Plimpton, J. Comp. Phys. 117, 1995), but is implemented within a code that was designed for execution on GPUs from the start (Anderson et al., J. Comp. Phys. 227, 2008). The software supports short-ranged pair force and bond force fields and achieves optimal GPU performance using an autotuning algorithm. We are able to demonstrate equivalent or superior scaling on up to 3,375 GPUs in Lennard-Jones and dissipative particle dynamics (DPD) simulations of up to 108 million particles. GPUDirect RDMA capabilities in recent GPU generations provide better performance in full double precision calculations. For a representative polymer physics application, HOOMD-blue 1.0 provides an effective GPU vs. CPU node speed-up of 12.5x.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Understanding shape entropy through local dense packing

Entropy drives the phase behavior of colloids ranging from dense suspensions of hard spheres or rods to dilute suspensions of hard spheres and depletants. Entropic ordering of anisotropic shapes into complex crystals, liquid crystals, and even quasicrystals has been demonstrated recently in computer simulations and experiments. The ordering of shapes appears to arise from the emergence of directional entropic forces (DEFs) that align neighboring particles, but these forces have been neither rigorously defined nor quantified in generic systems. Here, we show quantitatively that shape drives the phase behavior of systems of anisotropic particles upon crowding through DEFs. We define DEFs in generic systems, and compute them for several hard particle systems. We show that they are on the order of a few kT at the onset of ordering, placing DEFs on par with traditional depletion, van der Waals, and other intrinsic interactions. In experimental systems with these other interactions, we provide direct quantitative evidence that entropic effects of shape also contribute to self-assembly. We use DEFs to draw a distinction between self-assembly and packing behavior. We show that the mechanism that generates directional entropic forces is the maximization of entropy by optimizing local particle packing. We show that this mechanism occurs in a wide class of systems, and we treat, in a unified way, the entropy-driven phase behavior of arbitrary shapes incorporating the well-known works of Kirkwood, Onsager, and Asakura and Oosawa.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Complexity in surfaces of densest packings for families of polyhedra

Packings of hard polyhedra have been studied for centuries due to their mathematical aesthetic and more recently for their applications in fields such as nanoscience, granular and colloidal matter, and biology. In all these fields, particle shape is important for structure and properties, especially upon crowding. Here, we explore packing as a function of shape. By combining simulations and analytic calculations, we study three 2-parameter families of hard polyhedra and report an extensive and systematic analysis of the densest packings of more than 55,000 convex shapes. The three families have the symmetries of triangle groups (icosahedral, octahedral, tetrahedral) and interpolate between various symmetric solids (Platonic, Archimedean, Catalan). We find that optimal (maximum) packing density surfaces that reveal unexpected richness and complexity, containing as many as 130 different structures within a single family. Our results demonstrate the utility of thinking of shape not as a static property of an object in the context of packings, but rather as but one point in a higher dimensional shape space whose neighbors in that space may have identical or markedly different packings. Finally, we present and interpret our packing results in a consistent and generally applicable way by proposing a method to distinguish regions of packings and classify types of transitions between them.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Emergent collective phenomena in a mixture of hard shapes through active rotation

We investigate collective phenomena with rotationally driven spinners of concave shape. Each spinner experiences a constant internal torque in either a clockwise or counterclockwise direction. Although the spinners are modeled as hard, otherwise non-interacting rigid bodies, we find that their active motion induces an effective interaction that favors rotation in the same direction. With increasing density and activity, phase separation occurs via spinodal decomposition, as well as self-organization into rotating crystals. We observe the emergence of cooperative, super-diffusive motion along interfaces, which can transport inactive test particles. Our results demonstrate novel phase behavior of actively rotated particles that is not possible with linear propulsion or in non-driven, equilibrium systems of identical hard particles.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hard-disk equation of state: First-order liquid-hexatic transition in two dimensions with three simulation methods

We report large-scale computer simulations of the hard-disk system at high densities in the region of the melting transition. Our simulations reproduce the equation of state, previously obtained using the event-chain Monte Carlo algorithm, with a massively parallel implementation of the local Monte Carlo method and with event-driven molecular dynamics. We analyze the relative performance of these simulation methods to sample configuration space and approach equilibrium. Our results confirm the first-order nature of the melting phase transition in hard disks. Phase coexistence is visualized for individual configurations via the orientational order parameter field. The analysis of positional order confirms the existence of the hexatic phase.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Massively parallel Monte Carlo for many-particle simulations on GPUs

Current trends in parallel processors call for the design of efficient massively parallel algorithms for scientific computing. Parallel algorithms for Monte Carlo simulations of thermodynamic ensembles of particles have received little attention because of the inherent serial nature of the statistical sampling. In this paper, we present a massively parallel method that obeys detailed balance and implement it for a system of hard disks on the GPU. We reproduce results of serial high-precision Monte Carlo runs to verify the method. This is a good test case because the hard disk equation of state over the range where the liquid transforms into the solid is particularly sensitive to small deviations away from the balance conditions. On a Tesla K20, our GPU implementation executes over one billion trial moves per second, which is 148 times faster than on a single Intel Xeon E5540 CPU core, enables 27 times better performance per dollar, and cuts energy usage by a factor of 13. With this improved performance we are able to calculate the equation of state for systems of up to one million hard disks. These large system sizes are required in order to probe the nature of the melting transition, which has been debated for the last forty years. In this paper we present the details of our computational method, and discuss the thermodynamics of hard disks separately in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Design of Patchy Particles using Ternary Self-Assembled Monolayers

Recent simulations have studied the formation of patterns in a binary mixture of immiscible surfactants absorbed onto the surface of a spherical nanoparticle. The resulting patterns (Janus, spots and stripes) were in good agreement with experimental results. We perform dissipative particle dynamics (DPD) simulations to study the patterns obtained by adding a third surfactant to the monolayer as a guide towards increasing the richness and diversity of patchy particles synthesized this way. We predict a variety of new patterns that can be produced through different combinations of simple design elements, like nanocolloid size, degree of surfactant immiscibility, stoichiometry of the monolayer, and length difference between surfactants. In all cases, free energy minimization through conformational entropy maximization determines equilibrium pattern formation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimal Filling of Shapes

We present filling as a type of spatial subdivision problem similar to covering and packing. Filling addresses the optimal placement of overlapping objects lying entirely inside an arbitrary shape so as to cover the most interior volume. In n-dimensional space, if the objects are polydisperse n-balls, we show that solutions correspond to sets of maximal n-balls. For polygons, we provide a heuristic for finding solutions of maximal discs. We consider the properties of ideal distributions of N discs as N approaches infinity. We note an analogy with energy landscapes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimal Fillings - A new spatial subdivision problem related to packing and covering

We present filling as a new type of spatial subdivision problem that is related to covering and packing. Filling addresses the optimal placement of overlapping objects lying entirely inside an arbitrary shape so as to cover the most interior volume. In n-dimensional space, if the objects are polydisperse n-balls, we show that solutions correspond to sets of maximal n-balls and the solution space can reduced to the medial axis of a shape. We examine the structure of the solution space in two dimensions. For the filling of polygons, we provide detailed descriptions of a heuristic and a genetic algorithm for finding solutions of maximal discs. We also consider the properties of ideal distributions of N discs in polygons as N approaches infinity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Self Assembled Clusters of Spheres Related to Spherical Codes

We consider the thermodynamically driven self-assembly of spheres onto the surface of a central sphere. This assembly process forms self-limiting, or terminal, anisotropic clusters (N-clusters) with well defined structures. We use Brownian dynamics to model the assembly of N-clusters varying in size from two to twelve outer spheres, and free energy calculations to predict the expected cluster sizes and shapes as a function of temperature and inner particle diameter. We show that the arrangements of outer spheres at finite temperatures are related to spherical codes, an ideal mathematical sequence of points corresponding to densest possible sphere packings. We demonstrate that temperature and the ratio of the diameters of the inner and outer spheres dictate cluster morphology and dynamics. We find that some N-clusters exhibit collective particle rearrangements, and these collective modes are unique to a given cluster size N. We present a surprising result for the equilibrium structure of a 5-cluster, which prefers an asymmetric square pyramid arrangement over a more symmetric arrangement. Our results suggest a promising way to assemble anisotropic building blocks from constituent colloidal spheres.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Self-Assembly of Patchy Particles into Diamond Structures through Molecular Mimicry

Fabrication of diamond structures by self-assembly is a fundamental challenge in making three-dimensional photonic crystals. We simulate a system of model hard particles with attractive patches and show that they can self-assemble into a diamond structure from an initially disordered state. We quantify the extent to which the formation of the diamond structure can be facilitated by "seeding" the system with small diamond crystallites or by introducing a rotation interaction to mimic a carbon-carbon antibonding interaction. Our results suggest patchy particles may serve as colloidal "atoms" and "molecules" for the bottom-up self-assembly of three-dimensional crystals.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Structural diversity and the role of particle shape and dense fluid behavior in assemblies of hard polyhedra

A fundamental characteristic of matter is its ability to form ordered structures under the right thermodynamic conditions. Predicting these structures - and their properties - from the attributes of a material's building blocks is the holy grail of materials science. Here, we investigate the self-assembly of 145 hard convex polyhedra whose thermodynamic behavior arises solely from their anisotropic shape. Our results extend previous works on entropy-driven crystallization by demonstrating a remarkably high propensity for self-assembly and an unprecedented structural diversity, including some of the most complex crystalline phases yet observed in a non-atomic system. In addition to 22 Bravais and non-Bravais crystals, we report 66 plastic crystals (both Bravais and topologically close-packed), 21 liquid crystals (nematic, smectic, and columnar), and 44 glasses. We show that from simple measures of particle shape and local order in the disordered fluid, the class of ordered structure can be predicted.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Crystalline Assemblies and Densest Packings of a Family of Truncated Tetrahedra and the Role of Directional Entropic Forces

Polyhedra and their arrangements have intrigued humankind since the ancient Greeks and are today important motifs in condensed matter, with application to many classes of liquids and solids. Yet, little is known about the thermodynamically stable phases of polyhedrally-shaped building blocks, such as faceted nanoparticles and colloids. Although hard particles are known to organize due to entropy alone, and some unusual phases are reported in the literature, the role of entropic forces in connection with polyhedral shape is not well understood. Here, we study thermodynamic self-assembly of a family of truncated tetrahedra and report several atomic crystal isostructures, including diamond, β-tin, and high- pressure lithium, as the polyhedron shape varies from tetrahedral to octahedral. We compare our findings with the densest packings of the truncated tetrahedron family obtained by numerical compression and report a new space filling polyhedron, which has been overlooked in previous searches. Interestingly, the self-assembled structures differ from the densest packings. We show that the self-assembled crystal structures can be understood as a tendency for polyhedra to maximize face-to-face alignment, which can be generalized as directional entropic forces.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Facetted patchy particles through entropy-driven patterning of mixed ligand SAMS

We present a microscopic theory that describes the ordering of two distinct ligands on the surface of a faceted nanoparticle. The theory predicts that when one type of ligand is significantly bulkier than all others, the larger ligands preferentially align themselves along the edges and vertices of the nanoparticle. Monte Carlo simulations confirm these predictions. We show that the intrinsic conformational entropy of the ligands stabilizes this novel edge-aligned phase.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase Diagram of Hard Tetrahedra

Advancements in the synthesis of faceted nanoparticles and colloids have spurred interest in the phase behavior of polyhedral shapes. Regular tetrahedra have attracted particular attention because they prefer local symmetries that are incompatible with periodicity. Two dense phases of regular tetrahedra have been reported recently. The densest known tetrahedron packing is achieved in a crystal of triangular bipyramids (dimers) with packing density 4000/4671=85.63%. In simulation a dodecagonal quasicrystal is observed; its approximant, with periodic tiling (3.4.3^2.4), can be compressed to a packing fraction of 85.03%. Here, we show that the quasicrystal approximant is more stable than the dimer crystal for packing densities below 84% using Monte Carlo computer simulations and free energy calculations. To carry out the free energy calculations, we use a variation of the Frenkel-Ladd method for anisotropic shapes and thermodynamic integration. The enhanced stability of the approximant can be attributed to a network substructure, which maximizes the free volume (and hence the 'wiggle room') available to the particles and facilitates correlated motion of particles, which further contributes to entropy and leads to diffusion for packing densities below 65%. The existence of a solid-solid transition between structurally distinct phases not related by symmetry breaking -- the approximant and the dimer crystal-- is unusual for hard particle systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Self Assembly of Soft Matter Quasicrystals and Their Approximants

The surprising recent discoveries of quasicrystals and their approximants in soft matter systems poses the intriguing possibility that these structures can be realized in a broad range of nano- and micro-scale assemblies. It has been theorized that soft matter quasicrystals and approximants are largely entropically stabilized, but the thermodynamic mechanism underlying their formation remains elusive. Here, we use computer simulation and free energy calculations to demonstrate a simple design heuristic for assembling quasicrystals and approximants in soft matter systems. Our study builds on previous simulation studies of the self-assembly of dodecagonal quasicrystals and approximants in minimal systems of spherical particles with complex, highly-specific interaction potentials. We demonstrate an alternative entropy-based approach for assembling dodecagonal quasicrystals and approximants based solely on particle functionalization and shape, thereby recasting the interaction-potential-based assembly strategy in terms of simpler-to-achieve bonded and excluded-volume interactions. Here, spherical building blocks are functionalized with mobile surface entities to encourage the formation of structures with low surface contact area, including non-close-packed and polytetrahedral structures. The building blocks also possess shape polydispersity, where a subset of the building blocks deviate from the ideal spherical shape, discouraging the formation of close-packed crystals. We show that three different model systems with both of these features -- mobile surface entities and shape polydispersity -- consistently assemble quasicrystals and/or approximants. We argue that this design strategy can be widely exploited to assemble quasicrystals and approximants on the nano- and micro- scales. In addition, our results further elucidate the formation of soft matter quasicrystals in experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermal and Athermal Swarms of Self-Propelled Particles

Swarms of self-propelled particles exhibit complex behavior that can arise from simple models, with large changes in swarm behavior resulting from small changes in model parameters. We investigate the steady-state swarms formed by self-propelled Morse particles in three dimensions using molecular dynamics simulations optimized for GPUs. We find a variety of swarms of different overall shape assemble spontaneously and that for certain Morse potential parameters coexisting structures are observed. We report a rich "phase diagram" of athermal swarm structures observed across a broad range of interaction parameters. Unlike the structures formed in equilibrium self-assembly, we find that the probability of forming a self-propelled swarm can be biased by the choice of initial conditions. We investigate how thermal noise influences swarm formation and demonstrate ways it can be exploited to reconfigure one swarm into another. Our findings validate and extend previous observations of self-propelled Morse swarms and highlight open questions for predictive theories of nonequilibrium self-assembly.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dense crystalline dimer packings of regular tetrahedra

We present the densest known packing of regular tetrahedra with density phi = 4000/4671 = 0.856347... Like the recently discovered packings of Kallus et al. [arXiv:0910.5226] and Torquato-Jiao [arXiv:0912.4210], our packing is crystalline with a unit cell of four tetrahedra forming two triangular dipyramids (dimer clusters). We show that our packing has maximal density within a three-parameter family of dimer packings. Numerical compressions starting from random configurations suggest that the packing may be optimal at least for small cells with up to 16 tetrahedra and periodic boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Disordered, Quasicrystalline and Crystalline Phases of Densely Packed Tetrahedra

All hard, convex shapes are conjectured by Ulam to pack more densely than spheres, which have a maximum packing fraction of ϕ = π/\sqrt18 ~ 0.7405. For many shapes, simple lattice packings easily surpass this packing fraction. For regular tetrahedra, this conjecture was shown to be true only very recently; an ordered arrangement was obtained via geometric construction with ϕ = 0.7786, which was subsequently compressed numerically to ϕ = 0.7820. Here we show that tetrahedra pack much better than this, and in a completely unexpected way. Following a conceptually different approach, using thermodynamic computer simulations that allow the system to evolve naturally towards high-density states, we observe that a fluid of hard tetrahedra undergoes a first-order phase transition to a dodecagonal quasicrystal, which can be compressed to a packing fraction of ϕ = 0.8324. By compressing a crystalline approximant of the quasicrystal, the highest packing fraction we obtain is ϕ = 0.8503. If quasicrystal formation is suppressed, the system remains disordered, jams, and compresses to ϕ = 0.7858. Jamming and crystallization are both preceded by an entropy-driven transition from a simple fluid of independent tetrahedra to a complex fluid characterized by tetrahedra arranged in densely packed local motifs that form a percolating network at the transition. The quasicrystal that we report represents the first example of a quasicrystal formed from hard or non-spherical particles. Our results demonstrate that particle shape and entropy can produce highly complex, ordered structures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Harmonic Order Parameters for Characterizing Complex Particle Morphologies

Order parameters based on spherical harmonics and Fourier coefficients already play a significant role in condensed matter research in the context of systems of spherical or point particles. Here, we extend these types of order parameter to more complex shapes, such as those encountered in nanoscale self-assembly applications. To do so, we build on a powerful set of techniques that originate in the computer science field of "shape matching." We demonstrate how shape matching techniques can be applied to identify unknown structures and create highly-specialized \textit{ad hoc} order parameters. Additionally, we investigate the special symmetry properties of harmonic descriptors, and demonstrate how they can be exploited to provide optimal solutions to certain classes of problems. Our techniques can be applied to particle systems in general, both simulated and experimental, provided the particle positions are known.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

How Do Quasicrystals Grow?

Using molecular simulations, we show that the aperiodic growth of quasicrystals is controlled by the ability of the growing quasicrystal `nucleus' to incorporate kinetically trapped atoms into the solid phase with minimal rearrangement. In the system under investigation, which forms a dodecagonal quasicrystal, we show that this process occurs through the assimilation of stable icosahedral clusters by the growing quasicrystal. Our results demonstrate how local atomic interactions give rise to the long-range aperiodicity of quasicrystals.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Measurement of Growing Dynamical Lengthscales and Prediction of the Jamming Transition in a Granular Material

Supercooled liquids and dense colloids exhibit anomalous behaviour known as "spatially heterogeneous dynamics" (SHD), which becomes increasingly pronounced with approach to the glass transition. Recently, SHD has been observed in confined granular packings under slow shear near the onset of jamming, bolstering speculation that the two transitions are related. Here, we report measurements of SHD in a system of air-driven granular beads, as a function of both density and effective temperature. On approach to jamming, the dynamics become progressively slower and more spatially heterogeneous. The rapid growth of dynamical time and length scales characterizing the heterogeneities can be described both by mode-coupling theory and the Vogel-Tammann-Fulcher (VTF) equation, in analogy with glass-forming liquids. The value of the control variable at the VTF transition coincides with point-J, the random close-packed jamming density at which all motion ceases, indicating analogy with a zero temperature ideal glass transition. Our findings demonstrate further universality of the jamming concept and provide a significant step forward in the quest for a unified theory of jamming in disparate systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Pair Interaction Potentials of Colloids by Extrapolation of Confocal Microscopy Measurements of Collective Structure

A method for measuring the pair interaction potential between colloidal particles by extrapolation measurement of collective structure to infinite dilution is presented and explored using simulation and experiment. The method is particularly well suited to systems in which the colloid is fluorescent and refractive index matched with the solvent. The method involves characterizing the potential of mean force between colloidal particles in suspension by measurement of the radial distribution function using 3D direct visualization. The potentials of mean force are extrapolated to infinite dilution to yield an estimate of the pair interaction potential, $U(r)$. We use Monte Carlo (MC) simulation to test and establish our methodology as well as to explore the effects of polydispersity on the accuracy. We use poly-12-hydroxystearic acid-stabilized poly(methyl methacrylate) (PHSA-PMMA) particles dispersed in the solvent dioctyl phthalate (DOP) to test the method and assess its accuracy for three different repulsive systems for which the range has been manipulated by addition of electrolyte.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Crossover to Potential Energy Landscape Dominated Dynamics in a Model Glass-forming Liquid

An equilibrated model glass-forming liquid is studied by mapping successive configurations produced by molecular dynamics simulation onto a time series of inherent structures (local minima in the potential energy). Using this ``inherent dynamics'' approach we find direct numerical evidence for the long held view that below a crossover temperature, $T_x$, the liquid's dynamics can be separated into (i) vibrations around inherent structures and (ii) transitions between inherent structures (M. Goldstein, J. Chem. Phys. {\bf 51}, 3728 (1969)), i.e., the dynamics become ``dominated'' by the potential energy landscape. In agreement with previous proposals, we find that $T_x$ is within the vicinity of the mode-coupling critical temperature $T_c$. We further find that at the lowest temperature simulated (close to $T_x$), transitions between inherent structures involve cooperative, string like rearrangements of groups of particles moving distances substantially smaller than the average interparticle distance.

0 citations0 reviewsNo code linkedOpen access