Skip to main content

Thank you for visiting 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.

Thermally reconfigurable monoclinic nematic colloidal fluids


Fundamental relationships are believed to exist between the symmetries of building blocks and the condensed matter phases that they form1. For example, constituent molecular and colloidal rods and disks impart their uniaxial symmetry onto nematic liquid crystals, such as those used in displays1,2. Low-symmetry organizations could form in mixtures of rods and disks3,4,5, but entropy tends to phase-separate them at the molecular and colloidal scales, whereas strong elasticity-mediated interactions drive the formation of chains and crystals in nematic colloids6,7,8,9,10,11. To have a structure with few or no symmetry operations apart from trivial ones has so far been demonstrated to be a property of solids alone1, but not of their fully fluid condensed matter counterparts, even though such symmetries have been considered theoretically12,13,14,15 and observed in magnetic colloids16. Here we show that dispersing highly anisotropic charged colloidal disks in a nematic host composed of molecular rods provides a platform for observing many low-symmetry phases. Depending on the temperature, concentration and surface charge of the disks, we find nematic, smectic and columnar organizations with symmetries ranging from uniaxial1,2 to orthorhombic17,18,19,20,21 and monoclinic12,13,14,15. With increasing temperature, we observe unusual transitions from less- to more-ordered states and re-entrant22 phases. Most importantly, we demonstrate the presence of reconfigurable monoclinic colloidal nematic order, as well as the possibility of thermal and magnetic control of low-symmetry self-assembly2,23,24. Our experimental findings are supported by theoretical modelling of the colloidal interactions between disks in the nematic host and may provide a route towards realizing many low-symmetry condensed matter phases in systems with building blocks of dissimilar shapes and sizes, as well as their technological applications.

Access options

Rent or Buy article

Get time limited or full article access on ReadCube.


All prices are NET prices.

Fig. 1: Colloidal disks in a nematic host of molecular rods.
Fig. 2: Phase diagrams and characterization of phases.
Fig. 3: Colloidal orientational order within an orthorhombic nematic phase.
Fig. 4: Colloidal orientational order within a monoclinic nematic phase.
Fig. 5: Order within monoclinic smectic and columnar states.

Data availability

Additional relevant data, such as details of the experimental phase diagram, are provided in the Supplementary Information. Additional raw data that further support the findings of this study are available from the corresponding author on reasonable request. Source data are provided with this paper.

Code availability

