The following article is Open access

Copper Doped Zeolitic Imidazole Frameworks (ZIF-8): A New Generation of Single-Atom Catalyst for Oxygen Reduction Reaction in Alkaline Media

Published 6 August 2020 © 2020 The Author(s). Published on behalf of The Electrochemical Society by IOP Publishing Limited
, , Citation Anand Parkash 2020 J. Electrochem. Soc. 167 155504DOI 10.1149/1945-7111/abaaa5

1945-7111/167/15/155504

Abstract

Reasonable design of electrocatalyst based on abundant earth elements is of great importance for sustainable energy use. Transition metal nitrogen doped carbon (TM-N/C) materials demonstrate excellent catalytic activity. However, there is still a lack of comprehensive understanding of these TM-N/C-systems. Herein, we have developed an effective approach to develop highly active and stable Cu promoted N-doped carbon (Cu-N/C) catalysts. The impacts of copper doping and temperature of pyrolysis on catalytic performance have been studied. The Cu-NC-800 catalysts exhibited excellent catalytic activity and stability with an onset potential (Eo) of 0.99 V and half-wave potential (E1/2) of 0.85 V. It also shows strong the long-term stability. The synergistic effect between Cu (II)-N ligand and Cu0 nanoparticles is high, the active center is small, the transfer of mass is rapid, and the electrocatalytic efficiency is increased. The findings showed that the non-noble metal-based catalyst's bifunctional oxygen electrode activities prepared in this study were as high as those of commercial oxygen-based, noble metal catalysts.

Export citation and abstractBibTeXRIS

