Optimization principles and the figure of merit for triboelectric generators

See allHide authors and affiliations

Science Advances  15 Dec 2017:
Vol. 3, no. 12, eaap8576
DOI: 10.1126/sciadv.aap8576


Energy harvesting with triboelectric nanogenerators is a burgeoning field, with a growing portfolio of creative application schemes attracting much interest. Although power generation capabilities and its optimization are one of the most important subjects, a satisfactory elemental model that illustrates the basic principles and sets the optimization guideline remains elusive. We use a simple model to clarify how the energy generation mechanism is electrostatic induction but with a time-varying character that makes the optimal matching for power generation more restrictive. By combining multiple parameters into dimensionless variables, we pinpoint the optimum condition with only two independent parameters, leading to predictions of the maximum limit of power density, which allows us to derive the triboelectric material and device figure of merit. We reveal the importance of optimizing device capacitance, not only load resistance, and minimizing the impact of parasitic capacitance. Optimized capacitances can lead to an overall increase in power density of more than 10 times.


Rubbing two different materials is a charge-generating method as old as the discovery of electricity itself. Attempts to produce electricity from this triboelectric method has led to modern instruments such as the Van de Graaff generator (1), which is for sourcing an extremely high voltage. Today, with the renewed interest in energy harvesting for distributed power sources, triboelectric generators are gaining interest as an alternative source. Many variants of application schemes have been demonstrated or conceptualized recently (217) and described as triboelectric nanogenerators (18) in pioneering work, leading to deepened interest in understanding their power and efficiency (1925). Nevertheless, power optimization schemes are still premature where, often, only the load resistance is adjusted with lack of a general principle for full optimization. Furthermore, the impact of device imperfections (23) tends to be discussed in less detail despite its significant influence on device operation.

Here, we concisely model the contact-separation mode generator and reveal that the output power could be increased substantially by optimizing the capacitance in addition to the load resistance. Matching the resistance-capacitance (RC) product to the characteristic frequency of the triboelectric process is essential to use the mechanical motion and effectively convert it into electrical power. An added capacitor is necessary to stabilize the RC product because the total capacitance varies during the cycle. Adjusting resistance alone leads to poor matching conditions and, thus, low power.

We also model how parasitic capacitance, which is practically inevitable, affects the performance of the device. We show that proper adjustments to minimize the ineffectiveness caused by parasitic capacitance could also lead to a substantial increase in output power. We anticipate that, with the power improvement strategies learned from our model analysis, triboelectric generators could be a competitive energy-harvesting source.


Using the triboelectric effect

Triboelectricity is best known for forming charges on a material surface upon contact or rubbing. Although a pair of oppositely charged surfaces are essential for triboelectric power generation, the analysis described below shows that it is the motion of these electrostatic charges (not their production itself) that leads to useful electric power.

Any two different surfaces put into contact will exchange charge as determined by their Fermi level of electrons [see the study of Matsusaka et al. (26) for a review on detailed complexities]. The Fermi level relative to the vacuum level is known as the work function of a material. Materials with a high work function are more electronegative and easily accept electrons from a lower–work function (electropositive) material. Therefore, a larger difference in work function drives more charge exchange. As the charges exchange, they create an electrostatic potential difference, which tends to pull them together, thus creating a contact potential at the point of contact that exactly counteracts the work function difference.

Like a parallel plate capacitor, the electrostatic potential V = σ⋅d/ε between surface charges is the product of surface charge density σ and charge separation distance d (relative to permittivity ε). When the surfaces touch, charge separation distance is smallest and, thus, σ is largest. As the surfaces separate, the electric field increases rapidly, pushing the charges back to their original material. If the materials have sufficient conductivity, these charges flow through points still in contact during the separation process (Fig. 1). In insulators, where charges are not as free to move, some surface charge density remains after the surfaces are completely separated. As these fixed surface charges are pulled apart, mechanical work against the force of attraction is required. This work is what is converted into electrical energy observed at the load. It is the electrostatic induction of these residual surface charges that can account for all the electrical currents and voltages observed in an external circuit.

Fig. 1 Triboelectric charge generation.

Two materials with different work functions φw exchange electrons when brought into contact, generating surface charges. Rubbing is not essential as long as intimate contact is made. When the two materials separate, an electrostatic voltage is produced, which establishes a force to move the charges back if there is any conducting path. When at least one of the materials is not conductive at the surface, surface charges remain after the separation.

Device operation and its optimization

Because the interfacial charge transfer is not involved in the external circuit, we only need to consider the steady-state total surface charge σS (S is the surface area) that remains to analyze the generation cycle during the operation. Any charge loss due to leakage through the dielectric, although reducing efficiency, is replenished with the new contact at every cycle, and a steady state surface charge is established.

