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

    Prussian Blue Analogues in Aqueous Batteries and Desalination Batteries

    2021-10-21 03:30:42ChiweiXuZhengweiYangXikunZhangMaotingXiaHuihuiYanJingLiHaoxiangYuLiyuanZhangJieShu
    Nano-Micro Letters 2021年11期

    Chiwei Xu ,Zhengwei Yang ,Xikun Zhang ,Maoting Xia ,Huihui Yan ,Jing Li ,Haoxiang Yu,Liyuan Zhang,Jie Shu

    ABSTRACT In the applications of large-scale energy storage,aqueous batteries are considered as rivals for organic batteries due to their environmentally friendly and low-cost nature.However,carrier ions always exhibit huge hydrated radius in aqueous electrolyte,which brings difficulty to find suitable host materials that can achieve highly reversible insertion and extraction of cations.Owing to open threedimensional rigid framework and facile synthesis,Prussian blue analogues (PBAs)receive the most extensive attention among various host candidates in aqueous system.Herein,a comprehensive review on recent progresses of PBAs in aqueous batteries is presented.Based on the application in different aqueous systems,the relationship between electrochemical behaviors (redox potential,capacity,cycling stability and rate performance) and structural characteristics (preparation method,structure type,particle size,morphology,crystallinity,defect,metal atom in highspin state and chemical composition) is analyzed and summarized thoroughly.It can be concluded that the required type of PBAs is different for various carrier ions.In particular,the desalination batteries worked with the same mechanism as aqueous batteries are also discussed in detail to introduce the application of PBAs in aqueous systems comprehensively.This report can help the readers to understand the relationship between physical/chemical characteristics and electrochemical properties for PBAs and find a way to fabricate high-performance PBAs in aqueous batteries and desalination batteries.

    KEYWORDS Prussian blue analogues;Electrochemistry;Cathode;Aqueous batteries;Desalination batteries

    1 Introduction

    As is known to all,energy is the source of economic development and the material basis for human survival.Therefore,energy has always been a common concern of the world.Traditional energy sources have brought many changes to the development of human society and the world.However,the excessive use of traditional energy produces carbon dioxide and other pollutants,causing serious ecological problems.To solve these environmental problems,scientists are beginning to explore renewable energy sources to reduce environmental pollution and meet the needs of social development[1-5].However,these clean energy sources generally have the disadvantages of location dependence [6-10],high cost and low efficiency,which restrict their large-scale applications [11-16].In this case,a lot of attentions have been paid to the developments of rechargeable batteries [17-22],fuel cells [23-29] and supercapacitors [30-35],among which the rechargeable batteries have received the most attention.For rechargeable batteries,most of them are powered by organic electrolytes due to high voltage,high energy density and wide electrochemical window of organic electrolyte systems[36-41].However,the volatilization,toxicity and flammability of organic electrolytes have brought many environmental and safety problems to their widespread applications [42,43],especially in large-scale electrical storage applications.To get better batteries,much effort has been paid to aqueous electrolytes due to their advantages of high safety,environmental friendliness and low cost [44-46].

    Aside from electrolytes,the selection of electrode material also occupies an important position in the application of energy storage devices.In various types of electrode materials,tunnel-and layered-type transition metal oxides [47-53],polyanionic compounds [54-59],PBAs [60-70] have been widely studied.In recent researches,the most commonly used cathode materials are metal oxides because of their high theoretical capacity.However,this material has two disadvantages.One is that the synthesis method is complex,and the other is that many metal oxides are unable to resist the erosion of corrosive electrolytes,which are usually characterized by low rate capacity and poor structural stability.Polyanionic compounds as typical electrode materials possess high safety,low cost and environmental friendliness.However,their poor electrochemical properties restrict the practical application.PBAs are blessed with the advantages of acceptable theoretical capacity,high stability and low cost.In particular,many PBAs can be synthesized by simple methods,and nano-structured materials are usually obtained by the co-precipitation method.Therefore,these materials have great application prospects in rechargeable batteries.In the early eighteenth century,the PB was first synthesized in Berlin,but is initially used as a pigment.With the development of energy storage technology,the requirement of the low cost of the electrode materials promotes the application of PB in the battery industry.In 1978,Neff et al.first explored the reversible inserted/extracted electrochemical property of K+in K2FeIIFeII(CN)6and opened the door of PBAs as electrode material [71].Then,in 1999,the various types of PBAs used as electrode in organic Li+battery are reported [72,73].However,the extensive research of PBAs as electrode material in aqueous batteries began to rise after the twenty-first century.In 2012,Cui et al.proved that Li+,Na+,K+and NH4+exhibit reversible inserted/extracted ability in K0.9Cu1.3Fe(CN)6and K0.6Ni1.2Fe(CN)6,respectively,which inspires the research interest for PBAs in aqueous battery [74].Aside from these monovalent ions,they later demonstrated that the multivalent ions (Zn2+,Mg2+,Al3+,etc.) can also be reversibly inserted/extracted in PBAs [75,76].Especially for aqueous multivalent-ion batteries,Zhi et al.recently found that CoFe(CN)6shows the high capacity and long cycle performance in aqueous Zn2+batteries,which further promoted the research of PBA in aqueous multivalent ions batteries [77].Meanwhile,the application of PBAs in desalination battery as a kind of aqueous battery also has attracted more and more attention in recent years;particularly,it has great development potential in realizing the nonmembrane of desalination battery [78].To understand the development process of PBAs in aqueous batteries more intuitively,we also summarized other milestone applications of PBAs in aqueous batteries in the timeline in Fig.1.

    Fig.1 A summary of the key milestone in the development of PBAs used as electrode material in aqueous batteries

    Obviously,PBAs have been widely studied as battery electrode material for some years [71-81],so there are some relevant reviews on the application of PBAs in the battery field.For example,Li and co-workers introduced the application of PBAs in sodium-ion and potassium-ion batteries in detail,and some strategies for performance optimization are proposed [82].Srinivasan discussed the application of various electrode materials,including PBAs,in aqueous Zn2+and Al3+batteries [83].Ji and co-workers introduced the PBAs as electrodes for aqueous monovalent-ion batteries [84].Zhi and co-workers analyzed the PBAs as host materials for various carriers in both organic and aqueous electrolyte [85].However,none of these reviews systematically discusses the application of PBAs in aqueous batteries.In this review,we comprehensively introduce the application of PBAs in aqueous monovalent-ion batteries and multivalent-ion batteries.Impressively,the electrochemical properties of various PBAs and their structural characteristics are discussed.Meanwhile,various aqueous desalination batteries with unique characteristics are also introduced because they are working with the same mechanism as aqueous batteries.Therefore,this paper has a high reference value for future research on the PBAs used as electrode materials in aqueous batteries and desalination batteries.

    2 PBAs in Aqueous Batteries

    With a metal-organic open framework,PBAs possess a general composition of AxM[R(CN)6]·wH2O in which A is the alkali metals and M and R represent the transition metal.It is worth noting that most of the R-site in PBAs for the aqueous battery is replaced by Fe atom.Here,the classic PBAs structure is illustrated in Fig.2.In this facecentered cubic structure,the nitrogen and carbon atoms from the CN ligands combine with the M and R atom,respectively,creating open ionic channels and spacious interstitial spaces [86].As electrode material used in aqueous batteries,PBAs release and store energy through the insertion and extraction of host ions in PBAs.Here,AxM[R(CN)6]·wH2O is taken for example as a cathode.During the discharging process,host ions in the electrolyte are inserted into PBAs,which causes the R or M and R atoms to simultaneously undergo a reduction reaction to release energy.In the charging process,host ions are extracted from PBAs with R or M and R atoms simultaneous undergoing oxidation reaction to store energy.The charging and discharging process can be reflected by the following chemical reaction equation:

    Fig.2 a Schematic crystal structure of PBAs frameworks:An intact AxM[R(CN)6] framework without structural defects.b A defective AxM[R(CN)6] framework with vacancies and interstitial water existing in each unit cell

    Furthermore,Fig.3a also shows the insertion and extraction process of host ions in PBAs.For the aqueous battery,host ions will undergo a hydration process before entering into PBAs and then dehydrate after entering into PBAs.As shown in Fig.3b,the hydration ion radius will be larger than the original ion radius after hydration,which puts forward higher requirements for the structural stability of PBAs.

    Fig.3 a Schematic of host ion insertion and extraction in PBAs.b Histogram of ionic radius and hydration radius of various host ions

    For PBAs,it possesses the 3D cubic structure and tunable chemical compositions that give PBAs some advantages when inserting ions with large hydrated ionic radius.Firstly,there are many available “A” sites with large interstitial (~ 4.6 ? diameter) in the body center of PBAs lattice.Meanwhile,the PBAs possess a special ions channel with approximately 3.2 ? diameter in the (100) direction [87,88].The large insertion site and ion channel ensure the fast reaction kinetics of most host ions with small hydrated radius when they are inserted/extracted in/from PBAs,which usually brings well cycling stability to PBAs [89,90].Secondly,PBAs exhibit multivalent redox-active centers such as the M3+/M2+and Fe3+/Fe2+couples,which means that PBAs can provide more capacity.The theoretical capacities of different PBAs in characteristic aqueous systems are summarized in Fig.4.For FeFe-PB,CoFe-PBAs and MnFe-PBAs,their theoretical capacities are all over 150 mAh g-1.But the theoretical capacity utilization of NiFe-PBAs,CuFe-PBAs and ZnFe-PBAs is only 50% because Ni,Cu and Zn in PBAs are electrochemically inert.In addition to these structural advantages,PBAs also have some structural defects that need to be addressed.During the preparation process of PBAs through co-precipitation,it will inevitably bring some interstitial water and vacancies in PBAs crystals,and Fig.2 exhibits the crystal structure diagram of PBAs with vacancies and interstitial water.For interstitial water,it will cause the following two disadvantages to PBAs.First of all,the presence of interstitial water will occupy the insertion site of host ions,resulting in a capacity decrease of PBAs [91,92].Secondly,interstitial water in PBAs may enter into the electrolyte in the process of ion insertion and extraction,resulting in the decrease of electrolyte concentration that will aggravate the dissolution of PBAs structure.For Fe(CN)6vacancies,its random distribution can cause the distortion of PBAs lattice,thus reducing the structural stability.Especially for reduced-form PBAs,due to the oxidation reaction of transition metal elements,the structural instability is more serious.Furthermore,as shown in Fig.5a,in the process of host ions insertion/extraction in PBAs,its structure will change from cubic structure to monoclinic structure,which reduces the stability of PBAs.Generally,it is concluded that the host ions-deficient PBAs present the cubic phase,while host ions-rich PBAs usually present the monoclinic phase.Impressively,the asymmetrical occupation of electrons from transition metal atom in the degenerate orbitals will produce John-Teller effect in PBAs structure.This phenomenon will cause different shielding effects of electrons in different directions on central ions,thus distorting the geometric configuration of PBAs molecules,reducing the symmetry and degeneracy of the orbitals and finally reducing the energy of the system.Figure 5b exhibits the phase transitions between monoclinic,cubic and rhombohedral structures caused by John-Teller distortion.Despite the presence of these vacancies and interstitial water in PBAs,a large number of studies have shown that their content can be well controlled by a controlled crystallization method.Besides,the John-Teller distortion can also be inhibited by partial substitution of the transition metal attached to N by other metal ion.Thus,PBAs are suitable electrode materials for aqueous batteries.Herein,for the sake of better introducing the application of PBAs in aqueous batteries,we review the recent advances of PBAs in aqueous monovalent-ion batteries and multivalention batteries,respectively.

    Fig.4 Theoretical capacity of different PBAs in various aqueous batteries (the shaded part represents the un-activated capacity of the singleatom redox PBAs)

    Fig.5 Structure transformation of PBAs a between monoclinic and cubic structures and b among monoclinic,cubic and rhombohedral structures

    2.1 PBAs for Aqueous Monovalent-Ion Batteries

    The application of PBAs in aqueous monovalent-ion batteries can be divided into different fields:aqueous Na+,Li+,K+,H+and NH4+batteries.For aqueous Na+and K+batteries,PBAs usually show excellent cyclic stability due to the better structural compatibility and moderate hydrated radius for Na+and K+.Usually,these PBAs in aqueous Na+and K+batteries can maintain more than 80% of the initial capacity within 300 cycles.Although PBAs have poor cyclic stability in aqueous Li+batteries due to the large hydration ion radius of Li+(higher than Na+,K+,H+,and NH4+),mixed-ion batteries composed of Li+and other metal ions have great application prospects in the improvement in energy density.Besides,the excellent rate performance and cycling stability of PBAs in aqueous H+batteries and aqueous NH4+batteries have also attracted more and more attention in recent years due to the special transport mechanism of H+in aqueous electrolyte and the good structural compatibility of NH4+and PBAs.This section will comprehensively introduce the application of PBAs in these aqueous monovalent-ion batteries.

    2.1.1 PBAs for Aqueous Na+ Batteries

    For the PBAs in aqueous Na+batteries (ASIB),they can be divided into single-metal-atom redox PBAs (NiFe-PBAs,CuFe-PBAs and ZnFe-PBAs) and double-metal-atom redox PBAs (FeFe-PB,CoFe-PBAs and MnFe-PBAs).The singlemetal-atom redox PBAs mean that only Fe atom can undergo the redox reaction and provide capacity during the insertion and extraction of Na+.Therefore,as shown in Fig.4,the theoretical capacity utilization of the single-metal-atom redox PBAs is only 50%.The double-metal-atom redox PBAs means that both Fe atom and M atom can undergo redox reaction and provide capacity during the insertion and extraction of Na+.Herein,for better comparison,PBAs in ASIB are cataloged into two groups:PBAs with single-atom redox and PBAs with double-atom redox.

    2.1.1.1 PBAs with Single-Atom RedoxFrom the above introduction,when the Ni,Cu,and Zn are located at M-site in AxM[Fe(CN)6]y·nH2O,these types of PBAs belong to the single-metal-atom redox PBAs.The typical single-metal-atom redox nickel hexacyanoferrate (K0.6Ni1.2Fe(CN)6·3.6H2O,NiFe-PBA) was firstly synthesized by Cui and co-workers through spontaneous precipitation in aqueous solution as electrode materials for ASIB [43].As only Fe3+/Fe2+redox couple possess electrochemical activity,its discharge capacity only has 60 mAh g-1at the current rate of approximately 0.8 C.However,K0.6Ni1.2Fe(CN)6·3.6H2O exhibits excellent electrochemical stability with almost no capacity fading at 8.3 C after 5000 cycles due to the small structural changes and structural stress-strain in the Na+insertion/extraction process.But it is in a Na-deficient state,making it impossible to construct a practical ASIB with conventional Na-free anodes [93].Generally speaking,a Na insertion cathode should be designed in a Na-rich state (a discharged state),so as to act as a Na+reservoir to provide removable Na+for the Na+-deficient anode,thus enabling a rocking-chair Na-ion battery.To address this problem,a Na-rich Na1.94Ni1.03Fe(CN)6·4.8H2O was synthesized by Yang and co-workers and construct a practical ASIB with NaTi2(PO4)3[94].

    Later,Zhu and co-workers also synthesized Na-rich Na1.45Ni[Fe(CN)6]0.87·3.02H2O and construct a practical ASIB with NaTi2(PO4)3(Fig.6a).Besides,they find monoclinic Na1.45Ni[Fe(CN)6]0.87·3.02H2O has advantages in reversible capacity and cycle stability compared with cubic Na1.21Ni[Fe(CN)6]0.86·3.21H2O [95].The reasons can be concluded into three points.Firstly,the monoclinic crystal structure possesses more sodium that makes it contain fewer crystal defects in the lattice.Secondly,during the synthesis of monoclinic Na1.45Ni[Fe(CN)6]0.87·3.02H2O,the addition of chelating agents and surfactants reduces the crystallization rate and makes the crystal shape of monoclinic Na1.45Ni[Fe(CN)6]0.87·3.02H2O more uniform,thus increasing the contact area between the electrode and the electrolyte and shortening the migration path of ions.Thirdly,a slower reaction rate leads to the formation of higher crystallinity product with less interstitial water and vacancy,thus accelerating Na+insertion/extraction rate.Moreover,an ex situ X-ray photoelectron spectroscope (XPS) at different charge/discharge states in Fig.6b explored the valance state change of Fe and Ni atoms during Na+extraction/insertion.The result shows that the valance change only occurs in Fe atom,which fully explained the Ni atom in NiFe-PBAs with no electrochemical activity.Although NiFe-PBAs show acceptable performance in terms of electrochemical stability,its low redox potential is not suitable for aqueous energy storage applications.

    Fig.6 a Schematic illustration of Na1.45Ni[Fe(CN)6]0.87·3.02H2O/NaTi2(PO4)3 full cell.b Ex situ XPS of Fe 2p (left) and Ni 2p (right) spectra.Reproduced with permission from Ref.[95].Copyright 2020,Elsevier.c Schematic illustration of sodium-ion intercalation mechanism in Zn3[Fe(CN)6]2.d SEM images of the pristine Zn3[Fe(CN)6]2.Reproduced with permission from Ref.[100].Copyright 2018,Elsevier.e Schematic fabrication process of K2Zn3[Fe(CN)6]2·9H2O on carbon cloth.f SEM image of K2Zn3[Fe(CN)6]2·9H2O.Reproduced with permission from Ref.[101].Copyright 2019,Wiley-VCH

    For the redox potential (E) of host ions in PBAs,there is a qualitative relationship [96]:

    The value ofEredoxis related to the potential of the redox couple itself.Eion-hostis related to the interaction between host ions and PBAs framework,which usually increases as the radius of the ion enlarges.Eion-ionis related to the repulsion between host ions in PBAs.Since the redox couple of singlemetal-atom redox PBAs is only [Fe(CN)6]3-/[Fe(CN)6]4-,the redox potential will not differ much between them when the same ion is inserted in PBAs.However,the increase of electrons in the antibonding orbital (eg) of M atom in PBAs can cause [Fe(CN)6]3-/[Fe(CN)6]4-redox potential to move toward more positive values [97].Thus,the redox potentialvs.standard hydrogen electrode (SHE,all the potentials in this manuscript are reported versus the SHE) of these single-metal-atom redox PBAs follows the order:ZnFe-PBA > CuFe-PBA > NiFe-PBA.Obviously,CuFe-PBA and ZnFe-PBA are more suitable as cathode material compared with NiFe-PBA in terms of redox potential.The typical Narich Na2CuFe(CN)6with the insertion potential of 0.82 V was synthesized by Wu and co-workers,and a practical ASIB was constructed with NaTi2(PO4)3[98].Impressively,the Na2CuFe(CN)6-NaTi2(PO4)3full cell shows high operating voltage (1.4 V) due to the high redox potential of Na2CuFe(CN)6.Combining the high stability of NiFe-PBA with high redox potential of CuFe-PBA,Zhang and co-workers synthesized Ni-substituted copper hexacyanoferrates(Na2Cu1-xNix[Fe(CN)6]) as cathode for ASIB [99].For this Na2Cu1-xNix[Fe(CN)6] cathode,it can exhibit a discharge capacity of 56 mAh g-1with a capacity retention of 96%after 1000 cycles.Impressively,its redox potential can be adjusted from 0.6 to 1.0 V with the increased content of Cu.

    For ZnFe-PBAs,it has great potential in improving energy density due to its high redox potential.Recently,Liu and co-workers reported a full cell that uses Zn3[Fe(CN)6]2as cathode,NaTi2(PO4)3as anode,and NaClO4-H2O-polyethylene glycol as electrolyte [100].Figure 6c,d shows the schematic illustration of Na+intercalation mechanism in Zn3[Fe(CN)6]2and its basic morphology.For this battery,its energy density can reach up to 59 Wh kg-1due to the high output voltage of 1.6 V.Besides,the energy density of ZnFe-PBAs can be further improved by improving the synthesis method.Yao and co-workers adopt a two-step method of electrodeposition and water bath treatment grew K2Zn3[Fe(CN)6]2·9H2O on carbon cloth (CC@ZnFe-PBAs)[101].The schematic fabrication process is exhibited in Fig.6e,and the as-formed K2Zn3[Fe(CN)6]2·9H2O cubes are uniformly deposited on carbon cloth (Fig.6f).Since the CC@ZnFe-PBA electrode avoids the use of inactive substances (conductive additive,binder),it exhibits a high energy density of 14.3 mW cm-2and high output voltage of 1.6 V when it is combined with NaTi2(PO4)3to form a full cell.In general,the single-metal-atom redox PBAs inherit the advantage of high electrochemical stability that makes it possible to be used as ASIB electrode.However,considering the specific capacity,double-mental-atom redox PBAs are more competitive.

    2.1.1.2 PBAs with Double-Atom RedoxObviously,single-metal-atom redox PBAs have a disadvantage that their capacity is not high.The essential cause of this problem is that only one redox-active couple (Fe3+/Fe2+) releases its redox capacity.In contrast,double-metal-atom redox PBAs have the advantage in capacity because both Fe3+/Fe2+and M3+/M2+can provide capacity.As typical double-metalatom redox PBAs,NaFeFe(CN)6[102] and Co3[Fe(CN)6]2[103] are initially proposed as host materials for Na+storage.However,its specific capacity is approaching 70 mAh g-1which is much less than the theoretical capacity (Fig.4).Through reviewing the recent advances of PBAs in ASIB,there are four main reasons for the low capacity of doublemetal-atom redox PBAs.Firstly,the disadvantage of the synthesis method results in a lot of Fe(CN) vacancies.Secondly,more interstitial water in PBAs hinders Na+insertion.Thirdly,the low concentration of electrolyte results in a narrow voltage window,which affects the charge/discharge capacity.Finally,M atoms connecting with N atoms have weaker electrochemical activity.To a large extent,the above four factors limit the released capacity of double-metalatom redox PBAs.

    For the control of defects and the amount of interstitial water,Yang and co-workers synthesized Na1.33Fe[Fe(CN)6]0.82with low defect and water content by a multi-step crystallization method [104].The cubic morphology is displayed in Fig.7a.In comparison with hydrated FeFe-PBA,Na1.33Fe[Fe(CN)6]0.82exhibits a higher specific capacity of 125 mAh g-1at 2 C.Even at 20 C,it still can present a highly efficient Na storage reversibility with a desirable capacity of 102 mAh g-1(Fig.7b).Later,Yang adopted a similar crystallization method as shown in Fig.7c to achieve vacancy-free Na2CoFe(CN)6(Na1.85Co[Fe(CN)6]0.99·2.5H2O) [105],which exhibits perfect crystal structure without vacancy in Fig.7d.Because of the vacancy-free structure and two redox centers,the capacity of Na2CoFe(CN)6can approach up to 130 mAh g-1,which is a high value for aqueous battery in existing reports.

    Besides the defects and interstitial water of PBAs,electrolyte concentration is also one of the factors affecting the electrochemical performance of PBAs.Take MnFe-PBAs for example,Okada and co-workers synthesized Na0.13Mn[Fe(CN)6]0.81and explored the influence of electrolyte concentration on electrochemical performance[106].For Na0.13Mn[Fe(CN)6]0.81/NaTi2(PO4)3full battery in 1 M NaClO4and 17 M NaClO4,they find the electrochemical window in 17 M NaClO4is 2.8 V which is higher than the 1.9 V in 1 M NaClO4and the initial charge/discharge capacities can raise to 124/116 mAh g-1(Fig.7e).The reason is that the free water molecules in high-concentration electrolyte are less than dilute concentration electrolyte.Thus,the evolution of oxygen and hydrogen is inhibited largely,and the electrochemical window is enlarged.Besides,low-concentration electrolyte also has a certain influence on the structural stability of double-metal-atom redox PBAs.In dilute electrolytes,the oxygen evolution is more severe,which can increase the amount of OH-in the electrolyte.After that,Mn+dissolved in the electrolyte reacts with OH-to form precipitation and this irreversible reaction will aggravate the structural collapse of double-metal-atom redox PBAs.For Na0.13Mn[Fe(CN)6]0.81,the structural collapse can be described in the following chemical reaction:

    Fig.7 a TEM images of as-prepared Na1.33Fe[Fe(CN)6]0.82 nanocrystals.b Rate performance of Na1.33Fe[Fe(CN)6]0.82.Reproduced with permission from Ref.[104].Copyright 2015,Elsevier.c Schematic representation of the formation mechanism of Na2CoFe(CN)6 in a controlled crystallization reaction.d Crystal structure of Na2CoFe(CN)6.Reproduced with permission from Ref.[105].Copyright 2015,Wiley-VCH.e Cyclability dependence on the current density of Na2MnFe(CN)6 half cell with 17 M NaClO4 aqueous electrolyte.The first and second charge/discharge curves of Na0.13Mn[Fe(CN)6]0.81/NaTi2(PO4)3 full cell and NaTi2(PO4)3//Zn half cell f in 1 M NaClO4 electrolyte and g in 17 M NaClO4 electrolyte.Reproduced with permission from Ref.[106].Copyright 2017,Electrochemical Society.h Crystal structure of manganesecobalt hexacyanoferrate with rhombohedral lattice.Reproduced with permission from Ref.[107].Copyright 2016,The Royal Society of Chemistry

    In addition,according to the crystal field theory,the d orbital is subjected to the electrostatic repulsion of the octahedron crystal field and split into t2gand e2gorbitals.Compared with t2gorbitals,the variation of electrons in the e2gorbital has more influence on the crystal structure due to the orientation of the e2gorbitals along the axis of the bond [97].Because their crystal field fission energy is weaker,these transition metals attached to the N atom tend to gain and lose electrons in the e2gorbital when they undergo redox reactions.In contrast,the transition metal connected to C atom has strong crystal field fission energy,and the gain and loss of electrons occur in the t2gorbital.Therefore,the redox reaction of the transition metal connected with N atom has a great influence on the crystal structure,causing it to contribute less capacity.About the relationship of transition metal atoms connecting with N atoms and electrochemical property of doublemental-atom redox PBAs,Cui et al.conducted a systematic study through Na1.54Co[Fe(CN)6]0.86·γ0.14·2.16H2O and Na1.33Mn[Fe(CN)6]0.79·γ0.21·1.88H2O [107].Figure 7f exhibits the crystal structure of the above two PBAs.They prove that the transition metal connected with C atom plays an important role in the stability of the crystal structure and provides more capacities.On the contrary,the transition metal connected with N atom has worse stability in crystal structure due to the weak N-coordinated crystal field and thus provides fewer capacities.Obviously,it provides a new idea for future research on the mechanism of capacity decline of double-metal-atom redox PBAs.

    2.1.2 PBAs for Aqueous Li+ Batteries

    In aqueous battery system,the host ions inserted in the host material are accompanied by a hydration process in which the free water molecules in the electrolyte will combine with the host ion to form a hydrated ion.After hydration process,the hydrated ionic radius of Li+will get larger.Thus,the insertion performance of Li+in PBAs is usually poor because the large hydrated ionic radius of Li+has a severe impact on the structural stability of PBAs.Earlier research in the insertion of Li+in PBAs also confirmed the theory.In 2012,Cui and co-workers first reported two types of PBAs,K0.9Cu1.3Fe(CN)6and K0.6Ni1.2Fe(CN)6,and studied their ability to store Li+in aqueous electrolyte [74].Both PBAs can host Li+,but both exhibit poor cycling performance.Later,Lee and co-workers reduced the particle size of copper hexacyanoferrate to nanoscale for achieving Li storage [88].But the result is still disappointing.Figure 8a,b exhibits the insertion sites and mechanism of Li+in the copper hexacyanoferrate,respectively.It is apparent that the cause of poor Li+insertion ability in PBAs is that in aqueous electrolyte the radius of hydrated Li+is too large.Thus,it is difficult for PBAs to receive hydrated Li+intercalation in the structure and maintain structural stability during the process of Li+insertion and extraction.

    Although the Li+storage performance of PBAs in aqueous electrolyte is poor,the co-insertion of Li+with K+or Na+in PBAs has received special attention.In 2014,Liu and co-workers used K0.08Ni0.75Zn0.70Fe(CN)6as cathode and TiP2O7as anode explored the co-insertion performance of Li+/K+[108].Figure 8c exhibits the schematic of the battery.Different from PBAs as Li+host material,K0.08Ni0.75Zn0.70Fe(CN)6exhibits good cyclic performance for Li+/K+storage that can keep 63% of the original capacity after 50 cycles (Fig.8d).It is worth noting that this system delivers an average operating voltage of 1.25 V,which is much higher than most PBAs in aqueous electrolyte.However,they did not give systematic explanation for the increase of operating voltage in this work.Later,they reported zero-strain K0.04Cu1.47Fe(CN)6to co-insert Li+/K+and explored this reason,which also exhibits a high operating potential up to 1.39 V [109].Figure 8e exhibits the SEM image (top) and TEM image (below) of K0.04Cu1.47Fe(CN)6.They found that the binding energy of K+inserted in its corresponding redox sites is smaller than Li+.Thus,K0.04Cu1.47Fe(CN)6prefers to insert K+.When Li+and K+co-insert in PBAs,the presence of K+will affect the insertion kinetics of Li+that can shift the operation potential of PBAs to a higher value and this phenomenon can be seen in the galvanostatic profiles of the co-insertion of Li+/K+(Fig.8f).Obviously,the increase of cathode potential can result in the increase of battery operation voltage.For further improving the cycling performance of the mixed-ion battery,they later explored the co-insertion of Li+and Na+in Fe4[Fe(CN)6]3[110].Impressively,it exhibits a wonderful rate performance with no significant capacity decrease from 5 to 2 C and then retains this capacity up to 400 cycles(Fig.8g).In summary,the co-insertion of Li+with Na+or K+in PBAs has a broad prospect.However,there are few studies on the co-insertion of Li+with other ions at present,so the future study of Li+insertion behavior in PBAs can be put in this aspect.

    Fig.8 a Crystal structure of copper hexacyanoferrate.b Schematic of Li+intercalation in copper hexacyanoferrate.Reproduced with permission from Ref.[88].Copyright 2019,Electrochemical Society.c A schematic of K0.08Ni0.75Zn0.70Fe(CN)6/TiP2O7 battery.d Cycle life ofK0.08Ni0.75Zn0.70Fe(CN)6/TiP2O7 battery at 1 C.Reproduced with permission from Ref.[108].Copyright 2014,Wiley-VCH.e SEM and TEM images of K0.04Cu1.47Fe(CN)6.f Galvanostatic profiles of K0.04Cu1.47Fe(CN)6/TiP2O7 with co-insertion of Li+and K+.Reproduced with permission from Ref.[109].Copyright 2017,The Royal Society of Chemistry.g Rate capability and cycling stability of Fe4[Fe(CN)6]3/TiP2O7 full cell.Reproduced with permission from Ref.[110].Copyright 2013,American Chemical Society

    2.1.3 PBAs for Aqueous K+ Batteries

    Due to the small hydration ion radius of K+(Fig.3b),its transport speed in the ion channel will be relatively fast.Besides,the 8c site,as the insertion site of K+in PBAs,has the largest volume among all sites,which will also accelerate the transport of K+in PBAs [96].Thus,PBAs usually have an excellent rate performance in aqueous K+batteries.As the crystal structure displayed in Fig.9a,Cui and co-workers synthesized K0.6Ni1.2Fe(CN)6·3.6H2O by co-precipitation method to insert K+[43].Due to the low-strain characteristic of this K0.6Ni1.2Fe(CN)6·3.6H2O and the characteristic of K+inserted in PBAs,it can retain 66% of the initial capacity even at 41.7 C and its Coulombic efficiency can reach 99%at low rates.Soon after,Zhao and co-workers synthesized nanosized potassium-rich mesoporous K2NiFe(CN)6·1.2H2O(Fig.9b) through a diffuse ion-induced co-precipitation approach and constructed full cell with NaTi2(PO4)3[80].It exhibits ultrafast potassium-ion storage within 4.1 s at 500 C and 98.6% of capacity retention at 30 C after 5000 cycles (Fig.9c).Later,Cui and co-workers proved that K0.9Cu1.3Fe(CN)6(Fig.9d) also possesses outstanding rate and cycle performance when it hosts K+in aqueous electrolyte [74].Although these PBAs both exhibit well rate performance in aqueous K+batteries,they have a common shortcoming that the reversible capacity only has approximately 60 mAh g-1.To achieve high K-uptake capacity,Wang and co-workers synthesized high-potassium content K2FeII[FeII(CN)]6·2H2O nanocubes (Fig.9e,f) as cathode to store K+[7].Due to the presence of two redox centers,it can exhibit a high discharge capacity of 120 mAh g-1and approximately 85% capacity retention over 500 cycles at 21.4 C.Even at higher rates,it still has good electrochemical performance (Fig.9g).

    Obviously,the FeFe-PBAs used as the cathode for aqueous K+batteries possess huge application potential and thus received more attention in further optimizing its electrochemical performance.The most mainstream methods are controlling the size of particle and achieving higher crystallinity.The influence of size of particle on performance can be summarized as follows:Firstly,small crystals provide shorter diffusion paths and thus diffusion rates of alkali ions are improved.Secondly,the strain produced by volume changes during ion charge and discharge processes is minimized,which in turn improves rate capability and stability.Thirdly,the increase of the surface area of PBAs can effectively improve electron transport and structural resilience [111].Solanki and co-workers synthesized different-sized PBAs (KFeIIIFeII(CN)6with size of 50-75 nm,FeIIIFeIII(CN)6with size of 2-10 microns) as shown in Fig.9h and explored the effect of particle size on the properties of PBAs for inserting K+[112].The results show that PBAs with small particle size do have greater advantage in electrochemical performance.For obtaining small particle size,in addition to the widely used controlled crystallization method,Coronado and co-workers propose a new method that uses MoS2as reducing agent and platform to synthesize PBAs shells with small and uniform particle size (Fig.9i)and got excellent electrochemical performance [113].Therefore,in the future preparation of PBAs,supporting materials such as MoS2can be used to provide two-dimensional active carriers to achieve uniform nucleation of nanoparticles.In general,the PBAs with open framework possess huge progress potential in aqueous K+batteries and are waiting to be further explored.

    Fig.9 a Face-centered cubic structure of NiFe-PBA (K0.6Ni1.2Fe(CN)6·3.6H2O).Reproduced with permission from Ref.[43].Copyright 2011,American Chemical Society.b General view of the face-centered cubic crystal structure of K2NiFe(CN)6·1.2H2O.c Long-term cycling performance at 30 C of K2NiFe(CN)6·1.2H2O.Reproduced with permission from Ref.[80].Copyright 2018,Wiley-VCH.d Face-centered cubic structure of K0.9Cu1.3Fe(CN)6.Reproduced with permission from Ref.[74].Copyright 2012,Electrochemical Society.e TEM image of K2FeII[FeII(CN)]6·2H2O nanocubes.f Refined crystal structure of FeFe-PBA.g Long cycle performance of K2FeII[FeII(CN)]6·2H2O nanocubes at different current density.Reproduced with permission from Ref.[7].Copyright 2017,Wiley-VCH.h Crystal lattice unit cell of Prussian blue(left,KFeIIIFeII(CN)6) and crystal lattice unit cell of Prussian green (right,FeIIIFeIII(CN)6) (The black,blue,purple and brown balls represent C,N,K and Fe atoms,respectively).Reproduced with permission from Ref.[112].Copyright 2015,Elsevier.i TEM image showing a MoS2 flake uniformly decorated by Prussian blue nanoparticles to form K0.47Fe4III[FeII(CN)6]3.14@(MoSO1.7)0.44·18H2O.Reproduced with permission from Ref.[113].Copyright 2018,Wiley-VCH

    2.1.4 PBAs for Aqueous H+ Batteries

    At present,most of charge carriers of the battery under study are metal ions,and there are few researches on such carriers as H+.As shown in Fig.3b,H+has a radius nearly one-fifth that of Li+,and their hydration radius is also much smaller than that of commonly used metal charge carriers.Therefore,the insertion of H+in the electrode material will result in a small lattice strain due to its smaller ionic radius [114].In addition to smaller lattice strains,H+have another advantage as carriers.In aqueous electrolyte,there are a large number of hydrogen bonds between water molecules,which form a network of hydrogen bonds.As H+migrates,the H+displaces hydrogen ions from one of the water molecules in the hydrogen bond network,and the displaced hydrogen ions then displace hydrogen ions from the next water molecule through the hydrogen bond,thus transferring the H+to the electrode surface by continuous displacement.This mode of H+transport is called Grotthuss proton conduction,and it is this mode of transport that gives the H+faster migration dynamics.

    For PBAs,it will produce a lot of crystal water in the crystal lattice when using the co-precipitation method.In addition,due to some defects in PBAs,these defects will expose M atoms in PBAs outside the unit cell and combine with water molecules.Then,these water molecules bonded with M atoms can communicate crystal water in each lattice and form a hydrogen bond network in PBAs.Therefore,H+can also be transported by Grotthuss proton conduction in PBAs to achieve faster H+transport.Combined with the small ion radius of H+and the fast dynamics of Grotthuss proton conduction,PBAs as electrode for aqueous H+batteries will produce excellent cyclic stability and rate performance.For the application of PBAs in aqueous H+batteries,Ji et al.first use the Cu[Fe(CN)6]0.63·□0.37·3.4H2O as cathode and explore its electrochemical performance in 2.0 M H2SO4electrolyte [81].Impressively,the Cu[Fe(CN)6]0.63·□0.37·3.4H2O exhibits a capacity of 48 mAh g-1at 5000 C that is the half capacity in 1 C.In addition,it shows an unprecedented long cycle performance that can retain 60% of its initial capacity at 500 C after 730,000 cycles.Later,they find that the Ni[Fe(CN)6]0.60·□0.40·4.5H2O also exhibits excellent cycle and rate performance when it used as electrode for aqueous H+batteries [115].Furthermore,they find that Grotthuss proton conduction can still occur even at low temperature due to the unique chemical environment of crystal water in PBAs.At -40 °C,the Ni[Fe(CN)6]0.60·□0.40·4.5H2O still possesses well rate performance that can deliver~50% of theoretical capacity at 50 C.Therefore,the application of PBAs in aqueous H+batteries at low temperature also has a broad development prospect.

    2.1.5 PBAs for Aqueous NH4+ Batteries

    Apart from H+,NH4+as a non-metallic charge carrier also has received more and more attention in recent years.In 2011,Cui and co-workers explored the insertion performance of NH4+,Na+and K+in K0.9Cu1.3Fe(CN)6and K0.6Ni1.2Fe(CN)6[74].Due to the large ionic radius of NH4+,it exhibits a higher insertion potential.Furthermore,both of the two PBAs exhibit excellent cycle stability during the insertion of NH4+,which opens the door of the PBAs in aqueous NH4+batteries.The first rocking-chair NH4+battery was proposed by Ji and co-workers [79].For this battery,it uses (NH4)1.47Ni[Fe(CN)6]0.88as cathode and 3,4,9,10-perylenetetracarboxylic diimide (PTCDI) as anode.Due to the good structure stability of (NH4)1.47Ni[Fe(CN)6]0.88,it exhibits an initial capacity of 51.3 mAh g-1and retains 74% at 5 C after 2000 cycles.However,due to the limited stability of the anode,the cycling stability of this full cell is worse that can retain 67% at 3 C after 1000 cycles.Later,Huanget.al.proposed a novel aqueous Zn2+/NH4+hybrid battery that uses the K0.72Cu[Fe(CN)6]0.78·3.7H2O as cathode,zinc as anode in 1 M (NH4)2SO4/0.1 M ZnSO4electrolyte [116].It is noteworthy that the insertion/extraction of NH4+and deposition/plating of Zn2+occurs in cathode and anode,respectively,during the battery working.The K0.72Cu[Fe(CN)6]0.78·3.7H2O shows an excellent cyclic stability that can retain 78% at the current density of 1800 mA g-1after 1000 cycles.Impressively,due to the stability of the zinc anode in the electrolyte,the full cell shows a better cycle stability that can retain 76.5% at the current density of 1800 mA g-1after 1000 cycles.Thus,the zinc anode is worthy of consideration in future studies of PBAs in aqueous NH4+batteries.Aside from the CuFe-PBAs and NiFe-PBAs,the FeFe-PB (Fe[Fe(CN)6]0.88·□0.12·2.8H2O) in aqueous NH4+battery also has well electrochemical performance [117].Impressively,the Fe[Fe(CN)6]0.88·□0.12·2.8H2O exhibits a nearly zero strain characteristic when it hosts NH4+.Thus,its cyclic retention rate can reach up to 78%after 50,000 cycles.By increasing the agitation rate during the preparation of PBAs,Huang and co-workers synthesized NaFeIIIFeII(CN)6with ball-cutting morphology [118].Compared with the cubic morphology NaFeIIIFeII(CN)6,the ball-cutting morphology NaFeIIIFeII(CN)6improves the uniformity and reduces the particle size of NaFeIIIFeII(CN)6.So,it possesses an excellent cycle performance with no fading after 50,000 cycles.In conclusion,different types of PBAs both have good electrochemical performance in hosting NH4+and waiting to be further explored.

    2.2 PBAs for Aqueous Multivalent-Ion Batteries

    For the storage of multivalent ions,it mainly presents three challenges.Firstly,multivalent ion possesses strong hydration in the process of ion diffusion [119].After hydration,the hydrated ionic radius of multivalent ion is going to be particularly large and the specific change value is shown in Fig.3b.Secondly,the insertion of multivalent ion in electrode needs more electrons to retain electric neutrality,which means more valence changes will occur.Generally speaking,the large sudden change of valence will cause large local deformation to the crystal structure of host materials [120-122].Thirdly,the multivalent ion will be subject to stronger electrostatic repulsion in host materials that will result in the instability of crystal structure.Therefore,the selection of electrode materials for aqueous multivalent-ion batteries becomes particularly important.

    For PBAs electrode materials,it has a large ion insertion site which can host the multivalent ion with large hydrated ionic radius.Besides,the vacancy with large volume in PBAs structure also can be used as insertion site for multivalent ion in aqueous electrolyte.As shown in Fig.2,in the crystal structure of PBAs,the vacancy site often absorbs some coordinated water molecules.Therefore,compared with other insertion site,the vacancy site will allow the host ions adjacent to more coordinated water molecules within the structure,thus forming a hydration shell.This hydration shell can provide more electrostatic shielding to host ions and thus improve the thermodynamic stability of host ions at vacancy sites.It is this mechanism that allows vacancies to be used as insertion sites for multivalent ions.Meanwhile,the PBAs possess multiple redox sites,which make them able to accept more valence changes.Combined with these advantages,PBAs will be a highly potential electrode material for aqueous multivalent-ion batteries.Next,these sections will detail the applications of different PBAs in various aqueous multivalent-ion batteries.

    2.2.1 PBAs for Aqueous Mg2+ Batteries

    Magnesium is one of the largest reserves of metals in nature,which makes Mg2+batteries have a greater cost advantage over Li+batteries.In addition,aqueous Mg2+batteries also show great advantages in terms of safety and environmental friendliness.Due to its excellent performance in Na+and K+batteries,PBAs have been promoted to be applied in aqueous Mg2+batteries in recent years.In 2013,Cui et al.systematically explored the electrochemical behavior of Mg2+storage in K0.6Ni1.2Fe(CN)6cathode and proposed a novel solution to decrease capacity decay [75].For the insertion of Mg2+in K0.6Ni1.2Fe(CN)6,they found that the partial dissolution of K0.6Ni1.2Fe(CN)6in the electrolyte caused significant decay process at the beginning of the cycle and then stabilization.By adding a small amount of Ni2+to the electrolyte,not only the capacity is recovered but also the cycling stability is improved.The principle of this method can be attributed to two points:Firstly,in the process of discharge,Ni2+in the electrolyte is inserted into PBAs to improve the structural stability of PBAs;secondly,electrolyte with a certain concentration of Ni2+can inhibit the dissolution of high-spin Ni atom.Obviously,this simple method of improving PBAs cycling performance by changing electrolyte composition deserves more attention.For CuFe-PBAs,Yamada et al.reported the reversible insertion of Mg2+in K0.1Cu[Fe(CN)6]0.7·3.6H2O [123].Impressively,through ex situ M?ssbauer spectroscopy,they reveal that both Fe3+/Fe2+and CuII/CuIexperience solid-state redox when Mg2+inserts into K0.1Cu[Fe(CN)6]0.7·3.6H2O.About the PBAs in aqueous full Mg2+batteries,Xia and co-workers recently used Na1.4Ni1.3Fe(CN)6·5H2O cathode and polyimide anode to construct a full battery [124].Figure 10a,b shows the electrochemical redox mechanism (insertion and extraction)of polyimide with Li+,Na+and Mg2+.For this battery,it exhibits a high output voltage of 1.5 V and an excellent cycle life that can retain about 60% of its initial capacity after 5000 cycles,which is currently the best cycling performance of aqueous Mg2+battery.Obviously,it has a strong guiding value for the future study for the PBAs used in aqueous full Mg2+batteries.

    2.2.2 PBAs for Aqueous Zn2+ Batteries

    About the PBAs in aqueous Zn2+batteries,Liu et al.proved that K2Zn3[Fe(CN)6]2is more suitable to be used in aqueous Zn2+batteries than aqueous Na+and K+batteries because K2Zn3[Fe(CN)6]2has better structural stability in aqueous Zn2+electrolyte [22].Besides,K2Zn3[Fe(CN)6]2-based full battery was successfully assembled the Zn anode,achieving a high output voltage of 1.7 V and an energy density of 100 Wh kg-1.Later,they explored the effect of Zn3[Fe(CN)6]2with different morphologies for Zn2+insertion [17].Figure 10c shows the coordination environments for Zn and Fe atoms in two different crystal systems (cubic and rhombohedral),and corresponding SEM images of three Zn3[Fe(CN)6]2with different morphologies are shown in Fig.10d.And the experiment results show that,compared with truncated octahedral and octahedral particles,cube octahedral particles have more advantages in rate and cycling performance because the cube octahedral structure surface orientations is more beneficial to Zn2+diffusion and Zn3[Fe(CN)6]2structure stability.

    About the CuFe-PBAs in aqueous Zn2+battery,Wang and co-workers demonstrated for the first time that Zn2+can reversibly insert/extract into/from KCuFe(CN)6nanocube (Fig.10e) in 1 M ZnSO4electrolyte [125].However,its cycling stability is poor with 77% capacity retention after 20 cycles.For comparison,Mantia and co-workers reported that KCuFe(CN)6(K0.71Cu[Fe(CN)6]0.72·3.7H2O)can provide 90% of theoretical capacity with a capacity retention of 96.3% after 100 cycles at 1 C in a dilute electrolyte of 20 mM ZnSO4solution [126].According to the conventional view,the aging performance may be attributed to the destruction of crystal structure of PBAs upon Zn2+insertion/extraction.However,Mantia and co-workers failed to find Cu and Fe elements in the electrolyte through ICP-MS analysis [127].For this phenomenon,Lim and co-workers claimed that the inserted Zn2+not only existed in the interstitial lattice vacancy but also appeared in the pre-occupied lattice site during the cycling process [128].Zn atoms in the crystal structure will replace Cu atoms,leading to the formation of Cu(CN)2and CuZn(CN)4.Since Cu(CN)2and CuZn(CN)4do not have redox couple,they cannot provide capacity.Thus,this irreversible phase transition leads to the capacity fading.In particular,the lower the concentration of Zn2+in the solution,the more difficult this irreversible phase transition is to occur,which explains why KCuFe(CN)6has better performance in 20 mM ZnSO4electrolyte compared with 1 M ZnSO4electrolyte.In order to overcome this irreversible phase change,Mantia fabricated a series of copper zinc hexacyanoferrates(CuZnHCF) mixtures with different Cu/Zn ratios [129].It is found that the irreversible phase transition can be effectively avoided for CuZnHCF with Cu/Zn ratio of 93:7 (SEM image shown in Fig.10f) and the capacity fading problem can be solved to a large extent.Compared with the CuHCF without Zn element,it exhibits better capacity retention of 85.5% after 1000 cycles in 20 mM ZnSO4solution.Obviously,this method provides a new direction for optimizing PBAs in aqueous multivalention batteries.About the application of double-atom redox PBAs in aqueous Zn2+batteries,Zhi and co-workers synthesized KCoFe(CN)6and used it as cathode of aqueous Zn2+batteries [77].With two redox sites,this battery exhibits a high capacity of 173.4 mAh g-1at current density of 0.3 A g-1.Impressively,they adopt a sol-gel transition strategy constructing a solid-state cable-type battery,and Fig.10g schematically illustrates its configuration.Figure 10h,i shows the practical application of the solid-state cable-type battery in two different electronic devices,electronic hygrometer and LEDs,respectively.Obviously,combined with the lower electrode potential of zinc anode,the cable-type battery based on double-atom redox PBAs exhibits a huge potential in flexible energy storage that provides a new direction for the practical application of aqueous battery based on PBAs.

    2.2.3 PBAs for Aqueous Ca2+ Batteries

    Since hydrated Ca2+binds fewer water molecules on the surface [130],the investigations of aqueous Ca2+batteries have great development potential.About the PBAs used in aqueous Ca2+batteries,Yao and co-workers synthesized K0.02Cu[Fe(CN)6]0.66·3.7H2O as cathode material for Ca2+storage [131].Due to the low crystal structure strain (1.1%),the K0.02Cu[Fe(CN)6]0.66·3.7H2O exhibits an excellent long cycle performance that can retain 88% of its initial capacity at 5 C after 2000 cycles.In addition,the K0.02Cu[Fe(CN)6]0.66·3.7H2O also exhibits excellent rate performance in aqueous Ca2+batteries.Its capacity can be retained to 65% when the current density increases from 0.2 C to 20 C,while that of in aqueous Mg2+battery can only be maintained to 33%.This difference in rate performance can be summarized in two points:First,Ca2+with a smaller charge density will encounter less resistance when it is inserted/extracted into/from the PBAs.Second,hydrated Ca2+binds fewer water molecules on the surface,speeding up the transport of Ca2+in PBAs.Similarly,the K0.02Cu[Fe(CN)6]0.64·9H2O also demonstrated same excellent Ca2+storage performance in aqueous Ca2+battery,which can retain 94% of its initial capacity at 5 C after 1000 cycles [132].

    About the electrochemical performance of PBAs in aqueous Ca2+battery,the concentration of electrolyte also has influence.Jeong and co-workers studied the electrochemical properties of CuFe(CN)6in the different concentrations of electrolyte and found that in 1.0 mol dm-3Ca(NO3)2,the CuFe(CN)6delivers a 60 mAh g-1initial capacity and retains 50% at 10 C after 1000 cycles [133].But in 8.4 dm-3Ca(NO3)2,it exhibits a better storage performance that can deliver a 70 mAh g-1initial capacity and retains 95% at 10 C after 1000 cycles because the hydration number of Ca2+in the high-concentration electrolyte was smaller.Therefore,it is a simple and effective method to improve the electrochemical performance of PBAs in aqueous Ca2+batteries by changing the electrolyte concentration,which is worthy to be applied in future studies on aqueous Ca2+batteries.

    2.2.4 PBAs for Aqueous Al3+ Batteries

    For the insertion of Al3+in CuFe-PBAs,Cui et al.use K0.03Cu[Fe(CN)6]0.65·2.6H2O as host material to reversibly store Al3+[76].Figure 11a presents the lattice structure with all occupied atoms and possible insertion sites for hosting ions in this K0.03Cu[Fe(CN)6]0.65·2.6H2O.Furthermore,they confirmed the shielding effect of crystal water and hydration shell of Al3+on electrostatic repulsion between Al3+and K0.03Cu[Fe(CN)6]0.65·2.6H2O,which lays a foundation for the insertion of Al3+in PBAs.Later,Gao et al.used Al2(SO4)3solution as aqueous electrolyte to explore the Al3+reversible insertion characteristic in nanosized KCu[Fe(CN)6]·8H2O (Fig.11b) [10].As reported in Fig.11c,typical redox couples in CV curve are well indexed to the schematic insertion/extraction positions of Al3+in PBAs framework.For Al3+uptake,KCu[Fe(CN)6]·8H2O can release a discharge capacity of 46.9 mAh g-1at 400 mA g-1with retention of 54.9% after 1000 cycles.In addition to reducing electrostatic repulsion,PBAs with small particle size can also improve the migration rate of Al3+by reducing the diffusion distance.This is proved by Wills and Pang et al.via the investigation on CuFe(CN)6(Fig.11e) and K2CoFe(CN)6(Fig.11d,f),respectively,and got well performance [134,135].Although these two kinds of PBAs have wonderful Al3+insertion ability,the reason of capacity fading is not clear.Later,Wu and co-workers reported the working mechanism in PBAs-based Al3+batteries [136].They clarified that the reversible insertion/extraction of Al3+in K0.02Ni1.45[Fe(CN)6]·2.6H2O is related to a single-phase reaction without producing polycrystalline structure.In addition,they also found that the capacity fading is not mainly caused by the collapse of PBAs structure,but the transfer of Ni in K0.02Ni1.45[Fe(CN)6]·2.6H2O to Al foil anode,resulting in interface instability between anode and electrolyte.Obviously,this study provides a new idea for us to solve capacity fading problem of PBAs upon multivalent ions uptake.

    Fig.11 a Structure of K0.03Cu[Fe(CN)6]0.65·2.6H2O and corresponding insertion sites for trivalent ions.Reproduced with permission from Ref.[76].Copyright 2015,Wiley-VCH.b SEM and TEM images of as-prepared KCu[Fe(CN)6]·8H2O powders.c Typical CV curve of KCu[Fe(CN)6]·8H2O in Al2(SO4)3 solution and the schematic positions of Al3+ in the framework.Reproduced with permission from Ref.[10].Copyright 2015,The Royal Society of Chemistry.d TEM image of K2CoFe(CN)6.Reproduced with permission from Ref.[135].Copyright 2020,Elsevier.e Working process of aqueous CuFe(CN)6/TiO2 battery for Al3+ storage.Reproduced with permission from Ref.[134].Copyright 2018,Springer.f Crystal structure of K2CoFe(CN)6.Reproduced with permission from Ref.[135].Copyright 2020,Elsevier

    2.2.5 PBAs for Other Aqueous Multivalent Metal-Ion Batteries

    Besides these above multivalent metal ions,the insertion performance of Ba2+,Sr2+,Fe2+,Co2+,Ni2+,Cu2+,Pb2+,Y3+,La3+,Ce3+,Nd3+and Sm3+in PBAs also has been explored.For Ba2+and Sr2+,Cui and co-workers explored their insertion performance in NiFe(CN)6[75].When NiFe(CN)6inserting Sr2+,it exhibits an initial capacity of 40 mAh g-1and retains 64% at 5 C after 2000 cycles.Because Ba2+has smaller hydrated ionic radius,the insertion performance of NiFe(CN)6in aqueous Ba2+battery is better.It exhibits a 42 mAh g-1initial capacity and retains 93% at 5 C after 2000 cycles.About these two ions as host ion in PBAs,the studies are few.Thus,future studies can put into these aspects.Later,Cui et al.continue to explore the insertion performance of Co2+,Ni2+,Cu2+,Pb2+,Y3+,La3+,Ce3+,Nd3+and Sm3+in K0.03Cu[Fe(CN)6]0.65·2.6H2O [76].They find that K0.03Cu[Fe(CN)6]0.65·2.6H2O has a serious decay process when inserting Co2+,Ni2+and Cu2+resulting in poor cycling performance.The main reason for this phenomenon is that partial replacement of these ions with the highspin transition metal atom in K0.03Cu[Fe(CN)6]0.65·2.6H2O results in irreversible phase transition.Therefore,the main problem for the insertion of these ions into PBAs is to resolve this irreversible phase transition.It should be noted that the insertion and extraction of Fe2+in FeFe-PB do not produce this irreversible structural phase change.In 2019,Ji et al.used Fe[Fe(CN)6]0.73·□0.27·3.6H2O as cathode,Fe as anode and FeSO4solution as electrolyte to build a full cell[137].For this cell,it exhibits an excellent cycle stability that can retain 88% of its initial capacity at 1 C after 1000 cycles.

    In addition to these common metal ions,Cui et al.found that some lanthanide metal ions can also be reversibly inserted/extracted in/from PBAs such as La3+,Ce3+,Nd3+and Sm3+.For these ions inserted into K0.03Cu[Fe(CN)6]0.65·2.6H2O,they all show similar voltammogram,which means that they have the same insertion sites in K0.03Cu[Fe(CN)6]0.65·2.6H2O.Impressively,they both can deliver a capacity of approximately 60 mAh g-1at 0.2 C.Although these lanthanide metal ions can be reversibly inserted/extracted in/from PBAs,they are rarely studied due to their low natural reserves and high cost.It is worth noting that the insertion of Pb2+and Y3+possesses the fastest ion migration kinetics and most reversible electrochemistry among these uncommon divalent and trivalent metal ions,respectively [76].For the insertion of Pb2+in K0.03Cu[Fe(CN)6]0.65·2.6H2O,it exhibits an initial capacity of 55 mAh g-1and retains 73% at 5 C after 2000 cycles.When the K0.03Cu[Fe(CN)6]0.65·2.6H2O inserting Y3+,it exhibits an initial capacity of 40 mAh g-1and retains 87%at 5 C after 2000 cycles.Obviously,the insertion of these two ions in PBAs possesses great exploration value and is waiting to be studied.

    2.3 Relationship between Performance and Structure

    From the above introduction,it can be clearly found that for PBAs,their electrochemical performance (redox potential,capacity,cycling stability and rate performance) is closely related to its structure (metal atom in high-spin state,defect,particle size,morphology,etc.) and these electrochemical performances of various PBAs in aqueous battery are summarized in Table 1.From Table 1,it can be concluded that different types of metal in high-spin state lead to different redox potential and voltage platforms for PBAs.When Fe,Co and Mn are in high-spin state,this type of PBAs has two voltage platforms because both these metals and Fe in low spin have electrochemical activity.Conversely,when Ni,Cu and Zn are in high-spin state,this type of PBAs has only one voltage platform because the Ni,Cu and Zn are electrochemically inert.It is worth noting that Cu is sometimes electrochemically active in some acidic electrolytes and produces voltage platform.Meanwhile,in terms of redox potential,FeFe-PB,CoFe-PBAs and MnFe-PBAs are generally higher than NiFe-PBAs,CuFe-PBAs and ZnFe-PBAs.Although the redox potential of different PBAs is different in aqueous batteries,most of them are within the decomposition voltage of water (Fig.12).A small part of PBAs whose charging potential exceeds the oxygen evolution potential can be still used as electrode in aqueous battery by increasing electrolyte concentration or adjusting pH to inhibit the evolution of oxygen.

    Fig.12 Voltage working window of PBAs for various aqueous batteries

    Table 1 Electrochemical performance of representative PBAs in various aqueous batteries

    Table 1(continued)

    Table 1(continued)

    In addition to the redox potential,the metal atom in highspin state also has great influence on capacity and cycling stability.In order to more intuitively express the capacity and cyclic stability of different PBAs in aqueous battery,we summarize them in Fig.13.Obviously,compared with NiFe-PBAs,CuFe-PBAs and ZnFe-PBAs,FeFe-PB,CoFe-PBAs and MnFe-PBAs can provide more capacity because both metals in high-and low-spin state can occur redox reaction.But in terms of cyclic stability,NiFe-PBAs,CuFe-PBAs and ZnFe-PBAs are superior.In addition,defects in the PBAs structure will lead to the reduction of active sites in PBAs,resulting in the decrease of capacity.At the same time,the existence of defects leads to the increase of crystal water in PBAs and thus reduces the cyclic stability.Therefore,Fig.14 shows that the same-type PBAs have different crystal water content,resulting in different capacities and cycling performance.For the rate performance,high structural stability is beneficial to the rapid insertion and extraction of ions under high current density.Therefore,Table 1 shows that the capacity loss of NiFe-PBAs,CuFe-PBAs and ZnFe-PBAs is less when the current density increases.In addition,the morphology and particle size of PBAs crystals also have a great influence on the rate performance.Generally,PBAs with more uniform morphology and smaller particle size have better rate performance because they have more sufficient contact with the electrolyte and thus are more conducive to rapid ion transport.

    Fig.13 Cycling performance of different PBAs in various aqueous batteries

    Fig.14 a Principle of Na2Ni[Fe(CN)6]/Na2Fe[Fe(CN)6] rocking-chair desalination battery.b SEM image of Na2Ni[Fe(CN)6] particles.c SEM image of Na2Fe[Fe(CN)6] particles.d Galvanostatic cycling performance of the Na2Ni[Fe(CN)6]/Na2Fe[Fe(CN)6] cell in seawater.Reproduced with permission from Ref.[144].Copyright 2017,American Chemical Society.e Schematic of zinc/ferricyanide hybrid desalination flow battery operation during discharge (desalination) and charge (salination).Reproduced with permission from Ref.[145].Copyright 2018,American Chemical Society.f Scheme showing the operation of the Cu3[Fe(CN)6]2·nH2O/Bi desalination battery.Reproduced with permission from Ref.[78].Copyright 2019,American Chemical Society

    Through these relationships between PBAs performance and structure,we can more intuitively select the appropriate PBAs in each aqueous battery.As the capacity provided by NiFe-PBAs,CuFe-PBAs and ZnFe-PBAs in aqueous monovalent metal-ion batteries is limited,after improving the structural stability of FeFe-PB,CoFe-PBAs and MnFe-PBAs,they are more suitable for application in aqueous monovalent metal-ion batteries.For aqueous H+batteries,the high structural stability NiFe-PBAs are a better choice,because the rapid transfer of H+requires high structural stability for host materials.In aqueous NH4+batteries,due to the excellent structural compatibility of NH4+and FeFe-PB,the application of FeFe-PB has more prominent.For aqueous multivalent-ion batteries,NiFe-PBAs and CuFe-PBAs with higher structural stability are more suitable because of the large hydration ion radius of multivalent ions.Impressively,due to the excellent performance of CoFe-PBAs in aqueous Al3+batteries and Zn2+batteries recently,it is also worth trying in other aqueous multivalent-ion batteries.

    3 PBAs in Desalination Batteries

    To solve the shortage of freshwater source,various physical and electrochemical desalination ways have been studied in the past decades [138-140].Among physical desalination method,the reverse osmosis,also known as ultrafiltration,is the most representative.The reverse osmosis method is to separate freshwater from seawater by using a semi-permeable membrane which allows only solvent but not solute to permeate.The desalination process can be divided into two steps.The first step is to use a semi-permeable membrane to separate seawater and freshwater on both sides of the membrane,and then,the freshwater moves to the seawater side under the action of osmotic pressure.The second step is to use applied pressure to reverse osmosis the freshwater in the seawater side into the freshwater.Obviously,this physical desalination method has two major disadvantages:the use of high-pressure equipment and regular cleaning and replacement of semi-permeable membrane.In contrast,for electrochemical desalination,such as desalination battery,the desalination is done by the impurity ions hosted in the electrodes,so it does not require an external device.Even in some desalination batteries,semi-permeable membranes may not be used.It is worth noting that the semi-permeable membrane used in desalination battery can be mainly divided into anion exchange membrane(AEM) and cation exchange membrane (CEM).AEM is mainly composed of ion exchange resin and basic active group,and CEM is mainly composed of ion exchange resin and acid active group.Obviously,avoiding the use of semi-permeable membranes will greatly reduce the cost of desalination.Of course,because desalination battery needs to work under the condition of applied voltage,its energy consumption is higher than reverse osmosis.But in terms of simple equipment and environmentally friendly,desalination battery still has great potential.For desalination battery,it consists of cation and anion storage electrodes which remove salt ions from the water through a reversible electrochemical reaction.At present,the electrode materials commonly used in desalination batteries include metal oxide,polyanion polymer,PBAs,etc.Compared with other electrode materials,PBAs mainly have the following advantages as a desalination battery electrode material.Firstly,due to its large ion channel and insertion site,the PBAs can provide higher desalination rate and requires less energy consumption.Secondly,PBAs as desalination battery electrode will greatly reduce the cost due to the simple synthesis method of PBAs.Thirdly,PBAs have better compatibility with multivalent ions in seawater because the crystal water in PBAs structure has a certain shielding effect on the electrostatic repulsion of multivalent ions.In addition to these advantages,PBAs also have a structural disadvantage that needs to be noted.Since most seawater is alkaline,this will reduce the stability of PBAs in seawater and thus affect the desalination performance.However,such structural advantage can be well controlled by controlling PBAs vacancies and crystallinity,so the PBAs-based desalination batteries have received a large number of concerns in recent years.

    In desalination battery,the mass transfer resistance of ions in water inevitably leads to the degradation of battery performance,especially when the water concentration is very low.Different from high-concentration electrolytes that are usually defined as concentration higher than 0.1 M,the brackish water (about 500-2000 mg of solute per liter of solution)possesses a small number of accessible ions in the interface of electrode and electrolyte.Therefore,the mass transfer rate of ions in this kind of desalination batteries is slow,and the desalination process can only be guaranteed under the low current density,which is usually about 10 A m-2[141,142].In order to introduce the applications of PBAs in desalination batteries more comprehensively,herein the desalination batteries are divided into two groups:the desalination batteries for highly concentrated brine and the desalination batteries for brackish water.The recent advances of PBAs in these two desalination batteries are discussed,respectively,in the following sections.

    3.1 PBAs for Highly Concentrated Brine Desalination

    For traditional desalination battery like Ag/Na0.5MnO2battery,Na+and Cl-are captured in electrode materials through a complex four-step charge and discharge progress and the removal rate of Na+can reach 25% [143].During the process of discharging (step 1),Na+and Cl-are inserted into the host material by insertion reaction.The deionized water is then exchanged with the water to be treated (step 2).After the above steps,the battery in step 2 is charged so that the ions inserted in discharge reaction can be extracted from the host material (step 3).Then,exchanging the concentrated water with the water to be treated,and the battery continues to perform the discharging process (step 4).Obviously,this desalination process is complex.In order to simplify the desalination process and promote the removal rate of Na+,Yoon and co-workers used PBAs to design a simpler and more efficient desalination battery [144].Figure 14a represents the schematic diagram of PBAs-based desalination battery.The electrode materials used in this desalination battery are highly crystalline Na2Ni[Fe(CN)6] (NaNiHCF)and Na2Fe[Fe(CN)6] (NaFeHCF) for better cation capture.The morphology of NaNiHCF and NaFeHCF is revealed in Fig.14b,c,respectively.In this desalination battery,the water to be treated is divided into two parts by AEM.In fact,the motion of ions during the progress of charging and discharging shows that the desalination battery works in the same way as aqueous rocking-chair battery.Impressively,the use of AEM creates a charge difference between two different solutions of desalination battery,which leads to the concentration and dilution of solutions during charging and discharging without the need to replace the treated water during desalination process,thus simplifying the desalination process.It is noteworthy that NaNiHCF/NaFeHCF full cell demonstrates a good galvanostatic cycling performance with an initial specific capacity about 56 mAh g-1at 0.1 A m-2(Fig.14d).Meanwhile,this desalination battery possesses an excellent Na+removal efficiency up to 40%.For further improving the desalination performance and promoting its practical application,shortly afterward,Rivest and his co-workers presented a hybrid desalination flow battery that consists of a zinc anode and a PBA cathode [145].Impressively,this battery possesses a high operation cell voltage of 1.25 V,showing high round-trip efficiency and electrical storage capacity.Figure 14e represents the schematic diagram of hybrid desalination flow battery.Different from above reports of PBAs in desalination batteries,the cathode of this battery is mixed PBAs in liquid form,K4Fe(CN)6/K3Fe(CN)6. Na+removal depends on the reaction between K4Fe(CN)6/K3Fe(CN)6.Since the liquid-form PBAs possess larger surface,it is more favorable for hosting ions.Thus,impressive salt removal efficiency can be achieved at 85%that is at least two times higher than most desalination batteries reported till now.

    As is well known,the membranes play an important role in existing desalination technologies.For reverse osmosis,it requires a semi-permeable membrane for only transporting water molecules but blocking other ions to ensure that the desalination process takes place [146-148].For conventional desalination batteries,they also need membranes to simplify the desalination steps or improve the performance.Actually,the use of membrane is associated with membrane fouling and replacement cost.Thus,eliminating the need for membrane in highly concentration brine desalination batteries can greatly reduce the cost of desalinations and realize large-scale application of desalination batteries.Since the principle of desalination battery is performed by specific electrode reaction (such as Na+insertion in cathode and Cl-insertion in anode),it possesses the potential to avoid using semi-permeable membrane.Soon after,membranefree Na3Ti2(PO4)3/Bi desalination battery was put forward[149].Since Cl-removal from BiOCl requires a lower overpotential in acidic conditions than in neutral solutions,the anode reaction is more suitable to perform in acidic conditions.But under acidic conditions,the self-discharge phenomenon of Na3Ti2(PO4)3is serious and it will lead to the spontaneous oxidation of Na3Ti2(PO4)3,so it is not conducive to the desalination process of Na3Ti2(PO4)3[150].Thus,a membrane is still needed to separate the cathode and anode to allow them to work under different acidic conditions.As is well known,CuFe-PBA was demonstrated with excellent cyclability for charge and discharge reactions in acidic media.Therefore,replacing Na3Ti2(PO4)3with CuFe-PBA can achieve membrane-free desalination battery.On the basis of this view,Choi and co-workers constructed a novel desalination battery that uses Cu3[Fe(CN)6]2·nH2O as cathode and Bi as anode to enable membrane-free desalination,and Fig.14f represents the schematic diagram of this rocking-chair desalination battery [78].This work further promotes the application of PBAs in membrane-free desalination batteries.

    3.2 PBAs for Brackish Water Desalination

    In past studies,researchers have tried various ways to solve the problem of large mass transfer resistance of ions when PBAs-based desalination batteries were used in brackish water.One of the easiest is to control particle size of PBAs[151-157].Co-precipitation method is a widely accepted route to prepare PBAs.The reaction rate can be slowed down by coordination agents to get less defective PBAs with small size.Based on this view,Smith and co-workers prepared small-sized Na2NiFe(CN)6via co-precipitation [158].Since Na-rich NiFe-PBA does not require a pre-solinization process,they investigated a NaNiFe(CN)6/Na2NiFe(CN)6battery to treat brackish water as represented in Fig.15a.Figure 15b shows the salt adsorption capacity of NaNiFe(CN)6.For this desalination battery,it only exhibits a 3% sodium removal rate,because the system simply uses granularity control.Although this desalination battery is capable of treating brackish water,its desalination capacity is too poor to meet practical requirements.Later,Logan and co-workers investigated a NaCu[Fe(CN)6]/NaCu[Fe(CN)6] symmetrical cell to optimize the desalination performance in brackish water [159].For this desalination battery,it needs to operate at a low current density due to the influence of mass-transfer limitation between electrode and electrolyte.To increase desalination capability,different numbers of ion exchange membranes are added between electrodes to form a hybrid system.Thus,the removal rate of sodium ions is increased to 25%,which demonstrates the feasibility of this system in brackish water.Although these above two desalination batteries are capable of handling brackish water,the fundamental solution is to solve the problem of mass transfer resistance in electrolytes [160].Then,Yoon and co-workers proposed a multi-channel desalination battery (MC-DB) to solve this problem [161].Figure 15c exhibits the operation process based on NaNiFe(CN)6/Ag desalination battery.Compared with conventional desalination battery with only one feed channel,the MC-DB has two side feed channels and one middle feed channel.The two side feed channels are used for introducing a highly concentrated solution to electrode and one middle feed channel for water desalination.The most innovative part of the battery is that it can adjust the concentration of water around the PBA cathode and Ag anode to reduce the ion migration resistance.Impressively,the MC-DB system shows a desalination capacity of approximately 53 mg g-1.Particularly,the maximum salt removal rate can reach up to 0.0576 mg g-1s-1when it treats 10 mM NaCl solution.Therefore,the future application of PBAs-based desalination battery for treating brackish water can be further explored with this battery as a model.

    Fig.15 a Schematic of NaNiFe(CN)6/Na2NiFe(CN)6 desalination battery operation during discharge (desalination) and charge (salination).b Salt adsorption capacity of NaNiFe(CN)6.Reproduced with permission from Ref.[158].Copyright 2017,Elsevier.c Schematic of multi-channel NaNiFe(CN)6/Ag desalination battery operation during discharge (desalination) and charge (salination).Reproduced with permission from Ref.[161].Copyright 2019,American Chemical Society

    Table 2 summarizes the performance metrics of recently advanced PBAs as desalination battery electrode.As we can see,NaNiHCF/NaFeHCF desalination battery possesses excellent cycle performance with 100 reversible cycles.For desalination capability,the Zn/K3[Fe(CN)6] cell has the highest specific capacity with 81 and 820 mAh g-1(approximately 86% salt removal) based on the mass of cathode and anode,respectively.Therefore,PBAs-based desalination batteries inherit the potential to be explored.However,only few PBAs can be used in desalination batteries due to poor structural stability of PBAs,such as ZnFe-PBA and MnFe-PBA.Thus,future research can be focused on improving the stability of PBAs to promote more applications of PBAs in desalination batteries.More importantly,we should recognize that this electrochemical desalination technology is a bit too simple for large-scale application.Therefore,in the future exploration of electrochemical seawater desalination,we can develop a hybrid technology that can combine multiple desalination systems to achieve better results.

    Table 2 Summary of performance metrics of typical desalination batteries

    4 Conclusions and Outlook

    In this review,the recent progress of PBAs as electrode for aqueous monovalent-ion and multivalent-ion batteries is discussed.Meanwhile,the application of PBAs in aqueous desalination batteries is also introduced.PBAs have received a lot of attention for electrode materials in aqueous batteries due to their 3D rigid framework and facile synthesis.Among various host ions,Na+with small hydrated ionic radius and large energy density makes it more suitable for inserting in PBAs lattice.Although Li+is a superstar in non-aqueous batteries,it is unsuitable to be hosted in PBAs in aqueous electrolyte due to its large hydration radius,causing structural damage upon inserting/extracting into/from PBAs.However,co-insertion of Li+with other metal ions possesses great research significance due to synergistic effect.For K+,PBAs show excellent rate performance because of the small hydrated ionic radius of K+among monovalent ions.In addition,the application of PBAs in aqueous non-metallic ion battery also exhibits excellent performance.In aqueous H+batteries and NH4+batteries,PBAs exhibit excellent rate performance and cycling stability,respectively,due to the special transport mechanism of H+in aqueous electrolyte and the good structural compatibility among NH4+and PBAs.As for the application of PBAs in aqueous multivalent ions battery,the vacancies of PBAs can be used as insertion site to contain large volume of hydrated multivalent ions,and crystal water in PBAs also has a certain shielding effect on the electrostatic repulsion of multivalent ions in PBAs structure.Therefore,the application of PBAs in aqueous multivalent-ion battery also has application potential.Although PBAs have a good performance in the application of aqueous monovalent battery,it still has great room for improvement in capacity delivery,cycle stability,rate performance,etc.Future research can be considered in the following aspects.(1) The control of vacancies and interstitial water:Although co-precipitation,the most widely used preparation method of PBAs,is simple and efficient in operation,the faster crystallization rate will bring more vacancies and interstitial water.The presence of defects and interstitial water will lead to the reduction of redox sites in PBAs and the stability of PBAs structure,thus resulting in the reduction of capacity and the degradation of cycling performance.Therefore,during the preparation of PBAs,the reaction rate can be reduced by adding chelating agents or adopting multi-step crystallization to improve the crystallinity of PBAs.(2) The suppression of structural change:Generally,PBAs display a transition between one or two crystal types during the insertion and extraction of host ion.Due to the large insertion sites and ion channels of PBAs,such a degree of phase transition will not have a great impact on the structural stability for PBAs.However,for some PBAs such as MnFe-PBAs,the asymmetric distribution of electrons in the d orbitals of Mn3+will lead to the generation of John-Teller distortion,resulting in the phase transformation between three crystal types.This type of structural change has a great impact on the stability of PBAs structure,thus reducing its cyclic stability.In general,this distortion can be well controlled by Mn being partially replaced by Fe to reduce Mn density in PBAs.Therefore,the influence of such multiple-phase reactions on the structural stability of PBAs can be well solved by metal substitution.(3) The enhancement of electronic conductivity:In fact,PBAs as electrode material are poorly conductive.Therefore,it can be combined with various carbon materials or conductive polymers to improve its electronic conductivity.However,these conductive additive agents and binders do not provide capacity in the electrochemical reaction,thus reducing the energy density of PBAs.Actually,improving the conductivity of PBAs on the premise of ensuring the energy density can be solved by directly depositing PBAs onto carbon materials,which are also acted as current collector.This method not only avoids the use of binder and conductive agent,but also achieves the mixing of PBAs and conductive current collector,which has great development potential in improving the conductivity and energy density of the electrode.(iv) The optimization of electrolyte:As shown in Fig.13,the charging voltage of some PBAs has reached the theoretical decomposition voltage of water.Therefore,for the application of PBAs in aqueous batteries,it is very important to widen the electrochemical stability window of electrolyte.At present,the simplest and most effective method is to increase the solute content in the electrolyte and enhance the interaction between anions/cations and water molecules to reduce the content of free water molecules,so as to widen the electrochemical stability window of electrolyte.In addition,as shown in Eq.3,under alkaline conditions,the dissolution of PBAs will be accelerated,especially the dissolution of PBAs containing Co,V,etc.,which will produce toxic components in the electrolyte.Therefore,acidic electrolyte should be selected in general to improve the stability of PBAs and avoid the generation of toxic components.

    For the application of PBAs in aqueous multivalent-ion battery,besides the four points discussed above,there is another aspect that needs special attention.Although the crystal water in PBAs has a certain shielding effect on the electrostatic interaction between multivalent ions and PBAs,the presence of crystal water has negative impact on the structural stability and capacity delivery of PBAs.Therefore,it becomes important to explore the balance between the content of crystal water and the cycling life of PBAs.In addition,most of the current researches on the reaction mechanism of aqueous battery are based on ex situ characterization methods.Therefore,it is also important to use a variety of in situ characterization methods to explore the insertion/extraction mechanism of host ions.

    As desalination battery is a special kind of aqueous battery,the development of PBAs in conventional aqueous battery also promotes its application in desalination battery.Among various PBAs,Zn/K3[Fe(CN)6] shows the best performance,which widens the way for the future research of PBAs in desalination battery.Moreover,PBAs-based multi-channel desalination battery also shows excellent salt removal efficiency in brackish water through special structural design,which provides a clear direction for desalination battery in brackish water in future.In general,PBAs as a desalination battery electrode have the following advantages:(i) For PBAs,they have large ion channels and well electrochemical reversibility,which enable PBAs-based desalination battery to achieve better desalination performance with simplified desalination steps.(ii) PBAs have good chemical stability under acidic conditions,which enables it to form a practical desalination battery with traditional Cl-storage electrode and realize the membrane-free.(iii) Except for large ion channel,PBAs also have small particle size,which makes it suitable to purify brackish water in desalination battery.For the future application of PBAs in desalination batteries,the simplification of desalination steps and the development of membrane-free desalination batteries are the promising directions.Furthermore,in order to get better desalination performance,it is suggested to take the following methods:(i) preparing high-crystallinity,low-defect and small-sized PBAs,(ii) exploring PBAs with new structure,such as MnMn-PBA,CoCo-PBA and CrCr-PBA and their derivatives,and (iii) combining PBAs-based desalination battery with other systems to form hybrid desalination system.

    AcknowledgementsThis work was sponsored by NSAF Joint Fund (U1830106),Science and Technology Innovation 2025 Major Program of Ningbo (2018B10061) and K.C.Wong Magna Fund in Ningbo 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/.

    av电影中文网址| 成年美女黄网站色视频大全免费| 欧美精品亚洲一区二区| 亚洲精品av麻豆狂野| 一区在线观看完整版| av网站免费在线观看视频| 亚洲成人国产一区在线观看| 午夜激情av网站| 99国产精品99久久久久| 午夜福利视频在线观看免费| 精品少妇久久久久久888优播| 韩国精品一区二区三区| 欧美黄色片欧美黄色片| 欧美+亚洲+日韩+国产| 777久久人妻少妇嫩草av网站| 国产精品一区二区在线不卡| 欧美人与性动交α欧美精品济南到| 亚洲av电影在线进入| 新久久久久国产一级毛片| 成年人午夜在线观看视频| 久久久久久久精品精品| 法律面前人人平等表现在哪些方面 | 欧美精品高潮呻吟av久久| 丝袜脚勾引网站| 午夜日韩欧美国产| 青春草亚洲视频在线观看| 最黄视频免费看| 丝袜美腿诱惑在线| 日韩电影二区| 纯流量卡能插随身wifi吗| 国产精品99久久99久久久不卡| 夫妻午夜视频| 十分钟在线观看高清视频www| 一级,二级,三级黄色视频| 国产免费现黄频在线看| 午夜久久久在线观看| 伊人亚洲综合成人网| 欧美xxⅹ黑人| 日韩欧美一区二区三区在线观看 | 久久人人爽人人片av| 在线av久久热| 国产成人精品久久二区二区91| 搡老岳熟女国产| 欧美亚洲日本最大视频资源| 9热在线视频观看99| 亚洲欧美一区二区三区久久| 久久午夜综合久久蜜桃| 亚洲精品乱久久久久久| 别揉我奶头~嗯~啊~动态视频 | 午夜成年电影在线免费观看| 久久久久视频综合| 麻豆国产av国片精品| 欧美成狂野欧美在线观看| 国产精品 国内视频| 国产有黄有色有爽视频| 国产精品九九99| 男女床上黄色一级片免费看| 亚洲成人国产一区在线观看| 国产淫语在线视频| 亚洲少妇的诱惑av| 伊人久久大香线蕉亚洲五| 黄网站色视频无遮挡免费观看| 一本大道久久a久久精品| 午夜福利视频精品| 久久久久久久精品精品| 亚洲精品一区蜜桃| 一边摸一边抽搐一进一出视频| 成人18禁高潮啪啪吃奶动态图| 18禁裸乳无遮挡动漫免费视频| 日本av免费视频播放| 桃红色精品国产亚洲av| 18禁黄网站禁片午夜丰满| 国产97色在线日韩免费| 亚洲国产精品一区三区| 9191精品国产免费久久| 免费在线观看日本一区| 搡老乐熟女国产| 国产精品 欧美亚洲| 久久久久国产一级毛片高清牌| 中文欧美无线码| 精品久久久久久久毛片微露脸 | 国产有黄有色有爽视频| 亚洲精品美女久久av网站| 午夜福利,免费看| 欧美黄色片欧美黄色片| 久久毛片免费看一区二区三区| 新久久久久国产一级毛片| 97在线人人人人妻| 丁香六月天网| 在线观看人妻少妇| bbb黄色大片| 日韩有码中文字幕| 国产欧美日韩综合在线一区二区| 一区二区三区精品91| 国产黄频视频在线观看| 久久久国产成人免费| 在线精品无人区一区二区三| 一区福利在线观看| 国产主播在线观看一区二区| 在线永久观看黄色视频| 美女高潮喷水抽搐中文字幕| 久久久久精品人妻al黑| 午夜福利在线免费观看网站| 国产精品一区二区在线观看99| 久久天堂一区二区三区四区| 日韩欧美一区二区三区在线观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 日韩人妻精品一区2区三区| 亚洲美女黄色视频免费看| 久久精品久久久久久噜噜老黄| 久9热在线精品视频| 操出白浆在线播放| 久久国产精品大桥未久av| 久久99一区二区三区| 国产成人影院久久av| 欧美精品一区二区免费开放| 日本vs欧美在线观看视频| 亚洲免费av在线视频| 欧美大码av| 一区二区三区四区激情视频| 久久影院123| 我的亚洲天堂| 精品少妇内射三级| 丝袜脚勾引网站| 纯流量卡能插随身wifi吗| 日本91视频免费播放| 国产日韩欧美在线精品| 国产一卡二卡三卡精品| 国产在线视频一区二区| 少妇裸体淫交视频免费看高清 | 免费久久久久久久精品成人欧美视频| 国产极品粉嫩免费观看在线| 狂野欧美激情性xxxx| 成人三级做爰电影| 91国产中文字幕| 日本撒尿小便嘘嘘汇集6| 69精品国产乱码久久久| 欧美成狂野欧美在线观看| 99国产极品粉嫩在线观看| 亚洲成人免费电影在线观看| 美女视频免费永久观看网站| 男女午夜视频在线观看| 亚洲男人天堂网一区| 欧美xxⅹ黑人| 黄片大片在线免费观看| 女人精品久久久久毛片| 久久午夜综合久久蜜桃| 精品久久久精品久久久| 国产精品国产三级国产专区5o| 在线观看www视频免费| 精品亚洲成国产av| 国产一卡二卡三卡精品| 俄罗斯特黄特色一大片| 国产成人精品无人区| 大香蕉久久成人网| 各种免费的搞黄视频| 国产精品一区二区在线不卡| 大片免费播放器 马上看| 老熟妇乱子伦视频在线观看 | 两个人免费观看高清视频| 久久久水蜜桃国产精品网| 国产精品 国内视频| 人人妻人人澡人人看| 777米奇影视久久| 啦啦啦中文免费视频观看日本| 国产成人精品久久二区二区免费| 美女国产高潮福利片在线看| 女人被躁到高潮嗷嗷叫费观| 久久狼人影院| 久久久久精品人妻al黑| 久久99一区二区三区| 男人操女人黄网站| 爱豆传媒免费全集在线观看| 人人澡人人妻人| 亚洲一区中文字幕在线| 美女视频免费永久观看网站| 一区福利在线观看| 成人免费观看视频高清| 久久精品成人免费网站| 99热网站在线观看| 欧美大码av| www日本在线高清视频| 亚洲专区中文字幕在线| 蜜桃在线观看..| 国产有黄有色有爽视频| 欧美日韩一级在线毛片| 丝袜人妻中文字幕| 一进一出抽搐动态| 熟女少妇亚洲综合色aaa.| 久久青草综合色| 亚洲av片天天在线观看| 亚洲精品成人av观看孕妇| 丝瓜视频免费看黄片| 久久亚洲国产成人精品v| 中文字幕另类日韩欧美亚洲嫩草| 国产成人精品无人区| 国产欧美亚洲国产| 母亲3免费完整高清在线观看| 在线观看免费午夜福利视频| 国产精品 欧美亚洲| 91精品伊人久久大香线蕉| 国产精品一区二区精品视频观看| 日韩制服丝袜自拍偷拍| 一本一本久久a久久精品综合妖精| 俄罗斯特黄特色一大片| 99久久国产精品久久久| 欧美精品av麻豆av| 男男h啪啪无遮挡| 一二三四在线观看免费中文在| 国产亚洲精品第一综合不卡| 性高湖久久久久久久久免费观看| 国产福利在线免费观看视频| 国产一区二区三区av在线| 老司机福利观看| 精品少妇内射三级| 麻豆乱淫一区二区| 在线永久观看黄色视频| 国产成人精品在线电影| 国产一区二区在线观看av| 午夜免费鲁丝| 777久久人妻少妇嫩草av网站| 嫩草影视91久久| 午夜福利,免费看| 久久精品亚洲av国产电影网| 欧美精品av麻豆av| 欧美激情久久久久久爽电影 | 国产黄频视频在线观看| 亚洲第一欧美日韩一区二区三区 | 麻豆国产av国片精品| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看免费高清a一片| 一本大道久久a久久精品| 超碰97精品在线观看| 免费少妇av软件| 日韩视频一区二区在线观看| 最黄视频免费看| 久久久久久久精品精品| 久久国产精品影院| 日韩 欧美 亚洲 中文字幕| 国产亚洲欧美在线一区二区| 香蕉国产在线看| 国产高清国产精品国产三级| 欧美精品啪啪一区二区三区 | av欧美777| 精品国产乱码久久久久久男人| 成人手机av| 丰满人妻熟妇乱又伦精品不卡| 国产99久久九九免费精品| 国产精品99久久99久久久不卡| 黄色 视频免费看| 高清视频免费观看一区二区| 欧美精品人与动牲交sv欧美| 亚洲第一欧美日韩一区二区三区 | av福利片在线| 亚洲国产精品一区三区| 亚洲欧美一区二区三区久久| 12—13女人毛片做爰片一| 亚洲成人国产一区在线观看| av网站免费在线观看视频| 老司机福利观看| 亚洲欧美一区二区三区久久| 中文字幕制服av| 久久天堂一区二区三区四区| 亚洲一码二码三码区别大吗| 十分钟在线观看高清视频www| 精品一区二区三区av网在线观看 | 亚洲综合色网址| 一本综合久久免费| www.自偷自拍.com| 亚洲,欧美精品.| 电影成人av| 久久久国产一区二区| 国产在视频线精品| av线在线观看网站| 国产精品.久久久| 亚洲欧美成人综合另类久久久| 国产人伦9x9x在线观看| 老司机亚洲免费影院| 欧美精品高潮呻吟av久久| 国产有黄有色有爽视频| 日本wwww免费看| 啦啦啦中文免费视频观看日本| 精品福利观看| 久久狼人影院| 人成视频在线观看免费观看| 中文精品一卡2卡3卡4更新| 久久精品国产综合久久久| 伊人久久大香线蕉亚洲五| 欧美黑人欧美精品刺激| 欧美成狂野欧美在线观看| 制服人妻中文乱码| 午夜两性在线视频| 蜜桃在线观看..| 嫁个100分男人电影在线观看| 精品久久蜜臀av无| 国产精品熟女久久久久浪| 日本欧美视频一区| 久久人妻熟女aⅴ| 久久久久国产一级毛片高清牌| 欧美人与性动交α欧美软件| 嫩草影视91久久| 午夜福利免费观看在线| 老司机午夜十八禁免费视频| 真人做人爱边吃奶动态| 欧美+亚洲+日韩+国产| 国产精品九九99| 国产成人一区二区三区免费视频网站| 免费久久久久久久精品成人欧美视频| 丁香六月欧美| 99国产极品粉嫩在线观看| 国产免费福利视频在线观看| 国产精品亚洲av一区麻豆| 成年人午夜在线观看视频| 精品视频人人做人人爽| 中亚洲国语对白在线视频| av天堂久久9| 激情视频va一区二区三区| 高清视频免费观看一区二区| 制服诱惑二区| 脱女人内裤的视频| 自线自在国产av| 国产精品.久久久| 精品一区二区三区四区五区乱码| 欧美中文综合在线视频| 九色亚洲精品在线播放| 久久久精品免费免费高清| 日韩大码丰满熟妇| 淫妇啪啪啪对白视频 | 精品亚洲成国产av| 欧美精品人与动牲交sv欧美| 法律面前人人平等表现在哪些方面 | av电影中文网址| 91成人精品电影| 在线看a的网站| 亚洲精品乱久久久久久| 超碰成人久久| 久久狼人影院| 免费看十八禁软件| 亚洲成国产人片在线观看| 欧美日韩精品网址| 午夜福利免费观看在线| 久久ye,这里只有精品| 97在线人人人人妻| a在线观看视频网站| 女人爽到高潮嗷嗷叫在线视频| 国产不卡av网站在线观看| 日本wwww免费看| 国产视频一区二区在线看| 老司机影院毛片| 国产人伦9x9x在线观看| 国产成人欧美| 国产精品自产拍在线观看55亚洲 | 精品少妇一区二区三区视频日本电影| 久久久精品免费免费高清| 视频在线观看一区二区三区| 国产成人av激情在线播放| 狠狠精品人妻久久久久久综合| 欧美日本中文国产一区发布| 欧美日韩亚洲综合一区二区三区_| 在线观看免费午夜福利视频| 王馨瑶露胸无遮挡在线观看| av在线app专区| 永久免费av网站大全| 9色porny在线观看| 老熟妇仑乱视频hdxx| 人人妻人人澡人人爽人人夜夜| 香蕉国产在线看| 我要看黄色一级片免费的| 少妇被粗大的猛进出69影院| 亚洲美女黄色视频免费看| 免费一级毛片在线播放高清视频 | 免费不卡黄色视频| 国产老妇伦熟女老妇高清| 精品人妻在线不人妻| 久久久国产成人免费| 亚洲专区字幕在线| 午夜福利一区二区在线看| 欧美成狂野欧美在线观看| 色播在线永久视频| 久热这里只有精品99| 久久九九热精品免费| 欧美乱码精品一区二区三区| 午夜免费观看性视频| 美女主播在线视频| 夜夜骑夜夜射夜夜干| av视频免费观看在线观看| 亚洲天堂av无毛| 国产成人a∨麻豆精品| 十八禁高潮呻吟视频| 欧美精品一区二区大全| 中国国产av一级| 国产亚洲精品第一综合不卡| 日本a在线网址| 欧美国产精品一级二级三级| 飞空精品影院首页| 黄频高清免费视频| 黄片播放在线免费| 精品人妻一区二区三区麻豆| 啪啪无遮挡十八禁网站| 精品国产乱子伦一区二区三区 | a级毛片在线看网站| 久久人人爽av亚洲精品天堂| av天堂在线播放| 母亲3免费完整高清在线观看| 欧美精品av麻豆av| 伊人久久大香线蕉亚洲五| 纵有疾风起免费观看全集完整版| 国产片内射在线| 老熟女久久久| 麻豆av在线久日| 99久久综合免费| 国产主播在线观看一区二区| 男人操女人黄网站| 高清欧美精品videossex| 欧美老熟妇乱子伦牲交| 国产欧美日韩综合在线一区二区| 女人爽到高潮嗷嗷叫在线视频| 国产真人三级小视频在线观看| 国产精品久久久久久人妻精品电影 | 精品少妇黑人巨大在线播放| 99热网站在线观看| 少妇 在线观看| 天天添夜夜摸| 国产精品一区二区在线不卡| 亚洲精品一区蜜桃| 麻豆av在线久日| 一本久久精品| 男女床上黄色一级片免费看| 久久久久国产一级毛片高清牌| 97人妻天天添夜夜摸| videosex国产| av一本久久久久| 成人av一区二区三区在线看 | 国产精品免费视频内射| 色94色欧美一区二区| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩av久久| 一本一本久久a久久精品综合妖精| 久久精品亚洲av国产电影网| 亚洲色图综合在线观看| 国产在线观看jvid| 日本av免费视频播放| 老司机深夜福利视频在线观看 | 狂野欧美激情性xxxx| 欧美日韩黄片免| 亚洲精品中文字幕一二三四区 | 9色porny在线观看| 午夜激情久久久久久久| 亚洲国产欧美一区二区综合| 窝窝影院91人妻| 1024香蕉在线观看| 久久久久视频综合| 丝袜在线中文字幕| av超薄肉色丝袜交足视频| 日韩免费高清中文字幕av| 国产一区二区三区综合在线观看| 日韩大片免费观看网站| 黄网站色视频无遮挡免费观看| 亚洲av电影在线进入| 国产精品一区二区在线不卡| 日本91视频免费播放| 男人操女人黄网站| 在线观看舔阴道视频| 亚洲男人天堂网一区| 美女扒开内裤让男人捅视频| 精品亚洲乱码少妇综合久久| 久久久久国内视频| 无限看片的www在线观看| 国产不卡av网站在线观看| 99久久精品国产亚洲精品| 国产精品一二三区在线看| 国产亚洲精品久久久久5区| 丁香六月欧美| 国产91精品成人一区二区三区 | 这个男人来自地球电影免费观看| xxxhd国产人妻xxx| 巨乳人妻的诱惑在线观看| 午夜精品久久久久久毛片777| 亚洲精品国产av成人精品| 99精品欧美一区二区三区四区| 久久久久视频综合| av线在线观看网站| 91国产中文字幕| 亚洲专区国产一区二区| 大型av网站在线播放| 一本久久精品| www.999成人在线观看| 欧美中文综合在线视频| 无遮挡黄片免费观看| 精品国产乱码久久久久久小说| 在线观看舔阴道视频| 一二三四社区在线视频社区8| av又黄又爽大尺度在线免费看| 99热全是精品| 高清黄色对白视频在线免费看| 日本av手机在线免费观看| 巨乳人妻的诱惑在线观看| 亚洲国产av影院在线观看| 999久久久精品免费观看国产| 国产av精品麻豆| 国产1区2区3区精品| 亚洲成人国产一区在线观看| 国产男人的电影天堂91| 丁香六月天网| 亚洲欧美成人综合另类久久久| 成年人午夜在线观看视频| 亚洲人成77777在线视频| 看免费av毛片| 宅男免费午夜| 男女高潮啪啪啪动态图| 亚洲国产看品久久| 在线av久久热| 亚洲成人免费电影在线观看| 亚洲精品中文字幕一二三四区 | 国产一区二区 视频在线| 18禁观看日本| 欧美精品一区二区免费开放| 久热爱精品视频在线9| 亚洲av欧美aⅴ国产| 一区二区三区乱码不卡18| 一个人免费在线观看的高清视频 | 久久国产精品影院| 亚洲av日韩在线播放| 天天躁日日躁夜夜躁夜夜| 亚洲国产日韩一区二区| 9色porny在线观看| 国产又色又爽无遮挡免| 伦理电影免费视频| 在线永久观看黄色视频| www.精华液| 美女午夜性视频免费| 国产男人的电影天堂91| 波多野结衣av一区二区av| 女人精品久久久久毛片| 成人av一区二区三区在线看 | 亚洲专区字幕在线| 亚洲av电影在线进入| 亚洲中文日韩欧美视频| 成人免费观看视频高清| 国产av国产精品国产| 一二三四社区在线视频社区8| 亚洲精品国产色婷婷电影| 在线观看免费视频网站a站| 一级片'在线观看视频| 国产一区有黄有色的免费视频| 操出白浆在线播放| 12—13女人毛片做爰片一| 久久香蕉激情| 中文字幕制服av| 日韩有码中文字幕| 女人被躁到高潮嗷嗷叫费观| 欧美成狂野欧美在线观看| 激情视频va一区二区三区| 超色免费av| 亚洲三区欧美一区| 在线观看免费日韩欧美大片| 少妇 在线观看| 又大又爽又粗| 又紧又爽又黄一区二区| 亚洲精品久久久久久婷婷小说| 欧美 亚洲 国产 日韩一| 欧美日韩中文字幕国产精品一区二区三区 | 777米奇影视久久| 午夜影院在线不卡| 亚洲第一欧美日韩一区二区三区 | 男男h啪啪无遮挡| 大香蕉久久网| 美女福利国产在线| 午夜福利在线免费观看网站| 国产免费一区二区三区四区乱码| 啦啦啦 在线观看视频| 国产99久久九九免费精品| 19禁男女啪啪无遮挡网站| 亚洲色图综合在线观看| 高清视频免费观看一区二区| 一级毛片女人18水好多| 岛国毛片在线播放| 91精品伊人久久大香线蕉| 天天添夜夜摸| 亚洲va日本ⅴa欧美va伊人久久 | 日韩一卡2卡3卡4卡2021年| 亚洲精品一区蜜桃| 日本wwww免费看| 日本91视频免费播放| 国产男人的电影天堂91| 欧美变态另类bdsm刘玥| 国产日韩一区二区三区精品不卡| 久久久欧美国产精品| 成人三级做爰电影| 国产日韩欧美在线精品| 亚洲精华国产精华精| 亚洲欧美日韩高清在线视频 | 日本撒尿小便嘘嘘汇集6| 老司机在亚洲福利影院| 亚洲精品粉嫩美女一区| www.熟女人妻精品国产| 999久久久国产精品视频| 亚洲色图 男人天堂 中文字幕| 欧美成人午夜精品| 亚洲精品久久久久久婷婷小说| a 毛片基地| 一区二区三区精品91| 不卡一级毛片| 热99re8久久精品国产| 欧美精品高潮呻吟av久久| www.精华液| 亚洲第一青青草原| 一级毛片女人18水好多| 老司机影院毛片| 一边摸一边抽搐一进一出视频| 男女免费视频国产| 精品少妇一区二区三区视频日本电影| 纵有疾风起免费观看全集完整版| 色婷婷久久久亚洲欧美| 高清黄色对白视频在线免费看|