Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Nature volume  615, pages 418–424 (2023 )Cite this article

Chirality is a geometrical property described by continuous mathematical functions1,2,3,4,5. However, in chemical disciplines, chirality is often treated as a binary left or right characteristic of molecules rather than a continuity of chiral shapes. Although they are theoretically possible, a family of stable chemical structures with similar shapes and progressively tuneable chirality is yet unknown. Here we show that nanostructured microparticles with an anisotropic bowtie shape display chirality continuum and can be made with widely tuneable twist angle, pitch, width, thickness and length. The self-limited assembly of the bowties enables high synthetic reproducibility, size monodispersity and computational predictability of their geometries for different assembly conditions6. The bowtie nanoassemblies show several strong circular dichroism peaks originating from absorptive and scattering phenomena. Unlike classical chiral molecules, these particles show a continuum of chirality measures2 that correlate exponentially with the spectral positions of the circular dichroism peaks. Bowtie particles with variable polarization rotation were used to print photonically active metasurfaces with spectrally tuneable positive or negative polarization signatures for light detection and ranging (LIDAR) devices.

Mathematical definitions of mirror asymmetry1,2,3,4,5 recognize the continuity of chiral geometries that can be visualized by stretching macroscale helical springs to obtain coils of different lengths, that is, pitch. At smaller scales, continuously variable chirality can be observed for origami/kirigami sheets7,8,9, nanocomposites10,11 and polymeric solids whose shape and circular dichroism (CD) can be varied by external fields12. However, chirality in chemistry commonly manifests as a binary property: chiral molecules are either right- or left-handed and the descriptors for stereochemical configurations of enantiomers at molecular scale are correspondingly binary, for example d /l , R/S, M/P and Δ/Λ. The binary chirality of amino acids originates from a high-energy penalty for distorting the optical centre based on sp3-carbon atoms. The discreteness of chirality in liquid crystals13, macromolecular compounds14, helical polymers15 and nanoparticles (NPs)16,17 manifests in abrupt transitions between chiral phases with different crystallinity or particle shapes18. Energy penalties are relaxed for large flexible molecules19, supramolecular complexes16,20,21 and biopolymers18,22,23, but the restrictions on chiral shapes remain stringent. Unlike macroscale springs, the helical pitch across the variety of biomolecules varies little. For DNA, protein α-helices and β-sheets, pitch changes only in the ranges of 11–46, 2.3–5.5 and 7–8 Å, respectively; these narrow ranges are essential for the precise folding of biomolecules24.

The transition from discrete chiral phases and shapes to a palette of chemical compounds with continuously tuneable chirality would be transformative for the development of chiral photonics, chiral metamaterials, biochemical separations and chiral catalysis. The availability of continuously variable chiral compounds is essential for establishing fundamental correlations between chirality measures and chemical properties. For example, the attempts to correlate optical activity with a variety of chirality measures have largely failed25,26, but they could be possible for chiral nanostructures and their assemblies due to the differences in physics of chiroptical activity in them compared to molecules with binary chirality

Continuously variable chiral geometries become possible for nanostructured microparticles with bowtie shapes. They are hierarchically assembled from nanoribbons containing helical chains of cystine (CST, dipeptide of cysteine through the S–S bridge) interconnected by Cd2+ ions. The balance between short- and long-range interactions and defect tolerance of the electrostatically restricted assembly process27 enables the synthesis of bowties with widely tuneable pitch, width, thickness and length.

Mixing an aqueous solutions of Cd2+ and l - or d -cystine (Supplementary Fig. 1a,b) results in bowties (Methods) with pitch (p), width (w), thickness (t) and length (l) of 4.1, 1.3, 0.5 and 3.1 µm, respectively (Fig. 1a). The process is 100% enantioselective with l - and d -CST resulting in correspondingly left- and right-twisted bowties (Fig. 1a,c, and Supplementary Video 2); they are uniform in shape, size and handedness with standard deviation in p, w, t and l of 19.0, 13.6, 11.9 and 11.5%, respectively. Monodispersity of the bowtie particles indicates that they are formed in a self-limited assembly process (Supplementary Fig. 14)5,6. Nanostructured microparticles resembling a stack of flat nanoscale ‘pancakes’ are observed for rac-CST (Fig. 1b).

a–c, SEM images of bowties formed from Cd2+ and l -CST (left) (a), rac-CST (middle) (b) and d -CST (right) (c). Scale bars, 2 µm. d–f, CD (top) and extinction spectra (bottom) in THz (d), mid-IR (VCD) (e) and UV-visible-NIR (f) ranges for bowtie particles. g, Zeta potential of fully formed bowties at different ratios of l -CST and Cd2+ ions. h, Zeta potential of bowties formed from l -CST solution at different pHs. i–k, Synchrotron XRD pattern (i) and the calculated crystal structure of nanosheets (j) and nanoplatelets (k) forming the bowties. AU, arbitrary units. The unit cell is shown in Supplementary Fig. 6.

Electron microscopy, X-ray diffraction (XRD) and electron diffraction show levels of hierarchical organization within the cadmium cystinate bowties (Figs. 1i and 2a). Scanning electron microscopy (SEM) images for the terminal (Fig. 1a–c) and intermediate (Fig. 2a,b) stages of the synthesis show that bowties are structured as a stack of twisted nanoribbons 200–1,200 nm in length and 45 nm in thickness (Supplementary Video 1). The nanoribbons are assembled from nanoplatelets 50–200 nm in length (Fig. 2c,d) and a thickness of roughly 1.2 nm calculated from the solved crystal structure (Supplementary Table 1, Supplementary Fig. 6 and Fig. 1k). Structural data in Supplementary Fig. 12 indicate that the nanoplatelets are formed from nanoclusters 1–2 nm in thickness, which is consistent with the thickness of the nanoplatelets and previous studies28. The nanoribbons and bowties acquire progressively stronger twists as the enantiomeric excess (±χ) of l - or d -CST increases (Fig. 2f).

