
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

38556534
58054
10.1038/s41598-024-58054-2
Article
Morphogenesis of a chiral liquid crystalline droplet with topological reconnection and Lehmann rotation
Yoshioka Jun j-yoshi@fc.ritsumei.ac.jp

Ito Yuki
Fukao Koji
https://ror.org/0197nmd03 grid.262576.2 0000 0000 8863 9909 Department of Physical Sciences, Ritsumeikan University, 1-1-1 Noji-Higashi, Kusatsu, Shiga 525-8577 Japan
31 3 2024
31 3 2024
2024
14 759726 7 2023
25 3 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ 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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.
Morphogenesis is a hierarchical phenomenon that produces various macroscopic structures in living organisms, with high reproducibility. This study demonstrates that such structural formation can also be observed in a chiral liquid crystalline droplet under a temperature gradient. Through specific control of the temperature change process, we were able to switch the final structure obtained as a result of the formation via the appearance and reconnection of loop defects in the transient state during structure formation. Simultaneously, the existence of the gradient resulted in a characteristic rotational phenomenon called Lehmann rotation, which was prominently induced in the transient state. By demonstrating three-dimensional measurements of the flow field, we revealed the existence of Marangoni convection in the state. Consequently, it is indicated that the convection results in high-speed Lehmann rotation and large structural deformation with topological changes, thereby playing a significant role in the structure formation.

Subject terms

Liquid crystals
Fluid dynamics
JSPS KAKENHI22H01191 Yoshioka Jun issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Morphologies of living organisms are formed via ordering and disordering processes with various size scales, as observed in the segmentation and differentiation of eggs, or the metamorphosis of insects1,2. Morphogenesis is a highly controlled phenomenon that generates various macroscopic structures with a high reproducibility. A key factor in this phenomenon is the appearance of a concentration gradient, indicating that a potential gradient is required for morphogenesis. The present study demonstrates that such a structural formation process under a gradient can also be observed in liquid crystalline (LC) systems, which exist as components of living organisms.

An example of the macroscopic structure formed by LC has been observed in droplets composed of nematic (N), smectic, or cholesteric (Ch) phase3–20. Notably, the chiral LC system of the Ch droplets exhibits a wide range of structural diversity. Various structures with three-dimensional helices and point, line, and loop defects are formed, depending on the droplet size, helical pitch length, elastic constants, and anchoring constants3,11–20. However, even when these parameters, which are considered to determine an equilibrium state, are successfully controlled, the realised structure may not be uniquely determined in real systems. Reportedly, the internal director field in the droplet can show different states, although the experimental conditions are controlled to be identical13,15,18,20. One may consider that this structural discrepancy is attributed to variations in the structure formation process. Consequently, it is hypothesised that control of the process should result in complete structural control of the Ch droplet. In this study, we fabricated disk-shaped droplets sandwiched by substrates via a cooling process from an isotropic liquid (I) phase under a constant temperature gradient. The existence of this gradient enables gradual growth of the domain with LC order from the low-temperature side of the droplet21. This controlled growth should enable control of the final structure of the Ch droplets obtained by the cooling process.

The existence of a gradient, such as a thermal, electric, or chemical potential gradient in the Ch LC, may yield a rotation of the director field; this characteristic phenomenon is called Lehmann rotation22–25. The phenomenon was first discovered by Lehmann in the Ch LC under a temperature gradient in 1900 but was not reproduced for more than a century after the discovery. In this situation, as a reproduction of the Lehmann rotation, a steady rotation was observed in a Ch droplet under a temperature gradient by Oswald et al. in 200826. Subsequently, this phenomenon was reproduced by several research groups18,27–35. As the mechanism of the rotational phenomenon, two possibilities have been proposed. The first is based on Leslie’s theory, which predicts the existence of direct coupling between director rotation and the temperature gradient in Ch LC23,36. The other mechanism is attributed to Marangoni convection resulting from the surface tension gradient; the rotation is driven by the material flow in the twisted director field31,33. To date, the rotational phenomenon has mainly been observed in droplets dispersed in specific liquid solvents18,26–35. In contrast, we used a system without solvents; a droplet directly exposed to air was used in this study. However, despite this difference, Lehmann rotation should also be observed in our system if the rotation is driven by any of the aforementioned mechanisms.

In this study, we demonstrate the structural formation processes of N and Ch droplets with an air interface in the presence of a temperature gradient. These processes are accompanied by characteristic dynamics, such as structural deformation, Lehmann rotation, and Marangoni convection. A comparison of the observation results for the Ch and N droplets clarifies the contribution of the existence or nonexistence of chirality to the dynamics. Based on experimental analyses of the director and flow fields, we report a possible description of the structure formation in these droplets.

Results

Structure formation of nematic and cholesteric droplets

First, by creating an N droplet, we applied a temperature gradient (see Fig. 1(v), the Methods section and Supplementary Note 2). In this study, T0 is defined as the temperature at the droplet centre, and ΔT is the temperature difference between the upper and lower glass substrates (For more detail about the determination of T0 and ΔT, see Supplementary Note 3 and ref.37). Keeping ΔT constant, we changed T0 according to the chart in Fig. 1(w), which consists of two cooling and heating processes. Starting from the I phase, we decreased T0 and observed the textural changes in the droplet using polarised microscopy (POM). In the 1st cooling process, a cross-texture was observed in the coexisting state of the I and N phases (I + N phase), as shown in Fig. 1b–e. This indicates that the director field has a radial alignment and a point defect exists in the droplet18,23,38. During the cooling process, a point defect first appeared at the droplet centre (Fig. 1b) and then moved away from the centre at a lower temperature (Fig. 1d). Subsequently, the defect returned to the droplet centre (Fig. 1e), and this configuration was preserved until the entire droplet transited into the N phase. A centred cross texture was also observed after the transition (Fig. 1f); the radial director alignment was also formed in the N phase39, as well as in the I + N phase.Figure 1 Structural formation of an N droplet. Temperature difference between upper and lower glass substrates was ΔT = 10K, and heating/cooling rate was 0.2 K/min. (a–u) show POM images at each T0, defined as the temperature at the droplet centre. A and P in (a) indicate analyser and polariser respectively, and the white bar in (b) 100 μm. Schematic image of experimental system is shown in (v), and the heating and cooling chart is in (w). The chart consists of two cooling and heating processes, which are called 1st and 2st coolings, and 1st and 2st heatings, respectively. The corresponding movies are available in Supplementary Videos 1–3.

After the 1st cooling process, the temperature was increased. In this 1st heating process, the point defect was located at the droplet centre at a lower T0 of less than ~ 40°C and moved away from the centre at a higher T0 as shown in Fig. 1g–k, as well as in the 1st cooling process. The 1st heating process ended with the off-centre state of the defect (Fig. 1k), after which the temperature decreased. In this 2nd cooling process, the point defect returned to the droplet centre in the I + N phase, as shown in Fig. 1l–n. Finally, a cross texture similar to that observed in the 1st cooling process was formed (Fig. 1o). Subsequently, the temperature was increased. In this 2nd heating process, the movement of the point defect similar to that in the 1st heating was observed, as shown in Fig. 1g–j and p–s. Finally, the point defect returned to the droplet centre, and the droplet transitioned into the I phase (Fig. 1t and u).

By defining R as the droplet radius and Rc as the distance of the point defect from the droplet centre, we measured T0 dependence of RcRcRR as shown in Fig. 3a. The temperature region wherein RcRcRR shows a non-zero value depends on whether the temperature protocol is cooling or heating, suggesting the existence of thermal hysteresis. RcRcRR in the 2nd cooling process agrees with RcRcRR in the 1st heating when T0 is relatively high (T0>∼43∘C), and agrees with RcRcRR in the 1st cooling when T0 is relatively low (T0<∼43∘C). Owing to the agreement in the cooling processes, the final structures obtained in the N phase in the 1st and the 2nd coolings were identical, as shown in Fig. 1f and o.

Subsequently, we observed the Ch droplets in a manner similar to that described previously. In the 1st cooling process, the centred cross texture was first formed in the coexistence phase (I + Ch phase), as shown in Fig. 2b. The appearance of the texture indicates the formation of a point defect as well as an N droplet. The defect moved away from the centre during the cooling process, as shown in Fig. 2c. In contrast to N droplet, the point defect exhibited a clockwise circular motion (see Supplementary Video 4). As the temperature decreased further, the defect returned to the droplet centre, and the circular motion stopped (Fig. 2d). This configuration was preserved just before the entire droplet transitioned into the Ch phase. As the transition occurred, a comma-like texture formed near the droplet centre, as shown in Fig. 2f. In this study, we refer to a Ch droplet with this texture as the ‘Type-A’ state.Figure 2 Structural formation of a Ch droplet. Temperature difference between upper and lower glass substrates was ΔT = 10K, and heating/cooling rate was 0.2 K/min. (a–u) show POM images at each T0, defined as the temperature at the droplet centre. A and P in (a) indicate analyser and polariser respectively, and the white bar in (b) 100μm. Schematic image of experimental system is shown in (v), and the heating and cooling chart is in (w). The chart consists two cooling and heating processes, which are called 1st and 2st coolings, and 1st and 2st heatings, respectively. The corresponding movies are available in Supplementary Videos 4–6.

