Shape and structural relaxation of colloidal tactoids

Oct 19, 2023

Nature Communications volume 13, Article number: 2778 (2022) Cite this article

3042 Accesses

2 Citations

8 Altmetric

Metrics details

Facile geometric-structural response of liquid crystalline colloids to external fields enables many technological advances. However, the relaxation mechanisms for liquid crystalline colloids under mobile boundaries remain still unexplored. Here, by combining experiments, numerical simulations and theory, we describe the shape and structural relaxation of colloidal liquid crystalline micro-droplets, called tactoids, where amyloid fibrils and cellulose nanocrystals are used as model systems. We show that tactoids shape relaxation bears a universal single exponential decay signature and derive an analytic expression to predict this out of equilibrium process, which is governed by liquid crystalline anisotropic and isotropic contributions. The tactoids structural relaxation shows fundamentally different paths, with first- and second-order exponential decays, depending on the existence of splay/bend/twist orientation structures in the ground state. Our findings offer a comprehensive understanding on dynamic confinement effects in liquid crystalline colloidal systems and may set unexplored directions in the development of novel responsive materials.

Colloidal liquid crystals are a class of soft matter formed when shape-anisotropic nanoparticles are dispersed in a an isotropic fluid1. When confined to a finite volume, rod-like colloidal particles self-organize into various structures that are set by a delicate balance between anisotropic viscoelastic and surface properties2,3,4. The subtle balance between these contributions results in facile response to external fields1,5 such as mechanical, flow, electric, and magnetic, giving rise to many opportunities and technological applications. Examples include displays, spatial light modulator and tunable filters in medical devices and optics, liquid crystal biosensors for rapid diagnostics, and new functional material such as artificial muscles exploiting liquid crystalline anisotropic physical properties6,7,8. However, facile responsiveness to external fields (and disturbance) makes the colloidal liquid crystals very fragile to be studied experimentally under dynamical conditions5,9. In particular, relaxation of liquid crystalline droplets under mobile confinement is still poorly understood despite its central importance in a variety of phenomena in condensed matter physics. This includes particle packing10, self-assembly11, and relaxation of colloidal liquids12 with implications in the field of active nematic, e.g., living liquid crystals13, where the understanding of the hydrodynamics of the liquid crystals is critical14.

Liquid crystalline droplets, known as tactoids, are a particularly significant example of colloidal liquid crystals, since they consist of micro-confined liquid crystalline colloids with a self-selected shape/structure stemming out from the thermodynamically-driven phenomena from which they emerge, i.e., spontaneous liquid–liquid crystalline phase separation15,16,17,18,19,20,21,22,23. In stark contrast with spherical liquid crystalline emulsions, achieved commonly by emulsifying liquid crystals in another immiscible liquid (like water in oil)24,25, tactoids hold spindle-like, prolate, or oblate shapes with different nematic-cholesteric internal structures2,3,4,15,16,17,18,19,20,21,22,26, as a consequence of the coupling between the vanishingly small interfacial tension, the surface anchoring at the interface, the chirality of colloids and the anisotropic elastic properties3,4. These features make tactoids a very unique system with peculiar viscoelastic27,28,29 and boundary3,4 properties, thus adding theoretical challenges to the experimental ones when describing these complex colloidal systems under dynamical conditions. For instance, recent experiments suggest that the boundary has a significant impact on the local structure of colloids30,31,32,33 and on the equilibration pathways of structural relaxation of colloidal systems34; yet such understanding mainly comes from the examination of colloidal systems with static boundary conditions30,31,32,33,34,35. Moreover, one of the challenges of the current study is to disentangle the rate of self-assembly of the liquid crystalline tactoids from shape relaxation, providing insights on the kinetics of self-assembled complex colloidal systems.

Here we report the shape and structural relaxation dynamics of colloidal liquid crystalline tactoids. We use β-lactoglobulin amyloid fibrils and cellulose nanocrystals as model rod-like colloidal liquid crystalline systems. We disentangle shape and structural relaxation and show—by integrated experimental and numerical measurements— that the shape relaxation of the tactoids follows a single exponential decay that depends on the material isotropic and anisotropic properties and the size of the droplets. We develop a theoretical model to predict the shape of the tactoids out of the equilibrium state, by considering the Hamiltonian of the tactoids in presence of an external flow field. We also show that the structural relaxation of the tactoids follows different paths depending on the colloidal mesogens configurations at the ground state; homogenous and bipolar tactoids relax through a first-order exponential decay whereas cholesteric ones follow a second-order exponential decay. We use direct experimental measurements of the order parameter, supported by direct numerical simulation (DNS) of a complete structure-composition multiscale model, to discuss the nature of the structural relaxation of the liquid crystalline droplets and how it relates to the splay, bend and twist structures of tactoids at relaxed state. Our results offer original insights on the structural organization of colloidal suspensions out of equilibrium, under dynamic confined boundaries and evolving shapes.

The classical approach to study relaxation of droplets involves using the four-roll mill geometry developed originally by Taylor36. This approach, however, is not applicable to tactoids due to their labile nature, making prohibitive isolating a single tactoid in such a geometry. Thus, in our experiments we take advantage of a microfluidic system with contraction-abrupt expansions design37,38, allowing to elongate tactoids with different volumes and then let the elongated tactoids relax to the equilibrium state in the abrupt expansion zone (the details on the microfluidic channel are provided in Methods and Supplementary Note 1). To be able to form the tactoids in a microfluidic chip, we prepared the suspension of the liquid crystalline with a concentration that is set within the isotropic–nematic coexistence region. After injection of the suspension into the microfluidic system, tactoids with various volumes are formed inside the channel following nucleation and growth path. Series of experiments with β-lactoglobulin amyloid fibril and cellulose nanocrystals liquid crystalline droplets were performed and analyzed under crossed polarizers and LC (liquid crystal)-PolScope device, allowing capturing not only the shape of the tactoids but also their internal structure (Fig. 1), see Supplementary Movies 1–3. Regardless of the tactoids volume, all the tactoids at the initial state are extended and the director field inside the tactoids is aligned parallel to the long axis of the tactoid, known as homogenous configuration2,4. During the relaxation, an initially extended droplet with volume V (≈r2R with R the major and r the minor axes of tactoids) ~102 µm3 undergoes the shape relaxation while its structure remains unchanged as homogenous configuration (Fig. 1a). As explained in detail later, while its configuration remains unchanged as homogenous, the tactoid still undergoes structural relaxation. A tactoid with a larger volume, ~103 µm3, shows relaxation on both its shape and structure where the director field changes from homogenous structure to bipolar with a director field that smoothly follows the tactoid interface (Fig. 1b). For a tactoid with volume ~104 µm3 as shown in Fig. 1c, while its shape relaxes to a nearly—yet not perfectly—spherical shape, the structural relaxation takes place while changing the director field from homogenous to cholesteric configuration, that is easily distinguishable from its characteristic striped texture.

