Ionic

Jul 2, 2018 - Grupo Mesturas, Departamento de Física e Ciencias da Terra, Facultade de Ciencias, Universidade de A Coruña, Campus A. Zapateira s/n ...
0 downloads 0 Views 3MB Size
Article Cite This: J. Phys. Chem. C XXXX, XXX, XXX−XXX

pubs.acs.org/JPCC

Density Functional Study of Charge Transfer at the Graphene/Ionic Liquid Interface ́ ez-Gonzaĺ ez,† A. García-Fuente,‡ A. Vega,§ J. Carrete,∥ O. Cabeza,⊥ L. J. Gallego,† V. Gom and L. M. Varela*,†

Downloaded via TUFTS UNIV on July 4, 2018 at 05:12:47 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.



Grupo de Nanomateriais, Fotónica e Materia Branda, Departamento de Física de Partículas, Universidade de Santiago de Compostela, Campus Vida s/n, E-15782 Santiago de Compostela, Spain ‡ Departamento de Física, Universidad de Oviedo, E-33007 Oviedo, Spain § Departamento de Física Teórica, Atómica y Ó ptica, Universidad de Valladolid, E-47011 Valladolid, Spain ∥ Institute of Materials Chemistry, TU Wien, A-1060 Vienna, Austria ⊥ Grupo Mesturas, Departamento de Física e Ciencias da Terra, Facultade de Ciencias, Universidade de A Coruña, Campus A Zapateira s/n, E-15012 A Coruña, Spain ABSTRACT: We use density functional theory to analyze the charge transfer between lithium or magnesium cations and a graphene wall beyond the predictions of classical Marcus theory. To that end, metal atoms are placed in three different kinds of environments: (i) in a vacuum, (ii) among fluorine atoms to simulate a molten salt, and (iii) in an ionic liquid. We prove that a complete charge transfer to the electrodes takes place in all of the studied environments and that the charge transfer process starts at longer distances from the electrode in the ionic liquid. Vertical ionization potentials and vertical electron affinities are studied, and they confirm that the nanoconfined region close to the electrode is a favorable environment for electronic exchange. No significant difference between monovalent and divalent cations was found. Our results suggest a certain catalyzing effect of ionic liquids regarding metal-electrode charge transfer in these densely ionic environments. Moreover, they show that ionic liquids can actually enhance charge transfer to electrodes in electrochemical devices without significantly altering the nature of the process.



aprotic ILs in the bulk or near an interface33−39 have been reported, as well as different DFT studies, both in bulk and at interfaces.35,40−42 From these studies, as well as from others involving molecular cosolvents and other solutes (see ref 43 for a recent review on solvation in ILs), a detailed picture of a peculiar way of solvating metals in these media has emerged: metal cations are swollen into the polar nanoregions of these densely ionic solvents, forming a solid-like pseudolattice of metal-anion solvates or solvation complexes (according to the nanostructural solvation paradigm43), limiting the electrical conductivity and the ability of metal cations to reach the electrochemical interface. The formation of these metal-anion carrying species, quite different from those found in traditional solutions of metal cations in molecular solvents, is expected to have also important consequences in electron transfer in this highly ionic environments. A correct microscopic description of electron transfer in ILs is crucial to their proper use in electrochemical applications.

INTRODUCTION Over the past 2 decades, the field of ionic liquids (ILs) has enormously expanded due to the current and expected applications of these green, amphiphilically nanostructured “designer” solvents,1,2 whose properties can be specifically “tuned” for each application choosing the appropriate building blocks from a large number of known cations and anions.3 One of the most outstanding applications of ILs is as electrolytes in electrochemical devices, specifically in batteries, fuel cells, and supercapacitors, as has been repeatedly reported in recent years.4−7 This usage relies on the very large electrochemical stability of these solvents, which results in unrivaled electrochemical windows. Notwithstanding, due to this same property, ILs must be doped with some electrochemically active additive, usually metal salts, in order to get practical electrolytes for batteries. In the past few years, different studies on this subject have been reported for mixtures of ILs with alkaline8−29 and alkaline earth30−32 metals, including classical molecular dynamics (MD) and density functional theory (DFT) to dissect the microscopic mechanisms behind solvation, structure, and transport in these media. In particular, several MD analyses of mixtures of monovalent and divalent salts with protic and © XXXX American Chemical Society

Received: March 23, 2018 Revised: June 8, 2018

A

DOI: 10.1021/acs.jpcc.8b02795 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

