Next Article in Journal
MEG Effects on Hydrolysis of Polyamide 66/Glass Fiber Composites and Mechanical Property Changes
Next Article in Special Issue
A Molecular Electron Density Theory Study of the Chemoselectivity, Regioselectivity, and Diastereofacial Selectivity in the Synthesis of an Anticancer Spiroisoxazoline derived from α-Santonin
Previous Article in Journal
Transformation Products of Organic Contaminants and Residues—Overview of Current Simulation Methods
Previous Article in Special Issue
Understanding the Molecular Mechanism of the Rearrangement of Internal Nitronic Ester into Nitronorbornene in Light of the MEDT Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Li and Na Adsorption on Graphene and Graphene Oxide Examined by Density Functional Theory, Quantum Theory of Atoms in Molecules, and Electron Localization Function

1
Department of Physics and Astronomy, University of Texas Rio Grande Valley, Edinburg, TX 78539, USA
2
Achieve Early College High School, McAllen, TX 78501, USA
3
PSJA Thomas Jefferson T-STEM Early College HS, Pharr, TX 78577, USA
4
The Science Academy, Mercedes, TX 78570, USA
5
South Texas High School for Health Professions, Mercedes, TX 78570, USA
6
Department of Physics and Astronomy, University of Texas Rio Grande Valley, Edinburg, TX 78539, USA
*
Author to whom correspondence should be addressed.
Molecules 2019, 24(4), 754; https://doi.org/10.3390/molecules24040754
Submission received: 14 January 2019 / Revised: 8 February 2019 / Accepted: 12 February 2019 / Published: 19 February 2019
(This article belongs to the Special Issue The Molecular Electron Density Theory in Organic Chemistry)

Abstract

:
Adsorption of Li and Na on pristine and defective graphene and graphene oxide (GO) is studied using density functional theory (DFT) structural and electronic calculations, quantum theory of atoms in molecules (QTAIM), and electron localization function (ELF) analyses. DFT calculations show that Li and Na adsorptions on pristine graphene are not stable at all metal coverages examined here. However, the presence of defects on graphene support stabilizes both Li and Na adsorptions. Increased Li and Na coverages cause metal nucleation and weaken adsorption. Defective graphene is associated with the presence of band gaps and, thus, Li and Na adsorptions can be used to tune these gaps. Electronic calculations show that Li– and Na–graphene interactions are Coulombic: as Li and Na coverages increase, the metal valences partially hybridize with the graphene bands and weaken metal–graphene support interactions. However, for Li adsorption on single vacancy graphene, QTAIM, ELF, and overlap populations calculations show that the Li-C bond has some covalent character. The Li and Na adsorptions on GO are significantly stronger than on graphene and strengthen upon increased coverages. This is due to Li and Na forming bonds with both carbon and oxygen GO atoms. QTAIM and ELF are used to analyze the metal–C and metal–metal bonds (when metal nucleation is present). The Li and Na clusters may contain both covalent and metallic intra metal–metal bonds: This effect is related to the adsorption support selection. ELF bifurcation diagrams show individual metal–C and metal–metal interactions, as Li and Na are adsorbed on graphene and GO, at the metal coverages examined here.

Graphical Abstract

1. Introduction

Rechargeable lithium ion batteries (LIB) serve as energy storage devices for a variety of applications, such as portable electronics and hybrid and electric vehicles [1,2,3,4]. The LIB performance depends on the structure and composition of the anode and cathode electrodes [5]. Graphite has been used as anode material for LIB [6], with graphite’s theoretical capacity limit set at 372 mAh/g (intercalation compounds of LiC6 configurations) [7]. Current and future needs dictate the search for LIB materials with capacities significantly higher than graphite to satisfy portable electronic devices’ requirements [5]. The use of other carbon materials, such as graphene [8,9], C60, and carbon nanotubes, increases the LIB capacity to 540 mAh/g, 730 mAh/g, and 784 mAh/g, respectively [10]. Graphene is a monolayer honeycomb lattice of carbon atoms and is known for its high mechanical strength and high thermal and electrical conductivity [11,12,13]. Graphene B and N-doping further increase the LIB capacity to 1227 mAh/g and 872mAh/g, respectively [14].
Increased LIB graphene capacities relative to graphite seem indicative of increased Li atoms being adsorbed on graphene relative to graphite. However, the contrary has been observed via in situ Raman spectroscopy [15]. Although computational simulations support the adsorption of single Li atoms on graphene (one side [16,17,18] and both side adsorption [5,19]), this adsorption is less stable relative to bulk metallic Li. Lee et al. showed that Li adsorption on pristine graphene is not favorable relative to bulk metallic Li under any configuration and is only possible on graphene defects and edges [20]. Okamoto further supported Li adsorption on defects relative to edges [21]. This is due to the edges of the anode LIB materials being covered by low electron conductivity solid electrolyte interphase, thus affecting the reversibility of Li movement during charge and discharge processes. Li adsorption on defective graphene has been reported using computational calculations [21,22,23,24]. The presence of double carbon vacancies on N-doped graphene support can also improve LIB performance [25].
The main disadvantage of LIBs is Li cost, which has been rising steadily in recent years, with the exception of a price drop observed in the last quarter of 2018 [26]. Li is a very rare element (Li concentration is ~35 ppm in the upper continental crust [27]), and thus there are concerns that Li may not be sufficient for meeting future energy demands. LIB, as well as sodium ion batteries (NIB) [26,28], suffer from dendrite growth due to metal nucleation (i.e., clustering), which lowers their performance [1,29]. Li forms clusters due to the small Li diffusion barrier on graphene support [30]. Li cluster formation on pristine graphene (i.e., graphene without defects) and defective graphene [21] has been explored using density functional theory (DFT) [30,31]. Defective graphene is produced during graphene oxide’s (GO) thermal reduction process, as hydroxyl and epoxide O atoms are removed from the GO sheet, via oxidation and reduction [32,33]. Li6 clusters collapse during adsorption on 24 vacancies defective graphene [21].
Sodium ion batteries could serve as alternative batteries due to the abundance of Na. However, the NIB capacity of graphite is only 35 mAh/g, substantially lower relative to LIB [26]. This is attributed to NaC64 formation [34,35] due to the large Na+ size relative to Li+ (100 pm vs. 76 pm) and Na high ionization potential [21]. NIB capacities were recorded at 1450 mAh/g, for Na adsorption on defective graphene with divacancy [24]. DFT calculations on Na adsorbed on defective graphene have been reported [36]. Specifically, Na cluster formation for adsorption on pristine and defective graphene has been explored via DFT calculations by Liang et al. [37].
In this work, we examine Li and Na adsorptions on pristine and defective graphene, as well as on graphene oxide (GO), at various metal coverages, using periodic DFT, with and without long-range interactions. The metal-support interactions are elucidated the via changes in the adsorption and cohesive energies (Eads and Ec, respectively) per metal, the metal–Carbon plane distances (height, h), and the electronic charge transfers from the metals to graphene and GO supports. The Eads and the E c for the adsorption of n-metals (n = 1, 2, …) on graphene and GO are defined as
E a d s = ( E n · X / s u p p o r t n · E X E s u p p o r t ) / n
E c = ( E n · X / s u p p o r t E n · X E s u p p o r t ) / n
where X is the metal adsorbate and the support is either graphene or GO. The Eads is the minimum energy required to remove a single Li or Na atom during adsorption on graphene and GO, whereas Ec measures the cluster-support interaction. For n = 1, Eads = Ec. Single metals and their clusters form a bond with the support if Ec < 0. The metal-support interaction is considered stable if the |Eads| is larger than the metal binding energy of its bulk. The metal-support interactions are also examined via electronic structure analyses.
The local metal–C and the metal–metal interactions are examined by two independent analyses: The quantum theory of atoms in molecules (QTAIM) and the electron localization function (ELF). QTAIM was developed by Bader and co-workers, who note that orbitals, much less wave functions, are unphysical in nature; thus, chemical bonds are not quantum observables [38,39]. However, the electron densities ( ρ r ) and their gradients are quantum observables and can be used to characterize chemical systems [39,40]. ELF [41] is a measure of the probability of finding an electron in the vicinity of another same-spin electron. It reveals core, binding, and lone electron pairs in molecular and crystalline systems. Results from both QTAIM and ELF analysis are basis-set and method-independent, as long as a minimally adequate basis set is used in these calculations [42,43].

2. Computational Methods

2.1. Periodic Structure Modeling

Monolayer pristine graphene is modeled as a two-dimensional 4 × 4 hexagonal lattice with 32 carbon atoms in the unit cell. Figure 1 shows the DFT optimized structures for Lin and Nan (n = 1, 3, 5) adsorbed on pristine graphene, the Li and Na Ec values, the heights h for the closest metal to the graphene plane, and the smaller of the metal–metal distances per configuration. The Lin/graphene and Nan/graphene (n =1, 3, 5) unit cells are built as follows: for n = 1, single metal atoms (i.e., Li and Na) are initially placed distant from graphene, followed by geometry optimizations; for n = 3, two metal atoms are placed distant from the previously optimized Li/graphene and Na/graphene structures, followed by geometry optimizations. A similar approach is followed for the Li5/graphene and Na5/graphene cases. Figure 2 shows the optimized structures for Lin and Nan (n =1, 3, 5) adsorbed on single and double vacancy defective graphene, as well as the Ec and geometrical parameters, in a similar fashion as in Figure 1. Graphene defects are modeled by the use of “ghost” massless C atoms in the graphene lattice: These atoms share the same basis set as C, in agreement with our past report [44]. Here, the ghost atoms are placed in the center of the supercell and the defect is periodically repeated by the use of translational symmetry vectors [45]. Ghost atoms remain fixed in the locations occupied by the original C atoms. Figure 3 shows the optimized structures for Lin and Nan (n = 1, 3) adsorbed on graphene oxide (GO, O/C = 43.75%), with Ec and geometrical parameters. Here, GO supports contain epoxides and hydroxyls at 25% and 18.75%, respectively (percentages relative to C). Our GO structure is in agreement with the proposed non-hydrated GO structure from Zhou and Bongiorno [46] (26% and 18% for epoxides and hydroxyls, respectively), found by using X-ray photoelectron spectroscopy and DFT.

2.2. DFT Parameters

