LinC6m → Lin-1C6m + Li+(EC) (R).
Methods
We applied density functional theory (DFT) + effective screening medium (ESM) method3 combined with reference interaction site method (RISM);4 ESM-RISM calculation treats graphite surface (LixC6 slab + Li+) and liquid solution (1 M LiPF6 EC) by quantum mechanics and implicit solvation model, respectively.5 ESM-RISM calculations were performed on the configuration of a vacuum/slab/solvent system as shown in Figure 1, where the DFT slab domain is on the left-hand side, and the RISM solvents (EC, Li+, and PF6−) treated by the RISM equations are on the right-hand side. The Li insertion/desorption path (reaction R) was assumed as a combination of one or two straight lines, as indicated in Figure 1. The reacting Li was moved from the stable Li site in the graphite layer (r = 0) to the interior region of the solution. In the ESM-RISM calculations, the distribution function of the solvation is automatically formulated for the solvation/desolvation structures.
Results
We successfully obtained the energy landscapes of reaction R under the various constant chemical potentials of electron; a profile of the grand potential Ω along the reaction R path as indicated in Figure 1. The activation energy of reaction R at the equilibrium potential is approximately 0.6 eV, which is consistent with the experimental measurements.1,2 In the presentation, we will discuss changes of the activation energy and the surface charge density with changing the Li content x of LixC6 (x=0, 0.5, 1) on two differently tilted graphite (30 and 90 degrees).
References
- T. Abe, H. Fukuda, Y. Iriyama, Z. Ogumi, J. Electrochem. Soc. 151, A1120 (2004).
- K. Xu, A. von Cresce, U. Lee, Langmuir 26, 11538 (2010).
- M. Otani, O. Sugino, Phys. Rev. B 73, 115407 (2006).
- F. Hirata, P. J. Rossky, Chem. Phys. Lett. 83, 329, (1981).
- S. Nishihara, M. Otani, Phys. Rev. B 96, 115429 (2017).