paper. Their possible contributions are estimated to be well within the uncertainty of the model70 and significantly lower than, for instance, solvent effects. Therefore, the added computational cost would not be justified. The systems are placed in a box of X and Y dimensions matching the simulated graphene sheet (comprising 32 C atoms in simulations with IL and 72 C atoms in all others and assuming its planarity in order to focus in the charge transfer and avoid extra complexities) and Z dimension of 20 Å (25 Å in the case of the systems with IL), in order to have sufficient vacuum so as to avoid influence from periodic replicas. Monkhorst−Pack grids of 8 × 8 × 1, 2 × 2 × 1, and 6 × 6 × 1 were employed for the calculations involving boxes in vacuum, with fluorine, and with IL respectively. We base our initial configuration for each system on the results of MD simulations with large numbers of particles, which can capture the effect of medium- and long-range order: the initial simulation box for the system with IL was built from previously reported MD results of [BMIM][BF4] near a graphene wall,36 taking the average configuration of the IL pair when placed close to a neutral graphene sheet. We then place the salt ion at a predefined distance from the wall and run a self-consistent-field (SCF) DFT calculation to obtain converged wave functions for the Kohn−Sham electrons and the charge density. We sample at least 35 distances along the normal to the graphene surface (see Figure 1) for each studied configuration. The charges of all atoms were calculated from the charge density using the Bader criterion72−75 after each simulation. In the Bader method, space is partitioned in volumes divided by zero-flux surfaces, which guarantees that the sum of the Bader charges equals the total charge of the system. The geometries for each of the three different systems are shown in Figure 1. With respect to the configuration with only a graphene sheet plus a Li or Mg atom, the system with fluorine adds a F atom at a distance of 4 Å from the wall and at 1.77 Å from the line perpendicular to the wall along which the central atom is moved. Those at 1.77 Å result from an average of the distance between Li and F atoms in LiF and the sum of the ionic radii of Li+ and F−. When we simulate Mg instead of Li, we introduce an additional F atom in the system in order to preserve electroneutrality. The position of this new atom is symmetric with respect to the line where the magnesium atom is getting moved. Finally, to simulate the system with IL, a [BMIM]+ ion and a [BF4]− ion are placed at around 3.5 Å from the carbon sheet, qualitatively reproducing the structure of the IL near a carbon interface, which was calculated previously by means of classical MD simulations. The lateral distance between them also matches the value found in MD simulations, which is longer than the radii of the metal atoms.36,39 The orientation of the IL pair relative to the graphene surface is the average orientation/configuration extracted from MD simulations of [BMIM][BF4] near a graphene wall, given the low mobility of the species in the innermost layer of the electrical double layer in these systems. More details of the system with IL are shown in Figure 2.

Up to now, most of the theoretical efforts devoted to the description of this phenomenon have been done using MD simulations and in the framework of the classical Marcus theory,44 including some extensions to account for deviations from linear behavior,45,46 which focuses in the description of the influence of solvent fluctuations on the rate of electron transfer and whose validity in bulk ILs was originally proved by LyndenBell.47,48 The predictions of Marcus theory may show significant deviations if fluctuations around reactants and products are similar, as has been shown by MD studies based on DFT.46 These deviations associated with the existence of different solvation states for a redox species were further confirmed by Salanne and co-workers49 by means of classical MD simulations of the Fe2+/Fe3+ electron transfer reaction in ILs under nanoconfinement, in which the authors found two solvation states for Fe3+. Although classical simulations provide a valuable average picture of this phenomenon, a proper microscopic analysis requires a detailed description of the electronic densities of the reacting species throughout the process, so quantum chemical methods are needed. Several examples of studies on the subject of charge transfer using DFT can be found in the literature: variations of effective charge of a recoiling atom along crystallographic directions of a SiC interface,50 analysis of the net charge of Li and Si atoms in LixSi compounds,51 use of DFTcalculated charges to assess a theoretical transfer parameter,52 charge analysis in an electron donor−acceptor complex of PA and CA,53 and even a study of charge transfer between NO2 and graphene,54 among others. With respect to ILs in particular, the study of charge transfer using DFT has only been applied to a pure IL and a Li electrode by Yildirim et al.55 Other DFT studies of ILs concerning their interaction with Li metal have also been reported in the past decade, regarding their adsorption,56−58 electrochemical stabilities,59−61 and the chemical decomposition of the IL species.62−65 In this work we present, for the first time up to our knowledge, a DFT study of the electron transfer from a monovalent or a divalent salt cation (lithium and magnesium, respectively) to a graphene wall as a function of the distance of this atom to the interface. In order to mimic different environments, this analysis has been performed for three different situations: in vacuum, with added fluorine atoms to mimic a high-temperature molten salt, and in the presence of 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF4]), to study the transfer reaction in ILs. In all cases, highly electronegative fluorine and spherically symmetric fluorine-based anions (BF4−) have been chosen, since F will expectedly enhance the charge transfer. On the other hand, we use an imidazolium cation as a model IL, since this family is frequently chosen in electrochemical applications, and it has been widely studied and properly parametrized for its computational analysis. The structure of the paper is as follows. After this Introduction, the computational details are included, and in section 3, we report our results and discussion. Finally, in section 4, the conclusions of the paper are summarized.



RESULTS AND DISCUSSION Our simulations aim at qualitatively understanding the trends of electronic redistribution depending on the chemical environment that the salt cation may encounter, as well as how this environment influences the net charge transferred to the graphene sheet. Although the system is not fully relaxed in our vertical treatment of the process (the environment is not relaxed upon the movement of the salt cation through the selected path), the fact that the atomic and chemical environment and



COMPUTATIONAL METHOD For all DFT calculations in this paper, we have used the projector augmented wave (PAW) method under the Perdew− Burke−Ernzerhoff (PBE) form of the generalized gradient approximation (GGA) for exchange and correlation as implemented in the Vienna ab initio simulation package (VASP).66−69 Dispersion corrections are not used in this B

DOI: 10.1021/acs.jpcc.8b02795 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

