Computational self-assembly of a one-component icosahedral quasicrystal

Plastic Crystal-to-Crystal Transition of Janus Particles under Shear

Colloidal Janus spheres in the bulk typically spontaneously assemble into plastic crystalline phases, while particle orientations exhibit glasslike dynamics without long-range order. Through Brownian dynamics simulations, we demonstrate that shear can trigger a phase transition from an isotropic crystal with orientational disorder to an orientationally ordered crystal with lamellae along the shear direction. This nonequilibrium transition is accompanied with the orientational ordering following a nucleation and growth mechanism. By performing a phenomenological extension of free energy analysis, we reveal that the nucleation originates from the orientation fluctuations induced by shear. The growth of the orientationally crystalline cluster is examined to be disklike, captured by developing a lattice model with memoryless state functions. These findings bring new insights into the mechanisms for the ordering transition of anisotropic particles at nonequilibrium states.

Complex Crystals from Size-Disperse Spheres

Colloids are rarely perfectly uniform but follow a distribution of sizes, shapes, and charges. This dispersity can be inherent (static) or develop and change over time (dynamic). Despite a long history of research, the conditions under which nonuniform particles crystallize and which crystal forms is still not well understood. Here, we demonstrate that hard spheres with Gaussian radius distribution and dispersity up to 19% always crystallize if compressed slowly enough, and they do so in surprisingly complex ways. This result is obtained by accelerating event-driven simulations with particle swap moves for static dispersity and particle resize moves for dynamic dispersity. Above 6% dispersity, AB_{2} Laves, AB_{13}, and a region of Frank-Kasper phases are found. The Frank-Kasper region includes a quasicrystal approximant with Pearson symbol oS276. Our findings are relevant for ordering phenomena in soft matter and alloys.

Static and dynamic structures of liquid Ba8Ga16Sn30: a melt of the thermoelectric clathrate compounds

X-ray diffraction and inelastic x-ray scattering measurements of liquid Ba₈Ga₁₆Sn₃₀ have been carried out to investigate local structure and atomic dynamics in the liquid. The pair distribution function shows shorter and longer interatomic distances in the first coordination shell. The dynamic structure factor exhibits the inelastic excitations on both sides of the quasielastic central peak. The inelastic excitations disperse with increasing the momentum transfer, suggesting the longitudinal acoustic mode. We found a low energy excitation in addition to the longitudinal acoustic excitation in the dynamic structure factor and it reminds us a strong relationship with a rattling motion of a guest (Ba) atom in the solid state. The temperature dependence of the pair distribution function and the longitudinal acoustic excitation energy is very weak in a range from 600 to 900 °C. The result suggests that Ba and other atoms in the melt are located around minimum positions of the effective pair potential approximated as a harmonic one.

Growth of two-dimensional dodecagonal colloidal quasicrystals: Particles with isotropic pair interactions with two length scales vs. patchy colloids with preferred binding angles

We explore the growth of colloidal quasicrystals with dodecagonal symmetry in two dimensions by employing Brownian dynamics simulations. On the one hand, we study the growth behavior of structures obtained in a system of particles that interact according to an isotropic pair potential with two typical length scales. On the other hand, we consider patchy colloids that possess only one typical interaction length scale but prefer given binding angles. In case of the isotropic particles, we show that an imbalance in the competition between the two distances might lead to defects with wrong nearest-neighbor distances in the resulting structure. In contrast, during the growth of quasicrystals with patchy colloids such defects do not occur due to the lack of a second interaction length scale. However, as a downside, the diffusion of patchy particles along a surface typically is slower such that domains occur where the particles possess different phononic and phasonic offsets. Our results are important to understand how soft matter quasicrystals can be grown as perfectly as possible and to obtain a deeper insight into the mechanisms of the growth of quasicrystals in general.

Multiple-scale structures: from Faraday waves to soft-matter quasicrystals

