TMPA-VOH was synthesized via a hydrothermal method. 2 mmol of V₂O₅ and 2 mL of H₂O₂ (30 wt%) were added to 50 mL of deionized water. In a separate beaker, 1 mmol of trimethylphenylammonium chloride was dissolved in 30 mL of deionized water. The two solutions were then mixed and stirred for 30 min. The resulting solution was transferred to a 100 mL Teflon-lined autoclave and maintained at 120 °C for 6 h. Then the precipitates were collected by centrifugation and washed with deionized water and ethanol. The precipitates were first air-dried at 70 °C overnight and then vacuum dried at 120 °C for 2 h. For comparison, VOH was synthesized using the same procedure, excluding the addition of the organic salt and was dried using a freeze dryer.

Phase identification was carried out using a powder X-ray diffractometer (XRD, Bruker D8 Discover). Microstructural, morphological, and elemental analyses were conducted using FEI Sirion XL30 scanning electron microscope (SEM) equipped with the Oxford Instruments Energy Dispersive X-ray Spectrometer (EDS) system, and an FEI Tecnai G² F20 SuperTwin 200 keV Transmission Electron Microscope (TEM) equipped with an EDAX Elite-T EDS. The surface chemical composition was quantified using a Kratos Axis Ultra DLD X-ray photoelectron spectroscopy (XPS) with a monochromatic Al Kα x-ray source. Thermal stability and water content were measured using a dual thermogravimetric analysis/differential scanning calorimetry instrument (Mettler Toledo TGA/DSC 3 +). Functional groups and chemical bonds were identified by a Fourier transform infrared (FTIR) spectroscopy (Thermo Scientific Nicolet iS10) with laser excitation at 514 nm. Electron paramagnetic resonance (EPR) tests were conducted using a Bruker A300 spectroscopy.

The TMPA-VOH and VOH cathodes were prepared as follows: first, active materials, carbon black, and polyvinylidene fluoride were blended in a weight ratio of 7:2:1 and dispersed in N-Methyl-2-Pyrrolidone solvent to form a uniform slurry. Then, the slurry was coated onto a titanium foil and dried at 60 °C, followed by additional drying in a vacuum oven at 120 °C for 12 h. Electrochemical tests were conducted using CR2032 coin cells, with metallic zinc serving as the counter electrode, glass fiber as the separator, and a 3 M zinc trifluoromethanesulfonate (Zn(CF₃SO₃)₂) aqueous solution as the electrolyte. All cells were assembled in air and tested within a voltage range of 0.2-1.6 V (vs. Zn²⁺/Zn). Galvanostatic charge and discharge (GCD) tests were performed using the Neware Battery Testing Systems (CT-4008). Cyclic voltammetry (CV) and electrochemical impedance spectroscopy (EIS) data were collected on the Solartron SI 1287 Potentiostat and Galvanostat equipped with a 1260 Solartron Impedance Analyzer.

TMPA-VOH is synthesized using the same hydrothermal conditions as VOH, with the addition of trimethylphenylammonium chloride (TMPACl) salts at a V₂O₅/TMPACl ratio of 2. The formation of dark green precipitates implies a significant reduction of V⁵⁺ to V⁴⁺ [28]. Figure 1a illustrates the X-ray diffraction (XRD) patterns of TMPA-VOH and VOH. VOH exhibits a typical layered structure with characteristic peaks (001), (003), (004), and (005), well-matched to the V₂O₅·1.6H₂O pattern (JCPDS No. 40-1296) [14]. The main (001) peak at 2θ of 7.54° corresponds to an interlayer distance of approximately 11.7 Å in VOH. However, upon introduction of TMPA⁺ cations, the diffraction peaks of TMPA-VOH align with those of the monoclinic (Na, Ca)(V, Fe)₈O₂₀·nH₂O phase (Space group: C2/m, JCPDS No. 45-1363) [29]. As clearly depicted in Fig. 1a, all the (00l) peaks of TMPA-VOH shift to a lower degree compared to the standard (Na, Ca)(V, Fe)₈O₂₀·nH₂O phase. Meanwhile, all the other peaks such as (201), (110), (40 $\overline{1 }$), (11 $\overline{3 }$), (60 $\overline{1 }$), (020), (021), and (71 $\overline{1 }$) remain at the same positions, strongly indicating that TMPA⁺ has entered the interplanar space of the V-O structure. Compared to the standard (Na, Ca)(V, Fe)₈O₂₀·nH₂O pattern, the (00l) peak of TMPA-VOH shifts from 8.18° to 6.36°, corresponding to an expansion of lattice spacing from 10.8 to 13.9 Å. A schematic illustration of the TMPA-VOH structure is provided in Fig. S1. The [V₈O₂₀] layer is composed of four-fold [V₄O₁₀] chains, formed by (opposite) corner sharing [VO₆] octahedra (parallel to b). Those chains densely packed by sharing edges to form layers, with sufficient interlayer sites partially occupied by cations or water molecules [29].

