Showing 1–50 of 99 results for author: Glotzer, S C

Change Point Detection of Events in Molecular Simulations using dupin

Authors: Brandon L. Butler, Domagoj Fijan, Sharon C. Glotzer

Abstract: Particle tracking is commonly used to study time-dependent behavior in many different types of physical and chemical systems involving constituents that span many length scales, including atoms, molecules, nanoparticles, granular particles, and even larger objects. Behaviors of interest studied using particle tracking information include thermodynamic phase transitions, structural transitions, pro… ▽ More Particle tracking is commonly used to study time-dependent behavior in many different types of physical and chemical systems involving constituents that span many length scales, including atoms, molecules, nanoparticles, granular particles, and even larger objects. Behaviors of interest studied using particle tracking information include thermodynamic phase transitions, structural transitions, protein folding, and crystallization. A common challenge in studies of these systems involves change detection. Change point detection discerns when a temporal signal undergoes a change in distribution. These changes can be local or global, instantaneous or prolonged, obvious or subtle. Moreover, system-wide changes marking an interesting physical or chemical phenomenon (e.g. crystallization of a liquid) are often preceded by events (e.g. pre-nucleation clusters) that are localized and can occur anywhere at anytime in the system. For these reasons, detecting events in particle trajectories generated by molecular simulation is challenging and typically accomplished via ad hoc solutions unique to the behavior and system under study. Consequently, methods for event detection lack generality, and those used in one field are not easily used by scientists in other fields. Here we present a new Python-based tool, dupin, that allows for universal event detection from particle trajectory data irrespective of the system details. dupin works by creating a signal representing the simulation and partitioning the signal based on events (changes within the trajectory). This approach allows for studies where manual annotating of event boundaries would require a prohibitive amount of time. Furthermore, dupin can serve as a tool in automated and reproducible workflows. We demonstrate the application of dupin using two examples and discuss its applicability to a wider class of problems. △ Less

Submitted 23 January, 2024; v1 submitted 22 December, 2023; originally announced December 2023.

Layer-by-Layer Assembled Nanowire Networks Enable Graph Theoretical Design of Multifunctional Coatings

Authors: Wenbing Wu, Alain Kadar, Sang Hyun Lee, Bum Chul Park, Jeffery E. Raymond, Thomas K. Tsotsis, Carlos E. S. Cesnik, Sharon C. Glotzer, Valerie Goss, Nicholas A. Kotov

Abstract: Multifunctional coatings are central for information, biomedical, transportation and energy technologies. These coatings must possess hard-to-attain properties and be scalable, adaptable, and sustainable, which makes layer-by-layer assembly (LBL) of nanomaterials uniquely suitable for these technologies. What remains largely unexplored is that LBL enables computational methodologies for structural… ▽ More Multifunctional coatings are central for information, biomedical, transportation and energy technologies. These coatings must possess hard-to-attain properties and be scalable, adaptable, and sustainable, which makes layer-by-layer assembly (LBL) of nanomaterials uniquely suitable for these technologies. What remains largely unexplored is that LBL enables computational methodologies for structural design of these composites. Utilizing silver nanowires (NWs), we develop and validate a graph theoretical (GT) description of their LBL composites. GT successfully describes the multilayer structure with nonrandom disorder and enables simultaneous rapid assessment of several properties of electrical conductivity, electromagnetic transparency, and anisotropy. GT models for property assessment can be rapidly validated due to (1) quasi-2D confinement of NWs and (2) accurate microscopy data for stochastic organization of the NW networks. We finally show that spray-assisted LBL offers direct translation of the GT-based design of composite coatings to additive, scalable manufacturing of drone wings with straightforward extensions to other technologies. △ Less

Submitted 30 October, 2023; v1 submitted 23 October, 2023; originally announced October 2023.

Complex motion of steerable vesicular robots filled with active colloidal rods

Authors: Sophie Y. Lee, Philipp W. A. Schönhöfer, Sharon C. Glotzer

Abstract: While the collective motion of active particles has been studied extensively, effective strategies to navigate particle swarms without external guidance remain elusive. We introduce a method to control the trajectories of two-dimensional swarms of active rod-like particles by confining the particles to rigid bounding membranes (vesicles) with non-uniform curvature. We show that the propelling agen… ▽ More While the collective motion of active particles has been studied extensively, effective strategies to navigate particle swarms without external guidance remain elusive. We introduce a method to control the trajectories of two-dimensional swarms of active rod-like particles by confining the particles to rigid bounding membranes (vesicles) with non-uniform curvature. We show that the propelling agents spontaneously form clusters at the membrane wall and collectively propel the vesicle, turning it into an active superstructure. To further guide the motion of the superstructure, we add discontinuous features to the rigid membrane boundary in the form of a kinked tip, which acts as a steering component to direct the motion of the vesicle. We report that the system's geometrical and material properties, such as the aspect ratio and Peclet number of the active rods as well as the kink angle and flexibility of the membrane, determine the stacking of active particles close to the kinked confinement and induce a diverse set of dynamical behaviors of the superstructure, including linear and circular motion both in the direction of, and opposite to, the kink. From a systematic study of these various behaviors, we design vesicles with switchable and reversible locomotions by tuning the confinement parameters. The observed phenomena suggest a promising mechanism for particle transportation and could be used as a basic element to navigate active matter through complex and tortuous environments. △ Less

Submitted 24 August, 2023; originally announced August 2023.

Comments: 11 pages, 5 figures

Quasicrystalline Nanocrystal Superlattice with Partial Matching Rules

Authors: Xingchen Ye, Jun Chen, M. Eric Irrgang, Michael Engel, Angang Dong, Sharon C. Glotzer, Christopher B. Murray