Triboelectric generators have a generation mechanism involving a periodic mechanical motion that separates the triboelectrically established surface charge σS and induces a change in the electric potential between electrodes. The changed electrical potential induces a flow of charge (dQ/dt) through a load resistor, which could be used as electrical energy. Because mechanical energy is converted into electrical potential energy at the same time energy is being extracted (decreases the potential), the optimization of generation is inherently related to matching the frequency or time scales of these two processes. In short, the mechanical motion should “resonate” with the electrical energy removal. It is expected from this general understanding that the optimization of a device will involve tuning circuit elements (which changes the characteristic time of charge response) to match a specific mechanical motion or mode.

The most effective way to do mechanical work by separating charged flat surfaces is to separate them in the direction perpendicular to the surface, as depicted in Fig. 2A. For this process to be effective throughout a cycling motion, the frequency of the mechanical motion (ω) should match the characteristic frequency of the circuit (1/RCtotal). However, most types of triboelectric generators (including variants of the contact-separation mode and lateral sliding mode) have a time-varying capacitance; the dielectric layer and air gap sandwiched between metals altogether acts as a capacitor, and the mechanical motion changes the distance between the two surfaces and, thus, the total capacitance. This time-varying RC product is why the generated current is, in general, not in phase with the mechanical motion and, thus, not most effective [unless the generator is designed to have a constant capacity, such as the freestanding mode (27)].

Fig. 2 Principles for optimizing triboelectric generation.

(A) Schematic of a typical triboelectric generator (contact-separation mode). The dielectric maintains a surface charge throughout the cycle. When the moving electrode pulls away from the dielectric and enlarges the air gap (left), charge separation builds a potential that drives some of the charge to the stationary electrode. A load resistor is placed between the electrodes to produce work from this charge motion. The moving electrode cycles back and forth, producing an alternating current. The mechanical motion is characterized by an angular frequency ω. (B) A device circuit scheme where the device 1/RC is maintained close to ω. The device 1/Cdevice is large, making the relative fluctuation in 1/Ctotal = 1/Cair + 1/Cdevice smaller. A large load resistance is used to compensate for the large 1/Cdevice. (C) A device circuit scheme where overall matching is poor. 1/Cdevice is small, resulting in a large fluctuation in 1/RC. (D) The charge deposited in the stationary electrode. The good matching case (blue) resembles a sinusoid, whereas the poor matching case (red) is highly distorted. (E) The overall current decreases with larger 1/Cdevice because of the requirement of a larger load resistor. (F) The time-averaged power is optimized with a balance of both good matching and current. Power is generally higher with good matching rather than with high current.

The time variance of the RC product can be reduced to better match the mechanical driving motion over the entire cycle. We compare Fig. 2B, where the device has a thick dielectric layer (or large 1/Cdevice), to Fig. 2C, where the device has a thin dielectric layer (or small 1/Cdevice). For a given mechanical motion, the change in the overall RC product is much less when 1/Cdevice is large, which leads to better matching. As a result, the charge motion with respect to time (Fig. 2D) becomes a smoother function similar to that of the driving motion, from which mechanical energy is converted more effectively. The strategy for reducing the time variance of the RC product for better matching is applicable to any type of triboelectric generator with a varying capacitance.

The trade-off for using a large 1/Cdevice for good overall matching is that one needs a larger load resistance RL to match a similar 1/RC = ω. Figure 2E shows the decrease in current due to the requirement of a larger RL for matching.

As a result, power generation is optimized when both Cdevice and RL are balanced (Fig. 2F). As one departs from the optimum, it is generally advantageous to have better matching rather than higher currents (Fig. 2F), which highlights the importance of good matching. This concept of optimizing both R and C is in contrast to the dominant literature that optimizes only the load resistance to find maximum power while keeping the device at a large Cdevice (poor matching). Cdevice can be tuned not only by changing the dielectric thickness but also by adding a series capacitor to the circuit. Using a simple model, we demonstrate how this simple two-element optimization (RL and Cdevice) can lead to output power many times higher than conventional operating conditions.

Model for triboelectric power generation

We study the most basic form of triboelectric generation, the contact-separation mode, which could be the basis for analyzing other variations. As long as the generation involves a mechanical motion, which forms a steady-state generation cycle with a characteristic frequency, the scaling behavior should be generally applicable to other types of triboelectric generators.

Device parameters that could be experimentally controlled are numerous, including the maximum separation distance, mechanical frequency, dielectric surface area, load resistance, and capacitance. These parameters form a prohibitively large phase space with interdependencies, making it difficult to devise a general strategy and define a figure-of-merit for the maximum power density of triboelectric generation.