The custom codes used to generate plots and theoretical phase diagrams, as well as the Landau–de Gennes modelling and polarizing microscopy image simulation codes, are all provided in the Supplementary Information.


  1. 1.

    Chaikin, P. M. & Lubensky, T. C. Principles of Condensed Matter Physics (Cambridge Univ. Press, 1995).

  2. 2.

    de Gennes, P. G. & Prost, J. The Physics of Liquid Crystals 2nd edn (Clarendon Press, 1993).

  3. 3.

    Alben, R. Liquid crystal phase transitions in mixtures of rodlike and platelike molecules. J. Chem. Phys. 59, 4299–4304 (1973).

    CAS  ADS  Google Scholar 

  4. 4.

    van der Kooij, F. M. & Lekkerkerker, H. N. W. Liquid-crystal phases formed in mixed suspensions of rod- and platelike colloids. Langmuir 16, 10144–10149 (2000).

    Google Scholar 

  5. 5.

    Berardi, R. & Zannoni, C. Low–temperature biaxial nematic from rod and disc mesogen mixture. Soft Matter 8, 2017–2025 (2012).

    CAS  ADS  Google Scholar 

  6. 6.

    Poulin, P., Stark, H., Lubensky, T. C. & Weitz, D. A. Novel colloidal interactions in anisotropic fluids. Science 275, 1770–1773 (1997).

    CAS  PubMed  Google Scholar 

  7. 7.

    Muševič, I., Škarabot, M., Tkalec, U., Ravnik, M. & Žumer, S. Two-dimensional nematic colloidal crystals self-assembled by topological defects. Science 313, 954–958 (2006).

    PubMed  ADS  Google Scholar 

  8. 8.

    Silvestre, N. M., Patrício, P., Tasinkevych, M., Andrienko, D. & Telo da Gama, M. M. Colloidal discs in nematic liquid crystals. J. Phys. Condens. Matter 16, S1921–S1930 (2004).

    CAS  ADS  Google Scholar 

  9. 9.

    Smalyukh, I. I. Liquid crystal colloids. Annu. Rev. Condens. Matter Phys. 9, 207–226 (2018).

    ADS  Google Scholar 

  10. 10.

    Mundoor, H. et al. Electrostatically controlled surface boundary conditions in nematic liquid crystals and colloids. Sci. Adv. 5, eaax4257 (2019).

    CAS  PubMed  PubMed Central  ADS  Google Scholar 

  11. 11.

    Mertelj, A., Lisjak, D., Drofenik, M. & Čopič, M. Ferromagnetism in suspensions of magnetic platelets in liquid crystal. Nature 504, 237–241 (2013).

    CAS  PubMed  ADS  Google Scholar 

  12. 12.

    Lubensky, T. C. & Radzihovsky, L. Theory of bent-core liquid-crystal phases and phase transitions. Phys. Rev. E 66, 031704 (2002).

    MathSciNet  CAS  ADS  Google Scholar 

  13. 13.

    Luckhurst, G. R., Naemura, S., Sluckin, T. J., To, T. B. T. & Turzi, S. Molecular field theory for biaxial nematic liquid crystals composed of molecules with C2h point group symmetry. Phys. Rev. E 84, 011704 (2011).

    ADS  Google Scholar 

  14. 14.

    Mettout, B. Macroscopic and molecular symmetries of unconventional nematic phases. Phys. Rev. E 74, 041701 (2006).

    MathSciNet  CAS  ADS  Google Scholar 

  15. 15.

    Luckhurst, G. R. & Sluckin, T. J. Biaxial Nematic Liquid Crystals: Theory, Simulation and Experiment (Wiley, 2015).

  16. 16.

    Liu, Q., Ackerman, P. J., Lubensky, T. C. & Smalyukh, I. I. Biaxial ferromagnetic liquid crystal colloids. Proc. Natl Acad. Sci. USA 113, 10479–10484 (2016).

    CAS  PubMed  ADS  Google Scholar 

  17. 17.

    Freiser, M. J. Ordered states of a nematic liquid. Phys. Rev. Lett. 24, 1041–1043 (1970).

    CAS  ADS  Google Scholar 

  18. 18.

    Yu, L. J. & Saupe, A. Observation of a biaxial nematic phase in potassium laurate-1-decanol-water mixtures. Phys. Rev. Lett. 45, 1000–1003 (1980).

    CAS  ADS  Google Scholar 

  19. 19.

    Severing, K. & Saalwächter, K. Biaxial nematic phase in a thermotropic liquid-crystalline side-chain polymer. Phys. Rev. Lett. 92, 125501 (2004).

    PubMed  ADS  Google Scholar 

  20. 20.

    van den Pol, E., Petukhov, A. V., Thies-Weesie, D. M. E., Byelov, D. V. & Vroege, G. J. Experimental realization of biaxial liquid crystal phases in colloidal dispersions of boardlike particles. Phys. Rev. Lett. 103, 258301 (2009).

    PubMed  ADS  Google Scholar 

  21. 21.

    Mundoor, H., Park, S., Senyuk, B., Wensink, H. H. & Smalyukh, I. I. Hybrid molecular-colloidal liquid crystals. Science 360, 768–771 (2018).

    MathSciNet  CAS  PubMed  MATH  Google Scholar 

  22. 22.

    Cladis, P. E. New liquid-crystal phase diagram. Phys. Rev. Lett. 35, 48–51 (1975).

    CAS  ADS  Google Scholar 

  23. 23.

    Jákli, A., Krüerke, D., Sawade, H. & Heppke, G. Evidence for triclinic symmetry in smectic liquid crystals of bent-shape molecules. Phys. Rev. Lett. 86, 5715–5718 (2001).

    PubMed  ADS  Google Scholar 

  24. 24.

    Vaupotič, N. et al. Evidence for general tilt columnar liquid crystalline phase. Soft Matter 5, 2281–2285 (2009).

    ADS  Google Scholar 

  25. 25.

    Yang, J. et al. One-step hydrothermal synthesis of carboxyl-functionalised upconversion phosphors for bioapplications. Chemistry 18, 13642–13650 (2012).

    CAS  PubMed  Google Scholar 

  26. 26.

    Graf, C., Vossen, D. L. J., Imhof, A. & van Blaaderen, A. A. General method to coat colloidal particles with silica. Langmuir 19, 6693–6700 (2003).

    CAS  Google Scholar 

  27. 27.

    Jerôme, B. Surface effects and anchoring in liquid crystals. Rep. Prog. Phys. 54, 391 (1991).

    ADS  Google Scholar 

  28. 28.

    Helbing, D., Farkas, I. J. & Viscek, T. Freezing by heating in a driven mesoscopic system. Phys. Rev. Lett. 84, 1240–1243 (2000).

    CAS  PubMed  ADS  Google Scholar 

  29. 29.

    Mamlok, L., Malthete, J., Tinh, N. H., Destrade, C. & Levelut, A. M. Une nouvelle mésophase en colonnes. J. Phys. Lett. 43, 641–647 (1982).

    CAS  Google Scholar 

  30. 30.

    Lydon, J. Chromonic review. J. Mater. Chem. 20, 10071–10099 (2010).

    CAS  Google Scholar 

  31. 31.

    He, M. et al. Colloidal diamond. Nature 585, 524–529 (2020).

    CAS  PubMed  ADS  Google Scholar 

  32. 32.

    Moynihan, H. A. & O’Hare, I. P. Spectroscopic characterisation of the monoclinic and orthorhombic forms of paracetamol. Int. J. Pharm. 247, 179–185 (2002).

    CAS  PubMed  Google Scholar 

  33. 33.

    Senyuk, B., Glugla, D. & Smalyukh, I. I. Rotational and translational diffusion of anisotropic gold nanoparticles in liquid crystals controlled by varying surface anchoring. Phys. Rev. E 88, 062507 (2013).

    ADS  Google Scholar 

  34. 34.

    Loudet, C., Hanusse, P. & Poulin, P. Stokes drag on a sphere in a nematic liquid crystal. Science 306, 1525 (2004).

    CAS  PubMed  Google Scholar 

  35. 35.

    Conkey, D. B., Trivedi, R. P., Pavani, S. R. P., Smalyukh, I. I. & Piestun, R. Three-dimensional parallel particle manipulation and tracking by integrating holographic optical tweezers and engineered point spread functions. Opt. Express 19, 3835–3842 (2011).

    CAS  PubMed  ADS  Google Scholar 

  36. 36.

    Liu, Q., Yuan, Y. & Smalyukh, I. I. Electrically and optically tunable plasmonic guest-host liquid crystals with long-range ordered nanoparticles. Nano Lett. 14, 4071–4077 (2014).

    CAS  PubMed  ADS  Google Scholar 

  37. 37.

    Mundoor, H. & Smalyukh, I. I. Mesostructured composite materials with electrically tunable upconverting properties. Small 11, 5572–5580 (2015).

    CAS  PubMed  Google Scholar 

  38. 38.

    Evans, J. S., Beier, C. & Smalyukh, I. I. Alignment of high-aspect ratio colloidal gold nanoplatelets in nematic liquid crystals. J. Appl. Phys. 110, 033535 (2011).

    ADS  Google Scholar 

  39. 39.

    Shah, R. R. & Abbott, N. L. Coupling of the orientations of liquid crystals to electrical double layers formed by the dissociation of surface-immobilized salts. J. Phys. Chem. B 105, 4936–4950 (2001).

    CAS  Google Scholar 

  40. 40.

    Smalyukh, I. I., Kaputa, D. S., Kachynski, A. V., Kuzmin, A. N. & Prasad, P. N. Optical trapping of director structures and defects in liquid crystals using laser tweezers. Opt. Express 15, 4359 (2007).

    CAS  PubMed  ADS  Google Scholar 

  41. 41.

    Zhou, J. et al. Ultrasensitive polarized up-conversion of Tm3+−Yb3+ doped β-NaYF4 single nanorod. Nano Lett. 13, 2241–2246 (2013).

    CAS  PubMed  ADS  Google Scholar 

  42. 42.

    Yeh, P. & Gu, C. Optics of Liquid Crystal Displays (Wiley, 1999).

  43. 43.

    Weeks, E. R., Crocker, J. C., Levitt, A. C., Schofield, A. & Weitz, D. A. Three-dimensional direct imaging of structural relaxation near the colloidal glass transition. Science 287, 627–631 (2000).

    CAS  PubMed  ADS  Google Scholar 

  44. 44.

    Anderson, V. J. & Lekkerkerker, H. N. W. Insights into phase transition kinetics from colloid science. Nature 416, 811–815 (2002).

    CAS  PubMed  ADS  Google Scholar 

  45. 45.

    Gasser, U., Weeks, E. R., Schofield, A., Pusey, P. N. & Weitz, D. A. Real-space imaging of nucleation and growth in colloidal crystallization. Science 292, 258–262 (2001).

    CAS  PubMed  ADS  Google Scholar 

  46. 46.

    Besseling, T. H. et al. Determination of the positions and orientations of concentrated rodlike colloids from 3D microscopy data. J. Phys. Condens. Matter 27, 194109 (2015).

    CAS  PubMed  ADS  Google Scholar 

  47. 47.

    Hsiao, L. C. et al. Metastable orientational order of colloidal discoids. Nat. Commun. 6, 8507 (2015).

    CAS  PubMed  PubMed Central  ADS  Google Scholar 

  48. 48.

    Zheng, Z., Wang, F. & Han, Y. Glass transitions in quasi-two-dimensional suspensions of colloidal ellipsoids. Phys. Rev. Lett. 107, 065702 (2011).

    PubMed  ADS  Google Scholar 

  49. 49.

    Mundoor, H., Senyuk, B. & Smalyukh, I. I. Triclinic nematic colloidal crystals from competing elastic and electrostatic interactions. Science 352, 69–73 (2016).

    CAS  PubMed  ADS  Google Scholar 

  50. 50.

    Lee, T., Trivedi, R. P. & Smalyukh, I. I. Multimodal nonlinear optical polarizing microscopy of long-range molecular order in liquid crystals. Opt. Lett. 35, 3447–3449 (2010).

    CAS  PubMed  ADS  Google Scholar 

