Rectifying the output of vibrational piezoelectric energy harvester using quantum dots.

Piezoelectric energy harvester scavenges mechanical vibrations and generates electricity. Researchers have strived to optimize the electromechanical structures and to design necessary external power management circuits, aiming to deliver high power and rectified outputs ready for serving as batteries. Complex deformation of the mechanical structure results in charges with opposite polarities appearing on same surface, leading to current loss in the attached metal electrode. External power management circuits such as rectifiers comprise diodes that consume power and have undesirable forward bias. To address the above issues, we devise a novel integrated piezoelectric energy harvesting device that is structured by stacking a layer of quantum dots (QDs) and a layer of piezoelectric material. We find that the QD can rectify electrical charges generated from the piezoelectric material because of its adaptable conductance to the electrochemical potentials of both sides of the QDs layer, so that electrical current causing energy loss on the same surface of the piezoelectric material can be minimized. The QDs layer has the potential to replace external rectification circuits providing a much more compact and less power-consumption solution.

Piezoelectric effect is exhibited in crystals with no centre of symmetry. Under an external stress, a net polarization appears on the surface of the material, generating piezoelectric voltage, inducing electrical charges. Mechanical-to-electrical energy harvesting devices based on piezoelectric materials have been receiving much attention in recent decades1234. Moreover it has been successfully advanced to nanoscale energy devices recently56. Note that charge can only be continuously generated when there is a change in the applied force, which means that a static force only generates a fixed amount charges that will leak away quickly across the external circuit (internal screening effect is neglected). The time taken for leaking away all generated charges depends on the input impedance of the external circuit. For instance a 1 kΩ load resistance R and a 1 nF internal capacitance C of the piezoelectric transducer results in a time constant of R × C = 1 μs, which implies that after 1 μs all generated charges from a static mechanical force disappear. For this reason the device works mainly under AC conditions for energy harvesting purpose, in order to continuously generate electrical charges. This brought out the concept of ‘vibrational’ energy harvesting structures.

The optimal power generation is when the device operates at its resonant frequency7, which leads to essential additions of external rectifying circuits8910, causing unavoidable power consumptions. The vibrational energy harvesters usually appear in the form of a cantilever with one end fixed, the other suspended11121314, as mechanical cantilevers are able to generate larger strain, easy to fabricate, and having lower fundamental frequencies to match with ambient mechanical sources. Although tremendous efforts have been spent to increase the output power of the cantilever harvester, including some innovative structural designs151617. There is an ongoing issue of designing low power consumption AC-DC converters for rectifying the AC output of the harvester181920. Complex vibrations of the cantilever subjecting to arbitrary excitations from ambient sources sometimes result in both positive and negative charges accumulating on the same surface, subsequently introducing electrical current on the metal electrode attaching to that piezoelectric surface. This amount of power loss is not compensatable using the external circuits. We used segmentation protocol aiming for collecting positive and negative charges separately212223, essentially creating many sub-cells. However it will require more complex external circuits to handle these sub-cells, and the consequence will be more electrical power dissipated in the power management circuits.

Inspired from the electron quantum transport theory of conductors24, we propose to coat a layer of QDs onto the piezoelectric cantilever (schematic diagram is shown in Fig. 1a). The QD layer is able to function as an integrated filter or rectifier that has adaptable resistance under different chemical potentials of both sides of the QDs layer, which are the metal electrode and piezoelectric material respectively. Thus charges with different polarities and quantities cause variation of electrochemical potential of the piezoelectric surface attaching to the QDs layer, adjusting the resistance of the QDs layer. The theory and method underpinning this hypothesis is detailed in the next section, followed by the numerical case study. This new method is much more efficient than previous approaches of using external circuits or segmenting electrodes, since it does not need additional circuits consisting of active electronic components like transistors, diodes. This design seamlessly integrates power management circuits with the actual transducer, moreover it significantly reduces power loss on the same surface of the piezoelectric material in the cases of the harvester experiencing complex deformations.

Figure 1

Schematic view of the integrated energy harvester.

(a) 3D view of the stacked QD and piezoelectric layers. Circled diagrams show the principle of adaptable conductance of the QD subjecting to different surface charges, where red dot represents low conductance, green dot shows high conductance. (b) Energy diagram of the electrode-QD-piezoelectric interface. (c) equivalent circuit of the device.

Models and Methods

The multiphysical treatment for this new concept is split into three steps, namely mechanical-to-electrical (MTE), charge-to-electrochemical potential (CTE), and elastic tunnelling (ET) through the quantum dot. Euler–Bernoulli beam theory together with the constitutive piezoelectric equation will be used in the MTE stage, CTE simulation is based on semiconductor physics, and ET simulation is built on the quantum mechanics theory. This open-loop multi-physics simulation links the external applied mechanical excitation to the conductance of the QD, through which the rectifying function of this concept is demonstrated. Finally this integrated device is connected to a load resistor for demonstration in the real applications.

