Icosahedral quasicrystals (IQCs) are materials that exhibit long-range order but lack periodicity in any direction. Although IQCs were the first reported quasicrystals1, they have been experimentally observed only in metallic alloys2, not in other materials. By contrast, quasicrystals with other symmetries (particularly dodecagonal) have now been found in several soft-matter systems3,4,5. Here we introduce a class of IQCs built from model patchy colloids that could be realized experimentally using DNA origami particles. Our rational design strategy leads to systems that robustly assemble in simulations into a target IQC through directional bonding. This is illustrated for both body-centred and primitive IQCs, with the simplest systems involving just two particle types. The key design feature is the geometry of the interparticle interactions favouring the propagation of an icosahedral network of bonds, despite this leading to many particles not being fully bonded. As well as furnishing model systems in which to explore the fundamental physics of IQCs, our approach provides a potential route towards functional quasicrystalline materials.
This is a preview of subscription content
Subscription info for Chinese customers
We have a dedicated website for our Chinese customers. Please go to naturechina.com to subscribe to this journal.
Rent or Buy article
Get time limited or full article access on ReadCube.
All prices are NET prices.
Shechtman, D., Blech, I., Gratias, D. & Cahn, J. W. Metallic phase with long-range orientational order and no translational symmetry. Phys. Rev. Lett. 53, 1951–1953 (1984).
Takakura, H., Gómez, C. P., Yamamoto, A., Boissieu, M. D. & Tsai, A. P. Atomic structure of the binary icosahedral Yb–Cd quasicrystal. Nat. Mater. 6, 58–63 (2007).
Zeng, X. et al. Supramolecular dendritic quasicrystals. Nature 428, 157–160 (2004).
Haji-Akbari, A. et al. Disordered, quasicrystalline and crystalline phases of densely packed tetrahedra. Nature 462, 773–777 (2009).
Dotera, T. Quasicrystals in soft matter. Isr. J. Chem. 51, 1197–1205 (2011).
Dotera, T., Oshiro, T. & Ziherl, P. Mosaic two-length scale quasicrystals. Nature 506, 208–211 (2014).
Savitz, S., Babadi, M. & Lifshitz, R. Multiple-scale structures: from Faraday waves to soft-matter quasicrystals. IUCrJ 5, 247–268 (2018).
Subramanian, P., Archer, A. J., Knobloch, E. & Rucklidge, A. M. Three-dimensional icosahedral phase field quasicrystal. Phys. Rev. Lett. 117, 075501 (2016).
Engel, M., Damasceno, P. F., Phillips, C. L. & Glotzer, S. C. Computational self-assembly of a one-component icosahedral quasicrystal. Nat. Mater. 14, 109–116 (2015).
Damasceno, P. F., Glotzer, S. C. & Engel, M. Non-close-packed three-dimensional quasicrystals. J. Phys. Condens. Matter 29, 234005 (2017).
Reinhardt, A., Schreck, J. S., Romano, F. & Doye, J. P. K. Self-assembly of two-dimensional binary quasicrystals: a possible route to a DNA quasicrystal. J. Phys. Condens. Matter 29, 014006 (2017).
Sadoc, J.-F. & Mosseri, R. Geometric Frustration (Cambridge Univ. Press, 1999).
Doye, J. P. K. & Wales, D. J. Structural consequences of the range of the interatomic potential: a menagerie of clusters. J. Chem. Soc. Faraday Trans. 93, 4233–4243 (1997).
Dmitrienko, V. & Kléman, M. Tetrahedral structures with icosahedral order and their relation to quasi-crystals. Crystallogr. Rep. 46, 527–533 (2001).
Liu, L., Li, Z., Li, Y. & Mao, C. Rational design and self-assembly of two-dimensional dodecagonal quasicrystals. J. Am. Chem. Soc. 141, 4248–4251 (2019).
Tracey, D. F., Noya, E. G. & Doye, J. P. K. Programming patchy particles to form complex ordered structures. J. Chem. Phys. 151, 224506 (2019).
Janot, C. Quasicrystals: A Primer (Oxford Univ. Press, 2002).
Cademartiri, L. & Bishop, K. J. M. Programmable self-assembly. Nat. Mater. 14, 2–9 (2015).
Liu, W., Halverson, J., Tian, Y., Tkachenko, A. V. & Gang, O. Self-organized architectures from assorted DNA-framed nanoparticles. Nat. Chem. 8, 867–873 (2016).
Zhang, Z., Keys, A. S., Chen, T. & Glotzer, S. C. Self-assembly of patchy particles into diamond structures through molecular mimicry. Langmuir 21, 11547–11551 (2005).
Vega, C., Sanz, E., Abascal, J. L. F. & Noya, E. G. Determination of phase diagrams via computer simulation: methodology and applications to water, electrolytes and proteins. J. Phys. Condens. Matter 20, 153101 (2008).
Glotzer, S. C. & Solomon, M. J. Anisotropy of building blocks and their assembly into complex structures. Nat. Mater. 6, 557–562 (2007).
Li, W. et al. Colloidal molecules and patchy particles: complementary concepts, synthesis and self-assembly. Chem. Soc. Rev. 49, 1955–1976 (2020).
Xiong, Y. et al. Three-dimensional patterning of nanoparticles by molecular stamping. ACS Nano 14, 6823–6833 (2020).
Veneziano, R. et al. Designer nanoscale DNA assemblies programmed from the top down. Science 352, 1534 (2016).
Tian, Y. et al. Ordered three-dimensional nanomaterials using DNA-prescribed and valence-controlled material voxels. Nat. Mater. 19, 789–796 (2020).
Ma, N. et al. Directional assembly of nanoparticles by DNA shapes: towards designed architectures and functionality. Top. Curr. Chem. 378, 36 (2020).
Man, W., Megens, M., Steinhardt, P. J. & Chaikin, P. M. Experimental measurement of the photonic properties of icosahedral quasicrystals. Nature 436, 993–996 (2005).
Steurer, W. Boron-based quasicrystals with sevenfold symmetry. Phil. Mag. 87, 2707–2712 (2007).
Snodin, B. E. K. et al. Introducing improved structural properties and salt dependence into a coarse-grained model of DNA. J. Chem. Phys. 142, 234901 (2015).
Wilber, A. W., Doye, J. P. K., Louis, A. A. & Lewis, A. C. F. Monodisperse self-assembly in a model with protein-like interactions. J. Chem. Phys. 131, 175102 (2009).
Ben Zion, M. Y. et al. Self-assembled three-dimensional chiral colloidal architecture. Science 358, 633–636 (2017).
Anderson, J. A., Jankowski, E., Grubb, T. L., Engel, M. & Glotzer, S. C. Massively parallel Monte Carlo for many-particle simulations on GPUs. J. Comput. Phys. 254, 27–38 (2013).
Rapaport, D. C. The Art of Molecular Dynamics Simulation (Cambridge Univ. Press, 2004).
Engel, M., Umezaki, M., Trebin, H.-R. & Odagaki, T. Dynamics of particle flips in two-dimensional quasicrystals. Phys. Rev. B 82, 134206 (2010).
Rovigatti, L., Šulc, P., Reguly, I. Z. & Romano, F. A comparison between parallelization approaches in molecular dynamics simulations on GPUs. J. Comput. Chem. 36, 1–8 (2015).
Bai, X.-C., Martin, T. G., Scheres, S. H. W. & Dietz, H. Cryo-EM structure of a 3D DNA-origami object. Proc. Natl Acad. Sci. USA 109, 20012–20017 (2012).
Snodin, B. E. K., Schreck, J. S., Romano, F., Louis, A. A. & Doye, J. P. K. Coarse-grained modelling of the structural properties of DNA origami. Nucleic Acids Res. 47, 1585–1597 (2019).
Jun, H. et al. Automated sequence design of 3D polyhedral wireframe DNA origami with honeycomb edges. ACS Nano 13, 2083–2093 (2019).
Douglas, S. M. et al. Rapid prototyping of 3D DNA-origami shapes with caDNAno. Nucleic Acids Res. 37, 5001–5006 (2009).
Suma, A. et al. TacoxDNA: a user-friendly web server for simulations of complex DNA structures, from single strands to origami. J. Comput. Chem. 40, 2586–2595 (2019).
Dietz, H., Douglas, S. M. & Shih, W. M. Folding DNA into twisted and curved nanoscale shapes. Science 325, 725–730 (2009).
We are grateful for financial support from the Agencia Estatal de Investigación and the Fondo Europeo de Desarrollo Regional (FEDER) under grant numbers FIS2015-72946-EXP(AEI) and FIS2017-89361-C3-2-P(AEI/FEDER,UE) (to E.G.N.) and the Croucher Foundation (to C.K.W.). We acknowledge the use of the University of Oxford Advanced Research Computing (ARC) facility and the Cambridge Service for Data Driven Discovery (CSD3).
The authors declare no competing interests.
Peer review information Nature thanks Oleg Gang, Ron Lifshitz and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Projections of the BCI 5P (left) and 2P (right) IQCs along the five-fold axis, three-fold axis and two-fold axis (top to bottom). On the three-fold axis, self-similar triangular motifs follow a 1:τ:τ2:τ3 scaling that is a characteristic of IQCs.
Comparison of the diffraction patterns of the ideal body-centred IQC and its 3/2 approximant, with those of the assembled IQCs using the BCI 5P, 3P, 2P and 2P-mod models. All four assembled systems have diffraction patterns with clear five-fold symmetry confirming their IQC character. By contrast, the diffraction pattern of the approximant when viewed along the pseudo-five-fold axis shows clear deviations from fivefold symmetry. The diffraction patterns of the BCI 5P and 3P models most closely resemble those of the ideal IQC; the indexing of all peaks with a BCI indexing scheme is detailed in Supplementary Information section 4E. The diffraction patterns for the BCI 2P and 2P-mod models exhibit additional features, particularly along the two-fold and five-fold axes; these include both diffuse scattering and extra weak peaks that cannot be indexed in a BCI indexing scheme (see Supplementary Information).
a, Primitive IQC model with a dodecahedral occupation domain viewed along the five-, three- and two-fold rotational axes (left to right). In the top row, only red and those blue particles that are bonded to red particles are shown; these views clearly show that not all the clusters are perfect. In the bottom row, all particles are shown. The yellow particles form an additional shell around the icosahedral cluster depicted in Fig. 3a and all other matrix particles are coloured pink. b, Views of the 8/5 rational approximant along the psuedo-five-, three- and two-fold rotational axes (left to right).
Projections of the primitive icosahedral 3P WT (left) and 2P WT (right) models along the five-fold axis, three-fold axis and two-fold axis (top to bottom). On the three-fold axis, self-similar triangular motifs follow a 1:τ:τ2:τ3 scaling that is a characteristic of IQCs.
Comparison of the diffraction patterns of the ideal primitive IQC and its 8/5 approximant with those of the IQCs obtained from simulations. All four assembled systems have diffraction patterns with clear five-fold symmetry confirming their IQC character. By contrast, the diffraction pattern of the approximant when viewed along the pseudo-five-fold axis shows clear deviations from five-fold symmetry. Despite the assembled systems being simpler than the ideal IQC in terms of the number of environments, the diffraction patterns are very similar to that for the ideal IQC with the patterns along the two-fold axis showing the most differences. The primitive icosahedral 2P diffraction patterns could be satisfactorily indexed with a primitive icosahedral indexing scheme, but the better-resolved primitive icosahedral 3P patterns have additional weak peaks that require a face-centred icosahedral indexing scheme with a hypercubic lattice parameter that is double that of the primitive icosahedral indexing scheme, indicating some superstructural ordering (see Supplementary Information section 4E). The diffraction patterns for the systems with and without torsions show no clear differences, confirming their structural similarity.
The properties of the body-centred IQCs obtained with the 5P, 3P and 2P models are compared with those of the ideal body-centred IQC and its 3/2 rational approximant. a, PDFs considering only those particles within a sphere of radius 15σLJ centred at the centre of mass of the solid cluster. There is good coincidence between the broadened (due to thermal motion) peaks for the assembled IQCs and the distances present in the ideal IQC consistent with their similar structure. The PDFs for the BCI 5P and 3P systems are essentially the same, with the differences from the BCI 2P system still being relatively small. b, Radial density. The bulk densities of the ideal IQC, the BCI 5P and 3P systems and the 3/2 approximant are very similar, but the bulk density for the BCI 2P system is lower, as expected, due to the absence of the red particles. The radial density of the assembled systems decreases from its bulk value to zero over 10σLJ–15σLJ, reflecting the irregular nature of the surface of the assemblies. c, Mole fraction of each particle type in the different structures. The incorporation of the yellow and purple particles into the BCI 5P assembly is much lower than for the ideal IQC. The fraction of red and blue particles in the BCI 5P and 3P assemblies is also closer to that for the approximant than the ideal IQC. The BCI 2P assemblies have a substantially increased fraction of green particles. d–i, Probability of the number of dangling bonds for each particle type (d–h) and for all particle types (i). The assembled IQCs show a greater number of dangling bonds for the blue matrix particles than the ideal IQC, whereas the ideal IQC has a greater number of dangling bonds for the red and green particles that form the triacontahedral clusters.
The properties of the primitive IQCs obtained with the 3P and 2P models with and without torsions are compared with those of the ideal primitive IQC. a, PDF considering only those particles within a sphere of radius 15σLJ centred at the centre of mass of the solid cluster. The peaks in the PDFs are consistent with the distances in the ideal IQC, even though the assembled systems do not have particles corresponding to many of the environments in the ideal IQC. The 2P system does not have the peak at ~0.8σLJ that is due to the red–red bonds. b, Radial density. The bulk densities of the 3P systems are a bit lower than the ideal IQC with the 2P systems having even lower densities due to the absence of red particles. c, Mole fraction of each particle type in the different structures. Note that in the ideal structure the mole fractions do not add to 1, as there are additional local environments that are not included in this plot. Particularly noticeable is the much smaller number of red particles in the ideal IQC. d–g, Probability of the number of dangling bonds for each particle type (d–f) and for all particle types (g). The blue particles in the 3P systems have a particularly low probability of having no missing bonds, because when forming part of the inter-cluster matrix, they do not use their two patches that can only bond to the red particles that form the dodecahedral clusters; such particles contribute to the relatively large number of blue particles with two dangling bonds.
a–d, Contributions to the BOOD for the BCI 5P (a), BCI 2P (b), primitive icosahedral 3P WT (c) and primitive icosahedral 2P WT (d) quasicrystals from the different patch–patch bonds types in these systems. For the BCI systems, the green–green and purple–purple bonds are directed along the two-fold axes, the blue–green and blue–purple along the three-fold axes, the blue–red, red–green, blue–yellow and yellow–purple along the five-fold axes and the blue–blue can be directed along five-fold (bonds formed through patches 1 to 3) and three-fold (bonds formed by patch 4) axes. For the primitive icosahedral systems, the red–red bonds are directed along the two-fold axes, and the green–blue and red–blue along the three-fold axes. Thus, for the primitive icosahedral 2P WT system the total BOOD only exhibits spots directed along the three-fold axes, whereas for the primitive icosahedral 3P WT system there are also spots along the two-fold axes. Contributions to the BOODs calculated using a distance cutoff exhibit additional spots but always directed along the rotational axes of the Ih point group (Supplementary Fig. 7). e, f, The relationship between the symmetry of the patchy particles and the Ih point group. All particles can be oriented so that their patch vectors point exclusively along the rotational axes of Ih. This is illustrated for the patchy particles in the BCI 5P (e) and primitive icosahedral 3P (f) models by depicting the direction of the patch vectors Pi by circles (with the colours representing (one of) the particles with which they interact) on the surface of an icosahedron. The projections chosen are along the highest rotational axis of the particles. Edges on the back faces of the icosahedra are dashed and cyan. Similarly, for patch vectors on the back faces, the colour shade is lighter and ringed in cyan rather than black. The BOODs in a–d are fully consistent with the particles being oriented in this way with respect to the global icosahedral order of the IQCs. We also note that the preferred torsional angles of each patch–patch interaction are those that ensure the propagation of this global orientational order.
Van Hove autocorrelation functions for the simulated body-centred IQCs, evaluated after one million Monte Carlo cycles and considering only those particles within a radius of 20σLJ, which corresponds to the interior of the cluster where the radial density is constant (Extended Data Fig. 6b). The majority of the particle mobility is associated with the blue matrix particles, whereas the red and green particles that form the inner shell of the icosahedral clusters have very limited mobility. The yellow particles in the 5P system are the most mobile but are only present at very low mole fraction.
Van Hove autocorrelation function for the simulated primitive IQCs, evaluated after one million Monte Carlo cycles and considering only those particles within a radius of 20σLJ, which corresponds to the interior of the cluster where the radial density is constant (Extended Data Fig. 7b). Although the overall pattern shows clear five-fold symmetry, the peaks merge into each other more than for the BCI systems. This is partly because the size of the allowed hops is shorter. The green particles are the most mobile, and the red particles, which form the inner shell of the icosahedral clusters, are least mobile. It is also noticeable that without torsions, the 3P system has greater mobility, and that the 2P patterns are less well defined, with the motion of the blue particles close to isotropic.
About this article
Cite this article
Noya, E.G., Wong, C.K., Llombart, P. et al. How to design an icosahedral quasicrystal through directional bonding. Nature 596, 367–371 (2021). https://doi.org/10.1038/s41586-021-03700-2