For many years, quasicrystals were observed only as solid-state metallic alloys, yet current research is now actively exploring their formation in a variety of soft materials, including systems of macromolecules, nanoparticles and colloids. Much effort is being invested in understanding the thermodynamic properties of these soft-matter quasicrystals in order to predict and possibly control the structures that form, and hopefully to shed light on the broader yet unresolved general questions of quasicrystal formation and stability. Moreover, the ability to control the self-assembly of soft quasicrystals may contribute to the development of novel photonics or other applications based on self-assembled metamaterials. Here a path is followed, leading to quantitative stability predictions, that starts with a model developed two decades ago to treat the formation of multiple-scale quasiperiodic Faraday waves (standing wave patterns in vibrating fluid surfaces) and which was later mapped onto systems of soft particles, interacting via multiple-scale pair potentials. The article reviews, and substantially expands, the quantitative predictions of these models, while correcting a few discrepancies in earlier calculations, and presents new analytical methods for treating the models. In so doing, a number of new stable quasicrystalline structures are found with octagonal, octadecagonal and higher-order symmetries, some of which may, it is hoped, be observed in future experiments.

Micrometer-sized TPM emulsion droplets with surface-mobile binding groups

Colloids coated with lipid membranes have been widely employed for fundamental studies of lipid membrane processes, biotechnological applications such as drug delivery and biosensing, and more recently, for self-assembly. The latter has been made possible by inserting DNA oligomers with covalently linked hydrophobic anchors into the membrane. The lateral mobility of the DNA linkers on micrometer-sized droplets and solid particles has opened the door to creating structures with unprecedented structural flexibility. Here, we investigate micro-emulsions of TPM (3-(trimethoxysilyl)propyl methacrylate) as a platform for lipid monolayers and further functionalization with proteins and DNA oligonucleotides. TPM droplets can be produced with a narrow size distribution and are polymerizable, thus providing supports for model lipid membranes with controlled size and curvature. With fluorescence recovery after photobleaching, we observed that droplet-attached lipids, NeutrAvidin proteins, as well as DNA oligonucleotides all show mobility on the surface. We explored the assembly of micron-sized particles on TPM-droplets by exploiting either avidin-biotin interactions or double-stranded DNA with complementary single-stranded end groups. While the single molecules are mobile, the particles that are attached to them are not. We propose that this is caused by the heterogeneous nature of emulsified TPM, which forms an oligomer network that limits the collective motion of linkers, but allows the surface mobility of individual molecules.

Forming quasicrystals by monodisperse soft core particles

In traditional approaches to form quasicrystals, multiple competing length scales involved in particle size, shape, or interaction potential are believed to be necessary. It is unexpected that quasicrystals can be formed by monodisperse, isotropic particles interacting via a simple potential that does not contain explicit multiple length scales to stabilize quasicrystals. Here, we report the surprising finding of the formation of such quasicrystals in high-density systems of soft-core particles. Although there are length scales naturally introduced in our model systems, they do not establish the quasicrystalline order. In two dimensions, we find not only dodecagonal but also octagonal quasicrystals, which have not been found yet in soft quasicrystals. In such unexpected quasicrystals, particles tend to form pentagons, which are essential elements to develop the quasicrystalline order. Our findings thus pave an unexpected and simple way to form quasicrystals and pose a challenge for theoretical understanding of quasicrystals.

Quasicrystals with structural symmetries forbidden in crystals have been found in alloys or mono-component systems composed of anisotropic units. Zu et al. show a formation of two-dimensional quasicrystals in an isotropic soft solid with a spring-like potential, which challenges the existing theory.

Clathrate colloidal crystals

DNA-programmable assembly has been used to deliberately synthesize hundreds of different colloidal crystals spanning dozens of symmetries, but the complexity of the achieved structures has so far been limited to small unit cells. We assembled DNA-modified triangular bipyramids (~250-nanometer long edge, 177-nanometer short edge) into clathrate architectures. Electron microscopy images revealed that at least three different structures form as large single-domain architectures or as multidomain materials. Ordered assemblies, isostructural to clathrates, were identified with the help of molecular simulations and geometric analysis. These structures are the most sophisticated architectures made via programmable assembly, and their formation can be understood based on the shape of the nanoparticle building blocks and mode of DNA functionalization.

Novel Pyrochlorelike Crystal with a Photonic Band Gap Self-Assembled Using Colloids with a Simple Interaction Potential

Using computer simulations, we investigate the phase behavior of a system of particles interacting with a remarkably simple repulsive square-shoulder pair potential and report the formation of a novel (and stable) pyrochlorelike crystal phase. The lattice structure of the pyrochlorelike phase formed in our simulations possesses two inherent length scales corresponding to the inter- and intratetrahedral neighbors. We show that it can be used to fabricate a photonic crystal which displays complete photonic band gaps in both the direct and inverted dielectric structures.

Shaping micro-clusters via inverse jamming and topographic close-packing of microbombs

