## Abstract

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.

## INTRODUCTION

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 (*2*–*17*) and described as triboelectric nanogenerators (*18*) in pioneering work, leading to deepened interest in understanding their power and efficiency (*19*–*25*). 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.

## RESULTS

### 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.

### 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/*RC*_{total}). 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*)].

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/*C*_{device}), to Fig. 2C, where the device has a thin dielectric layer (or small 1/*C*_{device}). For a given mechanical motion, the change in the overall *RC* product is much less when 1/*C*_{device} 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/*C*_{device} for good overall matching is that one needs a larger load resistance *R*_{L} to match a similar 1/*RC* = ω. Figure 2E shows the decrease in current due to the requirement of a larger *R*_{L} for matching.

As a result, power generation is optimized when both *C*_{device} and *R*_{L} 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 *C*_{device} (poor matching). *C*_{device} 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 (*R*_{L} and *C*_{device}) 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 time(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)(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*/*x*_{max} 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 system(3)(4)where *R*_{L} is the load resistance, ε_{0} is the permittivity of vacuum, *C*_{device} 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 *C*_{air} = ε_{0}*S*/*x*_{max} 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 power(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 is pinpointed at one particular condition(6)

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 (7)which offers the definition for a “device figure of merit” for power density in contact-separation mode by setting to its maximum value(8)

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

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

Note that FOM_{mater} is a surface property dependent on which material it is brought into contact with. FOM_{mater} 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/m^{2}, *x*_{max} = 1 mm, and *w* = 2π × 60 Hz. Then, the maximum power density obtainable with a fully optimized generator is FOM_{device} = 124 mW/m^{2}. In a conventional device scheme where only the dielectric layer of thickness *d* is used without an additional series capacitor, 1/*C** = *d*/ε_{r}*x*_{max} (ε_{r} 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/m^{2}, more than five times lower than the FOM_{device}.

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/*C*_{device} = 1/*C*_{dielectric} + 1/*C*_{series}.

Partly similar conclusions regarding the relevant parameters have been reached in the literature (for example, regarding σ^{2} or ) (*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 FOM_{device} 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 *C*_{par}. This leakage affects the effectiveness (but not efficiency) of generating power, slowing down the energy conversion rate.

Parasitic leakage is characterized with the dimensionless parasitic capacitance(11)where *C*_{par} 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 *x*_{max} = 1 mm and *S* = 22 cm^{2}, yielding *C*_{air} = 20 pF. Then, a parasitic capacitance of 20 pF indicates , which reduces the maximum output power to 54 % of the FOM_{device}. The optimum *R** and 1/*C** also downshift with *C*_{par} (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 for a given *C*_{par}, one can either increase *S* or decrease *x*_{max}(Eq. 11); however, because *x*_{max} 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.

## DISCUSSION

### 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.

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 *x*_{max}. Niu *et al*. (*19*) increased *x*_{max} by five times, which increases the FOM_{device} by the same factor. However, the observed power increase was limited to a factor of 2.3. Note that for a given parasitic capacitance, scales with *x*_{max} (Eq. 11); that is, the influence of parasitic capacitance increases. As a result, the benefit of increased *x*_{max} is diminished.

Regarding the material figure of merit, several experimental reports have demonstrated the significance of σ^{2} (*28*–*31*). 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 .

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.

## CONCLUSIONS

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/*C*_{device} 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/*C*_{device}), 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.

## METHODS

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 MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/3/12/eaap8576/DC1

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.

## REFERENCES AND NOTES

**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.

- Copyright © 2017 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).