Download references


We thank M. Almansouri, B. Fleury, J. van de Lagemaat, Q. Liu, L. Longa, S. Park, B. Senyuk and J.-S. Tai for discussions and technical assistance. This research was supported by the US Department of Energy, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering, under contract DE-SC0019293 with the University of Colorado at Boulder.

Author information




H.M. synthesized disks and performed experiments and J.-S.W. performed Landau–de Gennes modelling of disk-induced distortions of the molecular order, under the supervision of I.I.S. H.H.W. designed the theoretical model and performed the numerical calculations of colloidal ordering of disks, with input from I.I.S. H.M. and I.I.S. analysed data. I.I.S. conceived, designed and directed the research, provided funding and wrote the manuscript, with input from all authors.

Corresponding author

Correspondence to Ivan I. Smalyukh.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature thanks Fulvio Bisi 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

Extended Data Fig. 1 Characterization of the building blocks of nematic colloids.

a, Three-photon excitation fluorescence spectra of 5CB obtained for polarizations (P) of excitation light parallel and perpendicular to nm. The inset shows the chemical structure and length of a 5CB molecule. b, Scanning electron micrographs of silica-coated disks. c, Transmission electron micrograph image of a single disk on a copper grid. d, A zoomed-in view of the region in the red square in c, showing the silica layer. e, Schematic of a disk with the NaYF4 core, silica layer and polyethylene glycol functionalization. f, Photon-upconverting luminescence spectrum of disks used in confocal imaging. gi, Simplified schematics of an effective building block of a molecular–colloidal LC formed by a disk (blue) in a fluid host of molecular rods (grey), illustrating its orthorhombic (g), monoclinic (h) and uniaxial (i) symmetries at different orientations. j, Orientational probability versus θ at different T values marked next to the distributions. The peaks of the curves correspond to θe, as depicted with a black dashed line at θe = 32° ± 1° for T = 33.0 ± 0.1 °C. The error in measuring θ is ±1°. k, Schematics of the main notations. l, Asymmetry of the colloidal surface anchoring potential Ue(θ) ≈ ε2(θ − θe)2 + ε3(θ − θe)3 + ε4(θ − θe)4, plotted for ε2 = 88, ε3 = −550 and ε4 = 1,608. m, Time-lapse bright-field micrographs of a disk in 5CB, showing its field-induced motion when subjected to an electric field EDC (black arrow) generated by a d.c. voltage of 5 V applied to in-plane electrodes 1 mm apart. The elapsed time and nm orientations are marked on the micrographs. n, Displacement versus time for disks with different charges in 5CB when subjected to 5 V between electrodes 1 mm apart. o, Histograms of the disk Brownian displacement probed along and perpendicular to nm using video microscopy, with the inset showing a bright-field micrograph of a disk in 5CB. Solid curves are Gaussian fits to experimental data. Errors in measuring the displacements in n, o are ±10 nm. a.u., arbitrary units.

