Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Interaction between functionalized graphene and sulfur compounds in a lithium–sulfur battery – a density functional theory investigation

Kimal Chandula Wasalathilake , Md Roknuzzaman , Kostya (Ken) Ostrikov , Godwin A. Ayoko and Cheng Yan *
School of Chemistry, Physics and Mechanical Engineering, Faculty of Science and Engineering, Queensland University of Technology (QUT), Brisbane, QLD 4001, Australia. E-mail: c2.yan@qut.edu.au

Received 21st October 2017 , Accepted 27th December 2017

First published on 9th January 2018


Abstract

Lithium–sulfur (Li–S) batteries are emerging as one of the promising candidates for next generation rechargeable batteries. However, dissolution of lithium polysulfides in the liquid electrolyte, low electrical conductivity of sulfur and large volume change during electrochemical cycling are the main technical challenges for practical applications. In this study, a systematic first-principles density functional theory calculation is adopted to understand the interactions between graphene and graphene with oxygen containing functional groups (hydroxyl, epoxy and carboxyl groups) and sulphur (S8) and long chain lithium polysulfides (Li2S8 and Li2S4). We find the adsorption is dominated by different mechanisms in sulphur and lithium polysulfides, i.e. van der Waals attraction and formation of coordinate covalent Li–O bonds. The adsorption strength is dependent on the inter-layer distance and electron rich functional groups. Through these mechanisms, sulphur and lithium polysulfides can be successfully retained in porous graphene, leading to improved conductivity and charge transfer in the cathode of Li–S batteries.


Introduction

Recently, there has been an increasing demand for higher energy density in energy storage systems such as batteries and supercapacitors for portable devices and electric vehicles. With an energy capacity of 1673 mA h g−1 and a specific energy of 2600 W h kg−1, Li–S batteries are emerging as one of the promising candidates for next generation rechargeable batteries.1–3 However, dissolution of lithium polysulfides in the liquid electrolyte, low electrical conductivity of sulfur and its final discharge products, and large volume change during electrochemical cycling are the main technical challenges.4,5 To solve these problems, recent efforts have been put into the design of nanostructured electrodes to improve their capacity and cycling performance.6,7 It has been confirmed that composite cathodes consisting of sulfur and nanostructured carbon materials such as meso-micro porous carbon,8 carbon spheres9,10 carbon nanotubes (CNTs),11 graphene,12–15 and graphene oxide,16,17 could mitigate the polysulfide shuttle via physical confinement of soluble polysulfides within the conductive carbonaceous structures. Among the various types of carbon materials, graphene has attracted much attention due to its good electrical conductivity and high surface area. Another alternative is reduced graphene oxide (rGO), a layered material with graphene domains and residual functional groups such as hydroxyl, epoxy and carboxyl groups.18,19 Consequently, attempts were made to design cathodes with porous graphene and rGO.20–23 Chen's group24 employed a three-dimensional rGO sponge to produce a sulfur nanogranular film-coated composite cathode with a reversible capacity of 1080 mA h g−1 at 0.1C rate and 86.2% capacity retention after 500 cycles at 1.0C rate.

In General, microporous anchoring materials (pore size < 2 nm) could successfully confine polysulfides to achieve good electro activity and long-term cyclability.25–27 However it was realized that anchoring materials with non-polar surfaces, alone cannot mitigate the polysulfide shuttle as they fail to make sufficient interaction with polar lithium polysulfides, and further surface modification is generally needed to chemically bind polysulfides onto the carbonaceous matrix. Recently, density functional theory (DFT) studies were carried out to investigate the discharge mechanisms in Li–S batteries,28,29 and several investigators modelled the interactions between defective graphene, heteroatom doped graphene, and lithium polysulfides.30–32 However, the anchoring mechanisms of porous graphene functionalized with different oxygen groups during lithiation process have not been well understood. In addition, it is unclear if polar groups can facilitate the adsorption of non-polar S8 as they interact with lithium polysulfides and if the functional groups can build up a barrier for electron transfer at the interface between S8 and the substrate.

In this paper, we report a systematic first-principle investigation on the anchoring effects of microporous epoxy-, hydroxyl- and carboxyl-functionalized graphene and elucidate the mechanisms responsive for interfacial interaction and electron transfer.

Computational methods

