ORIGINAL RESEARCH article
Sec. Plant Biophysics and Modeling
Evaluation of the surface free energy of plant surfaces: toward standardizing the procedure
- 1Forest Genetics and Ecophysiology Research Group, Plant Physiology and Anatomy Unit, School of Forest Engineering, Technical University of Madrid, Madrid, Spain
- 2Department of Applied Physics I, Faculty of Physics, Complutense University of Madrid, Madrid, Spain
- 3Madrid Institute for Advanced Studies of Water (IMDEA Water Institute), Madrid, Spain
Plant surfaces have been found to have a major chemical and physical heterogeneity and play a key protecting role against multiple stress factors. During the last decade, there is a raising interest in examining plant surface properties for the development of biomimetic materials. Contact angle measurement of different liquids is a common tool for characterizing synthetic materials, which is just beginning to be applied to plant surfaces. However, some studies performed with polymers and other materials showed that for the same surface, different surface free energy values may be obtained depending on the number and nature of the test liquids analyzed, materials' properties, and surface free energy calculation methods employed. For 3 rough and 3 rather smooth plant materials, we calculated their surface free energy using 2 or 3 test liquids and 3 different calculation methods. Regardless of the degree of surface roughness, the methods based on 2 test liquids often led to the under- or over-estimation of surface free energies as compared to the results derived from the 3-Liquids method. Given the major chemical and structural diversity of plant surfaces, it is concluded that 3 different liquids must be considered for characterizing materials of unknown physico-chemical properties, which may significantly differ in terms of polar and dispersive interactions. Since there are just few surface free energy data of plant surfaces with the aim of standardizing the calculation procedure and interpretation of the results among for instance, different species, organs, or phenological states, we suggest the use of 3 liquids and the mean surface tension values provided in this study.
Plant surfaces play a crucial role in protecting organs against an array of biotic and abiotic stress factors (Riederer, 2006). For instance, they may serve as barrier for the attack of pests and pathogens (Eigenbrode and Jetter, 2002; Serrano et al., 2014), high UV and visible radiation intensities (Reicosky and Hanover, 1978; Karabourniotis and Bormann, 1999), and chiefly the uncontrolled loss of water (Riederer and Schreiber, 2001). On the other hand, plant surface permeability to liquids may be a phenomenon of ecophysological and agronomic significance since it may enable the absorption of, for example, liquid water (Oliveira et al., 2005; Fernández et al., 2014a), fog (Limm and Dawson, 2010; Eller et al., 2013; Berry et al., 2014), dew (Konrad et al., 2012), foliar fertilizers (Fernández and Eichert, 2009; Fernández and Brown, 2013) or biostimulants (Saa et al., 2015). This will be first influenced by the wettability and adhesion or repellence of water drops onto plant surfaces (Fernández et al., 2014a).
Most aerial plant organs are protected by epidermal cells (including stomata and/or trichomes) of different shapes and sizes (Javelle et al., 2011; Fernández and Brown, 2013). The epidermis itself is covered with a cuticle, a layer which may be interpreted as a lipidized area of the epidermal cell wall (Guzmán et al., 2014). The outermost surface of the cuticle is the epicuticular wax layer, which may present a major degree of variability depending on the chemical and structural nature of epicuticular waxes (Barthlott et al., 1998; Jetter et al., 2006), and may also be affected by the prevailing environmental conditions during plant growth (Jeffree et al., 1975; Shepherd and Wynne Griffiths, 2006). Depending on the micro- and nano-structure of epidermal cells, epicuticular waxes and cuticular folds, plant surfaces may have different degrees of roughness and water drop repellence or adhesion (Koch and Barthlott, 2009; Fernández et al., 2014a). During the last decade and in light of the major plant surface variability observed by scanning electron microscopy (SEM) (Barthlott et al., 1998), there is a rising interest in examining plant surface topography and performance when in contact with water drops for the development of bio-mimetic materials (Bhushan and Jung, 2010). The existing epidermal micro-and nano-surface features together with surface chemical composition (Khayet and Fernández, 2012) will determine the degree of plant surface roughness and hydrophobicity (Fernández et al., 2014a).
Despite the fact that the mechanisms of foliar uptake of water and solutes by plant surfaces are still not fully understood (Fernández and Eichert, 2009; Burkhardt and Hunsche, 2013), several studies analyzed the interactions of plant surfaces with water (Brewer et al., 1991; Brewer and Smith, 1997; Hanba et al., 2004; Brewer and Nuñez, 2007; Dietz et al., 2007; Roth-Nebelsick et al., 2012; Rosado and Holder, 2013; Urrego-Pereira et al., 2013; Wang et al., 2014). Such investigations largely focussed on analysing the wettability and hydrophobicity of plant surfaces based on water contact angles. However, Fernández et al. (2011, 2014a) and Khayet and Fernández (2012) suggested the implementation of membrane science approaches for analysing the physico-chemical properties of plant surfaces and introduced surface free energy results derived from the measurement of contact angles of 3 liquids with different degrees of dispersive and non-dispersive interactions according to the method of Van Oss et al. (1987, 1988).
The surface tensions of the interfaces of a solid with liquids and the surrounding vapors are key parameters, which are commonly applied for characterizing the physico-chemical properties of surfaces in many areas of applied science and technology (Correia et al., 1997). The theory of interfacial tensions (γ) and surface free energy (SFE) was first introduced by Young (1805). Subsequently, several theoretical approaches were developed between the 1950's and the 1980's (e.g., Fox and Zisman, 1950; Girifalco and Good, 1957; Fowkes, 1964; Owens and Wendt, 1969; Wu, 1971; Good, 1977; Van Oss et al., 1986, 1988) with the aim of estimating solid surface tension (γs) or SFE which is the same) from contact angles of 2 or 3 liquids with different dispersive and non-dispersive components.
Despite the interest in estimating the SFE of different solids ranging from polymers and co-polymer blends (e.g., Correia et al., 1997; Adão et al., 1999), to polymeric membranes (Khayet et al., 2003, 2007), geological materials (Chibowski and Staszczuk, 1988; Helmy et al., 2004), or wood (Gindl et al., 2001), there is a fair degree of controversy concerning several associated aspects such as the liquid surface tension and contact angle determination methods, the use of 2, 3 or more test liquids and the different calculation methods reported (Morra, 1996; Della Volpe and Siboni, 1997; Jañczuk et al., 1999; Kwok and Neumann, 1999; Tretinnikov, 2000; Chibowski and Perea-Carpio, 2002; Della Volpe et al., 2004).
Since there are still few SFE data of plant materials (Gorb et al., 2004; Fernández et al., 2011, 2014a,b; Khayet and Fernández, 2012; Wang et al., 2014), and current contact angle measuring devices are available with software enabling: (i) the direct estimation of the SFE of solids and their components according to different methods, and (ii) the selection of different test liquid surface tension values, we carried out preliminary efforts to revise the existing SFE theory together with the most commonly used SFE determination methods, and to calculate the mean values for the total surface tension and surface tension components of water (W), glycerol (G) and diiodomethane (DM), which have been often used for such purpose. It must be stressed that SFE estimations of a solid surface by contact angle measurements of either 2 or 3 test liquids of the same nature should be similar or at least within the same range as reported in several studies performed with some synthetic materials (e.g., Adão et al., 1999; Khayet et al., 2003). However, the major structural (roughness) and chemical (polarity, apolarity and hydrogen-bonding interactions) variation observed among plant species, organs or organ parts made us develop trials to calculate the SFE of model plant surfaces when measured with the most commonly used SFE estimation methods, which employ either 2 or 3 test liquids. With the aim of choosing a procedure for determining the SFE of plant materials to be systematically used in future research attempts, we further calculated this parameter using contact angle measurements of either pairs of liquids (i.e., W-DM, G-DM, W-G) considering both the geometric and harmonic mean approaches, or 3 liquids (W-G-DM) applying the Lifshitz-van der Waals-acid-base or van Oss, Good and Chaudhury method (Van Oss et al., 1986, 1987, 1988). The contact angles of W, G and DM were determined on to 3 rough and 3 rather flat model plant surfaces and the SFE obtained by different methods were compared.
Materials and Methods
The plant materials analyzed as model correspond to intact, mature red ironbark (Eucalyptus sideroxylon A. Cunn. ex Woolls) leaves, juvenile blue gum eucalypt (Eucalyptus globulus Labill.) leaves, mature Chilean myrtle (Luma apiculata (DC.) Burret) leaves, mature rubber tree (Ficus elastica Roxb. ex Hornem) leaves, mature holm-oak (Quercus ilex subsp. ballota (Desf.) Samp.) leaves, and red bell peppers (Capsicum annum L. cv. “Genil”). Rubber tree and Chilean myrtle leaves and also pepper fruit surfaces are smooth while red iron bark, blue gum eucalypt and holm oak leaf surfaces are extremely rough. We examined the properties of the upper leaf side of red ironbark blue gum eucalypt, Chilean myrtle, and rubber tree. The lower leaf side of holm-oak and the pepper fruit surface were also analyzed.
Scanning Electron Microscopy
Gold sputtered plant materials were examined with a Hitachi S-3400 N (Tokyo, Japan) and a Philips XL30 (Eindhoven, The Netherlands) scanning electron microscope (SEM).
Contact Angle Measurements
The test liquids analyzed were double-distilled water (W), glycerol (G) and diiodomethane (DM, both 99% purity, Sigma-Aldrich). Advancing contact angles of drops of such liquids were measured at 20°C using a contact angle meter CAM 200 (KSV Instruments Ltd., Helsinki, Finland). Drops of each liquid were deposited on to intact plant surfaces (30 repetitions) with a manual dosing system holding a 1 mL syringe (0.5 mm diameter needle). Side view images of the drops were captured at a rate of 10 frames s−1. Contact angles were automatically calculated by fitting the captured drop shape to the one calculated from the Young–Laplace equation.
Surface Free Energy Calculation Methods
Since the total SFE, γs, cannot be calculated directly (Khayet et al., 2003), different indirect methods have been proposed which based on the contact angle measurement of either 2 or 3 different liquids deposited on to a solid surface (Jañczuk and Białopiotrowicz, 1989; Della Volpe and Siboni, 1997). Three calculation methods were used for estimating γs, its components and the resulting solubility parameter (δ) of different plant materials. Provided that we observed a significant dispersion of the total surface tension and surface tension component data provided in the literature for various liquids including W, G and DM, the mean values of all results gathered were considered in this study (Table 1).
Table 1. Total surface free energy or surface tension (γl) and its components (measured at 20°C) of water (W), Glycerol (G), and diiodomethane (DM).
Calculations Based on 2 Liquids
The 2-liquid based, geometric mean and harmonic mean methods of calculation of SFE are commonly used for characterizing synthetic surfaces (e.g., Correia et al., 1997; Adão et al., 1999) and have been used in this investigation.
Based on the theory of Girifalco and Good (1957), Fowkes (1962, 1964) proposed that γs (SFE of the solid to be analyzed) is the sum of contributions from the different intramolecular forces at the surface: those due to dispersion interactions, dipole-dipole interactions, dipole induced-dipole interactions, hydrogen bonding, Π-bonding, electrostatic interactions and acceptor-donor interactions.
While there is agreement in the nomenclature for the dispersive SFE component (γds), many authors referred to the remaining term (i.e., γnds), which obviously includes many types of forces (Fowkes, 1964), as polar component γps (e.g., Wu, 1971). To avoid any possible misinterpretation, we will subsequently refer to the γds and γnds components throughout the manuscript.
Geometric Mean Method
Based on Fowkes' assumption that dispersive interactions between the molecules of the solid and the liquid prevailed among other acting forces, and considering Young's equation, the following expression was suggested (Fowkes, 1964):
where, γl is the total surface tension of the liquid, γdl is the dispersive component of the liquid and θ is the contact angle measured between the liquid and the solid under study.
This equation was further expanded by Owens and Wendt (1969) to:
Where, γndl is the non-dispersive component of the liquid and γnds is the non-dispersive component of the solid.
Equation (3) is known as the Geometric mean, SFE calculation method, termed hereafter GM method.
Harmonic Mean Method
Wu (1971) proposed an alternative equation, based on the reciprocal mean and force additivity, which once applied to Young's equation for two liquids results in:
Equation (4) is known as the Harmonic mean, SFE calculation method, termed hereafter HM method.
where i denotes either the solid or the liquid phase and the acid-base component (γAB) breaks down into the electron-donor (γ−) and the electron-acceptor (γ+) interactions. For a solid-liquid system, the following expression is given (Van Oss et al., 1987, 1988):
where the 3 surface tension components (i.e., γLWs, γ+s and γ−s) together with the SFE of the solid measured (γs) can be obtained. This SFE calculation procedure is termed hereafter 3-L method.
Solubility Parameter Calculations
From the total surface tension and surface tension components estimated by the methods described above, the solubility parameter (δ) of plant surfaces was calculated using the following relation (Khayet et al., 2002):
where ec (MJ m−3) is the cohesive energy density, which is related to γs (mJ m−2) as follows:
Liquid Surface Tension Mean Values
The average values of the total surface tension γl and its components for the 3 liquids measured (W, G, and DM) were calculated together with their standard deviation using the data reported in the literature. For the GM and the HM, 2-liquids methods, the surface tension components are γdl and γndl while for the 3-L method they are γLW, γ−, and γ+ (Table 1). The estimation of γdl either via measuring contact angles on to a known solid or by liquid-liquid interfacial tension measurements may differ (Jañczuk et al., 1993). This was chiefly noticeable in the case of G.
Beginning with the 2-liquid approaches (i.e., the GM and HM methods), the mean γl, γdl, and γndl values of W and DM were obtained from: Fowkes (1964), Owens and Wendt (1969), Wu (1971, 1982), and Jañczuk and Białlopiotrowicz (1989). Glycerol surface values are the average of the results reported by Fowkes (1964), Dann (1970), Panzer (1973), Wu (1971, 1982), Jañczuk and Chibowski (1983) and Jañczuk et al. (1993).
The Liftshizt-van der Waals and acid-base components (including γ− and γ+) of W were taken as reference by Van Oss et al. (1987). The γLW, γ− and γ+ mean values of G and DM were calculated from the data reported by van Oss (1994) and Jañczuk et al. (1993, 1999). For G, it is remarkable the high degree of dispersion of the γ+ and chiefly the γ− components (see Table 1).
Plant Surface Topography
The adaxial leaf surface of red ironbark and blue gum eucalypt, and abaxial leaf side of holm oak were found to have a high degree of roughness (Figures 1A–C). This is due to the presence of milimetric trichomes on to the lower leaf side of holm-oak, wax nano-tubes on the surface of blue gum eucalypt leaves, and wax platelets and a complex micro-topography provided by the epidermal cells of red ironbark leaves. By contrast, the upper leaf surface of rubber tree and Chilean myrtle is rather flat alike that of pepper fruit (Figures 1D–F).
Figure 1. Plant surfaces analysed. (A) Upper leaf side of red ironbark (major micro- and nano roughness), (B) upper leaf side of juvenile blue gum (major micro- and nano roughness conferred by wax nano-tubes), (C) lower leaf side of holm oak (great micro-roughness conferred by the hairs/trichomes), (D) upper leaf side of rubber tree (rather smooth), (E) upper leaf side of Chilean myrtle (rather smooth), and (F) pepper fruit surface (rather smooth).
The contact angle values of W, G, and DM with the surface of the examined plant materials are shown in Table 2. In Figure 2 the contact angles of water (W), glycerol (G), and diiodomethane (DM) with one of the rough (blue gum leaf) and smooth (pepper fruit) surfaces analyzed are shown as an example. The structural complexity observed in the adaxial surfaces of red ironbark, blue gum eucalypt, and abaxial holm-oak leaf surfaces led to very high contact angles with polar liquids (i.e., W and G; Figures 2A,B). While the blue gum eucalypt leaf is wettable (θ < 90°) for DM (Figure 2C), which suggests the occurrence of chemical interactions between such liquid and epicuticular wax nano-tubes, the surfaces of red ironbark and holm oak leaves are unwettable (θ > 90°) for such largely apolar liquid.
Table 2. Contact angles of water (θw), glycerol (θg) and diiodomethane (θd) on the upper side of red ironbark, blue gum eucalypt, rubber tree and Chilean myrtle leaves, lower side of holm oak leaves, and pepper fruit surfaces.
Figure 2. Contact angles of water (A,D), glycerol (B,E) and diiomethane (C,F) on to blue gum (rough surface; A–C) and pepper fruit (smooth surface; D–F) surfaces, as an example.
Pepper fruit and rubber gum leaf surfaces are wetted by W, G and DM drops (see Figures 2D–F as an example). The surface of Chilean myrtle is however unwettable for W and G, but the contact angles measured (approximately 100°) are well below those determined for red ironbark, blue gum eucalypt, and holm oak leaves.
Surface Free Energy Calculation by Different Procedures
The SFE results calculated for the mean surface tension values (Table 1) and contact angles (Table 2) according to the 3-L, GM or HM methods are grouped in Tables 3–8 for the different plant materials assessed. It must be highlighted that the 3 model smooth surfaces analyzed had lower contact angles chiefly with polar liquids and higher SFE values, which were within a similar range (approximately from 31 to 38 mJ m−2) at least when determined by the 3-L method. The rough surfaces examined had very high contact angles with polar liquids (see Figure 2 as an example) and SFE varying from 5 to 19 mJ m−2 when calculated by the 3-Liquids method.
Table 3. Surface free energy (γs), solubility parameter (δ), surface free energycomponents, i.e., Lifschitz-van der Waals or dispersive component (γLW or ds) acid-base or non-dispersive component (γAB or nds), electron-donor (γ−s) and the electron-acceptor (γ+s) components), of adaxial red ironbark leaf surfaces calculated by the 3-liquids (3-L), geometric mean (GM) and harmonic mean (HM) methods.
Since we understand that the characterization of surfaces by measuring contact angles of 3 liquids with different dispersive and non-dispersive contributions (W, G, and DM) provides more information than comparing pairs of liquids (W-DM, G-DM, or W-G), we analyzed the SFE results taking the 3-L method as reference. Furthermore, since we used average surface tension and surface tension component values, which are comparable for all the SFE estimation approaches employed, it could be expected that, for the same plant material, at least some of the SFE results derived from pairs of liquids would be within the range obtained by the 3-L method. Based on the developed theory (Equations 3–6) it would be desirable that data estimated for γLWs and γds, γABs and γnds, and also γs and δ would be comparable among methods. For assessing the potential similarity between the results derived from 2 and 3 liquid calculation procedures, the percentage deviation from de 3-L method is provided in Tables 3–8.
The upper leaf surface of red iron bark, which is quite rough (Figure 1A), was found to have quite low total SFE (6.42 mJ m−2) and δ (5.01 MJ½ m−3/2) values which are however within the range determined for holm oak abaxial surfaces. The γLWs is lower than the γABs component chiefly due to the high γ−s value (Table 3). While selecting W-DM data and calculating the SFE by the HM method, γs and δ values do not deviate so much for the 3-L method results (−9.0 and 7.0% deviation, respectively). However, this approach led to a major over-estimation of the γLWs component. The remaining GM and HM results have a high degree of dispersion when compared to the 3-L method.
The adaxial side of juvenile blue gum leaves has also a high degree of roughness (Figure 1B). The total SFE calculated by the 3-L method is 19.07 mJ m−2, principally due to a large γLWs compared to the γABs component (Table 4). The SFE and δ results obtained with W-DM and the GM method had the lowest deviation (from 12 to 27% for all parameters) from the 3-L method. The remaining results obtained with the GM and HM methods led to significant SFE and δ under or over-estimations in relation to the 3-L method values.
Table 4. Surface free energy (γs) and its components (γLW or ds γAB or nds, γ+s and γ−s), and solubility parameter (δ) of adaxial blue gum eucalypt leaf surfaces calculated by the 3-L, GM, and HM methods.
The lower leaf surface of holm oak has a high degree of micro-roughness provided by pubescence (Figure 1C), and low total SFE (5. 97 mJ m−2) and δ (4.74 MJ½ m−3/2) values (Table 5). When considering W-DM data and using the HM method, γs and δ values presented a low deviation from the 3-L method (11.1 and 8.2%, respectively). However, the HM approach led to a major over-estimation of γds (189.1%) and an under-estimation (−80.0%) of the γnds components, as also recorded for red ironbark (Table 3). The remaining GM and HM values significantly deviated from the results obtained using the 3-L method (Table 5).
Table 5. Surface free energy (γs) and its components (γLW or ds γAB or nds, γ+s and γ−s), and solubility parameter (δ) of abaxial holm oak leaf surfaces calculated by the 3-L, GM and HM methods.
The upper leaf surface of rubber tree is rather smooth (Figure 1D) and has a higher total SFE (principally associated with a high γLWs and a significant γ−scomponent) and δ value (31.76 mJ m−2 and 16.60 MJ½ m−3/2, respectively). The lowest deviation from the results obtained by the 3-L method was obtained for W-DM and the GM approach (from 2.0 to 3.9% deviation for γds, γs, and δ), which however led to a 24.3% under-estimation of γnds (Table 6). Using W-DM and G-DM pairs and calculating γs by the GM or HM approaches led to results with low deviations from the 3-L method, with the exception of the γnds component which was under-estimated in the case of W-DM and largely G-DM using the GM method (−24.3 and −85.9%, respectively), and G-DM calculated with the HM method (−60.5% deviation).
Table 6. Surface free energy (γs) and its components (γLW or ds γAB or nds, γ+s and γ−s), and solubility parameter (δ) of adaxial rubber tree leaf surfaces calculated by the 3-L, GM and HM methods.
The rather flat upper leaf surface of Chilean myrtle (Figure 1E) has the highest total SFE (chiefly due to a high γLWs and a significant γ−s component) and δ values of all the plant materials analyzed (38.12 mJ m−2 and 19.04 MJ½ m−3/2, respectively). The lowest γds, γs, and δ deviations from the 3-L method were obtained for G-DM using the GM method, followed by W-DM calculated by the HM and GM approaches (Table 7). However, all of these results based on pairs of liquids led to a major under-estimation of the γnds component (ranging from −80 to −96% deviations).
Table 7. Surface free energy (γs) and its components (γLW or ds γAB or nds, γ+s and γ−s), and solubility parameter (δ) of adaxial Chilean myrtle leaf surfaces calculated by the 3-L, GM, and HM methods.
Finally, the rather smooth pepper fruit surface (Figure 1F) has also high γs (principally associated with a large γLWs component) and δ values (31.66 mJ m−2 and 16.56 MJ½ m−3/2, respectively). Using W-DM contact angles following the GM approach led to γs and δ data with low deviations from the 3-L method (with 3.1–7.4% under-estimations for γds, γs and δ and 11.4% over-estimation for γnds; Table 8). For W-G, and the GM method, γs and δ results slightly deviated from the 3-L method (8–18% over-estimations for γds, γs and δ and 23.2% under-estimation of γnds). Using G-DM and W-DM with the GM method, and even W-DM, G-DM and W-G with the HM approach also led to γds, γs and δ values with low deviations. However, such 2-liquid combinations led to a major over-estimation of the γnds component (from 40 to 100%) as compared to the results obtained by the 3-L method.
Table 8. Surface free energy (γs) and its components (γLW or ds γAB or nds, γ+s and γ−s), and solubility parameter (δ) of pepper fruit surfaces calculated by the 3-L, GM, and HM methods.
Calculations by the HM method led to some negative results for blue gum and holm oak leaves or to no results in the case of G-DM and G-W combinations for red iron bark and also for G-DM calculations for Chilean myrtle and holm oak. These results have no physical significance and are related to mathematical constraints associated with HM method equations.
The SFE of a material is a fundamental property that determines its surface and interfacial performance in processes like wetting and adhesion (Adão et al., 1999). Despite its importance for characterizing contact phenomena related to plant surfaces, such as the interaction of plant organs with surface deposited water (Brewer et al., 1991) or insects (Prüm et al., 2012), only few studies considered this parameter in a plant science context so far (Gorb et al., 2004; Fernández et al., 2011, 2014a,b; Khayet and Fernández, 2012; Wang et al., 2014). The estimation of the SFE and solubility parameter of plant materials may be, for example, useful for predicting water-plant surface interactions and the potential of plant organs to absorb water and agrochemical solutions (Fernández et al., 2011, 2014a,b).
Several SFE calculation methods were introduced during the last 50 years and the drawbacks of estimating SFE based on contact angle measurements by different procedures have been highlighted in some reports (e.g., Kwok and Neumann, 1999; Chibowski and Perea-Carpio, 2002; Della Volpe et al., 2004), which may be minimized if carefully considering the procedures, at least for comparing between different plant materials, which is the aim of this study.
Surface roughness can affect liquid contact angles and it is a matter of interest for many researchers dealing with solid-liquid interactions (Bhushan and Nosonovsky, 2010). According to the models of Wenzel (1936) and Cassie and Baxter (1944), there are two regimes of wetting of a rough surface: a homogeneous regime with a 2-phase solid–liquid interface, and a non-homogeneous regime with a 3-phase solid–water–air interface (with air pockets between the solid and the liquid).
When revising the existing SFE-related multi-disciplinary literature, we recognized that an array of test liquids with different total surface tension and surface tension components were available. We focused on water, glycerol and diiodomethane which are commonly used due to their distinct dispersive and non-dispersive components (e.g., Adão et al., 1999; Jañczuk et al., 1999). It is worth highlighting that for test liquids, we recognized that different surface tension values were used by authors when calculating SFE for instance, by the GM vs. the HM approach and that the method of determination of the dispersive component (e.g., with contact angles or interfacial tension measurements) may also have an influence on the surface tension values of liquids. Additionally, we actually noticed that the calculation methods were extremely sensitive to slight surface tension and surface tension component value modifications. Therefore, we decided to work with average surface tension and surface tension component values, which may be comparable among different SFE calculation methods as shown in Table 1. This choice enabled us to trust that for the same plant surface, the potential SFE differences are not due to the different surface tension values of test liquids, but rather to the particular calculation method in combination with the surface features of the materials analyzed. Thereby, we standardized the surface tension values of W, G and DM by using the mean values calculated after revising the exiting literature, and subsequently used them for estimating the SFE of plant materials.
Aware that any of the SFE calculation methods based on either 2 or 3 liquids analyzed produced certain SFE and SFE component values, we wondered if the results were meaningful and comparable between and within species. For verifying this aspect, 3 rough and 3 smother plant materials were selected as model species, and contact angles with W, G, and DM were carefully measured under similar conditions.
Since we understand that the characterization of a surface by measuring drops of 3 test liquids with different properties is more informative than using only 2 liquids having different dispersive and non-dispersive interactions, we took the results obtained using the 3-L method as reference and estimated the potential deviation of SFE, SFE components and solubility parameter (δ) values. This last parameter was theoretically calculated by Khayet and Fernández (2012) for model wax compounds, and provides insight into the combined effect of surface chemical composition (plant surfaces are generally covered with waxes having a theoretical δ ranging between 16 and 17 MJ½ m−3/2) and roughness. Looking at the mean surface tension components shown in Table 1, it can be derived that in a decreasing order, W has the highest non-dispersive component, followed by G and DM. On the contrary, DM has the highest dispersive (apolar) component, followed by G and W. It is hence easy to understand that drops of such liquids may perform differently depending on the chemical and structural nature of every different surface, as observed in this study. While measuring contact angles of drops of 2 liquids may save time when determining SFE of samples either by the GM or HM methods, we must be sure that results are meaningful and comparable with those obtained by the 3-L method.
We could not find a clear relationship between surface roughness and SFE calculation methods, since for example, for the blue gum eucalypt leaf which is quite rough, W-DM contact angles estimated by the GM method provided fair SFE and δ estimations. On the contrary, for the rather smooth Chilean myrtle leaf, none of the pairs of liquids combinations and 2-liquid calculation methods lead to values within the range obtained with the 3-L method, chiefly due to a severe under-estimation of γnds. Our observations agree with the comments of Kwok and Neumann (1999), who noted that “there are as yet no general criteria to answer the question of how smooth a solid surface has to be for surface roughness not to have an effect on the contact angle.”
The combination W-DM calculated with the GM method provided fair SFE and δ results for rubber tree, blue gum leaves and pepper fruit. Most of the 2-liquid combinations and calculation methods generally led to γnds under-estimations and sometimes γds over-estimations for red iron bark, holm oak and Chilean myrtle compared to the 3-L results. The HM method generally led to the highest deviations from the 3-L method as compared to the GM method. According to Wu (1971, 1982) the HM approach has been described to be suitable for analysing surfaces with low SFE and not for applying it to systems involving phases with widely different polarizabilities. Limitations associated with the use of the HM method for assessing the SFE of surfaces having some degree of polarity have been reported earlier (e.g., Wu, 1971; Adão et al., 1999).
The results obtained for blue gum and holm oak leaves, and pepper fruit surfaces with the 3-L method and mean liquid surface tension values shown in Table 1, are within the range reported by Khayet and Fernández (2012) and Fernández et al. (2014a). The most remarkable difference is that we have considered the γ+l of DM, which may be crucial for surfaces having significant polar components as noted by Tretinnikov (2000). The δ empirically determined for the smooth rubber tree and pepper fruit surfaces are coincident with the theoretical values estimated by Khayet and Fernández (2012), hence indicating a negligible effect of surface roughness and the prevalence of chemical interactions between such plant materials and liquid drops. On the contrary, the lower δ values determined for the rough iron bark, holm oak and to a lower degree, blue gum leaves reflect a significant effect of surface micro- and or-nano-structure on plant surface-liquid interactions.
An analysis of the total surface tension and surface tension component values of test liquids reported in the literature led us to select the mean values for water (W), glycerol (G), and diiodomethane (DM), and enabled the comparison between different surface free energy calculation methods. Despite it would be less time-consuming to be able to measure only 2 liquids for characterizing plant surfaces of unknown physical and chemical properties, we gained evidence that it is not possible to a priori select a combination of 2 liquids and a suitable calculation method (GM or HM methods), which may provide results comparable to those obtained using the Lifshitz-van der Waals-acid-base method. We hence suggest the standard use of the gamma values calculated for W, G and DM (Table 1) and the 3-L method for drawing coherent conclusions concerning the surface free energy of plant surfaces, related for instance, to different species, organs, developmental stages or growing conditions.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
VF is supported by a Ramón y Cajal contract (MINECO, Spain), co-financed by the European Social Fund. Thanks to Ricardo Fernández for his technical support for the development of this study.
Adão, M. H. V. C., Saramago, B. J. V., and Fernandes, A. C. (1999). Estimation of the surface properties of styrene-acrylonitrile random copolymers from contact angle measurements. J. Colloid Interface Sci. 217, 94–106. doi: 10.1006/jcis.1999.6279
Barthlott, W., Neinhuis, C., Cutler, D., Ditsch, F., Meusel, I., Theisen, I., et al. (1998). Classification and terminology of plant epicuticular waxes. Bot. J. Linn. Soc. 126, 237–260. doi: 10.1111/j.1095-8339.1998.tb02529.x
Berry, Z. C., Hughes, N. M., and Smith, W. K. (2014). Cloud immersion: an important water source for spruce and fir saplings in the southern Appalachian Mountains. Oecologia 174, 319–326. doi: 10.1007/s00442-013-2770-0
Bhushan, B., and Jung, Y. C. (2010). Natural and biomimetic artificial surfaces for superhydrophobicity, self- cleaning, low adhesion, and drag reduction. Prog. Mat. Sci. 56, 1–108. doi: 10.1016/j.pmatsci.2010.04.003
Brewer, C. A., Smith, W. K., and Vogelmann, T. C. (1991). Functional interaction between leaf trichomes, leaf wettability and the optical properties of water droplets. Plant Cell Environ. 14, 955–962. doi: 10.1111/j.1365-3040.1991.tb00965.x
Correia, N. T., Ramos, J. J. M., Saramago, B. J., and Calado, J. C. (1997). Estimation of the surface tension of a solid: application to a liquid crystalline polymer. J. Colloid Interface Sci. 189, 361–369. doi: 10.1006/jcis.1997.4857
Dann, J. R. (1970). Forces involved in the adhesive process: I. Critical surface tensions of polymeric solids as determined with polar liquids. J. Colloid Interface Sci. 32, 302–320. doi: 10.1016/0021-9797(70)90054-8
Della Volpe, C., Maniglio, D., Brugnara, M., Siboni, S., and Morra, M. (2004). The solid surface free energy calculation: I. In defense of the multicomponent approach. J. Colloid Interface Sci. 271, 434–453. doi: 10.1016/j.jcis.2003.09.049
Dietz, J., Leuschner, C., Hölscher, D., and Kreilein, H. (2007). Vertical patterns and duration of surface wetness in an old-growth tropical montane forest, Indonesia. Flora 202, 111. doi: 10.1016/j.flora.2006.03.004
Eller, C. B., Lima, A. L., and Oliveira, R. S. (2013). Foliar uptake of fog water and transport belowground alleviates drought effects in the cloud forest tree species, Drimys brasiliensis (Winteraceae). New Phytol. 199, 151–162. doi: 10.1111/nph.12248
Fernández, V., and Eichert, T. (2009). Uptake of hydrophilic solutes through plant leaves: current state of knowledge and perspectives of foliar fertilization. Crit. Rev. Plant Sci. 28, 36–68. doi: 10.1080/07352680902743069
Fernández, V., Guzmán, P., Peirce, C. A. E., McBeath, T. M., Khayet, M., and McLaughlin, M. J. (2014b). Effect of wheat phosphorus status on leaf surface properties and permeability to foliar applied phosphorus. Plant Soil 384, 7–20. doi: 10.1007/s11104-014-2052-6
Fernández, V., Khayet, M., Montero-Prado, P., Heredia-Guerrero, J. A., Liakopoulos, G., Karabourniotis, G., et al. (2011). New insights into the properties of pubescent surfaces: peach fruit as a model. Plant Physiol. 156, 2098–2108. doi: 10.1104/pp.111.176305
Fernández, V., Sancho-Knapik, D., Guzmán, P., Peguero-Pina, J. J., Gil, L., Karabourniotis, G., et al. (2014a). Wettability, polarity and water absorption of holm oak leaves: effect of leaf side and age. Plant Physiol. 166, 168–180. doi: 10.1104/pp.114.242040
Fowkes, F. M. (1962). Determination of interfacial tensions, contact angles, and dispersion forces in surfaces by assuming additivity of intermolecular interactions in surfaces. J. Phys. Chem. 66, 382–382. doi: 10.1021/j100808a524
Gindl, M., Sinn, G., Gindl, W., Reiterer, A., and Tschegg, S. (2001). A comparison of different methods to calculate the surface free energy of wood using contact angle measurements. Coll. Surf. A 181, 279–287. doi: 10.1016/S0927-7757(00)00795-0
Girifalco, L. A., and Good, R. J. (1957). A theory for the estimation of surface and interfacial energies. I. Derivation and application to interfacial tension. J. Phys. Chem. 61, 904–909. doi: 10.1021/j150553a013
Gorb, E., Kastner, V., Peressadko, A., Arzt, E., Gaume, L., Rowe, N., et al. (2004). Structure and properties of the glandular surface in the digestive zone of the pitcher in the carnivorous plant Nepenthes ventrata and its role in insect trapping and retention. J. Exp. Biol. 207, 2947–2963. doi: 10.1242/jeb.01128
Guzmán, P., Fernández, V., Graça, J., Cabral, V., Kayali, N., Khayet, M., et al. (2014). Chemical and structural analysis of Eucalyptus globulus and E. camaldulensis leaf cuticles: a lipidized cell wall region. Front. Plant Sci. 5:481. doi: 10.3389/fpls.2014.00481
Hanba, Y. T., Moriya, A., and Kimura, K. (2004). Effect of leaf surface wetness and wettability on photosynthesis in bean and pea. Plant Cell Environ. 27, 413–421. doi: 10.1046/j.1365-3040.2004.01154.x
Jañczuk, B., and Białlopiotrowicz, T. (1989). Surface free-energy components of liquids and low energy solids and contact angles. J. Colloid Interface Sci. 127, 189–204. doi: 10.1016/0021-9797(89)90019-2
Jañczuk, B., Białopiotrowicz, T., and Wójcik, W. (1989). The components of surface tension of liquids and their usefulness in determinations of surface free energy of solids. J. Colloid Interface Sci. 127, 59–66. doi: 10.1016/0021-9797(89)90007-6
Jañczuk, B., Wójcik, W., and Zdziennicka, A. (1993). Determination of the components of the surface tension of some liquids from interfacial liquid-liquid tension measurements. J. Colloid Interface Sci. 157, 384–393. doi: 10.1006/jcis.1993.1200
Javelle, M., Vernoud, V., Rogowsky, P. M., and Gwyneth, C. I. (2011). Epidermis: the formation and functions of a fundamental plant tissue. New Phytol. 189, 17–39. doi: 10.1111/j.1469-8137.2010.03514.x
Jetter, R., Kunst, L., and Samuels, A. L. (2006). “Composition of plant cuticular waxes,” in Biology of the Plant Cuticle, Annual Plant Reviews, Vol. 23, eds M. Riederer and C. Müller (Oxford: Blackwell), 145–181.
Karabourniotis, G., and Bormann, J. F. (1999). Penetration of UV-A UV-B and blue light through the leaf trichome of two xeromorphic plants, olive and oak, measured by optical fibre microprobes. Physiol. Plantarum 105, 655–661. doi: 10.1034/j.1399-3054.1999.105409.x
Khayet, M., Feng, C. Y., and Matsuura, T. (2003). Morphological study of fluorinated asymmetric polyetherimide ultrafiltration membranes by surface modifying macromolecules. J. Membr. Sci. 213, 159–180. doi: 10.1016/S0376-7388(02)00523-9
Khayet, M., and Fernández, V. (2012). Estimation of the solubility parameters of model plant surfaces and agrochemicals: a valuable tool for understanding plant surface interactions. Theor. Biol. Med. Model. 9, 45. doi: 10.1186/1742-4682-9-45
Khayet, M., Vazquez Alvarez, M., Khulbe, K. C., and Matsuura, T. (2007). Preferential surface segregation of homopolymer and copolymer blend films. Surf. Sci. 601, 885–895. doi: 10.1016/j.susc.2006.11.024
Konrad, W., Ebner, M., Traiser, C., and Roth-Nebelsick, A. (2012). Leaf surface wettability and implications for drop shedding and evaporation from forest canopies. Pure Appl. Geophys. 169, 835–845. doi: 10.1007/s00024-011-0330-2
Limm, E. B., and Dawson, T. E. (2010). Polystichum munitum (Dryopteridaceae) varies geographically in its capacity to absorb fog water by foliar uptake within the redwood forest ecosystem. Am. J. Bot. 97, 1121–1128. doi: 10.3732/ajb.1000081
Oliveira, R. S., Dawson, T. E., and Burgess, S. O. (2005). Evidence for direct water absorption by the shoot of the dessication-tolerant plant Vellozia flavicans in the savannas of central Brazil. J. Tropical Ecol. 21, 585–588. doi: 10.1017/S0266467405002658
Roth-Nebelsick, A., Ebner, M., Miranda, T., Gottschalk, V., Voigt, D., Gorb, S., et al. (2012). Leaf surface structures enable the endemic Namib desert grass Stipagrostis sabulicola to irrigate itself with fog water. J. R. Soc. Interface 9, 1965–1974. doi: 10.1098/rsif.2011.0847
Saa, S., Olivos-Del Rio, A., Castro, S., and Brown, P. H. (2015). Foliar application of microbial and plant based biostimulants increases growth and potassium uptake in almond (Prunus dulcis [Mill.] D. A. Webb). Front. Plant Sci. 6:87. doi: 10.3389/fpls.2015.00087
Tretinnikov, O. N. (2000). On neglecting the polar nature of halogenated hydrocarbons in the surface energy determination of polar solids from contact angle measurements. J. Colloid Interface Sci. 229, 644–647. doi: 10.1006/jcis.2000.7024
Urrego-Pereira, Y. F., Martínez-Cob, A., Fernández, V., and Cavero, J. (2013). Daytime sprinkler irrigation effects on net photosynthesis of maize and alfalfa. Agron. J. 105, 1515–1528. doi: 10.2134/agronj2013.0119
Van Oss, C. J., Good, R. J., and Chaudhury, M. K. (1986). The role of van der Waals forces and hydrogen bonds in “hydrophobic interactions” between biopolymers and low energy surfaces. J. Colloid Interface Sci. 111, 378–390. doi: 10.1016/0021-9797(86)90041-X
Wang, H., Shi, H., Li, Y., and Wang, Y. (2014). The effects of leaf roughness, surface free energy and work of Adhesion on leaf water drop adhesion. PLoS ONE 9:e107062. doi: 10.1371/journal.pone.0107062
Keywords: contact angles, cuticle, geometric mean, plant surfaces, surface free energy, three-liquids method
Citation: Fernández V and Khayet M (2015) Evaluation of the surface free energy of plant surfaces: toward standardizing the procedure. Front. Plant Sci. 6:510. doi: 10.3389/fpls.2015.00510
Received: 26 March 2015; Accepted: 25 June 2015;
Published: 07 July 2015.
Edited by:Michael Knoblauch, Washington State University, USA
Reviewed by:Kaare Hartvig Jensen, Harvard University, USA
Eric M. Kramer, Bard College at Simon's Rock, USA
Copyright © 2015 Fernández and Khayet. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Victoria Fernández, Forest Genetics and Ecophysiology Research Group, Plant Physiology and Anatomy Unit, School of Forest Engineering, Technical University of Madrid, Ciudad Universitaria s/n, 28040 Madrid, Spain, email@example.com; firstname.lastname@example.org;
Mohamed Khayet, Department of Applied Physics I, Faculty of Physics, Complutense University of Madrid, Avda. Complutense s/n, 28040 Madrid, Spain, email@example.com