Designing topographic clusters is of significant interest, yet it remains challenging as they often lack mobility or deformability. Here we exploit the huge volumetric expansion (up to 3000%) of a new type of building block, thermally expandable microbombs. They consist of a viscoelastic polymeric shell and a volatile gas core, which, within structural confinement, create micro-clusters via inverse jamming and topographical close-packing. Upon heating, microbombs anchored in rigid confinement underwent balloon-like blowing up, allowing for dense clusters via soft interplay between viscoelastic shells. Importantly, the confinement is unyielding against the internal pressure of the microbombs, thereby enabling self-assembled clusters, which can be coupled with topographic inscription to introduce structural hierarchy on the clusters. Our strategy provides densely packed yet ultralight clusters with a variety of complex shapes, cleavages, curvatures, and hierarchy. In turn, these clusters will enrich our ability to explore the assemblies of the ever-increasing range of microparticle systems.

Self-assembled systems are normally composed of incompressible building blocks, which constrain their space filling efficiency. Yu et al. show programmable, densely packed clusters using thermally expandable soft microparticles, whereby the self-assembling process is realized via a jamming transition.

Effective potentials induced by self-assembly of patchy particles

Effective colloid-colloid interactions can be tailored through the addition of a complex cosolute. Here we investigate the case of a cosolute made by self-assembling patchy particles. Depending on the valence, these particles can form either polymer chains or branched structures. We numerically calculate the effective potential Veff between two colloids immersed in a suspension of reversible patchy particles, exploring a wide region of the cosolute phase diagram and the role of valence. In addition to well-known excluded volume and depletion effects, we find that, under appropriate conditions, Veff is completely attractive but shows an oscillatory character. In the case of polymerizing cosolute, this results from the fact that chains are efficiently confined by the colloids through the onset of local order. This argument is then generalized to the case of particles with higher valence, under the condition that they are still able to maintain a fully bonded organization upon confinement. The resulting effective potentials are relevant for understanding the behavior of complex mixtures in crowded environments, but may also be exploited for tuning colloidal self-assembly at preferred target distances in order to build desired superstructures.

Universal self-assembly of one-component three-dimensional dodecagonal quasicrystals

Using molecular dynamics simulations, we study computational self-assembly of one-component three-dimensional dodecagonal (12-fold) quasicrystals in systems with two-length-scale potentials. Existing criteria for three-dimensional quasicrystal formation are quite complicated and rather inconvenient for particle simulations. So to localize numerically the quasicrystal phase, one should usually simulate over a wide range of system parameters. We show how to universally localize the parameter values at which dodecagonal quasicrystal order may appear for a given particle system. For that purpose, we use a criterion recently proposed for predicting decagonal quasicrystal formation in one-component two-length-scale systems. The criterion is based on two dimensionless effective parameters describing the fluid structure which are extracted from the radial distribution function. The proposed method allows reduction of the time spent for searching the parameters favoring a certain solid structure for a given system. We show that the method works well for dodecagonal quasicrystals; this result is verified on four systems with different potentials: the Dzugutov potential, the oscillating potential which mimics metal interactions, the repulsive shoulder potential describing effective interactions for the core/shell model of colloids and the embedded-atom model potential for aluminum. Our results suggest that the mechanism of dodecagonal quasicrystal formation is universal for both metallic and soft-matter systems and it is based on competition between interparticle scales.

Quasicrystalline and crystalline types of local protein order in capsids of small viruses

Like metal alloys and micellar systems in soft matter, the viral capsid structures can be of crystalline and quasicrystalline types. We reveal the local quasicrystalline order of proteins in small spherical viral capsids using their nets of dodecahedral type. We show that the structure of some of the viral shells is well described in terms of a chiral pentagonal tiling, whose nodes coincide with centers of mass of protein molecules. The chiral protein packing found in these capsids originates from the pentagonal Penrose tiling (PPT), due to a specific phason reconstruction needed to fit the protein order at the adjacent dodecahedron faces. Via examples of small spherical viral shells and geminate capsid of a Maize Streak virus, we discuss the benefits and shortcomings of the usage of a dodecahedral net in comparison to icosahedral one, which is commonly applied for the modeling of viral shells with a crystalline local order.

On the formation of stripe, sigma, and honeycomb phases in a core-corona system