The sequence of time-lapse images of relaxation of initially extended amyloid fibrils liquid crystalline droplets with different volumes. In each panel, the first row shows the experimental results taken with LC (liquid crystal)-PolScope device. The colormaps corresponding to experimental results denote the orientation of the director filed in the x–z plane; the second row demonstrates the numerical simulation results with color bar capturing the director field orientation with respect to z-axis. The tactoids are at the homogenous configuration at the initial state and upon relaxation, they hold different configurations. a An initially extended tactoid with volume 644 µm3 undergoes shape relaxation while its configuration remains unchanged at homogenous configuration. b An elongated tactoid with volume 2751 µm3 relaxes both its shape and structure recovering a bipolar configuration upon relaxation. c A droplet with volume 16,414 µm3, having larger volume compared to (a, b), relaxes to a cholesteric structure with three bands. Since colors and director lines show the same information due to the axisymmetric nature of the homogenous (a) and bipolar (b) tactoids, the lines are not shown for better readability (see Supplementary Movies 4–6 for the version with lines). Note that the brightness of the experimental images is varied for better visualization.

We additionally captured both the shape and structural relaxation of the tactoids using DNS (see Supplementary Movies 4–6), and found good agreement with our experimental results as shown in Fig. 1. Details of the DNS are given in the Methods section.

To characterize the underpinning physics of the relaxation of the liquid crystalline droplets, we first analyze the shape of the tactoids during the relaxation. We measured the long axis of tactoids at different time t and quantified the relaxation behavior, as suggested previously for homogenous droplets39, by \({{{{{\mathcal{R}}}}}}=\tfrac{R\left(t\right)\,-\,{R}_{{{{{{\rm{equil}}}}}}.}}{{R}_{{{{{{\rm{init}}}}}}.}\,-\,{R}_{{{{{{\rm{equil}}}}}}.}}\), where Requil., Rinit., and R(t) are the half-length of the long axis of tactoid at equilibrium, at the initial time, and at a given time, respectively (Fig. 2a–c). The value of \({{{{{\mathcal{R}}}}}}\) is one at time zero and zero when the tactoid reaches its equilibrium shape. We observed that \({{{{{\mathcal{R}}}}}}\) for all classes of tactoids (homogenous, bipolar, and cholesteric configurations) follows essentially a single exponential decay (Fig. 2a–c). Thus, to capture the shape relaxation behavior of the tactoids, we use \({{{{{\mathcal{R}}}}}}\,=\,{{{{{\rm{exp}}}}}}(-t/{\tau }_{{{{{{\rm{s}}}}}}})\), where the τs is the characteristic shape relaxation time of the tactoids. A good agreement between the fitting, using exponential decay \({{{{{\mathcal{R}}}}}}\,=\,{{{{{\rm{exp }}}}}}(-t/{\tau }_{{{{{{\rm{s}}}}}}})\), and the experimental data of the shape relaxation of the tactoids with different relaxed configurations, allows us to obtain the τs from the shape decay curve (Fig. 2a–c). Hence, a good collapse of the shape relaxation data of the tactoids with different volumes and relaxed configurations onto a single master curve is observed when time is rescaled with τs (Fig. 2d). We additionally performed experiments on the shape relaxation of the tactoids resulting from breakup events of initially extended tactoids in homogeneous configuration, and observed a single exponential decay trend similar to the initially extended tactoids relaxation behavior, see Supplementary Fig. 2.

a–c Evaluation of \({{{{{\mathcal{R}}}}}}\) (defined as \(\tfrac{R\left(t\right)\,-\,{R}_{{{{{{\rm{equil}}}}}}.}}{{R}_{{{{{{\rm{init}}}}}}.}\,-\,{R}_{{{{{{\rm{equil}}}}}}.}}\) where Requil., Rinit., and R(t) are the half-length of the long axis of tactoid at equilibrium, at the initial time, and at a given time t, respectively) with respect to time for tactoids that relax to homogenous, Requiv. (V1/3, with V the volume) = 9.4 µm (a), bipolar, Requiv. = 19.3 µm (b) and cholesteric, Requiv. = 27.8 µm (c) configurations at equilibrium. Symbols and black lines denote the experimental and numerical simulation results, respectively; colored lines show the fitting (\({{{{{\mathcal{R}}}}}}={{{{{\rm{exp }}}}}}\left(-\frac{t}{{\tau }_{{{{{{\rm{s}}}}}}}}\right)\)) that is used to obtain the characteristic shape relaxation time, τs. d Evaluation of \({{{{{\mathcal{R}}}}}}\) with respect to scaled time t⁄τs resulting in a universal curve, \({{{{{\mathcal{R}}}}}}={{{{{\rm{exp }}}}}}\left(-\frac{t}{{\tau }_{{{{{{\rm{s}}}}}}}}\right)\), for shape relaxation of the different classes of tactoids with various volumes and initial elongation values. e Circle, triangle, and square symbols denote homogenous, bipolar, and cholesteric tactoids, respectively. The error bars represent standard deviation. The developed theory, solid line, predicts the τs for different classes of BLG and SCNC liquid crystalline tactoids, confirming the generality of our approach to predict the bio-colloidal liquid crystalline tactoids relaxation behavior. Here, ω is the anchoring strength, ckB T is the thermal energy per unit volume of dispersion, K and K2 are the Frank elastic constants, ξ is the coherence length, M∅ is the mass mobility, MQ is the rotational mobility, γ is the interfacial tension, βμI is the effective viscosity, and b is a single constant pre-factor.

To calculate τs, we combine theories on simple droplet relaxation, dimensional analysis, and DNS that shows an excellent agreement with our experiments, see Figs. 1, 2a–c. As elaborated in Supplementary Note 2, we find:

where τs is expressed as the sum of two contributions: the first term is the liquid crystalline anisotropic contribution (τa) induced by the presence of colloidal mesogens accounting for orientational order, gradient elasticity, anisotropic viscoelasticity, rotational dissipation, and concentration gradients while the second term is the characteristic shape relaxation time of elongated isotropic tactoids (τi). We express τi, following the well-established relation of the characteristic shape relaxation time of elongated isotropic droplets37,40,41, as \({\tau }_{{{{{{\rm{i}}}}}}}\,=\,\frac{\beta {\mu }_{{{{{{\rm{I}}}}}}}{R}_{{{{{{\rm{equiv}}}}}}{{{{{\boldsymbol{.}}}}}}}}{\gamma }\), where γ is the interfacial tension, Requiv. = ((r2R)1/3) is the equivalent radius of tactoids, μI the viscosity of the medium phase -taken to be equal to viscosity of the isotropic phase- and \(\beta \,=\,\frac{\left(2\hat{\eta }\,+\,3\right)\left(19\hat{\eta }\,+\,16\right)}{40\left(\hat{\eta }\,+\,1\right)}\), where \(\hat{\eta }\,=\,\frac{{\mu }_{{{{{{\rm{N}}}}}}}}{{\mu }_{{{{{{\rm{I}}}}}}}}\) is the ratio of viscosities of the nematic phase μN and that of the isotropic medium, μI. We obtain τa, as noted in Supplementary Note 2, based on dimensional analysis and parametric studies through our validated DNS on material properties involved in the liquid crystalline self-assembly and tactoid size. In Eq. 1, ω is the anchoring strength, which takes into account the coefficient of the concentration-orientation gradient and its coupling in the numerical simulations42 (Supplementary Note 2), the term ckBT is thermal energy per unit volume of dispersion with c, kB, and T, the number density, Boltzman constant, and temperature, respectively. The term K is the Frank elastic constant for splay and bending (assumed to be equal) and K2 is the Frank twist elastic constant. It should be noted that, from both theory43,44 and experimental measurements45,46,47 on different systems of rigid rod-shaped liquid crystals, including filamentous colloids48 analogous to those studies here, the ratio of \(\frac{K}{{K}_{2}}\) in Eq. 1 is always greater than ½, and thus Eq. 1 only contains real arguments. The term ξ represents the coherence length that is an indicator of length over which long-range ordering takes place. The term \({M}_{{{\varnothing }}}\,\propto\, \frac{{{{{{\rm{ln}}}}}}\left(L/D\right)}{{{{{{{\rm{c}}}}}}\mu }_{{{{{{\rm{N}}}}}}}L}\) is the mass mobility and \({M}_{{{{{{\rm{Q}}}}}}}\,\propto\, \frac{2{{{{{\rm{ln}}}}}}\left(2L/D\right)\,-\,1}{{{{{{{\rm{c}}}}}}\mu }_{{{{{{\rm{N}}}}}}}{L}^{3}}\) is the rotational mobility49 (∝ stands for proportionality); here L and D, respectively, the length and diameter of the rod-like mesogen, assumed to be equal to the weighted mean length of the fibrils Lf,w and effective diameter Deff. proposed by Onsager50, respectively. We list all properties of the liquid crystals used in this study in Table 1 along with the details on calculations and measurements in Supplementary Notes 3, 4, supported by refs. 51,52,53,54,55,56,57,58,59. The relation for τa is formulated through a hybrid approach based on DNS results and dimensional analysis and can be turned to an equation by use of a single constant pre-factor (b = 54.0). The term b is thus an aggregated value reflecting pre-factors present in the proportionality terms such as mobilities. We compare the experimental results for τs with our prediction with a single fitting parameter b and find an excellent agreement (see Fig. 2e).

To test the generality of the present approach, in addition to BLG I in Table 1, we also perform experiments with one other system of amyloid fibrils (BLG II) having different length distribution compared to BLG I as a result different material properties60 and with sulfated cellulose nanocrystals (SCNC); see Table 1. The tactoids of both systems also illustrate single exponential decay during the shape relaxation like BLG I in Fig. 2a–c. The results for characteristic shape relaxation time of BLG II and SCNC show very good agreement with Eq. 1 prediction (Fig. 2e), suggesting that Eq. 1 is general enough to describe the relaxation behavior of most bio-colloidal liquid crystalline tactoids.

Having established a general picture on the dynamic shape relaxation of the tactoids, we go into modeling the deformation of the tactoids under external stresses. While this has been well documented for simple fluids40,61, for liquid crystalline tactoids the physics become complex due to the energy terms associated with the internal structure of the tactoids, their anisotropic viscoelasticity and the confining boundary features. We look at the deformation of the droplet under uniaxial flow field with extension rate given \(\dot{\varepsilon }_{xx}\,{{\mbox{=}}}\,\frac{\partial {u}_{x}}{\partial x}\), where the ux is the flow speed and x is the direction of the motion of the flow, although the approach is general and can be used to model tactoids deformation under any external stresses. We consider the energy gained by the tactoids under the external stresses imposed by the extensional flow field and incorporate it to free-energy landscape of the tactoids that is well described by a scaling form of Frank–Oseen elasticity theory2,4. In particular, the rate of energy \(\frac{{dE}}{{dt}}\) gained by the tactoids under any external normal stresses σ can be expressed as \(\frac{{dE}}{{dt}}\,=\,\int {{{{{\boldsymbol{\sigma }}}}}}.{{{{{{\bf{u}}}}}}}_{{{{{{\rm{i}}}}}}}\,{dS}\), where ui is the displacement of the interface of the tactoids. We consider the case where the tactoid is stretched due to the stresses applied by uniaxial extensional flow and in Supplementary Note 5 we calculate \(\frac{{dE}}{{dt}}\) to be

The free-energy landscape of the tactoids includes two energetic terms, the bulk elastic and surface free energies, associated with the tactoid at equilibrium. The total free energy of the tactoid FE is described in scaling form as2,4

where the first term accounts for the surface free energy of the tactoids that is due to the interfacial tension and the anchoring strength. The last two terms are the bulk elastic free energy of the tactoids where the second term accounts for splay and bending free energy and the third term is the twist elastic free energy. The term θ (=n ∙ ∇ × n with n the nematic director) is the twist term in the Frank–Oseen elasticity theory and the term q∞ (=2π/P∞ with P∞ the natural pitch of the system) is the chiral wave number. We measured and summarized the properties of the suspensions used in this study in Table 1 (see also Supplementary Note 4).

By energy conservation, the rate of the energy gained by the tactoids due to the normal stresses from the flow field must be equal to the rate of the energy changes in the free energy of the tactoids associated with their elastic/interfacial energy. Note that setting the two energies equal stands valid here since the process happens significantly faster than the rate at which the heat can flow out, thus separating the time scale for the transfer of the energy associated with structural changes from the time scale for heat dissipation. Additionally, as all three classes of the homogenous, bipolar, and cholesteric tactoids hold a homogenous configuration under extreme deformation, as can be seen in Fig. 1 and our recent study9, we ignore the second term in Eq. 3, implying, as usual for homogeneous tactoids, that the bulk elastic energy due to the splay and bending is zero. Additionally, the third term is eliminated as, in the homogenous configuration and at constant tactoid volume, it does not change under deformation, so the rate of the energy gained by this term becomes zero. All in all, by setting Eq. 2 equal to the time derivative of Eq. 3 yields:

giving the steady-state elongated shape of the tactoids under a given extensional flow field in terms of r as a function of \(\dot{\varepsilon }\) and V. There is no analytical solution for Eq. 4, but the numerical solution along with our experimental data are presented in Fig. 3. Here, to best match experimental observations, the second term is rescaled by a pre-factor of 0.14, which is fully justified by the use of a scaling form of the Frank–Oseen energy landscape. Our results suggest that the short axis of the tactoids r decreases as the extension rate increases, where r lines, corresponding to tactoids with different volumes, converge to a single universal curve at large values of extension rate (Fig. 3a). The most remarkable consequence of our analysis is that for high extension rates, r becomes independent of the volume, that is, the cross-section of the tactoid is simply ruled by extension rate, and that at identical extension rates, tactoids of different volumes V only differ by their long radius R, which is directly proportional to V. In the regime of low extension rate, the short axis of the tactoids becomes volume-dependent and increases logarithmically with the increase in the volume of the tactoids (Fig. 3b), which is nicely supported by the collected experimental dataset. Equation 4 can be used to predict the Rinit.(=V/r2) as the maximum tactoids deformation that can be reached under a given extensional flow rate. Note that in Fig. 3a, b the experiments of the large volume of the tactoids (e.g., V = 30,000 μm3) at the high extension rate is limited by the experimental setup. This is so as the length of the tactoids in the stretched forms becomes extremely high compared to the size of taken images, preventing to fully capture tactoids with large volume at high extension rate.