a,b, SEM image of a fully formed bowtie (a), composed of twisted segments that are assembled from nanoribbons (b). c, SEM image of nanoribbons composed of polydisperse nanoplatelets. d, TEM image of the platelet. e. cryo-SAED pattern of the nanosheet in d. f, SEM images of the bowtie particles obtained by mixing different ratios of l -CST and d -CST, as defined by the enantiomeric excess (χ) show the transition from left-handed to pancake to right-handed bowties. Potentials of mean force obtained for preferred directions of interaction between two basic clusters of two permutations of chirality in CST ligands. g,h, LCL–LCL (g) and LCD–LCD (h) are computed and used for Monte Carlo growth simulations. The coarse-grained particle shapes deduced from the atomic structure of the nanocluster with the hydrogen bonding interaction sites marked by the spheres in LCL and LCD type nanoclusters are shown. i,j, Magnified snapshots from Monte Carlo simulations of the fully formed petal shows mismatched domains with local crystalline order arranged to form a twisted nanoribbon (j, top) and pancake type (j, bottom) petal for LCL (i) and LCD (j) type nanoclusters, respectively. Colour coding in i is only for visualization purposes. Scale bars, a, 1 µm; b, 0.5 µm; c, 200 nm; d, 100 nm; e, 2 nm−1 and f, 2 µm.

Synchrotron XRD patterns helps in understanding how the atomic structure of bowties can accommodate their continuous changes in chirality while remaining chemically stable (Fig. 1a–c). Solving for the atomic structure followed by Rietveld refinement suggested several Cd-CST polymorphs (Supplementary Fig. 5–8 and Supplementary Table 1). Cumulative analysis of the XRD models points to Cd2CST2 nanoclusters (Fig. 1j) as the structural units of nanoribbons (Supplementary Fig. 8). These units produce helical chains (Fig. 1j,k and Supplementary Fig. 8) enabling partial swapping of l -CST for d -CST without drastic disruption of the short-range coordination bonds. The amine and carboxylic groups located on the exterior of the helical chains form hydrogen bonds with neighbouring chains stabilizing the structure and facilitating the nanosheet stacking (Figs. 1j,k and 2c,d). A peak overlap below 2.5 Å in the XRD data indicates imperfections in atomic organization of nanosheets that are necessary for the chemical structure to produce a multiplicity of twisted shapes.

The reasons for variable chirality are encoded in the nanoscale building blocks. First, flexible hydrogen bonds allow for variable bond angles17,18. Second, the ability of nanoribbons to ionize (Fig. 1g,h) leads to the long-range repulsive interactions between nanoscale building blocks that can be tuned over a wide range by changing the pH and ionic strength29. Because the nanoribbons twist, the electrostatic interactions are chiral, which reinforces the handedness of their assemblies. Geometric frustrations and related mechanical deformations are contributing factors, too. In case of  nanoribbons, they leed to preferential geometries at atomic, nanometer, and submicron scales17,28,30,31, which can be captured by microscale chirality vectors, Ξ (see below). 

Selected-area electron diffraction (SAED) from a single nanoribbon confirms the XRD structural analysis (Fig. 2e). Cryo-SAED patterns are consistent with corresponding distances between planes in the orthorhombic Cd2CST2 structure (Supplementary Fig. 11). The nanoribbons are polycrystalline with ordered domains of helical chains appearing as sharp spots in the SAED pattern (Fig. 2e and Supplementary Fig. 17). Some degree of polycrystallinity enables the nanoribbons to accommodate CST ligands of opposite chirality while retaining the physical integrity of the nanoassembly. ‘Imperfections’ and variability of nanoribbon lengths is essential for the formation of such complex particles with continuously variable twisted shape.

Extinction and CD spectroscopies from ultraviolet to terahertz (THz) frequencies demonstrate a rich set of mirror symmetric bands for the left- and right-handed bowties, whereas the ‘pancake stack’ assemblies are chiroptically silent (Fig. 1d,f). Visible light, Raman, Fourier transform infrared, mid-infrared and THz spectra (Fig. 1d–f) confirm the molecular structure of the bowties established by XRD and SAED. Peaks at 510 cm−1 in Raman spectra correspond to S–S disulfide stretches (Supplementary Fig. 2a). The strong and wide band from 3,000 to 3,750 cm−1 in Fourier transform infrared spectra confirms the multiplicity of hydrogen bonds stabilizing nanoribbons and their stacks. The CD spectrum has several positive and negative peaks denoted as P1 (−, 1,040 nm), P2 (+, 460 nm), P3 (−, 330 nm), P4 (+, 270 nm) and P5 (−, 235 nm) (Fig. 1f). P5 and P4 in the ultraviolet (UV) part of the spectrum are attributed to electronic transitions in the CST ligands. These two peaks are broadened and red-shifted compared to free CST due to coordination with Cd2+. The vibrational CD (VCD) spectra of bowties and free CST are very different between 1,550 and 1,650 cm−1 (Fig. 1e). The former is 100 times more intense and contains many split bands. Density functional theory calculations (Supplementary Figs. 9 and 10) identified these bands as two carbonyl stretching modes (out-of-phase and in-phase modes) coupled with H-N-H scissor oscillations. These delocalized modes with long-range couplings between the helical chains in the nanoribbons result in large absorption cross-section and drastic increase in the intensity of 1,610–1,620 cm−1 band. The non-equivalency of the O atoms bonded to a Cd atom removes the degeneracy of the out-of-phase and in-phase modes resulting in band splitting (Figs. 1j,k and Supplementary Fig. 6). THz CD shows the existence of chiral phonons propagating along the nanoribbons (Fig. 1d)32,33. These data also indicate that the bowties made with l - or d -CST are true multiscale enantiomers of each other for which the mirror asymmetry is present at the subnanometre, nanometre and micrometre scales (the multiscale chirality vector below).