Scanning electron microscopy (SEM) images (Fig. 1c) and transmission electron microscopy (TEM) images (Fig. S2a) reveal that TMPA-VOH exhibits a three-dimensional flower-like/urchin-like morphology, comprising numerous long and ultrathin VO_x nanobelts. These nanobelts nucleate from the core and grow radially outward in all directions, typically reaching several micrometers in length and less than 100 nm in width (Figs. 1d and S2b, c). This suggests a rapid diffusion-controlled crystal growth without secondary nucleation. The high-resolution transmission electron microscopy (HRTEM) image of TMPA-VOH (Figs. 1e and S2d) displays lattice fringes with an interlayer spacing of 1.42 nm, corresponding to the d-spacing of (001) planes from the XRD pattern (Fig. 1a), implying a preferential growth direction for the nanobelts. In contrast, VOH exhibits a more agglomerated and irregular morphology (Fig. S3), consisting of randomly stacked faceted flakes. This difference implies that the incorporation of TMPA⁺ cation or the presence of TMPACl induces a significant morphology change.

Possible explanations for this morphological change might be: the uniformly distributed TMPA⁺ cations in the V₂O₅ sol could (1) provide more initial nucleation sites, (2) alter the surface energy of nuclei and crystals by adsorbing TMPA⁺ cations onto the surface, (3) lower the critical energy barrier of initial nucleation and/or subsequent crystal growth, and (4) slow the supply or diffusion of vanadium oxide growing species to the growth front. These factors lead to the formation of a morphology far from the thermodynamic equilibrium shape. The 3D micro flower-like morphology with ultrathin nanobelts of TMPA-VOH provides a larger accessible surface area for ion diffusion and more open space for electrolyte penetration, reducing the diffusion pathways for guest ions. TEM energy dispersive spectrum (TEM-EDS) element mapping (Fig. 1f) and SEM-EDS (Fig. S4) confirm the uniform distribution of C, N, V, and O throughout the entire TMPA-VOH architecture, validating the successful introduction of TMPA⁺ into the V-O lattice. The molar ratio of N:V, calculated to be ~ 13.6% based on the SEM-EDS and TEM-EDS, gives an average stoichiometric formula of (TMPA)_1.08V₈O₂₀·nH₂O for TMPA-VOH.

The water content in TMPA-VOH and VOH is determined through thermogravimetric (TG) analysis. In Fig. 1b, the water loss in VOH is divided into: (1) a first stage (< 150 °C) corresponding to the loss of weakly bonded water molecules (reversibly adsorbed between the V-O layers), and (2) subsequent stages (150-450 °C) corresponding to the removal of more tightly bonded water within the layers [30]. The stoichiometric formula of the prepared VOH was determined to be V₂O₅·1.21H₂O, with a total mass loss of 10.4% between 100 to 450 °C. Notably, in TMPA-VOH, there is a negligible mass loss in the first stage (below 150 °C), indicating the absence of reversibly adsorbed water molecules between layers. It can be inferred that intercalating TMPA⁺ cations into the V-O ribbons is accompanied by the expulsion of interlayer water molecules. The primary weight loss of TMPA-VOH (~ 17%) occurs between 150 and 450 °C, attributed to the removal of strongly bonded crystal water and the combustion of organic components in air. Based on the N content from EDS data (N: V ~ 13.6%), the amount of TMPA⁺ cation and crystal water can be estimated, determining the formula of TMPA-VOH roughly to be (TMPA)_1.08V₈O₂₀·0.06H₂O.

The combined results of XRD, SEM/TEM and TG affirm that the intercalation of large, weakly polarized organic TMPA⁺ cations induces a phase change from V₂O₅·nH₂O to (TMPA)_1.08V₈O₂₀·nH₂O. This transformation is characterized by an interlayer spacing expansion from 11.7 to 13.9 Å, a morphology change from randomly stacked faceted flakes to 3D spherical micro-urchins, and the extrusion of weakly bonded water molecules between interlayers.

The chemical composition and valance state of TMPA-VOH are elucidated through XPS. The XPS survey spectra of TMPA-VOH (Fig. 2a) reveal signals from V, O, N, and C, confirming the presence of TMPA⁺ cations in the material. Notably, no trace of Cl signal is detected in the survey spectra (Fig. 2a) or the high-resolution Cl 1s spectra (Fig. S5), indicting the exclusive insertion of organic TMPA⁺ cations without Cl⁻ anions. In the high-resolution C 1s spectra of VOH (Fig. 2b), in addition to the main line at the binding energy of 285.0 eV (C-C/C-H bonding), two other weak but recognized peaks at 286.7 and 288.5 eV are observed, assigned to the C-O, and O-C = O components [31]. These detectable carbon components in VOH originate from surface contamination, commonly observed in samples exposed to the atmosphere [32,33]. In comparison to the C 1s of carbon contamination in VOH, the C 1s spectra of TMPA-VOH (Fig. 2b) reveal a strong resolved peak at 286.5 eV, corresponding to the C-N bonds from the amino group in TMPA⁺. The presence of C--N bonds or the amino group is also evident at 400.6 eV in the N 1s spectra (Fig. 2c), where the other peak at 402.9 eV could be attributed to N-O bonds formed between TMPA⁺ cations and oxygen in the vanadium oxide layer [34,35]. This further verifies the successful pre-intercalation of TMPA⁺ cations into the VOH lattice and the formation of chemical bonds between nitrogen atoms (from TMPA⁺) and oxygen atoms (on the interlayer surface of vanadium oxide).