In the 1st heating process, the comma-like texture was preserved just after the transition into the I + Ch phase, as shown in Fig. 2h, whereas at higher T0 the texture deformed into a cross-like texture, again indicating the formation of a point defect at the droplet centre (Fig. 2i). Subsequently, the defect moved away from the centre and exhibited a clockwise circular motion, as shown in Fig. 2j and Supplementary Video 5. The 1st heating process ended with an off-centre circulating state of the point defect (Fig. 2k). In the following 2nd cooling process, the cross texture gradually deformed into a comma-like texture, as shown in Fig. 2k–n, and it did not return to the droplet centre. The circular motion continued after the texture deformed into a comma-like shape, whereas its direction changed from clockwise to counterclockwise (see Supplementary Video 5). As the transition into the Ch phase occurred, the comma-like texture remained near the droplet edge, as shown in Fig. 2o; we refer to this state as ‘Type-B’, distinguished from Type-A (Fig. 2f). In contrast to N droplet, the final structure of Ch droplet depends on whether it is obtained during the 1st or 2nd cooling process. Since both type-A and -B states hold under the absence of the temperature gradient, they are considered as the stable or one of the metastable structures at equilibrium condition (see Supplementary Note 4). In the 2nd heating process, the comma-like texture gradually deformed into a cross texture, and the direction of the circular motion changed from counterclockwise to clockwise (Fig. 2p–s and Supplementary Video 6). Subsequently, the cross-texture defect moved to the droplet centre before the transition into the I phase (Fig. 2t and u).

T0 dependence of RcRcRR on the Ch droplets is shown in Fig. 3b. It demonstrates thermal hysteresis in the 1st cooling and heating processes, as well as in the N droplet. However, the RcRcRR in the 2nd cooling did not agree RcRcRR in the 1st heating or cooling process, and the point defect gradually changed into a structure with a comma-like texture, as described previously (Fig. 2k–n). After the structural change we could no longer define Rc, so that the data is absent in the region with comma-like texture in Fig. 3b. These hysteresis behaviours were not observed in the N droplet, and the observation results in Figs. 2 and 3b show that the structure formation in the Ch droplet followed different kinetic pathways in the 1st and the 2nd cooling processes, respectively. This difference in the pathway results in structural differences observed in the Ch phase (Type-A and Type-B states).Figure 3 Temperature dependences of RcRcRR in N (a) and Ch (b) droplets. R is defined as the droplet radius, and Rc as the distance between the point defect and the droplet centre; the definitions of them are illustrated in the insets of (a) and (b). The measurements of RcRcRR were performed under the temperature (T0) chart shown in Figs. 1(w) and 2(w), consisting of 1st cooling and heating, and 2nd cooling and heating processes, where ΔT was 10K. In (b), the data is not shown when T0 is lower than 43 °C in the 2nd cooling and heating processes. In this period, the point defect transformed into the structure with the comma-like texture (Fig. 2(l)–(n) and (q)–(s)), where Rc cannot be well-defined. (c) is the T0 dependence of steady-state solution of RcRcRR and the absolute value of τ, obtained by theoretical analysis. h and R were assumed to be 100 μm and 150 μm, respectively, and Vs was 1.0 μm/s. The rotational viscosities were assumed as γ1 = γ2 = 0.030 Pa·s40,41, and the elastic constants were K1 = K3 = 4.0 pN and K2 = 2.0 pN42.

The cross texture observed in the Ch droplets was twisted. In the 1st cooling process, this twist drastically progressed in the temperature range between T0 = 44.1 and 43.8 °C as shown in Fig. 4a. In addition, by observing the droplet without polarisers, we discovered that double-ring textures were formed just after the drastic progress of the twist, as shown in Fig. 4a–c and Supplementary Videos 4 and 7. The appearance of rings suggests the formation of double-loop defects, as reported in the literature13,16.Figure 4 Observation of the formation, deformation and reconnection of defect structures. (a), (d) and (g) are POM images, and the others are the observation results without polarisers. Red symbol in (a) indicates the direction of the temperature gradient. A and P in (a) indicate analyser and polariser respectively, and the white bar in (b) 100 μm. In (a, b, d, e, g and h), temperature difference in each figure is 0.1K and the cooling/heating rate 0.2K/min, so that the time difference in each figure is 30 s. (c and f) are displayed with higher time resolution (3 s each), showing the detailed change of the texture when the loop defect is formed and reconnected. During the 1st cooling process, the double loop defects, indicated by the aqua arrows in (b), are formed. In the end of the 1st heating process, the double loop defects are transformed into a single loop, as shown by the orange arrow in (e). In this instance the point defect (green arrow) is outside of the loop. After that, the point defect goes inside the single loop as indicated by the purple arrow. The corresponding movies are available in Supplementary Videos 4–8.

Reconnection of the defect structure was observed at the end of the 1st heating process. As shown in Fig. 4e and f, the double loops changed into a single loop in the temperature between T0 = 44.1 and 44.2 °C (see also Supplementary Videos 5 and 8). Immediately after reconnection, the point defect observed as a dark point was outside the single-loop defect, whereas the defect then entered the inside of the loop, as shown in Fig. 4e. Thus, a topological change in the defect structure was induced, resulting in a structural difference in the Ch phase (type-A and B states). The single loop was gradually elongated with a point defect inside the loop in the 2nd cooling process (Fig. 4h), and it never followed the structural change in the 1st heating or cooling processes.

Director rotation in cholesteric droplet under temperature gradient

Director rotation was driven in the Ch droplet not only under an unsteady state with a change in T0 but also under a steady state. Keeping both T0 and ΔT constant, we observed that the steady clockwise rotation of the director field was driven in both the Type-A and Type-B Ch states, as shown in Fig. 5a and f (see also Supplementary Video 9(a) and (f)). Since the structural change was not observed during the rotation, we should consider that these rotating droplets were in the steady state. The rotational speed was proportional to ΔT (see Supplementary Note 5), as is often reported for Ch droplets with Lehmann rotation18,26–28,35. This indicates that the rotations observed in the Ch phase were driven by the application of a temperature gradient.Figure 5 Director rotations under steady states. The time evolutions of POM images are shown. A and P in (a) indicate analyser and polariser respectively, and the white bar 100 μm. Red symbol in (f) indicates the direction of the temperature gradient. ΔT is set to be 10 K, and T0 is 35.0 °C in (a) and (f), 37.5 °C in (b) and (g), 38.5 °C in (c) and (h), 42.0 °C in (d) and (i) and 43.0 °C in (e) and (j), respectively. (a) and (f) are observed in the Ch phase, and the others are in the I + Ch phase. Time interval in each figure is 3000 s in (a) and (f), 300 s in (b), (c), (g) and (h) and 30 s in (d), (e), (i) and (j). Clockwise rotation was observed in (a), (d)–(f), (i) and (j), while counter-clockwise was in (b), (g) and (h). The corresponding movies are available in Supplementary Video 9.

Subsequently, keeping ΔT constant and increasing T0, we obtained the I + Ch phases for Type-A and Type-B, respectively. In these situations, we observed the droplets under constant T0 and ΔT. Just above the transition temperature from the Ch to the I + Ch phase, rotation was observed in both the Type-A and B as shown in Fig. 5b and g. Herein, both showed counter-clockwise rotation, reversed from the rotation observed in the Ch phase (see Supplementary Video 9(a), (b), (f), and (g)). In addition, the rotational speed in the I + Ch phase is one order of magnitude higher than that in the Ch phase (closed symbols with warm colours and open symbols with cold colours in Fig. 6a and b). As T0 increased further, in the case of Type-A, a four-fold symmetric texture with the centred point defect was formed and the rotation stopped, whereas in Type-B, the structure with the off-centred defect was preserved and the director field continued rotating, as shown in Fig. 5c and h (see also Supplementary Video 9(c) and (h)). With increasing T0, the rotational direction again reversed: in both cases of Type-A and B, the clockwise rotation was observed as shown in Fig. 5d, e, i and j and Supplementary Video 9(d), (e), (i) and (j); simultaneously, the rotational speed increased by one order of magnitude (open symbols with cold and warm colours in Fig. 6a and b). Heating the droplet further from this state, we observed that the four-fold symmetric texture shown in Fig. 2b was formed in both Types-A and B; in this state, rotation was never observed.Figure 6 Temperature dependence of rotational speed and state diagram. (a) and (b) show the dependences of the rotational speed on T0 and ΔT in the case of Type-A and B, respectively. The symbols with warm colours indicate the appearance of clockwise (CW) rotation, and cold colours indicate counter-clockwise (CCW). The data obtained in the Ch and I + Ch phases are shown by closed and open symbols, respectively. (c) and (d) are the state diagrams under constant T0 and ΔT in the case of Type-A and B Ch droplets, respectively. Based on the observation images under crossed polarizers and without polarisers, the structures were categorised into five groups as shown in (e) and (f). The structure in the Ch phase is categorised into Type-A and B, which are obtained at the end of the 1st and 2nd cooling processes, respectively. The I + Ch phase is categorised into three groups: DL, SL and P states. Only a point defect exists in the P state, while in the DL and SL states double and single loop defects are formed, respectively. The point defect can exist also in the DL and SL states, while it transforms into the structure with comma-shaped texture at lower temperature (Fig. 5 (b) and (g)–(i)). The observation was first performed under a constant ΔT in the Type-A or Type-B state, and the following observations were performed in stepwise increase of T0. Changing ΔT and repeating these observations, we obtained the results of (a), (b), (c) and (d). Here, the red and blue symbols indicate that CW and CCW rotations were observed respectively, and the black symbols indicate that no rotation was observed. In the I + Ch phase in (d), the direction of the rotation reverses at the broken line: this indicates that the direction is controlled by T0.

