Structural and electronic properties of Li_n, B_n, N_n and O_n (n=1-4) clustering on graphene: Density functional theory calculations with dispersive forces correction.
الموضوعات : Journal of Theoretical and Applied PhysicsMalika Gallouze 1 , Mahrez Drir 2 , Abdelhafid Kellou 3
1 - Laboratory of Theoretical Physics and Didactics (LPTD), Physics Faculty, University of Science and Technology-Houari Boumediene (U.S.T.H.B), Algiers, Algeria.
2 - Laboratory of Theoretical Physics and Didactics (LPTD), Physics Faculty, University of Science and Technology-Houari Boumediene (U.S.T.H.B), Algiers, Algeria.
3 - Laboratory of Theoretical Physics and Didactics (LPTD), Physics Faculty, University of Science and Technology-Houari Boumediene (U.S.T.H.B), Algiers, Algeria.
الکلمات المفتاحية: DFT, Adsorption, Cluster, Physical Properties, Van de Waals,
ملخص المقالة :
We have performed density functional theory calculations to study the adsorption of Li, B, N, and O clusters on graphene to clarify the interaction of these atoms with a graphene sheet. The stable structure, the adsorption energy, and the density of states (DOS) of atom-graphene systems are calculated. The obtained results of the structural, adsorption and interaction energies show that light atoms on graphene show physisorption/chemisorption mechanisms with and without clusteringoccurrence. More interestingly, the formation of diatomic molecules is noticed. For example, Li atoms prefer to be physisorbed where as B atoms prefer the clustering configuration on the graphene sheet. For N and O atoms, the results show that the formation of molecules on a graphene sheet is preferred. The densities of states calculations have been performed to study the electronic properties of the studied systems. The results show that electronic structures are affected by the increasing concentration and the nature of adsorbed light atoms.
Structural and electronic properties of Lin, Bn, Nn and On (n=1-4) clustering on graphene: Density Functional Theory calculations with dispersive forces correction.
M. Gallouze, M. Drir and A. Kellou
Laboratoire de Physique Théorique et de didactique (LPTD), Faculté de physique, Université des Sciences et de Technologie Houari Boumediene (U.S.T.H.B),
BP 32 El Alia, 16111 Bab-Ezzouar, Alger, Algérie
Abstract :
Density Functional Theory calculations are performed to study the adsorption of Li, B, N, and O clusters on graphene and to clarify the interaction of these atoms with a graphene sheet. The stable structures, the adsorption and interaction energies and the Density of States (DOS) of adatom-graphene systems are calculated. The obtained results show that light atoms on graphene are governed by physisorption/chemisorption mechanisms with and without clustering occurrence. More interestingly, the formation of diatomic molecules is noticed. For example, Li atoms prefer to be physisorbed whereas B atoms prefer the clustering configuration on the graphene sheet. For N and O atoms, the molecule formation is preferred to the clustering one. The electronic structure of studied systems is analyzed by performing densities of states calculations. The results show that electronic structures are affected by the increasing concentration and the nature of adsorbed light atoms.
Keywords: Graphene, Adsorption, Clustering, Electronic properties, Van de Waals, DFT.
Corresponding author: mgallouze@usthb.dz
1. Introduction:
Graphene, a single-atom thick carbon nanosheet, has emerged as a versatile material with distinctive properties originating from its two-dimensional (2D) honeycomb lattice [1-4]. Graphene is the basic material to form nanotubes, fullerenes and graphite. It attracts attention not only as part of nanocarbons but also for its own interesting electronic, mechanic, and optical properties [5], which make it a good candidate for many applications. Therefore, there have been many theoretical and experimental attempts to apply metal-embedded graphene to gas storage, sensing materials, catalyst, biodevices, low-dimensional transport and magnetic structures, and other applications [6-11]. This is due to the combination of the high chemical property of metal atoms and the large surface area of graphene. The graphene and nanoscale carbons are particularly suitable for the stabilization of metal particles while maintaining excellent mass transfer properties, which is important in catalytic studies. In electronic applications, materials with a tunable bandgap are very interesting. Therefore, several studies have reported possible graphene-like structures [12-13] by reducing the graphite thickness and achievement of monolayer and multi-Layer graphene sheets [14]. In other hand, doped-graphene systems show small bandgap by breaking the symmetry. In addition, their electronic structure and stability are highly affected by graphene interaction with atomic defects. For example, the doping graphene with N atoms is generally suggested for this kind of applications because of its high efficiency.
However, there is a major obstacle for metal-graphene complex since metal atoms prefer to aggregate without difficulty [15-16]. For example, in Li ion batteries [17], it was reported that the Li cluster formation on graphene declines the charge/discharge capacity of a battery when graphene is used as an anode material. Knowing that the interaction of atoms with graphene is a feasible way to open an energy gap for graphene, it is still difficult to understand the stability, the structural as well as the electronic and the magnetic properties of such systems. In spintronic devices, the application of graphene must include a ferromagnetic spin injector. It is then important to study the atom clustering on graphene knowing that clustering could greatly affect the long-term device stability. Recently, small ferromagnetic conductors have been used as spin injectors in graphene-based spintronics devices.
In this work, we study Li, B, N, and O clustering on graphene using the Density Functional Theory (DFT) and the Generalized Gradient Approximation (GGA). To take into account the Van der Waals (VdW) interaction effects, we have added Dispersive Force correction (DF) to the standard DFT functional. The paper is organized as follows. In Section 2, we present a brief description of the computational methods based on the Density Functional Theory (DFT). In Section 3, we present the obtained results of Li, B, N, and O clustering on graphene. The obtained results are discussed and compared with available data in the literature. The concluding remarks are given in Section 4.
2. Density Functional Theory (DFT) calculations
All calculations of the present study were performed using the framework of the Density Functional Theory (DFT) [18]. The Quantum Espresso program package was used to study the stability of nX clusters (X=Li, B, N, O; n =1-4) adsorption on a graphene sheet [19]. We used a plane-wave basis set for the expansion of the single particle Kohn–Sham wave functions. The Generalized Gradient Approximation (GGA) functional of Perdew, Burkeand Ernzerhof (PBE) [20] was employed in a spin-polarized scheme including the Ultra-Soft Pseudo-Potential (USPP) of Vanderbilt [21] with nonlinear core corrections. For the convergence, relaxation parameters were performed over the Monkhorst Pack meshes of 8x8x1 k-points along the three basic axes in reciprocal space. A plane-wave basis up to a kinetic cut-off of 400 eV was selected. All systems are allowed to fully relax using a Broyden–Fletcher–Goldforb–Shanno (BFGS) scheme. The convergence of the force between atoms is lower than 0.004eV/Å and the maximum displacement is 0.002Å. In our calculations, the Van der Waals (VdW) interactions between the adsorbed nX clusters and graphene are taken into account using the Dispersive Force (DF) correction (DF). Dispersion Forces arise from the mutual dynamic polarization of the electronic clouds of two or more atoms. London developed a formula to explain the intermolecular attractions that exist between nonpolar molecules [22]. DF correction models the long-range attractive potential of Dispersive Forces in terms of pair-wise atomic potential proportional to C6R-6, where R is the interatomic distance and C6 is the interaction parameter (0.75 in this study).
The graphene is described using the supercell approach where a single layer is separated from its image along the z-axis by 27 Å, which ensures a convergence precision of the total energy within10−5eV/atom. For this purpose, a 4x4 graphene supercell (32 C atoms) was considered for one, two, and three atoms defect adsorption. For the single atom, the adsorption was performed in three different sites: Top (T), Hollow (H) and Bridge (B), as shown in Fig. 1. The Top site is directly above a C atom; the Hollow site is at the center of a hexagon formed by six C atoms and the Bridge site is at the middle of the C-C bound. For two and three atoms adsorption, several configurations were considered taking into account previous equilibrium positions of one atom and only favorable configurations were considered. Accordingly, nX clusters-graphene systems were fully relaxed to equilibrium geometries.
3. Results and discussion
3.1 Structural and energetic properties
Equilibrium structures and equilibrium distances between nearest-neighbor adsorbed atoms and graphene atoms are reported in Fig. 2 and Fig. 3.
a. A single X atom adsorption on graphene
In the single atom case, the results show that Li and B atoms do not bound to a graphene surface as shown in Fig. 2. The B atom prefers the bridge site with the minimum adsorption energy of -1.188 eV and the equilibrium distance to the graphene surface of 1.69 Å. However, the Li atom prefers the hollow site with the minimum adsorption energy of -0.759 eV and the equilibrium distance to the graphene surface of 1.85 Å. For the graphene, interatomic distances between C atoms are not affected and remain unchanged in the presence of Li and B. Taking into account the above results, Li and B atoms prefer to be physisorbed. These results were confirmed by previous studies [23-26]. Several authors have found that the most favorable adsorption site for the Li atom is the hollow. The equilibrium distance to the graphene is 1.84 Å in ref [23], 1.71 Å in ref [24] and 1.67 Å in ref [25]. In addition, Li et al. [26] have indicated that the most favorable site for B is the bridge site at 1.69 Å from the graphene surface. However, the N and O atoms prefer to be bounded to two carbon atoms in the bridge site forming a triangle with two carbon atoms below, as shown in Fig. 2. The distance between N and O atoms and the graphene is 1.61 Å and 1.51Å, respectively. The presence of N and O has affected the planer configuration of C atoms at the adsorption site. Interatomic distances between N-C and O-C are 1.55 Å and 1.51 Å, respectively. These distances are smaller than Li and B ones. Accordingly, N and O have reduced the overlap of the electron cloud between two carbon atoms at the adsorption site. In addition, adsorption energies are predicted to be -3.564 eV and -3.597 eV for N and O, respectively. Consequently, from the perspective of bond length and the reported values of the adsorption energy, N and O are strongly bounded to two carbon atoms. These facts mean that N and O are chemisorbed.
The adsorption of N and O atoms has induced a deformation of the graphene lattice. However, Li and B are rather far from the graphene and the its lattice is not deformed. The results concerning the adsorption of O on graphene are in good agreement with ref [27].
b- The adsorption of nX atoms on graphene
By increasing their coverage, Li atoms have preferred the hollow site as show in Fig. 3. The relaxed configurations are represented with equilibrium atomic distances. As reported in Fig. 3, Li-graphene distances increase as the Li coverage increases. However, Li-Li distances remain unchanged (2.85 Å). This value is slightly smaller than the distance between the Li-Li distance in the bulk bcc-Li (3.01 Å [28]). This means that a formation of Li-cluster is not favored. The distance between Li atoms and graphene is smaller than Li-Li one. This means that the physorption of Li atoms on graphene is significant without deformation of the graphene sheet. Moreover, the adsorption site of the fourth Li atom is localized above the other three atoms. Li-Li distances on the same layer are identical to that of three Li atoms adsorbed on graphene with the distance between the fourth Li atom and others atoms is 3.05 Å. This behavior is the consequence of the Li clustering occurrence. This result is confirmed by G. Yang et al. [29] which have indicated that the stability of Li on graphene is due to the Li cluster configuration with increasing Li content.
For two B atoms, the first one is bounded to carbon atoms in bridge sites and the second one is vertically bounded to the first one. The B-B bond is 1.59 Å and the adsorption energy is -2.41 eV/B. For three and four B atoms, the equilibrium structure has indicated that two B atoms are bounded to the graphene sheet and other B atoms are bounded to each other in a triangular form. It can be concluded that the adsorption of B atoms on graphene has preferred the clustering configuration.
For two N atoms, the relaxed structure shows that two atoms are bounded to two surface carbon atoms while remaining bounded to each other with a bond length of 1.27 Å. The adsorption energy is -5.576 eV/N. Moreover, when three and four N atoms are adsorbed, two of them form a free dimer with an interatomic distance of 1.11 Å.
For two O atoms, the epoxy group is formed in the vicinity of each other. However, for three O atoms adsorbed on graphene, two atoms have formed a free dimer with a bond length of 1.25 Å. The remaining O atom is chemisorbed on the graphene sheet. When four O atoms are adsorbed, the relaxed structure has shown the formation of a free trimer (ozone) with a bond length of 1.31Å and an angle of 117.21°. This indicates that oxygen atoms have a strong tendency to form a molecular structure on graphene.
The bond length values of free 2N, 2O, and 3O are similar to experimental molecular counterparts which are 1.098 Å, 1.208 Å, and 1.28 Å [30], respectively. The adsorption of B, N, and O atoms has stretched the C-C bound to 1.52, 1.55 Å, and 1.51 Å, respectively. This means that the planar form of graphene was altered. As a consequence, the graphene sheet was distorted and atoms settle into a valley created by them. O atoms have formed an epoxy group as reported in available theoretical results [31, 32].
c- Formation, adsorption and interaction energies
To investigate the stability of nX atoms upon adsorption on a graphene sheet, we have calculated the adsorption energy of X atoms referenced to their bulk energies. We have also determined the interaction energy between the cluster and the graphene sheet. The interaction of nX atoms with a graphene sheet is considered stable if the adsorption energy absolute value is larger than the X bulk binding energy. Moreover, the nX atoms form a bond with the graphene sheet if the interaction energy is negative (Eint< 0).
The enthalpy DHf , the adsorption energy Eads and the interaction energy Eint are defined as follows:
Where Egra-nX is the total energy of graphene with nX atoms (X= Li, B, N, and O), and are the energy of X and C in bulk structure, respectively, n is the number of adsorbed atoms.
Where EX is the total energy of an isolated X atom, and Egra is the total energy of graphene sheet. The and are the sublattice energies of nX and graphene, respectively.
In Fig. 4, we have represented adsorption and interaction energies for nX atoms on a graphene sheet versus the number of adsorbed atoms. As it can be shown in Fig. 4(a), adsorption energies for any X concentration have negative values. It is interesting to notice that the absolute value increases with the number of adsorbed B and N atoms. However, adsorption energies of Li and O remain unchanged. With increasing adsorbed atoms, obtained results indicate that Li and O have preferred to be physisorbed. Clustering and molecular configurations are noticed for Li and O, respectively. These results are confirmed by configurations of Fig. 3. Moreover, adsorption energies of Li and B are greater than bulk binding energies of Li (-1.63 eV/Li) and B (-5.81 eV /B). For Li, obtained results are in good agreement with those of Dimakis et al. [33]. For O and N, adsorption energies are smaller than molecular binding energies of O (-2.60 eV/O) and N (-4.92 eV/N).
To investigate the interaction of adsorbed atoms and the graphene sheet, we have represented X-C interaction energies in Fig. 4(b). As shown, all interaction energies are attractive except for two and four N atoms. In these cases, interaction energies are repulsive. In addition, the Li-C bound is weaker when compared to other bounds. More details about formation energies and X-X interaction energies are given in Table 1. As indicated, all formation energies are negative indicating the stability of present systems. Also, all energies are low then the formation energy of monolayer of graphene which is about -0.55 eV /at. Thus, systems with adsorbed X atoms are more stable than the clean graphene and the stability increases with increasing number of X atoms.
For the X-X interaction, the Li-Li interaction is closer to the Li-C one. However, the Li cluster interacts weakly with the graphene sheet without any distortions (see Fig. 3). It can be deduced that the Li cluster is physisorbed on the graphene sheet, as concluded above. In contrast, the B-B interaction is weaker than the B-C one, which explains the strong distortion of the graphene sheet. This fact also leads to the formation of B clusters. Concerning N and O, N-N and O-O interactions are stronger than N-C and O-C ones, respectively. The N-N interaction for 2N and 4N atoms adsorbed on the graphene sheet is greater than the binding energy of N2 (4.92 eV/at). This confirms the occurrence of the molecular configuration of N and O on the graphene sheet. The detailed description of N and O adsorption is given below.
For 3N adsorbed on the graphene sheet, the N-N interaction is closer to the binding energy of N2. This indicates that one N atom is bounded to one C atom. Remaining two atoms have formed a molecule which interacts weakly with graphene. For 4N adsorbed atoms, the N-N interaction energy is close to the N-N one of 2N adsorbed atoms. Consequently, one N2 molecule is freely formed on the graphene sheet. For two O adsorbed atoms, the O-O interaction energy is close to zero, indicating equilibrium between the attractive and the repulsive interaction. These two atoms have formed an attractive interaction with C atoms (epoxy group), as shown in the Fig. 4(b). However, interaction energies of 3O and 4O have close values and are comparable to the O2 binding energy (2.6 eV/at). These results confirm the formation of O2 and O3(ozone) molecules, as shown in the Fig. 3.
3.2 The densities of states (DOSs)
In this section, we examine the effects of the nX atoms on electronic properties of the graphene system. The Fig. 5 shows the DOS of nLi-C, nB-C, nN-C, and nO-C (n=1-4). For comparison, the total DOS of pristine graphene is also represented. As shown in all curves, graphene is a zero-bandgap material. However, nX atoms adsorbed on graphene change completely the electronic character of the system. DOSs near the Fermi level are significantly altered involving a charge transfer between the cluster and pristine graphene. The zero-bandgap character has disappeared in all cases.
In nLi graphene systems, the total DOS represented in Fig.5(a) shows that the nLi adsorption leads to an upshift of the Fermi level into the conduction band due to a charge transfer from Li to graphene without any significant disruption of the p electronic system or graphene corrugation. All of these structures turn to a metallic character. These results are similar to those obtained by G. Yang et al. [29].
In the case of nB graphene systems, all total DOS represented in Fig.5(b) have the same behavior of Li-systems. Beyond the two atoms, the zero-bandgap character has disappeared and B-systems become metallic.
In the case of nN graphene systems, the total DOS represented in Fig.5(c) shows that the curves of N and 3N are similar. This is due to the formation of the N2 molecule. The Dirac cone is shifted toward a higher binding energy with respect to pristine graphene which has the Dirac point at the Fermi energy, Ef. The conduction band appears below Ef. The zero-bandgap character has disappeared and systems become metallic. The DOS of 2N and 4N atoms are similar and more electronic states are fulfilled. The form of the DOS at the Fermi level suggests that the adsorption of 2N and 4N has pseudo-finite bandgap (semiconducting character).
For nO, the total DOS represented in Fig. 5(d) show that the curves of O and 2O adsorbed on graphene are similar. In these cases, finite energy gap arises at K point which leads to the change of the electronic states from a zero bandgap to a finite one. This result is confirmed by experiments [27]. The systems with O2 or O3 molecules have a metallic character as found experimentally [34] for O2 adsorbed on the doped-graphene. It is important to notice the high conductivity character of 3O when compared to the 4O one.
To study the behavior of the atom cluster adsorbed on the graphene sheet, density of states for different concentrations are shown in Fig. 6. The results show that for all atoms and for each concentration, curves upshift in the conduction band where the metallic character is pronounced except for the case of one and two oxygen. In these cases, the curves are very distinct. As shown, the adsorption of one O atom leads to a pseudo gap formation, whereas the adsorption of 2O atoms leads to an insulating character.
4. Conclusion
We have performed density functional theory calculations to investigate the structural and electronic properties of Li, B, N, and O atoms adsorbed on a graphene sheet. The results have shown an interesting behavior of clustering or molecule formation on a graphene sheet. To confirm these trends, relaxed geometries, adsorption and interaction energies, and electronic structures were calculated.
In the case of a single atom, Li and B atoms have preferred to be physisorbed without deforming the graphene sheet. However, the N and O atoms have preferred to be chemisorbed inducing a deformation of the graphene sheet. By increasing the number of X atoms, Li has preferred to be physisorbed with a cluster formation with a number greater than three. B has preferred to be chemisorbed with the clustering configuration. For N and O, a strong tendency to form a molecular structure on graphene was observed. It is important to note that the adsorption of B, N, and O atoms modifies the planar form of graphene.
DOS calculations for all adsorbed atoms have shown that the zero-bandgap character was disappeared in all cases. The adsorption of Li and B atoms has induced the loss of the zero-bandgap character which a metallic character appearance. However, the adsorption of 2N and 4N leads to an insulating (semiconducting) character. In addition, systems with O2 or O3 molecules have shown a metallic character. It is important to notice the high conductivity character of 3O when compared to the 4O one. The adsorption of one O atom has led to a pseudo gap formation, whereas the adsorption of 2O atoms has led to an insulating behavior.
References:
[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, A. A. Firsov, Electric Field Effect in Atomically Thin Carbon Films. Science 306 (2004) 666–669.
[2] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, A. A. Firsov, Two-dimensional gas of massless Dirac fermions in graphene. Nature 438 (2005) 197–200.
[3] C. Lee, X. D. Wei, J. W. Kysar, J. Hone, Measurement of the Elastic Properties and Intrinsic Strength of Monolayer Graphene. Science 321 (2008) 385–388.
[4] T. Zhang, Q. Z. Xue, S. Zhang, M. D. Dong, Theoretical approaches to graphene and graphene-based materials. Nano Today. 7 (2012) 180–200.
[5] K. Nakada and A. Ishii, DFT Calculation for Adatom Adsorption on Graphene, Graphene Simulation. in: J.R. Gong (Eds.), Graphene simulation, Technology Inc., Croatia (2011) pp. 3-20.
[6] H. Lee, J. Ihm, M. L. Cohen, S. G. Louie, Calcium-decorated graphene-based nanostructures for hydrogen storage. Nano letters 10 (2010) 793-798.
[7] X. Y. Lianga, N. Dinga, S. Nga, C. L. Wua, Adsorption of gas molecules on Ga-doped graphene and effect of applied electric field: A DFT study », Applied Surface Science 411 (2017) 11–17.
[8] S. Shanmugam, S. Nachimuthu, V. Subramaniam, Enhanced Li+ ion adsorption on pristine and defected graphene via organic radical interaction - A DFT study, Physica B 611 (2021) 412700.
[9] G. Yu, Y. Xie, Q. Ge, Q. Dai, J. Xu, H. Cao, Mechanism of ozone adsorption and activation on B-, N-, P-, and Si-doped graphene: A DFT study, Chemical Engineering Journal 430 (2022) 133114.
[10] J. Huang, L. M. Zhang, B. Chen, N. Ji, F. H. Chen, Y. Zhang, Z. J. Zhang, The synthesis of perylene-coated graphene sheets decorated with Au nanoparticles and its electrocatalysis toward oxygen reduction. Nanoscale 2 (2010) 2733–2738.
[11] T. T. Baby, S. S. Aravind, T. Arockiadoss, R. B. Rakhi, S. Ramaprabhu, Metal decorated graphene nanosheets as immobilization matrix for amperometric glucose biosensor. Sens. Actuators B 145 (2010) 71–77.
[12] M. M. Sarmazdeh, R. T. Mendi, A. Zelati, A. Boochani and F. Nofeli, First-principles study of optical properties of InN nanosheet, International Journal of Modern Physics B 30 (19), 1650117 (2016).
[13] M. Yeganeh, F. Kafi, A. Boochani, Thermoelectric properties of InN graphene-like nanosheet: A first principle study, Superlattices and Microstructures 138 (2020) 106367.
[14] F. Lahourpour, A. Boochani, S. S. Parhizgar, S. M. Elahi, Structural, electronic and optical properties of graphene‑like nano‑layers MoX2(X:S,Se,Te): DFT study,Journal of Theoretical and Applied Physics 13 (3), 191-201 (2019),
[15] Q. Sun, Q. Wang, P. Jena, Y. Kawazoe, Clustering of Ti on a C60 Surface and Its Effect on Hydrogen Storage. J. Am. Chem. Soc. 127-42 (2005) 14582–14583.
[16] S. Li, P. Jena, Combinatorial Search for Optimal Hydrogen-Storage Nanomaterials Based on Polymers. Phys. Rev. Lett. 97 (2006) 056104.
[17] A. L. M. Reddy, A. Srivastava, S. R. Gowda, H. Gullapalli, M. Dubey, P. M. Ajayan, Synthesis of Nitrogen-Doped Graphene Films For Lithium Battery Application. ACS Nano, 4 (2010) 6337–6342.
[18] W. Kohn, and L. J. Sham, Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. A 140 (1965) 1133.
[19] S. Baroni, A. Dal Corso, S. de Gironcoli, P. Giannozzi, C. Cavazzoni, G. Ballabio, S. Scandolo, G. Chiarotti, P. Focher, A. Pasquarello, K. Laasonen, A. Trave, R. Car, N. Marzari, A. Kokalj, Journal of Physics: Condensed Matter 21 (2009) 39. http://www.pwscf.org/.
[20] J. P. Perdew, K. Burke, M. Ernzerhof, Generalized gradient approximation made simple. Phys Rev Lett 77 (1996) 3865-3868.
[21] D. Vanderbilt, Soft self-consistent pseudopotentials in a generalized eigenvalue formalism. Physical Review B 41 (1990) 7892-7895.
[22] V. Barone, M. Casarin, D. Forrer, M. Pavone, M. Sambi, A. Vittadini, Role and effective treatment of dispersive forces in materials: Polyethylene and graphite crystals as test cases. J. Comput. Chem. 30 (2009) 934-939.
[23] J. Zheng, Z. Ren, P. Guo, L. Fang and J. Fan, Diffusion of Li+ ion on graphene: A DFT study. Applied Surface Science 258 (2011) 1651-1655.
[24] K. T. Chan, J. B. Neaton, and M. L. Cohen, First-principles study of metal adatom adsorption on graphene. Phys Rev B 77 (2008) 235430-12.
[25] W. Liu, C. Zhang, M. Deng, S. Cai, The structural and electronic properties of metal atoms adsorbed on graphene. Physica E 93 (2017) 265-270.
[26] J. Li, X. Wanga, K. Liu, Y. Suna and L. Chena, High hydrogen-storage capacity of B-adsorbed graphene: First-principles calculation. Solid State Communications 152 (2012) 386-389.
[27] S. Horiuchi, T. Gotou, M. Fujiwara, R. Sotoaka, M. Hirata, K. Kimoto, T. Asakai, T. Yokosawai, Y. Matsui, K. Watanabe and M. Sekita, Carbon Nanofilm with a New Structure and Property. J. Appl. Phys. 42 (2003) L107.
[28] F. Tran, R. Laskowski, P. Blaha, K. Schwarz, Performance on molecules, surfaces and solids of the Wu-Cohen GGA exchange-correlation energy functional. Phys. Rev. B 75 (2007) 115131.
[29] G. yang, X. Fan, Z. Liang, Q. Xu and W. T. Zheng, Density functional theory study of Li binding to graphene. RSC Adv. 6 (2016) 26540-26545.
[30] L. E. Sutton, Tables of Interatomic Distances and Configuration in Molecules and Ions. London, Chemical Society, 1958.
[31] Z. Sljivancanin, A. S. Milosevic, Z. S. Popovic, F. R. Vukajlovic, Binding of atomic oxygen on graphene from small epoxy clusters to a fully oxidized surface. CARBON, 2013, 54, 482– 488.
[32] A. M. Suarez, L. R. Radovic, E. Bar-Ziv and J. O. Sofo, Gate-Voltage Control of Oxygen Diffusion on Graphene. Phys. Rev. Lett. 106 (2011) 146802.
[33] N. Dimakis, I. Salas, L. Gonzalez, O. Vadodaria, K. Ruiz and M. I. Bhatti, 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 (2019) 754-28.
[34] J. Dai, and J. Yuan, Adsorption of molecular oxygen on doped graphene: Atomic, electronic, and magnetic properties. Phys. Rev. B 81 (2010) 165414-7.
Table 1: The X-X interaction energy of adsorbed atoms on graphene, Eint (in eV/at), and the formation energy, DHf (in eV/at), of C32Xn.
X | Li | B | N | O | ||||
| Eint | DHf | Eint | DHf | Eint | DHf | Eint | DHf |
C32X1 | -0.00 | -0.558 | 0.0 | -0.571 | 0.0 | -0.644 | 0.0 | -0.645 |
C32X2 | -0.57 | -0.573 | -1.56 | -0.682 | -7.26 | -0.874 | -0.003 | -0.765 |
C32X3 | -0.67 | -0.603 | -3.25 | -0.880 | -5.26 | -1.124 | -2.80 | -0.900 |
C32X4 | -0.74 | -0.658 | -3.51 | -1.062 | -7.56 | -1.369 | -3.08 | -1.046 |
Figures captions
Figure 1: The three high symmetrical adsorption sites on graphene: Hollow (H), Bridge (B), and Top (T).
Figure 2: Equilibrium relaxed structures of one X atoms adsorbed on a 4×4 graphene sheet. X=Li, B, N and O.
Figure 3: Equilibrium relaxed structures of nX atoms (n=2-4) adsorbed on a 4×4 graphene sheet.
Figure 4: (a) The adsorption energy, Eads, for nX atoms on a graphene sheet and (b)X-C interaction energy,, as function of the number of adsorbed atoms, n.
Figure 5: The total DOS of nX atoms on a graphene sheet compared to pure graphene (closed line). The Fermi energies of all systems are shifted to 0 eV.
Figure 6: The comparison of the total DOS between the nX-C systems for different concentrations.
Figure 1
|
| ||
|
| ||
|
| ||
|
|
Figure 2
| B | N | O | ||||||
|
|
|
| ||||||
|
|
|
| ||||||
|
|
|
| ||||||
|
|
|
| ||||||
|
|
|
| ||||||
|