Restoration of thermally reduced graphene oxide by atomic-level selenium doping

The use of reduced graphene oxide (rGO) suffers from irreparable damage because of topological defects and residual heteroatoms, which degrade the inherent properties of graphene. To restore its electrical transport properties, charge-transfer chemical doping with d-electron-rich heteroatoms has been proposed. Herein, we report the effects of atomic-level selenium doping in rGO. Using first-principles calculations, we found that selenium atoms could be selectively bonded in particular locations, such as the pseudo-edge sites of hole-cluster defects in the basal plane and edge defect sites of graphene; however, we found that the intrinsic topological defects of the basal plane were unfavorable for bonding. Numerous selenium atoms were introduced on the fully amorphorized rGO surface, inducing a dramatic change of its electrical transport properties by electron doping. The large metallic regions formed by the selenium atoms on rGOs led to the enhancement of electrical conductivity by 210 S cm–1 at 300 K. Moreover, the temperature-dependent conductivities (σ)/σ20K of selenium-doped rGOs (Se-rGOs) were almost constant in the temperature range of 20–300 K, indicating that the carrier mobility of Se-rGOs becomes temperature-independent after selenium doping, similar to that of pure graphene. Atomic-thin layers of selenium can turn microscale graphene oxide sheets into enhanced energy-storage devices, reports a new study. Graphene's extraordinary conductivity has attracted interest for applications such as battery anodes, but manufacturers seeking sizeable quantities often use reduced graphene oxide – chemically synthesized carbon films containing oxygen and other defects that hinder charge transport. Hyoung-Joon Jin from Inha University in South Korea and colleagues now demonstrate that heating graphene oxide with elemental selenium returns metal-like conductivity to the carbon sheet through a process called surface transfer doping. Scanning electron microscopy and first-principles calculations revealed that selenium atoms attach to surface edge defect sites and transfer electrons to graphene oxide. These dopants enabled graphene oxide to achieve a similar battery capacity and longevity as pure graphene when incorporated in a prototype lithium-ion device. Selenium atoms were selectively introduced in particular locations such as the pseudo-edge sites of hole-cluster defects in the basal plane and edge defect sites of the fully amorphorized surface of reduced graphene oxide (rGO), inducing a dramatic change of electrical transport properties of rGO by electron doping.