CRYSTAL17 [47] is used to calculate the optimized geometries and electronic properties for Li and Na atoms adsorbed on graphene (with and without vacancies) and GO. An advantage of CRYSTAL is its treatment on two-dimensional slabs, without the need for artificial periodicity perpendicular to the surface, as is typically found in other DFT codes [48]. Periodic restricted DFT is used here, as in our past work [17,44,49,50,51]. The PBE0 non-empirical/parameter-free functional is used [52,53]. All calculations are treated using the D2 semiempirical correction by Grimme [54]. This correction improves the DFT functionals description of the long-range electron correlations, which are responsible for van der Waals interactions between the adsorbed metals and graphene/GO supports. Moreover, for selective cases, our calculations are repeated without employing D2 corrections. All atoms are described by all-electron basis sets. The C and H atoms are treated by the 3-ζ valence plus polarization basis set (pob-TZVP) from Peintinger et al. [55] and the O atoms by the split-valence 8-411G(2d1f) basis set from Mahmoud et al. [56]. These basis sets are optimized for crystalline calculations. The Li and Na are described by the atomic basis sets (6s3p1d) and (9s5p1d), respectively, and were used in our past work [17]. Calculations are at 0 K, as expected. Temperature and pressure effects will be described in a future work.
For geometry optimizations, Brillouin zone integrations (Monkhorst-Pack grid) [57], the Fermi energy, and the density matrix calculations (Gilat grid) [58,59] are performed on a 12 × 12 grid. The densities of states (DOS) and the electronic band structures are obtained using the denser 48 × 48 grid. The band structures are calculated on the reciprocal space Brillouin zone path Γ-Μ-Κ-M´-Γ, where Γ = (0, 0, 0), Μ = (1/2, 0, 0), K= (2/3, 1/3, 0), and M´ = (0, 1/2, 0). The Fermi surface is smeared with a Gaussian of 0.005 Hartrees for convergence purposes. Moreover, the SCF energy convergence is achieved using Anderson quadratic mixing [60], coupled with additional mixing of the occupied with the virtual orbitals. The SCF energy threshold value for our calculations is set at 10-9 Hartrees (default value is 10-7 Hartrees). A large integration grid is used: this is a pruned grid with 75 radial and 974 angular points. The Crystal Orbital Overlap Populations (COOP) [61], for selected pairs of atoms, are calculated directly by CRYSTAL17. Charge transfers are calculated using Mulliken population analysis [62] and the iterative Hirshfeld population [63], the latter being based on the atoms in molecules (AIM) approach [38,39]. Hirshfeld populations are less basis-set depended than Mulliken [64] and smaller in magnitude (relative to Mulliken) [65]. The absolute magnitude of the calculated charges obtained through populations analysis has little physical meaning: These calculations may depend on the method and basis sets used [65], and thus we only focus on the relative changes of these values.
The stabilities of the final geometry conformations are secured via post-geometry optimizations. When Gaussian basis sets are used, Eads and Ec calculations suffer from basis set superposition errors (BSSE) [66]. The BSSE errors are minimized using large basis sets and employing counterpoise corrections using “ghost” massless atoms in the fragment energy calculations of the metal–graphene structures. Here, we only report BSSE corrected Eads and Ec values. Calculated Eads and Ec values using the counterpoise correction are less negative in energy than corresponding uncorrected values [17,49,67].

2.3. QTAIM Parameters

The ρ r topology is performed by the TOPOND program [68], which is incorporated in CRYSTAL17. TOPOND calculates ρ r and its Laplacian ( 2 ρ r ) at critical points, for periodic and non-periodic systems. The critical points in the ρ r and 2 ρ r topological analyses fulfill the conditions ρ r r c p =0 and 2 ρ r r c p =0, respectively, where r c p is the vector from the atom to the critical point location. The bond critical points are denoted as (3, -1), where 3 stands for the three non-zero eigenvalues of the ρ r Hessian at r c p and −1 is the algebraic summation of the number of positive and negative eigenvalues. Bonds strengths are better described by (H/ ρ ) r relative to ρ r and 2 ρ r , at the bond critical point r c p , where H r = G r + V r , G r is the positive definite kinetic energy density and V r is the potential energy density [69]. Specifically, for covalent bonds, the bond strength and (H/ ρ ) r , at the bond critical point, are positively correlated, whereas the opposite is observed for ionic bonds. Moreover, for covalent and ionic bonds, at the bond critical point r c p , (H/ ρ ) r and ρ r are positively and negatively correlated, respectively [70,71]. Gatti stated that adoption of a single QTAIM criterion for bond assessment is challenging [70]. Therefore, we used ρ r , 2 ρ r , and (H/ ρ ) r at metal–C and metal–metal bond critical points r c p , to determine changes in the corresponding bond strengths via QTAIM.
Ionic, covalent, and metallic bonds are classified via changes in ρ r , 2 ρ r , (H/ ρ ) r , and (G/ ρ ) r , at bond critical points r c p [69]. Ionic bonds, such as Li-C and Na-C (vide infra), satisfy the conditions 2 ρ r c p > 0, (H/ ρ ) r c p > 0 , and (G/ ρ ) r c p 1 , whereas ρ r c p is small (at Li-C and Na-C bond critical points). Here, the metal–metal bonds appear as non-polar covalent or metallic bonds (vide infra). For non-polar covalent bonds, 2 ρ r c p 0, (H/ ρ ) r c p 0 , and (G/ ρ ) r < 1 . The C-C bonds have large ρ r c p values, whereas covalent Li-Li bonds, found in small Li clusters, have small ρ r c p values [72]. For metallic bonds, the above conditions are modified as follows: ρ r c p is small, 2 ρ r c p ~ 0, and (H/ ρ ) r c p < 0 . Metallicity ξ r = ρ r 2 ρ r (units of length2) serves as an additional criterion to classify non-polar covalent and metallic bonds [73,74]. For metallic bonds, ξ r > 1 , whereas negative metallicity values, due to 2 ρ r < 0 , are indicative of non-polar covalent bonds [73].

2.4. ELF Parameters

The electron localization function (ELF) is an alternative to QTAIM for analyzing chemical bonding [43,75,76,77]. ELF measures the electron localization in atomic and molecular systems and it was originally described by Becke and Edgecombe [41] as follows:
η r = 1 / 1 + χ r 2
where χ r is the relevant ELF kernel and 0 η r 1 [70]. The η r = 1 case corresponds to perfect localization, whereas the η r = 0.5 to the electron gas like pair probability. The ELF function η r reveals atomic shell structure and core, similar to 2 ρ r . However, 2 ρ r fails to reveal the shells structure for heavy atoms [78]. Kohout and Savin found that η r can resolve all shells for all atoms from Li to Sr [79]. Savin et al. [80] interpreted ELF as a local measure of the Pauli repulsion on the kinetic energy ρ r , and applied ELF to a variety of atoms, molecules, and solids. Specifically, η r   ~   1 is indicative of weaker Pauli repulsion than in the uniform electron gas. ELF gradient analysis reveals attractors (local maxima) and corresponding basins [43]. In QTAIM, basins are regions in space bounded by zero flux surfaces in the gradient of ρ r [81]. Basins are classified as core and valance and characterized by their synaptic order (i.e., the number of core basins that share a common boundary [82]). Only core basins contain a nucleus, whereas valence basins are always connected to at least one core basin. At low η r values only one surface is visualized, whereas at higher η r basins separate. Therefore, bifurcation (i.e., tree) diagrams [83], with the bifurcation points corresponding to η r minima, can be used to examine bond types and regions in space, where lone pair electrons are located. The first bifurcation point corresponds the separation of the core and the valance basins from the same atom. The lower the bifurcation point the more localized are the corresponding basins [43]. The ELF 3D isosurfaces are plotted using Chimera [84].

3. Results and Discussion

3.1. Structural Properties, Eads, and Ec for Li and Na Adsorbed on Graphene and Graphene Oxide

3.1.1. Adsorption on Pristine Graphene

At 1/32 monolayer (ML) coverage, Li and Na are adsorbed on the hollow site, in agreement with past reports [5,16,17,18,22]. At the higher coverages 3/32 ML and 5/32 ML, both metals form clusters (Li3, Li5, Na3, and Na5, Figure 1b,c,e,f), as expected [30,31,37]. Here, the effect of D2 corrections in our calculated Ec values is small. However, disagreements between the DFT calculations with and without D2 corrections are observed in the h values and the metal–metal closest distances (e.g., for Na5/graphene, |Δh| = 0.11 Å, Figure 1f). For Li3 and Li5 adsorptions on pristine graphene, our closest calculated Li-Li distances are 2.85 Å and 2.77 Å, when D2 corrections are included in our calculations (Figure 1b,c). These values are smaller than the bulk Li closest Li-Li distance of 3.04 Å [85]. Similarly, for the Na case, our closest Na-Na distances during adsorption on pristine graphene are less than the bulk Na closest Na-Na distance of 3.72 Å [86]. These are indicative of the presence of stronger Li and Na intra-cluster bonds relative to the ones in their metallic bulks.
Figure 4 shows the Eads for Li and Na adsorptions on graphene and GO. At the same metal coverages, Li adsorption on pristine graphene is significantly stronger relative to Na, as shown by Eads, Ec, and h values. This statement is extended for adsorption on defective graphene and GO, with the exception of Na3/GO calculations without D2 corrections (vide infra). Our Li and Na |Eads| values are smaller than the binding energies of bulk Li (1.63 eV/Li [87]) and Na (1.13 eV/Na [88]): These are indicative of unstable Li and Na adsorptions on pristine graphene, at the coverages of this work, in agreement Lee at al. [20] and Liang et al. [37]. Although our calculations are at 0 K, the above metal stability statements are in agreement with experimental observations at room temperature. For example, Na adsorption on pristine graphene has only been detected under Na/C64 configuration [89], where in this case Na |Eads| is larger than its bulk binding energy [37]. For adsorption on pristine graphene, the Ec becomes less negative along with increased metal coverage, thus weakening the metal-support interaction, as expected [17,37]. This statement is extended for adsorption on defective graphene. Table 1 shows our Eads values and from past reports, for Li and Na adsorbed on pristine and single vacancy graphene: our Eads values are in good agreement with past work.

3.1.2. Adsorption on Defective Graphene

Li and Na adsorptions on defective graphene are stable (Figure 4). Here, Li and Na |Eads| values at all coverages are larger than the binding energies of their metallic bulk, in agreement with past reports [20,37]. The presence of defects on graphene significantly strengthens Li and Na adsorption relative to their adsorption on pristine graphene (Figure 1 and Figure 2). Increased carbon vacancies decrease the metal-graphene interaction (Figure 2). For example, for Li/graphene at 1/32 ML coverage, increased C vacancies decrease |Ec| from 3.44 eV to 2.14 eV. At 3/32 ML coverage, the presence of a single C vacancy causes the Li clusters to collapse, whereas in all other cases Li and Na form clusters, in a similar fashion as during adsorption on pristine graphene.
The strong adsorption of Na on defective graphene (and GO, vide infra) is indicative of its possible use in NIB. Although, Na adsorption on graphene is not as strong as Li, its abundancy is a significant advantage over Li.