This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 License (CC BY, http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse of the work in any medium, provided the original work is properly cited.

Polymer electrolyte membrane fuel cell (PEMFC) has safety advantages, high power density, rapid start-up, low operating temperature, etc. The major problem of PEMFC in practical implementation stems from the sluggish kinetics of oxygen reducing reaction (ORR).15 Platinum (Pt) and Pt related metals for ORR are widely regarded as the strongest electrocatalysts.68 But the shortage of these costly metals has, however, seriously impeded their broad use.911 Hence, it is important to develop high-efficiency, low-cost multifunctional electrocatalysts, which can catalyze ORR.1215 To minimize the cost of a noble metal catalyst, the metal-free catalyst has been widely studied,16 and due to its low price and high stability,1720 porous carbon is the most promising catalyst.2123 Nevertheless, the electrocatalytic behavior of raw carbon (carbon nanotubes, graphene sheets, and nanostructured carbon) is difficult to equate with that of commercial precious metals/carbon.2430

The production of highly efficient electrocatalysts based on rich earth elements is important. In recent years, TM-N/C (Fe, Co, Ni) derived from the zeolite imidazolate (ZIF) containing nitrogen has a wide specific surface area, a tunable pore structure, and a rich nitrogen composition.3135 Fair modification of the coordination structure between the metal and nitrogen transition atoms will effectively adsorb and activate oxygen molecules, making the active site an excellent electrocatalyst for ORR.3640 However, to understand the process of ORR from a broader and deeper perspective, it is important to research other TM-N/C (such as Cu).4145 Theoretical estimation suggests that the activity of Cu based catalysts is greater than that of other transition metals (Fe, Co, Ni, etc).4650 Moreover, Cu's conductivity as the second-highest conductivity substance is just 6% lower than that of Ag, demonstrating the ability to facilitate the transfer of charges between the active center and the reactant.45,47 The key drawbacks of these copper macrocycles are reduced chemical and thermal stability, poor conductivity, and low activity to minimize oxygen.44,48 Researchers have recently discovered another way of directly breaking copper salts or nanostructures, carbon, and nitrogen precursors into an inert gas atmosphere to prepare copper nitrogen/carbon nanostructures, thus preparing hybrid nanocatalysts with improved stability, conductivity, and operation4.9,49

Herein, Cu promoted nitrogen-doped carbon Cu/ZIF-8 and pyrolyzed in the N2 atmosphere at different temperature (500 °C–900 °C) and termed as Cu–NC–X (X = 500 °C, 600 °C, 700 °C, 800 °C, and 900 °C). The Cu-NC-800 electrocatalyst demonstrated superior activity in 0.1 M KOH electrolyte, with an Eo of 0.99 V and E1/2 of 0.83 V compared to all prepared, as well as commercial catalysts, and close to the results recently reported. This also demonstrates excellent long-term stability. This impressive catalytic bifunctional performance can be attributed to the strong synergy between Cu (II)-N ligands and Cu0 nanoparticles, rich active centers, and rapid transfer of mass.

Experimental

Chemical reagents

Synthesis of CTAB-caped Cu NPs

Typically, copper nitrate trihydrate (Cu(NO3)2·3H2O, 1.5 g) and cetyltrimethylammonium bromide (CTAB, 106.6 mg) are typically added into methanol (50 ml). The mixture was stirred for 3 h at 75 °C. The Cu NPs were obtained by centrifugation, washed several times using methanol, then eventually vacuum-dried over 6 h at 80 °C.

Synthesis of ZIF-8

Typically, two solutions A and B are typically are synthesized. Solution A contains Zinc hexahydrate nitrate (Zn(NO3)2·6H2O, 0.9 g) dispersed in methanol (50 ml), and solution B contains 2-methylimidazole (11.0 g) dissolved in methanol (50 ml). The two solutions were then mixed under stirring for 0.5 h, and thus the resultant mixture then sterilized over 2 h at ambient temperature. The acquired purple precipitate was obtained through centrifugation, washed several times with methanol then vacuum-dried over 24 h at 80 °C.

Preparation of Cu-ZIF-8

For the preparation of Cu/ZIF-8, 40 ml ethanol, 0.210 g ZIF-8 was ultrasonically dispersed for three min, and then, after stirring, 10 ml ethanol containing 0.142 g and 0.210 g Cu NPS was added. Stir the mixture in a water bath at 400 rpm (80 °C) for 1 h. After the reaction was completed, the samples were centrifuged, washed three times with ethanol, and dried at 80 °C for 24 h. Prepared catalysts were named as Cu/ZIF-8. For the calcination, 250 mg of prepared catalysts was burned in the furnace at a flow rate of 20 ml min−1 of nitrogen for 2 h at different temperatures (500 °C–900 °C), at the rate of 5 °C min−1.

Results and Discussion

As shown in Scheme 1, the precursor Cu/ZIF-8 is pyrolyzed in an N2 atmosphere to synthesize the electrocatalyst Cu-NC. Briefly, CTAB caped nanoparticles were gradually dispersion on the surface of ZIF-8 and stirred at room temperature. The collected powder was then pyrolyzed to obtain the Cu-NC electrocatalyst directly at various temperatures in the N2 atmosphere, without any purification processes. The organic linker transforms into a carbon skeleton containing nitrogen during high-temperature pyrolysis, and the volatilization of Zn forms free nitrogen places. At the same time, the adsorbed copper ions are stabilized by combination with nitrogen and then reduced by the carbon around them.

Scheme 1. Refer to the following caption and surrounding text.

Scheme 1. Schematic illustration for the synthesis of Cu/ZIF-8.

Standard image High-resolution image

XRD characterization

The structural features of ZIF-8 and Cu doped ZIF-8 (Cu/ZIF-8) catalysts were explored by the X-ray diffraction (XRD) are shown in Fig. 2a. The synthesized Cu doped ZIF-8 morphology is like ZIF-8, suggesting the formation of a body-centered cubic crystal lattice symmetry. Also, the deposition of Cu nanoparticles to the reaction mixture of ZIF-8 crystals did not affect the frame structure of the ZIF-8 materials, and the structural integrity is maintained due to the close ionic size between Cu2+ (0.71 Å) and Zn2+ (0.74 Å) ions. Cu-NC and N/C show high absorption peaks of 23.4° and 43.5°, respectively corresponding to amorphous carbon and graphite carbon (Fig. 1b). For Cu-NC, the other two peaks centered at 43.2° and 50.4° correspond to metallic copper (111) and (200), suggesting the presence of zero-valent copper (Cu0). Due to the formation of graphite structure and amorphous carbon during pyrolysis, at 23° and 44°, respectively, there are two large peaks of Cu-NC and N–C, corresponding to carbon diffraction (002) and (101) (Fig. 1b). In particular, the two peaks at 43.4°, 50.4°, and 74.3° respectively correspond to diffraction Cu (111), (200) and (220), suggesting the presence of pure copper. The XRD analysis reveals that the copper addition has no apparent effect on the ZIF-8 structure. The X-ray diffraction spectrum does not show an apparent copper peak, which may be due to the relatively low copper content.

Figure 1. Refer to the following caption and surrounding text.

Figure 1. XRD patterns of pure (a) as well as calcinated ZIFs (b).

Standard image High-resolution image

Specific surface area characterization

The adsorption-desorption isotherm of N2 at 77 K was also measured to investigate the porosity characteristics of the sample (Fig. S1 is available online at stacks.iop.org/JES/167/155504/mmedia, Table I). Table I summarizes the characteristic parameters of samples, including an area of the BET surface and the total volume of the pores. The I-type isotherms of crystals with BET surface are 1437, 1230, 1177, 954, 765, 653, 632, and 621 m2 g−1, these are for ZIF-8, NC-900, Cu/ZIF-8, Cu/NC-500, Cu/NC-600, Cu/NC-700, Cu/NC-800, and Cu/NC-900, respectively, indicating microporous structure (Table I). The adsorption-desorption hysteresis curve of N2 near P/P0 = 1 is consistent with the porosity between particles, further demonstrating the dual microporosity and mesoporous existence of Cu/ZIF-8 crystal (Fig. S1). The BET specific surface area and pore volume of Cu-NC-900 nanocrystals are lower than other nanocrystals. Finally, the pore size distribution is centered at 5.63 nm, which is in accordance with the ZIF-8 structure model. Figure S1 reveals that the Cu-NC and N–C micropores are less than the Cu/ZIF-8 and ZIF-8 micropores, which could be caused by the collapse of certain structures.42 The heavy zinc evaporation and carbon gasification are mainly attributed to the N–C and Cu-NC mesopores and micropores, and the damage to the cubic structure is mainly attributed to certain macropores. Cu-NC has more mesopores and macropores relative to the matrix, whereas N–C has fewer mesopores and macropores.41 The particle size distribution (PSD) reveals that Cu-NC is a kind of hierarchical porous carbon with micropores, mesopores, and macropores in Fig. S1 and Table II. Three-dimensional porous hierarchical structure, combined with the synergistic effect of different pore sizes, may reduce resistance to mass transfer and reveal more ORR active sites.

Table I.  Structural properties of the materials analyzed by N2 adsorption-desorption isotherms.

Entry Catalysts SBET (m2 g−1)a) Pore Volume (cm3 g−1)b) Pore size (nm)c)
1 ZIF-8 1437 0.65 6.09
2 NC-900 1230 0.55 4.78
3 Cu/ZIF-8 1177 0.58 6.23
4 Cu/NC-500 954 0.18 6.11
5 Cu/NC-600 765 0.24 6.51
6 Cu/NC-700 653 0.19 7.23
7 Cu/NC-800 632 0.19 5.35
8 Cu/NC-900 621 0.19 5.63