INTRODUCTION
Graphene has attracted much attention in physics, chemistry and materials science because of its unique characteristics, such as high electrical conductivity (~10 6 S cm − 1 ), high intrinsic mobility (200 000 cm 2 V − 1 s − 1 ), large Young's modulus (~1.0 TPa), and exceptional thermal conductivity (~5000 W m − 1 K − 1 ). 1 − 3 However, the practical application of graphene is limited because of the difficulty of its large-scale production. Therefore, a chemical route has been established for producing the so-called graphene oxide (GO). 3− 5 GO possesses numerous topological defects and oxygen heteroatoms that distinguish it from pure graphene, and it cannot be restored to its initial state even via healing by various reduction processes. The large number of sp 3 sites and residual oxygen groups of reduced GO (rGO) negatively affect the intrinsic properties of graphene. 6−8 Hence, it is indispensable to endow rGO with properties similar to those of graphene.
Chemical doping with heteroatoms or molecules that can increase the density of free charge carriers in rGO is expected to enhance its electrical conductivity. 9−13 Charge-transfer chemical doping of graphene can be subdivided into surface transfer doping and substitutional doping. 11 In general, surface transfer doping does not destroy the graphene structure and is thus reversible. In contrast, substitutional doping refers to the permanent substitution of carbon atoms in the aromatic hexagonal lattice of graphene. Surface transfer doping occurs via charge transfer from the adsorbed dopant to graphene, and is closely related to their density of states (DOS). Therefore, adsorption of electron-withdrawing molecules leads to p-type doping, while the adsorption of electron-donating molecules induces n-type doping. Similarly, the p-or n-type doping can be induced by substitution with dopants having fewer or more valence electrons, respectively. A large number of studies have reported substitutional doping by nitrogen, 9,10,14,15 boron, 14−16 sulfur 17−20 and phosphorus 20−23 in carbon-based materials. These types of dopants affect the carbon crystalline structure by acting as defects, while providing the free charge carrier density and electroactive surface properties. In addition, surface transfer doping has been performed using p-type dopants such as NO 2 , Br 2 and I 2 vapors, 24−28 and n-type dopants such as K and Rb, 27−31 leading to greatly enhanced electrical properties. Nevertheless, since these dopants exhibit air-instability and toxicity, their use in practical applications is limited.
Selenium is a polyatomic, nonmetallic chalcogen with a large atom size (located in the period 4 along with potassium and bromine) and high polarizability, and is d-electron-rich (located in group 16). 32 Because of these characteristics, selenium is a promising potential dopant for enhancing electrical transport properties of rGOs. A number of studies have recently reported the catalytic effects of selenium-doped carbon materials in the oxygen reduction, 33,34 as well as electrooxidation of methanol 35 and ethanol. 36 Although these results showed the existence of direct chemical bonding between selenium and carbon, the exact configuration of selenium and its effects on the intrinsic properties of carbon materials are still veiled. Further applications of selenium-doped carbon materials are possible only after a deeper understanding of the role of selenium in carbon host structures.
In this study, we report atomic-level selenium doping and its effects on rGOs. Selenium-doped rGOs (Se-rGOs) were prepared by heating of GO with elemental selenium powder above the sublimation temperature (4685°C) of selenium. The d-electron-rich selenium atoms were homogeneously introduced on the entire rGO surface, leading to greatly improved electrical transport properties caused by the n-type doping. As a result, the carrier mobility of rGOs was dramatically restored by the introduced selenium, reaching a value comparable to that of graphene. This study not only provides a fundamental understanding of selenium doping of rGOs but also shows their potential for energy storage applications.

Preparation of Se-rGO
GO was prepared from natural graphite (Sigma-Aldrich, St Louis, MO, USA) using the Hummers method. 18 Aqueous GO suspensions were frozen in liquid nitrogen and freeze-dried at -50°C and 0.045 mbar for 72 h using a lyophilizer (LP3, Jouan, France). After lyophilization, low-density and loosely packed GO powders were obtained. A mixture of the as-obtained lyophilized GO powder (100 mg) and selenium powder (300 mg) was heated in a tubular furnace from room temperature to 800°C at a heating rate of 10°C min -1 , and the temperature was held for 2 h. An Ar flow rate of 200 ml min -1 was used. The obtained product (Se-rGO) was stored in a vacuum oven at 30°C. rGO was prepared using the same procedure without the elemental selenium powder.

Material characterization
The morphologies of Se-rGO were observed by atomic force microscopy (AFM, NT-MDT, Russia) with an NSG-10 cantilever (NT-MDT, Russia) in semi-contact operation mode and FE-transmission electron microscopy (TEM; JEM-ARM-200F, JEOL, USA). The chemical states of Se-rGO were investigated by X-ray photoelectron spectroscopy (XPS; PHI 5700 ESCA) using monochromated Al K α radiation (hn = 1486.6 eV). Elemental analysis was performed using an elemental analysis1112 instrument (CE instrument, Italy). Raman measurements were performed using an NTEGRA Spectra spectrometer (NT-MDT, Russia) equipped with a 473 nm (2.62 eV) laser in backscattering configuration. The spectral resolution was~2 cm -1 with a 600 groove per mm grating. A × 100 objective (NA = 0.9) provided a laser spot of ca. 330 nm diameter. The laser power was kept well below 0.3 mW for nondestructive Raman measurements.