The high-resolution V 2p_3/2 spectrum of TMPA-VOH (Fig. 2d) reveals a prominent peak at 517.9 eV and a shoulder at 516.2 eV, corresponding to V⁵⁺ and V⁴⁺ ions, respectively [31,36,37]. From the fitted peak area, the ratio of V⁴⁺/V⁵⁺ in TMPA-VOH is calculated to be 21.1%, resulting in an average vanadium valence state of + 4.83 (in agreement with the literature data of 4.8 [29]). In contrast, VOH has a V⁴⁺/V⁵⁺ ratio of 6.8%, estimating a vanadium valence state of + 4.94. The substantial increase in the proportion of V⁴⁺ ions after TMPA⁺ introduction promotes the hopping of unpaired electrons between V⁴⁺ and V⁵⁺, enhancing the electronic conductivity and electrochemical kinetics of TMPA-VOH. An intriguing feature from Fig. 2d is the full width at half maximum (FWHM) of the V 2p_3/2 peak in TMPA-VOH (1.23 eV), which is broader than that of VOH (1.14 eV). The V 2p width of vanadium oxides depends on the population of the d band; VO_x with more unpaired electrons exhibits a broader line [36,38,39]. The broadening of the V 2p line for TMPA-VOH suggests a higher concentration of electrons in the d-band, further illustrating the significant increase in V⁴⁺ ions after the introduction of TMPA⁺. The O 1s spectra of TMPA-VOH were deconvolved into three components at 530.5, 531.4, and 532.5 eV (Fig. 2d). The dominant peak at 530.5 eV arises from the V-O lattice [36,37]. The latter two components can be assigned to C-O and H-O groups, coming from the surface contamination on the sample [31,36], as also detected in the C 1s spectra (Fig. 2b) [37]. In the O 1s region of VOH, an additional peak at 533.7 eV is observed, related to the interlayer water molecules [31,38]. Notably, the corresponding peak is absent in TMPA-VOH, suggesting that the introduction of TMPA⁺ cations not only intercalate into the interlayers but also expels the interlayer water molecules, consistent with the TG results.

The valence band (VB) spectrum of TMPA-VOH and VOH is presented in Fig. 2e. Both samples exhibit a dominated band at ~ 6 eV (from O 2p contribution) and a low intensity band at ~ 1.5 eV (from V 3d contribution) [36]. The O 2p component is associated with the emission from the vanadyl oxygen (center peak), chain oxygen, and bridge oxygen in the VO_x structure, hybridized with V 3d. The appearance of the V 3d component, which will not be observed in the orthorhombic V₂O₅, confirms the presence of V⁴⁺, attributed to its higher electron density in the conduction band [36,40]. In comparison with VOH, the valance band of TMPA-VOH exhibits a slight shift to lower energy. This indicates increased V-O interactions in TMPA-VOH, possibly arising from the interaction of TMPA⁺ with the VO_x interlayer surface, resulting in a subtle change in V-O bond length or angle [41].

Figure 3 displays the FTIR spectra of TMPA-VOH and VOH, with corresponding frequencies and assignments listed in Table S1. In the TMPA-VOH spectra, distinctive peaks in the 1550-1200 cm⁻¹ region (Region 1 in Fig. 3a, b) and 1030-1180 cm⁻¹ region (Region 2 in Fig. 3a and Fig. S6) correspond to characteristic vibration modes of -N-CH₃, -CH₃, and the C-H ring in the TMPA⁺ salt [42], affirming the presence of TMPA⁺ in the TMPA-VOH lattice. In Fig. 3a, c, peaks observed at 742, 975, and 1007 cm⁻¹ in TMPA-VOH represent the characteristic bands of the vanadium oxide framework [43,44]. The 738 cm⁻¹ peak in VOH is related to the asymmetric stretching of V-O-V bridging bands, with the shift to a higher wavenumber in TMPA-VOH (742 cm⁻¹) indicating a stronger V-O-V bond in TMPA-VOH [43]. The symmetric stretching of V-O-V observed at 586 cm⁻¹ in VOH is not detected in TMPA-VOH [45]. In the enlarged spectra (Region 3 in Fig. 3a, c), the broad peak near ~ 1010 cm⁻¹ in VOH stems from the symmetric stretching of the terminal oxygen V = O bond. In TMPA-VOH, this V = O bond splits into two peaks at 1007 and 975cm⁻¹, possibly corresponding to the stretching vibrations of V⁵⁺ = O and V⁴⁺ = O, respectively [43,44,46,47]. The higher concentration of V⁴⁺ ions in TMPA-VOH (V⁴⁺/V⁵⁺ ratio of 21.1%) results in a significant difference between the two signals, leading to a pronounced splitting of the V = O peak. In contrast, VOH has a lower proportion of V⁴⁺ ions (V⁴⁺/V⁵⁺ ratio of 6.8%), making the splitting of the V = O peak less obvious. This phenomenon of separated V = O peaks has been observed in other systems, such as pyridine-V₂O₅·nH₂O and Nitrogrn-doped V₂O₅, where the authors attributed it to the stronger interactions between PyH⁺ and V = O and the V = O-H₄N⁺, respectively [48,49]. In the spectrum of VOH, peaks at 3565 and 1613 cm⁻¹ are related to the stretching of O-H and bending of H-O-H in the structural water molecules, respectively [43,47]. Notably, similar peaks are not well recognized in TMPA-VOH, suggesting the expulsion of most interlayer water upon TMPA⁺ pre-intercalation, in agreement with the TG results.