a)Specific surface area according to BET. b)t-Plot micropore volume. c)BJH Adsorption average pore diameter (4 V/A).

Table II.  Electrocatalytic parameters of prepared catalysts.

Catalysts Eo (E vs RHE) E1/2 (E vs RHE) n
ZIF-8 0.92 0.79 2.7
NC-900 0.92 0.80 3.3
Cu/ZIF-8 0.93 0.81 3.4
Cu-NC-500 0.95 0.82 3.4
Cu-NC-600 0.97 0.83 3.5
Cu-NC-700 0.98 0.84 3.9
Cu-NC-800 0.99 0.85 4.1
Cu-NC-900 0.98 0.83 3.3
Pt/C (20 wt.%) 0.97 0.84 3.8

FE-SEM, TEM and HRTEM characterization

SEM images of samples were measured to compare the morphology (Figs. 2c and S3–S9). These images clearly show that all samples show a well-defined dodecahedral structure with truncated rhomboid, which is a typical ZIF-8 morphology as reported. These findings suggest that the ZIF-8 surface morphology and skeleton structure did not change during Cu nanoparticle doping. EDX spectrum analyzes the composition of the elements of the samples, and the findings are shown in Figs. 2d and S3–S9. These EDX results agree well with the copper content in theoretical calculation. EDX showed that after carbonization, the accurate content of copper increased. It is because zinc evaporates during pyrolysis, and some small molecules such as water and nitrogen compounds are released. The copper content increased as a result. Furthermore, the mapping of the elements is shown in Figs. 2e and S3–S9. The element mappings show the presence of carbon, nitrogen, and copper in the samples. The element mapping images show uniform copper, carbon, and nitrogen dispersion. This indicates copper was successfully incorporated into the material. The transmission electron microscope (TEM) image of Cu-NC-800 (Fig. 2a) shows a porous structure, which increases the exposure of the active center and the speed of mass transfer. The high-resolution transmission electron microscope (HRTEM) image reveals a 0.285 nm d-spacing, which belongs to the (111) crystal surface of Cu(N3)2 (Fig. 2b). As shown in Fig. 2, the metal active sites are evenly distributed on Cu-N/C, which after high temperature pyrolysis shows the loose and porous characteristics, while the composite metal nanoparticles are coated in the carbon layer or at the top of the carbon, so the ingenious structure is conducive to improving the stability of the catalyst while avoiding the active sites directly exposed in the electrolyte.4345 At the same time, the carbon layer is coated with active metal material, which avoids the direct contact between the active site and the electrolyte and improves the stability of the catalyst.