As shown in Fig. 1a, a basic structure with two parallel graphene layers is constructed to simulate porous graphene. The interlayer distance is changed in the range of 7.5–20 Å to understand how the pore size affects the adsorption of S compounds. The smallest interlayer distance (pore size) is set to 7.5 Å to effectively accommodate the dimensions of cyclo-S8 (∼0.7 nm).33,34 To model oxygen functionalized graphene (OFG) which consists of hydroxyl, epoxy and carboxyl functional groups, three distinct microporous structures are built using the basic model shown in Fig. 1a with a interlayer distance of 12.5 Å. For simplicity, they are referred as hG, eG and cG, respectively (Fig. 1b). Based on literature,35,36 hydroxyl and epoxy groups are introduced on the basal plane and a carboxyl group is introduced at the edge of the plane corresponding to oxygen atomic concentration of 1.4–2.6%.
image file: c7ra11628d-f1.tif
Fig. 1 (a) Microporous graphene (b) microporous structures of hG, eG and cG respectively (c) geometrically optimized minimum energy structures of S8, Li2S8, Li2S6 and Li2S4 respectively (carbon, sulfur, oxygen, hydrogen and lithium atoms are represented by grey, yellow, red, white and purple, respectively).

Molecule configurations of S compounds, adsorption energies of S compounds to microporous structures (graphene and OFG), charge transfer from S compounds at different lithiation stages, and density of states near Fermi energy region are examined using DFT calculations with DMol3 package37,38 of Materials Studio 2016. Electron–electron exchange correlations are described by generalized gradient approximations (GGA) of the Perdew–Burke–Ernzerhof (PBE) functional.39 To consider the van der Waals interaction, a semi-empirical dispersion potential of DFT-D2 method of Grimme is used.40 Energy, maximum force and maximum displacement are set to 2.0 × 10−5 ha, 4 × 10−3 ha A−1 and 5 × 10−3 A, respectively. Self-consistent-field density is set to 1 × 10−5 eV. Double numerical plus polarization (DNP) is selected as the basis set and effective core potentials (ECP) are employed to describe the core electrons. The corresponding k point grid is generated by the Monkhrost–Pack technique41 for the Brillouin zone sampling and simulations are carried out using a 5 × 5 × 1 grid.

The binding energy (Ebind) of each Li2Sx molecules can be calculated as, Ebind = [ELi2Sx − (2ELi + ESx)]/2, where ELi2Sx, ELi and ESx are the energies of Li2Sx, Li atom and Sx molecule respectively. The adsorption energy (Eads) of each S species to various graphene structures are calculated according to, Eads = ETotal − (Egraphene + ES), where ETotal, Egraphene and ES are the total energies of the system, graphene structure and S containing cluster (S8, Li2S8 or Li2S4), respectively. According to the equation, more negative adsorption energy indicates stronger interaction between graphene structure and S cluster. Initially, sulfur species were placed with different orientations inside the graphene structure. Eight different potential configurations of each S compound were considered and after relaxation, corresponding adsorption values were compared to obtain the most stable configuration (Fig. S1). The energy difference between intact structure (Li2S8/Li2S4) and decomposed structure (Li + LiS8/Li + LiS4) is calculated according to, ΔE = E(Li2Sx + graphene) − E(Li + LiSx + graphene), where x = 4 and 8.

Results and discussion

Effect of microporous graphene on adsorption of S species

To find the structures with the minimum energy in S8, Li2S8, Li2S6 and Li2S4, geometry optimization was carried out to obtain possible linear or closed atomic arrangements.30,42,43Fig. 1c shows the stable ground state structures of S8 and lithium polysulfides. Their bond lengths are consistent with previous investigations.30 The calculated bond length, charge of each atom via the Mulliken population analysis (MPA) and the binding energies of lithium polysulfides are shown in Table 1.
Table 1 Average bond distances, average charges and binding energies
  Li2S8 Li2S6 Li2S4
Li–S (Å) 2.450 2.429 2.405
S–S (Å) 2.103 2.138 2.116
Charge on Li 0.44 0.44 0.46
Binding energy (eV) −3.39 −4.03 −3.50


