Introduction

The increasing demand for grid energy storage in modern society continues to propel the global battery industry1,2,3. The transition to TWh-scale energy storage presents new challenges in battery performance, especially battery cost, sustainability, and intrinsic safety4,5,6. Lithium-ion batteries, with their high energy density (>150 Wh kg−1), currently dominate the market of energy storage batteries. However, the low identified resources of lithium and the flammability and toxicity of the organic electrolytes pose a critical challenge to sustainability. Aqueous batteries, compared to commercial lithium-ion batteries, offer significant advantages in cost-effectiveness, recyclability and safety, positioning them as strong contenders in grid energy storage7,8. And zinc metal has emerged as an ideal candidate for aqueous anodes due to its low redox potential, high theoretical capacity, and abundant availability, making aqueous zinc metal batteries (ZMBs) a promising alternative for grid energy storage9,10,11. ZMBs require high reversibility to meet the extended operational demands of grid energy storage. To illustrate, the CE for the plating and stripping of zinc must surpass 99.97% to facilitate 5000 cycles at 60% end-of-life for a ZMB with an N/P ratio of 2. However, challenges like corrosion12 and dendrite formation13 hinder achieving high CE and long lifespan, significantly impeding ZMBs development.

Numerous efforts have been made to enhance the reversibility of the Zn anode, encompassing electrolyte engineering14,15,16,17,18, crystal-facet manipulation19,20,21,22,23, various coatings24,25,26,27,28, and substrate optimization29,30,31. These advancements have greatly improved the cycle life of the Zn anode. However, achieving a near-unity CE, especially at current densities of 1 mA cm−2 and below, remains a significant challenge17. The difficulty in achieving near-unity CE arises from corrosion and HER issues in aqueous electrolytes, which are inherently tied to the low potential of Zn metal and controlled by thermodynamics. To make matters worse, dendrite formation can greatly exacerbate these side reactions. As a result, satisfactory CE (>99.9%) have remained elusive32. This CE challenge is the same for other metal anodes33,34. Concentrated electrolytes have been proposed to extend the electrochemical stability window of aqueous electrolytes, enabling high CE of the Zn anode9,35,36. However, these highly concentrated electrolytes face challenges such as low ionic conductivity and high viscosity, posing barriers to commercialization.

Given the distinctive plating/stripping chemistry of metal anodes, a singular approach is typically insufficient to achieving highly reversible metal chemistry. An integration of multiple approaches is necessary, yet their interaction must be carefully regulated to prevent adverse effects. Different from highly concentrated electrolytes, commercial lithium-ion batteries are built on a cocktail strategy of various solvents and additives, each contributing to the overall electrolyte performance. In this study, we explore electrolyte design principles and apply this cocktail strategy to Zn metal chemistry towards high CE at moderate concentrations. Our approach involves a dual-salt strategy for the purpose of promoting homogeneous Zn deposition through anion-rich Zn2+ solvation structures. The other aspect of our cocktail strategy is to identify by theoretical screening appropriate organic co-solvents to inhibit corrosion and hydrogen evolution through reduced water activity. These anion-rich solvation structures and organic co-solvents facilitate the formation of a solid electrolyte interphase (SEI). The resulting hybrid dual-salt electrolyte (HDE) renders high reversibility for Zn stripping/plating, with over 99.95% CE at 1 mA cm−2 with a moderate concentration (3.5 m). Furthermore, this HDE imparts high cycling stability to Cu | |ZnI2 anode-free batteries, maintaining more than 95% capacity after 1000 cycles under practical conditions.

Results

Design principles of cocktail strategy