Mechanical strain to electrical charge

To begin with the analysis, we describe the integrated energy harvester as the form of a cantilever, which consists of a layer of piezoelectric material acting as the strain to electrical charge converter, and a layer of quantum dots that can be spin coated or printed. In the real scenario, metal electrodes are plated on both sides of the composite cantilever. Due to very thin electrode layer, its mechanical contribution is neglected in the analysis. The length L, width w, and thickness g of the piezoelectric layer are designated as 5 mm, 0.5 mm, and 50 μm for the later numerical demonstration. AlN is chosen as the piezoelectric material with density ρ of 3300 kg/m3 25, Poisson’s ratio ν and Young’s modulus Y of 350 GPa26. The quantum dots layer is so thin (normally <1 μm), and it can also be neglected in the mechanical analysis. The dynamic motion of the cantilever can be described using the mass-spring-damper equation27

where m, b, k1, k2 are total effective mass (here we attach an additional proof mass (100 × Lwgρ) at the free end of the cantilever to reduce the resonant frequency), damping coefficient, linear stiffness, and nonlinear stiffness respectively. z(t) is the tip displacement. The first derivative of the z with respect to the time t represents the velocity, and the second derivative denotes the acceleration. A is the amplitude of the external periodic driving force, and ω is the driving frequency. When ω = ω0 (ω0 being the natural frequency of the cantilever , here we only analyse the fundamental resonance), the deformation of the cantilever reaches to the maximum. Here we set the driving frequency same as the nature frequency of the cantilever to ensure an optimal match between the excitation and the beam. Under the external driving, the cantilever vibrates, the deflected cantilever can be treated as an arc. The curvature 1/r (r is the radius) of this arc is expressed with respect to the tip deflection z

The axial strain on the surface of the piezoelectric layer along the x-axis can be obtained in terms of curvature and tip deflection27

The axial stress σ in the piezoelectric layer is derived through the stress-strain relation, that is σ(t) = Yε(t). The piezoelectric constitutive equations relate four field variables stress components Σ, strain components S, electric field components E, and the electric displacement components D, which can be described as 2829. We shall focus on the generated electric displacement D from the mechanical stress Σ. Assuming no external electric field is applied to the cantilever, D = dΣ, d being the piezoelectric constant. The electrical charges generated from a tip deflection z(t) can therefore be obtained

where d31 is the piezoelectric constant, relating the strain in the ‘1’ direction (x-direction) to a generated field along the ‘3’ direction (z-direction). In the simulation, we take d31 = −1.9 pC/N30. A upward deflection (z > 0) corresponds to a compressive stress on the top surface, therefore inducing negative charges.

Charge to electrochemical potential

After obtaining piezoelectrically generated charges, electrochemical potentials (top and bottom surface) of the piezoelectric layer can be derived from the carrier density equation31, which is

where is the density of state for 3-dimensional bulk materials, which takes the form of  32, where ħ is the reduced Planck’s constant (1.05 × 10−34 Js), m* is the effective mass of electron m* = 0.33 m033, where m0 is the bare electron mass 9.11 × 10−31 Kg. n is the carrier density, which is varying upon deflection of the piezoelectric cantilever. The additional carrier density is calculated by Δn = Q(t)/(wLd), d denotes the depletion thickness on the surface of the piezoelectric, which can be approximated according to the Poisson equation . Here for simplicity it is taken as a constant 10 nm. is the Fermi distribution of electrons , where is the electron energy, μ stands for the chemical potential, k is the Boltzmann constant (8.62 × 10−5 eV/K−1), and T denotes the temperature. Iteration method/self-consistent method is used to obtain μ through the numerical approach34. In the work we calculate the Fermi distributions of electrons in conduction band f and holes in valence band f, where the electrochemical potentials μ and μ are electrochemical potentials for electrons and holes respectively. The effective mass of heavy hole is taken as 3.68 m033. Focusing on the interface between QD and piezoelectric material, accumulated electrons due to piezoelectric effect causes electrochemical potential move to the conduction band minimum μ. Accumulated holes at the interface tend to shift the chemical potential to the valence band maximum μ. The middle of the bandgap g is set to potential 0, μ = g/2, μ = −g/2. To solve the equation (5), the integral is calculated from conduction band minimum to a large value of or from to the valence band maximum, is large enough to arrive at an accurate result. Calculated chemical potential μ and μ are then mapped to the energy graph considering the bandgap g. For AlN, g = 6.2 eV35. The mapped electrochemical potential has been designated as the chemical potential of the piezoelectric interface μ = [μ + g/2, (n > n)] ∪ [−μ − g/2, (n < n)]. Figure 1b schematically displays the energy graph of the metal-QD-piezoelectric material interface. It is assumed that the piezoelectric material is intrinsic with the initial chemical potential in the middle of the bandgap, i.e. internal screening36 of piezoelectric polarization is minimal. This assumption is reasonable as the AlN films exhibit very low carrier concentration in the real case.

