1. Introduction
The discovery of the ferroelectric nematic phase [
1,
2,
3,
4,
5] was a groundbreaking event, since the combination of fluidity in ferroelectric fluid and its polar coupling to electric fields allow the observation of a whole new world of phenomena, which are rapidly becoming the focus of an increasing number of scientists, as demonstrated by the large number of peer review articles already published on the subject [
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16].
In this scenario, we recently performed experiments aimed at characterizing the behavior of sessile ferroelectric nematic liquid crystal droplets into contact with a ferroelectric solid substrate [
8]. We observed that upon entering the ferroelectric nematic phase, droplets experience an electromechanical instability that manifests itself through the sudden ejection of fluid jets, which branch into smaller streams and eventually form secondary tiny droplets [
8]. This behavior resembles the instability predicted by Lord Rayleigh in 1882, for charged conductive liquid droplets that are above the critical charge-to-volume ratio [
17]. Since ferroelectric nematic droplets are neutral, in our case this process occurs in the absence of free charges, but the required charging within the droplet arises from the intrinsic polarization of the ferroelectric liquid crystal via its contact with the ferroelectric substrate. The coupling between the polarization in the solid and fluid materials is mediated by the fringing field generated by the pyroelectric charging of the substrate. This polarization coupling induces the accumulation of surface charges on the droplet–air interface. As the droplet polarization grows by cooling the material, the local accumulation of polarization charges produces repulsive forces that overcome the surface tension. As this condition is met, the instability turns to an explosive runaway process, since the flow of the ferroelectric liquid crystal in the arising jet induces the orientational order of the nematic director along the jet direction. This transports polarization charge to the tip, thereby increasing the electrostatic repulsion.
The observed polarization-induced droplet instability crucially depends on the properties of the ferroelectric nematic phase and on the combination of polarization and fluidity unique to this system. Ferroelectric nematic liquid crystals offer for the first time the opportunity of studying the coupling between the polarization of a fluid and a solid ferroelectric, which understanding might provide the basis for novel electro-hydromechanical applications.
In this work, we aim at characterizing in detail the effects of the fringing field created by the ferroelectric substrate on the droplets’ instability. To this purpose, we did not work with sessile droplets but confined them between two ferroelectric substrates arranged in two different ways, to generate fringing fields with different profiles. Specifically, we realized sandwich cells where the substrates’ polarization vectors are parallel or antiparallel to each other, so to develop opposite or equal surface electrostatic charges at the interfaces with the liquid crystal droplet.
Noteworthy, conventional nematic liquid crystals in combination with ferroelectric substrates have been studied in several different configurations see for example [
18,
19]. The consequence of substrate charging on the liquid crystal’s average molecular orientation has been demonstrated; however, no effect comparable to those observed with ferroelectric nematic droplets has been reported.
2. Materials and Methods
The ferroelectric liquid crystal used in this work is 4-[(4-nitrophenoxy)carbonyl]phenyl2,4-dimethoxybenzoate (RM734). This compound was synthesized as described in [
1], and its structure and phase diagram have already been reported [
1]. The ferroelectric nematic phase (N
F) appears here through a weakly fist order phase transition upon cooling from the conventional higher temperature nematic (N) phase and is stable in the temperature range 133–80 °C [
1]. The value of the spontaneous polarization
P of RM734 exceeds 6 µC/cm
2 at the lowest temperature in the N
F phase; moreover,
P is locally collinear to the molecular director
n, defining the average orientation of the molecular axis, and is either parallel or antiparallel to it [
1].
As ferroelectric solid substrates, we used 900 µm thick z-cut undoped lithium niobate (LN) crystals (PI-Kem). Although the bulk spontaneous polarization of LN crystals along the [0001] z-axis is of the order of 70 µC/cm
2, because of very efficient compensation mechanisms at the z-cut surfaces, the surface charge at equilibrium is only of the order of about 10
−2 µC/cm
2 [
20]. When, however, temperature variations are induced in the crystal, the surface charge of LN can significantly increase thanks to the pyroelectric effect [
21,
22,
23], a transient phenomenon observable during and shortly after the temperature variation, and due to the slow free charge relaxation in LN. The pyroelectric coefficient of undoped LN is known to be of the order of 10
−4 C/m
2K at room temperature [
24], increases by one order of magnitude around 100 °C [
25]. Given the temperature used in our experiments, dictated by the RM734 phase diagram, we can thus expect an induced surface charge density of the order of 1 µC/cm
2, for temperature variations of a few degrees, ramped in a short time compared to the LN charge relaxation. To match the conditions of the previous work on sessile droplets [
8], LN crystals were used as bare substrates, with no coating applied.
To prepare the samples used in this work, a RM734 droplet was deposited on a LN substrate at a temperature T = 200 °C, which corresponds to the liquid crystal isotropic phase, and then covered by a second LN slab, previously heated at the same T. RM734 droplets were realized following two steps. First, a small amount of RM734 powder was deposited at room temperature on a clean glass slide and heated to 150 °C until it melted. To create the initial droplets, a cold stainless needle was dipped into the melt and retracted, causing the droplet on its tip to solidify upon contact with the surrounding air. To increase the size of the RM734 “pearl”, rapid (to avoid re-melting) successive dipping is performed. Then, the pearl was remelted into a droplet on the heated LN substrate. The two substrates were then stuck together by means of two 100 µm thick stripes of Kapton tape, which also define the cell thickness. The liquid crystal droplet confined between two solid substrates assumes the form of a capillary bridge (
Figure 1a). The entire cell was then transferred to a small, closed oven, suitable for the optical microscope, and cooled down to the N
F phase. In order to perform a systematic analysis of the instability events, we decreased the cell temperature by steps of 5 °C ramped in 60 s each, resulting in a cooling rate of 0.08 °C/s. Due to the structure of the oven used, the temperature variation is the same in the entire cell, which guarantees the same amount of pyroelectric charging on both the LN confining crystals.
The liquid crystal cells realized with two LN substrates are of two different kinds. Specifically, LN crystals were arranged so as to expose equally or oppositely charged surfaces at the interfaces with the liquid crystal, as sketched in
Figure 1a. The two configurations will be referred to as
np/
pn and
np/
np, where
n and
p stand for negative and positive, respectively. Noteworthy, the large thickness of the cells (100 µm) makes the droplets’ volume high enough to balance the additional friction due to confinement and the much lower thermal gradients on the LN substrates resulting from the closed arrangement.
The details of the electromechanical droplets’ instability have been analyzed by POM observations and videos recorded by a CCD camera operating at 25 frames per second.
3. Results and Discussion
Due to the finite size of LN crystals, their pyroelectric charging gives rise to an external fringing field. This field is a fraction
f of the internal field σ
LN/ε
0, with
f depending on the crystal shape and size, and of the order of
f ≈ 10
−3 in our experimental conditions [
8]. The fringing fields generated by the LN slabs combine in different ways in
np/
np and
np/
pn cells, which results in different profiles and values of the total field present within the region between the substrates. The total fringing field
Eff in the two cells, in the absence of the liquid crystalline fluid bridge, is reported in
Figure 1.
Figure 1a shows the normalized field E (intensity and lines) and its x- and z-components on the xz plane at y = 0 for the two configurations (top:
np/
np; bottom:
np/
pn). The
np/
pn arrangement gives rise to a total fringing field, which is mainly vertical and quite uniform along the cell thickness. The x-component is different from zero, although very weak, only close to the substrate’s corners. On the contrary, the
np/
pn configuration produces a total fringing field with a very weak vertical component, different from zero only at the corners, and an x-component that is uniform along the thickness, and increases from the center to the edges of the LN plates. The vertical and in-plane field components are reported in panels b (
np/
np) and c (
np/
pn) at three different positions along the cell thickness d: z = 0, z = d/2 and z = d. For the in-plane component
Et, the field lines are also shown. They are radially distributed on the xy planes. The absolute values of the field are also different in the two cells, being about 8.4 × 10
5 V/m for the
np/
np cells and 3.2 × 10
5 V/m for the
np/
pn ones (maximum values).
The specific profile of
Eff affects the features of the electromechanical instability of the N
F fluid bridge. In particular, (i) the number of instability events observed within the temperature window corresponding to the N
F phase, (ii) the number of ejected fluid jets for each of these events and (iii) the temperature at which these events start, depend on the specific substrates’ arrangement. This is shown in
Figure 2a, where the number N of observed fluid jets is reported as a function of the instability temperature T
i, defined as the temperature at which the instability event takes place. Different colors correspond to different kinds of cells. It is evident that the instability starts at higher temperatures in
np/
np cells, which thus exhibit the highest number of instability events. Moreover, these events are, on average, characterized by the highest number of ejected jets. Since a higher temperature corresponds to a lower value of the N
F polarization
P [
1], the results in
Figure 2a indicate that in
np/
np cells, the charging threshold is reached for lower values of
P than in the other kind of cells.
The
np/
np configuration, with its mainly vertical and uniform fringing field, is similar to the geometry used by Mathe et al. in [
26]. In their study, a ferroelectric fluid bridge was confined between two conductive glasses and an external electric field was applied perpendicular to the cell substrates. The instability caused by this field was interpreted as a kind of labyrinthine instability, already observed both in magnetic fluid exposed to a magnetic field orthogonal to the bounding plates and in dielectric fluids [
27,
28,
29]. In both cases, the presence of a gap separating the fluid from the magnetic poles or from the electrodes is necessary for the instability to take place. Such a gap allows the presence of a field component parallel to the plates. In our case, the gap is intrinsically present since LN is itself an insulator.
The presence of the ferroelectric fluid bridge changes the fringing field profile reported in
Figure 1. In
np/
np cells, this change consists in the appearance of an additional component of the field parallel to the bounding surfaces. It is known that the bulk polarization of the N
F liquid crystal spontaneously self organizes to minimize the internal and external electric fields. Generally,
P will end up parallel to the interfaces to avoid the accumulation of surface charge σ =
P.
u (where
u is the unit vector perpendicular to the surfaces). Additionally,
P will adopt bend deformations, which do not produce space charge, thus preventing nonzero
∇∙P as much as compatible with geometric constraints. In the presence of the fringing field, which in this geometry is mainly normal to the LN/N
F interface planes, the ferroelectric nematic becomes polarized. This happens through a small reorientation of
P by an angle such that it deposits polarization charge on the fluid bridge surfaces, canceling the internal field, a peculiarity of the N
F phase referred to as “fluid superscreening” [
8,
13]. This process leads to a mismatch between the field inside the bridge and the one outside. Indeed, in the first case, the potential difference is virtually negligible, being due to the field in the thin gap between the LN-charged surface and the liquid crystal interfacial layer; in the second case, the potential difference is determined by the field in air. This mismatch generates an additional in-plane component of the field, which can be arbitrarily large, depending on the thickness of the layer along which the potential difference changes from zero to the value dσ
LN/ε
0. The additional field component drives jet ejection from the charge accumulation sites. The result is the occurrence of several instability events, characterized by the ejection of a large number of jets from different portions of the capillary bridge at the interfaces with the bounding plates (
Figure 3a–d).
An additional feature of np/np cells is the temperature at which the instability starts, which is higher than that in np/pn cells, a behavior that we ascribe to the higher value of Eff components that characterize this specific substrates’ arrangement.
In
np/
pn cells, the fringing field has a lower absolute value and is mostly in plane, while the vertical component is weak and practically negligible at the center of the cell.
P and
Eff are thus both parallel to the LN/N
F interface planes, and the required in-plane component of the field is present from the beginning due to the specific LN substrate arrangement. In these conditions, the polarization charges that cancel the internal field are generated at the N
F/air lateral interfaces and are most probably due to twist distortions compatible with
P being parallel to the two bounding surfaces. Again, this creates an additional electric field that, in this case, is comparable in magnitude to the original value of the fringing field. The threshold charging that gives rise to the electromechanical instability is reached here at temperatures lower than those in the
np/
np cells, corresponding to higher values of
P, and on average produce the ejection of a lower number of fluid jets. We understand this phenomenon to be a result of the weaker fringing field that characterizes this configuration. Noteworthy, experiments on N
F sessile droplets deposited on glass surfaces with patterned electrodes showed that jet ejection preferentially occurs in regions where the fringing field is in the vertical direction [
26]. This is an additional indication that the in-plane electric field component, arising in these conditions due to the N
F superscreening, is higher than in any other configuration.
Once jets are formed in
np/
pn samples, they are, however, more disruptive compared to the other kind of cell. This is shown in
Figure 2b, where the variation of the average droplet radius after each instability event is reported as a function of T
i. This parameter exhibits a decreasing trend for both configurations, but such a decreasing trend is different in the two situations, being smooth for the
np/np cells and very steep for the
np/pn ones. In this latter case, the average radius decreases by more than 60% after the first instability events, indicating an extremely explosive and disruptive phenomenon.
We understand this behavior to be due to the radial shape of
Eff in the
np/
pn LN arrangement. Indeed, jets are polar fluid tubes carrying polarization charges on their tips, which thus keep on moving in the direction of the field, accelerated by the field itself. The polar nature of the ejected jets is well visualized in
Figure 4a, where jets parallel to the polarizer appear dark, while those at 45 deg. exhibit the maximum brightness, showing that the liquid crystal director
n is along the jet axis. Since the polarization vector is locally collinear to the optical axis
n,
Figure 4a indicates that
P is also along the jet’s axis, in the direction of flow. Noteworthy, POM analysis of the nematic director in conditions involving thick cells with LN substrates is cumbersome and made difficult by the formation of a thick fluid bridge and by the birefringence of the substrates. For this reason, the image in
Figure 4a is related to a RM734 N
F sessile droplet lying on a LN substrate. The results are, however, extendable to the double substrate configuration.
In the np/pn configuration, the in-plane field that accelerates the jet tips depends on the specific substrates’ arrangement and remains approximately unperturbed during jet motion. On the contrary, in the np/np case, the in-plane component of the field is generated by NF superscreening and depends on the position of the polarization charges. As such, it changes both value and direction as soon as jets protrude from the fluid bridge. In this case, there is not a “constant” radial field that moves the charged expelled fluid along a fixed direction, and jets do not show the collective disruptive motion observed in np/pn cells. The observed disruptive instability in np/pn cells might additionally be due to the acceleration experienced by the jets’ tips, which may generate large distortions on the fluid bridge portions close to the ejection sites. This, in turn, causes additional charge accumulation in a sort of self-sustained effect.
An example of the instability of
np/
pn cells is reported in
Figure 3e–h.
The evolution of the electromechanical instability of N
F fluid bridges is shown in more detail in
Figure 4b–d. In these images, the gap between the LN crystals is larger than in
Figure 3 (150 vs. 100 μm), which results in a more evident meniscus (thicker dark ring around the “droplet”).
Figure 4b shows the appearance of small spikes that protrude from several portions of the bridge rim at the onset of the instability event. This is followed by the formation of a single main jet (
Figure 4c) and by the ejection of a large portion of the total fluid mass (
Figure 4d). These events are accompanied by the decrease of the droplet diameter, which is quantified in
Figure 2b for both kinds of cells. This specific sequence of frames is related to the
np/
pn configuration. Noteworthy, large variations of the droplet texture can be observed by comparing the three images, which suggests that instability induces severe rearrangements of the N
F liquid crystal polarization. This reasonably gives rise to the additional accumulation of bound charges which, in turn, stimulates additional fluid ejection as a strategy to decrease electrostatic energy.
Noteworthy, measuring the position of the jet tips from the initial frames after the ejection, we noticed that the average jet speed for equal values of Ti is higher in np/pn cells compared to the others, in agreement with the notion that fluid motion is faster in this configuration.
In conclusion, we studied the behavior of ferroelectric liquid bridges confined between two solid ferroelectric substrates arranged in different configurations realized so that, once pyroelectrically charged, they generate fringing fields of different values and shapes. Our observations highlighted that the features of the liquid crystal instability are affected by the specific fringing field profile in a way dominated by the minimization of the electrostatic energy associated with the bulk polarization of the ferroelectric fluid.
Our results show that the electromechanical instability of ferroelectric droplets confined between two ferroelectric solid substrates can be controlled in terms of instability, temperature, number of ejected jets and violence by acting on the substrates’ arrangement. This might allow the realization of novel electro-hydrodynamic applications based on the electrostatic instability of polar liquids.