The irreversibility of Zn plating/stripping process is mainly attributed to dead Zn (QdeadM), corrosion (Qcor) and hydrogen evolution reaction (HER, QHER). The irreversible aspect of QHER intricately is associated with the thermodynamic properties of electrolyte (specifically, water activity) and catalytic behavior of the anode-electrolyte interface. Qcor shares a similar dependence on water activity and interfacial properties. QdeadM, on the other hand, arises from Zn dendrites and is exacerbated by Qcor and QHER, playing a major role in calendar aging37,38. The irreversibility associated with Qcor and QHER is also closely intertwined with the inhomogeneous deposition of Zn, highlighting that these irreversible components are interconnected. Hence, achieving a high CE of the Zn anode necessitates both homogeneous Zn deposition and a relatively ‘inert’ anode-electrolyte interface. This inert interface may exhibit suppressed kinetics, such as the formation of a passivated interphase, or hindered thermodynamics, such as the reduction of water activity.

In response to the above irreversible aspects, we can take two key steps (Fig. 1a): (1) fine-tuning the Zn2+ solvation structures to achieve uniform zinc deposition, (2) customizing the inner Helmholtz plane to reduce water activity. Furthermore, the formation of a robust SEI is highly related to how we design the solvation structures and co-solvents. Concerning Zn2+ solvation structures, contact-ion pairs (CIP) solvation structures were found conducive to promote more homogeneous metal deposition by reducing the exchange current density in alkali metal batteries39. High donor number anion depicts a strong affinity between the cation and anion and are more likely to form CIP solvation structures. As shown in Fig. 1b, acetate (OAc-) shows much higher donor number among common non-corrosive anions. As discussed in the Figs. S1S3, OAc- is more likely to form CIP structures than the widely used sulfate and triflate (OTf-) ions, promoting uniform Zn deposition. But the former has the drawback of lower solubility. In contrast, the low-polarity triflate possesses higher solubility and a more hydrophobic double-layer structure40, diminishing corrosion and HER. This inspires us to explore a dual-salt approach.

Fig. 1: Design principle and the realization of cocktail strategy for Zn metal anode.
figure 1

a Illustration of cocktail strategy and hybrid dual-salt electrolyte design. b, c Rational selection of dual-salts and co-solvent. d The average CEs of various hybrid dual-salt electrolytes based on 10th~100th cycles.

To construct an ‘inert’ anode-electrolyte interface, a co-solvent approach is deployed. We first screened some common organic co-solvents by theoretical calculations. The free-water density and H-bond number of water in inner Helmholtz plane are calculated using classic molecular dynamics simulations (Fig. S4). Triethyl phosphate (TEP) and dimethylformamide (DMF) are selected as rational co-solvent because they can induce the free-water-poor inner Helmholtz plane (within 5 Å from the electrode) and reduce the interfacial water activity (Fig. 1c). And it has been shown that both CIP structures and organic co-solvent will support the SEI formation41,42. Therefore, a cocktail strategy combining a dual-salt and organic co-solvent approach is expected to promote uniform Zn deposition while inhibiting corrosion and HER, as confirmed in Fig. 1d. It should be noted that approaches of individual components, such as co-solvent engineering using DMF43,44,45 or TEP46,47,48,49, have been extensively applied in the design of ZMBs electrolytes. However, these approaches have typically resulted in CE below 99.8%, even when incorporating a high content of organic co-solvents (Supplementary Table 1). The unique plating/stripping chemistry of metal anodes makes their requirements for the electrolyte more stringent. Likewise, the concept of the cocktail strategy maintains its consistency across both aqueous and nonaqueous systems, with the goal of achieving enhanced comprehensive performance through the integration of multiple effects. The key difference lies in the management of these effects: in aqueous systems, the regulation of water activity is paramount, whereas in nonaqueous systems, the formation and control of the SEI take precedence.

Reversibility of Zn metal anode