Electron paramagnetic resonance (EPR) spectroscopy was employed to provide direct evidence for the existence of paramagnetic V⁴⁺ species, whose outer shell electronic structure 3d¹ has one unpaired electron to interact with the magnetic field. The EPR spectra of TMPA-VOH exhibit a prominent peak centered at a g-factor of 1.9728 (with a peak-to-peak line width of ~ 200 G), confirming the existence of V⁴⁺ ions (Fig. S7) [50]. In Fig. 3d, both TMPA-VOH and VOH display a symmetric signal at a g value of 2.003, originating from unpaired electrons captured at the oxygen sites (oxygen free radicals ${{O}_{2}}^{\bullet -}$) [51]. The signal intensity of TMPA-VOH is notably stronger than that of VOH, indicating a higher content of unpaired electrons in TMPA-VOH and emphasizing the formation of more low-valance state V⁴⁺ ions following TMPA⁺ intercalation.

To evaluate the impacts of TMPA⁺ pre-intercalation on Zn²⁺ ion storage properties, we assembled coin cells using TMPA-VOH or VOH as the working electrode, metallic zinc foil as the counter electrode, glass fiber as a separator, and 3 M Zn(CF₃SO₃)₂ as the electrolyte. Figure 4a compares the 1st galvanostatic charge-discharge (GCD) curves of TMPA-VOH and VOH at a current density of 0.1 A g^-1. TMPA-VOH demonstrates a high open circuit voltage (OCV) of 1.58 V and an initial specific capacity of 451 mAh g^-1, surpassing that of VOH (OCV: 1.31 V, discharge capacity: 399 mAh g^-1). Additionally, the mid-point voltage difference of TMPA-VOH is 0.14 V, significantly lower than that of VOH (0.22 V), indicating a reduced discharge/charge overpotential for TMPA-VOH. The rate performance of two samples is illustrated in Fig. 4b. TMPA-VOH delivers discharge capacities of 403, 382, 355, 326, and 294 mAh g^-1 at current densities of 0.5, 1, 2, 4, and 8 A g^-1, respectively, demonstrating an outstanding rate performance. In contrast, the discharge capacity of VOH decreases from 354 to 184 mAh g^-1 with the current density increasing from 0.5 to 8 A g^-1, resulting in a capacity retention of 52%, much lower than that of TMPA-VOH (73%). The corresponding GCD curves of TMPA-VOH and VOH are depicted in Fig. 4c. The narrower capacity gap between each change in current density and the smaller voltage gap between charge/discharge plateaus (e.g., the mid-point voltage difference at 0.5 A g^-1 is 0.15 V for TMPA-VOH and 0.27 V for VOH) imply better tolerance to high rates and low polarization of TMPA-VOH, possibly arising from its promoted charge transfer and reaction kinetics.

A comparison of the cycling stability between TMPA-VOH and VOH is present in Fig. 4d. At a current density of 4 A g^-1, TMPA-VOH exhibits a superior initial discharge capacity of 329 mAh g^-1, outperforming VOH (264 mAh g^-1). After 2000 cycles, TMPA-VOH maintains a remarkable capacity of 290 mAh g^-1, exhibiting a high retention of 87% (of the highest capacity). In contrast, VOH experiences rapid capacity decay, retaining only 162 mAh g^-1 (61% retention) after cycling. Figure 4e compares the corresponding GCD profiles of the 1st and the 2000th cycle at 4 A g^-1. Following cycling, the degradation of VOH is evident not only in the reduced specific capacity but also in a significantly enlarged overpotential. Another notable observation from Fig. 4e is that the capacity decay in both samples primarily originates from the lower voltage plateaus (V⁴⁺/V³⁺), suggesting that the insertion of Zn²⁺ into the low-energy redox sites cannot be fully achieved. This difficulty arises from the gradual collapse of the V-O structure, driven by strong electrostatic interactions between divalent Zn²⁺ ions and the lattice. For TMPA-VOH, the introduction of weakly polarized organic cations serves to widen migration channels and weaken electrostatic interactions, promoting Zn²⁺ diffusion kinetics and preventing obvious structural deterioration and ensuring superior cycling stability.

Energy efficiency, a crucial metric calculated as the ratio of discharge specific energy to charge specific energy, provides a quantifiable measure of the energy consumed or released during the extraction or intercalation of working ions [52]. Charging and discharging at different rates not only impact specific capacity but also influence redox overpotential (Fig. 4c). Therefore, plotting energy efficiency as a function of rates offers a more insightful perspective on the practical rate capability of the electrode. In Fig. 4f, TMPA-VOH consistently demonstrates higher energy efficiency than VOH across all current densities. Notably, as the current density increases from 0.1 to 8 A g^-1, TMPA-VOH only experiences a moderate drop in energy efficiency, from 89% to 66%, corresponding to an impressive retention of 74%. In contrast, VOH shows a more substantial decay from 86% to 48%, with an energy efficiency retention of 56%.

