• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Aerophilic Triphase Interface Tuned byCarbon Dots Driving Durable andFlexible Rechargeable Zn?Air Batteries

    2023-02-20 01:09:54KuixingDingYuYeJiugangHuLimingZhaoWeiJinJiaLuoShanCaiBaichengWengGuoqiangZouHongshuaiHouXiaoboJi
    Nano-Micro Letters 2023年2期

    Kuixing Ding, Yu Ye, Jiugang Hu?, Liming Zhao, Wei Jin, Jia Luo, Shan Cai,Baicheng Weng, Guoqiang Zou, Hongshuai Hou?, Xiaobo Ji

    ABSTRACT Efficient bifunctional catalysts for oxygen reduction reaction (ORR) and oxygen evolution reaction (OER) are vital for rechargeable Zn-air batteries (ZABs). Herein, an oxygen-respirable sponge-like Co@C–O–Cs catalyst with oxygen-rich active sites was designed and constructed for both ORR and OER by a facile carbon dot-assisted strategy. The aerophilic triphase interface of Co@C–O–Cs cathode efficiently boosts oxygen diffusion and transfer. The theoretical calculations and experimental studies revealed that the Co–C–COC active sites can redistribute the local charge density and lower the reaction energy barrier. The Co@C–O–Cs catalyst displays superior bifunctional catalytic activities with a half-wave potential of 0.82 V for ORR and an ultralow overpotential of 294 mV at 10 mA cm?2 for OER. Moreover, it can drive the liquid ZABs with high peak power density(106.4 mW cm?2), specific capacity (720.7 mAh g?1), outstanding long-term cycle stability (over 750 cycles at 10 mA cm?2), and exhibits excellent feasibility in flexible all-solid-state ZABs. These findings provide new insights into the rational design of efficient bifunctional oxygen catalysts in rechargeable metal-air batteries.

    KEYWORDS Aerophilic triphase interface; Oxygen-rich active sites; O2 diffusion; Bifunctional oxygen catalyst; Flexible rechargeable Zn-air battery

    1 Introduction

    Nowadays, flexible electronics with portable and wearable characteristics have gradually attracted extensive attention in practical applications, such as foldable smartphones, roll-up displays, and smart bracelets [1–3]. To realize these novel energy storage devices, power sources with high energy density and excellent flexibility are essential [4, 5]. Flexible rechargeable zinc-air batteries (ZABs) are considered to be one of the most promising flexible/wearable energy storage devices owing to their high theoretical energy density (1086 Wh kg?1), low cost, remarkable safety, and environmental friendliness [6, 7]. Oxygen evolution reaction (OER) and oxygen reduction reaction (ORR) are essential for the charging and discharging processes, which determine the overall energy efficiency of ZABs. Since both OER and ORR undergo four electron transfer processes, their extremely slow proton coupling and electron transfer kinetics limit the energy conversion efficiency [8, 9]. So far, although precious Pt/RuO2are outstanding catalysts for ORR/OER [10, 11],their high price, single catalytic function, and poor stability hinder their large-scale application [12]. Therefore, exploring cost-effective and non-noble-metal bifunctional oxygen catalysts for ORR and OER will greatly facilitate the development of flexible rechargeable ZABs.

    In the past several years, various transition metal (Mn, Fe,Ni, and Co) catalysts, including transition metals [13, 14],oxides [15, 16], nitrides [17], phosphides [18], and chalcogenides [19, 20], have become extremely promising ORR and OER catalysts due to their abundant reserves and adjustable electronic states [21]. However, many transition metal catalysts have poor electrical conductivity and chemical durability, which seriously hinder their catalytic performance. It has been widely reported that the coupling between carbon materials and transition metal catalysts exhibits admirable activity and durability [22, 23]. In particular, the inevitable agglomeration of transition metal catalysts during hightemperature sintering can be effectively restrained through the confinement of carbon supports. Moreover, heteroatom (N, S, P, B, F, and Se) doping in carbon skeleton can activate the spin states and redistribute the charge of carbon atoms by adjusting the Fermi level ofsp2carbon[24], thereby improving the activities for ORR and OER.However, the influence of heteroatomic oxygen-containing species in the carbon-based catalysts on their electronic states and oxygen-catalytic activities is often ignored. Since the electronegativity of O (3.44) is greater than that of C atom (2.55), the doping of oxygen may regulate the electron delocalization of carbon skeleton and the adsorption of ORR/OER intermediates [25]. Moreover, oxygen doping can introduce abundant defection into the carbon skeleton as catalytic sites [26]. It is feasible to prepare bifunctional ORR/OER catalysts with high activity and excellent stability by coupling active metals and oxygen-doped carbon [27].

    The performance of bifunctional catalysts not only relies on the intrinsic active sites but also highly depends on the utilization efficiency of the catalysts at the solid–liquid–gas three-phase interface. Several strategies have been reported to improve the three-phase interface environment, including surface/interface engineering [28], anion/cation modification [29], and all-in-one design [30], and so on. Wen et al.fabricated Co and Fe atoms encapsulated in the N,S-codoped hollow carbon spheres, and found that the CoFe-SNC catalyst exhibits excellent bifunctional ORR/OER activity and impressive Zn-air battery performance due to the continuous oxygen mass transfer at the three-phase contact region [31].Tian et al. synthesized bifunctional B/N codoped carbon catalyst on carbon fibers for ORR/OER and found that the unique porous structure of BN-C/CF catalyst not only facilitates the exposure of more active sites, but also provides a specific channel for oxygen transport, making it exhibit remarkable capability for high-current and long-cycle Znair battery [32]. Zeng et al. designed NiCo2O4@N-OCNT composite catalysts by in situ encapsulation of NiCo2O4nanoparticles on nitrogen–oxygen double-doped carbon nanotubes. The unique tubular structure of NiCo2O4@N-OCNT is beneficial to repel the formation of bubbles and alleviate the pinning effect of the bubbles at the three-phase interface[33]. These studies demonstrate that optimizing the surface/interface configuration of bifunctional catalysts is vital to sufficient exposure of active sites and rapid oxygen transport of air–cathode for flexible rechargeable ZABs.

    Herein, guided by the density functional theory (DFT)calculations, an oxygen-respirable Co@C–O–Cs porous sponge with aerophilic triphase interface and oxygen-rich active sites was constructed with the assistance of carbon dots (CDs). Benefited from the high hydrophilicity and aerophilicity of Co@C–O–Cs at the three-phase interface, the O2diffusion and mass transfer were accelerated, thus boosting the ORR/OER kinetics in alkaline solution. Theoretical calculations and electrochemical experiments revealed that the oxygen-rich active sites modulate the local charge density and lower the reaction energy barrier of carbon sponge,thus enhancing the ORR/OER activity. As expected, the Co@C–O–Cs catalyst delivers a decent half-wave potential of 0.82 V for ORR, and a low overpotential of 294 mV at 10 mA cm?2for OER under an alkaline medium. The Co@C–O–Cs-based rechargeable liquid ZAB delivers superior performance with a high peak power density (106.4 mW cm?2), maximum specific capacity (720.7 mAh g?1), and outstanding long-term cycle stability (over 750 cycles) to Pt/C + RuO2-based ZAB. Furthermore, the flexible all-solidstate Co@C–O–Cs-based ZAB exhibits high power density and good flexibility at 0–180° bending.

    2 Experimental Section

    2.1 Materials

    All the reagents were analytical pure and used without any further treatment. Acetaldehyde (40% water solution) was supplied by Shanghai McLean Biochemical Technology Co.,Ltd. Sodium hydroxide was purchased from Tianjin Hengxing Chemical Reagent Co., Ltd. Hydrochloric acid and anhydrous ethanol were purchased from Chengdu Kelong Chemical Co., Ltd. Cobalt nitrate hexahydrate were provided by Sinopharm Chemical Reagent Co., Ltd. Nafion solution(5 wt% in alcohol and water) and ruthenium dioxide (RuO2,99.9%) were provided from Sigma-Aldrich. A commercial Pt/C catalyst (20 wt% Pt on the black carbon) was purchased from Alfa Aesar. High-purity water (18.2 MΩ cm) was used for all the experimental processes.

    2.2 Material Synthesis

    Carbon dots were synthesized by reacting sodium hydroxide with acetaldehyde at room temperature according to our previous work [34]. To obtain 3D porous sponge-like Co@C–O–C (denoted as Co@C–O–Cs), 1.5 g of cobalt nitrate and 2.0 g of CDs were dispersed in 50-mL absolute ethanol and magnetically stirred for 10 min to form a uniform dark brown solution. The solution was then placed in a round-bottomed flask and left to evaporate solvent completely, resulting in a dark brown powder. After loading into a ceramic boat, the dark brown powder was heated at 1000 °C under an Ar atmosphere for 2 h at a rate of 5 °C min?1, and then naturally cooled to 25 °C to obtain Co@C–O–Cs. The honeycomb-like Co@C–O–C (denoted as Co@C–O–Ch), vesicles-like Co@C–O–C (denoted as Co@C–O–Cv), and flake-like Co@C–O–C (denoted as Co@C–O–Cf) were also prepared by the same method except for the different calcination temperature (600, 800,1200 °C). The Co@C–O–Cs was washed with 5 mol L?1hydrochloric acid and distilled water to prepare the C–O–Cs control sample.

    2.3 Material Characterization

    The microstructure of samples was characterized by using scanning electron microscopy (SEM, Magellan 400) and transmission electron microscopy (TEM, JEOL JEM 2100F).Power X-ray diffraction (XRD) were recorded on a Rigaku Ultimate Iv diffractometer with Cu Kα radiation (λ= 1.54 ?)operating at 40 kV and 250 mA. Raman spectroscopy was performed using a Raman spectrometer (Renishaw, inVia Reflex) with 532 nm-wavelength incident laser light. The Brunauer–Emmett–Teller (BET) specific surface area and pore size distribution were measured by a specific surface area and pore size analyzer (BET ASAP2460) with N2as absorbate. The oxygen adsorption capacity of the prepared materials was determined from its O2adsorption isotherms,which were obtained using a gas analyzer (Belsorp-Mini II,Japan). Inert atmosphere TG experiments were conducted on Netzsch-STA449F5 under the N2atmosphere from 30 to 1200 °C with heating rate of 10 °C min?1. Oxidizing atmosphere TG experiments were conducted on Netzsch-STA449F5 under the oxygen atmosphere from 30 to 800 °C with heating rate of 10 °C min?1. X-ray photoelectron spectroscopy (XPS) was measured by Thermo Scientific EscaLab 250Xi. All XPS profiles were corrected using the C 1sline at 284.6 eV as an internal standard. C K-edge, Co L-edge,and O K-edge were measured at the beamline U19 of national synchrotron radiation laboratory (NSRL, Hefei) in the total electron yield mode by collecting the sample drain current under a vacuum better than 10?7Pa. The Co K-edge were collected at BL11B station in Shanghai Synchrotron Radiation Facility and 4B9A beamline at the Beijing Synchrotron Radiation Facility. The acquired EXAFS data were extracted and processed according to standard procedures using the ATHENA module implemented in the IFEFFIT software packages. Thek3-weighted EXAFS spectra were obtained by subtracting the post-edge background from the overall absorption and subsequently normalizing with respect to the edge-jump step.

    2.4 Electrochemical Measurements

    All of the ORR and OER performance tests were performed at 25 ± 0.5 °C using a typical three-electrode system on a CHI 760E electrochemical workstation (Chenhua, Shanghai). A rotating disk electrode (RDE) with a glassy carbon disk (5.0 mm diameter) was typically used as the substrate for measuring ORR and OER. Ag/AgCl (saturated with 3.5 mol L?1KCl) and a graphite rod were used as the reference and counter electrodes, respectively. And a 1.0 mol L?1KOH solution was employed as the electrolyte. All potentials were converted to the reversible hydrogen electrode (RHE)scale according to the Nernst equation:E(RHE) =E(Ag/AgCl) + 0.198 V + 0.059 × pH. The RDE or rotating ring disk electrode (RRDE) coated with catalyst inks was used as the working electrode to investigate the kinetics behaviors of Co@C–O–Cs. To prepare the catalyst inks, 5.0 mg catalyst powders were dispersed into a mixed aqueous solution containing ethanol (490 μL), H2O (490 μL), and Nafion solution (5 wt%, 20 μL) under sonication for 1 h. Then, the obtained homogeneous catalyst inks were dropped onto RDE or RRDE and dried at 25 ± 0.5 °C to afford a mass loading capacity of about 0.26 mg cm?2. For comparison,the catalyst inks of both benchmark Pt/C and RuO2with carbon black were also prepared with a loading amount of 0.26 mg cm?2. Prior to the ORR and OER catalytic activity experiments, the electrolyte was purged with high-purity O2for 30 min to obtain an O2-saturated solution.

    The ORR performance of the as-prepared samples was first investigated by CV in a N2/O2-saturated 1.0 mol L?1KOH solution at 25 ± 0.5 °C at a sweep rate of 50 mV s?1.LSV experiments were performed in a O2-saturated 1.0 mol L?1KOH solution at different rotating speeds(400–1600 rpm) with a scan rate of 5 mV s?1. For RDE measurements, the K–L Eqs. (1–3) were used to calculate the kinetic current density (Jk) and transferred electron number (n):

    whereIdandIrrepresent the disk and ring currents, respectively; andNis the current collection efficiency of the Pt ring, which was determined to be 0.37 [35]. The stability of the samples was examined by collecting current–time (it) chronoamperometric responses at 0.67 V (vs. RHE) in 1.0 mol L?1KOH for 12,000 s.

    For OER measurements, the LSV scan was conducted from 1.2 to 1.8 V (vs. RHE) at a scan rate of 5 mV s?1and corrected by 95%iRcompensation to compare the performance of the catalysts, whereiis the measured current andRis the compensated resistance between the working and reference electrodes. The overpotential was calculated usingη=E(vs. RHE) ? 1.23 V. The OER stability of the samples was examined based oni-tchronoamperometric responses at 1.8 V (vs.RHE) for 12,000 s. The ECSAs of the catalysts were compared using the double-layer capacitance(Cdl) in the non-Faradaic region from 1.20 to 1.30 Vvs.RHE. The scanning rates were set at 10, 20, 40, 60, 80, and 100 mV s?1. EIS experiments were performed using a Zahner electrochemical analyzer (Zennium, Germany) in the frequency range of 100 kHz to 10 mHz. The applied potential and excitation amplitude were set to 1.6 V (vs.RHE) and 5 mV, respectively.

    2.5 Assembly andMeasurements ofLiquid Zn?air Battery

    To evaluate the performance of liquid Zn-air batteries, a home-made two-electrode device was constructed. The air electrode was made of nickel foam with a gas diffusion layer on the air-facing side and a catalyst layer on the electrolytefacing side. The gas diffusion layer with an effective area of 1.0 cm2allowed O2diffusion from ambient air to the catalyst sites. The catalyst layer was made by drop-casting the catalyst ink onto the nickel foam with a loading capacity of 6.0 mg cm?2for all the catalysts. The anode was a polished Zn plate with a thickness of 1.0 mm (99.99% purity). The electrolyte solution contained 6.0 mol L?1KOH and 0.2 mol L?1Zn(Ac)2. A rechargeable battery fabricated using a mixture of commercial Pt/C (20%) and RuO2(99.95%) with a mass ratio of 1:1 was used for comparison. All battery performance tests were performed under an ambient atmosphere. The charge–discharge cycling tests of the rechargeable liquid ZABs were performed at 10 mA cm?2with 20 min per cycle (10 min for charging and 10 min for discharging)using the Land-CT2001A testing system.

    2.6 Assembly andMeasurements ofFlexible Solid?state Zn?air Battery

    The gel polymer electrolyte (GPE) was firstly prepared as follows: 2.0 g of PVA powder was added into 20 mL of deionized water, followed by magnetically stirring for 1 h at 90 °C. After the solution was transformed into transparent state, 2 mL of 18.0 mol L?1KOH solution containing 0.2 mol L?1Zn(Ac)2was added into the mixture. After stirred for 0.5 h, the gel was poured into a container and frozen at?20 °C for 2 h. The PVA gel electrolyte was obtained after thawing at room temperature. In a typical assembly of a flexible all-solid-state Zn-air battery, a clean Zn foil (0.08 mm thickness) was used as anode. The air electrode was made by dropping catalyst ink onto a flexible nickel foam substrate with a catalyst loading capacity of 1.5 mg cm?2. Then, the flexible solid-state Zn-air battery was assembled with air electrode and Zn foil placed on the two sides of the PVA gel electrolyte, and two pieces of acrylic tape were used to seal the device. The galvanostatic charge–discharge curves were recorded by chronopotentiometry at a current density of 1 mA cm?2with 20 min per cycle (10 min charge and 10 min discharge).

    2.7 DFT Calculations

    Computational calculations were performed based on DFT with generalized gradient approximations (GGA) [36, 37]implemented in the Vienna ab initio simulation package(VASP) [38, 39]. Thekpoints in the Brillouin zone were sampled with 3 × 2 × 1 unit cell for the structural optimization, with four lines being used along the Brillouin zone for the self-consistent calculations. The vacuum layer was set to 15 ? filled along thecaxis of the unit cell to avoid interactions between the periodic images. The pseudopotential generated in a projector augmented wave (PAW) [40] was used for the electron–ion interactions. In addition, the influence of spin polarization was considered in the calculation process. During geometric optimization, the bottom atomic layer was fixed, and all other atoms and adsorbates were fully relaxed until the maximum force on each atom was less than 0.02 eV ??1. The total energy convergence criterion between every two electronic steps was set to 10?5eV. All models were geometrically optimized prior to the singlepoint energy calculations.

    Changes in the free energies of the ORR and OER were obtained using the DFT calculations. The ORR could proceed incompletely through a two-electron pathway and reduce O2to H2O2, or completely progress via a four-electron route and reduce O2to H2O; the latter was adopted herein because of its efficiency and desirability. Therefore,the 4e?reaction mechanism about ORR could be presented as below:

    where * refers to active sites; (l) and (g) stand for liquid and gas phases, respectively; and *O, *OH, and *OOH are the adsorbed intermediates. As a reverse reaction of ORR, the mechanism of OER is described as below:

    where ?Eis the reaction energy; ?ZPE is the change of zero-point energy;T(298.15 K) is temperature; ?Sis the difference in entropy. The zero-point energies were calculated from the vibration frequencies [41]. The entropies were taken from standard tables for gas-phase molecules [42].Gas phase H2O at 0.035 bar and 300 K was used as the reference state because it equilibrated with liquid water at this situation.Note: the free energy of oxygen is not a correction in DFT calculation, so the free energy of oxygen is usually calculated from H2O and H2. The above adsorption energy values of ΔG*O, ΔG*OH, and ΔG*OOHhave been revised.

    3 Results andDiscussion

    3.1 Theoretical Design ofCo@C?O?Cs Catalyst

    The effect of different oxygen-containing active sites in the Co@C–O–Cs catalysts on their bifunctional ORR/OER performance was evaluated by DFT calculations. The Co–C–C = O, Co–C–COOH, and Co–C–COC models were constructed for theoretical analysis (Figs. 1a-c and S1).As shown in Fig. S2, the calculated band structures show that the gap between the conduction and valence band of the three theoretical models are negative, indicating that they have no obvious band gap and have strong electrical conductivity and metallic properties, which will facilitate electron transport. The electronic projected density of states(PDOS) of Co–C–C = O, Co–C–COOH, and Co–C–COC is shown in Fig. 1d–f. Compared with the Co–C–C = O and Co–C–COOH models, the Co–C–COC shows a slight increase in electron-occupied states near the Fermi level,indicating that it has a higher electrical conductivity and carrier density, which contribute to the enhancement of electron transfer capacity between the catalyst surface and the adsorbed intermediates. Combined with thed-band theory,the adsorption strength and stability of active sites are highly dependent on the electron filling degree of antibonding orbitals between the atomic orbitals of oxygen-containing intermediates and thed-electron state of metal [43, 44]. The obvious increase in thed-band center in Co–C–COC indicates that the presence of –C–O–C species may enhance the adsorption capacity of Co–C–COC sites for oxygen intermediates during ORR and OER processes. The charge density difference analysis of three oxygen-containing functional groups was further performed (Figs. 1g–i and S3).Compared with the Co–C–C = O and Co–C–COOH models, the Co–C–COC sites has a higher local charge density,indicating that the electrons are mainly concentrated on carbon atoms. These accumulated electrons can be transferred to oxygen intermediates when O2is adsorbed on the Co–C–COC active sites [45, 46]. The Co–C–COC model has a large overlap between O-2pstates and Co-3dstates (Fig.S4), further confirming the presence of abundant –C–O–C active sites can improve the oxygen adsorption capacity of Co@C–O–Cs catalyst, which is consistent with the above differential charge results. Moreover, the O2adsorption energy of the Co–C–COC active sites is ?1.073 eV,which is lower than that of the Co–C–C = O active sites of(?0.761 eV) and Co–C–COOH active sites (?0.839 eV)(Fig. 1j).

    Fig. 1 Side view of three theoretical models: a Co–C–C = O, b Co–C–COOH, and c Co–C–COC; Projected total state density of d Co–C–C = O, e Co–C–COOH, and f Co–C–COC; Side view of three theoretical models with a charge density difference of g Co–C–C = O, h Co–C–COOH, and i Co–C–COC, where the yellow and cyan represents electron depletion and accumulation, respectively; j O2 absorption energy values on Co–C–C = O, Co–C–COOH, and Co–C–COC models; k Gibbs free energy plots of ORR at U = 1.23 V, and l Gibbs free energy plots of OER at U = 1.23 V on Co–C–C = O, Co–C–COOH, and Co–C–COC models

    The ORR/OER pathways on the three kinds of oxygencontaining active sites were further disclosed. Both the ORR and OER processes consist of four adsorbed intermediates:*OO, *OOH, *O, and *OH (* refers to active sites) (Fig.S5) [47]. In all ORR processes, the adsorption free energies of all three models show a decreasing trend atU= 0 V (Fig.S6a), indicating that the reaction can proceed spontaneously.At the potentialU= 1.23 V, the overpotential of Co–C–COC active sites is 0.38 eV (Fig. 1k), which is lower than that of Co–C–C = O active sites (0.55 eV) and Co–C–COOH active sites (0.49 eV). These results indicate that the ORR activity of the Co@C–O–Cs catalyst is dominated by the strong electronic interaction between the metallic Co and oxygenrich carbon skeleton. For OER, the free energy of these four elementary reactions tends to rise at the potential ofU= 0 V(Fig. S6b). As shown in Fig. 1l, the Co–C–COC active sites exhibit an overpotential of 0.39 eV, which is smaller than that of Co–C–C = O active sites (0.52 eV) and Co–C–COOH active sites (0.45 eV). This low overpotential can effectively promote the kinetics of the OER process of Co@C–O–Cs catalyst. The lower ORR and OER reaction energy barriers of Co@C–O–Cs catalyst may be mainly attributed to the fast adsorption/desorption of oxygen-containing intermediates on the Co–C–COC active sites [48]. The proposed ORR/OER pathways on the Co@C–O–Cs catalyst is shown in Fig. S7. Therefore, the oxygen-containing active sites can effectively regulate the adsorption of intermediates on the Co@C–O–Cs catalyst surface, thus promoting the intrinsic activity of ORR/OER.

    3.2 Construction andCharacterization ofCo@C?O?Cs Catalyst

    Inspired by DFT calculations, a CDs-assisted synthesis strategy was proposed to construct a three-dimensional spongelike Co@C–O–Cs catalyst with encapsulated abundant Co NPs and oxygen-rich active sites, as shown in Fig. 2a. Initially, the carbon dots (CDs) with oxygen-rich functional groups were fabricated by condensation reaction between sodium hydroxide and acetaldehyde at 25 °C. Then, CDs and cobalt salt were incorporated into ethanol. Since the used CDs have been proven to be extremely well-dispersed in ethanol, a brownish-yellow solution with a homogeneous dispersion can be obtained after stirring, where partial Co2+may be coordinated with the oxygen-rich functional groups of CDs by the electrostatic interaction. When ethanol in the mixed solution was completely evaporated by rotary evaporation, a uniformly dispersed dark brown aggregate was obtained. Finally, the desired Co@C–O–Cs catalyst with abundant Co NPs and oxygen-containing active sites was obtained by annealing under the Ar atmosphere. During the heat treatment, the oxygen-rich CDs can be decomposed at high temperature to generate a large amount of gas and some carbon atoms evolved into graphitized carbon. Meanwhile,the cobalt ions can be reduced to Co NPs by the carbon.Under the catalysis role of Co NPs and the rapid expansion of gas, the aggregated zero-dimensional CDs can transform into three-dimensional (3D) porous carbon sponge. Co NPs can be highly dispersed within the carbon sponge during the rapid volume expansion process [49].

    Scanning electron microscopy (SEM) images of the asobtained Co@C–O–Ch, Co@C–O–Cv, Co@C–O–Cs, and Co@C–O–Cf samples at 600, 800, 1000, and 1200 °C are shown in Figs. 2b–c and S8. It can be seen that Co@C–O–Cs sample exhibits a 3D carbon sponge structure with large vesicle pores, while Co@C–O–Ch, Co@C–O–Cv, and Co@C–O–Cs samples exhibit honeycomb-like, vesicleslike, and flake-like structure, respectively. With the increase in carbonization temperature, the vesicle pores and fluffy degree of carbon sponge gradually become larger; however,the carbon skeleton is broken, resulting in the decrease in vesicle pores when the temperature reaches 1200 °C. The fluffy 3D oxygen-rich carbon sponge and the interdigitated macropore may facilitate the exposure of active sites and mass transport during the catalysis process [50]. From the transmission electron microscopy images of Figs. 2f and S9(a, d, g), it can be found that metallic Co NPs with an average size of ca. 20 nm are dispersed in the 3D porous carbon sponge, effectively increasing the number of active sites [51]. The dispersed Co NPs could result from the cooperative coordination between Co2+and C–O–C functional groups in CDs, which suppresses the agglomeration of metallic Co NPs during high-temperature carbonization.When the temperature is lower than 800 °C, the bivalent cobalt species cannot be completely reduced to metallic Co. For Co@C–O–C, the high-resolution TEM images in Figs. 2g and S9(b, e, h) show the lattice spacing of 0.20 and 0.34 nm, which corresponds to the (111) plane of Co NPs and the (002) plane of graphite phase, respectively [52]. The selected area electron diffraction pattern of Co@C–O–Cs catalyst (Fig. S10) verifies the (111), (200), and (220) lattice planes of Co NPs, respectively [53, 54]. The lattice of graphite phase indicates that the reduced cobalt NPs may act as a catalyst to boost the formation of fluffy porous carbon sponge and the graphitization of carbon. The element distribution images in Figs. 2h and S9 (c, f, i) show that the Co, O, and C elements are uniformly distributed in the Co@C–O–Cs catalyst.

    Fig. 2 a Schematic fabrication procedure of Co@C–O–Cs catalyst; b–e SEM images of Co@C–O–Ch, Co@C–O–Cv, Co@C–O–Cs, and Co@C–O–Cf; f–h TEM, HRTEM, and mapping image of Co@C–O–Cs catalyst; i XRD patterns, j oxygen adsorption curve, k contact angle of the bubble in KOH solution of Co@C–O–Ch, Co@C–O–Cv, Co@C–O–Cs, and Co@C–O–Cf; l Schematic aerophilicity of Co@C–O–Cs catalyst in KOH solution

    Based on thermogravimetric analysis under Ar atmosphere (Fig. S11), the precursor mass decreases rapidly when the temperature rises. The large mass loss before 500 °C should be caused by the decomposition of CDs, and then the slight mass loss may result from the reduction of Co2+and graphitization of carbon. The X-ray diffraction patterns in Fig. 2i show that the as-obtained Co@C–O–Ch,Co@C–O–Cv, Co@C–O–Cs, and Co@C–O–Cf have identical diffraction peaks at 44.2°, 51.5°, and 75.8°, which are well-indexed to the (111), (200), and (220) planes of Co(PDF#15–0806) [55, 56]. As the temperature increases from 600 to 1200 °C, the diffraction peaks of metallic Co become sharper, indicating an increase in crystallinity. The wide diffraction peak at 25° results from the (002) plane of graphitic carbon, indicating the formation of graphitic structure during the annealing process. According to the TG results of Co@C–O–Ch, Co@C–O–Cv, Co@C–O–Cs, and Co@C–O–Cf in air atmosphere (Fig. S12a), the determined Co contents of four samples are 20.5%, 17.0%, 20.0%, and 30.2%, respectively. X-ray diffraction pattern confirms that the residual oxidation product at 800 °C was Co3O4(Fig.S12b). The Raman peaks located at 1340 and 1580 cm?1belong to the D-band (structural defects) and the G-band(graphite structure), respectively [57, 58]. The integrated intensity ratio of the D to G peak (ID/IG), which is an indicator of the degree of disorder, was calculated from the Raman spectra of all samples. The results verify that theID/IGvalue of the Co@C–O–Cs is 0.84, which is lower than that of Co@C–O–Ch (0.99), Co@C–O–Cv (0.98),and Co@C–O–Cf (0.85), indicates that the graphitization degree of the catalysts increases as the carbonization temperature increases (Fig. S13). As shown in Fig. S14a,all of the isotherms belong to type-IV curve combined with the H4-type hysteresis loop, suggesting the existence of microporous and mesoporous in the Co@C–O–Ch,Co@C–O–Cv, Co@C–O–Cs, and Co@C–O–Cf [59]. The Brunauer–Emmett–Teller (BET) specific surface area of Co@C–O–Cs (203.7 m2g?1) is much larger than Co@C–O–Ch (127.3 m2g?1), Co@C–O–Cv (140.7 m2g?1),and Co@C–O–Cf (155.9 m2g?1), indicating that an increase in temperature can induce a larger specific surface area.When the temperature increases to 1200 °C, the carbon sponge structure collapses and the macropores are broken,resulting in the decrease in specific surface area. Moreover,it can be found from Fig. S14b that these samples possess a similar non-uniformly distributed nano-meso-micro-hierarchical porous structure. The hierarchical porous structure of Co@C–O–Cs may facilitate the mass diffusion of oxygen species at the three-phase interface during the catalytic process [60]. The ORR and OER involve a gas–liquid–solid tri-phase interface, the mass transport between the reactant and active sites is a vital factor to determine reaction rate.Therefore, oxygen adsorption isotherms were measured to evaluate the oxygen adsorption property of Co@C-O–Ch,Co@C–O–Cv, Co@C–O–Cs, and Co@C–O–Cf. As shown in Fig. 2j, with the increase in carbonization temperature, the oxygen adsorption capacity becomes larger. However, when the temperature reaches 1200 °C, the adsorption of oxygen decreases due to the collapse of the carbon sponge structure for Co@C–O–Cf. In Fig. S15, the maximum specific surface area of Co@C–O–Cs catalyst for oxygen adsorption can reach up to 2.3 m3g?1, which is larger than Co@C–O–Ch(1.1 m3g?1), Co@C–O–Cv (1.4 m3g?1), and Co@C–O–Cf(1.9 m3g?1), indicating the higher adsorption capacity of oxygen for Co@C–O–Cs catalyst [61, 62].

    The wettability in air and the adhesion behavior of oxygen bubbles in KOH solution on the as-fabricated catalysts were investigated by contact angle measurement.The contact angle value of KOH solution on the surface of Co@C–O–Cs catalyst in air is 18.2°, which is lower than that of Co@C–O–Ch (75.1°), Co@C–O–Cv (44.8°),and Co@C–O–Cf (34.6°) (Fig. S16). The good wettability and hydrophilicity of the oxygen-rich carbon sponge facilitate the contact between the active sites and electrolyte, thereby potentially improving the ORR/OER catalytic kinetics on the Co@C–O–Cs catalyst [63]. Meanwhile, as shown in Fig. 2k, the contact angle of oxygen bubbles in KOH solution on Co@C–O–Cs catalyst is 45.9°, which is smaller than that of Co@C–O–Ch (90.1°), Co@C–O–Cv(65.9°), and Co@C–O–Cf (60.5°), which strongly demonstrates that the Co@C–O–Cs catalyst has strong adhesion ability of oxygen bubbles, thus showing excellent aerophilicity of the Co@C–O–Cs catalyst surface. Benefiting from the lower surface roughness and loose sponge structure, the Co@C–O–Cs exhibits the high hydrophilicity in air and aerophilicity in KOH solution (Figs. 2l and S17). Based on the Cassie-Baxter wetting regime [64], the loose sponge structure of Co@C–O–Cs improves the adhesion to oxygen bubbles. This strong surface adhesion role can continuously accelerate the diffusion and transfer of O2, thus the large amount of reactants at the three-phase interface can boost the ORR/OER kinetics in the alkaline electrolyte [65].

    The near-surface chemical structure and components of the Co@C–O–C at different temperatures were characterized by X-ray photoelectron spectroscopy (XPS). The XPS full spectrum shows the presence of Co, C, and O elements in all samples (Fig. S18a–d). The high-resolution C 1sfor Co@C–O–C at different temperatures (Fig. S18e–h) can be divided into four peaks, including C–C (284.7 V), C–O(285.9 V), C = O (286.7 V), and O–C = O (287.7 V) bands[66], demonstrating the presence of oxygen species on the carbon skeleton. As shown in Figs. 3a and S18i-k, two pairs of 2p3/2/2p1/2doublets in the high-resolution Co 2pspectrum were deconvoluted into three major peaks, including Co0(779.5 and 795.7 eV), Co2+(781.8 and 797.4 eV), and shakeup satellite (786.6 and 804.9 eV) [67]. Notably, the peak corresponding to Co0is smaller than that of Co2+because small Co NPs are prone to surface oxidation, resulting in the presence of Co–O species on the catalyst surface.This phenomenon can be widely observed in related transition metal catalysts. The high-resolution O 1sspectrum(Figs. 3b and S18l-n) can be fitted with C–O–C (533.5 eV),C(O) –OH (535.1 eV), and C = O (532.1 eV) bands [68], further suggesting the formation of abundant oxygen-bearing components in the 3D porous carbon sponge.

    Fig. 3 a Co 2p and b O 1s high-resolution XPS spectra of Co@C–O–Cs; c Co L-edge and d O K-edge XANES spectra of Co@C–O–Cs and reference samples; e Co K-edge XANES spectra of Co@C–O–Cs and reference samples; f Fourier transform spectra of Co@C–O–Cs and reference samples; g WT-EXAFS spectra of Co@C–O–Cs and reference samples

    The soft/hard X-ray absorption spectroscopy analysis was performed to identify the local electronic structure and coordination environment of Co, C, and O atoms in the Co@C–O–Cs catalyst. The C K-edge near-edge X-ray absorption fine structure (NEXAFS) spectra of Co@C–O–Cs (Fig. S19) displays two typical spectroscopic peaks at 285.8 eV (C = C π*) and 292.3 eV (C–C σ*), implying a typicalsp2hybrid carbon structure [69]. Moreover,the peak at 288.9 eV is a fingerprint of thesp3carbon, evidencing the existence of C–O–C in the carbon lattice [70].The Co L-edge NEXAFS spectra of Co@C–O–Cs (Fig. 3c)show two prominent peaks, including Co L3at 781.1 eV and Co L2at 796.4 eV, indicating the 2pcore–hole spin–orbit interaction between Co and carbon skeleton [71]. The O K-edge NEXAFS spectra of Co@C–O–Cs in Fig. 3d exhibit three characteristic peaks, including π*(COO?) (530.2 eV),π*(C = O) (531.6 eV), and σ*(C–O–C) (541.7 eV), respectively [72, 73]. The Co K-edge X-ray absorption near edge structure spectrum of Co@C–O–Cs shows a pre-edge feature at 7708 eV, which is attributed to the 1s→ 3delectron transition (Fig. 3e). The peak intensity of Co@C–O–Cs (1s→ 4pz)is between that of Co foil and CoO reference samples, confirming the existence of interfacial interaction between the Co NPs and porous carbon sponge. Moreover, the near-edge absorption threshold of Co K-edge lies between those of Co foil (Co0) and CoO (Co2+), demonstrating the average valence state of Co is between 0 and + 2 mixed oxidation states [74]. The Fourier transform EXAFS curve of Co@C–O–Cs in Fig. 3f presents a peak at 2.1 ?, which is attributed to the Co–Co scattering paths [75]. Because of the high resolution of wavelet-transforming (WT)-EXAFS analysis, it has become an effective method to discriminate the coordination environment. Compared to the Co foil and CoO samples, the WT-EXAFS plot of Co@C–O–Cs shows two sole contour peak intensity maximums at around 7.1 ??1(Fig. 3g), which corresponds to the Co–Co coordination [76]. This phenomenon indicates the presence of Co–C–O–C coordination configuration in the Co@C–O–Cs catalyst.

    3.3 Electrochemical Evaluation

    The catalytic activities of Co@C–O–Cs catalysts and benchmarked commercial Pt/C toward ORR were evaluated by a typical three-electrode system in an O2-saturated 1.0 mol L?1KOH solution. Cyclic voltammetry profile of Co@C–O–Cs displays an obvious cathodic ORR peak at 0.77 V in the O2-saturated electrolyte, compared with that in the N2-saturated electrolyte (Fig. S20), implying effective catalytic activity for oxygen reduction. Moreover, linear sweep voltammetry (LSV) was conducted to disclose their kinetics and activities for ORR. The Co@C–O–Cs catalyst affords a remarkable ORR performance with a half-wave potential (E1/2) of 0.82 V and an onset and half-wave potential (Eonset) of 0.90 V (Figs. 4a and S21), which are comparable to the benchmarked Pt/C(E1/2= 0.86 V andEonset= 0.93 V), and obviously superior than those of C–O–Cs (E1/2= 0.64 V andEonset= 0.72 V),Co@C–O–Ch (E1/2= 0.71 V andEonset= 0.79 V),Co@C–O–Cv (E1/2= 0.73 V andEonset= 0.83 V),Co@C–O–Cf (E1/2= 0.75 V andEonset= 0.85 V), and the most of recently reported catalysts in Table S1. The corresponding Tafel plots of Co@C–O–Cs catalyst under alkaline conditions is 81.1 mV dec?1(Fig. 4b), which is lower than those of C–O–Cs (119.6 mV dec?1), Co@C–O–Ch (94.3 mV dec?1), Co@C–O–Cv (92.4 mV dec?1), Co@C–O–Cf(90.4 mV dec?1), and is comparable to that of the commercial Pt/C catalyst (79.6 mV dec?1), suggesting the faster ORR kinetics on Co@C–O–Cs catalyst. To investigate the ORR mechanism, the electron transfer number (n) was further determined by analyzing LSV curves at different rotating speeds using the Koutecky–Levich (K–L) equations (Figs. 4c and S22). The K–L plots have shown good linearity relationships (inset of Fig. 4c), confirming the firstorder kinetics process of the Co@C–O–Cs catalyst, which is dependent on the concentration of dissolved oxygen in the alkaline medium [77]. Moreover, the electron transfer number of Co@C–O–Cs is calculated as 3.8, which is close to 4.0 of benchmarked Pt/C, suggesting a four electron transfer pathway during the ORR process. The ORR pathway was further investigated by the rotating ring-disk electrode(RRDE) measurements (Figs. S23 and 4d). The average H2O2yield is 5.6% and the calculated n value is 3.82, which is in agreement with the results obtained from the K–L plot.Therefore, these results verify that the Co@C–O–Cs catalyst exhibits an efficient ORR kinetic process with four-electron transfer reaction mechanism. The stability of Co@C–O–Cs and Pt/C catalysts was assessed by chronoamperometric(CA) measurement within an O2-saturated 1.0 mol L?1KOH solution for 12,000 s. The current density of Co@C–O–Cs is maintained at 93.8% (Fig. S24a), which is higher than the 78.9% retention rate of Pt/C, demonstrating the high catalytic stability of the Co@C–O–Cs catalyst. In addition, a slight decay in theE1/2of Co@C–O–Cs (12.5 mV decay inE1/2) after 1000 cycles by comparing with that of Pt/C catalyst (18.6 mV decay inE1/2) (Fig. S24b), indicating excellent ORR durability of the Co@C–O–Cs catalyst because the Co NPs are encapsulated in the porous oxygen-rich carbon sponges.

    Fig. 4 a ORR LSV measurements in O2-saturated 1.0 mol L?1 KOH solution at 1600 rpm; b Tafel plots of ORR process; c LSV curves of Co@C–O–Cs catalyst under different rotation speeds (inset: K-L plots of the Co@C–O–Cs catalyst); d Transferred electron number and H2O2%of Co@C–O–Cs and Pt/C catalysts; e OER LSV measurements in O2-saturated 1.0 mol L?1 KOH solution; f Tafel plots of OER process; g Cdl curves of various catalysts; h In situ Raman spectra of Co@C–O–Cs catalyst in 1.0 mol L?1 KOH solution at different potentials for OER process; i Overall polarization curves of various catalysts

    The OER catalytic performance of Co@C–O–C was further evaluated by LSV in 1.0 mol L?1KOH solution. The Co@C–O–Cs presents a low overpotential value of 294 mV at a current density of 10 mA cm?2(Fig. 4e), which is obviously lower than C–O–Cs (450 mV), Co@C–O–Ch (380 mV),Co@C–O–Cv (340 mV), Co@C–O–Cf (320 mV), commercial RuO2catalyst (330 mV), and those of previously reported catalysts in Table S2. Moreover, the Tafel slope of Co@C–O–Cs is 83.1 mV dec?1, which is much lower than those of C–O–Cs (184.6 mV dec?1),Co@C–O–Ch (112.4 mV dec?1), Co@C–O–Cv (102.3 mV dec?1), Co@C–O–Cf (86.7 mV dec?1), and comparable to that of RuO2(100.3 mV dec?1) (Fig. 4f), indicating faster OER kinetics on Co@C–O–Cs catalyst [78]. From the Nyquist plots and impedance fitting data shown in Fig.S25 and Table S3, the Co@C–O–Cs catalyst exhibits the smallest charge-transfer resistance (Rct= 15.15 Ω) and the highest constant phase elements (CPE = 9.85 × 10–4mF),reflecting that more efficient electron transfer occurs at the Co@C–O–Cs/electrolyte interface during the catalytic reaction. After a 12,000 continuous reaction, the Co@C–O–Cs catalyst retains 93.7% of the initial current, which is much lower than the commercial RuO2catalyst (86.5%) (Fig. S26a). Additionally, a negligible activity loss of Co@C–O–Cs can be observed after 1000 potential cycles (~7.2 mV atEj=10), which is lower than that of commercial RuO2(19.4 mV) (Fig. S26b). These results suggest that the Co@C–O–Cs catalyst exhibits excellent OER electrochemical stability in alkaline solution. The electrochemically active surface area (ECSA)was obtained to demonstrate the origin of the intrinsic activity of the Co@C–O–Cs catalyst based on the electrochemical double layer capacitance (Cdl) at different scan rates (Fig. S27). The calculatedCdlof Co@C–O–Cs is 18.1 mF cm?2(Fig. 4g), which is significantly higher than those of C–O–Cs (1.3 mF cm?2), Co@C–O–Ch (7.6 mF cm?2), Co@C–O–Cv (12.1 mF cm?2), and Co@C–O–Cf(14.6 mF cm?2), resulting in the excellent mass transportation for the Co@C–O–Cs catalyst. To demonstrate the stability of the Co@C–O–Cs catalyst, the morphology,composition, and elemental electronic states of samples after OER stability tests were characterized, as shown in Fig. S28. It can be seen from Fig. S28a-c that the Co NPs are still distributed in the carbon sponge after cycling and no corrosion occurs. Furthermore, the XRD pattern of the Co@C–O–Cs catalyst after cycling shows the same characteristic peaks with before cycling (Fig. S28d), indicating that no phase transition occurs during cycling. The XPS results of Co@C–O–Cs catalyst (Fig. S28e–g) also show that the positions of the Co, C, and O peaks hardly change.These results suggest that the Co@C–O–Cs catalyst still remains its initial structure after OER stability test.

    The possible catalytic mechanism and surface structure changes of Co@C–O–Cs catalyst during OER under controlled overpotentials in 1.0 mol L?1KOH solution were measured by in situ Raman spectroscopy technique (Fig.S29). As shown in Fig. 4h, there is no feature peak of lattice vibration between 400 and 800 cm?1under the opencircuit condition. Moreover, no new peak can be observed when the voltage is lower than 1.5 V, indicating that OER does not proceed. When the potential rises from 1.5 to 1.8 V, a newly formed vibrational peak at 472 cm?1indicates the formation ofβ-Co(OH)2species [79]. Meanwhile,a more intense peak at 611 cm?1may be attributed to the formation of CoOOH intermediate under the harsh oxidative environment [80, 81]. The increased peak intensity of the CoOOH intermediate is larger than that ofβ-Co(OH)2due to the transformation ofβ-Co(OH)2into CoOOH intermediate. The formation of these intermediates over 1.5 V may result in the low OER overpotential of Co@C–O–Cs catalyst. These results demonstrate that the partial Co–CCOC active sites may be converted into CoOOH species under the high potential and strong alkali environments,which have the potential to boost OER and resist the destruction of Co@C–O–Cs catalyst. The potential gap(ΔE) between the OER overpotential at 10 mA cm?2and ORR half-wave potential (ΔE=Ej=10?E1/2) is generally used to evaluate the bifunctional activity and reversible electrochemical properties [82]. As shown in Figs. 4i and S30, the Co@C–O–Cs catalyst has a comparable ΔEvalue of 0.704 V, which is markedly lower than that of C–O–Cs(1.04 V), Co@C–O–Ch (0.90 V), Co@C–O–Cv (0.84 V),Co@C–O–Cf (0.80 V), Pt/C (1.03 V), commercial RuO2(0.98 V), and some bifunctional ORR/OER catalysts previously reported in Table S4. Therefore, the Co@C–O–Cs has great potential as a bifunctional catalyst as air cathode for rechargeable ZABs.

    3.4 Rechargeable ZABs

    Considering the excellent bifunctional activity of Co@C–O–Cs catalyst, a home-made rechargeable liquid ZAB was assembled. As shown in Figs. 5a and S31,the liquid ZAB consists of a polished Zn plate as anode,Co@C–O–Cs catalyst as cathode, 6.0 mol L?1KOH solution with 0.2 mol L?1zinc acetate (Zn(Ac)2) as electrolyte.The open-circuit voltage of Co@C–O–Cs-based ZAB is 1.426 V, which is closed to Pt/C + RuO2-based ZAB (1.40 V)(Fig. 5b). The discharge and charge curves in Fig. 5c show that under different current densities, the Co@C–O–Csbased liquid ZAB has the smaller voltage gaps than the control device based on Pt/C + RuO2. As shown in Fig. 5d, the Co@C–O–Cs-based liquid ZAB achieves a current density of 73.1 mA cm?2at 1.0 V and a peak power density of 106.4 mW cm?2at 164.5 mA cm?2, which are superior to those of Pt/C + RuO2-based ZAB (49.4 mA cm?2at 1.0 V, maximum 88.7 mW cm?2). Additionally, the specific capacity of the Co@C–O–Cs-based liquid ZAB can reach 720.7 mAh g?1,which is comparable to those of the Pt/C + RuO2-based ZAB(690.7 mAh g?1; Fig. 5e). A light-emitting diode (LED) with‘CSU’ (Central South University) screen can be powered by two liquid ZABs in series with the Co@C–O–Cs catalyst as the air cathode (inset of Fig. 5e), demonstrating its potential application in energy conversion devices. As increasing current density, the discharge potential plateau decreases and remains stable. Even at 25 mA cm?2, the Co@C–O–Csbased liquid ZAB still exhibits a high discharging potential of 1.19 V (Fig. 5f), suggesting its excellent high-rate performance. Furthermore, the cycling stabilities for Co@C–O–Cs and Pt/C + RuO2-based rechargeable liquid ZABs was investigated under conditions of charge/discharge at 10 mA cm?2(Fig. 5g). The Co@C–O–Cs-based liquid ZAB exhibits an ultra-stable operation window of over 750 cycles, which is much better than Pt/C + RuO2-based ZAB (300 cycles).Interestingly, the Co@C–O–Cs-based ZAB exhibits a small charge–discharge voltage gap of 0.79 V at 1.14 V discharge voltage and 1.93 V charge voltage, and the high round-trip efficiency is 59.1%. After 360 cycles, the charge/discharge voltage gap only varies to 0.87 V and round-trip efficiency only decreases to 55.6%. Even after 750 cycles, it still maintains good stability. Therefore, the Co@C–O–Cs-based rechargeable liquid ZAB exhibits the better stability performance than that of Pt/C + RuO2and some recently reported batteries based on non-precious metal catalysts in Table S5.

    Fig. 5 a Schematic diagram of rechargeable liquid ZAB with Co@C–O–Cs catalyst; b Open-circuit plots, c charge and discharge polarization curves, d discharge polarization curves and the corresponding power density curves, and e discharge polarization curves at the current density of 10 mA cm?2 of rechargeable liquid ZABs with Co@C–O–Cs and Pt/C + RuO2 catalysts (inset: photograph of a ‘CSU’ pattern formed by LEDs driven by two rechargeable liquid ZABs with Co@C–O–Cs catalyst); f Discharge curves at various current densities; g Discharge–charge cycling performance of the rechargeable liquid ZABs with Co@C–O–Cs and Pt/C + RuO2 catalysts at 10 mA cm?2

    Fig. 6 a Schematic diagram of rechargeable flexible all-solid-state ZAB with Co@C–O–Cs catalyst; b Open-circuit plots and c charge and discharge polarization curves and the corresponding power density curves of rechargeable flexible all-solid-state with Co@C–O–Cs catalyst; d Discharge curves of the battery at various current densities; e Cycling stability of rechargeable flexible all-solid-state ZAB with different bending angles; f Photographs of LED lightened by three flexible ZABs connected in series when bent at various angles and the three-series all-solidstate ZABs as a wearable bracelet to light up the LED

    With the increasing demand for flexible and portable electronic devices, all-solid-state ZABs are emerging as a promising alternative for efficient and sustainable power supply. Flexible all-solid-state ZABs were further fabricated with Co@C–O–Cs supported on nickel foam as air cathode, alkaline polyvinyl alcohol (PVA) gel as electrolyte,Zn foil as metal anode (Fig. 6a). The Co@C–O–Cs-based all-solid-state ZAB exhibits an open-circuit voltage of 1.43 V (Fig. 6b) and a high power density of 59.1 mW cm?2(Fig. 6c), which is comparable to those of recently reported bifunctional catalysts in Table S6. As shown in Fig. 6d, the discharge voltage gradually decreases from 1.28 to 1.03 V with the current density varying from 1 to 12 mA cm?2.Moreover, as shown in Fig. 6e, the flexible Co@C–O–Csbased ZAB can charge and discharge for a long time at a constant current density of 1.0 mA cm?2, even if the battery is bent at different angles (including 90° and 180°). Under the initial state (0°), the flexible all-solid-state battery exhibits an ultra-low voltage polarization of 0.50 V (charge and discharge voltages are 1.77 and 1.27 V, respectively), and the round-trip efficiency is 71.5%. A low-voltage polarization and high round-trip efficiency can be still remained at 90° (0.55 V, 69.4%) and 180° (0.56 V, 69.1%). Moreover,the voltage polarization rate and round-trip efficiency can maintain at 0.58 V and 68.5% after the battery was returned to the initial state, indicating the excellent flexibility of the Co@C–O–Cs-based all-solid-state ZAB. A working prototype that contains three series of flexible all-solid-state ZABs based on the Co@C–O–Cs catalyst was also established. As shown in Fig. 6f, it can successfully light up the LED with ‘CSU’ indicator even bending at different angles.Moreover, the three series of ZABs can be used as wearable bracelets to light up the LED with stable operating conditions. Compared with rechargeable liquid ZABs, this rechargeable all-solid-state ZAB exhibits more attractive features due to their small size, flexibility, high power density, and high safety.

    4 Conclusions

    In summary, the oxygen-respirable Co@C–O–Cs catalyst with oxygen-rich active sites were successfully constructed with CDs-assisted synthesis strategy as the efficient bifunctional ORR/OER catalysts. The hydrophilicity and aerophilicity of Co@C–O–Cs are beneficial to oxygen diffusion and mass transfer. The density functional theory calculations and experimental studies revealed that the Co–O–COC active sites can modulate the local charge density, lower the reaction energy barrier, and enhance the ORR/OER activity. In situ Raman spectroscopy revealed that partial Co–OCOC active site could transform to the Co–OOH intermediate during OER process to maintain high catalysis activity and stability. The Co@C–O–Cs catalyst delivers a halfwave potential of 0.82 V for ORR and a low overpotential of 294 mV at 10 mA cm?2for OER. The Co@C–O–Csbased ZAB displays a high discharge peak power density of 106.4 mW cm?2, high specific capacity of 720.7 mAh g?1, and outstanding cycling stability over 750 cycles,which outperforms the Pt/C + RuO2. Furthermore, the flexible Co@C–O–Cs-based all-solid-state ZAB exhibits high power density and good flexibility under 0–180° bending degrees. These findings provide new insights into the design of efficient bifunctional catalysts for versatile energy storage and conversion devices.

    AcknowledgementsThis work was financially supported by the National Key Research and Development Program of China(No. 2019YFC1907801), National Natural Science Foundation of China (No. 52174286), the Science and Technology Innovation Program of Hunan Province (2021RC3014), and Innovation-Driven Project of Central South University (No. 2020CX007). Moreover,the authors wish to thank facility support of the 4B9A beamline of Beijing Synchrotron Radiation Facility (BSRF), the BL11B station in Shanghai Synchrotron Radiation Facility (SSRF) and beamline U19 station in National Synchrotron Radiation Laboratory (NSRL,Hefei) for XAS measurement.

    FundingOpen access funding provided by Shanghai Jiao Tong University.

    Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing,adaptation, distribution and reproduction in any medium or format,as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.

    Supplementary InformationThe online version contains supplementary material available at https:// doi. org/ 10. 1007/s40820- 022- 00994-3.

    亚洲国产精品成人综合色| 国产v大片淫在线免费观看| 亚洲三级黄色毛片| av福利片在线观看| 亚洲综合色惰| 在线观看美女被高潮喷水网站 | 亚洲在线观看片| 国产精品1区2区在线观看.| 久久人人精品亚洲av| 国产亚洲欧美在线一区二区| 国产成人av教育| 免费人成在线观看视频色| 在线观看免费视频日本深夜| 亚洲欧美日韩无卡精品| 深爱激情五月婷婷| 国产v大片淫在线免费观看| 久久久久免费精品人妻一区二区| 国产男靠女视频免费网站| 色综合站精品国产| xxxwww97欧美| 老司机深夜福利视频在线观看| 国产精品久久久久久亚洲av鲁大| 听说在线观看完整版免费高清| 日韩 亚洲 欧美在线| 夜夜躁狠狠躁天天躁| 久久香蕉精品热| 国产蜜桃级精品一区二区三区| 欧美激情在线99| 国产精品伦人一区二区| 白带黄色成豆腐渣| 日本与韩国留学比较| 免费电影在线观看免费观看| 国内精品一区二区在线观看| av福利片在线观看| av中文乱码字幕在线| 国产蜜桃级精品一区二区三区| 成年女人永久免费观看视频| 色播亚洲综合网| 99热6这里只有精品| 欧美一级a爱片免费观看看| 欧美激情国产日韩精品一区| 91久久精品电影网| 久久亚洲精品不卡| 舔av片在线| 欧美乱妇无乱码| 一卡2卡三卡四卡精品乱码亚洲| 精品久久久久久久久久免费视频| 韩国av一区二区三区四区| 国产视频一区二区在线看| 村上凉子中文字幕在线| 亚洲av.av天堂| 女生性感内裤真人,穿戴方法视频| 一区二区三区四区激情视频 | 在线观看66精品国产| 两性午夜刺激爽爽歪歪视频在线观看| av福利片在线观看| 神马国产精品三级电影在线观看| 男女做爰动态图高潮gif福利片| 国产黄a三级三级三级人| 91久久精品国产一区二区成人| 久99久视频精品免费| aaaaa片日本免费| 欧美在线一区亚洲| 一本一本综合久久| 精品久久久久久成人av| 亚洲国产精品久久男人天堂| 欧美+日韩+精品| 亚洲 欧美 日韩 在线 免费| 97人妻精品一区二区三区麻豆| 国产一区二区在线av高清观看| 国模一区二区三区四区视频| 亚洲人成伊人成综合网2020| 国产又黄又爽又无遮挡在线| 日韩国内少妇激情av| 日韩免费av在线播放| 亚洲片人在线观看| av视频在线观看入口| 亚洲精品456在线播放app | 久久精品国产亚洲av天美| 丁香六月欧美| 看十八女毛片水多多多| 国产乱人视频| 亚洲片人在线观看| 亚洲无线在线观看| 欧美xxxx黑人xx丫x性爽| 免费在线观看亚洲国产| 99热这里只有是精品50| 成人特级黄色片久久久久久久| 久久国产乱子免费精品| 亚洲激情在线av| av视频在线观看入口| 国产人妻一区二区三区在| 欧美极品一区二区三区四区| 欧美最黄视频在线播放免费| 婷婷亚洲欧美| 国产爱豆传媒在线观看| 日日夜夜操网爽| 国产激情偷乱视频一区二区| 国产乱人伦免费视频| 波多野结衣高清无吗| 欧美午夜高清在线| 亚洲性夜色夜夜综合| 桃红色精品国产亚洲av| 久久精品国产99精品国产亚洲性色| 国产成人a区在线观看| 免费看光身美女| 国产高潮美女av| 九九久久精品国产亚洲av麻豆| 欧美性猛交╳xxx乱大交人| 美女cb高潮喷水在线观看| 国产av一区在线观看免费| 搡老熟女国产l中国老女人| 一区二区三区激情视频| av国产免费在线观看| 色哟哟·www| 国产亚洲精品久久久com| 老司机午夜福利在线观看视频| 99热精品在线国产| 国产精品98久久久久久宅男小说| 免费看a级黄色片| 亚洲av美国av| 91久久精品国产一区二区成人| 中文字幕精品亚洲无线码一区| 国产大屁股一区二区在线视频| 搞女人的毛片| 此物有八面人人有两片| 美女高潮喷水抽搐中文字幕| 亚洲欧美日韩卡通动漫| 免费在线观看成人毛片| 天堂网av新在线| 国产乱人伦免费视频| 亚洲美女黄片视频| 九色成人免费人妻av| 最近最新免费中文字幕在线| 久久99热6这里只有精品| 日韩大尺度精品在线看网址| 看十八女毛片水多多多| 亚洲乱码一区二区免费版| 久久精品综合一区二区三区| 十八禁人妻一区二区| 99久国产av精品| 亚洲成av人片在线播放无| 永久网站在线| 欧美黄色片欧美黄色片| 亚洲aⅴ乱码一区二区在线播放| 欧美色欧美亚洲另类二区| 淫妇啪啪啪对白视频| 成人特级黄色片久久久久久久| 搡女人真爽免费视频火全软件 | 欧美成人免费av一区二区三区| 夜夜躁狠狠躁天天躁| 亚洲欧美日韩东京热| 欧美+亚洲+日韩+国产| 亚洲综合色惰| 国产高清有码在线观看视频| 国产成人av教育| 九色国产91popny在线| 免费无遮挡裸体视频| 国产av在哪里看| 美女高潮喷水抽搐中文字幕| 99久久无色码亚洲精品果冻| 美女 人体艺术 gogo| 亚洲人成网站高清观看| 又黄又爽又免费观看的视频| 国产视频一区二区在线看| 国产欧美日韩精品一区二区| 少妇高潮的动态图| 日韩精品青青久久久久久| 亚洲,欧美精品.| 一区二区三区免费毛片| 老司机午夜十八禁免费视频| 看免费av毛片| 亚洲人成网站在线播放欧美日韩| av国产免费在线观看| 18禁在线播放成人免费| 午夜精品久久久久久毛片777| 国产 一区 欧美 日韩| 伦理电影大哥的女人| 精品人妻视频免费看| a级毛片a级免费在线| 精品人妻熟女av久视频| 国产精品电影一区二区三区| 午夜a级毛片| 久9热在线精品视频| 亚洲一区二区三区不卡视频| 国产精品免费一区二区三区在线| 中文字幕精品亚洲无线码一区| 精品福利观看| 美女xxoo啪啪120秒动态图 | 亚州av有码| 久久久久性生活片| 身体一侧抽搐| 嫩草影院入口| 中文字幕av在线有码专区| 国产野战对白在线观看| 免费大片18禁| 久久这里只有精品中国| 国产亚洲欧美在线一区二区| 97超级碰碰碰精品色视频在线观看| 91九色精品人成在线观看| 美女免费视频网站| 久久午夜福利片| 久久伊人香网站| 99国产精品一区二区三区| 在线观看一区二区三区| 91久久精品国产一区二区成人| 精品人妻一区二区三区麻豆 | 久久人人爽人人爽人人片va | 国产精品电影一区二区三区| 国产爱豆传媒在线观看| 亚洲av成人av| 一夜夜www| 中文字幕精品亚洲无线码一区| 成人欧美大片| 91久久精品电影网| 国产一区二区在线av高清观看| 最近在线观看免费完整版| 久久人妻av系列| 在线a可以看的网站| 国产午夜精品久久久久久一区二区三区 | 91麻豆精品激情在线观看国产| 一进一出抽搐动态| 伦理电影大哥的女人| 国产成人福利小说| 亚洲精品乱码久久久v下载方式| 国产精品一区二区性色av| 一夜夜www| 午夜精品一区二区三区免费看| 午夜福利18| 一本久久中文字幕| 久久精品国产亚洲av涩爱 | 国产免费男女视频| 国产成年人精品一区二区| 97超级碰碰碰精品色视频在线观看| 国产精品1区2区在线观看.| 琪琪午夜伦伦电影理论片6080| 中出人妻视频一区二区| av在线蜜桃| 性色av乱码一区二区三区2| av在线老鸭窝| 国内揄拍国产精品人妻在线| 欧美色欧美亚洲另类二区| 中文亚洲av片在线观看爽| 中文字幕熟女人妻在线| 国产综合懂色| 在线观看一区二区三区| 欧美丝袜亚洲另类 | 亚州av有码| 老司机福利观看| 欧美+亚洲+日韩+国产| 99久久无色码亚洲精品果冻| 国产成人啪精品午夜网站| 亚洲美女搞黄在线观看 | 久久久久久九九精品二区国产| 国产一区二区在线av高清观看| 国产美女午夜福利| 国产精品乱码一区二三区的特点| 国产伦精品一区二区三区视频9| 国产精品嫩草影院av在线观看 | 欧美一级a爱片免费观看看| 一个人看的www免费观看视频| 久久精品国产清高在天天线| 能在线免费观看的黄片| 免费一级毛片在线播放高清视频| 国产爱豆传媒在线观看| 丰满人妻一区二区三区视频av| 欧美性感艳星| 国内精品久久久久久久电影| 18+在线观看网站| 久久热精品热| 一区福利在线观看| 老鸭窝网址在线观看| 日本三级黄在线观看| 老熟妇乱子伦视频在线观看| eeuss影院久久| 亚洲精品影视一区二区三区av| 日韩欧美国产在线观看| 国产伦一二天堂av在线观看| 丁香欧美五月| 丝袜美腿在线中文| 国产探花极品一区二区| 欧美一区二区亚洲| 99久久99久久久精品蜜桃| 国产欧美日韩精品亚洲av| 高清在线国产一区| 草草在线视频免费看| 最新在线观看一区二区三区| 午夜福利在线在线| 99在线视频只有这里精品首页| 一级黄色大片毛片| 欧美日韩乱码在线| 成熟少妇高潮喷水视频| 欧美日韩瑟瑟在线播放| 91av网一区二区| 欧美在线一区亚洲| 亚洲精品乱码久久久v下载方式| 久久香蕉精品热| 亚洲av成人av| 国产精品美女特级片免费视频播放器| 免费av观看视频| 身体一侧抽搐| 最近视频中文字幕2019在线8| 又爽又黄无遮挡网站| 精品福利观看| 亚洲,欧美精品.| 可以在线观看的亚洲视频| 91在线观看av| 久久精品久久久久久噜噜老黄 | 欧美日韩中文字幕国产精品一区二区三区| 在线播放国产精品三级| 草草在线视频免费看| 老司机福利观看| 精品久久久久久,| 免费在线观看亚洲国产| 午夜免费成人在线视频| xxxwww97欧美| 久久久久亚洲av毛片大全| 99精品久久久久人妻精品| 午夜福利在线观看免费完整高清在 | 免费一级毛片在线播放高清视频| 女生性感内裤真人,穿戴方法视频| 国产成人a区在线观看| 香蕉av资源在线| 亚洲男人的天堂狠狠| 亚洲乱码一区二区免费版| 国产麻豆成人av免费视频| 每晚都被弄得嗷嗷叫到高潮| 午夜精品久久久久久毛片777| 欧美三级亚洲精品| 午夜福利视频1000在线观看| 亚洲电影在线观看av| 亚洲精华国产精华精| 日本熟妇午夜| 91久久精品国产一区二区成人| 亚洲va日本ⅴa欧美va伊人久久| a级毛片a级免费在线| av专区在线播放| 亚洲熟妇中文字幕五十中出| 欧美一级a爱片免费观看看| 黄片小视频在线播放| 国产精品永久免费网站| 久久精品国产清高在天天线| 久久人妻av系列| 亚洲av.av天堂| 国产av一区在线观看免费| 国产精品av视频在线免费观看| 色av中文字幕| 国产高清激情床上av| 内地一区二区视频在线| 色综合亚洲欧美另类图片| 欧美国产日韩亚洲一区| 免费av毛片视频| 亚洲美女搞黄在线观看 | 午夜影院日韩av| 免费在线观看亚洲国产| 久久人人精品亚洲av| 久久99热6这里只有精品| 午夜福利在线在线| 如何舔出高潮| avwww免费| 亚洲avbb在线观看| 草草在线视频免费看| 久久精品国产亚洲av天美| 搡女人真爽免费视频火全软件 | 日本 欧美在线| 日韩精品中文字幕看吧| 嫩草影院入口| 亚洲中文字幕日韩| 最近在线观看免费完整版| 人妻夜夜爽99麻豆av| 最近在线观看免费完整版| 三级国产精品欧美在线观看| 18禁黄网站禁片午夜丰满| 好男人电影高清在线观看| 久99久视频精品免费| 中文字幕av成人在线电影| 99精品久久久久人妻精品| 欧美一区二区精品小视频在线| 可以在线观看的亚洲视频| 日韩人妻高清精品专区| 亚洲成a人片在线一区二区| 国产三级在线视频| 精品日产1卡2卡| 最新中文字幕久久久久| 亚洲天堂国产精品一区在线| 一二三四社区在线视频社区8| 国产探花在线观看一区二区| 一个人免费在线观看的高清视频| 搞女人的毛片| 一个人观看的视频www高清免费观看| 激情在线观看视频在线高清| 在线观看午夜福利视频| 不卡一级毛片| 亚洲人与动物交配视频| 亚洲美女黄片视频| 精品久久久久久成人av| 亚洲,欧美精品.| 天堂网av新在线| 有码 亚洲区| 国产精品人妻久久久久久| 俄罗斯特黄特色一大片| 少妇熟女aⅴ在线视频| 制服丝袜大香蕉在线| 国产探花极品一区二区| 日本 欧美在线| 中文亚洲av片在线观看爽| 亚洲最大成人手机在线| 又粗又爽又猛毛片免费看| 国产精品美女特级片免费视频播放器| 国产野战对白在线观看| 亚洲久久久久久中文字幕| 国产主播在线观看一区二区| 九色国产91popny在线| 国产私拍福利视频在线观看| 亚洲男人的天堂狠狠| 成年女人看的毛片在线观看| 伦理电影大哥的女人| 国产精品一区二区三区四区久久| 亚洲男人的天堂狠狠| 日本精品一区二区三区蜜桃| 免费av观看视频| 亚洲精品色激情综合| 精品午夜福利视频在线观看一区| 日本免费一区二区三区高清不卡| 久久午夜福利片| 欧美激情国产日韩精品一区| 亚洲熟妇中文字幕五十中出| 国产成人欧美在线观看| 动漫黄色视频在线观看| 天美传媒精品一区二区| 看片在线看免费视频| 欧美黄色片欧美黄色片| 日韩欧美国产在线观看| 亚洲片人在线观看| 国产在线精品亚洲第一网站| 成年女人永久免费观看视频| 好看av亚洲va欧美ⅴa在| 51午夜福利影视在线观看| 国语自产精品视频在线第100页| 久久久久九九精品影院| 91午夜精品亚洲一区二区三区 | 中文在线观看免费www的网站| 国产高清视频在线观看网站| 国产精品人妻久久久久久| 精品免费久久久久久久清纯| 99视频精品全部免费 在线| 午夜精品久久久久久毛片777| 在线a可以看的网站| 禁无遮挡网站| 91午夜精品亚洲一区二区三区 | 天天一区二区日本电影三级| a在线观看视频网站| 国产精品综合久久久久久久免费| 色视频www国产| 男人舔女人下体高潮全视频| 亚洲欧美日韩高清在线视频| 最近视频中文字幕2019在线8| 欧美激情久久久久久爽电影| 午夜福利视频1000在线观看| 人人妻,人人澡人人爽秒播| 国产精品久久久久久久久免 | 老司机福利观看| 国产私拍福利视频在线观看| 波多野结衣巨乳人妻| 亚洲av第一区精品v没综合| .国产精品久久| 久久久国产成人免费| 日本与韩国留学比较| 别揉我奶头 嗯啊视频| 亚洲av成人av| 欧美绝顶高潮抽搐喷水| 午夜免费激情av| 中文字幕av成人在线电影| 国产v大片淫在线免费观看| 高清日韩中文字幕在线| 波多野结衣巨乳人妻| 欧美成人性av电影在线观看| 精品一区二区三区视频在线观看免费| 亚洲人成伊人成综合网2020| 男女做爰动态图高潮gif福利片| 免费观看精品视频网站| 性色av乱码一区二区三区2| 亚洲人成网站在线播放欧美日韩| 国产成人av教育| 少妇被粗大猛烈的视频| 国产成人a区在线观看| 欧美性猛交╳xxx乱大交人| 国产精品98久久久久久宅男小说| 熟女人妻精品中文字幕| 婷婷色综合大香蕉| 午夜精品一区二区三区免费看| 欧美性感艳星| 99视频精品全部免费 在线| 无人区码免费观看不卡| 搡女人真爽免费视频火全软件 | 久久热精品热| 色综合欧美亚洲国产小说| 丰满人妻一区二区三区视频av| 人人妻人人看人人澡| 啪啪无遮挡十八禁网站| 99久久久亚洲精品蜜臀av| 一级av片app| 宅男免费午夜| 日日摸夜夜添夜夜添av毛片 | 精品久久久久久久久久久久久| 国产一区二区亚洲精品在线观看| 黄色配什么色好看| 欧美成狂野欧美在线观看| 亚洲乱码一区二区免费版| 国产在视频线在精品| 90打野战视频偷拍视频| 亚洲精品一卡2卡三卡4卡5卡| 男插女下体视频免费在线播放| 中文亚洲av片在线观看爽| 给我免费播放毛片高清在线观看| 亚洲国产日韩欧美精品在线观看| 久久久久久久久久成人| 国产一区二区三区在线臀色熟女| 怎么达到女性高潮| 久久亚洲精品不卡| 精品人妻一区二区三区麻豆 | 可以在线观看的亚洲视频| 亚洲成人精品中文字幕电影| 能在线免费观看的黄片| 亚洲真实伦在线观看| 极品教师在线视频| 久久99热这里只有精品18| 又爽又黄a免费视频| 99久久精品热视频| 欧美国产日韩亚洲一区| 一进一出好大好爽视频| av黄色大香蕉| 91麻豆精品激情在线观看国产| 中文字幕人妻熟人妻熟丝袜美| 我的女老师完整版在线观看| 久久国产精品人妻蜜桃| 亚洲精品色激情综合| 免费在线观看日本一区| 国产精品三级大全| 国产亚洲精品综合一区在线观看| 嫩草影院精品99| 欧美绝顶高潮抽搐喷水| 亚洲av熟女| 变态另类成人亚洲欧美熟女| 最近最新免费中文字幕在线| 国产视频一区二区在线看| 国产精华一区二区三区| 我的老师免费观看完整版| 简卡轻食公司| 亚洲精品在线观看二区| 久久久国产成人精品二区| 国产一区二区在线观看日韩| www.熟女人妻精品国产| 99热只有精品国产| 中文字幕久久专区| 欧美性猛交黑人性爽| 亚洲,欧美,日韩| 啪啪无遮挡十八禁网站| 全区人妻精品视频| 国产真实乱freesex| 亚洲精品在线观看二区| 一级av片app| 久久这里只有精品中国| 可以在线观看的亚洲视频| 国产久久久一区二区三区| 99热只有精品国产| 国产在线男女| 亚洲av一区综合| 免费无遮挡裸体视频| 欧美成狂野欧美在线观看| 欧美又色又爽又黄视频| 欧美日韩亚洲国产一区二区在线观看| 美女黄网站色视频| 亚洲av中文字字幕乱码综合| 国产精品精品国产色婷婷| 色5月婷婷丁香| 一夜夜www| 亚洲电影在线观看av| 99精品久久久久人妻精品| 久久人人精品亚洲av| 别揉我奶头 嗯啊视频| 亚洲熟妇熟女久久| av专区在线播放| 免费观看人在逋| 国产乱人伦免费视频| 国产精品嫩草影院av在线观看 | a级一级毛片免费在线观看| 男女之事视频高清在线观看| 欧美不卡视频在线免费观看| 18美女黄网站色大片免费观看| 亚洲人成网站在线播放欧美日韩| 国产精品美女特级片免费视频播放器| 人人妻人人澡欧美一区二区| 婷婷丁香在线五月| 亚洲av.av天堂| 99久久精品一区二区三区| 免费观看人在逋| 国产亚洲精品久久久久久毛片| 国产成人aa在线观看| 3wmmmm亚洲av在线观看| 欧美一级a爱片免费观看看| 听说在线观看完整版免费高清| 精品福利观看| 99在线视频只有这里精品首页| 动漫黄色视频在线观看| 一级作爱视频免费观看| 亚洲真实伦在线观看| 校园春色视频在线观看| 宅男免费午夜| 国产亚洲精品av在线| 久久午夜福利片| 91字幕亚洲| 99热这里只有精品一区| 久久人人爽人人爽人人片va | 国产免费av片在线观看野外av|