By changing ΔT and repeating the above observation, we created state diagrams for the rotation, as shown in Fig. 6c and d. The red and blue symbols indicate that clockwise and counter-clockwise rotation was observed in the diagrams, respectively, and the black symbols indicate that no rotation was under the corresponding conditions. In addition, based on the structural difference of the droplet, the states in the I + Ch phase were categorised into the three groups: DL state with the double-loop defects, SL state with a single-loop defect, and the P state with only a point defect. In the regions of no rotation (black symbols), four-fold symmetric textures, as shown in Fig. 2b and d, were always observed. The appearance of this texture under the POM observation indicates the existence of the director field with circular symmetry18,23,38,39. In this situation, director rotation cannot be driven unless the director field is deformed. Thus, director rotation was always observed in this experiment unless a structure with circular symmetry was formed.

Notably, the rotational direction is determined by T0 in the I + Ch phase, as seen in the type-B state diagram of Fig. 6d. Herein, we considered that this is attributed to the peculiar property of the surface tension near the phase transition temperature, and to the fact that the direction of the Marangoni flow depends on the surface tension gradient43–45. In the N LC material of 7CB used in this study (see Methods section), it was found that the flow direction depends on the temperature owing to the drastic change in surface tension near the I–N transition point21,46. By measuring the surface tension of the Ch LC sample used in this study, we confirmed that its temperature dependence is similar to that of 7CB (see Supplementary Note 6). As the microscopic origin for the singular properties about the surface tension, the contribution of the orientational order parameter in the surface has been discussed in refs.47–50, while the origin remains to be solved. Anyway, we consider that the inversion of the director rotation in the I + Ch phase might be due to the inversion of the Marangoni convection. To validate this notion, we measured the flow field inside the droplet, as described in the next section.

Flow field in the droplet under temperature gradient

The flow field was measured using the fluorescence photobleaching method, as shown in Fig. 7. The use of a confocal microscope enabled us to three-dimensionally analyse the field. Assuming that the flow was circularly symmetric, we estimated the radial flow velocity vr for the measurement results obtained in each focal plane, as shown in Fig. 8 (see the Methods section and Supplementary Note 7). When T0 was relatively high in the I + Ch phase, an inward flow was observed near the high-temperature side of the substrate, whereas an outwards flow was observed near the low-temperature side, as shown in Figs. 7c, d, o, p, 8c, and d. These observations indicated the existence of convection inside the droplet, as shown in Fig. 7s and t. In contrast, when T0 was relatively low in the I + Ch phase, the measurement result near the high-temperature side of the substrate showed the existence of an outwards flow (Figs. 7b and 8b), in contrast to the inward flow observed at a higher T0 (Figs. 7c, d, 8c, and d). This inversion indicates that the direction of the convection was inverted depending on T0, and this dependence is consistent with the expectation of the direction of the Marangoni convection associated with director rotation in the previous section. Therefore, the measurement results obtained in this study support the notion that the director rotation in the I + Ch phase is driven by convection. Owing to the coupling between the director and the flow fields, the director rotation is induced as described in refs31,33.Figure 7 Measurement results of the flow field. The measurement was performed in the focal planes (I)–(IV) shown in (q). The planes (I) and (II) are 15 and 35μm far from the high-temperature side substrate respectively, and (III) and (IV) are 35 and 15 μm from the low-temperature side. For (I) and (II), the light was incident from the high-temperature side, and for (III) and (IV), the light was from the low-temperature side. The measurement was performed under four different temperatures T0 of 34.0, 38.0, 42.0 and 46.0 °C, and the temperature difference ΔT was set to be 12K. The droplet was in the Type-A Ch state when T0 = 34.0 °C, and I + Ch phase for otherwise, as illustrated in (q–t). In (a–p), the distribution of the flow velocity in each focal plane and T0 is drawn together with the photo-bleached pattern, which is normalised by the fluorescence image before the bleaching. The black bar in (a) indicates 100 μm.

Figure 8 Temperature dependence of radial flow velocity vr. The velocity was estimated by the measurement results shown in Fig. 7. The measurement was performed in the four focal planes of (I)–(IV) illustrated in Fig. 7q under four different temperatures T0 of 34.0, 38.0, 42.0 and 46.0 °C, as shown in (a–d). vr is plotted with respect to r/R, where r is radial coordinate, and R is the droplet radius. The positive value of vr indicates the outward flow, and negative vr indicates inward.

In the I + Ch phase with a low T0, flow was hardly observed near the low-temperature side of the substrate (Figs. 7n and 8b). This indicates that convection was localised near the high-temperature side of the substrate, as shown in Fig. 7r, in contrast to the convection induced at a higher T0 (Fig. 7s and t). This difference in the location of convection results in a large difference in the rotational speed of the director field of the I + Ch phase. As mentioned in the previous section, the speed decreased by one order of magnitude when T0 decreased in the phase (Fig. 6a and b). This drastic decrease can be attributed to the localisation of the Marangoni convection near the high-temperature region in the droplet: owing to the localisation, the flow generates the rotational torque only in the narrow region near the I-Ch interface when T0 is relatively low (see Fig. 7r), and the rotation is strongly suppressed by the dissipation produced in the region without flows.

Under the assumption that the rotation is induced by the Marangoni flow, we can give a rough estimation for the rotational speed. Considering the advection effect of the director field according to refs.31,33, the angular velocity of the rotation can be estimated as ω∼2πV2πVP0P0, where V is a characteristic flow speed and P0 is the helical pitch length. Assuming that P0 is ~ 50 µm (see Method section) and V is the order of 10–1–100 µm/s in the I + Ch phase from Fig. 8, ω can be estimated to be ~ 101–102 mrad/s. This is overlapping with the measured rotational speed in the I + Ch phase as shown in Fig. 6a and b, whereas ranges slightly higher than the measurement results. Since the factors inhibiting the rotation, such as viscous dissipation and the deformation of the director field, were neglected in the estimation, it necessarily results in an overestimation. Thus, it is sufficiently reasonable to consider that the rotation is mainly induced by the flow, also from the quantitative discussion about the rotational speed.

In the Ch phase, the flow was barely observed anywhere in the droplet, as shown in Figs. 7 a, e, i, m, and 8a. This is consistent with the rotational speed: the speed in the Ch phase was one or two orders of magnitude lower than that in the I + Ch phase (Fig. 6a and b). It is difficult to determine the mechanism of rotation in the Ch phase based on the results obtained in this study. One potential explanation for the mechanism is the direct coupling effect predicted by Leslie23,36. However, Marangoni convection remains a plausible candidate as well; the weak flow, which was not detected by our measurements owing to insufficient resolution, might result in slow rotation. As seen here, in the rotating Ch LC droplet under the temperature gradient, it is generally difficult to discuss in detail about the contribution of the Leslie’s thermomechanical coupling effect, separately from the possibility that the rotation is induced by other factors, such as the flow. How the existence or non-existence of the coupling effect should be verified is still a problem to be solved.

As described in this section, the flow property drastically changed, depending on the temperature. Unfortunately, we cannot clarify the reason why the drastic change was induced in the present situation. For the clarification, detailed measurement and discussion about the distribution of the temperature and the surface tension in the droplet and the energy dissipation due to the anisotropic viscosity of LC would be necessary.

Discussion

Based on the experimental results, we discuss the deformation process of the director fields induced by flow. First, we consider the temperature dependence of Rc, indicating the distance of the point defect from the droplet centre. Since RcRcRR shows a similar behaviour in both N and Ch droplets, except for the hysteresis property (Fig. 3a and b), we should consider that there is a common mechanism about the migration of the point defect apart from the droplet centre. For simplicity, here we consider the situation in N LC droplet. When T0 was relatively high in the coexistence phase, convection occurred in the entire droplet, as shown in Fig. 7s and t. In the similar way with refs.31,38,39,44,45, where the convection inside the droplet is discussed, we designed the trial functions for the flow velocity v = (vr, vϕ, vz) as shown below:1 vr=16Vsr1-r2R2z(h-z)(h-2z)h4vϕ=0vz=-16Vs1-2r2R2z2(h-z)2h4,