Figure 2. Refer to the following caption and surrounding text.

Figure 2. (a) TEM, (b) HRTEM, (c) FESEM images, (d) EDX spectra, (e) EDX mapping of images of Cu-NC-800.

Standard image High-resolution image

FTIR characterization

Figure S9a displays the FT-IR spectra of the ZIF-8 and the copper doped ZIF-8. A sharp peak at 423 cm−1 is caused by the Zn-N stretching mode, suggesting that the linker's zinc and the nitrogen atom are bound to form a porous ZIF coordination structure. The C–H vibration is assigned to the tensile frequency from 1100 to 1300 cm−1, while the C–N stretching vibration is assigned to the peaks at 1466 and 1470 cm−1. Additionally, the 1229 and 1231 cm−1 peaks confirmed C–N vibration in the imidazole ring. Cu-N bonding vibration of Cu doped ZIF-8 samples is within the 550–620 cm−1 range, which is consistent with the previous studies.

TGA characterization

ZIF-8 and Cu/ZIF-8 samples are just slightly different in the thermogravimetric analysis (TGA) (Fig. S9b). The TGA curve indicates only a very low weight loss of less than 0.5% from 20 to 200 °C, which corresponds to the removal from the cavity of guest molecules (methanol or Hmim and/or CO2). After guest-free Cu/Zn (MIM)2 crystal was formed at 350 °C, a long platform was observed, which indicated that the three-dimensional network of all samples had good thermal stability, like that found in the literature. Both curves indicate a step that can be assigned to the MIM linker decomposition above 350 °C (the start of exothermic decomposition).5,46 A sharp weight loss of 63%–64% is observed as the temperature increases from 400 °C to 500 °C, which is very consistent with the theoretical weight loss of 64%. The residue was XRD-doped zinc oxide doped by wurtzite copper. It should be noted that the thermal stability of Cu/ZIF-8 crystal obtained in this study is greater than that of copper-based MOFs (such as Cu (BDC)) used for catalytic use, which was decomposed at 200 °C. It should be noted that, during the heating, the zinc oxide content is still very low, and the connector is completely burnt out, suggesting that zinc oxide (ZnO) is the final calcined component of ZIF-8 nanocrystals.

XPS characterization

X-ray photoelectron spectroscopy (XPS) further analyzed the elementary structure and chemical state of the electrocatalysts Cu-NC-800. The XPS measurement continuum (Fig. S10) shows that Cu-NC-800 is composed of elements C, Cu, N, and O. The high-resolution N 1 s spectrum can be transformed into four peaks at 398.2 eV, 400.2 eV, 401.0 eV, and 403.3 eV, respectively, attributed to pyridine-N, pyrrole-N, graphite-N, and oxidation-N. The Cu0 is related to the main peak with a binding energy of 933.0 eV in the high-resolution Cu2p range. Cu (II) can be allocated with the shoulder peak at 934.6 eV and the satellite peak at 945.0 eV. Coordination with N atoms to form Cu (II)-N bond is suggested. Considering that Cu-NC 's ORR behavior is higher than N/C's, this means that Cu-NC's N content is rich and can play its role since the essence of nitrogen is not essential to ORR. The pyridine-N interaction with Cu (II) forms Cu (II)-N, while graphite-N affects the electronic and geometric carbon structure. It is reported that the Cu2+ center-coordinated MOF with nitrogen atoms such as pyrimidine and imidazole is an active catalyst for the aerobic oxidation of alkanes without any free radical initiator under mild reaction.2,7 We are continuing to research and find the high stability of ZIF-8 in this way. The goal of this work is to dope Cu nanoparticles in the ZIF-8 system and convert the inactive ZIF-8 into the active ZIF-8 to promote ORR without any free radical initiator/co-catalyst. The introduction of another metal Cu can increase the proportion of O–C–O/N–O and Cu oxygen bond, which greatly promotes its initial potential as ORR electrocatalyst.

