Density Functional Theory (DFT) Study on the Ternary ... - MDPI

5 downloads 0 Views 2MB Size Report
Dec 29, 2015 - Mami Mutoh 1, Shigeaki Abe 1,*, Teruo Kusaka 1, Mariko Nakamura 2, ... [email protected] (M.M.); kusaka.teruo@gmail.com (T.K.); ...
Article

Density Functional Theory (DFT) Study on the Ternary Interaction System of the Fluorinated Ethylene Carbonate, Li+ and Graphene Model Mami Mutoh 1 , Shigeaki Abe 1, *, Teruo Kusaka 1 , Mariko Nakamura 2 , Yasuhiro Yoshida 1 , Junichiro Iida 1 and Hiroto Tachikawa 3 Received: 23 September 2015; Accepted: 21 December 2015; Published: 29 December 2015 Academic Editors: James Babb and Hyun-Kyung Chung 1

2 3

*

Graduate School of Dental Medicine, Hokkaido University, Sapporo 060-8586, Japan; [email protected] (M.M.); [email protected] (T.K.); [email protected] (Y.Y.); [email protected] (J.I.) School of Health Science, Kyushu University of Health and Welfare, 1714-1 Yoshino-cho, Nobeoka, Miyazaki 882-8508, Japan; [email protected] Division of Materials Chemistry, Graduate School of Engineering, Hokkaido University, Sapporo 060-8628, Japan; [email protected] Correspondence: [email protected]

Abstract: The ternary interaction system composed of fluorinated ethylene carbonate, denoted by EC(F), lithium ion (Li+ ) and a model of nano-structured graphene has been investigated by means of the density functional theory (DFT) method. For comparison, fluorinated vinylene carbonate, denoted by VC(F), was also used. The model of graphene consisting of 14 benzene rings was examined as a nano-structured graphene. The effects of fluorine substitution on the electronic state and binding energy were investigated from a theoretical point of view. It was found that both EC(F) and VC(F) bind to a hexagonal site corresponding to the central benzene ring of the model of the graphene surface. The binding energies of Li+ EC(F) and Li+ VC(F) to the model of graphene decreased with increasing number of fluorine atoms (n). Keywords: density functional theory; fluorinated ethylene carbonate; nano-structured graphene model; lithium ion secondary battery; surface interaction

1. Introduction There are several interactions in a lithium ion secondary battery. Desolvation, Li+ -EC Ñ Li+ + EC, is an important process in the movement of the Li+ ion from the solvent to graphene. The weaker solvation of the Li+ ion by ethylene carbonate (EC) is a positive factor as a solvent. Ethylene carbonate (EC) and its related compounds are widely used in the high dielectric solvent in the lithium secondary battery, because EC has many desirable properties as an electrolyte solvent for lithium batteries. For example, EC has a high dielectric constant and low viscosity [1–4]. The interaction of the electrode (graphite) with the solvent molecule (EC) plays an important role in the redox process in electrochemistry. Recently, Zhu et al. (2014) [5] have synthesized a new family of polyfluoroalkyl-substituted ethylene carbonates, and they have been tested as additives in lithium-ion cells. The four additives studied were 4-(trifluoromethyl)-1,3-dioxolan-2-one (TFM-EC), 4-(perfluorobutyl)-1,3-dioxolan-2-one (PFB-EC), 4-(perfluorohexyl)-1,3-dioxolan-2-one (PFH-EC) and 4-(perfluorooctyl)-1,3-dioxolan-2-one (PFO-EC). Differential capacity plots from graphite-Li cells suggest that PFO-EC is involved in solid electrolyte interphase (SEI) formation. Raman data from anodes after cycling suggest that structural disordering of graphite is reduced by the addition of

Atoms 2016, 4, 4; doi:10.3390/atoms4010004

www.mdpi.com/journal/atoms

Atoms 2016, 4, 4

2 of 8

Atoms 2016, 4, 4  2 of 10  Atoms 2016, 4, 4  2 of 10  PFO-EC, which may explain the improved cell capacity retention. For comparison, fluorinated vinylene carbonate, denoted by VC(F), was also used. Figure 1 shows the typical interaction systems, carbonate,denoted by VC(F), was also used. Figure 1 shows the typical interaction systems, such as  carbonate,denoted by VC(F), was also used. Figure 1 shows the typical interaction systems, such as  such as solvent-graphene, solvent-lithium ion and lithium ion-graphene. solvent‐graphene, solvent‐lithium ion and lithium ion‐graphene.  solvent‐graphene, solvent‐lithium ion and lithium ion‐graphene. 

   Scheme 1. Structures of ethylene carbonate (EC) and vinylene carbonate (VC).  Scheme 1. Structures of ethylene carbonate (EC) and vinylene carbonate (VC). Scheme 1. Structures of ethylene carbonate (EC) and vinylene carbonate (VC). 

   Figure 1. Schematic illustration of the types of interactionsin the graphene-ethylene carbonate Figure 1. Schematic illustration of the types of interactionsin the graphene‐ethylene carbonate (EC)‐ Figure 1. Schematic illustration of the types of interactionsin the graphene‐ethylene carbonate (EC)‐ (EC)-lithium system. lithium system.  lithium system. 