where we used cylindrical coordinates whose origin is set at the centre of the bottom plane of the droplet, as the present droplet has rotational symmetry. R and h are the droplet radius and height, respectively. Vs is the characteristic velocity; in the central plane of the z=hh22, vz is − Vs and Vs at the centre (r = 0) and circumference (r = R), respectively. Equation (1) satisfies the incompressible (∇·v = 0) and non-slip conditions at the substrates (vr = 0 for z = 0, h)43. The influx and outflux are zero at the droplet surface (vr = 0 for r = R, vz = 0 for z = 0, h). vϕ is assumed to be zero, because no rotational flow is considered to be induced in the N LC droplet. The velocity field described by Eq. (1) is depicted in Fig. 9a and d.Figure 9 Analysis of flow and director fields in I + N phase. (a) and (d) represent the flow field depicted based on Eq. (1); (b), (c), (e) and (f) represent the director fields based on Eqs. (2) and (3). (a–c) show the fields at the plane including the central axis of the droplet, and (d–f) at the plane indicated by the broken line in (a–c). (g) and (i) are the light transmission intensity profiles calculated by the Jones matrix method under the assumption of Eqs. (2) and (3). P and A in (h) indicate polariser and analyser, respectively. RcRcRR is set to be 0 in (b), (e) and (g), and to be 0.7 in (c), (f) and (i). (h) and (j) are POM images obtained in the I + N phase.

The director should align along the z-axis at the substrate owing to the strong homeotropic anchoring condition (see Methods). The side surface of the air interface also exhibits homeotropic anchoring5, whereas the interface between the I and LC phases is planar51,52. Considering these boundary conditions, radial alignment of the director should be realised in the I + N phase, as shown in Fig. 9b and e. The trial functions for similar situations have already been reported in refs.38,39. Based on these reports, we assumed the director field in the N droplet to be2 nr′=R′+Rccosϕ′Rsinθnnϕ′=-Rcsinϕ′Rsinθn,nz=cosθn

where we used the cylindrical coordinates of (r′, ϕ′, z). The parameters θn and R′ are described as3 θn=π2zhLCr′R′∑k=021-r′R′kR′=-Rccosϕ′+R2-Rc2sin2ϕ′,

where hLC is the thickness of the LC phase in the droplet (see Fig. 9c) and the director field described by Eqs. (2) and (3) is depicted in Fig. 9b, c, e, and f. Equation (2) has a singular axis at r′=0, where the director is along the z-axis, and the coordinate origin is set to the point where the singular axis crosses the bottom plane of the droplet. Strictly speaking, the field does not include defects with a discontinuous director. However, the director varies drastically near the I-N interface in the singular axis; we can consider that the point defect exists substantially at (r, z) = (0, hLC). Under this consideration, the parameter Rc in Eqs. (2) and (3) indicates the distance of the defect from the central axis of the droplet, which is consistent with the definition of Rc described in the previous section (inset of Fig. 3a). Using the Jones matrix method, we numerically calculated the light transmission intensity profiles under crossed polarisers, as shown in Fig. 9g and i. In both Rc = 0 and Rc > 0, the profiles were well consistent with the POM images (Fig. 9g–j).

The relation between hLC and T0 can be estimated as follows:4 h-hLCh∼T0-TNIΔT+12,

where TNI is the I-N transition temperature. In Fig. 3a and b, Rc shows non-zero value when T0 ranges from 41 to 45 °C; in this range hLChLChh shows the values between 0.2 and 0.6. In this situation, the flow field in the N phase region is dominated by the outwards flow, as shown in Fig. 9d. Because of this field, a structure with a centred point defect is destabilised.

Using Eqs. (1) and (2), we deduced the time evolution of Rc. To adopt the Onsager variational principle, we calculated the Rayleighian ℜ composed of the dissipation function W and the time derivative of the free energy F (ℜ=W+F˙). In this case, W and F can be approximated as follows:5 W=γ1VsRhLCARcR+BRcR3R˙cR+γ1R2hLCXR˙cR2,F=K3hLCCRcR2+DRcR4,

where R˙c indicates the time derivative of Rc, and terms independent of Rc are neglected. The non-dimensional parameters A and B depend on hLChLChh, RRhh and γ2γ2γ1γ1, where γ1 and γ2 are the rotational viscosities. X is a positive constant and C and D depend on K1K1K3K3 and K2K2K3K3, where K1, K2 and K3 are elastic constants23,38,39,53,54 (for more details, see Supplementary Note 8).

Minimising ℜ with respect to R˙c according to the variational principle, we obtained the time evolution of Rc as follows:6 RcR=Rc0Rexp-tτ1+(γ1VsRB+4K3D)Rc02(γ1VsRA+2K3C)R21-exp-2tτ-12τ=2γ1R2Xγ1VsRA+2K3C,

where Rc0 is the initial value of Rc (Rc0 = Rc(t = 0)). In Eq. (6), as time t approaches infinity, Rc converges to zero or a non-zero value, depending on the sign of τ:7 RcR→0forτ>0RcR→-γ1VsRA+2K3Cγ1VsRB+4K3Dforτ<0.

Thus, the steady-state solution of Rc can exhibit a threshold behaviour. In addition, in both cases of Eq. (7), the time to reach steady state is characterised by the absolute value of τ defined in Eq. (6).

Using Eqs. (4), (6), and (7) and assuming that Vs is a constant of 1.0 μm/s, which is comparable to the measurement results of the flow field, we calculated the T0 dependence of the steady-state solution of RcRcRR and the absolute value of τ, as shown in Fig. 3c. RcRcRR shows a threshold behaviour, which is consistent with the behaviour observed in Fig. 3a and b. τ is of the order of 100 s and diverges near the threshold value. In Fig. 3a and b, temperature hysteresis of the order of 1 K is observed in RcRcRR. Considering that the cooling and heating rates are 0.2 K/min, we can estimate that the characteristic time for the hysteresis behaviour is ~ 300 s, which is comparable with τ of the characteristic time for the structural deformation. This indicates that the deformation progresses sufficiently slowly compared to the rate of change in the temperature; this slowness is considered to result in a temperature hysteresis.

The director field in the I + Ch phase is also deduced. The present situation is similar to the Ch LC confined to the sandwich cells with strong homeotropic anchoring as reported in refs.55–58. According to these literatures, when the helical pitch length P0 and the cell thickness d is comparable with each other, the director fields can show the two states called “translationally invariant configuration” (TIC) and “cholesteric fingers” (CF). The difference of these states is in the distribution of the helical axis: in TIC the axis is uniformly parallel to the cell depth direction, while in CF the axis also distributes in the plane parallel to the substrates in addition to the depth direction. Similar to this situation, in our system P0 and hLC are comparable with each other (see method section). Since the comma-like texture in relatively low T0 (Fig. 2g, h, m etc.) is similar to the pattern observed in CF, the helical axis parallel to the substrates is considered to be formed in this state. Here, we focus on the director field in relatively high T0, and discuss the field with a point defect as shown in Fig. 10. In this situation, a twisted cross-texture is often observed (e. g. Figure 2b and d), indicating that the director twists. We considered that this twist is attributed to the helical structure along the z-axis as similar to TIC, and that it gradually unwinds as r′ increases owing to homeotropic anchoring at the air interface (see Fig. 10a–f). Double- or single-loop defects were observed, in addition to helical structures (Fig. 4). For an intuitive understanding of the appearance of these defects, we considered the formation of cylindrical wall defects, which would actually deform into loops23,59. The wall resulted from the gap of the twist deformation along the z-axis; the twisting strength changed discontinuously at the wall (for more details, see Supplementary Note 9). Assuming that the director included double and single wall defects, we numerically calculated the light transmission intensity profiles under crossed polarisers, as shown in Fig. 10g and o, respectively. These profiles are well consistent with the POM images of Fig. 10h and p immediately after the double- and single-loop defects were formed, respectively. We consider that a field similar to that shown in Fig. 10 was formed in the I + Ch phase.Figure 10 Analysis of director field in I + Ch phase. (a–f) represent the director field with double wall defects, and (i–n) represent the field with a single wall. (e), (f), (m) and (n) show the fields at the plane including the central axis of the droplet. (a), (b), (i) and (j) show the director fields at the plane indicated by the orange broken line (I) in (e) and (m), and (c), (d), (k) and (l) at the plane of the purple broken line (II). (a), (c), and (e) are depicted based on Eqs. (3), (S14), (S15a) and (S15b) in Supplementary Note 9, and (i), (k) and (m) are based on Eqs. (3), (S14), (S17a) and (S17b). In the schematics of (b), (d), (f), (j), (l) and (n), the colour gradation of the ellipsoids from red to blue represents the director tilt from up to down against the plane. The green broken lines/circles in the schematics indicate the position of the wall defects. Actually, the walls assumed in these director fields are considered to transform into the loop defects. (g) and (o) are the light transmission intensity profiles calculated by the Jones matrix method; (g) is based on equations (3), (S14), (S15a) and (S15b), and (o) is on equations (3), (S14), (S17a) and (S17b). P and A in (g) and (o) indicate polariser and analyser, respectively. (h) and (o) are POM images obtained just after the loop defects are formed; double and single loops exist in (h) and (p), respectively.