Abstract: Expanding the library of self-assembled superstructures provides insight into the behavior of atomic crystals and supports the development of materials with mesoscale order. Here we build upon recent findings of soft matter quasicrystals and report a quasicrystalline binary nanocrystal superlattice that exhibits correlations in the form of partial matching rules reducing tiling disorder. We determ… ▽ More Expanding the library of self-assembled superstructures provides insight into the behavior of atomic crystals and supports the development of materials with mesoscale order. Here we build upon recent findings of soft matter quasicrystals and report a quasicrystalline binary nanocrystal superlattice that exhibits correlations in the form of partial matching rules reducing tiling disorder. We determine a three-dimensional structure model through electron tomography and direct imaging of surface topography. The 12-fold rotational symmetry of the quasicrystal is broken in sub-layers, forming a random tiling of rectangles, large triangles, and small triangles with 6-fold symmetry. We analyze the geometry of the experimental tiling and discuss factors relevant for the stabilization of the quasicrystal. Our joint experimental-computational study demonstrates the power of nanocrystal superlattice engineering and further narrows the gap between the richness of crystal structures found with atoms and in soft matter assemblies. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 18 pages, 4 figures

Journal ref: Nature Materials 16, 214-219 (2017)

Rationalizing Euclidean Assemblies of Hard Polyhedra from Tessellations in Curved Space

Authors: Philipp W. A. Schönhöfer, Kai Sun, Xiaoming Mao, Sharon C. Glotzer

Abstract: Entropic self-assembly is governed by the shape of the constituent particles, yet a priori prediction of crystal structures from particle shape alone is non-trivial for anything but the simplest of space-filling shapes. At the same time, most polyhedra are not space-filling due to geometric constraints, but these constraints can be relaxed or even eliminated by sufficiently curving space. We show… ▽ More Entropic self-assembly is governed by the shape of the constituent particles, yet a priori prediction of crystal structures from particle shape alone is non-trivial for anything but the simplest of space-filling shapes. At the same time, most polyhedra are not space-filling due to geometric constraints, but these constraints can be relaxed or even eliminated by sufficiently curving space. We show using Monte Carlo simulations that the majority of hard Platonic shapes self-assemble entropically into space-filling crystals when constrained to the surface volume of a 3-sphere. As we gradually decrease curvature to "flatten" space and compare the local morphologies of crystals assembling in curved and flat space, we show that the Euclidean assemblies can be categorized either as remnants of tessellations in curved space (tetrahedra and dodecahedra) or non-tessellation-based assemblies caused by large-scale geometric frustration (octahedra and icosahedra). △ Less

Submitted 12 May, 2023; originally announced May 2023.

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

Authors: Philipp W. A. Schönhöfer, Sharon C. Glotzer

Abstract: 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… ▽ More 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. △ Less

Submitted 24 February, 2022; v1 submitted 19 February, 2022; originally announced February 2022.

Comments: 13 pages, 8 figures

Tuning Stoichiometry to Promote Formation of Binary Colloidal Superlattices

Authors: R. Allen LaCour, Timothy C. Moore, Sharon C. Glotzer

Abstract: 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 pro… ▽ More 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. △ Less

Submitted 11 December, 2021; originally announced December 2021.

A Route to Hierarchical Assembly of Colloidal Diamond

Authors: Yuan Zhou, Rose K. Cersonsky, Sharon C. Glotzer

Abstract: Photonic crystals, appealing for their ability to control light, are constructed by periodic regions of different dielectric constants. Yet, the structural holy grail in photonic materials, diamond, remains challenging to synthesize at the colloidal length scale. Here we explore new ways to assemble diamond using modified gyrobifastigial (mGBF) nanoparticles, a shape that resembles two antialigned… ▽ More Photonic crystals, appealing for their ability to control light, are constructed by periodic regions of different dielectric constants. Yet, the structural holy grail in photonic materials, diamond, remains challenging to synthesize at the colloidal length scale. Here we explore new ways to assemble diamond using modified gyrobifastigial (mGBF) nanoparticles, a shape that resembles two antialigned triangular prisms. We investigate the parameter space that leads to the self-assembly of diamond, and we compare the likelihood of defects in diamond self-assembled via mGBF vs. the nanoparticle shape that is the current focus for assembling diamond, the truncated tetrahedra. We introduce a potential route for realizing mGBF particles by dimerizing triangular prisms using attractive patches, and we report the impact of this superstructure on the photonic properties. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: Originally submitted to Soft Matter October 3, 2021

Shape-Driven Caging Dynamics of Hard Polygons

Authors: Vyas Ramasubramani, Thi Vo, Joshua A. Anderson, Sharon C. Glotzer

Abstract: The fundamentals of Brownian motion have been largely well understood since the early 20th century, with most recent additions focusing on understanding anomalous diffusion via rescaling of drag coefficents. That focus emphasizes long-time dynamic behavior, but recent results indicate that additional, secondary modes are also present at short and intermediate times in fluids of anisotropic particl… ▽ More The fundamentals of Brownian motion have been largely well understood since the early 20th century, with most recent additions focusing on understanding anomalous diffusion via rescaling of drag coefficents. That focus emphasizes long-time dynamic behavior, but recent results indicate that additional, secondary modes are also present at short and intermediate times in fluids of anisotropic particles. Here, we study the dynamics of a representative family of nearly-hard n-gons. Using molecular dynamics simulations, we study a distinct form of caging only present in anisotropic systems. We show that this caging behavior emerges in the mean squared displacement of n-gons at intermediate particle volume fractions. We then develop an extended Langevin theory directly coupling translational and rotational motion as a function of the relative anisotropy for different n-gons that predicts the observed caging behavior. Extending our theory to incorporate secondary, off-phase cross-correlations between particles further enables the prediction of both translational and rotational relaxation times of the system. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Microscopic Theory of Entropic Bonding for Colloidal Crystal Prediction

Authors: Thi Vo, Sharon C. Glotzer