a The theory (lines) and the experimental data (symbols) predict that short axis of the tactoids r decreases as the extension rate \(\dot{\varepsilon }\) increases where r lines, corresponding to tactoids with different volumes V, converge to a single curve at large values of extension rate. b At a given extension rate, the short axis of the tactoids increases logarithmically with an increase in the volume of the tactoids.

The above treatment describes comprehensively the evolution of the confining boundaries of the tactoids, but provides no information on the evolution of their internal structure. We thus turn our attention to the study of the orientational order parameter of the director field during relaxation of the tactoids. In our experiments, we captured the relaxation of the tactoids by LC-PolScope allowing us to access the retardance images of the tactoids during relaxation and analyze the order parameter S, where S = \({{{{{\mathcalligra{r}}}}}}\)/d∆n0 with \({{{{{\mathcalligra{r}}}}}}\) the optical retardance value, d the thickness of the sample and ∆n0 the birefringence corresponding to a perfectly aligned nematic phase, i.e., when the order parameter is 162,63. We measured the retardance value \({{{{{\mathcalligra{r}}}}}}\) of every pixel within the tactoids, and accordingly, the d is calculated for every pixel assuming a spindle shape for the tactoids (see Supplementary Note 6). The exact value of the ∆n0 is unique for a given liquid crystalline system and is often challenging to obtain experimentally, thus here we present our calculation independent form ∆n0. We define \({{{{{\mathcal{S}}}}}}=\frac{S\left(t\right)\,-\,{S}_{{{{{{\rm{equil}}}}}}.}}{{S}_{{{{{{\rm{init}}}}}}.}\,-\,{S}_{{{{{{\rm{equil}}}}}}.}}\) to capture the structural relaxation of the tactoids, similar to the one used for the shape relaxation of the tactoid, and most importantly, fully independently of ∆n0. The experimental results of \({{{{{\mathcal{S}}}}}}\) obtained for tactoids with different relaxed configurations against scaled time are shown in Fig. 4a–c. The structural relaxation of the tactoids with small volumes that hold homogenous and bipolar at equilibrium follows a first-order exponential decay. However, for the larger extended droplets that relax to cholesteric structure, \({{{{{\mathcal{S}}}}}}\) shows non-monotonic behavior where initially \({{{{{\mathcal{S}}}}}}\) decreases until its minimum, then it starts to increase before reaching its equilibrium, indicating a second-order exponential decay (Fig. 4c). To illustrate such a undershoot behavior more clearly, while maintaining \({{{{{\mathcal{S}}}}}}\) as a positively-defined object, we use in this case \({{{{{\mathcal{S}}}}}}=\frac{S\left(t\right)\,-\,{S}_{{{{{{\rm{minimum}}}}}}}}{{S}_{{{{{{\rm{init}}}}}}.}\,-\,{S}_{{{{{{\rm{minimum}}}}}}}}\) and show the rate of the changes in \({{{{{\mathcal{S}}}}}}\) versus the scaled time in Fig. 4d–f (note that the two definitions of \({{{{{\mathcal{S}}}}}}=\frac{S\left(t\right)\,-\,{S}_{{{{{{\rm{minimum}}}}}}}}{{S}_{{{{{{\rm{init}}}}}}.}\,-\,{S}_{{{{{{\rm{minimum}}}}}}}}\) and \({{{{{\mathcal{S}}}}}}=\frac{S\left(t\right)\,-\,{S}_{{{{{{\rm{equil}}}}}}.}}{{S}_{{{{{{\rm{init}}}}}}.}\,-\,{S}_{{{{{{\rm{equil}}}}}}.}}\) are equivalent for homogeneous and bipolar tactoids following a monotonic decay). It is clear from Fig. 4 that numerical simulations capture the first-order exponential decay in the case of homogenous and bipolar tactoids. In the case of cholesteric tactoids, our numerical simulation results bear qualitatively similar behavior for the structural relaxation of the tactoids although simulations underestimate the equilibrium order parameter obtained experimentally for cholesteric tactoids structural relaxation (Fig. 4c). To illustrate this further, we show the rate of the changes in \({{{{{\mathcal{S}}}}}}\) versus the scaled time in Fig. 4d–f. Quantitatively, the tactoids with the homogenous and bipolar configuration at relaxed state follow \({{{{{\mathcal{S}}}}}}={{{{{\rm{exp }}}}}}(-t/{\tau }_{{{{{{\rm{c}}}}}}})\), where we experimentally find the characteristic structural or configurational relaxation time τc to be 24.4 s and 31.9 s for homogenous (V = 644 µm3) and bipolar (V = 2751 µm3) tactoids, respectively. In contrast, the tactoid with cholesteric configuration at equilibrium (V = 16,414 µm3), shows two characteristic configurational relaxation times captured by \({{{{{\mathcal{S}}}}}}={{{{{{\rm{c}}}}}}}_{1}{{{{{\rm{exp }}}}}}\left(-t/{\tau }_{{{{{{\rm{c}}}}}},1}\right)\,+\,(1\,-\,{{{{{\rm{c}}}}}}_{1}){{{{{\rm{exp }}}}}}(-t/{\tau }_{{{{{{\rm{c}}}}}},2})\), with c1 a constant equal to −1.6, and τc,1 and τc,2 found to be 105.9 s and 333.7 s, respectively.

a–c Evaluation of \({{{{{\mathcal{S}}}}}}\) (defined as \(\tfrac{S\left(t\right)\,-\,{S}_{{{{{{\rm{equil}}}}}}.}}{{S}_{{{{{{\rm{init}}}}}}.}\,-\,{S}_{{{{{{\rm{equil}}}}}}.}}\) for homogeneous and bipolar tactoids and \({{{{{\mathcal{S}}}}}}=\frac{S\left(t\right)\,-\,{S}_{{{{{{\rm{minimum}}}}}}}}{{S}_{{{{{{\rm{init}}}}}}.}\,-\,{S}_{{{{{{\rm{minimum}}}}}}}}\) for cholesteric tactoids, where Sequil., Sinit., and S(t) are order parameter values at equilibrium, at the initial time, and at a given time t, respectively) with respect to scaled time, \(\frac{t}{{\tau }_{{{{{{\rm{c}}}}}}}}\) with τc the characteristic configurational relaxation time, for tactoids that relax to homogenous (a), bipolar (b) and cholesteric (c) configurations at equilibrium. The experimental insets showing the retardance images taken with LC-PolScope along with numerical simulation results present the critical state of the relaxation for each class of the tactoids. Color bar denotes the order parameter values in numerical simulation insets. Note that the brightness of the experimental images is increased for better visualization. The symbols denote the experimental data, black solid lines are numerical simulation results. Colored and dashed black lines show the fitting that is used to obtain τc from experimental and numerical simulation results, respectively. d–f The changes in \(\frac{{{{{{\rm{d}}}}}}{{{{{\mathcal{S}}}}}}}{{dt}}\), obtained from the fitted lines in a–c, during relaxation for different classes of tactoids: homogenous (d), bipolar (e), and cholesteric (f) configurations. While homogenous and bipolar tactoids follow monotonic single exponential decay during relaxation \({{{{{\mathcal{S}}}}}}={{{{{\rm{exp }}}}}}\left(-\frac{t}{{\tau }_{{{{{{\rm{c}}}}}}}}\right)\), the cholesteric tactoids are characterized by a non-monotonic behavior of \({{{{{\mathcal{S}}}}}}\) during relaxation (see (c)), described by a second-order exponential decay, \({{{{{\mathcal{S}}}}}}\,=\,{c}_{1}{{{{{\rm{exp}}}}}}\left(-t/{\tau }_{{{{{{\rm{c}}}}}},1}\right)\,+\,(1\,-\,c_{1}){{{{{\rm{exp }}}}}}(-t/{\tau }_{{{{{{\rm{c}}}}}},2})\), where c1 is a constant.