Source data

Extended Data Fig. 2 Comparison and visualization of perturbed nematic order at different disk orientations.

a, Comparison of colloidal particles with tangential surface boundary conditions shaped as thin disks with different orientations (top) and a microring (bottom) in 5CB imaged by bright-field microscopy under similar conditions, with the latter inducing four boojum surface point defects visible as dark, light-scattering spots. No light-scattering singular defects are detected around the disks. bg, Numerical modelling of nm distortions around a disk at θe = 90°. bd, Visualization of energy-minimizing nm distortions around the disk, where the green isosurface-enclosed region indicates deviation of nm by >2° from its uniform far field and the black isosurface-enclosed region shows where Sm is reduced by >5% below its bulk equilibrium value. Colours of streamlines depict the opposite directions of local nm tilting relative to the midplane of the disk. eg, nm distortions near the same disk at 0 nm (e), 10 nm (f) and 50 nm (g) away from its midplane. hm, Numerical modelling of a disk-induced nm corona at θe = 75°. hj, Visualization of energy-minimizing nm distortions around the disk, where the green isosurface-enclosed region indicates deviation of nm by >2° from its uniform far-field background. The grey pane in i, j that intersects the green isosurfaces helps to show the monoclinic symmetry of nm distortions relative to the midplane of the disk. km, Isosurfaces enclosing small regions with Sm reduced by >5% below its bulk value (black surfaces) and colour-coded local nm distortions around the disk. Right-side insets in k depict details of nematic order perturbations at the disk edges. Field lines in l are shown 250 nm away from the disk centre and the single mirror-symmetry plane. nr, Visualizations similar to hm, but for θe = 30°. Note the opposite nm tilt at disk edges compared to that in hm, as highlighted by orange/blue colouring.