Abstract: Entropy alone can self-assemble hard particles into colloidal crystals of remarkable complexity whose structures are the same as atomic and molecular crystals, but with larger lattice spacings. Although particle-based molecular simulation is a powerful tool for predicting self-assembly by exploring phase space, it is not yet possible to predict colloidal crystal structures a priori from particle s… ▽ More Entropy alone can self-assemble hard particles into colloidal crystals of remarkable complexity whose structures are the same as atomic and molecular crystals, but with larger lattice spacings. Although particle-based molecular simulation is a powerful tool for predicting self-assembly by exploring phase space, it is not yet possible to predict colloidal crystal structures a priori from particle shape as we can for atomic crystals based on electronic valency. Here we present such a first principles theory. By directly calculating and minimizing excluded volume within the framework of statistical mechanics, we describe the directional forces that emerge between hard shapes in familiar terms used to describe chemical bonds. Our theory predicts thermodynamically preferred structures for four families of hard polyhedra that match, in every instance, previous simulation results. The success of this first principles approach to entropic colloidal crystal structure prediction furthers fundamental understanding of both entropically driven crystallization and conceptual pictures of bonding in matter. △ Less

Submitted 5 July, 2021; originally announced July 2021.

Formation of a single quasicrystal upon collision of multiple grains

Authors: Insung Han, Kelly L. Wang, Andrew T. Cadotte, Zhucong Xi, Hadi Parsamehr, Xianghui Xiao, Sharon C. Glotzer, Ashwin J. Shahani

Abstract: Quasicrystals exhibit long-range order but lack translational symmetry. When grown as single crystals, they possess distinctive and unusual properties owing to the absence of grain boundaries. Unfortunately, conventional methods such as bulk crystal growth or thin film deposition only allow us to synthesize either polycrystalline quasicrystals or quasicrystals that are at most a few centimeters in… ▽ More Quasicrystals exhibit long-range order but lack translational symmetry. When grown as single crystals, they possess distinctive and unusual properties owing to the absence of grain boundaries. Unfortunately, conventional methods such as bulk crystal growth or thin film deposition only allow us to synthesize either polycrystalline quasicrystals or quasicrystals that are at most a few centimeters in size. Here, we reveal through real-time and 3D imaging the formation of a single decagonal quasicrystal arising from a hard collision between multiple growing quasicrystals in an Al-Co-Ni liquid. Through corresponding molecular dynamics simulations, we examine the underlying kinetics of quasicrystal coalescence and investigate the effects of initial misorientation between the growing quasicrystalline grains on the formation of grain boundaries. At small misorientation, coalescence occurs following rigid rotation that is facilitated by phasons. Our joint experimental-computational discovery paves the way toward fabrication of single, large-scale quasicrystals for novel applications. △ Less

Submitted 20 September, 2021; v1 submitted 26 June, 2021; originally announced June 2021.

Randomness in self-assembled colloidal crystals can widen photonic band gaps through particle shape and internal structure

Authors: Duanduan Wan, Sharon C. Glotzer

Abstract: Using computer simulations, we explore how thermal noise-induced randomness in a self-assembled photonic crystal affects its photonic band gaps (PBGs). We consider a two-dimensional photonic crystal comprised of a self-assembled array of parallel dielectric hard rods of infinite length with circular or square cross section. We find the PBGs can exist over a large range of intermediate packing dens… ▽ More Using computer simulations, we explore how thermal noise-induced randomness in a self-assembled photonic crystal affects its photonic band gaps (PBGs). We consider a two-dimensional photonic crystal comprised of a self-assembled array of parallel dielectric hard rods of infinite length with circular or square cross section. We find the PBGs can exist over a large range of intermediate packing densities. Counterintuitively, the largest band gap does not always appear at the packing density where the crystal is most ordered, despite the randomness inherent in any self-assembled structure. For rods with square cross section at intermediate packing densities, we find that the transverse magnetic (TM) band gap of the self-assembled (i.e. thermal) system can be larger than that of identical rods arranged in a perfect square lattice. By considering hollow rods, we find the band gap of transverse electric (TE) modes can be substantially increased while that of TM modes show no obvious improvement over solid rods. Our study suggests that particle shape and internal structure can be used to engineer the PBG of a self-assembled system despite the positional and orientational randomness arising from thermal noise. △ Less

Submitted 10 October, 2020; originally announced October 2020.

Comments: 5 pages, 5 figures

Journal ref: Phys. Rev. Lett. 126, 208002 (2021)

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

Authors: Bryan VanSaders, Sharon C. Glotzer

Abstract: 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 interst… ▽ More 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. △ Less

Submitted 29 December, 2019; originally announced December 2019.

Looped liquid-liquid coexistence in protein crystallization

Authors: Jens Glaser, Sharon C. Glotzer

Abstract: In view of the notorious complexity of protein--protein interactions, simplified models of proteins treated as patchy particles offer a promising strategy to obtain insight into the mechanism of crystallization. Here we report liquid--liquid phase separation (LLPS) with a highly asymmetric coexistence region in a computational model of rubredoxin with real molecular shape. The coexistence region t… ▽ More In view of the notorious complexity of protein--protein interactions, simplified models of proteins treated as patchy particles offer a promising strategy to obtain insight into the mechanism of crystallization. Here we report liquid--liquid phase separation (LLPS) with a highly asymmetric coexistence region in a computational model of rubredoxin with real molecular shape. The coexistence region terminates in both an upper (UCST) and a lower (LCST) critical solution temperature, and the complex molecular shape explains the closed-loop behavior of the LLPS. △ Less

Submitted 15 October, 2019; originally announced October 2019.

Comments: 6 pages, 4 figures, 1 supplementary PDF

Energy Landscapes for Digital Alchemy

Authors: John W. R. Morgan, Sharon C. Glotzer

Abstract: We apply energy landscape methods to digital alchemy, defining a system in which the parameters of the potential are treated as degrees of freedom. Using geometrical optimisation, we locate minima and transition states on the landscape for small clusters. We show that it is easy to find the parameters that give the lowest energy minimum, and that the distribution of minima on the alchemical landsc… ▽ More We apply energy landscape methods to digital alchemy, defining a system in which the parameters of the potential are treated as degrees of freedom. Using geometrical optimisation, we locate minima and transition states on the landscape for small clusters. We show that it is easy to find the parameters that give the lowest energy minimum, and that the distribution of minima on the alchemical landscape is concentrated in particular areas. We also conclude that the alchemical landscape is more frustrated, in terms of competition between low energy structures separated by high barriers. Transition states on the alchemical landscape are classified by whether they become minima or transition states when the potential parameters are fixed. Those that become minima have a significant alchemical component, while those that remain as transition states can be characterised mainly in terms of atomic displacements. △ Less