3.1.3. Adsorption on Graphene Oxide

Li and Na adsorptions on GO are also stable (Figure 4). Moreover, contrary to adsorption on pristine and defective graphene, increased metal coverage strengthens the metal-GO interaction. This statement does not apply for the Na3/GO geometry calculated without D2 corrections. Figure 3d,e show striking differences between the optimized geometries for Na3/GO with and without D2 corrections. Therefore, for Na adsorption on GO, at coverages > 1/32 ML, it is important to include long-range interactions due to Na atoms’ large polarizations [37]. For Na3/GO, Na clustering is observed, in contrast to the Li3 case.

3.2. Charge transfers and COOP

Table 2 shows the calculated Li and Na Hirshfeld and Mulliken charges and the metal–C COOP values for the smaller metal–C distances per configuration, for Li and Na absorbed on graphene and GO, under the metal coverages examined here. Chan et al. stated that alkalis are very close to ideal ionic bonding and, thus, Li-C and Na-C bonds are expected to be ionic [16]. There is a positive relationship between the average per metal charge transferred to the support and the metal–C bond strength. Both Hirshfeld and Mulliken calculations show that charge is transferred from Li and Na towards the support (graphene or GO). This charge transfer is decreased along with increased metal coverage, in agreement with decreases in the |Ec| and, thus, the weakening of the metal-support interaction.
Dobrota et al. showed that Na adsorption on graphene-OH compounds leads to Na being neutral [96]. This is in contrast with our work, which shows that charge is transferred from Li and Na to GO (Table 2). However, the support of Dobrota et al. represents a reduced GO structure (rGO) with a significantly smaller O/C ratio than in our case.
The metal–C bond strengths can be estimated via changes in their COOP values. The majority of the metal–C COOP values in Table 1 are of the order expected for ionic bonding. However, there are some exceptions: The Li-C COOP value for Li adsorption on single C vacancy graphene (i.e., 0.138) is at the order expected for covalent bonding. Doan et al., using COOP, found that Li bonds covalently to the Nafion SO3 group [97]. Our COOP calculations indicate that metal–C bonds, for adsorption on single vacancy graphene, are covalent (Table 1). COOP calculations are basis set dependent [98]. Therefore, these bonds will be also examined via the method and basis set independent QTAIM and ELF analyses in the next sections.

3.3. DOS and electronic structures

Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 show the DOS and band structures for the configurations of Figure 1, Figure 2 and Figure 3, calculated with the inclusion of D2 corrections.

3.3.1. Adsorption on Pristine Graphene

The band structure calculations show that Li and Na adsorptions downshift the Dirac point (ED, DOS(ED) = 0) below the Fermi energy (Figure 5 and Figure 6). Recall that pristine graphene ED coincides with the Fermi energy and appears in the band structure as a K-point, where the Cpz valance (π) and conduction bands (π*) meet in a linear dispersion fashion [8]. Therefore, pristine graphene is a zero-gap semiconductor (i.e., semi-metal) [99]. The ED downshift, due to Li and Na adsorptions on pristine graphene, is indicative of n-doped graphene support via charge transfer from Li and Na. At low Li and Na coverages, the ED is intact (Figure 5d and Figure 6d) and the Li-2s and Na-3s valance orbitals appear in the band structures as flat bands near and above the Fermi energies. Moreover, the DOS spectra below the Fermi energy are unaffected of the adsorption (Figure 5a and Figure 6a), thus supporting a Columbic interaction between the metals and pristine graphene. The DOS and band structure calculations show that Li/graphene and Na/graphene are metallic.
At higher Li coverages, the ED is broken (Figure 5e,f), whereas for the Na case the ED is only broken at the 5/32 ML Na coverage (Figure 6f). However, only for the Li3/graphene case, a small band gap of 0.05 eV is observed in the ED region (Figure 5e,f and Figure 6f), turning Li3/graphene to a semiconductor. However, this effect may be an artifact of the calculations and it will be further examined in the future. All other metal-pristine graphene configurations are metallic. The flat bands below the Fermi energy in Figure 5e and Figure 6e and the wavy bands in Figure 5f and Figure 6f represent Li-2s, Na-3s, and Na-2p valence states. In the DOS spectrum, these bands appear as peaks near and below the Fermi energy (for Na5 case this peak is broadened) and thus, partially occupied. The partial occupancy of these states decreases the Coulombic interaction between the metals and pristine graphene, as the metal coverage increases. This statement also applies for adsorption on defective graphene.

3.3.2. Adsorption on Defective Graphene.

Defects on graphene upshift the ED above the Fermi energy and create band gaps. For example, the double vacancy graphene is a semiconductor [37]. This ED upshift is indicative of a p-doped graphene support [37]. At low Li and Na coverages and for Na3 adsorbed on double vacancy graphene, the ED appears above the Fermi energy (Figure 7d, Figure 8d, and Figure 10e). Liang et al. showed that doping defective graphene with Na ions downshifted ED, but still it remained above the Fermi energy. Therefore, Li and Na adsorptions on defective graphene decrease p-doping. Moreover, at low metal coverages, the observed band gaps further grow as vacancies increase. The DOS spectra (Figure 7a and Figure 8a) show minimal hybridization of the metal valance states with the graphene bands, which supports metal-graphene Coulombic interaction. However, we will show in the next section that the Li-C bond for Li adsorbed on single vacancy graphene has some covalent character.
With the exception of Na3 adsorbed on double vacancy graphene, increased Li and Na coverages on defective graphene push the ED below the Fermi energy, thus changing graphene to n-doped. For Li adsorption on single vacancy graphene the band gap grows due to Li clustering, whereas the opposite is observed for the Na case. Multiple band gaps are observed for select Li and Na clusters adsorbed on double vacancy graphene. Here, all metal-defective graphene configurations are semiconductors. Our Na band structure calculations are in agreement with Liang et al. [37].

3.3.3. Adsorption on GO

Graphene Oxide, at O/C = 43.75%, is an insulator. Its calculated band gap is 5.26 eV, when only epoxides are present in the GO structure (PBE0 calculations without D2 correction) [50]. This value is decreased to 4.40 eV, for GO that includes both epoxides and hydroxyls, under the same DFT calculations. The GO band structure calculations do not reveal ED (Figure 11). Li and Na adsorption on GO is associated with multiple band gaps below the Fermi energy. For example, the largest of these band gaps are 3.23 eV and 3.17 eV for Li/GO and Na/GO, respectively. The largest band gap width and the Li and Na coverages are negatively correlated (Figure 11). The GO DOS spectra below the Fermi energy show minimal contributions from the metals valance states. However, the Li/GO and Na/GO DOS spectra below the Fermi energy are different, which implies structural changes on GO support due to metal adsorption. This statement also applies to adsorptions at higher metal coverages.

3.4. QTAIM analysis

Figure 12 shows the ρ r and 2 ρ r contour plots for Li adsorbed on pristine graphene. The topology of the Laplacian is given by 2 ρ r , since charge is concentrated at 2 ρ r < 0. The 2 ρ r reveal altering shells of charge concentrations and depletions (valence shell charge concentration-VSCC [100]) and (3, -1) saddle points, whereas this information is absent in ρ r (Figure 12a). Aray and Rodriguez [101] stated that maxima are linked with other maxima by trajectories of 2 ρ r originating at (3, -1) saddle points. Evidently, the (3, -1) saddle points that connect Li and C are not shown in Figure 12.

3.4.1. The Metal–C Bond

Table 3 and Table 4 show the QTAIM properties ρ r , 2 ρ r , and (G/ ρ ) r at the Li-C and C-Na bond critical points, respectively (closest metal–C distances). Moreover, for adsorption on GO, the above QTAIM properties are also evaluated at the metal-O bond critical points. Table 3 shows that some Li-C bond critical points for Li3 adsorption on pristine graphene and GO are not revealed by TOPOND. Moreover, for the former case, the Li-C critical point at 2.28 Å appears as a (3, +1) point (i.e., charge minima). In all cases here, with the exception of Li adsorbed on single vacancy graphene, QTAIM shows that metal–C bonds are ionic. Recall that for ionic bonding 2 ρ r c p > 0, (H/ ρ ) r c p > 0 , and (G/ ρ ) r c p 1 , whereas ρ r c p is small. The Li-C bond, for Li adsorption on single vacancy graphene, is analyzed using the bond classification of Espinosa et al. [71]. Here (|V|/ G ) r c p < 1 is indicative of ionic bonding and (|V|/ G ) r c p = 1 is the lower limit of the transit region indicating incipient covalent bond formation. Recall that for this case, the Li-C COOP value is 0.138, and is the highest value observed in this work for metal–C bonds. Therefore, the presence of vacancies in the support may add a covalent character in the metal–C bonds.
QTAIM shows that for adsorption on GO, Li and Na form ionic bonds of similar strength with both C and O atoms. This explains the significantly stronger metal-GO interaction relative to adsorption on pristine graphene, at the same metal coverages. Figure 13 shows a negative correlation between (H/ ρ ) r and ρ r thus, verifying ionic metal–C bonding.

3.4.2. The Metal–metal Bond

Table 5 and Table 6 show the QTAIM properties ρ r , (H/ ρ ) r , (G/ ρ ) r , ξ r , and the sign of 2 ρ r and ELF η r values at the Li-Li and Na-Na bond critical points, respectively. Not all metal–metal bond critical points are revealed for both Li and Na clusters adsorbed on graphene and GO (Table 4 and Table 5). Moreover, the Li-Li bond critical points are not revealed for Li3 adsorbed on defective graphene and for Li3 and Li5 adsorbed on GO. Recall that for metallic bonds ξ > 1 , whereas ξ < 0 , due to ρ < 0 , is indicative of non-metallic non-polar covalent bonds. For Li3 adsorbed on pristine graphene, Na3 adsorbed on pristine and double vacancy graphene, and for Li5 and Na5 adsorbed on defective graphene all intra metal-metal bonds are metallic ( ξ > 1 , Table 3). However, Li5, Na3, Na5 show a mixture of covalent and metallic metal-metal bonds, when are adsorbed on other supports (e.g., Li5 on pristine graphene, Table 5). Therefore, the support selection during metal cluster adsorption could alter the metal-metal bond types.

3.5. ELF Analysis