Using a general modeling approach, the large number of device parameters can be combined into a few dimensionless quantities that guide the effectiveness of power generation. We find two, but only two, independently adjustable parameters—dimensionless resistance and capacitance—which can be optimized for the maximum power output. We show through the scaling properties of the dimensionless parameters that our results could be applied to both nanogenerators and larger-scale generators. Practical limitations causing reduced effectiveness or efficiency can be understood in terms of departure from the ideal model.

Theoretical maximum power and figure of merit

Our interest is in the elemental driving motion, which is a sinusoidal motion that changes the air gap x as a function of timeEmbedded Image(1)

We assume that the two electrodes (moving and stationary) of area S form an ideal capacitor where a charge σS gives a voltage V across the capacitor with capacitance C, where V = σS/C. Here, C is due to the dielectric layer and the air gap. Electrical power is extracted from a load resistance in series (Fig. 2).

During the initial contact (x = 0), opposite charges will be produced on the contacting surfaces. Because we are interested in the steady-state power production, we assume that a total surface charge of − Sσ remains constant on the dielectric layer surface and that it represents the steady-state surface charge. In the electrodes, a charge Q is transferred back and forth between the moving and stationary electrode (through the load resistor) due to the electrical potential change that accompanies the change in x, eventually reaching a periodic oscillation in steady state. In this ideal system, charge is not lost or dissipated, and triboelectric charge generation need not be considered at steady state.

Combining these principles results in a dimensionless differential equation describing the charge in the stationary electrode (see the Supplementary Materials)Embedded Image(2)

Here, θ = ωt is the phase of the driving motion; Q* = Q/Sσ is the dimensionless charge in the stationary electrode, which is a fraction of the total surface charge Sσ; and x* = x/xmax is the fractional air gap that sinusoidally varies between 0 and 1. The dimensionless capacitance C* and load resistance R* are the only two independent variables that characterize the systemEmbedded Image(3)Embedded Image(4)where RL is the load resistance, ε0 is the permittivity of vacuum, Cdevice is the capacitance of the entire device when the air gap is closed (that is, includes the dielectric layer and also any additional capacitor connected in series), and Cair = ε0S/xmax is the capacitance of the maximum air gap. The governing equation is identical when the moving electrode surface has an additional dielectric layer (19).

The power output is characterized by the averaged dimensionless powerEmbedded Image(5)where I* = dQ*/dθ.

By solving for the steady-state solution of Eq. 2 (see the Supplementary Materials), we find how the power output is optimized with respect to R* and C* (Fig. 3). It is seen that the maximum condition for Embedded Image is pinpointed at one particular conditionEmbedded Image(6)

Fig. 3 Analysis of the ideal model in steady state with dimensionless parameters.

(A) Power map for the two device parameters R* and 1/C*. The color represents the dimensionless power per cycle Embedded Image. Restricted optimization of power for a given 1/C* (solid line) and a given R* (dashed line) is shown, and the two coincide at the global optimum condition (green dot). Two example suboptimal conditions (red dot, poor matching case; blue dot, good matching case) are indicated. (B) V* (= I* R*), (C) I*, and (D) P* curves for the global optimum condition and two suboptimal conditions indicated in (A). It is seen that the poor matching condition (red lines) produces highly distorted curves with low power, whereas the good matching condition (blue lines) shows characteristics as good as the global optimum condition (green lines).

This condition, together with Eqs. 3 and 4, specifies how the real physical parameters should be optimized. For instance, for a lower-frequency driving motion (small ω), one would need either a larger load resistance or a larger device area; if two different sets of physical parameters yield the same R* and C*, they specify identical device operation conditions.

Real power density (time-averaged) can be calculated by scaling Embedded ImageEmbedded Image(7)which offers the definition for a “device figure of merit” for power density in contact-separation mode by setting Embedded Image to its maximum valueEmbedded Image(8)

FOMdevice is the maximum power density obtainable from the contact-separation generator. Note that Embedded Image is the average speed of the mechanical motion. It is seen that the only material-related parameter in the FOMdevice is the steady-state surface charge density, which leads to the definition of the “material figure of merit” for power densityEmbedded Image(9)

To keep the number more tractable, one could also use a normalized material figure of meritEmbedded Image(10)

Note that FOMmater is a surface property dependent on which material it is brought into contact with. FOMmater coincides to what has been established in the study of Zi et al. (21).