To interpret the physics behind these structural relaxations, we inspect the time scale related to bending, splay and twist terms in tactoids. Compared to the relaxation to homogeneous/bipolar tactoids, which involves nematic ordering with at best splay/bend relaxation, when the configuration relaxes to cholesteric tactoids, an additional twist relaxation takes place (see numerical simulation results in Supplementary Movies 4–6). This suggests that the second exponential decay associated with the cholesteric tactoids originates from the twist term. We suggest that in \({{{{{\mathcal{S}}}}}}={{{{{{\rm{c}}}}}}}_{1}{{{{{\rm{exp}}}}}}\left(-t/{\tau }_{{{{{{\rm{c}}}}}},1}\right)\,+\,(1\,-\,{{{{{\rm{c}}}}}}_{1}){{{{{\rm{exp }}}}}}(-t/{\tau }_{{{{{{\rm{c}}}}}},2})\), where τc,2 is significantly longer than τc,1, the second exponential decay τc,2 originates from the chiral twist term while τc,1 originates from simple nematic ordering. We base this statement on two grounds. First, we compare the length-scales of the nematic and cholesteric ordering. We consider the length scale of the cholesteric phase as the length which is required for the phase to form a single periodic twist, that is set by the inverse wave number (i.e., the pitch) which is in the order of 101 µm. In contrast, to form nematic ordering (characterized by splay and bending), the length scale is defined in the range of the length of the fibrils (mesogens) which is in the order of 10−1 µm. Thus, we argue that the larger twist length-scales compared to splay/bend imply longer relaxation times for the twisting deformation. Secondly, from the experiments of the relaxation of the cholesteric tactoids using a LC-Polscope, allowing us to capture the changes in the director field over time and follow the twist dynamics in the director field (see Supplementary Fig. 6), we find that change/rotation in the director field (the twist) takes place until the latest stages of the relaxation process and the twist changes are significant in the latest stage, suggesting again that the second exponential decay τc,2 originates from the twist re-arrangement. It is also worth mentioning that the structural relaxation time in tactoids is much lower than under static fixed boundary conditions. The structural relaxation time is of the order of hours for similar BLG and SCNC systems inside a capillary tube35 as opposed to tens/hundreds of seconds here. It is indeed known that boundary mobility increases the mobility of the mesogens thus promoting faster kinetics in structural relaxation12.

What is the configuration of the extended tactoid with a given initial volume after relaxation? We are able to predict the relaxed configuration of the tactoids using the theoretical modeling recently developed starting from a scaling form of Frank–Oseen elasticity theory4. According to this theory, the tactoids at equilibrium hold homogeneous configuration when (V/α)Homogenous < (K/γω)3, bipolar configuration when (K/γω)3 < (V/α)Bipolar < [1.7γ/(K2q∞2)]3, and cholesteric configuration when (V/α)Cholesteric > [1.7γ/(K2q∞2)]3. Approximating α equal to 3 for homogenous-bipolar and 1.5 for bipolar-cholesteric boundaries following ref. 4, we computed these threshold values for BLG I and found VHomogenous ≲ 800, 800 ≲ VBipolar ≲ 11,000, and VCholesteric ≳ 11,000 µm3. This confirms that the tactoids shown in Fig. 1 follow a relaxation path until equilibrium. Thus, knowing the initial volume of the tactoids, their configuration after relaxation can be predicted simply from the scaling form of Frank–Oseen elasticity theory and physical parameters of the system such as elastic constants, interfacial energy and anchoring strength4. We provide the nematic-cholesteric phase diagram of the tactoids collected from the samples in a cuvette at equilibrium showing tactoids configuration as a function of the volume in Supplementary Note 7, as a further demonstration that the initially stretched tactoids reach an equilibrium configuration after relaxation.

We have presented an integrated picture based on experiments, numerical simulations and theory, allowing a disentanglement and comprehensive description of the shape and structural relaxation of initially stretched colloidal liquid crystalline droplets. We have shown that these tactoids undergo relaxation on both shape and structure, when the external flow field maintaining them in a non-equilibrium state is released. Independently of the size, the shape relaxation of the tactoids is characterized by a single exponential decay, which is explained well by taking into account both isotropic and anisotropic features of the tactoids. In contrast, the structural relaxation follows different fates, with first- and second-order exponential decays, depending on the existence of splay, bend and twist contributions in the ground state, whose relative weight depend directly on the size of the tactoid. We have discussed the fundamental physical mechanisms behind the shape and structural relaxation and highlighted their interdependence. These results bring forward our understanding of dynamic processes in liquid crystalline systems based on filamentous colloids and introduce a combined experimental, theoretical and numerical formalism which can be extended to heterogeneous complex fluids, soft matter and biological colloids in general.

β-lactoglobulin was purified from whey protein following ref. 64 and dissolved in Milli-Q water at 2 wt%. The solution was cleared from aggregates by filtering using 0.45 µm Nylon syringe filter and pH of the solution was adjusted 2 by adding HCl. Later, the solution was heated for 5 hrs over a hot plate at 90 °C. Once the amyloid fibrils were prepared, we shortened the length of fibrils using the mechanical shear force method. Two sets of the solutions were prepared, see Supplementary Note 3 for details on length and height distributions. The solutions were dialyzed for 5 days using 100 kDa MWCO Spectra/Por dialysis membrane against pH 2 Milli-Q. The bath was changed every 24 h. To reach isotropic–nematic coexistence region concentration for the solution, suspensions up-concentrated using 6–8 kDa MWCO Spectra/Por 1 dialysis membrane against 6 wt% polyethylene glycol solution (mol wt: Mr ~ 35,000, Sigma Aldrich) in pH 2 milli-Q water. The solutions were kept in the fridge until phase separation happens, allowing reporting the isotropic and nematic phases concentrations in Table 1.