Elastic tunnelling

In the normal operation conditions of the integrated piezoelectric vibrational energy harvesting device, surrounding temperature is usually room temperature. The energy cost to add an electron to the dot is determined by the dot’s self-capacitance, that is e2/CΣ. Simple estimation of the self-capacitance is CΣ = e2/(8εε0R)37 where R is the radius of a 2D dot. For a typical GaAs dot with radius of 50 nm, CΣ is calculated to be 0.046 fF, resulting in e2/CΣ = 3.5 meV, much smaller than thermal energy at room temperature 300 K, kT = 25.9 meV. Hence Coulomb blockade effect can be ignored in this analysis. After μ is obtained, electrical current I flowing through the quantum dot is derived using the Landauer formula2438, that is

where f and f are electron Fermi distributions for the top electrode (note the QD layer is deposited on the top surface) and piezoelectric layer respectively, , . In the analysis, temperature T remains constant 300 K. Here we shall calculate the μ according to the previous section, and μ is assumed to be 0, aligning with the chemical potential of the initial state of the device. Without external mechanical excitation, the initial chemical potential of the piezoelectric layer is 0. represents the transmission function of the two interfaces for each incident electron energy . It is governed by the Lorentzian equation

where the γ represents the elastic tunnel coupling coefficients through electrode-quantum dot, and quantum dot-piezoelectric layer respectively. They are determined by 39, here δz denotes the distance between quantum dot edge and the electrode, quantum dot edge and piezoelectric layer. ξ is the characteristic length of the dot defining the exponential decay for the probability of finding electron outside the dot, i.e. decay of the wavefunction of the electron, which can normally be calculated from theory of electron passing barriers. Here in this analysis it is assumed that γ0 = kT. For the simulation we set δz = 20 nm, ξ = 2 nm, then we have γ = 0.45 × 10−4 kT. It is worth noted that γ can be adjusted by changing positions of the QDs to achieve various rectification effects. The energy level of the quantum dot is set to 1.1 eV (Here the fluctuation of the due to electrons flowing through the dot is not considered). Here we design the in-plane spacing between dots to be much larger than ξ, so that the current flow between dots can be neglected. The conductance of the quantum dots layer is

where the integral is calculated for the energy range of (−200 kT, 200 kT), which is sufficient to cover the varying potential of the piezoelectric material. Once G is obtained, further analysis on the integrated harvester connecting with an external load is discussed. The complete circuit consists of a load resistor and the integrated QD-piezoelectric energy harvester (Fig. 1c). It is assumed that the piezoelectric material has infinite resistance i.e. a pure capacitor. As the cantilever is vibrating under the external load, conductance/resistance of the quantum dots layer changes with the varying μ, only allowing electrons having energy within a certain band passing through. The next section demonstrates numerical analysis of the integrated harvester.

Results and Discussion

Mechanical analysis

A case study to quantitatively demonstrate this idea is detailed in this section. The numerical procedure begins with the mechanical analysis. In solving the equation (1), the amplitude of periodical driving force A is set to ma where m is the sum of the cantilever mass and the proof mass, a denotes the externally applied acceleration (10 m/s2). The linear stiffness of the cantilever can be calculated from k1 = (Ywg3)/(4L3)40, nonlinear stiffness k2 = 2k1, and the damping factor is usually very small for miniaturized devices. Here we use b = 2(k1/m)1/2 × 10−5. From the numerical solution of the equation (1), substitute tip deflection z(t) into the equations (2, 3, 4), generated charge Q(t) is calculated with respect to time t. Figure 2 depicts the result of z(t) and Q(t). It is seen that the cantilever oscillates periodically with a maximum peak displacement of 1.1 μm, generated charge oscillating in a similar manner with the peak amplitude of 3.66 × 10−12 C.

Figure 2

Vibration of the cantilever inducing electrical charges, frequency is calculated to be 335.9 Hz.

(a) Tip deflection vs. time. (b) Generated charge vs. time t.

Chemical potential