The most stable configuration of sulfur is the cycloocta-S which consists of covalently bonded 8 sulfur atoms in a crown formation with a calculated average S–S bond length of 2.092 Å. During the initial phase of the discharge process, S8 is reduced to S82−, S62−, S42− forming high-order lithium polysulfides and further gets reduced to S22− and S2− forming low-order lithium polysulfides at the latter stages.2 According to our calculations, as x in Li2Sx decreases from 8 to 4, the average Li–S distance decreases from 2.450 Å to 2.405 Å, implying that the interaction between Li and S atom strengthens as the x decreases. The average S–S bond lengths of Li2S8, Li2S6 and Li2S4 are 2.103, 2.138 and 2.118 Å, respectively. Our results revealed that S–S bond lengths of polysulfides are larger than that of S8 molecule which is 2.092 Å due to the less covalent nature of the S cluster within the lithium polysulfide. Due to the weak interaction between Li and S, high order polysulfides tend to ionize easily into Li and polysulfide ions in the electrolyte than low-order Li2S2 and Li2S. MPA is a qualitative method to understand how the charge is distributed among each atom. According to MPA, the average net charge of Li is 0.44, 0.44 and 0.46 in Li2S8, Li2S6 and Li2S4 respectively. Therefore, the interaction between Li and S in Li2S8 is the weakest due to the largest Li–S distance and the lowest average net charge of Li. Moreover, Li2S8 exhibits the lowest binding energy of −3.39 eV, suggesting that it is the most soluble compound among these lithium polysulfides during a charge/discharge cycle, consistent with previous experiments.44 As an intermediate which presents in both high and low voltage plateau regimes of charge/discharge cycle, Li2S4 plays a vital role in the Li–S redox reaction. Elemental S is the starting point of the multi-electron-transfer cathode reaction, and hence the retention of S8 is important for the long lasting performance of a Li–S cell. For these reasons, the interaction between S8, Li2S8 and Li2S4 and the porous graphene deserves further investigation.

As shown in Fig. 2a (Table S1), the variation of energy profiles indicates that the adsorption energy (Eads) increases as the pore size decreases. The adsorption energy of S compounds to a typical anchoring material of the cathode depends on (1) the chemical interaction between the lithium polysulfide and the anchoring material in which a covalent bond can be formed between the Li atom in Li2Sx and the functional group of the anchoring material and/or between the S atom of Li2Sx and the functional group of the anchoring material and (2) the physical van der Waals attraction. Unlithiated S8 is a non-polar molecule, and does not form any chemical interaction with graphene domains. Consequently, the interfacial interaction is mainly governed by the physical interaction. However, as the lithiation begins, apart from the physical interaction, a chemical interaction between Li atoms and the anchoring material is formed. In graphene, regardless of different lithiation stages the physical interaction overpowers the chemical interaction and there is a significant increase of adsorption strength as the pore size decreases from 10 Å to 7.5 Å. When the pore size is 7.5 Å, the adsorption energies are −1.55 eV, −1.50 eV and −1.22 eV for Li2S8, Li2S4 and S8, respectively.


image file: c7ra11628d-f2.tif
Fig. 2 (a) Adsorption energies of S8, Li2S8 and Li2S4 to microporous graphene. Deformation charge density and adsorption site of (b) S8, (c) Li2S8, and (d) Li2S4.

The effect of van der Waals interaction on the confinement of S species inside microporous graphene can be visualized from the deformation charge density analysis (Fig. 2b–d), in which the increase and decrease of total electron density are denoted in blue and green respectively. The electron density difference is expressed as, Δρ = ρTotal − (ρgraphene + ρS), where ρTotal, ρgraphene and ρS are the electron densities of the system, graphene structure and the S-containing cluster respectively. It can be seen that charge is transferred inside the S species and inside the graphene surfaces, but no apparent charge transfer occurs between them, suggesting no strong chemical interaction. However, by thoroughly comparing the deformation charge density of Li2S8 and Li2S4 with their adsorption sites (Fig. 2c–d), it can be seen that Li atoms in S complexes have slightly moved away from S and towards the graphene surface, indicating a slight chemical attraction between Li and C atoms which explains why Li2S8 and Li2S4 exhibit higher adsorption energy values than S8. Due to the fact that unlithiated S8 is confined inside a narrow pore which induces a strong physical interaction from both sides of the graphene surfaces, the adsorption value of −1.22 eV exhibited by microporous graphene with a pore size of 7.5 Å is much higher than values recorded for anchoring materials like V2O5, MoS2 and phosphorene.45,46

Influence of different functional groups of OFG towards adsorption of S species

The interactions between S8, Li2S4 and Li2S8 molecules with functional groups hG, eG and cG, were simulated and the adsorption energies are summarized in Fig. 3 (Table S2). Fig. 4a shows the final optimized structures of S8 on pure graphene, hG, eG and cG, respectively. Note that in all the cases, the highest adsorption energy was given when S8 adsorbed from the opposite direction of the functional group and parallel to the graphene surface. According to our calculations, it was found that when compared with graphene, hG, eG and cG exhibit almost similar adsorption energy to S8, illustrating that oxygen functional groups do not have a major influence on the adsorption of S8, as the adsorption strength is dominated by van der Waals attraction.
image file: c7ra11628d-f3.tif
Fig. 3 Adsorption values of (a) S8, (b) Li2S8, and (c) Li2S4 interacting with graphene, hG, eG and cG.

