Advanced
Mechanisms of Na adsorption on graphene and graphene oxide: density functional theory approach
Mechanisms of Na adsorption on graphene and graphene oxide: density functional theory approach
Carbon letters. 2015. Apr, 16(2): 116-120
Copyright © 2015, Korean Carbon Society
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : February 02, 2015
  • Accepted : March 03, 2015
  • Published : April 30, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Hye Sook Moon
Department of Organic Material Science and Engineering, Pusan National University, Busan 609-735, Korea
Ji Hye Lee
Department of Organic Material Science and Engineering, Pusan National University, Busan 609-735, Korea
Soonchul Kwon
School of Civil and Environmental Engineering, Georgia Institute of Technology, Atlanta, GA 30332, USA
Il Tae Kim
Department of Chemical and Biological Engineering, Gachon University, Seongnam 461-701, Korea
itkim@gachon.ac.kr
Seung Geol Lee
Department of Organic Material Science and Engineering, Pusan National University, Busan 609-735, Korea
itkim@gachon.ac.kr
Abstract
We investigated the adsorption of Na on graphene and graphene oxide, which are used as anode materials in sodium ion batteries, using density functional theory. The adsorption energy for Na on graphene was -0.507 eV at the hollow sites, implying that adsorption was favorable. In the case of graphene oxide, Na atoms were separately adsorbed on the epoxide and hydroxyl functional groups. The adsorption of Na on graphene oxide-epoxide (adsorption energy of -1.024 eV) was found to be stronger than the adsorption of Na on pristine graphene. However, the adsorption of Na on graphene oxide-hydroxyl resulted in the generation of NaOH as a by-product. Using density of states (DOS) calculations, we found that the DOS of the Na-adsorbed graphene was shifted down more than that of the Na-adsorbed graphene oxide-epoxide. In addition, the intensity of the DOS around the Fermi level for the Na-adsorbed graphene was higher than that for the Na-adsorbed graphene oxide-epoxide.
Keywords
1. Introduction
Lithium ion batteries (LIBs) are essential energy storage systems present in a number of devices used in our daily lives such as laptops and mobile phones. However, the developments made in energy storage technology have outpaced the increase in energy storage needs [1] . Energy storage technologies such as Li batteries are considered as alternatives to fossil fuels, which are rapidly depleting all over the world. However, because Li is a scarce resource on Earth, producing large-sized LIBs is limited by the exponential increase in its cost [2] . In this context, sodium ion batteries (SIBs) are excellent alternatives to LIBs because of their lower cost and the greater availability of Na compared to Li. In fact, Na is abundantly available and the de-insertion behavior of Na in SIBs is very similar to that of Li in LIBs [3] . Hard carbon materials that have graphite layers and amorphous regions are widely used as anodes in SIBs [4] . In addition, carbon black, amorphous carbon, mesoporous carbon, and carbon fiber are being experimentally tested as anode materials for SIBs [5] . Among the various carbon materials, graphene, which has excellent electrical characteristics and physical properties, is gaining attention as an anode material candidate for SIBs. In a recent study, David et al. [6] synthesized layered freestanding papers comprised of acid-exfoliated few-layer molybdenum disulfide (MoS 2 ) and reduced graphene oxide (rGO) flakes for use as self-standing flexible electrodes in SIBs. The electrode showed a good Na cycling ability with a stable charge capacity of approximately 230 mAh g -1 taking into account the total weight of the electrode and exhibited coulombic efficiencies up to approximately 99%. In addition, static uniaxial tensile tests performed on the crumpled composite papers showed high average strains to failure up to approximately 2%. Xie et al. [7] synthesized SnS 2 nanoplatelets at graphene nanocomposites using a morphology controlled hydrothermal method. When SnS 2 /graphene nanosheets were used as anode materials in SIBs, a high reversible specific Na-ion storage capacity of 725 mAh g -1 , stable cyclability, and an enhanced high-rate capability were observed.
Although such studies exist, the detailed mechanisms on Na adsorption on graphene and graphene oxide have not been fully elucidated yet. Therefore, in this work, we investigated the Na adsorption mechanisms on graphene-based anodes used in SIBs at the atomistic level by analyzing the adsorption energies and electronic properties using density functional theory (DFT).
2. Experimental
The DFT calculations were performed using the Vienna Ab-initio Simulation Package (VASP) [8 , 9] with the projector augmented wave (PAW) method [10] , and generalized gradient approximation (GGA) and Perdew-Burke-Ernzerhof (PBE) functionals [11 , 12] . The PBE functional has been used to successfully describe the interactions between an adsorbate and various surfaces [13 - 15] including sp 2 carbon-based materials [16 - 19] .
The energy cut-off was set at 500 eV, and the cell size was 9.840 × 9.840 × 15 Å. The calculations were performed under the periodic boundary conditions, with the vacuum separation distance set to 15 Å along the z-direction, to avoid interactions between the images. In the energy and electronic property calculations, k-point samplings for the Brillouin zone were performed using the 5 × 5 × 1 Monkhorst-Pack k-point mesh [20] . The atomic positions were optimized until the change in energy was less than 1 × 10 -6 eV per cell, and the force on each atom was less than 0.01 eV Å -1 . Spin-polarization and dipole corrections were also incorporated into the calculations. The adsorption energies of Na-adsorbed graphene and graphene oxide were defined by Eq. (1):
PPT Slide
Lager Image
, where E Graphene(oxide)+Na is the total energy of the Na-adsorbed graphene or graphene oxide; E Graphene(oxide) is the total energy of the graphene or graphene oxide, and E Na is the total energy of the isolated Na. Dispersion corrections [21] were also considered in all the calculations. Charge transfer was calculated using the density derived electrostatic and chemical (DDEC) method [22 , 23] . Graphene oxide was assumed to consist of epoxide and hydroxyl oxidation functional groups, and each group was introduced to the pristine graphene shown in Fig. 1 .
PPT Slide
Lager Image
Optimized structure of graphene oxide: (a) epoxide group and (b) hydroxyl group.
3. Results and Discussion
- 3.1. Adsorption energies and charge analysis
Na atoms were adsorbed onto hollow, bridge, and top sites of pristine graphene to determine the most stable position ( Fig. 2 ). The hollow site was the most stable position for Na atoms, and the adsorption energy was -0.507 eV. After Na adsorption, the optimized distance between the graphene and Na was 2.29 Å. The bridge sites were the second most stable adsorption sites for Na with an adsorption energy of -0.504 eV, whereas adsorption onto the top sites was the least stable with an adsorption energy of -0.493 eV. These results were confirmed by charge analysis showing that at all three sites, the charge was transferred from the Na to the graphene after adsorption, and the number of DDEC charges of the graphene at the hollow, bridge, and top sites was -0.657 e , -0.515 e , and -0.508 e , respectively. These results agreed well with the adsorption energy values. The adsorption energies, optimized distance between the Na atoms and the graphene, and charge distribution of the hollow, bridge, and top sites of the graphene are summarized in Fig. 3a and Table 1 .
PPT Slide
Lager Image
Na adsorption sites on the graphene: (a) hollow, (b) bridge, and (c) top.
PPT Slide
Lager Image
The optimized structures of the Na-adsorbed graphene and graphene oxide: (a) Na adsorption on graphene, (b) Na adsorption on graphene oxide-epoxide, and (c) initial and final structures of the Na-adsorbed graphene oxide-hydroxyl before and after Na adsorption.
Adsorption energy, distance between Na atoms and graphene, and charge distribution of the bridge, hollow, and top sites on the graphene
PPT Slide
Lager Image
Adsorption energy, distance between Na atoms and graphene, and charge distribution of the bridge, hollow, and top sites on the graphene
Graphene oxide was assumed to contain epoxide and hydroxyl functional groups, and Na was adsorbed onto each of these types of functional groups on the graphene oxide. First, Na was adsorbed vertically onto the oxygen atoms of the graphene oxide-epoxide (GO-epoxide). The adsorption energy for this case was -1.024 eV, which was approximately twice the value for the adsorption of Na onto the pristine graphene. We believe that the oxygen in the epoxide caused a doping effect, and subsequently, the adsorption energy was enhanced. In addition, the optimized distance between the graphene and Na was 2.068 Å. The optimized structure of the Na-adsorbed GO-epoxide and the charge transfer values for the graphene oxide and Na-adsorbed graphene oxide are summarized in Fig. 3b and Table 2 . The sum of the charges of the oxygen and carbon in the GO-epoxide was -0.287 e and +0.287 e before Na adsorption, respectively. However, after Na adsorption, the sum of the charges of the oxygen and carbon in the GO-epoxide was -0.672 e and -0.061 e , respectively. In other words, the oxygen and carbon in the GO-epoxide gained charges of -0.385 e and -0.348 e , respectively, as a result of the Na adsorption. Thus, we found that a total charge of +0.733 e was transferred from the Na to the oxygen and carbon in the graphene oxide during adsorption.
Charge distribution on the Na-adsorbed graphene oxide (epoxide and hydroxyl) compared with those on the pristine graphene oxide (epoxide and hydroxyl) and NaOH
PPT Slide
Lager Image
GO: graphene oxide.
Next, Na was adsorbed vertically onto the oxygen atom of the hydroxyl group in the graphene oxide-hydroxyl (GO-hydroxyl). Unlike graphene and GO-epoxide, the hydroxyl group dissociated as a result of Na adsorption and produced NaOH as a by-product ( Fig. 3c ). The optimized distance between the oxygen of the NaOH and graphene was 2.831 Å, and the adsorption energy between the NaOH and graphene was -0.402 eV. Kim et al. [24] showed that the epoxide groups are much more stable than the hydroxyl groups in graphene oxide by DFT with the GGA-PW91 functional. Because of the low binding energy, a hydroxyl group in the interior aromatic domain of the graphene oxide is unstable, and as a consequence, the hydroxyl group can be dissociated [25 , 26] . However, it should be noted that because the present investigation is based on DFT calculations that focus on the ground state of Na-graphene systems, the results obtained during the operation of SIBs may be different. Dissociated NaOH has a similar charge distribution as an isolated NaOH molecule. As a result, charges of -0.042 e were transferred from the NaOH to the graphene. The small amount of charge transferred implies that the NaOH and graphene molecules interacted via weak physical adsorption.
- 3.2. Electronic properties
We investigated the electronic properties of the graphene and graphene oxide by analyzing the density of states (DOS). The Fermi levels for all the DOS shifted to zero. The DOS of pristine and Na-adsorbed graphene are shown in Fig. 4a . After Na adsorption, however, the DOS of the Na-adsorbed graphene shifted down, and their intensity near the Fermi level increased, indicating that more electronic states can exist at the Fermi level. This could be attributed to the fact that electrons from Na are injected towards the graphene.
PPT Slide
Lager Image
The density of states (DOS) of the graphene and graphene oxide (GO) systems: (a) graphene and Na-adsorbed graphene, (b) GO-epoxide and Na-adsorbed GO-epoxide, (c) GO-hydroxyl and Na-adsorbed GO-hydroxyl, and (d) local DOS of the NaOH molecules and the carbon of the Na-adsorbed GO-hydroxyl.
Fig. 4b shows the DOS of the GO-epoxide and Na-adsorbed GO-epoxide. The Dirac point was not observed in the DOS of the GO-epoxide because the 2 p orbital of the oxygen hybridized with the π-orbital of the graphene, leading to a change in the Dirac point [27] . After Na adsorption, all the DOS were shifted down, and the intensity of DOS near the Fermi level increased, similar to the case of Na-adsorbed graphene. These results were also attributed to the injection of electrons from the Na to the GO-epoxide.
Fig. 4c shows the DOS of the pristine GO-hydroxyl and Na-adsorbed GO-hydroxyl. While the DOS of the GO-hydroxyl exhibited several peaks between -2 and 2 eV, the DOS of the Na-adsorbed GO-hydroxyl exhibited a broad line between -2 and 2 eV, similar to the case of the graphene. In addition, the intensity of the DOS near the Fermi level was significantly increased after Na adsorption. To investigate this in more detail, we calculated the local density of states (LDOS) of the carbon (graphene) after Na adsorption and the NaOH molecules ( Fig. 4d ). The LDOS of the carbon, shown in Fig. 4d , was similar to that of the graphene. Therefore, we concluded that the GO-hydroxyl reverted to the graphene, as a result of the dissociation of OH. In addition, a broad line between -2 and 2 eV in the DOS of the Na-adsorbed GO-hydroxyl was a result of the LDOS of the carbon. In the case of the LDOS of the NaOH molecules, the DOS of the NaOH exhibited a huge peak at the Fermi level, providing the energy levels of the molecular orbitals for the NaOH molecules rather than the band structure. Therefore, in Fig. 4c , the intensity of the DOS peak around the Fermi level was significantly increased after Na adsorption, which appears to be affected by the DOS of the dissociated NaOH.
As a result, the DOS of the Na-adsorbed graphene shifted down more than that of the GO-epoxide, and the intensity of the DOS around the Fermi level was higher than that of the Na-adsorbed GO-epoxide. Although the adsorption energy for Na onto GO-epoxide was higher than the corresponding value for graphene, after Na adsorption, the electronic properties of the graphene were better than those of the GO-epoxide because oxygen doping weakens the sp 2 - hybridization of the carbon atoms in the GO. In the case of the GO-hydroxyl, after Na adsorption, the OH group dissociated and produced NaOH molecules. In addition, the dissociated NaOH molecules were adsorbed onto the graphene by weak physical adsorption.
4. Conclusions
In this paper, we investigated the adsorption energies and electronic properties for the adsorption of Na onto graphene and graphene oxide, using the first-principles DFT method. Na adsorption was the most stable at the hollow sites of pristine graphene with an adsorption energy of -0.507 eV. In the case of the GO-epoxide, the Na adsorption energy was -1.024 eV, which was more negative than the adsorption energy for Na onto pristine graphene. In the case of the GO-hydroxyl, the hydroxyl group dissociated from the GO-hydroxyl, and NaOH was formed as a by-product. After Na adsorption, the DOS of both graphene and GO-epoxide shifted down, and the intensity of the DOS around the Fermi level increased due to electron injection in both cases. In addition, it was found that the DOS of the Na-adsorbed graphene shifted down more than that of the GO-epoxide, and the intensity of the DOS around the Fermi level for the Na-adsorbed graphene was higher than that for the Na-adsorbed GO-epoxide. Therefore, after Na adsorption, the electronic properties of the graphene were better than those of the GO-epoxide because of oxygen doping which weakened the sp 2 -hybridization of the carbon atoms in the GO.
Acknowledgements
This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (2014R1A1A1004096). This research was supported by the Gachon University research fund of 2014 (GCU-2014-0122).
References
Slater MD , Kim D , Lee E , Johnson CS 2013 Sodium-ion batteries Adv Funct Mater http://dx.doi.org/10.1002/adfm.201200691 23 947 -
Ellis BL , Nazar LF 2012 Sodium and sodium-ion energy storage batteries Curr Opin Solid State Mater Sci http://dx.doi.org/10.1016/j.cossms.2012.04.002 16 168 -
Palomares V , Casas-Cabanas M , Castillo-Martinez E , Han MH , Rojo T 2013 Update on Na-based battery materials. A growing research path Energy Environ Sci http://dx.doi.org/10.1039/C3ee41031e 6 2312 -
Kim SW , Seo DH , Ma XH , Ceder G , Kang K 2012 Electrode materials for rechargeable sodium-ion batteries: potential alternatives to current lithium-ion batteries Adv Energy Mater http://dx.doi.org/10.1002/aenm.201200026 2 710 -
Palomares V , Serras P , Villaluenga I , Hueso KB , Carretero-Gonzalez J , Rojo T 2012 Na-ion batteries, recent advances and present challenges to become low cost energy storage systems Energy Environ Sci http://dx.doi.org/10.1039/C2ee02781j 5 5884 -
David L , Bhandavat R , Singh G 2014 MoS2/graphene composite paper for sodium-ion battery electrodes Acs Nano http://dx.doi.org/10.1021/Nn406156b 8 1759 -
Xie XQ , Su DW , Chen SQ , Zhang JQ , Dou SX , Wang GX 2014 SnS2 nanoplatelet@graphene nanocomposites as high-capacity anode materials for sodium-ion batteries Chem Asian J http://dx.doi.org/10.1002/asia.201400018 9 1611 -
Kresse G , Furthmuller J 1996 Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set Phys Rev B http://dx.doi.org/10.1103/PhysRevB.54.11169 54 11169 -
Kresse G , Hafner J 1993 Abinitio molecular-dynamics for liquid-metals Phys Rev B http://dx.doi.org/10.1103/Phys-RevB.47.558 47 558 -
Kresse G , Joubert D 1999 From ultrasoft pseudopotentials to the projector augmented-wave method Phys Rev B http://dx.doi.org/10.1103/PhysRevB.59.1758 59 1758 -
Perdew JP , Burke K , Ernzerhof M 1996 Generalized gradient approximation made simple Phys Rev Lett http://dx.doi.org/10.1103/PhysRevLett.77.3865 77 3865 -
Perdew JP , Burke K , Wang Y 1996 Generalized gradient approximation for the exchange-correlation hole of a many-electron system Phys Rev B http://dx.doi.org/10.1103/Phys-RevB.54.16533 54 16533 -
Kwon S , Choi JI , Lee SG , Jang SS 2014 A density functional theory (DFT) study of CO2 adsorption on Mg-rich minerals by enhanced charge distribution Comput Mater Sci http://dx.doi.org/10.1016/j.commatsci.2014.07.042 95 181 -
Lee SG , Choi JI , Koh W , Jang SS 2013 Adsorption of beta-D-glucose and cellobiose on kaolinite surfaces: density functional theory (DFT) approach Appl Clay Sci http://dx.doi.org/10.1016/j.clay.2012.11.002 71 73 -
Kwon S , Lee SG , Chung E , Lee WR 2015 CO2 adsorption on H2O-saturated BaO(1 0 0) and induced barium surface dissociation Bull Korean Chem Soc 36 11 -
Koh W , Choi JI , Donaher K , Lee SG , Jang SS 2011 Mechanism of Li adsorption on carbon nanotube-fullerene hybrid system: a first-principles study ACS Appl Mater Interfaces http://dx.doi.org/10.1021/Am200018w 3 1186 -
Koh W , Choi JI , Jeong E , Lee SG , Jang SS 2014 Li adsorption on a Fullerene-Single wall carbon nanotube hybrid system: density functional theory approach Curr Appl Phys http://dx.doi.org/10.1016/j.cap.2014.09.031 14 1748 -
Koh W , Moon HS , Lee SG , Choi JI , Jang SS 2015 A first-principles study of lithium adsorption on a graphene-fullerene nanohybrid system ChemPhysChem http://dx.doi.org/10.1002/cphc.201402675 16 789 -
Koh W , Choi JI , Lee SG , Lee WR , Jang SS 2011 First-principles study of Li adsorption in a carbon nanotube-fullerene hybrid system Carbon http://dx.doi.org/10.1016/j.carbon.2010.09.022 49 286 -
Monkhorst HJ , Pack JD 1976 Special points for brillouin-zone integrations Phys Rev B http://dx.doi.org/10.1103/PhysRevB.13.5188 13 5188 -
Grimme S , Antony J , Ehrlich S , Krieg H 2010 A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu J Chem Phys http://dx.doi.org/10.1063/1.3382344 132 154104 -
Manz TA , Sholl DS 2012 Improved atoms-in-molecule charge partitioning functional for simultaneously reproducing the electrostatic potential and chemical states in periodic and nonperiodic materials J Chem Theory Comput http://dx.doi.org/10.1021/Ct3002199 8 2844 -
Manz TA , Sholl DS 2010 Chemically meaningful atomic charges that reproduce the electrostatic potential in periodic and nonperiodic materials J Chem Theory Comput http://dx.doi.org/10.1021/Ct100125x 6 2455 -
Kim MC , Hwang GS , Ruoff RS 2009 Epoxide reduction with hydrazine on graphene: a first principles study J Chem Phys http://dx.doi.org/10.1063/1.3197007 131 064704 -
Pei SF , Cheng HM 2012 The reduction of graphene oxide Carbon http://dx.doi.org/10.1016/j.carbon.2011.11.010 50 3210 -
Gao XF , Jang J , Nagase S 2010 Hydrazine and thermal reduction of graphene oxide: reaction mechanisms, product structures, and reaction design J Phys Chem C http://dx.doi.org/10.1021/Jp909284g 114 832 -
Yan JA , Chou MY 2010 Oxidation functional groups on graphene: structural and electronic properties Phys Rev B http://dx.doi.org/10.1103/Physrevb.82.125403 82 125403 -