Using Monte Carlo simulations and free-energy calculations, we investigate the phase behaviour of a two-dimensional core-corona system. We model this system as particles consisting of an impenetrable hard core of diameter σ_HD surrounded by a purely repulsive soft corona of diameter δ = 1.95σ_HD. At low densities, we observe the spontaneous formation of a phase with a stripe texture as well as a honeycomb-like phase driven by both energy and entropy considerations. At high densities, we find that a two-dimensional analogue of the periodic sigma phase, considered as an approximant of dodecagonal quasicrystals, is energetically stabilised with respect to two distinct dodecagonal quasicrystals, namely, a square-triangle tiling and a square-triangle-shield tiling. We also find the formation of stable hexagonal phases at three distinct density ranges, which are energetically driven, i.e. by minimising the overlap of coronas. Furthermore, our calculations show that the low-density dodecagonal quasicrystal that was previously reported by Dotera et al., [Nature, 2014, 506, 208] is kinetically formed in the coexistence region between the honeycomb and the medium-density hexagonal phase.

Non-close-packed three-dimensional quasicrystals

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.

The effect of temperature, interaction range, and pair potential on the formation of dodecagonal quasicrystals in core-corona systems

A two-dimensional dodecagonal quasicrystal was previously reported by Dotera et al (2014 Nature 506 208) in a system of particles interacting with a hard core of diameter σ and a repulsive square shoulder of diameter [Formula: see text]. In the current work, we examine the formation of this quasicrystal using bond orientational order parameters, correlation functions and tiling distributions. We find that this dodecagonal quasicrystal forms from a fluid phase. We further study the effect of the width of the repulsive shoulder by simulating the system over a range of values of δ. For the range of densities and temperatures considered, we observe the formation of the dodecagonal quasicrystal between [Formula: see text] and [Formula: see text]. We also study the effect of shape of the interaction potential by simulating the system using three other interaction potentials with two length scales, namely hard-core plus a linear ramp, modified exponential, or Buckingham (exp-6) potential. We observe the presence of the quasicrystal in all three systems. However, depending on the shape of the potential, the formation of the quasicrystal takes place at lower temperatures (or higher interaction strengths). Using free-energy calculations, we demonstrate that the quasicrystal is thermodynamically stable in the square-shoulder and linear-ramp system.

Quasicrystalline nanocrystal superlattice with partial matching rules

Expanding the library of self-assembled superstructures provides insight into the behaviour of atomic crystals and supports the development of materials with mesoscale order. Here we build on 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 sublayers, forming a random tiling of rectangles, large triangles and small triangles with 6-fold symmetry. We analyse 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.

Self-assembly of two-dimensional binary quasicrystals: a possible route to a DNA quasicrystal

We use Monte Carlo simulations and free-energy techniques to show that binary solutions of penta- and hexavalent two-dimensional patchy particles can form thermodynamically stable quasicrystals even at very narrow patch widths, provided their patch interactions are chosen in an appropriate way. Such patchy particles can be thought of as a coarse-grained representation of DNA multi-arm 'star' motifs, which can be chosen to bond with one another very specifically by tuning the DNA sequences of the protruding arms. We explore several possible design strategies and conclude that DNA star tiles that are designed to interact with one another in a specific but not overly constrained way could potentially be used to construct soft quasicrystals in experiment. We verify that such star tiles can form stable dodecagonal motifs using oxDNA, a realistic coarse-grained model of DNA.

Screening Libraries of Semifluorinated Arylene Bisimides to Discover and Predict Thermodynamically Controlled Helical Crystallization

Synthesis, structural, and retrostructural analysis of a library containing 16 self-assembling perylene (PBI), 1,6,7,12-tetrachloroperylene (Cl₄PBI), naphthalene (NBI), and pyromellitic (PMBI) bisimides functionalized with environmentally friendly AB₃ chiral racemic semifluorinated minidendrons at their imide groups via m = 0, 1, 2, and 3 methylene units is reported. These semifluorinated compounds melt at lower temperatures than homologous hydrogenated compounds, permitting screening of all their thermotropic phases via structural analysis to discover thermodynamically controlled helical crystallization from propeller-like, cogwheel, and tilted molecules as well as lamellar-like structures. Thermodynamically controlled helical crystallization was discovered for propeller-like PBI, Cl₄PBI and NBI with m = 0. Unexpectedly, assemblies of twisted Cl₄PBIs exhibit higher order than those of planar PBIs. PBI with m = 1, 2, and 3 form a thermodynamically controlled columnar hexagonal 2D lattice of tilted helical columns with intracolumnar order. PBI and Cl₄PBI with m = 1 crystallize via a recently discovered helical cogwheel mechanism, while NBI and PMBI with m = 1 form tilted helical columns. PBI, NBI and PMBI with m = 2 generate lamellar-like structures. 3D and 2D assemblies of PBI with m = 1, 2, and 3, NBI with m = 1 and PMBI with m = 2 exhibit 3.4 Å π-π stacking. The library approach applied here and in previous work enabled the discovery of six assemblies which self-organize via thermodynamic control into 3D and 2D periodic arrays, and provides molecular principles to predict the supramolecular structure of electronically active components.