The superior energy efficiency and its slower decline at higher charge/discharge rates in TMPA-VOH indicate lower overpotential and enhanced rate capability with TMPA⁺ pre-intercalation. This improved reaction kinetics in TMPA-VOH may stem from: (1) its larger diffusion channels, facilitating a lower energy barrier for intercalation ions to reach redox sites; (2) its higher content of low valance state vanadium, providing more unpaired electrons for the conductivity of electroactive materials. Figure 4g presents a Ragone plot comparing the specific energy vs. specific power of our Zn//TMPA-VOH and Zn//VOH cells with recently reported aqueous ZIBs. Remarkably, the Zn//TMPA-VOH cell achieves the highest specific energy of 316 Wh kg^-1 at a specific power of 71 W kg^-1 and a specific energy of 197 Wh kg^-1 at a high specific power of 4925 W kg^-1. These values not only surpass those of the Zn//VOH cell in our study, but also outperform/match many other guest species-intercalated vanadium oxides, including PANI-VOH [53], NaCa_0.6V₆O₁₆·3H₂O [54], ([N(CH₃)₄]_0.77,Zn_0.23)V₈O₂₀·3.8H₂O [55], N-doped V₂O₅ [49], Na_0.33V₂O₅ [56], and Ca_0.24V₂O₅⋅0.83 H₂O [57].

In Fig. 5a, a comparison of cyclic voltammetry (CV) curves at a scan rate of 0.1 mV s⁻¹ is presented for TMPA-VOH and VOH. VOH exhibits two distinct pairs of well-defined redox peaks at 0.98/1.06 and 0.45/0.65, representing the redox reactions between V⁵⁺/V⁴⁺ and V⁴⁺/V³⁺, respectively [13,14,22]. Surprisingly, TMPA-VOH displays five couples of redox peaks: cathodic peaks centered at 1.33, 0.98, 0.91, 0.61, and 0.43 V, along with respective anodic peaks at 1.35, 0.99, 1.06, 0.67, and 0.51 V. This observation is further supported by the differential capacity (dQ/dV) vs. voltage curves at 0.5 A g⁻¹ (Fig. S8), revealing five redox couples in TMPA-VOH and aligning well with the CV results. While multiple redox peaks have been reported in various vanadium oxides [58,59,60], the precise assignment of these peaks remains unclear in the literature.

Firstly, from the initial three CV cycles of TMPA-VOH at 0.1 mV s⁻¹ (Fig. 5b), it becomes evident that the five couples of redox peaks in TMPA-VOH remain consistently present and nearly overlapped throughout subsequent cycles. This result strongly suggests the good reversibility of these redox reactions, effectively ruling out the possibility of irreversible phase transitions or structural changes. Secondly, as illustrated in Fig. 5c, with an increase in the scan rate, two pairs of redox peaks around 0.5 V begin to merge into a single pair, so do the two pairs around 1 V. In other words, while TMPA-VOH exhibits five distinct redox pairs at low scan rates, it manifests only three redox pairs at high scan rate. Based on these observations, it is likely that at the small scan rate of 0.1 mV s⁻¹, both pairs of redox peaks at 0.98/0.99 V and 0.91/1.06 V can be assigned to the V⁵⁺/V⁴⁺ redox reaction, while both pairs at 0.61/0.67 V and 0.43/0.51 V represent the V⁴⁺/V³⁺ redox centers. This split of vanadium redox peaks could be attributed to the existence of two distinct vanadium coordination in TMPA-VOH. During discharging/charging, the oxidation state of vanadium undergoes transitions between 5 + , 4 + , and 3 + , coinciding with the transformation of the corresponding vanadium coordination polyhedra [58,61]. In TMPA-VOH, the existence of two distinct polyhedral coordination of vanadium is evident from the separate V = O peaks (V⁵⁺ = O and V⁴⁺ = O) in the FTIR spectra (Fig. 3c). Consequently, the transformation of these two types of vanadium coordination occurs at distinct redox voltages, and this becomes more apparent at low scan rates, where the redox reaction is relatively slow. Unlike VOH, where this redox peak split is undistinguishable in CV curves, consistent with the single or nearly merged V = O peaks in FTIR spectra. The additional tiny peak at 1.33/1.35 V suggests the emergence of a new redox site, with some attributing this to the interaction between the empty orbit of Zn²⁺ and the lone pair electrons of N [49].

Further insights can be derived from the CV curves at 0.1 mV s⁻¹ (Fig. 5a). The voltage gap between each redox pair of TMPA-VOH is smaller than that of VOH, as detailed in Table S2). This indicates reduced polarization and accelerated reaction kinetics in TMPA-VOH, aligning with its diminished overpotential as evident from the reduced mid-point voltage difference in the GCD curves. The area enclosed by the CV curves corresponds to the charge storage in the material. Notably, the peak area under the V⁵⁺/V⁴⁺ peaks is nearly the same for both samples. However, the peak area under the V⁴⁺/V³⁺ pairs is significantly larger in TMPA-VOH than in VOH. This indicates that TMPA-VOH possesses more accessible redox sites in its structure, facilitating the insertion of more Zn²⁺ ions to reduce vanadium ions to V³⁺. This finding is consistent with the higher specific capacity of TMPA-VOH obtained from the GCD tests.