Submitted 25 June, 2019; originally announced June 2019.

Topological Order in Densely Packed Anisotropic Colloids

Authors: William Zygmunt, Erin G. Teich, Greg van Anders, Sharon C. Glotzer

Abstract: The existence of topological order is frequently associated with strongly coupled quantum matter. Here, we demonstrate the existence of topological phases in classical systems of densely packed, hard, anisotropic polyhedrally shaped colloidal particles. We show that previously reported transitions in dense packings lead to the existence of topologically ordered thermodynamic phases, which we show… ▽ More The existence of topological order is frequently associated with strongly coupled quantum matter. Here, we demonstrate the existence of topological phases in classical systems of densely packed, hard, anisotropic polyhedrally shaped colloidal particles. We show that previously reported transitions in dense packings lead to the existence of topologically ordered thermodynamic phases, which we show are stable away from the dense packing limit. Our work expands the library of known topological phases, whose experimental realization could provide new means for constructing plasmonic materials that are robust in the presence of fluctuations. △ Less

Submitted 24 June, 2019; originally announced June 2019.

Comments: 5 Pages, 3 Figures

Journal ref: Phys. Rev. E 100, 032608 (2019)

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

Authors: Vyas Ramasubramani, Bradley D. Dice, Eric S. Harper, Matthew P. Spellings, Joshua A. Anderson, Sharon C. Glotzer

Abstract: 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… ▽ More 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. △ Less

Submitted 27 March, 2020; v1 submitted 14 June, 2019; originally announced June 2019.

Comments: This paper is accepted for publication in Computer Physics Communications. For the source code on Github, see https://github.com/glotzerlab/freud . For the documentation on ReadTheDocs, see https://freud.readthedocs.io/en/stable

FCC-to-BCC phase transitions in convex and concave hard particle systems

Authors: Duanduan Wan, Chrisy Xiyu Du, Greg van Anders, Sharon C. Glotzer

Abstract: Particle shape plays an important role in the phase behavior of colloidal self-assembly. Recent progress in particle synthesis has made particles of polyhedral shapes and dimpled spherical shapes available. Here using computer simulations of hard particle models, we study face-centered cubic to body-centered cubic (FCC-to-BCC) phase transitions in a convex 432 polyhedral shape family and a concave… ▽ More Particle shape plays an important role in the phase behavior of colloidal self-assembly. Recent progress in particle synthesis has made particles of polyhedral shapes and dimpled spherical shapes available. Here using computer simulations of hard particle models, we study face-centered cubic to body-centered cubic (FCC-to-BCC) phase transitions in a convex 432 polyhedral shape family and a concave dimpled sphere family. Particles in both families have four-, three-, and two-fold rotational symmetries. Via free energy calculations we find the FCC-to-BCC transitions in both families are first order. As a previous work reports the FCC-to-BCC phase transition is first order in a convex 332 family of hard polyhedra, our work provides additional insight into the FCC-to-BCC transition and how the convexity or concavity of particle shape affects phase transition pathways. △ Less

Submitted 28 January, 2019; originally announced January 2019.

Comments: 6 pages, 5 figures

Journal ref: J. Phys. Chem. B 123, 9038 (2019)

Pressure-Tunable Photonic Band Gaps in an Entropic Colloidal Crystal

Authors: Rose K. Cersonsky, Julia Dshemuchadse, James A. Antonaglia, Greg van Anders, Sharon C. Glotzer

Abstract: Materials adopting the diamond structure possess useful properties in atomic and colloidal systems, and are a popular target for synthesis in colloids where a photonic band gap is possible. The desirable photonic properties of the diamond structure pose an interesting opportunity for reconfigurable matter: can we create a colloidal crystal able to switch reversibly to and from the diamond structur… ▽ More Materials adopting the diamond structure possess useful properties in atomic and colloidal systems, and are a popular target for synthesis in colloids where a photonic band gap is possible. The desirable photonic properties of the diamond structure pose an interesting opportunity for reconfigurable matter: can we create a colloidal crystal able to switch reversibly to and from the diamond structure? Drawing inspiration from high-pressure transitions of diamond-forming atomic systems, we design a system of polyhedrally-shaped particles that transitions from diamond to a tetragonal diamond derivative upon a small pressure change. The transition can alternatively be triggered by changing the shape of the particle in-situ. We propose that the transition provides a reversible reconfiguration process for a potential new colloidal material, and draw parallels between this transition and phase behavior of the atomic transitions from which we take inspiration. △ Less

Submitted 16 March, 2018; originally announced March 2018.

Comments: 6 REVTeX pages, 4 figures

Journal ref: Phys. Rev. Materials 2, 125201 (2018)

Mapping Disorder in Entropically Ordered Crystals

Authors: James A. Antonaglia, Greg van Anders, Sharon C. Glotzer

Abstract: Systems of hard shapes crystallize due to entropy. How is entropy distributed among translational and rotational microscopic contributions? We answer this question by decomposing thermal fluctuation of crystals of hard hexagons into collective modes, a generalization and quantification of the Onsager picture of hard rod liquid crystals. We show that at densities both near densest packing and near… ▽ More Systems of hard shapes crystallize due to entropy. How is entropy distributed among translational and rotational microscopic contributions? We answer this question by decomposing thermal fluctuation of crystals of hard hexagons into collective modes, a generalization and quantification of the Onsager picture of hard rod liquid crystals. We show that at densities both near densest packing and near the solid-hexatic melting transition, solids of hard regular hexagons hold most of their entropy in translational degrees of freedom. △ Less

Submitted 15 March, 2018; originally announced March 2018.

Comments: 5 REVTeX pages, 3 figures

Efficient Phase Diagram Sampling by Active Learning