On the basis of the solved crystal structure and zeta-potential measurements (Fig. 1g,h and Supplementary Figs. 20–22), we constructed a coarse-grained version of the basic nanocluster (Supplementary Fig. 13) (Methods). A Monte Carlo growth algorithm was used to simulate bowtie assembly from nanoclusters with an average net charge of +1. Because the addition of each new nanocluster increases the total charge of the growing particle, the self-assembly process is electrostatically restricted32, with each extra nanocluster having to overcome an increasingly repulsive energy barrier. Computation of the pair-wise potentials (Fig. 2g,h) between l -homochiral and racemic versions of the nanoclusters, abbreviated to l -CSTs (LCL) and LCD, revealed growth directions along which the nanoclusters anisotropically interact through hydrogen bonds that reflect the structural hierarchy. Nanoclusters assemble into nanoplatelets that grow into nanoribbons that stack to form the bowties. As the number of building blocks in the assembly increases, the size of the ordered domains in the nanoribbons also increases (shown in colour in Fig. 2i and Supplementary Fig. 15). The growth of nanoribbons stops when the net charge repulsion on the particle becomes sufficient to prevent further attachment of nanoclusters (Supplementary Fig. 36), balancing out the attractive forces. LCL nanoclusters form a twisted platelet (Fig. 2j, top row), whereas LCD nanoclusters form a flat sheet (Fig. 2j, bottom row). The model combining both electrostatic and geometric reconstructions correctly predicts the relationship between the enantiomeric composition of the nanoclusters, χ, and the chirality of the bowties propagating up in scale as the particles evolve in size. The long-range electrostatic restrictions combined with efficacy of multiscale chirality transfer for different χ results in the uniform pitch and size of the bowties (Fig. 2f and Supplementary Fig. 31). The evidence of the elastic stress in the stacks contributing to the hierarchical twist can be appreciated by the partial unwinding of the nanoribbons on disassembly of the bowties (Supplementary Fig. 18).

Hierarchical assembly from chiral nanoclusters enables wide-range tuning of long- and short-range interactions by changing various physical–chemical parameters (Supplementary Fig. 25), which can be exemplified by three synthetic cases. In case 1, the increase of ionic strength results in the increase of l, w and t because the range of electrostatic interactions is reduced for high ionic strengths (Supplementary Fig. 20). In case 2, the zeta potential increase with pH results in the increase of l, w and t. De-protonation of the –COOH groups on the surface of nanoribbons and nanoplatelets results in the reduction of hydrogen bonding between the ribbons and, thus, in greater separation between the nanoribbons (Supplementary Figs. 21 and 22). In case 3, an increase of l -CST/Cd leads to an increase of l, w and t due to adsorption of the counter ions on the bowtie surface (Supplementary Fig. 16a). The increase of the CST concentration promotes (Fig. 3a and Supplementary Figs. 20 and 26) the close-range hydrogen bonding between the nanoclusters, resulting in thicker and wider nanoribbons (Fig. 3a and Supplementary Fig. 24).

a, Effect of NaCl, pH, excess charge (l -CST/Cd), solvent (1:1 water to solvent) on bowtie morphology (water, [Cd2+] = 1 mM, 4 mM for pH, [l -CST] = 4 mM) are shown. Scale bar, 1 µm.b, Continuously variable twist and sizes for bowtie particles obtained for different χ and [Cd2+]. Scale bar, 5 µm. c, SEM image of a typical bowtie assembly and the corresponding morphological parameters of length (l), width (w), thickness (t) and pitch (2wtan(90 − θ)) overlaid. Three-dimensional models (right) used in computations of OPD chirality measures. Scale bar, 2 µm. d, Scatter chart demonstrating the variability of the geometrical parameters (l, w, t) versus pitch (p) for the bowtie assemblies obtained in this study and the same parameters obtained in other studies16,36,37,38,39; bowties retain their shape over two orders of length scales. e, Variation of the morphological parameters with calculated OPD chirality measure show empirical correlations between the l, w, t, p and OPD.

The variability of the pitch serves as a demonstration of the tuneable chirality of these particles (Supplementary Fig. 33). The computational model describing the hierarchical self-assembly of the bowties from nanoclusters with intermediate formation of nanoplatelets and nanoribbons quantitatively predicts their geometry (Supplementary Figs. 31, 24 and 35). It also captures the variations of the geometrical parameters of the bowties in response to variations of ionic screening, solvents, enantiomeric excess, temperature and counter ions (Fig. 3d,e and Supplementary Figs. 25–33 and 39–43).

The chirality continuum can be demonstrated by calculating chirality measures, such as Hausdorff distance1, continuous symmetry measure2 or Osipov–Pickup–Dunmur (OPD) index (Supplementary Information)3. All of them are suitable for establishing the variable chirality of the bowties, but OPD has an additional useful feature: it changes sign when switching from left to right enantiomers. The transition from discrete to continuous chirality is evident from Fig. 3d,e because the pitch and OPD cover a wide range of values. For example, standard OPD varies gradually from −4.7 to +4.6, whereas other molecular or nanoscale systems have only discrete points (Fig. 3d).

Nanoassemblies have chirality at many levels, which is essential for establishing the properties of chemical structures displaying chirality continuum. Multiscale chirality can be expressed as a chirality vector Ξ(ξ1,ξ2,ξ3,ξ4...) where components ζi represent chirality measures at specific scales. For example, Ξbowties will comprise OPD for CST ligands, constituent nanoclusters, individual nanosheets (Fig. 1j,k) and bowties overall (Fig. 2a), resulting in Ξ(ξi) = Ξ(−1.1 × 10−3; 1.4 × 10−2; 6.81 × 10−4; 2.75 × 10−2); as all of them have different numbers of points, the point-scaled chirality measure OPDs instead of standard OPD was used. For some components, i.e.  ξ1, chirality continuum is not observed, but  for the level of fully assembled particles it is and is being considered henceforth.