Several theoretical works on the water in carbon nanotubes and on graphene have recently been  Several theoretical works on the water in carbon nanotubes and on graphene have recently Several theoretical works on the water in carbon nanotubes and on graphene have recently been  performed  to  feature  of  Werber  et  [6]  out  been performed to the  gaindynamical  the dynamical feature the interaction. Werber et al. (2003) [6] carried performed  to  gain  gain  the  dynamical  feature  of  the  the ofinteraction.  interaction.  Werber  et  al.  al.  (2003)  (2003)  [6]  carried  carried  out  a  a  theoretical  calculation  of  the  carbon‐water  interaction  system  and  showed  the  importance  of  out a theoretical calculation the carbon-water interaction system showed the importance of theoretical  calculation  of  the ofcarbon‐water  interaction  system  and  and showed  the  importance  of  the  the  relation  between  the  contact  angle  of  a  water‐carbon  the relation between contact angle a waterdroplet  dropleton  ongraphite  graphitechanges  changes and  and the  the water-carbon relation  between  the the contact  angle  of of a  water  water  droplet  on  graphite  changes  and  the  water‐carbon  interaction energy. Rana and Chandra (2007) [7] investigated the interaction of water containing a  interaction energy. Rana and Chandra (2007) [7] investigated the interaction of water containing a interaction energy. Rana and Chandra (2007) [7] investigated the interaction of water containing a  narrow carbon nanotube. The hydrogen bond distributions of water inside and also in the vicinity of  narrow carbon nanotube. The hydrogen bond distributions of water inside and also in the vicinity narrow carbon nanotube. The hydrogen bond distributions of water inside and also in the vicinity of  the  outer  surfaces  of  the  simulation.  The  of the outer surfaces thenanotube  nanotubewere  wereobtained  obtainedby  bya  molecular dynamics  dynamics (MD)  (MD) simulation. the  outer  surfaces  of of the  nanotube  were  obtained  by  a a molecular  molecular  dynamics  (MD)  simulation.  The The  solvation mechanism of lithium ions in pure ethylene carbonate (EC) solutions was studied in a wide  solvation mechanism of lithium ions in pure ethylene carbonate (EC) solutions was studied in a wide solvation mechanism of lithium ions in pure ethylene carbonate (EC) solutions was studied in a wide  concentration  range  by  different  techniques  and  temperatures  Thus,  the interaction  of  concentration techniques and temperatures [8]. [8].  Thus, the interaction of solvent with concentration range range by by different different  techniques  and  temperatures  [8].  Thus,  the interaction  of  solvent  solvent  + with++  with carbon materials has been investigated by several groups. However, the interaction of EC‐Li carbon materials has been investigated by several groups. However, the interaction of EC-Li with carbon materials has been investigated by several groups. However, the interaction of EC‐Li   with graphene is not clearly understood.  graphene is not clearly understood. with graphene is not clearly understood.  ++ the  In  our  previous  papers  [9–11],  we  investigated  the  interaction  systems  of  graphene‐Li In In  our our  previous previous  papers papers [9–11], [9–11], we we investigated investigated the the interaction interaction systems systems of of graphene-Li graphene‐Li+, ,,  the the  graphene‐Li  atom  and  graphene‐water  and  showed  that  the  water  molecule  can  bind  an  edge  graphene-Li atom and graphene-water and showed that the water molecule can bind an edge graphene‐Li  atom  and  graphene‐water  and  showed  that  the  water  molecule  can  bind  an  edge  of  of  graphene.  of graphene. graphene.  In  In  the  the  present  present  study,  study,  the  the  density  density  functional  functional  theory  theory  (DFT)  (DFT)  method  method  [12–14]  [12–14]  was  was  applied  applied  to  to  a  a  ternary interaction system composed of fluorinated ethylene carbonate, denoted by EC(F), lithium  ternary interaction system composed of fluorinated ethylene carbonate, denoted by EC(F), lithium  + ion (Li ion (Li+) and a nano‐structured graphene. For comparison, fluorinated vinylene carbonate,denoted  ) and a nano‐structured graphene. For comparison, fluorinated vinylene carbonate,denoted  by VC(F), was also used.  by VC(F), was also used. 