focus to the electronic configuration only when it matters the most. The results of these calculations for the system containing fluorine are shown in Figure 3, and the ones for the IL system are represented in Figure 4. The first relevant result to analyze is what happens in vacuum (only graphene and salt cation; dashed lines in Figures 3 and 4). Obviously, in this case all electronic charge lost by the atom is gained by the wall. Both the Li and Mg atoms are essentially neutral far away from the graphene sheet, but they start to ionize as they approach its surface, specifically at ∼3−4 Å, until they reach the expected ionization charge. However, while the Mg atom does not start transferring its charge until a distance of ∼4 Å, the Li atom is not entirely neutral even at long distances. This effect correlates with the lower electronegativity of lithium with respect to magnesium (0.98 and 1.31 in the Pauling scale, respectively), which makes electron transfer easier for the former. Another remarkable feature is that the Mg atom loses its two valence electrons in two separate steps, as the plateau around 2 Å suggests. The presence of fluorine atoms in the system drastically alters the shape of the charge curve. The fluorine atom remains ionized almost independently of the salt cation position (an unsurprising effect in view of its high electronegativity: 3.98 in the Pauling scale). There seem to be two distinct zones in the system with Li with regard to charge transfer and three in the system with Mg. In the Li−F system, the Li atom starts out ionized (with around +0.3e) when it is far away from graphene, and increases its charge exponentially until it reaches almost complete ionization at 4 Å, where the fluorine atom is placed. This charge is mainly transferred to the graphene wall, although some ends in the fluorine atom. From this point on, the charges of all species in the simulated system remain almost constant, as they reach an equilibrium point where Li+ is bonded with F−. In the Mg−F2 system, the behavior is similar to that in the monovalent case: fluorine atoms maintain their ionized state throughout all distances, and the salt cation increases its charge when approaching graphene, mainly by electron transfer to the wall but also to the fluorine atoms to a lesser extent. The main difference seems to be that Mg is almost fully ionized at a longer distance from the wall than Li. This effect is probably related to the fact that in this situation two fluorine atoms were introduced instead of one, so the opportunities to shed electrons are increased. Again, at distances to the wall shorter than 4 Å, the salt atom seems to be bonded with the F−, so the charges of all three species in the system remain approximately constant. However, when the salt cation is very close to the wall (≲2 Å), charge is transferred from the fluorine atoms to graphene. The preponderance of kinetic energy in the regime of high charge densities could be favoring a more uniform distribution of charge between graphene and the fluorine atom. However, it is clearly shown that the high electronegativity of the fluorine atom does not allow the graphene to reach negative enough charges (only around −0.1e/−0.25e in the monovalent/ divalent system), so although the salt cation atom is able to start losing electrons at longer distances from the wall than in vacuum, the transfer efficiency to graphene is much lower. This efficiency is obviously very important with regard to electrochemical applications. In order to get some qualitative insight into the energy barriers our salt cation must overcome when getting close to the graphene wall, we show in Figure 5 the total energy of the systems with magnesium as a function of the magnesium distance to the graphene sheet, choosing as the origin of energy

Figure 1. Starting configurations for all simulated systems: (a) graphene + Li or Mg cation system, (b) system with fluorine, and (c) system with IL, including the periodic replicas. The dashed red line represents the normal to the graphene surface, along which Li or Mg atoms are moved in order to get the starting geometry of the simulation for each distance to the wall. This line goes through the center of a hexagon of the graphene sheet (hollow site), which is the energetically favored site for the adsorption of groups I−III metal adatoms on graphene (see ref 71). Atom colors are the following: C, cyan; Li/Mg, pink; F, blue; N, purple; H, white; B, green.

interatomic distances match the optimal configuration given by previous MD simulations,36 whose results were benchmarked against experimental data, means that those environments will not differ much from snapshots of the ones occurring during the dynamical process. Of course, ab initio MD simulations could provide a more quantitative description of the systems investigated, but their complex nature makes the multiscale workflow combining classical MD and DFT a trade-off between accuracy and feasibility. Moreover, it is unrealistic to try to use DFT to simulate large boxes of solute for a long enough time to achieve finite-temperature equilibrated configurations close to the graphene interface with the same quality that MD can afford. We hence harness the strengths of each method by shifting our C

DOI: 10.1021/acs.jpcc.8b02795 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Figure 2. Simulated system of [BMIM][BF4] and Li/Mg and relevant distances between different species. Some replicas of the unit cell (which is shown in Figure 1c) are included for the sake of clarity. The blue arrow represents the direction along which the Li/Mg atom is being moved in order to get the starting geometry of the simulation for each distance to the graphene sheet.

Figure 3. Charges of the different species in the system with fluorine as a function of the (a) lithium or (b) magnesium distance to the graphene sheet. The dashed lines represent the charges of salt cation and graphene in the vacuum system. The red dotted line marks the position of the fluorine atom(s).

Figure 4. Charges of the different species in the system containing IL as a function of the (a) lithium or (b) magnesium distance to the graphene sheet. Dashed lines represent the charges of salt cation and graphene in the vacuum system. The positions of the other atoms of the system are represented in the background.

the configuration with the magnesium atom farthest from the graphene wall. The results are in good agreement with the

expected ones. In the system where only Mg and graphene are present, the energy of the system remains constant until the salt D

DOI: 10.1021/acs.jpcc.8b02795 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

2 Å. No qualitative differences between the monovalent and divalent salt cations are found in this regard. Furthermore, the efficiency of the charge transfer to the graphene wall is higher when IL is present than in vacuum or with fluorine atoms, without substantially altering the nature of its cations and anions (as the cation has the same charge at short and long distances of the salt cation to the wall and the anion maintains the same charge throughout the whole trajectory of the salt cation approaching the wall) in the whole process. Therefore, we could claim that ILs work as “catalyzers” of the charge transfer. We have also calculated the charges of all species in the system when the Mg atom is near (1.6 Å) and far away (8.0 Å) from the graphene sheet, when we add or remove an electron in the system. The results are shown in Figure 6. We can clearly see