Testing the correlation between the sign of peaks in the CD spectra (Figs.3e,4a and Supplementary Fig. 47d,e) is fundamentally important because of the long-standing controversy about the suitability of sign-changing pseudoscalar measures of chirality for the description of experimental data5,34. We found that the sign of OPD correlates with the experimental CD spectra, g-factor spectra and P3, P2 and P1 peaks (Fig. 4a and Supplementary Figs. 49 and 50)5,34. An exponential relationship between OPD and the spectral position of the P2 and P3 peaks was observed (Fig. 4c). We note that organic and biological chemistries do not provide, hitherto, examples of quantifiable property–chirality relationships, despite extensive discussions of their possibility4,5,25,35.

a, Normalized CD spectra for right-handed bowties of different OPD with peaks P1(−), P2(+), P3 (−) labelled. b, The corresponding extinction spectra. c, Dependences of spectral position of P1, P2, P3 with OPD made with varying [l -CST]/[Cd2+]. For peak P1, the dependence of its position on chirality measures is difficult to assess due to its breadth. d,e, Calculated g factor (d) and extinction (e) spectra are shown here for the right-handed twisted bowtie models with OPD values as in a and b. f, The g-factor spectra obtained for bowties prepared in an water-acetonitrile mixture with a peak at a 1,550 nm wavelength. g, Simulated g-factor spectra averaged for a bowtie model with l, w, t and θ as 5.3 ± 0.3 µm, 2.1 ± 0.1 µm, 0.7 ± 0.1 µm, 46 ± 5° along 25 different orientations with respect to an incident electromagnetic wave. h, Spatial distribution of differential (\({| E| }_{{\rm{LCP}}}^{2}-{| E<?RemoveMO1?>| }_{{\rm{RCP}}}^{2}\) ) field around a bowtie model described in g shows enhancement of field within the layers. i,j, Printed coatings of 1.5 × 1.5 cm2 l - and d -bowties on glass (i) and cloth (j) with the point-cloud of scattered signal on illumination with 1,550 nm wavelength chiral-LIDAR laser. Scale bar, 1 cm.

Experimental dependences of P2/P3 peaks versus OPD indicate that the physical laws governing chirality–property relationships for nanostructured microparticles and small molecules are different. We found that the computational models describing bowties as a stack of twisted ribbons accurately explain the signs and red-shifts of P3, P2 and P1 peaks across the entire parameter space (Fig. 4a–e and Supplementary Figs. 25, 27, 28, 39–43 and 45). Decomposition of the extinction spectra from right-handed bowties into scattering and absorptive components showed that the negative P3 and P1 peaks and positive P2 peak in Fig. 1f and Fig. 4a are from differential scattering (Supplementary Fig. 49). The computational predictions for peak position match experiments nearly perfectly (Fig. 3a,b and Supplementary Figs. 26, 27, 45 and 48–50). Electromagnetic field distributions indicate that P1 and P3 originate from a mix of dipolar and quadrupolar resonance scattering modes (Fig. 4a–h and Supplementary Fig. 54). Although the peak positions vary with particle size, the attribution of the peaks remains the same for all l, w, t and θ, as demonstrated by the calculations for the different models (Supplementary Figs. 49 and 50). Note that both scattering and absorption peaks of bowties involve the interaction of photons with the particle as a whole, not specific quantum states, which is essential for understanding the origin of observed chirality–property relationships (Fig. 4a–e). Optical properties of molecules originate from chiral shapes of molecular orbitals that may or may not correlate with the geometry of the sp3 carbon atom or Δ/Λ enantiomers. Therefore, the similarity of the physical dimensions of chiral objects in bowties and incident photons engenders the observed dependences (Fig. 4c) and, thus, the fundamental difference in chirality–property relationships in comparison with small molecules, polymers and so on.

The simplicity of bowtie synthesis and the predictive relationship between particle structure and optical properties opens the door to the engineering of particle-based metasurfaces. We demonstrate that bowties can be used as templates for inducing chiroptical emission from achiral perovskite nanocrystals (Supplementary Fig. 51). Coatings with circular polarization signatures (Supplementary Fig. 53) are needed for machine vision systems to increase the accuracy of object identification (for example, pedestrians). Dispersions of bowtie particles with strong circularly polarized P1 peak scattering in the near infrared (NIR) part of the spectrum were printed on glass slides and cotton fabric. A LIDAR operating at 1,550 nm with a circularly polarized beam (Supplementary Fig. 52) was used to acquire images of the printed objects following their decomposition into left- and right-handed components of back-scattered light originating from the P1 peak in CD spectra (Fig. 1f). The distinct contrast between prints made from the left- and right-handed bowties was observed (Fig. 4g and Supplementary Figs. 54–56). The layer of left-handed bowties can also be identified with high contrast on cotton cloth in LIDAR images (Supplementary Fig. 53d,e), which can be used in autonomous vehicles and other robots for the identification of specific materials and surface features exemplified by the polarization contrast from the micrometre-scale chirality of the fabrics and bowties.

Dispersions of bowties translate into an emerging opportunity for scalable manufacturing of printable metasurfaces. The observed continuous correlations between the position of CD peaks and OPD indicate chirality–property correlations with pseudoscalar chirality measures arise when simple chiral shapes continuously transition from left- to right-handed enantiomers through an achiral state of flat sheets by gradual twisting. Conversion of Cd-CST bowties into similarly shaped CdS particles and coating of bowties with a perovskite nanolayer (Supplementary Fig. 51) indicates that several other materials can be prepared with the same diversity of morphology (Supplementary Fig. 12). Besides machine vision, dipolar and quadrupolar photonic resonances enable chirality-based engineering of composite materials for chemical and biomedical technologies.