image file: c7ra11628d-f4.tif
Fig. 4 (a) S8 adsorption on graphene, hG, eG and cG. Oxygen and hydrogen atoms are represented by red and white respectively. (b) Electron transfer from S8 to different substrates, and (c) density of states near Fermi energy region for S8 on graphene, hG, eG and cG.

To further investigate the interaction between S8 and porous graphene, MPA is applied to determine how much charge transfer has occurred between them. Interestingly, 0.09 electrons have been transferred from S8 to the graphene surface but there is only a minor contribution to charge transfer between S8 and OFG where only 0.04, 0.03 and 0.06 electrons have been withdrawn from hG, eG and cG respectively (Fig. 4b). Fig. 4c shows the density of states (DOS) near Fermi energy (Ef) for the adsorption system of S8 in graphene and OFG, in which the Ef is set to 0 eV and represented by a vertical dotted line. For S8 adsorbed graphene system, the DOS value at Ef is found to be 4.33 electrons per eV and it reduces to 4.00 and 3.92 electrons per eV when S8 gets adsorbed to hG and cG, respectively (Table S3). Furthermore DOS value at Ef drops to 3.15 electrons per eV when eG adsorbs S8, suggesting that electron transfer at the interface is slightly hindered due to the oxygen functional groups on the graphene surface which then leads to building up of ohmic resistance at the S8/OFG interface. This phenomenon is also reported by Shiqi et al. where they discovered that S8 is decoupled from the graphitic surface due to the presence of Triton X-100, a non-ionic surfactant with a polyethylene oxide chain.47

As for the adsorption of Li2S8 molecule, cG exhibits the highest adsorption energy of −1.29 eV and all the oxygen functional groups show significantly higher adsorption energies compared to pristine graphene (Fig. 3b). When interacting with Li2S4, hG exhibits the highest adsorption energy of −1.43 eV and it is almost over 2 times higher than pristine graphene (Fig. 3c). Fig. 5b–i shows the optimized structures of Li2S8 and Li2S4 on hG, eG and cG, respectively. To further get an understanding of the magnitude of the Li–O interaction, we optimized the geometry of a Li2O crystal (Fig. S2) so that it could be used as a benchmark to compare the bonding nature of Li and O in the optimized models.48 The calculated Li–O distances of adsorption geometries are almost similar with the strong ionic bond length of Li2O molecule implying the existence of a strong attraction between Li and O in all the cases. hG had the shortest Li–O distance of 1.87 Å when it interacted with Li2S8 while the longest Li–O distance of 2.01 Å was given by eG when it interacted with Li2S4 (Table S4).


image file: c7ra11628d-f5.tif
Fig. 5 Mulliken charge distribution before and after adsorption of Li2S8 and Li2S4 of (a–c) hG, (d–f) eG and (g–i) cG respectively. Note that charge of each atom is shown in parentheses. The shortest distance between Li atom and O atom is also shown in angstroms (Å).

As O atom has a high electronegativity of 3.44 and Li atom has a low electronegativity of 0.98, the strong attraction between Li and O can be explained by the Lewis-acid base theory. Since epoxy, hydroxyl and carboxyl groups all consist of O atoms with lone electron pairs in their outer p orbitals; they act as electron pair donors (Lewis bases). These Lewis-base sites are attracted by the terminal Li atoms of Li2S8 and Li2S4 which act as strong Lewis acids according to the Lewis-acid base theory. Fig. 6 shows the deformation charge density corresponding to Li2S8 and Li2S4 adsorption sites. A significantly high electron density is visible around the lone pairs of the O atom strengthening the fact that extra pairs of electrons act as electron rich donor to interact with strong Lewis acid of Li ion to form a coordinate covalent bond. The strong attraction between negatively charged oxygen atom and the positively charged Li atom can be further illustrated in Fig. 5 where the Mulliken charge distribution of hG, eG and cG before and after adsorption of Li2S8 and Li2S4 are presented. Because of the polarization of the O atom by the terminal Li atom of Li2S8, 0.12, 0.13 and 0.12 electrons are withdrawn by the closest O atom of each eG, hG and cG respectively. The same phenomenon occurs when those substrates interact with Li2S4 where 0.06, 0.12 and 0.13 electrons have been transferred to the O atom. Therefore functionalized carbon materials consisting of highly electronegative atoms with lone electron pairs, are good candidates for the immobilization of high order lithium polysulfides due to the moderate adsorption ability. Recently, a composite cathode made out of chlorine-reinforced carbon nanofibers reported to have enhanced cycling performance, validating our theoretical predictions.49


image file: c7ra11628d-f6.tif
Fig. 6 Deformation charge density at Li2S8 and Li2S4 adsorption sites of (a and b) hG, (c and d) eG and (e and f) cG. (The increase/decrease of electron density is denoted by blue/green, respectively.)