Commencing with a plating/stripping test at 1 mA cm−2 and 1 mAh cm−2, we compare the CEs of three types of electrolytes, viz., 4 m ZnOTf2 (single salt electrolyte, SE), 4 m ZnOTf2 + 1 m ZnOAc2 (dual-salt electrolyte, DE), and 2.8 m ZnOTf2 + 0.7 m ZnOAc2 + 30 wt% DMF (hybrid dual-salt electrolyte, HDE). The concentrations of dual-salt and co-solvent were optimized through CE values (Fig. S5). The SE exhibits an average CE of 99.3% and short-circuit after less than 100 cycles (Fig. 2a). Issues related to side reactions and dead zinc has contributed to the low reversibility and CE (Figs. S6, S7), as well as poor calendar aging property (Fig. S8). Similarly, the Zn | |Zn symmetric cell with SE encounters quick short-circuit, primarily due to uneven Zn deposition (Fig. S9). Upon the introduction of the dual-salt electrolyte, the CE significantly improves, reaching an average of 99.8%. With further inclusion of co-solvents, DMF and TEP exhibit the highest CE among these tested co-solvents, proving the validity of theoretical screening (Fig. S10). The average CEs of HDE reaches 99.95% (based on an average of 1000 cycles) and stabilizes at approximately 99.97%. The other four individual cells utilizing HDE also exhibit similar average CEs, around 99.95% (Fig. S11). On the contrary, both the single salt electrolyte with ZnOAc2 and hybrid electrolyte with single salt of ZnOTf2 exhibit unsatisfactory CEs of <99.5% (Fig. S10). This achievement of HDE represents one of the highest reported CEs in the literature (Supplementary Table 1). Furthermore, both DE and HDE can enable stable cycling of Zn | |Zn symmetric cell for more than 2000 h (Fig. S9). HDE also facilitates stable cycling of zinc foils at a high depth of discharge (DOD) of 60% for over 600 h (Fig. S12).

Fig. 2: Coulombic efficiencies of different electrolytes.
figure 2

Zn plating/stripping CE with magnified views of selected cycles between 800th cycle to 9000th cycle in at (a) 1.0 mA cm−2 and 1.0 mAh cm−2 and (c) 0.5 mA cm−2 and 0.5 mAh cm−2. Zn plating/stripping potential profiles on a Cu electrode cycled in HDE at (b) 1.0 mA cm−2 and 1.0 mAh cm−2 and (d) 0.5 mA cm−2 and 0.5 mAh cm−2. e, f CE tests under ‘reservoir half-cell’ galvanostatic protocol. A fixed capacity of Zn (5 mA cm−2, 5 mAh cm−2) is firstly plated, followed by plating/stripping at a fixed DOD of 60% (200 cycles) and 80% (100 cycles).

In general, low current densities will amplify the effects of HER and corrosion. in our test at a current density of 0.5 mA cm−2, HDE still achieves average CEs of 99.90% and a stabilized CE of 99.93%. Furthermore, achieving high CEs at high areal capacities presents another challenge for practical applications. As shown in Fig. S13, our HDE can attain a CE of 99.9% within the initial 5 cycles and sustain this level of efficiency even at an areal capacity of 5 mAh cm−2 under a current density of 1 mA cm−2. Aurbach’s galvanostatic protocol was additionally employed to evaluate CE at high DOD. The CE approaches almost 100% at high DODs of 60% and 80%. In contrast, SE fails the test and DE shows a CE of only 99.7% (Fig. S14). Moreover, HDE offers notably low temperature ionic conductivity (3.2 mS cm−1 at −40 °C) and a wider electrochemical window (Fig. S15). HDE also shows a higher zinc ion transference number of 0.41, but it was still not sufficient (Fig. S16). It should be noted that this strategy is not applicable for harsh condition of very high current density and very high areal capacity (Fig. S13d). Fast deposition of zinc metal may no longer be dominated by CIP solvation structures, but rather by SSIP solvation structures with faster kinetics, thus leading to short circuits.

To demonstrate the adaptability of our cocktail strategy, we substituted the highly soluble triflate with high-concentrate ZnCl2 and obtained CEs exceeding 99.9% (Fig. S17). Moreover, we explored the substitution of the ZnOTf2 separately with the hydrophobic ZnTFSI2 and substitution with hydrophobic quaternary ammonium salts. The obtained CEs from these two electrolytes at 1 mA cm−2, 1 mAh cm−2 with a medium concentration exceed 99.9% (Fig. S18), validating the universality of our cocktail strategy.