Cadmium chloride (CdCl2), 99.99% trace metals basis, l -CST, ≥98% (TLC), crystalline, d -CST, 98%, sodium hydroxide (NaOH) pellets >97% were purchased from Sigma-Aldrich. De-ionized water (18.2 mΩ cm−1) was used for the preparation of stock solutions and aqueous dispersions. The 10 ml stock solutions of CdCl2 (0.1 M), l -CST (0.1 M), d -CST (0.1 M) and NaOH (2.5 M) were prepared by dissolving the required amounts in de-ionized water. pH of CST stock solution was raised to 11 by adding 1 ml of 2.5 NaOH in 10 ml of solution. An immediate change in appearance was observed by the formation of a clear solution. Remaining solids were dissolved by mild sonication for 10 s. Bowties were synthesized by the mixing stock solutions of CdCl2 and l -CST or d -CST in a 1:1 stoichiometric ratio in an aqueous media. Typical synthesis involved 20 µl of l -CST added to 960 µl of water followed by the addition of 20 µl of CdCl2. The solution mixture was shaken vigorously until the dispersion turned milky, which is an indicator of the bowties’ assembly process. After that, the mixture was kept still at room temperature for 15 min to ensure the completion of the self-assembly. The dispersion was subsequently centrifuged three times in de-ionized water at 6,000 rpm for 3 min. Final aqueous dispersion was stored at room temperature and used for further characterization and studies. Bowtie dispersions for coatings were prepared by mixing the freeze-dried powder in polyacrylic acid.

SEM samples were prepared by drop-casting 5 µl of aqueous dispersion on 1 × 1 cm silicon wafers (TedPella), following by drying at room temperature. This was followed by sputter coating a 5–10 nm film of gold on the wafer to avoid charging due to the electron beam. SEM measurements were performed in FEI Nova NanoLab Dual Beam SEM and FEI Helios Nanolab at 5 kV accelerating voltage and 0.4 nA beam current under secondary electron detection mode.

TEM samples were prepared by drop-casting 10 µl of aqueous solution on to a copper grid coated with holey carbon supported by a continuous carbon film (TedPella 01824). Bright-field TEM was performed on JEOL 2010 operating at 300 kV accelerating voltage and acquisition were done using Gatan OneView camera. Annular dark-field–scanning transmission electron microscopy (ADF–STEM) and STEM–energy dispersive X-ray spectroscopy measurements were performed on cold-FEG JEOL 3100R05 with Cs aberration correction operating at 300 keV. A HAADF detector was used to acquire Z-contrast images in which the intensity is proportional to the atomic number of the column over which the electron probe is placed. Diffraction experiments were acquired using Thermo Fisher Talos F200X operated at 200 keV equipped with a Gatan OneView camera. The Elsa (698) Gatan Cryo Holder cooled specimens down to roughly 93 K for low-temperature measurements.

The synthesized bowtie particles were dispersed in water and drop cast using a micropipet onto a 3 mm copper TEM grid dried at room temperature. The TEM grid was an ultrathin (3 nm) carbon film with a large hexagonal mesh (100) to provide high specimen tilts without beam shadowing (Electron Microscopy Sciences). The tomographic tilt series were acquired using a Thermo Fisher Talos F200X operated at 200 keV as a scanning transmission electron microscope with a 10.5 mrad semiconvergence angle using a 36 and 165 mrad inner and outer semicollection angles for the ADF detector. The right-handed bowtie particles were acquired over a tilt range of −72º to +73º and left-handed bowtie particles specimen were acquired over a tilt range of −75º to +71º both with a +1º tilt increment. At each angle, ADF images of size 1,024 × 1,024 pixels were recorded with a dwell time of 4 µs and pixel size of 4.94 nm. The tomograms were reconstructed with the additive simultaneous iterative reconstruction technique for 150 iterations. The three-dimensional reconstructions were visualized by tomviz40.

Synchrotron X-ray powder diffraction data (λ = 0.45192 Å) were collected at Beamline 17-BM at the Advanced Photon Source at the Argonne National laboratory. A freeze-dried powder sample of l -/d -/rac-bowties was loaded in a 1.0 mm diameter Kapton capillary tube, which was then measured using a VAREX XRD 4343CT amorphous silicon area detector in the Debye–Scherrer geometry.

Length (l), width (w), thickness (h) and angle of twist (θ) as defined in Fig. 3c were measured using Jann5s or the measure tool from MATLAB central. A minimum of ten measurements were taken for each parameter and standard deviation (σ) was calculated and represented as percentage variation by following the equation \({\rm{e}}{\rm{r}}{\rm{r}}{\rm{o}}{\rm{r}}=\frac{{\sigma }_{X}}{\mathop{x}\limits_{\_}}\times 100\) , where x = l, w, h, θ and x is the mean of corresponding values.

The CD and g-factor spectra were calculated using a software package (Lumerical Solutions Inc.; www.lumerical.com/tcad-products/fdtd/). We used total-field scattered-field sources that surrounded the structure being modelled. Circularly polarized light  was generated by positioning two total-field scattered-field sources along the same forward axis at a 90° angle and with a phase difference of either −90° for photons with left-handed polarization or 90° for photons with right-handed polarization. Two analysis groups using box power monitors monitored the absorption and scattering cross sections (extinction is the sum of absorption and scattering). The FDTD simulation region was defined by a larger box monitor with a stretched-coordinate perfectly matched layer and non-uniform mesh type. Frequency profile monitors were inserted in the total-field region to calculate electric field distribution in 3D. The refractive index for water is 1.33. Convergence tests with different mesh sizes were performed to determine the best balance between computational time restraints and simulation accuracy. Simulations were carried out for the bowties oriented with their long axis being parallel with the propagation direction (k vector) of photons. The 10 and 1 nm mesh sizes produced similar CD spectra; therefore, we used 10-nm mesh size for bowtie simulations presented in Supplementary Figs. 49, 50 and 54.