The stability of the defects is considered to depend on the parameter RcRcRR. When RcRcRR is small, the double loop defects are formed (Fig. 4a–c); on the other hand, as RcRcRR increases the defects are destabilised and transformed into the single loop (Fig. 4d–f). In addition, RcRcRR shows the temperature hysteresis as shown in Fig. 3b. Owing to this hysteresis and the RcRcRR-dependent stability property of the defects, the topologically different structures would be formed, depending on the temperature change process (Figs. 2 and 4).

In summary, we observed the formation of a director field in disk-shaped N and Ch droplets sandwiched by substrates under a temperature gradient. In both cases, during the coexistence phase, the field was largely deformed, showing hysteresis (Figs. 1, 2, 3). In the I + Ch phase, the deformation was accompanied by the rotation of the director field owing to the temperature gradient (Lehmann rotation). In addition, depending on the temperature change process, the formation and reconnection of loop defects were observed (Fig. 4). Owing to this topological variation, the final structures obtained in the Ch phase depend on the process.

Three-dimensional measurements of the flow field revealed the existence of the Marangoni convection. Convection was strongly induced in the entire droplet when the temperature was relatively high in the I + Ch phase, while it was suppressed and localised near the high-temperature region as the temperature decreased; in the Ch phase, the flow was hardly observed. We should consider that the Lehmann rotation in the I + Ch phase was mainly induced by the convection. This notion is consistent with the fact that the speed of the rotation in the I + Ch phase was higher by 1–2 orders of magnitude than the speed in the Ch phase, and that both direction of the convection and the rotation switched depending on the temperature (Figs. 5, 6, 7, 8). The appearance of the Marangoni convection also results in structural deformation with thermal hysteresis. By estimating the characteristic time of deformation and comparing it with the time of hysteresis, we found that both times exhibited the same order of magnitude. This indicates that the hysteresis behaviour was strongly associated with that the structural deformation was induced enough slowly, compared with the change rate of the temperature. This hysteresis property is considered to result in the topological variation of loop defects depending on the temperature change process.

As mentioned in the Introduction, the morphogenesis of living organisms proceeds in the presence of a potential gradient. To mimic this aspect, in this study, we focused on the structure formation of a Ch LC droplet under a temperature gradient and found that several nontrivial phenomena were induced. A remarkable fact is the simplicity of our system, in which only a typical Ch LC was used as the sample material. This suggests that a special gimmick would not be required to induce the phenomena: we expect that the Lehmann rotation and the topological change induced by flow are not specific phenomena induced in specific situations, but can be observed generally in the morphogenesis process of living organisms and soft matter systems, including the LC. Deoxyribonucleic acids and proteins, which strongly contributes to the morphogenesis, are chiral materials1. Since the morphogenesis is accompanied by the potential gradient, the rotation might be induced in these materials during the morphogenesis, as seen in the Lehmann rotation of Ch LC under the temperature gradient. In addition to the gradient, it has been reported that the flow is induced and plays significant roles in the morphogenesis60–62. Since it eventually needs the topological change as seen in segmentation and organogenesis, we expect that the flow might result also in the topological change of the living organisms, as observed in the LC droplet in this study.

Methods

Sample preparation

A rod-like molecule, 7CB (LCC Co., Ltd.), was used as the N LC. By adding the chiral dopant (S)-2-octyl 4-[4-(hexyloxy)benzoyloxy] benzoate (Tokyo Chemical Industry Co., Ltd.) to 7CB, we prepared Ch LC with a dopant concentration of 0.16 wt%. The phase sequences of the LC samples are N-42°C-I and Ch-42°C-I. The helical pitch length of the Ch LC sample was measured to be P0 ~ 50μm at 41 °C using the Cano-wedge method23. For fluorescence microscopy in the flow field measurement, the fluorescent dye poly[tris(2,5-bis(hexyloxy)-1,4-phenylenevinylene)-alt-(1,3-phenylenevinylene)] (Sigma-Aldrich Co., LLC) was added to the Ch LC sample.

To prepare sandwiched N LC and Ch LC droplets, we prepared cover glass substrates coated with CYTOP (Asahi Glass Co., Ltd.) for homeotropic alignment. By sandwiching the N LC and Ch LC samples with the substrates, we made the droplets whose side surfaces were exposed to air (see Figs. 1 and 2(v)). The distance between the glass substrates was maintained at 100 μm using polyimide film spacers, and droplets with diameters of 250–300 μm were used in this study (for more detail, see Supplementary Note 2).

Temperature control and polarised microscopy

A temperature gradient was applied to the sample in a homemade furnace, as described in the literature18. The furnace consisted of a homemade hot stage and a commercial hot/cool stage (Tokai Hit Co., Ltd.). These stages respectively control the temperatures of the upper and lower substrates, touching the droplet (see Figs. 1 or 2(v)). In this study, the positive direction of the z axis was defined by the direction of the temperature gradient. The droplets were observed by polarised optical microscopy (POM) using a commercial upright microscope (ECLIPSE LV100, Nikon Co.).

Surface tension measurement

Surface tension was measured using two methods. One is the pendant drop method, in which tension is determined by the size and shape of the pendant droplet46,63. The other one can be called ‘bubble method’, which is described in refs.21,64,65. We created a bubble using the LC samples and measured the pressure difference between the inside and outside of the bubble. Based on the Young–Laplace equation, we obtained the surface tension using the relationship between the pressure difference and radius of the bubble.

Flow-field measurement by fluorescence photobleaching method

The flow field was measured using the previously described fluorescence photobleaching method21,27,33. The fluorescent LED illumination system D-LEDI (Nikon Co., Ltd.) was used as the light source for photobleaching. First, the sample was bleached into a striped pattern using strong light illumination through a photomask with a periodic array of slits. The flow distribution was obtained from the time evolution of the fluorescence images after photobleaching. Images were obtained using a confocal microscope system constructed with the MAICO MEMS confocal unit (Hamamatsu Photonics Co., Ltd.) and an inverted microscope Ti2-U (Nikon Co., Ltd.). The flow measurement was performed in each focal plane, as shown in Fig. 7, to realise a three-dimensional measurement of the flow field in the droplet. In this study, by defining a striped pattern aligned along the y-axis, we obtained the distribution of the flow velocity component in the x-direction, vx.

Setting the coordinate origin as the droplet centre, we considered that the observed 2D flow was a linear combination of radial and rotational flows. By defining the velocity components of the radial and azimuthal flows as vr and vϕ, respectively, we estimated them from the distribution of vx obtained from the measurement. For simplicity, assuming that vr and vϕ are independent of the azimuthal coordinates ϕ, we can describe vx as8 vx(r,ϕ)=vr(r)cosϕ-vϕ(r)sinϕ.

From the ϕ dependence on vx, we estimated vr(r), as shown in Fig. 8 (for more details, see Supplementary Note 7).

Analytical calculation and optical simulation

Analytical calculations were performed using the commercial software Mathematica (Wolfram Research, Inc.), and three-dimensional drawings of the director and flow fields based on the trial functions were created using the freely available POV-Ray software. An optical simulation based on trial functions of the director field was performed using the Jones matrix method66 (for more details, see ref.18). The ordinary refractive index of the sample was assumed to be 1.567. For the birefringence index Δn, we deduced Δn of 7CB by measuring the transmission spectrum under crossed polarisers. The measurement was performed by a commercial spectroscope of BIM-6001A-06–Miniature Spectrometer (Brolight Technology Co., Ltd.), and the temperature was set to be 41 °C.

Supplementary Information

Supplementary Information 1.

Supplementary Video 1.

Supplementary Video 2.

Supplementary Video 3.

Supplementary Video 4.

Supplementary Video 5.

Supplementary Video 6.

Supplementary Video 7.

Supplementary Video 8.

Supplementary Video 9.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-58054-2.

Acknowledgements

This work was supported by JSPS KAKENHI Grant Number 22H01191.

Author contributions

J.Y. designed the study, and the experiments were performed by J.Y. and Y.I. J.Y. and K.F. analysed the data and prepared the manuscript. All the authors discussed the results and approved the final manuscript.

Data availability