Characterization of solvation structure and hydration structure

To investigate the origin of the high CEs in HDE, we conducted comprehensive characterization to understand the solvation structures and water environments in these electrolytes. Fourier transform infrared spectroscopy (FTIR) can reveal the evolution of water environments. As illustrated in Fig. 3a, the intensity of O-H stretching (3200–3600 cm−1) of H2O in SE at higher wavenumbers range decreases with the introduction of OAc- (Supplementary Table 2). Specifically, three hydrogen-bond states were identified in the water structure36,50: strong H-bond (3263.8 cm−1), weak H-bond (3247.3 cm−1), and none H-bond (3549.2 cm−1). The decrease in strong H-bond content in DE and HDE is related to the zinc ion concentration and the presence of the DMF co-solvent. The solvation sheath of zinc ions as well as DMF co-solvent disrupt the hydrogen-bonding network between water molecules, thereby reducing water reactivity. Furthermore, gas chromatography tests (Fig. S19) demonstrated that DE and HDE has effectively inhibited HER.

Fig. 3: Characterization of solvation structures and H2O structure.
figure 3

a FTIR spectra of various electrolytes at the wavenumber range of O-H stretching vibration. b, c Raman spectra in the wavenumber region of the SO3 stretching modes and symmetric C=O stretching modes, respectively. d Strong H-bond contents and CIP solvation structures contents of various electrolytes. e The XANES spectra and (f) EXAFS spectra in R space of various electrolytes. g MD simulation snapshot and typical solvation structures of HDE, and distribution percentage of possible main solvation sheath components of various electrolytes.

Raman spectroscopy was employed to reveal the solvation structures of Zn2+ (Fig. 3b–d). The −SO3 stretching modes (1000–1080 cm−1) of OTf- were utilized to identify OTf- involved solvation structures42. Limited reports are available regarding the OAc- involved solvation structure identification. A blue shift of νs(C = O) is observed with an increased content of OAc- (Fig. S20), indicating that solvated OAc- may introduce a blueshift in νs(C = O), similar to the trend of the −SO3 stretching modes of OTf-. This means that νs(C = O) can be used to identify OAc- involved solvation structures. Moreover, the Raman peaks were deconvoluted into solvent-separated ion pair (SSIP), contact ion pair (CIP), and aggregate (AGG). In SE, the proportion of CIP (1042.3 cm−1) is approximately 36.9% (Supplementary Table 3), suggesting that the dominant solvation structure contains (Zn-(H2O)6)2+. Intriguingly, the coordination number of OTf is slightly reduced after the introduction of acetate, with acetate entering in the coordination of Zn2+ in DE (64.5%, 1448.3 cm−1). This indicates that acetate is more competitive for Zn2+ solvation structures in the dual-salt42, replacing part of the OTf solvation, in line with the solvation energy trend (Fig. S21), viz., (Zn-OAc-(H2O)5)+ (20.5 eV) > ((Zn-OTf-(H2O)5)+ (19.8 eV). The introduction of a significant amount of acetate in the sheath also led to the appearance of AGG (12.4%, 1473.1 cm−1) in DE. In HDE, the introduction of DMF increased the coordination number of OTf. The case of OAc- becomes more complicated and it is difficult to obtain the exact proportions in HDE. Nonetheless, the proportion of OAc- is significantly decreased compared to DE (Fig. S20). The decrease in AGG is favorable for the increase in electrolyte ionic conductance, especially at low temperatures.

The solvation structures were further studied by X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) analysis of Zn K-edge (Fig. 3e, f). Compared with 1 m ZnOAc2 electrolyte, the 4 m ZnOTf2 exhibits lower edge energy in the XANES spectrum (Fig. 3e). This difference may be attributed to variations in the electron-donating ability between OTf- and OAc-. And the ESP analysis indicates the charge of Zn2+ in OTf- involved CIP (+1.116) is lower than that in OAc- involved CIP (+1.270). In the DE, the presence of more CIP structures leads to a shift of the K-edge energy to lower values. The introduction of DMF induces further alterations in the solvation structure. As a result, the K-edge shift to higher energies due to a reduction in the CIP content and the entry of DMF into the sheath. EXAFS spectra in R space confirm the rearrangement of the local solvation structure upon the introduction of DMF (Fig. 3f).

