Selenium cathodes have attracted considerable attention due to high electronic conductivity and volumetric capacity comparable to sulphur cathodes. However, practical development of lithium-selenium batteries has been hindered by the low selenium reaction activity with lithium, high volume changes and rapid capacity fading caused by the shuttle effect of polyselenides. Recently, single atom catalysts have attracted extensive interests in electrochemical energy conversion and storage because of unique electronic and structural properties, maximum atom-utilization efficiency, and outstanding catalytic performances. In this work, we developed a facile route to synthesize cobalt single atoms/nitrogen-doped hollow porous carbon (CoSA-HC). The cobalt single atoms can activate selenium reactivity and immobilize selenium and polyselenides. The as-prepared selenium-carbon (Se@CoSA-HC) cathodes deliver a high discharge capacity, a superior rate capability, and excellent cycling stability with a Coulombic efficiency of ~100%. This work could open an avenue for achieving long cycle life and high-power lithium-selenium batteries.
Rechargeable lithium-ion batteries (LIBs) are considered to be the promising candidates towards sustainable energy storage devices due to its long cycle life, high specific power and energy density1,2. However, the energy density of current LIBs can not meet the ever-increasing demands from many emerging applications such as electric vehicles3. On one hand, lithiumâ€“sulfur (Liâ€“S) batteries have attracted growing attention because of several advantages such as the natural abundance of sulfur, high specific energy density (2600â€‰Wâ€‰hâ€‰kgâˆ’1) and high theoretical capacity (1675â€‰mAâ€‰hâ€‰gâˆ’1)4. However, the development of Liâ€“S batteries still suffers from the inherent issues of low electronic conductivity of sulfur and the shuttle effect of polysulfides. As an element in the same group of sulfur in the periodic table, selenium owns similar chemical properties to sulfur and has been considered as an alternative cathode material for lithiumâ€“selenium battery because of its high theoretical volumetric capacity (3253â€‰mAâ€‰h cmâˆ’3)5,6. Additionally, the conductivity of Se (1â€‰Ã—â€‰10âˆ’3â€‰Sâ€‰mâˆ’1) is much higher than that of S (5â€‰Ã—â€‰10âˆ’30â€‰Sâ€‰mâˆ’1), which enables higher active material utilization and better rate capability7. However, the Se cathodes also have a dissolution issue associated with high-order lithium selenides (Li2Sex, xâ€‰>â€‰4) and large volume expansion during the charge/discharge process, resulting in a low Se utilization, inferior capacity and short cycle life7,8,9.
From the pioneering work by Amine et al.10,11, various strategies have been proposed to improve the electrochemical performance of selenium cathode. The most effective method is to incorporate Se particles with electronically conductive materials and encapsulate Se particles within a porous carbon matrix6,12,13. For Se/porous carbon composite materials, the charge transfer resistance can be decreased and the shuttle effect of polyselenides can be suppressed because of the high conductivity of carbon matrix and the strong affinity of porous carbon for Se particles7. Many porous carbon materials have been studied to construct Se/porous carbon composites for Liâ€“Se batteries, such as carbon nanospheres14,15, carbon nanofibers16, hierarchical porous carbon17 and porous hollow carbon bubbles18. However, high-power Liâ€“Se batteries with long cycling performance under high currents have never been reported due to the unsatisfactory performance of Se cathodes.
Single-atom catalysts (SACs) consist of isolated metal atoms dispersed or anchored on matrix materials. SACs attracted extensive attention due to their maximum atom utilization efficiency, homogenous active centres, and unique reaction mechanisms19,20,21. Currently, SACs have been successfully applied in batteries, including metal-air batteries and metal sulfur batteries22,23,24,25,26. In addition, the metal-organic framework (MOF)-derived SACs have been intensively investigated in the area of electrocatalysts because of their high electrical conductivity, superior activity, and maximum atomic utilization27,28,29,30. However, based on literature reviews, as a result of the great challenge of controllable synthesis atomic metals with a selenium host, there has been no report of single atoms in Liâ€“Se batteries, where single atoms can maximize the multi-functions of a selenium host to achieve high rate and cycling performance in a Liâ€“Se battery.
Herein, for the first time, we demonstrate that single-atom catalysts can enable highly effective cathodes for Liâ€“Se batteries with superior rate capability and outstanding long-term cycling performance. A facile and straightforward approach (Fig.Â 1) facilitates the delicate control of Zeolitic Imidazolate Framework (ZIF) particles deposited on the surface of polystyrene (PS) spheres. More importantly, the coreâ€“shell ZIF hybrid structure can be further converted into hollow structured carbon materials via a pyrolysis process with zinc evaporation. Through finely tuning the ratio between Zn and Co, we successfully prepared atomic cobalt electrocatalyst/nitrogen-doped hollow porous carbon (CoSA-HC), nitrogen-doped hollow porous carbon (HC) and cobalt nanoparticles/nitrogen-doped hollow porous carbon (CoNP-HC). In addition, by embedding Se in hollow structured carbon particles, carbon/selenium composites (Se@CoSA-HC) was obtained. When applied as cathode materials for Liâ€“Se batteries, the Se@CoSA-HC cathode exhibited a superior electrochemical performance, including a superior rate capability (311â€‰mAâ€‰hâ€‰gâˆ’1 at 50â€‰C) and excellent cycling stability (267â€‰mAâ€‰hâ€‰gâˆ’1 after 5000 cycles with a 0.0067% capacity decay per cycle at a current density of 50â€‰C) with the Coulombic efficiency of ~100%. This work reveals that the maximal utilization of cobalt single atoms can optimize the features of porous carbon materials towards the activation of selenium reactivity and immobilization of selenium and polyselenides. Our results demonstrate that the Se@CoSA-HC composite is a promising cathode material for lithiumâ€“selenium batteries with long cycle life and high-power.
The synthesis strategy for PS@ZIFs is schematically depicted in Fig.Â 1. After mixing the PVP-modified PS spheres with zinc nitrate, cobalt nitrate, and 2-methyl-imidazolate successively in methanol solution at room temperature, PS@ZIF particles and small ZnCo-ZIFs clusters are formed on the surface of the PS spheres. Three samples with Zn/Co molar ratios of 20:1, 21:0 and 17:4 have been obtained, designated as PS@ZIF-1, PS@ZIF-2, and PS@ZIF-3, respectively. The growth of ZIF particles on the surface of PS spheres is determined by the interfacial reaction between the substrate and surfactant. Polyvinylpyrrolidone (PVP) plays an important role in the homogeneous growth of ZIFs particles on the surface of PS spheres. To improve electrostatic forces and provide enough coordination sites to uniformly adsorb metal ions, the negatively charged PS spheres are firstly enriched with PVP molecules. The amide carbonyl groups of PVP can fully coordinate with metal ions through chemical bonds and make it possible to achieve the coating of ZIFs on the surfaces of the PS spheres31. As shown in Supplementary Fig.Â 1a, the particle size of the prepared PS@ZIF-1 is about 700â€‰nm and the particle size of the polyhedral-structured ZIFs on the surface of PS spheres is ~200â€‰nm. The XRD result for PS@ZIF-1 shown in Supplementary Fig.Â 2 agrees well with that of as-prepared bimetallic ZnCo-ZIFs32, confirming the formation of bimetallic ZIF structures. This growth strategy can be extended to the fabrication of other types of ZIFs with different dimensions and components. The SEM images and XRD patterns in Supplementary Fig.Â 3 show the successful preparation of rGO@ZIF and MnO2@ZIF. The as-prepared PS@ZIF materials were converted into carbonaceous nanocomposites via one-step pyrolysis. The PS template was removed in situ by evaporation to form a hollow morphology at 700â€‰Â°C under a N2 atmosphere. Owing to the composition and size of cobalt in the three products, the obtained carbon materials are denoted Co single atom/nitrogen-doped hollow porous carbon (CoSA-HC), nitrogen-doped hollow carbon (HC) and cobalt nanoparticle/nitrogen-doped hollow porous carbon (CoNP-HC), respectively. In the previous reports, hollow-structured materials have shown excellent electrochemical performance in energy storage due to large interior voids, high surface area, and shortened mass/charge transport lengths33,34,35,36,37. As shown in the SEM image (Supplementary Fig.Â 4), CoSA-HC particles with a decreased particle size of the polyhedral ZIF particle ~130â€‰nm were achieved. TEM images (Fig.Â 2a) and high angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) images (Fig.Â 2b) also confirmed the hollow nature of the as-prepared carbon materials and no cobalt aggregates are observed, suggesting that the cobalt nanoparticles are anchored uniformly on the carbon matrix. Energy-dispersive X-ray spectroscopy (EDS) images were obtained to identify the elemental distribution of cobalt, nitrogen, oxygen, and carbon. As shown in Fig.Â 2c, elements are homogenously distributed in the whole carbon framework. From aberration-corrected HAADF-STEM images in Fig.Â 2d and Supplementary Fig.Â 5, high-density bright dots (highlighted by red circles) can be detected, revealing the formation of atomic cobalt. The generation of atomically dispersed Co catalysts was attributed to the Co ions being reduced by carbonized organic linkers with the evaporation of elemental Zn during the calcination process28,38. The TEM images in Supplementary Fig.Â 6 and Supplementary Fig.Â 7 show that the generation of HC particles and CoNP-HC particles from pyrolyzed bimetallic ZnCo-ZIFs with Zn/Co molar ratios of 21:0 and 17:4. It is revealed that the distance of adjacent Co atoms can be expanded through the introduction of zinc species, achieving the controllable synthesis of different aggregation degrees of atomic Co.
X-ray diffraction (XRD) measurement was performed to investigate the phase structure of CoSA-HC, HC and CoNP-HC particles (Supplementary Fig.Â 8). The peaks around 26.8Â° and 44.1Â° in the XRD pattern of CoSA-HC and HC particles fit the (002) and (100) planes of graphitic carbon. After increasing the cobalt content in the CoNP-HC particles, newly formed peaks are observed at 44.6Â°, 51.9Â° and 76.7Â°, which can be ascribed to cubic cobalt metal (JCPDS 15-0806). The porosity of the CoSA-HC particles was then studied by nitrogen adsorption-desorption analysis. The CoSA-HC displays mixed type I and IV isotherms (Supplementary Fig.Â 9a). The sharp adsorption increase in the low-pressure region indicates that the existence of abundant micropores and the hysteresis loop in the medium-pressure region originates from the mesopores with the carbon framework. The pore size distribution of CoSA-HC particles was then analysed based on the isotherms (Supplementary Fig.Â 9b), which shows peaks centred at 1, 17 and 35â€‰nm. The Brunauerâ€“Emmettâ€“Teller (BET) surface area obtained on the basis of the aforementioned adsorption isotherm of CoSA-HC, HC and CoNP-HC are 221, 265 and 136â€‰m2â€‰gâˆ’1, respectively. The graphitization degree of a carbon material can be characterized by an ID/IG value (ID and IG represents the intensity of Raman D- and G- bands of carbon materials) in Raman spectra (Supplementary Fig.Â 10). The ratios for CoSA-HC, HC and CoNP-HC were 0.99, 1.00 and 1.01, respectively, indicating the similar graphitic structures in these three hybrids. This allows the exclusion of the influence of carbon supports on electrochemical performance, allowing us to concentrate instead on variances of Co particle size and configuration.
Selenium encapsulation into carbon particles was carried out at 300â€‰Â°C under Ar atmosphere to promote the infusion of Se into the microporous carbon clusters. The obtained selenium carbon composites are denoted Se@CoSA-HC, Se@HC and Se@CoNP-HC, respectively. To obtain a high selenium content in the composite materials, selenium powders were mixed with HC particles in a weight ratio of 1:3, respectively. SEM images (Supplementary Fig.Â 11a) and TEM images (Supplementary Fig.Â 11b) of the as-formed selenium composite material (Se@CoSA-HC) revealed the formation of a size-reduced structure. The HAADF-STEM and EDS images in Supplementary Fig.Â 11c suggest Se atoms are dispersed uniformly throughout the hollow carbon microstructure. The diffraction peaks of Se in Supplementary Fig.Â 12 correspond to trigonal crystalline (JCPDS 06-0362) and the XRD pattern of Se@CoSA-HC (Supplementary Fig.Â 12) indicated the successful incorporation of Se in the HC carbon matrix. After the selenization process, the crystallinity of selenium is low, which can be attributed to the transformation of trigonal Se to amorphous Se with the formation of low ordering of selenium and successful confinement in the carbon matrix7,8,9. The Se contents in Se@CoSA-HC, Se@HC, Se@CoNP-HC and Se@CoSA-HC with high Se loading were determined to be 57, 55, 57 and 73â€‰wt%, respectively, through thermogravimetric analysis (TGA) (Supplementary Fig.Â 13).
XPS measurements were also used to explore the chemical environment of CoSA-HC particles. The XPS high-resolution C 1â€‰s spectrum of CoSA-HC particles is shown in Supplementary Fig.Â 14 and can be divided into four individual peaks corresponding to Câ€“C (284.8â€‰eV), Câ€“N (286.2â€‰eV), Câ€“O (287.8â€‰eV) and C=O (289.3â€‰eV)39,40. In the high-resolution Co 2p spectrum (Fig.Â 3a), the deconvoluted peaks at 781.4â€‰eV and 796.6â€‰eV can be ascribed to Co 2p3/2 and Co 2p1/2 orbitals of Co2+ species and the peaks at 780.0â€‰eV and 795.5â€‰eV correspond to Co 2p3/2 and Co 2p1/2 orbitals of Co3+ species28,38. The peaks at 784.6 and 802.8â€‰eV are satellite peaks that can be ascribed to the shakeup excitation of the high-spin Co2+ ions36. The N 1â€‰s spectrum in Fig.Â 3b can be deconvoluted into four peaks located at 398.6, 399.3, 400.7 and 401.8â€‰eV, corresponding to pyridinic N, Coâ€“N, pyrrolic N, and quaternary-N, respectively41,42,43. X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure measurements (EXAFS) were conducted to investigate the chemical state and coordination environment of Co atoms in the CoSA-HC particles. As shown in the XANES spectra (Fig.Â 3c), the comparison of near-edge adsorption energy with reference cobalt phthalocyanine (CoPc) and Co foil implies that the Co single atoms in CoSA-HC particles are positively charged, which agrees well with previous results38,44. Furthermore, from the EXAFS spectra of the three samples (Fig.Â 3d), a main peak at ~1.48â€‰Ã… can be observed for CoPc, which is typically assigned to the Coâ€“N tetrahedral coordination (Coâ€“N4)44. The Coâ€“N coordination peak shifts to a low R-position at 1.43â€‰Ã… in the HC particles, revealing a slightly variation of Coâ€“N coordination. Compared with the spectra of HC particles and Co foil, no Co-Co peak around 2.1â€‰Ã… is observed in CoSA-HC particles, indicating atomically dispersed Co single atoms45,46. Furthermore, according to the fitting parameters given in Supplementary TableÂ 1, the Coâ€“N coordination number of CoSA-HC particles is 3.3, implying that the Coâ€“N interaction within HC particles are consisted of Coâ€“N tridentate (Coâ€“N3) and tetrahedral (Coâ€“N4) coordination. Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES) analysis of CoSA-HC particles indicates that the cobalt content was determined to be ~1.3%.
FigureÂ 4a reveals the cyclic voltammogram (CV) curves of the first three cycles of a Se@CoSA-HC electrodes at a scan rate of 0.1â€‰mV/s between 1.0â€‰V and 3.0â€‰V. During the first discharge process, two peaks at around 1.74 and 1.78â€‰V are observed, which may be related to the lithiation of different Se molecules in the Se@CoSA-HC composites47,48. However, these two peaks disappear in the second discharge process, simultaneous with the appearance of two reductions peak at around 1.81 and 2.02â€‰V, resulting from the electrochemical activation behaviour of the Se@CoSA-HC electrode during the lithiation process16,49. During the charge process, there is only one single anodic peak at 2.07â€‰V in all cycles, and this peak remains stable during following lithiation/delithiation cycles. The CV curves after the second cycle are overlapping, demonstrating the good electrochemical stability of the Se@CoSA-HC electrodes. As illustrated in Supplementary Fig.Â 15, the CV curves after the 10th and 100th cycle overlap very well at a higher scan rate of 2.0â€‰mV/s. The inset of Fig.Â 4a shows the second-cycle galvanostatic discharge-charge voltage profiles of Se@CoSA-HC at the current density of 0.1â€‰C. Both discharge and charge processes exhibit stable voltage plateaus at around 2.0â€‰V, which is consistent with the characteristic peaks from the CV curves. The specific capacities of Se@CoSA-HC electrodes at different current densities as a function of cycle number (from the second cycle) are presented in Supplementary Fig.Â 16a. The Se@CoSA-HC electrodes exhibit high specific capacities of 613, 579, 569, 548, 516, 467, 427, 385 and 311â€‰mAâ€‰hâ€‰gâˆ’1 at current densities of 0.1, 0.2, 0.5, 1, 2, 5, 10, 20 and 50â€‰C, respectively, which are higher than those of Se@HC and Se@CoNP-HC (Supplementary Fig.Â 16b). When the current rate was returned to 1â€‰C, a high discharge capacity of 537â€‰mAâ€‰hâ€‰gâˆ’1 could still be attained with almost no capacity degradation, suggesting excellent stability of Se@CoSA-HC cathodes. Supplementary FigureÂ 17 shows the discharge and charge profiles of Se@CoSA-HC, Se@HC and Se@CoNP-HC at various rates and cycles. It indicates that both discharge and charge processes exhibit stable voltage plateaus at around 2.0â€‰V. Compared with Se/rGO and Se/MnO2 electrodes (Supplementary Fig.Â 18), Se@CoSA-HC cathodes exhibit superior rate electrochemical performances mainly due to the enhanced conductive network, which was previously reported using rGO-based electrodes for Liâ€“Se batteries50. Supplementary FigureÂ 19 exhibits the rate capability of an unselenated CoSA-HC electrode at different current densities from 0.1 to 5â€‰C. It exhibits a much lower specific capacity of 43â€‰mAâ€‰hâ€‰gâˆ’1 at 0.1â€‰C. This indicates that the capacity contribution from the CoSA-HC matrix alone in the Se@CoSA-HC composite is negligible. This superior rate performance exceeds those of the reported Se-based materials, especially at high rates (Fig.Â 4b and Supplementary TableÂ 2)9,12,49,51,52.
FigureÂ 4c, d show the cycling stability of the Se@CoSA-HC cathode. After 100 cycles at current density 0.1â€‰C (Fig.Â 4c), the Se@CoSA-HC cathodes delivered a high reversible capacity of 563â€‰mAâ€‰hâ€‰gâˆ’1 with 94% capacity retention. FigureÂ 4c further shows the cycling performance of the Se@CoSA-HC at a current density of 0.5â€‰C. A large reversible specific capacity of 457â€‰mAâ€‰hâ€‰gâ€“1 was obtained after 1700 cycles. To highlight the role of atomic Co, the cycling stability of the Se@HC and Se@CoNP-HC cathodes were also tested, as presented in Supplementary Fig.Â 20, which indicates the better capacity retention of Se@CoSA-HC. Furthermore, Supplementary Fig.Â 21 shows the cycling performance of the Se@CoSA-HC cathode with a high areal loading of selenium about 5â€‰mgâ€‰cmâˆ’2 at 0.2â€‰C for 100 cycles, demonstrating stable cyclability and high capacities. When the selenium mass ratio was increased to 73% in the composite materials, the resultant Se@CoSA-HC cathodes with high Se loading delivered 242â€‰mAâ€‰h gâˆ’1 at 0.2â€‰C after 100 cycles and 220â€‰mAâ€‰hâ€‰gâˆ’1 at 0.5â€‰C after 100 cycles (Supplementary Fig.Â 22). As a highlight in Supplementary Fig.Â 23, the Se@CoSA-HC composite cathode operated for 1500 cycles at 5â€‰C with only 0.015% capacity decay per cycle from the 10th to the 1500th cycle, along with a Coulombic efficiency of nearly 100%, indicative of a quite stable prolonged cycle life. Moreover, to further demonstrate the long cycle life of Liâ€“Se battery, the long-term cycle life at a high rate of 20â€‰C (Supplementary Fig.Â 24) and 50â€‰C (Fig.Â 4d) was also tested. The Se@CoSA-HC electrodes delivered a remarkable capacity of 237â€‰mAâ€‰hâ€‰gâˆ’1 after 2500 cycles with 0.015% capacity decay per cycle at a current density of 20â€‰C and 267â€‰mAâ€‰h gâˆ’1 after 5000 cycles with 0.0067% capacity decay per cycle at a current density of 50â€‰C, and both had nearly 100% Coulombic efficiency. To the best of our knowledge, such an exceedingly good cycling stability of Liâ€“Se batteries especially at high current rates such as 50â€‰C for 5000 cycles has not been reported previously9,12,49,51,52.
In order to further clarify the electrochemical processes in Se@CoSA-HC electrodes after the first discharge-charge cycle at 3.0â€‰V, XPS measurements were carried out. Compared with bare CoSA-HC, there is a newly formed peak at 290.5â€‰eV in the C 1â€‰s XPS spectrum for Se@CoSA-HC electrodes as shown in Supplementary Fig.Â 25a. This peak was detected after the first discharge-charge cycle to 3.0â€‰V, indicating the formation of Câ€“Se bonds7,47,53. This peak can probably be ascribed to the strong interaction between the Se molecular and the carbon framework or the interaction between Se particles and a carbonyl group7. The peaks centred at 55.6 and 56.3â€‰eV in the Se 3d XPS spectrum in Supplementary Fig.Â 25b are attributed to Se 3d5/2 and Se 3d3/2, respectively. In addition, the broad peak located at 59.6â€‰eV is ascribed to Câ€“Se bonding7,47,53, which is consistent with the results from high-resolution C1s XPS spectra. Only one peak can be deconvoluted from the N 1s spectrum in Supplementary Fig.Â 25c. This is the peak located at 400.0â€‰eV, which corresponds to pyrrolic N.
To understand the excellent rate capability of Se@CoSA-HC in Liâ€“Se batteries, the electrochemical kinetics were investigated by CV measurements at various scan rates from 0.1â€‰mV/s to 0.5â€‰mV/s. As shown in Fig.Â 5a, there are two cathodic peaks (denoted as R1 and R2) related to two sequential reactions at 0.1â€‰mV/s. With increasing scan rate to 0.5â€‰mV/s, the peak R1 becomes dominant, while the intensity of the peak R2 is relatively weak. At a scan rate of 0.5â€‰mV/s, the R2 peak in the CV curve almost disappears. To further analyse the difference in kinetics, the current changes with regard to the scan rates could be analysed via as log(i)â€‰=â€‰log(a) + blog(v), where i and v are the peak current and scan rate, and a and b are derived parameters54. The reaction kinetics can be analysed via the b values. When the b-value reaches around 0.5, a diffusion-limited process is occurring during electrochemical reactions. The electrochemical behaviour is an interface-limited process and more capacitive, then b value is closer to 1, indicating faster kinetics processes55. As shown in Fig.Â 5b, the fitting curve is almost a linear relationship, and the peak R1 slope value, peak R2 slope value, peak O slope value are calculated to be 0.85, 0.85 and 0.84, respectively. This indicates that the conversion of Se into Li2Se is faster than that of Li2Se into Se and both conversions are far from being a diffusion-controlled process. FigureÂ 5c reveals that 85% of the total capacity is contributed by a capacitive process at a scan rate of 0.5â€‰mVâ€‰sâˆ’1. With an increase in the scan rate, the relative ratio of capacitive contribution to the total capacity gradually increases, as presented in Fig.Â 5d. It is believed that the capacitive behaviour is closely related to the degree of electrochemical dynamics or kinetics in battery electrodes56. Thus, the high ratio of capacitive-controlled contribution in battery electrodes is highly beneficial for fast transport of lithium ions, which would lead to superior electrochemical performance of rate capability and long-term cyclability.
To further evaluate the electrochemical processes and kinetics of Liâ€“Se reactions when using Se@CoSA-HC composite cathodes, EIS spectra were collected as a function of the state of discharge/charge. Supplementary Fig.Â 26a illustrates the typical discharge/charge profile of a Liâ€“Se battery using a Se@CoSA-HC composite cathode at 0.1â€‰C. The EIS spectra of the Se@CoSA-HC electrode at various depths (marked in Supplementary Fig.Â 26a) during the discharge/charge process are shown in Supplementary Fig.Â 26b. All the results exhibit two depressed or overlapping semicircles followed by a sloping line. For Liâ€“Se batteries based on ether-based electrolyte, the first semicircle and second semicircle in EIS spectra for Liâ€“Se batteries could be ascribed to the charge transfer resistance of Se/C electrode and accumulation of interfacial layer on cathode surface, respectively57. Combined with the XPS results in Supplementary Fig.Â 25, it is proposed that a stable layer formed at the surface of the Se@CoSA-HC particles during the first lithiation process. Despite the solubility of Se, lithium polyselenides and Li2Se in ether-based electrolytes11, side reactions have been restricted when using Se@CoSA-HC composites and the formed LixSe may be protected from further reaction by the stable layer. The equivalent circuit for fitting is shown in the inset of Supplementary Fig.Â 26c. In the equivalent circuits, R0 represents the impedance that is mainly derived from the resistance of the electrolyte, R1 is the charge transfer resistance at the conductive agent interface, and R2 is the resistance of interfacial layer57. CPE1 represents double-layer capacitance (Cdl), while CPE2 (Constant phase element) describes the space charge capacitance of the layer. W0 is the Warburg impedance corresponding to the polyselenide diffusion processes. The resistance values (R1 and R2) obtained from Supplementary Fig.Â S26b are summarized in Supplementary Fig.Â S26c. Supplementary Fig.Â S27 shows the typical Nyquist plots collected at different discharge-charge state and their fitted curves with two depressed or overlapping semicircles. The obtained resistance values are presented in Supplementary TableÂ 3. The R1 value is relatively stable throughout the whole cycle process, indicating the excellent charge transfer capability. The variation of the charge transfer resistance is attributed to the transformation of crown-like Se8 to Li2Sex during the discharge process, then maintaining amorphous chain-like Se molecules during the charge process. The formation of a stable layer and the reversible Se transformations are responsible for the superior electrochemical performance. The cell with the Se@CoSA-HC cathode after 1700 cycles was dissembled and the retrieved Se@CoSA-HC cathode materials were characterized by TEM. As shown in Supplementary Fig.Â 28a, it is observed that the morphology of the Se@CoSA-HC cathode has been well preserved, indicating the superior stability of the cathode structure. From STEM element mapping images in Supplementary Fig.Â 28c, d, energy-dispersive X-ray spectroscopy (EDS) images were obtained to identify the elemental distribution of cobalt and selenium, which are homogeneously distributed in the whole carbon framework. As shown in Supplementary Fig.Â 28b, a homogeneous thin layer with a thickness of ~20â€‰nm has been identified on the surface of the electrode after long-term cycling. These characterizations confirmed that the composition of the layer could be Li2Se2/Li2Se, which could enhance the cycling performance of the Liâ€“Se batteries. This result is consistent with the previously reported literature58. Furthermore, EIS spectra were recorded at open-circuit voltage before cycling and after 1st cycle, 2nd cycle, 5th cycle, 10th cycle and 50th cycle at 0.1â€‰C (Supplementary Fig.Â 29), 2â€‰C (Supplementary Fig.Â 30), 5â€‰C (Supplementary Fig.Â 31) and 20â€‰C (Supplementary Fig.Â 32). As shown in Supplementary Fig.Â 29â€“32, the charge-transfer resistance is very stable after charge/discharge at various current densities. These results further strengthen our claim that Se@CoSA-HC materials possess superior rate performance. Furthermore, it is also found that the diameter of the semicircle at 5â€‰C and 20â€‰C decreases with time, which may be attributed to the further conversion reaction with lithium ions and the stabilization of the layer with increasing cycling time59.
Owing to the electrocatalytic effect resulting from single cobalt atoms within CoSA-HC particles, the Se@CoSA-HC electrode shows the lowest overpotentials in charge profiles among the three as-synthesized materials (Fig.Â 6a). This is in accordance with the smaller value of voltage change Î”V (lowest voltage hysteresis) of test cells by using Se@CoSA-HC to produce discharge-charge curves under a current density of 0.1â€‰C compared with Se@HC and Se@CoNP-HC, as shown in Fig.Â 6b. Galvanostatic intermittent titration technique (GITT) testing and EIS spectra were employed to determine the catalyst effect of the single cobalt atoms in Liâ€“Se batteries. The GITT plots of Se@CoSA-HC and Se@HC (Supplementary Fig.Â 33a) and the diffusion coefficients DLi+ (Supplementary Fig.Â 33b) of the Se@CoSA-HC calculated from GITT plots were estimated to be 1.02â€‰Ã—â€‰10âˆ’13â€‰~â€‰1.7â€‰Ã—â€‰10âˆ’13 (Supplementary NoteÂ 1), which is almost ten times that of the Se@HC cathode (1.03â€‰Ã—â€‰10âˆ’14â€‰~â€‰8.35â€‰Ã—â€‰10âˆ’14 (Supplementary NoteÂ 1)) material, indicating the better electrochemical kinetics for Se@CoSA-HC cathode materials in Liâ€“Se batteries. According to the Butlerâ€“Volmer equation and EIS spectra (Supplementary Fig.Â 34), the exchange current density i0 of Se@CoSA-HC and Se@HC cathode is 1.14 and 0.76â€‰mAâ€‰cmâˆ’2 (Supplementary NoteÂ 2), respectively. To confirm the catalytic role and stability of single cobalt atoms after long-term cycling, the cell after cycled for 1700 cycles was dissembled and the Se@CoSA-HC cathode was retrieved and characterized. From high angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) images in Supplementary Fig.Â 35, high-density bright dots (highlighted by red circles) have been clearly detected, revealing that single cobalt atoms in the Se@CoSA-HC cathode are stable during long-term cycling. All these key kinetic parameters confirm the catalyst effect of the single cobalt atoms.
To further confirm the catalytic activity of single Co atoms after cycling, we conducted a visual observation on the suppression for the formation of lithium polyselenides during the cycling process via single Co atom catalysis. The Se@CoSA-HC cathode after cycled for 1700 cycles and the cycled bare Se@HC cathode after cycled 0.5â€‰C for 100 cycles were used for the visual observation. As shown in Supplementary Fig.Â 36a, H-type cell was used in which the cycled Se@CoSA-HC cathode and bare Se@HC cathode were used as the cathode. The electrolyte in the cell with cycled bare Se@HC cathode (Supplementary Fig.Â 36b, c) changed from colourless to yellow after the first discharge process, implying there was a dissolution of polyselenides in electrolyte60. It is clear that by using bare Se@HC cathode without single cobalt catalysts, lithium polyselenides were detached from the cathode and dissolved in the electrolyte, indicating that bare Se@HC cathode owns poor immobilization of polyselenides.
In contrast, the colour change of the electrolyte for the cycled Se@CoSA-HC cathode was not observed during three cycles. Supplementary Fig.Â 36d, e show that the cell cycled Se@CoSA-HC cathode presents transparent electrolyte without colour change, indicating the atomic cobalt electrocatalyst could effectively alleviate the dissolution of polyselenides, maximize polyselenides immobilization, electro-catalyse the transformation from polyselenides to Li2Se. Therefore, this visual experiment verifies that the single cobalt atoms play an important role that the transformation from polyselenides to Li2Se have been electro-catalysed, polyselenides immobilization has been maximized and the dissolution of polyselenides has been inhibited. This visual observation using cycled electrodes clearly confirmed that single atom cobalt catalysts remain active after long cycling.
To further understand the enhancement of reaction kinetics of charge/discharge of the Se@CoSA-HC cathodes, first-principles calculations were performed to investigate the different possible reactions of lithium polyselenides on nitrogen-doped carbon support (NC) as a reference and atomic Co/nitrogen-doped carbon supports (Coâ€“NC). As shown in Supplementary Fig.Â 37, two models of nitrogen-doped carbon without and with Co atoms were considered in our simulation. The reversible overall reaction for the formation of Li2Se originating from Se8 and Li was considered12. During discharge, the first step involves the reduction of Se8 and the generation of Li2Se8, followed by further reduction and disproportionation with the formation of three intermediate lithium polyselenides, namely, Li2Se6, Li2Se4, and Li2Se2, achieving the formation of Li2Se as the final product12. The Gibbs free energies were calculated for the above reactions on both NC and Coâ€“NC supports (Supplementary TableÂ 4). The optimized structures of the intermediates and their Gibbs free energy profiles are displayed in Fig.Â 6c. It can be observed that the transformations from Se8 to Li2Se6 are exothermic and the following three steps involving the conversion of Li2Se4, Li2Se2, and Li2Se are endothermic. The largest positive Gibbs free energy can be found in the conversion process from Li2Se2 to Li2Se, revealing its role as the rate-determining step in the whole discharge process. The Gibbs free energy of Coâ€“NC support (0.85â€‰eV) is much lower than that of NC support (0.96â€‰eV), indicating that the reduction of Se is thermodynamically more favourable on Coâ€“NC than on NC support. In the charging process, the transformation of Li2Se is the first step7. Through the climbing-image nudged elastic band method, the transformation energy and barrier of Li2Se were calculated to evaluate the delithiation reaction kinetics from Li2Se to selenium on the surfaces of Coâ€“NC and NC supports. FigureÂ 6d, e show the energy profiles for the transformation processes on both Coâ€“NC and NC supports. The calculated energy barriers of Li2Se transformation of Coâ€“NC supports (1.82â€‰eV) is smaller than that of NC (2.04â€‰eV), revealing that atomic cobalt nanoparticles are serving as active sites to enhance the phase transformation of Li2Se and the Se utilization in Liâ€“Se batteries.
In addition, from the DFT calculation results shown in Supplementary TableÂ 5, the Liâ€“Se bond length on the surface of Coâ€“NC support is elongated, revealing the interaction between Li atom and the rest of the molecule is weakened. This leads to the easier delithiation for lithium polyselenides, confirming the role of single Co atom catalysts. Therefore, we proposed a mechanism that single Co atoms can quickly catalyse the transformation from Li2Se2 into Li2Se during the discharging process and the transformation of Li2Se during the charging process. To demonstrate the mechanism, the schematic illustration of electrode reaction mechanisms for the Se@CoSA-HC cathodes is shown in Supplementary Fig.Â 38. The confined polyselenides within the carbon framework could be fully reduced into Li2Se by single atom Co catalysts, leading to high Se utilization. Therefore, the atomic cobalt on the Coâ€“NC support could effectively alleviate the dissolution of polyselenides, electro-catalyse the transformation from polyselenides to Li2Se and minimize the reaction energy barriers, leading to the full utilization of selenium, superior cycling ability, and reversible capability.
Until now, single-atom catalysts were achieved through various strategies, including wet impregnation and coprecipitation methods61, atomic layer deposition62, pyrolysis63 and photodeposition64. Among them, pyrolysis is a facile method to construct single-atom catalysts through thermal decomposition of suitable precursors. Herein, we have developed a simple synthesis method to achieve hollow structured particles with isolated, positively charged and highly dispersed single Co atoms through tuning the cobalt and zinc contents from bimetallic ZnCo-ZIFs precursors. The aberration-corrected HAADF-STEM image and EXAFS data provide strong evidence that single Co atoms are high dispersed within the CoSA-HC particles. It is also revealed from XANES results that the single Co atoms are positively charged and the EXAFS data show that isolated single Co atoms can be atomically anchored into the carbon frameworks through the formation of Coâ€“N3 and Coâ€“N4 coordination moieties within CoSA-HC particles. By using CoSA-HC particles, more accessible storage sites and larger electrode/electrolyte contact area are provided, and mass/charge transportation lengths are shortened through the formation of the hollow structure. The volume expansion during lithiation can be inhibited by the large internal void spaces.
Based on the density functional theory (DFT) calculations, it indicates that the reaction rate from the reduction of Li2Se2 into Li2Se is noticeably increased during the discharge process by using a single Co atom catalysts. From the charging process data, a smaller value than the calculated energy barriers for Li2Se transformation through the single Co atom catalysts can be obtained. It is proposed that the mechanism of single Co atom enhancement of Liâ€“Se batteries is that single Co atoms can quickly catalyse the transformation from Li2Se2 into Li2Se during the discharging process and the transformation of Li2Se during the charging process. Therefore, the atomic cobalt plays the key role in the alleviation of polyselenide dissolution, maximation of polyselenides immobilization and activation via strong electrocatalytic behaviour, achieving the best cycling performance in the field of Liâ€“Se batteries. Therefore, the Se@CoSA-HC composite is a promising candidate for long cycle life and high-power lithiumâ€“selenium batteries.
In summary, we developed a facile approach for synthesizing coreâ€“shell structured PS@ZIF materials which can be further converted into atomic Co electrocatalyst/nitrogen-doped hollow porous carbon (CoSA-HC) particles through one-step pyrolysis. To highlight the importance of CoSA-HC particles for energy applications, CoSA-HC particles were used for Liâ€“Se batteries. More specifically, Se@CoSA-HC cathodes delivered an excellent discharge capacity of 564â€‰mAâ€‰hâ€‰gâˆ’1 after 100 cycles at a current density of 0.1â€‰C and a superior rate capability of 385â€‰mAâ€‰hâ€‰gâˆ’1 and 311â€‰mAâ€‰hâ€‰gâˆ’1 at a current density of 20â€‰C and 50â€‰C. In addition, they show a superior reversible capacity of 457â€‰mAâ€‰hâ€‰gâˆ’1 at a current density of 0.5â€‰C after 1700 cycles with only 0.011% capacity decline per cycle and 267â€‰mAâ€‰hâ€‰gâˆ’1 after 5000 cycles at 50â€‰C with a 0.0067% capacity decay per cycle with Coulombic efficiency nearly 100%. These distinctive features of the HC particles are concluded as follows: (i) the atomic cobalt electrocatalyst could effectively alleviate the dissolution of polyselenides, electro-catalyse the transformation from polyselenides to Li2Se and minimize the adsorption energy barriers; (ii) the hollow structures provide more accessible selenium storage sites, larger electrode/electrolyte contact area, shortened mass/charge transport lengths; (iii) the large internal void spaces accommodate volume expansion during lithiation; (iv) the conductive carbon materials enhance the electrode conductivity and effectively confine the soluble polyselenides. This work provides an efficient route for the preparation of single-atom materials and paves a new strategy for developing high-power electrochemical energy storage devices.
The following provides information on chemicals used in this work: methanol (99%), styrene, potassium persulfate (KPS), polyvinylpyrrolidone (PVP, MW~50,000), ethanol (95â€“100%), zinc nitrate hexahydrate (Zn(NO3)2Â·6H2O), cobalt nitrate hexahydrate (Co(NO3)2Â·6H2O), 2-methylimidazole were purchased from Sigma-Aldrich and used as received without any further purification. Washing was achieved with ultrapure water and reagent grade ethanol where required. Ultrapure water was used for solution preparations.
The styrene was firstly treated by passing through a column containing neutral alumina to remove polymerization inhibitor and then subjected to vacuum distillation to obtain a purified styrene monomer. 1.1â€‰g of PVP was dissolved in 100â€‰mL deionized water, and then 13â€‰mL of the above treated styrene was added to the solution. After stirring for 30â€‰mins, an aqueous 20â€‰mL solution containing 0.3â€‰g of potassium persulfate (KPS) was then added to the above mixture. Then the mixture was heated to 70â€‰Â°C and remained at this temperature for 24â€‰h.
In a typical procedure, 0.15â€‰g of as-prepared PS spheres were fully ground and then dispersed into 90â€‰mL of methanol containing 1â€‰g of PVP. After ultrasonic dispersion and vigorous agitation for 3â€‰h, 2.125â€‰g of Zn(NO3)2Â·6H2O and 0.104â€‰g of Co(NO3)2Â·6H2O (the molar ratio of Zn2+/Co2+ was 20 were subsequently dissolved into the mixed solution and stirred for another 0.5â€‰h. Then, 4.926â€‰g of 2-methylimidazole dissolved in 90â€‰mL of methanol was quickly added into the above solution followed by vigorous stirring for 4â€‰h. Finally, the PS@ZIF-1 precursor were collected by centrifugation, washed with methanol for several times, and dried at 60â€‰Â°C overnight. Similarly, for the synthesis of PS@ZIF-2 and PS@ZIF-3, the procedures were carried out except that the Zn/Co molar ratio was 1:0 and 17:4, respectively.
PS@ZIF particles (PS@ZIF-1, PS@ZIF-2 and PS@ZIF-3) were carbonized in flowing N2 in a tube furnace using a heating rate of 5â€‰Â°Câ€‰minâˆ’1 up to 700â€‰Â°C and dwelling for 5â€‰h. After naturally cooling down to room temperature, the products were collected and marked as CoSA-HC, HC and CoNP-HC particles according to the Zn/Co molar ratio in the precursors, respectively.
Se powder (99.99%, Sigma-Aldrich) and the as-prepared CoSA-HC, HC, CoNP-HC particles with a weight ratio of 1:1 were mixed. Subsequently, the mixture was heated at 300â€‰Â°C for 12â€‰h with heating rate of 5â€‰Â°Câ€‰minâˆ’1 in a tubular furnace under argon atmosphere to achieve selenium carbon composites. The as-prepared materials were named Se@CoSA-HC, Se@HC and Se@CoNP-HC based on the different carbon precursors.
To achieve high selenium mass ratio in selenium carbon composite, Se powder and the as-prepared CoSA-HC particles with a weight ratio of 3:1 were mixed, followed by the same heat treatment as described above.
Synthesis of MnO2 nanowires: in a typical synthesis, 50â€‰mg of polyvinylpyrrolidone (PVP) and 40â€‰mL of 0.015â€‰M KMnO4 aqueous solution were mixed with magnetic stirring, and then the mixture was transferred into a 50â€‰mL Teflon-lined stainless autoclave. The autoclave was sealed and put in an electronic oven at 160â€‰Â°C for 9â€‰h and then naturally cooled down to room temperature. The precipitates were collected by filtration, washed with deionized water and absolute ethanol several times before drying at 60â€‰Â°C overnight.
Synthesis of reduced graphene oxide (rGO): Graphene oxide (3â€‰mg/ml) solution were dried via lyophilization and then was ground into powder. Reduced graphene oxide can be obtained through hydrothermal treatment of graphene oxide (0.4â€‰g) with 400â€‰ÂµL hydrazine hydrate at 95â€‰Â°C for 24â€‰h.
Synthesis of the MnO2@ZIF: in a typical procedure, 0.07â€‰g of as-prepared MnO2 nanowires were fully ground and then dispersed into 90â€‰mL of methanol containing 1â€‰g of PVP (K-30). After ultrasonic dispersion and vigorous agitation for 3â€‰h, 2.125â€‰g of Zn(NO3)2Â·6H2O and 0.104â€‰g of Co(NO3)2Â·6H2O (the molar ratio of Zn2+/Co2+ was 20 were subsequently dissolved into the mixed solution and stirred for another 0.5â€‰h. Then, 4.926â€‰g of 2-methylimidazole dissolved in 90â€‰mL of methanol was quickly added into the above solution followed by vigorous stirring for 4â€‰h. Finally, the MnO2@ZIF precursors were collected by centrifugation, washed with methanol for several times, and dried at 60â€‰Â°C overnight.
Synthesis of rGO@ZIF: in a typical procedure, 0.07â€‰g of as-prepared rGO were fully ground and then dispersed into 90â€‰mL of methanol containing 1â€‰g of PVP (K-30). After ultrasonic dispersion and vigorous agitation for 3â€‰h, 2.125â€‰g of Zn(NO3)2Â·6H2O and 0.104â€‰g of Co(NO3)2Â·6H2O (the molar ratio of Zn2+/Co2+ was 20 were subsequently dissolved into the mixed solution and stirred for another 0.5â€‰h. Then, 4.926â€‰g of 2-methylimidazole dissolved in 90â€‰mL of methanol was quickly added into the above solution followed by vigorous stirring for 4â€‰h. Finally, the rGO@ZIF precursors were collected by centrifugation, washed with methanol for several times, and dried at 60â€‰Â°C overnight.
Electrodes were prepared by mixing the selenium carbon composite, carbon black and polyvinylidene fluoride at a weight ratio of 8:1:1 in N-methyl-2-pyrrolidone solvent, The slurry was pasted onto aluminium foil and dried in a vacuum oven at 60â€‰Â°C for 12â€‰h. CR2032 coin cells were assembled in an argon-filled glove box (Mbraun, Unilab, Germany) and then used for electrochemical evaluation. The electrolyte contained 1â€‰wt% lithium nitrate (LiNO3) in 1,3-dioxolane and 1,2-dimethoxyethane (volume ratio 1:1). Porous polypropylene (Celgard 2300TM) was used for the separator membranes. About 20â€‰Î¼L electrolyte was added for each coin cell. The areal loading of selenium in the cathode is 0.8â€‰mgâ€‰cmâˆ’2. The cells were discharged and charged galvanostatically in the fixed voltage range 1.00â€“3.00â€‰V with current rates of 0.1 â€“50â€‰C rate using a NEWARE battery tester. Cyclic voltammetry (CV) and electrochemical impedance spectroscopy (EIS) were performed on a Biologic VMP3 electrochemical station. CV responses were also tested in voltage range 1.00â€‰V to 3.00â€‰V (vs. Li+/Li) at scan rate from 0.1 to 0.5â€‰mVâ€‰sâˆ’1. EIS (100â€‰kHz to 10â€‰mHz with an amplitude of 5â€‰mV). Analysis was then applied to evaluate electrochemical behaviors of the electrode.
The morphology and chemical composition of the as-prepared samples were observed by field emission scanning electron microscopy (FESEM, Zeiss Supra 55VP), transmission electron microscopy (TEM, Tecnai G2 F30 S-TWIN). X-ray diffraction (XRD) measurements were carried out by using a scanning step of 0.04Â° per second in the 2Î¸ range from 10Â° to 80Â° (Bruker D8 Discovery XRD). X-ray photoelectron spectroscopy (XPS) measurements were performed on an ESCALAB250Xi (Thermo Scientific, UK) equipped with mono-chromated Al K alpha (energy 1486.68â€‰eV). The BET specific surface area and single-point pore volume were obtained from nitrogen adsorption isotherms measured at âˆ’196â€‰Â°C using a nitrogen sorption instrument (ASAP 2460 Micropore Physisorption Analyzer). Prior to nitrogen adsorption measurements, the samples were degassed at 250â€‰Â°C overnight. Raman spectra were obtained from a Renishaw inVia Raman spectrometer system (Gloucestershire, UK) equipped with a Leica DMLB microscope (Wetzlar, Germany) and a Renishaw He-Ne laser source producing 17â€‰mW at 633â€‰nm. Thermogravimetric analysis (TGA) was performed with a 2960 SDT system. X-ray absorption fine structure (XAFS) measurements of the Co K-edge were performed at the 1W1B beamline of the Beijing Synchrotron Radiation Facility (BSRF) in transmission mode. The X-ray was monochromatized using a double-crystal Si (111) monochromator and the energy was calibrated by a cobalt metal foil for Co K-edge. Cobalt foil and CoPc were used as the reference substance. The XAFS data was analysed using the WinXAS3.1 program. Cobalt content in HC particles was determined by Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES).
All calculations were carried out by using the projector augmented wave method in the framework of DFT65, as implemented in the Vienna ab-initio Simulation Package (VASP). The generalized gradient approximation (GGA) and Perdewâ€“Burkeâ€“Ernzerhof (PBE) exchange functional65 were used. Structural relaxation calculations were performed by using the spin-polarized GGA method66. A supercell of graphene containing 4â€‰Ã—â€‰2(sqrt 3) unit cells was used to model N-G and Coâ€“N-G systems. After convergence tests, the plane-wave energy cutoff was set to 500â€‰eV, and the Monkhorstâ€“Pack method67 with 2â€‰Ã—â€‰3â€‰Ã—â€‰1 supercells was employed for the Brillouin zone sampling of N-G and Coâˆ’N-G systems. The convergence criterions of energy and force calculations were set to 10âˆ’5 eV/atom and 0.01â€‰eVâ€‰Ã…âˆ’1, respectively. The Gibbs free energies for electrochemical reduction of molecular Se8 and polyselenide Li2Sen on the N-G and Coâ€“N-G systems were calculated by the DFT energy difference between two different reduction steps68,69. The DFT-D2 empirical correction method was employed to describe van der Waals interactions. The barriers for Li2Se transformation on N/G and Coâ€“N/G were calculated with the nudged elastic band (NEB) method to evaluate its delithiation reaction kinetics.
The data that support the findings of this work are available from the corresponding author upon reasonable request.
Dunn, B., Kamath, H. & Tarascon, J.-M. Electrical energy storage for the grid: a battery of choices. Science 334, 928â€“935 (2011).
Posada, J. O. G. et al. Aqueous batteries as grid scale energy storage solutions. Renew. Sustain. Energy Rev. 68, 1174â€“1182 (2017).
Choi, J. W. & Aurbach, D. Promise and reality of post-lithium-ion batteries with high energy densities. Nat. Rev. Mater. 1, 16013 (2016).
Bruce, P. G., Freunberger, S. A., Hardwick, L. J. & Tarascon, J.-M. Liâ€“O2 and Liâ€“S batteries with high energy storage. Nat. Mater. 11, 19 (2011).
Yang, X., Wang, H., Yu, D. Y. W. & Rogach, A. L. Vacuum calcination induced conversion of selenium/carbon wires to tubes for high-performance sodiumâ€“selenium batteries. Adv. Funct. Mater. 28, 1706609 (2018).
Yao, Y. et al. CNT interwoven nitrogen and oxygen dual-doped porous carbon nanosheets as free-standing electrodes for high-performance Na-Se and K-Se flexible batteries. Adv. Mater. 30, 1805234 (2018).
Li, Z., Yuan, L., Yi, Z., Liu, Y. & Huang, Y. Confined selenium within porous carbon nanospheres as cathode for advanced Liâ€“Se batteries. Nano Energy 9, 229â€“236 (2014).
Yang, C.-P. et al. An advanced seleniumâ€“carbon cathode for rechargeable lithiumâ€“selenium batteries. Angew. Chem. Int. Ed. 52, 8363â€“8367 (2013).
Lee, J. T. et al. Micro- and mesoporous carbide-derived carbonâ€“selenium cathodes for high-performance lithium selenium batteries. Adv. Energy Mater. 5, 1400981 (2015).
Cui, Y. et al. (De)lithiation mechanism of Li/SeSx (xâ€‰=â€‰0â€“7) batteries determined by in situ synchrotron X-ray diffraction and X-ray absorption spectroscopy. J. Am. Chem. Soc. 135, 8047â€“8056 (2013).
Abouimrane, A. et al. A new class of lithium and sodium rechargeable batteries based on selenium and seleniumâ€“sulfur as a positive electrode. J. Am. Chem. Soc. 134, 4505â€“4508 (2012).
Luo, C. et al. Selenium@mesoporous carbon composite with superior lithium and sodium storage capacity. ACS Nano 7, 8003â€“8010 (2013).
Xin, S. et al. The electrochemistry with lithium versus sodium of selenium confined to slit micropores in carbon. Nano Lett. 16, 4560â€“4568 (2016).
Liu, T. et al. Selenium embedded in metalâ€“organic framework derived hollow hierarchical porous carbon spheres for advanced lithiumâ€“selenium batteries. ACS Appl. Mater. Interfaces 8, 16063â€“16070 (2016).
Tian, H. et al. Tunable porous carbon spheres for high-performance rechargeable batteries. J. Mater. Chem. A 6, 12816â€“12841 (2018).
Zeng, L. et al. A flexible porous carbon nanofibers-selenium cathode with superior electrochemical performance for both Li-Se and Na-Se batteries. Adv. Energy Mater. 5, 1401377 (2015).
Qu, Y. et al. Confining selenium in nitrogen-containing hierarchical porous carbon for high-rate rechargeable lithiumâ€“selenium batteries. J. Mater. Chem. A 2, 12255â€“12261 (2014).
Zhang, J. et al. Selenium/interconnected porous hollow carbon bubbles composites as the cathodes of Li-Se batteries with high performance. Nanoscale 6, 12952â€“12957 (2014).
Wang, A., Li, J. & Zhang, T. Heterogeneous single-atom catalysis. Nat. Rev. Chem. 2, 65â€“81 (2018).
Guo, X. et al. Direct, nonoxidative conversion of methane to ethylene, aromatics, and hydrogen. Science 344, 616â€“619 (2014).
Zhang, J. et al. Single platinum atoms immobilized on an MXene as an efficient catalyst for the hydrogen evolution reaction. Nat. Catal. 1, 985â€“992 (2018).
Yang, L., Shi, L., Wang, D., Lv, Y. & Cao, D. Single-atom cobalt electrocatalysts for foldable solid-state Zn-air battery. Nano Energy 50, 691â€“698 (2018).
Zang, W. et al. Single Co atoms anchored in porous N-doped carbon for efficient zincâˆ’air battery cathodes. ACS Catal. 8, 8961â€“8969 (2018).
Wang, J. et al. Single-atom catalyst boosts electrochemical conversion reactions in batteries. Energy Storage Mater. 18, 246â€“252 (2019).
Zhang, B.-W. et al. Atomic cobalt as an efficient electrocatalyst in sulfur cathodes for superior room-temperature sodium-sulfur batteries. Nat. Commun. 9, 4082 (2018).
Du, Z. et al. Cobalt in nitrogen-doped graphene as single-atom catalyst for high-sulfur content lithium-sulfur batteries. J. Am. Chem. Soc. 141, 3977â€“3985 (2019).
Qu, Y. et al. Direct transformation of bulk copper into copper single sites via emitting and trapping of atoms. Nat. Catal. 1, 781â€“786 (2018).
Wang, X. et al. Regulation of coordination number over single Co sites: triggering the efficient electroreduction of CO2. Angew. Chem. Int. Ed. 57, 1944â€“1948 (2018).
Li, J. et al. Atomically dispersed manganese catalysts for oxygen reduction in proton-exchange membrane fuel cells. Nat. Catal. 1, 935â€“945 (2018).
Chong, L. et al. Ultralow-loading platinum-cobalt fuel cell catalysts derived from imidazolate frameworks. Science 362, 1276â€“1281 (2018).
Zhong, H. et al. ZIF-8 derived graphene-based nitrogen-doped porous carbon sheets as highly efficient and durable oxygen reduction electrocatalysts. Angew. Chem. Int. Ed. 53, 14235â€“14239 (2014).
Chen, Y. et al. From bimetallic metal-organic framework to porous carbon: high surface area and multicomponent active dopants for excellent electrocatalysis. Adv. Mater. 27, 5010â€“5016 (2015).
Liu, D., Wan, J., Pang, G. & Tang, Z. Hollow metalâ€“organic-framework micro/nanostructures and their derivatives: emerging multifunctional materials. Adv. Mater. 31, 1803291 (2018).
Nai, J. & Lou, X. W. Hollow structures based on prussian blue and its analogs for electrochemical energy storage and conversion. Adv. Mater. 31, 1706825 (2018).
Wang, J., Cui, Y. & Wang, D. Design of hollow nanostructures for energy storage, conversion and production. Adv. Mater. 31, 1801993 (2018).
Tian, H. et al. Enhanced hydrogenation performance over hollow structured Co-CoOx@N-C capsules. Adv. Sci. 6, 1900807 (2019).
Tian, H., Liang, J. & Liu, J. Nanoengineering carbon spheres as nanoreactors for sustainable energy applications. Adv. Mater. 31, 1903886 (2019).
Yin, P. et al. Single cobalt atoms with precise N-coordination as superior oxygen reduction reaction catalysts. Angew. Chem. Int. Ed. 55, 10800â€“10805 (2016).
Yi, Z. et al. High-performance lithiumâ€“selenium batteries promoted by heteroatom-doped microporous carbon. J. Mater. Chem. A 3, 3059â€“3065 (2015).
Li, Z. & Yin, L. MOF-derived, N-doped, hierarchically porous carbon sponges as immobilizers to confine selenium as cathodes for Liâ€“Se batteries with superior storage capacity and perfect cycling stability. Nanoscale 7, 9597â€“9606 (2015).
Tian, H. et al. The development of yolkâ€“shell-structured Pd&ZnO@carbon submicroreactors with high selectivity and stability. Adv. Funct. Mater. 28, 1801737 (2018).
Wang, Q. et al. NiFe layered double hydroxide nanoparticles on Co,N-codoped carbon nanoframes as efficient bifunctional catalysts for rechargeable zincâ€“air batteries. Adv. Energy Mater. 7, 1700467 (2017).
Tian, H. et al. Metal-organic-framework-derived formation of Coâ€“N-doped carbon materials for efficient oxygen reduction reaction. J. Energy Chem. 40, 137â€“143 (2020).
Lu, Z. et al. An isolated zincâ€“cobalt atomic pair for highly active and durable oxygen reduction. Angew. Chem. Int. Ed. 58, 2622â€“2626 (2019).
Cheng, Q. et al. Single cobalt atom and N codoped carbon nanofibers as highly durable electrocatalyst for oxygen reduction reaction. ACS Catal. 7, 6864â€“6871 (2017).
Han, X. et al. Generation of nanoparticle, atomic-cluster, and single-atom cobalt catalysts from zeolitic imidazole frameworks by spatial isolation and their use in zincâ€“air batteries. Angew. Chem. Int. Ed. 131, 5413â€“5418 (2019).
Ye, H., Yin, Y.-X., Zhang, S.-F. & Guo, Y.-G. Advanced Seâ€“C nanocomposites: a bifunctional electrode material for both Liâ€“Se and Li-ion batteries. J. Mater. Chem. A 2, 13293â€“13298 (2014).
Li, X. et al. A new salt-baked approach for confining selenium in metal complex-derived porous carbon with superior lithium storage properties. Adv. Funct. Mater. 25, 5229â€“5238 (2015).
Liu, Y. et al. A selenium-confined microporous carbon cathode for ultrastable lithiumâ€“selenium batteries. J. Mater. Chem. A 2, 17735â€“17739 (2014).
Zeng, L. et al. Green synthesis of a Se/HPCFâ€“rGO composite for Liâ€“Se batteries with excellent long-term cycling performance. J. Mater. Chem. A 5, 22997â€“23005 (2017).
Zhou, Y., Li, Z. & Lu, Y.-C. A stable lithiumâ€“selenium interface via solid/liquid hybrid electrolytes: blocking polyselenides and suppressing lithium dendrite. Nano Energy 39, 554â€“561 (2017).
Han, K. et al. A free-standing and ultralong-life lithium-selenium battery cathode enabled by 3D mesoporous carbon/graphene hierarchical architecture. Adv. Funct. Mater. 25, 455â€“463 (2015).
Zhang, H., Yu, F., Kang, W. & Shen, Q. Encapsulating selenium into macro-/micro-porous biochar-based framework for high-performance lithium-selenium batteries. Carbon 95, 354â€“363 (2015).
Augustyn, V. et al. High-rate electrochemical energy storage through Li+ intercalation pseudocapacitance. Nat. Mater. 12, 518 (2013).
Ding, J., Zhou, H., Zhang, H., Tong, L. & Mitlin, D. Selenium impregnated monolithic carbons as free-standing cathodes for high volumetric energy lithium and sodium metal batteries. Adv. Energy Mater. 8, 1701918 (2018).
Chen, Z. et al. General synthesis of dual carbon-confined metal sulfides quantum dots toward high-performance anodes for sodium-ion batteries. Adv. Funct. Mater. 27, 1702046 (2017).
Wang, X. et al. New insight into the confinement effect of microporous carbon in Li/Se battery chemistry: a cathode with enhanced conductivity. Small 16, 2000266 (2020).
Liu, F. et al. A mixed lithium-ion conductive Li2S/Li2Se protection layer for stable lithium metal anode. Adv. Funct. Mater. 30, 2001607 (2020).
Yang, L. et al. Hierarchical MoO2/Mo2C/C hybrid nanowires as high-rate and long-life anodes for lithium-ion batteries. ACS Appl. Mater. Interfaces 8, 19987â€“19993 (2016).
Gu, X. et al. Highly reversible Liâ€“Se batteries with ultra-lightweight N, S-codoped graphene blocking layer. Nano-Micro Lett. 10, 59 (2018).
Lang, R. et al. Non defect-stabilized thermally stable single-atom catalyst. Nat. Commun. 10, 234 (2019).
Piernavieja-Hermida, M. et al. Towards ALD thin film stabilized single-atom Pd1 catalysts. Nanoscale 8, 15348â€“15356 (2016).
Han, A. et al. A polymer encapsulation strategy to synthesize porous nitrogen-doped carbon-nanosphere-supported metal isolated-single-atomic-site catalysts. Adv. Mater. 30, 1706508 (2018).
Liu, P. et al. Photochemical route for synthesizing atomically dispersed palladium catalysts. Science 352, 797â€“800 (2016).
Kohn, W. & Sham, L. J. Self-consistent equations including exchange and correlation effects. Phys. Rev. 140, A1133â€“A1138 (1965).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865â€“3868 (1996).
Monkhorst, H. J. & Pack, J. D. Special points for brillouin-zone integrations. Phys. Rev. B 13, 5188â€“5192 (1976).
Xu, Z. et al. First principles study of adsorption and oxidation mechanism of elemental mercury by HCl over MoS2 (100) surface. Chem. Eng. J. 308, 1225â€“1232 (2017).
Lv, X., Xu, Z., Li, J., Chen, J. & Liu, Q. Investigation of fluorine adsorption on nitrogen doped MgAl2O4 surface by first-principles. Appl. Surf. Sci. 376, 97â€“104 (2016).
Prof. Hao Liu thanks the support of the ARC Future Fellowship (FT180100705) and Discovery Project (DP180102297). This work was partially supported by the DNL Cooperation Fund, CAS (DNL180402).
The authors declare no competing interests.
Peer review information Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisherâ€™s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.