Cellulose nanocrystal suspensions were prepared by mixing freeze-dried cellulose nanocrystal (FPInnovations) in Milli-Q water. To make sure that cellulose nanocrystal dispersed well, the solution was ultra-sonicated for 120 s. This was followed by centrifugation for 20 min at 12,000 × g to remove aggregates. SCNC solution with concentration within the isotropic–nematic coexistence region was obtained by initially mixing 2.5 wt% freeze-dried cellulose nanocrystal in Milli-Q water.

To perform the AFM measurement, a droplet of diluted suspension (0.01 wt%) was deposited on freshly cleaved mica. After 2 min, the mica was rinsed with Milli-Q water and dried with air stream. The images of the sample were captured at ambient conditions with MultiMode VIII scanning probe microscope (Bruker) in tapping mode. The software FiberApp65 was used to analyze the images and measure the length and height distributions of fibrils.

We performed the experiments using an optical microscope Zeiss equipped with crossed polarizers and combined with LC-PolScope universal compensator. Under crossed polarizers, time-series images at the frame rate of 12 frames per minute were taken. The microfluidic channel was placed on the microscope in a way that the tactoids long axis held 45° angle with respect to one of the crossed polarizers. This allowed unambiguous measurement of the tactoids short and long axis during relaxation. To perform the measurements, MATLAB program and ImageJ software were used. Furthermore, we used optical microscopy combined with LC-PolScope universal compensator at time-series mode capturing 3 frames per minute. LC-PolScope images were used to analyze the internal structure of the tactoids. Additionally, LC-PolScope produces the retardance images giving the retardance value of the image, pixel by pixel, which were used to measure the order parameter.

The classical soft lithography approach was employed to make microfluidic systems66. We made PDMS by mixing polydimethylsiloxane (PDMS) monomer and curing agent (Dow Corning Slygard 184) with ratio 10 to 1. The plain glass slide (Corning 2947) was used as the base plate to attach the PDMS channel.

We used microfluidic system with rectangular cross-section with width of the channel at expansion zone we = 600 µm, the width of the contraction zone wc = 50 µm, and the height of the channel h = 100 µm (see Supplementary Information for the schematic of the microfluidic system).

We performed all the experiments at room condition. The equipment used to run experiments in microfluidics consists of a Harvard Apparatus syringe pump, 250 µl Hamilton syringe, flexible tubing with inner diameter 0.8 mm, and the needle with inner diameter 0.34 mm and outer diameter 0.64 mm.

Past studies23,28,67,68,69 have shown that the time-dependent Ginzburg–Landau model can capture the spatio-temporal coupled relaxation dynamics of shape and structure. Furthermore, our approach can result in high fidelity simulations capturing dynamics of spatio-temporal liquid crystalline self-assembly including breakage, coalescence, and defect evolution. However, these are out of scope of the current study. In the present study, we applied this modeling approach and we focused on the relaxation dynamics of the shape and structure in initially extended tactoids, which are isolated in isotropic phase. The implementation has been elaborated in detail in our previous works, see refs. 23,28,67. At the initial time, we consider an elongated tactoid according to the experimental observations. Thereafter, we let the elongated tactoids relax. Through relaxation, the total free energy is minimized according to the time-dependent Ginzburg–Landau model, by which the excess free energy such as surface tension and elasticity are relieved. The elongated tactoid self-selects the equilibrium shape and structure through a spontaneous thermodynamic-driven relaxation. Furthermore, the matrix surrounded the tactoid under relaxation is essentially kept at the isotropic concentration. Note that, in present work, we study three prime liquid crystalline configurations; homogenous nematic, bipolar nematic, uniaxial cholesteric, which are all fully rotationally symmetric70. Given this fact, to reduce computational costs, we rely on rectangular two-dimension simulations which provide a good 3D description since there is no need to discriminate between point, line and ring disclinations; see Supplementary Movies 4–6.

The data that support the findings of this study are available from the corresponding author upon request

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

Article ADS Google Scholar

Prinsen, P. & van der Schoot, P. Shape and director-field transformation of tactoids. Phys. Rev. E 68, 21701 (2003).

Article ADS CAS Google Scholar

Prinsen, P. & van der Schoot, P. Parity breaking in nematic tactoids. J. Phys. Condens. Matter 16, 8835 (2004).

Article ADS CAS Google Scholar

Nyström, G., Arcari, M. & Mezzenga, R. Confinement-induced liquid crystalline transitions in amyloid fibril cholesteric tactoids. Nat. Nanotech. 13, 330 (2018).

Article ADS CAS Google Scholar

Wang, P. X. & MacLachlan, M. J. Liquid crystalline tactoids: ordered structure, defective coalescence and evolution in confined geometries. Philos. Trans. A Math. Phys. Eng. Sci. 376, 20170042 (2018).

ADS PubMed Google Scholar

De Luca, A. et al. Optical nanotomography of anisotropic fluids. Nat. Phys. 4, 869–872 (2008).

Article CAS Google Scholar

Kim, J.-H., Yoneya, M. & Yokoyama, H. Tristable nematic liquid-crystal device using micropatterned surface alignment. Nature 420, 159–162 (2002).

Article ADS CAS PubMed Google Scholar

Woltman, S. J., Jay, G. D. & Crawford, G. P. Liquid-crystal materials find a new order in biomedical applications. Nat. Mater. 6, 929–938 (2007).

Article ADS CAS PubMed Google Scholar

Almohammadi, H., Bagnani, M. & Mezzenga, R. Flow-induced order–order transitions in amyloid fibril liquid crystalline tactoids. Nat. Commun. 11, 1–9 (2020).

Article CAS Google Scholar

Velev, O. D., Lenhoff, A. M. & Kaler, E. W. A class of microstructured particles through colloidal crystallization. Science 287, 2240–2243 (2000).

Article ADS CAS PubMed Google Scholar

Whitesides, G. M. & Grzybowski, B. Self-assembly at all scales. Science 295, 2418–2421 (2002).

Article ADS CAS PubMed Google Scholar

Hunter, G. L., Edmond, K. V. & Weeks, E. R. Boundary mobility controls glassiness in confined colloidal liquids. Phys. Rev. Lett. 112, 218302 (2014).

Article ADS CAS Google Scholar

Shuang, Z., Sokolov, A., Lavrentovich, O. D. & Aranson, I. S. Living liquid crystals. Proc. Natl Acad. Sci. USA 111, 1265–1270 (2014).

Article ADS CAS Google Scholar

Marchetti, M. C. et al. Hydrodynamics of soft active matter. Rev. Mod. Phys. 85, 1143–1189 (2013).

Article ADS CAS Google Scholar

Azzari, P., Bagnani, M. & Mezzenga, R. Liquid-liquid crystalline phase separation in biological filamentous colloids: nucleation, growth and order-order transitions of cholesteric tactoids. Soft Matter 17, 6627–6636 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Zocher, H. Über freiwillige Strukturbildung in Solen. (Eine neue Art anisotrop flüssiger Medien.). Z. f.ür. anorganische und Allg. Chem. 147, 19251470111 (1925).

Google Scholar