Though there is a relatively strong interaction between lithium polysulfides and the OFG, the polysulfide molecule itself remains intact without being dissociated into short chain lithium polysulfides upon adsorption. We observed that after adsorption, the internal Li–S bond length increased slightly from 2.45 Å to 2.48 Å in Li2S8 and 2.41 Å to 2.46 Å in Li2S4 (Table S4). Various layered metal oxides and sulfides were reported to induce a strong chemical bond with lithium polysulfides in terms of Li–O, Li–S or M–S bonds (M represents metal oxides/sulfides).50,51 However, recently it was discovered that strong interaction with polysulfides could interfere on S reduction reactions and sometimes weaken the Li–S bond causing the dissociation of the Li2Sx molecule. Such separation between Li and S atom eventually could lead up to the formation of Li+ and Sx2− ions and as a result sulfur could be dissolved in the electrolyte adversely affecting the performance of the Li–S cell.49,50,52 Furthermore, Yu et al. revealed the importance of employing a cathode material with a moderate binding capability which allows a small amount of polysulfide dissolution in order to improve the stability of the solid electrolyte interface (SEI) on the lithium anode.49

To further investigate the possibility of the decomposition of Li2Sx by the weakening of the Li–S bond due to the attraction of OFG, we considered the adsorption of the decomposed LiS8 and LiS4 structures along with an isolated Li atom and calculated the energy difference (ΔE) between the decomposed structure and the intact structure (Li2S8 and Li2S4). A negative value for ΔE indicates that the intact structure has lower energy than the decomposed structure and the intact structure is energetically stable. Fig. 7 shows the adsorbed Li + LiS8 and Li + LiS4 clusters on graphene, hG, eG and cG. According to ΔE values in all the cases, the intact structure proved to be energetically preferable over the decomposed structure. Therefore microporous graphene decorated with oxygen functional groups could be considered as an effective choice as a cathode material for Li–S batteries which strikes a balance between adsorption strength and intactness of high order lithium polysulfides.


image file: c7ra11628d-f7.tif
Fig. 7 Atomic configurations and ΔE values of Li + LiS8 and Li + LiS4 clusters adsorbed on graphene (a and e), hG (b and f), eG (c and g) and cG (d and h) respectively.

To gain further insights into the bond interaction between OFG and Li2S8/Li2S4, we analysed the atomic partial density of states (PDOS) near Ef for the adsorption systems (Fig. S3). The PDOS for Li-1 and O-1 atoms (indicated by black and blue colour) are seen to overlap at the upper part of the valence band (just below the Fermi level) suggesting that there is a hybridization between Li-2s and O-2p orbitals. This further ascertains the existence of a Li–O covalent bond which enhances the interaction between the substrate and the polysulfides and thereby mitigates the polysulfide shuttle effectively.

The DOS near Fermi energy for hG, eG and cG before and after the anchoring of S8, Li2S8 and Li2S4 are presented in Fig. S4, in which hG, eG and cG exhibit metallic nature with 3.98, 3.16 and 3.92 electrons per eV DOS values at Ef respectively at oxygen atomic concentrations of 1.4–2.6%. The good electrical conductivity of reduced graphene oxide has also been demonstrated by previous computational studies53,54 and, furthermore Stankovich et al.55 produced rGO nanosheets with a significantly high conductivity value of (∼2 × 102 S m−1) which closely approaches that of pristine graphene, even at an atomic C/O ratio of ∼10.

The strong affinity towards Li2S8/Li2S4 by the OFG is clearly observed from the DOS values of hG, eG and cG after anchoring of Li2S8 and Li2S4. When Li2S8 and Li2S4 are adsorbed, due to the newly formed Li–O covalent bond between Li atom of Li2S8/Li2S4 and O atom of the functional groups, more electrons are transferred to the adsorption sites and the DOS curve exhibits a higher value at Ef when compared with the S8 adsorbed system. Therefore polar groups maintain a strong interaction with the lithium polysulfides in the discharge/charge cycle of the Li–S battery and in the meantime improve the electrical conductivity of the cathode by facilitating the charge transfer at the interface.