Figure 5. Relative energies of the simulated systems with Mg, ΔE, as a function of the magnesium distance to the graphene sheet. Note that the origin of energy, ΔE = 0, has been chosen to be that of the configuration with the magnesium atom farthest from the graphene wall.

cation gets close to the wall. On the other hand, when we dope the system with two fluoride atoms, a minimum of energy is reached at 4 Å, where Mg forms bonded MgF2. However, when we have an IL added in the system, we recover the already known energy barriers that come from the anion and the cation, which were reported in previous works using classical MD simulations.39,76 These energy barriers are artificially high in these DFT calculations due to the lack of solvent relaxation. Finally, we analyze the electronic charge transfer in the IL system (shown in Figure 4). The main features of the behavior of the salt cation are qualitatively the same as in the previous fluorinated system, as there is a distance range far from the wall where the salt ion starts to ionize, and around 4.5 Å the charge stabilizes and remains approximately constant, mainly for Li. Again, the anion charge also remains constant independent of the position of the Li/Mg, as expected because of its nature as a heavily fluorinated compound. On the other hand, the IL cation is almost fully ionized when the salt atom is either very near or far away from the surface but displays a charge minimum when the Li/Mg atom is placed close to [BMIM]: the cation is half ionized when it is near the monovalent cation and almost fully neutral in the system with the divalent salt in the same position. This difference probably comes from the higher electronegativity of the divalent cation but also because electroneutrality has to be preserved, and the Mg atom can donate one more electron than the Li atom. If we had simulated this system with one extra anion, as we did in the system with fluorine, the cation charge curve should be approximately the same for both monovalent and divalent salt atoms. Moreover, it is clear that even when we approach bulk conditions (salt cation far away from the wall) the effect of the graphene wall over the IL anion is marked, as its charge is −0.98e while that of the IL cation is 0.87e, much more similar to the charge reported by Hollóczki et al. (approximately ±0.82) when the IL pair is completely isolated.35 In conclusion, clearly the main effect of the IL with respect to the vacuum system is to enhance the charge transfer from the corresponding salt cation: it starts at longer distances from the wall, reaching its fully ionized form at around 4.5 Å instead of at

Figure 6. Charges of the different species in the Mg + [BMIM][BF4] system as a function of the total charge of the system. The system with the Mg atom placed near the graphene sheet (1.6 Å) is represented with dots (a) while the system with this atom placed at 8 Å is represented with squares (b).

there is a difference in the behavior of the system’s components when the salt cation is near or far away from the wall, except for the anion: again, its charge remains practically immune to the addition or removal of one electron from the system, as it is fully ionized in all cases. An analogous behavior appears for the Mg atom when it is near the wall, as it is almost fully ionized and does not get influenced by a small variation of the number of electrons in the system, meaning that this geometric distribution forces the salt cation to give all its charge. At this point, it is natural to wonder where the excess/defect of charge in the system goes. It appears that both the IL cation and the wall distribute it among them, although not in the same proportion: carbon atoms in graphene seem to be more likely to lose some of their electrons when the total system charge is +1, while in the case of charge −1 the IL cation gets most of this electron excess. E

DOI: 10.1021/acs.jpcc.8b02795 J. Phys. Chem. C XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry C



On the other hand, when the salt cation remains far away from the graphene wall, the cation charge is almost constant, and only graphene and the salt cation change their charge in basically the same proportion. This is evidence that the IL pair is not affected by small charge variations in the system when it is in bulk phase (“isolated” from the salt cation), unlike in the situation where the salt cation is ionized near the IL pair, and hence the cation behaves as the predominant species that likes to lose/gain electrons. This was also observed in Figure 4, as in the distance range from ∼1.5 to 5 Å the IL cation can almost recover its electroneutrality because of the influence of the nearby salt cation. Finally, and in order to get some more insight into the ionic liquid system, we calculated the vertical electronic affinity (VEA) and the vertical ionization potential (VIP) of the system when the Mg atom is near and far away from the graphene sheet. These calculations are frequently done for clusters in DFT studies and are based on adding or removing an electron from the system without relaxing the structure.77,78 For clusters, these magnitudes are very similar to the corresponding ionization potentials and electronic affinities, so they tell us how difficult (energetically) it is to add or remove an electron from the system. In our system, we can qualitatively translate that to how difficult it is for an external potential (such as the one we should apply to a graphene electrode near our mixture of ionic liquid and salt) to be able to extract or supply electric energy (electrons) from or to the electrolyte. The results for the VIP and VEA for the system where the magnesium atom is close to the wall give values of 1.80 eV and −0.30 eV, respectively, and when this atom is placed far away from the graphene sheet, they change to 3.24 eV and −2.28 eV. This again shows that the system is more receptive to changes in its electronic environment when the metal ion is near the graphene wall, precisely what we want for the charge exchange to occur in a real battery.

Article

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