Atoms 2016, 4, 4

3 of 8

In the present study, the density functional theory (DFT) method [12–14] was applied to a ternary interaction system composed of fluorinated ethylene carbonate, denoted by EC(F), lithium ion (Li+ ) and a nano-structured graphene. For comparison, fluorinated vinylene carbonate, denoted by VC(F), was also used. 2. Method of Calculations 3 of 10  AsAtoms 2016, 4, 4  a graphene system, a carbon sheet composed of 14 benzene rings (C42 H16 ) was used as the nano-structured graphene in the present study. The edge carbon atom of graphene was terminated by 2. Method of Calculations  a hydrogen atom. First, each structure of graphene, EC and VC was individually optimized using the As a graphene system, a carbon sheet composed of 14 benzene rings (C42H16) was used as the  Handy and coworkers’ long-range corrected version of B3LYP using the Coulomb-attenuating method nano‐structured graphene in the present study. The edge carbon atom of graphene was terminated  (CAM-B3LYP) with a 6-31G(d) basis set [15]. Next, the graphene-EC-Li+ ternary interaction system by a hydrogen atom. First, each structure of graphene, EC and VC was individually optimized using  was fully optimized at the same level. All density functional theory (DFT) calculations were carried the Handy and coworkers’ long‐range corrected version of B3LYP using the Coulomb‐attenuating  + ternary interaction  method (CAM‐B3LYP) with a 6‐31G(d) basis set [15]. Next, the graphene‐EC‐Li out using the Gaussian 09 program package [16].

system was fully optimized at the same level. All density functional theory (DFT) calculations were 

3. Results carried out using the Gaussian 09 program package [16].  3. Results  3.1. Structures of EC(F) and VC(F)

Figure 2 shows the optimized structures of EC(F) and VC(F) used in the present study. The 3.1. Structures of EC(F) and VC(F)  notation of n in EC(Fn) and VC(Fn) means the number of fluorine atoms in EC and VC, respectively. Figure  2  shows  the  optimized  structures  of  EC(F)  and  VC(F)  used  in  the  present  study.  The  The number n in EC(Fn) and VC(Fn) was varied as (n = 0 to 4) and (n = 0 to 2), respectively. notation of n in EC(Fn) and VC(Fn) means the number of fluorine atoms in EC and VC, respectively.  The bond lengths of C=O carbonyl and the C–O single bond of EC(F0) were 1.190 and 1.355 Å, The number n in EC(Fn) and VC(Fn) was varied as (n = 0 to 4) and (n = 0 to 2), respectively.  The bond lengths of C=O carbonyl and the C–O single bond of EC(F0) were 1.190 and 1.355 Å,  respectively. These bond lengths were slightly changed by the F-substitution: 1.078 and 1.367 Å in EC(F4).respectively. These bond lengths were slightly changed by the F‐substitution: 1.078 and 1.367 Å in  A large, different point appeared as a distance of r(C–C): 1.628 Å in EC and 1.552 Å in EC(F4). EC(F4). A large, different point appeared as a distance of r(C–C): 1.628 Å in EC and 1.552 Å in EC(F4).  The F-substitution largely decreases the C–C single bond of EC. The F‐substitution largely decreases the C–C single bond of EC. 

  Figure Figure 2. Optimized structures of EC, fluorinated ethylene carbonate (EC(F)), vinylene carbonate (VC)  2. Optimized structures of EC, fluorinated ethylene carbonate (EC(F)), vinylene carbonate (VC) and fluorinated vinylene carbonate (VC(F)), calculated at the B3LYP/6‐31G(d) level. The notation of n  and fluorinated vinylene carbonate (VC(F)), calculated at the B3LYP/6-31G(d) level. The notation of n in EC(Fn) and VC(Fn) means the number of fluorine atoms. Bond lengths are in Å.  in EC(Fn) and VC(Fn) means the number of fluorine atoms. Bond lengths are in Å. The dipole moments of EC(F) and VC(F) are given in Table 1. The dipole moment of EC in n = 0  was  5.42  Debye,  which  in  reasonable  agreement  with  experimental  value  (4.9  Debye).  The  The dipole moments of is  EC(F) and VC(F) are given in the  Table 1. The dipole moment of EC in n = 0 dipole moment decreased significantly with increasing n. In the case of EC(F), the value (n = 4) was  was 5.42 Debye, which is in reasonable agreement with the experimental value (4.9 Debye). The 1.39 D, which is lower than that of n = 0 (5.42 D). This is due to the fact that the dipole moment of the  C=O carbonyl vanishes in the presence of the C‐F dipole. Thus, the effect of the halogen atom on the  electronic  state  is  very  important  in  ethylene  carbonate.  In  the  case  of  the  VC  system,  a  similar  diminishing was found after the F‐substitution. 