In a typical Li–S battery, the organic electrolyte consists of bis(trifluoromethane)sulfonimide lithium salt (LiTFSI) dissolved in a mixture of 1,3-dioxolane (DOL) and 1,2-dimethoxyethane (DME). To investigate the influence of the electrolyte towards the discharge/charge cycle of the battery, we calculated the adsorption energy values of DOL and DME towards S8, Li2S8 and Li2S4 (Fig. 8a) and compared them with the values we gained for graphene and OFG. For ease of comparison, we used the average adsorption values of hG, eG, and cG as the adsorption value for OFG in each lithiation stage. According to Fig. 8b, adsorption energies of graphene and OFG towards S8 are almost similar and they are significantly higher than DOL and DME, suggesting that van der Waals interaction is more than sufficient to retain S8 in the cathode. However Li2S8 and Li2S4 are barely adsorbed by graphene and there's a high possibility of these long chain polysulfides being dissolved in the liquid electrolyte. When compared with DOL and DME, adsorption energy of OFG towards Li2S8/Li2S4 is significantly high proving the fact that functional groups could confine lithium polysulfides inside the cathode and mitigate the shuttle effect.


image file: c7ra11628d-f8.tif
Fig. 8 (a) Interaction of S8, Li2S8 and Li2S4 with DOL and DME. (b) Adsorption energies of S8, Li2S8 and Li2S4 interacting with DOL, DME, graphene and OFG.

Based on the insights gained from our theoretical study, we can conclude that adsorption strength and electron transfer at the interface are two critical aspects which have to be considered when selecting an anchoring material for the cathode in Li–S batteries. At the initial stage of the discharge cycle, confinement of S8 is controlled by the van der Waals interaction, and it is highly unlikely for S8 being dissolved in the electrolyte due to the poor attraction of the solvents towards S8 molecules. A microporous substrate is more suitable to confine S8 as it can be confined inside a narrow pore which induces a strong physical interaction. Surface functionalization of graphene has no additional influence on attracting S8 to its surface since almost no chemical interaction can be formed between functional groups and non-polar S8. Furthermore functional groups adversely affect the electron transfer and increase the ohmic resistance at the substrate/S8 interface. However functional groups have a major impact on mitigating the shuttle effect and improving the capacity retention of Li–S cell by ensuring a strong interaction with Li2S8/Li2S4. On the other hand, the strong interaction with lithium polysulfides facilitates the interfacial electron transfer through the Li–O bond which can be explained by the Lewis-acid base theory.

Conclusions

In this work, we have investigated the interactions between graphene with oxygen containing functional groups (hydroxyl, epoxy and carboxyl groups) and sulphur (S8) and long chain lithium polysulfides (Li2S8 and Li2S4), via density functional theory computations. During the initial unlithiated stage, we find the interaction between sulphur and pristine graphene and OFG are dominated by van der Waals attraction. Although graphene and OFG exhibits almost the same adsorption capability towards non-polar S8, the functional groups develop a slight barrier for electron transfer at the interface and increase the ohmic resistance. The highest adsorption energy is observed when the distance between two graphene layers approaches to about 7.5 Å. During the lithiation stage, surface functionalization of graphene significantly enhances the interaction with Li2S8/Li2S4 by forming a coordinate covalent Li–O bond. Due to the covalent nature of the Li–O bond, polysulfides are well retained inside the cathode and it also improves the conductivity of the electrode upon the deposition of Li2S8 and Li2S4 by facilitating the interfacial charge transfer. Furthermore, our work explains the reason why porous graphene with oxygen functional groups is more effective as a cathode material compared to materials with too strong interactions which could cause destruction effects on the adsorbed lithium polysulfides. Due to the moderate binding affinity, OFG strikes a balance between adsorption strength and intactness of high order lithium polysulfides. Therefore based on our simulations, we suggest that microporous graphene decorated with hydroxyl, epoxy and carboxyl functional groups can successfully anchor sulfur and lithium polysulfides via moderate interactions, leading to improved conductivity and charge transfer in the cathode of Li–S batteries.

Conflicts of interest

The authors declare no conflict of interest.

Acknowledgements

The authors wish to acknowledge the High Performance Computing and Research Support (HPC) group in Queensland University of Technology for access to its computation resources. KW would like to acknowledge Queensland University of Technology for financial support via APA and IPRS scholarships. The financial support from an ARC Discovery Project (DP150101717) is appreciated.

