Yan
Feng
*a,
Yuliang
Zhang
a,
Guixiang
Du
a,
Jingbo
Zhang
a and
Xiaohui
Qu
*bc
aCollege of Chemistry, Tianjin Normal University, Tianjin 300387, P. R. China. E-mail: hxxyfy@mail.tjnu.edu.cn
bEnergy Storage and Distributed Resources Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA. E-mail: xqu1981@gmail.com
cCollege of Engineering Research, University of California, Berkeley, California 94720, USA
First published on 20th June 2018
Lithium–sulfur (Li–S) batteries, which are well-known and much studied rechargeable batteries, are very promising because of their low cost, environmental friendliness, very high specific capacity and superior energy density. However, applications of Li–S batteries have been obstructed by their fast capacity fading and low coulombic efficiency due to soluble polysulfide migration during charge/discharge cycling. Herein, we present a strategy utilizing cobalt metal–organic framework (CoMOF) with rough porous surface and defective structure as a host material for sulfur accommodation, which implements a CoMOF and S composite (CoMOF–S) as a cathode. Capacity retention of CoMOF–S cathode surpasses that of pure sulfur by 87.18% after 100 cycles at 0.1C, while achieving coulombic efficiency above 98% at a high rate of 2.0C. We reveal structural features by combining X-ray photoelectron spectroscopy (XPS) and density functional theory (DFT) calculations, which confirm the covalent bond connection between CoMOF and S. We attribute the excellent cycling performance to covalent bond immobilization of sulfur.
To address the abovementioned challenges, metal–organic frameworks (MOFs) have been used to accommodate S and reduce dissolution of Li polysulfides. MOFs are formed by linking organic and inorganic groups through covalent coordination linkages. In addition, other types of interactions such as H-bonding, π–π stacking and van der Waals forces may also play important roles in shaping the three-dimensional (3D) structures of the MOFs.9 Because of their tunable chemical composition and highly porous frameworks, in which fast mass transportation of related species is favorable, MOFs have been proved to be particularly suitable for electrochemical applications.10–12 A pioneering study focused on the utilization of MOFs as a confined matrix for sulfur impregnation was performed by Tarascon et al.13 The group employed a highly porous and hydrothermally stable Cr-MOF (MIL-100(Cr)) to encapsulate sulfur via melt diffusion process. Capacity retention was improved by using Cr-MOF compared with mesoporous carbon as host. Subsequently, Qian et al. utilized Cu-MOF (HKUST-1) as host material to trap sulfur and thus diminish the dissolution problem.14 Following these studies, many other MOFs, such as Zn-MOF (ZIF-8),15,16 Al-MOF (MIL-53, NH2-MIL-53),15 Cr-MOF (MIL-101),17 Ni-MOF,18 and Cu-MOF (MOF-525)19 have been used as sulfur confinements. Nevertheless, there is still a need to find suitable MOFs as sulfur hosts to further improve the electrochemical performance to meet the demands of high-capacity and high-power energy Li–S batteries. Moreover, physical and chemical interactions between MOF and S still need to be explored.
In this study, we report the preparation of cobalt MOF (CoMOF) with rough porous surface and defective structure as a matrix host to accommodate and confine sulfur. By sulfur-melt-diffusion method,20,21 sulfur effectively impregnates and absorbs into the CoMOF structure. Electrochemical performance results show that the CoMOF and S composite exhibits high specific capacity, improved cycle life, and high coulombic efficiency during high rate charge/discharge process. In addition, we investigated the chemical interaction/connection between the CoMOF host and sulfur both by X-ray photoelectron spectroscopic (XPS) analysis and density functional theory (DFT) calculations.
BE = EC 1s + Ecorr | (1) |
(2) |
The structures of all the species were fully optimized using TPSSh density functional in combination with the 6-31+G(d,p) basis set.24 A larger basis set 6-311++G(d,p) was used in the single point energy calculation and orbital energy analysis. Orbital energies were extracted from canonical molecular orbitals. All quantum mechanical calculations were carried out using Gaussian 09.25
Fig. 1 SEM images of CoMOF (a–c) and CoMOFS2 (d–f) materials; (g) elemental mapping images of CoMOFS2. |
Fig. 2a shows the XRD patterns of CoMOF and CoMOF–S composite materials. After sulfur impregnation at 155 °C, the diffraction peaks of crystalline sulfur are clearly seen in the CoMOF–S composite materials. However, intensity of crystalline CoMOF peaks is low, mainly because a large number of sulfur crystal diffraction peaks are covered up by the MOF structure. Sulfur ratios in CoMOF–S composite materials were determined by TGA under N2 atmosphere, as shown in Fig. 2b. Below 200 °C, the observed mass loss most likely corresponds to the loss of physically absorbed water on the CoMOF material. From 200 °C to nearly 300 °C, there is almost no weight loss for CoMOF, while pure sulfur exhibits one-step weight loss. This indicates that the weight loss at 200–300 °C can only be attributed to the content of sulfur in the composite materials. Based on these findings, the calculated weights of sulfur in CoMOF–S composites are 39.0 wt%, 56.6 wt%, and 65.6 wt%, corresponding to CoMOFS1, CoMOFS2, and CoMOFS3, respectively.
Porosity of CoMOF and CoMOFS2 was estimated by N2 adsorption–desorption analysis, as shown in Fig. 2c. CoMOF and CoMOFS2 exhibit type-IV isotherms with a hysteresis loop in the relative pressure range of 0.8 to 1.0, indicative of inhomogeneous mesoporous structure. The pore size distribution curves of CoMOF and CoMOFS2, as shown in Fig. 2d and calculated from the adsorption branch by the BJH method, display mesoporous pores with typical sizes of 3–10 nm (CoMOF) and 3–30 nm (CoMOFS2). Both the BET surface area (4.106 m2 g−1) and BJH surface area (13.077 m2 g−1) of CoMOFS2 became larger than those of pristine CoMOF (0.874 m2 g−1 and 3.095 m2 g−1, respectively), while the average pore diameter of CoMOFS2 (3.428 nm) became slightly smaller than that of CoMOF (3.823 nm). This phenomenon may be caused by sulfur impregnation. In the heating process at 155 °C, some part of the sulfur is embedded into the CoMOF matrix, resulting in the decrease in pore diameter of CoMOF; the rest of the sulfur may cover the surface of CoMOF crystals, leading to amplification of surface area.
Fig. 3a shows the cycle performances of CoMOF–S cathodes at charge/discharge current density of 0.1C with cut-off voltage of 1.5–3.0 V. CoMOF–S cathodes display much better cycling performance than the pure sulfur cathode; among them, the CoMOFS2 cathode exhibits the best cycle stability and reversible capacity. Initial delithiation capacity of CoMOFS2 is 1522 mA h g−1, which is higher than those of CoMOFS1 (1429.1 mA h g−1) and CoMOFS3 (1515.5 mA h g−1). The reversible capacity of CoMOFS2 is above 800 mA h g−1 at 0.1C rate and capacity decay of CoMOFS2 cathode is 34.66% lower than that of pure sulfur cathode under the same charge/discharge conditions. All the above results indicate that CoMOF–S with appropriate sulfur content (56.6 wt%) improves electrochemical performance of sulfur cathodes. However, all CoMOF–S cathodes suffer severe irreversible capacity fading in the first 20 cycles, which is typically observed in other types of cathode materials.26,27 After 20 cycles, their capacity decay drastically decreased. Cycling stability of the CoMOF–S cathodes is much higher than that of pure sulfur cathodes. Fig. 3b shows the specific capacity–potential curves of the CoMOFS2 cathode at the 1st, 5th, 10th, 25th, 50th, and 100th cycles. A typical discharge process with two plateaus at around 2.3 and 2.1 V is observed. The upper voltage plateau at 2.3 V corresponds to conversion of S8 to long chain polysulfides (Li2S8, Li2S6, or Li2S4), while the lower voltage plateau at 2.1 V reflects the further reduction of high-order polysulfides to low-order polysulfides (Li2Sn, n < 4) and finally to insoluble lithium sulfides (Li2S2/Li2S).28 The CoMOFS2 cathode manifests excellent high-rate performance (Fig. 3c and d), decreasing slowly from the reversible 1290 mA h g−1 at 0.1C to 970, 840, 710, and 610 mA h g−1 at 0.2C, 0.5C, 1.0C, and 2.0C, respectively. In addition, a capacity of 853 mA h g−1 is obtained when current is recovered from 2.0C to 0.1C, indicating a highly reversible low fading in capacity. From Fig. 2c, we can clearly see that both the rate performances and coulombic efficiency of CoMOFS2 are remarkably higher than those of the pure sulfur cathode at the same charge/discharge rate. The coulombic efficiency is ∼93% at low rate of 0.1C, which may be caused by the severe “shuttle effect”, where polysulfides diffuse through the electrolyte and migrate between the cathode and anode. However, the coulombic efficiency is ∼98% and ∼97% at high rate of 1.0C and 2.0C, respectively, which indicates that the CoMOF structure can stabilize some of the sulfur species and sulfur inside the MOF can be accessed and utilized efficiently at high current rate. In addition, the role of CoMOF structure is only as the sulfur host; it provides almost no capacity in the CoMOF–S cathode, which can be proven by the cycling performance of pristine CoMOF, as shown in Fig. S4 of the ESI.†
Cyclic voltammetry (CV) profiles of CoMOFS2 cathodes at a scan rate of 0.02 mV s−1 for the first three cycles are shown in Fig. 4a. The curve shapes of the first cycle are different from the subsequent cycles, which is attributed to the irreversible reaction involved in the first cycle. During the positive scan in the second scanning, two anodic peaks are observed at ∼2.5 and ∼2.7 V. The former peak indicates the oxidation reaction to form sulfur species is ∼2.5 V and the latter peak indicates some side reaction, such as oxidation of the carbon black at ∼2.7 V. During the negative scan, two cathodic peaks appear at ∼2.4 and ∼2.0 V, showing that the reduction reaction is around these two potentials. The peak at ∼2.4 V indicates the transformation of S to higher-order Li2Sx (4 ≤ x ≤ 8) and the peak at ∼2.0 V was caused by further reduction to higher-order Li2Sx (x ≤ 4) and finally to Li2S.29,30 After scanning downward to 1.5 V, one clear reduction peak at ∼1.6 V is found, which may be caused by reaction of LiNO3 in the electrolyte. The same situations are also found in the second and third cycles.
As shown in Fig. 4b, EIS was conducted to evaluate the electrochemical kinetic properties of the CoMOF–S cathodes and their equivalent circuit for EIS data fitting. The Nyquist plots for all the electrodes consist of one semicircle and a straight line, which are respectively associated with charge transfer and solid state diffusion process of lithium ions in the electrode and the Warburg region.31,32 EIS data were fitted to an electrochemical model using Zview module software, as shown in the inset in Fig. 4b. Rs at the high-frequency intercept on the real axis represents bulk resistance, and is dependent on electrolyte solution, active material, current collector and so on. Rp is the charge transfer resistance at the electrode–electrolyte interface, and is associated with the interface between electrolyte and active particles. CPE is the constant phase element and represents the conventional double-layer and passivation film capacitance. W is the finite length Warburg impedance with an open circuit terminus.33,34Table 1 shows the EIS fitting results of CoMOF–S cathodes. No distinct difference exists in the Rs values of the four samples, which proves that the cells have the same test conditions. Rp value of the CoMOFS2 cathode is the smallest among the CoMOF–S cathodes, which indicates the highest electrochemical activity and best utilization of electrodes among them. Comparing the Rp values in Table 1, CoMOFS2 cathode has the lowest Rp value, which indicates that lithium diffusion in CoMOFS2 cathode is the fastest among the CoMOF–S cathodes. This further demonstrates that the CoMOFS2 cathode has good lithium intercalation kinetics, corresponding to improved rate performance.
Parameter | CoMOFS1 | CoMOFS2 | CoMOFS3 |
---|---|---|---|
R s/Ω | 2.81 | 1.09 | 0.34 |
R p/Ω | 278.90 | 139.20 | 201.60 |
CPE-T/μF | 23.55 | 10.58 | 8.64 |
CPE-P/F | 0.67 | 0.70 | 0.73 |
W-R/Ω | 143.60 | 147.70 | 65.92 |
W-T/Ω | 0.39 | 0.38 | 0.11 |
W-P/Ω | 0.72 | 0.71 | 0.47 |
To further investigate the capacity decay during charge/discharge cycling, SEM images of CoMOFS2 cathode before cycling and after 10, 20, 50, and 100 cycles were recorded, as shown in Fig. 5. Before cycling, the CoMOFS2 particles are in tight contact with carbon black powder to form a homogenous and smooth electrode surface (Fig. 5a). With increasing cycling, the layer on the outside of CoMOFS2 cathode may undergo deposition of polysulfides (Fig. 5b–e). From the EIS analysis, as shown in Fig. 5f, we can draw the same conclusion. The single depressed semicircles in the high-to-medium frequency region can be ascribed to the charge-transfer resistance Rp, which represents the electrode reaction kinetics. After extended cycling from 1 to 20 cycles, the charge transfer resistance of CoMOFS2 decreases tremendously. This decrease is ascribed to loss of active materials upon cycling, which leads to capacity decay. With continued cycling, charge transfer resistance remains unchanged, which evidences the stability of the electrode after 20 cycles. This phenomenon matches well with the charge/discharge cycling performance shown in Fig. 3a.
To further identify the chemical status of the CoMOF–S material and reveal the reason for its improvement in electrochemical performance, X-ray photoelectron spectroscopy (XPS) measurements and density functional theory (DFT) calculations were performed. As shown in Fig. 6a, the C 1s spectrum of the CoMOFS2 material exhibits three peaks with binding energies of 284.5, 285.7, and 288.6 eV that correspond to carbon atoms in different functional groups of the MOF structure.35 Evidently, the fourth peak with binding energy of 286.5 eV is related to the structure of MOF as a result of sulfur binding, which is absent from the XPS of pristine CoMOF (Fig. 6b). To match the XPS peaks with the corresponding carbon atoms, DFT calculations were used to predict binding energy for each carbon. In agreement with four distinct XPS peaks, there are four unique carbon atoms, as shown in Fig. 6c. The ortho and piso aromatic carbons contribute two theoretical peaks at 284.7 and 285.2 eV, respectively, and are responsible for the 284.5 and 285.7 eV peaks observed experimentally. Ortho and piso aromatic carbons are also the primary carbon types in pristine CoMOF, which explains the duplication with pristine CoMOF for these two peaks. The other two carbons, namely, [OC–S]− and OC–O–S carbons, are associated with peaks at 286.8 and 288.7 eV, respectively. The OC–O–S carbon is responsible for the peak at 288.6 eV, which overlaps with the OC–O–H carbon peak in pristine CoMOF. The additional peak at 286.5 eV emitted by CoMOF–S is assigned to the [OC–S]− group, which deviates from theoretical value (286.8) by only 0.3 eV. Deviations of the DFT predictions for all CoMOF–S XPS peaks are within 0.5 eV. It is worth noting that the predicted accuracy is better for pristine CoMOF. Compared with the experimental values (284.7, 285.5 and 288.7 eV) in Fig. 6b, the maximum deviation from theoretical values (284.9, 285.3 and 288.4 eV) is only 0.3 eV. This excellent consistency is strong evidence that the proposed structure model in Fig. 5c exemplifies the connection mode between CoMOF and S. It is unambiguous that sulfur bonds with carbon to form strong covalent bonds.
We want to emphasize that we enumerated all possible structures for CoMOF and S composite, including covalent bond, hydrogen bond and van der Waals structures. The structure presented in Fig. 5c is the one with the best agreement with experiment. We computed the C 1s orbital energy using Koopmans' theorem22,36 at TPSSh/6-311++G(d,p) level,24 which was then corrected by an empirical parameter to obtain the XPS binding energy. The detailed computational methods are described in the methods section. As shown in Fig. 6c, we proposed two schemes for sulfur connecting to CoMOF. (1) Substitution for carboxylic –OH group, which changes a OC–O–H to a OC–S–H or a deprotonated [OC–S]− group. DFT simulation reveals that the energy of the deprotonated state is lower than that of the protonated state by 12.33 kJ mol−1. Therefore, the structure with [OC–S]− group is more energetically favorable and will be used in the XPS predictions. (2) Substitution for H atom of carboxyl group or H2O ligand. The combination of two substitutions results in a cross-link inside the pores of CoMOF. The cross-links tilt the carboxylic group plane slightly due to spatial constraints. In addition to direct O atom substitution, the structure distortion will also contribute to atomic charge redistribution. Natural bond orbital (NBO)37 population analysis shows that the atomic charge on carbon atom of OC–S–H is 0.85, which is close to that of the carbon of OC–O–H (0.82) in pristine CoMOF. However, the charge on carbon atom of [OC–S]− is sharply decreased to 0.28, which indicates higher electron density. Consequently, binding energy drops as a result of increased shielding of the nuclear charge. This can be seen more intuitively from the electrostatic potential in Fig. 6c. The [OC–S]− group resides in an area with significantly more negative potential than OC–S–H. As a consequence, it exerts less attraction to the electrons, which leads to lower binding energy. Both substitution schemes immobilize the sulfur atoms to the CoMOF pores by strong covalent bond. We expect that these schemes will facilitate sulfur confinement in the CoMOF framework.
Footnote |
† Electronic supplementary information (ESI) available: SEM images of pristine CoMOF after heating at 155 °C. XRD patterns of pristine CoMOF before and after heating at 155 °C. SEM images of CoMOFS1 and CoMOFS3. Cycle performance of pristine CoMOF (without S) cathode at a charge/discharge current density of 167.5 mA h g−1. See DOI: 10.1039/c8se00195b |
This journal is © The Royal Society of Chemistry 2018 |