Molecular dynamics (MD) simulations were performed to provide additional insights into the structure of these electrolytes (Figs. 3gS22S25). The evolution of water environments influences the order of the H-bond network, resulting in less ice-like water. The tetrahedral order parameter of water moves towards a more disordered state, impacting the low-temperature ionic conductance of the electrolyte. Consequently, HDE exhibits superior low-temperature ionic conductance. The solvation structures are also confirmed in the simulations. In SE, SSIP is the primary component. The SSIP content decreases to 50.4% in DE, with OAc- entering the solvation sheath from a larger proportion. A high proportion of AGG in DE adversely affects ionic conductance (Fig. S17). In HDE, SSIP content is further reduced, and the proportion of AGG in DE also decreases. The simulations also confirm that DMF enters the solvation sheath of Zn2+ with a coordination number of 0.19 (Fig. S24). All these results imply that the dual-salt introduces more anions into the Zn2+ primary solvation sheath, and the DMF co-solvent is also able to enter the solvation sheath.

Zn plating/stripping morphology and SEI chemistry

We assessed the morphology of plated Zn and Zn after soaking in different electrolytes. The current density reports an important role in deposition morphology. In DE and HDE, high current density leads to more uneven deposition morphology. Nevertheless, the elevated CIP content in DE and HDE leads to more uniform Zn nucleation and smoother deposition under plating currents of 1 and 5 mA cm−2 (Figs. S26, S27). Moreover, the HDE results Zn surface with more exposed (002)Zn planes, which also contributes to the higher stability of Zn metal anode (Fig. S28). This uniform Zn nucleation and deposition was further verified by chronoamperometry (CA) (Fig. S29). The variation of current at constant overpotential correlates with the nucleation process. In SE, the current continuously increases within 600 s, which represents an uneven nucleation process and rough three-dimensional propagation. HDE and DE, on the other hand, both exhibit two-dimensional like growth mode.

Immersion tests demonstrate that, compared to SE and DE, corrosion is significantly reduced in HDE (Figs. S30S32). Dendrite formation and corrosion byproduct are obvious in SE, but much less on DE and nearly absent in HDE (Figs. 4aS33, S34), highlighting the superior anti-corrosion property rendered by HDE. The water-deficient inner Helmholtz plane contributes to the corrosion-resistant interface (Figs. S35, S36). In the in-situ optical microscopy observation, the plated Zn becomes rough with dendrite formation after 20 minutes of plating with clear gas bubble formation in SE (Fig. 4d). Conversely, both DE and HDE electrolytes enable flat and compact Zn plating (Figs. 4eS37). Surface 3D reconstructions of plated Zn also confirm the much more homogeneous Zn deposition with DE and HDE (Fig. 4c).

Fig. 4: The morphology of plated Zn.
figure 4

SEM images of Zn anode after 50 cycles at 1 mA cm−2 and 1 mAh cm−2 in (a) SE and (b) HDE. Scale bar: 30 μm in a and 2 μm in (b) (magnified view). c 3D height images of plated Zn in SE and HDE with a capacity of 5 mAh cm−2 at 5 mA cm−2. In-situ optical microscopy images of Zn plating in (d) SE and (e) HDE at a current density of 5 mA cm−2.