References

  1. A. Manthiram, S.-H. Chung and C. Zu, Adv. Mater., 2015, 27, 1980–2006 CrossRef CAS .
  2. Z. Lin and C. Liang, J. Mater. Chem. A, 2015, 3, 936–958 RSC .
  3. A. Manthiram, Y. Fu, S.-H. Chung, C. Zu and Y.-S. Su, Chem. Rev., 2014, 114, 11751–11787 CrossRef CAS .
  4. D.-W. Wang, Q. Zeng, G. Zhou, L. Yin, F. Li, H.-M. Cheng, I. R. Gentle and G. Q. M. Lu, J. Mater. Chem. A, 2013, 1, 9382–9394 RSC .
  5. H. Kim, H.-D. Lim, J. Kim and K. Kang, J. Mater. Chem. A, 2014, 2, 33–47 RSC .
  6. M. A. Pope and I. A. Aksay, Adv. Energy Mater., 2015, 5, 1500124 CrossRef .
  7. L. Ma, K. E. Hendrickson, S. Wei and L. A. Archer, Nano Today, 2015, 10, 315–338 CrossRef CAS .
  8. X. Ji, K. T. Lee and L. F. Nazar, Nat. Mater., 2009, 8, 500–506 CrossRef CAS .
  9. J. Liu, T. Yang, D.-W. Wang, G. Q. Lu, D. Zhao and S. Z. Qiao, Nat. Commun., 2013, 4, 2798 CrossRef .
  10. R. Fang, S. Zhao, P. Hou, M. Cheng, S. Wang, H.-M. Cheng, C. Liu and F. Li, Adv. Mater., 2016, 28, 3374–3382 CrossRef CAS .
  11. X.-B. Cheng, J.-Q. Huang, Q. Zhang, H.-J. Peng, M.-Q. Zhao and F. Wei, Nano Energy, 2014, 4, 65–72 CrossRef CAS .
  12. L. Fei, X. Li, W. Bi, Z. Zhuo, W. Wei, L. Sun, W. Lu, X. Wu, K. Xie, C. Wu, H. L. W. Chan and Y. Wang, Adv. Mater., 2015, 27, 5936–5942 CrossRef CAS .
  13. J. Zhu, D. Yang, Z. Yin, Q. Yan and H. Zhang, Small, 2014, 10, 3480–3498 CrossRef CAS .
  14. M. Yu, R. Li, M. Wu and G. Shi, Energy Storage Materials, 2015, 1, 51–73 CrossRef .
  15. R. Chen, T. Zhao, J. Lu, F. Wu, L. Li, J. Chen, G. Tan, Y. Ye and K. Amine, Nano Lett., 2013, 13, 4642–4649 CrossRef CAS .
  16. L. Ji, M. Rao, H. Zheng, L. Zhang, Y. Li, W. Duan, J. Guo, E. J. Cairns and Y. Zhang, J. Am. Chem. Soc., 2011, 133, 18522–18525 CrossRef CAS .
  17. C. Zu and A. Manthiram, Adv. Energy Mater., 2013, 3, 1008–1012 CrossRef CAS .
  18. S. Pei and H.-M. Cheng, Carbon, 2012, 50, 3210–3228 CrossRef CAS .
  19. R. K. Singh, R. Kumar and D. P. Singh, RSC Adv., 2016, 6, 64993–65011 RSC .
  20. K. C. Wasalathilake, G. A. Ayoko and C. Yan, in Recent Advances in Graphene Research, InTech, 2016,  DOI:10.5772/63554 .
  21. X. Huang, B. Sun, K. Li, S. Chen and G. Wang, J. Mater. Chem. A, 2013, 1, 13484–13489 RSC .
  22. K. H. Kim, Y.-S. Jun, J. A. Gerbec, K. A. See, G. D. Stucky and H.-T. Jung, Carbon, 2014, 69, 543–551 CrossRef CAS .
  23. G. Zhou, L.-C. Yin, D.-W. Wang, L. Li, S. Pei, I. R. Gentle, F. Li and H.-M. Cheng, ACS Nano, 2013, 7, 5367–5375 CrossRef CAS .
  24. W. Ahn, M. H. Seo, Y.-S. Jun, D. U. Lee, F. M. Hassan, X. Wang, A. Yu and Z. Chen, ACS Appl. Mater. Interfaces, 2016, 8, 1984–1991 CrossRef CAS .
  25. Z. Li, X. Li, Y. Liao, X. Li and W. Li, J. Power Sources, 2016, 334, 23–30 CrossRef CAS .
  26. K. Yang, Q. Gao, Y. Tan, W. Tian, L. Zhu and C. Yang, Microporous Mesoporous Mater., 2015, 204, 235–241 CrossRef CAS .
  27. X. Gu, Y. Wang, C. Lai, J. Qiu, S. Li, Y. Hou, W. Martens, N. Mahmood and S. Zhang, Nano Res., 2015, 8, 129–139 CrossRef CAS .
  28. M. Vijayakumar, N. Govind, E. Walter, S. D. Burton, A. Shukla, A. Devaraj, J. Xiao, J. Liu, C. Wang, A. Karim and S. Thevuthasan, Phys. Chem. Chem. Phys., 2014, 16, 10923–10932 RSC .
  29. L. Wang, T. Zhang, S. Yang, F. Cheng, J. Liang and J. Chen, J. Energy Chem., 2013, 22, 72–77 CrossRef CAS .
  30. S. P. Jand, Y. Chen and P. Kaghazchi, J. Power Sources, 2016, 308, 166–171 CrossRef CAS .
  31. T.-Z. Hou, X. Chen, H.-J. Peng, J.-Q. Huang, B.-Q. Li, Q. Zhang and B. Li, Small, 2016, 12, 3283–3291 CrossRef CAS .
  32. H. Ting-Zheng, P. Hong-Jie, H. Jia-Qi, Z. Qiang and L. Bo, 2D Mater., 2015, 2, 014011 CrossRef .
  33. Y.-X. Yin, S. Xin, Y.-G. Guo and L.-J. Wan, Angew. Chem., Int. Ed., 2013, 52, 13186–13200 CrossRef CAS .
  34. H. Liao, H. Ding, B. Li, X. Ai and C. Wang, J. Mater. Chem. A, 2014, 2, 8854–8858 RSC .
  35. D. R. Dreyer, S. Park, C. W. Bielawski and R. S. Ruoff, Chem. Soc. Rev., 2010, 39, 228–240 RSC .
  36. A. S. Dobrota, I. A. Pasti, S. V. Mentus and N. V. Skorodumova, Phys. Chem. Chem. Phys., 2016, 18, 6580–6586 RSC .
  37. B. Delley, J. Chem. Phys., 1990, 92, 508–517 CrossRef CAS .
  38. B. Delley, J. Chem. Phys., 2000, 113, 7756–7764 CrossRef CAS .
  39. J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett., 1996, 77, 3865–3868 CrossRef CAS .
  40. E. R. McNellis, J. Meyer and K. Reuter, Phys. Rev. B: Condens. Matter Mater. Phys., 2009, 80, 205414 CrossRef .
  41. H. J. Monkhorst and J. D. Pack, Phys. Rev. B: Solid State, 1976, 13, 5188–5192 CrossRef .
  42. B. Wang, S. M. Alhassan and S. T. Pantelides, Phys. Rev. Appl., 2014, 2, 034004 CrossRef .
  43. Z. Ji, B. Han, Q. Li, C. Zhou, Q. Gao, K. Xia and J. Wu, J. Phys. Chem. C, 2015, 119, 20495–20502 CrossRef CAS .
  44. M. Hagen, P. Schiffels, M. Hammer, S. Dörfler, J. Tübke, M. J. Hoffmann, H. Althues and S. Kaskel, J. Electrochem. Soc., 2013, 160, A1205–A1214 CrossRef CAS .
  45. Q. Zhang, Y. Wang, Z. W. Seh, Z. Fu, R. Zhang and Y. Cui, Nano Lett., 2015, 15, 3780–3786 CrossRef CAS .
  46. J. Zhao, Y. Yang, R. S. Katiyar and Z. Chen, J. Mater. Chem. A, 2016, 4, 6124–6130 RSC .
  47. S. Li, T. Mou, G. Ren, J. Warzywoda, B. Wang and Z. Fan, ACS Energy Lett., 2016, 1, 481–489 CrossRef CAS .
  48. N. G. Szwacki and T. Szwacka, Basic Elements of Crystallography, Pan Stanford, 2nd edn, 2016 Search PubMed .
  49. L. Fan, H. L. Zhuang, K. Zhang, V. R. Cooper, Q. Li and Y. Lu, Adv. Sci., 2016, 3, 1600175 CrossRef .
  50. Y. Zhong, K. R. Yang, W. Liu, P. He, V. Batista and H. Wang, J. Phys. Chem. C, 2017, 121, 14222–14227 CrossRef CAS .
  51. X. Chen, H.-J. Peng, R. Zhang, T.-Z. Hou, J.-Q. Huang, B. Li and Q. Zhang, ACS Energy Lett., 2017, 2, 795–801 CrossRef CAS .
  52. Y. Zhao and J. Zhao, Appl. Surf. Sci., 2017, 412, 591–598 CrossRef CAS .
  53. D. W. Boukhvalov and M. I. Katsnelson, J. Am. Chem. Soc., 2008, 130, 10697–10701 CrossRef CAS .
  54. M. Barhoumi, D. Rocca, M. Said and S. Lebègue, J. Chem. Phys., 2017, 146, 044705 CrossRef CAS .
  55. S. Stankovich, D. A. Dikin, R. D. Piner, K. A. Kohlhaas, A. Kleinhammes, Y. Jia, Y. Wu, S. T. Nguyen and R. S. Ruoff, Carbon, 2007, 45, 1558–1565 CrossRef CAS .

Footnote

Electronic supplementary information (ESI) available. See DOI: 10.1039/c7ra11628d

This journal is © The Royal Society of Chemistry 2018