Computational details
All geometric optimizations and DOS calculations were conducted using Vienna ab initio simulation package software 37 of the density functional theory framework. We used a plane-wave basis set with a kinetic energy cutoff of 500 eV and a projector-augmented wave pseudopotential 38 provided by Vienna ab initio simulation package. A generalized gradient approximation suggested by Perdew-Burke-Ernzerhof 39 was adopted for the correction of exchange correlation energy. We constructed a model graphene structure containing 72 and 192 carbon atoms in a single layer for the calculation of the adsorption energies. To describe the edge site adsorption, a zigzag graphene nanoribbon structure was used, containing the same amount of carbon atoms as the graphene model. A 15-Å vacuum slab was introduced to avoid the interaction between graphene layers. The adsorption energy E ad was defined as where E C+Se , E C , E Se are the free energies of Se-adsorbed graphene, perfect graphene and the chemical potential of Se. For the DOS calculation, 18 k-points were uniformly distributed along the graphene nanoribbon, with the structures of the latter being fully relaxed until the system energies converge to o0.05 eV.

Characterization of electronic transport properties
Se-rGO and rGO were dispersed in DMF and deposited on a 300 nm SiO 2 /highly p-doped Si substrate. Electrode fabrication was performed by standard electron beam lithography, evaporation of Cr/Au (5/120 nm), and lift-off procedures. The temperature dependence of I-V characteristics was measured using four-probe configurations in a Janis cryogenic system with a current source (6221, Keithley) and nanovoltmeter (2182A, Keithley). The temperature dependence of the conductivity of Se-rGO and rGO was investigated between 20 and 300 K by four-probe measurements using a Keithley 6221/2182A Delta Mode System. Currents of 100 and 10 nA were applied for Se-rGO and rGO, respectively.

Electrochemical characterization
The electrochemical properties of Se-rGO and rGO as anode materials for Li-ion batteries were evaluated using a WonATech automatic battery cycler and CR2032-type coin cells. The working electrodes were prepared by mixing the active material (80 wt.%) with conductive carbon (10 wt.%) and polyvinylidene fluoride (10 wt.%) in N-methyl-2-pyrrolidone. The resulting slurries were uniformly applied onto Cu foil. The electrodes were dried at 120°C for 2 h and roll-pressed. The coin cells were assembled in a glove box filled with argon, employing a composite electrode with metallic lithium foil and 1 M LiPF 6 (Aldrich, 99.99%) dissolved in a solution of ethylene carbonate/dimethyl carbonate (1:1 v/v) as the electrolyte. The cells were galvanostatically cycled between 0.01 and 3.0 V vs Li/Li + at various current densities.
The air electrode of the Li-air battery was composed of catalyst materials (Se-rGO or rGO) and the polyvinylidene fluoride binder in an 8:2 ratio (w/w). The slurry was prepared using N-methyl-2-pyrrolidone and coated on the Ni mesh substrate. Li metal (3/8 inch diameter), a glass fiber separator (Whatman GF/D microfiber filter paper, 2.7 μm pore size), and the prepared air electrodes were assembled in a Swagelok-type cell. 1 M lithium bis(trifluoromethylsulfonyl)imide dissolved in tetraethylene glycol dimethylether was used as an electrolyte. The cells were operated at a constant outer pressure of pure oxygen at 770 Torr.