Ayers and Jenkins stated that metallic bonds rarely have ELF η r values larger than 0.6 [73]. Our η r values, for metallic metal–metal bonds at their bond critical points, are less than 0.43, which is within the Ayers and Jenkins ELF limit for metallic bonds. However, some Li-Li and Na-Na non-metallic covalent bonds have low η r values comparable to the ones for metallic bonds. For example, for Li5 adsorbed on pristine graphene the Li-Li bond with bond length 2.91 Å has a negative metallicity but its η r value is 0.39. A similar case appears for Na3/GO.
Figure 14 shows the ELF isosurfaces for Li5 adsorbed on pristine graphene at η r values 0.5 and 0.94. Recall that ELF analyzes chemical bonding by yielding basins of attractors. At η r = 0.5 the Li-Li valance basin (cyan surface) surrounds several Li atoms, whereas at η r = 0.94 it breaks in smaller isosurfaces. Figure 15 and Figure 16 show the bifurcation diagrams and corresponding isosurfaces for Li and Na adsorbed on pristine and single vacancy graphene, respectively. The values shown in bifurcation diagrams correspond to the η r value, where a separation between the basins occurs. We start our ELF analyzes with Li and Na adsorption on pristine graphene, at low metal coverages.
For Li/graphene, the bifurcation diagram reveals the Li and C core basins (C(Li) and C(C)) and the polysynaptic valence basin V(C,C) (Figure 15a). For Na/graphene, the bifurcation diagram additionally reveals a monosynaptic valance V(Na) (Figure 15g). For Li/graphene, the first bifurcation occurs at η r = 0.04, where C(Li) separates from the disynaptic basin. Here, the second bifurcation occurs at η r = 0.12, for the C(C) separation, whereas the V(C,C) separation appears at η r = 0.69 . Similar η r separation values are found for Na/graphene. Our ELF separation values are in agreement with those found by Matito et al. [81], who studied the bonding in methylalkalimetals (CH3M)n (M = Li, Na, and K; n = 1, 4). The η r C(Li) and C(Na) separation values are indicative of ionic bonding between Li and Na metals and C, in agreement with QTAIM and electronic structure analysis (vide supra). It is noteworthy that our C(Na) η r separation values are never higher than C(Li) values, in agreement with the weaker Na-graphene interaction relative to Li-graphene (Figure 15 and Figure 16). At higher metal coverages, ELF reveals monosynaptic V(Li) and V(Na) (Figure 15b) and polysynaptic V(Li,Li) and V(Na,Na) basins (Figure 15c).
At low metal coverages and for adsorption on single vacancy graphene, the first bifurcation point occurs at higher η r values relative to adsorption on pristine graphene. These higher η r separation values are indicative of stronger metal–C bonding and are in the covalent bond region (Figure 16a). COOP, QTAIM, and ELF agree that the Li-C bond, for adsorption on single vacancy graphene, has some covalent character. However, the Na-C covalent character, for Na adsorbed on single vacancy graphene, indicated by ELF analyses is not agreement with QTAIM, whereas it agrees with COOP. Disparagements between ELF and QTAIM have been reported [81].

4. Conclusions

Li and Na adsorptions on pristine and defective graphene and GO has been studied using periodic DFT structural and electronic calculations paired with COOP, QTAIM, and ELF analyzes. DFT optimized geometries show that increased metal coverage clusters Li and Na during adsorption on graphene, which lowers the efficiency of LIB and NIB. However, Li3 clusters collapse during adsorption on single vacancy graphene and GO. The Eads calculations show that Li and Na adsorption on pristine graphene is unstable at all the metal coverages examined here, whereas the presence of defects in the graphene support stabilizes Li and Na adsorptions. Increased metal coverage weakens metal-graphene and strengthen metal-GO interactions. Adsorption on GO is strong due to simultaneous metal bonding with C and O GO atoms.
DOS calculations show that, at low metal coverages, Li and Na adsorptions on pristine does not affect the DOS spectra below the Fermi energy. This effect is indicative of Coulombic interactions between the adsorbates and the support. These interactions weaken as metal coverages increase, due to the downshift of the metal valances below the Fermi energy. The band calculations show that the ED is pushed below the Fermi energy due to Li and Na adsorptions on pristine graphene and may be associated with small band gaps at higher metal coverages. The presence of defects on graphene is associated with ED shifting and bang gap openings.
The metal–C bonds are analyzed using COOP, QTAIM, and ELF. The COOP metal–C COOP values at the order expected for ionic bonding, for all cases except for adsorption on single vacancy graphene, at low metal coverages. QTAIM shows that only the Li-C bond contains covalent character for adsorption on single vacancy graphene, whereas in all other cases the metal–C bond is ionic. ELF partially agrees with QTAIM and shows that both Li-C and Na-C bonds are covalent for adsorption on single vacancy graphene. QTAIM and EFT are also used to study Li-Li and Na-Na bonds: Li and Na clusters may contain metallic and covalent metal–metal bonds depending on adsorption support (i.e., pristine and graphene and GO).

Author Contributions