A. García-Fuente: 0000-0002-4570-8315 J. Carrete: 0000-0003-0971-1098 L. M. Varela: 0000-0002-0569-0042 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The financial support of the Spanish Ministry of Economy and Competitiveness (Projects MAT2014-57943-C3-1-P, MAT2014-57943-C3-3-P, MAT2017-89239-C2-1-P, MAT2017-89239-C2-2-P, FIS2014-59279-P) is gratefully acknowledged. Moreover, this work was funded by the Xunta de Galicia (Grants AGRUP2015/11 and GRC ED431C 2016/ 001). All these research projects were partially supported by FEDER. V.G.-G. thanks the Spanish Ministry of Education for his FPU grant. Facilities provided by the Galician Supercomputing Centre (CESGA) are also acknowledged. Funding from the European Union (COST Actions CM1206 and MP1303) is also acknowledged.



REFERENCES

(1) Rogers, R. D.; Seddon, K. R. Ionic liquids−solvents of the future? Science 2003, 302, 792−793. (2) Seddon, K. R. Ionic liquids: a taste of the future. Nat. Mater. 2003, 2, 363−365. (3) Plechkova, N. V.; Seddon, K. R. Methods and Reagents for Green Chemistry; John Wiley & Sons, Inc., 2007; pp 103−130. (4) MacFarlane, D. R.; Forsyth, M.; Howlett, P. C.; Pringle, J. M.; Sun, J.; Annat, G.; Neil, W.; Izgorodina, E. I. Ionic liquids in electrochemical devices and processes: managing interfacial electrochemistry. Acc. Chem. Res. 2007, 40, 1165−1173. (5) Armand, M.; Endres, F.; MacFarlane, D. R.; Ohno, H.; Scrosati, B. Ionic-liquid materials for the electrochemical challenges of the future. Nat. Mater. 2009, 8, 621−629. (6) Lu, X.; Burrell, G.; Separovic, F.; Zhao, C. Electrochemistry of room temperature protic ionic liquids: a critical assessment for use as electrolytes in electrochemical applications. J. Phys. Chem. B 2012, 116, 9160−9170. (7) Fedorov, M. V.; Kornyshev, A. A. Ionic liquids at electrified interfaces. Chem. Rev. 2014, 114, 2978−3036. (8) Sakaebe, H.; Matsumoto, H. N-Methyl-N-propylpiperidinium bis(trifluoromethanesulfonyl)imide (PP13-TFSI) - novel electrolyte base for Li battery. Electrochem. Commun. 2003, 5, 594−598. (9) García, B.; Lavallée, S.; Perron, G.; Michot, C.; Armand, M. Room temperature molten salts as lithium battery electrolyte. Electrochim. Acta 2004, 49, 4583−4588. (10) Henderson, W. A.; Passerini, S. Phase behavior of ionic liquidLiX mixtures: pyrrolidinium cations and TFSI-anions. Chem. Mater. 2004, 16, 2881−2885. (11) Diaw, M.; Chagnes, A.; Carré, B.; Willmann, P.; Lemordant, D. Mixed ionic liquid as electrolyte for lithium batteries. J. Power Sources 2005, 146, 682−684. (12) Zhou, Q.; Henderson, W. A.; Appetecchi, G. B.; Passerini, S. Phase behavior and thermal properties of ternary ionic liquid-lithium salt (IL-IL-LiX) electrolytes. J. Phys. Chem. C 2010, 114, 6201−6204. (13) Lee, S. Y.; Yong, H. H.; Lee, Y. J.; Kim, S. K.; Ahn, S. Two-cation competition in ionic-liquid-modified electrolytes for lithium ion batteries. J. Phys. Chem. B 2005, 109, 13663−13667. (14) Nicotera, I.; Oliviero, C.; Henderson, W. A.; Appetecchi, G. B.; Passerini, S. NMR investigation of ionic liquid-LiX mixtures:



CONCLUSIONS We report the first, to our knowledge, DFT simulations of the process of electronic charge transfer between an alkaline (monovalent) and alkaline earth (divalent) salt cation (Li, Mg) to graphene electrodes in different environments (vacuum, a molten salt, and an IL). A complete charge transfer to the electrodes was found to take place in all the studied systems, but the IL seems to lead to a significant enhancement in the charge transfer process, which starts at longer distances than in the other environments, ultimately leading to the same degree of charge transfer to the electrode. Higher efficiency in the charge transfer process and higher vertical ionization potentials were found for the configuration with the metal closer to the electrode, and the opposite trend was found for the vertical electron affinities, confirming the favorable electronic exchange in the nanoconfined region close to the electrode. Finally, no relevant difference between monovalent and divalent cations was found. All our results suggest a certain catalyzing effect of ILs regarding metal/electrode charge transfer in these densely ionic environments. Our findings are promising for the future electrochemical applications of ILs: enhanced charge transfer to the electrode without an associated change in their behavior as solvents make them more attractive candidates for such purposes. It is, hence, of utmost importance to extend this research to more kinds of species and to larger systems that include more details about the equilibrium structure of the bulk phase in contact with the electrode. F