CIP structures in aqueous batteries, particularly those containing fluorine-containing anions, have been reported to promote dense SEI formation. This is consistent with organic co-solvents. We examined the surface structure of cycled Zn in HDE using X-ray photoelectron spectroscopy (XPS), time-of-flight secondary ion mass spectrometry (ToF-SIMS), and transmission electron microscopy (TEM). XPS tests show a predominance of oxides in the SEI (Figs. S38, S39), which is consistent with the results of the soaking tests. As shown in Fig. 5a, the peaks at approximately 688.6 eV and ~684 eV in the F 1 s spectrum are attributed to *CF3 and ZnF2, respectively. After Ar+ sputtering for 100 s, the F 1 s peak is predominated by ZnF2, indicating the presence of a ZnF2-contained SEI structure. Theoretical calculations (Fig. S40) and XPS tests (Fig. S41) reveal that ZnF2 is derived from solvated OTf- ion. The C 1 s profiles in Fig. 5b indicate the presence of some ZnCO3 (~289.1 eV) on the top surface of Zn, which can be attributed to the decomposition of acetate ions and DMF (Fig. S42). ToF-SIMS results further confirm the presence of ZnF2 in the SEI, while ZnCO3 is found in only small amounts at the top surface (Fig. 5c). TEM tests reveal that the SEI thickness ranges from approximately 50 to 100 nm and is rich in inorganic components (Fig. S43). This robust SEI plays a crucial role in passivating the Zn anode and suppressing side reactions (Fig. 5d).

Fig. 5: Characterizations of SEI on cycled Zn.
figure 5

Depth-resolved XPS of (a) C 1 s spectra, and (b) F1s after 50 cycles in HDE. c TOF-SIMS color mappings of ZnF, CO32−, and ZnO fragments. Test area is 50 µm × 50 µm. d Schematic presentation of proposed inorganic rich SEI.

Full-cell performance

We further examined the reversibility of full cells in HDE, including Zn | |ZnI2 cell, Zn | |vanadium oxide full cell, and Zn | |activated carbon (AC) cell. As illustrated in Fig. 6a, the SE supports only around 200 cycles of Zn | |ZnI2 cell at a NP ratio of ~0.8. In contrast, HDE can prolong it to more than 3000 cycles without evident capacity degradation, maintaining approximately 97% of its initial capacity (~142.7 mAh mg−1). Such near-unity CE of HDE can enable anode-free Cu | |ZnI2 batteries to cycle more than 1000 times with minimal capacity loss, preserving more than 95% of its initial capacity. Additionally, the anode-free Cu | |ZnI2 cell under practical conditions, high areal mass loading (>20 mg cm−2) and lean electrolyte (~10 µL mAh−1), retains over 80% of its initial capacity after 300 cycles at 0.15 A g−1 (Fig. 6c). Zn | |ZnI2 cell with NP ratio about 1.0 demonstrate a similar decay rate with this anode-free cell (Fig. S44), suggesting the iodine dominates the cell capacity degradation. When a CCCV charge protocol (1.6 V/50 mA g−1) was employed, the anode-free Cu | |ZnI2 cell can be stably cycled at a high mass loading for more than 1000 cycles without significant capacity degradation. This anode-free Cu | |ZnI2 cell exhibits good rate performance and low self-discharge characteristics (Figs. S45, S46). The specific energy density of this anode-free cell reaches 82 Wh kg−1 based on the total mass of electrodes. Further introduction of Br or Cl increases the energy density to 113 and 145 Wh kg−1, respectively (Fig. S47).

Fig. 6: Electrochemical performance of full cells.
figure 6

a Anode-less cell cycling with a ZnI2 cathode (mass loading of 8.2 mg cm−2) at 1 A g−1. b Cycling performance of anode-free Cu | |ZnI2 cell. c, d Anode-free Cu | |ZnI2 cells and Zn | |V2O5 cells cycling performance under practical conditions (low current density and high mass loading of active materials). The ZnI2 content in the CCCV test is 50 wt% and the rests are 43.2 wt%.

