Results and Discussion
Figure 1 presents the core design of a self-powered platform for metabolite biomarker analysis utilizing an LS-TENG sensing array assisting TES. This platform exploits the triboelectric signals produced by different analyte droplets sliding across the FEP film to construct distinctive triboelectric spectra, thereby enabling label-free identification of solutions with different chemical compositions. As depicted in
Figure 1a, the operational workflow begins when droplets containing specific chemical components—in this study, solutions of glucose, lactic acid, urea, uric acid, sodium chloride, potassium chloride, calcium chloride, and ammonium chloride solutions—traverse an inclined FEP surface. The LS-TENG sensing array can readily capture the interfical triboelectrification properties induced by electron transfer and ion adsoprtion and represent this process with specific temporal electrical signals, which are closely related to the analyte species. A high-precision multichannel electrometer (NI PXIe-8880) is utilized to continuously record the transferred electrostaic charge
Q from each sensing electrode as the target analyte droplet sequentially traverses the 16 electrically isolated copper electrodes. Due to inherent differences in interfacial triboelectrification, solutions of varying chemical composition generate distinct and characteristic spatial-temporal distribution curves of
Q, i.e., the representative TES for different analytes. These spectral profiles form the foundational dataset for subsequent chemical/biomolecular identification and quantification. The structural schematic of the LS-TENG sensing array for TES is presented in
Figure 1b. The device employs a multilayer architecture: a PMMA plate serves as the supporting substrate; a patterned 16-copper-electrode array fabricated on this substrate serves to electrostatically sense the charge variations induced by the passing droplets; and a top layer of FEP film, which constitutes the critical liquid–solid interface for the triboelectric charge generation. The interelectrode spacing is optimized to 0.5 cm to effectively minimize electrical crosstalk between adjacent electrodes (Figure S1). By processing the acquired triboelectric spectral data with a convolutional neural network (CNN), as outlined in
Figure 1c, the system achieves accurate classification of the eight distinct biofluid analytes. This integrated approach not only facilitates single-droplet diagnostics but also exhibits significant potential for extension to complex physiological biofluid analysis and for the development of scalable, self-powered biofluid sensing networks.
The biofluid analysis platform employs a FEP film as the triboelectric dielectric overlayer. This layer electrically isolates the droplet from the metal electrodes, so that the readout is governed primarily by electrostatic induction through the dielectric. When an analyte droplet slides on the FEP surface, charge separation at the liquid/FEP interface generates interfacial charges and a localized electric field. As the droplet advances along the array, the field distribution evolves with its position as it approaches, covers, and then departs from successive electrode regions. The resulting time-varying field couples capacitively to the underlying copper strip electrodes and drives transient redistribution of induced charge on each channel, producing electrode-resolved electrical signals. By recording these induced charge variations with the multichannel acquisition system and indexing them by electrode order, a spatiotemporally resolved response map can be constructed, defined here as the triboelectric spectrum (TES) of a single droplet transit.
Figure 2a shows a magnified view of a deionized (DI) water droplet sliding across the FEP-covered electrode array, highlighting the dielectric isolation and the liquid–solid interface that underpin the induction-based signal generation. The upper panel of
Figure 2a further illustrates the temporal evolution at defined electrode coordinates: as the DI water droplet sequentially traverses the copper strips, distinct pulse waveforms arise on individual channels, as exemplified by the signals from the first three electrodes. In this manner, continuous droplet motion is converted into a sequence of 16 spatially indexed pulses, establishing a comprehensive triboelectric profile for the fluid sample. The selection of FEP film as the triboelectric layer is primarily attributed to its superior surface properties. Its strong hydrophobicity facilitates rapid and stable droplet motion, while its pronounced electron affinity enhances interfacial charge transfer between the droplet and the surface, thereby amplifying the strength of the resulting triboelectric signal (Figure S3).
Figure 2b shows representative high-speed images of droplets sliding on the FEP film. The corresponding video (Video S1, Supporting Information) is provided to visually document the continuous sliding process along the device. The static contact angle on FEP is also characterized to be 96.75˚, serving as further confirmation of its excellent hydrophobicity and the controllability of droplet motion.
During the experiments, the analyte droplets (≈30 μL per drop) were dispensed at a constant flow rate from a grounded stainless-steel needle using a syringe pump. The needle was positioned approximately 1 cm above the FEP surface. As these droplets slid along the FEP film inclined at 50°, the transferred electrostatic charges
Q at different time points synchronously lead to the inducted charges on the underlying copper-electrode array and recorded in real time with the multichanel electrometer. Using DI water as a representative case, the entire detection time—from droplet release to complete traversal of the FEP surface—was completed in less than 1 second. The transferred charge at each electrode position was calculated from the acquired current signals. Based on 20 independent experimental trials of the liquid-solid sliding process (Figure S4), a representative TES was established as shown in
Figure 2c. The TES of DI water is dominated by three prominent characteristic peaks, located at approximately 6 cm, 11 cm, and 14 cm along the sliding path. These local maxima are not random artifacts but correspond to specific hydrodynamic stages of the droplet’s motion. The first peak (≈ 6 cm) arises from the initial capillary-inertial oscillation, where the droplet undergoes its first major deformation cycle upon acceleration, temporarily maximizing the contact area. Subsequently, the signal evolves into a broader, double-humped structure comprising the second (≈ 11 cm) and third (≈ 14 cm) peaks. This region represents the stable sliding phase, where the droplet achieves kinetic equilibrium. The global maximum at 14 cm specifically signifies the point of charge saturation, attributed to the optimal spreading of the droplet and the fully established electric double layer at the interface. The discernible dip between 11 cm and 14 cm likely reflects a periodic contraction of the droplet shape, further confirming the high sensitivity of the TES to dynamic fluid behaviors.
To elucidate the sensing mechanism underlying the electrode-resolved TES,
Figure 2d (I–IV) illustrates a representative process as the droplet (DI water) passes over the triboelectric sensing electrodes (taking the first four electrodes as the examples). Driven by gravity, the droplet slides along the FEP surface, successively transvers each electrode. Owing to the strong electron affinity of FEP, contact electrification occurs immediately upon intial contact, resulting in electron transfer from water molecules to the FEP surface. This processs renders the FEP negatively charged and the droplet positively charged. During this interaction, a fraction of water molecules ionize to form H
2O
+, which rapidly converts to H
3O
+. These cations migrate within the droplet and actively participate in the interfacial charge dynamics. When the droplet reaches the first copper electrode, the electric field established by the positively charged droplet and negatively charged FEP drives H
3O
+ migration toward the electrode surface. This ion migration, coupled with charge transfer through the electrode and external circuit, establishes a potential difference and results in a measurable transferred-charge signal. As the droplet advances to the second electrode, repeated contact-electrification effect and contact-separation process promote further charge partitioning, accumulation, and redistribution between adjacent electrodes. Each discrete contact–separation cycle between the droplet and one electrode modifies the overall charge distribution, yielding a corresponding transferred-charge signal on each successive electrode. To quantitatively characterize the charge-transfer dynamics during droplet motion, an electrostatic potential model was implemented using the COMSOL Multiphysics finite-element platform.
Figure 2e displays the simulated potential distributions, whch correspond to the experimental stages depicted in
Figure 2d I–IV. The close agreement between simulation and experiment is evident, with cross-sectional potential maps clearly delineating the dynamic evolution of potential differences across the electrode array.
Based on the foregoing discussion of the system design and signal acquisition, this section systematically elucidates the triboelectric response behaviors and governing mechanisms of representative cationic and non-ionic solutions on the LS-TENG based TES sensing platform.
Figure 3a–d dispaly the transferred-charge distributions—that is, the triboelectric signals induced at different positions as a droplet slides along the electrode array—for K⁺, NH₄⁺, Na⁺, and Ca²⁺ solutions at an identical concentration (5 mM). Although all four solutions share the same anion (Cl⁻), their triboelectric spectra exhibit pronounced differences, indicating that the physicochemical properties of the cation predominantly dictate charge-transfer behavior at the liquid–solid interface. During the triboelectrification process, immediate electron transfer occurs from water molecules to FEP surface upon initial droplet contact, rendering the FEP negatively charged and the droplet positively charged [
40,
41]. Subsequently, cations within the droplet are electrostatically attracted to and adsorb on the negatively charged FEP, significantly modulating charge induction [
38,
42]. The overall interfacial charge-transfer dynamics can be approximated by:
where
denotes the contribution from electron transfer, and
represents that contribution from ion adsorption and desorption [
43]. Consequently, differences in interfacial mobility and adsorption kinetics among various ions directly shape the resulatnt TES profiles.
Figure 3e presents a schematic illustration of this ion adsorption–desorption mechanism. As the droplet traverses the FEP surface, a portion of cations is captured by the FEP, while another fraction is carried away by the moving droplet, thereby generating localized charge variations. The magnitude and temporal characteristics of the charge transfer depend not only on cation mobility but also on the stability of the “ion atmosphere” surrounding each ion in proximity to the interface [
44]. In bulk solution, a cation is typically enveloped by a spherically symmetric, negatively charged cloud to maintain electroneutrality. However, near the negatively charged FEP surface, the interfacial electric field perturbs this arrangement, partially separating anions from cations and altering the migration pathway and velocity of the central ion. The continuous rupture and reformation of this ion atmosphere dynamically regulate the triboelectric signal output.
Figure 3f–g systematically analyze the influence of ionic structural parameters on the triboelectric response. The hydrodynamic radius, atomic mass, and valence of each ion collectively determine its interfacial migration efficiency, which in turn governs the onset time and amplitude of characteristic features in the TES. For instance, K
+ and NH
4+, which possess relatively small hydrodynamic radii (~1.296 Å and 1.25 Å, respectively), display higher interfacial mobility, and consequently exhibit earlier signal onset, with their characteristic peaks appearing closer to the initial electrodes. Although Ca²⁺ is divalent and experiences stronger electrostatic attraction toward the negatively charged surface, its more complex hydration structure substantially limits its migration speed, resulting in response signals that slightly lag behind those of K
+ and NH
4+. In contrast, Na⁺, with the largest hydrodynamic radius among the four ions (≈1.84 Å), encounters the greatest migrational resistance and thus requires more time to reach and interact effectively with the FEP interface; its corresponding charge peaks therefore emerge markedly later, reflecting the slowest overall response. These experimental observations indicate that triboelectric response is not solely determined by ionic valence states but is strongly governed by interfacial migration kinetics, particularly the decisive effect of hydrodynamic size on adsorption rates. Accordingly, the cation’s specific size, mass, and associated hydration properties cooperatively regulate its migration and adsorption dynamics at the FEP interface, ultimately dictating the peak positions, response speeds, and magnitudes of transferred charge observed in the triboelectric spectra across different ionic species.
For polar (non-ionic) solutions,
Figure 3h–k compares the triboelectric responses of glucose, urea, lactic acid, and uric acid under identical conditions (5 mM concentration with matched sliding profiles). All four analytes exhibit a characteristic signal profile comprising an initial rapid rise, a brief or extended plateau, and a subsequent slow decay. Crucially, the timing of the first detectable peak follows a distinct and reproducible sequence: glucose → urea → lactic acid → uric acid. This ordering is determined not by variations in signal amplitude but prmarily by the interplay between the adsorption kinetics of polar functional groups and the formation rate of electrical double-layer (EDL) at the interface. Upon contact and separation from the FEP surface, charge separation occurs at the liquid–solid interface. This process, involving selective adsorption of ions or molecular dipoles, imparts a negative charge to the FEP and a corresponding positive charge to the droplet. The result is a train of transferred-charge pulses that propagates sequentially along the electrode array [
45]. For glucose solutions, the abundance of hydroxyl (–OH) groups facilitates the formation of dense hydrogen-bond networks with interfacial water molecules, promoting rapid molecular orientation, This establishs an efficient, high-throughput interfacial charge-transfer pathway during the initial contact and early sliding phases, leading to the early emergence of the first peak. As adsorption progresses and the nascent EDL develops, charge injection is gradually attenuated, giving rise to the observed plateau and late-stage decay [
38,
46]. Urea, though highly polar due to its carbonyl (C=O) and amine (–NH₂) groups, must undergo a reorientation from its bulk-favored configuration to an interface-favored one to establish a stable hydrogen-bond network and achieve optimal dipole alignment. This requisite reorientation delays its initial response, shifting the first peak to a later position than that of glucose; the plateau region thereafter reflects a dynamic equilibrium between ongoing charge injection and increasing electrostatic screening from the EDL. Lactic acid partially dissociates in solution, thereby increasing the local ionic strength at the interface and accelerating EDL formation. This early development of the EDL, coupled with competitive adsorption at the surface, suppresses the initial electron/ion transfer kinetics, thereby shifting the first peak further downstream and limiting its overall growth. Uric acid, featuring multiple hydrogen-bond donor/acceptor sites and a tendency for intermolecular association, forms a particular dense interfacial adsorbate layer. The slower dynamics of interfacial reorganization and molecular desolvation, combined with stronger screening effects, result in the least efficient early-stage charge injection. Consequently, its first peak emerges at the latest position, accompanied by a more gradual signal transition (Figure S5). To capitalize on these distinct kinetic signatures in the subsequent machine-learning analysis, we standardized all sample concentrations to 5 mM. This normalization decouples the signal from amplitude–concentration linearity, ensuring that variations in ‘first-peak timing’ serve as robust indicators of intrinsic molecular composition and interfacial kinetics, rather than confounding concentration effects.
To quantitatively analyze concentration gradients from the characterized triboelectric signals, we employ three-dimensional (3D) waterfall plots to visualize the evolution of the droplet-generated response as a function of solute concentration during sliding. As shown in
Figure 4a, a tri-axial data framework is constructed, plotting the transferred charge
Q along the x-axis, the electrode position within the array along the y-axis, and the solution concentration along the z-axis. This approach generates a volumetric triboelectric spectral map that facilitates metabolite identification with concentration-resolved resolution. As the droplet traverses and sequentially activates the electrode array, the generated triboelectric response trajectory simultaneously captures information on both the interfacial charge transfer and the concentration-dependent interfacial chemistry. This dual encoding produces a distinctive “concentration-sensitive” triboelectric signature for analytical discrimination.
Figure 4b–e present 3D waterfall plots for four ionic solutions (NaCl, KCl, CaCl₂, NH₄Cl) across a range of concentrations, revealing a systematic, monotonic decrease in transferred charge
Q with increasing solute concentration. This trend indicates that solute concentration directly modulates the strength of interfacial triboelectrification, highlighting the close coupling between solute loading and electrical output. At low concentrations, cations (Na⁺, K⁺, Ca²⁺, NH₄⁺) can adsorb relatively uniformly onto the negatively charged FEP surface, establishing a robust interfacial electrical field that results in a larger measured
Q. However, as the concentration increases, surface adsorption gradually approaches saturation; the rate of increase in
Q progressively slows and eventually plateaus. This is because the available adsorption sites become largely occupied and charge-transfer efficiency reach its limit. With further elevation in concentration, a reduction in
Q is observed, which is consistent with the onset of electrostatic screening effects: a denser interfacial adsorption layer weakens the potential difference between the FEP and the droplet, thereby suppressing electron transfer and reducing the detectable charge signal.
In contrast to the consistent trend observed with ionic solutions, the non-ionic (polar) solutions—namely glucose, urea, lactic acid, and uric acid—exhibit distinct concentration-dependent behaviors (
Figure 4f–i). For glucose and urea, the transferred charge
Q increases with rising concentration. Although these molecules carry no net charge, their polar functional groups engage in strong interactions with the FEP interface, thereby facilitating electron transfer. At elevated concentrations, a greater number of molecules are recruited to participate in these interfacial processes, which enhances the overall triboelectric signal output. In comparison, lactic acid and uric acid demonstrate the inverse dependence:
Q decreases as concentration increases. This contrasting behavior can be attributed to acidity-induced variations in hydronium (H₃O⁺) concentration and the consequent redistribution of interfacial charge. Under conditions of increased concentration, the heightened presence of H₃O⁺ intensifies electrostatic screening and introduces competitive adsorption at the interface. These factors collectively reduce the efficiency of electron transfer, leading to the observed decline in measurable charge. Collectively, these results establish that the LS-TENG output is determined not simply by analyte concentration, but by a complex interplay of chemical composition, the kinetics of ion adsorption and desorption, and local gradients in H₃O⁺ concentration. These parameters jointly determine the amplitude and spatiotemporal evolution of the resultant triboelectric spectra.
In scenarios requiring home- or field-deployable rapid testing, the transferred-charge sequences generated by different metabolite-related solutions under identical conditions are often visually indistinguishable Subtle discriminative features are further obscured by non-stationary fluctuations stemming from droplet velocity perturbations, wetting hysteresis, contact-line pinning, and residual static charge. To circumvent the confounding effect of amplitude–concentration collinearity on classification, the concentration of all samples used for machine learning was standardized at 5 mM. This deliberate design choice focuses the model exclusively on fine-grained temporal structure, which is intrinsically governed by the chemical composition of each analyte.
Figure 5 systematically presents the complete data processing pipeline—from the raw LS-TENG time series to multi-class discrimination—under uniform concentration (5 mM) and identical operating conditions, along with the corresponding quantitative evaluation results. To avoid the confounding effect of amplitude–concentration collinearity on classification performance, all samples used for modeling and evaluation were standardized at this concentration. The overall workflow, depicted in
Figure 5a, initiates as a droplet sweeps sequentially across the over 16 equally sized copper electrodes beneath the FEP overlayer. Each contact-separation event induces a measurable transferred charge at the corresponding electrode location. These electrical signals are acquired in real time by an NI PXIe-8880 data acquisition system before being passed to subsequent preprocessing and event extraction modules. Preprocessing consists of baseline correction (to remove slow drift and DC offset), band-pass denoising and smoothing (to suppress both low-frequency undulations and high-frequency noise spikes), followed by precise peak detection and pulse segmentation using a threshold–slope joint criterion, with temporal alignment enforced by the known electrode sequence.
Figure 5b displays representative time-domain signal overlays for eight investigated solutions (NaCl, KCl, CaCl₂, NH₄Cl, glucose, urea, uric acid, lactic acid) under identical sliding dynamics. While signals within the same class demonstrate stable consistency in the timing of the first-peak, pulse duration, and overall envelope morphology, the most discriminative inter-class differences reside precisely in these subtle spatiotemporal patterns. This characteristic makes visual differentiation highly unreliable, thereby necessitating the adoption of a statistical learning methodology. To explicitly retain the spatial ordering information across electrodes, the segmented and temporally aligned signal trajectories are transformed into a 2D input sensor with dimensions “time × electrode-segment” (size 639×15×1). Feature learning is performed by a three-layer convolutional neural network, which applies convolutional filters and hierarchical pooling operations to this structured input. The detailed network architecture and the corresponding tensor size evolution are illustrated in
Figure 5c: Conv2D (32 filters, 3×3 kernel, same padding, ReLU) → MaxPool (2×2 window) yields 320×8×32; Conv2D (64 filters, 3×3 kernel, same padding, ReLU) → MaxPool (2×1 window) yields 160×8×64; Conv2D (128 filters, 3×3 kernel, same padding, ReLU) → MaxPool (2×1 window) yields 80×8×128; then Flatten (81920 units) → Dense (256 units, ReLU) → Dense (8 units, softmax) produces the final eight class labels. It should be noted that the “8” along the electrode axis in the diagram arises from the first pooling operation (15→8) and is not related to the number of input channels. Key hyperparameters—including kernel size, pooling window dimensions, stride values, and the width of the fully connected layer—were jointly optimized via Bayesian optimization (300 iterations) over a predefined search space. This approach was designed to balance the effective receptive field, total parameter count, and model generalization capability. The architecture of the final layer and the detailed parameter breakdown are provided in Figure S7, with the network containing a total of 21 066 504 trainable parameters. Using this optimized configuration, five-fold cross-validation was performed. The resulting confusion matrix (
Figure 5d) is predominantly diagonal, with misclassifications concentrated among the most spectrally similar class pairs. The training/validation accuracy curves rise synchronously and begin to plateau after approximately 30–40 epochs (
Figure 5e), demonstrating stable convergence without significant signs of overfitting or high variance. The overall mean classification accuracy reaches 95.3% (Figure S6). Furthermore, the receiver operating characteristic (ROC) curves (
Figure 5f), characterized by false-positive rate (FPR) versus true-positive rate (TPR), cluster tightly in the upper-left quadrant. This pattern signifies that the model achieves a high recall rate at a minimal false-positive cost, thereby facilitating application-specific threshold tuning for practical deployment. In sum, the integrated analytical framework—encompassing “uniform-concentration control → standardized preprocessing and pulse segmentation → a three-layer CNN with controlled pooling → Bayesian hyperparameter optimization → cross-validation-based performance assessment—can successfully extract generalizable compositional fingerprints from the TES-assisted LS-TENG transferred-charge time series, even in the presence of realistic perturbations and highly similar waveforms, achieving robust discrimination among the eight different solution classes.