Atoms 2016, 4, 4

4 of 8

dipole moment decreased significantly with increasing n. In the case of EC(F), the value (n = 4) was 1.39 D, which is lower than that of n = 0 (5.42 D). This is due to the fact that the dipole moment of the C=O carbonyl vanishes in the presence of the C-F dipole. Thus, the effect of the halogen atom on the electronic state is very important in ethylene carbonate. In the case of the VC system, a similar diminishing was found after the F-substitution. Table 1. Dipole moments of EC, EC(F), VC and VC(F) (in Debye). The number of fluorine atoms in EC is denoted by n. The experimental value of the dipole moment is given in parenthesis [17]. CAM, Atoms 2016, 4, 4  4 of 10  Coulomb-attenuating method. Table 1. Dipole moments of EC, EC(F), VC and VC(F) (in Debye). The number of fluorine atoms in  Dipole Moment EC is denoted by n. The experimental value of the dipole moment is given in parenthesis [17]. CAM,  Coulomb‐attenuating method.    n CAM-B3LYP/6-31G(d) CAM-B3LYP/6-311G(dp)

EC0 5.36 (4.9)Dipole Moment 5.42 EC(F1)1 4.76 4.80 CAM‐B3LYP/6‐31G(d) CAM‐B3LYP/6‐311G(dp)  n  EC(F2)2 3.61 3.50 EC0  5.36 (4.9)  5.42  EC(F3)3 2.87 2.75 EC(F1)1  4.76  4.80  EC(F4)4 1.67 1.39 EC(F2)2  3.61  3.50  VC0 4.71 4.69 EC(F3)3  2.87  2.75  VC(F1)1 3.70 3.59 EC(F4)4  1.67  1.39  VC(F2)2 2.61 2.40 VC0  4.71  4.69  VC(F1)1  3.70  3.59  VC(F2)2  In EC, the C–C bond is composed of 2.61  sp3 -sp3 carbons. On the2.40  other hand, the C–C bond in VC is

the sp2 -sp2 carbons. Hence, the C=C bond length of VC is significantly shorter than that of EC. The 3‐sp3 carbons. On the other hand, the C–C bond in VC is  In EC, the C–C bond is composed of sp C=C bond lengths of VC, VC(F1) and VC(F2) were 1.326, 1.323 and 1.321 Å, respectively. the sp2‐sp2 carbons. Hence, the C=C bond length of VC is significantly shorter than that of EC. The  C=C bond lengths of VC, VC(F1) and VC(F2) were 1.326, 1.323 and 1.321 Å, respectively.  3.2. Structures of Li+ -EC and Li+ -VC Complexes

+‐EC and Li+‐VC Complexes  3.2. Structures of Li The optimized structures of Li+ -EC and -VC complexes are given in Figure 3. The procedure of optimization was carried out without +symmetry restriction. The Li+ ion binds to the oxygen atom of The optimized structures of Li ‐EC and ‐VC complexes are given in Figure 3. The procedure of  + ion binds to the oxygen atom of  the C=O carbonyl of EC. The moiety of Li–O–C is almost linear in the Li+ EC complexes. The bond optimization was carried out without symmetry restriction. The Li + EC complexes. The bond  the C=O carbonyl of EC. The moiety of Li–O–C is almost linear in the Li distances of Li–O and C=O carbonyl were calculated to be 1.740 and 1.225 Å, respectively. The C–O distances of Li–O and C=O carbonyl were calculated to be 1.740 and 1.225 Å, respectively. The C–O  bond in the C=O carbonyl is elongated by the interaction with the Li+ ion. In the Li+ -VC system, + ion. In the Li+‐VC system, similar  similar bond in the C=O carbonyl is elongated by the interaction with the Li binding structures were obtained.

binding structures were obtained. 

 

Figure Figure 3. Optimized structures of lithium‐EC and ‐VC complexes calculated at the B3LYP/6‐31G(d)  3. Optimized structures of lithium-EC and -VC complexes calculated at the B3LYP/6-31G(d) level. Bond lengths are in Å.  level. Bond lengths are in Å. The distances of Li+ from the oxygen atom of C=O carbonyl are plotted in Figure 4 as a function  of the number of F atoms. The Li+‐EC complex has a distance of 1.740 Å. The distance increases with  increasing n. In n = 4, the distance is as large as 1.779 Å. 