Bawden, F. C., Pirie, N. W., Bernal, J. D. & Fankuchen, I. Liquid crystalline substances from virusinfected plants. Nature 138, 1051–1052 (1936).

Article ADS Google Scholar

Dogic, Z. Surface freezing and a two-step pathway of the isotropic-smectic phase transition in colloidal rods. Phys. Rev. Lett. 91, 165701 (2003).

Article ADS PubMed CAS Google Scholar

Oakes, P. W., Viamontes, J. & Tang, J. X. Growth of tactoidal droplets during the first-order isotropic to nematic phase transition of F-actin. Phys. Rev. E 75, 061902 (2007).

Article ADS CAS Google Scholar

Puech, N., Grelet, E., Poulin, P., Blanc, C. & Van Der Schoot, P. Nematic droplets in aqueous dispersions of carbon nanotubes. Phys. Rev. E 82, 020702 (2010).

Article ADS CAS Google Scholar

Bagnani, M., Azzari, P., Assenza, S. & Mezzenga, R. Six-fold director field configuration in amyloid nematic and cholesteric phases. Sci. Rep. 9, 1–9 (2019).

Article CAS Google Scholar

Nyström, G., Arcari, M., Adamcik, J., Usov, I. & Mezzenga, R. Nanocellulose fragmentation mechanisms and inversion of chirality from the single particle to the cholesteric phase. ACS Nano 12, 5141–5148 (2018).

Article PubMed CAS Google Scholar

Khadem, S. A. & Rey, A. D. Nucleation and growth of cholesteric collagen tactoids: A time-series statistical analysis based on integration of direct numerical simulation (DNS) and long short-term memory recurrent neural network (LSTM-RNN). J. Colloid Interface Sci. 582, 859–873 (2021).

Article ADS CAS PubMed Google Scholar

Li et al. Colloidal cholesteric liquid crystal in spherical confinement. Nat. Commun. 7, 12520 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

Parker, R. M. et al. Hierarchical self-assembly of cellulose nanocrystals in a confined geometry. ACS Nano 10, 8443–8449 (2016).

Article CAS PubMed PubMed Central Google Scholar

Verhoeff, A. A., Bakelaar, I. A., Otten, R. H. J., van der Schoot, P. P. A. M. & Lekkerkerker, H. N. W. Tactoids of plate-like particles: Size, shape, and director field. Langmuir 27, 116–125 (2011).

Article CAS PubMed Google Scholar

Rey, A. D. & Denn, M. M. Dynamical phenomena in liquid-crystalline materials. Ann. Rev. Fluid Mech. 34, 233–266 (2002).

Article ADS MathSciNet MATH Google Scholar

Rey, A. D. Liquid crystal models of biological materials and processes. Soft Matter 6, 3402–3429 (2010).

Article ADS CAS Google Scholar

Larson, R. G. The Structure and Rheology of Complex Fluids (Oxford University Press. 1999).

Gârlea, I. C. et al. Finite particle size drives defect-mediated domain structures in strongly confined colloidal liquid crystals. Nat. Commun. 7, 12112 (2016).

Article ADS PubMed PubMed Central CAS Google Scholar

Tortora, L. & Lavrentovich, O. D. Chiral symmetry breaking by spatial confinement in tactoidal droplets of lyotropic chromonic liquid crystal. Proc. Natl Acad. Sci. USA 108, 5163–5168 (2011).

Article ADS CAS PubMed PubMed Central Google Scholar

Dammone, O. J. et al. Confinement induced splay-to-bend transition of colloidal rods. Phys. Rev. Lett. 109, 108303 (2012).

Article ADS PubMed CAS Google Scholar

Erdmann, J. H., Žumer, S. & Doane, J. W. Configuration transition in a nematic liquid crystal confined to a small spherical cavity. Phys. Rev. Lett. 64, 1907 (1990).

Article ADS CAS PubMed Google Scholar

Gârlea, I. C. et al. Colloidal liquid crystals confined to synthetic tactoids. Sci. Rep. 9, 20391 (2019).

Article ADS PubMed PubMed Central CAS Google Scholar

Khadem, S. A., Bagnani, M., Mezzenga, R. & Rey, A. Relaxation dynamics in bio-colloidal cholesteric liquid crystals confined to cylindrical geometry. Nat. Commun. 11, 1–10 (2020).

Article CAS Google Scholar

Taylor, G. I. The formation of emulsions in definable fields of flow. Proc. R. Soc. Lond. A 146, 501–523 (1934).

Article ADS CAS Google Scholar

Cabral, J. T. & Hudson, S. D. Microfluidic approach for rapid multicomponent interfacial tensiometry. Lab. Chip 6, 427–436 (2006).

Article CAS PubMed Google Scholar

Brosseau, Q., Vrignon, J. & Baret, J.-C. Microfluidic dynamic interfacial tensiometry (μDIT). Soft Matter 10, 3066–3076 (2014).

Article ADS CAS PubMed Google Scholar

Weirich, K. L. et al. Liquid behavior of cross-linked actin bundles. Proc. Natl Acad. Sci. USA 114, 2131–2136 (2017).

Article ADS CAS PubMed PubMed Central Google Scholar

Rallison, J. M. The deformation of small viscous drops and bubbles in shear flows. Annu. Rev. Fluid Mech. 16, 45–66 (1984).

Article ADS MATH Google Scholar

Hudson, S. D., Cabral, J. T., Goodrum, W. J., Beers, K. & Amis, E. Microfluidic interfacial tensiometry. Appl. Phys. Lett. 87, 081905 (2005).

Article ADS CAS Google Scholar

Das, S. K. & Rey, A. D. Magnetic field-induced shape transitions in multiphase polymer-liquid crystal blends. Macromolecular theory and simulations. Macromol. Theory Simul. 15, 469–489 (2006).

Article ADS CAS Google Scholar

Odijk, T. Elastic constants of nematic solutions of rod-like and semi-flexible polymers. Liq. Cryst. 1, 553–559 (1986).

Article CAS Google Scholar

Ericksen, J. L. Inequalities in liquid crystal theory. Phys. Fluids 9, 1205 (1966).

Article ADS Google Scholar

Dietrich, C. F., Collings, P. J., Sottmann, T., Rudquist, P. & Giesselmann, F. Extremely small twist elastic constants in lyotropic nematic liquid crystals. Proc. Natl Acad. Sci. USA 117, 27238–27244 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Zhou, S. et al. Elasticity of lyotropic chromonic liquid crystals probed by director reorientation in a magnetic field. Phys. Rev. Lett. 109, 037801 (2012).

Article ADS PubMed CAS Google Scholar

Taratuta, V. G., Hurd, A. J. & Meyer, R. B. Light-scattering study of a polymer nematic liquid crystal. Phys. Rev. Lett. 55, 246 (1985).

Article ADS CAS PubMed Google Scholar

Bagnani, M., Azzari, P., De Michele, C., Arcari, M. & Mezzenga, R. Elastic constants of biological filamentous colloids: estimation and implications on nematic and cholesteric tactoid morphologies. Soft Matter 17, 2158–2169 (2021).