RESULTS AND DISCUSSION
Se-rGO samples have a thickness of~15 nm, a lateral size of several micrometers, and numerous ripples on their surface (Figure 1a), indicating morphologies similar to those of rGOs reduced at 800°C in the absence of elemental selenium (Supplementary Figure S1). TEM and electron energy loss spectroscopy mapping images of Se-rGOs shows that a large amount of selenium atoms were doped on the entire area of the carbon hosts structure of rGOs (Figure 1b-d). Moreover, scanning TEM (STEM) images of Se-rGOs and rGOs show clear differences, namely that Se-rGOs have a relatively rough surface, as indicated by the presence of numerous bright spots (Figure 1e and Supplementary Figure S2). The presence of selenium atoms on the rGO surface (with atomic selenium content of~13%) is also confirmed by energy-dispersive X-ray spectroscopy (EDS; Supplementary Figure S3). Nevertheless, the high-resolution TEM image (Figure 1f) of Se-rGO shows no aggregates or nanoparticles on its surface, indicating atomic-level selenium doping of rGO. Notably, the carbon structures of Se-rGOs and rGOs were highly amorphorized (Figure 1f and Supplementary Figure S4). Magnified high-resolution TEM images of Se-rGOs also show selenium atoms (Figure 1g). The large d-electron-rich selenium atoms (radius: 120 pm) are depicted as bright spots, which are observed on the highly amorphorized Se-rGO surfaces. The chemical states of Se-rGOs were characterized by XPS (Figure 2a-c). The C 1s spectrum shows several distinct peaks due to C-O and C-Se bonds centered at 285.8 eV and the ones due to C = O bonds centered at 290.5 eV. In addition, the spectrum shows the main C-C bond peak centered at 284.5 eV (Figure 2a), which is similar to the value observed for rGOs (Supplementary Figure S5a). Oxygenated functional groups were also found in the O 1s spectrum, as indicated by the peaks due to C-O and C = O bonds, centered at 533.3 and 530.5 eV, respectively (Figure 2b). The C/O atomic ratio of Se-rGOs was~10.9, which was higher than the corresponding value of 9.2 for rGOs, indicating the more reduced nature of Se-rGOs. Selenium atoms were mostly present in the form of C-Se-C bonds, as indicated by the peak centered at 56.3 eV, with peak due to the minor Se-O bonds located at 59.3 eV (Figure 2c), showing that selenium atoms were chemically doped on the rGO surface. The C/Se atomic ratio was calculated as~40.7, which is much lower than the value suggested by the energy-dispersive X-ray spectroscopy data. To determine the selenium content more accurately, elemental analysis was performed, and we obtained~15 wt.% selenium content for Se-rGOs, in agreement with the XPS results. In addition, Fourier transform infrared spectroscopy data supports the presence of the C-Se-C bonds (Supplementary Figure S5b). Nevertheless, it is notable that the Raman spectra showed similar D/G intensity ratios (I D /I G ), as well as integrated sp 3 and sp 2 area ratios (A sp3 /A sp2 ) for Se-rGO and rGO (Figure 2d, Supplementary Figures S6 and S7). The deconvoluted Raman spectra of Se-rGO and rGO and their integrated area ratios of sp 3 and sp 2 (A sp3 /A sp2 ) are shown in Supplementary Figures S6 and  S7. 40 The D band corresponds to the disorder in the A 1g breathing mode of the six-membered aromatic ring near the basal edge, while the G band reflects the hexagonal structure related to the E 2g vibration mode of the sp 2 -hybridized carbon atoms. Therefore, the similar I D /I G values mean that selenium doping slightly affects the hexagonal carbon structure, suggesting that the doping could be close to surface transfer doping.
More specific information on chemical bonding was obtained using density functional theory calculations (Figure 3). Various defect sites of rGOs including intrinsic topological defects, pseudo-edge sites of hole-cluster defects in the basal plane of graphene, and graphene edge defects are considered possible adsorption sites for selenium atoms (Figure 3a). As shown in Figure 3a, configurations 1-5 represent selenium adsorption on top of the carbon atoms (1), C-C bonds (2), carbon rings (3), at the pseudo-edge sites of hole-cluster defects in the graphene basal plane (4), and at the graphene edge defect sites (5; Supplementary Figure S8). The calculated selenium adsorption energies for each adsorption site are given in Figure 3b. Note that multiple points are described for configurations 1-3, since there are several distinctive local environments for selenium adsorption, composed of 5-, 6-and 7-membered carbon rings (see Supporting Information for details). The most stable adsorption sites in configurations 1-3, denoted by red circles in Figure 3b, were the carbon atoms 577 of configuration 1, 77 of configuration 2, and the fivemembered carbon ring of configuration 3. However, all selenium adsorption energies in configurations 1-3 were positive, indicating that binding of Se atoms at intrinsic topological defect sites such as monovacancies, di-vacancies and Stone-Wales defects is unfavorable. These results agree with the previous computational reports on substitutional sulfur doping of graphene 41 and surface sulfur doping of defected graphene, 18 both showing a positive binding energy for sulfur adsorption at the intrinsic defect sites of graphene. On the other hand, the binding energies of selenium atoms at both pseudo-edge hole-cluster defect sites in the basal plane and the edge defect sites of graphene are negative, indicating that Se atoms can be spontaneously bound only at the pseudo-edge and edge defect sites, which are relatively unstable. Moreover, the negative energy of selenium adsorption implies chemical bonding with graphene on an atomic level, not just physical adsorption. We also considered the possibility of selenium adsorption in the presence of oxygenated functional groups (Supplementary Figure S9). We find that selenium atoms are not likely to directly bind with oxygen; however, their individual adsorption is possible in near oxygenated functional groups, in agreement with the XPS result that showed Se-O bonding (Figure 2b). In addition, we calculated the DOS for H-and Se-terminated graphene nanoribbons to investigate the effect of Se doping on electronic properties (Figure 3c and d). The detailed model structures used for DOS calculations are depicted in Supplementary Figure S8e. As shown in Figure 3c and d, the Fermi level of Se-terminated graphene nanoribbons is 0.5 eV higher than that of H-terminated ones. As selenium has more valence electrons than carbon atoms, they are naturally transferred from selenium to rGO, hybridizing with the carbon atoms of rGO and occupying high-energy states, which leads to a positive shift of the Fermi level and makes Se-rGO an n-type conducting material. This effect is also reported for doping with K 42 and NH. 43 Since n-type doped graphene typically shows increased electrical conductivity due to the introduction of charge carriers 43 Se-rGOs can exhibit improved electrical conductivity compared to rGOs.
Temperature-dependent current-voltage (I-V) characterization of rGO and Se-rGO was performed (Figure 4a and b). For both samples, the I-V curves show highly symmetric and linear behavior in all temperature regions, with the slope of I-V curves increasing concomitantly with temperature. The curves indicate an increase in conductivity with temperature. Interestingly, the current of Se-rGO is two orders of magnitude higher than that of rGO, and the I-V curve slope for Se-rGO changes little, despite the temperature variation from 20 to 300 K. For detailed observations, the temperature-dependent conductivities (σ(T)) of rGO and Se-rGO were obtained as the temperature was continuously varied from 20 to 300 K (Figure 4c). The σ(T) of Se-rGO is larger than that of rGO by 1-2 orders of magnitude. The defects and functional groups cause an increase of dipole polarization, thus, reducing the conduction path length in GO. 44 For rGO, the dipole polarization is reduced by the decreased amount of defects and oxygenated functional groups, which increase Selenium-doped reduced graphene oxide YS Yun et al the probability of electron hopping. Hence, the electronic transport mechanism of rGO has been attributed to variable hopping. 45 The rGO follows the Efros-Shklovskii variable-range hopping (σ = σ 0 exp (-T 0 /T) 1/2 , where σ 0 and T 0 are 87.8 S cm -1 and 508.2 K, respectively), which is known for the transport properties of rGO with different fractions of sp 2 carbon (blue color in Figure 4c). 46 However, the electrical transport properties of Se-rGO were not fitted by the Efros-Shklovskii variable-range hopping model, but were well explained by the fluctuation-induced tunneling mechanism 47 (red color in Figure 4c). The σ(T) in the fluctuation-induced tunneling model is expressed as σ = σ ∞ exp(-T 1 /(T+T 0 )), where T 0 and T 1 are 200 and 96 K, respectively. T 0 is the onset temperature of thermally activated conduction, and T 1 is a measure of the energy required for an electron to cross the insulator gap between conductive regions. These values indicate that Se-rGO possesses large metallic regions separated by a relatively small insulating gap. These metallic regions created by selenium doping increase the conductivity to 210 S cm -1 at 300 K, a very large conductivity value compared to that of rGO (obtained using N 2 H 4 ,~2 S cm -1 ) 45 and is comparable to the ultrahigh bulk conductivity of GO papers reduced by metal iodides. 48 Note that the variation of σ(T) in Se-rGO is quite small (Figure 4d). The σ/σ 20 K value of rGO increases 41.5-fold as the temperature increases from 20 to 300 K, while that of Se-rGO stays almost constant (1.3-fold increase from 20 to 300 K). This means that the carrier mobility becomes temperature-independent because of selenium doping; the carrier mobility being similar to that of pure graphene. 49 Furthermore, the effects of selenium doping on the electrical transport properties of monolayer rGO were directly investigated by in situ measurements at 300 K (Supplementary Figure S10). The sample was loaded in a doping tube, with selenium powder placed on the bottom. The tube containing the sample was evacuated using a rotary pump (o3 × 10 -2 Torr; Supplementary Figure S10a). As the bottom of the tube was heated, the conductivity of rGO abruptly increased from 0.27 to 0.58 S cm -1 and then saturated (Supplementary Figure S10b). Moreover, we also studied the gate voltage (V G )dependent conductivity at 300 K, observing graphene-like field-effect behavior. The minimum conductivity (σ min ) of rGO was located near V G = 60 V (Supplementary Figure S10c). However, this value shifted toward negative V G after selenium doping, implying that this doping was n-type. Based on the electrical transport properties, we conclude that the chemisorbed selenium atoms on rGO result in electron doping of the latter, enhancing its electrical conductivity.
The superior electrical transport properties of Se-rGOs lead to significantly enhanced electrochemical properties of energy storage systems. As an anode material for Li-ion batteries, Se-rGOs show a much higher rate performance as well as an increased capacity, as compared to rGO (Figure 5a-c). The reversible capacity of Se-rGO equals 570 mAh g -1 at a current density of 1 C (372 mA g -1 ; Figure 5a), which is much higher than the value of 310 mAh g -1 , observed for rGO (Figure 5a). The enhanced capacity of Se-rGO suggests that selenium atoms contribute to the reversible Li-ion storage, as supported by the change of galvanostatic discharge/charge Selenium-doped reduced graphene oxide YS Yun et al profiles. The same test conducted at low temperature (0°C) shows a more distinct capacity gap between Se-rGO and rGO, as displayed in Supplementary Figure S11. In the low-temperature cell test, the reversible capacity of Se-rGO is 360 mAh g -1 , corresponding tõ 63% of the capacity at 25°C. In contrast, rGO shows a reversible capacity of~140 mAh g -1 , corresponding to~45% of that at 25°C. These results suggest that the enhanced electrical properties of Se-rGO affect its electrochemical performance at low temperature to a great extent. Highly stable capacities of 165 mAh g -1 could be obtained at a large current rate of 50 C, which is approximately seven times higher than the value of rGO anode materials (Figure 5b). In addition, the high rate performance of Se-rGO is significant compared to graphenebased electrode materials and nitrogen-doped graphene nanosheets for Li-ion batteries. 50,51 After 200 cycles with a consecutive increase of current density from 1 C (372 mA g -1 ) to 50 C (18,600 mA g -1 ), the specific capacity of Se-rGO returns to its initial value of ca. 460 mAh g -1 , indicating good reversibility. In addition, cycling stability is found to be maintained during 1000 repetitive cycles, showing a Coulombic efficiency of nearly 100% (Figure 5c).
We thought that Se-rGO could be a potential catalyst material for Li-O 2 batteries because of its high-energy density. However, we found that such batteries suffered from large polarization during charging, resulting in poor cycling stability (Figure 5d and e). The above behavior was due to the formation of non-conductive Li 2 O 2 during discharge. When a conventional ether-based electrolyte (that is, tetraethylene glycol dimethylether) was used for Li-O 2 batteries, amorphous Li 2 O 2 entangled on the carbon surface and bulk crystalline Li 2 O 2 with toroidal morphology were both observed simultaneously. The preference between amorphous and crystalline Li 2 O 2 is significantly affected by the properties of the electrolyte, current rates, and depth of discharge. 52,53 In the fully discharged state, as shown in Figure 5d, nucleated Li 2 O 2 has enough time to grow into large toroidal particles with high crystallinity. Based on theoretical studies of the charge transport properties of bulk Li 2 O 2 , its intrinsic ionic and electronic conductivities were calculated as~10 -19 S cm -1 , which is unacceptably low for charge conduction. 54 Therefore, despite the enhanced electrical conductivity, the charge profile of Se-rGO showed a large polarization similar to that of rGO because of the low conductivity of the large toroidal Li 2 O 2 particles. In this respect, a particular electrochemical protocol should be used to clearly demonstrate the effect of Se-rGO in Li-O 2 cells. Thus, we used capacity-limited discharge to produce surface-derived Li 2 O 2 . Since the surfacederived Li 2 O 2 directly comes into contact with the Se-rGO surface, the enhanced electrical properties of Se-rGO could highly affect the charge polarization of Li-O 2 cells. As shown in Figure 5e, the air cells tested under a protocol, wherein the discharge capacity was limited to 500 mA h g À1 , show an obvious difference in their charge profiles. Overall, the Se-rGO-based electrode showed a lower charge potential than rGO; in particular, the initial part (~200 mA h g -1 ) had a charge potential below 4.0 V. Based on this result, we could verify the effect of selenium doping in rGO. Additionally, when a soluble LiI catalyst was used, the Se-rGObased electrode showed an abrupt decrease of charge polarization near 0.5 V (Figure 5e), which was superior to the previous reports on rGO electrodes in terms of energy efficiency. [55][56][57] These results, thus, demonstrate the possibility of using Se-rGO-based electrodes for Li-O 2 batteries.