Atoms 2016, 4, 4

5 of 8

The distances of Li+ from the oxygen atom of C=O carbonyl are plotted in Figure 4 as a function of the number of F atoms. The Li+ -EC complex has a distance of 1.740 Å. The distance increases with increasing n. In n = 4, the distance is as large as 1.779 Å. Atoms 2016, 4, 4 

3.3. Binding Structures of

5 of 10 

Li+ -EC

Complexes on Graphene

3.3. Binding Structures of Li+‐EC Complexes on Graphene  Atoms 2016, 4, 4 

5 of 10 

The optimized structures of the complexes of EC(F) and EC(F4) with graphene are illustrated The optimized structures of the complexes of EC(F) and EC(F4) with graphene are illustrated in  3.3. Binding Structures of Li+‐EC Complexes on Graphene  in FigureFigure  5. In 5.  both complexes, thethe  C=O carbonyl towardthe  thecenter  center benzene In  both  complexes,  C=O  carbonyl orients orients  toward  of  of the the benzene  ring ring of  of The optimized structures of the complexes of EC(F) and EC(F4) with graphene are illustrated in  graphene. The bond lengths of the C=O carbonyl of EC and the height of the oxygen atom (h) of EC graphene. The bond lengths of the C=O carbonyl of EC and the height of the oxygen atom (h) of EC  Figure  5.  In  both  complexes,  the  C=O  carbonyl  orients  toward  the  center  of  the  benzene  ring  of  are 1.783 and 1.915 Å, respectively. In the case of EC(F4), the height becomes lower (h = 1.882 Å).  are 1.783 and 1.915 Å, respectively. In the case of EC(F4), the height becomes lower (h = 1.882 Å). graphene. The bond lengths of the C=O carbonyl of EC and the height of the oxygen atom (h) of EC  are 1.783 and 1.915 Å, respectively. In the case of EC(F4), the height becomes lower (h = 1.882 Å). 

    Figure 4. Bond lengths of the Li ‐oxygen atom of Li‐EC complexes plotted as a function of the number  + +

Figure 4. Bond lengths of the Li -oxygen atom of Li-EC complexes plotted as a function of the number +‐oxygen atom of Li‐EC complexes plotted as a function of the number  of F‐atoms. The values were calculated at the B3LYP/6‐31G(d) level.  Figure 4. Bond lengths of the Li of F-atoms. The values were calculated at the B3LYP/6-31G(d) level. of F‐atoms. The values were calculated at the B3LYP/6‐31G(d) level. 

+‐graphene ternary systems. Bond lengths are in Å. (A); Li +‐ +‐graphene ternary systems. Bond lengths are in Å. (A); Li +‐ Figure 5. Optimized structures of EC‐Li Figure 5. Optimized structures of EC‐Li

+ -graphene ternary systems. Bond lengths are in Å. (A); Li` -EC Figure 5.EC Complexeson Graphene, (B); Li Optimized structures of EC-Li+‐EC(F4) Complexeson Graphene.  EC Complexeson Graphene, (B); Li+‐EC(F4) Complexeson Graphene.  ` Complexes on Graphene, (B); Li -EC(F4) Complexes on Graphene.

Atoms 2016, 4, 4 Atoms 2016, 4, 4 

6 of 8 6 of 10 

Figure  6  shows  the  distance  the  oxygen  of  EC  the  graphene  as  a  Figure 6 shows the distance of theof  oxygen of EC from thefrom  graphene surface (h)surface  plotted(h)  as aplotted  function +‐VC are also  function of n. The distance decreases slightly from 1.915 to 1.883 Å. The values of Li of n. The distance decreases slightly from 1.915 to 1.883 Å. The values of Li+ -VC are also plotted as a plotted as a dashed line in Figure 6: 1.907 Å (n = 0), 1.899 Å (n = 1) and 1.893 Å (n = 2).  dashed line in Figure 6: 1.907 Å (n = 0), 1.899 Å (n = 1) and 1.893 Å (n = 2).

+‐EC (or Li+‐VC) from the graphene surface (h) plotted as  Figure 6. Distances of the lithium atom of Li Figure 6. Distances of the lithium atom of Li+ -EC (or Li+ -VC) from the graphene surface (h) plotted a function of n, the number of F‐atoms. The values were calculated at the B3LYP/6‐31G(d) level. (a)  as a function of n, the number of F-atoms. The values were calculated at the B3LYP/6-31G(d) level. + + + -VC. (a) Li Li+‐EC; (b) Li -EC; (b) Li‐VC. 