Electrochemical characterization

Firstly, cyclic voltammetry (CV) was used to calculate ORR activity in 0.1 M KOH electrolyte saturated with oxygen or nitrogen. If the electrolyte is saturated with O2, Cu-NC has an apparent cathodic reduction peak, while N2 has no noticeable oxidation-reduction peak, suggesting Cu-NC has ORR activity. The results of the Cu-NC electrocatalyst were reported with the linear sweep voltammetry (LSV) to evaluate the ORR operation (Figs. 3a, 3b, 3d). From all prepared catalysts, Cu-NC-800 exhibited the onset potential (Eo 0.99 V) and the halfwave potential (E1/2 0.85 V), superior than that of all prepared ZIF-8 (Eo 0.92 V, E1/2 0.79 V), NC-900 (Eo 0.92 V, E1/2 0.80 V), Cu/ZIF-8 (Eo 0.93 V, E1/2 0.81 V), Cu-NC-500 (Eo 0.95 V, E1/2 0.82 V), Cu-NC-600 (Eo 0.97 V, E1/2 0.83 V), Cu-NC-700 (Eo 0.98 V, E1/2 0.84 V), Cu-NC-900 (Eo 0.98 V, E1/2 0.83 V), and commercial Pt/C (Eo 0.97 V, E1/2 0.84 V). Cu/ZIF-8 and Cu-N/C displayed the better ORR activity compared to ZIF-8 and NC-900, suggesting that Cu species introduction led to ORR activity. Furthermore, the optimization experiment shows that the synthesis of Cu-NC at 800 °C is the optimal parameter for achieving the optimum balance between Cu(II)-N and Cu0 NPs. The excellent catalytic efficiency can be due to the following factors: the porous structure is conducive to the transfer of electron and exposure of rich active centers; the transport of electrolyte and gas diffusion; nitrogen doping in carbon matrix also promotes catalyst activity significantly, and the incorporation of Cu nanoparticles is the key to improving catalytic efficiency.4,16,24

Figure 3. Refer to the following caption and surrounding text.

Figure 3. (a) ORR curves of all prepared catalysts compared with commercial Pt/C in an O2-saturated 0.1 M KOH, (b) electron transfer number (n) of all prepared catalysts and commercial Pt/C, (c) ORR curves of Cu-N/C-800 at 400, 800, 1200, 1600 and 2000 rpm, (d) Number of electron transfers of Cu-N/C-800, (e) The ORR curves of Cu-N/C-800 in O2-saturated 0.1 M KOH before and after accelerated aging tests at a scan rate of 10 mV s−1 under 1600 rpm.

Standard image High-resolution image

A change of the LSV polarization curve with rotation speed was recorded to understand better the ORR mechanism (Figs. 3c, 3e). The measured transferring electron number (n) of ZIF-8 (2.7), NC-900 (3.3), Cu/ZIF-8 (3.4), Cu-NC-500 (3.4), Cu-NC-600 (3.5), Cu-NC-700 (3.9), Cu-NC-900 (4.1), Cu-NC-900 (3.3), and commercial Pt/C (3.8), respectively according to the slope of the Koutecky-Levich (K-L) curve (Figs. 3c, 3e, Table II). Results show that nearly all catalysts are conducive to the creation of nearly four-electron channels, which is the most powerful ORR process.10,21,47 This enhanced activity may be due to the presence of high content of graphite nitrogen can improve the conductivity of carbon materials.12,42