CV tests were conducted at various scan rates (ranging from 0.1 to 1.2 mV s⁻¹) to reveal the Zn²⁺ ion storage kinetics in TMPA-VOH (Fig. 5c). With increasing scan rates, the voltage difference between each redox pair expands, causing reduction peaks to shift to lower voltage and oxidation peaks to shift to higher voltage. In comparison to VOH (Fig. S9), TMPA-VOH exhibits a less significant enlargement in peak separation, indicating its reduced electrochemical polarization. b values, reflective of the dominant charge storage process, were derived from the log (peak current) vs. log (scan rate) plot (Fig. 5d, calculation details in Supplementary Information). For TMPA-VOH, the b values of peak O₁, R₁, O₂, R₂, O₃, and R₃ are 0.77, 0.96, 1.0, 0.84, 0.97, and 0.96, respectively. These values suggest that the Zn²⁺ storage in TMPA-VOH involves a combination of diffusion-controlled behaviors and surface-controlled capacitive behaviors [62].

As summarized in Fig. 5e, f, the contribution from each process can be quantitatively distinguished. Firstly, the capacitive contribution in TMPA-VOH at all sweep rates surpasses that in VOH, indicating an overall faster reaction kinetics of TMPA-VOH. This enhanced pseudocapacitive behavior may be attributed to the large surface area and the wide interlayer spacing of TMPA-VOH. Secondly, since capacitive charge storage is sweep rate independent (Fig. 5f), decreases in diffusion-controlled charge storage can serve as an effective indicator for evaluating the rate resistance of electrodes. As shown in Fig. 5e, the contribution ratio of the diffusion-controlled process in TMPA-VOH declines less (from 30 to 11% with scan rate increasing from 0.1 to 1.2 mV s⁻¹) than that of VOH (from 46 to 20%), indicating a superior rate capability of TMPA-VOH. Thirdly, from Fig. 5f, the specific capacitance/stored charge in TMPA-VOH (e.g., 1180 C g⁻¹ at 0.1 mV s⁻¹) surpasses that of VOH (925 C g⁻¹ at 0.1 mV s⁻¹), consistent with their specific capacify values from GCD tests. Additionally, as scan rates increase, the slower decay of specific capacitance in TMPA-VOH aligns with its improved rate performance (Fig. 4b). In summary, the kinetics analysis from CV data reveals that the charge storage in TMPA-VOH primarily results from the surface-controlled capacitive processes. The pre-intercalation of TMPA⁺ cations enhances pseudocapacitive behavior and results in a slower decay in diffusion-controlled charge storage with increasing scan rate.

Figure 6a compares the Nyquist plots of TMPA-VOH and VOH before and after CV tests. Both spectra consist of a semicircle in the high frequency region, associated with the charge transfer resistance (R_ct) in electrodes, and a sloped line at low frequency, related to the ion diffusion process. Before cycling, the R_ct of TMPA-VOH (28 Ω) is notably half that of the pristine VOH electrode (57 Ω), highlighting improved electronic conductivity resulting from organic cation pre-intercalation. This enhancement is likely attributed to the increased density of unpaired electrons in TMPA-VOH, supported by the significantly higher amount of V⁴⁺ in TMPA-VOH (21.1%) compared to VOH (6.8%). After cycling, R_ct decreases for both TMPA-VOH and VOH, reaching 13 and 26 Ω, respectively. The substantial reduction in R_ct after cycling may be ascribed to the small amount of residual Zn²⁺ in the structure, boosting the electrical conductivity of the cathode. Figure 6b illustrates the real part of impedance (Z’) as a function of frequency (ω^−1/2) in the low-frequency region. The line for the TMPA-VOH before cycling exhibits a smaller slope than VOH, indicating faster Zn²⁺ ion transport within TMPA-VOH. These EIS results underscore the active role of TMPA⁺ introduction in enhancing the electrochemical reaction kinetics, particularly by boosting charge transfer conductivity.