Authors: Chengyu Dai, Isaac R. Bruss, Sharon C. Glotzer

Abstract: We address the problem of efficient phase diagram sampling by adopting active learning techniques from machine learning, and achieve an 80% reduction in the sample size (number of sampled statepoints) needed to establish the phase boundary up to a given precision in example application. Traditionally, data is collected on a uniform grid of predetermined statepoints. This approach, also known as gr… ▽ More We address the problem of efficient phase diagram sampling by adopting active learning techniques from machine learning, and achieve an 80% reduction in the sample size (number of sampled statepoints) needed to establish the phase boundary up to a given precision in example application. Traditionally, data is collected on a uniform grid of predetermined statepoints. This approach, also known as grid search in the machine learning community, suffers from low efficiency by sampling statepoints that provide no information about the phase boundaries. We propose an active learning approach to overcome this deficiency by adaptively choosing the next most informative statepoint(s) every round. This is done by interpolating the sampled statepoints' phases by Gaussian Process regression. An acquisition function quantifies the informativeness of possible next statepoints, maximizing the information content in each subsequently sampled statepoint. We also generalize our approach with state-of-the-art batch sampling techniques to better utilize parallel computing resources. We demonstrate the usefulness of our approach in a few example simulations relevant to soft matter physics, although our algorithms are general. Our active learning enhanced phase diagram sampling method greatly accelerates research and opens up opportunities for extra-large scale exploration of a wide range of phase diagrams by simulations or experiments. △ Less

Submitted 8 March, 2018; originally announced March 2018.

Strain Fields in Repulsive Colloidal Crystals

Authors: Bryan VanSaders, Julia Dshemuchadse, Sharon C. Glotzer

Abstract: The concept of a local linear elastic strain field is commonly used in the metallurgical research community to approximate the collective effect of atomic displacements around crystalline defects. Here we show that the elastic strain field approximation is a useful tool in colloidal systems. For colloidal crystals with repulsive particle interaction potentials, given similar mechanical properties,… ▽ More The concept of a local linear elastic strain field is commonly used in the metallurgical research community to approximate the collective effect of atomic displacements around crystalline defects. Here we show that the elastic strain field approximation is a useful tool in colloidal systems. For colloidal crystals with repulsive particle interaction potentials, given similar mechanical properties, sharper potentials lead to: 1) free energies of deformation dominated by entropy, 2) lower variance in strain field fluctuations, 3) increased tension-compression asymmetry near dislocation core regions, and 4) smaller windows of applicability of the linear elastic approximation. We show that the window of linear behavior for entropic colloidal crystals is broadened for pressures at which the inter-particle separation sufficiently exceeds the range of steep repulsive interactions. △ Less

Submitted 20 February, 2018; originally announced February 2018.

Journal ref: Phys. Rev. Materials 2, 063604 (2018)

Predicting colloidal crystals from shapes via inverse design and machine learning

Authors: Yina Geng, Greg van Anders, Sharon C. Glotzer

Abstract: A fundamental challenge in materials design is linking building block attributes to crystal structure. Addressing this challenge is particularly difficult for systems that exhibit emergent order, such as entropy-stabilized colloidal crystals. We combine recently developed techniques in inverse design with machine learning to construct a model that correctly classifies the crystals of more than ten… ▽ More A fundamental challenge in materials design is linking building block attributes to crystal structure. Addressing this challenge is particularly difficult for systems that exhibit emergent order, such as entropy-stabilized colloidal crystals. We combine recently developed techniques in inverse design with machine learning to construct a model that correctly classifies the crystals of more than ten thousand polyhedral shapes into 13 different structures with a predictive accuracy of 96% using only two geometric shape measures. With three measures, 98% accuracy is achieved. We test our model on previously reported colloidal crystal structures for 71 symmetric polyhedra and obtain 92% accuracy. Our findings (1) demonstrate that entropic colloidal crystals are controlled by surprisingly few parameters, (2) provide a quantitative model to predict these crystals solely from the geometry of their building blocks, and (3) suggest a prediction paradigm that easily generalizes to other self-assembled materials. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: 4 figures

Shapes within shapes: how particles arrange inside a cavity

Authors: Duanduan Wan, Sharon C. Glotzer

Abstract: We calculate the configurational entropy of hard particles confined in a cavity using Monte Carlo integration. Multiple combinations of particle and cavity shapes are considered. For small numbers of particles $N$, we show that the entropy decreases monotonically with increasing cavity aspect ratio, regardless of particle shape. As $N$ increases, we find ordered regions of high and low particle de… ▽ More We calculate the configurational entropy of hard particles confined in a cavity using Monte Carlo integration. Multiple combinations of particle and cavity shapes are considered. For small numbers of particles $N$, we show that the entropy decreases monotonically with increasing cavity aspect ratio, regardless of particle shape. As $N$ increases, we find ordered regions of high and low particle density, with the highest density near the boundary for all particle and cavity shape combinations. Our findings provide insights relevant to engineering particles in confined spaces, entropic barriers, and systems with depletion interactions. △ Less

Submitted 6 January, 2018; originally announced January 2018.

Comments: 6 pages, 8 figures

Journal ref: Soft Matter, 14, 3012 (2018)

Relevance of Packing to Colloidal Self-Assembly

Authors: Rose Cersonsky, Greg van Anders, Paul M. Dodd, Sharon C. Glotzer

Abstract: 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… ▽ More 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 comprised 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 in the infinite pressure limit are imperfect versions of the ideal packing shape. △ Less

Submitted 6 December, 2017; originally announced December 2017.

Comments: 6 pages

Journal ref: PNAS 115, 1439-1444 (2018)

Engineering Entropy for the Inverse Design of Colloidal Crystals from Hard Shapes

Authors: Yina Geng, Greg van Anders, Paul M. Dodd, Julia Dshemuchadse, Sharon C. Glotzer