We next show an example of power density estimation using typical device parameters. We assume a generator with σ = 12 μC/m2, xmax = 1 mm, and w = 2π × 60 Hz. Then, the maximum power density obtainable with a fully optimized generator is FOMdevice = 124 mW/m2. In a conventional device scheme where only the dielectric layer of thickness d is used without an additional series capacitor, 1/C* = drxmaxr is relative permittivity). In this case, 1/C* is typically < 0.05 (that is, d < 150 μm for εr = 3), at which the maximum power density is only 23 mW/m2, more than five times lower than the FOMdevice.

The above example illustrates how the device capacitance has a significant impact on the output power, an aspect that has been mostly overlooked in the recent literature. Increasing 1/C* toward the global optimum can be done not only by using a thicker dielectric layer but also by adding a capacitor in series to the dielectric layer: 1/Cdevice = 1/Cdielectric + 1/Cseries.

Partly similar conclusions regarding the relevant parameters have been reached in the literature (for example, regarding σ2 or Embedded Image) (19, 21); however, earlier reports have based their model analysis on the transient behavior (in the transient state, the power output could be instantaneously high at nonoptimal conditions; see the Supplementary Materials); the literature tends to focus on a restricted optimization only based on load resistance (or sometimes no load conditions), leading to disparate conclusions regarding power output and optimization.

Zi et al. (21) offered one of the first dedicated discussions on the figure of merit for triboelectric nanogenerators by considering an upper bound for power output constructed with a combined cycle of open-circuit (potential buildup) and short-circuit (instantaneous discharge) processes. Open-circuit voltages and short-circuit charges were used to define a figure of merit. However, from the perspective of designing a device, it is more convenient to have a figure of merit that could be estimated just with the design parameters rather than with output characteristics. Furthermore, the practical output limit at a load resistor is much lower than what can be calculated from the instantaneous discharge process. The FOMdevice we propose is the actual maximum power density one could obtain at the load resistor under optimized capacitance and no parasitic reduction.

Power reduction due to parasitic capacitance

An important factor causing a deviation from the ideal case will be the presence of parasitic capacitance in the circuit. It is a small but unavoidable element that always exists in an electric circuit. In triboelectric generators, the parasitic capacitance can easily be comparable to that of the generator, significantly affecting the power output and the optimum device parameters (23). This leakage mechanism can be modeled with an additional capacitor in parallel to the device (Fig. 4). Then, instead of the charge produced from the triboelectric effect always going through the load resistor, part of it leaks into Cpar. This leakage affects the effectiveness (but not efficiency) of generating power, slowing down the energy conversion rate.

Fig. 4 The influence of parasitic capacitance.

(A) Decrease in the maximum obtainable power with increasing dimensionless parasitic capacitance Embedded Image. (B) Shift in the optimum R* (solid line) and 1/C* (dashed line) condition with the presence of Embedded Image.

Parasitic leakage is characterized with the dimensionless parasitic capacitanceEmbedded Image(11)where Cpar is the real parasitic capacitance.

By solving for the steady-state solution (see the Supplementary Materials), we find how the output power is reduced with increasing parasitic capacitance (Fig. 4A). For example, assume xmax = 1 mm and S = 22 cm2, yielding Cair = 20 pF. Then, a parasitic capacitance of 20 pF indicates Embedded Image, which reduces the maximum output power to 54 % of the FOMdevice. The optimum R* and 1/C* also downshift with Cpar (Fig. 4B).

Increasing the area of the metal-dielectric interface is the best way to minimize the impact of a given parasitic capacitance. To decrease Embedded Image for a given Cpar, one can either increase S or decrease xmax(Eq. 11); however, because xmax is directly proportional to the device figure of merit (Eq. 8), only increasing S is beneficial. Therefore, larger-area generators are much more robust to parasitic leakage.


Comparison to experiments

Our model shows good agreement with the experimental data reported by Niu et al. (19). We compare the power generation curve with respect to load resistance. For the data set where the device capacitance was not varied, the corresponding model curve is a vertical cross section of Fig. 3 (also see fig. S3). It is seen in Fig. 5A that the curve shapes from three different driving speeds all fit well to a fixed curve shape given by 1/C* = 0.07.

Fig. 5 Model fitting to experimental data from Niu et al. (19).

(A) Power change with respect to load resistance. Experimental data sets consist of driving conditions with average speeds of 0.08 m/s (black circles), 0.04 m/s (blue circles), and 0.02 m/s (red circles). Solid lines are from the ideal model by assuming a fixed 1/C* = 0.07, which determines the shape of the curve. The scaling of the model curve to real dimensions was determined by fitting due to the lack of experimental details in the study of Niu et al. (19). (B) The change of maximum power (left axis) and optimum resistance (for a fixed device capacitance; right axis) with respect to the average driving speed. The curves represent the expected scaling behavior of maximum power (Embedded Image, Eq. 7) and optimum load resistance (∝ 1/ω for fixed xmax;, Eq. 3) from the ideal model.