The Galvanostatic intermittent titration technique (GITT) was employed to compare the ion diffusion kinetics of TMPA-VOH and VOH at various charging/discharging stages. The calculated Zn²⁺ ion diffusion coefficients (D_Zn²⁺) during the 3rd GITT cycle at 50 mA g⁻¹ are presented in Fig. 6c (calculation details in Supplementary Information). TMPA-VOH exhibits D_Zn²⁺ values ranging from 9.1 × 10⁻¹¹ to 2.4 × 10⁻⁹ cm² s⁻¹ during discharging and 1.9 × 10⁻¹⁰ to 2.0 × 10⁻⁹ cm² s⁻¹ during charging. In contrast, VOH displays smaller D_Zn²⁺ values, approximately 8.6 × 10⁻¹² ~ 1.5 × 10⁻⁹ cm² s⁻¹ during discharging and 8.1 × 10⁻¹¹ ~ 1.0 × 10⁻⁹ cm² s⁻¹ during charging. Compared with VOH, the faster Zn²⁺ ion diffusion rate of TMPA-VOH throughout the entire intercalation/de-intercalation process indicates that the pre-insertion of TMPA⁺ into V-O layers effectively benefits the transport of charge carriers. As demonstrated in Fig. 6c, both systems exhibit relatively stable diffusivity at the beginning of the discharge process, followed by a rapid decrease after 1.0 V. Initially, the layered structure with large channels allows the “free” insertion of Zn²⁺ ions; afterwards, deep intercalation is hindered due to the increasing electrostatic interactions between the inserted Zn²⁺. However, this diffusivity decay is more significant in VOH than in TMPA-VOH, especially at lower voltage. This supports the previous statement derived from the CV curves that TMPA-VOH can accommodate more Zn²⁺ ions, inserting into the deeper lattice to reduce more V⁵⁺/V⁴⁺ to V³⁺. Another important factor, internal resistance, can be estimated from the IR-drop in the GITT results (Fig. 6d). Internal resistance represents the overall resistance in electrochemical reactions, originating from active materials, current collectors, conductive additives, contact between each element, and electrolyte conductivity. Throughout the charge/discharge process, the internal resistance values for TMPA-VOH are consistently lower than those of VOH (Fig. 6d), contributing to a reduced overpotential and higher energy efficiency for TMPA-VOH, as confirmed by previous GCD and CV results. During discharging (Fig. 6d), the internal resistance of VOH displays two steps of increase at ~ 1.2 and 0.6 V, coinciding with the two-step redox reactions upon the Zn²⁺ intercalation. The dramatic increase in resistance towards the end of discharge in VOH is a consequence of kinetic limitations, attributed to increased concentration polarization and electrochemical polarization. In contrast, the internal resistance of TMPA-VOH remains relatively stable throughout the entire discharge process, further indicating the reduced electrostatic interactions between the inserted Zn²⁺ and the lattice, facilitating the insertion of more Zn²⁺ ions into the TMPA-VOH lattice.

Benefiting from the 3D flower-like morphology, a wide interlayer spacing, and a large amount of V⁴⁺ with more unoccupied electrons, TMPA-VOH exhibits an accelerated ion diffusion rate and improved charge transfer kinetics. The low polarity of the pre-inserted organic ions minimizes interactions with the cycled Zn²⁺ ions, resulting in reduced polarization. Moreover, the pre-intercalation of organic cations has a significant impact on preventing the decay in diffusion rate and the increase in internal resistance upon deep discharging, enabling relatively stable kinetics throughout the entire discharging/charging process.

To elucidate the energy storage mechanism of the Zn//TMPA-VOH battery, the structural evolution, valance state and morphological changes of TMPA-VOH during charge/discharge were investigated through ex situ XRD, XPS and SEM analyses. Figure 7a displays the XRD pattern of TMPA-VOH electrodes in the pristine state, fully discharged state (0.2 V) and fully charged state (1.6 V) during the first cycle. The XRD pattern of the pristine electrode closely matches the as synthesized TMPA-VOH powder and the Ti current collector. Upon complete discharge to 0.2 V, there is no significant shift in the position of the (00l) peak, indicating minimal expansion or contraction of the interlayer distance during Zn²⁺ intercalation into the structure. This contrasts with some reported vanadium oxides that exhibit noticeable interlayer expansion, which was attributed to screened electrotactic repulsions between V-O layers caused by inserted Zn²⁺ [13]. Conversely, certain vanadium oxides demonstrate an obvious interlayer contraction explained by the electrostatic attractions between Zn²⁺ and the host structure, pulling the layers closer [34,55,63]. The negligible interlayer spacing variation in TMPA-VOH suggests alleviated electrostatic interactions between intercalated Zn²⁺ ions and the V-O lattice due to TMPA⁺ pre-insertion [27]. This reduction in electrostatic interactions is beneficial for Zn²⁺diffusion and structural stability during cycling.

It is noteworthy that no discernible new peaks are identified in the discharged electrode, which is different from observations in some studies where the generation of byproducts or secondary phases, such as Zn_x(CF₃SO₃)_y(OH)_2x-y·nH₂O or Zn₃(OH)₂V₂O₇·2H₂O, is commonly observed during the discharging process [53]. In the case of TMPA-VOH, characteristic peaks of Zn_x(CF₃SO₃)_y(OH)_2x-y·nH₂O or Zn₃(OH)₂V₂O₇·2H₂O phases cannot be clearly detected, possibly due to the low abundance of new phases generated during discharging, resulting in an XRD signal too weak to be identified [55]. Additionally, the main peaks corresponding to Zn₃(OH)₂V₂O₇·2H₂O (JCPDS No. 50-0570) appear at 12.3°, 20.9°, 30.1°, and 34.2° [64], closely located to the peaks of (002), (003), (40 $\overline{1 }$), and (11 $\overline{3 }$) planes of TMPA-VOH (at 12.6°, 19.2°, 30.8°, and 34.7°, respectively), making it even more challenging to differentiate the new phases from the main peaks. However, based on the XPS results, we speculate the formation or disappearance of these new phases upon discharging or charging (as discussed later). This phenomenon is similar to that observed in (Na,Mn)V₈O₂₀·nH₂O [65] and [N(CH₃)₄]_0.77,Zn_0.23) V₈O₂₀·3.8H₂O [55].