Extended Data Fig. 3 Uniaxial nematic dispersions and magnetic field alignment of disks.

ac, Uniaxial nematic phase characterized by upconverting luminescence confocal images (a, c) and numerical visualization (b) for in-plane (a) and out-of-plane (c) nm, showing the random orientations of the disk normals ωc in the plane perpendicular to nm. The orientation of the cross-sectional confocal image shown in a is depicted in b with a grey plane. d, Confocal image of a sample within a uniaxial–orthorhombic co-existence region with local correlations of ωc orientations. The disk number density is ρ = 0.16 μm−3 (a), ρ = 0.11 μm−3 (c) and ρ = 0.21 μm−3 (d). eh, Magnetic alignment of initially randomly oriented disk normals ωc (e, f) to point perpendicular to both nm and B at magnetic field amplitude B = 100 mT (g, h). Visualizations in f, h correspond to the images in e, g, respectively; ρ = 0.09 μm−3 and Z*e ≈ +80e in e, g. i, j, Confocal image (i) and numerical visualization (j) for a uniaxial nematic with disks oriented obliquely to nm at T = 31.5 ± 0.1 °C and ρ = 0.08 μm−3. k, l, Confocal image (k) and numerical visualization (l) of a uniaxial nematic at ρ = 0.19 μm−3 and T = 34.0 ± 0.1 °C.

Extended Data Fig. 4 Temperature-dependent orientations and magnetic control of disks.

ah, Upconversion-based luminescence confocal images (ac, eg) and corresponding numerical visualizations (d, h) of disk dispersions at ρ = 0.19 μm−3 with in-plane nm and ωc at T = 27.0 ± 0.1 °C (a), T = 28.0 ± 0.1 °C (b), T = 30.0 ± 0.1 °C (c), T = 31.0 ± 0.1 °C (e), T = 32.0 ± 0.1 °C (f) and T = 33.0 ± 0.1 °C (g), showing the temperature-dependent disk orientations in the presence of a magnetic field. With increasing T at a field of B = 100 mT applied along a normal to the sample, a domain with a magnetically induced unidirectional alignment of ωc (ad) transforms into a polydomain state at higher temperatures owing to the random ±θe tilting directionality of ωc relative to nm. Image planes are depicted in grey in fh. The orientations of nm, magnetic field and scale bars shown between a, b and e, f are the same for all confocal images presented. The disk surface charge is Z*e ≈ +80e.

Extended Data Fig. 5 Coarse-grained visualization of LC phases and background alignment.

ai, Simplified schematics of different hybrid phases of 5CB molecular rods (grey) and colloidal disks (blue) for the orthorhombic nematic (a), monoclinic nematic (b), monoclinic smectic (c), columnar (d), orthorhombic columnar nematic (e), uniaxial nematic (f), uniaxial nematic with colloidal columns (g), isotropic (h) and colloidal nematic (i) states. An elementary columnar cell in d is depicted in red, with details shown in Fig. 5h. j, Phase diagram upon variation of ρ, T and Z*e, with ‘Co-ex’ referring to co-existence regions. Disks aggregate at Z*e < +10e and form Wigner-type crystals at Z*e > +100e (ref. 10). T is measured with an error of ±0.1°; relative errors for ρ and Z*e are ±5% and ±1%, respectively. k, 3PEFPM image of a 5CB dispersion of disks at ρ = 0.32 μm−3 and T = 27.0 ± 0.1 °C for linear polarization of the excitation light parallel to nm. l, 3PEFPM intensity, averaged over the field of view, versus the angle between the linear polarization of the excitation light and nm (measured with an error of less than ±1°) for 5CB with colloidal disks and pure 5CB under the same conditions. m, n, Visualization of nm between two disks at θe = 75° (m) and θe = 30° (n), with black isosurfaces enclosing spatial regions in which Sm is reduced by >5% relative to its equilibrium bulk value. Regions of distorted nm are highlighted by coloured streamlines, with blue and yellow depicting opposite tilt directions.