High mass loading vanadium oxide cell also demonstrates stable cycling for 300 cycles under practical conditions (Figs. 6dS48). In addition, due to the high ionic conductivity and non-flammable nature of HDE (Fig. S49), the Zn | |V2O5 cell remains stable cycling from −40 °C to 50 °C (Figs. S50S52). We also validated the stable cycling performance of large anode-free Cu | |ZnI2 pouch cells (Fig. S53).

Discussion

We have demonstrated a cocktail strategy to design hybrid electrolytes for zinc metal batteries. Integrating a high donor number acetate salt and a highly soluble triflate salt results in a dual-salt electrolyte with a solvation structure rich in CIP and a hydrophobic IHP structure. The dual-salt configuration promotes uniform zinc deposition while suppressing corrosion and hydrogen evolution side reactions. The simultaneous addition of organic co-solvents, such as DMF or TEP, further reduces water activity and facilitates the formation of a robust SEI. As a result, this HDE renders highly reversible Zn plating and stripping, achieving a high CE of 99.95% at a current density of 1 mA cm−2. We have conducted extensive characterizations to delve into the fundamental properties of this hybrid dual-salt system. As proof-of-concept demonstration, we have implemented such hybrid system in an anode-free Cu | |ZnI2 full cell, which shows long-cycle stability under practical conditions.

The cocktail strategy presented here offers a rational and promising approach to achieving highly reversible aqueous metal batteries chemistry. For the case of MnO2 cathode (Fig. S54) which is chemically unstable, one may consider introducing additional effects into the cocktail strategy to achieve highly stable ZMBs with high specific energy density. Finally, similar to zinc, most metal-based batteries, such as aqueous manganese- and aluminum-metal batteries, also face challenges related to dendrite formation and side reactions, which will deteriorate battery performance and safety. By leveraging the distinct features of different salts and solvents, this cocktail strategy may guide the design of hybrid electrolyte systems to address these challenges.

Methods

Chemicals

Zinc triflate (ZnOTf2, >99%), Zn sulfate (ZnSO4*7H2O, AR), anhydrous zinc chloride (ZnCl2, >99%) and Zn acetate (ZnOAc2, AR) were purchased from Adamas, Sinopharm, Sigma-Aldrich and Macklin, respectively. Dimethylformamide (DMF, >99.8%) and triethyl phosphate (TEP, >98%) were purchased from BLDpharm and Sigma-Aldrich, respectively. Zn metal foils were purchased from Hebei Qingyuan Metal Material Technology Co., Ltd.

Preparation of electrolytes and electrodes

All the aqueous electrolytes were prepared by dissolving Zn salts into DI water. The hybrid electrolytes were prepared by dissolving various mass ratios (to the weight of water) of organic solvent in the aqueous electrolytes. Vanadium oxide nanobelt was synthesized based on previous report1. Typically, 3 g V2O5 ( > 98%, Sigma-Aldrich) powder was added into 100 mL 2 mol/L NaCl solution and stirred for 48 h at room temperature. The brown powder was filtered, washed and dried overnight at 60 °C. For ZnI2 cathodes preparation with medium mass loading, first, a mixture of activated carbon, super P conductive carbon and PVDF binder (weight ratio of 8:1:1) was coated on a carbon paper with a mass loading of ~12 mg cm−2. Then, an ethanol solution of ZnI2 (0.1 m) was drop coated on the activated carbon electrode and vacuum dried at 120 °C. For Zn | |ZnI2 coin cells, Cu foils pre-electrodeposited with Zn was used as anode, and 100 uL of electrolyte was added. For ZnI2 cathodes with high mass loading, a mixture of active material ZnI2, activated carbon, CNTs and PTFE binder was pressed to form freestanding electrodes in a weight ratio of 6:6:1:1. The operating voltage range for ZnI2-based full batteries was 0.6–1.6 V under room temperature conditions (~20 °C). Pouch cell shared the same method. For V2O5 cathodes preparation, a mixture of active material V2O5, CNTs and PTFE binder was pressed on a carbon paper in a weight ratio of 7:2:1.