Abstract: Throughout the physical sciences, entropy stands out as a pivotal but enigmatic concept that, in materials design, often takes a backseat to energy. Here, we demonstrate how to precisely engineer entropy to achieve desired colloidal crystals. We demonstrate the inverse design of hard particles that assemble six different target colloidal crystals due solely to entropy maximization. Our approach ef… ▽ More Throughout the physical sciences, entropy stands out as a pivotal but enigmatic concept that, in materials design, often takes a backseat to energy. Here, we demonstrate how to precisely engineer entropy to achieve desired colloidal crystals. We demonstrate the inverse design of hard particles that assemble six different target colloidal crystals due solely to entropy maximization. Our approach efficiently samples $10^8$ particle shapes from 88- and 192-dimensional design spaces to discover thermodynamically optimal shapes. We design particle shapes that self assemble known crystals with optimized thermodynamic stability, as well as new crystal structures with no known atomic or other equivalent. △ Less

Submitted 6 December, 2017; originally announced December 2017.

Comments: 6 pages

Journal ref: Science Advances 5, eaaw0514 (2019)

Phase separation of self-propelled ballistic particles

Authors: Isaac R. Bruss, Sharon C. Glotzer

Abstract: Self-propelled particles phase separate into coexisting dense and dilute regions above a critical density. The statistical nature of their stochastic motion lends itself to various theories that predict the onset of phase separation. However, these theories are ill equipped to describe such behavior when noise become negligible. To overcome this limitation, we present a predictive model that relie… ▽ More Self-propelled particles phase separate into coexisting dense and dilute regions above a critical density. The statistical nature of their stochastic motion lends itself to various theories that predict the onset of phase separation. However, these theories are ill equipped to describe such behavior when noise become negligible. To overcome this limitation, we present a predictive model that relies on two density-dependent timescales: $τ_F$, the mean time particles spend between collisions; and $τ_C$, the mean lifetime of a collision. We show that only when $τ_F < τ_C$ do collisions last long enough to develop a growing cluster and initiate phase separation. Using both analytical calculations and active particle simulations, we measure these timescales and determine the critical density for phase separation in both 2D and 3D. △ Less

Submitted 23 April, 2018; v1 submitted 5 December, 2017; originally announced December 2017.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. E 97, 042609 (2018)

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

Authors: Carl S. Adorf, James Antonaglia, Julia Dshemuchadse, Sharon C. Glotzer

Abstract: 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… ▽ More 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. △ Less

Submitted 21 November, 2018; v1 submitted 11 November, 2017; originally announced November 2017.

Comments: 11 pages, 6 figures

Journal ref: J. Chem. Phys. 149, 204102 (2018)

Machine learning for crystal identification and discovery

Authors: Matthew Spellings, Sharon C Glotzer

Abstract: As computers get faster, researchers -- not hardware or algorithms -- become the bottleneck in scientific discovery. Computational study of colloidal self-assembly is one area that is keenly affected: even after computers generate massive amounts of raw data, performing an exhaustive search to determine what (if any) ordered structures occur in a large parameter space of many simulations can be ex… ▽ More As computers get faster, researchers -- not hardware or algorithms -- become the bottleneck in scientific discovery. Computational study of colloidal self-assembly is one area that is keenly affected: even after computers generate massive amounts of raw data, performing an exhaustive search to determine what (if any) ordered structures occur in a large parameter space of many simulations can be excruciating. We demonstrate how machine learning can be applied to discover interesting areas of parameter space in colloidal self assembly. We create numerical fingerprints -- inspired by bond orientational order diagrams -- of structures found in self-assembly studies and use these descriptors to both find interesting regions in a phase diagram and identify characteristic local environments in simulations in an automated manner for simple and complex crystal structures. Utilizing these methods allows analysis methods to keep up with the data generation ability of modern high-throughput computing environments. △ Less

Submitted 26 March, 2018; v1 submitted 26 October, 2017; originally announced October 2017.

Comments: Fixed typo, added missing acknowledgment, added supplementary information

Non-Close-Packed Three-Dimensional Quasicrystals

Authors: Pablo F. Damasceno, Sharon C. Glotzer, Michael Engel

Abstract: Quasicrystals are frequently encountered in condensed matter. They are important candidates for equilibrium phases from the atomic scale to the nanoscale. Here, we investigate the computational self-assembly of four quasicrystals in a single model system of identical particles interacting with a tunable isotropic pair potential. We reproduce a known icosahedral quasicrystal and report a decagonal… ▽ More Quasicrystals are frequently encountered in condensed matter. They are important candidates for equilibrium phases from the atomic scale to the nanoscale. Here, we investigate the computational self-assembly of four quasicrystals in a single model system of identical particles interacting with a tunable isotropic pair potential. We reproduce a known icosahedral quasicrystal and report a decagonal quasicrystal, a dodecagonal quasicrystal, and an octagonal quasicrystal. The quasicrystals have low coordination number or occur in systems with mesoscale density variations. We also report a network gel phase. △ Less

Submitted 23 June, 2017; v1 submitted 28 February, 2017; originally announced February 2017.

Comments: 9 pages, 8 figures

Journal ref: J. Phys. Cond. Mat. 29, 234005 (2017)

Shape and symmetry determine two-dimensional melting transitions of hard regular polygons

Authors: Joshua A. Anderson, James Antonaglia, Jaime A. Millan, Michael Engel, Sharon C. Glotzer