CONCLUSION
In this work, we investigated the effects of selenium doping in rGO. Numerous selenium atoms were homogeneously doped on the top surface of rGOs; it caused little structural change, indicative of surface transfer doping. high-resolution TEM images showed no aggregates or nanoparticles on the surface of rGOs, indicating that selenium atoms were introduced in an atomic level. Based on firstprinciples calculations, we found the selenium atoms to be mainly bonded to the surface edge defect sites of highly amorphorized rGO, leading to the enhancement of electrical conductivity due to n-type doping. In contrast, the overall intrinsic topological defects of the graphene basal plane showed an unfavorable selenium atom binding energy. The doped selenium atoms changed the electrical transport properties of rGO. The properties were well fitted by the fluctuation-induced tunneling mechanism, and the results indicated the presence of a large metallic region separated by a relatively small insulating gap. The σ value of Se-rGO reached 210 S cm -1 at 300 K, while the σ/σ 20 K parameter was almost constant (1.3-fold increase from 20 to 300 K), this indicated that the carrier mobility of Se-rGO was similar to that of graphene. In addition, Se-rGO showed a superior electrochemical performance as an anode for Li-ion batteries and a catalyst for Li-O 2 batteries because of the effects of selenium doping in rGO. Figure 4 Temperature-dependent electron transport properties of rGO and Se-rGO. Temperature-dependent electrical transport properties of (a) rGO and (b) Se-rGO measured between 20 and 300 K in 20-K steps with inset optical images of the four-electrode configurations. The slope of the I-V curves increased concomitantly with temperature, indicating an increase of conductivity with temperature. (c) Temperature dependence of rGO and Se-rGO conductivities obtained from I-V curves and variable-range hopping (VRH) and fluctuation-induced tunneling (FIT) models. rGO and Se-rGO were well-fitted by VRH and FIT models, respectively. (d) Temperature vs σ/σ 20 K plots for rGO and Se-rGO. The above curve for Se-rGO is almost flat in the temperature range from 20 to 300 K (1.3-fold increase), indicating that the carrier mobility becomes temperature-independent due to selenium doping, implying mobility restoration. Figure 5 Electrochemical performances of rGO and Se-rGO as the anode and catalyst for Li-ion and Li-air batteries, respectively. (a) Galvanostatic discharge/ charge profiles of Se-rGO and rGO at a current density of 1 C (372 mA g -1 ), (b) rate performances of Se-rGO and rGO at consecutively increasing current densities from 1 C (372 mA g -1 ) to 50 C (18 600 mA g -1 ), with subsequent reversal to 1 C (372 mA g -1 ), and (c) cycling performance of Se-rGO as a Li-ion battery anode for 1000 repetitive cycles at a current density of 3 C (1116 mA g -1 ).