Pushing the limits of crystallography

A very serious concern of scientists dealing with crystal structure refinement, including theoretical research, pertains to the characteristic bias in calculated versus measured diffraction intensities, observed particularly in the weak reflection regime. This bias is here attributed to corrective factors for phonons and, even more distinctly, phasons, and credible proof supporting this assumption is given. The lack of a consistent theory of phasons in quasicrystals significantly contributes to this characteristic bias. It is shown that the most commonly used exponential Debye-Waller factor for phasons fails in the case of quasicrystals, and a novel method of calculating the correction factor within a statistical approach is proposed. The results obtained for model quasiperiodic systems show that phasonic perturbations can be successfully described and refinement fits of high quality are achievable. The standard Debye-Waller factor for phonons works equally well for periodic and quasiperiodic crystals, and it is only in the last steps of a refinement that different correction functions need to be applied to improve the fit quality.

Complex Arrangement of Orthogonal Nanoscale Columns via a Supramolecular Orientational Memory Effect

Memory effects, including shape, chirality, and liquid-crystallinity, have enabled macroscopic materials with novel functions. However, the generation of complex supramolecular nanosystems via memory effects has not yet been investigated. Here, we report a cyclotriveratrylene-crown (CTV) compound that self-assembles into supramolecular columns and spheres forming, respectively, hexagonal and cubic mesophases. Upon transition from one phase to the other, an epitaxial relationship holds, via an unprecedented supramolecular orientational memory effect. Specifically, the molecular orientation and columnar character of supramolecular packing is preserved in the cubic phase, providing an otherwise inaccessible structure comprising orthogonally oriented domains of supramolecular columns. The continuous columnar character of tetrahedrally distorted supramolecular spheres self-organized from the CTV derivative in the faces of the Pm3̅n lattice is the basis of this supramolecular orientational memory, which holds throughout cycling in temperature between the two phases. This concept is expected to be general for other combinations of periodic and quasiperiodic arrays generated from supramolecular spheres upon transition to supramolecular columns.

A smectic dodecagonal quasicrystal

We report a solid smectic phase that exhibits dodecagonal global order. It is composed of axially stacked hexagonally ordered particle layers, and its 12-fold rotational symmetry induced by the 30° rotation of adjacent layers with respect to each other. A quasicrystal was produced in a molecular-dynamics simulation of a single-component system of particles interacting via a spherically-symmetric potential. It was formed as a result of a first-order phase transition from an isotropic liquid state that occurred under constant-density cooling. This finding implies that a similarly structured quasicrystal can possibly be produced by the same class of systems as those forming smectic-B crystals. This quasicrystal can also be expected to arise in a system of spherically-shaped colloidal particles with appropriately tuned potential.

Molecular dynamics study of colloidal quasicrystals

Colloidal quasicrystals have received increased interest recently due to new insight in exploring their potential for photonic materials as well as for optical devices [Vardeny et al., Nat. Photonics, 2013, 7, 177]. Colloidal quasicrystals in aqueous solutions have been found in systems of micelles with impenetrable cores [Fischer et al., Proc. Natl. Acad. Sci. U. S. A., 2011, 108, 1810]. A simple model potential for micelle-micelle interaction is the step potential, which is infinite for core overlaps and constant for shell overlaps. Dotera et al. performed Monte Carlo simulations of the step potential model and found quasicrystals for specific values of the packing fraction η and the shell-core ratio λ [Dotera et al., Nature, 2014, 506, 208 ]. However, the overlap of real micelles causes repulsive forces, which increase with decreasing core distance. We consider this by introducing a novel model potential with repulsive forces depending on a third parameter α. In a systematic manner we study this more realistic potential with two-dimensional molecular dynamics simulations. For α = 0 the model is similar to the step potential model. For the first time, we provide a comprehensive overview of crystalline, quasicrystalline, and disordered structures as a function of η and λ. Simulations performed with α > 0 show the impact of the repulsive forces. We find that quasicrystalline structures at high densities vanish while new quasicrystalline structures appear at intermediate densities. Our results help to tailor colloidal systems for today's advanced applications in photonics and optical devices.