Conceptualization, N.D. and M.I.B.; methodology, N.D. and M.I.B.; validation, N.D. and M.I.B.; formal analysis, N.D., I.S., L.G., O.V., and K.R; investigation, N.D.; resources, N.D.; data curation, N.D., I.S., L.G., O.V., and K.R.; writing—original draft preparation, N.D., I.S., L.G, O.V., and K.R; writing—review and editing, N.D., I.S., L.G, and M.I.B; visualization, N.D., O.V., K.R; supervision, N.D.; project administration, N.D.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tarascon, J.-M.; Armand, M. Issues and challenges facing rechargeable lithium batteries. Nature 2001, 414, 359–367. [Google Scholar] [CrossRef] [PubMed]
  2. Goodenough, J.B.; Park, K.-S. The Li-Ion Rechargeable Battery: A Perspective. J. Am. Chem. Soc. 2013, 135, 1167–1176. [Google Scholar] [CrossRef] [PubMed]
  3. Bruce, P.G.; Scrosati, B.; Tarascon, J.-M. Nanomaterials for Rechargeable Lithium Batteries. Angew. Chem. Int. Ed. 2008, 47, 2930–2946. [Google Scholar] [CrossRef] [PubMed]
  4. Dahn, J.R.; Zheng, T.; Liu, Y.H.; Xue, J.S. Mechanisms for Lithium Insertion in Carbonaceous Materials. Science 1995, 270, 590–593. [Google Scholar] [CrossRef]
  5. Kheirabadi, N.; Shafiekhani, A. Graphene/Li-ion battery. J. Appl. Phys. 2012, 112, 124323–124325. [Google Scholar] [CrossRef]
  6. Yang, S.; Li, S.; Tang, S.; Dong, W.; Sun, W.; Shen, D.; Wang, M. Sodium adsorption and intercalation in bilayer graphene from density functional theory calculations. Theor. Chem. Acc. 2016, 135, 164. [Google Scholar] [CrossRef]
  7. Winter, M.; Besenhard, J.O.; Spahr, M.E.; Novak, P. Insertion Electrode Materials for Rechargeable Lithium Batteries. Adv. Mater. 1998, 10, 725–763. [Google Scholar] [CrossRef]
  8. Geim, A.K.; Novoselov, K.S. The rise of graphene. Nat. Mater. 2007, 6, 183–191. [Google Scholar] [CrossRef]
  9. Novoselov, K.S.; Geim, A.K.; Morozov, S.V.; Jiang, D.; Katsnelson, M.I.; Grigorieva, I.V.; Firsov, A.A. Electric field effect in anatomically thin carbon films. Science 2004, 306, 666–669. [Google Scholar] [CrossRef]
  10. Yoo, E.; Kim, J.; Hosono, E.; Zhou, H.-S.; Kudo, T.; Honma, I. Large Reversible Li Storage of Graphene Nanosheet Families for Use in Rechargeable Lithium Ion Batteries. Nano Lett. 2008, 8, 2277–2282. [Google Scholar] [CrossRef]
  11. Liu, Y.; Xie, B.; Zhang, Z.; Zheng, Q.; Xu, Z. Mechanical properties of graphene papers. J. Mech. Phys. Solids 2012, 60, 591–605. [Google Scholar] [CrossRef] [Green Version]
  12. Pop, E.; Varshney, V.; Roy, A. Thermal properties of graphene: Fundamentals and applications. MRS Bull. 2012, 37, 1273–1281. [Google Scholar] [CrossRef] [Green Version]
  13. Castro Neto, A.H.; Guinea, F.; Peres, N.M.R.; Novoselov, K.S.; Geim, A.K. The electronic properties of graphene. Rev. Mod. Phys. 2009, 81, 110–155. [Google Scholar] [CrossRef]
  14. Wu, M.; Cao, C.; Jiang, J.Z. Light non-metallic atom (B, N, O and F)-doped graphene: A first-principles study. Nanotechnology 2010, 21, 505202–505206. [Google Scholar] [CrossRef] [PubMed]
  15. Pollak, E.; Geng, B.; Jeon, K.-J.; Lucas, I.T.; Richardson, T.J.; Wang, F.; Kostecki, R. The Interaction of Li+ with Single-Layer and Few-Layer Graphene. Nano Lett. 2010, 10, 3386–3388. [Google Scholar] [CrossRef] [PubMed]
  16. Chan, K.T.; Neaton, J.B.; Cohen, M.L. First-Principles Study of Metal Adatom Adsorption on Graphene. Phys. Rev. B 2008, 77, 235430. [Google Scholar] [CrossRef]
  17. Dimakis, N.; Valdez, D.; Flor, F.A.; Salgado, A.; Adjibi, K.; Vargas, S.; Saenz, J. Density functional theory calculations on alkali and the alkaline Ca atoms adsorbed on graphene monolayers. Appl. Surf. Sci. 2017, 413, 197–208. [Google Scholar] [CrossRef]
  18. Nakada, K.; Ishii, A. DFT Calculation for adatom adsorption on Graphene. In Graphene Simulation; Gong, J.R., Ed.; InTech Inc.: Zagreb, Croatia, 2011; pp. 3–20. [Google Scholar]
  19. Medeiros, P.V.C.; de Brito Mota, F.; Mascarenhas, A.J.S.; de Castilho, C.M.C. Adsorption of Monovalent Metal Atoms on Graphene: A Theoretical Approach. Nanotechnology 2010, 21, 115701–115706. [Google Scholar] [CrossRef]
  20. Lee, E.; Persson, K.A. Li Absorption and Intercalation in Single Layer Graphene and Few Layer Graphene by First Principles. Nano Lett. 2012, 12, 4624–4628. [Google Scholar] [CrossRef]
  21. Okamoto, Y. Density Functional Theory Calculations of Lithium Adsorption and Insertion to Defect-Free and Defective Graphene. J. Phys. Chem. C 2016, 120, 14009–14014. [Google Scholar] [CrossRef]
  22. Yildirim, H.; Kinaci, A.; Zhao, Z.-J.; Chan, M.K.Y.; Greeley, J.P. First-Principles Analysis of Defect-Mediated Li Adsorption on Graphene. ACS Appl. Mater. Interfaces 2014, 6, 21141–21150. [Google Scholar] [CrossRef] [PubMed]
  23. Shiota, K.; Kawai, T. Li atom adsorption on graphene with various defects for large-capacity Li ion batteries: First-principles calculations. Jpn. J. Appl. Phys. 2017, 56, 06GE11. [Google Scholar] [CrossRef] [Green Version]
  24. Datta, D.; Li, J.; Koratkar, N.; Shenoy, V.B. Enhanced lithiation in defective graphene. Carbon 2014, 80, 305–310. [Google Scholar] [CrossRef]
  25. Yu, Y.-X. Can all nitrogen-doped defects improve the performance of graphene anode materials for lithium-ion batteries? Phys. Chem. Chem. Phys. 2013, 15, 16819–16827. [Google Scholar] [CrossRef] [PubMed]
  26. Slater, M.D.; Kim, D.; Lee, E.; Johnson, C.S. Sodium-ion batteries. Adv. Funct. Mater. 2013, 23, 947–958. [Google Scholar] [CrossRef]
  27. Teng, F.Z.; McDonough, W.F.; Rudnick, R.L.; Dalpé, C.; Tomascak, P.B.; Chappell, B.W.; Gao, S. Lithium isotopic composition and concentration of the upper continental crust. Geochim. Cosmochim. Acta 2004, 68, 4167–4178. [Google Scholar] [CrossRef] [Green Version]
  28. Palomares, V.; Serras, P.; Villaluenga, I.; Hueso, K.B.; Carretero-González, J.; Rojo, T. Na-ion batteries, recent advances and present challenges to become low cost energy storage systems. Energy Environ. Sci. 2012, 5, 5884–5901. [Google Scholar] [CrossRef]
  29. Harris, S.J.; Timmons, A.; Baker, D.R.; Monroe, C. Direct in situ Measurements of Li Transport in Li-Ion Battery Negative Electrodes. Chem. Phys. Lett. 2010, 485, 265–274. [Google Scholar] [CrossRef]
  30. Zhou, J.; Sun, Q.; Wang, Q.; Jena, P. Tailoring Li adsorption on graphene. Phys. Rev. B 2014, 90, 205427. [Google Scholar] [CrossRef]
  31. Liu, M.F.; Kutana, A.; Liu, Y.; Yakobson, B.I. First-Principles Studies of Li Nucleation on Graphene. J. Phys. Chem. Lett. 2014, 5, 1225–1229. [Google Scholar] [CrossRef]
  32. Gómez-Navarro, C.; Meyer, J.C.; Sundaram, R.S.; Chuvilin, A.; Kurasch, S.; Burghard, M.; Kern, K.; Kaiser, U.A. Atomic Structure of Reduced Graphene Oxide. Nano Lett. 2010, 10, 1144–1148. [Google Scholar] [CrossRef] [Green Version]
  33. Acik, M.; Lee, G.; Mattevi, C.; Pirkle, A.; Wallace, R.M.; Chhowalla, M.; Cho, K.; Chabal, Y. The Role of Oxygen during Thermal Reduction of Graphene Oxide Studied by Infrared Absorption Spectroscopy. J. Phys. Chem. C 2011, 115, 19761–19781. [Google Scholar] [CrossRef]
  34. Nobuhara, K.; Nakayama, H.; Nose, M.; Nakanishi, S.; Iba, H. First-Principles Study of Alkali Metal-Graphite Intercalation Compounds. J. Power Sources 2013, 243, 585–587. [Google Scholar] [CrossRef]
  35. Asher, R.C. A lamellar compound of sodium and graphite. J. Inorg. Nucl. Chem. 1959, 10, 238–249. [Google Scholar] [CrossRef]
  36. Datta, D.; Li, J.; Shenoy, V.B. Defective Graphene as a High-Capacity Anode Material for Na- and Ca-Ion Batteries. ACS Appl. Mater. Interfaces 2014, 6, 1788–1795. [Google Scholar] [CrossRef]
  37. Liang, Z.; Fan, X.; Zheng, W.; Singh, D.J. Adsorption and Formation of Small Na Clusters on Pristine and Double-Vacancy Graphene for Anodes of Na-Ion Batteries. ACS Appl. Mater. Interfaces 2017, 9, 17076–17084. [Google Scholar] [CrossRef] [PubMed]
  38. Bader, R.F.W. Atoms in Molecules. In Encyclopedia of Computational Chemistry; John Wiley and Sons: London, UK, 1998; Volume 1, pp. 64–86. [Google Scholar]
  39. Koritsanszky, T.S.; Coppens, P. Chemical Applications of X-ray Charge-Density Analysis. Chem. Rev. 2001, 101, 1583–1628. [Google Scholar] [CrossRef]
  40. Bader, R.F.W. Atoms in Molecules; Wiley Online Library: Hoboken, NJ, USA, 1990. [Google Scholar]
  41. Becke, A.D.; Edgecombe, K.E. A simple measure of electron localization in atomic and molecular systems. J. Chem. Phys. 1990, 92, 5397–5403. [Google Scholar] [CrossRef]
  42. Jabłoński, M.; Palusiak, M. Basis set and method dependence in atoms in molecules calculations. J. Phys. Chem. A 2010, 114, 2240–2244. [Google Scholar] [CrossRef]
  43. Fuentealba, P.; Chamorro, E.; Santos, J.C. Understanding and using the electron localization function. In Theoretical Aspects of Chemical Reactivity; Toro-Labbé, A., Ed.; Theoretical and Computational Chemistry Series; Elsevier: Amsterdam, The Netherlands, 2007; Volume 19, pp. 57–85. [Google Scholar]
  44. Gupta, S.; Carrizosa, S.B.; Jasisnki, J.; Dimakis, N. Charge Transfer Dynamical Processes at Graphene-Transition Metal Oxides/Electrolyte Interface for Energy Storage: Insights from In-Situ Raman Spectroelectrochemistry. AIP Adv.J. Raman Spectrosc. 2018, 8, 065225. [Google Scholar] [CrossRef]
  45. Mallia, G.; Orlando, R.; Roetti, C.; Ugliengo, P.; Dovesi, R. F center in LiF: A quantum mechanical ab initio investigation of the hyperfine interaction between the unpaired electron at the vacancy and its first seven neighbors. Phys. Rev. B 2001, 63, 235102–235107. [Google Scholar] [CrossRef]
  46. Zhou, S.; Bongiorno, A. Density Functional Theory Modeling of Multilayer “Epitaxial” Graphene Oxide. Acc. Chem. Res. 2014, 47, 3331–3339. [Google Scholar] [CrossRef] [PubMed]
  47. Dovesi, R.; Saunders, V.R.; Roetti, C.; Zicovich-Wilson, C.M.; Pascale, F.; Civalleri, B.; Doll, K.; Harrison, N.M.; Bush, I.J.; D’Arco, P.; et al. CRYSTAL17 User’s Manual; University Torino: Torino, Italy, 2017. [Google Scholar]
  48. Shi, H.; Eglitis, R.; Borstel, G. Ab initio calculations of the BaF2 bulk and surface F centres. J. Phys. Condens. Matter 2006, 18, 8367. [Google Scholar] [CrossRef]
  49. Dimakis, N.; Flor, F.A.; Salgado, A.; Adjibi, K.; Vargas, S.; Saenz, J. Density functional theory calculations on transition metal atoms adsorbed on graphene monolayers. Appl. Surf. Sci. 2017, 421, 252–259. [Google Scholar] [CrossRef]
  50. Gupta, S.; Dimakis, N. Computational Predictions of Electronic Properties of Graphene with Defects, Adsorbed Transition Metal-Oxides and Water using Density Functional Theory. Appl. Surf. Sci. 2019, 467, 760–772. [Google Scholar] [CrossRef]
  51. Gupta, S.; Meek, R.; Evans, B.; Dimakis, N. Graphene-based “hybrid” aerogels with carbon nanotubes: Mesoporous network–functionality promoted defect density and electrochemical activity correlations. J. Appl. Phys. 2018, 124, 124304–124317. [Google Scholar] [CrossRef]
  52. Ernzerhof, M.; Scuseria, G.E. Assessment of the Perdew–Burke–Ernzerhof exchange-correlation functional. J. Chem. Phys. 1999, 110, 5029–5036. [Google Scholar] [CrossRef]
  53. Adamo, C.; Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model. J. Chem. Phys. 1999, 110, 6158–6170. [Google Scholar] [CrossRef]
  54. Grimme, S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787–1799. [Google Scholar] [CrossRef]
  55. Peintinger, M.F.; Oliveira, D.V.; Bredow, T. Consistent Gaussian basis sets of triple-zeta valence with polarization quality for solid-state calculations. J. Comput. Chem. 2013, 34, 451–459. [Google Scholar] [CrossRef] [PubMed]
  56. Mahmoud, A.; Erba, A.; El-Kelany, K.E.; Rérat, M.; Orlando, R. Low-temperature phase of BaTiO3: Piezoelectric, dielectric, elastic, and photoelastic properties from ab initio simulations. Phys. Rev. B 2014, 89, 045103–045109. [Google Scholar] [CrossRef]
  57. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188. [Google Scholar] [CrossRef]
  58. Gilat, G.; Raubenheimer, L. Accurate numerical method for calculating frequency-distribution functions in solids. Phys. Rev. 1966, 144, 390. [Google Scholar] [CrossRef]
  59. Gilat, G. Analysis of methods for calculating spectral properties in solids. J. Comput. Phys. 1972, 10, 432–465. [Google Scholar] [CrossRef]
  60. Anderson, D.G. Iterative procedures for nonlinear integral equations. J. ACM 1965, 12, 547–560. [Google Scholar] [CrossRef]
  61. Hoffmann, R. A chemical and theoretical way to look at bonding on surfaces. Rev. Mod. Phys. 1988, 60, 601–628. [Google Scholar] [CrossRef]
  62. Mulliken, R.S. Electronic population analysis on LCAO-MO [linear combination of atomic orbital-molecular orbital] molecular wave functions. I. J. Chem. Phys. 1955, 23, 1833–1840. [Google Scholar] [CrossRef]
  63. Bultinck, P.; Van Alsenoy, C.; Ayers, P.W.; Carbó-Dorca, R. Critical analysis and extension of the Hirshfeld atoms in molecules. J. Chem. Phys. 2007, 126, 144111–144119. [Google Scholar] [CrossRef]
  64. Marenich, A.V.; Jerome, S.V.; Cramer, C.J.; Truhlar, D.G. Charge Model 5: An Extension of Hirshfeld Population Analysis for the Accurate Description of Molecular Interactions in Gaseous and Condensed Phases. J. Chem. Theory Comput. 2012, 8, 527–541. [Google Scholar] [CrossRef]
  65. Davidson, E.R.; Chakravorty, S. A test of the Hirshfeld definition of atomic charges and moments. Theor. Chim. Acta 1992, 83, 319–330. [Google Scholar] [CrossRef]
  66. Boys, S.F.; Bernardi, F. The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Mol. Phys. 1970, 19, 553–566. [Google Scholar] [CrossRef]
  67. Dimakis, N.; Flor, F.A.; Navarro, N.E.; Salgado, A.; Smotkin, E.S. Adsorption of Carbon Monoxide on Platinum–Ruthenium, Platinum–Osmium, Platinum–Ruthenium–Osmium, and Platinum–Ruthenium–Osmium–Iridium Alloys. J. Phys. Chem. C 2016, 120, 10427–10441. [Google Scholar] [CrossRef]
  68. Gatti, C. TOPOND-96: An electron density topological program for systems periodic in N (N=0-3) dimensions, User’s Manual. In User’s Manual; CNR-CSRSRC: Milano, Italy, 1996. [Google Scholar]
  69. Macchi, P.; Proserpio, D.M.; Sironi, A. Experimental electron density in a transition metal dimer: Metal–metal and metal–ligand bonds. J. Am. Chem. Soc. 1998, 120, 13429–13435. [Google Scholar] [CrossRef]
  70. Gatti, C. Chemical bonding in crystals: New directions. Z. Kristallogr. Cryst. Mater. 2005, 220, 399–457. [Google Scholar] [CrossRef]
  71. Espinosa, E.; Alkorta, I.; Elguero, J.; Molins, E. From weak to strong interactions: A comprehensive analysis of the topological and energetic properties of the electron density distribution involving X–H⋯F–Y systems. J. Chem. Phys. 2002, 117, 5529–5542. [Google Scholar] [CrossRef]
  72. Azizi, A.; Momen, R.; Xu, T.; Kirk, S.R.; Jenkins, S. Non-nuclear attractors in small charged lithium clusters, Limq (m = 2–5, q = ±1), with QTAIM and the Ehrenfest force partitioning. Phys. Chem. Chem. Phys. 2018, 20, 24695–24707. [Google Scholar] [CrossRef]
  73. Ayers, P.W.; Jenkins, S. Bond metallicity measures. Comput. Theor. Chem. 2015, 1053, 112–122. [Google Scholar] [CrossRef]
  74. Jenkins, S.; Ayers, P.W.; Kirk, S.R.; Mori-Sánchez, P.; Martín Pendás, A. Bond metallicity of materials from real space charge density distributions. Chem. Phys. Lett. 2009, 471, 174–177. [Google Scholar] [CrossRef]
  75. Centeno, S.A.; Vila, A.; Barro, L. Characterization of unprocessed historic platinum photographic papers by Raman, FTIR and XRF. Microchem. J. 2014, 114, 8–15. [Google Scholar] [CrossRef]
  76. Contreras-García, J.; Pendás, Á.M.; Silvi, B.; Recio, J.M. Useful applications of the electron localization function in high pressure crystal chemistry. J. Phys. Chem. Solids 2008, 69, 2204–2207. [Google Scholar] [CrossRef]
  77. Silvi, B.; Savin, A. Classification of chemical bonds based on topological analysis of electron localization functions. Nature 1994, 371, 683–686. [Google Scholar] [CrossRef]
  78. Simas, A.M.; Sagar, R.P.; Ku, A.C.T.; Smith, V.H.J. The radial charge distribution and the shell structure of atoms and ions. Can. J. Chem. 1988, 66, 1923–1930. [Google Scholar] [CrossRef] [Green Version]
  79. Kohout, M.; Savin, A. Atomic shell structure and electron numbers. Int. J. Quantum Chem. 1996, 60, 875–882. [Google Scholar] [CrossRef] [Green Version]
  80. Savin, A.; Nesper, R.; Wengert, S.; Fassler, T. ELF: The Electron Localization Function. Angew. Chem. Int. Ed. Engl. 1997, 36, 1808–1832. [Google Scholar] [CrossRef]
  81. Matito, E.; Poater, J.; Bickelhaupt, F.M.; Solà, M. Bonding in Methylalkalimetals (CH3M)n (M = Li, Na, K; n = 1, 4). Agreement and Divergences between AIM and ELF Analyses. J. Phys. Chem. B 2006, 110, 7189–7198. [Google Scholar] [CrossRef]
  82. Savin, A.; Silvi, B.; Colonna, F. Topological analysis of the electron localization function applied to delocalized bonds. Can. J. Chem. 1996, 74, 1088–1096. [Google Scholar] [CrossRef] [Green Version]
  83. Marx, D.; Savin, A. Topological Bifurcation Analysis: Electronic Structure of CH5+. Angew. Chem. Int. Ed. Engl. 1997, 36, 2077–2080. [Google Scholar] [CrossRef]
  84. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera—A visualization system for exploratory research and analysis. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef]
  85. Nadler, M.R.; Kempier, C.P. Crystallographic Data 186. Lithium. Anal. Chem. 1959, 31, 2109. [Google Scholar] [CrossRef]
  86. Aruja, E.; Perlitz, H. Neubestimmung der Gitterkonstante von Natrium. Z. Kristallogr. Cryst. Mater. 2015, 100, 195–200. [Google Scholar] [CrossRef]
  87. Kittel, C. Introduction to Solid State Physics, 8th ed.; Wiley: New York, NY, USA, 2005. [Google Scholar]
  88. Maezono, R.; Towler, M.D.; Lee, Y.; Needs, R.J. Quantum Monte Carlo Study of Sodium. Phys. Rev. B 2003, 68, 165103–165109. [Google Scholar] [CrossRef]
  89. Asher, R.C.; Wilson, S.A. Lamellar Compound of Sodium with Graphite. Nature 1958, 181, 409–410. [Google Scholar] [CrossRef]
  90. Pašti, I.A.; Jovanović, A.; Dobrota, A.S.; Mentus, S.V.; Johansson, B.J.; Skorodumova, N.V. Atomic adsorption on pristine graphene along the Periodic Table of Elements—From PBE to non-local functionals. Appl. Surf. Sci. 2018, 436, 433–440. [Google Scholar] [CrossRef]
  91. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [PubMed]
  92. Nakada, K.; Ishii, A. Migration of adatom adsorption on graphene using DFT calculation. Solid State Commun. 2011, 151, 13–16. [Google Scholar] [CrossRef]
  93. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104–154118. [Google Scholar] [CrossRef] [PubMed]
  94. Pašti, I.A.; Jovanović, A.; Dobrota, A.S.; Mentus, S.V.; Johansson, B.J.; Skorodumova, N.V. Atomic adsorption on graphene with a single vacancy: Systematic DFT study through the periodic table of elements. Phys. Chem. Chem. Phys. 2018, 20, 858–865. [Google Scholar] [CrossRef]
  95. Lee, K.; Murray, E.D.; Kong, L.; Lundqvist, B.I.; Langreth, D.C. Higher-accuracy van der Waals density functional. Phys. Rev. B 2010, 82, 081101–081104. [Google Scholar] [CrossRef]
  96. Dobrota, A.S.; Pašti, I.A.; Mentus, S.V.; Johansson, B.J.; Skorodumova, N.V. Functionalized graphene for sodium battery applications: The DFT insights. Electrochim. Acta 2017, 250, 185–195. [Google Scholar] [CrossRef]
  97. Doan, J.; Navarro, N.E.; Kumari, D.; Anderson, K.; Kingston, E.; Johnson, C.; Vong, A.; Dimakis, N.; Smotkin, E.S. Symmetry-based IR group modes as dynamic probes of Nafion ion exchange site structure. Polymer 2015, 73, 34–41. [Google Scholar] [CrossRef]
  98. Grechnev, A.; Ahuja, R.; Eriksson, O. Balanced crystal orbital overlap population—A tool for analysing chemical bonds in solids. J. Phys. Condens. Matter 2003, 15, 7751. [Google Scholar] [CrossRef]
  99. Rani, P.; Jindal, V.K. Designing band gap of graphene by B and N dopant atoms. RSC Adv. 2013, 3, 802–812. [Google Scholar] [CrossRef] [Green Version]
  100. Bader, R.F.W.; MacDougall, P.J.; Lau, C.D.H. Bonded and nonbonded charge concentrations and their relation to molecular geometry and reactivity. J. Am. Chem. Soc. 1984, 106, 1594–1605. [Google Scholar] [CrossRef]
  101. Aray, Y.; Rodriguez, J. Study of CO adsorption on the Fe (100) surface using the Laplacian of the electronic charge density. Surf. Sci. 1998, 405, L532–L541. [Google Scholar] [CrossRef]