Abstract: The melting transition of two-dimensional (2D) systems is a fundamental problem in condensed matter and statistical physics that has advanced significantly through the application of computational resources and algorithms. 2D systems present the opportunity for novel phases and phase transition scenarios not observed in 3D systems, but these phases depend sensitively on the system and thus predict… ▽ More The melting transition of two-dimensional (2D) systems is a fundamental problem in condensed matter and statistical physics that has advanced significantly through the application of computational resources and algorithms. 2D systems present the opportunity for novel phases and phase transition scenarios not observed in 3D systems, but these phases depend sensitively on the system and thus predicting how any given 2D system will behave remains a challenge. Here we report a comprehensive simulation study of the phase behavior near the melting transition of all hard regular polygons with $3\leq n\leq 14$ vertices using massively parallel Monte Carlo simulations of up to one million particles. By investigating this family of shapes, we show that the melting transition depends upon both particle shape and symmetry considerations, which together can predict which of three different melting scenarios will occur for a given $n$. We show that systems of polygons with as few as seven edges behave like hard disks; they melt continuously from a solid to a hexatic fluid and then undergo a first-order transition from the hexatic phase to the fluid phase. We show that this behavior, which holds for all $7\leq n\leq 14$, arises from weak entropic forces among the particles. Strong directional entropic forces align polygons with fewer than seven edges and impose local order in the fluid. These forces can enhance or suppress the discontinuous character of the transition depending on whether the local order in the fluid is compatible with the local order in the solid. As a result, systems of triangles, squares, and hexagons exhibit a KTHNY-type continuous transition between fluid and hexatic, tetratic, and hexatic phases, respectively, and a continuous transition from the appropriate "x"-atic to the solid. [abstract truncated due to arxiv length limitations]. △ Less

Submitted 22 December, 2016; v1 submitted 2 June, 2016; originally announced June 2016.

Journal ref: Phys. Rev. X 7, 021001 (2017)

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

Authors: Andrew T. Cadotte, Julia Dshemuchadse, Pablo F. Damasceno, Richmond S. Newman, Sharon C. Glotzer

Abstract: 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… ▽ More 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. △ Less

Submitted 23 May, 2016; originally announced May 2016.

Comments: 6 pages, 3 figures

Journal ref: Soft Matter 12, 7073-7078 (2016)

Shape Driven Solid--Solid Transitions in Colloids

Authors: Chrisy Xiyu Du, Greg van Anders, Richmond S. Newman, Sharon C. Glotzer

Abstract: Despite the fundamental importance of solid--solid transitions for metallurgy, ceramics, earth science, reconfigurable materials, and colloidal matter, the details of how materials transform between two solid structures are poorly understood. We introduce a class of {simple} model systems in which the direct control of local order, via colloid shape change, induces solid--solid phase transitions,… ▽ More Despite the fundamental importance of solid--solid transitions for metallurgy, ceramics, earth science, reconfigurable materials, and colloidal matter, the details of how materials transform between two solid structures are poorly understood. We introduce a class of {simple} model systems in which the direct control of local order, via colloid shape change, induces solid--solid phase transitions, and characterize how the transitions happen thermodynamically. We find that within a single shape family there are solid--solid transitions that can occur with or without a thermal activation barrier. Our results provide new means for the study of solid--solid phase transitions, and have implications for designing reconfigurable materials. △ Less

Submitted 12 March, 2017; v1 submitted 2 March, 2016; originally announced March 2016.

Comments: 8 REVTeX pages, 7 figures

Journal ref: PNAS 114, E3892--E3899 (2017)

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

Authors: Luis S. Froufe-Pérez, Michael Engel, Pablo F. Damasceno, Nicolas Muller, Jakub Haberko, Sharon C. Glotzer, Frank Scheffold

Abstract: 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 c… ▽ More 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. △ Less

Submitted 19 June, 2016; v1 submitted 2 February, 2016; originally announced February 2016.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. Lett. 117, 053902 (2016)

Using Depletion to Control Colloidal Crystal Assemblies of Hard Cuboctahedra

Authors: Andrew S. Karas, Jens Glaser, Sharon C. Glotzer

Abstract: 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 th… ▽ More 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. △ Less

Submitted 11 March, 2016; v1 submitted 14 October, 2015; originally announced October 2015.

Scalable Metropolis Monte Carlo for simulation of hard shapes

Authors: Joshua A. Anderson, M. Eric Irrgang, Sharon C. Glotzer

Abstract: 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.… ▽ More 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. △ Less

Submitted 15 September, 2015; originally announced September 2015.

Comments: 5 figures

A parallel algorithm for implicit depletant simulations

Authors: Jens Glaser, Andrew S. Karas, Sharon C. Glotzer

Abstract: 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 surroundin… ▽ More 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. △ Less

Submitted 27 August, 2015; originally announced August 2015.

Comments: 10 pages, 8 figures

Coarsening dynamics of binary liquids with active rotation

Authors: Syeda Sabrina, Matthew Spellings, Sharon C. Glotzer, Kyle J. M. Bishop

Abstract: 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 leve… ▽ More 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. △ Less

Submitted 23 July, 2015; originally announced July 2015.

Digital Alchemy for Materials Design: Colloids and Beyond

Authors: Greg van Anders, Daphne Klotsa, Andrew S. Karas, Paul M. Dodd, Sharon C. Glotzer

Abstract: 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 b… ▽ More 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. △ Less

Submitted 2 March, 2016; v1 submitted 17 July, 2015; originally announced July 2015.

Comments: 17 REVTeX pages, title fixed to match journal version

Journal ref: ACS Nano 9, 9542-9553 (2015)

Controlling Chirality of Entropic Crystals

Authors: Pablo F. Damasceno, Andrew S. Karas, Benjamin A. Schultz, Michael Engel, Sharon C. Glotzer

Abstract: 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… ▽ More 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. △ Less

Submitted 13 July, 2015; originally announced July 2015.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. Lett. 115, 158303 (2015)

Strong Orientational Coordinates and Orientational Order Parameters For Symmetric Objects

Authors: Amir Haji-Akbari, Sharon C. Glotzer

Abstract: 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… ▽ More 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. △ Less

Submitted 29 October, 2015; v1 submitted 8 July, 2015; originally announced July 2015.

Comments: 28 pages, 10 figures, 5 tables

Shape Allophiles Improve Entropic Assembly

Authors: Eric S. Harper, Ryan L. Marson, Joshua A. Anderson, Greg van Anders, Sharon C. Glotzer