Source data

Extended Data Fig. 6 Colloidal director rotation and nematic ordering versus temperature.

ai, Numerical visualizations (a, e, i) and upconversion-based luminescence confocal images (bd, fh) of disks in the monoclinic and orthorhombic nematic phases with out-of-plane nm at disk number density ρ = 0.31 μm−3 and at T = 27.0 ± 0.1 °C (b), T = 28.0 ± 0.1 °C (c), T = 30.0 ± 0.1 °C (d), T = 31.0 ± 0.1 °C (f), T = 32.5 ± 0.1 °C (g) and T = 33.0 ± 0.1 °C (h). Note the pretransitional smectic correlations in the vicinity of the smectic phase temperature range in f, g. nc tilts out of the sample plane owing to the temperature-dependent variation of θne, consistent with the elliptical cross-section of disks, revealing oblique orientations of ωc and nc relative to nm at T = 32.5 ± 0.1 °C (g) and T = 33.0 ± 0.1 °C (h). Grey planes in a, e, i depict the orientation of confocal images relative to the directors. Scale bars and nm orientation shown between panels b, c and f, g are the same for all images. j, fc(θφ) versus azimuthal angle φ of ωc at ρ = 0.34 μm−3 and different T values. k, fc(θφ) versus φ at T = 27.0 ± 0.1 °C for different ρ values. l, Δ (red squares) and S (black circles) versus ρ at T = 27.0 ± 0.1 °C. The disk surface charge is Z*e ≈ +80e.

Source data

Extended Data Fig. 7 Characterization of colloidal interactions within different phases.

ac, Distribution of the angle ϑr between nm and the centre-to-centre pair separation vector, defined in the right-side inset of a, in the monoclinic nematic at ρ = 0.31 μm−3 and T = 30.7 ± 0.1 °C (a), orthorhombic nematic at ρ = 0.34 μm−3 and T = 27.0 ± 0.1 °C (b) and uniaxial nematic at ρ = 0.1 μm−3 and T = 27 ± 0.1 °C (c). The error in measuring ϑr is ±2°. df, Scatter plots of the nearest-neighbour disk positions in the monoclinic (d), orthorhombic (e) and uniaxial (f) nematic phases at ρ and T values corresponding to ac. Anisotropic distributions in d and e correlate with the ordering of disks. The disk surface charge is Z*e ≈ +80e.

Source data

Extended Data Fig. 8 Monoclinic and magnetic-field-induced triclinic smectic order.

ad, Upconversion-based luminescence confocal microscopy images (a, b) and numerical visualizations (c, d) of disks in the monoclinic smectic phase at ρ ≈ 0.3 μm−3, T = 32.3 ± 0.1 °C and θne ≈ 45° (a, c) and for the magnetically induced triclinic smectic state under otherwise similar conditions (b, d). The smectic layer normal (yellow double arrow), nc and nm are initially within the image plane (a, c), but then nc is rotated out of the image plane by a magnetic field of B = 100 mT (b), thus switching the smectic order from monoclinic (a, c) to triclinic (b, d). This switching is driven by the tendency of ωc and nc to align perpendicular to B, where the disk normals also tend to reside on a cone of anchoring-defined easy-angle ωc orientations. Right-side insets of a, b show the layer normal, nm, nc and n and the applied magnetic field B for the monoclinic and the triclinic smectic states, respectively. The layer normal, nc and n are all mutually orthogonal and the layer normal, nc and nm are all co-planar in a but not in b, which lacks nontrivial symmetry operations and represents triclinic symmetry.

Extended Data Fig. 9 Nematic and columnar phases and their thermal melting.