Next the resulting charge values are input to the equation (5) to arrive at chemical potentials μ and μ, which are mapped to the band diagram and shown together with μ and μ in Fig. 3. It is seen from Fig. 3a that in the low density range, a small addition of the carrier density causes chemical potential to deviate to μ or μ rapidly. But the chemical potential remains steady in the large density range. That is due to the exponential nature of the Fermi-Dirac distribution of electrons in both conduction and valence bands. As the cantilever vibrates, the chemical potential on the top side of the piezoelectric layer (attaching to the QDs) μ with respect to time t is shown in the Fig. 3b. Compared with the results in Fig. 2b, μ curve displays a square wave instead of a sinusoidal, which coincides with the results of μ and μ versus carrier density.

Figure 3

Chemical potentials calculated from the increase of carrier density and polarity change of generated charges on the piezoelectric surface due to vibration of the piezoelectric cantilever.

(a) μ, μ, μ, and μ in relation to the carrier density. (b) μ vs. time t.

Analysis of conductance

Further on, the electrochemical potential μ is used to calculate the equation (8) for the purpose of obtaining conductance as the function of time t. Conductance G(t) is then obtained numerically for a single quantum dot. Suppose dots density is 2.78 × 1014 m−2, total conductance of the whole surface area of the cantilever G can be considered to be the sum of many single conductances G0 connected in parallel. Then it is G = ∑ G0 = 6.94 × 108 G0. The results are shown in the Fig. 4, where both R(t) (resistance) and conductance G(t) curves are displayed. Conductance varies from a very low value (1.74 × 10−8 S) to a high value (0.07 S), corresponding to resistance from 5.74 × 107 Ω to 14.5 Ω. A simple circuit consisting of the integrated device and a 1 kΩ resistor is sketched to demonstrate the rectification effect. From Fig. 5, clearly the rectification effect has been achieved and almost all positive part of the voltages has been blocked due to very large impedance of the QDs layer, attributed to its energy selection mechanism, where V1 and V2 refer to the equivalent circuit in Fig. 1c. Graphical demonstration of conductance vs. chemical potential (G-μ) curve with the varying δz for a single QD is shown in Fig. 6, which clearly displays a performance similar to that of diodes. As the current flows through the close-loop circuit consisting of a finite external load resistance, the electrochemical potential μ shifts to the neutral. Results from equation (5) (Fig. 3b) reveals that the variation of the μ is mainly determined by the polarity of the accumulated charges in the QD-Piezoelectric interface. A small quantity of charges will bring the μ close to μ or μ, and further increasing the quantity does not vary μ significantly (Fig. 3a). This makes the metal-QD-piezoelectric interface as a perfect rectifier compared to the classic diodes that require a forward bias.

Figure 4

Calculated conductance and resistance of the QD layer with respect to time.

(a) Conductance vs. t. (b) Resistance vs. t.

Figure 5

Simulation of the rectification effect of the integrated piezoelectric harvester.

(a) Voltage with respect to the time for piezoelectric layer. (b) Voltage vs. t for the integrated device.

Figure 6

Conductance vs. chemical potential of piezoelectric layer.

Complex mechanical vibration

Complex motions of the mechanical structure give charges with opposite signs on the the same surface, for example a 3rd resonant mode displayed in Fig. 7. The mode shape is obtained from Euler–Bernoulli beam theory on cantilever resonance41, the 3rd resonant mode is found by numerically solving the third root β3 of the nonlinear equation cosh(βL)cos(βL) + 1 = 0. d2 is for the curvature at x point of the cantilever. Subsequently the mode shape and charge distribution can be calculated. It is seen that both positive and negative charges are generated on the surface, causing charges cancellation. With the QDs layer, only one type of charges can flow through, reducing the current loss on electrode. Though this hypothesis is mainly based on the theoretical validation, experiments in previous published work42 on thermoelectric energy devices showed that using QDs in power generation devices is viable.

Figure 7

Third resonant mode shape of the cantilever and charge distribution.

(a) 3D view of the mode shape. (b) Charge distribution. Without QDs layer, the net value of charges is −2.6 × 10−11 C. With QDs layer, the net value of charges is −5.8 × 10−11 C. Cantilever in static and charge in neutral are also shown in the figure (semi-transparent plane).

Conclusion

An integrated piezoelectric energy harvesting concept comprising a quantum dots layer and a piezoelectric layer stacked together is reported in this article. The new structure enables rectification function to be achieved without additional external circuits, as well as avoidance of the screening effect on the surface of the piezoelectric material. Multiphysical simulations taking account of the mechanical-electrical, charge-chemical potential, and elastic tunnelling process validate the concept. The modelling results show clearly that the rectification effect has been achieved. It is worth mentioning that this configuration can be applied to other structures besides the cantilever shape.

Additional Information

How to cite this article: Li, L. Rectifying the output of vibrational piezoelectric energy harvester using quantum dots. Sci. Rep. 7, 44859; doi: 10.1038/srep44859 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.