Sample Availability: Input CRYSTAL17 files are available upon request.
Figure 1. Optimized geometries, heights h (blue arrows) for the metal (Li, Na) closest to graphene plane, the closest metal–metal distances (red arrows), and Ec for Lin (yellow spheres) (ac) and Nan (purple spheres) (df) (n = 1, 3, 5) adsorbed on pristine graphene (C, gray spheres). Values in parentheses include D2 correction. The thin black line is the unit cell boundary.
Figure 1. Optimized geometries, heights h (blue arrows) for the metal (Li, Na) closest to graphene plane, the closest metal–metal distances (red arrows), and Ec for Lin (yellow spheres) (ac) and Nan (purple spheres) (df) (n = 1, 3, 5) adsorbed on pristine graphene (C, gray spheres). Values in parentheses include D2 correction. The thin black line is the unit cell boundary.
Molecules 24 00754 g001
Figure 2. Optimized geometries, heights h (blue arrows) for the metal (Li, Na) closest to graphene plane, the closest metal–metal distances (red arrows), and Ec for Lin (yellow spheres) (af) and Nan (purple spheres) (gn) (n = 1, 3, 5) adsorbed on defective graphene (C, gray spheres, “ghost” C, blue spheres) with one and two C vacancies. All calculations include D2 correction. The thin black line is the unit cell boundary. The thick dashed black line is the carbon plane.
Figure 2. Optimized geometries, heights h (blue arrows) for the metal (Li, Na) closest to graphene plane, the closest metal–metal distances (red arrows), and Ec for Lin (yellow spheres) (af) and Nan (purple spheres) (gn) (n = 1, 3, 5) adsorbed on defective graphene (C, gray spheres, “ghost” C, blue spheres) with one and two C vacancies. All calculations include D2 correction. The thin black line is the unit cell boundary. The thick dashed black line is the carbon plane.
Molecules 24 00754 g002
Figure 3. Optimized geometries, heights h (blue arrows) for the metal (Li, Na) closest to the C plane, the closest metal–metal distances (red arrows), and Ec for (a) Lin (yellow spheres) and (b) Nan (purple spheres) (ce) (n = 1, 3) adsorbed on GO (C, gray spheres; O, red spheres; H white spheres). Values in parentheses include D2 correction. The thin black line is the unit cell boundary. The thick dashed black line is the carbon plane.
Figure 3. Optimized geometries, heights h (blue arrows) for the metal (Li, Na) closest to the C plane, the closest metal–metal distances (red arrows), and Ec for (a) Lin (yellow spheres) and (b) Nan (purple spheres) (ce) (n = 1, 3) adsorbed on GO (C, gray spheres; O, red spheres; H white spheres). Values in parentheses include D2 correction. The thin black line is the unit cell boundary. The thick dashed black line is the carbon plane.
Molecules 24 00754 g003
Figure 4. The Eads for (a) Li and (b) Na adsorptions on various supports. The horizontal dashed lines are the corresponding bulk metals Eads. Calculations include D2 correction.
Figure 4. The Eads for (a) Li and (b) Na adsorptions on various supports. The horizontal dashed lines are the corresponding bulk metals Eads. Calculations include D2 correction.
Molecules 24 00754 g004
Figure 5. Li adsorption on pristine graphene (ac) DOS spectra and (df) band structures. The red horizontal line is the Fermi energy.
Figure 5. Li adsorption on pristine graphene (ac) DOS spectra and (df) band structures. The red horizontal line is the Fermi energy.
Molecules 24 00754 g005
Figure 6. Na adsorption on pristine graphene (ac) DOS spectra and (df) band structures. The red horizontal line is the Fermi energy.
Figure 6. Na adsorption on pristine graphene (ac) DOS spectra and (df) band structures. The red horizontal line is the Fermi energy.
Molecules 24 00754 g006
Figure 7. Li adsorption on defective graphene with a single vacancy per 32 C (ac) DOS spectra and (df) band structures. The red horizontal line is the Fermi energy. Blue arrows show band gaps.
Figure 7. Li adsorption on defective graphene with a single vacancy per 32 C (ac) DOS spectra and (df) band structures. The red horizontal line is the Fermi energy. Blue arrows show band gaps.
Molecules 24 00754 g007
Figure 8. Li adsorption on defective graphene with double vacancies per 32 C (ac) DOS spectra and (df) band structures. The red horizontal line is the Fermi energy. Blue arrows show band gaps.
Figure 8. Li adsorption on defective graphene with double vacancies per 32 C (ac) DOS spectra and (df) band structures. The red horizontal line is the Fermi energy. Blue arrows show band gaps.
Molecules 24 00754 g008
Figure 9. Na adsorption on defective graphene with single vacancy per 32 C. (ac) DOS spectra and (df) band structures. The red horizontal line is the Fermi energy. Blue arrows show band gaps.
Figure 9. Na adsorption on defective graphene with single vacancy per 32 C. (ac) DOS spectra and (df) band structures. The red horizontal line is the Fermi energy. Blue arrows show band gaps.
Molecules 24 00754 g009
Figure 10. Na adsorption on defective graphene with double vacancies per 32 C (ac) DOS spectra and (df) band structures. The red horizontal line is the Fermi energy. Blue arrows show band gaps.
Figure 10. Na adsorption on defective graphene with double vacancies per 32 C (ac) DOS spectra and (df) band structures. The red horizontal line is the Fermi energy. Blue arrows show band gaps.
Molecules 24 00754 g010
Figure 11. DOS spectra and band structures for adsorption on GO (ad) Li and (eh) Na. The red horizontal line is the Fermi energy. Blue arrows show band gaps close to the Fermi energy.
Figure 11. DOS spectra and band structures for adsorption on GO (ad) Li and (eh) Na. The red horizontal line is the Fermi energy. Blue arrows show band gaps close to the Fermi energy.
Molecules 24 00754 g011
Figure 12. (a) Contour plots for Li/graphene electron density ρ r . (b) The negative of its Laplacian 2 ρ r , which reveals saddle points (3, -1). Colors as follows: Red: High ρ r ; Blue: Low ρ r .
Figure 12. (a) Contour plots for Li/graphene electron density ρ r . (b) The negative of its Laplacian 2 ρ r , which reveals saddle points (3, -1). Colors as follows: Red: High ρ r ; Blue: Low ρ r .
Molecules 24 00754 g012
Figure 13. QTAIM parameters vs. metal–C distance for Li and Na adsorbed on graphene and GO. (a) ρ r , (b) 2 ρ r , and (c) (H/ ρ ) r at metal–C bond critical points (cp).
Figure 13. QTAIM parameters vs. metal–C distance for Li and Na adsorbed on graphene and GO. (a) ρ r , (b) 2 ρ r , and (c) (H/ ρ ) r at metal–C bond critical points (cp).
Molecules 24 00754 g013
Figure 14. Electron localization isosurfaces for Li5 adsorbed on pristine graphene at ELF values. (a) 0.5 and (b) 0.94. Colors are as follows: C(Li) blue, C(C) red, V(Ci,Cj) green, and V(Li) and V(Li,Li) cyan. Calculations include D2 correction.
Figure 14. Electron localization isosurfaces for Li5 adsorbed on pristine graphene at ELF values. (a) 0.5 and (b) 0.94. Colors are as follows: C(Li) blue, C(C) red, V(Ci,Cj) green, and V(Li) and V(Li,Li) cyan. Calculations include D2 correction.
Molecules 24 00754 g014
Figure 15. (ac) Bifurcation diagrams for Lin and Nan n = 1, 3, 5 adsorbed on pristine graphene. Colors are as follows: blue and red refer to Li and Na adsorption, respectively and black refers to both metals. (di) corresponding ELF isosurfaces at ELF = 0.75. Colors are as follows: C(M) blue, C(C) red, V(Ci,Cj) green, and V(M) and V(M,M) cyan, M = Li, Na. Calculations include D2 correction.
Figure 15. (ac) Bifurcation diagrams for Lin and Nan n = 1, 3, 5 adsorbed on pristine graphene. Colors are as follows: blue and red refer to Li and Na adsorption, respectively and black refers to both metals. (di) corresponding ELF isosurfaces at ELF = 0.75. Colors are as follows: C(M) blue, C(C) red, V(Ci,Cj) green, and V(M) and V(M,M) cyan, M = Li, Na. Calculations include D2 correction.
Molecules 24 00754 g015
Figure 16. (ac) Bifurcation diagrams for Lin and Nan n = 1, 3, 5 adsorbed on single vacancy defective graphene. Colors are as follows: blue and red refer to Li and Na adsorption, respectively, and black refers to both metals. (di) Corresponding ELF isosurfaces at ELF = 0.75. Colors are as follows: C(M) blue, C(C) red, V(Ci,Cj) green, and V(M) and V(M,M) cyan, M = Li, Na. M* is the metal atom close to the defect. Calculations include D2 correction.
Figure 16. (ac) Bifurcation diagrams for Lin and Nan n = 1, 3, 5 adsorbed on single vacancy defective graphene. Colors are as follows: blue and red refer to Li and Na adsorption, respectively, and black refers to both metals. (di) Corresponding ELF isosurfaces at ELF = 0.75. Colors are as follows: C(M) blue, C(C) red, V(Ci,Cj) green, and V(M) and V(M,M) cyan, M = Li, Na. M* is the metal atom close to the defect. Calculations include D2 correction.
Molecules 24 00754 g016
Table 1. Calculated Eads for Li and Na adsorbed on graphene compared with past work. Values in parentheses include D2 correction. Unless otherwise stated, the graphene support contains 32 atoms per unit cell (4 × 4 supercell).
Table 1. Calculated Eads for Li and Na adsorbed on graphene compared with past work. Values in parentheses include D2 correction. Unless otherwise stated, the graphene support contains 32 atoms per unit cell (4 × 4 supercell).
SystemVacanciesCurrent WorkPast work
Eads (eV)Method
Li /Graphene0−1.34 (−1.33)−1.19 [90], −1.096 [16]PBE [91]
−1.36 [92]LDA
−1.56 [90]PBE+D2 [54,91]
−1.23 [90]PBE+D3 [91,93]
−1.05 [90]PBE+vdW-DF2 [91]
1(−3.44)−2.65 [94]PBE [91]
−2.94 [94]PBE+D2 [54,91]
−2.72 [94]PBE+D3 [91,93]
−2.59 [94]PBE+vdW-DF2 [91,95]
Na /Graphene0−0.59 (−0.59)−0.55 [90], −0.462 [16]PBE [91]
−0.72 [92] aLDA
−0.93 [90]PBE+D2 [54,91]
−0.64 [90]PBE+D3 [91,93]
−0.49 [90]PBE+vdW-DF2 [91,95]
−0.88 [37]PBE+D2 [91]
1(−2.35)−1.88 [94]PBE [91]
−2.20 [94]PBE+D2 [54,91]
−2.00 [94]PBE+D3 [91,93]
−1.74 [94]PBE+vdW-DF2 [91,95]
a 3 × 3 supercell.
Table 2. Li and Na Hirshfeld and Mulliken electronic charges (average values per metal) and metal–C COOP for the smaller metal–C distances, for Li and Na absorbed on graphene and GO. Values in parentheses include D2 correction.
Table 2. Li and Na Hirshfeld and Mulliken electronic charges (average values per metal) and metal–C COOP for the smaller metal–C distances, for Li and Na absorbed on graphene and GO. Values in parentheses include D2 correction.
MetalsSupportVacanciesMetal charge (e)Metal–C COOP a
HirshfeldMulliken
LiGraphene02.27 (2.27)2.66 (2.64)0.041 (0.040)
1(2.37)(2.73)(0.138) 3
2(2.31)(2.67)(0.039) 6
GO02.13 (2.14)2.43 (2.45)0.085 (0.082) 1
Li3Graphene02.98 (2.98)2.97 (2.97)0.024 (0.024)
1(2.64)(2.84)(0.036) 6
2(2.76)(2.84)(0.006) 4
GO02.12 (2.13)2.42 (2.43)0.088 (0.090) 1
Li5Graphene03.01 (3.02)2.99 (2.99)0.015 (0.022)
1(2.81)(2.91)(0.029) 6
2(2.84)(2.43)(0.036) 6
NaGraphene010.53 (10.53)10.81 (10.58)0.043 (0.042)
1(10.32)(10.76)(0.116) 3
2(10.21)(10.68)(0.043) 6
GO010.13 (10.13)10.53 (10.54)0.051(0.050) 1
Na3Graphene010.85 (10.85)10.97 (10.98)0.039 (0.039)
1(10.72)(10.97)(0.039) 6
2(10.85)(10.97)(0.035) 6
GO010.72 (10.52)10.84 (10.64)0.018 (0.032)
Na5Graphene010.95 (10.95)11.03 (11.04)0.031 (0.037) 6
1(10.85)(11.01)(0.039) 6
2(10.86)(11.00)(0.035) 6
a Superscript: Number of metal–C COOP values included in the calculations, i = 1−6.
Table 3. QTAIM parameters ρ r , 2 ρ r , (H/ ρ ) r , and (G/ ρ ) r at the Li-C bond critical points for the shortest Li-C distances, for Li adsorbed on graphene, with and without vacancies, and GO. For adsorption on GO, these properties are also calculated at the Li-O bond critical point. Values in parentheses include D2 correction.
Table 3. QTAIM parameters ρ r , 2 ρ r , (H/ ρ ) r , and (G/ ρ ) r at the Li-C bond critical points for the shortest Li-C distances, for Li adsorbed on graphene, with and without vacancies, and GO. For adsorption on GO, these properties are also calculated at the Li-O bond critical point. Values in parentheses include D2 correction.
MetalsSupportVacanciesDistances (Å)QTAIM properties (a.u.)
Li-C [Li-O] ρ r 2 ρ r ( H / ρ ) r ( G / ρ ) r
LiGraphene02.160.0250.1330.0951.224
(2.22)(0.023)(0.117)(0.101)(1.192)
1(1.98)(0.044)(0.197)(< 0.01)(1.111)
2(2.12)(0.026)(0.132)(0.094)(1.185)
GO02.330.0170.0790.1271.020
(2.26)(0.021)(0.093)(0.107)(1.000)
[(1.84)][(0.039)][(0.276)][(0.136)][(1.569)]
Li3Graphene02.2810.02210.13610.18811.3251
(2.30)(0.018)(0.093)(0.142)(1.147)
1(2.14)(0.024)(0.131)(0.118)(1.228)
2(2.12)(0.025)(0.128)(0.099)(1.162)
GO02.3620.0160.0780.1571.057
(2.30)(0.018)(0.089)(0.143)(1.068)
[(1.82)][(0.042)][(0.294)][(0.186)][(1.565)]
Li5Graphene02.240.0180.0920.1431.131
(2.28)(0.018)(0.095)(0.153)(1.155)
1(1.99)(0.043)(0.188)(< 0.01)(1.098)
2(2.10)(0.030)(0.153)(0.057)(1.209)
1 Values are for (3, +1) cp. 2 Bond critical point for dLi-C = 2.31 Å was not revealed.
Table 4. QTAIM parameters ρ r , 2 ρ r , (H/ ρ ) r , and (G/ ρ ) r at the Na-C bond critical points for the shortest Na-C distances, for Na adsorbed on graphene, with and without vacancies, and GO. For adsorption on GO, these properties are also calculated at the Li-O bond critical point. Values in parentheses include D2 correction.
Table 4. QTAIM parameters ρ r , 2 ρ r , (H/ ρ ) r , and (G/ ρ ) r at the Na-C bond critical points for the shortest Na-C distances, for Na adsorbed on graphene, with and without vacancies, and GO. For adsorption on GO, these properties are also calculated at the Li-O bond critical point. Values in parentheses include D2 correction.
MetalsSupportVacanciesDistances (Å)QTAIM properties (a.u.)
Na-C[Na-O] ρ r 2 ρ r ( H / ρ ) r ( G / ρ ) r
NaGraphene02.620.0150.0750.1471.105
(2.64)(0.015)(0.073)(0.149)(1.098)
1(2.47)(0.020)(0.092)(0.106)(1.054)
2(2.49)(0.020)(0.096)(0.108)(1.108)
GO02.540.0190.0870.1161.048
(2.50)(0.019)(0.092)(0.120)(1.074)
[(2.37)][(0.019)][(0.162)][(0.237)][(1.321)]
Na3Graphene02.500.0180.0880.1421.093
(2.59)(0.016)(0.078)(0.151)(1.085)
1(2.40)(0.029)(0.127)(0.047)(1.031)
2(2.63)(0.014)(0.072)(0.165)(1.083)
GO02.870.0090.0430.2771.005
(2.60)(0.016)(0.077)(0.155)(1.048)
[(2.27)][(0.023)][(0.152)][(0.237)][(1.426)]
Na5Graphene02.750.0120.0590.1891.027
(2.62)(0.014)(0.070)(0.179)(1.067)
1(2.40)(0.029)(0.126)(0.048)(1.035)
2(2.51)(0.020)(0.090)(0.110)(1.001)
Table 5. QTAIM parameters ρ r , (H/ ρ ) r , (G/ ρ ) r , and ξ r and ELF at all available Li-Li bond critical points, for Li adsorbed on graphene, with and without vacancies, and GO. Here, 2 ρ r < 0.01 . Values in parentheses include D2 correction.
Table 5. QTAIM parameters ρ r , (H/ ρ ) r , (G/ ρ ) r , and ξ r and ELF at all available Li-Li bond critical points, for Li adsorbed on graphene, with and without vacancies, and GO. Here, 2 ρ r < 0.01 . Values in parentheses include D2 correction.
SystemVacanciesLi-Li (Å)QTAIM properties (a.u.)ELF
ρ r ( H / ρ ) r ( G / ρ ) r ξ r η r
Li3/Graphene02.920.0131−0.1590.2333.370.31
2.920.0131−0.1580.2343.280.31
(2.85)(0.014)1(−0.166)(0.243)(3.25)(0.31)
(2.85)(0.014)1(−0.165)(0.243)(3.21)(0.31)
Li5/Graphene02.730.0112−0.1740.093−3.090.69
2.970.0111−0.1500.2283.200.28
3.040.0131−0.1530.2492.550.30
(2.77)(<0.01)2(−0.148)(0.044)(−2.40)(0.89)
(2.91)(0.010)2(−0.170)(0.164)(−40.60)(0.39)
(2.92)(0.012)1(−0.163)(0.211)(5.27)(0.34)
1(2.70)(<0.01)1(−0.115)(0.227)(2.24)(0.22)
(2.84)(0.011)1(−0.159)(0.215)(4.24)(0.32)
(3.18)(<0.01)1(−0.128)(0.246)(2.13)(0.21)
2(2.87)(0.010)1(−0.152)(0.214)(4.07)(0.29)
(3.37)(<0.01)1(−0.115)(0.251)(1.84)(0.18)
(3.49)(<0.01)1(−0.087)(0.233)(1.71)(0.11)
1 2 ρ r > 0 and 2 2 ρ r < 0.
Table 6. QTAIM parameters ρ r , (H/ ρ ) r , (G/ ρ ) r , and ξ r and ELF at all available Na-Na bond critical points, for Na adsorbed on graphene, with and without vacancies, and GO. Here, 2 ρ r < 0.01 . Values in parentheses include D2 correction.
Table 6. QTAIM parameters ρ r , (H/ ρ ) r , (G/ ρ ) r , and ξ r and ELF at all available Na-Na bond critical points, for Na adsorbed on graphene, with and without vacancies, and GO. Here, 2 ρ r < 0.01 . Values in parentheses include D2 correction.
SystemVacanciesNa-Na (Å)QTAIM properties (a.u)ELF
ρ r ( × 10 ) ( H / ρ ) r ( G / ρ ) r ξ r η r
Na3/Graphene03.350.091−0.1030.1863.020.31
3.370.091−0.0960.1972.500.28
(3.23)(0.10)1(−0.109)(0.205)(2.60)(0.30)
(3.24)(0.10)1(−0.108)(0.207)(2.51)(0.29)
1(3.82)3(0.04)2(−0.124)(0.054)(−3.59)(0.65)
2(3.35)(0.08)1(−0.099)(0.190)(2.75)(0.28)
(3.38)(0.08)1(−0.093)(0.198)(2.38)(0.26)
Na3/GO03.340.091−0.0970.1842.840.30
3.430.081−0.0880.1952.330.26
(3.19)0.052(−0.151)(0.137)(−17.51)(0.28)
Na5/Graphene03.660.062−0.1140.061-4.710.71
4.100.061−0.1040.11426.060.43
(3.54)0.072(−0.121)(0.029)(−2.71)(0.93)
(4.05)(0.07)1(−0.110)(0.141)(8.04)(0.34)
1(3.55)(0.06)1(−0.113)(0.115)(201.62)(0.41)
(3.63)(0.06)1(−0.089)(0.151)(4.017)(0.28)
2(3.34)(0.07)1(−0.084)(0.220)(1.84)(0.20)
(3.41)(0.07)1(−0.090)(0.182)(2.70)(0.26)
1 2 ρ r > 0 and 2 2 ρ r < 0. 3 Bond cp for dNa-Na = 3.77 Å was not revealed.