3.4. Binding Energies of Li+‐EC Complexes on Graphene  3.4. Binding Energies of Li+ -EC Complexes on Graphene The binding energies of EC and VC in several binding systems are calculated;these systems are  The binding energies of EC and VC in several binding systems are calculated;these systems are expressed by:  expressed by:  → EC(F)‐( Li EC(F) + Li ECpFq ` Li+ +Ñ ECpFq-pLi+ q+)(a)  paq + + +)‐graphene (b)  ECpFq-pLi q` graphene Ñ ECpFq-pLi+ q-graphene pbq ) + graphene → EC(F)‐( Li EC(F)‐( Li

VCpFq ` Li++Ñ VCpFq-pLi++q) (c)  pcq VC(F) + Li  → VC(F)‐( Li + VCpFq-pLi q +` graphene Ñ VCpFq-pLi+ q-graphene pdq +)‐graphene (d)  VC(F)‐( Li ) +graphene → VC(F)‐( Li

Reaction (a) means the binding reaction of lithium ion to solvent EC(F). The solvated lithium ion Reaction (a) means the binding reaction of lithium ion to solvent EC(F). The solvated lithium ion  is formed by Reaction (a). Reaction (b) is the binding reaction of the solvated lithium ion EC(F)-(Li+ ) to +)  is formed by Reaction (a). Reaction (b) is the binding reaction of the solvated lithium ion EC(F)‐(Li theto the graphene surface. Reactions (c) and (d) in VC(F) are the same as EC(F).  graphene surface. Reactions (c) and (d) in VC(F) are the same as EC(F). The reverse reaction of Reaction (a), i.e., desolvation of Li+ from+ from solvent EC(F), Li solvent EC(F), Li+ -EC Ñ+‐EC → Li Li+ + EC,+ +  The reverse reaction of Reaction (a), i.e., desolvation of Li is the important process process  in the movement of the Li+of  ion from tosolvent  graphene. The weaker EC, most is  the  most  important  in  the  movement  the  Li+  solvent ion  from  to  graphene.  The  + + ion by EC is a positive factor as a solvent. Namely, the weak solvation of  solvation of the Li ion by EC is a positive factor as a solvent. Namely, the weak solvation of Li+ is weaker solvation of the Li preferred as an effective solvent. Li+ is preferred as an effective solvent.  The calculated binding energies generated by Reactions (a) to (d) are plotted in Figure 7. The The calculated binding energies generated by Reactions (a) to (d) are plotted in Figure 7. The  binding energies of Li+ to+ to EC(F) (Reaction (a)) decreases with increasing n: 56.1 kcal/mol (n = 0) and  EC(F) (Reaction (a)) decreases with increasing n: 56.1 kcal/mol (n = 0) and binding energies of Li   39.6 kcal/mol (n = 4). The results suggest that EC(F4) is more preferred as the solvent than EC(H), 39.6 kcal/mol (n = 4). The results suggest that EC(F4) is more preferred as the solvent than EC(H),  + ion is the positive factor as the Li + solvent.  because the weaker solvation of the Li+ ion is the positive factor as the Li+ solvent. because the weaker solvation of the Li Furthermore, the binding energies for Reactions (b), (c) and (d) show the same tendency as Reaction (a).

Atoms 2016, 4, 4 Atoms 2016, 4, 4 

7 of 8 7 of 10 

Figure 7. Binding energies in several systems plotted as a function of n, the number of F‐atoms. (a)  Figure 7. Binding energies in several systems plotted as a function of n, the number of F-atoms. +‐graphene; (c) VC‐Li+; and (d) VC‐Li+‐graphene.  EC‐Li+; (b) EC‐Li (a) EC-Li+ ; (b) EC-Li+ -graphene; (c) VC-Li+ ; and (d) VC-Li+ -graphene.

Furthermore,  the  binding  energies  for  Reactions  (b),  (c)  and  (d)  show  the  same  tendency  as  4. Discussion Reaction (a).  The weaker solvation of Li+ ion is a positive factor as the solvent. In the present study, the effects of4. Discussion  F-substitution on the binding energies of EC to the Li+ ion and Li+ (EC) to the graphene surface have been investigated by means of the DFT method. The calculations showed that the binding energy The weaker solvation of Li+ ion is a positive factor as the solvent. In the present study, the effects  decreases with increasing n. F-substitution makes the interaction of graphene with EC(F) significantly of F‐substitution on the binding energies of EC to the Li+ ion and Li+(EC) to the graphene surface have  weaker due to the small dipole moment of EC(F) (1.4 D vs. 5.4 D). This result implies that desolvation been investigated by means of the DFT method. The calculations showed that the binding energy  of the Li+ ion from the F-substituted system, EC(F)-Li+ , occurs more easily than that of EC(H). Hence, decreases with increasing n. F‐substitution makes the interaction of graphene with EC(F) significantly  the ability of EC(F) is greater than that of EC(H). weaker due to the small dipole moment of EC(F) (1.4 D vs. 5.4 D). This result implies that desolvation  In the present study, a carbon sheet composed of 14 benzene rings (C42 H16 ) was used as a model of of the Li+ ion from the F‐substituted system, EC(F)‐Li+, occurs more easily than that of EC(H). Hence,  the graphene surface. This is a small number of rings for modeling a finite graphene surface. However, the ability of EC(F) is greater than that of EC(H).  recent theoretical calculations indicated that this size of graphene sheet can reasonably represent the In the present study, a carbon sheet composed of 14 benzene rings (C42H16) was used as a model  electronic states of the graphene system [18–21]. Therefore, the present calculations would provide of  the  graphene  surface.  This  is  a  small  number  of  rings  for  modeling  a  finite  graphene  surface.  useful information on the ternary interaction system. However,  recent  theoretical  calculations  indicated  that  this  size  of  graphene  sheet  can  reasonably  represent the electronic states of the graphene system [18–21]. Therefore, the present calculations would  Acknowledgments: The author acknowledges partial support from JSPS KAKENHI Grant Number 15K05371 provide useful information on the ternary interaction system.  and MEXT KAKENHI Grant Number 25108004. Author Contributions: Mami Mutoh, Teruo Kusaka and Mariko Nakamura carried out the DFT calculations of Acknowledgments: The author acknowledges partial support from JSPS KAKENHI Grant Number  Li` -EC Complexes on Graphene. Yasuhiro Yoshida, Junichiro Iida, and Hiroto Tachikawa conceived of the study, 15K05371 and MEXT KAKENHI Grant Number 25108004.  and participated in its design and coordination and helped to draft the manuscript.

Author  Mami  Mutoh,  Teruo ofKusaka  Conflicts ofContributions:  Interest: The authors declare no conflict interest.and  Mariko  Nakamura  carried  out  the  DFT  calculations  of  Li+‐EC  Complexeson  Graphene.  Yasuhiro  Yoshida,  Junichiro  Iida,  and  Hiroto  References Tachikawa  conceived  of  the  study,  and  participated  in  its  design  and  coordination  and  helped  to  1. draft the manuscript.  Bhatta, M.D.; Cho, M.; Cho, K. Interaction of Li+ ions with ethylene carbonate (EC): Density functional theory calculations. Appl. Surf. Sci. 2010, 257, 1463–1468. [CrossRef] Conflicts of Interest: The authors declare no conflict of interest.  2. Nakagawa, H.; Fujino, Y.; Kozono, S.; Katayama, Y.; Nukuda, T.; Sakaebe, H.; Matsumoto, H.; Tatsumi, K. Application of nonflammable electrolyte with room temperature ionic liquids (RTILs) for lithium-ion cells. References  J. Power Sources 2007, 174, 1021–1026. [CrossRef] 1. Bhatta,  M.D.;  Cho,  M.;  Cho,  K.  Interaction  of  Li+  ions  with  ethylene  carbonate  (EC):  Density  functional  Xing, L.; Li, W.; Wang, C.; Gu, F.; Xu, M.; Tan, C.; Yi, J. Theoretical Investigations on Oxidative Stability 3. theory calculations. Appl. Surf. Sci.2010, 257, 1463–1468.  of Solvents and Oxidative Decomposition Mechanism of Ethylene Carbonate for Lithium Ion Battery Use. 2. Nakagawa, H.; Fujino, Y.; Kozono, S.; Katayama, Y.; Nukuda, T.; Sakaebe, H.; Matsumoto, H.; Tatsumi, K.  J. Phys. Chem. B 2009, 113, 16596–16602. [CrossRef] [PubMed] Application of nonflammable electrolyte with room temperature ionic liquids (RTILs) for lithium‐ion cells.  J. Power Sources 2007, 174, 1021–1026.  3. Xing, L.; Li, W.; Wang, C.; Gu, F.; Xu, M.; Tan, C.; Yi, J. Theoretical Investigations on Oxidative Stability of  Solvents and Oxidative Decomposition Mechanism of Ethylene Carbonate for Lithium Ion Battery Use. J.  Phys. Chem. B 2009, 113, 16596–16602. 

Atoms 2016, 4, 4

4.

5.

6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16.

17. 18. 19. 20.

21.

8 of 8

Zhang, H.L.; Sun, C.H.; Li, F.; Liu, C.; Tan, J.H.; Cheng, M. New insight into the interaction between propylene carbonate-based electrolytes and graphite anode material for lithium ion batteries. J. Phys. Chem. C 2007, 111, 4740–4748. [CrossRef] Zhu, Y.; Casselman, M.D.; Li, Y.; Wei, A.; Abrahama, D.P. Perfluoroalkyl-substituted ethylene carbonates: Novel electrolyte additives for high-voltage lithium-ion batteries. J. Power Sources 2014, 246, 184–191. [CrossRef] Werber, T.; Walther, J.H.; Jaffe, R.L.; Halicioglu, T.; Koumoutsakos, P. On the water-carbon interaction for use in molecular dynamics simulations of graphite and carbon nanotubes. J. Phys. Chem. B 2003, 107, 1345–1352. Rana, M.; Chandra, A. Filled and empty states of carbon nanotubes in water: Dependence on nanotube diameter, wall thickness and dispersion interactions. J. Chem. Sci. 2007, 119, 367–376. [CrossRef] Castriota, M.; Cazzanelli, E.; Nicotera, I.; Coppola, L.; Oliviero, C.; Ranieri, G.A. Temperature dependence of lithium ion solvation in ethylene carbonate-LiClO4 solutions. J. Chem. Phys. 2003, 118, 5537–5541. [CrossRef] Tachikawa, H. A direct molecular orbital-molecular dynamics study on the diffusion of the Li ion on a fluorinated graphene surface. J. Phys. Chem. C 2008, 112, 10193–10199. [CrossRef] Abe, S.; Nagoya, Y.; Watari, F.; Tachikawa, H. Interaction of Water Molecules with Graphene: A Density Functional Theory and Molecular Dynamics Study. Jpn. J. Apl. Phys. 2010, 49, AH071–AH074. [CrossRef] Abe, S.; Nagoya, Y.; Watari, F.; Tachikawa, H. Structures and Electronic States of Water Molecules on Graphene Surface: A Density Functional Theory Study. Jpn. J. Apl. Phys. 2010, 49, GJ131–GJ133. [CrossRef] Amft, M.; Sanyal, B.; Eriksson, O.; Skorodumova, N.V. Small gold clusters on graphene, their mobility and clustering: A DFT study. J. Phys. Condens. Matter 2011, 23, 205301. [CrossRef] [PubMed] Ding, N.; Lu, X.Q.; Wu, C.M.L. Nitrated tyrosine adsorption on metal-doped graphene: A DFT study. Comput. Mater. Sci. 2012, 51, 141–145. [CrossRef] Lim, D.H.; Wilcox, J. DFT-Based Study on Oxygen Adsorption on Defective Graphene-Supported Pt Nanoparticles. J. Phys. Chem. 2011, 115, 22742–22747. [CrossRef] Yanai, T.; Tew, D.; Handy, N. A new hybrid exchange-correlation functional using the Coulomb-attenuating method (CAM-B3LYP). Chem. Phys. Lett. 2004, 393, 51–57. [CrossRef] Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Montgomery, J.A., Jr.; Vreven, T.; Kudin, K.N.; Burant, J.C.; et al. Gaussian 03, Revision B.04; Gaussian, Inc.: Pittsburgh, PA, USA, 2003. Payne, R.; Theodorou, I.E. Dielectric properties and relaxation in ethylene carbonate and propylene carbonate. J. Phys. Chem. 1972, 76, 2892–2900. [CrossRef] Mirzaei, M.; Ahangari, R.S. Formations of CNT modified 5-(halogen)uracil hybrids: DFT studies. Superlattice Microstruct. 2014, 65, 375–379. [CrossRef] Peyghan, A.A.; Noei, M.; Yourdkhani, S. Al-doped graphene-like BN nanosheet as a sensor for para-nitrophenol: DFT study. Superlattice Microstruct. 2013, 59, 115–122. [CrossRef] Jenness, G.R.; Karalti, O.; Jordan, K.D. Benchmark calculations of water-acene interaction energies: Extrapolation to the water-graphene limit and assessment of dispersion-corrected DFT methods. Phys. Chem. Chem. Phys. 2010, 12, 6375–6381. [CrossRef] [PubMed] Pisula, W.; Feng, X.; Mullen, K. Charge-Carrier Transporting Graphene-Type Molecules. Chem. Mater. 2011, 23, 554–567. [CrossRef] © 2015 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons by Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).