Abstract: We investigate a class of "shape allophiles" that fit together like puzzle pieces as a method to access and stabilize desired structures by controlling directional entropic forces. Squares are cut into rectangular halves, which are shaped in an allophilic manner with the goal of re-assembling the squares while self-assembling the square lattice. We examine the assembly characteristics of this syst… ▽ More We investigate a class of "shape allophiles" that fit together like puzzle pieces as a method to access and stabilize desired structures by controlling directional entropic forces. Squares are cut into rectangular halves, which are shaped in an allophilic manner with the goal of re-assembling the squares while self-assembling the square lattice. We examine the assembly characteristics of this system via the potential of mean force and torque, and the fraction of particles that entropically bind. We generalize our findings and apply them to self-assemble triangles into a square lattice via allophilic shaping. Through these studies we show how shape allophiles can be useful in assembling and stabilizing desired phases with appropriate allophilic design. △ Less

Submitted 25 June, 2015; originally announced June 2015.

Comments: 17 pages, 7 figures. Published as E. S. Harper, R. L. Marson, J. A. Anderson, G. van Anders, and S. C. Glotzer, Soft Matter, 2015, DOI: 10.1039/c5sm01351h

Journal ref: Soft Matter 11, 7250-7256 (2015)

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

Authors: N. Khalid Ahmed, Greg van Anders, Elizabeth R. Chen, Sharon C. Glotzer

Abstract: 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… ▽ More 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. △ Less

Submitted 4 October, 2016; v1 submitted 13 January, 2015; originally announced January 2015.

Comments: V2: 18 pages, 22 figures. Revised text with clarifications, corrections, and additional supporting data

Complexity in surfaces of densest packings for families of polyhedra

Authors: Elizabeth R. Chen, Daphne Klotsa, Michael Engel, Pablo F. Damasceno, Sharon C. Glotzer

Abstract: 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 analyt… ▽ More 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. △ Less

Submitted 10 September, 2013; originally announced September 2013.

Comments: 16 pages, 8 figures

Journal ref: Phys. Rev. X 4 (2014) 011024

Understanding shape entropy through local dense packing

Authors: Greg van Anders, Daphne Klotsa, N. Khalid Ahmed, Michael Engel, Sharon C. Glotzer

Abstract: 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 directio… ▽ More 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. △ Less

Submitted 28 October, 2014; v1 submitted 4 September, 2013; originally announced September 2013.

Comments: v3: 9+10 revtex pages, 10 figures, minor changes, changed title to match journal version

Journal ref: PNAS 111, E4812-E4821 (2014)

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

Authors: Nguyen H. P. Nguyen, Daphne Klotsa, Michael Engel, Sharon C. Glotzer

Abstract: 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 increasin… ▽ More 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. △ Less

Submitted 9 August, 2013; originally announced August 2013.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. Lett. 112 (2014) 075701

Entropically Patchy Particles: Engineering Valence through Shape Entropy

Authors: Greg van Anders, N. Khalid Ahmed, Ross Smith, Michael Engel, Sharon C. Glotzer

Abstract: 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… ▽ More 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. △ Less

Submitted 13 February, 2014; v1 submitted 28 April, 2013; originally announced April 2013.

Comments: v1: revtex, 6 pages, 5 figures. v2: completely rewritten, revtex, 11 pages, 7 figures. Figures are higher resolution than journal version

Journal ref: ACS Nano 8, 931-940 (2014)

Packing and Self-assembly of Truncated Triangular Bipyramids

Authors: Amir Haji-Akbari, Elizabeth R. Chen, Michael Engel, Sharon C Glotzer

Abstract: Motivated by breakthroughs in the synthesis of faceted nano- and colloidal particles, as well as theoretical and computational studies of their packings, we investigate a family of truncated triangular bipyramids. We report dense periodic packings with small unit cells that were obtained via numerical and analytical optimization. The maximal packing fraction $φ_{\max}$ changes continuously with th… ▽ More Motivated by breakthroughs in the synthesis of faceted nano- and colloidal particles, as well as theoretical and computational studies of their packings, we investigate a family of truncated triangular bipyramids. We report dense periodic packings with small unit cells that were obtained via numerical and analytical optimization. The maximal packing fraction $φ_{\max}$ changes continuously with the truncation parameter $t$. Eight distinct packings are identified based on discontinuities in the first and second derivatives of $φ_{\max}(t)$. These packings differ in the number of particles in the fundamental domain (unit cell) and the type of contacts between the particles. In particular, we report two packings with four particles in the unit cell for which both $φ_{\max}(t)$ and $φ'_{\max}(t)$ are continuous and the discontinuity occurs in the second derivative only. In the self-assembly simulations that we perform for larger boxes with 2048 particles, only one out of eight packings is found to assemble. In addition, the degenerate quasicrystal reported previously for triangular bipyramids without truncation [{http://prl.aps.org/abstract/PRL/v107/i21/e215702}{Haji-Akbari \emph{et al.}, \emph{Phys. Rev. Lett.} \textbf{107}: 215702 (2011)}] assembles for truncations as high as 0.45. The self-assembly propensities for the structures formed in the thermodynamic limit are explained using the isoperimetric quotient of the particles and the coordination number in the disordered fluid and in the assembled structure. △ Less

Submitted 14 August, 2013; v1 submitted 10 April, 2013; originally announced April 2013.

Comments: 14 pages, 12 figures

Journal ref: Phys. Rev. E 88, 012127 (2013)

Massively parallel Monte Carlo for many-particle simulations on GPUs

Authors: Joshua A. Anderson, Eric Jankowski, Thomas L. Grubb, Michael Engel, Sharon C. Glotzer

Abstract: 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 i… ▽ More 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. △ Less

Submitted 23 August, 2013; v1 submitted 7 November, 2012; originally announced November 2012.

Comments: 11 pages, 4 figures, 4 tables

Journal ref: Journal of Computational Physics 254, 27-38 (2013)

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

Authors: Michael Engel, Joshua A. Anderson, Sharon C. Glotzer, Masaharu Isobe, Etienne P. Bernard, Werner Krauth

Abstract: 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… ▽ More 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. △ Less

Submitted 23 April, 2013; v1 submitted 7 November, 2012; originally announced November 2012.

Comments: 9 pages, 8 figures, 2 tables

Journal ref: Phys. Rev. E 87, 042134 (2013)