Arxiv:2005.03854V1 [Cond-Mat.Mtrl-Sci] 8 May 2020 by Insulating fluorographene [46]
Total Page:16
File Type:pdf, Size:1020Kb
Tuning Hydrogen Adsorption and Electronic Properties from Graphene to Fluorographone Gabriel R. Schleder ,1, ∗ Enesio Marinho Jr ,1, y Douglas J. R. Baquiao ,1 Yuri M. Celaschi ,1 Felipe Gollino ,2 Gustavo M. Dalpian ,1 and Pedro A. S. Autreto 1, z 1Federal University of ABC (UFABC), 09210-580 Santo Andr´e,S~aoPaulo, Brazil 2S~aoCarlos Institute of Chemistry, University of S~aoPaulo, 13566-590 S~aoCarlos, S~aoPaulo, Brazil Graphene functionalization by hydrogen and fluorine has been proposed as a route to modulate its reactivity and electronic properties. However, until now, proposed systems present degradation and limited hydrogen adsorption capacity. In this study, combining first-principles calculations based on density functional theory (DFT) and reactive molecular dynamics, we analyze the tuning of hydrogen adsorption and electronic properties in fluorinated and hydrogenated monolayer graphenes. Our results indicate that fluorine adsorption promotes stronger carbon{hydrogen bonds. By changing the concentration of fluorine and hydrogen, charge density transfer and electronic properties such as the band gap and spin-splitting can be tailored, increasing their potential applicability for electronic and spintronic devices. Despite fluorine not affecting the total H incorporation, the ab initio molecular dynamics results suggest that 3% fluorinated graphene increases hydrogen anchoring, indicating the hydrogenated and fluorinated graphenes potential for hydrogen storage and related applications. I. INTRODUCTION of graphene-oxide in electrochemical studies and electro- analytical applications [32]. Graphene is a one-atom-thick layer of sp2 bonded car- Another important functionalization strategy used bon atoms, which forms a two dimensional (2D) honey- is adsorbing hydrogen and/or fluorine atoms [33{37]. comb lattice. This nanomaterial was first synthesized by Graphane is defined as a theoretical nonmagnetic semi- Novoselov and Geim [1], and since then has inspired the conductor composed by 100% hydrogenation of graphene, research and development of new 2D materials [2{17] for resulting in a CH stoichiometry [38]. In addition, gra- applications in nanoscale electronics. Due to its unique phone is the midpoint between graphene and graphane in and remarkable electronic [18], thermal [19] and mechan- which the graphene sheet is only one-side hydrogenated ical properties [20], graphene has been applied in several [38, 39]. Experimental synthesis of graphone requires a technological areas, such as flexible display screens [21], back-side stabilization in order to overcome the thermo- spintronic devices [22], solar cells [23], and electrochemi- dynamic and kinetic instabilities at room temperature cal energy storage [24]. [40]. Usually this stabilization is carried out by metal However, single-layer pristine graphene can be consid- surfaces such as Ir or Ni(111) which can also control the ered a zero-gap semiconductor (semimetal), which hin- hydrogen atoms desorption (dehydrogenation) [41, 42]. ders its application in technologies that require a semi- During the hydrogenation or fluorination processes, conducting band gap such as optical and electronic de- the planar sp2 structure of the pristine graphene changes vices. To overcome this issue, many experimental strate- to a tridimensional structure based also on sp3 hy- gies have proposed the functionalization of graphene by bridization [43]. These other configurations of the reactions with organic and inorganic molecules, as well as hydrogenated-graphene are classified according to the an- chemical modifications of graphene surface [25]. Doping gles and the atomic positions, and are so-called chair or of graphene by ad-atoms adsorption, introducing point boat configurations [34]. defects, applying external strains, or combinations of Recently, graphane was synthesized in a stoichiometric these, are some potential routes to open the band gap proportion similar to theoretical predictions, presenting in graphene [26{28]. intensive luminescence in the visible spectra and highly The adsorption of various chemical species on graphene hydrophobic behavior [44]. Besides optoelectronic appli- is a widely studied subject [25, 29]. Mohamed and co- cations [45], this can lead to chemically isolated struc- workers investigated the application of nitrogen-doped tures where conducting graphene devices are surrounded graphene as a suitable photocatalyst [30]. Poh et al. arXiv:2005.03854v1 [cond-mat.mtrl-sci] 8 May 2020 by insulating fluorographene [46]. proposed a method for the preparation of halogenated graphene by thermal exfoliation in gaseous halogen at- Fluorographane, which is the functionalized graphene mospheres [31]. Chen and co-workers reviewed the prop- by H and F incorporated in mixed atmospheres, was erties of oxidized graphene and suggested the application experimentally obtained [47] and showed heterogeneous electron transfer properties. The hydrophobic properties can be tuned under different atmosphere fluorine con- centrations. Fluorographane formation process and its ∗ [email protected] properties are still subject to investigation. It is already y [email protected] known that fluorine atoms do not tend to form clusters z [email protected] [37], but in high concentrations, they can cause dam- 2 a 10.3 Å b Fluorine Graphene 10.3 Å Hydrogen c Graphone d Fluorographene e Fluorographone FIG. 1. Configurations of graphone, fluorographene and fluorographone: (a) 4×4 fluorographone supercell scheme (solid line) and unit cell (dashed line); (b) pristine graphene sheet with fluorine and hydrogen atmospheres used in the molecular dynamics simulations; 4×4 supercells of (c) graphone, (d) fluorographene, and (e) fluorographone. Carbon, fluorine, and hydrogen atoms are represented in green, blue, and orange, respectively age to the graphene sheet [48]. The authors associate PBE approximation [56]. Although this functional ap- increased hydrogen incorporation to the presence of fluo- proximation is known to underestimate electronic band rine. Hydrofluorinated Janus graphene, a functionalized gaps by approximately 50% [57], energy differences tend graphene with H and F adsorbed on opposite sides, has to be well described presenting a weighted total mean ab- been proposed as a promising system in terms of piezo- solute deviation of 0.23 eV [58]. The kinetic energy cutoff electricity, reaching values one order of magnitude higher for the plane wave expansion was 400 eV. For the adsorp- than single element adsorption [49], with potential ad- tion study of H and/or F on graphene, we have consid- sorption properties [50]. ered 4×4 supercells containing 32 carbon atoms, shown To obtain physical insights regarding the challenge of in Figures 1a,c-e. Binding sites are determined by opti- reducing the material damage while further improving mization of atomic positions performed using conjugate the hydrogen storage capacity, in the present work we gradient algorithm until the Hellmann-Feynman forces investigate the hydrofluorinated Janus graphene, which are smaller than 0.01 eV A˚−1. We used 4×4 k-meshes. here we describe as fluorographone. Based on ab initio To avoid spurious interactions between graphene sheets, calculations and reactive molecular dynamics, we system- a 10 A˚ vacuum was added. atically studied the modulation of hydrogen adsorption We investigated the hydrogen and fluorine incorpora- and the corresponding variation in physical properties tion processes on graphene using fully atomistic reactive of the hydrogenated and/or fluorinated graphene mate- molecular dynamics (MD) simulations with the ReaxFF rials. We compared the electronic and binding proper- [59] reactive force field as implemented in lammps [60] ties among different H and F functionalized systems and package. ReaxFF force fields allow us to describe bond our findings suggest that simultaneous fluorine adsorp- creations and dissociations and has been applied with tion improves hydrogen incorporation capabilities. success to hydrogenation and fluorination of graphene [48, 61]. The simulated systems are composed of a 100×100 A˚2 II. COMPUTATIONAL DETAILS graphene sheet supported on a square border of 15 A˚ and hydrogen and fluorine atoms on the bottom and top re- We performed ab initio calculations [51] based on spin- gions, respectively, confined in the exposed graphene area polarized density functional theory [13, 52, 53] as im- consisting of 11,040 C atoms. In this way, the effective plemented in vasp package [54]. The interaction of va- area used for simulation of hydrogenation/fluorination lence electrons with ionic cores was described by projec- process is then the middle 85×85 A˚2 area of the graphene tor augmented wave (PAW) pseudopotentials [55]. For sheet, shown in Figure 1b. the exchange-correlation functional, we used the GGA- We have randomly placed hydrogen/fluorine atoms in 3 a b c d Adsorption site Adsorption site Adsorption site 0 1 2 3 4 5 6 7 8 0 1 2 3 4 0 1 2 3 4 5 6 7 8 9 Graphene without F -2 0 0 0 Graphene with 3% F -0.5 -2.5 -0.4 -1 -1 -3 -2 -0.8 -1.5 -3.5 -3 -1.2 -2 -4 -4.09 -1.36 -2.5 -4 -1.6 -4.5 -3 -2.86 Binding energy per H atom [eV] Binding energy -5 5 2 per H atom [eV] Binding energy -3.5 6 3 1 0 10 20 30 40 50 12 8 7 4 Hydrogen coverage [%] 12 4 35 7 6 3 Fluorine Hydrogen Hydrogen Fluorine FIG. 2. Binding energy for adsorbing a H into the upper surface of (a) single-fluorinated graphene, (b) single-hydrogenated graphene, (c) single-fluorinated and single-hydrogenated graphene. The structures represented below describe the adsorption site used to adsorb the probe H atom in each system. (d) Binding energies per H atom as function of hydrogen coverage in graphene without F or with 3% of F adsorbed. different concentrations on the lower/upper surface of III. RESULTS AND DISCUSSION the graphene sheet, respectively. For fluorine adsorp- tion analyses an initial concentration of 7.5% was used, First, we evaluate the most stable adsorption sites corresponding to 828 F atoms, and the F atoms were for hydrogen incorporation, and how this incorporation added under the graphene square area. This concentra- changes considering fluorine co-adsorption.