All data that support the findings in this study are available in the article and in Supplementary Information. Additional information is available from the corresponding author upon request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Bard J Morphogenesis: The Cellular and Molecular Processes of Developmental Anatomy 1990 Cambridge University Press
Bard, J. Morphogenesis: The Cellular and Molecular Processes of Developmental Anatomy (Cambridge University Press, 1990).
2. Truman JW The evolution of Insect metamorphosis Curr. Biol. 2019 29 R1252 1268 10.1016/j.cub.2019.10.009 31794762
Truman, J. W. The evolution of Insect metamorphosis. Curr. Biol. 29, R1252-1268 (2019).31794762 10.1016/j.cub.2019.10.009
3. Candau S Roy PL Debeauvais F Magnetic field effects in nematic and cholesteric droplets suspended in an isotropic liquid Mol. Cryst. Liq. Cryst. 1973 23 283 297 10.1080/15421407308083378
Candau, S., Roy, P. L. & Debeauvais, F. Magnetic field effects in nematic and cholesteric droplets suspended in an isotropic liquid. Mol. Cryst. Liq. Cryst. 23, 283–297 (1973).10.1080/15421407308083378
4. Drzaic PA New director alignment for droplets of nematic liquid crystal with low bend-to-splay ratio Mol. Cryst. Liq. Cryst. 1988 154 289 306
Drzaic, P. A. New director alignment for droplets of nematic liquid crystal with low bend-to-splay ratio. Mol. Cryst. Liq. Cryst. 154, 289–306 (1988).
5. Gupta VK Abbott NL Using droplets of nematic liquid crystal to probe the microscopic and mesoscopic structure of organic surfaces Langmuir 1999 15 7213 7223 10.1021/la981780g
Gupta, V. K. & Abbott, N. L. Using droplets of nematic liquid crystal to probe the microscopic and mesoscopic structure of organic surfaces. Langmuir 15, 7213–7223 (1999).10.1021/la981780g
6. Khullar S Zhou C Feng JJ Dynamic evolution of topological defects around drops and bubbles rising in a nematic liquid crystal Phys. Rev. Lett. 2007 99 237802 10.1103/PhysRevLett.99.237802 18233413
Khullar, S., Zhou, C. & Feng, J. J. Dynamic evolution of topological defects around drops and bubbles rising in a nematic liquid crystal. Phys. Rev. Lett. 99, 237802 (2007).18233413 10.1103/PhysRevLett.99.237802
7. Pairam E Stable nematic droplets with handles Proc. Natl. Acad. Sci. 2013 110 9295 9300 10.1073/pnas.1221380110 23690570
Pairam, E. et al. Stable nematic droplets with handles. Proc. Natl. Acad. Sci. 110, 9295–9300 (2013).23690570 10.1073/pnas.1221380110
8. Jiang J Yang D-K Bipolar to toroidal configuration transition in liquid crystal droplets Liq. Cryst. 2018 45 102 111 10.1080/02678292.2017.1301582
Jiang, J. & Yang, D.-K. Bipolar to toroidal configuration transition in liquid crystal droplets. Liq. Cryst. 45, 102–111 (2018).10.1080/02678292.2017.1301582
9. Yoshioka J Spherical-cap droplets of a photo-responsive bent liquid crystal dimer Soft Matter 2019 15 989 998 10.1039/C8SM01751D 30657150
Yoshioka, J. et al. Spherical-cap droplets of a photo-responsive bent liquid crystal dimer. Soft Matter 15, 989–998 (2019).30657150 10.1039/C8SM01751D
10. Peddireddy K Self-shaping liquid crystal droplets by balancing bulk elasticity and interfacial tension Proc. Natl. Acad. Sci. 2021 118 174118 10.1073/pnas.2011174118
Peddireddy, K. et al. Self-shaping liquid crystal droplets by balancing bulk elasticity and interfacial tension. Proc. Natl. Acad. Sci. 118, 174118 (2021).10.1073/pnas.2011174118
11. Xu F Crooker PP Chiral nematic droplets with parallel surface anchoring Phys. Rev. E 1997 56 6853 6860 10.1103/PhysRevE.56.6853
Xu, F. & Crooker, P. P. Chiral nematic droplets with parallel surface anchoring. Phys. Rev. E 56, 6853–6860 (1997).10.1103/PhysRevE.56.6853
12. Seč D Porenta T Ravnik M Žumer S Geometrical frustration of chiral ordering in cholesteric droplets Soft Matter 2012 8 11982 11988 10.1039/c2sm27048j
Seč, D., Porenta, T., Ravnik, M. & Žumer, S. Geometrical frustration of chiral ordering in cholesteric droplets. Soft Matter 8, 11982–11988 (2012).10.1039/c2sm27048j
13. Orlova T Aßhoff ST Yamaguchi T Katsons N Brasselet E Creation and manipulation of topological states in chiral nematic microspheres Nat. Commun. 2015 6 7603 10.1038/ncomms8603 26145716
Orlova, T., Aßhoff, S. T., Yamaguchi, T., Katsons, N. & Brasselet, E. Creation and manipulation of topological states in chiral nematic microspheres. Nat. Commun. 6, 7603 (2015).26145716 10.1038/ncomms8603
14. Yoshioka J Ito F Tabe Y Stability of a double twisted structure in spherical cholesteric droplets Soft Matter 2016 12 2400 2407 10.1039/C5SM02838H 26806762
Yoshioka, J., Ito, F. & Tabe, Y. Stability of a double twisted structure in spherical cholesteric droplets. Soft Matter 12, 2400–2407 (2016).26806762 10.1039/C5SM02838H
15. Posnjak G Čoper S Muševič I Hidden topological constellations and polyvalent charges in chiral nematic droplets Nat. Commun. 2017 8 14594 10.1038/ncomms14594 28220770
Posnjak, G., Čoper, S. & Muševič, I. Hidden topological constellations and polyvalent charges in chiral nematic droplets. Nat. Commun. 8, 14594 (2017).28220770 10.1038/ncomms14594
16. Krakhalev MN Bipolar configuration with twisted loop defect in chiral nematic droplets under homeotropic surface anchoring Sci. Rep. 2017 7 14582 10.1038/s41598-017-15049-6 29109533
Krakhalev, M. N. et al. Bipolar configuration with twisted loop defect in chiral nematic droplets under homeotropic surface anchoring. Sci. Rep. 7, 14582 (2017).29109533 10.1038/s41598-017-15049-6
17. Sleczkowski P Light-activated helical inversion in cholesteric liquid crystal microdroplets Proc. Natl. Acad. Sci. 2018 115 4334 4339 10.1073/pnas.1720742115 29626129
Sleczkowski, P. et al. Light-activated helical inversion in cholesteric liquid crystal microdroplets. Proc. Natl. Acad. Sci. 115, 4334–4339 (2018).29626129 10.1073/pnas.1720742115
18. Yoshioka J Araoka F Topology-dependent self-structure mediation and efficient energy conversion in heat-flux-driven rotors of cholesteric droplets Nat. Commun. 2018 9 432 10.1038/s41467-018-02910-z 29382841
Yoshioka, J. & Araoka, F. Topology-dependent self-structure mediation and efficient energy conversion in heat-flux-driven rotors of cholesteric droplets. Nat. Commun. 9, 432 (2018).29382841 10.1038/s41467-018-02910-z
19. Biagio RL Souza RT Evangelista LR Zola RS Frustrated structures and pattern formation after thermal quenches in cholesteric liquid crystal droplets J. Mater. Chem. C 2021 9 8623 8639 10.1039/D1TC02056K
Biagio, R. L., Souza, R. T., Evangelista, L. R. & Zola, R. S. Frustrated structures and pattern formation after thermal quenches in cholesteric liquid crystal droplets. J. Mater. Chem. C 9, 8623–8639 (2021).10.1039/D1TC02056K
20. Gardymova AP Krakhalev MN Rudyak VY Barbashov VA Zyryanov VY Polymer-dispersed cholesteric liquid crystal under homeotropic anchoring: Electrically induced structures with λ1/2-disclination Polymers 2022 14 1454 10.3390/polym14071454 35406327
Gardymova, A. P., Krakhalev, M. N., Rudyak, V. Y., Barbashov, V. A. & Zyryanov, V. Y. Polymer-dispersed cholesteric liquid crystal under homeotropic anchoring: Electrically induced structures with λ1/2-disclination. Polymers 14, 1454 (2022).35406327 10.3390/polym14071454
21. Yoshioka J Sakikawa T Ito Y Fukao K Marangoni convection driven by temperature gradient near an isotropic-nematic phase transition point Phys. Rev. E 2022 105 L012701 10.1103/PhysRevE.105.L012701 35193218
Yoshioka, J., Sakikawa, T., Ito, Y. & Fukao, K. Marangoni convection driven by temperature gradient near an isotropic-nematic phase transition point. Phys. Rev. E 105, L012701 (2022).35193218 10.1103/PhysRevE.105.L012701
22. Lehmann O Structure, system and magnetic behaviour of liquid crystals and their miscibility with the solid ones Ann. Phys. 1900 2 649 705 10.1002/andp.19003070802
Lehmann, O. Structure, system and magnetic behaviour of liquid crystals and their miscibility with the solid ones. Ann. Phys. 2, 649–705 (1900).10.1002/andp.19003070802
23. de Gennes PG Prost J The Physics of Liquid Crystals 1993 2 Clarendon Press
de Gennes, P. G. & Prost, J. The Physics of Liquid Crystals 2nd edn. (Clarendon Press, 1993).
24. Madhusudana NV Pratibha R An experimental investigation of electromechanical coupling in cholesteric liquid crystals Liq. Cryst. 1989 5 1827 1840 10.1080/02678298908045691
Madhusudana, N. V. & Pratibha, R. An experimental investigation of electromechanical coupling in cholesteric liquid crystals. Liq. Cryst. 5, 1827–1840 (1989).10.1080/02678298908045691
25. Tabe Y Yokoyama H Coherent collective precession of molecular rotors with chiral propellers Nat. Mater. 2003 2 806 809 10.1038/nmat1017 14625540
Tabe, Y. & Yokoyama, H. Coherent collective precession of molecular rotors with chiral propellers. Nat. Mater. 2, 806–809 (2003).14625540 10.1038/nmat1017
26. Oswald P Dequidt P Measurement of the continuous Lehmann rotation of cholesteric droplets subjected to a temperature gradient Phys. Rev. Lett. 2008 100 217802 10.1103/PhysRevLett.100.217802 18518635
Oswald, P. & Dequidt, P. Measurement of the continuous Lehmann rotation of cholesteric droplets subjected to a temperature gradient. Phys. Rev. Lett. 100, 217802 (2008).18518635 10.1103/PhysRevLett.100.217802
27. Yoshioka J Director/barycentric rotation in cholesteric droplets under temperature gradient Soft Matter 2014 10 5869 5877 10.1039/C4SM00670D 24866557
Yoshioka, J. et al. Director/barycentric rotation in cholesteric droplets under temperature gradient. Soft Matter 10, 5869–5877 (2014).24866557 10.1039/C4SM00670D
28. Yamamoto T Kuroda M Sano M Three-dimensional analysis of thermo-mechanically rotating cholesteric liquid crystal droplets under a temperature gradient EPL 2015 109 46001 10.1209/0295-5075/109/46001
Yamamoto, T., Kuroda, M. & Sano, M. Three-dimensional analysis of thermo-mechanically rotating cholesteric liquid crystal droplets under a temperature gradient. EPL 109, 46001 (2015).10.1209/0295-5075/109/46001
29. Ignés-Mullol J Poy G Oswald P Continuous rotation of achiral nematic liquid crystal droplets driven by heat flux Phys. Rev. Lett. 2016 117 057801 10.1103/PhysRevLett.117.057801 27517793
Ignés-Mullol, J., Poy, G. & Oswald, P. Continuous rotation of achiral nematic liquid crystal droplets driven by heat flux. Phys. Rev. Lett. 117, 057801 (2016).27517793 10.1103/PhysRevLett.117.057801
30. Bono S Maruyama Y Tabe Y Formation and dynamics of the aggregates of cholesteric double-twist cylinders Soft Matter 2018 14 9798 10.1039/C8SM01565A 30398276
Bono, S., Maruyama, Y. & Tabe, Y. Formation and dynamics of the aggregates of cholesteric double-twist cylinders. Soft Matter 14, 9798 (2018).30398276 10.1039/C8SM01565A
31. Oswald P Ignés-Mullol J Dequidt A Lehmann rotation of cholesteric droplets driven by Marangoni convection Soft Matter 2019 15 2591 2604 10.1039/C8SM02574F 30816902
Oswald, P., Ignés-Mullol, J. & Dequidt, A. Lehmann rotation of cholesteric droplets driven by Marangoni convection. Soft Matter 15, 2591–2604 (2019).30816902 10.1039/C8SM02574F
32. Nishiyama K Bono S Maruyama Y Tabe Y Direct observation of rigid-body rotation of cholesteric droplets subjected to a temperature gradient J. Phys. Soc. Jpn. 2019 88 063601 10.7566/JPSJ.88.063601
Nishiyama, K., Bono, S., Maruyama, Y. & Tabe, Y. Direct observation of rigid-body rotation of cholesteric droplets subjected to a temperature gradient. J. Phys. Soc. Jpn. 88, 063601 (2019).10.7566/JPSJ.88.063601
33. Yoshioka J Araoka F Differential rotation in cholesteric pillars under a temperature gradient Sci. Rep. 2020 10 17226 10.1038/s41598-020-73024-0 33057019
Yoshioka, J. & Araoka, F. Differential rotation in cholesteric pillars under a temperature gradient. Sci. Rep. 10, 17226 (2020).33057019 10.1038/s41598-020-73024-0
34. Kiang-ia J Anomalous Lehmann rotation of achiral nematic liquid crystal droplets trapped under linearly polarized optical tweezers Molecules 2021 26 4108 10.3390/molecules26144108 34299382
Kiang-ia, J. et al. Anomalous Lehmann rotation of achiral nematic liquid crystal droplets trapped under linearly polarized optical tweezers. Molecules 26, 4108 (2021).34299382 10.3390/molecules26144108
35. Takano S Bono S Tabe Y Heat-flux-driven rotation of cholesteric droplets dispersed in glycerol J. Phys. Soc. Jpn. 2023 92 024601 10.7566/JPSJ.92.024601
Takano, S., Bono, S. & Tabe, Y. Heat-flux-driven rotation of cholesteric droplets dispersed in glycerol. J. Phys. Soc. Jpn. 92, 024601 (2023).10.7566/JPSJ.92.024601
36. Leslie FM Some thermal effects in cholesteric liquid crystals Proc. R. Soc. Lond. A 1968 307 359 372 10.1098/rspa.1968.0195
Leslie, F. M. Some thermal effects in cholesteric liquid crystals. Proc. R. Soc. Lond. A 307, 359–372 (1968).10.1098/rspa.1968.0195
37. Jiji LM Danish-Yaadi AH Heat Conduction 2024 4 Springer Nature
Jiji, L. M. & Danish-Yaadi, A. H. Heat Conduction 4th edn. (Springer Nature, 2024).
38. Yoshioka J Fukao Self-excited oscillation of the director field in cholesteric liquid crystalline droplets under a temperature gradient J. Phys. Condens. Matter 2020 32 325102 10.1088/1361-648X/ab83b1
Yoshioka, J. & Fukao,. Self-excited oscillation of the director field in cholesteric liquid crystalline droplets under a temperature gradient. J. Phys. Condens. Matter 32, 325102 (2020).10.1088/1361-648X/ab83b1
39. Yoshioka J Fukao K Horizontal transportation of a Maltese cross pattern in nematic liquid crystalline droplets under a temperature gradient Phys. Rev. E 2019 99 022702 10.1103/PhysRevE.99.022702 30934222
Yoshioka, J. & Fukao, K. Horizontal transportation of a Maltese cross pattern in nematic liquid crystalline droplets under a temperature gradient. Phys. Rev. E 99, 022702 (2019).30934222 10.1103/PhysRevE.99.022702
40. Siedler LTS Hyde AJ Pethrick RA Leslie FM Zvetkow twist viscosity measurements of some nematic liquid crystals Mol. Cryst. Liq. Cryst. 1983 90 255 10.1080/00268948308072454
Siedler, L. T. S., Hyde, A. J., Pethrick, R. A. & Leslie, F. M. Zvetkow twist viscosity measurements of some nematic liquid crystals. Mol. Cryst. Liq. Cryst. 90, 255 (1983).10.1080/00268948308072454
41. Liu PY Jamieson AM Twist viscosity of mixtures of low molar mass nematics Rheol. Acta 2000 39 532 10.1007/s003970000118
Liu, P. Y. & Jamieson, A. M. Twist viscosity of mixtures of low molar mass nematics. Rheol. Acta 39, 532 (2000).10.1007/s003970000118
42. Karat PP Madhusudana NV Elasticity and orientational order in some 4’-n-Alkyl-4-cyanobiphenyls: Part II Mol. Cryst. Liq. Cryst. 1977 40 239 10.1080/15421407708084487
Karat, P. P. & Madhusudana, N. V. Elasticity and orientational order in some 4’-n-Alkyl-4-cyanobiphenyls: Part II. Mol. Cryst. Liq. Cryst. 40, 239 (1977).10.1080/15421407708084487
43. Guyon E Hulin J-P Petit L Mitescu CD Physical Hydrodynamics 2015 2 Oxford University Press
Guyon, E., Hulin, J.-P., Petit, L. & Mitescu, C. D. Physical Hydrodynamics 2nd edn. (Oxford University Press, 2015).
44. Levan MD Motion of a droplet with a Newtonian interface J. Colloid Interface Sci. 1981 83 11 17 10.1016/0021-9797(81)90003-5
Levan, M. D. Motion of a droplet with a Newtonian interface. J. Colloid Interface Sci. 83, 11–17 (1981).10.1016/0021-9797(81)90003-5
45. Herminghaus S Interfacial mechanisms in active emulsions Soft Matter 2014 10 7008 7022 10.1039/C4SM00550C 24924906
Herminghaus, S. et al. Interfacial mechanisms in active emulsions. Soft Matter 10, 7008–7022 (2014).24924906 10.1039/C4SM00550C
46. Tintaru M Moldovan R Beica T Frunza S Surface tension of some liquid crystals in the cyanobiphenyl series Liq. Cryst. 2001 28 793 797 10.1080/02678290010025459
Tintaru, M., Moldovan, R., Beica, T. & Frunza, S. Surface tension of some liquid crystals in the cyanobiphenyl series. Liq. Cryst. 28, 793–797 (2001).10.1080/02678290010025459
47. Kasten H Strobl G Nematic wetting at the free surface of 4-cyano-4’-n-alkylbiphenyls J. Chem. Phys. 1995 103 6768 6774 10.1063/1.470355
Kasten, H. & Strobl, G. Nematic wetting at the free surface of 4-cyano-4’-n-alkylbiphenyls. J. Chem. Phys. 103, 6768–6774 (1995).10.1063/1.470355
48. Martínez-Ratón Y Velasco E Somoza AM Mederos L Sluckin TJ Theoretical study of the anomalous surface tension properties of liquid crystals J. Chem. Phys. 1998 108 2583 2593 10.1063/1.475643
Martínez-Ratón, Y., Velasco, E., Somoza, A. M., Mederos, L. & Sluckin, T. J. Theoretical study of the anomalous surface tension properties of liquid crystals. J. Chem. Phys. 108, 2583–2593 (1998).10.1063/1.475643
49. Shimizu RN Demarquette NR Study of the surface and interfacial tensions in systems containing a low molar mass liquid crystal Liq. Cryst. 2001 28 1855 1862 10.1080/02678290110082356
Shimizu, R. N. & Demarquette, N. R. Study of the surface and interfacial tensions in systems containing a low molar mass liquid crystal. Liq. Cryst. 28, 1855–1862 (2001).10.1080/02678290110082356
50. Gomes LS Demarquette NR Influence of temperature on surface tension of liquid crystals in the cyanobiphenyl and cyano-oxybiphenyl series Mol. Cryst. Liq. Cryst. 2005 437 181 194 10.1080/15421400590954669
Gomes, L. S. & Demarquette, N. R. Influence of temperature on surface tension of liquid crystals in the cyanobiphenyl and cyano-oxybiphenyl series. Mol. Cryst. Liq. Cryst. 437, 181–194 (2005).10.1080/15421400590954669
51. Faetti S Palleschi V Nematic-isotropic interface of some members of the homologous series of 4-cyano-4'-(n-alkyl) biphenyl liquid crystals Phys. Rev. A 1984 30 3241 10.1103/PhysRevA.30.3241
Faetti, S. & Palleschi, V. Nematic-isotropic interface of some members of the homologous series of 4-cyano-4’-(n-alkyl) biphenyl liquid crystals. Phys. Rev. A 30, 3241 (1984).10.1103/PhysRevA.30.3241
52. Yokoyama H Kobayashi S Kamei H Measurement of director orientation at the nematic-isotropic interface using a substrate-nucleated nematic film Mol. Cryst. Liq. Cryst. 1984 107 311 10.1080/00268948408070444
Yokoyama, H., Kobayashi, S. & Kamei, H. Measurement of director orientation at the nematic-isotropic interface using a substrate-nucleated nematic film. Mol. Cryst. Liq. Cryst. 107, 311 (1984).10.1080/00268948408070444
53. Doi M Onsager’s variational principle in soft matter J. Phys. Condens. Matter 2011 23 284118 10.1088/0953-8984/23/28/284118 21709334
Doi, M. Onsager’s variational principle in soft matter. J. Phys. Condens. Matter 23, 284118 (2011).21709334 10.1088/0953-8984/23/28/284118
54. Doi M Soft Matter Physics 2013 Oxford University Press
Doi, M. Soft Matter Physics (Oxford University Press, 2013).
55. Oswald P Baudry J Pirkl S Static and dynamic properties of cholesteric fingers in electric field Phys. Rep. 2000 337 67 96 10.1016/S0370-1573(00)00056-9
Oswald, P., Baudry, J. & Pirkl, S. Static and dynamic properties of cholesteric fingers in electric field. Phys. Rep. 337, 67–96 (2000).10.1016/S0370-1573(00)00056-9
56. Smalyukh II Electric-field-induced nematic-cholesteric transition and three-dimensional director structures in homeotropic cells Phys. Rev. E 2005 72 061707 10.1103/PhysRevE.72.061707
Smalyukh, I. I. et al. Electric-field-induced nematic-cholesteric transition and three-dimensional director structures in homeotropic cells. Phys. Rev. E 72, 061707 (2005).10.1103/PhysRevE.72.061707
57. Varanytsia A Topology-commanded optical properties of bistable electric field-induced torons in cholesteric bubble domains Sci. Rep. 2017 7 16149 10.1038/s41598-017-16241-4 29170409
Varanytsia, A. et al. Topology-commanded optical properties of bistable electric field-induced torons in cholesteric bubble domains. Sci. Rep. 7, 16149 (2017).29170409 10.1038/s41598-017-16241-4
58. Echeverría-Alar S Emergence of disordered branching patterns in confined chiral nematic liquid crystals Proc. Natl. Acad. Sci. 2023 120 e2221000120 10.1073/pnas.2221000120 37027428
Echeverría-Alar, S. et al. Emergence of disordered branching patterns in confined chiral nematic liquid crystals. Proc. Natl. Acad. Sci. 120, e2221000120 (2023).37027428 10.1073/pnas.2221000120
59. Panasyuk G Allender DW Model for the director and electric field in liquid crystal cells having twist walls or disclination lines J. Appl. Phys. 2002 91 9603 9612 10.1063/1.1477613
Panasyuk, G. & Allender, D. W. Model for the director and electric field in liquid crystal cells having twist walls or disclination lines. J. Appl. Phys. 91, 9603–9612 (2002).10.1063/1.1477613
60. Helm CE Fleury ME Zisch AH Boschetti F Swartz MA Synergy between interstitial flow and VEGF directs capillary morphogenesis in vitro through a gradient amplification mechanism Proc. Natl. Acad. Sci. 2005 102 15779 15784 10.1073/pnas.0503681102 16249343
Helm, C. E., Fleury, M. E., Zisch, A. H., Boschetti, F. & Swartz, M. A. Synergy between interstitial flow and VEGF directs capillary morphogenesis in vitro through a gradient amplification mechanism. Proc. Natl. Acad. Sci. 102, 15779–15784 (2005).16249343 10.1073/pnas.0503681102
61. Howard J Grill SW Bois JS Turing’s next steps: The mechanochemical basis of morphogenesis Nat. Rev. Mol. Cell Biol. 2011 12 392 398 10.1038/nrm3120 21602907
Howard, J., Grill, S. W. & Bois, J. S. Turing’s next steps: The mechanochemical basis of morphogenesis. Nat. Rev. Mol. Cell Biol. 12, 392–398 (2011).21602907 10.1038/nrm3120
62. Alt S Ganguly P Salbreux G Vertex models: From cell mechanics to tissue morphogenesis Philos. Trans. R. Soc. B 2017 372 20150520 10.1098/rstb.2015.0520
Alt, S., Ganguly, P. & Salbreux, G. Vertex models: From cell mechanics to tissue morphogenesis. Philos. Trans. R. Soc. B 372, 20150520 (2017).10.1098/rstb.2015.0520
63. Matsubara H Murase M Mori YH Nagashima A Measurement of the surface tensions and the interfacial tensions of n-pentane-water and R113-water systems Int. J. Thermophys. 1988 9 409 424 10.1007/BF00513080
Matsubara, H., Murase, M., Mori, Y. H. & Nagashima, A. Measurement of the surface tensions and the interfacial tensions of n-pentane-water and R113-water systems. Int. J. Thermophys. 9, 409–424 (1988).10.1007/BF00513080
64. Stannarius R Cramer C Self-supporting bubbles of thermotropic smectic liquid crystals Europhys. Lett. 1998 42 43 10.1209/epl/i1998-00543-x
Stannarius, R. & Cramer, C. Self-supporting bubbles of thermotropic smectic liquid crystals. Europhys. Lett. 42, 43 (1998).10.1209/epl/i1998-00543-x
65. Schüring H Thieme C Stannarius R Surface tensions of smectic liquid crystals Liq. Cryst. 2001 28 241 10.1080/02678290010006270
Schüring, H., Thieme, C. & Stannarius, R. Surface tensions of smectic liquid crystals. Liq. Cryst. 28, 241 (2001).10.1080/02678290010006270
66. Jones RC A new calculus for the treatment of optical systems I. Description and discussion of the calculus J. Opt. Soc. Am. 1941 31 488 493 10.1364/JOSA.31.000488
Jones, R. C. A new calculus for the treatment of optical systems I. Description and discussion of the calculus. J. Opt. Soc. Am. 31, 488–493 (1941).10.1364/JOSA.31.000488
67. Li J Gauzia S Wu S-T High temperature-gradient refractive index liquid crystals Opt. Express 2004 12 2002 2010 10.1364/OPEX.12.002002 19475035
Li, J., Gauzia, S. & Wu, S.-T. High temperature-gradient refractive index liquid crystals. Opt. Express 12, 2002–2010 (2004).19475035 10.1364/OPEX.12.002002