Because the maximum power conditions for the three different driving speeds in the data set are expected to be described with identical dimensionless parameters (that is, a set of 1/C* and R*), the maximum power and optimal load resistance in real dimensions are expected to scale with respect to average speed as compared in Fig. 5B. The small deviations are thought to be associated with parasitic capacitance, which was not possible to be accounted for with the given information.

Parasitic capacitance seems to explain the power scaling observed with respect to an increase in xmax. Niu et al. (19) increased xmax by five times, which increases the FOMdevice by the same factor. However, the observed power increase was limited to a factor of 2.3. Note that for a given parasitic capacitance, Embedded Image scales with xmax (Eq. 11); that is, the influence of parasitic capacitance increases. As a result, the benefit of increased xmax is diminished.

Regarding the material figure of merit, several experimental reports have demonstrated the significance of σ2 (2831). For quantitative assessment of the scaling behavior with respect to σ2, power output at the load resistance should be measured while taking into account the 1/C* of the generator.

Conversion efficiency

Because the power production in a triboelectric generator is purely by electrostatic induction, the energy conversion efficiency of this process itself is 100%. It can be confirmed (see the Supplementary Materials) with our model that the power output is exactly equal to the net work done by the mechanical motion, which is true, in general, even for nonoptimal R* and 1/C* and a finite Embedded Image.

Loss of efficiency originates from other imperfections rather than from the electrostatic conversion process or purely capacitive leakage. For example, the surface charge on the dielectric layer could leak through the film, which is equivalent to resistive loss. This type of loss can be minimized by keeping the dielectric film resistance much larger than the load resistor. Frictional losses regarding the mechanical motion or the charge exchange process (during the contact of two surfaces) are other sources of losses that could be important.

Practical requirements in implementing the generator could also lead to reduced efficiency. The 100% conversion efficiency of the electrostatic induction is by considering the net mechanical work, including both positive (during separation of charged surfaces) and negative (while joining the charged surfaces back together) work input. In applications where negative work could not be meaningful, the relevant efficiency would have to be calculated with only the positive work input. This alternative definition gives 45 % efficiency for the maximum power condition (see the Supplementary Materials). Another example is when the load is not simply a load resistor. Then, there would be losses associated with power conditioning (32) to convert the output into a “useful” form.


Triboelectric generators make use of the surface charge maintained in an insulating material during the repeated contact and separation of a dissimilar material. Electrostatic induction due to the separation of two charged surfaces builds a potential that drives a current through the load.

Power generation is more effective when the time-varying RC product of the generator better matches the mechanical motion frequency during the cycle. A larger 1/Cdevice is beneficial for good matching but reduces the current level. Generally, the power output characteristic is more stable in the good matching regime (high 1/Cdevice), which could be achieved by connecting a series capacitor to the circuit.

The device figure of merit (Eq. 8) determines the maximum output power density obtainable from a generator at steady state. Here, the only material-related parameter is surface charge density, which defines the material figure of merit.

Parasitic capacitance reduces the power generation effectiveness. Enlarging the surface area of the dielectric-metal interface can minimize power reduction due to parasitic capacitances. The fundamentals of triboelectric generation presented in this paper could lead to designs of more sophisticated applications and help to establish triboelectric generators as an attractive energy-harvesting source.


Model differential equations were found by constructing circuit diagrams equivalent to the generator device. The governing equation was reduced to have a minimum number of independent variables equal to the degrees of freedom. Steady-state solutions were found analytically, which were numerically evaluated for model analysis. Details can be found in the Supplementary Materials.


Supplementary material for this article is available at

fig. S1. Circuit model diagrams of triboelectric generators.

fig. S2. Transient characteristics of the ideal model.

fig. S3. One-dimensional projections of the steady-state dimensionless power output.

fig. S4. The influence of parasitic capacitance on device characteristics.

fig. S5. Comparison of the output power and mechanical work input.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: Funding: This work was supported by the Solid State Solar Thermal Energy Conversion Center, an Energy Frontier Research Center, funded by the U.S. Department of Energy, Office of Science, Basic Energy Sciences (DE-SC0001299). Author contributions: J.P. performed the background experiments, which inspired the project, and conducted model analysis. S.D.K. constructed the dimensionless models and wrote the manuscript. G.J.S. supervised the project. All authors discussed the results throughout the project and edited the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Stay Connected to Science Advances

Navigate This Article