Article ADS CAS PubMed Google Scholar

Doi, M. & Edwards, S. F. The Theory of Polymer Dynamics. Vol. 73 (Clarendon Press, 1988).

Onsager, L. The effect of shape on the interaction of colloidal particles. Ann. NY Acad. Sci. 51, 627–659 (1949).

Article ADS CAS Google Scholar

Vroege, G. J. & Odijk, T. Elastic moduli of a nematic liquid-crystalline solution of polyelectrolytes. J. Chem. Phys. 87, 4223–4232 (1987).

Article ADS CAS Google Scholar

Mezzenga, R., Jung, J. M. & Adamcik, J. Effects of charge double layer and colloidal aggregation on the isotropic-nematic transition of protein fibers in water. Langmuir 26, 10401–10405 (2010).

Article CAS PubMed Google Scholar

Usov, I. et al. Understanding nanocellulose chirality and structure-properties relationship at the single fibril level. Nat. Commun. 6, 7564 (2015).

Article ADS CAS PubMed Google Scholar

Nyström, G., Fong, W. K. & Mezzenga, R. Ice-templated and cross-linked amyloid fibril aerogel scaffolds for cell growth. Biomacromolecules 18, 2858–2865 (2017).

Article PubMed CAS Google Scholar

Zhou, S., Cervenka, A. J. & Lavrentovich, O. D. Ionic-content dependence of viscoelasticity of the lyotropic chromonic liquid crystal sunset yellow. Phys. Rev. E 90, 042505 (2014).

Article ADS CAS Google Scholar

Adamcik, J. et al. Understanding amyloid aggregation by statistical analysis of atomic force microscopy images. Nat. Nanotechnol. 5, 423–428 (2010).

Article ADS CAS PubMed Google Scholar

van der Schoot, P. Remarks on the interfacial tension in colloidal systems. J. Phys. Chem. B 103, 8804–8808 (1999).

Article CAS Google Scholar

Koch, D. L. & Harlen, O. G. Interfacial tension at the boundary between nematic and isotropic phases of a hard rod solution. Macromolecules 32, 219–226 (1999).

Article ADS CAS Google Scholar

Milliken, W. J. & Leal, L. G. Deformation and breakup of viscoelastic drops in planar extensional flows. J. Non-Newton. Fluid Mech. 40, 355–379 (1991).

Article CAS Google Scholar

Bagnani, M., Nyström, G., De Michele, C. & Mezzenga, R. Amyloid fibrils length controls shape and structure of nematic and cholesteric tactoids. ACS Nano 13, 591–600 (2019).

Article CAS PubMed Google Scholar

Stone, H. A. Dynamics of drop deformation and breakup in viscous fluids. Annu. Rev. Fluid Mech. 26, 65–102 (1994).

Article ADS MathSciNet MATH Google Scholar

Oldenbourg, R., Salmon, E. & Tran, P. Birefringence of single and bundled microtubules. Biophys. J. 74, 645–654 (1998).

Article ADS CAS PubMed PubMed Central Google Scholar

Gentry, B., Smith, D. & Käs, J. Buckling-induced zebra stripe patterns in nematic F-actin. Phys. Rev. E 79, 031916 (2009).

Article ADS CAS Google Scholar

Vigolo, D. et al. Continuous isotropic-nematic transition in amyloid fibril suspensions driven by thermophoresis. Sci. Rep. 7, 1211 (2017).

Article ADS PubMed PubMed Central CAS Google Scholar

Usov, I. & Mezzenga, R. FiberApp: an open-source software for tracking and analyzing polymers, filaments, biomacromolecules, and fibrous objects. Macromolecules 48, 1269–1280 (2015).

Article ADS CAS Google Scholar

Duffy, D. C., McDonald, J. C., Schueller, O. J. & Whitesides, G. M. Rapid prototyping of microfluidic systems in poly (dimethylsiloxane). Anal. Chem. 70, 4974–4984 (1998).

Article CAS PubMed Google Scholar

Khadem, S. A. & Rey, A. D. Theoretical platform for liquid-crystalline self-assembly of collagen-based biomaterials. Front. Phys. 7, 88 (2019).

Article Google Scholar

Gutierrez, O. F. A. & Rey, A. D. Biological plywood film formation from para-nematic liquid crystalline organization. Soft Matter 13, 8076–8088 (2017).

Article ADS Google Scholar

Hohenberg, P. C. & Halperin, B. I. Theory of dynamic critical phenomena. Rev. Mod. Phys. 49, 435–479 (1977).

Article ADS CAS Google Scholar

Seč, D., Porenta, T., Ravnik, M. & Žumer, S. Geometrical frustration of chiral ordering in cholesteric droplets. Soft Matter 8, 11982–11988 (2012).

Article ADS CAS Google Scholar

Download references

We thank X.Cao (ETHZ) for help in microfluidic chips fabrications and Samuel Mathews for support to maintain in-house high-performance supercomputer. We thank Prof. Andrew de Mello (ETHZ) for granting access to his laboratory and Y.Yuan for helpful discussions. This work is supported by Sinergia grant no. CRSII5_189917 from the Swiss National Science foundation (R.M.).

These authors contributed equally: Hamed Almohammadi, Sayyed Ahmad Khadem.

Department of Health Sciences and Technology, ETH Zurich, Zurich, Switzerland

Hamed Almohammadi, Massimo Bagnani & Raffaele Mezzenga

Department of Chemical Engineering, McGill University, Montreal, QC, Canada

Sayyed Ahmad Khadem & Alejandro D. Rey

Quebec Centre for Advanced Materials, Canada (QCAM/CQMF), Montreal, QC, Canada

Sayyed Ahmad Khadem & Alejandro D. Rey

Department of Materials, ETH Zurich, Zurich, Switzerland

Raffaele Mezzenga

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

H.A. and R.M. conceived and initiated the project, designed the experiments, analyzed the data, developed theoretical modeling of the tactoids deformation, contributed to the theoretical analysis of relaxation time of the tactoids, and wrote the bulk of the paper. H.A. built the experimental apparatus and performed the experiments. M.B. contributed to the experiments and carried out the AFM measurements. S.A.K. and A.D.R. designed the simulations, analyzed the data, developed the theoretical analysis of relaxation time of the tactoids, and contributed to the writing of the paper. S.A.K. performed the simulations. A.D.R. and R.M. supervised the research. All authors discussed and edited the paper.

Correspondence to Raffaele Mezzenga.

The authors declare no competing interests.

Nature Communications thanks Mehrana Raeisin Nejad and the other anonymous reviewer(s) for their contribution to the peer review of this work. Peer review reports are available.

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

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Almohammadi, H., Khadem, S.A., Bagnani, M. et al. Shape and structural relaxation of colloidal tactoids. Nat Commun 13, 2778 (2022). https://doi.org/10.1038/s41467-022-30123-y

Download citation

Received: 10 December 2021

Accepted: 19 April 2022

Published: 19 May 2022

DOI: https://doi.org/10.1038/s41467-022-30123-y

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

Nature Communications (2023)

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.