DOI: 10.1021/acs.jpcc.8b02795 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C pyrrolidinium cations and TFSI-anions. J. Phys. Chem. B 2005, 109, 22814−22819. (15) Castriota, M.; Caruso, T.; Agostino, R. G.; Cazzanelli, E.; Henderson, W. A.; Passerini, S. Raman investigation of the ionic liquid n-methyl-n-propylpyrrolidinium bis(trifluoromethanesulfonyl)imide and its mixture with LiN(SO2CF3)2. J. Phys. Chem. A 2005, 109, 92−96. (16) Lassègues, J.; Grondin, J.; Talaga, D. Lithium solvation in bis(trifluoromethanesulfonyl)imide-based ionic liquids. Phys. Chem. Chem. Phys. 2006, 8, 5629−5632. (17) Markevich, E.; Baranchugov, V.; Aurbach, D. On the possibility of using ionic liquids as electrolyte solutions for rechargeable 5V Li ion batteries. Electrochem. Commun. 2006, 8, 1331−1334. (18) Xu, J.; Yang, J.; NuLi, Y.; Wang, J.; Zhang, Z. Additive-containing ionic liquid electrolytes for secondary lithium battery. J. Power Sources 2006, 160, 621−626. (19) Seki, S.; Kobayashi, Y.; Miyashiro, H.; Ohno, Y.; Usami, A.; Mita, Y.; Kihira, N.; Watanabe, M.; Terada, N. Lithium secondary batteries using modified-imidazolium room-temperature ionic liquid. J. Phys. Chem. B 2006, 110, 10228−10230. (20) Egashira, M.; Todo, H.; Yoshimoto, N.; Morita, M.; Yamaki, J. Functionalized imidazolium ionic liquids as electrolyte components of lithium batteries. J. Power Sources 2007, 174, 560−564. (21) Saito, Y.; Umecky, T.; Niwa, J.; Sakai, T.; Maeda, S. Existing condition and migration property of ions in lithium electrolytes with ionic liquid solvent. J. Phys. Chem. B 2007, 111, 11794−11802. (22) Borgel, V.; Markevich, E.; Aurbach, D.; Semrau, G.; Schmidt, M. On the application of ionic liquids for rechargeable Li batteries: high voltage systems. J. Power Sources 2009, 189, 331−336. (23) Lassègues, J.; Grondin, J.; Aupetit, C.; Johansson, P. Spectroscopic identification of the lithium ion transporting species in LiTFSI-doped ionic liquids. J. Phys. Chem. A 2009, 113, 305−314. (24) Nicolau, B. G.; Sturlaugson, A.; Fruchey, K.; Ribeiro, M. C. C.; Fayer, M. D. Room temperature ionic liquid-lithium salt mixtures: optical Kerr effect dynamical measurements. J. Phys. Chem. B 2010, 114, 8350−8356. (25) Zhou, Q.; Fitzgerald, K.; Boyle, P. D.; Henderson, W. A. Phase behavior and crystalline phases of ionic liquid-lithium salt mixtures with 1-alkyl-3-methylimidazolium salts. Chem. Mater. 2010, 22, 1203−1208. (26) Yoon, H.; Lane, G.; Shekibi, Y.; Howlett, P.; Forsyth, M.; Best, A.; MacFarlane, D. Lithium electrochemistry and cycling behaviour of ionic liquids using cyano based anions. Energy Environ. Sci. 2013, 6, 979−986. (27) Castiglione, F.; Famulari, A.; Raos, G.; Meille, S. V.; Mele, A.; Appetecchi, G. B.; Passerini, S. Pyrrolidinium-based ionic liquids doped with lithium salts: how does Li+ coordination affect its diffusivity? J. Phys. Chem. B 2014, 118, 13679−13688. (28) Aguilera, L.; Völkner, J.; Labrador, A.; Matic, A. The effect of lithium salt doping on the nanostructure of ionic liquids. Phys. Chem. Chem. Phys. 2015, 17, 27082−27087. (29) Lahiri, A.; Li, G.; Olschewski, M.; Endres, F. Influence of polar organic solvents in an ionic liquid containing lithium bis (fluorosulfonyl) amide: effect on the cation−anion interaction, lithium ion battery performance, and solid electrolyte interphase. ACS Appl. Mater. Interfaces 2016, 8, 34143−34150. (30) Cheek, G. T.; O’Grady, W. E.; El Abedin, S. Z.; Moustafa, E. M.; Endres, F. Studies on the electrodeposition of magnesium in ionic liquids. J. Electrochem. Soc. 2008, 155, D91−D95. (31) Giffin, G. A.; Moretti, A.; Jeong, S.; Passerini, S. Complex nature of ionic coordination in magnesium ionic liquid-based electrolytes: solvates with mobile Mg2+ cations. J. Phys. Chem. C 2014, 118, 9966− 9973. (32) Vardar, G.; Sleightholme, A. E. S.; Naruse, J.; Hiramatsu, H.; Siegel, D. J.; Monroe, C. W. Electrochemistry of magnesium electrolytes in ionic liquids for secondary batteries. ACS Appl. Mater. Interfaces 2014, 6, 18033−18039. (33) Méndez-Morales, T.; Carrete, J.; Bouzón-Capelo, S.; PérezRodríguez, M.; Cabeza, O.; Gallego, L. J.; Varela, L. M. MD simulations of the formation of stable clusters in mixtures of alkaline salts and