The stability of catalyst is also very important in its practical application. Therefore, we test the stability of catalyst by chronoamperometry and accelerated ageing method. Firstly, the current density curves of Cu/NC-800 catalysts were measured in 0.612 V fixed potential by chronoamperometry and 0.1 M KOH electrolyte saturated with O2 at 1600 rpm. Under the same conditions, the chronoamperometric curves of Pt/C (20 wt%) and N/C-900 catalyst with the second-best performance were compared. As shown in Fig. 4, that the current density of Cu-N/C-800 catalyst remained at 98.5% after 10 h long-term test at a constant potential, followed by the stability of N/C-900 catalyst under the same test conditions, the current density decreased to 74.6% of the initial current density, while the current density decreased to 83.3% of the initial current density at Pt/C (20 wt%). Using the accelerated aging method, in oxygen saturated 0.1 M KOH electrolyte, at a sweep rate of 50 mV s−1, CV cycle scanning was carried out in the potential range of −0.8 ∼ 0.2 V (vs SCE) for 1000 cycles, and the ORR performance curves in oxygen saturated 0.1 M KOH electrolyte before and after scanning were recorded respectively. Results, as shown in Fig.3e, after 1000 cycles of CV aging, the ORR properties of Cu-N/C-800 in 0.1 M KOH electrolyte saturated with O2 remained unchanged. However, after 1000 cycles of CV aging, the current density and half-wave potential of Pt/C (20 wt%) decreased significantly, which also reflected the excellent stability of Cu-N/C-800 catalyst and exceeded the commercial Pt/C (20 wt%). In commercial Pt/C catalyst, aggregation and migration of Pt nanoparticles are important factors for the degradation of Pt nanoparticles.10,15,39 The strong covalent bond between carbon and nitrogen, and the protection of copper's carbon matrix, will prevent the degradation of the Cu-NC catalyst. The results suggest high stability for the porous carbon Cu-NC. The high-performance electrochemical activity of Cu-NC is mainly due to three factors. Secondly, Cu-NC has distributed active sites with adequate uniformity. Many researchers suppose that pyridine nitrogen and graphite nitrogen are ideal for electrochemical behavior. It has now been found that oxygen molecules are readily adsorbed on carbon atoms alongside pyridine nitrogen, as these carbon atoms are the active centers of ORR.11,23 Therefore, we assume the contribution of graphite nitrogen to the electrochemical activity of the two materials may be similar. Hence, we conclude that Cu-NC's high activity is attributable to the higher pyridine nitrogen content than N-C-900, and due to the sacrifice of the Cu-ZIF-8 prototype, the active center is evenly distributed in the material. Second, Cu-NC has an acceptable porous layered structure, and Cu manages its high large surface area. Further ORR active sites are exposed due to the proper structure, and the resistance to mass transfer is significantly reduced. Cu-NC has adequate ORR active site access, short channel length and good mass transfer capacity, and has high electrochemical activity.9,14 Second, the copper inlay increases the conductivity of Cu-NC, which is useful for operating ORR.

Figure 4. Refer to the following caption and surrounding text.

Figure 4. The chronoamperometric response of N/C-900, Cu-NC-800, and Pt/C (20 wt.%) in an O2-saturated 0.1 M KOH for 10 h.

Standard image High-resolution image

Conclusions

In conclusion, we have developed a new method to prepare Cu-NC-X (X = 500 °C–900 °C) nanoparticles by pyrolysis of ZIF-8 doped Cu nanoparticles in N2 atmosphere. The catalyst Cu-NC-800 has an outstanding ORR activity compared to all prepared as well as commercial Pt/C catalyst. In addition, the catalyst has good selectivity for the reduction of O2 to OH by four-electron transfer. Cu-NC-800 exhibits higher stability compared to Pt/C as well as NC-900. The improvement of the ORR efficiency of Cu-NC is mainly due to the increase of the content of pyridine N, the ordered porous structure of Cu, N doping, and the improvement of the conductivity. Metal-N co-doping strategy will significantly improve the catalyst's electrocatalytic operation, which provides a valuable guideline for the design of a new non-noble N-doped carbon-based ORR catalyst. This work will promote the development of effective, low cost, nanostructured non-noble metal electrocatalysts for renewable energy conversion and storage.

Acknowledgments

This author wishes to express his sincere thanks to the Lab facilities provided for this work in the School of Chemistry and Chemical Engineering, Shaanxi Normal University, Xi'an, China.

10.1149/1945-7111/abaaa5
undefined