ac, Upconversion-based luminescence confocal micrographs (a, c) and numerical visualization (b) showing a uniaxial nematic with differently oriented columns at ρ = 0.08 μm−3 and T = 27.0 ± 0.1 °C (a, b) and the colloidal structure of the biphasic nematic–columnar co-existence region at ρ = 0.17 μm−3 and T = 27.0 ± 0.1 °C (c). dg, Monoclinic columnar phase with a two-dimensional oblique lattice characterized experimentally (f, g) and numerically (df). The inset in f shows the ordering of columns in a sample with out-of-plane nc and in-plane nm, revealing the oblique lattice of colloidal columns. Probing the distribution of the oblique lattice parameters defined in e across multiple sample regions yields angles of χ = 110° ± 5°; d1 = 1.67 ± 0.15 μm, d2 = 1.75 ± 0.15 μm. The angle between the diagonal of the oblique primitive cell and nm varies within 5°–7°. The measuring error is ±1° for angles and ±50 nm for d1 and d2. Coloured nm director streamlines in e depict nm distortions right above the disks. Orange/blue isosurfaces in f enclose regions in which local nm orientation departs from its uniform far-field background by >2°, with colours corresponding to the tilt directionality of the streamlines in e. g, Confocal cross-sections of the columnar phase showing disks at sample depths shifted by 2 μm sequentially, corresponding to the grey planes labelled in d; ρ = 0.32 μm−3 and T = 27.0 ± 0.1 °C. h, i, Confocal photon-upconverting images of disks at ρ = 0.19 μm−3 that show melting of the colloidal columns upon increasing the temperature from T = 27.0 ± 0.1 °C (h) to T = 33.0 ± 0.1 °C (i). j, k, Columns re-emerge upon decreasing T, as revealed by confocal images at T = 30.0 ± 0.1 °C (j) and T = 27.0 ± 0.1 °C (k). The disk charge is +20e.

Extended Data Fig. 10 Stability analysis of phases with partial positional order.

a, Schematic showing the principal reference frames, vectors and angles. b, Phase diagrams with coordinate axes θne versus ρ (left column) and T versus ρ (right column) of the molecular–colloidal LC for higher (0.002; top row) and lower (0.001; bottom row) values of the elastic energy Ξel (renormalized in units of thermal energy) that quantifies the relative strength of elastic attractions with respect to steric repulsion (Supplementary Information). Red-shaded regions represent the monoclinic smectic phase, consistent with the experimental diagram (Fig. 2a). For lower Ξel (weak elastic interactions), the phases with partial positional order appear at higher disk concentrations owing to stronger electrostatic repulsions between the disks. c, Theoretical estimation of the smectic layer spacing σsm versus T in the monoclinic smectic phase.

Supplementary information

Supplementary Information

This file contains the following sections: Supplementary Methods I: Details of theoretical modelling; Supplementary Methods II: Details of computer simulations; Supplementary Methods III: Parameter specifications for the experimental phase diagram; Supplementary References; and Appendix: Numerical codes.

Video 1: Comparison of confocal and brightfield imaging of a uniaxial LC

Confocal photon-upconverting (top left) and brightfield (top right) imaging videos show edge-on discs diffusing in a rectangular 0.2 × 2 mm capillary with out-of-plane nm. Confocal video tracks colloids through a series of depth-resolved scans, each revealing discs’ positions and orientations, but brightfield images show discs at different depths, impeding such characterisation. Elapsed times are marked; Z*≈ +80e. On the bottom-left diagram, a circle marks ρ = 0.1 µm−3 and T = 27.0 ± 0.1 °C. Bottom-right visualisation shows discs and the confocal scanning plane; green isosurfaces enclose regions where nm departs from its far-field orientation by >3°.

Video 2: Discs within an orthorhombic nematic

Time-sequential confocal scanning video shows orientations and diffusion of discs at T = 27.0 ± 0.1 °C, ρ = 0.34 µm−3 (first part of the video, out-of-plane nm) and ρ = 0.31µm−3 (second part, in-plane nm). Discs self-align with ncnm. Inset in the last part of the video is a zoomed-in view within the same sample. Elapsed time within each movie fragment (from different samples) is marked on the frames; ρ and T are shown by a filled circle in the bottom-right diagram. Top-right visualisation shows disc ordering and the scanning plane; green isosurfaces enclose regions where nm departs from its far-field by >3°. Z*e ≈ +80e.

Video 3: Discs in a monoclinic nematic LC

Confocal video shows diffusion and ordering of discs at different ρ and T (moving circle in the top-right diagram). Images have monodomain nm and polydomain (beginning of the movie) or monodomain nc oriented obliquely to nm, as shown by double arrows. Note the ellipsoidal appearance of discs in the frames where discs orient obliquely to the image plane (last part of the video). Elapsed time within each movie fragment (from different samples) is marked on the frames. Bottom-right visualisation of monoclinic order also shows cross-sectional planes corresponding to confocal images. Z*≈ +80e.