CD measurements in UV-visible and NIR range were performed using JASCO J-1700 equipped with one photomultiplier tube detector in 200–800 nm range and two InGaAs NIR detectors in 800–1,600 nm range. Typical scanning parameters were as follows: scanning speed, 500 nm min−1; data pitch, 0.1 nm; bandwidth, 5 nm (NIR bandwidth, 10 nm), digital integration time, 0.25 s and one accumulation. The anisotropy g factor was calculated according to the equation \(g=\frac{{\rm{C}}{\rm{D}}}{32,980\times \,{\rm{a}}{\rm{b}}{\rm{s}}}\) , where CD is the signal obtained from CD/DC channel and abs is the absorbance calculated from ABS and DC channels of the spectrophotometer. CD spectra were stopped at wavelengths in which the high-tension voltage was greater than 800 V to avoid artefacts during acquisition.

The terahertz (THz) time-domain polarimetry system based on three linear polarizers was used to measure THz absorption and THz CD spectra7. We used calculations of Stokes parameters based on the Ex and Ey, the electric field in the x and y directions, respectively. THz absorptions were extracted from the equations used in ref. 7 after retrieval of phase from the transmission data. To eliminate linear birefringence effects from samples, highly concentrated mixture of bowtie powders with mineral oil was used. The reference sample is a quartz sandwich cell filled with mineral oil32.

VCD measurements were performed on freeze-dried bowtie samples dispersed in heavy water (D2O) at 33 mg ml−1 concentration. A 100 µl drop was sandwiched between two BaF2 crystals separated by a 50 µm Teflon spacer. An MCT-V detector was used to acquire infrared and VCD data in the range of 2,000 to 850 cm−1 with a resolution of 4 cm−1 and a total of 100 and 500 accumulations, respectively. The sandwiched dispersion between BaF2 crystals was rotated along the a axis coinciding with the direction of the beam at a constant speed to avoid settling of particles. Corresponding infrared and VCD were plotted as A and \(\frac{\triangle \,A}{A}\) , respectively, with exclusion of the 1,300 to 1,100 cm−1 range that corresponds to strong absorption from D2O. The simulation of the VCD spectra was performed using the method described in the Supplementary Information.

We pick features consistent with XRD results to construct a coarse-grained building block for the nanoclusters. The key features we chose to adapt are (Supplementary Fig. 13) as follows. (1) Flat, rectangular prism-like building blocks with rough dimensions of 20 × 10 × 5 commensurate with XRD nanoclusters. (2) H-bonding sites decorate the two largest and two smallest facets to reflect the bonding between the –COOH and –NH2 groups. (3) The two intermediate facets show a directional interaction to capture the disulfide bridge formation intrinsic to CST. These interactions sites only interact with themselves. (4) Building blocks show a slight shear on one facet to enforce chirality of l -CST versus R- CST. (5) Isotropic charge–charge repulsion between Cd2+ atoms embedded inside the building block: each has a net charge of +1 to build in charge-driven, self-limiting behaviours. (6) All other surface sites not defined above show steric repulsion.

Potential of mean force (PMF) values between nanoclusters are computed by adding together all interactions as a function of the centre-to-centre distance between nanocluster pairs and performing weighted averaging. The resulting PMFs between LCL highlight directional growth features that are commensurate with the hierarchy of structural formation observed in experiments (Fig. 2g). PMFs computed for a mix of l -CST and r -CST (LCD) conversely suggest a reduction in hierarchical structural formation and chiral propagation (Fig. 2h). We then use the computed PMFs to perform Monte Carlo growth simulations. Specifically, the well depth of each computed PMF plus the net charge–charge repulsion between the next nanocluster and growing structure are used as the input energy in the Metropolis algorithm. This enables Boltzmann weighted selection of the direction in which to place the next building block in the growing structure. Visualization of the growth simulation results confirm the formation of bowtie and pancake morphologies (Fig. 2i,j).

The authors declare that the data supporting the findings of this study are available within the article and its supplementary information files. Source data are provided with this paper.

Buda, A. B. & Mislow, K. A Hausdorff chirality measure. J. Am. Chem. Soc. 114, 6006–6012 (1992).

Osipov, M. A., Pickup, B. T. & Dunmur, D. A. A new twist to molecular chirality: intrinsic chirality indices. Molec. Phys. 84, 1193–1206 (1995).

Article  CAS  ADS  Google Scholar 

Zabrodsky, H. & Avnir, D. Continuous symmetry measures. 4. Chirality. J. Am. Chem. Soc. 117, 462–473 (1995).

Rassat, A. & Fowler, P. W. Is there a ‘most chiral tetrahedron’? Chemistry 10, 6575–6580 (2004).

Article  CAS  PubMed  Google Scholar 

Harris, A. B., Kamien, R. D. & Lubensky, T. C. Molecular chirality and chiral parameters. Rev. Mod. Phys. 71, 1745–1757 (1999).

Article  CAS  ADS  Google Scholar 

Xia, Y. et al. Self-assembly of self-limiting monodisperse supraparticles from polydisperse nanoparticles. Nat. Nanotechnol. 6, 580–587 (2011).

Article  CAS  PubMed  ADS  Google Scholar 

Choi, W. J. et al. Terahertz circular dichroism spectroscopy of biomaterials enabled by kirigami polarization modulators. Nat. Mater. 18, 820–826 (2019).

Article  CAS  PubMed  ADS  Google Scholar 

