Diego
Díaz
a,
Diana
Garcia-Gonzalez
ab,
Pravash
Bista
a,
Stefan A. L.
Weber
ac,
Hans-Jürgen
Butt
a,
Amy
Stetten
*a and
Michael
Kappl
*a
aMax Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany. E-mail: kappl@mpip-mainz.mpg.de; stetten@mpip-mainz.mpg.de
bPhysics of Fluids Group, Max Planck Center Twente for Complex Fluid Dynamics, Department of Science and Technology, University of Twente, P.O. Box 217, Enschede 7500 AE, The Netherlands
cDepartment of Physics, Johannes Gutenberg University, Staudingerweg 10, 55128 Mainz, Germany
First published on 24th January 2022
When neutral water drops impact and rebound from superhydrophobic surfaces, they acquire a positive electrical charge. To measure the charge, we analyzed the trajectory of rebounding drops in an external electric field by high-speed video imaging. Although this charging phenomenon has been observed in the past, little is known about the controlling parameters for the amount of drop charging. Here we investigate the relative importance of five of these potential variables: impact speed, drop contact area, contact line retraction speed, drop size, and type of surface. We additionally apply our previously reported model for sliding drop electrification to the case of impacting drops, suggesting that the two cases contain the same charge separation mechanism at the contact line. Both our experimental results and our theoretical model indicate that maximum contact area is the dominant control parameter for charge separation.
It is known that impacting drops on hydrophobic surfaces may lead to charge separation and consequent electrification of the liquid and surface. A rebounding drop usually leaves the surface negatively charged, while the drop charges positively.22–25 The existence of an electrical double layer and the ions naturally present in water could explain the charge separation process. One example of charge separation occurs when drops slide down a hydrophobic tilted plane. After sliding, drops usually acquire a positive charge and deposit a surface charge of opposite sign on the solid.26–30 In a similar way, water drops that condense and coalesce on a superhydrophobic surface jump-off with a positive charge.31 Applications using such charging include designing surface charge gradients to transport drops on superamphiphobic surfaces.20
Although many such drop charging phenomena have been observed, charge separation processes that involve the movement of the three-phase contact line are not well understood due to difficulties in performing quantitative experiments. Furthermore, charge separation caused by flowing liquids involves processes far from equilibrium. As a result, it is difficult to derive a quantitative thermodynamic theory to describe them. Different fundamental effects may lead to net drop charging and it is challenging to recognize which ones contribute or dominate.
We will focus on the case of impacting drops that acquire some charge during contact with a surface. The charging of such drops has never been systematically quantified. The magnitude of drop charge after impact and the dependency on parameters like drop size, contact area, impact speed, and the retraction motion of the contact line is still unclear. Knowledge of such dependencies is a prerequisite for the understanding of the charging mechanisms and controlling the charge of impacting drops.
In this article, we address the issue of systematic measurement and modeling of impacting drop charge. We use an external electric field to reproducibly quantify the charge of rebounding drops on different superhydrophobic surfaces. When a water drop impacts these surfaces, it bounces off with a net charge and then deflects due to the electrostatic force. Our measurements allow us to isolate the influence of diverse drop parameters. We propose an extension of a previous theoretical model for electrification of sliding drops,30 through which we can quantitatively describe the charging process. We show that the maximum spreading contact area is the dominant parameter controlling the amount of charge acquired.
We characterized the surface topography of silicone nanofilaments and candle soot by scanning electron microscopy (Fig. 1). Candle soot consists of carbon particles of up to 40 nm, forming a loose fractal-like network (Fig. 1a). We coated the particles with a silica shell of 20 nm.32 The static contact angles of 4 μL drops on these surfaces is 162° ± 2. Silicone nanofilaments form a dense network on the glass substrate, showing a tangled, bent and hooked distribution, with diameters between 20 and 50 nm (Fig. 1b).
We carried out three types of experiments: (1) 4 μL drops of deionized water were released from heights between 0.4 and 5 cm onto all the surfaces mentioned above, (2) volumes between 3.5 and 13 μL onto SN-FDTS from 3 cm height, and (3) 6 μL drops released onto SN-PDMS brush and candle soot surfaces from 1.5 cm height. We recorded the impacting process using a high-speed camera in front of the sample (Photron FastCam Mini UX100, 5000 fps) to track the motion of the drop mass centre (Fig. S1, ESI†). Two copper plates were placed vertically over a Teflon plate to act as electrodes (Fig. 2a) and hence generate a horizontal electrical field E = V/d = 55 kV m−1, where V is the applied voltage and d the separation distance of the copper plates (18 mm). The sample surfaces were placed between the electrodes on the Teflon plate. The Teflon served as a non-conductive surface to isolate copper plates and sample. A grounded syringe pump generated neutral water drops. We shielded drop and needle with an aluminum tube to prevent polarization effects. We controlled the maximum spreading radius Rmax by the Weber number We = ρv02D0/γ, where ρ, v0, D0 and γ are the density, impact speed (0.2 m s−1 < v0 < 1 m s−1), diameter and surface tension of the drop. We will consider the impact speed in terms of We in the next sections. Sample surfaces were neutralized by an Ionizing Air Blower (Aerostat PC Ionizing Air Blower, Pennsylvania, USA) for 2 minutes before each drop. We observed complete rebounds with conservation of volume in the range of 1 < We < 13, and drop break up for We > 13. We can assure then that drops remain in the Cassie–Baxter state.33
The distribution of the electric field along several directions between our two metal plates was simulated using COMSOL Multiphysics 5.5 (Section III, ESI†). The simulations indicated that the electric field varies with height less than ± 5% around the average value (Fig. S4a, ESI†) in the range of the rebounding heights of our drops, which were always below the height of the plates. Within the range of the drop lateral movement (maximum of 2 mm) around the center, the lateral change in field (Fig. S4b and c, ESI†) is less than 1 kV m−1 (∼2% error), leading to an error of less than 10% in the charge calculation. Therefore, according these small and systematic deviations, we assume for our experiments an homogeneous electric field in all directions.
The homogeneity of the electric field implies that drops laterally move as , where Δx = |x2–x1|, with x1 = x(t1) as the lateral position of the drop mass centre at which the drop bounces from the surface and x2 = x(t2) is the lateral position where the drop impacts the surface again after bouncing (Fig. S1, ESI†). t = t2–t1 is the time interval between these events. Since the electrostatic force is FE = QE = max (Q is the drop charge and m is the drop mass), we calculated the drop charge as Q = 2ρVΔx/Et2. In our experiments, V is the initial and constant drop volume as we ensure that the drops rebound completely. For the case of drop break up, we calculated the charge of the primary drop subtracting the volume of the secondary ejected drops (Fig. S1g, ESI†). We only considered Weber numbers, where the volume of the main drop after break up is just slightly smaller (3.7–3.8 μL) than the original drop (4 μL). Only for higher Weber numbers, we observed also deflection of the secondary drops in the electric field, indicating partial charge transfer to the secondary drop as well. For the drop charge calculations, we did not measure directly the charge distribution in the drop. However, since the electric field is homogeneous, the exact charge distribution does not affect our measurements because the total force on the drop will equal to that on a point charge at the drop center.
The drag force can be neglected during the entire process and is defined by the Stokes' law FD = 6πrηv (r is the drop radius), η the viscosity of the air and v the drop velocity. For instance, a 4 μL drop impacting at We = 5, will experience a drag force of 1 nN when is laterally deflected by the electric field (see ESI,† Section IV) and an electrostatic force of FE = 9.81 μN. However, the drag force could become larger for the highly deformed drop shapes, where its velocity dependence would be in the worst case: , with CD the drag coefficient, A the cross section area and vx the lateral drop velocity. Since the lateral drop velocity is close to zero at the spreading/retraction phase and at the maximum deformation moments upon bouncing, the lateral drag force will remain negligible at this point. When drops move close to the maximum lateral speed measured in the experiments (0.03 m s−1, see Fig. S5, ESI†), then the lateral drag force reaches less than 16% of the electrostatic force in the worst case. Thus, the drag force does not affect significantly our experiments.
We performed control experiments assuring that drop charge is independent of the magnitude of the external electric field (Fig. 2e). This ruled out the possibility of induced electrification or dielectrophoretic effects. Thus, the acquired charge is the result of the interaction between the drop and the superhydrophobic surface. During a rebound without electric field, drops can exhibit a horizontal displacement of 0.03 ± 0.01 mm (We < 2). This displacement sets the limit of our detection method, which is equivalent to less than 3 pC.
Our detection technique of drop charge was calibrated using a similar method to Stetten et al.30 We detected the current signal of a rebounding drop touching an electrode probe prior impact and after the rebound from a superhydrophobic surface (Section II.I, ESI†). Simultaneously, we recorded the process by high-speed video imaging (Fig. S2a–d, ESI†). From the current signals (Fig. S2e and f, ESI†), we obtained values of drop charge very close to those obtained by our electric field method and not charge for the falling drop prior impact.
Ions, electrons, among other possibilities could be responsible for the positive charge of drops that interact with hydrophobic surfaces. Nevertheless, it has not been proven yet. We attribute the positive charge of rebounding drops to an accumulation of hydronium ions. When a hydrophobic surface interacts with water, it charges negatively most likely due to the absorption of hydroxyl ions.34–38 Accordingly, our superhydrophobic surfaces tend to absorb negative charges in form of hydroxyl ions from the impacting water drop. Hydronium ions form the complementary part of the electric double layer. When the drop rebounds positively charged, some of these ions accumulate in the drop.
Extracting Ṙ(t) from the retracting curves (Fig. 3a and Fig. S5, ESI†), we observed a speed fluctuating around a central value during the drop recoil (Fig. S6a, ESI†). Therefore, we define an average retraction speed ret = Rmax/tr as a characteristic speed of the retraction phase. Here, tr is the time of the entire retraction movement until drop detachment. This parameter is relevant because charge deposition is most likely the result of a non-equilibrium process. If we assume an electric double layer formed at the coating/liquid interface, its disruption and subsequent charge separation is most likely influenced by the speed of the retraction motion.
We found that the charge for 4 μL drops tends to increase with ret (Fig. 3b). The increasing
ret is a result of a further spreading as We is greater (Fig. S7, ESI†). Consequently, one might be tempted to conclude that the faster the contact line moves, the higher the charge. However, the opposite trend is seen when decreasing the drop volume (Fig. 3c). When drop volume decreases, drops spread at smaller Rmax and rebounds occur at shorter contact times40 and tr. The lower values of tr are sufficient to produce a faster retracting motion, so decreasing the volume increases
ret. And yet, decreasing drop volume leads to decreased drop charge. This means that there is no direct correlation between
ret and Q. Thus, we can exclude
ret as a dominant parameter in the charge separation process.
Similarly, we can use the same argument for an average spreading speed s = Rmax/ts, where ts is the time at which the drop reaches Rmax from the impact on the surface. Taking
s of the impact for the cases of Fig. 3b and c, we found that drop charge varies with
ret and
s following the same trend (Fig. S8, ESI†). Thus,
s is also excluded as a controlling parameter.
Drop charge increased linearly with the Weber number in the region where drop break up does not occur (Fig. 4a). We can write the relationship as Q = β(We–We0) with β ∈ [2,5]. The parameter β depends on the coating and drop volume, and We0 corresponds to the lowest Weber number to observe a complete rebound. We demonstrated this linear increase by varying drop speed at a constant volume of 4 μL. The linear trend is robust, as varying the drop volume while keeping constant its falling height reproduces the same behaviour (Fig. 4a, white triangles). Thus, with higher We, we find higher drop charge. The data in the shaded region of drop break up shows the existence of a saturation zone where the drop charge values tend to be unaltered with further increase of We. Besides, the difference in the saturation points of the surface coatings suggests that the maximum charge acquired by a drop depends on the surface chemistry of the impacted sample.
Higher We also brings with it a larger maximum contact area between the drop and surface, which could also be increased by increasing the drop volume. Then the maximum contact area arises as an implicit parameter influencing charging. Considering a symmetric spreading, we calculated the maximum contact area as Amax = πRmax2. We found that the drop charge increases with Amax (Fig. 4b) as well as with We and drop volume. However, note that for experiments with similar contact area and different Weber number and drop volume, the charge shows similar values (Fig. 4b, circular dashed lines). In particular, drops of 4 μL (We = 8.70, v0 = 0.75 m s−1) and 6 μL (We = 7.60, v0 = 0.55 m s−1) rebounding on SN-FDTS reach a maximum contact area of Amax ∼ 6 mm2, with charge values of ∼29.3 pC and ∼30.0 pC respectively. This suggests that We and drop volume are not the main control parameters of the charge separation process. In fact, Amax is a function of these two variables. Therefore, we can conclude that Amax is the governing parameter in the charging mechanism.
In our case, drops are impacting and then retracting radially off the surface. During the impacting process, the contact line moves in two phases: spreading and retraction. We consider the formation of an electrical double layer as the drop spreads on the surface. At the end of the spreading, the drop reaches the maximum contact diameter. This is comparable to a drop that wets a hydrophobic surface prior sliding, until it reaches its maximum width or contact diameter. The formation time of the electrical double layer is estimated to be 6 μs,30 which occurs during the wetting phase (2 ms, Fig. 3a). Accordingly, the retracting movement of an impacting drop should lead to charge separation by disruption of the electrical double layer in a similar way to a sliding contact line.
In this model, we assume that the charge is only deposited at the rear of the receding contact line, so we only need to consider the retracting motion. It could be, of course, that charges already deposit on the solid surface during the spreading, where the electrical double layer forms. However, since the drop is still in contact with the surface during spreading, counter charges are expected to exist close by and the sum of all charges in the electrical double layer is zero. So true charge separation is most likely not to occur until the drop retracts, which breaks the electrical equilibrium at the liquid/solid interface as the surface and liquid are no more in contact. As the drop retracts after reaching Rmax, it deposits charge along its outer radius. This causes that drops gain the opposite net charge like a charging capacitor. Since the drop has a limited electrical capacitance, charge separation will be highest during the initial phase of the retraction and will become more and more hindered with increasing drop charge. Therefore, the amount of charge left on the surface should decrease in density with decreasing radius. We assume that the infinitesimal charge collected as the contact line slides over an area, 2πr′dr′ (Fig. 5a), equals α times the charge in the electric double layer dQ(r′) = −σs(2πr′)dr′. Here, r′ is the retracted position of the drop on the surface.
![]() | ||
Fig. 5 (a) Schematic of charge deposition for a retracting drop. (b) Drop charge versus the retracted position r′ from eqn (1). The constants σ0S = 14 μC m−2 and λr = 0.6 mm were taking fitting eqn (1) into Fig. 4 at r′ = Rmax for SN-FDTS. |
Assuming that the drop starts neutrally charged at its maximum spreading radius Q(r′ = Rmax) = 0, (see ESI,† Section XI), we find the charge accumulated in the drop after it has retracted to radius r′:
![]() | (1) |
![]() | (2) |
The surface charge density left behind starts at its highest value, σ0 (at r′ = Rmax), and decreases with decreasing radius due to the increasing potential between drop and surface (Fig. S9, ESI†). For larger drops (Rmax = 10 mm), most of the charge left behind by a rebounding drop is concentrated at the outer ring of the deposition pattern (Fig. 6a). The drop fully charges before finishing its retraction, leaving a bare area in the centre of the deposition area because the drop capacitance limits further charge deposition on the surface. Smaller drops do not charge fully (Rmax < 5 mm) and there is charge left across the entire deposition area. This pattern is quite similar to that measured by Nie et al.41 as the charge pattern left behind by a drop pressed between two surfaces.
![]() | ||
Fig. 6 (a) Visualization of the predicted surface charge distribution left behind by three drops with the same volume but different maximum spreading areas. (b) Total drop charge (retracted to r′ = 0) versus maximum spreading area from eqn (3) for PFOTS coated surfaces (σ0S = 12 μC m−2; λr = 2.7 mm, black dotted line, from30) and our SN-FDTS surfaces for 4 μL (σ0S = 14 μC m−2; λr = 0.6 mm, brown dashed line). (c) Experimental data for drop charge versus maximum spreading area. The theoretical prediction for SN-FDTS (shown in (b), brown dashed line), SN and candle soot is plotted on the same axes. The parameters used for the fitting were: σ0S = 5 μC m−2 and λr = 0.52 mm for SN; σ0S = 17 μC m−2 and λr = 0.63 mm for candle soot. All the surface coatings have the fits starting from A0 as the minimum measured area. SN-FDTS data is fitted with an additional curve using A0 = 0.57 mm for a resting drop, as predicted by the model (brown dashed line). |
To compare our drop charge measurements with the proposed theory, we set r′ = 0 because the drop has fully retracted and when it leaves surface. Inserting r′ = 0 into eqn (1), we find the total charge on the drop as a function of Amax:
![]() | (3) |
Our experimental data in Fig. 6c are offset in terms of area compared to the theoretical prediction. Furthermore, charge dissipation processes not included in the theory could partially discharge the drop during extended contact times, reducing the net charge of the drop. This explains the difference in the slope when varying the drop volume compared to varying the impact speed. Since drops retract slower with an increase of volume, more time is available for discharging mechanisms. Another important issue that also explains the offset at low Amax is the surface structure of our substrates. The parameters σ0S and λr in30 were determined for flat fluorinated glass surfaces. These should have the same surface chemistry as our fluorinated superhydrophobic substrates. However, the nanoscale roughness of our substrates should lead to a much smaller effective liquid–solid contact area than the apparent contact area Amax, as the drops remain in the Cassie state during bouncing. In particular, we can estimate the effective contact area between the surface and drop using the Cassie–Baxter equation:33
cos![]() ![]() ![]() | (4) |
Another aspect that may play a role in the charging process is the relative humidity (RH) of the environment. We performed our experiments under lab ambient conditions, so we did not control RH (measured values were 30–40%). Recently, Sosa et al.42 showed that RH affects the charging of sliding drops only when RH > 70%. Thus, we assume that relative humidity does not affect the charging process in this work.
Our model, based on a radial solution of the slide electrification theory, is in good agreement with the experimental results. This supports the idea that a drop deposits charge along the outer ring of the impact area during the retraction phase. Impacting drops occur universally in nature and industry. Understanding and being able to model these charge separation phenomena has broad-ranging impacts. It could allow the optimization of energy-harvesting applications using rebounding drops. It could aid in the electrostatic guidance of charged bouncing drops or it could help us understand the fundamentals of tribocharging in liquids. Future studies for a better understanding of the underlying charge separation mechanism might involve use of other liquids and probing the influence of salts and pH.
Footnote |
† Electronic supplementary information (ESI) available: Sample preparation details, drop impact dynamics, full derivation of the theoretical model of charging for bouncing drops. See DOI: 10.1039/d1sm01725j |
This journal is © The Royal Society of Chemistry 2022 |