Share and Cite

MDPI and ACS Style

Dimakis, N.; Salas, I.; Gonzalez, L.; Vadodaria, O.; Ruiz, K.; Bhatti, M.I. Li and Na Adsorption on Graphene and Graphene Oxide Examined by Density Functional Theory, Quantum Theory of Atoms in Molecules, and Electron Localization Function. Molecules 2019, 24, 754. https://doi.org/10.3390/molecules24040754

AMA Style

Dimakis N, Salas I, Gonzalez L, Vadodaria O, Ruiz K, Bhatti MI. Li and Na Adsorption on Graphene and Graphene Oxide Examined by Density Functional Theory, Quantum Theory of Atoms in Molecules, and Electron Localization Function. Molecules. 2019; 24(4):754. https://doi.org/10.3390/molecules24040754

Chicago/Turabian Style

Dimakis, Nicholas, Isaiah Salas, Luis Gonzalez, Om Vadodaria, Korinna Ruiz, and Muhammad I. Bhatti. 2019. "Li and Na Adsorption on Graphene and Graphene Oxide Examined by Density Functional Theory, Quantum Theory of Atoms in Molecules, and Electron Localization Function" Molecules 24, no. 4: 754. https://doi.org/10.3390/molecules24040754

APA Style

Dimakis, N., Salas, I., Gonzalez, L., Vadodaria, O., Ruiz, K., & Bhatti, M. I. (2019). Li and Na Adsorption on Graphene and Graphene Oxide Examined by Density Functional Theory, Quantum Theory of Atoms in Molecules, and Electron Localization Function. Molecules, 24(4), 754. https://doi.org/10.3390/molecules24040754

Article Metrics

Back to TopTop