imidazolium-based ionic liquids. J. Phys. Chem. B 2013, 117, 3207− 3220. (34) Méndez-Morales, T.; Carrete, J.; Cabeza, O.; Russina, O.; Triolo, A.; Gallego, L. J.; Varela, L. M. Solvation of lithium salts in protic ionic liquids: a molecular dynamics study. J. Phys. Chem. B 2014, 118, 761− 770. (35) Hollóczki, O.; Malberg, F.; Welton, T.; Kirchner, B. On the origin of ionicity in ionic liquids. Ion pairing versus charge transfer. Phys. Chem. Chem. Phys. 2014, 16, 16880−16890. (36) Méndez-Morales, T.; Carrete, J.; Pérez-Rodríguez, M.; Cabeza, O.; Gallego, L. J.; Lynden-Bell, R. M.; Varela, L. M. Molecular dynamics simulations of the structure of the graphene-ionic liquid/alkali salt mixtures interface. Phys. Chem. Chem. Phys. 2014, 16, 13271−13278. (37) Gómez-González, V.; Docampo-Á lvarez, B.; Cabeza, O.; Fedorov, M. V.; Lynden-Bell, R. M.; Gallego, L. J.; Varela, L. M. Molecular dynamics simulations of the structure and single-particle dynamics of mixtures of divalent salts and ionic liquids. J. Chem. Phys. 2015, 143, 124507−124515. (38) Russina, O.; Caminiti, R.; Méndez-Morales, T.; Carrete, J.; Cabeza, O.; Gallego, L. J.; Varela, L. M.; Triolo, A. How does lithium nitrate dissolve in a protic ionic liquid? J. Mol. Liq. 2015, 205, 16−21. (39) Gómez-González, V.; Docampo-Á lvarez, B.; Méndez-Morales, T.; Cabeza, O.; Ivaništšev, V. B.; Fedorov, M. V.; Gallego, L. J.; Varela, L. M. Molecular dynamics simulation of the structure and interfacial free energy barriers of mixtures of ionic liquids and divalent salts near a graphene wall. Phys. Chem. Chem. Phys. 2017, 19, 846−853. (40) Angenendt, K.; Johansson, P. Ionic liquid based lithium battery electrolytes: charge carriers and interactions derived by density functional theory calculations. J. Phys. Chem. B 2011, 115, 7808−7813. (41) Haskins, J. B.; Bauschlicher, C. W., Jr; Lawson, J. W. Ab initio simulations and electronic structure of lithium-doped ionic liquids: structure, transport, and electrochemical stability. J. Phys. Chem. B 2015, 119, 14705−14719. (42) Gomez-Gonzalez, V.; Docampo-Alvarez, B.; Otero-Mato, J. M.; Cabeza, O.; Gallego, L. J.; Varela, L. M. Molecular dynamics simulations of the structure of mixtures of protic ionic liquids and monovalent and divalent salts at the electrochemical interface. Phys. Chem. Chem. Phys. 2018, 20, 12767−12776. (43) Varela, L.; Méndez-Morales, T.; Carrete, J.; Gómez-González, V.; Docampo-Á lvarez, B.; Gallego, L.; Cabeza, O.; Russina, O. Solvation of molecular cosolvents and inorganic salts in ionic liquids: a review of molecular dynamics simulations. J. Mol. Liq. 2015, 210, 178−188. (44) Marcus, R. A. On the theory of oxidation-reduction reactions involving electron transfer. I. J. Chem. Phys. 1956, 24, 966−978. (45) Small, D. W.; Matyushov, D. V.; Voth, G. A. The theory of electron transfer reactions: what may be missing? J. Am. Chem. Soc. 2003, 125, 7470−7478. + 2+ /Cuaq redox reaction exhibits strong (46) Blumberger, J. Cuaq nonlinear solvent response due to change in coordination number. J. Am. Chem. Soc. 2008, 130, 16065−16068. (47) Lynden-Bell, R. Can Marcus theory be applied to redox processes in ionic liquids? A comparative simulation study of dimethylimidazolium liquids and acetonitrile. J. Phys. Chem. B 2007, 111, 10800−10806. (48) Lynden-Bell, R. Does Marcus theory apply to redox processes in ionic liquids? A simulation study. Electrochem. Commun. 2007, 9, 1857− 1861. (49) Li, Z.; Jeanmairet, G.; Méndez-Morales, T.; Burbano, M.; Haefele, M.; Salanne, M. Confinement effects on an electron transfer reaction in nanoporous carbon electrodes. J. Phys. Chem. Lett. 2017, 8, 1925−1931. (50) Gao, F.; Xiao, H.; Zu, X.; Posselt, M.; Weber, W. J. Defectenhanced charge transfer by ion-solid interactions in SiC using largescale ab initio molecular dynamics simulations. Phys. Rev. Lett. 2009, 103, 027405. (51) Chevrier, V.; Zwanziger, J.; Dahn, J. First principles study of Li− Si crystalline phases: charge transfer, electronic structure, and lattice vibrations. J. Alloys Compd. 2010, 496, 25−36. (52) Kokalj, A. On the HSAB based estimate of charge transfer between adsorbates and metal surfaces. Chem. Phys. 2012, 393, 1−12. G