Nguyen, M.-K. & Kuzyk, A. Reconfigurable chiral plasmonics beyond single chiral centers. ACS Nano 13, 13615–13619 (2019).

Article  CAS  PubMed  Google Scholar 

Dey, S. et al. DNA origami. Nat. Rev. Methods Prim. 1, 13 (2021).

Probst, P. T. et al. Mechano-tunable chiral metasurfaces via colloidal assembly. Nat. Mater. 20, 1024–1028 (2021).

Article  CAS  PubMed  ADS  Google Scholar 

Kim, Y. et al. Reconfigurable chiroptical nanocomposites with chirality transfer from the macro- to the nanoscale. Nat. Mater. 15, 461–468 (2016).

Article  CAS  PubMed  ADS  Google Scholar 

Miyajima, D. et al. Liquid crystalline corannulene responsive to electric field. J. Am. Chem. Soc. 131, 44–45 (2009).

Article  CAS  PubMed  Google Scholar 

Ohta, E. et al. Redox-responsive molecular helices with highly condensed π-clouds. Nat. Chem. 3, 68–73 (2011).

Article  CAS  PubMed  Google Scholar 

Ma, W. et al. Chiral plasmonics of self-assembled nanorod dimers. Sci. Rep. 3, 1934 (2013).

Article  PubMed  PubMed Central  Google Scholar 

Yashima, E. et al. Supramolecular helical systems: helical assemblies of small molecules, foldamers, and polymers with chiral amplification and their functions. Chem. Rev. 116, 13752–13990 (2016).

Article  CAS  PubMed  Google Scholar 

Lee, H. E. et al. Amino-acid- and peptide-directed synthesis of chiral plasmonic gold nanoparticles. Nature 556, 360–364 (2018).

Article  CAS  PubMed  ADS  Google Scholar 

Jiang, W. et al. Emergence of complexity in hierarchically organized chiral particles. Science 368, 642–648 (2020).

Article  CAS  PubMed  ADS  Google Scholar 

Medina, D. D. & Mastai, Y. Chiral polymers and polymeric particles for enantioselective crystallization. Isr. J. Chem. 58, 1330–1337 (2018).

Gingras, M. One hundred years of helicene chemistry. Part 3: applications and properties of carbohelicenes. Chem. Soc. Rev. 42, 1051–1095 (2013).

Article  CAS  PubMed  Google Scholar 

Ma, W. et al. Chiral inorganic nanostructures. Chem. Rev. 117, 8041–8093 (2017).

Article  CAS  PubMed  Google Scholar 

Aida, T., Meijer, E. W. & Stupp, S. I. Functional supramolecular polymers. Science 335, 813–817 (2012).

Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

Ho, R.-M. et al. Block copolymers with a twist. J. Am. Chem. Soc. 131, 18533–18542 (2009).

Article  CAS  PubMed  Google Scholar 

Wang, H.-F. et al. Generalizing the effects of chirality on block copolymer assembly. Proc. Natl Acad. Sci. USA 116, 4080 LP–4084089 (2019).

Ido, S. et al. Beyond the helix pitch: direct visualization of native DNA in aqueous solution. ACS Nano 7, 1817–1822 (2013).

Article  CAS  PubMed  Google Scholar 

Petitjean, M. Chirality and symmetry measures: a transdisciplinary review. Entropy https://doi.org/10.3390/e5030271 (2003).

Buda, A. B., der Heyde, T. A. & Mislow, K. On quantifying chirality. Angew. Chemie Int. Ed. English 31, 989–1007 (1992).

Kotov, N. A., Fendler, J. H. & Dekany, I. Layer-by-layer self-assembly of polyelectrolyte- semiconductor nanoparticle composite films. J. Phys. Chem. 99, 13065–13069 (1995).

Tang, Z., Zhang, Z., Wang, Y., Glotzer, S. C. & Kotov, N. A. Self-assembly of CdTe nanocrystals into free-floating sheets. Science 314, 274–278 (2006).

Article  CAS  PubMed  ADS  Google Scholar 

Nguyen, T. D., Schultz, B. A., Kotov, N. A. & Glotzer, S. C. Generic, phenomenological, on-the-fly renormalized repulsion model for self-limited organization of terminal supraparticle assemblies. Proc. Natl Acad. Sci. USA 112, E3161–E3168 (2015).

Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

Gao, C. et al. Electrostatic shape control of a charged molecular membrane from ribbon to scroll. Proc. Natl Acad. Sci. USA 116, 22030–22036 (2019).

Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

Miller, J. M. et al. Conformational switching of chiral colloidal rafts regulates raft–raft attractions and repulsions. Proc. Natl Acad. Sci. USA 116, 15792–15801 (2019).

Article  CAS  PubMed  PubMed Central  ADS  Google Scholar 

Choi, W. J. et al. Chiral phonons in microcrystals and nanofibrils of biomolecules. Nat. Photonics. 16, 366-373 (2022).

Yeom, J. et al. Chiromagnetic nanoparticles and gels. Science 359, 309–314 (2018).

Millar, G., Weinberg, N. & Mislow, K. On the Osipov–Pickup–Dunmur chirality index: why pseudoscalar functions are generally unsuitable to quantify chirality. Mol. Phys. 103, 2769–2772 (2005).

Article  CAS  ADS  Google Scholar 

Pelayo, J. J., Whetten, R. L. & Garzón, I. L. Geometric quantification of chirality in ligand-protected metal clusters. J. Phys. Chem. C. 119, 28666–28678 (2015).

Shaw, L. A. et al. Scanning two-photon continuous flow lithography for synthesis of high-resolution 3D microparticles: erratum. Opt. Express 26, 14718 (2018).

Article  PubMed  ADS  Google Scholar 

Chang, Y.-H., Jang, J.-W., Chang, Y.-C., Lee, S.-H. & Siao, T.-F. Gold nanohelices: a new synthesis route, characterization, and plasmonic E-field enhancement. ACS Omega 5, 14860–14867 (2020).