Cell assembling

The Zn | |Cu cells were assembled with Zn foil, Cu foil, and glass fiber separator. Zn symmetric cells were fabricated using Zn foils and a glass fiber separator. If not specified otherwise, Zn foils with a thickness of 100 µm and a diameter of 12 mm were used in half-cell configurations. The electrolyte volume for each cell was maintained at 60 μL. The NP ratio of full cell is defined as the ratio of the capacity of the anode to the cathode. For the anode-free Cu | |ZnI2 cell, the NP ratio is 0. For the Zn | |ZnI2 cells in Fig. 6a, we first plated 1 mAh of zinc on the surface of the copper foil as the anode, so the NP ratio is close to 0.8. For Zn | |V2O5 cells, a 10 µm zinc foil is employed as anode with theoretical areal capacity of about 5.8 mAh cm−2, so the NP ratio is about to 1.7. The half cells were tested in an open environment, while the full cells were evaluated in an environmental chamber. Unless otherwise specified, both half cells and full cells were tested at approximately 20 °C using regular coin cells (2032 type, made of 316 stainless steel). The diameters of the cathodes, anodes, and separators were 12 mm, 12 mm, and 14 mm, respectively. The specific energy density was calculated using the formula: \({E}_{{sp}}=\frac{{E}_{{dcg}}}{{m}_{a}+{m}_{c}}\). Here Edcg is the discharge energy, and ma and mc are the total mass of the anode and cathode electrodes, respectively.

Electrochemical tests

Ionic conductivity measurements were performed using custom-made cells equipped with two platinum electrodes and calibrated using a standard conductivity solution (Mettler Toledo, 12.88 mS cm−1). Aurbach’s protocol was utilized for CE tests at high depth of discharge (DOD). Initially, a fixed capacity of Zn (5 mAh cm−2, denoted as Qr) was plated and stripped on a Cu foil at a current density of 5 mA cm−2. Subsequently, the same capacity of Zn was plated as a reservoir, followed by plating and stripping at a fixed capacity (Qc). After several cycles (n), the final stripping process was carried out to a predefined upper cut-off potential (0.5 V), corresponding to the remaining Zn (Qs). The average CE was calculated using: \({CE}=\frac{n*{Q}_{c}+{Q}_{s}}{n*{Q}_{c}+{Q}_{r}}\). The electrochemical performance of the cells was evaluated using the Neware (Shenzhen) battery testing system (CT4008) at room temperature. Chronoamperometry (CA) and Tafel tests were conducted on a CHI760E electrochemical workstation. EIS tests were conducted on a Gamry 1010E electrochemical workstation.

Characterizations

Raman spectroscopy was conducted using a Horiba HR Evolution microscope with an excitation wavelength of 532 nm. FT − IR spectroscopy was performed with an attenuated total reflection (ATR) accessory on an instrument from Thermo Fisher Scientific, Frontier. X-ray adsorption fine structure (XAFS) analysis was carried out in transmission mode at the Singapore Synchrotron Light Source, specifically at the XAFCA beamline, and at the 8-BM beamline of the National Synchrotron Light Source II. XPS analysis was performed using a Thermo Fisher Scientific Model K−Alpha spectrometer equipped with Al Kα radiation (1486.6 eV). All XPS spectra were calibrated with the binding energy in the C 1 s spectra set at 284.8 eV. High-resolution transmission electron microscopy (HRTEM) was conducted with a Talos F200X microscope (accelerating voltage 200.0 kV) to examine the solid-electrolyte interphase (SEI). XRD patterns were obtained using an advanced X-ray diffraction with Cu−Kα radiation (λ = 1.5418 nm). SEM images were acquired from a field emission scanning electron microscope (JSM − 7100 F). Ion distribution analysis was performed with time-of-flight secondary ion mass spectrometry (IONTOF GmbH, Germany measurements).