DOI: 10.1021/acs.jpcc.8b02795 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C (53) Shukla, M.; Srivastava, N.; Saha, S. Investigation of ground state charge transfer complex between paracetamol and p-chloranil through DFT and UV−visible studies. J. Mol. Struct. 2012, 1021, 153−157. (54) Leenaerts, O.; Partoens, B.; Peeters, F. Paramagnetic adsorbates on graphene: a charge transfer analysis. Appl. Phys. Lett. 2008, 92, 243125−243127. (55) Yildirim, H.; Haskins, J. B.; Bauschlicher, C. W., Jr; Lawson, J. W. Decomposition of ionic liquids at lithium interfaces. 1. Ab initio molecular dynamics simulations. J. Phys. Chem. C 2017, 121, 28214− 28234. (56) Valencia, H.; Kohyama, M.; Tanaka, S.; Matsumoto, H. Ab initio study of EMIM-BF4 crystal interaction with a Li (100) surface as a model for ionic liquid/Li interfaces in Li-ion batteries. J. Chem. Phys. 2009, 131, 244705−244715. (57) Buchner, F.; Forster-Tonigold, K.; Uhl, B.; Alwast, D.; Wagner, N.; Farkhondeh, H.; Groß, A.; Behm, R. J. Toward the microscopic identification of anions and cations at the ionic liquid/Ag(111) interface: a combined experimental and theoretical investigation. ACS Nano 2013, 7, 7773−7784. (58) Buchner, F.; Forster-Tonigold, K.; Bozorgchenani, M.; Gross, A.; Behm, R. J. Interaction of a self-assembled ionic liquid layer with graphite(0001): a combined experimental and theoretical study. J. Phys. Chem. Lett. 2016, 7, 226−233. (59) Koch, V.; Dominey, L.; Nanjundiah, C.; Ondrechen, M. J. The intrinsic anodic stability of several anions comprising solvent-free ionic liquids. J. Electrochem. Soc. 1996, 143, 798−803. (60) Ue, M.; Murakami, A.; Nakamura, S. A convenient method to estimate ion size for electrolyte materials design. J. Electrochem. Soc. 2002, 149−, A1385−A1388. (61) Vollmer, J. M.; Curtiss, L. A.; Vissers, D. R.; Amine, K. Reduction mechanisms of ethylene, propylene, and vinylethylene carbonates: a quantum chemical study. J. Electrochem. Soc. 2004, 151, A178−A183. (62) Budi, A.; Basile, A.; Opletal, G.; Hollenkamp, A. F.; Best, A. S.; Rees, R. J.; Bhatt, A. I.; O’Mullane, A. P.; Russo, S. P. Study of the initial stage of solid electrolyte interphase formation upon chemical reaction of lithium metal and n-methyl-n-propyl-pyrrolidinium-bis(fluorosulfonyl)imide. J. Phys. Chem. C 2012, 116, 19789−19797. (63) Camacho-Forero, L. E.; Smith, T. W.; Bertolini, S.; Balbuena, P. B. Reactivity at the lithium-metal anode surface of lithium-sulfur batteries. J. Phys. Chem. C 2015, 119, 26828−26839. (64) Sodeyama, K.; Yamada, Y.; Aikawa, K.; Yamada, A.; Tateyama, Y. Sacrificial anion reduction mechanism for electrochemical stability improvement in highly concentrated Li-salt electrolyte. J. Phys. Chem. C 2014, 118, 14091−14097. (65) Piper, D. M.; Evans, T.; Leung, K.; Watkins, T.; Olson, J.; Kim, S. C.; Han, S. S.; Bhat, V.; Oh, K. W.; Buttry, D. A.; Lee, S.-H.; et al. Stable silicon-ionic liquid interface for next-generation lithium-ion batteries. Nat. Commun. 2015, 6, 6230−6239. (66) Kresse, G.; Furhmuller, J. Software VASP, Vienna, 1999. Kresse, G.; Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B: Condens. Matter Mater. Phys. 1993, 47, R558−R561. (67) Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 54, 11169−11186. (68) Kresse, G.; Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 1996, 6, 15−50. (69) Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59, 1758−1775. (70) Lozano, A.; Escribano, B.; Akhmatskaya, E.; Carrasco, J. Assessment of van der Waals inclusive density functional theory methods for layered electroactive materials. Phys. Chem. Chem. Phys. 2017, 19, 10133−10139. (71) Chan, K. T.; Neaton, J. B.; Cohen, M. L. First-principles study of metal adatom adsorption on graphene. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 77, 235430.

(72) Henkelman, G.; Arnaldsson, A.; Jónsson, H. A fast and robust algorithm for Bader decomposition of charge density. Comput. Mater. Sci. 2006, 36, 354−360. (73) Sanville, E.; Kenny, S. D.; Smith, R.; Henkelman, G. Improved grid-based algorithm for Bader charge allocation. J. Comput. Chem. 2007, 28, 899−908. (74) Tang, W.; Sanville, E.; Henkelman, G. A grid-based Bader analysis algorithm without lattice bias. J. Phys.: Condens. Matter 2009, 21, 084204−084210. (75) Yu, M.; Trinkle, D. R. Accurate and efficient algorithm for Bader charge integration. J. Chem. Phys. 2011, 134, 064111−064118. (76) Ivaništšev, V.; Méndez-Morales, T.; Lynden-Bell, R. M.; Cabeza, O.; Gallego, L. J.; Varela, L. M.; Fedorov, M. V. Molecular origin of high free energy barriers for alkali metal ion transfer through ionic liquidgraphene electrode interfaces. Phys. Chem. Chem. Phys. 2016, 18, 1302− 1310. (77) Aguilera-Granja, F.; Gallego, L. Structural and electronic properties of Ni26−p X p clusters (X= Pd, Pt): a density-functionaltheoretic study. J. Appl. Phys. 2013, 114, 054311−054317. (78) Aguilera-Granja, F.; Carrete, J.; Vega, A.; Gallego, L. J. Structural, magnetic, and vibrational properties of stoichiometric clusters of CrN. Int. J. Quantum Chem. 2015, 115, 523−528.

H

DOI: 10.1021/acs.jpcc.8b02795 J. Phys. Chem. C XXXX, XXX, XXX−XXX