As shown in the Zn 2p XPS spectra (Fig. 7b), no signal is detected for the pristine electrode. When discharged to 0.2 V, strong Zn peaks clearly appear, confirming the intercalation of Zn²⁺ into TMPA-VOH. Upon charged to 1.6 V, the Zn signal is still present but with reduced intensity, suggesting the incomplete extraction of Zn²⁺ from the TMPA-VOH structure. These residual Zn²⁺ ions are believed to be trapped in the host lattice, contributing to enhanced electronic conductivity, as reflected by the reduced charge transfer resistance after cycling (as evidenced by EIS data). Notably, the discharged electrode shows two Zn components in the Zn 2p_3/2 region, a strong one at 1022.7 eV and a small side one at 1025.2 eV, which can be assigned to the intercalated Zn²⁺ in TMPA-VOH and the Zn precipitate (e.g., Zn_x(CF₃SO₃)_y(OH)_2x-y·nH₂O or Zn₃(OH)₂V₂O₇·2H₂O) generated on the surface [66,67]. After being fully charged, the latter peak disappears, with only the peak associate with intercalated Zn²⁺ showing at 1022.5 V, suggesting the reversible generation/dissolution of the zinc second phase.

The V 2p region is displayed in Fig. 7c. Consistent with the as-synthesized TMPA-VOH material (Fig. 2d), two peaks corresponding to V⁵⁺ (517.8 eV) and V⁴⁺ (516.1 eV) exist in the pristine electrode, with V⁵⁺ being the dominant one. After discharge, the signal of V⁴⁺ intensifies and becomes the major component, while that of V⁵⁺ decreases, indicating the reduction of vanadium ions accompanied with the intercalation of Zn²⁺. A shift in the binding energy of V in the discharged state could be attributed to the interactions between intercalated Zn²⁺ and the V-O layers [68]. A smaller shift of the V⁵⁺ peak in TMPA-VOH (0.1 eV) compared to other guest species pre-intercalated vanadium oxides, like Cr-VOH [69] or Polypyrrole-VOH [70], suggests reduced electrostatic interactions due to TMPA⁺ pre-intercalation. Interestingly, the V³⁺ signal was not observed in the discharged electrode. Although, based on the CV curves or the specific capacity of TMPA-VOH (451 mAh g^-1), there exists a valence conversion from + 4 to + 3 in vanadium. This phenomenon is widely observed for other vanadium-based cathode [22,71], as trivalent vanadium cations are very unstable and can be easily oxidized during sample handling (exposed to air) and ex situ XPS characterization. After complete charging, both the position and fraction of vanadium spectra return to the initial state, demonstrating the reversible redox reaction. The slightly higher amount of V⁴⁺ in the charged electrode than in the pristine electrode can be ascribed to the incomplete de-intercalation of Zn²⁺ in the first charging process, in agreement with the Zn 2p spectra analysis (Fig. 7b).

In Fig. 7d, the O 1s spectra of the pristine electrode reveal three distinct components at 530.4 V (attributed to lattice oxygen in the V-O matrix), 531.6 V (indicative of C-O from surface contaminations), and 532.8 eV (corresponding to surface O-H species). Upon complete discharge, an additional peak emerges at 533.8 eV, possibly associated with H₂O molecules. This may be attributed to the insertion of hydrated zinc ions or the formation of Zn_x(CF₃SO₃)_y(OH)_2x-y·nH₂O or Zn₃(OH)₂V₂O₇·2H₂O. The appearance of the H₂O peak, along with the strengthening of the -OH peak, indicates the possible formation of Zn_x(CF₃SO₃)_y(OH)_2x-y·nH₂O or Zn₃(OH)₂V₂O₇·2H₂O, consistent with the Zn spectra analysis. Upon charging, the diminishing -OH peak and H₂O peak suggest the decomposition of these newly formed phases. The -OH in the zinc complex arises from the dissociation of water molecules in the electrolyte, induced by the intercalation of H⁺ into the cathode [72]. Therefore, a Zn²⁺/H⁺ co-intercalation mechanism can be reasonably proposed for TMPA-VOH, aligning with many other vanadium oxide cathodes [34,55,65]. As demonstrated in the ex situ SEM images (Figs. 7d and S10), there is no significant morphological evolution of TMPA-VOH during charging and discharging. The homogeneous distribution of V, O, N, C, and Zn across the sample further confirms the insertion of Zn²⁺ into TMPA-VOH and the stability of TMPA⁺ cations during discharging/charging. The weak Zn signal detected at the fully charged state suggests the incomplete extraction of Zn²⁺ ions, which agrees with the XPS data. Based on the XRD, XPS, and SEM analyses, the energy storage mechanism of Zn//TMPA-VOH involves a Zn²⁺/H⁺ co-insertion. The pre-intercalation of TMPA⁺ cations significantly decrease the electrostatic interactions between inserted Zn²⁺ and the V-O lattice, resulting in negligible interlayer variation and structural evolution during discharging/charging. This phenomenon is responsible for the enhanced kinetics and stability observed during cycling.