Article  CAS  PubMed  PubMed Central  Google Scholar 

Kaschke, J. & Wegener, M. Gold triple-helix mid-infrared metamaterial by STED-inspired laser lithography. Opt. Lett. 40, 3986 (2015).

Article  CAS  PubMed  ADS  Google Scholar 

Yan, J. et al. Self-assembly of chiral nanoparticles into semiconductor helices with tunable near-infrared optical activity. Chem. Mater. 32, 476–488 (2020).

Levin, B. D. A. et al. Tutorial on the visualization of volumetric data using tomviz. Micros. Today 26, 12–17 (2018).

This project was supported by the Vannevar Bush DoD Fellowship to N.A.K. titled ‘Engineered Chiral Ceramics’ Office of Naval Research (ONR), grant no. N000141812876, National Science Foundation (NSF) project ‘Energy- and Cost-Efficient Manufacturing Employing Nanoparticles’, grant nos. NSF 1463474, ONR N000141812876; ONR COVID-19 Newton Award ‘Pathways to Complexity’ with ‘Imperfect’ Nanoparticles’, grant no. HQ00342010033, NSF 1566460 ‘Nanospiked Particles for Photocatalysis’; NSF grant no. 1538180 and NSF grant no. DMR-9871177 for funding the JEOL 2010F analytical electron microscope used in this work. N.A.K., P.K. M.V. and J.S. are also thankful to the Office of Naval Research for the support from a MURI project  Mechanical Metamaterials, ONR grant no. N00014-20-1-2479. Support for the Dual Source and Environmental X-ray Scattering Facility at the University of Pennsylvania was provided by the Laboratory for Structure and Matter, which is funded in part by NSF grant no. MRSEC 1720530. P.K., T.V., S.G., E.M., C.M. and N.A.K. were supported by Office of Naval Research Multidisciplinary University Research Initiative Award no. ONR N00014-18-1-2497; R.H. acknowledges support from Army Research Office (ARO), grant no. W911NF-22-1-0056. This research was supported in part through computational resources and services provided by Advanced Research Computing, a division of Information and Technology Services at the University of Michigan, Ann Arbor, USA. V.P.N. acknowledges funding from Unitatea Executiva pentru Finantarea Invatamantului Superior, a Cercetarii, Dezvoltarii si Inovarii (UEFISCDI), grant no. PN-III-P4-ID-PCE-2020-2783-P, contract no. PCE 2/2022. Initial experiments by J. Yan and Z. Mu are gratefully acknowledged39. We are also grateful to G. Grason, A. Hunter and B. Palen for helpful discussions.

Department of Chemical Engineering, University of Michigan, Ann Arbor, MI, USA

Prashant Kumar, Thi Vo, Minjeong Cha, Anastasia Visheratina, Ji-Young Kim, Alexander Simon, Daniel Katz, Won Jin Choi, Michael Veksler, Sharon Glotzer & Nicholas A. Kotov

Biointerfaces Institute, University of Michigan, Ann Arbor, MI, USA

Prashant Kumar, Thi Vo, Minjeong Cha, Anastasia Visheratina, Ji-Young Kim, Won Jin Choi, Sharon Glotzer & Nicholas A. Kotov

X-ray Science Division, Advanced Photon Source, Argonne National Laboratory, Lemont, IL, USA

Wenqian Xu & Si Chen

Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI, USA

Jonathan Schwartz, Robert Hovden, Sharon Glotzer & Nicholas A. Kotov

Pro-Vitam Ltd, Sfantu Gheorghe, Romania

Department of Chemistry, University of Pennsylvania, Philadelphia, PA, USA

Emanuele Marino & Christopher Murray

Department of Physics and Chemistry, University of Palermo, Palermo, Italy

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

N.A.K. designed and directed all aspects of the project. P.K. conducted the synthesis of libraries of bowties with assistance from A.S. and D.K. P.K. measured VCD and CD data. W.J.C. measured the terahertz CD spectroscopy data. P.K., W.X. and S.C. acquired, solved and refined the crystal structure from XRD data with input from N.A.K. P.K. conducted all the SEM imaging. J.S., P.K. and R.H. conducted the TEM and tomography measurements. T.V. and S.G. conducted the Monte Carlo growth simulations and theoretical model development for structure prediction with input from P.K. and N.A.K. J.-Y.K. and M.C. wrote the code for calculation of OPD indices and J.-Y.K. calculated all the values. A.V. built 3D models of bowties based on geometrical parameters observed in SEM images for OPD calculation. V.P.N. and P.K. conducted the VCD simulations. P.K. simulated the chiroptical properties with input from N.A.K. E.M. and C.M. measured single particle scattering data. M.V. developed perovskite coated bowties. M.C. and N.A.K. conducted the LIDAR measurement and analysis. All the authors discussed the results. P.K. and N.A.K. wrote the paper with written contributions from all the authors. N.A.K. put together chirality vectors Ξ.

Correspondence to Sharon Glotzer or Nicholas A. Kotov.

N.A.K. is a founder of a start-up company Photon Semantics developing LIDAR technologies for machine vision. The other authors declare no competing interests.

Nature thanks Tobias König and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

This file contains Supplementary Methods, Figs. 1–56, Tables 1–3 and References.

Python code for OPD calculations.

Slice and view of the bowtie particle.

Tomogram of the bowtie particle.

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Kumar, P., Vo, T., Cha, M. et al. Photonically active bowtie nanoassemblies with chirality continuum. Nature 615, 418–424 (2023). https://doi.org/10.1038/s41586-023-05733-1

DOI: https://doi.org/10.1038/s41586-023-05733-1

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Nature (Nature) ISSN 1476-4687 (online) ISSN 0028-0836 (print)

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.