Open Access Article
Bochao Xie
*,
Yingying Ma and
Yihuang Xie
School of Engineering & Applied Science, Yale University, New Haven 06250, USA. E-mail: bochao.xie@yale.edu
First published on 13th February 2026
The central open question in triboelectric nanogenerators is not how much power they can deliver, but what dynamic principle governs the conversion from motion to charge. We answer this by establishing a similarity modulation mechanism that extends the classical V–Q–x relation into a predictive dynamic law. A single similarity descriptor Sm = ω0T organizes the electromechanical response across device geometries and operating modes under a clearly stated condition where dielectric relaxation and charge decay are slow relative to motion. Within this conditional similarity, normalized outputs collapse onto universal curves that reveal three regimes with simple design rules: a low frequency growth region, a resonance near Sm ≈ 2π that captures impedance alignment, and a high Sm ceiling where average power becomes frequency independent due to incomplete charge refresh per cycle. The framework also ranks input waveforms and shows why compact top hat profiles outperform sinusoids at fixed amplitude, and it provides explicit guidance for load matching and storage-oriented operation without introducing new empirical parameters. This theory converts scattered models into a single map that explains published trends, exposes real limits, and supplies a compact rule set for rational TENG design and scalable energy management.
New conceptsWe demonstrate a similarity modulation mechanism (SMM) showing that triboelectric nanogenerators can be governed by one universal dynamic similarity control factor that compares how fast the device is mechanically driven with how fast charges and polarization can relax. Using this concept, we identify clear operating regimes that are easy to understand and apply: a regime where increasing driving speed boosts output, an optimum regime where energy transfer is most efficient, and a saturation regime where further speeding up no longer increases average power because interfacial charge cannot fully refresh within each cycle. This work adds new insight to nanoscience and nanotechnology by exposing the fundamental dynamic limit of contact electrification systems and enabling rational engineering of excitation waveforms, load matching, and power management strategies for scalable triboelectric energy harvesting and sensing. |
The dynamic behavior of TENGs defines the bridge between mechanical excitation and electrical response, determining not only the instantaneous performance but also the scaling limits and design logic of the entire technology. In principle, this dynamic coupling is encapsulated by the governing V–Q–x relationship, which links voltage (V), transferred charge (Q), and displacement (x) through a capacitor-like framework.11 Originating from the contact-separation model, this relationship has been extended to other operation modes12 and to structured systems such as grating-based generators.13,14 Later works constructed first-principles derivations from electrostatic or Maxwell formulations,15–22 reaffirming that triboelectric conversion is not an empirical phenomenon, but a self-consistent electromechanical process.
However, the V–Q–x formalism, despite its elegance, captures only the static or quasi-steady coupling among these variables. It does not yet provide a unified dynamic theory capable of describing how external excitations, structural parameters, and circuit conditions collectively shape system evolution. As a result, existing models often treat motion input, device structure, and load impedance as independent variables—an assumption that oversimplifies the inherently coupled nature of TENG operation.23–30 The absence of such a dynamic unification has constrained the field to empirical optimization rather than predictive design.
In this work, we address this missing foundation by introducing a similar theory based analytical framework that reformulates TENG dynamics into a dimensionless, self-scaling representation. This framework reveals that the entire electromechanical behavior of a TENG can be governed by a single similarity descriptor that couples the mechanical timescale of motion and the electrical relaxation of the circuit. Within this formalism, diverse structural modes and excitation forms collapse into a unified dynamic landscape, allowing direct interpretation of power scaling, motion regimes, and saturation behavior from a common theoretical perspective (Fig. 1).
= t/T,
= Q/q, and Ĉ = C/C0, where T represents the characteristic time scale, q denotes the characteristic charge, and C0 signifies the characteristic capacitance, we obtain the normalized governing equation:Here, the dimensionless functions are defined as ĝ−1 = C/C0 and ĵ = QSC/q, while the crucial similarity parameter emerges as Sm = ω0T = T/(RC0). The theoretical framework reveals that Sm serves as the fundamental similarity criterion governing TENG dynamics. While the structural functions ĝ = ĝ(
) and ĵ = ĵ(
) depend exclusively on dimensionless displacement
= x/D and define the specific geometric characteristics of each TENG configuration, Sm remains independent of displacement or time variables. This leads to our central conclusion: two structurally analogous TENGs experiencing kinematically similar motion will demonstrate completely scalable performance, provided they share identical Sm values. The physical significance of ω0 = 1/(RC0) corresponds to the characteristic frequency of the equivalent RC circuit formed by the initial TENG capacitance and external load, thereby providing concrete physical meaning to the similarity parameter. The general solution formalism for the system dynamics can be expressed as:
I = ω0qĝ(ĵ − ) |
This theoretical construct provides a comprehensive framework where analysing any specific TENG design simply requires incorporating its structural relations ĝ(
) and ĵ(
). For a lateral-sliding mode TENG, the structural characteristic relations adopt particularly elegant forms. The capacitance follows C = ε0W(L − x)/δ, while the short-circuit charge increases linearly as QSC = σxW. Here, the effective thickness is defined as δ = ε0d1/ε1 + ε0d2/ε2. The corresponding dimensionless variables become
= x/L, ĝ = (1 −
)−1, and ĵ =
, with the characteristic frequency given by ω0 = δ/(ε0WLR). This similarity framework underscores a fundamental distinction between TENGs and electromagnetic generators. While EMGs with analogous structures exhibit unconditional similarity due to the absence of governing dimensionless groups in their dynamic equations (as detailed in the SI), the existence of the Sm criterion in TENGs fundamentally explains their rich and complex dynamic behaviours. The presence of this scaling parameter enables both precise performance prediction and systematic device optimization across diverse configurations and operating conditions.
(
) possesses a unit period. The general solution for the dimensionless transferred charge can be expressed as:
=
− N ∈ [0,1) and N is a non-negative integer. This representation demonstrates that the charge dynamics over the semi-infinite time domain
∈ [0,+∞) are completely determined by their behavior within a single, finite period
∈ [0,1].
As time evolves (
→ +∞), the ratio of successive periods satisfies
(
+ 1)/
(
) → 1, confirming the asymptotic convergence to a steady-state periodic solution
S(
). This periodic limit cycle can be compactly defined for
∈ [0,1) as:
This establishes that despite the complex, state-dependent capacitance, TENGs under periodic excitation inevitably converge to a unique periodic operating regime.
, indicating that the dimensionless transferred charge approaches a constant value:Physically, this limit represents insufficient time for charge redistribution during rapid oscillations, causing the transferred charge to saturate. Consequently, the current and voltage become:
Remarkably, in this limit, the output current becomes independent of both the excitation frequency T−1 and the device's primary dimensions L and W, imposing a fundamental upper bound on power generation. This saturation behavior, absent in EMGs, reveals a unique scaling law for TENGs. Furthermore, the output voltage becomes load-independent, characterizing the TENG as an ideal voltage source in this operational regime.
This indicates complete charge transfer tracking the short-circuit value QSC. The resulting current and voltage scale as:
In this regime, both current and voltage exhibit conventional 1/T scaling, similar to EMGs. However, a fundamental distinction remains: the TENG operates as a current source, whereas the EMG behaves as a voltage source, reflecting their different energy conversion mechanisms.
The contrasting limiting behaviours produce a distinctive scaling law for the average output power
S. For EMGs, power universally scales as
∝ T−2. For TENGs, the scaling is regime-dependent:
As T → 0 (Sm → 0):
The power saturates to a finite constant, independent of both frequency and device dimensions.
As T → ∞ (Sm → ∞):
The power recovers the conventional T−2 scaling.
At finite periods, the average power
smoothly transitions between these asymptotic limits, creating a universal scaling curve that encapsulates the fundamental performance limits of TENGs across all operating conditions.
(
). This finding prompts a fundamental question: what form of
(
)—hereafter termed the shape function—maximizes the average output power (
ST→0) under this limiting condition?
The theoretical optimum is a top-hat function, given by:
m). This waveform yields the maximum achievable power,The optimal shape function for T is a top-hat form, fully determined by the structural characteristic ĝ(
) and the displacement limit
m. The resulting current waveform is synchronous with the motion input, suggesting the potential for a mechanical power management strategy.
We further employed numerical simulations to evaluate the performance of
opt at finite periods, comparing it against harmonic (
Har), triangular (
Tri), and two composite waveforms (
Haropt,
Triopt; see the SI). The corresponding
S–T curves are shown in Fig. 2. Evidently, the power output for
opt markedly surpasses that of other inputs. Furthermore, the composite waveforms
Haropt and
Triopt outperform their pure counterparts, indicating that intermittent and abrupt motion profiles enhance the average power output under identical motion periods.
It should be noted that the optimal function
opt is derived under the condition where period T approaches zero, and its advantages are maintained under relatively high frequency conditions.
Haropt and
Triopt over pure harmonic and triangular waves at high frequencies requires careful analysis. Since these composite waveforms feature steeper transitions, we develop the concept of motion input decomposition to compare waveforms with identical maximum velocities.
As shown in Fig. 3, we analyze a sequence of motions where motion 1 (a single continuous cycle) is progressively decomposed into multiple intermittent sub-cycles: motion 2 (two sub-cycles), motion 3 (three sub-cycles), motion 4 (four sub-cycles), motion 5 (six sub-cycles), and motion 6 (twelve sub-cycles). The mathematical form of these decomposed motions is given by:
) = Ni
for 0 ≤
≤ Ni−1, and Πi(
) represents a rectangular window function that is unity within each active sub-interval and zero elsewhere. The decomposition levels correspond to Ni = {1, 2, 3, 4, 6, 12} for i = 1–6, respectively.
To ensure kinematic consistency, we define the effective period Teff = T·Ni−1, representing the duration of each active motion segment. This normalization guarantees identical maximum velocities when comparing different waveforms at the same Teff. The power output analysis reveals a distinct optimum in the decomposition strategy. Motion 2 demonstrates enhanced power generation at shorter effective periods compared to the continuous motion 1. However, this trend reverses with excessive decomposition – motions 3 through 6 exhibit progressively reduced performance across all frequencies. This establishes motion 2 as the optimal decomposition level, whose waveform characteristics most closely resemble the beneficial properties of
Haropt, confirming that moderate intermittency rather than maximal decomposition maximizes power output.
For a specified motion profile, theoretical analysis confirms the existence of a unique dimensionless parameter
that maximizes power extraction. The corresponding optimal load resistance follows the scaling law:
Experimental validation across four distinct operational regimes—varying in device dimensions and excitation period—confirms this relationship. The power-resistance characteristics consistently peak at
despite parametric variations, demonstrating the universal criterion for impedance matching in TENG systems.
For harmonic excitation profiles, the optimal similarity parameter admits a closed-form expression:
The characteristically large magnitude of Ropt emerges naturally from this scaling relationship, accounting for the fundamental constant and typical displacement ranges in practical systems.
Further analysis reveals a profound physical interpretation: maximum power transfer occurs when the mechanical excitation frequency synchronizes with the system's electrical response. This resonant condition satisfies:
The general expression for optimal load resistance under resonant operation becomes:
. Experimental validation using published data confirms the predictive capability of this formulation, with the model achieving agreement within one order of magnitude across diverse experimental configurations. This resonant impedance matching framework provides a physical basis for TENG optimization, advancing beyond empirical approaches toward principled design methodologies.
The operation of EMGs is unconditionally similar, whereas the similarity of TENGs depends exclusively on the similarity number Sm = ω0T, where ω0 represents the characteristic frequency of the TENG structure and T is the motion input period. The operating states of TENG systems can be categorized based on the value of Sm into high-frequency (Sm ≪ 2π), medium-frequency (Sm ∼ 2π), and low-frequency (Sm ≫ 2π) states, each exhibiting distinct dynamic characteristics.
The output of an EMG depends solely on the instantaneous motion input, while the output of a TENG exhibits history-dependent behaviour. Under periodic motion excitation, the EMG output is inherently periodic, whereas the TENG output asymptotically approaches a periodic solution. In terms of power source characteristics, the EMG behaves as a voltage source, while the TENG exhibits voltage-source characteristics in the high-frequency state and transitions to current-source behaviour in the low-frequency state.
The energy output scaling reveals another fundamental difference: the average output power of an EMG follows an inverse-square relationship with the input period (
∝ T−2), while the TENG maintains this relationship only in the low-frequency state but saturates to a frequency-independent constant in the high-frequency state. Furthermore, in the high-frequency state, the TENG exhibits an optimal motion input in the form of a top-hat function, whose specific form is determined solely by displacement constraints and structural characteristics. The power output can be enhanced through appropriate motion decomposition and intermittent operation strategies, and a resonance state exists where the matching resistance corresponds to the condition where the time-averaged oscillation frequency equals the input motion frequency.
The distinctive dynamics of TENGs suggests that their design and application demand innovative approaches fundamentally different from those used for EMGs. The significantly greater complexity of TENG dynamics, particularly the strong coupling among motion input parameters, structural parameters, and load circuit parameters, indicates that advancing TENG technology requires collaborative efforts and deep integration across multiple disciplines, including materials science, structural design, power management, and dynamic analysis. These findings not only provide a theoretical foundation for TENG optimization but also open new research directions for overcoming fundamental limitations in mechanical energy harvesting systems.
The charge state evolves according to:
Initialization is defined by:
[0] = 0,
[0] = 1.
For determining steady-state periodic solutions, we derived a direct computational approach:
Additional data related to numerical codes and model validation can be provided by the corresponding author upon reasonable request.
| This journal is © The Royal Society of Chemistry 2026 |