Three-Dimensional Icosahedral Phase Field Quasicrystal

We investigate the formation and stability of icosahedral quasicrystalline structures using a dynamic phase field crystal model. Nonlinear interactions between density waves at two length scales stabilize three-dimensional quasicrystals. We determine the phase diagram and parameter values required for the quasicrystal to be the global minimum free energy state. We demonstrate that traits that promote the formation of two-dimensional quasicrystals are extant in three dimensions, and highlight the characteristics required for three-dimensional soft matter quasicrystal formation.

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

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

Cornucopia of Nanoscale Ordered Phases in Sphere-Forming Tetrablock Terpolymers

We report the phase behavior of a series of poly(styrene)-b-poly(isoprene)-b-poly(styrene)'-b-poly(ethylene oxide) (SIS'O) tetrablock terpolymers. This study was motivated by self-consistent field theory (SCFT) calculations that anticipate a rich array of sphere-forming morphologies with variations in the molecular symmetry parameter τ = NS/(NS + NS'), where N is the block degree of polymerization and the volume fraction of O is less than about 0.22. Eight SIS'O samples, with τ ranging from 0.21 to 0.73, were synthesized and investigated using small-angle X-ray scattering and transmission electron microscopy, yielding evidence of nine different spherical phases: hexagonal, FCC, HCP, BCC, rhombohedral (tentative), liquid-like packing, dodecagonal quasicrystal, and Frank-Kasper σ and A15 phases. At temperatures close to the order-disorder transition, these tetrablocks behave as pseudo-[SIS']-O diblocks and form equilibrium morphologies mediated by facile chain exchange between micelles. Transition from equilibrium to nonequilibrium behavior occurs at a temperature (Terg) several tens of degrees below the order-disorder transition temperature, speculated to be coincident with the loss of ergodicity, as chain exchange is arrested due to increased segregation strength between the core (O) and corona (SIS') blocks. Nonequilibrium ordered structures form when T < Terg; these are interpreted using SCFT calculations to elucidate the free energy landscape driving ordering in the S and I block matrix. These experiments demonstrate a profound dependence on phase stability with variations in τ and temperature, providing insights into the formation of ordered phase symmetry in this class of asymmetric multiblock polymers.

Emergence of Chirality from Isotropic Interactions of Three Length Scales

Chirality is known to play a pivotal role in determining material properties and functionalities. However, it remains a great challenge to understand and control the emergence of chirality and the related enantioselective process particularly when the building components of the system are achiral. Here we explore the generic mechanisms driving the formation of two-dimensional chiral structures in systems characterized by isotropic interactions and three competing length scales. We demonstrate that starting from isotropic and rotationally invariant interactions, a variety of chiral ordered patterns and superlattices with anisotropic but achiral units can self-assemble. The mechanisms for selecting specific states are related to the length-scale coupling and the selection of resonant density wave vectors. Sample phase diagrams and chiral elastic properties are identified. These findings provide a viable route for predicting chiral phases and selecting the desired handedness.

Equilibrium cluster fluids: pair interactions via inverse design

Inverse methods of statistical mechanics are becoming productive tools in the design of materials with specific microstructures or properties. While initial studies have focused on solid-state design targets (e.g., assembly of colloidal superlattices), one can alternatively design fluid states with desired morphologies. This work addresses the latter and demonstrates how a simple iterative Boltzmann inversion strategy can be used to determine the isotropic pair potential that reproduces the radial distribution function of a fluid of amorphous clusters with prescribed size. The inverse designed pair potential of this "ideal" cluster fluid, with its broad attractive well and narrow repulsive barrier at larger separations, is qualitatively different from the so-called SALR form most commonly associated with equilibrium cluster formation in colloids, which features short-range attractive (SA) and long-range repulsive (LR) contributions. These differences reflect alternative mechanisms for promoting cluster formation with an isotropic pair potential, and they in turn produce structured fluids with qualitatively different static and dynamic properties. Specifically, equilibrium simulations show that the amorphous clusters resulting from the inverse designed potentials display more uniformity in size and shape, and they also show greater spatial and temporal resolution than those resulting from SALR interactions.

Mathematical crystallography in the 21st century

As crystallography merges with materials science and engineering, mathematical crystallography is growing in new directions, including: Characterizing new materials with unusual properties; Imaging, including but not limited to diffraction; Exploring and exploiting superspaces; Mapping the aperiodic landscape, from chaos to classical periodicity and beyond; Re-modeling the structures of real crystals, both periodic and aperiodic; Modeling self-assembly and self-reorganization on the nanoscale. In short, it’s not (just) about space groups and tilings anymore.

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 a direct connection between building-block attributes and bulk structure 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 structure 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. Although our focus is on colloidal systems, our methods generalize to any building blocks with adjustable interactions.

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, thus 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.

Self-assembly of the decagonal quasicrystalline order in simple three-dimensional systems

Using molecular dynamics simulations we show that a one-component system can be driven to a three-dimensional decagonal (10-fold) quasicrystalline state just by purely repulsive, isotropic and monotonic interaction pair potential with two characteristic length scales; no attraction is needed. We found that self-assembly of a decagonal quasicrystal from a fluid can be predicted by two dimensionless effective parameters describing the fluid structure. We demonstrate stability of the results under changes of the potential by obtaining the decagonal order for three particle systems with different interaction potentials, both purely repulsive and attractive, but with the same values of the effective parameters. Our results suggest that soft matter quasicrystals with decagonal symmetry can be experimentally observed for the same systems demonstrating the dodecagonal order for an appropriate tuning of the effective parameters.

Experimental Observation of Quasicrystal Growth

The growth of an Al-Ni-Co decagonal quasicrystal was observed by in situ, high-temperature, high-resolution transmission electron microscopy. The tiling patterns extracted from a series of high-resolution transmission electron microscopy images were analyzed on the basis of the high-dimensional description of quasicrystalline structures. The analyses indicated that the growth proceeded with frequent error-and-repair processes. The final, grown structure showed nearly perfect quasicrystalline order. Our observations suggest that the repair process by phason relaxation, rather than local growth rule, plays an essential role in the construction of ideal quasicrystalline order in real materials.

Soft-core particles freezing to form a quasicrystal and a crystal-liquid phase

Systems of soft-core particles interacting via a two-scale potential are studied. The potential is responsible for peaks in the structure factor of the liquid state at two different but comparable length scales and a similar bimodal structure is evident in the dispersion relation. Dynamical density functional theory in two dimensions is used to identify two unusual states of this system: a crystal-liquid state, in which the majority of the particles are located on lattice sites but a minority remains free and so behaves like a liquid, and a 12-fold quasicrystalline state. Both are present even for deeply quenched liquids and are found in a regime in which the liquid is unstable with respect to modulations on the smaller scale only. As a result, the system initially evolves towards a small-scale crystal state; this state is not a minimum of the free energy, however, and so the system subsequently attempts to reorganize to generate the lower-energy larger-scale crystals. This dynamical process generates a disordered state with quasicrystalline domains and takes place even when this large scale is linearly stable, i.e., it is a nonlinear process. With controlled initial conditions, a perfect quasicrystal can form. The results are corroborated using Brownian dynamics simulations.

Effective substrate potentials with quasicrystalline symmetry depend on the size of the adsorbed particles

We explore the effective potential landscapes that extended particles experience when adsorbed on the surface of quasicrystals. Commonly, these are solids with long-ranged order but no translational symmetry. The effective potentials significantly depend on the size of the adsorbed particles. We show how changing the particle radius changes the so-called local isomorphism class of the effective quasicrystalline pattern. This means effective potentials for different particle sizes cannot directly be mapped onto each other. Our theoretical predictions are confirmed by Monte Carlo simulations. The results are important for colloidal particles with different sizes that are subjected to laser fields with quasicrystalline symmetry as well as for systems where extended molecules are deposited onto the surface of metallic quasicrystals.

Biomolecular engineering of virus-like particles aided by computational chemistry methods

Multi-scale investigation of VLP self-assembly aided by computational methods is facilitating the design, redesign, and modification of functionalized VLPs.

Atomic structure of icosahedral quasicrystals: stacking multiple quasi-unit cells

An effective tiling approach is proposed for the structural description of icosahedral quasicrystals based on the original substitution algorithm.