Video 4: Numerical visualisation of perturbed nematic order of the LC host around a disc

Video shows energy-minimising nm-distortions around a disc obtained from Landau-de Gennes numerical modelling. Green isosurface-enclosed regions (left) indicate deviations of nm by >2° from its uniform far-field and the black isosurface-enclosed regions (right) show where Sm is reduced by >5% below its bulk equilibrium value. Orientation of nm and the associated field lines (pink) stay vertical with changing θe; orientation of \({{\boldsymbol{\omega }}}_{{\rm{c}}}\) is shown with a white double arrow.

Video 5: Visualisation of orthorhombic, monoclinic and uniaxial colloidal nematic states

Video shows ordering of discs in nematic phases at different θne. Green isosurface-enclosed regions (left) indicate deviations of nm by >2° from its far-field; black isosurfaces enclose regions (right) where Sm reduces by >5% below its bulk equilibrium value. Pink lines/arrows indicate the far-field nm; blue and green double arrows depict nc and n, respectively. θe = 90° corresponds to the orthorhombic phase, 0° < θe < 90° to the monoclinic state and θe = 0° to the uniaxial phase. For simplicity, the same disc positions are used at different θne; no smectic correlations are considered at intermediate θne.

Video 6: Smectic order and magnetic control

First part of video shows smectic order at ρ and T indicated by a filled circle in the top-right diagram, with disc orientations/positions and directors depicted in the bottom right. The second part compares the monoclinic smectic (top-left) to the smectic with induced triclinic symmetry when nc tilts out of the sample plane in response to |B| = 100 mT (top-right) while discs satisfy the anchoring conditions and tend orienting with normals orthogonal to B; visualisations beneath show orientations of the layer normal (yellow double arrow), B and directors. Elapsed time within each movie fragment is marked; Z*≈ +80e.

Video 7: Uniaxial colloidal nematic order of discs within 5CB in nematic and isotropic states

The first part of the video shows the colloidal orientational order of discs dispersed within 5CB in uniaxial nematic state whereas the second half corresponds to 5CB being in the isotropic phase. The corresponding parameters ρ and T of the phase diagram and schematics of ordering are depicted in the top-right and bottom-right, respectively. The elapsed time is marked on the frames within each movie fragment. The disc surface charge is Z*≈ +80e.

Video 8: Out-of-equilibrium thermal reconfiguration

Video shows control of colloidal order upon heating from 27 °C to 32.2 °C and then cooling within 990s. Synchronously with video, black filled circle shows ρ and T (top-left diagram); bottom-left visualisation shows disc ordering. Reconfigurations emerge from tendency to form temperature-dependent equilibrium structures, including orthorhombic/monoclinic nematics and monoclinic smectic. Discs align with nc orthogonal to frames at T = 27.0 °C, then gradually tilt to have in-plane nc rotating to θne ≈ 45° relative to nm, exhibiting first nematic and then smectic monoclinic order; smectic layers dissociate and nc rotates back upon decreasing T. Elapsed time is marked on frames; ρ = 0.30 µm−3; Z*≈ +80e.

Video 9: Columnar orthorhombic nematic order

Video shows colloidal columns formed by discs, representing a columnar D2h nematic at T = 27.0 °C and ρ = 0.23 µm−3 (marked in the bottom-right diagram). Numerical visualisation (top right) shows colloidal structures, with coloured isosurfaces depicting regions where nm-departures away from its far-field background are >3° in opposite directions (different colours). Individual discs self-organise into columns as a result of competing electrostatic and anisotropic elastic interactions, which spontaneously align and define ncnm. Video was obtained by time-sequential confocal scanning in a 100-μm-thick glass cell ~10-μm away from one of the inner surfaces. Elapsed time is marked on the frames; Z*≈ +20e.

Video 10: Orthorhombic nematic order emerging from melting colloidal columns

Video of confocal (left) and brightfield (right) sequential images showing flow and Brownian motion of orientationally-ordered edge-on discs within the orthorhombic nematic phase formed by melting of colloidal columns at T = 34.0 °C and ρ = 0.29µm-3 (marked in the bottom-left diagram) in a rectangular glass capillary with a 200 μm by 2 mm cross-section. Disc order is shown in the bottom-right visualisation, with green isosurfaces enclosing regions where nm departs from its far-field orientation by >3°. Elapsed time is marked on confocal and brightfield video frames; Z*≈ +20e.

Source data

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Mundoor, H., Wu, JS., Wensink, H.H. et al. Thermally reconfigurable monoclinic nematic colloidal fluids. Nature 590, 268–274 (2021).

Download citation


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.


Quick links